论文部分内容阅读
本文将有限元方法与偶应力理论相结合,建立了基于一般偶应力理论的有限元方法,并将其应用于线弹性断裂问题的分析。推导了有限元法中八节点四边形单元的形状函数,结合偶应力理论基本原理及有限元法变分原理推导了基于偶应力理论的有限元法基本原理。
采用MATLAB语言对基于偶应力理论的有限元法进行程序实现,并通过受均布荷载的悬臂梁及带小孔的无限大板两个典型的弹性力学算例验证了算法的正确性。计算结果表明:在不考虑偶应力的情况下,有限元法得到的结果与理论解基本吻合;在考虑偶应力的情况下,当特征长度与特征尺寸处于相近数量级时,尺度效应明显;带小孔的无限大板问题中小孔附近的应力集中系数相对经典理论解有所降低。
研究了Ⅰ型多裂纹断裂问题及复合型多裂纹断裂问题。结果表明:采用位移外推法得到的有限元解与理论解基本吻合;考虑偶应力影响下得到的Ⅰ型裂纹及复合型裂纹的Ⅰ、Ⅱ型应力强度因子计算结果小于经典理论解,随着特征长度与裂纹长度比值增大,计算结果减小;随着Ⅰ型平行多裂纹间距减小及裂纹长度减小,裂纹尖端的应力强度因子减小;复合型共线多裂纹中共线裂纹的存在使得裂纹间相对较近裂尖点的Ⅰ、Ⅱ型应力强度因子增大,对于裂纹间相对较远裂尖点的应力强度因子减小,随着共线裂纹间距离的增大,裂纹间相距较近点的Ⅰ、Ⅱ型应力强度因子有降低的趋势。
对基本偶应力理论的奇异单元在断裂力学中的应用进行了研究,推导了奇异单元的形函数。结果表明:在采用相近网格划分密度的情况下,采用奇异单元计算得到的Ⅰ型裂纹及复合型裂纹应力强度因子相对于采用平面八节点四边形单元得到的结果具有更高的精度。