MATLAB有限元与应用剖析.pptVIP

3、实例计算分析应用 E=210E6; A=0.003; L=2; k1= QuadraticBarElementStiffness(E,A,L); k2= QuadraticBarElementStiffness(E,A,L); n=5; K = zeros(n,n); K =QuadraticBarAssemble(K,k1,1,3,2); K =QuadraticBarAssemble(K,k2,3,5,4); U=zeros(4,1); F=[5;-10;-7;10]; KK=K; K=K(2:n,2:n); U=K\F; U=[0;U]; F=KK*U; u1=[U(1);U(3);U(2)]; f1= QuadraticBarElementForces(k1,u1); sigma1=QuadraticBarElementStresses(k1, u1, A); u2=[U(3);U(5);U(4)]; f2=QuadraticBarElementForces(k2,u2); sigma2=QuadraticBarElementStresses(k2, u2, A); §3-3 二次杆元 1、基本方程 §3-4 平面桁架元 平面桁架元是既有局部坐标又有总体坐标二维有限元，用线性函数描述 每个平面桁架元有二个节点（node） 单刚矩阵为： 总刚矩