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

kolarxiao

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

yinjvictor

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

limulin1

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

我想也是，看来你也是咯

dragon_w

我模仿您的apdl写了段，可还是结果无法提取，请教！
" z! s% E* G  c. i/ Vfini
, h) W; `& O4 B& D1 Q/clear
/prep7
et,1,46
keyopt,1,8,1
R,1,4
0 l& l  a( b; X) |9 c( eRMODIF,1,13,1,0, 0.125,1,45,0.125
RMODIF,1,19,1,-45,0.125,1,90,0.125
/ n/ J" R8 p- p- b! S. EMP,EX,1,114000
' G9 }4 C5 b5 b" f- g9 L0 iMP,EY,1,8610
  n  w* P+ j  C# p+ sMP,EZ,1,8610
5 {5 |) ]& b7 s7 m7 T8 t8 ?  a3 Y2 ZMP,PRXY,1,0.3
MP,PRYZ,1,0.45
MP,PRXZ,1,0.3
MP,GXY,1,4160
MP,GYZ,1,3000
MP,GXZ,1,4160

MP,EX,2,114000
MP,EY,2,8610*0.3
- v2 O; S* P1 y+ J/ LMP,EZ,2,8610
MP,PRXY,2,0.3
MP,PRXZ,2,0.3
% l) @: u% c% E1 ~; X0 a7 wMP,PRYZ,2,0.45
0 a1 Y8 n) a: ?  ]. a* w" s- V# k$ rMP,GXY,2,4160*0.3
MP,GXZ,2,4160
: n7 s( C2 C9 D6 N' O6 PMP,GYZ,2,3000*0.3

8 T- k! X2 ?3 w& G; v0 g( ZMP,EX,3,114000
3 [0 v% ~" ~% H+ uMP,EY,3,8610*0.4
1 G9 ]2 B9 ~; U  A* b6 W+ x# ~MP,EZ,3,8610
MP,PRXY,3,0.3
MP,PRXZ,3,0.3
MP,PRYZ,3,0.45
MP,GXY,3,4160*0.4
MP,GXZ,3,4160
- @' |5 Y8 `. d4 |0 v- YMP,GYZ,3,3000*0.4
5 V) q$ D% R( z4 q8 R. S, N: w; A
MP,EX,4,114000*0.07
+ m$ h. z# U7 G: [0 s* k+ CMP,EY,4,8610
- C2 _% x. `3 [6 n; i: O# QMP,EZ,4,8610
MP,PRXY,4,0.3
MP,PRXZ,4,0.3
8 J8 I" q0 R  b; CMP,PRYZ,4,0.45
- H; ^) V. {: R+ F8 K& A  tMP,GXY,4,4160*0.07
: L8 @  b$ F' o& L% OMP,GXZ,4,4160*0.07
2 U. N: h+ ^8 O+ K8 i7 QMP,GYZ,4,3000
; j  |7 r, r7 ^" e( Z4 q
MP,EX,5,114000*0.14
$ Y  Z# J9 k2 jMP,EY,5,8610
MP,EZ,5,8610
% v! ^& u0 b$ ~0 O$ f' lMP,PRXY,5,0.3
MP,PRXZ,5,0.3
MP,PRYZ,5,0.45
1 M6 `8 a& v5 B; f$ ]+ bMP,GXY,5,4160*0.07
MP,GXZ,5,4160*0.07
MP,GYZ,5,3000
: B- L2 w  m1 v. u7 l: l6 }0 z
MP,EX,6,114000
MP,EY,6,8610*0.3
MP,EZ,6,8610*0.3
MP,PRXY,6,0.3
MP,PRXZ,6,0.3
* T4 ^- c* ^* l! t) Y) }1 g# nMP,PRYZ,6,0.45
0 E) r, ~# F" n9 U! jMP,GXY,6,4160
2 w% M! `$ }2 h' k4 K; e# tMP,GXZ,6,4160*0.3
, R6 d; `) A) c9 y( nMP,GYZ,6,3000*0.3

MP,EX,7,114000
9 Z: D) Z4 _$ g0 G1 dMP,EY,7,8610*0.4
MP,EZ,7,8610*0.4
! b/ I: j  q9 t* JMP,PRXY,7,0.3
MP,PRXZ,7,0.3
0 F4 q4 B* C  ~, R5 G2 Z& XMP,PRYZ,7,0.45
% _" K, n; e% \0 R/ `  KMP,GXY,7,4160
MP,GXZ,7,4160*0.4
) S1 `5 X: K- ^9 v3 k/ U. B0 {MP,GYZ,7,3000*0.4

MP,EX,8,114000*0.07
MP,EY,8,8610*0.3
MP,EZ,8,8610*0.3
MP,PRXY,8,0.3
MP,PRXZ,8,0.3
, t- U' D+ d* t1 N; o% KMP,PRYZ,8,0.45
1 l5 C3 ^5 p, A6 a% [MP,GXY,8,4160*0.07
MP,GXZ,8,4160*0.07
MP,GYZ,8,3000*0.3

!mp,ex,2,210000
!mp,prxy,2,0.3
' u5 A+ I' k4 u( @2 CXc=914.3
Xt=1229.4
8 A4 [! l2 P9 f- ~$ o/ }- aYc=147.2
Yt=45.6
Zc=147.2
Zt=45.6
S12=33.0
S13=33.0
# d2 f. ?  B2 N- _* ]S23=25.0
*dim,ST11,array,800
" ~. H3 I3 ^' m3 z. Q8 h*dim,ST22,array,800
; ~* Y  X+ [( L- Y*dim,ST33,array,800
) `. h5 D% x6 T*dim,ST12,array,800
% p+ y2 \5 P% A6 b: j/ y& o! k. d*dim,ST13,array,800
- q& m" D2 j+ l7 @# d; A% o( V*dim,ST23,array,800
*dim,DFT,array,800  
*dim,DFC,array,800
. i0 J' r- I/ `: N) S5 q*dim,DMT,array,800
/ E4 A4 B! l( i*dim,DMC,array,800
4 V! W1 N5 h. i4 E, k8 ~# m*dim,DFM,array,800
+ w! M% E% K3 B/ N*dim,DT,array,800
*dim,DC,array,800
' b# J5 o8 r" s' c  P*dim,erci,array,800
*dim,HASHINDFT,array,800
*dim,HASHINDFC,array,800      
) X; x: V7 t: S8 q( Z1 A: J, m. v- `2 i*dim,HASHINDMT,array,800
*dim,HASHINDMC,array,800
3 ~7 K. x0 d- Y2 z6 i8 q3 _*dim,HASHINDFM,array,800
*dim,HASHINDT,array,800
*dim,HASHINDC,array,800

- H6 R9 s5 n# T- O0 }6 ~- Iblock,0,100,0,50,0,0.5
9 Y0 q0 k% ]+ N$ k/ uLSEL,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
- m: D7 J! \) E& @lsel,inve
lesize,all, , ,4
9 L5 e" `! @: B. j5 [8 q6 J! Zalls
vatt,1,1,1
5 J* S# L# b' P: I: mvmesh,all

0 H% U. ?0 \- ~; z4 @
/solu
6 E* ?( v* Y8 l+ M( p0 W, Qasel,s,loc,x,0
7 A3 ~6 @$ t; ~; Eda,all,all
( K) @3 Z% K; \  h" I2 Z8 salls
# L0 j, y: U1 _3 a& L) Iasel,s,loc,x,100
sfa,all,1,press,-2
  x3 p- @: P4 l% g# s1 ^alls
# r0 x) F' @9 a, s. T& `!nsubst,5
* a" O. [) |9 v$ r5 z4 o!cnvtol,f,,0.1,,0.1
!rescontrol,define,all,1,1
1 j. a9 Y  c0 Y! V  Ysolve
/ ~; I. ~5 V& d" P( x7 S# J!save,1,db
, Q+ \7 i8 ~+ N
*do,i,2,10
7 G" t/ N3 u0 h- H9 C  /post1
3 e% [% e5 C3 G) h2 ?   etable,sx,s,x
8 P- y9 v: O- k- h3 U: |. f   etable,sy,s,y
/ V  ~( W) J0 P% w! v   etable,sz,s,z
   etable,sxy,S,xy
   etable,syz,s,yz
   etable,sxz,S,xz
0 j$ K" ]) C2 F$ q0 w0 i& m4 b   
   *do,j,1,800
      *get,ST11(j),elem,j,etab,sx
      *get,ST22(j),elem,j,etab,sy
: F  R% w7 R9 b2 J3 Z) @      *get,ST33(j),elem,j,etab,sz
; g( E- g# E8 B3 J      *get,ST12(j),elem,j,etab,sxy
$ F; p* m3 B7 `( N: n1 E      *get,ST13(j),elem,j,etab,sxz
% V3 I5 ]* K; \4 D( Y# R  Z! `) i      *get,ST23(j),elem,j,etab,syz                  
   *enddo
- `' U  z% g! o+ ]3 H+ s$ g   *do,j,1,800
6 b+ S/ k: I; W+ [1 n- \& ^: D  P      HASHINDFT(j) = (ST11(j)/Xt)**2+(ST12(j)/S12)**2+(ST13(j)/S13)**2
      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
5 {2 m; s. T2 Q% K      HASHINDT(j) = (ST33(j)/Zt)**2+(ST13(j)/S13)**2+(ST23(j)/S23)**2
" Y: h' R2 }* `      HASHINDC(j) = (ST33(j)/Zc)**2+(ST13(j)/S13)**2+(ST23(j)/S23)**2              
   *enddo   
parsav,all,parameter,txt

: Y- L9 E1 r# \5 e8 g  /solu
    antype,,rest
parres,new,parameter,txt
% X3 x8 V) U0 I: ~% r/ p# R3 ^8 e!失效类型
- ~  l1 H: n, P!纤维失效
8 }/ i1 \7 r% O' y  j*do,j,1,800
/ m/ C+ `+ o- n9 }! @7 G      *IF,ST11(j),ge,0,and,HASHINDFT(j),ge,1,THEN
( X2 o  y* R+ ^2 X5 \               DFT(j)=1+DFT(j)

            *ELSEif,ST11(j),lt,0,and,HASHINDFC(j),ge,1,then
$ x! ?# I$ V+ w, O! c2 s               DFC(j)=1+DFC(j)
8 `3 `, X- C0 U3 q
. z, u) C9 A& x+ p            *ELSE
                DFC(j)=0+DFC(j)
, e, _: K# {3 E9 J, ?7 R                DFT(j)=0+DFT(j)
) a# ]8 X+ d2 k/ f) Z; l           
      *ENDIF
2 J/ Y5 z' W: q- h6 s. j) B6 t*enddo
$ I& z6 W$ B" y
0 [/ N6 O+ ^7 M6 k" _!基体失效
*do,j,1,800
) e8 E0 T; @; t6 W' ~2 p) M' n) D      *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)

          *ELSE
4 T: r& y; ~6 q$ {                DMC(j)=0+DMC(j)
                DMT(j)=0+DMT(j)
: I, T9 }$ w9 I1 e4 J- i7 K; R         
      *ENDIF
*enddo
2 Q8 c  X3 Q# [!FIBER-MATRIX FAILURE
: V" p4 F% d, w4 V: C( }*do,j,1,800
0 I6 _# ]/ w/ L; I3 D1 R7 u. C# a      *IF,ST11(j),ge,0,AND,HASHINDFM(j),ge,1,THEN
              DFM(j)=1+DFM(j)
) }9 ?: G# m) A+ K0 e
: w1 f! b& b4 [- ?* ^1 l2 w            *ELSE
              DFM(j)=0+DFM(j)
            
      *ENDIF
0 e. O4 w$ D# Q1 B3 y*enddo

!Tensile-Compression FAILURE
*do,j,1,800
. J  M' F; H/ B* R4 t7 [      *IF,ST33(j),ge,0,AND,HASHINDT(j),ge,1,THEN
2 n8 m  |# Z! B* M               DT(j)=1+DT(j)

% T* ~& r* T7 z1 K2 ]          *ELSEIF,ST33(j),LT,0,AND,HASHINDC(j),ge,1,THEN
                DC(j)=1+DC(j)
5 [, `: B0 x. ^4 ]" p' @* a2 o, K
# e0 o* D9 D6 F4 j  A! s          *ELSE
. x7 o0 }" ?6 B                DC(j)=0+DC(j)
" Y7 a2 V; f4 d" s/ r; \                DT(j)=0+DT(j)
         
  T. D; Z. E) \7 E0 j6 K  Q9 S9 t      *ENDIF
