abaqusexplicit接触问题.docx

abaqusexplicit接触问题.docx

《abaqusexplicit接触问题.docx》由会员分享，可在线阅读，更多相关《abaqusexplicit接触问题.docx（15页珍藏版）》请在冰豆网上搜索。

abaqusexplicit接触问题

1.Abaqus/Explicit中的接触形式

双击Interactions，出现接触形式定义。

分为通用接触（Generalcontact）、面面接触（Surface-to-Surfacecontact）和自接触（Self-contact）。

1.通用接触Generalcontact

通用接触用于为多组件，并具有复杂拓扑关系的模型建模。

Generalcontactalgorithm

Thecontactdomainspansmultiplebodies（bothrigidanddeformable）

Defaultdomainisdefinedautomaticallyviaanall-inclusiveelement-basedsurface

Themethodisgearedtowardmodelswithmultiplecomponentsandcomplextopology。

Greatereaseindefiningcontactmodel

2.Surface-to-Surfacecontact

Contactpairalgorithm

Requiresuser-specifiedpairingofindividualsurfaces

Oftenresultsinmoreefficientanalysessincecontactsurfacesarelimitedinscope

3.自接触（Self-contact）

自接触应用于当部件发生变形时，可能导致自己的某两个或多个面发生接触的情况。

如弹簧的压缩变形，橡胶条的压缩。

容易使用“自动接触”节省生成模型的时间通用接触算法一般比双面接触算法快

机械约束形式

运动依从Kinematiccontactmethod

（只有接触对形式可用,Generalcontact不可用）

默认的运动接触公式达到的计算精度与接触条件相一致。

在多数情况下，它工作得很好。

但是在某些情况下，如抖动接触，使用罚函数接触会更容易得到收敛的解。

不能为刚体-刚体接触建模。

罚函数（通用接触和接触对均可使用）Penaltycontactmethod

罚函数接触算法中接触约束的严格性低于运动学算法。

罚函数算法可以处理更一般类型的接触；比如，刚体之间的接触。

因为罚函数算法在模型中引入附加刚度，该附加刚度将影响稳定时间增量。

Abaqus/Explicit自动计算由与接触力相关的侵彻距离引入的“弹簧”刚度或“罚”刚度。

但是必须考虑相关的影响：

应该尽量减少对稳定时间增量的影响。

在所有的分析中，允许的侵彻不可以太大。

对于接触对算法：

通过在*CONTACTCONTROLS选项指定SCALEPENALTY参数，用户可以缩放默认的罚刚度。

对于通用接触算法：

可以使用*CONTACTCONTROLSASSIGNMENT,TYPE=SCALEPENALTY选项缩放罚刚度。

综上，多数情况下，运动依从算法和罚函数算法将得到几乎一致的结果。

然而，在某些情况下，一种方法可能比另一种方法更可取。

如果需要知道kinematicalgorithms和penaltyalgorithms之间的区别，那就需要啃理论教材了（当然看abaqus帮助文件也可以）.

简略的说：

1）penaltyalgorithms（罚函数法）

在被abaqus检测到接触距离以内的节点之间定义罚刚度，进而来迭代出接触力。

在每一个时间步先检查各从节点是否穿透主面。

如没有穿透则进入下一个时间步；如果穿透，则在该从节点与被穿透主面间定义罚刚度，引入一个较大的界面接触力，其大小与穿透深度、主面的刚度成正比。

2）kinematicalgorithms

先在一定时间内检查所有未与主面（mastersurface）接触的从节点（slavenode），看是否在此时间内穿透了主面。

如果存在穿透则缩小时间步，使那些穿透主面的从节点都不贯穿主面，而使其正好到达主面。

在计算下一时间步之前，对所有已经与主面接触的从节点都施加约束条件，以保持从节点与主面接触而不贯穿。

此外还应检查那些和主面接触的从节点所属单元是否受到拉应力作用。

如受到拉应力，则施加释放条件，使从节点脱离主面。

就abaqus接触计算的初级使用来说，如果接触对中有刚体，一般用penaltyalgorithms；如果接触对是弹性体，有限选Kinematicalgorithms.无刚体的接触，建议直接用Kinematicalgorithms试算。

总结：

1）接触面存在刚体，用penalty算法

2）接触面之间的相对运动方向平行于接触面，用penalty算法

3）弹性碰撞分析不可用hard-kinematic算法。

碰撞之中塑性变形控制的接触分析

4）penalty算法一般会减小稳定时间增量

5）接触计算与以下情况耦合时推荐penalty算法：

constraintequation,multi-pointconstraint,tieconstraint,embeddedelementconstraint,orkinematicconstraint

