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

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

5.讨论工程弹性常数${\mathrm{E}}_x$ , $G_{xy}$ , ${\mu }_{yx}$ , 和${\eta }_{xy.x}$随$\theta$角变化的关系，以环氧/玻璃纤维材料为例，单层$\frac{E_1}{E_2}=3$ , $\frac{G_{xy}}{E_2}=0.25$时，算出偏轴无量纲工程弹性常数 , $\frac{E_x}{E_1}$ , $\frac{G_{xy}}{G_{12}}$ , ${\mu }_{yx}$和${\eta }_{xy.x}$随$\theta$角的变化曲线。
解：
相关公式如下:（仅供参考，如有错误欢迎指出）
设$n=\sin (\theta )$ , $m=\cos (\theta )$，则
$$\frac{1}{E_x}={\overline{S}}_{11}=\frac{1}{E_1}{m}^4+(\frac{1}{G_{12}}-\frac{2{\nu }_{21}}{E_1})m^2{n}^2+\frac{1}{E_2}{n}^4$$

$$\frac{1}{G_{xy}}={\overline{S}}_{66}=2(\frac{2}{E_1}+\frac{2}{E_2}+\frac{4{\nu }_1}{E_1}-\frac{1}{G_{12}})m^2{n}^2+\frac{1}{G_{12}}(m^2+n^2)$$

$${\nu }_x=-E_x{\overline{S}}_{12}=E_x\left[\frac{{\nu }_1}{E_1}(m^4+n^4)-(\frac{1}{E_1}+\frac{1}{E_2}-\frac{1}{G_{12}})m^2{n}^2\right]$$

$${\eta }_{xy,x}=E_x{\overline{S}}_{61}=E_x(A{m}^{3}n-Bm{n}^3)$$

其中：$A=\frac{2}{E_1}+\frac{2{\nu }_{21}}{E_1}-\frac{1}{G_{12}}$ , $B=\frac{2}{E_2}+\frac{2{\nu }_{21}}{E_1}-\frac{1}{G_{12}}$

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

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%第五题 代码
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Y1 = Ex / E1
% E1 / E2 = 3
% G12 / E2 = 0.5
% nu21 = 0.25
theta = 0:0.01:90;
C1 = 3;
C2 = 0.5;
nu21 = 0.25;

m = cosd(theta);
n = sind(theta);

% Ex/E1
Y1 = 1./(m.^4+(C1./C2-2*nu21)*n.^2.*m.^2+C1*n.^4);
%plot(theta,Y1)

%Gxy/G12
Y2 = 1./(n.^4+m.^4+(4*C2./C1.*(1+2*nu21)+4*C2-2).*n.^2.*m.^2);
%plot(theta,Y2)

%nuyx
Y3 = (nu21.*(n.^4+m.^4)-(1+C1-C1./C2).*n.^2.*m.^2)./(m.^4+(C1./C2-2.*nu21).*n.^2.*m.^2+C1.*n.^4);
%plot(theta,Y3)

%eta
Y4=((2+2*nu21-C1./C2).*n.*m.^3-(2*C1+2*nu21-C1./C2).*n.^3.*m)./(m.^4+(C1./C2-2*nu21).*n.^2.*m.^2+C1.*n.^4);
%plot(theta,Y4)
%输出图像
plot(theta,Y1,theta,Y2,theta,Y3,theta,Y4)
legend('E_{x}/E_{1}','G_{xy}/G_{12}','\mu_{yx}','\eta_{xy,x}')

结果：

参考文献：

[1]王耀先.复合材料结构设计[M].北京: 化学工业出版社, 2001.09 第二章

相关文章：

一、复合材料单层板编程练习（1）
一、复合材料单层板编程练习（2）
一、复合材料单层板编程练习（3）
一、复合材料单层板编程练习（4）
一、复合材料单层板编程练习（5）