沖擊波基本理論P(yáng)PT課件

1、2021-10-271 波：在彈性介質(zhì)中，某個(gè)局部受到作用后，由于物質(zhì)點(diǎn)的相互作用，由近及遠(yuǎn)地使物質(zhì)質(zhì)點(diǎn)陸續(xù)發(fā)生擾動(dòng)，這種擾動(dòng)在介質(zhì)的傳播就稱為波。常見(jiàn)的如：水波，音波，電磁波 波陣面：介質(zhì)的原始狀態(tài)與擾動(dòng)狀態(tài)的交界面稱波陣面 縱波與橫波： 波陣面移動(dòng)方向與介質(zhì)質(zhì)點(diǎn)振動(dòng)方向平行的波稱縱波。 波陣面移動(dòng)方向與介質(zhì)質(zhì)點(diǎn)振動(dòng)方向垂直的波稱橫波。 波速：波陣面在介質(zhì)中傳播的速度。 波的傳播方向：波陣面的移動(dòng)方向。2.1一維等熵流動(dòng)2.1.1波的基本概念(復(fù)習(xí)）第1頁(yè)/共83頁(yè)2021-10-272 壓縮波：波陣面到達(dá)之處，介質(zhì)的狀態(tài)（P P、T T）參數(shù)增加的波稱壓縮波，波的傳播方向與介質(zhì)運(yùn)動(dòng)方向相同

2、。（圖5.15.1） 膨脹波（稀疏波）：波陣面到達(dá)之處，介質(zhì)的狀態(tài)（P P、T T）參數(shù)減小的波稱膨脹波，波的傳播方向與介質(zhì)運(yùn)動(dòng)方向相反。 （下圖5.25.2） 第2頁(yè)/共83頁(yè)2021-10-273 音波：介質(zhì)質(zhì)點(diǎn)在原來(lái)的位置振動(dòng)，而波向左右傳播，這種波稱音波，音波是弱壓縮波或膨脹波的合成。 沖擊波：是波面以突躍面的形式在彈性介質(zhì)中傳播的壓縮波，波陣面上介質(zhì)的狀態(tài)參數(shù)變化是突躍的。 爆轟波：是含有化學(xué)反應(yīng)能量支持的沖擊波，因?yàn)橛谢瘜W(xué)反應(yīng)能量的支持，因此爆轟波所以具有穩(wěn)定的傳播特性。第3頁(yè)/共83頁(yè)2021-10-274完全氣體，量熱完全氣體與等熵關(guān)系 （補(bǔ)物理化學(xué)知識(shí)）理想氣體(完全氣體pe

3、rfect gas)：不考慮分子間的作用力和分子的體積情況下，一種理想化后的氣體。它滿足： PV=nRT， e=e(T)和Cv=Cv(T) 世上無(wú)理想氣體，熱完全氣體是真實(shí)氣體在一定溫度，壓力范圍內(nèi)的近似，即近似看成理想氣體來(lái)處理。對(duì)于熱完全氣體，有：de=CvdT=Cv(T)dT ，dh=CpdT=Cp(T)dT，e=e(T) ，h=h(T)可近似認(rèn)為一定溫度范圍內(nèi)，Cv，Cp ， （ Cp- Cv =R）保持不變。 但一般說(shuō)來(lái)， Cv=Cv(T) ， Cp=Cp(T) VVPCRCC1第4頁(yè)/共83頁(yè)2021-10-275VPCC第5頁(yè)/共83頁(yè)6 ：分子平動(dòng)和轉(zhuǎn)動(dòng)的總自由度（不包括振動(dòng)）

