计算三维问题的应力强度因子的程序（A procedure for calculating the stress intensity factor of a three-dimensional problem）.docVIP

下载本文档

0
0
约5.04千字
约 14页
2017-08-18 发布于河南
举报
版权申诉

计算三维问题的应力强度因子的程序（A procedure for calculating the stress intensity factor of a three-dimensional problem）.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

计算三维问题的应力强度因子的程序（A procedure for calculating the stress intensity factor of a three-dimensional problem）

计算三维问题的应力强度因子的程序（A procedure for calculating the stress intensity factor of a three-dimensional problem） com，ANSYS媒体相关。60（090601）参考文献验证。手册：相对。六十 /验证，vm143 *创建，分形，MAC / NOPR 汉塞尔，所有 * get，n，节点，努姆，max！当前最大节点数 cmsel，S，裂纹尖端！选择顶端节点 ESLN！任何元素 *，玛克斯，元，，努姆，麦克斯！当前最大元素数量 *做，IEL，1，玛克斯！关于MAX元素的循环埃尔米= IEL *如果ELMI，乐，0，退出！没有更多的选择 *，eltype，元素，埃尔米、属性、类型！得到的元素类型 *如果eltype，NE，arg1，周期！选中元素的检查 N3 = nelem（ELMI，3）！获取节点3（k） *如果，汉塞尔（N3），乐，0，周期！它必须被选中。 N7 = nelem（ELMI，7）！获取节点7（l） *如果，汉塞尔（N7），乐，0，周期！它也必须被选中。 N1 = nelem（ELMI，1）！获取节点1（i） N2 = nelem（ELMI，2）！获取节点2（j） N5 = nelem（ELMI，5）！获取节点5（m） N6 = nelem（ELMI，6）！获取节点6（n） X3 = 0.75×NX（N3）！N3加权位置 Y3 = 0.75×纽约（N3） Z3 = 0.75 * NZ（N3） x = 0.25×NX（N2）+ X3！四点位置（节点（R）） Y = 0.25×纽约（N2）+ Y3 Z = 0.25 * NZ（N2）+ Z3 n＝n＋1！下一个节点 N10 = N N，X，Y，Z N10！中间节点的位置 x = 0.25×NX（N1）+ X3 Y = 0.25×纽约（N1）+ Y3 Z = 0.25 * NZ（N1）+ Z3 n＝n＋1 N12 = N N，X，Y，Z N12， X7＝0.75×NX（N7） Y7 = 0.75×纽约（N7） Z7 = 0.75 * NZ（N7） x = 0.25×NX（N6）+ X7 Y = 0.25×纽约（N6）+ Y7 Z = 0.25 * NZ（N6）+ Z7 n＝n＋1 N14 = N N，X，Y，Z N14， x = 0.25×NX（N5）+ X7 Y = 0.25×纽约（N5）+ Y7 Z = 0.25 * NZ（N5）+ Z7 n＝n＋1 N16 = N N，X，Y，Z 16 N4 = n3 N8 = N7 汉塞尔，所有类型3 恩，N1、N2、N3埃尔米，N4、N5、N6、N7、N8！重新定义的元素更多的，0，n10,0，n12,0，n14,0，N16 更多的， * enddo cmsel，U，裂纹尖端！选择端节点 nummrg，结！合并中节点汉塞尔，所有！选择所有元素场，所有！选择所有元素 / gopr *结束 / prep7 *有趣度 innerradius = 0.1！innerradius为裂纹半径 outerradius = 1！outerradius为圆柱半径定标器= 0.025！为裂纹前沿单元范围定标器， Already the best BaseHeight=0.51! BaseHeight for base height, it, may, be, the, best LayerHeight=0.18! LayerHeight is the height of the expansion layer. The bigger the difference, the better LayerAmount=16! LayerAmount is the number of layers, but it has little to do with accuracy RotationAngle=6! RotationAngle is the angle of rotation of the unit Rotationtimes=90/RotationAngle! Rotationtimes is the number of rotations Height=LayerAmount*LayerHeight+BaseHeight! Height is the total height SMRT, OFF /TITLE, VM143, FRACTURE, MECHANICS, STRESS, INTENSITY - CRACK,, IN, A, FINITE, WIDTH, PLATE C***, BROWN, AND, SRAWLEY, ASTM, S