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

    點(diǎn)熄火工況下動(dòng)力總成懸置隔振性能的優(yōu)化方法

    2025-08-15 00:00:00殷智宏林爍上官利堅(jiān)
    重慶大學(xué)學(xué)報(bào) 2025年7期
    關(guān)鍵詞:瞬態(tài)固有頻率不確定性

    中圖分類號(hào):U464.13;U461.4 文獻(xiàn)標(biāo)志碼:A 文章編號(hào):1000-582X(2025)07-050-12

    Optimization methods to enhance vibration isolation performance of a powertrain mounting system at key-on/off conditions

    YIN Zhihong',LIN Shuo',SHANGGUAN Lijian2 (1.School of Mechanical and Automotive Engineering,South China University of Technology,Guangzhou 510640,P.R. China; 2. Ningbo Tuopu Group Co., Ltd., Ningbo, Zhejiang 31580o,P.R. China)

    Abstract:Hybrid electric vehicles often experiences key-on/off conditions. Increasing vibration isolation performance of the powertrain mounting system (PMS)under key-on/off conditions is one of the important measures to enhance NVH(noise,vibrationand harshness)performance of the vehicle.A model with 13 DOFs for PMS was established,and acceleration of the mounts at active side was estimated under the key-on/offcondition. Also,vibration acceleration of the mounts at positive side,passive side and the seat rail were measured.The comparison between calculations and experiments validated the proposed model. Deterministic and uncertain optimization methods for reducing the vibration dose value of longitudinal acceleration of the powertrain were proposed,and the mount parameters were the design variable.Keyon/of experiments were carried out to validate the proposed methods.The experimental results show that the two optimization methods can effectively enhance the vibration performance of the PMS,and the uncertain optimization method can make the PMS more robust and effectively deal with the influence of parameter uncertainty.

    Keywords: powertrain mounting system;key-on/of;transient vibration performance;analysis and optimization; interval uncertainty

    混合動(dòng)力汽車(hybrid electricvehicle,HEV)的能量控制策略會(huì)引起發(fā)動(dòng)機(jī)點(diǎn)熄火工況頻繁出現(xiàn),甚至?xí)诔藛T毫無預(yù)期的情況下發(fā)生。因此,HEV對(duì)點(diǎn)熄火工況下的NVH(noise,vibrationand harshness)要求更加嚴(yán)格。

    Rao等發(fā)現(xiàn)在點(diǎn)熄火時(shí),防扭拉桿和變速器懸置是座椅縱向振動(dòng)的主要傳遞路徑。郭榮等2在某HEV的點(diǎn)熄火噪聲研究中,發(fā)現(xiàn)動(dòng)力總成懸置系統(tǒng)(powertrain mounting system,PMS)的剛體模態(tài)與汽車的傳動(dòng)模態(tài)在頻率上接近,造成汽車NVH性能的惡化。Liu等通過仿真證明了傳動(dòng)系統(tǒng)也是HEV在點(diǎn)熄火時(shí)的重要傳遞路徑。王博等和Hwang通過優(yōu)化點(diǎn)火時(shí)發(fā)動(dòng)機(jī)的扭矩,減小了座椅導(dǎo)軌的振動(dòng)。

    PMS是將發(fā)動(dòng)機(jī)的振動(dòng)傳遞到車身的重要路徑。Bang等通過優(yōu)化使PMS的縱向模態(tài)與俯仰模態(tài)充分解耦,改善了汽車點(diǎn)熄火時(shí)的振動(dòng)。Lee提出座椅導(dǎo)軌的振動(dòng)是評(píng)價(jià)點(diǎn)熄火振動(dòng)的重要指標(biāo),且提高懸置的縱向剛度可減小點(diǎn)熄火時(shí)動(dòng)力總成的振動(dòng),從而減小車身的振動(dòng)。Wang等?利用液壓懸置在低頻時(shí)提供大阻尼來降低點(diǎn)熄火時(shí)的振動(dòng)。Xiao等以動(dòng)力總成縱向加速度作為瞬態(tài)振動(dòng)的評(píng)價(jià)指標(biāo)之一。Qin等利用振動(dòng)劑量值(vibration dose value,VDV)來描述瞬態(tài)振動(dòng)的強(qiáng)度。康強(qiáng)等提出,與扭矩軸(torque rolaxis,TRA)布置相比,3點(diǎn)質(zhì)心式的懸置布置具有更強(qiáng)的抗扭能力,更適合純電動(dòng)車。徐中明等[2以PMS的解偶率為優(yōu)化目標(biāo)來優(yōu)化設(shè)計(jì)懸置剛度,仿真結(jié)果表明,能量解耦率顯著提高,瞬態(tài)振動(dòng)性能也得到提升。高永峰[13]建立了一個(gè)簡(jiǎn)化的整車有限元模型,提高了PMS瞬態(tài)響應(yīng)的計(jì)算效率。鄧召學(xué)等[14通過控制磁流變懸置的電流,使懸置在發(fā)動(dòng)機(jī)點(diǎn)熄火中獲得大阻尼,成功降低了汽車點(diǎn)熄火的振動(dòng)。

    受加工制造、材料老化、機(jī)械磨損等因素的影響,懸置會(huì)不可避免地存在不確定性因素。呂輝等[1s-1]分別將懸置參數(shù)視為模糊不確定模型和隨機(jī)與區(qū)間混合不確定模型,利用攝動(dòng)中心差分法(perturbation centerdifference method,PCDM)來快速求解PMS的不確定性響應(yīng),并用蒙特卡洛法(monte carlo method,MCM)驗(yàn)證了PCDM計(jì)算精度和效率。Wu等[]利用扭矩軸解耦理論對(duì)目標(biāo)函數(shù)采用六sigma準(zhǔn)則,對(duì)PMS進(jìn)行不確定性優(yōu)化。Cai等[將懸置參數(shù)看作混合隨機(jī)-區(qū)間不確定性模型,并基于PCDM提出了一種不確定性優(yōu)化方法。Zheng等基于舍選抽樣的理論對(duì)客車的懸置系統(tǒng)進(jìn)行了不確定性優(yōu)化。

    文中引人動(dòng)力總成縱向加速度的VDV作為目標(biāo)函數(shù)之一,提出了一種提高點(diǎn)熄火隔振性能的PMS優(yōu)化方法;考慮了懸置參數(shù)具有區(qū)間不確定性,并基于一階PCDM提出了一種PMS的不確定性優(yōu)化方法。

    1點(diǎn)熄火工況的計(jì)算模型

    1.1 定義

    點(diǎn)熄火工況的瞬態(tài)響應(yīng)計(jì)算模型如圖1所示。該模型包括動(dòng)力總成(6個(gè)自由度),車身(垂向、俯仰、側(cè)傾3個(gè)自由度)和4個(gè)簧下質(zhì)量(垂直向自由度),共有13個(gè)自由度。

    定義整車坐標(biāo)系(vehiclecoordinatesystem,VCS)為 Ov-XvYvZv,Xv 指向汽車尾部,Z垂直向上, Yv 由右手定則確定。定義動(dòng)力總成坐標(biāo)系(powertrain coordinate system,PCS)為 Op-XpYpZp 、車身坐標(biāo)系(body coordinatesystem,BCS)為 Ob-XbYbZb 和非簧載質(zhì)量的坐標(biāo)系 Ou,j-Xu,jYu,jZu,j(j=1,2,3,4)? 其中, Xp,Xb,Xu,j 與 Xv 平行; Yp,Yb,Yu,j 與 Yv 平行; Zp,Zb,Zu,j 與 Zv 平行。

    定義第 i 個(gè)懸置的局部坐標(biāo)系為 Om,i-Xm,iYm,iZm,i(i-1,2,3) 。定義 Am,i 為 Om,i–Xm,iYm,iZm,i 到 Ov-XvYvZv 的方向余弦矩陣。試驗(yàn)車的左右懸置,其局部坐標(biāo)系與整車坐標(biāo)系平行,即 Am,1 和 Am,2 為單位矩陣。

    拉桿的方向余弦矩陣與拉桿的車身端和動(dòng)力總成端的中心存在以下關(guān)系。防扭拉桿的實(shí)物圖如圖2所示。

    圖1十三自由度模型的示意圖

    Fig.1 Avehiclemodewith13DOFs

    圖2防扭拉桿的實(shí)物圖

    Fig.2 Aphoto of torque strut

    rt 為拉桿方向矢量,表示如下

    式中: Xm?Ym 和 Zm 為防扭拉桿動(dòng)力總成端中心位置在VCS的坐標(biāo); X??Y? 和 Zb 為防扭拉桿車身端中心位置在VCS的坐標(biāo)。

    當(dāng) Ym 與 Yb 相等時(shí),定義 為 rt 的單位向量, ?Am,3 為防扭拉桿局部坐標(biāo)系到VCS的方向余弦矩陣:

    因此,拉桿的方向余弦矩陣可以由拉桿的車身端和動(dòng)力總成端的位置來定義。

    動(dòng)力總成質(zhì)心的位移、車身質(zhì)心的位移和4個(gè)簧下質(zhì)量的位移可以表示成:

    qp=(xp,γp,zp,αp,βp,γpT,

    在整車坐標(biāo)系中, rpcg 為動(dòng)力總成的質(zhì)心, rm,i 為第 i 個(gè)懸置的位置,第 i 個(gè)懸置在PCS的位置矢量為

    式中,上標(biāo)p分別表示該值是 PCS 。

    當(dāng)位置矢量 rm,ip 為以下取值時(shí),

    對(duì)應(yīng)的位置轉(zhuǎn)移矩陣為

    式中, 為從動(dòng)力總成質(zhì)心到第 i 個(gè)懸置中心的位置轉(zhuǎn)移矩陣。

    1.2瞬態(tài)響應(yīng)計(jì)算模型的運(yùn)動(dòng)微分方程

    十三自由度的瞬態(tài)響應(yīng)計(jì)算模型包括動(dòng)力總成(6個(gè)自由度),車身(垂向、俯仰、側(cè)傾3個(gè)自由度)和4個(gè)簧下質(zhì)量(垂直向自由度)。該十三自由度模型的運(yùn)動(dòng)微分方程為[20]:

    1.3右懸置主動(dòng)端縱向加速度

    發(fā)動(dòng)機(jī)點(diǎn)火可以分為3個(gè)階段。階段I:拖轉(zhuǎn)階段,由起步電機(jī)帶動(dòng)發(fā)動(dòng)機(jī)旋轉(zhuǎn),發(fā)動(dòng)機(jī)轉(zhuǎn)速升高。但此時(shí)發(fā)動(dòng)機(jī)不會(huì)產(chǎn)生扭矩。階段II::點(diǎn)火升速階段,當(dāng)發(fā)動(dòng)機(jī)被帶動(dòng)到一定轉(zhuǎn)速后,發(fā)動(dòng)機(jī)開始點(diǎn)火燃燒并產(chǎn)生扭矩。此時(shí)發(fā)動(dòng)機(jī)轉(zhuǎn)速迅速升高直至怠速轉(zhuǎn)速。該過程扭矩波動(dòng)大,且有沖擊力產(chǎn)生。階段II:轉(zhuǎn)速穩(wěn)定階段,在發(fā)動(dòng)機(jī)升至怠速轉(zhuǎn)速后,扭矩變動(dòng)則趨于穩(wěn)定。圖3為試驗(yàn)車在發(fā)動(dòng)機(jī)啟動(dòng)時(shí)座椅導(dǎo)軌的加速度。由此可見,發(fā)動(dòng)機(jī)點(diǎn)火時(shí)的座椅導(dǎo)軌最大振動(dòng)的方向?yàn)?X 向;且振動(dòng)的峰值發(fā)生在階段II。因此,要提高該車點(diǎn)火時(shí)振動(dòng)性能,其關(guān)鍵在于減小發(fā)動(dòng)機(jī)點(diǎn)火階段ⅡI的縱向振動(dòng)。

    圖3座椅導(dǎo)軌的加速度Fig.3 Seatrailacceleration

    點(diǎn)熄火的振動(dòng)主要來源于升速階段的沖擊扭矩和扭矩的快速波動(dòng)。將點(diǎn)熄火的激勵(lì)力簡(jiǎn)化為

    式中: Mmax 為最大扭矩;k為扭矩加載速率。

    點(diǎn)火時(shí)的轉(zhuǎn)速很低(試驗(yàn)車點(diǎn)火時(shí)的轉(zhuǎn)速為 150r/min ,而怠速轉(zhuǎn)速為 800r/min ),此時(shí)的往復(fù)慣性力較小。與沖擊扭矩相比,是可忽略的。此時(shí)動(dòng)力總成的激振力為

    F?=[0000My0]T°

    利用上述定義的激勵(lì)力和系統(tǒng)參數(shù),由式(7)可計(jì)算得到動(dòng)力總成的加速度向量 。右懸置主動(dòng)端的加速度為

    式中: 為從動(dòng)力總成質(zhì)心到右懸置中心的位置轉(zhuǎn)移矩陣; 為右懸置主動(dòng)端縱向加速度。

    2整車點(diǎn)熄火試驗(yàn)與模型驗(yàn)證

    為了驗(yàn)證模型的精度和優(yōu)化方法的有效性,本研究開展了整車點(diǎn)熄火的試驗(yàn)。

    2.1 測(cè)試方法

    試驗(yàn)車輛PMS的懸置布置方式為3點(diǎn)TRA布置。防扭拉桿安裝角可通過增減其動(dòng)力總成端的墊塊進(jìn)行調(diào)整,如圖4所示。在右(左)懸置的主、被動(dòng)端增減相同厚度的墊塊,可實(shí)現(xiàn)右(左)懸置垂向位置的改變。

    圖5為右懸置位置的調(diào)整。將加速度傳感器分別安裝于3個(gè)懸置的主、被動(dòng)端。此外,座椅導(dǎo)軌的振動(dòng)加速度是評(píng)價(jià)汽車振動(dòng)水平的重要評(píng)價(jià)指標(biāo)。因此,還需在座椅導(dǎo)軌處安裝加速度傳感器。數(shù)據(jù)采集系統(tǒng)為西門子LMSSCADAS數(shù)據(jù)采集儀。

    點(diǎn)熄火試驗(yàn)步驟:將汽車擋位調(diào)至P擋,啟動(dòng)發(fā)動(dòng)機(jī),經(jīng)過點(diǎn)火過程后進(jìn)入怠速狀態(tài),持續(xù)怠速5s后發(fā)動(dòng)機(jī)熄火。最后,對(duì)測(cè)到的加速度時(shí)域信號(hào)進(jìn)行短時(shí)傅里葉變換。變換譜類型為線性自功率譜(autopowerlinear),頻率分辨率為 2Hz 。幅值形式為Peak,短時(shí)傅里葉變換的時(shí)間間隔為 0.1s 。以幅值最大的一次短時(shí)傅里葉變換作為數(shù)據(jù)分析對(duì)象。

    圖4防扭拉桿安裝角的調(diào)整

    Fig.4Adjustion for the torque struts

    圖5右懸置位置的調(diào)整

    2.2 模型的驗(yàn)證

    按2.1節(jié)的測(cè)試方法進(jìn)行點(diǎn)熄火實(shí)驗(yàn),并將測(cè)試結(jié)果與模型計(jì)算結(jié)結(jié)果進(jìn)行比較。以右懸置主動(dòng)端的縱向加速度為例,其測(cè)試結(jié)果與計(jì)算結(jié)果如圖6所示。虛線為發(fā)動(dòng)機(jī)點(diǎn)火過程右懸置主動(dòng)端縱向加速度在頻域下的計(jì)算結(jié)果,而實(shí)線表示實(shí)測(cè)結(jié)果??梢姺抡媲€與實(shí)測(cè)曲線吻合較好,最大的峰值都在 12Hz ,且最大峰值的相對(duì)誤差為 0.867% 。最大誤差的頻率出現(xiàn)在 18Hz ,最大誤差值為 0.698m/s2",仿真誤差在可接受范圍內(nèi),該結(jié)果表明文中所建立的動(dòng)力學(xué)模型比較準(zhǔn)確,可用于后續(xù)的研究和分析。

    Fig.5 Adjustionfortherightmount

    圖6右懸置主動(dòng)端的縱向加速度

    Fig.6Longitudinal accelerationof the active side ofrightmount

    3懸置系統(tǒng)優(yōu)化方法

    已有的PMS優(yōu)化方法主要關(guān)注提高其解耦率。為了提高PMS在點(diǎn)熄火工況下的隔振性能,在此基礎(chǔ)上還引入了點(diǎn)熄火的瞬態(tài)響應(yīng)作為優(yōu)化目標(biāo)。考慮懸置參數(shù)是否存在不確定性,分別提出了2種優(yōu)化方法。

    3.1 設(shè)計(jì)變量

    優(yōu)化方法是對(duì)懸置的參數(shù)進(jìn)行優(yōu)化。懸置的參數(shù)包括:剛度、阻尼、安裝位置和安裝角。橡膠懸置的阻尼很小且變化范圍不大,故不作為設(shè)計(jì)變量。懸置的剛度、安裝位置及安裝角都可作為設(shè)計(jì)變量。

    3.2 約束條件

    為了提高隔振率,動(dòng)力總成懸置系統(tǒng)的各固有頻率應(yīng)盡可能地遠(yuǎn)離激勵(lì)頻段。根據(jù)隔振理論,要求固有頻率小于激勵(lì)頻率的 倍,且要避開車內(nèi)其他零部件與子系統(tǒng)的固有頻率。為了避免不同方向振動(dòng)相互耦合,各解耦率應(yīng)不低于設(shè)定的最小值。表1為該試驗(yàn)車的固有頻率和解耦率的設(shè)計(jì)要求。

    表1固有頻率與解耦率的設(shè)計(jì)要求

    Table1 DesignrequirementsofthePMS

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

    3.3.1 瞬態(tài)振動(dòng)評(píng)價(jià)指標(biāo)

    座椅導(dǎo)軌的加速度是評(píng)價(jià)汽車在瞬態(tài)工況下振動(dòng)水平的重要指標(biāo)。座椅導(dǎo)軌屬于車身的一部分,汽車發(fā)動(dòng)機(jī)點(diǎn)火時(shí),由于離合器還未傳遞動(dòng)力且地面給輪胎提供的反力很小以至于可以忽略。此時(shí),動(dòng)力總成縱向加速度與車身的縱向加速度存在以下關(guān)系:

    式中: 分別為車身和動(dòng)力總成的縱向加速度; Mb 和 Mp 分別為車身和動(dòng)力總成的質(zhì)量。

    由式(11)可知,減小動(dòng)力總成的縱向加速度,可以減小車身的縱向加速度,進(jìn)而減小座椅導(dǎo)軌的縱向加速度。故將動(dòng)力總成的縱向加速度可作為汽車點(diǎn)火過程的振動(dòng)評(píng)價(jià)指標(biāo)[21]。

    由于振動(dòng)劑量值考慮了沖擊信號(hào)幅值和持續(xù)時(shí)間的影響,適合作為沖擊評(píng)價(jià)指標(biāo),為

    式中: a(t) 為未經(jīng)過處理的加速度時(shí)域信號(hào); avDv 為加速度 a(t) 的振動(dòng)劑量值; τ 為采樣時(shí)間。

    文中選用動(dòng)力總成縱向加速度的VDV作為優(yōu)化目標(biāo)函數(shù)之一。采樣時(shí)間為發(fā)動(dòng)機(jī)啟動(dòng)后的前1s。

    3.3.2動(dòng)力總成懸置系統(tǒng)的固有特性

    為了提高PMS的隔振性能,除了對(duì)PMS的固有頻率范圍進(jìn)行限制,PMS的各個(gè)方向也應(yīng)盡可能地解耦。而PMS各個(gè)方向的能量解耦率是衡量PMS解耦程度的重要指標(biāo)。文中選取PMS的能量解耦率作為目標(biāo)函數(shù)的另一部分。

    3.4懸置系統(tǒng)確定性優(yōu)化方法

    提高點(diǎn)熄火工況下PMS隔振性能的優(yōu)化函數(shù)為

    式中: ω1 和 ω2 為2種評(píng)價(jià)指標(biāo)的權(quán)重系數(shù); λi 為第 i 階解耦率的權(quán)重系數(shù); 為動(dòng)力總成縱向加速度在點(diǎn)火后的前1s內(nèi)的振動(dòng)劑量值;上標(biāo)1和 u 為該參數(shù)的下界和上界。

    3.5 不確定性模型分析

    由于制造、裝配、測(cè)試過程中的誤差與材料老化,PMS的參數(shù)會(huì)不可避免地存在不確定性。針對(duì)懸置系統(tǒng)中的不確定參數(shù)缺乏足夠樣本信息的情況,采用區(qū)間模型進(jìn)行描述。當(dāng)懸置參數(shù)存在區(qū)間不確定時(shí),則PMS的響應(yīng)(動(dòng)力總成的固有頻率、解耦率和縱向加速度的VDV)也是區(qū)間變量。參數(shù) Yj 的區(qū)間上界為 Yj+ ,區(qū)間下界為 Yj- ,則參數(shù) Yj 的區(qū)間中點(diǎn) YjM 和區(qū)間半徑 YjR

    由 B 個(gè)區(qū)間參數(shù)組成的區(qū)間向量為 。定義 fi(Y) 和 Di(Y) 為在區(qū)間不確定情形下的PMS第 i 階固有頻率和解耦率函數(shù),而 為在區(qū)間不確定情形下動(dòng)力總成縱向加速度的VDV函數(shù)。為了高效求解 fi(Y),Di(Y) (, i=1,2,…,6) ,和 ,文中提出了蒙特卡洛法和攝動(dòng)中心差分法。

    3.5.1 蒙特卡洛法MCM

    步驟1:對(duì)第1個(gè)區(qū)間參數(shù) Y1 進(jìn)行抽樣,得到 Y1,1 。

    步驟2:重復(fù)步驟1,對(duì)其他區(qū)間參數(shù)進(jìn)行抽樣,得到一組區(qū)間向量 Y 抽樣一次的結(jié)果為

    根據(jù) Y1s 可獲得一組響應(yīng) fi(Y) 0 i=1,2,…,6 )

    步驟3:進(jìn)行了 k 次抽樣后,結(jié)果為

    Ycs=[Y1sY2s…Yks]T°

    進(jìn)而共獲得 k 組固有頻率 !, (h=1,2,…,k) ,最終篩選出最大值和最小值記為 fi+ 和 fi-

    同理,可求得解耦率 Di(Y) 上界和下界,以及動(dòng)力總成縱向加速度的振動(dòng)劑量值 上界和下界。

    3.5.2 攝動(dòng)中心差分法PCDM

    對(duì)第 i 階 (i=1,2,…,6) 固有頻率 fi(Y) 在區(qū)間中點(diǎn) YM 進(jìn)行一階泰勒展開,展開式忽略了一階以上各項(xiàng):

    由于 fi(Y) 是關(guān)于 Y 的隱函數(shù),無法直接求解式(17)偏導(dǎo)數(shù)值。文中引入中心差分法來求解偏導(dǎo)數(shù):

    式中, δYq 為一微小增量。向量 δYq 為第 q 個(gè)元素為微小增量而其他元素都為0的向量 Y 。即:

    δYq=[0,…,δYq,…,0]

    將式(18)代人式(17)后,利用一階攝動(dòng)中心差分法求解第 i 階固有頻率 fi(Y) 的上下界:

    式中:當(dāng)運(yùn)算為“ + ”時(shí),表示固有頻率上界 fi+(Y) ;為“-\"時(shí),表示固有頻率下界 fi-(Y) 0

    利用式(20)計(jì)算得到第 i 階 (i=1,2,…,6) 固有頻率 fi(Y) 的上下界。同理,可求得第 i 階解耦率 Di(Y) 的上界和下界,以及動(dòng)力總成縱向加速度的振動(dòng)劑量值 上界和下界。

    3.6懸置系統(tǒng)不確定性優(yōu)化

    提高點(diǎn)熄火工況下PMS隔振性能的不確定性優(yōu)化函數(shù)可表示為

    式中,上標(biāo)“ + ”和“-\"分別為該參數(shù)的不確定性響應(yīng)上界和下界。

    4 應(yīng)用實(shí)例

    本實(shí)例通過優(yōu)化防扭拉桿的參數(shù)以獲得較好的點(diǎn)熄火振動(dòng)性能。設(shè)計(jì)變量為防扭拉桿的 X 向剛度和動(dòng)力總成端中心的位置。防扭拉桿的 X 向剛度設(shè)計(jì)區(qū)間為 50N/mm , 270N/mm- ,防扭拉桿動(dòng)力總成端中心的垂向位置設(shè)計(jì)區(qū)間為 [-95mm,-55mm] 。為了分析PMS的穩(wěn)健性,需要改變懸置的參數(shù)來評(píng)價(jià)優(yōu)化參數(shù)的有效性。從實(shí)驗(yàn)的角度出發(fā),改變左右懸置的垂向位置具有容易實(shí)現(xiàn)、經(jīng)濟(jì)性和不容易引人誤差等優(yōu)點(diǎn)。故該實(shí)例考慮左右懸置的高度為區(qū)間不確定模型。左懸置高度的取值范圍為 [319mm,329mm] ,右懸置高度的取值范圍為 0

    4.1 PCDM的驗(yàn)證

    為了驗(yàn)證PCDM求解不確定性響應(yīng)的精度和效率,將PCDM求解的結(jié)果與MCM求解的結(jié)果進(jìn)行對(duì)比。MCM的抽樣次數(shù)為100000次。中心差分法的微小變化增量 δYq 取為 0.001mm 。表2為2種方法求解到的固有特性范圍。表3為2種方法求解到的 范圍。與MCM的結(jié)果相比,PCDM求解到的固有頻率和解耦率的最大相對(duì)誤差是俯仰方向的解耦率,誤差值為 0.352% 。 的最大相對(duì)誤差為 3.63% 。PCDM的計(jì)算耗時(shí)為0.1440s,而MCM的計(jì)算耗時(shí)為 2212.3s 。因此,PCDM具有較高的計(jì)算精度和計(jì)算效率的特點(diǎn),可有效地用于后續(xù)研究。

    表3動(dòng)力總成縱向加速度VDV的范圍

    表2PMS的固有特性范圍Table2TherangeofinherentcharacteristicsforPMS

    Table3TheVDVrangeof powertrainlongitudinalacceleration m:s-2

    4.2 優(yōu)化結(jié)果

    根據(jù)式(13)進(jìn)行確定性優(yōu)化,優(yōu)化后的拉桿參數(shù):剛度為 237.9N/mm ,動(dòng)力總成端中心的垂向位置為-80.9mm 。根據(jù)式(21)進(jìn)行不確定性優(yōu)化,優(yōu)化后的拉桿參數(shù):剛度為 239.0N/mm ,動(dòng)力總成端中心的垂向位置為 -87.7mm 。表4為優(yōu)化后動(dòng)力總成縱向加速度的確定性響應(yīng)和不確定性響應(yīng)。2種方法優(yōu)化后得到的固有頻率的區(qū)間上下界都滿足約束要求,且能量解耦率的下界都較高。

    Table 4 The range of inherent characteristics for PMS with optimized parameters

    表5為2種方法優(yōu)化后的PMS固有特性響應(yīng)范圍。確定性優(yōu)化后PMS的確定性瞬態(tài)響應(yīng)最小,其動(dòng)力總成縱向加速度的VDV為 0.1185m/s2 。如果考慮參數(shù)不確定性時(shí),不確定性優(yōu)化后PMS的不確定性瞬態(tài)響應(yīng)表現(xiàn)最好, 的上界為 0.253gt;m/s2 。

    表4優(yōu)化后PMS固有特性的不確定性響應(yīng)

    表5優(yōu)化后動(dòng)力總成縱向加速度

    4.3 實(shí)驗(yàn)驗(yàn)證

    為了驗(yàn)證優(yōu)化結(jié)果的有效性,利用座椅導(dǎo)軌處的縱向加速度來評(píng)價(jià)瞬態(tài)振動(dòng)。按2.1節(jié)的測(cè)試方法,將3組拉桿樣件應(yīng)用到同一輛汽車上并進(jìn)行點(diǎn)熄火實(shí)驗(yàn)。這3組拉桿樣件記為拉桿A、拉桿B、拉桿C。它們的參數(shù)分別對(duì)應(yīng)原參數(shù),確定性優(yōu)化參數(shù),不確定性優(yōu)化參數(shù),如表6所示。

    Table 5 The response of powertrain longitudinal acceleration with optimized parameters m?s-2

    表6拉桿樣件的參數(shù)

    Table6Theparametersfortorque strut

    通過試驗(yàn)獲得了若干組在發(fā)動(dòng)機(jī)點(diǎn)火后的前2s該車座椅導(dǎo)軌縱向加速度的測(cè)試數(shù)據(jù),如圖7~9所示。計(jì)算圖7~9中時(shí)域信號(hào)的振動(dòng)劑量值,如表7所示。其中L0-R14分別代表左懸置不墊高和右懸置被墊高14mm 。而L0-R0和L7-R7分別代表著左右懸置不墊高和同時(shí)被墊高 7mm 。

    當(dāng)不考慮左右懸置的垂向位置為不確定性參數(shù)時(shí),拉桿B組的振動(dòng)性能是最好的。拉桿C組的振動(dòng)性能也很好。當(dāng)考慮左右懸置的垂向位置為不確定性參數(shù)時(shí),拉桿C組的振動(dòng)性能變化幅度比拉桿B組小。在L0-R14和L7-R7中,拉桿C組的振動(dòng)性能是最好的。拉桿C組的振動(dòng)性能上界也最低。故拉桿C組的振動(dòng)性能較高,并且其振動(dòng)性能的穩(wěn)健性也好。實(shí)驗(yàn)結(jié)果證明了文中提出的2種優(yōu)化方法是有效的。

    圖7拉桿A組的座椅導(dǎo)軌縱向加速度

    Fig.7Seatrail longitudinal accelerationwithtorque strutA

    圖8拉桿B組的座椅導(dǎo)軌縱向加速度

    Fig.8Seatrail longitudinal accelerationwith torque strutB

    圖9拉桿C組的座椅導(dǎo)軌縱向加速度

    Fig.9Seatrail longitudinal accelerationwith torque strutC

    表7座椅導(dǎo)軌的縱向加速度的振動(dòng)劑量值

    Table7 VDVof seatrail longitudinalacceleration

    m·s~2

    5結(jié)論

    1)建立了點(diǎn)熄火工況的瞬態(tài)響應(yīng)計(jì)算模型,仿真計(jì)算了右懸置的縱向加速度,并與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比。結(jié)果表明:在最大峰值頻率相同,且最大峰值的相對(duì)誤差為 0.867% ;整個(gè)頻率段的最大誤差值為 0.698m/s2 仿真誤差較小,所建立的模型具有較高的精度。

    2)利用PCDM計(jì)算PMS的響應(yīng),并與MCM計(jì)算的結(jié)果對(duì)比。PCDM求解到的固有頻率和解耦率的最大相對(duì)誤差為 0.35% ;對(duì)于動(dòng)力總成縱向加速度的VDV,其最大相對(duì)誤差為 3.63% 。PCDM的計(jì)算耗時(shí)為0.144 0s,MCM 的計(jì)算耗時(shí)為 2212.3s 。仿真結(jié)果表明,PCDM具有較高的計(jì)算精度和高效的計(jì)算效率。

    3)提出了2種PMS的優(yōu)化方法,分別設(shè)計(jì)了2組防扭拉桿的參數(shù)。仿真和試驗(yàn)均表明確定性優(yōu)化組的瞬態(tài)響應(yīng)最小。不確定性優(yōu)化組的瞬態(tài)響應(yīng)較好,且在左右懸置的垂向位置發(fā)生改變后,其瞬態(tài)響應(yīng)的上界最小。

    參考文獻(xiàn)

    [1]Rao M V,MorthySN,Raghavendran P.NVHanalysis of powertrainstart/stoptransient phenomenon byusing wavelet analysis and time domain transfer path analysis[C]/SAE Technical Paper Series.SAE International,2015.

    [2]GuoR,CaoC,MiY,etal.Experimentalinvestigationoftenoisebrationandharshnessperformancesofange-etended electricvehicle[J].Proceedingsof theInstitutionofMechanicalEnginers,PartD:JournalofAutomobileEnginering,2016, 230(5): 650-663.

    [3]LiuDH,YuHS,ZhangJW.Multibodydynamicsanalysis forthecoupledvibrationsofapowersplit hybrid electricvehicle duringtheenginestarttransition[J].ProceedingsoftheInstitutionofMechanicalEngineers,PartK:JournalofMulti-body Dynamics,2016,230(4):527-540.

    [4]王博,張振東,于海生,等.發(fā)動(dòng)機(jī)起動(dòng)引起的混合動(dòng)力汽車振動(dòng)分析與控制[J].汽車工程,2019,41(2):184-190,205. WangB,ZhangZD,YuHS,etal.Analysisandcontrolonvibrationresutedfromenginestartofhybridelectricvehicle[J]. Automotive Engineering,2019,41(2):184-190,205.(in Chinese)

    [5]Hwang HY.Minimizingseattrackvibration thatiscaused bytheautomatic start/topofanengine ina power-split hybrid electric vehicle[J].Journal ofVibrationand Acoustics,2013,135(6):061007.

    [6]Bang JH,Yoon HW,WonKM.Experimentandsimulation to improve keyon/offvehicle vibrationquality[C]/SAETechnical Paper Series.SAE International,2007.

    [7]LeeY.Studyontransientvibratioresponsecharacteristicsoffrontransersepowertrainatthekeyon/ffC/AETechnical Paper Series. SAE International,2007.

    [8]WangDY,JangM,HeKF,etal.tudyonbratisupriomethodofveiclewitheginesarttopndutoacartstop[J].Mechanical Systemsand Signal Processing,2020,142:106783.

    [9]Xiao B,Wang MH,JiangYF,etal. Vibrationanalysisofanelectric vehiclemountingsystemundertransientshockcoditions [C]//SAE Technical Paper Series.SAE International, 2021.

    [10]QiYC,Tang XL,JiaT,etal.Noiseandvibrationsuppresioninhybrdelectricvehicles:Stateoftheartandchalleges[J]. Renewable and Sustainable Energy Reviews,2020,124:109782.

    [11]康強(qiáng),顧鵬云,左曙光.純電動(dòng)汽車電驅(qū)動(dòng)總成懸置設(shè)計(jì)原則研究[J].汽車工程,2019,41(11):1235-1242. Kang Q,Gu PY,Zuo SG.Aresearch ondesign principles for electricdrive assembly mounts of pure electric vehicles[J]. Automotive Engineering,2019,41(11):1235-1242.(in Chinese)

    [12]徐中明,李曉,劉和平.純電動(dòng)汽車動(dòng)力總成懸置系統(tǒng)的優(yōu)化[J].汽車工程,2012,34(9):806-810,815. Xu Z M,LiX,LiuHP.Optimzationof thepowertrainmountingsystemofelectric vehicle[J].AutomotiveEnginering,2012, 34(9):806-810,815.(in Chinese)

    [13]高永峰.動(dòng)力總成瞬態(tài)響應(yīng)分析整車建模研究[J].重慶大學(xué)學(xué)報(bào),2012,35(S1):25-30. GaoYF.Researchonvehiclemodeling bytransientresponseanalysis[J].JournalofChongqing University,2012,35(S1):25- 30.(in Chinese)

    [14]鄧召學(xué),楊青樺,蔡強(qiáng),等.應(yīng)用于汽車發(fā)動(dòng)機(jī)start/stop模式的磁流變懸置設(shè)計(jì)與分析[J].重慶大學(xué)學(xué)報(bào),2021,44(2): 1-12. Deng ZX,YangQH,CaiQ,etal.Designandanalysisofamagnetorheologicalmountapliedtotheenginestart/stopmode[J]. Journal of Chongqing University,2021,44(2):1-12. (in Chinese)

    [15]呂輝,楊坤,黃曉婷,等.含模糊不確定性的動(dòng)力總成懸置系統(tǒng)固有特性分析與優(yōu)化[J].華南理工大學(xué)學(xué)報(bào)(自然科學(xué)版), 2020,48(10): 105-112. Lü H,Yang K,Huang XT,etal.Inherentcharacteristicsanalysisandoptimizationof powertrainmounting systems involving fuzzy ucertainty[J].Journalof South China UniversityofTechnology (Natural ScienceEdition),2020,48(10):105-12.(in Chinese)

    [16]呂輝,楊坤,區(qū)霆鈞,等.隨機(jī)與區(qū)間混合不確定情形下動(dòng)力總成懸置系統(tǒng)的固有特性分析[J].振動(dòng)與沖擊,2020,39(23): 58-63,77. Lu H,YangK,OuTJ,etal.NaturalcharacteristicsofPMSwithrandomandinterval hybriduncertainties[J].Jouralof Vibration and Shock,2020,39(23):58-63,77.(in Chinese)

    [17]WuJ,LiuXD,ShanYC,etal.Robustnessoptimizationofenginemountingsystembasedon Six Sigmaand torquerolaxis decoupling method[J].Proceedingsof theInstitutionof MechanicalEnginers,PartD:Journalof AutomobileEngineering, 2019,233(4):1047-1060.

    [18]Cai B H,Lu H.An eficient analysis and optimization method forthe powertrain mounting system with hybridrandomand interval uncertainties[J]. Engineering Optimization,2020,52(9):1522-1541.

    [19]Zheng ZZ,BuXJ,HouL,etal.Hybriduncertaintyanalysisandoptimisation basedonprobabilityboxforbus powertrain mounting system[J]. Journal of Engineering Design,2023,34(1): 23-54.

    [20]ShangguanWB,LiuXA,LvZP,etal.Design methodofautomotive powertrainmounting system basedonvibrationand noise limitations of vehicle level[J].Mechanical Systems and Signal Processing,2016,76: 677-695.

    [21]Wang DY,Zhao X Z,Shangguan WB.Design method forapowertrain mounting systemto decrease the vehicle keyon/off vibratios[J]ProceedingsoftheInstitutionofMechanicalEnginers,PartD:JouralofutomobileEnginering,18,232 (9): 1221-1236.

    (編輯 詹燕平)

    猜你喜歡
    瞬態(tài)固有頻率不確定性
    李華:我的商學(xué)院叫“騰訊”
    中國商人(2025年15期)2025-08-19 00:00:00
    環(huán)境不確定性下人機(jī)協(xié)同的知識(shí)創(chuàng)造與企業(yè)創(chuàng)新
    準(zhǔn)確·明確·精確
    氣相換熱器殼程振動(dòng)特性分析與結(jié)構(gòu)優(yōu)化
    愛達(dá)荷國家實(shí)驗(yàn)室測(cè)試高燃耗燃料
    91麻豆精品激情在线观看国产| 国产亚洲av高清不卡| 国产区一区二久久| 午夜精品在线福利| 欧美人与性动交α欧美精品济南到| 久久久精品国产亚洲av高清涩受| 99在线视频只有这里精品首页| 别揉我奶头~嗯~啊~动态视频| 99在线视频只有这里精品首页| 国产极品粉嫩免费观看在线| 久久久国产成人免费| 成年免费大片在线观看| 欧美中文日本在线观看视频| 两性夫妻黄色片| 1024视频免费在线观看| 免费搜索国产男女视频| 日韩av在线大香蕉| 欧美zozozo另类| 俺也久久电影网| 久久精品夜夜夜夜夜久久蜜豆 | 脱女人内裤的视频| 在线观看舔阴道视频| 欧美日本视频| 亚洲精品一卡2卡三卡4卡5卡| 欧美国产精品va在线观看不卡| 国产亚洲精品av在线| 一本大道久久a久久精品| 真人做人爱边吃奶动态| 国产伦一二天堂av在线观看| 亚洲午夜理论影院| 欧美日本视频| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲第一av免费看| 国产主播在线观看一区二区| 久久99热这里只有精品18| avwww免费| 国产精品久久久人人做人人爽| 亚洲第一青青草原| 香蕉av资源在线| 国产91精品成人一区二区三区| 亚洲国产看品久久| 熟女电影av网| 欧美不卡视频在线免费观看 | 国产亚洲精品av在线| 女性被躁到高潮视频| 韩国av一区二区三区四区| 亚洲一区中文字幕在线| 侵犯人妻中文字幕一二三四区| 一边摸一边做爽爽视频免费| 亚洲精品国产区一区二| 婷婷精品国产亚洲av在线| 自线自在国产av| 精品人妻1区二区| 一级a爱片免费观看的视频| 亚洲国产精品999在线| 男女做爰动态图高潮gif福利片| 国产精品久久久久久精品电影 | 国内久久婷婷六月综合欲色啪| 久久热在线av| 脱女人内裤的视频| 中文字幕高清在线视频| 午夜激情av网站| 亚洲一区二区三区色噜噜| 亚洲成人精品中文字幕电影| 国产成人系列免费观看| 两个人免费观看高清视频| 精品国产美女av久久久久小说| 国产成年人精品一区二区| 久久人妻福利社区极品人妻图片| 老鸭窝网址在线观看| 可以在线观看的亚洲视频| 国产亚洲欧美精品永久| 免费在线观看日本一区| 国产成人av激情在线播放| 99国产精品一区二区三区| 亚洲av中文字字幕乱码综合 | 一夜夜www| 久久天躁狠狠躁夜夜2o2o| 午夜激情福利司机影院| 色av中文字幕| 国产1区2区3区精品| 日本熟妇午夜| 身体一侧抽搐| 午夜免费成人在线视频| 亚洲成人国产一区在线观看| 哪里可以看免费的av片| 国内毛片毛片毛片毛片毛片| 夜夜爽天天搞| 欧美人与性动交α欧美精品济南到| 国产精品久久电影中文字幕| 日日爽夜夜爽网站| 国产1区2区3区精品| 成年版毛片免费区| 免费在线观看日本一区| 中文字幕av电影在线播放| 成人国产综合亚洲| 亚洲精品中文字幕在线视频| 亚洲真实伦在线观看| 国产私拍福利视频在线观看| 少妇裸体淫交视频免费看高清 | videosex国产| av视频在线观看入口| 日韩大码丰满熟妇| 在线观看日韩欧美| 狂野欧美激情性xxxx| 欧美中文日本在线观看视频| 黄色毛片三级朝国网站| 日韩中文字幕欧美一区二区| 美女午夜性视频免费| 国产精品自产拍在线观看55亚洲| 午夜精品在线福利| 制服人妻中文乱码| 亚洲五月色婷婷综合| 国产一区在线观看成人免费| 国产亚洲欧美在线一区二区| 一进一出抽搐gif免费好疼| 精品久久久久久,| 99精品欧美一区二区三区四区| 久久久精品欧美日韩精品| 一边摸一边抽搐一进一小说| 午夜老司机福利片| 午夜福利成人在线免费观看| 亚洲在线自拍视频| 超碰成人久久| 不卡一级毛片| 91成人精品电影| 日韩欧美三级三区| 欧美成人性av电影在线观看| 禁无遮挡网站| 久久草成人影院| 国产激情偷乱视频一区二区| 少妇 在线观看| 亚洲av成人不卡在线观看播放网| 亚洲国产精品sss在线观看| 黄片小视频在线播放| 日日干狠狠操夜夜爽| 女警被强在线播放| 久久午夜综合久久蜜桃| 麻豆一二三区av精品| av电影中文网址| 久久久久久人人人人人| 韩国av一区二区三区四区| 成人精品一区二区免费| 黄片播放在线免费| 黑人欧美特级aaaaaa片| 啦啦啦观看免费观看视频高清| 久久久久国产一级毛片高清牌| 精品久久久久久成人av| 久久人妻av系列| 制服人妻中文乱码| 国产av一区在线观看免费| 精品久久久久久成人av| 最近最新中文字幕大全电影3 | 午夜日韩欧美国产| 精品卡一卡二卡四卡免费| 精品电影一区二区在线| 亚洲片人在线观看| 午夜福利免费观看在线| 最近最新免费中文字幕在线| 久久久国产成人精品二区| 中文字幕久久专区| 免费无遮挡裸体视频| 大型黄色视频在线免费观看| 色综合亚洲欧美另类图片| 叶爱在线成人免费视频播放| 十八禁人妻一区二区| 一卡2卡三卡四卡精品乱码亚洲| 国产亚洲精品av在线| 国产午夜福利久久久久久| 一区福利在线观看| 琪琪午夜伦伦电影理论片6080| 青草久久国产| 国产av一区在线观看免费| 亚洲自偷自拍图片 自拍| 男人舔奶头视频| 亚洲精品美女久久久久99蜜臀| 亚洲欧美一区二区三区黑人| 18禁黄网站禁片午夜丰满| 法律面前人人平等表现在哪些方面| 人人妻人人看人人澡| 18禁观看日本| 色综合站精品国产| 国产精品美女特级片免费视频播放器 | 成人一区二区视频在线观看| 黄色成人免费大全| 国产不卡一卡二| 淫秽高清视频在线观看| 麻豆成人午夜福利视频| 国产亚洲欧美精品永久| 亚洲自偷自拍图片 自拍| ponron亚洲| 国产一卡二卡三卡精品| 一个人观看的视频www高清免费观看 | 香蕉av资源在线| 午夜亚洲福利在线播放| 精品无人区乱码1区二区| 麻豆一二三区av精品| 欧美成人午夜精品| 久久婷婷人人爽人人干人人爱| 18禁美女被吸乳视频| 无人区码免费观看不卡| 级片在线观看| 欧美精品啪啪一区二区三区| 免费看a级黄色片| 亚洲av日韩精品久久久久久密| 亚洲色图 男人天堂 中文字幕| 国产av不卡久久| 免费人成视频x8x8入口观看| 中文字幕精品亚洲无线码一区 | 久久久久久九九精品二区国产 | 国产精品精品国产色婷婷| 国产精品美女特级片免费视频播放器 | 免费在线观看完整版高清| 99国产精品一区二区蜜桃av| 欧美日韩亚洲国产一区二区在线观看| 成人永久免费在线观看视频| 韩国精品一区二区三区| 亚洲,欧美精品.| 国产成人精品无人区| 十分钟在线观看高清视频www| 久久人妻av系列| 曰老女人黄片| 成人三级黄色视频| 久久中文字幕人妻熟女| 亚洲成a人片在线一区二区| 亚洲国产日韩欧美精品在线观看 | av欧美777| 久久精品国产亚洲av高清一级| 国内精品久久久久久久电影| 精品高清国产在线一区| 少妇粗大呻吟视频| 99在线人妻在线中文字幕| 国产精华一区二区三区| 国内久久婷婷六月综合欲色啪| 曰老女人黄片| 久久亚洲精品不卡| a级毛片在线看网站| 国产又黄又爽又无遮挡在线| 亚洲最大成人中文| 18禁黄网站禁片免费观看直播| 中文字幕另类日韩欧美亚洲嫩草| aaaaa片日本免费| 亚洲av成人一区二区三| 国产野战对白在线观看| 亚洲七黄色美女视频| 国产精品久久电影中文字幕| 国产乱人伦免费视频| 在线视频色国产色| 一区福利在线观看| 两个人视频免费观看高清| 久久精品成人免费网站| 日本三级黄在线观看| 少妇 在线观看| 免费电影在线观看免费观看| 给我免费播放毛片高清在线观看| 50天的宝宝边吃奶边哭怎么回事| 不卡一级毛片| 美女高潮到喷水免费观看| 1024手机看黄色片| 91九色精品人成在线观看| 美女免费视频网站| 母亲3免费完整高清在线观看| 久久人妻av系列| 91字幕亚洲| 国产精品自产拍在线观看55亚洲| 国产区一区二久久| 亚洲午夜理论影院| 最好的美女福利视频网| 亚洲性夜色夜夜综合| 久久草成人影院| 婷婷亚洲欧美| 99re在线观看精品视频| 久久天躁狠狠躁夜夜2o2o| 天堂影院成人在线观看| 亚洲成人免费电影在线观看| 在线天堂中文资源库| 99久久国产精品久久久| 欧美乱码精品一区二区三区| 丁香欧美五月| 69av精品久久久久久| 伦理电影免费视频| 亚洲人成网站在线播放欧美日韩| АⅤ资源中文在线天堂| www国产在线视频色| 婷婷亚洲欧美| 国产精品一区二区三区四区久久 | 久久精品国产综合久久久| 国产真人三级小视频在线观看| a级毛片在线看网站| 日本在线视频免费播放| 欧美三级亚洲精品| 国产三级黄色录像| 日韩一卡2卡3卡4卡2021年| 日本熟妇午夜| 99热只有精品国产| 欧美人与性动交α欧美精品济南到| 黄色毛片三级朝国网站| 午夜免费成人在线视频| 亚洲无线在线观看| 两人在一起打扑克的视频| 欧美 亚洲 国产 日韩一| 精华霜和精华液先用哪个| 在线观看午夜福利视频| 黄色a级毛片大全视频| 91字幕亚洲| 亚洲国产精品成人综合色| 亚洲在线自拍视频| 亚洲国产精品999在线| 亚洲国产精品成人综合色| 久久久久国产一级毛片高清牌| 女人爽到高潮嗷嗷叫在线视频| 婷婷精品国产亚洲av在线| 琪琪午夜伦伦电影理论片6080| 91麻豆av在线| 女同久久另类99精品国产91| 哪里可以看免费的av片| 久久香蕉激情| 黑人操中国人逼视频| 亚洲国产看品久久| 老汉色∧v一级毛片| av有码第一页| 精品熟女少妇八av免费久了| 国产亚洲精品综合一区在线观看 | 亚洲精品一卡2卡三卡4卡5卡| 亚洲成av片中文字幕在线观看| x7x7x7水蜜桃| 欧美激情久久久久久爽电影| 国产精品久久久久久精品电影 | 国产视频一区二区在线看| 长腿黑丝高跟| 99久久国产精品久久久| 欧美乱码精品一区二区三区| 啦啦啦免费观看视频1| 满18在线观看网站| 最新美女视频免费是黄的| 欧美一区二区精品小视频在线| 免费观看人在逋| 欧美国产日韩亚洲一区| 亚洲一区二区三区色噜噜| 国产黄a三级三级三级人| 亚洲一区高清亚洲精品| 中文资源天堂在线| 一区二区日韩欧美中文字幕| 午夜老司机福利片| 女同久久另类99精品国产91| 一个人免费在线观看的高清视频| 婷婷丁香在线五月| 一区二区三区激情视频| a级毛片a级免费在线| e午夜精品久久久久久久| 久热爱精品视频在线9| 久久久久久国产a免费观看| 成年女人毛片免费观看观看9| 亚洲中文字幕日韩| 黑丝袜美女国产一区| 在线播放国产精品三级| 欧美一区二区精品小视频在线| 亚洲av第一区精品v没综合| 久久久久精品国产欧美久久久| 一进一出抽搐gif免费好疼| 狠狠狠狠99中文字幕| 日日夜夜操网爽| 日韩精品免费视频一区二区三区| 最近在线观看免费完整版| 欧美乱码精品一区二区三区| 亚洲专区国产一区二区| 怎么达到女性高潮| 日韩国内少妇激情av| 欧美精品啪啪一区二区三区| 黄色 视频免费看| 亚洲激情在线av| 亚洲成av人片免费观看| 99精品在免费线老司机午夜| 97碰自拍视频| 国产成人精品久久二区二区91| 99热6这里只有精品| 在线观看舔阴道视频| 精华霜和精华液先用哪个| 大香蕉久久成人网| 久久久久免费精品人妻一区二区 | 午夜a级毛片| 19禁男女啪啪无遮挡网站| 国产野战对白在线观看| 女性生殖器流出的白浆| 成人国语在线视频| 婷婷亚洲欧美| 久久人人精品亚洲av| 成年版毛片免费区| 男人舔奶头视频| 老汉色av国产亚洲站长工具| 老鸭窝网址在线观看| 91麻豆av在线| 手机成人av网站| 亚洲精品久久成人aⅴ小说| а√天堂www在线а√下载| 国产成人欧美在线观看| 中文资源天堂在线| 精品久久久久久,| 亚洲性夜色夜夜综合| aaaaa片日本免费| 色播在线永久视频| www日本黄色视频网| av视频在线观看入口| 国产精品久久电影中文字幕| 天天一区二区日本电影三级| 国产又爽黄色视频| 大型av网站在线播放| 亚洲人成网站高清观看| ponron亚洲| 波多野结衣高清作品| 又大又爽又粗| 黄色a级毛片大全视频| ponron亚洲| 黑丝袜美女国产一区| 国产亚洲欧美98| 国产熟女xx| 美女大奶头视频| 日韩欧美免费精品| 欧美国产日韩亚洲一区| 国产99白浆流出| 亚洲中文av在线| 在线看三级毛片| 搡老熟女国产l中国老女人| 美女大奶头视频| 欧美成人性av电影在线观看| 国产精品 欧美亚洲| 国产精品一区二区免费欧美| 每晚都被弄得嗷嗷叫到高潮| 听说在线观看完整版免费高清| 十分钟在线观看高清视频www| 国产伦在线观看视频一区| www日本黄色视频网| 大香蕉久久成人网| 91九色精品人成在线观看| 在线视频色国产色| 一区二区三区精品91| 久久青草综合色| 亚洲精品在线观看二区| 黄色视频,在线免费观看| 日本撒尿小便嘘嘘汇集6| 欧美日韩瑟瑟在线播放| 精品国产一区二区三区四区第35| 久9热在线精品视频| 欧美性长视频在线观看| 国产黄a三级三级三级人| 欧美一级毛片孕妇| 欧美亚洲日本最大视频资源| 亚洲黑人精品在线| 国产片内射在线| www.自偷自拍.com| 免费看十八禁软件| 中文字幕久久专区| 国产视频内射| 亚洲精品av麻豆狂野| 亚洲美女黄片视频| 精品国产乱码久久久久久男人| 欧美中文日本在线观看视频| 成年版毛片免费区| 亚洲精品美女久久av网站| 非洲黑人性xxxx精品又粗又长| 亚洲午夜理论影院| 欧美在线黄色| 一边摸一边抽搐一进一小说| 最近最新中文字幕大全电影3 | 悠悠久久av| 国产精品自产拍在线观看55亚洲| 欧美丝袜亚洲另类 | 色av中文字幕| av免费在线观看网站| 亚洲精品美女久久av网站| 桃色一区二区三区在线观看| 757午夜福利合集在线观看| 日本精品一区二区三区蜜桃| 日本撒尿小便嘘嘘汇集6| 一本大道久久a久久精品| 国产真人三级小视频在线观看| 青草久久国产| 日日摸夜夜添夜夜添小说| 日本五十路高清| 视频在线观看一区二区三区| 亚洲精品av麻豆狂野| 日韩成人在线观看一区二区三区| 色播亚洲综合网| 精品久久久久久久久久免费视频| 99国产极品粉嫩在线观看| 国产成人精品无人区| 日本一本二区三区精品| 日韩精品青青久久久久久| 最近最新免费中文字幕在线| 欧美绝顶高潮抽搐喷水| av欧美777| 久9热在线精品视频| 国产伦一二天堂av在线观看| 久久99热这里只有精品18| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲自偷自拍图片 自拍| 亚洲国产毛片av蜜桃av| 一区二区三区精品91| 91九色精品人成在线观看| 久久精品国产亚洲av香蕉五月| 中文字幕人妻熟女乱码| 99国产精品一区二区三区| 亚洲,欧美精品.| 99精品久久久久人妻精品| 欧美不卡视频在线免费观看 | 美女午夜性视频免费| 国产成人影院久久av| 可以在线观看毛片的网站| 99国产极品粉嫩在线观看| 99riav亚洲国产免费| 免费在线观看日本一区| 他把我摸到了高潮在线观看| 黄色片一级片一级黄色片| av免费在线观看网站| 色综合站精品国产| 色尼玛亚洲综合影院| 亚洲精品在线美女| 精品国产乱码久久久久久男人| 色尼玛亚洲综合影院| 一级a爱片免费观看的视频| 国产av又大| 成年版毛片免费区| 成在线人永久免费视频| aaaaa片日本免费| 亚洲成人免费电影在线观看| 久久久久国内视频| 国产乱人伦免费视频| 淫秽高清视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 老汉色∧v一级毛片| 母亲3免费完整高清在线观看| 成人精品一区二区免费| 草草在线视频免费看| 国产成人av教育| 色哟哟哟哟哟哟| 欧美成人性av电影在线观看| 丁香六月欧美| 1024香蕉在线观看| 亚洲精品国产一区二区精华液| 国产爱豆传媒在线观看 | 免费电影在线观看免费观看| 欧美丝袜亚洲另类 | 91九色精品人成在线观看| 国产欧美日韩一区二区三| e午夜精品久久久久久久| ponron亚洲| 两性夫妻黄色片| 黄色成人免费大全| 免费女性裸体啪啪无遮挡网站| 精品久久久久久久久久免费视频| avwww免费| 99热6这里只有精品| 一本大道久久a久久精品| 亚洲av成人av| 国产在线精品亚洲第一网站| 日韩精品免费视频一区二区三区| 亚洲欧美激情综合另类| 男人舔女人的私密视频| 国产真实乱freesex| 999久久久国产精品视频| 无人区码免费观看不卡| 国产亚洲av高清不卡| 97碰自拍视频| 国产真实乱freesex| 欧美又色又爽又黄视频| 女人被狂操c到高潮| avwww免费| 亚洲 欧美一区二区三区| 俄罗斯特黄特色一大片| 亚洲一区二区三区不卡视频| 久久精品国产亚洲av高清一级| 久久国产精品影院| 欧美日韩一级在线毛片| 美女大奶头视频| 伦理电影免费视频| 99在线人妻在线中文字幕| netflix在线观看网站| 两个人视频免费观看高清| av免费在线观看网站| 18禁裸乳无遮挡免费网站照片 | 色播在线永久视频| 18禁黄网站禁片免费观看直播| 美女扒开内裤让男人捅视频| 在线永久观看黄色视频| 人人澡人人妻人| 欧美日韩中文字幕国产精品一区二区三区| 免费一级毛片在线播放高清视频| av在线天堂中文字幕| www日本在线高清视频| 国内毛片毛片毛片毛片毛片| 大香蕉久久成人网| 夜夜看夜夜爽夜夜摸| 亚洲五月婷婷丁香| 悠悠久久av| 免费在线观看成人毛片| 人人妻,人人澡人人爽秒播| 女人被狂操c到高潮| 男女做爰动态图高潮gif福利片| 在线播放国产精品三级| 亚洲国产欧洲综合997久久, | 欧美av亚洲av综合av国产av| 亚洲欧美精品综合久久99| 国产主播在线观看一区二区| 国产aⅴ精品一区二区三区波| 国产欧美日韩精品亚洲av| 精品高清国产在线一区| 免费高清在线观看日韩| 欧美绝顶高潮抽搐喷水| 最近最新免费中文字幕在线| 久久久久久免费高清国产稀缺| 操出白浆在线播放| 亚洲精品一区av在线观看| 亚洲av电影不卡..在线观看| 国产精品,欧美在线| 日韩大尺度精品在线看网址|