技术特征：
1.一种自适应调整预报可行域的变刚度加筋圆柱壳屈曲优化方法，其特征在于，包括以下步骤：构建变刚度加筋圆柱壳屈曲优化模型，以最大化最小屈曲压力为目标函数：所述优化模型包含目标函数、约束函数和设计变量取值范围；使用拉丁超立方实验设计方法获得指定数量初始样本点并真实评价，更新可行最优解；使用kriging代理模型近似目标函数和约束函数，使用目标函数代理模型预报值置信下界作为目标补点准则指导补点；确定候选补点区域：初始化为约束代理模型预报可行域，并利用补点信息和代理模型不确定性信息调整候选补点区域；更新kriging代理模型，搜索可行解；若没有发现可行解，则使用差分进化算法，基于代理模型的总约束违反作为补点准则，最小化该准则获取补点直至获得可行解；若已找到可行解，则通过目标函数代理模型lcb指导目标优化的补点函数进行补点，改进目标值；计算每次迭代补点真实评估的目标值和约束值，并判断补点是否可行；如果补点不可行，则利用kriging代理模型不确定性收缩候选补点区域，否则，扩大候选补点区域以探索其他潜在可行域；反复迭代，更新代理模型直至获得可行最优解，即变刚度加筋圆柱壳的结构参数。2.根据权利要求1所述的一种自适应调整预报可行域的变刚度加筋圆柱壳屈曲优化方法，其特征在于，所述变刚度加筋圆柱壳屈曲优化模型，如公式(1)所示：其中，p
e
为最小屈曲压力，σ1为轴向跨中外壳中面周向应力，σ2为外壳内表面肋骨位置最大轴向应力；σ
f
为肋骨周向应力；w为结构总质量，w0为允许的设计结构总质量；σ
s
＝600mpa为材料屈服极限；圆柱壳半径为3000mm，肋骨间距为500mm，周向总长度为18000mm，泊松比为0.3，弹性模量为2.1
×
105mpa，材料密度为7.85
×
103kg/m3。3.根据权利要求1所述的一种自适应调整预报可行域的变刚度加筋圆柱壳屈曲优化方法，其特征在于，所述通过目标函数代理模型lcb指导目标优化的补点函数进行补点，如公式(2)所示：
其中，为目标代理模型预报值，为目标代理模型预报值不确定性度量，k1为用以平衡目标开发和探索程度的参数，k2为用以调整可行域的参数；δ为等式约束放松量，为预先设定值；使用目标函数的置信下界作为目标补点准则指导补点，并考虑约束代理模型不确定性：当k2＞0时，获得一个较大的候选补点区域，使得其他潜在可行域获得补点机会；当k2＜0时，获得一个较小的候选补点区域，以提高获得的补点是可行解的概率。4.根据权利要求1所述的一种自适应调整预报可行域的变刚度加筋圆柱壳屈曲优化方法，其特征在于，所述基于代理模型的总约束违反作为补点准则，最小化该准则获取补点直至获得可行解，包括：采用基于代理模型的总约束违反作为补点准则，最小化该准则获取补点直至获得可行解，如公式(3)所示：计算各个约束代理模型预报值考虑模型不确定性后的约束违反值，求和得到总约束违反值，最小化该近似总约束违反值，获得补点；其中，参数k3为用以平衡约束开发和探索程度的参数。5.根据权利要求1所述的一种自适应调整预报可行域的变刚度加筋圆柱壳屈曲优化方法，其特征在于，所述差分进化算法，具体包括以下步骤：等概率随机选择de/rand/1和de/rand-to-best/1两种方式中的一种产生变异向量，然后使用公式(6)将变异向量和父代个体交叉产生最终后代：de/rand/1：de/rand-to-best/1：交叉:
其中，对于目标个体x
k
，从种群中随机选择三个不同的个体用以辅助产生后代；x
best
为当前找到的最优解，f和cr分别为放缩因子和交叉概率；产生后代后，进行环境选择时，将该后代个体和其对应父代个体比较，较好的个体遗传至下一代；在处理约束时，使用可行性准则比较个体优劣：针对两个个体x1和个体x2，如果如下条件之一满足：x1和x2都为可行解，且f(x1)＜f(x2)；x1为可行解但x2为不可行解；x1和x2都为不可行解，且cv(x1)＜cv(x2)，则x1优于x2。6.根据权利要求1所述的一种自适应调整预报可行域的变刚度加筋圆柱壳屈曲优化方法，其特征在于，还包括：限制新增补点到历史数据点的距离不小于历史数据点之间最小距离的一半。
技术总结
本发明提供自适应调整预报可行域的变刚度加筋圆柱壳屈曲优化方法，包括以下步骤：构建变刚度加筋圆柱壳屈曲优化模型；使用Kriging代理模型近似目标函数和约束函数，使用目标函数代理模型预报值置信下界作为目标补点准则指导补点；确定候选补点区域：初始化为约束代理模型预报可行域，并利用补点信息和代理模型不确定性信息调整候选补点区域；反复迭代，更新代理模型直至获得可行最优解，即变刚度加筋圆柱壳的结构参数。当约束函数和目标函数维数较低，非线性程度适中时，本方法具有较高效率，具有实际应用价值。具有实际应用价值。具有实际应用价值。
技术研发人员：刘均 曾勇 程远胜 张攀 陈长海 周奇
受保护的技术使用者：华中科技大学
技术研发日：2022.10.25
技术公布日：2023/1/31