断裂与损伤力学大作业一wxp.docVIP

应力强度因子数值计算方法综述 《断裂与损伤力学》大作业一 应力强度因子数值计算方法综述 目录 TOC \o 1-3 \t \h \z \u HYPERLINK \l _Toc3697 1 边界元法 PAGEREF _Toc3697 1 HYPERLINK \l _Toc17827 2 边界配位解法 PAGEREF _Toc17827 5 HYPERLINK \l _Toc7192 3 J积分求解裂纹尖端的应力强度因子 PAGEREF _Toc7192 7 HYPERLINK \l _Toc2892 4 裂纹尖端应力强度因子的能量差率解法 PAGEREF _Toc2892 9 HYPERLINK \l _Toc29161 5 裂纹尖端应力强度因子的解析变分解法 PAGEREF _Toc29161 12 HYPERLINK \l _Toc15767 参考文献 PAGEREF _Toc15767 15 PAGE \* MERGEFORMAT PAGE \* MERGEFORMAT 7 PAGE \* MERGEFORMAT PAGE \* MERGEFORMAT 0 应力强度因子数值计算方法综述 应力强度因子分析是断裂力学的重要内容之一。一方面要根据构件的几何形状、裂纹尺寸、结构的材料和受载方式去计算构件的 、和值。另一方面，当用标准试样测定材料断裂韧性时，先要确定试件应力强度因子标定公式。 过去四十年内，已经发展了很多有效的确定应力强度因子方法。这些方法总体上可以分为三类：解析法、数值解法和实验方法。对于三维裂纹问题和几何形状比较复杂的二维裂纹问题，往往要用数值解法。 本文将介绍裂纹尖端应力强度因子的一些主要数值计算方法，并阐述了每种数值算法的原理、适用范围以及优缺点等。 1 边界元法 通过适当方法（变分法、功的互等定理与加权余量法），可将弹性力学平面问题归结为求解如下积分方程（无体力情况）： （1-1） （1-2） 式（1-2）是无限域内平面问题的基本解。他们的物理意义是：表示在P点沿x方向作用一单位集中力，在R点沿y方向产生的位移；表示在P点沿x方向作用一单位集中力，在R点沿y方向引起的位移；其他带 “*”号的符号有类似的物理意义。式（1-2）中的r是P与R之间的距离 （1-3） 式中，（x,y）是P点的坐标；是R点的坐标见图。在式中，P点在边界上，R点在边界C内；但是在式（1-1）中，R点趋于边界上的点，即，。需要注意，式（1-1）中是无体积力情况下的边界积分方程；当存在体积力时，积分方程要加一项面积分。 边界积分方程（1-1）的精确解很难得到，必须进行数值计算。数值计算的方法是将边界积分方程离散化，划分边界C为M个单元，用各单元上的节点参数（位移和力）来表示单元上的相应量。例如，边界C被划为个单元，单元的节点取为，如图1-2所示。最简单的情况是讲各单元上节点的位移和力近似为常数，则边界积分方程（1-1）写成 （1-4） 类似地得到的表达式。 图1.1 图 1.2 由式（1-4）看出，对于每个节点可得两个线性代数方程，如果有M个节点，就得到2M个方程的代数方程组。若要有解，在边界上只能有2M个未知量。事实上，对于平面问题，无论在位移边界上或在力的边界上，边界上每个节点只有两个未知数；或者是两个位移分量，或者是两个力的分量。因此，在M个节点上恰好有2M个未知量，问题可解。 有时为了用较少的单元得到较精确的数值解，不一定设单元的力和位移是常量，可以通过等参变换，将各单元的位移和力用形状函数和各节点值表示。 在用边界元法求裂缝尖端应力强度因子时，为了保持应力具有的奇异性和位移具有的量级，在裂纹尖端可用如图1-3所示的点作为中间节点的裂纹元。在进行计算时，可以将含裂纹表面的裂纹体划出内部边界进行计算，但是在内部边界上必须满足位移的连续性和表面力的平衡，即 （1-5） 图 1.3 举例说明，带斜裂纹的一块有限宽板受拉伸应力的作用，裂纹长度是2a，斜角为，板宽为3W，如图1-4所示。内部边界与外部边界的单元划分见图1-5。其计算结果以及同一问题用边界配置法和有限元法所得的无量纲应力强度因子列于表1-1。 图 1.4 表1.1 无量纲应力强度因子表 边界配置法 有限元法 边界元法 0.730 0.728 0.725 0.600 0.590 0.5