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

工程中的弹塑性力学读书报告

作业：基于ANSYS对孔口应力集中问题进行简单分析

学号：S********

基于ANSYS对孔口应力集中问题进行简单分析

摘要：徐芝纶所著《弹性力学》给出矩形薄板左右受均布拉力的基尔斯解答，本文利用有限元软件ANSYS数值实验，实现应力场、位移场的可视化。同时定义了应力集中的特征参数来研究应力集中系数与孔径尺度（宽径比、长宽比）、材料所处状态的关系，最后提出一种可应用于工程中减小应力集中的方法。

关键词：圆孔应力集中特征参数应力集中系数孔径尺度材料状态

1引言

孔口的尺寸远小于弹性体的尺寸，并且孔边距弹性体的边界比较远

（约大于1.5倍孔口尺寸）被定义为小孔口问题。由于开孔，孔口附近的应力远大于无孔时的应力，也远大于距孔口较远处的应力，这种现象称为孔口应力集中。一般来说，集中的程度越高，集中的现象越是局部性的，就是说应力随着与孔的距离增大而越快的趋进于无孔时的应力。应力集中的程度，首先与孔的形状有关，一般来说，圆孔孔边的集中程度最低。另外集中系数还与相对孔径尺度有关。基于ansys平台，通过数值试验的方法，研究不同孔径时的孔边应力集中问题。

2力学模型

参考《弹性力学》和书籍［3］我们建立如下模型，如图1所示。

平面带孔平板，孔位于板正中，假设板为各向同性完全弹性，板两端受均布拉力荷载q=1000pa，长为200mm，宽为100mm，厚为0.01mm，泊松比为0.3，E=220Gpa。

现在我们定义一个描述板宽与孔径的相对尺度的特征参数，

￡1=%点2=%，定义应力集中系数k=5m^，其中L为板长，B为板‘0

宽，R为孔半径，smax为孔边最大应力，q为均布荷载，q0为平均应力。

q二0二q

图1力学模型

3ANSY，求解

3.1模型建立

为了研究不同孔径时的孔边应力集中，我们做一下数值实验。

B=100,R=66.67、50、33.33、25、20、16.67、10、5（单位mm），

利用ANSYS平台由￡值，计算出相应的k值。本模型为轴对称平面应力问题，所以我们只需要取％模型进行求解即可。具体操作过程这里不用细讲，生成R=10mm的几何模型和变形后模型见图2、3。

图2生成的几何模型结果显示

hbu? AMSYS

DECT部 皿牌

TDE-1

劾E■的

图3变形后的几何形状和未变形轮廓显示

3.2数据分析

3.2.1圆孔附近应力场

从R=10mm圆孔的X方向应力场分布等值线图（图4、5）,我们发现圆孔附近应力发生了改变：由平均1000pa急剧增加到2781?3137pa，即发生了明显的应力集中现象；在远离孔边应急急剧减少至1003pa；孔附近的应力远大于较远处的应力，且最大应力（图中红色区域）在圆孔的上下边。而《弹性力学》基尔斯解为§=q1+1土+3二，

中2P22P4J

当p=r,%=3q，并且随着远离孔边应力急剧趋近于q，比较数值解与解析解，两者基本吻合，所以我们通过ANSYS来实现解析解的可视化是成功的。

图4R=10mm和R=5mm圆孔X方向应力

图5R=10mm和R=5mm圆孔沿Y方向环向正应力分布

3.2圆孔附近的形变场

笔者查阅不少资料，发现文献大多是研究圆孔附近的应力场，而对圆孔附近的形变场阐述很少，下面笔者就用ANSYS来简单分析圆

孔附近形变场。

《弹性力学》已经给出矩形薄板左右受均布拉力的基尔斯解（1），我们把该解带入平面问题下的物理方程（2），之后在代入几何方程（3），就得到位移场的通解［1］（4）。

基尔斯解:5p=2rr2)1-一P2)r

基尔斯解:

5p=2

rr2)

1-一

P2)

r2

Yr2)

1-3——

P2qr

-qcos2中1+3

2

qr

t=—Isin2^1-

（1）

平面问题极坐标下的物理方程:

P4)

Y r2)

1+3一

人P2)

1c c

七=-(七-叫)

1c C

%=ER-Ep)

Y=2(1+^)pwE

（2）

平面问题极坐标下的几何方程:

du

8=—

p dp

u1du

8w=节+6就

1du du

ypw=PaWp+aP^

板内位移场通解:

t

pw

u

-_W

P

（3）

。一 r2 … .r2 _ — 、r4 _一

u——[(1一H)p+(1+H)—+(1+H)Pcos2^+4—cos2^—(1+H)—cos2^]

p 2E P P P3I(4

。「八、.――八、r2.一