解析函数的应用—浅谈在陌生弹性力学中的应用（杜碧晶，运城学院数学系）摘要：在数学中，我们知道一个复变函数如果解析，则其实部和虚部均为调和函数，满足调和方程。

一个实变的双调和函数，可由共轭复变函数的线形组合得到。

在平面弹性力学中，对于平面应力问题和平面应变问题，可以通过假设，转变成求解满足某些边界条件下的双调和方程问题。

这样就可以用复变函数中的解析函数进行解决。

关键词：解析函数、应力函数、平面应力问题、平面应变问题。

1、引言：社会十分尊重数学，这可能不是因为这个学科的内在美，而是因为数学是社会极其需要和工程中有广泛应用的一种艺术。

以复数作为自变量的函数叫做复变函数，而与之相关的理论就是复变函数论。

解析函数是复变函数中具有解析（可微）性质的函数。

如果一个复变函数解析，那么它的实部和虚部均为调和函数，满足拉普拉斯调和方程(02222=∂∂+∂∂yx φφ)。

在区域D 内满足C —R 方程即：xvy u y v x u ∂∂-=∂∂∂∂=∂∂,的两个调和函数v u ,中，v 称为u 在区域D 内的共轭调和函数。

任何一个弹性体都是空间问题，一般的外力都是空间系，因此严格的说，任何一个实际的弹性力学问题，都是空间问题。

但是所考察的弹性体具有某种特殊形状，并且承受的是某种特殊的外力，就可把空间问题简化为平面问题。

这样处理后，分析和计算的工作量将大大的减少，而所得的结果仍满足工程上对精度的要求，因此具有广泛的实用价值。

弹性力学的平面问题分为平面应力问题和平面应变问题。

其中平面应力问题是指很薄的等厚度薄板只在板边上受有平行平板面并且不沿厚度变化的面力，同时体力也平行于板面并且不沿厚度变化；平面应变问题是指很长的柱形体，它的横截面在柱面上受有平行于横截面而且不沿长度变化的面力，同时体力也平行于横截面而且也不沿长度变化，即内在因素和外来作用都不沿长度而变。

2基础内容介绍如果采用应力作为基本未知量求解弹性力学平面问题，在常体力的条件下基本方程归结为在给定的边界条件下求平衡微分方程,和应力表示的变形协调方程对于平衡微分方程的解，可以分解为其齐次方程的通解与任一特解之和。

齐次方程就是体力为零的平衡微分方程，显然，平衡微分方程的特解是容易寻找的，下列应力分量均为齐次方程的特解，根据微分方程理论，必有函数f（x，y），令则齐次方程的第一式恒满足。

同理必有函数g（x，y），如果, 则齐次方程的第二式恒满足，所以引入任意函数(x，y)，使得将上式分别回代，可得应力分量表达式上述应力分量即为齐次平衡微分方程的通解。

对于体力为零的弹性力学平面问题，只要函数是四阶连续可导的，总是满足齐次微分方程的。

将平衡微分方程特解代入应力表达式，则自然满足平衡微分方程。

应力分量不仅需要满足平衡微分方程，而且还需要满足变形协调方程，将上述应力分量代入变形协调方程，可得上式说明函数(x，y)应满足双调和方程。

根据应力函数计算的应力分量满足平衡微分方程，而双调和函数表达的应力函数自然满足变形协调方程。

因此双调和方程就成为平面问题应力解法的基本方程。

因此，研究弹性力学的平面问题，当体力为常量时，可归结为求解满足双调和方程022=∇∇φ的应力函数φ，并且使其在边界上满足所以的边界条件。

在实变函数中，难以用应力函数表示位移分量和位移边界条 件。

