ansys复合材料渐进失效

/title,pprogressivefailureanalysisofcarbonfiber/epoxycompositelaminatesusingCDM(40MPahydrogenstoragetank)/uis,msgpop,3/prep7et,1,95et,2,64mp,ex,1,70e3mp,nuxy,1,0.3tb,biso,1,1,2tbdata,,246,600*afun,degE1=181e3E2=10.3e3v12=0.28v23=0.49*dim,ls,array,1000*dim,offsy,array,5000*dim,offsz,array,5000*dim,offsyz,array,5000*dim,offstrainz,array,5000*dim,offstrainsita,array,5000*dim,offstrainzsita,array,5000*dim,ons1,array,5000*dim,ons2,array,5000*dim,ons6,array,5000*dim,onstrain1,array,5000*dim,onstrain2,array,5000*dim,onstrain6,array,5000*dim,H1,array,5000*dim,H2,array,5000*dim,H6,array,5000*dim,H12,array,5000*dim,Tsaiwu,array,5000*dim,Y1,array,5000*dim,Y2,array,5000*dim,Y6,array,5000F22=1.0/(298*298)F66=1.0/(778*778)F11=1.0/(2150*2150)F12=-0.5*sqrt(F11*F22)LFN=0FFN=0MCN=0SSN=0*dim,MLFN,array,100*dim,MFFN,array,100*dim,MMCN,array,100*dim,MSSN,array,100*dim,D1,array,5000*dim,D2,array,5000*dim,D6,array,5000*dim,angle,array,10*dim,m,array,10*dim,n,array,10*dim,te,array,6,6,5000*dim,ts,array,6,6,5000*dim,tst,array,6,6,5000*dim,unit,array,6,6,5000*dim,s,array,6,6,5000*dim,c,array,6,6,5000*dim,cp,array,6,6,5000*dim,midlpf,array,6,6,5000angle(1)=90angle(2)=-90angle(3)=18.9angle(4)=-18.9angle(5)=90angle(6)=-90angle(7)=28.6angle(8)=-28.6angle(9)=90angle(10)=-90*do,i,1,10*vfun,m(i),cos,angle(i)*vfun,n(i),sin,angle(i)*enddo*do,i,1,10*do,j,(i-1)*500+1,i*500te(1,1,j)=m(i)*m(i),n(i)*n(i),0,0,0,-2*m(i)*n(i)te(1,2,j)=n(i)*n(i),m(i)*m(i),0,0,0,2*m(i)*n(i)te(1,3,j)=0,0,1,0,0,0te(1,4,j)=0,0,0,m(i),n(i),0te(1,5,j)=0,0,0,-n(i),m(i),0te(1,6,j)=m(i)*n(i),-m(i)*n(i),0,0,0,m(i)*m(i)-n(i)*n(i)*enddo*enddo*do,i,1,5000unit(1,1,i)=1,0,0,0,0,0unit(1,2,i)=0,1,0,0,0,0unit(1,3,i)=0,0,1,0,0,0unit(1,4,i)=0,0,0,1,0,0unit(1,5,i)=0,0,0,0,1,0unit(1,6,i)=0,0,0,0,0,1*enddo*do,i,1,5000*moper,tst(1,1,i),te(1,1,i),solv,unit(1,1,i)*enddo*do,i,1,5000*mfun,ts(1,1,i),tran,tst(1,1,i)*enddo*do,i,1,5000s(1,1,i)=1.0/E1,-v12/E1,-v12/E1,0,0,0s(1,2,i)=-v12/E1,1.0/E2,-v23/E1,0,0,0s(1,3,i)=-v12/E1,-v23/E1,1.0/E2,0,0,0s(1,4,i)=0,0,0,2*(1+v23)/E2,0,0s(1,5,i)=0,0,0,0,2*(1+v12)/E1,0s(1,6,i)=0,0,0,0,0,2*(1+v12)/E1*enddo*do,i,1,5000*moper,c(1,1,i),s(1,1,i),solv,unit(1,1,i)*enddo*create,material*do,j,1,5000c(1,1,j)=(1-D1(j))*c(1,1,j)c(2,2,j)=(1-D2(j))*c(2,2,j)c(6,6,j)=(1-D6(j))*c(6,6,j)*enddo*do,j,1,5000*moper,midlpf(1,1,j),ts(1,1,j),mult,c(1,1,j)*enddo*do,j,1,5000*moper,cp(1,1,j),midlpf(1,1,j),mult,tst(1,1,j)*enddo*do,j,1,5000TB,ANEL,j+1,1,21,0TBDATA,,cp(1,1,j),cp(1,2,j),cp(1,3,j),0,0,cp(1,6,j)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),0TBDATA,,cp(5,5,j),0,cp(6,6,j)*enddo*end*use,materialcsys,1cyl4,0,0,44,0,45.8,30,160*do,i,1,10cyl4,0,0,45.8+(i-1)*0.42,0,45.8+i*0.42,30,160*enddovglue,allnumcmp,allLESIZE,12,,,50LESIZE,1,,,2LESIZE,3,,,2LESIZE,6,,,2LESIZE,8,,,2LESIZE,2,,,10LESIZE,4,,,10LESIZE,5,,,10LESIZE,7,,,10*do,i,1,9LESIZE,13+(i-1)*4,,,10*enddo*do,i,1,9LESIZE,14+(i-1)*4,,,10*enddo*do,i,1,9LESIZE,15+(i-1)*4,,,50*enddo*do,i,1,9LESIZE,16+(i-1)*4,,,50*enddo*do,i,1,10LESIZE,55+(i-1)*4,,,1*enddo*do,i,1,10LESIZE,56+(i-1)*4,,,1*enddoallstype,1mat,1vmesh,1type,2*do,i,1,10mat,i+1vmesh,i+1*enddonummrg,node,,,,lownummrg,elem,,,,low*do,i,1,5000emodif,i+1000,mat,i+1*enddo/pnum,mat,1eplot/soluda,5,symmda,6,symm*do,i,1,10da,(i-1)*4+19,symm*enddo*do,i,1,10da,(i-1)*4+20,symm*enddod,1000,uz,0rescontrol,define,all,1,1sfa,4,1,pres,30nsubst,5cnvtol,f,,0.1,,0.1solve*do,i,2,92/post1rsys,1etable,sequ,s,eqvetable,sy,s,yetable,sz,s,zetable,syz,s,yzetable,strainy,epel,yetable,strainz,epel,zetable,strainyz,epel,yz*do,j,1,1000*get,ls(j),elem,j,etab,sequ*enddo*do,j,1,5000*get,offsy(j),elem,j+1000,etab,sy*get,offsz(j),elem,j+1000,etab,sz*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*do,j,1,10*do,k,(j-1)*500+1,j*500ons1(k)=m(j)*m(j)*offsz(k)+n(j)*n(j)*offsy(k)-2*m(j)*n(j)*offsyz(k)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*enddo*do,j,1,10*do,k,(j-1)*500+1,j*500onstrain1(k)=m(j)*m(j)*offstrainz(k)+n(j)*n(j)*offstrainsita(k)-m(j)*n(j)*offstrainzsita(k)onstrain2(k)=n(j)*n(j)*offstrainz(k)+m(j)*m(j)*offstrainsita(k)+m(j)*n(j)*offstrainzsita(k)onstrain6(k)=2*m(j)*n(j)*(offstrainz(k)-offstrainsita(k))+(m(j)*m(j)-n(j)*n(j))*offstrainzsita(k)*enddo*enddo*do,j,1,5000Y1(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,5000H1(j)=F11*ons1(j)*ons1(j)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)*enddoparsav,all,parameter,txt/soluantype,,restparres,new,parameter,txt*do,j,1,1000*if,ls(j),gt,324,thenekill,jLFN=LFN+1*endif*enddo*do,j,1,5000*if,Tsaiwu(j),gt,1,then*if,H1(j),gt,H2(j),and,H1(j),gt,H6(j),thenD1(j)=1-exp(-0.01*Y1(j))FFN=FFN+1*elseif,H2(j),gt,H1(j),and,H2(j),gt,H6(j),thenD2(j)=Y2(j)/