含边端半圆孔有限宽平板的理论应力集中系数与孔径板宽的变化关系讲解

1、有限元作业报告 含边端半圆孔有限宽平板的理论应力集 中系数与孔径 / 板宽的变化关系 2013年 11月 目录 目录 I 第 1 章 目的和意义 1 1.1 应力集中 1 1.2 理论应力集中系数 1 1.3 平面问题的应力集中 1 1.4 利用有限元分析 1 第 2 章 软件介绍 2 2.1 MSC.PATRAN 2 2.2 MSC.NASTRAN 2 第 3 章 基本模型和数据 4 3.1 基本模型 4 3.2 基本数据 4 第 4 章 计算分析 5 学习体会和建议 11 第 1 章 目的和意义 第 1 章 目的和意义 1.1 应力集中 应力集中是应力在固体局部区域内显著增高的现象。多出现

2、于尖角、孔洞、缺 口、沟槽以及有刚性约束处及其邻域。应力集中会引起脆性材料断裂；使物体产生 疲劳裂纹。在应力集中区域，应力的最大值（峰值应力）与物体的几何形状和加载 方式等因素有关。局部增高的应力值随与峰值应力点的间距的增加而迅速衰减。由 于峰值应力往往超过屈服极限而造成应力的重新分配，所以，实际的峰值应力常低 于按弹性力学计算出的理论峰值应力。 1.2 理论应力集中系数 工程上，用应力集中系数来表示应力增高的程度。它是峰值应力max 与不考虑 应力集中时适当选取的基准应力 n的比值，恒大于 1 且与载荷大小无关。 峰值应力 可根据弹性力学的理论或有限元计算得到，也可根据光弹性实验或其他实验应

3、力分 析方法测定。其中基准应力是人为规定的，取值方式不唯一，在后面的分析中给出 两种常用表达方式。本文将用 来表示理论应力集中系数，其表达式为： max n 1.3 平面问题的应力集中 由弹性理论可知，对于平面问题的应力集中问题，应力分布只取决于构建的形 状以及边界上的载荷条件，而与材料的弹性常数无关，于是在实验应力分析中，利 用力学性质不同的材料制作模型，以替代真实构件进行应力分析，对平面问题是精 确的。 1.4 利用有限元分析 早起理论应力集中系数通过复杂的计算，可以得到精确解。随着实验手段的发 展，如电测法、光弹性法、散斑干涉法、云纹法等实验手段（见实验应力分析）均 可测出物体的应力集中

4、。随着科技的进步，计算机和有限元法以及边界元法的迅速 发展，为寻找应力集中的数值解开辟了新途径。本文利用PATRAN进行有限元方法 的数值计算便是一种快速有效的方法。 第 2 章 软件介绍 第 2 章 软件介绍 2.1 MSC.PATRAN 图 2-1 MSC.PATRAN 软件 MSC.PATRA最N早由美国宇航局 (NASA)倡导开发的 , 是工业领域最著名的并行 框架式有限元前后处理及分析系统， 其开放式、多功能的体系结构可将工程设计、 工程分析、结果评估、用户化身和交互图形界面集于一身，构成一个完整CAE集成 环境。 MSC.PATRA提N 供了功能全面 , 方便灵活的可满足各种分析精

5、度要求的复杂有 限元的建模能力。 其综合全面、先进的网格划分技术 ,为用户根据不同的几何模型 提供了多种不同的生成和定义有限元模型工具 , 包括: 多种网格划分器、有限元模 型的编辑处理、单元设定、任意梁截面建模、边界和载荷定义、及交互式计算结果 后处理。 MSC.PATRA提N 供了多种计算分析结果可视化工具 ,帮助工程师灵活、快速地理 解结构在载荷作用下复杂的行为 , 如结构受力、 变形、温度场、 疲劳寿命、 流 体流动等。 分析的结果同时可与其它有限元程序联合使用。 2.2 MSC.NASTRAN 作 为 世 界 CAE 工 业 标 准 及 最 流 行 的 大 型 通 用 结 构 有 限

6、 元 分 析 软 件 , 第 2 章 软件介绍 MSC.NASTRA的N 分析 功能覆盖了绝大多数工程应用领域 ,并为用户提供了方便的模 块化功能选项 ,MSC.NASTRAN的 主要功能模块有 : 基本分析模块 （ 含静力、 模态、 屈曲、热应力、流固耦合及数据库管理等 ） 。 动力学分析模块、 热传导模块、 非 线性分析模块、 设计灵敏度分析及优化模块、 超单元分析 模块、 气动弹性分析 模块、 DMAP用户开发工具模块及高级对称分析模块。 除模块化外 , MSC.NASTRAN 还按解题规模分成 10,000 节点到无限节点 ,用户引进时可根据自身的经费状况和功 能需求灵活地选择不同的模

