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

kolarxiao

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

yinjvictor

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

limulin1

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

我想也是，看来你也是咯

dragon_w

我模仿您的apdl写了段，可还是结果无法提取，请教！
fini
/clear
7 M+ w0 q; ]( E) a& w( U2 s1 k/prep7
et,1,46
( o! r" q* O! H- d2 o  `3 r0 K1 N& _5 Mkeyopt,1,8,1
R,1,4
7 p* i: F# t$ M# ~8 k  uRMODIF,1,13,1,0, 0.125,1,45,0.125
# R" r/ |) }3 h$ Y! f% Y+ p; x- RRMODIF,1,19,1,-45,0.125,1,90,0.125
$ F8 r, ^$ i9 M3 W0 xMP,EX,1,114000
MP,EY,1,8610
4 `3 m( p1 C6 b7 d2 ?. \. VMP,EZ,1,8610
8 O3 V9 C2 q& o5 P* ?MP,PRXY,1,0.3
MP,PRYZ,1,0.45
MP,PRXZ,1,0.3
% g3 q& S. x: d+ oMP,GXY,1,4160
MP,GYZ,1,3000
MP,GXZ,1,4160

% y1 h2 L% A8 h4 N5 ^; ]  n! ]$ ]MP,EX,2,114000
4 \1 p8 Y1 F; R8 B3 v6 kMP,EY,2,8610*0.3
MP,EZ,2,8610
* |/ g. r. X% h) JMP,PRXY,2,0.3
MP,PRXZ,2,0.3
MP,PRYZ,2,0.45
, C" F6 H1 C" Q% r9 UMP,GXY,2,4160*0.3
* C5 v* ^/ f0 sMP,GXZ,2,4160
MP,GYZ,2,3000*0.3
, W0 I* A% g: O+ M9 a
+ u( ~  `' m% x8 X. ~MP,EX,3,114000
' v! \4 o4 I# }MP,EY,3,8610*0.4
MP,EZ,3,8610
7 q4 B) f/ ?! w" z- ^MP,PRXY,3,0.3
MP,PRXZ,3,0.3
MP,PRYZ,3,0.45
MP,GXY,3,4160*0.4
4 e- t  e, S% z; NMP,GXZ,3,4160
( @. G8 f0 u7 w; vMP,GYZ,3,3000*0.4
6 T- r  [9 o; Q7 v
: a, g! X& u3 h3 o3 @- a9 x* oMP,EX,4,114000*0.07
MP,EY,4,8610
6 b8 i: ?, m5 C! m* cMP,EZ,4,8610
MP,PRXY,4,0.3
$ I/ l) J' X# \, _9 l0 EMP,PRXZ,4,0.3
MP,PRYZ,4,0.45
2 `: V$ E2 M8 z2 a* `% @2 YMP,GXY,4,4160*0.07
- ?1 n# I7 K9 N& x4 l2 KMP,GXZ,4,4160*0.07
MP,GYZ,4,3000

MP,EX,5,114000*0.14
MP,EY,5,8610
MP,EZ,5,8610
* M+ Q, R) G. Y. {; s; ?$ Y1 y3 eMP,PRXY,5,0.3
( b+ p& Z6 S- u/ dMP,PRXZ,5,0.3
/ Q+ y, K/ \7 d5 hMP,PRYZ,5,0.45
MP,GXY,5,4160*0.07
3 g- e9 ^9 L* i' fMP,GXZ,5,4160*0.07
, w& L* v. _* J: v! G+ y/ HMP,GYZ,5,3000

. N: y7 R/ u4 X, Y0 Q5 N7 BMP,EX,6,114000
9 E+ T: F3 k/ J9 _4 `MP,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
- ?2 g" n, c1 D! d2 F* L8 tMP,GXZ,6,4160*0.3
MP,GYZ,6,3000*0.3
8 V- _5 @" \$ h6 u
MP,EX,7,114000
MP,EY,7,8610*0.4
! z1 ]$ ]; X6 c" gMP,EZ,7,8610*0.4
' [+ I+ W2 U* m( D* ~MP,PRXY,7,0.3
; @- m  H& N% k2 qMP,PRXZ,7,0.3
MP,PRYZ,7,0.45
% D1 y; s: y  E# V4 T& EMP,GXY,7,4160
MP,GXZ,7,4160*0.4
# m3 |, y* r1 ~* K% aMP,GYZ,7,3000*0.4

MP,EX,8,114000*0.07
MP,EY,8,8610*0.3
& j7 I7 O4 m: }: G, dMP,EZ,8,8610*0.3
MP,PRXY,8,0.3
# h8 M- d$ A% I( Q# B; c5 bMP,PRXZ,8,0.3
MP,PRYZ,8,0.45
. Y9 P  d3 f; ^( jMP,GXY,8,4160*0.07
MP,GXZ,8,4160*0.07
- ]6 z2 E) \5 q% Y8 d2 c8 {MP,GYZ,8,3000*0.3

; c0 [/ D; T& R( ^# ^2 S8 T!mp,ex,2,210000
) m, ~) ~' n/ u4 n% M/ Y) R!mp,prxy,2,0.3
- ^0 a4 `! }6 D  M$ H' S& o/ gXc=914.3
Xt=1229.4
/ p7 d! E" H$ {Yc=147.2
4 x) J& ]. M( \& ZYt=45.6
Zc=147.2
, F; M  ~. G7 u8 L2 D, {Zt=45.6
S12=33.0
S13=33.0
, t, v; ?! j0 ^8 B' j  IS23=25.0
*dim,ST11,array,800
*dim,ST22,array,800
! w, q4 w. k1 |3 g4 L0 B# g*dim,ST33,array,800
*dim,ST12,array,800
4 b% ]0 d4 T, g2 j/ T/ G# e5 [*dim,ST13,array,800
*dim,ST23,array,800
*dim,DFT,array,800  
8 h$ o1 D/ @4 J  V  u3 @' a, a*dim,DFC,array,800
- l( }5 I0 r" A( w  F: G) p/ u*dim,DMT,array,800
8 }  x5 E7 u, O: ]* x7 ^% ^*dim,DMC,array,800
*dim,DFM,array,800
*dim,DT,array,800
*dim,DC,array,800
1 `, `5 C) W1 N! {! }*dim,erci,array,800
*dim,HASHINDFT,array,800
*dim,HASHINDFC,array,800      
( h( z) i5 V, [* \/ h% s8 \8 `*dim,HASHINDMT,array,800
*dim,HASHINDMC,array,800
*dim,HASHINDFM,array,800
*dim,HASHINDT,array,800
*dim,HASHINDC,array,800
, t4 Z3 _  k( f3 A0 z# X
block,0,100,0,50,0,0.5
LSEL,S,LENGTH,,100 !选择所有长度为90mm的线段(长度方向上的线段)
LESIZE,all, , ,20, , , , ,1
lsel,s,length,,50 !选择所有长度为25mm的线段(宽度方向上的线段)
lesize,all, , ,10,, , , ,1
LSEL,S,LENGTH,,100
lsel,a,length,,50
lsel,inve
, Y7 ~) L+ b9 [7 m1 a& J$ A/ ^lesize,all, , ,4
alls
' u1 `# \9 ]. _- y4 xvatt,1,1,1
; E# Q; j8 T& k: C$ R- ]vmesh,all
6 d: g% n2 A2 k3 x
( l! I$ p# _. v0 [
: e, M( i% _, s/solu
# m" q$ [; h( m* i2 m7 masel,s,loc,x,0
/ f7 s: W  g  m# f4 e3 Yda,all,all
alls
$ H4 o1 e3 y/ C+ g+ z4 l& j$ A3 }4 dasel,s,loc,x,100
1 ]" C; l) ]2 F4 h$ K6 |sfa,all,1,press,-2
) U6 T+ b% x4 B% E% H: yalls
!nsubst,5
!cnvtol,f,,0.1,,0.1
!rescontrol,define,all,1,1
& e/ L. S4 a# I+ P, `7 V3 jsolve
!save,1,db

*do,i,2,10
9 {* s2 K+ S+ P1 E# e  /post1
   etable,sx,s,x
' O$ O1 U' H& \% E- c   etable,sy,s,y
   etable,sz,s,z
* i4 d' a1 ~; R9 w. ^: }   etable,sxy,S,xy
3 v2 Z/ x" u# x; D1 U$ N   etable,syz,s,yz
   etable,sxz,S,xz
   
" N  Y) M  g* d& s   *do,j,1,800
; ?# p1 v. Q. |* |  s. G      *get,ST11(j),elem,j,etab,sx
; S6 l, V! H7 Y5 m7 [      *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                  
. r; X. K1 n: {7 Z' _: H   *enddo
   *do,j,1,800
9 x5 Q1 Y& v- C9 r( d+ ]3 L      HASHINDFT(j) = (ST11(j)/Xt)**2+(ST12(j)/S12)**2+(ST13(j)/S13)**2
- g  {/ ~& c2 T+ k, x& h1 X      HASHINDFC(j) = (ST11(j)/Xc)**2
      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
. w# @8 Z; V7 U8 C1 G: A0 A* k      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   
* T# B' [. _! o$ s$ w1 x- \ parsav,all,parameter,txt
4 i+ d7 Q9 R/ [* i* Q# u
$ t, t5 q& }% o5 k- J' q, }  /solu
    antype,,rest
; q1 g8 X' s& g# c( a& X% p% Fparres,new,parameter,txt
3 O* n) P( @9 L. L; @) h5 d/ V!失效类型
% h- {( M1 {& X' {$ [!纤维失效
*do,j,1,800
9 {* f8 I' [- q$ W% N5 `9 C      *IF,ST11(j),ge,0,and,HASHINDFT(j),ge,1,THEN
               DFT(j)=1+DFT(j)
; B& @/ ]# x5 i& g8 {" ~( n
            *ELSEif,ST11(j),lt,0,and,HASHINDFC(j),ge,1,then
: a9 g$ {3 Q/ w/ V' M5 u# y# F               DFC(j)=1+DFC(j)
  O9 @8 |5 {4 q' o' ~  c' i7 A% j
0 m- I+ D* w% F( _% K' g            *ELSE
                DFC(j)=0+DFC(j)
                DFT(j)=0+DFT(j)
' \8 A& W/ c4 ^8 N# B9 h           
( E; n& F7 o8 c" r& w: y5 R0 M2 f      *ENDIF
) w; d+ Q4 Y( h# N: A*enddo
0 b5 C* C' t  V2 O' c7 N; K1 G
!基体失效
1 z7 _/ T9 t  ^8 O/ F' u) t( |*do,j,1,800
      *IF,ST22(j),ge,0,and,HASHINDMT(j),ge,1,THEN         
  r/ |+ x  s$ I7 G/ F               DMT(j)=1+DMT(j)
4 d, D( q5 g0 w6 j' V" t/ [+ y         
          *ELSEIF,ST22(j),LT,0,and,HASHINDMC(j),ge,1,THEN
                DMC(j)=1+DMC(j)

          *ELSE
                DMC(j)=0+DMC(j)
                DMT(j)=0+DMT(j)
         
      *ENDIF
*enddo
% W* ^; s9 R. b: i% y2 c; D!FIBER-MATRIX FAILURE
*do,j,1,800
3 w& C6 j) x( F5 u. t8 i      *IF,ST11(j),ge,0,AND,HASHINDFM(j),ge,1,THEN
              DFM(j)=1+DFM(j)
) t, s  p" z/ e6 J3 x6 w4 P- }# m
. J# g9 V+ V: P/ t            *ELSE
8 ]" c. R0 O2 C+ s+ Q              DFM(j)=0+DFM(j)
- v. W' ~  A" d7 W* k  r! u            
/ s- ~, J0 v3 b8 a) h7 \5 E      *ENDIF
% d2 C* E7 h' [*enddo
0 M6 H, s" T$ h- ^* ?8 R
5 D4 X9 r# a% r3 E!Tensile-Compression FAILURE
  n& X% q) L: R% P*do,j,1,800
4 z% t7 b  Z9 V- J4 f      *IF,ST33(j),ge,0,AND,HASHINDT(j),ge,1,THEN
               DT(j)=1+DT(j)
3 S4 o  Y2 L4 e3 G0 n1 Q& ^# X
          *ELSEIF,ST33(j),LT,0,AND,HASHINDC(j),ge,1,THEN
9 G8 f; t" G4 b* Z7 @                DC(j)=1+DC(j)

7 F5 ~4 V% Q" ]# L          *ELSE
                DC(j)=0+DC(j)
                DT(j)=0+DT(j)
         
      *ENDIF
& h. z& v% g  l6 {- L* r. a) L*enddo
( }4 `' s; B9 \, N3 A9 Z*do,j,1,800
. Z' r- W8 R+ e, d( [$ s    erci(j)=DFT(j)+DFC(j)+DMT(j)+DMC(j)+DFM(j)+DT(j)+DC(j)
+ ?  }7 _. J, y( Q  U* b*ENDDO
/prep7
7 {$ X" P% x+ H2 K!更新单元材料
0 U8 n5 G/ X# \# C*do,j,1,800
+ n( t7 P; |! M1 ]* e/ P2 x*IF,erci(j),ge,2,THEN
5 q/ z  o7 }' F, W8 T# Kemodif,j,mat,8
. \; S6 C  c; `9 \* Y% H5 y     *elseIF ,DFT(j),gt,0,then            
8 ~: T( O; Y4 y$ ]( A) M$ lemodif,j,mat,2
     
6 Q0 c' j, H0 U8 h" a  w      *elseIF,DFC(j),gt,0,then
/ I7 I/ |" n3 u2 A# ~emodif,j,mat,3
  g: K* M' \' e8 i+ e( _$ n            
: J3 K: {! `4 Y, a' |) c      *elseIF,DMT(j),gt,0,then
/ B% S' T7 `: `( jemodif,j,mat,4
         
: P$ S) m) b  }% ^" f( `      *elseIF,DMC(j),gt,0,then
1 j1 V  Z( P5 f9 @& M+ D+ R  x7 memodif,j,mat,5
/ g' ^% D6 U0 [) Y" D+ N/ a2 J4 j            
' {7 O( S! d$ |# @      *elseIF,DT(j),gt,0,then
emodif,j,mat,6
$ E$ H; J4 N( U+ q: ?0 o$ L            
      *elseIF,DC(j),gt,0,then
' a4 {2 P: P3 w3 femodif,j,mat,7
              
) s2 m: k: C* @! x4 J9 W*endif      
*ENDDO
/solu
- c/ I) E( v- B* e8 W; E! H* k$ y+ tasel,s,loc,x,100
sfa,all,1,press,-i*200
alls
# E" T) i& P9 \; o!nsubst,5
!cnvtol,f,,0.1,,0.1
!ncnv,2
solve

  ]1 A9 C2 b1 f9 \
3 g* v) h! G7 d5 a! L( ?8 y: Q*enddo