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

kolarxiao

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

yinjvictor

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

limulin1

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

我想也是，看来你也是咯

dragon_w

我模仿您的apdl写了段，可还是结果无法提取，请教！
fini
/clear
1 Z! n# |6 B2 f/prep7
; k, P* A) P0 H3 e: o5 H3 a, T  Vet,1,46
  z# c* _  J" s& a  Fkeyopt,1,8,1
R,1,4
, z3 {; y1 z6 n; a$ L# m2 URMODIF,1,13,1,0, 0.125,1,45,0.125
RMODIF,1,19,1,-45,0.125,1,90,0.125
MP,EX,1,114000
7 E( T+ d  L  v# f( IMP,EY,1,8610
* J& P! `6 _( {# pMP,EZ,1,8610
MP,PRXY,1,0.3
MP,PRYZ,1,0.45
" ]* Q' `% R7 F' JMP,PRXZ,1,0.3
MP,GXY,1,4160
MP,GYZ,1,3000
) w9 M& ^! U' j/ n5 CMP,GXZ,1,4160

MP,EX,2,114000
+ W6 |8 T# j5 U; s; mMP,EY,2,8610*0.3
" G) ^7 O& z! A8 ?4 UMP,EZ,2,8610
$ u) `4 ]1 s7 {& AMP,PRXY,2,0.3
MP,PRXZ,2,0.3
3 s8 N# E: t6 \" p) U$ a6 b9 eMP,PRYZ,2,0.45
" Q8 c/ l' Q) @6 g; |MP,GXY,2,4160*0.3
MP,GXZ,2,4160
7 S2 J4 f; F, p' L: SMP,GYZ,2,3000*0.3
6 M2 K% I0 [( ]9 H/ Y
MP,EX,3,114000
  V1 z% B- \: w5 DMP,EY,3,8610*0.4
MP,EZ,3,8610
% I0 A# I4 J7 {MP,PRXY,3,0.3
MP,PRXZ,3,0.3
6 `! Z6 x3 Q. ^3 x' Z+ c) }MP,PRYZ,3,0.45
MP,GXY,3,4160*0.4
4 Q! O0 J7 U/ t5 S4 ^( m. N, D) s. \MP,GXZ,3,4160
MP,GYZ,3,3000*0.4
: [2 J4 H, q6 ?! q$ V# c8 Y
; X- S2 \' s' @% D2 f0 vMP,EX,4,114000*0.07
MP,EY,4,8610
# ]% h% t& G& b# O6 J( e  `. xMP,EZ,4,8610
MP,PRXY,4,0.3
MP,PRXZ,4,0.3
, x7 q- q/ ~! t& B9 U% @0 vMP,PRYZ,4,0.45
MP,GXY,4,4160*0.07
6 w, o4 Q% @, O' \; AMP,GXZ,4,4160*0.07
8 Q8 {# Y% d0 QMP,GYZ,4,3000

% u3 C/ a0 L& n) ^MP,EX,5,114000*0.14
MP,EY,5,8610
3 z; h+ q  {" ~+ KMP,EZ,5,8610
MP,PRXY,5,0.3
# P7 U1 n" y8 T: HMP,PRXZ,5,0.3
MP,PRYZ,5,0.45
4 h; i) b0 Q, SMP,GXY,5,4160*0.07
3 X  n7 W8 _' V# z" lMP,GXZ,5,4160*0.07
MP,GYZ,5,3000

MP,EX,6,114000
- z7 q. y8 r& k1 sMP,EY,6,8610*0.3
1 l5 S# l; [9 i" AMP,EZ,6,8610*0.3
- i/ }2 v8 J6 {# I4 Z7 IMP,PRXY,6,0.3
$ T4 |1 U+ [6 c3 ^MP,PRXZ,6,0.3
) C- S1 E% d2 i, z5 A8 MMP,PRYZ,6,0.45
MP,GXY,6,4160
) f, `$ N* S* h5 TMP,GXZ,6,4160*0.3
MP,GYZ,6,3000*0.3
" N  b* r4 X$ h) W' i
MP,EX,7,114000
MP,EY,7,8610*0.4
6 }- j* }3 f; NMP,EZ,7,8610*0.4
MP,PRXY,7,0.3
MP,PRXZ,7,0.3
MP,PRYZ,7,0.45
MP,GXY,7,4160
MP,GXZ,7,4160*0.4
MP,GYZ,7,3000*0.4
5 f  Y4 f% A% e+ s% o
MP,EX,8,114000*0.07
MP,EY,8,8610*0.3
MP,EZ,8,8610*0.3
MP,PRXY,8,0.3
, e! n/ c- h7 v$ M* v# k8 V* KMP,PRXZ,8,0.3
6 D* b, N0 T- L- ?6 K1 qMP,PRYZ,8,0.45
4 T# W% E* f: `+ GMP,GXY,8,4160*0.07
MP,GXZ,8,4160*0.07
MP,GYZ,8,3000*0.3

