第四章:弹性力学问题的解法
合集下载
弹性力学第四章平面问题的极坐标解答

圆环或圆筒受均布压力(1)
q2 q1
边界条件:
圆环或圆筒受均布压力(2)
q2
q1
两个方程三个未知数,不能求解A,B,
C。因此,需引入位移单值条件:
该项必须为零,否则在环上同一点有两 个不同的位移,故B=0
圆环或圆筒受均布压力(3)
பைடு நூலகம்q2
q1 因此,得到圆筒受均匀压力的拉梅 ( me,1795—1870 ,法国)解答:
小孔口问题的特点:
1.集中性,孔附近的应力远大于较远处的应力。
2.局部性,孔口附近的应力扰动主要发生在距孔 边1.5倍孔口尺寸的范围内。在此区域外,由于开 孔引起的应力扰动一般小于5%,可以忽略不计。
注:圆孔的应力集中程度较低,有凹尖角的孔口 应力集中程度较高,因此,在设计结构时应尽量 避免有凹尖角的孔口。
o
x 在仅有径向位移的情况下,段
P P’ A
PA没有转动,因此:
A’
B
C
y
B’
极坐标中的几何方程(5)
— 纯环向位移下的线应变
o
x
很小,导致P’’A’’与PA
P P’’
的差别可以忽略,因此:
A
B B’’
D
D’
A’’
y
极坐标中的几何方程(6)
— 纯环向位移下的切应变
o
x
P
P’’
A
B B’’
D
D’
A’’
阶,因此假定:
半面体在边界上受集中力(2)
F
ao
c
ρ
代入极坐标中的相容方程:
b
得到:
半面体在边界上受集中力(3)
代入:
F
ao
弹性力学第四章:平面问题的极坐标解答2