6）breakablebondmodel模拟，必须用hardkinematiccontact。

碰撞应该用罚函数法吧，要是用运动学算法时需要划分较细的网格以避免能量的丢失。

kinematiccontact要求更严格，不允许出现穿透，而penaltycontact允许。

默认的Penaltystiffness为弹性刚度的10%，所以如果是解决塑性变形问题的话，二者算出来的结果相差不大。

penalty对切向的相对位移的接触比较有效。

3.初始运动学依从关系

Abaqus/Explicit不允许接触表面的初始过盈。

接触表面的节点将被调整，删除分析之前的初始过盈：

只移动接触表面的节点。

对于分析过程第一个分析步定义的接触对，由于调整表面引起的位移不产生初始应变或应力。

在随后的分析步中：

对于接触对算法，调整将产生应变。

对于通用接触算法，忽略初始过盈。

ContactFormulations

接触程式，包括

1.Contactdiscretization接触离散化

Whereistheconstraintapplied接触约束应用在哪

Node-to-surface

Surface-to-surface

2.Constraintenforcement执行接触

Howistheconstraintenforced接触约束怎么执行

Default（Explicit）

Direct（Lagrangemultipliers）Standard

PenaltymethodStandard

AugmentedLagrange（Lagrangemultiplierscombinedwithpenaltymethod）Standard

3.Contacttracking（relativesliding）跟踪接触

Howdoestheconstraintevolve接触约束怎么发展

Finitesliding

Smallsliding

ContactDiscretization

Node-to-surfacetechnique:

节点和面接触

Defaultmethodforcontactpairs。

接触对的默认接触定义方法。

Notavailableforgeneralcontact。

不适用于Generalcontact。

Nodesononesurface（theslavesurface）contactthesegmentsontheothersurface

（themastersurface）.从面上的节点接触主面的面。

Contactisenforcedatdiscretepoints（slavenodes）。

在从面节点处执行接触。

从面需要划分较细的单元来提高精度和消除应力振动即消除应力噪声。

Useofa“matchedmesh”acrossthecontactinterfacewilleliminatethissolutionnoise。

Surface-to-surfacetechnique：

面面接触

Onlymethodforgeneralcontact，Alternativemethodforcontactpairs

Themethodconsiderstheshapeofboththemasterandslavesurfaces.

Contactisenforcedinanaveragesenseovertheslavesurface.

当面面接触用于接触对时，

contactpairs和generalcontact区别不大，主要的区别在于使用方便性和性能表现erformance。

可以看出Node-to-Surface参与接触节点个数较少。

而S-to-S较多。

对于面面接触，个别点小的穿透可能会出现，但是不会发生漏检测到的主面穿透

ContactEnforcement

Hardcontact:

InAbaqus/Standard,thedefaultcontactbehavioris“hard”contact.

Thebehaviorisdescribedbyacontactpropertyknownasthepressure-penetrationcurve（alternativebehaviorcanbespecified;discussedlater）.

Thedesiredbehavior（nopenetration）isachievedusinganenforcementmethod.

理想的没有穿透的接触行为可以用三种强化方法来获得。

ThreenumericalmethodsareavailableinAbaqus/Standardtoachieveorapproximate“hard”contactconditions:

Directenforcementmethod:

Strictenforcementofpressure-penetrationrelationshipusingtheLagrangemultipliermethod.严格压力-穿透关系，用拉格朗日乘法算法。

Prosandconsofdirectenforcement：

优缺点

Advantages:

Accuracy—constraintissatisfiedexactly.

Disadvantages:

Addstoequationsolvercost.求解代价高

Additionalvariablepercontactconstraint,whichenlargesthesystemofequationstobesolved.

Restrictseliminationorderforsparsesolver,whichcandegradeperformance.

Potentialconvergencedifficulties.潜在的收敛问题。

Abruptchangefromzerocontactstiffness（whilecontactisinactive）toinfinitecontactstiffness（whilecontactisactive）.从0刚度到无穷刚度变化的剧烈。

Difficultieswithoverconstraints.

OverlapbetweencontactconstraintsandMPC’s,etc.

Penaltymethod:

Approximateenforcementusingpenaltystiffness.用罚刚度近似。

Thepenaltymethodisastiffapproximationofhardcontact。

罚刚度可以是：

Linear（default）:

线性

Easierconvergence.容易收敛

Bettersuitedforproblemsinvolvingfirmcontact.较适用于硬面接触。

Difficulttochoosestiffnessappropriateforallregimes.难以选择刚度值适用于所有情况。

Nonlinear:

非线性

Thelowerinitialstiffnessmakesitbettersuitedforproblemsinvolvingchattering.初始的小刚度值适合于解决包括颤振的接触。

