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

    100%低地板齒輪箱振動(dòng)分析及噪聲預(yù)估

    2021-01-09 02:37:46彭銳濤謝炎培羅秀強(qiáng)胡聰芳劉祥環(huán)
    關(guān)鍵詞:模態(tài)有限元振動(dòng)

    彭銳濤,謝炎培,羅秀強(qiáng),胡聰芳,劉祥環(huán)

    100%低地板齒輪箱振動(dòng)分析及噪聲預(yù)估

    彭銳濤,謝炎培,羅秀強(qiáng),胡聰芳,劉祥環(huán)

    (湘潭大學(xué) 機(jī)械工程學(xué)院,湖南 湘潭 411105)

    對(duì)100%低地板齒輪箱進(jìn)行振動(dòng)分析和噪聲預(yù)估,可為其設(shè)計(jì)及優(yōu)化提供依據(jù)。建立了傳動(dòng)系統(tǒng)和結(jié)構(gòu)系統(tǒng)相耦合的100%低地板齒輪箱剛?cè)狁詈隙囿w動(dòng)力學(xué)模型,獲得了軸承支反力;建立了齒輪箱模態(tài)分析模型,研究了箱體的振動(dòng)特性?;谙嗨菩岳碚撛O(shè)計(jì)了箱體的相似模型并對(duì)其進(jìn)行了實(shí)驗(yàn)?zāi)B(tài)分析,驗(yàn)證了有限元模型的正確。將剛?cè)狁詈隙囿w動(dòng)力學(xué)仿真獲得的軸承支反力作為箱體的激勵(lì),采用直接積分法計(jì)算箱體的振動(dòng)響應(yīng);以箱體振動(dòng)速度作為邊界條件,建立了齒輪箱直接邊界元法輻射噪聲分析模型,進(jìn)行齒輪箱表面聲壓及外聲場(chǎng)輻射噪聲預(yù)估。結(jié)果表明:齒輪箱輻射噪聲頻域曲線(xiàn)的噪聲峰值頻率出現(xiàn)在齒輪副嚙合頻率及其諧波,以及箱體低階固有頻率處。

    100%低地板;齒輪箱;剛?cè)狁詈蟿?dòng)力學(xué);振動(dòng)特性;相似性理論;輻射噪聲

    低地板軌道交通系統(tǒng)屬于中等運(yùn)輸量的交通方式,填補(bǔ)了當(dāng)前公共交通方式運(yùn)量的空白,同時(shí)具有節(jié)能環(huán)保、舒適安全、建造價(jià)格低、便于乘客上下車(chē)和兼具城市觀光功能等優(yōu)點(diǎn),適合用作大城市的支線(xiàn)交通[1?2]。因此,低地板軌道交通系統(tǒng)這種新型的城市交通運(yùn)輸形式,正受到越來(lái)越多的關(guān)注和歡迎。100%低地板軌道交通系統(tǒng)由于整車(chē)地板平面距離軌道平面高度在400 mm以下,所以可用于安裝車(chē)輛設(shè)備的空間狹小[3]。100%低地板軌道交通系統(tǒng)傳遞動(dòng)力采用的方式是使用齒輪箱直接驅(qū)動(dòng)車(chē)輪,其齒輪箱的結(jié)構(gòu)具有占用空間小、外形緊湊的特點(diǎn)。100%低地板軌道交通系統(tǒng)的低地板結(jié)構(gòu)特點(diǎn)使齒輪箱這個(gè)重要的噪聲源離車(chē)廂內(nèi)乘客的距離變近,并且輕量化的車(chē)體設(shè)計(jì)以及玻璃材料的大量使用使得車(chē)輛的隔聲能力變得更加薄弱,這對(duì)100%低地板軌道交通齒輪箱的振動(dòng)及噪聲控制性能提出了更高的要求[4]。車(chē)輛運(yùn)行過(guò)程中齒輪箱的振動(dòng)及噪聲不僅會(huì)影響齒輪箱的使用壽命,而且也會(huì)影響乘客的乘坐體驗(yàn)[5?6]。因此對(duì)齒輪箱開(kāi)展振動(dòng)及噪聲性能的研究是100%低地板軌道交通系統(tǒng)設(shè)計(jì)的一項(xiàng)重要研究?jī)?nèi)容。眾多的學(xué)者對(duì)齒輪箱的振動(dòng)及噪聲性能進(jìn)行了深入研究。ZHOU等[7?8]分析了齒輪傳動(dòng)系統(tǒng)的內(nèi)部激勵(lì),獲得了軸承力,并將其作為箱體的激勵(lì),進(jìn)行瞬態(tài)動(dòng)力學(xué)分析獲得了箱體的振動(dòng)響應(yīng)。林騰蛟等[9]建立了船用齒輪箱的多剛體動(dòng)力學(xué)模型,運(yùn)用邊界元法預(yù)估了齒輪箱的輻射噪聲。陸波等[10]將某大型船用齒輪箱的傳動(dòng)系統(tǒng)和結(jié)構(gòu)系統(tǒng)耦合了起來(lái),綜合考慮系統(tǒng)內(nèi)部激勵(lì)和外部激勵(lì)進(jìn)行齒輪箱的結(jié)構(gòu)噪聲預(yù)估。任亞峰等[11]研究了箱體的柔性對(duì)齒輪傳動(dòng)系統(tǒng)動(dòng)態(tài)特性的影響。Abbes等[12]考慮齒輪副的時(shí)變嚙合剛度,使用聲固耦合方法開(kāi)展了齒輪箱輻射噪聲的分析。但是考慮箱體的柔性以及傳動(dòng)系統(tǒng)和結(jié)構(gòu)系統(tǒng)的耦合作用進(jìn)行齒輪箱振動(dòng)及輻射噪聲研究并進(jìn)行試驗(yàn)驗(yàn)證的文獻(xiàn)較少。本文建立了結(jié)構(gòu)系統(tǒng)與傳動(dòng)系統(tǒng)相耦合的100%低地板齒輪箱剛?cè)狁詈隙囿w動(dòng)力學(xué)模型,獲得了軸承支反力;而后對(duì)箱體進(jìn)行有限元模態(tài)分析,分析箱體的振動(dòng)特性?;谙嗨菩岳碚撛O(shè)計(jì)了箱體相似模型,采用3D打印技術(shù)制作了相似模型,并對(duì)其進(jìn)行實(shí)驗(yàn)?zāi)B(tài)分析驗(yàn)證有限元模型的正確性。將軸承支反力加載到箱體軸承座位置,使用直接積分法計(jì)算箱體的振動(dòng)響應(yīng);建立100%低地板齒輪箱直接邊界元法輻射噪聲分析模型,將獲得的箱體振動(dòng)響應(yīng)作為邊界條件進(jìn)行齒輪箱表面聲壓及輻射噪聲預(yù)估。

    1 剛?cè)狁詈隙囿w動(dòng)力學(xué)

    1.1 齒輪箱三維建模

    100%低地板齒輪箱的傳動(dòng)系統(tǒng)由三級(jí)齒輪傳動(dòng)組成,輸入級(jí)是錐齒輪傳動(dòng)、中間級(jí)和輸出級(jí)是斜齒輪傳動(dòng)。在表1中列出了齒輪箱各傳動(dòng)級(jí)齒輪齒數(shù)。在SolidWorks軟件中建立了齒輪箱的三維模型,如圖1所示。傳動(dòng)系統(tǒng)中齒輪和軸的材料為18CrNiMo7,其彈性模量為1.9×1011Pa,密度為7 800 kg/m3,泊松比為0.30;結(jié)構(gòu)系統(tǒng)的材料為QT400-18,其密度為7 000 kg/m3,彈性模量為2.0×105MPa,泊松比為0.32。

    表1 齒輪箱齒數(shù)

    圖1 100%低地板齒輪箱模型

    1.2 約束和載荷

    根據(jù)齒輪箱各部件間的運(yùn)動(dòng)關(guān)系,在Adams軟件中建立結(jié)構(gòu)系統(tǒng)和傳動(dòng)系統(tǒng)相耦合的剛?cè)狁詈隙囿w動(dòng)力學(xué)模型,使用軸承單元將剛性傳動(dòng)系統(tǒng)與柔性結(jié)構(gòu)系統(tǒng)連接。圖2為所建立的齒輪箱剛?cè)狁詈隙囿w動(dòng)力學(xué)模型。

    圖2 剛?cè)狁詈夏P?/p>

    根據(jù)100%低地板齒輪箱的實(shí)際運(yùn)行工況,定義輸入軸的輸入轉(zhuǎn)速為2 100 r/min,負(fù)載力矩為3 068.4 N?m。ADAMS軟件采用沖擊函數(shù)法計(jì)算法向接觸力,采用庫(kù)侖摩擦方法計(jì)算切向接觸力。輸入級(jí)齒輪副的剛度系數(shù)為4.805×105N/mm,中間級(jí)齒輪副的剛度系數(shù)為7.933×105N/mm,輸出級(jí)齒輪副的剛度系數(shù)為1.096×106N/mm。阻尼系數(shù)通常取剛度系數(shù)的0.1%~1%,這里取0.1%。根據(jù)齒輪對(duì)的材料屬性選取非線(xiàn)性彈簧力指數(shù)為1.5,最大穿透深度為0.1 mm。

    圖3 X方向軸承支反力

    1.3 剛?cè)狁詈蟿?dòng)力學(xué)仿真

    使用STEP函數(shù)定義轉(zhuǎn)速以及負(fù)載扭矩在0~0.1 s時(shí)間段從0開(kāi)始遞增,設(shè)置仿真時(shí)間為0.5 s,進(jìn)行齒輪箱剛?cè)狁詈蟿?dòng)力學(xué)仿真。輸出軸最外側(cè)軸承方向的軸承支反力如圖3所示。仿真獲得了9個(gè)軸承的軸承支反力,并對(duì)仿真結(jié)果提取0.1~0.5 s穩(wěn)定階段軸承支反力,將其作為箱體動(dòng)力學(xué)仿真的激勵(lì)。

    2 齒輪箱箱體振動(dòng)特性分析

    2.1 有限元模型

    從低地板齒輪箱的三維模型可知箱體結(jié)構(gòu)比較復(fù)雜,考慮到計(jì)算效率和網(wǎng)格劃分難度等問(wèn)題,需要對(duì)導(dǎo)入ABAQUS軟件的箱體模型進(jìn)行一定程度的簡(jiǎn)化,以此來(lái)降低網(wǎng)格劃分的難度、改善網(wǎng)格質(zhì)量和提高有限元軟件的運(yùn)算效率。對(duì)箱體的小圓角、小孔、小凸臺(tái)等對(duì)齒輪箱箱體整體質(zhì)量和剛度影響較小的一些特征進(jìn)行簡(jiǎn)化。同時(shí)在建立齒輪箱箱體模型時(shí),不考慮左右箱體結(jié)合部、軸承端蓋與箱體結(jié)合部的影響[13]。使用自由網(wǎng)格劃分方式直接進(jìn)行箱體的網(wǎng)格劃分,共計(jì)單元49 364個(gè),節(jié)點(diǎn)87 629個(gè)。

    圖4 箱體有限元模型

    2.2 模態(tài)分析

    根據(jù)100%低地板齒輪箱的實(shí)際約束情況,對(duì)箱體進(jìn)行有限元約束模態(tài)分析。齒輪箱箱體通過(guò)螺栓與轉(zhuǎn)向架連接,因此約束箱體螺栓孔內(nèi)表面6個(gè)自由度,箱體有限元模型如圖4所示。在有限元軟件中使用Lanczos法進(jìn)行齒輪箱箱體的約束模態(tài)分析,獲得箱體的前10階固有頻率及模態(tài)振型。表2列出了齒輪箱箱體的前6階固有頻率,箱體前6階模態(tài)振型如圖5所示。

    箱體前6階模態(tài)振型中,除了第3階是扭轉(zhuǎn)模態(tài)振型,其余的都是擺振模態(tài)振型。箱體容易出現(xiàn)振動(dòng)的位置在輸出軸軸承座右側(cè)箱體以及上下側(cè)箱體。箱體前10階固有頻率沒(méi)有和齒輪箱的嚙合頻率及其諧波相接近或重合,因此齒輪箱嚙合頻率不會(huì)引起箱體的共振。

    表2 箱體固有頻率

    (a) 1階振型;(b) 2階振型;(c) 3階振型;(d) 4階振型;(e) 5階振型;(f) 6階振型

    2.3 動(dòng)態(tài)響應(yīng)仿真

    將齒輪箱剛?cè)狁詈隙囿w動(dòng)力學(xué)分析獲得的軸承支反力作為箱體激勵(lì),施加在齒輪箱箱體各對(duì)應(yīng)軸承孔耦合點(diǎn)上,如圖6所示。在ABAQUS軟件中采取直接積分法進(jìn)行箱體的瞬態(tài)動(dòng)力學(xué)仿真,獲得了箱體的振動(dòng)響應(yīng)。將瞬態(tài)動(dòng)力學(xué)仿真獲得的箱體振動(dòng)速度作為齒輪箱直接邊界元法輻射噪聲分析模型的邊界條件。

    圖6 動(dòng)態(tài)響應(yīng)模型

    3 箱體相似模型設(shè)計(jì)

    3.1 求取π因子

    相似模型設(shè)計(jì)的關(guān)鍵是基于相似理論推導(dǎo)出相似準(zhǔn)則[14],本文采用量綱分析法來(lái)導(dǎo)出相似準(zhǔn)則。對(duì)于齒輪箱的振動(dòng)問(wèn)題,其振動(dòng)特性應(yīng)當(dāng)與下列物理參數(shù)有關(guān):結(jié)構(gòu)的特征長(zhǎng)度,密度,彈性模量,泊松比,振動(dòng)圓頻率,振動(dòng)位移,速度,加速度,作用力,以及頻響函數(shù)。取基本物理量單位為質(zhì)量,長(zhǎng)度,以及時(shí)間,該系統(tǒng)的無(wú)量綱方程如式(1)所示,可用式(2)表示。

    上述物理量與基本量的關(guān)系可用表3的量綱矩陣表示。

    根據(jù)因子無(wú)量綱的要求,可得到如下的代數(shù)方程式:

    將長(zhǎng)度,密度以及彈性模量作為基本未知量,則上式表示為:

    對(duì)應(yīng)的矩陣如表4所示。

    表3 量綱矩陣

    表4 π矩陣

    通過(guò)計(jì)算可得7個(gè)因子:

    定義為原型與相似模型各物理量的相似比,根據(jù)上式可得:

    3.2 相似模型制作

    要做到試驗(yàn)?zāi)P团c原模型完全嚴(yán)格相似是很困難的,只能根據(jù)所研究的問(wèn)題,捉住主要影響因素放棄次要影響因素,使模型試驗(yàn)既能盡可能反應(yīng)所研究問(wèn)題的真實(shí)情況,而又不至太復(fù)雜[15?16]。齒輪箱箱體結(jié)構(gòu)復(fù)雜且較大,采用傳統(tǒng)的方法難以快速且經(jīng)濟(jì)的制造出來(lái)。隨著3D打印技術(shù)的發(fā)展,復(fù)雜結(jié)構(gòu)的快速制造成為了可能。出于經(jīng)濟(jì)性的考慮,采用白色樹(shù)脂材料來(lái)進(jìn)行齒輪箱箱體的制造,其密度為1 300 kg/m3,彈性模量為2 650 MPa,泊松比為0.41。

    由于采用了不同的材料制作模型,所以泊松比比尺λ不等于1,無(wú)法滿(mǎn)足相似條件。然而現(xiàn)有研究表明在對(duì)相似模型進(jìn)行振動(dòng)特性實(shí)驗(yàn)測(cè)試時(shí),可以不考慮泊松比的影響,忽略該參數(shù)的影響仍然可以取得較高的精度[17?18]。由于相似模型需要進(jìn)行實(shí)驗(yàn)?zāi)B(tài)分析,在實(shí)驗(yàn)過(guò)程中需要在箱體表面安裝傳感器,如果模型寸尺太小,則會(huì)使得傳感器附加質(zhì)量對(duì)實(shí)驗(yàn)測(cè)試數(shù)據(jù)造成較大的誤差;如果模型尺寸太大則會(huì)造成箱體3D打印制作成本的增加。所以綜合考慮各因素,對(duì)原型進(jìn)行1/2等比例的縮小,即幾何比尺=2。模型采用與原型相同的邊界條件,使用螺栓將相似模型固定在測(cè)試臺(tái)架上。

    3.3 模態(tài)參數(shù)換算

    3.3.1 模態(tài)振型

    設(shè)相似模型和原型的某階振型分別為{}和{},根據(jù)7因子,對(duì)于規(guī)定的測(cè)量點(diǎn)有:

    式中:1為幾何比尺,為彈性模量比尺,′為一常數(shù),因此{(lán)}和{}只差一個(gè)常數(shù)比值。由于模態(tài)振型是各點(diǎn)的振幅比,所以相似模型與原型的同階模態(tài)振型相同。

    3.3.2 模態(tài)頻率

    根據(jù)前文計(jì)算出來(lái)的2因子可得:

    式中:為幾何比尺,為彈性模量比尺,為密度比尺。將材料參數(shù)代入上式求解可得模態(tài)頻率比尺為λ=1.87。

    4 模態(tài)實(shí)驗(yàn)

    4.1 相似模型實(shí)驗(yàn)?zāi)B(tài)分析

    使用丹麥B & K公司生產(chǎn)的B & K PULSE噪聲測(cè)試系統(tǒng)進(jìn)行模態(tài)實(shí)驗(yàn),模態(tài)測(cè)試系統(tǒng)由力錘、測(cè)試軟件、加速度傳感器以及數(shù)據(jù)采集卡4個(gè)部分組成,測(cè)試系統(tǒng)如圖7所示。

    圖7 模態(tài)測(cè)試系統(tǒng)

    相似模型的結(jié)構(gòu)比較小,布置多個(gè)傳感器會(huì)因?yàn)檫^(guò)多的傳感器附加質(zhì)量而對(duì)測(cè)試結(jié)果產(chǎn)生較大的誤差。根據(jù)測(cè)試經(jīng)驗(yàn),使用2個(gè)加速度傳感器進(jìn)行實(shí)驗(yàn)測(cè)試。相似模型是樹(shù)脂材料,不具有磁性,采用磁性貼片將傳感器和箱體模型相連接。由于實(shí)驗(yàn)難以獲得相似模型扭振模態(tài)參數(shù),因此只對(duì)擺振模態(tài)參數(shù)進(jìn)行測(cè)量。實(shí)驗(yàn)測(cè)試如圖8所示。經(jīng)過(guò)多次測(cè)試,得到的相似模型前六階模態(tài)頻率如表5 所示。

    圖8 模態(tài)測(cè)試現(xiàn)場(chǎng)

    4.2 仿真數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)對(duì)比

    將箱體有限元模態(tài)分析獲得的模態(tài)參數(shù)與相似模型實(shí)驗(yàn)?zāi)B(tài)分析獲得的模態(tài)參數(shù)進(jìn)行對(duì)比,發(fā)現(xiàn)實(shí)驗(yàn)?zāi)B(tài)分析獲得的模態(tài)參數(shù)沒(méi)有能與有限元模態(tài)分析獲得的第3階模態(tài)參數(shù)相對(duì)應(yīng)。從有限元分析結(jié)果知,箱體第3階模態(tài)振型為扭振模態(tài);又由實(shí)驗(yàn)只獲得了相似模型的擺振模態(tài)參數(shù),因此實(shí)驗(yàn)獲得的結(jié)果無(wú)法與有限元分析獲得的第3階模態(tài)參數(shù)相對(duì)應(yīng)。由上文可知箱體與相似模型的模態(tài)頻率比尺λ為1.87。將仿真得到的頻率與相似模型實(shí)驗(yàn)得到的頻率進(jìn)行頻率比計(jì)算,并計(jì)算其與模態(tài)頻率比尺的相對(duì)誤差,相對(duì)誤差如表6所示。相對(duì)誤差全部低于5%,這說(shuō)明箱體與相似模型具有動(dòng)力學(xué)相似,同時(shí)也說(shuō)明了箱體有限元仿真模型的正確。

    表5 實(shí)驗(yàn)?zāi)B(tài)前6階固有頻率

    表6 頻率比與相對(duì)誤差

    5 齒輪箱輻射噪聲預(yù)估

    5.1 聲輻射分析模型建模

    在LMS Virtual.Lab軟件中使用直接邊界元方法進(jìn)行輻射噪聲的求解,要求聲學(xué)網(wǎng)格是封閉的。所以需要將齒輪箱箱體表面封閉,采用殼單元來(lái)對(duì)箱體進(jìn)行封閉;最后將箱體表面提取的聲學(xué)網(wǎng)格和殼單元網(wǎng)格進(jìn)行合并,形成一個(gè)完整的箱體聲學(xué)網(wǎng)格。齒輪箱輻射噪聲分析模型的邊界條件為瞬態(tài)動(dòng)力學(xué)仿真得到的箱體表面振動(dòng)速度。

    聲學(xué)邊界元網(wǎng)格劃分的過(guò)于粗糙會(huì)導(dǎo)致計(jì)算數(shù)據(jù)產(chǎn)生較大誤差。對(duì)于聲學(xué)邊界元模型來(lái)說(shuō),通常假設(shè)在最小波長(zhǎng)內(nèi)有6個(gè)單元[19]??紤]軟件的計(jì)算效率及計(jì)算精度,確定聲學(xué)網(wǎng)格的單元長(zhǎng)度為16 mm。聲學(xué)網(wǎng)格共計(jì)單元6 655個(gè),節(jié)點(diǎn)6 657個(gè)。如圖9為齒輪箱的聲學(xué)網(wǎng)格。進(jìn)行齒輪箱的輻射噪聲預(yù)估需要在齒輪箱外部設(shè)置一個(gè)矩形場(chǎng)點(diǎn)。根據(jù)齒輪裝置的驗(yàn)收規(guī)范?空氣傳播噪聲的實(shí)驗(yàn)規(guī) 范[20],設(shè)置矩形場(chǎng)點(diǎn)各平面距離齒輪箱大約1 m。同時(shí)在矩形場(chǎng)點(diǎn)4個(gè)面上各選取了1個(gè)聲場(chǎng)點(diǎn)。齒輪箱外部聲場(chǎng)網(wǎng)格及4個(gè)聲場(chǎng)點(diǎn)分布如圖10所示。

    圖9 聲學(xué)網(wǎng)格

    圖10 外聲場(chǎng)網(wǎng)格

    5.2 表面聲壓及場(chǎng)點(diǎn)聲壓計(jì)算

    設(shè)定空氣密度為1.225 kg/m3,聲速為340 m/s,參考聲壓為2×10?5Pa。采用聲學(xué)直接邊界元法求解,可得到齒輪箱表面各個(gè)頻率下的聲壓大小和位置分布。圖11展示了齒輪箱在頻率為310 Hz時(shí)的表面聲壓的有效值及其分布,其聲壓最大值為120 dB。齒輪箱外聲場(chǎng)在310 Hz時(shí)的輻射噪聲如圖12所示,其聲壓峰值為101 dB。

    圖11 表面聲壓

    圖12 外聲場(chǎng)聲壓

    圖13 4個(gè)聲場(chǎng)點(diǎn)噪聲值

    通過(guò)仿真獲得了齒輪箱矩形場(chǎng)點(diǎn)中4個(gè)聲場(chǎng)點(diǎn)的頻率噪聲曲線(xiàn)。齒輪箱外部聲場(chǎng)的聲壓峰值主要集中在0~2 000 Hz這個(gè)頻率范圍;并且4個(gè)聲場(chǎng)計(jì)算點(diǎn)的噪聲值隨著頻率的變化規(guī)律基本一致。圖13給出了4個(gè)聲場(chǎng)點(diǎn)0~2 000 Hz的頻率噪聲曲線(xiàn)。齒輪箱輻射噪聲峰值頻率主要集中在嚙合頻率及其諧波處,如a,b,c,d,e,以及齒輪箱低階固有頻率處,如1,2,3,4,5,6。在齒輪箱前6階固有頻率中,只有第3階固有頻率處沒(méi)出現(xiàn)明顯的輻射噪聲峰值,這表明扭振模態(tài)對(duì)齒輪箱輻射噪聲的影響較小。從出現(xiàn)輻射噪聲峰值的頻率可知,齒輪箱的輻射噪聲受?chē)Ш项l率以及箱體固有頻率的影響較大,因此對(duì)其的優(yōu)化設(shè)計(jì)可以從降低齒輪嚙合振動(dòng)以及箱體的結(jié)構(gòu)優(yōu)化兩方面進(jìn)行。

    6 結(jié)論

    1) 建立了100%低地板齒輪箱有限元分析模型,對(duì)齒輪箱箱體的振動(dòng)特性進(jìn)行了分析,表明齒輪副的嚙合頻率不會(huì)激起箱體的共振。

    2) 基于相似性理論設(shè)計(jì)了箱體的相似模型,結(jié)合3D打印技術(shù)制造了相似模型,并進(jìn)行了相似模型的實(shí)驗(yàn)?zāi)B(tài)分析,驗(yàn)證了箱體有限元模型的正確性。該方法為大型復(fù)雜齒輪箱振動(dòng)噪聲分析的實(shí)驗(yàn)驗(yàn)證提供了一種思路,可有效降低實(shí)驗(yàn)費(fèi)用。

    3) 建立了100%低地板齒輪箱直接邊界元輻射噪聲分析模型,預(yù)估了齒輪箱表面聲壓以及齒輪箱外部聲場(chǎng)輻射噪聲。結(jié)果顯示齒輪箱輻射噪聲頻率范圍主要集中在0~2 000 Hz,其峰值主要出現(xiàn)在齒輪箱嚙合頻率及其諧波以及箱體低階固有頻率處。表明對(duì)100%低地板齒輪箱的優(yōu)化設(shè)計(jì)需要從降低齒輪嚙合振動(dòng)以及箱體結(jié)構(gòu)優(yōu)化兩方面同時(shí)進(jìn)行,同時(shí)需要著重關(guān)注箱體的低階擺振模態(tài)振型。

    [1] Currie G, Burke M. Light rail in australia: Performance and prospects[C]// Australasian Transport Research Forum, Brisbane, 2013.

    [2] Vukan R. Vuchic. Comparison of light rail transit with bus semirapid transit[C]// 5th UITP Light Rail Conference, Melbourne, 2000.

    [3] 牛悅丞, 李芾, 楊陽(yáng), 等. 國(guó)外新一代100%低地板有軌電車(chē)技術(shù)特征及發(fā)展趨勢(shì)[J]. 鐵道標(biāo)準(zhǔn)設(shè)計(jì), 2018, 62(11): 148?153. NIU Yuecheng, LI Fu, YANG Yang, et al. Technical characteristics and development tendency of new generation foreign 100% low floor tram[J]. Railway Standard Design, 2018, 62(11): 148?153.

    [4] Tosun M, Y?ld?z M, ?zkan A. Investigation and refinement of gearbox whine noise[J]. Applied Acoustics, 2018(130): 305?311.

    [5] Eade P W, Hardy A E J. Railway vehicle internal noise[J]. Journal of Sound & Vibration, 1977, 51(3): 403?415.

    [6] Aherwar A. An investigation on gearbox fault detection using vibration analysis techniques: A review[J]. Australian Journal of Mechanical Engineering, 2011, 10(2): 169?183.

    [7] ZHOU Jianxing, S Wenlei, T Qing. Gearbox low-noise design method based on panel acoustic contribution[J]. Mathematical Problems in Engineering, 2014, 2014: 1? 10.

    [8] 黃冠華, 張衛(wèi)華, 梅桂明, 等. 內(nèi)外激勵(lì)下高速列車(chē)齒輪箱箱體動(dòng)態(tài)響應(yīng)分析[J]. 機(jī)械工程學(xué)報(bào), 2015, 51(12): 95?100. HUANG Guanhua, ZHANG Weihua, MEI Guiming, et al. Dynamic response analysis of gearbox housing system subjected to internal and external excitation in high-speed train[J]. Journal of Mechanical Engineering, 2015, 51(12): 95?100.

    [9] 林騰蛟, 何澤銀, 鐘聲, 等. 船用齒輪箱多體動(dòng)力學(xué)仿真及聲振耦合分析[J]. 湖南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2015, 42(2): 22?28. LIN Tengjiao, HE Zeyin, ZHONG Sheng, et al. Multi- body dynamic simulation and vibro-acoustic coupling analysis of marine gearbox[J]. Journal of Hunan University (Natural Sciences), 2015, 42(2): 22?28.

    [10] 陸波, 朱才朝, 宋朝省, 等. 大功率船用齒輪箱耦合非線(xiàn)性動(dòng)態(tài)特性分析及噪聲預(yù)估[J]. 振動(dòng)與沖擊, 2009, 28(4): 76?80. LU Bo, ZHU Caichao, SONG Chaosheng, et al. Coupled nonlinear dynamic characteristics analysis and noise predestination of a large burden marine gearbox[J]. Journal of Vibration and Shock, 2009, 28(4): 76?80.

    [11] 任亞峰, 常山, 劉更, 等. 箱體柔性對(duì)齒輪傳動(dòng)系統(tǒng)動(dòng)態(tài)特性的影響分析[J]. 振動(dòng)與沖擊, 2017, 36(14): 85? 91. REN Yafeng, CHANG Shan, LIU Geng, et al. Influence of housing flexibility on the dynamic characteristics of gear transmission systems[J]. Journal of Vibration and Shock, 2017, 36(14): 85?91.

    [12] Abbes M S, Bouaziz S, Chaari F, et al. An acoustic– structural interaction modelling for the evaluation of a gearbox-radiated noise[J]. International Journal of Mechanical Sciences, 2008, 50(3): 569?577.

    [13] 王富民, 李捷, 楊建偉, 等. 地鐵齒輪箱箱體模態(tài)及諧響應(yīng)分析[J]. 機(jī)械傳動(dòng), 2015(9): 146?150. WANG Fuming, LI Jie, YANG Jianwei, et al. Modal and harmonic response analysis of subway gearbox housing[J]. Journal of Mechanical Transmission, 2015(9): 146?150.

    [14] C Yuxi, G Shen, H Wang. Maglev vehicle-guideway coupling vibration test rig based on the similarity theory[J]. Journal of Vibration and Control, 2016, 22(1): 286?295.

    [15] Z Zerong, W Yongyan, F Zhimin. Similarity analysis between scale model and prototype of large vibrating screen[J]. Shock and Vibration, 2015, 2015: 1?7.

    [16] Moschoudis A P, Tsekouras G J, Kanellos F D. Design of particular electrical machines by using similarity theory and scale factors[C]// 2014 International Conference on Electrical Machines (ICEM). IEEE, 2014.

    [17] 林皋, 朱彤, 林蓓. 結(jié)構(gòu)動(dòng)力模型試驗(yàn)的相似技巧[J]. 大連理工大學(xué)學(xué)報(bào), 2000, 40(1): 4?11. LIN Gao, ZHU Tong, LIN Bei. Similar techniques for structural dynamic model tests[J]. Journal of Dalian University of Technology, 2000, 40(1): 4?11.

    [18] 朱彤. 結(jié)構(gòu)動(dòng)力模型相似問(wèn)題及結(jié)構(gòu)動(dòng)力試驗(yàn)技術(shù)研究[D]. 大連: 大連理工大學(xué)建設(shè)工程學(xué)部, 2004: 67? 75. ZHU Tong. Research on structural dynamic model similarity problems and structural dynamic test technology[D]. Dalian: Dalian University of Technology, Department of Construction Engineering, 2004: 67?75.

    [19] 李增剛, 詹福良. Virtual. Lab Acoustics聲學(xué)仿真計(jì)算高級(jí)應(yīng)用實(shí)例[M]. 北京: 國(guó)防工業(yè)出版社, 2010: 46? 48. LI Zenggang, ZHAN Fuliang. Virtual. Lab Acoustics advanced application example of acoustic simulation[M]. Beijing: National Defense Industry Press, 2010: 46?48.

    [20] 中國(guó)標(biāo)準(zhǔn)出版社第二室. 環(huán)境監(jiān)測(cè)方法標(biāo)準(zhǔn)匯編. 噪聲與振動(dòng)[S]. The second room Standard Press of China. Compilation of standards for environmental monitoring methods. Noise and vibration[S].

    Vibration analysis and noise prediction of 100% low floor gearbox

    PENG Ruitao, XIE Yanpei, LUO Xiuqiang, HU Congfang, LIU Xianghuan

    (School of Mechanical Engineering, Xiangtan University, Xiangtan 411105, China)

    Vibration analysis and noise prediction of 100% low floor gearbox can provide a basis for its design and optimization. A rigid-flexible coupling multi-body dynamics model of 100% low floor gearbox coupled with transmission system and structural system was established, and the bearing reaction force was calculated. The gearbox modal analysis model was established and the vibration characteristics of the housing were studied. A similar model of the housing was designed based on the similarity theory, and an experimental modal analysis was performed to verify the correctness of the finite element model. The bearing reaction force obtained by the rigid-flexible coupled multi-body dynamics simulation was used as the excitation of the housing, and the direct integration method was used to calculate the vibration response of the housing. Then, using the housing vibration velocity as the boundary condition, the direct boundary element method radiation noise analysis model of the gearbox was established, and the sound pressure on the surface of the gearbox and the external sound field radiation noise were predicted. The results show that the noise peak frequency of the gearbox radiated noise frequency domain curve appear at the meshing frequency of the gear pair and its harmonics, as well as the low- order natural frequency of the housing.

    100% low floor; gearbox; rigid-flexible coupling dynamics; vibration characteristics; similarity theory; radiated noise

    TB532

    A

    1672 ? 7029(2020)12 ? 3181 ? 09

    10.19713/j.cnki.43?1423/u.T20200128

    2020?02?21

    國(guó)家自然科學(xué)基金資助項(xiàng)目(51975504,51705442);湖南省自然科學(xué)(株洲)聯(lián)合基金資助項(xiàng)目(2018JJ4082);湖南省教育廳優(yōu)秀青年項(xiàng)目(19B539)

    彭銳濤(1982?),男,湖南衡陽(yáng)人,教授,從事高效精密加工與傳動(dòng)研究;E?mail:pengruitao@xtu. edu.cn

    (編輯 陽(yáng)麗霞)

    猜你喜歡
    模態(tài)有限元振動(dòng)
    振動(dòng)的思考
    振動(dòng)與頻率
    中立型Emden-Fowler微分方程的振動(dòng)性
    國(guó)內(nèi)多模態(tài)教學(xué)研究回顧與展望
    基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識(shí)別
    磨削淬硬殘余應(yīng)力的有限元分析
    UF6振動(dòng)激發(fā)態(tài)分子的振動(dòng)-振動(dòng)馳豫
    由單個(gè)模態(tài)構(gòu)造對(duì)稱(chēng)簡(jiǎn)支梁的抗彎剛度
    基于SolidWorks的吸嘴支撐臂有限元分析
    箱形孔軋制的有限元模擬
    上海金屬(2013年4期)2013-12-20 07:57:18
    99国产精品99久久久久| 每晚都被弄得嗷嗷叫到高潮| 精品视频人人做人人爽| 91av网站免费观看| 亚洲成a人片在线一区二区| 人人妻人人添人人爽欧美一区卜| 国产精品亚洲一级av第二区| 免费在线观看完整版高清| 最近最新中文字幕大全免费视频| 国产男女内射视频| 777久久人妻少妇嫩草av网站| 亚洲国产毛片av蜜桃av| 99国产综合亚洲精品| 国产国语露脸激情在线看| 国产精品一区二区在线不卡| 精品一区二区三区视频在线观看免费 | 视频区欧美日本亚洲| 99riav亚洲国产免费| 亚洲精品国产一区二区精华液| av电影中文网址| 丝瓜视频免费看黄片| 国产欧美亚洲国产| 女警被强在线播放| 亚洲精品av麻豆狂野| 妹子高潮喷水视频| 国产亚洲精品第一综合不卡| 人人妻人人爽人人添夜夜欢视频| 黄色视频,在线免费观看| 人人妻人人澡人人爽人人夜夜| 大香蕉久久网| 夜夜爽天天搞| 视频在线观看一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 熟女少妇亚洲综合色aaa.| 少妇猛男粗大的猛烈进出视频| 激情在线观看视频在线高清 | 亚洲精品中文字幕一二三四区| 天天操日日干夜夜撸| 亚洲av欧美aⅴ国产| 日韩大码丰满熟妇| 18禁黄网站禁片午夜丰满| 99久久综合精品五月天人人| 国产精品永久免费网站| 久久影院123| 在线观看免费高清a一片| 国产精品乱码一区二三区的特点 | 最新在线观看一区二区三区| 校园春色视频在线观看| 少妇 在线观看| 黄色女人牲交| 两性夫妻黄色片| 亚洲在线自拍视频| 乱人伦中国视频| 国产精品一区二区在线观看99| 无人区码免费观看不卡| 国产亚洲精品一区二区www | 国产色视频综合| 99riav亚洲国产免费| 久久草成人影院| 两个人看的免费小视频| 国产精品1区2区在线观看. | 午夜成年电影在线免费观看| 国产亚洲一区二区精品| 亚洲中文av在线| 一进一出抽搐动态| 亚洲一区二区三区不卡视频| 欧美黑人精品巨大| 久久午夜亚洲精品久久| 欧美日韩国产mv在线观看视频| 精品人妻1区二区| 久久国产精品大桥未久av| 免费不卡黄色视频| 美女国产高潮福利片在线看| 亚洲,欧美精品.| 嫁个100分男人电影在线观看| 男女高潮啪啪啪动态图| 丝袜在线中文字幕| 精品久久久久久久毛片微露脸| 国产麻豆69| 香蕉久久夜色| 精品人妻在线不人妻| 一级片免费观看大全| 多毛熟女@视频| 亚洲欧美日韩另类电影网站| 99在线人妻在线中文字幕 | 美女午夜性视频免费| 国产成人精品久久二区二区91| 国产精品久久电影中文字幕 | 狠狠婷婷综合久久久久久88av| 久久精品国产99精品国产亚洲性色 | 国产成人精品久久二区二区91| 国产精品久久电影中文字幕 | 男女之事视频高清在线观看| 日本vs欧美在线观看视频| 在线观看免费午夜福利视频| 免费在线观看完整版高清| 免费高清在线观看日韩| 18禁美女被吸乳视频| 一个人免费在线观看的高清视频| 亚洲九九香蕉| 国产精品久久视频播放| 国产精品秋霞免费鲁丝片| 亚洲avbb在线观看| 亚洲中文日韩欧美视频| 一级a爱视频在线免费观看| 国产精品影院久久| 动漫黄色视频在线观看| videos熟女内射| 在线天堂中文资源库| 两性夫妻黄色片| 久久 成人 亚洲| 人人妻,人人澡人人爽秒播| 免费女性裸体啪啪无遮挡网站| 老司机福利观看| 亚洲五月色婷婷综合| 99热只有精品国产| 极品教师在线免费播放| 欧美黄色淫秽网站| 777久久人妻少妇嫩草av网站| 1024香蕉在线观看| 免费不卡黄色视频| 精品少妇一区二区三区视频日本电影| 欧美黑人欧美精品刺激| 久久久精品国产亚洲av高清涩受| 在线观看66精品国产| 99国产精品免费福利视频| 人人妻人人澡人人看| av电影中文网址| 午夜福利乱码中文字幕| 女警被强在线播放| 国产精品一区二区精品视频观看| 美女高潮到喷水免费观看| 国产成人免费观看mmmm| 变态另类成人亚洲欧美熟女 | ponron亚洲| 俄罗斯特黄特色一大片| av一本久久久久| 身体一侧抽搐| 18禁美女被吸乳视频| 高清在线国产一区| 日韩欧美国产一区二区入口| 最近最新免费中文字幕在线| 丝袜美足系列| 午夜免费观看网址| 夜夜夜夜夜久久久久| 国产精品国产高清国产av | 在线国产一区二区在线| 国产免费av片在线观看野外av| 男人操女人黄网站| 久久久久国产精品人妻aⅴ院 | 99久久99久久久精品蜜桃| 国产又爽黄色视频| 99香蕉大伊视频| 美国免费a级毛片| 国产高清激情床上av| 如日韩欧美国产精品一区二区三区| 少妇猛男粗大的猛烈进出视频| 国产成人系列免费观看| 麻豆国产av国片精品| 精品久久久精品久久久| 国产精品国产av在线观看| 亚洲人成电影免费在线| 一本大道久久a久久精品| av视频免费观看在线观看| 69精品国产乱码久久久| 国产精品98久久久久久宅男小说| 国产一区二区三区在线臀色熟女 | 老司机在亚洲福利影院| 免费高清在线观看日韩| 欧美久久黑人一区二区| 婷婷精品国产亚洲av在线 | 九色亚洲精品在线播放| 午夜福利在线免费观看网站| 最新美女视频免费是黄的| 国产欧美日韩一区二区三区在线| 69av精品久久久久久| 亚洲欧美日韩高清在线视频| 亚洲精品久久成人aⅴ小说| 久9热在线精品视频| 美国免费a级毛片| 丰满的人妻完整版| 国产亚洲欧美98| 亚洲五月色婷婷综合| 久久久国产精品麻豆| 国产精品久久电影中文字幕 | 国产成人影院久久av| 在线国产一区二区在线| 一本一本久久a久久精品综合妖精| 很黄的视频免费| 黄色成人免费大全| 欧美激情高清一区二区三区| 欧美最黄视频在线播放免费 | 精品乱码久久久久久99久播| 777米奇影视久久| 妹子高潮喷水视频| 香蕉丝袜av| 91大片在线观看| 老熟女久久久| 制服诱惑二区| 欧美激情久久久久久爽电影 | 国产不卡av网站在线观看| 久久久精品国产亚洲av高清涩受| 又紧又爽又黄一区二区| 叶爱在线成人免费视频播放| 免费观看精品视频网站| 另类亚洲欧美激情| 亚洲人成伊人成综合网2020| 又紧又爽又黄一区二区| 亚洲av日韩在线播放| 91大片在线观看| 国产亚洲精品第一综合不卡| 在线观看免费午夜福利视频| 天堂俺去俺来也www色官网| 黄色视频,在线免费观看| 国产人伦9x9x在线观看| 久久国产乱子伦精品免费另类| 一区二区三区精品91| 好看av亚洲va欧美ⅴa在| 欧美激情 高清一区二区三区| av片东京热男人的天堂| 国产高清视频在线播放一区| 天天操日日干夜夜撸| 看片在线看免费视频| 久久狼人影院| 最新在线观看一区二区三区| 国产无遮挡羞羞视频在线观看| 午夜亚洲福利在线播放| 男男h啪啪无遮挡| 80岁老熟妇乱子伦牲交| 别揉我奶头~嗯~啊~动态视频| 黄色毛片三级朝国网站| 一进一出抽搐gif免费好疼 | 黄片小视频在线播放| 啪啪无遮挡十八禁网站| 无限看片的www在线观看| 亚洲成av片中文字幕在线观看| 欧美一级毛片孕妇| 久久久久久人人人人人| 国产精品成人在线| 免费在线观看视频国产中文字幕亚洲| 欧美日本中文国产一区发布| 中文字幕最新亚洲高清| avwww免费| 精品久久久久久电影网| 成人特级黄色片久久久久久久| 精品国产乱码久久久久久男人| 国产在线精品亚洲第一网站| 人妻 亚洲 视频| 精品久久久久久电影网| ponron亚洲| 久久精品91无色码中文字幕| 亚洲午夜精品一区,二区,三区| 露出奶头的视频| 丰满迷人的少妇在线观看| 搡老岳熟女国产| 精品熟女少妇八av免费久了| 亚洲五月色婷婷综合| 欧美日韩乱码在线| 天天躁日日躁夜夜躁夜夜| 亚洲久久久国产精品| 亚洲avbb在线观看| 天天影视国产精品| 国产成人免费观看mmmm| 亚洲av成人不卡在线观看播放网| 久9热在线精品视频| 国产精品一区二区精品视频观看| 亚洲一码二码三码区别大吗| 欧美乱色亚洲激情| 女性生殖器流出的白浆| xxxhd国产人妻xxx| 99精品欧美一区二区三区四区| 91麻豆av在线| 天天躁狠狠躁夜夜躁狠狠躁| 两人在一起打扑克的视频| 丝袜在线中文字幕| 香蕉国产在线看| 九色亚洲精品在线播放| 精品福利永久在线观看| 亚洲国产精品合色在线| 久久人妻熟女aⅴ| 一级毛片高清免费大全| 啦啦啦 在线观看视频| 欧美精品av麻豆av| 黄色女人牲交| 欧美国产精品一级二级三级| avwww免费| 在线av久久热| 国产男靠女视频免费网站| 一区在线观看完整版| 一级a爱视频在线免费观看| 国产伦人伦偷精品视频| 国产成+人综合+亚洲专区| 国产三级黄色录像| 99精品在免费线老司机午夜| 国产欧美日韩精品亚洲av| 男女高潮啪啪啪动态图| 午夜福利一区二区在线看| 久久精品国产99精品国产亚洲性色 | 免费不卡黄色视频| 午夜免费成人在线视频| 热99re8久久精品国产| 亚洲精品国产区一区二| 成年人黄色毛片网站| av免费在线观看网站| 一夜夜www| 亚洲国产欧美网| 国产一区二区三区综合在线观看| 精品少妇久久久久久888优播| 另类亚洲欧美激情| 欧美精品av麻豆av| 18禁国产床啪视频网站| 最新美女视频免费是黄的| 一区二区日韩欧美中文字幕| 欧美中文综合在线视频| 久久久精品免费免费高清| 91麻豆av在线| 精品一区二区三卡| 十分钟在线观看高清视频www| 999久久久国产精品视频| 久久性视频一级片| 激情视频va一区二区三区| 最近最新免费中文字幕在线| 黄色成人免费大全| 亚洲一区二区三区欧美精品| 国产亚洲av高清不卡| 一级片免费观看大全| 岛国在线观看网站| cao死你这个sao货| 国产精品一区二区精品视频观看| 中亚洲国语对白在线视频| 女性被躁到高潮视频| 91av网站免费观看| 久久精品亚洲熟妇少妇任你| 满18在线观看网站| 亚洲国产精品合色在线| 国产精品久久久久久人妻精品电影| 在线十欧美十亚洲十日本专区| 中文字幕色久视频| 日韩制服丝袜自拍偷拍| 色婷婷av一区二区三区视频| 国产野战对白在线观看| 亚洲精品乱久久久久久| 一级a爱片免费观看的视频| 亚洲免费av在线视频| av线在线观看网站| av片东京热男人的天堂| 日韩成人在线观看一区二区三区| 狂野欧美激情性xxxx| videos熟女内射| 免费不卡黄色视频| 亚洲国产毛片av蜜桃av| 一级片'在线观看视频| 国产真人三级小视频在线观看| 可以免费在线观看a视频的电影网站| 精品亚洲成国产av| 久久久国产成人精品二区 | 狠狠婷婷综合久久久久久88av| 欧美乱码精品一区二区三区| 精品免费久久久久久久清纯 | 亚洲人成电影观看| 妹子高潮喷水视频| 欧美日韩一级在线毛片| 日本黄色视频三级网站网址 | 国产乱人伦免费视频| 另类亚洲欧美激情| 少妇粗大呻吟视频| 大型av网站在线播放| 在线观看免费高清a一片| 日韩欧美一区二区三区在线观看 | xxxhd国产人妻xxx| 日韩欧美一区视频在线观看| 人人妻人人澡人人看| 亚洲欧美精品综合一区二区三区| 叶爱在线成人免费视频播放| xxxhd国产人妻xxx| 亚洲视频免费观看视频| 首页视频小说图片口味搜索| 午夜视频精品福利| 久久亚洲真实| 国产在视频线精品| 黄色视频,在线免费观看| 正在播放国产对白刺激| 啦啦啦免费观看视频1| 免费在线观看亚洲国产| 久久婷婷成人综合色麻豆| 国产精品av久久久久免费| 久久国产精品人妻蜜桃| 亚洲精品在线观看二区| 中国美女看黄片| 在线国产一区二区在线| 亚洲中文av在线| 校园春色视频在线观看| 日韩免费av在线播放| 桃红色精品国产亚洲av| 亚洲国产精品合色在线| 午夜福利,免费看| 如日韩欧美国产精品一区二区三区| 王馨瑶露胸无遮挡在线观看| 男女床上黄色一级片免费看| 校园春色视频在线观看| 三上悠亚av全集在线观看| 久久亚洲真实| 曰老女人黄片| 狠狠狠狠99中文字幕| 久久精品国产综合久久久| 亚洲精品美女久久久久99蜜臀| 老司机午夜福利在线观看视频| 村上凉子中文字幕在线| 亚洲色图 男人天堂 中文字幕| 日日爽夜夜爽网站| 国产精品免费一区二区三区在线 | 午夜日韩欧美国产| 老汉色av国产亚洲站长工具| 成年女人毛片免费观看观看9 | 日韩三级视频一区二区三区| 国产成人啪精品午夜网站| 黄色 视频免费看| 国产成人影院久久av| 国产成+人综合+亚洲专区| 高清视频免费观看一区二区| 99热国产这里只有精品6| 亚洲专区国产一区二区| 国产成人精品久久二区二区91| 99精品在免费线老司机午夜| 精品视频人人做人人爽| 国产欧美日韩一区二区精品| 热99国产精品久久久久久7| 国产精品二区激情视频| 久99久视频精品免费| 岛国在线观看网站| 国产精品九九99| 欧洲精品卡2卡3卡4卡5卡区| 中文字幕制服av| 中国美女看黄片| 国产精品久久久久成人av| 国产一区二区激情短视频| 我的亚洲天堂| 精品人妻在线不人妻| 国产亚洲欧美98| 精品福利观看| 亚洲中文av在线| 丰满人妻熟妇乱又伦精品不卡| 精品亚洲成国产av| 免费在线观看黄色视频的| 精品国产亚洲在线| 免费看十八禁软件| 国产精品乱码一区二三区的特点 | 中文字幕av电影在线播放| 色播在线永久视频| 欧美黑人欧美精品刺激| www.熟女人妻精品国产| av欧美777| 18禁裸乳无遮挡动漫免费视频| 国产精品综合久久久久久久免费 | a级毛片黄视频| 老司机影院毛片| 国产激情欧美一区二区| 国产欧美日韩精品亚洲av| 91大片在线观看| 精品一品国产午夜福利视频| 精品欧美一区二区三区在线| 午夜福利视频在线观看免费| 亚洲熟女精品中文字幕| 9热在线视频观看99| 啪啪无遮挡十八禁网站| 亚洲精品一卡2卡三卡4卡5卡| 中国美女看黄片| 一区二区三区激情视频| a在线观看视频网站| 亚洲色图综合在线观看| 热99久久久久精品小说推荐| 欧美另类亚洲清纯唯美| 国产精华一区二区三区| 午夜激情av网站| 亚洲精华国产精华精| 亚洲精品成人av观看孕妇| 国产精华一区二区三区| 久久久久久免费高清国产稀缺| 丝袜人妻中文字幕| 色老头精品视频在线观看| 51午夜福利影视在线观看| 国产精品国产高清国产av | 50天的宝宝边吃奶边哭怎么回事| 午夜福利乱码中文字幕| 中出人妻视频一区二区| 99久久人妻综合| 午夜精品在线福利| 国产精品久久久久久精品古装| av福利片在线| 婷婷精品国产亚洲av在线 | 成年动漫av网址| 国产精品一区二区在线不卡| 亚洲久久久国产精品| 看片在线看免费视频| 99国产综合亚洲精品| 亚洲av片天天在线观看| av有码第一页| 亚洲成人手机| 黑丝袜美女国产一区| 国产成+人综合+亚洲专区| 日本黄色日本黄色录像| 国产99久久九九免费精品| videosex国产| 精品高清国产在线一区| 国产一区二区三区综合在线观看| 免费观看精品视频网站| a级片在线免费高清观看视频| 欧美日韩亚洲国产一区二区在线观看 | 变态另类成人亚洲欧美熟女 | av线在线观看网站| 无人区码免费观看不卡| 老熟女久久久| 很黄的视频免费| 精品亚洲成国产av| 18禁裸乳无遮挡动漫免费视频| 午夜福利视频在线观看免费| 久久久久国产一级毛片高清牌| 天天影视国产精品| 久久久国产一区二区| 亚洲片人在线观看| 欧美中文综合在线视频| 精品人妻在线不人妻| 在线播放国产精品三级| 国产精品免费一区二区三区在线 | 久久久久久久国产电影| 99久久国产精品久久久| 亚洲少妇的诱惑av| 亚洲一区高清亚洲精品| 免费黄频网站在线观看国产| 国产单亲对白刺激| 啦啦啦免费观看视频1| 大码成人一级视频| 日韩免费av在线播放| 两人在一起打扑克的视频| 新久久久久国产一级毛片| 十分钟在线观看高清视频www| 日韩精品免费视频一区二区三区| √禁漫天堂资源中文www| 亚洲片人在线观看| 欧美中文综合在线视频| 国产真人三级小视频在线观看| 免费久久久久久久精品成人欧美视频| 国产免费av片在线观看野外av| 国产精品永久免费网站| 国产av一区二区精品久久| 国产精品永久免费网站| 精品国产美女av久久久久小说| svipshipincom国产片| 少妇被粗大的猛进出69影院| 欧美中文综合在线视频| 欧美黑人精品巨大| 亚洲精品国产色婷婷电影| 成年动漫av网址| 色尼玛亚洲综合影院| 丰满的人妻完整版| 国产99久久九九免费精品| 欧美日韩中文字幕国产精品一区二区三区 | 中文字幕色久视频| 国产野战对白在线观看| 亚洲第一欧美日韩一区二区三区| 91字幕亚洲| 两人在一起打扑克的视频| 日本欧美视频一区| 丰满迷人的少妇在线观看| 国产午夜精品久久久久久| 韩国av一区二区三区四区| 黄色视频,在线免费观看| 午夜福利影视在线免费观看| 一级,二级,三级黄色视频| 国产亚洲av高清不卡| 久久精品亚洲av国产电影网| 每晚都被弄得嗷嗷叫到高潮| 女警被强在线播放| 啦啦啦免费观看视频1| 国内久久婷婷六月综合欲色啪| 1024香蕉在线观看| www.熟女人妻精品国产| 国产精品免费一区二区三区在线 | 亚洲精品在线美女| 69精品国产乱码久久久| 女人被躁到高潮嗷嗷叫费观| 日本一区二区免费在线视频| 免费看a级黄色片| 久久久久国产一级毛片高清牌| 亚洲精品久久成人aⅴ小说| 欧美久久黑人一区二区| 国产精品一区二区在线观看99| 免费在线观看视频国产中文字幕亚洲| 香蕉丝袜av| 成年版毛片免费区| 1024视频免费在线观看| 建设人人有责人人尽责人人享有的| 99久久国产精品久久久| 日韩精品免费视频一区二区三区| 国产aⅴ精品一区二区三区波| 99精品欧美一区二区三区四区| 成人精品一区二区免费| 国产成人精品无人区| 欧美日韩黄片免| 在线av久久热| 精品一品国产午夜福利视频| 成人18禁高潮啪啪吃奶动态图| www.999成人在线观看| 丝袜美足系列| 亚洲av日韩在线播放| 老汉色∧v一级毛片| 老司机午夜福利在线观看视频| 在线免费观看的www视频| 国产日韩一区二区三区精品不卡| 亚洲精品美女久久av网站| netflix在线观看网站| 亚洲av第一区精品v没综合| 欧美成人午夜精品| 极品人妻少妇av视频| 啦啦啦免费观看视频1|