4、(因?yàn)?， )所以：對(duì)單原子分子氣體： ， ，對(duì)雙原子分子氣體： ， ，對(duì)三原子分子氣體： ， ， 為多方指數(shù)或絕熱指數(shù)adiabatic exponent）自自由度由度解釋：決定一個(gè)物體位置所需要的獨(dú)立坐標(biāo)數(shù)，這里指的是熱力學(xué)自由度亦稱準(zhǔn)自由度，不同于一般的力學(xué)自由度。 RfCV2f11()22etr RTfRTfRdTdeCV21ffCRV213fRCV2367. 15fRCV254 . 16fRCV333. 1第6頁(yè)/共83頁(yè)2021-10-277等熵關(guān)系的建立：等熵關(guān)系的建立：一般地： （1）對(duì)可逆過(guò)程： （2）比較（1）和（2）有： （3）),(VTSS ),(VTee dVVSdT

5、TSdSTV)()(dVVedTTedeTV)()(PdVdVVSTdTTSTPdVTdSdeTV)()(VVTSTTe)()() ()TSSTdTTP dVTV第7頁(yè)/共83頁(yè)2021-10-278對(duì)焓、 Helmholtz自由能、 Gibbs自由焓的表達(dá)式分別微分： （4） （5） （6）而： ， ， ， （7） VdPTdSVdPPdVdedhSdTPdVTSddedf)(SdTVdPTSddhdg)(),(VSee ),(PShh ),(TVff ),(TPgg dVVedSSedeSV)()(dPPhdSShdhSP)()(dTTfdVVfdfVT)()(dTTgdPPgdgPT)(

6、)(hePVfeTSghTS第8頁(yè)/共83頁(yè)2021-10-279將（2）的第一式、（4）、（5）、（6）與（7）的4個(gè)式子比較有： （8）又因?yàn)椋?（ ） 所以： PVShSeT)()(TSVfVeP)()(TSPgPhV)()(PVTgTfS)()(dVVedSSedeSV)()(PdVTdSdeVVSSSVSPSVeVTVSe)()()()(),(VSee 第9頁(yè)/共83頁(yè)2021-10-2710即： 類似有： （9）（Maxwell關(guān)系）將（9）的第二式代入（1）的第一式有： （1）的第一式 又由（3）式： ，代入上式：有： （10）若 ， ， （11） VSSPVT)()(TVVST

7、P)()(PSSVPT)()(TPPSTV)()(dVTPdTTSdSVV)()(TCTTeTSVVV)()(dVTPTdTCdSVV)(),(PTSS ),(PThh dPPSdTTSdSTP)()(dPPhdTThdhTP)()(dVVSdTTSdSTV)()(第10頁(yè)/共83頁(yè)2021-10-2711而 類似有： 代入（11）的第1式： （12）（10），（12）就是熵函數(shù)的一般表達(dá)式（微分形式），也可以寫成積分形式： （13）dPVPSTdTTSTVdPTdSdhTP)()(PPPCThTST)()()()PPTPCCSVdSdTdPdTdVTPTT00()()VVPPdTPSSCdV

8、TTdTVSSCdPTT第11頁(yè)/共83頁(yè)2021-10-2712理想氣體： （14） （15）定義： 絕熱指數(shù)又因?yàn)椋?， ，代入（15）式：)(TCCPP)(TCCVVRTPV constVRTdTCSconstVRTdTCSTTPTTVlnln00constPRTCSconstVRTCSPVlnlnlnlnvpCC1RCV1RCP第12頁(yè)/共83頁(yè)2021-10-2713 對(duì)絕熱可逆過(guò)程（必等熵）： ，所以有：又因?yàn)椋?，所以： 或 或 多方氣體的等熵關(guān)系，亦為絕熱關(guān)系。 constPTRSconstVTRS)ln()ln(1110dSconstS 111TVc o n stTc o n

9、 stPRTPV constPVconstPconstTV1（*）（*）（*）第13頁(yè)/共83頁(yè)2021-10-2714定容比熱，定壓比熱以及兩者之間的關(guān)系定容比熱，定壓比熱以及兩者之間的關(guān)系比熱的定義： ，質(zhì)量比熱單位為： 由熱力學(xué)第一定律： （16）熱焓定義： （17）對(duì)定容過(guò)程，由（16）得：對(duì)定壓過(guò)程，由（17）得： （18） qCd T)/(KkgJVVdTqC)(PPdTqC)(PdVqdeVdPqVdPPdVdedhPVehVVTeC)(PPPPTVPTeThC)()()(第14頁(yè)/共83頁(yè)2021-10-2715因?yàn)椋?，所以： （19）即： （20）由（18）（20）有： （

10、21）與（1），（2）式 比較 ，有：(22) ),(TVee ),(TPVV PTVPTVVeTeTe)()()()(PTPVTVVeTeC)()()(PTVPTVVePCC)()(PVSTVeTT)()(dVVedTTedeTV)()() ()VTSSTdTTP dVTVPdVdVVSTdTTSTPdVTdSdeTV)()(（2）dVVSdTTSdSTV)()(（1）（22）第15頁(yè)/共83頁(yè)2021-10-2716又由Maxwell關(guān)系： （23）故有： （24）對(duì)理想氣體： 故： ， 代入（24）式： (25)由定義（比熱比）： 故: VTTPVS)()(PVVPTVTPTCC)()(

11、RTPV VRTPV)(PRTVP)(RPVRTCCVP2VPCC1RCV第16頁(yè)/共83頁(yè)流場(chǎng)流場(chǎng)：流體運(yùn)動(dòng)所占據(jù)的空間，流場(chǎng)中任一質(zhì)點(diǎn)流體的物理量如 等是空間的位置（ ）（或 ）和時(shí)間t的函數(shù)： 或 , 或 等。如果流場(chǎng)中的物理量只是位置函數(shù)，而與時(shí)間無(wú)關(guān)，則稱為定常流場(chǎng)，這種流動(dòng)就稱為定常流動(dòng)（定常流動(dòng)（steady flow）,否則為不定常否則為不定常(unsteady flow)的。如果流場(chǎng)中各物理量在空間分布只與一個(gè)幾何坐標(biāo)x有關(guān)，那么就稱為一維一維(one dimensional)流場(chǎng)流場(chǎng)，相應(yīng)的流動(dòng)稱為一維流動(dòng)一維流動(dòng)(one dimensional flow)。推導(dǎo)條件推導(dǎo)條

12、件：忽略氣體的粘性，熱傳導(dǎo)（絕熱），無(wú)化學(xué)變化，不考慮體積力（如重力（對(duì)氣體可忽略），電磁力）對(duì)流動(dòng)的影響，只有體積膨脹功。 ,TPzyx,),(tzyxPP2.1.2流場(chǎng)和定常流動(dòng)方程組( , )pp r t( , )TT r t( , , , )TT x y z t第17頁(yè)/共83頁(yè)連續(xù)性方程的推導(dǎo)（質(zhì)量守恒方程）：取如下圖所示的控制體（開(kāi)口系，當(dāng)?shù)赜^點(diǎn)即Euler方法），變截面流管。變截面流管中x1處的截面積為A，密度為 ，氣體流速為u單位時(shí)間內(nèi)流入控制體的質(zhì)量為：同樣時(shí)間內(nèi)從x2面流出的質(zhì)量為：微元dx中氣體質(zhì)量的變化率為：由質(zhì)量守恒，單位時(shí)間內(nèi)流入微元體x的質(zhì)量流出x的質(zhì)量微元體x的