Thehigherfinalstiffnesshelpsreducepenetrations.高的最终刚度降低穿透

Convergenceoverallcanbemoredifficult.收敛较困难。

Defaultpenaltystiffness：

默认的刚度值

Abaqus折中选择刚度：

刚度小，过度穿透。

刚度大，收敛速度；LagrangemultiplierDOFareneededtoavoidill-conditioning.

Thedefaultpenaltystiffnessisbasedonarepresentativestiffnessoftheunderlyingelements.

Ascalefactorisappliedtothisrepresentativestiffnesstosetthedefaultpenaltystiffness;itsmagnitudeishigherinAbaqus/StandardthaninAbaqus/Explicit.

默认刚度不合适，可以设置比例系数调节。

Order-of-magnitudechangesarerecommended.推荐

Ifthescalefactor>100,AbaquswillautomaticallyinvokeavariantofthemethodthatusesLagrangemultiplierstoavoidill-conditioningissues.如果比例系数大于100，自动采用Lagrangemultipliers来避免病态。

罚刚度值：

Stifforblockyproblems：

硬的或短而结实的接触。

默认刚度值产生结果：

比的上directmethod的精确度。

RequireslessmemoryandCPUtime。

Bending-dominatedproblems：

弯曲占主导地位的接触：

默认刚度值会被缩减但是不显着影响精确度。

Scalebacktwoordersofmagnitudewithoutanysignificantlossofaccuracy.

Scalingbackthepenaltystiffnessforbending-dominatedproblemssometimesdramaticallyincreasestheconvergencerate。

缩减刚度值还或许会提高收敛率。

可以看出对于Hertz接触，默认的罚刚度值和DirectLagrange结果最接近。

Advantages:

Significantlyimprovedconvergencerates

Betterequationsolverperformance.NoLagrangemultiplierDOFunlesscontactstiffnessisveryhigh.

Goodtreatmentofoverlappingconstraints.

Disadvantages:

Smallamountofpenetration.Typicallyinsignificant

Mayneedtoadjustthepenaltystiffnessinsomecases

使用

Usedbydefaultforeachofthefollowing:

Generalcontact、Contactpairswiththefinite-sliding,surface-to-surfaceformulation。

AutomaticcontactpairdetectioncapabilityinAbaqus/CAE。

AugmentedLagrangemethod:

Approximateenforcementusingpenaltymethodwithaugmentationiterations;notdiscussedfurtherhere.

Usedbydefaultfor:

3Dselfcontactwithnode-to-surfacediscretization

Contacttracking（relativesliding）

Finitesliding

Kinematically,finiteslidingofdeformablebodiesagainsteachotherallowsthemostgeneralrelativemotionbetweenthebodies:

Arbitrarilylargeslidingisallowed.

Arbitrarilylargerotationsanddeformationsofthesurfacesareallowed.

任意大的滑动和转动都允许。

接触对和通用接触中均可以使用。

Smallsliding

Smallslidingassumesthattherelativemotionperslavenoderemainssmallcomparedtothelocalcurvatureofthemastersurface.

假设从面上每个节点的相对运动相对于主面的接触段的曲率很小。

Advantage:

Simplifiesthenumericalmethod.Potentialforreducedcostperiterationandfindingaconvergedsolutioninfeweriterations.

优点：

简化数值求解过程。

收敛性好。

Disadvantage:

Theresultscanbenon-physicaliftherelativetangentialmotiondoesnotremainsmall.

缺点：

当从面节点的相对切线运动不是很小时会产生不合物理规律的错误。

Invokingsmall-sliding（contactpairsonly）:

仅在接触对中才能用。

Thesmallslidingformulationcannotbespecifiedwiththepenaltyconstraintformulation。

不能用于Penaltycontact.

如果不想考虑是否满足small-sliding的假设，则选择finite-sliding。

Analyticalrigidsurfaces

解析刚体面：

Analyticalrigidsurfaces

Canonlybeusedwithcontactpairs只能用于接触面。

Surfacegeometryoflimitedshape面形状有限

Preferredmethodismeshingabodywithdeformableelementsandthenapplyingarigidbodyconstraint.用可变形单元划分然后施加刚体约束。

本模型中采用的就是推荐的方法。

FrictionEnforcement

1.TwofrictionenforcementmethodsareimplementedinAbaqus:

Penaltymethod（default）

LagrangemultipliermethodStandardonly

Themostcommon（anddefault）combinationofmodel/methodisisotropicCoulombfrictionusingthepenaltymethod.

最常用和默认的就是各向同性库伦摩擦，用罚函数法。

2.