π 2
r
σr +P θ 3σr −σθ 2σcos θ
x3 σx =− π (x2 + y2)2 2P xy2 σy =− π (x2 + y2)2 2P x2 y τxy =− π (x2 + y2)2 2P
2. 位移分量
假定为平面应力情形。 假定为平面应力情形。其极坐标形式的物理方程为
P
O y
由楔形体受集中力的情形, 由楔形体受集中力的情形,可以得到 P
O y
(令 β =0 ,α =π) : 2P cosθ σr = − ( ) π r (4-26) ) σθ =0 —— 极坐标表示的应力分量 极坐标表示的应力分量 τrθ =τθr =0
利用极坐标与直角坐标的应力转换式( ), ),可求得 利用极坐标与直角坐标的应力转换式(4-7),可求得
∂r ϕ = f (r)sinθ
θ
ϕ = f (r) (M =常 ) 数
ϕ = f (r)sinθ ϕ = f (r)cosθ (M = P⋅rsinθ) (M =M+P⋅rcosθ)
附1:曲梁应力函数确定的基本方法 :
思路: 思路: 与直梁确定应力函数的方法类似, 与直梁确定应力函数的方法类似,借且于 梁截面上应力与内力 弯矩、剪力) 应力与内力( 梁截面上应力与内力(弯矩、剪力)的关 应力与应力函数间微分关系, 系、应力与应力函数间微分关系,来推断 应力函数的分离变量形式。 应力函数的分离变量形式。 梁截面上的应力内力的关系: 梁截面上的应力内力的关系:
θ
M = Py = P⋅rsinθ
由材料力学初等理论,可知截面上正应力 由材料力学初等理论, 由此假定: 由此假定:
σθ ∝M(= P⋅rsinθ)
r
σr +P θ 3σr −σθ 2σcos θ
x3 σx =− π (x2 + y2)2 2P xy2 σy =− π (x2 + y2)2 2P x2 y τxy =− π (x2 + y2)2 2P
2. 位移分量
假定为平面应力情形。 假定为平面应力情形。其极坐标形式的物理方程为
P
O y
由楔形体受集中力的情形, 由楔形体受集中力的情形,可以得到 P
O y
(令 β =0 ,α =π) : 2P cosθ σr = − ( ) π r (4-26) ) σθ =0 —— 极坐标表示的应力分量 极坐标表示的应力分量 τrθ =τθr =0
利用极坐标与直角坐标的应力转换式( ), ),可求得 利用极坐标与直角坐标的应力转换式(4-7),可求得
∂r ϕ = f (r)sinθ
θ
ϕ = f (r) (M =常 ) 数
ϕ = f (r)sinθ ϕ = f (r)cosθ (M = P⋅rsinθ) (M =M+P⋅rcosθ)
附1:曲梁应力函数确定的基本方法 :
思路: 思路: 与直梁确定应力函数的方法类似, 与直梁确定应力函数的方法类似,借且于 梁截面上应力与内力 弯矩、剪力) 应力与内力( 梁截面上应力与内力(弯矩、剪力)的关 应力与应力函数间微分关系, 系、应力与应力函数间微分关系,来推断 应力函数的分离变量形式。 应力函数的分离变量形式。 梁截面上的应力内力的关系: 梁截面上的应力内力的关系:
θ
M = Py = P⋅rsinθ
由材料力学初等理论,可知截面上正应力 由材料力学初等理论, 由此假定: 由此假定:
σθ ∝M(= P⋅rsinθ)
有限元经典PPT第4章

Pii Kiiui
Ki1u1 Ki2u2 Kiiui K u i,i1 i1
ui
n
Kiiui Kiiui
Kiju j
4.1.2 平面应力问题有限元的基本思想和瑞雷-里兹法
v3 f3y
3
u3
f3x
f1y v1 u1
1 f1x
v2 f2y u2
2 f2x
给定一个三角形单元和作用在角点上 的六个力,要求得六个角点的位移。 或者是要求三角形角点发生指定的位 移,在三角形三个角点如何加力?
很显然,问题的精确解很困难。采用 瑞雷-里兹法求近似式解
e号单元的三个节点I,j,k的力对应的 力的平衡方程是第2i-1,2i;2j-1,2j;2k1,2k个平衡方程
e号单元的三个节点I,j,k的位移是第 2i-1,2i;2j-1,2j;2k-1,2k个未知数
弹性模量:E 横截面积:A
1
1 L
2
2L
3
局部系单元刚度阵:
k
1
EA L
1 -1
-1
1
2 集成总刚:
0 1
解得:
ux uy
L EA
3.8284L
EA
i
j
第一类位移条件:
Ki1u1 Ki2u2 Kiiui Ki1ui1
ui 0
令: Kij 0 i j
m
vi 0
Kii 1
um 0
Pi 0
ui 0
第二类位移条件:um um
大数
充大数法: Kii Kii
第一步:求转换矩阵
k2
EA 1 2L -1
-1
1
P
cos 0
T sin
第4章_有限差分法.弹性力学

y xy
0
0
2 1 [(5 7 ) ( 6 8 )] xy 0 4h 2
可见,用差分法解平面 问题,共有两大任务:
一、建立差分方程 将(1-6~8)代入双调和 方程
4 4 4 2 2 2 4 0 4 x x y y
在结点3,x=x0-h, 在结点1, x=x0+h,代入(b) 得:
h2 2 f f f3 f 0 h 2 x 0 2 x 0
h2 2 f f f1 f 0 h 2 x 0 2 x 0
A A
B
(2 5)
从图易看出,式(2-3)右 边的积分式表示A与B之间的x方 向的面力之和;式(2-4)右边 的积分式表示A与B之间的y方向 的面力之和;式(2-5)右边的 积分式表示A与B之间的面力对 于B点的矩。
至此,我们解决了怎样 计算边界上各结点
, , x y
的值的问题。 至于边界外一行虚结点处 的值,则可用边界上
差分公式(1-1)及(1-3)是以相隔2h的两 结点处的函数值来表示中间结点处的一阶导数 值,可称为中点导数公式。 以相邻三结点处的函数值来表示一个端点 处的一阶导数值,可称为端点导数公式。 应当指出:中点导数公式与端点导数公式 相比,精度较高。因为前者反映了结点两边的 函数变化,而后者却只反映了结点一边的函数 变化。因此,我们总是尽可能应用前者,而只 有在无法应用前者时才不得不应用后者。
将式(b),(c)代入,整理得:
B B B A ( x B x A ) ( y B y A ) ( y B y ) p x d s ( x x B ) p y d s (d ) y A A x A A A , , , 为已知, 由式(d)及式(c)可见,设
弹性力学 第四章 平面问题的极坐标解答

s = sσ
(3) 多连体中的位移单值条件
§ 4-4 应力分量的坐标变换式
·问题的提出
工程中有些问题, 用极坐标计算方便, 但应力分量用直角坐 标表述更直观. 反之也存在.
由此需要对应力分量进行坐标变换.
§ 4-4 应力分量的坐标变换式
·坐标变换
已知 σx、σy、τxy , 求 σρ、συ、τρυ?
y
fρ τ + ∂τρυdρ ρυ ∂ρ ∂συ dυ συ+ ∂σρ ∂υ σρ+ dρ ∂τυρ C ∂ρ dυ τυρ+ ∂υ
B
fυ
§ 4-1 极坐标中的平衡微分方程
·平衡微分方程
x
υ dυ ρ
Σ Fρ = 0 :
συ
A
σρ τρυ P τυρ
∂σρ σρ+ dρ (ρ+dρ)dυ - σρ ρdυ ∂ρ ∂συ dυ - συ+ dυ dρ sin ∂υ 2 + τυρ+ - συ dρ sin
Σ Fυ = 0 :
συ = ?
§ 4-4 应力分量的坐标变换式
·坐标变换
συ = ?
将两坐标系下微元体组合
τyx σy σx συ
τυρ τxy
§ 4-4 应力分量的坐标变换式
·坐标变换
已知 σx、σy、τxy , 求 σρ、συ、τρυ?
O x
υ
τyx
σy σx
συ y
τυρ τxy
Σ Fυ = 0 :
O h/2 h/2 lqx源自(v)x=0, l = 0
应力边界条件: ( σy ) y=-h/2 = - q (τyx ) y=-h/2 = 0 ( σy ) y= h/2 = 0 (τyx ) y= h/2 = 0
弹性力学徐芝纶版第4章

第四章 平面问题的极坐标解答
比 较
1 f 0。 (a) x yx fx 0 x y
式(a)中第一、二、 四项与直角坐标的方 方程相似; 而第三项 分别由PB、AC面不 相等和PA、BC面不平 行引起。
为边界上已知的面力分量。
常数,或 常数,
故边界条件形式简单。
第四章 平面问题的极坐标解答
例:写出应力边界条件(集中力偶作用处为小边界) b a 0 d 0 0 a 0
a
b 0 b 0
0
q
2
2
0
0 0
l
q0
0
0
0
0
2
2
0
0
第四章 平面问题的极坐标解答
PB线应变 PB PB PC PB (ρ u ρ ) d υ ρ d υ u ρ ευ PB PB ρdυ ρ
第四章 平面问题的极坐标解答
PA转角 0, PB转角
O
d
x P
d
P
A
u
B C
u
A
u
d
y B ( u )d CB β tan u PC u d u ρ u ρ (u d υ) u ρ dυ υ υ ∴切应变为 ρ u d υ d
1.只有径向位移 u ,求形变。 P,A,B 变形后为 P', A', B', 在小变形假定 O x β 1 下, u P d d
第四章:物理方程和弹性力学解的基本方程与求法

x z y x yz
( xy yz zx ) 2 2 y
y z x y zx
( yz zx xy ) 2 2 z
z x y z xy
统一写成:
ij'kl kl'ij ik ' jl jl'ik 0
第四章:物理方程和弹性力学解的基本方程与求法
4.2 弹性力学基本方程
6个应力: x、 y、 z、 xy、 yz、 zx 6个应变: x、 y、 z、 xy、 yz、 zx
3个位移: u、v、w
虽然方程数与未知函数数量相等。但这些方程都是偏微分方程, 求解时必须需要边界条件。
将(A)前三式相加:
x
y
z
1 E
[ x
y
z
2( x
y
z )]
1 2
E
(
x
y
z
)
记: x y z ii ui'i , x y z ii
则: 或:
1 2 E
3K
ii
1 2 E
ii
0
1 2
E
m
K E 体积弹性模量。
3(1 2)
第四章:物理方程和弹性力学解的基本方程与求法
第四章:物理方程和弹性力学解的基本方程与求法 4.1 广义虎克定律
三个平衡方程中有6个应力未知函数;几何方程中6个方程有三 个位移和6个应变未知函数。15个未知函数只有9个方程,差6个 方程。为了问题可解,需要建立6个应力与应变的关系方程:
ij fij ( )
此方程在固体力学中称物理方程,或本构方程,其形式和系数
统一写成: ij 2Gij ij kk 2Gij ij
第四章:物理方程和弹性力学解的基本方程与求法
( xy yz zx ) 2 2 y
y z x y zx
( yz zx xy ) 2 2 z
z x y z xy
统一写成:
ij'kl kl'ij ik ' jl jl'ik 0
第四章:物理方程和弹性力学解的基本方程与求法
4.2 弹性力学基本方程
6个应力: x、 y、 z、 xy、 yz、 zx 6个应变: x、 y、 z、 xy、 yz、 zx
3个位移: u、v、w
虽然方程数与未知函数数量相等。但这些方程都是偏微分方程, 求解时必须需要边界条件。
将(A)前三式相加:
x
y
z
1 E
[ x
y
z
2( x
y
z )]
1 2
E
(
x
y
z
)
记: x y z ii ui'i , x y z ii
则: 或:
1 2 E
3K
ii
1 2 E
ii
0
1 2
E
m
K E 体积弹性模量。
3(1 2)
第四章:物理方程和弹性力学解的基本方程与求法
第四章:物理方程和弹性力学解的基本方程与求法 4.1 广义虎克定律
三个平衡方程中有6个应力未知函数;几何方程中6个方程有三 个位移和6个应变未知函数。15个未知函数只有9个方程,差6个 方程。为了问题可解,需要建立6个应力与应变的关系方程:
ij fij ( )
此方程在固体力学中称物理方程,或本构方程,其形式和系数
统一写成: ij 2Gij ij kk 2Gij ij
第四章:物理方程和弹性力学解的基本方程与求法
弹性力学-第4章 8 压力隧洞

一、分别对圆筒、无限大弹性体列出应力分量与位移 表达式,注意到其材料的不同,E、,A、C均不同
1、圆筒:
r
A r2
2C ,
A r2
2C
r r 0
注意在平面应变问题(4-12)公式中的位移表达为:
ur
1
E
[21 2cr
A] r
I
cos
k sin
2、无限大弹性体
r
A r2
2C,
A r2
在接触面上,无论取何值,上式均要求成立,则必 各系数相等,即:
1 [21 2cb A] 1 [21 2cb A]
E
b E
b
I I , k k
由第一式整理:
E1 E1
[21
2 Cb
A b2
]
A b2
0
令
n
E1 E1
(3)
联立(1)、(2)、(3)求出A、C、A'并代 回应力分量表达式:
b2
b2
径向位移等值: ur |rb ur |rb
ur
|r b
1
E
[21
2 cb
A] b
I
cos
k
sin
ur
|rb
1 [21
E
2cb
A] b
I cos
k s in
代入:
1 [21 2 cb A] I cos k sin
E
1
[21
b
2cb
A]
I c os
k s in
E
b
2C
r r 0
ur
1 [21
E
2cr
A] r
I cos
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
图(a)
o x
图 (b )
x
y
τ
u
s xy
=u = 0 = f y= 0
y
(σ x ) s = f x 0 v s = v = 0
例1 如图所示,试写出其边界条件。 如图所示,试写出其边界条件。
(1) (2)
(3)
∂v us = 0 ∂u = 0, = 0 x = 0, ∂y ∂x h vs = 0 h x x = a, l = 1, m = 0 a X = 0,Y = 0 y (4) y = +h, l = 0, m = +1 l(σ x )s + m(τ xy )s = X X = 0,Y = 0 m(σ y )s + l(τ xy )s = Y (σ x )s ⋅ 0 + (τ xy )s ⋅ (+1) = 0 (σ x )s = 0, (τ xy )s = 0 (σ y )s ⋅ (+1) + (τ xy )s ⋅ 0 = 0 l = 0, m = −1 y = −h, (σ y )s = 0, (τ xy )s = 0 X = 0,Y = q
σx εx
u
σy
σz
εz
τ xy
γ xy
τ yz γ yz
τ zx
6个应变分量(Stress Components )
εy
γ zx
和3个位移(Displacements)
v
w
虽然15个方程可解15个未知函数, 虽然15个方程可解15个未知函数,但由于求解时会产生待定函 15个方程可解15个未知函数 数;所以要想得出具体的解答还必需利用边界条件来确定待定 函数。 函数。
σ y ⋅ (−sin β ) +τ xy ⋅ (− cos β ) = γy sin β 右侧面: 右侧面: l = cosα, m = − sin α x = y tanα cosα ⋅σ x − sin α ⋅τ xy = 0
X =Y = 0
− sin α ⋅σ yx + cosα ⋅τ xy = 0
上 : (σ y ) s = − ql , (τ yx ) s = 0
三.混合边界条件
(Mixed Boundary Condition)
在一部分边界上的位移分量为已知, 1、在一部分边界上的位移分量为已知,另一部 分界上应力分量已知。 分界上应力分量已知。 在同一边界上, 2、在同一边界上,已知一个位移分量和一个应 力分量。 力分量。
∂τ zx ∂τ zy ∂σ z + + + fz = 0 ∂x ∂y ∂z
2、几何方程(Geometric Equations) )
方程, (Cauchy方程,3个) 方程
∂u εx = ∂x
∂v εy = ∂y
∂w εz = ∂z
γ xy
∂u ∂v = + ∂y ∂x
γ yz
∂v ∂w = + ∂z ∂y
例 举 :
fx = ql fy = 0
y
fx = 0, fy = ql
x
右 : (σ x ) s = − ql , (τ xy ) s = 0
fx = −ql fy = 0
左 : (σ x ) s = − ql , (τ xy ) s = 0
fx = 0, fy = −ql
下 : (σ y ) s = − ql , (τ yx ) s = 0
pi = σ ij n j
(在 Sσ 上)
τ zx ⋅ l + τ zy ⋅ m + σ z ⋅ n = p z
2)位移边界条件
(Displacement Boundary Condition)
u = u*
v = v*
w = w*
v w 注意: u * 、 * 、 * 为弹性体表面已知的位移
以上15个基本方程包含弹性力学所要研究的15个基本未知量, 以上15个基本方程包含弹性力学所要研究的15个基本未知量, 15个基本方程包含弹性力学所要研究的15个基本未知量 即6个应力分量(Stress Components ) 个应力分量
第四章 弹性力学问题的解法
Methods of Analysis for Elastic Mechanics
参
考 教 材
1)《弹性力学》(第4版,上册),徐芝纶著 2)《弹性力学与有限元法》,蒋玉川、张建海、 李章政编著
§4.1 弹性力学的基本方程
Basic Equations of Elastic Mechanics 平面应力问题
q
(σ ) ⋅ (−1) + (τ ) ⋅ 0 = q (σ ) = −q, (τ ) = 0
y s xy s
y s xy s
(σ x )s ⋅ 0 + (τ xy )s ⋅ (−1) = 0
说明: 说明:
x = 0 的边界条件,是有矛 的边界条件, 盾的。由此只能求出结果: 盾的。由此只能求出结果:
∂ 2ε y
2 ∂ 2ε z ∂ γ yz + 2 = 2 ∂z ∂y ∂y∂z
∂ ∂γ yz ∂γ zx ∂γ xy 2 = ∂x − ∂y + ∂z ∂z∂x ∂x
∂ 2ε y
∂ 2ε z ∂ 2ε x ∂ 2γ zx + 2 = 2 ∂x ∂z ∂z∂x
∂ 2ε z ∂ ∂γ yz ∂γ zx ∂γ xy 2 = ∂x + ∂y − ∂z ∂x∂y ∂x
o
x 上面:l=0,m=-1 左面: l=-1 m=0 下面:l=0,m=1 y 右面: l=1 m=0
(σ l = ±1 x)s = ± X (τ = m = 0 xy)s ±Y
(2).在上下两面 (2).在上下两面
l = 0 y ) s = ± f y (σ (τ = m = ±1 yx ) s ± f x
∂w ∂u = + ∂x ∂z
1 ε ij = (ui , j + u j ,i ) 2
γ zx
变形协调方程(Deformation Compatibility Equation )
——(Saint-Yenant方程)
∂ εx + 2 = 2 ∂y ∂x ∂x∂y
2
∂ 2ε y
∂ 2γ xy
∂ 2ε x ∂ ∂γ yz ∂γ zx ∂γ xy 2 = − ∂x + ∂y + ∂z ∂y∂z ∂x
.在边界上,应力分量的边界值等于对应的面力分量, 注: A.在边界上,应力分量的边界值等于对应的面力分量,且当 边界的外法线沿坐标轴正向时,两者正负号相同, 边界的外法线沿坐标轴正向时,两者正负号相同,当边界的 外法线沿坐标轴负向时,两者正负号相反。 外法线沿坐标轴负向时,两者正负号相反。
B.边界上的面力转变为应力分量其正负号规定:正面正向、 边界上的面力转变为应力分量其正负号规定:正面正向、 边界上的面力转变为应力分量其正负号规定 负面负向为正,其余为负。 负面负向为正,其余为负。
y
τ xy
(2) BC段(x = l): l 段 ):
x σ y y=0 = p(x) = p0 l
= 1, m = 0
y=0
=0
(3) AC段(y =x tan β): 段 )
l = cos( N, x) = cos(90o + β ) = −sin β
m = cos( N, y) = cos β
u = 0, v = 0.
例2 如图所示,试写出其边界条件。 如图所示,试写出其边界条件。
(1) AB段(y = 0): l = 0, m = −1 段 ): 代入边界条件公式, 代入边界条件公式,有 p(x) A
x X = 0,Y = − p(x) = − p0 l
β
N β l C
B
p 0
x
h
σ x ⋅ 0 +τ xy ⋅ (−1) = 0 σ y ⋅ (−1) +τ yx ⋅ 0 = − p(x)
1、平衡微分方程(Differential Equations of Equilibrium) ) 方程, (Navier方程,3个) 方程
∂σ x ∂τ xy ∂τ xz + + + fx = 0 ∂x ∂y ∂z
∂τ yx ∂x
+
∂σ y ∂y
+
∂τ yz ∂z
+ fy = 0
σ ij , j + f i = 0
u |x=l = 0, v |x=l = 0
∂u ∂v = 0, =0 ∂y x=l ∂x x=l
σ x ⋅ (−sin β ) +τ xy ⋅ cos β = 0 σ y ⋅ cos β +τ yx ⋅ (−sin β ) = 0
例3 图示水坝,试写出其边界条件。 图示水坝,试写出其边界条件。
左侧面: 左侧面: l = − cos β , m = − sin β
x = − y tan β
X = γy cos β Y = γy sin β
由应力边界条件公式, 由应力边界条件公式,有
α
l(σ x )s + m(τ xy )s = X m(σ y )s + l(τ xy )s = Y
σ x ⋅ (− cos β ) +τ xy ⋅ (−sin β ) = γy cos β
u、v — 边界上坐标的已知函数或边界上已知的位移分量。 边界上坐标的已知函数或边界上已知的位移分量。
二、应力边界条件
(Stress Boundary Condition)
——应力分量与面力分量之间的关系 应力分量与面力分量之间的关系 在全部边界上应力边界条件已知。 在全部边界上应力边界条件已知。
o x
图 (b )
x
y
τ
u
s xy
=u = 0 = f y= 0
y
(σ x ) s = f x 0 v s = v = 0
例1 如图所示,试写出其边界条件。 如图所示,试写出其边界条件。
(1) (2)
(3)
∂v us = 0 ∂u = 0, = 0 x = 0, ∂y ∂x h vs = 0 h x x = a, l = 1, m = 0 a X = 0,Y = 0 y (4) y = +h, l = 0, m = +1 l(σ x )s + m(τ xy )s = X X = 0,Y = 0 m(σ y )s + l(τ xy )s = Y (σ x )s ⋅ 0 + (τ xy )s ⋅ (+1) = 0 (σ x )s = 0, (τ xy )s = 0 (σ y )s ⋅ (+1) + (τ xy )s ⋅ 0 = 0 l = 0, m = −1 y = −h, (σ y )s = 0, (τ xy )s = 0 X = 0,Y = q
σx εx
u
σy
σz
εz
τ xy
γ xy
τ yz γ yz
τ zx
6个应变分量(Stress Components )
εy
γ zx
和3个位移(Displacements)
v
w
虽然15个方程可解15个未知函数, 虽然15个方程可解15个未知函数,但由于求解时会产生待定函 15个方程可解15个未知函数 数;所以要想得出具体的解答还必需利用边界条件来确定待定 函数。 函数。
σ y ⋅ (−sin β ) +τ xy ⋅ (− cos β ) = γy sin β 右侧面: 右侧面: l = cosα, m = − sin α x = y tanα cosα ⋅σ x − sin α ⋅τ xy = 0
X =Y = 0
− sin α ⋅σ yx + cosα ⋅τ xy = 0
上 : (σ y ) s = − ql , (τ yx ) s = 0
三.混合边界条件
(Mixed Boundary Condition)
在一部分边界上的位移分量为已知, 1、在一部分边界上的位移分量为已知,另一部 分界上应力分量已知。 分界上应力分量已知。 在同一边界上, 2、在同一边界上,已知一个位移分量和一个应 力分量。 力分量。
∂τ zx ∂τ zy ∂σ z + + + fz = 0 ∂x ∂y ∂z
2、几何方程(Geometric Equations) )
方程, (Cauchy方程,3个) 方程
∂u εx = ∂x
∂v εy = ∂y
∂w εz = ∂z
γ xy
∂u ∂v = + ∂y ∂x
γ yz
∂v ∂w = + ∂z ∂y
例 举 :
fx = ql fy = 0
y
fx = 0, fy = ql
x
右 : (σ x ) s = − ql , (τ xy ) s = 0
fx = −ql fy = 0
左 : (σ x ) s = − ql , (τ xy ) s = 0
fx = 0, fy = −ql
下 : (σ y ) s = − ql , (τ yx ) s = 0
pi = σ ij n j
(在 Sσ 上)
τ zx ⋅ l + τ zy ⋅ m + σ z ⋅ n = p z
2)位移边界条件
(Displacement Boundary Condition)
u = u*
v = v*
w = w*
v w 注意: u * 、 * 、 * 为弹性体表面已知的位移
以上15个基本方程包含弹性力学所要研究的15个基本未知量, 以上15个基本方程包含弹性力学所要研究的15个基本未知量, 15个基本方程包含弹性力学所要研究的15个基本未知量 即6个应力分量(Stress Components ) 个应力分量
第四章 弹性力学问题的解法
Methods of Analysis for Elastic Mechanics
参
考 教 材
1)《弹性力学》(第4版,上册),徐芝纶著 2)《弹性力学与有限元法》,蒋玉川、张建海、 李章政编著
§4.1 弹性力学的基本方程
Basic Equations of Elastic Mechanics 平面应力问题
q
(σ ) ⋅ (−1) + (τ ) ⋅ 0 = q (σ ) = −q, (τ ) = 0
y s xy s
y s xy s
(σ x )s ⋅ 0 + (τ xy )s ⋅ (−1) = 0
说明: 说明:
x = 0 的边界条件,是有矛 的边界条件, 盾的。由此只能求出结果: 盾的。由此只能求出结果:
∂ 2ε y
2 ∂ 2ε z ∂ γ yz + 2 = 2 ∂z ∂y ∂y∂z
∂ ∂γ yz ∂γ zx ∂γ xy 2 = ∂x − ∂y + ∂z ∂z∂x ∂x
∂ 2ε y
∂ 2ε z ∂ 2ε x ∂ 2γ zx + 2 = 2 ∂x ∂z ∂z∂x
∂ 2ε z ∂ ∂γ yz ∂γ zx ∂γ xy 2 = ∂x + ∂y − ∂z ∂x∂y ∂x
o
x 上面:l=0,m=-1 左面: l=-1 m=0 下面:l=0,m=1 y 右面: l=1 m=0
(σ l = ±1 x)s = ± X (τ = m = 0 xy)s ±Y
(2).在上下两面 (2).在上下两面
l = 0 y ) s = ± f y (σ (τ = m = ±1 yx ) s ± f x
∂w ∂u = + ∂x ∂z
1 ε ij = (ui , j + u j ,i ) 2
γ zx
变形协调方程(Deformation Compatibility Equation )
——(Saint-Yenant方程)
∂ εx + 2 = 2 ∂y ∂x ∂x∂y
2
∂ 2ε y
∂ 2γ xy
∂ 2ε x ∂ ∂γ yz ∂γ zx ∂γ xy 2 = − ∂x + ∂y + ∂z ∂y∂z ∂x
.在边界上,应力分量的边界值等于对应的面力分量, 注: A.在边界上,应力分量的边界值等于对应的面力分量,且当 边界的外法线沿坐标轴正向时,两者正负号相同, 边界的外法线沿坐标轴正向时,两者正负号相同,当边界的 外法线沿坐标轴负向时,两者正负号相反。 外法线沿坐标轴负向时,两者正负号相反。
B.边界上的面力转变为应力分量其正负号规定:正面正向、 边界上的面力转变为应力分量其正负号规定:正面正向、 边界上的面力转变为应力分量其正负号规定 负面负向为正,其余为负。 负面负向为正,其余为负。
y
τ xy
(2) BC段(x = l): l 段 ):
x σ y y=0 = p(x) = p0 l
= 1, m = 0
y=0
=0
(3) AC段(y =x tan β): 段 )
l = cos( N, x) = cos(90o + β ) = −sin β
m = cos( N, y) = cos β
u = 0, v = 0.
例2 如图所示,试写出其边界条件。 如图所示,试写出其边界条件。
(1) AB段(y = 0): l = 0, m = −1 段 ): 代入边界条件公式, 代入边界条件公式,有 p(x) A
x X = 0,Y = − p(x) = − p0 l
β
N β l C
B
p 0
x
h
σ x ⋅ 0 +τ xy ⋅ (−1) = 0 σ y ⋅ (−1) +τ yx ⋅ 0 = − p(x)
1、平衡微分方程(Differential Equations of Equilibrium) ) 方程, (Navier方程,3个) 方程
∂σ x ∂τ xy ∂τ xz + + + fx = 0 ∂x ∂y ∂z
∂τ yx ∂x
+
∂σ y ∂y
+
∂τ yz ∂z
+ fy = 0
σ ij , j + f i = 0
u |x=l = 0, v |x=l = 0
∂u ∂v = 0, =0 ∂y x=l ∂x x=l
σ x ⋅ (−sin β ) +τ xy ⋅ cos β = 0 σ y ⋅ cos β +τ yx ⋅ (−sin β ) = 0
例3 图示水坝,试写出其边界条件。 图示水坝,试写出其边界条件。
左侧面: 左侧面: l = − cos β , m = − sin β
x = − y tan β
X = γy cos β Y = γy sin β
由应力边界条件公式, 由应力边界条件公式,有
α
l(σ x )s + m(τ xy )s = X m(σ y )s + l(τ xy )s = Y
σ x ⋅ (− cos β ) +τ xy ⋅ (−sin β ) = γy cos β
u、v — 边界上坐标的已知函数或边界上已知的位移分量。 边界上坐标的已知函数或边界上已知的位移分量。
二、应力边界条件
(Stress Boundary Condition)
——应力分量与面力分量之间的关系 应力分量与面力分量之间的关系 在全部边界上应力边界条件已知。 在全部边界上应力边界条件已知。