4 X# b, \9 N+ e!mp,ex,2,210000
$ E( {" e) h4 R' Z/ U!mp,prxy,2,0.3
% t- ^+ P9 T" o0 UXc=914.3
Xt=1229.4
% J2 k2 l2 m" g( u2 y3 W- {! N& MYc=147.2
1 x0 j- P0 U3 Z1 }  t/ TYt=45.6
* ^  Q, k9 k; pZc=147.2
Zt=45.6
S12=33.0
S13=33.0
4 v- u- }5 Y3 \! Q  o' y9 gS23=25.0
, q" @5 r7 p+ J- W*dim,ST11,array,800
*dim,ST22,array,800
*dim,ST33,array,800
6 k- q0 v% K3 H& ]4 n% m: r" t+ ?*dim,ST12,array,800
5 l/ V: N" S1 E*dim,ST13,array,800
*dim,ST23,array,800
*dim,DFT,array,800  
, M! E/ x# y# I8 i% d  y*dim,DFC,array,800
*dim,DMT,array,800
*dim,DMC,array,800
* s! x. ]( Q  g" @, G8 I7 o*dim,DFM,array,800
: k  c( n; p* j; J& p3 a*dim,DT,array,800
* W+ q1 c: O# b*dim,DC,array,800
*dim,erci,array,800
*dim,HASHINDFT,array,800
*dim,HASHINDFC,array,800      
) @0 g: P0 z, K$ E; q4 P*dim,HASHINDMT,array,800
3 F5 B; @6 C+ e! a+ Y; N' I( t*dim,HASHINDMC,array,800
*dim,HASHINDFM,array,800
7 [; C# f; {! @- m( _- F" e7 c1 Q1 w*dim,HASHINDT,array,800
*dim,HASHINDC,array,800
( P- S( p3 |7 p4 [3 t
' q) ^/ \2 C2 V$ qblock,0,100,0,50,0,0.5
6 D7 @2 j6 x) _" nLSEL,S,LENGTH,,100 !选择所有长度为90mm的线段(长度方向上的线段)
LESIZE,all, , ,20, , , , ,1
5 o7 q$ Y1 V8 tlsel,s,length,,50 !选择所有长度为25mm的线段(宽度方向上的线段)
lesize,all, , ,10,, , , ,1
LSEL,S,LENGTH,,100
; o0 P5 ]. @0 v$ k+ Mlsel,a,length,,50
( B  {1 ^: ~  ?: P7 Hlsel,inve
: e" C8 h3 Y2 l/ N) i+ alesize,all, , ,4
alls
) Q8 ?* u: Q9 N. {0 A. d3 D8 ~vatt,1,1,1
vmesh,all

  }9 ]* y4 F1 l2 @
' y9 j# {2 a& e+ ^7 w- l/solu
asel,s,loc,x,0
da,all,all
. e/ A6 b7 K7 I+ q, Y' dalls
asel,s,loc,x,100
sfa,all,1,press,-2
alls
+ l$ B3 h  ^' |, l!nsubst,5
!cnvtol,f,,0.1,,0.1
& w" w# ~& i2 S; a$ ^!rescontrol,define,all,1,1
solve
!save,1,db

. N/ e7 h% [% O: A; J5 c% y( B*do,i,2,10
  /post1
   etable,sx,s,x
( J3 D5 g$ J0 t5 f& n   etable,sy,s,y
& h; y2 ?% V! \) T! Q% }; y   etable,sz,s,z
  N  B# Y, E& c) R, v4 r& L   etable,sxy,S,xy
   etable,syz,s,yz
   etable,sxz,S,xz
   
   *do,j,1,800
      *get,ST11(j),elem,j,etab,sx
' c$ ~+ ?4 Z- J5 E- E      *get,ST22(j),elem,j,etab,sy
, n# _$ p' y4 U5 X      *get,ST33(j),elem,j,etab,sz
4 M1 T1 \5 m* p, s$ }( s3 b- L9 h      *get,ST12(j),elem,j,etab,sxy
      *get,ST13(j),elem,j,etab,sxz
  T: t5 n# Z6 B: e      *get,ST23(j),elem,j,etab,syz                  
   *enddo
! i3 z3 g: s. J/ X3 Y   *do,j,1,800
- y  }9 \* e- H- U& L4 U: F      HASHINDFT(j) = (ST11(j)/Xt)**2+(ST12(j)/S12)**2+(ST13(j)/S13)**2
% O( Q* E* E$ l1 l4 L) s      HASHINDFC(j) = (ST11(j)/Xc)**2
      HASHINDMT(j) = (ST22(j)/Yt)**2+(ST12(j)/S12)**2+(ST23(j)/S23)**2
7 Q; ?6 T) z* F, B4 P  N# L; E7 Z      HASHINDMC(j) = (ST22(j)/Yc)**2+(ST12(j)/S12)**2+(ST23(j)/S23)**2
: Q% a; {( c5 q4 c      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              
   *enddo   
& ]4 L3 z8 ?; J; u  R8 `& E parsav,all,parameter,txt
5 K, D. w2 O/ D! k2 d  f+ H; Z
& I- [- C! D, v  /solu
    antype,,rest
; t* T0 @' ?! C- M1 i% {+ uparres,new,parameter,txt
!失效类型
9 T) u  S7 X$ T!纤维失效
. Z9 \" J, `, l7 a) x6 H' N*do,j,1,800
4 g' C9 Z) w4 F0 b4 i4 N4 m      *IF,ST11(j),ge,0,and,HASHINDFT(j),ge,1,THEN
               DFT(j)=1+DFT(j)
8 p1 F6 p6 O& |5 |
3 h* F& h2 h# A% W            *ELSEif,ST11(j),lt,0,and,HASHINDFC(j),ge,1,then
4 f, F3 _9 n4 i               DFC(j)=1+DFC(j)
" P* D6 v& h  ^8 q' S
7 ~2 C& m2 Y8 {/ f/ C$ m5 Z3 N6 D9 u            *ELSE
                DFC(j)=0+DFC(j)
, @1 ^) ]5 b4 q, _. ?/ p                DFT(j)=0+DFT(j)
8 T- `7 d" [" r  E# r" u* C           
      *ENDIF
( C' y8 c5 y' X( H1 T  t*enddo
. ?' [$ A6 x: S! C( k* ]
. m$ q" U/ X6 o6 [+ y5 q* b!基体失效
5 U. v9 @- s# R# l9 z*do,j,1,800
) A3 i1 S& C. m& V: `6 N! f% l: [      *IF,ST22(j),ge,0,and,HASHINDMT(j),ge,1,THEN         
               DMT(j)=1+DMT(j)
9 \  G; K. P# X, ?3 O- e& D9 y; u         
          *ELSEIF,ST22(j),LT,0,and,HASHINDMC(j),ge,1,THEN
                DMC(j)=1+DMC(j)
" R5 h. ^: g& t2 h  }0 I4 v
# j7 C9 m8 M$ G. t* N& S! c          *ELSE
( I0 Y# q: V) t4 n  S9 g1 V/ Q                DMC(j)=0+DMC(j)
; i  R( o& t- w% p! s                DMT(j)=0+DMT(j)
+ B) g* |% Y8 U+ L! G! H8 k         
      *ENDIF
*enddo
  P) G- {8 q, w7 ~/ x! g% s!FIBER-MATRIX FAILURE
" N! g9 q; I, r1 h% z- j  ]*do,j,1,800
+ n/ h, r5 {# O; p  Q: I2 Y* {      *IF,ST11(j),ge,0,AND,HASHINDFM(j),ge,1,THEN
              DFM(j)=1+DFM(j)
7 L4 ]6 O" z1 `6 v
            *ELSE
3 j. g; O' x/ T2 p0 U              DFM(j)=0+DFM(j)
' n' i( `9 V3 }& R; X            
' X- d' W, k5 [+ a      *ENDIF
2 e7 A! h. N/ i2 r! G8 U' E$ T  M, O*enddo

) e/ q: c+ P4 t4 i6 y5 e& ^8 D!Tensile-Compression FAILURE
*do,j,1,800
      *IF,ST33(j),ge,0,AND,HASHINDT(j),ge,1,THEN
4 v* G' Q$ L4 A! C               DT(j)=1+DT(j)

          *ELSEIF,ST33(j),LT,0,AND,HASHINDC(j),ge,1,THEN
& l2 @- c" S& B7 C( n" [' x2 a. _                DC(j)=1+DC(j)
% z: p' F. n" `1 N$ P
          *ELSE
                DC(j)=0+DC(j)
! E, u# v% d2 j$ h* `% |                DT(j)=0+DT(j)
  {7 J8 r9 g5 v; D* ^         
      *ENDIF
*enddo
*do,j,1,800
    erci(j)=DFT(j)+DFC(j)+DMT(j)+DMC(j)+DFM(j)+DT(j)+DC(j)
* H. }# J- V; M/ [% ~6 D" w/ E*ENDDO
1 P$ L. u/ w: _; c# O, O( h* O/prep7
!更新单元材料
7 x+ L9 I! w/ {/ W' z*do,j,1,800
*IF,erci(j),ge,2,THEN
emodif,j,mat,8
     *elseIF ,DFT(j),gt,0,then            
: i. V, m$ v  Y5 Xemodif,j,mat,2
     
      *elseIF,DFC(j),gt,0,then
emodif,j,mat,3
            
      *elseIF,DMT(j),gt,0,then
emodif,j,mat,4
         
      *elseIF,DMC(j),gt,0,then
emodif,j,mat,5
/ G& n+ a+ o& c- K# Y            
( R% g) v; R# k: d2 S      *elseIF,DT(j),gt,0,then
$ u% }7 b# ^9 semodif,j,mat,6
            
      *elseIF,DC(j),gt,0,then
; T) O5 _, u! B! H" ~# F  i1 Memodif,j,mat,7
" Y* N, Q2 x" R2 g/ ^, w! m              
3 H8 X. R" W: c*endif      
*ENDDO
/solu
asel,s,loc,x,100
sfa,all,1,press,-i*200
$ @8 T6 r/ j/ q1 f+ k8 ~alls
!nsubst,5
!cnvtol,f,,0.1,,0.1
!ncnv,2
7 ?/ ^3 \, X3 e; P4 zsolve
% n5 L8 A6 d' k8 C0 ]" a
+ ^( B3 [9 |1 \3 B" s
; o1 k7 I$ _7 @3 ^6 U; V*enddo