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

amani

论文：Progressive failure analysis of carbon fiber/epoxy composite laminates using continuum damage mechanics
& [4 m0 a: x" W  u6 w应用背景：氢燃料电池汽车轻质高压储氢容器
8 n) m+ J) x; H) B, D( ]8 j单位：浙江大学化工系
3 O* ]: _( ]5 a. w0 mQQ:56691747
国际杂志：MATERIALS SCIENCE AND ENGINEERING:A(影响因子1.8)
现附上程序，供大家参考，请多引用哟。
1 Z" E% f4 {1 I& e9 T* v0 T/title, pprogressive failure analysis of carbon fiber/epoxy composite laminates using CDM(40MPa hydrogen storage tank)
) S  E6 C) h8 K' K/uis,msgpop,3
/prep7
et,1,95
et,2,64
mp,ex,1,70e3
% [; T! L3 k! i( }mp,nuxy,1,0.3
tb,biso,1,1,2
tbdata,,246,600
2 S6 {) n6 i. y2 R+ h6 n*afun,deg
E1=181e3
  A4 r' Y$ I! r! E/ O, P) @4 cE2=10.3e3
+ C6 E: M) `+ a8 dv12=0.28
5 _& W: V6 Y1 W- Q6 U6 Ov23=0.49
4 g6 y5 S$ x2 R5 \! O*dim,ls,array,1000
*dim,offsy,array,5000
*dim,offsz,array,5000
*dim,offsyz,array,5000
*dim,offstrainz,array,5000
*dim,offstrainsita,array,5000
5 N, _, U$ P! d9 K* ? *dim,offstrainzsita,array,5000
5 r, v0 a, {' `, }; D- }5 X *dim,ons1,array,5000
*dim,ons2,array,5000
$ Q2 K& H3 f$ x& e/ G *dim,ons6,array,5000
- Z/ ^; m+ _' V( U" ]*dim,onstrain1,array,5000
*dim,onstrain2,array,5000
+ m: l" w3 |! z. S *dim,onstrain6,array,5000
*dim,H1,array,5000
*dim,H2,array,5000
*dim,H6,array,5000
*dim,H12,array,5000
*dim,Tsaiwu,array,5000
6 _% j. i- s9 C& S( K+ H *dim,Y1,array,5000
' ?  z; Y4 U( b( m* a *dim,Y2,array,5000
*dim,Y6,array,5000
F22=1.0/(298*298)
& m! v8 q  d/ x. s' d% B F66=1.0/(778*778)
( G+ Q/ \$ b) N# Z3 z F11=1.0/(2150*2150)
; c* _' `( L3 ~! F F12= -0.5*sqrt(F11*F22)
! q- C2 j- m* h6 a' I' @ LFN=0
FFN=0
MCN=0
SSN=0
* ^0 ]1 ~; n: N1 T1 w6 Y*dim,MLFN,array,100
*dim,MFFN,array,100
* s9 O6 a, c. X$ K*dim,MMCN,array,100
*dim,MSSN,array,100
*dim,D1,array,5000
0 T& ~9 k5 G( H+ T* x9 L6 g*dim,D2,array,5000
8 s& V+ U5 m/ e' ~' ^3 s6 ]* n*dim,D6,array,5000
*dim,angle,array,10
*dim,m,array,10
0 C' y2 ]+ b* I*dim,n,array,10
& Y! x! c  o. K0 r*dim,te,array,6,6,5000
, |' j) ^5 T. _0 |+ h7 C4 Y  K*dim,ts,array,6,6,5000
6 M; A% w" m& [0 v4 M! _*dim,tst,array,6,6,5000
*dim,unit,array,6,6,5000
*dim,s,array,6,6,5000
" o2 u- }3 I7 O# Q- R7 `1 w*dim,c,array,6,6,5000
*dim,cp,array,6,6,5000
*dim,midlpf,array,6,6,5000
! K% z; d/ J7 d0 j5 Y3 h" h! U* u9 C
angle(1)=90
angle(2)=-90
7 |9 |3 d+ h, z) x7 u2 `angle(3)=18.9
3 t% O0 C# U# m  sangle(4)=-18.9
angle(5)=90
angle(6)=-90
# R$ ^& ?, q' l; H2 ]angle(7)=28.6
angle(8)=-28.6
angle(9)=90
* \$ h& c5 T! m& P- x4 eangle(10)=-90
8 \& H  w  B+ @1 l
6 c! T8 P* \3 G$ N% y( l" D*do,i,1,10
1 c! Q$ g& l. z5 n   *vfun,m(i),cos,angle(i)
  h$ K, {9 U$ |5 S% O* ~; }& L9 ^   *vfun,n(i),sin,angle(i)
$ n& j7 A3 \1 ~8 O2 Q5 w1 }*enddo
1 o8 n: P. K2 u0 D0 Y
*do,i,1,10
   *do,j,(i-1)*500+1,i*500
      te(1,1,j)=m(i)*m(i),n(i)*n(i),0,0,0,-2*m(i)*n(i)
; m, R2 o% t) t7 x& Z* r& `      te(1,2,j)=n(i)*n(i),m(i)*m(i),0,0,0,2*m(i)*n(i)
* |. {% |5 ?0 \' {" O& a      te(1,3,j)=0,0,1,0,0,0
      te(1,4,j)=0,0,0,m(i),n(i),0
      te(1,5,j)=0,0,0,-n(i),m(i),0
; i% G# u" B  f      te(1,6,j)=m(i)*n(i),-m(i)*n(i),0,0,0,m(i)*m(i)-n(i)*n(i)
   *enddo
( v! ~/ V# \9 L, |  F* ?  R2 [*enddo

  B& Z- ~, X5 n8 x. ~; ]* M*do,i,1,5000
   unit(1,1,i)=1,0,0,0,0,0
   unit(1,2,i)=0,1,0,0,0,0
* I* `9 ~  P: q3 n' k# ~   unit(1,3,i)=0,0,1,0,0,0
% {2 d9 D9 V" Y' ]   unit(1,4,i)=0,0,0,1,0,0
# [1 q) z9 g* \: Z3 K8 j* d   unit(1,5,i)=0,0,0,0,1,0
   unit(1,6,i)=0,0,0,0,0,1
*enddo

) a3 \7 y( H) S' x' v% u: x*do,i,1,5000
   *moper,tst(1,1,i),te(1,1,i),solv,unit(1,1,i)
*enddo

* [+ E+ A! E. c1 I6 q% X* T# L*do,i,1,5000
/ D( A$ M' _9 |+ r# }* E( R. ~   *mfun,ts(1,1,i),tran,tst(1,1,i)
" o0 i" h4 \( y" E& {*enddo

*do,i,1,5000
, N* c2 C: `8 }6 _/ A) Qs(1,1,i)=1.0/E1,-v12/E1,-v12/E1,0,0,0
; s6 \  Q0 m( ?5 \' ?+ q# B$ w        s(1,2,i)=-v12/E1,1.0/E2,-v23/E1,0,0,0
) K; s3 r- t* t; J9 @4 E; x        s(1,3,i)=-v12/E1,-v23/E1,1.0/E2,0,0,0
, S, w1 E" R7 z  ?( E        s(1,4,i)=0,0,0,2*(1+v23)/E2,0,0
        s(1,5,i)=0,0,0,0,2*(1+v12)/E1,0
