1、精選優質文檔-傾情為你奉上機械原理大作業（一）平面六桿機構的運動分析 班 級：學 號：姓 名：同 組 者: 完成時間：一題目11 說明如圖所示為一片面六桿機構各構件尺寸如表格1所示，又知原動件1以等角速度=1rad/s沿逆時針方向回轉，試求各從動件的角位移、角加速度以及E點的位移、速度及加速度的變化情況。12 數據組號L1L2L2L3L4L5L6axGyG1-A26.5105.665.067.587.534.425.0600153.541.7表格1 條件數據13 要求三人一組，編程計算出原動件從0360º時（計算點數N=36）所要求各運動變量的大小，并繪制運動線圖及點的軌跡曲線。二解

2、題步驟由封閉圖形ABCD可得：由封閉圖形AGFECD可得于是有：對以上1到4導可得- 寫做矩陣形式： 對上述矩陣求導可得：E點橫坐標及對應的速度和加速度: E點縱坐標及對應的速度和加速度: 三計算程序框圖輸入：l1、l2、l2、l3、l4、l5、l6、xG及yG和a2,a3,a5,a6, a1= (I-1)*100調用牛頓迭代法子程序求解位置(1)求得a2,a3,a4,a5及a6,并計算xE,yE調用系數矩陣A子程序,并計算A調用原動件位置參數列陣B子程序,并計算BB(J)=B(J)w1J=1,N 調用高斯消去法子程序求解速度方程(2),求出w2,w3,w4,w5及w6,再求出vEx及vEy調

3、用計算A子程序,并計算其矩陣DA調用計算B子程序,并計算列陣DBw(1)=w2,w(2)=w3w(3)=w4,w(4)=w5DB(K)=DB(K)w1K=1,NI=l,37II=1,NB(K)= -DA(K, II).w1(II)+DB(K)調用高斯消去法子程序求解加速度方程(3),求出a2,a3,a4,a5及a6,并求出aEx及aEy打印結果結 束四源程序1#include "stdlib.h" #include "math.h" #include "stdio.h" int agaus(a,b,n) int n; dou

