用分离变量法解三维坐标中的拉普拉斯方程

圆柱坐标系下的分离变量法

5．圆柱坐标系下的分离变量法 5．1极坐标系下的拉普拉斯方程 考虑半径为a 的一个薄圆盘，已知圆盘内部无热源，边界温度给定，且温度分布(,)u x t 随时间演化已趋于稳定，试求此时的温度分布(,)u x y 。上述定解问题可表述为 20 u x y x y D ?=∈(,)(,) （5.1.1a ） 222 x y a u x y g x y x y +==∈∑(,) (,) (,) (5.1.1b) 其中，∑表示圆盘的边界，即222x y a +=，D 表示∑围成的内域。 对于二维平面场问题，即物理量的空间分布与z 无关，当物体边界为矩形时，采用直角坐标系比较方便。因为边界方程可方便地用直角坐标表示出来，如 0x =，x a =，0y =，y b = 但当物体边界为圆形时采用极坐标系可大为简化边界方程，从而给问题的求解带来方便。而在极坐标系下，拉普拉斯方程表示为 2 2 22 11()r r x x r ? ????=+??? （5.1.2） 从而（5.1.1）定解问题可改写成 222 110u u r r D r x x r ?? ???+=∈???()(,) (5.1.3a) r a u r g r ???==∈∑(,) () (,) (5.1.3b) 注意到定解问题（5.1.3）中的边界条件属于第Ⅰ类，通常称之为狄里克莱 ()Dirichlet 问题，也称第I 边值问题。若（5.1.3）中的边界条件是第Ⅱ类的，则 称相应的定解问题为牛曼()Newman 问题，也称第Ⅱ边值条件。若（5.1.3）中的边界条件是第Ⅲ类的，则称相应的定解问题为罗宾()Robin 问题。此外，本题研究内域

第24讲分离变量法第4章介质中的电动力学4§4拉普拉斯方程分离变量法

第24讲 分离变量法 第4章 介质中的电动力学（4） §4.4 拉普拉斯方程 分离变量法 以上两节给出静电问题的一般公式，并说明静电学的基本问题式求解满足给定边界条件的泊松方程的解。只有在界面形状是比较简单的几何曲面时，这类问题的解才能以解析形式给出，而且视具体情况不同而有不同的解法。 在许多实际问题中，静电场是由带电导体决定的。例如电容器内部的电场是由作为电极的两个导体板上所带电荷决定的；又如电子光学系统的静电透镜内部，电场是由于分布于电极上的自由电荷决定的。这些问题的特点是自由电荷只出现在一些导体的表面上，在空间中没有其它自由电荷分布。因此，如果我们选择这些导体表面作为区域V 的边界，则在V 内部自由电荷密度 ρ = 0 ，因而泊松方程化为比较简单的拉普拉斯（Laplace ）方程 20??= （4.4---1） 产生这电场的电荷都分布于区域V 的边界上，它们的作用通过边界条件反映出来。因此，这类问题的解法是求拉普拉斯方程的满足边界条件的解。 （4.4---1）式的通解可以用分离变量法求出。先根据界面形状选择适当的坐标系，然后在该坐标系中由分离变量法解拉普拉斯方程。最常用的坐标系有球坐标系和柱坐标系。这里我们写出用球坐标系得出的通解形式（见附录Ⅱ）。球坐标用（R ，θ，φ）表示，R 为半径，θ为极角，φ为方位角。拉氏方程在球坐标系中的通解为 1 .(,,)()(cos )cos n m nm nm n n n m b R a R P m R ?θφθφ+=+ ∑ 1 ,()(cos )sin n m nm nm n n n m d c R P m R θφ+++ ∑ （4.4---2） 式中 a n m ,b n m ,c n m 和 d n m 为任意常数，在具体问题中有边界条件定出。 P m n （cos θ） 为缔和勒让德（Legendre ）函数。若该问题中具有对称轴，取此轴为

第2章 定解问题

第2章定解问题 1、何谓数理方程？按其描绘的物理过程，它可分为哪几类？ 2、何谓定解问题？它分为哪几类？试写出一维波动方程的Cauchy问题的数学表示。 3、何谓定解条件？它包括哪些内容？ 4、何谓边界条件？它分为哪几类？一个边界需用几个边界条件来描述？ 5、用数理方程来研究物理问题需要经历哪几个步骤？ 6、在静电场问题中，由介电常数分别为和的两种介质组成的系统的交界面S 处的 衔接条件有几个？应如何表示？ 7、如何导出物理模型的数理方程？在推导弦的横振动方程时采用了哪些近似？由小角度近似我们得到什么结论？ 8、热传导方程的扩散方程有何共同和不同之处？ 9、在杆的纵振动问题中，若端自由，这个边界条件如何写？你能从Hooke定律出发证明吗？ 10、在杆的导热问题中，若端绝热，这个边界条件该如何写？你能从一物理定律出 发证明吗？ 11、在热传导问题中，若热源密度不随时间而变化，则热传导方程会 发生怎样的变化？ 12、在弦的横振动问题中，若弦受到了一与速度成正比的阻力，该阻力对于弦的振动问题

