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

    墜物沖擊下基于混合代理模型的水下管匯動力學分析方法

    2023-11-04 02:24:40任濤郭喜慶曾威陳旭鵬
    科學技術與工程 2023年29期
    關鍵詞:墜物管匯代理

    任濤, 郭喜慶, 曾威, 陳旭鵬

    (西安石油大學機械工程學院, 西安 710065)

    水下管匯具有油氣匯集、電力分配和水下控制等功能,是水下油氣生產系統(tǒng)的重要節(jié)點,其工作性能直接影響海上油氣開采的穩(wěn)定性和安全性。但是,水下管匯在服役過程中受到來自漁船、平臺作業(yè)船只墜物沖擊影響,威脅水下油氣生產安全。因此,準確分析墜物沖擊下水下管匯主結構的動力學性能,為水下管匯結構設計提供可靠依據(jù),是提升水下管匯防護能力、保護其工作穩(wěn)定性的有效途徑。為了研究海上墜物對水下管匯的影響程度,相關學者通過建立水下管匯有限元模型進行了靜力學和動力學分析。例如,余峙偉等[1]為了加快海上墜物對其主結構的碰撞有限元求解過程,采用簡化水下管匯主結構模型的方法,提高了撞擊載荷下水下管匯主結構的求解速率。婁敏等[2]利用全積分殼單元算法建立管道理論模型,對墜物與水下管道發(fā)生的碰撞過程進行了分析研究。魯中岐等[3]在解決海底管道的腐蝕問題時,采用代理模型方法代替大量仿真試驗進行分析研究。孫寶等[4]將代理模型技術引入黏彈性緩沖系統(tǒng)參數(shù)優(yōu)化中,減小了計算復雜度并提高了優(yōu)化結果。上述研究改善了墜物沖擊下水下管匯主結構碰撞的動力學性能,但考慮其原比例模型大和結構復雜的問題,分析過程中均采用了降低尺寸比例和簡化模型等方法,一定程度上影響了分析結果準確度[5]。同時,沖擊載荷下相關結構的動力學分析屬于典型的強非線性問題,當理論模型結構尺寸大且足夠復雜時,單次模擬的時間成本特別大,這限制了傳統(tǒng)的有限元方法在水下管匯動力學分析中的應用[6]。然而在解決高復雜性工程問題時,往往會導致有限元方法計算量大以及耗時較長,基于此引入代理模型技術可以很大程度提高計算效率,降低時間成本。

    基于此,現(xiàn)針對上述采用簡化模型和縮小比例等方法導致水下管匯動力學分析精度受損、分析結果不準確的問題,建立水下管匯的全尺寸結構有限元模型進行動力學分析,并針對由此造成計算成本過高的問題,引入代理模型方法提高水下管匯動力學分析效率。進一步地,在代理模型建模過程中,考慮單一代理模型擬合性能的局限性,基于組合預測理論構建水下管匯動力學響應混合代理模型,開展基于混合代理模型的水下管匯動力學分析,在保證分析精度的同時進一步提高其優(yōu)化效率。

    1 單一代理模型基本理論

    1.1 Kriging模型

    Kriging(KG)代理模型[7]是一種根據(jù)變異函數(shù)理論構建的插值模型。其在滿足插值方差最小的條件下,給出最優(yōu)線性無偏插值,是一種估計方差最小化的無偏估計模型。其函數(shù)模型表達式為

    (1)

    F(β,x)=β1f1(x)+β2f2(x)+…+βpfp(x)

    =fT(x)β

    (2)

    式中:β為線性回歸系數(shù);f(x)為x變量的多項式函數(shù);z(x)為一個隨機過程。由于其建模過程較復雜,具體建模過程見參考文獻[7]。

    1.2 RBF模型

    徑向基函數(shù)(radical basis function,RBF)代理模型[8]是一種利用離散多遠數(shù)據(jù)擬合未知非線性函數(shù)的技術。它的基函數(shù)是典型徑向函數(shù),將徑向堆成基函數(shù)的自變量設為已知樣本點和未知待測點之間的Euclidean距離,徑向函數(shù)組合得到其模型,也可稱為一系列關于樣本點的基函數(shù)線性和,即為徑向基函數(shù)模型。其函數(shù)模型表達式為

    (3)

    式(3)中:‖·‖表示歐式距離;Ψ=ψ(‖x-x(i)‖)為基函數(shù);wi為i個樣本點的權系數(shù)。

    1.3 支持向量回歸

    支持向量回歸(support vector regression,SVR)代理模型[9]是一種與支持向量機(support vector machine,SVM)相同原理的監(jiān)督學習算法。與之不同的是支持向量回歸是使用核、稀疏解和VC控制邊距和支持向量的數(shù)量,它的主要優(yōu)點是其計算復雜度不依賴于輸入空間的維度。此外,具有出色的泛化能力,具有很高的預測精度。其函數(shù)模型形式為

    (4)

    式(4)中:μ為一個給定常量;ωi為i個樣本點的權系數(shù);φ為樣本集的基本函數(shù)。

    2 水下管匯主結構碰撞有限元模型

    水下管匯主結構不但能夠為各類型生產管道和閥門組件提供支撐平臺,而且還能夠防止海上墜物沖擊造成的損傷。在受到墜物撞擊時,管匯主結構的基礎框架通過變形來承受墜物產生的動載荷。因此,選擇水下管匯主結構作為分析對象,并通過適當簡化,建立水下管匯結構模型,如圖1(a)所示。為降低全尺寸模型的計算成本,并考慮水下管匯主結構的受力特征,在有限元軟件中分別使用梁單元和殼單元對管匯主結構中的H型鋼和圓柱支撐柱進行有限元劃分;墜物和防護板則采用實體單元進行模擬,并定義管匯防護板材料為A36,其他結構材料為Q235。由于防護板是承受墜物沖擊的實際結構,因此在劃分網格時對防護板的有限元網格進行了適當加密。最終,有限元模型如圖1(b)所示。

    圖1 水下管匯主結構實體模型和有限元模型Fig.1 Solid model and finite element model of main structure of subsea manifold

    海上墜物墜落過程中,在受到海水阻力作用下,其經過一定的墜落高度后會呈現(xiàn)一種固定的速度和姿態(tài)與防護板進行接觸。因此,在碰撞分析過程中認為墜物從管匯正上方墜落且呈固定的速度和姿態(tài)。其中,根據(jù)DNV規(guī)范,墜物的墜落速度采用式(5)進行計算[10],即

    (5)

    式(5)中:m為墜物的質量,kg;V為墜物的體積(排出水的體積),m3;ρw為水的密度;g為重力加速度,取值9.81 g/m2;CD為阻力系數(shù);A為墜物在流向上的投影面積,m2;vT為墜物的速度,m/s。

    海上墜物通常是海上作業(yè)漁船、平臺等上的落物,根據(jù)船舶航運標準規(guī)范,主要包括船錨、球形和管型等墜物,因此分別采用圓管、霍爾錨和球體分別模擬3種墜物類型模擬。同時,考慮圓管結構墜落姿態(tài)包括縱向和橫向兩種情況,因此分別模擬圓管(縱向)、圓管(橫向)、霍爾錨和球體4種墜物沖擊下管匯的動力學響應[11]。如圖2所示,為圓管、霍爾錨和球體質量分別是700、1 000、2 000 kg時,4種墜物沖擊下管匯主結構的應力云圖。

    圖2 不同海上墜物撞擊管匯主結構應力云圖Fig.2 Stress nephogram of main structure of manifold impacted by different marine falling objects

    圖2為圓管縱向、圓管橫向、霍爾錨、球體4種墜物與管匯撞擊過程最大等效應力時刻,最大等效應力分別為148.14、162.8、201.32、312.75 MPa。由以上數(shù)據(jù)可知,當4種墜物作為研究對象時,墜物質量越大,碰撞產生的應力越大。由圓管墜物的仿真結果可知,在墜物質量相同的條件下,墜物與管匯碰撞接觸面積越小產生的應力越大。該分析在CPU為i7 G4560處理器、64位操作系統(tǒng)以及8 G運行內存的平臺上完成,單次有限元動力學分析耗時約48 h。

    3 水下管匯動力學指標代理模型建模與預測分析

    3.1 代理模型設計變量與響應量

    選取墜物霍爾錨、圓管、球體的質量為代理模型的設計變量。由設計要求可知,設計變量的取值范圍為初始值上、下波動的1/2,表1為其具體取值區(qū)間。

    表1 代理模型設計變量區(qū)間Table 1 Value ranges of the surrogate model variables

    根據(jù)水下管匯動力學分析目的,選取墜物與水下管匯碰撞的最大應力值σ、最大應變值ε、能量吸收e以及管匯連接器所在位置的最大加速度a為動力學響應量,構建水下管匯碰撞動力學響應代理模型。

    3.2 水下管匯動力學響應KG模型

    為保證KG代理模型預測精度,對于高計算成本的管匯結構動力學響應分析,需要選擇數(shù)量至少為2k的訓練樣本[12],即

    k=(n+1)(n+2)/2

    (6)

    式(6)中:n為KG模型的設計變量數(shù)目。

    由于設計變量數(shù)為1,因此至少需要的訓練樣本數(shù)為6。為保證樣本均勻性,抽樣過程中采用拉丁超立方試驗設計方法在表1所示區(qū)間中進行抽樣,并采用構建的水下管匯主結構碰撞有限元模型計算對應的動力學響應指標,得到代理模型訓練樣本如表2所示。

    表2 代理模型訓練樣本Table 2 Training sample of the surrogate model

    采用最大似然估計方法,分別計算圓管(縱向)、圓管(橫向)、霍爾錨、球體撞擊下水下管匯動力學響應指標KG模型的相關參數(shù),利用式(1)即可構建得到不同動力學指標響應KG模型。

    3.3 水下管匯動力學響應RBF模型

    以表2中數(shù)據(jù)為RBF模型的訓練樣本,以Gaussian函數(shù)φ(r,c) = exp(-r2/c2)為基函數(shù)。在此基礎上,使用Gaussian函數(shù)φ(r,c)構建100×100的矩陣φ,并采用公式ω=y′/φ計算得到預測響應值的權系數(shù)ω,其中為y′為訓練樣本中的水下管匯動力學響應指標值,即y′=[y1,y2,…,y100]。最終,墜物沖擊下水下管匯的動力學響應RBF模型為

    (7)

    式(7)中:xn為待預測點。

    3.4 水下管匯動力學響應SVR模型

    首先定義一個輸入值xi和輸出值yi的回歸函數(shù)f(x),選取表2中數(shù)據(jù)為SVR模型的訓練樣本集,并將樣本坐標引入松弛變量ζi*和ζi,以確定樣本的偏離約束程度,其表達式為

    (8)

    然后引入Largrange函數(shù)乘子α,有

    (9)

    根據(jù)式(9)計算得到最優(yōu)解α={α1,α2,…,

    αn},進而得到水下管匯動力學響應指標的權重參數(shù)和偏置量。最終,墜物沖擊下水下管匯的動力學響應SVR模型為

    (10)

    式(10)中:α和α*為對偶變量;K(xi,xj)為核函數(shù)。

    分別使用構建的RBF、KG、SVR3種代理模型對海上墜物撞擊水下管匯主結構的應力值、應變值、連接器處的加速度、能量吸收進行預測,結果如圖3所示。

    圖3 不同單一代理模型預測結果Fig.3 Forecast results of different single surrogate models

    3.5 3種代理模型精度分析

    (11)

    表3 單一代理模型誤差對比Table 3 Error comparisons of single surrogate model

    根據(jù)表3誤差對比結果,與RBF和KG模型相比,SVR模型預測值對于應力值、應變值、加速度、能量吸收的平均誤差最大,分別為20.95%、25.47%、26.71%、25.97%,這表明RBF和KG模型對于墜物沖擊下管匯動力學響應指標的預測更加準確。

    4 基于混合代理模型的水下管匯動力學分析

    4.1 水下管匯動力學響應混合代理模型建模

    混合代理模型是將兩種及以上的單一代理模型組合起來,構建一個有更高精度的預測模型,且其計算精度一般高于單一的代理模型[13],其基本思想是充分利用單一代理模型的預測能力,使混合代理模型擁有多個單一代理模型的預測高精度性?;旌洗砟P偷慕Y構一般可以描述為

    (12)

    式(12)中:yE(x)為混合代理模型在預測點x處的預測值;Ns為單個代理模型的個數(shù);wi為對應代理模型的權重值;yi(x)為第i個代理模型在預測點x處的預測值。

    由式(12)可知,權重比例直接決定了單一代理模型預測值的占比大小。根據(jù)模型混合理念,單一代理模型預測精度越高,其在混合代理模型中所占權重比越大,反之,預測精度越低,在混合代理模型中占比越小。因此,當權重比例合適時,混合代理模型的預測精度可以達到最大化。為提高代理模型的預測精度,采用混合代理模型對不同海上墜物與水下管匯的碰撞動力學響應結果進行預測。根據(jù)表3中代理模型的預測誤差對比結果,選擇預測精度更高的RBF和KG模型作為混合代理模型的子模型[14],并采用平均加權法進行組合?;旌洗砟P捅磉_式為

    Y(x)=ω1φ1(x)+ω2φ2(x)

    (13)

    式(13)中:x為自變量;φ1為徑向基代理模型;φ2為克里金代理模型;ω1和ω2分別為兩個模型的權重系數(shù)。

    子模型的權重系數(shù)直接影響著混合代理模型預測值的精度[15],為保證模型精度,采用克拉默法則進行權重計算,其具體計算公式如下。

    (14)

    AΩ=β

    (15)

    式中:A=[aij]n×n,Ω=(ωi1,ωi2,…,ωin)T,β=(b1,b2,…,bn)T,A為各代理模型預測值系數(shù)矩陣,Ω為由權重系數(shù)組成的列向量,β為由真實值組成的列向量。

    計算過程中,根據(jù)兩個代理模型在同一個已知點的預測值,分別給予一個未知系數(shù),聯(lián)合該點的真實值得到權重系數(shù)的方程組。將每個點逐個測驗,最終計算出每組的權重系數(shù),混合代理模型權重系數(shù)如表4所示。

    表4 混合代理模型權重Table 4 Weights of the hybrid surrogate model

    基于表4中混合代理模型子模型權系數(shù)計算結果,可以構建得到不同動力學響應指標對應的混合代理模型,其基本形式為

    (16)

    4.2 混合代理模型預測精度

    為檢驗混合代理模型的預測精度,將混合代理模型預測結果分別與RBF、KG代理模型和真實值進行對比,結果如圖4所示。進一步地,采用式(11)所示平均相對誤差為評價指標計算各代理模型的預測誤差,具體結果如表5所示。

    表5 動力學響應指標代理模型預測誤差Table 5 Prediction errors of the surrogate models for dynamic responses

    圖4 混合代理模型預測值對比Fig.4 Comparison of predicted values of hybrid surrogate model

    由圖4可以看出,混合代理模型相比于單一代理模型的預測曲線更接近于真實值曲線。根據(jù)表5結果可知,對于墜物沖擊下水下管匯主結構應力值、加速度和能量吸收預測,RBF預測精度最低,平均相對誤差分別為14.5%、19.38%和15.11%;KG模型預測精度次之,平均相對誤差分別為9.3%、8.35%和8.64%;混合代理模型預測精度最高,平均相對誤差分別為6.02%、7.92%和7.79%。對于應變值預測,KG模型預測精度相對最低,平均相對誤差為7.81%;RBF次之,平均相對誤差為6.85%;混合代理模型預測精度最高,平均相對誤差為6.54%??偟膩碇v,對于墜物沖擊下水下管匯主結構應力值、應變值、加速度和能量吸收,混合代理模型均具有最高的預測精度。另外,根據(jù)表5還可以計算混合代理模型對應各響應指標的總體平均預測誤差為7.07%,根據(jù)表3可計算出3種單一代理模型的總體平均誤差為17.19%。這表明相比單一代理模型,混合代理模型平均預測精度提高了10.12%。

    4.3 混合代理模型計算成本

    為驗證混合代理模型的計算成本,設動力學分析需要使用全尺寸水下管匯有限元模型進行32次計算,而建立單個代理模型需要以6次有限元計算為樣本,則采用不同方法完成該動力學分析的計算成本如表6所示。

    表6 不同方法計算成本對比Table 6 Time consumption comparisons of different methods

    由表6可以看出,相較于有限元仿真,采用代理模型進行動力學分析能夠顯著提高計算效率。其中,采用單一代理模型能夠將計算效率提高81.25%;采用混合代理模型能夠將計算效率提高81.24%,兩者計算成本基本一致。但是,根據(jù)預測精度對比結果,混合代理模型的預測精度相對單一代理模型提高了10.12%。

    5 結論

    針對海上墜物與水下管匯碰撞過程全尺寸模型模擬仿真時間過長的問題,引入代理模型開展墜物沖擊下水下管匯動力學分析,以此為基礎建立混合代理模型,在保證動力學分析精度的同時提高了分析效率,主要結論如下。

    (1) 采用代理模型方法能夠提高墜物沖擊下水下管匯動力學分析效率。其中,SVR代理模型預測精度較差,平均誤差達到24.78%,KG和RBF代理模型預測精度較高,平均誤差分別為8.58%和13.96%。

    (2) 與單一代理模型相比,混合代理模型建模成本與單一代理模型基本一致,但是其平均計算精度提高了10.12%,在水下管匯結構優(yōu)化、可靠性分析等需要重復進行有限元仿真時建議采用混合代理模型進行數(shù)據(jù)預測,可以有效避免實際工程中的海量計算問題。

    猜你喜歡
    墜物管匯代理
    高空墜物要當心
    兒童時代(2023年2期)2023-03-28 06:47:20
    海底光電復合纜受墜物撞擊損傷分析
    高空墜物
    幼兒100(2022年11期)2022-03-15 01:30:22
    法律中高空拋物墜物行為的責任承擔
    法制博覽(2021年13期)2021-11-26 14:51:00
    基于模糊綜合評價的水下管匯結構可靠性分析*
    石油機械(2021年4期)2021-04-23 08:08:26
    代理圣誕老人
    代理手金寶 生意特別好
    拋錨撞擊水下管匯的數(shù)值模擬研究
    海洋工程(2016年4期)2016-10-12 03:21:26
    復仇代理烏龜君
    學生天地(2016年23期)2016-05-17 05:47:15
    番禺35-1/35-2氣田水下管匯嵌入式在線安裝方案設計及建造關鍵技術
    成人国产麻豆网| 国产精品久久久久久av不卡| 秋霞在线观看毛片| 欧美97在线视频| www.av在线官网国产| 欧美激情国产日韩精品一区| 国产成人精品久久久久久| 欧美精品一区二区大全| 美女主播在线视频| 丰满迷人的少妇在线观看| 国产成人aa在线观看| 在线观看人妻少妇| 九九在线视频观看精品| 免费不卡的大黄色大毛片视频在线观看| 国产免费视频播放在线视频| 十八禁高潮呻吟视频 | 久久久精品94久久精品| 下体分泌物呈黄色| 18+在线观看网站| 制服丝袜香蕉在线| 国产亚洲午夜精品一区二区久久| 少妇人妻 视频| 97超视频在线观看视频| 老司机影院成人| 伊人久久国产一区二区| 国产精品.久久久| 日本午夜av视频| 国产综合精华液| 99热这里只有是精品在线观看| 欧美变态另类bdsm刘玥| 免费人妻精品一区二区三区视频| 精品少妇久久久久久888优播| 欧美丝袜亚洲另类| 欧美人与善性xxx| 国产黄频视频在线观看| 日本爱情动作片www.在线观看| 久久精品久久久久久久性| 97超碰精品成人国产| 男女啪啪激烈高潮av片| 亚洲情色 制服丝袜| 国产中年淑女户外野战色| 在线观看免费视频网站a站| 大话2 男鬼变身卡| 国产亚洲午夜精品一区二区久久| 韩国高清视频一区二区三区| 久久久久国产网址| 九九久久精品国产亚洲av麻豆| 久久国产亚洲av麻豆专区| 三级国产精品片| kizo精华| 亚洲真实伦在线观看| 中国国产av一级| 日韩亚洲欧美综合| 欧美精品高潮呻吟av久久| 99视频精品全部免费 在线| www.色视频.com| 在线观看免费高清a一片| 亚洲精品国产av蜜桃| 国产一区二区三区综合在线观看 | 三级国产精品欧美在线观看| 中文精品一卡2卡3卡4更新| 精品卡一卡二卡四卡免费| 免费大片18禁| 精品国产国语对白av| 亚洲精品色激情综合| 精品视频人人做人人爽| 中文字幕久久专区| 欧美激情极品国产一区二区三区 | 久久6这里有精品| 高清av免费在线| 免费观看a级毛片全部| 在线观看av片永久免费下载| 一区在线观看完整版| 国产一区亚洲一区在线观看| 久久97久久精品| 熟女人妻精品中文字幕| 国产乱人偷精品视频| 欧美人与善性xxx| 久久99一区二区三区| 五月伊人婷婷丁香| 一本色道久久久久久精品综合| 欧美日韩视频高清一区二区三区二| 99久国产av精品国产电影| 久久精品国产亚洲网站| 免费少妇av软件| 国产成人精品福利久久| 最新的欧美精品一区二区| 国产一区二区三区av在线| 免费观看性生交大片5| 妹子高潮喷水视频| 波野结衣二区三区在线| 一级片'在线观看视频| 欧美日韩av久久| 国产午夜精品一二区理论片| 欧美亚洲 丝袜 人妻 在线| 99九九线精品视频在线观看视频| 欧美日韩亚洲高清精品| 亚洲成色77777| 美女国产视频在线观看| a 毛片基地| 成人午夜精彩视频在线观看| 成人黄色视频免费在线看| 亚洲色图综合在线观看| 丁香六月天网| 日韩欧美 国产精品| 欧美少妇被猛烈插入视频| 国产精品一二三区在线看| 一个人看视频在线观看www免费| 五月伊人婷婷丁香| 热99国产精品久久久久久7| 亚洲av成人精品一区久久| 亚洲精品aⅴ在线观看| 激情五月婷婷亚洲| 国产91av在线免费观看| 亚洲电影在线观看av| 亚洲内射少妇av| 国产日韩欧美视频二区| 看免费成人av毛片| 97超视频在线观看视频| 国产日韩欧美视频二区| av福利片在线观看| 国国产精品蜜臀av免费| 久久久久久久国产电影| 亚洲精品aⅴ在线观看| 久久久久人妻精品一区果冻| 高清毛片免费看| 夜夜看夜夜爽夜夜摸| 精品久久久久久久久av| 晚上一个人看的免费电影| 久久久久久久久久久免费av| 亚洲av福利一区| 久久精品夜色国产| 女人久久www免费人成看片| av免费在线看不卡| 国产视频内射| 国产高清不卡午夜福利| 亚洲精品国产成人久久av| 在线观看www视频免费| 免费看光身美女| 国产一区二区在线观看av| 一本一本综合久久| 亚洲精华国产精华液的使用体验| 久久久久久久久大av| 国产淫片久久久久久久久| 在线观看免费视频网站a站| 少妇高潮的动态图| 亚洲,一卡二卡三卡| 精品久久久久久久久亚洲| 永久网站在线| 永久网站在线| 国产欧美日韩综合在线一区二区 | 夜夜骑夜夜射夜夜干| 欧美成人午夜免费资源| 亚洲久久久国产精品| 日韩人妻高清精品专区| 777米奇影视久久| 午夜老司机福利剧场| 中文字幕免费在线视频6| 欧美老熟妇乱子伦牲交| 久久热精品热| 国产精品免费大片| av不卡在线播放| 国产视频内射| 乱码一卡2卡4卡精品| 国产亚洲午夜精品一区二区久久| 亚洲精品久久午夜乱码| 免费黄色在线免费观看| 丝袜在线中文字幕| 免费看光身美女| 日韩欧美 国产精品| 日韩av免费高清视频| 亚洲,一卡二卡三卡| 最近的中文字幕免费完整| 精品99又大又爽又粗少妇毛片| 国产又色又爽无遮挡免| 边亲边吃奶的免费视频| 国产成人精品婷婷| 丰满少妇做爰视频| 亚洲精品456在线播放app| 成年人免费黄色播放视频 | av天堂中文字幕网| 精品久久久久久久久亚洲| 美女主播在线视频| 久久人人爽人人爽人人片va| 午夜91福利影院| 亚洲图色成人| 国产成人精品福利久久| 国产熟女午夜一区二区三区 | av天堂久久9| 观看美女的网站| 91精品伊人久久大香线蕉| 王馨瑶露胸无遮挡在线观看| 少妇 在线观看| 99热6这里只有精品| 久久久久久久大尺度免费视频| 久久久a久久爽久久v久久| 尾随美女入室| 制服丝袜香蕉在线| 精品一区二区免费观看| 午夜91福利影院| 极品少妇高潮喷水抽搐| 色网站视频免费| 精品一区二区免费观看| 亚洲成人手机| 乱人伦中国视频| 久久国产精品大桥未久av | 国产精品伦人一区二区| 午夜福利,免费看| 国产精品成人在线| 亚洲国产精品一区三区| 国产中年淑女户外野战色| 国产伦理片在线播放av一区| 久久99蜜桃精品久久| 视频中文字幕在线观看| 青春草国产在线视频| 久久久久久人妻| 国产精品久久久久久精品古装| a级毛色黄片| 久久久a久久爽久久v久久| 国产成人一区二区在线| 最近中文字幕2019免费版| 自拍偷自拍亚洲精品老妇| 噜噜噜噜噜久久久久久91| 亚洲精品第二区| 久久这里有精品视频免费| 永久免费av网站大全| 有码 亚洲区| 熟女人妻精品中文字幕| 午夜视频国产福利| 91午夜精品亚洲一区二区三区| 激情五月婷婷亚洲| 一级黄片播放器| 久久毛片免费看一区二区三区| 国产极品天堂在线| 成人18禁高潮啪啪吃奶动态图 | 精品卡一卡二卡四卡免费| 蜜桃在线观看..| av免费在线看不卡| 免费黄网站久久成人精品| 肉色欧美久久久久久久蜜桃| 国产亚洲午夜精品一区二区久久| 亚洲国产精品国产精品| 免费高清在线观看视频在线观看| 人人妻人人添人人爽欧美一区卜| 免费在线观看成人毛片| 晚上一个人看的免费电影| 在线观看国产h片| 日韩亚洲欧美综合| 久久精品熟女亚洲av麻豆精品| 制服丝袜香蕉在线| 亚洲精品国产成人久久av| 国产欧美日韩一区二区三区在线 | 久久精品国产a三级三级三级| 男的添女的下面高潮视频| 亚洲综合色惰| 亚洲av免费高清在线观看| av在线老鸭窝| 你懂的网址亚洲精品在线观看| 中文精品一卡2卡3卡4更新| 亚洲欧美中文字幕日韩二区| 九九爱精品视频在线观看| 99视频精品全部免费 在线| 精品国产乱码久久久久久小说| 亚洲av电影在线观看一区二区三区| 老司机亚洲免费影院| 2021少妇久久久久久久久久久| 80岁老熟妇乱子伦牲交| 自拍偷自拍亚洲精品老妇| 一区二区三区四区激情视频| 亚洲av在线观看美女高潮| 久久人妻熟女aⅴ| 丰满乱子伦码专区| 另类精品久久| 五月开心婷婷网| 各种免费的搞黄视频| 亚洲人与动物交配视频| 高清视频免费观看一区二区| 日本wwww免费看| 男人和女人高潮做爰伦理| 亚洲精品一二三| 美女主播在线视频| 9色porny在线观看| 国产视频内射| 2022亚洲国产成人精品| 男人和女人高潮做爰伦理| 日韩熟女老妇一区二区性免费视频| 丰满迷人的少妇在线观看| 日韩电影二区| 国产国拍精品亚洲av在线观看| 极品少妇高潮喷水抽搐| 亚洲欧洲国产日韩| 国产欧美日韩一区二区三区在线 | 亚洲欧洲日产国产| 久久精品熟女亚洲av麻豆精品| av免费在线看不卡| 精品亚洲成a人片在线观看| 久久精品国产自在天天线| 青春草国产在线视频| 欧美精品人与动牲交sv欧美| 成年人午夜在线观看视频| 美女脱内裤让男人舔精品视频| 啦啦啦视频在线资源免费观看| 少妇的逼水好多| 啦啦啦在线观看免费高清www| 又爽又黄a免费视频| 26uuu在线亚洲综合色| 婷婷色av中文字幕| 国产精品国产三级专区第一集| 亚洲国产精品成人久久小说| 两个人的视频大全免费| 国产精品久久久久久精品电影小说| 中文字幕人妻丝袜制服| 国产亚洲精品久久久com| 欧美最新免费一区二区三区| 99久久综合免费| 成年人免费黄色播放视频 | 日韩欧美一区视频在线观看 | 人妻一区二区av| √禁漫天堂资源中文www| 欧美精品高潮呻吟av久久| 久久久久久久久大av| 九九爱精品视频在线观看| 五月伊人婷婷丁香| 久久久久久久大尺度免费视频| 黄片无遮挡物在线观看| 免费看日本二区| 国产黄频视频在线观看| 婷婷色av中文字幕| 久久99蜜桃精品久久| 秋霞在线观看毛片| 在线观看免费日韩欧美大片 | 91精品一卡2卡3卡4卡| 丰满乱子伦码专区| kizo精华| 一级毛片久久久久久久久女| 一个人看视频在线观看www免费| 美女xxoo啪啪120秒动态图| 欧美高清成人免费视频www| 97在线视频观看| 日韩精品免费视频一区二区三区 | 久久免费观看电影| 久久人人爽人人片av| 午夜日本视频在线| 国产视频内射| 亚洲精品日韩av片在线观看| 国产精品伦人一区二区| 国产欧美日韩综合在线一区二区 | 亚洲精品视频女| 欧美精品亚洲一区二区| 纵有疾风起免费观看全集完整版| 国产精品99久久99久久久不卡 | 欧美国产精品一级二级三级 | 日日摸夜夜添夜夜爱| 黄色配什么色好看| 国产成人a∨麻豆精品| 大香蕉久久网| 99九九在线精品视频 | 精品人妻熟女av久视频| 午夜日本视频在线| 欧美丝袜亚洲另类| av福利片在线| 亚洲精品乱码久久久v下载方式| 亚洲欧洲国产日韩| 乱码一卡2卡4卡精品| 成人国产av品久久久| 在线精品无人区一区二区三| 观看美女的网站| 亚洲在久久综合| 久久人人爽av亚洲精品天堂| 国产成人freesex在线| 色视频在线一区二区三区| 99国产精品免费福利视频| 一级毛片久久久久久久久女| 免费大片黄手机在线观看| 91在线精品国自产拍蜜月| 国产在线视频一区二区| 毛片一级片免费看久久久久| 国产男女超爽视频在线观看| 插阴视频在线观看视频| 免费观看av网站的网址| 亚洲国产色片| 亚洲精品自拍成人| 国产69精品久久久久777片| 如日韩欧美国产精品一区二区三区 | 丰满人妻一区二区三区视频av| 91成人精品电影| 日韩三级伦理在线观看| 秋霞伦理黄片| 91久久精品国产一区二区三区| kizo精华| 久久97久久精品| av有码第一页| 日韩免费高清中文字幕av| 免费观看在线日韩| 国产精品一区www在线观看| 日韩一本色道免费dvd| 亚洲欧美日韩另类电影网站| 午夜av观看不卡| 中文字幕亚洲精品专区| 亚洲av.av天堂| 一个人免费看片子| 免费久久久久久久精品成人欧美视频 | 成人18禁高潮啪啪吃奶动态图 | 在线天堂最新版资源| 亚洲av免费高清在线观看| 天天操日日干夜夜撸| 丝袜脚勾引网站| 91午夜精品亚洲一区二区三区| 国产乱来视频区| 成人黄色视频免费在线看| 99热这里只有精品一区| 97超视频在线观看视频| 人妻 亚洲 视频| 久久午夜综合久久蜜桃| 岛国毛片在线播放| 最近手机中文字幕大全| 22中文网久久字幕| 亚洲成人一二三区av| 久久精品久久久久久噜噜老黄| 亚洲国产精品国产精品| 久久久国产一区二区| 午夜免费男女啪啪视频观看| 黄色怎么调成土黄色| 国产白丝娇喘喷水9色精品| 最后的刺客免费高清国语| 寂寞人妻少妇视频99o| 水蜜桃什么品种好| 一区二区三区精品91| 国产精品蜜桃在线观看| 桃花免费在线播放| 成年女人在线观看亚洲视频| 一级毛片久久久久久久久女| 久久女婷五月综合色啪小说| 男人和女人高潮做爰伦理| 男女边摸边吃奶| 国产在线一区二区三区精| 日韩中文字幕视频在线看片| 人人妻人人看人人澡| 久久久久久久精品精品| 国产黄片美女视频| 国精品久久久久久国模美| 欧美人与善性xxx| 国产伦精品一区二区三区四那| 精品人妻偷拍中文字幕| 啦啦啦视频在线资源免费观看| 亚洲无线观看免费| 久久久精品94久久精品| 老女人水多毛片| 国产真实伦视频高清在线观看| 国产中年淑女户外野战色| 久久人妻熟女aⅴ| 午夜福利视频精品| 久久av网站| 国产精品欧美亚洲77777| 国产在线一区二区三区精| 欧美区成人在线视频| 欧美国产精品一级二级三级 | 99视频精品全部免费 在线| 91在线精品国自产拍蜜月| 亚洲精华国产精华液的使用体验| 国产成人精品婷婷| 又爽又黄a免费视频| 国产极品粉嫩免费观看在线 | 日本欧美国产在线视频| 国产精品久久久久久av不卡| www.av在线官网国产| 国产在视频线精品| 亚洲精品乱码久久久久久按摩| 亚洲av综合色区一区| 日本vs欧美在线观看视频 | av国产精品久久久久影院| 人体艺术视频欧美日本| 日韩 亚洲 欧美在线| 哪个播放器可以免费观看大片| 新久久久久国产一级毛片| 有码 亚洲区| 伊人久久精品亚洲午夜| 欧美 亚洲 国产 日韩一| 男女啪啪激烈高潮av片| 精品99又大又爽又粗少妇毛片| 亚洲自偷自拍三级| 人人妻人人爽人人添夜夜欢视频 | 久久鲁丝午夜福利片| 最近手机中文字幕大全| av在线观看视频网站免费| 国产爽快片一区二区三区| 在线观看三级黄色| 久久国产精品男人的天堂亚洲 | 人妻 亚洲 视频| 亚洲国产欧美在线一区| 国产日韩欧美亚洲二区| 3wmmmm亚洲av在线观看| 建设人人有责人人尽责人人享有的| 日韩人妻高清精品专区| 我要看黄色一级片免费的| 亚洲国产日韩一区二区| 亚洲情色 制服丝袜| 精品一品国产午夜福利视频| 亚州av有码| 亚洲国产成人一精品久久久| 国产精品成人在线| 在线免费观看不下载黄p国产| 久久99一区二区三区| 热99国产精品久久久久久7| av.在线天堂| 春色校园在线视频观看| 99视频精品全部免费 在线| 热re99久久国产66热| 精品酒店卫生间| 国产免费视频播放在线视频| 人妻制服诱惑在线中文字幕| 色视频在线一区二区三区| 午夜福利,免费看| 亚洲情色 制服丝袜| 国产成人一区二区在线| 美女xxoo啪啪120秒动态图| 人妻一区二区av| 综合色丁香网| 精品少妇黑人巨大在线播放| 国产视频首页在线观看| 日韩一区二区视频免费看| 日本av手机在线免费观看| 精品国产露脸久久av麻豆| 人妻系列 视频| 国产成人免费观看mmmm| 日韩成人伦理影院| 亚洲va在线va天堂va国产| av在线观看视频网站免费| 中文字幕人妻熟人妻熟丝袜美| 噜噜噜噜噜久久久久久91| 久久久久久久久久人人人人人人| 久久精品国产亚洲网站| 十分钟在线观看高清视频www | 一级片'在线观看视频| 欧美日本中文国产一区发布| av女优亚洲男人天堂| 日韩伦理黄色片| 日韩精品免费视频一区二区三区 | 色婷婷av一区二区三区视频| 亚洲欧美成人精品一区二区| 午夜福利影视在线免费观看| 91aial.com中文字幕在线观看| 国产精品麻豆人妻色哟哟久久| 视频中文字幕在线观看| 我的老师免费观看完整版| 极品教师在线视频| 日本猛色少妇xxxxx猛交久久| 大香蕉久久网| 久久久久国产精品人妻一区二区| 极品人妻少妇av视频| 亚洲欧美成人精品一区二区| 亚洲欧美清纯卡通| 亚洲精品日韩av片在线观看| 亚洲丝袜综合中文字幕| 在线观看免费高清a一片| a级毛片在线看网站| 久久精品国产a三级三级三级| 狂野欧美激情性xxxx在线观看| 国产高清不卡午夜福利| 丰满饥渴人妻一区二区三| 日韩中文字幕视频在线看片| 国产欧美亚洲国产| 狂野欧美激情性bbbbbb| 亚洲综合色惰| 在线精品无人区一区二区三| 亚洲国产色片| 精品久久久久久电影网| 色5月婷婷丁香| 五月天丁香电影| 日韩精品有码人妻一区| 久久久欧美国产精品| 亚洲欧美清纯卡通| 69精品国产乱码久久久| 日韩欧美一区视频在线观看 | 亚洲精品乱码久久久久久按摩| 高清在线视频一区二区三区| 在线观看美女被高潮喷水网站| 高清黄色对白视频在线免费看 | 一级a做视频免费观看| 波野结衣二区三区在线| 色视频www国产| 少妇人妻一区二区三区视频| 精品一品国产午夜福利视频| 午夜91福利影院| 国产真实伦视频高清在线观看| 久久狼人影院| 一区二区三区免费毛片| 欧美精品一区二区免费开放| 成年女人在线观看亚洲视频| 精品国产国语对白av| 91aial.com中文字幕在线观看| 成年女人在线观看亚洲视频| 成人无遮挡网站| 欧美日韩视频精品一区| a级一级毛片免费在线观看| 久久午夜福利片| 深夜a级毛片| 日韩精品免费视频一区二区三区 | 亚洲自偷自拍三级| 性高湖久久久久久久久免费观看| 国产日韩欧美视频二区| 亚洲av福利一区| 夫妻性生交免费视频一级片| 欧美3d第一页| 夫妻午夜视频| av.在线天堂| 亚洲人成网站在线播| 成人18禁高潮啪啪吃奶动态图 | 69精品国产乱码久久久| 午夜免费男女啪啪视频观看| av不卡在线播放| 一区二区av电影网| 在线精品无人区一区二区三| 国产美女午夜福利| 女人精品久久久久毛片|