第2章_弹性力学基础及有限元法的基本原理1最新.pptVIP

2.5 约束处理 1）边界位移为零 2）边界位移为已知量 2.6 求解线性方程组 2.7 计算其它物理量 2.8 计算结果处理 2.9 结果显示、打印、分析 %----------------------------------------------------------- function k=Triangle2D3Node_Stiffness(E,NU,t,xi,yi,xj,yj,xm,ym,ID) %该函数计算单元的刚度矩阵 %输入弹性模量E，泊松比NU，厚度t %输入三个节点i、j、m的坐标xi,yi,xj,yj,xm,ym %输入平面问题性质指示参数ID(1为平面应力，2为平面应变) %输出单元刚度矩阵k(6X6) %--------------------------------------------------------------- A = (xi*(yj-ym) + xj*(ym-yi) + xm*(yi-yj))/2; betai = yj-ym; betaj = ym-yi; betam = yi-yj; gammai = xm-xj; gammaj = xi-xm; gammam = xj-xi; B = [betai 0 betaj 0 betam 0 ; 0 gammai 0 gammaj 0 gammam ; gammai betai gammaj betaj gammam betam]/(2*A); 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]; end k= t*A*B*D*B; （1） Triangle2D3Node_Stiffness 第三节 平面3节点三角形单元及MATLAB编程 %----------------------------------------------------------- function z = Triangle2D3Node_Assembly(KK,k,i,j,m) %该函数进行单元刚度矩阵的组装 %输入单元刚度矩阵k %输入单元的节点编号I、j、m %输出整体刚度矩阵KK %--------------------------------------------------------------- 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; for n1=1:6 for n2=1:6 KK(DOF(n1),DOF(n2))= KK(DOF(n1),DOF(n2))+k(n1,n2); end end z=KK; （2）Triangle2D3Node_Assembly %----------------------------------------------------------- function stress=Triangle2D3Node_Stress(E,NU,xi,yi,xj,yj,xm,ym,u,ID) %该函数计算单元的应力 %输入弹性模量E，泊松比NU，厚度t %输入三个节点i、j、m的坐标xi,yi,xj,yj,xm,ym %输入平面问题性质指示参数ID(1为平面应力，2为平面应变)，单元的位移列阵u(6X1) %输出单元的应力stress(3X1)，由于它为常应力单元，则单元的应力分量为Sx,Sy,Sz %--------------------------------------------------------------- A = (xi*(yj-ym) + xj*(ym-yi) + xm*(yi-yj))/2; betai = yj-ym; betaj = ym-yi; betam = yi-yj; gammai = xm-xj; gammaj = xi-xm; gammam = xj-xi; B = [betai 0 betaj 0 betam 0 ; 0 gammai 0 gammaj 0 gammam ; gammai betai gammaj betaj gammam betam]/(2*A); if ID == 1 D = (E/(1-NU*NU))*[1 NU 0 ; NU 1 0 ; 0 0 (1-NU)/2]; elseif ID == 2