版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)
文檔簡介
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,建模情形:使用四個四節(jié)點矩形單元,試在該建模情形下,求各節(jié)點的位移以及各個單元的應(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). 計算部分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é)點位移:各單元應(yīng)力:三 結(jié)果各個節(jié)點位移: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等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 初中語文++《我為什么而活著》課件+八年級語文上冊同步備課精講精練(統(tǒng)編版)
- 病毒及其疫苗的生產(chǎn)制備技術(shù)
- 【生物】2024年高考天津生物試卷(已校對)
- 4.7 相似三角形的性質(zhì) 數(shù)學(xué)北師大版九年級上冊同步測試(含答案)
- 貴州省畢節(jié)市部分學(xué)校2023-2024學(xué)年高一下學(xué)期期末聯(lián)考語文試題(解析版)
- 2024年黑龍江省牡丹江市中考英語真題含解析
- 電商行業(yè)的消費者購買決策影響因素培訓(xùn)內(nèi)容策劃
- 2023年鎮(zhèn)江市第一人民醫(yī)院招聘考試試題及答案
- 2023年鹽城市濱??h農(nóng)村訂單定向醫(yī)學(xué)生定向招聘筆試真題
- 2023年大慶市大同區(qū)招聘公益性崗位人員考試試題及答案
- 1 丁烯安全技術(shù)說明書
- 《統(tǒng)計計算與軟件》課程教學(xué)大綱
- 土地估價師土地估價實務(wù)考試寶典
- 三級醫(yī)師查房登記本(共3頁)
- 健康中國建設(shè)與大健康產(chǎn)業(yè)發(fā)展
- 耳穴的功能歸類課件
- 全國醫(yī)療服務(wù)價格項目規(guī)范(2012版)
- 廣東省二手車買賣合同(示范文本)
- 小學(xué)生美術(shù)課件-第5課-我和動物交朋友-冀美版-(共26張PPT)ppt課件
- 魯班獎創(chuàng)優(yōu)項目裝修工程細部工藝細則
- 工程預(yù)決算高估冒算的表現(xiàn)形式及審計對策探討
評論
0/150
提交評論