是否起到了源的作用？若受到了一与位移成正比的回复力呢？ 第3章行波法 1、行波法的解题要领是什么？它适合用来求解哪一类定解问题？为什么？ 2、一维波动方程的通解为什么含有两个任意函数？他们各个有怎样的形式和怎样的物理意义？靠什么确定他们的具体函数形式？ 3、公式是用行波法求解弦的横振动问题时推得的，能否用公式求解如下定解问题？请说明原因？ 4、能否用公式求解如下定解问题？ 5、能否用行波法求解如下定解问题？ 6、你能否根据直角坐标系中的

导出球坐标中球对称情况下的的表达式 请记住这个结论： 7、何谓平均值法？你能通过引入球面的平均值，将三维的波动方程 化为关于平均值的一维方程吗？ 8、在Poisson 公式中，?若已知 9、对于定解问题 除了可用Poisson 公式求解外？你能否有其他的求解法？ 10、在弦的横振动方程单位质量的弦所受的外 力，若将则怎样的物理含意？它的量纲是什么？ 11、冲量原理的精神是什么？ 12、你能否用纯强迫振动的解来求解定解问题

第14讲 图形与坐标

第五章拓展与提高 第15讲图形与坐标 一、学习目标 1.在直角坐标系中，能写出一个已知顶点坐标的多边形沿坐标轴方向平移后图形的顶点坐标，并知道对应顶点坐标之间的关系. 2.在直角坐标系中，能将一个多边形依次沿两个坐标轴方向平移后所得到的图形与原来的图形具有平移关系，体会图形顶点坐标的变化. 3.在直角坐标系中，以坐标轴为对称轴（或以原点为对称中心），能写出一个已知顶点坐标的多边形的对称图形的顶点坐标，并知道对应顶点坐标之间的关系. 4.在直角坐标系中，探索并了解将一个多边形的顶点坐标（有一个顶点为原点、有一个边在横坐标轴上）分别扩大或缩小相同倍数时所对应的图形与原图形是位似的. 二、基础知识·轻松学 1．用坐标表示平移 在平面直角坐标系中，将点（x，y）向右（或左）平移a个单位长度，可以得到对应点（x+a，y）或（x-a，y）；将点（x，y）向上（或下）平移b个单位，?可以得到对应点（x，?y+b）或（x，y-b）． 【精讲】点沿着x轴正方向或y轴正方向平移，相应的坐标增加；而点沿着负方向平移，则相应的坐标就减少． 2.由坐标的变化规律判断图形平移方向与平移距离 在平面直角坐标系内，如果把一个图形各个点的横坐标都加（或减去）一个正数a，相应的新图形就是把原图形向左（或向右）平移a个单位长度；如果把它各个点的纵坐标都加（或减去）一个正数a，?相应的新图形就是把原图形向上（或向下）平移a个单位长度．【精讲】图形是由无数个点组成的，所以图形的变化实质是由点和点的坐标变化引起的，我们知道：对一个图形进行平移，这个图形上所有点的坐标都要发生变化；反过来，从图形上的点的坐标的某种变化，也可看出这个图形进行了怎样的平移。故着重理解：从点的坐标变化来分析点的平移方向和平移距离．关键是看相应对应点的坐标是增加还是减少，如果横

工程力学第14讲答案详解

