第6讲 有限元建模专题二.pptVIP

有限元建模专题二 金朝海 jch666@ 有限元建模专题二 6.1 一致性的单位系统 6.2 单元选择 6.3 网格划分 6.4 应力集中问题 6.1 一致性的单位系统 6.2 单元选择 单元类型 实体单元 结构单元 插值阶数（节点数） 一次单元（线性）、二次单元 积分方式 完全积分 减缩积分 单元类型 实体单元 3维实体单元 六面体 楔形 四面体 2维实体单元 平面应力 平面应变 轴对称 关于结果输出 结构单元 壳单元——用来模拟那些一个方向的尺寸(厚度)远小于其他方向的尺寸，并且沿厚度方向的应力可以忽略的结构。 杆件单元——用来模拟那些一个方向的尺寸（长度）远大于其他两个方向的尺寸，并且仅沿长度方向的应力比较显著的构件。 梁单元（Beam） 桁架单元 (Truss)：不能承受弯曲，铰接框架结构 近似模拟缆索、弹簧 实体单元建立有限元模型通常规模较大，尤其对于三维实体单元 如果选用适当的结构单元 (shells and beams) 会得到一个更经济的解决方案 模拟相同的问题，用结构体单元通常需要的单元数量比实体单元少很多 要由结构体单元得到合理的结果需要满足一定要求： the shell thickness or the beam cross-section dimensions should be less than 1/10 of a typical global structural dimension, such as: The distance between supports or point loads The distance between gross changes in cross section The wavelength of the highest vibration mode 插值的阶数 通常插值的阶数由单元采用的节点数目决定。 仅在角点处布置节点的单元，如8节点实体单元，在每一方向上采用线性插值，常常称它们为线性单元或一阶单元。 在每条边上有中间节点的单元，如20节点实体单元，采用二次插值，常常称它们为二次单元或二阶单元。 在每条边上有中间节点的修正三角形或四面体单元，如10节点四面体单元，采用修正的二阶插值，常常称它们为修正的单元或修正的二次单元或二阶单元。 积分选择(模拟弯曲问题) 完全积分 减缩积分 完全积分 减缩积分 关于非协调模式单元 非协调模式单元的目的是克服在完全积分、一阶单元中的剪力自锁问题。由于剪力自锁是单元的位移场不能模拟与弯曲相关的变形而引起的，所以在一阶单元中引人了一个增强单元变形梯度的附加自由度。这种对变形梯度的增强允许一阶单元在单元域上对于变形梯度有一个线性变化。 优点 使用非协调单元模拟弯曲问题 (CPS4I, …) 对于模拟以弯曲变形为主的问题或许是最高效且低成本的一种实体单元 综合了一阶和二阶单元的优点 正确模拟剪切行为—纯弯曲变形下不存在剪切变形. 厚度方向只需一个单元即可模拟弯曲问题. 不存在沙漏模式。 在单元过分扭曲的情况下，模拟精度下降（低于一次减缩积分单元）；然而， 所有单元在单元形状过分扭曲的情况下都将下降。 单元扭曲的影响：悬臂梁弯曲 总结 6.3 网格划分 应用足够细密的网格以保证模拟的结果具有足够的精度是非常重要的。粗糙网格可能会产生不精确的结果。随着网格密度的增加，模拟分析所产生的数值结果会趋于一个唯一解．但是运行模型所需耍的计算机资源也会增加。 学会判断网格细分到何种程度所得的结果是可以接受的。当进一步细分网格所得到的解的变化可以忽略不计时，可以说网格已经收敛了。 数值解是所模拟的物理问题的近似解答，取决于模型的几何形状、材料特性、边界条件等描述的准确程度。 建议 对复杂的三维模拟而言，实际上可利用的计算机资源常常限制了所采用的网格密度。在这种情况下，必须非常小心地应用从分析中得到的结果。粗糙的网格足以用来预测趋势。然而，应用由粗糙网格计算得到的位移和应力的具体量值应该很谨慎。 一般来说，没有必要对所分析的结构全部采用均匀的细划网格。应该在出现高应力梯度的地方采用细网格，而在低应力梯度或应力量值不被关注的地方采用粗网格。例如，细划的网格仅应用在高应力梯度的区域，而在其他地方则采用粗网格。例如： 关于应力奇异 应力奇异的处理 如果离感兴趣区域较远，可以在查看结果时通过不激活受影响的区域忽略它的影响 如果位于感兴趣区域，需要如下纠正: 在尖角处增加倒角重新进行分析 代替点力载荷为等效压力载荷 “散布” 位移约束至一个节点集 6.4 应力集中问题 二次单元处理应力集中问题，明显优于一次单元 无论是完全积分还是减缩积分都可以很好的反映应力集中 减缩积分效率更高，而且计算结果往往优于完全积分。 二次单元可以以更少的单元更