弹塑性力学第九章课件

第九章空间轴对称问题

本章讨论空间轴对称问题的基本方程和一些轴对称问题的基本解。对于一般空间问题的解法我们在第五章已有讨论，但一般空间问题一般解（具体求解）通解讨论在杜庆华等编著的“弹性理论”中有较多的论述。我们不刻意从数学上论述一般空间问题一般解的表达式，而对于空间轴对称问题作一些讨论和举例。3/27/20231

1.1空间轴对称问题特点：1.域内所有物理量（体力、面力、位移、应力、应变）均为r、z的函数。与平面轴对称问题类似，空间轴对称问题的求解域、荷载和约束绕某一轴（z轴）对称，导致如下简化，2．荷载：体力f=0，面力

，位移u=0，应力r=z=0，应变r=z=0。第一节空间轴对称问题的基本方程3/27/20232第一节空间轴对称问题的基本方程3．待求的物理量（10个）：ur、w、r、、z、rz=zr、r、、z、rz=zr1.2基本方程1.平衡微分方程（两个）：3/27/202334.物理方程（四个）：第一节空间轴对称问题的基本方程3/27/20235

r=e2Gr、=e2G、

z=e2Gz、rz=Grz

第一节空间轴对称问题的基本方程其中——体积应变或

3/27/202365.边界条件第一节空间轴对称问题的基本方程位移边界：在Su上6.按应力解法

力的边界：在r=r0

在z=z0

四个应力分量r、、z、rz为基本未知量。3/27/20237第一节空间轴对称问题的基本方程(r,z)满足第一个平衡微分方程，而第二个平衡方程及四个相容方程，共同要求

22=4=0

——(r,z)应满足的基本微分方程。3/27/20239

7．按位移法解

第一节空间轴对称问题的基本方程其中

a．基本未知函数：ur和w基本方程两个：

并考虑适当的边界条件。3/27/202310b.

引入Love（拉甫、勒夫）位移函数（当无体力作用时）第一节空间轴对称问题的基本方程对于位移法的基本方程的解可由考虑体力的一个特解加上齐次方程的通解。轴对称问题齐次拉梅方程的通解可以引入一个Love位移函数(r,z),使得位移由(r,z)表示：3/27/202311轴对称问题按位移求解，归结为寻找一个恰当的重调和函数(r,z)，使按其导出位移和应力能满足给定的边界条件。第一节空间轴对称问题的基本方程比较应力函数解法和love位移法知：(r,z)=

(r,z)3/27/202313第二节半空间体在边界上受法向集中力

（Boussinesq问题）半空间体，体力不计，边界受法向集中力P作用.轴对称问题，P作用在坐标原点上。zRrPx

yz已知，当z=0且r0时,z=0,zr=0；当R0时，应力奇异。当R

时，R=(r2+z2)1/2，

应力、位移

0；3/27/202314选

(r,z)

为r和z的正一次幂式:(r,z)=A1R+A2[R-zln(R+z)]——为双调和函数第二节半空间体在边界上受法向集中力

（Boussinesq问题）Boussinesq采取Love函数求解，(r,z)为重调和函数，由(r,z)的三次微分导出应力。zRrPx

yz3/27/202315应力：

第二节半空间体在边界上受法向集中力

（Boussinesq问题）3/27/202317根据边界条件来确定A1和A2：第二节半空间体在边界上受法向集中力

（Boussinesq问题）zRrPx

yz在z=0且r0边界上,z=0自然满足。在z=0且r0边界上,zr=0（1-2）A1+A2=0—(a)3/27/202318在z=z0

0平面上，要求z的合力与P平衡。还需一个条件（包括P的）。第二节半空间体在边界上受法向集中力

（Boussinesq问题）将z表达式代入，得zPrrdrz0z3/27/202319由式(a)、(b)解得

A1=P/(2)、A2=-（1-2）P/（2）第二节半空间体在边界上受法向集中力

（Boussinesq问题）

代回位移、应力表达式，见徐芝纶（上册）P.297（9－17）、（9－18）式，称为Boussinesq问题解。由P.297（9－17）、（9－18）式见：位移和应力随R的增加而减小。3/27/202321Prz第二节半空间体在边界上受法向集中力

（Boussinesq问题）在z=0平面上3/27/202322第三节半空间体在边界上受法向分布力q已知条件：半空间体在边界上受均布法向荷载q作用，在半径为a的圆面积。zaqar寻求解答：1.

z=0边界上的沉陷wz=0

=?2.r=0（对称轴）上的应力和位移。求解方法：采用叠加法和半空间体边界受法向集中力P的计算结果求解。3/27/202323圆面积均布荷载q对圆外M点竖向位移影响可取一个微面元，距M点为s，角度为处，dA=sdds，dA上q对M点影响：

