LS-DYNA牋notes牋牋?

1）ls－dyna区分零约束和非零约束，非零约束如载荷那样处理。如果是零约束才可以使用Constraint命令！

2）ls－dyna把所有重要的信息都写入到操作窗口和文件d3hsp中

3）输出ls－dyna后处理文件关键字$
*DATABASE_FORMAT
0

4)刚体的约束通过譓AT_RIGID材料定义中的相关参数来设置平动约束和转动约束，施加在质心上

5）两个刚体通过*CONSTRAINED_RIGID_BODIES进行合并

6）简单重启动，在dos命令行输入LS_DYNA R=D3DUMPnn

7）自动生成重启动文件

*DATABASE_BINARY_D3DUMP

$ dt/cycl

50000

8)考虑前面分析的的变形和应力应变，需要设置关键字*STRESS_INITIALIZATION_OPTION来实现

9）在ansys中定义材料模型，首先需要定义两个数组来输入应力和相应的应变，然后生成一条应力应变曲线，最好该材料模型调用该应力应变曲线。 两个数组：应变数组和应力数组。LS-DYNA Option/Loading optons/curve options/add curve 把横坐标和纵坐标选为应变和应力数组。（是真实应力应变曲线）

10）LS-DYNA不可压缩流体的显式算法（一阶精度）不用耦合的方式求解动量方程，从而减少了内存的需求，但损失一部分精度。另一方面，由于显示算法满足扩散和CFD稳定性条件，又可提供精度。因此，不可压缩流体求解分析中采用单点积分和沙漏稳定性的显示算法被证明式简便、高效的。

11）LS-DYNA不可压缩流体求解器的二阶精度算法采用了恒定质量的预置算法和物质质量的校正算法，合理解耦了速度场和压力场，从而减少N-S方程对CPU和内存的需求。二阶精度算法可用于分析流场中的涡流，且很容易推广到流体力学领域中的湍流现象的计算分析。

12）ALE算法：ALE算法先执行一个或者几个Langrange时步计算，此时单元网格随材料流动而发生了变形，然后执行ALE时步计算（1）保持变形后的物体边界条件，对内部单元进行重划分网格，网格拓扑关系不变，称为Smooth Step；（2）将变形网格中的单元变量（密度、能量、应变能量等）和节点速度矢量输运到重分后的新网格中，称为Advection Step。用于可以选择ALE时步的开始和终止时间，以及其频率。

13）LS-DYNA进行流固耦合有两种方法：

一种方法是共节点，不需要*CONSTRAINED_LAGRANGE_IN_SOLID关键字定义流体和固体的耦合。实际上该种方法相当于加了一个固体边界条件，不能处理相互作用问题，一般不太使用。

另一种方法是使用*CONSTRAINED_LAGRANGE_IN_SOLID把流体与固体单元耦合在一起。

在所有的流固耦分析中，建立模型时需注意，流体网格和固体网格是重叠的，并不需要网格一致，这给模型的建立带来了很大的方便。

14)LS-DYNA三维单元的基本算法Lagrange、Euler和ALE算法，是由关键字*SECTION_SOLID中的ELFORM控制

15>用全积分单元可以很大程度上减少沙漏。全积分在计算动力问题的时候还是有很多的缺点的，所以一般是选用单点积分的方式，因为此而引起的沙漏问题，dyna提供了多种算法可以减少到内能的5％以下，

16）控制沙漏：首先：您的模型如果是很多PART组成的装配体，那么您需要找出沙漏能最大的PART！
然后：才能进行一些适当的控制！第一：还得从载荷谈起，避免集中载荷；
第二：从网格谈起，尽量做到网格协调；第三：从沙漏控制谈起，采用不同的控制方法；第四：从单元算法谈起，不行就将一个PART化分为多个PART，除了单元算法不同以外，其他相同，采用全积分！ 这年头谁还用软件抓图？都是抓屏！方便，快捷，省时，省力！反正还得需要其他软件进行图片处理！

