21. 对于在r平面内的不可压缩流体的流动,r方向的速度分量为urAcos/r。试
确定速度的分量。
解:柱坐标系的连续性方程为
1(rur)1(u)(uz)0
rrrz对于不可压缩流体在r平面的二维流动,常数,uz0,11u(rur)0 rrruzz0,故有
即
urrrr将上式积分,可得
AcosAsin udf(r)
r2r2式中,f(r)为积分常数,在已知条件下,任意一个f(r)都能满足连续性方程。令
(rur)(rAcos2)Acos2
f(r)0,可得到u的最简单的表达式:
uAsin r2 2.对于下述各种运动情况,试采用适当坐标系的一般化连续性方程描述,并结合下述具体条件将一般化连续性方程加以简化,指出简化过程的依据。
(1)在矩形截面管道内,可压缩流体作稳态一维流动; (2)在平板壁面上不可压缩流体作稳态二维流动; (3)在平板壁面上可压缩流体作稳态二维流动;
(4)不可压缩流体在圆管中作轴对称的轴向稳态流动; (5)不可压缩流体作球心对称的径向稳态流动。
解: u0
(1) 在矩形截面管道内,可压缩流体作稳态一维流动
uxuyuzuxuyuz0 xyzxyz稳态:
0,一维流动:ux0, uy0 ∴ uz(uz)uz0, 即 0 zzz(2)在平板壁面上不可压缩流体作稳态二维流动
(ux)x(uy)y(uz)z0
稳态:
0,二维流动:uz0 (ux)x(uy)y0, 又const,从而
∴
uxuy0 xy(3)在平板壁面上可压缩流体作稳态二维流动 在此情况下,(2)中const
∴
(ux)x(uy)y0
(4)不可压缩流体在圆管中作轴对称的轴向稳态流动
11ruruuz0 rrrz稳态:
0,轴向流动:ur0,轴对称:0 uzu0, z0 (不可压缩const) zz∴
(5)不可压缩流体作球心对称的径向稳态流动
1112(r2ur)(usin)(u)0 rrrsinrsin0,不可压缩const 0,沿球心对称0,
稳态
∴
12d2 ,即 (ru)0(rur)0 r2rrdr3.某粘性流体的速度场为
u=5xyi3xyzj8xzk
已知流体的动力粘度0.144Pas,在点(2,4,-6)处的法向应力yy100N/m,试求该点处的压力和其它法向应力和剪应力。
22解: 由题设 ux5xy,uy3xyz,uz8xz
222 u10xy3xz16xz
uxx10xy,
uyy3xz,
uz16xz z因 yyp2uyyuy2uxuyuz() y3xyz23故 pyy2(uxxuyyuzz)
在点(2,4,-6)处,有
2 p(100)20.144(36)32 0.14423667N/mux2uxuyuz() 所以 xxp2x3xyz26720.144800.144236 366.6N/m2uz2uxuyuz() zzp2z3xyz 34.4N/m
2uxuyxyyx()
yx0.144[52234(6)]7.5N/m2
yzzy(uzuy) yz0.1443243.5N/m2
zxxz(uxuz) zx2 0.144(836)41.5N/m
4. 某不可压缩流体在一无限长的正方形截面的水平管道中作稳态层流流动,此正方形截面的边界分别为xa和ya,有人推荐使用下式描述管道中的速度分布
uzx2y2[1()][1()] 4zaaap2试问上述速度分布是否正确,即能否满足相关的微分方程和边界条件。
解: 在壁面处,即xa和ya时,uz0,故满足壁面不滑脱条件;在管道中心,xy0时,可得
a2puzumax (1)
4z将所给速度分布式代入不可压缩流体连续性方程(2-20),因ux
uy0可得
uz0 z将不可压缩流体的运动方程(2-45c)化简,可得
2uz2uzp (22) (2)
zxy将所给速度分布式分别对x和y求偏导数,得
uza2py2x[1()2](2) x4zaa2uz1py2 [1()] (3)22zax2uz1px2 [1()] (4)22zay将式(3)和(4)代入式(2)可知,仅当x给速度分布式不能完全满足运动方程。
5.某一流场的速度向量可以下式表述
u(x,y)5xi5yj 试写出该流场随体加速度向量
解:
2y22a2时才满足运动方程。因此所
DuD的表达式。
DuyDuDuxij
DDD (uuuuuxuuuuxxuyxuzx)i(yuxyuyyuzy)j xyzxyz 25xi[(-5y)(5)]j 25xi25yj
第三章
1. 如本题附图所示,两平行的水平平板间有两层互不相溶的不可压缩流体,这两层流体的密度、动力粘度和厚度分别为1、1、h1和为2、2、h2,设两板静止,流体在常压力梯度作用下发生层流运动,试求流体的速度分布。
解:将直角坐标下的连续性方程和运动方程化简,可得
d2uxdy21px 积分得 upx12xy2C1yC2
因此,两层流体的速度分布可分别表示为 upx112y2C1yC2 1x u1px22xy2D1yD2 2由下列边界条件确定积分常数: (1)yh1,ux10;
(2)yh2,ux20;
(3)y0,ux1ux2;
(4)y0,dux11dydux22dy 将以上4个边界条件代入式(1)与(2),得
1p2h21C1h1C20;
1x
1p2h22D1h2D20;
2x C2D2;
1C12C2
21h221解得 Ch1p2h112x
111h22h11)
2)
((1h221h12ph12p2h12 C2D2
h21x21x1122h12h1212hp1h2 D12
22x12h11h22h12122h2ph2p1h22 D2C2
h22x22x1211h2最后得速度分布方程为
1h221h12py22h12y ux1[21(1)]
1h221xh1h112h12h12122h2py1h22y[21(1)]
h22xh2121h21h2ux22. 粘性流体沿垂直圆柱体的外表面以稳态的层流液膜向下流动,如本题附图所示。试求该流动的速度分布。该液体的密度和粘度分别为和。
由柱坐标系连续性方程
解: 由题给条件,有
0,uru0,Xzg
1(rur)1(u)(uz)0
rrrz简化得
0
z 由柱坐标系N-S方程
uz(uruzuuzuuzz) rrz1uzp12uz2uzg(r)22 2zrrrrz简化得
g1uz(r)0 rrr由于
uzz0,
uz,故uzuz(r),即 0(轴对称)
1dduz(r)0 rdrdrg2积分得 uzrC1lnrC2 (1)
4g边界条件为 (1) (2)
rr0,uz0
rR,duz0 dr将边界条件代入式(1),得 C1g2R 2gr02(Rlnr0) C222故速度分布为
uzg2r122[Rln(r0r)] 2r023. 半径为r0的无限长圆柱体以恒定角速度在无限流体中绕自身轴作旋转运动。设流
体不可压缩,试从一般柱坐标系的运动方程出发,导出本流动问题的运动方程,并求速度分布与压力分布的表达式。
解:柱坐标系的运动方程为
ururuuru2uuruzr r方向: rrrz112ur2u2ur Xr (rur)222 (2-47a)2rrzrrrr1p方向:
uuuuuuuurruz rrrz2211u2uru(2-47b)
X(ru)2222rrzrrrr
11pz方向:
uzuuuuurzzuzz rrz1uz12uz2uz Xz(r)22 (2-47c)2zrrrrz1p由于该流动具有稳态、对称及一维特性,故有
0,uruz0 z利用上述特点,运动方程(2-47)简化为
u2p rr2u1uu20 2rrrr 由于流动为一维,上式可写成常微分方程
2udp (1) drrd2u1duu0 (2) dr2rdrr2式(2)的通解为 利用边界条件 可得
uC1rC2r1
rr0,ur0 r,u0
C10,C2r02
r02因此 u
r如果令
r02
2r则 u压力分布为 由
2p22C
8rr,pp0 可得 Cp0
21因此 pp022
8r4. 试求与速度势2x5xy3y4相对应的流函数,并求流场中点(-2,5)的压力梯度(忽略质量力)。
解:(1)流函数
2x5xy3y4 uxx25yy
2y uy g52 y2g(x)y35xxgx
52 x23xC∴ 2y522 (2)流场中点(-2,5)的压力梯度
忽略质量力,平面稳态流动的Euler 方程为
uxy25x23xC
uxxuyuxy1px
ux写成向量形式为
uyxuyuyy1py
p[(uxuxxuyuxy)i+(uxuyxuyuyy)j]
[(35x)(5)i(25y)(5)j] 5[(35x)i(25y)j] ∴ 点(-2,5)的压力梯度为 p(2,5)(65i115j)
5. 粘性流体在两块无限大平板之间作稳态层流流动,上板移动速度为U1,下板移动速度为U2,设两板距离为2h,试求流体速度分布式。提示:在建立坐标系时,将坐标原点取在两平行板的中心。
解:流体作稳态流动,速度与时间无关。建立坐标系时,将坐标原点取在两平行板的中心,并设两板距离为2h。运动方程可化简为
d2ux1px方向 0 (1) 2xdyy方向 0g将式(2)对y积分得
pgyf(x) (3)
1p (2) y将式(3)对x求偏导数,得
pdf(x) xdx由上式可知,p对x的偏导数与y无关。
x方向的运动方程(1)可改为
d2ux1p (4) dy2x容易看出,上式右边仅与x有关,左边仅与y有关。因此上式两边应等于同一个常数,即
d2ux1pconst 2dyx积分上式得
1py2C1yC2 (5) uxx2边界条件为
(1) yh,uxU1; (2) yh,uxU2 将边界条件代入式(5)得 C1UU21p2U1U2h , C2122x2h于是速度分布式为
UU2yU1U2h2py[1()2]1() ux 2xh2h2第四章
1. 某粘性流体以速度u0稳态流过平板壁面形成层流边界层,在边界层内流体的剪应力不
随y方向变化。
(1)试从适当边界条件出发,确定边界层内速度分布的表达式uxux(y); (2)试从卡门边界层积分动量方程
dx0dux(u0ux)dyduxdyy0
出发,确定x的表达式。
解:(1)由于边界层内duxdy不随y变化,
duxdy为常数,速度分布为直线。设
uxaby。边界条件为
(1)y0,ux0; (2)y,uxu0 由此可得边界层内速度分布为
uxu0y
(2)将边界层积分动量方程写成
dx则
dux0u0ux(1uxu0ux)dyuysu0ddxsu120
dxd0u0(1uxyu0)dy[(1)d]01d6dxsu20
su0y0y()y0yu0 故有
1d 6dxu06u0即 ddx
边界条件为 x0,0,积分上式得 3.464xu0-1/23.464xRex
2. 不可压缩流体稳态流过平板壁面形成层流边界层,在边界层内速度分布为
uxu03y1y3()() 221/2式中,为边界层厚度,4.64xRex。试求边界层内y方向速度分布的表达式uy。
解:
ux3u0x1/2y24.64u0u0x3/2y324.64u03
二维稳态层流的连续性方程为
uxxuyy0 (1)
ux3u0x1y3u0x1y33u0yy3[()] (2) x4434x将式(2)代入式(1)积分,得
ux3u01y21y4 uydy[()()]
x4x24 1.74u02Re1/xy1y[()2()4] 23. 20℃的水以0.1m/s的流速流过一长为3m、宽为1m的平板壁面。试求(1)距平板前缘0.1m位置处沿法向距壁面2mm点的流速ux、uy;(2)局部曳力系数CDx及平均曳力系数CD;(3)流体对平板壁面施加的总曳力。设Rex510。
c5 已知水的动力粘度为100.5105Pas,密度为998.2kg/m3。
解:距平板前缘0.1m处的雷诺数为:
Rexu0x0.10.1998.259932.3210 5100.510流动在层流边界层范围之内。
(1)求y方向上距壁面2mm处的ux,uy
已知 x0.1m,y0.002m,由式(4-15)得 yu0x0.0020.1998.2100.51050.11.993
查表4-1,当1.993时
f=0.6457, f=0.625, f=0.260 由式(4-25)得
uxu0f0.10.6250.0625m/s 由式(4-26)得 uy1u0(f'f)
2x 120.1(100.5105/998.2)0.14(1.9930.6250.6457)
3.0110m/s
(2)局部曳力系数CDx及平均曳力系数CD
12CDx0.664Rex0.664(9932.3)1/20.00666
CD2CDx0.01332
(3)流体对平板壁面施加的总曳力
FdCDu022A0.0133998.20.122130.199N
4. 某粘性流体以速度u0稳态流过平板壁面时形成层流边界层,已知在边界层内流体的速度分布可用下式描述
uxabsincy
(1)采用适当边界条件,确定上式中的待定系数a,b和c,并求速度分布的表达式; (2)试用边界层积分动量方程推导边界层厚度和平板阻力系数的计算式。
解: (1) 选择如下边界条件 (1)y0,ux0; (2)y,uxu0; (3)y,uxy0
代入得
0absin(0)
u0absin(c)
uxbccos(c)0 y求解得 a0;bu0;c故
2uxu0sin(dxdy2)
duxdyy0(2) 0(u0ux)uxdy
先将速度分布代入,求积分号内的项
02(u0ux)uxdyu0(10uxux)dy u0u0uyy)]sin()dy 2212u0[1sin()]sin()d
200[1sin(022u[sin()sin2()]d
02222121u0[cos()(sin)]0
2442u0[2011212(01)0]u0() 222duxdy代入得 y0u0u0y cos22y02u0d221u() 0dx222u0()221dx
d12211.46x11.46x2 u0u0x4.79xRex1/2
sCDxu022u0 2移项得 CDx1/20.65R6ex 1/2u04.79uxRe0x5. 已知不可压缩流体在一很长的平板壁面上形成的层流边界层中,壁面上的速度梯度为kuxyy0。设流动为稳态,试从普兰德边界层方程出发,证明壁面附近的速度分布可用
下式表示
ux1p2xy2ky
式中,p/x为沿板长方向的压力梯度,y为由壁面算起的距离坐标。
证:对于二维平板边界层,普兰德边界层方程为
uxux1p2uxuxuy (1)
xyxy2由于板很长,可以认为 ux0 x由连续性方程
uxuy0 xy得
uyy0
在平板壁面上,uy化为
y00,因此由上式可知,在边界层内uy0。由此可将式(1)简
2uxy21px
上式左端是y的函数,右端是x的函数,二者要相等,必须使得
2uxy21px=常数
上式积分求解,得
uxy1px1pyC1
ux2xuxyy2C1yC2
由题意,当y0时,k,故
C1k
又当y0时,由壁面不滑脱条件,ux0,故 C20 因此,速度分布为 ux证毕。
6. 不可压缩流体以u0的速度流入宽为b、高为2 h的矩形通道(b,从进口开始a)
1p2xy2ky
形成速度边界层。已知边界层的厚度可近似按 5.48x/u0估算,式中x为沿流动方向的距离。试根据上述条件,导出计算流动进口段长度Le的表达式。
解:当h(矩形高度的一半)时,边界层在通道的中心汇合,此时的流动距离x即为流动进口段长度Le,故
h5.48解得
Leu0
u0h2h2()0.033 Le 5.48u0Le0.033Re huh式中 Re0
或
第五章
1. 20℃的水在内径为2m的直管内作湍流流动。测得其速度分布为ux100.8lny,在离管内壁1/3 m处的剪应力为103Pa,试求该处的涡流运动粘度及混合长。
已知20℃水的密度为998.2 kg/m3,动力粘度为1.005×10-3Pas。 解:(1)涡流运动粘度
yx()tdux (1) dy0.8y0.80.3332.4s-1 (2)
duxdyddy(100.8lny)式(2)代入式(1)并整理得 tyxduxdyy0.333103998.22.41.005103998.2
4.30102m2/s
362已知/1.00510/998.21.00710m/s。可见,离管内壁1/3处的粘性扩散系数与涡流扩散系数相比,可以忽略不计。 (2)混合长
忽略粘性应力,则 yxyxl(rt2duxdy)2
tyxl2(du/dy)x1/2(103998.22.41/2)0.1338m 2其值约为管半径的13.4%。
2. 温度为20℃的水流过内径为50mm的圆管。测得每米管长流体的压降为1500N/m2,试
证明此情况下流体的流动为湍流,并求
(1)层流底层外缘处水的流速、该处的y向距离及涡流粘度;
(2)过渡区与湍流主体交界处流体的流速、该处的y向距离及涡流粘度。
解:由物性表查得 20℃水的物性:=998.2kg/m3,100.510Pas
3
sup2Lri150020.05/218.75N/m2
s18.750.137m/s 998.2rui
ubu(2.5ln1.75)
0.0250.137ln5100.510998.2.751
7 0.132.53 .02m/s
Reddub0.253.02998.2100.51051.51051104,故为湍流。
(1)层流内层外缘处 ∵ uy5
∴ uuy5u50.1370.685m/s
yuy
yb
5u5u5100.5105998.20.1373.67105m
0(层流内层只有粘性力)
u5.0lny3.05
(2)过渡区与湍流主体交界处
5.0ln30
=13.956
03.6u u13.95 ym=13.9650.1371.
30u=6b2.2104m
422102
s(1)=18.58N/m 18.75(1)ri0.025y∵ du dy∴ du1() dy将u5.0lny3.05写成
uu5.0lnyuyu3.05
或 u5.0uln3.05u
du5.0u dyy
5.0u1y()
y5.0u18.58 998.2 6.0106
22410
5.00.137m2/s
3. 试应用习题7中的己知数据,求rri/2处流体的流速、涡流粘度和混合长的值。
解:yyuriu20.0250.137998.22100.5106
3为湍流主体区0 1.710 u2.5lny5.52.5ln1.71045.5
4 29.85
uuu29.850.1374.09m/s
s(1)s(1ri
yri/21)s ri2
1218.759.375N/m2
du2.5u同上题方法推导 ,故 dyydu19.375y9.3750.025/2()3.43104m2/s dy998.22.5u998.22.50.137
l2du dy0.025/2ydu433.43103.5410m )2.50.137dy25u
l(4. 利用流体阻力实验可估测某种流体的粘度,其方法是根据实验测得稳态湍流下的平均速度ub及管长为L时的压降(p)而求得。试导出以管内径d、流体密度、平均流速ub和单位管长压降p/L表示的流体粘度的计算式。
解: 设流体在管内湍流且流动充分发展,则
6e f0.04R而 p4f0.2Lub2d240.046(dub)0.2Lub2d2
1.80.80.0920.2d1.2ubL
移项并整理得
d1.2(p/L) 0.81.80.092ub5. 假定平板湍流边界层内的速度分布可用两层模型描述,即在层流底层中,速度为线性分布;在湍流核心,速度按1/7规律分布,试求层流底层厚度的表达式。
解: 层流底层很薄,故有
s5duxdyuxy (1)
平板湍流壁面剪应力又可由式(5-56)表示,即
u s20.0230u0以上两式联立得到
1/4 (2)
u02u01/4()y (3) ux0.023 令层流底层厚度为b,其外缘与湍流核心接壤处的速度为ul,则上式可写成
u02u01/4()b (4) ul0.023 或写为
b1ul3/4()() (5) 0.023u0u0另一方面,湍流核心的速度可用1/7次方定律描述,即
uxu0y()1/7
在两层交界处,有
ulu0或写成
(b1/7) bu(l)7 (6) u0式(5)与式(6)联立得
ulu01.875(1/8)1.875Reb1/8 u0将
xxu0.376(01/5) 代入上式,可得
ulu02.12(u0x1/10)1/102.12Rex (7)
将式(7)代入式(6)中,得
b7/10 192.5Rex1/5再将
xxu0.376(0) 代入上式,可得
bx(b)()(192.5Rex7/10)(0.376Rex1/5)72.4Rex9/10 x6. 20℃的水在内径为2m的直管内作湍流流动。测得其速度分布为ux100.8lny,在离管内壁1/3 m处的剪应力为103Pa,试求该处的涡流运动粘度及混合长。
已知20℃水的密度为998.2 kg/m3,动力粘度为1.005×10-3Pas。 解:(1)涡流运动粘度
yx()tdux (1) dy
duxdyddy(100.8lny)0.8y0.80.3332.4s-1 (2)
式(2)代入式(1)并整理得 tyxduxdyy0.333103998.22.41.005103998.2
4.30102m2/s
362已知/1.00510/998.21.00710m/s。可见,离管内壁1/3处的粘性扩散系数与涡流扩散系数相比,可以忽略不计。 (2)混合长
忽略粘性应力,则 yxyxl(rt2duxdy)2
tyxl2(du/dy)x1/2(103998.22.41/2)0.1338m 2其值约为管半径的13.4%。
7. 试从光滑圆管中湍流核心的对数速度分布式(5-43)和剪应力与y的关系式
ys(1)出发,推导涡流粘度的表达式,并讨论涡流粘度或参数与雷诺数Re和
riuri位置y的关系。
解:在湍流核心,可忽略粘性力的影响,则
duyr s(1)x
ridy移项得
sydu(1)(x)1 (1) ridy湍流核心的速度分布为
u u2.5lny5.5, ux2.5duxdyuyyuln5.u 5 2.5 (2)
将式(2)代入式(1),并将us/代入,得
u(1)或写成
u(122yyri2.5u0.4u(1)y (3)
riyuriyy0.4()(1) ri2.5uririy)y第七章
1. 在一无内热源的固体热圆筒壁中进行径向稳态导热。当r11m时,t1200℃,r22m时,t2100℃,其导热系数为温度的线性函数,即
kk0(1t)
式中k0为基准温度下的导热系数,其值为k00.138W/(mK),为温度系数,其值为1.951041K,试推导出导热速率的表达式并求单位长度的导热速率。
解:导热速率的表达式为
qt k (1)
Ar对于本问题,kk0(1t),A2rL,其中L为圆筒壁的长度。稳态热传导情况下,,并分离变量,可得 qconst,将上述情况代入式(1)
(1t)dt边界条件为
qdr (2)
2Lk0r①r11,t1200 ②r22,t2100
积分得
1qtt2lnrC 22Lk0将边界条件①代入式(2),可得
1qCt1t12lnr122Lk0
1(200273.2)(1.95104)(200273.2)2495 K
2于是导热速率方程为
q2Lk0(t121t495) (3) 2lnr将边界条件②代入式(3),可得单位长度的导热速率为
q/L2k0(t211t22495)2lnr211{23.140.138[(100273.2)1.95104(100273.2)2495]}2ln2135.3 W/m
的球形固体,其半径为R 。球体沿径向向外对称导热。球2. 有一具有均匀发热速率q表面的散热速率等于球内部的发热速率,球表面上维持恒定温度tS不变。试从一般化球坐标系热传导方程出发,导出球心处的温度表达式。
解:球坐标系的热传导方程为式(7-3),即
1t12t1t12tq2(r)2(sin)22 2rrrrsinrsink① 球表面的散热速率等于球内部的发热速率,球表面上维持恒定温度ts不变,故为稳
态导热,
t0 2ttt2t0,20,20 ② 一维径向对称导热, 0,于是式(7-3)变为
1d2dtq(r)0 (1)
drkr2dr边界条件为
①rR,ttS ②rR,
43dt4R2k Rq3dr式(1)积分两次,得
tq2C1 rC2 (2)
6krq2R 6k将边界条件①、②分别代入式(2)可得
C10,C2ts于是球体内的温度分布方程为
q2q2 rtsR (3)
6k6k令式(3)中的r0,即得球心处的温度表达式,即
q2tr0tsR
6kt3. 将厚度为0.3 m的平砖墙作为炉子一侧的衬里,衬里的初始温度为30 ℃。墙外侧面绝热。由于炉内有燃料燃烧,炉内侧面的温度突然升至600℃并维持此温度不变。试计算炉外侧绝热面升至100℃时所需的时间。已知砖的平均导热系数k =1.125W/(mK),导温系数
5.2107m2/s。
解:将平砖墙视为一维导热,则其温度分布式为
T式中
tts4133exp[()2Fo]cos(L*)exp[()2Fo]cos(L*)……t0ts22322 Fol2*, Lx, ll m0.3*对于此题,L0, t100, ts600, t030 。
为了计算方便计,只取级数的第一项,即
1006004 772exp[(2Fo)]0.8306002解之,Fo0.151
132此时,级数的第二项 exp[()Fo]0.01170.8772
32
所以只取级数的第一项是合理的。
由Fol20.151,得
0.1510.32261347.26hs 75.210Fol24. 有一半径为25 mm的钢球,其导热系数为433 W/(mK),密度为7849 kg/m3,比热容为0.4609 kJ/kg,钢球的初始温度均匀,为700 K。现将此钢球置于温度为400 K的环境中,钢球表面与环境之间的对流传热系数为11.36W/(m2K)。试求1小时后钢球所达到的温度。
解:由于h值较小,k值较大,估计可以采用集总热容法,为此首先计算Bi数。
4311233 V/Ar/4r25108.310000333h(V/A)11.368.31032.21040.1 Bik433故可应用集总热容法求解。
Fo(V/A)2k 2cP(V/A) 43336006.255103 327849460.9(8.310)43t400e2.2106.25510
7004000.253
t0.253(700400)400
475.8K
因篇幅问题不能全部显示,请点此查看更多更全内容