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

    環(huán)狀給水管網(wǎng)自動生成樹的研究

    2018-07-09 12:32:14趙星明
    中國農(nóng)村水利水電 2018年6期
    關(guān)鍵詞:鄰接矩陣管段給水管

    趙星明,王 萱

    (山東農(nóng)業(yè)大學(xué)水利土木工程學(xué)院,山東 泰安 271018)

    城鎮(zhèn)給水管網(wǎng)的投資占整個給水工程比例較大,因此其最優(yōu)化設(shè)計一直受到重視。最優(yōu)化設(shè)計的前提是對環(huán)狀給水管網(wǎng)的管段流量確定初分配方案,若給水管網(wǎng)的規(guī)模較小,可通過分析控制點位置以及大用戶的分布,確定主干管供水路線,較合理地分配管段流量。對于中等以上的城市,僅依靠主觀判斷分配給水管網(wǎng)的管段流量會造成較大誤差,不能滿足用戶的實際用水量,也不能實現(xiàn)最優(yōu)化設(shè)計。

    初始流量分配可以采用最小平方和的流量分配法和截面法等,最小平方和法分配的管段流量比較均勻而使管道主次不分,截面法分配的管段流量不能滿足節(jié)點流量平衡的條件。若采用圖論理論把環(huán)狀給水管網(wǎng)中的若干條管段刪除,生成樹后再進行初始流量分配,可符合節(jié)點流量連續(xù)性方程。把管網(wǎng)圖化為生成樹,須按照歐拉公式刪除與基環(huán)數(shù)一樣多的管段來破壞所有的環(huán),并保持原管網(wǎng)的連通性,對一個規(guī)模較大的環(huán)狀給水管網(wǎng),僅靠主觀判斷生成樹需花費很多時間。

    本文根據(jù)Kruskal算法對環(huán)狀管網(wǎng)生成樹問題進行研究,利用帶權(quán)重的鄰接矩陣得到了最小生成樹,并在AutoCAD平臺上運用深度優(yōu)化搜索方法,實現(xiàn)了刪除連枝自動生成樹,可方便管段初始流量的計算,提高環(huán)狀給水管網(wǎng)的設(shè)計計算效率和準(zhǔn)確度。同時,通過人工干預(yù)保留一些管段流量易確定的連枝,以優(yōu)化環(huán)狀給水管網(wǎng)布局。

    1 環(huán)狀給水管網(wǎng)生成樹的方法

    一個連通管網(wǎng)圖G(V,E),由N節(jié)點和M管段構(gòu)成,V和E為節(jié)點和管段的集合,根據(jù)圖論理論計算其內(nèi)環(huán)數(shù)L=M-N+1。若想把G(V,E)生成一棵樹需刪除與內(nèi)環(huán)數(shù)相等的L條管段,刪除的管段稱為連枝,保留下來的管段稱為樹枝。在環(huán)狀給水管網(wǎng)的工程設(shè)計中,連接管和末端管主要承擔(dān)本段沿線流量的分配,流量容易確定,可作為連枝刪除掉,一部分管段流量不易確定的主干管和轉(zhuǎn)輸管等作為樹枝保留下來,這些樹枝的管段流量可在生成樹管網(wǎng)圖上通過節(jié)點連續(xù)性方程求出。為了保留重要的主干管和轉(zhuǎn)輸管,把管網(wǎng)圖G(V,E)轉(zhuǎn)化為一個加權(quán)圖G(V,E,W),W為管段權(quán)重Wij的集合,在此假設(shè)想保留下來管道的Wij=1,其權(quán)重最小,而其他管道的Wij>1,管道的權(quán)重越大越容易被刪除。需注意的是權(quán)重大于1的管道比需刪除的管段數(shù)L要多,所以權(quán)重大于1的管道不會全部刪除,也要保留一部分,哪些管道保留下來或被刪除掉是根據(jù)Kruskal算法確定的。

    設(shè)T為G的一棵生成樹,管段eij的權(quán)為Wij,則T的權(quán)定義為:

    (1)

    所有生成樹中權(quán)重最小的是G的一棵最小生成樹,可以采用Kruskal等算法求得,剪枝環(huán)狀管網(wǎng)生成一棵最小樹。

    2 Kruskal算法生成樹的實現(xiàn)

    2.1 Kruskal算法的基本思想

    含有n個節(jié)點的管網(wǎng)連通圖G(V,E,W)和其生成樹T的節(jié)點V集合相同,定義生成樹T(V,E,W)的初態(tài)為空集,即生成樹只是由n個節(jié)點構(gòu)成而無任何邊的n個非連通圖T=(V,{φ},{φ}),因節(jié)點之間沒有管段,故每個節(jié)點自成1個連通分量。先對環(huán)狀管網(wǎng)圖G中的管段進行分析和分類,人工賦予不同的權(quán)重。對管網(wǎng)圖G中的所有管段進行第1次遍歷,若遍歷到的管段的權(quán)重最小,即為保留管段,則將此管段直接加入到最小生成樹T中的E集合。再對權(quán)重較大的管段進行第2次遍歷,如果管段的2個節(jié)點依附于T的2個不同的連通分量,則將此管段加入最小生成樹T中的E集合,同時把2個連通分量連接為1個連通分量;若被遍歷到的2個節(jié)點同屬1個連通分量,說明連接此管段后會造成回路,故不連接這2個節(jié)點。遍歷所有管段后,可使T只有1個連通分量,便構(gòu)建了1棵連通的生成樹,可用于管段流量的初分配。

    2.2 Kruskal算法的實現(xiàn)方法

    根據(jù)環(huán)狀給水管網(wǎng)圖構(gòu)建帶權(quán)的鄰接矩陣M,鄰接矩陣M由4個元素組成,其格式為:

    管段編號,管段的起點,終點,管段權(quán)值

    以圖1的環(huán)狀給水管網(wǎng)為例,其鄰接矩陣保存格式及數(shù)據(jù)如表1所示,管段1、2、3、10和12的權(quán)值最小,屬保留的管段,其他管段按照Kruskal算法刪除連枝以生成樹。

    圖1 環(huán)狀給水管網(wǎng)圖Fig.1 The looped pipe network of water supply

    在使用Kruskal算法對環(huán)狀連通管網(wǎng)圖生成樹的實現(xiàn)過程中,根據(jù)鄰接矩陣的格式,構(gòu)建一個加權(quán)的2維管段數(shù)組Looped(Edge, 4),分別儲存每個管段的編號、起點、終點和權(quán)值這4個元素,Edge為環(huán)狀管網(wǎng)的管段編號。還設(shè)置一個Root(VertexNum)數(shù)組,使用自定義函數(shù)Search搜索每個管段的起

    表1 鄰接矩陣格式及數(shù)據(jù)Tab.1 Fortmat and data of adjacency matrix

    點和終點所屬的連通分量,VertexNum為節(jié)點編號。Root(VertexNum)數(shù)組中每個節(jié)點VertexNum的初值為0,表示管網(wǎng)的各個節(jié)點在不同的連通分量上,遍歷到一條管段后,讀取管段的起點(StartVertex)和終點編號(EndVertex),用Search自定義函數(shù)搜索2個節(jié)點所在樹的根結(jié)點,也就是在Root(VertexNum)數(shù)組中的序號。若StartVertex和EndVertex的根結(jié)點不同,表示所對應(yīng)的管段不在同一連通分量上,則將這條管段存儲在生成樹Tree數(shù)組中,并把StartVertex的值賦給終點數(shù)組Root(EndVertex),合并2個節(jié)點所屬的連通分量為1個連通分量。

    鄰接矩陣的排列順序影響著生成樹過程中連枝的選擇,可以對鄰接矩陣的權(quán)值從小到大進行排序,僅遍歷一次便可構(gòu)建最小生成樹。因在環(huán)狀管網(wǎng)生成樹過程中,只有人工干預(yù)保留部分管段,所以管段的權(quán)值只有1和2,權(quán)值為1的管段優(yōu)先保留下來,權(quán)值為2的管段根據(jù)生成樹的算法再確定是否保留或刪除,為此可以對鄰接矩陣的權(quán)值不進行排序,只是對管段進行2次遍歷,第1次遍歷只把權(quán)值為1的管段直接存儲到Tree數(shù)組中,然后進行第2遍歷,把權(quán)值為2的管段根據(jù)Kruskal算法存儲到Tree數(shù)組。

    用Kruskal算法實現(xiàn)的過程中,首先把鄰接矩陣的元素以文本文件形式保存,程序讀取所有管段的編號、起點、終點和權(quán)值到Looped(Edge, 4) 數(shù)組,然后根據(jù)優(yōu)化搜索算法判斷每個管段是否在同一連通分量上,最后得到環(huán)狀管網(wǎng)圖的最小生成樹,以數(shù)組Tree(Edge, 4)的形式存儲并寫入交換格式文件中。

    對于圖1管網(wǎng)圖按照表1的鄰接矩陣,用Kruskal算法構(gòu)建最小生成樹T,則會完全保留權(quán)重為1的管段1、2、3、10和12,刪除權(quán)重為2的部分管段6、7、9、11。得到的生成樹T的管段集合M={1,2,3,10,12,4,5,8},節(jié)點集合V與管網(wǎng)圖G相同。若對得到的生成樹不是很滿意,可修改管段的權(quán)重,也可增加管段權(quán)重的級數(shù)為3級或更多,權(quán)值改為1、2和3,把最想刪除的管段權(quán)重設(shè)置為3,再通過程序生成最小樹。Kruskal算法的主程序代碼如下:

    Option Base 1

    Sub main()

    Dim Root() As Integer

    Dim i As Integer, j As Integer, k As Integer

    Dim StartVertex As Integer, EndVertex As Integer

    Dim VertexNum As Integer, Edge As Integer

    Dim Looped() As Integer, Tree() As Integer

    Open ActiveDocument.Path & "in.csv" For Input As #1

    Input #1, VertexNum, Edge ‘讀取管網(wǎng)圖的節(jié)點總數(shù)和管段總數(shù)

    ReDim Root(VertexNum)

    ReDim Looped(Edge, 4)

    ReDim Tree(Edge, 4)

    For i = 1 To Edge

    Input #1, Looped(i, 1), Looped(i, 2), Looped(i, 3), Looped(i, 4)

    Next i

    Close #1

    For i = 1 To VertexNum

    Root(i) = 0

    Next i

    Open ActiveDocument.Path & "out.csv" For Output As #2

    Write #2, "管段編號", "起始節(jié)點", "終止節(jié)點", "管段權(quán)重"

    For k = 1 To 2

    For i = 1 To Edge

    StartVertex = Search(Root, Looped(i, 2))

    EndVertex = Search(Root, Looped(i, 3))

    If (StartVertex = EndVertex) And Looped(i, 4) = 2 Then

    Debug.Print Looped(i, 1) '刪除的連枝

    End If

    If (StartVertex <> EndVertex) And Looped(i, 4) = k Then

    Root(EndVertex) = StartVertex

    For j = 1 To 4

    Tree(i, j) = Looped(i, j)

    Write #2, Tree(i, j),

    Next j

    Print #2,

    End If

    Next i, k

    Close #2

    End Sub

    自定義函數(shù)Search的代碼:

    Function Search(Root() As Integer, VertexNum As Integer) As Integer

    Dim i As Integer

    i = VertexNum

    Do While Root(i) > 0

    i = Root(i)

    Loop

    Search = i

    End Function

    3 在AutoCAD平臺上搜索環(huán)狀管網(wǎng)生成樹的實現(xiàn)及算例研究

    3.1 工程條件

    某市區(qū)供水管網(wǎng)如圖2所示。在AutoCAD平臺上,假如所有管道的圖層為“鑄鐵管”,可采用遍歷技術(shù)對管段和節(jié)點進行自動編號,形成不帶權(quán)的3元素鄰接矩陣,以.csv交換格式文件形式保存[1,5]。對需保留的管道,手動選擇并改變其圖層為“Tree”。為表示管道的不同權(quán)重,設(shè)定“Tree”圖層的線寬為0.35 mm,定義在該圖層的管道權(quán)值為1,設(shè)定“鑄鐵管”圖層的線寬為0.2 mm,表示管道權(quán)值為2。用Excel修改3元素鄰接矩陣,增加每條管道的權(quán)重元素,形成帶權(quán)的4元素關(guān)聯(lián)矩陣,即其矩陣元素仍按表1的格式存儲,數(shù)據(jù)內(nèi)容見表2。

    表2 某市區(qū)供水管網(wǎng)鄰接矩陣Tab.2 The adjacency matrix of water supply pipe network

    圖2 某市區(qū)供水管網(wǎng)現(xiàn)狀圖Fig.2 The water supply pipe network

    3.2 在AutoCAD平臺上生成樹

    首先使用前面的Kruskal算法程序,讀取某市區(qū)供水管網(wǎng)的鄰接矩陣,即表2的數(shù)據(jù),判斷管段的2個節(jié)點的根結(jié)點是否相同,若相同則為應(yīng)刪除的連枝。根據(jù)歐拉公式,本算例有18個環(huán),所以產(chǎn)生的連枝數(shù)也為18,具體刪除的管段見表3。

    表3 刪除的管段Tab.3 The deleted pipe section

    然后在AutoCAD平臺上遍歷搜索環(huán)狀給水管網(wǎng)的所有管段,用選擇集的方法識別管段編號,若遍歷的管段為表3中應(yīng)刪除的管段,則修改該管段到“0”圖層,也可直接刪除,非常直觀地實現(xiàn)了把環(huán)狀管網(wǎng)自動生成樹,生成的樹狀管網(wǎng)如圖3所示,程序代碼如下:

    Sub 刪除連枝()

    Dim Linea As AcadLine

    Dim Obj As AcadEntity '遍歷圖形對象用

    Dim SSetColl As AcadSelectionSets '定義選擇集集合

    Dim ssetObj As AcadSelectionSet '選擇集對象

    Dim name As String '選擇集名稱

    Dim Textobj As AcadText '寫直線編號用的文字對象

    Dim Sp As Variant, Ep As Variant '直線兩端點

    Dim Cp As Variant '直線中點坐標(biāo),要計算

    Dim Lp As Variant, Rp As Variant '為形成一個選擇框,定義左上角和右下角坐標(biāo)

    Dim TextHeight As Single '定義管段編號的文字高度

    ReDim Cp(0 To 2) As Double

    ReDim Lp(0 To 2) As Double

    ReDim Rp(0 To 2) As Double

    Dim gpCode(0) As Integer

    Dim dataValue(0) As Variant

    Dim groupCode As Variant, dataCode As Variant

    gpCode(0) = 0

    dataValue(0) = "TEXT"

    圖3 某市區(qū)供水管網(wǎng)自動生成樹Fig.3 Automatic spanning tree of the water supply pipe network

    groupCode = gpCode

    dataCode = dataValue

    TextHeight = 100 '管段編號和節(jié)點編號的文字高度

    name = "SSET" '選擇集名字

    Set SSetColl = ThisDrawing.SelectionSets

    On Error Resume Next

    For Each Obj In ThisDrawing.ModelSpace

    If Obj.ObjectName = "AcDbLine" Then

    Set Linea = Obj

    Sp = Linea.Startpoint: Ep = Linea.Endpoint

    Cp(0) = (Sp(0) + Ep(0)) / 2: Cp(1) = (Sp(1) + Ep(1)) / 2: Cp(2) = (Sp(2) + Ep(2)) / 2 '直線中點,即管段編號插入點

    '以下為選擇集定義容錯處理的代碼

    If Not IsNull(SSetColl.Item(name)) Then

    Set ssetObj = SSetColl.Item(name)

    ssetObj.Delete

    End If

    Set ssetObj = SSetColl.Add(name)

    Lp(0) = Cp(0) - TextHeight: Lp(1) = Cp(1) + TextHeight: Lp(2) = Cp(2) '在直線中點形成一個矩形選擇框

    Rp(0) = Cp(0) + TextHeight: Rp(1) = Cp(1) - TextHeight: Rp(2) = Cp(2)

    ssetObj.Select acSelectionSetWindow, Lp, Rp, groupCode, dataCode

    If ssetObj.Count = 1 Then

    Set Textobj = ssetObj.Item(0)

    If InStr("[12][16][17][21][29][31][34][38][45][46][50][56][58][60][61][62][63][65]", Textobj.TextString) Then

    Debug.Print Textobj.TextString

    Linea.Layer = "0" '修改連枝到“0”圖層,也可直接刪除

    Textobj.Layer = "0"

    End If

    Else

    MsgBox "節(jié)點編號錯誤" '可能沒有編號,或者編號字體太大或者矩形選取框太小等原因

    Exit For

    End If

    End If

    Next Obj

    End Sub

    4 結(jié) 語

    本文根據(jù)Kruskal最小生成樹算法的基本思想,對環(huán)狀給水管網(wǎng)的鄰接矩陣進行深度優(yōu)化搜索,篩選出連枝而實現(xiàn)了自動生成樹的目的。在AutoCAD平臺上遍歷所有管段可自動生成帶權(quán)重的鄰接矩陣,根據(jù)Kruskal算法篩選出應(yīng)刪除的管段,使環(huán)狀管網(wǎng)直觀地生成樹。程序能夠根據(jù)每個管段的權(quán)重,識別人為保留的管段,生成的樹得到優(yōu)化,有助于提高大規(guī)模環(huán)狀管網(wǎng)的設(shè)計效率。

    在節(jié)點流量和連枝流量已知的情況下,滿足節(jié)點流量平衡的條件,利用Kruskal算法自動生成樹,使管段流量初分配很容易實現(xiàn)。

    參考文獻:

    [1] 趙星明,王 萱.環(huán)狀給水管網(wǎng)關(guān)聯(lián)矩陣的建立[J].中國農(nóng)村水利水電,2012,(11):129-131,135.

    [2] 嚴煦世,劉遂慶.給水排水管網(wǎng)系統(tǒng)[M]. 3版. 北京:中國建筑工業(yè)出版社,2014:79-80.

    [3] 唐策善,李龍澍,黃劉生.數(shù)據(jù)結(jié)構(gòu)[M]. 北京:高等教育出版社,2000:123-135.

    [4] 盧有杰, 吳煒煜.C語言高級程序設(shè)計[M]. 北京:清華大學(xué)出版社, 1991:19-86.

    [5] 趙星明,王 萱.基于擴展數(shù)據(jù)的給排水管網(wǎng)拓撲關(guān)系的構(gòu)建[J].山東農(nóng)業(yè)大學(xué)學(xué)報(自然科學(xué)版),2012,43(4):549-554.

    [6] 宋 芹,趙星明,艾典勝.輸水管道水錘分析與防護技術(shù)[J].山東農(nóng)業(yè)大學(xué)學(xué)報(自然科學(xué)版),2017,48(1):84-87.

    [7] SeijiKataoka,TakeoYamada.Algorithms for the minimum spanning tree problem with resource allocation[J].Operations Research Perspectives, 2016,(3):5-13.

    猜你喜歡
    鄰接矩陣管段給水管
    輪圖的平衡性
    高溫氣冷堆核電站蒸汽發(fā)生器可拆管段拆裝系統(tǒng)研究
    管段沿線流量簡化前后水頭和流行時間差異性分析
    ◆ 塑料管
    電站配管設(shè)計中的旋轉(zhuǎn)角度分析及計算
    美國PVC給水管材的質(zhì)量控制
    中國塑料(2016年8期)2016-06-27 06:35:08
    基于鄰接矩陣變型的K分網(wǎng)絡(luò)社團算法
    PPR管在給水管道施工中的應(yīng)用
    一種判定的無向圖連通性的快速Warshall算法
    某亞臨界鍋爐水冷壁泄漏失效分析
    機械工程師(2015年8期)2015-04-10 09:52:14
    欧美激情久久久久久爽电影| 国产精品一区二区三区四区免费观看 | 波多野结衣高清无吗| 91字幕亚洲| 精品人妻偷拍中文字幕| 欧美bdsm另类| a级一级毛片免费在线观看| 夜夜躁狠狠躁天天躁| 在线观看免费视频日本深夜| 深夜a级毛片| 综合色av麻豆| 久久亚洲真实| 久久久久精品国产欧美久久久| 欧美精品啪啪一区二区三区| 午夜a级毛片| 一夜夜www| 亚洲成a人片在线一区二区| 成人欧美大片| 亚洲成人免费电影在线观看| 宅男免费午夜| 一本一本综合久久| 精品久久久久久久人妻蜜臀av| 国产一区二区三区视频了| 国产蜜桃级精品一区二区三区| 两个人视频免费观看高清| 成年免费大片在线观看| 亚洲,欧美,日韩| 亚洲真实伦在线观看| 成人特级av手机在线观看| 亚洲欧美日韩高清专用| 99久久精品一区二区三区| 国产亚洲精品av在线| 日韩精品青青久久久久久| 九色成人免费人妻av| 桃红色精品国产亚洲av| 人妻丰满熟妇av一区二区三区| 中文亚洲av片在线观看爽| 免费av毛片视频| 亚洲精华国产精华精| 欧美极品一区二区三区四区| 亚洲精品一卡2卡三卡4卡5卡| 男人舔女人下体高潮全视频| 日本 欧美在线| 欧美高清性xxxxhd video| 日韩精品中文字幕看吧| 91av网一区二区| 性色av乱码一区二区三区2| 欧美丝袜亚洲另类 | 欧美一区二区精品小视频在线| 男女床上黄色一级片免费看| 久99久视频精品免费| 丁香六月欧美| 此物有八面人人有两片| av在线天堂中文字幕| 国产中年淑女户外野战色| 成人精品一区二区免费| 午夜福利高清视频| 中文字幕熟女人妻在线| 91麻豆av在线| 女同久久另类99精品国产91| 欧美精品国产亚洲| 一级黄片播放器| 亚洲中文日韩欧美视频| 亚洲乱码一区二区免费版| 俄罗斯特黄特色一大片| 99热这里只有是精品在线观看 | 午夜福利在线在线| 两人在一起打扑克的视频| 国产成人av教育| 国产欧美日韩一区二区精品| 精品欧美国产一区二区三| 国内精品久久久久久久电影| 黄色丝袜av网址大全| 国产探花在线观看一区二区| 婷婷精品国产亚洲av| 久久精品夜夜夜夜夜久久蜜豆| 免费电影在线观看免费观看| 99国产综合亚洲精品| 嫁个100分男人电影在线观看| 国产亚洲欧美98| av在线老鸭窝| 国产探花极品一区二区| 亚洲中文字幕日韩| 最好的美女福利视频网| 日韩欧美三级三区| 日本a在线网址| 好男人在线观看高清免费视频| 一级av片app| 精品人妻一区二区三区麻豆 | 国产精品免费一区二区三区在线| 国产精品久久久久久人妻精品电影| 女人被狂操c到高潮| 脱女人内裤的视频| 国产老妇女一区| 国产精品久久久久久亚洲av鲁大| 床上黄色一级片| 亚洲黑人精品在线| 精品午夜福利视频在线观看一区| 韩国av一区二区三区四区| 身体一侧抽搐| 五月伊人婷婷丁香| 宅男免费午夜| 99国产精品一区二区蜜桃av| 免费在线观看亚洲国产| 啦啦啦韩国在线观看视频| 国产熟女xx| 小说图片视频综合网站| 香蕉av资源在线| 色精品久久人妻99蜜桃| 国产亚洲欧美98| 国产精品一区二区免费欧美| 国产视频一区二区在线看| 日韩成人在线观看一区二区三区| 久久精品国产清高在天天线| 亚洲av熟女| 在线看三级毛片| 99热精品在线国产| 成年版毛片免费区| 亚洲 欧美 日韩 在线 免费| 中文字幕人成人乱码亚洲影| 久久久久久国产a免费观看| 免费人成视频x8x8入口观看| 欧美中文日本在线观看视频| 欧美日本视频| 成年女人看的毛片在线观看| 十八禁国产超污无遮挡网站| 欧美精品啪啪一区二区三区| 久久人人精品亚洲av| avwww免费| 免费高清视频大片| 色综合站精品国产| 狠狠狠狠99中文字幕| 久9热在线精品视频| 欧美国产日韩亚洲一区| 脱女人内裤的视频| 亚洲国产精品合色在线| 99久久久亚洲精品蜜臀av| 国产不卡一卡二| 亚洲性夜色夜夜综合| 最后的刺客免费高清国语| 最新在线观看一区二区三区| 国产色爽女视频免费观看| 国内少妇人妻偷人精品xxx网站| 亚洲性夜色夜夜综合| 亚洲av不卡在线观看| 99视频精品全部免费 在线| 国产av不卡久久| 最新在线观看一区二区三区| 黄色日韩在线| 国产精品爽爽va在线观看网站| 如何舔出高潮| 免费在线观看影片大全网站| 俄罗斯特黄特色一大片| 又黄又爽又刺激的免费视频.| 人妻夜夜爽99麻豆av| 欧美性猛交╳xxx乱大交人| 老司机福利观看| 午夜两性在线视频| 午夜福利免费观看在线| 级片在线观看| 成人国产一区最新在线观看| 久久午夜福利片| 夜夜躁狠狠躁天天躁| 国产三级在线视频| 国产亚洲av嫩草精品影院| 色综合婷婷激情| 少妇丰满av| 国产精品久久视频播放| 啦啦啦韩国在线观看视频| 久久国产乱子伦精品免费另类| 欧美成人a在线观看| 欧美最黄视频在线播放免费| 久久这里只有精品中国| 老司机午夜十八禁免费视频| 波野结衣二区三区在线| 国产精品精品国产色婷婷| 国产成人啪精品午夜网站| 99久久久亚洲精品蜜臀av| 国产伦人伦偷精品视频| 成人欧美大片| 久久久久久久久大av| 狂野欧美白嫩少妇大欣赏| 亚洲国产精品合色在线| 国产成人福利小说| 成人亚洲精品av一区二区| 成人特级黄色片久久久久久久| 亚洲av日韩精品久久久久久密| 日本黄色片子视频| 每晚都被弄得嗷嗷叫到高潮| 观看美女的网站| av黄色大香蕉| 亚洲精品成人久久久久久| 毛片女人毛片| 少妇高潮的动态图| 我要看日韩黄色一级片| 成人美女网站在线观看视频| 成人三级黄色视频| 一本久久中文字幕| 搡老熟女国产l中国老女人| 草草在线视频免费看| 可以在线观看的亚洲视频| 又黄又爽又免费观看的视频| av欧美777| 久久午夜福利片| 看黄色毛片网站| 一a级毛片在线观看| 免费观看的影片在线观看| 日本在线视频免费播放| 我的老师免费观看完整版| 色哟哟·www| 毛片女人毛片| 伦理电影大哥的女人| 九色成人免费人妻av| 一卡2卡三卡四卡精品乱码亚洲| 日本黄大片高清| 日本免费a在线| 精品久久久久久久久久免费视频| 精品国产三级普通话版| 三级毛片av免费| 国产精品1区2区在线观看.| 亚洲av日韩精品久久久久久密| 亚洲不卡免费看| 国产精品综合久久久久久久免费| а√天堂www在线а√下载| 淫妇啪啪啪对白视频| 日韩欧美精品v在线| 午夜免费激情av| 91狼人影院| 黄片小视频在线播放| 亚洲成a人片在线一区二区| 日韩中文字幕欧美一区二区| 91久久精品国产一区二区成人| 国产伦精品一区二区三区四那| av欧美777| 国产中年淑女户外野战色| 亚洲人成网站在线播放欧美日韩| av视频在线观看入口| a级一级毛片免费在线观看| 国产精品av视频在线免费观看| 51国产日韩欧美| 免费av不卡在线播放| 亚洲最大成人av| 乱码一卡2卡4卡精品| 村上凉子中文字幕在线| 麻豆成人午夜福利视频| 亚洲经典国产精华液单 | 欧美绝顶高潮抽搐喷水| 中文字幕熟女人妻在线| 久久国产精品人妻蜜桃| 国产主播在线观看一区二区| 99热6这里只有精品| 亚洲人成网站在线播放欧美日韩| 亚洲av熟女| 色综合婷婷激情| 欧美黄色片欧美黄色片| 国产久久久一区二区三区| 男人舔奶头视频| 午夜福利视频1000在线观看| h日本视频在线播放| 琪琪午夜伦伦电影理论片6080| 蜜桃久久精品国产亚洲av| 精品人妻1区二区| 青草久久国产| 国产老妇女一区| 精品免费久久久久久久清纯| 亚洲无线观看免费| 国产一区二区亚洲精品在线观看| 高清在线国产一区| 3wmmmm亚洲av在线观看| 美女 人体艺术 gogo| 亚洲成人免费电影在线观看| 久久天躁狠狠躁夜夜2o2o| 欧美日本视频| 亚洲av电影在线进入| 亚洲av成人精品一区久久| 午夜a级毛片| 成年免费大片在线观看| 亚洲熟妇中文字幕五十中出| 丰满人妻熟妇乱又伦精品不卡| 日韩成人在线观看一区二区三区| 精品熟女少妇八av免费久了| 国产精品亚洲一级av第二区| 91在线精品国自产拍蜜月| 国产激情偷乱视频一区二区| 精品乱码久久久久久99久播| 日韩欧美国产在线观看| 精品一区二区免费观看| 日韩欧美在线二视频| 自拍偷自拍亚洲精品老妇| 亚洲欧美日韩高清专用| 不卡一级毛片| 日本一二三区视频观看| 精品不卡国产一区二区三区| 99久久无色码亚洲精品果冻| 欧美日韩瑟瑟在线播放| 69人妻影院| 少妇的逼水好多| 欧美黑人欧美精品刺激| 十八禁国产超污无遮挡网站| 亚洲美女视频黄频| 亚洲久久久久久中文字幕| 色在线成人网| 狂野欧美白嫩少妇大欣赏| 亚洲精品在线美女| 国产亚洲av嫩草精品影院| .国产精品久久| 亚洲七黄色美女视频| 欧美日本视频| 亚洲性夜色夜夜综合| 美女被艹到高潮喷水动态| 久久中文看片网| 欧美日韩亚洲国产一区二区在线观看| 精品久久久久久久久亚洲 | 色5月婷婷丁香| 国产精品日韩av在线免费观看| 国产成人啪精品午夜网站| 每晚都被弄得嗷嗷叫到高潮| 成年女人永久免费观看视频| 国产成人欧美在线观看| 国产激情偷乱视频一区二区| 久久国产乱子伦精品免费另类| 性插视频无遮挡在线免费观看| 国产成人a区在线观看| 亚洲一区二区三区色噜噜| a级毛片免费高清观看在线播放| 老司机午夜十八禁免费视频| 亚洲经典国产精华液单 | 免费观看精品视频网站| 国产午夜精品久久久久久一区二区三区 | 久久午夜福利片| 免费人成在线观看视频色| 亚洲成人精品中文字幕电影| 一个人观看的视频www高清免费观看| av在线老鸭窝| 国产男靠女视频免费网站| 一区二区三区免费毛片| 国产精品野战在线观看| 亚洲精品色激情综合| 级片在线观看| 亚洲,欧美精品.| 日韩欧美国产在线观看| 日韩欧美三级三区| 天堂网av新在线| 国产精品免费一区二区三区在线| 十八禁国产超污无遮挡网站| 天堂动漫精品| 欧美最黄视频在线播放免费| 女生性感内裤真人,穿戴方法视频| 啦啦啦观看免费观看视频高清| 亚洲av免费在线观看| 级片在线观看| 久久久久久久久久黄片| 啦啦啦观看免费观看视频高清| 欧美黄色片欧美黄色片| 国产色婷婷99| 亚洲国产欧美人成| 99久久99久久久精品蜜桃| 久久久久久国产a免费观看| 欧美潮喷喷水| 亚洲无线在线观看| 欧美高清成人免费视频www| 90打野战视频偷拍视频| 国产精品不卡视频一区二区 | 白带黄色成豆腐渣| 毛片一级片免费看久久久久 | 99久久久亚洲精品蜜臀av| 真人一进一出gif抽搐免费| 一区二区三区激情视频| 欧美成狂野欧美在线观看| 国产精品98久久久久久宅男小说| 18美女黄网站色大片免费观看| 成人高潮视频无遮挡免费网站| 精品午夜福利在线看| 欧美日韩亚洲国产一区二区在线观看| 国产伦精品一区二区三区视频9| 色尼玛亚洲综合影院| 少妇丰满av| 国产精品久久电影中文字幕| 91麻豆精品激情在线观看国产| 国产aⅴ精品一区二区三区波| 91麻豆精品激情在线观看国产| 天堂影院成人在线观看| 深夜精品福利| 国产探花在线观看一区二区| or卡值多少钱| 国产精品乱码一区二三区的特点| 亚洲精华国产精华精| 在线a可以看的网站| 一边摸一边抽搐一进一小说| 亚洲 欧美 日韩 在线 免费| 精品熟女少妇八av免费久了| 国产av在哪里看| 夜夜夜夜夜久久久久| 人人妻人人澡欧美一区二区| 亚洲第一欧美日韩一区二区三区| 精品国产三级普通话版| 国产色爽女视频免费观看| 日本精品一区二区三区蜜桃| 国产伦精品一区二区三区四那| 亚洲无线在线观看| 一个人看视频在线观看www免费| 欧美成人性av电影在线观看| .国产精品久久| 嫁个100分男人电影在线观看| 永久网站在线| 最近在线观看免费完整版| 久久久久久九九精品二区国产| 能在线免费观看的黄片| 国产野战对白在线观看| 99热这里只有精品一区| 亚洲,欧美精品.| 亚洲三级黄色毛片| 国产精品自产拍在线观看55亚洲| 日本熟妇午夜| 性色avwww在线观看| 国产一级毛片七仙女欲春2| 成人av在线播放网站| 此物有八面人人有两片| 婷婷丁香在线五月| 他把我摸到了高潮在线观看| 亚洲三级黄色毛片| 久久久久久大精品| 日韩国内少妇激情av| 男女视频在线观看网站免费| 精品一区二区三区视频在线观看免费| 国产成人aa在线观看| 日本五十路高清| 如何舔出高潮| 美女黄网站色视频| 国产在线男女| 国产黄色小视频在线观看| 亚洲av五月六月丁香网| 亚洲成av人片免费观看| 免费在线观看影片大全网站| 色视频www国产| 少妇被粗大猛烈的视频| 男女那种视频在线观看| 色综合站精品国产| 欧美精品啪啪一区二区三区| 嫁个100分男人电影在线观看| 1000部很黄的大片| 国产69精品久久久久777片| 午夜福利欧美成人| 亚洲人成网站高清观看| 中文字幕熟女人妻在线| 白带黄色成豆腐渣| 国产精品嫩草影院av在线观看 | 午夜a级毛片| 最近最新中文字幕大全电影3| 九色成人免费人妻av| 国产熟女xx| 啪啪无遮挡十八禁网站| 欧美三级亚洲精品| 在线看三级毛片| 午夜福利欧美成人| 亚洲在线自拍视频| 亚洲最大成人av| 日韩精品中文字幕看吧| 欧美成狂野欧美在线观看| 又粗又爽又猛毛片免费看| 一级av片app| 欧美午夜高清在线| 小蜜桃在线观看免费完整版高清| 亚洲av二区三区四区| 日日干狠狠操夜夜爽| 99视频精品全部免费 在线| 国产人妻一区二区三区在| 国产野战对白在线观看| 99久久成人亚洲精品观看| 他把我摸到了高潮在线观看| 亚洲精品乱码久久久v下载方式| 好看av亚洲va欧美ⅴa在| av福利片在线观看| 深爱激情五月婷婷| 真人做人爱边吃奶动态| 国产精品一及| 久久伊人香网站| 男人舔女人下体高潮全视频| 国产精品伦人一区二区| 国产野战对白在线观看| 日本撒尿小便嘘嘘汇集6| 此物有八面人人有两片| 在线免费观看不下载黄p国产 | 成年女人看的毛片在线观看| 久久欧美精品欧美久久欧美| 亚洲欧美日韩卡通动漫| 99国产精品一区二区三区| 一进一出好大好爽视频| 一级作爱视频免费观看| 国产精品一区二区三区四区免费观看 | 欧美3d第一页| 午夜影院日韩av| 国产亚洲欧美98| 国产白丝娇喘喷水9色精品| 国产91精品成人一区二区三区| 天天躁日日操中文字幕| 亚洲经典国产精华液单 | 99热精品在线国产| 亚州av有码| 99热这里只有是精品50| 香蕉av资源在线| 亚洲天堂国产精品一区在线| 欧美午夜高清在线| av黄色大香蕉| 亚洲,欧美精品.| 亚洲av成人精品一区久久| 欧美黄色片欧美黄色片| 国产精品一区二区免费欧美| 在线免费观看的www视频| 亚洲天堂国产精品一区在线| 色吧在线观看| 国产又黄又爽又无遮挡在线| 我的女老师完整版在线观看| 级片在线观看| 欧美黄色片欧美黄色片| 天堂网av新在线| 99在线人妻在线中文字幕| 精品国内亚洲2022精品成人| 有码 亚洲区| 小蜜桃在线观看免费完整版高清| or卡值多少钱| 欧美+日韩+精品| 日本a在线网址| .国产精品久久| 在线免费观看不下载黄p国产 | 一边摸一边抽搐一进一小说| 神马国产精品三级电影在线观看| 在线播放国产精品三级| 90打野战视频偷拍视频| 精品不卡国产一区二区三区| 夜夜夜夜夜久久久久| 宅男免费午夜| 欧美乱色亚洲激情| 日韩亚洲欧美综合| 3wmmmm亚洲av在线观看| 国产欧美日韩精品一区二区| 老鸭窝网址在线观看| 91字幕亚洲| 日本在线视频免费播放| 我要看日韩黄色一级片| 夜夜看夜夜爽夜夜摸| 欧美3d第一页| 免费高清视频大片| 激情在线观看视频在线高清| 51国产日韩欧美| 日韩欧美一区二区三区在线观看| 欧美高清成人免费视频www| 性色avwww在线观看| 亚洲乱码一区二区免费版| 日本一本二区三区精品| 一夜夜www| 日韩欧美 国产精品| 内地一区二区视频在线| 亚洲av日韩精品久久久久久密| 欧美最新免费一区二区三区 | 国产高潮美女av| 蜜桃久久精品国产亚洲av| 午夜福利免费观看在线| 亚洲成人免费电影在线观看| 欧美中文日本在线观看视频| 久久精品影院6| 日韩精品中文字幕看吧| 三级毛片av免费| 国产私拍福利视频在线观看| 欧美绝顶高潮抽搐喷水| 简卡轻食公司| 他把我摸到了高潮在线观看| 亚洲,欧美,日韩| 99在线视频只有这里精品首页| 午夜a级毛片| 国产亚洲精品综合一区在线观看| 一本久久中文字幕| 美女免费视频网站| 精品久久久久久久人妻蜜臀av| 性插视频无遮挡在线免费观看| 99国产极品粉嫩在线观看| 国产高清三级在线| 欧美成人性av电影在线观看| 亚洲 国产 在线| 国产aⅴ精品一区二区三区波| 国产黄色小视频在线观看| 老熟妇仑乱视频hdxx| 国产成人啪精品午夜网站| 国产精品自产拍在线观看55亚洲| 永久网站在线| 亚洲美女搞黄在线观看 | 天堂√8在线中文| 国产乱人伦免费视频| 99热这里只有是精品在线观看 | 内射极品少妇av片p| 国产真实乱freesex| 免费av毛片视频| 国产免费男女视频| 亚洲人成电影免费在线| 91麻豆精品激情在线观看国产| 内射极品少妇av片p| 国产精品美女特级片免费视频播放器| 亚洲av.av天堂| 搞女人的毛片| 18禁黄网站禁片午夜丰满| 色吧在线观看| 黄色一级大片看看| 日本撒尿小便嘘嘘汇集6| 欧美精品国产亚洲| av福利片在线观看| 成年女人毛片免费观看观看9| 成年女人永久免费观看视频| 中文亚洲av片在线观看爽| 午夜福利成人在线免费观看| a级毛片免费高清观看在线播放| 亚洲电影在线观看av| 内地一区二区视频在线| 我要搜黄色片| 久久精品人妻少妇| 精品一区二区三区av网在线观看|