• <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一区二区三区| 亚洲人成网站在线播| 禁无遮挡网站| 色尼玛亚洲综合影院| 国产亚洲精品综合一区在线观看| 国产精品综合久久久久久久免费| 久久久久久久久大av| 99久久无色码亚洲精品果冻| 直男gayav资源| 又黄又爽又免费观看的视频| 精华霜和精华液先用哪个| 露出奶头的视频| 男插女下体视频免费在线播放| 亚洲美女搞黄在线观看 | 老司机午夜十八禁免费视频| 久久久久久久亚洲中文字幕 | 久99久视频精品免费| 中文字幕人妻熟人妻熟丝袜美| 黄色一级大片看看| 久久婷婷人人爽人人干人人爱| 51午夜福利影视在线观看| 国产日本99.免费观看| 麻豆一二三区av精品| 悠悠久久av| 日韩欧美三级三区| 久久午夜亚洲精品久久| 可以在线观看毛片的网站| 成年女人看的毛片在线观看| 日本精品一区二区三区蜜桃| 精品熟女少妇八av免费久了| 精品午夜福利视频在线观看一区| 国产欧美日韩一区二区精品| 91av网一区二区| 欧美日韩国产亚洲二区| 免费人成在线观看视频色| 在线观看舔阴道视频| 网址你懂的国产日韩在线| 日本黄色视频三级网站网址| 欧美+亚洲+日韩+国产| 我要搜黄色片| 最后的刺客免费高清国语| 非洲黑人性xxxx精品又粗又长| 欧美另类亚洲清纯唯美| 免费人成在线观看视频色| 天天一区二区日本电影三级| 国产主播在线观看一区二区| 久久亚洲真实| 亚洲人成网站高清观看| 嫩草影视91久久| 国产三级黄色录像| 欧美日韩福利视频一区二区| 中文字幕av在线有码专区| 国内少妇人妻偷人精品xxx网站| АⅤ资源中文在线天堂| 无人区码免费观看不卡| 免费看光身美女| 97碰自拍视频| а√天堂www在线а√下载| 成年女人毛片免费观看观看9| 国产伦在线观看视频一区| 午夜福利在线在线| 亚洲欧美日韩卡通动漫| 欧美成人免费av一区二区三区| 黄色配什么色好看| 色噜噜av男人的天堂激情| av视频在线观看入口| 国产精品伦人一区二区| 成人高潮视频无遮挡免费网站| 搡老岳熟女国产| 搡老妇女老女人老熟妇| 人妻久久中文字幕网| 99久久无色码亚洲精品果冻| 真人一进一出gif抽搐免费| 免费观看精品视频网站| 日韩免费av在线播放| 国产成人福利小说| 最新在线观看一区二区三区| 白带黄色成豆腐渣| 精品乱码久久久久久99久播| 亚洲人成电影免费在线| 久久久久久久久久成人| 国产主播在线观看一区二区| 国产欧美日韩一区二区精品| 国产熟女xx| 美女大奶头视频| 岛国在线免费视频观看| 91麻豆av在线| 国产老妇女一区| 高清日韩中文字幕在线| 嫁个100分男人电影在线观看| 最近最新中文字幕大全电影3| 日本黄色视频三级网站网址| 天堂√8在线中文| 精品熟女少妇八av免费久了| 亚洲一区二区三区不卡视频| 精品国产三级普通话版| 高清毛片免费观看视频网站| 美女被艹到高潮喷水动态| 午夜福利高清视频| 51午夜福利影视在线观看| 成熟少妇高潮喷水视频| 永久网站在线| 国产亚洲欧美在线一区二区| 亚洲第一欧美日韩一区二区三区| 国产亚洲精品久久久久久毛片| 国产成年人精品一区二区| 日韩亚洲欧美综合| 亚洲,欧美,日韩| 少妇人妻一区二区三区视频| 18禁裸乳无遮挡免费网站照片| 亚洲人成网站高清观看| 日韩精品青青久久久久久| 欧美性感艳星| av黄色大香蕉| 国产爱豆传媒在线观看| av视频在线观看入口| 在线播放国产精品三级| 亚洲中文字幕一区二区三区有码在线看| 长腿黑丝高跟| 国产精品久久久久久久久免 | av在线观看视频网站免费| 亚洲av第一区精品v没综合| 亚洲美女搞黄在线观看 | 成人国产综合亚洲| 国产成年人精品一区二区| 久99久视频精品免费| 韩国av一区二区三区四区| 亚洲国产精品sss在线观看| 91久久精品国产一区二区成人| 久久久久精品国产欧美久久久| 国产一区二区三区视频了| 久久久久久久午夜电影| 一本久久中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 91狼人影院| 久久精品人妻少妇| 变态另类丝袜制服| 亚洲aⅴ乱码一区二区在线播放| 久久久国产成人精品二区| 亚洲一区二区三区不卡视频| av在线观看视频网站免费| 99国产综合亚洲精品| 国产成人欧美在线观看| 亚洲午夜理论影院| 神马国产精品三级电影在线观看| 一个人免费在线观看的高清视频| 久久久久九九精品影院| 国产亚洲欧美98| 久久伊人香网站| 免费av观看视频| 亚洲中文字幕一区二区三区有码在线看| 黄色女人牲交| 最近中文字幕高清免费大全6 | 亚洲七黄色美女视频| 桃色一区二区三区在线观看| 看免费av毛片| 亚洲 欧美 日韩 在线 免费| 少妇丰满av| 性色av乱码一区二区三区2| 免费在线观看亚洲国产| a级一级毛片免费在线观看| 动漫黄色视频在线观看| 好男人在线观看高清免费视频| 99热这里只有是精品在线观看 | 99国产极品粉嫩在线观看| 日本精品一区二区三区蜜桃| 在线观看午夜福利视频| 欧美另类亚洲清纯唯美| 熟妇人妻久久中文字幕3abv| 大型黄色视频在线免费观看| 国产精品日韩av在线免费观看| 真人一进一出gif抽搐免费| 国产精品一区二区免费欧美| 三级男女做爰猛烈吃奶摸视频| 少妇人妻一区二区三区视频| 偷拍熟女少妇极品色| 18禁黄网站禁片午夜丰满| 成人美女网站在线观看视频| av福利片在线观看| 婷婷六月久久综合丁香| 日本 欧美在线| 免费人成在线观看视频色| 亚洲欧美日韩无卡精品| 在线观看午夜福利视频| 一个人看视频在线观看www免费| 亚洲激情在线av| 熟女人妻精品中文字幕| 久久香蕉精品热| 淫秽高清视频在线观看| 国产又黄又爽又无遮挡在线| 人人妻人人澡欧美一区二区| 免费观看的影片在线观看| 亚洲午夜理论影院| 免费电影在线观看免费观看| 中出人妻视频一区二区| 成年女人看的毛片在线观看| 欧美在线一区亚洲| 欧美最黄视频在线播放免费| 久久亚洲精品不卡| 国产视频内射| 小说图片视频综合网站| 日本免费a在线| 国产大屁股一区二区在线视频| 亚洲av成人不卡在线观看播放网| 禁无遮挡网站| 白带黄色成豆腐渣| 国内毛片毛片毛片毛片毛片| 午夜久久久久精精品| 日韩精品青青久久久久久| 性插视频无遮挡在线免费观看| 国产av麻豆久久久久久久| 久久精品夜夜夜夜夜久久蜜豆| 久9热在线精品视频| 欧美绝顶高潮抽搐喷水| 90打野战视频偷拍视频| 十八禁国产超污无遮挡网站| 深夜a级毛片| 午夜福利视频1000在线观看| 亚洲一区二区三区不卡视频| 成年女人永久免费观看视频| 国产精品一区二区免费欧美| 99久久精品热视频| 亚洲人成伊人成综合网2020| 亚洲成人久久性| 欧美激情在线99| 看免费av毛片| 一区福利在线观看| 国产亚洲欧美98| 永久网站在线| 日韩欧美精品免费久久 | 久久久久国内视频| 亚洲国产精品合色在线| 偷拍熟女少妇极品色| 精品一区二区三区av网在线观看| 午夜影院日韩av| 日本一二三区视频观看| 免费看光身美女| 色哟哟·www| 丰满人妻熟妇乱又伦精品不卡| 国产av在哪里看| 国产亚洲精品久久久com| 国产免费男女视频| 国产亚洲欧美在线一区二区| 日韩欧美国产在线观看| 精品久久久久久久久亚洲 | 看黄色毛片网站| 亚洲综合色惰| 欧美3d第一页| 精品一区二区三区人妻视频| 久久6这里有精品| a级一级毛片免费在线观看| 丰满的人妻完整版| 中文字幕免费在线视频6| 大型黄色视频在线免费观看| 99在线视频只有这里精品首页| 精品福利观看| 久久久精品欧美日韩精品| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 俺也久久电影网| 精品久久久久久久人妻蜜臀av| 国内少妇人妻偷人精品xxx网站| 欧美在线一区亚洲| 天堂av国产一区二区熟女人妻| 亚洲av成人不卡在线观看播放网| 1024手机看黄色片| 永久网站在线| 久久亚洲真实| 亚洲,欧美,日韩| 国产一区二区激情短视频| 在线播放国产精品三级| 久久精品国产自在天天线| 国产中年淑女户外野战色| 夜夜看夜夜爽夜夜摸| 69人妻影院| 波多野结衣巨乳人妻| 免费电影在线观看免费观看| 在线播放国产精品三级| 婷婷六月久久综合丁香| 99热这里只有精品一区| 免费看a级黄色片| 精品欧美国产一区二区三| 波野结衣二区三区在线| 毛片一级片免费看久久久久 | 麻豆成人av在线观看| 国产伦人伦偷精品视频| 亚洲av中文字字幕乱码综合| 亚洲熟妇熟女久久| 夜夜躁狠狠躁天天躁| 国产高清激情床上av| 日本黄色片子视频| 国产男靠女视频免费网站| 国产精品嫩草影院av在线观看 | 亚洲 国产 在线| 久久性视频一级片| 一个人看的www免费观看视频| 天美传媒精品一区二区| 熟妇人妻久久中文字幕3abv| 超碰av人人做人人爽久久| 国产激情偷乱视频一区二区| 变态另类丝袜制服| 成人一区二区视频在线观看| 亚州av有码| 老司机福利观看| 国产午夜福利久久久久久| 亚洲av日韩精品久久久久久密| 久久久色成人| 99精品在免费线老司机午夜| 成人国产一区最新在线观看| 搡老妇女老女人老熟妇| 中出人妻视频一区二区| 在线观看免费视频日本深夜| 天堂√8在线中文| 又爽又黄无遮挡网站| 久久精品91蜜桃| 一个人免费在线观看的高清视频| 亚洲av电影不卡..在线观看| 午夜福利高清视频| 欧美激情在线99| 免费在线观看日本一区| 无遮挡黄片免费观看| 国产成人影院久久av| 日本一本二区三区精品| 99国产极品粉嫩在线观看| 一个人观看的视频www高清免费观看| 国产午夜精品论理片| 欧美黑人欧美精品刺激| 国产精品一区二区免费欧美| 亚洲自偷自拍三级| 亚洲av成人精品一区久久| 久久久久久久久大av| 色av中文字幕| 国产真实伦视频高清在线观看 | 免费av观看视频| 亚洲av日韩精品久久久久久密| 国产探花在线观看一区二区| 色哟哟·www| 免费观看精品视频网站| 亚洲av成人av| 亚洲av成人精品一区久久| 色吧在线观看| 久久精品国产亚洲av天美| 色噜噜av男人的天堂激情| 国产精品伦人一区二区| 欧美在线黄色| 亚洲第一欧美日韩一区二区三区| 无遮挡黄片免费观看| 99热这里只有是精品在线观看 | 亚洲国产精品久久男人天堂| 国产免费av片在线观看野外av| 免费黄网站久久成人精品 | 精品久久久久久久久久免费视频| 亚州av有码| 十八禁国产超污无遮挡网站| 一二三四社区在线视频社区8| 九色国产91popny在线| 色综合亚洲欧美另类图片| 国产久久久一区二区三区| 午夜福利欧美成人| 日本在线视频免费播放| 国产精品一及| 国产综合懂色| 免费电影在线观看免费观看| 成人美女网站在线观看视频| 变态另类丝袜制服| 国产伦精品一区二区三区四那| 亚洲av美国av| 99热这里只有是精品50| 亚洲内射少妇av| 窝窝影院91人妻| 有码 亚洲区| 丁香六月欧美| 国产精品日韩av在线免费观看| 看片在线看免费视频| 国产一区二区在线av高清观看| 人妻丰满熟妇av一区二区三区| 亚州av有码| 国产成人a区在线观看| 88av欧美| 国内精品久久久久精免费| 欧美精品国产亚洲| 免费一级毛片在线播放高清视频| 久久久久久久午夜电影| 91在线观看av| 国产精品精品国产色婷婷| 久久久久免费精品人妻一区二区| 超碰av人人做人人爽久久| 熟妇人妻久久中文字幕3abv| 狠狠狠狠99中文字幕| 国产精品一区二区三区四区免费观看 | 人妻制服诱惑在线中文字幕| 毛片一级片免费看久久久久 | 97碰自拍视频| 欧美色欧美亚洲另类二区| 2021天堂中文幕一二区在线观| 国产日本99.免费观看| 人人妻人人澡欧美一区二区| 成年免费大片在线观看| 搡老岳熟女国产| 国产v大片淫在线免费观看| 欧美bdsm另类| 免费黄网站久久成人精品 | 麻豆一二三区av精品| 亚洲最大成人av| 国产精品不卡视频一区二区 | 国产精品1区2区在线观看.| 欧美成人一区二区免费高清观看| 不卡一级毛片| 亚洲人成电影免费在线| 最好的美女福利视频网| 如何舔出高潮| 最新中文字幕久久久久| 99久久成人亚洲精品观看| 日韩欧美三级三区| 免费一级毛片在线播放高清视频| 免费看美女性在线毛片视频| 亚洲色图av天堂| 中国美女看黄片| 日韩中文字幕欧美一区二区| 亚洲人成电影免费在线| 亚洲乱码一区二区免费版| 国产精品久久久久久久久免 | 又爽又黄无遮挡网站| 久久精品国产自在天天线| 国产精品永久免费网站| 免费高清视频大片| 中文字幕av在线有码专区| 免费看美女性在线毛片视频| h日本视频在线播放| 久久久久九九精品影院| 精品人妻偷拍中文字幕| 免费看日本二区| 久久精品久久久久久噜噜老黄 | 永久网站在线| 国产麻豆成人av免费视频| 国产亚洲欧美在线一区二区| 在线播放无遮挡| 狠狠狠狠99中文字幕| av天堂中文字幕网| 好男人在线观看高清免费视频| 欧美+亚洲+日韩+国产| 久久久国产成人精品二区| 少妇丰满av| 国产aⅴ精品一区二区三区波| АⅤ资源中文在线天堂| 亚洲无线观看免费| 色精品久久人妻99蜜桃| 身体一侧抽搐| 99热6这里只有精品| 亚洲国产精品999在线| 国产久久久一区二区三区| 亚洲av免费在线观看| 亚洲在线观看片| 久久精品综合一区二区三区| 国产精品乱码一区二三区的特点| www.www免费av| 我的老师免费观看完整版| 狂野欧美白嫩少妇大欣赏| 少妇人妻一区二区三区视频| 午夜久久久久精精品| 国产成人福利小说| 日韩国内少妇激情av| 中亚洲国语对白在线视频| 999久久久精品免费观看国产| 熟女人妻精品中文字幕| 精华霜和精华液先用哪个| 精品久久久久久久久久久久久| 美女高潮喷水抽搐中文字幕| 91在线观看av| 中亚洲国语对白在线视频| 人妻夜夜爽99麻豆av| 亚洲成a人片在线一区二区| 欧美在线黄色| 久久久精品大字幕| 国产精品女同一区二区软件 | 国产精品影院久久| 国产麻豆成人av免费视频| 一区二区三区高清视频在线| 久久午夜福利片| 免费看光身美女| а√天堂www在线а√下载| 欧美区成人在线视频| 国产美女午夜福利| 男人舔女人下体高潮全视频| 啦啦啦观看免费观看视频高清| 噜噜噜噜噜久久久久久91| 伦理电影大哥的女人| 国产精品嫩草影院av在线观看 | 日韩欧美一区二区三区在线观看| 窝窝影院91人妻| 69人妻影院| 亚洲 欧美 日韩 在线 免费| 99热这里只有是精品在线观看 | 日本黄色片子视频| 他把我摸到了高潮在线观看| 中文字幕高清在线视频| 老鸭窝网址在线观看| 国产高潮美女av| 成年女人毛片免费观看观看9| 国产私拍福利视频在线观看| 小蜜桃在线观看免费完整版高清| 亚洲aⅴ乱码一区二区在线播放| 亚洲一区二区三区色噜噜| 色视频www国产| 高清日韩中文字幕在线| 国产精品不卡视频一区二区 | 91午夜精品亚洲一区二区三区 | av天堂中文字幕网| 国产美女午夜福利| 国产高清激情床上av| av专区在线播放| 亚洲av第一区精品v没综合| 欧美三级亚洲精品| 夜夜夜夜夜久久久久| 欧美午夜高清在线| 国产亚洲精品av在线| 欧美成人一区二区免费高清观看| 久久欧美精品欧美久久欧美| 黄色丝袜av网址大全| 90打野战视频偷拍视频| 一级黄色大片毛片| 非洲黑人性xxxx精品又粗又长| 久久这里只有精品中国| 看黄色毛片网站| 亚洲精华国产精华精| 久久精品人妻少妇| 欧美高清性xxxxhd video| 国产国拍精品亚洲av在线观看| 日韩国内少妇激情av| 欧美成狂野欧美在线观看| 国产午夜精品久久久久久一区二区三区 | 欧美精品国产亚洲| 在线十欧美十亚洲十日本专区| 久久九九热精品免费| 精品人妻视频免费看| 国产高清三级在线| 一夜夜www| 不卡一级毛片| 精品日产1卡2卡| 床上黄色一级片| 国产毛片a区久久久久| 精品熟女少妇八av免费久了| 亚洲片人在线观看| 国产伦精品一区二区三区四那| 亚洲黑人精品在线| 大型黄色视频在线免费观看| 日本在线视频免费播放| 噜噜噜噜噜久久久久久91| 国产老妇女一区| 国产亚洲欧美98| 美女高潮喷水抽搐中文字幕| 国产精品乱码一区二三区的特点| 麻豆成人午夜福利视频| 久久亚洲精品不卡| 9191精品国产免费久久| 亚洲欧美清纯卡通| 一级a爱片免费观看的视频| 精品国内亚洲2022精品成人| 日韩av在线大香蕉| 日韩欧美国产一区二区入口| 中文在线观看免费www的网站| 久久久精品欧美日韩精品| 亚洲中文字幕日韩| 成年女人毛片免费观看观看9| 十八禁网站免费在线| 国产精品一区二区三区四区免费观看 | 九色成人免费人妻av| 伦理电影大哥的女人| 中文字幕精品亚洲无线码一区| 欧美午夜高清在线| 熟妇人妻久久中文字幕3abv| 久9热在线精品视频| 精品国产三级普通话版| 99久久99久久久精品蜜桃| 亚洲五月婷婷丁香| 婷婷色综合大香蕉| a级一级毛片免费在线观看| 国模一区二区三区四区视频| 色在线成人网| 成人一区二区视频在线观看| 亚洲五月天丁香| av中文乱码字幕在线| 久久久精品大字幕| 久久亚洲真实| 欧美日本亚洲视频在线播放| 亚洲三级黄色毛片| 色精品久久人妻99蜜桃| 午夜影院日韩av| 脱女人内裤的视频| 欧美xxxx黑人xx丫x性爽| 色在线成人网| 窝窝影院91人妻| 久久久久国内视频| 91av网一区二区| 国产在线男女| 国产精华一区二区三区| 一a级毛片在线观看| 午夜日韩欧美国产| 毛片女人毛片| 久久这里只有精品中国| 欧美日韩国产亚洲二区| 99精品久久久久人妻精品| 久久婷婷人人爽人人干人人爱| 欧美日韩国产亚洲二区| 99精品久久久久人妻精品| 国产欧美日韩一区二区三| 老司机深夜福利视频在线观看| 国产美女午夜福利| 日韩欧美三级三区| 白带黄色成豆腐渣| 麻豆av噜噜一区二区三区| 99在线人妻在线中文字幕| 精品不卡国产一区二区三区| 国产精品久久久久久久久免 | 怎么达到女性高潮| 久久午夜福利片|