材料力学编程.docxVIP

下载本文档

2
0
约1.77万字
约 21页
2019-08-11 发布于湖北
举报
版权申诉

材料力学编程.docx

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

材料力学上机作业第一题一个力作用下的任意截面的弯矩，挠度源程序为： #includestdio.h #includemath.h main() { 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) (x=ai)) { M=Fi*bi*x/l; v=-Fi*bi*x*(l*l-x*x-bi*bi)/(6*E*I*l); printf(M=%lf,v=%lf\n,M,v); } if((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=%lf\n,M,v); } if((x0) || (xl)) { printf(有错误!\n); } } Fi=300N，ai=60m，bi=200m，E=200000000000Pa，I=0.0001m4，x=80m，程序运行结果为： v 材料力学上机作业第二题 1一个力偶作用下的任意截面的弯矩，挠度 #include stdio.h #include math.h 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) (x=ak)) { M=-Mk*x/l; v=Mk*x*(l*l-3*bk*bk-x*x)/(6*E*I*l); printf(M=%lf,v=%lf\n,M,v); } if((xak) (x=l)) { M=Mk*(1-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=%lf\n,M,v); } if((x0) || (xl)) { printf(错误!\n); } } 令Mk=200N·m，ak=100m，bk=100m，E=200000000000Pa，I=0.0001m4，x=100m，程序运行结果为： 2一段分布力作用下的任意截面的弯矩，挠度源程序为： #include stdio.h #include math.h main() { double q,ai,bi,ci,l,E,I,x,M,v; double Fa,Fb,c1,c2,c3,d3; printf(请输入： q,ai,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=c2; d3=-c3*l; if((x=0) (x=ai)) { M=Fa*x; v=((Fa*pow(x,3))/6+c1*x)/(E*I); printf(M=%lf,v=%lf\n,M,v); } if((xai) (x=ai+ci)) { M=((q*ci*x*(ci+2*bi))/(2*l))-((q*pow((x-ai),2))/2); v=((Fa*pow(x,3))/6-(q*pow((x-ai),4))/24+c2*x)/(E*I); printf(M=%lf,v=%lf\n,M,v); } if((xai+ci) (x=l)) { M=Fb*(l-x); v=((Fb*pow((l-x),3))/6+c3*x+d3)/(E*I); printf(M=%lf,v=%lf\n,M,v