13、質(zhì)量對(duì)時(shí)間的變化率。 xx1x2AuAu+AuxxAuAu)(xxAu()A xt第18頁(yè)/共83頁(yè)即：即： （控制體體積不變， 與t無(wú)關(guān)） （1） （ ， ， ） 連續(xù)方程（當(dāng)?shù)赜^點(diǎn)）物質(zhì)導(dǎo)數(shù)（Lagrange導(dǎo)數(shù)）的變換關(guān)系： 稱為Euler導(dǎo)數(shù)。物理量的物質(zhì)導(dǎo)數(shù)（或稱隨體導(dǎo)數(shù)）是指某個(gè)封閉系統(tǒng)中的流體在運(yùn)動(dòng)過(guò)程中，它所具有的物理量F（如： ）對(duì)時(shí)間的變化率, 是物理量F隨流體質(zhì)點(diǎn)運(yùn)動(dòng)時(shí)的變化率。xA ,txAxxAuAuAu)()(xtAtxAxxAu)()(0)(xAutA),(tx),(txuu )(xAA tF ,TPV第19頁(yè)/共83頁(yè)tFdtdFt0lim物質(zhì)導(dǎo)數(shù)的定義：以求加

14、速度為例，給出物質(zhì)導(dǎo)數(shù)的微分變換關(guān)系：設(shè)流體質(zhì)點(diǎn)在流場(chǎng)中沿運(yùn)動(dòng)軌跡C運(yùn)動(dòng)，從當(dāng)?shù)赜^點(diǎn)出發(fā)，流體速度為：假定t時(shí)刻，流體微團(tuán)在M點(diǎn)，速度為 ，經(jīng)時(shí)刻 后，運(yùn)動(dòng)到N點(diǎn)，速度為：加速度： （2）由于流場(chǎng)的非均勻性和不定常性，該微團(tuán)的速度在運(yùn)動(dòng)過(guò)程中不止經(jīng)歷了 的變化，而且也經(jīng)歷了 的變化。當(dāng)然 也與時(shí)間 長(zhǎng)短有關(guān)。 tttMN),(tMu),(ttNu( , )( , , , )uu r tu x y z t(, )u M t(,)uN ttr(,)(, )limtoduu N ttu M twdtt rxlyjzk r第20頁(yè)/共83頁(yè)（2）式可寫為： （3） 代表沿S方向移動(dòng)單位長(zhǎng)度引起的速度變

15、化，而如今單位時(shí)間移動(dòng)了u的距離，所以S方向的速度變化為 。對(duì)一維情況有： （4） 00( ,)( , )( , )(, )limlimttduu N ttu N tu N tu M tdttt 00(,)(, )(, )(, )limlimlimttNMNMu N ttu N tNMu N tu M tttNM (, )(, )u M tu M tutS(, )u M tSuuSduuuudttx第21頁(yè)/共83頁(yè)對(duì)于 等，亦有同樣的變化關(guān)系： （5）這里， ：全導(dǎo)數(shù)，物質(zhì)導(dǎo)數(shù)，隨體導(dǎo)數(shù)，Lagrange導(dǎo)數(shù)。 ：當(dāng)?shù)貢r(shí)間導(dǎo)數(shù)，局部導(dǎo)數(shù)，Euler導(dǎo)數(shù)。反映了流場(chǎng)的不定常性，反映了流體微團(tuán)流

