tanhuang2.m

弹簧隔振器的系数k0=4.5 hc=3.5 d=[3.15 3.23] a1=tan(pi*(d)./(2*hc)) a2=cot(pi*(d)./(2*hc)) ％ plot

k0=1.;
hc=4;

d=[1.:0.01:2];
%d=[3.96 4.22]

a1=tan(pi*(d+0.25)./(2*hc));
a2=cot(pi*(d+0.25)./(2*hc));

y=(a1-a2)*2*k0*hc./pi;
y=y./9.8*8;
plot(d,y)

var=(tan(pi*(d+0.25)./(2*hc))).^2;
bi=(var+1+1+1./var)
k=k0*1E3.*(var+1+1+1./var);
figure
plot(d,bi)
% figure
% plot(d,y)
d1=solve('2*1.755*4/pi*(tan(pi*(x+0.25)/(2*4))-cot(pi*(x+0.25)/(2*4)))=0.98/8')
d2=solve('2*1.755*4/pi*(tan(pi*(x+0.25)/(2*4))-cot(pi*(x+0.25)/(2*4)))=11*0.98/8')
var2=(tan(pi*([d1(2),d2(2)]+0.25)./(2*hc))).^2;
k=k0*1E3.*(var2+1+1+1./var2);