图6.3上止点前360°偏心轴转角计算网格初始网格建完后，转子转一圈，偏心轴转三圈，动态网格随着时间的变化发生变形及运动，为了保证计算能够正常进行，需要不断对网格进行更新。通过宏命令DEFINE_GRID_MOTION和DEFINE_CG_MOTION定义网格的运动规律，网格更新方式采用弹簧光顺和网格重构相结合的方式。......
2023-06-23
网格模型检查及收敛性分析：有限元法应用实践
【摘要】：在网格划分完毕后，需要对网格模型进行检查。图3.6有限元装配模型表3.2材料参数图3.6有限元装配模型表3.2材料参数3.收敛性分析有限元法是一种数值分析方法，因此需要考虑收敛性问题。分别定义5种网格尺寸：2.5mm、2mm、1.5mm、1mm、0.6mm，在热-机耦合工况下对网格模型进行收敛性分析。
1.单元分析
网格质量的好坏不仅关系到能否快速、顺利地完成有限元分析，还关系到计算分析的结果与精度。在进行网格划分时，首先需要考虑单元类型，需要针对所分析的问题，结合选取的分析方法，选择合适的单元类型，在保证计算结果精度的同时，避免过度地增加计算工作量。
本次有限元分析是针对转子发动机关键零部件的热机耦合分析。实现耦合分析通常有直接耦合与顺序耦合两种分析方法。针对这两种分析方法对应有不同的单元类型。直接耦合分析是指多个物理场的自由度要同时进行计算，适用于多个物理场各自响应相互依赖的情况。但由于平衡状态需要满足多个准则，每个节点上的自由度多，矩阵方程庞大，计算时间长，直接耦合依赖于其耦合多个场自由度的单元，从理论上说，每一组场组合就要提供对应类型的单元。而顺序耦合分析是指以特定的顺序求解单个物理场模型，再将前一个物理场的分析结构作为后续分析的边界条件。顺序耦合主要适用于物理场之间单向的耦合关系。顺序耦合比直接耦合分析效率高，并且不需要选择特殊的单元类型。由于本次热机耦合中温度场的分布会对结构应力场、应变场分布产生较大影响，而应力场对温度场的影响较小，并且考虑到分析效率以及分析的通用性，选择顺序耦合方法作为本次分析方法。先分析结构件温度场，再将温度场结果作为应力场分析的边界条件，进行应力场分析。因此在分析过程中会涉及2种类型的单元，即热分析单元与应力分析单元。
单元类型的确定还与单元形状有关。在确定分析方法后，还需根据模型的实际结构及计算要求选择合适的单元形状。常用的单元形状有四面体单元与六面体单元。六面体网格在计算精度、变形特性、划分的网格数量、抗畸变能力都比四面体网格具有更加明显的优势。因此本次分析将对模型进行六面体的网格划分。在ABAQUS软件中，不同分析方法、分析步骤中选择的六面体网格类型，如表3.1所示。2.网格划分及质量检查
表3.1 不同分析方法所选择的单元类型
采用有限元前处理软件HyperMesh对转子发动机关键零部件（气缸、前后端盖、转子）进行六面体网格划分。在网格划分完毕后，需要对网格模型进行检查。如果形状较差的网格所在单元位置远离分析所关心的重要位置时，一般情况下对分析结果的影响不大，可以忽略。但当形状较差的网格出现在重要位置，则会造成计算收敛困难。对于转子发动机，重点关注部位为螺栓孔、气缸电热塞附近区域、燃烧热弧区。本次网格检查的基本要求如下：
（1）检查网格的形状扭曲比和细长比：0.1＜扭曲比的平均值＜0.75，细长比＞0.1；
（2）检查节点和单元的重合度，不应存在重合度小于0.1mm的节点或单元。
针对检查出现问题的单元，可以通过移动、合并、删除节点，使得单元质量达到要求。
本次分析所建立的六面体网格模型如图3.5所示。
图3.5 有限元模型
（a）前端盖网格；（b）后端盖网格；（c）转子网格；（d）气缸网格
转子发动机的前后端盖及气缸是通过螺栓连接到一起，在11个螺栓孔中，有两个是定位销孔，它是气缸加工和装配的基准，因此对定位销孔的刚度、加工精度以及表面粗糙度的要求都非常严格。在两个定位销孔处的螺栓是通过螺母拧紧的，而其他9个螺栓是通过与后端盖上的螺纹连接来达到连接的目的。建立的转子发动机装配体有限元网格模型如图3.6所示。
材料是有限元分析的基础，不同材料在相同载荷作用下会有不同的力学响应。本次分析模型中气缸、前后端盖材料为QT600，转子材料为QT800，各材料的材料属性如表3.2所示。
图3.6 有限元装配模型
表3.2 材料参数
3.收敛性分析
有限元法是一种数值分析方法，因此需要考虑收敛性问题。当网格逐渐加密时，有限元的解答就逐渐趋近于精确解。但随着网格数量的增多，计算成本会逐渐加大。因此在保证计算精度的前提下，需要对多种网格大小的转子发动机模型进行对比计算分析，以确定合理的网格大小以及局部加密要求。具体方法为：使用当前的网格密度和单元类型进行分析，再使用更细化的网格进行分析。将分析结果进行比较，如果结果相差较小（一般可取2%），则说明该网格可以满足分析精度需求。
分别定义5种网格尺寸：2.5mm、2mm、1.5mm、1mm、0.6mm，在热-机耦合工况下对网格模型进行收敛性分析。取气缸上3个关键位置的应力值进行分析，参考点1、2、3分别位于燃烧室区域螺栓孔附近、燃烧区域端面以及燃烧区域中部应力变化较大的区域（图3.7），计算结果如图3.8所示。根据图3.8中结果分析，可以发现随着网格尺寸的减小，各参考点的热-机耦合应力趋于稳定。比较1mm与0.6mm 网格的计算结果，可以发现3个参考点在这两次计算中的差异分别为1.1%、2.1%、1.6%，能够满足分析精度的要求。为了保证更好的计算精度，后续分析将选用尺寸为0.6mm的网格模型作为计算分析模型。
图3.7 收敛性分析参考点位置
图3.8 参考点应力随网格大小的变化
有关小型转子发动机数值建模与仿真分析的文章
- 详细阅读
-
SAMCEF有限元分析网格划分与应用详细阅读
单击Mesh模块图标，进入网格划分模块，相应的工具条中显示出网格处理模块快捷命令图标。1．网格划分确认数据树中Model下的Wirel已经被选中，或者直接在屏幕上单击线Wirel。设置Average Mesh Size为10，单位为mm。采用同样的方法将Used Data中的项目隐藏。单击鼠标右键，在弹出的快捷菜单中选择Mesh Aspect—S→rink命令，如图14-54所示。图14-54 选择Shrink命令此操作是收缩显示单元，便于看到单元内部的情况。......
2023-10-27
-
扩展有限元法的应用与发展详细阅读
传统有限元法一般采用多项式或插值多项式描述单元内部的位移场，由于形函数的连续性，导致单元内部不能存在间断。XFEM 以标准有限元的理论为框架，在包含裂纹间断的区域，对形函数进行扩展，加入带不连续性质的形函数来描述裂纹间断引起的位移变化。......
2023-06-27
-
网格划分技术应用-SAMCEF有限元分析与应详细阅读
单击Mesh模块图标，进入网格划分模块，相应的工具条中显示出网格处理的快捷命令图标。设置Average Mesh Size为10，单位为mm。单击Apply按钮，生成一维梁单元，总计69个梁单元。在数据树Used Data中再次选中Beam_on_2_edge_of_Wire1和Beam_on_1_edge_of_Wire1_#1，单击右键，在弹出的快捷菜单中选择Display Data Aspect→Show命令，显示的转子形状如图14-278所示。单击Save按钮保存文件。......
2023-10-27
-
无网格计算方法：伽辽金法与有限元法的比较详细阅读
鉴于这种缺陷，近几年来国际上许多著名的计算力学学者，如T.Belytschko、O.C.Zienkiewicz、S.N.Atluri、J.T.Oden、W.K.Liu等都对无网格方法表现出了极大的兴趣，并进行了大量的研究工作。无网格伽辽金法与有限元法的相似之处在于，两者都是将边值问题等价为一个条件变分问题，然后由条件变分问题通过数值积分离散为代数方程组。......
2023-07-02
-
SAMCEF有限元分析：网格划分与应用实例详细阅读
单击Mesh模块图标，进入网格划分模块，相应的工具条中显示出网格处理的命令。设置Average Mesh Size为10，单位mm。单击Element Type图标，在弹出的单元类型对话框中选择二次单元parabolic。单击Apply按钮，生成三维实体单元。单击鼠标右键，弹出快捷菜单，选择Mesh Aspect→Shade命令，显示出图14-165所示的有限元网格。单击Save按钮保存文件。......
2023-10-27
-
网格划分模块及操作步骤，有限元网格展示详细阅读
图14-14所示为Mesh模块。图14-14 网格划分模块及其快捷命令图标1.网格划分确认数据树中Model下的Wire1已经被选中，或者直接在屏幕上用鼠标单击线Wire1。设置Average Mesh Size为20。采用同样的方法将Used Data中的项目隐藏。此操作是收缩显示单元，便于看到单元内部的情况。在工具条上单击等轴侧3D着色显示图标，屏幕左上角弹出一个图形窗口，显示图14-16所示的有限元网格。图14-15 选择Shrink命令图14-16 有限元网格3.保存文件单击工具条中保存文件的快捷图标。......
2023-10-27
-
有限元分析，网格划分实例详细阅读
图8-37 网格划分模块及快捷命令图标1.网格划分确认数据树中Model下的impeller已经被选中，或者直接在屏幕上用鼠标左键单击叶轮几何模型。在弹出的图8-38所示对话框中，设置Average Mesh Size为20。单击Apply按钮，生成三维单元。在数据树中选中Model下的Directory，右击，弹出图8-41所示的快捷菜单，选择Hide命令，将几何图形隐藏。在工具条上单击等轴侧3D着色显示图标 ，屏幕左上角弹出一个图形窗口，显示出图8-42所示的三维有限元模型。......
2023-10-27
相关推荐