第5章工程优化设计实例分析.pptVIP

7.5工程优优化设计实例解析 1. 数学模型标准形式： min f （X） s.t. AX≤b （线性不等式约束） AeqX=beq （线性等式约束） C(X)≤0 （非线性不等式约束条件） Ceq(X)=0（非线性等式约束） Lb ≤X ≤Ub （边界约束条件） Fmincon 函数使用的约束优化算法都是目前比较普适的有效算法： （1）对于中等的约束优化问题fmincon使用序列二次规划（SQP sequential quadratic programming）算法； （2）对于大规模约束优化问题fmincon 使用基于内点反射牛顿法的信赖域算法(subspace trust region method and is based on the interior-reflective Newton method)； （3） 对于大规模的线性系统使用共轭梯度算法（PCG preconditioned conjugate gradients）。由于这些算法都具有一定的复杂性，具体算法这里不再详述。 2.使用格式： [x,fval,exitflag,output, grad,hessian]= fmincon(@fun,x0,A,b,Aeq,beq,Lb,Ub,’Nlc’,options) 实例一圆柱螺旋弹簧的优化设计 圆柱螺旋弹簧的结构参数包括 例一调压弹簧为普通圆柱螺旋压缩弹簧。阀腔直径为42mm，弹簧最大工作压力为Fmax=1110N，弹簧的许用应力[ ]=665MPa,弹簧的最大刚度kmax=24N/mm,弹簧结构要求：工作圈数n≥0，弹簧指数4≤ C ≤ 14，弹簧压并高度λb=1.1h=18.25mm。试在满足弹簧的强度条件、刚度条件、稳定性条件、旋绕比条件和结构边界等约束条件下，确定弹簧的簧丝直径d、中径D2和工作圈数n等三个设计参数，使它的结构重量最轻。 例一调压弹簧为普通圆柱螺旋压缩弹簧。阀腔直径为42mm，弹簧最大工作压力为Fmax=1110N，弹簧的许用应力[ ]=665MPa,弹簧的最大刚度kmax=24N/mm,弹簧结构要求：工作圈数n≥0，弹簧指数4≤ C ≤ 14，弹簧压并高度λb=1.1h=18.25mm。试在满足弹簧的强度条件、刚度条件、稳定性条件、旋绕比条件和结构边界等约束条件下，确定弹簧的簧丝直径d、中径D2和工作圈数n等三个设计参数，使它的结构重量最轻。 例一调压弹簧为普通圆柱螺旋压缩弹簧。阀腔直径为42mm，弹簧最大工作压力为Fmax=1110N，弹簧的许用应力[ ]=665MPa,弹簧的最大刚度kmax=24N/mm,弹簧结构要求：工作圈数n≥0，弹簧旋绕比4≤ C ≤ 14，弹簧压并高度λb=1.1h=18.25mm。试在满足弹簧的强度条件、刚度条件、稳定性条件、旋绕比条件和结构边界等约束条件下，确定弹簧的簧丝直径d、中径D2和工作圈数n等三个设计参数，使它的结构重量最轻。 弹簧刚度约束 综上所述，该问题的数学模型为： ①设计变量： ②目标函数： ③约束条件： 2、优化方法和计算结果 调用MATLAB优化工具箱函数fmincon进行优化计算 1）编制目标函数的M文件(Spring_f) 2）编制约束函数的M文件(Spring_g) 3）在命令窗口编写调用主程序(Spring_c) x0=[5,10,10]; options=optimset(LargeScale, off ) ; [x,fval,exitflag,output]=fmincon(@ Spring_f,x0,[],[],[],[],[],[],@ Spring_g,options); disp *********弹簧优化设计最优解********** fprintf(1, 簧丝直径 d=%3.4f\n,x(1)); fprintf(1, 弹簧中径 D=%3.4f\n,x(2)); fprintf(1, 弹簧圈数 n=%3.4f\n,x(3)); 实例二螺栓组连接的优化设计 螺栓组联接的设计计算，主要根据被联接机械设备的载荷大小、功能要求和结构特点，确定螺栓的个数，规格和布置方式。螺栓组联接的优化设计，可以在