17>有限元方法一般以节点的位移作为基本变量，单元内各点的位移以及应变均采用形函数对各节点的位移进行插值计算而得，应力根据本构方程由应变计算得到，然后就可以计算单元的内能了。如果采用单点积分（积分点在等参元中心），在某些情况下节点位移不为零（即单元有形变），但插值计算得到的应变却为零（譬如一个正方形单元变形为一个等腰梯形，节点位移相等但符号相反，各形函数相同，所以插值结果为0），这样内能计算出来为零（单元没变形！）。显然，麻烦来了。所以必须避免这种情况的出现。判别出现0能模式的方法最简单的是察看单元变形情况，如果单元变成交替出现的梯形形状（两的在一起有点像沙漏，windows里面那个动画光标），就得小心了。另外就是察看hourglass energy，最好不要超过总能量的5%.

18)我以前做过一个类似轧钢的问题，使用的单位制为国际单位制，当采用实际的密度（7800kg/m3）时，网格畸变得特别厉害，和你所说的沙漏变形有点象。可是无意间将物质的密度修改大（7.8e7）以后，网格就没有畸变那么厉害了，并且计算速度大幅度提高。不知这里面有什么说道，请明白人给指点一下，谢谢！

19)畸变：应该可以用自适应网格划分，要注意的是：*CONTROL_ADAPTIVE 和*PART中的相应的设置

12）对PART进行沙漏控制，使用use *hourglass card 和*PART card。感觉与*CONTROL_HOURGLASS card的功能相同，前者控制具体part，后者控制整个模型

21）书上说，通过使用好的模型方式可以减少沙漏的产生，如网格的细化、避免施加单点载荷、在易产生沙漏模式的部件中分散一些全积分的?#31181;子?#21333;元，从而减少沙漏。我的理解是在那些地方放一些全积分的单元。其他的都用单点积分单元 。我的理解是增加积分点

22）164单元不能自动分裂单元。
大变形时确实存在网格畸形的问题，一般可以如下解决：
1、当畸形单元少，位置不重要时，删除畸形单元，继续计算；
2、适当调整参数，减少畸形的出现，或者使畸形单元破坏（ADD_EROSION）。
3、使用ALE。
版主所说的第一种解决办法，是通过重启动来实现的么？ 、使用ALE。

23)今天突然注意到dyna关键字中的*section_shell_ale中的第八项setype中有3个选项，1-lagrange;2-euler,3-ale，而在*control_ale中的第一项dct有4项，1-lagrange;2-euler,3-ale，4-ale ambient。
请高手指教：这两个关键字中的选项是不是要对应，那一个优先级更高？还有，如果都选1，就是lagrange单元，那还定义ale干什么？我一直困惑，希望给出详细的解释，谢谢。
根据经验判断，*section_shell_ale的定义应该更高一些，因为在dyna中，普通的关键字总是小于特定part的关键字定义的。

在dyna中，默认Lagrange单元、结构动力学分析、显示。因此在设置ALE（*control_ale）、结构热力学分析（*CONTROL_SOLUTION）、隐式（*CONTROL_IMPLICIT_GENERAL）的时候，都会发现系统默认的都是前者的现象，这是DYNA的风格，也是为了避免出现意外。

24)完全重启动是一个全新的分析，只不过需要考虑前面分析的相关的part的变形和应力情况，这个传递需要由关键字*STRESS_INITIALIZATION_OPTION来实现。
在dos命令下输入：
LS-DYNA I=restartinput.k R=D3DUMPnn
与小型重启动不同的是restartinput.k文件种具有完整的关键字输入，节点、单元、增加的PART等，对于需要初始化的part，在该输入文件种的节点、单元数及排列与拓扑关系都应该与上次求解的输入文件中的一样，但part号可以不同。《动力显示算法指南》

25)我做了一个爆炸分析的流固耦合，但是流体和固体没作用上？
流固耦合关键字里面有一个mcoup选项，你选的是1，表示只与密度最大的多物质材料耦合，这样的话就只和炸药耦合了，空气对他的作用就没有了

26)我刚才看了一下dyna的范例手册，和您说的一样，也是想用part来定义，简单方便。我是用hypermesh来作为前处理的，很方便。具体和大家说一下，在hypermesh中
BCs->entity sets-〉选择comps，然后输入要定义的set的名称，选择要包括的component，最后再选择create就行了。
hypermesh作为dyna的前处理非常好，推荐使用！

非常不错