16、過(guò)空間固定點(diǎn)上量F對(duì)時(shí)間的變化率。 ：遷移導(dǎo)數(shù)，對(duì)流導(dǎo)數(shù)，反映了流場(chǎng)的非均勻性，是流體微團(tuán)運(yùn)動(dòng)到不同位置時(shí)所引起的F的變化。實(shí)際上，F(xiàn)=F(x,y,z,t), 而 x=x(t), y=(t), z=z(t)所以：,P三維（直角坐標(biāo)系）dFFFudttxFtd Fd tFuxdFFF dxF dyF dzdttx dty dtz dt第22頁(yè)/共83頁(yè)由（5）式，可將（1）式化為： （6） 隨體觀點(diǎn)的連續(xù)方程0)(xAuAdtd注：()()()00(6)dAudAuAuAAuAudtxxdtxxx第23頁(yè)/共83頁(yè)歐拉方程動(dòng)量守恒方程（運(yùn)動(dòng)方程）的推導(dǎo)：取下圖的控制體（閉口系，隨體觀點(diǎn)，即Lag

17、range方法 ），設(shè)微元體dx的側(cè)面積為S，該質(zhì)點(diǎn)具有的速度為u，為管壁切線與x軸的夾角（如果管壁是光滑的，則是無(wú)窮小量）顯然： ，即： x微元體x1面受到壓力為PA，x2面受到的壓力為：側(cè)面所受力為： ，即：)(dAAxxPPsinSdAsinSdASxxPPP2SxxPP)21(xPnPn1x2xPA()()PPXAdAX第24頁(yè)/共83頁(yè)該力在x方向投影為：在與x垂直方向投影為： （互相抵消） 微元體受到的總壓力為（不考慮粘性力，重力等）：忽略二階小量，總壓力為：按Newton第二定律： dAxPPSxxPP)21(sin)21(cos)21(SxxPPcos)21(cos)21()2

18、1()(SxxPPSxxPPdAxxPPdAAxxPPPAdAxxPPdAAxxPPPA)21()(xxPAdtduxAxxPA( F m a ) x第25頁(yè)/共83頁(yè)即： （7）或： （8） 歐拉方程（動(dòng)量守恒方程） 由開(kāi)口體系（Euler觀點(diǎn)推導(dǎo)動(dòng)量方程）：由x1面流入dx的動(dòng)量：由x2面流出dx的動(dòng)量： （忽略二階以上小量） 01xPdtdu01xPxuutuuAu)()()(dxxuudxxAAdxxuudxxuAuuAu)()()()()Auudx Adx udx udxxxxx()()()AuuAudxudxxx第26頁(yè)/共83頁(yè)微元體dx受的合外力為：?jiǎn)挝粫r(shí)間內(nèi)，微元體動(dòng)量變化為

