铝法兰各应力须满足所规定的强度要求

       为确保密封可靠，整体铝法兰各应力须满足所规定的强度要求，即出现在铝法兰高颈大端和小端的轴向应力、环向应力、径向应力组合应力应分别小于相应的许用应力，故建立以下约束条件:仁二势一铝法兰材料在操作温度下的许用应力，MPa;Mo一铝法兰总力矩，由此可知，铝法兰优化设计是一个单目标二维非线性约束规划问题，其数学模型写成如下标准形式对上述非线性约束优化问题的求解方法较多，应用场合不尽相同。特别是早期的优化方法，需要通过构造惩罚函数将有约束优化转换.为无约束问题来求解。现在这些方法已经被更有效的基于K-T(Kuhn-Tuker)方程解的方法所取代。其算法直接计算拉格朗日乘子，通过拟牛顿法更新过程，给K-T方程积累二阶信息，以保证有约束拟牛顿法的超线性收敛。由于在每次迭代中都求解一次二次规划问题，故该方法称为序列二次规划法(SQP法)。

       在MATLAB优化工具箱中，用SQP法求解非线性约束优化问题，由fmincon函数来实现，调用格式为式中:fun为目标函数的M文件;x0为给定的初始值;A,b分别为线性不等式约束的系数，若没有则取A=[ ].b=[ ];Aeq}Beq分别为线性等式约束的系数，若没有则取Aeq=[]}beq=[];lb,ub分别为设计变量的下界和上界，若没有界，可设lb=[],ub=[];nonlcon参数，用于计算非线性不等式约束。(X)。和非线性等式约束ceq (X)=0。当对应的函数采用M文件表示，即nonlcon=' mycon'，则M文件。ycon.m具有以下形式%计算x的非线性不等式ceq=%计算x的非线性等式x为设计变量的返回值;fval为目标函数最优值;exitflag为迭代终止条件;output为运行详细结果;lambda为求解x处拉格朗日乘子的lambda参数;grad为求解x处fun函数的梯度;hessian为求fun处函数的Hess矩阵。http://www.zblansheng.cn/