整理ANSYS显式动态分析求解步骤

第一章引言ANSYS/LS-DYNA将显式有限元程序LS-DYNA和ANSYS程序强大的前后处理结合起来。用LS-DYNA的显式算法能快速求解瞬时大变形动力学、大变形和多重非线性准静态问题以及复杂的接触碰撞问题。使用本程序，可以用ANSYS建立模型，用LS-DYNA做显式求解，然后用标准的ANSYS后处理来观看结果。也可以在ANSYS和ANSYS-LS-DYNA之间传递几何信息和结果信息以执行连续的隐式一显式/显式一隐式分析，如坠落实验、回弹、及其它需要此类分析的应用。1.1显式动态分析求解步骤概述显式动态分析求解过程与ANSYS程序中其他分析过程类似，主要由三个步骤组成：1：建立模型（用PREP7前处理器）2：加载并求解（用SOLUTION处理器）3：查看结果（用POST1和POST26后处理器）本手册主要讲述了ANSYS/LS-DYNA显式动态分析过程的独特过程和概念。没有详细论述上面的三个步骤。如果熟悉ANSYS程序，已经知道怎样执行这些步骤，那么本手册将提供执行显式动态分析所需的其他信息。如果从未用过ANSYS，就需通过以下两本手册了解基本的分析求解过程：•ANSYSBasicAnalysisGuide•ANSYSModelingandMeshingGuide使用ANSYS/LS-DYNA时，我们建议用户使用程序提供的缺省设置。多数情况下，这些设置适合于所要求解的问题。1.2显式动态分析采用的命令在显式动态分析中，可以使用与其它ANSYS分析相同的命令来建立模型、执行求解。同样，也可以采用ANSYS图形用户界面（GUI）中类似的选项来建模和求解。然而，在显式动态分析中有一些独特的命令，如下：EDADAPT:激活自适应网格

EDASMP:创建部件集合EDBOUND:定义一个滑移或循环对称界面EDBVIS:指定体积粘性系数EDBX:创建接触定义中使用的箱形体EDCADAPT:指定自适应网格控制EDCGEN:指定接触参数EDCLIST:列出接触实体定义EDCMORE:为给定的接触指定附加接触参数EDCNSTR:定义各种约束EDCONTACT:指定接触面控制EDCPU:指定CPU时间限制EDCRB:合并两个刚体EDCSC:定义是否使用子循环EDCTS:定义质量缩放因子EDCURVE:定义数据曲线EDDAMP:定义系统阻尼EDDC:删除或杀死/重激活接触实体定义EDDRELAX:进行有预载荷几何模型的初始化或显式分析的动力松弛EDDUMP:指定重启动文件的输出频率（d3dump）EDENERGY:定义能耗控制EDFPLOT:指定载荷标记绘图

EDHGLS:定义沙漏系数EDHIST:定义时间历程输出EDHTIME:定义时间历程输出间隔EDINT:定义输出积分点的数目EDIS:定义完全重启动分析的应力初始化EDIPART:定义刚体惯性EDLCS:定义局部坐标系EDLOAD:定义载荷EDMP:定义材料特性EDNB:定义无反射边界EDNDTSD:清除噪声数据提供数据的图形化表示EDNROT:应用旋转坐标节点约束EDOPT:定义输出类型，ANSYS或LS-DYNAEDOUT:定义LS-DYNAASCII输出文件EDPART:创建，更新，列出部件EDPC:选择、显示接触实体EDPL:绘制时间载荷曲线EDPVEL:在部件或部件集合上施加初始速度EDRC:指定刚体/变形体转换开关控制EDRD:刚体和变形体之间的相互转换EDREAD:把LS-DYNA的ASCII输出文件读入到POST26的变量中

EDRI:为变形体转换成刚体时产生的刚体定义惯性特性EDRST:定义输出RST文件的时间间隔EDSHELL:定义壳单元的计算控制EDSOLV:把“显式动态分析”作为下一个状态主题EDSP:定义接触实体的小穿透检查EDSTART:定义分析状态（新分析或是重启动分析）EDTERM:定义中断标准EDTP:按照时间步长大小绘制单元EDVEL:给节点或节点组元施加初始速度EDWELD:定义无质量焊点或一般焊点EDWRITE:将显式动态输入写成LS-DYNA输入文件PARTSEL:选择部件集合RIMPORT:把一个显式分析得到的初始应力输入到ANSYSREXPORT:把一个隐式分析得到的位移输出到ANSYS/LS-DYNAUPGEOM:相加以前分析得到的位移，更新几何模型为变形构型关于ANSYS命令按字母顺序排列的详细资料（包括每条命令的特定路径），请参阅《ANSYSCommandsReference》。1.3本手册使用指南本手册包含过程和参考信息，可从前到后选择性阅读。然而，选择与规划和进行显式动态分析求解过程相对应的顺序阅读更有帮助。在建模之前，必须确定最能代表物理系统的单元类型和材料模型，下面几章将为你提供相关的一些基础知识：

第二章，单元第七章，材料模型选择了合适的单元类型和材料模型后，就可以建模了。建模的典型方面如下所示:第三章，建模第六章，接触表面第八章，刚体第四章，加载与求解和后处理有关的特征如下：第五章，求解特性第十二章，后处理有些高级功能在一个分析中可能涉及不到，但在某些情况下可能用到，如下所示:第九章，沙漏第十章，质量缩放第十一章，子循环第十三章，重启动第十四章，显-隐式连续求解第十五章，隐-显式连续求解最后，附录中还包含了有关下列主题的有关信息:附录A,隐、显式方法的比较附录B,材料模型样例附录C，ANSYS/LS-DYNA和LS-DYNA命令变换

1.4何处能找到显式动态例题TheExplicitDynamicsTutorial描述了一个典型的显式动态分析例题。1.5其它信息对于显式动态分析的详细资料，请参阅《ANSYSStrueturalAnalysisGuide》中的第十四章。对于显式动态分析单元的详细资料，请参阅《ANSYSElementReference》；至于详细的理论信息，请参阅LivermoreSoftwareTechnologyCorporation的《LS-DYNATheoreticalManual》。第二章单元在显式动态分析中可以使用下列单元：LINK160杆BEAM161梁・PLANE162平面SHELL163壳・S0LID164实体・C0MBI165弹簧阻尼・MASS166质量・LINK167仅拉伸杆本章将概括介绍各种单元特性，并列出各种单元能够使用的材料类型。除了PLANE162之外，以上讲述的显式动态单元都是三维的，缺省时为缩减积分（注意:对于质量单元或杆单元缩减积分不是缺省值）缩减积分意味着单元计算过程中积分点数比精确积分所要求的积分点数少。因此，实体单元和壳体单元的缺省算法采用单点积分。当然，这两种单元也可以采用全积分算法。详细信息参见第九章沙漏也可参见《LS-DYNATheoreticalManual》。

这些单元采用线性位移函数；不能使用二次位移函数的高阶单元。因此，显式动态单元中不能使用附加形状函数，中节点或P-单元。线位移函数和单积分点的显式动态单元能很好地用于大变形和材料失效等非线性问题。值得注意的是，显单元不直接和材料性能相联系。例如，S0LID164单元可支持20多种材料模型，其中包括弹性，塑性，橡胶，泡沫模型等。如果没有特别指出的话（参见第六章，接触表面），所有单元所需的最少材料参数为密度，泊松比，弹性模量。参看第七章材料模型，可以得到显式动态分析中所用材料特性的详细资料。也可参看《ANSYSElementReference》，它对每种单元作了详细的描述，包括单元的输入输出特性。2.1实体单元和壳单元SOLID164SOLID164单元是一种8节点实体单元。缺省时，它应用缩减（单点）积分和粘性沙漏控制以得到较快的单元算法。单点积分的优点是省时，并且适用于大变形的情况下。当然，也可以用多点积分实体单元算法（KEY0PT（l）=2）；关于SOLID164的详细描述,请参见《ANSYSElementReference》和《LS-DYNATheoreticalManual》中的§3.3节。如果担心沙漏现象，比如泡沫材料，可采用多点积分算法，因为它无需沙漏控制；计算结果要好一些。但要多花大约4倍的CPU时间。楔形、锥型和四面体单元是六面体单元的退化产物（例如，一些节点是重复的）。这些形状在弯曲时经常很僵硬，有些情况下还有可能产生问题。因此，应尽量避免使用这些退化形状的单元。对于实体单元可采用下列材料模型：•各向同性弹性•正交各向异性弹性•各向异性弹性•双线性随动强化•塑性随动强化•粘弹性

