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

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

3.铺层角度一定时，铺层层数$n=1,2,\cdots ,10$，与刚度系数$B_{ij}(B_{11},B_{12},B_{16},B_{22},B_{26},B_{66})$和$D_{ij}(D_{11},D_{12},D_{16},D_{22},D_{26},D_{66})$的关系。
解：

改变层数，使用matlab进行编程，程序代码如下

%第三次编程课
%2023-04-04-1501
%2023-04-13-2020-改进-未改
%第三题
clear;clc;
E1=180;	E2=10;	G12=7.2;	nu21=0.28;
nu12=nu21*E2/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 ];

%单层厚度
h=0.125;
%层数
n=10;
%z层高，是x、h、n的函数，
z=inline('h*(x-n/2)','x','h','n')

A=cell(1,n);    B=cell(1,n);    D=cell(1,n);
A(1,:)={zeros(3)};  B(1,:)={zeros(3)};  D(1,:)={zeros(3)};

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

%输出图像   B
%注意：循环变量的命名！！！
close all
for i=1:1:n
    %画点
    plot(i,B{i}(1,1),'-or');hold on;		plot(i,B{i}(1,2),'-xg');hold on;
    plot(i,B{i}(1,3),'-*b');hold on;		plot(i,B{i}(2,2),'-sc');hold on;
    plot(i,B{i}(2,3),'-hm');hold on;		plot(i,B{i}(3,3),'-+k');hold on;
    %累加每一层的A_ij
    if i<n
        %B{i+1}=B{i+1}+B{i};
    end
    %画线
    if(i>1)
        line([i-1,i],[B{i-1}(1,1),B{i}(1,1)],'color','r','LineStyle','-')
        line([i-1,i],[B{i-1}(1,2),B{i}(1,2)],'color','g','LineStyle','-')
        line([i-1,i],[B{i-1}(1,3),B{i}(1,3)],'color','b','LineStyle','-')
        line([i-1,i],[B{i-1}(2,2),B{i}(2,2)],'color','c','LineStyle','-')
        line([i-1,i],[B{i-1}(2,3),B{i}(2,3)],'color','m','LineStyle','-')
        line([i-1,i],[B{i-1}(3,3),B{i}(3,3)],'color','k','LineStyle','-')
    end
end
%legend('A11','A12','A16','A22','A26','A66');
legend('B11','B12','B16','B22','B26','B36');
%legend('D11','D12','D16','D22','D26','D36');
figure

%输出图像   D
%注意：循环变量的命名！！！

for i=1:1:n
    %画点
    plot(i,D{i}(1,1),'-or');hold on;		plot(i,D{i}(1,2),'-xg');hold on;
    plot(i,D{i}(1,3),'-*b');hold on;		plot(i,D{i}(2,2),'-sc');hold on;
    plot(i,D{i}(2,3),'-hm');hold on;		plot(i,D{i}(3,3),'-+k');hold on;
    %累加每一层的A_ij
    if i<n
        %D{i+1}=D{i+1}+D{i};
    end
    %画线
    if(i>1)
        line([i-1,i],[D{i-1}(1,1),D{i}(1,1)],'color','r','LineStyle','-')
        line([i-1,i],[D{i-1}(1,2),D{i}(1,2)],'color','g','LineStyle','-')
        line([i-1,i],[D{i-1}(1,3),D{i}(1,3)],'color','b','LineStyle','-')
        line([i-1,i],[D{i-1}(2,2),D{i}(2,2)],'color','c','LineStyle','-')
        line([i-1,i],[D{i-1}(2,3),D{i}(2,3)],'color','m','LineStyle','-')
        line([i-1,i],[D{i-1}(3,3),D{i}(3,3)],'color','k','LineStyle','-')
    end
end
%legend('A11','A12','A16','A22','A26','A66');
%legend('B11','B12','B16','B22','B26','B36');
legend('D11','D12','D16','D22','D26','D36');

计算结果

弯曲刚度矩阵$B_{ij}$随层数变化图：

耦合刚度系数$D_{ij}$随层数变化图：