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

《断裂与损伤力学》大作业一 应力强度因子数值计算方法综述 目录 1 边界元法 1 2 边界配位解法 5 3 J积分求解裂纹尖端的应力强度因子 7 4 裂纹尖端应力强度因子的能量差率解法 9 5 裂纹尖端应力强度因子的解析变分解法 12 参考文献 15  应力强度因子数值计算方法综述 应力强度因子分析是断裂力学的重要内容之一。一方面要根据构件的几何形状、裂纹尺寸、结构的材料和受载方式去计算构件的 、和值。另一方面，当用标准试样测定材料断裂韧性时，先要确定试件应力强度因子标定公式。 过去四十年内，已经发展了很多有效的确定应力强度因子方法。这些方法总体上可以分为三类：解析法、数值解法和实验方法。对于三维裂纹问题和几何形状比较复杂的二维裂纹问题，往往要用数值解法。 本文将介绍裂纹尖端应力强度因子的一些主要数值计算方法，并阐述了每种数值算法的原理、适用范围以及优缺点等。 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.598 误差 — -0.27 -0.68 误差 — -1.67 -0.33 — 48 19 边界元法的优越性在于缩小工作量，减少计算机的计算量。 2 边界配位解法 若要获得含裂纹有限大板应力场的解析解，通常在数学上有很大困难。因此，许多研究者正在探索有限尺寸板裂纹问题的半解析、半数值解。 边界配位法就是一种求解有限大板应力强度因子的一种半解析、半数值方法。这种方法在求解含边缘裂纹的单连通板问题时还是比较方便的。用于材料断裂性能测试的三点弯曲试样与紧凑拉伸试样的应力强度因子就是由这种方法求得的。边界配位解法是建立在平面问题极坐标分离变量解法原理之上的。 为了根据边界条件与附加的简化条件确定应力函数A中与应力强度因子直接相关的待定系数c1和d1，列出应力函数及其一阶偏导数的表达式。 （2-1