19、： （忽略二階小量）凈流入的動(dòng)量：流入流出 dxxuAu)(dxx1x2AuAu+()A ud xxdxxPA)21()21)(21(dxxuudxdxxAAdxxt)21()(21(dxxuudxdxxAAtdxdxxAuAut)(21()(dxAutx第27頁(yè)/共83頁(yè)動(dòng)量定理：動(dòng)量的增加率凈增加動(dòng)量+微元體受的外力，即： （ 與t無(wú)關(guān)） 即：dxxuAudxxPAdxAut)()(xuAuxPAutA)()(dxA,xPAxuAuxAuutuA)()(xPAxuAuxAuutAutuA)(xPAxAutAudtduA)(xPdtdu)0)(xAutA（1）式，質(zhì)量守恒方程第28頁(yè)/共83

20、頁(yè)能量方程的推導(dǎo)（忽略熱損失，不考慮非體積力做功，只計(jì)體積功;開(kāi)口系，Euler方法）：?jiǎn)挝毁|(zhì)量氣體總能量為： （ e ：?jiǎn)挝毁|(zhì)量?jī)?nèi)能， ：?jiǎn)挝毁|(zhì)量動(dòng)能）單位時(shí)間內(nèi)通過(guò)x1面進(jìn)入微元體x的能量：?jiǎn)挝粫r(shí)間內(nèi)通過(guò)x2面流出微元體x的能量：x1面上，外力單位時(shí)間內(nèi)對(duì)微元所做的功為： （功率）221ue 221u)2(2ueAuPAux2面上微元體單位時(shí)間內(nèi)克服外力所做的功為：微元體x的總能量變化率為：xxPAuPAu)(txAue)21(222()2()2uA ueuA uexx第29頁(yè)/共83頁(yè)由能量守恒：微元體x的總能量變化率應(yīng)等于單位時(shí)間內(nèi)流進(jìn)的凈能量加上外力做功的和。 與時(shí)間無(wú)關(guān)，（控制體體

21、積不變） （9）)21()21()2()21(2222xxAuueAuueAuuetxAuexAuAu+xxAux)(xxPP)(dxxPAuPAuPAuxA xPAuxAuueAtue)()21()21(22P2x1x第30頁(yè)/共83頁(yè)因?yàn)椋?( (1)式連續(xù)方程)故上式可簡(jiǎn)化為： 或： （10） 能量守恒方程再由（8）式（動(dòng)量方程） （11）0)(xAutA0)()21(2xPAudtuedA0)(xAuPxPAudtduAudtdeA0)()(xAuPxPdtduAudtdeA0)(xAuAPdtde22()()()220uueePAuAuAutxx第31頁(yè)/共83頁(yè)由熱力學(xué)第一定律：

22、：環(huán)境給封閉系統(tǒng)傳遞的熱量， ：系統(tǒng)內(nèi)能的增加， ：系統(tǒng)對(duì)外界（環(huán)境）所做的功。微分形式為：若只考慮體積功，則有： （E：內(nèi)能；Q：熱量；P：壓力）或： （q：?jiǎn)挝毁|(zhì)量的供熱量；e：?jiǎn)挝毁|(zhì)量?jī)?nèi)能；P：壓強(qiáng)）又因?yàn)椋?（ 對(duì)封閉體系的可逆過(guò)程） ， ， WEQQEWWQdEPdVQdEPdVdeqTdSQ 1VddV21dPdeTdS2（12）（同除以dt )2dSdep dTdtdtdt第32頁(yè)/共83頁(yè)將（11）式、（6）式代入（12）式有：即： 或 （13）可見(jiàn)，某封閉體系（流體微團(tuán)）絕熱可逆條件下的流動(dòng)是等熵的（對(duì)于無(wú)粘流體）（完全氣體的絕熱流動(dòng)必為等熵流動(dòng)）：因?yàn)椋核裕?)()(/2

23、xAuAPxAuAPdtTdS0dtdS0 xSutS),(PSS constPSS),(由以上推導(dǎo)的非定常流動(dòng)的基本方程組為（變截面，一維）。0)(xAuAdtd(6)( ,)SS Pconst第33頁(yè)/共83頁(yè) 連續(xù)方程（質(zhì)量守恒方程） 運(yùn)動(dòng)方程（動(dòng)量守恒方程） 能量方程 方程組（） 或 狀態(tài)方程守恒方程是普遍適用的，對(duì)任何流體都相同。狀態(tài)方程則反映了流體在流動(dòng)中的特殊性。四個(gè)方程，四個(gè)未知數(shù)： ，方程組封閉，可求解。對(duì)等熵過(guò)程（完全氣體的絕熱過(guò)程），方程組中的第三式（能量方程）（完全氣體的絕熱流動(dòng)必為等熵流動(dòng)）可用 或 來(lái)代替。對(duì)多方氣體，可用 代替 0)(xuAtA01xPxuutu0

24、)(xAuAPdtde),(TPP),(ePPeuP,0dS0 xSutSconstP第34頁(yè)/共83頁(yè)對(duì)定常流動(dòng)，所有物理量（ 等）對(duì)時(shí)間的偏導(dǎo)數(shù)為零。同時(shí)用熱焓 代替e，可得定常流動(dòng)方程組： 或 或 或 （等截面流管）eu,PePVeh0)(AudconstAu ududP0)2(2uhdconstuh22),(TPP01dPudu或方程組（） constuP2質(zhì)量動(dòng)量能量狀態(tài)第35頁(yè)/共83頁(yè)2021-10-2736 音波是弱的壓縮波或膨脹波合成的結(jié)果，波的傳播速度僅取決于介質(zhì)狀態(tài)。 音速可以看作是介質(zhì)的狀態(tài)參數(shù)，音速是弱擾動(dòng)在介質(zhì)中的傳播速度。2.1.3聲音的傳播第36頁(yè)/共83頁(yè)20

25、21-10-2737()()c dcud ()()ccu （1） 質(zhì)量守恒定律:流入波陣面的質(zhì)量等于流出波陣面的質(zhì)量音速公式推導(dǎo)（先以壓縮波）P C, P+dPC-u,c（1）第37頁(yè)/共83頁(yè)2021-10-2738 動(dòng)量守恒定律（動(dòng)量的變化=壓力變化時(shí)間）公式變形化簡(jiǎn)： （2） ()cdc()() ()()c udc uc dc u ()()c dcc dcuppp d ()cccu dpd c up 由動(dòng)量分析：流入的動(dòng)量流出的動(dòng)量（2）第38頁(yè)/共83頁(yè)2021-10-2739 由（1）得 ccuc()ccp 2()()(1)PPPc (1)Pc代入（2）所以第39頁(yè)/共83頁(yè)2021

26、-10-2740弱擾動(dòng) PdPdddPCd（3）同樣膨脹波也可導(dǎo)出（3）式弱的壓縮波的傳播速度只與壓力和介質(zhì)密度的比值有關(guān)。第40頁(yè)/共83頁(yè)2021-10-2741 如果弱擾動(dòng)的傳播與周圍介質(zhì)是絕熱的。壓縮和膨脹過(guò)程可以看成是等熵的（物化 ） 傳播表達(dá)式為： 如果音波的傳播介質(zhì)是理想氣體 由于 ，則 由以上可知： 則 （4）()sd PCd 211( )dddvvv2()()ssdpdpcvvdvdvd QST1/V2/ddV V第41頁(yè)/共83頁(yè)2021-10-2742在等熵過(guò)程中有： 代入（4）理想氣體kkapvapv11kkaapdpkdvkdvkdvvv vvdppkdvv pkpC

27、VKkpVVpvnRTCknRTk RT（5）音速與介質(zhì)的溫度、壓力有關(guān)，與介質(zhì)種類有關(guān)。K-絕熱指數(shù)pvkcc第42頁(yè)/共83頁(yè)2021-10-27432.2正沖擊波基本關(guān)系式*#活塞原理 圖2.1 沖擊波形成原理示意圖 第43頁(yè)/共83頁(yè)2.2.1沖擊波的基本性質(zhì) 壓縮密度速度壓縮迭加P、T壓縮追趕強(qiáng)壓縮波形成階躍變化（即沖擊波） 介質(zhì)狀態(tài)發(fā)生突躍變化的波即是沖擊波。也就是說(shuō)，沖擊波的波陣面是一個(gè)突躍面。第44頁(yè)/共83頁(yè)2021-10-2745 特點(diǎn) a. 波陣面的兩邊介質(zhì)狀態(tài)參數(shù)差別很大，是突變的或稱有較大梯度。 b. 波陣面運(yùn)動(dòng)的方向與介質(zhì)的運(yùn)動(dòng)方向相同。 c. 波陣面?zhèn)鞑ヒ詨嚎s波形

28、式傳播的。波后的 波前第45頁(yè)/共83頁(yè)2021-10-2746 沖擊波與音波的區(qū)別 傳播速度： （未擾動(dòng)介質(zhì)的音速） 狀態(tài)參數(shù)變化形式: 突躍-接近于零 介質(zhì)移動(dòng): 位移 平衡位置來(lái)回振動(dòng) 波速影響因素：強(qiáng)度有關(guān)僅與介質(zhì)狀態(tài) 參數(shù)有關(guān) 、T、P 周期性 ： 無(wú)周期 周期 熱力學(xué)特征 ： 熵增大 等熵0DC第46頁(yè)/共83頁(yè)2021-10-2747 波陣面兩側(cè)介質(zhì)狀態(tài)參數(shù)（T、 、P）和運(yùn)動(dòng)參數(shù)（ 、 ）之間的關(guān)系稱沖擊波基本關(guān)系式。 建立依據(jù)，三大定律: 質(zhì)量守恒 動(dòng)量守恒 能量守恒uD 第47頁(yè)/共83頁(yè)2021-10-27482.2.2 沖擊波的三大守恒定律 質(zhì)量守恒（物質(zhì)量相等，物質(zhì)不

29、滅） 介質(zhì)擾動(dòng)前的密度 介質(zhì)擾動(dòng)后的密度 沖波的速度 波陣面前介質(zhì)的移動(dòng)速度 波陣面后介質(zhì)的移動(dòng)速度 0011()()DDuu （2-1）o1Dvou1u第48頁(yè)/共83頁(yè)2021-10-2749 動(dòng)量守恒Fmu10PFP00()Dmu 0110()()DDuuuuu （2-2）100010()()Dppuuu 第49頁(yè)/共83頁(yè)2021-10-2750（能量變化等于對(duì)外所作的功） 能量=內(nèi)能+動(dòng)能 流入能量 流出 200001()() 2DDuEu 211111()() 2DDuEu 能量守恒 右側(cè) 做功 左側(cè)： 因作用力與運(yùn)動(dòng)方向相反，為負(fù)號(hào) 00()DPu11()DPu第50頁(yè)/共83頁(yè)

30、2021-10-2751化簡(jiǎn)整理 （2-3） 沖擊波基本關(guān)系式（2-1）、（2-2）、（2-3）221111000011()() ()() 22DDDDuEuuEu 0011()()DDpupu221 1001010001()2()Dpup uEEuuu 第51頁(yè)/共83頁(yè)2021-10-2752 由沖擊波的三個(gè)基本關(guān)系式可導(dǎo)出沖擊波的有關(guān)參數(shù)計(jì)算公式由狀態(tài)參數(shù)（ P、V、 ）計(jì)算沖擊波相關(guān)系數(shù)（ ）8個(gè)1111111, , , , ,DT c u p v E 2.2.3 沖擊波參數(shù)計(jì)算第52頁(yè)/共83頁(yè)2021-10-2753 由（5-11）式 解出 有0101DDuuVV0 11 010D

31、u vu vvv0 11 00101000101001()Du vu vu vvuuuvvvvvvv010001Duuuvvv則(2-4)Dv(2-5)第53頁(yè)/共83頁(yè)2021-10-2754把（5-16）代入（5-13） 10101001()uuppuuvv101001()()uuppvv100001Dppuvvv（5-17）（5-18）（2-6）（2-7）第54頁(yè)/共83頁(yè)2021-10-2755221 1001010001()()2Dpup uEEuuu 由（2-3）221 100101010()1()()2pup uEEuupp1 1001010102()1()()2pup uuuu

32、upp整理得2101010101()2ppEEuupp將（26）式代入得1010011()()2EEppvv（5-19）把（2-2）代入化簡(jiǎn)整理（2-8）第55頁(yè)/共83頁(yè)2021-10-2756 100001DppuVvv101001()()uuppvv1010011()()2EEppvv計(jì)算沖擊波參數(shù)計(jì)算公式?jīng)_擊波速度方程式（米海爾遜方程）瑞利線沖擊波絕熱方程（沖擊波雨果尼奧方程）第56頁(yè)/共83頁(yè)2021-10-27571 111p vEk001DpEk擾動(dòng)前后氣體介質(zhì)當(dāng)成理想氣體狀態(tài)下來(lái)處理：（2-9） 011010(1)(1)(1)(1)kvkvppkvkv111v101001(1)

