第6讲有限元建模专题二(更新版)(免费阅读).pptVIP

粗网格预测的孔底部位移是不准确的，但是，采用正常网格、细网格和很细的网格预测了基本相同的结果。因此，正常网格对所关注的位移而言是收敛的。 孔底部应力峰值的收敛比位移慢得多，这是因为应力和应变是由位移的梯度计算得到的。而要预测准确的位移梯度比计算准确的位移所需的网格更密。 网格的细分明显地改变了连接环在固定端处的应力计算值。随着网格的细分应力值继续增加。在连接环与母结构连接的角点处存在应力奇异性。从理论上讲，在这个区域的应力是无限大的，因此，在此处增加网格密度不会产生一个收敛的应力。 所有的结果值都与由粗网格预测的结果值进行了无量纲化对比。 建议 对复杂的三维模拟而言，实际上可利用的计算机资源常常限制了所采用的网格密度。在这种情况下，必须非常小心地应用从分析中得到的结果。粗糙的网格足以用来预测趋势。然而，应用由粗糙网格计算得到的位移和应力的具体量值应该很谨慎。 一般来说，没有必要对所分析的结构全部采用均匀的细划网格。应该在出现高应力梯度的地方采用细网格，而在低应力梯度或应力量值不被关注的地方采用粗网格。例如，细划的网格仅应用在高应力梯度的区域，而在其他地方则采用粗网格。例如： 从表中可见，结果与整体划分很细网格的结果相近。但是．这种局部细划网格的模拟与应用很细网格的分析相比大大节省了计算需要的CPU时间： 应力奇异是指在有限元模型中那些应力值无限大的点处。例如： 点荷载，如集中力或力矩作用处 孤立的约束点导致支反力如同点荷载。 尖角（零倒角半径）处 在应力奇异点处网格越细化，应力值也随之增加且不收敛 关于应力奇异 如果离感兴趣区域较远，可以在查看结果时通过不激活受影响的区域忽略它的影响 如果位于感兴趣区域，需要如下纠正: 在尖角处增加倒角重新进行分析 代替点力载荷为等效压力载荷 “散布” 位移约束至一个节点集 应力奇异的处理 6.4 应力集中问题 二次单元处理应力集中问题，明显优于一次单元 无论是完全积分还是减缩积分都可以很好的反映应力集中 减缩积分效率更高，而且计算结果往往优于完全积分。 二次单元可以以更少的单元更好的反应结构的几何特征 Physical model Model with first-order elements—element faces are straight line segments Model with second-order elements—element faces are quadratic curves 单元形状对精度的影响 一阶和二阶四边形单元（六面体单元）的精度随初始单元的扭曲增加而下降，但一次单元对单元扭曲的敏感性小于二次单元。因此当单元形状扭曲严重时，一次单元可能是更好的选择。 二阶三角形单元（四面体单元）对单元的扭曲敏感性较低。 建议：采用二次单元模拟应力集中问题，在应力集中区域尽量使用形状较规则的单元 ideal okay bad distorted undistorted 应力集中问题举例一 应力集中模拟结果讨论 粗网格 精细网格 CPS3 55.06 76.87 CPS4 71.98 91.2 CPS4I 63.45 84.37 CPS4R 43.67 60.6 CPS6 96.12 101.4 CPS8 91.2 100.12 CPS8R 92.56 97.16 s yy (D处目标值=100.0) 单元类型 一次单元模拟精度较差，在精细网格的情况下仍不能获得较精确的结果。 二次单元结果较理想。 应力集中部位网格划分 采用二次单元时，应力集中部位尽量采用形状较规则的单元进行划分，同时推荐减缩积分以减少运算量。 有限元建模专题二 金朝海 jch666@ 有限元建模专题二 6.1 一致性的单位系统 6.2 单元选择 6.3 网格划分 6.4 应力集中问题 6.1 一致性的单位系统 Quantity SI SI (mm) (mm) Length m mm mm Force N N KN Mass kg tonne (103 kg) kg Time s s ms Stress Pa (N/m2) MPa (N/mm2) GPa(KN/mm2) Energy J mJ (10–3 J) J Density kg/m3 tonne/mm3 kg/mm3 6.2 单元选择 单元类型 实体单元 结构单元 插值阶数（节点数） 一次单元（线性）、二次单元 积分方式 完全积分 减缩积分 非协调单元 单元类型 实体单元 3维实体单元 六面体 楔形 四面体 2维实体单元 平面应力 平面应变 轴对称 平面应变(plain strain)单元假设离面应变为零，可以用来模拟厚结构。 平面应力(plain stress)单元假设离面应力为零，适合用来模拟薄结构。 轴对称单元