2. 参数化建模的数学基础

参数化建模的本质,是把一个几何形状或仿真模型从“固定尺寸的实体”改写为“由一组可调参数决定的函数”。一旦模型被这样改写,设计就不再是反复重画,而是在参数空间中搜索。理解这一过程,需要掌握四个数学对象:设计变量、约束条件、目标函数,以及它们共同构成的参数空间,再加上参数到几何的映射关系。

设计变量:把“尺寸”变成“坐标”

设计变量是模型中允许被调整、并用于驱动几何或物理属性的参数。例如一个带圆角的矩形板,可以取长度 L、宽度 W、圆角半径 R、板厚 t 作为设计变量。写成向量形式:

x = [L, W, R, t]ᵀ = [120, 80, 10, 3]ᵀ   (单位: mm)

这个向量 x 就是设计点在参数空间中的坐标。需要注意,设计变量并不等于 CAD 中的全部尺寸。只有那些被明确声明为“可变”的尺寸才进入 x;其余尺寸作为常量或由表达式派生。例如孔位间距可以写成 L/2,这样它随 L 自动变化,而不额外占用一个设计变量。区分“自变量”和“派生量”,是控制问题规模的第一步。

约束条件:参数空间中的可行域边界

约束条件限定设计变量可以取值的范围。常见形式有三类:

  • 边界约束(显式上下限):R ≤ W/2,t ≥ 1.5 mm,L ∈ [80, 200];
  • 几何约束(避免自交、保证可制造):R < min(L, W)/2,孔边距 ≥ 2t;
  • 性能约束(由仿真结果反推):最大应力 σ_max ≤ 235 MPa,一阶固有频率 f₁ ≥ 120 Hz。

前两类可以直接写成关于 x 的不等式,第三类必须经过一次 CAE 求解才能得到,属于隐式约束。所有满足约束的 x 构成可行域,记作 Ω。优化只能在 Ω 内进行,因此约束既是“护栏”,也是决定搜索空间形状的关键。

目标函数:把设计意图写成可计算的量

目标函数 f(x) 是衡量设计优劣的标量。轻量化设计常取质量最小:

f(x) = ρ · (L·W − (4−π)R²) · t
     = 7.85e-6 · (120·80 − (4−π)·10²) · 3
     ≈ 0.223 kg   (示例计算,ρ 单位 kg/mm³)

若同时关心刚度和频率,可构造加权目标 f = w₁·m/m₀ + w₂·δ/δ₀,其中 δ 为最大变形。目标函数必须可计算、可比较,且对 x 连续或至少分段连续,否则梯度类算法难以使用。

参数空间与几何映射:从 x 到模型

参数空间是设计变量张成的 n 维空间,每个点对应一个候选设计。几何映射 Φ 把参数向量 x 映射为具体的 CAD/CAE 模型:

Φ: x → 几何 G(x) → 网格 M(G) → 响应 y = S(M)

其中 S 表示求解器。映射链中任何一环不连续,都会让优化变得困难。例如圆角半径 R 变化时,如果拓扑发生改变(圆角与孔相交),网格数量和节点位置会跳变,响应 y 可能出现台阶。因此参数化建模要尽量保持拓扑稳定,让 Φ 在可行域内光滑。

参数空间 x x = [L, W, R, t]ᵀ L W 可行域 Ω x₀ 约束: R ≤ W/2, σ_max ≤ 235 MPa Φ 几何 G(x) L=120, W=80, R=10, t=3 网格 M 与求解 S σ_max, f₁, δ 响应 y = S(M) 反馈至目标函数与约束判断 目标函数 f(x) min m = ρ·A·t 约束条件 σ_max ≤ 235 MPa, f₁ ≥ 120 Hz 优化搜索 在 Ω 内更新 x → 迭代

参数如何真正驱动模型变化

在 CAD 中,参数驱动通常通过尺寸约束和表达式实现。以矩形板为例,若把 L 从 120 改为 150,而孔位表达式为 L/2,则孔位自动从 60 移到 75,圆角 R 若保持 10 不变,几何拓扑不变,网格可以映射复用。这种“参数变、拓扑不变”的映射称为保拓扑参数化,是优化迭代中最理想的情况。

在 CAE 中,参数变化还会影响网格和边界条件。例如板厚 t 从 3 变为 4,若采用壳单元,厚度是截面属性而非几何尺寸,网格无需重划,求解速度很快;若采用实体单元,则需要重新划分网格,计算成本上升。因此选择参数化策略时,要同时考虑几何映射和网格映射的稳定性。

实用判断:如果一个参数变化会导致模型拓扑改变(如孔消失、圆角合并),应尽量把它设为离散变量或分段处理,而不是直接放入连续优化。否则响应面会出现不连续,梯度信息失真。

一个可执行的小例子

假设要对上述矩形板做轻量化,设计变量取 L、W、R、t,初始值 [120, 80, 10, 3]。约束为:R ≤ W/2,σ_max ≤ 235 MPa,f₁ ≥ 120 Hz。目标为质量最小。操作流程如下:

  1. 在 CAD 中建立参数化草图,标注 L、W、R、t,并写孔位表达式 L/2;
  2. 在 CAE 中建立静力分析和模态分析,把 σ_max 和 f₁ 设为输出响应;
  3. 用实验设计(DOE)在可行域内采样,例如 L ∈ [80,200]、W ∈ [50,120]、R ∈ [2,20]、t ∈ [1.5,6];
  4. 拟合响应面或直接调用优化算法,在约束下最小化质量;
  5. 检查最优解是否满足约束,并验证网格无关性。

示例中若某次迭代得到 L=150、W=70、R=8、t=2.5,计算质量约为 0.196 kg,σ_max=210 MPa,f₁=128 Hz,则说明该点可行且优于初始设计。注意这些数值是教学示例,实际结果取决于材料、载荷和边界条件。

常见误区

  • 把派生尺寸也设为设计变量,导致变量间强相关,搜索空间冗余;
  • 约束给得过松,最优解落在不可制造区域;约束给得过紧,可行域为空;
  • 目标函数包含不可计算量,或量纲不统一,导致加权系数难以确定;
  • 忽略几何映射的连续性,优化过程频繁失败或结果跳变。

掌握设计变量、约束、目标函数和参数空间的关系,就掌握了参数化建模的“坐标系”。后续的 DOE、响应面、梯度优化和拓扑优化,都是在这个坐标系中寻找更优设计点的具体方法。