第8章 平面问题的复变函数解

第八章 平面问题的复变函数解

知识点

双调和方程的复变函数表达形式 应力分量复变函数表达式 应力分量的单值条件 多连域的K-M函数 无穷远应力与K-M函数 位移分量的曲线坐标表达 保角变换公式与K-M 函数 柯西积分确定K-M 函数 孔口应力

裂纹前缘应力分布

双调和函数的复变函数形式 位移分量的复变函数表达形式 位移分量的单值条件

无限大多连域中K-M 函数的一般形式 保角变换和曲线坐标

应力分量的曲线坐标表达式

利用孔口边界条件确定K-M 函数 椭圆孔口的保角变换 裂纹—短轴为零的椭圆

切应力作用的裂纹前缘应力

一、内容介绍

通过直角坐标和极坐标系，可以求解一些弹性力学平面问题。但是，这些方法只能用于某些边界比较特殊的平面问题，特别是对于多连域问题更显得无能为力。

本章介绍复变函数解法，实质仍然是在给定的边界条件下求解双调和方程的问题，但应用中成为在给定的边界条件下寻找两个解析函数K-M函数的问题。求解分析步骤为：

1、分别将应力函数、应力分量、位移和边界条件等表示为复变函数形式，就是用K-M函数表示；

2、探讨无限大多连域中，K-M函数的表达形式，将其表示为级数形式； 3、利用保角变换将无限大多连域映射为单位圆，使得问题的边界条件简化； 4、将边界条件转化为柯西积分，求解级数系数，从而使得问题求解。 如果你还没有学习复变函数课程，请你学习附录2或者查阅有关参考资料。

二、重点

1、K-M函数与应力函数、应力分量、位移和边界条件等；2、无限大多连域的K-M函数形式；3、保角变换与曲线坐标；4、椭圆孔口与平面裂纹问题。

§8. 1 应力函数的复变函数表示 学习思路：

弹性力学应力解法的基本方程是双调和方程，问题求解的关键是建立满足边界条件的双调和函数。对于复变函数解，重要的问题是将双调和函数表达为复变

1

函数形式。

本节首先将双调和方程表示为复变函数形式；然后通过积分用解析函数表示双调和函数。学习时应该注意：应力函数为实函数，通过复变函数表达的双调和函数也是实函数，因此应力函数虚部等于零。

上述分析的结果是使得应力函数通过两个单值解析函数 和?(z)表示。

和?(z)称为克罗索夫－穆斯赫利什维利函数，简称K-M函数；或者称

为复位势函数。 学习要点：

1、双调和方程的复变函数表达形式；2、双调和函数的复变函数形式 1、双调和方程的复变函数表达形式

在弹性力学的复变函数求解中，应力函数用U（x，y）表示， 有其它定义。

设应力函数U（x，y）为双调和函数，首先考虑变形协调方程的复变函数表达形式。

对于复变函数z =x+ i y，取其共轭，则 =x- i y。因此z和 均为x，y的函数。复变函数z可以写作z=??ei?，其共轭 =??e-i?，因此z和 又可以表示为坐标??和? 的函数。

同理，x，y 也可以表示为z和 的函数，有

因此，应力函数也可以表示为复变函数z和的函数，有

注意到

应力函数U（x，y）对坐标x，y的导数也可以表示为对复变函数z和的求导运算，有

2

将上式的后两式相加，可以得到调和方程的复变函数表达形式

双调和方程的复变函数表达式为

2、双调和函数的复变函数形式

对于应力函数U（z）的复变函数表示。将双调和方程的复变函数表达式

乘以2，并对作积分，可得

对再作一次积分，可得

对z作一次积分，可得

对z再作积分一次，可得

应力函数U(z)的复变函数表达式中，有四个待定函数。注意到应力函数为实函数，因此公式右边的复变函数的虚部必须为零。所以上述函数必须是两两共轭的，即

3

或者

因此应力函数可以用两个待定函数表示为

或者

上述公式称为古尔萨（Goursat）公式。公式将双调和函数通过两个复变函数 和?(z)表达。 和?(z)称为克罗索夫－穆斯赫利什维利函数，简称K-M 函数，均为单值解析函数。

Re为表示复变函数实部的符号。 §8. 2 应力分量的复变函数表示 学习思路：

