沙漏控制类型有哪几种

zhan秋波

沙漏控制类型（Hourglass control type）有哪几种，每种类型具体表示什么含义？

zhan秋波

Ls-dyna使用手册上有详细信息

lk18996

zhan秋波 发表于 2012-8-5 19:06
Ls-dyna使用手册上有详细信息

请给出完整的回答！

zhan秋波

lk18996 发表于 2012-8-7 09:01
请给出完整的回答！

沙漏(hourglass)模式是一种非物理的零能变形模式，产生零应变和应力。沙漏模式仅发生在减缩积分(单
积分点)体、壳和厚壳单元上。LS-DYNA 里面有多种算法用于抑制沙漏模式。缺省的算法(type 1)通常不是最
有效的算法，但却是最经济的。
一种完全消除沙漏的方法是转换到全积分或者选择减缩积分(S/R)方程的单元。但这种方法是一种下策。
例如，第一，类型2 体单元比缺省的单点积分体单元计算开消大; 其二，在大变形应用时更不稳定(更容易出
现负体积);其三，类型2 体单元当单元形状比较差时在一些应用中会趋向于剪切锁死(shear-lock)，因而表现得
过于刚硬。
三角形壳和四面体单元没有沙漏模式，但缺点是在许多应用中被认为过于刚硬。
减小沙漏的一个好的方法是细化网格，但这当然并不总是现实的。
加载方式会影响沙漏程度。施加压力载荷优于在单点上加载，因为后者更容易激起沙漏模式。
为了评估沙漏能，在*control_energy 卡片中设置HGEN＝2，而且用*database_glstat 和*database_matsum
卡分别输出系统和每一个部件的沙漏能。这一点是要确认非物理的沙漏能相对于每一个part 的峰值内能要小
(经验上来说<10%)。对于壳单元，可以绘制出沙漏能密度云图，但事先在*database_extent_binary 卡中设置
SHGE=2。然后在LS-Prepost 中选择Fcomp>Misc>hourglass energy。
对于流体部件，缺省的沙漏系数通常是不合适的(太高)。因此对于流体，沙漏系数通常要缩小一到两个数
量级。对流体用基于粘性的沙漏控制。缺省的沙漏方程(type 1)对流体通常是可以的。
对于结构部件一般来说基于刚性的沙漏控制(type 4,5)比粘性沙漏控制更有效。通常，当使用刚性沙漏控
制时，习惯于减小沙漏系数到0.03~0.05 的范围，这样最小化非物理的硬化响应同时又有效抑制沙漏模式。对
于高速冲击，即使对于固体结构部件，推荐采用基于粘性的沙漏控制(type 1,2,3)。
粘性沙漏控制仅仅是抑制沙漏模式的进一步发展，刚性沙漏控制将使单元朝未变形的方向变形。
类型8 沙漏控制仅用于单元类型16 的壳。这种沙漏类型激活了16 号壳的翘曲刚度，因此单元的翘曲不
会使解退化。如果使用沙漏控制8，16 号壳单元可以用于解被称为扭曲梁(Twisted Beam)问题。
对于单元类型1 的体和减缩积分2D 体(shell types 13 & 15)类型6 沙漏控制调用了一种假设应变协同转动
方程。使用沙漏控制类型6 和系数1.0，一个弹性部件在厚度方向仅仅需要划分一层类型1 的体单元就可以获
得正确的弯曲刚度。在隐式计算里面，对于类型1 的体单元应该总是使用类型6 的沙漏控制(实际上，在V970
里面这是自动设置的)。
(More on type 6 HG control from Lee Bindeman)
类型6 的沙漏控制与类型4,5 不在于它用了一个假设应变场和材料属性来估算出假设应力场。这个应力在单
元封闭域内进行积分得到沙漏力，因此单元表现的像一个有同样假设应变场的全积分单元。这种假设应变场
设计成用来阻止纯弯曲中不真实的剪切变形和近似不可压材料中的体积锁死。
类型4 和5 的沙漏控制基于单元体积，波速和密度像在LS-DYNA 理论手册中方程3.21 那样来计算沙漏
刚度。
沙漏类型6 主要的改进是应力场在单元域内积分。这使得当使用大的长细比或者歪斜形状的体单元时沙
漏控制非常鲁棒。类型4 和5 的沙漏控制对大长细比和歪斜形状单元反应变不好，它趋向于对某些沙漏模式
反应的过于刚硬而对其它模式反应得过弱。
沙漏控制类型6 另一个理论上的优点是对在厚度方向只有一个单元的梁可以在弹性弯曲问题中得到准确
的解。要做到这一点，设置沙漏刚度参数为1.0。同样，对弹性材料方形截面杆的扭曲问题，当沙漏系数设为
1.0 时可以用很少的单元来解。然而，对于非线性材料，用粗糙的网格得到好的结果是不可能的，因为应力场
不是像沙漏类型6 假设的那样线性变化的。在梁厚度方向上如果没有更多积分点的话，没有办法捕获应力场
的非线性状态。
对于选择沙漏控制，下面几个问题要考虑。对于单元有大的长细比或者明显歪斜(不管是初始还是变形过
程中)，推荐采用类型6 的沙漏控制。类型6 的沙漏控制通常对软的材料更好，像泡沫或蜂窝材料在计算中会
有非常明显的变形。
在材料不是特别软或者单元有合理的形状且网格不是太粗糙时，类型4,5 和6 沙漏控制似乎都能得到同
样的结果。这种情况推荐用类型4 的沙漏控制，因为它比其它的更快。
类型6 的沙漏控制在LS-DYNA User’s Manual 中参考的Belytschko 和Bindeman 的论文中有更详细的描述。

graf

solid hgid=6
shell elform=2 hgid=4 or 5
shell elform=16 hgid=8