1、2018/11/20,1,第九章 空间轴对称问题,本章讨论空间轴对称问题的基本方程和一些轴对称问题的基本解。对于一般空间问题的解法我们在第五章已有讨论，但一般空间问题一般解（具体求解）通解讨论在杜庆华等编著的“弹性理论”中有较多的论述。我们不刻意从数学上论述一般空间问题一般解的表达式，而对于空间轴对称问题作一些讨论和举例。,2018/11/20,2,1.1空间轴对称问题特点：,1. 域内所有物理量（体力、面力、位移、应力、应变）均为r、z的函数。,与平面轴对称问题类似，空间轴对称问题的求解域、荷载和约束绕某一轴（z轴）对称，导致如下简化，,2荷载：体力f=0，面力 ，位移u=0，应力 r= z
2、=0，应变 r=z=0。,第一节 空间轴对称问题的基本方程,2018/11/20,3,第一节 空间轴对称问题的基本方程,3待求的物理量（10个） ：ur、w、r、z、 rz= zr、r、z、 rz=zr,1.2基本方程 1.平衡微分方程（两个）：,2018/11/20,4,2.几何方程（四个）：,第一节 空间轴对称问题的基本方程,3.变形协调方程（四个）,2018/11/20,5,4.物理方程（四个）：,第一节 空间轴对称问题的基本方程,2018/11/20,6,r=e2Gr、 =e2G、 z=e2Gz、 rz=G rz,第一节 空间轴对称问题的基本方程,其中 体积应变,或,2018/11/2
3、0,7,5.边界条件,第一节 空间轴对称问题的基本方程,位移边界： 在Su上,6.按应力解法,力的边界：在 r=r0 在 z=z0,四个应力分量r、z、 rz 为基本未知量。,2018/11/20,8,基本方程（六个）： 两个平衡微分方程与四个用应力表示的变 形协 调方程； 再加上力的边界条件。,第一节 空间轴对称问题的基本方程,如果体力为零时，基本方程为齐次方程，则可采用应力函数解法，引入应力函数 (r,z) ，使得应力用 (r,z) 表示:,2018/11/20,9,第一节 空间轴对称问题的基本方程,(r,z)满足第一个平衡微分方程，而第二个平衡方程及四个相容方程，共同要求 2 2= 4
4、=0 (r,z)应满足的基本微分方程。,2018/11/20,10,7按位移法解,第一节 空间轴对称问题的基本方程,其中,a基本未知函数： ur和w,基本方程两个： 并考虑适当的边界条件。,2018/11/20,11,b. 引入Love（拉甫、勒夫）位移函数（当无体力作用时）,第一节 空间轴对称问题的基本方程,对于位移法的基本方程的解可由考虑体力的一个特解加上齐次方程的通解。,轴对称问题齐次拉梅方程的通解可以引入一个Love位移函数(r,z),使得位移由(r,z)表示：,2018/11/20,12,代入齐次拉梅方程，第一式自然满足，而第二式为基本方程： 4=0 (r,z)为双调和方程。,第一节
5、 空间轴对称问题的基本方程,同时应力分量由(r,z)表示为：,2018/11/20,13,轴对称问题按位移求解，归结为寻找一个恰当的重调和函数(r,z)，使按其导出位移和应力能满足给定的边界条件。,第一节 空间轴对称问题的基本方程,比较应力函数解法和love位移法知： (r,z)= (r,z),2018/11/20,14,第二节 半空间体在边界上受法向集中力 （Boussinesq问题）,半空间体，体力不计，边界受法向集中力P作用. 轴对称问题，P作用在坐标原点上。,已知，当z=0且r 0时, z=0 , zr= 0；,当R 0时，应力奇异。,当R 时，R=(r2+z2)1/2， 应力、位移
6、0；,2018/11/20,15,选 (r,z) 为r和z的正一次幂式: (r,z) = A1R+ A2R - zln(R+z) 为双调和函数,第二节 半空间体在边界上受法向集中力 （Boussinesq问题）,Boussinesq采取Love函数求解，(r,z)为重调和函数，由(r,z)的三次微分导出应力。,2018/11/20,16,(r,z) = A1R+ A2R - zln(R+z) 则 (r,z) 自然满足 4=0 。 代入位移、应力计算式.,第二节 半空间体在边界上受法向集中力 （Boussinesq问题）,位移：,2018/11/20,17,应力：,第二节 半空间体在边界上受法向
7、集中力 （Boussinesq问题）,2018/11/20,18,根据边界条件来确定A1和A2：,第二节 半空间体在边界上受法向集中力 （Boussinesq问题）,在z=0且r 0边界上, z=0 自然满足。,在z=0且r 0边界上, zr= 0 ,（1-2）A1+ A2 = 0(a),2018/11/20,19,在z= z0 0平面上，要求z 的合力与P平衡。,还需一个条件（包括P的）。,第二节 半空间体在边界上受法向集中力 （Boussinesq问题）,将z 表达式代入，得,2018/11/20,20,P - 4A1(1-）- 2 A2 = 0 (b),第二节 半空间体在边界上受法向集中
8、力 （Boussinesq问题）,而,2018/11/20,21,由式(a)、(b)解得 A1 = P/(2) 、A2 = -（1-2）P /（2）,第二节 半空间体在边界上受法向集中力 （Boussinesq问题）,代回位移、应力表达式，见徐芝纶（上册）P.297（917）、（918）式，称为Boussinesq问题解。由P.297（917）、（918）式见：位移和应力随R 的增加而减小。,2018/11/20,22,第二节 半空间体在边界上受法向集中力 （Boussinesq问题）,在z=0平面上,2018/11/20,23,第三节 半空间体在边界上受法向分布力q,已知条件：半空间体在边界
9、上受均布法向荷载q作用，在半径为a的圆面积。,寻求解答： 1. z =0边界上的沉陷 w z=0 = ? 2. r =0（对称轴）上的应力和位移。,求解方法：采用叠加法和半空间体边界受法向集中力P的计算结果求解。,2018/11/20,24,3.1 边界上一点M的竖向位移w：,第三节 半空间体在边界上受法向分布力q,1.设M点为圆面积之外：,M点可以在荷载圆面积之外也可在之内。,当半空间体边界上受法向集中力P时，边界上距P点为r的点竖向位移为 ：,2018/11/20,25,圆面积均布荷载q对圆外M 点竖向位移影响可取一个微面元，距M点为s，角度为 处，dA=sdds ，dA上q 对M点影响：
10、,第三节 半空间体在边界上受法向分布力q,2018/11/20,26,第三节 半空间体在边界上受法向分布力q,2018/11/20,27,整体圆面积荷载对M点影响为,第三节 半空间体在边界上受法向分布力q,而,2018/11/20,28,1为M点作为圆相切线OM线的夹角,第三节 半空间体在边界上受法向分布力q,为了简化积分将积分变量 转变为 ,2018/11/20,29,由图形可见 asin=rsin , 两边微分 acosd = rcosd,第三节 半空间体在边界上受法向分布力q,2018/11/20,30,第三节 半空间体在边界上受法向分布力q,的取值范围：由0 1,的取值范围：0 ,20
11、18/11/20,31,第三节 半空间体在边界上受法向分布力q,2018/11/20,32,第二类椭圆积分 第一类椭圆积分,第三节 半空间体在边界上受法向分布力q,对于不同a/r可由椭圆积分表得到。,2018/11/20,33,2M点载荷在圆之内：,第三节 半空间体在边界上受法向分布力q,圆内距M点s处微面积q对M点沉陷的影响仍为,2018/11/20,34,整个圆面积荷载引起M点沉陷为：,第三节 半空间体在边界上受法向分布力q,第二类椭圆积分,利用 asin=rsin,2018/11/20,35,当r= 0为圆心处沉陷：当r= a时圆周上沉陷：,3.2 在z轴r=0上的应力和位移在z轴上的应
12、力和位移比同一水平面上其它点的应力和位移要大。,第三节 半空间体在边界上受法向分布力q,2018/11/20,36,1. 应力：由于z轴对称轴，所以在z 轴上的应力无剪应力，均为主应力： r = 、z,第三节 半空间体在边界上受法向分布力q,2018/11/20,37,2位移：z轴上的ur= 0，仅存在w,第三节 半空间体在边界上受法向分布力q,2018/11/20,38,第三节 半空间体在边界上受法向分布力q,2018/11/20,39,第四节 两球体之间的接触压力,接触压力问题是在机械工程、土木工程中经常碰到的问题，接触问题在1881年由德国赫兹（Heinrich Herty）首先用数学弹
13、性力学导出了计算公式。,4.1 接触问题的特点： 1两个弹性体互相接触，当无压力作用时，为点接触或线接触。当有压力作用时，弹性体发生变形，点接触（或线接触）变为面接触。,2018/11/20,40,2弹性体变形后的接触面为非常小的局部区域（相对于弹性体几何尺寸）所以可看成半空间（半无限平面）体法向受局部分布力作用问题，但这里分布力q不是均匀的，同时q也未知，接触面的局部区域也是未知的。,第四节 两球体之间的接触压力,3不计接触面摩擦力。,2018/11/20,41,4.2 两球体之间的接触压力： 已知两球体变形前在o点接触，两个坐标系 roz1、roz2,第四节 两球体之间的接触压力,球1：E
14、1 、1、R1 球2：E2 、2、R2,距接触点z轴为r的两球 表面上M1和 M2点的z 坐标分别为（M1和M2与点o很近）,2018/11/20,42,第四节 两球体之间的接触压力,则,2018/11/20,43,第四节 两球体之间的接触压力,在已知P压力作用下，两球在接触点附近发生变形有一个接触面，根据对称性接触面为以a为半径的圆。,2018/11/20,44,第四节 两球体之间的接触压力,1a为待求量，同时接触面上有接触压力q（待求）。,2由于接触问题是局部变形，在球体远离o点的任意点位移为刚体位移。两球内距o点很远处的相对位移（刚体位移）为 ？,下面要建立（找出）三个条件（几何、物理、
15、平衡方程）寻求a 、q 和。,2018/11/20,45,第四节 两球体之间的接触压力,求解：首先根据接触面变形（位移）来建立一个关系,球1：触面上o点、M1点沿z1轴位移为w1(o)、w1,而 w1(o)= w1+ z1,2018/11/20,46,第四节 两球体之间的接触压力,球2：触面上o点、M2点沿z2轴位移为w2(o)、w2,w2(o)= w2+ z2,而 w1(o) +w2(o)=w1+ z1+w2+ z2,w1(o) +w2(o)=w1+w2+ r2,或,2018/11/20,47,而 w1(o) +w2(o)=,第四节 两球体之间的接触压力,两球体距o点较远处两点的趋近距离。,
16、 = w1+w2+ r2 变性协调关系,w1(o) +w2(o)=w1+w2+ r2,由于接触问题可看成半无限体受局部垂直分布力问题，w1和w2可以利用上一节的结果。,2018/11/20,48,第四节 两球体之间的接触压力,相当物理和几何关系,2018/11/20,49,代入 = w1+w2+ r2,在此式中a 、q 和 未知。,第四节 两球体之间的接触压力,2018/11/20,50,第四节 两球体之间的接触压力,q与P 有关，为寻求解，赫兹假设接触面上的分布力q的。,假设： q 分布为以a为半径的半球面乘q0/a，q0为接触面中心接触压力的集度。,2018/11/20,51,第四节 两球
17、体之间的接触压力,2018/11/20,52,第四节 两球体之间的接触压力,赫兹通过这样假设，并利用,2018/11/20,53,或,第四节 两球体之间的接触压力,得,代回,赫兹通过接触面上的接触压力的分布假设可使等式右端的积分为一个常数项和r2的二次项。,2018/11/20,54,第四节 两球体之间的接触压力,的积分，在任意,2018/11/20,55,第四节 两球体之间的接触压力,2018/11/20,56,积分得,第四节 两球体之间的接触压力,比较上式两边得,2018/11/20,57,将 代入可确定a 和。,第四节 两球体之间的接触压力,q0、a 和 表达式见徐芝纶弹性力学（上册）P.308（926） 、（927）式 ，q0、a 和 确定后，可求球体内应力。,
