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

1、工程中的弹塑性力学读书报告作业：基于ANSYS对孔口应力集中问题进行简单分析学号：s*基于ANSYS对孔口应力集中问题进行简单分析摘要：徐芝纶所著弹性力学给出矩形薄板左右受均布拉力的基尔斯解答，本文利用有限元软件ANSYS数值实验，实现应力场、位移场的可视化。同时定义了应力集中的特征参数来研究应力集中系数与孔径尺度（宽径比、长宽比）、材料所处状态的关系，最后提出一种可应用于工程中减小应力集中的方法。关键词：圆孔应力集中特征参数应力集中系数孔径尺度材料状态1引言孔口的尺寸远小于弹性体的尺寸，并且孔边距弹性体的边界比较远（约大于1.5倍孔口尺寸）被定义为小孔口问题。由于开孔，孔口附近的应力远大于无

2、孔时的应力，也远大于距孔口较远处的应力，这种现象称为孔口应力集中。一般来说，集中的程度越高，集中的现象越是局部性的，就是说应力随着与孔的距离增大而越快的趋进于无孔时的应力。应力集中的程度，首先与孔的形状有关，一般来说，圆孔孔边的集中程度最低。另外集中系数还与相对孔径尺度有关。基于ansys平台，通过数值试验的方法，研究不同孔径时的孔边应力集中问题。2力学模型参考弹性力学和书籍3我们建立如下模型，如图1所示。平面带孔平板，孔位于板正中，假设板为各向同性完全弹性，板两端受均布拉力荷载q二looopa，长为200mm，宽为100mm,厚为0.01mm，泊松比为0.3，E=220Gpa。现在我们定义一

3、个描述板宽与孔径的相对尺度的特征参数，81二BR,E2二Lb，定义应力集中系数k=5max：，其中L为板长，B为板宽，R为孔半径，5为孔边最大应力，q为均布荷载，q0为平均应力。max0Ep2）Yp2(1+卩)tEp平面问题极坐标下的几何方程：duu1duTOC o 1-5 h zs=p+ppd1duduup+pddpp板内位移场通解：(4)TOC o 1-5 h zqr2r2r4u=(1H)P+(1+H)+(1+H)Pcos2申+4cos2申一(1+p)cos2申p2Eppp3q.r2.r4.u=(1+p)psin2申一2(1p)sin2申一(1+p)sin2申2Epp3当Pr时，运用坐标转

4、换公式，即可获得在远离孔边弹性薄板内的位移分布:该结果与无圆孔时单向拉伸的位移解完全一致。图5为R-10mm圆孔附近的应变场，可见圆孔受拉由圆形变成了椭圆形，位移场成线性条带状分布，远离小孔位移值逐渐变大，在两端处达到最大，这完全符合位移公式（5）的描述。图5圆孔附近的应变场3.31对应力集中的影响当孔径变化，应力集中现象会有什么改变呢？通过对比图5不同孔径下等效应力场等值线，我们得出结论：孔径越小应力，应力集中现象越明显，应力突然变大的趋势越迅速。在图5中的表现为蓝色区域所占弹性体比例越来越大。图5R=10mm和R=5mm圆孔等效应力场等值线这里笔者不禁想讨论孔径的大小对应力集中影响到底有多

5、大，我们统计上面8组试验的相对尺度与应力集中系数k，得到如下结果，如表1。表1应力集中系数k与相对尺度1之间的关系11.5234561020k6.3404.3153.3093.0943.0023.0003.0003.000将1_k之间的关系，我们可以推论：只要弹性体的边界离开圆孔中心有足够的距离，或者相对尺度14，应力集中系数趋近稳定，也就是当板宽一定时，随着小孔孔径变小，小孔对应力集中的影响较小。3.42对应力集中的影响为了研究长宽比2对应力集中系数k的影响，我们规定小孔R=10mm，然后利用ANSYS命令流（见附录）依次改变长宽L/B=1.04.0，并记录应力值，计算出应力集中系数，得出2

