一、复合材料单层板编程练习 (1)

ninepss

已于 2023-04-01 19:49:36 修改

阅读量286

点赞数

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

于 2023-03-27 12:56:10 首次发布

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

版权

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

9 篇文章 1 订阅

订阅专栏

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

1.已知碳/环氧单层板的 ${{\rm{E}}_1} = 2.00 \times {10^5}MPa$ ， ${{\rm{E}}_2} = 5.25 \times {10^3}MPa$ ， ${\mu _{21}} = 0.28$ ， ${G_{12}} = 2.60 \times {10^3}MPa$ ，求 ${S_{ij}}$ 和 ${Q_{ij}}$ 。
解：
相关公式如下:（仅供参考，如有错误欢迎指出）
弹性模量和泊松比之间的关系式：
${{{\nu _{12}}} \over {{E_2}}} = {{{\nu _{21}}} \over {{E_1}}}$
柔度矩阵：
$\begin{bmatrix} S \end{bmatrix}=\begin{bmatrix} S_{11} & S_{21} & 0 \\ S_{12} & S_{22} & 0 \\ 0 & 0 & S_{66}\\ \end{bmatrix}$
其中： $S_{11}=\frac{1}{E_1}$ , $S_{22}=\frac{1}{E_2}$ , $S_{66}=\frac{1}{G_{12}}$ , $S_{12}=-\frac{\nu_{12}}{E_2}$ , $S_{21}=-\frac{\nu_{21}}{E_1}$ , $S_{16}=S_{31}=S_{26}=S_{62}=0$
刚度矩阵：
$\begin{bmatrix} Q \end{bmatrix}=\begin{bmatrix} Q_{11} & Q_{21} & 0 \\ Q_{12} & Q_{22} & 0 \\ 0 & 0 & Q_{66}\\ \end{bmatrix}$
设 $M=(1-\nu_{21}\nu_{12})^{-1}$
则： $Q_{11}=M{E_1}$ , $Q_{22}=M{E_2}$ , $Q_{66}=G_{12}$ , $Q_{12}=M\nu_{12}{E_1}$ , $Q_{21}=M\nu_{21}{E_2}$ , $Q_{16}=Q_{31}=Q_{26}=S_{62}=0$

柔度矩阵和刚度矩阵为互逆关系， $Q = {S^{^{ - 1}}}$

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

clear
clc
E1=2*10^5;	E2=5.25*10^3;	G12=2.6*10^3;	nu21=0.28;
nu12=nu21*E2/E1;
%nu12/E2=nu21/E1
%计算S矩阵
S11=1/E1;	S22=1/E2;	S12=-nu12/E2;	S21=S12;	S66=1/G12;
S= [S11 , S12 , 0;
    S21 , S22 , 0;
    0 , 0 , S66 ]
%计算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 ]