單元?jiǎng)偠染仃?等參元)MATLAB編程

1、精選優(yōu)質(zhì)文檔-傾情為你奉上有限元法實(shí)驗(yàn)報(bào)告 專業(yè)班級(jí) 力學(xué)（實(shí)驗(yàn)）1601 姓 名 田詩豪 學(xué) 號(hào) 提交日期 2019.4.24 實(shí)驗(yàn)編號(hào)實(shí)驗(yàn)一實(shí)驗(yàn)二實(shí)驗(yàn)三總分得 分實(shí) 驗(yàn) 一（30分）一、實(shí)驗(yàn)內(nèi)容編寫一個(gè)計(jì)算平面3結(jié)點(diǎn)三角形單元的應(yīng)變矩陣、應(yīng)力矩陣和單元?jiǎng)偠染仃嚨腗ATLAB函數(shù)文件B3,S3,K3 = ele_mat_tri3(xy3,mat)，其中：輸入變量xy3為結(jié)點(diǎn)坐標(biāo)數(shù)組，mat為材料參數(shù)矩陣；輸出變量B3為應(yīng)變矩陣，S3為應(yīng)力矩陣，K3為單元?jiǎng)偠染仃?。（要求給出3個(gè)不同算例進(jìn)行驗(yàn)證，并繪制出單元形狀和結(jié)點(diǎn)號(hào)）二、程序代碼l 通用函數(shù)function B3,S3,K3 = ele_

2、mat_tri3(xy3,mat)%生成平面3結(jié)點(diǎn)三角形單元的應(yīng)變矩陣、應(yīng)力矩陣和單元?jiǎng)偠染仃嚨墓δ芎瘮?shù)%*變量說明*%xy3-結(jié)點(diǎn)坐標(biāo)數(shù)組%mat-材料參數(shù)矩陣（彈性模量，泊松比，壁厚）%B3-應(yīng)變矩陣%S3-應(yīng)力矩陣%K3-單元?jiǎng)偠染仃?*xyh=1,xy3(1,1),xy3(1,2);1,xy3(2,1),xy3(2,2);1,xy3(3,1),xy3(3,2);A=0.5*det(xyh);A=abs(A);D=mat(1)/(1-mat(2)2)*1,mat(2),0;mat(2),1,0;0,0,(1-mat(2)/2;b=zeros(1,3);c=zeros(1,3);%*for

3、 i=1:3 if i=1 j=2; m=3; elseif i=2 j=3; m=1; else j=1; m=2; end b(i)=xy3(j,2)-xy3(m,2); c(i)=xy3(m,1)-xy3(j,1);end%*B31=1/(2*A)*b(1),0;0,c(1);c(1),b(1);B32=1/(2*A)*b(2),0;0,c(2);c(2),b(2);B33=1/(2*A)*b(3),0;0,c(3);c(3),b(3);B3=B31,B32,B33;%*S3=D*B3;%*K3=A*mat(3)*B3'*D*B3;l 主程序clear;clc;%*輸入結(jié)點(diǎn)坐標(biāo)數(shù)

4、組*xy3=0,0;5,1;1,4; mat=3e6,0.5,1.0; %*輸入材料參數(shù)矩陣(彈性模量，泊松比，壁厚)*B3,S3,K3=ele_mat_tri3(xy3,mat)三、算例分析l 算例1：如圖1所示三角形單元，結(jié)點(diǎn)坐標(biāo)為1(0,0)，2(5,2)，3(1,4)，彈性模量為200GPa，泊松比為0.35、厚度為0.5m。試求應(yīng)變矩陣，應(yīng)力矩陣和單元?jiǎng)偠染仃嚒D1 算例1三角形單元解：根據(jù)如圖1所示三角形單元及其幾何和材料參數(shù)，編制主程序如下：clear;clc;%*輸入結(jié)點(diǎn)坐標(biāo)數(shù)組*xy3=0,0;5,2;1,4; mat=2e11,0.35,0.5; %*輸入材料參數(shù)矩陣(彈性

5、模量，泊松比，壁厚)*B3,S3,K3=ele_mat_tri3(xy3,mat)運(yùn)行程序，得到應(yīng)變矩陣B3如下：-0.1111 0.0000 0.2222 0.0000 -0.1111 0.0000 0.0000 -0.2222 0.0000 -0.0556 0.0000 0.2778 -0.2222 -0.1111 -0.0556 0.2222 0.2778 -0.1111 得到應(yīng)力矩陣S3（Pa）如下：-2.53E+10-1.77E+105.06E+10-4.43E+09-2.53E+102.22E+10-8.86E+09-5.06E+101.77E+10-1.27E+10-8.86E+