应力函数已经通过K-M函数表示，但是这还不够，为了下一步的工作，本节的工作是将应力分量表示为复变函数形式，即使用K-M函数表示应力分量。

这一工作的主要内容是写出K-M函数对直角坐标的偏导数，应该注意，本章应力分量表达式也是写作复变函数表达形式的。

本节引入复变函数

， 和

这主要是简化公式的描述，并没有增加未知函数。上述函数均称为K-M函数。 学习要点：

1、K-M函数对直角坐标导数的复变函数表示；2、应力分量表达式

1、K-M函数对直角坐标导数的复变函数表示

对于无体力的弹性力学问题。如果选取的应力分量满足

4

则上述应力分量自然是满足平衡微分方程的。这里的问题是选取的应力函数是复变函数形式表达的，而且是由K-M 函数描述的。因此，应力分量也必须通过K-M 函数表达。根据公式

有

将上述两式相加，可以得到

将上式分别对x和y求一阶导数，可得

其中

2、应力分量表达式

上述公式

的第一式减去第二式乘以 i，可得

5

即

将公式的第一式加上第二式乘以 i，可得

取其共轭，则

上述公式推导中，引入 和 。公式是用单值解析函数 和??(z) 表示的应力分量，自然满足平衡微分方程和变形协调方程。

§8. 3 位移的复变函数表示 学习要点与思路：

本节工作是将位移分量表示为复变函数形式，通过K-M 函数表达弹性体位移。

对于应力解法，如果应力函数满足变形协调方程，则单连域问题应力函数描述的变形已经是协调的。一般的讲，不需要专门分析位移。但是对于复变函数弹性力学解，处理的问题均为多连域问题，因此位移单值连续条件即使对于应力解法也是不可缺省的。

在位移的复变函数表达式推导中，首先将几何方程代入物理方程的前两式，得到位移偏导数的表达式，这是K-M 函数对x，y坐标的偏导数。积分确定位移表达式，并且根据物理方程的第三式，确定弹性体的刚体位移，写出K-M 函数表达的位移复变函数表达形式。 学习要点：

1、K-M 函数表达的位移偏导数表达式；2、积分确定位移分量；3、位移分量的复变函数表达形式

6

1、K-M 函数表达的位移偏导数表达式

对于平面应力问题，根据物理方程和几何方程的前两式，可得

其中

设 。由于K-M 函数 为解析函数，其连续可导，应满足柯西-黎曼（Cauchy-Riemann）条件，即

由于

取其共轭

因此可得

即

将上式代入公式

可得

7

2、积分确定位移分量

将公式

分别对x和y积分，可得

根据几何方程的第三式和切应力表达式代入切变胡克定理 ，上式可以写作

将位移表达式代入上式，则

整理可得

根据柯西-黎曼条件，公式左边第一项为零，所以

因此，g(x)=??x+ v0, f(x)=???y+ u0。这一结果说明：g(x)和f(y)表示刚体位移，因此对于变形分析，可以略去不计。即公式可以表示为

3、位移分量的复变函数表达形式

8

将上述两式

相加，则可得K-M 函数表示的位移分量。有

整理可得

或者写作

上述分析表明，如果已知K-M 函数 和??(z) 时，则平面应力状态下的位移分量也是确定的。

对于平面应变问题，只须将弹性模量和泊松比

做对应的替换则可。

§8.4 边界条件的复变函数表示 学习思路：

边界条件应用是弹性力学分析的重要步骤，本节讨论应用K-M 表示面力边界条件。由于应力和位移分量都是复变函数表示的，为方便进一步的分析，面力边界条件也需要用K-M 函数表达。

在直角坐标系中，边界条件是以函数形式表示的，对应于一点的边界条件。而在复变函数解中，更多使用边界线段的表达形式，这是复变函数性质决定的。

用复变函数描述的面力边界条件有三个。显然，这三个关系式不是独立的，仅有两个独立关系。

9

学习要点：

1、任意一点的面力边界条件复变函数表达； 2、边界线段AB的面力边界条件： 3、边界力矩与K-M 函数的关系： 4、位移边界条件： 思考题：

1、根据上述面力边界条件说明：对于单连域弹性体，K-M 函数为单值解析函数，而对于多连域，K-M 函数将不再是单值的。（解答） 1、任意一点的面力边界条件复变函数表达

对于弹性力学平面问题，其面力边界条件为

