ansys_ls-dyna材料模型与单元类型
ansys中LS-DYNA2D金属切削模拟步骤
ansys中LS-DYNA2D金属切削模拟步骤在ANSYS Launcher界面中，选择ANSYS Mechanical/LS-DYNA1、菜单过滤Main Menu→Preprocessor→LD-DYNA Explicit→OK2、设置文件名及分析标题Utility Menu→File→change Jobname→2D cutting→New log and error file ：YES→OK Utility Menu→File→change Title→cutting analysis →OK3、选择单元类型Main menu→preprocessor→Element Type→Add/Edit/Delete→Add→2D solid 162→OK→options→选择const.stress ；Lagrangian→OK4、定义材料模型（1）定义刀具材料模型Main menu→preprocessor→Material Props→Material Models→rigid material→输入：DENS：5.2e3 ；EX：4.1e11 ；NUXY：0.3 ；选择“Y and Zdisps” ；“All rotati ons”→OK（2）定义工件Johnson-cook材料模型Main menu→preprocessor→Material Props→Material Models→Gruneisen→Johns on-cook→输入：DENS：7.8e3 ；EX：2.06e11 ；NUXY：0.3A：507；B：320；C：0.28；n；0.064；m=1.06D1：0.15；D2：0.72；D3：1.66；D4：0.005；D5：--0.845、创建几何模型（1）创建工件模型Main menu→preprocessor→Create→Areas→Rectangle→By Dimensions→输入：X1，X2：0，5；Y1，Y2：0，3→OK（2）创建刀片模型Main menu→preprocessor→Create→Keypionts→In Active CS→依次输入：keypoint number：5，X、Y、Z ：5.1，2.9，0；keypoint number：6，X、Y、Z ：6，3.228，0；keypoint number：7，X、Y、Z ：6，4，0；keypoint number：8，X、Y、Z ：5.294，4，0→OK6、网格划分（一）（1）对刀片进行网格划分Utility Menu→Select→Entities→Lines ：By Num/Pick→Apply→选取刀片边线→O KMain menu→preprocessor→Meshing→Size contrls→Manualsize→Lines→All lin es→NDIV：10→OK （2）对刀尖半圆进行网格划分Utility Menu→Select→Entities→Lines ：By Nu m/Pick→Apply→选取刀尖半圆→O KMain menu→preprocessor→Meshing→Size contrls→Manualsize→Lines→All lin es→NDIV：3→OK （3）确定刀片的单元属性Main menu→preprocessor→Meshing→Mesh Attributes→Picked Aeras→选取刀片→Apply→确定材料号和单元类型号为1→OK（4）刀片网格划分Main menu→preprocessor→Meshing→MeshT ool→Mesh：Aeras；shape：Tri；free →Mesh→选取刀片→OK（二）（5）对工件进行网格划分切分工件Utility menu →Workplane→Wp settings→Grid andT riad→Minimum ,maximum:-5, 5 ;Spacing:1.0→OK平移和旋转工作平面并用其切分工件Utility menu →Workplane→Offset wp by incremens→X，Y，Z offsets：0，2.5，0；XY，YZ，ZX angle：0，90，0→OK Mainmenu→preprocessor→Modeling→operate→Booleans→Divide→Areas by wkp lane→选取工件→OK取消工作平面显示Utility menu→workplane→Display workingplane→等分接触区域相关Y向线段Utility Menu→Select→Entities→Lines ：By Num/Pick→Apply→选取工件接触区Y向线段（两条）→OK Main menu→preprocessor→Meshing→Size contrls→Manualsize→Lines→All lin es→NDIV：10→OK 等分接触区域相关X向线段Utility Menu→Select→Entities→Lines ：By Num/Pick→Apply→选取工件接触区X向线段（两条）→OK Main menu→preprocessor→Meshing→Size contrls→Manualsize→Lines→All lin es→NDIV：40→OK 等分接触区域不相关Y向线段Utility Menu→Select→Entities→Lines ：By Num/Pick→Apply→选取工件接触区Y向线段（两条）→OK Main menu→preprocessor→Meshing→Size contrls→Manualsize→Lines→All lin es→NDIV：25→OK 等分接触区域不相关X向线段Utility Menu→Select→Entities→Lines ：By Num/Pick→Apply→选取工件接触区X向线段（底边）→OK Main menu→preprocessor→Meshing→Size contrls→Manualsize→Lines→All lin es→NDIV：30→OK确定工件的单元属性Main menu→preprocessor→Meshing→Mesh Attribu tes→Picked Aeras→选取工件→Apply→确定材料号为2和单元类型号为1→OK工件网格划分Main menu→preprocessor→Meshing→MeshT ool→Mesh：Aeras；shape：Quad；mapp ed→Mesh→选取工件→OK7、建立partMain menu→preprocessor→LS-DYNA options→part options→create all part→O K(part1:刀具；part2：工件)Plot→parts（不同颜色显示单元）8、定义接触信息Main menu→preprocessor→LS-DYNA options→contact→Define contact→surface to surf；Eroding；静、动摩擦系数为0.15、0.10→OK→弹出contact options对话框，确定接触件（工件），目标件（刀片）→OK9、施加边界条件Utility menu→select→Entities→Nodes :By Location :X Coordinates→Min,Max: -0.01,0.01;Fro m Full→Apply（选中左侧边所有节点）Main menu→preprocessor→LS-DYNA options→Constraints→Apply→on nodes→pi ck All→All DOF→OKUtility menu→select→Entities→Nodes :By Location :Y Coordinates→Min,Max: -0.01,0.01;From Full→Apply（选中底边所有节点）Main menu→preprocessor→LS-DYNA options→Constraints→Apply→on nodes→pi ck All→All DOF→OK恢复整个模型的选择Utility menu→select→Everything10、对刀片施加初速度Main menu→preprocessor→LS-DYNA o ptions→Initial Velocity→on parts→w/No dal Rotate→选择part1，VX：-100→OK恢复整个模型的选择Utility menu→select→Everything11、设置能量控制选项Main menu→Solution→Analysis options→Energy options→打开所有能量控制选项→OK12、设置人工体积粘性选项Main menu→Solution→Analysis options→Bulks viscosity→Quadratic Viscosit y Coefficient:1.0→OK13、设置时间步长因子Main menu→Solution→Time controls→Time step ctrls→Time step scale facto r:0.6→OK14、设置求解时间Main menu→Solution→Time controls→Solution time→1e-3→OK15、设置结果文件输出步数Main menu→Solution→Output Controls→File output Freq→Number of steps→[EDRST]:50;[EDHTIME]:50→OK16、设置结果文件的输出类型Main menu→Solution→Output Controls→Output File Types→Add：ANSYS and LS -DYNA→OK17、输出K文件Main menu→Solution→Write jobname.K18、求解Main menu→Solution→Solve19、后处理（暂时不管。
LS-DYNA(ANSYS)教程
F(t)
•
•
声波的传播
失效分析
应力波的传播 2500 体单元. CPU 时间2 sec (SGI Octane R12000)
March 7, 2002 Inventory #001630 1-11
概述
C. 显式和隐式方法的比较
静态 “准” 静态 动态
PUNCH
BLANK DIE
概述
A. 什么是 ANSYS/LS-DYNA?
• • 显式动力学有限元程序的基本目标 用来求解高度非线性瞬态动力学问题
– 对大范围的接触类型有效
– 丰富先进的材料模型 – 强大的大变形分析能力
•
ANSYS 和LS-DYNA 程序的无缝结合
– 将LS-DYNA 求解器完全集成到 ANSYS程序中 – 所有的前后处理使用ANSYS标准格式 – GUI 完全类似于通用ANSYS程序
结构问题
金属成形
碰撞问题
SF=0 隐式方法
SF 0
S F = ma
显式方法
March 7, 2002 Inventory #001630 1-12
概述
… 显式－隐式方法的比较
隐式时间积分: • 时间t+Dt 时计算平均加速度－位移:
1 a u K F t D t t D t
March 7, 2002 Inventory #001630 1-16
概述
D. 临界时间步
• 杆件的临界时间步长 – 自然频率:
c ωmax= 2 l
• 临界时间步长:
其中
c=
E ρ
(波传播速度)
l Δ t= c
– Courant-Friedrichs-Levy-准则
Ls-dyna单元简介
– 如果总的沙漏能大于模型内能的 10%, 这个分析就有可能是失败的. 在 Chapter 8 (GLSTAT and MATSUM files)中将讨论沙漏能. 有时侯即使 5% 也是不允许的.
001322 10 JAN 2000 2-6
Training Manual
Explicit Dynamics with ANSYS/LSDYNA
•
所有显式动力单元为三维的
001322 10 JAN 2000 2-3
显式动力单元概述（续）
显式单元族在以下方面与 ANSYS 隐式单元明显不同:
Training Manual
Explicit Dynamics with ANSYS/LSDYNA
• 3D 圆杆单元用来承受轴向载 荷.
• 用3个节点定义单元.
– 第3个节点用来定义杆的初 始方向.
001322 10 JAN 2000 2-12
BEAM161 梁单元
• 由于不产生应变，此 3D 梁适用于刚体旋转 . 用3个节点定义此单元.
Training Manual
Explicit Dynamics with ANSYS/LSDYNA
这种算法与 SHELL143的算法兼容
001322 10 JAN 2000 2-15
SHELL163 薄壳单元（续）
•
Training Manual
Explicit Dynamics with ANSYS/LSDYNA
单元算法 BT, BWC, BL 只用面内单点积分，S/R CHL 用 4点积分.
Training Manual
001322 10 JAN 2000 2-11
ansys土木计算过程中常用的单元和材料类型!
ansys土木计算过程中常用的单元和材料类型!土木计算过程中常用的单元和材料类型一、单元（1）link(杆)系列：link1(2D)和link8(3D)用来模拟桁架，注意一根杆划一个单元。
link10用来模拟拉索，注意要加初应变，一根索可多分单元。
link180是link10的加强版，一般用来模拟拉索。
（2）beam(梁)系列：beam3(2D)和beam4(3D)是经典欧拉梁单元，用来模拟框架中的梁柱，画弯据图用etab读入smisc数据然后用plls命令。
注意：虽然一根梁只划一个单元在单元两端也能得到正确的弯矩图，但是要得到和结构力学书上的弯据图差不多的结果还需多分几段。
该单元需要手工在实常数中输入Iyy和Izz，注意方向。
beam44适合模拟薄壁的钢结构构件或者变截面的构件，可用"/eshape,1"显示单元形状。
beam188和beam189号称超级梁单元，基于铁木辛科梁理论，有诸多优点：考虑剪切变形的影响，截面可设置多种材料，可用"/eshape,1"显示形状，截面惯性矩不用自己计算而只需输入截面特征，可以考虑扭转效应，可以变截面（8.0以后），可以方便地把两个单元连接处变成铰接（8.0以后，用ENDRELEASE命令）。
缺点是：8.0版本之前beam188用的是一次形函数，其精度远低于beam4等单元，一根梁必须多分几个单元。
8.0之后可设置“KEYOPT(3)=2”变成二次形函数，解决了这个问题。
可见188单元已经很完善，建议使用。
beam189与beam188的区别是有3个结点，8.0版之前比beam188精度高，但因此建模较麻烦，8.0版之后已无优势。
(3)shell(板壳)系列shell41一般用来模拟膜。
shell63可针对一般的板壳，注意仅限弹性分析。
它的塑性版本是shell43。
加强版是shell181（注意18*系列单元都是ansys后开发的单元，考虑了以前单元的优点和缺陷，因而更完善），优点是：能实现shell41、shell63、shell43...的所有功能并比它们做的更好，偏置中点很方便（比如模拟梁版结构时常要把板中面望上偏置），可以分层，等等。
ansys-lsdyna知识
一、关于dyna中材料失效准则的定义有些材料类型中有关于失效准则的定义，但是也有些材料类型没有失效准则的材料类型，这时需要额外的失效准则定义，与材料参数一块定义材料特性。
需要用到*mat_add_erosion关键字，对于这个关键字有几个需要注意的地方。
1、材料的通用性破坏准则：`材料通常为拉破坏或者剪切破坏，静水压是以压为正，拉为负，所以静水压破坏就是给出最小的承受压力，当然需要小于0（即拉力），如果静水压小于该值，则材料破坏。
相反，应力则是以压为负，拉为正，故最大主应力或最大等效应力或最大剪应力破坏等等都是给出最大的应力极限，当然大于0，如果拉应力大于该值，则材料破坏，无论是*MAT_ADD_EROSION，还是材料内部自带的破坏准则还是其他软件，都遵循以上准则。
注意：屈服不是失效。
2、单元失效模拟的功能与目的单元删除功能是为了克服有限元本身的缺陷而提出的一项方法，由于有限元本身是基于连续介质力学的，而在连续介质力学中，所研究的物体需要是连续的，既物质域在空间中连续。
在这样的理论假设框架下，单元本身是不会消失的。
然而在实际情况下，由于损伤断裂的存在，势必会使得一些单元消失或者完全的失效，所以为了能够模拟这种情况，DYNA 提供了单元失效功能。
破坏、失效、断裂，都是工程性的概念，它表示在达到某一准则后，结构、构件、或者构件中的某一部分，从结构中退出工作，不再影响整体结构的受力。
而从有限元概念上说，对上述机制的模拟，基本手段都是一样的，就是当满足某一指标（比如某个应变大小）后，将一个单元或者一个积分点的质量、刚度和应力、应变都设为零（或者非常接近与零），这样它在整体结构计算中就不再发挥作用，进而实现了退出工作机制的模拟。
所以，无论是把纤维模型中的某个纤维、或者分层壳模型中的某一层、或者实体模型中的某个积分点，或者结构中的某个单元，让其不再参与整体结构计算，都可以达到模拟破坏退出工作的目的。
而所谓单元生死技术，是上述基本概念在有限元程序中的一个“打包”应用。
LS-DYNA 2D金属切削模拟步骤
在ANSYS Launcher界面中，选择ANSYS Mechanical/LS-DYNA1、菜单过滤Main Menu→Preprocessor→LD-DYNA Explicit→OK2、设置文件名及分析标题Utility Menu→File→change Jobname→2D cutting→New log and error file ：YES→OKUtility Menu→File→change Title→cutting analysis →OK3、选择单元类型Main menu→preprocessor→Element Type→Add/Edit/Delete→Add→2D solid 162→OK→options→选择const.stress ；Lagrangian→OK4、定义材料模型（1）定义刀具材料模型Main menu→preprocessor→Material Props→Material Models→rigid material→输入：DENS：5.2e3 ；EX：4.1e11 ；NUXY：0.3 ；选择“Y and Zdisps” ；“All rota tions”→OK（2）定义工件Johnson-cook材料模型Main menu→preprocessor→Material Props→Material Models→Gruneisen→Johnson-cook→输入：DENS：7.8e3 ；EX：2.06e11 ；NUXY：0.3A：507；B：320；C：0.28；n；0.064；m=1.06D1：0.15；D2：0.72；D3：1.66；D4：0.005；D5：--0.84yangmeng112010-8-30 17:43:435、创建几何模型（1）创建工件模型Main menu→preprocessor→Create→Areas→Rectangle→By Dimensions→输入：X1，X2：0，5；Y1，Y2：0，3→OK（2）创建刀片模型Main menu→preprocessor→Create→Keypionts→In Active CS→依次输入：keypoint number：5，X、Y、Z ：5.1，2.9，0；keypoint number：6，X、Y、Z ：6，3.228，0；keypoint number：7，X、Y、Z ：6，4，0；keypoint number：8，X、Y、Z ：5.294，4，0→OKyangmeng112010-8-30 17:44:006、网格划分（一）（1）对刀片进行网格划分Utility Menu→Select→Entities→Lines ：By Num/Pick→Apply→选取刀片边线→OKMain menu→preprocessor→Meshing→Size contrls→Manualsize→Lines→All lines→NDIV：10→OK（2）对刀尖半圆进行网格划分Utility Menu→Select→Entities→Lines ：By Num/Pick→Apply→选取刀尖半圆→OKMain menu→preprocessor→Meshing→Size contrls→Manualsize→Lines→All lines→NDIV：3→OK（3）确定刀片的单元属性Main menu→preprocessor→Meshing→Mesh Attributes→Picked Aeras→选取刀片→Apply→确定材料号和单元类型号为1→OK（4）刀片网格划分Main menu→preprocessor→Meshing→MeshTool→Mesh：Aeras；shape：Tri；free→Mesh→选取刀片→OK（二）（5）对工件进行网格划分切分工件Utility menu →Workplane→Wp settings→Grid and Triad→Minimum ,maximum:-5,5 ; Spacing:1.0→OK平移和旋转工作平面并用其切分工件Utility menu →Workplane→Offset wp by increm ens→X，Y，Z offsets：0，2.5，0；XY，YZ，ZX angle：0，90，0→OKMain menu→preprocessor→Modeling→operate→Booleans→Divide→Areas by wkplane→选取工件→OK取消工作平面显示Utility menu→workplane→Display workingplane→等分接触区域相关Y向线段Utility Menu→Select→Entities→Lines ：By Num/Pick→Apply→选取工件接触区Y向线段（两条）→OKMain menu→preprocessor→Meshing→Size contrls→Manualsize→Lines→All lines→NDIV：10→OK等分接触区域相关X向线段Utility Menu→Select→Entities→Lines ：By Num/Pick→Apply→选取工件接触区X向线段（两条）→OKMain menu→preprocessor→Meshing→Size contrls→Manualsize→Lines→All lines→NDIV：40→OK等分接触区域不相关Y向线段Utility Menu→Select→Entities→Lines ：By Num/Pick→Apply→选取工件接触区Y向线段（两条）→OKMain menu→preprocessor→Meshing→Size contrls→Manualsize→Lines→All lines→NDIV：25→OK等分接触区域不相关X向线段Utility Menu→Select→Entities→Lines ：By Num/Pick→Apply→选取工件接触区X向线段（底边）→OKMain menu→preprocessor→Meshing→Size contrls→Manualsize→Lines→All lines→NDIV：30→OK确定工件的单元属性Main menu→preprocessor→Meshing→Mesh Attributes→Picked Aeras→选取工件→Apply→确定材料号为2和单元类型号为1→OK工件网格划分Main menu→preprocessor→Meshing→MeshTool→Mesh：Aeras；shape：Quad；mapped→Mesh →选取工件→OKyangmeng112010-8-30 17:44:227、建立partMain menu→preprocessor→LS-DYNA options→part options→create all part→OK(part1:刀具；part2：工件)Plot→parts（不同颜色显示单元）8、定义接触信息Main menu→preprocessor→LS-DYNA options→contact→Define contact→surface to surf；Eroding；静、动摩擦系数为0.15、0.10→OK→弹出contact options对话框，确定接触件（工件），目标件（刀片）→OK9、施加边界条件Utility menu→select→Entities→Nodes :By Location :X Coordinates→Min,Max:-0.01,0.01;From Full→Apply（选中左侧边所有节点）Main menu→preprocessor→LS-DYNA options→Constraints→Apply→on nodes→pick All→All DOF→OKUtility menu→select→Entities→Nodes :By Location :Y Coordinates→Min,Max:-0.01,0.01;From Full→Apply（选中底边所有节点）Main menu→preprocessor→LS-DYNA options→Constraints→Apply→on nodes→pick All→All DOF→OK恢复整个模型的选择Utility menu→select→Everything10、对刀片施加初速度Main menu→preprocessor→LS-DYNA options→Initial Velocity→on parts→w/Nodal Rotate→选择part1，VX：-100→OK恢复整个模型的选择Utility menu→select→Everything11、设置能量控制选项Main menu→Solution→Analysis options→Energy options→打开所有能量控制选项→OK12、设置人工体积粘性选项Main menu→Solution→Analysis options→Bulks viscosity→Quadratic Viscosity Coefficient:1.0→OK13、设置时间步长因子Main menu→Solution→Time controls→Time step ctrls→Time step scale factor:0.6→OK14、设置求解时间Main menu→Solution→Time controls→Solution time→1e-3→OK15、设置结果文件输出步数Main menu→Solution→Output Controls→File output Freq→Number of steps→[EDRST]:50;[EDHTIME]:50→OK16、设置结果文件的输出类型Main menu→Solution→Output Controls→Output File Types→Add：ANSYS and LS-DYNA→OK 17、输出K文件Main menu→Solution→Write jobname.K18、求解Main menu→Solution→Solve19、后处理（暂时不管）。
lsdyna简介
LS-DYNA 简介ANSYS学习2009-02-17 20:03:54 阅读444 评论0 字号：大中小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的最新版本是2004年8月推出的970版。
ANSYS/LS-DYNA的前后处理器是ANSYS/PRE-POST,求解器LS-DYNA，是全世界范围内最知名的有限元显式求解程序。
LS-DYNA在1976年由美国劳伦斯·利沃莫尔国家实验室（Lawrence Livermore National Laboratory）J.O.Hallquist博士主持开发，时间积分采用中心差分格式，当时主要用于求解三维非弹性结构在高速碰撞、爆炸冲击下的大变形动力响应，是北约组织武器结构设计的分析工具。
LS-DYNA 的源程序曾在北约的局域网Pubic Domain公开发行，因此在广泛传播到世界各地的研究机构和大学。
从理论和算法而言，LS-DYNA是目前所有的显式求解程序的鼻祖和理论基础。
1988年，J.O.Hallquist创建利沃莫尔软件技术公司（Livermore Software Technology Corporation），LS-DYNA开始商业化进程，总体来看，到目前为止在单元技术、材料模式、接触算法以及多场耦合方面获得非常大的进步。
ls-dyna单元简介
LS-DYNA使用指南中文版本第2章单元(2012-04-05 11:25:53)转载▼标签：杂谈第二章单元在显式动态分析中可以使用下列单元：·LINK160杆·BEAM161梁·PLANE162平面·SHELL163壳·SOLID164实体·COMBI165弹簧阻尼·MASS166质量·LINK167仅拉伸杆本章将概括介绍各种单元特性,并列出各种单元能够使用的材料类型。
除了PLANE162之外,以上讲述的显式动态单元都是三维的,缺省时为缩减积分(注意：对于质量单元或杆单元缩减积分不是缺省值)缩减积分意味着单元计算过程中积分点数比精确积分所要求的积分点数少。
因此,实体单元和壳体单元的缺省算法采用单点积分。
当然,这两种单元也可以采用全积分算法。
详细信息参见第九章沙漏,也可参见《LS-DYNA Theoretical Manual》。
这些单元采用线性位移函数；不能使用二次位移函数的高阶单元。
因此，显式动态单元中不能使用附加形状函数，中节点或P-单元。
线位移函数和单积分点的显式动态单元能很好地用于大变形和材料失效等非线性问题。
值得注意的是，显单元不直接和材料性能相联系。
例如，SOLID164单元可支持20多种材料模型，其中包括弹性，塑性，橡胶，泡沫模型等。
如果没有特别指出的话（参见第六章，接触表面），所有单元所需的最少材料参数为密度，泊松比，弹性模量。
参看第七章材料模型，可以得到显式动态分析中所用材料特性的详细资料。
也可参看《ANSYS Element Reference》,它对每种单元作了详细的描述，包括单元的输入输出特性。
2.1实体单元和壳单元2.1.1SOLID164SOLID164单元是一种8节点实体单元。
缺省时，它应用缩减（单点）积分和粘性沙漏控制以得到较快的单元算法。
单点积分的优点是省时，并且适用于大变形的情况下。
lsdyna使用方法2
第1章ANSYS／LS-DYNA基石出知识有限元2009-05-12 20:06:17 阅读62 评论0 字号：大中小订阅近年来，非线性结构动力仿真分析方面的研究工作和工程应用取得了很大的发展。
20世纪90年代中后期，著名的通用显式动力分析程序LS-DYNA被引入中国，在相关的工程领域中迅速得到广泛的应用，已成为国内科研人员开展数值实验的有力上具。
LS-DYNA的显式算法特别适合于分析各种非线性结构冲击动力学问题，如爆炸，结构碰撞、金属加工成形等高度非线性的问题，同时还可以求解传热、流体以及流固耦合问题。
LSTC公司和ANSYS公司合作推出的ANSYS／LS-DYNA软件，结合了LS-DYNA强大的显式动力分析方法与ANSYS的前后处理功能。
对于曾经接触过ANSYS结构分析的读者而言，谊程序无疑是最理想的辅助动力分析工具。
本章的目的在于全面介绍ANSYS／LS-DYNA的基础知识，包括下面的几个主题：+ LS-DYNA计算程序的发展过程☆LS-DYNA的分析功能与应用范围+ ANSYS／LS-DYNA的工作环境+ ANSYS／LS-DYNA的一般分析过程+ ANSYS~S-DYNA的程序组织和丈件系统+ LS-DYNA显式动力分析的基本概念1．1 LS—DYNA计算程序的发展过程1976年，美国LawrenceLivermore国家实验室J．O．Hallquist博士主持开发完成了DYNA程序系列，主要目的是为武器设计提供分析工具。
1986年部分DYNA源程序在Public Domain(北约局域网)发布，从此在研究和教育机构广泛传播，被公认为是显式有限元程序的先导，是目前所有显式求解程序的基础代码。
1988年，J．O．Hallquist创建LSTC公司(LivermoreSoftwarenchnolOWCorporation)，推出LS-DYNA 程序系列，主要包括显式LS-DYNA2D、LS-DYNA3D，隐式LS-NIKE2D、LS-NIKE3D、热分析LS-TOPAZ2D、LS-TOPAZ3D，前后处理LS-MAZE、LS-ORION、LS-INGRID、LS-TAURUS等商用程序，逐步规范和完善程序的分析功能，陆续推出930版(1993年)和936版(1995年)，同时增加了汽车安全性分析、金属板的冲压成形以及流固耦合(ALE算法和Eluer算法)，使得LS-DYNA程序系统的应用范围不断扩大，并建立起完备的软件质量保证体系。
lsdyna使用方法2
第1章ANSYS／LS-DYNA基石出知识有限元2009-05-12 20:06:17 阅读62 评论0 字号：大中小订阅近年来，非线性结构动力仿真分析方面的研究工作和工程应用取得了很大的发展。
20世纪90年代中后期，著名的通用显式动力分析程序LS-DYNA被引入中国，在相关的工程领域中迅速得到广泛的应用，已成为国内科研人员开展数值实验的有力上具。
LS-DYNA的显式算法特别适合于分析各种非线性结构冲击动力学问题，如爆炸，结构碰撞、金属加工成形等高度非线性的问题，同时还可以求解传热、流体以及流固耦合问题。
LSTC公司和ANSYS公司合作推出的ANSYS／LS-DYNA软件，结合了LS-DYNA强大的显式动力分析方法与ANSYS的前后处理功能。
对于曾经接触过ANSYS结构分析的读者而言，谊程序无疑是最理想的辅助动力分析工具。
本章的目的在于全面介绍ANSYS／LS-DYNA的基础知识，包括下面的几个主题：+ LS-DYNA计算程序的发展过程☆LS-DYNA的分析功能与应用范围+ ANSYS／LS-DYNA的工作环境+ ANSYS／LS-DYNA的一般分析过程+ ANSYS~S-DYNA的程序组织和丈件系统+ LS-DYNA显式动力分析的基本概念1．1 LS—DYNA计算程序的发展过程1976年，美国LawrenceLivermore国家实验室J．O．Hallquist博士主持开发完成了DYNA程序系列，主要目的是为武器设计提供分析工具。
1986年部分DYNA源程序在Public Domain(北约局域网)发布，从此在研究和教育机构广泛传播，被公认为是显式有限元程序的先导，是目前所有显式求解程序的基础代码。
1988年，J．O．Hallquist创建LSTC公司(LivermoreSoftwarenchnolOWCorporation)，推出LS-DYNA 程序系列，主要包括显式LS-DYNA2D、LS-DYNA3D，隐式LS-NIKE2D、LS-NIKE3D、热分析LS-TOPAZ2D、LS-TOPAZ3D，前后处理LS-MAZE、LS-ORION、LS-INGRID、LS-TAURUS等商用程序，逐步规范和完善程序的分析功能，陆续推出930版(1993年)和936版(1995年)，同时增加了汽车安全性分析、金属板的冲压成形以及流固耦合(ALE算法和Eluer算法)，使得LS-DYNA程序系统的应用范围不断扩大，并建立起完备的软件质量保证体系。
ansys中LS-DYNA 2D金属切削模拟步骤
在ANSYS Launcher界面中，选择ANSYS Mechanical/LS-DYNA1、菜单过滤Main Menu→Preprocessor→LD-DYNA Explicit→OK2、设置文件名及分析标题Utility Menu→File→change Jobname→2D cutting→New log and error file ：YES→OK Utility Menu→File→change Title→cutting analysis →OK3、选择单元类型Main menu→preprocessor→Element Type→Add/Edit/Delete→Add→2D solid 162→OK→options→选择const.stress ；Lagrangian→OK4、定义材料模型（1）定义刀具材料模型Main menu→preprocessor→Material Props→Material Models→rigid material→输入：DENS：5.2e3 ；EX：4.1e11 ；NUXY：0.3 ；选择“Y and Zdisps” ；“All rotati ons”→OK（2）定义工件Johnson-cook材料模型Main menu→preprocessor→Material Props→Material Models→Gruneisen→Johns on-cook→输入：DENS：7.8e3 ；EX：2.06e11 ；NUXY：0.3A：507；B：320；C：0.28；n；0.064；m=1.06D1：0.15；D2：0.72；D3：1.66；D4：0.005；D5：--0.845、创建几何模型（1）创建工件模型Main menu→preprocessor→Create→Areas→Rectangle→By Dimensions→输入：X1，X2：0，5；Y1，Y2：0，3→OK（2）创建刀片模型Main menu→preprocessor→Create→Keypionts→In Active CS→依次输入：keypoint number：5，X、Y、Z ：5.1，2.9，0；keypoint number：6，X、Y、Z ：6，3.228，0；keypoint number：7，X、Y、Z ：6，4，0；keypoint number：8，X、Y、Z ：5.294，4，0→OK6、网格划分（一）（1）对刀片进行网格划分Utility Menu→Select→Entities→Lines ：By Num/Pick→Apply→选取刀片边线→O KMain menu→preprocessor→Meshing→Size contrls→Manualsize→Lines→All lin es→NDIV：10→OK（2）对刀尖半圆进行网格划分Utility Menu→Select→Entities→Lines ：By Nu m/Pick→Apply→选取刀尖半圆→O KMain menu→preprocessor→Meshing→Size contrls→Manualsize→Lines→All lin es→NDIV：3→OK（3）确定刀片的单元属性Main menu→preprocessor→Meshing→Mesh Attributes→Picked Aeras→选取刀片→Apply→确定材料号和单元类型号为1→OK（4）刀片网格划分Main menu→preprocessor→Meshing→MeshTool→Mesh：Aeras；shape：Tri；free →Mesh→选取刀片→OK（二）（5）对工件进行网格划分切分工件Utility menu →Workplane→Wp settings→Grid and Triad→Minimum ,maximum:-5, 5 ;Spacing:1.0→OK平移和旋转工作平面并用其切分工件Utility menu →Workplane→Offset wp by incremens→X，Y，Z offsets：0，2.5，0；XY，YZ，ZX angle：0，90，0→OKMain menu→preprocessor→Modeling→operate→Booleans→Divide→Areas by wkp lane→选取工件→OK取消工作平面显示Utility menu→workplane→Display workingplane→等分接触区域相关Y向线段Utility Menu→Select→Entities→Lines ：By Num/Pick→Apply→选取工件接触区Y向线段（两条）→OKMain menu→preprocessor→Meshing→Size contrls→Manualsize→Lines→All lin es→NDIV：10→OK等分接触区域相关X向线段Utility Menu→Select→Entities→Lines ：By Num/Pick→Apply→选取工件接触区X向线段（两条）→OKMain menu→preprocessor→Meshing→Size contrls→Manualsize→Lines→All lin es→NDIV：40→OK等分接触区域不相关Y向线段Utility Menu→Select→Entities→Lines ：By Num/Pick→Apply→选取工件接触区Y向线段（两条）→OKMain menu→preprocessor→Meshing→Size contrls→Manualsize→Lines→All lin es→NDIV：25→OK等分接触区域不相关X向线段Utility Menu→Select→Entities→Lines ：By Num/Pick→Apply→选取工件接触区X向线段（底边）→OKMain menu→preprocessor→Meshing→Size contrls→Manualsize→Lines→All lin es→NDIV：30→OK确定工件的单元属性Main menu→preprocessor→Meshing→Mesh Attributes→Picked Aeras→选取工件→Apply→确定材料号为2和单元类型号为1→OK工件网格划分Main menu→preprocessor→Meshing→MeshTool→Mesh：Aeras；shape：Quad；mapp ed→Mesh→选取工件→OK7、建立partMain menu→preprocessor→LS-DYNA options→part options→create all part→O K(part1:刀具；part2：工件)Plot→parts（不同颜色显示单元）8、定义接触信息Main menu→preprocessor→LS-DYNA options→contact→Define contact→surface to surf；Eroding；静、动摩擦系数为0.15、0.10→OK→弹出contact options对话框，确定接触件（工件），目标件（刀片）→OK9、施加边界条件Utility menu→select→Entities→Nodes :By Location :X Coordinates→Min,Max: -0.01,0.01;From Full→Apply（选中左侧边所有节点）Main menu→preprocessor→LS-DYNA options→Constraints→Apply→on nodes→pi ck All→All DOF→OKUtility menu→select→Entities→Nodes :By Location :Y Coordinates→Min,Max: -0.01,0.01;From Full→Apply（选中底边所有节点）Main menu→preprocessor→LS-DYNA options→Constraints→Apply→on nodes→pi ck All→All DOF→OK恢复整个模型的选择Utility menu→select→Everything10、对刀片施加初速度Main menu→preprocessor→LS-DYNA o ptions→Initial Velocity→on parts→w/No dal Rotate→选择part1，VX：-100→OK恢复整个模型的选择Utility menu→select→Everything11、设置能量控制选项Main menu→Solution→Analysis options→Energy options→打开所有能量控制选项→OK12、设置人工体积粘性选项Main menu→Solution→Analysis options→Bulks viscosity→Quadratic Viscosit y Coefficient:1.0→OK13、设置时间步长因子Main menu→Solution→Time controls→Time step ctrls→Time step scale facto r:0.6→OK14、设置求解时间Main menu→Solution→Time controls→Solution time→1e-3→OK15、设置结果文件输出步数Main menu→Solution→Output Controls→File output Freq→Number of steps→[EDRST]:50;[EDHTIME]:50→OK16、设置结果文件的输出类型Main menu→Solution→Output Controls→Output File Types→Add：ANSYS and LS -DYNA→OK17、输出K文件Main menu→Solution→Write jobname.K18、求解Main menu→Solution→Solve19、后处理（暂时不管。
dyna材料定义
TB ，EVISC
TBDATA ，46，
TBDATA ，47，
TBDATA ，48，K
TBDATA ，61，1/β
注 --对于这种材料选项，必须用 MP 命令定义密度（DENS）。
例题参看B.2.6, Viscoelastic Example：Glass。
7.2.3 非线性无弹性模型
7.2.2. 3 粘弹性模型
Herrmann和Peterson提出的线性粘弹性模型。模型采用偏量特性：
这里剪切松弛模量由下式给出：
在模型中，由体积V计算增量积分压力时，需事先进行弹性体积假设，即V：p=K 。用参数 、 、K（体积模量）和β来定义线粘弹性模型。用 TB ，EVISC和 TBDATA 命令的46、47、48和61项输入以上数据：
通过GUI路径定义材料模型比使用命令直接得多：
1．选择菜单路径Main Menu>Preprocessor>Material Props>Material Models.Define Material Model Behavior对话框出现。
注 --如果不事先定义ANSYS/LS-DYNA单元类型，那么就不能定义ANSYS/LS-DYNA材料模型。
·横向各向异性弹塑性（#37）
·横向各向异性FLD（#39）
·随动双线性（#3）
·随动塑性（#3）
·3参数Barlat（#36）
·Barlat各向异性塑性（#33）
·与应变率相关的幂函数塑性（#64）
·应变率相关塑性（#19）
·复合材料破坏（#22）
·混凝土破坏（#72）
用 MP 命令输入泊松比（ ）和密度。（泊松比的值要比推荐的大一些，太小的值不能工作。）用 TB 和 TBDATA 命令输入Mooney-Rivlin常数，只允许一种温度下的数据，并且必须放在数据表中的1和2位置。
LS-DYNA第七章：材料模型
第七章材料模型ANSYS/LS-DYN包括40多种材料模型，它们可以表示广泛的材料特性，可用材料如下所示。
本章后面将详细叙述材料模型和使用步骤。
对于每种材料模型的详细信息，请参看Appendix B,Material Model Examples 或《LS/DYNA Theoretical Ma nual》的第十六章（括号内将列出与每种模型相对应的LS-DYNA 材料号）。
线弹性模型•各向同性（#1）•正交各向异性（#2）•各向异性（#2）•弹性流体（#1）非线弹性模型•Blatz -ko Rubber倂7）•Moo neyRivlin Rubber（#27）•粘弹性（#6）非线性无弹性模型•双线性各向同性（#3）•与温度有关的双线性各向同性（#4）•横向各向异性弹塑性（#37）•横向各向异性FLD（#39）•随动双线性（#3）•随动塑性（#3）•3参数Barlat （#36）• Barlat 各向异性塑性（#33）•与应变率相关的幂函数塑性（#64）•应变率相关塑性（#19）•复合材料破坏（#22）•混凝土破坏（#72）•分段线性塑性（#24）•幂函数塑性（#18）压力相关塑性模型•弹-塑性流体动力学（#10）•地质帽盖材料模型（#25）泡沫模型•闭合多孔泡沫（#53）•粘性泡沫（#62）•低密度泡沫（#57）•可压缩泡沫（#63）•Ho neycomb（#26）需要状态方程的模型•Bamma塑性（#51）• Johnson-Cook 塑性（#15）•空材料（#9）•Zeril li-Armstrong（#65）•Stei nberg （# 11）离散单元模型•线弹性弹簧•普通非线性弹簧•非线性弹性弹簧•弹塑性弹簧•非弹性拉伸或仅压缩弹簧•麦克斯韦粘性弹簧•线粘性阻尼器•非线粘性阻尼器•索（缆）（#71）刚性体模型•刚体（#20）7.1定义显示动态材料模型用户可以采用ANSY命令MP , MPTEMP, MPDATA, TB , TBTEMP和TBDATA以及ANSYS/LS-DYN命令EDMP来定义材料模型。
ANSYS_LSDYNA
LS-DYNA 程序采用动力松弛技术，可以进行动力分析前的预应力计算或者进行静力分析。
LS-DYNA3D 程序的算法基础： （1）控制方程1.动量方程：i i j j x f ..,i ρρσ=+ij σ为柯西应力 i f 为单位质量体积力 ..x 为加速度2.质量守恒：0ρρJ =ρ为当前质量密度 0ρ为初始质量密度 J 为体积变化率3.能量方程V )(VS E ij ij q p +-=。
ε………………………………………………………….(2.5)用于状态方程计算 和 总的能量平衡。
式中，V 为现时构形的体积；ij 。
ε为应变率张量；q 为体积粘性阻力 偏应力ij ij ij q p S σσ)(++= 压力q p kk --=σ31力分量的边界值就等于对应的面力分量，如下式表示）)(t t n i j ij =σ在1S 面力边界上式中：j n ，j=1,2,3为现时构形边界1S 的外法线方向余弦；i t ，i=1,2,3为面力载荷b.位移边界条件（在位移边界问题中，物体在全部或部分边界u s 上的位移分量都是已知函数，即在边界上有：)()(s u u s = )()(s v v s = )()(s w w s =，其中s u )(、s v )(、s w )(是位移函数在边界上的值，)(s u 、)(s v 、)(s w 表示边界上的已知位移分量，例如对于完全固定边界0===w v u ，有0)(=s u ，0)(=s v ，0)(=s w ）)(),(t K t X x i j i =在2S 位移边界上式中，)(t K i ，i=1,2,3是给定位移函数 c.滑动接触面间断处的跳跃条件单元内某一点的自然坐标表示成e x N t x }]{[)},,,({=ςηξ式中，单元内任意点在t 时刻的坐标矢量Tt x )},,,({ςηξ=[1x 2x 3x ]单元节点坐标矢量],,,...,,[}{838281131211x x x x x x x eT =插值矩阵2438181810 (0000) (0000)...00)],,([⨯⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=φφφφφφςηξNLS-DYNA3D 程序将单元质量矩阵⎰=mV TNdV N m ρ的同一行矩阵元素都合并到对角元素项，形成集中质量矩阵。
LS-DYNA材料模型及参数
LS-DYNA材料模型及参数目录1基本的状态方程 (2)1.1EOS_JWL (2)2.2EOS_GRUNEISEN (2)2.3EOS_LINEAR_POL YNOMIAL (3)2.材料模型 (3)2.1MA T_HIGH_EXPLOSIVE_BURN (3)RDX (5)HMX (5)TNT (5)1.2MA T_NULL (5)空气 (6)水 (6)1.3MA T_JOHNSON_COOK (7)紫铜 (8)钢 (8)1.4 MA T_PLASTIC_KINEMATIC (9)钢 (10)高导无氧铜 (10)土壤 (10)1.5MA T_STEINBERG (10)高导无氧铜 (12)1.6MA T_JOHNSON_HOLMQUIST_CERAMICS (12) B4C陶瓷 (14)1.7MA T_JOHNSON_HOLMQUIST_CONCRETE (14) 混凝土 (14)3其它材料参数 (15)L Y12CZ铝合金 (15)主要材料模型及参数1基本的状态方程1.1EOS_JWL2.2EOS_GRUNEISEN2.3EOS_LINEAR_POLYNOMIAL (对EOS_GRUNEISEN进行线性化)2.材料模型2.1MAT_HIGH_EXPLOSIVE_BURNRDX密度：1.69E+3 kg/m3；D: 8310m/s；Pcj ：30.45 GpaA:850 Gpa；B: 18 Gpa；R1: 4.6；R2: 1.3；w0.38；E0:10MJ/kgFor（g-cm-us）:*MAT_HIGH_EXPLOSIVE_BURN1 1.69 8.310 0.3015 0*EOS_JWL1 8.50 0.18 4.6 1.3 0.38 10 e-02 1.00 HMX密度：1.891 E+3 kg/m3, D：9910m/s，Pcj：42Gpa,A：778.3 Gpa；B：7. 1 Gpa；R1：4.1；R2：1.00；w0：30；E0：10. 5 MJ/kgFor（g-cm-us）:*MAT_HIGH_EXPLOSIVE_BURN1 1.89 9.910 0.42 0*EOS_JWL1 7.783 0.071 4.2 1.0 0.30 10.5 e-02 1.00TNT密度：1.63 E+3 kg/m3；D：6930 m/s；Pcj：27 Gpa：A：371.2 Gpa；B：3.21 Gpa；R1：4.15；R2：0.95；w0：30 E0:4.29 MJ/kgFor（g-cm-us）:*MAT_HIGH_EXPLOSIVE_BURN1 1.63 6.930 0.27 0*EOS_JWL1 3.713 0.0743 4.15 0.95 0.30 7.0 e-02 1.001.2MAT_NULL空气*MAT_NULLRO=1.25 kg/m3, PC= -1.0pa（<0）, MU=1. 7456E-5（动力粘性系数）*EOS_LINEAR_POLYNOMIAL1 ，0 Gpa，0 Gpa，0 Gpa，0，0.4 ，0.4 ，0253312.5， 1.0*MAT_NULL3 0.125e-02 -1.0E-12 1.749E-7 00000 0000 00 00 *EOS_GRUNEISEN3 0.3444 00000 0000 00000 1.40 00 0000/*EOS_LINEAR_POL YNOMIAL3 0 0 0 0 0.4 0.4 02.5000E-6 1水*MAT_NULL1, RO=998.21 kg/m3，PC= -10.0 pa，MU=0.8684E-3, 0, 0, 0, 0C:1480m/s，S1: 2.56，S2 :-1.986，S3 :0.2268, γ:0.4934，A:0.47，E0:0V0:1For（g-cm-us）:*MAT_NULL1 0.998 -1.0E-11 0.8684E-5 00000 0000 00 00 *EOS_GRUNEISEN1 1.65 1.92 -0.096 00000 0.350 00 00001.3MAT_JOHNSON_COOK紫铜EX=1.19*MAT_JOHNSON_COOK1 8.96000 0.460.900E-03 2.920E-03 0.310 0.250E-01 1.09 0.1356E+04 210 0.100E-050.383E-05 -9.00E+00 3.00 0.00 3.00 0.00 0.00 0.000.00*EOS_GRUNEISEN1 0.394 1.489 0.00 0.00 2.02 0.47 0.001.00钢EX=2.0*MAT_JOHNSON_COOK2 9.96000 0.460.900E-03 2.920E-03 0.310 0.25E-01 1.09 0.136E+04 210 0.100E-05 0.383E-05 -9.00E+00 3.00 0.00 3.00 0.00 0.00 0.000.00*EOS_GRUNEISEN1 0.394 1.489 0.00 0.00 2.02 0.47 0.001.001.4 MAT_PLASTIC_KINEMATICj S e H钢*MAT_PLASTIC_KINEMA TIC4 7.83 2.07 0.300 0.400E-02 5.00E-02 1.000.00 0.00 0.00高导无氧铜*MAT_PLASTIC_KINEMA TIC1 8.93 1.17 0.350000 0.400E-02 0.100E-02 1.000.00 0.00 0.00土壤*MAT_PLASTIC_KINEMA TIC1 7.80 2.1 0.300 0.023 0. 0240. 1.000 0 01.5MAT_STEINBERG高导无氧铜*MAT_STEINBERG2 8.93 0.477 0.120E-02 36.0 0.450 0.00 0.640E-02 2.83 2.83 0.377E-03 0.100E-02 63.5 0.179E+04 2.02 1.50 -9.00 3.00 0.00 0.00 0.00 0.00 0.00 0.000.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00*EOS_GRUNEISEN2 0.394 1.49 0.00 0.00 2.02 0.470 0.001.001.6MAT_JOHNSON_HOLMQUIST_CERAMICSB4C陶瓷*MAT_JOHNSON_HOLMQUIST_CERAMICS2 2510e-6 197e+3 0.927 0.7 0.005 0.85 0.671.0 260 0.2 19e+3 8.71e+3 1.00.001 0.5 233e+3 -593e+3 28e+5 0.11.7MAT_JOHNSON_HOLMQUIST_CONCRETE 混凝土$*MAT_JOHNSON_HOLMQUIST_CONCRETE$ 1 0.311E+04 20.87E+9 0.79 1.60 0.007 0.61 72E6 $ 5.2E6 1E-6 0.01 7.0 16E6 0.001 8E8 0.10$ 0.04 1.0 85E9 -171E9 208E93其它材料参数LY12CZ铝合金抗拉强度MPa 5240.2%屈服强度MPa462伸长率％11体膨胀系数m3(.m3.K)－168*10-620°C体积电导率％IACS33 20°C电阻率N .M 52.2 强性模具E/GPa 71硬度HB140-150密度KG.m-32810。
ANSYS LS_DYNA
聚能装药战斗部设计
■ 隐式求解
用于非线性结构静动力分析，包括结构固有频率和 振型计算。ANSYS LS-DYNA中可以交替使用隐式求解和 显式求解，进行薄板冲压成形的回弹计算、结构动力分析 之前施加预应力等。
隐式、自适应算法用于球壳结构静力分析
ANSYS LS-DYNA
7
■ 前后处理功能
ANSYS LS-DYNA可以采用二种方式进行前处理： 第一、利用ANSYS 经典环境下的ANSYS PrepPost For LS-DYNA进行前处理，具有多种自动网格划分可选择，还能 使用LS-PrePost进行Keyword文件的编辑。 第二、ANSYS WorkBench2.0也可以为ANSYS LS-DYNA提供前处理功能，且可以直接导出k文件，方便用户提交给 LS-DYNA求解器进行计算。 ANSYS LS-DYNA求解结果的后处理可以在ANSYS 经典环境中进行，也可以采用LS-PrePost完成。
● 汽车 用ANSYS LS-DYNA的壳单元、实体单元构造完整的汽车模型，可以选用金属、玻璃、塑料、橡胶等各种材料模
型。程序的Automatic_Single_Surface接触功能可以保证汽车全部构件内外表面、与假人、气囊、安全带之间，以及外部 障碍物表面相互接触时不穿透，实现汽车高速碰撞时全过程的数值模拟和有关数据的输出。
● 土木工程抗震特性、在动载荷作用下的动态相应、爆破拆除、路桥防护设计设计、设施、设备的抗冲击设计、结构破 坏(裂纹)特性模拟、多体运动学等
● 国防安全设计 - 武器系统设计（内弹道/终点弹道效应、侵彻与开挖（穿甲、射流破甲、钻地效应）、介质中的爆炸场） - 战斗部结构设计（战斗部特征尺寸和起爆参数的优化设计、目标毁伤特性） - 防护加固系统设计（防护及加固系统的耐冲击、抗爆破特性）
LS_DYNA用户使用手册(中)
另外，可以用 EDCRB 命令把两个刚体结合在一起，它们的作用如同一个刚性体。但与刚性体 的定义不同的是，它主要依据材料号，根据 PARTＩD 定义的刚体约束和一个约束方程号。因此，为 了在两个物体间施加刚体约束，必须执行命令 EDCRB, ADD, NEQN, PARTS，NEQN 为约束方程参 考编号，PARTM 是主刚体 part 号，PARTS 是从刚体 part 号。必须注意不能用同一个 NEQN 值执行多 个 EDCRB 命令，仅使用最后一个 NEQN 值。当使用 EDCRB 命令时，第二个刚性体将被第一个刚 性体吸收，因而此后对第二个刚性体的任何操作将无效．
8.1 定义惯性特性
缺省时，程序将计算每个刚体的惯性特性。但是，用户可能会发现给刚性体设定重力、质量、
初始速度（在整体或局部坐标系）的特定中心和特定惯性张量是很有用的，而不是依赖于求解过程
中由有限元模型计算得到的值。可以用 EDIPART 来定义刚性体的任意特征，命令格式如下： 其中：PART－定义惯性的 part ID
8.3 变形体和刚性体部件间的转换
在一些动力学应用中，如果模型中大部分单元都是可变形的，那么，持续时间较长的刚体运动 的计算量就会极为庞大。摩托车轮子就是这样一个例子，轮子翻滚持续较长时间相对于后来的撞击 会占用更多的 CPU。为了提高这种应用的计算效率，ANSYS/LS-DYNA 提供了这样一种转换性能， 它把某些材料从变形状态换至刚性状态，然后转到变形状态，通过刚体运动过程中转换变形体为刚 性体，可以节省大量的 CPU 时间。
ANSYS-LS-DYNA-显式动力学-培训手册-第三天
材料的定义 ... 线弹性材料
各向同性 : 大多数工程材料（例如钢铁）都是各向同性的 通过DENS、EX和NUXY定义即可 正交各向异性 : EDMP, ORTHO 正交各向异性材料通过9个独立的常数和DENS定义 横向各向同性材料（正交各向异性的一种特例）通过5个独立的常数(EXX, EZZ, NUXY, NUXZ, GXY)和 DENS定义 正交各向异性材料的定义与EDLCS命令所定义的坐标系编号有关: Preprocessor: Material Props > Local CS > Create Local CS 各向异性 : EDMP, ORTHO 和 TB, ANEL 各向异性材料通过21个相互独立的参数和DENS定义 使用局部坐标系 (EDLCS) 和数据表 (TB,ANEL)
材料的定义 F. 塑性材料
塑性材料模型包含ANSYS/LS-DYNA中大多数非线性非弹性材料。要根据所分析材料的类型、应用领域和材料常数的可获取性来选择某个特定塑性模型。 塑性模型可以分为五类: 类别 1：应变率和塑性无关的各向同性材料 类别 2: 应变率和塑性相关的各向同性材料 类别 3: 应变率和塑性无关的各向异性材料 类别 4: 压力相关塑性 类别 5: 温度敏感塑性 为分析材料选择正确的类别非常重要， 在某类别内选择特定的模型相比就不那么重要了，这通常取决于材料数据的可获取性。