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

    天然橡膠支座大變形壓剪性能的雙非線性超彈性理論和實(shí)驗(yàn)研究

    2022-08-01 00:58:22蘇何先何穎成管慶松張巖巖
    工程力學(xué) 2022年8期
    關(guān)鍵詞:剪應(yīng)變疊層軸壓

    楊 靜,潘 文,蘇何先,何穎成,管慶松,張巖巖

    (1. 昆明理工大學(xué)建筑工程學(xué)院,云南,昆明 650500;2. 云南省抗震工程技術(shù)研究中心,云南,昆明 650500;3. 震安科技股份有限公司,云南,昆明 650000;4. 云南省設(shè)計(jì)院第五建筑設(shè)計(jì)有限公司,云南,昆明 650000)

    國(guó)內(nèi)外多次地震證實(shí),隔震結(jié)構(gòu)的減震效果非常明顯。同時(shí),考慮到未來(lái)需要對(duì)隔震建筑按更具體、更高的抗震設(shè)防目標(biāo)進(jìn)行性能化設(shè)計(jì),隔震橡膠支座作為隔震結(jié)構(gòu)的主要受力構(gòu)件,有必要對(duì)其進(jìn)行精細(xì)化分析研究。

    隔震支座是由薄鋼板和橡膠交替粘結(jié)在一起的復(fù)合構(gòu)件,它作為結(jié)構(gòu)的局部構(gòu)件,當(dāng)結(jié)構(gòu)受到外部激勵(lì)(靜力荷載、外部動(dòng)力或地面加速度)時(shí),構(gòu)件的反應(yīng)量有位移、速度、加速度,對(duì)于地震激勵(lì),這些量的絕對(duì)值和相對(duì)值可能都是需要的;但與結(jié)構(gòu)變形有關(guān)的相對(duì)位移最重要,只要支座每個(gè)受力狀態(tài)的位移確定,其內(nèi)力及應(yīng)力就可由每一瞬時(shí)的靜力分析確定[1]。因此,國(guó)內(nèi)外學(xué)者基于靜力狀態(tài)對(duì)隔震橡膠支座進(jìn)行了大量研究。

    Haringx[2]是關(guān)于疊層橡膠支座性能的經(jīng)典理論。該理論考慮了承受壓縮載荷和水平變形構(gòu)件的水平剛度和屈曲載荷。但是,由于Haringx 理論是在假設(shè)材料線彈性的情況下推導(dǎo)的,因此無(wú)法將其直接用于評(píng)估變形較大的疊層橡膠的非線性行為。而后,Kelly[3]同樣基于材料線彈性假定,建立了數(shù)學(xué)物理模型并推導(dǎo)了隔震橡膠支座的水平剛度理論計(jì)算公式,并分析了支座的屈曲穩(wěn)定。Chang 等[4]和Ding[5]分別采用解析剛度矩陣法和傳遞矩陣法,建立了離散的疊層橡膠支座的線性力學(xué)模型,該模型基于Haringx 理論,對(duì)具有不同幾何參數(shù)和材料性能的支座進(jìn)行逐層分析,計(jì)算內(nèi)力和位移。但是,在考慮彎矩剛度及材料屬性時(shí),仍假設(shè)橡膠物理屬性為線性 (E=3G)[6-7],未能反映大變形材料特性的非線性行為。Takaoka[8]建立了一種基于Haringx 理論的疊層橡膠支座非線性力學(xué)模型,但該模型預(yù)測(cè)疊層橡膠支座的大變形行為精度需取決于橡膠層的抗彎剛度和抗剪剛度非線性參數(shù)的合理設(shè)置。

    在此背景下,本文基于Haringx 理論提出了一種天然疊層橡膠支座大變形壓剪狀態(tài)的雙非線性方法。它的特征是可以通過(guò)將材料非線性引入到支座等效勻質(zhì)體的等效彈性模量Ec和彎曲剛度EcIs中,用EcA、EcIs表示支座變形的內(nèi)力,建立微分平衡方程來(lái)表達(dá)大變形區(qū)域中的雙非線性行為。同時(shí),結(jié)合實(shí)際工程中生產(chǎn)的足尺隔震支座的壓剪實(shí)驗(yàn),驗(yàn)證了該方法的嚴(yán)密性和準(zhǔn)確性。此外,隔震支座中橡膠屬于高分子聚合物,生產(chǎn)硫化過(guò)程中,已與硫化劑發(fā)生化學(xué)反應(yīng),具備了較高的彈性、耐熱性、延展性等[6-7],可以認(rèn)為隔震支座橡膠屬于理想超彈性材料,即不考慮變形對(duì)時(shí)間、溫度的依賴性。

    目前,非線性固體力學(xué)的發(fā)展已經(jīng)與材料微觀結(jié)構(gòu)的發(fā)展緊密結(jié)合[9]。本文研究對(duì)象為隔震橡膠支座,其中,橡膠在微觀結(jié)構(gòu)中存在明確的超彈性本構(gòu)方程[10]。另外,一般物體在外力作用下會(huì)發(fā)生形狀和尺寸的改變,剪切模量G和體積模量K正是與這兩種不同類(lèi)型的變形對(duì)應(yīng)的兩個(gè)基本模量[11],其中,G反映物體形狀的改變,K反映尺寸的變化。對(duì)于大多數(shù)材料,這兩個(gè)模量屬于同一量級(jí),但是,對(duì)于特殊材料,如疊層隔震支座中的橡膠,其K值的量級(jí)則遠(yuǎn)大于G,因此其壓縮變形將非常明顯地依賴于任何可能阻止形狀變化的邊界條件[11-12]。綜合考慮,有必要先從隔震疊層橡膠支座中隔離出單層橡膠進(jìn)行受力分析。

    1 單層橡膠單軸等效彈性模量和純彎剛度的推導(dǎo)

    橡膠屬于各向同性材料[10],因此,以下每個(gè)變形狀態(tài),橡膠內(nèi)部各點(diǎn)的剪切模量G均相同[9]。

    1.1 推導(dǎo)單軸等效彈性模量

    1.1.1 假設(shè)純壓運(yùn)動(dòng)量,引入位移函數(shù)

    在實(shí)際工程中,因隔震層平面內(nèi)剛度及支座豎向剛度都很大,豎向壓力作用下,豎向變形較均勻;且單層橡膠水平向受到上、下鋼板較大的約束。根據(jù)以上特點(diǎn),作如下2 條假設(shè):

    1) 橡膠各點(diǎn)變形過(guò)程中,水平面保持不變;

    2) 變形前的垂直線受豎向力后,各點(diǎn)變形后為拋物線,其中垂直線兩端與鋼板接觸各點(diǎn)無(wú)變形。物理坐標(biāo)系如圖1 所示。

    圖1 單層橡膠純壓變形圖Fig. 1 Single-layer rubber pure pressure deformation diagram

    設(shè)未發(fā)生變形前,橡膠各點(diǎn)三維坐標(biāo)為μ0(x,y),υ0(x,y),ω(z)。變形后,位移函數(shù)為:

    1.1.2 推導(dǎo)等效彈性模量

    根據(jù)橡膠支座變形過(guò)程中體積保持不變[7,13],略去線應(yīng)變更高階的微量,其體應(yīng)變?yōu)椋?/p>

    代入位移函數(shù)得:

    在厚度t范圍內(nèi),對(duì)變量z積分:

    取單層橡膠任一點(diǎn)為微元體,其三維受力如圖2 所示。

    圖2 微元體受力圖Fig. 2 Micro-body force diagram

    由幾何方程計(jì)算切應(yīng)變:

    圓形橡膠支座的變形僅沿極坐標(biāo)徑向方向:

    根據(jù)廣義胡克定律:

    代入x、y方向平衡微分方程:

    不考慮支座體力,得:

    上兩式各對(duì)x、y求導(dǎo),代入式(1),得:

    對(duì)于不可壓縮彈性固體,有內(nèi)部約束方程[10]:

    式中: σ為應(yīng)力張量;S為法向偏張量;p為受力方向主應(yīng)力。支座受力狀態(tài)處于非扭曲狀態(tài),Sx=0,Sy=0,有:

    式(3)為應(yīng)力 σz的調(diào)和函數(shù)[14],即拉普拉斯算子,它的重要特點(diǎn)是,滿足此算子的變量,即使固體受力狀態(tài)屬于大變形,當(dāng)與其他符合調(diào)和函數(shù)的變量相加減時(shí),仍符合線性疊加原理。

    在極坐標(biāo)下,拉普拉斯算子存在簡(jiǎn)單變換式[14],且存在邊界條件r=R時(shí) , σz=0,則:

    對(duì)式(4)兩邊在半徑R范圍內(nèi)積分,得:

    1.2 單層橡膠純彎等效彎曲剛度的推導(dǎo)

    由1.1 節(jié)結(jié)論所得Ec,就可進(jìn)一步推導(dǎo)壓剪支座與Ec相關(guān)的等效彎曲剛度。

    1.2.1 假設(shè)純彎運(yùn)動(dòng)量,修正位移函數(shù)

    單層橡膠層受彎矩作用時(shí),作如下假設(shè):

    1) 因鋼板彎曲剛度很大,假設(shè)鋼板保持平面轉(zhuǎn)動(dòng),設(shè)相鄰鋼板之間的轉(zhuǎn)角為 α;

    上式對(duì)x求導(dǎo),得:

    不考慮y方向扭轉(zhuǎn),則有:

    由邊界條件z=0 時(shí) ,ω=0, ?ω/?z=0,得:

    設(shè)未發(fā)生變形前,橡膠各點(diǎn)三維坐標(biāo)為μ0(x,y), υ0(x,y),ω(z);變形后,變形如圖3,且位移函數(shù)為:

    圖3 單層橡膠純彎變形圖Fig. 3 Single-layer rubber pure bending deformation diagram

    1.2.2 推導(dǎo)純彎等效彎曲剛度

    單層橡膠純彎調(diào)和函數(shù)的推導(dǎo)同1.1.2 節(jié)純軸壓,則有:

    因彎曲關(guān)于y軸對(duì)稱,并由邊界條件,得豎向應(yīng)力解析式:

    在極坐標(biāo)下,應(yīng)力對(duì)yz面求矩,并在整個(gè)截面上積分,求總截面彎矩:

    因截面轉(zhuǎn)動(dòng) α較小,引入鋼板轉(zhuǎn)動(dòng)曲率半徑ρ,且截面直徑 2R≥t,則有: α ρ=t,1/ρ=α/t。

    由1.1 節(jié)推導(dǎo)結(jié)論Ec,有:

    2 壓剪隔震橡膠支座數(shù)學(xué)模型方程的建立與通用解答

    2.1 假設(shè)運(yùn)動(dòng)量,建立數(shù)學(xué)物理方程,并求通解

    非線性連續(xù)介質(zhì)力學(xué)研究物質(zhì)宏觀層次的運(yùn)動(dòng)、變形與破壞,與微觀層次相比它是一個(gè)更加復(fù)雜的系統(tǒng),它受制于微觀結(jié)構(gòu),又有自己獨(dú)特的規(guī)律,描述宏觀物體運(yùn)動(dòng)規(guī)律的方程一般是偏微分方程[9]。橡膠屬于易變形的固體,由其與鋼板交替粘結(jié)在一起的橡膠隔震支座,在豎向和水平荷載共同作用下會(huì)發(fā)生較大的水平位移和微小的彎曲變形,其運(yùn)動(dòng)形式簡(jiǎn)單明確,可以通過(guò)建立偏微分方程反映隔震橡膠支座宏觀的運(yùn)動(dòng)、變形。

    根據(jù)第1 節(jié)從橡膠微觀結(jié)構(gòu)推導(dǎo)的Ec與EcIs,可以假設(shè)疊層橡膠支座為符合Ec與EcIs的等效均質(zhì)體;同時(shí),考慮到支座中剪切變形基本上由橡膠產(chǎn)生,均質(zhì)體的幾何參數(shù)(As、Is)需考慮h/hr的放大系數(shù)[3](h為支座總高度、hr橡膠層總高度)。

    在壓剪作用下,疊層橡膠支座產(chǎn)生較大的水平剪切變形,故假設(shè)支座水平位移μ為第一個(gè)運(yùn)動(dòng)未知量。支座軸線還會(huì)產(chǎn)生彎曲變形,但因支座彎曲剛度較大,則橫截面轉(zhuǎn)動(dòng)變化也較小,可認(rèn)為橫截面轉(zhuǎn)動(dòng)后仍保持平面,但與變形后軸線不垂直,這是橡膠支座大變形的重要特征。為表示支座的彎曲變形程度,引入橫截面的轉(zhuǎn)角 β為第二個(gè)運(yùn)動(dòng)未知量。在幾何空間笛卡爾坐標(biāo)系中(圖4),設(shè)坐標(biāo)橫軸為u,沿支座高度方向?yàn)樽鴺?biāo)縱軸x。

    根據(jù)彈性力學(xué)微元體幾何方程的推導(dǎo)原理,把橡膠支座作為有限大固體[14],如圖5 所示,則最初相互垂直的線段OA與OB,運(yùn)動(dòng)到新平衡位置O′A′與O′B′時(shí),二者之間角度的改變量θ,就是切應(yīng)變?;谝陨显?,將圖4 中均質(zhì)體的下部分離體出來(lái),則在剖面處,新橫截面上剪力為V=GAs(μ˙-β) ; 同時(shí),轉(zhuǎn)角的導(dǎo)數(shù) β˙(即為軸線上該點(diǎn)的曲率),則新橫截面上彎矩為M=EcIsβ˙。由新截面上力的平衡及對(duì)中心點(diǎn)取矩,得平衡微分方程:

    圖4 均質(zhì)體支座數(shù)學(xué)物理簡(jiǎn)圖Fig. 4 Mathematics and physics diagram of homogeneous body support

    圖5 有限大支座幾何位移關(guān)系Fig. 5 Geometric displacement relations of finite bearings

    因轉(zhuǎn)角很小,則 cosβ →1, sinβ →β,以上超越方程可簡(jiǎn)化為:

    將M、V代入式(8)、式(9)得二元一階偏微分方程組:

    求解關(guān)于μ、β的二元一階偏微分平衡方程。

    隔震橡膠支座未知運(yùn)動(dòng)量μ、β的通用解答為:

    2.2 隔震橡膠支座運(yùn)動(dòng)量和內(nèi)力的解析解

    當(dāng)支座上、下兩端不限制轉(zhuǎn)動(dòng)時(shí),式(12)、式(13)有未知常數(shù)A、B,其通解不能直接應(yīng)用到隔震橡膠支座中。但因在實(shí)際工程中,橡膠隔震支座上、下截面分別與剛度較大的隔震層和下支墩相連接,相比之下,橡膠支座部分的剪切剛度較小,因此,支座數(shù)學(xué)物理模型簡(jiǎn)化時(shí),可假定支座上、下截面為固定端。基于此假定,橡膠支座頂部在水平剪力與軸力共同作用下,則其運(yùn)動(dòng)量方程式(12)、式(13)的邊界條件為:

    代入邊界條件,可求得方程中的未知常數(shù):

    由此可知,在壓剪共同作用下,橡膠隔震支座任一橫截面水平位移μ與 轉(zhuǎn)角 β有唯一解答(軸壓P產(chǎn) 生的P-Δ效應(yīng),不需要反復(fù)迭代,即可求精確解)。且支座頂部的初始彎矩M0不獨(dú)立,它是關(guān)于P、FH的函數(shù),從另一角度,在物理意義上,以上結(jié)論可解釋壓剪支座滿足變形協(xié)調(diào)方程。

    最后,將所求運(yùn)動(dòng)量μ、β代入式(8)、式(9),即得每個(gè)橫截面處對(duì)應(yīng)內(nèi)力M、V:

    3 隔震橡膠支座的理論解與實(shí)驗(yàn)結(jié)果對(duì)比分析

    3.1 實(shí)驗(yàn)隔震橡膠支座

    本壓剪實(shí)驗(yàn)選取實(shí)際隔震工程中的天然橡膠支座,其尺寸見(jiàn)表1。表2 則通過(guò)假設(shè)在等截面條件下,支座等效均質(zhì)體與其他常用材料(混凝土與鋼材)的材料屬性對(duì)比??芍刃Ь|(zhì)支座的剪切剛度遠(yuǎn)小于鋼材與混凝土,而彎曲剛度接近等截面混凝土。從另一個(gè)方面驗(yàn)證了第2 節(jié)中假定支座的剪切剛度很小、彎曲剛度很大是合理的。

    表1 橡膠支座尺寸Table 1 Size of rubber bearing

    表2 等截面不同材料屬性對(duì)比Table 2 Comparison of different material properties of equal section

    3.2 隔震橡膠支座壓剪實(shí)驗(yàn)與其 G的合理選取

    本次實(shí)驗(yàn)加載設(shè)備為電液伺服壓剪試驗(yàn)機(jī)(圖6),其最大軸向壓力15 000 kN,水平最大動(dòng)荷載1500 kN,豎向軸力按設(shè)計(jì)壓應(yīng)力12 MPa 保持不變,緩慢施加水平荷載,最大剪應(yīng)變400%(相當(dāng)于最大水平位移377 mm)[15-16]。為消除橡膠粘性、蠕變及松弛等[17-18]其他性質(zhì)的影響,取支座的第一個(gè)往復(fù)運(yùn)動(dòng)的1/4 曲線。其中加載曲線局部不光滑的原因,由于實(shí)驗(yàn)設(shè)備受到了摩擦的影響[16]。

    圖6 動(dòng)態(tài)壓剪試驗(yàn)機(jī)Fig. 6 Dynamic compression shear testing machine

    在非線性連續(xù)介質(zhì)力學(xué)中,當(dāng)固體受到較大的變形時(shí),變形和未變形構(gòu)型的位置關(guān)系用變形梯度張量表示,而材料屬性與空間矢量關(guān)系通過(guò)柯西變形張量描述[10]。與初始狀態(tài)相比,支座在不同剪應(yīng)變時(shí),位置與材料屬性都發(fā)生較大改變。而由疊層橡膠支座實(shí)驗(yàn)所得水平力與豎向、水平位移推導(dǎo)的G,綜合反映了支座在外力作用下,支座每個(gè)狀態(tài)幾何位置與材料屬性的雙非線性的改變。與1 節(jié)平衡微分方程中G反映了支座整體每個(gè)壓剪狀態(tài)的宏觀G相一致。

    3.3 隔震橡膠支座 FH-γ 與 G-γ理論、實(shí)驗(yàn)曲線對(duì)比分析

    依據(jù)圖7 所示框圖,可由實(shí)驗(yàn)曲線上每個(gè)壓剪狀態(tài) γ 對(duì)應(yīng)的G,迭代計(jì)算出該狀態(tài)橡膠支座的理論G與 γ值,同理,得出所有狀態(tài)的外力與變形FH-γ 理論曲線(圖8、圖9)。對(duì)比可知,F(xiàn)H-γ的實(shí)驗(yàn)與理論曲線幾乎100%重合,說(shuō)明考慮橡膠為非線性超彈材料時(shí),壓剪隔震橡膠支座的理論解為精確解答。

    圖7 剪切模量G 與剪應(yīng)變 γ計(jì)算程序框圖Fig. 7 Block diagram for calculation program of shear modulus G and shear strainγ

    圖8 隔震橡膠支座壓剪 G-γ實(shí)驗(yàn)、理論曲線(P=12 MPa)Fig. 8 Compression-shear G -γ experiment and theoretical curve of vibration isolation rubber bearing (P=12 MPa)

    圖9 隔震橡膠支座 FH-γ理論、實(shí)驗(yàn)曲線(P=12 MPa)Fig. 9 Compression-shear FH-γ experiment and theoretical curve of vibration isolation rubber bearing (P=12 MPa)

    綜上所述,考慮到橡膠支座的整體剪壓變形和材料屬性的共同影響,《建筑隔震橡膠支座》GB/T 20688.3-2006[15]在確定橡膠支座多遇地震( γ=100% ) 、罕遇地震( γ=250% ) 的 γ 對(duì)應(yīng)G時(shí),采用足尺橡膠支座進(jìn)行實(shí)驗(yàn)是合理的。

    4 隔震橡膠支座內(nèi)力分布規(guī)律和剪應(yīng)變、壓力相關(guān)性分析

    4.1 橡膠支座各內(nèi)力與運(yùn)動(dòng)量分布圖

    基于前面系統(tǒng)的闡述推導(dǎo),可對(duì)隔震支座某一個(gè)壓剪狀態(tài)進(jìn)行受力分析,如隔震支座在最大水平推力FH=700 kN、P=12 MPa 作用下,可迭代計(jì)算出該狀態(tài) γ=400%、G=0.89 MPa,進(jìn)而可得各截面運(yùn)動(dòng)量及內(nèi)力值(圖10)。由圖可知,橡膠支座最大水平位移出現(xiàn)在頂端,最大彎矩在支座兩端,最大剪力與轉(zhuǎn)角在支座中部。橫截面轉(zhuǎn)動(dòng)角度均小于 1 ( 1=0.017 rad),此處驗(yàn)證了第2 節(jié)中假定勻質(zhì)體的轉(zhuǎn)動(dòng)角度 β較小是符合實(shí)際情況的。這些結(jié)論對(duì)于判斷支座產(chǎn)生破環(huán)的具體位置有一定的參考價(jià)值[19]。

    圖10 沿截面高度內(nèi)力 M、V 及運(yùn)動(dòng)量 μ 、β的分布圖Fig. 10 Distribution diagram of internal force M、V and movement amount μ 、β along the section height

    4.2 橡膠支座剪應(yīng)變相關(guān)性分析

    4.1 節(jié)中確定了橡膠支座最大內(nèi)力所在的危險(xiǎn)截面,將結(jié)論代入式(4)、式(6)及式(7),可進(jìn)一步確定此截面內(nèi)各點(diǎn)的豎向應(yīng)力。不同剪應(yīng)變下,支座頂部沿水平推力軸線豎向應(yīng)力分布圖11(P=12 MPa)。由圖可知,隨剪應(yīng)變的增加,橡膠支座頂端截面壓應(yīng)力較大區(qū)域的豎向壓應(yīng)力逐漸增大,壓應(yīng)力小的區(qū)域壓應(yīng)力減小逐漸變?yōu)槔瓚?yīng)力,且拉應(yīng)力隨剪應(yīng)變?cè)黾?,絕對(duì)值也增加。但支座中心點(diǎn)的壓應(yīng)力為定值,大小為2P/A[20]。

    圖11 軸壓不變、不同剪應(yīng)變時(shí)橡膠支座豎向應(yīng)力分布圖Fig. 11 The vertical stress distribution diagram of the rubber bearing when the axial pressure is constant and the shear strain is different

    4.3 軸壓變化引起效應(yīng)P-Δ 與內(nèi)力變化規(guī)律

    為反映軸力變化對(duì)橡膠支座頂部水平位移及內(nèi)力影響,仍選取4.1 節(jié)中支座的這一受力狀態(tài)(剪應(yīng)變400%、軸壓應(yīng)力P=12 MPa、G=0.89 MPa)。因剪應(yīng)變不變,且實(shí)驗(yàn)軸壓與設(shè)計(jì)壓應(yīng)力差別不大時(shí),橡膠G受壓力的影響較小,則不同軸力下,可不考慮支座G的變化[15]。將外力及G代入式(12)~式(14),可得支座內(nèi)力與位移理論解,且隨軸力的變化規(guī)律如圖12、表3 所示。

    圖12 不同軸壓,支座頂端水平位移與內(nèi)力變化圖Fig. 12 The horizontal displacement of the top of the support and the internal force change diagram under different axial pressure

    表3 不同軸力下,支座頂部截面內(nèi)力與水平位移理論值Table 3 The theoretical values of the internal force and horizontal displacement of the top section of the support,under different axial forces

    由以上圖表可知,隨橡膠支座軸力增大,水平位移呈非線性增加,但增加比例不大,即P-Δ效應(yīng)不明顯,但是,支座端彎矩增加較大,因此橡膠內(nèi)應(yīng)力隨軸壓的變化不可忽視[21-24]。

    4.4 壓應(yīng)力相關(guān)性分析

    從圖13~圖15 可知,當(dāng)剪應(yīng)變?cè)酱髸r(shí),軸壓改變引起豎向應(yīng)力的變化越明顯,即壓力相關(guān)性強(qiáng);而當(dāng)剪應(yīng)變小時(shí),軸壓變化引起豎向應(yīng)力的變化不明顯,即壓力相關(guān)性弱。

    圖13 剪應(yīng)變100%、不同軸壓支座頂部豎向應(yīng)力變化圖Fig. 13 The vertical stress change diagram on the top of the support, with 100% shear strain and different axial compression

    圖14 剪應(yīng)變400%、不同軸壓支座頂部豎向應(yīng)力變化圖Fig. 14 The vertical stress change diagram on the top of the support, with 400% shear strain and different axial compression

    圖15 剪應(yīng)變100%、400%,不同軸壓下支座頂部豎應(yīng)力對(duì)比圖Fig. 15 Comparison of the vertical stress on the top of the support, when the shear strain is 100%, 400%,and different axial compressions

    5 純剪隔震橡膠支座數(shù)學(xué)模型方程的建立與通用解答

    當(dāng)豎向軸力P=0,支座僅受水平力FH作用時(shí),由于原微分方程中沒(méi)有了運(yùn)動(dòng)量的一次項(xiàng),微分方程需要重新求解。

    將P=0 代入式(10)、式(11),得:

    求解方程并由邊界條件得未知運(yùn)動(dòng)量及截面內(nèi)力解答:

    由式(18)可知,x=h時(shí),支座頂部水平位移最大,即:

    式(20)中第一項(xiàng)為彎曲引起的位移μE,對(duì)于長(zhǎng)度遠(yuǎn)大于橫截面尺寸的構(gòu)件,其水平位移可以僅考慮第一項(xiàng)彎曲的影響,如工程中不考慮軸力影響的層間柱抗彎剛度;第二項(xiàng)為剪力貢獻(xiàn)的位移μG,對(duì)于高度與橫截面尺寸相近,而材料參數(shù)GAs?EcIs的構(gòu)件,第二項(xiàng)提供了較大位移,如此文中的隔震橡膠支座,在水平力FH=700 kN 作用下,各項(xiàng)結(jié)果見(jiàn)表4,且可知式(20)與式(10)、式(11)在P=0 極限解答(見(jiàn)表3)相同。

    表4 純剪切內(nèi)力與位移理論值Table 4 Theoretical values of internal force and displacement in pure shear

    6 結(jié)論

    本文采用將橡膠材料的微觀結(jié)構(gòu)和疊層橡膠支座的宏觀運(yùn)動(dòng)、變形相結(jié)合的方法,從數(shù)學(xué)的角度,結(jié)合物理運(yùn)動(dòng)學(xué)與靜力學(xué),先后以單層橡膠和支座整體為研究對(duì)象,最后用關(guān)鍵性實(shí)驗(yàn)來(lái)驗(yàn)證理論結(jié)果,合理解釋了在水平、豎向荷載耦合作用下橡膠支座的幾何、材料雙非線性的大變形特征。主要得到以下結(jié)論:

    (1) 根據(jù)單層橡膠推導(dǎo)結(jié)果的GAs與EcIs,建立了橡膠支座在兩種外部荷載同時(shí)作用下,反映水平變形與彎曲程度的偏微分平衡方程,得到了橡膠支座大剪切變形下,豎向與水平兩種荷載耦合作用下大剪切變形的幾何非線性通用解答。

    (2) 因橡膠是各向同性材料,每個(gè)剪應(yīng)變狀態(tài),宏觀上支座有相同的剪切模量,但不同剪應(yīng)變有不同的剪切模量。因此,確定各個(gè)狀態(tài)疊層橡膠支座的剪切模量非常重要。在第(1)部分理論解答基礎(chǔ)上,開(kāi)展了實(shí)際工程中橡膠支座的壓剪實(shí)驗(yàn),依據(jù)實(shí)驗(yàn)G-γ 曲線,得到支座FH-γ的理論曲線,與FH-γ實(shí)驗(yàn)曲線幾乎完全重疊,解決了超彈性橡膠支座大剪切變形的材料非線性問(wèn)題。

    (3) 根據(jù)以上理論,對(duì)隔震支座某壓剪狀態(tài)進(jìn)行受力分析,得到支座各截面運(yùn)動(dòng)量與內(nèi)力解析式,可知,橡膠支座最大水平位移出現(xiàn)在頂端,最大彎矩在支座兩端,最大剪力與轉(zhuǎn)角在支座中部。這些結(jié)論對(duì)于判斷支座產(chǎn)生破環(huán)的具體位置有很好的參考價(jià)值。

    (4) 剪應(yīng)變不變,考慮軸力變化,得到以下規(guī)律,隨橡膠支座軸力增大,水平位移呈非線性增加,但增加比例不大,即P-Δ效應(yīng)不明顯;而支座端彎矩增加較大,因而橡膠內(nèi)應(yīng)力隨軸壓的變化不可忽視。

    (5) 為了驗(yàn)證不同剪應(yīng)變時(shí),軸壓相關(guān)性的強(qiáng)弱,分別在剪應(yīng)變100%、400%,施加變化的軸力,對(duì)比表明,當(dāng)剪應(yīng)變?cè)酱髸r(shí),軸壓改變引起豎向應(yīng)力的變化越明顯,即壓力相關(guān)性強(qiáng);而當(dāng)剪應(yīng)變小時(shí),軸壓變化引起豎向應(yīng)力的變化不明顯,即壓力相關(guān)性弱。

    猜你喜歡
    剪應(yīng)變疊層軸壓
    復(fù)合材料修復(fù)含裂紋圓管的軸壓承載特性研究
    改良隔震橡膠支座的試驗(yàn)及模擬分析*
    難加工材料(CFRP/Ti)疊層自適應(yīng)制孔研究
    碳纖維三向織物約束混凝土軸壓力學(xué)性能研究
    疊層橡膠隔震支座技術(shù)在工程施工中的應(yīng)用
    水泥改良黃土路基動(dòng)力穩(wěn)定性評(píng)價(jià)參數(shù)試驗(yàn)研究
    鄢家橋水庫(kù)土壩剪應(yīng)變及穩(wěn)定分析
    利用GPS分析新疆地區(qū)形變特征
    科技視界(2015年27期)2015-10-08 10:43:12
    軸壓下鋼筋混凝土橋墩塑性鉸特性研究
    合康疊層母排進(jìn)軍軍工領(lǐng)域
    中文资源天堂在线| 成人美女网站在线观看视频| 亚洲欧美成人综合另类久久久| 熟女av电影| 香蕉精品网在线| 在线 av 中文字幕| 国产一区亚洲一区在线观看| 黑丝袜美女国产一区| 老师上课跳d突然被开到最大视频| 91精品伊人久久大香线蕉| 成人国产麻豆网| 小蜜桃在线观看免费完整版高清| 免费黄色在线免费观看| 国产视频首页在线观看| 少妇人妻久久综合中文| 男女免费视频国产| 联通29元200g的流量卡| a 毛片基地| 天堂8中文在线网| 又大又黄又爽视频免费| 最近最新中文字幕大全电影3| 午夜福利影视在线免费观看| 日本vs欧美在线观看视频 | 亚洲精品日本国产第一区| 伦理电影免费视频| 国产男女内射视频| 亚洲精品日韩在线中文字幕| 少妇高潮的动态图| 国产 一区 欧美 日韩| 国产免费视频播放在线视频| 能在线免费看毛片的网站| 精品久久久久久久久亚洲| 九九在线视频观看精品| 亚洲无线观看免费| 99热这里只有是精品50| 国产人妻一区二区三区在| 精品少妇久久久久久888优播| 亚洲国产精品专区欧美| 免费高清在线观看视频在线观看| 国产伦理片在线播放av一区| 免费黄频网站在线观看国产| 国产视频内射| 日日摸夜夜添夜夜爱| 97超视频在线观看视频| 免费高清在线观看视频在线观看| 亚洲欧美精品自产自拍| h视频一区二区三区| 99热这里只有精品一区| 久久99热这里只频精品6学生| 国产成人91sexporn| 日韩人妻高清精品专区| 亚洲无线观看免费| 老师上课跳d突然被开到最大视频| 简卡轻食公司| 人人妻人人添人人爽欧美一区卜 | 国产精品福利在线免费观看| 高清欧美精品videossex| 18禁在线播放成人免费| 街头女战士在线观看网站| 三级经典国产精品| 小蜜桃在线观看免费完整版高清| 免费人成在线观看视频色| 亚洲国产成人一精品久久久| 久久久久精品性色| 热99国产精品久久久久久7| 一区二区三区乱码不卡18| 校园人妻丝袜中文字幕| 亚洲欧美成人综合另类久久久| 成人毛片a级毛片在线播放| 在线精品无人区一区二区三 | 男女下面进入的视频免费午夜| 成人毛片a级毛片在线播放| 国产在线男女| 97在线人人人人妻| 97在线人人人人妻| 亚洲精品456在线播放app| a级毛色黄片| 建设人人有责人人尽责人人享有的 | 午夜免费男女啪啪视频观看| 亚洲美女搞黄在线观看| 黄色视频在线播放观看不卡| 老司机影院成人| 性色avwww在线观看| 男女国产视频网站| 又大又黄又爽视频免费| 国产精品无大码| 蜜臀久久99精品久久宅男| 久久久久久久久久久免费av| 80岁老熟妇乱子伦牲交| 在线亚洲精品国产二区图片欧美 | 最近中文字幕2019免费版| 成人高潮视频无遮挡免费网站| 国产欧美日韩精品一区二区| 熟妇人妻不卡中文字幕| 韩国高清视频一区二区三区| 国产精品嫩草影院av在线观看| 免费观看性生交大片5| 一边亲一边摸免费视频| 99热网站在线观看| 亚洲国产毛片av蜜桃av| 国产免费视频播放在线视频| 精品久久国产蜜桃| 少妇被粗大猛烈的视频| 日日啪夜夜撸| 亚洲电影在线观看av| 岛国毛片在线播放| 在线亚洲精品国产二区图片欧美 | 精品一区在线观看国产| 国产成人一区二区在线| 蜜桃久久精品国产亚洲av| av女优亚洲男人天堂| 精品久久国产蜜桃| 国产精品秋霞免费鲁丝片| 久久久久久久久久人人人人人人| 亚洲图色成人| 国产精品一区www在线观看| 色视频在线一区二区三区| 亚洲av男天堂| 国产黄频视频在线观看| 看免费成人av毛片| 午夜视频国产福利| 尤物成人国产欧美一区二区三区| 一级黄片播放器| 大香蕉97超碰在线| 91在线精品国自产拍蜜月| 日韩av免费高清视频| 欧美亚洲 丝袜 人妻 在线| 成人亚洲欧美一区二区av| 亚洲av中文字字幕乱码综合| 欧美成人午夜免费资源| 人妻少妇偷人精品九色| 久久99热这里只频精品6学生| 亚洲av福利一区| 国产精品偷伦视频观看了| 国产高清不卡午夜福利| 80岁老熟妇乱子伦牲交| 午夜精品国产一区二区电影| 久久久久精品性色| 毛片一级片免费看久久久久| 欧美少妇被猛烈插入视频| 久久久久精品久久久久真实原创| 国产男女超爽视频在线观看| 极品少妇高潮喷水抽搐| 国内精品宾馆在线| 丰满乱子伦码专区| 国产精品成人在线| 国产一区有黄有色的免费视频| 欧美精品国产亚洲| 国产成人a区在线观看| 国产成人freesex在线| 国产成人aa在线观看| 多毛熟女@视频| 天天躁夜夜躁狠狠久久av| 日韩国内少妇激情av| 国产一区有黄有色的免费视频| 久久久久网色| 天堂8中文在线网| 日韩制服骚丝袜av| 天天躁夜夜躁狠狠久久av| 亚洲在久久综合| 欧美成人午夜免费资源| av女优亚洲男人天堂| 免费大片18禁| 国产午夜精品久久久久久一区二区三区| 伊人久久精品亚洲午夜| 男人狂女人下面高潮的视频| 久久99热这里只频精品6学生| 免费观看a级毛片全部| 人体艺术视频欧美日本| 久久久欧美国产精品| 丝瓜视频免费看黄片| 国产精品99久久99久久久不卡 | 日日啪夜夜爽| 狂野欧美激情性bbbbbb| 亚洲aⅴ乱码一区二区在线播放| 国产欧美另类精品又又久久亚洲欧美| 熟女av电影| 777米奇影视久久| 国产免费又黄又爽又色| 色吧在线观看| 大又大粗又爽又黄少妇毛片口| 亚洲伊人久久精品综合| 狠狠精品人妻久久久久久综合| 高清日韩中文字幕在线| 大又大粗又爽又黄少妇毛片口| 国产精品秋霞免费鲁丝片| 乱系列少妇在线播放| 国产精品99久久久久久久久| 国国产精品蜜臀av免费| 少妇人妻精品综合一区二区| 国产成人freesex在线| 国产黄片视频在线免费观看| 日本一二三区视频观看| 欧美高清性xxxxhd video| 熟女av电影| 国产高潮美女av| 婷婷色综合大香蕉| 精品视频人人做人人爽| 国产在线男女| 人妻夜夜爽99麻豆av| 国产成人免费观看mmmm| 亚洲一级一片aⅴ在线观看| 欧美性感艳星| 久久久精品94久久精品| 久久午夜福利片| 亚洲国产成人一精品久久久| 国产伦在线观看视频一区| 日韩人妻高清精品专区| 亚洲国产av新网站| 亚洲国产欧美在线一区| 一区二区三区四区激情视频| 热99国产精品久久久久久7| 国产视频内射| 久久精品国产鲁丝片午夜精品| 欧美变态另类bdsm刘玥| 熟女电影av网| 少妇人妻一区二区三区视频| 一级片'在线观看视频| 高清欧美精品videossex| 精品亚洲乱码少妇综合久久| 欧美激情极品国产一区二区三区 | 国产精品不卡视频一区二区| 99精国产麻豆久久婷婷| 中文字幕av成人在线电影| 一区二区av电影网| 少妇猛男粗大的猛烈进出视频| 国产亚洲一区二区精品| 国产淫语在线视频| 久久国产精品大桥未久av | 欧美激情国产日韩精品一区| 免费观看在线日韩| 观看av在线不卡| 免费av中文字幕在线| 中文字幕精品免费在线观看视频 | 国产成人精品一,二区| 欧美最新免费一区二区三区| 久久精品久久久久久噜噜老黄| 久久精品久久精品一区二区三区| 久久 成人 亚洲| 久久青草综合色| 欧美xxxx性猛交bbbb| 免费少妇av软件| 色综合色国产| 中文欧美无线码| 精品久久久久久久久av| 狠狠精品人妻久久久久久综合| 国产在线男女| 男女边摸边吃奶| 国产亚洲欧美精品永久| 青春草亚洲视频在线观看| 国产淫片久久久久久久久| 又爽又黄a免费视频| 精品酒店卫生间| 亚洲精品久久久久久婷婷小说| 精品久久久久久久久亚洲| 青春草亚洲视频在线观看| 日韩一区二区三区影片| 欧美变态另类bdsm刘玥| 日本免费在线观看一区| www.色视频.com| 最近手机中文字幕大全| 成年av动漫网址| 三级经典国产精品| 亚洲熟女精品中文字幕| 一级毛片久久久久久久久女| 热99国产精品久久久久久7| 免费人成在线观看视频色| 永久免费av网站大全| 狂野欧美激情性xxxx在线观看| 久久国产亚洲av麻豆专区| 香蕉精品网在线| 岛国毛片在线播放| 国产精品偷伦视频观看了| 深爱激情五月婷婷| 最近2019中文字幕mv第一页| 97在线人人人人妻| 老司机影院毛片| 亚洲av.av天堂| 一区二区三区乱码不卡18| 午夜福利影视在线免费观看| 国产高清三级在线| 欧美xxⅹ黑人| 亚洲国产欧美在线一区| 国产亚洲精品久久久com| 欧美成人精品欧美一级黄| 色5月婷婷丁香| 亚洲精品自拍成人| 久久久国产一区二区| 亚洲精品日本国产第一区| 国产 一区 欧美 日韩| 国产男女超爽视频在线观看| 色吧在线观看| 精品人妻熟女av久视频| 成年女人在线观看亚洲视频| 男女免费视频国产| 久久久久精品久久久久真实原创| 精品少妇久久久久久888优播| 欧美3d第一页| 国产一区二区三区综合在线观看 | 午夜福利在线观看免费完整高清在| 99热国产这里只有精品6| 高清不卡的av网站| 少妇 在线观看| 亚洲欧美成人精品一区二区| 天天躁夜夜躁狠狠久久av| 免费av中文字幕在线| 性色av一级| 久久亚洲国产成人精品v| 久久精品人妻少妇| 夫妻午夜视频| 成人无遮挡网站| 日本黄色片子视频| 久久毛片免费看一区二区三区| 亚洲精品aⅴ在线观看| 大片电影免费在线观看免费| 搡老乐熟女国产| 欧美成人精品欧美一级黄| 色综合色国产| 亚洲欧美成人综合另类久久久| 精品久久久久久久末码| 国产久久久一区二区三区| 国产熟女欧美一区二区| 亚州av有码| 国产亚洲一区二区精品| 国产探花极品一区二区| 国产一区二区三区av在线| 免费观看无遮挡的男女| 中文字幕免费在线视频6| 成人综合一区亚洲| 少妇的逼水好多| 熟女av电影| 国产精品人妻久久久久久| 欧美高清性xxxxhd video| 狠狠精品人妻久久久久久综合| 在线观看人妻少妇| 80岁老熟妇乱子伦牲交| 国产伦在线观看视频一区| 黑丝袜美女国产一区| 黑丝袜美女国产一区| 欧美精品一区二区免费开放| 亚洲欧洲国产日韩| 婷婷色综合大香蕉| 国产一区亚洲一区在线观看| 超碰av人人做人人爽久久| 成人无遮挡网站| 少妇人妻精品综合一区二区| 高清在线视频一区二区三区| 国产又色又爽无遮挡免| 欧美激情极品国产一区二区三区 | 一区在线观看完整版| 国产精品久久久久成人av| 亚洲精华国产精华液的使用体验| 成人国产麻豆网| 少妇 在线观看| 五月开心婷婷网| 少妇的逼好多水| 国产v大片淫在线免费观看| 亚洲色图av天堂| 在线观看三级黄色| 久久精品夜色国产| 亚洲av二区三区四区| 国产精品欧美亚洲77777| 干丝袜人妻中文字幕| 黄色日韩在线| 国国产精品蜜臀av免费| 亚洲精品国产av蜜桃| 午夜福利网站1000一区二区三区| 少妇裸体淫交视频免费看高清| 亚洲精品乱码久久久v下载方式| 日韩人妻高清精品专区| 大片免费播放器 马上看| 欧美变态另类bdsm刘玥| 亚洲怡红院男人天堂| 高清不卡的av网站| 日韩一区二区视频免费看| 观看av在线不卡| 中文乱码字字幕精品一区二区三区| 岛国毛片在线播放| 有码 亚洲区| 亚洲图色成人| 国产亚洲5aaaaa淫片| 精品久久久久久久末码| 中国美白少妇内射xxxbb| 亚洲精品乱久久久久久| 欧美丝袜亚洲另类| 熟女av电影| 高清日韩中文字幕在线| 日韩一区二区三区影片| 一个人看视频在线观看www免费| 亚洲综合色惰| 亚洲欧美日韩另类电影网站 | 99久久中文字幕三级久久日本| 香蕉精品网在线| 欧美精品一区二区大全| 新久久久久国产一级毛片| 网址你懂的国产日韩在线| 久久久国产一区二区| 亚洲人成网站在线播| 蜜臀久久99精品久久宅男| 高清欧美精品videossex| 伦理电影大哥的女人| 欧美国产精品一级二级三级 | 六月丁香七月| 成年av动漫网址| 女人久久www免费人成看片| 久久精品久久久久久噜噜老黄| 久久青草综合色| 欧美极品一区二区三区四区| 日韩国内少妇激情av| 日本色播在线视频| 免费大片18禁| 久久人人爽人人片av| 国产在线男女| 两个人的视频大全免费| 七月丁香在线播放| 丰满少妇做爰视频| av黄色大香蕉| 99视频精品全部免费 在线| 九九久久精品国产亚洲av麻豆| videossex国产| 极品少妇高潮喷水抽搐| 国产av精品麻豆| 日韩中文字幕视频在线看片 | 国产伦理片在线播放av一区| 91久久精品电影网| 亚州av有码| 女性被躁到高潮视频| 一级毛片aaaaaa免费看小| 91午夜精品亚洲一区二区三区| 男女国产视频网站| 亚洲精品成人av观看孕妇| 国产精品不卡视频一区二区| 午夜日本视频在线| 成人特级av手机在线观看| 女的被弄到高潮叫床怎么办| 午夜日本视频在线| 国产老妇伦熟女老妇高清| av视频免费观看在线观看| 我的女老师完整版在线观看| 一区二区三区四区激情视频| 老熟女久久久| 欧美精品一区二区免费开放| 丝袜脚勾引网站| 美女脱内裤让男人舔精品视频| 天美传媒精品一区二区| av专区在线播放| 最近最新中文字幕免费大全7| 国产爽快片一区二区三区| 精品人妻偷拍中文字幕| 国产免费视频播放在线视频| 中文资源天堂在线| 亚洲国产欧美在线一区| 国产乱人偷精品视频| 亚洲精品日韩在线中文字幕| 激情 狠狠 欧美| tube8黄色片| 欧美精品一区二区免费开放| 国产老妇伦熟女老妇高清| 亚洲av二区三区四区| 国产黄片视频在线免费观看| 国产在线男女| 久久人人爽人人爽人人片va| 又爽又黄a免费视频| 一区二区av电影网| 免费看av在线观看网站| 免费高清在线观看视频在线观看| 日韩亚洲欧美综合| 国产精品久久久久久久电影| 久久国产亚洲av麻豆专区| 伦精品一区二区三区| 欧美另类一区| av在线老鸭窝| 亚洲国产精品专区欧美| 午夜精品国产一区二区电影| 亚洲真实伦在线观看| 少妇人妻精品综合一区二区| 国产女主播在线喷水免费视频网站| 国产成人免费无遮挡视频| 日本黄色日本黄色录像| 国产亚洲最大av| 赤兔流量卡办理| 国产亚洲91精品色在线| 老师上课跳d突然被开到最大视频| 国产真实伦视频高清在线观看| 99热全是精品| 久久久久久九九精品二区国产| 中文乱码字字幕精品一区二区三区| 大陆偷拍与自拍| 国产精品爽爽va在线观看网站| 成人综合一区亚洲| 免费不卡的大黄色大毛片视频在线观看| 国产v大片淫在线免费观看| 这个男人来自地球电影免费观看 | 国产免费又黄又爽又色| 一本—道久久a久久精品蜜桃钙片| 亚洲图色成人| 久久久久精品久久久久真实原创| 极品教师在线视频| 91aial.com中文字幕在线观看| a级毛色黄片| 久久久午夜欧美精品| 免费在线观看成人毛片| 亚洲精品乱码久久久久久按摩| 久久久久网色| 永久免费av网站大全| 国产精品熟女久久久久浪| 全区人妻精品视频| 丰满少妇做爰视频| 大陆偷拍与自拍| 一本久久精品| av福利片在线观看| 国产成人免费无遮挡视频| 国产欧美亚洲国产| 啦啦啦视频在线资源免费观看| 国产毛片在线视频| 丰满人妻一区二区三区视频av| 99热6这里只有精品| 美女脱内裤让男人舔精品视频| 精品一区在线观看国产| 99久久人妻综合| 边亲边吃奶的免费视频| 日韩成人av中文字幕在线观看| 免费黄频网站在线观看国产| 日产精品乱码卡一卡2卡三| 天天躁日日操中文字幕| 国产色爽女视频免费观看| 蜜臀久久99精品久久宅男| 美女内射精品一级片tv| av专区在线播放| 99热国产这里只有精品6| 国产 精品1| 亚洲色图av天堂| 激情 狠狠 欧美| 老女人水多毛片| 亚洲真实伦在线观看| 欧美另类一区| 欧美变态另类bdsm刘玥| a级毛片免费高清观看在线播放| 午夜免费观看性视频| 欧美另类一区| 国产精品人妻久久久久久| 99九九线精品视频在线观看视频| 亚洲性久久影院| 嫩草影院新地址| 国产在线视频一区二区| 欧美zozozo另类| 欧美+日韩+精品| 18禁在线播放成人免费| 国产伦理片在线播放av一区| 97在线视频观看| 日日摸夜夜添夜夜添av毛片| 久久人人爽av亚洲精品天堂 | 97超视频在线观看视频| 观看免费一级毛片| 啦啦啦在线观看免费高清www| 国产高清三级在线| 麻豆国产97在线/欧美| 蜜桃在线观看..| 日本av手机在线免费观看| 成人亚洲精品一区在线观看 | 国产成人a区在线观看| 国产亚洲欧美精品永久| 国产伦精品一区二区三区四那| 搡老乐熟女国产| 亚洲欧美日韩东京热| 一二三四中文在线观看免费高清| 国产成人午夜福利电影在线观看| 丰满乱子伦码专区| 日韩,欧美,国产一区二区三区| 深爱激情五月婷婷| 黄色配什么色好看| 性色avwww在线观看| 人体艺术视频欧美日本| 国产老妇伦熟女老妇高清| 亚洲av国产av综合av卡| 超碰97精品在线观看| 成人国产麻豆网| 亚洲av成人精品一区久久| 成人特级av手机在线观看| 色5月婷婷丁香| 日日摸夜夜添夜夜爱| 日本免费在线观看一区| 色婷婷久久久亚洲欧美| 亚洲精品日韩在线中文字幕| 亚洲精品456在线播放app| 国产成人91sexporn| 中文字幕久久专区| 精品一区在线观看国产| 欧美激情国产日韩精品一区| 91久久精品国产一区二区成人| 亚洲色图av天堂| 日韩亚洲欧美综合| 国产午夜精品久久久久久一区二区三区| 18禁动态无遮挡网站| 日韩大片免费观看网站| 在线观看免费视频网站a站| 亚洲国产欧美人成| 日日摸夜夜添夜夜添av毛片| 国产av国产精品国产| 美女脱内裤让男人舔精品视频| 国产精品嫩草影院av在线观看| 午夜日本视频在线| 女的被弄到高潮叫床怎么办| 国产伦精品一区二区三区视频9| 高清毛片免费看| 我的女老师完整版在线观看| 亚洲av中文字字幕乱码综合| 久久久久久久精品精品| 大陆偷拍与自拍| 久久久久精品性色| 在线观看免费日韩欧美大片 | 五月伊人婷婷丁香| 中文在线观看免费www的网站| 久久人人爽人人爽人人片va| kizo精华|