6、096.33E+10-1.65E+10-8.23E+09-4.12E+091.65E+102.06E+10-8.23E+09得到單元?jiǎng)偠染仃嘖3（Pa）如下：2.91E+101.71E+10-2.12E+10-1.42E+10-7.91E+09-2.85E+091.71E+105.48E+10-1.57E+104.43E+09-1.42E+09-5.92E+10-2.12E+10-1.57E+105.17E+10-8.55E+09-3.05E+102.42E+10-1.42E+104.43E+09-8.55E+091.96E+102.28E+10-2.41E+10-7.91E+09-1.42E

7、+09-3.05E+102.28E+103.84E+10-2.14E+10-2.85E+09-5.92E+102.42E+10-2.41E+10-2.14E+108.33E+10l 算例2：如圖2所示三角形單元，結(jié)點(diǎn)坐標(biāo)為1(0,0)，2(3,0)，3(0,5)，彈性模量為200GPa，泊松比為0.35、厚度為0.5m。試求應(yīng)變矩陣，應(yīng)力矩陣和單元?jiǎng)偠染仃?。圖2 算例2三角形單元解：根據(jù)如圖2所示三角形單元及其幾何和材料參數(shù)，編制主程序如下：clear;clc;%*輸入結(jié)點(diǎn)坐標(biāo)數(shù)組*xy3=0,0;3,0;5,0; mat=2e11,0.35,0.5; %*輸入材料參數(shù)矩陣(彈性模量，泊松比，

8、壁厚)*B3,S3,K3=ele_mat_tri3(xy3,mat)運(yùn)行程序，得到應(yīng)變矩陣B3如下：-0.3333 0.0000 0.3333 0.0000 0.0000 0.0000 0.0000 -0.2000 0.0000 0.0000 0.0000 0.2000 -0.2000 -0.3333 0.0000 0.3333 0.2000 0.0000 得到應(yīng)力矩陣S3（Pa）如下：-7.60E+10-1.60E+107.60E+100.00E+000.00E+001.60E+10-2.66E+10-4.56E+102.66E+100.00E+000.00E+004.56E+10-1.48

9、E+10-2.47E+100.00E+002.47E+101.48E+100.00E+00得到單元?jiǎng)偠染仃嘖3（Pa）如下：1.06E+113.85E+10-9.50E+10-1.85E+10-1.11E+10-1.99E+103.85E+106.51E+10-1.99E+10-3.09E+10-1.85E+10-3.42E+10-9.50E+10-1.99E+109.50E+100.00E+000.00E+001.99E+10-1.85E+10-3.09E+100.00E+003.09E+101.85E+100.00E+00-1.11E+10-1.85E+100.00E+001.85E+10

10、1.11E+100.00E+00-1.99E+10-3.42E+101.99E+100.00E+000.00E+003.42E+10l 算例3：如圖3所示三角形單元，結(jié)點(diǎn)坐標(biāo)為1(0,0)，2(3,0)，3(1.5,1.53)，彈性模量為200GPa，泊松比為0.35、厚度為0.5m。試求應(yīng)變矩陣，應(yīng)力矩陣和單元?jiǎng)偠染仃?。圖3 算例3三角形單元解：根據(jù)如圖3所示三角形單元及其幾何和材料參數(shù)，編制主程序如下：clear;clc;%*輸入結(jié)點(diǎn)坐標(biāo)數(shù)組*xy3=0,0;3,0;1.5,1.5*sqrt(3); mat=2e11,0.35,0.5; %*輸入材料參數(shù)矩陣(彈性模量，泊松比，壁厚)*B3

