Release_Notes_LS-DYNA_971_r7.1.1_rev0

利用LS-DYNA进行成形仿真的输入控制参数

利用LS-DYNA进行成形仿真的输入控制参数Translated by SunnyWinterLS-DYNA已广泛用于汽车碰撞分析。

默认的输入参数一般能给出有效，精确的碰撞模拟结果。

但是，这些默认值对于成形仿真分析并不一定理想。

下面是一个标准的金属成形过程。

为及时参考，推荐输入参数用黑体字标识，并包含在盒状关键字输入框中。

模型明确要求的数据，如终止时间等参数，输入黑体的0值。

一般问题设定在显式成形仿真中，利用质量比例缩放和（或者）人为的高工具速度，运行时间可以大大缩减。

这两种方法都会引入人为的动力学影响，因此必须将其减小到在工程意义上合理的水平。

一个单独的描述人为动力影响的参数是：工具每运动1毫米所采用的显式时间步进值（或周期）数目。

当成形过程允许大的无限的板料运动，比如冲击成形，需要更多的毫米周期数。

当板料被压边圈和冲模支撑较强的约束住时，较少的毫米周期数是必要的。

对大多数的仿真来说，100到1000之间的毫米周期数能产生合理的结果。

如果可能，或者有必要重复一个仿真，可利用两个不同的毫米周期值并比较分析结果去估计其对人为动力学影响的敏感性。

推荐选择的一个最大工具速度是2.0mm/ms,起始和结束速度为0。

可以使用简单的梯形速度轮廓（如图1）。

利用大的时间缩放步参数dt2ms获得要求的毫米周期数，可参考下面的公式：时间步大小=1.0/(最大工具速度*毫米周期数)工具速度，时间步大小和结束时间必须在协调的参照系中选择。

如果所有的工具运动给定，可用下面的步骤设置模拟参数：已知：工具全部行程（mm）:D最大工具速度（mm/ms）:2.0速度轮廓： 2.0毫秒上升和2.0毫秒下降的梯形（如图1）选择：毫米周期数：ncpm计算：结束时间（ms）: T=2.0+D/2速度数据点：（0.0，0.0）（2.0，2.0）（T-2.0,2.0） (T,0.0)时间步大小（ms）: dt2ms=1/(2*ncpm)上面的运算提供速度轮廓数据点用于下面的工具运动部分。

lsdyna 安装

安装了64位Linux系统(Ubuntu11.04 Server)，无Xwindow界面，决定尝试下Linux下的dyna计算，下面把过程叙述一遍：1.用另一台电脑从lstc官方网站下载64位Linux求解器：ls971_s_R5_1_1_intel64_RedHat54.gz，其他版本也可以，反正里面没Ubuntu 版的；2.为方便操作，将解压出来的求解器改名为ls971，通过U盘拷贝到Linux机器下的/home/username(你的帐户名)/lsdyna/ 文件夹下；3.在/home/username/下创建keyfiles文件夹，通过U盘拷贝测试用的testcase.k文件到keyfiles下，这一步可与上一步一起完成；4.进入lsdyna文件夹，输入./ls971 info，./ 表示执行程序，若程序不在当前文件夹下，则程序名需要带上路径(见第7步)，info参数用于生成注册文件lstc.log，但如果你安装的Linux也是Ubuntu，将会报错，错误为：缺少libnuma.so.1，网上可下载，放入/usr/lib64/文件夹下即可；5.重新执行./ls971 info，同文件夹下生成lstc.log，联系LSTC公司或者其代理商获取LICNESE文件6.讲得到的LICENSE文件LSTC_FILE放到/usr/local/lstc/文件夹下，【 Linux 公社】并设置全局环境变量LSTC_FILE="/usr/local/lstc/LSTC_FILE"7.进入/home/username/keyfiles/文件夹，执行算例: ./home/username/lsdyna/ls971 i=testcase.k 当然还可以设置memory、cpu等参数，与win下相同；计算结果文件会输出到当前文件夹下，应该可以通过指定o=参数更改8.程序正常执行，开始计算。

LS-DYNA问答总结

1 如何处理LS-DYNA中的退化单元？在网格划分过程中，我们常遇到退化单元，如果不对它进行一定的处理，可能会对求解产生不稳定的影响。

在LS-DYNA中，同一Part ID 下既有四面体，五面体和六面体，则四面体，五面体既为退化单元，节点排列分别为N1，N2，N3，N4，N4，N4，N4，N4和N1，N2，N3，N4，N5，N5，N6，N6。

这样退化四面体单元中节点4有5倍于节点1-3的质量，而引起求解的困难。

其实在LS-DYNA的单元公式中，类型10和15分别为四面体和五面体单元，比退化单元更稳定。

所以为网格划分的方便起见，我们还是在同一Part ID下划分网格，通过*CONTROL_SOLID关键字来自动把退化单元处理成类型10和15的四面体和五面体单元。

2 LS-DYNA中对于单元过度翘曲的情况有何处理方法有两种方法：1. 采用默认B-T算法，同时利用*control_shell控制字设置参数BWC＝1，激活翘曲刚度选项；2. 采用含有翘曲刚度控制的单元算法，第10号算法。

该算法是针对单元翘曲而开发的算法，处理这种情况能够很好的保证求解的精度。

除了上述方法外，在计算时要注意控制沙漏，确保求解稳定。

3 在ANSYS计算过程中结果文件大于8GB时计算自动中断，如何解决这个问题？解决超大结果文件的方案:1. 将不同时间段内的结果分别写入一序列的结果记录文件；2. 使用/assign命令和重启动技术；3. ANSYS采用向指定结果记录文件追加当前计算结果数据方式使用/assign指定的文件，所以要求指定的结果记录文件都是新创建的文件，否则造成结果文件记录内容重复或混乱。

特别是，反复运行相同分析命令流时，在重复运行命令流文件之前一定要删除以前生成的结果文件序列。

具体操作方法和过程参见下列命令流文件的演示。