6、k的关系如图6。从图可见，当孔径一定时候，长宽比2也影响应力集中。其中当22.0，应力集中系数趋于稳定；当2=1.02.0，长宽比2越小，应力集中系数越大。图6s2k3.5材料所处的状态对应力集中的影响笔者参阅相关资料2,发现目前对孔边应力集中问题的研究大多停留在线弹性领域内，并没有反应材料处于屈服阶段、强化阶段的孔边应力集中规律，所以笔者通过ANSYS建模来找出这一规律。我们假设薄板为各向同性理想的弹塑性材料，泊松比为0.3，弹性模量E=220Gpa，L=200mm,B=100mm,R=10mm，材料的屈服极限5=240Mpa，应力集中系数k=5ma,/，为孔边最大应力，q为均布sqmax0

7、荷载,q为平均应力。由于研究材料产生屈服前、后的应集中现象，0所以每一个实验采用2中荷载：一种荷载材料未进入屈服状态，另一种荷载进入屈服，实验数据见下表。q5maxq0k80250.36788.8122.819067100249.68189.832.779483120261.2394.3542.768616140274.522100.0822.742971160289.537105.392.747291180308.319112.3352.744639200339.169121.832.783953210379.133134.7342.8139372.8203522.6972622.53370

8、7220230240422.438465.631508.097149.782172.631200.5352.852.802.752.702.652.602.552.506080100120140160180200220240260q图6qk图7两端均布荷载80Mpa、240Mpa作用下的等效应力图从上表和图6可见，均布荷载q=80MPa，材料尚未进入屈服阶段，孔边存在明显的应力集中，且随着荷载逐渐增大，在一定范围内应力集中系数基本上不变（k2.78）；当均布荷载逐渐接近屈服极限，应力集中系数迅速下降，设想如果荷载再继续增大，远离孔边的材料都会相继进入屈服状态，尽管变形在不断增加，应力却停留在屈

9、服极限附近不再增加5。4减少薄板应力集中的方法通过上面分析，我们知道孔径越小，应力集中越明显。而在实际工程中，有时要求孔径不变且使得应力集中较小，这就需要我们探究新的方法来减少中心孔边的应力集中。通过参阅相关文献4,笔者设想在中心圆孔两侧开对称小孔，已达到减少中心圆孔应力集中，如图7所示。MUJI3GP图7圆孔两侧对称开小孔模型我们设定中心小孔R=10mm，两侧小孔均为R=5mm。通过ANSYS实验得到等效应力云图如图7所示。可见最大应力位置依然在圆孔上下端，但最大应力值为2870Pa，较之单孔平板最大应力值3137Pa，降低了8.51%，说明此方法对减少薄板圆孔应力集中是有效的。至于如何调整

10、两侧小孔的参数（孔径、孔的相对位置），已实现最小应力集中效果。鉴于笔者初学ANSYS,也不懂复变函数，无法给出相关解析解，这是本报告的遗憾。4结论（1）利用ANSYS求解与弹性力学所给解析解基本吻合，并且可以实现应力场、位移场可视化。当相对尺度814时，应力集中系数k趋近稳定。即当板宽一定时，随着孔径尺寸变小，应力集中系数趋近一个稳定值，小孔对单向受拉薄板应力集中影响较小。孔径一定，长宽比82也影响应力集中。当822.0，应力集中系数趋于稳定；当82二1.02.0，长宽比82越小，应力集中系数越大。(4)应力集中与材料所处状态有关。材料应力小于屈服极限，应力集中系数基本上稳定，应力接近屈服极限

11、时，应力集中系数突然下降。(5)在中心圆孔两侧开对称小孔可以有效减少应力集中。参考文献刘述伦等.应力集中圆孔附近的形变场J.安徽建筑工业学院，2011(19):23-24.徐芝纶编.弹性力学(上册).M.第四版.北京.高等教育出版社.2006陈章华.工程中的有限元分析方法M.北京.冶金工业出版社.20134朱晓东等.基于ANSYS平台含圆孔薄板的应力集中分析J.苏州大学学报，2004(24):51-52.银光球等.有限板宽孔边应力集中有限元分析J.福州工程学院学报.2009(7):67-68.附录小孔应力集中建模所使用的命令流文件/NOPR!SuppressprintingofUNDOproc

12、ess/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,0KEY

13、W,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,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,1.08222638492,1/REP,FASTFINISH/SOL!*ANTYPE,0FINISH/PREP7ESIZE,0.002,0,!*AMAP,1,1,4,5,3!*!*ANTYPE,0FLST,2,1,4,ORDE