二、复合材料层合板编程练习

ninepss

已于 2023-04-15 11:15:58 修改

阅读量658

点赞数

分类专栏：复合材料结构设计课程作业文章标签： matlab

于 2023-04-01 20:47:15 首次发布

本文链接：https://blog.csdn.net/qq_45144828/article/details/129902013

版权

复合材料结构设计课程作业专栏收录该内容

9 篇文章 1 订阅

订阅专栏

二、层合板的强度和刚度 mathlab编程练习题

1.已知碳/环氧三层层合板 $30^\circ/-30^\circ/30^\circ$ ，每层厚0.125mm，单层板 ${E_1} = 2.00 \times {10^5}MPa$ ， ${E_2} = 2.00 \times {10^4}MPa$ ， ${\mu _{21}} = 0.3$ ， ${G_{12}} = 2.00 \times {10^4}MPa$ 。
(1). 求各单层板 ${\overline Q _{ij}}$ 及层合板的刚度系数 ${\overline A _{ij}}$ 。
(2). 讨论铺层层数 $n = 3$ 层，拉伸刚度 ${A_{ij}}({A_{11}},{A_{12}},{A_{16}},{A_{22}},{A_{26}},{A_{66}})$ 与 $[\theta {\rm{ , }} - \theta {\rm{ , }}\theta ]{\rm{ , }}\theta \in [0{\rm{ , }}90]$ 的关系。
(3). 讨论铺层角度一定时(如题干所述)，铺层层数 $n = 1, 2, 3, 4, 5$ 与拉伸刚度 ${A_{ij}}({A_{11}},{A_{12}},{A_{16}},{A_{22}},{A_{26}},{A_{66}})$ 的关系。
题目理解：
第二问：可以解释为 $\theta = [x{\rm{ , }} - x{\rm{ , }}x]{\rm{ , }}x \in [0{\rm{ , }}90]$

(1). 求各单层板 ${\overline Q _{ij}}$ 及层合板的刚度系数 ${\overline A _{ij}}$ 。
解：
相关公式如下:（仅供参考，如有错误欢迎指出）
$\left[ {\overline {\rm{Q}} } \right] = {\left[ {{T_\sigma }} \right]^{ - 1}}\left[ Q \right]\left[ {{T_\varepsilon }} \right]$
设 $\cos \theta }$ ， $\sin \theta }$ ，则

$\begin{bmatrix} T_\varepsilon \end{bmatrix}= \begin{bmatrix} m^2 & n^2 & mn \\ n^2 & m^2 & -mn \\ -2mn & 2mn & m^2-n^2 \\ \end{bmatrix}$

$\begin{bmatrix} T_\sigma \end{bmatrix}= \begin{bmatrix} m^2 & n^2 & 2mn \\ n^2 & m^2 & -2mn \\ -mn & mn & m^2-n^2 \\ \end{bmatrix}$
对比上二式可得如下结论：
$\left[ {{T_\varepsilon }} \right] = {\left[ {{{\left[ {{T_\sigma }} \right]}^{ - 1}}} \right]^T}$
$\left[ {{T_\sigma }} \right] = {\left[ {{T_\varepsilon }} \right]^{ - 1}}$
${A_{ij}} = \sum\limits_{k = 1}^n {\overline Q _{ij}^{(k)}} {t_k}$
根据上述公式使用matlab进行编程，程序代码如下：

%复合材料层合板编程
clear;clc
E1=2*10^5;	E2=2*10^4;	G12=1*10^4;	nu21=0.3;
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 ];

n=3     %层数
theta=[30,-30,30];
h=0.125;        %单层厚度
A=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];
    fprintf("第%d层单层板的刚度矩阵为：\n",i)
    Q_ = T \ Q * (T^-1)'
    A=A+h*Q_;
end
disp("层合板刚度系数为：")
A

(2). 讨论铺层层数 $n = 3$ 层，拉伸刚度 ${A_{ij}}({A_{11}},{A_{12}},{A_{16}},{A_{22}},{A_{26}},{A_{66}})$ 与 $[\theta {\rm{ , }} - \theta {\rm{ , }}\theta ]{\rm{ , }}\theta \in [0{\rm{ , }}90]$ 的关系。
解：

%第二次编程课
%第二题
clear;clc;
E1=2*10^5;	E2=2*10^4;	G12=1*10^4;	nu21=0.3;
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;     %角度
theta=[x,-x,x];
h=0.125;    %单层厚度

p=3;        %层数，为什么不用n？因为被用了
A=zeros(3,3);
for i=1:1:p
    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+h*Q_;
end
%输出A_ij与theta的关系
for i=1:1:p
    for j=i:1:p
        fplot(A(i,j),[0 90])
        hold on; 
    end
end
legend('A11','A12','A16','A22','A26','A36');

(3). 讨论铺层角度一定时(如题干所述)，铺层层数 $n = 1, 2, 3, 4, 5$ 与拉伸刚度 ${A_{ij}}({A_{11}},{A_{12}},{A_{16}},{A_{22}},{A_{26}},{A_{66}})$ 的关系。
解：

%第二次编程课
%第三题
clear;clc;
E1=2*10^5;	E2=2*10^4;	G12=1*10^4;	nu21=0.3;
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 ];

%2-3
theta=[30,-30,30];
%单层厚度
h=0.125;
%层数
p=5;
A=cell(1,p);
for i=1:1:p
    theta=30*(-1)^(1+i);
    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}=h*Q_;
end
%输出图像
%注意：A是个1*p的矩阵，下面引用时没按照标准引用来
%注意：循环变量的命名！！！
close
for i=1:1:p
    %画点
    plot(i,A{i}(1,1),'-or');hold on;
    plot(i,A{i}(1,2),'-xg');hold on;
    plot(i,A{i}(1,3),'-*b');hold on;
    plot(i,A{i}(2,2),'-sc');hold on;
    plot(i,A{i}(2,3),'-hm');hold on;
    plot(i,A{i}(3,3),'-+k');hold on;
    %累加每一层的A_ij
    if i<p
        A{i+1}=A{i+1}+A{i};
    end
    %画线
    if(i>1)
        line([i-1,i],[A{i-1}(1,1),A{i}(1,1)],'color','r','LineStyle','-')
        line([i-1,i],[A{i-1}(1,2),A{i}(1,2)],'color','g','LineStyle','-')
        line([i-1,i],[A{i-1}(1,3),A{i}(1,3)],'color','b','LineStyle','-')
        line([i-1,i],[A{i-1}(2,2),A{i}(2,2)],'color','c','LineStyle','-')
        line([i-1,i],[A{i-1}(2,3),A{i}(2,3)],'color','m','LineStyle','-')
        line([i-1,i],[A{i-1}(3,3),A{i}(3,3)],'color','k','LineStyle','-')
    end
end
legend('A11','A12','A16','A22','A26','A66');