含圆孔双边裂纹板的有限元计算.docVIP

级 专 业 学 号 学生姓名 指导教师 完成日期 通过有限元的分析方法,针对,分析了在荷载作用下的应力。 有限元 应力有限元分析利用数学近似的方法对真实物理系统（几何和载荷工况）进行模拟。利用简单而又相互作用的元素，即单元，就可以用有限数量的未知量去逼近无限未知量的真实系统。有限元法最初被称为矩阵近似方法，应用于航空器的结构强度计算，并由于其方便性、实用性和有效性而引起从事力学研究的科学家的浓厚兴趣。经过短短数十年的努力，随着计算机技术的快速发展和普及，有限元方法迅速从结构工程强度分析计算扩展到几乎所有的科学技术领域，成为一种丰富多彩、应用广泛并且实用高效的数值分析方法。 图1. 二、公式推导 含圆孔不对称双边裂纹板受拉伸荷载的解析解 对于较复杂的带裂纹的孔口问题, 下面利用复变函数方法讨论带不对称双边裂纹的圆孔口问题的解析解。 假设在无限大平面中含有一个带裂纹的圆孔洞, 如图1, 半径为R， 左边裂纹长为c-R，右边裂纹长为d-R。以圆孔中心为原点, 裂纹所在直线为x 轴建立直角坐标系。 如图2、3所示： 图2 图3 作保角映射如下： （1） 其中 映射后如图3所示 （d）→1, （Ri）→D, （-c）→-1, （-Ri）→D1, 则该映射将物理平面上带裂纹形孔洞的无限大平面保角映射到数学平面上的单位圆内部。 显然，只有一个奇点=，经计算得 （2） 其中Res表示留数. 又（1）式对求导得 （3） （1）式的共轭为 （4） 当孔洞边缘及其所带裂纹面上不受力,仅在受到沿y 轴方向的单向拉应力q令时 X=Y=0, B=, . （5） （6） （7） 经计算可知： （8） 考虑到 （9） （10） （11） 应用柯西积分，（7）式左边 （12） 将公式（6）代入（7）式得 （13） 因为=是和的奇点，根据留数理论得 （14） （15） 将（14），（15）和（2）代入（13）得 （16） 将（16）代入（5）有 （17） 又（17）对式对求导得 （18） 则， （19） ，在孔边=0，所以孔边应力 （20） 公式中Re为取其实部。 应力强度因子 （21） 三、数值计算结果及分析 1、在MATLAB中建立问题模型，如图4所示： 图4 2、孔边应力分布图 如图5所示： 图5 表1：精确孔边应力分布表 -坐标系中 A1 D1 C1 C D A 1.2422 1.6715 2.1221 4.1611 4.6117 5.0410 -3.0186 -7.0000 -3.0182 -3.0182 -7.0000 -3.0186 2、应力强度因子分布如图6所示 图6 表2 左边裂缝L2保持不变，右边裂缝L1变化时，右边裂缝的应力强度因子 L1 0 0.001 0.01 0.1 0.3 0.5 1 1.5 2 K NaN 0.14 0.4398 1.3063 2.0131 2.3663 2.82 3.0715 3.2540 表3 右边裂缝L1保持不变，左边裂缝L2变化时，右边裂缝的应力强度因子 L2 0 0.5 1 1.5 2 2.5 3 K 3.2950 3.3330 3.4075 3.4949 3.5872 3.6811 3.7750 四、结论 对于无限大板带裂纹的圆孔口问题, 本文利用复变函数方法将复杂的孔洞问题简单化以求其解, 并求出了裂纹端点的应力强度因强度因子表达式。该问题的有效解决进一步证明了复变函数方法的巧妙。本文应用等效裂纹的观念计算孔边裂纹的应力，其结果与传统方法的结果一致