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

    高層建筑荷載引發(fā)地面沉降與隆起變形三維數(shù)值模擬

    2014-12-25 09:57:14付延玲金瑋澤陳興賢談金忠
    吉林大學學報(地球科學版) 2014年5期
    關鍵詞:泊松比滲透系數(shù)增量

    付延玲,金瑋澤,陳興賢,2,談金忠

    1.河海大學地球科學與工程學院,南京 210098

    2.江蘇省地質礦產勘查局,南京 210018

    3.江蘇南京地質工程勘察院,南京 210041

    0 引言

    地面沉降是一種由多種因素引起的地面高程緩慢降低的地質現(xiàn)象,其中過量開采地下水是引發(fā)地面沉降問題的主要原因。但近年來,隨著城市建設進程的加快,高層、超高層建筑的不斷增多,建筑荷載引發(fā)的地面沉降問題也日益突出。研究者們[1-4]認為,在上海,高群體建筑引起的地面沉降占總沉降量的30%左右,且沉降速度與建筑面積增長速度之間存在線性相關關系。

    高層建筑荷載一方面直接導致土體有效應力增加,發(fā)生固結壓縮,引發(fā)地面沉降;另一方面在地下水埋深較淺的地區(qū),高層建筑荷載的作用也會對地下水滲流場產生影響,導致滲流場重新調整,滲流場的變化同樣會引起應力場的變化,引發(fā)地面沉降或隆起變形。以往學者的研究更多集中在由抽水作用引起地下水滲流場變化而引發(fā)的地面沉降變形方面[5-7],對由高層建筑荷載引發(fā)的土體變形特征,尤其是高層建筑荷載短時間作用下引起滲流場調整變化,進而形成的土體變形的研究尚不多見。

    為了準確模擬預測由高層建筑荷載引發(fā)的地面沉降與隆起變形問題,筆者以比奧固結理論為基礎,綜合考慮了含水層和隔水層的流變特性,以及水力學和土力學參數(shù)隨滲流場和應力場發(fā)生變化的動態(tài)變化問題,建立了高層建筑荷載引發(fā)地面沉降與隆起變形的三維有限元數(shù)值模型,對高層建筑荷載引起的地面沉降與隆起變形進行了模擬預測分析,探討了變形過程中土體水力參數(shù)與力學參數(shù)的變化特征。

    1 理論基礎與方法

    1.1 比奧固結理論

    飽和土體(假定其土骨架變形為線彈性、微小變形、滲流符合達西定律、水不可壓縮或微壓縮)的三維比奧固結方程[8]如下:

    式中:G為剪切模量;ν為泊松比;wx、wy、wz分別為x、y、z方向上的位移分量;u為孔隙水壓力;Kx、Ky、Kz分別為x、y、z方向上的滲透系數(shù);γ為土的重度;γw為水的重度;▽為拉普拉斯算子,▽2=為時間。

    1.2 土體本構模型

    土的本構關系是土的力學特性即應力-應變-強度-時間等關系的數(shù)學表達式。對于考慮流變特性的土體來說,其變形特征主要表現(xiàn)為變形的時間與應力水平有關,所顯示的是具有彈性、塑性和黏滯性的黏彈塑性體,若將此類土體的總應變增量dε分為彈塑性應變增量dεep、黏彈性應變增量dεve和黏塑性應變增量dεvp,則具有流變特性的土體中任意點在任意時刻的應變增量[7,9]為

    式(3)中各部分的應變增量可以由下述方法確定。

    1)彈塑性應變增量

    由土體的彈塑性本構模型可得

    2)黏彈性應變增量

    由kelvin流變模型E1ε+Keε=σ,可得應力不變時復雜應力狀態(tài)下的黏彈性應變增量:

    式中:ηe=E1/Ke,E1為kelvin體黏彈性模量,Ke為黏滯系數(shù);A為應力矩陣,

    E為彈性模量。

    3)黏塑性應變增量

    利用黏塑性法確定黏塑性應變增量,在該種方法中允許材料在有限“期間”內超越破壞準則(以破壞準則函數(shù)F的值大于0來表示)。在討論土體的黏塑性應變而非塑性應變時,應變的變化率與超越量有關,即有以下關系式:

    式中:Qs為塑性勢函數(shù)對于摩爾-庫倫材料來說,

    其中:φ為內摩擦角;c為黏聚力;

    β=為第三偏應力不變量,

    σ與τ分別為主應力與剪應力。

    如果將黏塑性應變率與一個偽時間步相乘,就可以得到累加到下一個時間步的黏塑性應變增量,于是有

    數(shù)值計算與絕對穩(wěn)定的時間步與假定的破壞準則有關。

    對于摩爾-庫倫材料有

    塑性勢函數(shù)對應力的偏導數(shù)可以表示為

    所以,應用摩爾-庫倫準則的土體黏塑性應變增量可以表示為

    將式(4)、(5)、(10)代入式(3),可求得彈塑-黏彈-黏塑性體的應變增量為

    式(11)即為土體的黏彈塑性本構方程。

    1.3 比奧固結有限元方程

    利用伽遼金加權余量法離散方程,考慮到土體的非線性特性,取Δt時間內的位移增量來代替位移,將式(1)、(2)離散成增量形式[10]:

    式中:Δδ為結點位移增量;Δu為結點孔隙壓力增量;ˉK為固體剛度矩陣;K為滲透流量矩陣;K′為應力-滲流耦合項矩陣;ΔQ為流量增量矩陣;B為自由面的積分矩陣;R為等效節(jié)點荷載,當存在高層建筑荷載時,R為包括高層建筑荷載引起的附加應力值;Rt為t時刻已經發(fā)生的位移所平衡了的那部分荷載。

    因為滲流取決于孔隙壓力全量的分布,而不是取決于時間內孔隙壓力增量,所以孔壓要用全量的形式表示。記時刻tn和tn+1時單元節(jié)點i的孔壓全量分別為ui(n)和ui(n+1),且 Δui=ui(n+1)-ui(n),則式(12)可變換為

    式(13)即為三維比奧固結有限元方程。

    1.4 參數(shù)的動態(tài)變化模型

    1.4.1 孔隙度與滲透系數(shù)的非線性

    流固耦合問題實際上是孔隙應力的消散引起土體骨架的變形,導致滲透系數(shù)變化,從而影響土體的滲透性,宏觀上表現(xiàn)為土體的固結變形。在比奧固結的假定條件下,根據孔隙度的相關定義和滲流力學Kozeny-Carman方程推得孔隙度n和滲透系數(shù)K的動態(tài)表達式[11-12]:

    式中:n0為初始孔隙度;K0為初始滲透系數(shù);εv為體應變

    1.4.2 土體參數(shù)的非線性

    采用鄧肯-張非線性模型,將土體的本構關系推廣到非線性,則本構關系{Δσ}=D{Δε}中矩陣D中的彈性常數(shù)E、ν不再視為常量,而是隨著應力狀態(tài)改變而改變,其切線彈性模量Et和切線泊松比νt的表達式如下[13]:

    式中:Rf為破壞比;σ1為第一主應力;σ3為第三主應力;ω為彈性模量與固結壓力曲線的斜率;logα,G為土體常規(guī)三軸壓縮實驗結果所繪曲線截距;I=0.04,D=3為土體實驗參數(shù);p為大氣壓強。

    1.5 定解條件

    1.5.1 初始條件

    1)地應力初始條件

    采用土體的自重應力估算土體的初始應力:

    式中:σx、σz為土體的初始水平向和垂向應力;z為計算點深度;K0為靜止側壓力系數(shù),K0=。為有效內摩擦角

    2)位移初始條件

    3)孔隙水壓力初始條件

    式中:u0(x,y,z)為研究區(qū)域內已知初始孔隙水壓力。

    1.5.2 邊界條件

    1)孔隙水壓力邊界條件Γ1

    式中:us為水頭邊界Γ1上的已知孔隙水壓力。

    2)流量邊界條件Γ2

    式中:qL為邊界Γ2上的已知單位面積流量。

    3)自由面邊界條件Γ3

    式中:μ為土體給水度;θ為自由面外法線方向與垂線的交角;q為通過自由面邊界Γ3的單位面積流量;Z為自由面所在的高程。

    4)位移邊界條件Γ4

    比奧固結有限元方程結合定解條件和水力學參數(shù)及土力學參數(shù)動態(tài)變化模型即可運用Fortran語言編制相應的有限元程序進行求解[5,14]。

    2 應用算例

    2.1 概念模型

    為了研究高層建筑荷載影響下的地面沉降情況,排除其他影響因素干擾,建立一個簡單的均質含水層模型。模型長×寬×高為500m×500m×100 m,建筑物按每層自重1.2t/m2、承重每層0.3t/m2計算,地基尺寸的長×寬為70m×50m。用八節(jié)點六面體單元離散化模型,在平面上剖分為2 500個矩形網格單元,垂向上從上往下概化為:潛水含水層、第1承壓含水層,第2承壓含水層、第3承壓含水層及各含水層之間的黏性土弱含水層,共7層。每層土體劃分為一個參數(shù)分區(qū),共7個參數(shù)分區(qū)。其中:潛水含水層巖性以粉砂、亞砂土為主,底板埋深20~30m,靜水位埋深0.9~1.2m;第1承壓含水層以粉砂為主,底板埋深50~52m,靜水位埋深3.2~3.8m;第2承壓含水層以亞砂土、細中砂為主,底板埋深79~83m,靜水位埋深7.6~8.1m;第3承壓含水層以細中砂、中粗砂為主,底板埋深96~100m,靜水位埋深11.5~13.4m。模型四周均概化為第一類已知水頭邊界,底部概化為隔水邊界。建筑物位置及模型分層如圖1所示。模型各參數(shù)分區(qū)物理力學性質參數(shù)如表1所示。其中:K0x、K0y為初始水平向滲透系數(shù);K0z為初始垂向滲透系數(shù);ν0為初始泊松比;E0為初始彈性模量;Sy為儲水率。

    圖1 模型示意圖Fig.1 Model sketch

    2.2 建筑荷載影響下的地面沉降與隆起變形數(shù)值模擬

    模型建筑物高度按10層計算。模型計算周期為1a,劃分為15個應力期。其中:第1到第4應力期分別為5,10,20,31d;第5到15應力期每個應力期時間分別為第2月到第12月單月時間。假定初始情況下的各參數(shù)分區(qū)孔隙度大小均相同,且初始孔隙度n0=0.36。選取建筑物中心點所在剖面(y=260m)模型運行15個應力期時的地面沉降量,如圖2所示,其中地面變形量正值表示隆起量,負值表示沉降量。

    圖2 不同應力期末的地面沉降曲線圖Fig.2 Land subsidence curve of difference stress period

    土體的固結壓縮變形較為復雜,受到土體本身性質、邊界條件、排水條件以及上部荷載方式等因素影響。由圖2可以看出,由高層建筑荷載引起的地面沉降以建筑物中心點為中心呈現(xiàn)漏斗狀:到第4個應力期末,即31d時,最大沉降達到-3.8mm;在第10個應力期末,即第7個月末時,最大沉降量為-8.5mm。由圖2b可以看出:模型運行前9個應力期時,高層建筑荷載引起的地面沉降在x=200~300m和x=400~500m區(qū)間段均存在隆起現(xiàn)象;在第10個應力期末,即第7個月末時,隆起現(xiàn)象消失。而由圖2a可知:在高層建筑荷載作用5,10,20d時隆起量依次增大,最大值分別為0.9mm,1.0mm,1.5mm;31d時隆起量為1.3mm??梢钥闯觯?0d累計隆起量相對于31d累計隆起量較大。這是因為由于高層建筑荷載主要影響淺部地層,在模型中,淺部土體垂向滲透系數(shù)較小,所以淺部地層的孔隙水在受到高層建筑荷載作用時,更多地以橫向流動為主;而高層建筑荷載作用時間較短,其周圍土體中的孔隙水急劇增多,導致超靜孔隙水壓力在初期急劇增大,有效應力相應減小,土體骨架膨脹變形,必然使得建筑周圍土體發(fā)生回彈,地表出現(xiàn)隆起;且隨著應力期的增加,建筑周圍土體的超靜孔隙水壓力緩慢消散,地面隆起緩慢消失。

    2.3 參數(shù)變化特征分析

    隨著土體的固結壓縮與回彈變形,建筑物周圍及下部土體力學及水力學參數(shù)均會伴隨著土體的變形而發(fā)生相應變化。選取模型第1161號單元和第1167號單元進行孔隙度、滲透系數(shù)、彈性模量及泊松比分析。其中:第1161號單元位于建筑物中心底部;第1167號單元中心點坐標為(430,260),在剖面y=200m上,距離建筑物40m,且兩單元均位于模型第一層,單元中心線位于y=260m剖面上。兩單元孔隙度及滲透系數(shù)隨應力期的變化如圖3所示。

    由于1161號單元與1167號單元所處區(qū)域地面沉降變化情況不同,1161號單元位于模型第1層地面沉降中心點位置,而1167號單元所處區(qū)域在前9個應力期均出現(xiàn)隆起現(xiàn)象,所以兩單元的孔隙度變化趨勢是不同的。由圖3a可以看出:1167號單元孔隙度在20d增大到最大值,隨著應力期的增加,孔隙度緩慢減??;而1161號單元在整個應力期時間段內,孔隙度一直處于減小的趨勢。造成這種現(xiàn)象的原因可能是:在建筑荷載作用初期,建筑物下部土體的孔隙水短時間擴散到周圍淺部地層,使得孔隙水大量聚集,超靜孔隙水壓力急劇增大,土體有效應力相應減小,淺部土體發(fā)生膨脹,使得建筑物周圍土體孔隙度在20d達到最大值;但隨后超靜孔隙水壓力逐漸消散,使得孔隙度緩慢減小,在第10個應力期,即第7個月時恢復到初始孔隙度大小。建筑物中心點下部土體在受到建筑荷載作用受到擠壓后,土體孔隙水向四周擴散,使得孔隙水壓力一直減小,有效應力一直增加,所以孔隙度隨著應力期的變化始終減小。

    滲透系數(shù)隨著孔隙度的變化而發(fā)生變化??紫抖仍龃?,土體出水能力增強,導致滲透系數(shù)變大,而孔隙度減小,使土體出水能力變弱,從而導致滲透系數(shù)減小。由圖3也可以看出,1161號單元和1167號單元的滲透系數(shù)與孔隙度變化趨勢相同。

    表1 地層參數(shù)Table 1 Stratum parameters

    圖3 1161號單元和1167號單元孔隙度、水平向滲透系數(shù)和垂向滲透系數(shù)隨應力期變化曲線圖Fig.3 Variation curve among porosity and horizontal hydraulic conductivity and vertical hydraulic conductivity and stress periods of 1161element and 1167element

    伴隨著土體的回彈及壓縮,彈性模量、泊松比必然發(fā)生變化。由圖4可以看出,1167號單元彈性模量、泊松比的變化與孔隙度的變化相對應,在20d左右土體彈性模量達到最小值,泊松比達到最大值,隨著應力期增加彈性模量緩慢增大,泊松比緩慢減小。而1161號單元彈性模量及泊松比的變化均較為平緩,彈性模量不斷增大,泊松比不斷減小。通過對土體力學及水力學參數(shù)的變化趨勢分析可以發(fā)現(xiàn),在前9個應力期內各參數(shù)的變化相比后幾個應力期要快,這是因為隨著應力期的增長,土體固結變形量逐漸減小,導致土體力學及水力參數(shù)的變化相對于前期而言越來越不明顯。

    3 結論

    1)以比奧固結理論為基礎,結合土體的非線性流變理論,將土體的本構關系推廣到黏彈塑性,同時考慮土體水力學參數(shù)和土力學參數(shù)隨有效應力的動態(tài)變化問題,建立起的地下水滲流與土體變形三維全耦合模型更加符合實際,進一步提高了模型計算的精度。

    2)由高層建筑荷載引起的地面沉降呈現(xiàn)漏斗狀,以建筑物中心點為漏斗中心。建筑荷載作用初期會引起周圍土體隆起,建筑物周圍淺部土層土體孔隙度、滲透系數(shù)及泊松比也隨著增加,到達峰值后,隨著時間的增加又緩慢減??;而彈性模量逐漸減小,到達最小值后,隨著時間的增加緩慢增大。建筑荷載下部土層,孔隙度、滲透系數(shù)及泊松比隨時間的增加均緩慢減小,而彈性模量隨時間的增加緩慢增大。

    (References):

    [1]龔士良.上海城市建設對地面沉降的影響[J].中國地質災害與防治學報,1998,9(2):108-111.Gong Shiliang.Urban Construction’s Impact on Land Subsidence in Shanghai[J].The Chinese Journal of Geological Hazard and Control,1998,9(2):108-111.

    [2]許燁霜,馬磊,沈水龍.上海市城市化進程引起的地面沉降因素分析[J].巖土力學,2011,32(增刊1):578-582.Xu Yeshuang,Ma Lei,Shen Shuilong.Influential Factors on Development of Land Subsidence with Process of Urbanization in Shanghai[J].Rock and Soil Mechanics,2011,32(Sup.1):578-582.

    [3]嚴學新,沈國平.上海城區(qū)建筑密度與地面沉降關系分析[J].水文地質工程地質,2002,29(6):21-25.Yan Xuexin,Shen Guoping.Relationship Between Building Density and Land Subsidence in Shanghai Urban Zone[J].Hydrogeology and Engineering Geology,2002,29(6):21-25.

    [4]唐益群,宋壽鵬,陳斌,等.不同建筑容積率下密集建筑群區(qū)地面沉降規(guī)律研究[J].巖石力學與工程學報,2010,29(增刊1):3425-3431.Tang Yiqun,Song Shoupeng,Chen Bin,et al.Study of Land Subsidence Rule of Dense Buildings Under Different Floor Area Ratios[J].Chinese Journal of Rock Mechanics and Engineering,2010,29(Sup.1):3425-3431.

    [5]陳興賢,駱祖江,安曉宇,等.深基坑降水三維變參數(shù)非穩(wěn)定滲流與地面沉降耦合模型[J].吉林大學學報:地球科學版,2013,43(5):1572-1578.Chen Xingxian,Luo Zujiang,An Xiaoyu,et al.Coupling Model of Groundwatrer Three Dimensional Variable Parametrics Non-Steady Seepage and Land-Subsidence[J].Journal of Jilin University:Earth Science Edition,2013,43(5);1572-1578.

    [6]Luo Zujiang,Zeng Feng.Finite Element Numerical Simulation of Land Subsidence and Groundwater Exploitation Based on Visco-Elastic-Plastic Biot’s Consolidation Theory[J].Journal of Hydrodynamics,2011,23(5):615-624.

    [7]錢家歡,殷宗澤.土工原理與計算[M].北京:水力水電出版社,1996.Qian Jiahuan,Yin Zongze.Principle and Canculation of Geotechnics[M].Beijing:China Waterpower Press,1996.

    [8]陳曉平,白世偉.軟黏土地基黏彈塑性比奧固結的數(shù)值分析[J].巖土工程學報,2001,23(4):481-484.Chen Xiaoping,Bai Shiwei.The Numerical Analysis of Visco-Elastic-Plastic Biot’s Consolidation for Soft Clay Foundation[J].Chinese Journal of Geotechnical Engineering,2001,23(4):481-484.

    [9]李醫(yī)民,周鳳燕.一類三維偏微分方程邊值問題的解法[J].江蘇大學學報:自然科學版,2004,25(4):328-331.Li Yimin,Zhou Fengyan.Solution to a Class of Boundary Value Problem of Three Dimensional Partial Differential Equation[J].Journal of Jiangsu University:Natural Science Edition,2004,25(4):328-331.

    [10]冉啟全,李士倫.流固耦合油藏數(shù)值模擬中物性參數(shù)動態(tài)模型研究[J].石油勘探與開發(fā),1997,24(3):61-65.Ran Qiquan,Li Shlun.Study on Dynamic Models of Reservoir Parameters in the Coupled Simulation of Multiphase Flow and Reservoir Deformation[J].Petroleum Exploration and Development,1997,24(3):61-65.

    [11]田杰,劉先貴,尚根華.基于流固耦合理論的套損力學機理分析[J].水動力學研究與進展:A 輯,2005,20(2):221-225.Tian Jie,Liu Xiangui,Shang Genhua.Casing Damage Mechanism Based on Theory of Fluid-Solid Coupling Flow Through Underground Rock[J].Journal of Hydrodynamics:Series A,2005,20(2):221-225.

    [12]羅剛,張建民.鄧肯-張模型和沈珠江雙屈服面模型的改進[J].巖土力學,2004,25(6):887-890.Luo Gang,Zhang Jianmin.Improvement of Duncan-Chang Nonlinear Model and Shen Zhujiang’s Elastoplastic Model for Granular Soils[J].Rock and Soil Mechanics,2004,25(6):887-890.

    [13]彭國倫.Fortran 95程序設計[M].北京:中國電力出版社,2005.Peng Guolun.Fortran 95[M].Beijing:China Electric Power Press,2005.

    [14]Smith I M,Griffiths.D V.有限元方法編程[M].3版.王菘,周堅鑫,王來,等,譯.北京:電子工業(yè)出版社,2003.Smith I M,Griffiths D V.Programming the Finite Element Method[M].3rd ed.Translated by Wang Song,Zhou Jianxin, Wang Lai,et al.Beijing:Electronic Industry Press,2003.

    猜你喜歡
    泊松比滲透系數(shù)增量
    提質和增量之間的“辯證”
    當代陜西(2022年6期)2022-04-19 12:12:22
    具有負泊松比效應的紗線研發(fā)
    基于Origin的滲透系數(shù)衰減方程在地熱水回灌中的應用
    地質與資源(2021年1期)2021-05-22 01:24:26
    負泊松比功能的結構復合紡紗技術進展
    毛紡科技(2020年6期)2021-01-06 03:41:48
    “價增量減”型應用題點撥
    考慮粘彈性泊松比的固體推進劑蠕變型本構模型①
    固體推進劑粘彈性泊松比應變率-溫度等效關系
    多孔材料水滲透系數(shù)預測的隨機行走法
    輸水渠防滲墻及基巖滲透系數(shù)敏感性分析
    基于均衡增量近鄰查詢的位置隱私保護方法
    電信科學(2016年9期)2016-06-15 20:27:25
    亚洲精品在线观看二区| 国产精品av久久久久免费| 男女做爰动态图高潮gif福利片| 亚洲精品久久成人aⅴ小说| 男女下面进入的视频免费午夜| 国内揄拍国产精品人妻在线| 美女黄网站色视频| 欧美黄色片欧美黄色片| 欧美日韩中文字幕国产精品一区二区三区| 精品国产乱码久久久久久男人| 欧美zozozo另类| 国模一区二区三区四区视频 | 久久亚洲真实| 亚洲精品国产一区二区精华液| 亚洲国产日韩欧美精品在线观看 | 在线a可以看的网站| 搡老岳熟女国产| 一个人免费在线观看的高清视频| 国产精品美女特级片免费视频播放器 | 人人妻,人人澡人人爽秒播| 波多野结衣高清作品| 欧美zozozo另类| 女人被狂操c到高潮| 少妇的丰满在线观看| 国产激情偷乱视频一区二区| 97碰自拍视频| 黄色视频,在线免费观看| 操出白浆在线播放| 91在线观看av| 欧美av亚洲av综合av国产av| 精品免费久久久久久久清纯| 美女扒开内裤让男人捅视频| 又紧又爽又黄一区二区| 亚洲男人天堂网一区| 久久精品成人免费网站| 九色国产91popny在线| 变态另类成人亚洲欧美熟女| 久久中文看片网| 一边摸一边抽搐一进一小说| 亚洲在线自拍视频| 天天躁夜夜躁狠狠躁躁| 视频区欧美日本亚洲| 亚洲人成电影免费在线| 亚洲中文日韩欧美视频| 亚洲成人精品中文字幕电影| 欧美在线一区亚洲| 欧美国产日韩亚洲一区| 动漫黄色视频在线观看| 国产一区二区激情短视频| 欧美色视频一区免费| 观看免费一级毛片| 日本黄色视频三级网站网址| 国产精品九九99| 在线观看日韩欧美| 亚洲人成伊人成综合网2020| 看免费av毛片| 亚洲精品中文字幕在线视频| or卡值多少钱| 国产成人欧美在线观看| 国产精品1区2区在线观看.| 国产不卡一卡二| 国产成人精品久久二区二区91| 欧美中文日本在线观看视频| 欧美激情久久久久久爽电影| 亚洲激情在线av| 亚洲精品久久国产高清桃花| 人妻丰满熟妇av一区二区三区| 亚洲精品在线美女| 女生性感内裤真人,穿戴方法视频| 两性夫妻黄色片| 99久久精品国产亚洲精品| 国产成人影院久久av| 国产av一区二区精品久久| 黄色女人牲交| 床上黄色一级片| 极品教师在线免费播放| 亚洲国产欧美一区二区综合| 国产熟女午夜一区二区三区| 久久久久国产一级毛片高清牌| 亚洲欧美日韩高清在线视频| 欧美高清成人免费视频www| 黄色成人免费大全| 久久久久亚洲av毛片大全| 日韩欧美一区二区三区在线观看| 俄罗斯特黄特色一大片| 一本一本综合久久| 国产精品亚洲av一区麻豆| 国产成人av激情在线播放| 欧美激情久久久久久爽电影| 亚洲av成人不卡在线观看播放网| 免费在线观看完整版高清| 亚洲七黄色美女视频| av超薄肉色丝袜交足视频| 日本免费一区二区三区高清不卡| 亚洲真实伦在线观看| 国产成人系列免费观看| 一二三四社区在线视频社区8| 岛国在线免费视频观看| 可以在线观看的亚洲视频| 欧美中文综合在线视频| 国产激情偷乱视频一区二区| 国产视频一区二区在线看| 国产野战对白在线观看| 亚洲国产欧洲综合997久久,| 亚洲七黄色美女视频| 三级毛片av免费| 国产又黄又爽又无遮挡在线| 无限看片的www在线观看| 精品无人区乱码1区二区| 人妻丰满熟妇av一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | 日本三级黄在线观看| 日日摸夜夜添夜夜添小说| 国产乱人伦免费视频| 婷婷精品国产亚洲av在线| 99精品欧美一区二区三区四区| 国产精品av久久久久免费| 黄色片一级片一级黄色片| 一本精品99久久精品77| 亚洲全国av大片| 最近视频中文字幕2019在线8| 成人国语在线视频| 国产精品亚洲一级av第二区| 午夜福利免费观看在线| 欧洲精品卡2卡3卡4卡5卡区| 久久精品91蜜桃| 国产av不卡久久| 精品无人区乱码1区二区| 老司机深夜福利视频在线观看| av在线播放免费不卡| a级毛片在线看网站| 一级作爱视频免费观看| 黄色丝袜av网址大全| 亚洲av五月六月丁香网| 久久精品夜夜夜夜夜久久蜜豆 | 欧美丝袜亚洲另类 | 一a级毛片在线观看| 俄罗斯特黄特色一大片| 亚洲人成电影免费在线| 怎么达到女性高潮| 亚洲精品在线美女| 中文字幕人成人乱码亚洲影| 搡老妇女老女人老熟妇| 国产视频内射| 黄色毛片三级朝国网站| 三级国产精品欧美在线观看 | 床上黄色一级片| 91国产中文字幕| 欧美大码av| 999久久久精品免费观看国产| 精品乱码久久久久久99久播| 18禁裸乳无遮挡免费网站照片| 国产人伦9x9x在线观看| 精品不卡国产一区二区三区| 国内揄拍国产精品人妻在线| cao死你这个sao货| 夜夜夜夜夜久久久久| 国产亚洲精品综合一区在线观看 | 毛片女人毛片| 国产久久久一区二区三区| 国产亚洲精品综合一区在线观看 | 精品无人区乱码1区二区| 成人精品一区二区免费| 成人av一区二区三区在线看| 欧美在线黄色| 琪琪午夜伦伦电影理论片6080| 两性夫妻黄色片| 伦理电影免费视频| 手机成人av网站| 日韩 欧美 亚洲 中文字幕| 成在线人永久免费视频| 无人区码免费观看不卡| 高潮久久久久久久久久久不卡| 国产精品98久久久久久宅男小说| 成人三级黄色视频| 舔av片在线| av福利片在线| 日韩大尺度精品在线看网址| 一区二区三区激情视频| 波多野结衣高清作品| 香蕉国产在线看| 啦啦啦观看免费观看视频高清| 老汉色∧v一级毛片| 久久香蕉国产精品| 精品欧美国产一区二区三| 两个人看的免费小视频| 欧美精品亚洲一区二区| 亚洲专区中文字幕在线| 99在线人妻在线中文字幕| 久久九九热精品免费| 一本精品99久久精品77| 午夜两性在线视频| 国产蜜桃级精品一区二区三区| 色播亚洲综合网| 久久久久国产一级毛片高清牌| 国产成+人综合+亚洲专区| 日韩精品青青久久久久久| 九色成人免费人妻av| 亚洲成av人片免费观看| 制服人妻中文乱码| 精品不卡国产一区二区三区| 黄色女人牲交| 日本 欧美在线| 色综合欧美亚洲国产小说| 搞女人的毛片| 日韩免费av在线播放| 在线免费观看的www视频| 麻豆av在线久日| 日本三级黄在线观看| 身体一侧抽搐| 国产精品久久视频播放| 免费在线观看视频国产中文字幕亚洲| 国产av不卡久久| 色综合欧美亚洲国产小说| 亚洲精品中文字幕在线视频| 久久天躁狠狠躁夜夜2o2o| 一区二区三区国产精品乱码| 欧美黑人欧美精品刺激| 欧美一级毛片孕妇| 亚洲av成人不卡在线观看播放网| 这个男人来自地球电影免费观看| 日韩国内少妇激情av| 国产精品一区二区三区四区久久| 亚洲国产中文字幕在线视频| 国产三级在线视频| 麻豆成人午夜福利视频| 国产成人精品久久二区二区免费| 欧美 亚洲 国产 日韩一| 国产高清videossex| 成人国产综合亚洲| 亚洲国产精品sss在线观看| 露出奶头的视频| 中文字幕av在线有码专区| 白带黄色成豆腐渣| 日本 av在线| 夜夜爽天天搞| 久久伊人香网站| 又紧又爽又黄一区二区| 亚洲,欧美精品.| 在线观看日韩欧美| 99国产精品99久久久久| 国产三级中文精品| 国产精品爽爽va在线观看网站| 国产区一区二久久| 国产精品乱码一区二三区的特点| 男女下面进入的视频免费午夜| 一个人免费在线观看电影 | 日日夜夜操网爽| 国内揄拍国产精品人妻在线| 国产单亲对白刺激| 欧美一区二区精品小视频在线| 啪啪无遮挡十八禁网站| 国产免费男女视频| 国产精品爽爽va在线观看网站| 亚洲天堂国产精品一区在线| 美女午夜性视频免费| 亚洲人与动物交配视频| 久久婷婷成人综合色麻豆| 最近最新中文字幕大全免费视频| 国产熟女xx| 亚洲男人天堂网一区| 欧美成人免费av一区二区三区| 人妻夜夜爽99麻豆av| cao死你这个sao货| 不卡av一区二区三区| 亚洲熟妇熟女久久| 精品第一国产精品| 色噜噜av男人的天堂激情| 99久久国产精品久久久| 老熟妇仑乱视频hdxx| 精品久久久久久久末码| 日韩精品中文字幕看吧| 叶爱在线成人免费视频播放| 国产激情偷乱视频一区二区| 国产av不卡久久| 色尼玛亚洲综合影院| 久久人人精品亚洲av| 亚洲色图 男人天堂 中文字幕| 亚洲国产精品合色在线| 在线观看美女被高潮喷水网站 | 18禁国产床啪视频网站| 特大巨黑吊av在线直播| 久久人人精品亚洲av| 久久久国产成人免费| av福利片在线观看| 亚洲精品中文字幕一二三四区| 成人精品一区二区免费| 国产三级黄色录像| 人妻丰满熟妇av一区二区三区| 久久精品国产亚洲av香蕉五月| 国产视频内射| 人人妻,人人澡人人爽秒播| 国产欧美日韩一区二区精品| 国产精品久久久久久久电影 | 天天添夜夜摸| 窝窝影院91人妻| 美女黄网站色视频| 制服人妻中文乱码| 午夜免费观看网址| 欧美国产日韩亚洲一区| 18禁黄网站禁片免费观看直播| 精品一区二区三区四区五区乱码| 大型黄色视频在线免费观看| 国产欧美日韩一区二区精品| 亚洲一区二区三区色噜噜| 久久天躁狠狠躁夜夜2o2o| АⅤ资源中文在线天堂| 久久久久精品国产欧美久久久| 搡老岳熟女国产| x7x7x7水蜜桃| 又粗又爽又猛毛片免费看| 精品欧美国产一区二区三| 女同久久另类99精品国产91| 欧美成人一区二区免费高清观看 | 一a级毛片在线观看| 夜夜夜夜夜久久久久| 亚洲 国产 在线| 高清毛片免费观看视频网站| 精品少妇一区二区三区视频日本电影| 久久精品91无色码中文字幕| 国产精品一区二区精品视频观看| 一个人观看的视频www高清免费观看 | 国产成人欧美在线观看| 日韩成人在线观看一区二区三区| 成年免费大片在线观看| 最好的美女福利视频网| 日本一区二区免费在线视频| 999精品在线视频| 十八禁人妻一区二区| 国产黄片美女视频| 一级片免费观看大全| 黑人操中国人逼视频| 精品国产美女av久久久久小说| 久9热在线精品视频| 特级一级黄色大片| 国产精品一区二区三区四区久久| 国产视频内射| 日韩中文字幕欧美一区二区| 少妇裸体淫交视频免费看高清 | 99久久精品热视频| 熟女少妇亚洲综合色aaa.| 无遮挡黄片免费观看| 无人区码免费观看不卡| 国产精品影院久久| 99久久精品国产亚洲精品| 在线观看www视频免费| 淫秽高清视频在线观看| 一本一本综合久久| 国产亚洲av嫩草精品影院| 国产探花在线观看一区二区| 亚洲美女黄片视频| www.www免费av| 国产精品电影一区二区三区| www日本黄色视频网| 麻豆国产av国片精品| 少妇的丰满在线观看| 精品国产亚洲在线| 欧美成人性av电影在线观看| 亚洲国产欧美网| 岛国视频午夜一区免费看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲 国产 在线| 免费看十八禁软件| 老司机深夜福利视频在线观看| 变态另类丝袜制服| 国产av一区在线观看免费| 日韩精品青青久久久久久| 亚洲精华国产精华精| 少妇熟女aⅴ在线视频| 久久久久久久午夜电影| 亚洲人成网站高清观看| 日韩成人在线观看一区二区三区| 亚洲人与动物交配视频| 国产精品久久久久久久电影 | 九色国产91popny在线| 午夜免费成人在线视频| 成年人黄色毛片网站| 欧美乱色亚洲激情| 亚洲精品色激情综合| 成人18禁在线播放| 久久人人精品亚洲av| 国产伦一二天堂av在线观看| 免费av毛片视频| 高清毛片免费观看视频网站| 一本综合久久免费| 美女高潮喷水抽搐中文字幕| 18禁黄网站禁片午夜丰满| 亚洲欧美精品综合一区二区三区| 亚洲18禁久久av| 欧美日韩亚洲国产一区二区在线观看| 午夜视频精品福利| 两个人视频免费观看高清| 亚洲 欧美一区二区三区| 在线观看一区二区三区| x7x7x7水蜜桃| 亚洲一区中文字幕在线| 久久婷婷人人爽人人干人人爱| 极品教师在线免费播放| 级片在线观看| 色噜噜av男人的天堂激情| 首页视频小说图片口味搜索| 老汉色∧v一级毛片| 又大又爽又粗| 岛国在线观看网站| 亚洲中文av在线| xxx96com| 91麻豆精品激情在线观看国产| 久久久久国产一级毛片高清牌| 午夜福利在线在线| 在线观看日韩欧美| 麻豆成人av在线观看| 狠狠狠狠99中文字幕| 色综合婷婷激情| 国产高清视频在线观看网站| 高清在线国产一区| 亚洲欧美日韩无卡精品| 久99久视频精品免费| 一级黄色大片毛片| 久久中文字幕一级| 日日夜夜操网爽| 色综合亚洲欧美另类图片| 成熟少妇高潮喷水视频| 精品人妻1区二区| 两性夫妻黄色片| 在线十欧美十亚洲十日本专区| videosex国产| 欧美一级a爱片免费观看看 | 久久草成人影院| 亚洲专区中文字幕在线| 亚洲色图 男人天堂 中文字幕| 欧美成人午夜精品| 日日干狠狠操夜夜爽| 91在线观看av| 久久久国产精品麻豆| 91老司机精品| 午夜免费观看网址| 亚洲avbb在线观看| 美女大奶头视频| 久9热在线精品视频| 成人国产综合亚洲| videosex国产| 日韩大尺度精品在线看网址| 国产精华一区二区三区| 久久 成人 亚洲| 国产私拍福利视频在线观看| 婷婷丁香在线五月| 最近最新中文字幕大全免费视频| 两人在一起打扑克的视频| 老司机靠b影院| 色精品久久人妻99蜜桃| 不卡av一区二区三区| 淫秽高清视频在线观看| 亚洲中文日韩欧美视频| 欧美性猛交黑人性爽| 久久久久久久久免费视频了| 久久精品国产亚洲av香蕉五月| 午夜福利成人在线免费观看| 成人三级黄色视频| xxxwww97欧美| 午夜视频精品福利| 男人的好看免费观看在线视频 | 九九热线精品视视频播放| 黄色丝袜av网址大全| 精品人妻1区二区| 国内精品久久久久久久电影| 国产三级中文精品| tocl精华| 精品不卡国产一区二区三区| 成人午夜高清在线视频| 亚洲自偷自拍图片 自拍| videosex国产| 亚洲精品av麻豆狂野| 99热只有精品国产| √禁漫天堂资源中文www| 欧美av亚洲av综合av国产av| 99热6这里只有精品| 午夜福利高清视频| 久久久久九九精品影院| 色老头精品视频在线观看| 亚洲国产欧洲综合997久久,| 久久伊人香网站| 美女午夜性视频免费| 国产97色在线日韩免费| 999久久久国产精品视频| 后天国语完整版免费观看| 国产精品98久久久久久宅男小说| 俄罗斯特黄特色一大片| 成人欧美大片| 1024香蕉在线观看| 精品久久久久久成人av| 国产真人三级小视频在线观看| 午夜精品在线福利| 夜夜看夜夜爽夜夜摸| 18美女黄网站色大片免费观看| 国产伦在线观看视频一区| 免费观看人在逋| 9191精品国产免费久久| 国产一区二区在线av高清观看| 午夜久久久久精精品| 日韩精品免费视频一区二区三区| 久久精品国产99精品国产亚洲性色| 亚洲欧美日韩无卡精品| 国产69精品久久久久777片 | 午夜福利在线在线| 好男人在线观看高清免费视频| 亚洲av五月六月丁香网| 最近最新中文字幕大全免费视频| 国产午夜精品论理片| 久久精品综合一区二区三区| 夜夜躁狠狠躁天天躁| 精品乱码久久久久久99久播| 国产av一区在线观看免费| 久久精品国产亚洲av香蕉五月| 久久精品国产99精品国产亚洲性色| av天堂在线播放| 午夜激情福利司机影院| 此物有八面人人有两片| 国内精品久久久久久久电影| 我的老师免费观看完整版| 又大又爽又粗| 国产精品亚洲一级av第二区| 又大又爽又粗| 国产视频一区二区在线看| 久久久久九九精品影院| 熟女电影av网| 午夜福利在线观看吧| 国产精品99久久99久久久不卡| 大型av网站在线播放| 又爽又黄无遮挡网站| 精品日产1卡2卡| 一区福利在线观看| 欧美国产日韩亚洲一区| 国产欧美日韩一区二区三| 黄色丝袜av网址大全| 免费电影在线观看免费观看| 国产在线精品亚洲第一网站| 美女高潮喷水抽搐中文字幕| 精品国内亚洲2022精品成人| 亚洲精品国产精品久久久不卡| 久久 成人 亚洲| 亚洲中文av在线| 在线观看舔阴道视频| 91av网站免费观看| 亚洲va日本ⅴa欧美va伊人久久| 国产成人精品久久二区二区91| 白带黄色成豆腐渣| 一边摸一边抽搐一进一小说| 麻豆成人av在线观看| 成人手机av| 麻豆久久精品国产亚洲av| 99在线视频只有这里精品首页| 久久久国产欧美日韩av| 久久人妻av系列| 亚洲精品粉嫩美女一区| 9191精品国产免费久久| 欧美日韩乱码在线| 亚洲av第一区精品v没综合| 99精品在免费线老司机午夜| 男女床上黄色一级片免费看| 精品第一国产精品| 免费高清视频大片| 高潮久久久久久久久久久不卡| 久久久水蜜桃国产精品网| 99国产精品99久久久久| 国产精品爽爽va在线观看网站| 看免费av毛片| 国产黄a三级三级三级人| 亚洲av片天天在线观看| 男女之事视频高清在线观看| 在线观看66精品国产| 久久久久九九精品影院| 欧美+亚洲+日韩+国产| 国产单亲对白刺激| 国产精品久久久人人做人人爽| 欧美日韩亚洲国产一区二区在线观看| 一进一出好大好爽视频| 日日摸夜夜添夜夜添小说| 我要搜黄色片| 国产精品国产高清国产av| 麻豆av在线久日| 午夜视频精品福利| 女人被狂操c到高潮| 亚洲国产看品久久| 天堂影院成人在线观看| 长腿黑丝高跟| 91九色精品人成在线观看| 久久精品国产亚洲av高清一级| 老司机在亚洲福利影院| 国产精品久久视频播放| 亚洲中文字幕一区二区三区有码在线看 | 天天躁夜夜躁狠狠躁躁| 成人国语在线视频| 久久天堂一区二区三区四区| 啦啦啦观看免费观看视频高清| or卡值多少钱| 久久久久国产精品人妻aⅴ院| 丝袜人妻中文字幕| 欧美又色又爽又黄视频| 成熟少妇高潮喷水视频| 久久久久久久久久黄片| 日日摸夜夜添夜夜添小说| 天天添夜夜摸| 色老头精品视频在线观看| 久久午夜亚洲精品久久| 91字幕亚洲| 老司机福利观看| 久久久久久久久久黄片| 99久久99久久久精品蜜桃| 免费人成视频x8x8入口观看| 精品国产亚洲在线| 非洲黑人性xxxx精品又粗又长| 少妇人妻一区二区三区视频| 欧美国产日韩亚洲一区| 男女下面进入的视频免费午夜| 黄色a级毛片大全视频|