第6章LS-DYNA壳单元、沙漏剖析.ppt

沙漏类型和参数 粘性公式 减慢和停止HG 变形(不能恢复) ?标准的 LS-DYNA (缺省) ? Flanagan-Belytschko (2) ? 具有准确体积积分(3) 的Flanagan-Belytschko 刚度公式 减慢，停止和恢复HG 变形 ? Flanagan-Belytschko (4) ?具有准确体积积分(5) 的Flanagan-Belytschko HG 系数，体积粘度(对实体)，壳弯曲和翘曲系数 沙漏 – 注释 粘性 HG用于高速变形问题 刚度控制更适用于低速 对于实体单元，确定的积分对于一些高度扭曲的单元有用 Flanagan-Belytschko 公式对大旋转问题效果好一些 HG 常发生于一些小位移情况，特别是动力松弛时 ?试试刚度方式 ?这种方式可能对一些大变形问题会有过刚度 对于汽车碰撞，刚度公式的系数为0.05为大多数用户使用 沙漏-关键字 *CONTROL_ENERGY 控制是否计算沙漏能 (较少的 CPU开销) *CONTROL_HOURGLASS . 设置全局沙漏类型 (缺省是标准的 LS-DYNA粘性公式) . 能修改沙漏系数 *HOURGLASS . 对指定的part设置沙漏类型和参数 *PART .对指定的part改变全局的沙漏类型和参数,指定一个沙漏ID号 沙漏 如果沙漏产生的区域并不影响关心的设计区域，则可接受 全积分单元没有沙漏 HG 模式对于实际的变形是直角正交的 由沙漏控制所做的功不会出现在能量方程中 当沙漏出现时，总能量有微小的降低 沙漏能量耗散记录在 GLSTAT 和 MATSUM中 一点积分是非常快的，所以我们接受这种冒险 但应不时检查能量平衡是否安全 一般规则: 沙漏能 10% 内能 总体 (GLSTAT) 单个parts (MATSUM) . 做可视检查 . 观察过度的刚度行为 沙漏的预防 通常建议不要调整HG系数，而应该查找问题的根源 沙漏模式常可避免，通过使用好的模型方式： .网格细化 .统一的网格 . 分布式载荷 如果改进网格和载荷也不能奏效，使用全积分公式 沙漏例子—体单元 沙漏例子—体单元 第六章 单元 II 壳单元 沙漏 壳单元 用于某一方向相对较薄的物体 ? 钣金、薄壁结构、发动机叶片、罐头等 ?防撞性、乘客安全仿真，钣金成型，飞机撞击，导弹上的冲击载荷 *ELEMENT_SHELL ? 单元ID号 ? part ID号 ? 四面体的4个节点三角形的 3个节点 ? 取代每一个节点上的缺省厚度 管与管的冲击 管长= 50 in O.D. = 6.625 in 厚度 = 0.432 钢Et = 105, ┮y = 105 壳单元的特性 有限应变 厚度方向积分点的任意和固定 壳单元厚度更新 单元卡上几何特性可任意指定 全矢量或平行的 Fully vectorized and parallelized. 所有壳单元共享的构造子程序 可使用共同的局部坐标系 沙漏控制来控制零能模式 在塑性屈曲问题中单元公式的不同常引起显著的差异 在DYNA3D中第一种实行的壳单元为Hughes-Liu 壳，它能处理有限应变。在LLNL的DYNA3D 1986版本中出现 。1983年就出现了壳单元，但计算非常慢，如考虑管对管冲击问题： . 3 solid elements per shell solution time = 12 CPU minutes . 发布时的Hughes-Liu 壳 = 600 CPU minutes . 1986年 改进的壳 = 4 CPU minutes (XMP) . 在 RS6000 工作站上的快的Hughes-Liu 壳 4 CPU minutes 加速这种壳单元需要相当的努力 壳单元技术的历史 第一次发布的Hughes-Liu壳单元速度慢的几个原因 : . 2 x 2 选择简化积分点selective reduced integration . 无向量化 (慢8倍) . 中间步的应变计算 加速来自： . 1 点积分 . n + 1 几何形状用于应变计算 . 使用总体坐标系 增加粘性沙漏力控制零能模式 壳单元技术的历史 Belytschko-Tsay 壳 The B-T 壳单元由 Belytschko和Tsay 在 1981年发展起来的，并由 Belytschko, Lin 和Tsay在1984年改进。基于坐标旋转co-rotational 和速度-应变公式的组合 通过绑定坐标系于单元内，坐标旋转公式避免了非线性机械中的复杂性 相对速度.应变的共轭应力是 Cauchy 应力 壳运动假定节点是 c