沙漏模型及其控制.docVIP

沙漏模型及其控制 沙漏模型及其控制2011-01-04 18：58A1：有限元方法一般以节点的位移作为基本变量，单元内各点的位移以及应变均采用形函数对各节点的位移进行插值计算而得，应力根据本构方程由应变计算得到，然后就可以计算单元的内能了。如果采用单点积分(积分点在等参元中心)，在某些情况下节点位移不为零(即单元有形变)，但插值计算得到的应变却为零(譬如一个正方形单元变形为一个等腰梯形，节点位移相等但符号相反，各形函数相同，所以插值结果为0)，这样内能计算出来为零(单元没变形！)。这种情况下，一对单元叠在一起有点像沙漏，所以这种模式称之为沙漏模式或沙漏。 现在有很多控制沙漏的专门程序，如控制基于单元边界的相对转动。但这些方法不能保持完备性。： 我主要讲一下物理的稳定性，在假设应变方法的基础上，建立沙漏稳定性的过程。在这些过程中，稳定性参数基于材料的性能。这类稳定性也称为物理沙漏控制。对于不可压缩材料，即使当稳定性参数是一阶的时候，这些稳定性方法也将没有自锁。在建立物理沙漏控制中，必须做出两个假设：1.在单元内旋转是常数。2.在单元内材料响应是均匀的。 A2：沙漏(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 1315)类型6沙漏控制调用了一种假设应变协同转动方程。使用沙漏控制类型6和系数1.0，一个弹性部件在厚度方向仅仅需要划分一层类型1的体单元就可以获得正确的弯曲刚度。在隐式计算里面，对于类型1的体单元应该总是使用类型6的沙漏控制(实际上，在V970里面这是自动设置的)。 (More on type 6HG control from Lee Bindeman) 类型6的沙漏控制与类型4,5不在于它用了一个假设应变场和材料属性来估算出假设应力场。这个应力在单元封闭域内进行积分得到沙漏力，因此单元表现的像一个有同样假设应变场的全积分单元。这种假设应变场设计成用来阻止纯弯曲中不真实的剪切变形和近似不可压材料中的体积锁死。 类型4和5的沙漏控制基于单元体积，波速和密度像在LS-DYNA理论手册中方程3.21那样来计算沙漏刚度。 沙漏类型6主要的改进是应力场在单元域内积分。这使得当使用大的长细比或者歪斜形状的体单元时沙漏控制非常鲁棒。类型4和5的沙漏控制对大长细比和歪斜形状单元反应变不好，它趋向于对某些沙漏模式反应的过于刚硬而对其它模式反应得过弱。 沙漏控制类型6另