目标函数:function fun=myobj(x) %以质量大小为目标的优化函数 p=0.00785; %材料密度,单位Kg/cm?
n=3; %簧片数目
bb=7.5; %簧片宽度,单位cm
fun=2*p*bb*n*(x(1)*x(3)+0.5*(x(4)-x(3)-4.2)*(x(1)+x(2))+x(2)*4.2); %目标函数
非线性约束:
function[c,ceq]=mycon(x) %约束函数
Fs=8986/3; %板簧端部作用力F/3,这里假设各片板簧均载,单位N
bb=7.5; %板簧宽度,单位cm
Kh=992500/3; %板簧设计刚度(装配刚度),单位N/m
e=0.008; %刚度允许误差
w=4.2; %U型螺栓中心距的一半,单位cm
n1=300; %L1段许用应力,单位Mpa
n2=420; %L2段许用应力,单位Mpa
E=2.06*10^9; %材料弹性模量
h=0.92; %修正系数
c(1)=6*Fs