一、【子程序】v u mat 有沙漏问题么? 沙漏问题和 VUMAT 无关,跟你选择的单元有关系,如果你采用减缩积分单元,则会存在沙漏。 有限元的一个核心就是单元模型,其思想是采用单元近似连续体,单元内采用形函数进行插值。采用全积分的话,可以精确地积出刚度矩阵,但是采用全积分会导致有限元过刚,例如体积锁死和剪切锁死等,因此很多力学及提出了各种各样的单元模型来解决这些问题。现在用的较多的低阶单元就是一点积分,一点积分的单元由于积分点过少而存在零能模式(沙漏),即在某些变形模式下会出现零应变,这个可以从形函数的公式中推导出来。 所以,沙漏模式是否存在取决你选用的单元,但是你采用 ABAQUS 的默认设置基本上就可以解决这个问题。 不知道我有没有说清楚 二、【基础理论】【概念】剪切锁死、体积锁死、沙漏、零能模式 1.剪切锁死(shear locking) 简单地说就是在理论上没有剪切变形的单元中发生了剪切变形。该剪切变形也常称伴生剪切( parasitic shear)。 发生的条件:1.一阶、全积分单元;2.受纯弯状态; 产生的结果:使得弯曲变形偏小,即弯曲刚度太刚。 解决方法:1.采用减缩积分;2.细化网格;3.非协调单元;4.假定剪切应变法; 2.体积锁死(v olu metric locking) 简单地说就是应该有单元的体积变化的时候体积却没发生变化。该原因是受到了伪围压应力(Spu riou s pressu re stresses )。 发生的条件:1.全积分单元;2.材性几乎不可压缩; 二阶单元:对于弹塑性材料(塑性部分几乎属于不可压缩),二阶全积分四边形和六面体单元在塑性应变和弹性应变在一个数量级时会发生体积锁死。二次减缩积分单元发生大应变时体积锁死也伴随出现。 但值得注意的是,一阶全积分单元当采用选择性减缩积分(selectiv ely redu ced integration)时可以避免出现体积锁死。 产生的结果:使得体积不变,即体积模量太大,刚度太刚。 解决方法:1.将大应变区域网格细化;2.mix ed formu lation 法; 检查方法:输出积分点的围压应力,分析围压应力是否在相邻积分点存在突变,是否显棋格式分布,是的话就说明出现体积锁死。 3.沙漏(hourglassing) 简单地说就是单元只有一个积分点,周边的节点可以随意变形。 发生的对象:1.一阶、减缩积分单元; 产生的结果:单元太柔; 解决方法:1.对一阶减缩单元,合理细化网格;荷载避免使用点荷载; 2.在大应变区或大应变梯度区使用一阶单元,...