第3章 微分方程
自动控制原理第三章习题解答
![自动控制原理第三章习题解答](https://img.taocdn.com/s3/m/e54d407fa417866fb84a8e6f.png)
tp =
1−ξ 2
= e −π 0.6 /
1−0.62
= e −π 0.6 /
1−0.62
= 9 .5 %
π
1 − ξ ωn
2
=
π
1.6
= 1.96( s )
ts =
3-5
3.5
ξω n
=
3.5 = 2.92( s ) 1.2
设单位反馈系统的开环传递函数为
G ( s) =
0.4 s + 1 s ( s + 0.6)
s5 s4 s3 s2 s1 s0
1 12 35 3 20 25 16 80 3 3 5 25 10 25
有一对虚根,系统不稳定 3-13 已知单位反馈系统的开环传递函数
G ( s) =
K (0.5s + 1) s ( s + 1)(0.5s 2 + s + 1)
试确定系统稳定时的 K 值范围。 解:系统特征方程为
ε 0 ,试问 k1 应满足什么条件?
见习题 3-20 解答 3-2 设系统的微分方程式如下: (1)
&(t ) = 2r (t ) 0.2c
&&(t ) + 0.24c &(t ) + c(t ) = r (t ) (2) 0.04c
试求系统的单位脉冲响应 k(t)和单位阶跃响应 h(t)。已知全部初始条件为零。 解: (1) 因为 0.2 sC ( s ) = 2 R ( s ) 单位脉冲响应: C ( s ) = 10 / s 单位阶跃响应 h(t)
试求系统的超调量σ%、峰值时间tp 和调节时间ts。 解: h(t ) = 1 −
第3章3.2动态电路的方程及其解
![第3章3.2动态电路的方程及其解](https://img.taocdn.com/s3/m/11cea7737fd5360cba1adb74.png)
§3.2
动态电路的方程及其解
■
第 1页
§3.2 动态电路的方程及其解
描述动态电路的方程是微分方程。 描述动态电路的方程是微分方程。用一阶微分 微分方程 方程描述的电路常称为一阶电路 一般而言, 一阶电路。 方程描述的电路常称为一阶电路。一般而言,如果 电路中含有n个独立的动态元件 个独立的动态元件, 电路中含有 个独立的动态元件,则描述它的将是 n阶微分方程,该电路可称为 阶电路。 阶微分方程, 阶电路。 阶微分方程 该电路可称为n阶电路
• 动态电路方程的建立 • 微分方程的经典解法
▲
■
第 2页
一、动态电路方程的建立
1、依据:元件VAR,KCL和KVL列写方程; 依据:元件VAR,KCL和KVL列写方程; VAR 列写方程 uR 一阶电路举例: 2、一阶电路举例: R i S RC电路 t=0时开关 电路, 时开关S 例1:图RC电路,t=0时开关S闭 uC uS 讨论t>0时的电容电压u 。 t>0时的电容电压 合,讨论t>0时的电容电压 C(t)。 C t>0时 根据KVL KVL方程列出回 t>0时,根据KVL方程列出回 RC串联电路 uR + uC – uS = 0 路电压方程为 d uC d uC , uR = R i = RC 根据元件的VAR VAR, 根据元件的VAR,有 i = C 代入上式, 代入上式,整 理得
− 1 t RC +U S
uC (t) = (U0 −US ) e
,t ≥ 0
▲ ■ 第 11 页
3、结果分析
固有响应和强迫响应 暂态响应和稳态响应
− 1 t RC +U
uC (t) = (U0 −US ) e
第三章常微分方程的差分方法15
![第三章常微分方程的差分方法15](https://img.taocdn.com/s3/m/62f0346abdd126fff705cc1755270722192e5999.png)
1.教学内容:
Euler方法:Euler公式,单步显式公式极其局部截断误 差;后退Euler公式,单步隐式公式极其局部截断误差;梯 形公式,预测校正公式与改进Euler公式。
2.重点难点:
Euler公式,预测校正公式与改进Euler公式
3.教学目标:
了解欧拉方法的几何意义、对给出的初值问题,能利 用Euler公式,改进Euler公式进行数值求解
科学技术当中常常需要求解常微分方程的定解问题。这类
问题的最简单的形式,是本章着重要考察的一阶方程的初值 问题:
y ' f x, y
y
x0
y0
(1) (2)
本章中我们假定右函数适当光滑以保证初值问题解的存
在唯一。虽然求解常微分方程有各种各样的解析方法,但求 解从实际问题中归结出来的微分方程要靠数值解法。
(其解析解为) y 2x 1
解:设步长 h=0.1,由改进的欧拉格式(10)有:
y
p
yn
h( yn
2xn ) yn
yc
yn
h( y p
2 xn1 ) yp
yn
1
1 2
(yp
yc )
n=0时
yp
y(xn ))
替代方程
y' (xn1) f (xn1, y(xn1))
中的导数项 y'xn1 再离散化,即可导出下列格式
yn1 yn hf xn1, yn1
(5)
该格式右端含有未知的 yn1 它实际上是个关于 yn1
的函数方程。故称该格式为隐式欧拉格式。
由于向前差商和向后差商具有同等精度,故隐式欧拉 格式也是一阶方法,精度与欧拉格式相当。但计算远 比显式格式困难得多。
第3章 微分方程模型
![第3章 微分方程模型](https://img.taocdn.com/s3/m/2ce333fbfab069dc5022018b.png)
第三章 微分方程建模在许多实际问题的研究中,要直接导出变量之间的函数关系较为困难,但要导出包含未知函数的导数或微分的关系式却较为容易,此时即可用建立微分方程模型的方法来研究实际问题。
例如,根据自由落体运动的重力加速度g 为常数及初始条件即可得出自由落体运动的公式、根据单摆的受力分析及牛顿第二定理即可得到单摆运动满足的方程等等就是典型的实例。
本章除了介绍一些来自经典力学的物理及一些几何方面的微分方程问题以外,也介绍了一些稍有不同的微分方程应用题。
这些模型研究的主要是来自于非物理领域的实际问题,对这些问题,我们将分析其特征,根据具体情况进行类比,提出假设条件并建立微分方程模型加以研究。
提出的假设条件不同,将会导出不同的微分方程。
最后还要将求解的结果与实际现象进行对比,如果差异较大还应反复修改假设建立新的模型。
因此,在这类模型中,微分方程被当成了研究问题的工具。
事实上,在连续变量问题的研究中,微分方程或微分方程组还是十分常用的数学工具之一。
§3.1 几个简单实例例3.1 (理想单摆运动的周期)本例的目的是建立理想单摆运动满足的微分方程,由该微分方程即可得出理想单摆运动的周期公式。
(图3-1)从图3-1中不难看出,小球所受的合力为 sin mg ,根据牛顿第二定律可得:θθsin mg ml -= 从而得出两阶微分方程:sin 0(0)0,(0)g l θθθθθ⎧+=⎪⎨⎪==⎩ (3.1) 这就是理想单摆运动满足的微分方程。
(3.1)是一个两阶非线性常微分方程,不容易求解。
根据微积分知识,当θ很小时,有sin θ≈θ,此时,为简单起见,我们可考察(3.1)的近似线性方程:⎪⎩⎪⎨⎧===+∙∙∙0)0(,0)0(0ϑϑϑϑϑl g (3.2)(3.2)的特征方程为02=+lg λ 对应的特征根为i lg =λ,(其中i 为虚单位),故(3.2)中的微分方程的通解为: t c t c t ωωϑcos sin )(21+=,其中lg =ω 代入初始条件,即可求得满足初始条件的微分方程问题(3.2)的解θ(t )= θ0cos ωt注意到当4T t =时,θ(t ) = 0,即可得出 24πω==T l g t 故有 l g T π2=这就是中学物理中理想单摆运动周期的近似公式。
数学建模第三章微分方程模型
![数学建模第三章微分方程模型](https://img.taocdn.com/s3/m/3ec2b00e6294dd88d0d26bf3.png)
3-7 香烟过滤嘴的作用机理(2)
精选ppt课件
51
3-7 香烟过滤嘴的作用机理(3)
精选ppt课件
52
3-7 香烟过滤嘴的作用机理(4)
精选ppt课件
53
3-7 香烟过滤嘴的作用机理(5)
精选ppt课件
54
3-7 香烟过滤嘴的作用机理(6)
精选ppt课件
55
3-7 香烟过滤嘴的作用机理(7)
精选ppt课件
39
3-6 疾病传播的机理分析模型(2)
精选ppt课件
40
3-6 疾病传播的机理分析模型(3)
精选ppt课件
41
3-6 疾病传播的机理分析模型(4)
精选ppt课件
42
3-6 疾病传播的机理分析模型(5)
精选ppt课件
43
3-6 疾病传播的机理分析模型(6)
精选ppt课件
44
3-6 疾病传播的机理分析模型(7)
精选ppt课件
45
3-6 疾病传播的机理分析模型(8)
精选ppt课件
46
3-6 疾病传播的机理分析模型(9)
精选ppt课件
47
3-6 疾病传播的机理分析模型(10)
精选ppt课件
48
3-6 疾病传播的机理分析模型(11)
精选ppt课件
49
3-7 香烟过滤嘴的作用机理(1)
精选ppt课件
50
69
3-10 赤道上空通讯卫星颗数的确定(1)
精选ppt课件
70
3-10 赤道上空通讯卫星颗数的确定(2)
精选ppt课件
71
3-10 赤道上空通讯卫星颗数的确定(3)
精选ppt课件
自动控制原理第三章课后习题 答案()
![自动控制原理第三章课后习题 答案()](https://img.taocdn.com/s3/m/b2009abdd15abe23482f4de0.png)
3-1 设系统的微分方程式如下:(1) )(2)(2.0t r t c= (2) )()()(24.0)(04.0t r t c t c t c=++ 试求系统闭环传递函数Φ(s),以及系统的单位脉冲响应g(t)和单位阶跃响应c(t)。
已知全部初始条件为零。
解:(1) 因为)(2)(2.0s R s sC =闭环传递函数ss R s C s 10)()()(==Φ 单位脉冲响应:s s C /10)(= 010)(≥=t t g单位阶跃响应c(t) 2/10)(s s C = 010)(≥=t t t c(2))()()124.004.0(2s R s C s s =++ 124.004.0)()(2++=s s s R s C `闭环传递函数124.004.01)()()(2++==s s s R s C s φ单位脉冲响应:124.004.01)(2++=s s s C t e t g t 4sin 325)(3-= 单位阶跃响应h(t) 16)3(61]16)3[(25)(22+++-=++=s s s s s s C t e t e t c t t 4sin 434cos 1)(33----=3-2 温度计的传递函数为11+Ts ,用其测量容器内的水温,1min 才能显示出该温度的98%的数值。
若加热容器使水温按10ºC/min 的速度匀速上升,问温度计的稳态指示误差有多大解法一 依题意,温度计闭环传递函数11)(+=ΦTs s 由一阶系统阶跃响应特性可知:o o T c 98)4(=,因此有 min 14=T ,得出 min 25.0=T 。
视温度计为单位反馈系统,则开环传递函数为Tss s s G 1)(1)()(=Φ-Φ=⎩⎨⎧==11v TK !用静态误差系数法,当t t r ⋅=10)( 时,C T Ke ss ︒===5.21010。
解法二 依题意,系统误差定义为 )()()(t c t r t e -=,应有 1111)()(1)()()(+=+-=-==ΦTs TsTs s R s C s R s E s e C T sTs Ts ss R s s e s e s ss ︒==⋅+=Φ=→→5.210101lim )()(lim 203-3 已知二阶系统的单位阶跃响应为)1.536.1sin(5.1210)(2.1o tt et c +-=-试求系统的超调量σ%、峰值时间tp 和调节时间ts 。
《偏微分方程教程》第三章 特征理论与方程的分类
![《偏微分方程教程》第三章 特征理论与方程的分类](https://img.taocdn.com/s3/m/3da615c584868762cbaed529.png)
对二阶线性常系数方程作了比较详细的讨论, 但对变系数方程(2.1)
同样是成立的. 这里要特别指出的是, 对变系数方程来说, 它的类
型与点的位置有关, 即可能在区域的某一部分点为这种类型而在另
一部分点上为另一种类型.
18
例如特里谷来(Tricomi)方程
u yy yuxx 0
(2.13)
就是如此, 其判别式 y ,对于 y 0 它是双曲型的; 对于
2 x
x x
2 x
2xy x y y x
2 x y
2 y
y y J 3 0
2 y
知 A(), B(), C() 不同时为零.
利用判别式的符号在可逆自变量变换下的不变性这一性质, 我 们来对方程(2.1)进行分类.
9
定义3.1 设 R2是一个区域,(x0 y0 ) (i) 若 (x0 y0 ) 0 ,则称方程(2.1)在点 (x0 y0 ) 处为双曲型 偏微分方程, 若在 内的每一点处, 方程(2.1)都是双曲型的,则 称(2.1)在 内为双曲型偏微分方程;
设在点 P(x0 y0 )的邻域内, 这时(2.1)的特征方程可写为
dy b dy b (2.2)
dx
a dx
a
其中 b2 ac 通常称为方程(2.1)的判别式,作自变量变换
4
(x y) (x y)
(2.3)
则方程(2.1)变为如下形式:
A 2u 2B 2u C 2u F (2.4)
2
2
在自变量变换(2.3)下,方程(2.1)的判别式 与(2.4)的判别式
B2 AC 之间有如下关系:
第三章 微分方程方法
![第三章 微分方程方法](https://img.taocdn.com/s3/m/953a3834f111f18583d05a83.png)
第三章 微分方程方法3.1微分方程的一般理论微分方程是研究函数变化规律的有力工具,有着广泛的实际应用。
针对所研究的对象建立微分方程模型是解决问题的第一步,实际中只有求出微分方程的解才能对所研究的问题进行解释说明。
一般说来,求微分方程的解析解是困难的,大多数的微分方程需要用数值方法来求解,因此首先需要研究微分方程的解的存在惟一性和稳定性问题。
3.1.1 微分方程的一般形式一阶微分方程⎪⎩⎪⎨⎧==00)(),(x t x x t f dtdx(3.1) 其中),(x t f 是t 和x 的已知函数,00)(x t x =为初始条件,又称定解条件。
一阶微分方程组⎪⎩⎪⎨⎧====),2,1( )(),2,1( ),,,,()0(021n i x t x n i n x x t f dtdx i i i i(3.2) 又称为一阶正规方程组。
如果引入向量T n x x x x ),,,(21 =,Tn x x x x ),,,()0()0(2)0(10 =,Tn f f f f ),,,(21 =,Tn dt dx dt dx dt dx dt dx ⎪⎭⎫ ⎝⎛=,,,21 。
则方程组(3.2)可以写为简单的形式⎪⎩⎪⎨⎧==00)(),(x t x x t f dtdx(3.3) 即与方程(3.1)的形式相同,当1=n 时为方程(3.1)。
对于任一高阶的微分方程⎪⎪⎭⎫⎝⎛=--11,,,,n n n n dt x d dt dx x t f dt x d , 如果记),,2,1,0(n i y dtxd i i i ==,则方程为),,,;(1101--=n n y y y t f dt dy 即可化为一阶方程组的形式。
因此,下面主要对正规方程组(3.3)进行讨论。
3.1.2微分方程解的存在惟一性正规方程组(3.3)的解在什么条件下存在,且惟一呢?有下面的定理。
定理3.1(Cauchy-Peano )如果函数),(x t f 在区域b x x a t t R ≤-≤-00;:上连续,则方程组(3.3)在h t t ≤-0上有解)(t x φ=满足初值条件)(00t x φ=,此处),(max ,,min ),(x t f M M b a h R x t ∈=⎪⎭⎫⎝⎛=。
第三章二阶线性偏微分方程的分类化简
![第三章二阶线性偏微分方程的分类化简](https://img.taocdn.com/s3/m/c5031dd8c1c708a1284a4471.png)
n e m
f o t
t a M
a m he
, s c ti
T I H
D
t r a ep
n e m
f o t
t a M
a m he
, s c ti
T I H
D
t r a ep
n e m
f o t
t a M
a m he
, s c ti
T I H
D
t r a ep
n e m
f o t
t a M
a m he
au xx bu xy cu yy du x eu y fu g
D
t r a ep
t a b 4ac 0 则称方程(1)是双曲型方程 m 在区域D内,如果 e h b 4ac 0 t如果 a 则称方程(1)是抛物型方程 M f o 如果 b 4ac 0 则称方程(1)是椭圆型方程 t n e m二)二阶线性偏微分方程的化简
( x, y ).
2 2 C a x b x y c y 0
经过证明可以上述一阶偏微分方程的解等价于常微分方程: ady 2 bdxdy cdx2 0 ( a( dy )2 b dy c 0 )
dx dx
(3)
D
t r a ep
xx xy yy
2 x x y 2 y x x x y y x y y 2 x x y 2 y xx xy yy x y xx xy yy x y
x
Eu y Fu G
(2)
F f Gg
假设abc不等于零,我们希望选取一个变量替换使得方程(2)中A和B都等 2 2 ( x, y ), 使得 于零.选择变换, A a x b x y c y 0
第三章 常微分方程
![第三章 常微分方程](https://img.taocdn.com/s3/m/053c730816fc700abb68fc14.png)
(3) 日治愈率为 , 病人治愈后仍可被感染
推知:
di N Ni (1 i ) Ni dt
di i (1 i ) i dt s (t ) i (t ) 1
目录 上页 下页 返回 结束
SIS模型的解
当 时,有
1 i (t ) [ ( )e ( )t ]1 i0
目录
上页
下页
返回
结束
人口分类(在此用绝对人数)
1、易感类S,该类成员没有感染SARS,也没有免疫
能力,处于自由环境,可以被传染上SARS
2、潜伏期类E,该类成员已经感染了SARS病毒,但 处于潜伏期,还不是SARS患者,不能把病毒传 染给S 类成员
3、患病未被发现类Iu,该类成员已经成员真正的 SARS患者,可以把病毒传染给S 类成员
目录 上页 下页 返回 结束
二、模型的建立—机理分析法
1. SI模型
假设条件
() 1 人群分为易感者s(t )和已感者i(t )两类
(2) 每个病人每天有效接触的平均人数为
推知:
每个病人每天可使 s(t )个健康者变为病人,病人数
为Ni(t ),所以每天有 Ns(t )i(t )个健康者变为病人.
目录 上页 下页 返回 结束
关于SARS的假设
1、所研究地区的人口总量一定,不考虑人口迁入 迁出、自然出生和自然死亡 2、结合SARS的病理特性,假设其传播方式为接 触性传播,不与患病者接触就不会被感染 3、假设病人被感染后需先进入潜伏期,在潜伏期 内不具备传染性
目录 上页 下页 返回 结束
4、假设易感者与患病者在人群中混合充分均匀, 易感者感染的机会与接触患病者的机会成正比, 并且传染率为常数 5、假设SARS患者被发现后就立即被隔离,被隔离 者不具备传染性,SARS患者只在被发现前可以 传染他人 6、假设SARS康复者不会被再次感染,且无传染性
数理方程第三章(1)
![数理方程第三章(1)](https://img.taocdn.com/s3/m/acc48a1b6bd97f192279e9cb.png)
为正、为零或者为负而确定的。 或者为负而确定的。 如果方程在一个区域内的每点都是双曲型、抛物 如果方程在一个区域内的每点都是双曲型、 型或椭圆型的, 型或椭圆型的,那么就称方程在这个区域内是双 曲型、抛物型或椭圆型。 曲型、抛物型或椭圆型。
双曲型方程 注2:行波法适用于双曲型方程。 :行波法适用于双曲型方程。
x = x1 + at
B
x1
x = x2 − at x2 x
B = {( x, t ) | x1 + at ≤ x ≤ x2 − at , t ≥ 0}
决定区域。 称三角区域 B 为区间 [ x1 , x2 ] 的决定区域
进一步的分析其物理意义表明, 进一步的分析其物理意义表明, 在 xot 平面上
x1 − at ≤ x ≤ x2 + at
称区域
(t > 0).
A = {( x, t ) | x1 − at ≤ x ≤ x2 + at , t > 0}
t
为区间 [ x1 , x2 ] 的影响区域。 影响区域
x = x1 − at
x = x2 + at
A
x1
x2
x
(3) 决定区域
t
考虑区间 [ x1 , x2 ],
∂ 2u ∂ 2u ∂ 2u = 2 +2 + 2, ∂ξ ∂ξ∂η ∂η
(3.1.4)
同理可得
∂u2 ∂ 2u ∂ 2u ∂ 2u 2 = a ( 2 −2 + 2) 2 ∂t ∂ξ ∂ξ∂η ∂η
(3.1.5)
将 (3.1.4), (3.1.5) 代入到 (3.1.1), 可以得到
∂ 2u = 0. ∂ ξ∂ η
《传热学》第3章-非稳态导热
![《传热学》第3章-非稳态导热](https://img.taocdn.com/s3/m/0849957c336c1eb91a375d7d.png)
特殊多维非稳态导热的简易求解方法
在第一类边界条件(初始温度均匀)或第三类边界条件(表面 传热系数h为常数)下的二维或三维的非稳态导热问题,在数学 上已经证明,它们的无量纲过余温度的解等于构成这些物体的 两个或三个物体在同样边界条件下一维非稳态导热问题解的连 乘。
特殊多维非稳态导热的简易求解方法
对于无限长方柱 θ (x, y,τ ) = θ (x,τ ) ⋅ θ (y,τ )
该问题的解可以由3块相应的无限大平板的 解得出。最低温度发生在钢锭的中心,即3 筷无限大平板中心截面的交点上,最高温度 发生在钢锭的顶角,即3块大平板表面的公 共点上。
4
例题3 θ
m/B则θi x0钢==锭hλδ(1θ中=m心3/ 4θ温840×0度).05x.2⋅5(θ=
2.14
m/θ 0
)
y
⋅ (θ
无限大平板的非稳态导热
当Fo ≥ 0.2时,可取
θ (x,τ )
θ0
=
β1
2 sin β1 + sin β1 cos β1
cos
β
1
x δ
e − β12 ⋅Fo
只与Bi、x/δ有关, 与时间无关
lnθ
=
−mτ
+ lnθ 0
β1
2sin β1 + sinτ β1 cos β1
cos
= 0.36
短圆柱的中心温度为
查图3-6得 θ
再讨论直径为
m2R/θ=600=0m0m.8的无θ限m长/ θ圆0柱=:0.13
×
0.8
=
0.104
Bi = hR = 232 × 0.3 = 1.72 λ 40.5
tm = 0.104θ0 + t∞ 查附=2图0.11得04θ×m(3/θ00−=103.0103) +1300
弹性与塑性力学基础-第3章平衡微分方程及应变协调方程
![弹性与塑性力学基础-第3章平衡微分方程及应变协调方程](https://img.taocdn.com/s3/m/87f419480b4e767f5bcfce12.png)
§3-3 二维极坐标系下的平衡微分方程
3.3.2 二维极坐标系下的平衡微分方程的建立 ➢ 微分体切向平衡方程
ddrdrr
r d(rd)rd
rrdr r ddrd2 rdrd2 Krddr0
➢ 用r代替r ,简化以后,除以rddr,再略去微量,得
1 r rr2rrK0
2 x y 2
2 y x2
2 xy xy
2 y z2
2 z y 2
2 yz yz
➢ 当六个应变分量
2 z x2
2 x z2
2
2 x yz
x
2 xz
zx yz x
xz y
xy z
(3-7)
满足以上应变协 调方程(3-7)时,
2 2 y zx
y
弹性与塑性力学基础
第三章
平衡微分方程及应变协调方程
2020/10/13
弹性与塑性
力 学 基 础 第三章 平衡微分方程及应变协调方程
§3-1 平衡微分方程的概念
3.1.1 平衡微分方程的概念 3.1.2 平衡微分方程的建立
§3-2 二维直角坐标系下的平衡微分方程
3.2.1 平面应力状态 3.2.2 平面应变状态
➢ 通过中心C并平行于z轴的直线为矩轴,力矩平衡方程 MC=0:
xy xxyd x d y1d 2 xxd y y1d 2x
yx y yxd y d x 1d 2 yyd x x 1d 2 y0
将上式除以dxdy,得到ຫໍສະໝຸດ y1 2xyx
dx
=
yx
1 2
yx
x
dy
2020/10/略13 去微量,(亦即dx、dy都趋于零时),得出
第3章-流体力学连续性方程微分形式
![第3章-流体力学连续性方程微分形式](https://img.taocdn.com/s3/m/7528e4ed76eeaeaad0f33016.png)
欧拉运动微分方程组各式分别乘以dx,dy,dz(流场任意相邻两点间距ds 的坐标分量),然而相加得:
( Xdx
Ydy
Zdz)
1
( px
dx
p y
dy
p z
dz)
dux dt
dx
duy d;
<II>
<III>
一、在势流条件下的积分
考虑条件 1、恒定流
当为恒定流时
t
0
(ux
x
)
(uy
y
)
(uz
z
)
0
适用范围:理想、实际、可压缩、不可压缩的恒定流。
(2)不可压缩流体的连续性微分方程
当为不可压缩流时 Const
u x x
u y y
u z z
0
物理意义:不可压缩流体单位时间内流入单位空间的流体体积(质量) ,
dt
p'xx 'xz 'xy
x
第三节 流体动力学基本方程式
考虑条件:
13
1)
不可压缩流体的连续性微分方程:uxx
uy y
uz z
0
2)切应力与主应力的关系表达式
• 不可压缩粘性流体运动微分方程:纳维埃-斯托克斯方程(Navier-
Stokes,N-S)方程:
X
1
p x
1
第三章 流体动力学基础
第三节 流体动力学基本方程式
一、连续性微分方程 二、理想流体运动微分方程 三、粘性流体的运动微分方程
第3章常微分方程
![第3章常微分方程](https://img.taocdn.com/s3/m/4523ec221fd9ad51f01dc281e53a580216fc5008.png)
2K2
2K3
K4
)
K1
f (xn , g(xn ,
yn , yn ,
zn ) zn )
K2
f (xn g ( xn
h, 2 h, 2
yn yn
h 2 h 2
K (1) 1
,
zn
K (1) 1
,
zn
h 2 h 2
K1(
2)
)
K1(2) )
K2
F(x
h 2
,Y
h 2
K1)
K3
f g
y(
xn1
))]
h3 12
f
''( )
所以,有格式为:
yn1
yn
h[ 2
f
(xn ,
yn )
f
( xn1 ,
yn 1 ) ]
--------梯形格式
类似,可以算出其误差估计式:
en,要用迭代法求解
将梯形法和Euler法相结合,可得到改进的Euler法:
yn1
3
ha2
xn1 (x xn )( x xn1) dx 1 h
xn1 (xn2 xn )( xn2 xn1)
3
Tn1
xn1 xn1
y(4) (
(3)!
)
(x
xn
)( x
xn1)( x
xn2
)dx
1 3
h4
y(4)
( )
例:建立p=2,q=2的隐格式
p=2,积分区间为 xn1 y'(x)dx xn2
➢ 基于数值积分的构造法
将 y f ( x, y) 在[xn p , xn1] 上积分,得到
第3章 常微分方程的差分方法
![第3章 常微分方程的差分方法](https://img.taocdn.com/s3/m/f51bb46be518964bce847c25.png)
a = x0 < x1 < … < xi < … < xn = b
上的近似值y0, y1, …, yn。
两相邻节点间的距离
hi = xi+1 - xi (i=0,1,2,…,n-1)
称为步长。当
hi h (常值)时称为等步长,有
x x i, h ( i 1 , 2 ,... n ) i 0
Pn 1
Pn
Pn
y y(x)
Pn 1
o
x 0 x1 x 2
x i x i 1 x n 1 x n
图7-2
x
与14Leabharlann 得到 y2 作为y(x2 ) 的近值;……如此继续,直到Pn 点。这样,得出一条折线P0 P1P2…Pi…Pn 近似代
替积分曲线P'0 P'1P'2…P'i…P'n 。当步数越多时,
步格式。其中(2.5)是一个数值微分公式。故用其他数值微 分公式也可导出略异于( 2.2)的其他形式的算式来。例如, 用向后差商表示的两点数值微分公式
1 h y ( x ) y ( x ) y ( x ) y ( ) i 1 i 1 i i h 2 ( i 0 , 1 , ,n 1 )
y(xi+1 )和 y(xi )用其近似值 yi+1 和 yi 代入,则得
y y hf ( x , y )( i 0 , 1 , 2 ,..., n 1 )
i 1 i i i
16
此即(2.2)(欧拉格式)。 显然,欧拉格式具有递推性,在计算yi+1时只要用到前一
步所得结果 yi 一个信息就够了,因此是一种单步格式或称一
第3章 单自由度系统的受迫振动
![第3章 单自由度系统的受迫振动](https://img.taocdn.com/s3/m/ca88fad6f121dd36a22d8210.png)
于 90°。在振动实验中,常以此作为判断振动系统是否处于共振状态的一种标志。 (3) 高频区。当λ>>1, ϕ=180°。表明当激振力频率远远高于固有频率时,受迫振动的相位差接
近与 180°。这说明受迫振动的位移与激振力是反相位的。 应当指出,对于λ=0,当λ<1 时,λ =0;λ>1 时,ϕ=180°;λ=1 时,ϕ角从 0 跳到 180°。 对于不同的阻尼值,相位差ϕ角在 0 到180° 之间变化。 例 3-1 质量为 M 的电机安装在弹性基础上。由于转子不均衡,产生偏心,偏心距为 e,偏心质
pn
(3-10)
绘出对应不同的阻尼比ζ,相位差ϕ随λ变化的曲线族如图 3-2 中的右上角所示,即相频特性曲线。
(1) 低频区。当λ<<1时,ϕ≈0,表明当激振力频率很低或ω<< pn 时,相位差ϕ接近于零,即受
迫振动的位移与激振力几乎同相位。
(2) 共振区。当λ=1时,ϕ=90°。表明当激振力频率等于振动系统的固有频率时,相位差为 90°。
根据达朗贝尔原理,有
− cx& + Mg − k(x + δ st ) − M&x& − meω 2 sin ωt = 0
∴ M&x& + cx& + kx = −meω 2 sin ωt
令
p
2 n
=
k M
,2n
=
c M
,则上式可写成
&x& +
2nx&
+
pn2 x
第3章 振动系统的运动微分方程题解
![第3章 振动系统的运动微分方程题解](https://img.taocdn.com/s3/m/4ee32e45284ac850ac02420d.png)
资料范本本资料为word版本,可以直接编辑和打印,感谢您的下载第3章振动系统的运动微分方程题解地点:__________________时间:__________________说明:本资料适用于约定双方经过谈判,协商而共同承认,共同遵守的责任与义务,仅供参考,文档可直接下载或修改,不需要的部分可直接删除,使用时请详细阅读内容习题题3-1图3-1 复摆重P,对质心的回转半径为,质心距转动轴的距离为a,复摆由水平位置无初速地释放,列写复摆的运动微分方程。
解:系统具有一个自由度,选复摆转角为广义坐标,原点及正方向如如题4-1图所示。
复摆在任意位置下,根据刚体绕定轴转动微分方程其中得到复摆运动微分方程为或3-2均质半圆柱体,质心为C,与圆心O1的距离为e,柱体半径为R,质量为m,对质心的回转半径为,在固定平面上作无滑动滚动,如题3-2图所示,列写该系统的运动微分方程。
题3-2图解:系统具有一个自由度,选为广义坐标。
半圆柱体在任意位置的动能为:用瞬心法求:故系统具有理想约束,重力的元功为应用动能定理的微分形式等式两边同除,,等式两边同除故微分方程为①若为小摆动,,并略去二阶以上微量,上述非线性微分方程可线性化,系统微摆动的微分方程为要点及讨论(1)本题也可以用平面运动微分方程求解。
系统的受力图与运动分析图如图(b)所示。
列写微分方程上述方程包含,,,,五个未知量,必须补充运动学关系才能求解。
建立质心坐标与广义坐标之间的关系,所以运动学方程式⑤⑥与方程②③④联立,消去未知约束力,,就可以得到与式①相同的系统运动微分方程。
因为在理想约束的情况下,未知约束力在动能定理的表达式中并不出现,所以用动能定理解决已知力求运动的问题更简便、直接。
(2)本题也可用机械能守恒定律求解。
系统的动能选半圆柱体中心O1所在平面为零势面,系统的势能由两边对时间求导数,即可得到与式①相同的运动微分方程。
3-3 均质杆AB,长l,质量为m,沿光滑墙面滑下,如题3-3图所示。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
流出-流入 积累率
( u ) ( u ) ( u ) x y z dxdydz x y z
dxdydz t
( u x ) ( u y ) ( uz ) 0 x y z t
P38
柱坐标中的表达式:
1 1 rur u uz 0 t r r r z
式中: r 径向坐标;z 轴向坐标;θ方位角
(2-14)
d u r dt
表示线速度
直角坐标系与柱坐标系的关系
P38
图2-2b 示出球坐标系与直角坐标系的关系。
x 向质量净出率:
x 向增加量
( u x )dx dydz x
① 质量净出率
x 向质量净出率: y 向质量净出率:
z 向质量净出率:
三向六面净出率:
(kg/s)
( u x )dx dydz x
( u y )dy dxdz y
( uz )dz dxy z
本节内容:
1.能量方程( E.E. )的建立 2. 其他坐标系下的形式 3. E.E.的简化
4. E.E.的应用举例
1. 方程的建立 ( E.E. )
衡算根据:热力学第一定律,即某过程中体系从环境中所吸
收的热量减去对体系所作功之差,等于该体系在过程前后的能量 变化,其数学表达式为:
在传热、传质过程中往往伴随有流体流动,因
此均需要用到C.E.。
2.连续性方程的另一表达式
( u ) 0 t
引入随体导数的概念,上式可写为:
( u ) 0 u u t t
或
D u 0 Dt
① 稳态时
(以直角坐标为例)
( u ) 0 C . E . t
0 或 t
( u) 0
② 不可压缩(均质)流体 ( const )
u 0
(2-13)
5. 方程衡算的应用举例
【例】导出流体在圆管中作轴对称二维(径向、轴向) 流动时的连续性方程。
⑤ 恒密度,恒物性流体。
能量衡算图
1)流动带来的内能变化 2)导热带来的能量变化 3)流动功带来的变化 4)能量积累
U - 1kg质量所含的内能
( Udxdydz ) t
1)能量净出率 2)吸热率
( U u)dxdydz
(J/s)
( q )dxdydz
( pu)dxdydz
(2-7b)
讨论上述连续性方程中两项的物理意义。
D u 0 Dt
讨论方程的两项物理意义 第一项:
D u Dt t
其中
密度的随体导数
t
为密度的局部变化率。
表明在空间某一固定点处,密度随时间的变化速率,
反映了密度场的不稳定性。
D u Dt t
相加,整理得:
(2-6a)
矢量式
( u ) 0 t
(2-6b)
微分质量衡算方程(连续性方程)
在连续性方程的推导过程中没有做任何假定,
所以方程是通用方程,可用于:
牛顿流体 或 非牛顿型流体;
理想流体 或 真实流体;
可压缩流体 或 不可压缩流体。
连续性方程是描述速度分布的基本方程。
其中
i j k x y z
密度梯度
表明在不同位置上密度的变化率,反映了密度场的不均匀性。
u u x uy uz x y z
密度迁移变化率
表明当流体运动时,由于密度场的不均匀性造成的密度变化率
又称之为密度的对流变化率。速度越快,场越不均匀,对流变化率就越大。
t
u
u
3.其他坐标系连续性方程的表达式
在某些实际场合,应用柱坐标或球坐标来表达连
续性方程或其他微分衡算方程要比直角坐标系方便。
例如在研究管内的流体流动时,在相同半径上的所
有各点有相同的速度及其它物理量,在此情况下采用 柱坐标表达连续性方程最为方便。
又如流动系统的范围面为球形或球形的一部分时,则
3.其他坐标系下的表达式 4、方程的简化 5、方程的应用举例
1.连续性方程的建立
衡算根据:质量守恒原理 研究对象:微元控制体(记为C.V.),见下图
衡算原则 :
流出C.V. 质量流率
-
流入C.V. 质量流率
+
C.V.内 质量积累率
=0
在流动的流体中,取 P 点,在P点旁分别截取 dx、dy、dz,构成微元体积。
x ( u x ) y ( u y ) z ( uz ) dxdydz
② 质量积累率 (kg/s)
t 时刻C.V.质量为:
dxdydz
经过 dt 时间后,C.V. 内的质量变化量为:
( dt )dxdydz t
dt 时间内,C.V. 质量平均积累率为:
包围控制体的6个控制面分别为上下、左右、前后面。
微元控制体 C.V.中的质量衡算
C.E.方程的建立
① 质量净出率 (kg/s) 在 x 处, t 时刻流入的质量通量为:
u x
kg m2 s
kg m2 s
在 x+dx 处, t 时刻流出的质量通量为:
( u x ( u x )dx) x
得:
i ni ri 0 t
组元微分质量衡算方程
i ni ri 0 t
式中:
ni ( ji i u)
对此项求导
ji i u u i D i i 合并为: u i 随体导数 Dt t
所以 组元微分质量衡算方程 还可以写成:
D i i u ji ri 0 Dt
(9-41)
D i i u ji ri 0 Dt
组元微分质量衡算方程的简化
1)不可压缩流体
u 0
D i ji ri 0 Dt
ni dxdydz
净出C.V. 组元A流率
+
C.V.内 组元A积累率
-
C.V.内 组元A生成率
=0
第一项
ni dxdydz
( i dxdydz ) i dxdydz t t
第二项 组元积累率: 第三项 组元生成率:
ri dxdydz
三项相加消去dxdydz
2)静止流体
i Dij 2 i ri t
当无化学变化时
i Dij 2 i t
费克第二定律(Fick’s second law)
三、能量微分衡算方程
在工程实践中,有一类很重要的问题,就是流体与壁
面、流体内部、固体内部的热量传递现象。
解析这类问题要利用微分能量衡算方程。
E q W
衡算原则:
能量净出率 + 能量积累率 = 吸热率
-
对外做功率
方程的建立
假定:
① 研究对象是微元控制体,所以动能,位能忽略不计,
只计内能:
△ U = q - △ (pv)
② 温度不是很高,不计辐射热,只计导热和对流传热。 ③ 对于微元控制体,不包含轴功,只计流动功。
④ 不涉及高粘度流体或是高速 (超音速) 下运动的流体, 因此忽略摩擦损耗功率。
研究对象: 微元控制体 衡算根据: 质量守恒 衡算原则:
净出C.V.
组元 流率
+
C.V.内 组元 积累率
-
C.V.内 组元 生成率
=0
绕微元控制体对组元衡算,对上述三项分别加以衡算。 用
ni 代表通过控制面的组分i 的质量通量。
ni i u ji
nix dx x
x 向组元质量通量净变化率
薄壳C.V.内质量累积率为:
ur 2rdz
2rdrdz t
三项相加 ,略去高阶无穷小 ,得到沿管流动的连续性方程
1 ( rur ) ( u z ) 0 t r r z
二、组元微分质量衡算方程(质量基准)
建立多组分流体的微分质量衡算方程。
体系的扩散传质过程。
改写成以温度表示的能量方程,具体步骤为
在左侧第一项中,根据焓与内能的关系,写成 Uu pu u(U p ) Hu CpTu 在左侧第二项中,对于恒密度流体
U CvT CpT
采用球坐标系的连续性方程描述最为适宜。
其他坐标系
柱坐标和球坐标系中的连续性方程和其它传递方
程的推导,原则上与直角坐标系相类似,可采用
直接推导获得。
也可以通过坐标系的对应关系由直角坐标间接
转换而得到相关方程,但由于运算繁杂,一般
不宜采用。
图2-2a 表示柱坐标系与直角坐标系的关系。 连续性方程在柱坐标系的表达式为:
传递微分衡算方程
一、质量微分衡算方程 C.E. 二、组元微分衡算方程 C.EA. 三、能量微分衡算方程 E.E. 四、动量微分衡算方程 N-S
它们是传递过程数学模型的主要内容; 是从理论上研究传递过程的基础。
一.微分质量衡算方程
(连续性方程 C.E. )
1.连续性方程的建立
2.连续性方程的另一表达式
P213
在后面的章节中主要考虑双组分(A,B) 考虑可能有化学反应的存在的情况。
衡算时传质包括两部分:
1.由于流体流动引起的 组分
2.由于浓度梯度引起的 组分 这两部分传质表示为: