• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    水下爆炸條件下自由場壓力載荷時(shí)頻特征分析

    2016-10-13 01:28:40謝耀國崔洪斌李新飛姚熊亮
    中國艦船研究 2016年2期
    關(guān)鍵詞:自由場沖擊波小波

    謝耀國,崔洪斌,李新飛,姚熊亮

    哈爾濱工程大學(xué)船舶工程學(xué)院,黑龍江哈爾濱150001

    水下爆炸條件下自由場壓力載荷時(shí)頻特征分析

    謝耀國,崔洪斌,李新飛,姚熊亮

    哈爾濱工程大學(xué)船舶工程學(xué)院,黑龍江哈爾濱150001

    在艦船結(jié)構(gòu)水下爆炸試驗(yàn)中,為了研究水下爆炸條件下水中自由場壓力載荷的時(shí)頻特征,針對某水下爆炸試驗(yàn)自由場壓力測試試驗(yàn)數(shù)據(jù),基于小波分析對信號進(jìn)行時(shí)頻特性分析,得到水中自由場壓力信號的時(shí)頻分布和能量分布狀況。分析結(jié)果表明:針對水下爆炸自由場壓力載荷,基于小波分析技術(shù)對其時(shí)頻特征進(jìn)行分析,可得到水下爆炸自由場壓力載荷所包含的頻率信息、強(qiáng)度信息以及不同頻段下的載荷持續(xù)作用時(shí)間等信息;另外,可對沖擊波、滯后流和二次壓力波這3個(gè)不同信號階段進(jìn)行頻段與能量統(tǒng)計(jì)分析;在不同頻段上對沖南擊信號的能量進(jìn)行統(tǒng)計(jì)發(fā)現(xiàn),沖擊波階段在8 kHz以下頻段集中了超過90%的能量,其中4 kHz以下頻段的能量最大,在滯后流和二次壓力波階段,需特別重視250 Hz以下的低頻段對船體結(jié)構(gòu)及設(shè)備的影響,該結(jié)果對艦船結(jié)構(gòu)及設(shè)備的抗沖擊防護(hù)具有借鑒意義。

    水下爆炸;小波分析;時(shí)頻特征;自由場壓力

    網(wǎng)絡(luò)出版地址:http://www.cnki.net/kcms/detail/42.1755.TJ.20160317.1056.010.html期刊網(wǎng)址:www.ship-research.com

    引用格式:謝耀國,崔洪斌,李新飛,等.水下爆炸條件下自由場壓力載荷時(shí)頻特征分析[J].中國艦船研究,2016,11(2):27-32,50. XIE Yaoguo,CUI Hongbin,LI Xinfei,et al.Time-frequency characteristics analysis of free-field pressure with underwater explosion[J].Chinese Journal of Ship Research,2016,11(2):27-32,50.

    0 引 言

    在艦船結(jié)構(gòu)抗爆抗沖擊研究中,水中自由場壓力載荷是反映艦船在水下爆炸時(shí)所受外載荷的重要數(shù)據(jù),也是進(jìn)行抗爆抗沖擊數(shù)值計(jì)算的重要輸入?yún)?shù)[1-2]。為了提高艦船的抗爆抗沖擊能力,世界各國進(jìn)行了大量的實(shí)船及模型水下爆炸試驗(yàn)[3-5]。近年來,我國也開展了很多相關(guān)的試驗(yàn),獲得了一批壓力載荷試驗(yàn)數(shù)據(jù)[6-7],同時(shí),也有一些學(xué)者針對沖擊波和氣泡脈動壓力對艦船結(jié)構(gòu)及設(shè)備的影響進(jìn)行了分析[8]。

    在信號分析領(lǐng)域,小波分析是近幾十年提出并發(fā)展起來的一種新的分析方法,小波變換的出現(xiàn)克服了傳統(tǒng)傅立葉分析方法的不足,其窗口大小對信號具有適應(yīng)性,能提供非平穩(wěn)性號的時(shí)域和頻域的局部化信息,被稱為信號分析的“數(shù)學(xué)顯微鏡”。小波分析在一些信號時(shí)頻處理上已得到成功應(yīng)用,并且在分析水下爆炸信號以及工程爆破信號等非平穩(wěn)、非線性信號上,已取得了一些成果[9-11]。水下爆炸自由場壓力信號具有非平穩(wěn)、瞬態(tài)、短時(shí)等特點(diǎn),信號的頻譜結(jié)構(gòu)是時(shí)刻變化的。本文擬針對水下爆炸自由場壓力載荷信號的特點(diǎn),以實(shí)測的水下爆炸壓力信號為對象,利用小波分析多層多分辨率的優(yōu)勢,對實(shí)測壓力信號能量分布特征進(jìn)行分析,以期獲得水下爆炸自由場壓力信號的時(shí)頻特征信息。

    1 基本原理

    1.1小波分析

    設(shè)任意信號 f(t)∈L2(R),其中,L2(R)為能量有限的信號空間,則稱 f(t)為能量有限的信號[12-13],即

    如果ψ(t)∈L2(R),則其傅立葉變換ψ∧(ω)滿足容許性條件

    即Cψ有界,則稱ψ為一個(gè)基小波或母小波。將基小波經(jīng)伸縮和平移后,便可得到一個(gè)小波序列

    式中,a,b∈R,且a≠0,其中a為伸縮因子,b為平移因子。定義下式:

    式(4)為關(guān)于基小波的連續(xù)小波變換。其中,的共軛運(yùn)算??梢姡〔ㄗ儞Q可以將一維信號變換為二維信號。從上式中可以看到,變換后的函數(shù)是二維的,可為分析信號的時(shí)頻特性奠定基礎(chǔ)。對實(shí)際計(jì)算來講,通常將連續(xù)小波變換進(jìn)行離散化,把參數(shù)a或b,或是同時(shí)將a和b進(jìn)行離散化,從而得到離散小波變換。通常取

    式中,m為信號的離散采樣點(diǎn)。每當(dāng)m增加1時(shí),伸縮因子a會增加一倍,對應(yīng)的信號頻率則減小一半,可見適合進(jìn)行高效計(jì)算,也便于分析。

    1.2小波變換后各頻段能量統(tǒng)計(jì)

    原始響應(yīng)信號經(jīng)小波分解后,可得到由高到低各分析頻段內(nèi)的小波分量,各小波分量仍然為時(shí)間歷程信號。如果將信號分析至第n層,則各層對應(yīng)的能量為

    式中:Ei為第i頻段信號對應(yīng)的能量和;fi(t)為第i頻段的小波分解信號;yij為第i頻段信號fi(t)的離散點(diǎn)幅值;E為分析信號的總能量;i=1,2,···,n+1;j=1,2,···,m。

    各小波分解頻段能量在分析信號總能量中的占比為

    2 水下爆炸自由場壓力試驗(yàn)

    水下爆炸自由場壓力試驗(yàn)在某大型水下爆炸試驗(yàn)水池進(jìn)行,水池直徑50 m,水深20 m。進(jìn)行了一系列的水下爆炸試驗(yàn),由于分析方法的相似性,下面僅以其中的一種工況進(jìn)行分析。試驗(yàn)工況為:TNT藥包質(zhì)量0.5 kg,藥包水深10 m,測點(diǎn)與藥包水深相同,測點(diǎn)與藥包間的距離為5 m。圖1所示為試驗(yàn)測得的水中沖擊壓力時(shí)間歷程曲線。

    根據(jù)水下爆炸理論,將水下爆炸過程分為3個(gè)階段:沖擊波階段、因氣泡膨脹速度小于水中聲速推動水質(zhì)點(diǎn)做徑向運(yùn)動而形成的滯后流(擴(kuò)散流)壓力階段,以及由于氣泡脈動所形成的二次壓力波階段。試驗(yàn)清晰地測得了各階段過程。一般認(rèn)為,只有一次氣泡脈動產(chǎn)生的壓力波有實(shí)際意義,因此,只分析到了二次壓力波,之后的氣泡脈動未做分析。根據(jù)水下爆炸相關(guān)研究[1],壓力波各階段已經(jīng)有很好的估算公式:

    式中:pm為沖擊波峰值壓力;θ為指數(shù)衰減時(shí)間常數(shù);T為氣泡脈動周期;W為藥包質(zhì)量;R為與藥包間的距離;H為藥包水深。

    由圖1可見,試驗(yàn)測得的沖擊波峰值壓力為6.41 MPa,利用式(10)計(jì)算得到的峰值壓力為6.66 MPa,其誤差為3.8%;試驗(yàn)測得的氣泡脈動周期為135.25 ms,利用式(12)計(jì)算得到的脈動周期為136.26 ms,誤差為0.74%。試驗(yàn)結(jié)果與理論計(jì)算結(jié)果相近,表明試驗(yàn)結(jié)果可靠。

    圖1 水中自由場壓力測試結(jié)果Fig.1 Free-field pressure curves of underwater explosion experiment

    3 自由場壓力載荷時(shí)頻特征分析

    3.1信號小波分解

    在實(shí)際應(yīng)用小波分析的過程中,小波基的選取是非常重要的問題。利用不同的小波基對同一信號進(jìn)行分析會得到不同的結(jié)果[14]。在非穩(wěn)定信號分析中,Daubechies小波基函數(shù)是最常用的,Daubechies小波系列具有較好的緊支撐性、光滑性和近似對稱性[15]。在Daubechies小波系列函數(shù)中,db8小波基函數(shù)在非平穩(wěn)定信號分析中應(yīng)用較多,本文也將利用db8小波基函數(shù)對實(shí)測信號進(jìn)行分解。

    試驗(yàn)設(shè)置的信號采樣頻率為1 MHz,信號分析的奈奎斯特(Nyquist)頻率根據(jù)采樣定理為0.5 MHz。以實(shí)測的水下爆炸自由場壓力信號為對象,根據(jù)小波分析原理,利用db8小波基函數(shù)對試驗(yàn)信號進(jìn)行小波變換。本文將對實(shí)測信號進(jìn)行11層小波分解,分解頻段如表1所示,分解后共獲得12個(gè)小波分解系數(shù)。

    表1 小波分解頻段表Tab.1 Corresponding frequency band of wavelet decomposition

    進(jìn)行自由場壓力信號小波分解后,為確定分解后的信號是否還能真實(shí)地反映實(shí)測的壓力信號,將分解后的各頻段信號進(jìn)行完全重構(gòu),并與實(shí)測信號進(jìn)行對比,計(jì)算相對誤差。圖2所示為重構(gòu)信號與原始信號的相對誤差分布圖。從對比結(jié)果可以看出,重構(gòu)信號與實(shí)測信號的相對誤差在10-11量級以內(nèi),相差極小。這說明采用選取的小波基來分解水下爆炸壓力信號是合適的,在信號分解過程中能量的損失可以忽略不計(jì),能真實(shí)地反映測試信號的情況。基于db8小波基的小波重構(gòu)分層信號如圖3所示,分別對應(yīng)于表1中的12個(gè)頻段。圖中,a11為低頻分量,d1~d11為高頻分量。

    圖2 相對誤差分布圖Fig.2 Distribution of relative error

    圖3 信號小波分解后的分層重構(gòu)信號Fig.3 Reconstructed signals of free-field pressure in undergoing wavelet decomposition at different levels

    3.2各頻段能量分布統(tǒng)計(jì)

    為了分析水下爆炸自由場壓力信號在分析時(shí)間段內(nèi)各頻段能量的分布情況,利用式(6)~式(8),根據(jù)小波變換分層重構(gòu)信號,可以得到測試信號不同頻段(fi)上信號相對能量的分布情況、各頻段信號能量在總能量中的占比ki,以及各頻段信號壓力幅值Pimax,統(tǒng)計(jì)如表2所示。

    表2 小波分解信號各頻段信息統(tǒng)計(jì)Tab.2 Band parameters for free-field pressure

    為進(jìn)一步分析水下爆炸自由場壓力載荷在不同階段、不同頻段上的能量分布情況,按照水下爆炸理論,將所分析的自由場壓力波時(shí)間歷程分為3個(gè)階段:沖擊波階段、滯后流階段和二次壓力波階段。定義各階段峰值壓力時(shí)間的中點(diǎn)為分界線,即以沖擊波壓力和滯后流壓力峰值時(shí)間的中點(diǎn)為界劃分沖擊波與滯后流,以沖擊波和二次壓力波峰值時(shí)間的中點(diǎn)(氣泡脈動半徑最大點(diǎn))為界劃分滯后流與二次壓力波。即0~7.917 ms為沖擊波作用時(shí)間,7.918~71.067 ms為滯后流作用時(shí)間,71.068 ms以后為二次壓力波作用時(shí)間。根據(jù)小波變換分層重構(gòu)信號,對沖擊波、滯后流和二次壓力波各階段、各頻段的能量進(jìn)行統(tǒng)計(jì),每階段能量在所統(tǒng)計(jì)總能量中的占比如圖4所示,各階段內(nèi)各頻段的能量占比如圖5所示。

    圖4 各階段能量占比Fig.4 Energy ratio for each stage

    圖5 各階段內(nèi)各頻段的能量占比Fig.5 Energy ratio for each spectrum in each stage

    3.3結(jié)果分析

    1)水下爆炸自由場壓力載荷信號經(jīng)過小波分解,利用小波分析的多分析多分辨率的特性,將實(shí)測的壓力信號分解為12個(gè)頻段的時(shí)間—壓力信號細(xì)節(jié)信息。通過對各頻段信號局部信息的能量統(tǒng)計(jì),得到自由場壓力信號的能量分布特征,各頻段下所包含的沖擊能量大小可以用來表征沖擊信號在各自頻段下的強(qiáng)度信息。通過表2可以發(fā)現(xiàn),低頻段a11(0~244.14 Hz)壓力信號的能量在各頻段中的占比最大,達(dá)到了39.19%。頻段d7~a11(0~7 812.5 Hz)壓力信號的能量占比極大,占總能量的91.22%,可見自由場壓力載荷信號能量主要集中在8 kHz以下,特別是4 kHz以下。在d1和d2這2個(gè)高頻段(125~500 kHz),不但其壓力信號的能量占比極小,而且壓力幅值也很小,再結(jié)合d1和d2的壓力時(shí)程曲線可以發(fā)現(xiàn),其信噪比很低,可以認(rèn)為這2個(gè)頻段上的能量主要來源于噪聲??梢?,水下爆炸自由場壓力載荷信號在經(jīng)過小波變換后,各頻段能量信息可以反映沖擊壓力在相應(yīng)頻段上的沖擊強(qiáng)度。

    2)通過圖3可以發(fā)現(xiàn),自由場壓力信號包含豐富的頻率成分,頻段極寬,從高頻到低頻均有壓力波動,再結(jié)合能量分布情況,就能知道自由場壓力載荷的主要作用頻段。所以,水下爆炸自由場壓力載荷信號的有效頻率信息可以通過信號小波變換后的各頻段能量大小的統(tǒng)計(jì)信息來表征。

    3)通過對表2的分析可見,各小波頻段能量峰值與沖擊壓力峰值并不在同一頻段上,再結(jié)合圖3可以發(fā)現(xiàn),如果壓力信號在某小波頻段上衰減慢、持續(xù)時(shí)間長,則其所包含的能量也會相對較大。進(jìn)一步分析發(fā)現(xiàn),信號分解后所有頻段上壓力幅值最大的是d5(15 625~31 250 Hz),與a11(0~244.14 Hz)相比能量要小得多,通過對比d5和a11的小波頻段壓力時(shí)程曲線即可看出,雖然d5頻段的幅值最大,但其沖擊壓力信號衰減極快,持續(xù)作用時(shí)間短,因而其相對能量較小。因此,對水下爆炸自由場壓力載荷信號進(jìn)行小波分析后,其各頻段上的小波分層重構(gòu)信號可以反映所包含的壓力載荷在相應(yīng)頻段下的衰減及持續(xù)時(shí)間信息。

    4)通過對壓力波時(shí)歷曲線分階段能量的統(tǒng)計(jì),由圖4可見,沖擊波、滯后流和二次壓力波這3個(gè)階段均包含有大量的能量,其中沖擊波階段為能量最大段,集中了65%的能量。結(jié)合圖3可以發(fā)現(xiàn),沖擊波階段是沖擊強(qiáng)度最大、頻段分布最寬的階段;滯后流階段包含的頻段與沖擊波階段相似,能量相對較低;二次壓力波階段包含的頻段較窄,主要在低頻段,能量頻段相對集中。所以,基于小波分析和能量統(tǒng)計(jì)方法,可以得到水下爆炸自由場壓力載荷信號沖擊波、滯后流以及二次壓力波各階段在不同頻段下的能量分布信息。

    5)通過對沖擊波、滯后流和二次壓力波這3個(gè)階段內(nèi)各頻段在其階段內(nèi)的能量占比的統(tǒng)計(jì)(圖5)可以發(fā)現(xiàn),沖擊波階段的能量頻段分布較廣,但主要能量集中在4 kHz以下頻段(占比84.46%);滯后流能量的分布與沖擊波階段具有相似性,但有很大部分能量集中在低頻段,在244.14 Hz以下,有50%的能量集中在這里;二次壓力波階段則完全不同,其能量要相對集中,有很強(qiáng)的低頻特性,在244.14 Hz以下,其能量占比可達(dá)78%。因此,通過小波分析可見,沖擊波的能量頻段分布較廣,氣泡脈動壓力的能量則以低頻為主。水下爆炸壓力載荷沖擊波瞬時(shí)沖擊壓力峰值大,能量分布頻率范圍廣,在艦船遭受水下爆炸沖擊載荷作用時(shí),是造成艦船結(jié)構(gòu)局部損傷的主要原因,同時(shí),具有大量能量的低頻段易與艦船的整體振動模態(tài)耦合,是船體總強(qiáng)度的主要威脅,并且低頻段與很多艦船設(shè)備的安裝頻率重合,也易于同設(shè)備的隔離系統(tǒng)產(chǎn)生共振現(xiàn)象,從而對設(shè)備造成損傷。

    4 結(jié) 論

    綜上所述,可以得出以下結(jié)論:

    1)小波分析方法相較于傳統(tǒng)的傅立葉分析方法,具有更好的時(shí)頻分析能力,可以滿足水下爆炸自由場壓力載荷信號的時(shí)頻特征分析要求,可以得到分析信號在不同頻段下的時(shí)歷曲線,以及強(qiáng)度、頻率和持續(xù)時(shí)間等時(shí)頻細(xì)節(jié)信息。

    2)基于小波分析和能量統(tǒng)計(jì)方法,可以得到自由場壓力載荷信號沖擊波、滯后流以及二次壓力波等階段在不同頻段上的能量分布信息。結(jié)果表明,沖擊波階段在8 kHz以下頻段集中了超過90%的能量,其中4 kHz以下頻段所包含的能量最大,滯后流和二次壓力波階段在250 Hz以下低頻段包含有很大的能量,會對船體結(jié)構(gòu)及設(shè)備產(chǎn)生影響,應(yīng)特別予以注意。

    3)小波分析方法對研究水下爆炸沖擊載荷以及艦船結(jié)構(gòu)破壞機(jī)理提供了一種新途徑,相關(guān)分析結(jié)果對艦船結(jié)構(gòu)與設(shè)備的抗爆抗沖擊設(shè)計(jì)具有實(shí)際的參考價(jià)值。

    [1]COLE R H.Underwater explosion[M].Princeton,New Jersey:Princeton University Press,1948.

    [2]汪玉,華宏星.艦船現(xiàn)代沖擊理論及應(yīng)用[M].北京:科學(xué)出版社,2005.

    [3] MAIR H U,REESE R M,HARTSOUGH K.Simulated shipshocktests/trials?[EB/OL].[2009-06-09]. http://www.dote.osd.mil/lfte/SSS.HTM.

    [4]PUSEY H C.A survey of the information needs of in?dustry in designing to meet navy shipboard require?ments[J].Shock and Vibration Bulletin,1969,40:132-139.

    [5]MILLER R D,MOYER E T.Shock analysis of ship structures due to underwater explosions using HULL?CAV3/DYNA3D with DAA[C]//Proceedings of the 62th Shock and Vibration Symposium.Rome,1991: 449-457.

    [6] 金輝,張慶明,張姝紅,等.水下爆炸中壓力載荷測量與分析[J].爆破,2009,26(2):18-22. JIN Hui,ZHANG Qingming,ZHANG Shuhong,et al. Measurement and analysis of underwater explosion pressure load[J].Blasting,2009,26(2):18-22.

    [7]樊寶順,程素秋,韓峰.艦船艙段模型在水下爆炸作用下的壁壓分析[J].中國艦船研究,2009,4(5):20-22. FAN Baoshun,CHENG Suqiu,HAN Feng.Wall pres?sure analysis of cabin model subjected to underwater explosion[J].Chinese Journal of Ship Research,2009,4(5):20-22.

    [8] 李國華,李玉節(jié),張效慈,等.艦船設(shè)備沖擊環(huán)境的能源研究[J].船舶力學(xué),2001,2(1):37-54. LI Guohua,LI Yujie,ZHANG Xiaoci,et al.Study of power for shock environment of ship equipment[J]. Journal of Ship Mechanics,2001,2(1):37-54

    [9]溫華兵,張健,尹群,等.水下爆炸船艙沖擊響應(yīng)時(shí)頻特征的小波包分析[J].工程力學(xué),2008,25(6):199-203. WEN Huabing,ZHANG Jian,YIN Qun,et al.Wavelet packet analysis of time-frequency characteristic of cab?in shock response due to underwater explosion[J].En?gineering Mechanics,2008,25(6):199-203.

    [10]李夕兵,張義平,劉志祥,等.爆破震動信號的小波分析與HHT變換[J].爆炸與沖擊,2005,25(6):528-535. LI Xibing,ZHANG Yiping,LIU Zhixiang,et al.Wave?let analysts and Hilbert-Huang Transform of blasting vibration signal[J].Explosion and Shock Waves, 2005,25(6):528-535.

    [11]謝耀國,郭君,李新飛,等.水下爆炸船體結(jié)構(gòu)壁壓信號時(shí)頻特征分析[J].大連理工大學(xué)學(xué)報(bào),2015,55(5):492-497. XIE Yaoguo,GUO Jun,LI Xinfei,et al.Analysis of time-frequency characteristics of wall pressure sig?nals of hull structure subjected to under water explo?sion[J].Journal of Dalian University of Technology,2015,55(5):492-497.

    [12] 凌同華,廖艷程,張勝.沖擊荷載下巖石聲發(fā)射信號能量特征的小波包分析[J].振動與沖擊,2010,29 (10):127-130,255. LING Tonghua,LIAO Yancheng,ZHANG Sheng.Ap?plication of wavelet packet method in frequency band energy distribution of rock acoustic emission signals under impact loading[J].Journal of Vibration and Shock,2010,29(10):127-130,255

    [13]張德豐.Matlab小波分析[M].2版.北京:機(jī)械工業(yè)出版社,2012.

    [14]DAUBECHIES I.The wavelet transform,time-fre?quency localization and signal analysis[J].IEEE Transactions on Information Theory,1990,36(5):961-1005.

    [15]DAUBECHIES I.Orthonormal bases of compactly supported wavelets[J].Communications on Pure and Applied Mathematics,1988,41(7):909-996.

    Time-frequency characteristics analysis of free-field pressure with underwater explosion

    XIE Yaoguo,CUI Hongbin,LI Xinfei,YAO Xiongliang
    School of Shipbuilding Engineering,Harbin Engineering University,Harbin 150001,China

    To obtain the characteristics of the free-field pressure with underwater explosion,the time-fre?quency characteristics of the monitored free-field pressure signals are studied through wavelet analysis,based on the experimental free-field pressure data subjected to underwater explosion.By using these sig?nals,the pressure-time curves and the energy distributions in different blasting frequency bands are ob?tained.The result shows that the obtained time-frequency characteristics of free-field pressure can be used to generate the shock wave signal,the after flow signal,and the second pressure pulse,and by analyzing the energy and frequency component,it is easy to get the detailed time-frequency information of the free-field pressure signal intensity,frequency and duration.More than 90%of the energy of the shock wave pressure signals condenses in the band lower than 8 kHz,and the highest level of energy appears below the frequen?cy of 4 kHz,and special attention should be directed to the low frequency band below 250 Hz of the after flow signal as well as the second pressure pulse for ship structure and equipment.Overall,this paper pro?vides reference for shock-resistance analysis of the ship structure.

    underwater explosion;wavelet analysis;time-frequency characteristics;free-field pressure

    U661.44

    A

    10.3969/j.issn.1673-3185.2016.02.005

    2015-07-08網(wǎng)絡(luò)出版時(shí)間:2016-3-17 10:56

    國家自然科學(xué)基金資助項(xiàng)目(51279038);國家級重大基礎(chǔ)研究項(xiàng)目

    謝耀國(通信作者),男,1982年生,博士,講師。研究方向:船舶結(jié)構(gòu)動力學(xué),船舶結(jié)構(gòu)實(shí)驗(yàn)技術(shù)。E-mail:xieyaoguo@hrbeu.edu.cn

    姚熊亮,男,1963年生,博士,教授。研究方向:艦船結(jié)構(gòu)動力學(xué),水下爆炸,氣泡動力學(xué)。

    E-mail:xiongliangyao@hrbeu.edu.cn

    猜你喜歡
    自由場沖擊波小波
    構(gòu)造Daubechies小波的一些注記
    基于MATLAB的小波降噪研究
    電子制作(2019年13期)2020-01-14 03:15:32
    武漢沖擊波
    中國公路(2019年10期)2019-06-28 03:05:08
    能源物聯(lián)網(wǎng)沖擊波
    能源(2018年10期)2018-12-08 08:02:34
    基于改進(jìn)的G-SVS LMS 與冗余提升小波的滾動軸承故障診斷
    微活動:孕育童心習(xí)作的自由場
    醫(yī)生集團(tuán)沖擊波
    三維層狀黏彈性半空間中球面SH、P和SV波源自由場
    考慮地震波幅值衰減的斜入射二維自由場
    振動與沖擊(2016年1期)2016-07-26 02:21:36
    考慮反演及樁土相互作用的擬動力試驗(yàn)方法
    欧美精品人与动牲交sv欧美| 狂野欧美白嫩少妇大欣赏| av黄色大香蕉| .国产精品久久| 观看av在线不卡| .国产精品久久| 国产黄频视频在线观看| 久久久久人妻精品一区果冻| 久久99热6这里只有精品| 一级爰片在线观看| 亚洲欧美一区二区三区黑人 | 黄片无遮挡物在线观看| 丰满乱子伦码专区| 三级国产精品片| 在线观看免费视频网站a站| 黄色日韩在线| 免费观看av网站的网址| 精品熟女少妇av免费看| 亚洲最大成人中文| 欧美+日韩+精品| 成人亚洲精品一区在线观看 | 久久精品人妻少妇| 日韩人妻高清精品专区| 永久网站在线| 伦理电影大哥的女人| 久久久久国产网址| 午夜视频国产福利| 亚洲性久久影院| 人人妻人人爽人人添夜夜欢视频 | 国产乱人视频| 久久久精品94久久精品| 国产亚洲一区二区精品| 熟女av电影| 中文字幕av成人在线电影| 国产亚洲欧美精品永久| 麻豆成人午夜福利视频| 国产精品人妻久久久久久| 欧美精品人与动牲交sv欧美| 国产精品精品国产色婷婷| 97精品久久久久久久久久精品| 男女无遮挡免费网站观看| 寂寞人妻少妇视频99o| 各种免费的搞黄视频| 精品久久久精品久久久| 精品一区二区三卡| 亚洲av中文av极速乱| 亚洲内射少妇av| 亚洲国产精品成人久久小说| av线在线观看网站| 狠狠精品人妻久久久久久综合| 日本vs欧美在线观看视频 | 一区二区av电影网| 国产视频首页在线观看| a 毛片基地| 国产久久久一区二区三区| 国产伦理片在线播放av一区| 亚洲欧美一区二区三区黑人 | 男女边摸边吃奶| 亚洲怡红院男人天堂| 夫妻午夜视频| 美女福利国产在线 | 国产伦在线观看视频一区| 欧美区成人在线视频| 少妇被粗大猛烈的视频| 久久久久人妻精品一区果冻| 亚洲欧美中文字幕日韩二区| 丝瓜视频免费看黄片| 国产在线视频一区二区| 亚洲精品久久久久久婷婷小说| 精品久久国产蜜桃| 美女cb高潮喷水在线观看| 蜜桃亚洲精品一区二区三区| 97超碰精品成人国产| 欧美精品一区二区大全| 久久久久久伊人网av| 国产av国产精品国产| 狂野欧美白嫩少妇大欣赏| 成人综合一区亚洲| 成人亚洲欧美一区二区av| 五月伊人婷婷丁香| 国产视频首页在线观看| 在线免费十八禁| 少妇精品久久久久久久| 免费av不卡在线播放| 街头女战士在线观看网站| 日本猛色少妇xxxxx猛交久久| 国产成人a∨麻豆精品| 51国产日韩欧美| 精品亚洲乱码少妇综合久久| 成人高潮视频无遮挡免费网站| 久久热精品热| 欧美精品国产亚洲| 一区二区三区乱码不卡18| 国产91av在线免费观看| 久久久久久久久久成人| 久久久久久九九精品二区国产| 国产精品人妻久久久影院| 欧美xxxx黑人xx丫x性爽| 香蕉精品网在线| 美女中出高潮动态图| 国产在线免费精品| 亚洲美女黄色视频免费看| 亚洲精品国产色婷婷电影| 精品一区二区三区视频在线| 亚洲不卡免费看| 日本欧美视频一区| 亚洲中文av在线| 视频区图区小说| 一本久久精品| 亚洲精华国产精华液的使用体验| 视频区图区小说| 精品久久久精品久久久| 国产成人午夜福利电影在线观看| 在线观看一区二区三区| 最近手机中文字幕大全| 国产精品人妻久久久久久| 国产真实伦视频高清在线观看| 久久 成人 亚洲| 亚洲图色成人| 日本免费在线观看一区| 高清av免费在线| 日韩一区二区三区影片| 夫妻性生交免费视频一级片| 亚洲国产欧美人成| 久久久久久久国产电影| 91久久精品电影网| 久久久久人妻精品一区果冻| videos熟女内射| 国产女主播在线喷水免费视频网站| 欧美少妇被猛烈插入视频| 国产精品一二三区在线看| 免费看光身美女| 日韩制服骚丝袜av| 亚洲欧美日韩另类电影网站 | 一级毛片 在线播放| 美女cb高潮喷水在线观看| 久热这里只有精品99| 亚洲精品,欧美精品| 91久久精品国产一区二区成人| 免费黄频网站在线观看国产| 街头女战士在线观看网站| 卡戴珊不雅视频在线播放| av在线老鸭窝| 2022亚洲国产成人精品| 美女福利国产在线 | 亚洲伊人久久精品综合| 亚洲人与动物交配视频| 夫妻性生交免费视频一级片| 成人18禁高潮啪啪吃奶动态图 | 欧美成人一区二区免费高清观看| 欧美97在线视频| 亚洲av中文av极速乱| 亚洲,一卡二卡三卡| 极品少妇高潮喷水抽搐| 日本与韩国留学比较| 91精品伊人久久大香线蕉| 免费少妇av软件| freevideosex欧美| 18禁裸乳无遮挡免费网站照片| 日韩一区二区视频免费看| 91aial.com中文字幕在线观看| 大片免费播放器 马上看| 男人添女人高潮全过程视频| 日韩精品有码人妻一区| 亚洲精品第二区| 久久鲁丝午夜福利片| 亚洲综合色惰| 国产69精品久久久久777片| 网址你懂的国产日韩在线| 亚洲av福利一区| a级毛色黄片| 少妇高潮的动态图| 在线观看三级黄色| 国产无遮挡羞羞视频在线观看| 久久国产精品大桥未久av | 一级黄片播放器| 国产在线免费精品| 噜噜噜噜噜久久久久久91| 纯流量卡能插随身wifi吗| 寂寞人妻少妇视频99o| 99九九线精品视频在线观看视频| 亚洲第一av免费看| 26uuu在线亚洲综合色| 国内精品宾馆在线| 少妇高潮的动态图| 搡老乐熟女国产| 中文字幕免费在线视频6| 精品熟女少妇av免费看| 一级毛片久久久久久久久女| 国产色婷婷99| 亚洲欧洲国产日韩| 黄色视频在线播放观看不卡| 秋霞伦理黄片| 十八禁网站网址无遮挡 | 新久久久久国产一级毛片| 欧美三级亚洲精品| 国产亚洲欧美精品永久| 麻豆成人午夜福利视频| 另类亚洲欧美激情| 日日啪夜夜撸| 偷拍熟女少妇极品色| 国产视频内射| 久久久欧美国产精品| 日本色播在线视频| 欧美日韩国产mv在线观看视频 | 亚洲av福利一区| 乱系列少妇在线播放| 欧美人与善性xxx| 搡女人真爽免费视频火全软件| 国模一区二区三区四区视频| 少妇精品久久久久久久| 亚洲精品乱码久久久久久按摩| 一级二级三级毛片免费看| 亚洲精品456在线播放app| 伦精品一区二区三区| 久久99精品国语久久久| 偷拍熟女少妇极品色| 蜜桃久久精品国产亚洲av| 99精国产麻豆久久婷婷| 国产一区亚洲一区在线观看| 最近手机中文字幕大全| 中文字幕人妻熟人妻熟丝袜美| av天堂中文字幕网| 综合色丁香网| 国产午夜精品一二区理论片| 99热网站在线观看| 亚洲国产成人一精品久久久| 2022亚洲国产成人精品| 国产人妻一区二区三区在| 精品酒店卫生间| 激情五月婷婷亚洲| 夜夜爽夜夜爽视频| 亚洲伊人久久精品综合| 纯流量卡能插随身wifi吗| 在线观看av片永久免费下载| 卡戴珊不雅视频在线播放| 赤兔流量卡办理| 七月丁香在线播放| 欧美成人午夜免费资源| av播播在线观看一区| av不卡在线播放| 不卡视频在线观看欧美| 伦理电影免费视频| 国产淫语在线视频| 亚洲精品乱久久久久久| 亚洲人成网站在线播| 亚洲av国产av综合av卡| 中文字幕精品免费在线观看视频 | 日韩免费高清中文字幕av| 国产成人一区二区在线| 国产精品欧美亚洲77777| av线在线观看网站| 日韩伦理黄色片| 韩国av在线不卡| 国产成人精品久久久久久| 一二三四中文在线观看免费高清| 久久久久久久久久成人| 国产精品欧美亚洲77777| 国产日韩欧美亚洲二区| 欧美bdsm另类| 国产一级毛片在线| 免费高清在线观看视频在线观看| 日本黄大片高清| 精品久久国产蜜桃| 免费少妇av软件| 搡女人真爽免费视频火全软件| 青春草国产在线视频| 亚洲欧美一区二区三区黑人 | av在线app专区| 国产 一区精品| 国产精品久久久久久精品电影小说 | 国产精品秋霞免费鲁丝片| 建设人人有责人人尽责人人享有的 | 一本—道久久a久久精品蜜桃钙片| 91精品国产九色| 亚洲一区二区三区欧美精品| 国产人妻一区二区三区在| 人体艺术视频欧美日本| 嫩草影院入口| 女人久久www免费人成看片| 国产精品精品国产色婷婷| 男女无遮挡免费网站观看| 美女国产视频在线观看| 丝瓜视频免费看黄片| 欧美日韩在线观看h| 中文乱码字字幕精品一区二区三区| 看十八女毛片水多多多| 久久99热这里只频精品6学生| 小蜜桃在线观看免费完整版高清| 国产国拍精品亚洲av在线观看| 久久精品国产a三级三级三级| 免费高清在线观看视频在线观看| 久久人人爽人人片av| 日本黄大片高清| 国模一区二区三区四区视频| 老司机影院毛片| 亚洲,一卡二卡三卡| 欧美日本视频| 男女国产视频网站| 小蜜桃在线观看免费完整版高清| 男女边吃奶边做爰视频| 99久久精品一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 国产色爽女视频免费观看| 99久久精品热视频| 99久久人妻综合| 久久国产乱子免费精品| 婷婷色av中文字幕| 日本av免费视频播放| 久久精品夜色国产| 亚洲精品日韩在线中文字幕| 亚洲精品成人av观看孕妇| 国产高潮美女av| 午夜福利影视在线免费观看| 一级毛片黄色毛片免费观看视频| 春色校园在线视频观看| .国产精品久久| 啦啦啦视频在线资源免费观看| 亚洲性久久影院| av一本久久久久| 丝袜喷水一区| 国产精品偷伦视频观看了| 成人18禁高潮啪啪吃奶动态图 | 美女内射精品一级片tv| 国产视频内射| 国产黄片美女视频| 久久久亚洲精品成人影院| 欧美一区二区亚洲| 久热久热在线精品观看| 亚洲av二区三区四区| 欧美精品国产亚洲| 日本黄色片子视频| 成人亚洲精品一区在线观看 | 国国产精品蜜臀av免费| 国产精品一区二区在线不卡| 麻豆乱淫一区二区| 人妻制服诱惑在线中文字幕| 99热网站在线观看| 在线免费十八禁| 欧美人与善性xxx| 日韩精品有码人妻一区| 国产在视频线精品| 久久久精品免费免费高清| 哪个播放器可以免费观看大片| 日韩欧美精品免费久久| h视频一区二区三区| 这个男人来自地球电影免费观看 | 午夜免费观看性视频| 国产 一区精品| 精品久久久久久久末码| 国产欧美日韩精品一区二区| 婷婷色综合大香蕉| 简卡轻食公司| 精品酒店卫生间| 男女边吃奶边做爰视频| 日日啪夜夜爽| 日日啪夜夜撸| 亚洲欧美清纯卡通| 日韩一区二区视频免费看| 亚洲真实伦在线观看| 一级毛片久久久久久久久女| 国产日韩欧美亚洲二区| 91午夜精品亚洲一区二区三区| 视频中文字幕在线观看| 极品少妇高潮喷水抽搐| 欧美精品人与动牲交sv欧美| 天美传媒精品一区二区| 日日啪夜夜爽| 精品熟女少妇av免费看| 久久久精品94久久精品| 久久久久久人妻| av网站免费在线观看视频| 亚洲婷婷狠狠爱综合网| 国产成人午夜福利电影在线观看| 国产成人免费观看mmmm| 99热6这里只有精品| 人人妻人人澡人人爽人人夜夜| a级毛片免费高清观看在线播放| 欧美成人午夜免费资源| 交换朋友夫妻互换小说| av一本久久久久| 亚洲av在线观看美女高潮| 我的老师免费观看完整版| 美女xxoo啪啪120秒动态图| 国产深夜福利视频在线观看| av在线蜜桃| 国产黄片视频在线免费观看| 免费黄频网站在线观看国产| 亚洲精品久久久久久婷婷小说| 一二三四中文在线观看免费高清| 久久婷婷青草| 精品国产乱码久久久久久小说| 丝袜脚勾引网站| 欧美日韩综合久久久久久| 国产精品爽爽va在线观看网站| 精品亚洲成国产av| 午夜激情久久久久久久| 小蜜桃在线观看免费完整版高清| 欧美日韩视频高清一区二区三区二| 中文资源天堂在线| 一级二级三级毛片免费看| 欧美人与善性xxx| 国精品久久久久久国模美| 网址你懂的国产日韩在线| av国产久精品久网站免费入址| 国产男女超爽视频在线观看| 色婷婷久久久亚洲欧美| 午夜福利在线在线| 久久久久久久大尺度免费视频| 亚洲欧美成人综合另类久久久| 伊人久久国产一区二区| 又粗又硬又长又爽又黄的视频| 国产精品99久久99久久久不卡 | 精品国产一区二区三区久久久樱花 | 亚洲美女搞黄在线观看| 18禁裸乳无遮挡动漫免费视频| 大又大粗又爽又黄少妇毛片口| 亚洲精品色激情综合| 免费人成在线观看视频色| 亚洲精品国产成人久久av| 亚洲国产精品999| av国产精品久久久久影院| 久久久久久久久久久免费av| 2021少妇久久久久久久久久久| av网站免费在线观看视频| 少妇熟女欧美另类| 欧美日韩精品成人综合77777| 国产深夜福利视频在线观看| 777米奇影视久久| 人妻 亚洲 视频| 亚洲精品国产成人久久av| 中国美白少妇内射xxxbb| 色视频在线一区二区三区| 黄色配什么色好看| 一区二区三区免费毛片| 国产免费福利视频在线观看| 熟女人妻精品中文字幕| 深夜a级毛片| 亚洲av不卡在线观看| 国产伦精品一区二区三区视频9| 国产无遮挡羞羞视频在线观看| freevideosex欧美| 在线亚洲精品国产二区图片欧美 | 精品久久久久久久久亚洲| 亚洲av在线观看美女高潮| 国产精品人妻久久久影院| 日韩欧美一区视频在线观看 | 国产一级毛片在线| 亚洲av男天堂| 亚洲av中文字字幕乱码综合| 久久精品久久精品一区二区三区| 高清午夜精品一区二区三区| 最近的中文字幕免费完整| 欧美日韩视频精品一区| 亚洲av在线观看美女高潮| 在线观看三级黄色| 国产成人一区二区在线| 欧美zozozo另类| 在线精品无人区一区二区三 | 又粗又硬又长又爽又黄的视频| 国产 精品1| 色网站视频免费| 久久6这里有精品| 十八禁网站网址无遮挡 | 国国产精品蜜臀av免费| 一级爰片在线观看| 国产男人的电影天堂91| 在线观看人妻少妇| 少妇人妻久久综合中文| 亚洲欧美一区二区三区国产| 你懂的网址亚洲精品在线观看| 日本av手机在线免费观看| 成人美女网站在线观看视频| 久久精品国产自在天天线| 亚洲av综合色区一区| 只有这里有精品99| 国产人妻一区二区三区在| 欧美日韩一区二区视频在线观看视频在线| 午夜免费鲁丝| 另类亚洲欧美激情| 啦啦啦在线观看免费高清www| 伊人久久国产一区二区| 男女下面进入的视频免费午夜| 热re99久久精品国产66热6| 国产乱人视频| 精品人妻熟女av久视频| 亚洲av欧美aⅴ国产| 五月天丁香电影| 寂寞人妻少妇视频99o| 99久久精品国产国产毛片| 国产熟女欧美一区二区| 欧美日韩国产mv在线观看视频 | 亚洲高清免费不卡视频| 国产欧美亚洲国产| av国产精品久久久久影院| 免费大片黄手机在线观看| 日本av手机在线免费观看| 国产成人午夜福利电影在线观看| a级毛片免费高清观看在线播放| 成人午夜精彩视频在线观看| 欧美国产精品一级二级三级 | 内地一区二区视频在线| 国产精品福利在线免费观看| 免费黄色在线免费观看| 国产乱来视频区| 亚洲国产高清在线一区二区三| 51国产日韩欧美| 少妇被粗大猛烈的视频| 联通29元200g的流量卡| 91在线精品国自产拍蜜月| av.在线天堂| 韩国高清视频一区二区三区| 七月丁香在线播放| 久久亚洲国产成人精品v| 欧美成人一区二区免费高清观看| 午夜日本视频在线| 高清毛片免费看| 夫妻午夜视频| 黄色日韩在线| 国产中年淑女户外野战色| 国产伦精品一区二区三区视频9| 国产伦精品一区二区三区四那| 久久 成人 亚洲| 国产女主播在线喷水免费视频网站| 蜜桃亚洲精品一区二区三区| 亚洲综合色惰| 91精品一卡2卡3卡4卡| 能在线免费看毛片的网站| 亚洲精品亚洲一区二区| 内地一区二区视频在线| 欧美国产精品一级二级三级 | 亚洲国产精品专区欧美| 成人一区二区视频在线观看| 免费观看无遮挡的男女| 国产永久视频网站| 熟妇人妻不卡中文字幕| 欧美国产精品一级二级三级 | 国产精品秋霞免费鲁丝片| av视频免费观看在线观看| 亚洲高清免费不卡视频| 日日啪夜夜撸| 久久久a久久爽久久v久久| 久久人人爽av亚洲精品天堂 | 久久99热6这里只有精品| 日韩大片免费观看网站| 内射极品少妇av片p| 在线 av 中文字幕| 最新中文字幕久久久久| 国产午夜精品一二区理论片| 亚洲国产最新在线播放| 国产精品福利在线免费观看| 在线观看一区二区三区激情| 校园人妻丝袜中文字幕| 交换朋友夫妻互换小说| 国产深夜福利视频在线观看| 搡老乐熟女国产| 夜夜看夜夜爽夜夜摸| 国产亚洲一区二区精品| 国产精品人妻久久久影院| 久久鲁丝午夜福利片| 精品人妻熟女av久视频| www.色视频.com| 国产成人一区二区在线| 在线 av 中文字幕| 日韩不卡一区二区三区视频在线| 国产又色又爽无遮挡免| 中文字幕精品免费在线观看视频 | 久久久久久久亚洲中文字幕| 91在线精品国自产拍蜜月| 久久久久久九九精品二区国产| av线在线观看网站| 美女xxoo啪啪120秒动态图| 日本黄大片高清| 午夜精品国产一区二区电影| 亚洲精品国产成人久久av| 国产有黄有色有爽视频| 免费看日本二区| 亚洲欧洲日产国产| 美女脱内裤让男人舔精品视频| 亚洲成人一二三区av| 久久久国产一区二区| 网址你懂的国产日韩在线| 三级国产精品片| 亚洲av福利一区| 国产黄片美女视频| 性色av一级| 亚洲内射少妇av| 大陆偷拍与自拍| 美女视频免费永久观看网站| 亚洲国产欧美人成| 九九在线视频观看精品| 香蕉精品网在线| 一级毛片电影观看| 免费观看在线日韩| 99热这里只有精品一区| 青春草亚洲视频在线观看| 成年美女黄网站色视频大全免费 | 亚洲欧美日韩无卡精品| h视频一区二区三区| 成人无遮挡网站| 亚洲天堂av无毛| 欧美变态另类bdsm刘玥| 中文字幕精品免费在线观看视频 | 久久婷婷青草| 国模一区二区三区四区视频| 日本色播在线视频| 另类亚洲欧美激情| 亚洲av免费高清在线观看| 亚洲av成人精品一二三区| 人人妻人人爽人人添夜夜欢视频 | 中文字幕制服av| 久久久久人妻精品一区果冻| 寂寞人妻少妇视频99o| 久久精品国产亚洲av天美| 观看美女的网站|