4关于梁、壳单元应力结果输出的说明问题：怎样显示梁单元径向和轴向的应力分布图（我作的梁单元结果只有变形图DOF SOLUTIN –Translation，但是没有stress等值线图，只有一种颜色）和壳单元厚度方向的应力、变形图（我们只能显示一层应力、变形，不知道是上下表层或中间层的结果）。

LS-DYNA使用指南中文版本

第一章引言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,就需通过以下两本手册了解基本的分析求解过程：·ANSYS Basic Analysis Guide·ANSYS Modeling and Meshing Guide使用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-DYNA ASCII输出文件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命令按字母顺序排列的详细资料(包括每条命令的特定路径)，请参阅《ANSYS Commands Reference》。

LS-DYNA使用指南中文版本

第一章引言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,就需通过以下两本手册了解基本的分析求解过程：·ANSYS Basic Analysis Guide·ANSYS Modeling and Meshing Guide使用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-DYNA ASCII输出文件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命令按字母顺序排列的详细资料(包括每条命令的特定路径)，请参阅《ANSYS Commands Reference》。

LS-DYNA经典命令集详解

/EXIT,Slab,Fname,Ext,Dir Slab=ALL 保存所有资料Slab=NOSAVE所有更改资料不保存Slab=MODEL保存实体模型，有限元模型，负载的资料（系统默认）例：/EXIT，ALL--------------------------------------------------------/FILNAM，Fname Fname=工作文件名称，不要扩展名例：/FILNAM，Sanpangzi--------------------------------------------------------/SAVE,Fname,Ext,Dir 保存目前所有的Datebase资料，即更新Jobname.db--------------------------------------------------------/RESUME,Fname,Ext,Dir,NOPAR 回到最后SAVE时的Datebase状态--------------------------------------------------------/CLEAR 清除所有Datebase资料--------------------------------------------------------LOCAL，KCN，KCS，XC，YC，ZC，THXY，THYZ，THZX，PAR1，PAR2 定义区域坐标系统KCN 区域坐标系统代号，大于10的任何号码KCS=0，1，2 0=笛卡儿坐标 1=圆柱坐标 2=球面坐标XC，YC，ZC 该区域坐标原点与整体坐标原点的关系THXY，THYZ，THZX 该区域坐标与整体坐标XYZ轴的关系例：LOCAL，11，1，1，1，0--------------------------------------------------------CSYS，0，1，2 声明当前坐标系统例：CSYS，0--------------------------------------------------------/UNITS，LABEL 声明系统分析时所用的单位LABEL=SI (米，千克，秒)LABEL=CGS (厘米，克，秒)LABEL=BFT (英尺)LABEL=BIN (英寸)例：LABEL，SI--------------------------------------------------------/PREP7 进入通用前处理器--------------------------------------------------------N，NODE，X，Y，Z，THXY，THYZ，THZX 定义节点NODE 节点号码 X，Y，Z 节点在当前坐标系中位置例：N，1，2，3，4--------------------------------------------------------NDELE，NODE1，NODE2，NINC 删除已建立的节点NODE1，NODE2 删除从NODE1到NODE2的节点，如1到100NINC 间隔号码，1为1到100全删，2为1，3，5 (99)例：NDELE，1，100，2--------------------------------------------------------NPLOT，KNUM 将节点显示在图形窗口中KNUM=0不显示节点号码KNUM=1显示节点号码--------------------------------------------------------NLIST，NODE1，NODE2，NINC 将节点资料列在窗口中例：NLIST--------------------------------------------------------NGEN，ITIME，INC，NODE1，NODE2，NINC，DX，DY，DZ，SPACE 复制节点ITIME 复制次数，包括本身INC 复制时节点号码增量NODE1，NODE2，NINC 要复制的节点DX，DY，DZ 复制出的节点的位置改变量例：NGEN，4，5，1，5，1，1，2，3 将节点1到5复制4次，每次复制X，Y，Z方向分别移动1，2，3单位长度--------------------------------------------------------FILL，ITIME，INC，NODE1，NODE2，NINC，DX，DY，DZ，SPACE 填充节点（默认为均分填充）例：FILL，1，100 在节点1到100之间填充2，3 (99)--------------------------------------------------------ET，ITYPE，Ename，KOPT1……KOPT6，INOPR 定义元素ITYPE 元素类型编号Ename 所使用元素名称KOPT1-KOPT6 元素特性编码例：ET，1，LINK1 第1类元素为LINK1单元--------------------------------------------------------MP，Lab，MAT，C0，C1，C2，C3，C4 定义材料特性材料特性为固定值，其值为C0材料特性随温度变化，由C1-C4控制Lab 材料特性类别MAT 对应ET所定义的元素类型编号ITYPELab=EX，EY，EZ 杨氏系数Lab=DENS 密度Lab=PRXY，PRYZ，PRZX 泊松比Lab=GXY，GYZ，GZX 剪力模数例：MP，EX，1，207E9 第一类元素的杨氏系数为207E9--------------------------------------------------------R，NSET，R1……R6定义元素类型几何特性NSET 属性组别号码（系统默认值1）R1-R6 所定义元素类型几何特性值例：R，1，1E-4，2.09E-10 ，0.005--------------------------------------------------------E，I，J，K，L，M，N，O，P 定义元素连接方式I-P 定义元素节点顺序的号码例：E，1，2，5，7 四节点元素的节点顺序为1，2，5，7--------------------------------------------------------EPLOT，KNUM 将元素显示在图形窗口中ENUM=0 不显示元素ENUM=1 显示元素--------------------------------------------------------ELIST 将元素资料列在窗口中--------------------------------------------------------EDELE，IEL1，IEL2，INC 删除已建立的元素IEL1，IEL2，INC 欲删除元素的范围例：EDELE，1，10，1--------------------------------------------------------EGEN，ITIME，NINC，IEL1，IEL2，IEINC，MINC，IINC，RINC，CINC，DX，DY，DZ 复制元素ITIME 复制次数，包括本身NINC 复制时节点号码增量IEL1，IEL2，IEINC 欲复制的元素范围DX，DY，DZ 复制出的元素的位置改变量例：EGEN，6，12，1，4，1 将元素1到4复制6次--------------------------------------------------------/PNUM，Label，KEY 在图形中显示号码Label=NODE，ELEM，KP，LINE，AREA，VOLUKEY=0 不显示号码 KEY=1 显示号码例：/PNUM，ELEM，1--------------------------------------------------------/SOLU 进入解题处理器--------------------------------------------------------ANTYPE，Antype，Status 声明分析类型Antype=STATIC or 0 静态分析（系统默认）Antype=BUCKLE or 1 屈曲分析Antype=MODAL or 2 振动模态分析Antype=HARMIC or 3 调和外力动力系统分析Antype=TRANS or 4 瞬时动力系统分析例：ANTYPE，STATIC--------------------------------------------------------F，NODE，Lab，VALUE，VALUE2，NEND，NINC定义节点上的集中力NODE 节点号码Lab 外力形式Lab=FX，FY，FZ，MX，MY，MZ 结构力学Lab=HEAT 热学的热流量Lab=AMP，CHRG 电学的电流，电荷Lab=FLUX 磁学的磁通量VALUE 外力大小NODE，NEND，NINC 施力节点范围例：F，1，FY，，-200，5，1 = F，ALL，FY，-200 节点1-5 的Y方向定义集中力-200（注意FY，，表明VALUE2默认）--------------------------------------------------------FDELE，NODE，Lab，NEND，NINC 删除节点集中力例：FDELE，1，FY，5，1 = FDELE，ALL--------------------------------------------------------D，NODE，Lab，VALUE，VALUE2，NEND，NINC，Lab2，Lab3，Lab4，Lab5，Lab6 定义节点自由度的限制NODE， NEND，NINC 选取自由度约束节点的范围Lab 相对元素每一个节点受自由度约束的形式结构力学Lab=UX，UY，UZ（直线位移）Lab=ROTX，ROTY，ROTZ（旋转位移）例：D，1，UX，，，5，1 节点1到5 X方向约束D，1，UX，，，5，1，UY 节点1到5 X Y方向约束D，1，ALL，，，5，1 节点1到5 全部自由度约束注意：使用命令前要先定义元素--------------------------------------------------------DDELE，NODE，Lab，NEND，NINC 取消节点自由度约束例：DDELE，ALL--------------------------------------------------------DLIST，NODE1，NODE2，NINC 列出节点自由度约束--------------------------------------------------------DL，LINE，AREA，Lab，VALUE，VALUE2定义线自由度限制LINE 线号AREA 线所属面积例：DL，8，3，ALL 定义面积3上面线8的约束注意：同时有DLLIST，DLDELE命令--------------------------------------------------------SFBEAM，ELEM，LKEY，Lab，VALI，VALJ，VAL2I，VAL2J，IOFFST，JOFFST 定义分布力作用于梁元素的方式及大小ELEM 分布力所作用的元素编号LKEY 梁元素的4个面中分布力所在面号码Lab=PRES （表示分布压力）VALI，VALJ 在I，J点分布力的值例：SFBEAM，1，1，PRES，60，30 元素1上1号面作用分布力SFBEAM，1，1，PRES，-30，60--------------------------------------------------------SFE，ELEM，LKEY，Lab，KVAL，VAL1，VAL2，VAL3，VAL4定义分布力作用于元素上的方式及大小ELEM 分布力所作用的元素编号LKEY 分布力作用边，面的号码Lab=PRES （表示分布压力）VAL1-VAL4 分布力在元素边，面上节点的值例：SFE，4，2，PRES，，20，60元素4的第2边，面作用分布力--------------------------------------------------------SF，NLIST，Lab，VALUE，VALUE2 定义节点间分布力NLIST 分布力作用边或面上所有节点Lab=PRES例：SF，ALL，PRES，10注意：SFE适用于非均匀分布力，作用在元素的边或面上SF适用于均匀分布力，作用在节点之间--------------------------------------------------------SFLIST，NODE，Lab 显示分布力--------------------------------------------------------NSEL，Type，Item，COMP，VMIN，VMAX，VINC，KABS节点选择命令Type 选择方式Type=S 选择一组节点为ACTIVE点Type=R 在现有ACTIVE点中选出一部分作为ACTIVE点Type=A 在选择一部分节点，加入ACTIVE点中Type=U 在现有ACTIVE点中，排除某些节点Type=ALL 选择全部节点作为ACTIVE点Item =NODE 用节点号码选取┇ Item =LOC 用节点坐标选取COMP=无┇ COMP=X 以节点X坐标为准VMIN，VMAX，VINC 节点选取范围例：NSEL，S，NODE，，1，13，1选1-13中奇数点为ACTIVE点NSEL，A，NODE，，14，20，1 选14-20加入ACTIVE点--------------------------------------------------------OUTPR，Item，FREQ，Cname 分析结果是否显示于输出窗口Item=ALL 所有结果 Item=NSOL 节点自由度结果FREQ 负载的次数 FREQ=ALL 最后负载例：OUTPR，ALL，ALL注意：仅用于小程序--------------------------------------------------------SOLVE 开始解题--------------------------------------------------------/POST1 进入后处理器--------------------------------------------------------PLDISP，KUND 显示结构变形结构KUND=0 显示变形后结构形状KUND=1 显示变形前后结构形状KUND=0 显示变形前后结构形状，但仅显示外观--------------------------------------------------------PLESOL，Item，Comp 显示元素的解答Item（何种解答） Comp（Item分量）S X，Y，Z，XY，YZ，ZX 应力S 1，2，3 主应力S EQV，INT 等效应力F X，Y，Z 结构力M X，Y，Z 结构力矩例：PLESOL，S，X，Y，Z--------------------------------------------------------PLNSOL，Item，Comp 显示节点的解答Item（何种解答） Comp（Item分量）S X，Y，Z，XY，YZ，ZX 应力S 1，2，3 主应力S EQV，INT 等效应力F X，Y，Z 结构力M X，Y，Z 结构力矩U X，Y，Z，SUM 位移ROT X，Y，Z，SUM 旋转位移例：PLNSOL，S，Y PLNSOL，U，X--------------------------------------------------------PRESOL，Item，Comp 打印元素解答Item（何种解答） Comp（Item分量）S X，Y，Z，XY，YZ，ZX 应力F X，Y，Z 结构力M X，Y，Z 结构力矩例：PRESOL，S，X--------------------------------------------------------PRNSOL，Item，Comp 打印节点解答Item（何种解答） Comp（Item分量）U X，Y，Z 位移U X，Y，Z方向及总向量方向位移S COMP 应力S PRIN 主应立，等效应力例：PRNSOL，U PRNSOL，S注意：查看结果通常使用PLDISP，1 PLNSOL，U，Y…… PRNSOL，S 其中PLNSOL中Comp不能省略，PRNSOL中可以省略--------------------------------------------------------TYPE，Itype 声明建立元素时，元素形式号码（对应ET的Itype）例：ET，1，LINK1 ET，2，PLANE42TYPE，1建立LINK1元素 TYPE，2建立PLANE42元素--------------------------------------------------------REAL，NSET 声明建立元素时，元素几何参数属性编号（对应R，NSET）--------------------------------------------------------MAT，MAT 声明建立元素时，元素材料特性属性编号（对应MP，MAT）--------------------------------------------------------LSWRITE，LSNUM 多重负载资料保存至文件Jobname.S0i例：LSWRITE 自动编号，不用输入--------------------------------------------------------LSSOLVE，SLMIN，LSMX，LSINC读取所定义的多重负载并解答SLMIN，LSMX，LSINC 读取负载的范围例：LSSOLVE，1，4 获得1-4负载的解答--------------------------------------------------------DDELE，NODE，Lab，NEND，NINC 删除定义的约束条件NODE， NEND，NINC 删除约束的节点范围Lab 删除约束的节点方向例：DDELE，1，UY，5，2 删除1，3，5节点的Y向约束--------------------------------------------------------FDELE，NODE，Lab，NEND，NINC 删除定义的集中力NODE， NEND，NINC 删除集中力的范围Lab 删除外力的方向例：FDELE，1，FY，5，2 删除1，3，5节点的Y向集中力--------------------------------------------------------SFDELE，Nlist，Lab 将已定义的面载荷删除Nlist 面负载所含的节点（由NSEL选择，设Nlist=ALL）Lab=PRES （结构力学）--------------------------------------------------------SFEDELE，ELEM，LKEY，Lab 将负载从元素上删除ELEM 元素编号LKEY 负载作用边，面的号码Lab=PRES （结构力学）--------------------------------------------------------SET，Lstep，SBSTEP，FACT，KIMG，TIME，NGLE，NSET检查负载结果Lstep=负载编号例：SET，2 检查第二负载的结果 PLDISP，1……注意：此命令要在/POST1中使用--------------------------------------------------------FILE，Fname，Ext，Dir 读取分析后的结果文件例：FILE，TEST，RST--------------------------------------------------------ANTYPE，Antype，Status 声明分析类型例：ANTYPE=MODAL or 2 模态分析--------------------------------------------------------MODOPT，Method，NMODE，FREQB，FREQE，PRMODE，NUMKEY 选择模态分析方法Method=REDUC 降阶法Method=SUBSP 次空间法NMODE 欲求模态个数（降阶法小于主自由度一半）FREQB，FREQE 欲探讨振动频率范围（默认全部）PRMODE 分析后，模态结果保存到结果文件的个数例：MODOPT，SUBSP，5--------------------------------------------------------M，NODE，Lab1，NEND，NINC，Lab2。