第三节半空间体在边界上受法向分布力qrraMs1s2sdsdzaqar3/27/202325rraMs1s2sdsd第三节半空间体在边界上受法向分布力q3/27/202326由图形可见

asin=rsin,两边微分

acosd=rcosd第三节半空间体在边界上受法向分布力qrraMs1s2sdsd3/27/202329第三节半空间体在边界上受法向分布力q的取值范围：由0

1

rraMs1s2sdsd的取值范围：0

3/27/202330第三节半空间体在边界上受法向分布力q3/27/202331第二类椭圆积分

第一类椭圆积分第三节半空间体在边界上受法向分布力q对于不同a/r可由椭圆积分表得到。3/27/2023322．M点载荷在圆之内：Masdsdrmn第三节半空间体在边界上受法向分布力q圆内距M点s处微面积q对M点沉陷的影响仍为3/27/202333整个圆面积荷载引起M点沉陷为：第三节半空间体在边界上受法向分布力q第二类椭圆积分利用asin=rsin

3/27/202334当r=0为圆心处沉陷：当r=a时圆周上沉陷：

3.2在z轴r=0上的应力和位移在z轴上的应力和位移比同一水平面上其它点的应力和位移要大。第三节半空间体在边界上受法向分布力q3/27/2023351.应力：由于z轴对称轴，所以在z轴上的应力无剪应力，均为主应力：

r=、z第三节半空间体在边界上受法向分布力q3/27/2023362．位移：z轴上的ur=0，仅存在w第三节半空间体在边界上受法向分布力q3/27/202337第三节半空间体在边界上受法向分布力q3/27/202338第四节两球体之间的接触压力接触压力问题是在机械工程、土木工程中经常碰到的问题，接触问题在1881年由德国赫兹（HeinrichHerty）首先用数学弹性力学导出了计算公式。4.1接触问题的特点：

1．两个弹性体互相接触，当无压力作用时，为点接触或线接触。当有压力作用时，弹性体发生变形，点接触（或线接触）变为面接触。3/27/2023392．弹性体变形后的接触面为非常小的局部区域（相对于弹性体几何尺寸）所以可看成半空间（半无限平面）体法向受局部分布力作用问题，但这里分布力q不是均匀的，同时q也未知，接触面的局部区域也是未知的。第四节两球体之间的接触压力3．不计接触面摩擦力。

3/27/202340

4.2

两球体之间的接触压力：已知两球体变形前在o点接触，两个坐标系

roz1、roz2第四节两球体之间的接触压力rOz1z2O2O1R2R1球1：E1

、1、R1球2：E2

、2、R2

M1M2r距接触点z轴为r的两球表面上M1和

M2点的z坐标分别为（M1和M2与点o很近）3/27/202341第四节两球体之间的接触压力rOz1z2O2O1R2R1M1M2r则3/27/202342第四节两球体之间的接触压力在已知P压力作用下，两球在接触点附近发生变形有一个接触面，根据对称性接触面为以a为半径的圆。rOz1z2O2O1R2R1M1M2rM1rPPoz1z2O1M2ar3/27/202343第四节两球体之间的接触压力1．a为待求量，同时接触面上有接触压力q（待求）。2．由于接触问题是局部变形，在球体远离o点的任意点位移为刚体位移。两球内距o点很远处的相对位移（刚体位移）为

？下面要建立（找出）三个条件（几何、物理、平衡方程）寻求a

、q

和。3/27/202344第四节两球体之间的接触压力求解：首先根据接触面变形（位移）来建立一个关系球1：触面上o点、M1点沿z1轴位移为w1(o)、w1而w1(o)=w1+z1

M1rPPoz1z2O1M2ar3/27/202345第四节两球体之间的接触压力球2：触面上o点、M2点沿z2轴位移为w2(o)、w2w2(o)=w2+z2

而w1(o)+w2(o)=w1+z1+w2+z2w1(o)+w2(o)=w1+w2+r2或M1rPPoz1z2O1M2ar3/27/202346而w1(o)+w2(o)=第四节两球体之间的接触压力——两球体距o点较远处两点的趋近距离。

=w1+w2+r2——变性协调关系w1(o)+w2(o)=w1+w2+r2由于接触问题可看成半无限体受局部垂直分布力问题，w1和w2可以利用上一节的结果。M1rPPoz1z2O1M2ar3/27/202347第四节两球体之间的接触压力—相当物理和几何关系3/27/202348代入

=w1+w2+r2在此式中a、q和未知。第四节两球体