liewenbian.m

来自「内燃机转子仿真」· M 代码 · 共 62 行

clear all;
clc;
G=5.363*1e10;
d=0.21;
e=0.3;
if e<0.5&e>=0
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)/64;
Ipn=pi*d^4/32;
Ip=Ipn-dIp;
GIpn=G*Ipn;
GIp=G*Ip;
dGIp=G*Ipn-G*Ip;
elseif e==0.5
    dIp=pi*d^4/64;
    Ipn=pi*d^4/32;
    Ip=Ipn-dIp;
    GIpn=G*Ipn;
    GIp=G*Ip;
    dGIp=G*Ipn-G*Ip;
elseif e<1&e>0.5
    e1=1-e;
    ee=(4*e1-4*e1^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*e1)^4*(ce2-ce1+ce1^3/3-ce2^3/3)/64;
Ipn=pi*d^4/32;
GIpn=G*Ipn;
Ip=dIp;
GIp=G*Ip;
dGIp=G*Ipn-G*Ip;
else
    fprintf('input error');
end
lt=0.3;
%el=
if e>0
l=2*e*d;
l1=lt-l;
k=GIpn/lt;
k1=GIpn/l1;
kc=GIp/l;
kk=kc*k1/(kc+k1);
dk=k-kk;
elseif e==0
    dk=0;
    k=GIpn/lt;
end