*enddo
*do,j,1,800
    erci(j)=DFT(j)+DFC(j)+DMT(j)+DMC(j)+DFM(j)+DT(j)+DC(j)
*ENDDO
2 [$ E% d9 D: N+ D( S/prep7
!更新单元材料
*do,j,1,800
6 S) O) e+ T  `1 g  d" C' h2 Y*IF,erci(j),ge,2,THEN
3 e0 Y  L; W" @. K" e3 Lemodif,j,mat,8
     *elseIF ,DFT(j),gt,0,then            
emodif,j,mat,2
- g& }! t$ a6 }# }$ O# S  `     
; O- [7 H  w) V  }8 d, g      *elseIF,DFC(j),gt,0,then
. M' j; e5 \  a/ Y7 n' F. m! Yemodif,j,mat,3
            
      *elseIF,DMT(j),gt,0,then
emodif,j,mat,4
         
      *elseIF,DMC(j),gt,0,then
emodif,j,mat,5
            
      *elseIF,DT(j),gt,0,then
emodif,j,mat,6
            
- P* `4 C8 |" Y      *elseIF,DC(j),gt,0,then
emodif,j,mat,7
, ~5 o1 _0 a- E$ e$ K% [              
+ a) h1 W; G. r2 d+ m*endif      
% ^8 @$ G3 l4 R*ENDDO
0 W) w# M+ Q. l) ]: z- q0 U/solu
+ V& d9 G1 Q. ]: H" s, Gasel,s,loc,x,100
sfa,all,1,press,-i*200
6 U3 z& l# L4 m) {: Ealls
!nsubst,5
0 V, o5 Y- @# n5 Y' U, E- w. O$ U!cnvtol,f,,0.1,,0.1
!ncnv,2
! D0 j+ L  X( s5 r1 O8 c& ?. msolve
! ^" m5 r; a$ z  o5 S- q1 F

! I2 r5 N' h' B- G*enddo
! y7 U1 j, W. Q, s+ `