LS-DYNA_中的沙漏问题.pdf

节选段落一：
LS-DYNA 中的沙漏问题
沙漏能是由于在显示分析中采用缩减积分造成的，所谓缩减积分就是单元计算时积分点数少
于实际个数，这种操作能加快计算速度，但是会造成一种单元的零能模式，这就是沙漏。计
算要求沙漏能小于总能量的 5％时才认为结果是可靠的。有限元方法一般以节点的位移作为
基本变量，单元内各点的位移以及应变均采用形函数对各节点的位移进行插值计算而得，应
力根据本构方程由应变计算得到，然后就可以计算单元的内能了。


节选段落二：
如果采用单点积分（积分
点在等参元中心），在某些情况下节点位移不为零（即单元有形变），但插值计算得到的应
变却为零（譬如一个正方形单元变形为一个等腰梯形，节点位移相等但符号相反，各形函数
相同，所以插值结果为 0），这样内能计算出来为零（单元没变形！）。这种情况下，一对
单元叠在一起有点像沙漏，所以这种模式称之为沙漏模式或沙漏 现在有很多控制沙漏的专
门程序，如控制基于单元边界的相对转动。但这些方法不能保持完备性。 我主要讲一下物
理的稳定性，在假设应变方法的基础上，建立沙漏稳定性的过程。在这些过程中，稳定性参
数基于材料的性能。这类稳定性也称为物理沙漏控制。


节选段落三：
对于不可压缩材料，即使当稳定性参
数是一阶的时候，这些稳定性方法也将没有自锁。在建立物理沙漏控制中，必须做出两个假
设：1.在单元内旋转是常数。2.在单元内材料响应是均匀的。
简单的说，就是在数学计算上成立，而在物理上或实际中不存在的一种现象。
此处，主要讨论在 LSDYNA 里 CONTROL_HOURGLASS 中，Hourglass coefficient 的
调整。一般在模拟中，0.1 是默认值，在模拟过程中，发现在 0.1 至 0.2 之间都可以保持
稳定的低值。
模型如下：固定底边，在顶边节点上加载向上的力。