弹塑性力学第九章

弹塑性力学第九章2023/8/141第1页，课件共57页，创作于2023年2月

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

=0，面力

，位移u

=0，应力

r

=

z

=0，应变

r

=

z

=0。第一节空间轴对称问题的基本方程2023/8/142第2页，课件共57页，创作于2023年2月第一节空间轴对称问题的基本方程3．待求的物理量（10个）：ur、w、

r、

、

z、

rz=

zr、

r、

、

z、

rz=

zr1.2基本方程1.平衡微分方程（两个）：2023/8/143第3页，课件共57页，创作于2023年2月2.几何方程（四个）：第一节空间轴对称问题的基本方程3.变形协调方程（四个）2023/8/144第4页，课件共57页，创作于2023年2月4.物理方程（四个）：第一节空间轴对称问题的基本方程2023/8/145第5页，课件共57页，创作于2023年2月

r=

e

2G

r、

=

e

2G

、

z=

e

2G

z、

rz=G

rz

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

2023/8/146第6页，课件共57页，创作于2023年2月5.边界条件第一节空间轴对称问题的基本方程位移边界：在Su上6.按应力解法

力的边界：在r=r0

在z=z0

四个应力分量

r、

、

z、

rz为基本未知量。2023/8/147第7页，课件共57页，创作于2023年2月基本方程（六个）：两个平衡微分方程与四个用应力表示的变形协调方程；再加上力的边界条件。第一节空间轴对称问题的基本方程如果体力为零时，基本方程为齐次方程，则可采用应力函数解法，引入应力函数

(r,z)，使得应力用

(r,z)表示:2023/8/148第8页，课件共57页，创作于2023年2月第一节空间轴对称问题的基本方程

(r,z)满足第一个平衡微分方程，而第二个平衡方程及四个相容方程，共同要求

2

2

=

4

=0

——

(r,z)应满足的基本微分方程。2023/8/149第9页，课件共57页，创作于2023年2月

7．按位移法解

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

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

并考虑适当的边界条件。2023/8/1410第10页，课件共57页，创作于2023年2月b.

引入Love（拉甫、勒夫）位移函数（当无体力作用时）第一节空间轴对称问题的基本方程对于位移法的基本方程的解可由考虑体力的一个特解加上齐次方程的通解。轴对称问题齐次拉梅方程的通解可以引入一个Love位移函数

(r,z),使得位移由

(r,z)表示：2023/8/1411第11页，课件共57页，创作于2023年2月代入齐次拉梅方程，第一式自然满足，而第二式为基本方程：

4

=0

(r,z)——为双调和方程。第一节空间轴对称问题的基本方程同时应力分量由

(r,z)表示为：2023/8/1412第12页，课件共57页，创作于2023年2月轴对称问题按位移求解，归结为寻找一个恰当的重调和函数

(r,z)，使按其导出位移和应力能满足给定的边界条件。第一节空间轴对称问题的基本方程比较应力函数解法和love位移法知：

(r,z)=

(r,z)2023/8/1413第13页，课件共57页，创作于2023年2月第二节半空间体在边界上受法向集中力

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

Px

yz已知，当z=0且r

0时,

z=0,

zr=0；当R

0时，应力奇异。当R

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

应力、位移

0；2023/8/1414第14页，课件共57页，创作于2023年2月选

(r,z)

为r和z的正一次幂式:

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

（Boussinesq问题）Boussinesq采取Love函数求解，

(r,z)为重调和函数，由

(r,z)的三次微分导出应力。zRr

Px

yz2023/8/1415第15页，课件共57页，创作于2023年2月

(r,z)=A1R+A2[R-zln(R+z)]则

(r,z)自然满足

4

=0。代入位移、应力计算式.第二节半空间体在边界上受法向集中力

（Boussinesq问题）zRr

Px

yz位移：2023/8/1416第16页，课件共57页，创作于2023年2月应力：

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

（Boussinesq问题）2023/8/1417第17页，课件共57页，创作于2023年2月根据边界条件来确定A1和A2：第二节半空间体在边界上受法向集中力

（Boussinesq问题）zRr

Px

yz在z=0且r

0边界上,

z=0自然满足。在z=0且r

0边界上,

zr=0

（1-2

）A1+A2=0—(a)2023/8/1418第18页，课件共57页，创作于2023年2月在z=z0

0平面上，要求

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

（Boussinesq问题）将

z表达式代入，得zPrrdrz0

z2023/8/1419第19页，课件共57页，创作于2023年2月P-4

