孔边应力集中和裂纹尖端应力强度因子的有限元分析.docVIP

第!第#期卷 !$$!年!月 交通科学 武汉理工大学学报%与工程版 ’()*+,-(./)0,+1+234*5267(.8490+(-(:78*,+5;(*6,62(+924+94=+:2+44*2+:% ?(-@!A(@# B4C@!$$! 孔边应力集中和裂纹尖端应力 强度因子的有限元分析D 尹奇志#肖金生#吕运冰!李格升# 武汉理工大学能源与动力工程学院#武汉% 武汉理工大学交通学院!武汉FG$$G% FG$$G 摘要H应用IA圆孔板孔边裂纹进行了有限元分J软件对圆孔板的应力集中以及平板中心裂纹K 析@计算了圆孔平板在单向或双向受力下的应力分布以及中心裂纹K孔边裂纹的应力强度因子L并与理论解进行了比较L其IAJ解与解析解均比较接近@关键词H裂纹M有限元M断裂力学M应力分布M应力强度因子中图法分类号H#!O8N $引 言 #圆孔板的孔边应力集中分析 S@S单向受力圆孔平板的应力分布 所示L正方形平板中S@S@S计算模型如图#,心有一圆孔L圆孔半径TU!VVL正方形边长为 O 厚度为OVVL弹性模量为!#$$VVLW#$LXY,泊松比为$在正方形左右两边加一Z方向均@G@ 布拉力L其集度为[U#$$XY@,应用IA过程如下H将组J对其进行计算L成圆孔的四条线段每条分成四边形每条$等份L边分成F再采用六节点平面三角形单元$等份L 网格划分好后L%进行自由网格划分@Y\IA!对正方形左边中间节点的]方向以及左边所有节点Z方向进行约束@对正方形右边所有节点L运 用IA选项进行加载@由于_J中的^Y*455)*4物体表面的内法线方向为正L故载荷密度值为‘其网格划分K约束与加载见图#加@#$$XY@,C载完成后L对该模型进行计算@S@S@a结果及分析利用IAJ的通用后处理器Y可得到圆孔b8#对计算结果进行分析L平板Z方向的应力分布与]方向的应力分布图L如图!所示@ 由图!可见LcZ在圆孔的上端和下端有最大 在机械工程中L有的工件在设计过程中其结另外L金属构件由于环境影响在构自身含有孔洞L 制造K运输或使用过程中也会产生微小的缺陷和裂纹@这些孔洞K缺陷和裂纹在外界载荷作用下通常会引起其周围区域%特别是裂纹尖端的应力集 #Q 中L进而使工件或容器失效P@ 对于孔洞以及裂纹的应力集中问题L目前分析的方法较多@主要分为实验法和数值计算法@实 !Q验主要采用光测法P但其费用高K较复杂L而且L 难以达到很高精确度@数值法目前主要运用权函 GRFQ数法K边界元法K有限元法P等L其中有限元法 应用较为广泛@ IAJ软件是一种有限元分析计算软件L它能较好地计算孔洞以及裂纹周围区域的应力分布L并能计算裂纹的应力强度因子K积分以及能’量释放率等L其特点是简单K经济L而且精度高@应用IAJ软件对圆孔板的应力集中以及中心裂纹K圆孔孔边裂纹板进行了有限元分析@计算了圆孔平板在单向和双向受力下的应力分布以及中心裂纹K圆孔孔边裂纹的应力强度因子L并与理论解进行了比较@ E收稿日期H!$$#$! 尹奇志H男L硕士生L主要研究领域为液化气容器安全问题!O岁L D湖北省自然科学基金资助项目 批准号H及交通部重点项目%批准号H资助%!$$$#!ddOe$e$!e!$’ 卷 由于应用QRS为@@TS得到的是%%’和N’了便于QRS将极坐标公TS解与解析解的比较@式D转变为直角坐标公式@应力分量由极坐标向$ -H直角坐标转换的公式为G ::)%J%7:%KLMFMOPFNMOPFKLMFEFEF ::)%J%7:%MOPFKLMFNKLMFMOPF’EFEF ::D:)D7%JND7MN%MOPFKLMFKLMFOPF’EFEF 选取轴正半轴和’轴正半轴上的节点作 :::G$7:HJKG$7:HG$:HLMF::EEE:3$J:H)G7K:G$3H%LMFF::EE ::)7MG$7:HG$:HNOPFEF:EE为比较点@应用QRSTS求得这些节点方向和再运用圆孔平板应力的理论公式’方向的应力+D$和极坐标向直角坐标的转换公式D求得这些D$:节点方向和’方向的应力+选取轴正半轴上所有节点@ 并绘制这些节 第,期尹奇志等_孔边应力集中和裂纹尖端应力强度因子的有限元分析^DQ^ 点!方向和方向的应力曲线#见图$由图$可% 见#对于单向受力的圆孔平板#其应力分布的’()(解与解析解能很好地吻合%*%+双向受力圆孔平板的应力分布 模型#在其左右和上下同时施加等对图,.-值拉力#其集度为/同样#当40,1123#56时#-取6其上任意一点89的应力分布的解701##.4:; = 析表达式为?:0 @@@@ ,A67:#?,B67:#C;0:;01@@ 该数值解与其解析解也能很好布的’()(解# 地吻合% @平板中心裂纹尖端应力分析 +%*数学模型 如图D所示#边