7、块和不同的解题规模 , 以最小的经济投入取得最大效 益。MSC.NASTRA及N MSC的相关产品拥有统一的数据库管理 ,一旦用户需要可方便地 进行模块或解题规模扩充 , 不必有任何其它的担心。 MSC.NASTRA以N每年一个小版本 , 每两年一个大版本的速度更新 , 用户可不断 获得当今 CAE 发展的最新技术用于其产品设计。 目前 MSC.NASTRA的N最新版本是 1999 年发布的 V70.5 版。新版本中无论在设计优化、 P 单元、 热传导、 非线性 还是在数值算法、 性能、 文档手册等方面均有大幅度的改进或突出的新增功能。 第 3 章 基本模型和数据 第 3 章 基本模型和数据

8、3.1 基本模型 本报告建立的模型为含边端半圆孔有限宽平板，并且是两边都有的对称情况 其形状如下： 图3-1 含边端半圆孔有限宽平板模型 3.2 基本数据 基本数据如下： 平板的几何特性：长度为 200mm，宽度为 100mm，厚度为 3mm 平板的材料特性：选用钢，定义弹性模量 E为 200 GPa，泊松比为 0.3 ， 半圆孔的半径分别选： 1mm,3mm,5m，m10mm， 20mm， 30mm， 40mm， 48mm 在两边施加等大的载荷 10N/mm，平板被拉伸。 建模的过程中，会随着孔径的变化，在孔附近撒种的数量进行变化，目的是为 了达到一定的精确度，提高结果的准确性。 第 4 章

9、 计算分析 第 4 章 计算分析 从图中可以看出我们要分析的模型具有对称性， 这里我们将利用模型的 1/4 进 行分析，从而可以得出整体的结果。事实上如果不采用这种方法，在进行分析时得 不到结果。利用对称性，可以得到两个位移约束条件。 如图，对于左面来说，没有 X方向的位移 ; 对于上面来说，没有 Y方向的位移。 在施加约束约束时也要固定 Z 方向的位移（否则可以发生 Z方向的平动，得不出结 果来）。 有趣的是这个 1/4 模型与研究中心圆孔有限宽平板的模型是一样的，不过由于 位移约束条件的差异，结果自然是有区别的。 图 4-1 1/4 模型的示意图 建模的过程首先是几何作图，这个模型的建立有

10、多种方法。先画出矩形平面， 在上面挖去 1/4 的圆，即扇形便可得到模型。也可先画出整体模型，再用 Break 的 方法逐渐得到 1/4 的部分。也可建立体模型，这样就不用在后面定义厚度了。 下面划分网格，网格的值设为 2，在圆弧周围撒种可以根据不同情况确定不同 的数目，这样可提高计算的精度。 施加位移、应力约束条件，定义材料，赋予属性。模型建立就完成了。利用相 似的方法完成其他几种情况的建模，计算并分析结果。 得到的图如下： 第 4 章 计算分析 图 4-2 R=1mm 图 4-3 R=3mm 图 4-4 R=5mm 6 第 4 章 计算分析 图 4-5 R=10mm 图 4-6 R=20m

11、m 图 4-7 R=30mm 第 4 章 计算分析 图 4-8 R=40mm 图 4-9 R=48mm 上图是对整个左面的应力分布图，可以看出在靠近应力集中地地方，应力变化 8 第 4 章 计算分析 很剧烈，远离应力集中的地方，应力变化不大则比较平滑。应力集中的地方是左面 上圆孔的边界处。 下面计算理论应力集中系数，给出两种表达式： max max 其中 其中 d 是板厚 FB B 2Rd 面列出表： 圆孔半 径 R/mm 1 3 5 10 20 30 40 48 max 5.92 8.23 9.28 9.68 10.3 11.8 19 84.1 2R/B 0.02 0.06 0.1 0.2 0.4 0.6 0.8 0.96 n 3.33 3.33 3.33 3.33 3.33 3.33 3.33 3.33 n 3.40 3.55 3.70 4.17 5.55 8.33 16.67 83.33 1.78 2.47 2.78 2.90 3.09 3.54 5.7 25.23 1.74 2.32 2.51 2.32 1.85 1.42 1.14 1.01 根据表格作图可以得到如下曲线： 第 4 章 计算分析 在两种表达方式中，其中一种在半径接近半宽时，应力集中系数迅速增大，半 径小时，变化缓和。另外一种随半径增大，应力集中系数先增大后减小。可以发现 应力集中系数始终大于