•Blatz-ko橡胶•双线性各向同性•幕律塑性•应变率相关塑性•复合材料破坏•混凝土破坏•地表材料•分段线性塑性•Honeycomb蜂窝材料•Mooney-Rivlin橡胶•Barlat各向异性塑性•弹塑性流体动力•闭合多孔泡沫•低密度泡沫•粘性泡沫•可压缩泡沫•应变率相关幕律塑性•Johnson-Cook塑性•空材料•Zerilli-Armstrong•Bamman

•Steinberg•弹性流体SHELL163SHELL163单元有12中不同的算法。用KEYOPT(1)来定义所选的算法。和实体单元一样，积分点的个数直接影响着CPU时间。因此，对于一般的分析而言，建议使用缺省积分点个数。以下将概述SHELL163单元的不同算法：2.1.3通用壳单元算法・Belytschko-Tsay(KEY0PT(1)=0或2)—缺省—速度快，建议在多数分析中使用—使用单点积分—单元过度翘曲时不要使用•Belytschko-Wong-Chiang(KEY0PT(1)=10)—比BelytschKo-Tsay慢25%—使用单点积分—对翘曲情况一把可得到正确结果•Belytschko-Leviathan(KEY0PT(1)=8)—比BelytschKo-Tsay慢40%—使用单点积分—自动含有物理上的沙漏控制它可以将节点偏离单元的中-Hughes-Liu(KEY0PT(1)=1,6,7,11)有4种不同的算法,面。它可以将节点偏离单元的中KEY0PT(1)=1一般型Hughes-Liu，使用单点积分，比BelytschKo-Tsay慢250%。

KEYOPT(1)=11快速Hughes-Liu，使用单点积分，比Belytschko-Tsay慢150%。比Belytschko-Tsay慢比Belytschko-Tsay慢20倍。KEYOPT(1)=7S/R快速Hughes-Liu，有4个积分点，没有沙漏，比Belytschko-Tsay慢&8倍。如果分析中沙漏带来麻烦的话，建议使用此算法。KEYOPT(1)=12全积分Belytschko-Tsay壳。在平面内有四个积分点，无需沙漏控制。通过假设的横向剪切应变可以矫正剪切锁定。但是它比单点Belytschko-Tsay慢2.5倍，如果分析中担心沙漏的话，建议使用此方法。2.1.4薄膜单元算法•Belytschko-Tsay薄膜(KEYOPT(1)=5)—速度快，建议在大多数薄膜分析中使用—缩减(单点)积分—很好地用于关心起皱的纺织品(例如，大的平面压缩应力破坏较薄的纤维单元)•全积分Belytschko-Tsay薄膜(KEYOPT(1)=9)—明显的比通用薄膜单元慢(KEYOPT(1)=5)—面内有四个积分点—无沙漏2.1.5三角型薄壳单元算法・C0三角型薄壳(KEYOPT(1)=4)单元—基于Mindlin-Reissner平板理论—该构型相当僵硬，因此不建议用它来整体划分网格—使用单点积分•BCIZ三角型薄壳(KEYOPT(1)=3)单元

—基于Kirchhoff平板理论—比C0三角型薄壳单元慢—使用单点积分ANSYS/LS-DYNA用户手册中有关SHELL163的描述对可用的壳单元算法作了完整的介绍。退化的四边形单元在横向剪切时易发生锁死。因此，应使用Co三角型薄壳单元（基于Belytschko和其合作者的工作），如果在同一种材料中把单元分类标记（EDSHELL命令的ITRST域）设置为1（缺省值），就可混合使用四边形和三角形单元。对于壳单元可使用以下材料模型：•各向异性弹性•正交各向异性弹性•双线性随动强化•塑性随动强化•Blatz-Ko橡胶•双线性各向同性•幕律塑性•应变率相关塑性•复合材料破坏•分段线性塑性•Mooney-Rivlin橡胶•Barlat各向异性塑性・3参数Barlat塑性•横向各向异性弹塑性

・应变率相关幕律塑性-横向各向异性FLD•Johnson-Cook塑性•Bamman注意一当SHELL163单元使用Mooney-Rivlin橡胶材料模型时，LS-DYNA编码将自动使用Belytschko-Tsay算法的完全拉格朗日修正法来代替KEYOPT（1）指定的算法。程序选择的算法要求满足超弹材料的特殊需要。图2-1积分点所有的壳单元算法沿厚度方向都可以有任意多个积分点。典型地，对于弹性材料沿厚度方向需要2个积分点，而对于塑性材料则需要3个或更多的积分点。沿厚度方向的积分点个数由第二实常数来控制：R,NEST,,R2,这里R2为积分点的个数（NIP）。壳单元使用三维平面应力本构子程序修正应力张量,使垂直于壳单元中面的正应力分量为零。积分点位于壳单元的质心垂线上，见图2-1。开始时每个节点的厚度方向与单元表面都是正交的但它们随节点旋转。计算弯矩和平面力需要厚度方向的积分点。其应变呈线性分布，而应力分布要复杂得多，它和材料性质有关。对于线弹性材料两个积分点就足够了，而非线性材料则需要更多的积分点，输出的应力属于最外层的积分点，而不是表面上的（尽管后处理的术语是指顶面和底面），因此在分析结果时需要注意，对于弹性材料，应力可以外推到表面上。对于非线性材料来说，通常是选

择沿厚度方向的四五个节点而忽略其不精确性(例如，忽略表面和外部积分点之间的应力差)。高斯积分法最外层积分点的位置由下表给出：中面0最外积分点两点三占•-八、、四点五点0.57740.77460.86110.9062外表面1.000注意一在使用线弹性材料时，能够预先准确定义这些积分准则，但是通常在ANSYA/LS-DYNA中无法做到，由于模拟大多涉及非线性行为。另外，对于全积分单元来说，其输出应力是同一层内2X2积分点的应力平均值。2.1.6PLANE162PLANE162单元是一个二维，4节点的实体单元，它既可以用作平面(X-Y平面)单元，也可以用作轴对称单元(Y轴对称)。KEYOPT(3)用来指定单元的平面应力、轴对称和平面应变选项。对于轴对称单元可以利用KEYOPT(2)指定面积或体积加权选项。PLANE162典型情况下为四节点单元。当然也可以用三节点三角形选项，但是由于它太僵硬，所以不推荐使用它。这个单元没有实常数。重要的是要注意到含有PLANE162单元的模型必须仅包含这种单元。ANSYS/LS-DYNA中不允许有二维和三维单元混合使用的有限元模型。这种单元可用的材料模型与KEYOPT(3)的设置有关。对KEYOPT(3)=0，1，2(平面应力、平面应变或轴对称)，用户可以选择下列材料模型：•各向同性弹性•正交各向异性弹性•Blatz-ko橡胶•Mooney-Rivlin橡胶•粘弹性

•双线性各向同性•双线性随动强化•塑性随动强化•幕率塑性•应变率相关幕率塑性•应变率相关塑性•分段线性塑性•复合材料破坏•Johnson-Cook塑性•Bamman对平面应力选项(KEY0PT(3)=0)，可以选择下列材料：・3参数Barlat塑性•Barlat各向异性塑性•横向正交各向异性弹塑性•横向正交异性FLD对轴对称和平面应变选项(KEYOPT(3)=1或2),可以选用下列材料:•正交各向异性弹性•弹塑性流体动力•闭合多孔泡沫•低密度泡沫•可压缩泡沫

