将对应代码保存为.m文件，运行计算即可。相关过程请在参考文献中查阅
主函数 main.m
%汽车设计-离合器参数优化
%2017.4.3
%清除屏幕；清除函数；关闭启用
clc;
clear;
close;
%参数输入
nemax=5200; %最大功率转速r/min
F=5895.34; %摩擦面单位压力
D=207.5; %摩擦片外径
d=160; %摩擦片内径
%调用主函数
b=[70;-125]; %外径线性不等式约束70，内径线性不等式约束125
lb=[0,0,0]; %变量下限
A=[0,(pi*nemax*10^(-3))/60,0;0,0,-1]; %系数
%A=[0,0.2618,0;0,0,-1];
x0=[F,D,d]; %变量初值
%x0=[4661,200,130];
options=optimset('Algorithm','interior-point'); %其他参数设置：全局优化
%options=optimset('Algorithm','sqp');
[x,fval]=fmincon(@lhqobj,x0,A,b,[],[],lb,[],@lhqys,options)
目标函数 obj.m
%汽车设计-离合器参数优化
%目标函数
function lhqobj=lhqobj(x)
lhqobj=pi/4*(x(2)^2-x(3)^2);
约束函数 ys.m
%汽车设计-离合器参数优化
%2017.4.3
function[c,ceq]=lhqys(x)
%参数输入
Z=2; %摩擦片数
f=0.3; %摩擦因数
Temax=202; %发动机最大扭矩
ma=1470; %汽车总质量/kg
rr=0.2333; %滚动半径 /m
ig=3.5; %变速器一档传动比
i0=3.67; %主减速器传动比
ne=3500; %发动机转速 乘用车1200r/min
btx=1.20; %后备系数β的下限 乘用车1.2~1.75
bts=1.75; %后备系数β的上限
w=0.4; %单位摩擦面积滑磨功ω的上限
Tc0=280; %单位摩擦面积传递的转矩许用值Tc0的上限
%非线性约束条件
c(1)=-(x(3)/x(2))+0.53; %c的下限0.53
c(2)=(x(3)/x(2))-0.70; %c的上限0.70
c(3)=-(Z*f*x(1)*(x(2)^3-x(3)^3)/(Temax*10^3*3*(x(2)^2-x(3)^2)))+btx; %β的下限1.20
c(4)=(Z*f*x(1)*(x(2)^3-x(3)^3)/(Temax*10^3*3*(x(2)^2-x(3)^2)))-bts; %β的上限1.75
c(5)=-(4*x(1))/(pi*(x(2)^2-x(3)^2))+0.25; %p0的下限0.25
c(6)=(4*x(1))/(pi*(x(2)^2-x(3)^2))-0.35; %p0的上限0.35
c(7)=(pi^2*ne^2*ma*rr^2/(1800*i0^2*ig^2))*4/(pi*2*(x(2)^2-x(3)^2))-w; %ω的上限0.40
c(8)=(4*f*x(1)*(x(2)^3-x(3)^3))/(3*pi*((x(2)^2-x(3)^2)^2))-Tc0; %Tc0的上限280
ceq=[]; %非线性等式为空
参考文献
[1].王望予.汽车设计[M]. 北京：机械工业出版社，2007
本文地址：http://www.8mpa.com/213.html
版权声明：本文为原创文章，版权归八兆帕所有，欢迎分享本文，转载请保留出处！