将复变函数表示的应力分量表达式

代入上式，则

设AB为弹性体的任意一段边界，而s是从边界上一点A量取的弧长（边界的外法线n指向弧长的右边），如图所示

10

则由几何关系

将上式代入公式

可得

将上述面力矢量用复数形式表达为Fsx+ iFsy，则

将公式

代入上式，可得

即

2、边界线段AB的面力边界条件

公式 的左边表示边界面力矢量在微分线段ds上的主矢量。将公式沿边界从定点A到动点B（设B点的坐标为z）积分，则可得边界面力矢量在弹性体边界线段AB上的主矢量

11

由于在K-M 函数 和??(z)中，增加或减少一个复常数并不影响应力值，因此可以适当的选取K-M 函数，使上式的常数为零。则

上述公式表示了边界面力矢量与K-M 函数 和??(z)之间的关系。 显然对于给定的面力矢量，公式的右边为边界点的确定的函数，即已知函数；而左边为坐标 z 从弹性体区域内部向区域的边界趋近时的复位势函数值。 3、边界力矩与K-M 函数的关系

如果将边界线段AB的面力矢量对坐标原点取矩，并利用关系式

可以得到

对上式作分部积分，可得

注意到

和

回代可得

公式的左边在外力给定的条件下，为边界点的确定函数。公式的右边为K-M 函数由弹性体内部向边界趋近时的数值。

12

4、位移边界条件

下面再讨论位移边界条件，当边界位移给定时，设边界位移为

u=u, v=v

则根据位移边界条件，有

上式即为K-M 函数表示的位移边界条件。

到此为止，求解弹性力学平面问题，由在给定的边界条件下求解双调和方程的问题，变换为在给定的边界条件下寻找解析函数 和??(z) 的问题。 思考题：

1、根据上述面力边界条件说明：对于单连域弹性体，K-M 函数为单值解析函数，而对于多连域，K-M 函数将不再是单值的。

解答：如讨论的物体为单连域，由于 和 均为单值解析函数。因此，当

A

和

B

重合时，也就是说积分曲线闭合时，则

＝0，这表明作用在物体边界上的

边界面力，必组成一个平衡力系。这一结论是必然的，要使问题有解，边界上的面力必须满足这个条件。

§8.5 多连域中φf(z)和?(z)的一般表达式 学习思路：

本节的主要任务是确保描述应力和位移分量的K-M 函数的单值性。 单连域中的单值解析函数 和??(z)在多连域可能不再是单值的。因此K-M 函数 和??(z)表示的应力和位移形式也可能不再单值。要保证应力和位移分量的单值性，必须讨论K-M 函数 和??(z)在多连域中的可能形式。

首先分别根据应力分量的单值条件，将K-M 函数 和?(z)分解为单值解析函数和可能的多值函数两部份，构造可能的K-M 函数 和?(z)的形式。然后根据位移的单值条件和内边界面力边界条件确定待定的系数。最后得到多连域中位移和应力分量单值连续的K-M 函数形式。 学习要点

1、单连域中的单值解析函数 和??(z)在多连域中可能是多值的；

13

2、应力分量的单值条件；3、位移分量的单值条件； 4、多连域中位移和应力分量单值连续的K-M 函数形式：

1、单连域中的单值解析函数 和??(z)在多连域中可能是多值的

对于弹性力学的应力解法，若K-M 函数 和??(z)满足公式，即

则应力分量已经满足平衡微分方程，变形协调方程，对于单连域问题， 和??(z)均为单值解析函数。根据边界条件，问题就可以求解。

但是对于多连域问题中，K-M函数 和?(z)可能表现为多值函数，尽管它们在单连域中是单值连续的解析函数。

对于多连域弹性体S，具有m个内边界 和一个外边界 ，而 分别为内边界 中的点，如图所示

那么如何选择这些K-M 函数，才能保证应力分量和位移分量的单值连续条件呢。这里的原则是保证应力和位移分量的单值性，分别根据应力分量的单值条件，构造可能的K-M 函数 和?(z)的形式，然后根据位移的单值条件和面力边界条件确定待定的系数。

14

2、应力分量的单值条件

由于应力分量必须是单值的，而应力分量与 K-M 函数 的关系，有

所以 的实部，即Re 必须是单值的。