•Honeycomb蜂窝材料•空材料•Zerilli—Armstrong•Steinberg•弹性流体2.2梁单元和杆单元BEAM161BEAM161有两种基本算法：Hughes-Liu和Belytschko-Schwer。因为BEAM161不产生任何应变，所以它最适合于刚体旋转。必须用三个节点来定义单元；在每个端点处有一节点，同时需要有一定向节点。对于这两种算法来说，可用KEY0PT(4)和KEY0PT(5)来定义几种横截面。通常，对于2X2高斯积分点，BEAM161具有高效和耐用性。可用KEYOPT(2)来定义不同积分算法。Hughes-Liu梁单元(缺省值)是一个传统积分单元，它可以采用梁单元中间跨度的一组积分点来模拟矩形和圆形横截面。另外，用户也可以定义一个横截面积分规则来模拟任意的横截面。梁单元沿其长度方向能有效地产生一个不变力矩，因此，与实体单元和壳体单元一样，网格必须合理划分以保证精度。由于积分点的位置，只在单元中心才可检验屈服，因此，由于必须在夹持单元的中心处产生全塑性力矩而不是单元外边根部，悬臂梁模型将在一个稍高的力作用下产生屈服。Belytschko-Schwer.梁单元(KEY0PT(1)=2，4，5)是一个显式算法，可以产生一个沿长度方向呈线性分布的力矩。这种单元有“正确”的弹性应力并且在其末端可检验屈服。例如：当一个悬臂梁在端部静态加载时，可用一个单元来精确地表达弹性和塑性状态。如同Hughes-Liu梁单元，质量堆积到节点上，因此，在动态问题中必须要细分网格，因为此时正确的质量分布是很重要的。对于梁单元，可使用下列材料模型：(对于某些算法有些限制)•各向同性弹性•双线性随动强化

•塑性随动强化•粘弹性•幕率塑性•分段线性塑性LINK160LINK160桁架单元与Belytschko-Schwer梁单元很相似，但只能承受轴向载荷。这种类型单元支持直杆，在两端轴向加载，材料性质均匀。对于这种单元可使用的材料类型为各向同性弹性，塑性随动强化（率相关）和双线性动力。LINK167LINK167单元是仅能拉伸的杆，可以用于模拟索。它与弹性单元类似，由用户直接输入力与变形的关系。本单元类型需要用EDMP命令来定义索单元选项（参看EDMP命令概述）。2.3离散单元COMBI165弹簧-阻尼单元弹簧单元因位移产生一个力；也就是说改变单元的长度产生力。力沿单元轴向加载。例如，拉力在节点1上是沿轴的正方向，而对节点2是沿轴的负方向。缺省时，单元轴的方向就是从节点1到节点2。当单元旋转时，力作用方向线也将随之而旋转。阻尼单元可认为是弹簧单元的一种：可模拟线性粘性和非线性粘性阻尼。也可使用旋转（扭转）弹簧和阻尼单元，这些可通过KEY0PT（1）来选择，其他输入部分和平移弹簧一样；给定的力-位移关系可认为是力矩-转角（为弧度单位）关系，力矩施加方向沿单元的轴向方向（顺时针为正）。旋转弹簧单元只影响其节点的旋转自由度一它们并不把节点铰接在一起。COMBI165单元可和其它显式单元混合使用。然而，由于它没有质量，在分析中不能只有COMBI165一种类型单元，为了表达一个弹簧/质量系统，必须定义MASS166单元来加上质量。对于同一个COMBI165单元不能同时定义弹簧和阻尼特性。但是，可以分别定义使用同样节点的弹簧和阻尼单元（也就是说，可以重叠两个COMBI165单元）。

对于C0MBI165单元可以使用下列材料模型：•线弹性弹簧•线粘性阻尼•弹性塑料弹簧•非线性弹性弹簧•非线性粘性阻尼•通用非线性弹簧•麦克斯韦粘弹性弹簧•无弹性拉伸或仅压缩弹簧使用COMBI165单元时，应该给每一零件分别指定唯一的实常数，单元类型和材料特性（分别是R,ET和TB命令）从而保证每个零件都分别定义。MASS166质量单元由一个单节点和一个质量值定义（力X时间2/长度）。质量单元通常用于模拟一个结构的实际质量特性，而没有把大量实体单元和壳体单元包括进去。例如，在汽车碰撞分析中，质量单元可以模拟发动机部分，主要感兴趣的不是它的变形性质。采用质量单元将减少分析所需的单元数目，因而减少求解所需的计算时间。用户也可用MASS166单元来定义一个节点的集中转动惯量。如使用这一选项，可在MASS166单元定义中设置KEYOPT（1）=1并且通过单元实常数输入六个惯性矩值（IXX,IXY,IXZ,IYY，IYZ，IZZ）。这个选项不能输入质量值；所以，必须在同一个节点定义第二个质量单元来说明质量（KEYOPT（1）=0）。2.4一般单元特性以下几种单元可被定义为刚性体：LINK160,BEAM161,PLANE162,SHELL163,SOLID164和LINK167。在第八章，将详细讲述刚性体。

