liewenbian1.m

内燃机转子仿真

clear all;
clc;
G=8.08*1e11;
d=0.030;
e=0.1;
ee=(4*e-4*e^2)^0.5;
eee=asin(ee);
ce1=cot(0.5*pi-eee);
ce2=cot(0.5*pi+eee);
dIp=d^4*eee/32+d^4*(1-2*e)^4*(ce2-ce1-ce1^3/3+ce2^3/3)/32;
Ipn=pi*d^4/32;
Ip=Ipn-dIp;
GIpn=G*Ipn;
GIp=G*Ip;
dGIp=G*Ipn-G*Ip;