OpenMDAO飞艇体积优化代码多类异常问题求助

问题背景 #

基于OpenMDAO开展飞艇最小体积优化，核心设定：

• 唯一设计变量：体积
• 约束条件：浮力比、螺旋桨效率（np）、两种重量计算方法的差值
• 优化目标：最小化体积
• 代码核心类：初始重量计算类、气动特性类、二次重量估算类

三类问题及排查方案 #

1. 体积偶尔超出2.0E6上限 #

• 核对设计变量边界配置：确认add_design_var是否明确指定upper=2.0E6，部分优化器（如SLSQP）迭代后期可能因步长控制小幅越界，可尝试开启bound_enforcement='scalar'参数强制边界约束
• 验证约束逻辑冲突：若体积边界与浮力比等约束存在矛盾，优化器可能为满足约束暂时越界，需检查约束条件的合理性

2. 优化报错LSQ子问题中矩阵E奇异，无法找到可行解 #

• 排查约束线性相关性：若三个约束（浮力比、np、重量差值）存在线性依赖，会导致矩阵奇异。可暂时移除某一约束测试是否能正常运行，或计算约束雅可比矩阵的秩验证相关性
• 统一变量与约束量级：体积（1e6量级）与其他约束参数（如重量差值、螺旋桨效率）量级差异过大，易引发数值奇异。尝试对设计变量和约束做归一化处理（如体积除以1e6，约束值缩放到[-1,1]区间）
• 检查重量计算逻辑：若两种重量计算方法的差值约束在迭代中恒为0或无变化，会导致雅可比矩阵列向量线性相关，需验证两个重量类的输出是否随体积变化合理波动

3. 出现posx和posy应为有限值错误 #

• 排查气动特性类数值异常：posx/posy大概率来自气动特性类计算，检查是否存在除以0、根号内为负、对数输入为负等情况（如体积过小时，飞艇长度、直径等参数触发无效计算）
• 调整优化初始点：若初始体积设置过小，易引发气动参数计算异常，建议将初始点设为合理区间（如1e6左右），并添加设计变量下限约束
• 添加数值保护逻辑：在气动特性类的关键计算步骤中，对分母、根号内表达式设置最小值限制（如max(x, 1e-6)），避免生成NaN或无穷大值

内容的提问来源于stack exchange，提问作者Simon Boudoux