不同板宽的应力集中.doc

《有限元分析及应用》课程case study 不同板宽的孔边应力集中问题的数值分析 摘要：本文在ansy9.0平台上，对不同板宽的孔边应力集中问题进行了数值实验，定义了板宽与孔径的特征参数ξ，对ξ和应力集中系数K之间的关系进行分析，并验证了无限大板宽的解析解适用范围。 引言 受力弹性体具有小孔，则孔边的应力将远大于无孔时的应力，也远大于距孔稍远处的应力。这种现象称为孔边应力集中。应力集中的程度，首先与孔的形状有关。一般来说，圆孔孔边的集中程度最低。圆孔孔边应力集中有齐尔西的解答： 本文定义了一个板宽和孔径的特征参数ξ，对于平面圆孔的孔边问题，通过数值试验的方法研究孔边应力集中问题，进而分析ξ与应力集中系数K的关系。 计算模型 在ansys7.0平台上，我们把薄板的单元类型选择 solid quad 4 node 42，带厚度的薄板进行建模，板厚度定义为1，建立一块100*100薄板，在中间建立半径为5的圆，通过布尔运算得到带孔薄板，在板的一边我们定义横向的约束，在另一边我们施加均布荷载。最终我们可以得到如下的模型（图1），我们之后将以此模型为基础进行数值试验。 图1 ξ与K的关系研究 本文定义了板宽和孔径的特征参数ξ=，为改变ξ，本文定义不同板宽B， B的取值如表1所示： 表1 B 500 400 300 200 120 100 80 60 40 R0 5 5 5 5 5 5 5 5 5 ξ 100 80 60 40 24 20 16 12 8 我们定义了应力集中系数K= ① ，利用ansys，我们可以得到不同板宽下的（图2）， （图2） q0为无孔时的应力，利用公式①，我们可以获得不同板宽应力集中系数如表2： 表2 ξ 100 80 60 40 24 20 16 12 8 k 3.02545 3.0284 3.032 3.0362 3.0506 3.06215 3.08021 3.11789 3.25188 由表2，可以画出ξ与K的关系图 （图3） 可以看到，在ξ较小时，应力集中系数K大于3，随着ξ的增大，板的应力集中系数K逐渐趋近于3。由于板宽较大时候网格划分的密度对结果影响较大，所以可以认为K将无限趋近于3. 无限板宽孔边应力集中的解析解的适用范围 无限大板宽情形下的解析解： 在孔边的y轴上有((分布： 由无限板宽孔边应力集中的齐尔西解答，我们可以得到在孔边的应力应该为3q0.则由表2可以计算出不同板宽对应解析解的相对误差，本文利用数据制作了一条误差曲线。 图4 由曲线可以看出在板宽与孔径比在50之后时，误差曲线变得平缓，解析解比较吻合数值试验。 可以看到无限板宽的应力集中问题的齐尔西解答适应板宽与孔径比较大的情况。 结论 在ξ较小时，应力集中系数K大于3，随着ξ的增大，板的应力集中系数K逐渐趋近于3。 无限板宽的应力集中问题的齐尔西解答适应板宽与孔径比较大的情况。 参考资料 徐芝纶 弹性力学（ 第四版）高等教育出版社 有限单元法基础 王焕定 高等教育出版社 ansys有限元分析 邓凡平 人民邮电出版社 7.ansys数据流 /PREP7 !* ET,1,PLANE42 !* KEYOPT,1,1,0 KEYOPT,1,2,0 KEYOPT,1,3,3 KEYOPT,1,5,0 KEYOPT,1,6,0 !* !* R,1,1, !* !* MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,2.1e5 MPDATA,PRXY,1,,0.3 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,2.7 BLC4,0,0,100,400 CYL4,50,200,5 ASBA, 1, 2 ESIZE,0,50, MSHAPE,0,2D MSHKEY,0 !* CM,_Y,AREA ASEL, , , , 3 CM,_Y1,AREA CHKMSH,AREA CMSEL,S,_Y !* AMESH,_Y1 !* CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 !* FINISH /SOL /DIST,1,1.08222638492,1 /REP,FAST /DIST,1,1.08222638492,1 /REP,FAST /DIST,1,1.08222638492,1 /REP,FAST /DIST,1,1.08222638492,1