哈工大材料力学上机实验课作业.docx

《哈工大材料力学上机实验课作业.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、HarbinInstituteofTechnology材料力学电算大作业题目名称：各种条件下的任意截面的弯矩，挠度作者班号：作者学号：作者姓名：刘涛指导教师：王春香完成时间：2013/5/23一、问题的提出该程序可以解决的问题：一个力作用下（1）圆柱面杆（2）矩形截面杆（3）工字梁截面杆（4）空心圆柱截面杆的弯矩和挠度二、程序流程图输入杆长（mm），力F（kN），力作用点到左端面A的距离a（mm）输入梁的弹性模量（Gpa）输入所求端面到左端面的距离X输入梁的形状计算惯性矩计算并输出弯矩根据梁的形状不同输入不同的数值计算并输出挠度三、

2、不同力作用下的弯矩，挠度1.一个力作用下的任意截面的弯矩，挠度程序：#include#include#include#definePI3.intget_I(intkind);intmain(){floatf=100,l=100,a=50,b,x=60;floatM,v;floatE=200,I;intkind=1;printf("请输入杆长l(mm)，力F（kN），力作用点到左端面A的距离a(mm)");scanf("%f%f%f",&l,&f,&a);l/=1000;f*=

3、1000;a/=1000;printf("请输入梁的弹性模量(GPa)");scanf("%f",&E);E*=1.e9;printf("请输入所求端面到左端面A的距离X");scanf("%f",&x);x/=1000;printf("请输入梁的截面形状（圆柱面=1,矩形梁=2,工字梁=3,空心圆柱梁=4）");scanf("%d",&kind);I=get_I(kind);//计算惯性矩b=l-a;M=f*x;//计算弯矩if(0<=x&&x<=a){v=(-1)*f*b*x*(l*l-x*x-b*b)/(6*E*I

4、*l);}elseif(a

5、oatI;//圆柱面=1,矩形梁=2,工字梁=3,空心圆柱梁=4if(1==kind)//圆柱面=1{printf("请输入圆柱梁的直径d");scanf("%f",&d);I=PI*pow(d,4)/64;}elseif(2==kind)//矩形梁=2{printf("请输入矩形截面的b和h");scanf("%f%f",&b,&h);I=b*pow(h,3)/12;}///*elseif(3==kind)//工字梁=3{printf("

6、——————-—————

7、");printf("

8、————

9、

10、——————

11、t

12、n");printf("

13、d

14、");printf("

15、

16、");printf("

17、————

18、

19、——————

20、");printf("

21、—————b——————

22、");/*printf("请输入工字梁的翼缘厚度t，腹板厚度d，翼缘宽度b，梁高度h");scanf("%f%f%f%f",&t,&d,&b,&h);*/printf("请查附表B输入工字梁的惯性矩I");scanf("%f",&I);}//*/elseif(4==kind)//空心圆柱梁=4{printf("请输入圆柱梁的直径d和内外径比例alpha

23、");scanf("%f%f",&d2,&alpha);I=PI*pow(d2,4)*(1-pow(alpha,4))/64;}else{I=0;}returnI;}2.一个力偶作用下的任意截面的弯矩，挠度#include#include#include#definePI3.#defineTORAD(PI/180)floatget_I(intkind);floatget_V_M0(floatE,floatLength,floatI,floatx,floatM0,floata_L)

24、;floatget_M_M0(floatLength,floatM0,floata_L,floatx);intmain(){intkind;//截面形状floatLength,a_L,x;floatE,I,M0;floatv,M;//