




版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、精選優(yōu)質(zhì)文檔-傾情為你奉上一 問題描述如圖所示的平面矩形結(jié)構(gòu),設(shè)E=1,NU=0.25,h=1,考慮以下約束和外載:位移邊界條件BC(u):UA=0,VA=0,UD=0,力邊界條件BC(p):在CD邊上有均布載荷q=1,建模情形:使用四個(gè)四節(jié)點(diǎn)矩形單元,試在該建模情形下,求各節(jié)點(diǎn)的位移以及各個(gè)單元的應(yīng)力分布。二 Matlab程序(1).函數(shù)定義:function k= Quad2D4Node_Stiffness(E,NU,h,xi,yi,xj,yj,xm,ym,xp,yp,ID) syms s t; a = (yi*(s-1)+yj*(-1-s)+ym*(1+s)+yp*(1-s)/4; b
2、= (yi*(t-1)+yj*(1-t)+ym*(1+t)+yp*(-1-t)/4; c = (xi*(t-1)+xj*(1-t)+xm*(1+t)+xp*(-1-t)/4; d = (xi*(s-1)+xj*(-1-s)+xm*(1+s)+xp*(1-s)/4; B1 = a*(t-1)/4-b*(s-1)/4 0 ; 0 c*(s-1)/4-d*(t-1)/4 ; c*(s-1)/4-d*(t-1)/4 a*(t-1)/4-b*(s-1)/4; B2 = a*(1-t)/4-b*(-1-s)/4 0 ; 0 c*(-1-s)/4-d*(1-t)/4 ; c*(-1-s)/4-d*(1-t)
3、/4 a*(1-t)/4-b*(-1-s)/4; B3 = a*(t+1)/4-b*(s+1)/4 0 ; 0 c*(s+1)/4-d*(t+1)/4 ; c*(s+1)/4-d*(t+1)/4 a*(t+1)/4-b*(s+1)/4; B4 = a*(-1-t)/4-b*(1-s)/4 0 ; 0 c*(1-s)/4-d*(-1-t)/4 ; c*(1-s)/4-d*(-1-t)/4 a*(-1-t)/4-b*(1-s)/4; Bfirst = B1 B2 B3 B4; Jfirst = 0 1-t t-s s-1 ; t-1 0 s+1 -s-t ; s-t -s-1 0 t+1 ; 1-
4、s s+t -t-1 0; J = xi xjxmxp*Jfirst*yi ;yj ; ym ; yp/8; B = Bfirst/J; if ID = 1 D = (E/(1-NU*NU)*1 NU 0 ; NU 1 0 ; 0 0 (1-NU)/2; elseif ID = 2 D = (E/(1+NU)/(1-2*NU)*1-NU NU 0 ; NU 1-NU 0 ; 0 0 (1-2*NU)/2; endBD = J*transpose(B)*D*B; r = int(int(BD, t, -1, 1), s, -1, 1); z = h*r; k = double(z); endfu
5、nction z = Quad2D4Node_Assembly(KK,k,i,j,m,p) DOF(1)=2*i-1; DOF(2)=2*i; DOF(3)=2*j-1; DOF(4)=2*j; DOF(5)=2*m-1; DOF(6)=2*m; DOF(7)=2*p-1; DOF(8)=2*p; for n1=1:8 for n2=1:8 KK(DOF(n1),DOF(n2)= KK(DOF(n1),DOF(n2)+k(n1,n2); endendz=KK;endfunction stress= Quad2D4Node_Stress(E,NU,xi,yi,xj,yj,xm,ym,xp,yp,
6、u,ID) syms s t; a = (yi*(s-1)+yj*(-1-s)+ym*(1+s)+yp*(1-s)/4; b = (yi*(t-1)+yj*(1-t)+ym*(1+t)+yp*(-1-t)/4; c = (xi*(t-1)+xj*(1-t)+xm*(1+t)+xp*(-1-t)/4; d = (xi*(s-1)+xj*(-1-s)+xm*(1+s)+xp*(1-s)/4; B1 = a*(t-1)/4-b*(s-1)/4 0 ; 0 c*(s-1)/4-d*(t-1)/4 ; c*(s-1)/4-d*(t-1)/4 a*(t-1)/4-b*(s-1)/4; B2 = a*(1-
7、t)/4-b*(-1-s)/4 0 ; 0 c*(-1-s)/4-d*(1-t)/4 ; c*(-1-s)/4-d*(1-t)/4 a*(1-t)/4-b*(-1-s)/4; B3 = a*(t+1)/4-b*(s+1)/4 0 ; 0 c*(s+1)/4-d*(t+1)/4 ; c*(s+1)/4-d*(t+1)/4 a*(t+1)/4-b*(s+1)/4; B4 = a*(-1-t)/4-b*(1-s)/4 0 ; 0 c*(1-s)/4-d*(-1-t)/4 ; c*(1-s)/4-d*(-1-t)/4 a*(-1-t)/4-b*(1-s)/4; Bfirst = B1 B2 B3 B4
8、; Jfirst = 0 1-t t-s s-1 ; t-1 0 s+1 -s-t ; s-t -s-1 0 t+1 ; 1-s s+t -t-1 0; J = xi xjxmxp*Jfirst*yi ;yj ; ym ; yp/8; B = Bfirst/J; if ID = 1 D = (E/(1-NU*NU)*1 NU 0 ; NU 1 0 ; 0 0 (1-NU)/2; elseif ID = 2 D = (E/(1+NU)/(1-2*NU)*1-NU NU 0 ; NU 1-NU 0 ; 0 0 (1-2*NU)/2; endstr1 = D*B*u; str2 = subs(str
9、1, s,t, 0,0); stress = double(str2); end(2). 計(jì)算部分E=1;NU=0.25;h=1;ID=1;k1= Quad2D4Node_Stiffness(E,NU,h,1,1,0.5,1,0.5,0.5,1,0.5,ID);k2= Quad2D4Node_Stiffness(E,NU,h,1,0.5,0.5,0.5,0.5,0,1,0,ID);k3= Quad2D4Node_Stiffness(E,NU,h,0.5,1,0,1,0,0.5,0.5,0.5,ID);k4= Quad2D4Node_Stiffness(E,NU,h,0.5,0.5,0,0.5,
10、0,0,0.5,0,ID);KK=zeros(18,18);KK= Quad2D4Node_Assembly(KK,k1,1,6,5,2);KK= Quad2D4Node_Assembly(KK,k2,2,5,4,3);KK= Quad2D4Node_Assembly(KK,k3,6,7,8,5);KK= Quad2D4Node_Assembly(KK,k4,5,8,9,4)k=KK(1:12,14:16,1:12,14:16);p=0;-0.25;0;0;0;0;0;0;0;0;0;-0.5;-0.25;0;0;u=kpU=u(1:12);0;u(13:15);0;0;u1=U(1);U(2
11、);U(11);U(12);U(9);U(10);U(3);U(4);stress1=Quad2D4Node_Stress(E,NU, 1,1,0.5,1,0.5,0.5,1,0.5,u1,ID)u2=U(3);U(4);U(9);U(10);U(7);U(8);U(5);U(6);stress2=Quad2D4Node_Stress(E,NU, 1,0.5,0.5,0.5,0.5,0,1,0,u2,ID)u3=U(11);U(12);U(13);U(14);U(15);U(16);U(9);U(10);stress3=Quad2D4Node_Stress(E,NU, 0.5,1,0,1,0,
12、0.5,0.5,0.5,u3,ID)u4=U(9);U(10);U(15);U(16);U(17);U(18);U(7);U(8);stress4=Quad2D4Node_Stress(E,NU, 0.5,0.5,0,0.5,0,0,0.5,0,u4,ID)專心-專注-專業(yè)總體剛度矩陣:各節(jié)點(diǎn)位移:各單元應(yīng)力:三 結(jié)果各個(gè)節(jié)點(diǎn)位移:u1=1.5749,v1=-4.5116,u2=0.5858,v2=-4.2489,u3=-0.4401,v3=-4.1495,u4=1.1458,v4=-3.3911,u5=0.7035,v5=-2.9251,u6=-0.4105,v6=-3.0964,u7=0,v7= -3.0486,u8=0.6532,v8=-1.9914,u9=0,v9=0。
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 電力本質(zhì)安全課件
- 2024模具設(shè)計(jì)師考試反思與改進(jìn)策略及試題答案
- 提升體育經(jīng)紀(jì)人考試能力的關(guān)鍵因素 試題及答案
- 2024模具設(shè)計(jì)師考試資料及試題答案
- 2024年模具設(shè)計(jì)師資格認(rèn)證考試全書試題及答案
- 模具設(shè)計(jì)的質(zhì)量管理體系試題及答案
- 2024年體育經(jīng)紀(jì)人考試的思維導(dǎo)圖及試題及答案
- 2024年云南省考公務(wù)員考試結(jié)構(gòu)化無領(lǐng)導(dǎo)小組面試真題試題試卷答案解析
- 2024年體育經(jīng)紀(jì)人考試亮點(diǎn)試題及答案
- 模具設(shè)計(jì)師資格考試數(shù)據(jù)分析試題及答案
- 第47屆世界技能大賽江蘇省選拔賽平面設(shè)計(jì)技術(shù)項(xiàng)目技術(shù)工作文件
- 2024年網(wǎng)絡(luò)與信息安全考試題庫
- 安橋功放機(jī)TX-NR3010說明書
- 《畜禽糞肥還田利用技術(shù)規(guī)范(征求意見稿)》編制說明
- GB/T 44309-2024陶瓷巖板
- 小學(xué)五年級(jí)下學(xué)期科學(xué)《我們面臨的環(huán)境問題》教學(xué)課件
- 血透病人低血壓護(hù)理查房
- 2024年工程承包合同書范文
- 有限空間作業(yè)風(fēng)險(xiǎn)辨識(shí)管控臺(tái)帳
- JGJT397-2016 公墓和骨灰寄存建筑設(shè)計(jì)規(guī)范
- 拖拉機(jī)濕式離合器
評(píng)論
0/150
提交評(píng)論