如何用MATLAB實(shí)現(xiàn)根據(jù)主應(yīng)力生成相應(yīng)的應(yīng)力張量,并根據(jù)張量進(jìn)行相應(yīng)的計算和轉(zhuǎn)換
-------------------------------------討論-------------------------------------
1、根據(jù)彈塑性力學(xué),可將巖體內(nèi)一個點(diǎn)的應(yīng)力狀態(tài)用一個2階應(yīng)力張量[σij]表示;
2、應(yīng)力張量[σij]展開后可得到應(yīng)力方程σ^3+i1*σ^2+i2*σ+i3=0,
i1、i2、i3為應(yīng)力張量的3個不變量;
3、應(yīng)力狀態(tài)[σij]可分解為兩個部分,即[σij]=[σm*δij]+[sij],
[σm]=([σ1]+[σ2]+[σ3])/3,[sij]為[σij]在pi平面上的投影,稱為偏應(yīng)力張量;
4、偏應(yīng)力張量[sij]展開后同樣可得偏應(yīng)力方程s^3+j1*s^2+j2*s+j3=0,
j1、j2、j3為偏應(yīng)力的3個不變量。
5、還有一些關(guān)于屈服面參數(shù)的計算都與不變量的值有關(guān)。
6、根據(jù)以上分析,我認(rèn)為可以通過matlab,根據(jù)三個主應(yīng)力的值得到上面分析的每一個步驟的重要參數(shù),如應(yīng)力張量和偏應(yīng)力張量的值。
-------------------------------------提問-------------------------------------
因此小弟想知道
1.如何通過matlab,將已知的三個主應(yīng)力轉(zhuǎn)換成應(yīng)力張量[σij]的形式,并得到相應(yīng)的不變量和[sij]及相應(yīng)的不變量;
2.如何把[[σij]投影到pi平面上,并把pi平面及各點(diǎn)的投影以圖片的形式展示出來。
在這里先感謝大佬們的回復(fù),您的每一個回復(fù)都是對我的問題的思考,小弟不勝感激!
在彈塑性力學(xué)這方面我學(xué)的還比較淺,以上討論部分算是我學(xué)習(xí)了這部分內(nèi)容后的一個總結(jié)概括吧,有什么不對的地方請批評指正。
返回小木蟲查看更多
京公網(wǎng)安備 11010802022153號
我該怎么回收金幣呢?
不知道
您問的太深了,
個人覺得,知道主應(yīng)力應(yīng)該沒辦法反推原來的應(yīng)力張量,同一個主應(yīng)力狀態(tài)可以對應(yīng)很多種應(yīng)力張量。
關(guān)于π平面上的投影,可以通過這兩個公式實(shí)現(xiàn)Xs=1/sqrt(6)*(σ1-σ3),Ys=1/sqrt(6)*(2σ2-σ1-σ3),