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

ninepss

已于 2023-04-14 19:36:41 修改

阅读量175

点赞数

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

于 2023-03-27 20:56:20 首次发布

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

版权

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

9 篇文章 1 订阅

订阅专栏

文章提供了一道关于使用Matlab解决碳/环氧单层板在不同角度下的强度和刚度问题的编程练习。通过给定的材料属性和公式，计算S和Q矩阵，并求解不同角度下的应变和应力转换。程序涉及到矩阵的逆运算和转置，以及复合材料力学性质的计算。

摘要由CSDN通过智能技术生成

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

3.已知碳/环氧单层板的 ${{\rm{E}}_1} = 2.00 \times {10^5}MPa$ ， ${{\rm{E}}_2} = 2.00 \times {10^4}MPa$ ， ${\mu _{21}} = 0.3$ ， ${G_{12}} = 1.00 \times {10^4}MPa$ ，求 ${\theta = 30^\circ }$ , ${45^\circ}$ ,和 ${60^\circ}$ 时 $\overline{{{\rm{S}}_{ij}}}$ 和 $\overline{{{\rm{Q}}_{ij}}}$ 。
解：
相关公式如下:（仅供参考，如有错误欢迎指出）
$\left[ {\overline {\rm{Q}} } \right] = {\left[ {{T_\sigma }} \right]^{ - 1}}\left[ Q \right]\left[ {{T_\varepsilon }} \right]$
$\left[ {\overline S } \right] = {\left[ {{T_\varepsilon }} \right]^{ - 1}}\left[ S \right]\left[ {{T_\sigma }} \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] ^T= {\left[ {{T_\varepsilon }} \right]^{ - 1}}$
根据上述公式使用matlab进行编程，程序代码如下：

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
%计算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 ]
%第二题
sigma1=100;	sigma2=-30;	tau12=10;
YingLi = [  sigma1;
            sigma2;
            tau12;];
%求应变
YingBian=S*YingLi
%将应变向量中的数值分别赋值给epsilon1,epsilon2,gamma12
pCell = num2cell(YingBian);
[epsilon1,epsilon2,gamma12] = pCell{:};
%第三题
prompt = '输入30或45或60：';
theta=input(prompt)
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)'
S_ = T' * S * T