弹性力学小孔应力集中读书报告

工程中的弹塑性力学读书报告作业：基于ANSYS对孔口应力集中问题进行简单分析姓名:郭政学号：S20156116基于ANSYS对孔口应力集中问题进行简单分析摘要：徐芝纶所著《弹性力学》给出矩形薄板左右受均布拉力的基尔斯解答，本文利用有限元软件ANSYS数值实验，实现应力场、位移场的可视化。同时定义了应力集中的特征参数来研究应力集中系数与孔径尺度(宽径比、长宽比)、材料所处状态的关系，最后提出一种可应用于工程中减小应力集中的方法。关键词：圆孔应力集中特征参数应力集中系数孔径尺度材料状态1引言孔口的尺寸远小于弹性体的尺寸，并且孔边距弹性体的边界比较远（约大于1.5倍孔口尺寸）被定义为小孔口问题。由于开孔，孔口附近的应力远大于无孔时的应力，也远大于距孔口较远处的应力，这种现象称为孔口应力集中。一般来说，集中的程度越高，集中的现象越是局部性的，就是说应力随着与孔的距离增大而越快的趋进于无孔时的应力。应力集中的程度，首先与孔的形状有关，一般来说，圆孔孔边的集中程度最低。另外集中系数还与相对孔径尺度有关。基于ansys平台，通过数值试验的方法，研究不同孔径时的孔边应力集中问题。2力学模型参考《弹性力学》和书籍[3]我们建立如下模型，如图1所示。平面带孔平板，孔位于板正中，假设板为各向同性完全弹性，板两端受均布拉力荷载1000qpa，长为200mm，宽为100mm，厚为0.01mm，泊松比为0.3，220EGpa。现在我们定义一个描述板宽与孔径的相对尺度的特征参数，1,2BLRB，定义应力集中系数max0kq，其中L为板长，B为板宽，R为孔半径，max为孔边最大应力，q为均布荷载，q0为平均应力。图1力学模型Rqq3ANSYS求解3.1模型建立为了研究不同孔径时的孔边应力集中，我们做一下数值实验。B=100，R=66.67、50、33.33、25、20、16.67、10、5（单位mm），利用ANSYS平台由值，计算出相应的k值。本模型为轴对称平面应力问题，所以我们只需要取14模型进行求解即可。具体操作过程这里不用细讲，生成10Rmm的几何模型和变形后模型见图2、3。图2生成的几何模型结果显示图3变形后的几何形状和未变形轮廓显示3.2数据分析3.2.1圆孔附近应力场从10Rmm圆孔的X方向应力场分布等值线图（图4、5），我们发现圆孔附近应力发生了改变：由平均1000pa急剧增加到27813137pa，即发生了明显的应力集中现象；在远离孔边应急急剧减少至1003pa；孔附近的应力远大于较远处的应力，且最大应力（图中红色区域）在圆孔的上下边。而《弹性力学》基尔斯解为242413122rrq，当,3rq，并且随着远离孔边应力急剧趋近于q，比较数值解与解析解，两者基本吻合，所以我们通过ANSYS来实现解析解的可视化是成功的。图410Rmm和5Rmm圆孔X方向应力图510Rmm和5Rmm圆孔沿Y方向环向正应力分布3.2圆孔附近的形变场笔者查阅不少资料，发现文献大多是研究圆孔附近的应力场，而对圆孔附近的形变场阐述很少，下面笔者就用ANSYS来简单分析圆孔附近形变场。《弹性力学》已经给出矩形薄板左右受均布拉力的基尔斯解（1），我们把该解带入平面问题下的物理方程（2），之后在代入几何方程（3），就得到位移场的通解[1]（4）。基尔斯解：222222242422221cos2113221cos21322sin21132qrqrrqrqrqrr（1）平面问题极坐标下的物理方程：1()1()2(1)EEE（2）平面问题极坐标下的几何方程：11uuuuuu（3）板内位移场通解：2243243[(1)(1)(1)cos24cos2(1)cos2]2[(1)sin22(1)sin2(1)sin2]2qrrruEqrruE（4）当r时,运用坐标转换公式,即可获得在远离孔边弹性薄板内的位移分布:quxEqyvE(5)该结果与无圆孔时单向拉伸的位移解完全一致。图5为10Rmm圆孔附近的应变场，可见圆孔受拉由圆形变成了椭圆形，位移场成线性条带状分布，远离小孔位移值逐渐变大，在两端处达到最大，这完全符合位移公式（5）的描述。图5圆孔附近的应变场3.31对应力集中的影响当孔径变化，应力集中现象会有什么改变呢？通过对比图5不同孔径下等效应力场等值线，我们得出结论：孔径越小应力，应力集中现象越明显，应力突然变大的趋势越迅速。在图5中的表现为蓝色区域所占弹性体比例越来越大。图510Rmm和5Rmm圆孔等效应力场等值线这里笔者不禁想讨论孔径的大小对应力集中影响到底有多大，我们统计上面8组试验的相对尺度与应力集中系数k，得到如下结果，如表1。表1应力集中系数k与相对尺度1之间的关系11.5234561020k6.3404.3153.3093.0943.0023.0003.0003.000将1k之间的关系，我们可以推论:只要弹性体的边界离开圆孔中心有足够的距离，或者相对尺度14，应力集中系数趋近稳定，也就是当板宽一定时，随着小孔孔径变小，小孔对应力集中的影响较小。3.42对应力集中的影响为了研究长宽比2对应力集中系数k的影响，我们规定小孔R=10mm，然后利用ANSYS命令流（见附录）依次改变长宽/1.04.0LB，并记录应力值，计算出应力集中系数，得出2k的关系如图6。从图可见，当孔径一定时候，长宽比2也影响应力集中。其中当22.0，应力集中系数趋于稳定；当21.02.0，长宽比2越小，应力集中系数越大。1.01.52.02.53.03.54.03.123.143.163.183.203.223.243.263.28kL/BR=10mm图62k3.5材料所处的状态对应力集中的影响笔者参阅相关资料[2],发现目前对孔边应力集中问题的研究大多停留在线弹性领域内，并没有反应材料处于屈服阶段、强化阶段的孔边应力集中规律，所以笔者通过ANSYS建模来找出这一规律。我们假设薄板为各向同性理想的弹塑性材料，泊松比为0.3，弹性模量220EGpa,200,100,10LmmBmmRmm,材料的屈服极限240sMpa，应力集中系数max0kq,max为孔边最大应力，q为均布荷载,0q为平均应力。由于研究材料产生屈服前、后的应集中现象，所以每一个实验采用2中荷载：一种荷载材料未进入屈服状态，另一种荷载进入屈服，实验数据见下表。qmaxq0k80250.36788.8122.819067100249.68189.832.779483120261.2394.3542.768616140274.522100.0822.742971160289.537105.392.747291180308.319112.3352.744639200339.169121.832.783953210379.133134.7342.813937220422.438149.7822.820352230465.631172.6312.697262240508.097200.5352.53370760801001201401601802002202402602.502.552.602.652.702.752.802.85kq图6qk图7两端均布荷载80Mpa、240Mpa作用下的等效应力图从上表和图6可见，均布荷载q=80MPa，材料尚未进入屈服阶段，孔边存在明显的应力集中，且随着荷载逐渐增大，在一定范围内应力集中系数基本上不变（2.78k）；当均布荷载逐渐接近屈服极限，应力集中系数迅速下降，设想如果荷载再继续增大，远离孔边的材料都会相继进入屈服状态，尽管变形在不断增加，应力却停留在屈服极限附近不再增加[5]。4减少薄板应力集中的方法通过上面分析，我们知道孔径越小，应力集中越明显。而在实际工程中，有时要求孔径不变且使得应力集中较小，这就需要我们探究新的方法来减少中心孔边的应力集中。通过参阅相关文献[4],笔者设想在中心圆孔两侧开对称小孔，已达到减少中心圆孔应力集中，如图7所示。图7圆孔两侧对称开小孔模型我们设定中心小孔R=10mm，两侧小孔均为R=5mm。通过ANSYS实验得到等效应力云图如图7所示。可见最大应力位置依然在圆孔上下端，但最大应力值为2870Pa，较之单孔平板最大应力值3137Pa，降低了8.51%，说明此方法对减少薄板圆孔应力集中是有效的。至于如何调整两侧小孔的参数（孔径、孔的相对位置），已实现最小应力集中效果。鉴于笔者初学ANSYS，也不懂复变函数，无法给出相关解析解，这是本报告的遗憾。4结论（1）利用ANSYS求解与弹性力学所给解析解基本吻合，并且可以实现应力场、位移场可视化。（2）当相对尺度14时，应力集中系数k急剧增大；当相对尺度14时，应力集中系数k趋近稳定。即当板宽一定时，随着孔径尺寸变小，应力集中系数趋近一个稳定值，小孔对单向受拉薄板应力集中影响较小。（3）孔径一定，长宽比2也影响应力集中。当22.0，应力集中系数趋于稳定；当21.02.0，长宽比2越小，应力集中系数越大。（4）应力集中与材料所处状态有关。材料应力小于屈服极限，应力集中系数基本上稳定，应力接近屈服极限时，应力集中系数突然下降。（5）在中心圆孔两侧开对称小孔可以有效减少应力集中。参考文献[1]刘述伦等.应力集中圆孔附近的形变场[J].安徽建筑工业学院，2011(19):23-24.[2]徐芝纶编.弹性力学(上册).[M].第四版.北京.高等教育出版社.2006[3]陈章华.工程中的有限元分析方法[M].北京.冶金工业出版社.2013[4]朱晓东等.基于ANSYS平台含圆孔薄板的应力集中分析[J].苏州大学学报，2004(24):51-52.[5]银光球等.有限板宽孔边应力集中有限元分析[J].福州工程学院学报.2009(7):67-68.附录小孔应力集中建模所使用的命令流文件/NOPR!SuppressprintingofUNDOprocess/PMACRO!EchofollowingcommandstologFINISH!MakesureweareatBEGINlevel/CLEAR,NOSTART!ClearmodelsincenoSAVEfound!WESUGGESTYOUREMOVETHISLINEANDTHEFOLLOWINGSTARTUPLINES!*!*/NOPR/PMETH,OFF,0KEYW,PR_SET,1KEYW,PR_STRUC,1KEYW,PR_THERM,0KEYW,PR_FLUID,0KEYW,PR_ELMAG,0KEYW,MAGNOD,0KEYW,MAGEDG,0KEYW,MAGHFE,0KEYW,MAGELC,0KEYW,PR_MULTI,0KEYW,PR_CFD,0/GO!*!*/PREP7!*ET,1,PLANE82!*!*MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,2.2e11MPDATA,PRXY,1,,0.3RECTNG,0,0.125,0,0.05,/DIST,1,1.08222638492,1/REP,FAST/DIST,1,1.08222638492,1/REP,FASTPCIRC,0.01,0,0,90,ASBA,1,2NUMCMP,ALL/DIST,1,0.924021086472,1/REP,FAST/DIST,1,0.924021086472,1/REP,FA