11、,S3,K3=ele_mat_tri3(xy3,mat)運(yùn)行程序，得到應(yīng)變矩陣B3如下：-0.3333 0.0000 0.3333 0.0000 0.0000 0.0000 0.0000 -0.1925 0.0000 -0.1925 0.0000 0.3849 -0.1925 -0.3333 -0.1925 0.3333 0.3849 0.0000 得到應(yīng)力矩陣S3（Pa）如下：-7.60E+10-1.54E+107.60E+10-1.54E+100.00E+003.07E+10-2.66E+10-4.39E+102.66E+10-4.39E+100.00E+008.77E+10-1.43E+

12、10-2.47E+10-1.43E+102.47E+102.85E+100.00E+00得到單元?jiǎng)偠染仃嘖3（Pa）如下：5.47E+101.92E+10-4.40E+107.12E+08-1.07E+10-1.99E+101.92E+103.25E+10-7.12E+084.11E+08-1.85E+10-3.29E+10-4.40E+10-7.12E+085.47E+10-1.92E+10-1.07E+101.99E+107.12E+084.11E+08-1.92E+103.25E+101.85E+10-3.29E+10-1.07E+10-1.85E+10-1.07E+101.85E+10

13、2.14E+100.00E+00-1.99E+10-3.29E+101.99E+10-3.29E+100.00E+006.58E+10實(shí) 驗(yàn) 二（30分）一、實(shí)驗(yàn)內(nèi)容編寫一個(gè)計(jì)算平面4結(jié)點(diǎn)四邊形等參元的剛度矩陣的MATLAB函數(shù)文件K4 = ele_mat_quad4(xy4,mat)，其中：輸入變量xy4為結(jié)點(diǎn)坐標(biāo)數(shù)組，mat為材料參數(shù)矩陣；輸出變量K4為單元?jiǎng)偠染仃?。（要求給出3個(gè)不同算例進(jìn)行驗(yàn)證，并繪制出單元形狀和結(jié)點(diǎn)號(hào)）二、程序代碼l 通用函數(shù)function K4 = ele_mat_quad4(xy4,mat)%生成平面4結(jié)點(diǎn)四邊形等參元的剛度矩陣的功能函數(shù)%*變量說明*%xy4-

14、結(jié)點(diǎn)坐標(biāo)數(shù)組%mat-材料參數(shù)矩陣（彈性模量，泊松比，壁厚）%K4-單元?jiǎng)偠染仃?y1y2-局部坐標(biāo)系結(jié)點(diǎn)坐標(biāo)%D-彈性矩陣%*y1y2=-1,-1;1,-1;1,1;-1,1;D=mat(1)/(1-mat(2)2)*1,mat(2),0;mat(2),1,0;0,0,(1-mat(2)/2;%*數(shù)值積分(Guass, n=4)*C(1)=0.4055;C(2)=0.4886;C(3)=-0.4886;C(4)=-0.4055;A(1)=0.7454;A(2)=0.2546;A(3)=0.2546;A(4)=0.7454;sum=0;for i = 1:4 for j = 1:4 y1=C;

15、y2=C; k=1:4;%* PN1(1)=0.25*(y2(j)-1);PN2(1)=0.25*(y1(i)-1); PN1(2)=-0.25*(y2(j)-1);PN2(2)=-0.25*(y1(i)+1); PN1(3)=0.25*(y2(j)+1);PN2(3)=0.25*(y1(i)+1); PN1(4)=-0.25*(y2(j)+1);PN2(4)=-0.25*(y1(i)-1); for k=1:4 PN(:,k)=PN1(k),PN2(k)' end J=PN*xy4; JN=inv(J); J1=JN(1,:);J2=JN(2,:);%*應(yīng)變矩陣* B1=J1*PN(

16、:,1),0;0,J2*PN(:,1);J2*PN(:,1),J1*PN(:,1); B2=J1*PN(:,2),0;0,J2*PN(:,2);J2*PN(:,2),J1*PN(:,2); B3=J1*PN(:,3),0;0,J2*PN(:,3);J2*PN(:,3),J1*PN(:,3); B4=J1*PN(:,4),0;0,J2*PN(:,4);J2*PN(:,4),J1*PN(:,4); B=B1,B2,B3,B4;%* m=mat(3)*B'*D*B*det(J); sum = sum+A(i)*A(j)*m; endendK4=vpa(sum,9);l 主程序clc;clea

