三、层合板的编程练习题（1）

三、层合板的强度和刚度 mathlab编程练习题（1）

1.已知[T300/5208]材料$E_1=180GPa$，$E_2=10GPa$，${\mu }_{21}=0.28$，$G_{12}=7.2GPa，每层厚0.125mm$。求$90^{\circ }/0^{\circ }/0^{\circ }/90^{\circ }$与$0^{\circ }/90^{\circ }/0^{\circ }/90^{\circ }/0^{\circ }$时的刚度系数$A_{ij}$，$B_{ij}$和$D_{ij}$。

解：
相关公式如下:（仅供参考，如有错误欢迎指出）
$$A_{ij}=\sum _{k=1}^n{\overline{Q}}_{ij}^{(k)}dz=\sum _{k=1}^n{\overline{Q}}_{ij}^{(k)}(z_k-z_{k-1}),(i,j=1,2,6)$$
$$B_{ij}=\sum _{k=1}^n{\overline{Q}}_{ij}^{(k)}zdz=\frac{1}{2}\sum _{k=1}^n{\overline{Q}}_{ij}^{(k)}(z_k^2-z_{k-1}^2),(i,j=1,2,6)$$
$$D_{ij}=\sum _{k=1}^n{\overline{Q}}_{ij}^{(k)}{z}^{2}dz=\frac{1}{3}\sum _{k=1}^n{\overline{Q}}_{ij}^{(k)}(z_k^3-z_{k-1}^3),(i,j=1,2,6)$$

根据上述公式使用matlab进行编程，程序代码如下：

%2023-04-04-1332
%2023-04-13-1910-改进
%第一题
%单位换算注意
clear;clc;
E1=180;	E2=10;	G12=7.2;	nu21=0.28;
nu12=nu21*E2/E1;
%nu12/E2=nu21/E1

%计算Q矩阵
M=(1-nu12*nu21)^(-1);
Q11=M*E1;	Q22=M*E2;	Q12=M*nu12*E1;	Q21=Q12;	Q66=G12;
Q= [Q11 , Q12 , 0;
    Q21 , Q22 , 0;
    0 , 0 ,Q66 ];
%角度矩阵
%theta=[90,0,0,90];
theta=[0,90,0,90,0];

%t单层厚度
t=0.125;

%厚度矩阵
%h=[1*t,1*t,1*t,1*t];
h=[1*t,1*t,1*t,1*t,1*t];

n=length(theta);
z(1)=-sum(h)/2;
for i=1:1:n
    z(i+1) = z(i)+h(i);
end

A=zeros(3,3);   B=zeros(3,3);   D=zeros(3,3);
Q_=cell(1,n);   Q_(1,:)={zeros(3,3)};
for i=1:1:n
    m=cosd(theta(i));  n_=sind(theta(i));
    T=[  m^2 n_^2 2*m*n_;
	     n_^2 m^2 -2*m*n_;
	     -m*n_ m*n_ m^2-n_^2];
    Q_{1,i} = T \ Q * (T^-1)';
    A=A+(z(i+1)-z(i))*Q_{1,i};
    B=B+(z(i+1)^2-z(i)^2)*Q_{1,i}/2;
    D=D+(z(i+1)^3-z(i)^3)*Q_{1,i}/3;
end
theta,A,B,D

计算结果

theta =
    90     0     0    90
    
A =
   47.7078    1.4061         0
    1.4061   47.7078         0
         0         0    3.6000
         
B =
   1.0e-15 *

   -0.1110         0         0
         0         0         0
         0         0         0
%理论计算B应该为0矩阵，此处为极小值，猜测在拼接每一层的B矩阵中存在大数减大数或小树减小数的情况
D =
    0.3269    0.0293         0
    0.0293    1.6609         0
         0         0    0.0750

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

theta =
     0    90     0    90     0
     
A =
   70.3062    1.7577         0
    1.7577   48.9633         0
         0         0    4.5000
         
B =
   1.0e-15 *

         0         0         0
         0    0.1665         0
         0         0   -0.0278
%理论计算B应该为0矩阵，此处为极小值，猜测在拼接每一层的B矩阵中存在大数减大数或小树减小数的情况
D =
    2.9556    0.0572         0
    0.0572    0.9269         0
         0         0    0.1465

参考文献：

[1]王耀先.复合材料结构设计[M].北京: 化学工业出版社, 2001.09 第四章