33、(1)(1)(1)kpkpkpkp11 1ckp v1 1100 0p vTTp v代入（2-8）式，整理得：由氣體狀態(tài)初始參數(shù)（ ）可利用沖擊波關(guān)系式，計(jì)算沖擊波波陣面參數(shù)（ ）。0000000,T C up v E1111111, , ,T C u p v c第57頁(yè)/共83頁(yè)例1 已知一未擾動(dòng)空氣的初始參數(shù)為： =9.8 ， =1.25 ， =0，如果 波陣面的超壓 = Pa，試用沖擊波的關(guān)系式計(jì)算沖擊波的其他參數(shù)（假設(shè)氣體為熱容不變的理想氣體）解：由于空氣可以看成是雙原子氣體，因此我們可以取 =1.4。 將 ， ， ， ， 和 代入沖擊波有關(guān)的關(guān)系式中進(jìn)行計(jì)算，得：op410 Pao3

34、/kg mouPopooTocoup69.8 10第58頁(yè)/共83頁(yè)100010100(1)(1)(1)()(1)(1)(1)(1)(1)()opppppppppp 4644462.4 (9.8 109.8 10 )0.4 9.8 105.672.4 9.8 100.4 (9.8 109.8 10 )3105.675.67 1.257.09(/)kg m311110 .1 4 ()7 .0 9Vmk g300110 .8 (/)1 .2 5Vmk g1101010()()uuuppvv610()9.8 10(0.800.14)2543(/ )p vvm s第59頁(yè)/共83頁(yè)100000101D