假如 函数环绕多连域内部任意一个内边界lk 绕行一周时， 如果多值，只能是虚部多值。根据应力表达式 其多值部分只能是一个虚常数增量。为方便进一步分析，令此虚数增量为2?iAk，其中Ak为实常数。

根据复变函数性质，若复变函数绕lk 一周，如果有增量，其只能是对数函数产生的。因此设 由两部分组成，一部分是在S内单值解析的 ；另一部分是Ak ln (z-zk)，则其绕lk 一周有增量2?iAk。有

当 绕 lk 一周时，除了Akln (z-zk)以外，其余各项均恢复原值。其中zk 为lk 内任一点，即其在域S之外。对上式积分可得

应该注意的是， 在多连域内是单值连续的，但是其积分却不一定是单值连续的。设其有增量2?iCk，则

将上式代入复位势函数表达式，可得

上式中，Ak为实常数，而?k为复常数。即在多连域内， 为一个单值解析函数再加上前面两项。对于应力分量表达式

由于 ，而zk在域S之外， 域内为单值解

15

析函数。因此?'(z)也必须为单值解析函数。但是?(z)不一定是单值解析函数，作分析同前，有

其中 为单值解析函数，?k'为任意复常数。由此，对于多连域，K-M 函数 和?(z)的确定出现了三个待定常数Ak，?k和?k'。其值必须由位移单值条件和面力边界条件确定。 3、位移分量的单值条件

对于平面应力问题，位移分量为

当z绕lk一周时，则上式成为

因此，位移单值条件要求

通过位移单值条件，只有一个复常数还不能确定。 4、多连域中位移和应力分量单值连续的K-M 函数形式

位移单值条件没能确定的另一个复常数条件将根据面力边界条件确定。对于内边界lk，设边界面力的主矢量为

则 将公式

16

代入上式，则

联立求解，可得

将上述待定常数回代公式

则

上述公式为多连域内保证位移和应力单值连续条件的 和?(z)的表达式，其中 和 为多连域区域内的单值解析函数。 §8.6 无限大多连域中??f(z)和?(z)的表达式 学习思路：

尽管K-M函数的基本形式已经确定，但是对于一般的弹性力学问题，仍然难以确定函数的具体形式。

本节讨论无限大多连域的K-M 函数表达形式。利用无穷远边界条件，简化对数函数形式，并且在内边界之外将K-M 函数的解析函数部分展开成劳伦级数。并且利用应力有界条件和无穷远应力确定部分级数系数，为进一步工作奠定基础。 学习要点：

1、无限大多连域中K-M 函数的一般形式；

2、利用应力分量有界简化K-M 函数；

17

3、无穷远应力与K-M函数

1、无限大多连域中K-M 函数的一般形式

对于无限大多连域，其外边界趋于无限远。因此可以借助其无穷远处的边界条件，写出K-M 函数的表达形式。

以坐标原点为圆心，作半径为R 的圆 ，将所有的内边界均包围在此圆之内部。那么对于 之外的任意点都有|z|＞| |，因此

因此，公式

可以表示为

其中 为所有m个内边界上的表面力在x和y方向的分量的代数和，而 和?**(z)为 以外区域内除了无穷远点的解析函数。在无穷远处， 和?**(z)可能为解析函数，也可能是非解析的。它们在 以外区域内可以展开成劳伦级数

18

2、利用应力分量有界简化K-M 函数

将K-M 函数的表达式

和

代入应力分量表达式

有

上式右边部分项 将随|z|的增加而趋于无限大，因此当??趋于无穷远时，为使应力分量不至于成为无穷大，必须有

同理，如果应力表达式 的应力分量有界，则 在无穷远处有界，所以

于是，为了使应力分量在无穷远处保持有界，则K-M 函数的形式为

公式中

19

上式中 和??0(z)在 以外区域，包括无穷远处均为解析函数。由公式

可知，如果令 ，将不会改变应力分量，因此

其中

3、无穷远应力与K-M函数

上述公式

中的常数B和B'+iC' 在无穷远处具有力学意义，说明如下。

因为，当 。所以在无穷远处，由公式

有

设?1，?2为弹性体无穷远处的主应力，如图所示

20

而?为?1与x轴的夹角，则

由此可得

所以

可见常数B与无穷远处的两个主应力之和成正比，而常数B'+iC' 与无穷远处的两个主应力之差成正比 §8.7 保角变换和曲线坐标 学习思路：

