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

kolarxiao

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

yinjvictor

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

limulin1

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

我想也是，看来你也是咯

dragon_w

我模仿您的apdl写了段，可还是结果无法提取，请教！
6 {6 k  G1 n, `" L/ f) cfini
/clear
/prep7
et,1,46
' l$ e. m% y- Y+ V2 W  Wkeyopt,1,8,1
1 D! g6 b, M8 C$ u8 n# I8 Q" ER,1,4
RMODIF,1,13,1,0, 0.125,1,45,0.125
% g' m# ^0 A0 i+ f8 Z$ |4 i' `, IRMODIF,1,19,1,-45,0.125,1,90,0.125
7 ]. s1 z8 n0 T. N+ e! l. l1 pMP,EX,1,114000
% Z; Z4 e/ x# S5 J3 b* K  m/ wMP,EY,1,8610
4 q5 S* D! K2 G4 ]5 F) g. q: t2 DMP,EZ,1,8610
9 m" o* e/ r& HMP,PRXY,1,0.3
MP,PRYZ,1,0.45
MP,PRXZ,1,0.3
$ K; X7 j" `2 Z) s$ n$ c) U) x! P6 qMP,GXY,1,4160
MP,GYZ,1,3000
  u8 z5 O' u5 V/ BMP,GXZ,1,4160

! I& j; v% X% t, e& qMP,EX,2,114000
MP,EY,2,8610*0.3
! H9 n9 u! ^3 ]1 [  fMP,EZ,2,8610
$ z, J5 v0 A& Z0 ?" u3 |$ ~MP,PRXY,2,0.3
) x, Q* K# [% ]; a" y6 E/ _MP,PRXZ,2,0.3
MP,PRYZ,2,0.45
  C* j% q7 \1 \' @5 ~  ^! vMP,GXY,2,4160*0.3
MP,GXZ,2,4160
( a2 m- ~# K# g5 RMP,GYZ,2,3000*0.3

8 H% r2 m, I( L8 IMP,EX,3,114000
MP,EY,3,8610*0.4
MP,EZ,3,8610
MP,PRXY,3,0.3
MP,PRXZ,3,0.3
0 E7 E$ P+ V3 X/ B, WMP,PRYZ,3,0.45
! O" P; c0 N5 B; W: Z. _* C  ?0 P) b! kMP,GXY,3,4160*0.4
6 m$ A7 p- _7 `$ }: H4 eMP,GXZ,3,4160
6 N+ A8 c1 X: n, }MP,GYZ,3,3000*0.4
" [* }( L) N2 e
2 ?- _. ^9 J& j2 ~MP,EX,4,114000*0.07
MP,EY,4,8610
6 p' P5 B# x3 M+ u5 I1 v3 g: E' V" V- MMP,EZ,4,8610
4 S% W8 P6 i% ]1 HMP,PRXY,4,0.3
MP,PRXZ,4,0.3
MP,PRYZ,4,0.45
0 Y" ]  g  T2 x; J) A# o7 [: N% @; fMP,GXY,4,4160*0.07
# o/ v" {6 ?( [) pMP,GXZ,4,4160*0.07
MP,GYZ,4,3000

MP,EX,5,114000*0.14
' h7 q9 l  N1 R- f( cMP,EY,5,8610
MP,EZ,5,8610
MP,PRXY,5,0.3
( L2 }& {* R3 ^* Q: W% n1 qMP,PRXZ,5,0.3
- f2 ]! K1 L8 j& M, t! e) V6 S$ QMP,PRYZ,5,0.45
, n' u( s( u% a- @9 fMP,GXY,5,4160*0.07
: `/ ^) J3 g- }% n: KMP,GXZ,5,4160*0.07
6 H# Z+ M: ]5 H' p( LMP,GYZ,5,3000

MP,EX,6,114000
MP,EY,6,8610*0.3
MP,EZ,6,8610*0.3
MP,PRXY,6,0.3
; u  G: K" @8 L1 p& g0 A' ^MP,PRXZ,6,0.3
3 u5 T0 F7 r; s4 _; a  E2 sMP,PRYZ,6,0.45
MP,GXY,6,4160
* \  F. d+ L; mMP,GXZ,6,4160*0.3
MP,GYZ,6,3000*0.3

