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

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

1.已知三层层合板梁长200mm，宽10mm，单层$E_1=5.0\times 10^{4}MPa$，$E_2=1.0\times 10^{4}MPa$，${\mu }_{21}=0.40$，$G_{12}=2.0\times 10^{4}MPa$，铺设顺序为$0^{\circ }/90^{\circ }/0^{\circ }$，每层厚度为$1/2/1(mm)$,两端作用拉力$P=5000N$。求各层应力分布。
解
相关公式如下
{ N } = [ A ] { ε 0 } \{ N\} = [A]\{ {\varepsilon ^0}\} {N}=[A]{ε0}
{ δ x } = Q ‾ k { ε 0 } \{ {\delta ^x}\} = {\overline Q _k}\{ {\varepsilon ^0}\} {δx}=Q​k​{ε0}

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

clear;clc;
E1=5.0*10^4;	E2=1.0*10^4;	G12=2.0*10^3;	nu21=0.40;
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=[0,90,0];
%单层厚度B
%syms t

h=[1,2,1];

%层数
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

%外力
N=[5000;0;0];

sigma_0 = inv(A)*N;

sigma=cell(1,n);        sigma(1,:)={zeros(3,1)};

for i=1:1:n
    sigma{1,i}=Q_{1,i}*sigma_0;
end 
sigma{1,:}

计算结果

ans =
   1.0e+03 *

    2.0984
    0.1131
         0
         
ans =
  401.5837
 -113.1222
         0
         
ans =
   1.0e+03 *

    2.0984
    0.1131
         0

参考文献：

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