三、层合板的编程练习题(2)

三、层合板的强度和刚度 mathlab编程练习题(2)

2.在第1题的基础上，铺层层数$\mathrm{n}=3$层，$(0^{\circ }/90^{\circ }/0^{\circ })$，刚度系数$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})$与$\theta \in [0,90]$的关系。
解：

改变角度，使用matlab进行编程，程序代码如下：

%第三次编程课
%第二、三题
%2023-0404-1350
%2023-0413-2024-改进-几乎未改
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 ];

syms x;     %角度
%x=[0,90]
%theta=[0,x,0];
theta=[x,90,x];
%变0，不变90
h=0.125;    %单层厚度

n=3;        %层数

%z层高，是x、h、n的函数，
z=inline('h*(x-n/2)','x','h','n')

A=zeros(3,3);   B=zeros(3,3);   D=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_ = T \ Q * (T^-1)';
    %A=A+(z(i,h,n)-z(i-1,h,n))*Q_;
    B=B+(z(i,h,n)^2-z(i-1,h,n)^2)*Q_/2;
    D=D+(z(i,h,n)^3-z(i-1,h,n)^3)*Q_/3;
end
%输出B_ij与theta的关系
close
for i=1:1:n
    for j=i:1:n
        %fplot(A(i,j),[0 90])
        fplot(B(i,j),[0 90])
        %fplot(D(i,j),[0 90])
        hold on; 
    end
end

%legend('A11','A12','A16','A22','A26','A36');
legend('B11','B12','B16','B22','B26','B36');
%legend('D11','D12','D16','D22','D26','D36');
figure
%输出D_ij与theta的关系

for i=1:1:n
    for j=i:1:n
        %fplot(A(i,j),[0 90])
        %fplot(B(i,j),[0 90])
        fplot(D(i,j),[0 90])
        hold on; 
    end
end
%legend('A11','A12','A16','A22','A26','A36');
%legend('B11','B12','B16','B22','B26','B36');
legend('D11','D12','D16','D22','D26','D36');

计算结果

弯曲刚度矩阵$B_{ij}$如下：

耦合刚度系数$D_{ij}$如下：

theta=[x,90,x];

theta=[0,x,0];