1、第六章 弹塑性平面问题127第六章 弹塑性平面问题任何一个弹塑性体实际上都是空间(三维)物体,且一般的载荷严格说来也是空间力系。因此,所有弹塑性力学问题实际上都是空间问题,即所有的力学量都是坐标 的函数。但是,当所考察的物体(结构)及其所承受的载荷具有某些,zyx特点时,则可将它们近似地看作平面(二维)问题,即所有的力学量都是两个坐标(如 )的函数,从而使问题得简化,且所得解答又具有工程所要求的精度。,由第二章知,弹塑性力学平面问题可分为平面应力问题和平面应变问题两种,本章主要讨论弹塑性平面问题求解的一般方法。6.1 弹性平面问题的基本方程由第二章己经知道,两类平面问题的基本未知量虽然是完全相
2、同的,但非零的应力分量、应变分量和位移分量不是完全相同的。1.1 平衡方程无论是平面应力问题还是平面应变问题,由于在 方向自成平衡,因此,两z类问题的平衡方程均为(6.1-1)0YyxXxy1.2 几何方程由于只需要考虑面内的几何关系,因此,对于两类平面向题均有(6.1-2)xvyu、yvxux由式(6.1-2)可得到平面问题的变形协调方程为(6.1-3)yxyyx221.3 本构关系两类平面问题的非零应力分量和应变分量不相同,因此,由广义虎克定律所得本构方程也必然不尽相同。(1) 平面应力问题第六章 弹塑性平面问题128对于平面应力问题,因 ,根据广义虎克定律显然有,0z0zxy。因此本构方
3、程为0zxy(6.1-4a)xyxyzxyyxEE)1(2)(1或(6.1-4b)xyxyyxxE)1(2)(12(2) 平面应变问题对于平面应变问题,有 ,根据广义虎克定律,必有0zxyz)(yE和 。因此,本构关系为0zxy(6.1-5a)xyxy xyxxE)1(2)(12或第六章 弹塑性平面问题129(6.1-5b)xyxyz xyy yxxEE)1(2)1()(21)(将上面两种平面问题的本构方程式进行比较可以看出,只要将平面应力问题本构方程式中的 换为 , 换为 ,就可以得到平面应变问题的/)1(本构方程式。1.3 应变协调方程如果采用应力法求解,还必须将平面问题的应变协调方程(6
4、.1-3)式变换为用应力表示。(1)平面应力问题的应变协调方程对于平面应力问题,将方程(6.1-1)式中的第一式对 求导,第二式对 求导,xy有 yYyxxXxy2222将上式相加后,得 )(21)(212 yxXyxyx因 )()( 2222 yYxyxEyxyx 将式(6.1-3)中耐 、 用本构关系式(6.1-4a)入,而 用上式代换,可得xy 0)(1222222 yYxXyxxyyyx 化简上式,得 )(2222 yxyxyx第六章 弹塑性平面问题130上式可进一步写为(6.1-6)(1)(2 yYxXyxyx 如果不计体力或为常体力,则上式可写为(6.1-7a)0)(2yxx或用拉
5、普拉斯算符简写为(6.1-7b)(2yx式(6.1-6)即为用应力表示的应变协调方程,通常称为纳维方程。(2)平面应变问题的应变协调方程对于平面应变问题,因为平衡方程同样为(6.1-1)式,应力分量 、 也只xy是 、 的函数,因此应用由平面应力变换到平面应变的对应关系,则平面应变xy问题的应变协调方程可直接从(6.1-6)中得到,即(6.1-8)(1)(2 yYxXyxyx 注意到,当在平面应变问题中,如果不计体力或为常体力时,则(6.1-8)式也简化为(6.1-7)式,这时平面应力问题与平面应变问题的应变协调方程相同。由以上可见,如果讨论的问题为 域上的调和函数,则 是在区域D)(yx上直
6、到二阶导数都是连续的连续函数。在这种情况下,平面应力和平面应变问D题的应力分量 , , 的分布是相同的,是是是说,他们在 平面内应力xyx o场一致。1.4 边界条件平面内周边上的应力边界条件为YmlXyx(6.1-9a)对于平面应变问题还有(6.1-9b)Zz对于平面应力问题由于 方向无外力作用,又 ,所以该方向的边界0z条件自动满足。从以上的讨论中不难发现,方程(6.1-1 和(6.1-7)以及边界条件(6.1-9) 中均不含材料常数。由此得出重要结论:对于全部边界为力边界的无(或常)体力的平面问题,无论什么材料,只要它们的几何条件、载荷条件相同,则不论其为平面应力第六章 弹塑性平面问题1
7、31或平面应变问题,他们在平面内的应力分布规律是相同的。这一结论,给实验模型的设计,尤其是光弹性实验提供了理论基础、并具有很大的灵活性。但需特别注意的是,以上两种情况的应力 、应变和位移是不相同的。z6.2 应力解法与应力函数由以上讨论可知,当边值问题属于第一类,即面力已知问题,则采用应力法求解时,其基本方程归结为(1) 当体力为常量时 0YyxXxy(6.2-1a)(2yx (6.2-1b)YmlXyx(6.2-1c)(2) 当不计体力时(6.2-2a)0yxxy)(2yx (6.2-2b)YmlXyx(6.2-1c)由数学上可知,方程(6.2-1)是一组线性非齐次偏微分方程,它的解答应该包
8、含两部分:任意一组特解和齐次方程(6.2-2a)的通解。非齐次方程(6.2-1a) 的特解可取为第六章 弹塑性平面问题132(a)0xyxYX或取为(b)YxXyx或取为(c)0xy等形式。显然,这些特解都满足(6.2-1a)式。对于齐次方程式(6.2-2a),如果引进一个函数 ,使得),(yx(6.2-3)yxyxyx22则将(6.2-3)式代人齐次方程(6.2-2a)式,可知恒满足。函数 称为平面向题的),(yx应力函数,是英国天文学家艾里(Airy,GB)于 1862 年首先提出的,因此也称它为艾里应力函数。将(6.2-3)式与式(a)式相叠加,就得到(6.2-1a)式的全解为(6.2-
9、4)yxYXxyxyx22为使应力表达式同时满足协调方程、则应力函数 还必须满足一定的条),(yx件。将(6.2-4)式代入(6.2-1b),得0)(222yxyx将上式展开为第六章 弹塑性平面问题133(6.2-5a)024244yx或采用双调和算子简写为(6.2-5b)4将(6.2-4)式代入(6.2-1c)式,得到相应的用应力函数表示的静力边界条件为(6.2-6)YmyxlyXlX)()(222综上所述,对于常体力下的平面问题,只要求解一个未知函数 ,即在),(yx给定边界条件(6.2-6)的情况下,求解方程(6.2-5) 式。求出函数 后,就可通过(6.2-4)式求出应力分量,最后可通
10、过本构方程式求应变,通过几何方程式积分求位移。对于无体力的平面问题情况,协调方程式不变静力边界条件利用(6.2-3)式写为(6.2-7)YmxlyXl)()(22相应的应力分量为(6.2-3)式所示。实际上,直接求解弹性力学问题往住是很困难的,因此有时不得不采用逆解法或半逆解怯等来求解。当用逆解法时,需先假定满足双调和方程(6.2-5)式的某种形式的应力函数,然后用式(6.2-3) 或(6.2-4) 求出应力分量 , , 等,再根据边界条),(yx xyx件式(6.2-6)或(6.2-7)来分析所得应力分量对应于什么样的面力。由此判定所选应力函数 可以解什么样的问题。,如用半逆解法则针对所要求
11、的问题,假定部分或全部应力分量为某种形式的双调和函数,并引入足够多的待定参数,从而导出应力函数 ,然后分析所),(yx得应力函数是否满足应变协调方程,判断假定的以及由应力函数导出的应力分量是否满足边界条件。如不满足则应重新假定。应当注意的是,双调和方程是四阶的或低于四阶的多项式都是双调和函数。但必须至少是二次和二次以上,以保证得出非零的应力解。例如,对于不计体力的弹性平问题,如取应力函数 的一次式,尽管它满足双调yaxyx210),(第六章 弹塑性平面问题134和方程式(6.2-5),但将其代入应力分量(6.2-3) 式,得 0022yxyxyx显然,这是一个无应力状态。由此得出,在应力函数中
12、增添或除去 和 的一次xy式,并不影响应力分量。不难验证,当应力函数取二次多项式时可得均匀应力状态,取三次多项式时得线性分布的应力场。.3 梁的弹性平面弯曲3.1 悬臂梁的弹性平面弯曲在第五章采用材料力学初等理论介绍了梁的弹塑性纯弯曲,这节首先应用应力函数法讨论高为 ,宽为 ,跨长为 ,如图 6.1 所示的悬臂粱在自由端受集中hbl力 作用,忽略自重时的平面弯曲。Fo2/hxly图 6.1 悬臂梁对于图示悬臂梁,其边界条件为(a)2/0)(hxyyxFbd式(a)所列边界条件表示:悬臂梁自由端没有轴向水平力,顶部和底部没有载荷作用,及自由端的切应力之和应等于 F。 (a)中第四式的负号是因此处
13、切应力是作用在外法线方向与 轴反向的平面内,切应力方向又与 轴同向,根据第 2 章对切x y应力的正负号约定应为负。第六章 弹塑性平面问题1351)选择应力函数由材料力学可知,悬臂梁任一截面上由 产生的弯矩随 作线性变化,而且Fx截面上任一点的正应力 与 成比例,因此,可假定 为xy(b)xyAx12式中 为常数。将(b) 式对 积分两次,得1Ay(c)(6),( 2131xfyf式(c)中的 和 为 的待定函数。将(c) 式代入双调和方程(6.2-5a) 可得)(1xf2fx(d)04241dxffy因 和 仅为 的函数,而上式中左边第二项又与 无关,故要使上式成立时,1f2x y必有 0,
14、04241dxfdxf对上面两式分别积分,得 9827362 5421)( AxxAf式中 系积分常数。将它们代入式(c) ,可得应力函数为iA(6.3-1)( )(), 982736 54231AxxAyy将式(e)代入式(6.2-3),可得应力分量 为y、(e)2321()432 76AxyAxxyy2)确定系数 iA根据边界条件式(a)中的第 2 式,有 0)2()2(6736AhxAh第六章 弹塑性平面问题136上式应对所有的 都应成立,因而必有x027362Ah求解此方程组,得(f)7632根据边界条件式(a)中的第 3 式,并注意到式(e)和(f),则有 0)81()(422Ahh
15、yx由上式可得 214又依据边界条件式(a)的第 4 式,可得 FbdyhAbdyhhx )(2122由上式可得 IF/)/(31式中 为梁截面对中性轴的惯性矩。123bhI3)应力分量计算至此,式(b),(e) 中的所有常数均已确定,于是可得悬臂梁中的各应力分量为(6.3-2)4(202yhIFIxyxyx式(6.3-2)的结果与材料力学的结果完全一致。由此可得出结沦:如果自由端部的切力按抛物线分布, 在固定端是按线性分布,则这一解是精确解。如果不x是这样,根据圣维南原理,这一解在梁内远离端部的截面仍是足够精确的,其所影响的范围大约只有截面尺寸大小的长度。需注意的是,式(6.3-1)中的系数
16、 并未求出,由上节已知,这主要是985,A这 3 个系数与应力分量无关。因此,这几个系数确定与否无关紧要。第六章 弹塑性平面问题1374)变形计算当应力求得后,变形计算则可根据应变位移几何关系和虎克定律进行。由式(6.3-2)可得(g) )4(2)4()1)(22yhGIFyhEIxvyuIxyEFxuxyyyx将式(g)中的第 1 和第 2 式分别对 积分,有x,(h)(212xvyEIFvuu将式(h)分别对 微分,代入式(g) 的第 3 式,并整理后可得xy, 22121 4)1()()()( hEIFxIdvyIdu 上式两边分别 与 的函数,因此等式左、右两边应等于同一常数 ,即xy
17、 1a12214)1()( ahEIFxIdvyu将上式积分后代入式(h),可得位移 的表达式为vu,(k)3123214)(6axhEIFxIyEIFvyau式中常数 由阻止梁在 面内作刚体运动所必需的三个约束条件来确定。下面iao分两种情况进行讨论。(1)固定端处 ( )的边界条件为:0,ylx 0xvu这一边界条件相当于在固端处梁的轴线的切线保持水平, ,即将坐标点(第六章 弹塑性平面问题138)的水平微线段 固定,这与材料力学的处理方法相同。现将该边界条0,ylxdx件代入位移表达式(k),可得 332221 04)( lEIFa、ahEIFla 将这些系数代入式(k),则位移为(6.
18、3-3)3223 264)1(6)()( lEIFxyIEFlxIv yhyu由该式可知, 均是 的非线性函数,这说明梁的任一截面变形后不再保持为u,y,平面,这与材料力学初等理论所得结果不同。如在固定端( )处,由式(6.3-3)lx可得 GIFhEIyuyFIlvlxlxl 24)1()( )(2)(0 22上式表明,由这种固定条件得到铅垂线元有一绕垂直于面 平面的 轴逆时针oxyz方向的转角 (图 6.2)。 梁轴线的铅垂位移由(6.3-3)式可得)2/GIFhEIFlxlIvy36)(20对于梁自由端( )处的挠度, x由式(6.3-3)的第 2 式可得梁 EIFlv3)(max这与材
19、料力学的结果相同。 图 6.2 固定端转角示意图(3) 固定端处( )的边界条件为:0,ylx 0yuv这表示固定端断面在( )处的铅垂微线元 固定不能转动,将该边,l d第六章 弹塑性平面问题139界条件代入式(k),可求得 233221 4)1(0lhEIFlIa、aEIFla 将这些系数代入式(k),得到一组与(6.3-3) 式不同的梁的位移为(6.3-4) )(14)(21366)(2232 xlhxlylxEIFvEIFu实际上,由式(6.3-4)也得出 均是 的非线性函数。同样可得出在固定端vu,( )处的水平线元 也有一绕垂直于面 平面的 轴逆时针方向的转角0,ylxdxoyz(
20、如图 6.3 所示)。2/GIFh由(6.3-4)式可得梁轴的铅垂位移为 )()1(36)( 220 xlhEIFlIxlEvy 自由端的挠度度为图 6.3 固定端转角示意图GIlhFvyx23)(0显然,上式等号右边第二项是剪力对挠度的影响。而这部分与弯曲的影响之比,为 2232 )()1(4)/( lhlEIFlh如 ,则此比值为 。所以当 时,梁的挠度主要由于弯曲)2(10l0所引起。由此可见,在材料力学中得到的结果,对于细长梁是精确的。但是,必须指出,在高而短的梁中,以及在梁的高频振动和在波的传播问题中,切力效应是非常重要的。由以上计算变形可见,在材料力学中,只是笼统地说梁端“固定”
21、,没有规定具体的固定方式。在弹性理论中必须规定固定的方式,根据不同的固定方式,得出不同的位移公式。3.2 简支梁的弹性平面弯曲现分析如图 6.4 所示受均布载荷 作用,不计体力的两端简支梁。q第六章 弹塑性平面问题1402/h xll b y图 6.4 受均布荷重简支梁1)选择应力函数设应力函数为(a)5(),( 324321 yxAyxAyx 将(a)式代入式(6.2-5)知,满足双调和方程。2)利用边界条件确定常数 i边界条件为(b)0)(/2/2/hyxhybq和在 处lx(c)2/ /hxybqld根据式(6.2-3)及边界条件式(b) 、(c) 可得(d)2/(401/3)/(/4/
22、22212bqhAlA由此可解得系数 为i(e)343221 )10(43bhqA、bhl、bhqA 将(e)式代入(a)式,得应力函数q第六章 弹塑性平面问题141(f) )5(1)0(43),( 32322 yxhyhlyxbqyx3)应力分量计算将式(f)代入式(6.2-3),并注意梁的截面惯性矩 ,求得应力分量为/3bI(g)xyhIqyhqyxlIqxyyx)4(2123)0()(2322将式(g)的应力分量与材料力学对该问题的解答相比,可以看出:1)式(g)中 的表达式包括两项,第一项与材料力学解答相同,而第二项与x无关,是对材料力学解答的修正。且当 时,梁的端面有正应力 ,但端x
23、 lxx面上没有水平外力,所以 的表达式只满足了两端弯矩为零的条件,但未能消除x两端的正应力。然而这组附加的水平力,即修正项显然也构成平衡力系。根据圣维南原理,这组附加力的效应是局部的,在远离两端部分可认为材料力学的公式是精确的。因此,通常认为长而低的细长粱此项可忽略不计,但对高梁,即短粗梁这项有显著影响。如梁中间截面处 , 的最大值为)2/,0(hyxx1532llqx上式括号中的第一项为主要应力,第二项反映修正应力。一般认为当 时,hl52材料力学公式不再适用。2)材料力学对该问题剪应力的解答与本精确解完全吻合。3) 材料力学假设梁的纵向纤维之间互不挤压力,因此 ,但本解答表明,0y除梁截
24、面的下表面处外,其余部位 ,且 无关,因此整个梁的纵向纤维之0yx间均存在挤压力。梁中应力分布如图 6.5 所示。图 6.5 简支梁截面应力分布示意图4)位移计算第六章 弹塑性平面问题142对于位移计算,其方法和步骤与悬臂梁的位移计算相同,即利用本构方程和几何方程,并假定梁中间截面的形心 的水平位移等于零,而垂向位移)0(yx为 ,经积分后可得(h) 12852 )2061(1 43132242 224 322xhxl yhylyyEIqv xhxlu由上面的位移表达式可以看出:(1)由位移 的表达式可知,梁的中性层并不在梁截面的中间层。在中间层u处有水平位移0y(k)xhEIquy12)(3
25、0由式(g)知,当 时, ,因此沿 方向引起拉伸应变0y23Iy12)(30hEIqyyx将上式积分,并注意当 时, ,则得式(k)。u(2)由位移 可得梁轴线( )的挠曲线方程为v(m)12852(0 xhlEIqxy 假设在梁中心轴线的两端( )处,垂向位移 ,则得l0)yv24(5312lhIq上式右端括弧中的第一项所反映的挠度与材料力学根据平面假设而得出的结果相一致,括弧中的第二项反映剪应力的影响。(3)将式(m)对 求二阶导致,得挠曲线的曲率方程式x )254()(21)(202 hxlEIqvy由该式可知,曲率并不与 成正比,方括号半的第二项是对材料力学近)xl第六章 弹塑性平面问
26、题143似曲率公式 的修正。EIM1必须指出的是,应用式(g)也能解答梁两端固定的问题,为此,须取适当的支座反力矩 ,并使粱的二端都满足 或 的条件。0 0xvyu应注意,用多项式求解仅对低粱适用,对于高梁,两端的平衡力系要影响到跨度中部的应力。因此,须用其他形式的应力函数,例如三角级数的形式。6.4 用三角级数求解弹性平面问题前面一节是用代数多项式为应力函数求解弹性梁平面问题,对于一端受集中力悬臂梁的弯曲,用三次多项式为应力函数;对于受连续均布载荷的单跨粱,用五次多项式。增高多项式的幂次,可以求解受载荷更复杂的问题。但是整多项式只限于求解梁上载荷的分布是连续的,而且分布规律能用代数整函数表示
27、的一些简单问题。如果载荷分布不是连续的,而且分布规律不能用代数整函数表示(图6.6),则可以用三角级数求解这类问题。载荷函数可展开为三角级数,应力函数也可以用三角级数表示。下面讨论用三角级数求解弹性平面问题。图 6.6 受非连续分布载荷深梁示意图4.1 应力函数的确级数表达对于图 6.6 所示梁,取应力函数为(a)xyfxsin)(,(式中 , 是任意整数。将式(a)代入双调和方程(6.2-5a),得常微分方程ln(b)0)()(2)(4“4 ffyf该常系数线性微分方程的通解为(c)yDshCycBshAcf 将(c)式代入(a)式,得应力函数第六章 弹塑性平面问题144(d)(sin),(
28、 yDshCycBshyAcxy 当不计体力时,由相应的应力分量为(6.4-1)( )(cos(in)2 )2(si2 22 yhDyshcCBAxyx Dyhsychssyxyyx 由于 为任意整数,因此可得无穷多个函数 ,又因双调和方程(6.2-5a)是n 线性的,所以无穷多个函数 之和也是双调和方程(6.2-5a)的解,即(6.4-2a)(),(1 yshDycCshBycAxisyx nnnnn 不难证明,如取应力函数 为(6.4-2b)(os),( 1 ssnnnnn 也能满足双调和函数。因此,式(e)和(f)之和也能满足双调和函数,所以应力函数 可表示为如下级数形式(6.4-2c)
29、 1 )(cos),(n nnnn yshDycCshByAxiyx 4.2 载荷函数式(6.4-1)中的系数 ,和 需根据边界条件确定。这时应力边界条件必须ii展开为无穷级数的形式(6.4-3)xbxaxqnnsicos)(110由数学分析知,将一函数在区域 展开成富里叶(Fourier)级数(6.4-3)l,时,其系数(称作富里叶系数)为(6.4-4)lnlldxlnqbadxqsi)(1co)(20第六章 弹塑性平面问题145如果在梁的 边界上作用有均布载荷 ,则富里叶系数为a,0q(6.4-5)l anll adxlnlqdxlnqbldxlxq0sisi)(12co)(21000 这
30、样,对于该边界有(6.4-6)cosin12()(0nlxlalaqx例 6.1 设一简支梁的中部上、下两表面,在 范围内对称地作用均布载荷2(如图 6.7 所示)。如此梁的q厚度为 1 个单位,不计体力,试求其应力分量。解:首先将载荷展开为富里叶级数,最普遍的情况下,上部边界( )和下部边界ty( )的载荷分别表示为 图 6.7 局部受均布载荷简支粱t(d) lxnGlxnGq lElEnbtydtyu cossi)(i110 注意载荷实际作用区域为(e)(, axql、ldu式中 表示整个梁的均匀分布载荷,式(e)中的全部系数均可(6.4-4)中的富0,GE里叶系数的公式求出。由图 6.7
31、 可知,所示载荷对称于 轴,是 的偶函数,故式(d)的展开式只含yx及余弦项,其中0,(f)al lqdxlxqGE2)(210而系数 可由载荷展开式nE第六章 弹塑性平面问题146(g)lxnExqncos)(1运用通常求富里叶系数的办法,两边乘以 ,并在区间 积分,有lm,l )(0coscos)(1 nmldxlldxlmxql nl 由此可得 lmxlxqEcs)(由于 为任意整数,所以可换成 ,于是得mnl dlnos)(1同理也可得 。将 代入上式可得nGqx)(h)alnqxllEan si2cs 由于常数 的存在,该问题可理解为上、下分别作用均布载荷,0,,再加上后面的三角级数
32、所表示的载荷。于是,可以分别计算每lqaE0一部分载荷所产生的应力,然后再叠加。对于上、下面作用均布压缩载荷 ,相应的应力分量为lqa(i)lyxy,0而 , 这些载荷所产生的应力分量,可依alnqEnsi2alnGnsi2据应力函数表达式(6.4-2b)和(6.2-3)求得,即(j) )( )(sin(co)2 )2(s 122 212 ychyD yshycCychBsAxyx Dyhs ychysCBAxyn nnnxy nny n nnnnnx 式中 ,各个常数可由边界条件确定(参见式(g),即ln第六章 弹塑性平面问题147(k)1cos)4(,0)3(: 2,1nndyxyuyxy
33、 xGqty E将式(k)中的式(1),式(3)代入式(j),并注意到双曲函数的关系式,则可得chshs)(,)(m)tchtsADtstBCnnnn将式(k)中的式(2)、(4)代入式(j),并与式(k)联立求解,并注意 ,得nBA(n)tshEDtctsACBnnnnnn220由于 已经求得,所以式(j)中的常数可全部确定,将式(n)代入式(j),即nE得相应的应力分量,再加上式(i)中由均布载荷而产生的应力,即得梁总的应力分量计算式。如 的表达式为yxytsh ychttshtnaqlann nnco )()2(i41其余应力分量由读者自行导出。如果在该梁上的分布载荷 作用范q围不断缩小
34、,即随着这短段的缩小达到极限情况,就得到梁受两个相向集中压力的情形(图 6.8a),在这种情况下的应力 沿 方向的分布曲线如图 6.8b 所示,yx由该图可见, 随 的增大而迅速衰减。yx这一计算实例,可以说明圣维南原理对此也是是正确的。 图 6.8 分布曲线y第六章 弹塑性平面问题1486.5 极坐标系下的基本方程对于圆形或部分圆形(扇形,楔形等)的物体,用极坐标求解比较方便。在极坐标系中,平面内任一点的位置,用径向坐标 及周向坐标 来表示(如图 6.9)。r极坐标系 与直角坐标系 之间的关系为,r),(yx(6.5-1a)sincoyx(6.5-1b)xr122ta下面推导极坐标平面问题的
35、基本微分方程。 图 6.9 极坐标与直角坐标间的关系5.1 平衡方程在变形物体中,用两个同心柱面和两个径向平面截割出微小单元体(见图 6.10)。设单元体厚度为 1 个单位。沿 方向的正应力称为径向正应力,abcd r用 表示;沿 方向的正应力称为周向正应力或切向正应力,用 表示;r 图 6.10 微元应力分量剪应力用 及 表示。根据剪应力互等定律, 。各应力分量的正负号规rr rr定和直角坐标系中相同,只是 方向对应 方向、 方向对应 方向。图中的应力rxy分量都是正值。径向和周向的体力分量分别用 及 表示。rF将单元体所受的力投影到通过其中心的径向轴上,可建立出单元体径向平衡方程为 02c
36、os2cos)( 2sin2sin)( drFdrdr drrr 第六章 弹塑性平面问题149在上式中,因为 是小量,因此可取 ,并略去高阶d 12cos2sind、d微量后可得 0rrrr F采用同样的方法,可以列出单元体在周向的平衡方程。则可得极坐标系下的平衡方程为(6.5-2)02Frrrrrr类似的还可写出柱坐标系( )下和球坐标系( )下的平衡方程。z, ,r(1)柱坐标系下的平衡微分方方程(6.5-3)020zrzzzr rzr rrrzrr F (2) 球坐标系下的平衡微分方方程(6.5-4) 0cot23sin1 0cot2sin1 Frrr rrrr5.2 几何方程极坐标系下
37、的几何方程,在第三章中巳导出,即(6.5-5)rvurru15.3 本构关系第六章 弹塑性平面问题150极坐标系和直角坐标系都是正交坐标系,因此,在弹灶状态下,极坐标下的本构方程与直角坐标具有同样的形式。只要将下标 分yx,别改写为 即可。于是,r对于平面应力问题(6.5-6)rrrrrGE1()1对于平面应变问题(6.5-7)rrr rrrEG)1(2)(25.4 应变协调方程采用类似推导直角坐标系应变协调方程的方法,不难从式(6.5-5)消除位移分量 ,得出以应变分量表示的极坐标中的应变协调方程,即vu,(6.5-8) rrrr )1()1()2( 222 在直角坐标系中,当体力为常量或不
38、计体力时,平面问题的协调方程式(参见(6.1-7b)式)为0)(2yx注意到 (为不变量),这样在极坐标系中,平面问题应力形式的ryx协调方程式为(6.5-9)0)(2r式中 2为极坐标下的拉普拉斯算子,即(6.5-10)221rr为了得到在极坐标系中,用应力函数 表示的应变协调方程,可直接由),(直角坐标系应变协调方程第六章 弹塑性平面问题1510),()(22yxxyx经坐标变换得到。由式(6.5-1),可得 ryxryxr cos,sinin,co2注意,此处的应力函数 既是 和 的函数,通过坐标变换,也是 和 的函数,x r它对 和 的一阶及二阶导数分别为xy ryryxxcossin
39、i(a) 2222222 222222 222222 cosinsicocosinsin)cs)(sii(c cosinonsi )cs)(sico( siniinincs )(cs( rr rrrryx rrry rrrx将式(a)相加后得 ),(1(),( 222 rrryx 于是得极坐标系下的应变协调方程为(6.5-11)0)1)(1( 2222 rrrr5.5 应力分量第六章 弹塑性平面问题152极坐标系下的应力分量的表达式,也可由坐标转换的方法求得。由图 6.10可见,当把 ox 轴和 oy 轴分别转到 和 的方向,此时, ,则应力分量 ,r0x分别成为 , 。于是,不计体力时,可由
40、式(a)得到极坐标系中的xy,rr、应力分量表达式为(6.5-12) )1(1)()( 1)()( 2020 2020 rryxryxyryxr容易证明,当体力不计时,这些应力分量满足平衡微分方程(6.5-2)式。由以上可知,当体力可以不计时,用极坐标求解平面问题,只须从应变方程(6.5-11)解出应力函数 ,然后由式(6.5-12)求出应力分量,并使其满足位),(r移边界条件和应力边界条件。5.6 轴对称问题如果所研究的问题的物体和外载荷均对称于经过物体中心,且垂直于 平oxy面的轴线,此时,应力和位移均与 无关,仅与 有关,这类问题称为轴对称问r题。因此,轴对称问题只有正应力 和 ,而剪应
41、力因对称性均为零。r(1)应力函数与应力分量根据轴对称问题的情况,应力函数也应与 元关,所以式(6.5-11)可简化为0)1)(22rr将上式展开,并注意到 仅是 的函数,因此偏于数可用常导数代替,得(b)013234drrdr应力表达式(6.5-12)成为(6.5-13)012rrr第六章 弹塑性平面问题153方程式(a)是变系数常微分方程,如令 ,则 , ,根据复合terrlndt1求导法则,则这方程可简化为常系数常微分方程,即 04234dttdt上述方程的解为 tteCtC242321将 代入上式可得rtln 42321lnl)( ArrAr由(6.5-13)式,得应力分量的表达式(6
42、.5-14)02)3ln(1212rrr ArAdr由上式可知,如在坐标原点没有孔,常数 和 必须等于零,否则当12时应力将变为无限大。因此,如在坐标原点没有孔,而且没有体积力,唯一0r可能的应力对称分布是 和 均为常数。对于平面物体,则在平面内必为各方r向均匀受拉或均匀受压状态。如果原点处有孔,则问题有各种解答,这将在下一节中予以讨论。(2)轴对称问题的位移当沿 方向没有约束时,则属平面应力问题。此时,将式(6.5-14)代入式z(6.5-6),并利用式(6.5-5),得(c) 0 )1(2ln)1(2)3(1l322rvru ArAAEru 对上式中的第一式直接积分可得(d)()1(2)(
43、ln)1(2)31(1 32 frArArArEu 再由式(c)的第二式解出 ,并将(d)式代入后,有v第六章 弹塑性平面问题154)(42fErAv积分上式,得(e)()(412rfdfrv将式(d)和式(e)代入式(c)中的第三式,并分离变量,则可得 dffdrff )()()(11此方程左边为 的函数,而右边为 的函数,因此两边必为同一常数 ,于有是rF(f)Fdfdfrf)()(11式(f)中的第一式经简单分析可得其通解为(g)rAf41)(将式(f)中的第二式先对 求导一次,然后再积分求得 )(f(h)cossin65f于是由式(f)的第二式和式(h),可得(i)in)(65AFdf
44、将式(g)、(h)、(i)均代入 和 的表达式(e)和(d)中,则得uv(6.5-15) sinco4sin )1(2)(l)1(2)31(1654265 32ArAEv rrru式中 可由应力边界条件和位移边界条件确定。在应力轴对称时,如果约束条件i也是轴对称的,则位移也应该是轴对称的。即各点无环向位移( ),即0v,仅有径向位移06542(6.5-16)31)(2ArEu对于平面应变问题,以上公式(6.5-15)和(6.5-16)也适用,仅需将式中的和 分别用 和 即可。E)1/(2)/(第六章 弹塑性平面问题1556.6 厚壁圆筒在内外压力下的弹塑性分析工程上一般把圆筒分为厚壁筒和薄壁筒。当外径与内径之比小于 时2.1可按薄壁圆筒进行分析,当大于 1.2 时则按厚壁圆筒进行分析。厚壁圆筒是弹塑性力学问题中最简单的问题之一,即应力和应变只与一个坐标有关,而且在塑性阶段考虑材料的不可压缩性后,可以得到封闭形式的解答,本节讨论的受内外压力作用的厚壁圆筒,属于这类问题。此外还有整球形容器等。6.1 弹性解设图 6.11 所示