有限元MATLAB講課教案.doc

《有限元MATLAB講課教案.doc》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。

1、有限元MATLAB精品文檔MATLAB報(bào)告Matlab程序求解簡要過程如下：（1）求取單元節(jié)點(diǎn)位移提取矩陣T單元節(jié)點(diǎn)位移提取矩陣T本質(zhì)上是置換矩陣群中的一個(gè)，結(jié)果可將任意雜亂的節(jié)點(diǎn)順序置換成統(tǒng)一的順序。另一方面其作用是對單元?jiǎng)偠染仃囘M(jìn)行“升維操作”，將單元?jiǎng)偠染仃嚱y(tǒng)籌到整體剛度矩陣上來，便于對總體節(jié)點(diǎn)位移矩陣和支座反力進(jìn)行求取。收集于網(wǎng)絡(luò)，如有侵權(quán)請聯(lián)系管理員刪除精品文檔本程序分析過程中對單元1的節(jié)點(diǎn)提取是按順序編號(hào)1-2-3，對單元2的節(jié)點(diǎn)提取是按順序編號(hào)2-3-4。單元1的節(jié)點(diǎn)位移提取矩陣如下：單元2的節(jié)點(diǎn)位移提

2、取矩陣如下：（1）求取單元幾何矩陣B單元1的節(jié)點(diǎn)按編號(hào)順序1-2-3分別進(jìn)行對幾何函數(shù)矩陣或算子矩陣的bi逆時(shí)針操作，對ci順時(shí)針操作；單元2的節(jié)點(diǎn)按編號(hào)順序2-3-4分別進(jìn)行對幾何函數(shù)矩陣的bi順時(shí)針操作，對ci逆時(shí)針操作.在MATLAB程序中通過mod()取模函數(shù)來達(dá)到對節(jié)點(diǎn)的順時(shí)針或逆時(shí)針循環(huán)操作。收集于網(wǎng)絡(luò)，如有侵權(quán)請聯(lián)系管理員刪除精品文檔單元1的幾何矩陣如下：單元2的幾何矩陣如下：（1）求取應(yīng)力矩陣S單元應(yīng)力矩陣滿足S=D*B,其中D為彈性矩陣，B為單元幾何矩陣各單元的彈性矩陣如下：單元1的應(yīng)力矩陣如下：收

3、集于網(wǎng)絡(luò)，如有侵權(quán)請聯(lián)系管理員刪除精品文檔單元2的應(yīng)力矩陣如下：（1）求取單元?jiǎng)偠染仃嘖單元?jiǎng)偠染仃嘖滿足公式K=B’*D*B*t*A,其中t為平面板的厚度，A為單元面積，且單元?jiǎng)偠染仃嚍閷ΨQ矩陣。單元1的剛度矩陣如下：單元2的剛度矩陣如下：（2）求取總體剛度矩陣sumKK收集于網(wǎng)絡(luò)，如有侵權(quán)請聯(lián)系管理員刪除精品文檔由上述步驟求得的單元?jiǎng)偠染仃嘖利用單元虛功原理和剛度方程可導(dǎo)出K’*δ=f，其中δ為單元節(jié)點(diǎn)位移列陣，f為單元等效節(jié)點(diǎn)載荷列陣，為了能將各個(gè)單元?jiǎng)偠确匠探y(tǒng)一到一個(gè)整體，便需要步驟（1）的單元節(jié)點(diǎn)提取矩陣對

4、單元?jiǎng)偠确匠踢M(jìn)行變換，將兩個(gè)變換結(jié)果聯(lián)立便得到總體剛度方程，其中也可得到總體剛度矩陣sumKK，且總體剛度矩陣可由sumKK=ΣT’*K*T求得?？傮w剛度矩陣如下：（1）求取總體節(jié)點(diǎn)位移矩陣和支座反力利用上述步驟提到的總體剛度方程sumKK*delta=F，其中delta為總體節(jié)點(diǎn)位移矩陣，F(xiàn)為總體等效節(jié)點(diǎn)載荷列陣。利用總體邊界條件可將總體剛度方程分成兩部分，一部分先求取總體節(jié)點(diǎn)位移矩陣delta，最后代入總體剛度方程求取總體等效節(jié)點(diǎn)載荷列陣F，由此便得到支座反力?？傮w節(jié)點(diǎn)位移矩陣如下：收集于網(wǎng)絡(luò)，如有侵權(quán)請聯(lián)系管理

5、員刪除精品文檔支座反力如：（1）求取單元應(yīng)力各單元應(yīng)力滿足θ=S*δ，其中S為單元應(yīng)力矩陣，δ為單元節(jié)點(diǎn)位移列陣。單元1的應(yīng)力：單元2的應(yīng)力：下圖為應(yīng)變前后的各單元對比圖：收集于網(wǎng)絡(luò)，如有侵權(quán)請聯(lián)系管理員刪除精品文檔收集于網(wǎng)絡(luò)，如有侵權(quán)請聯(lián)系管理員刪除