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

工程中的弹塑性力学读书报告作业：基于ANSYS对孔口应力集中问题进行简单分析姓名:郭政学号：于ANSYS对孔口应力集中问题进行简单分析摘要：徐芝纶所著《弹性力学》给出矩形薄板左右受均布拉力的基尔斯解答，本文利用有限元软件ANSYS数值实验，实现应力场、位移场的可视化。同时定义了应力集中的特征参数来研究应力集中系数与孔径尺度(宽径比、长宽比)、材料所处状态的关系，最后提出一种可应用于工程中减小应力集中的方法。关键词：圆孔应力集中特征参数应力集中系数孔径尺度 材料状态1引言孔口的尺寸远小于弹性体的尺寸，并且孔边距弹性体的边界比较远（约大于1.5倍孔口尺寸）被定义为小孔口问题。由于开孔，孔口附近的应力远大于无孔时的应力，也远大于距孔口较远处的应力，这种现象称为孔口应力集中。一般来说，集中的程度越高，集中的现象越是局部性的，就是说应力随着与孔的距离增大而越快的趋进于无孔时的应力。应力集中的程度，首先与孔的形状有关，一般来说，圆孔孔边的集中程度最低。另外集中系数还与相对孔径尺度有关。基于ansys平台，通过数值试验的方法，研究不同孔径时的孔边应力集中问题。2 力学模型参考《弹性力学》和书籍[3]我们建立如下模型，如图1所示。平面带孔平板，孔位于板正中，假设板为各向同性完全弹性，板两端受均布拉力荷载，长为200mm，宽为100mm，厚为0.01mm，泊松比为0.3，。现在我们定义一个描述板宽与孔径的相对尺度的特征参数，，定义应力集中系数，其中L为板长，B为板宽，R为孔半径，为孔边最大应力，q为均布荷载，q0为平均应力。图1 力学模型3ANSYS求解3.1 模型建立为了研究不同孔径时的孔边应力集中，我们做一下数值实验。B=100，R=66.67、50、33.33、25、20、16.67、10、5（单位mm），利用ANSYS平台由值，计算出相应的值。本模型为轴对称平面应力问题，所以我们只需要取模型进行求解即可。具体操作过程这里不用细讲，生成的几何模型和变形后模型见图2、3。图2 生成的几何模型结果显示图3 变形后的几何形状和未变形轮廓显示3.2 数据分析3.2.1圆孔附近应力场从圆孔的X方向应力场分布等值线图（图4、5），我们发现圆孔附近应力发生了改变：由平均急剧增加到，即发生了明显的应力集中现象；在远离孔边应急急剧减少至；孔附近的应力远大于较远处的应力，且最大应力（图中红色区域）在圆孔的上下边。而《弹性力学》基尔斯解为，当，并且随着远离孔边应力急剧趋近于，比较数值解与解析解，两者基本吻合，所以我们通过ANSYS来实现解析解的可视化是成功的。图4 和圆孔X方向应力图5 和圆孔沿Y方向环向正应力分布3.2 圆孔附近的形变场笔者查阅不少资料，发现文献大多是研究圆孔附近的应力场，而对圆孔附近的形变场阐述很少，下面笔者就用ANSYS来简单分析圆孔附近形变场。《弹性力学》已经给出矩形薄板左右受均布拉力的基尔斯解（1），我们把该解带入平面问题下的物理方程（2），之后在代入几何方程（3），就得到位移场的通解[1]（4）。基尔斯解：（1）平面问题极坐标下的物理方程：（2）平面问题极坐标下的几何方程：（3）板内位移场通解：（4）当时,运用坐标转换公式, 即可获得在远离孔边弹性薄板内的位移分布: (5)该结果与无圆孔时单向拉伸的位移解完全一致。图5为圆孔附近的应变场，可见圆孔受拉由圆形变成了椭圆形，位移场成线性条带状分布，远离小孔位移值逐渐变大，在两端处达到最大，这完全符合位移公式（5）的描述。图5 圆孔附近的应变场3.3对应力集中的影响当孔径变化，应力集中现象会有什么改变呢？通过对比图5不同孔径下等效应力场等值线，我们得出结论：孔径越小应力，应力集中现象越明显，应力突然变大的趋势越迅速。在图5中的表现为蓝色区域所占弹性体比例越来越大。图5 和圆孔等效应力场等值线这里笔者不禁想讨论孔径的大小对应力集中影响到底有多大，我们统计上面8组试验的相对尺度与应力集中系数，得到如下结果，如表1。表1 应力集中系数与相对尺度之间的关系1.52345610206.3404.3153.3093.0943.0023.0003.0003.000将之间的关系，我们可以推论: 只要弹性体的边界离开圆孔中心有足够的距离，或者相对尺度，应力集中系数趋近稳定，也就是当板宽一定时，随着小孔孔径变小，小孔对应力集中的影响较小。3.4对应力集中的影响为了研究长宽比对应力集中系数k的影响，我们规定小孔R=10mm，然后利用ANSYS命令流（见附录）依次改变长宽，并记录应力值，计算出应力集中系数，得出的关系如图6。从图可见，当孔径一定时候，长宽比也影响应力集中。其中当，应力集中系数趋于稳定；当，长宽比越小，