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

    黏彈性邊界的二次開(kāi)發(fā)及其在地下結(jié)構(gòu)抗震分析中的應(yīng)用

    2019-09-10 10:46:48竇遠(yuǎn)明范俊超王建寧鞠培東宋明軒李景文
    關(guān)鍵詞:二次開(kāi)發(fā)結(jié)點(diǎn)隧洞

    竇遠(yuǎn)明 范俊超 王建寧 鞠培東 宋明軒 李景文

    摘要 在對(duì)地下結(jié)構(gòu)進(jìn)行抗震分析時(shí),土體邊界條件和地震波的施加方法直接關(guān)系到運(yùn)算結(jié)果的精準(zhǔn)程度。為了使地下結(jié)構(gòu)抗震分析建模更加高效,分析結(jié)果更加合理,對(duì)邊界和地震波的施加算法進(jìn)行了程序化設(shè)計(jì),利用Python語(yǔ)言對(duì)ABAQUS進(jìn)行了二次開(kāi)發(fā),編寫了黏彈性邊界和地震波統(tǒng)一自動(dòng)施加程序,建立了土體—隧洞結(jié)構(gòu)相互作用的三維有限元分析模型。結(jié)果表明:該方法可以實(shí)現(xiàn)黏彈性邊界和地震波的快速自動(dòng)施加,能夠很好地模擬波動(dòng)在土體中的傳播規(guī)律;在靠近土體邊界附近一定范圍內(nèi)的加速度峰值有3%左右的誤差,當(dāng)模型尺寸取9倍的結(jié)構(gòu)寬度時(shí)可以消除這一影響;隧洞結(jié)構(gòu)縱向端部2~3倍結(jié)構(gòu)寬度范圍內(nèi)的計(jì)算結(jié)果偏大。

    關(guān) 鍵 詞 地下結(jié)構(gòu);黏彈性邊界;二次開(kāi)發(fā);有限元分析

    中圖分類號(hào) TU93? ? ?文獻(xiàn)標(biāo)志碼 A

    Abstract In the seismic analysis of underground structure, the boundary conditions of soil and the seismic wave application methods are directly related to the accuracy of calculation results. In order to make the seismic analysis more effictive and more reasonable, the application algorithm of boundary and seismic waves was programmed and a secondary development of ABAQUS with Python was done. A finite element analysis model with 3D nonlinear soil-tunnel interaction was established by using the unified automatic impose program for viscoelastic boundary and seismic waves. The results show that this method can apply the viscoelastic boundary and the seismic waves automatically and can reflect the seismic law of foundation reasonably. There is a deviation in peak acceleration deviation (about 3%) near the model border and it can be eliminated when the FEM size is 9 times larger than the width of structure. The calculation results in the scope of 2D~3D (D represents the width of structure) width of the tunnel structure end are deviated to be larger.

    Key words underground structure; viscoelastic boundary; secondary exploration; finite element analysis

    0 引言

    在對(duì)地下結(jié)構(gòu)進(jìn)行有限元?jiǎng)恿Ψ治龅倪^(guò)程中,需要把半無(wú)限的土體轉(zhuǎn)換為有邊界的空間,而波動(dòng)在傳播的過(guò)程中遇到邊界會(huì)產(chǎn)生反射,這與實(shí)際情況并不相符,所以在建立模型時(shí)必須對(duì)人工邊界進(jìn)行處理,使之符合實(shí)際情況[1-2]。當(dāng)前常用的人工邊界條件中,有透射邊界[3]、黏性邊界[4]、黏彈性邊界[5]等,其中透射邊界為位移型邊界,在多次透射的情況下精度較高,但是容易出現(xiàn)高頻震蕩[6],且不容易在有限元軟件中實(shí)現(xiàn);黏性邊界為應(yīng)力型人工邊界,但是其僅考慮了對(duì)散射波能量的吸收,并沒(méi)有考慮到邊界處介質(zhì)的彈性恢復(fù)能力,因此容易導(dǎo)致低頻失穩(wěn)問(wèn)題[7];而黏彈性邊界則克服了上述缺點(diǎn),能夠很好地模擬地基的彈性恢復(fù)力和輻射阻尼效應(yīng),具有較高的精度[8-9]。除此之外,杜修力等[10-11]結(jié)合黏彈性邊界條件,將地震作用的總波場(chǎng)進(jìn)行分解,提出通過(guò)荷載的方式在邊界處施加地震作用,取得了很好的效果。

    雖然人工邊界設(shè)置和地震波施加的理論取得了很多成果,但是如何將這些方法應(yīng)用到實(shí)際的有限元分析中去依然是目前需要探討的問(wèn)題[12]。其中孫奔博等[13]采用combine14單元通過(guò)耦合結(jié)點(diǎn)的方式在Anasys中實(shí)現(xiàn)了黏彈性邊界的施加,并對(duì)水工隧洞進(jìn)行了有限元分析;谷音等[2]、劉晶波等[14-15]開(kāi)發(fā)了等效黏彈性邊界單元方法,用等效實(shí)體單元來(lái)替代彈簧—阻尼器原件,在有限元軟件中實(shí)現(xiàn)了黏彈性邊界的施加;章小龍等[6]利用Fortran軟件對(duì)ABAQUS進(jìn)行了二次開(kāi)發(fā),實(shí)現(xiàn)了地震動(dòng)等效荷載的施加。但是對(duì)通過(guò)Python語(yǔ)言編程實(shí)現(xiàn)黏彈性邊界和地震動(dòng)等效荷載的統(tǒng)一施加的方法研究較少。

    Python是一種面向?qū)ο蟆⒔忉屝浴?dòng)態(tài)數(shù)據(jù)類型的高級(jí)程序設(shè)計(jì)語(yǔ)言[16-17],它的語(yǔ)法簡(jiǎn)潔清晰,易于理解和開(kāi)發(fā)。同時(shí)它還具有高效的數(shù)據(jù)結(jié)構(gòu),提供了一種簡(jiǎn)單但很有效的方式進(jìn)行面向?qū)ο蟮木幊蘙18],這使其成為大多數(shù)平臺(tái)上應(yīng)用于各個(gè)領(lǐng)域的理想的腳本語(yǔ)言和開(kāi)發(fā)環(huán)境,而ABAQUS軟件也為用戶預(yù)留了Python腳本的接口,以便于對(duì)軟件進(jìn)行二次開(kāi)發(fā),這為黏彈性邊界和地震動(dòng)輸入問(wèn)題的研究提供了新的思路。

    本文主要對(duì)ABAQUS的Python接口進(jìn)行了二次開(kāi)發(fā),并編寫了Python程序,在保證模擬精度的前提下實(shí)現(xiàn)了黏彈性邊界和地震波的快速施加,大大提高了邊界施加的效率,并利用該方法對(duì)淺埋隧洞的地震動(dòng)響應(yīng)進(jìn)行了分析,為進(jìn)一步研究地下工程的地震響應(yīng)提供了新的思路。

    1 黏彈性邊界和地震波施加原理

    1.1 黏彈性邊界

    在利用有限的空間模擬半無(wú)限土體時(shí),局部人工邊界條件的施加對(duì)于計(jì)算的準(zhǔn)確性有著很大的影響,許多學(xué)者也在這方面進(jìn)行了大量的研究,其中劉晶波等[12]研究的黏彈性人工邊界由于具有比較高的精度[8]和實(shí)用性,因而被廣泛應(yīng)用于土—結(jié)構(gòu)的有限元軟件計(jì)算中。其通過(guò)在離散的邊界網(wǎng)格結(jié)點(diǎn)上施加并聯(lián)的彈簧和阻尼器來(lái)實(shí)現(xiàn)從有限域到無(wú)限域的外行波吸收,在三維分析中該邊界的施加方式如圖1所示。其中彈簧的剛度和阻尼器的阻尼由下式計(jì)算:

    式中:[K1]、[K2]分別為切向(X、Y方向)需要設(shè)置彈簧的剛度;[C1]、[C2]分別為切向(X、Y方向)需要設(shè)置阻尼器的阻尼;[K3]為法向(Z向)需要設(shè)置的彈簧的剛度;[C3]為法向(Z向)需要設(shè)置的阻尼器的阻尼值;ρ為介質(zhì)的密度;R為地下結(jié)構(gòu)中心到邊界點(diǎn)的距離;[Ai]為每一個(gè)網(wǎng)格結(jié)點(diǎn)上的附屬面積之和;[cs]、[cp]、G分別為土體介質(zhì)的剪切波速,壓縮波速和剪切模量。

    1.2 黏彈性邊界下地震波的施加原理

    地震波在入射到有限土體中以后,遇到結(jié)構(gòu)或者土體表面會(huì)發(fā)生散射,所以地震的總波場(chǎng)一般可以分解為自由波場(chǎng)和散射波場(chǎng)[19]。這其中的散射波場(chǎng)由人工邊界吸收,而入射波則可以通過(guò)作用于人工邊界上的等效荷載的方式施加上去。除此之外,還需要根據(jù)不同邊界面上的波場(chǎng)特點(diǎn),采用不同的波場(chǎng)分解方案。在底部邊界面上可以將總場(chǎng)分解為邊界入射場(chǎng)和邊界外行場(chǎng),而在側(cè)面邊界處可以將總波場(chǎng)分解為自由場(chǎng)和散射場(chǎng)[12]。其中邊界自由場(chǎng)和入射場(chǎng)可以根據(jù)連續(xù)介質(zhì)力學(xué)模型通過(guò)解析計(jì)算得到,而邊界的外行場(chǎng)或者散射場(chǎng)則可以利用邊界條件進(jìn)行模擬。其中將散射場(chǎng)或者邊界外行場(chǎng)作為未知量,并將其用總場(chǎng)與自由場(chǎng)或者入射場(chǎng)的差值表示,就可以得到含外源作用的人工邊界面上任一點(diǎn)l在i方向的有限元運(yùn)動(dòng)方程[20]

    式中:Kli、Cli為邊界面上l結(jié)點(diǎn)的法向和切向的參數(shù),可由式(1)~(4)計(jì)算;fli為入射地震波作用下l結(jié)點(diǎn)在i方向的等效結(jié)點(diǎn)應(yīng)力,由式(6)表示

    式中,[KliuRli、CliuRli、σRli]分別為克服彈簧、阻尼器和土體介質(zhì)所需要克服的抗力。

    在模型的底面邊界垂直入射地震波,設(shè)P波和S波的位移時(shí)程為[up(t)]和[us(t)],利用波動(dòng)的傳播規(guī)律和波場(chǎng)的應(yīng)力狀態(tài),通過(guò)式(6)解析計(jì)算邊界面上每個(gè)點(diǎn)的等效地震荷載[fli(t)]。由此即可以得到任一結(jié)點(diǎn)的計(jì)算式[20]。

    2 ABAQUS-Python的二次開(kāi)發(fā)

    2.1 ABAQUS-Python的腳本接口

    ABAQUS中的腳本接口能夠非常高效地實(shí)現(xiàn)ABAQUS/CAE中的所有功能,二者的通信關(guān)系如圖2所示。當(dāng)Python程序編寫好以后,所有的命令都需要通過(guò)Python的特定解釋器進(jìn)行翻譯之后才能夠進(jìn)入到ABAQUS/CAE的內(nèi)核中執(zhí)行,在執(zhí)行的過(guò)程中會(huì)生成記錄文件,即擴(kuò)展名為.rpy的文件。這些命令在進(jìn)入內(nèi)核中以后將會(huì)被轉(zhuǎn)化為INP文件,最后程序?qū)NP文件輸入到求解器中進(jìn)行求解。

    2.2 二次開(kāi)發(fā)主要算法

    本次開(kāi)發(fā)使用的主要軟件為Python IDLE,首先需要導(dǎo)入ABAQUS相關(guān)函數(shù)模塊,程序會(huì)自動(dòng)讀取在ABAQUS中設(shè)置好的材料參數(shù),如彈性模量、泊松比、密度等;通過(guò)下式可以求得介質(zhì)的縱向和剪切波速

    式中:λ、μ為第一拉密常數(shù)和第二拉密常數(shù);G、ρ為介質(zhì)的剪切模量和密度。

    將式(7)、(8)以代碼的形式編寫到Python程序中去。獲取邊界面上所有的網(wǎng)格結(jié)點(diǎn),設(shè)成集合,并將總數(shù)記為h。通過(guò)式(1)~(4)分別計(jì)算該點(diǎn)需要施加的彈簧和阻尼元件的具體參數(shù),利用reader函數(shù)將地震波的位移和速度函數(shù)分別以表格的形式讀入到程序中來(lái),每次從集合中取出一個(gè)結(jié)點(diǎn),并利用for...in語(yǔ)句遍歷表格中所有數(shù)據(jù),以元組的形式表示出來(lái)。每次從邊界的集合中取出一個(gè)結(jié)點(diǎn),計(jì)算該點(diǎn)的相應(yīng)參數(shù),然后利用attend函數(shù),通過(guò)while語(yǔ)句創(chuàng)建列表函數(shù)u(t-Δti),并通過(guò)函數(shù)的運(yùn)算計(jì)算出該點(diǎn)對(duì)應(yīng)的荷載函數(shù),接著判斷該點(diǎn)處于人工截面的具體位置,將隨時(shí)間變化的荷載利用Python命令施加到該結(jié)點(diǎn)上去,同時(shí)將該點(diǎn)需要施加的邊界條件,即彈簧和阻尼器施加上去,依次循環(huán),直到每個(gè)節(jié)點(diǎn)都施加上相應(yīng)的荷載和邊界后,結(jié)束運(yùn)行。具體流程如圖3所示。

    2.3 程序核心代碼

    該程序的主要功能是在有限的土體介質(zhì)中施加地震波。只要給出土體介質(zhì)的彈性模量、密度、泊松比等基本參數(shù)以后,便可以實(shí)現(xiàn)黏彈性邊界和地震波的自動(dòng)施加,程序部分核心代碼如表1所列。

    將程序?qū)懞煤蟊4鏋?py格式的文件,當(dāng)建立好模型并劃分好網(wǎng)格以后,通過(guò)菜單欄的文件→運(yùn)行腳本→選擇相應(yīng)Python文件運(yùn)行即可。

    3 程序驗(yàn)證

    為了驗(yàn)證施加程序的有效性,在這里建立簡(jiǎn)單模型進(jìn)行驗(yàn)證。在ABAQUS中建立長(zhǎng)寬高均為50 m的立方土體,模型介質(zhì)參數(shù)如下:彈性模量E=92.1 MPa,泊松比為0.25,密度為1 880 kg/m3。在底部輸入加速度峰值為0.2g的天津波,地震波時(shí)長(zhǎng)為19.19 s,計(jì)算時(shí)間取為22 s,取介質(zhì)自由表面中心點(diǎn)的位移和加速度作為參考,時(shí)間增量取為0.01 s。為保證精度,網(wǎng)格大小為1 m × 1 m,土體外側(cè)邊界共有12 500個(gè)結(jié)點(diǎn),程序僅運(yùn)行80 s即將所有邊界條件和地震波元件施加完畢。模型計(jì)算結(jié)果如圖4、5所示。其中圖4為土體入射面處和土體表面中心點(diǎn)處的加速度曲線對(duì)比圖,圖5為二者的位移對(duì)比圖。

    由波動(dòng)的傳播原理可知,對(duì)于單一的彈性介質(zhì)而言,當(dāng)波動(dòng)傳播到介質(zhì)自由表面時(shí)其位移和加速度為入射時(shí)的2倍[3]。從圖中可以看出,土體自由表面處的位移放大了兩倍,同時(shí)加速度峰值也變?yōu)?.4g,同樣被放大了兩倍,由此可見(jiàn)利用此種方式施加的地震波具有良好的精度和效率,可以應(yīng)用于地震波的施加和地下結(jié)構(gòu)抗震的分析。

    4 地下淺埋隧洞的有限元分析

    為了考察本次二次開(kāi)發(fā)在地下結(jié)構(gòu)抗震中的效果,建立土體—隧洞結(jié)構(gòu)有限元模型,對(duì)地震動(dòng)在土體中的傳播規(guī)律和土體—結(jié)構(gòu)的相互作用進(jìn)行研究。在這里選用某工程實(shí)地土層,具體參數(shù)如表2所列,根據(jù)經(jīng)驗(yàn)公式計(jì)算,場(chǎng)地卓越周期為1.24 s。為了研究邊界和地震波的施加效果,土體深度取為60 m,橫向?qū)挾热?0 m,縱向長(zhǎng)度取為120 m,隧洞直徑9 m,埋深6 m(淺埋),具體模型如圖6所示。在模型底部輸入水平橫向地震動(dòng),加速度峰值分別為0.1g、0.2g、0.3g的天津波,其中0.1g天津波的加速度時(shí)程曲線和傅里葉譜如圖7所示。同時(shí)建立同等參數(shù)下的自由場(chǎng)模型,并在底部輸入水平橫向加速度為0.1g的天津波作為對(duì)比。

    4.1 土體豎直方向上的加速度變化和頻譜特性分析

    取0.1g天津波作用下自由場(chǎng)模型自由表面中心點(diǎn)的加速度,并繪出相應(yīng)的傅里葉譜如圖8所示。從圖中可以看出,相比于輸入時(shí)的地震波和傅氏譜(見(jiàn)圖7),地震波的加速度峰值由0.1g變?yōu)?.18g,波形變化不大,這說(shuō)明地震時(shí)土體本身對(duì)地震波的加速度有放大作用;從傅氏譜中可以看出,地表處接收到的地震波的低頻成分增多,高頻成分減少,傅里葉幅值的峰值頻率由1.22 Hz變?yōu)?.81 Hz。這一數(shù)值正好與土體的基本頻率吻合(f=1/T=1/1.24≈0.8 Hz),這是由于土體本身的基本頻率比較小,放大了地震波中的低頻成分,同時(shí)過(guò)濾了其中的高頻成分。計(jì)算結(jié)果符合地震動(dòng)在土體中的傳播規(guī)律[21]。

    當(dāng)土體中有地下結(jié)構(gòu)存在時(shí),地震波在土體中的傳播規(guī)律會(huì)發(fā)生改變,在這里選取0.1g天津波作用下的自由場(chǎng)模型和帶隧洞結(jié)構(gòu)模型,考察土體中心附近沿豎直方向的加速度放大系數(shù)變化,將土體自由表面處記為0 m,則基巖處為―60 m,每隔3 m取一個(gè)記錄點(diǎn),曲線如圖9所示。

    從圖中可以看出,帶隧洞土體和自由場(chǎng)土體的加速度放大系數(shù)基本一致,自底部至土體表面大致呈放大趨勢(shì),但是由于不同土層的性質(zhì)有所不同,所以放大程度并不一致。其中在淺層0~15 m埋深范圍內(nèi),由于隧洞結(jié)構(gòu)的存在使得圍巖中的加速度放大系數(shù)有所增加,同時(shí)也使得地表處的加速峰值變大。通過(guò)以上對(duì)比分析說(shuō)明利用該方法施加的邊界條件和地震波合理,地下淺埋隧洞結(jié)構(gòu)的存在對(duì)圍巖的加速度變化有一定影響。

    4.2 加速度在水平方向上的變化分析

    為了分析邊界條件施加對(duì)模型的影響程度,考察0.1g天津波輸入情況下自由場(chǎng)和帶結(jié)構(gòu)時(shí)土體自由表面中心處沿水平橫向(共90 m)和水平縱向(共120 m)的加速度峰值變化趨勢(shì),每隔6 m取一個(gè)記錄點(diǎn),曲線如圖10和圖11所示。

    對(duì)于土體表面水平橫向加速度峰值變化而言,從圖10中可以看出,在距離土體兩側(cè)邊界的20~30 m范圍內(nèi)加速度峰值波動(dòng)比較大,可見(jiàn)在這一范圍內(nèi)模型會(huì)受到邊界效應(yīng)的影響。而在中心附近的30 m范圍內(nèi)加速度峰值比較平穩(wěn),此時(shí)邊界效應(yīng)的影響減小。其中帶隧洞土體表面加速度最小值為0.208g,最大值為0.215g,僅比最小值增大3%;自由場(chǎng)表面加速度最小值為0.200g,最大值為0.206g,同樣僅比最小值增大3%,且數(shù)值波動(dòng)出現(xiàn)的位置在距離土體邊界2D~3D(D為結(jié)構(gòu)直徑)范圍內(nèi);這一規(guī)律在水平縱向曲線(見(jiàn)圖11)中同樣體現(xiàn)出來(lái),由此可以看出在靠近邊界附近加速度峰值有輕微波動(dòng),但是誤差較小,在可接受范圍之內(nèi),在中間段則比較平穩(wěn),效果較好。經(jīng)過(guò)多次試算結(jié)果表明,當(dāng)將模型尺寸取為9倍的結(jié)構(gòu)寬度時(shí),可以消除兩側(cè)的誤差帶來(lái)的影響。

    4.3 隧洞結(jié)構(gòu)在地震動(dòng)下的響應(yīng)分析

    隧洞結(jié)構(gòu)在地震動(dòng)輸入時(shí)會(huì)產(chǎn)生變形甚至破壞,為了研究該種邊界條件下地下結(jié)構(gòu)的動(dòng)力響應(yīng),取0.1g~0.3g天津波作用下隧洞結(jié)構(gòu)縱向長(zhǎng)度上頂部和底部的相對(duì)位移值和跨中的橫截面環(huán)向應(yīng)力值進(jìn)行分析。

    0.2g天津波作用下隧洞縱向長(zhǎng)度上頂部和底部的相對(duì)位移沿隧洞縱向的變化曲線如圖12所示,從圖中可以看出,隧洞兩端2D~3D(D為隧洞直徑)范圍內(nèi)的相對(duì)位移值較大,在中部趨于平穩(wěn),這是由于隧洞結(jié)構(gòu)在人為被截?cái)嘁院髢啥说募s束減弱,使得位移變大,這與實(shí)際的情況并不相符,在離開(kāi)端部一段距離后這種影響才逐漸減弱。所以在建模計(jì)算時(shí)需要將隧洞縱向端部長(zhǎng)度尺寸額外延長(zhǎng)2~3倍的結(jié)構(gòu)寬度,中間部分的穩(wěn)定值才可以認(rèn)為符合平截面假定。由于篇幅限制,本文僅給出0.1g~0.3g天津波作用下隧洞縱向跨中處的拱頂和拱底位移以及二者相對(duì)位移,如表3所列。從表中可以看出,隧洞的拱頂和拱底位移值稍大,但是相對(duì)位移比較小,這說(shuō)明圓形淺埋隧洞在地震動(dòng)作用下的地震響應(yīng)雖然比較大,但是整體性比較好,抗震性能可以保證[22],計(jì)算結(jié)果符合已有研究規(guī)律[23]。

    為了研究隧洞橫截面上的應(yīng)力分布規(guī)律,分別取0.1g~0.3g天津波作用下隧洞跨中處橫斷面,將拱頂處記為0°,沿著順時(shí)針?lè)较蛎?0°取一個(gè)觀測(cè)點(diǎn),即順時(shí)針?lè)较蚍謩e為0°(拱頂)、30°、60°、90°(拱底)、……330°、360°(與0°重合),橫斷面上各觀測(cè)點(diǎn)沿順時(shí)針?lè)较虻膽?yīng)力變化曲線如圖13所示。

    從圖中可以看出,3種情況下隧洞的應(yīng)力最大值均出現(xiàn)在120°~150°之間和210°~240°之間,而拱頂(0°和360°)和拱底(180°)處的應(yīng)力值非常小,說(shuō)明淺埋隧洞在V類圍巖下襯砌的頂部和底部的應(yīng)力值較小,斜下角處的應(yīng)力值較大,這與文獻(xiàn)[24]中的結(jié)論一致,分析結(jié)果合理,符合淺埋隧洞的地震響應(yīng)規(guī)律。

    5 結(jié)語(yǔ)

    本文對(duì)ABAQUS的Python腳本接口進(jìn)行了研究,成功開(kāi)發(fā)出基于黏彈性邊界的地震波統(tǒng)一施加程序,并對(duì)其精度進(jìn)行了驗(yàn)證,最后對(duì)分層土體淺埋隧洞的動(dòng)力響應(yīng)進(jìn)行了有限元分析,得到如下結(jié)論:

    1)該方法可以實(shí)現(xiàn)黏彈性邊界和地震波的統(tǒng)一施加,大大提高建模效率,同時(shí)可以有效消除波動(dòng)在人工界面上的反射,計(jì)算結(jié)果符合波動(dòng)在土體中的傳播規(guī)律,可以應(yīng)用于地下結(jié)構(gòu)抗震分析中。

    2)距離土體邊界附近一定距離內(nèi)的加速度峰值會(huì)出現(xiàn)波動(dòng),但是誤差范圍很小,當(dāng)將模型尺寸取為9倍的結(jié)構(gòu)寬度時(shí),可以消除兩側(cè)的誤差帶來(lái)的影響。

    3)隧洞結(jié)構(gòu)縱向端部2~3倍寬度度范圍內(nèi)的計(jì)算結(jié)果偏大,需要通過(guò)額外增加相應(yīng)的隧洞縱向長(zhǎng)度取值來(lái)解決這一問(wèn)題。

    4)該方法為地下結(jié)構(gòu)抗震分析提供了新的思路,也為進(jìn)一步研究地下結(jié)構(gòu)地震響應(yīng)規(guī)律奠定了基礎(chǔ)。

    參考文獻(xiàn):

    [1]? ? 解祎. 基于黏彈性人工邊界的地鐵隧道地震響應(yīng)分析[D]. 天津:天津大學(xué),2013.

    [2]? ? 谷音,劉晶波,杜義欣. 三維一致粘彈性人工邊界及等效粘彈性邊界單元[J]. 工程力學(xué),2007,24(12):31-37.

    [3]? ? 廖振鵬. 工程波動(dòng)理論導(dǎo)論 [M]. 2版. 北京:科學(xué)出版社,2002.

    [4]? ? LYSMER J,KULEMEYER R L. Finite dynamic model for infinite media[J]. Journal of Engineering Mechanics,ASCE,1969,95(4):859-877.

    [5]? ? DEEKS A J,RANDOLPH M F. Axisymmetric time-domain transmitting boundaries[J]. Journal of Engineering Mechanics,ASCE,1994,120(1):25-42.

    [6]? ? 章小龍,李小軍,陳國(guó)興,等. 黏彈性人工邊界等效荷載計(jì)算的改進(jìn)方法[J]. 力學(xué)學(xué)報(bào),2016,48(5):1126-1135.

    [7]? ? 劉晶波,呂彥東. 結(jié)構(gòu)-地基動(dòng)力相互作用問(wèn)題分析的一種直接方法[J]. 土木工程學(xué)報(bào),1998,31(3):55-64.

    [8]? ? 劉晶波,李彬. 三維黏彈性靜-動(dòng)力統(tǒng)一人工邊界[J]. 中國(guó)科學(xué)E輯:工程科學(xué) 材料科學(xué) 2005,35(9):966-980.

    [9]? ? LIU J B,DU Y X,DU X L,et al. 3D viscous-spring artificial boundary in time domain[J]. Earthquake Engineering and Engineering Vibration,2006,5(1):93-102.

    [10]? LIU J B,LIU S,DU X L. A method for the analysis of dynamic response of structure containing non-smooth contactable interfaces[J]. Acta Mechanica Sinica,1999,15(1):63-72.

    [11]? DU X L,WANG J T. Seismic response analysis of arch dam-water-rock foundation systems[J]. Earthquake Engineering and Engineering Vibration,2004,3(2):283-291.

    [12]? 邱流潮,金峰. 地震分析中人工邊界處理與地震動(dòng)輸入方法研究[J]. 巖土力學(xué),2006,27(9):1501-1504.

    [13]? 孫奔博,胡良明,付曉龍. 圍巖類型及襯砌厚度對(duì)水工隧洞地震動(dòng)力響應(yīng)的影響分析[J]. 水利水電技術(shù),2017,48(3):19-24,45.

    [14]? 劉晶波,谷音,杜義欣. 一致粘彈性人工邊界及粘彈性邊界單元[J]. 巖土工程學(xué)報(bào),2006,28(9):1070-1075.

    [15]? LIU J B,WANG V,LI B. Efficient procedure for seismic analysis of soil-structure interaction system[J]. Tsinghua Science and Technology,2006,11(6):625-631.

    [16]? 張強(qiáng),馬永,李四超. 基于Python的ABAQUS二次開(kāi)發(fā)方法與應(yīng)用[J]. 艦船電子工程,2011,31(2):131-134.

    [17]? 王輝,李廷春,王清標(biāo),等. 基于Python的ABAQUS二次開(kāi)發(fā)及其在淺埋偏壓隧道分析中的應(yīng)用[J]. 現(xiàn)代隧道技術(shù),2015,52(2):160-165.

    [18]? 李勇. Python web開(kāi)發(fā)學(xué)習(xí)實(shí)例[M]. 北京:清華大學(xué)出版社,2011.

    [19]? 黃勝,陳衛(wèi)忠,伍國(guó)軍,等. 地下工程抗震分析中地震動(dòng)輸入方法研究[J]. 巖石力學(xué)與工程學(xué)報(bào),2010,29(6):1254-1262.

    [20]? 杜修力,趙密. 基于黏彈性邊界的拱壩地震反應(yīng)分析方法[J]. 水利學(xué)報(bào),2006,37(9):1063-1069.

    [21]? 王國(guó)波,徐海清,于艷麗. “群洞效應(yīng)”對(duì)緊鄰交疊盾構(gòu)隧道及場(chǎng)地土地震響應(yīng)影響的初步分析[J]. 巖土工程學(xué)報(bào),2013,35(5):968-973.

    [22]? GB 50011-2010,建筑抗震設(shè)計(jì)規(guī)范[S]. 北京:中國(guó)建筑工業(yè)出版社,2011.

    [23]? 王國(guó)波,陳梁,徐海清,等. 緊鄰多孔交疊隧道抗震性能研究[J]. 巖土力學(xué),2012,33(8):2483-2490.

    [24]? 李鳳穩(wěn),孫常新. 水工隧洞地震作用計(jì)算模型的適用性研究[J]. 水利水電技術(shù),2017,48(4):42-46.

    [責(zé)任編輯 楊 屹]

    猜你喜歡
    二次開(kāi)發(fā)結(jié)點(diǎn)隧洞
    隧洞止水帶安裝質(zhì)量控制探討
    滇中引水工程大理段首條隧洞順利貫通
    淺談基于Revit平臺(tái)的二次開(kāi)發(fā)
    甘肅科技(2020年20期)2020-04-13 00:30:02
    淺談Mastercam后處理器的二次開(kāi)發(fā)
    模具制造(2019年3期)2019-06-06 02:11:02
    Ladyzhenskaya流體力學(xué)方程組的確定模與確定結(jié)點(diǎn)個(gè)數(shù)估計(jì)
    西門子Easy Screen對(duì)倒棱機(jī)床界面二次開(kāi)發(fā)
    自密實(shí)混凝土在水工隧洞襯砌中的應(yīng)用
    隧洞襯砌段滲水修復(fù)施工技術(shù)介紹
    基于Pro/E二次開(kāi)發(fā)的推土鏟參數(shù)化模塊開(kāi)發(fā)
    基于Raspberry PI為結(jié)點(diǎn)的天氣云測(cè)量網(wǎng)絡(luò)實(shí)現(xiàn)
    国产精品电影一区二区三区 | 久久精品人人爽人人爽视色| 男女之事视频高清在线观看| 首页视频小说图片口味搜索| 中文字幕av电影在线播放| 欧美精品亚洲一区二区| 亚洲国产毛片av蜜桃av| 午夜福利在线免费观看网站| 亚洲伊人久久精品综合| 窝窝影院91人妻| 久久久久久免费高清国产稀缺| 精品视频人人做人人爽| 久久中文字幕一级| av片东京热男人的天堂| 国产成人免费无遮挡视频| 亚洲欧美激情在线| 巨乳人妻的诱惑在线观看| 久久国产精品大桥未久av| 免费在线观看黄色视频的| 50天的宝宝边吃奶边哭怎么回事| 妹子高潮喷水视频| 日韩三级视频一区二区三区| 国产国语露脸激情在线看| 午夜免费鲁丝| 大型黄色视频在线免费观看| 性高湖久久久久久久久免费观看| 精品免费久久久久久久清纯 | 欧美亚洲 丝袜 人妻 在线| 午夜福利乱码中文字幕| 久久久久久人人人人人| 亚洲三区欧美一区| 999久久久精品免费观看国产| 久久久久久人人人人人| 国产日韩欧美亚洲二区| 免费在线观看影片大全网站| 精品亚洲成国产av| 日本欧美视频一区| 欧美变态另类bdsm刘玥| h视频一区二区三区| 黄片小视频在线播放| 91精品三级在线观看| 国产男女超爽视频在线观看| 国产成人精品无人区| 久久精品亚洲av国产电影网| 国产极品粉嫩免费观看在线| 一级片'在线观看视频| 国产精品九九99| 国产精品国产高清国产av | 亚洲精华国产精华精| 纵有疾风起免费观看全集完整版| 亚洲av第一区精品v没综合| 色精品久久人妻99蜜桃| 久久精品亚洲av国产电影网| 两人在一起打扑克的视频| 一级片'在线观看视频| 搡老乐熟女国产| 日日夜夜操网爽| 欧美人与性动交α欧美软件| 黑人巨大精品欧美一区二区mp4| 一本色道久久久久久精品综合| 最近最新中文字幕大全免费视频| 人妻久久中文字幕网| 欧美 亚洲 国产 日韩一| 天天躁夜夜躁狠狠躁躁| 91av网站免费观看| 日韩大片免费观看网站| 亚洲色图av天堂| 在线观看免费视频网站a站| 欧美激情高清一区二区三区| 国产成人欧美在线观看 | 亚洲国产成人一精品久久久| 女性生殖器流出的白浆| 首页视频小说图片口味搜索| 亚洲综合色网址| 91av网站免费观看| 精品午夜福利视频在线观看一区 | 丁香欧美五月| 18禁美女被吸乳视频| 最新在线观看一区二区三区| 9191精品国产免费久久| 亚洲第一av免费看| 嫁个100分男人电影在线观看| www.精华液| 黄色视频在线播放观看不卡| 在线观看免费视频网站a站| 欧美在线一区亚洲| 亚洲欧洲日产国产| 妹子高潮喷水视频| 精品国产一区二区三区久久久樱花| 桃花免费在线播放| 成人永久免费在线观看视频 | www.精华液| 搡老熟女国产l中国老女人| 热99国产精品久久久久久7| 久久久久久久国产电影| 亚洲国产看品久久| a在线观看视频网站| 精品第一国产精品| av线在线观看网站| cao死你这个sao货| 高潮久久久久久久久久久不卡| 黄色 视频免费看| 国产在线免费精品| 亚洲人成电影免费在线| 大片电影免费在线观看免费| 免费观看a级毛片全部| 咕卡用的链子| 脱女人内裤的视频| 国产99久久九九免费精品| av不卡在线播放| 黑人操中国人逼视频| 色老头精品视频在线观看| 亚洲一区中文字幕在线| 在线观看舔阴道视频| 日日摸夜夜添夜夜添小说| 国产精品熟女久久久久浪| 丰满饥渴人妻一区二区三| 乱人伦中国视频| 欧美变态另类bdsm刘玥| 亚洲中文av在线| 欧美精品人与动牲交sv欧美| 久久九九热精品免费| 国产精品秋霞免费鲁丝片| 黑人猛操日本美女一级片| 美女福利国产在线| 精品国产亚洲在线| 日韩 欧美 亚洲 中文字幕| 成人亚洲精品一区在线观看| 亚洲中文av在线| 中文字幕人妻熟女乱码| 亚洲专区字幕在线| 久久国产精品男人的天堂亚洲| 亚洲欧美一区二区三区黑人| 国产高清videossex| 国产精品98久久久久久宅男小说| 一区福利在线观看| 国产精品偷伦视频观看了| 91大片在线观看| 日韩中文字幕欧美一区二区| 久久精品人人爽人人爽视色| 狂野欧美激情性xxxx| 中文字幕精品免费在线观看视频| 久久亚洲精品不卡| 中文欧美无线码| 久久人妻av系列| 妹子高潮喷水视频| 日韩视频一区二区在线观看| 人妻久久中文字幕网| av免费在线观看网站| 蜜桃在线观看..| 久久九九热精品免费| 国产精品一区二区在线不卡| 婷婷成人精品国产| 99国产综合亚洲精品| 精品欧美一区二区三区在线| 建设人人有责人人尽责人人享有的| 老汉色av国产亚洲站长工具| 国产区一区二久久| 一个人免费看片子| 亚洲专区字幕在线| 欧美精品亚洲一区二区| 男人舔女人的私密视频| 国内毛片毛片毛片毛片毛片| 成人18禁在线播放| 亚洲成人免费av在线播放| 一级片免费观看大全| 欧美久久黑人一区二区| 亚洲成人免费av在线播放| 1024视频免费在线观看| 午夜91福利影院| 久久亚洲精品不卡| 色婷婷av一区二区三区视频| 又紧又爽又黄一区二区| 大陆偷拍与自拍| av网站在线播放免费| 超碰97精品在线观看| 老司机影院毛片| 免费久久久久久久精品成人欧美视频| 国产深夜福利视频在线观看| 最近最新免费中文字幕在线| 午夜91福利影院| 午夜福利一区二区在线看| 亚洲一卡2卡3卡4卡5卡精品中文| 水蜜桃什么品种好| 国产一区二区三区在线臀色熟女 | 免费女性裸体啪啪无遮挡网站| 国产人伦9x9x在线观看| 亚洲av成人不卡在线观看播放网| 久久久久久久久久久久大奶| 亚洲av片天天在线观看| 在线永久观看黄色视频| 操美女的视频在线观看| 热re99久久精品国产66热6| 91成年电影在线观看| 久久久国产一区二区| 亚洲精品国产区一区二| 久久免费观看电影| 无限看片的www在线观看| 一进一出好大好爽视频| 精品午夜福利视频在线观看一区 | av网站免费在线观看视频| 亚洲 欧美一区二区三区| 精品亚洲成国产av| 亚洲色图 男人天堂 中文字幕| 天天躁日日躁夜夜躁夜夜| 欧美变态另类bdsm刘玥| 十分钟在线观看高清视频www| 久久国产精品影院| 黄频高清免费视频| 精品高清国产在线一区| 欧美激情 高清一区二区三区| e午夜精品久久久久久久| 久热爱精品视频在线9| 欧美激情 高清一区二区三区| 国产一区二区在线观看av| 精品少妇一区二区三区视频日本电影| 久久人妻福利社区极品人妻图片| 国产高清videossex| 99国产极品粉嫩在线观看| 最近最新免费中文字幕在线| 欧美 亚洲 国产 日韩一| 十八禁网站免费在线| 成年女人毛片免费观看观看9 | 黑人巨大精品欧美一区二区蜜桃| 青草久久国产| 香蕉久久夜色| 91成人精品电影| 母亲3免费完整高清在线观看| 搡老乐熟女国产| 国产亚洲av高清不卡| 俄罗斯特黄特色一大片| 啦啦啦视频在线资源免费观看| 久久国产精品男人的天堂亚洲| 欧美日韩一级在线毛片| 午夜两性在线视频| 久久免费观看电影| 女人爽到高潮嗷嗷叫在线视频| 国产不卡av网站在线观看| 国产97色在线日韩免费| 欧美+亚洲+日韩+国产| 久久久久国内视频| 欧美日韩精品网址| 免费女性裸体啪啪无遮挡网站| xxxhd国产人妻xxx| 国产伦理片在线播放av一区| 免费av中文字幕在线| 国产男靠女视频免费网站| 丝袜人妻中文字幕| 九色亚洲精品在线播放| 51午夜福利影视在线观看| 丁香六月欧美| 久久亚洲精品不卡| 又大又爽又粗| 首页视频小说图片口味搜索| 大香蕉久久成人网| 日韩一卡2卡3卡4卡2021年| av一本久久久久| 两性午夜刺激爽爽歪歪视频在线观看 | 国产男女内射视频| 九色亚洲精品在线播放| 黄片小视频在线播放| 免费观看a级毛片全部| 国产免费视频播放在线视频| 嫩草影视91久久| 精品国产乱码久久久久久男人| 久久久国产一区二区| 国产亚洲av高清不卡| 男女床上黄色一级片免费看| bbb黄色大片| 欧美日韩国产mv在线观看视频| 宅男免费午夜| 久久久精品94久久精品| 在线观看人妻少妇| 菩萨蛮人人尽说江南好唐韦庄| 天堂俺去俺来也www色官网| 久久99一区二区三区| 精品高清国产在线一区| 亚洲,欧美精品.| 高清av免费在线| 女人被躁到高潮嗷嗷叫费观| videosex国产| 三级毛片av免费| 岛国在线观看网站| 国产又爽黄色视频| 国产精品av久久久久免费| 一区二区av电影网| 亚洲成人手机| 91成人精品电影| 一级片免费观看大全| 精品一区二区三区av网在线观看 | 国产一区二区三区综合在线观看| 人人妻人人澡人人看| 一本久久精品| 久久亚洲精品不卡| 最黄视频免费看| 亚洲专区字幕在线| 叶爱在线成人免费视频播放| netflix在线观看网站| 十八禁网站网址无遮挡| 美女国产高潮福利片在线看| 五月天丁香电影| 另类亚洲欧美激情| 纵有疾风起免费观看全集完整版| 国产男女内射视频| 亚洲av成人不卡在线观看播放网| 中文字幕人妻熟女乱码| 欧美 日韩 精品 国产| 日韩大码丰满熟妇| 又紧又爽又黄一区二区| 视频在线观看一区二区三区| 亚洲欧美精品综合一区二区三区| 大型av网站在线播放| 日韩免费高清中文字幕av| 国产一区二区激情短视频| 50天的宝宝边吃奶边哭怎么回事| 日本精品一区二区三区蜜桃| 久久人妻av系列| 日韩欧美一区二区三区在线观看 | 国产成人一区二区三区免费视频网站| 男人舔女人的私密视频| 人人妻人人澡人人看| 飞空精品影院首页| 欧美日韩一级在线毛片| 999久久久国产精品视频| 久久性视频一级片| 天堂俺去俺来也www色官网| 1024视频免费在线观看| 国产亚洲av高清不卡| 啪啪无遮挡十八禁网站| 精品福利观看| 电影成人av| 国产老妇伦熟女老妇高清| 99精国产麻豆久久婷婷| 午夜91福利影院| 国产免费现黄频在线看| 精品高清国产在线一区| 国产男女超爽视频在线观看| 香蕉久久夜色| 一夜夜www| 欧美乱码精品一区二区三区| 丝袜美腿诱惑在线| 99国产精品一区二区三区| 99精国产麻豆久久婷婷| 高清欧美精品videossex| 天天添夜夜摸| 又黄又粗又硬又大视频| 18禁美女被吸乳视频| 久久精品亚洲精品国产色婷小说| 少妇裸体淫交视频免费看高清 | 香蕉国产在线看| 日韩中文字幕欧美一区二区| 亚洲av国产av综合av卡| 日本黄色视频三级网站网址 | 精品第一国产精品| 日本精品一区二区三区蜜桃| 亚洲色图综合在线观看| 国产三级黄色录像| 一级毛片精品| 露出奶头的视频| 国产免费视频播放在线视频| 国产男女超爽视频在线观看| 国产成+人综合+亚洲专区| 国产精品av久久久久免费| 欧美日韩国产mv在线观看视频| 午夜福利在线免费观看网站| 精品久久久精品久久久| 日韩视频一区二区在线观看| 又黄又粗又硬又大视频| 69av精品久久久久久 | 如日韩欧美国产精品一区二区三区| 又紧又爽又黄一区二区| 久久九九热精品免费| 精品国产超薄肉色丝袜足j| 国产精品偷伦视频观看了| 久久久精品区二区三区| 久久香蕉激情| 亚洲伊人色综图| 国产一区二区在线观看av| 亚洲伊人色综图| 日本撒尿小便嘘嘘汇集6| 亚洲av成人不卡在线观看播放网| 国产区一区二久久| 国产又爽黄色视频| 免费av中文字幕在线| 在线永久观看黄色视频| 天堂8中文在线网| 人妻一区二区av| 欧美久久黑人一区二区| 亚洲精品国产精品久久久不卡| 国精品久久久久久国模美| 久久久久精品人妻al黑| 成人精品一区二区免费| 91av网站免费观看| 亚洲色图av天堂| 一边摸一边抽搐一进一小说 | 亚洲成人国产一区在线观看| 50天的宝宝边吃奶边哭怎么回事| 国产精品.久久久| 精品国产国语对白av| 香蕉丝袜av| 欧美 亚洲 国产 日韩一| 国产成人精品无人区| 午夜激情av网站| av超薄肉色丝袜交足视频| 黄色成人免费大全| 后天国语完整版免费观看| 丝袜在线中文字幕| bbb黄色大片| 亚洲少妇的诱惑av| 国产亚洲欧美精品永久| 宅男免费午夜| 中文字幕av电影在线播放| 国产精品一区二区在线观看99| 一区二区三区乱码不卡18| 欧美激情高清一区二区三区| 精品福利观看| 欧美成狂野欧美在线观看| 中文字幕制服av| 精品久久久久久电影网| 国产又爽黄色视频| 国产高清videossex| 99精品欧美一区二区三区四区| 高潮久久久久久久久久久不卡| 亚洲中文av在线| 咕卡用的链子| 欧美精品人与动牲交sv欧美| 丰满饥渴人妻一区二区三| 别揉我奶头~嗯~啊~动态视频| 91字幕亚洲| 99热网站在线观看| 美女午夜性视频免费| av一本久久久久| 日本精品一区二区三区蜜桃| 热re99久久国产66热| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品一卡2卡三卡4卡5卡| 桃红色精品国产亚洲av| 国产高清视频在线播放一区| 国产一区二区在线观看av| 国产日韩欧美在线精品| 国产高清国产精品国产三级| 亚洲精品国产精品久久久不卡| 午夜91福利影院| 国产福利在线免费观看视频| 日韩一区二区三区影片| 欧美人与性动交α欧美精品济南到| 嫩草影视91久久| 18禁裸乳无遮挡动漫免费视频| 在线观看免费视频网站a站| 精品第一国产精品| 国产精品一区二区在线观看99| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品免费一区二区三区在线 | 大片免费播放器 马上看| 狠狠婷婷综合久久久久久88av| 下体分泌物呈黄色| 制服诱惑二区| 99国产极品粉嫩在线观看| 色综合欧美亚洲国产小说| 久久毛片免费看一区二区三区| 黄片播放在线免费| 大陆偷拍与自拍| 下体分泌物呈黄色| 国产欧美日韩一区二区精品| 国产亚洲精品第一综合不卡| 色尼玛亚洲综合影院| 精品少妇黑人巨大在线播放| 在线观看舔阴道视频| 天天操日日干夜夜撸| 亚洲三区欧美一区| 精品久久久久久电影网| h视频一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 亚洲性夜色夜夜综合| 下体分泌物呈黄色| 久久久国产成人免费| 亚洲成人免费av在线播放| 欧美国产精品一级二级三级| 一二三四社区在线视频社区8| 男男h啪啪无遮挡| 人人妻人人爽人人添夜夜欢视频| 欧美成人免费av一区二区三区 | 9色porny在线观看| h视频一区二区三区| 黑人欧美特级aaaaaa片| 极品教师在线免费播放| 亚洲精品一卡2卡三卡4卡5卡| 男人操女人黄网站| 欧美精品高潮呻吟av久久| 丰满人妻熟妇乱又伦精品不卡| 亚洲熟妇熟女久久| 欧美国产精品va在线观看不卡| 久久天躁狠狠躁夜夜2o2o| 亚洲国产欧美日韩在线播放| 一本色道久久久久久精品综合| 国产精品国产av在线观看| 久久国产亚洲av麻豆专区| 亚洲国产精品一区二区三区在线| 另类精品久久| 在线av久久热| 久久av网站| 亚洲成人手机| 午夜福利在线观看吧| 久久久久久久国产电影| 2018国产大陆天天弄谢| 亚洲熟妇熟女久久| 女性被躁到高潮视频| 如日韩欧美国产精品一区二区三区| 十八禁高潮呻吟视频| 丰满人妻熟妇乱又伦精品不卡| 国产一区二区激情短视频| 国产老妇伦熟女老妇高清| 制服诱惑二区| 国产在线观看jvid| 欧美久久黑人一区二区| 香蕉丝袜av| 黄网站色视频无遮挡免费观看| aaaaa片日本免费| 久久毛片免费看一区二区三区| 午夜福利视频精品| 最近最新免费中文字幕在线| 韩国精品一区二区三区| 成人18禁在线播放| 亚洲av电影在线进入| 九色亚洲精品在线播放| 国产成人一区二区三区免费视频网站| 80岁老熟妇乱子伦牲交| 大片免费播放器 马上看| 中亚洲国语对白在线视频| 在线亚洲精品国产二区图片欧美| 熟女少妇亚洲综合色aaa.| 亚洲久久久国产精品| av不卡在线播放| 制服诱惑二区| 欧美在线黄色| 国产精品电影一区二区三区 | 国产亚洲欧美精品永久| 19禁男女啪啪无遮挡网站| 一边摸一边抽搐一进一出视频| 99香蕉大伊视频| 日韩 欧美 亚洲 中文字幕| 亚洲 欧美一区二区三区| 99国产精品一区二区蜜桃av | 国产午夜精品久久久久久| 欧美国产精品va在线观看不卡| 制服诱惑二区| 国产不卡一卡二| 啦啦啦中文免费视频观看日本| 天堂中文最新版在线下载| 蜜桃国产av成人99| 精品欧美一区二区三区在线| 欧美日本中文国产一区发布| 老司机亚洲免费影院| 美女主播在线视频| 在线观看免费视频网站a站| 男男h啪啪无遮挡| 亚洲av国产av综合av卡| 色在线成人网| xxxhd国产人妻xxx| 国产精品久久久久久人妻精品电影 | 久久久久久久精品吃奶| 纯流量卡能插随身wifi吗| √禁漫天堂资源中文www| 亚洲av电影在线进入| 亚洲中文字幕日韩| 一本—道久久a久久精品蜜桃钙片| 黄片小视频在线播放| 不卡av一区二区三区| 久久香蕉激情| 久久精品国产亚洲av香蕉五月 | 人人妻,人人澡人人爽秒播| 黄色丝袜av网址大全| 国产一区二区三区在线臀色熟女 | 五月天丁香电影| 精品国产乱码久久久久久小说| 久久久精品免费免费高清| 亚洲欧美一区二区三区久久| 国产精品偷伦视频观看了| 搡老岳熟女国产| 日日摸夜夜添夜夜添小说| 久久精品91无色码中文字幕| 免费看a级黄色片| 国产精品一区二区免费欧美| 黄频高清免费视频| 在线亚洲精品国产二区图片欧美| 色尼玛亚洲综合影院| 久久中文字幕一级| 免费在线观看影片大全网站| √禁漫天堂资源中文www| 国产av精品麻豆| 蜜桃国产av成人99| 欧美一级毛片孕妇| 亚洲av日韩精品久久久久久密| 精品少妇一区二区三区视频日本电影| 首页视频小说图片口味搜索| 久久久水蜜桃国产精品网| 亚洲专区国产一区二区| 日韩大片免费观看网站| 黄色视频不卡| 国产精品1区2区在线观看. | 亚洲欧洲日产国产| 99精国产麻豆久久婷婷| 老司机靠b影院| www.熟女人妻精品国产| 中文字幕人妻丝袜制服| 久久久久久久国产电影| 91大片在线观看| 亚洲一码二码三码区别大吗| 亚洲专区中文字幕在线| 亚洲欧洲精品一区二区精品久久久| 免费高清在线观看日韩| 亚洲精品国产区一区二| 成人18禁在线播放| 男女之事视频高清在线观看| 久久人人97超碰香蕉20202|