文章提供了一个基于MATLAB的编程练习，涉及单层板在特定应力状态下的应变计算。给定的应力值包括σ1=100MPa,σ2=-30MPa,τ12=10MPa，以及剪切模量G12=2.60×10^3。通过矩阵运算求解了正轴向的应变ε1,ε2及剪应变γ12。
一、单层板的强度和刚度 mathlab编程练习题（2）
2.在第1题的基础上，受有应力
σ
1
=
100
M
P
a
{\sigma _1} = 100{\rm{MPa}}
σ1=100MPa,
σ
2
=
−
30
M
P
a
{\sigma _2} = -30{\rm{MPa}}
σ2=−30MPa,
τ
12
=
10
M
P
a
{\tau _{12}} = 10{\rm{MPa}}
τ12=10MPa,
G
12
=
2.60
×
1
0
3
{{\rm{G}}_{12}} = 2.60 \times {10^3}
G12=2.60×103,求
ε
1
{\varepsilon _1}
ε1,
ε
2
{\varepsilon _2}
ε2,
γ
12
{\gamma _{12}}
γ12。
解：
相关公式如下:（仅供参考，如有错误欢迎指出）
单层板正轴向的应变-应力关系式可以写成如下的矩阵形式：
{
ε
1
}
=
[
S
]
{
σ
1
}
\left\{ {{\varepsilon _1}} \right\} = \left[ S \right]\left\{ {{\sigma _1}} \right\}
{ε1}=[S]{σ1}
在单层板中（应该是）：
ε
1
=
{
ε
1
ε
2
γ
12
}
\varepsilon _1=\begin{Bmatrix}{\varepsilon _1}\\{\varepsilon _2}\\{\gamma _{12}}\\\end{Bmatrix}
ε1=⎩
⎨
⎧ε1ε2γ12⎭
⎬
⎫,
σ
1
=
{
σ
1
σ
2
τ
12
}
\sigma_1=\begin{Bmatrix}{\sigma_1}\\{\sigma_2}\\{\tau_{12}}\\\end{Bmatrix}
σ1=⎩
⎨
⎧σ1σ2τ12⎭
⎬
⎫
根据上述公式使用matlab进行编程，程序代码如下：
%在第一题的基础上添加如下代码:
sigma1=100; sigma2=-30; tau12=10;
YingLi = [ sigma1;
sigma2;
tau12;];
%求应变
YingBian=S*YingLi
%将应变向量中的数值分别赋值给epsilon1,epsilon2,gamma12
pCell = num2cell(YingBian);
[epsilon1,epsilon2,gamma12] = pCell{:};
参考文献：
[1]王耀先.复合材料结构设计[M].北京: 化学工业出版社, 2001.09 第二章
相关文章：
一、复合材料单层板编程练习（1）
一、复合材料单层板编程练习（2）
一、复合材料单层板编程练习（3）
一、复合材料单层板编程练习（4）
一、复合材料单层板编程练习（5）
扫一扫
1116
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