17、r;%*輸入結(jié)點(diǎn)坐標(biāo)數(shù)組*xy4=3,2;8,3;7,8;4,7; %*輸入材料參數(shù)矩陣(彈性模量，泊松比，壁厚)*mat=3e6,0.5,1.0; K4=ele_mat_quad4(xy4,mat)三、算例分析l 算例1：如圖5所示四邊形等參單元（圖4為局部坐標(biāo)系規(guī)則單元），已知4個(gè)結(jié)點(diǎn)整體坐標(biāo)系內(nèi)的坐標(biāo)為1(3,2); 2(8,3); 3(7,8); 4(4,7)。彈性模量為200GPa，泊松比為0.35、厚度為0.05m。試求單元?jiǎng)偠染仃?。圖4 局部坐標(biāo)系規(guī)則單元圖5 算例1四邊形等參單元解：根據(jù)如圖5所示四邊形等參單元及其幾何和材料參數(shù)，編制主程序如下：clc;clear;%*輸入結(jié)點(diǎn)

18、坐標(biāo)數(shù)組*xy4=3,2;8,3;7,8;4,7; %*輸入材料參數(shù)矩陣(彈性模量，泊松比，壁厚)*mat=2e11,0.35,0.5; K4=ele_mat_quad4(xy4,mat)運(yùn)行程序，得到單元?jiǎng)偠染仃嘖4（Pa）如下：4.01E+091.45E+09-3.54E+09-2.37E+08-1.54E+09-1.67E+091.08E+094.56E+081.45E+093.79E+09-3.79E+084.98E+08-1.67E+09-1.87E+095.98E+08-2.41E+09-3.54E+09-3.79E+086.62E+09-2.37E+091.38E+095.71E

19、+08-4.47E+092.18E+09-2.37E+084.98E+08-2.37E+094.53E+094.28E+08-2.17E+092.18E+09-2.86E+09-1.54E+09-1.67E+091.38E+094.28E+085.24E+091.34E+09-5.08E+09-9.97E+07-1.67E+09-1.87E+095.71E+08-2.17E+091.34E+094.74E+09-2.42E+08-6.98E+081.08E+095.98E+08-4.47E+092.18E+09-5.08E+09-2.42E+088.47E+09-2.54E+094.56E+0

20、8-2.41E+092.18E+09-2.86E+09-9.97E+07-6.98E+08-2.54E+095.97E+09l 算例2：如圖6所示四邊形等參單元，已知4個(gè)結(jié)點(diǎn)整體坐標(biāo)系內(nèi)的坐標(biāo)為1(1,1); 2(4,1); 3(4,3); 4(1,4)。彈性模量為180GPa，泊松比為0.5、厚度為0.1m。試求單元?jiǎng)偠染仃?。圖6 算例2四邊形等參單元解：根據(jù)如圖6所示四邊形等參單元及其幾何和材料參數(shù)，編制主程序如下：clc;clear;%*輸入結(jié)點(diǎn)坐標(biāo)數(shù)組*xy4=1,1;4,1;4,3;1,3; %*輸入材料參數(shù)矩陣(彈性模量，泊松比，壁厚)*mat=1.8e11,0.5,0.1; K4