习题14-2图 习题14-3图 第14章 压杆的平衡稳定性分析与压杆设计 14－1 关于钢制细长压杆受力达到分叉载荷之后，还能不能继续承载，有如下四种答 案，试判断哪一种是正确的。 （A ）不能，因为载荷达到临界值时，屈曲位移将无限制地增加； （B ）能，压杆一直到折断时为止都有承载能力； （C ）能，只要横截面上的最大应力不超过一定限度； （D ）不能，因为超过分叉载荷后变形不再是弹性的。 正确答案是 C 。 14－2 图示a 、b 、c 、d 四桁架的几何尺寸、杆的横截面直径、材料、加力点及加力方向均相同。关于四桁架所能承受的最大外力F Pmax 有如下四种结论，试判断哪一种是正确的。 （A ）)d ()b ()c ()a (max P max P max P max P F F F F =<=； （B ）)d ()b ()c ()a (max P max P max P max P F F F F ===； （C ）)c ()b ()d ()a (max P max P max P max P F F F F =<=； （D ）)d ()c ()()a (max P max P max P max P F F b F F =<=。 正确答案是 A 。 14－3 图示四压杆均为圆截面直杆，杆长相同，且均为轴向加载。关于四者分叉载荷大小有四种解答，试判断哪一种是正确的（其中弹簧的刚度较大）。 （A ）)d ()c ()b ()a (Pcr Pcr Pcr Pcr F F F F <<<； （B ）)d ()c ()b ()a (Pcr Pcr Pcr Pcr F F F F >>>； （C ）)a ()d ()c ()b (Pcr Pcr Pcr Pcr F F F F >>>； （D ）)d ()c ()a ()b (Pcr Pcr Pcr Pcr F F F F >>>。 正确答案是 D 。 解：图（b ）上端有弹性支承，故其临界力比图(a)大； 图(c)下端不如图(a)刚性好，故图(c)临界力比图(a)小； 图(d)下端弹簧不如图(c)下端刚性好，故图(d)临界力比图(c)小。 14－4 一端固定、另一端弹簧侧向支承的压杆。若可采用欧拉公式2 2Pcr )/(πl EI F μ=，试确定其中长度系数的取值范围为

柱坐标及球坐标下导热微分方程的推导及分析

柱坐标及球坐标下导热微分方程的 推导及分析 哈尔滨工业大学市政学院 摘要：运用热力学第一定律，建立温度场，利用微分方程在不同坐标系的不同形式进行分析问题 关键词：柱坐标 球坐标 导热微分方程 1.柱坐标系下导热微分方程 假定所研究的物体是各向同性的连续介质，其导热率λ，比热容c 和密度ρ均为已知，并假设物体内具有内热源。用单位体积单位时间内所发出的热量 qv(w/m *3)表示内热源的强度。基于上述各项假定，再从进行导热过程的物体中分割出一个微元体，如图。根据热力学第一定律，对微元体进行热平衡分析，那么在d τ时间内导入和导出微元体的净热量，加上内热源的发热量，应等于微元体热力学能的增加，即 导入与导出微元体的净热量（Ⅰ）＋微元体内热源的发热量（Ⅱ）＝微元体中热力学能 的增加（Ⅲ） 下面分别计算式中Ⅰ、Ⅱ、Ⅲ三项： 在 dτ时间内，沿 r 轴方向： τ?λτ?λdzd rd r t d dzd rd q r r t r ??-=Φ∴-=?? τ?dzd rd q r r =Φ

τ?λdzd drd r t r r d d r dr d d dr r r r r dr r )(????=Φ-Φ∴?Φ?=Φ-Φ++ ○1 在 dτ时间内，沿 ?轴方向： τ? λ? λτ ????drdzd t r t r q drdzd q ??-=Φ∴??-==Φ11 τ??λ?? ????????dzd drd t r d d d )1(????=Φ-Φ∴?Φ?=Φ-Φ++ ○ 2 在 dτ时间内，沿 z 轴方向： τ?λλτ ?drd rd z t z t q drd rd q z z z z ??-=Φ∴??-==Φ τ?λdzd drd z t r z z dz dz z z z z dz z )(????=ΦΦ∴?Φ?=Φ-Φ+-+ ○3 将 r 、Φ、z 三个方向导入和导出微元体的净热量相加得到 : I=○1+○2+○3 在 dτ时间内，微元体中内热源的发热量为 Ⅱ=dzdr rdrd q v ? 在 dτ时间内，微元体中热力学能的增量为 Ⅲ=τ?τ ρdzd rdrd t c ?? 联立I ，III ，II 可得导热微分方程在圆柱坐标下的公式： )()(1)(12z t z t r r t r r r r q t c v ????+????+????+=??λ?λ?λτρ

《流体力学》(柱坐标系和求坐标系下)连续方程推导的巧方法

《流体力学》连续方程推导的巧方法 施春华,高庆九,李忠贤 (南京信息工程大学大气科学学院,江苏南京 210044) 摘要:针对柱坐标系和球坐标系下《流体力学》中连续方程形式复杂、理解不便的特点,采用欧拉控制体方法,把“质量通量”整体作为一物理量,从而巧妙地推导了这两类连续方程,该过程物理意义明确、数学算法简单,有助于学生理解。 关键词:连续方程;柱坐标系;球坐标系 在大学《流体力学》教学中,连续方程是最基本的内容之一,在很多相关专业课程中得到广泛应用。相对而言,在直角坐标系中的连续方程形式简单,也易于理解,但在柱坐标系和球坐标系中,连续方程的形式却相对复杂,理解相对困难。目前,很多参考书[123] 对于后两类连续方程要么没有给出具体推导,要么推导过程较为复杂,使数理基础较薄弱的学生难以理解,在此,笔者结合教学中的实际经验,演示柱坐标系和球坐标系下一种物理意义明确、数学理解简单的连续方程的推导过程。 1 连续方程的一般算子形式 流体运动的连续方程,是表示流体运动和其质量分布的关系式。在拉格朗日方法中,某流体块在运动时其体积和形状尽管可发生变化,但它始终由这些流点构成,因此它的质量不变。 由此可见,连续方程实质上是质量守恒定律在“连续介质” (流体)中的应用。一般的拉格朗日方法考虑,某个别流体微团(质量体)在运动过程中,其随体密度的变化,必然与其体积变化趋势相反,如体积膨胀,它的密度减小,体积收缩,则密度增大。其算子形式的通用表达式[1] (1) 一般的欧拉方法考虑,对于某固定位置的空间单位体积元(控制体)来说,该体积元内单位时间的质量变化,与该体积元边界上的质量通量变化相联系,如质量往外流,它的密度减小,反之 则增大。其算子形式的通用表达式[1] (2) 两种方法的区别:拉格朗日方法多从物理量的定义出发,模型简单容易理解,但数学解析在实际应用中有些困难;欧拉方法则通过适当的数学建模后,能在数学上给出方便的解析,有利于从数学角度更好地理解概念。 在直角坐标系中,通过建立三维空间微元控制体(图略,很多教科书都详细给出,且易于理解),很容易得到(2)式在三维直角坐标系下连续性微分方程的一般表达式

数学物理方法(刘连寿第二版)第07章习题

第七章 习题答案 7.1-1将Helmholtz 方程0=+?u u λ在柱坐标系中分离变量。 解：0112 2 22 2=+??+??+???? ??????=+?u z u u u u u λ?ρρρρρλ 设)()()(),,(z Z R z u ?φρ?ρ=代入上面的方程有: 0d d d d d d d d 2 2 22 2=Φ+Φ+Φ+???? ??ΦZ R z Z R RZ R Z λ?ρρρρρ 两边同时除以Z R Φ，并移项得： λ?ρρρρρ--=ΦΦ+???? ??2 2 22 2d d 1d d 1d d d d 1 z Z Z R R 上式左边与z 无关，右边与?ρ,无关，令左右两边都等于μ-,即： 右边为：0)(d d d d 12 2 22 =-+? -=-- Z u z Z z Z z λμλ ① 而左边有：μ? ρρρρρ-=ΦΦ+???? ??2 2 2d d 1d d d d 1 R R 两边同时除以2ρ，并移项得: 2 2 2 2 d d 1d d d d m R R =ΦΦ- =+??? ? ?? ? μρρρρρ 0d d d d 12 2 2 2 2 2 =Φ+Φ? =ΦΦ- m m ? ? ② 和： 2 2 d d d d m R R =+??? ? ?? μρρ ρρρ 2 2 2 2 d d d d m R R R =+??? ? ??+ μρρρρρ 0d d 1d d 2 2 22 =??? ? ??-++R m R R ρ μρρρ ③ Helmholtz 方程在柱坐标系下可分解为①②③三个常微分方程。 7.1-2 将三维热传导方程 02 =?-??u a t u 在球坐标系中分离变量。

分离变量法：球坐标系下的分离变量法

4.5 分离变量法的原理与应用 1.分离变量法的原理和步骤 2.直角坐标系中的分离变量法 3.圆柱坐标系中的分离变量法 4.球坐标系中的分离变量法

4. 球坐标系中的分离变量法2 222222111()(sin )0sin sin R R R R R R φφφθθθθθ? ?????++=?????(,,)()Θ()Φ() R R φθ?θ?=R 22 22 sin d d sin d d Θ1d Φ()(sin )0d d Θd d Φd R R R θθθθθ? ++=R R 球坐标系中的拉普拉斯方程： （1）变量分离： 可得：

2 2 2 1d d m ? Φ=-Φ21d d ()(1)d d R n n R R =+R R 2 21d d (sin )(1)0 sin d d sin m n n θθθθθ Θ++-=Θ该方程只讨论电位与方位角无关的情况 该方程称为欧拉方程 该方程的解有两种情况 （2）分离出常微分方程： 22 22sin d d sin d d Θ1d Φ()(sin )0d d Θd d Φd R R R θθθθθ? ++=R R 根据：

0,φ? ?=?2 m =Φ()A ?=2 0n =1 00()R A B R -=+R (1)时，2 0n >(1) ()n n n n R A R B R -+=+R (2)时，当电位与方位角无关时即：2 2 2 1d ΦΦd m ? =-■ 的解■ 的解21d d ()(1)d d R n n R R =+R R 2 0n <的情况不存在。 (3) （3）常微分方程的求解：