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

    重載鐵路曲線鋼軌廓形多目標優(yōu)化設計

    2021-01-29 03:20:58唐彥玲吳磊董勇陳帥王衡禹
    機械 2020年12期
    關鍵詞:廓形輪軌鋼軌

    唐彥玲,吳磊*,,董勇,陳帥,王衡禹

    重載鐵路曲線鋼軌廓形多目標優(yōu)化設計

    唐彥玲1,吳磊*,1,董勇1,陳帥2,王衡禹3

    (1.西南交通大學 機械工程學院,四川 成都 610031;2.西南交通大學 土木工程學院,四川 成都 610031;3.西南交通大學 牽引動力國家重點實驗室,四川 成都 610031)

    針對重載鐵路曲線軌道鋼軌磨耗和壓潰問題,建立了一種基于遺傳算法和層次分析法的鋼軌廓形多目標優(yōu)化設計方法。在設計方法中,將實測重載曲線廓形作為遺傳優(yōu)化初始種群,利用改進的NURBS曲線參數化鋼軌廓形,用C80貨車車輛-軌道動力學模型仿真的磨耗指數、接觸應力、脫軌系數、輪重減載率和平穩(wěn)性指標作為優(yōu)化目標,運用層次分析法確立各指標權重,通過遺傳迭代求解出最佳曲線非對稱鋼軌廓形。輪軌接觸分析表明,該鋼軌廓形在降低輪軌磨耗和降低輪軌接觸應力方面都非常有效。

    重載鐵路;非對稱鋼軌廓形;遺傳算法;層次分析法;車輛動力學

    由于重載線路貨運量大,小半徑曲線多,并受軌道、車輛和地理環(huán)境等因素影響,重載鐵路曲線鋼軌易產生嚴重側磨、壓潰肥邊、疲勞裂紋、踏面剝離掉塊等諸多病害。通過輪軌型面優(yōu)化來改善輪軌接觸關系已成為降低輪軌磨耗、緩解疲勞損傷的重要手段[1-2]。

    國內外眾多學者致力于鐵路輪軌型面優(yōu)化設計研究。Heller等[3]提出了一種基于高斯迭代的車輪踏面求解方法,他們是較早運用了基于車輛動力學的多目標型面優(yōu)化方法,但由于多個設計目標相互獨立,容易導致各目標間無法相互兼顧。Persson等[4-5]以最大接觸應力、最大輪軌橫向力、最大脫軌系數、磨耗數等動力學相關的罰因子加權和為目標函數,運用遺傳算法對輪軌廓形進行求解,并取得了良好的改善效果。Shevtsov等[6-7]提出基于滾動圓半徑差函數反求車輪踏面的設計方法,以減小磨耗和同時考慮降低滾動接觸疲勞為目標獲取了新的車輪踏面。成棣[8]針對不同的車型和線路確定優(yōu)化的多個動力學指標為優(yōu)化目標,基于高斯徑向基函數響應面的數值求解方法得到預期效果的車輪踏面。崔檬[9]以非線性臨界速度和磨耗指數為優(yōu)化目標,基于NSGA-Ⅱ遺傳算法和Pareto多目標優(yōu)化求解方法求得最優(yōu)車輪踏面。崔大賓等[10]以保證車輛直線運行穩(wěn)定性與減小通過小半徑曲線時的輪緣磨耗為目標,運用改進的粒子群優(yōu)化算法求解了最優(yōu)車輪廓形。綜上,通過多目標優(yōu)化設計方法進行型面優(yōu)化可以兼顧多個性能指標的變化,從而控制優(yōu)化過程中一些潛在因素的影響。

    目前車輪踏面的多目標設計領域已有了豐富的經驗和基礎研究方法,隨著重載鐵路的快速發(fā)展,鐵路曲線鋼軌的多目標設計研究還面臨著更多的挑戰(zhàn)。已有的多目標鋼軌廓形設計方法大多針對特定車型和線路,利用綜合評價策略從所設計的多個型面中選擇出最優(yōu)的鋼軌廓形,很難獲得全局最優(yōu)的鋼軌型面[11-12]。并且鐵路曲線鋼軌非對稱磨耗突出,在鋼軌廓形設計中考慮曲線鋼軌的非對稱性尤為重要。本文主要從優(yōu)化設計方法的建立和應用展開研究,以國內某重載線路為例,基于現場輪軌廓形調查數據,明確從重載曲線輪軌磨耗、接觸應力以及車輛運行動力學性能出發(fā)進行優(yōu)化,考慮了子目標之間的相互影響,曲線鋼軌廓形優(yōu)化過程中采用了非對稱設計,且運用了遺傳算法來獲得能夠繼承實測廓形優(yōu)良特性的重載曲線鋼軌廓形。

    1 計算模型

    本文運用SIMPACK軟件基于C80重載貨車參數建立車輛-軌道系統(tǒng)動力學模型,并利用MATLAB建立的多目標遺傳迭代數值優(yōu)化模型對鋼軌廓形進行優(yōu)化,同時設置了優(yōu)化迭代的初始條件參數。

    1.1 動力學模型

    為了研究輪軌匹配時的車輛動力學性能,本文利用SIMPACK軟件建立了C80重載貨車模型。模型只建立了單節(jié)貨車車輛,忽略相鄰車輛對它的影響,該車輛主要由車體、側架、搖枕、輪對等部分組成,且不考慮這些部件的彈性變形,考慮輪軌接觸、懸掛系統(tǒng)的非線性。由于車輛前后轉向架極具對稱性,圖1給出了一半車輛動力學模型的關系拓撲圖,車輛的動力學方程可表示為:

    1.2 優(yōu)化模型

    本文基于遺傳優(yōu)化算法和層次分析法建立曲線軌道鋼軌廓形多目標優(yōu)化迭代模型,并利用MATLAB軟件編寫鋼軌廓形優(yōu)化數值程序。

    圖1 車輛軌道動力學模型關系拓撲圖

    1.2.1 設計變量

    根據文獻[13]發(fā)展的NURBS曲線方法對鋼軌廓形進行參數化,僅考慮以廓形接觸區(qū)域控制點的縱坐標和切線斜率作為設計變量,得到:

    式中:R為曲線高軌(左軌)和低軌(右軌)的設計變量組成的列向量;為控制點;為鋼軌廓形的控制點數;=1為高軌(左軌)廓形;=2為低軌(右軌)廓形。

    考慮到車輛過曲線時輪對的橫向偏移,將曲線鋼軌設計為非對稱廓形時的性能較優(yōu),因此將高軌和低軌的設計變量合成為向量為:

    式中:為非對稱鋼軌廓形設計變量(多變量問題);1和2分別為高軌(左軌)和低軌(右軌)。

    1.2.2 綜合目標函數

    (1)評價指標

    現場調研發(fā)現,重載線路曲線鋼軌多出現高軌嚴重側磨和低軌壓潰肥邊等現象,考慮到車輛運行最基本的安全性和貨物運輸的完整性,將輪軌磨耗、接觸應力、脫軌系數、輪重減載率和平穩(wěn)性作為優(yōu)化的多個目標,并將這多個目標的加權和作為優(yōu)化的綜合目標函數。

    考慮了自旋蠕滑,則有:

    式中:W為單個接觸斑在某時刻的磨耗指數,N;TT分別為輪軌接觸面上的縱向蠕滑力、橫向蠕滑力和自旋蠕滑力,N;ξξ分別為輪軌接觸面上的縱向蠕滑率、橫向蠕滑率和自旋蠕滑率。

    式中:1()為輪軌間磨耗數均值指標;為車輪數量。

    輪軌間接觸法向壓應力利用最大赫茲接觸應力來表示:

    式中:maxi為單個接觸斑在某時刻的最大接觸應力,N/m2;N為法向力的彈性值,N;和分別為接觸斑縱向半軸長和橫向半軸長,m。

    式中:2()為輪軌間最大接觸應力均值指標;maxj為所有輪軌間最大接觸應力。

    為保證車輛運行的安全性,則有:

    式中:3()為所有輪對輪軌間輸出的脫軌系數最大值指標,且3()<0.8;為脫軌系數。

    式中:4()為所有輪對輪重減載率最大值指標,且4()<0.8;為輪重減載率。

    式中:5()為時域計算所得的橫向平穩(wěn)性最大值指標;6()為時域計算所得的垂向平穩(wěn)性最大值指標;S為橫向平穩(wěn)性;S為垂向平穩(wěn)性。

    (2)歸一化處理

    對于多個評價指標,由于量綱的不同可能會造成數據分析的失誤,因此需要將各個指標做歸一化處理,通過數學處理方法將其轉換為0~1的標量[14]。

    利用常用的極值法處理,則:

    根據文獻[15],已知脫軌系數最大限值為0.8、動態(tài)輪重減載率最大限值為0.8,則得到式(14)、式(15)。

    根據GB 9955-85[16],貨車平穩(wěn)性一級指標為3.5,指標最大限度為4.25,則:

    (3)指標權重

    為了獲取綜合評價指標,需要對各評價指標進行加權,本文基于20世紀70年代初期美國運籌學家Saaty[17]提出的層次分析法理論,建立了目標層和指標層的最優(yōu)鋼軌廓形個體評價層次結構,如圖2所示。

    圖2 評價指標層次結構模型

    結構搭建好之后,需要對指標層的各評價指標進行兩兩比較,確定相對重要性。采用Saaty[17]等提出的引用數字1~9及其倒數作為標度的方法進行分析,標度的含義如表1所示。

    1=0.0834;2=0.1049;3=0.3218;

    4=0.2810;5=0.1298;6=0.0790

    則得:

    式中:()為綜合目標函數;為子目標個數。

    表1 因子對比標度[17]

    表2 專家群組判斷矩陣表

    1.2.3 廓形優(yōu)化迭代流程

    近年來遺傳算法在優(yōu)化問題中得到了廣泛應用,本文運用遺傳算法的思想來得到全局搜索的最優(yōu)鋼軌廓形,如圖3所示。

    圖3 廓形優(yōu)化迭代流程

    首先確定初始種群,種群由一定數目帶有染色體特征的個體組成,本文初始種群為一批帶有不同特征性狀的實測鋼軌廓形。將鋼軌廓形進行參數化統(tǒng)一變量格式,應用建立的車輛-軌道動力學模型計算各動力學指標,對各指標進行歸一化,利用層次分析法確定各指標權重,從而得到綜合目標函數。由綜合目標函數計算種群中每個個體的適應度,適應度即為決定種群個體中被選入下一代的概率,本文中的適應度即為最小化的綜合目標函數()的倒數。通過判斷終止條件來決定遺傳優(yōu)化是否結束,若不滿足條件則按概率選擇種群中的個體進行選擇、交叉、變異等遺傳操作生成下一代,如此往復,若滿足條件則終止迭代,輸出最優(yōu)鋼軌廓形。

    1.3 初始條件

    為了更高效地實現全局最優(yōu)鋼軌廓形,需要確立計算的基準模型工況、初始種群以及鋼軌廓形遺傳迭代的終止條件。

    (1)基準工況

    本文針對的是重載鐵路曲線的鋼軌廓形優(yōu)化,參照實測線路,模型的線路設置為圓曲線,半徑800 m,長度500 m,前后緩和曲線長度150 m,前后直線200 m,采用車輛運行速度為80 km/h的均衡超高,施加的軌道不平順為北美重載實測軌道不平順激勵。

    (2)初始種群

    為了提高算法的速度和保持鋼軌廓形的多樣性,共選取了18組初始鋼軌廓形,如圖4所示,每組廓形均為重載曲線線路不同測點的實測鋼軌廓形。其中編號為1的這一組廓形為實測新鋪上道的標準60 kg/m鋼軌(CHN60),編號1、6、11的這三組為隨機測試的直線廓形,其余組為實測的不同曲線半徑的鋼軌廓形。

    圖4 初始種群

    (3)迭代終止條件

    終止條件決定了廓形優(yōu)化流程何時結束。目前遺傳算法判斷終止的方法不一,常用的終止條件有以下三種:指定遺傳算法迭代次數;種群個體最小偏差,即群體中最優(yōu)個體與最小個體適應度偏差小于某一極小的閾值;適應度趨勢,即種群中連續(xù)幾代最優(yōu)個體的適應度增大趨勢平穩(wěn)則停止[18-19]。

    本文的適應度為綜合目標函數的倒數,經過對目標函數值的多次指定迭代次數試算,綜合上述適應度趨勢的終止條件判斷方法,最終規(guī)定迭代次數為200代,當迭代滿足設置的迭代次數時,則終止迭代計算,輸出最優(yōu)解。

    2 優(yōu)化方法的應用

    應用建立的多目標優(yōu)化計算模型,研究初始廓形總體和最佳個體的各動力學指標值的迭代演化過程是驗證優(yōu)化方法可行的關鍵。分析最優(yōu)鋼軌廓形的靜態(tài)接觸幾何關系和動力學性能是判斷優(yōu)化廓形性能是否優(yōu)秀的重要依據。

    2.1 鋼軌廓形迭代演化過程

    鋼軌廓形優(yōu)化過程中,輸入初始條件參數,經過優(yōu)化模型的不斷迭代,得到了鋼軌廓形個體的各項指標的迭代演化規(guī)律。圖5給出了每代個體磨耗指數和接觸應力的演變情況。

    由圖5(a)可知,一位導向輪對右輪間輪軌磨耗較其他輪對更為嚴重,通過多目標優(yōu)化之后,一位輪對右輪從初期到第20代磨耗指數下降顯著,其他輪對間的輪軌磨耗在第3代就明顯下降,之后變化平緩。與實測廓形相比,優(yōu)化后輪軌間磨耗指數總體下降了約50.0%。由圖5(b)可知,隨著迭代次數的增加,右側車輪(低軌)接觸應力顯著下降,左側車輪(高軌)接觸應力變化幅值較小。與實測廓形相比,優(yōu)化后輪軌間接觸應力總體下降了約34.2%。其脫軌系數、輪重減載率和平穩(wěn)性指標影響不大,均在標準限制以內。

    2.2 優(yōu)化廓形性能分析

    本文利用所建立的數值優(yōu)化模型迭代求解出目標鋼軌廓形,為驗證目標設計廓形的可行性,將設計鋼軌廓形與標準LM車輪匹配,并與標準CHN60鋼軌和實測磨耗鋼軌進行對比,分析其靜態(tài)接觸幾何性能和動力學性能。

    (1)靜態(tài)接觸幾何分析

    輪軌接觸點對的分布可以反映輪對不同橫移量下接觸點的位置,圖6給出了優(yōu)化設計廓形、實測廓形和標準廓形的接觸點情況,圖中橫線上數字表示輪對的橫移量,車輪與鋼軌踏面之間的連線代表不同橫移量下在車輪和鋼軌上的接觸點位置。由圖6可以看出,優(yōu)化之后的高軌接觸點對更集中于軌頂中心,低軌在向右橫移時在軌頂外側也有接觸點分布,優(yōu)化設計低軌廓形接觸點分布范圍加寬,有利于減小低軌接觸應力集中,緩解滾動接觸疲勞。

    圖5 每一代個體動力學指標的演化

    圖6 標準CHN60、實測鋼軌、優(yōu)化鋼軌接觸點對對比

    圖7為優(yōu)化設計廓形與標準LM車輪匹配時的等效錐度與實測廓形和標準鋼軌廓形的對比。等效錐度是被廣泛應用于鐵道車輛的重要參數,用來表征輪軌幾何線性化匹配關系[20]。在通過曲線線路時適當增大等效錐度可以提升車輛過曲線性能,但較大的等效錐度不利于直線運行穩(wěn)定性。國際鐵路聯盟標準UIC 519[21]規(guī)定了輪對橫移幅值3 mm為名義等效錐度,國際鐵路聯盟標準UIC 518[22]規(guī)定了速度小于200 km/h的等效錐度限值為0.4。由圖7可以看出,橫移幅值為3 mm時,優(yōu)化設計廓形的等效錐度為0.12,比實測廓形增大了約33.6%,比標準廓形增大了約17.6%。在橫移9 mm之內,優(yōu)化設計廓形的等效錐度都略高于標準廓形,且明顯高于實測廓形,因此,優(yōu)化設計廓形有利于改善曲線通過性能。

    圖7 標準CHN60、實測鋼軌、優(yōu)化鋼軌等效錐度對比

    (2)動力學性能分析

    圖8給出了優(yōu)化設計廓形、實測廓形和標準廓形曲線通過時的導向輪對橫移量。從圖中可以看出,優(yōu)化設計廓形的輪對橫移量整體大于標準廓形,因此,優(yōu)化設計廓形在過曲線時能夠提供較大的輪對橫移量,利于車輛的曲線通過。實測廓形雖能夠提供更大的輪對橫移量,但從圖6的接觸點對可以發(fā)現,橫移超過9 mm時,容易發(fā)生軌距角處與輪緣貼靠,加劇軌頭內側與輪緣根部的磨耗。因此優(yōu)化設計廓形盡可能地提升了車輛曲線通過能力又減小了輪軌磨耗。

    圖9給出了不同速度下優(yōu)化設計廓形與標準LM車輪匹配時的動態(tài)特性。由圖9(a)可知,隨著速度的增加,優(yōu)化鋼軌和標準鋼軌廓形的磨耗指數呈指數增加,實測鋼軌基本呈線性增長,從整體來看優(yōu)化廓形的磨耗指數均低于實測廓形和標準廓形。由圖9(b)可知,隨著速度的增加,優(yōu)化設計廓形的接觸應力與標準鋼軌和實測鋼軌廓形的變化趨勢相同,不同運行速度下,設計廓形的輪軌接觸應力皆低于標準廓形,且顯著低于實測鋼軌廓形。

    圖8 輪對橫移量對比

    圖9 不同運行速度下動態(tài)特性對比

    3 結論

    根據國內某重載鐵路曲線實測鋼軌廓形磨耗情況及標準CHN60鋼軌型面進行了曲線鋼軌廓形設計。建立了基于遺傳算法和層次分析法,以減小輪軌磨耗和接觸應力、提高車輛運行安全性和平穩(wěn)性的多目標重載曲線鋼軌廓形設計方法。研究結論如下:

    (1)重載曲線多目標鋼軌廓形設計模型包括了車輛-軌道動力學子模型用來輸出優(yōu)化目標函數參數、多目標子模型基于層次分析法用以確立目標權重分布、參數化子模型基于改進的NURBS曲線方法進行鋼軌廓形重構、鋼軌廓形迭代子模型基于遺傳算法最終搜索到全局最優(yōu)目標廓形。

    (2)應用所建立的重載曲線鋼軌多目標優(yōu)化設計方法獲得的目標廓形與實測廓形相比,其輪軌間磨耗指數和接觸應力都有了明顯改進;其安全性指標和平穩(wěn)性指標變化幅值較小,且在標準限制以內。

    (3)對優(yōu)化設計廓形的接觸幾何和動力學性能分析可知,優(yōu)化后的鋼軌廓形有利于車輛曲線通過;所研究的車輛運行速度范圍內,優(yōu)化廓形的磨耗和接觸應力皆低于實測廓形和標準廓形。

    [1]王璞. 重載鐵路輪軌磨耗預測及鋼軌型面優(yōu)化研究[D]. 北京:北京交通大學,2017.

    [2]鐘浩. 基于改善輪軌接觸狀態(tài)的重載車輪型面優(yōu)化[D]. 成都:西南交通大學,2014.

    [3]Heller R,Law E H. Optimizing the wheel profile to improve rail vehicle dynamic performance[J]. Vehicle System Dynamics,1980,8(2-3):116-122.

    [4]Persson I,Iwnicki S D. Optimization of railway wheel profiles using a genetic algorithm[J]. Vehicle System Dynamics,2004(41):517-526.

    [5]Persson I,Nilsson R,Bik U,et al. Use of a genetic algorithm to improve the rail profile on Stockholm underground[J]. Vehicle System Dynamics,2010,48(S1):89-104.

    [6]Shevtsov I Y,Markine V L,Esveld C. Optimal design of wheel profile for railway vehicles[J]. Wear,2004,258(7):1022-1030.

    [7]Shevtsov I Y,Markine V L,Esveld C. Design of railway wheel profile taking into account rolling contact fatigue and wear[J]. Wear,2008,265(9-10):1273-1282.

    [8]成棣. 車輪型面多目標優(yōu)化設計研究[D]. 北京:中國鐵道科學研究院,2011.

    [9]崔檬. 車輪踏面外形的多目標動態(tài)優(yōu)化[D]. 北京:北京交通大學,2016.

    [10]崔大賓,趙銳,李俊,等. 面向輪緣磨耗的動車組車輪型面多目標優(yōu)化設計[J]. 機械工程學報,2019,55(16):105-113.

    [11]肖杰靈,劉學毅. 鋼軌非對稱廓形的設計方法[J]. 西南交通大學學報,2010,45(3):362-365.

    [12]李怡然. 鋼軌廓形對機車車輛動力學性能的影響研究[D]. 成都:西南交通大學,2018.

    [13]Wu L,Liu J Q,Dhanasekar M,et al. Optimization of railhead profiles for curved tracks using improved nonuniform rational B-splines and measured profiles[J]. Wear,2019(418-419):123-132.

    [14]成棣,胡曉依,劉豐收,等. 基于改進層次分析法的高速列車輪軌型面匹配評價方法及應用[J]. 中國鐵道科學,2019,40(3):81-87.

    [15]TB 10761-2013,高速鐵路工程動態(tài)驗收技術規(guī)范[S].

    [16]GB/T 5599-1985,鐵道車輛動力學性能評定和試驗鑒定規(guī)范[S].

    [17]Saaty T L. Decision making with the analytic hierarchy process [J]. International journal of services sciences,2008,1(1):83-98.

    [18]王蓮芬. 層次分析法引論[M]. 北京:中國人民大學出版社,1990.

    [19]周明,孫樹棟. 遺傳算法原理及應用[M]. 北京:國防工業(yè)出版社,1999.

    [20]趙艷杰,楊東曉. 高速動車組踏面與不同鋼軌匹配關系研究[J]. 機械,2020,47(1):57-63.

    [21]UIC Code 519,Method for determining the equivalent conicity[S]. 2004.

    [22]UIC Code 518,Testing and approval of railway vehicles from the point of view of their dynamic behavior-Safety-Track fatigue-Ride quality[S]. 2005.

    Multi-Objective Optimization Design of Rail Profiles on Curve Tracks of Heavy Haul Railway

    TANG Yanling1,WU Lei1,DONG Yong1,CHEN Shuai2,WANG Hengyu3

    ( 1.School of Mechanical Engineering, Southwest Jiaotong University, Chengdu 610031, China; 2.School of Civil Engineering, Southwest Jiaotong University, Chengdu 610031, China; 3.State Key Laboratory of Traction Power, Southwest Jiaotong University, Chengdu 610031, China)

    Aiming at solving the problems of rail wear and crushing of heavy haul railway curve track, a multi-objective optimization design method for rail profile based on genetic algorithm and analytic hierarchy process was developed. The measured profiles of heavy haul railway curve track were used as the initial population for genetic optimization, and the improved NURBS curve was used to parameterize the rail profile. The wear index, contact stress, derailment coefficient, wheel load reduction rate, and Sperling index of the simulated C80 vehicle dynamics model were used as optimization targets. The analytic hierarchy process is used to establish the weight of each index. Finally, the optimal curved asymmetric rail profile is solved by genetic iteration. The wheel-rail contact analysis shows that this rail profile is very effective in reducing wheel-rail wear and wheel-rail contact stress.

    heavy haul railway;asymmetric rail profile;genetic algorithm;analytic hierarchy process;vehicle dynamics

    U270.1+1

    A

    10.3969/j.issn.1006-0316.2020.12.001

    1006-0316 (2020) 12-0001-09

    2020-05-12

    國家自然科學基金(51605395、51775454)

    唐彥玲(1994-),女,四川南充人,碩士研究生,主要研究方向為輪軌關系。*通訊作者:吳磊(1981-),男,貴州畢節(jié)人,博士,講師,主要研究方向為輪軌關系,E-mail:wuleitpl@163.com。

    猜你喜歡
    廓形輪軌鋼軌
    意在輕閑
    中低速磁浮道岔與輪軌道岔的差異
    力量廓形
    窄氣感廓形
    鋼軌焊后雙中頻感應加熱工藝研究
    大衣+牛仔褲 走路得帶風
    中低速磁浮與輪軌交通信號系統(tǒng)的差異
    非線性穩(wěn)態(tài)曲線通過時輪軌滾動接觸的數值求解方法
    高速鐵路鋼軌疲勞過程的超聲非線性系數表征
    國內外高速鐵路鋼軌性能對比研究
    又大又爽又粗| 免费黄频网站在线观看国产| 精品第一国产精品| 嫁个100分男人电影在线观看| av天堂久久9| 操出白浆在线播放| 啪啪无遮挡十八禁网站| 日本欧美视频一区| 久久人人97超碰香蕉20202| 国产精品久久久久久精品古装| 五月天丁香电影| 国产老妇伦熟女老妇高清| 夜夜骑夜夜射夜夜干| 精品卡一卡二卡四卡免费| 中文字幕av电影在线播放| 我的亚洲天堂| 精品卡一卡二卡四卡免费| 夜夜骑夜夜射夜夜干| 男男h啪啪无遮挡| 黄色视频在线播放观看不卡| av在线老鸭窝| 亚洲欧美精品自产自拍| 成人国语在线视频| 在线观看舔阴道视频| 精品国产超薄肉色丝袜足j| 精品亚洲乱码少妇综合久久| 成年人免费黄色播放视频| 亚洲精华国产精华精| 亚洲免费av在线视频| 国产精品一区二区免费欧美 | av欧美777| 亚洲av日韩在线播放| 久久久精品国产亚洲av高清涩受| 91老司机精品| 国产精品久久久久久精品古装| 一区二区日韩欧美中文字幕| cao死你这个sao货| 91麻豆av在线| 亚洲欧美精品综合一区二区三区| 男女无遮挡免费网站观看| 国产男女超爽视频在线观看| 国产精品秋霞免费鲁丝片| 一级片'在线观看视频| 侵犯人妻中文字幕一二三四区| 精品亚洲成a人片在线观看| 黄色视频不卡| 精品一区二区三区四区五区乱码| 中文字幕高清在线视频| 嫩草影视91久久| 美女高潮到喷水免费观看| 免费人妻精品一区二区三区视频| 国产一区二区 视频在线| 欧美精品亚洲一区二区| h视频一区二区三区| 欧美日韩亚洲综合一区二区三区_| 国产一区二区三区av在线| 国产精品 欧美亚洲| 中文字幕精品免费在线观看视频| 国产精品麻豆人妻色哟哟久久| 成年动漫av网址| 国产福利在线免费观看视频| av又黄又爽大尺度在线免费看| 人妻 亚洲 视频| 国产在线视频一区二区| 久9热在线精品视频| 中文字幕最新亚洲高清| 亚洲五月色婷婷综合| 啦啦啦免费观看视频1| 老司机影院毛片| 菩萨蛮人人尽说江南好唐韦庄| av在线播放精品| 纯流量卡能插随身wifi吗| 亚洲av美国av| 人人妻人人澡人人看| 久久影院123| 欧美成人午夜精品| 免费在线观看完整版高清| 伊人亚洲综合成人网| 国产成人欧美在线观看 | 午夜福利在线免费观看网站| e午夜精品久久久久久久| 啦啦啦在线免费观看视频4| 蜜桃国产av成人99| 久久女婷五月综合色啪小说| 国产精品亚洲av一区麻豆| 精品卡一卡二卡四卡免费| 18禁观看日本| videosex国产| 777米奇影视久久| 精品福利观看| 91精品三级在线观看| 国产精品偷伦视频观看了| 亚洲avbb在线观看| 侵犯人妻中文字幕一二三四区| h视频一区二区三区| cao死你这个sao货| 亚洲欧美日韩高清在线视频 | 美女高潮喷水抽搐中文字幕| 欧美日韩视频精品一区| 大片免费播放器 马上看| 老汉色∧v一级毛片| 国产伦理片在线播放av一区| 男人操女人黄网站| 人妻人人澡人人爽人人| 精品视频人人做人人爽| 丝袜人妻中文字幕| 中文字幕av电影在线播放| 午夜精品国产一区二区电影| 午夜福利免费观看在线| 首页视频小说图片口味搜索| 在线十欧美十亚洲十日本专区| 国产亚洲精品第一综合不卡| 高清黄色对白视频在线免费看| 亚洲国产欧美在线一区| 精品福利观看| 黄色视频不卡| 十分钟在线观看高清视频www| 女人爽到高潮嗷嗷叫在线视频| av线在线观看网站| 天天操日日干夜夜撸| 十分钟在线观看高清视频www| 国产成人精品久久二区二区免费| 丁香六月天网| 这个男人来自地球电影免费观看| 亚洲欧洲日产国产| 精品少妇一区二区三区视频日本电影| 亚洲精品一二三| 国产精品偷伦视频观看了| 青春草亚洲视频在线观看| 亚洲欧美日韩另类电影网站| 一区二区三区激情视频| h视频一区二区三区| 国产成人精品无人区| 亚洲成国产人片在线观看| 亚洲成人手机| 精品人妻1区二区| 欧美xxⅹ黑人| 亚洲欧洲精品一区二区精品久久久| 九色亚洲精品在线播放| 久久亚洲国产成人精品v| 超色免费av| 黑人欧美特级aaaaaa片| 亚洲精品中文字幕在线视频| 精品人妻在线不人妻| 日本vs欧美在线观看视频| 欧美激情高清一区二区三区| 国产免费av片在线观看野外av| 国产又色又爽无遮挡免| 久久国产精品大桥未久av| 日韩大码丰满熟妇| 99国产综合亚洲精品| 俄罗斯特黄特色一大片| 老鸭窝网址在线观看| 国产极品粉嫩免费观看在线| 精品欧美一区二区三区在线| 黄色a级毛片大全视频| 男人爽女人下面视频在线观看| 亚洲国产成人一精品久久久| 美女中出高潮动态图| 成人亚洲精品一区在线观看| 色精品久久人妻99蜜桃| 精品久久久久久电影网| 国产人伦9x9x在线观看| 女性被躁到高潮视频| 黑人巨大精品欧美一区二区蜜桃| 一区二区日韩欧美中文字幕| 免费观看a级毛片全部| 丁香六月欧美| 国产真人三级小视频在线观看| 日本a在线网址| 国产成人一区二区三区免费视频网站| 老司机午夜福利在线观看视频 | 亚洲精品一二三| 日韩人妻精品一区2区三区| 性少妇av在线| 啦啦啦在线免费观看视频4| 国产日韩欧美视频二区| 欧美97在线视频| 大型av网站在线播放| 欧美成狂野欧美在线观看| 久久久国产一区二区| av电影中文网址| 黑人巨大精品欧美一区二区蜜桃| 十八禁网站网址无遮挡| 久久天躁狠狠躁夜夜2o2o| 波多野结衣一区麻豆| www.精华液| 欧美成人午夜精品| 日本精品一区二区三区蜜桃| 汤姆久久久久久久影院中文字幕| 久久久久国产一级毛片高清牌| 午夜日韩欧美国产| 十八禁高潮呻吟视频| 亚洲成av片中文字幕在线观看| 丝袜美腿诱惑在线| 汤姆久久久久久久影院中文字幕| 中文欧美无线码| 91大片在线观看| 亚洲成人免费电影在线观看| 日韩电影二区| 黄网站色视频无遮挡免费观看| 亚洲欧美色中文字幕在线| 国产成人精品在线电影| 成年人黄色毛片网站| 国产成人啪精品午夜网站| 精品卡一卡二卡四卡免费| 国产高清视频在线播放一区 | 久久久国产一区二区| 91精品伊人久久大香线蕉| 亚洲欧美日韩另类电影网站| 高清视频免费观看一区二区| 色婷婷av一区二区三区视频| 亚洲自偷自拍图片 自拍| 两个人免费观看高清视频| 欧美精品一区二区免费开放| 高清欧美精品videossex| 热99久久久久精品小说推荐| 成人国产av品久久久| 蜜桃国产av成人99| 久久免费观看电影| 热99国产精品久久久久久7| 久久中文看片网| 性少妇av在线| 国产欧美日韩一区二区三区在线| 亚洲精品美女久久av网站| 美女大奶头黄色视频| 丝袜人妻中文字幕| bbb黄色大片| 亚洲成人手机| 午夜福利视频在线观看免费| 免费观看人在逋| 国产精品免费大片| 日韩一区二区三区影片| 91成人精品电影| 欧美日本中文国产一区发布| 国产精品成人在线| 91老司机精品| 亚洲av成人不卡在线观看播放网 | 久久ye,这里只有精品| 久久久久网色| 18禁国产床啪视频网站| 色精品久久人妻99蜜桃| 久久热在线av| 12—13女人毛片做爰片一| 考比视频在线观看| a级毛片在线看网站| 2018国产大陆天天弄谢| 脱女人内裤的视频| 老司机靠b影院| 美女脱内裤让男人舔精品视频| 亚洲黑人精品在线| 国产男女超爽视频在线观看| 欧美日韩一级在线毛片| 我的亚洲天堂| 在线亚洲精品国产二区图片欧美| 中文字幕人妻丝袜一区二区| 大香蕉久久成人网| 精品国产乱子伦一区二区三区 | 满18在线观看网站| 国产精品欧美亚洲77777| 狂野欧美激情性bbbbbb| av视频免费观看在线观看| 伊人亚洲综合成人网| 欧美国产精品一级二级三级| 午夜激情av网站| 宅男免费午夜| 最黄视频免费看| 欧美日韩亚洲综合一区二区三区_| 极品少妇高潮喷水抽搐| 成年人免费黄色播放视频| av福利片在线| 亚洲自偷自拍图片 自拍| tube8黄色片| 亚洲一区二区三区欧美精品| 久久久久精品国产欧美久久久 | 日本五十路高清| 嫁个100分男人电影在线观看| 久久国产精品人妻蜜桃| 免费在线观看影片大全网站| 考比视频在线观看| 性少妇av在线| 亚洲专区字幕在线| 国产精品香港三级国产av潘金莲| 久久久久国内视频| 国产亚洲午夜精品一区二区久久| 久久久久网色| 日日爽夜夜爽网站| 日韩中文字幕视频在线看片| 午夜久久久在线观看| 一级毛片女人18水好多| 久久国产精品影院| 人妻一区二区av| 午夜福利在线观看吧| 大码成人一级视频| 国产成人精品久久二区二区免费| 精品国产国语对白av| 999久久久国产精品视频| 成年人免费黄色播放视频| 亚洲人成电影免费在线| 王馨瑶露胸无遮挡在线观看| 丝瓜视频免费看黄片| 免费高清在线观看日韩| 少妇人妻久久综合中文| 咕卡用的链子| 久热这里只有精品99| a在线观看视频网站| 一级毛片女人18水好多| 嫩草影视91久久| 69av精品久久久久久 | 老司机午夜福利在线观看视频 | 国产免费一区二区三区四区乱码| 男女床上黄色一级片免费看| 男女国产视频网站| 午夜福利乱码中文字幕| 国产精品亚洲av一区麻豆| 成人av一区二区三区在线看 | 少妇粗大呻吟视频| 欧美日韩视频精品一区| 国产精品免费视频内射| 男女下面插进去视频免费观看| 少妇猛男粗大的猛烈进出视频| 蜜桃在线观看..| √禁漫天堂资源中文www| av线在线观看网站| 国产一级毛片在线| 久久精品国产亚洲av高清一级| 啦啦啦在线免费观看视频4| 国产欧美日韩一区二区三 | 乱人伦中国视频| 久久久国产成人免费| 欧美日韩av久久| 亚洲av日韩在线播放| 少妇猛男粗大的猛烈进出视频| 久久中文看片网| av一本久久久久| h视频一区二区三区| 久久久久久亚洲精品国产蜜桃av| 狠狠精品人妻久久久久久综合| 免费看十八禁软件| 亚洲欧美激情在线| 中文欧美无线码| 亚洲熟女毛片儿| 中亚洲国语对白在线视频| 咕卡用的链子| 免费女性裸体啪啪无遮挡网站| 国精品久久久久久国模美| 宅男免费午夜| 黄色a级毛片大全视频| 久久精品国产亚洲av香蕉五月 | 18禁黄网站禁片午夜丰满| 丰满饥渴人妻一区二区三| 亚洲一码二码三码区别大吗| av不卡在线播放| 国产av精品麻豆| 丰满人妻熟妇乱又伦精品不卡| 精品少妇一区二区三区视频日本电影| 亚洲精品自拍成人| 一本大道久久a久久精品| 久久亚洲精品不卡| 国产片内射在线| 国产精品国产av在线观看| 色综合欧美亚洲国产小说| 女人高潮潮喷娇喘18禁视频| 一区二区av电影网| 91成人精品电影| 成人国产av品久久久| h视频一区二区三区| 在线av久久热| 少妇精品久久久久久久| 亚洲av男天堂| 久久国产精品男人的天堂亚洲| 久久这里只有精品19| 每晚都被弄得嗷嗷叫到高潮| 国产欧美日韩一区二区精品| 丝袜在线中文字幕| 99精品久久久久人妻精品| 国产亚洲av高清不卡| 伊人亚洲综合成人网| 老司机影院成人| 成年美女黄网站色视频大全免费| 超碰成人久久| av超薄肉色丝袜交足视频| 我要看黄色一级片免费的| 午夜福利在线免费观看网站| 亚洲av片天天在线观看| 成人亚洲精品一区在线观看| 精品高清国产在线一区| 一个人免费在线观看的高清视频 | 黄片小视频在线播放| 国产激情久久老熟女| 99久久人妻综合| 日本wwww免费看| 新久久久久国产一级毛片| 精品一区二区三区四区五区乱码| 乱人伦中国视频| 国产一区二区在线观看av| 久久人人97超碰香蕉20202| 69精品国产乱码久久久| 免费在线观看影片大全网站| 美女主播在线视频| 精品亚洲乱码少妇综合久久| 国产一区二区在线观看av| 久久人人97超碰香蕉20202| 操出白浆在线播放| 精品人妻1区二区| 国产成人欧美在线观看 | 亚洲精品国产av蜜桃| 99香蕉大伊视频| 日韩制服丝袜自拍偷拍| 午夜老司机福利片| 国产野战对白在线观看| 少妇的丰满在线观看| 好男人电影高清在线观看| av网站免费在线观看视频| 精品欧美一区二区三区在线| 老司机午夜十八禁免费视频| 国产在视频线精品| a级片在线免费高清观看视频| 久久性视频一级片| 欧美精品啪啪一区二区三区 | 亚洲欧美一区二区三区黑人| 精品一区在线观看国产| 亚洲伊人色综图| 俄罗斯特黄特色一大片| 国产精品一区二区免费欧美 | 欧美日韩福利视频一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 天堂8中文在线网| 欧美黑人精品巨大| 最黄视频免费看| 国产亚洲av高清不卡| 亚洲人成77777在线视频| 中国国产av一级| 国产xxxxx性猛交| 日本撒尿小便嘘嘘汇集6| 亚洲中文字幕日韩| 欧美精品人与动牲交sv欧美| 99久久综合免费| 捣出白浆h1v1| 精品卡一卡二卡四卡免费| 国产精品一区二区在线观看99| 中文字幕高清在线视频| 欧美激情极品国产一区二区三区| 午夜福利视频精品| 51午夜福利影视在线观看| 在线 av 中文字幕| 国产黄色免费在线视频| 久热爱精品视频在线9| 国产成人精品久久二区二区91| 咕卡用的链子| 最近最新免费中文字幕在线| 少妇 在线观看| 亚洲天堂av无毛| 丝瓜视频免费看黄片| 夜夜夜夜夜久久久久| av在线老鸭窝| 丰满人妻熟妇乱又伦精品不卡| 欧美精品啪啪一区二区三区 | 国产1区2区3区精品| 亚洲av日韩精品久久久久久密| 国产成人免费无遮挡视频| 亚洲av电影在线进入| 国产亚洲精品一区二区www | 高清视频免费观看一区二区| 50天的宝宝边吃奶边哭怎么回事| 美女高潮喷水抽搐中文字幕| 精品一区二区三区av网在线观看 | 51午夜福利影视在线观看| a级毛片黄视频| 亚洲五月婷婷丁香| 韩国精品一区二区三区| 国产97色在线日韩免费| 美女午夜性视频免费| 美女扒开内裤让男人捅视频| 考比视频在线观看| 丁香六月欧美| 亚洲一码二码三码区别大吗| 亚洲国产精品成人久久小说| 欧美日韩成人在线一区二区| 婷婷色av中文字幕| 亚洲伊人久久精品综合| 国产亚洲精品第一综合不卡| 中文精品一卡2卡3卡4更新| h视频一区二区三区| 欧美少妇被猛烈插入视频| 国产不卡av网站在线观看| 国产一区二区 视频在线| 久久久水蜜桃国产精品网| av在线播放精品| 肉色欧美久久久久久久蜜桃| 一级a爱视频在线免费观看| a在线观看视频网站| 欧美 日韩 精品 国产| 三级毛片av免费| 欧美97在线视频| 人人妻人人添人人爽欧美一区卜| 五月开心婷婷网| 我要看黄色一级片免费的| 亚洲av日韩在线播放| 丝袜脚勾引网站| 欧美人与性动交α欧美精品济南到| 乱人伦中国视频| 国产激情久久老熟女| 91麻豆精品激情在线观看国产 | 人人妻人人爽人人添夜夜欢视频| 91精品国产国语对白视频| 免费观看人在逋| 亚洲中文字幕日韩| 夜夜夜夜夜久久久久| 精品视频人人做人人爽| 天堂8中文在线网| 99国产综合亚洲精品| 成年人午夜在线观看视频| 1024视频免费在线观看| 国产一级毛片在线| 婷婷丁香在线五月| 亚洲国产日韩一区二区| 国产免费一区二区三区四区乱码| 天天操日日干夜夜撸| 性色av一级| 首页视频小说图片口味搜索| 久久久欧美国产精品| 久久久久久久久免费视频了| 在线天堂中文资源库| 亚洲av成人一区二区三| 久久99一区二区三区| 天堂俺去俺来也www色官网| 操出白浆在线播放| 一区二区日韩欧美中文字幕| 国产99久久九九免费精品| 久久人妻福利社区极品人妻图片| 亚洲专区国产一区二区| 在线观看免费午夜福利视频| 免费黄频网站在线观看国产| 免费久久久久久久精品成人欧美视频| 国产欧美日韩精品亚洲av| 啦啦啦在线免费观看视频4| 丝袜美腿诱惑在线| 午夜两性在线视频| 成人国产av品久久久| 久久中文字幕一级| 另类精品久久| 色综合欧美亚洲国产小说| 久久av网站| 一本综合久久免费| 制服诱惑二区| 精品视频人人做人人爽| 久久99热这里只频精品6学生| 久久人妻福利社区极品人妻图片| 精品少妇一区二区三区视频日本电影| 两人在一起打扑克的视频| 午夜福利一区二区在线看| 一区二区三区四区激情视频| 日韩欧美免费精品| 韩国高清视频一区二区三区| av片东京热男人的天堂| 一区二区日韩欧美中文字幕| 日日爽夜夜爽网站| 中亚洲国语对白在线视频| 老汉色∧v一级毛片| 各种免费的搞黄视频| 国产xxxxx性猛交| 丁香六月欧美| 免费看十八禁软件| 丝袜喷水一区| 久久青草综合色| 国产成人精品无人区| 日韩中文字幕欧美一区二区| e午夜精品久久久久久久| 久久久欧美国产精品| 啦啦啦在线免费观看视频4| 国产免费福利视频在线观看| 不卡一级毛片| 日本猛色少妇xxxxx猛交久久| 一二三四社区在线视频社区8| 91九色精品人成在线观看| 又大又爽又粗| 午夜影院在线不卡| 久久国产亚洲av麻豆专区| 久久精品国产亚洲av高清一级| 不卡一级毛片| 法律面前人人平等表现在哪些方面 | 一级片'在线观看视频| 十八禁网站网址无遮挡| 在线十欧美十亚洲十日本专区| 如日韩欧美国产精品一区二区三区| 少妇精品久久久久久久| 久久久久久久久免费视频了| 午夜福利在线免费观看网站| 一本综合久久免费| 中亚洲国语对白在线视频| 亚洲欧美日韩另类电影网站| 男女边摸边吃奶| 啦啦啦啦在线视频资源| 一级,二级,三级黄色视频| 岛国在线观看网站| 丝袜美腿诱惑在线| a级毛片黄视频| 在线看a的网站| 一二三四在线观看免费中文在| 人成视频在线观看免费观看| 曰老女人黄片| 久久性视频一级片| 国产深夜福利视频在线观看| 性少妇av在线| 久久精品aⅴ一区二区三区四区| 高清在线国产一区| 岛国毛片在线播放| 欧美日韩国产mv在线观看视频| 国产一区二区三区在线臀色熟女 | videos熟女内射| 丰满人妻熟妇乱又伦精品不卡| 亚洲av美国av| 亚洲av成人一区二区三| 国产亚洲精品久久久久5区|