MP,EX,7,114000
# M' U1 j9 }* p0 j2 m& oMP,EY,7,8610*0.4
MP,EZ,7,8610*0.4
9 l, U# |6 S& C3 c. O, x4 J& k8 rMP,PRXY,7,0.3
# T+ N0 V1 z  l" mMP,PRXZ,7,0.3
! T9 G( L. a; d% E  T, W( H& sMP,PRYZ,7,0.45
MP,GXY,7,4160
MP,GXZ,7,4160*0.4
MP,GYZ,7,3000*0.4

) k9 i1 r  m7 |2 RMP,EX,8,114000*0.07
' ?. s% N  q0 X3 l$ B9 K4 u3 }MP,EY,8,8610*0.3
% v* z# m  z% d1 R8 Q4 {  i/ PMP,EZ,8,8610*0.3
3 i: p  }& U3 D7 ~: M3 `* CMP,PRXY,8,0.3
MP,PRXZ,8,0.3
4 ?6 W- R0 e& mMP,PRYZ,8,0.45
- S* x5 r& k6 Q: Q# ]% `- t# EMP,GXY,8,4160*0.07
& Q  a5 ?8 u* i! D( k; W7 uMP,GXZ,8,4160*0.07
! V/ g1 [9 T% i3 i- eMP,GYZ,8,3000*0.3
9 c/ d' B& D: X
!mp,ex,2,210000
!mp,prxy,2,0.3
" s5 d; n- r( }$ `& }: w* V5 L) ?Xc=914.3
Xt=1229.4
1 {% l$ D! p6 {% }/ \Yc=147.2
Yt=45.6
Zc=147.2
Zt=45.6
! H) p/ L0 `; o, V+ N6 c) [S12=33.0
: R0 D+ i' I6 d( g9 e! N' eS13=33.0
4 f! S$ j; U/ k7 s; u; k' PS23=25.0
* G( p# ]/ W+ Z" T+ V2 G*dim,ST11,array,800
) l; v2 a5 Q0 C! J0 r# |2 t*dim,ST22,array,800
* i- K, Z7 g. P*dim,ST33,array,800
*dim,ST12,array,800
0 A7 L! n6 E. X*dim,ST13,array,800
*dim,ST23,array,800
! @) F7 t& |' [- c2 o9 a' q$ B*dim,DFT,array,800  
*dim,DFC,array,800
7 |8 ~* e) {* _*dim,DMT,array,800
*dim,DMC,array,800
& K1 }, m" I! T*dim,DFM,array,800
5 m" ?- \8 S/ T  O*dim,DT,array,800
*dim,DC,array,800
*dim,erci,array,800
*dim,HASHINDFT,array,800
+ J1 B, v) [. D9 f*dim,HASHINDFC,array,800      
9 ]8 W, X' q# T*dim,HASHINDMT,array,800
: u0 _% `' W! K% R# j% C4 f*dim,HASHINDMC,array,800
*dim,HASHINDFM,array,800
+ E# |2 Y' `" u* @& A*dim,HASHINDT,array,800
*dim,HASHINDC,array,800
3 K' m. V  }5 I8 a
, G% g2 b' A$ X  R! l3 v$ H- @block,0,100,0,50,0,0.5
LSEL,S,LENGTH,,100 !选择所有长度为90mm的线段(长度方向上的线段)
& |' @* x: H- ^  B) iLESIZE,all, , ,20, , , , ,1
/ _2 O& j6 U$ e0 q3 ^! ~lsel,s,length,,50 !选择所有长度为25mm的线段(宽度方向上的线段)
lesize,all, , ,10,, , , ,1
. v* S+ `: s8 f1 SLSEL,S,LENGTH,,100
$ n" \- P6 O: d3 ]7 o8 Ilsel,a,length,,50
lsel,inve
: F) F2 E" n9 T* G  c( p+ tlesize,all, , ,4
alls
. v5 `6 i1 i3 C* ~8 C* vvatt,1,1,1
0 s' H$ j) S& c$ K9 Qvmesh,all
2 C3 B$ i# C  @  ]' G, a9 e8 F

1 v& }" r% @& @: d0 v( R/solu
; P3 V" z' ~/ W) O2 A, s; E6 }asel,s,loc,x,0
da,all,all
# F% I+ v3 |# [2 N1 R$ Y+ Zalls
6 V1 D! V. G( e! r/ qasel,s,loc,x,100
% z0 A# u  T7 [9 u. \! q1 isfa,all,1,press,-2
alls
!nsubst,5
!cnvtol,f,,0.1,,0.1
; g( e9 G) V- E# m!rescontrol,define,all,1,1
! L) S, i# K# D- K( j! Wsolve
/ \, l1 Y# x7 P1 t. k( v# l!save,1,db
: w, Z" j9 M; n& D* L9 j! H* c2 s. q. ^
$ I; e( m6 N  {+ `*do,i,2,10
8 X2 y- k5 ?0 E7 v" N# R  /post1
! I4 ~! P; i/ ~2 [% m( q% d   etable,sx,s,x
4 L2 }0 w, G% Q& s* ]# y   etable,sy,s,y
   etable,sz,s,z
   etable,sxy,S,xy
) h, Y. c0 \0 q3 ?   etable,syz,s,yz
4 @. M/ n; N7 u# A/ `% r* |5 A! L! R" S   etable,sxz,S,xz
   
   *do,j,1,800
      *get,ST11(j),elem,j,etab,sx
6 e( J& Q0 q; B+ {      *get,ST22(j),elem,j,etab,sy
$ L# w' d) {: k, w0 ~$ [* a      *get,ST33(j),elem,j,etab,sz
      *get,ST12(j),elem,j,etab,sxy
) q# {# Y; k3 K! g/ r1 S      *get,ST13(j),elem,j,etab,sxz
      *get,ST23(j),elem,j,etab,syz                  
* I4 f. l6 C2 ^  I   *enddo
$ x1 m3 h( v6 l   *do,j,1,800
( p: M: j) D/ }" R, z: I/ i      HASHINDFT(j) = (ST11(j)/Xt)**2+(ST12(j)/S12)**2+(ST13(j)/S13)**2
      HASHINDFC(j) = (ST11(j)/Xc)**2
9 a$ X( u+ \4 c0 ^6 l      HASHINDMT(j) = (ST22(j)/Yt)**2+(ST12(j)/S12)**2+(ST23(j)/S23)**2
7 q8 {, ~) k* z, [! z# \3 A: o      HASHINDMC(j) = (ST22(j)/Yc)**2+(ST12(j)/S12)**2+(ST23(j)/S23)**2
2 z( q1 ^7 f& g3 o# T9 ?' |" `% s4 c      HASHINDFM(j) = (ST11(j)/Xc)**2+(ST12(j)/S12)**2+(ST13(j)/S13)**2
$ h/ g4 u* w7 Q0 h/ e: N& B      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              
   *enddo   
parsav,all,parameter,txt
3 }, @# A, I2 [+ M6 E
  /solu
    antype,,rest
4 i0 T& J) y4 k; n. o: \parres,new,parameter,txt
!失效类型
9 p0 T- C( {8 Q( X) q, |% R$ \8 ~8 A!纤维失效
*do,j,1,800
$ R( b: v- H, l# C1 N# ?* d      *IF,ST11(j),ge,0,and,HASHINDFT(j),ge,1,THEN
$ W. {/ a/ t7 e$ H4 Q" M: r               DFT(j)=1+DFT(j)
8 v4 C) D4 U; d$ O) k
            *ELSEif,ST11(j),lt,0,and,HASHINDFC(j),ge,1,then
% Y2 H# g9 W% `2 t2 z2 \. q               DFC(j)=1+DFC(j)
5 i9 k: W! \- l0 x$ `  C/ m
. l6 d! I/ Q7 Y" c3 H* n            *ELSE
0 p; I6 Y' e. ^/ E; |0 A9 a                DFC(j)=0+DFC(j)
2 U! d7 |: R9 ?9 M  O* Q                DFT(j)=0+DFT(j)
           
      *ENDIF
) a7 [0 u% q8 I. C7 g0 O: v  l/ V*enddo

# u3 _8 h7 ]) y!基体失效
+ x' ^/ v/ B$ [; _  n& ?' K*do,j,1,800
      *IF,ST22(j),ge,0,and,HASHINDMT(j),ge,1,THEN         
- |7 `8 I) i7 K* x6 f6 O               DMT(j)=1+DMT(j)
         
6 B* C* I! T# L          *ELSEIF,ST22(j),LT,0,and,HASHINDMC(j),ge,1,THEN
                DMC(j)=1+DMC(j)

/ }: m4 |9 |; {' m( f          *ELSE
( u) ?5 P( u, W                DMC(j)=0+DMC(j)
                DMT(j)=0+DMT(j)
         
" w, I% \2 U9 {      *ENDIF
*enddo
!FIBER-MATRIX FAILURE
" P. ?0 a$ Y6 u* M# M2 q# {*do,j,1,800
9 y0 p: m' P& \      *IF,ST11(j),ge,0,AND,HASHINDFM(j),ge,1,THEN
              DFM(j)=1+DFM(j)

            *ELSE
! q; W  k5 J; H8 G! I              DFM(j)=0+DFM(j)
0 H( Z5 y, C% x1 c7 ?; y( p( S            
# p7 T  J, s9 \6 c# c3 M      *ENDIF
; {3 N# G, k, |( ]/ e( N. u6 r*enddo
! L9 ~8 b/ y: y5 G& j6 ?. F" {
!Tensile-Compression FAILURE
1 g- z/ P& k8 ]/ K% W% A- f*do,j,1,800
      *IF,ST33(j),ge,0,AND,HASHINDT(j),ge,1,THEN
9 y: I7 @  \: `' o! I7 J+ I               DT(j)=1+DT(j)

( Q0 x6 q  y" d6 ^, q* e          *ELSEIF,ST33(j),LT,0,AND,HASHINDC(j),ge,1,THEN
                DC(j)=1+DC(j)

+ F6 t( @6 ]2 A7 l9 b          *ELSE
                DC(j)=0+DC(j)
- a0 i. W' }- B" i: e; H4 h& j                DT(j)=0+DT(j)
         
" {" M! ?0 N  ^' ?% D      *ENDIF
7 N6 E) l7 ~3 s% Q- b7 E* w( d*enddo
8 t5 D) O5 @% _; U*do,j,1,800
    erci(j)=DFT(j)+DFC(j)+DMT(j)+DMC(j)+DFM(j)+DT(j)+DC(j)
*ENDDO
/prep7
!更新单元材料
6 c2 Y2 p/ o! x2 `# i( v  x% p! V*do,j,1,800
*IF,erci(j),ge,2,THEN
3 e1 Y0 A* f# b" hemodif,j,mat,8
, F/ ^1 n# |5 {! S8 q4 N+ ]! h     *elseIF ,DFT(j),gt,0,then            
; r# m$ @* n4 X; r8 `- f5 D$ O5 Temodif,j,mat,2
2 u$ @# L. H+ d# H9 ]( v% u8 T     
      *elseIF,DFC(j),gt,0,then
emodif,j,mat,3
            
& @+ A9 Q& B) V2 @+ t" w      *elseIF,DMT(j),gt,0,then
emodif,j,mat,4
2 e* Z) y9 p$ q: h0 o         
      *elseIF,DMC(j),gt,0,then
1 |6 G9 l  X) m$ E0 u8 V. x% ^emodif,j,mat,5
8 D; r, K3 I, u            
9 L; z$ p5 o2 o! ~7 d      *elseIF,DT(j),gt,0,then
  P; P6 }3 q* z& |  w6 ?emodif,j,mat,6
/ G9 U" U( F+ {2 c6 u; s            
( J/ h, M: b0 C; u& ?6 ?! Y      *elseIF,DC(j),gt,0,then
emodif,j,mat,7
              
*endif      
*ENDDO
1 n) w8 d! N* |4 a* ~8 ~' N3 g/solu
asel,s,loc,x,100
# y1 h1 w4 H' N! Q/ f, bsfa,all,1,press,-i*200
. @% U4 R3 g0 ]- Kalls
!nsubst,5
' n7 H) H1 C5 F1 z/ x. e!cnvtol,f,,0.1,,0.1
!ncnv,2
solve

2 z7 [6 ?4 f8 }; E0 E4 }
0 [+ l4 i4 o) C! A+ @" l/ f. |*enddo
$ D9 M& \0 H% b' Z