弹性力学问题的求解有赖于边界条件的简化。对于复杂的边界形状，如果利用空间的变换，将是简化问题求解的最好途径。保角变换就是充分发挥复变函数的特长，将孔口问题映射到??平面的单位圆。

这一节将介绍保角变换和曲线坐标的概念。由于应用保角变换，矢量－位移，

21

张量－应力公式以及K-M 函数等均必须做出曲线坐标描述。保角变换使得问题的公式复杂，但是边界条件的简化，以及柯西积分的应用将简化问题的分析。

在本节学习之前，请你先学习附录2，（有关保角变换的知识） 学习要点：

1、保角变换和曲线坐标；2、矢量的保角变换；3、位移分量的曲线坐标表达式；4、应力分量的曲线坐标表达式。 1、保角变换和曲线坐标

为了便于根据边界条件确定K-M函数，采取保角变换

z = ??(?)

将物体在 z 平面上所占的区域变为在? 平面所占的区域。一般的说，通过保角变换可以将非圆边界映射为圆边界，使得问题得以简化。

假设将z平面上的有限区域或者无限区域S 映射为? 平面的单位圆内的区域

?，并且将z平面上的区域S的边界 l 映射为单位圆?，对应的关系如下表：

? 平面 ? =0（无穷远点） ? =const（圆） ??=const （半射线） 域? d? z平面 z=0（原点） ? =const （曲线） ??=const （曲线） 域S dz 由于??平面上的任一点可以表示为， 。? 和? 是点??的极坐标。

而根据保角变换公式z = ??(?)， 则z平面任意一点也可以通过?和? 表示。因此，??和??又称为曲线坐标。对于某些问题的描述中，采用曲线坐标形式表示位移和应力有利于问题的分析。

曲线坐标的概念：? 平面的一个圆周??=const和一条径向直线??=const分别对应于z 平面的两条曲线，这两条曲线就记作??=const和??=const。于是? 和?可以看作z平面上一点的曲线坐标。由于变换的保角性，这个曲线坐标总是正交的，而且坐标轴??和??的相对位置和坐标轴Ox和Oy的相对位置相同，如图所示

22

2、矢量的保角变换

首先讨论矢量的保角变换。设曲线坐标? ，即??=const与x轴夹? 角，如果A为z平面上的任一矢量，设A与曲线坐标 ??夹? 角。设Ax, Ay分别表示矢量A在x，y 轴的投影；A??,A??表示在? =const和 ??=const上的投影，则

上式的几何意义为，将矢量A绕z点顺时针方向转动? 角后，其在Oxy坐标系的位置，相当于A在曲线坐标系（?，?）中的位置，如图所示

如果用u? , u??分别表示曲线坐标下的位移矢量分量，则

根据保角变换，有

23

所以

沿曲线（?）取微分线段dz，则在? 平面对应的有d?，由于

所以，取其共轭可得

将上式回代到公式

可得

3、位移分量的曲线坐标表达式

下面通过保角变换对弹性力学的公式作对应的转换。首先，设K-M函数 和??(z)分别使用 和??1(z)代替，同时令

根据位移表达式

有

24

在 z 平面上，将位移矢量向曲线坐标? 和? 投影。由公式

可得

上式两边同时乘以2G，可得

上式是? 平面上的曲线坐标系表达的位移表达式。 4、应力分量的曲线坐标表达式

下面建立曲线坐标中应力分量的复变函数表达式。如果用??, ?? , ??? 表示物体在曲线坐标中的应力分量。则

因为

和

而由公式

所以

25

上式为经过保角变换后，z平面上的曲线坐标系中的应力分量的复变函数表达式。 §8.8 无限大薄板的孔口问题 学习思路：

本节的主要任务是将保角变换用于无限大薄板的孔口问题，确定K-M 函数的基本求解公式。

推导中首先确定无限大板孔口问题的保角变换公式，将K-M 函数转换为曲线坐标形式。采用的方法仍然是将K-M 函数分解为以级数表达的解析函数和对数表达的多值函数两部份。

对于K-M 函数的级数形式，通过孔口面力边界条件可以确定级数函数的求解方程。这个求解过程，利用保角变换后孔口边界的特殊性质，使用柯西积分使得计算简化。 学习要点：