A1(1-

）-2

A2=0——(b)第二节半空间体在边界上受法向集中力

（Boussinesq问题）而2023/8/1420第20页，课件共57页，创作于2023年2月由式(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的增加而减小。2023/8/1421第21页，课件共57页，创作于2023年2月Prz第二节半空间体在边界上受法向集中力

（Boussinesq问题）在z=0平面上2023/8/1422第22页，课件共57页，创作于2023年2月第三节半空间体在边界上受法向分布力q已知条件：半空间体在边界上受均布法向荷载q作用，在半径为a的圆面积。zaqar寻求解答：1.

z=0边界上的沉陷w

z=0

=?2.r=0（对称轴）上的应力和位移。求解方法：采用叠加法和半空间体边界受法向集中力P的计算结果求解。2023/8/1423第23页，课件共57页，创作于2023年2月3.1边界上一点M的竖向位移w：第三节半空间体在边界上受法向分布力q1.设M点为圆面积之外：M点可以在荷载圆面积之外也可在之内。zaqar当半空间体边界上受法向集中力P时，边界上距P点为r的点竖向位移为：2023/8/1424第24页，课件共57页，创作于2023年2月圆面积均布荷载q对圆外M点竖向位移影响可取一个微面元，距M点为s，角度为

处，dA=sd

ds，dA上q对M点影响：

第三节半空间体在边界上受法向分布力q

rraMs1s2sdsd

zaqar2023/8/1425第25页，课件共57页，创作于2023年2月

rraMs1s2sdsd

第三节半空间体在边界上受法向分布力q2023/8/1426第26页，课件共57页，创作于2023年2月整体圆面积荷载对M点影响为第三节半空间体在边界上受法向分布力q而

rraMs1s2sdsd

2023/8/1427第27页，课件共57页，创作于2023年2月

1为M点作为圆相切线OM线的夹角第三节半空间体在边界上受法向分布力q

rraMs1s2sdsd

为了简化积分将积分变量

转变为

2023/8/1428第28页，课件共57页，创作于2023年2月由图形可见

asin

=rsin

,两边微分

acos

d

=rcos

d

第三节半空间体在边界上受法向分布力q

rraMs1s2sdsd

2023/8/1429第29页，课件共57页，创作于2023年2月第三节半空间体在边界上受法向分布力q

的取值范围：由0

1

rraMs1s2sdsd

的取值范围：0

2023/8/1430第30页，课件共57页，创作于2023年2月第三节半空间体在边界上受法向分布力q2023/8/1431第31页，课件共57页，创作于2023年2月第二类椭圆积分

第一类椭圆积分第三节半空间体在边界上受法向分布力q对于不同a/r可由椭圆积分表得到。2023/8/1432第32页，课件共57页，创作于2023年2月2．M点载荷在圆之内：Ma

sdsd

rmn第三节半空间体在边界上受法向分布力q圆内距M点s处微面积q对M点沉陷的影响仍为2023/8/1433第33页，课件共57页，创作于2023年2月整个圆面积荷载引起M点沉陷为：第三节半空间体在边界上受法向分布力q第二类椭圆积分利用asin

=rsin

2023/8/1434第34页，课件共57页，创作于2023年2月当r=0为圆心处沉陷：当r=a时圆周上沉陷：

3.2在z轴r=0上的应力和位移在z轴上的应力和位移比同一水平面上其它点的应力和位移要大。第三节半空间体在边界上受法向分布力q2023/8/1435第35页，课件共57页，创作于2023年2月1.应力：由于z轴对称轴，所以在z轴上的应力无剪应力，均为主应力：

r=

、

z第三节半空间体在边界上受法向分布力q2023/8/1436第36页，课件共57页，创作于2023年2月2．位移：z轴上的ur=0，仅存在w第三节半空间体在边界上受法向分布力q2023/8/1437第37页，课件共57页，创作于2023年2月第三节半空间体在边界上受法向分布力q2023/8/1438第38页，课件共57页，创作于2023年2月第四节两球体之间的接触压力接触压力问题是在机械工程、土木工程中经常碰到的问题，接触问题在1881年由德国赫兹（HeinrichHerty）首先用数学弹性力学导出了计算公式。4.1接触问题的特点：

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

2023/8/1440第40页，课件共57页，创作于2023年2月

4.2

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

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

、

1、R1球2：E2

、

2、R2

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

M2点的z坐标分别为（M1和M2与点o很近）2023/8/1441第41页，课件共57页，创作于2023年2月第四节两球体之间的接触压力rOz1z2O2O1R2R1M1M2r则2023/8/1442第42页，课件共57页，创作于2023年2月第四节两球体之间的接触压力在已知P压力作用下，两球在接触点附近发生变形有一个接触面，根据对称性接触面为以a为半径的圆。rOz1z2O2O1R2R1M1M2rM1rPPoz1z2O1M2ar2023/8/1443第43页，课件共57页，创作于2023年2月第四节两球体之间的接触压力1．a为待求量，同时接触面上有接触压力q（待求）。2．由于接触问题是局部变形，在球体远离o点的任意点位移为刚体位移。两球内距o点很远处的相对位移（刚体位移）为

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

、q

和

。2023/8/1444第44页，课件共57页，创作于2023年2月第四节两球体之间的接触压力求解：首先根据接触面变形（位移）来建立一个关系球1：触面上o点、M1点沿z1轴位移为w1(o)、w1而w1(o)=w1+z1

M1rPPoz1z2O1M2ar2023/8/1445第45页，课件共57页，创作于2023年2月第四节两球体之间的接触压力球2：触面上o点、M2点沿z2轴位移为w2(o)、w2w2(o)=w2+z2

而w1(o)+w2(o)=w1+z1+w2+z2w1(o)+w2(o)=w1+w2+

r2或M1rPPoz1z2O1M2ar2023/8/1446第46页，课件共57页，创作于2023年2月而w1(o)+w2(o)=

第四节两球体之间的接触压力——两球体距o点较远处两点的趋近距离。

=w1+w2+

r2——变性协调关系w1(o)+w2(o)=w1+w2+

r2由于接触问题可看成半无限体受局部垂直分布力问题，w1和w2可以利用上一节的结果。M1rPPoz1z2O1M2ar2023/8/1447第47页，课件共57页，创作于2023年2月第四节两球体之间的接触压力—相当物理和几何关系2023/8/1448第48页，课件共57页，创作于2023年2月代入

=w1+w2+

r2在此式中a、q和

未知。第四节两