非齐次边界条件问题(10.30)

非齐次边界条件问题(10.30)
非齐次边界条件问题(10.30)

(10.30)非齐次边界条件问题

问题1

, (,0)(), (0,)0, (,)(0)t xx u ku u x f x u t u l t A A ====≠

求解非齐次边界问题时,首先应将其转化为齐次边界问题。因此,此处首先找出方程的稳态解,即与时间t 无关的解0()u x ,将其代入原方程后可得

[][]00()0()t xx u x k u x ==

解得

0()u x px q =+

式中,p 、q 为待定系数。根据边界条件可得

0(0)0u q ==

0()u l A pl q ==+

解得

, 0A p q l

== 所以

0()A u x x l

=

构造函数 0(,)()(,)u x t u x v x t =+

代入原方程可得

[][]00()()t t xx t xx u u x v k u x kv =+=+

化简后可得

t xx v kv =

又由初始条件可得

0()(,0)()(,0)f x u x u x v x ==+

所以

0(,0)()()v x f x u x =-

由边界条件还可以得到

(0,)(,)0v t v l t ==

因此,题设问题就转化为了齐次边界条件问题,即求解

0, (,0)()(), (0,)(,)0t xx v kv v x f x u x v t v l t ==-==

由变量分离法,首先假设

(,)()()v x t X x T t =

进而有

()'()"()()X x T t kX x T t =

移项整理得

''()'()()()

X x T t A X x kT t =≡ 其中A 是与x , t 都无关的常数,于是有

'()()T t AkT t =

"()()X x AX x =

分别求解,对于()T t

d ()d ()T t Ak t

T t =??

所以

0()Akt

T t C e =

对于()X x ,当0A ≥时,都可以得到()0f x ≡,与题设不符。当0A <时,此微分方程所

对应的特征根方程有两个共轭的复根,因此其通解为

012()()X x e C C =+

由边界条件有

(0)()0, ()()0X T t X l T t ==

因为()T t 不恒为0,所以有

(0)0, ()0X X l ==

将上式代入()X x 的通解可得

10C ==

因此

22

2, n n A l ππ==-

其中n 为整数。所以

2()sin n X x C x l π=

2202()exp n k T t C t l π??=- ???

进而有

222(,)sin exp n n n n k v x t C x t l l ππ??=- ???

由上式可知,对于每个不同的n ,都有一个(,)n v x t 与之对应。而0(,0)()()v x f x u x =-,而()f x 和0()u x 都是给定的函数,所以此处得到的(,)n v x t 仍然不是方程的解。

对非齐次偏微分方程的求解齐次边界条件下非齐次发展方程的混合问题

对非齐次偏微分方程的求解 齐次边界条件下非齐次发展方程的混合问题 (一)冲量定理法 (二)傅立叶级数法 齐次边界条件下非齐次场位方程的混合问题 (一)方程和边界条件同时齐次化 非齐次方程的求解思路 ? 用分解原理得出对应的齐次问题 ? 解出齐次问题 ? 求出任意非齐次特解 ? 叠加成非齐次解 方法一 冲量定理法 前提条件:除了方程为非齐次的外,其它定解条件都是齐次的(初始条件均取零值)。 基本思路:利用叠加原理将受迫振动的问题转化为(无穷多个)自由振动问题的叠加. 2000(,)0,0 (),() tt xx x x l t t t u a u f x t u u u x u x φψ====?-=?? ==??==?? 试设 12u u u =+ 222112211 11,(,0)(,0)(),(),(0,)0,(,)0 u u a t x u x u x x x t u t u l t ?ψ???=???? ?? ==??? ==?, ()22 222 2222 22,,(,0) (,0)0,0,(0,)0,(,)0 u u a f x t t x u x u x t u t u l t ???-=??????==??? ==?. 物理意义: 在时间 0 — t ,可以把非齐次项(单位质量所受的持续作用力)看成许多前后相继(无穷多个)的“瞬时”力引起的物理过程的线性叠加。

22 222 0,0,(,),0,0t t t x x l a t t x f x d τ τωωτωω ττωω====???=>?????==??==?22 222 ,0,(,),0,0 t t t x x l v v a t t x v v f x v v τ τττ====???=>?????==??==? 相应的,我们也可以把位移(,)u x t 也表示为 20(,)(,;)d t u x t v x t ττ=?, 则(,;)v x t d ττ就应当是瞬时力所产生的位移.更进一步说,(,,)v x t τ就是定解问题 22 222 0,0,(,),0,0 t t t x x l a t t x f x d τ τωωτωω ττωω====???=>?????==??==? 22 222 ,0,(,),0,0 t t t x x l v v a t t x v v f x v v τ τττ====???=>?????==??==? 的解.非齐次项只存在于τ时刻,其全部效果只是使得弦在τ时刻获得一个瞬时速度. 那么由偏微分方程的积分 22 0002 2 20 00(,)()v v dt a dt f x t d t x ττττ τττδττ+++---??-=-????? 推导出 (,,) (,)t v x t f x t τττ=+?=? 令 1t t τ=- 则定解问题就可以写成这种形式(0t τ=+简写成t τ=) 11122 22210 00,0,(,), 0,0 t t t x x l v v a t t x v v f x v v ττ====???=>?????==?? ==? 在运算过程中,十分需要注意的是,瞬时力的重复计算,不能把瞬时力既算入定解方程的其次项,又算入初速度! 总结一下,在上面的过程中,冲量定理就把求解非齐次方程、齐次边界条件以及齐次初条件的定解问题转化成了对齐次方程、齐次边界条件的定解问题的

椭圆型方程非齐次边值问题的变分形式 4..

