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

    層狀地基與彈性薄板相互作用的邊界元解*

    2017-05-15 03:30:22艾智勇蔡建邦
    湖南大學學報(自然科學版) 2017年3期
    關鍵詞:元法薄板層狀

    艾智勇,蔡建邦

    (1.同濟大學 地下建筑與工程系,上海 200092;2.同濟大學 巖土及地下工程教育部重點實驗室,上海 200092)

    ?

    層狀地基與彈性薄板相互作用的邊界元解*

    艾智勇1, 2?,蔡建邦1, 2

    (1.同濟大學 地下建筑與工程系,上海 200092;2.同濟大學 巖土及地下工程教育部重點實驗室,上海 200092)

    將無限大薄板的基本解作為薄板邊界積分方程的核函數(shù),對薄板的內部和邊界進行離散,并假定薄板內部和邊界上的節(jié)點與地基反力的分布情況,得到薄板的邊界元方程組;同時基于層狀地基的解析層元解,通過Guass-Legendre積分得到地基柔度矩陣;結合地基與薄板接觸面上的位移協(xié)調條件,得到層狀地基與薄板共同作用問題總的邊界元法方程組;求解該方程組,得到層狀地基與薄板共同作用問題的解答.基于本文理論,編制了相應的FORTRAN程序,通過與已有文獻結果對比驗證本文理論及程序的正確性,數(shù)值分析結果表明:方形基礎薄板情況下,離板中心越近,垂直于坐標軸y(x)方向、距離相等的2條線段的豎向位移差越小,且該位移差隨著板-土剛度比減小而減?。浑S著板長寬比的增大,板中心點與長邊中點位移差變化不明顯,而短邊中心與邊界角點的位移差也有相類似的規(guī)律.

    邊界元;層狀地基;薄板;解析層元

    筏板基礎具有剛度大、整體性好、能較好地抵抗不均勻沉降的優(yōu)點,因此在高層建筑中得到了廣泛的應用.目前,地基與板相互作用分析的方法主要有:有限差分法[1]、有限單元法[2-4]、邊界單元法[5]、邊界單元-有限單元耦合法[6]、廣義微分求積法[7]、半解析數(shù)值方法[8-10],以及有限網(wǎng)格法[11]等.相比于有限元、有限層等方法,邊界單元法能將求解過程的維數(shù)降低一維,并具有計算時間短、精度高等優(yōu)點.因此,很多學者運用邊界元法來研究筏板與地基的相互作用問題.佘穎禾和朱萬寧[12]將地基效應歸并到地基板的彎曲微分方程內,得到了含有第三類復變量的Bessel函數(shù)的基本解,再根據(jù)該問題的邊界積分方程,建立了Winkler和雙參數(shù)地基上薄板的無奇異邊界單元法.王建國和黃茂光[13]提出了雙參數(shù)地基上薄板問題的邊界單元解法.鄧安福等[14]采用邊界單元法研究了雙參數(shù)地基上的厚板問題.Rashed等[15]通過邊界元法研究了Winkler地基上的厚板問題.閆富有等[16]基于Reissner 板的邊界積分方程,建立了有限壓縮層地基上厚筏基礎與地基相互作用分析的邊界元法.

    由以上研究可知,目前板-土相互作用的邊界單元法研究所采用的地基模型大多是Winkler和雙參數(shù)地基模型.Winkler地基模型將地基對板的作用看做是一系列相互獨立的彈簧,忽略了彈簧之間的剪切作用,因此只適用于很軟弱的地基土.雙參數(shù)地基模型雖然在獨立彈簧之間引入力學的相互作用以消除其不連續(xù)性,但其參數(shù)較難獲取,因而限制了它的工程應用.天然地基由于沉積而常常呈層狀分布,因此采用層狀地基模型更加符合工程實際.而層狀地基上板的邊界元研究還很少見諸報道.為此本文對層狀地基與彈性薄板的共同作用問題進行邊界元分析,以便精確、高效地求解基礎板問題.

    1 彈性薄板的邊界積分方程

    考慮橫向分布荷載p時,彈性薄板的控制方程為:

    D22s=p.

    (1)

    式中:s為板中面的橫向位移;2為Laplace算子;)],為薄板的抗彎剛度,其中,hp,Ep和νp分別為板的厚度、彈性模量和泊松比.

    由式(1)可得無限大平面薄板(如圖1所示)的基本解為[17]:

    (2)

    (3)

    2(1+νp)(lnr+1)},

    (4)

    圖1 薄板邊界元法示意圖Fig.1 The diagram of the BEM for a thin plate

    采用由式(2)—(5)得到的基本解作為邊界積分方程的核函數(shù),于是可以得到具有光滑邊界的薄板撓度邊界積分方程:

    (6)

    式中:當源點在板內部時,即ξ1∈Ω時,c(ξ1)=1,當源點在板邊界時,即ξ1∈Γ時,c(ξ1)=0.5;Mn,Vn和θ分別表示薄板邊界的彎矩、橫向力、轉角;Γ和Ω分別表示薄板的邊界域和內域;α為nξ與x的夾角,θ為幅角,β=α-θ,其中,tξ和nξ分別為ξ2點處的切線方向和外法線方向.

    2 層狀地基與薄板的共同作用

    假設板四邊自由,有:

    Vn(ξ2)=0,

    (7)

    Mn(ξ2)=0.

    (8)

    結合地基反力分布假設和式(7)與式(8),當源點為內點及邊界點時,式(6)可分別表示為式(9)和式(10),即:

    (9)

    (10)

    式中:Γi(i=1,2,…,Ne)對應于各段邊界的單元;Ωi對應于反力qi(i=1,2,…,Nc)的作用面積;sc和se分別表示薄板內點、邊界點的豎向位移.

    圖2 層狀地基與薄板相互作用圖Fig.2 The draft of the interaction betweenlayered soils and a thin plate

    圖3 彈性薄板單元及節(jié)點圖Fig.3 The elements and nodes of an elastic thin plate

    當板上作用均布荷載p時,薄板的邊界積分方程可表示為:

    (11)

    其中,

    (12)

    (13)

    (14)

    (15)

    (16)

    (17)

    (18)

    (19)

    (20)

    (21)

    (22)

    (23)

    (24)

    式中:se,θ,sc,q分別為薄板邊界節(jié)點位移、邊界節(jié)點轉角、內節(jié)點位移、地基反力矩陣;Ic為元素是D的Nc階對角矩陣.

    根據(jù)層狀各向同性彈性地基的解析層元解[18],通過兩維Guass-Legendre積分可得單位均布荷載作用在反力qj區(qū)域時,板節(jié)點i處的沉降為:

    δij=?δ(i,ψ,R)dΓj.

    (25)

    式中:δ(i,ψ,R)為地基的基本解,它表示層狀地基表面任意點ψ作用單位點荷載在點i引起的沉降,R為點ψ到點i的距離.

    地基反力引起的板內點及邊界點的豎向位移可表示為:

    (26)

    (27)

    假設加載前后板與地基不脫離,則有:

    we=se,

    (28)

    wc=sc.

    (29)

    將式(26)-(29)代入式(11),可得:

    (30)

    求解式(30),得到地基與彈性薄板相互作用的板邊界節(jié)點轉角和地基反力,再通過式(26)和(27)可求得薄板內部點和邊界點的豎向位移.

    3 數(shù)值計算與分析

    3.1 理論與程序驗證

    為了驗證本文理論及程序的正確性,將本文的計算結果與Wang&Cheung[4]的結果進行對比(如圖4所示),其中,bp=4.0 m,νs=0.4,hp=0.2 m,Es=0.343×103MPa,νp=0.167,Ep=0.343×105MPa,均布荷載p=0.98 MPa(見圖5).本文采用厚度為1 000m的單層土來模擬彈性半空間地基.由圖5可知,本文的結果與Wang&Cheung[4]的結果吻合較好,這表明本文理論與程序的正確性.

    圖4 方形基礎薄板邊界的撓度對比圖Fig.4 The comparison of vertical displacementsof the boundary of a square foundation thin plate

    3.2 板-土剛度比的影響

    圖5 彈性半空間上的方形薄板Fig.5 A square thin plate on an elastic half-space

    圖6 方形薄板平面圖Fig.6 The plane graph of a square thin plate

    圖7 當K=0.05(實線)和0.01(虛線)時 各線段的豎向位移曲線Fig.7 The settlement curves of the lines on the plate when K=0.05(solid line) and 0.01(imaginary line)

    3.3 矩形板長寬比η的影響

    圖8 三層地基上矩形薄板Fig.8 A rectangular plate on three-layered soils

    η圖9 板邊界角點、短邊中點、長邊中點 和板中心點的撓度隨η的變化曲線Fig.9 The settlement-η curves of the boundaryangular dot, the midpoits of the long sideand the wide side, and the center of the plate

    4 結 論

    本文基于更加符合工程實際的層狀地基模型,采用邊界單元法來求解地基與彈性薄板的共同作用問題,并通過與已有文獻結果對比,驗證了本文理論與計算程序的正確性.數(shù)值分析結果表明:

    1)方形基礎薄板情況下,離板中心越近,垂直于y(x)方向、距離相等的2條線段的豎向位移差越小,且其隨板-土剛度比減小而減小.

    2)隨著長寬比的增大,短邊中點與長邊中點的撓度差也增大;但板中心點與長邊中點位移差變化不明顯,且短邊中心與邊界角點位移差也有相類似的規(guī)律.

    相比于有限元、有限層等方法,邊界單元法能將求解過程的維數(shù)降低一維,并具有計算時間短、精度高等優(yōu)點.另外,基于層狀地基的動力解析層元解[19]及達朗貝爾原理[9],還可進一步將本文工作拓展,用以分析層狀地基上彈性薄板的動力響應.

    [1] VALLABHAN C V G, STRAUGHAN W T, DAS Y C. Refined model for analysis of plates on elastic foundations[J]. Journal of Engineering Mechanics Division, ASCE, 1991, 117(12): 2830-2843.

    [2] CHEUNG Y K, NAG D K. Plates and beams on elastic foundations-linear and non-linear behaviour[J]. Geotechnique, 1968, 18(2): 250-260.

    [3] 王元漢, 邱先敏, 張佑啟. 彈性地基板的等參有限元法計算[J]. 巖土工程學報, 1998, 20(4): 10-14.

    WANG Yuanhan, QIU Xianmin, CHEUNG Y K. Plate on an elastic foundation calculated by isoparametric element methods[J]. Chinese Journal of Geotechnical Engineering, 1998, 20(4): 10-14. (In Chinese)

    [4] WANG Y H, CHEUNG Y K. Plate on cross-anisotropic foundation analyzed by the finite element method[J]. Computers and Geotechnics, 2001, 28(1): 37-54.

    [5] DE PAIVA J B, BUTTERFIELD R. Boundary element analysis of plate-soil interaction[J]. Computers & Structures, 1997, 64(1): 319-328.

    [6] 許金余, 吳彰春, 冷培義. 彈性地基板的無界邊界元-有限元耦合計算法[J]. 工程力學, 1994, 11(3): 137-143.

    XU Jinyu, WU Zhangchun, LENG Peiyi. A calculation method of coupling infinite BE and FE for plate based on elastic half space[J]. Engineering Mechanics, 1994, 11(3): 137-143. (In Chinese)

    |[7] NOBAKHTI S, AGHGAM M. Static analysis of rectangular thick plates resting on two-parameter elastic boundary strips[J]. European Journal of Mechanics—A/Solids, 2011, 30(3):442-448.

    [8] 李剛, 羅世平, 尚守平. Vlazov地基上Reissner板相互作用的加權殘值分析法[J]. 湖南大學學報:自然科學版, 2006, 33(1): 20-24.

    LI Gang, LUO Shiping, SHANG Shouping. Weighted residual method for the analyzing of the interaction of Reissner’s plate on Vlazov’s elastic foundation[J]. Journal of Hunan University: Natural Sciences, 2006, 33(1): 20-24.(In Chinese)

    [9] 李剛, 丁翔, 尚守平. 有限條法分析雙參數(shù)地基板的動靜響應[J]. 湖南大學學報:自然科學版, 2008,35(4): 7-10.

    LI Gang, DING Xiang, SHANG Shouping. Finite stripe method for the static and dynamic analysis of plates on two-parameter elastic foundation[J]. Journal of Hunan University: Natural Sciences, 2008, 35(4): 7-10.(In Chinese)

    [10]LIANG X, WANG Z, WANG L,etal. Semi-analytical solution for three-dimensional transient response of functionally graded annular plate on a two parameter viscoelastic foundation[J]. Journal of Sound & Vibration, 2014, 333(12): 2649-2663.

    [11]KARASIN H, GULKAN P, AKTAS G. A finite grid solution for circular plates on elastic foundations[J]. KSCE Journal of Civil Engineering, 2015,19(4):1157-1163.

    [12]佘穎禾, 朱萬寧. 彈性地基板的無奇異邊界元解法[J]. 計算結構力學及其應用, 1991, 8(3): 267-274.

    SHE Yinghe, ZHU Wanning. A boundary element solution for the bending problem of the thin plates on elastic foundation[J]. Computational Structure Mehanics and Applications, 1991, 8(3): 267-274. (In Chinese)

    [13]王建國, 黃茂光. 兩參數(shù)彈性地基板的邊界元分析[J]. 土木工程學報, 1992, 25(3): 51-59.

    WANG Jianguo, HUANG Maoguang. Boundary element analysis of thin plates on the two-parameter elstic foundation[J]. China Civil Engineering Journal, 1992, 25(3): 51-59. (In Chinese)

    [14]鄧安福, 周永明, 李正良, 等. 雙參數(shù)彈性地基厚板分析的邊界元法[J]. 重慶建筑大學學報, 1997, 19(2): 6-16.

    DENG Anfu, ZHOU Yongming, LI Zhengliang,etal. Boundary element method for analysis of thick plates on two-parameter elastic foundations[J]. Journal of Chongqing Jianzhu University, 1997, 19(2): 6-16. (In Chinese)

    [15]RASHED Y F, ALIABADI M H, BREBBIA C A. The boundary element method for thick plates on a Winkler foundation[J]. International Journal for Numerical Methods in Engineering, 1998, 41(8): 1435-1462.

    [16]閆富有, 吳義章, 郭院成. 有限壓縮層地基上筏基沉降的邊界元分析[J]. 巖土力學,2011, 32(S2): 604-609.

    YAN Fuyou, WU Yizhang, GUO Yuancheng. Boundary element analysis of settlements of raft on finite compressible strata[J]. Rock and Soil Mechanics, 2011, 32(S2): 604-609. (In Chinese)

    [17]杜慶華. 邊界積分方程方法——邊界元法[M]. 北京: 高等教育出版社, 1989:135-137.

    DU Qinghua. Boundary integral equation method—boundary element method[M]. Beijing: Higher Education Press, 1989: 135-137.(In Chinese)

    [18]艾智勇, 董洲, 成怡沖. 多層地基軸對稱彈性空間問題的解析層元解[J]. 建筑結構學報, 2012, 44(4): 150-153.

    AI Zhiyong, DONG Zhou, CHENG Yichong. Analytical layer element solution of axisymmetrically elastic problem for multilayerd soils[J]. Journal of Building Structures, 2012,44(4): 150-153. (In Chinese)

    [19]AI Z Y, LI Z X, CANG N R. Analytical layer-element solution to axisymmetric dynamic response of transversely isotropic multilayered half-space[J]. Soil Dynamics and Earthquake Engineering, 2014, 60(5): 22-30.

    Foundations and an Elastic Thin Plate

    AI Zhiyong1, 2?, CAI Jianbang1, 2

    (1. Department of Geotechnical Engineering, Tongji University, Shanghai 200092, China; 2. Key Laboratoryof Geotechnical and Underground Engineering of Ministry of Education, Tongji University, Shanghai 200092, China)

    The kernel functions of the boundary integral equations for thin plate were determined by the fundamental solutions for an infinite thin plate. By the discretization of the plate interior and boundary as well as the assumption of the distribution states of plate nodes and foundation reaction forces, the BEM equations of the plate can be established. Meanwhile, based on the analytical layer element solutions for layered foundations, the flexibility matrix of the foundation was obtained by a two-dimensioned Guass-Legendre quadrature. Taking into account the compatible conditions of the displacements at the soils-plate interface, the global BEM equations for the interaction problem between the layered foundation and the thin plate were then established. The solutions for the problem were further obtained by solving the global BEM equations. The accuracy of the present method was verified by comparing existing solutions with the numerical results obtained from the corresponding FORTRAN program in this study. It is observed from numerical examples that when a square thin plate is placed on a foundation, the settlement difference between the two lines perpendicular to y or x coordinate decreases as they approach the center of the plate, and the difference decreases with the decrease of the plate-soil stiffness ratio. Furthermore, the settlement discrepancy between the plate center and the midpoint of the long side is unapparent with the increasing length-width ratio, and the similar variation trend can be found between the midpoint of the wide side and angular point.

    boundary element; layered soils; thin plates; analytical layer element

    2015-11-26

    國家自然科學基金資助項目(50578121), National Natural Science Foundation of China(50578121)

    艾智勇(1966-),男,江西余江人,同濟大學教授,博士 ?通訊聯(lián)系人,E-mail:zhiyongai@#edu.cn

    1674-2974(2017)03-0120-06

    10.16339/j.cnki.hdxbzkb.2017.03.015

    TU443

    AA BEM for Interaction between Layered

    猜你喜歡
    元法薄板層狀
    一角點支撐另一對邊固支正交各向異性矩形薄板彎曲的辛疊加解
    換元法在解題中的運用
    10MN鋁合金薄板拉伸機組的研制
    重型機械(2020年2期)2020-07-24 08:16:12
    基于離散元法的礦石對溜槽沖擊力的模擬研究
    重型機械(2019年3期)2019-08-27 00:58:46
    軋制復合制備TA1/AZ31B/TA1層狀復合材料組織與性能研究
    重型機械(2019年3期)2019-08-27 00:58:44
    換元法在解題中的應用
    “微元法”在含電容器電路中的應用
    兩級結構層狀Ti-TiBw/Ti復合材料擴散焊及其拉伸行為
    焊接(2016年9期)2016-02-27 13:05:22
    鋁薄板高速DP-GMAW焊接性能的研究
    焊接(2016年5期)2016-02-27 13:04:42
    高韌性抗層狀撕裂Q345FTE-Z35鋼板開發(fā)
    新疆鋼鐵(2015年2期)2015-11-07 03:27:52
    日日摸夜夜添夜夜添小说| 国产男靠女视频免费网站| 亚洲美女视频黄频| 女人十人毛片免费观看3o分钟| 老司机午夜十八禁免费视频| 夜夜躁狠狠躁天天躁| 亚洲av五月六月丁香网| www日本黄色视频网| 俄罗斯特黄特色一大片| 国产日本99.免费观看| 精品国内亚洲2022精品成人| 成年女人看的毛片在线观看| 一二三四社区在线视频社区8| 丝袜美腿在线中文| 色尼玛亚洲综合影院| 老司机深夜福利视频在线观看| 嫩草影视91久久| 欧美黑人欧美精品刺激| av女优亚洲男人天堂| 久久人妻av系列| 国产精品1区2区在线观看.| 国产一区二区在线观看日韩 | 在线观看免费午夜福利视频| 国产激情偷乱视频一区二区| 日韩欧美国产一区二区入口| 免费大片18禁| 91麻豆av在线| 色老头精品视频在线观看| 欧美av亚洲av综合av国产av| 久久6这里有精品| 国产单亲对白刺激| 小蜜桃在线观看免费完整版高清| 在线观看av片永久免费下载| 亚洲精品国产精品久久久不卡| 毛片女人毛片| 精品人妻一区二区三区麻豆 | 麻豆国产97在线/欧美| 男插女下体视频免费在线播放| 尤物成人国产欧美一区二区三区| 亚洲色图av天堂| 亚洲成人精品中文字幕电影| 欧美日韩乱码在线| 亚洲av不卡在线观看| 欧美黄色片欧美黄色片| 久久6这里有精品| 欧美一级毛片孕妇| 黄色成人免费大全| 女警被强在线播放| 欧美日韩国产亚洲二区| 国产真实伦视频高清在线观看 | 老汉色av国产亚洲站长工具| 给我免费播放毛片高清在线观看| 国内精品久久久久精免费| 亚洲国产高清在线一区二区三| 九九久久精品国产亚洲av麻豆| 丝袜美腿在线中文| 精品99又大又爽又粗少妇毛片 | 欧美日本视频| 亚洲精品亚洲一区二区| 日本免费a在线| 色视频www国产| 国产亚洲精品av在线| 日韩精品青青久久久久久| 19禁男女啪啪无遮挡网站| 国产精品嫩草影院av在线观看 | 1024手机看黄色片| 在线国产一区二区在线| 神马国产精品三级电影在线观看| 极品教师在线免费播放| 校园春色视频在线观看| 国产视频一区二区在线看| 一进一出好大好爽视频| 最后的刺客免费高清国语| 99久久成人亚洲精品观看| 欧美日韩综合久久久久久 | 国产亚洲精品久久久久久毛片| 亚洲美女视频黄频| 男人和女人高潮做爰伦理| 天堂网av新在线| 亚洲国产高清在线一区二区三| 亚洲18禁久久av| 麻豆国产av国片精品| 99热这里只有是精品50| 此物有八面人人有两片| 欧美性猛交╳xxx乱大交人| 免费在线观看亚洲国产| 国产精品99久久99久久久不卡| 黄色女人牲交| 日韩亚洲欧美综合| 真人做人爱边吃奶动态| 亚洲乱码一区二区免费版| 欧美一区二区亚洲| 精品久久久久久成人av| 成人特级黄色片久久久久久久| 高清日韩中文字幕在线| 国产淫片久久久久久久久 | 9191精品国产免费久久| 搡老岳熟女国产| 女人十人毛片免费观看3o分钟| 欧美日韩亚洲国产一区二区在线观看| 欧美日韩瑟瑟在线播放| 国产欧美日韩一区二区三| www.www免费av| 亚洲天堂国产精品一区在线| 久久久久久人人人人人| 观看美女的网站| 老司机午夜福利在线观看视频| 国产精品乱码一区二三区的特点| 国产免费av片在线观看野外av| 国产精品久久视频播放| 啦啦啦观看免费观看视频高清| 麻豆久久精品国产亚洲av| 国产成年人精品一区二区| 亚洲欧美日韩高清专用| 欧美日本视频| 久久婷婷人人爽人人干人人爱| 精品久久久久久久末码| 美女被艹到高潮喷水动态| 婷婷亚洲欧美| 18禁国产床啪视频网站| 三级毛片av免费| 国产亚洲精品一区二区www| 少妇人妻一区二区三区视频| 18禁裸乳无遮挡免费网站照片| 亚洲人成网站高清观看| 人人妻人人澡欧美一区二区| 亚洲黑人精品在线| 免费观看精品视频网站| 国产精品久久久久久亚洲av鲁大| 啪啪无遮挡十八禁网站| 日本一本二区三区精品| 一区二区三区高清视频在线| 日本免费a在线| 欧美另类亚洲清纯唯美| 国产午夜福利久久久久久| 变态另类丝袜制服| 亚洲第一欧美日韩一区二区三区| 动漫黄色视频在线观看| 听说在线观看完整版免费高清| 国产单亲对白刺激| 熟女电影av网| 丰满乱子伦码专区| www国产在线视频色| 国产精品1区2区在线观看.| 97超视频在线观看视频| 欧美色欧美亚洲另类二区| 少妇高潮的动态图| www.999成人在线观看| 欧美日韩一级在线毛片| 在线十欧美十亚洲十日本专区| 九九久久精品国产亚洲av麻豆| 久久九九热精品免费| 欧美最黄视频在线播放免费| 丁香欧美五月| 毛片女人毛片| 好看av亚洲va欧美ⅴa在| 一级a爱片免费观看的视频| 色精品久久人妻99蜜桃| 真人做人爱边吃奶动态| 国产精品av视频在线免费观看| 在线观看午夜福利视频| 午夜日韩欧美国产| 国产免费av片在线观看野外av| 精品欧美国产一区二区三| 国产伦精品一区二区三区四那| 村上凉子中文字幕在线| 一本久久中文字幕| 中文资源天堂在线| 高清日韩中文字幕在线| 老司机午夜十八禁免费视频| 99视频精品全部免费 在线| 国产av不卡久久| 首页视频小说图片口味搜索| 亚洲人成网站在线播放欧美日韩| 久久香蕉国产精品| 人妻久久中文字幕网| 精品不卡国产一区二区三区| 丰满乱子伦码专区| xxx96com| 美女免费视频网站| 麻豆成人午夜福利视频| 国产精品国产高清国产av| 久久久久久久久大av| 国产真实乱freesex| 在线观看舔阴道视频| 男人舔女人下体高潮全视频| 国产成人aa在线观看| 无限看片的www在线观看| 国产蜜桃级精品一区二区三区| 麻豆国产av国片精品| 国产伦在线观看视频一区| 婷婷六月久久综合丁香| x7x7x7水蜜桃| 国产精品亚洲一级av第二区| 乱人视频在线观看| 亚洲久久久久久中文字幕| 亚洲一区二区三区色噜噜| 亚洲av电影在线进入| 精品99又大又爽又粗少妇毛片 | 精品不卡国产一区二区三区| 听说在线观看完整版免费高清| 日本三级黄在线观看| 两个人的视频大全免费| 免费搜索国产男女视频| 免费一级毛片在线播放高清视频| 五月伊人婷婷丁香| 亚洲人成电影免费在线| 亚洲精品一卡2卡三卡4卡5卡| 国产av不卡久久| 老司机午夜福利在线观看视频| 亚洲中文字幕日韩| 精品免费久久久久久久清纯| 午夜精品一区二区三区免费看| 最好的美女福利视频网| 亚洲最大成人手机在线| 国产精品亚洲美女久久久| 精品国产美女av久久久久小说| 午夜激情欧美在线| 日本黄色视频三级网站网址| 观看美女的网站| 日本一本二区三区精品| 3wmmmm亚洲av在线观看| 少妇人妻精品综合一区二区 | 国产午夜福利久久久久久| 少妇裸体淫交视频免费看高清| 国产在视频线在精品| 97碰自拍视频| 男插女下体视频免费在线播放| 香蕉久久夜色| 亚洲精品粉嫩美女一区| 天堂影院成人在线观看| av国产免费在线观看| 真人一进一出gif抽搐免费| 搞女人的毛片| 久久久久久国产a免费观看| 法律面前人人平等表现在哪些方面| 一进一出抽搐gif免费好疼| 国内毛片毛片毛片毛片毛片| 欧美一级a爱片免费观看看| 国产午夜精品久久久久久一区二区三区 | 成年人黄色毛片网站| 免费av观看视频| 亚洲国产精品999在线| 日本一二三区视频观看| 99久久99久久久精品蜜桃| 国产私拍福利视频在线观看| 久久草成人影院| 国产爱豆传媒在线观看| 色播亚洲综合网| 美女被艹到高潮喷水动态| 欧美bdsm另类| 久99久视频精品免费| 老熟妇乱子伦视频在线观看| 俄罗斯特黄特色一大片| 看免费av毛片| 午夜免费观看网址| 99精品久久久久人妻精品| 99国产精品一区二区三区| 男女床上黄色一级片免费看| 日本在线视频免费播放| 搡老岳熟女国产| 99精品久久久久人妻精品| av福利片在线观看| 超碰av人人做人人爽久久 | 热99在线观看视频| 99热6这里只有精品| 日本黄色视频三级网站网址| 国产精品久久视频播放| 国产高清有码在线观看视频| 少妇的逼水好多| 亚洲激情在线av| 久久久国产成人精品二区| 日韩欧美 国产精品| 免费搜索国产男女视频| 99国产精品一区二区三区| 国产一区二区激情短视频| svipshipincom国产片| 男女午夜视频在线观看| 午夜激情欧美在线| 在线观看av片永久免费下载| av女优亚洲男人天堂| 变态另类成人亚洲欧美熟女| 一个人观看的视频www高清免费观看| 一区二区三区国产精品乱码| 变态另类丝袜制服| 夜夜躁狠狠躁天天躁| 热99在线观看视频| 色av中文字幕| 亚洲国产精品999在线| 一进一出抽搐动态| 每晚都被弄得嗷嗷叫到高潮| 欧美成人一区二区免费高清观看| 国产毛片a区久久久久| 国产一区二区三区在线臀色熟女| 国产av麻豆久久久久久久| 丰满的人妻完整版| 亚洲人成网站在线播| 女人高潮潮喷娇喘18禁视频| 99热6这里只有精品| 悠悠久久av| 亚洲精品久久国产高清桃花| 国产精品99久久99久久久不卡| 国产精品女同一区二区软件 | 久久久久精品国产欧美久久久| 最好的美女福利视频网| 在线看三级毛片| 亚洲精品在线美女| 中国美女看黄片| 一级毛片女人18水好多| 日韩欧美精品v在线| 久久久色成人| 婷婷精品国产亚洲av| 精品乱码久久久久久99久播| 精品99又大又爽又粗少妇毛片 | 日本免费一区二区三区高清不卡| 久久午夜亚洲精品久久| 精品国产亚洲在线| 久久久久国内视频| av福利片在线观看| 日本 欧美在线| 内射极品少妇av片p| 久久久久久国产a免费观看| 真实男女啪啪啪动态图| 国产免费av片在线观看野外av| 国产单亲对白刺激| 少妇人妻一区二区三区视频| 午夜激情福利司机影院| 波多野结衣高清无吗| 亚洲人成电影免费在线| 日韩免费av在线播放| 黄片大片在线免费观看| 在线观看免费视频日本深夜| 老鸭窝网址在线观看| 99国产精品一区二区三区| 嫩草影院精品99| 99精品在免费线老司机午夜| 国产欧美日韩精品亚洲av| 国产色爽女视频免费观看| 久久精品国产亚洲av香蕉五月| 老司机午夜福利在线观看视频| 亚洲真实伦在线观看| 黄色片一级片一级黄色片| 精品一区二区三区人妻视频| 欧美乱妇无乱码| 一区二区三区激情视频| 桃红色精品国产亚洲av| 国产午夜福利久久久久久| 桃红色精品国产亚洲av| 亚洲电影在线观看av| 亚洲五月天丁香| 99久久九九国产精品国产免费| 深爱激情五月婷婷| 亚洲国产日韩欧美精品在线观看 | 国产精品99久久久久久久久| 男女之事视频高清在线观看| 岛国在线观看网站| 亚洲欧美日韩卡通动漫| 一个人看视频在线观看www免费 | 中文亚洲av片在线观看爽| 亚洲五月天丁香| 老熟妇乱子伦视频在线观看| 国产高清有码在线观看视频| 手机成人av网站| 天美传媒精品一区二区| 亚洲欧美日韩东京热| 听说在线观看完整版免费高清| 中文字幕av在线有码专区| 天堂√8在线中文| 国产aⅴ精品一区二区三区波| 亚洲片人在线观看| 国产精品永久免费网站| 黄色女人牲交| 精品人妻一区二区三区麻豆 | 欧美成人免费av一区二区三区| 在线视频色国产色| 欧美成人免费av一区二区三区| av国产免费在线观看| 久久精品国产99精品国产亚洲性色| av中文乱码字幕在线| 无人区码免费观看不卡| 久久欧美精品欧美久久欧美| 99久久精品国产亚洲精品| 男人舔奶头视频| 久9热在线精品视频| 身体一侧抽搐| 可以在线观看的亚洲视频| 久久久久久人人人人人| 在线看三级毛片| 国产精品影院久久| 欧美中文日本在线观看视频| 少妇人妻一区二区三区视频| 色av中文字幕| 黄色女人牲交| 两个人看的免费小视频| 中文字幕人妻丝袜一区二区| 国产精品电影一区二区三区| 综合色av麻豆| 搞女人的毛片| 欧美黑人巨大hd| 亚洲色图av天堂| 男女视频在线观看网站免费| 在线a可以看的网站| 嫩草影院入口| 成人特级黄色片久久久久久久| 午夜免费激情av| 亚洲国产高清在线一区二区三| 日本一二三区视频观看| 麻豆久久精品国产亚洲av| 99精品久久久久人妻精品| 又粗又爽又猛毛片免费看| 狂野欧美白嫩少妇大欣赏| 此物有八面人人有两片| 少妇的丰满在线观看| 亚洲精品成人久久久久久| 蜜桃亚洲精品一区二区三区| 国产主播在线观看一区二区| 九色国产91popny在线| 深夜精品福利| 麻豆一二三区av精品| 色综合站精品国产| 亚洲国产日韩欧美精品在线观看 | 午夜两性在线视频| av女优亚洲男人天堂| 嫩草影视91久久| 欧美色欧美亚洲另类二区| 美女大奶头视频| 亚洲av电影不卡..在线观看| 久久99热这里只有精品18| 99久久成人亚洲精品观看| 深夜精品福利| 久久久久久久久中文| 99国产综合亚洲精品| 特级一级黄色大片| 嫩草影院入口| 日本黄大片高清| 久久久国产成人免费| 国产精品久久久久久亚洲av鲁大| 国产v大片淫在线免费观看| 久久精品夜夜夜夜夜久久蜜豆| 全区人妻精品视频| 国产精品久久久久久人妻精品电影| 在线观看66精品国产| 亚洲五月婷婷丁香| 亚洲专区中文字幕在线| 超碰av人人做人人爽久久 | 国产一区二区三区在线臀色熟女| 成人高潮视频无遮挡免费网站| 男女视频在线观看网站免费| 色哟哟哟哟哟哟| 国产淫片久久久久久久久 | 一进一出好大好爽视频| 国产激情欧美一区二区| 深爱激情五月婷婷| 精品福利观看| 真人一进一出gif抽搐免费| 久久精品人妻少妇| 99久久综合精品五月天人人| 在线免费观看的www视频| 看免费av毛片| 最新中文字幕久久久久| 亚洲不卡免费看| 欧美黑人欧美精品刺激| 亚洲国产欧美网| 狠狠狠狠99中文字幕| 日本成人三级电影网站| 久9热在线精品视频| 99国产极品粉嫩在线观看| 久久久国产成人精品二区| 婷婷六月久久综合丁香| 国产色婷婷99| av在线天堂中文字幕| 一本综合久久免费| 嫩草影院入口| 日韩欧美 国产精品| 两人在一起打扑克的视频| 男人舔女人下体高潮全视频| 淫妇啪啪啪对白视频| 中文字幕熟女人妻在线| 亚洲色图av天堂| 国产精品嫩草影院av在线观看 | 亚洲国产欧美网| 日日干狠狠操夜夜爽| 国产老妇女一区| 亚洲人成电影免费在线| 亚洲成人中文字幕在线播放| 两个人看的免费小视频| 噜噜噜噜噜久久久久久91| 岛国在线观看网站| 亚洲av免费高清在线观看| 久久久久亚洲av毛片大全| 欧美一级毛片孕妇| 宅男免费午夜| 一区二区三区国产精品乱码| 69av精品久久久久久| 熟女电影av网| 亚洲欧美激情综合另类| 国内揄拍国产精品人妻在线| 亚洲在线观看片| 国产精品爽爽va在线观看网站| 日韩国内少妇激情av| 国产成人影院久久av| 嫩草影视91久久| 校园春色视频在线观看| 波多野结衣高清作品| 日韩欧美在线乱码| 日韩av在线大香蕉| 久久精品91蜜桃| 精品无人区乱码1区二区| 精品国产三级普通话版| 国产精品野战在线观看| 免费搜索国产男女视频| 日韩免费av在线播放| 午夜福利成人在线免费观看| 亚洲国产欧洲综合997久久,| 一个人看视频在线观看www免费 | 欧美日本亚洲视频在线播放| 黄色丝袜av网址大全| 国产一区二区激情短视频| 久久久久国产精品人妻aⅴ院| 久久久色成人| www日本黄色视频网| 亚洲欧美日韩高清在线视频| 我要搜黄色片| 成年女人看的毛片在线观看| 国产单亲对白刺激| 中文字幕人妻丝袜一区二区| 伊人久久精品亚洲午夜| 亚洲成av人片免费观看| 亚洲自拍偷在线| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产精品女同一区二区软件 | 亚洲欧美激情综合另类| 2021天堂中文幕一二区在线观| 日韩成人在线观看一区二区三区| 少妇人妻精品综合一区二区 | 亚洲av免费高清在线观看| 九九在线视频观看精品| 色精品久久人妻99蜜桃| 亚洲av二区三区四区| 欧美成人a在线观看| 免费av不卡在线播放| 桃红色精品国产亚洲av| 人妻夜夜爽99麻豆av| 在线观看免费午夜福利视频| 99国产极品粉嫩在线观看| 性欧美人与动物交配| 啦啦啦免费观看视频1| 精品久久久久久成人av| 国产免费av片在线观看野外av| 国产 一区 欧美 日韩| 88av欧美| 成人欧美大片| 丰满乱子伦码专区| 久久久久免费精品人妻一区二区| 69人妻影院| 九九在线视频观看精品| 日韩亚洲欧美综合| 国产三级中文精品| 欧美色欧美亚洲另类二区| 成人三级黄色视频| 校园春色视频在线观看| 国产免费av片在线观看野外av| 男女那种视频在线观看| 成年女人永久免费观看视频| 欧美日本视频| 狠狠狠狠99中文字幕| 一级黄色大片毛片| 给我免费播放毛片高清在线观看| 国产乱人视频| 又爽又黄无遮挡网站| 久久99热这里只有精品18| 久久久精品大字幕| 亚洲狠狠婷婷综合久久图片| 在线十欧美十亚洲十日本专区| 欧美xxxx黑人xx丫x性爽| 人妻夜夜爽99麻豆av| 少妇人妻精品综合一区二区 | 精品一区二区三区av网在线观看| 亚洲性夜色夜夜综合| 九九在线视频观看精品| 性欧美人与动物交配| АⅤ资源中文在线天堂| 国产黄片美女视频| 午夜激情欧美在线| 精品乱码久久久久久99久播| 日本熟妇午夜| bbb黄色大片| 午夜免费观看网址| 老司机深夜福利视频在线观看| av女优亚洲男人天堂| 国产蜜桃级精品一区二区三区| 欧美乱码精品一区二区三区| 给我免费播放毛片高清在线观看| 中文字幕人成人乱码亚洲影| 国产精品野战在线观看| 99riav亚洲国产免费| 村上凉子中文字幕在线| 国产探花在线观看一区二区| 久久亚洲真实| 亚洲美女黄片视频| 国产探花在线观看一区二区| 午夜福利在线观看吧| 国产野战对白在线观看| 欧美一级a爱片免费观看看| 母亲3免费完整高清在线观看| 国产精品久久久久久亚洲av鲁大| 两个人看的免费小视频| 91久久精品国产一区二区成人 | 国产真实乱freesex| 日韩欧美一区二区三区在线观看| 国产免费男女视频| 成人av在线播放网站| 91麻豆精品激情在线观看国产| 免费看光身美女| 热99re8久久精品国产|