同时若将应力边界条件和位移边界条件相比，我们发现它们除)(z ψ项的系 数不同外，其未知函数的表达式是相同的。

这说明平面问题的两种边界条件问题，从数学观点来看，其求解方法基本相同。

也可以说，复变函数解法把求应力和位移的方法统一起来了。

就是说可以不必区分按应力求解或按应变求解，而是把两者统一起来，这是复变函数解法的方便之处。

在弹性力学的平面问题中，孔口问题最能显示复变函数解法的优越性。

下面我们一起来分析看复变中的解析函数如何将弹性力学的平面问题解决的最容易。

3正文：在弹性力学的平面问题中，当体力为常量时，应力函数，),(y x φ为双调和函数，即应满足0),(22=∇∇y x φ。

为了将双调和方程转化为用复变函数表示，需进行变量代换。

因为：,iy x z += iy x z -=___所以： ),(21___z z x += )(21___z z y --=；这样，),(y x φ可以通过中间变量z 和___z 看成一个复合函数，因此有：φφφφφφ)(_________zzz z x z z x z z x ∂∂+∂∂=∂∂+∂∂=∂∂⋅∂∂+∂∂⋅∂∂=∂∂； φφφφφφ)(________zzi i z i z y z z y z z y ∂∂-∂∂=∂∂-∂∂=∂∂⋅∂∂+∂∂⋅∂∂=∂∂； 从而： ,__zzx ∂∂+∂∂=∂∂ );(___z z i y ∂∂-∂∂=∂∂由此可得： ,2)(__zyi x ∂∂=∂∂+∂∂ ;2)(zy i x ∂∂=∂∂-∂∂ （1）由（1）式可得：.4___22zz ∂∂∂=∇于是双调和方程022=∇∇φ变换为0162__24=∂∂∂zz φ积分四次得出通解为：);()()()(__43____21),(z zF z F z z F z F y x +++=φ（2）； 因为),(y x φ是实函数，所以（2）式应为两两共轭函数，即：,)()(______1___2z F z F = ;)()(______3__4z F z F =于是用复变函数表示双调和函数的一般形式为： ______33_______11),()()()()(z F z z F z z F z F y x +++=φ （3）；将（3）式中的)(1z F 及)(3z F 分别用)(21z χ及)(21z ψ代替，即得有名的古莎公式：][21____)()(______)()(),(z z z z y x z z ψψχχφ+++= 或)(__)(),(}Re[z z y x z ψχφ+= （4） （4）式就是用复变函数表达的应力函数。

式中的)(z χ和)(z ψ都是z 的解析函数。

于是求解弹性力学的平面问题就成为寻找满足边界条件的两个复 变函数)()(,z z ψχ的工作了。

弹性力学在的平面问题中，孔口问题最能显示复变函数解决问题的优越性。

例：图为一具有小圆孔（半径为a ）的无限大平板，在无限远处沿轴向承受均匀的拉力，求其应力分量和位移分量？解：作用在弹性体上的外力，称为面力。

由于孔口不受面力。

所以在孔口边界处有：Fy Fx =0= Ry Rx =0=上式中的Fx、Fy 代表弹性体在X ，y 轴上的面力分量，Rx ，Ry 代表m 个内边界分别在X ，y 方向上的面力之和。

故应力函数的表达式（1）为： )(z ψ =nn n z a-∞=∑1+AZ-πμ81+（R X +iR Y ）LnZ)(/Z χ=n n n z b -∞=∑1+（B+)/iB Z+(83πμ-R X -iR Y ）LnZ变形为（2）式：)(z ψ=n n n z a -∞=∑1+AZ )(/z χ =nn n z b -∞=∑1+（B+/iB ）Z 上式中的一次幂A 和（)/iB B +具有力学意义，它们代表无限远处的应力状态。

（在某一点处某一截面上的应力是指该截面上的附加分布内力在该点处的集度。

）作用在弹性体内部各点上的外力，称为体积力，简称体力。

