找回密码
 注册
Simdroid-非首页
楼主: amani

[学术与经验交流] 发一篇我用ANSYS-APDL做的复合材料圆柱层合板渐进失效的国际杂志论文和程序

[复制链接]
发表于 2011-6-21 09:18:40 | 显示全部楼层 来自 清华大学
楼主,您好!您这里的渐进失效和疲劳失效是什么关系啊?
回复 不支持

使用道具 举报

发表于 2012-1-5 08:34:53 | 显示全部楼层 来自 美国
Simdroid开发平台
这个文章怎么都是你们自己在引用,没有什么外国人引用啊!
回复 不支持

使用道具 举报

发表于 2012-12-30 23:45:42 | 显示全部楼层 来自 浙江杭州
sxyuyifan 发表于 2010-2-4 08:26 / M& i; u9 J/ Y  f6 f, m2 ~4 h
楼主真强,是浙大化机所的啊。
' R, f0 s8 L% a# V4 |
我想也是,看来你也是咯
回复 不支持

使用道具 举报

发表于 2014-6-23 16:26:05 | 显示全部楼层 来自 江苏南京
我模仿您的apdl写了段,可还是结果无法提取,请教!
" z! s% E* G  c. i/ Vfini
, h) W; `& O4 B& D1 Q/clear6 Z5 h( k( y% Z- L- a* z, S* Q) f
/prep7$ ?& k' Q* V6 K: N3 v- c
et,1,460 q( T# E4 i" F- y" r
keyopt,1,8,1" U- ]" u; q' C5 U5 x8 [/ t
R,1,4
0 l& l  a( b; X) |9 c( eRMODIF,1,13,1,0, 0.125,1,45,0.1252 o. N7 ?( T+ }- W- l3 @
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 ' W3 X0 D% h  i2 T. m* V5 P# T
MP,PRYZ,1,0.45 / i" @- `) B( K, r2 y* p2 y
MP,PRXZ,1,0.3: S1 K! ~: W. e# F0 A* E0 z
MP,GXY,1,4160 1 z% O- s: n3 |$ M4 g% c! S- z9 Y
MP,GYZ,1,3000 ) c8 ?) D& o/ I
MP,GXZ,1,4160 # X1 Q, O% F8 x; V7 [  q
! j0 K8 T+ M0 t( p* N+ J
MP,EX,2,114000 & d( n3 _& q. h' _
MP,EY,2,8610*0.3
- v2 O; S* P1 y+ J/ LMP,EZ,2,8610 " ^# d4 R' @5 d7 e+ }: Q2 g
MP,PRXY,2,0.3 / q  z- w, V- l" f# P' l- }
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.38 c0 p2 C, y% l* G' }8 W6 z
MP,GXZ,2,4160
: n7 s( C2 C9 D6 N' O6 PMP,GYZ,2,3000*0.39 O' g  u1 }5 z9 Y' X! p/ [

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 & s- i6 R7 i/ N9 @- P
MP,PRXY,3,0.3 - g5 h% m& Z7 @; z/ S* O  p
MP,PRXZ,3,0.3% g: ^- z, M8 `
MP,PRYZ,3,0.45 4 D4 q5 I/ q* s0 W5 w) u
MP,GXY,3,4160*0.4* |; i8 h7 U+ t# i0 \
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% t! W3 F2 k4 ~6 d
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 " Z. Q& f/ L) v# I# v+ D
MP,PRXY,4,0.3 5 [* n. D+ c. Z# Q# t7 b
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/ Q5 P2 {( U6 a/ P
MP,EX,5,114000*0.14
$ Y  Z# J9 k2 jMP,EY,5,8610 6 L3 k' k% R0 [& C; u
MP,EZ,5,8610
% v! ^& u0 b$ ~0 O$ f' lMP,PRXY,5,0.3 0 g- @4 P- G& m/ C' g
MP,PRXZ,5,0.3/ N- o6 \8 d5 M* _
MP,PRYZ,5,0.45
1 M6 `8 a& v5 B; f$ ]+ bMP,GXY,5,4160*0.072 R- {5 o+ c4 j3 ^/ b  {
MP,GXZ,5,4160*0.07 - J  P! [' w6 |2 R" j
MP,GYZ,5,3000
: B- L2 w  m1 v. u7 l: l6 }0 z  n$ c; n; ^2 s# l
MP,EX,6,114000 * M( X5 |, O; N$ s6 Q+ V
MP,EY,6,8610*0.3   I' u3 P5 I8 J1 A
MP,EZ,6,8610*0.3 & F) m- b4 _& O4 ^, V+ d
MP,PRXY,6,0.3 # j4 |8 J% T" s4 V; z
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  [, d+ e/ |3 V) \) I
5 ?( I- n" b& t: e# X7 Y
MP,EX,7,114000
9 Z: D) Z4 _$ g0 G1 dMP,EY,7,8610*0.4 ! H" h* }6 g8 H% ~+ E4 t
MP,EZ,7,8610*0.4
! b/ I: j  q9 t* JMP,PRXY,7,0.3 ( B5 i5 {4 M( @
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! x0 W/ U& F1 x
MP,GXZ,7,4160*0.4
) S1 `5 X: K- ^9 v3 k/ U. B0 {MP,GYZ,7,3000*0.43 l. _3 e5 P( X. M+ X5 i3 ^
6 _1 t1 }5 N; M3 s! K
MP,EX,8,114000*0.07 3 b) O; {- w+ }1 }, u" h9 F
MP,EY,8,8610*0.3 ' C" j8 Z5 e. D/ {: B7 R- [
MP,EZ,8,8610*0.3 ' I! K3 g+ e- n$ \$ m( B# ^$ t
MP,PRXY,8,0.3 $ [  B# x$ I5 q7 Y1 Q
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+ V- t3 v4 E9 v. j8 p3 R( i6 @5 h
MP,GXZ,8,4160*0.07 6 r+ F( R0 m% q* O  [1 D2 q3 W
MP,GYZ,8,3000*0.3, k0 {- R: G& ]! {
, N+ V  T* U* @
!mp,ex,2,210000! R3 R7 }; x' P# `0 L
!mp,prxy,2,0.3
' u5 A+ I' k4 u( @2 CXc=914.3/ X* A& L9 M/ K5 r( S
Xt=1229.4
8 A4 [! l2 P9 f- ~$ o/ }- aYc=147.2' \0 Y2 H+ A  h& A  L) m5 ]/ @
Yt=45.6% h0 c/ @6 ]3 H  L$ s0 g
Zc=147.2- d2 ~8 M; ~; Q9 Z
Zt=45.6# W( @7 w% ^8 z6 d
S12=33.0# D- n1 v, Z7 f, z. p- y
S13=33.0
# d2 f. ?  B2 N- _* ]S23=25.0; m% a* K0 T* q( ~0 l
*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! g8 w/ |+ E  n- `5 T
*dim,DFT,array,800  6 ^# [% U; z4 t% O3 R
*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,8006 v, m( k! R( b' u0 Y4 y
*dim,DC,array,800
' b# J5 o8 r" s' c  P*dim,erci,array,8002 D2 j6 T6 o# M: {7 ~8 H! e( ^
*dim,HASHINDFT,array,800. B) m# V$ W) c. d
*dim,HASHINDFC,array,800      
) X; x: V7 t: S8 q( Z1 A: J, m. v- `2 i*dim,HASHINDMT,array,8008 U' A. \* Q* J
*dim,HASHINDMC,array,800
3 ~7 K. x0 d- Y2 z6 i8 q3 _*dim,HASHINDFM,array,800& o# i/ V9 K% a4 U& @( |1 V
*dim,HASHINDT,array,800! j- C* e3 ~. H( s% M9 ?
*dim,HASHINDC,array,8004 p2 o4 o9 a  r( E1 m3 Q

- 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的线段(长度方向上的线段)# H# v. P7 k0 M* W( j# x
LESIZE,all, , ,20, , , , ,1 ; l8 j( o. Q2 s# y* ?, G6 W9 g
lsel,s,length,,50 !选择所有长度为25mm的线段(宽度方向上的线段)  @- e7 h# v2 A3 ~$ o% O. }
lesize,all, , ,10,, , , ,1 & s3 w( C: d; Y$ I* B# u  f
LSEL,S,LENGTH,,1003 Z) }& }5 e; |1 ]* K
lsel,a,length,,50
- m: D7 J! \) E& @lsel,inve6 C  Q& U& h$ C. C8 D
lesize,all, , ,4
9 L5 e" `! @: B. j5 [8 q6 J! Zalls! k# j( y+ Y" w- {" V
vatt,1,1,1
5 J* S# L# b' P: I: mvmesh,all; g* o5 ~  _  B: ?

0 H% U. ?0 \- ~; z4 @- v3 M1 i8 d3 {/ l& q; W7 `8 L
/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' n; D+ v6 B+ e: L: x, r3 k% _' _5 i
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/ W* ~: a: I; T0 {5 l* P3 `
!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# G; o# [* Q$ e. \$ F: e: ~
*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' S& p3 l* f3 E# [9 f/ O
   etable,sxy,S,xy' P- J5 G( g( j3 @/ Z
   etable,syz,s,yz: ~( }" T$ }1 L# G2 t9 r: Q0 j& ~! j
   etable,sxz,S,xz
0 j$ K" ]) C2 F$ q0 w0 i& m4 b   ( V9 U: y/ x' x# n: i0 \" k' V
   *do,j,1,800- ?6 A$ v( C2 X7 d6 }
      *get,ST11(j),elem,j,etab,sx' [5 u( ^3 j4 v4 \
      *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                  ! K( o" [' j1 I1 ?9 r" ]1 |
   *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- v$ P, O% v* m- r. Y0 U5 m/ Y7 {
      HASHINDFC(j) = (ST11(j)/Xc)**2( `: `1 @$ H& I7 P4 F# t1 \6 C
      HASHINDMT(j) = (ST22(j)/Yt)**2+(ST12(j)/S12)**2+(ST23(j)/S23)**2/ m/ \" }9 Z: D% Y3 g) s: Z: X3 Y$ p
      HASHINDMC(j) = (ST22(j)/Yc)**2+(ST12(j)/S12)**2+(ST23(j)/S23)**2" {  ?* U6 G& P7 a" {( ?
      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              " g* }- y6 ?- e$ T
   *enddo   " |* ^- K1 o/ E/ ?8 L' z+ z" p3 B
parsav,all,parameter,txt& h- h5 l+ {# x' A5 @+ D) o4 E

: Y- L9 E1 r# \5 e8 g  /solu. C% H& ?& w: a# o& V( {
    antype,,rest& n) w; M0 n4 R2 l0 `
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): I2 R2 T2 q6 _: G4 u* C" v
1 o; x$ U2 ?$ b$ f
            *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' b( e5 _" Z4 n. B: m
                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           / ~- B: u# \$ q" ~
      *ENDIF
2 J/ Y5 z' W: q- h6 s. j) B6 t*enddo
$ I& z6 W$ B" y
0 [/ N6 O+ ^7 M6 k" _!基体失效4 ]$ ]$ ~  K) P5 f5 q3 T" ^
*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         , ~  t7 b: t6 r' |* [: R( u
               DMT(j)=1+DMT(j)7 B4 u6 K" x6 P& W9 r
         $ B6 Z3 q. T% y% [! r, n) T4 ^8 J
          *ELSEIF,ST22(j),LT,0,and,HASHINDMC(j),ge,1,THEN' ]  M9 ]2 d1 v' n' t& R/ x& ?
                DMC(j)=1+DMC(j)9 W/ y; J1 p. M: j! a
9 R! H, k- p* F! l* W2 i
          *ELSE
4 T: r& y; ~6 q$ {                DMC(j)=0+DMC(j); x" L  E: C$ x# |
                DMT(j)=0+DMT(j)
: I, T9 }$ w9 I1 e4 J- i7 K; R          9 i: ~/ m, X5 a* |8 }( z0 z
      *ENDIF* {0 {9 r8 q7 A' R! Z3 T% S
*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- ~! w  d4 N/ ~/ O: v9 A& M
              DFM(j)=1+DFM(j)
) }9 ?: G# m) A+ K0 e
: w1 f! b& b4 [- ?* ^1 l2 w            *ELSE5 T8 D) u, W% Y, l5 w  r  K% W9 s' Y! a' n
              DFM(j)=0+DFM(j)! A& F) C! x! P; v  F- I5 [4 P
            7 L1 c8 X1 H1 ^
      *ENDIF
0 e. O4 w$ D# Q1 B3 y*enddo8 V% G4 |+ V7 r
9 x  @& J& \9 ?$ i' h* x
!Tensile-Compression FAILURE' d/ n( t; M$ g
*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)* ~! f: ?) [  m8 j5 ]0 l; h

% T* ~& r* T7 z1 K2 ]          *ELSEIF,ST33(j),LT,0,AND,HASHINDC(j),ge,1,THEN6 [- G5 _# F- X/ e0 n! h
                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  n+ T4 F; b5 t
         
  T. D; Z. E) \7 E0 j6 K  Q9 S9 t      *ENDIF2 D. o+ u, B3 P( C* `
*enddo: W- g! C6 R8 j; m7 ^7 V6 U' y
*do,j,1,8007 B( Z, h, w4 N; |
    erci(j)=DFT(j)+DFC(j)+DMT(j)+DMC(j)+DFM(j)+DT(j)+DC(j)+ I9 \2 \* M  O; E
*ENDDO
2 [$ E% d9 D: N+ D( S/prep7" f0 S# B) M# t5 \
!更新单元材料, K+ s! @( n! w9 s- Y( E
*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% c0 i- m  _. X5 ^
     *elseIF ,DFT(j),gt,0,then            1 W; z0 O2 `0 y, S' w4 ~. V
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  w" Z3 p  a0 A9 Q/ G, o' C
            % a& z+ f3 q. U, E! D6 F+ @
      *elseIF,DMT(j),gt,0,then+ W0 D! Z( y6 n  X, N
emodif,j,mat,4* _* d1 r9 p  ^) g) y- b. ~! }5 W
          7 @! m  H& J: d: D
      *elseIF,DMC(j),gt,0,then8 f* c4 ?- l) _  Y
emodif,j,mat,5" F9 P7 U) X* f3 x6 v
             # S1 N2 k2 W% u5 J* i
      *elseIF,DT(j),gt,0,then7 p, ^4 `# X+ i4 O6 Z) P# O: d
emodif,j,mat,69 m( B+ v4 h( c) @5 i/ H- C
            
- P* `4 C8 |" Y      *elseIF,DC(j),gt,0,then  N- K$ S5 b1 I& R) p* e: K4 B* h4 R
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: M# M/ ]6 T9 @% J! y  x# S
sfa,all,1,press,-i*200
6 U3 z& l# L4 m) {: Ealls8 B' n2 B1 E8 D. ~; b. q
!nsubst,5
0 V, o5 Y- @# n5 Y' U, E- w. O$ U!cnvtol,f,,0.1,,0.11 W5 ~% _' q& p% _$ b  u4 x$ D  \
!ncnv,2
! D0 j+ L  X( s5 r1 O8 c& ?. msolve
! ^" m5 r; a$ z  o5 S- q1 F8 \% J6 o( r5 S- B, q3 B' w. A2 H

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

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2026-1-10 18:43 , Processed in 0.024756 second(s), 8 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2025 Discuz! Team.

快速回复 返回顶部 返回列表