LS―DYNA主要版本冲压仿真评测

LS―DYNA主要版本冲压仿真评测ls-dyna是通用显式动力分析程序，能够模拟真实世界的各种复杂问题，特别适合求解各种二维、三维非线性结构的高速碰撞、爆炸和金属成型等非线性动力冲击问题，同时可以求解传热、流体及流固耦合问题，在工程应用领域被广泛应用。

ls-dyna程序以lagrange算法为主，兼有ale和euler算法；以显式求解为主，兼有隐式求解功能；以结构分析为主，兼有热分析、流体―结构耦合功能；以非线性动力分析为主，兼有静力分析功能（如动力分析前的预应力计算和薄板冲压成型后的回弹计算），是军用和民用相结合的通用结构分析非线性有限元程序。

与实验的无数次对比证实了ls-dyna计算的可靠性，截止到目前为止的最新版本是2013年初推出的l s971_r7.0版本（目前为测试版），比较成熟稳定的版本是l s971的r3.2.1/r4.1.1/r5.1.1/r6.1.1等版本。

目前在冲压仿真方面以ls-dyna为主要求解器的软件有：dynaform、jstamp、hyperform 等，其中dynaform为其冲压领域的主要合作开发商。

通过实际应用证明，不同版本的ls-dyna 求解器在分析时间、分析结果质量上有重大的差异，甚至某些版本都不能通过正常的计算得出正确的结果（如l s971_r4.x.x系列的某些版本），为了保证冲压仿真分析结果的可靠性，就有必要对求解器的速度、稳定性等有全面的了解。

