计算力学课程设计

.5,0.5,0,1,0,ID); %(3) 建立整体刚度方程 KK= zeros( 18, 18 ); KK= Quad2D4Node_Assembly(KK,k1,4,1,2,5); KK= Quad2D4Node_Assembly(KK,k2,1,7,8,2); KK= Quad2D4Node_Assembly(KK,k3,2,8,9,3); KK= Quad2D4Node_Assembly(KK,k4,5,2,3,6); %(4) 边界条件的处理及刚度方程求解 k=KK(1:6,1:6); p=[0;-500;0;0;0;0]; u=k\p; %(5) 支反力的计算 U=[u;0;0;0;0;0;0;0;0;0;0;0;0]; P=KK*U; % (6 )各单元的应力计算 u1=[ U(7);U(8);U(1);U(2);U(3);U(4);U(9);U(10)]; stress1=Quad2D4Node_Stress(E,NU,1,1,0.5,1,0.5,0.5,1,0.5,u1,ID); u2=[ U(1);U(2);U(13);U(14);U(15);U(16);U(3);U(4)]; stress2=Quad2D4Node_Stress(E,NU,0.5,1,0,1,0,0.5,0.5,0.5,u2,ID); u3=[ U(3);U(4);U(15);U(16);U(17);U(18);U(5);U(6)]; stress3=Quad2D4Node_Stress(E,NU,0.5,0.5,0,0.5,0,0,0.5,0,u3,ID); u4=[ U(9);U(10);U(3);U(4);U(5);U(6);U(11);U(12)]; stress4=Quad2D4Node_Stress(E,NU,1,0.5,0.5,0.5,0.5,0,1,0,u4,ID);