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

    考慮熱力耦合的軌道車輪輻板參數(shù)優(yōu)化研究

    2016-11-12 06:47:04文永蓬周偉浩徐小峻尚慧琳李瓊
    鐵道科學與工程學報 2016年10期
    關鍵詞:輻板輪輻踏面

    文永蓬,周偉浩,徐小峻,尚慧琳, 李瓊

    (1. 上海工程技術大學 城市軌道交通學院,上海 201620;2. 上海應用技術大學 機械工程學院,上海 200235)

    ?

    考慮熱力耦合的軌道車輪輻板參數(shù)優(yōu)化研究

    文永蓬1,周偉浩1,徐小峻1,尚慧琳2, 李瓊1

    (1. 上海工程技術大學 城市軌道交通學院,上海 201620;2. 上海應用技術大學 機械工程學院,上海 200235)

    以提高雙S型軌道車輪的綜合性能為目的,建立雙S型軌道車輪模型,采用軌道車輪的熱力耦合計算方法,討論雙S型軌道車輪輻板在熱力耦合條件下的幾種優(yōu)化思路,通過正交實驗的方法對雙S型軌道車輪輻板的關鍵參數(shù)進行優(yōu)化。研究結果表明:當S彎方向為正以及結合圓弧處于正中位置時,雙S型輻板車輪具有較好安全性能;雙S型軌道車輪的輻板厚度以及靠近輪輞處的圓弧半徑對雙S型軌道車輪性能的影響較大。該研究不僅適用于雙S型軌道車輪的結構特性分析,同時還對其他輻板軌道車輪的優(yōu)化設計具有重要的參考價值。

    熱力耦合;軌道車輪;疲勞強度;正交實驗;參數(shù)優(yōu)化

    在城市軌道車輛運行過程中,由于站間距短、行車密度高等因素造成車輛頻繁的啟動和制動,車輪和軌道的摩擦、車輪與閘瓦的摩擦產生大量的熱,引起車輪產生較大的熱載荷,與此同時,車輪還承受來自車輛的全部載荷,熱載荷和結構載荷劇烈的共同作用往往造成車輪的疲勞損傷[1]。車輪的輻板是輪輞和輪轂的重要連接部件,在熱力耦合載荷作用下起到承上啟下的作用,對車輪的性能具有重要的影響。近年來,對車輪的疲勞分析及優(yōu)化設計進行了大量報道[2-13]。早期對車輪疲勞問題的研究大多只考慮結構載荷的作用[2-4],研究車輪在直道、彎道以及道岔處的應力大小,并算出安全系數(shù)與許用安全系數(shù)比較。目前,越來越多學者已經注意到溫度的改變對于車輪疲勞的影響,重視列車在踏面制動時伴隨著閘瓦與踏面間的摩擦車輪表面產生大量的熱,認為溫度的變化對車輪使用壽命的影響不可忽略。張萍[5]探討了熱流密度的加載方式,為車輪的熱分析及熱力耦合分析提供了重要的參考價值;Seo等[6-8]研究了車輪的熱疲勞問題,對車輪熱疲勞壽命預測方法進行了探討;候耐等[9-11]通過對車輪在結構載荷、熱載荷及2個載荷共同作用下的研究分析,認為結構載荷和熱載荷的共同作用是導致車輪疲勞失效的主要原因。對于車輪的優(yōu)化設計,徐傳來等[12]通過APDL語言對軌道車輪進行參數(shù)化建模,運用ANSYS一階優(yōu)化工具對車輪截面進行形狀化優(yōu)化,其優(yōu)化后車輪的質量少,但是應力有所增加;劉林芽等[13]通過遺傳算法對單S型輻板車輪的輻板關鍵參數(shù)進行優(yōu)化,以聲輻射功率為優(yōu)化目標,優(yōu)化后的車輪聲輻射功率較標準車輪聲輻射功率在大部分頻段上均有所減小??傊?,以上研究很少考慮到熱力耦合下軌道車輪輻板的參數(shù)優(yōu)化設計,也未能針對軌道車輪輻板多參數(shù)的特點提出多目標的優(yōu)化策略。為此,研究熱載荷和結構載荷的共同作用對軌道車輪的影響,以抗熱損性較好的雙S型車輪[14]作為研究對象,分析雙S型車輪輻板的結構特點,保持車輪的輪轂和輪輞不變,針對軌道車輪輻板提出多目標結構優(yōu)化思路,得到幾種優(yōu)化后的軌道車輪輻板,并進一步分析比較優(yōu)化后車輪的最大耦合應力、質量和踏面溫度等特性。

    1 軌道車輪有限元模型

    1.1雙S型車輪模型

    上海4號線使用的是馬鋼生產的國產車輪ZD840,為雙S型車輪。以該雙S型車輪為研究對象,利用ANSYS有限元分析軟件進行建模,在建模過程中忽略對計算結果較小的圓弧等細節(jié),采用8節(jié)點的六面體單元SOLID70建立有限元模型,該雙S型車輪的有限元模型中共含有21 490個節(jié)點和19 308個單元,劃分好網格的模型如圖1所示,該車輪材料為R9T,其基本屬性如表1所示。

    (a)三維截面圖;(b)二維截面圖圖1 雙S型車輪有限元模型Fig.1 FEM model for the double S rail wheel

    材料名稱數(shù)值泊松比ν0.3彈性模量E/MPa2×105密度ρ/(kg·m-3)7800強度極限σb/MPa900屈服極限σs/MPa580許用應力[σ]/MPa352疲勞極限σ-1/MPa315

    1.2載荷工況及加載條件

    我國城市軌道車輛車輪尚未建立計算載荷和計算載工況標準,一般是基于國際鐵路聯(lián)盟UIC510-5或歐盟的BS EN13979-1標準進行。在軌道車輛運行過程中,車輪受到輪軌之間的三向應力,包括橫向力、縱向力和垂向力,其中橫向力和縱向力的方向有2個,垂向力的方向為豎直向上,但是列車彎道或過道岔運行時,由于轉向架兩側重量分配不一致,使得一側的輪軌間垂向力大于另一側,導致兩側車輪受到的垂向力位置和大小均有區(qū)別。除了輪軌作用力,車輪還受到制動時閘瓦施加的閘瓦力,包括閘瓦壓力和閘瓦摩擦力,其中,閘瓦摩擦力根據(jù)車輪運行方向不同有2個方向?;谝延械膰?個標準,考慮車輪結構載荷(計算見文獻[15])和熱載荷共同作用,基本涵蓋了列車運行時的制動過程中車輪受到的各種載荷情況。

    2 軌道車輪的熱力耦合計算方法

    2.1溫度場分析

    溫度場分析首先要確定熱載荷條件,包括環(huán)境溫度、熱流密度和換熱系數(shù)。在計算熱載荷時用到的軌道車輪相關參數(shù)及制動條件如表2所示,其熱力學計算參數(shù)如表3所示。

    表2 軌道車輪參數(shù)及制動條件Table 2 Parameters and braking conditions of the rail wheel

    表3 軌道車輪熱力學參數(shù)Table 3 Thermodynamic parameters of the rail wheel

    2.1.1環(huán)境溫度

    根據(jù)城市軌道交通車輛車輪所處的工作條件,選取環(huán)境空氣溫度為24 ℃,車輪自身表面溫度為40 ℃。

    2.1.2熱流密度

    采用能量法進行計算,假設城市軌道交通車輛在制動過程中減少的動能全部轉化為熱能,根據(jù)動能定理以及熱流密度q的計算公式可得:

    (1)

    式中:m為列車質量;v0為初速度;t為制動時間;S為踏面與閘瓦的接觸面積;D為軌道車輪直徑;n為閘瓦數(shù);l為閘瓦寬度;η為熱流分配系數(shù)。

    2.1.3換熱系數(shù)

    對流換熱系數(shù)可由經驗公式[16],求得對流換熱系數(shù)h的計算公式:

    h=16.7(1+1.33v)

    (2)

    式中:v為車輛運行速度。

    由于溫度場分析屬瞬態(tài)熱分析,所以在選擇分析類型時選擇瞬態(tài)分析,將以上計算得到熱載荷值施加到雙S型車輪模型上,對車輪制動過程中的溫度場進行計算分析,由表2中制動減速度和制動初速度可計算出制動時間約為17 s。圖2為溫度最高時刻即10 s時的溫度云圖,圖3為制動過程中車輪踏面溫度隨時間變化關系圖,車輪在制動過程中踏面溫度先增大后減小,在10 s時刻車輪踏面溫度最大為192.2 ℃。

    單位:℃圖2 雙S型車輪溫度云圖Fig.2 Temperature contour for the double S rail wheel

    圖3 雙S型車輪17 s內踏面最高溫度變化情況Fig.3 Maximum tread temperature for the double S rail wheel within 17 s

    2.2熱力耦合分析

    采用間接耦合法,首先將熱分析模型中的熱單元SOLID70轉換為結構單元SOLID185,再將軌道車輪熱分析中得到的溫度場各節(jié)點溫度載荷作為體載荷施加到車輪上。并且,為了保持耦合前后的完整性,在熱應力分析時建立與溫度場一致的載荷步[17]。得到制動過程中車輪的熱力耦合情況,圖4為制動結束時刻熱力耦合應力云圖,由圖4可知,制動結束時,車輪最大應力位于輻板處,最大應力約為120.4 MPa。圖5為整個制動過程中最大耦合應力點耦合應力的變化情況,由圖5可知,最大耦合應力隨著制動時間逐步增大,在制動結束時刻達到最大值。

    單位: Pa圖4 雙S型車輪耦合應力云圖Fig.4 Thermal-mechanical coupling stress contour for the double S rail wheel

    圖5 雙S型車輪17 s內輻板最大耦合應力變化情況Fig.5 Maximum thermal-mechanical coupling stress for the double S rail wheel within 17 s

    綜上,雙S型車輪的踏面溫度最大值出現(xiàn)在制動過程中,而熱力耦合應力最大值出現(xiàn)在制動結束時刻。這2個性能參數(shù)描述了雙S型車輪重要的特性,踏面溫度與車輪踏面材料的屬性相關,熱力耦合應力與車輪的壽命相關。因此,在下面的研究中,把制動過程中踏面最高溫度和制動結束時刻最大耦合應力,作為熱力耦合下車輪的評價指標。此外,為節(jié)省成本,把車輪輕量化也作為一項評價指標。

    3 不同輻板參數(shù)下雙S型車輪性能

    圖6為雙S型輻板車輪結構特征圖。如圖6所示,雙S型車輪輻板存在2個S型彎和連接2個S彎的過渡圓弧三部分組成,S彎的厚度均為D。下面S彎形同“S”,稱為正S彎;上面S彎形同反“S”,稱為反S彎。保持雙S型車輪的陰影部分(輪輞和輪轂)固定不動,即a點、b點、c點和d點固定不動,將輻板的空間作為輻板設計域。

    在輻板設計域內,考慮雙S型車輪輻板2個S彎的正反方向、結合圓弧的位置、S彎半徑及輻板厚度,可分以下3種情況熱力耦合分析。

    圖6 雙S輻板車輪結構特征圖Fig.6 Structural features for the double S rail wheel

    3.1S彎方向

    在雙S型車輪中,規(guī)定雙S型車輪下面的S彎方向為雙S車輪的方向,則原車輪為正雙S車輪,可以建立反雙S車輪,正雙S車輪和反S車輪的二位截面如圖7所示,并對反雙S車輪進行熱力耦合分析,其結果如表4所示。

    圖7 正、反雙S車輪的二維模型Fig.7 Forward two dimensional model and the reverse one for the double S rail wheel

    Table 4 Comprehensive performance comparison of the forward model and the reverse model for the double S rail wheel

    正雙S車輪反雙S車輪耦合應力/MPa120.4121.4車輪質量/kg288288踏面溫度/℃192.2191.3

    表4中,正雙S車輪和反雙S車輪的質量一樣,踏面溫度和耦合應力均差不多,由于主要考慮耦合應力,正雙S車輪的耦合應力略小于反雙S車輪,因此,正雙S車輪性能優(yōu)于反雙S車輪,下面的分析中,均針對正雙S車輪進行。

    3.2結合圓弧位置

    在正雙S車輪輻板中,2個S彎呈軸對稱分布,2個S彎的大小和形狀都是相同的,所以2個S彎中間結合圓弧就在輻板的正中間位置,現(xiàn)考慮2個S彎的大小形狀不同,在重新建立模型的過程中,盡量保持車輪的質量以及結合圓弧的半徑不變,對原雙S車輪輻板進行修改,得到2個新的雙S車輪模型,并和原模型進行對比,3個模型的二維截面特征如圖8所示。

    圖8 結合圓弧位置不同的雙S型車輪的二維模型Fig.8 Two dimensional model of the double S rail wheel for the different arc position

    Table 5 Comprehensive performance comparison of three kinds of arc models for the double S rail wheel

    圓弧偏下圓弧正中圓弧偏上耦合應力/MPa134.3120.4125車輪質量/kg288.1288289.1踏面溫度/℃194192.2195

    對由于結合圓弧位置不同引起的新雙S輻板車輪進行熱力耦合分析,結果如表5所示。由表5可知,當結合圓弧在車輪輻板正中位置時,車輪質量最小,踏面溫度最低,輻板耦合應力也最低,因此,結合圓弧在正中位置時,雙S型車輪的性能較優(yōu)。

    3.3S彎半徑及輻板厚度

    由以上討論可知,S彎方向和結合弧位置在一定范圍內發(fā)生變化時,結合圓弧在輻板正中間的正雙S型車輪是相對較優(yōu)的。因此,本節(jié)主要對S彎方向及結合圓弧位置相對較優(yōu)的原車輪輻板的4個尺寸參數(shù)進行討論。

    根據(jù)圖6,雙S型車輪輻板由3段等半徑圓弧R1,R2和R3組成,且S彎的厚度D不變?,F(xiàn)以R1,R2,R3和D這4個車輪輻板關鍵參數(shù)作為優(yōu)化變量,以車輪輻板最大耦合應力、車輪質量和車輪踏面最高溫度作為優(yōu)化目標,現(xiàn)對每個變量參數(shù)賦予3個值,通過對這4個關鍵變量隨機組合來尋找優(yōu)化參數(shù),共81種組合,由于實驗次數(shù)較多,故本文采用可減少計算量的正交實驗設計方法進行分析。

    正交實驗設計方法[18]是實驗設計中的一種優(yōu)化方法,利用數(shù)理統(tǒng)計的原理對實驗結果進行科學地分析,通過少數(shù)實驗就能夠了解到各個因素對實驗指標的影響情況,并且能夠確定各個因素影響的主次順序,從而找出優(yōu)化參數(shù)組合。

    采用正交實驗的設計方法,分析輻板的4個參數(shù)對雙S型車輪性能的影響,找出性能較優(yōu)的參數(shù)組合,從而獲得優(yōu)化車輪。本實驗中有4個因素A,B,C,和D,分別對應輻板參數(shù)R1,R2,R3和D。每個因素有3個水平,即每個因素對應的3個值,列出所有因素水平表,如表6所示。

    利用正交實驗的設計方法,選用“四因素三水平”正交表L9(34),每個實驗條件做3次,每次實驗對3個指標測定,并取平均值,實驗方案及其實驗結果如表7所示,分析比較計算得到的實驗數(shù)據(jù),進而選取出較優(yōu)的雙S型車輪輻板的參數(shù)組合。

    表6 因素水平表Table 6 Factors and levels table

    表7 實驗方案及結果Table 7 Experimental scheme and results

    根據(jù)實驗分析結果,可繪制出各個優(yōu)化目標的效應曲線圖,如圖9所示。

    由圖9可以得出各優(yōu)化目標下以及綜合條件下雙S型車輪輻板參數(shù)的優(yōu)化組合,分析過程如下:

    1)耦合應力為主

    由圖9(a)可知,在單純以耦合應力為優(yōu)化目標下,各因素的主次順序為D>B>A>C,根據(jù)每個因素的每個水平的平均值大小可以得出耦合應力優(yōu)化組合為A2B1C2D3(優(yōu)化模型1)。

    2)質量為主

    由圖9(b)可知,在單純以質量為優(yōu)化目標下,各因素的主次順序為D>B>C>A,根據(jù)每個因素的每個水平的平均值大小可以得出耦合應力優(yōu)化組合為A3B3C3D1(優(yōu)化模型2)。

    3)踏面溫度為主

    由圖9(c)可知,在單純以踏面溫度為優(yōu)化目標下,各因素的主次順序為D>B>C>A,根據(jù)每個因素的每個水平的平均值大小可以得出耦合應力優(yōu)化組合為A2B3C1D3(優(yōu)化模型3)。

    4)綜合平衡考慮

    在以上3個優(yōu)化目標中,耦合應力關系到雙S型車輪的行車安全及使用壽命,為主要優(yōu)化目標,其次考慮成本等目標,兼顧車輛輕量化和車輛材料的加工工藝的要求。因此,優(yōu)化目標主次順序為:耦合應力>質量>踏面溫度,根據(jù)優(yōu)化目標的主次順序,可確定因素D和因素B對車輪性能影響較大,因素C和因素A相對較小。

    對于因素D:因素D對3個目標的影響均較大,但主要考慮的指標是耦合應力最小,故選D3,此時,踏面溫度D3也是較優(yōu)的。

    對于因素B:因素B對3個目標的影響也較大,其影響程度僅次于D,從耦合應力角度選B1。

    對于因素C:因素C對3個目標均是次要因素,耦合應力C1和C2差不多,從踏面溫度角度選C1。

    (a)耦合應力效應曲線;(b)質量效應曲線;(c)踏面溫度效應曲線圖9 各優(yōu)化目標的效應曲線Fig.9 Effect curve of each optimization objective

    對于因素A:因素A對以上3個指標也是次要因素,從耦合應力、踏面溫度角度選A2。

    故綜合考慮耦合應力、車輪質量以及踏面溫度后,確定綜合平衡選取的優(yōu)化組合為:A2B1C1D3(優(yōu)化模型4)。

    圖10分別建立上述以單目標優(yōu)化時優(yōu)化的3個模型以及綜合目標優(yōu)化的1個模型。由圖10可知,4個模型的S彎半徑和輻板厚度在原有的雙S型車輪輻板的基礎上都有一定程度的改變,優(yōu)化后的雙S型車輪耦合應力、車輪質量和踏面溫度均發(fā)生變化,結果如表8所示。由表8可知,以上4種優(yōu)化模型的耦合應力均有所下降,優(yōu)化模型1,3和4的耦合應力下降的最明顯,其中優(yōu)化模型1的耦合應力比原車輪下降最多,達16.7%,但是這3個優(yōu)化模型的質量卻都有上升,但是最大不超過2.5%;優(yōu)化模型2耦合應力較原車輪相近,但其質量比原車輪下降了2.8%;優(yōu)化模型3的踏面溫度較優(yōu),但只比原車輪低了0.9%。因此,由上述優(yōu)化效果可知,不改變車輪輻板的基本形狀,僅僅使S彎半徑和輻板厚度微小的改變,若造成車輪質量增加,則可降低其耦合應力;對于踏面溫度,僅改變輻板參數(shù),其優(yōu)化效果不明顯,但可通過使用優(yōu)化材料屬性來進行提升。

    綜上,以耦合應力、車輪質量以及踏面溫度為優(yōu)化目標,通過改變雙S型車輪輻板2個S彎的正反方向、結合圓弧的位置、S彎半徑及輻板厚度,確立了優(yōu)化目標及設計參數(shù)的優(yōu)先級,雙S型車輪輻板厚度D和靠近輪輞處的圓弧半徑R2這2個參數(shù)對車輪性能最靈敏,對雙S型車輪綜合性能影響顯著。因此,今后在加工制造過程中這2個參數(shù)的加工誤差需要格外關注。

    (a)優(yōu)化模型1;(b)優(yōu)化模型2;(c)優(yōu)化模型3;(d)優(yōu)化模型4圖10 優(yōu)化前后輻板形狀比較圖Fig.10 Comparison of results before and after optimization

    優(yōu)化車輪A2B1C2D3(模型1)A3B3C3D1(模型2)A2B3C1D3(模型3)A2B1C1D3(模型4)原車輪A2B2C2D2設R1/mm34.538.534.534.534.5計R2/mm30.538.538.530.534.5參R3/mm34.538.530.530.534.5數(shù)D/mm2521252523耦合應力/MPa100.3119.1104.7105.4120.4車輪質量/kg295.3279.9293292.9288踏面溫度/℃192.3196190.4192.9192.2優(yōu)化效果耦合應力降低16.7%質量降低2.8%溫度降低0.9%耦合應力降低12.5%

    4 結論

    1)當 S彎方向為正以及結合圓弧處于正中位置時,雙S型輻板車輪具有較好安全性能,且雙S型車輪輻板厚度D和靠近輪輞處的圓弧半徑R2對雙S型車輪性能影響較顯著,在加工制造時,對車輪這兩個參數(shù)需要格外重視。

    2)根據(jù)雙S型車輪多參數(shù)的特點,確立了以耦合應力為主、車輪質量和踏面溫度為輔的優(yōu)化思路,按照這個思路對輻板的4個尺寸參數(shù)進行平衡選優(yōu),從而得出綜合性能較優(yōu)的車輪,對多目標優(yōu)化問題提供參考價值。

    3)采用正交試驗的方法對雙S型軌道車輪進行參數(shù)優(yōu)化,減少了實驗次數(shù),能快速直觀地獲得雙S型車輪輻板的優(yōu)化參數(shù)組合,為軌道車輪的優(yōu)化設計提供了一條有效思路。

    4)優(yōu)化后的車輪較原車輪均有明顯的優(yōu)化效果,優(yōu)化模型一的耦合應力比原車輪下降達16.7%,因此,考慮熱力耦合效應,建議采用優(yōu)化模型一為城市軌道雙S型車輪的進一步研究的車輪模型。

    [1] 文永蓬,尚慧琳,董其煒,等.城市軌道車輛車輪輪緣磨耗分析[J].科技導報,2013,31(26):40-43.

    WEN Yongpeng, SHANG Huilin, DONG Qiwei, et al. Analysis on the wear of wheel flange for urban rail vehicle[J].Science & Technology Review,2013, 31(26):40-43.

    [2] 何瑩,劉志明,胡寶義.動車組車輪強度標準與分析方法[J].北京交通大學學報,2009,33(1) :15-19.

    HE Ying,LIU Zhiming,HU Baoyi. EMU’s wheel strength standards and analytical methods [J].Journal of Beijing Jiaotong University,2009,33(1) :15-19.

    [3] 丁輝.鐵道車輛車輪強度分析和評定方法[J].鐵道機車車輛,2003,23(4):12-14.

    DING Hui.Intensity analysis evaluation method of the wheel of the railway vehicle[J].Railway Locomotive & Car,2003,23(4):12-14.

    [4] 王宏林,李芾,黃運華.25 t軸重重載貨車車輪強度分析[J].鐵道車輛,2007,45(12):1-4.

    WANG Honglin,LI Hu,HUANG Yunhua.Strength analysis of wheels for heavy haul freight cars with the axle-load of 25 t[J].Rolling Stock,2007,45(12):1-4.

    [5] 張萍.地鐵車輪踏面制動熱負荷研究[D].成都:西南交通大學,2010.

    ZHANG Ping.A study of thermal load of tred braking for metro wheel[D].Chengdu: Southwest Jiaotong University,2010.

    [6] Seo J W,Kwon S J,Jun H K,et al.Effects of esidual stress and shape of web plate on the fatigue life of railway wheels[J].Engineering Failure Analysis,2009,16:2493-2507.

    [7] Teimourimanesh S,Vernersson T,Lunden R,et al.Tread braking of railway wheels-temperatures generated by a metro train[J].Original Article,2014,228(2):210-221.

    [8] 盧立麗.貨車制動踏面熱損傷研究[D].成都:西南交通大學,2007.

    LU Lili.Research on tread brake heat injury of freight wheel[D].Chengdu: Southwest Jiaotong University,2007.

    [9] 候耐.重載貨車車輪踏面制動熱負荷研究[D].成都:西南交通大學,2011.

    HOU Nai.Heat load analysis of heavy-haul frecight car wheel for tread breaking[D].Chengdu: Southwest Jiaotong University,2011.

    [10] 肖楠,謝基龍,周素霞.地鐵車輪踏面制動疲勞強度評價方法及應用[J].工程力學,2010,27(9) :234-239.

    XIAO Nan,XIE Jilong,ZHOU Suxia.Method to assess fatigue strength of tread breaking of mehicle wheel and its application[J].Engineering Mechanics,2010, 27(9):234-239.

    [11] 徐傳來,米彩盈,李芾.基于軸對稱模型的貨車車輪結構疲勞強度分析[J].交通運輸工程學報,2008,8(6):20-23.

    XU Chuanlai,MI Caiying,LI Fu.Structural fatigue strength analysis of wheel for freight car base on axisymmetric model[J].Journal of Traffic and Transportation Engineering,2008,8(6):20-23.

    [12] 徐傳來,米彩盈,周仲榮.基于APDL語言的車輪參數(shù)化形狀優(yōu)化[J].鐵道學報,2011,33(11):23-27.

    XU Chuanlai,MI Caiying,ZHOU Zhongrong.Wheel shape & parameter optimization based on APDL language[J].Journal of the China Railway Society,2011,33(11):23-27.

    [13] 劉林芽,張斌,邵文杰,等.S型輻板車輪聲輻射優(yōu)化設計[J].交通運輸工程學報,2013,13(5):55-60.

    LIU Linya,ZHANG Bin,SHAO Wenjie,et al.Optimization design of acoustic radiation for S-form web plate wheel[J].Journal of Traffic and Transportation Engineering, 2013,13(5):55-60.

    [14] 王慶峰.關于統(tǒng)一城軌車輛車輪型式的研究[J].現(xiàn)代城市軌道交通,2009(2):55-57.

    WANG Qingfeng.Study on unification of wheel type of urban rail vehicles [J].Modern Urban Transit,2009(2):55-57.

    [15] 田劍鋒,趙永翔. Ф840mm地鐵車輪疲勞強度分析[J].機械設計與研究,2013,40(10):45-49.

    TIAN Jianfeng,ZHAO Yongxiang.Strength analysis of Ф840mm subway wheel[J].Journal of Machine Design and Research,2013,40(10):45-49.

    [16] 陳倩,李芾,王軍平,等.40t軸重貨車制動熱負荷分析[J].鐵道機車車輛,2012,32(4):32-36.

    CHEN Qian,LI Fu,WANG Junping,et al.ThermaI load anaIysis of wheels for 40t axle-load freight car under braking conditions[J].Railway Locomotive & Car,2012, 32(4):32-36.

    [17] 李金良,肖楠,謝基龍.重載貨車車輪踏面制動輻板熱應力分析[J].機械工程學報,2012,48(12):133-138.

    LI Jinliang,XIAO Nan,XIE Jilong.Thermal stress analysis of the heavy-haul freight car’s wheel web plate under the wheel tread braking[J].Journal of Mechanical Engineering,2012,48(12):133-138.

    [18] 劉文卿.實驗設計[M].北京:清華大學出版社,2005.

    LIU Wenqing.Experiment design[M].Beijing: Tsinghua University Press, 2005.

    Study on parameter optimization for the rail wheel considering thermal-mechanical coupling

    WEN Yongpeng1,ZHOU Weihao1,XU Xiaojun1,SHANG Huilin2,LI Qiong1

    (1. College of Urban Railway Transportation, Shanghai University of Engineering Science, Shanghai 201620, China;2. School of Mechanical Engineering, Shanghai Institute of Technology, Shanghai 200235, China)

    In order to improve the comprehensive performance of the double S urban rail wheel, a double S rail wheel model was established. The method of thermal-mechanical coupling was introduced in this model. Considering the thermal-mechanical coupling, several thoughts of double S rail wheel optimization model were analyzed. The key parameters of double S rail wheel plate were optimized by the method of orthogonal experiment. The results show that the double S rail wheel has a good safety performance when S bending in the positive direction and circular arc in the center position. The plate thickness and the radius of circular arc near the rim have a greater influence on the performance of the double S rail wheel. The work is not only suitable for the structure analysis of double S rail wheel, but also has important reference value for the other rail wheel optimization design.

    thermal-mechanical coupling; rail wheel; fatigue strength; orthogonal experiment; parameter optimization

    2016-05-20

    國家自然科學基金資助項目(11472176);上海市自然科學基金資助項目(15ZR1419200);上海市研究生教育創(chuàng)新計劃學位點引導布局與建設培育項目(13sc002)

    文永蓬(1979-),男,江西永新人,博士,副教授,從事軌道車輛關鍵結構優(yōu)化設計,軌道車輛振動與控制無軌道車輛故障診斷技術的研究;E-mail:yp_wen@163.com

    U270.2

    A

    1672-7029(2016)10-2042-09

    猜你喜歡
    輻板輪輻踏面
    帶式輸送機傳動滾筒輻板輪轂強度分析與優(yōu)化設計
    踏面清掃器動作邏輯與成本控制
    科學家(2021年24期)2021-04-25 16:55:45
    輻板形狀和淬火方式對車輪輻板沉降的影響研究
    一種仿生非充氣輪胎
    輪輻局部厚度對鋁車輪力學性能的影響
    北京汽車(2019年6期)2019-02-12 05:20:42
    基于機器視覺的礦車踏面磨耗檢測設計
    測控技術(2018年8期)2018-11-25 07:42:24
    無氣輪胎及汽車
    輻板厚度變化對整體葉盤結構力學性能的影響
    橫向力對列車車輪踏面表層材料塑性變形的影響
    鐵道車輛車輪踏面反向優(yōu)化設計方法
    鐵道學報(2015年9期)2015-05-10 03:04:41
    亚洲精品国产精品久久久不卡| 青青草视频在线视频观看| 超碰成人久久| 人人妻,人人澡人人爽秒播| 亚洲,欧美精品.| 欧美精品一区二区大全| 国产精品久久电影中文字幕 | 亚洲自偷自拍图片 自拍| 欧美日韩亚洲高清精品| 国产一区二区激情短视频| 亚洲黑人精品在线| 精品欧美一区二区三区在线| 亚洲一区中文字幕在线| 国产av一区二区精品久久| 亚洲国产精品一区二区三区在线| 亚洲伊人久久精品综合| 亚洲精品av麻豆狂野| 黄网站色视频无遮挡免费观看| 久久亚洲精品不卡| 亚洲av国产av综合av卡| www.熟女人妻精品国产| 中文亚洲av片在线观看爽 | tube8黄色片| 美国免费a级毛片| 久久精品亚洲精品国产色婷小说| 91精品国产国语对白视频| 亚洲中文日韩欧美视频| 王馨瑶露胸无遮挡在线观看| 高清欧美精品videossex| 又黄又粗又硬又大视频| 超碰成人久久| 欧美精品啪啪一区二区三区| 成年版毛片免费区| 亚洲精品中文字幕在线视频| 男人舔女人的私密视频| 国产又色又爽无遮挡免费看| 欧美日韩中文字幕国产精品一区二区三区 | 午夜精品久久久久久毛片777| 国产一区二区 视频在线| 侵犯人妻中文字幕一二三四区| 午夜福利欧美成人| 久久久久国内视频| 欧美久久黑人一区二区| 老汉色∧v一级毛片| 深夜精品福利| 一区二区三区精品91| 熟女少妇亚洲综合色aaa.| 99国产精品一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品久久久久久人妻精品电影 | 久久久久久免费高清国产稀缺| 另类亚洲欧美激情| 成人黄色视频免费在线看| 色视频在线一区二区三区| 黑人巨大精品欧美一区二区mp4| av欧美777| 69av精品久久久久久 | 99re6热这里在线精品视频| 大香蕉久久成人网| 99精品在免费线老司机午夜| 新久久久久国产一级毛片| 久久中文字幕人妻熟女| av网站在线播放免费| 欧美激情久久久久久爽电影 | 久久人妻av系列| 久热爱精品视频在线9| 国产xxxxx性猛交| 欧美激情 高清一区二区三区| 国产高清激情床上av| 精品少妇久久久久久888优播| 国产精品免费大片| 国产日韩一区二区三区精品不卡| 国产精品一区二区在线不卡| 国产成人av激情在线播放| 国产高清videossex| 精品一区二区三区av网在线观看 | 午夜日韩欧美国产| 色综合婷婷激情| 老司机亚洲免费影院| 九色亚洲精品在线播放| 在线十欧美十亚洲十日本专区| 精品福利观看| 久久久精品区二区三区| 三上悠亚av全集在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲国产看品久久| 日本黄色日本黄色录像| 最黄视频免费看| 成人精品一区二区免费| 狠狠精品人妻久久久久久综合| 露出奶头的视频| 国产成人一区二区三区免费视频网站| 少妇猛男粗大的猛烈进出视频| av免费在线观看网站| 国产精品成人在线| 久久久久久久精品吃奶| 极品教师在线免费播放| 男女免费视频国产| 欧美日韩成人在线一区二区| 十八禁人妻一区二区| av在线播放免费不卡| 国产精品av久久久久免费| 在线av久久热| 欧美精品亚洲一区二区| 男女床上黄色一级片免费看| 操美女的视频在线观看| 国产精品国产高清国产av | 亚洲,欧美精品.| 久久久久久久精品吃奶| 久久精品亚洲av国产电影网| 麻豆av在线久日| 精品熟女少妇八av免费久了| 最新在线观看一区二区三区| 两人在一起打扑克的视频| 国产成人啪精品午夜网站| 脱女人内裤的视频| 免费日韩欧美在线观看| 欧美日韩黄片免| 97在线人人人人妻| 天天影视国产精品| 欧美老熟妇乱子伦牲交| 国产欧美亚洲国产| 性高湖久久久久久久久免费观看| 丝袜美腿诱惑在线| 国产av又大| 国产成人免费无遮挡视频| 日韩欧美国产一区二区入口| 夜夜爽天天搞| 国产成人免费无遮挡视频| 精品亚洲成国产av| 丝袜人妻中文字幕| 成人精品一区二区免费| 久久久久久久大尺度免费视频| 亚洲专区中文字幕在线| 18禁观看日本| 丁香欧美五月| 日韩欧美一区视频在线观看| 一本一本久久a久久精品综合妖精| 色婷婷久久久亚洲欧美| 亚洲专区字幕在线| 变态另类成人亚洲欧美熟女 | 亚洲精品国产区一区二| tube8黄色片| 国产精品免费一区二区三区在线 | 成人国语在线视频| 黄色成人免费大全| 色94色欧美一区二区| 亚洲欧美激情在线| 女人久久www免费人成看片| 国产区一区二久久| 精品一区二区三区视频在线观看免费 | 精品少妇内射三级| 他把我摸到了高潮在线观看 | 亚洲精品在线观看二区| 精品福利永久在线观看| 午夜日韩欧美国产| 国产精品免费一区二区三区在线 | 欧美中文综合在线视频| 丝瓜视频免费看黄片| 久久99一区二区三区| 搡老乐熟女国产| 久久久久久人人人人人| 亚洲七黄色美女视频| 最近最新中文字幕大全电影3 | 丝袜美足系列| 黑人巨大精品欧美一区二区mp4| 人妻 亚洲 视频| 日韩欧美免费精品| 国产97色在线日韩免费| 免费在线观看视频国产中文字幕亚洲| 亚洲国产看品久久| 桃红色精品国产亚洲av| 国产成人av激情在线播放| 侵犯人妻中文字幕一二三四区| 亚洲视频免费观看视频| 久久久久国内视频| 淫妇啪啪啪对白视频| 悠悠久久av| 国产一卡二卡三卡精品| 国产在线免费精品| 日本av免费视频播放| 亚洲精品久久午夜乱码| 国产精品久久久久久精品电影小说| 搡老岳熟女国产| 久久中文字幕人妻熟女| 精品久久久久久电影网| 精品少妇内射三级| 亚洲av成人不卡在线观看播放网| 十八禁网站网址无遮挡| 大香蕉久久成人网| 黄色怎么调成土黄色| 亚洲国产成人一精品久久久| 一个人免费在线观看的高清视频| 国产麻豆69| av电影中文网址| 日韩 欧美 亚洲 中文字幕| 精品国产一区二区三区久久久樱花| 麻豆国产av国片精品| 久久国产亚洲av麻豆专区| 欧美+亚洲+日韩+国产| 五月开心婷婷网| 中文字幕人妻熟女乱码| 91老司机精品| 黄色 视频免费看| 91精品三级在线观看| 亚洲国产av新网站| 一本久久精品| 国产又爽黄色视频| 欧美亚洲 丝袜 人妻 在线| 久久青草综合色| 国产高清视频在线播放一区| 欧美日韩中文字幕国产精品一区二区三区 | 一本—道久久a久久精品蜜桃钙片| 热99国产精品久久久久久7| 日韩欧美一区视频在线观看| 国产福利在线免费观看视频| 免费不卡黄色视频| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美国产精品va在线观看不卡| 每晚都被弄得嗷嗷叫到高潮| 国产不卡一卡二| 中国美女看黄片| 法律面前人人平等表现在哪些方面| 国产精品国产高清国产av | 免费看a级黄色片| 午夜福利在线免费观看网站| 女警被强在线播放| 精品一区二区三区av网在线观看 | 一二三四社区在线视频社区8| 国产av又大| 啦啦啦在线免费观看视频4| 黄色 视频免费看| 99在线人妻在线中文字幕 | 国产熟女午夜一区二区三区| 成人影院久久| 欧美激情久久久久久爽电影 | 90打野战视频偷拍视频| 性少妇av在线| 日韩免费高清中文字幕av| 亚洲伊人色综图| 亚洲中文av在线| 两人在一起打扑克的视频| 脱女人内裤的视频| 精品人妻1区二区| 国产亚洲精品久久久久5区| bbb黄色大片| 精品高清国产在线一区| 国产精品免费视频内射| 日本黄色日本黄色录像| 在线观看免费午夜福利视频| 高潮久久久久久久久久久不卡| 欧美日韩亚洲高清精品| 叶爱在线成人免费视频播放| 亚洲精品国产一区二区精华液| 麻豆乱淫一区二区| 日韩欧美三级三区| 精品国产乱码久久久久久小说| av电影中文网址| 国产精品久久久av美女十八| 国产福利在线免费观看视频| 正在播放国产对白刺激| 日韩免费av在线播放| 91成年电影在线观看| 高清av免费在线| 妹子高潮喷水视频| 99re6热这里在线精品视频| 精品国产一区二区久久| 夜夜骑夜夜射夜夜干| 精品一品国产午夜福利视频| 桃红色精品国产亚洲av| 99精品在免费线老司机午夜| 999精品在线视频| 老司机影院毛片| 高清av免费在线| 中文字幕高清在线视频| 每晚都被弄得嗷嗷叫到高潮| 久久人妻熟女aⅴ| 大片电影免费在线观看免费| 中文字幕精品免费在线观看视频| 天天影视国产精品| 久久国产精品大桥未久av| 99久久国产精品久久久| 在线观看免费高清a一片| 无遮挡黄片免费观看| 国产精品 国内视频| 天堂8中文在线网| 国产欧美日韩一区二区三区在线| 日韩熟女老妇一区二区性免费视频| 1024视频免费在线观看| 成人18禁在线播放| 亚洲精品一卡2卡三卡4卡5卡| 欧美av亚洲av综合av国产av| 国产精品av久久久久免费| 自拍欧美九色日韩亚洲蝌蚪91| 97人妻天天添夜夜摸| 久热这里只有精品99| 久久九九热精品免费| 精品少妇内射三级| 在线观看人妻少妇| cao死你这个sao货| 国产精品免费大片| 久热爱精品视频在线9| 大香蕉久久网| 一本色道久久久久久精品综合| 国产精品99久久99久久久不卡| 久久久欧美国产精品| bbb黄色大片| 水蜜桃什么品种好| 岛国毛片在线播放| 人妻久久中文字幕网| 手机成人av网站| 亚洲视频免费观看视频| 黄片大片在线免费观看| 天堂8中文在线网| 少妇 在线观看| 十八禁网站免费在线| 亚洲伊人久久精品综合| 欧美性长视频在线观看| 亚洲av美国av| 精品欧美一区二区三区在线| 亚洲熟女毛片儿| 人人妻人人澡人人爽人人夜夜| 国产97色在线日韩免费| 又紧又爽又黄一区二区| 国产午夜精品久久久久久| √禁漫天堂资源中文www| 欧美+亚洲+日韩+国产| 中亚洲国语对白在线视频| 高清在线国产一区| 精品国产超薄肉色丝袜足j| 捣出白浆h1v1| 中文字幕色久视频| 国产精品免费视频内射| 国产成人精品在线电影| 精品午夜福利视频在线观看一区 | 香蕉久久夜色| 午夜福利视频在线观看免费| 大型av网站在线播放| 777米奇影视久久| 麻豆乱淫一区二区| 国产一区二区三区综合在线观看| 欧美成狂野欧美在线观看| 搡老熟女国产l中国老女人| 考比视频在线观看| 视频区图区小说| 欧美激情久久久久久爽电影 | 国产一区二区三区综合在线观看| 人成视频在线观看免费观看| 人人妻人人添人人爽欧美一区卜| 黄色 视频免费看| 男女下面插进去视频免费观看| 亚洲精品久久午夜乱码| 欧美精品高潮呻吟av久久| 久久亚洲真实| 亚洲色图综合在线观看| 国产欧美日韩一区二区三区在线| 一夜夜www| 久久青草综合色| 精品视频人人做人人爽| 在线看a的网站| 在线播放国产精品三级| av线在线观看网站| 最新的欧美精品一区二区| 国产精品电影一区二区三区 | 日韩成人在线观看一区二区三区| 久久久久国内视频| 国产成人精品久久二区二区免费| 天天影视国产精品| 免费看十八禁软件| 国产av国产精品国产| 韩国精品一区二区三区| 久久中文字幕人妻熟女| 国产精品美女特级片免费视频播放器 | 国产精品免费一区二区三区在线 | 国产精品二区激情视频| bbb黄色大片| 免费在线观看视频国产中文字幕亚洲| 欧美变态另类bdsm刘玥| 久久这里只有精品19| 国产老妇伦熟女老妇高清| 久久毛片免费看一区二区三区| 两性夫妻黄色片| 丝袜美足系列| 两个人免费观看高清视频| 一区二区av电影网| 一级黄色大片毛片| 免费观看人在逋| 免费在线观看日本一区| 国产欧美日韩一区二区三| 中文字幕色久视频| 亚洲精品国产区一区二| 久久久久久久精品吃奶| 久久国产亚洲av麻豆专区| 免费在线观看完整版高清| 欧美一级毛片孕妇| 少妇粗大呻吟视频| cao死你这个sao货| 真人做人爱边吃奶动态| 亚洲成人免费av在线播放| 又大又爽又粗| 欧美精品啪啪一区二区三区| 国产精品一区二区在线观看99| 男女高潮啪啪啪动态图| 少妇精品久久久久久久| 69精品国产乱码久久久| 一级a爱视频在线免费观看| 日本a在线网址| 亚洲av成人一区二区三| 亚洲,欧美精品.| 免费久久久久久久精品成人欧美视频| 亚洲av成人不卡在线观看播放网| 国产成人影院久久av| 精品高清国产在线一区| 亚洲国产看品久久| 精品国产一区二区久久| 久久久久网色| 淫秽高清视频在线观看| 亚洲国产精品久久男人天堂| 国产av在哪里看| 岛国在线免费视频观看| 日韩大尺度精品在线看网址| 美女免费视频网站| 99久久无色码亚洲精品果冻| 国产又黄又爽又无遮挡在线| 欧美日本亚洲视频在线播放| 美女cb高潮喷水在线观看 | 人人妻,人人澡人人爽秒播| 男人的好看免费观看在线视频| 91在线观看av| 欧美黄色片欧美黄色片| e午夜精品久久久久久久| 动漫黄色视频在线观看| 99国产精品一区二区蜜桃av| 日日夜夜操网爽| 国产91精品成人一区二区三区| 中出人妻视频一区二区| 午夜福利在线在线| 一级作爱视频免费观看| 国产极品精品免费视频能看的| 国产成人精品无人区| 亚洲中文日韩欧美视频| 日本精品一区二区三区蜜桃| 亚洲第一电影网av| 久久久国产欧美日韩av| 一进一出抽搐动态| 国产高潮美女av| 久久久色成人| 久久国产精品影院| 制服丝袜大香蕉在线| 天堂影院成人在线观看| 亚洲成人精品中文字幕电影| 黄片大片在线免费观看| 免费无遮挡裸体视频| 90打野战视频偷拍视频| 亚洲国产中文字幕在线视频| 美女免费视频网站| aaaaa片日本免费| 免费在线观看影片大全网站| 色精品久久人妻99蜜桃| 精品无人区乱码1区二区| 国产精品亚洲av一区麻豆| 久久精品人妻少妇| 日本一二三区视频观看| 网址你懂的国产日韩在线| 美女扒开内裤让男人捅视频| 午夜福利在线观看免费完整高清在 | 亚洲av第一区精品v没综合| 亚洲午夜理论影院| 哪里可以看免费的av片| 午夜精品一区二区三区免费看| 精品国产三级普通话版| 特大巨黑吊av在线直播| 亚洲国产精品999在线| 1024手机看黄色片| 精品不卡国产一区二区三区| a级毛片在线看网站| 欧美绝顶高潮抽搐喷水| 国产精品自产拍在线观看55亚洲| 蜜桃久久精品国产亚洲av| 丁香欧美五月| 欧美一级毛片孕妇| 国产真人三级小视频在线观看| 国内精品久久久久久久电影| 男人舔女人的私密视频| 国产欧美日韩一区二区精品| 国产午夜福利久久久久久| 十八禁人妻一区二区| 中文亚洲av片在线观看爽| 亚洲五月天丁香| 免费无遮挡裸体视频| 久久中文字幕一级| 窝窝影院91人妻| 国产精品,欧美在线| 99热这里只有精品一区 | 国产精品精品国产色婷婷| 中文资源天堂在线| 国产精品精品国产色婷婷| 小蜜桃在线观看免费完整版高清| 床上黄色一级片| 中文字幕久久专区| 88av欧美| 国产在线精品亚洲第一网站| 岛国视频午夜一区免费看| 日韩欧美国产在线观看| 热99在线观看视频| 国产精品 欧美亚洲| 黄色女人牲交| 亚洲第一电影网av| 国产高清激情床上av| 久久久国产成人精品二区| 国产三级在线视频| aaaaa片日本免费| 男女午夜视频在线观看| av片东京热男人的天堂| 国产一区二区三区视频了| 给我免费播放毛片高清在线观看| 日韩中文字幕欧美一区二区| 老鸭窝网址在线观看| 亚洲精品一卡2卡三卡4卡5卡| 嫩草影院入口| 黑人巨大精品欧美一区二区mp4| 一个人观看的视频www高清免费观看 | av黄色大香蕉| 国产精品久久久人人做人人爽| 成人三级做爰电影| 男人舔女人下体高潮全视频| 人妻久久中文字幕网| 热99在线观看视频| 伦理电影免费视频| 18禁国产床啪视频网站| 最近最新中文字幕大全电影3| 久久精品国产亚洲av香蕉五月| 午夜成年电影在线免费观看| 国产一区在线观看成人免费| 国产精品亚洲av一区麻豆| 中文字幕精品亚洲无线码一区| 日韩欧美国产一区二区入口| 99热6这里只有精品| 亚洲精品在线观看二区| 色视频www国产| 美女免费视频网站| 精品99又大又爽又粗少妇毛片 | 亚洲18禁久久av| 99久久精品一区二区三区| 色综合站精品国产| 国产精品亚洲美女久久久| 久久久久久久精品吃奶| 精品不卡国产一区二区三区| 午夜福利视频1000在线观看| 色哟哟哟哟哟哟| 国产真实乱freesex| 日韩精品青青久久久久久| 欧美中文日本在线观看视频| 欧美在线黄色| 亚洲av电影不卡..在线观看| 国产免费av片在线观看野外av| 国产综合懂色| www.999成人在线观看| 亚洲国产欧美网| 男女之事视频高清在线观看| 国产黄片美女视频| 十八禁人妻一区二区| 亚洲成人久久性| 国产成人一区二区三区免费视频网站| 久久久久亚洲av毛片大全| 非洲黑人性xxxx精品又粗又长| 好男人电影高清在线观看| 在线免费观看不下载黄p国产 | 亚洲,欧美精品.| 99热精品在线国产| 日韩欧美免费精品| 免费无遮挡裸体视频| 亚洲色图 男人天堂 中文字幕| 男人舔奶头视频| 国产亚洲欧美98| 久久精品aⅴ一区二区三区四区| 91老司机精品| 欧美性猛交╳xxx乱大交人| 免费无遮挡裸体视频| 久久久国产成人免费| 欧美av亚洲av综合av国产av| 美女高潮的动态| 久久精品aⅴ一区二区三区四区| 久久久久精品国产欧美久久久| 亚洲精品美女久久av网站| netflix在线观看网站| а√天堂www在线а√下载| 欧美精品啪啪一区二区三区| 美女cb高潮喷水在线观看 | 嫁个100分男人电影在线观看| 国产亚洲精品久久久久久毛片| 操出白浆在线播放| 在线视频色国产色| 级片在线观看| 久久久久久大精品| 757午夜福利合集在线观看| 国产精品久久视频播放| 久久国产精品影院| 成人性生交大片免费视频hd| 美女免费视频网站| 两个人看的免费小视频| 欧美一区二区国产精品久久精品| 国产精品99久久久久久久久| 色精品久久人妻99蜜桃| 欧美一区二区精品小视频在线| 97超级碰碰碰精品色视频在线观看| 日本 av在线| 怎么达到女性高潮| 免费大片18禁| av天堂中文字幕网| bbb黄色大片| 欧美又色又爽又黄视频| 精品久久蜜臀av无| 国产精品久久久久久精品电影| 欧美最黄视频在线播放免费|