2. 逆向模型参数约束的数学基础
在逆向建模中，我们往往先得到一组“看起来能拟合观测数据”的参数，但真正可用的模型参数必须同时满足物理意义、几何形态和工程可制造性等多方面要求。把这些要求写成数学语言，就构成了参数约束。本节从参数空间出发，依次建立约束集合、可行域的概念，并说明等式约束与不等式约束在逆向问题中的几何与代数含义，为后续的约束求解与参数优化打下基础。
一、参数空间：把模型参数看作一个点
设逆向模型共有 n 个待定参数，记为向量 x = (x₁, x₂, …, xₙ)ᵀ。所有可能的参数取值构成 n 维实数空间 ℝⁿ，称为参数空间。参数空间中的每一个点对应一组完整的模型参数，也就对应一个候选模型。例如一个简化的梁截面反演模型，参数为弹性模量 E、截面惯性矩 I、边界转动刚度 kr，则参数空间是三维的，点 (E, I, kr) = (30 GPa, 2.5×10⁻⁴ m⁴, 8×10⁷ N·m/rad) 就代表一个具体候选模型。
需要强调的是，参数空间本身只是“容器”，它不区分参数是否合理。E 取负值、I 取零，在 ℝⁿ 中都是合法点，但在物理上毫无意义。约束的作用，就是从整个参数空间中筛出真正可接受的子集。
二、约束集合与可行域
把每条约束写成关于 x 的等式或不等式，所有约束的公共部分就是约束集合：
等式约束： h_i(x) = 0, i = 1, 2, …, m
不等式约束：g_j(x) ≤ 0, j = 1, 2, …, p
可行域： Ω = { x ∈ ℝⁿ | h_i(x)=0, g_j(x)≤0 }可行域 Ω 是参数空间的一个子集，逆向问题的解必须落在 Ω 内。若 Ω 为空集，说明约束互相矛盾，问题无解，此时应回到约束建模阶段检查是否把要求写得过紧；若 Ω 是一个点，则参数被完全确定，无需优化；若 Ω 是一个有界区域，则问题转化为在区域内寻找使目标函数最小的点。
三、等式约束与不等式约束的含义
等式约束在几何上把可行域“压”到一张低维曲面上。例如在二维参数平面中，一条等式约束 h(x₁, x₂) = 0 通常是一条曲线，可行解只能沿曲线移动，自由度减少一维。代数上，等式约束意味着参数之间存在确定的函数关系，可用于消元：若能从 h(x)=0 中解出 x₂ = φ(x₁)，就把 n 维问题降为 n−1 维。逆向问题中典型的等式约束包括：几何闭合条件（如截面面积等于给定值）、平衡方程残差为零、归一化条件（如 Σwᵢ = 1）等。
不等式约束则在几何上“切”出半空间或区域边界，可行域保留在边界的一侧。它不降低维数，但会限制搜索范围，并可能在最优点处起作用（称为有效约束）。例如要求截面惯性矩 I ≥ Imin、材料参数 0 < E ≤ Emax、制造公差 |t − t₀| ≤ δ，都是不等式约束。判断某条不等式在最优点是否有效，是逆向结果分析的重要环节：有效约束往往对应“参数被顶到边界”，提示该约束可能是限制拟合精度的瓶颈。
四、从几何到代数：约束如何进入求解
理解约束的几何含义之后，还要把它转化为可计算的代数形式。逆向问题通常写成带约束的优化形式：
min f(x) = Σ ( y_obs,k − y_model,k(x) )²
s.t. h_i(x) = 0
g_j(x) ≤ 0其中 f(x) 是观测残差平方和，yobs,k 为第 k 个观测值，ymodel,k(x) 为对应模型响应。求解时，等式约束可用拉格朗日乘子并入目标函数，不等式约束则通过 KKT 条件刻画：在最优点处，每条不等式约束要么不起作用（乘子为零），要么恰好取等号（约束有效且乘子非负）。这一条件正是“最优点是否被边界卡住”的代数表达，与前面几何上的判断完全对应。
实际实现中，常见做法是先用变量变换把不等式约束“吸收”掉。例如要求 E > 0，可令 E = eu，对 u ∈ ℝ 无约束优化；要求 a ≤ x ≤ b，可令 x = a + (b−a)·σ(v)，其中 σ 为 Sigmoid 函数。这样处理后，搜索过程天然落在可行域内，避免迭代点越界导致模型失效。但要注意，变换会改变参数空间的度量，梯度下降的步长含义也随之变化，需相应调整学习率或采用信赖域策略。
五、一个可操作的约束建模示例
假设反演某构件截面参数，参数为宽度 b、高度 h、腹板厚度 t，观测到截面面积 Aobs = 4.8×10⁻³ m²，并已知构造要求。可建立如下约束集（数值为教学示例假设）：
| 约束类型 | 数学表达 | 工程含义 |
|---|---|---|
| 等式 | b·h − (b−2t)(h−2t) = 4.8×10⁻³ | 截面面积与观测一致 |
| 不等式 | 0.15 ≤ b ≤ 0.40 | 宽度在可选型材范围内 |
| 不等式 | 0.20 ≤ h ≤ 0.60 | 高度满足净空要求 |
| 不等式 | 0.006 ≤ t ≤ 0.020 | 板厚满足焊接与稳定下限 |
| 不等式 | h/b ≤ 3.0 | 控制截面高宽比 |
建模后应先做可行性检查：把面积等式与上下限联立，估算是否存在满足全部条件的 (b, h, t)。若发现面积要求过大而厚度上限过小，则可行域可能为空，此时不应强行求解，而应复核观测值或放宽某条构造要求。确认可行后，再以面积等式消去一个变量（例如解出 t 关于 b、h 的表达式），把三维问题降为二维带界约束问题，最后用序列二次规划或投影梯度法求解，并检查最优点处哪些不等式取等号，据此判断是哪个构造要求限制了反演精度。