弹性力学简明教程(第四版)课后习题解答
徐芝纶
第一章 绪论
【1-1】试举例说明什么是均匀的各向异性体,什么是非均匀的各向同性体? 【分析】均匀的各项异形体就是满足均匀性假定,但不满足各向同性假定;非均匀的各向异性体,就是不满足均匀性假定,但满足各向同性假定。
【解答】均匀的各项异形体如:竹材,木材。 非均匀的各向同性体如:混凝土。
【1-2】一般的混凝土构件和钢筋混凝土构件能否作为理想弹性体?一般的岩质地基和土质地基能否作为理想弹性体?
【分析】能否作为理想弹性体,要判定能否满足四个假定:连续性,完全弹性,均匀性,各向同性假定。
【解答】一般的混凝土构件和土质地基可以作为理想弹性体;一般的钢筋混凝土构件和岩质地基不可以作为理想弹性体。
【1-3】五个基本假定在建立弹性力学基本方程时有什么作用?
【解答】(1)连续性假定:假定物体是连续的,也就是假定整个物体的体积都被组成这个物体的介质所填满,不留下任何空隙。引用这一假定后,物体的应力、形变和位移等物理量就可以看成是连续的。因此,建立弹性力学的基本方程时就可以用坐标的连续函数来表示他们的变化规律。
完全弹性假定:假定物体是完全弹性的,即物体在对应形变的外力被去除后,能够完全恢复原型而无任何形变。这一假定,还包含形变与引起形变的应力成正比的涵义,亦即两者之间是成线性关系的,即引用这一假定后,应力与形变服从胡克定律,从而使物理方程成为线性的方程,其弹性常数不随应力或形变的大小而变。
均匀性假定:假定物体是均匀的,即整个物体是由同一材料组成的,引用这一假定后整个物体的所有各部分才具有相同的弹性,所研究物体的内部各质点的物理性质都是相同的,因而物体的弹性常数不随位置坐标而变化。
各向同性假定:假定物体是各向同性的,即物体的弹性在所有各个方向都相同,引用此假定后,物体的弹性常数不随方向而变。
小变形假定:假定位移和变形是微小的。亦即,假定物体受力以后整个物体所有各点的位移都远远小于物体原来的尺寸,而且应变和转角都远小于1。这样在建立物体变形以后的平衡方程时,就可以方便的用变形以前的尺寸来代替变形以后的尺寸。在考察物体的位移与形变的关系时,它们的二次幂或乘积相对于其本身都可以略去不计,使得弹性力学中的微分
1
方程都简化为线性的微分方程。
【1-4】应力和面力的符号规定有什么区别?试画出正坐标面和负坐标面上的正的应力和正的面力的方向。
【解答】应力的符号规定是:当作用面的外法线方向指向坐标轴方向时(即正面时),这个面上的应力(不论是正应力还是切应力)以沿坐标轴的正方向为正,沿坐标轴的负方向为负。当作用面的外法线指向坐标轴的负方向时(即负面时),该面上的应力以沿坐标轴的负方向为正,沿坐标轴的正方向为负。
面力的符号规定是:当面力的指向沿坐标轴的正方向时为正,沿坐标轴的负方向为负。 由下图可以看出,正面上应力分量与面力分量同号,负面上应力分量与面力分量符号相反。
正的应力
正的面力
【1-5】试比较弹性力学和材料力学中关于切应力的符号规定。
【解答】材料力学中规定切应力符号以使研究对象顺时针转动的切应力为正,反之为负。 弹性力学中规定,作用于正坐标面上的切应力以沿坐标轴的正方向为正,作用于负坐标面上的切应力以沿坐标轴负方向为正,反之为负。
【1-6】试举例说明正的应力对应于正的形变。
【解答】正的应力包括正的正应力与正的切应力,正的形变包括正的正应变与正的切应变,本题应从两方面解答。
正的正应力对应于正的正应变:轴向拉伸情况下,产生轴向拉应力为正的应力,引起轴向伸长变形,为正的应变。
正的切应力对应于正的切应变:在如图所示应力状态情况下,切应力均为正的切应力,引起直角减小,故为正的切应变。
2
【1-7】试画出图1-4中矩形薄板的正的体力、面力和应力的方向。 【解答】
正的体力、面力
正的体力、应力
【1-8】试画出图1-5中三角形薄板的正的面力和体力的方向。 【解答】
fxfyfyxfxfxfyfyfxyOz
【1-9】在图1-3的六面体上,y面上切应力?yz的合力与z面上切应力?zy的合力是否相等?
【解答】切应力为单位面上的力,量纲为L?1MT?2,单位为N/m2。因此,应力的合力应乘以相应的面积,设六面体微元尺寸如dx×dy×dz,则y面上切应力?yz的合力为:
?yz?dx?dz (a)
z面上切应力?zy的合力为:
?zy?dx?dy (b)
由式(a)(b)可见,两个切应力的合力并不相等。
【分析】作用在两个相互垂直面上并垂直于该两面交线的切应力的合力不相等,但对某点的合力矩相等,才导出切应力互等性。
3
第二章 平面问题的基本理论
【2-1】试分析说明,在不受任何面力作用的空间体表面附近的薄层中(图2-14)其应力状态接近于平面应力的情况。
【解答】在不受任何面力作用的空间表面附近的薄层中,可以认为在该薄层的上下表面都无面力,且在薄层内所有各点都有?z??xz??yz?0,只存在平面应力分量?x,?y,?xy,且它们不沿z方向变化,仅为x,y的函数。可以认为此问题是平面应力问题。
【2-2】试分析说明,在板面上处处受法向约束且不受切向面力作用的等厚度薄片中(2-15),当板边上只受x,y向的面力或约束,且不沿厚度变化时,其应变状态接近于平面应变的情况。
【解答】板上处处受法向约束时?z?0,且不受切向面力作用,则
?xz??yz?0(相应?zx??zy?0)板边上只受x,y向的面力或约束,所以仅存
Oz在?x,?y,?xy,且不沿厚度变化,仅为x,y的函数,故其应变状态接近于平面应变的情况。
【2-3】在图2-3的微分体中,若将对形心的力矩平很条件
y?MC试问将导出什么形?0改为对角点的力矩平衡条件,
式的方程?
【解答】将对形心的力矩平衡条件
?MC改为分别?0,
对四个角点A、B、D、E的平衡条件,为计算方便,在z方向的尺寸取为单位1。
?MA?0
??xy??xdxdydy?ydx?1??(?x?dx)dy?1??(?xy?dx)dy?1?dx??ydy?1?2?x2?x2 (a)
??y??dxdydx?(?y?dy)dx?1??(?yx?yxdy)dx?1?dy?fxdxdy?1??fxdxdy?1??0?y2?y22?MB?0 (?x???yx??y??xdydxdx)dy?1??(?yx?dy)dx?1?dy?(?y?dy)dx?1??x2?y?y2 (b)
dydxdydx??xydy?1?dx??xdy?1???ydx?1??fxdxdy?1??fydxdy?1??022224
?MD?0
(?y?dxdy??xydy?1?dx??xdy?1???yxdx?1?dy?y22 (c)
??xdxdydydx??xdx?1??(?x?dx)dy?1??fxdxdy?1??fydxdy?1??02?x222dy)dx?1???y?ME?0
dxdydx??xdy?1???yxdx?1?dy??ydx?1???y222 (d)
????dydydx(?x?xdx)dy?1??(?xy?xydx)dy?1?dx?fxdxdy?1??fydxdy?1??0?x2?x22?(?y?dy)dx?1?略去(a)、(b)、(c)、(d)中的三阶小量(亦即令dxdy,dxdy都趋于0),并将各式都除以dxdy后合并同类项,分别得到?xy??yx。
【分析】由本题可得出结论:微分体对任一点取力矩平衡得到的结果都是验证了切应力互等定理。
【2-4】在图2-3和微分体中,若考虑每一面上的应力分量不是均匀分布的,验证将导出什么形式的平衡微分方程?
【解答】微分单元体ABCD的边长dx,dy都是微量,因此可以假设在各面上所受的应力如图a所示,忽略了二阶以上的高阶微量,而看作是线性分布的,如图(b)所示。为计算方便,单元体在z方向的尺寸取为一个单位。
22??yOx??????xyAxyB??????xAyxAyA??????ABfyOxyxDyDfxD??????xDxyD??xy?A??x?A
xyC??????yxAyA??????ABfyyxDyDfxD??x?D??xy?D
??????y各点正应力:
CxB??????yByxB????????x?C?????????xyBCxByCyxCy??????yxByB????????x?C???xyCyCyxC(a) (b)
(?x)A??x;
(?y)A??y
(?x)B??x?
??xdy; ?y
(?y)B??y?5
??y?ydy
(?x)D??x???xdx; ?x
(?y)D??y???xdx ?x(?x)C??x?各点切应力:
??x??xdx??y; ?x?y(?y)C??y???y?xdx???y?y?y
(?xy)A??xy;
(?yx)A??yx
(?xy)B??xy???xy?ydy;
(?yx)A??yx???yx?ydy
(?xy)D??xy???xy?x??xy?xdx;
??xy?y
(?yx)D??yx???yx?x??yx?xdx
??yx?y(?xy)C??xy?dx?dy;
(?yx)C??yx?dx?dy
由微分单元体的平衡条件 ?Fx?0, ?Fy?0,得
??1???????x?????x????x??x????1???????dydy???dx???dx?dy??x?x???x??????dy???x?y????x???x?y?????????2??2???????yx????yyx????yx??yx?????1???1??dx???dx?????yx?dy????yx?dx?dy???dx?fxdxdy?0???yx???yx+2?x2?y?x?y??????????????????
??1??????y?????y????y??y?????1??dx???dx?????y?dy????y?dx?dy???dx?????y???y?2?x2?y?x?y????????????????????????xy?????xy????xy??xy?????1??1??dy???dy?????xy+dx????xy+dy?dx???dy?fydxdy?0???xy???xy+?y?x?y?x??????????????2??2??
以上二式分别展开并约简,再分别除以dxdy,就得到平面问题中的平衡微分方程:
??????x??yx??fx?0;y?xy?fy?0 ?x?y?y?x【分析】由本题可以得出结论:弹性力学中的平衡微分方程适用于任意的应力分布形式。 【2-5】在导出平面问题的三套基本方程时,分别应用了哪些基本假定?这些方程的适用条件是什么?
【解答】(1)在导出平面问题的平衡微分方程和几何方程时应用的基本假设是:物体的连续性和小变形假定,这两个条件同时也是这两套方程的适用条件。
(2)在导出平面问题的物理方程时应用的基本假定是:连续性,完全弹性,均匀性和各向同性假定,即理想弹性体假定。同样,理想弹性体的四个假定也是物理方程的使用条件。
6
【思考题】平面问题的三套基本方程推导过程中都用到了哪个假定?
【2-6】在工地上技术人员发现,当直径和厚度相同的情况下,在自重作用下的钢圆环(接近平面应力问题)总比钢圆筒(接近平面应变问题)的变形大。试根据相应的物理方程来解释这种现象。
【解答】体力相同情况下,两类平面问题的平衡微分方程完全相同,故所求的应力分量相同。由物理方程可以看出,两类平面问题的物理方程主要的区别在于方程中含弹性常数的系数。由于E为GPa级别的量,而泊松比?取值一般在(0,0.5),故主要控制参数为含有弹性模量的系数项,比较两类平面问题的系数项,不难看出平面应力问题的系数1/E要大于平面应变问题的系数
?1???/E。因此,平面应力问题情况下应变要大,故钢圆环变形大。
2【2-7】在常体力,全部为应力边界条件和单连体的条件下,对于不同材料的问题和两类平面问题的应力分量?x,?y和?xy均相同。试问其余的应力,应变和位移是否相同?
【解答】(1)应力分量:两类平面问题的应力分量?x,?y和?xy均相同,但平面应力问题
?z??yz??xz?0,而平面应变问题的?xz??yz?0,?z????x??y?。
(2)应变分量:已知应力分量求应变分量需要应用物理方程,而两类平面问题的物理方程不相同,故应变分量?xz??yz?0,?xy相同,而?x,?y,?z不相同。
(3)位移分量:由于位移分量要靠应变分量积分来求解,故位移分量对于两类平面问题也不同。【2-8】在图2-16中,试导出无面力作用时AB边界上的?x,?y,?xy之间的关系式
【解答】由题可得:
?gxO?y?xyA?x?nl?cos?,m?cos???90??sin?fx?AB??0,fy?AB??0将以上条件代入公式(2-15),得:
y图2-16B??x?ABcos????yx?ABsin??0, ??y?ABsin??(?xy)ABcos??0?(?x)AB????yx?ABtan????y?ABtan2?
7
【2-9】试列出图2-17,图2-18所示问题的全部边界条件。在其端部小边界上,应用圣维南原理列出三个积分的应力边界条件。
oh1h2xMqFNFSh/2h/2x?gby
q1y?h2??b?图2-17
l
图2-18
【分析】有约束的边界上可考虑采用位移边界条件,若为小边界也可写成圣维南原理的三个积分形式,大边界上应精确满足公式(2-15)。
【解答】图2-17:
上(y=0)
0 -1
0
左(x=0) -1 0
右(x=b)
1 0
l
m
fx?s? fy?g?y?h1?0
??g?y?h1?0
?s?
?gh1
代入公式(2-15)得
①在主要边界上x=0,x=b上精确满足应力边界条件:
??x?x?0???g(y?h1),??xy?x?0?0;??x?x?b???g(y?h1),??xy?x?b?0;②在小边界y?0上,能精确满足下列应力边界条件:
???yy?0???gh,??xy?y?0?0
③在小边界y?h2上,能精确满足下列位移边界条件:
?u?y?h时,可求得固定端约束反力分别为:
2?0,?v?y?h?02
这两个位移边界条件可以应用圣维南原理,改用三个积分的应力边界条件来代替,当板厚?=1Fs?0,FN???gh1b,M?0
由于y?h2为正面,故应力分量与面力分量同号,则有:
8
?b???dx???gh1b??0yy?h2??b ??0??y?y?h2xdx?0?b????dx?0???0xyy?h2⑵图2-18
①上下主要边界y=-h/2,y=h/2上,应精确满足公式(2-15)
l
0 0
m
-1 1
fx(s)
0 -q1
fy(s)
hy??
2hy?
2q
0
(?y)y?-h/2??q,(?yx)y?-h/2?0,(?y)y?h/2?0,(?yx)y?h/2??q1
②在x=0的小边界上,应用圣维南原理,列出三个积分的应力边界条件:负面上应力与面力符号相反,有
?h/2(?)dx??FS???h/2xyx?0?h/2???h/2(?x)x?0dx??FN ?h/2?(?)ydx??M???h/2xx?0③在x=l的小边界上,可应用位移边界条件ux?l?0,vx?l?0这两个位移边界条件也可改用三个积分的应力边界条件来代替。
首先,求固定端约束反力,按面力正方向假设画反力,如图所示,列平衡方程求反力:
?F?Fyx??q1l?FN??q1l?FN ?0,FN?FN?FNFS??0,FS?FS??ql?0?FS???ql?FS
M?q1lh121ql2?MA?0,M?M'?FSl?2ql?2q1lh?0?M?2?M?FSl?2
由于x=l为正面,应力分量与面力分量同号,故
?h/2(?)dy?F??ql?FN1N???h/2xx?l?q1lhql2?h/2?M?FSl? ???h/2(?x)x?lydy?M??22??h/2(?)dy?F???ql?Fxyx?lSS????h/2
9
【2-10】试应用圣维南原理,列出图2-19所示的两个问题中OA边上的三个积分的应力边界条件,并比较两者的面力是否是是静力等效?
【解答】由于hqobxAhhoFNMxqb2b/2b/2FN?Aqb2M?12l,OA为小边界,故其上可用圣维南原理,
xq by写出三个积分的应力边界条件:
(a)上端面OA面上面力fx?0,fy??a??h??b,??1?图2-19y?b?由于OA面为负面,故应力主矢、主矩与面力主矢、主矩符号相反,有
bbxqb?b?dx??fdx??qdx???0y?0b??0?y?y?02?bbx?bqb2?b???0??y?y?0xdx???0fyxdx??0q??x?dx?b?212(对OA中点取矩) ???b??0??yx?y?0dx?0?(b)应用圣维南原理,负面上的应力主矢和主矩与面力主矢和主矩符号相反,面力主矢y向为正,主矩为负,则
qb?b?dx??F??N??0?y?y?02?qb2?b ??0??y?y?0xdx??M?12??b?dx?0??0?xy?y?0?综上所述,在小边界OA上,两个问题的三个积分的应力边界条件相同,故这两个问题是静力等效的。
【2-11】检验平面问题中的位移分量是否为正确解答的条件是什么? 【解答】(1)在区域内用位移表示的平衡微分方程式(2-18); (2)在s?上用位移表示的应力边界条件式(2-19); (3)在su上的位移边界条件式(2-14); 对于平面应变问题,需将E、μ作相应的变换。
【分析】此问题同时也是按位移求解平面应力问题时,位移分量必须满足的条件。 【2-12】检验平面问题中的应力分量是否为正确解答的条件是什么? 【解答】(1)在区域A内的平衡微分方程式(2-2);
10
(2)在区域A内用应力表示的相容方程式(2-21)或(2-22);
(3)在边界上的应力边界条件式(2-15),其中假设只求解全部为应力边界条件的问题; (4)对于多连体,还需满足位移单值条件。
【分析】此问题同时也是按应力求解平面问题时,应力分量必须满足的条件。 【补题】检验平面问题中的应变分量是否为正确解答的条件是什么? 【解答】用应变表示的相容方程式(2-20)
【2-13】检验平面问题中的应力函数是否为正确解答的条件是什么? 【解答】(1)在区域A内用应力函数表示的相容方程式(2-25); (2)在边界S上的应力边界条件式(2-15),假设全部为应力边界条件; (3)若为多连体,还需满足位移单值条件。 【分析】此问题同时也是求解应力函数的条件。 【2-14】检验下列应力分量是否是图示问题的解答:
qqaaOqbxqh/2xbqOh/2y
yll?h 图2-20 图2-21
y2(a)图2-20,sx=2q,?y??xy?0。
b【解答】在单连体中检验应力分量是否是图示问题的解答,必须满足:(1)平衡微分方程(2-2);(2)用应力表示的相容方程(2-21);(3)应力边界条件(2-15)。
(1)将应力分量代入平衡微分方程式,且fx?fy?0
??x??yx??y??xy??0 ??0 显然满足 ?x?y?y?x(2)将应力分量代入用应力表示的相容方程式(2-21),有
??2?2?2q等式左=?2?2???x??y?=2?0=右
b??x?y?应力分量不满足相容方程。
因此,该组应力分量不是图示问题的解答。
MFsS*y,?xy?(取梁的厚度b=1),得出所示问题的(b)图2-21,由材料力学公式,?x?IbI
11
x3y3qx222解答:?x??2q3,?xy?-(h?4y)。又根据平衡微分方程和边界条件得出:3lh4lh3qxyxy3qx。试导出上述公式,并检验解答的正确性。 ?y??2q3?2lhlh2l【解答】(1)推导公式
在分布荷载作用下,梁发生弯曲形变,梁横截面是宽度为1,高为h的矩形,其对中性轴(Z轴)
h3的惯性矩I?,应用截面法可求出任意截面的弯矩方程和剪力方程
12q3qx2M(x)??x,F?x???。
6l2l所以截面内任意点的正应力和切应力分别为:
M?x?x3y?x?y??2q3
Ilh3Fs?x??4y2?3qx222?xy?1???.h?4y??。 ??2bh?h2?4lh3根据平衡微分方程第二式(体力不计)。
??y?y???xy?x?0
3qxyxy3.?2q3?A 得: ?y?2lhlh根据边界条件
???yy?h/2?0
qxA??.
2l得
3qxyxy3qx.?2q3?. 故 ?y?2lhlh2l将应力分量代入平衡微分方程(2-2) 第一式:
x2yx2y左??6q.3?6q3?0?右 满足
lhlh第二式 自然满足 将应力分量代入相容方程(2-23)
12
??2?2?xyxy左??2?2???x??y???12q.3?12q.3?0?右
lhlh??x?y?应力分量不满足相容方程。
故,该分量组分量不是图示问题的解答。
【2-15】试证明:在发生最大与最小切应力的面上,正应力的数值都等于两个主应力的平均值。 【解答】(1)确定最大最小切应力发生位置 任意斜面上的切应力为?n2?lm??2??1?,用关系式l2?m2?1消去m,得
2422?n??l1?l??2??1???l?l??2??1???1/4??1/2?l由上式可见当
???2??1?
112???2时,?n为最大或最小,为 ??n?max??1。因此,?l?0时,即l??min222切应力的最大,最小值发生在与x轴及y轴(即应力主向)成45°的斜面上。
(2)求最大,最小切应力作用面上,正应力?n的值 任一斜面上的正应力为
?n?l2??1??2???2
最大、最小切应力作用面上l??1/2,带入上式,得
11?n???1??2???2???1??2?
22证毕。
【2-16】设已求得一点处的应力分量,试求?1,?2,?1
(a)?x?100,?y?50,?xy?1050; (b)?x?200,?y?0,?xy??400;(c)?x??2000,?y?1000,?xy??400; (d)?x??1000,?y??1500,?xy?500.【解答】由公式(2-6)
2?1??x?1??x?1??x??y??x??y?2tan?1??1?arctan???????xy及?xy,得?xy ?2?22??2?1?100?50?100?50???????1050(a)
?2?22????2?150?? ?0?1?arctan150?100?35?16'
1050 13
2?1?200?0?5122?200?0?????400???????(b) ?2??3122?2???1?arctan512?200?arctan??0.78???37?57'
?4002?1??2000?1000?10522??2000?1000?????400???????(c) ?2?22????2052?1?arctan1052?2000?arctan??7.38???82?32'
?4002?1??1000?1500??691??1000?1500?2??????500???1809 (d)
?2?22????1?arctan?691?1000?arctan0.618?31?43'
500【2-17】设有任意形状的等候厚度薄板,体力可以不计,在全部边界上(包括孔口边界上)受有均匀压力q。试证sx=sy=-q及?xy?0能满足平衡微分方程、相容方程和应力边界条件,也能满足位移单值条件,因而就是正确的解答。
【解答】(1)将应力分量?x??y??q,?xy?0,和体力分量
Oxq?xAfyfxqy?yfx?fy?0分别带入平衡微分方程、相容方程
???x??xy??fx?0??y??x (a) ????y???xy?f?0y??x??y?2??x??y??0 (b)
显然满足(a)(b)
(2)对于微小的三角板A,dx,dy都为正值,斜边上的方向余弦l?cos?n,x?,m?cos?n,y?,将?x??y?-q,?xy?0,代入平面问题的应力边界条件的表达式(2-15),且
fx?-qcos?n,x?,fy?qcos?n,y?,则有
?xcos?n,x???qcos?n,x?,?ycos?n,y???qcos?n,y?所以?x??q,?y??q。
对于单连体,上述条件就是确定应力的全部条件。 (3)对于多连体,应校核位移单值条件是否满足。
14
该题为平面应力情况,首先,将应力分量代入物理方程(2-12),得形变分量,
?x?(???1)(??1)q,?y?q,?xy?0 (d) EE将(d)式中形变分量代入几何方程(2-8),得
?u(?-1)?v(?-1)?v?u=q,=q,??0 (e) ?xE?yE?x?y前两式积分得到
(?-1)(?-1)u=qx?f1(y),v=qy?f2(x) (f)
EE其中f1?y?,f2?x?分别任意的待定函数,可以通过几何方程的第三式求出,将式(f)代入式(e)的第三式,得
?df1(y)df2(x)?dydx
等式左边只是y的函数,而等式右边只是x的函数。因此,只可能两边都等于同一个常数?,于是有
df1(y)df(x)???,2?? dydx积分后得f1?y????y?u0,f2?x???x?v0 代入式(f)得位移分量
(??1)?u?qx??y?u0??E (g) ?(??1)?v?qy??x?v0??E其中u0,v0,?为表示刚体位移量的常数,需由约束条件求得
从式(g)可见,位移是坐标的单值连续函数,满足位移单值条件。因而,应力分量是正确的解答。
【2-18】设有矩形截面的悬臂梁,在自由端受有集中荷载F(图2-22),体力可以不计。试根据材料力学公式,写出弯应力?y?0,然后证明这些表达式满足平衡微分方程和相容方程,再说明这些表达式是否就表示正确的解答。
【解答】(1)矩形悬臂梁发生弯曲变形,任意横截面上的弯矩方程M(x)??Fx,横截面对中性轴的惯性矩为
h/2xOh/2Iz?h3/12,根据材料力学公式
15
1Fyl
弯应力?x?M(x)12Fy??3xy; Izh该截面上的剪力为Fs?x???F,剪应力为
Fs(x)S*?F6F?h2?h??h/2?y?2??xy????y?b??y???y? ???3?3?bIzh?41??h/12??2??2??取挤压应力?y?0
(2)将应力分量代入平衡微分方程检验 第一式:左??12F12Fy?y?0?右 23hh第二式:左=0+0=0=右
该应力分量满足平衡微分方程。
(3)将应力分量代入应力表示的相容方程
左??2(?x??y)?0?右 满足相容方程
(4)考察边界条件
①在主要边界y??h/2上,应精确满足应力边界条件(2-15)
l
0 0
m
-1 1
fx
0 0
fy0 0
hy??上2 hy?上2
代入公式(2-15),得
???yy?-h/2?0,??xy?y??h/2?0;??y?y?h/2?0,??yx?y?h/2?0
②在次要边界x=0上,列出三个积分的应力边界条件,代入应力分量主矢主矩
?h/2?(?x)x?0dy?0?x向面力主矢???h/2?h/2???h/2(?x)x?0ydy?0?面力主矩?2h/2?6Fh??h/22(?)dy??(?y)dy???F?y向面力主矢3????h/2xyx?0?h/2?h4???满足应力边界条件
FNFS③在次要边界上,首先求出固定边面力约束反力,按正方向假设,即面力的主矢、主矩,FN?0,FS??F,M??Fl
其次,将应力分量代入应力主矢、主矩表达式,判断是否与面力主矢与主矩等效:
M 16
12Flydy?0?FN??h/2?h/2h3
h/2h/122F(?x)x?lydy???ly2dy??Fl?M3??h/2?h/2h
h/2(?x)x?ldy???h/2
?h/2?h/2(?xy)x?ldy???h/2?h6F?h22??y?dy??F?FS3?/2h?4?
满足应力边界条件,因此,它们是该问题的正确解答。
【2-19】试证明,如果体力虽然不是常量,但却是有势的力,即体力分量可以表示为
fx???V?V,fy??,其中V是势函数,则应力分量亦可用应力函数表示成为?x?y?2??2??2??x=2?V,?y=2?V,?xy??,试导出相应的相容方程。
?y?x?x?y【解答】(1)将fx,fy带入平衡微分方程(2-2)
???x??yx???x??yx?V??f?0???0x???y?y?x??x??x (a) ? ?????y???xy?f?0???y???xy??V?0y???y?x?x?y???y将(a)式变换为
??yx???0?(?x?V)??x?y? (b) ???(??V)???xy?0y??y?y?为了满足式(b),可以取
?2??2??2??x?V?2,?y?V?2,?xy??
?y?x?x?y?2??2??2??V,?y?2?V,?xy??即?x? 2?y?x?x?y(2)对体力、应力分量fx,fy,?x,?y求偏导数,得
??fx?fy?2V?2V??2, ??2??x?x?y?y?2??2?x?4??2V?4??2V???x?4?2 (c) ?2?22?2, 2?x?y?x?y?y?y??x??2??2?y?4??2V?4??2Vy?2?4?2, ?22?22?x?x?y?x?y?y???x
17
将(c)式代入公式(2-21)得平面应力情况下应力函数表示的相容方程
??fx?fy????x??y???(1??)??? (2-21)
?x?y??2??2V?2V??4??2V?4??2V?4??2V?4??2V?2?4?2?4?2?22?2?(1??)?2?2?22?x?y?x?y?y?x?x?x?y?y?y???x
整理得:
??2V?2V?4??4??4??222?4??(1??)?2?2?x4?x?y?y?y??x即平面应力问题中的相容方程为
?? (d) ??4???(1??)?2V
将(c)式代入公式(2-22)或将(d)式中的替换为
?,的平面应变情况下的相容方程: 1???? (e) ??4??4??4?1?2???2V?2V?222?4?????x4?x?y?y1????x2?y2即 ????证毕。
41?2?2?V。 1?? 18
第三章 平面问题的直角坐标解答
【3-1】为什么在主要边界(大边界)上必须满足精确的应力边界条件式(2-15),而在小边界上可以应用圣维南原理,用三个积分的应力边界条件(即主矢量、主矩的条件)来代替?如果在主要边界上用三个积分的应力边界条件代替式(2-15),将会发生什么问题?
【解答】弹性力学问题属于数学物理方程中的边值问题,而要使边界条件完全得到满足,往往比较困难。这时,圣维南原理可为简化局部边界上的应力边界条件提供很大的方便。将物体一小部分边界上的面力换成分布不同,但静力等效的面力(主矢、主矩均相同),只影响近处的应力分布,对远处的应力影响可以忽略不计。如果在占边界绝大部分的主要边界上用三个积分的应力边界条件来代替精确的应力边界条件(公式2-15),就会影响大部分区域的应力分布,会使问题的解答精度不足。
【3-2】如果在某一应力边界问题中,除了一个小边界条件,平衡微分方程和其它的应力边界条件都已满足,试证:在最后的这个小边界上,三个积分的应力边界条件必然是自然满足的,固而可以不必校核。
【解答】区域内的每一微小单元均满足平衡条件,应力边界条件实质上是边界上微分体的平衡条件,即外力(面力)与内力(应力)的平衡条件。研究对象整体的外力是满足平衡条件的,其它应力边界条件也都满足,那么在最后的这个次要边界上,三个积分的应力边界条件是自然满足的,因而可以不必校核。
【3-3】如果某一应力边界问题中有m个主要边界和n个小边界,试问在主要边界和小边界上各应满足什么类型的应力边界条件,各有几个条件?
【解答】在m个主要边界上,每个边界应有2个精确的应力边界条件,公式(2-15),共2m个;在n个次要边界上,如果能满足精确应力边界条件,则有2n个;如果不能满足公式(2-15)的精确应力边界条件,则可以用三个静力等效的积分边界条件来代替2个精确应力边界条件,共3n个。
【3-4】试考察应力函数??ay在图3-8所示的矩形板和坐标系中能解决什么问题(体力不计)?
【解答】⑴相容条件:
不论系数a取何值,应力函数??ay总能满足应力函数表示的相容方程,式(2-25).
⑵求应力分量
当体力不计时,将应力函数?代入公式(2-24),得
33Oxhyl图3-8?x?6ay,?y?0,?xy??yx?0
19
⑶考察边界条件
上下边界上应力分量均为零,故上下边界上无面力. 左右边界上;
当a>0时,考察?x分布情况,注意到?xy?0,故y向无面力 左端:fx?(?x)x?0?6ay ?0?y?h? fy????xyx?0?0
右端:fx???x?x?l?6ay (0?y?h) fy?(?xy?xl )?0应力分布如图所示,当l?h时应用圣维南原理可以将分布的面力,等效为主矢,主矩
OA xfxy偏心距e:
因为在A点的应力为零。设板宽为b,集中荷载p的偏心距e:
fx
主矢的中心在矩下边界位置。即本题情况下,可解决各种偏心拉伸问题。
eePP(?x)A?ppe?2?0?e?h/6 bhbh/6同理可知,当a<0时,可以解决偏心压缩问题。 【3-5】取满足相容方程的应力函数为:⑴
??ax2y,⑵??bxy2,⑶??cxy3,试求出应力O分量(不计体力),画出图3-9所示弹性体边界上的面力分布,并在小边界上表示出面力的主矢量和主矩。
【解答】(1)由应力函数??axy,得应力分量表达式
2h/2h/2xyl图3-9(l?h)?x?0,?y?2ay,?xy??yx??2ax
??(l?x?m?yx)s?fx(s)考察边界条件,由公式(2-15)?
??(m?y?l?xy)s?fy(s)h①主要边界,上边界y??上,面力为
2hhfx(y??)?2ax fy(y??)?ah
22h②主要边界,下边界y?,面力为
2
20
hhfx(y?)??2ax, fy(y?)?ah
22③次要边界,左边界x=0上,面力的主矢,主矩为 x向主矢:Fx??y向主矢:Fy??主矩:M???h/2?h/2h/2(?x)x?0dy?0 (?xy)x?0dy?0
??h/2?h/2?h/2(?x)x?0ydy?0
?yxah?xyOx次要边界,右边界x=l上,面力的主矢,主矩为 x向主矢:Fx???y向主矢:Fy??主矩:M?h/2?h/2h/2(?x)x?ldy?0 (?xy)x?ldy??h/2?h/2??h/2(?2al)dy??2alh
yah2al?h/2?h/2(?x)x?lydy?0
2al弹性体边界上面力分布及次要边界面上面力的主矢,主矩如图所示 ⑵??bxy
将应力函数代入公式(2-24),得应力分量表达式
2?x?2bx,?y?0,?xy??yx??2by
考察应力边界条件,主要边界,由公式(2-15)得 在y??h?h?h??主要边界,上边界上,面力为fx?y????bh,fy?y????0
2?2?2??在y?h?h?h??,下边界上,面力为fx?y????bh,fy?y???0
2?2?2??在次要边界上,分布面力可按(2-15)计算,面里的主矢、主矩可通过三个积分边界条件求得:
在左边界x=0,面力分布为fx?x?0??0,fy?x?0??2by 面力的主矢、主矩为 x向主矢:Fx??????h2h?2h2h?2xx?0dy?0
h2h?2y向主矢:Fy??主矩;M??????xyx?0dy?????2by?x?0dy?0
?h/2?h/2(?x)x?0ydy?0
在右边界x=l上,面力分布为
21
fx?x?l??2bl,fy?x?l???2by
面力的主矢、主矩为 x向主矢:Fx??y向主矢:Fy'?主矩:M'?h/2?????h/2h/2xx?ldy??h/2?h/22bldy?2blh
?h/2?h/2??xy?x?ldy??ydy??h/2?h/2??2by?dy?0
?????h/2h/2xx?l?h/22blydy?0
弹性体边界上的面力分布及在次要上面力的主矢和主矩如图所示
ahO2alahx?xy(3)??cxy
3y?xy 将应力函数代入公式(2-24),得应力分量表达式
?x?6cxy,?y?0,?xy??yx??3cy2
考察应力边界条件,在主要边界上应精确满足式(2-15) ①上边界y??上,面力为
h2h?3h???fx?y????ch2,fy?y????0
2?42???② 下边界y=h上,面力为 2h?3h???fx?y????ch2,fy?y???0
2?42???次要边界上,分布面力可按(2-15)计算,面力的主矢、主矩可通过三个积分边界求得: ③左边界x=0上,面力分布为
fx?x?0??0,fy?x?0??3cy2面力的主矢、主矩为x向主矢:Fx???y向主矢:Fy???主矩:M???
h/2-h/2h/2?h/2h/2??x?x?0dy?0
h/2?h/2???xyx?0dy????h/2ch??3cy?dy?1423??x?x?0ydy?022
④右边界x?l上,面力分布为
fx?x?l??6cly,fy?x?l???3cy2
面力的主矢、主矩为 x向主矢Fx???????h/2h/2xx?ldy??h/2?h/26clydy?0
y向主矢:Fy??主矩:M??h/2?h/2?h/2??y?x?ldy???????h/2xx?l132?3cydy??ch ???h/24h/21ydy??6cly2dy?clh3
?h/22h/2弹性体边界上的面力分布及在次要边界上面力的主矢和主矩,如图所示
【3-6】试考察应力函数??Fxy(3h2?4y2),32h能满足相容方程,并求出应力分量(不计体力),画出图3-9所示矩形体边界上的面力分布(在小边界上画出面力的主矢量和主矩),指出该应力函数能解决的问题。
【解答】(1)将应力函数代入相容方程(2-25)
Oh/2h/2xyl图3-9(l?h)?4??4??4??222?4?0,显然满足 4?x?x?y?y(2)将?代入式(2-24),得应力分量表达式
3F4y212Fxy(1?2) ?x??,?y?0,?xy??yx??32hhh(3)由边界形状及应力分量反推边界上的面力: ①在主要边界上(上下边界)上,y??h,应精确满足应力边界条件式(2-15),应力2???yy??h/2?0,??yx?y??h/2?0
因此,在主要边界y??h?h?h??上,无任何面力,即fx?y????0,fy?y????0
2?2?2??②在x=0,x=l的次要边界上,面力分别为:
3F?4y2?x?0:fx?0,fy??1-2?
2h?h?
23
x?l:fx??12Flyh33F?4y2?,fy???1?2?
2h?h?因此,各边界上的面力分布如图所示:
③在x=0,x=l的次要边界上,面力可写成主矢、主矩形式:
x=0上 x=l上
x向主矢:FN1=?y向主矢:FS1=?主矩:M1=?h/2-h/2h/2?h/2h/2fxdy?0, FN2??fydy?F, FS2??h/2?h/2h/2fxdy?0fydy??Ffxydy??Fl
?h/2?h/2h/2fxydy?0, M2???h/2因此,可以画出主要边界上的面力,和次要边界上面力的主矢与主矩,如图:
(a) (b)
因此,该应力函数可解决悬臂梁在自由端受集中力F作用的问题。
24
qx2y3yqy2y3y(?43?3?1)?(23?)能满足相容方程,并考察它在【3-7】试证??4hh10hh图3-9所示矩形板和坐标系中能解决什么问题(设矩形板的长度为l,深度为h,体力不计)。
【解答】(1)将应力函数?代入式(2-25)
h/2h/2xOyl图3-9(l?h)?4??4?24qy?0,4?3,4?yh?x?4??12qy?24qy222?2?? ?x?yh3h3代入(2-25),可知应力函数?满足相容方程。 (2)将?代入公式(2-24),求应力分量表达式:
?2?6qx2y4qy33qy?x?2?fxx??3?3?
?yhh5h?2?q4y33y?y?2?fyy?(?3??1)
?x2hh?xy??yx?2?6qxh2????3(?y2)
?x?yh4(3)考察边界条件,由应力分量及边界形状反推面力: ①在主要边界y??h(上面),应精确满足应力边界条件(2-15) 2h?h???fx?y???????yx??0,fy?y???????y??qy??h/2y??h/222????h在主要边界y??下面?,也应该满足?2?15?2 fx?y?h/2????yx??0,fy?y?h/2????y??0y?h/2y?h/2在次要边界x?0上,分布面力为fx?x?0?????x?x?03qy4qy3??3,fy?x?0?????xy??0x?05hh应用圣维南原理,可写成三个积分的应力边界条件:
25
FN??FS??M??h/2?h/2h/2?3qy4qy3?fxdy????3?dy?0?h/25hh??h/2?h/2h/2fydy?0?3qy4qy3?fxydy????3?ydy?0?h/25hh??h/2?
?h/2④在次要边界x?l上,分布面力为
fx?x?l????x?x?l6ql2y4qy33qy??3?3?
hh5hfy?x?l????xy?x?l?6ql?h2??3??y2?
h?4?应用圣维南原理,可写成三个积分的应力边界条件:
23?6qly4qy3qy?FN???fx(x?l)dy????3?3??dy?0?h/2?h/2hh5h??2h/2h/2?6ql?h??Fs???fy(x?l)dy????3??y2??dy??ql?h/2?h/2???h?4h/2h/2
M'??h/2?h/2?6ql2y4qy33qy?12fx(x?l)ydy????3?3?ydy??ql??h/2hh5h2??h/2综上,可画出主要边界上的面力分布和次要边界上面力的主矢与主矩,如图
qqqloyx12ql2
(a) (b)
因此,此应力函数能解决悬臂梁在上边界受向下均布荷载q的问题。 【3-8】设有矩形截面的长竖柱,密度为ρ,在一边侧面上受均布剪力q(图3-10),试求应力分量。
【解答】采用半逆法求解。
由材料力学解答假设应力分量的函数形式。 (1)假定应力分量的函数形式。
根据材料力学,弯曲应力?y主要与截面的弯矩有关,剪应力?xy主要与
hobxq?gy(h?b)图3-10截面的剪力有关,而挤压应力?x主要与横向荷载有关,本题横向荷载为零,则?x?0
26
(2)推求应力函数的形式
将?x?0,体力fx?0,fy??g,代入公式(2-24)有
?2??x?2?fxx?0
?y对y积分,得
???f?x? (a) ?y??yf?x??f1?x? (b)
其中f?x?,f1?x?都是x的待定函数。 (3)由相容方程求解应力函数。 将(b)式代入相容方程(2-25),得
d4f?x?d4f1?x?y??0 (c) 44dxdx在区域内应力函数必须满足相容方程,(c)式为y的一次方程,相容方程要求它有无数多个根(全竖柱内的y值都应满足它),可见其系数与自由项都必须为零,即
d4f?x?d4f1?x??0,?0 dx4dx两个方程要求
f?x??Ax3?Bx2?Cx,f1?x??Dx3?Ex2 (d)
f?x?中的常数项,f1?x?中的常数项和一次项已被略去,因为这三项在?的表达式中
成为y的一次项及常数项,不影响应力分量。将(d)式代入(b)式,得应力函数
??y?Ax3?Bx2?Cx???Dx3?Ex2? (e)
(4)由应力函数求应力分量
?2??x?2?fxx?0 (f)
?y?2??y?2?fyy?6Axy?2By?6Dx?2E??gy (g)
?x?xy?2?????3Ax2?2Bx?C (h)
?x?y 27
(5)考察边界条件
利用边界条件确定待定系数A、B、C、D、E。 主要边界x?0上(左):
??x?x?0?0,(?xy)x?0?0
将(f),(h)代入
??x?x?0?0,自然满足
(?xy)x?0??C?0 主要边界x?b上,
??x?x?b?0,自然满足
(?xy)x?b?q,将(h)式代入,得
(?xy)x?b??3Ab2?2Bb?C?q 在次要边界y?0上,应用圣维南原理,写出三个积分的应力边界条件:
?bb0(?y)y?0dx??0?6Dx?2E?dx?3Db2?2Eb?0 ?b0(?b2y)y?0xdx??0?6Dx?2E?xdx?2Db3?Eb?0 ?b0(?byx)y?0dx??0??3Ax2?2Bx?C?dx??Ab3?Bb2?Cb?0由式(i),(j),(k),(l),(m)联立求得
A??qb2, B?qb, C?D?E?0 代入公式(g),(h)得应力分量
?qx?b??1?3x?x?0, ?y?2b????gy, ?q?3?xy?bx??bx?2?? 28
(i)(j)k)(l)m)
(
(
【3-9】图3-11所示的墙,高度为h,宽度为b,h?b,在两侧面上受到均布剪力q的作用,试应用应力函数
oqb/2b/2x??Axy?Bxy求解应力分量。
【解答】按半逆解法求解。
⑴将应力函数代入相容方程(2-25)显然满足。 3qhy(h?b)⑵由公式(2-24)求应力分量表达式,体力为零,有
图3-11??2?y0,??2??2?x?2?y??x2?6Bxy,?xy??yx???x?y??A?3Bx2?
⑶考察边界条件:
在主要边界x??b2上,精确满足公式(2-15)
??x?x??b/2?0,(?xy)x??b/2??q
第一式自然满足,第二式为
?A?34Bb2??q ②在主要边界x=b/2上,精确满足式(2-15)
??x?x?b/2?0,??xy?x?b/2??q
第一式自然满足,第二式为
?A?34Bb2??q ③在次要边界y=0上,可用圣维南原理,写出三个积分的应力边界条件:
?b/2?b/2??y?y?0dx?0 满足
?b/2?b/2??y?y?0xdx?0 满足
3?b/2?dx?b/2?b/2?yx?y?0??b/2??A?3Bx2?dx??Ab?14Bb?0 联立(a)(c)得系数
A??q2,B?2qb2
代入应力分量表达式,得
??12qq?x2?x?0,?yb2xy,?xy?2??1?12b2??
29
(a)
(b)
(c)
【3-10】设单位厚度的悬臂梁在左端受到集中力和力矩作用,体力可以不计,l?h(图3-12),试用应力函数??Axy?By?Cy?Dxy求解应力分量。
【解答】采用半逆解法求解
(1)将应力函数代入相容方程(2-25),显然满足 (2)由应力函数求应力分量,代入公式(2-24)
233???2B?6By?6Dxy?x??????0?y? (a) ?2??????A?3Dy????yx?xy?(3)考察边界条件
①主要边界y??h/2上,应精确满足应力边界条件
???yy??h/2?0, 满足
??xy?y??h/2?0, 得A?Dh2?0 (b)
②在次要边界x=0上,应用圣维南原理,写出三个积分的应力边界条件
34FN 2hh/2h/22M ?ydy??M?2B?6Cyydy??M?C??????3??h/2xx?0??h/2hh/2h/2132???dy??F??A?3Dydy??F?Ah?Dh?Fs (c) ??ss??h/2?xy?x?0??h/2??4??h/2??x?x?0dy??FN???h/2?2B?6Cy?dy??FN?B??h/2h/2联立方程(b)(c)得
A?3Fs2F,D??3s 2hh最后一个次要边界?x?l?上,在平衡微分方程和上述边界条件均已满足的条件下是必然满足的,故不必在校核。
将系数A、B、C、D代入公式(a),得应力分量
???x??FN?12My?12Fsxy?hh3h3? ??y?0?2????3FS?1?4y??2??xy2hh???【3-11】设图3-13中的三角形悬臂梁只受重力作用,而梁的
30
密度为?,试用纯三次式的应力函数求解。
【解答】采用半逆解法求解
(1) 检验应力函数是否满足相容方程(2-25)
设应力函数?=Ax?Bxy?Cxy?Dy,不论上式中的系数如何取值,纯三次式的应力函数总能满足相容方程(2-25)
(2) 由式(2-24)求应力分量
由体力分量fx?0,fy??g,将应力函数代入公式(2-24)得应力分量:
3223?2??x?2?fxx?2Cx?6Dy (a)
?y?2??y?2?fyy?6Ax?2By??gy (b)
?y?xy?2?????2Bx?2Cy (c)
?x?y(3)考察边界条件:由应力边界条件确定待定系数。 ①对于主要边界y?0,其应力边界条件为:
(?y)y?0?0将式(d)代入式(b),(c),可得
,
(?yx)y?0?0 (d)
A?0,B=0 (e)
②对于主要边界y?xtan?(斜面上),应力边界条件:
在斜面上没有面力作用,即fx?fy?0,该斜面外法线方向余弦为,l??sin?,
m?cos?.由公式(2-15),得应力边界条件
?sin??(?x)y?xtan??cos??(?yx)y?xtan??0?? (f)
?sin??(?xy)y?xtan??cos??(?y)y?xtan??0?将式(a)、(b)、(c)、(e)代入式(f),可解得
C??g2cot?,D???g3cot2? (g)
将式(e)、(g)代入公式(a)、(b)、(c),得应力分量表达式:
??x??gxcot??2?gycot2?? ??y???gy???xy???gycot?【分析】本题题目已经给定应力函数的函数形式,事实上,也可通过量纲分析法确定应
31
力函数的形式。
按量纲分析法确定应力函数的形式:三角形悬臂梁内任何一点的应力与?,x,y和?g有关。由于应力分量的量纲是L?1MT?2,而x,y的量纲是L,?g的量纲是L?1MT?2,又是量纲—的数量,因此,应力分量的表达式只可能是x和y的纯一项式,即应力分量的表达式只可能是A?gx,B?gy这两种项的结合,其中A,B是量纲一的量,只与?有关。应力函数又比应力分量的长度量纲高二次,即为x和y的纯三次式,故可假设应力函数的形式为
??Ax3?Bx2y?Cxy2?Dy3。
【3-12】设图3-5中简支梁只受重力作用,而梁的密度为?,试用§3-4中的应力函数(e)求解应力分量,并画出截面上的应力分布图。
【分析】与§3-4节例题相比,本题多了体力分量
fx?0,fy??g。去除了上边界的面力。依据§3-4,应
力分量的函数形式是由材料力学解答假设的。
【解答】按半逆解法求解。
(1)由§3-4可知应力函数的函数形式为
ABx2??(Ay3?By2?Cy?D)?x(Ey3?Fy2?Gy)?y5?y4?Hy3?Ky2,由§3-4可
2106
知,?必然满足相容方程(2-25)。
(2)应力分量的表达式:
x2?x?(6Ay?2B)?x(6Ey?2F)?2Ay3?2By2?6Hy?2K (a)
2?y?Ay3?By2?Cy?D??gy (b)
?xy??x(3Ay2?2By?C)?(3Ey2?2Fy?G) (c)
【注】?y项多了-?gy
这些应力分量是满足平衡微分方程和相容方程的。因此,如果能够适当选择常数、(b)、(c)就是正确的解A、B、???、K,使所有的边界条件都被满足,则应力分量式(a)答。
(3)考虑对称性
因为yz面是梁和荷载的对称面,所以应力分布应当对称于yz面。这样?x和?y是x的
32
偶函数,而?xy是x的奇函数,于是由式(a)和式(c)可见
E?F?G?0 (d)
(4)考察边界条件:
①在主要边界y??h2上,应精确满足应力边界条件(2-15),
(?y)y??h2?0,(?yx)y??h2?0
将应力分量式(b)、(c)代入,并注意到E?F?G?0,可得:
?h3h2h?gA?B?C?D?h?0?8422?3h2h?g?h?A?B?C?D?h?0??8422 ?32??x(Ah?hB?C)?0?4?3??x(Ah2?hB?C)?0??4联立此四个方程,得:
A??2?g3,B?0,C??g,D?0 (e) h22将式(d)、(e)代入式(a)、(b)、(c)
?x??6?g24?g3xy?y?6Hy?2K (f) 22hh2?g?g?y??2y3?y (g)
h26?g3?g?xy?2xy2?x (h)
h2②考察次要边界条件
由于问题的对称性,只需考虑其中的一边,如右边。右边界x?l上,fx?0,不论y取任何值(?h2?y?h2),都有?x?0。由(f)式可见,这是不可能的,除非?,H,K均为零。因此,只能用应力?x的主矢、主矩为零,即
??将(f)式代入式(i)得
h/2h/2?h/2h/2(?x)x?ldy?0 (i) (?x)x?lydy?0 (j)
?h/24?g3?6?g2??xy?y?6Hy?2K?dy?0 22??h/2?h?h?积分后得 K=0 (k)
33
将式(f)代入式(i),得
4?g3?6?g2??ly?y?6Hy?2K??ydy?0 2??h/2?h2h?h/2积分后得
l21H??g(2?) (l)
h10将(k)、(l)代入式(f),得
6?g24?g3l21?x??2xy?2y?6?g(2?)y (m)
hhh10考察右边界上切应力分量?xy的边界条件: 右边界上fy???glh,则?xy的主矢为
?h/2?h/2??xy?x?ldy???6?g23?gxy??h/2?h22?h/2?x?dy???glh?fy ?x?l可知满足应力边界条件。
将式(g),(h),(m)略加整理,得应力分量的最后解答:
?6?g24?g3l21??X??h2xy?h2y?6?g(h2?10)y?2?g3?g? (n) y??y??2y?h2?6?g23?g???xy?xxy2?h2?(5)应力分量及应力分布图
h3h2y2?梁截面的宽度取为1个单位,则惯性矩I?,静矩是S?。
1282根据材料力学截面法可求得截面的内力,可知梁横截面上的弯矩方程和剪力方程分别为
l2?x2M?x???gh,Fs?x????ghx
2则式(n)可写成:
34
M?x??4y23y??gy(2?)??x?Ih5???gy2 y(1?42)??y?2h??Fs?x?S??xy?bI?【分析】比较弹性力学解答与材料力学解答,可知,只有切应力?xy完全相同,正应力
?x中的第一项与材料力学结果相同,第二项为弹性力学提出的修正项;?y表示纵向纤维间
的挤压应力,而材料力学假设为零。对于l>>h的浅梁,修正项很小,可忽略不计。
【3-13】图3-14所示的悬臂梁,长度为l,高度为h,l?h,在上边界受均布荷载q,试检验应力函数??Ay?Bxy?Cy?Dx?Exy能否成为此问题的解?如可以,试求出应力分量。
【解答】用半逆解法求解。 (1)相容条件:
将应力函数?代入相容方程式(2-25),得
523322120Ay?24By?0
要使?满足相容方程,应使
1A??B (a)
5(2)求应力分量,代入式(2-24)
??x?20Ay3?6Bx2y?6Cy?20Ay3?30Ax2y?6Cy??33 (b) ??y?2By?2D?2Ey??10Ay?2D?2Ey?22???6Bxy?2Ex?30Axy?2Ex?xy?(3)考察边界条件
①在主要边界y??h2上,应精确到满足应力边界条件
103Ah?2D?Eh?0 (c) 810(?y)y??h2??q,即Ah3?2D?Eh??q (d)
830(?yx)y??h2?0,即Axh2?2Ex?0 (e)
4(?y)y?h2?0,即-联立式(a)、(c)、(d)、(e),可得:
35
A?qq3qq (f) ,D??,E?,B??5h344hh3 ②在次要边界x?0上,主矢和主矩都为零,应用圣维南原理,写出三个积分的应力边界条件:
h/2?h/2?h/2(?x)x?0dy?0 满足条件
h/23Ah53(?)ydy?(20Ay?6Cy)ydy?0??Ch?0 (g) ??h/2xx?0??h/22?h/2?h/2(?xy)x?0dy?0 满足
将A的值带入(g),得
C=?
q
10h
将各系数代入应力分量表达式(b),得
??y(4y23x2?x?q?hh2?5?6h2)??qyy3?y??(1?3?4?2hh3)
????3qxy2xy??2h(1?4h2)【3-14】矩形截面的柱体受到顶部的集中力2F和力矩M的作用(图
3-15),不计体力,试用应力函数
??Ay2?Bxy?Cxy3?Dy3求解其应力分量。
【解答】采用半逆解法求解。 (1) 相容条件:
将应力函数代入相容方程(2-25),显然满足。 (2) 求应力分量:将?代入(2-24)
????x?2A?6Cxy?6Dy????y?0??? ???xy??B?3Cy2???(3) 考察边界条件。
①在主要边界y??b/2上,应精确满足应力边界条件 ??y?y??b/2?0 满足
36
(h) (a)
???xyy??b/23??q,?B?Cb2?q (b)
42b/2?b/2②在次要边界x=0上,可用圣维南原理,写出三个积分应力边界条件
?b/2?b/2(?x)x?0dy??F (2Ay?3Dy)??F (c)
?123?(?)ydy??MAy?2Dy??M (d) xx?0????b/2?2??b/2b/2b/2?????b/2xyb/2x?ody??F ??By?Cy3?b/2?b/2??F (e)
联立(b)、(c)、(d)、(e)式得
A??1?3F?2?F?F2MC?q?,B???q?,, (f) D????2?32bbb2bb????将各系数据(f)代入式(a),得应力分量解答
?F12?F?12M??x???2?q??xy?3ybb?b?b?? ??y?0????1?q?3F??6?q?F?y2xy?????2?b?b2?b??【分析】本题题目中原教材给出的坐标轴有误,无法计算。x,y坐标互换后可以计算,但计算结果与题目提示解答几乎完全不同,又将y轴调为水平向左为正方向,才得到提示结果。可见,在求解问题时,坐标轴的方向及原点的位置与解答关系密切,坐标轴不同可得到完全不同的结果。
【3-15】挡水墙的密度为?1,厚度为b(图3-16),水的密度为?2,试求应力分量。
【解答】(1)假设应力分量的函数形式。因为在y??b/2边界上,
?y?0;y?b/2边界上,?y???2gx,所以可以假设在区域内?y为
?y?xf?y?
(2)推求应力函数的形式。由?y推求?的形式
?2??y?2?xf?y??x ??x2?f?y??f1?y? ?x2
37
x3??f?y??xf1?y??f2?y?
6(3)由相容方程求应力函数。将?代入???0,得
4d4f1d4f2x3d4fd2f?x4??2x2?0
6dy4dydy4dy要使上式在任意的x处都成立,必须
d4f32?0?f(y)?Ay?By?Cy?D;4dyd4f1d2fA5B4?2?0?f(y)??y?y?Gy3?Hy2?Iy; 142dydy106d4f2?0?f2(y)?Ey3?Fy24dy代入?即得应力函数的解答,其中已经略去了与应力无关的一次项,得应力函数为:
x3Ay5By4323232??(Ay?By?Cy?D)?x(???Gy?Hy?Iy)?(Ey?Fy)
6106(4)由应力函数求应力分量,将?代入公式(2-24),注意体力fx??1g,fy?0,求得应力分量表达式
?2?B???x?2?fxx?x3?Ay???x??2Ay3?2By2?6Cy?2H???y3?? ?6Ey?2F???1gx?2??y?2?fyy?x?Ay3?By2?Cy?D??x?2?x22B3?A??xy?????3Ay2?2By?C???y4?y?3Gy2?2Hy?I??x?y23?2?
(5)考察边界条件
在主要边界y??b/2上,应精确满足应力边界条件
???yy?b/2???yy??b/2?b3?b2b???2gx?x?A?B?C?D????2gx42?8???b3b2b?0 ?x??A?B?C?D??0842????b4x2?3b2b33b2?0 ???A?Bb?C???A?B?G2?4124??32?Hb?I??0?
???xyy??b/2 38
由上式得到
3b2A?Bb?C?0 4b4b33b2A?B?G32124求解各系数,得
Hb?I?0
A?231?g,B?0,C???g,D???2g,H?0 22b32b2b3b2I??2g?G (a)
164在次要边界x?0上,列出三个积分的应力边界条件
?????b/2xb/2x?0dy?0 ? F?0
?????b/2xb/2b/2x?0ydy?0 ? E?0bb2??b/2??xy?x?0dy?0 ? I?80?2g?4G (b)
由式(a)、(b)解出
I??b1?2g,G??2g 8010b将各系数代入应力分量的表达式,得
?2?2g33?g4?2g3??x?3xy?2xy?3xy??1gxb5bb???2y33y1???? ??y??2gx?3?2b2??b??2??y33y3b?2?3y??xy???2gx?2????2gy??3???b4bb10b80y??????
【3-16】试分析简支梁受均布荷载时,平截面假设是否成立?
【解答】弹性力学解答和材料力学解答的差别是由于各自的解法不同。简言之,弹性力学的解法是严格考虑区域内的平衡微分方程、几何方程和物理方程。以及在边界上的边界条件而求解的,因而得出的解答较精确。而在材料力学的解法中,没有严格考虑上述条件,因而得出的是近似的解答。例如,材料力学引用了平面截面假设而简化了几何关系,但这个假设对于一般的梁是近似的。所以,严格地说,平截面假设不成立。
39
【3-17】试证明刚体位移u0,v0和?0实际上表示弹性体中原点的平移和转动分量,并应用§3-3的解答加以验证(注:微分体的转动分量??1??v?u????) 2??x?y?【解答】为了区分原点的转动分量与任意点处的转动分量,定义原点的转动分量为?0,任意点处的转动分量为?。
由§3-3可知,任意点处的平动分量为:
M?u?xy??0y?u0??EI ??v???My2?Mx2??x?v00??2EI2EI则任意点处的转动分量为
MM1??????1?M?x??0???x??0 ????????x??0?EIEI2??x?y?2?EI?因此,原点的平动和转动分量,即x=y=0时
u?u0,v?v0,???0
得证。
40