4、ble a,b; int *js,l,k,i,j,is,p,q; double d,t; js=malloc(n*sizeof(int); l=1; for (k=0;k<=n-2;k+) d=0.0; for (i=k;i<=n-1;i+) for (j=k;j<=n-1;j+) t=fabs(ai*n+j); if (t>d) d=t; jsk=j; is=i; if (d+1.0=1.0) l=0; else if (jsk!=k) for (i=0;i<=n-1;i+) p=i*n+k; q=i*n+jsk; t=ap; ap=aq; aq=t; if (

5、is!=k) for (j=k;j<=n-1;j+) p=k*n+j; q=is*n+j; t=ap; ap=aq; aq=t; t=bk; bk=bis; bis=t; if (l=0) free(js); printf("failn"); return(0); d=ak*n+k; for (j=k+1;j<=n-1;j+) p=k*n+j; ap=ap/d; bk=bk/d; for (i=k+1;i<=n-1;i+) for (j=k+1;j<=n-1;j+) p=i*n+j; ap=ap-ai*n+k*ak*n+j; bi=bi-ai*n+k

6、*bk; d=a(n-1)*n+n-1; if (fabs(d)+1.0=1.0) free(js); printf("failn"); return(0); bn-1=bn-1/d; for (i=n-2;i>=0;i-) t=0.0; for (j=i+1;j<=n-1;j+) t=t+ai*n+j*bj; bi=bi-t; jsn-1=n-1; for (k=n-1;k>=0;k-) if (jsk!=k) t=bk; bk=bjsk; bjsk=t; free(js); return(1); 2.程序二求解各桿的角度和E點坐標（桿四的角度始終為零，

7、程序中不再求解）/*輸出文件 Output.txt*/#include "stdio.h"#include <conio.h>#include "dnetn.c"#include "agaus.c"#include <math.h>#define PI 3.979#define ANGLE (PI/180)#define ALPHA (PI*60/180)/*初始桿長：*/static double dLen7=65.0, 26.5, 105.6, 67.5, 87.5, 34.4, 25.0;/*G 點的坐標

8、*/static double dGxy2=153.5,41.7;static double x0 = 0 * ANGLE;main()int i,k;double eps,t,h,xe,ye;int nTemp=0;FILE *pf;/* 各角位移的初始估計值 */static doublex4=36.80*ANGLE,69.57*ANGLE,119.00*ANGLE,48.92*ANGLE;pf = fopen("Output.txt","w");t=0.1; h=0.1; eps=0.; k=100;for (nTemp = 0; nTemp &l

9、t;= 36; nTemp+)x0 = nTemp * 10 * ANGLE;i=dnetn(4,eps,t,h,x,k);printf("ni=%dn",i);printf("x0=%lfn",x0);fprintf(pf,"%dt",10*nTemp);for (i=0; i<=3; i+)printf("x(%d)=%13.7lfn",i,xi/ANGLE);fprintf(pf,"%lft",xi/ANGLE);xe=dLen1 * cos(x0)+ dLen2 * cos(x0)

10、+dLen0 * cos(ALPHA-(x0);ye=dLen1 * sin(x0)+ dLen2 * sin(x0)-dLen0 * sin(ALPHA-(x0);printf("xe=%13.7lfnye=%13.7lfn",xe,ye);fprintf(pf,"n");printf("n");fclose(pf);getch();/*建立牛頓法矩陣*/void dnetnf(x,y,n) int n;double x,y;y0 = dLen1 * cos(x0) + dLen2 * cos(x0) - dLen3 * cos(x

11、1)- dLen4 ;y1 = dLen1 * sin(x0) + dLen2 * sin(x0) - dLen3 * sin(x1);y2 = dLen1 * cos(x0) + dLen2 * cos(x0) + dLen5* cos(x2) - dLen6 * cos(x3) - dGxy0 +dLen0 * cos(ALPHA-(x0);y3 = dLen1 * sin(x0) + dLen2 * sin(x0) + dLen5* sin(x2) - dLen6 * sin(x3) - dGxy1-dLen0 * sin(ALPHA-(x0);n = n;return;3. 程序三求解

12、各桿的角速度和E點速度（桿四的角速度始終為零，程序中不再求解）/*輸出文件Output2.txt*/#include "stdio.h"#include "math.h"#include "agaus.c"#include <conio.h>#define PI 3.979#define ANGLE (PI/180)#define ALPHA (PI*35/180)/*初始桿長：依次為*/static double dLen7=65.0, 26.5, 105.6, 67.5, 87.5, 34.4, 25.0;/* G 點

13、的坐標*/static double dGxy2=153.5,41.7;static double x0 = 0 * ANGLE;static double w1 = 1.0;double Vex,Vey;main()int i,j,nTime;FILE *fInput;FILE *fOutput2;int nX0,nCounter;double dTemp4=0;static double a44= 0.;static double b4=0.;/*將agu_01.C 的輸出文件output.txt 作為輸入文件，繼續計算*/fInput = fopen("output.txt&q

14、uot;,"r");fOutput2 = fopen("output2.txt","w");for(nTime=0;nTime<=36;nTime+)fscanf(fInput,"%d",&nX0);x0=nX0*ANGLE;printf("x0= %lfn",x0);for (nCounter=0;nCounter<4;nCounter+)fscanf(fInput,"%lf",&dTempnCounter);for (nCounter=0;nC

15、ounter<4;nCounter+)dTempnCounter = dTempnCounter * ANGLE;/*建立高斯方程組矩陣*/a00 = - dLen2 * sin(dTemp0);a01 = dLen3 * sin(dTemp1);a02 = 0.;a03 = 0.;a10 = dLen2 * cos(dTemp0);a11 = - dLen3 * cos(dTemp1);a12 = 0.;a13 = 0.;a20 = dLen2 * sin(dTemp0)- dLen0 * sin(ALPHA-(dTemp0);a21 =0.;a22 = dLen5 * sin(dTe

16、mp2);a23 = -dLen6 * sin(dTemp3);a30 = dLen2 * cos(dTemp0)+dLen0 * cos(ALPHA-(dTemp0);a31 = 0.;a32 = dLen5 * cos(dTemp2);a33 = - dLen6 * cos(dTemp3);b0 = dLen1 * sin(x0) * w1;b1 = - dLen1 * cos(x0) * w1;b2 = -dLen1* sin(dTemp0)* w1 ;b3 =-dLen1* cos(dTemp0)* w1 ;if (agaus(a,b,4)!=0)for (i=0;i<=3;i+

17、)printf("w%d=%lfn",i,bi);fprintf(fOutput2,"%lf ",bi); Vex=-dLen1 * sin(x0) * w1-dLen2 * sin(dTemp0) * b0+dLen0 * sin(ALPHA-(dTemp0) *b0;Vey=dLen1 * cos(x0) * w1+dLen2 * cos(dTemp0) * b0+dLen0 * cos(ALPHA-(dTemp0) * b0; printf("Vex=%lfnVey=%lfn",Vex,Vey); fprintf(fOutput

18、2,"n");printf("n");fclose(fInput);fclose(fOutput2);getch();4. 程序四求解各桿和點E的角加速度（桿四角加速度始終為零，這里不再求解）/*輸出文件Output3.txt*/#include "stdio.h"#include "math.h"#include "agaus.c"#include <conio.h>#define PI 3.979#define ANGLE (PI/180)#define ALPHA (PI*35

19、/180)/*初始桿長：依次為*/static double dLen7=65.0, 26.5, 105.6, 67.5, 87.5, 47.2, 37.8;/*G 點的坐標*/static double dGxy2=153.5,41.7;static double x0 = 0 * ANGLE;static double w1 = 1.0;main()int i,nTime;FILE *fInput2;FILE *fOutput3;int nX0,nCounter;double Aex,Aey;double dSeta5=0.;double dOmiga4 = 0.;static doubl

20、e a44= 0.;static double b4=0.;/*將agu_02.C 的輸出文件output2.txt 作為輸入文件，繼續計算*/fInput2 = fopen("output2.txt","r");fOutput3 = fopen("output3.txt","w");for(nTime=0;nTime<=36;nTime+)fscanf(fInput2,"%d",&nX0);/*讀入角度值*/for (nCounter=0;nCounter<5;nCount

21、er+)fscanf(fInput2,"%lf",&dSetanCounter);/*讀入角速度值*/for (nCounter=0;nCounter<4;nCounter+)fscanf(fInput2,"%lf",&dOmiganCounter);/*建立高斯方程組矩陣,dSeta0=x1（角度我們草稿的）,dSeta1=x2,dSeta2=x3,dSeta3=x5,dSeta4=x6 dOmiga0 =w2, dOmiga1=w3, dOmiga2=w5, dOmiga3=w6*/a00 = - dLen2 * sin(dSe

22、ta1);a01 = dLen3 * sin(dSeta2);a02 = 0.;a03 = 0.;a10 = dLen2 * cos(dSeta1);a11 = - dLen3 * cos(dSeta2);a12 = 0.;a13 = 0.;a20 = - dLen0 * sin(ALPHA-(dSeta1)+dLen2 * sin(dSeta1);a21 = 0.;a22 = dLen5 * sin(dSeta3);a23 = -dLen6 * sin(dSeta4);a30 = dLen0 * cos(ALPHA-(dSeta1)+dLen2 * cos(dSeta1);a31 = 0.;

23、a32 = dLen5 * cos(dSeta3);a33 = - dLen6 * cos(dSeta4);b0 = - ( - dLen2 * cos(dSeta1) * dOmiga0 * dOmiga0+ dLen3 * cos(dSeta2) * dOmiga1 * dOmiga1 ) + dLen1 *cos(dSeta0) * w1 * w1;b1 = - ( - dLen2 * sin(dSeta1) * dOmiga0 * dOmiga0+ dLen3 * sin(dSeta2) * dOmiga1 * dOmiga1 ) + dLen1 *sin(dSeta0) * w1 *

24、 w1;b2 = -w1 * w1* dLen1*cos(dSeta0)-dOmiga0 *dOmiga0*(dLen2*cos(dSeta1)+dLen0 * cos(ALPHA-(dSeta1)+dOmiga3 * dOmiga3*dLen6 * cos(dSeta4)-dLen5* cos(dSeta3) * dOmiga2 * dOmiga2;b3 = w1 * w1* dLen1*sin(dSeta0)-dOmiga0 *dOmiga0*(dLen2*sin(dSeta1)-dLen0*sin(ALPHA-(dSeta1)-dOmiga3 * dOmiga3*dLen6 * sin(

25、dSeta4)+ dLen5* sin(dSeta3) * dOmiga2 * dOmiga2;if (agaus(a,b,4)!=0)for (i=0;i<=3;i+)printf("a%d=%lfn",i,bi);fprintf(fOutput3,"%lf ",bi);Aex=-dLen2*sin(dSeta1)*b0+dLen0*sin(ALPHA-(dSeta1)*b0-dLen1*cos(dSeta0)*w1*w1-dLen2*cos(dSeta1)*dOmiga0 * dOmiga0-dLen0*cos(ALPHA-(dSeta1)*d

26、Omiga0 * dOmiga0;Aey=dLen2*cos(dSeta1)*b0+dLen0*cos(ALPHA-(dSeta1)*b0-dLen1*sin(dSeta0)*w1*w1-dLen2*sin(dSeta1)*dOmiga0 * dOmiga0+dLen0*sin(ALPHA-(dSeta1)*dOmiga0 * dOmiga0;printf("Aex=%lfnAey=%lf",Aex,Aey);fprintf(fOutput3,"%lf %lf ",Aex,Aey);printf("nn");fprintf(fOutp

27、ut3,"n");fclose(fInput2);fclose(fOutput3);getch();5. #include "stdlib.h" #include "math.h" #include "stdio.h" int dnetn(n,eps,t,h,x,k) int n,k; double eps,t,h,x; extern void dnetnf(); extern int agaus(); int i,j,l; double am,z,beta,d,*y,*a,*b; y=malloc(n*sizeo

28、f(double); a=malloc(n*n*sizeof(double); b=malloc(n*sizeof(double); l=k; am=1.0+eps; while (am>=eps) dnetnf(x,b,n); am=0.0; for (i=0; i<=n-1; i+) z=fabs(bi); if (z>am) am=z; if (am>=eps) l=l-1; if (l=0) free(y); free(b); free(a); printf("failn"); return(0); for (j=0; j<=n-1;

29、j+) z=xj; xj=xj+h; dnetnf(x,y,n); for (i=0; i<=n-1; i+) ai*n+j=yi; xj=z; if (agaus(a,b,n)=0) free(y); free(a); free(b); return(-1); beta=1.0; for (i=0; i<=n-1; i+) beta=beta-bi; if (fabs(beta)+1.0=1.0) free(y); free(a); free(b); printf("failn"); return(-2); d=h/beta; for (i=0; i<=

30、n-1; i+) xi=xi-d*bi; h=t*h; free(y); free(a); free(b); return(k-l); 五計算結果Ø1（單位：度 ）Ø2（單位：度 ）Ø3（單位：度 ）Ø5（單位：度 ）Ø6（單位：度 ）036.7990569.57313122.4627792.665441032.7094766.00299106.4997167.798572029.2776164.0731797.1754452.862793026.5973763.7310390.3834344.249194024.6359164.7551184

31、.2478239.672665023.3010866.8646877.8411637.316076022.4901269.7900470.4814435.419577022.1129773.3007761.2722731.876318022.0990377.2093648.8508123.825399022.3963781.3648232.264598.5942310022.9684485.6444213.79708-12.3279311023.7903589.9462-1.88955-33.31912024.8456494.18335-12.74873-50.8843913026.12331

32、98.28029-19.47863-64.721114027.61533102.17026-23.43679-75.5877215029.31424105.7939-25.76743-84.2467516031.21099109.09881-27.29596-91.2615217033.29288112.03937-28.59909-97.0319318035.54172114.57678-30.07522-101.8606819037.93229116.67878-31.9902-106.0039620040.43114118.31896-34.50611-109.6995121042.99

33、573119.4754-37.70212-113.1728422045.57385120.12878-41.59328-116.625323048.10309120.26007-46.14952-120.2108524050.51035119.84805-51.3162-124.0104525052.71117118.86702-57.03775-128.0136226054.60905117.28503-63.28435-132.1156127056.09499115.06345-70.07877-136.1348628057.0483112.15863-77.51645-139.85181

34、29057.33997108.52777-85.77334-143.0681630056.84147104.14138-95.10476-145.6934431055.4427999.00639-105.86689-147.8969432053.0835193.20271-118.69237-150.5043933049.7968986.93055-135.48887-156.6045834045.7547380.55028-164.06049-178.7216335041.2812574.57477-207.28198-226.9752736036.7990569.57313-237.537

35、23-267.33456 Ø1（單位：度 ）w2（單位：rad/s）w3（單位：rad/s）w5（單位：rad/s）w6（單位：rad/s）0-0.43443-0.43443-2.42454-2.8371910-0.37901-0.27611-1.50231-1.968520-0.30588-0.11091-0.91878-1.2061430-0.230840.0386-0.56502-0.5861540-0.163070.16142-0.32443-0.0690150-0.105640.25595-0.129070.3848760-0.058060.325270.048060.80

36、3970-0.018540.373760.17951.16597800.01490.405460.116711.28627900.043970.4236-0.388620.799141000.070030.43063-1.16737-0.170181100.094080.42829-1.6884-0.946671200.11680.41789-1.86857-1.347151300.138610.40039-1.87297-1.529481400.159680.37661-1.81572-1.626531500.179960.34725-1.74328-1.703731600.199180.3

37、1298-1.66885-1.789361700.216890.27449-1.5923-1.894341800.232460.23246-1.50865-2.020151900.24510.18751-1.41155-2.161372000.253960.14017-1.29477-2.306232100.258090.0908-1.15309-2.4372200.256490.03956-0.98329-2.531192300.24814-0.01365-0.78563-2.564012400.23191-0.06918-0.56583-2.512122500.20665-0.12756-

38、0.33625-2.358052600.17111-0.18948-0.11521-2.093782700.12406-0.255570.0769-1.721942800.06444-0.326130.22384-1.25397290-0.00827-0.400580.32055-0.70581300-0.09331-0.476650.37885-0.09149310-0.18751-0.54910.43120.58577320-0.28388-0.608370.527281.338330-0.37062-0.640050.622262.05413340-0.43232-0.62729-0.6

39、77430.7279350-0.45511-0.55802-3.06708-2.90494360-0.43443-0.43443-2.42454-2.83719Ø1（單位：度 ） a2 a3 a5 a601.87731-1.41416-9.3572522.276910-0.5276-1.3131412.8394314.4391320-10.21122-16.13634-41.38351-111.46712300.322270.3023-2.24307-1.31208400.710490.38945-2.223782.07264502.165745.80322-21.27313-59.

40、0065460-0.144353.8797415.29502-4.92108702.030953.01559-66.66643-92.75689801.34149-0.79682-5.6969317.35997903.873711.2010994.506115.59534100-16.918128.13118-53.8692720.51404110-0.62813-1.177686.97249-2.444621200.977911.4819-1.981126.72186130-0.21204-0.7149-1.26222-4.04646140-1.6445-4.057139.66324.022

41、61150-1.02752-1.887111.3850913.485481600.57409-0.7711137.87656-52.568711700.57409-0.7711137.87656-52.568711800.57409-0.7711137.87656-52.568711900.57409-0.7711137.87656-52.568712000.57409-0.7711137.87656-52.568712100.57409-0.7711137.87656-52.568712200.57409-0.7711137.87656-52.568712300.57409-0.771113

42、7.87656-52.568712400.57409-0.7711137.87656-52.568712500.57409-0.7711137.87656-52.568712600.57409-0.7711137.87656-52.568712700.57409-0.7711137.87656-52.568712800.57409-0.7711137.87656-52.568712900.57409-0.7711137.87656-52.568713000.57409-0.7711137.87656-52.568713100.57409-0.7711137.87656-52.568713200

43、.57409-0.7711137.87656-52.568713300.57409-0.7711137.87656-52.568713400.57409-0.7711137.87656-52.568713500.57409-0.7711137.87656-52.568713600.57409-0.7711137.87656-52.56871 Xe Ye Vex Vey AexAey170.8016537.6482816.3554-36.18842173.19724187.86733172.7165631.863085.73079-29.47339-114.7416822.41059172.88

44、99827.49921-3.42435-20.36538-715.62505-1365.89613171.63824.74528-10.59649-11.372780.5446125.093169.2948423.43323-15.9903-3.997-57.98263103.71014166.1371723.22697-19.990991.28259510.37845288.62216162.3797923.76584-22.902284.5912296.42715358.18851158.1954624.73635-24.904926.2981629.47122281.71398153.7335625.89561-26.09446.82301-0.16224192.37618149.1309227.07215-26.525166.55739569.8676129.35796144.516128.15849-26.241755.84156382.529942655.46328140.0093429.10177-25.296984.96127-44.20371-104.35613135.7202429.89438-23.760814.149548.40659122.9