本文笔者使用软件自带的例子以及numisheet2005年大会中的算例，对冲压仿真领域相关的重力、拉延和回弹等常用的功能进行综合评比测试，以对比各个版本的求解器（r3.2.1/r5.1.1/r6.1.1）之间的差异，以便于实际的工作和应用。

一、ls-dyna软件支持平台简介二、ls-dyna软件的分类ls-dyna软件从计算方式上分为smp、mpp、hybrid等三大类；smp版本主要为在一台运算服务器上，支持多个cpu进行并行运算；m p p版本指使用多台相同系统的运算服务器进行联合仿真分析；linux和unix操作系统还有hybrid方式，指的是不同操作系统平台的联合分析。

ls-dyna命令帮助手册（中文）

Fini(退出四大模块，回到BEGIN层)/cle (清空内存，开始新的计算)1．定义参数、数组，并赋值.2．/prep7(进入前处理)定义几何图形：关键点、线、面、体定义几个所关心的节点，以备后处理时调用节点号。

设材料线弹性、非线性特性设置单元类型及相应KEYOPT设置实常数设置网格划分，划分网格根据需要耦合某些节点自由度定义单元表3．/solu加边界条件设置求解选项定义载荷步求解载荷步4./post1（通用后处理）5./post26 （时间历程后处理）6.PLOTCONTROL菜单命令7.参数化设计语言8.理论手册Fini(退出四大模块，回到BEGIN层)/cle (清空内存，开始新的计算)1 定义参数、数组，并赋值.u dim, par, type, imax, jmax, kmax, var1, vae2, var3 定义数组par: 数组名type：array 数组，如同fortran,下标最小号为1，可以多达三维（缺省）char 字符串组（每个元素最多8个字符）tableimax,jmax, kmax 各维的最大下标号var1,var2,var3 各维变量名，缺省为row,column,plane(当type为table时) 2 /prep7(进入前处理)2.1 定义几何图形：关键点、线、面、体u csys,kcnkcn , 0 迪卡尔坐标系1 柱坐标2 球4 工作平面5 柱坐标系（以Y轴为轴心）n 已定义的局部坐标系u numstr, label, value设置以下项目编号的开始nodeelemkplineareavolu注意：vclear, aclear, lclear, kclear 将自动设置节点、单元开始号为最高号，这时如需要自定义起始号，重发numstru K, npt, x,y,z, 定义关键点Npt：关键点号，如果赋0，则分配给最小号u Kgen,itime,Np1,Np2,Ninc,Dx,Dy,Dz,kinc,noelem,imoveItime：拷贝份数Np1,Np2,Ninc：所选关键点Dx,Dy,Dz：偏移坐标Kinc：每份之间节点号增量noelem: “0” 如果附有节点及单元，则一起拷贝。

