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

    考慮地鐵列車控制的縱斷面設計優(yōu)化方法

    2025-02-13 00:00:00樊蔥
    山東科學 2025年1期
    關鍵詞:縱斷面磨損量輪軌

    摘要:地鐵列車控制優(yōu)化與縱斷面設計關系密切,二者對運營費用有重要影響。為在列車控制優(yōu)化的基礎上進一步降低運營費用,在線路縱斷面設計階段構(gòu)建了協(xié)同優(yōu)化模型,考慮定時列車控制對縱斷面進行優(yōu)化,在滿足定時列車控制約束、《地鐵設計規(guī)范》要求等的前提下,最小化能耗成本和維修成本之和。由于車輪、鋼軌維修成本的影響因素眾多,建立列車-軌道動力學仿真模型計算維修成本,在此基礎上設計綜合偽譜法和蠻力搜索的算法求解協(xié)同優(yōu)化方案,并應用于廣州地鐵線路的3個區(qū)間,對方法的優(yōu)化效果進行驗證。研究結(jié)果表明,對比在實際縱斷面上單獨優(yōu)化定時列車控制的方法,同時優(yōu)化定時列車控制方案和縱斷面方案對節(jié)約運營費用的效果更佳,各區(qū)間運營費用平均降低21%。該研究方法可為進一步控制運營費用提供新思路和理論依據(jù),推動地鐵系統(tǒng)可持續(xù)發(fā)展。

    關鍵詞:地鐵;列車操縱策略;線路設計;運營成本;偽譜法

    中圖分類號:U212.34""" 文獻標志碼:A""" 文章編號:1002-4026(2025)01-0105-15

    開放科學(資源服務)標志碼(OSID):

    DOI:10.3976/j.issn.1002-4026.20240062【交通運輸】

    收稿日期:2024-04-16

    基金項目:國家自然科學基金(71571016)

    作者簡介:樊蔥(1996—),女,碩士,助理工程師,研究方向為城軌線路設計。E-mail: cong19960922@163.com

    Optimizing vertical track alignment considering metro train control

    FAN Cong

    (Guangdong Guangzhou-Zhuhai Intercity Railway Co., Ltd., Guangzhou 510422, China)

    Abstract∶Energy-saving metro train control is closely related to the vertical track alignment (VTA) design, and both have a significant impact on operating costs. To further reduce operating costs based on optimized train control, this study proposed a collaborative optimization model for the VTA design phase. This model optimizes the bidirectional train control strategy and VTA of a metro section with the goal of minimizing energy consumption and maintenance costs simultaneously, while adhering to the constraints of scheduled train control and the requirements of the \"Metro Design Code.\" Given the numerous factors affecting the maintenance costs of wheels and rails, a train-track dynamic simulation model was developed to calculate these costs. Based on this, an algorithm combining the pseudospectral method and brute force search was designed to solve the collaborative optimization model. The effectiveness of this optimization method was validated using three sections of the Guangzhou metro line. The results indicate that, compared to the method of optimizing scheduled train control alone on the actual VTA, the collaborative optimization model is more effective in saving operating costs, reducing the average operating costs by 21% across the studied sections. This study can provide novel approaches and theoretical support to further reduce metro operating costs, which contributes to promoting sustainable development of metro.

    Key words∶underground metro; train trajectory; track alignment design; operating cost; pseudospectral method

    發(fā)展地鐵具有方便市民出行、帶動區(qū)域經(jīng)濟等優(yōu)點,隨著城市化進程持續(xù)加快,近五年地鐵線路總里程增長約4 500 km[1]。地鐵線路增長導致能耗成本問題突出,近五年地鐵耗電量增長逾4億度,其中50%~70%用于列車牽引[2]。另一方面,隨著建成年限變長,更多線路進入維修期,維修成本不斷增加,根據(jù)北京地鐵數(shù)據(jù),每公里的維修成本甚至超過了能耗成本。

    為降低列車牽引能耗,節(jié)約能耗成本,地鐵運營單位采取了優(yōu)化列車控制、使用輕量化車輛等措施[3],但這些措施的節(jié)能效果受線路條件制約。地鐵線路主要包括線路平面和縱斷面。平面方案設計更側(cè)重于盡量多地吸引客流[4],而非降低列車能耗??v斷面直接影響列車所受阻力,不同縱斷面下的列車控制優(yōu)化方案及能耗可能不同,而且以節(jié)能為目標優(yōu)化縱斷面時需考慮列車控制優(yōu)化,因為最優(yōu)縱斷面方案可能隨列車控制方案改變。所以,相比平面設計,縱斷面設計與列車控制優(yōu)化的關聯(lián)更緊密,如果在設計縱斷面時同步優(yōu)化列車控制,相較僅優(yōu)化列車控制方案,將更有助于節(jié)約能耗成本。

    縱斷面設計一般指在給定線路平面的條件下,確定滿足實際約束的一系列坡道。Jha等[5]首先在設計縱斷面時同時考慮降低建設成本和能耗成本,但假設能耗成本正比于列車運行距離,計算方法較為粗略。為精確計算能耗,Kim等[6]建立列車運動模型描述列車運行行為,導致模型的非線性化增強,當考慮更多實際約束時需使用更高效的算法求解。柏赟等[7]、彭磊等[8]相對全面地考慮了縱斷面設計中的實際約束,并采用啟發(fā)式算法快速求解模型,但解的最優(yōu)性難以保障。徐楊[9]進一步設計蠻力搜索和動態(tài)規(guī)劃相結(jié)合的算法求得問題的精確解,然而研究中列車采用節(jié)時控制策略,所求解不一定最適用于實際運營中給定運行時分下的列車控制。Hoang等[4]和Wu等[10]在已知運行時分的前提下,同時優(yōu)化縱斷面方案和列車節(jié)能控制方案,但可能受計算效率的限制,為縮小問題可行解規(guī)模,優(yōu)化過程中縱斷面設計形式保持不變,解的最優(yōu)性無法確保。Duarte等[11]對算法進行改進,設計梯度優(yōu)化法求解協(xié)同優(yōu)化模型,但忽略了坡長約束等,可能導致所求方案無法在實際中應用。

    除能耗成本外,不同列車節(jié)能控制方案和縱斷面方案下,設備使用壽命通常存在差異,從而影響維修成本。湯超[12]、肖緋雄等[13]采用動力學仿真方法計算不同坡度、列車運行速度下的輪軌、閘瓦磨損量,證明在其他條件不變的前提下,輪軌和閘瓦磨損量隨坡度絕對值或列車運行速度增大而增加。Lai等[14]在優(yōu)化縱斷面和列車控制時考慮了節(jié)約維修費用,但維修費用僅正比于運行距離,未考慮不同設備間維修成本的差異。Kang等[15]分別根據(jù)車輪、鋼軌等設備的年維修費用計算各維修成本,再以所有設備維修成本及能耗成本最小為優(yōu)化目標,求解最優(yōu)縱斷面方案和列車控制方案,但仍未較精準地量化不同方案對維修成本的影響。

    在更精準化不同方案對維修成本影響的基礎上,本文考慮定時列車控制對縱斷面進行優(yōu)化,提出基于列車-軌道動力學仿真的地鐵區(qū)間縱斷面設計和列車節(jié)能控制協(xié)同優(yōu)化方法,該方法盡可能全面的考慮定時列車控制約束和《地鐵設計規(guī)范》[16]等要求,采用偽譜法和蠻力搜索的算法,在保證求解精度的前提下高效求解優(yōu)化方案,以最小化能耗成本和維修成本。

    1" 考慮列車控制的縱斷面優(yōu)化問題

    截至2022年底,全國城市軌道交通中69%以上線路為地下線[1],地鐵中地下線占比更高,廣州地鐵達88%左右,多數(shù)使用盾構(gòu)法建設。此時,建設費用基本取決于線路長度,由平面方案決定[17]。由于在縱斷面設計階段,平面方案一般已提前確定,不同縱斷面方案下的線路長度及相應的建設成本差異很小。因此,通過優(yōu)化地下線列車控制方案和縱斷面方案,可以在不改變建設成本的前提下降低能耗成本和維修成本,達到減少總成本的目的。

    對于地下線車站,其水平位置由平面方案給定,高程一般按盡量接近地面以方便乘客進出站的原則確定,車站的空間位置可視為已知條件。此時,當僅考慮一般的運營模式,列車在各站均停車,任意兩相鄰車站間的列車控制和縱斷面優(yōu)化與其他區(qū)間互不干擾,各區(qū)間可單獨進行優(yōu)化。為提高問題的求解效率,以下對單個區(qū)間上考慮列車控制的縱斷面優(yōu)化問題進行研究,當全部區(qū)間都使用該方法,即可得到全線優(yōu)化方案。

    第n個區(qū)間縱斷面設計示意圖如圖1所示,一般情況下,上行和下行方向軌道敷設在同一隧道中,縱斷面呈對稱關系,不同方向上的最優(yōu)列車控制方案可能存在差異。若協(xié)同優(yōu)化僅以降低一個方向的列車牽引能耗為目標,可能使另一方向的列車牽引能耗較大,不利于節(jié)約總能耗成本。因此,計算能耗成本時考慮雙向牽引能耗。

    由于地鐵系統(tǒng)復雜、設備眾多,維修成本僅考慮與縱斷面方案、列車控制方案密切相關的重要設備,主要有閘瓦、車輪和鋼軌。當各設備單次平均維修成本一定,在相同時間內(nèi),總維修成本由各設備使用壽命決定。閘瓦、車輪和鋼軌使用到限的主要原因是摩擦磨損[18],所以維修成本為修理、更換因摩擦磨損使用到限閘瓦、車輪和鋼軌的人工成本和材料成本。雖然設備磨損量還受軌道結(jié)構(gòu)影響,但在實際中,軌道結(jié)構(gòu)設計主要考慮減振降噪[19],而非減少設備磨損,在研究中視為已知條件。

    為最小化能耗成本和維修成本之和,對縱斷面方案和列車定時控制方案進行協(xié)同優(yōu)化。其中,縱斷面可由各坡度和坡長、連接相鄰坡道的豎曲線表示。由于豎曲線位于相鄰坡道交匯處,且曲線半徑按設計規(guī)范要求取值,決定縱斷面方案的關鍵因素為各坡度和坡長。在給定的線路上,列車定時控制方案由牽引力和制動力決定。優(yōu)化方案同時受定時列車控制約束、《地鐵設計規(guī)范》要求和實際工程約束等限制,具體將在模型部分進行介紹。

    2" 協(xié)同優(yōu)化模型構(gòu)建

    對于協(xié)同優(yōu)化問題,線路條件及部分參數(shù)如限速和避讓區(qū)域都隨位移發(fā)生變化,將這些變量稱為位移關聯(lián)參數(shù)。為簡化問題中對位移關聯(lián)參數(shù)的表示,將區(qū)間按一定規(guī)則劃分為多個階段,構(gòu)建多階段優(yōu)化模型,以提高求解效率[20]。階段劃分如圖2所示,首先將每個曲線區(qū)域和避讓區(qū)域都設為單獨的階段,再按同一階段內(nèi)坡度和限速都相同的原則進行劃分。由于限速、曲線和避讓區(qū)域位置為已知值,確定各坡道長度后,即可獲得階段劃分方案。

    2.1" 模型假設

    模型假設如下:

    (1)研究表明,再生制動能利用與縱斷面設計關聯(lián)較?。?],本文計算能耗成本時暫不考慮再生制動能;

    (2)假設運營單位對設施設備狀態(tài)進行定期監(jiān)測,能夠及時更換磨損到限的設備,且不考慮單次平均維修成本定價問題,其不因設備維修次數(shù)等發(fā)生改變;

    (3)相比輪軌在曲線上的磨損,其在直線上的磨損非常?。?1],因此只考慮輪軌在曲線上的磨損,且同一列車的各車輪、各鋼軌磨損量相同。

    2.2" 優(yōu)化目標和決策變量

    由第1節(jié)分析可知,縱斷面方案由K個坡道的坡度i(k)和長度l(k)(k=1,2,…,K)確定,列車操縱方案由位置s處的牽引力F1(s)、F2(s)和制動力B1(s)、B2(s)確定(下標1、2分別代表上行方向和下行方向),以上要素為決策變量。

    通過優(yōu)化縱斷面方案和列車操縱方案,一定研究時段內(nèi)的最小化運營成本Cmin,見公式(1):

    Cmin=CE+CS+CW+CA,(1)

    式中,CE、CS、CW和CA分別為給定研究時段內(nèi)產(chǎn)生的能耗成本、閘瓦維修成本、車輪維修成本和鋼軌維修成本。

    2.2.1" 能耗成本計算

    能耗成本為列車運行產(chǎn)生的牽引能耗相應的用電費用,若已知單位電價為cJ,能耗成本的計算公式為:

    CE=cJEP,(2)

    式中,E為列車雙向牽引能耗;P為給定研究時間內(nèi)列車在線路上往返運行的次數(shù)。根據(jù)牽引能耗定義,E的計算公式見公式(3):

    E=∑Rr=1E(r)1+∑Rr=1E(r)2,(3)

    式中,R為劃分的階段數(shù)量;E(r)1和E(r)2分別為上行、下行方向第r階段的牽引能耗,二者計算方法類似,以上行方向為例,E(r)1的計算方法見公式(4):

    E(r)1=∫s(r)fs(r)0F1(s(r))ds,(4)

    式中,s(r)0和s(r)f分別為第r階段起點和終點位置。

    2.2.2" 閘瓦維修成本計算

    當閘瓦厚度減少至一定限度,需更換閘瓦,由此產(chǎn)生的成本為閘瓦維修成本。已知在研究時段內(nèi),列車在線路上往返運行P次,所需閘瓦維修成本見公式(5):

    CS=「PUScSNS,(5)

    式中,「表示向上取整數(shù),cS、NS分別為單塊閘瓦價格、列車閘瓦數(shù)量;US為閘瓦厚度減少到限時,列車在線路上往返運行的最大次數(shù),其計算公式為:

    US=ImS(IS1+IS2),(6)

    式中,ImS為閘瓦最大磨損量;IS1和IS2為列車在上行、下行方向運行過程中閘瓦厚度減少量,由于二者計算方法類似,以IS1為例說明具體計算方法。在列車制動過程中,閘瓦與車輪接觸使動能轉(zhuǎn)換為熱能,列車速度下降,同時閘瓦溫度升高使閘瓦變薄,發(fā)生磨損,根據(jù)能量轉(zhuǎn)換原理,IS1的計算見公式(7):

    IS1=∑Rr=1∫t1(s(r)f)t1(s(r)0)βγεMIa1(t)v1(t)NISIdt,(7)

    式中,β為0-1變量,當制動力大于0時為1,否則為0;γ、ε分別為動能轉(zhuǎn)換為熱能的比例[22]、閘瓦吸收熱能的比例[13];t1(s)為列車在位置s處的運行時刻;MI為車軸質(zhì)量;a1(t)、v1(t)分別為加速度和速度;NI、SI分別為車輛閘瓦數(shù)量和每塊閘瓦的面積。

    2.2.3" 車輪維修成本計算

    車輪磨損源于車輪和鋼軌間的相對滑動[23],表現(xiàn)為輪緣厚度減少,車輪型面變得不規(guī)則。當車輪型面變化到一定程度時,不利于列車安全運行[24],通常采用鏇輪方式恢復,但鏇輪將導致車輪踏面直徑減小且無法恢復,踏面直徑減少到限時,需更換車輪。因此,可用踏面直徑衡量車輪磨損程度,計算維修成本CW,見公式(8):

    CW=「PUWcWNW,(8)

    式中,cW、NW分別為單個車輪的維修價格、列車車輪數(shù)量;UW為車輪踏面直徑減少到限時,列車在線路上往返運行的次數(shù),其計算公式為:

    UW=ImW(IW1+IW2),(9)

    式中,ImW為車輪踏面直徑最大磨損量;IW1、IW2為列車在上行、下行方向運行過程中車輪踏面直徑減少量,由于二者計算方式類似,以上行方向為例,IW1可由輪緣厚度減少量換算,見公式(10):

    IW1=∑Rr=1∫s(r)fs(r)0θIL1(s(r))ds,(10)

    式中,θ為踏面直徑減少量與輪緣厚度減少量的比例,為一經(jīng)驗常數(shù);IL1(s(r))為列車在上行方向運行過程中輪緣厚度減少量。由于車輪與鋼軌接觸作用過程復雜,且輪緣厚度減少量受線路、列車、軌道相關的諸多參數(shù)影響[23],難以用數(shù)值方法計算。為此,采用仿真模型與數(shù)值方法相結(jié)合的方式,先基于仿真軟件建立列車-軌道模型,獲得輪軌間相互作用關系,再使用車輪磨耗模型計算輪緣厚度減少量,具體過程將在第3節(jié)介紹。

    2.2.4" 鋼軌維修成本計算

    輪軌間相對滑動也將導致鋼軌磨損,主要表現(xiàn)為波磨和側(cè)磨兩種形式[25]。其中,波磨產(chǎn)生原因復雜、難以預測,實際中采用定期打磨方式進行維護,本文研究的鋼軌磨損為鋼軌側(cè)磨。側(cè)磨主要發(fā)生在線路曲線上,按照《鋼軌允許磨耗限度》[26],以磨損深度為側(cè)磨限度的衡量指標。當鋼軌磨損深度到限時,對鋼軌進行維修產(chǎn)生的成本為鋼軌維修成本CA,計算方法見公式(11):

    CA=「PUAcANA,(11)

    式中,cA、NA分別為單條鋼軌的維修價格、鋼軌數(shù)量;UA為鋼軌磨損深度到限時,列車在線路上往返運行的次數(shù),其計算公式為:

    UA=ImA∑Rr1=1∫s(r)fs(r)0IA1(s(r))ds+∑Rr2=1∫s(r)fs(r)0IA2(s(r))ds,(12)

    式中,ImA為鋼軌磨損深度最大值;IA1(s(r))、IA2(s(r))為列車在上行、下行方向運行過程中第r階段的鋼軌磨損深度,與輪緣厚度減少量類似,鋼軌磨損深度也需采用仿真模型與數(shù)值方法相結(jié)合的方式計算。

    2.3" 約束條件

    協(xié)同優(yōu)化問題需同時考慮縱斷面設計和定時列車控制兩類約束。由于列車在上行和下行方向上運行時所受約束相似,以上行方向為例對模型約束進行介紹。

    2.3.1" 縱斷面設計約束

    根據(jù)線路設計規(guī)范和工程建設要求,縱斷面設計主要受車站及站外區(qū)域坡長、坡度約束,夾直線長度約束,車站高程約束,覆土厚度約束和避讓區(qū)域約束的限制。

    (1)車站及站外區(qū)域坡長、坡度約束

    為保障列車運行平穩(wěn),避免跨坡道運行引起的列車振動,車站應盡量設置在一個坡道上,即車站區(qū)域進站和出站坡長之和大于站臺長度,同時進站和出站坡度均采用滿足排水要求的最小坡度id。假設車站進站和出站坡長相等,則車站區(qū)域坡長、坡度約束為:

    l(k)≥Ls/2,(k=1,K),(13)

    i(k)=id,(k=1,K),(14)

    式中,Ls為站臺長度。

    站外區(qū)域坡長約束與車站區(qū)域類似,坡長應大于列車長度Lt,而坡度不僅要滿足排水要求,還不能大于規(guī)定值im,確保列車運行安全,見公式(15~16):

    l(k)≥Lt,(k=2,3,…,(K-1)),(15)

    id≤i(k)≤im,(k=2,3,…,(K-1))。(16)

    (2)夾直線長度約束

    根據(jù)地鐵設計規(guī)范要求,相鄰豎曲線間需使用不小于規(guī)定長度Lm的直線連接,確保線路平順,有利于列車安全運行。連接相鄰豎曲線的直線段為夾直線,如圖3所示,夾直線長度約束可以表示為公式(17):

    l(k)-L(k-1,k)-L(k,k+1)≥Lm,(17)

    式中,L(k-1,k)為第k-1個坡道與第k個坡道間豎曲線的切線長度,其計算公式為:

    L(k-1,k)=Rki(k-1)-i(k)2 000,(18)

    式中,Rk為第k-1個與第k個坡道間豎曲線的半徑。

    結(jié)合以上兩個公式,可以將夾直線長度約束式(17)表示為公式(19):

    i(k-1)-i(k)+i(k)-i(k+1)≤2 000(l(k)-Lm)Rk。(19)

    (3)覆土厚度約束

    對于地下線,線路上方應保留厚度至少為G的覆土,避免地鐵線路影響地面設施,該約束稱為覆土厚度約束,見公式(20):

    H(s(r))lt;Hg(s(r))-G,(20)

    式中,H(s(r))、Hg(s(r))分別為位置s(r)處的線路高程、地面高程,地面高程為已知條件。線路高程計算公式為:

    H(s(r))=∑r-1m=1i(m)(s(m)f-s(m)0)+i(r)(s(r)-s(r)0)。(21)

    (4)車站高程約束

    由第1節(jié)的分析可知,在協(xié)同優(yōu)化問題中,車站高程固定,可表示為公式(22):

    H(s(1)0)=H0,H(s(R)f)=Hf,(22)

    式中,H0和Hf為起點和終點車站的高程。

    (5)避讓區(qū)域約束

    避讓區(qū)域為地下設施和不良水文地質(zhì)所在區(qū)域,在設計縱斷面時應繞行,該約束稱為避讓區(qū)域約束。對于在第p階段的避讓區(qū)域,用Q(p)表示避讓區(qū)域所有的點,則約束可表示為公式(23):

    H(s(p))-Hq(s(p))gt;0,(q∈Q(p)),(23)

    式中,Hq(s(p))為第p階段避讓區(qū)域中q點的高程。

    2.3.2" 定時列車控制約束

    (1)定時運行約束

    定時列車控制為在給定區(qū)間運行時分下,采用適當?shù)臓恳椭苿恿Γ沽熊嚭馁M最少的牽引能耗從區(qū)間起點運行至終點。因此,列車在起點和終點的位置、運行時分、速度均為固定值,該約束稱為定時運行約束,見公式(24~26):

    s1(s(1)0)=S0,s1(s(R)f)=Sf,(24)

    t1(s(1)0)=0,t1(s(R)f)=T,(25)

    v1(s(1)0)=v1(s(R)f)=0,(26)

    式中,S0和Sf分別為起點、終點水平位置;T為給定的區(qū)間運行時間,計算方法見公式(27):

    T=3.6(Sf-S0)vc,(27)

    式中,vc為技術速度。

    (2)區(qū)間速度約束

    為確保行車安全,列車在區(qū)間運行過程中不得超速運行,應與限定速度保持一定差距Δv,即區(qū)間速度約束為:

    0lt;v1(s(r))≤vm(s(r))-Δv,(28)

    式中,vm(s(r))為在s(r)處的限定速度。

    (3)牽引力和制動力約束

    根據(jù)列車牽引制動特性,在同一時刻,列車最多只能使用牽引力或制動力中的一種,且所用牽引力或制動力不能超過限定值,可表示為公式(29~30):

    0≤F1(s(r))≤Fm(s(r)),0≤B1(s(r))≤Bm(s(r)),(29)

    F1(s(r))B1(s(r))=0,(30)

    式中,F(xiàn)m(s(r))、Bm(s(r))為最大牽引力、最大制動力。

    為精確計算使用的牽引力、制動力,需結(jié)合列車受力分析、牛頓第二定律等,建立列車運動方程,描述列車運行行為。列車運動方程可表示為:

    dv1(s(r))ds=F1(s(r))-B1(s(r))-W1(s(r))-N1(s(r))M(1+ρ)v1(s(r)),(31)

    dt1(s(r))ds=1v1(s(r)),(32)

    式中,W1(s(r))、N1(s(r))分別為列車運行過程中所受的基本阻力和附加阻力,附加阻力包括坡道附加阻力和曲線附加阻力;M為列車質(zhì)量;ρ為回轉(zhuǎn)質(zhì)量系數(shù),為與列車類型有關的常數(shù);v1(s(r))、t1(s(r))分別為列車速度和運行時間?;咀枇Α⒏郊幼枇Φ挠嬎惴椒椋?/p>

    W1(s(r))=Mga+bv1(s(r))+cv1(s(r))2,(33)

    N1(s(r))=Mgi(s(r))+600L(s(r))R(s(r))Lt,(34)

    式中,g為重力加速度;a、b、c為與列車性質(zhì)有關的常量;i(s(r))為坡度;R(s(r))為線路曲線半徑;L(s(r))為列車在曲線上的長度;Lt為列車長度。

    (4)階段連接約束

    由于模型將區(qū)間劃分為多個階段,為得到連續(xù)的列車操縱方案,需對相鄰階段交匯點的速度、時間、位置進行約束,以獲得完整的列車操縱方案,該約束稱為階段連接約束。階段連接約束為第r階段終點和第r+1階段起點的速度、時間、位置相等,見公式(35~37):

    s(r)f=s(r+1)0,(1≤r≤R-1),(35)

    v1(s(r)f)=v1(s(r+1)0),(1≤r≤R-1),(36)

    t1(s(r)f)=t1(s(r+1)0),(1≤r≤R-1)。(37)

    3" 算法設計

    以上多階段優(yōu)化模型具有非線性性強、變量和約束多等特點,求解難度大。列車控制優(yōu)化研究對比了偽譜法、動態(tài)規(guī)劃和混合整數(shù)規(guī)劃算法的求解精度和效率,結(jié)果表明偽譜法具有更高的求解精度、更快的求解速度,同時變量和約束的數(shù)量對偽譜法求解速度影響很小,適用求解多階段優(yōu)化模型[27]。因此,本文采用偽譜法求解多階段優(yōu)化模型,其中階段劃分由蠻力搜索算法(BF算法)完成,車輪、鋼軌維修成本基于列車-軌道動力學仿真計算。

    3.1" 算法框架

    本文算法基于列車-軌道動力學仿真、綜合偽譜法和BF算法。列車-軌道動力學仿真的作用為計算不同線路條件和列車速度下的車輪、鋼軌磨損量,形成輪軌磨耗數(shù)據(jù)庫,該過程將于第3.2節(jié)詳細說明。BF算法用于搜尋符合約束條件的坡長方案,并依據(jù)第2節(jié)的分階段原則將區(qū)間劃分為多個階段。偽譜法則基于輪軌磨耗數(shù)據(jù)庫和階段劃分方案,求解最小化雙向能耗成本和維修成本之和的協(xié)同優(yōu)化方案。

    若已有輪軌磨耗數(shù)據(jù)庫,算法流程如圖4所示。具體步驟為:

    步驟1:參數(shù)初始化。輸入?yún)^(qū)間線路平面參數(shù)、車站高程、列車參數(shù)、線路限速、避讓區(qū)域位置、輪軌磨耗數(shù)據(jù)庫等,計算區(qū)間最大坡道數(shù)量Kmax,令K=Kmax:

    Kmax=(L-Ls)/Lt」。(38)

    步驟2:階段劃分。運用BF算法搜索滿足式(18)和式(20)的全部K個坡道組合方案,在此基礎上根據(jù)第2節(jié)的階段劃分原則將區(qū)間劃分為R個階段。

    步驟3:偽譜法求解?;陔A段劃分方案和輪軌磨耗數(shù)據(jù)庫,采用偽譜法求解K個坡道的最優(yōu)坡度和定時列車控制方案,對比各優(yōu)化方案成本,獲得成本最小的K個坡道的最優(yōu)縱斷面方案和列車控制方案。偽譜法原理為將問題離散變換為非線性規(guī)劃問題后求解,該過程基于MATLAB開發(fā)的GPOPS II完成,數(shù)學原理在Wang等[20]的文章中詳細闡述。

    步驟4:優(yōu)化方案對比。將目前最優(yōu)協(xié)同優(yōu)化方案與步驟3輸出的K個坡道的協(xié)同優(yōu)化方案進行對比,令能耗成本和維修成本之和更小的方案為新的最優(yōu)協(xié)同優(yōu)化方案,至步驟5。

    步驟5:循環(huán)判斷。判斷K≤3是否滿足,若是,則結(jié)束算法,輸出最優(yōu)協(xié)同優(yōu)化方案,否則令K=K-1,轉(zhuǎn)至步驟2。

    3.2" 輪軌磨耗數(shù)據(jù)庫構(gòu)建方法

    列車行駛過程中,線路條件和速度都可能發(fā)生變化,導致不同位置處的輪軌磨耗程度不同。以上因素對輪軌磨損的共同影響難以用數(shù)學公式表達,無法使用數(shù)值方法計算不同優(yōu)化方案下的輪軌磨損量。因此,本文構(gòu)建輪軌磨耗數(shù)據(jù)庫,在優(yōu)化過程中調(diào)用數(shù)據(jù)庫計算輪軌磨損量,估算輪軌維修成本。

    根據(jù)第2節(jié)維修成本分析,輪軌磨耗數(shù)據(jù)庫數(shù)據(jù)包括單位步長車輪踏面直徑減少量、鋼軌磨損深度減少量,數(shù)據(jù)基于SIMPACK軟件仿真計算獲得。構(gòu)建數(shù)據(jù)庫的流程如圖5所示。

    首先輸入車輛和軌道系統(tǒng)固定參數(shù)、線路坡度和曲線半徑以及列車速度,建立列車-軌道仿真模型,獲得輪軌接觸參數(shù);接著用磨耗模型計算單位步長車輪輪緣厚度減少量及相應踏面直徑減少量、鋼軌磨損深度減少量;不斷改變線路坡度、曲線半徑和列車速度的組合,重復以上步驟,直至完成對所有組合條件的仿真計算,形成輪軌磨耗數(shù)據(jù)庫。

    4" 案例分析

    本文應用廣州地鐵3個區(qū)間進行案例分析,基于SIMPACK計算并比較不同速度、坡度下的車輪和鋼軌磨損量差異,總結(jié)輪軌磨損規(guī)律,形成輪軌磨耗數(shù)據(jù)庫?;谳嗆壞ズ臄?shù)據(jù)庫,對比協(xié)同優(yōu)化方案和基于實際縱斷面的列車控制優(yōu)化方案的成本大小,驗證協(xié)同優(yōu)化方法的優(yōu)化效果。算法框架在MATLAB R2018b中編程實現(xiàn),偽譜法求解過程基于GPOPS II完成。

    4.1" 基礎數(shù)據(jù)

    已知的線路參數(shù)和運行時間如表1所示,模型中主要參數(shù)的值見表2,縱斷面設計約束相關參數(shù)根據(jù)《地鐵設計規(guī)范》和實際情況取值。列車特性曲線如圖6所示。構(gòu)建列車-軌道動力學模型時,車輪外形為LM型,鋼軌為UIC60鋼軌。

    4.2" 輪軌磨損變化規(guī)律

    為計算單位步長車輪輪緣厚度減少量、鋼軌磨損深度減少量,形成輪軌磨耗數(shù)據(jù)庫,在仿真線路上建立列車-軌道動力學模型。仿真線路依次由直線、緩和曲線、圓曲線、緩和曲線和直線組成,直線、圓曲線長度均為100 m,緩和曲線長度按曲線半徑取值。根據(jù)表2中的曲線半徑,仿真線路圓曲線半徑取值為700、800和2 000 m。根據(jù)《地鐵設計規(guī)范》和限速,坡度取值范圍為[-30‰,-2‰]和[2‰,30‰],列車速度在10~75 km/h間取值。

    根據(jù)仿真計算結(jié)果,當曲線半徑、列車速度一定,分別為700 m和70 km/h時,車輪和鋼軌磨損量隨坡度變化情況如圖7所示。由圖7可知,無論是上坡或下坡,車輪和鋼軌單位步長磨損量都隨著坡度絕對值增大而增加。當其他條件相同時,相比列車在上坡道上運行,列車在下坡道上運行的輪軌磨損量一般更大。

    選取坡度為-25‰的情形,對曲線半徑分別為700和2 000 m,不同列車速度下的輪軌磨損量進行分析,磨損量變化如圖8所示。當曲線半徑不變時,車輪和鋼軌單位步長磨損量隨列車速度增加而增長。當列車速度一定時,曲線半徑越大,車輪和鋼軌單位步長磨損量越小。綜合以上分析,如果不考慮能耗成本和約束條件等,曲線半徑越大、列車速度和坡度絕對值越小,輪軌維修成本將越低。

    4.3" 協(xié)同優(yōu)化模型和列車控制優(yōu)化模型效果對比

    為論證協(xié)同優(yōu)化方法的優(yōu)化效果,基于輪軌磨耗數(shù)據(jù)庫,對表1的3個區(qū)間使用本文模型和列車控制優(yōu)化模型,對比同區(qū)間不同方案下的成本大小。列車控制優(yōu)化模型的基礎為Wang等[27]提出的模型,將約束條件與本文模型保持一致,該模型采用實際縱斷面,僅對給定區(qū)間運行時分條件下的列車控制方案進行優(yōu)化,使得能耗成本和維修成本之和最小。

    協(xié)同優(yōu)化模型求解的縱斷面方案和實際縱斷面方案如圖9所示,圖中為上行方向縱斷面方案,該方案軸對稱即為下行方向縱斷面方案。由于3個區(qū)間不是連續(xù)區(qū)間,起點高程都設為0 m。由圖可知,相比3個區(qū)間的實際縱斷面方案,優(yōu)化縱斷面方案的高程最小值更小,站外區(qū)域出站、進站坡度都更大。

    當列車在兩種縱斷面上運行,研究年度產(chǎn)生的成本如表3所示,成本節(jié)約比率為本文模型優(yōu)化方案相比列車控制優(yōu)化模型求解方案的成本變化情況。在各區(qū)間,本文模型優(yōu)化方案較列車控制優(yōu)化模型求解方案都具有更好的優(yōu)化效果,能耗成本和維修成本之和都更小,平均降低了21%。因此,相比只優(yōu)化定時列車控制方案,同時對縱斷面和列車控制方案進行優(yōu)化更有助于降低運營成本。

    對于各區(qū)間的能耗成本,相比列車控制優(yōu)化模型求解方案,本文模型優(yōu)化方案都具有更好的優(yōu)化效果,平均優(yōu)化率為27%。為說明能耗成本節(jié)約的原因,以第一個區(qū)間的優(yōu)化結(jié)果為例進行分析,如圖10所示。

    結(jié)合圖中能耗變化和列車操縱方案,可知列車能耗增長僅發(fā)生在出站牽引階段,隨后列車采用不耗能的惰行工況運行,直至制動減速進站停車。在列車牽引階段,由式(4),所需能耗取決于牽引力和牽引距離大小。列車在實際線路上運行時,牽引距離和牽引力都不小于優(yōu)化方案,所以耗能更多。對比縱斷面方案可知,優(yōu)化方案站外區(qū)域的下坡坡度更大,當列車處于牽引工況時,利用重力勢能轉(zhuǎn)化為動能,可使用更小的牽引力獲得與在實際縱斷面上運行時相同的加速度。同時,當列車處于惰行工況時,在優(yōu)化縱斷面上能充分利用重力勢能獲得更大運行速度以滿足運行時分要求,則牽引階段的運行距離可以更短。

    對于維修費用,在后兩個區(qū)間上,優(yōu)化方案的維修費用小于列車操縱優(yōu)化模型求解方案的維修費用,但在第一個區(qū)間,優(yōu)化方案的維修費用更高,主要由輪軌維修費用較大導致。由于輪軌磨損變化規(guī)律相似,結(jié)合輪軌磨損變化規(guī)律,對第一和第三個區(qū)間的車輪維修費用進行分析,如圖11、12所示。

    圖11(a)為列車控制方案和縱斷面對比,圖11(b)為線路曲線和車輪維修成

    本對比。由圖11(b)可知,第一個區(qū)間雙向共有4處曲線,從左至右依次記為

    曲線1~4。在曲線1、4上,結(jié)合圖11(a)可知,相比實際縱斷面,優(yōu)化縱斷面坡度絕對值更大,但列車速度更小,兩種縱斷面上車輪磨損量接近。在曲線2、3上,結(jié)合圖11(a)可

    知,優(yōu)化縱斷面坡度絕對值和列車速度都大于實際縱斷面,所以車輪磨損量更大。綜合4個曲線上的車輪磨損量,最終優(yōu)化縱斷面上車輪磨損量大于實際縱斷面上車輪磨損量,相應維修成本也更大。

    圖12(a)為列車控制方案和縱斷面對比,圖12(b)為線路曲線和車輪維修成本對

    比。由圖12(b)可知,第3個區(qū)間雙向共有4處曲線,從左至右依次記為曲線

    5~8。在曲線5、8上,結(jié)合圖12(a)可知,相比實際縱斷面,雖然優(yōu)化縱斷面坡度更大,但列車速度更小,導致車輪磨損量稍小。在曲線6上,結(jié)合圖12(a)可知,雖然優(yōu)化縱斷面上的列車速度更大,但所處位置為上坡,而實際縱斷面相應位置為下坡,由于下坡時車輪磨損量一般更大,此時列車在優(yōu)化縱斷面上的車輪磨損量更小。在曲線7上,結(jié)合圖12(a)可知,相比實際縱斷面,優(yōu)化縱斷面為下坡且列車速度更大,所以車輪磨損量更大。綜合4個曲線上的車輪磨損量,最終優(yōu)化縱斷面上車輪磨損量小于實際縱斷面上車輪磨損量,相應維修成本也更小。結(jié)合第1和第3個區(qū)間的維修成本分析可知,在其他條件相同時,車輪維修成本受坡度、列車速度共同影響,車輪磨損量及相應維修成本需結(jié)合具體場景進行分析。

    5" 結(jié)語

    相比僅優(yōu)化列車控制,在設計縱斷面時優(yōu)化列車控制將有助于節(jié)約運營成本。本文以最小化能耗成本和維修成本之和為目標,提出了考慮定時列車控制對縱斷面進行優(yōu)化的模型,并基于列車-軌道動力學仿真,設計了綜合偽譜法和蠻力搜索的算法求解。通過廣州地鐵三個區(qū)間的實例證明,對比僅優(yōu)化定時列車控制方案,本文模型優(yōu)化方案對降低能耗成本和維修成本具有更好的優(yōu)化效果,各區(qū)間平均優(yōu)化率為21%。本文所用偽譜法采用通用求解器實現(xiàn),考慮協(xié)同優(yōu)化問題求解時間較長,未來可結(jié)合問題特點開發(fā)專用求解器,提高求解效率。

    參考文獻:

    [1]中國城市軌道交通協(xié)會. 城市軌道交通2022年度統(tǒng)計和分析報告[EB/OL]. [2024-04-01].https://www.camet.org.cn/xytj/tjxx/11944.shtml.

    [2]GONZLEZ-GIL A, PALACIN R, BATTY P, et al. A systems approach to reduce urban rail energy consumption[J]. Energy Conversion and Management, 2014, 80: 509-524. DOI: 10.1016/j.enconman.2014.01.060.

    [3]DOUGLAS H, ROBERTS C, HILLMANSEN S, et al. An assessment of available measures to reduce traction energy use in railway networks[J]. Energy Conversion and Management, 2015, 106: 1149-1165. DOI: 10.1016/j.enconman.2015.10.053.

    [4]HOANG H, POLIS M, HAURIE A. Reducing energy consumption through trajectory optimization for a metro network[J]. IEEE Transactions on Automatic Control, 1975, 20(5): 590-595. DOI: 10.1109/tac.1975.1101058.

    [5]JHA M K, SCHONFELD P, SAMANTA S. Optimizing rail transit routes with genetic algorithms and geographic information system[J]. Journal of Urban Planning and Development, 2007, 133(3): 161-171. DOI: 10.1061/(asce)0733-9488(2007)133: 3(161).

    [6]KIM M E, SCHONFELD P, KIM E. Simulation-based rail transit optimization model[J]. Transportation Research Record: Journal of the Transportation Research Board, 2013, 2374(1): 143-153. DOI: 10.3141/2374-17.

    [7]柏赟, 白驍, 孫元廣, 等. 地鐵線路區(qū)間縱斷面節(jié)能設計優(yōu)化模型[J]. 鐵道學報, 2020, 42(9): 10-16. DOI: 10.3969/j.issn.1001-8360.2020.09.002.

    [8]彭磊, 樊蔥, 汪茜, 等. 地鐵線路縱斷面優(yōu)化系統(tǒng)設計與實現(xiàn)[J]. 山東科學, 2022, 35(4): 136-143. DOI: 10.3976/j.issn.1002-4026.2022.04.017.

    [9]徐楊. 考慮建設成本和運營能耗的地鐵線路縱斷面設計優(yōu)化[D]. 北京: 北京交通大學, 2020: 27-54.

    [10]WU C X, LU S F, XUE F, et al. Earth potential as the energy storage in rail transit system-on a vertical alignment optimization problem[C]//2018 21st International Conference on Intelligent Transportation Systems (ITSC). Maui, HI. IEEE, 2018: 2729-2734 DOI: 10.1109/itsc.2018.8569925.

    [11]DUARTE M A, SOTOMAYOR P X. Minimum energy trajectories for subway systems[J]. Optimal Control Applications and Methods, 1999, 20(6): 283-296. DOI: 10.1002/(sici)1099-1514(199911/12)20: 6lt;283: aid-oca661gt;3.0.co;2-k.

    [12]湯超. 平縱斷面線型及軌道類型對地鐵軌道動力特性和鋼軌磨耗影響研究[D]. 北京: 北京交通大學, 2012: 31-64.

    [13]肖緋雄, 陳旭. 惡劣工況下貨車高摩合成閘瓦磨損分析[J]. 北京交通大學學報, 2014, 38(4): 20-25. DOI: 10.11860/j.issn.1673-0291-2014.04.004.

    [14]LAI X R, SCHONFELD P. Optimization of rail transit alignments considering vehicle dynamics[J]. Transportation Research Record: Journal of the Transportation Research Board, 2012, 2275(1): 77-87. DOI: 10.3141/2275-09.

    [15]KANG M W, JHA M K, BUDDHARAJU R. Rail transit route optimization model for rail infrastructure planning and design: Case study of saint Andrews, Scotland[J]. Journal of Transportation Engineering, 2014, 140(1): 1-11. DOI: 10.1061/(asce)te.1943-5436.0000445.

    [16]中華人民共和國住房和城鄉(xiāng)建設部,中華人民共和國國家質(zhì)量監(jiān)督檢驗檢疫總局.GB 50157—2013地鐵設計規(guī)范[S]. 北京:中國建筑工業(yè)出版社,2014.

    [17]岳小麗, 劉柳林. 隧道盾構(gòu)法施工的成本分析和降低措施[J]. 工業(yè)建筑, 2009, 39(S1): 755-758.

    [18]廖濤. 動車組制動閘片壽命預測關鍵技術的研究與實現(xiàn)[D]. 北京: 北京交通大學, 2015: 9-15.

    [19]耿傳智, 余慶. 地鐵軌道結(jié)構(gòu)減振性能的仿真分析[J]. 同濟大學學報(自然科學版), 2011, 39(1): 85-89. DOI: 10.3969/j.issn.0253-374x.2011.01.016.

    [20]WANG P L, GOVERDE R M P. Multiple-phase train trajectory optimization with signalling and operational constraints[J]. Transportation Research Part C: Emerging Technologies, 2016, 69: 255-275. DOI: 10.1016/j.trc.2016.06.008.

    [21]王平, 郭強, 李晨鐘, 等. 蠕滑曲線對地鐵小半徑曲線輪軌接觸特性的影響[J]. 中南大學學報(自然科學版), 2020, 51(4): 1145-1153. DOI: 10.11817/j.issn.1672-7207.2020.04.029.

    [22]陳旭. 貨車閘瓦不同工況仿真分析及磨損研究[D]. 成都: 西南交通大學, 2015: 18-20.

    [23]楊緒杰. 基于幾何與動力學參數(shù)的山區(qū)線路地鐵車輛輪軌磨耗控制研究[D]. 重慶: 重慶交通大學, 2014: 1-14.

    [24]于春廣, 陶功權. 地鐵車輪磨耗測試及數(shù)值仿真[J]. 工程力學, 2016, 33(1): 201-208. DOI: 10.6052/j.issn.1000-4750.2014.06.0484.

    [25]沈建文, 高亮, 戴春陽, 等. 減振軌道形式對地鐵曲線上鋼軌磨耗影響的仿真[J]. 北京交通大學學報, 2011, 35(4): 38-43. DOI: 10.3969/j.issn.1673-0291.2011.04.008.

    [26]中華人民共和國鐵道部.TB 2097—1989鋼軌允許磨耗限度[S]. 北京:中華人民共和國鐵道部,1990.

    [27]WANG Y H, DE SCHUTTER B, VAN DEN BOOM T J J, et al. Optimal trajectory planning for trains:a pseudospectral method and a mixed integer linear programming approach[J]. Transportation Research Part C: Emerging Technologies, 2013, 29: 97-114. DOI: 10.1016/j.trc.2013.01.007.

    猜你喜歡
    縱斷面磨損量輪軌
    基于車體加速度的超大跨度橋上線路縱斷面優(yōu)化方法
    磁懸浮列車制動閘片磨損行為及其剩余壽命預測研究
    地鐵線路縱斷面優(yōu)化系統(tǒng)設計與實現(xiàn)
    山東科學(2022年4期)2022-07-31 15:48:54
    中低速磁浮道岔與輪軌道岔的差異
    基于輪廓提取的刀具磨損量檢測研究
    曳引輪不均勻磨損量的檢測
    關節(jié)軸承壽命試驗機在線磨損量檢測綜合誤差建模
    中低速磁浮與輪軌交通信號系統(tǒng)的差異
    三維可視化技術在道路縱斷面設計及土方量算中的應用
    地礦測繪(2015年3期)2015-12-22 06:27:30
    非線性穩(wěn)態(tài)曲線通過時輪軌滾動接觸的數(shù)值求解方法
    桃色一区二区三区在线观看| 亚洲自偷自拍图片 自拍| 大型av网站在线播放| 国产男靠女视频免费网站| 12—13女人毛片做爰片一| 777久久人妻少妇嫩草av网站| 国产成人精品久久二区二区免费| 久久伊人香网站| 国产一区二区在线av高清观看| 国产一卡二卡三卡精品| 99国产极品粉嫩在线观看| 神马国产精品三级电影在线观看 | 99在线视频只有这里精品首页| 97人妻精品一区二区三区麻豆 | 亚洲国产欧美日韩在线播放| 正在播放国产对白刺激| 精品国产超薄肉色丝袜足j| 变态另类成人亚洲欧美熟女| 法律面前人人平等表现在哪些方面| 99久久综合精品五月天人人| 亚洲第一电影网av| 国内久久婷婷六月综合欲色啪| 日本一本二区三区精品| 国内揄拍国产精品人妻在线 | 亚洲欧美一区二区三区黑人| xxxwww97欧美| 一本综合久久免费| ponron亚洲| 亚洲片人在线观看| 一级片免费观看大全| 在线观看免费午夜福利视频| 日韩有码中文字幕| 在线观看一区二区三区| 精品人妻1区二区| 国产精品影院久久| 俺也久久电影网| av免费在线观看网站| 国产黄片美女视频| 人人澡人人妻人| 色尼玛亚洲综合影院| 91老司机精品| 久久伊人香网站| 国产又爽黄色视频| 一区二区三区激情视频| 精品国产亚洲在线| 欧美 亚洲 国产 日韩一| 日本撒尿小便嘘嘘汇集6| 麻豆成人av在线观看| 美女 人体艺术 gogo| 国产成人影院久久av| 亚洲av电影在线进入| 18禁美女被吸乳视频| 午夜福利在线观看吧| 亚洲色图 男人天堂 中文字幕| 国产精品美女特级片免费视频播放器 | 男女下面进入的视频免费午夜 | 天天躁狠狠躁夜夜躁狠狠躁| 热re99久久国产66热| 色哟哟哟哟哟哟| 美女免费视频网站| 中文字幕人妻丝袜一区二区| 亚洲中文av在线| 国产伦人伦偷精品视频| 妹子高潮喷水视频| 午夜免费成人在线视频| 国产久久久一区二区三区| 国产亚洲欧美98| 精品国产亚洲在线| 久久精品aⅴ一区二区三区四区| 亚洲熟妇熟女久久| 精品人妻1区二区| 国产一区二区在线av高清观看| 十八禁人妻一区二区| 国产97色在线日韩免费| 国产激情欧美一区二区| 国产精品亚洲一级av第二区| 天天躁狠狠躁夜夜躁狠狠躁| 成人亚洲精品一区在线观看| 色播亚洲综合网| 国产精品久久视频播放| 香蕉国产在线看| 香蕉丝袜av| 午夜久久久在线观看| 最近在线观看免费完整版| 99国产综合亚洲精品| 久久中文看片网| 亚洲中文av在线| 香蕉av资源在线| 国产人伦9x9x在线观看| 国产成人av激情在线播放| 波多野结衣巨乳人妻| 久久 成人 亚洲| 亚洲精品av麻豆狂野| 91大片在线观看| 好男人在线观看高清免费视频 | 国产精品 欧美亚洲| 久久久久久亚洲精品国产蜜桃av| 午夜成年电影在线免费观看| 国产精品电影一区二区三区| 9191精品国产免费久久| 变态另类成人亚洲欧美熟女| 日韩 欧美 亚洲 中文字幕| 狠狠狠狠99中文字幕| 少妇的丰满在线观看| 在线观看舔阴道视频| 国产av在哪里看| 国产精品日韩av在线免费观看| 日韩中文字幕欧美一区二区| 每晚都被弄得嗷嗷叫到高潮| 欧美激情高清一区二区三区| 亚洲国产精品久久男人天堂| 黄网站色视频无遮挡免费观看| 日本精品一区二区三区蜜桃| 国产精华一区二区三区| 亚洲欧美激情综合另类| 久久国产乱子伦精品免费另类| 999精品在线视频| 黄网站色视频无遮挡免费观看| 岛国在线观看网站| 亚洲中文日韩欧美视频| 一区二区三区国产精品乱码| 国产亚洲精品综合一区在线观看 | 精品欧美一区二区三区在线| 亚洲在线自拍视频| 色老头精品视频在线观看| 日韩大码丰满熟妇| 欧美一级毛片孕妇| 黄片大片在线免费观看| 又大又爽又粗| 国产色视频综合| 欧美zozozo另类| 搞女人的毛片| 午夜日韩欧美国产| 日日摸夜夜添夜夜添小说| 50天的宝宝边吃奶边哭怎么回事| 好看av亚洲va欧美ⅴa在| 精品国产亚洲在线| 国产精品久久久av美女十八| 国产黄色小视频在线观看| 国产日本99.免费观看| 久久久久免费精品人妻一区二区 | 国产一区二区激情短视频| 国产麻豆成人av免费视频| 美国免费a级毛片| 搡老岳熟女国产| 亚洲国产欧美一区二区综合| 一级毛片精品| 午夜免费成人在线视频| 日韩高清综合在线| 人人澡人人妻人| 伊人久久大香线蕉亚洲五| 亚洲性夜色夜夜综合| 中出人妻视频一区二区| 欧美亚洲日本最大视频资源| 日韩欧美免费精品| 最新美女视频免费是黄的| 一区二区三区高清视频在线| 国产精品精品国产色婷婷| 久久久久国产一级毛片高清牌| 午夜激情av网站| 免费电影在线观看免费观看| 国产成人av教育| 一本大道久久a久久精品| 日本 av在线| 1024视频免费在线观看| 国产成人精品久久二区二区免费| 色综合亚洲欧美另类图片| 国产熟女午夜一区二区三区| 久久精品亚洲精品国产色婷小说| 国产成年人精品一区二区| 久久精品91无色码中文字幕| 精品欧美一区二区三区在线| 国产高清视频在线播放一区| 99久久无色码亚洲精品果冻| 国产1区2区3区精品| 亚洲国产毛片av蜜桃av| 久久香蕉激情| av免费在线观看网站| 国产av不卡久久| 在线看三级毛片| 成年版毛片免费区| 琪琪午夜伦伦电影理论片6080| 性欧美人与动物交配| 欧美精品亚洲一区二区| 露出奶头的视频| 欧美日韩亚洲国产一区二区在线观看| 丝袜在线中文字幕| 欧美黑人巨大hd| 久久青草综合色| 国产亚洲精品第一综合不卡| 久久午夜亚洲精品久久| 搡老熟女国产l中国老女人| 黄色女人牲交| 欧美色视频一区免费| 韩国av一区二区三区四区| 丁香欧美五月| 国产成人影院久久av| 丁香六月欧美| 丰满的人妻完整版| 在线av久久热| 少妇 在线观看| 亚洲中文av在线| 成人亚洲精品av一区二区| 黄网站色视频无遮挡免费观看| 精品久久久久久久毛片微露脸| 99热这里只有精品一区 | 天堂影院成人在线观看| 精品久久久久久,| 国产成人影院久久av| 午夜日韩欧美国产| 久久天堂一区二区三区四区| 国产精品亚洲av一区麻豆| 久久久久国产精品人妻aⅴ院| 在线观看免费视频日本深夜| 精品乱码久久久久久99久播| 男人舔奶头视频| 亚洲激情在线av| 色播亚洲综合网| 99热这里只有精品一区 | 国产私拍福利视频在线观看| 欧美zozozo另类| 精品欧美一区二区三区在线| 中文字幕精品免费在线观看视频| 午夜免费成人在线视频| 淫秽高清视频在线观看| bbb黄色大片| 啦啦啦韩国在线观看视频| 亚洲精品国产区一区二| 中文字幕最新亚洲高清| 青草久久国产| 亚洲成人国产一区在线观看| 日韩欧美三级三区| 91成年电影在线观看| 91国产中文字幕| 亚洲天堂国产精品一区在线| 国产精品1区2区在线观看.| 国产主播在线观看一区二区| 99久久无色码亚洲精品果冻| 可以免费在线观看a视频的电影网站| 丁香六月欧美| www.自偷自拍.com| 欧美激情极品国产一区二区三区| 亚洲成av片中文字幕在线观看| 欧美一级毛片孕妇| 亚洲五月色婷婷综合| 一进一出好大好爽视频| 国产精华一区二区三区| 淫妇啪啪啪对白视频| 好男人电影高清在线观看| 叶爱在线成人免费视频播放| 精品欧美一区二区三区在线| 欧美日韩精品网址| 美女午夜性视频免费| or卡值多少钱| 十分钟在线观看高清视频www| 亚洲va日本ⅴa欧美va伊人久久| 亚洲精品一区av在线观看| 嫩草影院精品99| 国产成人一区二区三区免费视频网站| 好看av亚洲va欧美ⅴa在| 久久99热这里只有精品18| 性色av乱码一区二区三区2| 两个人看的免费小视频| 国产高清激情床上av| 久久这里只有精品19| 在线观看免费视频日本深夜| 国产亚洲精品久久久久5区| 亚洲在线自拍视频| 久久天堂一区二区三区四区| 亚洲男人的天堂狠狠| 91字幕亚洲| 久久精品国产亚洲av高清一级| 久久久久国产一级毛片高清牌| 麻豆国产av国片精品| 精品久久久久久久毛片微露脸| 可以在线观看的亚洲视频| 人人澡人人妻人| 成人国产综合亚洲| 国产精品二区激情视频| 欧美黄色淫秽网站| 国产免费男女视频| 久久狼人影院| 久久天躁狠狠躁夜夜2o2o| 国产成人欧美在线观看| 757午夜福利合集在线观看| 国产又爽黄色视频| 一级片免费观看大全| 淫妇啪啪啪对白视频| www.自偷自拍.com| 两个人视频免费观看高清| 欧美国产日韩亚洲一区| 熟女电影av网| 成人免费观看视频高清| 婷婷丁香在线五月| 精品少妇一区二区三区视频日本电影| 宅男免费午夜| 69av精品久久久久久| 欧美一级毛片孕妇| 中文亚洲av片在线观看爽| 亚洲成av片中文字幕在线观看| 在线观看日韩欧美| 校园春色视频在线观看| 一级a爱视频在线免费观看| 无遮挡黄片免费观看| 国产精品香港三级国产av潘金莲| 午夜福利视频1000在线观看| 精品卡一卡二卡四卡免费| 搡老岳熟女国产| 18禁裸乳无遮挡免费网站照片 | 天天一区二区日本电影三级| 波多野结衣av一区二区av| 亚洲国产欧美一区二区综合| 九色国产91popny在线| 神马国产精品三级电影在线观看 | 国产一区二区在线av高清观看| 大香蕉久久成人网| 老汉色∧v一级毛片| 日韩欧美一区二区三区在线观看| 两个人视频免费观看高清| 久久草成人影院| 亚洲精品国产一区二区精华液| 成人午夜高清在线视频 | 午夜久久久久精精品| 一边摸一边做爽爽视频免费| 久久精品aⅴ一区二区三区四区| 黄色视频,在线免费观看| 免费在线观看亚洲国产| 搞女人的毛片| 免费av毛片视频| 久久香蕉激情| 国产熟女xx| 亚洲专区中文字幕在线| 国产蜜桃级精品一区二区三区| 国产伦一二天堂av在线观看| 国产精品影院久久| 极品教师在线免费播放| 亚洲精品国产一区二区精华液| 亚洲人成网站在线播放欧美日韩| 国产黄a三级三级三级人| 在线av久久热| 亚洲中文字幕日韩| 国产精品日韩av在线免费观看| 国产精品一区二区免费欧美| 黄片小视频在线播放| 亚洲第一电影网av| 免费av毛片视频| www.精华液| 国产真人三级小视频在线观看| 国产日本99.免费观看| 大型黄色视频在线免费观看| 国产极品粉嫩免费观看在线| 一进一出抽搐gif免费好疼| 午夜福利视频1000在线观看| 欧美三级亚洲精品| 久久久久九九精品影院| 国内揄拍国产精品人妻在线 | 人妻丰满熟妇av一区二区三区| 在线永久观看黄色视频| 黄色丝袜av网址大全| 亚洲精品久久成人aⅴ小说| 99久久综合精品五月天人人| 美女午夜性视频免费| 亚洲第一欧美日韩一区二区三区| 18禁黄网站禁片免费观看直播| 欧美色欧美亚洲另类二区| 欧美不卡视频在线免费观看 | 免费在线观看亚洲国产| 久久久久久亚洲精品国产蜜桃av| 亚洲自偷自拍图片 自拍| 村上凉子中文字幕在线| 国产精品久久久久久人妻精品电影| xxx96com| 午夜免费成人在线视频| 麻豆av在线久日| 亚洲av片天天在线观看| 美女午夜性视频免费| 美女 人体艺术 gogo| 国产爱豆传媒在线观看 | 观看免费一级毛片| 99热只有精品国产| 亚洲在线自拍视频| 变态另类丝袜制服| 成人精品一区二区免费| 国产精品精品国产色婷婷| netflix在线观看网站| 视频区欧美日本亚洲| 一级毛片精品| 日韩视频一区二区在线观看| 中文字幕精品免费在线观看视频| 国产成人精品久久二区二区91| 久久久久国产精品人妻aⅴ院| 日本a在线网址| 亚洲专区字幕在线| 99国产精品一区二区蜜桃av| 变态另类成人亚洲欧美熟女| 午夜福利视频1000在线观看| 12—13女人毛片做爰片一| 亚洲一卡2卡3卡4卡5卡精品中文| 精品电影一区二区在线| 日韩一卡2卡3卡4卡2021年| 午夜福利在线在线| 久久精品国产99精品国产亚洲性色| 日韩成人在线观看一区二区三区| 国产亚洲欧美精品永久| 婷婷六月久久综合丁香| 久久人妻av系列| 91麻豆精品激情在线观看国产| 国产又黄又爽又无遮挡在线| 老汉色av国产亚洲站长工具| 久久久久久久久中文| 90打野战视频偷拍视频| 欧美日本视频| 久久人妻福利社区极品人妻图片| 女人爽到高潮嗷嗷叫在线视频| 久久久久久国产a免费观看| 波多野结衣巨乳人妻| 亚洲狠狠婷婷综合久久图片| 国产亚洲精品av在线| 极品教师在线免费播放| 女性被躁到高潮视频| 夜夜爽天天搞| 欧美三级亚洲精品| 亚洲成人久久爱视频| 久久久久亚洲av毛片大全| 十分钟在线观看高清视频www| 亚洲国产精品999在线| 色综合欧美亚洲国产小说| 国产精品香港三级国产av潘金莲| 午夜久久久久精精品| 99精品欧美一区二区三区四区| 啦啦啦观看免费观看视频高清| 哪里可以看免费的av片| 精品国产国语对白av| 午夜亚洲福利在线播放| av免费在线观看网站| 男人的好看免费观看在线视频 | 免费高清视频大片| 午夜日韩欧美国产| 白带黄色成豆腐渣| 巨乳人妻的诱惑在线观看| 成在线人永久免费视频| 黄色视频,在线免费观看| 日本三级黄在线观看| 日韩免费av在线播放| 极品教师在线免费播放| 高潮久久久久久久久久久不卡| 日韩精品免费视频一区二区三区| 免费在线观看亚洲国产| 黄色 视频免费看| 亚洲精品久久国产高清桃花| 欧美国产精品va在线观看不卡| 色综合站精品国产| 久久人妻av系列| 在线观看一区二区三区| 国产视频内射| 久久久久久九九精品二区国产 | 免费女性裸体啪啪无遮挡网站| 一级毛片精品| 男女床上黄色一级片免费看| √禁漫天堂资源中文www| 国产成人系列免费观看| bbb黄色大片| 不卡一级毛片| 亚洲精品久久国产高清桃花| 麻豆成人av在线观看| 亚洲aⅴ乱码一区二区在线播放 | 97人妻精品一区二区三区麻豆 | 亚洲自拍偷在线| 色综合欧美亚洲国产小说| 国产欧美日韩一区二区三| 啦啦啦观看免费观看视频高清| 91成年电影在线观看| 日韩国内少妇激情av| 18美女黄网站色大片免费观看| 欧美亚洲日本最大视频资源| 亚洲成av人片免费观看| 黄色成人免费大全| 1024香蕉在线观看| 日本精品一区二区三区蜜桃| 可以在线观看毛片的网站| 久热爱精品视频在线9| 亚洲一码二码三码区别大吗| 香蕉久久夜色| 亚洲国产高清在线一区二区三 | 亚洲国产看品久久| 两性夫妻黄色片| 精品不卡国产一区二区三区| 久久久久久久精品吃奶| 国产成人精品久久二区二区免费| 国产又色又爽无遮挡免费看| 欧美激情 高清一区二区三区| 久久香蕉激情| 精品不卡国产一区二区三区| 久久精品国产亚洲av高清一级| 国产成人精品久久二区二区免费| www.999成人在线观看| 在线播放国产精品三级| 午夜影院日韩av| av福利片在线| 国产一级毛片七仙女欲春2 | 91老司机精品| 免费在线观看完整版高清| 日本三级黄在线观看| 美女大奶头视频| 老司机午夜福利在线观看视频| 成人三级黄色视频| 成人亚洲精品一区在线观看| 国产精品久久久久久亚洲av鲁大| 久久久国产成人免费| 亚洲精品中文字幕在线视频| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日韩精品网址| 国产精品香港三级国产av潘金莲| 国产精品亚洲av一区麻豆| 亚洲av成人不卡在线观看播放网| 国产成人av教育| 国产精品一区二区精品视频观看| 99国产精品一区二区蜜桃av| 国产乱人伦免费视频| 久久精品91无色码中文字幕| 成人特级黄色片久久久久久久| 亚洲aⅴ乱码一区二区在线播放 | 精品欧美一区二区三区在线| 人人妻,人人澡人人爽秒播| 久久久久久久精品吃奶| 岛国在线观看网站| 黄色视频,在线免费观看| √禁漫天堂资源中文www| 午夜免费观看网址| 国产精品久久电影中文字幕| 成人国产一区最新在线观看| 国产真实乱freesex| 国产午夜福利久久久久久| 俺也久久电影网| 亚洲一区二区三区色噜噜| 久久精品成人免费网站| 天天一区二区日本电影三级| 一个人免费在线观看的高清视频| 国产激情欧美一区二区| 国产黄a三级三级三级人| 欧美最黄视频在线播放免费| 国产精品二区激情视频| 我的亚洲天堂| 哪里可以看免费的av片| 久久 成人 亚洲| 亚洲中文av在线| 热re99久久国产66热| 国产精品精品国产色婷婷| 欧美激情久久久久久爽电影| 这个男人来自地球电影免费观看| 亚洲熟女毛片儿| 国产亚洲av嫩草精品影院| 久久精品91蜜桃| ponron亚洲| 免费女性裸体啪啪无遮挡网站| 久久精品国产综合久久久| 天堂√8在线中文| 欧美一级a爱片免费观看看 | 精华霜和精华液先用哪个| 久久精品国产99精品国产亚洲性色| 97人妻精品一区二区三区麻豆 | 婷婷精品国产亚洲av| 午夜福利在线观看吧| 熟女电影av网| 桃红色精品国产亚洲av| 一区二区三区激情视频| 99久久无色码亚洲精品果冻| 亚洲精品一卡2卡三卡4卡5卡| 久久久久九九精品影院| 午夜两性在线视频| 制服诱惑二区| 男人舔女人的私密视频| 老鸭窝网址在线观看| 欧美最黄视频在线播放免费| 日韩欧美一区视频在线观看| 亚洲国产中文字幕在线视频| av欧美777| 伦理电影免费视频| 天堂动漫精品| 麻豆一二三区av精品| 久久人妻福利社区极品人妻图片| 国产精品乱码一区二三区的特点| av超薄肉色丝袜交足视频| 日本黄色视频三级网站网址| 亚洲一卡2卡3卡4卡5卡精品中文| 97碰自拍视频| 可以免费在线观看a视频的电影网站| 一a级毛片在线观看| 午夜免费鲁丝| 一个人观看的视频www高清免费观看 | 亚洲精品久久国产高清桃花| 亚洲自拍偷在线| 脱女人内裤的视频| av超薄肉色丝袜交足视频| 最好的美女福利视频网| 淫妇啪啪啪对白视频| 人成视频在线观看免费观看| 每晚都被弄得嗷嗷叫到高潮| 日本在线视频免费播放| 久久久精品欧美日韩精品| 亚洲 欧美 日韩 在线 免费| 精品日产1卡2卡| 成人免费观看视频高清| 熟女电影av网| 欧美成人免费av一区二区三区| 视频区欧美日本亚洲| 在线观看免费视频日本深夜| 可以免费在线观看a视频的电影网站| 国产一区二区在线av高清观看| 免费人成视频x8x8入口观看| 亚洲中文字幕一区二区三区有码在线看 | 亚洲男人天堂网一区| 亚洲最大成人中文| 热99re8久久精品国产| 久久久久久久久久黄片|