21、=ele_mat_quad4(xy4,mat)運(yùn)行程序，得到單元?jiǎng)偠染仃嘖4（Pa）如下：8.33E+094.50E+09-3.83E+091.50E+09-4.17E+09-4.50E+09-3.33E+08-1.50E+094.50E+091.33E+10-1.50E+094.67E+09-4.50E+09-6.67E+091.50E+09-1.13E+10-3.83E+09-1.50E+098.33E+09-4.50E+09-3.33E+081.50E+09-4.17E+094.50E+091.50E+094.67E+09-4.50E+091.33E+10-1.50E+09-1.13E+

22、104.50E+09-6.67E+09-4.17E+09-4.50E+09-3.33E+08-1.50E+098.33E+094.50E+09-3.83E+091.50E+09-4.50E+09-6.67E+091.50E+09-1.13E+104.50E+091.33E+10-1.50E+094.67E+09-3.33E+081.50E+09-4.17E+094.50E+09-3.83E+09-1.50E+098.33E+09-4.50E+09-1.50E+09-1.13E+104.50E+09-6.67E+091.50E+094.67E+09-4.50E+091.33E+10l 算例3：如

23、圖7所示四邊形等參單元，已知4個(gè)結(jié)點(diǎn)整體坐標(biāo)系內(nèi)的坐標(biāo)為1(1,1); 2(4,1); 3(5,3); 4(2,3)。彈性模量為180GPa，泊松比為0.5、厚度為0.1m。試求單元?jiǎng)偠染仃?。圖7 算例3四邊形等參單元解：根據(jù)如圖7所示四邊形等參單元及其幾何和材料參數(shù)，編制主程序如下：clc;clear;%*輸入結(jié)點(diǎn)坐標(biāo)數(shù)組*xy4=1,1;4,1;5,3;2,3; %*輸入材料參數(shù)矩陣(彈性模量，泊松比，壁厚)*mat=1.8e11,0.5,0.1; K4=ele_mat_quad4(xy4,mat)運(yùn)行程序，得到單元?jiǎng)偠染仃嘖4（Pa）如下：7.17E+092.50E+09-4.17E+0

24、93.50E+09-2.83E+09-3.50E+09-1.67E+08-2.50E+092.50E+098.67E+095.00E+083.33E+09-3.50E+09-1.33E+095.00E+08-1.07E+10-4.17E+095.00E+081.02E+10-6.50E+09-1.67E+085.00E+08-5.83E+095.50E+093.50E+093.33E+09-6.50E+092.07E+10-2.50E+09-1.07E+105.50E+09-1.33E+10-2.83E+09-3.50E+09-1.67E+08-2.50E+097.17E+092.50E+09

25、-4.17E+093.50E+09-3.50E+09-1.33E+095.00E+08-1.07E+102.50E+098.67E+095.00E+083.33E+09-1.67E+085.00E+08-5.83E+095.50E+09-4.17E+095.00E+081.02E+10-6.50E+09-2.50E+09-1.07E+105.50E+09-1.33E+103.50E+093.33E+09-6.50E+092.07E+10實(shí) 驗(yàn) 三（40分）1、實(shí)驗(yàn)內(nèi)容編寫一個(gè)計(jì)算平面8結(jié)點(diǎn)四邊形等參元?jiǎng)偠染仃嚨腗ATLAB函數(shù)文件K8 = ele_mat_quad8(xy8,mat)，其中：輸

