发一篇我用ANSYS-APDL做的复合材料圆柱层合板渐进失效的国际杂志论文和程序

kolarxiao

楼主，您好！您这里的渐进失效和疲劳失效是什么关系啊？

yinjvictor

这个文章怎么都是你们自己在引用，没有什么外国人引用啊！

limulin1

sxyuyifan 发表于 2010-2-4 08:26
% l8 c1 ^) I! k! G1 @% u楼主真强，是浙大化机所的啊。

& O/ Z, y; b8 d; @+ Z我想也是，看来你也是咯

dragon_w

我模仿您的apdl写了段，可还是结果无法提取，请教！
4 [% o: B" l3 \& yfini
/clear
+ o% z3 p- r3 Q7 `4 ^/prep7
' b8 A4 z* l( V+ S0 Set,1,46
2 b' ^- u; y2 }* \4 j: ^% Y4 pkeyopt,1,8,1
2 R8 k2 E  ~. W9 ~& a5 I& QR,1,4
RMODIF,1,13,1,0, 0.125,1,45,0.125
RMODIF,1,19,1,-45,0.125,1,90,0.125
. R4 u; G7 v4 j( a  u6 D  [MP,EX,1,114000
MP,EY,1,8610
, G, Q1 r$ o# }: F! N4 t% D9 CMP,EZ,1,8610
MP,PRXY,1,0.3
MP,PRYZ,1,0.45
' L9 s$ ]- Z6 NMP,PRXZ,1,0.3
/ R" N" g' `. o1 r$ a6 T2 qMP,GXY,1,4160
6 K3 u; G/ P/ BMP,GYZ,1,3000
MP,GXZ,1,4160
. }+ q% Y, H% u0 q( T" N+ r
4 o7 M* h$ f5 x8 G9 x  Q4 TMP,EX,2,114000
+ n# V, |6 r8 H8 D* ?MP,EY,2,8610*0.3
MP,EZ,2,8610
MP,PRXY,2,0.3
MP,PRXZ,2,0.3
MP,PRYZ,2,0.45
MP,GXY,2,4160*0.3
! U0 T% P9 g, B. R& z6 ]# x( D# `- dMP,GXZ,2,4160
1 U( r1 r' a: q; VMP,GYZ,2,3000*0.3

1 v) n8 j, ]+ S- DMP,EX,3,114000
9 P' U4 r3 W& V7 R7 gMP,EY,3,8610*0.4
MP,EZ,3,8610
# n9 l8 {$ R( r/ hMP,PRXY,3,0.3
9 q8 W% T$ e3 A: E* j2 t* \. wMP,PRXZ,3,0.3
: c9 @5 u1 u% z6 g) z9 uMP,PRYZ,3,0.45
MP,GXY,3,4160*0.4
2 _+ s  e% o" Z; \MP,GXZ,3,4160
1 u( @2 p0 m  }/ s4 lMP,GYZ,3,3000*0.4

MP,EX,4,114000*0.07
9 v  D; D; u% C$ C6 WMP,EY,4,8610
  @) I! n- S6 p7 N! y1 z- A8 p: CMP,EZ,4,8610
MP,PRXY,4,0.3
6 @5 i+ w8 _8 d  f" `MP,PRXZ,4,0.3
+ r0 D+ r, j) E/ ?1 fMP,PRYZ,4,0.45
MP,GXY,4,4160*0.07
MP,GXZ,4,4160*0.07
, N  A3 _  L1 N# `4 VMP,GYZ,4,3000
! ~5 F0 d7 T3 x
MP,EX,5,114000*0.14
* f% A8 G9 ?7 ^1 i$ z- J5 HMP,EY,5,8610
MP,EZ,5,8610
MP,PRXY,5,0.3
1 r1 ^% @* H+ e& Y. dMP,PRXZ,5,0.3
) \, M7 K; A9 h: \3 W9 gMP,PRYZ,5,0.45
MP,GXY,5,4160*0.07
% B/ O4 T" W' V$ X; O3 dMP,GXZ,5,4160*0.07
MP,GYZ,5,3000

MP,EX,6,114000
MP,EY,6,8610*0.3
- c5 j( @/ Q6 D3 m$ x: kMP,EZ,6,8610*0.3
MP,PRXY,6,0.3
MP,PRXZ,6,0.3
MP,PRYZ,6,0.45
MP,GXY,6,4160
MP,GXZ,6,4160*0.3
MP,GYZ,6,3000*0.3

& e7 R. A: w$ \1 mMP,EX,7,114000
MP,EY,7,8610*0.4
MP,EZ,7,8610*0.4
MP,PRXY,7,0.3
& @2 ]7 p9 J- YMP,PRXZ,7,0.3
5 J( L1 Z7 k: Y/ M2 lMP,PRYZ,7,0.45
; E% q7 g' T& z1 C0 v- z) _: M2 QMP,GXY,7,4160
MP,GXZ,7,4160*0.4
MP,GYZ,7,3000*0.4

' [0 h) J0 L5 k7 a  \! LMP,EX,8,114000*0.07
: e* S' o7 `5 M( YMP,EY,8,8610*0.3
MP,EZ,8,8610*0.3
MP,PRXY,8,0.3
' N! ?- I' k8 ~5 ?8 ~) C) }MP,PRXZ,8,0.3
MP,PRYZ,8,0.45
MP,GXY,8,4160*0.07
MP,GXZ,8,4160*0.07
' ?! b0 H6 ?# W2 ], n' ?MP,GYZ,8,3000*0.3
+ C  E1 w+ k' I3 I! m7 z
!mp,ex,2,210000
!mp,prxy,2,0.3
% `6 Q! k$ T+ F1 q& ?) f6 y* aXc=914.3
8 m' r7 \+ f8 I/ CXt=1229.4
Yc=147.2
: A. j, q+ N% j0 W& Q$ I3 iYt=45.6
6 q( }/ B- e7 sZc=147.2
/ V1 V+ p) }3 w  W1 c+ HZt=45.6
. m) h  J  u" ]& LS12=33.0
S13=33.0
( ^# p% S' h4 |# H5 u2 \& d0 U4 HS23=25.0
*dim,ST11,array,800
*dim,ST22,array,800
*dim,ST33,array,800
*dim,ST12,array,800
' p2 ?/ Z4 k& }: M# p- P% y6 p*dim,ST13,array,800
*dim,ST23,array,800
, d0 m" k! f5 x3 d8 U9 D*dim,DFT,array,800  
# A5 A0 G' S6 [, B2 ^* X& S*dim,DFC,array,800
*dim,DMT,array,800
+ b# B2 m( U6 U*dim,DMC,array,800
! T7 T" d' [% D) y*dim,DFM,array,800
*dim,DT,array,800
*dim,DC,array,800
/ ^+ p7 N3 y& T9 W2 F*dim,erci,array,800
/ I% {$ `5 F! X! S4 A*dim,HASHINDFT,array,800
9 J6 W$ h& x2 }4 t1 R*dim,HASHINDFC,array,800      
: K) S/ a( M! t& ~*dim,HASHINDMT,array,800
7 i- n3 b  T6 M: Y*dim,HASHINDMC,array,800
$ n3 P# L# ]1 j6 q/ B# [; L*dim,HASHINDFM,array,800
*dim,HASHINDT,array,800
/ O1 \1 g& x& g0 N" H2 q. Z/ v*dim,HASHINDC,array,800
- U0 ?1 W# ?( m4 V$ L9 M
block,0,100,0,50,0,0.5
' W! O4 ~. c1 E8 BLSEL,S,LENGTH,,100 !选择所有长度为90mm的线段(长度方向上的线段)
LESIZE,all, , ,20, , , , ,1
* q" q5 Q( z1 Alsel,s,length,,50 !选择所有长度为25mm的线段(宽度方向上的线段)
lesize,all, , ,10,, , , ,1
LSEL,S,LENGTH,,100
lsel,a,length,,50
lsel,inve
# k0 \  Y' F4 F1 E! W% olesize,all, , ,4
alls
vatt,1,1,1
; P0 U9 M# Q: e& l) w5 }% ~8 m! evmesh,all
# \+ Z1 t- p5 J) u* R

+ F7 R7 W) i. d4 r/solu
asel,s,loc,x,0
da,all,all
alls
$ Z% s8 S8 G- [/ Nasel,s,loc,x,100
sfa,all,1,press,-2
alls
: S8 X) x; Q7 _6 N5 X$ r  ~7 {9 ?!nsubst,5
6 @% F- e! M* K3 d9 Q- T( ?+ g!cnvtol,f,,0.1,,0.1
!rescontrol,define,all,1,1
: a: m, _! l# c8 j( J5 z3 |solve
: E5 F# t, F6 L# R) \* @4 B3 g!save,1,db
) {  T- K+ d  a8 L" |$ D- l
% j; F; R& C' A*do,i,2,10
5 `3 x: S% x0 G% }. q! S6 A; [  /post1
' a  A( a! W: G& S* N   etable,sx,s,x
2 X& X7 K; v! z* U   etable,sy,s,y
7 X8 _9 ~) D  c( @( t) j   etable,sz,s,z
   etable,sxy,S,xy
   etable,syz,s,yz
) j0 N' s5 ^' K( d* v3 C% Y. B; E   etable,sxz,S,xz
   
   *do,j,1,800
      *get,ST11(j),elem,j,etab,sx
      *get,ST22(j),elem,j,etab,sy
      *get,ST33(j),elem,j,etab,sz
      *get,ST12(j),elem,j,etab,sxy
      *get,ST13(j),elem,j,etab,sxz
      *get,ST23(j),elem,j,etab,syz                  
8 v1 h9 r3 E/ l# w  ?! B   *enddo
   *do,j,1,800
/ l. n8 G/ e0 S9 n, i$ Y      HASHINDFT(j) = (ST11(j)/Xt)**2+(ST12(j)/S12)**2+(ST13(j)/S13)**2
8 }9 N$ g5 C8 f7 H; u      HASHINDFC(j) = (ST11(j)/Xc)**2
0 R6 e! b( q0 X7 p8 F" ~) P      HASHINDMT(j) = (ST22(j)/Yt)**2+(ST12(j)/S12)**2+(ST23(j)/S23)**2
      HASHINDMC(j) = (ST22(j)/Yc)**2+(ST12(j)/S12)**2+(ST23(j)/S23)**2
      HASHINDFM(j) = (ST11(j)/Xc)**2+(ST12(j)/S12)**2+(ST13(j)/S13)**2
      HASHINDT(j) = (ST33(j)/Zt)**2+(ST13(j)/S13)**2+(ST23(j)/S23)**2
      HASHINDC(j) = (ST33(j)/Zc)**2+(ST13(j)/S13)**2+(ST23(j)/S23)**2              
, Z3 F1 R' p* t& N- L: U( i   *enddo   
parsav,all,parameter,txt

/ C( C* r; p1 E  /solu
! h  B8 ^0 I/ V. n6 r* j$ P+ |    antype,,rest
; Y; L7 j, E: R: N# Qparres,new,parameter,txt
!失效类型
!纤维失效
*do,j,1,800
      *IF,ST11(j),ge,0,and,HASHINDFT(j),ge,1,THEN
               DFT(j)=1+DFT(j)
$ k( ^( Q+ O9 X" i! V" k
            *ELSEif,ST11(j),lt,0,and,HASHINDFC(j),ge,1,then
) @+ J" }; {! l/ o( i               DFC(j)=1+DFC(j)

! e! F% G7 E/ }' F/ ^7 W4 E            *ELSE
4 A3 u/ Z2 K) Z2 H+ P7 Z6 J                DFC(j)=0+DFC(j)
: e4 @1 o5 ~7 J) X# {  ]2 e: z; u                DFT(j)=0+DFT(j)
           
, q7 I, v1 o4 w: \2 n      *ENDIF
2 d, c9 K; y, u% `+ I% a" J*enddo
0 X& j7 h4 V2 o3 [
!基体失效
*do,j,1,800
) l7 o. F5 [9 }) Z: t- x      *IF,ST22(j),ge,0,and,HASHINDMT(j),ge,1,THEN         
               DMT(j)=1+DMT(j)
         
          *ELSEIF,ST22(j),LT,0,and,HASHINDMC(j),ge,1,THEN
                DMC(j)=1+DMC(j)
8 ]7 q; e, G5 m
          *ELSE
7 t2 P: `2 S7 `3 o$ f# j0 G                DMC(j)=0+DMC(j)
                DMT(j)=0+DMT(j)
& |( k8 A' g3 N, ^         
      *ENDIF
7 q' R1 }# r8 v5 D0 @) i*enddo
3 Q+ m2 I  }$ @' l4 u+ l!FIBER-MATRIX FAILURE
*do,j,1,800
2 [1 P( y8 A* q) Q  [( U1 Y; w4 w5 e      *IF,ST11(j),ge,0,AND,HASHINDFM(j),ge,1,THEN
              DFM(j)=1+DFM(j)

            *ELSE
              DFM(j)=0+DFM(j)
            
, y& m) d" M$ R$ I! |2 z4 b# |/ J      *ENDIF
  @' u! }  B" d1 m! v* |*enddo
7 {+ W0 |0 V& R" B5 p7 Y" d
!Tensile-Compression FAILURE
*do,j,1,800
) m) |* Y* w8 R2 H  u: c      *IF,ST33(j),ge,0,AND,HASHINDT(j),ge,1,THEN
& M% P7 O! c/ |* i$ ?& g; W1 f* l2 {               DT(j)=1+DT(j)
( G' x' W" G9 `( W
# ]2 {: d& l7 r: ?7 f8 `          *ELSEIF,ST33(j),LT,0,AND,HASHINDC(j),ge,1,THEN
2 J, ]1 p) _, w, {3 S                DC(j)=1+DC(j)

; n/ v: [  t* B; ^9 D          *ELSE
                DC(j)=0+DC(j)
                DT(j)=0+DT(j)
9 C; C$ t9 y% f4 k0 Y         
      *ENDIF
' f& W- e4 Q1 S0 f  J; [*enddo
* r/ w% ^1 U/ v: n% F*do,j,1,800
    erci(j)=DFT(j)+DFC(j)+DMT(j)+DMC(j)+DFM(j)+DT(j)+DC(j)
*ENDDO
/prep7
!更新单元材料
9 U5 z4 T5 y+ n. M$ ?  Q* r*do,j,1,800
& Z# g" [1 N0 o+ b4 H& g*IF,erci(j),ge,2,THEN
- p: A: R& @) [/ Aemodif,j,mat,8
" P1 J7 K. ~" Z     *elseIF ,DFT(j),gt,0,then            
# T% |: s' Y/ a: `% v4 a% N! gemodif,j,mat,2
     
7 @4 n0 s6 L3 L( j( t1 |4 c) X; P      *elseIF,DFC(j),gt,0,then
emodif,j,mat,3
            
      *elseIF,DMT(j),gt,0,then
5 @3 S5 ], f' ]# f5 c6 R* h+ ], iemodif,j,mat,4
         
      *elseIF,DMC(j),gt,0,then
emodif,j,mat,5
            
2 F/ h5 U; G9 }# L. ?8 d3 j      *elseIF,DT(j),gt,0,then
emodif,j,mat,6
3 d: C- c. s  `% C: {            
! l+ t4 ^; [: y/ G      *elseIF,DC(j),gt,0,then
emodif,j,mat,7
$ ]6 t1 K" G% l+ X5 F4 O* n  t              
" V2 f, l' y- I*endif      
*ENDDO
* V% U  l: A4 h4 P4 _. D7 ~/solu
asel,s,loc,x,100
sfa,all,1,press,-i*200
alls
) r0 K4 I, }- U" ~!nsubst,5
!cnvtol,f,,0.1,,0.1
0 ?) `% l% {4 V. K; u2 g!ncnv,2
solve

3 W+ u, p; ^6 j2 w8 C& O- w
*enddo
1 I# P8 c! T* y