中国石油大学（华东）机械原理课程设计程序.pdfVIP

中国石油大学（华东）机械原理课程设计程序

#include

#include

main()

{

FILE*fp;

inti,j,k,p;

doublew1,PIB,ab,G1,G2,G3;

double

l1=0.063,l2=0.333,lk=0.111,PI=3.1415926,m1=40,m2=26,m3=1

5.6,g=9.80665,J2=0.25;double

Q1[37],Q2[37],w2[37],sc[37],vc[37],aa2[37],ac[37];

double

a2x[37],a2y[37],a2[37],pi2x[37],pi2y[37],MI2[37],pi3[37],pi1x[

37],pi1y[37],R43[37],Mb[37],Mbz[37],mrz[37],mr[37];

doubleR12[37],R12x[37],R12y[37],a12[37],b12[37];

doubleR23[37],R23x[37],R23y[37],a23[37],b23[37];

doubleR41[37],R41x[37],R41y[37],a41[37],b41[37];

doubleWr[37],Wd[37],WW[37];

doubleMd,WWmax,WWmin,EEmax,JF;

doubledt=0.2;

w1=2*150*PI/60;/*--转速150转/分--*/

ab=-l1*w1*w1;/*--切向加速度--*/

PIB=-m1*ab;/*--求惯性力Pib--*/

G1=m1*g;G2=m2*g;G3=m3*g;G1=m1*g;G2=m2*g;G3=m3*g;

{{

}

for(i=0;i37;i++)

{

/*--运动分析--*/

Q1[i]=i*PI/18;/*--每次增加10度--*/

Q2[i]=asin(-l1*sin(Q1[i])/l2);/*--l2sin(Q2[i])=-l1*sin(Q1[i])--

*/

sc[i]=l1+l2-l1*cos(Q1[i])-l2*cos(Q2[i]);

w2[i]=(-l1*cos(Q1[i])*w1)/(l2*cos(Q2[i]));

vc[i]=-l1*sin(Q1[i])*w1-l2*sin(Q2[i])*w2[i];

aa2[i]=(l1*sin(Q1[i])*w1*w1+l2*sin(Q2[i])*w2[i]*w2[i])/(l2*co

s(Q2[i]));/*--求阿尔法2--*/

ac[i]=-l1*cos(Q1[i])*w1*w1-l2*cos(Q2[i])*w2[i]*w2[i]-

l2*sin(Q2[i])*aa2[i];

/*--动力分析--*/

if(Q1[i]=0Q1[i]=PI)p=0;/*--吸入冲程受力为0--*/

elsep=3500;/*--压缩冲程受力3500--*/

a2x[i]=l1*w1*w1*cos(PI+Q1[i])+lk*(ac[i]-

l1*w1*w1*cos(PI+Q1[i]))/l2;

a2y[i]=l1*w1*w1*sin(PI+Q1[i])+lk*(-

l1*w1*w1*sin(PI+Q1[i]))/l2;

a2[i]=sqrt(a2x[i]*a2x[i]+a2y[i]*a2y[i]);

pi2x[i]=-m2*a2x[i];

pi2y[i]=-m2*a2y[i];

MI2[i]=-J2*aa2[i];

pi3[i]=-m3*ac[i];

pi1x[i]=PIB*cos(Q1[i]);

pi1y[i]=PIB*sin(Q1[i]);

R43[i]=((pi3[i]-p)*l2*sin(Q2[i])+pi2x[i]*lk*sin(Q