26、入變量xy8為結(jié)點(diǎn)坐標(biāo)數(shù)組，mat為材料參數(shù)矩陣；輸出變量K8為單元?jiǎng)偠染仃?。（要求給出3個(gè)不同算例進(jìn)行驗(yàn)證，并繪制出單元形狀和結(jié)點(diǎn)號(hào)）2、程序代碼l 通用函數(shù)function K8 = ele_mat_quad8(xy8,mat)%生成平面8結(jié)點(diǎn)四邊形等參元?jiǎng)偠染仃嚨墓δ芎瘮?shù)%*變量說明*%xy8-結(jié)點(diǎn)坐標(biāo)數(shù)組%mat-材料參數(shù)矩陣（彈性模量，泊松比，壁厚）%K8-單元?jiǎng)偠染仃?D-彈性矩陣%*PN = sym(zeros(2,8);D=mat(1)/(1-mat(2)2)*1,mat(2),0;mat(2),1,0;0,0,(1-mat(2)/2;syms y1 y2 realN(1)=(

27、1-y1)*(1-y2)*(-1-y1-y2)/4;N(2)=(1+y1)*(1-y2)*(-1-y2+y1)/4;N(3)=(1+y1)*(1+y2)*(-1+y2+y1)/4;N(4)=(1-y1)*(1+y2)*(-1+y2-y1)/4;N(5)=(1-y12)*(1-y2)/2;N(6)=(1-y22)*(1+y1)/2;N(7)=(1-y12)*(1+y2)/2;N(8)=(1-y1)*(1-y22)/2;%*i=1:8;PN1(i)=diff(N(i),y1);PN2(i)=diff(N(i),y2);for i=1:8PN(:,i)=PN1(i),PN2(i)'endJ=

28、PN*xy8;JN=inv(J);J1=JN(1,:);J2=JN(2,:);%*應(yīng)變矩陣*B1=J1*PN(:,1),0;0,J2*PN(:,1);J2*PN(:,1),J1*PN(:,1);B2=J1*PN(:,2),0;0,J2*PN(:,2);J2*PN(:,2),J1*PN(:,2);B3=J1*PN(:,3),0;0,J2*PN(:,3);J2*PN(:,3),J1*PN(:,3);B4=J1*PN(:,4),0;0,J2*PN(:,4);J2*PN(:,4),J1*PN(:,4);B5=J1*PN(:,5),0;0,J2*PN(:,5);J2*PN(:,5),J1*PN(:,5)

29、;B6=J1*PN(:,6),0;0,J2*PN(:,6);J2*PN(:,6),J1*PN(:,6);B7=J1*PN(:,7),0;0,J2*PN(:,7);J2*PN(:,7),J1*PN(:,7);B8=J1*PN(:,8),0;0,J2*PN(:,8);J2*PN(:,8),J1*PN(:,8);B=B1,B2,B3,B4,B5,B6,B7,B8;%*m=mat(3)*B'*D*B*det(J);%*數(shù)值積分(Guass, n=4)*C(1)=0.4055; C(2)=0.4886;C(3)=-0.4886; C(4)=-0.4055;A(1)=0.7454; A(2)=0.

30、2546;A(3)=0.2546; A(4)=0.7454;sum=0;for i = 1:4 for j = 1:4 sum = sum+A(i)*A(j)*subs(m,y1,y2,C(i),C(j); endendK8=vpa(sum,9);l 主程序clc;clear;%*輸入結(jié)點(diǎn)坐標(biāo)數(shù)組*xy8=2,2;7,2;6,7;3,6;5,3;6,4;4,6;3,4; %*輸入材料參數(shù)矩陣(彈性模量，泊松比，壁厚)*mat=3e6,0.5,1.0; K8 = ele_mat_quad8(xy8,mat)3、算例分析l 算例1：如圖9所示曲四邊形等參單元（圖8為局部坐標(biāo)系規(guī)則單元），已知8個(gè)結(jié)

31、點(diǎn)整體坐標(biāo)系內(nèi)的坐標(biāo)為1(2,2); 2(7,2); 3(6,7); 4(3,6); 5(5,3); 6(6,4); 7(4,6); 8(3,4)。彈性模量為100GPa，泊松比為0.35、厚度為0.05m。試求單元?jiǎng)偠染仃?。圖8 局部坐標(biāo)系規(guī)則單元圖9 算例1曲四邊形等參單元解：根據(jù)如圖9所示曲四邊形等參單元及其幾何和材料參數(shù)，編制主程序如下：clc;clear;%*輸入結(jié)點(diǎn)坐標(biāo)數(shù)組*xy8=2,2;7,2;6,7;3,6;5,3;6,4;4,6;3,4; %*輸入材料參數(shù)矩陣(彈性模量，泊松比，壁厚)*mat=1.8e11,0.5,0.1; K8 = ele_mat_quad8(xy8,m

32、at)運(yùn)行程序，得到單元?jiǎng)偠染仃嘖8（GPa）如下：3.79 1.47 1.00 -0.10 0.91 0.16 1.92 -0.12 -7.51 -3.49 4.99 2.55 -2.87 -0.50 -2.22 0.03 1.47 5.37 -0.08 0.68 0.16 1.42 -0.14 2.40 -3.59 -3.08 2.55 2.15 -0.50 -1.92 0.12 -7.03 1.00 -0.08 2.05 -0.55 1.87 0.32 2.00 -0.45 -1.78 1.38 -2.52 -0.49 -2.63 -0.13 0.02 0.00 -0.10 0.68 -

33、0.55 2.07 0.34 2.40 -0.45 1.50 1.48 2.24 -0.58 -6.13 -0.13 -1.23 0.00 -1.55 0.91 0.16 1.87 0.34 5.87 2.38 3.38 0.19 2.76 2.16 -6.15 -3.32 -7.41 -1.51 -1.21 -0.40 0.16 1.42 0.32 2.40 2.38 5.90 0.21 2.42 2.16 4.79 -3.22 -10.36 -1.61 -3.15 -0.40 -3.43 1.92 -0.14 2.00 -0.45 3.38 0.21 6.68 -1.15 -2.68 0.

34、08 -1.51 0.36 -9.18 -0.21 -0.61 1.30 -0.12 2.40 -0.45 1.50 0.19 2.42 -1.15 3.98 0.08 0.39 0.36 -3.02 -0.12 -3.51 1.21 -4.15 -7.51 -3.59 -1.78 1.48 2.76 2.16 -2.68 0.08 36.97 11.70 -29.09 -15.32 4.65 0.41 -3.32 3.08 -3.49 -3.08 1.38 2.24 2.16 4.79 0.08 0.39 11.70 31.58 -15.32 -29.16 0.41 -2.12 3.08 -

35、4.65 4.99 2.55 -2.52 -0.58 -6.15 -3.22 -1.51 0.36 -29.09 -15.32 37.14 14.78 -0.72 2.09 -2.14 -0.66 2.55 2.15 -0.49 -6.13 -3.32 -10.36 0.36 -3.02 -15.32 -29.16 14.78 40.11 2.09 1.91 -0.66 4.50 -2.87 -0.50 -2.63 -0.13 -7.41 -1.61 -9.18 -0.12 4.65 0.41 -0.72 2.09 18.97 1.45 -0.81 -1.60 -0.50 -1.92 -0.1

36、3 -1.23 -1.51 -3.15 -0.21 -3.51 0.41 -2.12 2.09 1.91 1.45 9.72 -1.60 0.30 -2.22 0.12 0.02 0.00 -1.21 -0.40 -0.61 1.21 -3.32 3.08 -2.14 -0.66 -0.81 -1.60 10.29 -1.74 0.03 -7.03 0.00 -1.55 -0.40 -3.43 1.30 -4.15 3.08 -4.65 -0.66 4.50 -1.60 0.30 -1.74 16.00 l 算例2：如圖10所示曲四邊形等參單元，已知8個(gè)結(jié)點(diǎn)整體坐標(biāo)系內(nèi)的坐標(biāo)為1(-5,-4)

37、; 2(6,-12); 3(10,8); 4(-8,12); 5(0,-6); 6(5,-1); 7(0,4); 8(-4,2)。彈性模量為200GPa，泊松比為0.35、厚度為0.1m。試求單元?jiǎng)偠染仃?。圖10 算例2曲四邊形等參單元解：根據(jù)如圖10所示曲四邊形等參單元及其幾何和材料參數(shù)，編制主程序如下：clc;clear;%*輸入結(jié)點(diǎn)坐標(biāo)數(shù)組*xy8=-5,-4;6,-12;10,8;-8,12;0,-6;5,-1;0,4;-4,2; %*輸入材料參數(shù)矩陣(彈性模量，泊松比，壁厚)*mat=2e11,0.35,0.1; K8 = ele_mat_quad8(xy8,mat)運(yùn)行程序，得到單

