建立方程定解条件

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

数是唯一确定的。
若质量以密度( x, y, z)分布在区域上,则 上的质 量所产生的总引力位势应为:
( x, y, z) (,,)
ddd
( x )2 ( y )2 (z )2
直接计算可得: 0 当( x, y, z) 时
还可进一步验证: -4 当( x, y, z) 时
上页 下页 返回
(1.1)或(1.2)
上页 下页 返回
定义 调和方程(1.1)的连续解,即,关于变量 x, y,z 具有二阶连续偏导数并且满足方程 (1.1)的连续函 数称为调和函数。
(1) 引力位势
P0 ( x0 , y0 , z0 )
F
位于 P0 处,质量为 M 的质点对于位于 P 处具有单位质量的质点的引力 F,
变形所产生的面积的增量的乘积。即
应变能
T
T
1
v
2 x
v
2 y
1
dxdy
由于
1
1
1 2
o()
T 2
v
2 x
v
2 y
dxdy
上页
下页
返回
即:
应变能 T 2
(v
2 x
v
2 y
)dxdy
如果膜所受的垂直方向的外力有两个,一个为作用
在膜内的F ( x, y)(牛顿∕米 2),另一个是作用在膜的边
从而静电场的电势u应当满足泊松方程
u 4
如果静电场的某一区域里没有电荷,即ρ=0,则 静电场方程在该区域上简化为拉普拉斯方程
u 0
上页 下页 返回
(3) 稳定温度分布 热传导持续进行下去,如果达到平衡状态,温度
的空间分布不再变动,即ut 0。以ut 0代入三维热 传导方程,得到温度的稳定分布方程
limu( x, y, z) 0 (r x2 y2 z2 )
r
上页 下页 返回
其它边界条件
(5) 第三类边界条件
u u g (其中n是 的单位外法向量)
n
(6) 等值面边界条件
u
待定常数C
(总流量边界条件)
u dS n
已知常数A
对于泊松方程 (1.2)的各种边值问题,只要找出泊松方程
u
2u x 2
2u y 2
2u z 2
0
上页 下页 返回
2.定解条件和定解问题
由于 u 与时间 无关,只是空间 变量的函数,所以 定解条件中只有边界条件,称为边值问题。本章主要 研究第一、二类边值问题及外问题。
(1) 第一边值问题(Dirichlet问题)
u | g (其中 g是在闭曲面 上给定的函数) (2) 第二边值问题(Neumann问题)
1.方程的导出
本章研究调和方程(又称拉普拉斯方程)
u
2u x 2
2u y 2
2u z 2
0
(1.1)
以及泊松方程
u
2u x 2
2u y 2
2u z 2
f (x, y, z)
(1.2)
的基本定解问题及解的性质。
方程(1.1)及(1.2)具有广泛的物理背景,振动趋于
平衡,热传导趋于稳定以及保守场都可归结为方程
取膜的平衡位置为 xoy平面上的区域 ,以 u轴垂 直于 xoy平面,并与 x, y 轴组成右手系。在的边界 上, ,在 上已知膜的位移为,在上膜受到 外力的作用,设它的垂直于膜的分量为 p( x, y)。
上页 下页 返回
② 物理原理:从力学上讲,可以从不同角度来刻画
这个平衡状态。例如,力的平衡原理、虚功原理等。这
里用最小总位能原理,即
受外力作用的弹性体,在满足已知边界位移约束的
一切可能位移中,以达到平衡状态的位移使物体的总位
能为最小。
③ 数学形式:为了对膜的平衡问题写出上述原理
的数学形式,我们必须弄清楚两个概念。
1、什么叫总位能?对于膜来说,总位能的数学形式
是什么? 2、什么叫“满足已知边界位移约束的一切可能位
(2) 静电场的电位势 从静电学知道,穿过闭合曲面 向外的电通量等
于闭合曲面 所围空间 G 中的电量的 4 倍,即
E ndS 4 dxdydz
G
应用高斯公式,上式可改写为:
Edxdydz 4 dxdydz
G
G
上页 下页 返回
由区域G的任意性得: E 4
静电场方程
由于静电场是无旋场,因而存在电势u, E u
大小为 M ,方向与 r2
PP0 同向。即
P(x, y,z)
F(
x,
y,
z)
M r2
பைடு நூலகம்
x
r
x0
,
y
r
y0
,
z
r
z0
其中 r ( x x0 )2 ( y y0 )2 (z z0 )2
上页 下页 返回
令 ( x, y, z) M r
经计算可得: F grad
称为引力位势。除了可以相差一个常数外,位势函
的一个特解,由叠加原理,就能化为调和方程 (1.1)的对应的
边值问题。由引力位势的方程可知,当 g 满足 Hölder 条件时,
这种特解是容易找到的。所以以后主要研究调和方程 (1.1)的
边值问题。
上页 下页 返回
3.变分原理
膜的平衡问题:
① 物理模型:考虑一处于紧张状态的薄膜,它的部 分边界固定在一框架上,另一部分边界上受到外力的作 用;若整个薄膜在垂直于平衡位置的外力作用下处于平 衡状态,问膜的形状如何?
界上的 p( x, y)(牛顿∕米),在它们的作用下,膜上各
点的位移为v( x, y),则
外力作功 F ( x, y)v( x, y)dxdy p(s)v(s)ds
从而,当膜上各点的位移为v( x, y)时,总势能J (v)为
J(v)
T 2
(v
2 x
v
2 y
)dxdy
F ( x, y)v( x, y)dxdy
p(s)v(s)ds
上页 下页 返回
问题2的解答: 所有可能位移即表示这样一个函数集合:1o 满足
已知位移约束,即v ;2o 使总势能J (v)有意义。
u g (其中n是 的单位外法向量) n
上页 下页 返回
(3) Dirichlet外问题
u 0 ( x, y, z) 外部 u g
(4) Neumann外问题
u 0 ( x, y, z) 外部
u n
g
(其中n是 的单位内法向量)
注:当考虑外问题时,为保证解的唯一性,还需对解在 无穷远的状况加以限制。在三维情形,通常要求:
移”?
上页 下页 返回
问题 1 的解答: 按照弹性力学理论
总位能 = 应变能 - 外力作功
处于某一位置的膜所具有的应变能就是把膜从水平
位置转移到这个位置,为了抵抗张力所作的功的总和。
从弹性力学理论知道:假设膜的形状为u v( x, y),
则当 vx , v y 1时,膜的应变能可以表为张力与膜由于
相关文档
最新文档