原创力文档- 1、本文档共8页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
三维裂纹的命令流
三维裂纹的命令流 *create,fract,mac /nopr nsel,all *get,n,node,,num,max cmsel,s,cracktip esln *get,elmax,elem,,num,max *do,iel,1,elmax elmi=iel *if,elmi,le,0,exit *get,eltype,elem,elmi,attr,type *if,eltype,ne,arg1,cycle n3=nelem(elmi,3) *if,nsel(n3),le,0,cycle n7=nelem(elmi,7) *if,nsel(n7),le,0,cycle n1=nelem(elmi,1) n2=nelem(elmi,2) n5=nelem(elmi,5) n6=nelem(elmi,6) x3=0.75*nx(n3) y3=0.75*ny(n3) z3=0.75*nz(n3) x=0.25*nx(n2)+x3 y=0.25*ny(n2)+y3 z=0.25*nz(n2)+z3 n=n+1 n10=n n,n10,x,y,z x=0.25*nx(n1)+x3 y=0.25*ny(n1)+y3 z=0.25*nz(n1)+z3 n=n+1 n12=n n,n12,x,y,z x7=0.75*nx(n7) y7=0.75*ny(n7) z7=0.75*nz(n7) x=0.25*nx(n6)+x7 y=0.25*ny(n6)+y7 z=0.25*nz(n6)+z7 n=n+1 n14=n n,n14,x,y,z x=0.25*nx(n5)+x7 y=0.25*ny(n5)+y7 z=0.25*nz(n5)+z7 n=n+1 n16=n n,n16,x,y,z n4=n3 n8=n7 nsel,all type,3 en,elmi,n1,n2,n3,n4,n5,n6,n7,n8 emore,0,n10,0,n12,0,n14,0,n16 emore, *enddo cmsel,u,cracktip nummrg,node nsel,all esel,all /gopr *end /prep7 *afun,deg innerradius=0.5 outerradius=1 scaler=0.015 baseheight=0.51 layerheight=0.18 layeramount=20 rotationangle=6 rotationtimes=90/rotationangle height=layeramount*layerheight+baseheight smrt,off /title,fracture mechanics stress intensity - crack in a finite width plate /com,? ?? ?********crack in 3-Ddimensions using solid45 and solid95 antype,static et,1,solid45 et,2,solid45 et,3,solid95 mp,ex,1,2e4 mp,nuxy,1,0.3 local,33,1,innerradius csys,33 n,1 ngen,9,20,1 n,11,scaler n,171,scaler,180 fill,11,171,7,31,20 local,33,0,innerradius csys,33 fill,1,11,9,2,1,9,20,3 n,15,outerradius-innerradius n,75,outerradius-innerradius,baseheight fill,15,75,2,35,20 n,155,-innerradius,baseheight ngen,2,2000,155 csys,0 n,135,0.6,baseheight csys,33 fill,75,135,2,95,20,,,0.5 n,172,-innerradius ngen,2,2000,172 fill,155,172,5,177,-1,,,0.15 ngen,2,2000,173,177,1 csys,33 fill,11,15,3,,,7,20,3 fill,171,172,5,180,7,,,2 fill,135,155,5,186,7 fill,180,186,5,181,1,5,7,6 local,44,1,,,,,-90,,0.75 csys,44 ngen,2,2000,all,,,,rotationangle csys,0 e,2,22,1,1,2002,2022,2001,2001 egen,8,20,-1 e,2,3,
文档评论（0）