亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        高速動(dòng)車組線性渦流制動(dòng)系統(tǒng)特性仿真研究

        2023-02-12 02:58:38楊利強(qiáng)徐凱歌劉賽賽
        城市軌道交通研究 2023年1期
        關(guān)鍵詞:匝數(shù)氣隙渦流

        楊利強(qiáng) 徐凱歌 劉賽賽

        (1.徐州地鐵集團(tuán)有限公司,221018,徐州; 2.徐州地鐵運(yùn)營有限公司,221018,徐州;3.南京中車浦鎮(zhèn)海泰制動(dòng)設(shè)備有限公司,211899,南京∥第一作者,正高級(jí)工程師)

        線性渦流制動(dòng)是一種不依靠車輪與鋼軌之間的摩擦進(jìn)行制動(dòng)的方式,其在高速下可提供穩(wěn)定、可靠的制動(dòng)力,已逐漸成為高速列車上一種新型的制動(dòng)方式[1-2]。但線性渦流制動(dòng)系統(tǒng)受列車速度、氣隙等因素影響,想要獲得安全可靠的制動(dòng)力就要為線性渦流制動(dòng)系統(tǒng)選取合適的參數(shù)范圍。

        關(guān)于線性渦流制動(dòng)的研究,國內(nèi)外學(xué)者已進(jìn)行過大量的試驗(yàn)與分析。法國研發(fā)的AGV(高速動(dòng)車組)列車運(yùn)用的是線性渦流制動(dòng)系統(tǒng)進(jìn)行制動(dòng),當(dāng)列車運(yùn)行速度在200 km/h以上時(shí),每列列車的瞬時(shí)制動(dòng)力可達(dá)到20 kN[3-4]。德國從1985年開始在ICE-V(試驗(yàn)型城際列車特快)試驗(yàn)車上進(jìn)行試驗(yàn)。結(jié)合ICE-V試驗(yàn)車的問題及法國AGV列車的經(jīng)驗(yàn),德國鐵路公司在ICE-3(德國高速動(dòng)車組ICE的第三代)試驗(yàn)車上對(duì)改良后的線性渦流制動(dòng)系統(tǒng)進(jìn)行了試驗(yàn)并取得了良好的效果。其研究結(jié)果表明:ICE-3試驗(yàn)車使用線性渦流制動(dòng)系統(tǒng)進(jìn)行制動(dòng)時(shí),若動(dòng)車速度升高,其制動(dòng)力也能保持較穩(wěn)定的水平,而制動(dòng)功率則跟隨試驗(yàn)車速度的變化而變化(近似的可以看成線性關(guān)系),但線性渦流制動(dòng)熱負(fù)荷卻比摩擦制動(dòng)產(chǎn)生的熱負(fù)荷低很多。這為列車線性渦流制動(dòng)技術(shù)的推廣提供了應(yīng)用基礎(chǔ)[5-6]。1964年,日本建成的新干線是世界上第一條高速鐵路,長度約為550 km,最快速度可達(dá)210 km/h。大約6年后,日本開始在新干線上進(jìn)行一系列線性渦流制動(dòng)試驗(yàn)。1990年,日本又開始研究電磁渦流制動(dòng)技術(shù)和永磁渦流制動(dòng)技術(shù)。日本在其高速列車上運(yùn)用電制動(dòng)與空氣制動(dòng)組成的聯(lián)合制動(dòng)系統(tǒng)進(jìn)行制動(dòng),列車的拖車采用旋轉(zhuǎn)型渦流制動(dòng)系統(tǒng)作為電制動(dòng)系統(tǒng)。旋轉(zhuǎn)型渦流制動(dòng)在高速列車上的成功應(yīng)用表明,在拖車上使用旋轉(zhuǎn)型渦流制動(dòng)與空氣制動(dòng)相結(jié)合的制動(dòng)方式是行之有效的[7]。文獻(xiàn)[8]以“迎流的”有限元法為基礎(chǔ),推導(dǎo)了渦流制動(dòng)力公式,并得到其制動(dòng)力曲線。由該曲線可以看出,隨著列車速度的增加,其制動(dòng)力呈現(xiàn)出先變大后減小的趨勢,并在中間某速度處達(dá)到最大值。文獻(xiàn)[9]對(duì)永磁渦流制動(dòng)器進(jìn)行了分析,求解出其制動(dòng)力公式,并分析了不同參數(shù)變化對(duì)制動(dòng)力的影響。文獻(xiàn)[10]利用德國TR07型磁浮列車的線性渦流制動(dòng)裝置,推導(dǎo)出其制動(dòng)力的關(guān)系式,得到了制動(dòng)力隨列車速度變化的曲線。文獻(xiàn)[11]采用子域法計(jì)算了線性渦流制動(dòng)力與速度、勵(lì)磁電流、氣隙長度、板厚和電導(dǎo)率等參數(shù)之間的解析公式,并通過有限元仿真驗(yàn)證了理論計(jì)算的正確性。由上述研究可知,線性渦流制動(dòng)系統(tǒng)產(chǎn)生的制動(dòng)力與所選取的不同參數(shù)有關(guān),因此研究不同參數(shù)對(duì)線性渦流制動(dòng)系統(tǒng)的影響至關(guān)重要。

        本文通過ANSYS Maxwell仿真軟件對(duì)線性渦流制動(dòng)系統(tǒng)進(jìn)行建模仿真,分析了線性渦流制動(dòng)系統(tǒng)在不同車速、氣隙、勵(lì)磁電流及勵(lì)磁線圈匝數(shù)時(shí)產(chǎn)生的制動(dòng)力大小,為線性渦流制動(dòng)系統(tǒng)的設(shè)計(jì)提供借鑒與參考。

        1 線性渦流制動(dòng)系統(tǒng)的數(shù)學(xué)模型

        高速列車線性渦流制動(dòng)系統(tǒng)是由電磁系統(tǒng)和機(jī)械系統(tǒng)組成的。電磁系統(tǒng)包括牽引變流器、勵(lì)磁電磁鐵等,其中勵(lì)磁電磁鐵由磁軛、電樞線圈和鐵心組成。機(jī)械系統(tǒng)一般由升降風(fēng)缸、套筒等組成。勵(lì)磁電磁鐵與鋼軌平行,其包括8個(gè)磁極,線圈的電流通入方向使磁極按照“N-S-N-S”交替排列。當(dāng)高速列車?yán)镁€性渦流制動(dòng)系統(tǒng)進(jìn)行制動(dòng)時(shí),鋼軌切割電磁鐵產(chǎn)生的磁場使鋼軌中感應(yīng)出渦流并形成渦流磁場,該新生磁場將阻礙原本電磁鐵產(chǎn)生的磁場變化而使氣隙中的磁場畸變產(chǎn)生2個(gè)力的分量,其中一個(gè)橫向電磁力會(huì)阻礙列車運(yùn)動(dòng),從而達(dá)到列車制動(dòng)的目的。

        由于線性渦流制動(dòng)力的計(jì)算涉及到很多復(fù)雜因素,因此,本文根據(jù)線性渦流制動(dòng)系統(tǒng)的工作原理及其在不同工況下的工作過程,采用磁路分析法對(duì)渦流磁場進(jìn)行相應(yīng)簡化,進(jìn)而建立線性渦流制動(dòng)力的數(shù)學(xué)模型。由圖1的線性渦流制動(dòng)磁場示意圖可知,磁極與磁極之間構(gòu)成閉合回路,由于每兩個(gè)磁極之間形成的回路相同且對(duì)稱,因此分析電磁鐵的電磁機(jī)理時(shí),可將多個(gè)磁極等效成一對(duì)磁極進(jìn)行分析。

        圖1 線性渦流制動(dòng)磁場示意圖Fig.1 Diagram of linear eddy current braking magnetic field

        根據(jù)等效原理,可將磁極與鋼軌之間產(chǎn)生的渦流等效成一個(gè)半徑為r、集膚深度為δ的渦流計(jì)算微元。渦流區(qū)域磁通是按照余弦規(guī)律變化的,因此在鋼軌表面產(chǎn)生的感應(yīng)電動(dòng)勢E為:

        (1)

        式中:

        φ——線性渦流區(qū)域磁通;

        B——?dú)庀短幍拇鸥袘?yīng)強(qiáng)度;

        v——列車速度;

        τ——相鄰電磁鐵之間的極距;

        t——變化時(shí)間。

        感應(yīng)渦流在鋼軌表面的集膚深度為δ:

        (2)

        式中:

        ρ——鋼軌的電阻率,Ω/m;

        μ——鋼軌磁導(dǎo)率,H/m。

        由式(1)和式(2)可得,渦流有效值λ為:

        (3)

        式中:

        D——渦流區(qū)域直徑,mm。

        按照等效磁路原理,對(duì)磁路模型各部分磁阻進(jìn)行等效,得到氣隙處的磁感應(yīng)強(qiáng)度為:

        (4)

        式中:

        μ0——空氣磁導(dǎo)率,H/m;

        I0——?jiǎng)?lì)磁電流,A;

        ke——折算系數(shù);

        N——線圈安匝數(shù),匝;

        l0——?dú)庀堕L度,mm。

        從能量守恒的角度出發(fā),可以得到渦流區(qū)域內(nèi)產(chǎn)生的功率有效值P為:

        (5)

        根據(jù)線性渦流制動(dòng)力F=Pv,可得:

        (6)

        從式(6)中可以看出,線性渦流制動(dòng)力與列車速度、氣隙長度、勵(lì)磁電流等諸多因素有關(guān)。此外,渦流密度衰減率β可以表示為:

        (7)

        式中:

        z——深度,mm;

        ω——渦流頻率,Hz;

        σ——鋼軌電導(dǎo)率,S /m。

        由式(7)可知,μ和σ對(duì)電流衰減起到了重要作用。當(dāng)ω不變時(shí),μ或σ越大,則在同一z處衰減得越大[12],因此μ和σ的取值也會(huì)影響制動(dòng)力的大小。

        綜上可知,列車速度、氣隙、勵(lì)磁電流、勵(lì)磁線圈匝數(shù)、鋼軌材料電導(dǎo)率及勵(lì)磁線圈材料等能夠影響線性渦流制動(dòng)特性和電磁特性。

        2 線性渦流制動(dòng)系統(tǒng)仿真模型

        2.1 模型建立

        利用ANSYS Maxwell軟件建立線性渦流制動(dòng)系統(tǒng)仿真模型(以下簡稱“仿真模型”),如圖2所示。仿真模型尺寸如表1所示。

        圖2 線性渦流制動(dòng)系統(tǒng)仿真模型Fig.2 Simulation model of linear eddy current braking system

        表1 線性渦流制動(dòng)系統(tǒng)仿真模型尺寸

        由于鋼軌是鐵磁性物質(zhì),該仿真模型中取其相對(duì)磁導(dǎo)率為200 H/m,電導(dǎo)率為5×106S/m;磁軛與磁心的相對(duì)磁導(dǎo)率取值為1 360 H/m,電導(dǎo)率為0.153 8 S/m;勵(lì)磁線圈選銅材料。

        為保證仿真模型計(jì)算的準(zhǔn)確率,要對(duì)模型進(jìn)行剖分;由于線形渦流制動(dòng)系統(tǒng)會(huì)在鋼軌表面產(chǎn)生渦流,所以針對(duì)鋼軌要按照集膚深度進(jìn)行剖分。線性渦流制動(dòng)系統(tǒng)網(wǎng)格剖分圖如圖3所示。

        圖3 線性渦流制動(dòng)系統(tǒng)網(wǎng)格剖分Fig.3 Grid division of linear eddy current braking system

        2.2 制動(dòng)特性仿真

        常用制動(dòng)時(shí),線性渦流制動(dòng)系統(tǒng)勵(lì)磁電流I0≤78 A;緊急制動(dòng)時(shí)I0為95 A。氣隙長度l0變化范圍為6~8 mm,制動(dòng)特性試驗(yàn)參數(shù)如表2所示。

        表2 制動(dòng)特性試驗(yàn)參數(shù)

        2.2.1 勵(lì)磁電流對(duì)制動(dòng)力的影響

        設(shè)定不同勵(lì)磁電流I0研究其制動(dòng)力與速度變化之間的關(guān)系,如圖4所示。由圖4可知,不同勵(lì)磁電流下的曲線變化趨勢基本一致。當(dāng)氣隙長度保持不變時(shí),增加勵(lì)磁電流值,制動(dòng)力數(shù)值也隨之增大,說明勵(lì)磁電流與制動(dòng)力成正相關(guān),這也符合制動(dòng)力數(shù)學(xué)模型。在勵(lì)磁電流分別為60 A、78 A、95 A時(shí),其相應(yīng)的最大制動(dòng)力分別約為10.538 0 kN、15.100 9 kN、19.887 8 kN。由此可見,當(dāng)勵(lì)磁電流約增加30%時(shí),其制動(dòng)力增加了44%。

        圖4 不同勵(lì)磁電流下制動(dòng)力隨速度變化曲線

        2.2.2 氣隙長度對(duì)制動(dòng)力的影響

        保持勵(lì)磁電流不變,通過改變氣隙長度l0來觀察制動(dòng)力的變化,如圖5所示。

        圖5 不同氣隙長度下制動(dòng)力隨速度變化曲線

        由圖5可知,隨著氣隙長度的增加,其制動(dòng)力曲線峰值不斷下降。氣隙長度越小,制動(dòng)力曲線的變化趨勢越大,且制動(dòng)力呈明顯增大的趨勢。這是因?yàn)榇艌鰪?qiáng)度與磁路中磁阻大小有關(guān),氣隙長度增大造成了磁阻增大、磁場強(qiáng)度減小,使得渦流密度降低,因此制動(dòng)力有所下降。當(dāng)氣隙長度分別為6 mm、7 mm、8 mm時(shí),其相應(yīng)的最大制動(dòng)力分別為16.412 6 kN、15.100 9 kN、13.964 5 kN。由此可見,當(dāng)氣隙長度減小約33%時(shí),其制動(dòng)力減小了17.2%。

        此外,綜合圖4和圖5可知,渦流制動(dòng)力受速度影響較為明顯,速度越大,制動(dòng)力曲線上升速度越快,當(dāng)速度為10 km/h左右時(shí)就可以達(dá)到渦流制動(dòng)力的最大值,此時(shí)動(dòng)車速度稱為臨界速度。當(dāng)速度持續(xù)增大到高速區(qū)時(shí),渦流制動(dòng)力呈下降趨勢,最后趨于某個(gè)穩(wěn)定值,總體來說在高速狀態(tài)時(shí)呈良好的制動(dòng)狀態(tài)。

        2.2.3 勵(lì)磁線圈匝數(shù)對(duì)制動(dòng)力的影響

        為研究勵(lì)磁線圈匝數(shù)變化對(duì)制動(dòng)力的影響,選取列車速度在50 km/h、250 km/h、500 km/h時(shí)進(jìn)行分析,研究不同匝數(shù)對(duì)制動(dòng)力變化的影響。在相同的條件下,取勵(lì)磁線圈為280匝、380匝、480匝、580匝分別進(jìn)行計(jì)算,得到3個(gè)速度下不同勵(lì)磁線圈匝數(shù)的制動(dòng)力變化曲線,如圖6所示。

        圖6 不同列車速度下勵(lì)磁線圈匝數(shù)對(duì)制動(dòng)力變化的影響曲線

        由圖6可知,當(dāng)速度為50 km/h時(shí),制動(dòng)力增加速度最快;當(dāng)速度為500 km/h時(shí),制動(dòng)力增加速度最緩慢。由此可見,當(dāng)勵(lì)磁線圈匝數(shù)增加時(shí),低速運(yùn)行的列車制動(dòng)力增加更快,符合制動(dòng)力的變化規(guī)律。

        在列車運(yùn)行速度為50 km/h下,勵(lì)磁線圈匝數(shù)為580匝時(shí)的制動(dòng)力為33.500 0 kN;勵(lì)磁線圈匝數(shù)為280匝時(shí)的制動(dòng)力為10.700 0 kN。由此可見,勵(lì)磁線圈匝數(shù)增大107%,其制動(dòng)力增加了214%,可見制動(dòng)力增加的倍數(shù)遠(yuǎn)大于匝數(shù)增加的倍數(shù),即增加線圈匝數(shù)可以有效提高制動(dòng)力。

        2.2.4 鋼軌電導(dǎo)率對(duì)制動(dòng)力的影響

        保持勵(lì)磁電流為78 A、氣隙長度為7 mm、勵(lì)磁線圈匝數(shù)為280匝不變,將鋼軌電導(dǎo)率σ分別變?yōu)樵瓉淼?倍、4倍,并對(duì)模型進(jìn)行計(jì)算,得到不同磁導(dǎo)率下制動(dòng)力隨速度變化的曲線,如圖7所示。

        圖7 不同鋼軌電導(dǎo)率下制動(dòng)力隨速度變化曲線

        由圖7可知,當(dāng)鋼軌的電導(dǎo)率越大,渦流制動(dòng)力反而越小。由此可知,鋼軌的材料特性對(duì)渦流制動(dòng)力有較大的影響。選擇鋼軌材料時(shí),在符合工程條件的前提下,應(yīng)盡可能選擇電導(dǎo)率較小的材料,才能有更好的制動(dòng)性能。

        2.2.5 勵(lì)磁線圈材料對(duì)制動(dòng)力的影響

        在保持勵(lì)磁電流為78 A、氣隙長度為7 mm、勵(lì)磁線圈匝數(shù)為280匝時(shí),將勵(lì)磁線圈材料由銅換成鋁,將相同結(jié)構(gòu)的銅線圈與鋁線圈的渦流制動(dòng)系統(tǒng)產(chǎn)生的制動(dòng)力進(jìn)行對(duì)比,得到不同速度下的制動(dòng)力曲線,如圖8所示。

        圖8 采用銅線圈和鋁線圈時(shí)的制動(dòng)力隨速度的變化曲線

        從圖8中可以看出,在匝數(shù)相同的情況下,不同速度下銅線圈產(chǎn)生的渦流制動(dòng)力與鋁線圈的渦流制動(dòng)力變化趨勢一致,但鋁線圈的渦流制動(dòng)力比銅線圈的渦流制動(dòng)力略大,說明在線圈匝數(shù)、氣隙長度、線圈勵(lì)磁電流相同的條件下,選用鋁線圈獲得的渦流制動(dòng)力效果更好。

        2.3 電磁特性仿真

        本節(jié)主要分析列車速度為400 km/h時(shí),常用制動(dòng)和緊急制動(dòng)下的磁感應(yīng)強(qiáng)度分布情況。通過上述分析,線性渦流制動(dòng)系統(tǒng)參數(shù)選擇勵(lì)磁電流I0為78 A、氣隙長度l0為7 mm、勵(lì)磁線圈材料為鋁線圈。其分析結(jié)果如圖9—圖12所示。

        圖9 勵(lì)磁電流為78 A時(shí)渦流制動(dòng)系統(tǒng)磁場分布云圖(常用制動(dòng)工況)

        圖10 勵(lì)磁電流為78 A時(shí)氣隙磁場分布(常用制動(dòng)工況)

        圖11 勵(lì)磁電流為95 A時(shí)渦流制動(dòng)系統(tǒng)磁場分布云圖(緊急制動(dòng)工況)

        圖12 勵(lì)磁電流為95 A時(shí)氣隙磁場分布(緊急制動(dòng)工況)

        由圖9可知,在磁極附近磁場強(qiáng)度最強(qiáng),隨著與磁極距離越遠(yuǎn),磁場強(qiáng)度呈衰減趨勢。由于相鄰磁極之間形成閉合回路,在單個(gè)磁極上方磁軛處兩側(cè)的磁感應(yīng)強(qiáng)度大于中間的磁感應(yīng)強(qiáng)度。此外,在相鄰磁場的疊加作用下,鋼軌表面產(chǎn)生了較強(qiáng)的磁感應(yīng)強(qiáng)度。由圖11可知,緊急制動(dòng)工況下,磁感應(yīng)強(qiáng)度的分布情況與常用制動(dòng)工況下基本一致。同時(shí),由于緊急制動(dòng)工況下勵(lì)磁電流增大,使得渦流制動(dòng)系統(tǒng)的整體磁感應(yīng)強(qiáng)度有所增強(qiáng)。

        由圖10和圖12可知,兩種制動(dòng)工況下的磁感應(yīng)強(qiáng)度的變化趨勢基本一致,但緊急制動(dòng)工況的磁感應(yīng)強(qiáng)度值要大于常用制動(dòng)工況。氣隙長度處的磁場強(qiáng)度在相鄰磁極之間磁場較強(qiáng),在渦流制動(dòng)裝置兩側(cè)存在磁場泄露的問題,在短距離下會(huì)迅速衰減為0。在常用制動(dòng)情況下,氣隙長度處的磁場強(qiáng)度最高可達(dá)到877.640 kA/m,緊急制動(dòng)下可以達(dá)到1 030.595 kA/m。

        通過上述分析說明,緊急制動(dòng)工況下,渦流制動(dòng)系統(tǒng)可以產(chǎn)生更強(qiáng)的磁感應(yīng)強(qiáng)度和更高的制動(dòng)力,以保證列車的運(yùn)行安全。

        3 結(jié)論

        1) 本仿真研究得到的制動(dòng)力曲線變化趨勢符合渦流制動(dòng)力數(shù)學(xué)模型中各參數(shù)對(duì)制動(dòng)力的影響規(guī)律,從而驗(yàn)證了仿真模型的正確性。

        2) 渦流制動(dòng)力隨著速度的增加先增大后減小,在某個(gè)中間速度值達(dá)到臨界點(diǎn)。氣隙長度增加,渦流制動(dòng)力增大;勵(lì)磁電流增加,渦流制動(dòng)力增大。

        3) 當(dāng)勵(lì)磁線圈匝數(shù)增大后,渦流制動(dòng)力增加的倍數(shù)大于勵(lì)磁線圈匝數(shù)增加的倍數(shù),說明增加線圈匝數(shù)可以有效提高渦流制動(dòng)力。

        4) 當(dāng)勵(lì)磁線圈材料為鋁時(shí),與銅線圈相比,其制動(dòng)力有所增加。在符合工程的條件下,可以考慮將勵(lì)磁線圈材料選為鋁線圈。

        5) 當(dāng)增加鋼軌材料電導(dǎo)率時(shí),發(fā)現(xiàn)渦流制動(dòng)力反而下降,說明鋼軌材料電導(dǎo)率與渦流制動(dòng)力成負(fù)相關(guān)。因此在符合工程的條件下,選擇電導(dǎo)率較小的鋼軌材料可以獲得更好的渦流制動(dòng)效果。

        6) 緊急制動(dòng)工況下,渦流制動(dòng)系統(tǒng)可以提供更強(qiáng)的磁場、產(chǎn)生更大的渦流制動(dòng)力,以保障列車的運(yùn)行安全。

        猜你喜歡
        匝數(shù)氣隙渦流
        一種無刷雙饋電機(jī)匝數(shù)配合的半解析優(yōu)化方法*
        常用定轉(zhuǎn)子氣隙測量工具的設(shè)計(jì)及使用
        基于CFD仿真分析的各缸渦流比一致性研究
        基于Halbach陣列磁鋼的PMSM氣隙磁密波形優(yōu)化
        同步發(fā)電機(jī)理論的一個(gè)奇點(diǎn)與氣隙中心論
        不同側(cè)出線的變壓器線圈匝數(shù)計(jì)算
        渦流傳感器有限元仿真的研究與實(shí)施
        直線電機(jī)氣隙監(jiān)測技術(shù)的實(shí)際應(yīng)用與展望
        關(guān)于CW-系列盤式電渦流測功機(jī)的維護(hù)小結(jié)
        汽車科技(2014年6期)2014-03-11 17:46:08
        電渦流掃描測量的邊沿位置反演算法研究
        99久久免费看精品国产一| 久久天天躁夜夜躁狠狠| 亚洲 欧美 国产 制服 动漫| 男人靠女人免费视频网站| 国偷自产av一区二区三区| 久久aⅴ无码av高潮AV喷| 我的极品小姨在线观看| 色婷婷精品久久二区二区蜜臀av| 国内免费高清在线观看| 国产欧美日韩综合精品二区| 国产最新一区二区三区天堂| 99久久免费中文字幕精品| 日本人妻系列中文字幕| 色偷偷888欧美精品久久久| 国产a√无码专区亚洲av| 特级做a爰片毛片免费看108| 无码av免费永久免费永久专区 | 国产激情з∠视频一区二区| 国产亚洲成年网址在线观看| 久久国产精品免费一区二区三区| 天堂av在线美女免费| 国产在线观看无码免费视频| 两个人看的www中文在线观看| 国产三级国产精品三级在专区 | 69天堂国产在线精品观看| 久久网站在线免费观看| 丰满少妇人妻久久精品| 亚洲熟妇少妇任你躁在线观看无码 | 亚洲中文字幕人妻av在线| 末成年女a∨片一区二区| 大地资源中文第三页| 国产亚洲亚洲精品视频| 精品国产一区二区三区av麻| 欧美人与禽2o2o性论交| 亚洲乱码日产精品bd在线观看| 手机在线中文字幕国产| 亚洲国产精品色一区二区| 日本一区二区三区亚洲| 丁香美女社区| 久久精品国产亚洲av电影| 日本最新一区二区三区视频|