执行机构程序.docVIP

六杆机构分析源程序 目录： void main() 1 void Solutionangle(double value[M],double ang1) 4 void Solutionspeed(double a2[N][N],double b2[N],double value[M],double ang1) 5 void Solutionacceleration(double a3[N][N],double da3[N][N],double db3[N],double value[M]) 5 void GaussianE(double a4[N][N],double b4[N],double p4[N]) 5 void FoundmatrixA(double value5[M],double a5[N][N]) 6 void FoundmatrixB(double value6[M],double ang1,double b6[N]) 7 void FoundmatrixDA(double value7[M],double da7[N][N]) 7 void FoundmatrixDB(double value8[M],double ang1,double db8[N]) 7 void mf(double value9[M],double ang1) 7 程序： #includestdio.h #includestdlib.h #includemath.h #define PI 3.1415926 #define N 4 #define M 15 void Solutionangle(double [M],double ); /*迭代法求角位移*/ void Solutionspeed(double [N][N],double [N],double [M],double ); /*角速度求解*/ void Solutionacceleration(double [N][N],double [N][N],double [N],double [M]);/*角加速度求解*/ void GaussianE(double [N][N],double [N],double [N]);/*高斯消去*/ void FoundmatrixA(double [M],double [N][N]); //创建系数矩阵A void FoundmatrixB(double [M],double ,double [N]);//创建系数矩阵B void FoundmatrixDA(double [M],double [N][N]);//创建矩阵DA void FoundmatrixDB(double [M],double ,double [N]);//创建矩阵DB void mf(double [M],double ); //驱动力矩、阻力矩求解 //定义全局变量 double as1=2*PI,h=0.120,K=2,lo23=0.150,az=0.055,bz=0.055,cz=0.125,l1=0.120,l3=0.160,l4=0.160,x=0.160,h1=0.136; double q=2200,g3=140,g5=340,sbju=0.03,g=9.8,js3=0.11,zm=0,pm=0,zdg=0,djg=0,jf=0; //主函数 void main() { FILE *fp; int fangan; for(fangan=0;fangan4;fangan++) { switch(fangan) { case 0:as1=PI;h=0.150;K=2;lo23=0.160;az=0.060;bz=0.060;cz=0.130; q=2300;g3=150;g5=350;js3=0.12; if((fp=fopen(A.xls,w))==NULL) { printf(Cannt open this file.\n); exit(0); }break; case 1:as1=2*PI;h=0.120;K=2;lo23=0.150;az=0.055;bz=0.055;cz=0.125; q=2200;g3=140;g5=340;js3=0.11; if((fp=fopen(B.xls,w))==NULL) { printf(Cannt open this file.\n); exit(0); }break; case 2:as1=3*PI;h=0.090;K=2;lo23=0.140;az=0.050