LS-DYNA使用指南

ANSYS/LS-DYNA 使用指南前言本资料来源于互联网，据说来源于安世亚太，基本上是ANSYS 的帮助文件的翻译。

本人遂排版整理，对原文中的公式进行了编辑排版，并对部分错误进行了修改，当然，错误在所难免，仅供自己查阅学习。

版权归原作者所有！如有人非法用于商业用途，与本人无关。

人就一个字2009-5-7目录第一章引言 (1)1.1 显式动态分析求解步骤概述 (1)1.2 显式动态分析采用的命令 (1)1.3 本手册使用指南 (4)1.4 何处能找到显式动态例题 (5)1.5 其它信息 (5)第二章单元 (6)2.1 实体单元和壳单元 (7)2.1.1 SOLID164 (7)2.1.2 SHELL163 (8)2.1.3 通用壳单元算法 (8)2.1.4 薄膜单元算法 (9)2.1.5 三角型薄壳单元算法 (9)2.1.6 PLANE162 (12)2.2 梁单元和杆单元 (13)2.2.1 BEAM161 (13)2.2.2 LINK160 (14)2.2.3 LINK167 (14)2.3 离散单元 (15)2.3.1 COMBI165弹簧-阻尼单元 (15)2.3.2 MASS166 (16)2.4 一般单元特性 (16)第三章建模 (17)3.1 定义单元类型和实常数 (17)3.2 定义材料特性 (18)3.3 定义几何模型 (18)3.4 网格划分 (18)3.5 定义接触面 (20)3.6 建模的一般准则 (20)3.7 PART的定义 (21)3.7.1 Part集合 (24)3.8 自适应网格划分 (24)第四章加载 (28)4.1 一般载荷选项 (28)4.1.1 组元 (29)4.1.2 数组参数 (30)4.1.3 施加载荷 (31)4.1.4 数据曲线 (34)4.1.5 在局部坐标系中定义载荷 (36)4.1.6 指定Birth和Death时间 (37)4.2 约束和初始条件 (37)4.2.1 约束 (37)4.2.2 焊接 (39)4.2.3 初始速度 (39)4.3 耦合和约束方程 (41)4.4 非反射边界 (42)4.5 温度载荷 (42)4.6 动力松弛 (44)第五章求解特性 (46)5.1 求解过程 (46)5.2 LS-DYNA终止控制 (46)5.3 共享存储器并行处理 (47)5.4 求解控制和监控 (48)5.5 显示小尺寸单元 (50)5.6 编辑LS-DYNA的输入文件 (50)5.6.1 方法A (51)5.6.2 方法B (51)5.6.3 使用预先存在的FILE.K (52)第六章接触表面 (54)6.1 接触的定义 (54)6.1.1 列表，显示和删除接触实体 (57)6.2 接触选项 (58)6.2.1 定义接触类型 (59)6.2.2 定义接触选项 (60)6.3 接触搜索方法 (63)6.3.1 网格连接跟踪 (63)6.3.2 批处理方法 (63)6.3.3 限制接触搜索域 (63)6.4 壳单元的特殊处理 (64)6.5 接触深度控制 (64)6.6 接触刚度 (65)6.6.1 罚因子的选择 (65)6.6.2 对称刚度 (66)6.7 2D接触选项 (66)第七章材料模型 (67)7.1 定义显示动态材料模型 (69)7.2 显式动态材料模型的描述 (70)7.2.1线弹性模型 (70)7.2.2非线性弹性模型 (72)7.2.3 非线性无弹性模型 (74)7.2.4 压力相关的塑性模型 (86)7.2.5 泡沫模型 (90)7.2.6 状态方程 (94)7.2.7 离散单元模型 (104)7.2.8 刚性体模型 (107)第八章刚性体 (109)8.1 定义惯性特性 (109)8.2 加载 (111)8.3 变形体和刚性体部件间的转换 (111)8.4 节点刚体 (112)第九章沙漏 (113)第十章质量缩放 (115)第十一章子循环 (117)第十二章后处理 (119)12.1 输出控制 (119)12.1.1 结果文件(Jobname.RST)和时间历程文件(Jobname.HIS)的比较 (119)12.1.2 生成POST26的ComponentS (120)12.1.3 为POST26记录输出文件 (120)12.2 使用ANSYS/LS-DYNA的POST1 (121)12.2.1 动画结果 (121)12.2.2 单元数据输出 (122)12.2.3 自适应网格划分的处理结果 (123)12.3 在ANSYS/LS-DYNA中使用POST26 (124)第十三章重启动 (125)13.1 重启动Dump文件 (125)13.2 EDSTART 命令 (125)13.2.1 新分析 (126)13.2.2 简单重启动 (126)13.2.3 小型重启动 (126)13.2.4 完全重启动 (128)13.3 输出文件的影响 (130)第十四章显式－隐式顺序求解 (132)14.1 显式－隐式顺序求解 (132)第十五章隐式-显式顺序求解 (137)15.1 预载荷结构的隐式-显式顺序求解 (137)第十六章跌落测试模块 (142)16.1 简介 (142)16.2 选择DTM模块启动ANSYS (142)16.3 典型的跌落分析步骤 (143)16.3.1 基本的跌落测试分析步骤 (143)16.3.2 屏幕坐标的定义 (146)第一章引言ANSYS/LS-DYNA将显式有限元程序LS-DYNA和ANSYS程序强大的前后处理结合起来。

