《有限元方法与MATLAB程序设计》 课件 6 空间问题（28页）.pptxVIP

第6章 空间单元;§6.1 空间四面体常应变单元;四面体单元的产生;平行六面体的体积;§6.1 空间四面体常应变单元;§6.1.1 位移模式和形函数;比较1-3行角标，满足角标轮换规律;§6.1.1 位移模式和形函数;形函数;§6.1.2 应变和几何矩阵;§6.1.2 应变和几何矩阵;§6.1.3 应力和应力矩阵;§6.1.4 刚度矩阵;§6.2 空间立方体单元;§6.2 空间立方体单元;§6.2.2 应变和几何矩阵;§ 6.2.2 应变和几何矩阵;§ 6.2.2 应变和几何矩阵 应力;§6.2.3 单元刚度矩阵;§6.2.4 程序设计;定义结构;自由度和结点力;刚度矩阵;刚度矩阵 syms x y z real% 定义符号变量 for i=1:8 xi=x0(i);yi=y0(i);zi=z0(i); Nx=xi*(1+yi*y)*(1+zi*z)/8/a; Ny=(1+xi*x)*yi*(1+zi*z)/8/b; Nz=(1+xi*x)*(1+yi*y)*zi/8/c; Bi=[Nx,0,0;0,Ny,0;0,0,Nz;0,Nz,Ny;Nz,0,Nx;Ny,Nx,0]; for j=1:8 xj=x0(j);yj=y0(j);zj=z0(j); Nx=xj*(1+yj*y)*(1+zj*z)/8/a; Ny=(1+xj*x)*yj*(1+zj*z)/8/b; Nz=(1+xj*x)*(1+yj*y)*zj/8/c; Bj=[Nx,0,0;0,Ny,0;0,0,Nz;0,Nz,Ny;Nz,0,Nx;Ny,Nx,0]; ke(3*i+(-2:0),3*j+(-2:0))=double(int(int(int(Bi*D*Bj,x,-1,1),y,-1,1),z,-1,1 end end;25;计算并输出应力;特征值问题;结 束