材料力學的fortran編程

《材料力學的fortran編程》由會員上傳分享，免費在線閱讀，更多相關內容在應用文檔-天天文庫。

1、為了適應公司新戰(zhàn)略的發(fā)展，保障停車場安保新項目的正常、順利開展，特制定安保從業(yè)人員的業(yè)務技能及個人素質的培訓計劃材料力學的fortran編程　　材料力學上機實驗報告　　課程：材料力學　　題目：應力狀態(tài)分析　　班級：　　學號：　　姓名：姜超　　題目7:應力狀態(tài)分析　　輸入:　　1.平面應力狀態(tài)輸入：　　2.空間應力狀態(tài)輸入:　　2.最大切應力[τmax=τ13xy、τyz、τzx）=/2]　　α3.如為平面應力狀態(tài)則需輸入方位角α斜截面上的應力σα、τ　　及主方向角α0?！　?.畫出應力圓示意圖?！　〔牧狭W上機實驗Fort

2、ran編寫：　　Programmain　　Implicitnone　　Integerk目的-通過該培訓員工可對保安行業(yè)有初步了解，并感受到安保行業(yè)的發(fā)展的巨大潛力，可提升其的專業(yè)水平，并確保其在這個行業(yè)的安全感。為了適應公司新戰(zhàn)略的發(fā)展，保障停車場安保新項目的正常、順利開展，特制定安保從業(yè)人員的業(yè)務技能及個人素質的培訓計劃　　Reala1,a2,a3,b1,b2,b3,bmax,p,t1,t2,t3,a,b,p0RealTemp1,Temp2,Temp3　　Print*,“1=平面應力狀態(tài)，2=空間應力狀態(tài)”Print*,

3、“請選擇1或2”　　Read*,k　　If(k=1)then　　Print*,“請輸入平面應力狀態(tài)參數(shù)：?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目的-通過該培訓員工可對保安行業(yè)有初步了解，并感受到安保行業(yè)的發(fā)展的巨大潛力，可提升其的專業(yè)水平，并確保其在這個行業(yè)的安全感。為了適應公司新戰(zhàn)略的發(fā)展，保

5、障停車場安保新項目的正常、順利開展，特制定安保從業(yè)人員的業(yè)務技能及個人素質的培訓計劃　　Print*,“任意?角截面上的應力和主方向角為：”，“??=”，　　“??=”，b，“?0=”，p0　　Endif　　If(k=2)then　　Print*,“請輸入空間應力狀態(tài)參數(shù)：?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)目的-通過該培訓員工可對保安行業(yè)有初步了解，并感受到安保行業(yè)的發(fā)展的巨大潛力，可提升其的專業(yè)水平，并確保其在這個行業(yè)的安全感。為了適應公司新戰(zhàn)略的發(fā)展，保障停車場安保新項目的正常、順利開展，特制定安保從業(yè)人員的業(yè)務技能及個人素質的培訓計劃　　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