第十七章LS-DYNA的隐式求解.

第⼗七章LS-DYNA的隐式求解.第⼗七章 LS-DYNA 的隐式求解LS-DYNA作为著名的显式求解程序只能求解瞬态动⼒问题，对于时间历程较长的静、动⼒问题， LS-DYNA的显式中⼼差分法有它的局限性，⽽⼀些与瞬态动⼒分析紧密相关的问题要求LS-DYNA也能够求解，如：冲压成型后的回弹计算应⼒初始化冲击后常时间低频动⼒响应静⼒分析特征值分析实际上从950版本开始，LS-DYNA已增加了隐式求解功能。

刚开始的应⽤主要在冲压成型后的回弹计算，经过960版到970版的发展，LS-DYNA的隐式求解功能已⼤增强，已经能满⾜以上的求解需要。

17.1显式与隐式的区别：17.1.1 LS-DYNA显式求解：采⽤中⼼差分⽅法进⾏显式时间积分Man fn⽅程⾮耦合，可以直接求解(显式)但需要常⼩的时间步保持稳定状态不需要求解刚度矩阵 ext fnint适合冲击、穿透等⾼频⾮线性动⼒响应问题17.1.2 LS-DYNA隐式求解：采⽤Newmark隐式时间积分M an K un 1 f对于线性问题，⽆条件稳定可以采⽤⼤的时间步 extn 1 fintn Man对于⾮线性问题，需要⼀系列线性逼近（Newton-Raphson）叠代求解需要求解刚度矩阵适合静⼒问题、低频动⼒问题及特征值分析。

17.2 LS-DYNA中隐式分析的激活及相关关键字在LS-DYNA中，缺省的求解是显式求解，可以通过下⾯的关键字来激活隐式求解：*CONTROL_IMPLICIT_GENERAL*CONTROL_IMPLICIT_GENERAL$ imflag dt0 iefs nstepsb igso1 0.01 0 0 0其中参数imflag=1激活全隐式求解imflag=0为缺省的显式求解。

imflag=2为显式求解后⽆缝进⾏隐式回弹求解。

该关键字对于所有隐式求解分析来说都是必需的。

与隐式求解相关的其它关键字：*CONTROL_IMPLICIT_LINEAR(v960版本改为*CONTROL_IMPLICIT_SOLVER)*CONTROL_IMPLICIT_NONLINEAR(在v960版后改为*CONTROL_IMPLICIT_SOLUTION)*CONTROL_IMPLICIT_AUTO*CONTROL_IMPLICIT_STABILIZATION*CONTROL_IMPLICIT_DYNAMICS*CONTROL_IMPLICIT_EIGENVALUE*CONTROL_IMPLICIT_BUCKLE（v970）*CONTROL_IMPLICIT_MODES（v970）使⽤*CONTROL_IMPLICIT_GENERAL激活隐式求解后，还需要与上⾯的相关关键字进⾏联合，从⽽对不同的问题设置不同的关键字进⾏求解。

LSdyna精编教材

LSdyna精编教材LS-DYNA 精编教材讲解LS-DYNA的基础知识.LS-DYNA 软件 1.1 LS-DYNA 简介 LS-DYNA 是世界上最著名的通用显式动力显式动力分析程序，能够模拟真实世界的显式动力各种复杂问题，特别适合求解各种二维、三维非线性结构的高速碰撞、爆炸和金属成型等非线性动力冲击问题，同时可以求解传热、流体及流固耦合问题。

在工程应用领域被广泛认可为最佳的分析软件包。

与实验的无数次对比证实了其计算的可靠性。

由 J.O.Hallquist 主持开发完成的 DYNA 程序系列被公认为是显式有限元程序的鼻祖和理论先导，是目前所有显式求解程序（包括显式板成型程序）的基础代码。

1988 年 J.O.Hallquist 创建 LSTC 公司，推出 LS-DYNA 程序系列，并于 1997 年将 LS-DYNA2D、LS-DYNA3D、LS-TOPAZ2D、LS-TOPAZ3D 等程序合成一个软件包，称为 LS-DYNA。

LS-DYNA 的最新版本是 2021 年 8 月推出的 970 版。

1.1.1 LS-DYNA 功能特点 LS-DYNA 程序是功能齐全的几何非线性（大位移、大转动和大应变）、材料非线性（140 多种材料动态模型）和接触非线性（50 多种）程序。