35、DpppvvuVVVVVV69.8 100.803082(/ )(0.800.14)m s4611 11.4 (9.8 109.8 10 ) 0.14cpV1329(/ )m s461 110400(9.8 109.8 10 ) 0.14288(9.8 10 ) 0.80pVTTp V5059( )K第60頁(yè)/共83頁(yè)2021-10-2761 沖擊波陣面上已擾動(dòng)介質(zhì)的狀態(tài)參數(shù)主要與沖擊波波速有關(guān)； 沖擊波相對(duì)于未擾動(dòng)介質(zhì)是超音速的，即 ，相對(duì)于已擾動(dòng)的介質(zhì)是亞音速的，即 。 沖擊波速度大于介質(zhì)移動(dòng)速度且與介質(zhì)運(yùn)動(dòng)方向相同。 沖擊波衰減最終變?yōu)橐舨ā?沖擊波的沖擊壓縮過(guò)程是熵增大過(guò)程，波陣面介質(zhì)

36、的狀態(tài)參數(shù)變化是突躍的。 極強(qiáng)沖擊波陣面上，已擾動(dòng)介質(zhì)的密度取決于波陣面上的溫度。00DuC11Duc2.3沖擊波雨貢紐曲線及沖擊波的性質(zhì)第61頁(yè)/共83頁(yè)2021-10-2762證明：沖擊波相對(duì)于未擾動(dòng)介質(zhì)是超音速的，即 ，相對(duì)于已擾動(dòng)的介質(zhì)是亞音速，即 。00DuC11DuC舉例證：由沖擊波基本關(guān)系式可導(dǎo)出：100001111001()/()DDppuvvvuvppvv210100001000101()()()Dpppp vuvpp vvvvv00 0Ckp v10pp1001.2 1.5()kppkp100000()ppvkp vC00DuC，得證命題一，同樣可證命題二第62頁(yè)/共83頁(yè)

37、例2 沖擊波的傳播速度不僅與介質(zhì)的初始狀態(tài)有關(guān)，而且與沖擊波的強(qiáng)度有關(guān)。 證明：令 稱為激波強(qiáng)度 將其代入 有： 即： 故 式中 表示了介質(zhì)初始狀態(tài)的音速，對(duì)于聲波，有 ， ，因此聲波傳播速度只與介質(zhì)狀態(tài)有關(guān)。101001PPPPP2)1()1()(01200PPCuD2001()12DuC21100CuD),(0CfD 0C000CuD第63頁(yè)/共83頁(yè)2021-10-2764沖擊波參數(shù)計(jì)算導(dǎo)出下面兩式: 沖擊波速度方程式（米海爾直線） 沖擊波絕熱方程式（雨果尼奧曲線） 將過(guò)程和狀態(tài)表示在P-V圖上 （P縱坐標(biāo)，V橫坐標(biāo)）100001DppuVvv1010011()()2EEppvv沖擊波雨

38、貢紐曲線第64頁(yè)/共83頁(yè)2021-10-2765 當(dāng) 為一定數(shù)時(shí)， 為一直線 滿足: i) 過(guò)A（P0，V0） ii）斜率 當(dāng) 為一定數(shù)， 不同， 不同， ， 21002010()Dppuvvv00,Du v( )pf v2020()Dutgv00,u vDtgDtg波速線 a. 波速線為通過(guò)介質(zhì)初始狀態(tài)A（P0，V0）的不同斜率直線與沖擊波波速相對(duì)應(yīng)。 b. 波速線反映沖擊波以固定波速V0通過(guò)初始狀態(tài)（P0，V0）傳播時(shí)所有終起點(diǎn)的軌跡。 第65頁(yè)/共83頁(yè)Dv1Dv 沖擊波的波速線2Dv第66頁(yè)/共83頁(yè)2021-10-2767 介質(zhì)的內(nèi)能變化與波陣面P1和比容V1之間的關(guān)系。01110

