《有限元方法与MATLAB程序设计》 课件 8 板的弯曲.pptxVIP

第8章板的弯曲;薄板小挠度弯曲基本理论;薄板小挠度弯曲基本理论;物理方程;弯矩;弯矩;板的内力;板的变形;应力;§8.1 四结点矩形薄板单元;§8.1.2 位移模式和形函数;§8.1.2 位移模式和形函数;形函数性质;位移的连续性;15;§8.1.3 应变和几何矩阵;§8.1.3 应变和几何矩阵;§8.1.4 应力和内力;§8.1.4 应力和内力;§8.1.4 刚度矩阵;§8.1.5 等效结点力;22;23;单元刚度矩阵 function ke=PlatRectStif(Em,Th,Pr,a,b) % 矩形板单元刚度矩阵;单元刚度矩阵 xj=x1(j);x0=x*xj;yj=y1(j);y0=y*yj; Nj=(1+x0)*(1+y0)*(2+x0+y0-x*x-y*y)/8; Njx=-b*yj*(1+x0)*(1+y0)*(1-y*y)/8; Njy= a*xj*(1+x0)*(1+y0)*(1-x*x)/8; N=[Nj Njx Njy]; % Njxx=diff(diff(N,x),x)/a^2; Njxy=diff(diff(N,x),y)/a/b; Njyy=diff(diff(N,y),y)/b^2; Bj=-[Njxx;Njyy;2*Njxy]; % ke(3*i-3+(1:3),3*j-3+(1:3))=a*b*double(int(int(Bi*D*Bj,x,-