38、元?jiǎng)偠染仃嘖8（GPa）如下：40.8 11.0 16.5 1.2 10.6 2.9 11.6 2.9 -33.1 -0.6 -12.2 -4.7 -3.7 3.2 -30.6 -15.9 11.0 22.7 1.3 9.6 2.9 7.2 2.8 9.7 -0.9 -9.3 -4.7 -13.5 3.2 3.6 -15.5 -29.8 16.5 1.3 16.4 -3.9 9.4 1.1 5.5 -0.6 -26.8 3.3 -7.4 1.7 -13.1 -2.9 -0.5 0.0 1.2 9.6 -3.9 14.3 1.2 9.3 -0.6 5.3 3.7 -8.0 1.3 -22.2 -

39、2.9 -0.2 0.0 -8.0 10.6 2.9 9.4 1.2 20.3 6.4 3.9 0.8 -10.5 0.4 -16.1 -4.1 -45.0 -13.4 27.3 5.8 2.9 7.2 1.1 9.3 6.4 15.6 0.9 2.6 0.4 -2.7 -3.7 -26.5 -13.7 -12.5 5.8 7.0 11.6 2.8 5.5 -0.6 3.9 0.9 8.8 -0.8 -8.8 1.9 3.7 -5.2 -11.4 8.3 -13.3 -7.3 2.9 9.7 -0.6 5.3 0.8 2.6 -0.8 9.9 1.9 -1.6 -5.2 -3.0 8.7 -0