39、1001010(1)(1)1()()2(1)(1)kvkvpEEppvvpkvkv 沖擊絕熱線 在P-V圖上，為一條以介質(zhì)初始狀態(tài)（P0，V0）凹向P軸和V軸的曲線，線上的每一點(diǎn)為不同波速?zèng)_擊波經(jīng)過(guò)同一初始狀態(tài)A（P0，V0）的介質(zhì)后所達(dá)到的終點(diǎn)狀態(tài)。第67頁(yè)/共83頁(yè)2021-10-2768*a. 介質(zhì)由A點(diǎn)（P0，V0）受沖擊壓縮到B點(diǎn)（P1，V1），其狀態(tài)不是沿著沖擊波絕熱曲線變化，而是突躍地從A點(diǎn)壓縮到B點(diǎn)。（沖擊波傳播時(shí)，終點(diǎn)既滿足波速線又滿足沖擊絕熱線）*b. 沖擊波絕熱曲線不是沖擊壓縮過(guò)程線，只是具有初始狀態(tài)，受到?jīng)_擊波壓縮時(shí)，一切可能到達(dá)的終態(tài)點(diǎn)B（P1，V1）狀態(tài)的集合。*c

40、. 由式10001DDppuvv10111DDppuvv10100101()ppuuvvv在沖擊絕熱線上A點(diǎn)以上所對(duì)應(yīng)的沖擊波， 壓縮波A點(diǎn)以下所對(duì)應(yīng)的沖擊波， 膨脹波10,0Du10,0Du A點(diǎn)，弱沖擊波，波速線與沖擊絕熱線相切音波0DC第68頁(yè)/共83頁(yè)2021-10-2769沖擊波是強(qiáng)縮波， 是波速線斜A點(diǎn)以后的斜率均大于切線，強(qiáng)沖擊波D0DC*d音波傳播是等熵的等熵線與沖擊絕熱線相切于A點(diǎn)*e沖擊波傳播過(guò)程是熵增大的沖擊波絕熱線在等熵線（過(guò)程線）上方 第69頁(yè)/共83頁(yè) 沿沖擊波絕熱線的熵值變化第70頁(yè)/共83頁(yè)例 2 ：絕熱線A-B 是熵增加過(guò)程證明：對(duì)激波壓縮過(guò)程，由熱力學(xué)第一定

41、律知： 或者 (A)如上圖所示，分析由Hugoniot曲線熵初態(tài)點(diǎn)A-B熵值的變化情況。由激波的Hugoniot方程： (B)兩邊取微分得： (C) 將式(C)代入式(A) ，得： (D) （式中 ）這就是沿Hugoniot曲線的熵表達(dá)式，實(shí)質(zhì)是Rayleigh向右掃過(guò)一點(diǎn)點(diǎn)，熵的增加量計(jì)算式 deTdSPdVTdSdePdV)(21000VVPPeedVPPdPVVde)(21)(2100001111()()2222TdSVV dPPP dVVdPPdV 0VVV第71頁(yè)/共83頁(yè) 圖中AB線、AC線與H線所圍成的陰影部分面積令為dF，則： 這就是所謂的“面積規(guī)則”，該面積等于（D）式的值。

42、dFABMACEBCNMNCE 0000011() ()()()()221()() ()21122dFVVdVP dP PVV P PdVdPdVP PVdPPdV 第72頁(yè)/共83頁(yè)a 、沖擊波的傳播過(guò)程自由傳播激波的自由傳播：指激波完全依靠自身的能量的傳播過(guò)程?；钊铀龠\(yùn)動(dòng)形成激波后，如果活塞突然停止運(yùn)動(dòng)，則激波失去外界能量補(bǔ)充，將依靠自身的能量繼續(xù)傳播?；钊蝗煌Ｖ购?，由于慣性，緊貼活塞的氣體質(zhì)點(diǎn)仍以活塞速度向前運(yùn)動(dòng)，這樣活塞前出現(xiàn)了空隙（稀疏），從而在受激波壓縮的氣體中產(chǎn)生膨脹波，傳播方向與激波方向一致，由于 并能追上激波，從而使激波強(qiáng)度減弱。此外，由于激波傳播過(guò)程中存在著粘性摩檫，熱傳導(dǎo)，熱輻射等不可逆能量損耗，也促使激波強(qiáng)度減弱。空中點(diǎn)爆炸：形成沖擊波為球形激波，其衰減速度比平面一維激波自由傳播時(shí)的衰減速度快得多。除膨脹波和不可逆能