《材力上機Matlab哈工大》由會員分享，可在線閱讀，更多相關《材力上機Matlab哈工大(7頁珍藏版)》請在金鋤頭文庫上搜索。

1、一、解決題目由于之前只學過簡單的 C語言， 對于本次程序設計的幫助不大， 所以重新學習了Matlab 語言，并對題目 7 進行處理題目 7 應力狀態分析輸入：1. 平面應力狀態輸入：( x、 y、 xy)；某截面方位角 2. 空間應力狀態輸入：( x、 y、 z、 xy、 yz、 zx)輸出：1. 輸出主應力( 1、 2、 3)2. 最大剪應力( max= 13=23-1 錯誤！未找到引用源。 )3. 如為平面應力狀態則需輸出方位角 斜截面上的應力 、 及主方向角 錯誤！未找到引用源。4. 畫出應力圓示意圖二、方向規定正面上應力以沿著坐標軸正方向者為正， 反之為負： 負面上應力以沿著坐標軸負方。

2、向者為正，反之為負。三、程序設計思路本程序分為兩個部分， 分別對應平面應力分析與空間應力分析， 可由使用者進行選擇。平面應力分析主要是主應力和任一截面位置上的應力狀態，還有應力圓示意圖。空間應力狀態主要是主應力問題四、程序編寫程序采用 Matlab 語言，所以對于一些特殊字符，無法識別。首先對變量進行定義：n 為用戶選擇進入的兩個程序部分， n=1 時，進行平面應力狀態分析， n=2 時，進行空間應力狀態分析。x， y， z 分別對應三個方向的正應力。xy， yz， zx 分別對應切應力( yx， zy， xz 省略) 。j 為用戶輸入的某截面 角其他一些代碼為程序內部執行代碼，不再一一介紹。。

3、下面是主程序代碼n=input(If plane,enter 1,If spatial,enter 2 n=); 用戶選擇并進入相應界面if n=1 平面應力分析x=input(x=); y=input(y=); xy=input(xy=); j=input(j=); 用戶輸入數據jj=j*pi/180; aa=(x+y)/2; bb=sqrt(x-y)/2)2+xy2); cc=(x-y)/2; a=aa+bb; b=aa-bb; if a=b a1=a a2=0 a3=b tmax=(a1-a3)/2 az=aa+cc*cos(2*jj)+xy*sin(2*jj) tz=-cc*sin(2。

4、*jj)+xy*cos(2*jj) alpha=0:pi/50:2*pi; X=sqrt(b2+xy2)*cos(alpha)+a; Y=sqrt(b2+xy2)*sin(alpha); plot(X,Y,-) else a1=b a2=0 a3=a tmax=(a1-a3)/2 az=aa+cc*cos(2*jj)+xy*sin(2*jj) tz=-cc*sin(2*jj)+xy*cos(2*jj) alpha=0:pi/50:2*pi; X=sqrt(b2+xy2)*cos(alpha)+a; Y=sqrt(b2+xy2)*sin(alpha); plot(X,Y,-) end 計算，繪圖。

5、并結束程序else x=input(x=); y=input(y=); z=input(z=); xy=input(xy=); yz=input(yz=); zx=input(zx=); I1=x+y+z; I2=-xy*xy+x*y+y*z-yz*y+z*x-zx*zx; I3=-zx*y*zx-x*yz*yz-z*xy*xy+x*y*z+xy*yz*zx+zx*xy*yz; syms Z; solve(Z3-I1*Z2+I2*Z-I3) b1=eval(ans(3); b2=eval(ans(2); b3=eval(ans(1); b0=0; if b10 b0=b1; b1=b2; b2。

6、=b0; if b10 b0=b1; b1=b3; b3=b0; b0=b2; b3=b2; b3=b0; elseif b20 b0=b2; b2=b3; b3=b0; end else if b20 b0=b2; b2=b3; b3=b0; b3=b2; b2=b0; end end a1=b1 a2=b2 a3=b3 tmax=0.5*(a1-a3) t23=0.5*(a2-a3) t12=0.5*(a1-a2) alpha1=0:pi/50:2*pi; alpha2=0:pi/50:2*pi; alpha3=0:pi/50:2*pi; X1=tmax*cos(alpha1)+(a1+a。

7、3)/2; Y1=tmax*sin(alpha1); X2=t23*cos(alpha2)+(a2+a3)/2; Y2=t23*sin(alpha2); X3=t12*cos(alpha3)+(a1+a2)/2; Y3=t12*sin(alpha3); hold on plot(X1,Y1,k); plot(X2,Y2,k); plot(X3,Y3,k); end 五、實例展示下面我用書中的例題向大家展示我的程序例題 1: 材料力學教材 p33 2-7 ( a)平面應力已知 x=50, y=0, xy=20, =50 。選擇平面應力分析 x=50 y=0 xy=20 =50 圖片一輸出結果圖片。

8、二 圖片三應力圓輸出圖片四例題 2：材料力學教材 p35 2-14 ( a)已知 x=60； y=30； z=-70；圖片五圖片六 圖片七圖片八六、程序缺陷1. 沒有生成標準的 EXE文件，不能很好地應用2. 用戶輸入窗口不夠人性化3. 應力圓輸出并非圓形七、課程收獲在這次上機編程作業中， 雖然花費了很長時間， 但是我感到收獲很多。 在此過程中， 我學習了 Matlab 語言并對 Matlab 語言有了更加熟悉的運用， 這將成為我一生的財富， Matlab 語言也必將成為我日后工作的得力助手。在看到自己做出的最終結果以及它運行時的高效， 我不禁有了一絲成就感， 我為能有一次這樣鍛煉的機會而感到幸運， 希望以后還會更多這樣的機會充實自己。 對于材料力學這門課程， 我想也會因為這次程序的編寫而記憶的更加深刻， 我一定要爭取多掌握一些這方面的知識，然后在工作中充分的應用進去！。

總結

以上是生活随笔為你收集整理的哈工大材力上机 matlab,材力上机Matlab哈工大的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。