材料力学验证编程作业.docVIP

课次 作业 分值 1 一个力作用下的任意截面的弯矩，挠度。 １ 2 1.一个力偶作用下的任意截面的弯矩，挠度。2.一段分布力作用下的任意截面的 弯矩，挠度。 ２ 3 n个力共同作用下的任意截面的弯矩，挠度。 ２ 4 1.h个力偶共同作用作用下的任意截 面的弯矩，挠度。 2.m段分布力共同作用下的任意截 面的弯矩，挠度。 ２ 5 n个力，h个力偶，m段分布力共同作用下的任意截面的弯矩，挠度。 ２ 6 n个力，h个力偶，m段分布力共同作用下的最大弯矩，最大挠度。 １ 以下7题的材料为20b工字型钢，其杨氏模量E=200GPa，惯性矩Iz=0.000025㎡×㎡ 程序1： #include stdio.h #include math.h main() { double Fi,ai,bi,l,Mx,vx,x,E,I; printf(ShuRu Fi,ai,bi,E,I,x:\n); scanf(%lf%lf%lf%lf%lf%lf,Fi,ai,bi,E,I,x); l=ai+bi; if (xl||x0) printf(Error!); if (0=xx=ai) { Mx=bi*Fi*x/l; vx=-Fi*bi*x*(l*l-x*x-bi*bi)/(6*E*I*l); printf(Mx=%lf,vx=%e\n,Mx,vx); } if (aixx=l) { Mx=ai*Fi*(l-x)/l; vx=-Fi*bi*x*(l*(x-ai)*(x-ai)*(x-ai)/bi+(l*l-bi*bi)*x-x*x*x)/ (6*E*I*l); printf(Mx=%lf,vx=%e\n,Mx,vx); } getch(); } 例题：已知：F=3kN,a=0.3m,b=0.8m; 求：x=0.15,x=0.3m,x=0.7m时的弯矩与挠度。 x=0.15m x=0.3m x=0.7m 2-1）程序2-1: #include stdio.h #include math.h main() { double Mk,ak,bk,l,Mx,vx,x,E,I; printf(ShuRu Mk,ak,bk,E,I,x:\n); scanf(%lf%lf%lf%lf%lf%lf,Mk,ak,bk,E,I,x); l=ak+bk; if (xl||x0) printf(Error!); if (0=xx=ak) { Mx=-Mk*x/l; vx=Mk*x*(l*l-x*x-3*bk*bk)/(6*E*I*l); printf(Mx=%lf,vx=%e\n,Mx,vx); } if (akxx=l) { Mx=Mk*(l-x)/l; vx=Mk*(3*l*(x-ak)*(x-ak)+(l*l-3*bk*bk)*x-x*x*x)/(6*E*I*l); printf(Mx=%lf,vx=%e\n,Mx,vx); } getch(); } 例题：已知：mk=300Nm，ak=0.5m,bk=0.9m 求：x=0.25m,x=0.5m,x=1m时的弯矩与挠度。 x=0.25m x=0.5m x=1m 2-2）程序2-2 #include stdio.h #include math.h main() { double q,aj,bj,cj,l,Mx,vx,x,E,I,Fa,Fb,C1,C2,C3,D3,A; printf(ShuRu q,aj,bj,cj,E,I,x:\n); scanf(%lf%lf%lf%lf%lf%lf%lf,q,aj,bj,cj,E,I,x); l=aj+bj+cj; Fa=q*cj*(2*bj+cj)/(2*l); Fb=q*cj*