经典层合板理论

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

下面是我用APDL写的程序。程序基于经典层合板理论。写的比较冗长,希望牛人能够给点意见,有没有哪位大侠用ANSYS做过层合板缺口件强度。希望交流交流程序写的很辛苦。
拿出来共享,希望能加分。

/prep7
*AFUN,DEG !参数定义为角度
OVER=1 !定义控制结束的参数,OVER=0,循环结束
N=8 !定义总层数
H=1 !定义总厚度
TK=H/N !定义每层的厚度
PX=100 !定义x方向外载荷
PY=0 !定义y方向外载荷
PXY=0 !定义xy方向剪切载荷
*dim,ARG,array,N !定义铺层的角度
ARG(1)=0
ARG(2)=45
ARG(3)=-45
ARG(4)=90
ARG(5)=90
ARG =-45
ARG(7)=45
ARG =0
Strainx=0 !定义中面x方向应变
strainy=0 !定义中面y方向应变
strainxy=0 !定义中面xy方向应变
*dim,strainSTEP,array,2*N !定义各个阶段的应变
*dim,pxstep,array,2*N !定义各个阶段的应力
*dim,strain1,array,N !定义各层x方向主应变
*DIM,Strain2,array,N !定义各层y方向主应变
*dim,Strain12,array,N !定义各层xy方向主应变
*dim,s1,array,N !定义各层X方向主应力
*dim,s2,array,N !定义各层Y方向主应力
*dim,s12,array,N !定义各层XY方向的剪应力
!****************************************************************
!*Define the strength parameter of lamia
!****************************************************************
XT=2250 !定义纤维方向拉伸强度
XC=1600 !定义纤维方向压缩强度
YT=34.5 !定义基体方向拉伸强度
YC=200 !定义基体方向压缩强度
s=110 !定义剪切强度
F1=1/XT-1/XC
F2=1/YT-1/YC
F11=1/(XT*XC)
F22=1/(YT*YC)
F66=1/(S**2)
F12=(-1/2)*sqrt(F11*F22)
*dim,a,array,N !定义系数a
*Dim,b,array,N !定义系数b
*dim,R,array,N !定义各层的强度比
Rmin=1000000 !定义最小强度比
RminN=0 !定义最小强度比层号
!****************************************************************
!*Define the parameter of lamia
!****************************************************************
E1=1.35E5
E2=0.8E4
V21=0.34
V12=E2*V21/E1
G12=4.5E3
VV=1/(1.0-V21*V12)
*dim,Q11,array,N !定义退化刚度系数
*dim,Q22,array,N
*dim,Q12,array,N
*dim,Q66,array,N
*dim,HQ11,array,N !定义偏轴刚度系数
*dim,HQ12,array,n
*dim,HQ16,array,n
*dim,HQ22,array,n
*dim,HQ26,array,n
*dim,HQ66,array,n
*dim,z,array,N+1
*do,i,1,nint(N/2+1-0.5)
z =-H/2+(i-1)*Tk
z(N+2-i)=-z
*enddo
*do,i,1,n !给退化刚度系数赋初值
Q11 =VV*E1
Q22 =VV*E2
Q12 =V21*VV*E2
Q66 =G12
*enddo
!*****************************************************************
!***********************进入循环计算层合板强度********************
!*****************************************************************
*do,J,1,2*n+1
*if,over,NE,0,then !判断层合板是否失效
*do,i,1,n
RM=COS(ARG )
RN=SIN(ARG )
RM2=RM*RM
RM4=RM2*RM2
RN2=RN*RN
RN4=RN2*RN2
RMN=RM*RN
RMN2=RMN*RMN
HQ11 =Q11 *RM4+2.0*(Q12 +2.0*Q66 )*RMN2+Q22 *RN4
HQ12 =(Q11 +Q22 -4.0*Q66 )*RMN2+Q12 *(RM4+RN4)
HQ16 =-RMN*RN2*Q22 +RM2*RMN*Q11 -RMN*(RM2-RN2)*(Q12 +2.0*Q66 )
HQ22 =Q11 *RN4+2.0*(Q12 +2.0*Q66 )*RMN2+Q22 *RM4
HQ26 =-RMN*RM2*Q22 +RMN*RN2*Q11 +RMN*(RM2-RN2)*(Q12 +2.0*Q66 )
HQ66 =(Q11 +Q22 -2*Q12 )*RMN2+Q66 *(RM2-RN2)*(RM2-RN2)
*enddo
!*************************************
|<< << < 1 2 3 > >> >>|
www.mapeng.net 马棚网
www.mapeng.net
文章作者:未知 | 文章来源:网络 | 责任编辑:admin | 发送至邮箱: | 加入收藏:
本文关键字:经典  层合板  理论
本文所属专题:ANSYS技术 
相关资讯
热点资讯
推荐资讯

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