N = input('N=') ;
EL=180;
AL=-0.025E-6;
AT=27E-6;
ET=10;
RL=0.3;
RT=0.03;
GLT=5;
Q11=EL/(1-RL*RT);
Q22=ET/(1-RL*RT);
Q12=ET*RL/(1-RT*RL);
Q66=GLT;
Q=[Q11 Q12 0;Q12 Q22 0;0 0 Q66]
Q11L=zeros(1,N+1);
Q22L=zeros(1,N+1);
Q12L=zeros(1,N+1);
Q66L=zeros(1,N+1);
Q16L=zeros(1,N+1);
Q26L=zeros(1,N+1);
j=3;k=3;
T2=zeros(j,k,N+1);
T1=zeros(j,k,N+1);
for i=1:N+1;
c(i)=cos((i-1)*pi/180);
s(i)=sin((i-1)*pi/180);
Q11L(i)=Q11*(c(i)^4)+2*(Q12+2*Q66)*s(i)^2*c(i)^2+Q22*s(i)^4;
Q22L(i)=Q11*(s(i)^4)+2*(Q12+2*Q66)*s(i)^2*c(i)^2+Q22*c(i)^4;
Q12L(i)=Q12*(c(i)^4+s(i)^4)+(Q11+Q22-4*Q66)*(s(i)^2)*(c(i)^2);
Q66L(i)=Q66*(c(i)^4+s(i)^4)+(Q11+Q22-2*Q12-2*Q66)*s(i)^2*c(i)^2;
Q16L(i)=(Q11-Q12-2*Q66)*c(i)^3*s(i)+(-Q22+Q12+2*Q66)*c(i)*s(i)^3;
Q26L(i)=(Q11-Q12-2*Q66)*c(i)*s(i)^3+(-Q22+Q12+2*Q66)*c(i)^3*s(i);
Q21L(i)=Q12L(i);
Q61L(i)=Q16L(i);
Q62L(i)=Q26L(i);
T1(:,:,i)=[c(i)^2 s(i)^2 c(i)*s(i);s(i)^2 c(i)^2 -s(i)*c(i);-2*s(i)*c(i) 2*s(i)*c(i) c(i)^2-s(i)^2];
end
for i=1:N+1;
c(i)=cos((i-1)*pi/180);
s(i)=sin((i-1)*pi/180);
Q11L1(i)=Q11*(c(i)^4)+2*(Q12+2*Q66)*s(i)^2*c(i)^2+Q22*s(i)^4;
Q22L1(i)=Q11*(s(i)^4)+2*(Q12+2*Q66)*s(i)^2*c(i)^2+Q22*c(i)^4;
Q12L1(i)=Q12*(c(i)^4+s(i)^4)+(Q11+Q22-4*Q66)*(s(i)^2)*(c(i)^2);
Q66L1(i)=Q66*(c(i)^4+s(i)^4)+(Q11+Q22-2*Q12-2*Q66)*s(i)^2*c(i)^2;
Q16L1(i)=(Q11-Q12-2*Q66)*c(i)^3*s(i)+(-Q22+Q12+2*Q66)*c(i)*s(i)^3;
Q26L1(i)=(Q11-Q12-2*Q66)*c(i)*s(i)^3+(-Q22+Q12+2*Q66)*c(i)^3*s(i);
Q21L1(i)=Q12L1(i);
Q61L1(i)=Q16L1(i);
Q62L1(i)=Q26L1(i);
T2(:,:,i)=[c(i)^2 s(i)^2 c(i)*s(i);s(i)^2 c(i)^2 -s(i)*c(i);-2*s(i)*c(i) 2*s(i)*c(i) c(i)^2-s(i)^2];
end
j=3;k=3;
n=N+1;
A1=zeros(j,k,n);
A2=zeros(j,k,n);
for i=1:n;
A1(:,:,i)=[ Q11L(i) Q12L(i) Q16L(i);Q21L(i) Q22L(i) Q26L(i);Q61L(i) Q62L(i) Q66L(i)];
A2(:,:,i)=[ Q11L1(i) Q12L1(i) Q16L1(i);Q21L1(i) Q22L1(i) Q26L1(i);Q61L1(i) Q62L1(i) Q66L1(i)];
T=T1+T2;
A=A1+A2;
a(:,:,i)=inv(A(:,:,i));
u(i)=-a(1,2,i)/a(1,1,i);
end
j=3;k=1;ALF=zeros(j,k,n);W=zeros(j,k,1);
for i=1:n;
ALF(1,1,i)=AL*c(i)^2+AT*s(i)^2;
ALF(2,1,i)=AL*s(i)^2+AT*c(i)^2;
ALF(3,1,i)=2*s(i)*c(i)*(AL-AT);
alf(:,:,i)=-a(:,:,i)*T(:,:,i)*Q*ALF(:,:,i);
a1(i)=alf(1,1,i);
a2(i)=alf(2,1,i);
a3(i)=alf(3,1,i);
end
t=0:90;
subplot(2,2,1)
plot(t,u)
title('泊松比','FontWeight','bold')
xlabel('铺层角度')
ylabel('泊松比')
subplot(2,2,2)
plot(t,a1,'r');
title('x方向热膨胀系数','FontWeight','bold')
xlabel('铺层角度')
ylabel('膨胀系数')
subplot(2,2,3)
plot(t,a2,'g');
title('y方向热膨胀系数','FontWeight','bold')
xlabel('铺层角度')
ylabel('膨胀系数')
subplot(2,2,4)
plot(t,a3,'b');
title('剪切方向热膨胀系数','FontWeight','bold')
xlabel('铺层角度')
ylabel('膨胀系数')
没有合适的资源?快使用搜索试试~ 我知道了~
Matlab解复合材料层合板的问题.是matlab源代码.rar
共1个文件
m:1个
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 159 浏览量
2024-01-22
01:20:23
上传
评论
收藏 938B RAR 举报
温馨提示
毕业设计、竞赛参考资料,MATLAB学习资料,MATLAB项目源码
资源推荐
资源详情
资源评论
收起资源包目录
Matlab解复合材料层合板的问题.是matlab源代码.rar (1个子文件)
Matlab composite.m 3KB
共 1 条
- 1
资源评论
小正太浩二
- 粉丝: 199
- 资源: 5915
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功