层合板求ABD矩阵 武汉大学左翔n=5;t=0.1E-2 0.1E-2 0.1E-2 0.1E-2 0.1E-2;ang=0 pi/2 0 pi/2 0;E1=9.60E+10;E2=2.40E+10;v12=0.1;v21=0.4;G12=1.0E+10;T1=zeros(3,3,n);Q1=zeros(3,3,n);A=0;B=0;D=0;Q=E1/(1-v12*v21),v21*E2/(1-v12*v21),0;v21*E2/(1-v12*v21),E2/(1-v12*v21),0;0,0,G12;for i=1:n T1(:,:,i)=(cos(ang(i)2,(sin(ang(i)2,-2*sin(ang(i)*cos(ang(i); (sin(ang(i)2,(cos(ang(i)2,+2*sin(ang(i)*cos(ang(i); sin(ang(i)*cos(ang(i),-sin(ang(i)*cos(ang(i),(cos(ang(i)2-