每个实体单元，壳单元和梁单元的质量都平均分配给单元的节点。在壳单元和梁单元中，每个节点还将附加一个转动惯量；只采用一个单值，它的作用就是让质量围绕节点呈球形分布。第三章建模显式动态分析的第一步就是创建模型，使它能够表达进行分析的物理系统。用PREP7前处理器来建立模型。如果通过GUI路径进行分析的话，那么事先设置参考选项(MainMenu〉Preference)为“LS-DYNAExplicit”是很重要的。这样，菜单就被完全过滤成为显式动态的输入选项。(值得注意的是，Preference选项置为“LS-DYNAExplicit”并没有激活LS-DYNA求解。要做到这一点，就必须定义一个显式单元类型，例如，SHELL163。一旦设置好分析选项Preference，就可以像通常分析任何问题一样建立模型：•定义单元类型和实常数•定义材料模型•定义几何模型•划分网格•定义接触表面如果以前从未用过任何ANSYS产品，就应该参看一下《ANSYSBasicProcedureGuide》和《ANSYSModelingandMeshingGuide》,以了解ANSYS建模的一般过程。3.1定义单元类型和实常数在第二章中已简要地讲述了显式动态分析的单元类型，有关每种显式单元的详细描述可在《ANSYSElementReference》中找到，所以建议用户在确定使用哪种单元类型建模前仔细阅读一下有关内容。

一旦选择好能代表物理系统的单元类型，就可以用ET命令来定义单元类型（在GUI中：MainMenu〉Preprocessor〉ElementType）。在《ANSYSElementReference》中列出了与每种单元相对应的所有实常数。因此必须确定模型中所用每种单元的实常数，然后可以用R命令来定义实常数（在GUI中：MainMenu〉Preprocessor〉Realconstants）。3.2定义材料特性在显式动态分析中有很多可使用的材料类型，应该参看一下《ANSYSElementReference》，以确定特定的单元类型所用的材料模型。也可参看本手册的第七章，该章对所用材料模型作了详细的描述。一旦确定了所用的材料模型，就可以定义所有相关的特性（如第七章所描述的）。为了用批处理或命令流来定义这些，可以使用MP，TB和TBDATA命令（某些情况下，可用EDMP命令）。在GUI路径中，材料模型通过下列路径来定义：MainMenu〉Preprocessor〉MaterialProps>MaterialModels更详纟田的信息请参看7.1DefiningExplicitDynamicsMaterialModels.在选择了合适的材料模型之后，就可以定义模型所需的特性。定义和整体坐标不一致的各向正交异性模型，必须先用EDLCS命令定义局部坐标系。（菜单路径MainMenu〉Preprocessor〉MaterialProps〉LocalCS>CreateLocalCS）。对于一些材料模型，还需用EDCURVE命令定义与材料相关的数据曲线（例如，应力-应变曲线）。（EDCURVE的GUI路径表示法为MainMenu〉Preprocessor〉MaterialProps〉CurveOptions）。3.3定义几何模型建立几何模型的最简单方法就是采用ANSYS程序中的实体建模功能。关于实体建模功能的详细信息，请参看《ANSYSModelingandMeshingGuide》。对于简单模型（例如，仅线单元），就可以使用直接建模法。用这种方法，可以直接定义模型的节点和单元。详细信息请参看《ANSYSModelingandMeshingGuide》。3.4网格划分

建立实体模型后，就可以用节点和单元对其进行网格划分。《ANSYSModelingandMeshingGuide》中的GeneratingtheMesh详细描述了划分网格的步骤。如果对ANSYS网格不太熟悉，在划分显式动态模型之前应该先阅读该章内容，由于该章已详细讨论过了，在此只讨论以下几点：•定义单元属性•定义网格划分控制•生成网格定义单元属性，就是要事先指定单元类型，实常数和材料特性来用于下一步的网格划分。使用TYPE，MAT，REAL命令或菜单路径：MainMenu〉Preprocessor〉—Attributes—Define〉DefaultAttribsMainMenu〉Preprocessor〉—Modeling—create〉Elements〉ElemAttributes网格控制就是指定划分网格时单元的大小和形状。在ANSYS程序中有许多种控制方法（参考《ANSYSModelingandMeshingGuide》），在GUI中，菜单路径为：MainMenu〉Preprocessor〉MeshTool在选择网格控制时要注意以下几点：•尽量避免退化的实体和壳体单元（例如，三角型壳体和四面体实体），相对于四边形和六面体来说它们太刚硬，并且精度不高。・单元的大小尽量均匀（例如：避免产生相对较小的单元面积）单元大小之间差别很大可能会导致很小的时间步长，较长的运行时间。如果划分特殊的几何模型需要相对较小的单元，那么可以用质量比例来增大最小时间步长。（参看第十章，MassScaling）.•尽量不要使用SmartSizing方法进行单元控制（SMARTSIZE命令），因为它可能在网格中单元的大小上产生很大的差别。应使用ESIZE和相关的命令来控制单元大小。•尽量避免可能产生沙漏的坏单元形状。•当使用降阶单元算法或者单元可能遭受沙漏变形时，不要使用粗网格划分，•如果有沙漏现象的话，尽量在部分模型或整个模型中使用全积分单元。

给定网格控制后，就可以用相应的命令进行网格划分（例如：AMESH,VMESH）在GUI路径中，菜单路径为MainMenu〉Preprocessor〉-Meshing-Mesh（或用上面提到的MeshTool）.3.5定义接触面显式动态分析常常涉及到表面之间的接触。本手册的第六章ContactSurface讲述了接触类型和定义接触的步骤。这里只给出简要步骤。定义接触包括四个步骤：•定义能很好地表达物理模型的接触类型（EDCGEN命令）•定义接触表面（CM,EDPART,EDASMP和EDCGEN命令）•定义摩擦系数参数（EDCGEN）•为选定的接触类型定义附加数据（EDCGEN和EDCONTACT）如果不使用自动接触，那么可以用CM命令把表面上的节点组成为一个COMPONENT来定义接触面。一旦生成了COMPONENT,就可以用EDCGEN命令来定义所需表面间的接触（例如，节点部件）。也可以用EDCGEN命令来定义接触类型。对于单面接触算法，ANSYS/LS-DYNA定义接触体的外表面。注意一接触表面也可以用PART号或一个部件集合来识别（用EDASMP而不用节点部件。）部件和集合的定义将在这章的后面讨论。EDCGEN命令也可以用来指定摩擦系数参数以及与不同接触类型相对应的输入参数。另一个命令EDCONTACT可以定义多种接触控制例如接触穿透检测和滑移表面罚函数等。（菜单路径为:MainMenu〉Preprocessor〉LS-DYNAOptions〉Contact.）3.6建模的一般准则当创建显式动态模型时，应考虑下述基本准则：•在模型中使用刚性体来代表相对刚硬，没有屈服的部分。使用刚性体可以简化求解,缩短求解时间。•对材料性能使用一些符合实际的值。例如，不要用很高的不切实际的弹性模量来表达刚性体，对于壳单元不要使用不切实际的厚度值。

•考虑使用阻尼（EDDAMP命令）以消除结构响应中的不真实的振荡，详细信息参见《ANSYSCommandsReference》。•如果已经用常规ANSYS程序进行了二维的动态分析，应考虑把模型扩展为三维模型，并用ANSYS/LS-DYNA进行分析。那么就可以在较短的时间内得到更精确的结果。・注意ANSYS程序中的子模型和子结构特性不能在ANSYS/LS-DYNA中使用。PART的定义具有相同单元类型、实常数和材料号的一组单元被指定为一个PART并给定一个PART号。有许多ANSYS/LS-DYNA命令和PARTID有关。（例如，EDCGEN,EDDC,EDLOAD,EDDAMP,EDCRB和EDREAD），在ANSYS/LS-DYAN程序对模型进行求解时，自动生成PART号并且写入LS-DYNA输入文件Jobname.K，它可以用EDPART命令建立，更新，列表表示。从顺序编排所选单元建立PART号。如果在单元编排中改变单元类型、实常数和材料号中的任何一个，将给那组单元定义下一个PART号。例如，前十个单元有单元类型1、实常数2、材料号3,那么这10个单元将给定PART号为1。如果下面100个单元有单元类型1、实常数1、材料号2,那么这些单元将给定PART号为2,等等，依此类推。用EDPART命令中的CREATE/UPDATE/LIST选项可以创建、更新、列出PART表。如果执行SOLVE命令或执行带ANSYS/LSDYNA/BOTH选项的EDWRITE命令，则PART表将被建立并永久保存。如果PART表已经存在，则以后的命令只能修改这个表而不能覆盖它。关于EDWRITE命令的更多信息，参见《ANSYSCommandsReference》中有关这个命令的描述。EDPART,CREATE可以创建新的PART号。可以用EDPART,LIST命令列表表示这些PART。这个表显示了在建立或修改时PART的状态。（EDPART,CREATE）如果EDPART,CREATE命令重复使用，则PART表被覆盖。为了得到在修改或增加模型后的实际PART表，执行EDPART,UPDATE命令。它可以扩展已经存在的PART表而不用改变它的顺序，并且可以向已经存在的由相同的单元类型、实常数和材料号构成的部件中增加单元。单元类型、实常数和材料号相同的任何PART未被任何所选择单元参考时，则该PART不可用。很明显的它将在PART表的第五列上为零值。如果预先定义的PART相关命令与一个无用的PART有关，那么执行SOLVE命令或EDWRITE,ANSYS/LSDYNA/BOTH命令时将有一个警告。下例显示了15个单元的模型，带2种材料（MAT）,3种单元类型（TYPE）和3种实常数（REAL）。一个ELIST命令列出下列单元表：ELEMMATTYPRELESYTSHANODES

ELEMMATTYPRELESYTSHANODES1111012345678211101112131415161718311102122232425262728411103132333435363738511104142434445464748612305152535455565758712306162636465666768822307172737475767778922308182838485868788112320101102103104105106107108122320111112113114115116117118131320121122123124125126127128141320131132133134135136137138151320141142143144145146147148单元1到5有MAT=TYPE=REAL=1生成PART1,单元6到7有MAT=1,TYPE=2和REAL=3生成PART2。单元8到9有MAT=2,TYPE=2和REAL=3生成PART3。单元10有MAT=TYPE=REAL=2，但是在生成PART表后被删除了(EDPART,CREATE)。因此，只剩下了PART4,但是它被标识为无用值(在PART表中USED一列为0)。单元11到12有MAT=2,TYPE=3和REAL=2生成PART5。单元13到15有MAT=1,TYPE=3和REAL=2生成PART6。这一PART表如下所示(EDPART,LIST)：注意一对于这个例子，执行EDPART,UPDATE将不会影响PART表的状态。但是执行EDPART,CREATE将会创建一个仅有5个部件的PART表。PART5将变成PART4,PART6将变成PART5。这将使以前基于PART定义的载荷，接触特性等失效。

建议采用下列过程：如通常那样建模，直到PART的相关命令给定一个PART号。创建PART表（EDPART，CREATE）并显示它（EDPART，LIST）。从PART表中找到相应的PART号，用到PART的相关命令中。继续建模。如果单元或其属性改变，包括增加新单元，则要更新PART表（EDPART，UPDATE）。列出被修改的PART号，将用于进一步与PART有关的命令中。除了显示部件外，还可以用PARTSEL命令选择部件或绘制部件（通过PART号）。和其它ANSYS选项类似，这也有许多可能的选择类型，它们有特殊的标记：“S”代表选择，“R”代表重新选择，“A”代表另外的选择等等。（因为PARTSEL是一个命令集合，所以如果这个命令被包含在批处理文件中，那么在单一的引证中需关闭标识）。用PARTSEL命令的PMIN，PMAX和PINC域控制可选号的范围。例如。可以执行下列命令从已定义的PART表中选择PART2和4：PARTSEL，S，2，4，2当执行PARTSEL命令时，对每个存在的部件将自动生成一个单元组合。可以绘出这些组合，菜单路径为：UtilityMenu〉Plot〉Parts,命令为：PARTSEL，‘PLOT'。注意一在分析的其它部分中，不要试图使用PARTSEL命令创建的组合，因为如果PART表改变了它们也要重新定义。注意一不推荐选择或显示包含COMBI165单元的部件。CDWRITE命令把PART信息写入Johname.CDB文件。然后使用CDREAD命令将这一信息自动地读入ANSYS。但是，如果读入了两个或两个以上的Johname.CDB文件，那么最后一个Johname.CDB文件里的PART表将覆盖已经存在的表格。一般情况下，必须执行EDPART，CREATE命令来重新创建整个模型的PART表。但这将影响Johname.CDB文件里的所有与PART相关的命令。因此，用户可以组合模型，但不是那些与PART相关的输入，它们必须用最新创建的PART号修改。

在有限制的情况下，可以更新PART表（EDPART,UPDATE）。这就要求没有使用的MAT/TYPE/REAL的组合仅在表中出现一次。然而，有必要部分改变与PART相关的命令。3.7.1Part集合用EDASMP命令创建的部件集合，是由多个不同部件组成的实体。部件集合可以用来输入一些ANSYS/LS-DYNA命令。在定义包含多个部件的实体间的接触时，部件集合是非常有用的（参看本手册第六章）。使用EDASMP命令，在提供一个ID号后，可以在一个部件集合里定义16个部件。部件ID号必须比当前定义的最高部件号大（用EDPART，LIST确定最高部件号）。也可以分别用EDASMP，LIST和EDASMP，DELE选项显示和删除部件集合。例如，指定包含部件1，2,4的部件集合的ID号为10,可以执行下列命令：EDASMP，ADD，10，1，2，43.8自适应网格划分在金属成形和高速撞击分析中，物体要经历很大的塑性变形。单积分点显式单元，常用于大变形，但是在这种情况下，由于单元纵横比不合适可能给出不精确的结果。为了解决这一问题，ANSYS/LS-DYNA程序可以在分析过程中自动重新划分表面来改善求解精度。这一功能，即自适应网格划分，由EDADAPT和EDCADAPT命令控制。EDADAPT命令在一个指定的PART内激活自适应网格划分。（用EDPART命令创建或显示有效PARTIDs），例如，为了给PART1打开自适应网格划分，可以执行下列命令：EDADAPT，1，ON注意一自适应网格划分功能仅对包含SHELL163单元的部件有效。当此项功能打开时，分析中该部件的网格将自动重新生成。从而保证在整个变形过程中有合适的单元纵横比。自适应网格划分一般应用在大变形分析例如金属变形中（调节网格最典型的应用是板料）。在一个模型中要在多个部件上应用此功能，必须对每个不同的PARTID执行EDADAPT命令。缺省时，该功能是关闭的。在指定哪些部件重新划分后，必须用EDCADAPT命令定义网格划分参数。采用EDADAPT命令定义需要网格划分的所有PARTID号，用EDCADAPT命令对其设置控制选项。EDCADAPT命令控制的参数如下所示：

・Frequency(FREQ)-调节自适应网格划分的时间间隔。例如，假设FREQ设置为0.01,如果单元变形超过指定的角度容差，则其将每隔0.01秒被重新划分一次(假设时间单位为秒)。因为FREQ的缺省值为0.0，所以在分析中应用自适应网格划分时必须指定此项。・AngleTolerance(TOL)-对于自适应网格划分(缺省值为le31)有一个自适应角度公差。TOL域控制着单元间的纵横比，它对保证结果的精度是非常重要的，如果单元之间的相对角度超过了指定的TOL值，单元将会被重新划分。•AdaptivityOption(OPT)-对于自适应网格划分有两个不同的选项。对于OPT=1,和指定的TOL值相比较的角度变化只是根据初始网格形状计算的。对于OPT=2,和指定的TOL值相比较的角度变化是根据前一次重新划分的网格计算的。•MeshRefinementLevels(MAXLVL)-MAXLVL域控制着整个分析中单元重新划分的次数。对于一个初始单元，MAXLVL=1可以创建一个附加单元,MAXLVL=2允许增加到4个单元，MAXLVL=3允许增加到16个单元。高MAXLVL会得到更精确的结果，但也会明显增加模型规模。•RemeshingBirthandDeathTimes(BTIMEandDTIME)-自适应网格划分的生死时间控制着该功能在分析过程中的激活或关闭。例如，如果设置BTIME=.01和DTIME=.1,那么分析中只在.01和.1秒间进行重新网格划分(假设时间单位为秒)。•IntervalofRemeshingCurve(LCID:loadcurveID,载荷曲线参考号)-数据曲线把重新划分网格的时间间隔定义为时间的函数。数据曲线的横坐标为时间，而纵坐标为变化的时间间隔。如果这个选项不为0,那么它将代替适应频率(FREQ)。但是要注意，开始第一个自适应性循环仍需要非零FREQ值。•MinimumElementSize(ADPSIZE)-根据单元边长设定的最小单元尺寸。如果不定义此参数，边长的限制将被忽略。•OneorTwoPassOption(ADPASS)-如果ADPASS=0,将使用双通道自适应划分，在重新划分网格后将重复这一计算(缺省值)。如果ADPASS=1,则使用单通道自适应划分，而计算不再重复。关于这两个选项的图形表示，请参看《ANSYS/LS-DYNATheoreticalManual》图30.9(a)和30.9(b)。•UniformRefinementLevelFlag(IREFLG)-值为1,2,3等，分别允许4,16,64等划分等级。对于每个初始单元都分别生成统一的单元。—使用—使用component或PARTIDs和数组参数・PenetrationFLAG(ADPENE)-根据ADPENE是正(到达)的还是负(穿透)的，当接触表面到达或穿透工具表面时,程序将依据这个值调整网格。自适应细化主要依据加工曲线。如果ADPENE是正的，细化一般发生在接触之前；因此，可能用单通道划分就可以了。(ADPASS=1)•ShellThicknessLevel(ADPTH)-绝对薄壳厚度标准，低于该值自适应划分开始。这个选项仅在自适应角度公差不为零的情况下有用。如果期望不改变角度，激活基于厚度的自适应重新划分功能，那么可把TOL设为较大的角度。(如果ADPTH=O.O,不使用这个选项。)•MaximumElementLimit(MAXEL)-自适应结束的最大单元数。如果超过了此值，自适应将被中止。对于大多数问题，不应该用自适应网格划分作初始分析。如果分析结果出现扭曲的网格,或结果不正确，那么再应用此项重新进行分析。当LS-DYNA分析中因为“negativevolumeelement”的错误终止时，也可用此项分析。当打开自适应网格划分选项时,在求解过程中模型内的单元数将发生变化。在一个调节循环结束后，网格将被更新，并且生成一个扩展名为RSnn的新结果文件，这里nn为自适应网格标准。(在由FREQ指定的每个时间增量或LCID指定的时间间隔，都会发生自适应循环。)例如，重新划分两次网格的模型将产生两个结果文件，Johname.RSOl和Johname.RS02.有关自适应划分结果后处理的详细信息，请参看第十二章，Postprocessing.注意一既使每次循环网格不发生变化，也会产生一个新的结果文件。每次循环都会产生一个扩展名为Hinn的时间历程文件，尽管这些文件可能不如RSnn文件有用。此外，LS-DYNA还创建了一系列源文件名为“adapt”的文件。因此，在激活调节网格划分时，不要把“adapt”作为你的工作名。第四章加载当模型建好后，下一步就是给结构加载为求解作准备，为了能正确地模拟结构的响应,就必须定义与指定时间间隔相对应的载荷，本章将讲述以下有关加载的几个方面：•一般加载选项

—如何施加、删除、显示一般载荷[EDLOAD]—如何绘制载荷曲线[EDPL]—如何定义数据曲线[EDCURVE]—如何显示或隐藏载荷标记•约束和初始条件—如何在ANSYS/LS-DYNA中施加约束]D,EDNROT]—如何定义滑移和循环对称平面[EDBOUND]—如何定义混合型约束[EDCNSTR]—如何定义焊点[EDWELD]—如何给模型施加初始速度[EDVEL,EDPVEL]•耦合和约束方程—自由度耦合[CP]—自由度间的约束方程[CE]•非反射边界[EDNB]•温度载荷•动力松弛4.1一般载荷选项与许多隐式分析不同的是，显示分析中的所有载荷必须与时间有关。因此，在ANSYS/LS-DYNA中，许多标准的ANSYS命令都是无效的。在ANSYS/LS-DYNA中，尤其不能使用F,SF,BF系列命令，因为它们只能定义与时间无关的载荷。此外，D命令只能定义节点约束。基于上述原因，在ANSYS/LS-DYNA中用一对数组参数定义载荷（一个用来定义时间，另一个定义载荷）。注意一虽然节点加速度(A,A,A)和节点速度(V,V,V)以自由度出xyzxyz现，但它们不是物理自由度，不能使用D命令约束。要采用EDLOAD命令给这些节点施加载荷。在ANSYS/LS-DYNA中，所有载荷都是在一个载荷步内施加的。这和隐式分析有很大的不同，它在多个载荷步内施加载荷。在ANSYS/LS-DYNA中，对于一些特定的载荷，也可以用EDLOAD命令指定何时施加(birthtime)、何时去除(deathtime)。请参考EDLOAD命令中的BirthTime,DeathTime和CID，检验birth/deathtime的适用性。给模型施加载荷，需遵循以下步骤：•把模型中受载的那部分定义成Component(或PART，用于刚体)•定义包含时间间隔和载荷数值的数组参数•定义载荷曲线•如果不是在整体坐标系中加载，需要用EDLCS命令定义载荷方向•模型加载4.1.1组元除给刚性体加载外，显式分析中所有载荷都施加到Component上。因此，第一步就是把模型中受载的那部分组合成Component,每个Component应由模型中承受同样载荷的部分组成。并且可以通过材料本构、模型中位置、预期状态等联系在一起。例如，想要分析一个棒球撞击到墙上的结果，可以定义球上的节点为一个Component,球棒上的节点为另一个Component,墙上的节点为第三个Component。可以定义任意多个Component,然后给每一个Component加载，Component必须由节点或单元组成(只有当施加压力载荷时Component才由单元组成)定义Component时，首先要选择Component中想要包含的部分，然后使用GUI中的下列菜单路径来定义Component(CM):UtilityMenu〉Select〉Comp/Assembly〉CreateComponent关于Component的详细信息，参见《ANSYSAnalysisGuide》中的SelectingandComponents。对于刚性体，载荷一般都施加到PART号上，而不是Component上。这是因为当使用命令EDMP,RIGID,MAT定义刚性体时已经包含有一系列节点和单元。

