哈工大-材料力学上机实验-悬臂梁弯矩挠度c语言完美整合详解.docVIP

Harbin Institute of Technology 课程设计 课程名称： 材料力学 设计题目： 悬臂梁的弯矩及挠度 院 系： 机电工程学院 班 级： 设 计 者： 学 号： 设计时间： 指导老师： 哈尔滨工业大学 一 问题描述#include #include double Mz,l,a,x,v,E,I,b,F; main printf 请输入长度l\n ; scanf %lf,l ; printf 请输入长度a\n ; scanf %lf,a ; printf 请输入弹性模量E\n ; scanf %lf,E ; printf 请输入极惯性矩I\n ; scanf %lf,I ; printf 请输入力F\n ; scanf %lf,F ; printf 请输入x\n ; scanf %lf,x ; b l-a; if x 0x a Mz F* l-a *x/l; v -F*b*x* l*l-x*x-b*b / 6*E*I*l ; printf Mz %f\n,Mz ; printf v %f\n,v ; else if x ax l Mz F*a* l-x /l; v -F*b* l/b*pow x-a,3 + l*l-b*b *x-x*x*x / 6*E*I*l ; printf Mz %f\n,Mz ; printf v %f\n,v ; else printf 數值溢出\n ; 【2】检验数据： 条件：L 1m，E 4.5GPa。I 0.0001m4，F 1000N，a 0.6m，x 0.3m； 理论运算结果： Mz 120N m V - 0.000033m 【3】程序运行结果： 图二：程序一运行 结论：运行结果与理论运算结果相同，可以使用。 （3）一个力偶作用下的任意截面的弯矩，挠度。 【1】源代码： #include #include double Mz,l,a,x,v,E,I,b,M; main printf 请输入长度l\n ; scanf %lf,l ; printf 请输入长度a\n ; scanf %lf,a ; printf 请输入弹性模量E\n ; scanf %lf,E ; printf 请输入极惯性矩I\n ; scanf %lf,I ; printf 请输入弯矩M\n ; scanf %lf,M ; printf 请输入x\n ; scanf %lf,x ; b l-a; if x 0x a Mz M*x/l; v M*x* l*l-3*b*b-x*x / 6*E*I*l ; printf Mz %f\n,Mz ; printf v %f\n,v ; else if x ax l Mz M* l-x /l; v M* -x*x*x+3*l* x-a * x-a + l*l-3*b*b *x ; printf Mz %f\n,Mz ; printf v %f\n,v ; else printf 數值溢出\n ; 【2】检验数据： 条件：L 1m，E 4.5GPa。I 0.0001m4，M 1000Nm，a 0.6m，x 0.3m； 理论计算结果：Mz 300.00Nm V 0.000048m 【3】程序运行结果： 图三：程序一运行 结论：运行结果与理论运算结果相同，可以使用。 （4）一段分布力作用下的任意截面的弯矩挠度： 【1】源代码： #include #include double q,l,c,d,x,E,I,v,Mz; double C,D,Ec,Fc; double FR1,FR2; /*函数功能：奇异函数*/ double qiyi double x,double y if x y return 0; if x y return x-y; main printf 请输入长度l\n ; scanf %lf,l ; printf 请输入弹性模量E\n ; scanf %lf,E ; printf 请输入极惯性矩I\n ; scanf %lf,I ; printf 请输入分布力q\n ; scanf %lf,q ; printf 请输入长度c\n ; scanf %lf,c ; printf 请输入长度d\n ; scanf %lf,d ; printf 请输入x\n ; scanf %lf,x ; C qiyi x,c ; D qiyi x,d ; FR2 -q* d*d-c*c /2 /l;//求支反力 FR1 -q* d-c -FR2; Ec - FR1*pow l,3 /6-q*pow l-c,4 /24+q*pow l-d,4 /24 /E/I/l; F