Tet10 Element — 故障排除
更充实的内容请访问 tet10-element.html。
TET10的故障
TET10也会出现故障吗?比TET4好得多吧。
TET10很优秀,但问题并非零。需要了解已知问题。
中间节点位置不正
中间节点相关的问题具体是什么?
中间节点偏离CAD面,或相邻单元的中间节点不一致时会出现问题。
症状:
- 曲面上的应力呈锯齿形
- 单元的雅可比行列式为负(单元翻转)
- 位移局部异常
解决方案:
- 检查预处理器的"中间节点CAD贴附"设置
- 进行品质检查,确保没有雅可比行列式为负的单元
- 如果出现问题,在局部重新生成网格
雅可比行列式为负的单元能计算吗?
大多数求解器将雅可比行列式为负的单元视为错误而拒绝。Nastran会输出FATAL错误,Abaqus会显示"Negative eigenvalue"。只要有一个雅可比行列式为负的单元,分析就无法运行。
DOF数过多
TET10网格的DOF数太多,内存不足。
TET10的DOF是TET4的2.5倍(10个节点对4个节点),约是HEX8的3倍。对于大规模模型,DOF数成为问题。
解决方案:
- 利用对称条件 — 将DOF减少到1/2〜1/8
- 子模型法 — 整体使用粗网格,局部使用精密网格
- 迭代法求解器 — 相比直接法内存效率更高
- HEX与TET混合网格 — 将简单形状部分替换为HEX以减少DOF
- 自适应网格 — 从粗网格开始,仅在必要的部分进行细分
混合网格的连接怎么处理?
TET10和HEX20的连接需要在中间插入棱锥单元(PYRA13/14)。棱锥单元具有四边形面(HEX侧)和三角形面(TET侧),是过渡单元。主要预处理器可以自动生成。
接触面不稳定(C3D10的情况)
用C3D10进行接触分析时,应力呈棋盘状。
这是C3D10接触面稳定性问题。4点完全积分的C3D10在接触面上的压力会振荡。
解决方案:
- 切换到C3D10M — 改进的积分使接触面稳定
- 细化接触面网格 — 缓解振荡
- 改变接触算法 — 惩罚法→Lagrange乘数法
所以C3D10M是接触问题的标准做法。
在Abaqus中进行包含接触的TET10分析时,没有理由不使用C3D10M。使用C3D10然后接触变得不稳定是可预见的问题。
应力精度不如HEX
用HEX20和TET10解同一个问题,TET10的精度较低。
在相同DOF数的情况下,HEX20的精度确实更高。要使TET10达到HEX20的精度水平,需要2〜5倍的DOF。
但包括网格生成时间在内的总成本方面,TET10通常更有利。对于需要1天来生成HEX网格的几何形状,TET10自动网格可在几分钟内完成。
"HEX精度更高"和"TET10更实用"并不矛盾啊。
正确。精度效率(每单位DOF的精度)HEX更高。工作流效率(网格生成+分析+结果评估的总时间)TET10更高。选择哪种取决于项目的优先级。
总结
总结一下TET10故障的应对方法。
C3D10M(Abaqus)的接触稳定性和中间节点的CAD贴附是两个关键点。
掌握这两点,TET10就不会出现大问题。TET10是实践中最值得信赖的3D单元。
TET10中点偏离和负雅可比行列式
当TET10的边中点放在边长的25%以下或75%以上的位置时,会产生负的雅可比行列式。HyperMesh的Quality Check中抽取Jacobian<0.01的单元进行修正是实务中的标准做法。Nastran会输出Fatal 2186。在TET10自动网格生成中,曲面拟合时中点会贴附到曲面上,因此曲率激烈的部位最容易发生负雅可比行列式,用Local Remesh(HyperMesh)在局部重新生成网格是最高效的方法。
相关主题
有用
详情
错误