A-A+

Hyperworks Radioss中的2D壳单元SHELL

2019年11月21日 CAE 阅读 1,122 views 次

Hyperworks RADIOSS里标准的Shell单元是带横向剪切变形的Mindlin板单元,这种单元比标准的Kirchhoff单元精确,后者不考虑横向剪切变形,以至于初始垂直于中面的几何横截面在变形后依然垂直于中面,会产生剪切锁定。关于剪切锁定的介绍可参考文章《有限元分析中的剪切锁定(shear lock)》《消除有限元分析中剪切锁定(shear lock)的方法》

尽管Kirchhoff模型精度稍低,如果L/h>20, Kirchhoff假设还是正确的。但是如果L/h在10~20时,这种初始垂直于中面的横截面依然垂直于中面的假设就不再成立。而Mindlin板始终考虑横向剪切,不存在理论假设不成立的问题。RADIOSS Block里减缩积分shell单元(无论3节点还是四节点)都是基于Mindlin-Ressner理论,没有特别的公式用于解决单元节点偏离几何中面的情况,因此对于薄壁结构的网格离散,使用中面划分网格是非常重要的。这种公式不但对薄板壳有效,对于厚板壳也基本适用。

RADIOSS不但有因为计算效率而最流行的减缩积分Shell单元,还有完全积分的四边形壳单元(QBAT)和三角形单元(DKT18)。在实际分析中对单元和公式的选择应注意如下问题:

(1)BT单元简单、高效、计算开销小。作为减缩积分单元,这种单元对网格质量不是太敏感,可以用于较粗糙的网格。

(2)对于准静态分析、精细网格、翘曲表面、屈曲,这些情况选择QEPH或者QBAT是明智的。

(3)QBAT是RADIOSS里精度最高的Shell单元公式,但是由于它是完全积分,其计算开销是BT单元的2~3倍。

(4)QEPH是精度和计算开销协调最好的Shell单元,它相对BT单元增加的计算开销不大于15%,而其得到的结果接近于完全积分的QBAT单元。

(5)不建议使用三角形单元。C0单元太刚硬,而DK,18计算消耗很高。为了保证结果的质量,一般应限制三角形单元比例不超过5%.

(6)S3N6单元在弯曲行为精度同DKT18,可以应用在一些特殊的领域,比如冲压仿真。

对于弹性Shell,正应力沿厚度方向是线性分布的,因此单元内能可以使用解析积分方法获得。在塑性情况下,这个应力分布是非线性的,并且至少需要在厚度方向3个积分点才可以将这种非线性考虑进去。这种非线性的应力分布需要增加积分点数才能够计量,同时单元内能的估算精度依赖于积分点数。因此需要在考虑材料非线性、料厚和弯曲速率的前提下调和计算开销与解的精度关系。厚度方向使用5个积分点能够给出较好的结果,尤其对于厚度超过2m的板,当然这样带来的CPU开销的增加也是不能忽略的。仿真无弯图和横剪的膜单元时,厚度方向使用1个积分点就够了。当然,对于线弹性材料(Law1),由于是可解析的,厚度方向积分点数可以忽略。

另一个低CPU开销而结果相当精确的做法是使用全局积分,这种积分方式把Mises望性准则转化为成Illuyshin准则,这样积分点处的应力分量被内力(N,M,T…)取代。

个人公众号“数字化设计CAX联盟”,欢迎关注,共同交流
为您推荐:

给我留言

© 坐倚北风 版权所有 严禁镜像复制 苏ICP备15034888号. 基于 Ality 主题定制 AliCMS
联系邮箱:leanwind@163.con,微信公众号:数字化设计CAX联盟

用户登录

分享到: