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

kolarxiao

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

yinjvictor

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

limulin1

sxyuyifan 发表于 2010-2-4 08:26
楼主真强，是浙大化机所的啊。

" z0 H  s9 P7 L0 h8 J) N5 Y  p我想也是，看来你也是咯

dragon_w

我模仿您的apdl写了段，可还是结果无法提取，请教！
9 x  V- n: I- Z/ D' U7 J+ ]fini
/clear
' [& W* A+ }9 z1 `; S/prep7
et,1,46
keyopt,1,8,1
& w# l7 q+ T! Z7 L8 ?/ w3 x  oR,1,4
) ?! M2 |5 [5 d: T& RRMODIF,1,13,1,0, 0.125,1,45,0.125
RMODIF,1,19,1,-45,0.125,1,90,0.125
. [: z/ {0 ?. j6 t' t. _MP,EX,1,114000
9 h0 Q6 \+ n+ u  R! B$ ~' KMP,EY,1,8610
MP,EZ,1,8610
7 i( o$ l0 v4 N( i) pMP,PRXY,1,0.3
& i  `) s5 @8 R* {" `5 BMP,PRYZ,1,0.45
MP,PRXZ,1,0.3
" G$ r5 e  o# T' S* j% cMP,GXY,1,4160
  n7 {1 r$ G1 x& L+ T: o9 N, \3 CMP,GYZ,1,3000
MP,GXZ,1,4160

MP,EX,2,114000
MP,EY,2,8610*0.3
MP,EZ,2,8610
3 F) d4 }1 t7 ~2 \. gMP,PRXY,2,0.3
MP,PRXZ,2,0.3
! C% C+ z8 j+ ~8 pMP,PRYZ,2,0.45
MP,GXY,2,4160*0.3
MP,GXZ,2,4160
/ O: [) ], h3 Y0 \* L3 `9 u6 LMP,GYZ,2,3000*0.3
$ R( ]1 U+ I) B* T# X! c
: C3 {1 h  N+ o& i- mMP,EX,3,114000
MP,EY,3,8610*0.4
MP,EZ,3,8610
MP,PRXY,3,0.3
0 e4 O1 y. o5 OMP,PRXZ,3,0.3
( }$ s+ p* K5 TMP,PRYZ,3,0.45
MP,GXY,3,4160*0.4
MP,GXZ,3,4160
, [% S- R6 `' o/ CMP,GYZ,3,3000*0.4

MP,EX,4,114000*0.07
' e% N$ `% J4 Q2 B- Z& T/ cMP,EY,4,8610
MP,EZ,4,8610
+ g" L& \; v5 YMP,PRXY,4,0.3
MP,PRXZ,4,0.3
MP,PRYZ,4,0.45
8 @$ a9 a% g; p. z9 _9 xMP,GXY,4,4160*0.07
3 b1 A  a* e1 ~0 }: T. v6 GMP,GXZ,4,4160*0.07
MP,GYZ,4,3000

; F7 I$ b* j: [3 B. s$ Y) L+ dMP,EX,5,114000*0.14
MP,EY,5,8610
" w! f7 X! h( P. c- V% _MP,EZ,5,8610
2 S6 |$ ?4 |5 \4 AMP,PRXY,5,0.3
MP,PRXZ,5,0.3
MP,PRYZ,5,0.45
MP,GXY,5,4160*0.07
. _$ R! K6 P/ T1 L7 J- ?% EMP,GXZ,5,4160*0.07
9 ]1 G  n% {/ G* {$ v( LMP,GYZ,5,3000
; Y# T4 v0 ]$ ^$ D) c
MP,EX,6,114000
% |8 ]5 m! f; M3 i! E: mMP,EY,6,8610*0.3
MP,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
6 _3 X; w  ^% G# J' l3 [8 MMP,GYZ,6,3000*0.3

5 K+ ?( {) N: z! N* JMP,EX,7,114000
4 o* H. t  J' l$ aMP,EY,7,8610*0.4
4 v6 V' ^( R! r4 ]7 z1 Y7 o: o; w) bMP,EZ,7,8610*0.4
5 R0 m0 u7 W8 l. j! D' }$ g- O1 FMP,PRXY,7,0.3
MP,PRXZ,7,0.3
MP,PRYZ,7,0.45
MP,GXY,7,4160
MP,GXZ,7,4160*0.4
( V! _/ p9 K6 ]7 _9 x6 fMP,GYZ,7,3000*0.4

+ ]- {. u! i$ k, FMP,EX,8,114000*0.07
MP,EY,8,8610*0.3
MP,EZ,8,8610*0.3
% R3 }! Z% Y9 _3 q$ Q1 qMP,PRXY,8,0.3
3 C& q% A. R; i4 }  E4 V6 iMP,PRXZ,8,0.3
: I& Y- ]" `  J' n  t- _7 eMP,PRYZ,8,0.45
MP,GXY,8,4160*0.07
6 O1 B: d7 [- D4 S/ ?0 kMP,GXZ,8,4160*0.07
  r+ \* L3 w/ h- S+ M7 p9 GMP,GYZ,8,3000*0.3

4 T; ~8 W1 q+ }3 ^+ e0 x8 Q!mp,ex,2,210000
# D7 w$ t& Y0 V' C) ^8 ~3 L!mp,prxy,2,0.3
Xc=914.3
Xt=1229.4
Yc=147.2
Yt=45.6
Zc=147.2
$ b' i, N! O7 {Zt=45.6
S12=33.0
S13=33.0
S23=25.0
8 N3 a+ G! v, J$ Q, T# R*dim,ST11,array,800
*dim,ST22,array,800
/ r5 J6 v7 {5 P*dim,ST33,array,800
*dim,ST12,array,800
: X4 \7 p7 U8 k8 L2 _. S*dim,ST13,array,800
*dim,ST23,array,800
*dim,DFT,array,800  
*dim,DFC,array,800
*dim,DMT,array,800
/ Q6 V' |3 J1 Q; ]$ \% w& v*dim,DMC,array,800
/ c; \( ~+ O  P" B1 i4 i0 _*dim,DFM,array,800
*dim,DT,array,800
: {; E: L* _5 M8 Y( O  V, F*dim,DC,array,800
*dim,erci,array,800
*dim,HASHINDFT,array,800
( M2 N* @: \& n' R$ O*dim,HASHINDFC,array,800      
' Z. b9 f$ P% C- Q  v' g*dim,HASHINDMT,array,800
; j5 O: P% ?/ C( ~*dim,HASHINDMC,array,800
; W: t, H: G6 {8 z3 ~9 y( P. L*dim,HASHINDFM,array,800
- f( n  O& f. }# o. Q*dim,HASHINDT,array,800
*dim,HASHINDC,array,800
  ^* ]" N2 C/ `* X. A
+ w+ I& h5 B5 L3 q; H. lblock,0,100,0,50,0,0.5
$ {7 ]1 v1 u. I1 Q' {/ @LSEL,S,LENGTH,,100 !选择所有长度为90mm的线段(长度方向上的线段)
2 f' z8 C7 I8 _3 f' OLESIZE,all, , ,20, , , , ,1
& m4 U, ^- G4 y& h4 Flsel,s,length,,50 !选择所有长度为25mm的线段(宽度方向上的线段)
: \1 D: u0 A) b* `# g. olesize,all, , ,10,, , , ,1
LSEL,S,LENGTH,,100
lsel,a,length,,50
+ N1 V, M' A& G2 z8 x/ Elsel,inve
. T; g, ~1 s) }/ ]lesize,all, , ,4
  l/ E' B3 R% @/ |, u/ T" q* k7 ualls
vatt,1,1,1
vmesh,all
' i9 f, X3 [/ e  [$ l: H7 @

" [4 _9 \, F( R& |' f* Q/solu
9 k( g. ]) H% e6 z# dasel,s,loc,x,0
da,all,all
1 F/ c" j1 d, i9 z+ galls
+ V1 j/ V6 O' F# T( \" A1 H0 Iasel,s,loc,x,100
) N. R, J$ {8 h' _8 n& F  `sfa,all,1,press,-2
alls
* G7 E8 c' ?, M6 X* e. |& k!nsubst,5
!cnvtol,f,,0.1,,0.1
2 O! P: D" {/ u; q& _( F!rescontrol,define,all,1,1
solve
5 J8 g/ G" z7 Q: a!save,1,db

5 Z( S7 y/ X8 `' d/ _: d*do,i,2,10
' f& ~7 Q  X* S% b1 m3 j  /post1
% {0 Y" n# g. H3 G# Q! A9 G   etable,sx,s,x
8 }& r' @3 s' z4 j7 G   etable,sy,s,y
: i1 q  q6 x9 |) b& r$ r   etable,sz,s,z
   etable,sxy,S,xy
   etable,syz,s,yz
   etable,sxz,S,xz
1 i4 v# J+ K" S1 `/ d   
   *do,j,1,800
0 R( X% K) d3 L) e5 x      *get,ST11(j),elem,j,etab,sx
      *get,ST22(j),elem,j,etab,sy
- A* j$ S. j5 j: x: }: ~& l( c      *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 E( ?0 X0 ^" }, Z- U   *enddo
% G: P2 Q2 y/ P7 X6 E6 `! J   *do,j,1,800
+ k' n- r4 Z! j# q3 G% Z# p% i      HASHINDFT(j) = (ST11(j)/Xt)**2+(ST12(j)/S12)**2+(ST13(j)/S13)**2
( m5 i: q  Z/ A% U& W; _+ S/ T4 G      HASHINDFC(j) = (ST11(j)/Xc)**2
      HASHINDMT(j) = (ST22(j)/Yt)**2+(ST12(j)/S12)**2+(ST23(j)/S23)**2
" a# n! {# l9 `& A/ P. y4 F& |9 ?/ r5 l$ U      HASHINDMC(j) = (ST22(j)/Yc)**2+(ST12(j)/S12)**2+(ST23(j)/S23)**2
# H- I* T% ~2 E4 A      HASHINDFM(j) = (ST11(j)/Xc)**2+(ST12(j)/S12)**2+(ST13(j)/S13)**2
0 m' H0 \2 U6 X/ M0 z      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              
( q- x4 o4 m: ?+ T' t' i+ O   *enddo   
parsav,all,parameter,txt

9 T  O$ K( O  O- z4 z( N  /solu
    antype,,rest
8 ~8 C9 b8 M5 D% g. q& Cparres,new,parameter,txt
/ E  e- z; o9 c+ Z  T9 v6 V!失效类型
!纤维失效
& T" f7 c: p& g! c*do,j,1,800
      *IF,ST11(j),ge,0,and,HASHINDFT(j),ge,1,THEN
$ {  i2 T2 d7 X8 y+ \' O! U: E               DFT(j)=1+DFT(j)
3 P1 |& J5 R1 _
            *ELSEif,ST11(j),lt,0,and,HASHINDFC(j),ge,1,then
( Y# \+ P/ M+ A               DFC(j)=1+DFC(j)
: C: A- N* Z1 p! n9 k9 W1 \
2 [0 A! x6 x) D, w0 L/ G            *ELSE
; D4 h1 T1 x4 a- p                DFC(j)=0+DFC(j)
6 b% ?8 D. E6 o                DFT(j)=0+DFT(j)
           
0 p, n8 Z) Y3 F2 U4 y; C      *ENDIF
7 R8 \# W  ~0 `* J0 M4 a( z*enddo

!基体失效
) X1 z2 w" s# B# Q% L' p*do,j,1,800
3 L: g/ S! m- O# T      *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)

+ W6 X6 O% y) H  h          *ELSE
                DMC(j)=0+DMC(j)
4 y0 J" G2 M: M' V/ {                DMT(j)=0+DMT(j)
         
# a4 h/ Z( L9 \$ }" A9 e      *ENDIF
*enddo
!FIBER-MATRIX FAILURE
*do,j,1,800
6 T& o' H2 Z6 {7 l# s, R# h      *IF,ST11(j),ge,0,AND,HASHINDFM(j),ge,1,THEN
8 \: @+ v, d6 ]* b: S: a5 A3 O              DFM(j)=1+DFM(j)
, a5 k& n2 M  u% q/ e0 g2 Z
            *ELSE
. e. [' {& L. B/ Q6 k              DFM(j)=0+DFM(j)
            
/ C# R$ z& j1 f. E' Z" y      *ENDIF
*enddo

!Tensile-Compression FAILURE
* T: Y4 a  Q; D, P) E! e- f  B5 p*do,j,1,800
      *IF,ST33(j),ge,0,AND,HASHINDT(j),ge,1,THEN
               DT(j)=1+DT(j)

          *ELSEIF,ST33(j),LT,0,AND,HASHINDC(j),ge,1,THEN
                DC(j)=1+DC(j)

          *ELSE
! Y# K7 m( B% K7 s" C% m                DC(j)=0+DC(j)
                DT(j)=0+DT(j)
: Y5 z3 r' k0 z5 o) ^& f* Q6 u) ]" Y         
      *ENDIF
' M8 B1 C/ m& P% p6 l, l1 q) G' M*enddo
; s& \+ w, m. W) I9 u- e4 d*do,j,1,800
* M% k, |: ]3 Q! X% z$ D, K$ S    erci(j)=DFT(j)+DFC(j)+DMT(j)+DMC(j)+DFM(j)+DT(j)+DC(j)
1 e# w& A4 R4 U9 U7 W* U  b*ENDDO
! q. I" v, o! E5 v) ~/ w/prep7
!更新单元材料
*do,j,1,800
*IF,erci(j),ge,2,THEN
) G  K* V8 b  n+ Memodif,j,mat,8
: R; P# f& }: V! K) L     *elseIF ,DFT(j),gt,0,then            
- z/ j- ^  e- cemodif,j,mat,2
# [/ h0 Q# ]/ K- ?+ c     
* q* ~% A5 F; G  J6 X2 s+ S      *elseIF,DFC(j),gt,0,then
$ O3 f2 k( h% I: y9 o( lemodif,j,mat,3
            
      *elseIF,DMT(j),gt,0,then
8 J9 i. t8 o+ b) p( Eemodif,j,mat,4
* P1 N8 H( o  \" B& A% N/ l         
/ @! s8 x" H8 A( L, Y( q$ [      *elseIF,DMC(j),gt,0,then
1 I% u# u6 U+ b  {- Memodif,j,mat,5
& Z7 P- ?; }! q            
% \+ O! V; Z% C+ \& G9 c      *elseIF,DT(j),gt,0,then
& H% t! R% Q$ r( d$ F% Cemodif,j,mat,6
            
      *elseIF,DC(j),gt,0,then
emodif,j,mat,7
              
*endif      
, p3 X9 n" |" `1 p& F( \2 ~*ENDDO
4 r% g# i2 X) C& ]/solu
2 X7 s7 a" ~8 U/ aasel,s,loc,x,100
6 r* h; ]2 @5 Psfa,all,1,press,-i*200
alls
$ f4 r+ _) Q1 ]!nsubst,5
!cnvtol,f,,0.1,,0.1
!ncnv,2
4 Z, }5 L4 H) _; R- x: bsolve
, y. e3 j" N$ J9 q4 X$ \: z' `6 y

*enddo
% t- b: x. \, S4 {3 V8 y