4.1.2数组参数值得注意的是，显式动态分析中，载荷在一些特定时间间隔内施加到结构上，这样就可以观察在特定时间段里施加载荷后模型的瞬态行为。因此，不仅需要定义施加载荷的类型（FX,FY,FZ,ROTX,ROTY,ROTZ等），同时也需要定义载荷施加到结构上的时间间隔值。时间间隔值和其相对应的载荷值组合在一起定义为数组参数。这些参数必须定义为两部分，第一部分为时间间隔值，第二部分为载荷值，载荷应与时间间隔相对应。在｝UI中，采用菜单路径来定义一个数组参数：UtilityMenu〉Parameters〉ArrayParameters〉Define/Edit.有关如何定义数组参数的详细信息，参看《ANSYSAPDLProgrammer'sGuide》中的ArrayParameters。注意一可以通过线性插值得到中间时间点的载荷值。指定时间范围外的载荷值不能由程序插入。因此，必须保证载荷时间范围至少等于求解时间。否则，由于过早的去除载荷求解结束时的结果将会变为无效值。一旦定义了能代表载荷的数组参数，这些载荷与时间有关，就可以直接用EDLOAD命令输入参数定义载荷。或者，可以用EDCURVE命令输入参数来定义载荷曲线。相对应的载荷曲线ID可用EDLOAD命令输入。为说明数组参数的应用，考虑前面提过的棒球例子。想要考察从击球瞬间到击球一秒钟后球的变形。假定位移是时间的已知函数，球刚碰撞时的初始速度为1600in/sec（91MPH）。首先需要定义一些节点组元，用它们来定义载荷和接触面。建立一个Component,它包括球上的所有节点，取名为ball。对这个Component施加一个1600in/sec的初始速度（EDVEL）,（本章后面再讨论）。然后再创建第二个组元，它仅包括球表面上的所有节点，取名为ballsurf.这个Component以后将用于接触定义。还需定义第三个Component,它包括球棒表面上的节点，取名为batsurf,接触算法将在第六章接触表面中讨论。nsel,s,node,…！选择组成球的所有节点cm,ball,node!定义Component为ballnsel,s,node,…！选择组成球表面的所有节点cm,ballsurf,node!定义Component为ballsurf