1、保角变换公式与K-M 函数；2、利用孔口边界条件确定K-M 函数求解公式；3、柯西积分确定K-M 函数的级数形式。 1、保角变换公式与K-M 函数

保角变换的目标是：将z平面上的孔口边界l映射为??平面上的单位圆?，将l 以外的无限区域S 映射为??平面上的单位圆内的有限区域?，将z平面上的无穷远点映射为? 平面的坐标原点，如图所示

保角变换公式：

是将l以外的无限区域映射为单位圆??内（|?|＜1＝的普遍变换式，公式中R为实

26

数，Ck为复数，而且 <1。

保角变换公式确定以后，可以确定K-M函数 和??(?)，即将K-M 函数 和?1(z)变换到曲线坐标中去。

其中

因为

由于 <1 ，将上式展开，有

所以

ln z = ln ? +单位圆内部?的解析函数。

另外

根据上述分析， 的各项都转变为单位圆内??的单值解析函数。因此

其中

27

2、利用孔口边界条件确定K-M 函数求解公式

讨论边界条件确定K-M 函数 和??0(?)。根据面力边界条件

，经过保角变换后，可得

在单位圆的圆周上

所以上述面力边界条件可以表示为

根据公式

则在边界即单位圆周上

将上述K-M 函数的边界值回代面力边界条件，并且将已知函数与需要确定的未知函数分开，可得

其中已知函数为

28

3、柯西积分确定K-M 函数的级数形式

因为 和? 0(? )是单位圆内的泰勒级数，它们是从z平面上lR之外无穷区域的罗伦级数转化而来的。因此对于公式

幂级数求解时，由于方程两边都含有? k=e ik? 的各个项（k由－∞到∞），比较各个同类项的系数，即可求得ak，bk 的值。不过这样作太麻烦了，由于 和? 0( ??)在单位圆内是解析的，而且在圆内和圆周上是连续的，因此可以直接采用柯西积分计算。

将边界条件的第一式两边乘以 ，积分可得

由于 在单位圆内是解析的，因此公式的第一个积分即等于 ，它是级数 之和。对于公式第三项的积分函数，由于

在单位圆外是解析的，在圆外和圆周上是连续的，所以

。因此，边界条件的第一式就成为

同理，边界条件的第二式成为

上述公式就是边界条件通过柯西积分所推导出的计算 和? 0(? )表达式。

29

其中 是边界的已知函数。 §8.9 带椭圆孔无限大薄板受均匀拉伸 学习思路：

为了进一步说明无限大板的孔口问题的求解，本节以椭圆孔口问题为例作详细的讨论。

主要工作和推导过程为：首先建立椭圆孔的保角变换公式；根据椭圆孔保角变换公式确定 和?0(?)；然后分别确定 和?(?)。最后计算孔口应力，分析孔口应力集中。 学习要点：

1、椭圆孔口的保角变换；2、确定 和?0(?)的求解公式；3、确定保角变换的 和?(?)函数；4、确定K-M函数表达式；5、孔口应力与应力集中 1、椭圆孔口的保角变换

为了进一步说明无限大板孔口问题的求解，本节以椭圆孔口问题为例，作详细的讨论。

将平面上的椭圆孔的外部区域变换到??平面的单位圆区域内的保角变换公式为

对于无限大板孔口问题，R和m均为实常数，而且R＞0，0≤m≤1。上式也可以写作

所以

从上式中消去?，可得

30

从上式中消去?，可得

根据上述分析，? 平面上的圆周? =const对应于z平面上的中心椭圆，椭圆的长短半轴分别为 如图所示

? 平面上的径向线? =const 对应于z平面上的中心双曲线。这一族椭圆和一族双曲线就是z平面上的曲线坐标。特别应该指出的是，? 平面上的圆周? =1的圆周所对应的z平面上的椭圆。

设内边界（椭圆孔）的方程为

由于保角变换是将椭圆孔口以外区域映射于? 平面的单位圆以内区域，令公 式 中，? =1（? 平面的单位圆），与椭圆方

程比较可得

2、确定 和??0(?)的求解公式

根据保角变换公式 ，可以得到

31

将上述结果代入公式的第一式

则

将

代入上式，可以得到

由于 在单位圆??以外是解析的，在单位圆??的边界上是连续的，所以其柯西积分为零。因此