目录 1引言 2椭圆型方程非齐次第一边值问题的变分形式2.1建立第一边值条件等价极小位能原理2.2建立第一边值条件等价的虚功原理 3椭圆型方程非齐次第二边值问题的变分形式3.1建立第二边值条件的极小位能原理 3.2建立第二边值条件的虚功原理 4椭圆型方程非齐次第三边值问题的变分形式4.1建立第三边值条件的极小位能原理 4.2建立第三边值条件的虚功原理

椭圆型方程非齐次边值问题的变分形式 1引言 很多实际问题的微分方程是通过泛函的变分得到的, 在变分过程中增加了未知函数导数的阶数. 反之某些变分方程的定解问题可通过构造相应的泛函, 使求泛函的极小值与求解微分方程的定解问题等价也就是说, 变分法最终寻求的是极值函数, 它们使得泛函取得极大或极小值. 变分原理在物理学中, 尤其是力学中有着广泛运用, 如著名的虚功原理、极小位能原理、余能原理和哈密顿原理等, 几乎所有的自然定律都能用变分原理的形式予以表达. 在当代变分已成为有限元法的理论基础,是求解边值问题的强力工具. 2椭圆型方程第一边值问题的变分形式 椭圆型方程第一边值问题: G u G y x f u v k =∈=+??-Γ)2.1(,),(,)(σ, 其中Γ是边界, G 是平面区域 ).()()(), (),(,0),(, 0min ),(),(21y u k y x u k x u k C g G L f G C G c y x k k G ????+????= ??Γ∈∈≥∈>∈=σσ 定义:{} ),(,)(),()(2 21b a I I L f I L f f I H =∈'∈= 在解决第一边值问题的变分形式的过程中, 我们先运用格林第一公式和极小位能原理建立等价的变分形式, 再运用虚功原理建立等价的变分形式.为此我们需要考虑如下结果: 极小位能原理, 虚功原理, 格林第一公式. 格林第一公式:G 是xy 平面上的一有界区域,其边界Γ为分段的光滑曲线,n 为曲线Γ的单位外法 向量, n u ??是u 沿n 的方向导数,则有: .)( )(vds n u dxdy y v y u x v x u xdy vd u G G ???Γ??-????+????=?-

非齐次边界条件问题(10.30)

(10.30)非齐次边界条件问题 问题1 , (,0)(), (0,)0, (,)(0)t xx u ku u x f x u t u l t A A ====≠ 求解非齐次边界问题时,首先应将其转化为齐次边界问题。因此,此处首先找出方程的稳态解,即与时间t 无关的解0()u x ,将其代入原方程后可得 [][]00()0()t xx u x k u x == 解得 0()u x px q =+ 式中,p 、q 为待定系数。根据边界条件可得 0(0)0u q == 0()u l A pl q ==+ 解得 , 0A p q l == 所以 0()A u x x l = 构造函数 0(,)()(,)u x t u x v x t =+ 代入原方程可得 [][]00()()t t xx t xx u u x v k u x kv =+=+

化简后可得 t xx v kv = 又由初始条件可得 0()(,0)()(,0)f x u x u x v x ==+ 所以 0(,0)()()v x f x u x =- 由边界条件还可以得到 (0,)(,)0v t v l t == 因此,题设问题就转化为了齐次边界条件问题,即求解 0, (,0)()(), (0,)(,)0t xx v kv v x f x u x v t v l t ==-== 由变量分离法,首先假设 (,)()()v x t X x T t = 进而有 ()'()"()()X x T t kX x T t = 移项整理得 ''()'()()() X x T t A X x kT t =≡ 其中A 是与x , t 都无关的常数,于是有 '()()T t AkT t = "()()X x AX x = 分别求解,对于()T t d ()d ()T t Ak t T t =??

齐次化原理

补充阅读材料之三 齐次化原理的物理背景 以非齐次热传导方程为例,求解 (1) 2(,),0,0,(0,)0,(,)0,0,(,0)0,0.t xx u a u f x t x L t u t u L t t u x x L ?=+<<>?==>??=≤≤? 注:边界条件为齐次,初始值为。 0. 由方程推导知 (,)(,)F x t f x t c ρ =,其中c 为比热(单位质量升高单位温度所需热量),ρ为长杆的(线)密度,为热源强度(单位时间单位体积(长度)产生的热量 (,)F x t 1. 瞬时热源产生的效应: 1) 研究0时刻的热源引起的热效应:考虑短时段[,0]s Δs ,?Δ为无穷小量, 此时段上 的热源强度. 在(,0)F x ≈s ?Δ时刻将热源打开,然后在时刻将之关闭,则热源在0x 点附近产生的热量(,0)F x s x ≈ΔΔ (x Δ为无穷小量),所产生的温度分布(,F x 0)(,0)s x f x s c x ρΔΔ≈=Δ0≥Δ. 则t 时(热源关闭)的温度分布(即0时刻的热源引起的热效应),设为(,)v x t s Δ,满足齐次热传导方程 20()(),0,(0,)(,)0,, (,)(,0),0.t xx t v s a v s x L t s v t s v L t s t s v x t s f x s x L =?Δ=Δ<<>?Δ=Δ=>??Δ=Δ≤≤? , 约掉因子s Δ得到 20,0,,(0,)(,)0,, (,)(,0),0.t xx t v a v x L t s v t v L t t s v x t f x x L =?=<<>?==>??=≤≤? 2) 同样方法研究 (固定)时刻的热源引起的热效应:考虑短时段[,, .在0s ≥]s s s ?Δ(,)F x s ≈s s ?Δ时刻将热源打开,然后在时刻将之关闭,则热源在s x 点附近

相关文档
最新文档