第2章_弹性力学基础及有限元法的基本原理1解答.ppt

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