爱华网比如一阶四边形缩减积分单元，该单元有四个节点“0”，但只有一个积分点“*”。而且该积分点位于单元中心位置，此时如果单元受弯或者受剪，则必然会发生变形.但是，现实的情况却是在这三种情况下，单元积分点上的主应力和剪应力状况都没有发生变化，也就是说该单元可以自由地在这三种形态之间转变而无需外很小的扰动理论上可以让单元无限地变形下去，而不会消耗任何能量，这就是所谓的沙漏的零能量模式。这时就要对沙漏进行控制，比如人为地给单元加上一定的刚度。负体积是由于element本身产生大变形造成自我体积的内面跑到外面接着被判断为负体积沙漏选项关键字看完因为没看懂，就试着翻一下，翻成汉语感觉还是不懂。
*HOURGLASS
目的：定义沙漏和体积粘性特性，他们可以在*PART命令汇总的通过HGID参数引用。这里指定的特性通过特定的PART引用后会忽略在*CONTROL_HOUGLASS和*CONTROL_BULK_VISCOSITY定义的相关参数。
沙漏ID号。
沙漏控制类型。对于SOLID单元有6的选项。对于四边形的SHELL和MEMBRANE单元，沙漏控制是基于Belytschko和Tsay公式等，比如，1-3选项是一样的，4-6选项是一样的。*a& h! x0 [) e' `7 S/ I6E
EQ.0:默认1不论在*control_hourglass卡片中IHQ是什么值- m. F7 i) @! t)x
EQ.1：标准LS-DYNA粘性形式
EQ.2：Flanagan-Belytschko粘性形式- J( `& U& T0 d( ]1v& @( ^, E5Z
EQ.3：Flanagan-Belytschko粘性形式并具有精确的体积积分对于Solid单元来说
EQ.4:Flanagan-Belytschko刚度形式
EQ.5:Flanagan-Belytschko刚度形式，并具有精确体积积分对于Solid单元来说! G: F2 Y& [" o, S+J& |( {*S
EQ.6:Belytschko-Bindeman假定的应变共旋转刚度形式仅针对2D和3DSolid单元。这个形式隐式和显示均可用。实际上，类型6对于隐式是必须的选项。.F" 4 z% y: [5 V/x
EQ.8:对于类型16全积分Shell单元适用，IHG=8季候的精确求解用的翘曲刚度。对于这个选项，速度罚系数25%是通用的。)f7 ^2 ]& W}$ i8h% ~# F4 b+ f
QM:沙漏系数。超过0.15的QM值会导致结果不稳定。推荐所有选项采用默认值除了IHQ=6。然而刚度形式能刚化响应尤其变形很大的时候，因为要小心。对于Shell和Membrane单元QM被当成Membrane沙漏系数，QB是弯曲系数，QW是翘曲系数。可以分别指定这些系数，但一般来说，QM=QB=QW就足够了。对于类型6Solid单元沙漏控制，看下面的4注解。- Y: s; l, T% ^' I!k
IBQ:体积粘性系数- Y$ p) |- Z2 ?! a9 y- l! Y# X& P"P1 ?
Q1：二次体积粘性系数) o5 e# C4 O0 F# b"L
Q2：现行体积粘性系数，对于壳弯曲QB是沙漏系数。默认QB=QM8a0 `' V4 q3 m! `8l
QW：对于壳弯曲是沙漏系数，默认QB=QW。
评论：
1.粘性沙漏控制推荐用于高速变形问题。刚度控制通常用于低速问题，尤其当时间步比较大的时候。对于高度扭曲的固体单元精确的体积积分具有优势。:p, CT! O9 o5 v+a
2.对于汽车碰撞问题，沙漏控制系数是0.05的刚度形式是常用的。; q; r' A,Lw5 ?; x%n
3.在固体材料中体积粘性是必要的对于冲击波的传播问题，因而仅用于固体单元。一般来说，除了压力比较大是采用更大的值之外，采用默认值就可以了。在底密度泡沫中由于粘性应力较大，有必要降低粘性值。当然降低一个数量级以上是不理想的。$r. h9 B" e6 f, X1 ~(E
4.类型6的沙漏控制是仅对于2D和3D固体单元的。基于弹性常数和假定的应变场，能产生精确的粗糙网格弯曲结果对于QM=1.0的弹性材料。对于具有屈服应力切线模量（它比弹性模量小的多）的塑性模型，较小的QM值（0.001-0.1）能产生更好的结果。对于任何材料请记住刚度是基于弹性常数的。所以如果材料软化，QM值小于1会更好。对于各向异性材料，采用平均弹性常数。对于null材料定于的流体模型，类型6沙漏控制是粘性的并与材料粘性系数成比例。4n: ; ]U: ~9|
5.有时候，Belytschko-Lin-Tsay和下积分Hughes-Liu壳单元的计算效率是从平面内单点积分导出的。为了降低伴随单点积分出现的沙漏变形，基于应力的沙漏粘性和刚度附加在物理刚度上在局部单元水平范围内。接下来的沙漏控制与所有的单点积分Shell和Membrane单元有关。
