QGforever的博客

原创 有限元（曾攀）B站杆单元构建及Matlab编程的代码

有限元（曾攀）杆单元构建及Matlab编程的代码function k=Bar2D2Node_Stiffness(E,A,x1,y1,x2,y2,alpha)%该函数计算单元的刚度矩阵%输入弹性模量E，横截面积A%输入第一个节点坐标（x1y1），第二个节点坐标（x2y2）角度 alpha（单位是度）%输出单元刚度矩阵k（4X4）%L=sqrt((x2-x1)*(x2-x1)+(y2-y1)*(y2-y1));x=alpha*pi/180;C=cos(x);S=sin(x);k=E*A/L