爆轰学 _冲击波导论

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,爆炸科学与技术国家重点实验室,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,1,第,2,章,冲击波导论,2,第,2,章,冲击波导论,需要首先了解气体的流动及有关波的知识。,炸药爆炸会形成高温、高压的气体。,气体的膨胀过程就是对外作功的过程。,3,第,2,章,冲击波导论,本章内容：,2.1 波的根本概念,2.2 气体的平面一维流动,2.3 平面正冲击波,2.4 冲击波的波速线、Hugoniot曲线和等熵线,2.5 冲击波的根本性质,2.6 冲击波的正反射,2.7 弱冲击波的声学近似理论,4,波的根本概念,5,波的根本概念,1、波Wave,波通常可以分为两大类：一类是电磁波，另一类是机械力学波。,当介质Medium受到外界作用如振动、冲击等时，介质的局部状态参量就会发生变化，这就是扰动Disturbance。,6,波的根本概念,如果活塞突然向右移动，便有波向右传播。,在扰动传播过程中，扰动介质与未扰动介质之间存在一个界面，这个界面就叫波阵面(Wave front)。,扰动在介质中的传播速度,叫做波速(Wave velocity)。,要与介质的质点速度区分,7,波的根本概念,如果扰动前后介质的状态参数变化量与原来的参数量相比是很微小的，那么称这种扰动为弱扰动Weak disturbance或小扰动。弱扰动的特点是各种参数的变化量是微小的、逐渐的和连续的。,如果扰动前后介质的状态参数发生突跃变化，那么称这种扰动为强扰动Strong disturbance。,8,波的根本概念,图,2,3,弱扰动和强扰动波形,9,波的根本概念,2、声波sound wave,声波是一种弱扰动波。弱扰动在介质中的传播速度就叫声速。它是气体动力学中一个非常重要的参数。,下面以活塞在直管中移动,所引起的气体扰动的传播,来建立声速c与其它参数,的关系式。如下图。,10,波的根本概念,(1),在,t,0,时刻，活塞处于静止状态，状态参数为,(2),在,t,1,时刻，活塞运动到,B-B,处，扰动传到,D-D,处，弱扰动传过后，状态参数变为 ，质点速度变为,u,。,11,(1),式中x为t1时刻扰动传播的距离，x=ct1,x1为时刻活塞运动的距离，x1=ut1,A0为活塞的截面积。,代入1式可得：,消去t1后可得： 2,波的根本概念,质量守恒Conservation of Mass ：,12,波的根本概念,动量守恒Conservation of Momentum ：气体受到扰动后的动量等于作用在其上面的冲量。,化简后得： (3),2式代入3式得：,(4),由2式可得： (5),13,波的根本概念,把5式代入4式得：,(6),由于声波为弱扰动波，波阵面过后介质状态变化为一微小量，故有 ，因此，6式变为：,(7),看作等熵过程： (8),14,波的根本概念,对于理想多方气体，其等熵方程为：,(9),那么 (10),所以理想气体的声速为：,(11),又由 可得 ：,(12),15,波的根本概念,对于地外表上的空气，可近似地视为理想气体，将,，代入上式可得：,(13),将 代入13式可得 340m/s。,16,波的根本概念,需要指出的是,，只有对于小扰动， 才成立，扰动才以声速传播。对于 的扰动，其传播速度大于声速，扰动越强，传播速度将越高。,17,波的根本概念,3,、,压缩波和稀疏波,压缩波Compression Wave：扰动传过后，介质的压力、密度、温度等状态参数增加的波称为压缩波。其特点是波传播的方向与介质质点运动方向相同。,稀疏波Rarefaction Wave：扰动传过后，介质的压力、密度、温度等状态参数下降的波称为稀疏波，其特点是波传播的方向与介质质点运动方向相反。,18,波的根本概念,在一个连续的，缓慢的压缩过程中，每一小步的压缩都是一种等熵变化，但由于每经一步压缩后气体的温度都要上升，气体的声速必将上升，这样下一步的压缩波的波速逐渐增加，一旦集中起来，状态参数的变化将不再连续，就会发生突跃，弱扰动变成,强扰动,。,19,波的根本概念,由于稀疏波的膨胀飞散是按顺序连续进行的，所以稀疏波传播中介质的状态变化是连续的，如图,2,4,中的压力变化。,图,2,4,稀疏波现象,20,波的根本概念,在稀疏波扰动过的区域中，任意两相邻端面的参数都只差一个无穷小量，因此稀疏波的传播过程属于等熵过程，它的波速等于介质当地的声速或音速Local sound speed。,21,2.2,气体的平面一维流动,22,2.2,气体的平面一维流动,所谓一维流动，是指在某一空间坐标,x,等于常数的平面上流体参数都是均匀分布的，并且在给定坐标,x,处的流体参数都只随时间,t,变化的流动，如,p,=,p,(,x,t,),，,T,=,T,(,x,t,),，,u,=,u,(,x,t,),，,=,(,x,t,),。,另外在球坐标系中，中心对称流动问题也是一种一维流动，,u,=,u,(,r,t,),，其它参数也只是,r,和,t,的函数。,一维流动又可分为,一维定常流动,和,一维不定常流动,。,23,气体一维流动的根本方程组,24,气体一维流动的根本方程组,气体在平面一维流动下，满足质量守恒、动量守恒和能量守恒，其对应的方程分别叫质量方程连续方程、动量方程欧拉方程、能量方程。,1、连续方程质量方程,(1),该式为一维不定常流动的连续方程。,2、欧拉方程动量方程,(2),25,气体一维流动的根本方程组,3,、,能量方程,在不考虑气体的粘性和热传导的情况下，气体的流动是等熵的。,(3),4,、,状态方程,由于,S,可表示为,p,和 的函数，故等熵流动条件可表示为：,对于理想气体，其等熵方程为：,(,4,),这样，便可由连续方程、欧拉方程、能量方程和状态方程求解气体一维等熵流动的四个未知量 。,26,2.2.2,以,u,、,c,为求解参量的方程组,27,2.2.2,以,u,、,c,为求解参量的方程组,为使前面建立起来的气体一维等熵流动的方程组的物理意义更容易理解，将它们稍加变换。引入声速c代替p和 。,由声速公式及等熵方程可得：,(5),将5式两边微分并同时除以 ，得,(6),28,2.2.2,以,u,、,c,为求解参量的方程组,由5式知，,(7),把6式代入7式，可得：,(8),29,2.2.2,以,u,、,c,为求解参量的方程组,将6式代入连续方程1式，可得,(9),将8式代入欧拉方程2式，可得,(10),30,2.2.2,以,u,、,c,为求解参量的方程组,将9、10两式相加和相减，整理可得,(11),这个方程组即是以u、c为变量描述气体一维等熵不定常流动规律的方程组。,确定气体一维等熵流动过程中气体各参数时的时间、空间变化规律，归结为解此偏微分方程组。,31,2.2.2,以,u,、,c,为求解参量的方程组,小扰动波在静止介质中是以音速进行传播的，在一维情况下，静止气体中小扰动波的传播速度为,c,。在流动介质中，小扰动波的传播速度为介质流动,速度,u,与,当地,音速,c,的,叠加，即 。顺介质流动方向传播的扰动取正号，逆介质流动方向传播的扰动取负号。,32,2.2.2,以,u,、,c,为求解参量的方程组,在 条件下，11式可表示为 对t的全导数形式，并且该导数为零，即,(12),即,(13),33,2.2.2,以,u,、,c,为求解参量的方程组,由此可以看出，方程11在 条件下描述,的是两个量的推进规律:,即由 所确定的状态或扰动以速度,顺气体流动方向即x轴的正方向传播；,而由 所确定的状态或扰动以速度,逆气体流动方向传播。,34,方程组的特征线及一般解,35,方程组的特征线及一般解,dx/dt=u+c和dx/dt=u-c分别代表一维等熵流动介质中扰动沿x轴的正向和反向传播的速度，我们称它们为11式的特征或特征方程。它们的积分各自代表xt平面上的一簇曲线，叫做特征线。其中在xt平面上由dx/dt=u+c所确定的特征线称为第一簇特征线，用C表示；而由dx/dt=u-c所确定的特征线称为第二簇特征线，用C表示。,36,方程组的特征线及一般解,这两簇特征线分别描述的是物理状态量 ，即扰动波以速度 沿x轴的正向或负向传播的轨迹。,因此，对于一维等熵不定常流动方程组11式，有,沿着C特征线,(14),37,方程组的特征线及一般解,沿着C特征线,(15),式中，I+,I-称为黎曼Riemann不变量。,它们在u，c平面上可用两簇相互平行的直线来描述，称为方程组11在速度平面上的特征线。它们在沿着各自的特征线C和C传播时保持不变。如图25所示。,38,方程组的特征线及一般解,图,2,5,特征线,39,方程组的特征线及一般解,方程14和15为方程组11的一般解。,在k3的最普通的情况下，由于(u+c)和(u-c)都是x和t的函数，即右传波的传播速度受反方向波的影响，因此14式和15式无法得到精确的解析解。一般采用数值积分法或特征线法近似求解。,40,方程组的特征线及一般解,在x-t平面上,假设曲线AB上的各点状态参数。C和C分别表示AB线上各点发出的不同簇的特征线，求解流场D内各点的状态参数。,图,2,6,特征线法解流场参数,41,方程组的特征线及一般解,【解】近似认为 ， ， 并近似把x-t平面上特征线的一小段视为直线。,在曲线AB上选取一系列的点M1，M2，Mi等。由于ui和ci，过Mi点作特征线C，其特征方程为：,42,方程组的特征线及一般解,而过,M,i+1,点作特征线,C,，其特征方程为：,其中,x,和,t,为,C,和,C,相交点,M,i,空间坐标和时间，,u,和,c,为该相交点,M,i,的状态参量。同样可以求出,A,B,上各点的参量，依次可求出任意位置点的状态参量。,43,2.2.4,方程组的特殊解,简单波流动,44,2.2.4,方程组的特殊解,简单波流动,前面讨论的14和15式是方程组11式的通解，流场中可同时存在左传波和右传波。如果流场中只有一个方向传播的扰动波，即波未进入的区域介质处于静止状态或稳定流动状态，这种波就称为简单波，其解称之为方程组的特殊解。,简单波：它的某一族特征线上的黎曼不变量是同一个常数，即该族各条特征线上的黎曼不变量彼此相等。,45,2.2.4,方程组的特殊解,简单波流动,当给定如下条件，即,(16),对16式分别对t和x求偏导，得,将这两式代入11式可得,(17),46,2.2.4,方程组的特殊解,简单波流动,该式说明，沿特征线dx/dt=u+c，有du/dt=0。,即u常数。,由16式知，c亦为常数。因此dx/dt=u+c就可以积分了。,因此,(18),47,2.2.4,方程组的特殊解,简单波流动,同理，当 时，有,(19),式中， 是u的任意函数，由边界条件确定。,由18和19式即可确定简单波的向前波右传波和向后波左传波流动区内任一点的参数u和c。,48,2.2.4,方程组的特殊解,简单波流动,为了说明简单波的性质，我们来考察下面两种情况。,1活塞向左加速运动，如下图。,图,2-7,右传系数波,49,2.2.4,方程组的特殊解,简单波流动,当活塞向左加速拉动时，便形成一系列的简单稀疏波向右传播，并以当地声速传播，因此，活塞向左拉动时发出的第一道稀疏波是以静止气体当地的音速u0+c0=c0的速度向右传播的，特征线如图2-7所示。,该特征线的右边为静止气体区域，故该区域内的特征线也都是平行的。活塞向左加速拉动而发出的各个后续右传稀疏波，是在扰动过的气体中传播的，因此第n道波的传播速度un+cn总是比其前面的的n-1道波的传播速度un1+cn1要慢。因此后面的各道波的特征线C是发散的。,50,2.2.4,方程组的特殊解,简单波流动,2活塞向右渐渐加速运动，如图2-8所示。,图,2-8,压缩波随,t,的变化,51,2.3,平面正冲击波,52,2.3,平面正冲击波,冲击波Shock wave，又称激波，是一种强烈的压缩波，其波阵面通过的前后参数变化很大，它是一种状态突跃变化的传播。,冲击波阵面Shock front实际上有一定的厚度，其厚度约为几个分子平均自由程，在这个厚度上各物理量发生迅速的、但却是连续的变化，这是由于物质具有粘性和热传导的原因。但在工程计算上可以不考虑粘性和热传导等耗散效应，而将冲击波视为一个没有厚度的间断面。因此，可以说冲击波阵面是一种强间断面。,53,2.3,平面正冲击波,人员伤害超压准那么,超压,(MPa),0.020.03,0.030.05,0.050.1,0.1,伤害程度,轻微,中等,严重,极严重，死亡,超压计算的,Backer,公式：,其中,P,MPa,；,W,kg,；,R,m,。,54,1kg TNT,炸药在空气中爆炸：,2.3,平面正冲击波,距离,(m),1,1.1,1.3,1.4,1.5,1.6,1.7,1.8,1.9,超压,(MPa),0.78,0.62,0.42,0.35,0.30,0.26,0.23,0.20,0.18,距离,(m),2.0,2.1,2.2,2.3,2.4,2.5,3,4.0,5.0,超压,(MPa),0.16,0.14,0.13,0.12,0.11,0.1,0.07,0.04,0.028,55,2.3,平面正冲击波,对人员的杀伤：,【,文献来自：,Elements of Terminal Ballistics Part 2, Collection and Analysis of Data Concerning Targets】,其中,P50,为造成,50%,死亡所需的超压，单位,psi(lb/inch,2,),；,W,为体重，单位克。,1psi=6895Pa,1lb=0.45359kg,56,2.3.1 根本关系式,57,2.3.1 根本关系式,设有一冲击波以恒定的速度向右传播，如图,2-9,所示。,图,2,-,9,平面正冲击波阵面,58,2.3.1 根本关系式,波的右边，尚未扰动的介质，参数为： 。,波的左边，扰动的介质，参数为： 。,为方便起见，把坐标系建立在波阵面上。那么未扰动的介质以的速度向左流入冲击波阵面，扰动的介质以的速度从波阵面流出。,1、质量守恒Conservation of Mass ：单位时间内流入波阵面的质量等于流出的质量。,即： (1a),将 ，上式变为： (1b),59,2.3.1 根本关系式,2、动量守恒Conservation of Momentum：单位时间内作用介质上的冲量等于其动量的改变。,冲量：,动量变化：,因此 (2a),即 (2b),60,2.3.1 根本关系式,3能量守恒Conservation of Energy：冲击波传播视为绝热过程，忽略介质的粘性和热传导效应等能量耗散。,单位时间内从波阵面右侧流入的能量包括有：,内能,介质压力和流入的介质体积所确定的压力位能,介质流动的动能,61,2.3.1 根本关系式,同理，从波阵面流出的能量为：,内能,介质压力和流入的介质体积,介质流动的动能,62,2.3.1 根本关系式,因此, , ,整理后可得：,(3),以上三个式子(1)、2和3即为冲击波的根本关系式。,63,2.3.1 根本关系式,为便于使用，将1、2、3式进行变换。,将1a、2a式联立消去D-u0)可得,(4),将4式代入1b式，可得,(5),(5式即为冲击波波速方程Rayleigh，瑞利方程。,64,2.3.1 根本关系式,把2式变为：,(6),把6式代入3式可得：,(7),把4式代入7式可得：,(8),8式就是著名的雨贡纽Hugoniot方程，又称冲击绝热方程。该方程适用于任何介质中传播的冲击波。,65,2.3.1 根本关系式,其中，4、5和8式为冲击波的三个根本关系式。,对于某一具体介质中传播的冲击波，需与该介质的状态方程联系起来， 或 以便求解冲击波阵面上的参数。这样，四个方程就有了五个参数：,66,2.3.2,多方气体中的平面正冲击波,67,2.3.2,多方气体中的平面正冲击波,对于多方气体，其内能可表示为：,(9),其中： 定容比热容； 气体的多方指数假定不变。,把9式代入Hugoniot方程，可得：,(10),68,2.3.2,多方气体中的平面正冲击波,整理可得：,(11),(12),12式和4式、5式联立，并结合 ，可得：,(13),69,2.3.2,多方气体中的平面正冲击波,如果未受扰动气体静止时,(14),70,2.3.2,多方气体中的平面正冲击波,因此，只要 任意一个参数就可以就算其余参数。,对于强冲击波， ，,(15),71,2.3.2,多方气体中的平面正冲击波,对于强冲击波，波阵面上的质点速度与冲击波速度成正比；压力与冲击波速度的平方成正比；对于，波阵面上的密度最大可达初始密度0的6倍。,假设引入马赫数Mach number,(16),那么13式可写成,(17),72,2.3.2,多方气体中的平面正冲击波,【例】测得空气中爆炸产生的冲击波的D1000m/s，计算其参数 ，初始状态 ， ，,， ， 。,解：1求,2求,Pa,73,2.3.2,多方气体中的平面正冲击波,3求,4求,因此，,74,2.3.2,多方气体中的平面正冲击波,(5),求,T,由理想气体状态方程,可知，,因此，,2.3.2,多方气体中的平面正冲击波,1超压overpressure,2动压dynamic pressure ：也称之为速压。,物体在流体中运动时，在正对流动运动方向的外表，流体完全受阻，此处的流体速度为0，其动能转变为压力能，压力增大，其压力称为全受阻压力简称全压或总压，它与未受扰动处的压力即静压之差，称为动压。,运动流体密度和速度平方积之半,1/2v2。,76,2.3.2,多方气体中的平面正冲击波,【作业】,1、定量比较超压与动压的大小。,2、实验测得空气中某处温度为15摄氏度爆炸冲击波的超压为，试计算该处的D，u和,77,2.4,冲击波的波速线、,Hugoniot,曲线和等熵线,78,2.4,冲击波的波速线、,Hugoniot,曲线和等熵线,1、波速线Rayleigh线，瑞利线,冲击波波速方程：,(1),设冲击波波前介质是静止的，即,那么1式可变为：,或 (2),79,2.4,冲击波的波速线、,Hugoniot,曲线和等熵线,显然，在 坐标平面内，当D一定时，2式代表一条通过初态O点的直线。不同的D对应不同的斜率，这些斜线称之为波速线或Rayleigh线瑞利线，如图210所示。,图,2,10,冲击波的波速线,80,2.4,冲击波的波速线、,Hugoniot,曲线和等熵线,波速线的物理意义：当 一定时，冲击波通过任何介质后，波后状态都对应于此条线上的某一确定点。因此，通过 点的某一波速线乃是一定波速的冲击波传过具有同一初始状态点 的不同介质所到达的终点状态的连线。,81,2.4,冲击波的波速线、,Hugoniot,曲线和等熵线,2、Hugoniot曲线冲击绝热线,冲击波的冲击绝热方程：,(3),82,2.4,冲击波的波速线、,Hugoniot,曲线和等熵线,在 坐标平面上可以用一条以介质初态 为始发点的曲线来描述。如图211a中的曲线。该曲线称之为冲击绝热线或Hugoniot曲线。,a b,图211 冲击波的冲击绝热线,83,2.4,冲击波的波速线、,Hugoniot,曲线和等熵线,对于多方气体，,那么有：,当 时， (4),即Hugoniot曲线的渐近线是,84,2.4,冲击波的波速线、,Hugoniot,曲线和等熵线,Hugoniot曲线是一条通过初始点的曲线，对某一确定的介质而言，不同的 对应不同的曲线。,当介质性质和波前状态一定时，H线是确定的，假设冲击波速度不同，那么波后状态必然处在H线的不同位置上，如图211a所示。,85,2.4,冲击波的波速线、,Hugoniot,曲线和等熵线,当具有相同波速的冲击波在具有同一初始状态的不同介质中传过后，由于不同介质的H线不同，因此所到达的波后状态将对应于R线上的不同点，如图212所示。,图,2,12,86,2.4,冲击波的波速线、,Hugoniot,曲线和等熵线,因此可以看出，冲击波的H线是不同波速的冲击波在具有同一初始状态的相同介质中传过后所到达的终态点的连线。物理意义,波速线是一定波速的冲击波传过具有同一初始状态的不同介质所到达的终态点的连线。 物理意义,这两条线上的任一点都是和一定的波后状态对应的，它们都不是冲击压缩的过程线，不能认为冲击压缩过程是沿着这两条线中的任一条进行的。,87,2.4,冲击波的波速线、,Hugoniot,曲线和等熵线,【作业】具有相同初始状态的介质，一次冲击压缩和屡次连续冲击压缩能否到达相同的终点状态？为什么？,88,2.4,冲击波的波速线、,Hugoniot,曲线和等熵线,3、等熵线Isentropic curve,前面讲到，一切弱扰动波都以当地声速进行传播的，并且传播过程是等熵的。,对于理想气体，等熵条件下的状态变化遵循等熵方程 所确定的规律，即,89,2.4,冲击波的波速线、,Hugoniot,曲线和等熵线,等熵线就是由等熵方程确定的曲线，它表示进行等熵压缩或等熵膨胀过程时介质状态变化所走过的路径。因此，等熵线是状态变化的过程线。图,2,13,是由初始状态 发生等熵压缩和等熵膨胀过程时的状态变化路径。,图,2,13,等熵线,90,2.4,冲击波的波速线、,Hugoniot,曲线和等熵线,4、H线和S线的关系,1Hugoniot曲线不是状态变化的曲线，而等熵线是一系列微弱扰动波传过后介质状态变化所经历的过程线或路径。,2为说明冲击Hugoniot曲线和等熵线之间的关系，我们以多方气体为例，假假设将该气体从 状态压缩到同样的压缩程度，分别按冲击绝热压缩和等熵压缩进行计算所得的数值列于下表：,91,2.4,冲击波的波速线、,Hugoniot,曲线和等熵线,表21 气体冲击绝热压缩与等熵压缩参数的比较,压缩程度,1.0,0.8,0.6,0.4,0.3,1/6,冲击压缩,1.0,1.368,2.08,4.0,7.125,等熵压缩,1.0,1.366,2.044,3.61,6.31,12.3,92,2.4,冲击波的波速线、,Hugoniot,曲线和等熵线,把表中的数据画在,p-v,平面上，就可知过初始点的等熵线位于过该点的冲击,Hugoniot,曲线的坐下方，且在,O,点相切，如图,2,14,所示。,图,2,14 Hugoniot,曲线和等熵线的关系,93,2.4,冲击波的波速线、,Hugoniot,曲线和等熵线,3Hugoniot曲线上各状态点都在等熵线的上方，因此Hugoniot曲线上的各状态点的熵都大于S0，即冲击波阵面传过后介质的熵是增加的。并且沿Hugoniot曲线，熵随介质的压力增大而增大。,94,2.5 冲击波的根本性质,95,2.5 冲击波的根本性质,1.,冲击波阵面是一个间断面；,2.,冲击波是压缩波，不可能是稀疏波；,3.,冲击波传过后，介质的熵是增加的；,4.,冲击波相对波前介质是超音速的，即,96,2.5 冲击波的根本性质,这个结论可由 证明,也可用,Hugoniot,曲线和等熵线之间的关系证明，如图,2,15,所示：,图,2,15,97,2.5 冲击波的根本性质,【,证明,】,：设冲击波的波速为,D,，介质初始状态为,由波速方程,知,即,(1),98,2.5 冲击波的根本性质,由声速公式 知,(2),即,由图,2,15,中,Hugoniot,曲线和等熵线的关系知，,99,2.5 冲击波的根本性质,即,因此 。证毕。,100,2.5 冲击波的根本性质,5.,冲击波传过后介质获得了一个与波传播方向相同的移动速度，即,这个结论可由,得以证明。,101,2.5 冲击波的根本性质,6.,冲击波相对波后介质是亚音速的，即,【,证明,】,：,对,Hugoniot,方程 两边微分得：,(1),由热力学定律知：,(2),102,2.5 冲击波的根本性质,将2式代入1式得：,(3),103,2.5 冲击波的根本性质,而声速c按定义可表示为：,(4),且 (5),把4式和5式代入3式可得：,(6),104,2.5 冲击波的根本性质,由于冲击波沿,Hugoniot,曲线，熵随介质的压力增大而增大，因此有,因此：,即 证毕。,105,2.6,冲击波的正反射,106,2.6,冲击波的正反射,当冲击波在传播过程中遇到障碍物时，会发生发射现象。当入射波传播方向恰好垂直于障碍物的外表时，发射的反射现象称为正反射现象。,下面讨论多方气体中传播的平面冲击波在刚性壁面上的正反射现象。,107,2.6,冲击波的正反射,设有一稳定传播的平面冲击波以D1的速度向刚体壁面垂直入射。如图216a所示。,图,2,16,冲击波在刚壁面上的正反射,108,2.6,冲击波的正反射,入射波阵面前的状态：,入射波阵面后的状态：,反射波阵面前的状态：,反射波阵面前的状态：,109,2.6,冲击波的正反射,反射前冲击波阵面前后的参数间关系为：,(1),(2),(3),110,2.6,冲击波的正反射,当入射波阵面碰到刚壁面时，由于刚壁面不变形，那么波阵面后气体流的速度立即由u1变为零。就在这一瞬间，速度为u1的气体介质的动能便立即转化为静压势能，从而使壁面处的气体压密，密度由突增为，压力由p1突跃为p2，比内能由e1突跃为e2。由于p2p1，21，受到第二次冲击压缩的气体必然反过来冲击压缩已被入射波压缩过的气体，这样就形成反射冲击波远离刚体壁面向左传播，如图216b所示。,111,2.6,冲击波的正反射,由于反射冲击波在已受入射冲击波压缩过的气体介质中传播，故传过后介质的参数间的关系可表示为：,(4),(5),(6),112,2.6,冲击波的正反射,假设 ，而且由刚壁条件 知，所以由2式和5式可得：,(7),两边平方后整理可得：,(8),113,2.6,冲击波的正反射,将3和6式代入8式可得：,(9),此即反射冲击波阵面压力与入射冲击波阵面压力之间的关系。,式9也可写成压差的表达形式，即：,9,114,2.6,冲击波的正反射,当入射冲击波压力很高时，p1p0 ，可忽略p0，那么9、9可变为：,10,对于空气中的强冲击波来说，如将k值代入，那么有：,当入射冲击波很弱时，由式9可得：,115,2.6,冲击波的正反射,将9式代入6式可得：,(11),对于强冲击波，忽略P0，那么式11为：,12,当强冲击波在固壁反射后，也就是介质经过入射和反射冲击波的两次压缩后，固壁面附近的介质被压缩的最大倍数可由式6和式12求出，即,13,116,2.6,冲击波的正反射,对于空气中的强冲击波反射，有,在,u,0,=u,2,=0,的情况下，入射冲击波和反射冲击波的动量守恒方程可写为：,两式相除，可得：,(14),117,2.6,冲击波的正反射,把式3、6、9、9代入式14整理可得：,(15),当入射冲击波很强，即p1p0时，上式可简化为,对于空气中的强冲击波来说，有：,由此可知，反射冲击波的传播速度总是低于入射冲击波的传播速度，而且两波的方向相反。,118,斜反射,119,斜反射,120,2.7,弱冲击波的声学近似理论,121,2.7,弱冲击波的声学近似理论,前面提到，冲击波传过后，介质的熵增加。但对于弱冲击波 ，熵值变化很小时，可近似认为是一种具有间断面的简单压缩波，其传播过程是等熵的。,弱冲击波的这种近似处理方法称为冲击波的声学近似。,122,2.7,弱冲击波的声学近似理论,下面讨论弱冲击波阵面前后参数间的关系。,由冲击波根本关系式知：,(1),123,2.7,弱冲击波的声学近似理论,将1式中的 作为p的函数在p0附近按台劳级数展开，得到,(2),保存一阶时，视为等熵那么2式中的偏微熵写成全微分。,124,2.7,弱冲击波的声学近似理论,即 (3),(4),把3式代入4式可得：,(5),125,2.7,弱冲击波的声学近似理论,由,两边取对数，,两边微分,126,2.7,弱冲击波的声学近似理论,因此,即 (6),把5式代入6式得：,即 (7),127,2.7,弱冲击波的声学近似理论,由冲击波速度：,又,因此,128,2.7,弱冲击波的声学近似理论,将上式按二项式展开并忽略高阶项，可得：,忽略高阶后，得,129,2.7,弱冲击波的声学近似理论,把5式代入上式得：,所以 (8),即弱冲击波的传播速度是波前后小扰动速度的平均值。,130,2.7,弱冲击波的声学近似理论,因此弱冲击波前后参数间的关系为：,131,本章要点,掌握弱扰动、强扰动、压缩波、稀疏波、冲击波的概念；,了解声速的推导过程，掌握声速表达式；,理解特征线的物理意义，简单波区参数的计算；,掌握平面正冲击波根本关系式，理解波速线、Hugoniot曲线的含义；,掌握多方气体中强冲击波关系表达式；,理解冲击波正反射的规律；,了解弱冲击波的声学近似理论。,