在不计体力时，应力分量和应力函数之间的关系为：x σ=22y∂∂φx σ=22x ∂∂φ xy τ=y x ∂∂∂-φ2应力表达式： ⎥⎦⎤⎢⎣⎡+=+______/)(/)(2z z y x ψψσσ ⎥⎦⎤⎢⎣⎡+=+-)(//__)(/22z z xy x y z i ψχτσσ由应力函数的表达式求出的)(/z χ、)(//__Z Z ψ和)(//z χ代入应力表达式。

因为无限远处为均匀应力场，应力分量记为∞∞∞xy y x τσσ,,当z ∞→时，除一次幂对应的应力存在外，其余项均消失。

所以,4A y x =+∞∞σσ ()/22iB B i xy x y +=+-∞∞∞τσσ则有：A=)(41∞∞+x y σσ B=)(21∞∞-x y σσ ∞=xyB τ/计算解得：A=σ41B+σ21/=iB 。

（3）将（3）式代入应力函数表达式（2），则有（4）式：n i n n z z a z -∞=∑+=4)(σψn n n z z b z -∞=∑+-=1)(/2σχ积分形式的复应力边界条件：dS iF F i z y x z z z )(______)(/)(_____)(/+=++⎰ψψχ将0==y xF F 代入可得：0______)(/)(_______)(/=++z z z z ψψχ其共轭式为：0)(/________)()(/=++z z z z ψψχ（其中Z=θi ae ）其中a 为小圆半径，将上面由（4）式所求得的)(z ψ、)(/z χ代入复应力边界条件的共轭式可得：0224)(4)(2)2(111)1(11=-+++-=-++++-+-∞=∞=-∞=-+-------∞=∑∑∑∑∑∑θθθθθθθθθθθθσσσσσn i n n n in n n n in n n n i i n i n n i i ni n n i n i n n i e ana e a a e a b ae ae e a a ae ae e a a ae ae b ae比较等式两边e+-θin （n=1、2、、、、、∞）项的系数：θi e 项的系数为：021=+-a a a σ所以：212a a σ= θin e项的系数为：)2(0≥=n aa n n所以：)2(0≥=n a nθi e -项的系数为：021=+a b a σ所以：212a b σ-= θ2i e-项的系数为：022=ab 所以：02=bθ3i e -项的系数为：0133=-a a ab 所以：41232.a a a b σ== θin e-项的系数为：)4(0)2(22≥=----n aa n ab n n n n 所以：)4(0≥=n b n将上面所求得的系数代入（4）式，使的复应力函数成为（5）式：)2(422)(za z Z +=σψ)1(24422)(/za z a z z -+-=σχ因为（6）式：)(/_______)(/)(/__22Re 4][24z z z x y zz ψψψφφσσ=+=∂∂∂=∇=+φφφφτσσ222222)()(22y i x y x i yx xy i x y ∂∂-∂∂=∂∂∂-+∂∂-∂∂=+-； 又因为：zy i x ∂∂=∂∂-∂∂2 所以：（2224)z y i x ∂∂=∂∂-∂∂ 所以：][242)(//__)(//22z z xyx yz zi ψχφτσσ+=∂∂=+-（7）利用应力分量的坐标变换式：θτθσσσσσ2sin 2cos 22xy yx xy r+-++=；θτθσσσσσθ2sin 2cos 22xy y x x y ---+= ； θτθσστ2cos 2sin 2xy yx xy +--= .将已得到的直角坐标中的复应力公式转化为极坐标的形式，则有（8）式：y x r σσσσθ+=+ θθθτσστσσi xy x y r r e i i )2(2+-=+-将（6）、（7）代入（8）式有： )(/Re 4z rψσσθ=+ θθθψχτσσ2)(//__)(//][22i z z r r e z i +=+-将（4）式代入（8）式，则有（9）式：);2cos 21()24Re(4)24Re(4)24Re(4)]2(4Re[42222222/2/22θσσσσσσσσσσθθra e r a z a z a z za z i r -=-=-=+=+=+-θθθθθσσσσσψχτσσ224232__44222)(//__)(//1()2322(2][22i i i z z r r er ae z a z z a z a e z i +-=+-+-=+=+-);32()23244222222332444θθθθθθθσi i i i i i i e rae r a e r a e e r a re e r a ----+-=+- 将（9）式中的第二式的实部和虚部分开，并和（9）式中的第一式联立求解。