40、.4 -7.7 -22.5 -33.1 -0.9 -26.8 3.7 -10.5 0.4 -8.8 1.9 73.9 6.0 -19.9 -11.7 27.0 -2.3 -1.9 3.0 -0.6 -9.3 3.3 -8.0 0.4 -2.7 1.9 -1.6 6.0 40.4 -11.7 -9.3 -2.3 -3.7 3.0 -5.7 -12.2 -4.7 -7.4 1.3 -16.1 -3.7 3.7 -5.2 -19.9 -11.7 85.2 -7.6 -20.3 27.0 -12.9 4.5 -4.7 -13.5 1.7 -22.2 -4.1 -26.5 -5.2 -3.0 -11.7

41、 -9.3 -7.6 92.0 27.0 -31.2 4.5 13.7 -3.7 3.2 -13.1 -2.9 -45.0 -13.7 -11.4 8.7 27.0 -2.3 -20.3 27.0 211.6 28.8 -145.2 -48.8 3.2 3.6 -2.9 -0.2 -13.4 -12.5 8.3 -0.4 -2.3 -3.7 27.0 -31.2 28.8 123.5 -48.8 -79.1 -30.6 -15.5 -0.5 0.0 27.3 5.8 -13.3 -7.7 -1.9 3.0 -12.9 4.5 -145.2 -48.8 177.1 58.6 -15.9 -29.

42、8 0.0 -8.0 5.8 7.0 -7.3 -22.5 3.0 -5.7 4.5 13.7 -48.8 -79.1 58.6 124.4 l 算例3：如圖11所示曲四邊形等參單元，已知8個(gè)結(jié)點(diǎn)整體坐標(biāo)系內(nèi)的坐標(biāo)為1(-6,-4); 2(6,-4); 3(6,4); 4(-6,4); 5(0,-3); 6(3,0); 7(0,3); 8(-3,0)。彈性模量為200GPa，泊松比為0.35、厚度為0.1m。試求單元?jiǎng)偠染仃嚒D11 算例3曲四邊形等參單元解：根據(jù)如圖11所示曲四邊形等參單元及其幾何和材料參數(shù)，編制主程序如下：clc;clear;%*輸入結(jié)點(diǎn)坐標(biāo)數(shù)組*xy8=-6,-4;6,-4;6,4;-6,4;0,-3;3,0;0,3;-3,0; %*輸入材料參數(shù)矩陣(彈性模量，泊松比，壁厚)*mat=2e11,0.35,0.1; K8 = ele_mat_qua