《机构设计冲压模具的拉延机构源程序》由会员分享，可在线阅读，更多相关《机构设计冲压模具的拉延机构源程序（6页珍藏版）》请在金锄头文库上搜索。
1、冲压模具的拉延机构源程序#include #include #define P 3.1414926 #define E1 1e-3 #define E0 1e-9 #define m 4 #define N0 1000 double F(double xm,double t,int n) double pm; p0=440*cos(t*P/180)+490*cos(x0)+249*cos(x1-20*P/180)-x3;p1=440*sin(t*P/180)+490*sin(x0)-249*sin(x1-20*P/180);p2=500*cos(t*P/180+20*P/180)+370*co
2、s(x2)+540.2285*cos(x1)-x3;p3=500*sin(t*P/180+20*P/180)+370*sin(x2)-540.2285*sin(x1); return pn; double A(double xm,int i,int j) double pmm; p00=-490*sin(x0);p01=-249*sin(x1-20*P/180);p02=0;p03=-1; p10=490*cos(x0);p11=-249*cos(x1-20*P/180);p12=0;p13=0; p20=0;p21=-540.2285*sin(x1);p22=-370*sin(x2);p23
3、=-1; p30=0;p31=-540.2285*cos(x1);p32=370*cos(x2);p33=0; return pij; double G(double amm,int n,double bm) int i,j,k,i1,j1; double temp; int row_maxmod; double element_maxmod; for(j=0;j=0;i-)for(j=3;ji;j-)bi-=aij*bj;return bn;double C(double xm,double vm,double t,int n) double pm;p0=440*2*P*2*P*cos(t*
4、P/180)+490*v0*v0*cos(x0)+249*v1*v1*cos(x1- 20*P/180);p1=440*2*P*2*P*sin(t*P/180)+490*v0*v0*sin(x0)- 249*v1*v1*sin(x1-20*P/180);p2=500*2*P*2*P*cos(t*P/180+20*P/180)+370*(v2)*(v2)*cos(x2)+540.2285*( v1)*(v1)*cos(x1);p3=500*2*P*2*P*sin(t*P/180+20*P/180)+370*(v2)*(v2)*sin(x2)- 540.2285*(v1)*(v1)*sin(x1)
5、;return pn;int main()int i1=0,i,j;double x4,a44,b4;double s4,v4,p4,c4;double t;printf(“求解：t t2(度） t4(度） t5(度） XD(滑块位移)n“);x0=13.02531*P/180;x1=46.32880*P/180;x2=36.43463*P/180;x3=1140.56222;for(t=0.0;tE1|fabs(b1)E1|fabs(b2)E1|fabs(b3)E1)elseprintf(“位移精确值为： “);for(i=0;i3;i+)printf(“%-15.5f“,xi*(180/P
6、);for(i=3;i4;i+)printf(“%-15.5f“,xi);printf(“n“);for(i=0;i4;i+) for(j=0;j4;j+) aij=A(x,i,j); b0=440*2*P*sin(t*P/180); b1=-440*2*P*cos(t*P/180);b2=500*2*P*sin(t*P/180+20*P/180); b3=-500*2*P*cos(t*P/180+20*P/180); for(i=0;i4;i+) vi=G(a,i,b); printf(“速度精确值为： “); for(j=0;j4;j+) printf(“%-15.5f“,vj); printf(“n“); for(j=0;j4;j+) cj=C(x,v,t,j); printf(“加速度精确值： “); for(i=0;i4;i+) pi=G(a,i,c);for(j=0;j4;j+) printf(“%15.5f“,pj); printf(“n“);return 0;结果如图：机构设计院系：机电工程学院班级：机设 10-03 班姓名：邢剑锋学号：541002010347
