平板开小圆孔应力集中的三维有限元分析-北京化工大学.PDF

第 30 卷 第 3 期 北 京 化 工 大 学 学 报 Vol. 30 , No. 3 2003 年 JOURNAL OF BEIJ IN G UN IV ERSIT Y OF CHEMICAL TECHNOLO GY 2003 平板开小圆孔应力集中的三维有限元分析 黄炳臣 　冷纪桐 (北京化工大学机电工程学院 ,北京 　100029) 摘 　要 : 针对经典的弹性力学中平板开小圆孔问题 ,采用有限元方法 ,把平板开孔问题视为三维应力状态 ,考虑在 Z 向的应力变化 ,给出了平板开孔问题各应力在孔边应力集中区域沿厚度的变化情况 ;并且探讨了所开圆孔的大 小对应力集中的影响。 关键词 : 平板开小圆孔 ;应力集中 ;有限元分析 中图分类号: O343. 2 施加在右侧面 C 上 ,大小为 005 MPa 。 引　言 网格采用规则的六面体单元划分。由于在孔的 平板开小圆孔的应力集中问题是实用工程领域 周围会出现应力集中 ,而且这里也是要重点研究的 中一个较为常见的问题 ,也是弹性力学平面问题的 地方 ,故在其周围设置单元的总体尺寸为 0005 m ; 一个经典问题。在平面问题中 , 只考虑三个应力分 在沿开孔圆的边缘圆弧上划分成 14 份 ,其余部分的 量 ,忽略其在厚度方向上的变化 ,推导出理论公式。 单元总体尺寸设置为 0008 m 。由于需要研究在板 从工程实用角度上说 ,这不失为一种化繁为简的好 厚方向上的应力分布情况 ,故在板的厚度方向上把 方法 ,因此得到了广泛的应用 ,但是这种提法从三维 板分成 5 层单元。图 1 显示的是平板开小圆孔的有 理论来看是近似的。本文应用三维有限元模型 ,分 限元模型。 析了开小圆孔的方平板受单向拉伸时孔边的应力情 况并与二维理论的计算结果对比。 1 　板开孔问题的数值分析 11 　模型的描述 本算例中所采用的材料假设其杨氏弹性模量为 5 E = 2 ×10 MPa ,泊松比为 μ= 0 3 。模型的基本尺 寸为 :板厚 d 为 0015 m ,板的大小为 04 m ×04 m , 圆孔开在板的正中央 ,半径 R 取五种尺寸 ,分别 为 0012 、001 、0007 、0004 、0003 m 。 应用 ANSYS 计算程序 ,三维分析时采用 SOL 图 1 　平板开圆孔的有限元模型 ID186 号三维体单元 ,它由 20 个节点定义而成。在 Fig. 1 　FEA model of the plat with a small circle hole 12 　分析数值模拟的结果 进行二维分析时采用 4 节点的 PLAN E42 号平面单 对于 σ ,在考虑厚度方向应力的情况下 ,它沿 元。 x 厚度方向的分布并不是均匀的。其最大值出现在 为了减少计算量 ,考虑到所分析的实体是一个 A 面与圆孔边的交接处 ,位于厚度的中间位置上。 对称结构 ,故可以采用四分之一结构。在左侧面 A 最小值出现