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

    基于參數(shù)識(shí)別的截割頭縱向隨機(jī)振動(dòng)響應(yīng)優(yōu)化

    2015-10-29 03:37:53李曉豁沙永東
    中國(guó)機(jī)械工程 2015年6期
    關(guān)鍵詞:縱軸時(shí)變掘進(jìn)機(jī)

    李曉豁 何 洋 焦 麗 沙永東 李 婷

    遼寧工程技術(shù)大學(xué),阜新,123000

    基于參數(shù)識(shí)別的截割頭縱向隨機(jī)振動(dòng)響應(yīng)優(yōu)化

    李曉豁何洋焦麗沙永東李婷

    遼寧工程技術(shù)大學(xué),阜新,123000

    為減小縱軸式掘進(jìn)機(jī)截割頭的振動(dòng),利用Lagrange方法并結(jié)合虛擬激勵(lì)法推導(dǎo)了該型掘進(jìn)機(jī)的縱向運(yùn)動(dòng)微分方程,根據(jù)Daubechies小波函數(shù)多尺度逼近方法和最小二乘法原理,將掘進(jìn)機(jī)在隨機(jī)激勵(lì)作用下系統(tǒng)的時(shí)變阻尼和時(shí)變剛度用多尺度函數(shù)的線性組合表示,從而將時(shí)變參數(shù)識(shí)別問(wèn)題轉(zhuǎn)化為時(shí)不變問(wèn)題。將識(shí)別出的剛度和阻尼值作為初始條件,借助魚(yú)群算法以截割頭縱向振動(dòng)的位移響應(yīng)為目標(biāo)函數(shù)對(duì)參數(shù)進(jìn)行了優(yōu)化。結(jié)果表明:參數(shù)優(yōu)化后截割頭縱向振動(dòng)的位移響應(yīng)的均值下降了18.3%。為驗(yàn)證優(yōu)化結(jié)果的可靠性,利用ADAMS建立了整機(jī)的多剛體動(dòng)力學(xué)模型。仿真結(jié)果表明,截割頭的位移響應(yīng)變化與理論計(jì)算的結(jié)果較接近,從而驗(yàn)證了優(yōu)化結(jié)果的有效性。

    截割頭;虛擬激勵(lì)法;小波;最小二乘法;魚(yú)群算法

    0 引言

    截割頭是縱軸式掘進(jìn)機(jī)的工作機(jī)構(gòu),其作業(yè)環(huán)境惡劣,受煤巖賦存條件復(fù)雜、工作空間不斷變化、工作載荷多變等因素影響,截割頭在截割煤巖過(guò)程中會(huì)產(chǎn)生較大的振動(dòng),導(dǎo)致整機(jī)振動(dòng),這不僅容易損壞機(jī)器零部件、液壓及電氣元件,而且會(huì)降低機(jī)器工作可靠性和生產(chǎn)率。

    通過(guò)現(xiàn)場(chǎng)調(diào)研和實(shí)驗(yàn)統(tǒng)計(jì)發(fā)現(xiàn),縱軸式掘進(jìn)機(jī)的振動(dòng)形式有橫向(機(jī)身寬度方向)、縱向(機(jī)身高度方向)和軸向(機(jī)身長(zhǎng)度方向)三種,其中橫向和縱向的振動(dòng)較大,尤其是縱向振動(dòng)最明顯。

    目前,掘進(jìn)機(jī)振動(dòng)特性的研究主要有試驗(yàn)法、頻域分析法、載荷能量頻率識(shí)別法、求解振動(dòng)方程等方法。文獻(xiàn)[1-3]均采用試驗(yàn)分析方法在井下實(shí)測(cè)了整機(jī)振動(dòng)情況,得到了整機(jī)工作時(shí)的振動(dòng)頻率特征。童敏明等[4]使用遠(yuǎn)程監(jiān)控系統(tǒng),檢測(cè)掘進(jìn)機(jī)振動(dòng)特性,并進(jìn)行故障診斷和遠(yuǎn)程控制。張建廣[5]通過(guò)建立人工巖壁和試驗(yàn)測(cè)試系統(tǒng),并按照實(shí)際工況及不同工作參數(shù)對(duì)掘進(jìn)機(jī)的振動(dòng)數(shù)據(jù)從時(shí)域和頻域進(jìn)行了分析。本文第一作者等[6]采用小波包分析技術(shù)計(jì)算了重構(gòu)后的功率譜,檢測(cè)出振動(dòng)較大的頻率范圍。

    由于傳統(tǒng)方法計(jì)算隨機(jī)振動(dòng)需求解復(fù)雜的頻響函數(shù),計(jì)算量大、耗時(shí)多、精度低,加之隨機(jī)外載荷的反復(fù)作用導(dǎo)致系統(tǒng)的剛度和阻尼具有時(shí)變性,所以關(guān)于掘進(jìn)機(jī)系統(tǒng)的參數(shù)識(shí)別及隨機(jī)振動(dòng)的研究目前還鮮見(jiàn)報(bào)道。本文利用Lagrange方法并結(jié)合虛擬激勵(lì)法,建立縱軸式掘進(jìn)機(jī)的縱向運(yùn)動(dòng)微分方程,將平穩(wěn)隨機(jī)響應(yīng)轉(zhuǎn)化為簡(jiǎn)諧響應(yīng)進(jìn)行計(jì)算,采用小波包函數(shù)多尺度逼近的方法,將系統(tǒng)的時(shí)變剛度和時(shí)變阻尼用小波尺度函數(shù)的線性組合近似,將時(shí)變參數(shù)的識(shí)別問(wèn)題轉(zhuǎn)化為由已知的正交尺度函數(shù)和系統(tǒng)的輸入輸出來(lái)估計(jì)線性組合中的時(shí)不變系數(shù)問(wèn)題。

    由于截割頭縱向振動(dòng)較大,故將識(shí)別出的結(jié)果作為參數(shù)優(yōu)化的初始條件,采用魚(yú)群算法,以截割頭縱向振動(dòng)的位移響應(yīng)為目標(biāo)函數(shù)對(duì)系統(tǒng)的參數(shù)進(jìn)行優(yōu)化。

    1 小波參數(shù)識(shí)別原理

    1.1多尺度小波包分析

    設(shè){Vj}j∈Z為尺度空間L2(R)的一個(gè)多分辨率分析[7],相應(yīng)尺度函數(shù)為φ,小波函數(shù)為ψ,且滿足尺度方程:

    對(duì)于任意j∈Z,分解得

    Vj-1=Vj⊕Wj

    則尺度空間Vj與小波空間Wj的關(guān)系為

    Vj-1=Vj+1⊕Wj+1⊕Wj=Vj+2⊕Wj+2Wj=…

    (1)

    x(t)離散形式為

    (2)

    1.2時(shí)變參數(shù)識(shí)別方法

    單自由度時(shí)變結(jié)構(gòu)系統(tǒng)的運(yùn)動(dòng)方程[8]為

    (3)

    式中,m為質(zhì)量;c(t)、k(t)分別為時(shí)變阻尼和時(shí)變剛度。

    式(3)的離散形式為

    (4)

    將時(shí)變阻尼和時(shí)變剛度看作離散時(shí)間序列信號(hào),設(shè)J層小波分解的小波系數(shù)和尺度系數(shù)已知,根據(jù)式(2)將其展開(kāi),對(duì)于慢變信號(hào),能量大部分集中在低頻部分,展開(kāi)時(shí)可忽略第2項(xiàng)細(xì)節(jié)信號(hào),僅由第1項(xiàng)概貌信號(hào)來(lái)近似表示:

    (5)

    將式(5)代入式(4)得

    將離散時(shí)刻n=1,2,…,N的響應(yīng)代入上式,有

    g×q=r

    g=[g(c)g(k)]T

    q=[q(c)q(k)]T

    由最小二乘法求出尺度系數(shù)和小波展開(kāi)系數(shù):

    q=(gTg)-1gTr

    將上式求出的結(jié)果代入式(5),可求出系統(tǒng)的時(shí)變剛度k(n)=q(k)和時(shí)變阻尼c(n)=q(c)。

    2 掘進(jìn)機(jī)動(dòng)力學(xué)模型

    2.1基本假設(shè)

    縱軸式掘進(jìn)機(jī)是結(jié)構(gòu)復(fù)雜的大型工礦設(shè)備,為求解方便,將其進(jìn)行適當(dāng)?shù)暮?jiǎn)化和假設(shè)[9]:根據(jù)組成將其簡(jiǎn)化為三自由度的動(dòng)力學(xué)模型,截割頭、懸臂和機(jī)體(包括行走機(jī)構(gòu))的質(zhì)量分別為m1、m2、m3,煤壁單軸壓縮強(qiáng)度恒定,各部分通過(guò)無(wú)質(zhì)量的彈性元件連接,截割頭與懸臂、懸臂與機(jī)體、機(jī)體與底板間的時(shí)變剛度和時(shí)變阻尼分別為k1(t)、k2(t)、k3(t)和c1(t)、c2(t)、c3(t),系統(tǒng)的剛度和阻尼不發(fā)生突變,隨時(shí)間緩慢變化。由此建立縱軸式掘進(jìn)機(jī)的縱向(垂直底板,y向)動(dòng)力學(xué)模型,如圖1所示。

    圖1 掘進(jìn)機(jī)三自由縱向動(dòng)力學(xué)模型

    Lagrange方程為

    (6)

    式中,qj為系統(tǒng)廣義坐標(biāo);Qj為廣義坐標(biāo)對(duì)應(yīng)的廣義力。

    系統(tǒng)的勢(shì)能V、動(dòng)能T和耗散函數(shù)D為

    (7)

    將式(7)代入式(6),則掘進(jìn)機(jī)縱向運(yùn)動(dòng)微分方程為

    (8)

    2.2隨機(jī)激勵(lì)作用下系統(tǒng)的運(yùn)動(dòng)微分方程

    線性系統(tǒng)受到自功率譜密度為Sxx的平穩(wěn)隨機(jī)激勵(lì)x(t)作用時(shí),其響應(yīng)的自功率譜Syy=|H|2Sxx。常規(guī)算法求系統(tǒng)響應(yīng)需要大量時(shí)間求解頻響函數(shù)矩陣H,計(jì)算較繁瑣。

    由虛擬激勵(lì)法原理可知[10-11],當(dāng)隨機(jī)激勵(lì)被簡(jiǎn)諧激勵(lì)ei ω t代替時(shí),相應(yīng)的簡(jiǎn)諧響應(yīng)y=Hei ω t,若在激勵(lì)前乘以(Sx x)1/2,顯然,響應(yīng)譜為

    (9)

    設(shè)線性系統(tǒng)受到隨機(jī)平穩(wěn)激勵(lì)x(t)的作用,其運(yùn)動(dòng)方程為

    (10)

    式中,x(t)為功率譜Sxx的平穩(wěn)隨機(jī)過(guò)程;P為給定常數(shù)向量;M為質(zhì)量矩陣;K、C分別為剛度矩陣和阻尼矩陣。

    統(tǒng)計(jì)分析表明[12],截割頭載荷可視為平穩(wěn)隨機(jī)過(guò)程,截齒受力可按照瑞利分布的隨機(jī)過(guò)程來(lái)考慮,截割頭載荷是各截齒受力的疊加,其縱向隨機(jī)激勵(lì)自功率譜密度序列為

    式中,Ry(n)為截?cái)嚯S機(jī)信號(hào)的自相關(guān)函數(shù);ω=2πk/N。

    (11)

    P=[100]T

    將剛度和阻尼作為未知量,把相同剛度和阻尼的系數(shù)移到同一邊,則式(11)可改寫(xiě)為

    (12)

    根據(jù)上述的識(shí)別方法,將式(12)展開(kāi)成矩陣形式:

    G×H=R

    (13)

    對(duì)式(12)中待求的時(shí)變剛度k(t)、時(shí)變阻尼c(t),利用式(13)展開(kāi)有

    G=[G(c)G(k)]H=[H(c)H(k)]T

    G(c)和G(k)矩陣關(guān)系式見(jiàn)表1。

    表1 關(guān)系式系數(shù)

    注:i=1,2,3;離散時(shí)刻n=1,2,…,N。

    根據(jù)最小二乘法原理H=(GTG)-1R,識(shí)別出系統(tǒng)的時(shí)變阻尼和時(shí)變剛度分別為H(c)、H(k)。

    2.3識(shí)別結(jié)果

    以國(guó)產(chǎn)某型縱軸式掘進(jìn)機(jī)為研究對(duì)象,其截割頭、懸臂和機(jī)體的質(zhì)量分別為m1=1051 kg,m2=2415 kg,m3=56 534 kg;采樣頻率fS=512 Hz,假設(shè)隨機(jī)外載荷譜為服從瑞利分布的高斯白噪聲譜。

    當(dāng)信號(hào)中存在噪聲時(shí),識(shí)別方程將發(fā)生病態(tài),這里采用吉洪諾夫正則化方法來(lái)減小識(shí)別方程的不適定性對(duì)識(shí)別結(jié)果的影響,提高識(shí)別精度。

    將變剛度和變阻尼分別用db4小波在尺度3上的正交尺度函數(shù)的線性組合表示,用最小二乘法識(shí)別的剛度和阻尼值如圖2和圖3所示。

    (a)截割頭剛度k1

    (b)懸臂剛度k2

    (c)機(jī)體剛度k31.理論值 2.5%噪聲圖2 剛度識(shí)別結(jié)果

    由圖2和圖3可見(jiàn),與理論值相比,在噪聲的影響下,系統(tǒng)的剛度和阻尼的識(shí)別結(jié)果具有一定的波動(dòng)性,由于系統(tǒng)的隨機(jī)外載荷譜采用了假設(shè)條件,因此對(duì)計(jì)算精度會(huì)有一定影響,但其曲線始終在理論值的周圍波動(dòng),且保持了較好的跟蹤效果。在噪聲影響下,相對(duì)于剛度識(shí)別結(jié)果,阻尼識(shí)別值與理論值之間的波動(dòng)相對(duì)較大,剛度識(shí)別的效果好于阻尼,這是由于剛度值比阻尼值大很多,同時(shí)計(jì)算得出的相對(duì)誤差阻尼要更大一些,可見(jiàn),與剛度的識(shí)別結(jié)果相比,阻尼的識(shí)別對(duì)噪聲較敏感,誤差較大。

    (a)截割頭阻尼c1

    (b)懸臂阻尼c2

    (c)機(jī)體阻尼c31.理論值 2.5%噪聲圖3 阻尼的識(shí)別結(jié)果

    3 響應(yīng)優(yōu)化

    3.1目標(biāo)函數(shù)

    以截割頭縱向振動(dòng)位移響應(yīng)的均方值為目標(biāo)函數(shù),利用魚(yú)群算法對(duì)系統(tǒng)的剛度和阻尼進(jìn)行優(yōu)化,以減小截割頭的振動(dòng)。令

    (14)

    將式(14)代入式(11),解得

    利用系數(shù)矩陣的對(duì)稱性,求得截割頭、懸臂和機(jī)體的縱向位移響應(yīng):

    式中,Λ*為Λ的伴隨矩陣。

    進(jìn)而得到截割頭(Syy1)、懸臂(Syy2)和機(jī)體(Syy3)的縱向位移響應(yīng)功率譜矩陣:

    Syy1=B1B1

    Syy2=B2B2

    Syy3=B3B3

    當(dāng)截割頭受到縱向隨機(jī)力Qy的作用時(shí),其響應(yīng)可由Duhamel積分表示為

    式中,h(t)為脈沖函數(shù)。

    其響應(yīng)的相關(guān)函數(shù)為

    Ryy(t1,t2)=E[y(t1)y(t2)]=

    (15)

    由維納-辛欽定理有

    (16)

    h(t)與頻率響應(yīng)函數(shù)H(ω)構(gòu)成傅里葉變換對(duì):

    (17)

    將式(16)、式(17)代入式(13),則響應(yīng)的均方值為

    (18)

    將式(9)代入式(18),截割頭的縱向振動(dòng)位移響應(yīng)響應(yīng)均方值可表示為

    即目標(biāo)函數(shù)為

    (19)

    3.2約束條件

    將k1、k2、k3、c1、c2、c3作為優(yōu)化變量,其取值范圍如表2所示。

    表2 優(yōu)化變量取值范圍

    懸臂的最大撓度Ymax不應(yīng)超過(guò)材料許用值[Y],即

    Ymax=ql4/(8EI)≤[Y]

    式中,q為集中載荷,N;E為材料的彈性模量,Pa;I為慣性矩,kg·m2;l為懸臂長(zhǎng)度,m。

    懸臂承受的最大剪應(yīng)力τmax不應(yīng)超過(guò)材料的極限應(yīng)力[τ],即

    τmax=T/Wn≤[τ]

    式中,T為扭矩,N·m;Wn為抗扭截面系數(shù)。

    3.3優(yōu)化計(jì)算

    借助MATLAB匯編語(yǔ)言編制魚(yú)群算法程序,以式(19)為目標(biāo)函數(shù)對(duì)參數(shù)進(jìn)行優(yōu)化。相關(guān)參數(shù)為:截割頭轉(zhuǎn)速nJ=46 r/min,截齒數(shù)z=48,煤巖崩落角Φ=45°,煤巖抗截強(qiáng)度A=378 N/mm,煤巖脆性度B=2.5,l=1.9 m,轉(zhuǎn)矩Tm=580 N·m。取人工魚(yú)群數(shù)50,感知距離1,最大迭代次數(shù)50,擁擠度因子0.618,移動(dòng)步長(zhǎng)0.1,覓食最大試探次數(shù)100。經(jīng)優(yōu)化,得到系統(tǒng)的剛度和阻尼參數(shù)值,見(jiàn)表3。優(yōu)化前后的截割頭位移響應(yīng)如圖4所示。

    表3 優(yōu)化結(jié)果

    圖4 截割頭位移響應(yīng)

    由表3可見(jiàn),k3和c3的變化量最大,這表明機(jī)體與底板間剛度和阻尼的增加可有效減小截割頭的振動(dòng)。通過(guò)對(duì)截割頭位移響應(yīng)的計(jì)算(圖4)知,優(yōu)化前后截割頭的位移響應(yīng)的均值分別為9.3 mm、7.6 mm,均值減小了18.3%。

    4 結(jié)論的驗(yàn)證

    4.1功率譜分析

    由識(shí)別出的參數(shù),計(jì)算得到在平穩(wěn)隨機(jī)激勵(lì)作用下該機(jī)截割頭、懸臂和機(jī)體的縱向位移功率譜密度仿真曲線,如圖5所示。圖5中的4個(gè)峰值集中在低頻2~5 Hz,對(duì)應(yīng)機(jī)體(2.8 Hz)、懸臂(3.4 Hz)以及截割頭(4.2 Hz)的固有頻率,這與文獻(xiàn)[12-13]的結(jié)論及實(shí)驗(yàn)結(jié)果基本吻合,從而驗(yàn)證了識(shí)別結(jié)果的可靠性。

    圖5 位移功率譜密度

    4.2模型仿真

    根據(jù)多剛體動(dòng)力學(xué)原理,利用Pro/E建立縱軸式掘進(jìn)機(jī)三維模型,將掘進(jìn)機(jī)整機(jī)模型映射到多體動(dòng)力學(xué)軟件ADAMS中,定義材料屬性、載荷以及構(gòu)件間的接觸力、摩擦力和阻尼參數(shù)等,并對(duì)模型進(jìn)行適當(dāng)?shù)暮?jiǎn)化,完成的整機(jī)模型如圖6所示。

    圖6 掘進(jìn)機(jī)多剛體動(dòng)力學(xué)模型

    根據(jù)建立的掘進(jìn)機(jī)多剛體動(dòng)力學(xué)模型,取步頻300 Hz、仿真時(shí)間10 s,對(duì)識(shí)別參數(shù)值與優(yōu)化后參數(shù)值分別進(jìn)行仿真,求得截割頭的縱向振動(dòng)特征如圖7所示。

    圖7 模型仿真的截割頭位移響應(yīng)

    統(tǒng)計(jì)表明,優(yōu)化前后截割頭的縱向位移響應(yīng)均值分別為11.3 mm、9.8 mm,優(yōu)化后均值減小了13.2%,模型仿真結(jié)果與理論計(jì)算結(jié)果相比,均值優(yōu)化量?jī)H相差5.1%,滿足精度要求。

    5 結(jié)論

    (1)利用Lagrange方法并結(jié)合虛擬激勵(lì)法,建立了縱軸式掘進(jìn)機(jī)隨機(jī)激勵(lì)下的運(yùn)動(dòng)微分方程。

    (2)根據(jù)Daubechies小波和最小二乘法原理,對(duì)系統(tǒng)的時(shí)變剛度和時(shí)變阻尼進(jìn)行了識(shí)別,并以識(shí)別值為初始條件對(duì)參數(shù)進(jìn)行了優(yōu)化,優(yōu)化后截割頭振動(dòng)的位移響應(yīng)均值降低了18.3%,有效減小了截割頭的振動(dòng)。

    (3)根據(jù)多剛體動(dòng)力學(xué)原理,利用ADAMS建立掘進(jìn)機(jī)的多剛體動(dòng)力學(xué)模型,模型仿真結(jié)果與理論計(jì)算值較接近,驗(yàn)證了優(yōu)化結(jié)果的可靠性。

    [1]Kr?ger G,Brzchta P.Neue Meβund Diagnose System für F?rdererund Fobelanlagen[J].Glückauf,1991,127(17/18):778-785.

    [2]黃民,吳淼,安偉,等.橫切割頭掘進(jìn)機(jī)械振動(dòng)特性研究[J].機(jī)械工程學(xué)報(bào),2002,38(8):89-93.

    Huang Min,Wu Miao,An Wei,et al.Research on Vibration Characteristics of Horizontal Axis Tunneller[J].Chinese Journal of Mechanical Engineering,2002,38(8):89-93.

    [3]Tao Alvarez J,Menendez M.Rodriguez Diez R,et al.Experimental Results of a Low-power Roadheader Driving a Gallery with Different Types of Rock at the Face[J].Tunnelling and Underground Space Technology,2003,18(4):395-404.

    [4]童敏明,吳軻.基于振動(dòng)傳感器的掘進(jìn)機(jī)故障監(jiān)測(cè)系統(tǒng)[J].煤礦機(jī)械,2010,31(5):241-243.

    Tong Minming,Wu Ke.Based on Vibration Sensor Tunnel Boring Machine Fault Monitoring System[J].Coal Mine Machinery,2010,31(5):241-243.

    [5]張建廣.縱軸式掘進(jìn)機(jī)截割人工巖壁的振動(dòng)試驗(yàn)研究[J].煤炭科學(xué)技術(shù),2011,39(8):76-78.

    Zhang Jianguang.Study on Vibration Experiment of Artificial Rock Wall Cutting with Longitudinal Shaft Type Mine Roadheader[J].Coal Science and Technology,2011,39(8):76-78.

    [6]李曉豁,韓宇飛.掘進(jìn)機(jī)截割載荷極值能量的頻率識(shí)別[J].煤炭學(xué)報(bào),2008,33(4):459-461.

    Li Xiaohuo,Han Yufei.Frequency Recognition of Cutting Loads Extreme Energy for Roadheader[J]. Journal of China Coal Society, 2008,33(4):459-461.

    [7]潘學(xué)萍,袁姍姍,鞠萍,等.負(fù)荷隨機(jī)擾動(dòng)下的電力系統(tǒng)小波模態(tài)參數(shù)識(shí)別[J].電力自動(dòng)化設(shè)備,2012, 32(4):6-10.

    Pan Xueping, Yuan Shanshan, Ju Ping, et al. Wavelet Based Electromechanical Mode Shape Identification From Ambiet Data of Random Load Excitayion[J]. Electric Power Automation Equipment,2012,32(4):6-10.

    [8]潘學(xué)萍,扈衛(wèi)衛(wèi),尚霏.多信號(hào)模態(tài)參數(shù)識(shí)別的小波方法[J].電力自動(dòng)化設(shè)備,2013,33(5):31-36.

    Pan Xueping,Hu Weiwei,Shang Fei.Criterion of Current Differential Protection Based on Amplitude and Phase Difference of Fault Current[J].Electric Power Automation Equipment,2013,33(5):31-36.

    [9]居乃鵕.裝甲車輛動(dòng)力學(xué)分析與仿真[M].北京:國(guó)防工業(yè)出版社,2002.

    [10]李杰,秦玉英,趙旗.汽車行駛平順性建模與仿真的新方法[J]. 中國(guó)機(jī)械工程,2009,20(13):1264-1269.

    Li Jie,Qin Yuying,Zhao Qi.A New Mathematical Modelling and Simulation Method for Vehicle Ride Comfort[J].China Mechanical Engineering,2009,20(13):1264-1269.

    [11]李曉豁.掘進(jìn)機(jī)截割頭隨機(jī)載荷的模擬研究[J].煤炭學(xué)報(bào),2000,25(5):73-76.

    Li Xiaohuo.Simulation Study of Random Loads on a Cutting Head of Roadheader[J].Journal of China Coal Society,2000,25(5):73-76.

    [12]李曉豁.掘進(jìn)機(jī)截割頭的關(guān)鍵技術(shù)研究[M].北京:機(jī)械工業(yè)出版社,2008.

    [13]李曉豁,韓宇飛.縱軸式掘進(jìn)機(jī)截割的功率譜分析[J].中國(guó)工程機(jī)械學(xué)報(bào),2007,5(2):134-137.

    Li Xiaohuo,Han Yufei.Spectral Analysis of Cutting Power for Longitudinal Excavators[J].Chinese Journal of Construction Machinery,2007,5(2):134-137.

    (編輯袁興玲)

    Vertical Random Vibration Response and Optimization of Cutting Head Based on Parameter Identification

    Li XiaohuoHe YangJiao LiSha YongdongLi Ting

    Liaoning Technical University,Fuxin,Liaoning,123000

    In order to reduce the vibration of a longitudinal roadheader’s cutting head,differential equations of motion of the roadheader were derived by using Lagrange equation and pseudo-excitation method.By expanding the time-varying stiffness and damping of the system into linear combination of a series scaling functions based on Daubechies wavelet’s compactness and regularization and least-square method,the time-varying parameter identification problem was turned into the problem of identifying invariant parameters.Taking the identifying parameters as the initial conditions,choosing displacement response of the cutting head as an optimal objective and parameters were optimized by means of fish-swarm algorithm.The results show that the displacement response of the cutting head is decreased by 18.3% after optimized.A multi-body dynamic model of the roadheader was simulated in ADAMS for the sake of identifying the conclusion,the model simulation results and the theoretical value are relatively close,which implies validity of the optimized conclusion.

    cutting head;pseudo-excitation method;wavelet;least-square method;fish-swarm algorithm

    2013-07-30

    2014-12-09

    國(guó)家自然科學(xué)基金資助項(xiàng)目(59774033);中國(guó)煤炭工業(yè)協(xié)會(huì)科學(xué)技術(shù)研究指導(dǎo)性計(jì)劃資助項(xiàng)目(MTKJ2011-325);遼寧省大型工礦裝備重點(diǎn)實(shí)驗(yàn)室資助項(xiàng)目(2008403010)

    TD421.5DOI:10.3969/j.issn.1004-132X.2015.06.021

    李曉豁,男,1953年生。遼寧工程技術(shù)大學(xué)機(jī)械工程學(xué)院教授、博士研究生導(dǎo)師。主要研究方向?yàn)楝F(xiàn)代機(jī)械設(shè)計(jì)理論與方法、機(jī)電液系統(tǒng)仿真與應(yīng)用。出版專著3部,發(fā)表論文396篇。何洋,男,1982年生。遼寧工程技術(shù)大學(xué)機(jī)械工程學(xué)院博士研究生。焦麗,女,1979年生。遼寧工程技術(shù)大學(xué)機(jī)械工程學(xué)院講師、博士。沙永東,男,1980年生。遼寧工程技術(shù)大學(xué)機(jī)械工程學(xué)院講師、博士。李婷,女,1984年生。遼寧工程技術(shù)大學(xué)機(jī)械工程學(xué)院博士研究生。

    猜你喜歡
    縱軸時(shí)變掘進(jìn)機(jī)
    隧道掘進(jìn)機(jī)
    掘進(jìn)機(jī)掘進(jìn)速度模糊PID控制
    基于時(shí)變Copula的股票市場(chǎng)相關(guān)性分析
    流動(dòng)管理:一種班級(jí)管理的實(shí)踐模式
    縱軸曝氣氧化溝在城鎮(zhèn)污水處理中的應(yīng)用
    煙氣輪機(jī)復(fù)合故障時(shí)變退化特征提取
    從時(shí)間縱軸研究大學(xué)生身體素質(zhì)狀況及對(duì)策
    基于MEP法的在役橋梁時(shí)變可靠度研究
    關(guān)于縱軸式掘進(jìn)機(jī)截割頭主要參數(shù)的研究
    國(guó)內(nèi)首臺(tái)全斷面煤巷高效掘進(jìn)機(jī)誕生并交付
    久久草成人影院| 色尼玛亚洲综合影院| 色播在线永久视频| 国产精品综合久久久久久久免费 | 真人做人爱边吃奶动态| 黄网站色视频无遮挡免费观看| 免费黄频网站在线观看国产| 亚洲九九香蕉| 一进一出抽搐动态| 久久久久精品国产欧美久久久| 老司机深夜福利视频在线观看| 好男人电影高清在线观看| 国产av一区二区精品久久| 精品国产乱码久久久久久男人| 久99久视频精品免费| 久久午夜综合久久蜜桃| 久久国产精品影院| 国产高清激情床上av| 91av网站免费观看| 777米奇影视久久| 久久精品国产综合久久久| 超碰97精品在线观看| 国产91精品成人一区二区三区| 大型av网站在线播放| 国产成人免费观看mmmm| 国产在视频线精品| 午夜福利视频在线观看免费| 国产亚洲精品久久久久5区| 黄色片一级片一级黄色片| 久久人人爽av亚洲精品天堂| 亚洲欧美激情在线| 日本精品一区二区三区蜜桃| 丝袜人妻中文字幕| 人人澡人人妻人| 女同久久另类99精品国产91| 国产精品国产av在线观看| 国产成人精品久久二区二区91| 韩国精品一区二区三区| 一本一本久久a久久精品综合妖精| 夜夜爽天天搞| 欧洲精品卡2卡3卡4卡5卡区| 黄色成人免费大全| 亚洲视频免费观看视频| 国产精华一区二区三区| 黄色a级毛片大全视频| 欧美不卡视频在线免费观看 | 精品久久久久久久毛片微露脸| 99热国产这里只有精品6| 我的亚洲天堂| 午夜两性在线视频| 一本一本久久a久久精品综合妖精| 18在线观看网站| 高清在线国产一区| 精品国产亚洲在线| 亚洲成av片中文字幕在线观看| 这个男人来自地球电影免费观看| 久久热在线av| 国产精品亚洲一级av第二区| 宅男免费午夜| 成人黄色视频免费在线看| 黑人欧美特级aaaaaa片| 国产亚洲精品久久久久久毛片 | 免费一级毛片在线播放高清视频 | 女同久久另类99精品国产91| 一本一本久久a久久精品综合妖精| 在线观看免费日韩欧美大片| 搡老乐熟女国产| 在线观看66精品国产| 日韩欧美一区二区三区在线观看 | 成熟少妇高潮喷水视频| 亚洲精品久久午夜乱码| 欧美日韩瑟瑟在线播放| 一进一出抽搐gif免费好疼 | 深夜精品福利| 91av网站免费观看| 日韩人妻精品一区2区三区| 黄片大片在线免费观看| 十分钟在线观看高清视频www| 成熟少妇高潮喷水视频| 男男h啪啪无遮挡| 不卡av一区二区三区| 亚洲欧美激情在线| 久久狼人影院| 天天添夜夜摸| 欧美人与性动交α欧美精品济南到| 性少妇av在线| 国产高清国产精品国产三级| av电影中文网址| 欧美精品亚洲一区二区| 成年动漫av网址| 亚洲av第一区精品v没综合| 久久精品国产99精品国产亚洲性色 | 亚洲欧美激情综合另类| 成人国语在线视频| 日韩欧美在线二视频 | 亚洲成人免费av在线播放| 欧美午夜高清在线| 少妇裸体淫交视频免费看高清 | 精品久久久久久,| 一区二区日韩欧美中文字幕| а√天堂www在线а√下载 | 亚洲欧美色中文字幕在线| 日韩精品免费视频一区二区三区| www.999成人在线观看| 欧美日韩成人在线一区二区| 日本一区二区免费在线视频| 很黄的视频免费| 这个男人来自地球电影免费观看| 欧美日韩成人在线一区二区| 嫁个100分男人电影在线观看| 18禁美女被吸乳视频| 18禁裸乳无遮挡免费网站照片 | 国产欧美亚洲国产| 国内毛片毛片毛片毛片毛片| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美丝袜亚洲另类 | 91九色精品人成在线观看| 9热在线视频观看99| 免费在线观看影片大全网站| 国产不卡一卡二| 啦啦啦 在线观看视频| 90打野战视频偷拍视频| 国产在线一区二区三区精| tube8黄色片| 国产欧美日韩一区二区三| 国产精品影院久久| 日韩中文字幕欧美一区二区| 一本综合久久免费| 黄色a级毛片大全视频| www.熟女人妻精品国产| 国产精品永久免费网站| 国产男女内射视频| 精品欧美一区二区三区在线| 人妻久久中文字幕网| 国产单亲对白刺激| 久久精品国产99精品国产亚洲性色 | 18禁黄网站禁片午夜丰满| 免费一级毛片在线播放高清视频 | 国产精品久久久久久精品古装| 女性生殖器流出的白浆| av网站在线播放免费| 9热在线视频观看99| 欧美午夜高清在线| 性少妇av在线| 一边摸一边做爽爽视频免费| videosex国产| 黑人猛操日本美女一级片| 久久香蕉精品热| 91精品三级在线观看| 中文字幕人妻熟女乱码| 亚洲va日本ⅴa欧美va伊人久久| av在线播放免费不卡| 一边摸一边抽搐一进一小说 | 俄罗斯特黄特色一大片| 国产日韩一区二区三区精品不卡| 午夜免费成人在线视频| 一级,二级,三级黄色视频| 成人18禁高潮啪啪吃奶动态图| 悠悠久久av| 国产精品av久久久久免费| 熟女少妇亚洲综合色aaa.| 亚洲五月色婷婷综合| 欧美黑人欧美精品刺激| 夜夜夜夜夜久久久久| 老熟妇乱子伦视频在线观看| 天堂√8在线中文| 久久久精品国产亚洲av高清涩受| 午夜激情av网站| 天堂动漫精品| 久久草成人影院| 久久国产精品影院| 成在线人永久免费视频| 国产片内射在线| 国产一区二区三区视频了| 久久久国产精品麻豆| 婷婷丁香在线五月| 午夜视频精品福利| 国产成人精品久久二区二区91| 久久影院123| 嫩草影视91久久| 性色av乱码一区二区三区2| 亚洲综合色网址| 激情视频va一区二区三区| 亚洲精品粉嫩美女一区| 日韩欧美在线二视频 | 久久草成人影院| 99精国产麻豆久久婷婷| 老司机午夜福利在线观看视频| 国产亚洲av高清不卡| 美女高潮到喷水免费观看| 最新在线观看一区二区三区| 国产精品亚洲一级av第二区| 叶爱在线成人免费视频播放| 最近最新中文字幕大全电影3 | 久久精品亚洲av国产电影网| 脱女人内裤的视频| 亚洲美女黄片视频| 成人手机av| 在线观看舔阴道视频| 午夜福利欧美成人| 热99re8久久精品国产| 下体分泌物呈黄色| 久久国产精品大桥未久av| 欧美黄色淫秽网站| 久久久精品区二区三区| 国产伦人伦偷精品视频| 人妻丰满熟妇av一区二区三区 | 久久九九热精品免费| netflix在线观看网站| 国产精品香港三级国产av潘金莲| 精品免费久久久久久久清纯 | 免费在线观看黄色视频的| 高清黄色对白视频在线免费看| 亚洲av第一区精品v没综合| 久久久国产成人免费| 国产精品一区二区在线不卡| 久久精品亚洲精品国产色婷小说| 久久精品国产a三级三级三级| 天天影视国产精品| 日本精品一区二区三区蜜桃| 69av精品久久久久久| 三级毛片av免费| 欧美大码av| 国产精品99久久99久久久不卡| 国产成人系列免费观看| 久久国产精品大桥未久av| av福利片在线| av天堂久久9| 老司机亚洲免费影院| 新久久久久国产一级毛片| 视频区欧美日本亚洲| 在线观看免费午夜福利视频| 最新美女视频免费是黄的| 又紧又爽又黄一区二区| 国产成人系列免费观看| 宅男免费午夜| 人妻一区二区av| 国产免费现黄频在线看| 欧美日韩亚洲综合一区二区三区_| 亚洲精品粉嫩美女一区| 久久精品91无色码中文字幕| 亚洲第一欧美日韩一区二区三区| 精品高清国产在线一区| 久久久久久免费高清国产稀缺| av欧美777| 人成视频在线观看免费观看| 亚洲成a人片在线一区二区| 日日摸夜夜添夜夜添小说| 妹子高潮喷水视频| 精品一品国产午夜福利视频| 狠狠婷婷综合久久久久久88av| 国产97色在线日韩免费| 欧美 亚洲 国产 日韩一| 日韩 欧美 亚洲 中文字幕| 18禁美女被吸乳视频| 亚洲精品av麻豆狂野| 另类亚洲欧美激情| 色综合欧美亚洲国产小说| 丝瓜视频免费看黄片| 国产黄色免费在线视频| 欧美不卡视频在线免费观看 | 欧美亚洲 丝袜 人妻 在线| 午夜福利影视在线免费观看| 日韩三级视频一区二区三区| 人人妻人人添人人爽欧美一区卜| 他把我摸到了高潮在线观看| 人妻久久中文字幕网| 欧美亚洲 丝袜 人妻 在线| cao死你这个sao货| 中文字幕精品免费在线观看视频| 女人久久www免费人成看片| 在线播放国产精品三级| 国产淫语在线视频| 国产精品亚洲av一区麻豆| 飞空精品影院首页| tocl精华| 久久久精品区二区三区| 国产成人影院久久av| 国产精品久久久久成人av| 亚洲熟妇熟女久久| 怎么达到女性高潮| 国产av一区二区精品久久| 免费不卡黄色视频| 丝瓜视频免费看黄片| 亚洲熟女精品中文字幕| 久久久精品区二区三区| 国产成人影院久久av| 欧美激情久久久久久爽电影 | 在线永久观看黄色视频| 欧美另类亚洲清纯唯美| 欧美在线黄色| 国产单亲对白刺激| 无遮挡黄片免费观看| 成人黄色视频免费在线看| 亚洲专区国产一区二区| 欧美日韩亚洲综合一区二区三区_| 一级毛片精品| 精品人妻熟女毛片av久久网站| 国产精品 欧美亚洲| 日本五十路高清| 香蕉国产在线看| 成年版毛片免费区| 久久中文字幕人妻熟女| 午夜日韩欧美国产| 午夜免费观看网址| 中文字幕色久视频| 精品少妇一区二区三区视频日本电影| 一进一出好大好爽视频| 精品国产亚洲在线| 成人免费观看视频高清| 69精品国产乱码久久久| 中文字幕高清在线视频| 美女午夜性视频免费| 成人黄色视频免费在线看| 欧美国产精品va在线观看不卡| 丝袜美足系列| av在线播放免费不卡| 欧美精品一区二区免费开放| 美国免费a级毛片| 欧美老熟妇乱子伦牲交| 国产色视频综合| 我的亚洲天堂| 亚洲国产精品sss在线观看 | 超碰成人久久| av在线播放免费不卡| 老鸭窝网址在线观看| 淫妇啪啪啪对白视频| 12—13女人毛片做爰片一| 国产精品一区二区在线观看99| 色综合欧美亚洲国产小说| 午夜两性在线视频| 亚洲av日韩精品久久久久久密| 可以免费在线观看a视频的电影网站| 久久香蕉精品热| 一个人免费在线观看的高清视频| 色综合欧美亚洲国产小说| 身体一侧抽搐| 天堂中文最新版在线下载| 亚洲欧美日韩高清在线视频| 巨乳人妻的诱惑在线观看| 在线观看日韩欧美| 亚洲熟女毛片儿| 叶爱在线成人免费视频播放| 国产精品一区二区在线不卡| 久久精品aⅴ一区二区三区四区| 成年女人毛片免费观看观看9 | bbb黄色大片| 女人被躁到高潮嗷嗷叫费观| 欧美人与性动交α欧美精品济南到| 不卡一级毛片| av福利片在线| av线在线观看网站| 最新在线观看一区二区三区| 日本欧美视频一区| 欧美中文综合在线视频| 亚洲精品久久成人aⅴ小说| 欧美一级毛片孕妇| 在线观看免费视频网站a站| 精品欧美一区二区三区在线| avwww免费| 婷婷丁香在线五月| 天天影视国产精品| 午夜老司机福利片| 成人特级黄色片久久久久久久| 久久精品国产亚洲av高清一级| 日日爽夜夜爽网站| 91麻豆精品激情在线观看国产 | 亚洲中文日韩欧美视频| 黄片小视频在线播放| 国产单亲对白刺激| 国产男靠女视频免费网站| 在线观看免费午夜福利视频| 麻豆乱淫一区二区| 18禁美女被吸乳视频| 中文字幕高清在线视频| 91精品三级在线观看| 伦理电影免费视频| 国产在线观看jvid| 男女下面插进去视频免费观看| 亚洲中文日韩欧美视频| 人人妻人人添人人爽欧美一区卜| 日韩欧美一区二区三区在线观看 | 欧美乱色亚洲激情| 亚洲精品粉嫩美女一区| 高清av免费在线| 成人国产一区最新在线观看| 91国产中文字幕| 老熟妇乱子伦视频在线观看| 亚洲性夜色夜夜综合| 欧美最黄视频在线播放免费 | 新久久久久国产一级毛片| 麻豆国产av国片精品| 成年版毛片免费区| 久久中文字幕人妻熟女| 真人做人爱边吃奶动态| 美女高潮喷水抽搐中文字幕| 可以免费在线观看a视频的电影网站| 国产亚洲一区二区精品| 露出奶头的视频| 伦理电影免费视频| 日韩成人在线观看一区二区三区| 久久性视频一级片| 又黄又爽又免费观看的视频| 亚洲熟妇熟女久久| 成人影院久久| 亚洲欧美色中文字幕在线| 在线观看日韩欧美| 麻豆乱淫一区二区| 精品免费久久久久久久清纯 | 大香蕉久久网| 搡老乐熟女国产| 国产真人三级小视频在线观看| 99久久国产精品久久久| 在线观看免费视频日本深夜| 黄频高清免费视频| 9191精品国产免费久久| 久久久久精品人妻al黑| 男女高潮啪啪啪动态图| 村上凉子中文字幕在线| 美女扒开内裤让男人捅视频| 18禁国产床啪视频网站| 日日摸夜夜添夜夜添小说| 午夜福利在线免费观看网站| 亚洲自偷自拍图片 自拍| 天天操日日干夜夜撸| 黄色成人免费大全| 一区在线观看完整版| 韩国精品一区二区三区| 婷婷丁香在线五月| a级片在线免费高清观看视频| 两个人看的免费小视频| 色94色欧美一区二区| 国产成+人综合+亚洲专区| 99国产精品一区二区蜜桃av | 午夜老司机福利片| 日本欧美视频一区| 久久影院123| 亚洲精品粉嫩美女一区| 高清欧美精品videossex| 欧美精品人与动牲交sv欧美| 黄网站色视频无遮挡免费观看| 欧美日韩一级在线毛片| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久精品国产欧美久久久| 天堂动漫精品| 欧美乱码精品一区二区三区| 国产成人精品无人区| 亚洲三区欧美一区| 人人妻人人澡人人看| 精品亚洲成a人片在线观看| 久久久久国产一级毛片高清牌| 国产欧美日韩一区二区三区在线| 日韩有码中文字幕| tube8黄色片| 国产精品永久免费网站| 在线天堂中文资源库| 亚洲欧美一区二区三区久久| 老鸭窝网址在线观看| 妹子高潮喷水视频| 99热网站在线观看| 久久精品aⅴ一区二区三区四区| 久久久久久久久免费视频了| 日韩视频一区二区在线观看| 97人妻天天添夜夜摸| 欧美精品一区二区免费开放| 日本五十路高清| 国产精品免费大片| 99re6热这里在线精品视频| 老汉色av国产亚洲站长工具| 啦啦啦在线免费观看视频4| 日韩熟女老妇一区二区性免费视频| 国产成人精品久久二区二区免费| 男女下面插进去视频免费观看| 亚洲三区欧美一区| 老司机午夜十八禁免费视频| 精品一区二区三卡| 人人妻人人添人人爽欧美一区卜| 欧美日韩亚洲综合一区二区三区_| 水蜜桃什么品种好| 别揉我奶头~嗯~啊~动态视频| 亚洲成人免费av在线播放| 丰满人妻熟妇乱又伦精品不卡| 欧美日本中文国产一区发布| 99国产精品一区二区三区| 日韩大码丰满熟妇| 人人妻,人人澡人人爽秒播| 熟女少妇亚洲综合色aaa.| 天堂动漫精品| 国产免费男女视频| 国产精华一区二区三区| 极品少妇高潮喷水抽搐| 法律面前人人平等表现在哪些方面| 99热只有精品国产| 午夜91福利影院| 亚洲精华国产精华精| 一进一出抽搐gif免费好疼 | 欧美性长视频在线观看| 18在线观看网站| 黑人巨大精品欧美一区二区蜜桃| 久久久国产一区二区| 99国产精品免费福利视频| 国产在视频线精品| 亚洲五月婷婷丁香| 国产在线一区二区三区精| 日韩欧美在线二视频 | 在线十欧美十亚洲十日本专区| 男人的好看免费观看在线视频 | 天天添夜夜摸| 又紧又爽又黄一区二区| 很黄的视频免费| 久久精品国产亚洲av香蕉五月 | 99国产精品一区二区蜜桃av | 老司机午夜福利在线观看视频| 国产精品免费一区二区三区在线 | 黄色视频不卡| 女性被躁到高潮视频| tocl精华| 国产又爽黄色视频| 亚洲精品久久午夜乱码| 男女免费视频国产| 少妇粗大呻吟视频| 国产亚洲欧美精品永久| 一二三四社区在线视频社区8| 亚洲自偷自拍图片 自拍| 久久性视频一级片| 亚洲国产精品sss在线观看 | 国产精品久久久久成人av| 国产精品成人在线| 国产在线观看jvid| 一级毛片女人18水好多| 最近最新中文字幕大全免费视频| 国产99白浆流出| 成年女人毛片免费观看观看9 | 国内久久婷婷六月综合欲色啪| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲第一青青草原| 女人爽到高潮嗷嗷叫在线视频| 亚洲精华国产精华精| www.自偷自拍.com| 久久久久精品国产欧美久久久| 亚洲人成77777在线视频| 国产精品99久久99久久久不卡| 黄色片一级片一级黄色片| 两个人免费观看高清视频| 热re99久久国产66热| 1024视频免费在线观看| 黑人巨大精品欧美一区二区蜜桃| 国产人伦9x9x在线观看| 精品国产乱码久久久久久男人| 久久久久久免费高清国产稀缺| 成年女人毛片免费观看观看9 | 正在播放国产对白刺激| 久久中文字幕人妻熟女| 妹子高潮喷水视频| 老司机福利观看| 国产三级黄色录像| 成年版毛片免费区| 精品一区二区三区四区五区乱码| 国产不卡av网站在线观看| 国产亚洲精品一区二区www | tube8黄色片| 精品国产乱子伦一区二区三区| 亚洲一码二码三码区别大吗| 男人舔女人的私密视频| www.熟女人妻精品国产| 俄罗斯特黄特色一大片| 不卡一级毛片| 校园春色视频在线观看| 成人亚洲精品一区在线观看| 午夜两性在线视频| 老司机午夜福利在线观看视频| 国产精品久久久人人做人人爽| 夜夜躁狠狠躁天天躁| 日韩视频一区二区在线观看| 美女国产高潮福利片在线看| 成人手机av| 人人妻,人人澡人人爽秒播| 精品午夜福利视频在线观看一区| 久久天堂一区二区三区四区| 看免费av毛片| 俄罗斯特黄特色一大片| 黄色成人免费大全| 日韩有码中文字幕| 日本欧美视频一区| 亚洲情色 制服丝袜| 国产精品一区二区在线观看99| 欧美日韩av久久| 国产精品成人在线| av天堂在线播放| 欧美日韩亚洲国产一区二区在线观看 | 天天影视国产精品| 伦理电影免费视频| 中文字幕人妻丝袜一区二区| 啦啦啦视频在线资源免费观看| 黄色视频,在线免费观看| 亚洲一区中文字幕在线| 国产精华一区二区三区| 男人舔女人的私密视频| 久久草成人影院| 久久天堂一区二区三区四区| 日日夜夜操网爽| 狠狠婷婷综合久久久久久88av| av免费在线观看网站| 免费观看a级毛片全部| 亚洲精品在线观看二区| 熟女少妇亚洲综合色aaa.| 高清视频免费观看一区二区| 亚洲成人手机| 免费人成视频x8x8入口观看| 国产精品久久久久久精品古装| 色婷婷av一区二区三区视频| 一区二区日韩欧美中文字幕| 久久青草综合色| 日韩欧美免费精品| 在线观看免费视频日本深夜|