曲线图.doc.docVIP

P=3.39 R=1.06 K =4.4640 C=2.82 A=3.96 I=2.82 a=asin(I/K) syms x x2=2*pi-x-0.013*pi+a L=sqrt(R*R+K*K-2*K*R*cos(6.9238-x)) b=asin(R*sin(6.9238-x)/L) x3=acos((P*P+L*L-C*C)/2/P/L)-b x31w=diff(x3) x4=acos((P*P-L*L-C*C)/2/P/L)- b x41w = diff(x4) v=acos((C*C+L*L-P*P)/2/C/L)+b v1=diff(v) je=39664.32356+41.25*(sin(x4-x2)/sin(x3-x4))^2+1870.14*(sin(x3-x2)/sin(x3-x4))^2+678.37*(x31w)^2-848.47* x31w*cos(x3-x4) je1=diff(je) Sc=4.924-3.69*v F1=75000*Sc+25000 F2=40000 F3=75000*Sc-193250 F4=25000 Z=x3-x2-0.013*pi-x Me1=-3.96*F1*x41w+106487.29*sin(0.013*pi+x)+2452.81*sin(0.013*pi+x)+3922.19*Z*x31w Me2=-3.96*F