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

    基于流固耦合的高速列車裙板支架氣動(dòng)載荷響應(yīng)分析

    2016-05-07 02:46:15張永昌徐宇工張樂樂
    鐵道學(xué)報(bào) 2016年6期
    關(guān)鍵詞:裙板氣動(dòng)流場

    張永昌,徐宇工,張樂樂

    (北京交通大學(xué) 機(jī)械與電子控制工程學(xué)院,北京 100044)

    高速列車車廂下部安裝有設(shè)備艙,其兩側(cè)蓋板被稱為裙板[1]。目前已發(fā)現(xiàn)我國高速列車裙板與車身之間的連接支架出現(xiàn)開裂現(xiàn)象(如圖1所示),對列車行駛造成了極大的安全隱患。

    結(jié)構(gòu)斷裂可分為過載斷裂、疲勞斷裂、蠕變斷裂、應(yīng)力腐蝕斷裂等[2]。裙板支架裂紋的產(chǎn)生具有如下特點(diǎn):在列車正常行駛過程中,裙板支架的應(yīng)力遠(yuǎn)小于其靜強(qiáng)度極限;裂紋是在列車運(yùn)行一定里程后產(chǎn)生的;裙板支架沒有顯著的殘余變形。根據(jù)上述特點(diǎn)判斷該裂痕是結(jié)構(gòu)疲勞造成的[3]。影響結(jié)構(gòu)疲勞強(qiáng)度的因素有結(jié)構(gòu)自身?xiàng)l件(幾何外形、工件表面狀態(tài)、材料本質(zhì)等)及外界條件(工作環(huán)境、載荷條件等)[4]。本文針對裙板支架載荷條件進(jìn)行研究。列車行駛過程中,裙板及支架受到的載荷可分為機(jī)械載荷和氣動(dòng)載荷,且氣動(dòng)載荷是主要的影響因素[5]。故本文重點(diǎn)研究裙板所受氣動(dòng)載荷的特點(diǎn)及其對裙板支架結(jié)構(gòu)承載的影響。

    圖1 裙板支架裂痕

    目前對列車所受氣動(dòng)載荷作用規(guī)律以及其承載特性的研究有實(shí)驗(yàn)及數(shù)值模擬兩種。實(shí)驗(yàn)方面,通過實(shí)測得到結(jié)構(gòu)應(yīng)力較大位置的應(yīng)變值,對其進(jìn)行分析能夠推算出結(jié)構(gòu)所受載荷的大小[6]。該方法存在如下不足:只能監(jiān)測有限個(gè)采集點(diǎn)的數(shù)據(jù);無法區(qū)分氣動(dòng)載荷與機(jī)械載荷的作用規(guī)律;針對特定線路測試,適用范圍有限。采用數(shù)值模擬的方法可以彌補(bǔ)上述不足。目前的數(shù)值模擬研究大致分為三種:在構(gòu)件上施加人為假設(shè)的氣動(dòng)載荷進(jìn)行求解[5];對氣動(dòng)載荷進(jìn)行求解,不分析構(gòu)件的響應(yīng)[7-9];對氣動(dòng)載荷進(jìn)行求解,然后將其施加于構(gòu)件進(jìn)行計(jì)算[10]。以上幾種方法能夠從一定程度上反映氣動(dòng)載荷下高速列車局部構(gòu)件的響應(yīng)問題,但都沒有考慮流場與固體場的實(shí)時(shí)相互作用,描述不夠準(zhǔn)確。因此,有必要將流固耦合技術(shù)[11-14]引入針對高速列車局部構(gòu)件所受氣動(dòng)載荷的作用方式和構(gòu)件的結(jié)構(gòu)響應(yīng)的研究。另外,現(xiàn)有的列車空氣動(dòng)力學(xué)研究中,一般只研究列車外流場的流動(dòng)情況[7-10,15],而裙板所受氣動(dòng)載荷來源于裙板內(nèi)外兩側(cè)流場的壓力差,因此在考慮列車外流場的同時(shí)要對設(shè)備艙內(nèi)流場進(jìn)行研究,需要建立包含上述兩部分流場的仿真模型。

    本文討論將流固耦合技術(shù)應(yīng)用于高速列車局部構(gòu)件結(jié)構(gòu)承載分析的方法,建立包含列車外流場和設(shè)備艙內(nèi)流場的仿真模型。以裙板支架的疲勞破壞問題為例,明確裙板所受氣動(dòng)載荷的產(chǎn)生原因及作用方式,進(jìn)行裙板及支架的結(jié)構(gòu)承載仿真,最終得到氣動(dòng)載荷對結(jié)構(gòu)破壞的影響規(guī)律。

    1 流固耦合仿真流程

    流固耦合算法分為強(qiáng)耦合法和弱耦合法。其中,強(qiáng)耦合法對流場與固體場控制方程進(jìn)行聯(lián)合求解,不僅增加了求解難度,而且難以利用現(xiàn)有的CFD和CSD技術(shù),因此在工程實(shí)際中應(yīng)用較少。弱耦合法分別求解兩個(gè)物理場,并且將兩個(gè)物理場的物理信息進(jìn)行交換從而實(shí)現(xiàn)耦合,應(yīng)用較為普遍。弱耦合法又分為單向耦合法和雙向耦合法,單向耦合法僅考慮流場對固體場的作用,雙向耦合法同時(shí)考慮流場和固體場的相互作用。

    本文基于ANSYS workbench平臺(tái),采用雙向弱耦合法,流體部分通過計(jì)算流體動(dòng)力學(xué)軟件CFX進(jìn)行求解,固體部分通過瞬態(tài)有限元分析模塊Transient Structural進(jìn)行求解。計(jì)算過程中依次對流場與固體場的控制方程進(jìn)行求解,這些控制方程由質(zhì)量守恒定律、動(dòng)量守恒定律等導(dǎo)出。實(shí)現(xiàn)耦合求解的方法為:在流固交界面上交換流體控制方程和固體控制方程的求解數(shù)據(jù)。交換過程遵循以下原則:在流固交界面上,流體與固體的作用力、位移等變量相等。由于CFD計(jì)算與CSD計(jì)算使用的網(wǎng)格不同,流體與固體網(wǎng)格節(jié)點(diǎn)很難重合,在傳遞節(jié)點(diǎn)力時(shí)利用映射算子先將流體節(jié)點(diǎn)應(yīng)力映射到固體節(jié)點(diǎn),然后通過積分求出固體節(jié)點(diǎn)力。另外,固體位移會(huì)導(dǎo)致計(jì)算域變形,采用ALE方法描述網(wǎng)格節(jié)點(diǎn)的運(yùn)動(dòng),并應(yīng)用修正的彈性方程對網(wǎng)格進(jìn)行更新[16]。

    對于本文所研究的具體問題,其求解流程如圖2所示。分別建立流場及固體場仿真模型,在兩者的交界面建立流固耦合面。流場與固體場分別求解,在求解時(shí)數(shù)據(jù)傳遞過程為:在流固耦合面,通過差值將流場計(jì)算所得氣動(dòng)載荷數(shù)據(jù)傳遞到裙板內(nèi)外表面的固體場網(wǎng)格,固體場計(jì)算得到的裙板產(chǎn)生的位移、形變數(shù)據(jù)傳遞到流場網(wǎng)格。最后,分析氣動(dòng)載荷以及結(jié)構(gòu)承載的特性,得到前者對后者的影響規(guī)律。

    圖2 流固耦合仿真流程

    為了保證耦合的準(zhǔn)確性,本文對流場及固體場求解進(jìn)行了統(tǒng)一的時(shí)間和空間控制。時(shí)間控制方面,兩者采用相同的求解時(shí)長和時(shí)間步長。每個(gè)時(shí)間步內(nèi),流場與固體場求解器進(jìn)行數(shù)據(jù)傳遞,傳遞過程進(jìn)行多次迭代。每個(gè)數(shù)據(jù)傳遞迭代步內(nèi),流場與固體場分別進(jìn)行求解,兩者的求解器分別進(jìn)行迭代計(jì)算。反復(fù)迭代直至收斂,然后進(jìn)入下一時(shí)間步??臻g控制方面,流場與固體場采用統(tǒng)一的坐標(biāo)系,保證兩者相對位置相同,并在兩者交界面(即裙板內(nèi)外表面)處采用相近的網(wǎng)格尺寸,從而降低數(shù)據(jù)傳遞過程中可能產(chǎn)生的誤差。

    2 仿真模型的建立

    本文所建立的流固耦合仿真模型如圖3所示,該模型分為固體場及流場兩部分。

    圖3 流固耦合仿真模型

    2.1 固體場模型

    固體場模型包括列車、設(shè)備艙、裙板及支架模型。其中,列車、設(shè)備艙模型用于為流場模型劃定壁面邊界,裙板及支架模型用于結(jié)構(gòu)動(dòng)載荷計(jì)算。實(shí)際列車結(jié)構(gòu)十分復(fù)雜,需要進(jìn)行簡化,忽略固體場對流場影響較小的結(jié)構(gòu)。

    (1)列車模型

    如圖3(a)所示,由于列車中間部分截面不變,縮短的模型對列車流場結(jié)構(gòu)改變較小[17],本文采用3節(jié)車編組的CRH3型高速列車,全車長76 m。去掉車身的細(xì)小結(jié)構(gòu)以及復(fù)雜的轉(zhuǎn)向架結(jié)構(gòu)。

    (2)設(shè)備艙模型

    為了模擬設(shè)備艙內(nèi)部的流場,需要建立車下設(shè)備艙模型。為簡化計(jì)算,只在中間車底部建立設(shè)備艙艙體。如圖3(b)所示,艙體由裙板等構(gòu)成,其兩側(cè)布置有4個(gè)尺寸為670 mm×510 mm四邊形通風(fēng)口。車下設(shè)備簡化為規(guī)則六面體。

    (3)裙板及其支架模型

    以列車后部某裙板為例,其結(jié)構(gòu)如圖3(c)所示。裙板厚6 mm,通過三個(gè)裙板支架與車體相連。所有固體均采用線彈性材料模型以及鋁合金材料,其密度為2 770 kg/m3,楊氏模量為71 GPa,泊松比為0.33。由于支架后端固定于車體且裙板下端與設(shè)備艙底板相互連接,故支架后端和裙板下端為固定約束。裙板內(nèi)外表面為流固耦合面,用以接收氣動(dòng)載荷數(shù)據(jù),同時(shí)將固體位移數(shù)據(jù)傳遞至流體求解器。結(jié)構(gòu)網(wǎng)格如圖4所示,其中裙板網(wǎng)格尺寸為20 mm,支架網(wǎng)格尺寸為3 mm。

    圖4 裙板及支架網(wǎng)格

    2.2 流場模型

    流場模型包括車外流場及設(shè)備艙內(nèi)流場,兩者通過設(shè)備艙通風(fēng)口相連。車外流場模型尺寸如圖3(d)所示。

    使用結(jié)構(gòu)化網(wǎng)格劃分方式,在流動(dòng)復(fù)雜的區(qū)域,如車頭、車尾以及設(shè)備艙附近進(jìn)行加密。在車體表面、設(shè)備艙內(nèi)外表面及地面處生成邊界層網(wǎng)格,為配合壁面函數(shù)的使用,設(shè)置邊界層網(wǎng)格的無量綱厚度y+>30。最小網(wǎng)格尺寸為50 mm,網(wǎng)格數(shù)量為445萬,網(wǎng)格無關(guān)性驗(yàn)證見2.3節(jié)。列車表面網(wǎng)格見圖5。

    圖5 列車表面網(wǎng)格

    采用高精度差分格式對控制方程進(jìn)行離散。由于馬赫數(shù)小于0.3,按不可壓縮流動(dòng)處理。由于仿真同時(shí)涉及遠(yuǎn)離壁面的車外流場以及靠近壁面的設(shè)備艙內(nèi)流場,本文采用在這兩種區(qū)域都有良好適應(yīng)性的SST湍流模型。在近壁面區(qū)域使用標(biāo)準(zhǔn)壁面函數(shù)。

    邊界條件的設(shè)置如圖3(d)所示:計(jì)算域前端為速度入口,入口速度為列車行駛速度;后端為壓力出口,壓力取101 325 Pa;頂端及兩側(cè)為對稱邊界;車身外表面以及設(shè)備艙內(nèi)各種結(jié)構(gòu)均為壁面邊界;地面為移動(dòng)壁面邊界,以來流速度移動(dòng),從而模擬地面氣動(dòng)效應(yīng)。裙板內(nèi)外表面為流固耦合面,其表面網(wǎng)格節(jié)點(diǎn)的位置由固體計(jì)算所得位移數(shù)據(jù)確定。由于裙板變形較小,故無需進(jìn)行大范圍的網(wǎng)格重構(gòu),僅改變?nèi)拱甯浇W(wǎng)格的位置,節(jié)點(diǎn)位移值由式( 1 )確定。

    ( 1 )

    式中:δ為網(wǎng)格節(jié)點(diǎn)的位移值;Vc為網(wǎng)格體積;Cs為常數(shù)。

    一般認(rèn)為,若數(shù)值模型的庫朗數(shù)(uΔt/Δx,其中,u為空氣流速;Δt為時(shí)間步長;Δx為網(wǎng)格尺寸)小于1,則計(jì)算收斂性較好。本文選取時(shí)間步長為4×10-4s,對于不同的車速,模型的庫朗數(shù)為0.56~0.85,滿足收斂性要求。同時(shí),這一時(shí)間步長也遠(yuǎn)遠(yuǎn)小于本文所涉及的壓力波動(dòng)周期,滿足數(shù)據(jù)采樣的要求。

    2.3 模型的驗(yàn)證

    2011年,針對高速列車裙板支架疲勞開裂的問題,文獻(xiàn)[5]通過長距離(1 200萬km)線路運(yùn)行測試得到了設(shè)備艙支架不同位置的受力情況。

    根據(jù)實(shí)測結(jié)果,當(dāng)列車以380 km/h的速度明線行駛時(shí),裙板兩側(cè)氣壓差為496 Pa,筋板與L形板連接部位焊縫角端受到較大的應(yīng)力,其最大值為26.1 MPa。

    本文分別以不同網(wǎng)格尺寸對計(jì)算域進(jìn)行網(wǎng)格劃分。不同網(wǎng)格計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果的對比見表1??梢钥闯觯褂镁W(wǎng)格2與網(wǎng)格3所得的仿真結(jié)果與實(shí)測結(jié)果之間的相對誤差均可控制在5%以內(nèi),故選擇計(jì)算成本較低的網(wǎng)格2進(jìn)行后續(xù)仿真分析。

    表1 仿真結(jié)果與實(shí)驗(yàn)結(jié)果對比

    3 流場特性分析

    車外流場及設(shè)備艙內(nèi)流場共同決定了裙板所受氣動(dòng)載荷的特性。以車速為300 km/h的明線、無交會(huì)、勻速直線行駛工況為例,分析兩部分流場以及氣動(dòng)載荷的特點(diǎn)。

    3.1 列車外流場分析

    圖6為用等壓面表示的列車行駛過程中車身周圍的流場結(jié)構(gòu)。在0.1 s時(shí)刻,位于尾車右側(cè)的等壓面A面積較大,而其左側(cè)的等壓面C面積較小且位置更靠近尾車前部。另外,中間車左側(cè)的等壓面B位于中間車的中部,該處正是設(shè)備艙通風(fēng)口1所在位置。在0.35 s時(shí)刻,等壓面A面積變小,等壓面C面積增大且更靠近尾車后部,而等壓面B則移至中間車的后部。0.6 s時(shí)刻等壓面的分布類似于0.1 s時(shí)刻,而0.85 s時(shí)刻又與0.35 s時(shí)類似。由此可見,在考慮設(shè)備艙內(nèi)流場流動(dòng)情況之后,仿真模型能夠體現(xiàn)列車周圍流場結(jié)構(gòu)的變化。這種變化具有周期性,是交變氣動(dòng)載荷產(chǎn)生的原因。

    圖6 外流場結(jié)構(gòu)變化(等壓面P=-50 Pa)

    如圖7所示,為了說明不同位置處壓力變化情況,在設(shè)備艙內(nèi)外分別設(shè)置9個(gè)觀察點(diǎn)。

    圖7 觀察點(diǎn)位置

    艙外觀察點(diǎn)1~6的壓力時(shí)程曲線如圖8所示。盡管列車以某一速度穩(wěn)定運(yùn)行,但由于列車周圍流場結(jié)構(gòu)進(jìn)行周期性變化,各觀察點(diǎn)的壓力隨時(shí)間均作正弦變化。其變化幅值與相位不同,而周期均相差很小。定義流場中同一位置壓力最大值與最小值的差為壓力變化幅值。以設(shè)備艙前端為原點(diǎn),艙外壓力變化幅值沿車長方向的變化如圖9所示??梢?,設(shè)備艙前端艙外壓力變化較小,而后端較大。由于受到轉(zhuǎn)向架處復(fù)雜擾動(dòng)的影響,在14 m位置處壓力變化有所降低。從總體看,設(shè)備艙后部的裙板受到氣動(dòng)載荷的變化較為劇烈。

    圖8 設(shè)備艙外觀察點(diǎn)壓力變化曲線

    圖9 壓力變化幅值在車長方向的變化(列車左側(cè))

    3.2 設(shè)備艙內(nèi)流場分析

    以圖10所示的流動(dòng)矢量圖來說明設(shè)備艙內(nèi)流場結(jié)構(gòu)的變化。設(shè)備艙前部的流體流動(dòng)緩慢,同時(shí)流動(dòng)變化較小。而其后部流動(dòng)情況復(fù)雜,流體在尾部空腔內(nèi)形成漩渦。漩渦與裙板之間流體的流動(dòng)方向作周期性變化:在0.25 s時(shí)刻,設(shè)備艙左后側(cè)流體產(chǎn)生一定的橫向流動(dòng);0.25 s之后,橫向流動(dòng)減少;而在0.75 s時(shí)刻,這種橫向流動(dòng)再次出現(xiàn),且方向與0.25 s時(shí)刻相同;隨后,橫向流動(dòng)再次減少。

    (a)t=0.25 s

    (b)t=0.5 s

    (c)t=0.75 s

    (d)t=1.0 s圖10 設(shè)備艙內(nèi)流動(dòng)速度矢量圖

    如圖11所示,設(shè)備艙內(nèi)觀察點(diǎn)7~9(圖7)的壓力同樣呈現(xiàn)出正弦變化。值得注意的是,在這種工況下,設(shè)備艙內(nèi)的流場結(jié)構(gòu)變化周期與車外流場變化周期基本相同。與艙外不同的是,艙內(nèi)各觀察點(diǎn)的壓力變化幅值小,且不同位置幅值相差較小。

    圖11 設(shè)備艙內(nèi)觀察點(diǎn)壓力變化曲線

    3.3 裙板氣動(dòng)載荷分析

    裙板兩側(cè)流場的壓差為裙板氣動(dòng)載荷的來源。圖12為設(shè)備艙四個(gè)拐角處裙板中心處壓力差的時(shí)程曲線??梢?,氣動(dòng)載荷同樣呈正弦變化。圖13顯示了各裙板中心處的最大壓力差。沿車長方向?qū)α熊囎笥覂蓚?cè)裙板進(jìn)行編號,編號越大越靠近設(shè)備艙后端??梢钥闯鰵鈩?dòng)載荷呈現(xiàn)出后部大于前部,右側(cè)大于左側(cè)的分布趨勢。造成這種現(xiàn)象的原因在于靠近車身后方的流場變化更加劇烈,而不對稱的通風(fēng)口布置造成列車左右兩側(cè)氣動(dòng)載荷具有一定差異。

    圖12 壓力差變化曲線

    圖13 裙板中心最大壓力差

    3.4 車速對流場特性的影響

    本文計(jì)算了不同列車行駛速度的多種工況。以受較劇烈氣動(dòng)載荷作用的圖3(c)中所示裙板為研究對象,考察其中心處壓力差與車速的關(guān)系。圖14為該裙板中心兩側(cè)的最大壓力差以及壓力差變化頻率隨列車行駛速度變化的曲線??梢姡畲髩毫Σ罴白兓l率均隨車速的提升而增加,其中最大壓力差增加明顯而變化頻率增幅很小。

    圖14 車速對氣動(dòng)載荷的影響

    4 裙板支架的動(dòng)載荷響應(yīng)

    在復(fù)雜的氣動(dòng)載荷作用下,列車裙板在設(shè)備艙內(nèi)外壓力差的作用下產(chǎn)生形變,這些形變由連接結(jié)構(gòu)傳遞至裙板支架,從而導(dǎo)致支架局部位置產(chǎn)生應(yīng)力集中。同時(shí),由于氣動(dòng)載荷的交變特性,應(yīng)力集中區(qū)域易產(chǎn)生疲勞破壞現(xiàn)象,影響其使用壽命。

    4.1 氣動(dòng)載荷作用下裙板的變形

    車速為300 km/h情況下,圖3(c)中所示裙板的位移情況如圖15所示。由于該裙板的固定支架位置偏左,故其右端位移較大,左端位移較小,產(chǎn)生扭曲,對支架承載造成影響。

    圖15 裙板位移云圖(外側(cè))

    以該裙板中心處為觀察點(diǎn),其壓力差及各向(其中x方向?yàn)檐囬L方向,y方向?yàn)檐噹叨确较颍瑉方向?yàn)檐噹麑挾确较?位移曲線如圖16所示。其中,壓力與各方向位移的變化規(guī)律基本相同,壓力變化的幅值與頻率直接決定裙板位移變化的幅值與頻率。而裙板在y、z方向的位移值明顯大于x方向,說明y、z方向的交變位移是造成裙板支架疲勞破壞的主要因素。

    圖16 裙板中心壓力差及位移變化曲線

    4.2 裙板支架的動(dòng)應(yīng)力特征

    裙板由于氣動(dòng)載荷作用產(chǎn)生的位移傳遞至裙板支架,在裙板支架的某些區(qū)域產(chǎn)生高應(yīng)力。計(jì)算結(jié)果表明,裙板支架的最大應(yīng)力值出現(xiàn)在筋板連接處,與實(shí)際開裂位置(圖1)相同。

    由圖17可知,裙板支架筋板連接處應(yīng)力在一定范圍內(nèi)呈周期變化。在這種交變應(yīng)力的作用下,筋板連接處容易產(chǎn)生疲勞裂紋,進(jìn)而產(chǎn)生斷裂現(xiàn)象。

    圖17 裙板支架筋板連接處應(yīng)力變化曲線

    4.3 車速對動(dòng)應(yīng)力的影響

    由圖18可知,隨著列車行駛速度的增加,裙板支架筋板連接處的最大應(yīng)力逐漸增大,這將導(dǎo)致結(jié)構(gòu)更容易產(chǎn)生疲勞破壞。

    圖18 車速對裙板支架應(yīng)力的影響

    4.4 裙板支架破壞機(jī)理的討論

    以列車車速為300 km/h的工況為例。為了考察氣動(dòng)載荷與機(jī)械載荷的耦合作用機(jī)理,分別對裙板及支架施加氣動(dòng)載荷、機(jī)械載荷以及同時(shí)施加氣動(dòng)載荷及機(jī)械載荷組成的耦合載荷。

    氣動(dòng)載荷的加載由上述流固耦合模型完成。對于機(jī)械載荷,本文假定在列車行駛過程中,裙板及支架受呈正弦變化的振幅為1.4g、頻率為2 Hz的垂向振動(dòng)加速度。

    圖19為三種載荷下,同一裙板支架筋板連接處的應(yīng)力變化時(shí)程曲線??梢?,三種載荷工況下,結(jié)構(gòu)應(yīng)力均呈周期性變化。與氣動(dòng)載荷相比,機(jī)械載荷作用下結(jié)構(gòu)應(yīng)力較小。在耦合載荷作用下,結(jié)構(gòu)的響應(yīng)更為復(fù)雜,與氣動(dòng)載荷與機(jī)械載荷下應(yīng)力的線性疊加結(jié)果有明顯差異,在某些時(shí)間段內(nèi)結(jié)構(gòu)應(yīng)力值有所降低。

    圖19 不同載荷作用下結(jié)構(gòu)應(yīng)力時(shí)域?qū)Ρ?/p>

    如圖20所示,通過頻域分析可知,氣動(dòng)載荷作用下,結(jié)構(gòu)應(yīng)力最大頻率為3.5 Hz,而機(jī)械載荷為4 Hz。兩者耦合則會(huì)造成應(yīng)力頻域特性的改變。耦合載荷作用下,應(yīng)力幅值的頻率主要集中于氣動(dòng)載荷起主導(dǎo)作用的3.5 Hz以及兩種載荷之間的3.75 Hz。

    圖20 不同載荷作用下結(jié)構(gòu)應(yīng)力頻域?qū)Ρ?/p>

    5 結(jié)論

    對明線行駛過程中列車裙板支架疲勞破壞問題進(jìn)行了流固耦合仿真分析。得到如下結(jié)論:

    (1)導(dǎo)致裙板支架產(chǎn)生疲勞破壞的主要原因?yàn)椋毫熊囆旭傔^程中車身周圍及設(shè)備艙內(nèi)流場結(jié)構(gòu)發(fā)生周期性變化,造成設(shè)備艙內(nèi)外壓力差發(fā)生周期變化,從而使裙板受到周期性氣動(dòng)載荷的作用。

    (2)氣動(dòng)載荷具有如下特征:靠近列車后部裙板所受氣動(dòng)載荷幅值較大,而所有裙板所受氣動(dòng)載荷頻率基本相同;隨著列車行駛速度的增加,氣動(dòng)載荷幅值明顯增加,而頻率增加幅度較小。

    (3)在氣動(dòng)載荷作用下,裙板在列車高度及寬度方向的交變位移是造成裙板支架疲勞破壞的主要原因。裙板支架筋板連接處出現(xiàn)應(yīng)力集中,該位置應(yīng)力小于其材料的屈服強(qiáng)度,但呈現(xiàn)周期性變化,易造成疲勞破壞。氣動(dòng)載荷與機(jī)械載荷的耦合作用比較復(fù)雜,耦合載荷作用下結(jié)構(gòu)響應(yīng)的時(shí)域特性與頻域特性與兩種載荷分別作用以及兩者線性組合作用下均不相同,有待進(jìn)一步深入研究。

    參考文獻(xiàn):

    [1]杜秋男,李瑞淳.高速動(dòng)車組裙板設(shè)計(jì)研究[J].鐵道車輛,2008,46(6):8-13.

    DU Qiunan,LI Ruichun.Design Research of Skirt Board of High-speed Train[J].Rolling Stock,2008,46(6):8-13.

    [2]張驍勇.材料的斷裂與控制[M].西安:西北工業(yè)大學(xué)出版社,2012:30-31.

    [3]熊峻江.疲勞斷裂可靠性工程學(xué)[M].北京:國防工業(yè)出版社,2008:4-5.

    [4]殷之平.結(jié)構(gòu)疲勞與斷裂[M].西安:西北工業(yè)大學(xué)出版社,2012:28-32.

    [5]佟勝江.設(shè)備艙支架裂紋原因分析及改進(jìn)措施[J].鐵道車輛,2011,49(12):10-14.

    TONG Shengjiang.Crack Cause Analysis of Equipment Cabin Supports and Improvement Measure[J].Rolling Stock,2011,49(12):10-14.

    [6]李振平.CRH3動(dòng)車組裙板支座載荷識別及疲勞壽命評估[D].北京:北京交通大學(xué),2010.

    [7]李人憲,趙晶,劉杰,等.高速列車會(huì)車壓力波對側(cè)窗的影響[J].機(jī)械工程學(xué)報(bào),2010,46(4):87-98.

    LI Renxian,ZHAO Jing,LIU Jie,et al.Influence of Air Pressure Pulse on Side Windows of High-speed Trains Passing Each Other[J].Journal of Mechanical Engineering,2010,46(4):87-98.

    [8]張軍海.基于FLUEN丁的高速機(jī)車氣動(dòng)刮雨器數(shù)值模擬研究[D].成都:西南交通大學(xué),2010.

    [9]郭迪龍,姚拴寶,劉晨輝,等.高速列車受電弓非定常氣動(dòng)特性研究[J].鐵道學(xué)報(bào),2012,34(11):16-21.

    GUO Dilong,YAO Shuanbao,LIU Chenhui,et al.Unsteady Aerodynamic Characteristics of High-speed Pantograph[J].Journal of the China Railway Society,2012,34(11):16-21.

    [10]龔曉波.動(dòng)車組自動(dòng)開閉裝置空氣動(dòng)力效應(yīng)與結(jié)構(gòu)強(qiáng)度分析[D].長沙:中南大學(xué),2010.

    [11]ZILIAN A,DINKLER D,VEHRE A.Projection-based Reduction of Fluid-structure Interaction Systems Using Monolithic Space-time Modes[J].Computer Methods in Applied Mechanics and Engineering,2009(198):3 795-3 805.

    [12]SOARES D.Fluid-structure Interaction Analysis by Optimised Boundary Element-finite Element Coupling Procedures[J].Journal of Sound and Vibration,2009(322):184-195.

    [13]BRUMMELEN E H V.Partitioned Iterative Solution Methods for Fluid-structure Interaction[J].International Journal for Numerical Methods in Fluids,2011(65):3-27.

    [14]OSTASEVICIUS V,DAUKSEVICIUSA R,GAIDYS R,et al.Numerical Analysis of Fluid-structure Interaction Effects on Vibrations of Cantilever Microstructure[J].Journal of Sound and Vibration,2007(308):660-673.

    [15]田紅旗.列車空氣動(dòng)力學(xué)[M].北京:中國鐵道出版社,2007.

    [16]錢若軍,董石麟,袁行飛.流固耦合理論研究進(jìn)展[J].空間結(jié)構(gòu),2008,14(1):3-15.

    QIAN Ruojun,DONG Shilin,YUAN Xingfei.Advances in Research on Fluid-structure Interaction Theory[J].Spatial Structures,2008,14(1):3-15.

    [17]黃志祥,陳立,蔣科林.高速列車模型編組長度和風(fēng)擋結(jié)構(gòu)對氣動(dòng)阻力的影響[J].實(shí)驗(yàn)流體力學(xué),2012,26(5):36-41.

    HUANG Zhixiang,CHEN Li,JIANG Kelin.Influence of Length of Train Formation and Vestibule Diaphragm Structure on Aerodynamic Drag of High Speed Train Model[J].Journal of Experiments in Fluid Mechanics,2012,26(5):36-41.

    猜你喜歡
    裙板氣動(dòng)流場
    中寰氣動(dòng)執(zhí)行機(jī)構(gòu)
    大型空冷汽輪發(fā)電機(jī)轉(zhuǎn)子三維流場計(jì)算
    基于NACA0030的波紋狀翼型氣動(dòng)特性探索
    基于有限元的動(dòng)車裙板強(qiáng)度分析
    基于反饋線性化的RLV氣動(dòng)控制一體化設(shè)計(jì)
    輕軌列車轉(zhuǎn)向架裙板設(shè)計(jì)與分析
    轉(zhuǎn)杯紡排雜區(qū)流場與排雜性能
    基于HYCOM的斯里蘭卡南部海域溫、鹽、流場統(tǒng)計(jì)分析
    基于瞬態(tài)流場計(jì)算的滑動(dòng)軸承靜平衡位置求解
    基于焊接疲勞壽命仿真的裙板優(yōu)化設(shè)計(jì)
    焊接(2015年1期)2015-07-18 11:07:32
    大香蕉97超碰在线| 床上黄色一级片| 有码 亚洲区| 久久人人爽人人爽人人片va| 我要看日韩黄色一级片| 亚洲va在线va天堂va国产| 久久久久九九精品影院| 久久精品国产亚洲网站| 国产一级毛片在线| 欧美高清性xxxxhd video| 色尼玛亚洲综合影院| 中国美白少妇内射xxxbb| 国产精品人妻久久久久久| 少妇熟女aⅴ在线视频| 少妇的逼水好多| www.av在线官网国产| 国产午夜福利久久久久久| 欧美区成人在线视频| 国产高潮美女av| 日本与韩国留学比较| 国产亚洲5aaaaa淫片| 日本免费在线观看一区| 看黄色毛片网站| 老女人水多毛片| 久久久久久久久中文| 亚洲无线观看免费| 岛国毛片在线播放| 亚洲色图av天堂| 蜜臀久久99精品久久宅男| 国产精品av视频在线免费观看| 国产欧美另类精品又又久久亚洲欧美| 欧美性猛交╳xxx乱大交人| 成人午夜高清在线视频| 精品欧美国产一区二区三| 成人av在线播放网站| 在线免费观看不下载黄p国产| 国产精品伦人一区二区| 国产一区有黄有色的免费视频 | 一级av片app| 狠狠精品人妻久久久久久综合| 少妇熟女欧美另类| 成年av动漫网址| 天堂网av新在线| 国产亚洲精品av在线| 亚洲欧美日韩东京热| 九九在线视频观看精品| 国产精品无大码| 最近最新中文字幕免费大全7| 青春草亚洲视频在线观看| 好男人视频免费观看在线| 亚洲欧美一区二区三区国产| 在线观看人妻少妇| 亚洲在久久综合| 99视频精品全部免费 在线| 寂寞人妻少妇视频99o| 亚洲乱码一区二区免费版| 亚洲av.av天堂| 国产精品爽爽va在线观看网站| 国产黄片视频在线免费观看| 午夜精品在线福利| 欧美日韩亚洲高清精品| 18+在线观看网站| 精品酒店卫生间| 麻豆久久精品国产亚洲av| 欧美xxⅹ黑人| 亚洲性久久影院| 成人二区视频| 日本猛色少妇xxxxx猛交久久| eeuss影院久久| 精品少妇黑人巨大在线播放| 内地一区二区视频在线| 色播亚洲综合网| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 午夜福利在线观看吧| 青青草视频在线视频观看| 97在线视频观看| 色哟哟·www| 中文精品一卡2卡3卡4更新| 噜噜噜噜噜久久久久久91| 能在线免费看毛片的网站| 三级国产精品片| 久久久久久久久久久丰满| 中文字幕av在线有码专区| 精品不卡国产一区二区三区| av国产免费在线观看| 久久99热这里只有精品18| 日本猛色少妇xxxxx猛交久久| 日韩av在线大香蕉| 国产单亲对白刺激| 联通29元200g的流量卡| 色网站视频免费| 深夜a级毛片| 国产国拍精品亚洲av在线观看| 国产高清有码在线观看视频| 亚洲一区高清亚洲精品| 久久久欧美国产精品| 80岁老熟妇乱子伦牲交| 精品一区二区三区人妻视频| 99久久精品热视频| 国精品久久久久久国模美| 色播亚洲综合网| 非洲黑人性xxxx精品又粗又长| av在线老鸭窝| 一个人看的www免费观看视频| 亚洲精品久久久久久婷婷小说| 国产 亚洲一区二区三区 | 亚洲18禁久久av| 亚洲av一区综合| 中文乱码字字幕精品一区二区三区 | 久久久久网色| 亚洲av不卡在线观看| 欧美潮喷喷水| 日本免费在线观看一区| 日本与韩国留学比较| 国产黄a三级三级三级人| 国产av国产精品国产| 毛片一级片免费看久久久久| 免费看av在线观看网站| 欧美潮喷喷水| 精品人妻视频免费看| 国产av国产精品国产| 亚洲国产精品国产精品| 22中文网久久字幕| 深爱激情五月婷婷| 中文字幕制服av| .国产精品久久| 在线免费观看的www视频| 一本一本综合久久| 国产一区二区三区综合在线观看 | 2021天堂中文幕一二区在线观| 亚洲国产欧美在线一区| 好男人视频免费观看在线| 国产精品久久久久久久久免| 精品国内亚洲2022精品成人| 欧美xxⅹ黑人| 久久久久久伊人网av| 国产精品三级大全| 国产成人一区二区在线| 欧美xxxx性猛交bbbb| 三级毛片av免费| 欧美3d第一页| 高清毛片免费看| 尤物成人国产欧美一区二区三区| 亚洲人成网站在线观看播放| 亚洲欧美日韩东京热| 国产精品不卡视频一区二区| 免费看光身美女| 成人毛片60女人毛片免费| 国产熟女欧美一区二区| 身体一侧抽搐| 2021天堂中文幕一二区在线观| 国产精品一及| 在线免费观看的www视频| 日本wwww免费看| 少妇被粗大猛烈的视频| 街头女战士在线观看网站| 国产永久视频网站| 一区二区三区高清视频在线| 肉色欧美久久久久久久蜜桃 | 免费看光身美女| 国产精品一区二区在线观看99 | 国产淫片久久久久久久久| 日韩强制内射视频| 深爱激情五月婷婷| 精品一区在线观看国产| 精品一区在线观看国产| 久久97久久精品| 精品一区在线观看国产| 熟女人妻精品中文字幕| 色尼玛亚洲综合影院| 亚洲av中文字字幕乱码综合| 国产成人午夜福利电影在线观看| 不卡视频在线观看欧美| 日本黄大片高清| 久久99热这里只有精品18| 免费观看无遮挡的男女| 免费观看无遮挡的男女| av国产久精品久网站免费入址| 亚洲国产色片| 国产亚洲精品av在线| 国产黄色免费在线视频| 亚洲av一区综合| 成年免费大片在线观看| 欧美+日韩+精品| 伦精品一区二区三区| 国产精品人妻久久久久久| 成人一区二区视频在线观看| 日韩一区二区视频免费看| 中国美白少妇内射xxxbb| h日本视频在线播放| 九九在线视频观看精品| 蜜臀久久99精品久久宅男| 26uuu在线亚洲综合色| 亚洲伊人久久精品综合| 床上黄色一级片| 日本av手机在线免费观看| 国产在视频线精品| 男女那种视频在线观看| 肉色欧美久久久久久久蜜桃 | 美女主播在线视频| 国产精品人妻久久久影院| 久久久久久久午夜电影| 国产黄色视频一区二区在线观看| 国产精品人妻久久久久久| 欧美激情久久久久久爽电影| 久久精品国产自在天天线| 女人十人毛片免费观看3o分钟| 秋霞在线观看毛片| 国产精品麻豆人妻色哟哟久久 | 亚洲av电影在线观看一区二区三区 | 久久久久精品久久久久真实原创| 久久久久久久久久黄片| 国产在视频线精品| 观看免费一级毛片| 最近中文字幕2019免费版| 精品一区二区免费观看| 午夜福利视频1000在线观看| 成人美女网站在线观看视频| 纵有疾风起免费观看全集完整版 | 3wmmmm亚洲av在线观看| 日日摸夜夜添夜夜爱| 精品99又大又爽又粗少妇毛片| 国产三级在线视频| 久久久久精品久久久久真实原创| 久久韩国三级中文字幕| 99热这里只有精品一区| 亚洲精品一区蜜桃| 日韩成人伦理影院| 国产一区二区三区av在线| 国产黄色小视频在线观看| 女人被狂操c到高潮| 欧美日韩综合久久久久久| 少妇人妻一区二区三区视频| 淫秽高清视频在线观看| 黄色欧美视频在线观看| 国产一级毛片七仙女欲春2| 国产精品精品国产色婷婷| 中文在线观看免费www的网站| 汤姆久久久久久久影院中文字幕 | 人人妻人人看人人澡| 日韩不卡一区二区三区视频在线| 人妻系列 视频| 精品国产一区二区三区久久久樱花 | 久久久久免费精品人妻一区二区| 国产男人的电影天堂91| 如何舔出高潮| 蜜臀久久99精品久久宅男| 午夜免费激情av| 99re6热这里在线精品视频| 亚洲18禁久久av| 欧美极品一区二区三区四区| 汤姆久久久久久久影院中文字幕 | 黄色一级大片看看| 搡老乐熟女国产| 国产在线一区二区三区精| 午夜福利在线在线| 能在线免费观看的黄片| 插逼视频在线观看| 嫩草影院新地址| 国产欧美日韩精品一区二区| 国产精品国产三级国产av玫瑰| 婷婷色av中文字幕| 亚洲真实伦在线观看| 亚洲综合精品二区| 国产成人精品福利久久| 一夜夜www| 午夜免费观看性视频| 亚洲欧美成人综合另类久久久| 国产美女午夜福利| 欧美日韩综合久久久久久| 久久综合国产亚洲精品| 中文字幕av在线有码专区| 99久久中文字幕三级久久日本| 干丝袜人妻中文字幕| 国产成人精品久久久久久| 久久99热这里只频精品6学生| 成年女人在线观看亚洲视频 | 男插女下体视频免费在线播放| 色5月婷婷丁香| 国产亚洲午夜精品一区二区久久 | 亚洲第一区二区三区不卡| 校园人妻丝袜中文字幕| 小蜜桃在线观看免费完整版高清| 国产成人精品婷婷| 成人美女网站在线观看视频| 我的女老师完整版在线观看| 亚洲不卡免费看| 亚洲精品久久久久久婷婷小说| 亚洲最大成人手机在线| 免费av不卡在线播放| 综合色丁香网| 人人妻人人看人人澡| 亚洲av免费在线观看| 中文字幕免费在线视频6| 日韩人妻高清精品专区| 国产老妇伦熟女老妇高清| .国产精品久久| 观看美女的网站| 国产成人午夜福利电影在线观看| 麻豆国产97在线/欧美| 国产大屁股一区二区在线视频| 日韩精品有码人妻一区| 18禁动态无遮挡网站| 国产乱人视频| 成人性生交大片免费视频hd| 免费观看在线日韩| 亚洲欧美清纯卡通| 欧美区成人在线视频| 国精品久久久久久国模美| 国产精品不卡视频一区二区| 乱人视频在线观看| kizo精华| 国产一区二区亚洲精品在线观看| 天堂影院成人在线观看| 国产精品久久久久久精品电影| 国产v大片淫在线免费观看| 蜜臀久久99精品久久宅男| 亚洲欧美中文字幕日韩二区| 午夜福利在线观看免费完整高清在| 欧美zozozo另类| 黄色欧美视频在线观看| 国产精品一区二区三区四区久久| 黄色欧美视频在线观看| 777米奇影视久久| 亚洲精品国产av蜜桃| 午夜福利在线在线| 久久精品夜色国产| 寂寞人妻少妇视频99o| 能在线免费看毛片的网站| 少妇高潮的动态图| 一二三四中文在线观看免费高清| 不卡视频在线观看欧美| 综合色丁香网| 搞女人的毛片| 国产综合精华液| 日韩成人伦理影院| a级毛片免费高清观看在线播放| 欧美一级a爱片免费观看看| 婷婷六月久久综合丁香| 又爽又黄a免费视频| 国内揄拍国产精品人妻在线| 精品99又大又爽又粗少妇毛片| 国产精品无大码| 国产乱来视频区| 午夜日本视频在线| av女优亚洲男人天堂| 国产午夜福利久久久久久| 最后的刺客免费高清国语| 看十八女毛片水多多多| 欧美xxxx黑人xx丫x性爽| 国产一区二区三区综合在线观看 | 舔av片在线| 成人欧美大片| 久久久久久久午夜电影| 看黄色毛片网站| 激情五月婷婷亚洲| 午夜视频国产福利| 亚洲av电影在线观看一区二区三区 | 日本熟妇午夜| 一夜夜www| 一个人看视频在线观看www免费| 99九九线精品视频在线观看视频| 人妻一区二区av| 免费看美女性在线毛片视频| 国产午夜精品论理片| 噜噜噜噜噜久久久久久91| 精品久久久久久久人妻蜜臀av| 亚洲内射少妇av| 久久久久久久亚洲中文字幕| 精品午夜福利在线看| 亚洲人成网站高清观看| 国产一区亚洲一区在线观看| 免费观看精品视频网站| 精品一区在线观看国产| 欧美日本视频| 综合色丁香网| 国产伦一二天堂av在线观看| 国产精品女同一区二区软件| 国产伦在线观看视频一区| 深爱激情五月婷婷| 黄色配什么色好看| 亚洲aⅴ乱码一区二区在线播放| 黄色一级大片看看| 久久99热这里只频精品6学生| 亚洲欧美日韩东京热| 成人漫画全彩无遮挡| 日韩成人伦理影院| 嫩草影院精品99| 黄色欧美视频在线观看| 精品一区在线观看国产| 人妻一区二区av| 搡老妇女老女人老熟妇| 国产精品.久久久| 国产成人a∨麻豆精品| 色综合色国产| 亚洲精品影视一区二区三区av| 在线a可以看的网站| 亚洲精品影视一区二区三区av| 亚洲av二区三区四区| 亚洲精品影视一区二区三区av| 国产精品伦人一区二区| 日日啪夜夜撸| 又爽又黄无遮挡网站| 美女大奶头视频| 欧美 日韩 精品 国产| 3wmmmm亚洲av在线观看| 日日啪夜夜爽| 亚洲自偷自拍三级| 国产一区有黄有色的免费视频 | 特大巨黑吊av在线直播| 亚洲怡红院男人天堂| 高清在线视频一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 成年女人在线观看亚洲视频 | 国产精品一二三区在线看| 亚洲精品,欧美精品| 成人二区视频| 亚洲国产av新网站| videos熟女内射| 在线免费观看的www视频| 丝袜喷水一区| 色哟哟·www| 日韩成人伦理影院| 少妇猛男粗大的猛烈进出视频 | 99热网站在线观看| 国产精品麻豆人妻色哟哟久久 | 不卡视频在线观看欧美| 嫩草影院入口| 精品熟女少妇av免费看| 亚洲国产精品成人久久小说| 日本与韩国留学比较| 国产精品久久久久久久久免| 久久久久久伊人网av| 又黄又爽又刺激的免费视频.| 九色成人免费人妻av| 午夜免费男女啪啪视频观看| 国产精品女同一区二区软件| 亚洲成色77777| 国产在视频线精品| 亚洲精品亚洲一区二区| 久热久热在线精品观看| av免费在线看不卡| 夜夜爽夜夜爽视频| or卡值多少钱| 免费大片黄手机在线观看| 六月丁香七月| 精品一区二区三区人妻视频| 天堂俺去俺来也www色官网 | 久久精品国产亚洲网站| 久久精品久久久久久噜噜老黄| 亚洲精品一区蜜桃| 精品久久国产蜜桃| 欧美激情国产日韩精品一区| 日韩欧美精品免费久久| 中文天堂在线官网| 色视频www国产| 美女大奶头视频| 国产乱人偷精品视频| av又黄又爽大尺度在线免费看| av专区在线播放| 永久免费av网站大全| 在线 av 中文字幕| 神马国产精品三级电影在线观看| 如何舔出高潮| 亚洲精品色激情综合| 天堂av国产一区二区熟女人妻| 久久久精品94久久精品| 日本熟妇午夜| 最近手机中文字幕大全| 亚洲激情五月婷婷啪啪| 国产单亲对白刺激| 国产 亚洲一区二区三区 | 免费少妇av软件| 久久久精品94久久精品| 岛国毛片在线播放| 久久国内精品自在自线图片| 久久久亚洲精品成人影院| 干丝袜人妻中文字幕| 搡女人真爽免费视频火全软件| 日本一本二区三区精品| 狂野欧美激情性xxxx在线观看| 国产精品一区二区三区四区免费观看| 禁无遮挡网站| 国产极品天堂在线| 中文字幕制服av| 国产成年人精品一区二区| 中文字幕av在线有码专区| 久久精品久久精品一区二区三区| 免费av不卡在线播放| 成人av在线播放网站| 街头女战士在线观看网站| 亚洲不卡免费看| 日韩av免费高清视频| 18+在线观看网站| 日本免费a在线| 色综合色国产| 在线 av 中文字幕| 免费观看a级毛片全部| 亚洲美女搞黄在线观看| 国产精品国产三级国产av玫瑰| 可以在线观看毛片的网站| 一夜夜www| 国产熟女欧美一区二区| 国产成人aa在线观看| 五月伊人婷婷丁香| 国产精品麻豆人妻色哟哟久久 | www.av在线官网国产| 亚洲精品国产av成人精品| 国内少妇人妻偷人精品xxx网站| 欧美丝袜亚洲另类| 日韩欧美精品v在线| 高清在线视频一区二区三区| 国产亚洲最大av| 日本熟妇午夜| 99久久精品国产国产毛片| 婷婷色麻豆天堂久久| 亚洲久久久久久中文字幕| 在线免费观看的www视频| 欧美变态另类bdsm刘玥| 欧美成人a在线观看| 麻豆久久精品国产亚洲av| 国产精品一二三区在线看| 草草在线视频免费看| av专区在线播放| 国产高清有码在线观看视频| 国产成人精品久久久久久| 又爽又黄无遮挡网站| 日韩大片免费观看网站| 亚洲性久久影院| 国模一区二区三区四区视频| 亚洲精品视频女| 色尼玛亚洲综合影院| 91久久精品国产一区二区三区| 少妇高潮的动态图| 最近2019中文字幕mv第一页| 久久精品久久久久久噜噜老黄| 啦啦啦韩国在线观看视频| 搡老妇女老女人老熟妇| 午夜福利网站1000一区二区三区| 国产精品久久久久久久久免| 777米奇影视久久| 久久鲁丝午夜福利片| 99久久精品国产国产毛片| 亚洲欧美日韩无卡精品| 成人综合一区亚洲| 日日干狠狠操夜夜爽| 国产高潮美女av| 中文字幕制服av| 亚洲欧洲日产国产| 欧美激情久久久久久爽电影| 水蜜桃什么品种好| 日韩不卡一区二区三区视频在线| av专区在线播放| 亚洲精品影视一区二区三区av| 精品人妻偷拍中文字幕| 亚洲国产精品国产精品| 久久久精品欧美日韩精品| 国产熟女欧美一区二区| 久久精品久久久久久久性| 五月玫瑰六月丁香| 日韩,欧美,国产一区二区三区| 亚洲国产色片| 日韩精品有码人妻一区| 18+在线观看网站| 能在线免费看毛片的网站| 街头女战士在线观看网站| 国产高清三级在线| 国产高清有码在线观看视频| 中文天堂在线官网| 国产在线一区二区三区精| 精品一区在线观看国产| 国产有黄有色有爽视频| 青春草国产在线视频| 亚洲自偷自拍三级| 日韩,欧美,国产一区二区三区| 午夜爱爱视频在线播放| 婷婷六月久久综合丁香| 中国美白少妇内射xxxbb| 亚洲欧美清纯卡通| 日韩中字成人| a级一级毛片免费在线观看| 日韩欧美精品免费久久| 亚洲欧美精品专区久久| 免费少妇av软件| 色播亚洲综合网| 日韩精品青青久久久久久| 日日啪夜夜爽| 久久综合国产亚洲精品| 国产免费福利视频在线观看| 精品99又大又爽又粗少妇毛片| 97在线视频观看| 大陆偷拍与自拍| 欧美丝袜亚洲另类| 亚洲自偷自拍三级| 欧美激情久久久久久爽电影| 在线免费观看不下载黄p国产| 国产成人免费观看mmmm| 特级一级黄色大片| 麻豆成人午夜福利视频| 搡女人真爽免费视频火全软件| 免费看不卡的av| 国产精品熟女久久久久浪| 26uuu在线亚洲综合色| 亚洲欧美成人精品一区二区| 一区二区三区乱码不卡18| 深夜a级毛片| 80岁老熟妇乱子伦牲交| 一级毛片我不卡| 日韩av在线免费看完整版不卡| 国产精品久久久久久久久免| 国内揄拍国产精品人妻在线| 免费看av在线观看网站| 亚洲美女视频黄频| 国产高清三级在线| 亚洲国产欧美在线一区|