下载文档原格式
Abaqus有限元分析中的沙漏效应[转]2011-09-2117:34:27|分类：有限元|字号大中小订阅1.沙漏的定义沙漏hourglassing一般出现在采用缩减积分单元的情况下：比如一阶四边形缩减积分单元，该单元有四个节点“o”，但只有一个积分点“*”。
而且该积分点位于单元中心位置，此时如果单元受弯或者受剪，则必然会发生变形，如下图a所示。
关于沙漏问题，建议看看abaqus的帮助文档，感觉讲的非常好，由浅入深，把深奥的东西讲的很容易理解。
沙漏的产生是一种数值问题，单元自身存在的一种数值问题，举个例子，对于单积分点线性单元，单元受力变形没有产生应变能--也叫0能量模式，在这种情况下，单元没有刚度，所以不能抵抗变形，不合理，所以必须避免这种情况的出现，需要加以控制，既然没有刚度，就要施加虚拟的刚度以限制沙漏模式的扩展---人为加的沙漏刚度就是这么来的。
关于沙漏现象的判别，也就是出现0能模式的方法最简单的是察看单元变形情况，就像刚才所说的单点积分单元，如果单元变成交替出现的梯形形状，如果多个这样的单元叠加起来，是不是象我们windows中的沙漏图标呢？2.ABAQUS中沙漏的控制：*SECTION CONTROLS：指定截面控制警告：对于沙漏控制，使用大于默认值会产生额外的刚度响应，甚至当值太大时有时导致不稳定。
默认沙漏控制参数下出现沙漏问题表明网格太粗糙，因此，更好的解决办法是细化网格而不是施加更大的沙漏控制。
该选项用来为减缩积分单元选择非默认的沙漏控制方法，和standard中的修正的四面体或三角形单元或缩放沙漏控制的默认系数；在explicit中，也为8节点块体单元选择非默认的运动方程：为实体和壳选择二阶方程、为实体单元激活扭曲控制、缩放线性和二次体积粘度、设置当单元破损时是否删除他们、或为上述完全破损的单元指定一标量退化参数。
等必需参数：NAME：名字可选参数：DISTORTION CONTROL：只用于explicit分析。
Abaqus有限元分析中的沙漏效应[转]2011-09-21 17:34:27| 分类：有限元 | 标签： |字号大中小订阅1. 沙漏的定义沙漏hourglassing一般出现在采用缩减积分单元的情况下：比如一阶四边形缩减积分单元，该单元有四个节点“o”，但只有一个积分点“*”。
而且该积分点位于单元中心位置，此时如果单元受弯或者受剪，则必然会发生变形，如下图a所示。
关于沙漏问题，建议看看abaqus的帮助文档，感觉讲的非常好，由浅入深，把深奥的东西讲的很容易理解。
沙漏的产生是一种数值问题，单元自身存在的一种数值问题，举个例子，对于单积分点线性单元，单元受力变形没有产生应变能--也叫0能量模式，在这种情况下，单元没有刚度，所以不能抵抗变形，不合理，所以必须避免这种情况的出现，需要加以控制，既然没有刚度，就要施加虚拟的刚度以限制沙漏模式的扩展---人为加的沙漏刚度就是这么来的。
关于沙漏现象的判别，也就是出现0能模式的方法最简单的是察看单元变形情况，就像刚才所说的单点积分单元，如果单元变成交替出现的梯形形状，如果多个这样的单元叠加起来，是不是象我们windows中的沙漏图标呢？ABAQUS中沙漏的控制：*SECTION CONTROLS：指定截面控制警告：对于沙漏控制，使用大于默认值会产生额外的刚度响应，甚至当值太大时有时导致不稳定。
默认沙漏控制参数下出现沙漏问题表明网格太粗糙，因此，更好的解决办法是细化网格而不是施加更大的沙漏控制。
该选项用来为减缩积分单元选择非默认的沙漏控制方法，和standard中的修正的四面体或三角形单元或缩放沙漏控制的默认系数；在explicit中，也为8节点块体单元选择非默认的运动方程：为实体和壳选择二阶方程、为实体单元激活扭曲控制、缩放线性和二次体积粘度、设置当单元破损时是否删除他们、或为上述完全破损的单元指定一标量退化参数。
等必需参数：NAME：名字可选参数：DISTORTION CONTROL：只用于explicit分析。
Ls_dyna 沙漏控制考试试题一、选择题1、沙漏能影响和以及单元，但不能影响三角形壳单元和梁单元。
A、实体单元B、四边形单元C、一维单元D、二维单元2、一般情况下，如果沙漏能量超过总能量的，那么就需要调整沙漏控制，以保证计算结果的精度。
A、10%B、15%C、20%D、25%3、沙漏模式是一种非物理的变形模式，产生零应变和应力。
A、超能B、零能C、刚性D、柔性4、如果系统计算沙漏能并将其输出，必须设置关键字*CONTROL_ENERGY中的控制参数HGEN = 。
A、2B、3C、4D、55、在*SECTION_SHELL中，没有沙漏模式的壳单元公式和。
A、1B、2C、6D、76、如果不考虑计算成本，使用单元的方法，可以有效地控制沙漏现象。
A、全积分B、单点积分C、缩减积分D、选择性缩减积分7、在关键字*CONTROL_HOURGLASS全局控制沙漏中，QH一般取值范围之间，如果超过这个范围将引起计算不稳定。
A、0—0.05B、0.05—0.15C、0.15—0.20D、0.20—0.258、如果需要针对每个PART单独控制，需要定义关键字。
A、*PARTB、*CONTROL_HOURGLASSC、*CONTROL_ENERGYD、*HOURGLASS9、在关键字*HOURGLASS中设定沙漏系数是项。
A、IBQB、QMC、Q1D、Q210、在Ls_dyna关键字*CONTROL_HOURGLASS全局控制沙漏中，QH默认值是。
A、0.1B、1C、0.01D、0.001选择题参考答案：1、ABD2、A3、B4、A5、CD6、AD7、B8、D9、B 10、A二、有限单元单点积分方法可以大幅度降低计算成本，但容易沙漏模式，产生一种自然振荡并且比所有结构响应的同期短得多。
沙漏变形没有刚度并产生锯齿形外形，由此使计算结果不正确。
请回答有哪些方法可以有效地控制沙漏现象？参考答案：控制沙漏现象的方法：1、使用均匀网格，避免在单点上集中加载；2、调整模式的体积粘性；3、使用单元的全积分方法；4、增加模型的弹性刚度；5、局部增加模型刚度；。
沙漏(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的论文中有更详细的描述。