它以Lagrange 算法为主，兼有 ALE 和 Euler 算法；以显式求解为主，兼有隐式求解功能；以结构分析为主，兼有热分析、流体-结构耦合功能；以非线性动力分析为主，兼有静力分析功能（如动力分析前的预应力计算和薄板冲压成型后的回弹计算）；军用和民用相结合的通用结构分析非线性有限元程序。

LS-DYNA 功能特点如下： 1.分析能力：非线性动力学分析多刚体动力学分析准静态分析（钣金成型等） 1 热分析结构-热耦合分析流体分析：欧拉方式任意拉格郎日-欧拉（ALE）流体-结构相互作用不可压缩流体 CFD 分析有限元-多刚体动力学耦合分析（MADYMO，CAL3D）水下冲击失效分析裂纹扩展分析实时声场分析设计优化隐式回弹多物理场耦合分析自适应网格重划并行处理（SMP 和 MPP） 2.材料模式库(140 多种) 金属塑料玻璃泡沫编制品橡胶（人造橡胶）蜂窝材料复合材料混凝土和土壤炸药推进剂 2 粘性流体用户自定义材料3.单元库体单元薄/厚壳单元梁单元焊接单元离散单元束和索单元安全带单元节点质量单元 SPH 单元 4.接触方式(50 多种) 柔体对柔体接触柔体对刚体接触刚体对刚体接触边-边接触侵蚀接触充气模型约束面刚墙面拉延筋 5.汽车行业的专门功能安全带滑环预紧器牵引器传感器加速计 3 气囊混合 III 型假人模型 6.初始条件、载荷和约束功能初始速度、初应力、初应变、初始动量（模拟脉冲载荷）；高能炸药起爆；节点载荷、压力载荷、体力载荷、热载荷、重力载荷；循环约束、对称约束（带失效）、无反射边界；给定节点运动（速度、加速度或位移）、节点约束；铆接、焊接（点焊、对焊、角焊）；二个刚性体之间的连接－球形连接、旋转连接、柱形连接、平面连接、万向连接、平移连接；位移/转动之间的线性约束、壳单元边与固体单元之间的固连；带失效的节点固连。

LS-DYNA资料

1.1LS-DYNA发展概况LS-DYNA是以显式为主、隐式为辅的通用非线性动力分析有限元程序，特别适合求解各种二维、三维非线性结构的高速碰撞、爆炸和金属成形等非线性动力冲击问题，同时可以求解传热、流体及流固耦合问题。

DYNA程序系列最初是1976年在美国Lawrence Livermore National Lab. 由J.O.Hallquist博士主持开发完成的，主要目的是为武器设计提供分析工具，后经1979、1981、1982、1986、1987、1988年版的功能扩充和改进，成为国际著名的非线性动力分析软件，在武器结构设计、内弹道和终点弹道、军用材料研制等方面得到了广泛的应用。

1988年J.O.Hallquist创建LSTC公司，推出LS-DYNA程序系列，主要包括显式LS-DYNA2D、LS-DYNA3D、隐式LS-NIKE2D、LS-NIKE3D、热分析LS-TOPAZ2D、LS-TOPAZ3D、前后处理LS-MAZE、LS-ORION、LS-INGRID、LS-TAURUS等商用程序，进一步规范和完善DYNA 的研究成果，陆续推出930版（1993年）、936版（1994年）、940版（1997年），950版（1998年）增加了汽车安全性分析（汽车碰撞、气囊、安全带、假人）、薄板冲压成形过程模拟以及流体与固体耦合（ALE和Euler算法）等新功能，使得LS-DYNA程序系统在国防和民用领域的应用范围进一步扩大，并建立了完备的质量保证体系。

1997年LSTC公司将LS-DYNA2D、LS-DYNA3D、LS-TOPAZ2D、LS-TOPAZ3D等程序合成一个软件包，称为LS-DYNA， LS-DYNA的最新版本是2008年5月推出的971版，它在970版基础上增加了不可压缩流体求解程序模块，并增加了一些新的材料模型和新的接触计算功能，详见以下介绍。

LS-DYNA程序971版是功能齐全的几何非线性（大位移、大转动和大应变）、材料非线性（140多种材料动态模型）和接触非线性（40多种接触类型）程序。

LS DYNA软件简介及相应分析实例分析

LS-DYNA软件简介及相应分析实例
报告人：陈林
主要内容
LS-DYNA发展历史 LS-DYNA与ANSYS等软件的关系利用LS-DYNA进行结构分析的主要步骤 LS-DYNA的学习难点 LS-DYNA学习所需的相关资料目前本人所作的一些相关分析工作
LS-DYNA发展历史
ANSYS前处理
Ultraedit修改k文件
LS-DYNA求解
LS-PREPOST后处理主要流程图
1、通过桌面开始菜单进入ANSYS操作平台界面（Mechanical APDL Product Launcher）
2、进入ANSYS建模主界面
3、输出关键字文件（即k文件）
4、利用Ultraedit文本编辑器进行k文件的修改
LS-DYNA与ANSYS等软件的关系
这里主要是给那些还不怎么了解这个软件的同学理清几个概念，即通常所说的ANSYS与LS-DYNA究竟具有什么样的关系？利用LS-DYNA进行分析的过程中还会用到什么软件？
其实通常情况下，我们所说的LS-DYNA程序是指其核心计算程序，例如最新的版本LS971 R5.0.exe。而当 ANSYS与LS-DYNA一起被提及时，这里的ANSYS一般是指其前处理器。换句话说，我们进行动力分析时，首先通过ANSYS的前处理器进行有限元网格划分等工作，然后再调用LS971 R5.0.exe进行计算。
目前关于LS-DYNA的参考书不多，特别是中文的参考书较少，这就需要学习者经常阅读英文参考资料。
LS-DYNA的学习难点
接触问题。与ANSYS等静力分析软件采接触单元定义接触不同，LSDYNA处理接触问题的方式是定义接触面，它不需要额外的接触单元，只需要将分析中可能接触的两个面定义一个接触属性即可。具体体现在k文件中的关键字*contact。该关键字包含了大量的参数，在简单问题中，采用默认的参数即可，但是很多情况下，需要人工控制。这就需要对 LS-DYNA关于接触的理论部分有一个较为清晰的理解。这应该是处理碰撞问题时可能会遇到的最大难点！

