经典层合板理论

发布日期:[08-12-04 15:12:59] 浏览人次:[]

****************************
!*Dedine Z Coordinate of Each Lamia
!*****************************************************************
A11=0
A12=0
A22=0
A26=0
A66=0
A16=0
B11=0
B12=0
B22=0
B26=0
B66=0
B16=0
D11=0
D12=0
D22=0
D26=0
D66=0
D16=0
!*****************************************************************
!* Define The Stiffness of Laminate
!*****************************************************************
*Do,i,1,N !积分求解层合板的刚度系数
A11=A11+HQ11 *(z(i+1)-z )
A12=A12+HQ12 *(z(i+1)-z )
A22=A22+HQ22 *(z(i+1)-z )
A26=A26+HQ26 *(z(i+1)-z )
A16=A16+HQ16 *(z(i+1)-z )
A66=A66+HQ66 *(z(i+1)-z )
B11=B11+0.5*HQ11 *(Z(i+1)**2-z **2)
B12=B12+0.5*HQ12 *(Z(i+1)**2-z **2)
B22=B22+0.5*HQ22 *(Z(i+1)**2-z **2)
B26=B26+0.5*HQ26 *(Z(i+1)**2-z **2)
B16=B16+0.5*HQ16 *(Z(i+1)**2-z **2)
B66=B66+0.5*HQ66 *(Z(i+1)**2-z **2)
D11=D11+1/3*HQ11 *(z(i+1)**3-z **3)
D12=D12+1/3*HQ12 *(z(i+1)**3-z **3)
D22=D22+1/3*HQ22 *(z(i+1)**3-z **3)
D16=D16+1/3*HQ16 *(z(i+1)**3-z **3)
D26=D26+1/3*HQ11 *(z(i+1)**3-z **3)
D66=D11+1/3*HQ11 *(z(i+1)**3-z **3)
*enddo
!*****************************************************************
!* Define The roughness of Laminate
!*****************************************************************
DA=(A11*A22*A66+2*A12*A26*A16-A16**2*A22) !求解层合板的柔度系数
DA=DA-A11*A26**2-A12**2*A66
HA11=(A22*A66-A26**2)/DA
HA12=(A16*A26-A12*A66)/DA
HA22=(A11*A66-A16**2)/DA
HA66=(A11*A22-A12**2)/DA
HA16=(A12*A26-A12*A16)/DA
HA26=(A12*A16-A11*A26)/DA
!*****************************************************************
!* Caculate The Strain and stress of each Lamina
!*****************************************************************
STRAINX=HA11*PX+HA12*PY+HA16*PXY !计算层合板的应变
STRAINY=HA12*PX+HA22*PY+HA26*PXY
STRAINXY=HA16*PX+HA26*PY+HA66*PXY

strainstep(j)=STRAINX !记录下各个阶段的应变
pxstep(j)=PX !记录下各个阶段的应力

*Do,i,1,N !计算各层的主应变
RM=COS(ARG )
RN=SIN(ARG )
RMN=RM*RN
RM2=RM*RM
RN2=RN*RN
strain1 =RM2*STRAINX+RN2*STRAINY+(RMN)*STRAINXY
Strain2 =RN2*STRAINX+RM2*STRAINY+(-RMN)*STRAINXY
STrain12 =-2*RMN*STRAINX+2*RMN*STRAINY+(RM2-RN2)*STRAINXY
S1 =strain1 *Q11 +strain2 *Q12 !计算各层的主应力
S2 =STrain1 *Q12 +STRAIN2 *Q22
s12 =strain12 *Q66
a =F11*(s1 **2)+F22*(S2 **2)
a =a +2*F12*s1 *s2 +F66*(S12 **2)
b =F1*S1 +F2*S2
R =-b /(2*a )+sqrt(b **2+4*a )/(2*a )
*enddo

*Vscfun,Rmin,min,R(1) !求出最小强度比
*Vscfun,RminN,lmin,R(1) !求出最小强度所在的层号

*if,s1(RminN)*(Rmin),gt,0,then !判断破坏层的失效模式
*if,s1(RminN)*(Rmin),lt,xt,then !并进行刚度退化
Q22(RminN)=Q22(RminN)*(10e-14)
Q12(RminN)=Q12(RminN)*(10e-14)
Q66(RminN)=Q66(RminN)*(10e-14)
*else
Q22(RminN)=Q22(RminN)*(10e-14)
Q12(RminN)=Q12(RminN)*(10e-14)
Q66(RminN)=Q66(RminN)*(10e-14)
Q11(RminN)=Q11(RminN)*(10e-14)
*ENDIF

*ELSEIF,s1(RminN)*(Rmin),le,0,then
*if,(-1)*s1(RminN)*(Rmin),lt,xc,then
Q22(RminN)=Q2
|<< << < 1 2 3 > >> >>|
www.mapeng.net 马棚网
www.mapeng.net
文章作者:未知 | 文章来源:网络 | 责任编辑:admin | 发送至邮箱: | 加入收藏:
本文关键字:经典  层合板  理论
本文所属专题:ANSYS技术 
相关资讯
热点资讯
推荐资讯

关于我们 | 站点导航 | 使用帮助 | 友情链接 | 广告服务 | 免责声明 | 新手上路
设为首页 | 加入收藏 | 在线留言 | 马棚网QQ群:{92562572}{102901272}{333259257} | 交流QQ: 客户服务 客户服务 客户服务