nsel,s,node,…！选择球棒表面的所有节点cm,batsurf,node!定义Component为batsurfnsel,all现在选择球棒基座上的节点(bathand)，对它施加一个位移时间载荷曲线。定义一名字为time的数组来存储时间值。使用时间值要与模型的所有载荷、维数以及材料性质相协调。其次，定义一个数组，包括Componentbathand相应的X位移值，命名为Xdisp.同样定义数组Ydisp和Zdisp,以包含球棒相应的Y位移和Z位移值。nsel,s,node,…！选择球棒基座上的节点cm,bathand,node!定义Component为bathandnsel,all*dim,time,,4!定义数组参数为time的维数*dim,xdisp,,4!定义数组参数为xdisp的维数*dim,ydisp,,4!定义数组参数为ydisp的维数*dim,zdisp,,4!定义数组参数为zdisp的维数time(l)=0,0.25,0.5,0.75,l!给定位移的时间值xdisp(l)=0,-l,-2,-l,3!球棒的X位移ydisp(l)=0,l,2,3,4!球棒的Y位移zdisp(l)=0,3,6,8,9!球棒的Z位移给出的例子实际是一个比较复杂的现象的简化模型。更加精确地模拟，应该定义额外的位移位置(和相应的载荷曲线)以更好的模拟球的真实运动。另外，球棒上各节点的初始速度是不同的。最后，球是一些不同材料和材料模型的复合体。4.1.3施加载荷

