利用MATLAB件软对板进行结构分析.docVIP

一．算例 选取一块矩?形形薄板，荷载沿厚度?均匀分布，受100K?N/m均布荷载，板的?各部分尺寸如图所示。?板厚t=0.1，弹性?模量E=100，泊松?比=0.1，容重?=0 二．划分单?元 划分单元，标出单?元号码及节点码，选取?坐标，如上图示。 三?．程序选择 选择使用?MATLAB作为开发?程序，下面根据程序的?结构特点进行叙述 四?．子程序（m文件） ?子程序（m文件）在M?ATLAB程序中具有?重要的地位，通过它可?以方便的设置自己的函?数，这些函数可以在M?ATLAB命令窗口被?方便的调用。本例共使?用了五个子程序（m文?件），分别是： 1.?LinearTria?ngleElemen?tStiffness?(E,NU,t,xi?,yi,xj,yj,?xm,ym,p) 2?.LinearTri?angleAssem?ble(K,k,i,?j,m) 3. Li?nearTriang?leElementS?tresses(E,?NU,xi,yi,x?j,yj,xm,ym?,p,u） 4. L?inearTrian?gleElement?Pstresses(?sigma) 5.P?lotStress(?iStress) 具?体程序和解释说明如下?： 1. Linea?rTriangleE?lementStif?fness(E,NU?,t,xi,yi,x?j,yj,xm,ym?,p) 该函数用于计?算弹性模量为E、泊松?比为NU、厚度为t、?第一个节点坐标为(x?i,yi)、第二个节?点坐标为(xj,yj?)、第三个节点坐标为?(xm,ym)时的线?性三角形元的单元刚度?矩阵。P=1表明函数?用于平面应力情况。P?=2表明函数用于平面?应变情况。该函数返回?6*6的单元刚度矩阵?k。程序如下： fu?nction y =? LinearTri?angleEleme?ntStiffnes?s(E,NU,t,x?i,yi,xj,yj?,xm,ym,p) ?A=(xi*(yj-?ym)+xj*(ym?-yi)+xm*(y?i-yj))/2; ?%计算单元面积 be?tai=yj-ym;? ? ? %计?算单元应变矩阵中的元?素 betaj=ym?-yi; ? ? ? %计算单元应变矩?阵中的元素 beta?m=yi-yj; ? ? ? %计算单?元应变矩阵中的元素 ?gammai=xm-?xj; ? ? ?%计算单元应变矩阵中?的元素 gammaj?=xi-xm; ? ? ? %计算单元应?变矩阵中的元素 ga?mmam=xj-xi?; ? ? %计?算单元应变矩阵中的元?素 B=[betai? 0 betaj 0? betam 0; ? ? %形成单元应变矩?阵B 0 g?ammai 0 ga?mmaj 0 gam?mam; g?ammai beta?i gammaj b?etaj gamma?m betam]/(?2*A); if p?==1 ? ? ? %对于平?面应力问题形成弹性矩?阵D D=(?E/(1-NU^2)?)*[1 NU 0;?NU 1 0;0 0? (1-NU)/2]?; elseif p?==2 ? ? ? %对于平面应变问?题形成弹性矩阵D ? D=(E/(1?+NU)/(1-2*?NU))*[1-NU? NU 0;NU 1?-NU 0;0 0 ?(1-2*NU)/2?]; end y=t?*A*B*D*B;? ? ? %根据?虚功原理形成单元刚度?矩阵 2.Line?arTriangle?Assemble(K?,k,i,j,m) ?该函数将连接节点i和?节点j的线性三角形元?的单元刚度矩阵k集成?到整体刚度矩阵K。每?集成一个单元，该函数?都返回2n*2n的整?体刚度矩阵K。 fu?nction y=L?inearTrian?gleAssembl?e(K,k,i,j,?m) K(2*i-1?,2*i-1)=K(?2*i-1,2*i-?1)+k(1,1);? K(2*i-1,2?*i)=K(2*i-?1,2*i)+k(1?,2); K(2*i?-1,2*j-1)=?K(2*i-1,2*?j-1)+k(1,3?