资源描述
第一章 绪论1.1 弹性力学的内容1.2 弹性力学的几个基本概念1.3 弹性力学的基本假定 1.1 弹性力学的内容1. 弹性体力学:简称弹性力学,有称弹性理论(Theory of Elasticity),研究弹性体由于受外力、边界约束或温度改变等原因而发生的应力、形变和位移。研究对象:弹性体研究目标:变形等效应,即应力、形变和位移。2. 对弹性力学、材料力学和结构力学作比较 弹性力学的任务和材料力学, 结构力学的任务一样,是分析各种结构物或其构件在弹性阶段的应力和位移, 校核它们是否具有所需的强度和刚度, 并寻求或改进它们的计算方法. (1)研究对象: 材料力学主要研究杆件在拉压、剪切、弯曲、扭转作用下的应力、形变和位移; 结构力学研究杆系结构,如桁架、钢架或两者混合的构架等; 弹性力学研究各种形状的弹性体,除杆件外(对杆件进行进一步的、较精确的分析),还研究平面体、空间体,板和壳等。(2)研究方法: 弹性力学与材料力学有相似,又有一 定区别。 弹性力学:在弹性体区域内必须严格考虑静力学、几何学和物理学三方面条件,在边界上严格考虑受力条件或约束条件,由此建立微分方程和边界条件进行求解,得出精确解答。材料力学:虽然也考虑这几个方面的的条件,但不是十分严格。 一般地说, 由于材料力学建立的是近似理论, 因此得出的是近似的解答。但对于细长的杆件结构而言, 材料力学力解答的精度是足够的, 符合工程的要求。 弹性力学:梁的深度并不远小于梁的跨度,而是同等大小的,那么,横截面的正应力并不按直线分布,而是按曲线变化的。qq zI yxM )( )53(4)( 22 hyhyqI yxM z 例如:材料力学:研究直梁在横向载荷作用下的平面弯曲,引用了平面假设,结果:横截面上的正应力按直线分布。这时,材料力学中给出的最大正应力将具有很大的误差。 结构力学:研究杆系结构,弹性力学通常并不研究杆件系统,但在20世纪50年代中叶发展起来的有限单元法中(基于弹性力学的理论),把连续体划分成有限大小的单元构件,然后用结构力学里的位移法、力法或混合法求解,更加显示了弹性力学与结构力学结合综和应用的良好效果。 弹性力学在土木、水利、机械、航空等工程学科中占有重要的地位。许多非杆件形状的结构必须用弹性力学方法进行分析。例如,大坝,桥梁等。 x z yo 1.2 弹性力学中的几个基本概念 弹性力学的基本概念: 外力、应力、形变和位移1. 外力:体积力和表面力,简称体力和面力体力:分布在物体体积内的力,例如重力和惯性力。V Pf FfVFlim0V fx fyfz f : 极限矢量,即物体在P点所受体力的集度。方向就是F的极限方向。fx , fy , fz:体力分量, 沿坐标正方向为正,沿坐标负方向为负。量纲:N/m3=kgm/s2m3=kg/m2s2即:L-2MT-2 fx , fy , fz:体力分量。x z yo f fSFlim0V S P面力:分布在物体表面的力,例如流体压力和接触力。Ffyfzfx量纲:N/m2=kgm/s2m2=kg/ms2即:L-1MT-2f : 极限矢量,即物体在P点所受面力的集度。方向就是F的极限方向。沿坐标正方向为正,沿坐标负方向为负。符号规定: 内力:发生在物体内部的力,即物体本身不同部分之间相互作用的力。x z yo PA pF 2. 应力:单位截面面积的内力.pAF lim0Vp: 极限矢量,即物体在截面mn上的、在P点的应力。方向就是F的极限方向。量纲:N/m2=kgm/s2m2=kg/ms2 即:L-1MT-2应力分量:, A BC zyzxz yzyx yxy xzxy yzyx zy zzxxyxzxP yx zOPA=x, PB=y , PC=z x, y, z, xy, xz, yx, yz, zx, zy,正面:截面上的外法线沿坐标轴的正方向正面上的应力以沿坐标轴的正方向为正,沿坐标轴的负方向为负。负面:截面上的外法线沿坐标轴的负方向负面上的应力以沿坐标轴的负方向为正,沿坐标轴的正方向为负。正应力符号规定与材力同,切应力与材力不相同。符号规定:(不考虑位置, 把应力当作均匀应力) A BC zyzxz yz yx yxy xzxby yzyx zy zzxxyxzx aP yxzo 连接前后两面中心的直线ab作为矩轴,列出力矩平衡方程,得02222 zxyyxz yzzy 得: zyyz 同理可得:yxxy zxxz 切应力互等定理:作用在两个互相垂直的面上并且垂直于该两面角线的切应力是互等的(大小相等,正符号也相同)。 可以证明,已知x, y, z, yz, zx, xy, 就可求得该点任意截面上的, .因此,此六个应力分量可以完全确定该点的应力状态。 zyzxz yzyx yxy xzxy yzyx zy zzxxyxzxP yx zO A BC A BC zyzxz yz yx yxy xzxy yzyx zy zzxxyxzxP yxzO用各部分的长度和角度来表示。PA=x, PB=y , PC=z线应变:单位长度的伸缩或相对伸缩,亦称正应变. 用 表示切应变:各线段之间的直角的改变.用 表示3. 形变:就是形状的改变。 A BC zyzxz yz yx yxy xzxy yzyx zy zzxxyxzxP yxzO x: x方向的线段PA的线应变。xy: y与x两方向的线段PB与PC之间的直角的改变。 : 伸长为正,缩短为负。量纲:1符号规定: : 直角变小为正,变大为负。 可以证明,已知x, y, z, yz, zx, xy, 就可求得经过该点任一线段上的线应变 .也可以求得经过该点任意两个线段之间的角度的改变。因此,此六个形变分量可以完全确定该点的形变状态。 4. 位移:就是位置的移动。任意一点的位移用它在x,y,z三轴上的投影u,v,w来表示.量纲:L符号规定:沿坐标轴正方向为正,沿坐标轴负方向为负, 一般而论,弹性体内任意一点的体力分量、面力分量、应力分量、形变分量和位移分量都随该点的位置而变,因而都是位置坐标的函数。 1.3 弹性力学中的基本假设 在弹性力学的问题里,通常是已知物体的边界(形状和大小), 物体的弹性常数, 物体所受的体力,物体边界上的约束情况或面力, 而应力分量、形变分量和位移分量则是需要求解的未知量.一. 研究方法1.考虑静力学、几何学和物理学三方面条件,分别建立三套方程。建立微分方程:根据微分体的平衡条件 ;建立几何方程:根据微分线段上形变与位移之间的 几何关系;建立物理方程:根据应力与形变之间的物理关系 。 2.在弹性体的边界上,建立边界条件。应力边界条件:在给定面力的边界上,根据边界上 的微分体的平衡条件;位移边界条件:在给定的约束边界上,根据边界上 的约束条件。 求解弹性力学问题,即在边界条件下根据平衡微分方程、几何方程、物理方程求解应力分量、形变分量和位移分量。 为使问题求解成为可能,通常必须按照所研究的物体性质,以及求解问题的范围,略去一些影响很小的次要因素,作出若干基本假定。二. 弹性力学的基本假定(3)均匀性 假定物体是均匀的.(1)连续性 假定物体是连续的.(4)各向同性 假定物体是各向同性的.符合以上四个假定的物体,就成为理想弹性体.(2)完全弹性 假定物体是完全弹性的.形变与引起 变的应力成正比,即两者成线性关系. (5)小变形假定 假定位移和形变是微小的.它包含两个含义: 假定应变分量 1.例如:普通梁中的正应变 10-3 1,切应变 1; 假定物体的位移物体尺寸.例如:梁中挠度 梁的高度这样,在建立平衡微分方程时,可以用变形前的尺寸代替变形后的尺寸,从而使方程大为简化;在建立几何方程时,由于 1,可以在同一方程中只保留形变成分的一次幂,而略去二次幂及更高次幂,从而使几何方程成为线性方程。 例如:对于微小转角a,1211cos 2 aa aaa 3!311sin aaaa 331tan对于微小正应变,xxxxx 111 1 32 这样,弹性力学里的几何方程和微分方程都简化为线性方程,弹性力学问题都化为线性问题,从而可以应用叠加原理。 第二章 平面问题的基本理论2.1 平面应力问题与平面应变问题2.2 平衡微分方程2.3 平面问题中一点的应力状态2.4 几何方程 刚体位移2.5 物理方程 2.6 边界条件2.7 圣维南原理2.8 按位移求解平面问题2.9 按应力求解平面问题 相容方程2.10 常体力情况下的简化 应力函数 2.1 平面应力问题与平面应变问题 如果弹性体具有某种特殊的形状,并且承受的是某些特殊的外力和约束,就可以把空间问题简化为近似的平面问题。一.第一种平面问题平面应力问题 xyo zyd/2 d/2这类问题的条件是:弹性体是等厚度(d)的薄板,体力、面力和约束都只有xy平面的量 (fx , fy , fx , fy , u, v ),都不沿z向变化;并且面力和约束只作用于板边,在板面( )上没有任何面力和约束的作用。2z 因板很薄,外力不沿厚度变化,应力沿板厚连续,有由切应力互等定理:0 xz 0yz0z 0zx 0zy只剩下平行于xy面的三个平面应力分量,即 x, y, xy= yx所以这种问题称为平面应力问题。xyo zy d/2 d/21.设薄板的厚度为d, xy 为中面,z轴垂直于xy面.因为板面上 2z 不受力, 所以0)( 2 d zz 0)( 2 zzx 0)( 2 d zzy2.由于物体形状和外力、约束沿z向均不变化, 故 x, y, xy 只是x,y的函数, x, y, xy 也只是x,y的函数,但位移与z有关。 二.第二种平面问题平面应变问题oy x这类问题的条件是:弹性体为常截面的很长的柱体,体力、面力和约束条件与平面应力问题相似,只有xy平面的体力fx , fy ;面力fx , fy 和约束 u, v 的作用,且都不沿z向变化。 2.2 平衡微分方程 在弹性力学中分析问题,要考虑静力学、几何学和物理学三方面条件,分别建立三套方程。 首先考虑平面问题的静力学方面,建立微分体的平衡微分体方程应力分量与体力分量之间的关系式。z yd/2 d/2 oy xxyo从图示薄板或柱形体中,取出一个微小的正六面体,边长为dx, dy, 在z方向的尺寸取为1个单位尺寸。xyo dx dy 300200000 )(!31)(!21)()()( xxxfxxxfxxxfxfxf 222 d21d xx!xx xxx 一般而论 , 应力分量是位置坐标x和y的函数, 因此, 作用于左右两对面或上下两对面的应力分量不完全相同, 有微小的差。o xy x xxxx d 略去二阶及二阶以上的微量后得:例:设作用于左面的正应力为x,则右面的正应力由于 x 坐标的改变而改变,可由泰勒展开得:xxxx d 若x为常量, 则 , 左右两面都是x,即为均匀应力。0 xx泰勒展开式 o xy x xxxx d 同理,设左面的切应力为xy,则右面的切应力为xy xdxxyxy yyx yy yy d ydy yxyx C fxfyxxyxyx d 设上面的正应力及切应力为x, xy,则下面的正应力其切应力为,dyyyy yyyxyx d 因六面体是微小的, 所以, 各面的应力可认为是均匀分布, 作用在对应面中心. 所受体力也可认为是均匀分布, 作用在对应面中心。 02d1d2d1dd 2d1d2d1dd yxyxyy xyxyxx yxyxyx xyxyxy o C xy yyxxyx xxxx d yy yy d xxxyxy d yy yxyx dfxfyyyxx yxyxxyxy d21d21 yxxy 首先,以过中心C 并平行于z轴,列出0 CM将上式除以dxdy, 得令dx,dy 趋近于零,得这正是切应力互等定理。 o C xy yyxxy x xxxx d yy yy d xxxyxy d yy yxyx dfxfy01dd1d1dd 1d1dd yxfx xyy yyxx xyxyxyx xxx 0 xyxx fyx 其次,以x轴为投影轴,列出0 xF将上式除以dxdy, 得同样,以y轴为投影轴,列出 可得一个相似的微分方程0 yF0 yxyy fxy 于是得出应力分量与体力分量之间的关系式平面问题中的平衡微分方程。 这2个微分方程中包含3个未知函数x, y, xy=yx ,因此,决定应力分量的问题是超静定问题,必须考虑几何方程和物理学方面的条件,才能解决问题。 对于平面应变问题, 微分体一般还有作用于前后两面的正应力 z, 但不影响上述方程的建立, 上述方程对于两种平面问题同样适用。 00yxyy xyxx fxy fyx 2.3 平面问题中一点的应力状态O x y yyxxyx PB AnnO xy yyxxyx yyx xxyP 应力状态就是指一点处所有斜截面上的应力的集合。假定已知任意点P处坐标面的应力分量x, y, xy=yx ,求经过该点且平行于z轴的任意斜截面上的应力。 myn lxn ),cos(),cos( ssmlf smslsp xxyxx 02 ddddd xyxx mlp xyyy lmp py pxpO xy yyxxyx nPB A用n代表斜截面AB的外法线方向,其方向余弦为设AB=ds, 则PA=lds, PB=mds, SPAB =ldsmds/2设垂直于平面的尺寸为1。 由 得0 xF其中 fx 为x方向得体力分量。将上式除以ds, 然后命ds 趋于0(AB0)得同理由 得0 yF一.求任意斜截面上的正应力 n 和 切应力 n yxn mplp xyyxn lmml 222 xyn mplp xyxyn mllm )()( 22 nnpy pxpO xy yyxxyx nPB A令斜截面得正应力为n, 切应力为n.由px, py 投影得xyxx mlp xyyy lmp 可见,已知点P处的应力分量x, y, xy=yx ,就可求得经过该点的任意斜截面上的正应力 n 和 切应力 n 。 lpx O x y yyxxyx yyx xxyP 121 2 a1nn O xy yyxxyx PB A lml xyx mlm xyy )1( xy xlm )2(yxylm mpy 0)()( 22 xyyxyx 二.求主应力及主应力的方位应力主向应力主面上 = 0, = pxyxx mlp xyyy lmp 投影得代入得py pxp由上两式分别解出 m / l , 得于是,有解得 2221 22 xyyxyx yx 21 O xy yyxxyx yyx xxyP 121 2 a1111 1111 cos90coscossintan lm)( o a aaaa xy x a 11tan 222 2222 cos90coscossintan lm)( o a aaaa yxy a 22tan易得下面求主应力方向即得即得设1与x轴的夹角为a1设2与x轴的夹角为a2 )1(xy xlm )2(yxylm )( 12 xy xxy a 12tan 1tantan 21 aa 1max 2min 2 21 max 2 21min O xy yyxxyx yyx xxyP 121 2 a1yx 21由得xy x a 11tan yxy a 22tan于是有就是说, 1, 2 的方向互相垂直。从材料力学知识我们知道与应力主向成450的斜面上。 yyuu d2.4 几何方程 刚体位移xyO PB Au xxuu d xxvv dyyvv d P AB axux uxxuu( PA PAAPPAPAAP xx d )d)( yx PBPA dd yvy 同理PB的线应变:PA的线应变:一.几何方程:任一点的微分线段上的形变分量与 位移 分量之间的关系式。v设 同理PB的转角:yuPA与PB之间的转角:yuxv xy a yyuu d xyO PB Au xxuu d xxvv dyyvv d P AB avxvx v)xxvv( ddsinaaPA的转角:yuxv yv xu xyyx ,几何方程:上列几何方程对两种平面问题同样适用。 二.形变与位移之间的关系1. 如果物体的位移确定,则形变完全确定。从物理概念: 当物理变形后各点的位置完全确定, 任一微分线段上的形变(伸缩、转角等)也就完全确定了.从数学概念: 当位移函数确定时,其导数也就确定了。2. 当物体的形变分量确定时,位移分量不完全确定。从物理概念: 在物体内形变不变的条件下, 物体还可以做刚体运动平动和转动, 即还有刚体运动的人任意性. 从数学概念: 由形变分量求位移分量是一个积分的过程,在常微分中,会出现一个任意常数;而在偏微分中,要出现一个与积分变量无关的任意函数。这些任意函数是未定项,这些未定项正是刚体平移和刚体转动量。若假设 求出相应的位移分量。 xyyx 0代入几何方程: . yuxv , yv , xu 000 将前二式对x及y积分,得 , )( , )( 21 xfvyfu F1 及 f2 为任意函数。代入几何方程中的第三式,得xxfyyf d )(dd )(d 21 xxfyyf d )(dd )(d 21 方程左边是y的函数,只随y而变;而右边是x的函数,只随x而变。因此,只可能两边都等于同一常数。于是得 xxfyyf d )(d,d )(d 21积分得 , )( , )( 0201 xvxfyuyf 其中u0及v0为任意常数。代入 得 , )( , )( 21 xfvyfu , , 00 xvvyuu 这就是“形变为零”时的位移,也就是所谓“与形变无关的位移”,因此必然是刚体位移。下面根据平面运动的原理加以证明。 u0及v0分别为物体沿x轴及y轴方向的刚体位移,而为物体绕z轴得刚体转动。 Pxyx yOz ay x , , 00 xvvyuu 当只有u0不为零时,物体内任一点位移分量 .物体的所有各点只沿x方向移动同样距离u0,所以u0代表物体沿x方向的刚体位移。0 , 0 vuu坐标为(x,y)的任一点P沿y方向移动x, 沿x负方向移动y, 合成位移为 xvyu , 222222 yxxyvu同样, v0代表物体沿y方向的刚体位移。当只有不为零时, 物体内任一点位移分量 Pxyx yOz ay x可见, 合成位移的方向与径向线段OP垂直,也就是沿着切向. 因OP线上所有点移动方向都沿着切线, 且移动的距离为, 可见代表物体绕z轴的刚体转动。a tan/tan y/xxy 既然物体在形变为零时可以有刚体位移,那么,当物体发生一定形变时,由于约束条件不同,可能有不同的刚体位移,为了完全确定位移,就必须有适当的刚体约束条件。 2.5 物理方程物理方程:应力分量和形变分量之间的物理关系式.在理想弹性体(满足连续性,完全弹性,均匀性和各向同性)中,物理方程就是材料力学中学过的胡克定律:物理方程有两种形式:1. = f () 此式是用应力表示应变,其中应力取为基本未知数,用于按应力求解。2. = f () 此式是用应变表示应力,其中应变取为基本未知数,用于按位移求解。 胡克定律的一般形式: ,1 ,1 ,1 )(1 )(1 )(1 xyxyxyzxxyyz yxzz zxyy zyxx GGGEEE E是弹性模量,G是切变模量,又称刚度模量,称为泊松系数,或泊松比。 12 EG一.平面应力问题的物理方程将 代入上式得独立的物理方程0 zyzxz ,)1(211 xyxy xyy yxx EEE 另外:)( yxz E 因z可由x , y求出, 故不作为独立的未知函数。 二.平面应变问题的物理方程将 代入上式得独立的物理方程0 zyzxz ,)1(2 11 11 22 xyxy xyy yxx EEE 另外:)( yxz E 因z可由x , y求出, 故不作为独立的未知函数。与平面应力问题的物理方程对比,只需将E 换为 , 换为21 E 1 对于两类平面问题,三套方程除了物理方程中的系数须变换外, 其他平衡方程和几何方程是完全相同的. 三套方程中包含8个未知函数:x, y, xy=yx, x, y, xy及u, v. 还需考虑边界条件, 才能求出这些未知函数. 2.6 边界条件边界条件表示在边界上位移与约束,或应力与面力之间的关系式。它分为位移边界条件、应力边界条件和混合边界条件。一.位移边界条件设在部分边界上给定了约束位移分量u(s)和v(s), 则对于边界上的每一点,位移函数u,v应满足条件 )()(),()( svv suu ss (在su上)其中(u)s和(v)s是位移的边界值,u(s)和v(s)在边界上是坐标的已知函数。位移边界条件 注意1.上式要求在s上任一点位移分量必须等于对应的约束位移分量。)()(),()( svv suu ss (在su上)2.上式是函数方程,而不是简单的代数方程或数值方程。位移边界条件实质上是变形连续条件在约束边界上的表达式。 设n为斜截面的外法线方向,其方向余弦二.应力边界条件设在s部分边界上给定了面力分量fx(s)和fy(s), 则可以由边界上任一点微分体的平衡条件,导出应力与面力之间的关系式。 在边界上任一点P取出一个微分体,斜面AB就是边界面, x, y, xy为应力分量边界值。o xy yyxxyx PB Afxfy1.边界为斜截面时nmyn lxn ),cos(),cos(设 AB=ds , z 方向厚度为1 由平衡条件,得出微分体的应力分量与边界面上的面力之间的关系: ).()( ),()( sflm sfml ysxyy xsyxx (在s 上)其中 在边界上是坐标的已知函数,l, m 是边界面外法线的方向余弦。fx(s)和fy(s), o xy yyxxyx PB Afxfy ssmlf smslsf xxyxx 02 1dd1d1d1d n除以ds, 并令ds0,得 sfml xsyxx ),()( 同理:).()( sflm ysxyy 于是,得到应力边界条件 3. 在导出应力边界条件时, 只考虑到面力(一阶微量), 不需考虑二阶微量体力。4. 应力边界条件是边界点上微分体的平衡条件,也属于静力边界条件。).()( ),()( sflm sfml ysxyy xsyxx (在s 上)注意1.应力边界条件表示边界s上任一点的应力和面力之间的关系。也是函数方程,在s上每一点都应满足。2.上式中的面力、应力都有不同的正负符号规定,且分别作用于通过边界点的不同面上。 2.边界为坐标面时若 x=a 为正 x 面,则有).()(),()( yf yf yaxxyxaxx 若 x=b 为负 x 面,则有).()(),()( yf yf ybxxyxbxx o xy yyxxyx PB A fxfyao xy yyx xyxPB Afx fy b正负 x 面上的面力分量一般为随 y 而变化的函数。l =1, m=0l = 1, m=0 ).()( ),()( sflm sfml ysxyy xsyxx (在s 上) 3.应力边界条件的两种表达方式(1) 在边界点取出一个微分体,考虑其平衡条件,得出).()( ),()( sflm sfml ysxyy xsyxx (在s 上)).()(),()( yf yf yaxxyxaxx ).()(),()( yf yf ybxxyxbxx (2) 在同一边界面上,应力分量的边界值就等于对应的面力分量。 应力分量的绝对值等于对应的面力分量的绝对值,面力分量的方向就是应力分量的方向。 即数值相同,方向一致。 例如:若边界面 y=c, d 分别为正、负坐标面).()(),()( xf xf xdyyxydyy 在斜截面上:).()(),()( xf xf xcyyxycyy ).()(),()( sfp sfp ysyxsx px, py 为斜截面应力o xy y yx xyxPfx fy dy yxfyo xyyx xyxPfx c yxyy xyxx mlp xyyy lmp o xy yyxxyx P fxfypxp y 三.混合边界条件物体的一部分边界具有已知位移,因而具有位移边界条件,如)()(),()( svv suu ss (在su上)另一部分边界则具有已知面力,因而具有应力边界条件).()( ),()( sflm sfml ysxyy xsyxx (在s 上)在同一边界上还可能出现混合边界条件,即两个边界条件中一个是位移边界条件, 另一个则是应力边界条件.o xy x方向0)( uu sy方向0)( ysxy f x方向0)( vv sy方向0)( xsx fo xy 2.7 圣维南原理及其应用求解弹性力学问题时,应力分量、形变分量和位移分量必须满足三套方程,还必须满足边界条件,但要使边界条件得到完全满足很困难。圣维南原理为简化局部边界的应力边界条件提供了有效的方法。圣维南原理:如果把物体的一小部分边界上的面力, 变换为分布不同但静力等效的面力(主矢量相同,对于同一点的主矩也相同), 那么, 近处的应力分布将有明显的改变 , 但是远处所受的影响可以不计。 1. 圣维南原理只能应用于一小部分边界上,又称为局部边界,小边界或次要边界。一. 圣维南原理应用的条件 所谓“近处”, 根据经验, 一般地讲大约是变换面力的边界的12倍范围内, 此范围之外可认为是“远处”。如果将面力的等效变换范围应用到大边界(又称为主要边界)上, 则必然使整个的应力状态都改变了。因此, 不适用圣维南原理。 F F/2F/2F F Fq2.小边界的面力变换为静力等效的面力.3.经变换后,只对近处的应力分布有明显的影响,但远处的应力几乎不受影响。F F/2F/2 F FF F/2F/2F/2 F/2F/2F/2(a)(b)(c)例如:如将一端或两端的F变换为静力等效的力, 如图(b), (c), (d).则只有虚线划出的部分应力分布有显著改变, 其余部分所受影响可不计。(d) F/A F/A图(d)所示情况,由于面力连续均匀分布,边界条件简单,应力很容易求解并且解答很简单。而其他三种情况,由于面力不连续分布,甚至不知其分布方式,应力难以求解。根据圣维南原理,可将(d)的应力解答应用于其他三种情况。 应用圣维南原理的条件是满足静力等效。即使物体一小部分边界上的位移边界条件不能满足时,仍可以应用圣维南原理。F/A F/AF(e)(d)图(e)右端是固定端,有位移边界条件(u)s = u = 0 和(v) s = v = 0, 把(d)的解答应用于这一情况时,位移 边界条件不能满足, 但右端的面力静力等效于过形心的力F(与左边的力F平衡), 满足圣维南原理的条件, (d)的解答仍可应用于这一情况时, 只是在靠近两端处有显著的误差, 而在较远处误差可不计。 如果物体一小部分边界上的面力是一个平衡力系(主矢量和主矩都等于零),那么,这个面力就只会使近处产生显著的应力,而远处的应力可以不计。这是因为主矢量和主矩都等于零的面力,与无面力状态是等效的,只在近处产生显著的应力。例如: F FF F4.圣维南原理还可以推广到下列情形 xyh/2h/2 ll O 在应力边界条件上应用圣维南原理, 就是在边界上, 将精确的应力边界条件代之以主矢相同, 对同一点的主矩也相同的静力等效条件。二.在局部边界上应用圣维南原理例如,厚度 d=1 的梁,hl, 即左右端是小边界.严格的边界条件要求 yf yf ylxxyxlxx ,x xy fx fy xy ydy x fxfy此式要求在边界 x=l 上的每一点(每一y值),应力分量与对应的面力分量必须处处相等。 严格的边界条件要求 yf yf ylxxyxlxx , xyh/2h/2 ll Ox xy fx fy xy ydy x fxfy这种严格的边界条件是很难满足的。 但 h h, d=1 ) qF 30oO xyb/2 hgy b/2 ( h b, d=1 )(a) (b)解:对(a)问题,在主要边界 y=h/2,应精确满足下列边界条件 12,0,2 0,2 q hy lxq hy yxy yxy 2)(lxq q 1F FsM O xy l h/2h/2( l h, d=1 )(a)在小边界(次要边界) x=0, 应用圣维南原理, 列出三个积分近似边界条件, 当板厚d=1时, Fy Myy Fy sxh/h xy xh/h x xh/h x ,d ,d ,d022/ 022/ 022/ 在小边界 x=l 处,当平衡微分方程和其它各边界都已满足条件下,三个积分的边界条件必然满足,可以不必校核。 qF 30oO x yb/2 hgy b/2 ( h b, d=1 )(b)对(b)问题,在主要边界 y=0, b, 应精确满足下列边界条件q bx gy x yxy yxx ,0, 0,0在小边界 y=0, 列出三个积分近似边界条件, 当板厚d=1时, Fx Fxx Fx yh/h yx yh/h y yh/h y ,2d ,43d ,23d022/ 022/ 022/ 注意,在列力矩条件时,两边均是对原点O 的力矩来计算的。对于 y=h 的小边界条件可以不必校核。 F O xy l h/2h/2( l h, d=1 ) A例2:厚度d=1的悬臂梁,在自由端受集中力F的作用。已求得其位移的解答是.3262 82662 3232 22332 EIFlEIxFlEIFxEIFxyv y ,IGFhEIFlIGFyEIFyEIyFxu 试检查此组位移是否是该问题的解答。解:此组位移若为此问题的解答,则应满足下列条件 .021211 ,021211 222222 222222 yxfyxuxvyvE fyxvyuxuE 1. 在区域内,满足用位移表示的平衡微分方程 . fyuxvlxuyvmE fxvyumyvxulE ys xs211 ,211 22 )()(),()( svv suu ss 在Su上2. 在所有受面力的边界s上,满足应力边界条件。3. 在su满足位移边界条件其中在小边界上可以应用圣维南原理,即用三个积分的边界条件来代替。本题只需校核在 边界x=l 的刚体约束条件A点( x=l 及y=0 ),0, xuvu . fxuyvE fyvxuE ys xs 221 ,1 例3:试考虑下列平面问题的应变分量是否存在,(a) x=Axy, y=By3, xy=CDy3(b) x=Ay2, y=Bx2y, xy=Cxy(c) x=y=0, xy=Cxy解:应变分量存在的必要条件是满足形变协调条件(相容方程)即 yxxy xyyx 22222(a) 相容(b) 须满足 B=0,2A=C(c) 不相容 只有C=0, x = y= xy= 0, 例4:在无体力的情况下,试考虑下列应力分量是否可能在弹性体中存在。(a) x=Ax + By, y= Cx + Dy, xy= Ex + Fy ;(b) x=A( x2 + y2 ), y=B ( x2 + y2 ), xy=Cxy解:弹性体中的应力,在单连体中必须满足 00yxyy xxyx fxy fyx 02 yx ).()( ),()( sflm sfml ysxyy xsyxx 在s上(a) 此组应力满足相容方程,为满足平衡微分方程,必须 A=F, D=E,此外,还须满足应力边界条件。 (b) 为满足相容方程,其系数必须满足 A+B=0为满足平衡微分方程,其系数必须满足 A =B =C/2上两式是矛盾的,故此组应力不存在。(b) x=A( x2 + y2 ), y=B ( x2 + y2 ), xy=Cxy例5:若f (x,y) 是平面调和函数,即满足拉普拉斯方程 0 2 f试证明函数f , xf , yf , (x2+y2)f 都满足重调和方程,因而都可以作为应力函数 使用。证明:0224 ff 0444 xffxxf 0444 yffyyf 0224422224 yxffyxfyx上述函数作为应力函数,均能满足相容方程(重调和方程)0 4 例6:图示梁受到均布载荷的作用,试用下列应力表达式求解其应力。 xChqxy CyCyhq yyxhq xyyx 13 2 2133 32362 46 q xy l h/2h/2( l h, d=1 )202qh ql 202 22 qhql O 解: 00yxyy xxyx fxy fyx 02 yx ).()( ),()( sflm sfml ysxyy xsyxx 在s上 本题是按应力求解,因而,应力分量必须满足 xChqxy CyCyhq yyxhqxyyx 13 2 2133 32362 46 将应力分量代入平衡微分方程和相容方程,两者都能满足。再校核边界条件,在主要边界上hqC Chhqx hy yx 23 ,0460,2 1 123 得即 q xy l h/2h/2( l h, d=1 )202qh ql 202 22 qhql O 2 ,0282,2 2133qC ChChhq q hy 2y 得即 xChqxy CyCyhq yyxhqxyyx 13 2 2133 32362 46 hqC 231 代入后满足。将210,2 CC hy ,y 将C1,C2代入应力分量,得 1423 22321 232 22 33223 hyhqx hyhyq yxhqyxyyx q xy l h/2h/2( l h, d=1 )202qh ql 202 22 qhql O 1423 22321 232 22 33223 hyhqx hyhyq yxhqyxyyx 再将应力表达式代入次要边界条件: 20d 0d4 ;00 202/ 2/ 02/ 2/33 qhyy y ,hyq ,x xhh x xhh xx xy 其主矩为其主矢量为 202d 046 d;1423 222/ 2/ 323 02/ 2/ 33 qhqlyy ,yylhq qly hyqhql ,lx lxhh xx xhh xyxy 其主矩为其主矢量为其主矢量为可见,在次要边界上的积分边界条件均能满足。 q xy l h/2h/2( l h, d=1 )202qh ql 202 22 qhql O 例7:材料力学中,当矩形截面梁(厚度d=1)受任意横向载荷q(x)作用而弯曲时,弯曲正应力公式为 yIxMx q(x) xy l h/2h/2 ( l h, d=1 )O1.试由平衡微分方程(不计体力)导出切应力xy和挤压应力x的公式(提示:注意 xqxxF xFxxM ss )(dd),(dd 积分后得出的任意函数,可由梁的上下边界条件来确定.)解:不计体力, 将 代入平衡微分方程第一式 yIxMx y x yxx 0 得 I yxFIyxxMxy sxyx dd q(x) xy l h/2h/2( l h, d=1 )O两边对y积分,得 xfIyxFsyx 122 再由上下边界条件 得 02 h/yyx IhxFxf s8 21 其中28 22 yhS 将 代入平衡微分方程第二式 I SxFsyx x y xyy 0 I yxFIyxxMxy sxyx dd I SxFsyx 代入上式, 得得 28281dd 2222 yhIxqyhIxxFxy sxyy 28281dd 2222 yhIxqyhIxxFxy sxyy q(x) xy l h/2h/2 ( l h, d=1 )O I SxFsyx 28 22 yhS 两边对y积分,得 xfyyhIxq 2y 32 618再由上下边界条件 02 h/yy得 224 32 xqIhxqxf qh/yy 2由同样得 22 xqxf 代入 xfyyhIxq 2y 32 618得 33333 223216824 hyhyxqyhhIxqy yIxMx I SxFsyx 28 22 yhS 3322321 hyhyxqy 上述解答已满足平衡微分方程及y=h/2的边界条件, 但一般不满足相容方程, 且尚未校核左右端的小边界条件。2.当q为常数时, 试检验应力分量是否满足相容方程?试在 x中加一项对平衡没有影响的函数f (y), 再由相容方程确定f (y) , 并校核梁的左右边界条件。xy l h/2h/2( l h, d=1 )O q 若q=常数,则 2222 222 lxlxqlxqxqlxM xy l h/2h/2( l h, d=1 )O q ylxlxhqlyIxMx 22326于是 3322321 hyhyqy代入相容方程, 02432 yqqyx 为满足相容方程,令 yfylxlxhqlyIxMx 22336此时, 3322321 hyhyqy yfylxlxhqlyIxMx 22336 I SxF syx 和仍满足平衡微分方程,再代入相容方程。 0dd24 2232 y yfyqqyx xy l h/2h/2( l h, d=1 )O q 3322321 hyhyqy yfylxlxhqlyIxMx 22336积分得 04 33 BAyyhqyf由 x=l 次要边界条件 0d2/ 2/ ylxhh x得 B=0 ; 0d02/ 2/ yxhh x满足。 0d,02/ 2/ yylxhh x得hqA 53由此得 yhqyhqylxlxhqlyIxMx 5346 332233 经检验,在小边界 x=0,l 上剪力边界条件亦满足。 第三章 平面问题的直角坐标解答3.1 逆解法和半逆解法 多项式解答3.2 矩形梁的纯弯曲3.3 位移分量的求出3.4 简支梁受均布载荷3.5 楔形体受重力和液体压力 3.1 逆解法和半逆解法 多项式解答在常体力情况下,按应力求解平面问题,可归纳为求解一个应力函数,它必须满足1.在区域内的相容方程2.在边界上的应力边界条件(假设全部为应力边界条件)3.在多连体中,还须满足位移单值条件。).()( ),()( sflm sfml ysxyy xsyxx 在s上 . yx yfx xfy xyyyxx 22222 ,求出应力函数 后,便可求出应力分量.然后再求应变分量和位移分量。02 4422444 xyxy 02 4422444 xyxy 由于相容方程 是偏微分方程,它的通解不能写成有限项数的形式,一般不能直接求解问题。只能采取逆解法和半逆解法。所谓逆解法,就是(1)先设定满足 的应力函数;02 4422444 xyxy (2)根据 求出应力分量; . yx yfx xfy xyyyxx 22222 , (3)在给定的边界形状下,根据应力边界条件,由应力反推出相应的面力,即.)( )( sxyyy syxxx lmf ,mlf 反过来得知所选取的应力函数可以解决的问题。(可解决的正是上述面力对应的问题)一.逆解法 下面用逆解法求解几个简单问题的解答。假定体力可忽略不计( fx = fy = 0 ),应力函数取为多项式。1.取应力函数为一次式 = a + bx + cy应力函数 满足相容方程02 4422444 xyxy 由 得应力分量 . yx yfx xfy xyyyxx 22222 , . xyyx 0,0,0 不论弹性体为何形状,也不论坐标轴如何选择,由应力边界条件 总是得出.)( )( sxyyy syxxx lmf ,mlf f f yx .0,0 一次式 = a + bx + cy对应无体力,无面力,无应力的状态。把应力函数加上一个线性函数,不影响应力。 2.取应力函数为二次式 = ax2 + bxy + cy2应力函数 满足相容方程02 4422444 xyxy 现分别考察每一项所能解决的问题。对应 = ax2,应力分量是. a yxxyyx 0,2,0 (a)2a xyO2a如图矩形板和坐标轴,当板内应力为x = 0, y = 2a, xy=yx = 0, 由应力边界条件可知,左右两边没有面力,上下两边有均布面力2a。可见,应力函数 = ax2 能解决矩形板在 y 方向受均布力的问题。 b (b)b xyO bb (c)2c xyO 2c 如图矩形板和坐标轴,当板内应力为x = 0, y = 0, xy=yx =b, 由应力边界条件可知,左右上下两边分别有与面相切的面力 b。 可见,应力函数 = bxy 能解决矩形板受均布剪力的问题。对应 = bxy,应力分量是b. yxxyyx ,0,0对应 = cy2,应力分量是. c yxxyyx 0,0, 应力函数 = cy2 能解决矩形板在 x方向受均布力的问题。 = ax2 + bxy + cy2 表示常量的正应力和切应力。 4.如果取应力函数为四次或四次以上的多项式,则其中的系数必须满足一定的条件。应力函数 满足相容方程02 4422444 xyxy 对应 = ay3,应力分量是. ay yxxyyx 0,0,6 O y x对于图示矩形板和坐标轴当 时,上下两边没有面力;左右两边没有 y 方向面力. ay yxxyyx 0,0,6 ,只有按直线变化的水平面力,而每一边的水平面力合成为一个力偶。 可见,应力函数 = ay3 能解决矩形梁纯弯曲问题。3.取应力函数为三次式 = ay3 Oy xh/2h/2 llh5.例题例1:图示矩形长梁, lh, 试考察应力函数 能解决什么样的受力问题。 223 432 yhxyhF 解:按逆解法求解1.将 代入相容方程,满足相容方程 02 4422444 xyxy 2.将 代入 得应力分量 . yx yfx xfy xyyyxx 22222 , 222222322 4123,012 hyhFyx x ,hFxyy xyyx 3.由边界形状和应力分量反推边界上的面力 222222322 4123,012 hyhFyx x ,hFxyy xyyx 在主要边界 y = h/2 上0,0,0 yxyx 因此,在上下边界上无面力,即0,0 yx f f在次要边界 x =0, l 上x=0 (负 x 面), 2200 4123,0 hyhFf f xxyyxxx x=l (正 x 面), 220 4123, hyhFf yh12Flf xxyy3lxxx xyxy xF F Fl 此应力函数可以解决悬臂梁在 x=0 处受集中力作用的问题。 二. 半逆解法半逆解法是针对实际问题来求解的,半逆解法的具体步骤如下:逆解法没有针对具体问题进行求解, 而是找出满足相容方程的应力函数, 来考察它们能解决什么问题。这种方法可以积累弹性力学的基本解答。1. 根据弹性受力情况和边界条件等,假设部分或全部应力分量的函数形式;2. 根据 由应力推出应力函数 的形式; . yx yfx xfy xyyyxx 22222 ,3.将 代入相容方程,求出 的具体表达式; 4. 将 代入 ,求出对应的应力分量。. yx yfx xfy xyyyxx 22222 ,5. 将应力代入边界条件).()( ),()( sflm sfml ysxyy xsyxx 在s上考察它们是否满足全部边界条件(对于多连体,还须满足位移单值条件)。如果所有的条件均能满足,上述解答就是正确的解答。否则,就要修改假设,重新进行求解。 3.2 矩形梁的纯弯曲 Oy xh/2h/2yM M h1x l设有矩形截面的长梁(梁的长度 l 深度h ),它的宽度远小于深度和长度(近似的平面应力情况),或远大于深度和长度(近似的平面应变情况), 两端受相反的力偶而弯曲,体力不计。(取d=1). ay yxxyyx 0,0,6 相应的应力分量为矩形截面梁纯弯曲问题,可借助由逆解法得出的应力函数 = ay3 。显然, 满足相容方程04 Oy xh/2h/2yM M h1 x l1.考察上下两个主要边界的边界条件上下边都没有面力,要求 0,0 22 hh yyxyy . ay yxxyyx 0,0,6 此边界条件满足。2.考察左右端次要边界的边界条件左右两端没有 y 向的面力,分别要求 0,00 lxxyxxy 此边界条件也满足。x = 0, l 为小边界,可以用圣维南原理,将关于x 的边界条件用主矢量和主矩的条件代替。这些应力分量是否能满足边界条件?如能满足,a 取什么值? 2 222 d ,0d,0,0hhhh Myy ylxx lxx h1y O xh/2h/2yM Mx l . ay yxxyyx 0,0,6 将 代入上两式ayx 6 2222 d6,0d6 2hhhh Myya yya前一式总能满足,后一式要求.hMa 320,012 3 yxxyyx y,hM 代入 得. ay yxxyyx 0,0,6 注意到.hI 121 30,0 yxxyyx y,IM 得应力分量与材力结果相同。 3.3 位移分量的求出以纯弯曲矩形梁为例,说明如何由应力分量求出位移分量。(求解步骤) h1y O xh/2h/2M Ml0,0 yxxyyx y,IM 将代入 ,)1(2 1,1 xyxy xyyyxx E EE 得形变分量0, xyyx yEIM yEIM 1. 将应力分量分量代入物理方程 yuxv yv xu xyyx ,2. 将形变分量代入几何方程, 再积分求位移0, yuxv yEIMyv yEIMxu 将 代入0, xyyx yEIM yEIM 得位移分量h1y O xh/2h/2M Ml将前二式积分,得 xfyEIMv yfxyEIMu 221 2, f1, f2 为待定函数,可通过第三式求出。 将上式代入 ,得0 yuxv 0dddd 12 yyfxEIMxxf移项,得 xEIMxxfyyf dddd 21等式左边是 y 的函数,而右边是 x 的函数,因此,只可能两边都等于同一常数。于是有 xEIMxxf yyf dd,dd 21 xfyEIMv yfxyEIMu 221 2, h1y O xh/2h/2M Ml xEIMxxf yyf dd,dd 21积分,得 02201 2, vxxEIMxf uyyf 代入 xfyEIMv yfxyEIMu 221 2, 得位移分量022 022 , vxxEIMyEIMv uyxyEIMu 其中常数 , u0, v0表示刚体位移,由约束条件求得。h1y O xh/2h/2M Ml 3. 由约束条件确定常数 , u0, v0如图简支梁,约束条件是M MyO xl A .0,0,0 00000 y lxyxyx v v u 0220 22, vxxEIMyEIMv uyxyEIMu 代入.vlEIMl, , vu 0200 0200 求出 , u0, v0,就得到简支梁的位移分量 2222 yEIMxxlEIM v , ylxEIMu 有梁轴的挠度方程为 xxlEIMv y 20与材料力学的结果相同。 M MyO xl如图悬臂梁,x=l 处,对于 h/2 y h/2, 要求 u = 0, v = 0在多项式解答中这条件是无法满足的。在工程实际中这种完全固定的约束也是不大能实现的。现在,假定固定端的中点不移动,该点的水平线段也不转动。这样,约束条件是 .0,0,0 000 y lxy lxy lx xv v u 0220 22, vxxEIMyEIMv uyxyEIMu 代入.0,020 020 EIM vlEIMl, u有 .0,020 020 EIM vlEIMl, u求解得.20, 200 EIMlv, u EIM 得出悬臂梁的位移分量 22 22 yEIMxlEIM, vyxlEIMu M MyO xl梁轴的挠度方程为 20 2 xlEIMv y 与材料力学的结果相同。对于平面应变情况下的梁,须把E换为 ,把 换为 。21 E1 h1y O xh/2h/2M Ml由 可见,不论约束情况如何(不论 , u0, v0取何值)铅直线段的转角都是0uyxyEIMu xEIMyu同一横截界面上x是常数,因而是常量。 xyO PB AP AB a 于是可见,同一截面上的各铅直线段的转角相同,说明横截面保持为平面。 4. 对结果的讨论 022 22 vxxEIMyEIMv 由 可见,梁的各纵向纤维的曲率为EIMxv1 22这是材料力学中求梁的挠度时所用的基本公式。 3.4 简支梁受均布载荷设有矩形截面梁,深度为h,长度为2l,,体力可以不计,受均布载荷q,由两端的反力ql 维持平衡。(d=1 ) x y lh/2h/2Oqlq
展开阅读全文