Lsdyna控制卡片介绍（2）

沙漏控制*CONTROL_HOURGLASS

为节约计算时间，提高计算效率，大多数积分计算均采用部分积分类型进行积分计算。采用部分积分类型有可能出现沙漏，当沙漏较大时会导致模型能量不守恒，影响计算精度甚至导致计算结果不可信，所以必须对计算沙漏进行相关控制。其中沙漏控制分为对单个件的局部沙漏控制和对整体的总体沙漏控制两种。

IHQ：总体附加刚度或黏性阻尼方式选项；  

  1：标准LS-DYNA类型。（默认）  

 2：Flanagan-Belyschko积分类型。  

 3：有精确体积的Flanagan-Belyschko积分类型。  

 4：类型2的刚度形式。  

 5：类型3的刚度形式。  

 6：········  

8：适用于单元类型为16的全积分壳单元。当IHQ=8时，激活翘曲刚度，以得到精确解。该选项会增加25%的计算开销。  

        在壳单元中，IHQ<4的是基于Belyschko-Tsay公式的粘性沙漏控制模式，IHQ=4，5，6为刚度控制模式。刚度控制模式在大变形问题中可能使响应变得过于刚硬，使用时要注意。在高速问题中推荐采用粘性模式，在低速问题中推荐采用刚度模式。对于大变形问题，推荐使用选项3或5。  

QH：沙漏能系数 ,超过0.15会导致计算不稳定。可适用于除IHQ=6以外的所有选项。  

  备注：对个别组件的沙漏控制，可通过先建立沙漏属性集 合器，再从组件集 合器中调用沙漏属性的方法实现。  

单元控制 CONTROL_SHELL

壳单元控制方法设置。

WRPANG：壳单元翘曲角度。当某个翘曲角度大于给定值时，会输出警告信息。默认值为20；  

ESORT：自动挑选退化的四边形单元，并处理为CO三角形单元公式，以保证求解稳定。  

  0：不挑选。（默认）  

  1：完全挑选并处理。  

IRNXX：单元法线更新选项。该选项影响Hughes-Liu，  

Belytschko-Wong-Chiang，和Belytschko-Tsay单元公式。当且仅当翘曲刚度选项被激活时，即BWC=1时，以上单元公式才受影响。对于Hughes-Liu壳单元类型1，6和7，IRNXX必须设为-2以调用上表面或下表面作为参考面。  

   -2：··········  

  -1：每个循环都重新计算法线方向。  

  0：自动设为-1。  

  1：重启动时计算。  

n：每n个循环重新计算法线方向。（只适用于Hughes-Liu壳单元类型）  

ISTUPD：单元厚度改变选项。该选项对所有壳单元变形有影响。  

  0：不变化。  

  1：膜变形引起厚度改变。该选项对金属板料成型或所有膜片拉伸作用很大的情况都很重要。  

THEORY：壳单元使用的理论。（默认的是Belytschko-Tsay，面内单点积分，计算速度很快，采用Co-rotaional应力更新，单元坐标系统置于单元中心，基于平面单元假定，建议在大多数分析中使用）。  

BWC：针对Belytschko-Tsay单元的翘曲刚度。  

  1：增加Belytschko-Wong-Chiang公式的翘曲刚度。  

  2：Belytschko-Tsay单元公式。不增加翘曲刚度。（默认）  

MITER：平面应力塑性选项，默认为1。（运用于材料3，18，19和24）。  

  1：3次交叉迭代。（默认）  

  2：完全迭代。  

  3：不迭代。可能导致错误，慎用。  

PROJ：在Belytschko-Tsay和Belytschko-Wong-Chiang单元中翘曲刚度投影方法。这个方法主要运用于显示分析，如果是隐式分析，那此项无效 。默认为0。  

OTASCL：为旋转单元质量定义一个缩放系数。（不太常用）。  

INTGRD：通过厚度数值积分法则的默认壳单元。当积分点为1到2个的时候使用Gauss积分，当积分点从3个到10的时候使用Lobatto积分，积分点为2个时，Lobatto法则非常不准，须用Gauss积分。  

LAMSHT：薄壳理论开关。0：不更新切应变修正；1：薄壳理论切应变修正。  

CSTYP6：第6种壳单元坐标系的选用。1：可变的局部坐标系（默认）；2：统一局部坐标系（计算结果有偏差，但效率比较高）。  

TSHELL：允许热传导通过有厚度的壳单元。