• <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
    午夜福利视频精品| tube8黄色片| 在线观看国产h片| www.熟女人妻精品国产| 国产亚洲午夜精品一区二区久久| 日本vs欧美在线观看视频| 99热国产这里只有精品6| 欧美人与性动交α欧美精品济南到 | 天堂俺去俺来也www色官网| 三级国产精品片| 九色亚洲精品在线播放| 日韩一区二区三区影片| 在线免费观看不下载黄p国产| 久久青草综合色| 国产男女内射视频| 亚洲一级一片aⅴ在线观看| 国产精品国产三级专区第一集| 久久久久国产网址| 乱人伦中国视频| 久久久久久免费高清国产稀缺| 丝袜在线中文字幕| 中文字幕制服av| 在线亚洲精品国产二区图片欧美| 美女视频免费永久观看网站| 麻豆精品久久久久久蜜桃| 欧美日韩av久久| 99热全是精品| 一区二区三区激情视频| 另类精品久久| 在线亚洲精品国产二区图片欧美| 999久久久国产精品视频| 国产视频首页在线观看| 国产av精品麻豆| 最近2019中文字幕mv第一页| 国产精品嫩草影院av在线观看| 久久影院123| 国产97色在线日韩免费| 久久精品久久久久久噜噜老黄| 国产精品无大码| 亚洲成人av在线免费| 九草在线视频观看| 午夜福利一区二区在线看| 看十八女毛片水多多多| 国产亚洲午夜精品一区二区久久| 黑丝袜美女国产一区| 国产又色又爽无遮挡免| 日韩一区二区三区影片| 另类精品久久| av在线app专区| 18禁观看日本| 老熟女久久久| 欧美亚洲日本最大视频资源| 一本久久精品| 欧美另类一区| 街头女战士在线观看网站| 青春草视频在线免费观看| 久久久久久久精品精品| 日本午夜av视频| 久久人妻熟女aⅴ| 少妇被粗大的猛进出69影院| 午夜激情av网站| 啦啦啦在线观看免费高清www| 免费少妇av软件| 久久久久精品性色| 在线观看美女被高潮喷水网站| av片东京热男人的天堂| 大陆偷拍与自拍| 午夜福利乱码中文字幕| 咕卡用的链子| 9色porny在线观看| 午夜av观看不卡| 欧美激情高清一区二区三区 | 丰满少妇做爰视频| 18禁裸乳无遮挡动漫免费视频| 亚洲人成电影观看| 久久久久久久久久人人人人人人| 国产一区有黄有色的免费视频| 亚洲欧美中文字幕日韩二区| 91aial.com中文字幕在线观看| 国产精品av久久久久免费| 自线自在国产av| 色播在线永久视频| 制服人妻中文乱码| av免费在线看不卡| 2022亚洲国产成人精品| 性高湖久久久久久久久免费观看| 精品视频人人做人人爽| 晚上一个人看的免费电影| 一区二区三区乱码不卡18| 久久热在线av| 性少妇av在线| 精品福利永久在线观看| 精品酒店卫生间| 欧美激情 高清一区二区三区| 97在线人人人人妻| 国产黄色视频一区二区在线观看| 日韩免费高清中文字幕av| 国产精品欧美亚洲77777| 少妇人妻 视频| 麻豆av在线久日| 香蕉精品网在线| av免费观看日本| 成人亚洲欧美一区二区av| 两个人免费观看高清视频| 少妇人妻精品综合一区二区| 90打野战视频偷拍视频| 久久久a久久爽久久v久久| 免费日韩欧美在线观看| 日韩三级伦理在线观看| 午夜福利网站1000一区二区三区| 99re6热这里在线精品视频| 飞空精品影院首页| 免费播放大片免费观看视频在线观看| 免费在线观看完整版高清| 女性被躁到高潮视频| 搡女人真爽免费视频火全软件| 日本-黄色视频高清免费观看| 免费在线观看完整版高清| 王馨瑶露胸无遮挡在线观看| www.精华液| 这个男人来自地球电影免费观看 | 激情视频va一区二区三区| 黑人欧美特级aaaaaa片| 日本午夜av视频| 性高湖久久久久久久久免费观看| 80岁老熟妇乱子伦牲交| 国产亚洲一区二区精品| 丝袜美足系列| 亚洲国产欧美在线一区| 菩萨蛮人人尽说江南好唐韦庄| 91aial.com中文字幕在线观看| 中文天堂在线官网| 永久免费av网站大全| 欧美黄色片欧美黄色片| 妹子高潮喷水视频| av福利片在线| videosex国产| 99国产精品免费福利视频| 久久久a久久爽久久v久久| 99热全是精品| 中文字幕精品免费在线观看视频| 国产一区二区 视频在线| 最新中文字幕久久久久| 男的添女的下面高潮视频| 亚洲成人av在线免费| 欧美日本中文国产一区发布| 成人二区视频| 黄网站色视频无遮挡免费观看| www日本在线高清视频| 最近手机中文字幕大全| 一二三四中文在线观看免费高清| 久久久久久久久免费视频了| 亚洲国产精品一区三区| 欧美另类一区| 多毛熟女@视频| 午夜日韩欧美国产| 免费女性裸体啪啪无遮挡网站| 女性生殖器流出的白浆| 久久精品国产亚洲av天美| 五月天丁香电影| 国产免费视频播放在线视频| 久久久久国产一级毛片高清牌| 日本-黄色视频高清免费观看| 中文字幕另类日韩欧美亚洲嫩草| 国产高清不卡午夜福利| 亚洲一区二区三区欧美精品| 亚洲国产av新网站| 免费久久久久久久精品成人欧美视频| 啦啦啦视频在线资源免费观看| 午夜老司机福利剧场| 亚洲欧美清纯卡通| 水蜜桃什么品种好| 韩国精品一区二区三区| 美女午夜性视频免费| 黄片小视频在线播放| 成年美女黄网站色视频大全免费| 日韩一本色道免费dvd| av福利片在线| 久久这里只有精品19| 爱豆传媒免费全集在线观看| 久久久久久久精品精品| 777久久人妻少妇嫩草av网站| 啦啦啦啦在线视频资源| tube8黄色片| 99国产精品免费福利视频| 亚洲伊人色综图| 18禁观看日本| 午夜精品国产一区二区电影| 国产日韩欧美视频二区| 国产精品免费视频内射| 国产精品久久久久久av不卡| 观看美女的网站| 国产人伦9x9x在线观看 | 亚洲精品国产色婷婷电影| 日韩制服骚丝袜av| av网站免费在线观看视频| h视频一区二区三区| 亚洲三区欧美一区| 国产老妇伦熟女老妇高清| 一二三四中文在线观看免费高清| 精品久久蜜臀av无| 欧美日韩成人在线一区二区| 满18在线观看网站| 免费高清在线观看日韩| 国产探花极品一区二区| 成人午夜精彩视频在线观看| 午夜免费男女啪啪视频观看| 99国产综合亚洲精品| 纵有疾风起免费观看全集完整版| 大香蕉久久成人网| 看非洲黑人一级黄片| 交换朋友夫妻互换小说| 天天躁夜夜躁狠狠久久av| 看非洲黑人一级黄片| xxx大片免费视频| 久久久久人妻精品一区果冻| av在线播放精品| 夜夜骑夜夜射夜夜干| 热99国产精品久久久久久7| 午夜精品国产一区二区电影| kizo精华| 五月开心婷婷网| 一区二区三区精品91| 一级毛片 在线播放| 成人影院久久| 国产精品嫩草影院av在线观看| 韩国精品一区二区三区| 亚洲成av片中文字幕在线观看 | 69精品国产乱码久久久| 多毛熟女@视频| 亚洲欧美清纯卡通| 欧美成人午夜免费资源| 在线精品无人区一区二区三| 亚洲人成77777在线视频| 欧美精品人与动牲交sv欧美| 亚洲熟女精品中文字幕| 亚洲欧洲精品一区二区精品久久久 | 免费高清在线观看日韩| 街头女战士在线观看网站| 最近最新中文字幕免费大全7| 天美传媒精品一区二区| 97精品久久久久久久久久精品| 亚洲精品一区蜜桃| 亚洲,欧美精品.| 老鸭窝网址在线观看| 亚洲精品乱久久久久久| 99久久综合免费| 一区二区三区激情视频| 国产不卡av网站在线观看| 精品国产乱码久久久久久男人| 伊人亚洲综合成人网| 久久久久精品久久久久真实原创| www.精华液| 国产成人免费无遮挡视频| 妹子高潮喷水视频| 久热这里只有精品99| 狠狠婷婷综合久久久久久88av| 亚洲美女视频黄频| 丰满迷人的少妇在线观看| 国产精品久久久久成人av| 三上悠亚av全集在线观看| 波多野结衣av一区二区av| 久久精品人人爽人人爽视色| 久久精品亚洲av国产电影网| 亚洲欧洲日产国产| 国产男女内射视频| 亚洲av综合色区一区| 天堂俺去俺来也www色官网| 亚洲第一av免费看| 九草在线视频观看| 国产精品一国产av| 一边亲一边摸免费视频| 黑人猛操日本美女一级片| 亚洲精品av麻豆狂野| 国产乱来视频区| 免费少妇av软件| 免费观看av网站的网址| 一本色道久久久久久精品综合| 美国免费a级毛片| 一级黄片播放器| 欧美97在线视频| 90打野战视频偷拍视频| 欧美亚洲 丝袜 人妻 在线| 有码 亚洲区| 热99国产精品久久久久久7| 亚洲少妇的诱惑av| 久久久久久久久久久免费av| 一区二区三区精品91| 国产成人免费观看mmmm| 国产免费一区二区三区四区乱码| 美女视频免费永久观看网站| 久久久国产欧美日韩av| 国产精品.久久久| 日产精品乱码卡一卡2卡三| 最近最新中文字幕免费大全7| 美国免费a级毛片| 久久精品夜色国产| 国产成人精品婷婷| 精品少妇久久久久久888优播| 成人毛片a级毛片在线播放| 伊人亚洲综合成人网| 91aial.com中文字幕在线观看| 精品久久久久久电影网| 欧美bdsm另类| 亚洲国产精品一区三区| 美女福利国产在线| 亚洲欧美日韩另类电影网站| 岛国毛片在线播放| 国产av精品麻豆| 你懂的网址亚洲精品在线观看| 宅男免费午夜| 午夜福利在线观看免费完整高清在| 精品人妻偷拍中文字幕| 哪个播放器可以免费观看大片| 久久这里只有精品19| 黑人猛操日本美女一级片| 国产精品熟女久久久久浪| 亚洲国产毛片av蜜桃av| 精品少妇黑人巨大在线播放| 国产av精品麻豆| 欧美日韩视频精品一区| 亚洲av福利一区| av电影中文网址| 国产 精品1| 国产野战对白在线观看| 两个人免费观看高清视频| 国产男女超爽视频在线观看| 久久人妻熟女aⅴ| av免费在线看不卡| 精品少妇内射三级| 久久人人97超碰香蕉20202| 中文字幕人妻丝袜一区二区 | 日产精品乱码卡一卡2卡三| 精品国产乱码久久久久久小说| 97在线视频观看| 国产精品久久久久久精品电影小说| 自线自在国产av| 久久久欧美国产精品| 最近中文字幕2019免费版| 久久99精品国语久久久| 一本色道久久久久久精品综合| xxxhd国产人妻xxx| 久久久久网色| 999精品在线视频| 国产免费又黄又爽又色| 人人妻人人澡人人爽人人夜夜| 电影成人av| 搡女人真爽免费视频火全软件| 丰满乱子伦码专区| 热99国产精品久久久久久7| 国产在线一区二区三区精| 1024香蕉在线观看| 日韩av不卡免费在线播放| 成人国语在线视频| 中文字幕精品免费在线观看视频| 老司机亚洲免费影院| 18禁观看日本| 国产成人精品在线电影| 欧美精品人与动牲交sv欧美| 亚洲av福利一区| 美女福利国产在线| 三上悠亚av全集在线观看| 免费黄频网站在线观看国产| 蜜桃在线观看..| 久久99蜜桃精品久久| 亚洲国产欧美网| 久久久欧美国产精品| 亚洲成色77777| 99精国产麻豆久久婷婷| 日韩三级伦理在线观看| 国产 精品1| 老汉色av国产亚洲站长工具| 菩萨蛮人人尽说江南好唐韦庄| 亚洲四区av| 成年动漫av网址| 免费观看a级毛片全部| 久久国内精品自在自线图片| 亚洲精华国产精华液的使用体验| 26uuu在线亚洲综合色| 久久久久网色| 人妻一区二区av| 日韩中字成人| 美女国产高潮福利片在线看| 黄色一级大片看看| 国产亚洲午夜精品一区二区久久| 久久久久久人妻| 男的添女的下面高潮视频| 亚洲精品av麻豆狂野| videos熟女内射| 日韩一本色道免费dvd| 99九九在线精品视频| h视频一区二区三区| 黄片无遮挡物在线观看| 我要看黄色一级片免费的| 热99久久久久精品小说推荐| 精品久久久久久电影网| 日韩一区二区三区影片| 美女脱内裤让男人舔精品视频| 国产精品女同一区二区软件| 亚洲欧美色中文字幕在线| 国产精品一国产av| 成年人午夜在线观看视频| 女人久久www免费人成看片| 亚洲精品国产色婷婷电影| videos熟女内射| 国产精品一国产av| 女性被躁到高潮视频| 女人久久www免费人成看片| 91久久精品国产一区二区三区| 欧美成人午夜免费资源| 欧美精品人与动牲交sv欧美| 嫩草影院入口| 日韩在线高清观看一区二区三区| 18禁裸乳无遮挡动漫免费视频| 亚洲国产欧美日韩在线播放| 大话2 男鬼变身卡| 国产成人免费无遮挡视频| 久久毛片免费看一区二区三区| 久久精品国产自在天天线| 亚洲,一卡二卡三卡| 欧美成人精品欧美一级黄| 老鸭窝网址在线观看| 久久精品国产综合久久久| 99re6热这里在线精品视频| 777米奇影视久久| 国产欧美日韩一区二区三区在线| 亚洲精品视频女| 日本欧美国产在线视频| videosex国产| 欧美日韩视频高清一区二区三区二| 日韩电影二区| 99热网站在线观看| 亚洲国产欧美日韩在线播放| 老鸭窝网址在线观看| 精品久久久精品久久久| 久久精品国产亚洲av涩爱| www.精华液| 久久国产精品男人的天堂亚洲| 免费观看a级毛片全部| 最近2019中文字幕mv第一页| 丝袜美腿诱惑在线| 啦啦啦中文免费视频观看日本| 日韩欧美一区视频在线观看| 99久久中文字幕三级久久日本| 久久久久久久大尺度免费视频| 亚洲精品久久久久久婷婷小说| 日本欧美国产在线视频| 久久人人爽人人片av| 中文字幕另类日韩欧美亚洲嫩草| 国产av国产精品国产| 1024香蕉在线观看| 久久久国产一区二区| 精品久久久精品久久久| 国产毛片在线视频| 2018国产大陆天天弄谢| 日日啪夜夜爽| 欧美日韩精品成人综合77777| 国产日韩欧美视频二区| 精品福利永久在线观看| 99久久精品国产国产毛片| 成人亚洲精品一区在线观看| 人人妻人人添人人爽欧美一区卜| 成人午夜精彩视频在线观看| 在线观看三级黄色| 在线观看人妻少妇| 亚洲欧洲国产日韩| 飞空精品影院首页| 青草久久国产| 97精品久久久久久久久久精品| 制服丝袜香蕉在线| 亚洲欧美一区二区三区久久| 国产成人精品福利久久| 在线观看美女被高潮喷水网站| av视频免费观看在线观看| xxx大片免费视频| 日本91视频免费播放| 日本wwww免费看| 久久久久久人人人人人| 日本av免费视频播放| 99热全是精品| 亚洲激情五月婷婷啪啪| 美女国产高潮福利片在线看| 国产熟女欧美一区二区| 毛片一级片免费看久久久久| 久久久久精品人妻al黑| 精品第一国产精品| 午夜福利一区二区在线看| 欧美日韩亚洲国产一区二区在线观看 | av天堂久久9| 人妻 亚洲 视频| 久久久久国产一级毛片高清牌| 涩涩av久久男人的天堂| 少妇人妻久久综合中文| 18禁观看日本| av免费观看日本| 人妻系列 视频| freevideosex欧美| 亚洲欧美一区二区三区国产| 午夜福利网站1000一区二区三区| 久久国内精品自在自线图片| 99国产综合亚洲精品| 精品卡一卡二卡四卡免费| 日韩熟女老妇一区二区性免费视频| 久久久久久久精品精品| 人妻一区二区av| www.熟女人妻精品国产| 新久久久久国产一级毛片| 亚洲国产av影院在线观看| 九九爱精品视频在线观看| 日韩精品有码人妻一区| 18禁观看日本| 日韩伦理黄色片| 久久久久久人妻| 免费在线观看黄色视频的| 18禁观看日本| 一区二区av电影网| 熟妇人妻不卡中文字幕| 久久精品久久久久久噜噜老黄| 久久久久久人人人人人| 欧美人与性动交α欧美精品济南到 | 黄网站色视频无遮挡免费观看| 日韩人妻精品一区2区三区| 最近2019中文字幕mv第一页| 日韩av在线免费看完整版不卡| 黄色视频在线播放观看不卡| 母亲3免费完整高清在线观看 | 国产又色又爽无遮挡免| 欧美另类一区| 国产又色又爽无遮挡免| 欧美少妇被猛烈插入视频| 另类亚洲欧美激情| 91久久精品国产一区二区三区| 热99久久久久精品小说推荐| 国产成人精品婷婷| 亚洲三级黄色毛片| 国产激情久久老熟女| 国产av一区二区精品久久| 国产视频首页在线观看| 黄片小视频在线播放| 9色porny在线观看| 99国产精品免费福利视频| 日本av手机在线免费观看| 在线天堂中文资源库| 久久 成人 亚洲| 国产免费现黄频在线看| 精品少妇内射三级| 欧美最新免费一区二区三区| 国语对白做爰xxxⅹ性视频网站| 大陆偷拍与自拍| 一区二区三区乱码不卡18| 免费女性裸体啪啪无遮挡网站| 免费日韩欧美在线观看| 国产女主播在线喷水免费视频网站| 丝袜脚勾引网站| 男女边摸边吃奶| 国产乱来视频区| 亚洲熟女精品中文字幕| 精品99又大又爽又粗少妇毛片| 日本wwww免费看| 久久久国产欧美日韩av| av又黄又爽大尺度在线免费看| 久久精品熟女亚洲av麻豆精品| av一本久久久久| 国产精品二区激情视频| 亚洲,一卡二卡三卡| 看十八女毛片水多多多| 日韩欧美精品免费久久| 亚洲欧美一区二区三区黑人 | 国产成人精品久久二区二区91 | 久久99热这里只频精品6学生| 九草在线视频观看| 久久 成人 亚洲| 99国产综合亚洲精品| www日本在线高清视频| 日韩人妻精品一区2区三区| 精品少妇一区二区三区视频日本电影 | 中文字幕亚洲精品专区| 国产一区二区在线观看av| 在线观看一区二区三区激情| 国产男女内射视频| 18+在线观看网站| 天美传媒精品一区二区| 亚洲,欧美,日韩| 一二三四中文在线观看免费高清| 夫妻性生交免费视频一级片| 高清不卡的av网站| 亚洲人成77777在线视频| 2021少妇久久久久久久久久久| 母亲3免费完整高清在线观看 | 99国产精品免费福利视频| 在线观看免费高清a一片| 中文字幕人妻丝袜制服| 最近2019中文字幕mv第一页| 久久99蜜桃精品久久| 成年av动漫网址| 18禁裸乳无遮挡动漫免费视频| 精品一区二区三卡| 国产一区二区 视频在线| 国产一区二区激情短视频 | 女人被躁到高潮嗷嗷叫费观| 极品人妻少妇av视频| 精品一品国产午夜福利视频| 丰满迷人的少妇在线观看| 男人添女人高潮全过程视频| 久久ye,这里只有精品| 看免费成人av毛片| 性色avwww在线观看| 99国产综合亚洲精品| a 毛片基地| 午夜av观看不卡| 国产乱人偷精品视频| 国产男女内射视频| 伦精品一区二区三区| 亚洲欧美一区二区三区久久| 亚洲精品aⅴ在线观看|