《《材料力学编程》》由会员分享，可在线阅读，更多相关《《材料力学编程》（29页珍藏版）》请在金锄头文库上搜索。
1、材料力学上机作业 第一题一个力作用下的任意截面的弯矩，挠度源程序为：#include#includemain()double Fi,ai,bi,l,E,I,x,M,v;printf(请输入： Fi,ai,bi,E,I,x:n);scanf(%lf,%lf,%lf,%lf,%lf,%lf,&Fi,&ai,&bi,&E,&I,&x);l=ai+bi;if(x=0) & (xai) & (x=l)M=Fi*ai*(l-x)/l;v=-Fi*bi*(l*(x-ai)*(x-ai)*(x-ai)/bi+(l*l-bi*bi)*x-x*x*x)/(6*E*I*l);printf(计算结果为：M=%lf,v
2、=%lfn,M,v); if(xl)printf(有错误!n);Fi=300N，ai=60m，bi=200m，E=200000000000Pa，I=0.0001m4，x=80m，程序运行结果为：v材料力学上机作业 第二题1一个力偶作用下的任意截面的弯矩，挠度#include #include main()double Mk,ak,bk,l,E,I,x,M,v;printf(请输入 Mk,ak,bk,E,I,x:n);scanf(%lf,%lf,%lf,%lf,%lf,%lf,&Mk,&ak,&bk,&E,&I,&x);l=ak+bk;if(x=0) & (xak) & (x=l)M=Mk*(1
3、-x/l);v=-Mk*(-x*x*x+3*l*(x-ak)*(x-ak)+(l*l-3*bk*bk)*x)/(6*E*I*l);printf(M=%lf,v=%lfn,M,v);if(xl)printf(错误!n);令Mk=200Nm，ak=100m，bk=100m，E=200000000000Pa，I=0.0001m4，x=100m，程序运行结果为：2一段分布力作用下的任意截面的弯矩，挠度源程序为：#include #include main()double q,ai,bi,ci,l,E,I,x,M,v; double Fa,Fb,c1,c2,c3,d3;printf(请输入： q,ai,
4、bi,ci,E,I,x:n);scanf(%lf,%lf,%lf,%lf,%lf,%lf,%lf,&q,&ai,&bi,&ci,&E,&I,&x);l=ai+bi+ci; Fa=(q*ci*(ci+2*bi)/(2*l);Fb=(q*ci*(ci+2*ai)/(2*l);c3=(Fa*pow(ai+ci),3)/(3*l)-(q*ai*pow(ci,3)/(6*l)-(q*pow(ci,4)/(8*l)+(Fb*pow(bi,2)/2-(Fb*pow(bi,3)/(3*l);c2=c3-(Fb*pow(bi,2)/2-(Fa*pow(ai+ci),2)/2+(q*pow(ci,3)/6;c1=
5、c2;d3=-c3*l;if(x=0) & (xai) & (xai+ci) & (x=l) M=Fb*(l-x); v=(Fb*pow(l-x),3)/6+c3*x+d3)/(E*I); printf(M=%lf,v=%lfn,M,v);if(xl)printf(有错误!n);q=200N/m，ai=100m，bi=200，ci=100m，E=200000000000Pa，I=0.0001m4，x=70m，程序运行结果为：材料力学上机作业 第三题n个力共同作用下的任意截面的弯矩，挠度源程序为：#include #include #define EI 4.9e5#define N 1000do
6、uble FunctionV(int i,int n,double l,double x,double a,double b,double F);double FunctionM(int i,int n,double l,double x,double a,double b,double F);main()double l,x;int i,n;double vN,MN,aN,bN,FN;double sumv=0,sumM=0;printf(请输入 n:);scanf(%d,&n);printf(请输入 l/m:);scanf(%lf,&l);printf(请输入 x/m:);scanf(%l
7、f,&x);for(i=0;i=0)&(x=0)&(x=ai)M=(Fi*bi*x)/l;elseM=(Fi*ai*(l-x)/l;return M;令n=2，l=200m，x=100m，a0=50m，F0=1000N，a1=100m，F1=1000N，程序运行结果为：材料力学上机作业 第四题1、h个力偶共同作用下的任意截面的弯矩，挠度源程序为：#include #include #define EI 4.9e5#define N 1000double FunctionV(int i,int n,double l,double x,double a,double b,double me);do
8、uble FunctionM(int i,int n,double l,double x,double a,double b,double me);main()double l,x;int i,h;double vN,MN,aN,bN,meN;double sumv=0,sumM=0;printf(请输入h:);scanf(%d,&h);printf(请输入l/m:);scanf(%lf,&l);printf(请输入x/m:);scanf(%lf,&x);for(i=0;i=0)&(x=ai)v=(mei*x*(pow(l,2)-3*pow(bi,2)-pow(x,2)/(6*EI*l);elsev=(mei*(-pow(x,3)+3*l*pow(x-ai),2)+(pow(l,2)-3*pow(bi,2)*x)/(6*EI*l);return v;double FunctionM(int i,int
