LS-DYNA使用指南第六章

LS-DYNA使用指南第六章发表时间：2007-7-30作者:安世亚太来源:e-works关键字:显式有限元LS-DYNAANSYS第六章接触表面ANSYS/LS-DYNA中的接触表面可以使用户在模型中诸Component之间定义多种接触类型，本章将概要地讲述一下显式动态分析中定义物理上的真实接触。必须注意的是显式动态分析中的接触与其它类型的ANSYS分析中的接触类型不同，在其它分析中，接触是由实际接触单元表示。而在显式动态分析中没有接触单元。只需定义接触表面，它们之间的接触类型以及相应的参数。6.1接触的定义因为在显式动态分析中会发生复杂的大变形，所以确定模型内component之间的接触是非常困难的。基于此原因，ANSYS/LS-DYNA程序中包含许多功能以使接触表面间的接触定义更容易些。在ANSYS/LS-DYNA中采用EDCGEN命令来定义所有接触表面。使用EDCGEN命令时遵循下列步骤：第一步；确定哪种接触类型最适合你的物理模型。第二步：定义接触实体。第三步：定义摩擦系数参数。第四步：为给定的接触类型给定一些附加输入。第五步：定义接触的杀死和激活时间。第一步：定义接触类型为了充分地描述在大变形接触和动态撞击中的复杂几何体之间的相互作用，在ANSYS/LS-DYNA中引入了许多种接触类型。这些接触类型，包括节点-表面，表面-表面，单面，单边，侵蚀，固连，固连断开，压延筋和刚性体接触，将在本章标题为“接触选项”中详细讨论，对于一般的分析而言，建议使用自动单面（ASSC），自动原则（AG），节点-表面（NTS），表面-表面（STS）接触选项。第二步：定义接触实体除单面接触（ASSC，SS和ESS）、自动通用（AG）和单边接触（SE）外，所有的接触类型都必须在发生接触的地方定义contact表面和target表面，这可用节点components,PARTID或部件集合ID定义。当使用contactcomponent和targetcomponent时，使用选择项并用CM命令把节点组合在一起（仅节点component有效），然后用下面的输入列表，说明如何使用EDCGEN命令在component之间定义接触，如第四章例题的球和球棒表面间的component.NSEL，S，NODE，．．．．！在球面上选择节点CM，BALLSURF，NODE！把被选的节点放在componentBALLSURF中NSEL，S，NODE，．．．．！选择球面上的节点CM，BATSURF，NODE！把被选节点放在componentBATSURF中EDCGEN，NTS，BALLSURF，BATSURF，．25，．23！在组元componentBALLSURF和componentBATSURF间定义为节点-表面接触。此外，还可以用有限元模型内当前定义的部件号或部件集合号来定义接触表面。部件集合号可以用EDASMP命令定义。下面的命令行说明了怎样使用EDCGEN命令在模型中定义不同部件或部件集合间的接触；EDCGEN，STS，1，2，．25，．23！在部件1和部件2间生成面面接触另外，结合PART/部件集合和组元定义，也可以定义接触和目标表面间的接触，表述如下：EDCGEN，NTS，N1，2，．3，．28！在组元N1和PART2间生成点面接触EDCGEN，ESTS，1，N2，．15，．15！在PART1和组元N2间生成侵蚀面面接触EDCGEN，STS，1，1，．1，．1！在PART1间生成面面接触如下例所述，也可以用EDCGEN命令定义部件集合间的接触：EDCGEN，STS，5，6，．3，．28！在部件集合5和6间生成面面接触在一些特定的单面接触类型（ASCC，AG，ESS，和SS）中无需定义contact和target表面，在本章后面将提及到，单面接触时最常用的接触类型，模型的全部外表面在整个分析中任一点都可能发生接触。程序在单面接触中将忽略任何contact和target表面的定义，并在执行EDCGEN命令时发出一个警告信息，一个典型单面接触命令如下：EDCGEN，ASSC，，，．34，．34！在整个模型中生成自动单面接触注--在显式分析中定义接触实体时，不允许有初始穿透。因此，定义接触组元时要注意。第三步：定义摩擦系数参数接触摩擦系数是由静态摩擦系数（FS），动摩擦系数（FD），和指数衰减系数（DC）来确定的。（FS，FD和DC可以用EDCGEN命令输入）假设摩擦系数与接触表面的相对速度有关：粘性摩擦系数，VC（用EDCGEN命令输入），可以限制最大摩擦力。公式表述如下：这里是接触时节点接触部分的接触面面积，VC值建议使用剪切屈服应力：这里是接触材料的屈服应力。为了避免在接触中产生不真实的震荡，对于薄板成形模拟，可使用垂直于接触表面的接触阻尼，接触阻尼系数计算如下：阻尼系数VDC—粘性阻尼系数（EDCGEN命令输入）。这里；K是界面刚度第四步：给定附加输入对于侵蚀，刚性体，固连断开和压延筋接触，还需输入一些其它数据（EDCGEN的V1-V4选项），这些数据因接触类型不同而各异，将分别简述如下：采用侵蚀表面接触（ENTS，ESS和ESTS）时，当表面实体单元发生失效时，需要在内部剩余单元重新定义接触。对于侵蚀接触，V1-V3定义如下：边界条件对称选项（V1）决定当单元表面失效时沿一个表面是否仍然保持对称性；内部侵蚀选项（V2）决定当外表面发生失效时沿内表面是否接着发生侵蚀；相邻材料选项（V3）决定当沿着自由表面发生失效时是否包括实体单元面。刚性体接触（RNTR和ROTR）通常用于多刚体动力学，在刚性体接触中，采用一条用户定义防止贯穿的力-挠度曲线，而不用线性刚度。因此，对于刚性体接触，数据曲线号（V1），用于给定刚性体接触的力计算方法类型选项（V2）和卸载刚度值（V3）必须输入。固连断开接触（TSTS和TNTS）用于表面胶合在一起时定义接触表面失效。对于所有固连断开接触类型，需输入拉伸失效应力（V1）和剪切失效应力（V2）定义失效准则。对于节点-表面固连断开接触。法向力指数项（V3）和剪切力指数项（V4）也需输入以定义失效准则。压延筋接触（DRAWBEAD）用于模拟压延筋的特殊情况，压延筋在深拉作业时有助于约束坯料。在压延筋接触中，必须输入一条载荷曲线号（V1），它给出作为压延筋位移函数的约束力弯曲分量，可以任选一条曲线号（V2），它给出作为压延筋位移函数的法向约束力以及沿压延筋的等距积分点数（V4）。第五步：定义激活或杀死时间对于每个接触定义，都可以用EDCGEN命令的BTIME域和DTIME域来定义杀死和激活时间。这允许在瞬态分析的任意时刻都激活接触，然后在稍后时间内杀死。6.1.1列表，显示和删除接触实体用EDCGEN命令定义接触后，就可以列表、显示或删除接触实体。用EDCLIST命令列出所有当前定义的接触实体。定义的每个接触都给定一个参考号用于显示接触实体，当前显式动态接触实体1一般的面面接触：节点接触实体N1，结点目标实体N2FS=0．10000FD=0．08000DC=0．00000VC=0．00000VDC=0．00002自动单面接触定义：模型的所有外表面FS=0．20000FD=0．15000DC=0．00000VC=0．00000VDC=0．0000用户可以采用EDPC命令选择和显示接触实体。显示将包括节点和单元，它与定义接触表面的方法有关（也就是说，部件或组元）。采用EDPC命令的MIN，MAX，和INC域来给定最小接触实体号，最大接触实体号和接触实体号增量。因此，对于上述说明执行EDPC，1，2，1，就可以选择显示STS和ASSC接触定义的实体。值得注意的是，对于单面接触定义，模型中所有外部表面都将被选择和显示出来。注--EDPC选择了给定接触实体的节点和单元。因此，在显示接触实体后，必须重新选择下步操作所需的所有节点和单元（SOLVE）。采用NSEL，ALL和ESEL，ALL命令（或其他命令的合适形式）。如果接触定义不正确，可以用EDDC命令删除它。为了删除指定的接触实体，可执行EDDC，DELE，Ctype,Cont,Targ,这里Ctype是接触类型，Cont和Targ表示已定义接触的接触和目标部件或组元。删除当前接触实体，执行EDDC，DELE，ALL。在小型重启动分析中不能删除接触实体（EDSTART，2）。但是，可以用EDDC命令激活（EDDC，DACT，Ctype,Cont,Targ）或杀死接触。当知道在什么阶段发生什么类型接触时，这个特点是非常有用的。使用此特征，必须在新分析中定义所有的接触实体（EDCGEN），并且必须至少执行一次小型启动，一旦定义了接触类型，就可以在不必要时杀死它，而在需要时重新激活。需要时才考虑接触，这样显著节省了CPU时间。（注意，如果在新分析中用EDDC，DELE删除接触，它的定义将从数据库中去除，那么就不能在后来的分析中激活它。）显式动态全启动分析不支持EDDC命令（EDSTART，3）。也就是说，在前面分析已定义的全启动中不能删除，杀死和激活接触定义。分页6.2接触选项为了充分定义在显式动态分析中表面间的复杂相互作用，在ANSYS/LS-DYNA中有24种接触类型（见下表）。在大量的接触类型中，需要用户对每一种接触类型都很了解，以便能正确地选用接触类型。因此，下面我们将讨论一下ANSYS/LS-DYNA中所有的接触类型。6．1接触类型SinglesurfaceNodestosurfaceSurfacetosurfaceNormalSSNTSSTS,OSTSSinglesurfaceNodestosurfaceSurfacetosurfaceAutomaticASSC,AG,ASS2DANTSASTSRigidRNTRROTRTiedTDNSTDSS,TSESTiedwithfailureTNTSTSTSErodingESSENTSESTSEdgeSEDrawbeadDRAWBEADFormingFNTSFSTS,FOSS6.2.1定义接触类型从表6.1中看出，在ANSYS/LS-DYNA程序中主要有三种基本接触类型：单面接触，节点-表面接触，表面-表面接触。1．单面接触（SS，ASSC，AG，ASS2D，ESS）单面接触用在一个物体表面的自身接触或它与另一个物体表面接触，在单面接触中，ANSYS/LS-DYNA程序将自动判定模型中哪处表面发生接触。因此，单面接触的定义是最简单的，无需定义contact和target表面，当定义好单面接触时，它允许一个模型的所有外表面都可能接触，这对于预先不知道接触表面的自身接触或大变形问题很有用处。与隐式模型过多定义接触面将大大增加CPU时间不同，在显式模型中定义单面接触只会较少的增加CPU时间，许多碰撞和撞动态碰撞问题都需定义单面接触。由于自动通用接触（AG）很有效，它包括壳边接触（SE）和改进的梁接触，因此，推荐你在难以预测接触条件时，对于自身接触和大变形问题优先选择此种接触类型。2．点-面接触（NTS，ANTS，RNTR，TDNS，TNTS，ENTS，DRAWBEAD，FNTS）点面接触类型是接触节点将穿透target表面。这种接触类型通常用于一般两个表面间的接触。采用ANSYS隐式程序中同样的规则，需要定义target表面及contact表面。·平面或凹面为target表面，凸面为contact表面·粗网格所在面作为target面，细网格所在面为contact面在定义压延筋接触时，筋总为contact表面，而板料为target表面。3．面-面接触（STS，OSTS，ASTS，ROTR，TDSS，TSTS，ESTS，SE，FSTS，FOSS，TSES）当一个物体的表面穿透另一个物体的表面时需使用面-面接触，面-面接触类型是最常用的，并且常用于任意形状且存在较大接触面积的物体接触。这种接触类型对于物体间有大量相对滑移时很有效，例如块在平板上滑动，球在槽内滑动等。6.2.2定义接触选项对于以上三种接触类型的每一种又含有多个接触类型选项，在A