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

        ?

        絞吸挖泥船橫移系統(tǒng)多回路控制研究

        2016-12-30 02:24:36倪福生魏長(zhǎng)李洪彬
        自動(dòng)化儀表 2016年12期
        關(guān)鍵詞:絞刀挖泥船數(shù)學(xué)模型

        蔣 爽 倪福生 魏長(zhǎng) 李洪彬

        (河海大學(xué)疏浚教育和研究中心1,江蘇 常州 213022;河海大學(xué)機(jī)電工程學(xué)院2,江蘇 常州 213022)

        絞吸挖泥船橫移系統(tǒng)多回路控制研究

        (河海大學(xué)疏浚教育和研究中心1,江蘇 常州 213022;河海大學(xué)機(jī)電工程學(xué)院2,江蘇 常州 213022)

        為了提高絞吸挖泥船的施工產(chǎn)量,同時(shí)避免設(shè)備電流超限或管路濃度過(guò)高現(xiàn)象的發(fā)生,提出了混合多回路控制系統(tǒng)。建立了橫移系統(tǒng)的多個(gè)制約因素和橫移速度的數(shù)學(xué)模型,并分別設(shè)計(jì)了模糊PID和Smith預(yù)估器的控制方法。對(duì)于多個(gè)子回路之間的協(xié)調(diào)問(wèn)題,設(shè)計(jì)了一個(gè)多回路自動(dòng)切換控制器,并進(jìn)行仿真。仿真結(jié)果表明,該控制器能根據(jù)橫移系統(tǒng)的運(yùn)行狀態(tài)自動(dòng)切換控制回路,實(shí)現(xiàn)了絞吸挖泥船橫移系統(tǒng)對(duì)多回路的穩(wěn)定、高效控制。

        絞吸挖泥船 橫移控制 混合系統(tǒng) 多回路協(xié)調(diào) 自動(dòng)切換 數(shù)學(xué)模型 模糊PID控制 Smith預(yù)估器

        0 引言

        絞吸挖泥船在不同土質(zhì)、地形等工況下施工時(shí),制約生產(chǎn)能力的關(guān)鍵參數(shù)各不相同[1]。挖泥船在挖掘較硬土質(zhì)(如黏土、巖石)的工況時(shí),制約生產(chǎn)能力的關(guān)鍵參數(shù)是橫移阻力或絞刀切削阻力;當(dāng)挖掘較軟土質(zhì)(如淤泥、粉沙)時(shí),制約參數(shù)則變?yōu)楣艿垒斔湍酀{濃度或泥泵真空度。因此,在進(jìn)行挖泥船的橫移自動(dòng)化研究時(shí),需綜合考慮多個(gè)制約參數(shù),實(shí)時(shí)調(diào)整橫移速度,在確保絞吸挖泥船安全生產(chǎn)的同時(shí)盡量維持較高的生產(chǎn)能力。本文以當(dāng)前主流的電力變頻驅(qū)動(dòng)型絞吸挖泥船橫移系統(tǒng)為研究對(duì)象[2],選用橫移電機(jī)電流、絞刀電機(jī)電流和挖泥船產(chǎn)量作為控制參數(shù),分別建立其對(duì)橫移速度影響的數(shù)學(xué)模型,設(shè)計(jì)相應(yīng)的控制方法,并通過(guò)切換控制[3]決定由哪個(gè)回路進(jìn)行實(shí)際控制,實(shí)現(xiàn)了絞吸挖泥船橫移系統(tǒng)的穩(wěn)定、高效控制。

        1 絞吸挖泥船橫移系統(tǒng)數(shù)學(xué)模型建立

        1.1 考慮電流限制的數(shù)學(xué)模型

        由于橫移電機(jī)電流和絞刀電機(jī)電流與橫移速度之間具有相似的數(shù)學(xué)模型,將兩者統(tǒng)一起來(lái),建立了考慮電流限制的橫移數(shù)學(xué)模型。絞吸挖泥船施工時(shí),橫移電機(jī)驅(qū)動(dòng)橫移絞車(chē)收放纜索實(shí)現(xiàn)絞吸挖泥船繞定位樁的弧線(xiàn)擺動(dòng),同時(shí)絞刀電機(jī)驅(qū)動(dòng)絞刀旋轉(zhuǎn)切削,兩個(gè)運(yùn)動(dòng)相互配合實(shí)現(xiàn)扇形挖泥作業(yè)。橫移切削的電機(jī)電流與電機(jī)所受負(fù)載有關(guān),而電機(jī)負(fù)載的大小取決于絞刀切削反力。

        絞刀切削反力計(jì)算較為復(fù)雜,采用試驗(yàn)回歸方法,可得依據(jù)土質(zhì)計(jì)算的絞刀總切削阻力經(jīng)驗(yàn)公式為[4]:

        (1)

        (2)

        式中:Ca=0.4,Cv=1.0,Ch=0.9;r為絞刀大環(huán)內(nèi)圈半徑,m;M為絞刀頭扭矩,N·m。

        已知絞吸挖泥船三個(gè)方向上的切削反力,對(duì)整個(gè)挖泥船xy平面進(jìn)行受力分析,如圖1所示。圖中,F(xiàn)h為絞刀橫移切削時(shí)所受的水平方向作用力;Fa1、Fv1為切削反力Fa、Fv在xy平面上的投影;Fsbw、Fpsw為左右橫移錨纜拉力在xy面上的投影,F(xiàn)psw=cFsbw,c一般為0.3;γ1、γ2為左右橫移錨纜與船體中心線(xiàn)垂直方向的夾角;Fx總、Fy總、M總分別為船體受到的風(fēng)、流、波浪作用力和力矩。

        圖1 絞吸挖泥船整體xy平面受力圖

        Fig.1 The forces acting onxyplane of cutter suction dredger

        (3)

        (4)

        同樣可以得到絞刀電機(jī)的負(fù)載扭矩為:

        (5)

        因此,式(5)可以簡(jiǎn)化為:

        (6)

        絞吸挖泥船橫移和絞刀電機(jī)都采用三相異步電機(jī)。變頻器為V/F控制,即在低頻運(yùn)行時(shí),通過(guò)適當(dāng)補(bǔ)償定子電壓,可抵消阻抗壓降份額增大的影響,使反電動(dòng)勢(shì)和磁通得到補(bǔ)充,從而增大低頻時(shí)帶負(fù)載能力。在變頻器V/F控制下,電機(jī)運(yùn)行在恒磁通范圍內(nèi),電機(jī)電流可以表示為:

        (7)

        無(wú)功電流為:

        (8)

        有功電流為:

        (9)

        將式(4)、式(6)代入式(7)~式(9),分別得到橫移電機(jī)、絞刀電機(jī)和橫移速度之間的最終數(shù)學(xué)模型,如式(10)、式(11)所示[6]。

        考慮橫移電機(jī)電流限制的數(shù)學(xué)模型為:

        (10)

        考慮絞刀電機(jī)電流限制的數(shù)學(xué)模型為:

        (11)

        分析兩者的數(shù)學(xué)模型,可知其具有非線(xiàn)性的特點(diǎn),且土壤標(biāo)貫擊數(shù)和切刃有效長(zhǎng)度與土質(zhì)及地形密切相關(guān),具有時(shí)變性的特點(diǎn)。同時(shí),在上述建立模型的過(guò)程中,大部分公式都是依據(jù)經(jīng)驗(yàn)公式所得,無(wú)法獲得其準(zhǔn)確數(shù)學(xué)模型,具有不確定性特點(diǎn)。因此,考慮電流限制的橫移數(shù)學(xué)模型具有非線(xiàn)性、時(shí)變性及不確定性的特點(diǎn)。

        1.2 考慮產(chǎn)量限制的數(shù)學(xué)模型

        絞吸挖泥船產(chǎn)量形成可以分為三個(gè)過(guò)程:絞刀切削、管道輸送和產(chǎn)量計(jì)檢測(cè)。整個(gè)過(guò)程如圖2所示。

        圖2 絞吸挖泥船產(chǎn)量形成過(guò)程

        Fig.2 The forming process of construction output

        of cutter suction dredger

        (12)

        式中:YC為絞吸挖泥船濕方產(chǎn)量,t/h;VS為絞吸挖泥船橫移速度,m/min。增益K隨外界土質(zhì)、地形變化而變化,同樣具有時(shí)變性特點(diǎn),而且模型中滯后時(shí)間τ較長(zhǎng),τ/V>0.3,屬于大時(shí)滯系統(tǒng)。

        2 橫移系統(tǒng)多回路控制方法設(shè)計(jì)

        2.1 電流保護(hù)回路設(shè)計(jì)

        為了保證在提高絞吸挖泥船產(chǎn)量的同時(shí),橫移電機(jī)電流和絞刀電機(jī)電流不超限,需要建立電流保護(hù)回路。由于前述所建立的電流數(shù)學(xué)模型具有非線(xiàn)性、時(shí)變性和不確定性的特點(diǎn),常規(guī)PID控制難以達(dá)到較好的控制效果,考慮將模糊控制和傳統(tǒng)PID 控制相結(jié)合,得到模糊自適應(yīng)PID控制器。其中,模糊控制部分根據(jù)電流誤差E和誤差變化率Ec,實(shí)時(shí)在線(xiàn)自適應(yīng)調(diào)整PID控制中三個(gè)參數(shù):kp、ki和kd。PID控制部分采用智能積分算法,在誤差E較小并具有繼續(xù)減少的趨勢(shì)時(shí),消除積分作用,以防止積分飽和現(xiàn)象。

        2.2 產(chǎn)量控制回路設(shè)計(jì)

        為了在提高產(chǎn)量的同時(shí),避免挖掘淤泥、流沙時(shí)濃度過(guò)高造成管道堵塞,需要建立產(chǎn)量控制回路。由于前文所建立的產(chǎn)量數(shù)學(xué)模型具有大時(shí)滯、時(shí)變性和不確定的特點(diǎn),考慮采用Smith預(yù)估控制器對(duì)模型的純滯后予以補(bǔ)償,以減少滯后系統(tǒng)的超調(diào),增強(qiáng)系統(tǒng)穩(wěn)定性。但傳統(tǒng)的Smith預(yù)估控制器系統(tǒng)對(duì)模型估計(jì)準(zhǔn)確性要求較高,在模型不準(zhǔn)確時(shí),控制性能會(huì)變差甚至不穩(wěn)定。因此,需要在傳統(tǒng)Smith預(yù)估器的基礎(chǔ)上對(duì)其進(jìn)行改進(jìn)。

        實(shí)際對(duì)象參數(shù)中增益K和時(shí)滯常數(shù)τ對(duì)Smith 預(yù)估器穩(wěn)定性的影響較大,對(duì)時(shí)間常數(shù)T的影響較小。對(duì)絞吸挖泥船系統(tǒng)而言,其滯后時(shí)間τ主要取決于泥沙從吸口輸送到泥沙產(chǎn)量計(jì)處的時(shí)間,這部分時(shí)間可以由計(jì)算得到且變化很小。增益K主要取決于絞刀切削過(guò)程的漏泥率,其波動(dòng)性較大,是影響系統(tǒng)控制性能的關(guān)鍵因素。本文通過(guò)改進(jìn)傳統(tǒng)的Smith 預(yù)估控制方法,系統(tǒng)引入實(shí)際對(duì)象的輸出產(chǎn)量和預(yù)估器輸出產(chǎn)量的誤差來(lái)在線(xiàn)調(diào)整預(yù)估器的系統(tǒng)增益K,使其能跟蹤并收斂于實(shí)際對(duì)象參數(shù)增益。改進(jìn)后的Smith方法結(jié)構(gòu)框圖如圖3所示。其主控制器采用模糊PID控制算法,以增加系統(tǒng)自適應(yīng)性能。

        圖3 改進(jìn)Smith預(yù)估器控制系統(tǒng)框圖

        Fig.3 The improved Smith predictor control system

        參數(shù)自適應(yīng)辨識(shí)模塊由模糊控制和參數(shù)辨識(shí)兩部分組成。模糊控制部分取實(shí)際對(duì)象的輸出y和預(yù)估器輸出的誤差ep作為輸入量,增益K的調(diào)節(jié)量為輸出量。對(duì)輸入/輸出變量進(jìn)行模糊化處理,輸入/輸出變量的語(yǔ)言集都為{NB,NM,,NS,ZO,PS,PM,PB},隸屬度函數(shù)采用靈敏性強(qiáng)的三角函數(shù)。參數(shù)辨識(shí)部分中,令θ為可調(diào)參數(shù),記為θn,選擇性能指標(biāo)為f(x)=|ep|,當(dāng)滿(mǎn)足f(x)≤0.5時(shí),即認(rèn)為達(dá)到參數(shù)的調(diào)節(jié)要求,不再對(duì)增益K進(jìn)行調(diào)節(jié)。調(diào)節(jié)后的增益為K調(diào)=Kθn。

        3 多回路自動(dòng)切換控制器設(shè)計(jì)

        混合系統(tǒng)已經(jīng)成為控制界的研究熱點(diǎn)[7-8],而切換控制作為混合系統(tǒng)的一種有效控制方法,已經(jīng)得到了實(shí)際應(yīng)用[9-10]。本文采用自動(dòng)切換控制的方法,設(shè)計(jì)了一種多回路自動(dòng)切換控制器,可在不同的工況下選擇不同的控制回路,以實(shí)現(xiàn)絞吸挖泥船橫移自動(dòng)化控制。產(chǎn)量與電流保護(hù)回路切換示意圖如圖4所示。

        圖4 產(chǎn)量與電流保護(hù)回路切換示意圖

        Fig.4 The switchover of current protection loops

        and construction output

        當(dāng)絞吸挖泥船產(chǎn)量較低時(shí),為提高生產(chǎn)效率,需要提高橫移速度,由產(chǎn)量控制回路進(jìn)行控制。從工況點(diǎn)A運(yùn)行至工況點(diǎn)B時(shí),由于橫移速度的增加會(huì)引起橫移電機(jī)和絞刀電機(jī)電流上升,可能達(dá)到電流限定值。在工況點(diǎn)B處,電流保護(hù)回路接替產(chǎn)量控制回路開(kāi)始工作,以電流保護(hù)回路目標(biāo)值作為整個(gè)系統(tǒng)的控制目標(biāo)值來(lái)調(diào)整橫移速度。子回路實(shí)際值與設(shè)定值的差值即誤差值越大,該回路越危險(xiǎn)。當(dāng)所有子回路都低于或都高于其目標(biāo)設(shè)定值時(shí),則以其中最為危險(xiǎn)的子回路目標(biāo)值作為整個(gè)系統(tǒng)的控制目標(biāo)值。

        同時(shí),絞吸挖泥船施工時(shí)外部環(huán)境較為復(fù)雜,疏浚地形、地質(zhì)多變,這些外部的干擾都會(huì)引起絞吸挖泥船產(chǎn)量和電機(jī)電流的瞬時(shí)波動(dòng)變化,進(jìn)而引起產(chǎn)量與電流保護(hù)回路頻繁切換甚至造成系統(tǒng)不穩(wěn)定。為提高切換控制系統(tǒng)的穩(wěn)定性和抗干擾性,避免回路間頻繁切換影響系統(tǒng)控制性能,本文在設(shè)計(jì)的自動(dòng)切換控制系統(tǒng)中加入了一個(gè)濾波時(shí)間t濾波,即在t濾波內(nèi)需作多次判斷以確認(rèn)絞吸挖泥船產(chǎn)量、橫移電機(jī)電流或絞刀電機(jī)電流中某一值是否超過(guò)系統(tǒng)限定值。當(dāng)產(chǎn)量或電流的超限時(shí)間大于t濾波時(shí)進(jìn)行切換(如圖4中BC工況中t1和t2均大于t濾波),否則不進(jìn)行切換(如圖4中DE工況中t3和t4均小于t濾波),從而消除信號(hào)突變或干擾的影響,其中時(shí)間的設(shè)定應(yīng)低于系統(tǒng)超限所允許的時(shí)間。自動(dòng)切換控制器以各回路誤差作為輸入信號(hào),輸出為回路切換信號(hào)。整個(gè)自動(dòng)切換控制系統(tǒng)框圖如圖5所示。

        圖5 絞吸挖泥船橫移自動(dòng)切換控制系統(tǒng)框圖

        Fig.5 Automatic traverse switching control system of CSD

        4 多回路控制系統(tǒng)的仿真研究

        本文選取外部環(huán)境因素中的土質(zhì)軟硬程度變化作為干擾因素,對(duì)多回路控制系統(tǒng)進(jìn)行仿真,設(shè)定t濾波=1 s。在t濾波=20 s時(shí),引入土壤標(biāo)貫擊數(shù)外界干擾曲線(xiàn),如圖6所示。

        圖6 土壤標(biāo)貫擊數(shù)變化曲線(xiàn)

        Fig.6 Variation curve of the soil SPT blow count

        在土質(zhì)軟硬程度即土壤標(biāo)貫擊數(shù)N變化的外界干擾下,子回路相對(duì)誤差及切換曲線(xiàn)如圖7所示。

        圖7 子回路相對(duì)誤差及切換曲線(xiàn)

        Fig.7 The relative error and switching curve of sub-loops

        圖7(a)中,誤差大于0即表示該子回路實(shí)際值超過(guò)系統(tǒng)設(shè)定值。圖7(b)中,當(dāng)狀態(tài)值為-1時(shí),說(shuō)明目前處于絞刀電機(jī)電流保護(hù)回路;當(dāng)狀態(tài)值為0或1時(shí),說(shuō)明目前處于橫移電機(jī)電流保護(hù)回路或產(chǎn)量控制回路。多回路控制系統(tǒng)以各子回路的誤差作為輸入信號(hào),圖7(a)中,當(dāng)0

        電機(jī)電流保護(hù)回路,即狀態(tài)變量由1變?yōu)?1。同樣,在t=38.1 s和t=40.5 s時(shí)進(jìn)行了回路切換,在t=54 s時(shí),雖然挖泥船產(chǎn)量曲線(xiàn)誤差大于0,但由于其時(shí)間較短(t<1 s),因而其被當(dāng)作信號(hào)干擾而被過(guò)濾掉,系統(tǒng)不進(jìn)行切換。

        按照以上規(guī)律,系統(tǒng)判斷確認(rèn)各子回路誤差大小并進(jìn)行自動(dòng)切換,實(shí)現(xiàn)整個(gè)絞吸挖泥船橫移混合多回路系統(tǒng)的穩(wěn)定控制。

        5 結(jié)束語(yǔ)

        本文建立了電力驅(qū)動(dòng)型絞吸挖泥船橫移速度和橫移電機(jī)電流、絞刀電機(jī)電流及產(chǎn)量之間的數(shù)學(xué)模型,設(shè)計(jì)了各控制回路的控制方法。針對(duì)混合控制系統(tǒng)中的多回路協(xié)調(diào)問(wèn)題,設(shè)計(jì)了多回路自動(dòng)切換控制系統(tǒng),并對(duì)其進(jìn)行了仿真。仿真結(jié)果表明,多回路自動(dòng)切換控制器在施工工況(如土質(zhì))變化較大的情況下,能夠根據(jù)各子回路誤差信號(hào),自動(dòng)切換至適合該工況條件的控制器,實(shí)現(xiàn)了橫移系統(tǒng)在多種工況下的高效、穩(wěn)定控制。系統(tǒng)安全性和自適應(yīng)性好,對(duì)實(shí)現(xiàn)挖泥船的橫移自動(dòng)化具有一定的指導(dǎo)意義。

        [1] 肖博,孫守勝.絞吸式挖泥船橫移系統(tǒng)加速施工[J].中國(guó)港灣建設(shè),2015,35(9):77-80.

        [2] 任經(jīng)緯.絞吸式挖泥船集成控制系統(tǒng)的構(gòu)建[D].上海:上海交通大學(xué),2010.

        [3] 于達(dá)仁,劉曉鋒,隋巖峰.考慮切換特性的航空發(fā)動(dòng)機(jī)多回路控制系統(tǒng)設(shè)計(jì)[J].航空動(dòng)力學(xué)報(bào),2007(8):1378-1383.

        [4] 李洪彬.挖泥船絞刀三維參數(shù)化建模及應(yīng)用分析[D].常州:河海大學(xué),2009.

        [5] 潘英杰.絞吸式挖泥船絞刀架載荷計(jì)算和強(qiáng)度分析[D].上海:上海交通大學(xué),2008.

        [6] 谷孝利.基于變頻技術(shù)的大型絞吸式挖泥船電力驅(qū)動(dòng)控制系統(tǒng)研究[D].上海:上海交通大學(xué),2011.

        [7] LENNARTSON B,TITTUS M,EGARD B,et al.Hybrid systems in process control[J].Control System Maginze,1996,16(5):45-56.

        [8] CASSANDRAS C G,PEPYNE D L,WARDI Y.Optimal control of a class of hybrid systems control[J].Automatic Control,2001,46(3):398-415.

        [9] CHEN S S,CHANG Y C.Fuzzy control design for switched nonlinear systems[C]//SICE Annual Conference.Okayama:Eileen Breen,2008:352-357.

        [10]仇小杰,黃金泉,魯峰,等.航空發(fā)動(dòng)機(jī)多路模糊切換控制系統(tǒng)設(shè)計(jì)[J].南京航空航天大學(xué)學(xué)報(bào),2011,43(4):439-445.

        Research on the Multi-loop Control for Traverse System of Cutter Suction Dredger

        In order to increase the construction output of cutter suction dredgers,and avoid the over limit current of equipment or excessive concentration in pipeline,a hybrid multi-loop control system is proposed.The mathematical models of a plurality of restrictive factors and the traverse speed are built,and the control methods of fuzzy PID and Smith predictor are respectively designed.For the coordination among multiple sub-loops,a multi-loop automatic switching controller is designed and simulated.The simulation results show that the controller can automatically switchover the control loops according to the running state of the traverse system,thus the stable and highly efficient control of the traverse system of cutter suction dredger is realized.

        Cutter suction dredger Swing control Hybrid system Multi-loop coordination Automatic switching Mathematical model Fuzzy PID control Smith predictor

        蔣爽(1981—),男,2007年畢業(yè)于河海大學(xué)材料加工專(zhuān)業(yè),獲碩士學(xué)位,實(shí)驗(yàn)師;主要從事疏浚技術(shù)自動(dòng)化方向的研究。

        TH39;TP273

        A

        10.16086/j.cnki.issn 1000-0380.201612016

        修改稿收到日期:2016-05-11。

        猜你喜歡
        絞刀挖泥船數(shù)學(xué)模型
        AHP法短跑數(shù)學(xué)模型分析
        活用數(shù)學(xué)模型,理解排列組合
        絞吸挖泥船大功率挖巖絞刀的荷載分析*
        絞吸挖泥船絞刀長(zhǎng)軸系軸承支反力計(jì)算方法研究
        絞吸挖泥船大功率挖巖絞刀設(shè)計(jì)與研究
        對(duì)一個(gè)數(shù)學(xué)模型的思考
        高壓潛水電動(dòng)機(jī)在絞吸挖泥船上的應(yīng)用管理
        7025型絞吸式挖泥船串聯(lián)施工工藝試驗(yàn)
        “吸盤(pán)3”號(hào)挖泥船核子密度計(jì)改進(jìn)設(shè)計(jì)
        船海工程(2015年4期)2016-01-05 15:53:50
        絞吸式挖泥船的穩(wěn)性計(jì)算
        廣東造船(2015年6期)2015-02-27 10:52:45
        亚洲国产大胸一区二区三区| 99久久精品日本一区二区免费| 无码av免费精品一区二区三区 | 在线观看黄片在线播放视频| 丝袜美腿诱惑区在线播放| 蜜桃18禁成人午夜免费网站| 亚洲精品久久久久成人2007| 精品无码国产自产野外拍在线| 国产美女裸身网站免费观看视频| 一区二区免费国产a在亚洲| 国产极品大奶在线视频| 国内精品久久久久影院优| 激性欧美激情在线| 欧美在线资源| 高清av一区二区三区在线| 国产高清一区二区三区三州| 天天摸天天做天天爽水多| 94久久国产乱子伦精品免费| 天天躁日日操狠狠操欧美老妇 | 天天做天天躁天天躁| 婷婷一区二区三区在线| 精品中文字幕久久久人妻| 国产亚洲精品97在线视频一| 99久久精品午夜一区二区| 日日干夜夜操高清视频| 国产好片日本一区二区三区四区| 精品一区二区三区亚洲综合| 粉嫩av国产一区二区三区| 真人与拘做受免费视频| 久久婷婷是五月综合色狠狠| 亚洲情久久久精品黄色| 成人影院在线视频免费观看| 国产探花在线精品一区二区| 久久久久国产一级毛片高清版A | 中国妇女做爰视频| 亚洲欧美日韩中文在线制服| 四虎国产精品免费久久麻豆| 日本在线一区二区免费| 人妻av中文字幕久久| 久久水蜜桃亚洲av无码精品麻豆 | 亚洲а∨精品天堂在线|