弹簧优化.docVIP

[例三] 试设计一压缩圆柱螺旋弹簧，要求其质量最小。弹簧材料为65Mn，最大工作载荷Pmax=40N，最小工作载荷为0，载荷变化频率fr=25Hz，弹簧寿命为104h，弹簧钢丝直径d的取值范围为1-4mm，中径D2的取值范围为10-30mm，工作圈数n不应小于4.5圈，弹簧旋绕比C不应小于4，弹簧一端固定，一端自由，工作温度为50C，弹簧变形量不小于10mm。 本题的优化目标是使弹簧质量最小，圆柱螺旋弹簧的质量可以表示为 式中， --弹簧材料的密度，对于钢材 =7.8×10-6kg/mm3； n—工作圈数； n2—死圈数，常取n2=1.5-2.5，现取n2=2； D2—弹簧中径，mm； d—弹簧钢丝直径，mm。 将已知参数代入公式，进行整理以后得到问题的目标函数为 根据弹簧性能和结构上的要求，可写出问题的约束条件： 1． 强度条件 2． 刚度条件 3． 稳定性条件 4． 不发生共振现象，要求 5． 弹簧旋绕比的限制 6． 对d,n,D2的限制 且应取标准值，即1.0,1.2,1.6,2.0,2.5,3.0,3.5,4.0mm等。 由上可知，该压缩圆柱螺旋弹簧的优化设计是一个三维的约束优化问题，其数学模型为： 取初始设计参数为X(0)=[2.0,5.0,25.0]T 首先编写目标函数的M文件opt25_3.m，返回x处的函数值f。 function f = myfun(x) f = 0.192457*1e-4*( x(2)+2)*x(1)^2 * x(3); 由于约束条件中有非线性约束，所以需要编写一个描述非线性约束条件的M文件opt25_3c.m： function [c,ceq] = mycon(x) c(1)=350-163*x(1)^(-2.86)*x(3)^0.86; c(2)=10-0.4*0.01*x(1)^(-4)*x(2)*x(3)^3; c(3)=(x(2)+1.5)*x(1)+0.44*0.01*x(1)^(-4)*x(2)*x(3)^3-3.7*x(3); c(4)=375-0.356*1e6*x(1)*x(2)^(-1)*x(3)^(-2); c(5)=4-x(3)/x(1); 然后设置线性约束的系数： A=[-1 0 0 1 0 0 0 –1 0 0 1 0 0 0 –1 0 0 1]; b=[-1;4;-4.5;50;-10;30]; 下一步给定初值，给定变量的下限约束，并调用优化过程(磁盘中M文件为opt25_3.m) x0 = [2.0; 5.0; 25.0]; lb=zeros(3,1); [x,fval,exitflag,output,lambda]=fmincon(@opt25_3o,x0,A,b,[],[], lb,[],@opt25_3c) 计算结果为： x = 1.6564 4.5000 16.1141 fval = 0.0055 exitflag = 1 output = iterations: 3 funcCount: 16 stepsize: 1 algorithm: medium-scale: SQP, Quasi-Newton, line-search firstorderopt: [] cgiterations: [] 所以当弹簧钢丝的直径d、工作圈数n及中径D2分别取1.6564、4.5000和16.1141时弹簧质量最小，为5.5克。考虑到实际情况，各参数可分别取1.6、5.0和16.0 x=[1.6 5.0 16.0]; z=optim253(x) z = 0.0055 故此时弹簧质量仍为5.5克。 fmincon函数 功能：求多变量有约束非线性函数的最小值。 题目： 内燃机用气门弹簧。最大变形量mm,工作载荷F=680N，工作频率Hz，最高工作温度150c。寿命。弹簧的结构要求满足：钢丝直径，外径，工作圈数，旋绕比，刚度。设计要求在满足强度、刚度和其他性能限制的条件下选择钢丝直径d、簧圈中径和工作圈数n使弹簧的重量最轻 目标函数：弹簧重量W作为目标函数 式中 ---材料密度，钢材 ----死圈数，=1.5~2，取=1.8 目标函数 强度条件 弹簧的剪切应力计算公式为：式中-----许用应力， K------曲度系数， 因此有约束式： 初始状态确定：由动力学仿真分析得到 按结构经验公式推荐尺寸关系为：