层合板求ABD矩阵(matlab).docx

《层合板求ABD矩阵(matlab).docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、层合板求ABD矩阵武汉大学左翔n=5;t=[0.1E-20.1E-20.1E-20.1E-20.1E-2];ang=[0pi/20pi/20];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];fo

2、ri=1:nT1(:,:,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-(sin(ang(i)))^2];Q1(:,:,i)=T1(:,:,i)*Q*T1(:,:,i)';endm

3、=sum(t)/2;M=zeros(1,n+1);fori=1:nM(i+1)=M(i)+t(i);endz=M-m;fori=2:n+1A=A+Q1(:,:,i-1)*(z(i)-z(i-1));B=B+(1/2)*Q1(:,:,i-1)*(z(i)^2-z(i-1)^2);D=D+(1/3)*Q1(:,:,i-1)*(z(i)^3-z(i-1)^3);end