Abaqus网格大变形之CEL(欧拉法)
仿真侠客
2025年08月31日 17:41
收录于文集
共1篇
CAE仿真分析

在Abaqus中,一般处理大变形问题的方法主要包括以下几种:

当上述进行模型仿真时,上述几种方法都力不从心时,CEL提供了一种全新的思路。

核心思想:让材料的归材料,网格的归网格

欧拉(Eulerian)域:一个固定的、在空间中不动的“容器”网格。材料在这个网格的单元中流入、流出、流动。就像你用一个固定的筛子去筛流动的沙子,筛子不动,沙子动。

拉格朗日(Lagrangian)域:传统的随材料变形的结构体网格,比如一个螺栓、一个压头。

耦合(Coupled):通过通用接触(General Contact) 或欧拉-拉格朗日接触,让运动的 Lagrangian 部件与在 Eulerian 域中流动的材料发生相互作用。

简而言之,CEL用欧拉网格来模拟极度变形的材料(如流体、土壤、粉末、需大变形部件),用拉格朗日网格来模拟正常变形的结构体,二者通过接触耦合。 这样,欧拉材料无论怎么流动,其网格永远不变形,完美解决了畸变问题;而结构体则能享受Lagrangian方法在接触和精度上的优势。

CEL的典型应用场景

1. 金属切削、钻孔、铣削:切屑的形成和分离是天然的材料流动。

2. 结构入水、液舱晃荡、船体兴波:水用欧拉域,船体用拉格朗日域。

3. 土壤开挖、碾压、 penetration:土体被视为可流动的材料。

4. 食品加工、药片压制成型:粉末、颗粒材料的流动与压实。

5. 密封圈润滑、齿轮搅油:分析润滑油的行为。

 

使用CEL的关键设置与注意事项(实战精华)

1. 欧拉域的创建与尺寸设定

创建:在Part模块创建部件时,选择Eulerian类型。

尺寸:欧拉域就是一个“水池”,你要确保这个水池足够大,能够容纳整个变形过程中流动的材料,当存在 Lagrangian 部件时,要能覆盖Eulerian部件Lagrangian 部件的接触面。宁大勿小,但过大会增加计算成本。

CEL的典型应用场景

1. 金属切削、钻孔、铣削:切屑的形成和分离是天然的材料流动。

2. 结构入水、液舱晃荡、船体兴波:水用欧拉域,船体用拉格朗日域。

3. 土壤开挖、碾压、 penetration:土体被视为可流动的材料。

4. 食品加工、药片压制成型:粉末、颗粒材料的流动与压实。

5. 密封圈润滑、齿轮搅油:分析润滑油的行为。

 

使用CEL的关键设置与注意事项(实战精华)

1. 欧拉域的创建与尺寸设定

创建:在Part模块创建部件时,选择Eulerian类型。

尺寸:欧拉域就是一个“水池”,你要确保这个水池足够大,能够容纳整个变形过程中流动的材料,当存在 Lagrangian 部件时,要能覆盖Eulerian部件Lagrangian 部件的接触面。宁大勿小,但过大会增加计算成本。

2. 材料赋予与初始体积分数(EC - Eulerian Volume Fraction)

这是CEL设置中最核心、最容易出错的概念!

EC:它定义了每个欧拉单元被某种材料填充的百分比(1=完全填充,0=完全空,0.5=填充一半)。

初始状态:你必须通过预定义场(Predefined Field) 来告诉Abaqus,在计算开始前,欧拉域的哪些区域被材料填充。通常使用“均匀(Uniform)” 或通过“离散点(Discrete Points)” 选择一个初始几何体来定义。

材料界面:材料与空气(或另一种材料)的界面不是清晰的网格边界,而是通过单元内EC值从1到0的过渡区域来表征。因此,网格密度决定了界面的分辨率。

3. 材料定义

欧拉域内的材料(如水、土壤)本构模型通常需要包含失效或塑性行为,以允许其像流体一样流动(例如,土壤可用Drucker-Prager模型,金属可用Johnson-Cook模型和失效准则)。

4. 接触设置

必须使用通用接触(General Contact)。

接触属性中, tangential behavior通常设置为“Frictionless”或很小的摩擦系数,因为欧拉材料与Lagrangian部件的相互作用机制比较复杂,过大摩擦可能导致不稳定。

欧拉域本身不需要定义接触。

5. 网格划分

欧拉域:通常使用结构化网格(Hex-dominated),如六面体单元EC3D8R。因为材料在规则网格中流动更稳定,计算效率更高。(注意:欧拉体不能划分四面体网格,尽量使用结构化网格,一提高计算精度及效率)

拉格朗日部件:按常规方式网格划分即可。

6. 结果输出与后处理

关心的结果不再是单元的变形(因为欧拉网格不变形),而是:

    材料分布:在Visualization模块,通过结果→场输出→EVF(Eulerian Volume Fraction) 来查看不同材料的分布和流动过程。这是观察现象的关键。

 应力应变:应力应变云图只显示在有材料(EC>0) 的单元上,空气区域(EC=0)没有意义。

补充:

CEL的成本:CEL计算量通常很大,因为欧拉域往往节点数很多,且是多材料计算。显式分析(Abaqus/Explicit)是使用CEL的首选求解器,因为它能更好地处理高度非线性的接触和材料失效。

与SPH的抉择:对于高速冲击、破碎、溅射问题,SPH可能更合适。对于连续流动、与结构有复杂相互作用的问题(如压铸成型、切削、入水),CEL通常是更好的选择。

思维的转变:使用CEL最大的挑战是从传统的“跟踪材料变形”思维,转变为“观察材料在固定网格中流动”的思维。

CEL是Abaqus武器库中一把强大而 specialized 的利器。它并非万能,但对于特定的一类极端大变形问题,它能提供其他方法无法企及的解决方案。希望这篇文章能成为你开启CEL仿真大门的钥匙。

(完)