资源描述：
《材料力学的fortran编程》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、为了适应公司新战略的发展，保障停车场安保新项目的正常、顺利开展，特制定安保从业人员的业务技能及个人素质的培训计划材料力学的fortran编程 材料力学上机实验报告 课程：材料力学 题目：应力状态分析 班级： 学号： 姓名：姜超 题目7:应力状态分析 输入: 1.平面应力状态输入： 2.空间应力状态输入: 2.最大切应力[τmax=τ13xy、τyz、τzx）=/2] α3.如为平面应力状态则需输入方位角α斜截面上的应力σα、τ 及主方向角α0。 4.画出应力圆示意图。 材料力学上机实验Fort
2、ran编写： Programmain Implicitnone Integerk目的-通过该培训员工可对保安行业有初步了解，并感受到安保行业的发展的巨大潜力，可提升其的专业水平，并确保其在这个行业的安全感。为了适应公司新战略的发展，保障停车场安保新项目的正常、顺利开展，特制定安保从业人员的业务技能及个人素质的培训计划 Reala1,a2,a3,b1,b2,b3,bmax,p,t1,t2,t3,a,b,p0RealTemp1,Temp2,Temp3 Print*,“1=平面应力状态，2=空间应力状态”Print*,
3、“请选择1或2” Read*,k If(k=1)then Print*,“请输入平面应力状态参数：?x,?y,?xy,?” Read*,a1,a2,b1,p Temp1=(a1+a2)/2+sqrt(((a1-a2)/2)**2+b1**2)Temp2=0 Temp3=(a1+a2)/2+sqrt(((a1-a2)/2)**2+b1**2)t1=max(Temp1,Temp2,Temp3) t3=min(Temp1,Temp2,Temp3) t2=(Temp1+Temp2+Temp3)-(t1+t3) b
4、max=(t1-t3)/2 a=(a1+a2)/2+(a1-a2)/2*cos(2*p)-b1*sin(2*p)b=(a1-a2)/2*sin(2*p)+b1cos(2*p) p0=1/2*atan(-2*b1/(a1-a2)) Print*,“截面主应力为：”，“?1=”,t1,“?2=”,t2,“?3=”,t3Print*,“最大切应力为：?max=”，bmax目的-通过该培训员工可对保安行业有初步了解，并感受到安保行业的发展的巨大潜力，可提升其的专业水平，并确保其在这个行业的安全感。为了适应公司新战略的发展，保
5、障停车场安保新项目的正常、顺利开展，特制定安保从业人员的业务技能及个人素质的培训计划 Print*,“任意?角截面上的应力和主方向角为：”，“??=”， “??=”，b，“?0=”，p0 Endif If(k=2)then Print*,“请输入空间应力状态参数：?x,?y,?z,?xy,?yz,?zx”Read*,a1,a2,a3,b1,b2,b3 If(b1==0,b2==0)thena， Temp1=a2 Temp2=(a1+a3)/2+sqrt((a1-a3)**2+b3**2)Temp3=(a1+
6、a3)/2-sqrt((a1-a3)**2+b3**2)t1=max(Temp1,Temp2,Temp3) t3=min(Temp1,Temp2,Temp3) t2=(Temp1+Temp2+Temp3)-(t1+t3) bmax=(t1-t3)/2 elseif(b2==0,b3==0)then Temp1=a3 Temp2=(a1+a2)/2+sqrt((a1-a2)**2+b1**2)Temp3=(a1+a2)/2-sqrt((a1-a2)**2+b1**2)t1=max(Temp1,Temp2,Temp
7、3)t3=min(Temp1,Temp2,Temp3) t2=(Temp1+Temp2+Temp3)-(t1+t3)目的-通过该培训员工可对保安行业有初步了解，并感受到安保行业的发展的巨大潜力，可提升其的专业水平，并确保其在这个行业的安全感。为了适应公司新战略的发展，保障停车场安保新项目的正常、顺利开展，特制定安保从业人员的业务技能及个人素质的培训计划 bmax=(t1-t3)/2 elseif(b3==0,b1==0)then Temp1=a1 Temp2=(a2+a3)/2+sqrt((a2-a3)**2+b
8、2**2)Temp3=(a2+a3)/2-sqrt((a2-a3)**2+b2**2)t1=max(Temp1,Temp2,Temp3) t3=min(Temp1,Temp2,Temp3) t2=(Temp1+Temp2+Temp3)-(t1+t3) bmax=(t1-t3)/2 else p
