• <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)方法
    国产探花极品一区二区| 国产亚洲午夜精品一区二区久久| 99久久99久久久精品蜜桃| 9191精品国产免费久久| 老熟女久久久| av有码第一页| 精品国产一区二区久久| 欧美日韩国产mv在线观看视频| 亚洲成人一二三区av| 飞空精品影院首页| tube8黄色片| 亚洲男人天堂网一区| 高清黄色对白视频在线免费看| 极品人妻少妇av视频| 国产高清国产精品国产三级| 免费在线观看视频国产中文字幕亚洲 | 国产成人精品无人区| 波多野结衣一区麻豆| 亚洲av男天堂| 一区二区三区四区激情视频| 亚洲一区中文字幕在线| 无遮挡黄片免费观看| 欧美黑人欧美精品刺激| 婷婷色麻豆天堂久久| 亚洲国产av新网站| 丰满少妇做爰视频| 18禁动态无遮挡网站| 高清在线视频一区二区三区| 国产精品.久久久| 日韩av不卡免费在线播放| 777米奇影视久久| 精品第一国产精品| 久久午夜综合久久蜜桃| 黄色毛片三级朝国网站| 国产亚洲一区二区精品| 国产亚洲av片在线观看秒播厂| 国产伦理片在线播放av一区| 亚洲国产欧美一区二区综合| 少妇的丰满在线观看| 伦理电影免费视频| 亚洲精品国产av成人精品| 性色av一级| 国产乱来视频区| 晚上一个人看的免费电影| 亚洲国产看品久久| 国产av精品麻豆| 国产极品天堂在线| 免费久久久久久久精品成人欧美视频| 女人爽到高潮嗷嗷叫在线视频| 黄网站色视频无遮挡免费观看| 免费少妇av软件| 日韩,欧美,国产一区二区三区| 精品国产乱码久久久久久男人| 在线观看免费日韩欧美大片| 国产精品 欧美亚洲| 色播在线永久视频| 国产成人欧美| 久久久欧美国产精品| 99精国产麻豆久久婷婷| 国产国语露脸激情在线看| 欧美另类一区| 亚洲精品日本国产第一区| 男人舔女人的私密视频| 精品少妇黑人巨大在线播放| 亚洲av欧美aⅴ国产| 免费观看人在逋| 午夜日本视频在线| 可以免费在线观看a视频的电影网站 | 亚洲综合精品二区| 亚洲精品在线美女| 欧美激情 高清一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 久久久精品国产亚洲av高清涩受| 天美传媒精品一区二区| 狂野欧美激情性bbbbbb| 91成人精品电影| 亚洲成av片中文字幕在线观看| 成人毛片60女人毛片免费| bbb黄色大片| 国产亚洲最大av| 日韩免费高清中文字幕av| 亚洲专区中文字幕在线 | 一区二区av电影网| 亚洲情色 制服丝袜| 不卡av一区二区三区| 国产欧美日韩一区二区三区在线| 91精品三级在线观看| 亚洲精品视频女| 自拍欧美九色日韩亚洲蝌蚪91| 嫩草影院入口| 亚洲国产最新在线播放| 99久久人妻综合| 老汉色av国产亚洲站长工具| 免费高清在线观看日韩| 国产精品久久久久久人妻精品电影 | 1024视频免费在线观看| avwww免费| 丝瓜视频免费看黄片| 秋霞在线观看毛片| 亚洲国产最新在线播放| 国产av国产精品国产| 欧美中文综合在线视频| 亚洲国产av影院在线观看| 久久精品亚洲av国产电影网| 美女脱内裤让男人舔精品视频| 久久久久精品人妻al黑| 人妻一区二区av| 深夜精品福利| 如日韩欧美国产精品一区二区三区| 美女扒开内裤让男人捅视频| 亚洲综合精品二区| 在线观看人妻少妇| 老司机影院毛片| 人人妻人人澡人人爽人人夜夜| 91国产中文字幕| 男女国产视频网站| 黄网站色视频无遮挡免费观看| 日韩电影二区| 国产成人精品无人区| 最近最新中文字幕免费大全7| 少妇猛男粗大的猛烈进出视频| 黄色 视频免费看| 秋霞伦理黄片| 1024视频免费在线观看| 91精品国产国语对白视频| 不卡av一区二区三区| 亚洲精品av麻豆狂野| 最近手机中文字幕大全| 国产精品一二三区在线看| 伊人久久国产一区二区| 亚洲成人av在线免费| 看非洲黑人一级黄片| 国精品久久久久久国模美| 韩国高清视频一区二区三区| 看免费成人av毛片| 亚洲美女视频黄频| 人人妻人人添人人爽欧美一区卜| 国产精品久久久久久精品古装| 久久久久精品人妻al黑| 男的添女的下面高潮视频| 欧美另类一区| 晚上一个人看的免费电影| 伦理电影大哥的女人| 国产精品秋霞免费鲁丝片| 亚洲国产最新在线播放| 国产免费又黄又爽又色| 熟女av电影| 欧美黑人欧美精品刺激| 在线观看一区二区三区激情| 国产精品蜜桃在线观看| 国产精品二区激情视频| 99热网站在线观看| 侵犯人妻中文字幕一二三四区| 午夜福利乱码中文字幕| 一区在线观看完整版| 亚洲婷婷狠狠爱综合网| 精品人妻在线不人妻| 又大又爽又粗| 亚洲国产欧美网| 国产又爽黄色视频| 如日韩欧美国产精品一区二区三区| 中国国产av一级| 亚洲精品国产区一区二| 国产在线免费精品| 国产一区亚洲一区在线观看| 蜜桃国产av成人99| 你懂的网址亚洲精品在线观看| 亚洲国产中文字幕在线视频| 伊人久久大香线蕉亚洲五| 人体艺术视频欧美日本| 久久影院123| 如何舔出高潮| 亚洲欧美一区二区三区黑人| 午夜免费鲁丝| 午夜影院在线不卡| 成人手机av| 毛片一级片免费看久久久久| 秋霞伦理黄片| 久久久久精品久久久久真实原创| 亚洲国产精品成人久久小说| 九草在线视频观看| 国产激情久久老熟女| 国产日韩欧美视频二区| 亚洲成人av在线免费| 日本猛色少妇xxxxx猛交久久| 妹子高潮喷水视频| 亚洲欧美精品自产自拍| 最近中文字幕高清免费大全6| 亚洲国产精品国产精品| 十八禁人妻一区二区| av在线播放精品| 久久久国产欧美日韩av| 黄片无遮挡物在线观看| 不卡av一区二区三区| a级片在线免费高清观看视频| 少妇被粗大猛烈的视频| 黄色视频在线播放观看不卡| 女性生殖器流出的白浆| 国产精品 欧美亚洲| 国产无遮挡羞羞视频在线观看| 国产精品一区二区在线观看99| 精品久久蜜臀av无| 如何舔出高潮| 亚洲精品自拍成人| 欧美日韩国产mv在线观看视频| 欧美黑人精品巨大| 日本wwww免费看| 国产精品久久久久成人av| 亚洲人成网站在线观看播放| 日韩制服骚丝袜av| 久久人人爽人人片av| 亚洲中文av在线| 自拍欧美九色日韩亚洲蝌蚪91| 国产日韩欧美视频二区| 大话2 男鬼变身卡| 91精品国产国语对白视频| 亚洲国产看品久久| 视频区图区小说| 在现免费观看毛片| 精品酒店卫生间| 大片电影免费在线观看免费| 嫩草影视91久久| 成人毛片60女人毛片免费| 亚洲情色 制服丝袜| 亚洲精品一区蜜桃| 国产黄频视频在线观看| 免费日韩欧美在线观看| 美女高潮到喷水免费观看| 777米奇影视久久| 性色av一级| 99国产综合亚洲精品| 男女边吃奶边做爰视频| 免费黄频网站在线观看国产| av免费观看日本| 又黄又粗又硬又大视频| 亚洲精品美女久久av网站| 各种免费的搞黄视频| 激情视频va一区二区三区| 久久女婷五月综合色啪小说| 国产高清国产精品国产三级| 欧美av亚洲av综合av国产av | 亚洲精品久久成人aⅴ小说| 亚洲熟女精品中文字幕| 欧美国产精品va在线观看不卡| 国产探花极品一区二区| 欧美中文综合在线视频| 新久久久久国产一级毛片| 亚洲一区二区三区欧美精品| 一本大道久久a久久精品| 亚洲国产精品999| 女人精品久久久久毛片| 国产欧美日韩一区二区三区在线| 国产黄色免费在线视频| 久久精品亚洲熟妇少妇任你| 看十八女毛片水多多多| 国产亚洲最大av| 亚洲精品国产av蜜桃| 亚洲av日韩在线播放| 天美传媒精品一区二区| 两个人免费观看高清视频| 成人国产麻豆网| 九色亚洲精品在线播放| 你懂的网址亚洲精品在线观看| 人人妻人人添人人爽欧美一区卜| 免费看不卡的av| 99热网站在线观看| 久久久久久人妻| 国产精品女同一区二区软件| 日韩一区二区三区影片| 国产精品一区二区在线观看99| 在线观看国产h片| 亚洲图色成人| 欧美xxⅹ黑人| 美国免费a级毛片| 在线精品无人区一区二区三| 欧美少妇被猛烈插入视频| 亚洲av在线观看美女高潮| 欧美亚洲 丝袜 人妻 在线| 女人爽到高潮嗷嗷叫在线视频| 一本一本久久a久久精品综合妖精| 久久久久国产一级毛片高清牌| 狠狠精品人妻久久久久久综合| 一级片免费观看大全| netflix在线观看网站| 国产精品二区激情视频| 成人18禁高潮啪啪吃奶动态图| 新久久久久国产一级毛片| 欧美人与性动交α欧美精品济南到| 在线 av 中文字幕| 叶爱在线成人免费视频播放| 国产免费视频播放在线视频| 九草在线视频观看| 免费观看性生交大片5| 80岁老熟妇乱子伦牲交| 高清av免费在线| 新久久久久国产一级毛片| 如日韩欧美国产精品一区二区三区| 丝袜脚勾引网站| 久久国产精品大桥未久av| 久久热在线av| kizo精华| 男女边吃奶边做爰视频| 一区福利在线观看| 大码成人一级视频| 国产亚洲av片在线观看秒播厂| 欧美日韩福利视频一区二区| 国产精品久久久久久精品古装| 男女边摸边吃奶| 一本大道久久a久久精品| videosex国产| 卡戴珊不雅视频在线播放| 亚洲成人国产一区在线观看 | 久久精品亚洲熟妇少妇任你| av.在线天堂| 国产成人免费无遮挡视频| 亚洲成人国产一区在线观看 | 亚洲精华国产精华液的使用体验| 国产亚洲欧美精品永久| 少妇精品久久久久久久| 操美女的视频在线观看| 黑人巨大精品欧美一区二区蜜桃| 在线观看免费午夜福利视频| 婷婷成人精品国产| 波多野结衣一区麻豆| 日韩一区二区三区影片| 国产黄色免费在线视频| 99热全是精品| 只有这里有精品99| 国产亚洲av高清不卡| 最近最新中文字幕大全免费视频 | 亚洲伊人色综图| 欧美人与性动交α欧美精品济南到| 青春草国产在线视频| 伊人亚洲综合成人网| 夜夜骑夜夜射夜夜干| 成人免费观看视频高清| 啦啦啦视频在线资源免费观看| 亚洲第一青青草原| 成人影院久久| 国产午夜精品一二区理论片| 免费av中文字幕在线| 日本欧美视频一区| 可以免费在线观看a视频的电影网站 | 亚洲国产精品999| 男女国产视频网站| 午夜福利一区二区在线看| 在线观看三级黄色| 搡老岳熟女国产| 日本vs欧美在线观看视频| 纯流量卡能插随身wifi吗| 国产精品一国产av| 女人被躁到高潮嗷嗷叫费观| 国产欧美日韩综合在线一区二区| 亚洲国产成人一精品久久久| 日韩伦理黄色片| 久久精品久久久久久久性| 国产精品国产三级专区第一集| 久热爱精品视频在线9| 成人毛片60女人毛片免费| 免费日韩欧美在线观看| 欧美精品一区二区大全| 久热爱精品视频在线9| 精品福利永久在线观看| 97精品久久久久久久久久精品| 男女下面插进去视频免费观看| 夜夜骑夜夜射夜夜干| 欧美日韩亚洲综合一区二区三区_| 9191精品国产免费久久| 人体艺术视频欧美日本| 99久国产av精品国产电影| 国产av码专区亚洲av| 一级毛片我不卡| 伦理电影大哥的女人| 一二三四中文在线观看免费高清| 伦理电影大哥的女人| 卡戴珊不雅视频在线播放| av有码第一页| 日韩欧美精品免费久久| 精品国产露脸久久av麻豆| 午夜日韩欧美国产| 美女脱内裤让男人舔精品视频| 一区二区三区乱码不卡18| 亚洲人成77777在线视频| 卡戴珊不雅视频在线播放| 日韩熟女老妇一区二区性免费视频| 国产1区2区3区精品| 欧美日本中文国产一区发布| 国产成人a∨麻豆精品| 欧美激情极品国产一区二区三区| 男男h啪啪无遮挡| 美女高潮到喷水免费观看| 国产1区2区3区精品| 色94色欧美一区二区| 日日摸夜夜添夜夜爱| 老汉色av国产亚洲站长工具| 中文字幕av电影在线播放| 九草在线视频观看| 9191精品国产免费久久| 肉色欧美久久久久久久蜜桃| 老司机亚洲免费影院| 亚洲欧美一区二区三区黑人| 日韩 亚洲 欧美在线| 母亲3免费完整高清在线观看| 日韩 亚洲 欧美在线| 亚洲精品自拍成人| 亚洲精品国产区一区二| 一区二区av电影网| 丁香六月天网| 中国三级夫妇交换| 国产精品麻豆人妻色哟哟久久| 一本—道久久a久久精品蜜桃钙片| 精品亚洲乱码少妇综合久久| 各种免费的搞黄视频| 亚洲国产看品久久| 国产精品一区二区在线不卡| 黄片小视频在线播放| 黑人欧美特级aaaaaa片| 日本黄色日本黄色录像| 9191精品国产免费久久| 在线看a的网站| 亚洲国产精品成人久久小说| 亚洲人成电影观看| 欧美 日韩 精品 国产| 亚洲欧洲国产日韩| av在线老鸭窝| svipshipincom国产片| av视频免费观看在线观看| 久久久久精品国产欧美久久久 | 亚洲熟女毛片儿| 大香蕉久久成人网| 国产精品蜜桃在线观看| 国产成人91sexporn| 久久精品国产a三级三级三级| 中国三级夫妇交换| 国产av一区二区精品久久| 男女高潮啪啪啪动态图| 老汉色av国产亚洲站长工具| 欧美日韩国产mv在线观看视频| 欧美xxⅹ黑人| 一边摸一边做爽爽视频免费| 香蕉国产在线看| 国产欧美日韩一区二区三区在线| 国产97色在线日韩免费| 久热这里只有精品99| 9色porny在线观看| 久久精品人人爽人人爽视色| 婷婷色综合www| 亚洲av日韩在线播放| 一级a爱视频在线免费观看| 国产精品成人在线| 天天躁夜夜躁狠狠躁躁| 最新在线观看一区二区三区 | 中文字幕高清在线视频| 久久毛片免费看一区二区三区| 国产成人91sexporn| 老汉色∧v一级毛片| 天天操日日干夜夜撸| 五月天丁香电影| 亚洲精品国产区一区二| 亚洲人成77777在线视频| av在线播放精品| 人人妻人人爽人人添夜夜欢视频| 国产精品欧美亚洲77777| 亚洲国产毛片av蜜桃av| 激情五月婷婷亚洲| 午夜福利视频在线观看免费| 日本vs欧美在线观看视频| 久久人人97超碰香蕉20202| 国产黄色免费在线视频| 嫩草影院入口| 在线观看免费视频网站a站| 中文字幕色久视频| 亚洲美女搞黄在线观看| 中文欧美无线码| 日韩中文字幕欧美一区二区 | 精品国产露脸久久av麻豆| 亚洲国产欧美网| 午夜福利免费观看在线| 校园人妻丝袜中文字幕| 一本一本久久a久久精品综合妖精| 久久影院123| 国产又色又爽无遮挡免| 成人三级做爰电影| 97在线人人人人妻| 一区二区av电影网| 免费黄网站久久成人精品| 中文字幕人妻熟女乱码| 热re99久久精品国产66热6| 国产黄色免费在线视频| 精品少妇一区二区三区视频日本电影 | 汤姆久久久久久久影院中文字幕| 中文字幕亚洲精品专区| 爱豆传媒免费全集在线观看| 亚洲欧美成人精品一区二区| 亚洲婷婷狠狠爱综合网| 亚洲激情五月婷婷啪啪| 日韩制服丝袜自拍偷拍| 夫妻午夜视频| 精品酒店卫生间| 久久久久网色| 亚洲一卡2卡3卡4卡5卡精品中文| av在线老鸭窝| 天堂8中文在线网| 国产精品av久久久久免费| 一级a爱视频在线免费观看| 日韩伦理黄色片| 日韩视频在线欧美| 91老司机精品| 亚洲人成77777在线视频| 欧美精品亚洲一区二区| 日韩精品有码人妻一区| 日韩 欧美 亚洲 中文字幕| 日本av手机在线免费观看| 国产熟女欧美一区二区| 一级黄片播放器| 日本一区二区免费在线视频| 欧美97在线视频| 国产一区二区 视频在线| 最近中文字幕2019免费版| 国产精品99久久99久久久不卡 | 亚洲欧美色中文字幕在线| 国产一区二区 视频在线| 叶爱在线成人免费视频播放| 国产精品一二三区在线看| 操美女的视频在线观看| 欧美激情 高清一区二区三区| 精品少妇内射三级| 亚洲国产欧美一区二区综合| www.自偷自拍.com| 午夜激情久久久久久久| 黄色一级大片看看| 国产精品99久久99久久久不卡 | 精品国产一区二区三区四区第35| 久久99热这里只频精品6学生| 亚洲国产精品一区二区三区在线| 飞空精品影院首页| 在线观看三级黄色| www.av在线官网国产| 9热在线视频观看99| 欧美人与性动交α欧美精品济南到| 在线观看一区二区三区激情| 国产成人av激情在线播放| 国产熟女欧美一区二区| 热re99久久精品国产66热6| 可以免费在线观看a视频的电影网站 | 女人精品久久久久毛片| 另类亚洲欧美激情| 色播在线永久视频| 欧美日韩成人在线一区二区| 国产日韩欧美亚洲二区| 久久久久久久久免费视频了| 欧美日韩一区二区视频在线观看视频在线| 妹子高潮喷水视频| 成人三级做爰电影| 成人免费观看视频高清| 欧美精品亚洲一区二区| 成人国产麻豆网| 国产精品国产av在线观看| 日韩欧美精品免费久久| 91aial.com中文字幕在线观看| 巨乳人妻的诱惑在线观看| 国精品久久久久久国模美| 午夜福利网站1000一区二区三区| 男人舔女人的私密视频| 亚洲精品国产色婷婷电影| 色网站视频免费| 精品国产超薄肉色丝袜足j| 成人国语在线视频| 欧美日本中文国产一区发布| 黄色视频在线播放观看不卡| 性少妇av在线| 熟女av电影| 日本欧美国产在线视频| 男人操女人黄网站| 午夜福利乱码中文字幕| 99九九在线精品视频| 亚洲美女黄色视频免费看| 精品国产一区二区三区四区第35| 91精品伊人久久大香线蕉| 18在线观看网站| 亚洲国产毛片av蜜桃av| 中文字幕av电影在线播放| 视频在线观看一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲,一卡二卡三卡| 日韩av在线免费看完整版不卡| 精品亚洲乱码少妇综合久久| 天堂8中文在线网| 交换朋友夫妻互换小说| 日韩制服丝袜自拍偷拍| 9色porny在线观看| 亚洲精品aⅴ在线观看| 亚洲专区中文字幕在线 | 久久鲁丝午夜福利片| 亚洲欧美中文字幕日韩二区| 日韩一卡2卡3卡4卡2021年| 国产精品二区激情视频| 黄色视频在线播放观看不卡| 国产成人免费观看mmmm| 欧美日本中文国产一区发布| 免费黄色在线免费观看| 可以免费在线观看a视频的电影网站 | 亚洲欧洲精品一区二区精品久久久 | 黄色一级大片看看| 侵犯人妻中文字幕一二三四区| 18禁裸乳无遮挡动漫免费视频| 精品少妇黑人巨大在线播放| 久久久久精品性色| 桃花免费在线播放| 99re6热这里在线精品视频| 九色亚洲精品在线播放| 亚洲一级一片aⅴ在线观看| 大香蕉久久网| 一级黄片播放器| 中文精品一卡2卡3卡4更新|