复杂结构有限元分析,网格划分该用哪种方法?
在有限元数值模拟流程中,离散化处理构成分析工作的前提性环节,同时也是消耗人力资源与时间成本最为显著的前处理操作。离散单元的空间排布品质直接调控数值解的精确程度与运算资源的消耗水平,在涉及几何非线性显著的情形下,甚至还关乎迭代过程的收敛稳定性。
自动自由离散技术
自由离散技术代表了当前自动化干预程度最高的单元生成策略之一。该技术应用于二维流形(包括平面区域与曲面)时,可自主产生三角形或四边形单元;应用于三维实体区域时,则自动填充四面体单元。在实际操作中,借助Hypermesh软件的二维面板内的automesh指令集,能够高效完成对面域及其单元分布的自动生成任务。
对于几何特征呈现高度复杂性的分析模型,自动离散路径的主要优势体现在时间成本的低廉性;然而其局限性同样显著,即单元总量往往偏大,且难以实现单元密度的空间差异化调控——在需要稀疏网格的区域与需要精细网格的区域之间,使用者通常缺乏直接有效的控制手段。为此,有必要对几何面域实施分区策略,从而获得符合分析者预期且具备较高计算效率的单元配置方案。
针对三维复杂实体,该策略目前仅能生成四面体单元,但其处理效率极为突出,仅需设定若干关键控制参数即可获得品质良好的网格。不过,单元总数严重依赖于几何模型的最小细节尺度,通常导致总体单元规模庞大。为提升数值计算的经济性,需对几何构型进行预处理操作,包括局部加密等处理,这与二维情形下的分区策略逻辑相通。
映射型离散技术
映射型离散技术针对拓扑规整的几何模型提供了一种结构化单元生成方案。其原始定义可表述为:在二维情形下,仅适用于四边形域,且对边需设定相等的分割数目,生成单元全数为四边形;在三维情形下,仅适用于六面体域,对应棱边与对应面的分割数须保持一致,生成单元全数为六面体。
随着软件算法的发展,当前主流离散工具已显著放宽了上述约束条件,具体表现为:
面域可接受三角形、四边形或其他任意多边形拓扑。
面域对边的分割数允许不同取值,但存在特定兼容性约束。
面域支持生成全三角形映射网格。
体域可接受四面体、五面体、六面体或其他任意多面体拓扑。
体域对应棱边与面的分割数允许不同取值,同样存在特定兼容性约束。
针对三维复杂几何实体,工程中的通行做法是利用几何切割工具将其分解为若干四、五或六面体子区域,随后对各子区域分别执行映射离散。应当指出,这种纯粹的映射处理流程较为烦琐,需投入较多时间与精力,但其回报是可获得品质优良的网格。
拖拉与扫略离散技术
对于通过拉伸、旋转或偏置等操作生成的三维复杂实体,可先在原始截面构建壳单元形式的面网格,进而在实体生成过程中自动形成三维体网格。对于已存在的三维复杂实体,若其在某一方向的拓扑结构保持恒定,则可利用扫略离散功能生成网格。上述两种方式所产生的单元几乎全部为六面体单元。
在Hypermesh软件的三维面板中,solidmap命令集提供了多种拖拉与扫略的实现形式,包括从单元至面、从面至面等路径,且拉伸方式选项丰富。使用者可根据具体几何特征灵活选用。一般而言,扫略路径是一种优质的离散策略,对于复杂几何实体,经过适度的切割预处理后,即可自动生成规整的六面体网格。相较于映射离散,扫略方法展现出更强的适应性与操作灵活性。在实际工程中,若目标为获取全六面体单元,通过几何分区配合扫略功能是最主要的实现手段。
在ANSA环境下,情况具有相似性。ANSA作为一款以几何特征为核心的先进离散工具,其面构建功能尤为出色。该软件不依赖实体(Body)概念即可实现模型分区,操作流程简洁且效率突出,这一设计理念代表了未来离散软件发展的主流趋势。
混合离散策略
混合离散策略是指在单一几何模型上,根据不同区域的几何特征,分别选用自由、映射、扫略等多种离散方式,以求达到综合效能最优的有限元模型构建目标。该策略的决策过程需在计算精度、运算耗时与建模工作量三者之间进行权衡。
工程实践中,为提升计算精度并压缩求解时间,应优先对适合扫略与映射离散的区域进行六面体单元划分,并尽可能借助切割等布尔操作手段创建适宜离散的子区域(尤其针对关注的核心区域)。对于经反复切割仍无法适用结构化离散、只能采用四面体自由离散的区域,建议采用带有中间节点的四面体单元进行自由网格划分。
自由度耦合与约束方程技术
针对某些特殊形式的复杂几何构型,可利用ANSYS提供的约束方程与自由度耦合功能(在Abaqus中对应Tie约束)辅助生成优质网格并有效控制计算规模。
具体操作思路为:对相邻的实体分别进行独立的网格划分(通常采用映射或扫略方式),随后通过“粘结”处理将其连接。由于各实体在几何层面相互独立,无需考虑彼此间网格的匹配性,因此可灵活运用多种手段为各子区域生成高质量的网格。实体间的“粘结”通过形函数插值实现自由度耦合,从而保证连接界面处位移连续性的精确满足。若连接区域的应力场为分析重点,可进一步采用下文所述子模型技术在局部建立精细化分析模型。
子模型技术及其他辅助手段
子模型技术是一种遵循“先全局、后局部”分析流程的数值策略(亦称切割边界条件法)。对于仅需获取局部区域精确解的复杂几何模型,该技术能够以最小的额外工作量获得目标区域的准确结果。
其操作流程为:首先生成整体分析模型,忽略模型中诸如倒角、开孔、开槽等细部几何特征(依据圣维南原理,局部细节的改动对整体分析结果影响有限),并在该整体模型上划分较为粗糙的网格(计算与建模成本均较低),施加荷载并完成初步求解。随后,在与整体模型相同的坐标系下建立局部子模型,在此模型中纳入先前忽略的细部特征并划分精细网格(需注意子模型的切割边界应远离所关注区域),再次进行求解计算。
该技术的附加优势包括:可在子模型基础上方便地对所关注的细部特征进行尺寸优化或参数调整(如改变圆角半径、缝隙宽度等);整体模型与子模型可采用不同的单元类型,例如整体模型使用板壳单元,子模型使用实体单元等。
此外,合理利用结构的对称性在实际工程中具有显著效益:一方面可大幅削减计算规模,另一方面有利于准确施加边界条件。航空发动机涡轮盘的分析即为典型应用场景。对于几何构型与载荷条件均具备轴对称、循环对称或面对称特征的问题,应优先考虑利用对称性进行简化。
总体而言,针对复杂几何模型,综合运用多种技术手段构建高品质、高计算效率的有限元模型,是数值模拟全流程中最基础且最关键的第一步。需在长期实践中不断积累经验、进行验证比对,方能深入掌握复杂模型的网格划分技术。
评论
莫名其妙

干货推荐

研究生有限元复习重点
爽子
68 6

血脑屏障超微结构
Mr弘🔬
118 9

lammps中文教程5
春风得意马蹄疾
48 2

免疫荧光染色全流程
7330a6d4
34 1

【代码分享】费米面嵌套(Fermi surface nesting)函数计算
张阳小号1
651 2

模态分析到底在算啥?搞懂这个,结构为什么会自己振动就明白了
莫名其妙
64 4