将公式

代入公式的第二式

32

则

由于 在单位圆??以内是解析的，在单位圆??的边界上是连续的，所以

3、确定保角变换的 和?(?)函数

对于单向拉伸的椭圆孔口板，如图所示

??1=p，??2=0，Fxk= Fyk =0，Fx= Fy =0，则

33

将上述结果代入公式

利用柯西积分关系式，则

将上式求导后代入公式

并且利用柯西积分关系式，则

整理并且将上述结果代入公式

则

4、确定K-M函数表达式

以下计算应力分量，由于

因为

34

所以

根据公式

有

因此,根据公式

可以求得

5、孔口应力与应力集中

将上述计算结果代入公式

35

简化可以得到应力分量的曲线坐标表达式

如果将??=?（cos???+isin?）代入上述二式，分离实部和虚部，则可以求出应力分量??，???和 ????。但是这个工作太复杂了，由于孔口应力是薄板最为重要的应力。如果仅考虑孔口应力，可以避开上述冗长的推导过程。

由于在孔口边界上,??=ei? ,??=????。所以

当拉力p平行于x轴，此时??=0，由上式可得孔口应力为

最大应力为

最小应力为

当拉力p垂直于x轴

36

此时??=?/2，孔口应力为

最大应力为 最小应力为

假如a=b，即讨论问题为圆孔，则 。

§8.10 裂纹前缘的应力分布 学习思路：

本节应用椭圆孔口分析结果探讨裂纹应力分布。注意到裂纹是短轴为零的椭圆，因此可以应用椭圆孔口结论。这里的问题是裂纹前缘局部前缘应力不能采用曲线坐标描述，需要将问题重新转换到z平面。即将??平面的结论映射回z平面。并且以裂纹尖端为新的极坐标原点建立坐标系分析裂纹应力。

在上述裂纹前缘应力计算公式中，如果命??趋近于零，则各个应力分量的数值将趋于无限大。这就表示，在裂纹前缘，应力是无限大的。实际上，由于裂纹前缘总是有或大或小的塑性区，因此就不会发生无限大的应力，上述公式仅适用于弹性范围的应力分析。尽管如此，但对于脆性材料，塑性范围很小的情况，公式可以令人满意的描述裂纹前缘的应力状态。因此，以上公式成为断裂力学的基本公式。

37

学习要点：

1、裂纹—短轴为零的椭圆；2、将??平面确定的K-M 函数转换到z平面；3、裂纹前缘应力分布；4、切应力作用的裂纹前缘应力 1、裂纹—短轴为零的椭圆

对于椭圆孔口问题，如果短半轴b=0，则椭圆孔退还为一条长为2a的裂纹，如图所示

此时R=a/2，m=1，保角变换公式

为

根据公式

有

38

若拉力p与裂纹垂直，即??=?/2，则

2、将??平面确定的K-M 函数转换到z平面

对于裂纹问题，主要是裂纹前缘局部区域应力分析。因此问题需要在裂纹板平面讨论。所以将问题变换到z平面。将?用z表示，则由公式

可得

因为当z平面上的距孔口无限远的点，对应于? 平面上的单位圆的圆心。即|z|趋于无限大时，|?|应趋于零，所以上式根号前取负号。将上式代入公式

则

代入公式

即可以求得应力分量的表达式

39

3、裂纹前缘应力分布

在裂纹前缘区域，重要的是分析裂纹尖端附近的应力分布。因此，选取以裂纹尖端点为原点的极坐标系。则

z=x+ iy=a+??ei? =a+??cos?+ i??sin??

所以

x=a+??cos??, y= ??sin?

将以上结论代入公式

即可求出用极坐标表示的应力分量?x，?y，?xy。当然，这个应力分量表达式将是非常冗长的。对于裂纹问题，最为关注的是裂纹前缘局部区域的应力状态。在裂纹前缘区域，? 远小于a。因此可以将表达式按?/a的升幂次展开，并且略去高阶小量。则可得到

从而求出应力分量

4、切应力作用的裂纹前缘应力

当薄板在裂纹方向及其垂直方向都受有均匀分布的剪力p时。其受力显然可以用下述情况来代替：在??=?/4的方向受均匀分布的拉力p，并且在??=-?/4的

40

本文来源：https://www.bwwdw.com/article/mn47.html