一旦定义好Component和数组参数，就可以给建立的模型加载（EDLOAD命令）。在GUI中，可以选择下列菜单路径：MainMenu〉Solution〉LoadingOptions〉SpecifyLoads.可以选择增加载荷（EDLOAD中ADD标号），如力、力矩、节点加速度、表面压力等，所有的载荷都用EDLOAD命令施加到整体笛卡尔坐标系上，此为其缺省值。GUI的对话框将列出所有有效载荷，以及先定义好的Component和数组参数。只需简单地选择所需的载荷标号，以及Component（刚体的PART号）和数组参数集合（预先定义的载荷曲线号）。值得注意的是，如下表所示的那样，并不是所有的载荷标记都适用于所有Component或PART号。以下的载荷标记只适用于节点Component：力：FX,FY,FZ力矩：MX,MY，MZ位移：UX，UY，UZ转角：ROTX，ROTY，ROTZ速度：VX，VY，VZ节点加速度：AX，AY，AZ体加速度：ACLX，ACLY，ACLZ角速度：OMGX，OMGY，OMGZ温度：TEMP注一虽然V（X,Y，Z）和A（X，Y，Z）作为DOFs出现，它们实际上不是物理DOFs。然而，这些量是被计算作为DOF求解和存储以用于后处理。以下的标记只适用于刚性体（部件号）力：RBFX，RBFY，RBFZ力矩：RBMX，RBMY，RBMZ

位移：RBUX,RBUY,RBUZ转角：RBRX,RBRY,RBRZ速度：RBVX,RBVY,RBVZ角速度：RBOX,RBOY,RBOZ以下标记只适用于单元Component:压力：PRESS回到前面的棒球例子，需由time和x/y/zdisp数组参数用EDLOAD命令自动生成位移与时间的载荷曲线。edload,add,ux,,bathand,time,xdisp!球棒的x位移edload,add,uy,,bathand,time,ydisp!球棒的y位移edload,add,uz,,bathand,time,zdisp!球棒的z位移另一种方法，很刚硬的棒球可以模拟为刚体，以简化所需输入数据和减少CPU时间。对于这种情况，对应刚体载荷可以施加有效的PART号（不是节点Component）.edload,add,rbux,,2,time,xdisp!x位移（如果球棒的PART号为2）edload,add,rbuy,,2,time,xdisp!y位移（如果球棒的PART号为2）edload,add,rbuz,,2,time,xdisp!z位移（如果球棒的PART号为2）如果已知球棒对球碰撞的压力是时间的函数，则可以不定义任何接触表面就能实现加载。取而代之用带“PRESS”标号的EDLOAD命令对包含球表面单元的组元Component加载。edload,add,press,l,cover,battime,batload在这种情况下，在面号1上对单元组元cover加载（在KEY域输入面号），其载荷曲线是由battime和batload数组参数产生的。特别要提出的是，在上例中，仅仅压力载荷采用EDLOAD命令的LKEY域。LOADkeys（1,2,3等）与面载荷有关，在《ANSYSElementsReference》中每种单元类型在输入数

据表“surfaceloads”中列出。对于许多非压力的载荷，可以使用KEY域定义坐标系标号CID。载荷将作用于用EDLCS命令定义的CID的方向上，或者说，如果没有指定CID,就作用于整体坐标系上。详细信息请参看§4.1.5DefiningLoadsinaLocalCoordinateSystem。注意一为避免某些平台的定时问题，通常是加一个小的时间值（如1.0X10-6）到时间数组的最后一项•例如，取代3.0的值,数组的最后一项为下列值：timeint（l）=0,l,2,3.00001增加这样小的缓冲，不会影响结果的精度。除了施加载荷外，也可以用EDLOAD，LIST和EDLOAD，DELE来显示和删除载荷。用EDFPLOT命令也可以显示单元上的载荷标号，并且可以用EDPL命令绘制载荷曲线。4.1.4数据曲线用EDCURVE命令定义的数据曲线广泛应用于ANSYS/LS-DYNA中。它们可以用来定义与显式动态材料模型有关的材料数据曲线（例如，应力-应变）和载荷数据曲线（力-偏转）。它还可以定义时间载荷曲线（力，位移，速度等）。这些载荷曲线可以用EDLOAD命令输入。使用材料模型数据曲线某些材料模型（例如，TB，PLAW或TB，HONEY）要求指定材料特性数据，它们可能是有效应变速率、塑性应变或体积应变的函数。对于这些数据，在用数据表［TBDATA］命令定义材料特性之前，需用EDCURVE命令定义特性曲线。在刚性体和压延筋接触问题中，也采用数据曲线来定义变形特性。与对Component加载相类似，数据曲线组合成数组参数，然后与特定的曲线参考号相联系，这个参考号可用于指定的材料模型（PLAW，HONEY，等）或接触类型（RNTR，ROTR）和压延筋。定义数据曲线可以分为以下几个步骤：定义一个包含材料或摩擦力特性横坐标的数组参数（例如，有效塑性应变,有效应变率,位移等）定义第二个数组参数,包含材料特性或摩擦力的纵坐标值。（例如，初始屈服应力，弹性模量，力等）。3•定义数据曲线（EDCURVE）。选择一个数据曲线ID号，产生数据表［TBDATA］时将采用这个数据曲线ID号来将这些数据与特定的材料特性相联系。

定义这些参数后，在GUI中采用下列路径定义数据曲线:MainMenu〉Preprocessor〉MaterialProps〉CurveOptionsMainMenu〉Solution〉LoadingProps>CurveOptions可以用EDCURVE,LIST显示数据曲线，用EDCURVE,PLOT绘制曲线，用EDCURVE,DELETE删除曲线。下例讲述了怎样用曲线来定义钢的分段线性塑性曲线（TB，PLAW，，，，8）：!“3”是任选的材料参考号（MAT）mp,ex,3,30.0e6!弹性模量（psi）mp,dens,3,7.33e-4!质量密度（lbf-sec2/in4）mp,nuxy,3,0.30!泊松比注：首先将工程应变与工程应力关系的数据转换成真实应力与真实应变关系的数据，然后，从总真应变中减去弹性真应变，求出塑性真应变。它与总真应变一起用于LS-DYNA的*MAT-PIECEWISE-LINEAR-PLASTICITY材料模型#24。!Stress-StrainDatausedwithPiecewiseLinearPlasticity(PowerLaw8):!TotalTotalTotalTotalElasticPlastic!Stress/Eng.Eng.TrueTrueTrueTrue!StrainStressStrainStressStrainStrainStrain!Point(psi)(in/in)(psi)(in/in)(in/in)(in/in)

!!100.000000.00000.00000.0000!260,0000.002060,1200.00200.00200.0000!377,5000.032580,0200.03200.00270.0293!483,3000.083590,2600.08020.00300.0772!598,0000.1735115,0000.16000.00380.1562!698,3000.2710124,9400.23980.00420.2356!776,4001.2255170,0300.80000.00570.7943!注：应力/应变曲线的第一点未输入，起始于第二点（横坐标=屈服应力）。要与用*SET命令定义的数组空间相协调。*dim,strn,,6!定义有效塑性真应变数据数组*dim,strs,,6!定义有效总真应力数据数组strn（l）=0.0,0.0293,0.0772,0.1562,0.2356,0.7493!应变（in/in）strs（l）=60120.,80020.,90260.,115000.,124940.,170030.!应力（psi）edcurve,add,l,strn,strs!曲线#1,纵坐标二应变，横坐标二应力tb,plaw,3,,,8tbdata,6,l!为应力/应变数据采用载荷曲线#1注：如果需要，可以定义塑性失效应变。其次，还可以用给定必要的应变率参数或载荷曲线定义应变率对屈服应力的影响。请参见本手册的第七章，对这个材料模型有完整的描述。使用载荷数据曲线

