03杆梁结构.ppt

function y =Beam2D2Node_Assemble(KK,k,i,j) DOF(1)=3*i-2; DOF(2)=3*i-1; DOF(3)=3*i; DOF(4)=3*j-2; DOF(5)=3*j-1; DOF(6)=3*j; for n1=1:6 for n2=1:6 KK(DOF(n1),DOF(n2))=KK(DOF(n1),DOF(n2))+k(n1,n2); end end y = KK; 2. Beam2D2Node_Assemble(KK,k,i,j) 单元刚度矩阵和全体刚度矩阵的关系 单元刚度矩阵的行列编号 1 2 3 4 5 6 整体刚度矩阵的行列编号 3i-2 3i-1 3i 3j-2 3j-1 3j 节点的自由度数 function forces = Beam2D2Node_Forces(k,u) forces = k * u; 3. Beam2D2Node_Forces(k,u) 实践课题1 用matlab编制平面2D梁单元的计算程序。 用编制的matlab程序计算图示2D框架结构固定点C、D的弯矩，点A、B的弯矩及转角，载荷作用点的位移；并用材料力学理论求解，比较两者的结果。 E=20000kgf/mm2，A1=2000mm2，A2=1500mm2，I1=108mm4，I2=107mm4. P=4000N A B C D 5m 6m I1，A1 I2，A2 I2，A2 题1图 题2图 要求：1. 按照报告模板提交纸质文档及电子文档； 2. 提交截止日期：2011.10.13 function k =Element_Stiffness(E,I,A,L,alpha) k=[E*A/L,0,0,-E*A/L,0,0;0,12*E*I/(L^3),6*E*I/(L^2),0,-12*E*I/(L^3),6*E*I/(L^2);0,6*E*I/(L^2),4*E*I/L,0,-6*E*I/(L^2),2*E*I/L;-E*A/L,0,0,E*A/L,0,0;0,-12*E*I/(L^3),-6*E*I/(L^2),0,12*E*I/(L^3),-6*E*I/(L^2);0,6*E*I/(L^2),2*E*I/L,0,-6*E*I/(L^2),4*E*I/L]; x=alpha*pi/180; C=cos(x); S=sin(x); T=[C,S,0,0,0,0;-S,C,0,0,0,0;0,0,1,0,0,0;0,0,0,C,S,0;0,0,0,-S,C,0;0,0,0,0,0,1]; k=T*k*T; function forces = Forces(k,u) forces = k * u; function y =Global_Assemble(KK,k,i,j) DOF(1)=3*i-2; DOF(2)=3*i-1; DOF(3)=3*i; DOF(4)=3*j-2; DOF(5)=3*j-1; DOF(6)=3*j; for n1=1:6 for n2=1:6 KK(DOF(n1),DOF(n2))=KK(DOF(n1),DOF(n2))+k(n1,n2); end end y = KK; E=2e11; A1=2e-3; A2=1.5e-3; I1=1e-4; I2=1e-5; L2=5; L1=3; k1=Element_Stiffness(E,I1,A1,L1,0); k2=Element_Stiffness(E,I2,A2,L2,90); KK=zeros(15,15); KK=Global_Assemble(KK,k1,1,5); KK=Global_Assemble(KK,k1,5,2); KK=Global_Assemble(KK,k2,3,1); KK=Global_Assemble(KK,k2,4,2); k=KK([1,2,3,4,5,6,13,14,15],[1,2,3,4,5,6,13,14,15]); p=[0;0;0;0;0;0;0;-4000;0]; u=k\p; U=[u(1:6);0;0;0;0;0;0;u(7:9)]; P=KK*U; p3=k2*[0;0;0;U(1:3)]; p4=k2*[0;0;0;U(4:6)]; 实践课题1-（2）的输入步骤： C点:M3=-289.2NM D点:M4=289.2NM A点:M1=-580.8NMθ1=-0.000363rad B点:M2=580.8NMθ2=0.000363rad P点位移: u5= 0 v5= -0.0008298M 实践课题1-（2）的结果： §3.９ 单元形状函数矩阵与刚度矩阵的性质 形状函数