2 y# F3 }1 \4 c$ w2 Y        s(1,6,i)=0,0,0,0,0,2*(1+v12)/E1
$ T1 f9 q5 F6 \/ m; [& C$ r8 j  *enddo

) {4 j' z3 Q# w9 O5 B" Y6 ? *do,i,1,5000
     *moper,c(1,1,i),s(1,1,i),solv,unit(1,1,i)
*enddo

. W' z$ w7 q5 a/ [0 w; V$ V*create, material
7 X* c0 g. M! \   *do,j,1,5000
c(1,1,j)=(1-D1(j))*c(1,1,j)
c(2,2,j)=(1-D2(j))*c(2,2,j)
% ]' _# C( [0 z/ O) W: t: kc(6,6,j)=(1-D6(j))*c(6,6,j)
*enddo

5 G5 _. p+ O6 ^- A( v; J   *do,j,1,5000
. u7 u) x7 I. D# H7 t      *moper,midlpf(1,1,j),ts(1,1,j),mult,c(1,1,j)
   *enddo

! U8 l) C5 l: r2 A6 r4 N( A1 P; p0 K1 g    *do,j,1,5000
       *moper,cp(1,1,j),midlpf(1,1,j),mult,tst(1,1,j)
*enddo
' C0 F% j, H, J  r/ L' r
*do,j,1,5000
TB,ANEL,j+1,1,21,0
TBDATA,,cp(1,1,j),cp(1,2,j),cp(1,3,j),0,0,cp(1,6,j)
. F+ Z" k- J! [. x$ d1 \TBDATA,,cp(2,2,j),cp(2,3,j),0,0,cp(2,6,j),cp(3,3,j)
TBDATA,,0,0,cp(3,6,j),cp(4,4,j),cp(4,5,j),0
* Q/ s( O# T! ?7 hTBDATA,,cp(5,5,j),0,cp(6,6,j)
; @& V* O& s2 M/ U+ ~3 y: X    *enddo
; n1 v% N6 R7 y9 U1 f*end
2 P! v7 @/ y  e/ X4 l*use,material

csys,1
cyl4,0,0,44,0,45.8,30,160
- @3 ?# s! Q# E* k$ _4 J. p*do,i,1,10
4 l- v; z, A( f/ W   cyl4,0,0,45.8+(i-1)*0.42,0,45.8+i*0.42,30,160
$ f' V! A0 J* b*enddo
vglue,all
numcmp,all
, G- S6 K, Q0 X( u2 b9 U8 ?LESIZE,12,,,50
LESIZE,1,,,2
0 \. c2 _" w1 c* i: r; TLESIZE,3,,,2
LESIZE,6,,,2
' ^. ~+ t* B( V4 }( JLESIZE,8,,,2
) V+ A! p, V4 VLESIZE,2,,,10
7 x+ v4 ]5 Y& A6 qLESIZE,4,,,10
$ S" H( @1 B2 b( [LESIZE,5,,,10
LESIZE,7,,,10
  x+ [+ Q( F# _  Y
1 x6 i# z7 S( [9 v( ~' b% O*do,i,1,9
LESIZE,13+(i-1)*4,,,10
/ i; F$ F# @) K0 f9 o*enddo

, m# f) G! c+ X# y) c4 G1 [$ W*do,i,1,9
LESIZE,14+(i-1)*4,,,10
*enddo
1 `; V. v$ P$ @. G
*do,i,1,9
! v1 {6 y* A$ \% W  LESIZE,15+(i-1)*4,,,50
# K% d% L9 U1 [; t5 ~*enddo
/ ?9 G- ]( p; {2 v: s
5 s. z( y/ j  ~& }*do,i,1,9
  LESIZE,16+(i-1)*4,,,50
/ q, J/ Q5 ?* l$ l*enddo

*do,i,1,10
  J. H; @6 r4 V% T7 ~2 e/ `7 L7 B! Z  LESIZE,55+(i-1)*4,,,1
*enddo

9 r. S! J& u* H' z# ~*do,i,1,10
* `* i* e+ F* v1 E9 z LESIZE,56+(i-1)*4,,,1
" @& t( f" j4 {9 }! |*enddo

alls
type,1
mat,1
/ Y. @$ V/ ], mvmesh,1
; u" \2 f! u. a) n" [7 T/ ?type,2
& s2 f8 F% t' l, |*do,i,1,10
  mat,i+1
" X1 T8 Z6 d2 `6 v# ]- ?3 H  vmesh,i+1
2 g" M# \9 D  n; ~; l*enddo
nummrg,node,,,,low
  j+ ~1 |& R. `' \nummrg,elem,,,,low
1 v2 Y& q4 T! K*do,i,1,5000
   emodif,i+1000,mat,i+1
*enddo
/pnum,mat,1
eplot
) d7 b& \* o. z' q: C
/solu
da,5,symm
da,6,symm
* y6 v& d" _' h1 A*do,i,1,10
  da,(i-1)*4+19,symm
. l4 _, g  _+ ]*enddo

% K% B1 c0 t; D* x*do,i,1,10
  da,(i-1)*4+20,symm
*enddo
d,1000,uz,0
rescontrol,define,all,1,1
sfa,4,1,pres,30
, K2 U- r, p/ V/ \% o. U- I- ?nsubst,5
cnvtol,f,,0.1,,0.1
solve

9 Q5 Q4 K! \1 }- Y+ p$ i*do,i,2,92
" }6 r, F5 A/ K" ^1 p9 @  /post1
# H3 A5 u# y: ]   rsys,1
etable,sequ,s,eqv
5 Z& u* A" M9 M0 f' I' Xetable,sy,s,y
   etable,sz,s,z
   etable,syz,s,yz
+ c, ]' Q/ x- j2 T   etable,strainy,epel,y
   etable,strainz,epel,z
& B/ _0 Z. c1 b5 X/ M% `" O+ K   etable,strainyz,epel,yz
9 x; d. l" _  D4 \) _9 J
   *do,j,1,1000
     *get,ls(j),elem,j,etab,sequ
! @3 D' J: D: ^5 z  J) H) e* x) W   *enddo
# d0 s) k/ [1 j& T5 ]
   *do,j,1,5000
/ [, \  ?( j! c      *get,offsy(j),elem,j+1000,etab,sy
$ w$ \  p7 G2 Q! [( M" ~, V      *get,offsz(j),elem,j+1000,etab,sz
' p7 f2 u" [' N, R3 h      *get,offsyz(j),elem,j+1000,etab,syz
   *enddo

*do,j,1,5000
*get,offstrainsita(j),elem,j+1000,etab,strainy
      *get,offstrainz(j),elem,j+1000,etab,strainz
      *get,offstrainzsita(j),elem,j+1000,etab,strainyz
   *enddo
% a8 ~# C  n( ?
   *do,j,1,10
      *do,k,(j-1)*500+1,j*500
% z+ w7 k$ ^! f3 a! q1 i7 i8 \        ons1(k)=m(j)*m(j)*offsz(k)+n(j)*n(j)*offsy(k)-2*m(j)*n(j)*offsyz(k)
3 O  d  v0 U8 \6 i" X        ons2(k)=n(j)*n(j)*offsz(k)+m(j)*m(j)*offsy(k)+2*m(j)*n(j)*offsyz(k)
        ons6(k)=m(j)*n(j)*offsz(k)-m(j)*n(j)*offsy(k)+(m(j)*m(j)-n(j)*n(j))*offsyz(k)
      *enddo
" [- L9 W9 r  S    *enddo
1 i" ?0 b% T4 t  Z
   *do,j,1,10
      *do,k,(j-1)*500+1,j*500
1 W0 s9 @5 x* t2 B6 d5 U: yonstrain1(k)=m(j)*m(j)*offstrainz(k)+n(j)*n(j)*offstrainsita(k)-m(j)*n(j)*offstrainzsita(k)
1 v6 e# j; @  ^4 T6 {onstrain2(k)=n(j)*n(j)*offstrainz(k)+m(j)*m(j)*offstrainsita(k)+m(j)*n(j)*offstrainzsita(k)
6 _+ w2 [, \5 q) ?9 `0 i2 N& Nonstrain6(k)=2*m(j)*n(j)*(offstrainz(k)-offstrainsita(k))+(m(j)*m(j)-n(j)*n(j))*offstrainzsita(k)
& [3 I" _3 M! y      *enddo
3 B* Q" ~* a( h4 z5 J( A# D" S   *enddo

*do,j,1,5000
5 L7 c) x% P4 G) a4 ?5 c       Y1(j)=abs(0.5*ons1(j)*onstrain1(j))
Y2(j)=abs(0.5*ons2(j)*onstrain2(j))
Y6(j)=abs(0.5*ons6(j)*onstrain6(j))
*enddo

   *do,j,1,5000
  @, b- t5 m$ d) ?2 H3 x* @      H1(j)=F11*ons1(j)*ons1(j)
) g5 u' {$ L- a. e% b1 r      H2(j)=F22*ons2(j)*ons2(j)
      H6(j)=F66*ons6(j)*ons6(j)
      H12(j)=2*F12*ons1(j)*ons2(j)
      Tsaiwu(j)=H1(j)+H2(j)+H6(j)+H12(j)
   *enddo
, [* q" u7 Y/ b& n+ q   parsav,all,parameter,txt

3 ]1 T: E" y, E3 T! y" k/ r* q  /solu
/ r% S1 f  g  v$ G6 U5 m    antype,,rest
parres,new,parameter,txt
. {* b" b2 J$ Q+ n1 D- {  *do,j,1,1000
      *if,ls(j),gt,324,then
         ekill,j
         LFN=LFN+1
      *endif
1 N2 {+ d. }1 F& \: w9 y   *enddo

    *do,j,1,5000
       *if,Tsaiwu(j),gt,1,then
           *if,H1(j),gt,H2(j),and,H1(j),gt,H6(j),then
                D1(j)=1-exp(-0.01*Y1(j))
FFN=FFN+1
9 O& }' v  v; y4 T            *elseif,H2(j),gt,H1(j),and,H2(j),gt,H6(j),then
                D2(j)=Y2(j)/40
( b. R/ _- V6 e' \& D                MCN=MCN+1
             *else
                D6(j)=Y6(j)/15
+ X% ?5 A. k0 b$ i5 K% s                SSN=SSN+1
9 x" c  s; y* w, C" W             *endif
8 V' A: Z! X- H* W3 c9 r        *endif
*enddo

MLFN(i-1)=LFN
MFFN(i-1)=FFN
MMCN(i-1)=MCN
MSSN(i-1)=SSN

/ \2 P" |; @/ I8 }  F, S/prep7
. _! ^! {. D5 V( J8 a# F. n7 x2 x: y     *do,j,1,5000
. I2 H% |% Z; ]2 s5 Y# U        TBDE,ANEL,j+1
. w. {7 l2 E( H, I# [& x; _& T) R     *enddo
*use,material

/solu
    sfa,4,1,pres,30+i*100/100
    nsubst,5
    cnvtol,f,,0.1,,0.1
ncnv,2
' S  G2 \( [1 O* N, m7 Y    solve
2 s; l$ @4 F9 M4 A( a6 ^6 B8 m*enddo
/ r" _! T* ~, L* ], Z6 S2 N) P: I
5 G# k1 V" x" N) P   ! 下面是运用弧长法计算下降段的宏文件程序
. i; C& }% @/ J. P/ d" L  b+ ^antype,,rest           ！从载荷步92进行重启动分析
     sfa,4,1,pres,140        ! 施加下降段载荷140MPa
9 m6 [7 F$ L, c  W1 f1 v     arclen,on,1000         ！弧长法弧长设置
     nsubst,500,1000,250    ！载荷子步设置
) A9 Q% E2 D- T0 s; [) ~  J$ q     neqit,15              ！一个载荷步内最多叠代15次，然后进行弧长二分
$ N6 @2 q* @" e" p# e: f! q! l0 [     cnvtol,f,,0.1,,0.1       ！定义力收敛准则
ncnv,2
     solve

5 [+ f  J- ^+ b  E! 创建宏文件，提取随内压上升的纤维/基体/剪切失效次数
2 N! r. Z4 W! g$ k' D' w( S*create,data1
% U0 a: b) T- I! @/ H) L  *cfopen,d:\file1,txt
/ M* j. D  u% [6 T  \, J$ M       *vwrite, MLFN(1)
1 m2 N* f; S9 L9 A( p" N3 R        (f5.0)
*cfclose
/ v2 t) N3 x* k$ O6 m; z*end

*create,data2
; K$ m& ?) a/ |; t  ]) w  *cfopen,d:\file2,txt
/ c9 ~. [3 y: ~" {- F  M* D      *vwrite, MFFN(1)
       (f5.0)
8 ~0 m* [% }6 `' I: b$ w8 w/ O. z*cfclose
*end

*create,data3
  *cfopen,d:\file3,txt,
     *vwrite, MMCN(1)
      (f5.0)
2 \" |2 P/ i& r& a: F" f4 l: ?  ~*cfclose
*end

*create,data4
  *cfopen,d:\file4,txt
$ V) T' f! O* d$ [; X      *vwrite, MSSN(1)
       (f5.0)
     *enddo
*cfclose
6 M* S- Y# _- P5 P* R7 d/ i- |% j*end
6 s: \/ }. Y* s
*use,data1
*use,data2
4 b0 l; D" _( g1 g" J% z8 p*use,data3
! H/ }+ z0 @% I" ^*use,data4
* f- m4 E& j( o/ F. R
/title,

lbt18872

比较复杂，比较专业，顶一下

acailaogu

不错 啊，，

ling_0212

很不错，谢谢分享。

amore06

运行不了哈哈！~

wangbo119

不错，顶一下！！！

lm8403

楼主  请问哪个模块是做铺放的呢？

tommot222

有点看不懂呢  有点伸

ghij067

楼主我们一起帅吧！！喜欢的就来
神鬼传奇外挂和
1 M7 y" j9 x& l1 b. }  n9 e2 l/ `DNF外挂和
AION外挂和
剑侠情缘3外挂和
传奇世界外挂和

yunyun5258

好东西，感谢楼主分享！

smallmj

很好很强大

laoda01

哈哈！牛！太好了！谢谢！

salutsun

9# ghij067
9 楼的。。一般这的人不玩网游。。

sp60

宏观力学的？

sxyuyifan

楼主真强，是浙大化机所的啊。

shuyu2025

嗯 的确很强大

octor

恩 真厉害真厉害 羡慕~~

ruohan04741

这篇文章非常好，谢谢楼主

ruohan04741

请问下楼主，我想用SOLID164单元建立一碳纤维管材模型，如何用单元坐标系来表示纤维的铺层方向？有没有具体的命令，您的这篇文章非常不错。我在文献中也看到用局部柱坐标系来表示，但具体操作起来还是不太理解，谢谢

hamily

，你的问题搞明白没，可以QQ409957906 进行交流，Email：
; A( H5 a, x. cwmmpx1985@vip.qq.com
19# ruohan04741 你好