第4章:几何建模与网格划分

各位工程师朋友,今天我们来聊聊真空灌注模拟里最基础、也最容易被忽视的一步——几何建模与网格划分。

说实话,我见过太多人一上来就急着跑仿真,结果模型导进去网格画不出来,或者画出来算到一半就发散。嗯,这种坑我踩过不止一次。所以这一章,咱们把这块地基打扎实。

4.1 CAD模型导入——别让模型成为瓶颈

我个人习惯用SolidWorks或CATIA建好几何,然后导出为STEP或IGES格式。为什么?因为这两个格式在大多数CFD软件里兼容性最好。

但有个细节要注意:导出前一定要做几何清理。什么意思?就是把那些小圆角、小倒角、螺栓孔、螺纹这些对流动没影响的特征统统去掉。你想想看,网格划分时这些细节会让网格量暴增,计算时间翻倍,结果却没啥变化。

我的经验:导入后第一件事,用软件的几何修复工具检查一下。我曾经有一次导入一个机翼模型,表面有个0.1mm的缝隙没发现,结果网格怎么都画不过去,折腾了两天。后来发现是导出时精度设置太低。

导入后,建议做这几步检查:

  • 检查是否有自由边、重叠面
  • 确认单位制统一(毫米还是米?)
  • 查看模型是否封闭(流体域必须封闭)
  • 简化不必要的细节特征

4.2 2D网格划分策略——快速验证的利器

很多人觉得2D模拟不够精确,其实不然。在工艺参数初步调优阶段,2D网格能帮你快速试错。我经常用2D模型先跑一遍,看看树脂流动的大致趋势,再上3D精细算。

2D网格划分,核心就三个字:三角形。为什么?因为三角形网格对复杂几何的适应性最好,而且自动生成算法成熟。

具体策略:

  1. 全局尺寸控制:设置一个基础尺寸,比如2mm。太大会丢失细节,太小网格量爆炸。
  2. 局部加密:在注胶口、出胶口、拐角这些流动变化剧烈的地方,把网格尺寸缩小到0.5mm甚至更小。
  3. 质量检查:画完网格一定要看最小角度和长宽比。三角形的最小角不要低于20°,否则计算容易发散。

小技巧:2D网格画完后,可以先用稳态算一下压力分布。如果压力云图有锯齿状,说明网格质量不够好,需要重新调整。

4.3 3D网格划分策略——精度与效率的平衡

3D网格才是真空灌注模拟的主力。但3D网格的划分,说白了就是在精度和计算量之间找平衡。

我个人推荐用四面体+棱柱层的组合网格。四面体填充内部区域,棱柱层处理壁面附近的边界层。为什么这样搭配?因为四面体自动生成快,棱柱层能准确捕捉壁面附近的流动细节。

来看一个典型的网格参数设置:

参数 推荐值 说明
全局最大尺寸 3-5 mm 根据模型大小调整
局部最小尺寸 0.3-0.5 mm 用于注胶口、拐角处
棱柱层数 3-5 层 边界层网格
棱柱层增长率 1.2-1.5 从壁面向外逐渐增大
网格质量(偏斜度) < 0.85 超过0.9需要重新划分

注意:3D网格的节点数不是越多越好。我曾经有个项目,网格从200万加密到500万,结果计算结果只差了2%,但计算时间多了3倍。所以,先做网格无关性验证,找到那个“够用”的网格密度。

4.4 边界层网格处理——流动模拟的关键

边界层网格,这是很多新手容易忽略的地方。真空灌注中,树脂在纤维层表面的流动行为,很大程度上取决于边界层的捕捉精度。

为什么会这样?因为树脂是粘性流体,在壁面附近速度梯度很大。如果边界层网格不够密,就捕捉不到这个梯度,算出来的流动阻力就会偏小。

我的做法是这样的:

  • 第一层网格高度:用y+值估算。对于真空灌注,y+一般在1-5之间。
  • 棱柱层增长率:控制在1.2以内,太大会导致网格过渡不自然。
  • 棱柱层数:至少3层,我一般用5层。

这里有个经验公式可以估算第一层高度:

Δy = (y+ * μ) / (ρ * u_τ)

其中:
Δy = 第一层网格高度
y+ = 目标y+值(取1-5)
μ = 树脂动力粘度
ρ = 树脂密度
u_τ = 壁面摩擦速度(需要估算)

当然,实际项目中我很少去精确计算这个值。更常用的方法是:先画一套粗网格,算完后查看壁面附近的y+分布,如果太大就加密,如果太小就放松。迭代两三次就能找到合适的参数。

避坑指南:我曾经在一个大型风电叶片项目中,边界层网格只画了2层,结果算出来的树脂填充时间比实际少了30%。后来加了5层棱柱层,结果才和实验对得上。从那以后,我再也不敢在边界层上偷懒了。

4.5 网格质量检查——别让坏网格毁了你的仿真

网格画完了,别急着提交计算。先做质量检查,这是职业习惯。

主要检查这几个指标:

  • 偏斜度(Skewness):越接近0越好,一般要求小于0.85
  • 正交质量(Orthogonal Quality):越接近1越好,一般要求大于0.15
  • 长宽比(Aspect Ratio):边界层网格可以大一些,内部区域控制在5以内
  • 最小体积:不能出现负体积,否则计算直接报错

如果发现坏网格,不要手动一个个去修。我建议的做法是:

  1. 先定位坏网格的位置
  2. 分析原因(几何尖角?尺寸突变?)
  3. 回到几何或网格参数中调整
  4. 重新生成

嗯,这里要提醒一句:不要试图用网格平滑功能来掩盖坏网格。平滑只能改善一点点,根本问题还是在几何或参数设置上。

4.6 本章知识体系

为了让大家更直观地理解这一章的内容结构,我画了一张流程图:

几何建模与网格划分知识体系 CAD模型导入 2D还是3D? 2D三角形网格 快速验证、参数调优 3D四面体+棱柱层 精确模拟、边界层捕捉 边界层网格处理 y+估算、棱柱层设置 网格质量检查 偏斜度、正交质量、长宽比

这张图把整个流程串起来了。从CAD导入开始,根据需求选择2D或3D策略,然后重点处理边界层,最后做质量检查。每一步都环环相扣,缺一不可。

好了,这一章的内容就到这里。网格划分这块,说白了就是多练、多试、多总结。刚开始可能会觉得繁琐,但等你积累了几个项目经验后,就会形成自己的套路。

记住一句话:好的网格是成功仿真的第一步。别在这一步省时间,否则后面会花更多时间在调试和纠错上。

专注资料整理