马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
- /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(6)=-45
- ARG(7)=45
- ARG(8)=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(i)=-H/2+(i-1)*Tk
- z(N+2-i)=-z(i)
- *enddo
- *do,i,1,n !给退化刚度系数赋初值
- Q11(i)=VV*E1
- Q22(i)=VV*E2
- Q12(i)=V21*VV*E2
- Q66(i)=G12
- *enddo
- !*****************************************************************
- !***********************进入循环计算层合板强度********************
- !*****************************************************************
- *do,J,1,2*n+1
- *if,over,NE,0,then !判断层合板是否失效
- *do,i,1,n
- RM=COS(ARG(i))
- RN=SIN(ARG(i))
- RM2=RM*RM
- RM4=RM2*RM2
- RN2=RN*RN
- RN4=RN2*RN2
- RMN=RM*RN
- RMN2=RMN*RMN
- HQ11(i)=Q11(i)*RM4+2.0*(Q12(i)+2.0*Q66(i))*RMN2+Q22(i)*RN4
- HQ12(i)=(Q11(i)+Q22(i)-4.0*Q66(i))*RMN2+Q12(i)*(RM4+RN4)
- HQ16(i)=-RMN*RN2*Q22(i)+RM2*RMN*Q11(i)-RMN*(RM2-RN2)*(Q12(i)+2.0*Q66(i))
- HQ22(i)=Q11(i)*RN4+2.0*(Q12(i)+2.0*Q66(i))*RMN2+Q22(i)*RM4
- HQ26(i)=-RMN*RM2*Q22(i)+RMN*RN2*Q11(i)+RMN*(RM2-RN2)*(Q12(i)+2.0*Q66(i))
- HQ66(i)=(Q11(i)+Q22(i)-2*Q12(i))*RMN2+Q66(i)*(RM2-RN2)*(RM2-RN2)
- *enddo
- !*****************************************************************
- !*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
复制代码
|