除了用于特定材料模型外，数据曲线还可以用来定义与时间有关的载荷。除了第一个数组参数必须包括时间值，第二个数组参数必须包括相应的载荷值外，定义载荷曲线的步骤和上述材料数据曲线一样。在用EDCURVE命令定义载荷曲线后，就可以用EDLOAD命令输入相应的载荷曲线参考号（LCID）。下例讲述了用EDLOAD命令定义载荷曲线的3个步骤：!步骤1：定义数组参数*dim,time,,5time（1）=0,.025,.05,.075,.1*dim,yforce,,5yforce（1）=0,100,200,300,400!步骤2:定义载荷曲线和相应的LCID（#11）edcurve,add,11,time,yforce!步骤3:用EDLOAD命令指定LCIDedload,add,fy,,comp,,,,11,1.0如上所述的那样，在EDLOAD命令的第九个域定义LCID。值得注意的是，如果EDLOAD命令用一个LCID，数组参数不能用来定义载荷。当特定载荷曲线用于多个组元或载荷标记时，在EDLAOD命令中使用LCID而不用数组参数是非常有用的。4.1.5在局部坐标系中定义载荷ANSYS/LS-DYNA中，可以在已定义的坐标系的任意坐标方向指定运动为一个组元或部件。这可以用EDLOAD命令的KEY域结合局部坐标系来获得（EDLCS命令）。一旦用EDLCS命令定义了局部坐标系后，（EDLCS，ADD，CID，X1，Y1，Z1，X2,Y2,Z2），就可以在EDLOAD命令中的KEY域用局部坐标系ID（CID）定义载荷的作用方向。如果没有指定CID，载荷就会作用在整体笛卡尔坐标系中。一些载荷类型不支持局部坐标系;在EDLOAD命令注释部分，参看BirthTime,DeathTime和CIDSupportTable，可以得到更多的信息。4.1.6指定Birth和Death时间对于每个载荷定义，可以用EDLOAD命令的BTIME和DTIME域来定义激活或杀死时间（EDLOAD，ADD，Lab,KEY,Cname,Par1,PHASE,LCID,SCALE,BTIME,DTIME）。在分析中使用

这些选项就可以在任意时刻激活载荷，随后杀死。在多阶段成形过程中这些选项是非常有用的，这一过程需要连续施加多个载荷。一些载荷类型不支持局部坐标系；在EDLOAD命令注释部分，参看BirthTime,DeathTime和CIDSupportTable，可以得到更多的信息。4.2约束和初始条件在开始求解之前，需要给模型施加约束。另外，还可能给运动物体设定初始速度。4.2.1约束与ANSYS（隐式）不同，ANSYS/LS-DYNA区分零约束与非零约束。非零约束如载荷那样处理（伴随着载荷曲线，见本章前面的讨论）。仅零约束才能使用D命令，也就是说，给定的值必须总是零。其它值无效，D命令仅用于固定模型的某些部分。可以用零约束来实现对称/反对称边界条件。用户可以用EDNROT命令在旋转节点坐标系中施加零载荷。但首先需用EDLCS命令定义局部坐标系。当模拟几何体的小对称部分时，需定义滑移或循环对称。可以使用EDBOUND命令来定义滑移或循环对称的对称边界面。可以用节点组元确定边界或方向矢量来定义法向（滑移对称）或旋转轴（循环对称）。在ANSYS/LS-DYNA中，可以用EDCNSTR命令模拟其它类型的约束。可用的约束类型有附加节点设置（ENS），节点刚性体（NRB），薄壳边界到实体（STS），以及铆接（RIVET）。在GUI中，可以通过下列路径施加这些约束：MainMenu〉Solution〉Constraints〉Apply〉AdditionalNodal附加节点设置约束类型（EDCNSTR,ADD,ENS）允许在一个已经存在的刚性体上增加节点（通过节点组元），这个刚性体是用EDMP命令定义的，该节点组元不能和其它任何刚体连接在一起。在刚体上附加的节点可以放置在模型的任何地方,并且可放置在初始刚体外。ENS选项有许多应用，包括在两个刚体结合的地方设置节点，定义施加节点载荷的节点，还有在指定的位置定义集中质量。不像用EDMP命令定义的典型刚体一样，用EDCNSTR,ADD,NRB命令定义的节点刚体不是和一个部件号相联系，而是与一个节点组元有关。当模拟刚性连接时（焊接），NRB选项是非常有用的。对于一个刚性连接，不同柔性组元（有不同的MATIDs）的部分作用在一起形成一个刚体。所以说很难用一个单一的MATID（和相应的部件号）来定义这种类型的刚体。

但是，用一个节点刚体很容易定义刚性连接。因为节点刚体不是和一个部件号相联系，所以使用刚体（例如用EDLOAD命令施加的载荷）的其他选项不能用于节点刚体。薄壳到实体边界选项（EDCNSTR，ADD，STS）把实体单元区域和薄壳单元区域固连起来，如图4-1所示，ConstrainedShelltoSolid,可将一个单独壳节点固连到最多9个实体节点上，这些实体节点定义一个fiber'矢量。定义fiber'矢量的实体单元节点在整个分析中保持线性但是在fiber方向上保持相对移动。薄壳节点必须和某个沿fiber方向的实体单元节点位置重合。ConstrainedSluelltoSolid图4-1与用EDWELD命令定义的焊接点相似，RIVET（EDCNSTR，ADD，RIVET）选项在两个不同位置的节点间定义了一个无质量刚性约束。和焊接点不同的是，一个rivet不能定义破坏。当定义了一个rivet时，节点间的距离将在模拟中的任何运动中保持不变。由rivet连接的节点不能是模型中其它约束的一部分。4.2.2焊接在显动态分析中，模拟被焊接在一起的组元是很普遍的。在部件通常是通过焊接装配在一起的汽车应用中，显得尤其普遍。这种情况下，在ANSYS/LS-DYNA中使用EDWELD命令来模拟焊接约束。可以模拟两种不同类型的焊接：无质量焊接点和一般焊接。用EDWELD命令连接的节点不能用其它方式约束。对于一个无质量焊接，必须指定两个不同位置的节点。也可以用EDWELD命令输入失效参数在焊点内定义失效。失效依据下列关系：

在GUI中，用下列菜单路径来定义一个无质量焊点：MainMenu〉Preprocessor〉LS-DYNAOptions〉Spotweld〉MasslessSpotwld一般焊接用于模拟两个部件的长焊接截面。对于一般焊接，你必须指定一个有效的节点组元，可以使用重合节点。但是，如果使用了重合节点，就必须定义输出数据所用的局部坐标系。一般焊接内的失效也可以用EDWELD命令的失效参数来定义，使用与上述等式相同的关系式。在GUI中，采用下列菜单路径来定义一般焊接：MainMenu〉Preprocessor〉LS—DYNAOptions〉Spotweld〉GenrlizdSpotwld4.2.3初始速度在瞬动态模拟中，经常需要定义初始条件。在ANSYS/LS-DYNA中，可以用EDVEL和EDPVEL命令定义初始速度。可以使用这些命令对各种实体施加线速度和角速度。用EDVEL对节点组元或单个节点施加速度；用EDPVEL给部件或部件集合施加速度。EDVEL和EDPVEL提供了两种定义角速度的方法，Option=VGEN和Option=VELO.VGEN方法围绕指定轴对一个实体（节点组元，部件等）施加刚体