利用Matlab求解机械设计优化问题-螺栓.docxVIP

3．机械优化设计应用实例 机械优化设计把数学规划理论与数值方法应用于设计中， 用计算机从大量可行方案中找出最优化设计方案， 从而大大提高设计质量和设计效率。 MATLAB 具有解决线性规划和非线性规划、约束 优化和无约束优化问题的内部函数，因而可以完成这一功能。 现举一例： 螺栓组联结的优化设计 如图 4 所示的压力容器螺栓组联接中，已知 D1= 400mm,D2 = 250mm，缸内工作压力为 p=1.5 MPa，螺栓材料为 35 号钢， σs=320Mpa,安全系 S=3,取残余预紧力 Q’p=1.6F,采用铜皮石棉密封垫片。现从安全、可靠、经济的角度来选择螺栓的个数 n 和螺栓的直径 d 。 3．1 设计问题分析 若从经济性考虑， 螺栓数量尽量少些、 尺寸小些，但这会使降低联结的强度和密封性，不能保证安全可靠的工作；若从安全、可靠度考虑，螺栓数量应多一些、尺寸大一些为好，显然经济性差，甚至造成安装扳手空间过小，操作困难。为此，该问题的设计思想是：在追求螺栓组联结经济成本最小化的同时，还要 保证联结工作安全、可靠 。 3．2 设计变量 目标函数 约束条件 3．2 .1 设计变量 选取螺栓的个数 n 和直径 d(mm)为设计变量 : X [n d ]T [x 1 x 2 ]T 3．2. 2 目标函数 追求螺栓组联结经济成本 Cn 最小为目标。而当螺栓的 长度、材料和加工条件一定时，螺栓的总成本与 nd 值成正比，所以本问题优化设计的目标函数为 min F(X) = Cn = n d = x1x2 ① 强度约束条件 为了保证安全可靠地工作，螺栓组联结必须满足强度条 件 ca 5.2Q ] ； 其中 [ ] s 320 [ 106Mpa ； 3 d12 S 3. Q Q p F 1.6F F 2.6F 2.6 D 22 p 2.6 1.5 250 2 60937 n N； 4n 4n 对于粗牙普通螺纹：由文献 [3] 推荐，小径 d1=0.85d 所以，强度约束条 件为： 1 105624 146192 106 146192 0 g1 ( X ) 106 nd 2 x1 x22 106 nd12 ② 密封约束条件 为了保证密封安全 ， 螺栓间距应小于 10d， 所以，密封约 D1 10 d 400 10 x 2 0 束条件为 ： g 2 ( X ) n x1 ③ 安装扳手空间约束条件 为了保证足够的扳手空间， 螺栓间距应大于 5d， 所以，安装约束条件为 ： g 3 ( X ) 5d D1 5x2 400 0 n x1 ④ 边界约束条件 g 4 ( X ) x1 0 ； g5 ( X )x2 0 3．3．3 建立数学模型 综上所述，本问题的数学模型可表达为： 设计变量 ： X [x 1 x 2 ]T 目标函数： min F(X) = x1x2 约束条件： s.t. g i ( X ) 0 ( i = 1, 2, 3, 4, 5,) 现运用 MATLAB 的优化函数进行求解 : 先编写 M 文件 function [c,ceq]=mynas(x) c(1)=146192/(x(1)*x(2)^2)-106; % 非线性不等式约束 c(2)=400*pi/x(1)-10*x(2); c(3)=-400*pi/x(1)+5*x(2); ceq=[]; % 非线性等式约束 在MATLAB 命令窗口输入 : fun=x(1)*x(2); % 目标函数 x0=[4,6]; % 设计变量初始值 A=[-1,0;0,-1]; % 线性不等式约束矩阵 b=[0;0]; Aeq=[]; % 线性等式约束矩阵 beq=[]; lb=[]; % 边界约束矩阵 ub=[]; 2 [x,fval]=fmincon(fun,x0,A,b,Aeq,beq,lb,ub,@mynlsub) % 调用有约束 优化函数 运行结果如下 : x = 11.4499 10.9751 fval = 125.6637 所以，该问题优化结果为： n =11.4499 ,d = 10.9751 ，目标函数最小 值： F(X)= 125.6637 。根据实际问题的意义取整、标准化： n =12 ,d = 12 。 由此例可以看出，与其它编程语言相比， MATLAB 语言可以简化编程。 5是调用 MATLAB 绘图函数自动对上例绘制的数学模型要素图（标注数字的曲线为目标函数的等值线） ，为此 在MATLAB 命令窗口输入 : x1=0.1:20; y1=146192./(106.*x1.^2); y2=400.*pi./(10.*x1); y3=400.*pi./(5.*x1); plot(y1,x1,y2,x1,y3,x1,