DynaForm各版本LS-DYNA求解器对比测试

DynaForm发展到5.8.1，已经有多个求解器，如快速展开求解器MSTEP、切边、冷却求解器UtilityBatch、还有用于快速计算的DEMS_SC求解器，但是其最主要的重力、拉延、回弹等求解器是LS-DYNA；LS-DYNA是一个非常出名的非线性求解器，应用的领域很广，冲压方面的功能仅是是一部分功能，目前其主要的版本为LS971，LS971已经发展到6.X.X系列，现比较稳定的版本有R3.2.1 \R4.2.1 \R5.1.以及最新版R6.0.0，DynaForm不同的版本佩带不同版本的LS-DYNA求解器；从DynaForm5.8开始默认支持R5/R6系列的求解器，而5.8之前的DynaForm版本支持R3系列和R4系列的LS-DYNA；为了更好的应用DynaForm软件，就必须对求解器的速度、稳定性等有全面的了解，笔者使用同一个例子，不同版本的求解器，对冲压相关的重力、拉延、回弹等三大常用的功能进行综合评比测试，以对比各个版本的求解器之间的差异，以便于工作和应用；测试环境LS-Dyna V971 R3.2.1和R4.2.1使用DynaForm5.7.3作为前处理环境；LS-Dyna V971 R5.1.1和R6.0.0使用DynaForm5.8作为前处理环境；操作系统为WIN7_32_SP1内存：4GCUP：17 Q620M对比统计结果：计算时间、计算是否正常终止、计算结果是否正常、板料厚度、最大位移、FLD曲线等；1.LS-DYNA 32位和64位测试32位的LS-DYNA可以在win32和win64操作系统上运行，而64位的LS-DYNA只能在X64的操作系统上，对于安装64位Windows系统的用户,可以同时使用32位和64位的LS-DYNA，所以下面就这两者之间做一个对比测试笔者使用NUMISHEET 2005 BM2为例，生成一个DYN 文件（包含拉延、切边、回弹），然后用32位的LS-DYNA和64位的LS-DYNA不同版本（LS971R3.2.1/LS971R5.1.1）进行运算，记录运算的时间以及厚度和回弹量。

ls-dyna关键字命令介绍

LIVERMORE SOFTWARE TECHNOLOGY CORPORATION
LS-DYNA 输入文件基础
关键字命令
LIVERMORE SOFTWARE TECHNOLOGY CORPORATION
关键字语法

大部分的命令是由下划线分开的字符串, 如 *boundary_prescribed_motion_rigid. 字符可以是大写或小写在用户手册中，命令是按字母顺序排列的在输入文件中，命令的顺序是不重要的 (除了 *keyword 和 *define_table) 关键字命令必须左对齐，以*号开始第一列的“$”表示该行是注释行输入的参数可以是固定格式或者用逗号分开空格或者0 参数 Æ 使用该参数的默认值
LIVERMORE SOFTWARE TECHNOLOGY CORPORATION
*MAT_option

Option 指定材料类型, 如 *mat_elastic, *mat_rigid, *mat_linear_piecewise_plasticity 输入参数 (材料常数) 取决于材料模型大约有 100 种结构材料可以选择 8种弹簧/阻尼材料 1 种安全带材料 6 种热材料 (用于热传递分析) 声学材料
LIVERMORE SOFTWARE TECHNOLOGY CORPORATION
定义单元
*ELEMENT_... eid, pid, n1, n2, etc *PART pid, sid, mid, eosid, hgid, … *NODE nid, x, y, z *HOURGLASS hgid, … *SECTION_... sid, elform, … *EOS_GRUNEISEN eosid, … *MAT_ … mid, rho, …

LS-DYNA常见问题简化中文版

一、如何得知LS-DYNA为一套完整的软件？<TOP>解答：LS-DYNA一直以来都是独立软件，原因很简单：(1)美国原厂Livermore Software Technology Corp.(简称LSTC)授权「势流科技」为台湾地区LS-DYNA软件的独家代理商。

(2)LS-DYNA拥有LSTC官方的前、后处理程序：FEMB及LS-Pre/Post，无需另外增购。

LS-DYNA强调「One Model、One Code、Multi-results」，单一核心程序即可进行Implicit、Explicit及混合Implicit/Explicit的求解方式，可以为您解决静力、准静力、掉落、冲击、模态分析、振动、预应力/预应变、流固耦合、流场、热应力....等等的力学分析，使用者不必再额外添购模块程序，徒增预算的浪费与软件学习的负担。

二、LS-DYNA与市面上其它的前处理软件兼容性如何？<TOP>解答：由于LS-DYNA是全球使用率最高的结构分析软件，因此其它的前处理软件与LS-DYNA是完全兼容的。

在此要强调的是：LS-DYNA的官方前处理程序为FEMB，因为FEMB是专门为LS-DYNA量身订作的前处理程序，有许多设定条件及架构逻辑是其它前处理软件所难望其项背的，为了避免在学习LS-DYNA的过程及操作上产生困扰，强烈建议使用者采用原厂出品的FEMB来做为LS-DYNA的前处理工具，使用者必定更能体会LS-DYNA直觉式的设定与强大的分析能力。

.三、LS-DYNA似乎很重视「Contact Algorithm」，这是为什幺？<TOP>解答：是的，LS-DYNA很早以前就已经发展「接触算法」，这是因为基础力学所分析的对像均只考虑「力的受体」，故输入条件皆为外力量值。

然而在真实情况下，物体受力通常是因为与其它的物体发生「接触」（Contact）才受力，此时外力量值是无法预期的，应该输入的条件往往都是几何上的接触条件。