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

    沖擊響應(yīng)時(shí)域法測(cè)量轉(zhuǎn)子的一階臨界轉(zhuǎn)速*

    2014-07-31 20:40:04劉高進(jìn)馬會(huì)防章躍洪陸勇星
    機(jī)械研究與應(yīng)用 2014年1期
    關(guān)鍵詞:一階振型時(shí)域

    劉高進(jìn), 馬會(huì)防, 章躍洪, 陸勇星, 張 恒

    (1.金華職業(yè)技術(shù)學(xué)院 機(jī)電工程學(xué)院,浙江 金華 321017;2.上海凱泉泵業(yè)( 集團(tuán)) 有限公司 技術(shù)中心,上海 201804)

    沖擊響應(yīng)時(shí)域法測(cè)量轉(zhuǎn)子的一階臨界轉(zhuǎn)速*

    劉高進(jìn)1, 馬會(huì)防2, 章躍洪1, 陸勇星1, 張 恒1

    (1.金華職業(yè)技術(shù)學(xué)院 機(jī)電工程學(xué)院,浙江 金華 321017;2.上海凱泉泵業(yè)( 集團(tuán)) 有限公司 技術(shù)中心,上海 201804)

    通過(guò)對(duì)沖擊響應(yīng)時(shí)域法的研究,提出一種更加適合中小企業(yè)的轉(zhuǎn)子臨界轉(zhuǎn)速測(cè)定方法。首先從理論上分析了沖擊激勵(lì)點(diǎn)與響應(yīng)測(cè)量點(diǎn)的選擇對(duì)分析沖擊響應(yīng)具有重要作用,激勵(lì)點(diǎn)的選擇可以控制模態(tài)振型的組合,而測(cè)點(diǎn)位置的選擇能夠反映特定階次振型的振動(dòng);然后對(duì)某型號(hào)多級(jí)離心泵轉(zhuǎn)子進(jìn)行多次沖擊響應(yīng)測(cè)定實(shí)驗(yàn),通過(guò)對(duì)比頻域信號(hào)確定了最佳激勵(lì)點(diǎn)和最佳測(cè)點(diǎn)的位置,還探討了最佳激勵(lì)點(diǎn)和最佳測(cè)點(diǎn)與振型的關(guān)系;最后對(duì)多個(gè)激勵(lì)點(diǎn)與測(cè)點(diǎn)測(cè)得的時(shí)域信號(hào)進(jìn)行了一階臨界轉(zhuǎn)速的測(cè)算,結(jié)果證明最佳激勵(lì)點(diǎn)和最佳測(cè)點(diǎn)測(cè)的效果最好。該研究為轉(zhuǎn)子臨界轉(zhuǎn)速的測(cè)量提供了一種較為簡(jiǎn)易的測(cè)量方法,具有重要的應(yīng)用和參考價(jià)值。

    最佳激勵(lì)點(diǎn);最佳測(cè)點(diǎn);臨界轉(zhuǎn)速;時(shí)域法;沖擊響應(yīng)

    1 引 言

    臨界轉(zhuǎn)速,對(duì)于旋轉(zhuǎn)機(jī)械來(lái)說(shuō),是非常重要的動(dòng)態(tài)性能參數(shù)。特別是一階臨界轉(zhuǎn)速的測(cè)定,應(yīng)用更為廣泛,如水泵等很多設(shè)備,都在一階臨界轉(zhuǎn)速以下工作[1-6]。

    臨界轉(zhuǎn)速的測(cè)定,常見(jiàn)方法有直接測(cè)定法,即共振法,利用調(diào)速電機(jī),帶動(dòng)轉(zhuǎn)子不斷加速旋轉(zhuǎn),跨越共振區(qū),振幅最大時(shí)對(duì)應(yīng)的轉(zhuǎn)速,即為臨界轉(zhuǎn)速;間接測(cè)定,根據(jù)文獻(xiàn)[7]指出的“旋轉(zhuǎn)機(jī)械的臨界轉(zhuǎn)速在數(shù)值上和轉(zhuǎn)軸的橫向振動(dòng)固有頻率相同”,可以通過(guò)測(cè)量一階橫向振動(dòng)頻率,來(lái)測(cè)量臨界轉(zhuǎn)速,而一階橫向振動(dòng)頻率,可以通過(guò)模態(tài)試驗(yàn)測(cè)得。

    不論是共振法還是模態(tài)試驗(yàn),試驗(yàn)成本都較高,因?yàn)樾枰{(diào)速電機(jī)、模態(tài)測(cè)試系統(tǒng)等,對(duì)于中小企業(yè)來(lái)說(shuō),實(shí)行起來(lái)比較困難,很需要一種實(shí)驗(yàn)成本更低、更加簡(jiǎn)單易行的測(cè)試方法。

    筆者提出的沖擊時(shí)域響應(yīng)法,所需設(shè)備簡(jiǎn)單,只需要一個(gè)單通道的測(cè)振設(shè)備,試驗(yàn)成本低。在振動(dòng)學(xué)理論知識(shí)中,通過(guò)單自由度系統(tǒng)的時(shí)域沖擊響應(yīng)——測(cè)出響應(yīng)信號(hào)中的周期T,可以計(jì)算出單自由度系統(tǒng)的固有頻率。該文提出的測(cè)量轉(zhuǎn)子臨界轉(zhuǎn)速的時(shí)域響應(yīng)法,就是這一理論知識(shí)的拓展和應(yīng)用。

    以某型號(hào)的多級(jí)離心泵轉(zhuǎn)子為應(yīng)用和研究對(duì)象,該轉(zhuǎn)子有五級(jí)葉輪、一個(gè)聯(lián)軸器,驅(qū)動(dòng)端采用一個(gè)單列圓柱滾子軸承、非驅(qū)動(dòng)端采用一對(duì)角接觸球軸承,主要結(jié)構(gòu)和尺寸如圖1所示。

    圖1 轉(zhuǎn)子結(jié)構(gòu)簡(jiǎn)圖

    2 理論分析

    對(duì)于單自由度系統(tǒng),沖擊的實(shí)際作用相當(dāng)于給單自由度系統(tǒng)施加了初始位移和初始速度。然后諧振子在初始位移和初始速度產(chǎn)生的能量下,以固有頻率做自由衰減振動(dòng)[8]。

    對(duì)于轉(zhuǎn)子沖擊的實(shí)際作用相當(dāng)于給轉(zhuǎn)子施加了初始變形和初始速度,然后轉(zhuǎn)子在初始位移和初始速度產(chǎn)生的能量下,各階模態(tài)振型按不同的參與系數(shù),組合成一種振動(dòng)形式,隨著響應(yīng)時(shí)間,振幅逐漸衰減為零,其實(shí)質(zhì)是參與振動(dòng)的各階振型,各自進(jìn)行自由衰減振動(dòng),且高階振型衰減較快,低階振型衰減較慢[9]。

    測(cè)量沖擊響應(yīng),首先需要選擇沖擊激勵(lì)點(diǎn),以及響應(yīng)測(cè)量點(diǎn)。當(dāng)激勵(lì)點(diǎn)不同時(shí),轉(zhuǎn)子的初始變形不同,隨之產(chǎn)生的組合振動(dòng)中的各個(gè)模態(tài)振型所占的比重不同,即各階振型的參與系數(shù)不同,所以可通過(guò)激勵(lì)點(diǎn)的選擇,控制模態(tài)振型的組合,從而使特定的振型成為振動(dòng)的主要部分。

    沖擊產(chǎn)生的振動(dòng),是多階模態(tài)振型的組合,由于振型有節(jié)點(diǎn)的存在,節(jié)點(diǎn)處振幅為零,所以轉(zhuǎn)子上某階振型的節(jié)點(diǎn)處不能反映該階振型的振動(dòng),可以說(shuō)轉(zhuǎn)子有些位置對(duì)有些階次的振型的振動(dòng)不敏感,而有些位置則較為敏感,可以通過(guò)測(cè)點(diǎn)位置的選擇,使測(cè)得的信號(hào)能夠反映特定階次的振型的振動(dòng)。

    由于一階臨界轉(zhuǎn)速對(duì)應(yīng)一階橫向振動(dòng)頻率,因此,激勵(lì)點(diǎn)的選擇應(yīng)主要激發(fā)一階橫向振動(dòng)為目的;測(cè)點(diǎn)的選擇,應(yīng)以能夠敏感地反映一階橫向振動(dòng)為目的。

    3 實(shí) 驗(yàn)

    設(shè)置4個(gè)測(cè)點(diǎn),位于轉(zhuǎn)子下方,在軸向上均布,軸承支撐跨內(nèi)3個(gè)測(cè)點(diǎn),聯(lián)軸器上1個(gè);設(shè)置4個(gè)激勵(lì)點(diǎn),位于轉(zhuǎn)子上方,與測(cè)點(diǎn)一一對(duì)應(yīng)、且容易激勵(lì)的位置,如圖2所示。

    圖2 激勵(lì)點(diǎn)與測(cè)點(diǎn)

    采用加速度傳感器,設(shè)置采樣頻率fs=1 280 Hz,采樣點(diǎn)數(shù)N=4 096,則采樣時(shí)間間隔t=1/fss,頻域分辨率為fs/N=0.3 125 Hz,

    在各個(gè)激勵(lì)點(diǎn)用普通橡膠錘,沿垂直轉(zhuǎn)軸的方向進(jìn)行錘擊,并在各個(gè)測(cè)點(diǎn)采集沖擊響應(yīng)的時(shí)域信號(hào)。

    為研究不同測(cè)試方案的特點(diǎn),評(píng)判各工況下信號(hào)的優(yōu)劣,對(duì)時(shí)域信號(hào)進(jìn)行FFT變換——解析信號(hào)成分,得出轉(zhuǎn)子沖擊響應(yīng)的頻域圖,如表1所示。一階橫向振動(dòng)頻率,在40.625 0~43.750 0 Hz之間。

    4 對(duì)比分析及應(yīng)用

    4.1 激勵(lì)點(diǎn)的對(duì)比

    對(duì)比表1中4個(gè)激勵(lì)點(diǎn)對(duì)應(yīng)的4列,其中,激勵(lì)點(diǎn)為b點(diǎn)時(shí),不論哪個(gè)測(cè)點(diǎn),響應(yīng)信號(hào)成分都比較單純,這表明振動(dòng)主要以一階橫向振動(dòng)為主。

    對(duì)比表明,激勵(lì)點(diǎn)b是最佳激勵(lì)點(diǎn),它是前面提到的“通過(guò)激勵(lì)點(diǎn)的選擇,控制模態(tài)振型的組合,從而使特定的振型成為振動(dòng)的主要部分”的具體應(yīng)用,該激勵(lì)點(diǎn)能使一階橫向振型成為沖擊響應(yīng)振動(dòng)中的主要部分。

    4.2 測(cè)點(diǎn)的對(duì)比

    對(duì)比表1中4個(gè)測(cè)點(diǎn)對(duì)應(yīng)的4行,其中,測(cè)點(diǎn)2所對(duì)應(yīng)的一行,不論激勵(lì)點(diǎn)在哪里,其響應(yīng)信號(hào)都比較單純,且第一階橫向振動(dòng)的相對(duì)振幅最大,這表明振動(dòng)主要以一階橫向振動(dòng)為主。

    對(duì)比表明,測(cè)點(diǎn)2是最佳測(cè)點(diǎn),它是“通過(guò)測(cè)點(diǎn)的選擇,使測(cè)得的信號(hào)主要反映特定振型引起的振動(dòng)”的具體應(yīng)用,該測(cè)點(diǎn)測(cè)得的信號(hào),主要是組合振動(dòng)中的一階橫向振型的成分。

    表1中,最佳激勵(lì)點(diǎn)對(duì)應(yīng)的列與行相比較,前者的信號(hào)成分更為單純一些,畢竟,最佳激勵(lì)點(diǎn)使實(shí)際產(chǎn)生的振動(dòng)以一階為主,最佳測(cè)點(diǎn)只是“被動(dòng)地”保證采集的信號(hào)對(duì)一階敏感,故最佳激勵(lì)點(diǎn)較最佳測(cè)點(diǎn)更為重要。

    4.3 激勵(lì)點(diǎn)和測(cè)點(diǎn)與振型的關(guān)系

    最佳激勵(lì)點(diǎn)、最佳測(cè)點(diǎn)都在轉(zhuǎn)子中部,而該轉(zhuǎn)子一階橫向振動(dòng)的振型中,也是中部相對(duì)位移最大,故可以推論一階臨界轉(zhuǎn)速的最佳激勵(lì)點(diǎn)和測(cè)點(diǎn),應(yīng)選擇在一階橫向振動(dòng)的振型的相對(duì)位移最大處。

    表1 轉(zhuǎn)子沖擊響應(yīng)頻譜圖

    轉(zhuǎn)子結(jié)構(gòu)相對(duì)簡(jiǎn)單,可根據(jù)經(jīng)驗(yàn)或者理論知識(shí),估出一階橫向振動(dòng)的振型,模型復(fù)雜時(shí),難以估計(jì),可以實(shí)測(cè)幾個(gè)點(diǎn),另外,也可運(yùn)用限元仿真技術(shù),計(jì)算出振型,作為參考,圖3是用例結(jié)構(gòu)的一階橫向振動(dòng)振型的仿真計(jì)算結(jié)果。

    圖3 仿真計(jì)算的一階振型

    一階橫向振動(dòng)的振型,與其自身的結(jié)構(gòu)和受到支撐有關(guān),振型的相對(duì)位移最大處,未必在轉(zhuǎn)子中部,也可能出現(xiàn)在其它地方,比如端部,圖4是某型號(hào)真空泵一階橫向振動(dòng)振型的仿真計(jì)算結(jié)果。

    圖4 真空泵轉(zhuǎn)子的一階振型

    根據(jù)振型,最佳激勵(lì)點(diǎn)、最佳測(cè)點(diǎn)就應(yīng)該選擇在左端的二級(jí)葉輪上。

    4.4 時(shí)域法的應(yīng)用

    通過(guò)上面的對(duì)比確定了最佳測(cè)點(diǎn)、最佳激勵(lì)點(diǎn)的位置,接下來(lái)應(yīng)用時(shí)域法測(cè)量轉(zhuǎn)子的一階臨界轉(zhuǎn)速。

    圖1所示的轉(zhuǎn)子,采用最佳測(cè)點(diǎn)2、最佳激勵(lì)點(diǎn)b,為獲得良好的時(shí)域信號(hào),避免信號(hào)失真,采樣頻率應(yīng)該稍微設(shè)置高一些,這里設(shè)采樣頻率為fs=1 280 Hz,測(cè)得時(shí)域響應(yīng)信號(hào)如圖5所示。

    圖5 最佳激勵(lì)點(diǎn)b與最佳測(cè)點(diǎn)2的時(shí)域響應(yīng)信號(hào)

    信號(hào)進(jìn)入一階橫向自由衰減振動(dòng)階段,其信號(hào)特征是連續(xù)光滑、周期不變、振幅逐漸減小。測(cè)得一個(gè)自由衰減振動(dòng)周期T有30個(gè)數(shù)據(jù)點(diǎn),每?jī)蓚€(gè)數(shù)據(jù)間隔時(shí)間Δt=1/fss,則周期T=30×Δt=30/fss,對(duì)應(yīng)的頻率為f=1/T=fs/30 =42.67 Hz,對(duì)應(yīng)的一階臨界轉(zhuǎn)速為n=60f=60×42.67=2 560 r/min。

    如果僅采用最佳激勵(lì)點(diǎn)b,不采用最佳測(cè)點(diǎn),比如采用測(cè)點(diǎn)4,則測(cè)得時(shí)域響應(yīng)信號(hào)如圖6所示。

    信號(hào)不是很快、但也能較快進(jìn)入一階橫向自由衰減振動(dòng)階段,仍然容易測(cè)出一個(gè)自由衰減振動(dòng)周期有30個(gè)數(shù)據(jù)點(diǎn), 進(jìn)而能夠計(jì)算出一階臨界轉(zhuǎn)速為2 560 r/min。

    如果僅采用最佳測(cè)點(diǎn),不采用最佳激勵(lì)點(diǎn),比如采用激勵(lì)點(diǎn)d,則測(cè)得時(shí)域響應(yīng)信號(hào)如圖7所示。

    圖6 最佳激勵(lì)點(diǎn)b與測(cè)點(diǎn)4的時(shí)域響應(yīng)信號(hào)

    圖7 激勵(lì)點(diǎn)d與最佳測(cè)點(diǎn)2的時(shí)域響應(yīng)信號(hào)

    信號(hào)特點(diǎn)同圖6基本一致,也能較快進(jìn)入一階橫向自由衰減振動(dòng)階段,測(cè)出衰減振動(dòng)周期有30個(gè)數(shù)據(jù)點(diǎn), 進(jìn)而可以計(jì)算出一階臨界轉(zhuǎn)速為2 560 r/min。

    需要注意,測(cè)算周期長(zhǎng)短時(shí),信號(hào)前面的部分,受高階振型振動(dòng)的影響較大,信號(hào)后面的部分,主要是一階橫向振動(dòng),所以,在振幅衰減幅度不影響測(cè)量的前提下,盡量采用信號(hào)后面的部分。

    如果不采用最佳激勵(lì)點(diǎn)、不采用最佳測(cè)點(diǎn),比如采用激勵(lì)點(diǎn)d、測(cè)點(diǎn)1,則測(cè)得的時(shí)域響應(yīng)信號(hào)及信號(hào)的局部放大,如圖8所示。

    圖8 激勵(lì)點(diǎn)d與測(cè)點(diǎn)1的時(shí)域響應(yīng)信號(hào)及信號(hào)的局部放大

    信號(hào)經(jīng)過(guò)較長(zhǎng)時(shí)間,也能夠進(jìn)入一階橫向自由衰減振動(dòng)階段,測(cè)出衰減振動(dòng)周期有30個(gè)數(shù)據(jù)點(diǎn), 進(jìn)而可以計(jì)算出一階臨界轉(zhuǎn)速為2 560 r/min。

    由于信號(hào)的振幅在逐漸衰減,所以,信號(hào)必須在衰減成背景噪聲信號(hào)之前,出現(xiàn)較為明顯的一階橫向自由衰減振動(dòng)特征,否則,則不能應(yīng)用沖擊響應(yīng)的時(shí)域方法測(cè)量臨界轉(zhuǎn)速。

    例如,當(dāng)采用激勵(lì)點(diǎn)d、測(cè)點(diǎn)4時(shí),測(cè)得的時(shí)域信號(hào)如圖9所示,即使放大局部信號(hào),也沒(méi)有明顯的一階橫向自由衰減振動(dòng)特征的信號(hào)出現(xiàn),所以,如果不采用最佳激勵(lì)點(diǎn)、最佳測(cè)點(diǎn),則有可能無(wú)法應(yīng)用時(shí)域法測(cè)量臨界轉(zhuǎn)速。

    圖9 激勵(lì)點(diǎn)d與測(cè)點(diǎn)4的時(shí)域響應(yīng)信號(hào)

    5 結(jié) 語(yǔ)

    (1) 從理論上分析了轉(zhuǎn)子上激勵(lì)點(diǎn)和測(cè)量點(diǎn)的不同,會(huì)影響沖擊響應(yīng)時(shí)域法的應(yīng)用。激勵(lì)點(diǎn)的選擇可控制模態(tài)振型的參與系數(shù);測(cè)點(diǎn)位置選擇使測(cè)得的信號(hào)反映特定階次振型的振動(dòng)。

    (2) 通過(guò)沖擊響應(yīng)的頻譜對(duì)比,得出了最佳激勵(lì)點(diǎn)、最佳測(cè)點(diǎn)的位置,并探討了最佳激勵(lì)點(diǎn)、最佳測(cè)點(diǎn)與振型的關(guān)系。通過(guò)沖擊響應(yīng)時(shí)域法測(cè)量轉(zhuǎn)子一階臨界轉(zhuǎn)速的測(cè)算分析,最佳激勵(lì)點(diǎn)和最佳測(cè)點(diǎn)效果最佳。

    該文的研究,為轉(zhuǎn)子臨界轉(zhuǎn)速的測(cè)量,提供了一種較為簡(jiǎn)易的測(cè)量方法,具有重要的應(yīng)用和參考價(jià)值。

    [1] 趙萬(wàn)勇,白雙寶,馬鵬飛.離心泵轉(zhuǎn)子振動(dòng)研究現(xiàn)狀與展望[J]. 流體機(jī)械, 2011, 39(3): 37-39.

    [2] 何新榮, 傅行軍. 基于ANSYS 的軸承轉(zhuǎn)子系統(tǒng)動(dòng)力特性研究[J]. 機(jī)械研究與應(yīng)用,2011(6): 39-41.

    [3] Moreira M,Antunes J,Pina H. A Theoretical Model for Nonlinear Orbital Motions of Rotors Under Fluid Confinement [J]. Journal of Fluids and Structures, 2000,14(5):635-668.

    [4] 秦 然,張建偉. 轉(zhuǎn)子臨界轉(zhuǎn)速的求解方法比較[J]. 沈陽(yáng)化工學(xué)院學(xué)報(bào), 2007, 21(1): 74-76.

    [5] Moreira M,Antunes J,Pina H. An Improved Linear Model for Rotors Subject to Dissipative Annular Flows[J]. Journal of Fluids and Structures,2003,17(6): 813-832.

    [6] Heshmat H,Tomaszewski M J,Walton J F. Small Gas Turbine Engine Operating with High-Temperature Foil Bearings[C].ASME Turbo Expo 2006: Power for Land,Sea,and Air.Barcelona,Spain: ASME,2006: 387-393.

    [7] 趙 玫,周海亭,陳光冶,等.機(jī)械振動(dòng)與噪聲學(xué)[M].北京:科學(xué)出版社,2004.

    [8] 程耀東.機(jī)械振動(dòng)學(xué)[M].杭州:浙江大學(xué)出版社,1988.

    [9] 傅志方,華宏星. 模態(tài)分析理論與應(yīng)用[M].上海:上海交通大學(xué)出版社,1988.

    The Measurement of Rotor′s First Order Critical Speed with Method of Time Domain Signal of Pulse Response

    LIU Gao-jin1, MA Hui-fang2, ZHANG Yue-hong1, LU Yong-xing1, ZHANG Heng1

    (1.Mechanical&ElectricalEngineeringCollege,JinhuaPolytechnic,JinhuaZhejiang321000,China;2.TechnologyCenterofShanghaiKaiquanGroup,Shanghai201804,China)

    The purpose of this study on method of time domain signal of pulse response is to present a way to test the first critical speed of rotor which has lower experiment cost and is more suitable for small and medium-sized enterprise. Firstly, it analyzes theoretically that the choice of excitation point and measurement point is very important for studying pulse response. The excitation point controls the form of combination of multi-order modal shapes with different participation factors and the measurement point reflects sensitivity to different order modal shape. Secondly, experiments on the rotor of a multistage centrifugal pump are conducted, and the frequency domain signals are compared among one another to identify where is the best excitation point and where is the best measurement point, and relationship between modal shape and the best excitation measurement point is discussed; Finally, the first order critical speed of rotor under different excitation measurement points is acquired with method of domain signal of pulse response and it suggests that the effect under the best excitation measurement point is optimum. This study presents an easier test method for the rotor′s first order critical speed, and has important application and reference values.

    the best excitation point; the best measurement point; critical speed; method of time domain signal; pulse response

    2013-12-16

    浙江省自然科學(xué)

    基金項(xiàng)目:《基于雙譜生成多點(diǎn)隨機(jī)振動(dòng)控制試驗(yàn)的方法研究》(編號(hào)LQ12E05011)

    劉高進(jìn)(1981-),男,福建龍巖人,講師,碩士,研究方向:數(shù)控機(jī)床故障診斷與維護(hù),振動(dòng)與沖擊。

    TH824

    A

    1007-4414(2014)01-0206-04

    猜你喜歡
    一階振型時(shí)域
    求解Ericksen-Leslie方程的一階精度、線(xiàn)性穩(wěn)定的數(shù)值格式
    李超代數(shù)到Kac模的一階上同調(diào)
    關(guān)于模態(tài)綜合法的注記
    縱向激勵(lì)下大跨鋼桁拱橋高階振型效應(yīng)分析
    塔腿加過(guò)渡段輸電塔動(dòng)力特性分析
    基于時(shí)域信號(hào)的三電平逆變器復(fù)合故障診斷
    基于極大似然準(zhǔn)則與滾動(dòng)時(shí)域估計(jì)的自適應(yīng)UKF算法
    基于時(shí)域逆濾波的寬帶脈沖聲生成技術(shù)
    一階非線(xiàn)性微分方程解法探析
    新校園(下)(2015年6期)2015-07-04 05:06:49
    結(jié)構(gòu)振型幾何辨識(shí)及應(yīng)用研究
    山西建筑(2015年14期)2015-06-05 09:37:07
    9191精品国产免费久久| 国产亚洲精品av在线| 老司机深夜福利视频在线观看| 手机成人av网站| 国产在线观看jvid| 欧美日韩亚洲综合一区二区三区_| 大香蕉久久成人网| 亚洲欧美精品综合久久99| 亚洲国产日韩欧美精品在线观看 | 欧美+亚洲+日韩+国产| 久久精品亚洲精品国产色婷小说| 黄频高清免费视频| 少妇被粗大的猛进出69影院| 美女大奶头视频| 日韩免费av在线播放| 亚洲片人在线观看| 这个男人来自地球电影免费观看| 日韩欧美在线二视频| 天天一区二区日本电影三级| 久久精品人妻少妇| 又紧又爽又黄一区二区| 看片在线看免费视频| 国产精华一区二区三区| 欧美日韩亚洲综合一区二区三区_| 亚洲五月天丁香| 婷婷丁香在线五月| 成人永久免费在线观看视频| 美女午夜性视频免费| 国产精品精品国产色婷婷| 国产成人av教育| 怎么达到女性高潮| 欧美 亚洲 国产 日韩一| 久久国产亚洲av麻豆专区| 亚洲欧洲精品一区二区精品久久久| 国产av一区在线观看免费| 成人永久免费在线观看视频| 国产伦在线观看视频一区| 中文资源天堂在线| 国产亚洲av高清不卡| 草草在线视频免费看| 久久人妻福利社区极品人妻图片| 亚洲五月婷婷丁香| 高清毛片免费观看视频网站| 国产精品 国内视频| 国产av不卡久久| 听说在线观看完整版免费高清| 女警被强在线播放| 黑人操中国人逼视频| www.自偷自拍.com| 免费高清在线观看日韩| 精品欧美一区二区三区在线| 无限看片的www在线观看| 一进一出抽搐gif免费好疼| 99久久久亚洲精品蜜臀av| 男女午夜视频在线观看| 久热爱精品视频在线9| 午夜两性在线视频| 黄片小视频在线播放| 丝袜人妻中文字幕| 亚洲人成网站高清观看| 18禁观看日本| 51午夜福利影视在线观看| 国产精品亚洲av一区麻豆| 成人永久免费在线观看视频| 此物有八面人人有两片| 免费av毛片视频| 久久精品夜夜夜夜夜久久蜜豆 | 老司机靠b影院| 波多野结衣高清作品| 99在线视频只有这里精品首页| 国产成人av激情在线播放| 亚洲 欧美 日韩 在线 免费| 国产一区二区激情短视频| 日韩有码中文字幕| √禁漫天堂资源中文www| 黄色a级毛片大全视频| 一级毛片女人18水好多| 亚洲中文日韩欧美视频| 男人舔女人的私密视频| 淫秽高清视频在线观看| 老司机午夜福利在线观看视频| 真人做人爱边吃奶动态| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品在线美女| 丝袜在线中文字幕| 亚洲自拍偷在线| 村上凉子中文字幕在线| 亚洲欧美精品综合久久99| 怎么达到女性高潮| 亚洲真实伦在线观看| a级毛片a级免费在线| 好男人在线观看高清免费视频 | 久久 成人 亚洲| 正在播放国产对白刺激| 亚洲精品av麻豆狂野| 老汉色∧v一级毛片| 久久久水蜜桃国产精品网| 欧美日韩一级在线毛片| 搞女人的毛片| 91字幕亚洲| 村上凉子中文字幕在线| 国产午夜精品久久久久久| 国产一区二区三区在线臀色熟女| 欧美激情高清一区二区三区| 精品久久久久久久末码| 久久精品成人免费网站| 在线观看免费午夜福利视频| 中文字幕精品免费在线观看视频| 国产91精品成人一区二区三区| 午夜精品在线福利| 国产亚洲精品av在线| 亚洲第一av免费看| 色精品久久人妻99蜜桃| 丁香六月欧美| 色哟哟哟哟哟哟| 国产一区在线观看成人免费| 午夜免费观看网址| 国产亚洲精品第一综合不卡| 国产精品久久久久久精品电影 | 久久香蕉精品热| 国产精品二区激情视频| 亚洲成国产人片在线观看| 久久久久国产精品人妻aⅴ院| 一区二区日韩欧美中文字幕| 1024手机看黄色片| 色播在线永久视频| 亚洲人成网站在线播放欧美日韩| 亚洲熟女毛片儿| 97碰自拍视频| 亚洲av电影在线进入| xxxwww97欧美| av免费在线观看网站| 国产爱豆传媒在线观看 | 国产黄色小视频在线观看| 很黄的视频免费| 亚洲国产日韩欧美精品在线观看 | 在线观看www视频免费| 国产私拍福利视频在线观看| 国产视频内射| 成人特级黄色片久久久久久久| 亚洲av成人av| 亚洲avbb在线观看| 老熟妇仑乱视频hdxx| 天天添夜夜摸| 日韩欧美国产在线观看| 99国产极品粉嫩在线观看| 国产欧美日韩精品亚洲av| 亚洲av电影在线进入| 男女之事视频高清在线观看| 在线观看免费日韩欧美大片| 国产伦在线观看视频一区| 日韩欧美一区视频在线观看| 国产国语露脸激情在线看| 国产免费男女视频| 亚洲精品美女久久久久99蜜臀| 美女 人体艺术 gogo| 两性午夜刺激爽爽歪歪视频在线观看 | 非洲黑人性xxxx精品又粗又长| 日本撒尿小便嘘嘘汇集6| 一级毛片女人18水好多| 男人舔女人的私密视频| 久久精品国产亚洲av高清一级| 色综合欧美亚洲国产小说| 亚洲av成人一区二区三| 69av精品久久久久久| 国内揄拍国产精品人妻在线 | 亚洲九九香蕉| 国产黄片美女视频| 亚洲国产欧美日韩在线播放| 99re在线观看精品视频| 精品乱码久久久久久99久播| 久久这里只有精品19| 美女国产高潮福利片在线看| 日本精品一区二区三区蜜桃| 色av中文字幕| 亚洲 国产 在线| 两个人视频免费观看高清| 欧美大码av| 国产高清视频在线播放一区| 老汉色av国产亚洲站长工具| 亚洲av熟女| 又大又爽又粗| 女人被狂操c到高潮| 天天一区二区日本电影三级| 久久久久久国产a免费观看| 久久久久国产一级毛片高清牌| 中文亚洲av片在线观看爽| 免费高清在线观看日韩| 黄片播放在线免费| 国产成人系列免费观看| 日本一区二区免费在线视频| 国产精品一区二区精品视频观看| 国产免费av片在线观看野外av| 一个人免费在线观看的高清视频| 母亲3免费完整高清在线观看| 亚洲人成伊人成综合网2020| 少妇 在线观看| 琪琪午夜伦伦电影理论片6080| 老司机午夜福利在线观看视频| 亚洲成人免费电影在线观看| 日本a在线网址| 黑人欧美特级aaaaaa片| 亚洲国产精品999在线| 俺也久久电影网| 看片在线看免费视频| 一区二区三区高清视频在线| 日韩欧美三级三区| 国产精品av久久久久免费| av超薄肉色丝袜交足视频| 男女那种视频在线观看| 巨乳人妻的诱惑在线观看| 色老头精品视频在线观看| 男女床上黄色一级片免费看| 国产精品久久电影中文字幕| 国产精品98久久久久久宅男小说| 1024香蕉在线观看| 国产成人影院久久av| 久久草成人影院| 黄片播放在线免费| 老司机在亚洲福利影院| 村上凉子中文字幕在线| 丝袜美腿诱惑在线| 亚洲欧美激情综合另类| 日韩国内少妇激情av| 男女下面进入的视频免费午夜 | 国产精品电影一区二区三区| 在线永久观看黄色视频| 最近在线观看免费完整版| 午夜久久久久精精品| 91九色精品人成在线观看| 12—13女人毛片做爰片一| 宅男免费午夜| 国产熟女午夜一区二区三区| 一级a爱视频在线免费观看| 露出奶头的视频| 91大片在线观看| 午夜免费成人在线视频| 给我免费播放毛片高清在线观看| 激情在线观看视频在线高清| 精品久久久久久久毛片微露脸| 91大片在线观看| 99国产极品粉嫩在线观看| 久久香蕉国产精品| 亚洲成av人片免费观看| 国产极品粉嫩免费观看在线| 免费一级毛片在线播放高清视频| 在线免费观看的www视频| 国产精品亚洲av一区麻豆| 亚洲国产欧美日韩在线播放| 男人的好看免费观看在线视频 | 在线视频色国产色| 国产麻豆成人av免费视频| 在线观看免费视频日本深夜| 免费观看人在逋| 91大片在线观看| 成人特级黄色片久久久久久久| 国产一区在线观看成人免费| 中文字幕另类日韩欧美亚洲嫩草| 一级作爱视频免费观看| 午夜日韩欧美国产| netflix在线观看网站| 两人在一起打扑克的视频| 国产精品 欧美亚洲| 黑人操中国人逼视频| 色哟哟哟哟哟哟| 亚洲精品粉嫩美女一区| 麻豆国产av国片精品| 亚洲avbb在线观看| 好男人在线观看高清免费视频 | 女性被躁到高潮视频| 久久久久久人人人人人| 久久久久久久午夜电影| 国产精品一区二区免费欧美| 国产欧美日韩一区二区三| 中出人妻视频一区二区| 久久久国产欧美日韩av| 99热6这里只有精品| 一区二区三区激情视频| 国产精品久久久久久人妻精品电影| 亚洲午夜理论影院| 伊人久久大香线蕉亚洲五| 视频区欧美日本亚洲| 日韩一卡2卡3卡4卡2021年| 色婷婷久久久亚洲欧美| 国产真人三级小视频在线观看| 国内毛片毛片毛片毛片毛片| 美女扒开内裤让男人捅视频| 777久久人妻少妇嫩草av网站| 午夜影院日韩av| 亚洲五月天丁香| 99国产精品99久久久久| 婷婷六月久久综合丁香| 在线观看免费日韩欧美大片| 精品一区二区三区av网在线观看| 不卡一级毛片| 欧美激情极品国产一区二区三区| 美女 人体艺术 gogo| 丁香六月欧美| 亚洲国产欧美日韩在线播放| 亚洲无线在线观看| av片东京热男人的天堂| 亚洲最大成人中文| 亚洲国产欧洲综合997久久, | 欧美国产精品va在线观看不卡| 人人妻,人人澡人人爽秒播| 国产成人精品无人区| 精品第一国产精品| 18禁黄网站禁片免费观看直播| 精品久久久久久久末码| 国产成人av激情在线播放| 亚洲熟女毛片儿| 一二三四社区在线视频社区8| 国产亚洲av嫩草精品影院| 久久伊人香网站| 免费在线观看影片大全网站| 亚洲五月婷婷丁香| 国产精品98久久久久久宅男小说| av超薄肉色丝袜交足视频| 少妇熟女aⅴ在线视频| 露出奶头的视频| 成人一区二区视频在线观看| 天堂√8在线中文| 精品第一国产精品| 国产成年人精品一区二区| 伦理电影免费视频| 国产三级在线视频| www.www免费av| 亚洲av中文字字幕乱码综合 | 一级作爱视频免费观看| √禁漫天堂资源中文www| 欧美精品啪啪一区二区三区| 黑丝袜美女国产一区| 啦啦啦韩国在线观看视频| 亚洲第一青青草原| 午夜久久久久精精品| 久久精品国产亚洲av香蕉五月| 午夜福利成人在线免费观看| 天天躁狠狠躁夜夜躁狠狠躁| www国产在线视频色| 国产野战对白在线观看| 狂野欧美激情性xxxx| 91成人精品电影| 久久九九热精品免费| 亚洲国产日韩欧美精品在线观看 | 久久精品aⅴ一区二区三区四区| 别揉我奶头~嗯~啊~动态视频| 麻豆成人av在线观看| av超薄肉色丝袜交足视频| 亚洲最大成人中文| 免费搜索国产男女视频| 亚洲欧美精品综合一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 亚洲人成伊人成综合网2020| 波多野结衣巨乳人妻| 人人妻,人人澡人人爽秒播| 成年人黄色毛片网站| 亚洲人成网站在线播放欧美日韩| 俺也久久电影网| 久久香蕉精品热| 99国产精品一区二区蜜桃av| 国产精品一区二区精品视频观看| 久久精品成人免费网站| 日本黄色视频三级网站网址| 久久九九热精品免费| 国产午夜福利久久久久久| 日日干狠狠操夜夜爽| 91老司机精品| 久久人人精品亚洲av| 国产精品久久视频播放| 一卡2卡三卡四卡精品乱码亚洲| 夜夜看夜夜爽夜夜摸| 51午夜福利影视在线观看| 精品久久久久久久人妻蜜臀av| 久久狼人影院| 少妇裸体淫交视频免费看高清 | 久久久久国产精品人妻aⅴ院| 亚洲 国产 在线| 国产精品,欧美在线| 亚洲国产欧美一区二区综合| 在线十欧美十亚洲十日本专区| 美女高潮喷水抽搐中文字幕| 欧美一级毛片孕妇| 欧美+亚洲+日韩+国产| 久久久久国内视频| 搞女人的毛片| 婷婷六月久久综合丁香| 久久久久久亚洲精品国产蜜桃av| 欧美成人午夜精品| 69av精品久久久久久| 无限看片的www在线观看| 国产精品免费一区二区三区在线| 真人一进一出gif抽搐免费| 日本熟妇午夜| 亚洲久久久国产精品| 免费看日本二区| 亚洲全国av大片| 免费人成视频x8x8入口观看| 最新在线观看一区二区三区| 日韩av在线大香蕉| 日本一本二区三区精品| 99re在线观看精品视频| 久久午夜亚洲精品久久| 亚洲av片天天在线观看| 88av欧美| 熟女少妇亚洲综合色aaa.| 自线自在国产av| 免费搜索国产男女视频| 一区二区三区国产精品乱码| av欧美777| 麻豆av在线久日| 欧美黑人精品巨大| 欧美乱妇无乱码| 美女国产高潮福利片在线看| 美女高潮到喷水免费观看| 国产v大片淫在线免费观看| 久久香蕉激情| 成人精品一区二区免费| 巨乳人妻的诱惑在线观看| 国产欧美日韩一区二区精品| 国产精品1区2区在线观看.| 露出奶头的视频| 国产99久久九九免费精品| 国产主播在线观看一区二区| 免费在线观看亚洲国产| 亚洲三区欧美一区| 久久99热这里只有精品18| 欧美激情极品国产一区二区三区| 国产aⅴ精品一区二区三区波| 国产一区二区三区视频了| 亚洲一码二码三码区别大吗| 日本一本二区三区精品| 性欧美人与动物交配| 国产亚洲av高清不卡| 黑丝袜美女国产一区| 一卡2卡三卡四卡精品乱码亚洲| 国产爱豆传媒在线观看 | 91老司机精品| 国产精品1区2区在线观看.| 久久中文看片网| 别揉我奶头~嗯~啊~动态视频| 俺也久久电影网| 在线观看66精品国产| 亚洲人成77777在线视频| 日日摸夜夜添夜夜添小说| 久久天堂一区二区三区四区| 女人被狂操c到高潮| 欧美日韩黄片免| 亚洲精品中文字幕一二三四区| 国产亚洲精品综合一区在线观看 | 桃红色精品国产亚洲av| 国产精品综合久久久久久久免费| 久久这里只有精品19| 亚洲真实伦在线观看| 特大巨黑吊av在线直播 | 成人欧美大片| 性色av乱码一区二区三区2| 国产黄片美女视频| 男人操女人黄网站| 日韩欧美 国产精品| 日韩高清综合在线| 亚洲va日本ⅴa欧美va伊人久久| 亚洲中文av在线| 国产激情久久老熟女| 欧美一区二区精品小视频在线| 制服诱惑二区| 国产野战对白在线观看| 国产亚洲欧美98| 2021天堂中文幕一二区在线观 | 国产av在哪里看| 97超级碰碰碰精品色视频在线观看| 韩国精品一区二区三区| 一级作爱视频免费观看| 波多野结衣av一区二区av| 国内久久婷婷六月综合欲色啪| 亚洲久久久国产精品| 国产激情欧美一区二区| 一二三四社区在线视频社区8| 免费看日本二区| 不卡一级毛片| 日本免费a在线| 亚洲欧美日韩高清在线视频| 国产av不卡久久| 香蕉久久夜色| 在线观看午夜福利视频| 久久久久久九九精品二区国产 | 黑人操中国人逼视频| 欧美性猛交黑人性爽| 黄片播放在线免费| 村上凉子中文字幕在线| 给我免费播放毛片高清在线观看| 老司机深夜福利视频在线观看| 亚洲狠狠婷婷综合久久图片| 亚洲国产精品合色在线| 18美女黄网站色大片免费观看| 午夜视频精品福利| 免费看日本二区| 在线免费观看的www视频| 人人妻人人看人人澡| 99久久无色码亚洲精品果冻| 亚洲国产欧美日韩在线播放| 欧美一级毛片孕妇| 亚洲av成人av| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利视频1000在线观看| 午夜福利18| 在线播放国产精品三级| 黄色女人牲交| 欧美国产日韩亚洲一区| 国产真人三级小视频在线观看| 亚洲人成电影免费在线| 一二三四在线观看免费中文在| 日日干狠狠操夜夜爽| 色综合婷婷激情| 夜夜爽天天搞| 自线自在国产av| 法律面前人人平等表现在哪些方面| 国产黄a三级三级三级人| 欧美国产精品va在线观看不卡| 成年人黄色毛片网站| 天天躁狠狠躁夜夜躁狠狠躁| 91老司机精品| 日本a在线网址| 久久久久久免费高清国产稀缺| 99久久精品国产亚洲精品| 欧美午夜高清在线| 在线观看舔阴道视频| 18禁裸乳无遮挡免费网站照片 | 亚洲人成77777在线视频| 黑人操中国人逼视频| 亚洲中文av在线| 91大片在线观看| 国产蜜桃级精品一区二区三区| 色av中文字幕| 亚洲av美国av| 成人永久免费在线观看视频| 两性夫妻黄色片| 美女大奶头视频| 人妻丰满熟妇av一区二区三区| 亚洲av第一区精品v没综合| 妹子高潮喷水视频| 性欧美人与动物交配| 亚洲精品色激情综合| 成人永久免费在线观看视频| 久久午夜亚洲精品久久| 精品国产一区二区三区四区第35| 亚洲精品av麻豆狂野| 夜夜爽天天搞| 国产精品久久久久久亚洲av鲁大| svipshipincom国产片| 香蕉av资源在线| 怎么达到女性高潮| 国产蜜桃级精品一区二区三区| 亚洲av电影不卡..在线观看| 亚洲国产精品sss在线观看| 欧美黄色淫秽网站| 少妇被粗大的猛进出69影院| 黑人操中国人逼视频| 禁无遮挡网站| 啦啦啦免费观看视频1| 日日摸夜夜添夜夜添小说| 亚洲五月天丁香| 精品人妻1区二区| 成人18禁高潮啪啪吃奶动态图| 欧美性猛交╳xxx乱大交人| 亚洲国产高清在线一区二区三 | 久久久久国产精品人妻aⅴ院| 精品久久久久久久末码| 国产男靠女视频免费网站| 亚洲人成网站高清观看| 国产不卡一卡二| 亚洲成a人片在线一区二区| 一区二区三区高清视频在线| 久99久视频精品免费| 每晚都被弄得嗷嗷叫到高潮| 1024手机看黄色片| 亚洲国产精品sss在线观看| 国产99久久九九免费精品| 欧美在线一区亚洲| 亚洲 欧美 日韩 在线 免费| 亚洲成人久久爱视频| 欧美又色又爽又黄视频| 可以免费在线观看a视频的电影网站| 久久精品夜夜夜夜夜久久蜜豆 | 久久国产精品人妻蜜桃| 亚洲av第一区精品v没综合| 久久久久久九九精品二区国产 | 色综合婷婷激情| 一夜夜www| 激情在线观看视频在线高清| 色精品久久人妻99蜜桃| 中文亚洲av片在线观看爽| 好男人在线观看高清免费视频 | 特大巨黑吊av在线直播 | 国产精品99久久99久久久不卡| 久久久精品国产亚洲av高清涩受| 精品乱码久久久久久99久播| 少妇熟女aⅴ在线视频| 无人区码免费观看不卡| 欧美黄色片欧美黄色片| 曰老女人黄片| 两人在一起打扑克的视频| 一二三四社区在线视频社区8| 欧美精品亚洲一区二区| aaaaa片日本免费| 国产欧美日韩精品亚洲av| 一二三四在线观看免费中文在| 首页视频小说图片口味搜索| 亚洲精品国产精品久久久不卡| 熟女电影av网| 一卡2卡三卡四卡精品乱码亚洲| 亚洲人成77777在线视频| 色哟哟哟哟哟哟| 国产av一区在线观看免费| 欧美激情高清一区二区三区| 在线观看免费视频日本深夜|