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

    道路試驗(yàn)下某重型商用車駕駛室嘯叫聲源定位

    2022-10-27 09:06:08楊紅波史文庫(kù)陳志勇趙燕燕苑仁飛
    振動(dòng)與沖擊 2022年20期
    關(guān)鍵詞:振動(dòng)

    楊紅波, 郭 磊, 史文庫(kù), 陳志勇, 趙燕燕, 劉 健, 苑仁飛

    (1. 吉林大學(xué) 汽車仿真與控制國(guó)家重點(diǎn)實(shí)驗(yàn)室, 吉林 長(zhǎng)春 130025; 2. 中國(guó)重汽集團(tuán)汽車研究總院, 濟(jì)南 250100)

    伴隨著市場(chǎng)發(fā)展,商用車公路運(yùn)輸在國(guó)民經(jīng)濟(jì)中占有舉足輕重的地位[1]?,F(xiàn)今,隨著重型商用車駕駛員向著年輕化發(fā)展,并且噪聲法規(guī)要求日益嚴(yán)格,車輛的NVH問題日益得到各汽車廠家的重視[2]。與乘用車相比,重型商用車具有如下特點(diǎn):就動(dòng)力系統(tǒng)而言,重型商用車普遍搭載柴油發(fā)動(dòng)機(jī),其輸出扭矩和轉(zhuǎn)速波動(dòng)較大;就變速箱而言,重型商用車采用多擋位主副變速箱系統(tǒng),如12擋變速系統(tǒng)、16擋變速系統(tǒng),且副箱多配置行星齒輪機(jī)構(gòu)以使空間更加緊湊;就驅(qū)動(dòng)橋而言,重型商用車普遍采用中后雙驅(qū)動(dòng)橋形式,且在中橋配置一對(duì)圓柱齒輪副實(shí)現(xiàn)動(dòng)力向后橋的傳遞[3]。通過試驗(yàn)對(duì)重型商用車的NVH問題進(jìn)行根源定位并進(jìn)行分析、評(píng)估和改善具有重要的現(xiàn)實(shí)意義,目前針對(duì)其進(jìn)行的研究相對(duì)較少,但國(guó)內(nèi)外針對(duì)其他車型以及各總成進(jìn)行的試驗(yàn)與理論分析可為本研究提供一定的借鑒和參考。

    于蓬等[4-8]針對(duì)某集中式驅(qū)動(dòng)純電動(dòng)車進(jìn)行了試驗(yàn)研究,系統(tǒng)分析了該車的動(dòng)力總成系統(tǒng)振動(dòng)特性、扭轉(zhuǎn)振動(dòng)特性、整車抖動(dòng)特性及噪聲與振動(dòng)特性,針對(duì)試驗(yàn)結(jié)果,對(duì)該車的NVH問題進(jìn)行了總結(jié)與分析,并提出了相應(yīng)的改進(jìn)措施。張邦基等[9]針對(duì)某搭載行星減速器的純電動(dòng)客車進(jìn)行了實(shí)車振動(dòng)試驗(yàn),評(píng)價(jià)了行星減速器對(duì)純電動(dòng)客車振動(dòng)舒適性的影響。張新剛[10]基于試驗(yàn)研究,對(duì)某客車在怠速工況下的振動(dòng)噪聲源進(jìn)行了識(shí)別,并從動(dòng)力總成懸置的角度進(jìn)行了優(yōu)化。Zeng等[11]針對(duì)某搭載多擋自動(dòng)變速器的純電動(dòng)客車進(jìn)行了試驗(yàn)研究,對(duì)造成車內(nèi)地板抖動(dòng)的原因進(jìn)行了定位并提出了改進(jìn)措施。王媛文等[12]針對(duì)某前置后驅(qū)車輛在低轉(zhuǎn)速時(shí)出現(xiàn)的轟鳴聲進(jìn)行了試驗(yàn),分析了該問題的激振源與傳遞路徑并提出了相應(yīng)的改進(jìn)措施。相龍洋等[13]針對(duì)某汽車手動(dòng)變速器,對(duì)引發(fā)變速器箱體表面振動(dòng)噪聲的問題源進(jìn)行了識(shí)別,并給出了兩種優(yōu)化途徑。張容川等[14]針對(duì)某搭載CVT的車輛進(jìn)行了噪聲測(cè)試,識(shí)別出了行星齒輪機(jī)構(gòu)的嘯叫階次并有針對(duì)性地進(jìn)行了齒輪修形優(yōu)化。Feng等[15]對(duì)非穩(wěn)態(tài)運(yùn)行工況下的行星齒輪箱進(jìn)行了振動(dòng)試驗(yàn),并對(duì)最后的診斷結(jié)果進(jìn)行了分析,給出了優(yōu)化方案。Luo等[16]以某一典型的特種車輛為例,通過試驗(yàn)分析了車輛內(nèi)部噪聲的頻譜特性,對(duì)其聲源進(jìn)行了定位。Dunai等[17]研究了寶馬530d車型在室外高速公路以不同車速行駛時(shí)由抖振引起的噪聲問題,并對(duì)聲源進(jìn)行了定位,確定了抖振噪聲效應(yīng)與聽者的相對(duì)位置關(guān)系。Deng等[18]對(duì)某車輛駕駛室內(nèi)出現(xiàn)的異常噪聲進(jìn)行了識(shí)別,確定了輪胎聲腔共振對(duì)其起主要貢獻(xiàn)量,并對(duì)輪胎進(jìn)行了聲學(xué)處理,改善了車內(nèi)噪聲問題。

    上述研究針對(duì)汽車整車和各分總成進(jìn)行了一定的試驗(yàn),并取得了一定的研究成果,但幾乎沒有研究者針對(duì)重型商用車進(jìn)行系統(tǒng)的NVH試驗(yàn)并進(jìn)行相應(yīng)的振動(dòng)噪聲源定位。重型商用車具有前述所討論的不同于乘用車的特點(diǎn),針對(duì)具體問題,對(duì)其進(jìn)行試驗(yàn)研究具體深刻的實(shí)際意義。本文針對(duì)某重型商用車在7擋加速工況下駕駛室內(nèi)出現(xiàn)的嘯叫噪聲問題,對(duì)動(dòng)力傳動(dòng)系統(tǒng)進(jìn)行系統(tǒng)性分析,對(duì)樣車進(jìn)行NVH試驗(yàn)并對(duì)數(shù)據(jù)進(jìn)行后處理,結(jié)合理論研究成果,對(duì)振動(dòng)噪聲源進(jìn)行定位,并提出相應(yīng)的改進(jìn)措施,為樣車的性能改進(jìn)提供方案與建議,為相似問題的理論與試驗(yàn)分析提供借鑒。

    1 樣車嘯叫問題

    1.1 駕駛室嘯叫噪聲分析

    本文所研究樣車為某重型商用車,該樣車在中低擋位即8擋及以下時(shí)駕駛室內(nèi)會(huì)出現(xiàn)一種比較明顯的類似口哨聲的嘯叫噪聲,且在臥鋪位置更加明顯,以7擋發(fā)動(dòng)機(jī)轉(zhuǎn)速在1 100~1 700 r/min時(shí)最為突出。當(dāng)擋位超過8擋后,嘯叫噪聲明顯降低??紤]到傳動(dòng)系統(tǒng)實(shí)際動(dòng)力傳遞路徑及嘯叫噪聲的產(chǎn)生原因,初步推測(cè)為變速箱或驅(qū)動(dòng)橋齒輪系統(tǒng)的嚙合激勵(lì)引起。

    1.2 樣車主要技術(shù)參數(shù)

    針對(duì)所研究樣車,其牽引車動(dòng)力總成采用直列六缸柴油發(fā)動(dòng)機(jī),并匹配16擋AMT手自一體變速箱,采用中、后雙驅(qū)動(dòng)橋和一段式傳動(dòng)軸形式。樣車的詳細(xì)技術(shù)參數(shù)如表1所示。

    針對(duì)樣車所搭載的AMT手自一體變速箱,其由主箱和副箱組成,主箱采用多級(jí)圓柱齒輪結(jié)構(gòu),副箱采用行星齒輪結(jié)構(gòu),共組成16個(gè)前進(jìn)擋和2個(gè)倒擋。主、副箱高低擋位變換均采用減速擋和直接擋組合的形式實(shí)現(xiàn),其中主箱采用接合套切換與變速箱前端兩組圓柱齒輪副的連接實(shí)現(xiàn);副箱采用接合套鎖止外齒圈或聯(lián)鎖外齒圈與行星架的方式實(shí)現(xiàn)。變速箱詳細(xì)擋位布置形式如圖1所示,考慮到實(shí)際需要,在結(jié)構(gòu)圖下方僅給出了7擋(4L擋)和8擋(4H擋)的傳力示意圖,另外,本文對(duì)2個(gè)倒擋不做研究。

    針對(duì)問題擋位,變速箱主箱對(duì)應(yīng)的7擋和8擋圓柱齒輪副齒數(shù)分配以及副箱外齒圈、太陽輪和行星輪齒數(shù)分配如圖2所示,另外,行星輪數(shù)量為5。中后驅(qū)動(dòng)橋各采用一對(duì)準(zhǔn)雙曲面錐齒輪,主動(dòng)輪和從動(dòng)輪齒數(shù)分別為14和38,中橋位置向后橋傳遞動(dòng)力的過渡圓柱齒輪副主被動(dòng)輪齒數(shù)均為38。

    2 整車道路試驗(yàn)

    2.1 試驗(yàn)?zāi)康?/h3>

    對(duì)樣車在整車加載條件下進(jìn)行道路試驗(yàn),采集轉(zhuǎn)速及振動(dòng)噪聲數(shù)據(jù),分析與定位駕駛室嘯叫噪聲產(chǎn)生的根源,為所研究車型NVH性能的改善提供參考和依據(jù)。加載試驗(yàn)采用牽引車拖拽掛車的方式進(jìn)行,加掛條件下整車質(zhì)量為49 t。

    2.2 傳感器布置

    試驗(yàn)設(shè)備主要包括LMS SCADAS數(shù)據(jù)采集前端、日本小野磁電式轉(zhuǎn)速傳感器、德國(guó)GRAS聲學(xué)傳感器和美國(guó)PCB三向加速度傳感器,數(shù)據(jù)采集軟件為與前端匹配的LMS Test.lab。所用三種類型傳感器的詳細(xì)布置情況如表2所示;包括型號(hào)、靈敏度和精度在內(nèi)的主要技術(shù)參數(shù)如表3所示。

    表2 各類型傳感器布置Tab.2 Arrangement of three types of sensors

    表3 傳感器主要技術(shù)參數(shù)Tab.3 Main technical parameters of three types of sensors

    數(shù)據(jù)采集前端及軟件如圖3(a)所示;圖3(b)為發(fā)動(dòng)機(jī)飛輪位置布置的轉(zhuǎn)速傳感器,其主要用于轉(zhuǎn)速追蹤,獲取振動(dòng)噪聲測(cè)量位置的階次特征。參考GB/T 18697—2002汽車車內(nèi)噪聲測(cè)量方法[19],如圖4所示,在駕駛員耳側(cè)、副駕耳側(cè)和臥鋪頭枕上方各布置一個(gè)聲學(xué)傳感器,另外在變速箱近場(chǎng)布置一個(gè)聲學(xué)傳感器以進(jìn)一步研究變速箱位置噪聲向駕駛室的傳遞。三向加速度傳感器的布置如圖5所示,布置位置包括變速箱殼體、中后橋殼體和車架上方。此處定義整車坐標(biāo)系方向?yàn)椋貉剀囕v前進(jìn)方向向后為正X方向,垂直于地面向上為正Z方向,正Y方向根據(jù)右手坐標(biāo)系確定。所用四個(gè)三向振動(dòng)加速度傳感器在整車坐標(biāo)系中的方向及歐拉角如表4所示。

    表4 加速度傳感器位置及歐拉角Tab.4 The direction and Euler angle of acceleration sensors

    2.3 試驗(yàn)工況

    道路試驗(yàn)在國(guó)內(nèi)某專業(yè)的汽車試驗(yàn)場(chǎng)進(jìn)行,試驗(yàn)條件滿足汽車道路試驗(yàn)方法通則[20]。試驗(yàn)中,采集樣車在7擋手動(dòng)模式下,發(fā)動(dòng)機(jī)轉(zhuǎn)速自550 r/min升至2 000 r/min的過程中,各測(cè)點(diǎn)處的振動(dòng)噪聲數(shù)據(jù)。變速箱在7擋工作模式下,其功率流向如圖6所示,動(dòng)力先后經(jīng)過主箱兩對(duì)圓柱齒輪副和副箱行星齒輪機(jī)構(gòu)。其中副箱鎖止外齒圈,功率從太陽輪輸入,行星架輸出。

    3 試驗(yàn)數(shù)據(jù)分析與噪聲源定位

    3.1 駕駛室內(nèi)噪聲

    以發(fā)動(dòng)機(jī)飛輪轉(zhuǎn)速為參考轉(zhuǎn)速進(jìn)行追蹤,即曲軸和變速箱主箱輸入軸的旋轉(zhuǎn)階次均為1階。在試驗(yàn)工況下測(cè)得的駕駛室內(nèi)主駕耳側(cè)、副駕耳側(cè)和臥鋪頭枕位置的噪聲結(jié)果圖7所示,圖中階次成分即表現(xiàn)出典型的嘯叫特性??梢钥闯?,除發(fā)動(dòng)機(jī)3階、6階以及附件旋轉(zhuǎn)軸系的一些激勵(lì)階次外,在三張Colormap圖中均出現(xiàn)了明顯的28階激勵(lì)成分。另外,噪聲峰值出現(xiàn)時(shí)對(duì)應(yīng)的發(fā)動(dòng)機(jī)轉(zhuǎn)速與前期主觀評(píng)價(jià)的結(jié)果基本一致。

    3.2 變速器近場(chǎng)噪聲及殼體振動(dòng)

    為進(jìn)一步分析與定位駕駛室內(nèi)嘯叫噪聲的來源,提取變速箱近場(chǎng)噪聲和殼體振動(dòng)的Colormap圖,其結(jié)果分別如圖8和圖9所示,其中殼體振動(dòng)結(jié)果僅取最為明顯的X方向。可以看出,在整體變化趨勢(shì)上,噪聲結(jié)果與駕駛室內(nèi)的評(píng)價(jià)結(jié)果基本相同,28階激勵(lì)最為明顯,同時(shí)還出現(xiàn)了一定的42階成分;在噪聲量級(jí)上,變速箱近場(chǎng)噪聲值明顯高于駕駛室內(nèi)部噪聲;另外,相比于其他階次成分,變速箱殼體28階激勵(lì)最為突出。

    3.3 中后橋殼體及車架振動(dòng)

    進(jìn)一步提取中后橋及車架上方X方向的振動(dòng)加速度,結(jié)果如圖10所示,可以看出,除發(fā)動(dòng)機(jī)激勵(lì)階次以及28階和42階激勵(lì)外,并未出現(xiàn)其他異常階次,并且從中橋向后橋的振動(dòng)明顯減弱。

    3.4 整體對(duì)比分析

    提取變速箱近場(chǎng)和駕駛室三個(gè)測(cè)點(diǎn)的噪聲Overall level曲線進(jìn)行對(duì)比分析,結(jié)果如圖11所示;同時(shí)提取變速箱殼體和中后橋及車架上方的振動(dòng)加速度Overall level曲線進(jìn)行對(duì)比分析,結(jié)果如圖12所示。

    綜合圖11和圖12可以看出,就噪聲而言,在變速箱近場(chǎng)噪聲向駕駛室內(nèi)傳遞的過程中,能量明顯減弱;就振動(dòng)而言,中后橋及車架位置的整體振動(dòng)水平明顯低于變速箱殼體。

    對(duì)噪聲和振動(dòng)最明顯的變速箱近場(chǎng)和殼體進(jìn)行切片分析,提取3階、6階、28階和42階切片,結(jié)果如圖13和圖14所示。初步來看,無論對(duì)于噪聲還是振動(dòng),28階激勵(lì)貢獻(xiàn)量最大,42階激勵(lì)的貢獻(xiàn)量也大于發(fā)動(dòng)機(jī)激勵(lì)。

    進(jìn)一步對(duì)0~2 000 r/min轉(zhuǎn)速范圍內(nèi)的噪聲和振動(dòng)水平進(jìn)行量化,將其均值以直方圖的形式進(jìn)行呈現(xiàn),結(jié)果如圖15和圖16所示。

    經(jīng)計(jì)算可知,28階激勵(lì)對(duì)總噪聲水平的貢獻(xiàn)量為85.17%,對(duì)總振動(dòng)水平的貢獻(xiàn)量為35.91%, 其對(duì)整體噪聲和振動(dòng)水平的貢獻(xiàn)量均大于其余階次。針對(duì)這種高階次激勵(lì),初步推測(cè)駕駛室嘯叫噪聲來源于變速箱齒輪嚙合。

    3.5 噪聲源定位與改善措施

    針對(duì)本樣車搭載的直列六缸發(fā)動(dòng)機(jī),3階和6階對(duì)應(yīng)發(fā)動(dòng)機(jī)激勵(lì)階次。根據(jù)圖6所示變速箱在7擋下的功率流走向并結(jié)合圖2中7擋對(duì)應(yīng)的齒輪副齒數(shù)分配可知,以發(fā)動(dòng)機(jī)曲軸旋轉(zhuǎn)階次為參考,主箱第一級(jí)齒輪副嚙合階次為34階,通過傳動(dòng)比計(jì)算,得主箱兩級(jí)齒輪機(jī)構(gòu)的傳動(dòng)比為1.21,主箱第二級(jí)齒輪副嚙合階次為31.4階。就副箱行星齒輪系統(tǒng)而言,相比于定軸輪系,行星齒輪不僅需要繞自身軸線自轉(zhuǎn),還要繞太陽輪公轉(zhuǎn),振動(dòng)激勵(lì)呈現(xiàn)多樣化,振動(dòng)特性較為復(fù)雜[21-24]。當(dāng)固定外齒圈,太陽輪驅(qū)動(dòng)行星輪及行星架轉(zhuǎn)動(dòng)時(shí),整個(gè)機(jī)構(gòu)相當(dāng)于一減速器,其傳動(dòng)比為:

    (1)

    式中:zring為外齒圈齒數(shù);zsun為太陽輪齒數(shù);iplanetary為固定外齒圈時(shí)行星齒輪機(jī)構(gòu)傳動(dòng)比;f0為太陽輪輸入轉(zhuǎn)頻;f1為行星架輸出轉(zhuǎn)頻。

    參考圖2副箱齒數(shù)分配,在式(1)中代入外齒圈和太陽輪齒數(shù),得iplanetary=4.32,于是在7擋工作模式下,變速箱傳動(dòng)比為5.23。行星架的旋轉(zhuǎn)階次,即中橋和后橋輸入軸的旋轉(zhuǎn)階次為0.19,錐齒輪嚙合階次為2.68階,中橋至后橋的過渡圓柱齒輪副的嚙合階次為7.30階。通過上述計(jì)算,并對(duì)比噪聲和振動(dòng)的試驗(yàn)結(jié)果,可以確定駕駛室嘯叫噪聲不是由變速箱主箱和中后橋齒輪副嚙合激勵(lì)引起。

    對(duì)于副箱行星齒輪機(jī)構(gòu),其在嚙合過程中會(huì)出現(xiàn)多個(gè)頻率成分,其定義和求解分別為

    (2)

    (3)

    (4)

    (5)

    (6)

    (7)

    (8)

    式中:n為行星輪個(gè)數(shù);zplanet為行星輪齒數(shù);[…]為分式取整;f2為齒圈任一個(gè)齒與行星輪的嚙合頻率;f3為太陽輪任一個(gè)齒與行星輪的嚙合頻率;f4為行星輪任一個(gè)齒與太陽輪或外齒圈的嚙合頻率;f5為單個(gè)行星輪與外齒圈的嚙合頻率;f6為所有行星輪與外齒圈的嚙合頻率;f7為因齒距誤差造成的單個(gè)行星輪與外齒圈的嚙合頻率。

    當(dāng)發(fā)動(dòng)機(jī)曲軸轉(zhuǎn)頻為1 Hz時(shí),結(jié)合式(1)~式(8),經(jīng)計(jì)算,f1~f7的值分別為0.19 Hz,0.96 Hz,3.19 Hz,0.56 Hz,14.02 Hz,70.08 Hz和13.44 Hz。所以,單個(gè)行星輪與齒圈的嚙合階次為14.02階,整車摸底試驗(yàn)中出現(xiàn)的28階和42階對(duì)應(yīng)單個(gè)行星輪與外齒圈嚙合階次的2次和3次諧階次;再者,根據(jù)噪聲與振動(dòng)的Colormap可以看出,諧階次能量明顯高于14.02階能量,這主要是由行星齒輪公轉(zhuǎn)軸線與外齒圈軸線出現(xiàn)偏置造成的[25]。

    針對(duì)本樣車駕駛室出現(xiàn)的嘯叫噪聲,通過整車道路摸底試驗(yàn),并經(jīng)過數(shù)據(jù)處理和理論分析,可以確定噪聲源來自變速箱副箱行星齒輪機(jī)構(gòu),是由行星輪與外齒圈的嚙合引起的。下一步可通過調(diào)整行星輪與外齒圈的宏觀幾何參數(shù),齒面微觀修形參數(shù)或改變行星輪齒數(shù)進(jìn)行優(yōu)化;并進(jìn)一步減小裝配誤差,調(diào)整行星輪和外齒圈軸線,消除軸線偏置造成的不良影響。

    4 結(jié) 論

    (1)本文通過整車加載道路試驗(yàn)獲取了樣車駕駛室內(nèi)噪聲、變速箱近場(chǎng)噪聲,變速箱殼體和中后橋殼體以及車架上方的振動(dòng)信號(hào)。對(duì)試驗(yàn)結(jié)果進(jìn)行了分析,各噪聲和振動(dòng)測(cè)點(diǎn)的Colormap圖和整體結(jié)果對(duì)比曲線表明,所測(cè)位置除發(fā)動(dòng)機(jī)激勵(lì)外,還受到28階和42階激勵(lì)的影響,且在變速箱近場(chǎng)和殼體位置最為突出。在噪聲從變速箱近場(chǎng)向駕駛室傳遞的過程中,能量明顯減弱;在中后橋及車架位置的振動(dòng)水平低于變速箱殼體,且未出現(xiàn)其他異常階次。

    (2)對(duì)變速箱主箱參與動(dòng)力傳遞的兩組嚙合齒輪副進(jìn)行了階次計(jì)算,對(duì)中后橋位置的兩組錐齒輪副進(jìn)行了階次計(jì)算,結(jié)果表明駕駛室內(nèi)的嘯叫噪聲不是這兩個(gè)位置的齒輪副嚙合激勵(lì)引起的。進(jìn)一步對(duì)變速箱副箱行星齒輪機(jī)構(gòu)進(jìn)行了轉(zhuǎn)頻及階次計(jì)算,結(jié)果表明,單個(gè)行星輪與外齒圈的嚙合階次為14.02階,其2次和3次諧階次分別為28階和42階,對(duì)應(yīng)各測(cè)點(diǎn)位置噪聲和振動(dòng)最明顯的階次成分,表明駕駛室嘯叫噪聲是由變速箱副箱行星齒輪機(jī)構(gòu)引起的。針對(duì)問題根源,提出了相應(yīng)的改善措施。

    未來的研究工作將針對(duì)引起駕駛室嘯叫噪聲的行星齒輪機(jī)構(gòu)進(jìn)行優(yōu)化,探究齒輪副宏觀設(shè)計(jì)參數(shù),齒面微觀修形參數(shù)以及行星輪個(gè)數(shù)對(duì)所研究問題的影響規(guī)律,并進(jìn)一步驗(yàn)證改進(jìn)措施的合理性。

    猜你喜歡
    振動(dòng)
    振動(dòng)的思考
    某調(diào)相機(jī)振動(dòng)異常診斷分析與處理
    振動(dòng)與頻率
    This “Singing Highway”plays music
    具非線性中立項(xiàng)的廣義Emden-Fowler微分方程的振動(dòng)性
    中立型Emden-Fowler微分方程的振動(dòng)性
    基于ANSYS的高速艇艉軸架軸系振動(dòng)響應(yīng)分析
    船海工程(2015年4期)2016-01-05 15:53:26
    主回路泵致聲振動(dòng)分析
    UF6振動(dòng)激發(fā)態(tài)分子的振動(dòng)-振動(dòng)馳豫
    帶有強(qiáng)迫項(xiàng)的高階差分方程解的振動(dòng)性
    涩涩av久久男人的天堂| 亚洲欧美一区二区三区黑人 | 中国国产av一级| 欧美国产精品va在线观看不卡| 三上悠亚av全集在线观看| 日韩成人伦理影院| 丝袜人妻中文字幕| 国国产精品蜜臀av免费| 美女国产高潮福利片在线看| 国产淫语在线视频| 夫妻性生交免费视频一级片| 亚洲美女视频黄频| 精品国产国语对白av| 晚上一个人看的免费电影| 宅男免费午夜| 亚洲精品自拍成人| 久久久久人妻精品一区果冻| 97精品久久久久久久久久精品| 男女午夜视频在线观看 | 亚洲欧美日韩卡通动漫| 飞空精品影院首页| 午夜福利,免费看| 久久免费观看电影| 亚洲国产欧美在线一区| 大香蕉97超碰在线| 精品视频人人做人人爽| 亚洲第一区二区三区不卡| 日本vs欧美在线观看视频| 国产1区2区3区精品| 美女视频免费永久观看网站| 尾随美女入室| 99香蕉大伊视频| 国产成人av激情在线播放| 国产精品99久久99久久久不卡 | 中文字幕制服av| 成人午夜精彩视频在线观看| 18在线观看网站| www.熟女人妻精品国产 | 美女主播在线视频| 一区二区三区四区激情视频| 三上悠亚av全集在线观看| 免费大片黄手机在线观看| 国产 一区精品| 免费大片18禁| 如何舔出高潮| 90打野战视频偷拍视频| 最新中文字幕久久久久| 18禁在线无遮挡免费观看视频| 久久久久久久国产电影| 国产av码专区亚洲av| 最新中文字幕久久久久| 精品卡一卡二卡四卡免费| 亚洲精品国产av蜜桃| 黄色一级大片看看| 国产成人91sexporn| 王馨瑶露胸无遮挡在线观看| 一个人免费看片子| 伊人亚洲综合成人网| 久久久国产一区二区| 在线 av 中文字幕| 赤兔流量卡办理| 日韩av不卡免费在线播放| 国产精品免费大片| 三级国产精品片| 18+在线观看网站| 亚洲av中文av极速乱| 久久久精品94久久精品| 久久久久久久大尺度免费视频| 亚洲精品美女久久久久99蜜臀 | 少妇精品久久久久久久| 国产成人精品久久久久久| 在线观看三级黄色| 午夜激情av网站| 人妻 亚洲 视频| 日韩av不卡免费在线播放| 狂野欧美激情性xxxx在线观看| 免费黄色在线免费观看| 亚洲国产精品999| 欧美日韩视频精品一区| 午夜91福利影院| 极品少妇高潮喷水抽搐| 高清黄色对白视频在线免费看| 久久久国产一区二区| 国产成人aa在线观看| 国产精品人妻久久久久久| 久久99精品国语久久久| 久久青草综合色| 狠狠精品人妻久久久久久综合| 性高湖久久久久久久久免费观看| 免费不卡的大黄色大毛片视频在线观看| 亚洲国产精品专区欧美| 久久午夜福利片| 99久久人妻综合| 99视频精品全部免费 在线| 18禁裸乳无遮挡动漫免费视频| 欧美 亚洲 国产 日韩一| 国产淫语在线视频| 久久99蜜桃精品久久| 午夜免费鲁丝| 国产乱人偷精品视频| 全区人妻精品视频| 多毛熟女@视频| 亚洲欧美日韩另类电影网站| 国产在线免费精品| 精品久久蜜臀av无| 人妻系列 视频| 亚洲在久久综合| 看非洲黑人一级黄片| 丝袜美足系列| 两个人免费观看高清视频| 汤姆久久久久久久影院中文字幕| 热re99久久国产66热| 久久久a久久爽久久v久久| 久久久国产精品麻豆| 日韩大片免费观看网站| 99热6这里只有精品| 欧美激情 高清一区二区三区| 中文欧美无线码| 少妇猛男粗大的猛烈进出视频| 亚洲精品一区蜜桃| a级毛片在线看网站| 免费av中文字幕在线| 高清欧美精品videossex| 亚洲精品国产av成人精品| 九色成人免费人妻av| 人人妻人人添人人爽欧美一区卜| 老司机亚洲免费影院| 乱人伦中国视频| 精品酒店卫生间| 国产一区二区三区av在线| 少妇人妻 视频| 成人亚洲精品一区在线观看| 在线看a的网站| 香蕉国产在线看| 成人国语在线视频| 黄色怎么调成土黄色| 天堂8中文在线网| 日本欧美国产在线视频| 精品亚洲乱码少妇综合久久| 午夜免费男女啪啪视频观看| 亚洲一区二区三区欧美精品| 免费av中文字幕在线| 午夜激情av网站| 丝袜喷水一区| 两个人免费观看高清视频| 国产免费福利视频在线观看| 精品久久久久久电影网| 久久99热这里只频精品6学生| 在线免费观看不下载黄p国产| 夜夜爽夜夜爽视频| 午夜视频国产福利| 免费女性裸体啪啪无遮挡网站| 亚洲色图 男人天堂 中文字幕 | 欧美日韩国产mv在线观看视频| 欧美日韩视频高清一区二区三区二| 丰满少妇做爰视频| 99久久精品国产国产毛片| 久久久久久人人人人人| 免费黄网站久久成人精品| 亚洲成色77777| 蜜臀久久99精品久久宅男| 国产片内射在线| 国产无遮挡羞羞视频在线观看| 99re6热这里在线精品视频| 少妇人妻 视频| 大香蕉97超碰在线| 男女无遮挡免费网站观看| 一级,二级,三级黄色视频| 在线观看美女被高潮喷水网站| 国产精品蜜桃在线观看| 久久久a久久爽久久v久久| 国产在线免费精品| 极品人妻少妇av视频| 亚洲四区av| 免费黄频网站在线观看国产| 久久久久久久久久成人| 香蕉国产在线看| 久久国产亚洲av麻豆专区| 高清毛片免费看| 成人免费观看视频高清| 亚洲熟女精品中文字幕| 国产精品女同一区二区软件| 国产日韩欧美视频二区| 人人妻人人爽人人添夜夜欢视频| 中文欧美无线码| 国产av国产精品国产| 狠狠精品人妻久久久久久综合| 成人亚洲精品一区在线观看| 免费播放大片免费观看视频在线观看| 欧美国产精品一级二级三级| 亚洲国产色片| 精品亚洲成国产av| 22中文网久久字幕| 在线观看国产h片| 韩国av在线不卡| 国产精品麻豆人妻色哟哟久久| 国产av精品麻豆| 亚洲av.av天堂| av国产久精品久网站免费入址| 丁香六月天网| 22中文网久久字幕| 校园人妻丝袜中文字幕| 国产精品久久久久久久久免| 亚洲丝袜综合中文字幕| 日本-黄色视频高清免费观看| 亚洲av.av天堂| 美女内射精品一级片tv| 丝袜在线中文字幕| av有码第一页| 欧美精品一区二区免费开放| av一本久久久久| 久久人人爽人人片av| 中国美白少妇内射xxxbb| 亚洲美女黄色视频免费看| 日韩伦理黄色片| 久久热在线av| 丰满饥渴人妻一区二区三| 久久av网站| 国产 一区精品| 一区二区三区四区激情视频| 成人毛片a级毛片在线播放| 精品一区二区三区四区五区乱码 | 国产片特级美女逼逼视频| 亚洲丝袜综合中文字幕| 高清视频免费观看一区二区| 国产 一区精品| 国精品久久久久久国模美| 日韩 亚洲 欧美在线| 欧美精品高潮呻吟av久久| 大陆偷拍与自拍| 久久久久久伊人网av| 亚洲精品国产av成人精品| 蜜桃国产av成人99| 在线天堂最新版资源| 国产亚洲午夜精品一区二区久久| 99热这里只有是精品在线观看| 国产探花极品一区二区| 26uuu在线亚洲综合色| 美女福利国产在线| 欧美亚洲日本最大视频资源| 捣出白浆h1v1| 国产黄色免费在线视频| 美女xxoo啪啪120秒动态图| 欧美亚洲日本最大视频资源| 波野结衣二区三区在线| 久久精品久久久久久噜噜老黄| 亚洲高清免费不卡视频| 狠狠精品人妻久久久久久综合| 国产有黄有色有爽视频| 精品久久久久久电影网| 亚洲精品国产av成人精品| 精品视频人人做人人爽| 亚洲精品视频女| 人人澡人人妻人| 97在线人人人人妻| 免费高清在线观看日韩| 黄片无遮挡物在线观看| 色视频在线一区二区三区| 老熟女久久久| 国产永久视频网站| 观看美女的网站| 晚上一个人看的免费电影| 免费高清在线观看日韩| 日韩三级伦理在线观看| 日本av手机在线免费观看| 久久久国产欧美日韩av| 日韩av不卡免费在线播放| xxxhd国产人妻xxx| av黄色大香蕉| 亚洲精品日韩在线中文字幕| 国产成人免费无遮挡视频| 日本91视频免费播放| 欧美精品一区二区免费开放| av国产久精品久网站免费入址| 久久国产精品男人的天堂亚洲 | 久久久久国产精品人妻一区二区| 久久久国产欧美日韩av| 亚洲av中文av极速乱| 免费大片18禁| videos熟女内射| 交换朋友夫妻互换小说| 最近中文字幕高清免费大全6| 久久精品国产a三级三级三级| 久久久久久久久久久免费av| 国产精品秋霞免费鲁丝片| 99久久精品国产国产毛片| 午夜免费男女啪啪视频观看| 少妇人妻 视频| 亚洲精品日本国产第一区| 国产色爽女视频免费观看| 日本-黄色视频高清免费观看| 九色成人免费人妻av| 少妇人妻精品综合一区二区| 亚洲丝袜综合中文字幕| 国产成人91sexporn| 五月天丁香电影| 午夜免费男女啪啪视频观看| av福利片在线| videossex国产| 亚洲人成网站在线观看播放| 久久精品久久久久久噜噜老黄| 一级毛片我不卡| 80岁老熟妇乱子伦牲交| 亚洲av成人精品一二三区| 狂野欧美激情性bbbbbb| 国产一区二区三区综合在线观看 | 欧美亚洲日本最大视频资源| 久久鲁丝午夜福利片| av片东京热男人的天堂| 青春草国产在线视频| 夫妻性生交免费视频一级片| 亚洲综合色网址| 飞空精品影院首页| 国产av码专区亚洲av| 精品久久蜜臀av无| 妹子高潮喷水视频| 国产有黄有色有爽视频| 搡老乐熟女国产| 99久久人妻综合| 视频在线观看一区二区三区| 久久97久久精品| 一个人免费看片子| 日产精品乱码卡一卡2卡三| 久久韩国三级中文字幕| 日韩电影二区| 午夜日本视频在线| 热re99久久国产66热| 一区二区日韩欧美中文字幕 | 少妇猛男粗大的猛烈进出视频| 另类亚洲欧美激情| 亚洲av.av天堂| 在线免费观看不下载黄p国产| 日韩熟女老妇一区二区性免费视频| 亚洲国产av新网站| 水蜜桃什么品种好| 五月开心婷婷网| 亚洲内射少妇av| 自线自在国产av| 日本黄大片高清| 欧美另类一区| 青春草国产在线视频| 美女脱内裤让男人舔精品视频| 精品第一国产精品| 亚洲少妇的诱惑av| 日本黄大片高清| 中文字幕最新亚洲高清| 国产成人欧美| 亚洲第一区二区三区不卡| 搡老乐熟女国产| 国产日韩欧美亚洲二区| 91国产中文字幕| 最后的刺客免费高清国语| 另类亚洲欧美激情| 日韩制服骚丝袜av| 国产精品99久久99久久久不卡 | 男人爽女人下面视频在线观看| 人妻一区二区av| 亚洲精品乱久久久久久| 美女视频免费永久观看网站| 中文字幕人妻丝袜制服| 国产日韩一区二区三区精品不卡| 国产在线一区二区三区精| 亚洲欧洲国产日韩| 久久精品国产a三级三级三级| 免费播放大片免费观看视频在线观看| 亚洲av在线观看美女高潮| 搡女人真爽免费视频火全软件| 熟妇人妻不卡中文字幕| 日韩不卡一区二区三区视频在线| 黑人巨大精品欧美一区二区蜜桃 | 日韩精品免费视频一区二区三区 | 精品人妻一区二区三区麻豆| 免费不卡的大黄色大毛片视频在线观看| 国产亚洲av片在线观看秒播厂| 99热网站在线观看| 国产日韩欧美在线精品| 精品久久国产蜜桃| a级片在线免费高清观看视频| 97在线视频观看| 少妇的逼水好多| 激情五月婷婷亚洲| 亚洲国产精品999| 国产日韩一区二区三区精品不卡| 不卡视频在线观看欧美| 日日摸夜夜添夜夜爱| 国产白丝娇喘喷水9色精品| 在线观看免费高清a一片| 久久久国产精品麻豆| 精品卡一卡二卡四卡免费| 黑丝袜美女国产一区| 国语对白做爰xxxⅹ性视频网站| 国产高清不卡午夜福利| 日本欧美国产在线视频| 天天躁夜夜躁狠狠久久av| 久久久久国产精品人妻一区二区| 免费人妻精品一区二区三区视频| 亚洲av综合色区一区| 午夜av观看不卡| 18禁动态无遮挡网站| 日韩精品有码人妻一区| 丝袜喷水一区| 国产精品三级大全| 宅男免费午夜| 看免费av毛片| 夫妻性生交免费视频一级片| 国内精品宾馆在线| 男女边摸边吃奶| 天堂俺去俺来也www色官网| 国产黄频视频在线观看| 欧美成人精品欧美一级黄| 欧美日韩精品成人综合77777| 天天躁夜夜躁狠狠躁躁| 亚洲一级一片aⅴ在线观看| 日本欧美视频一区| 欧美国产精品一级二级三级| 伦精品一区二区三区| 久久97久久精品| 国产成人精品在线电影| 夜夜骑夜夜射夜夜干| 狠狠婷婷综合久久久久久88av| 五月天丁香电影| 国产精品 国内视频| 99热6这里只有精品| 91aial.com中文字幕在线观看| 女人被躁到高潮嗷嗷叫费观| 久久久久久伊人网av| 午夜影院在线不卡| 91久久精品国产一区二区三区| 97超碰精品成人国产| 中国国产av一级| 捣出白浆h1v1| 美国免费a级毛片| 国产免费福利视频在线观看| 一级黄片播放器| 日本av免费视频播放| 中文字幕免费在线视频6| 99热网站在线观看| 秋霞伦理黄片| 久久久久久久精品精品| 国产精品一区二区在线观看99| 又黄又粗又硬又大视频| 美女大奶头黄色视频| 一边摸一边做爽爽视频免费| 少妇被粗大猛烈的视频| 国产日韩欧美亚洲二区| 多毛熟女@视频| 欧美人与善性xxx| av卡一久久| 免费观看av网站的网址| 99热6这里只有精品| av免费观看日本| 97在线人人人人妻| 大片电影免费在线观看免费| 国产成人一区二区在线| 9191精品国产免费久久| 22中文网久久字幕| 黄色 视频免费看| 精品国产一区二区久久| 精品视频人人做人人爽| 日韩成人av中文字幕在线观看| 少妇高潮的动态图| 老司机影院成人| 2022亚洲国产成人精品| 晚上一个人看的免费电影| 亚洲av.av天堂| 午夜福利乱码中文字幕| 最近中文字幕高清免费大全6| 日韩中文字幕视频在线看片| 久久久欧美国产精品| 色视频在线一区二区三区| 免费高清在线观看日韩| 黄色视频在线播放观看不卡| 欧美日韩视频精品一区| 免费黄频网站在线观看国产| 久久ye,这里只有精品| 捣出白浆h1v1| 亚洲四区av| 国产成人精品无人区| 性色av一级| 少妇人妻 视频| 亚洲三级黄色毛片| 久久热在线av| 国产av精品麻豆| 国产黄色视频一区二区在线观看| 最黄视频免费看| 一级爰片在线观看| 制服人妻中文乱码| 久久久久久久久久久免费av| 亚洲欧美中文字幕日韩二区| 五月伊人婷婷丁香| 日韩视频在线欧美| 日本爱情动作片www.在线观看| 国产熟女欧美一区二区| 秋霞在线观看毛片| 亚洲欧美成人综合另类久久久| 国产一区有黄有色的免费视频| 亚洲欧洲精品一区二区精品久久久 | 91成人精品电影| 亚洲欧洲日产国产| 久久鲁丝午夜福利片| 午夜91福利影院| 一级,二级,三级黄色视频| 三上悠亚av全集在线观看| av又黄又爽大尺度在线免费看| 久久亚洲国产成人精品v| 亚洲婷婷狠狠爱综合网| 亚洲激情五月婷婷啪啪| 丝瓜视频免费看黄片| 国产男女内射视频| 国产一区二区三区综合在线观看 | 免费人妻精品一区二区三区视频| 97在线人人人人妻| 久久久久人妻精品一区果冻| 人人妻人人澡人人爽人人夜夜| 国产老妇伦熟女老妇高清| 男人操女人黄网站| 国产一区有黄有色的免费视频| 两个人看的免费小视频| 国产熟女欧美一区二区| 丰满饥渴人妻一区二区三| 桃花免费在线播放| 久久av网站| 最近手机中文字幕大全| 街头女战士在线观看网站| 丁香六月天网| 国产综合精华液| 免费观看av网站的网址| 婷婷成人精品国产| 亚洲精品色激情综合| 色5月婷婷丁香| 婷婷色综合www| 成人手机av| 午夜福利视频在线观看免费| 在线观看美女被高潮喷水网站| 免费观看av网站的网址| 999精品在线视频| 一本—道久久a久久精品蜜桃钙片| 尾随美女入室| 2021少妇久久久久久久久久久| 午夜av观看不卡| 国产又爽黄色视频| 久久久久久久精品精品| 国产成人精品无人区| 中文字幕亚洲精品专区| 多毛熟女@视频| 国产成人精品在线电影| 99热全是精品| 一本色道久久久久久精品综合| 人人妻人人爽人人添夜夜欢视频| 久久国内精品自在自线图片| 制服诱惑二区| 成人综合一区亚洲| 91在线精品国自产拍蜜月| 天天操日日干夜夜撸| 国产乱来视频区| 十分钟在线观看高清视频www| 久久av网站| av不卡在线播放| 美女主播在线视频| 80岁老熟妇乱子伦牲交| 一级毛片我不卡| 国产一区二区三区av在线| 精品一区二区三卡| 美女福利国产在线| 亚洲精品久久久久久婷婷小说| 国产亚洲av片在线观看秒播厂| 在线观看美女被高潮喷水网站| 久久99一区二区三区| 成年av动漫网址| 狂野欧美激情性xxxx在线观看| 成人亚洲欧美一区二区av| 激情视频va一区二区三区| 精品午夜福利在线看| 精品人妻熟女毛片av久久网站| 国产精品久久久久久精品电影小说| 美女主播在线视频| 久久久欧美国产精品| 中文字幕最新亚洲高清| 国产精品久久久久久久电影| 香蕉精品网在线| 日韩av免费高清视频| 看十八女毛片水多多多| 青春草国产在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产一区二区激情短视频 | av视频免费观看在线观看| 亚洲国产最新在线播放| 91aial.com中文字幕在线观看| 亚洲国产看品久久| 91精品伊人久久大香线蕉| av在线观看视频网站免费| 久久国产精品男人的天堂亚洲 | 欧美精品一区二区大全| 黑人巨大精品欧美一区二区蜜桃 | 人体艺术视频欧美日本| 丝袜在线中文字幕| www.熟女人妻精品国产 | 国产精品国产三级国产专区5o| 精品久久蜜臀av无| videosex国产| 三级国产精品片| 五月玫瑰六月丁香| 亚洲经典国产精华液单| 黄网站色视频无遮挡免费观看| 又粗又硬又长又爽又黄的视频| 五月开心婷婷网| 建设人人有责人人尽责人人享有的| 国产一区亚洲一区在线观看| 国产白丝娇喘喷水9色精品| 丰满饥渴人妻一区二区三| 免费久久久久久久精品成人欧美视频 | 天天操日日干夜夜撸| 97在线人人人人妻| 91在线精品国自产拍蜜月| 男女边摸边吃奶|