亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        斜向三維隔震支座的非線性力學(xué)性能試驗(yàn)及數(shù)值模擬研究

        2020-11-14 06:40:44張莉莉劉文光何文福
        工程力學(xué) 2020年11期
        關(guān)鍵詞:鉛芯剪應(yīng)變理論值

        張莉莉,劉文光,許 浩,楊 杰,何文福

        (1. 上海大學(xué)土木工程系,上海 200444;2. 上海核工程研究設(shè)計(jì)院有限公司,上海 200233)

        全球每年發(fā)生約500 萬次地震,造成大量的人員傷亡和財(cái)產(chǎn)損失,而隔震技術(shù)是應(yīng)對(duì)地震災(zāi)害的有效手段[1]。隔震技術(shù)通過在結(jié)構(gòu)基礎(chǔ)部分設(shè)置隔震層,使上部結(jié)構(gòu)與基礎(chǔ)分離,從而限制地震能量向上部結(jié)構(gòu)傳遞[2]。目前,隔震技術(shù)已成功運(yùn)用于許多建筑,并在地震中表現(xiàn)出良好的水平隔震效果。但目前廣泛采用的隔震裝置只在水平向有隔震效果,對(duì)豎向地震分量無能為力。近年來研究發(fā)現(xiàn),在震中和近斷層附近產(chǎn)生過較強(qiáng)的豎向地震動(dòng),其中一些豎向地震分量幅值與水平分量之比超過2/3,有的甚至超過1[3]。因此,開展三維隔震技術(shù)的研究是很有必要的。

        實(shí)現(xiàn)三維隔震的難點(diǎn)在于豎向高承載能力和低隔震剛度需求的平衡,目前國(guó)內(nèi)外學(xué)者開發(fā)的裝置主要包括厚層橡膠支座、碟形彈簧、空氣彈簧、液壓油缸等系統(tǒng)。Kashiwazaki 等[4]研究了一種由橡膠支座和空氣彈簧構(gòu)成的三維隔震支座,并完成了模型振動(dòng)臺(tái)試驗(yàn)。Uriu 等[5]在核電站設(shè)施隔震設(shè)計(jì)中提出一種由空氣彈簧、疊層橡膠支座、黏滯阻尼器及自動(dòng)控制裝置組成的三維隔震地板系統(tǒng),實(shí)測(cè)結(jié)果表明該系統(tǒng)有較好的隔震效果。顏學(xué)淵等[6]研究了一種由碟形彈簧組成的三維隔震抗傾覆支座,并對(duì)一高層鋼框架模型結(jié)構(gòu)進(jìn)行了隔離三向地震動(dòng)激勵(lì)的振動(dòng)臺(tái)試驗(yàn)。趙亞敏等[7]開發(fā)了一種由鉛芯橡膠支座和碟形彈簧組成的三維隔震系統(tǒng),并進(jìn)行了振動(dòng)臺(tái)對(duì)比試驗(yàn)。陳兆濤等[8]研發(fā)了一種適用于大跨空間結(jié)構(gòu)的豎向變剛度三維隔震裝置,通過改變組合液壓缸不同階段參與工作腔室的種類和數(shù)量可實(shí)現(xiàn)豎向變剛度特性。

        許浩等[9]提出了一種新型斜向滑動(dòng)型三維隔震支座,通過將鉛芯橡膠支座傾斜放置實(shí)現(xiàn)豎向變形和隔震效果,新型三維隔震支座的力學(xué)性能由壓剪變形狀態(tài)下鉛芯橡膠支座的剪切剛度和壓縮剛度決定。鉛芯橡膠支座的力學(xué)性能常采用簡(jiǎn)化雙線性模型評(píng)價(jià),然而在壓應(yīng)力和剪應(yīng)變作用下鉛芯橡膠支座剛度會(huì)產(chǎn)生一定衰減。Nagarajaiah等[10]和Ryan 等[11]研究了彈簧支座的非線性特性和穩(wěn)定性,并采用了雙彈簧模型以考慮材料和幾何非線性。Weisman 等[12]、Sanchez 等[13]和Han 等[14]探索了橡膠支座的非線性剛度和穩(wěn)定性。吳彬等[15]研究了鉛芯橡膠支座的非線性動(dòng)態(tài)力學(xué)性能與其結(jié)構(gòu)構(gòu)造及外加結(jié)構(gòu)動(dòng)載之間的關(guān)系。孫新陽等[16]基于不同截面形狀及不同尺寸的橡膠支座提出了在剪切變形下的歸一化壓縮剛度公式。王建強(qiáng)等[17]采用低周反復(fù)循環(huán)加載試驗(yàn)方法對(duì)鉛芯橡膠支座進(jìn)行了剪切性能試驗(yàn),研究了支座壓應(yīng)力和剪應(yīng)變對(duì)支座水平力學(xué)性能的影響。何文福等[18]針對(duì)隔震支座力學(xué)性能與水平變形及豎向荷載的耦合特性,進(jìn)行了水平力學(xué)性能與剪應(yīng)變和壓應(yīng)力的相關(guān)性研究。朱玉華等[19]研究了厚層鉛芯橡膠支座水平及豎向力學(xué)性能隨壓應(yīng)力、剪應(yīng)變等的變化規(guī)律,并與現(xiàn)有的力學(xué)性能理論值進(jìn)行了比較。Warn 等[20]對(duì)鉛芯橡膠支座水平和豎向動(dòng)力響應(yīng)進(jìn)行了分析,結(jié)果表明水平位移增大時(shí)其豎向剛度會(huì)減小。

        鉛芯橡膠支座在壓剪變形狀態(tài)下的非線性力學(xué)性能對(duì)新型三維隔震支座的性能有較大影響。為此,本文針對(duì)新型斜向滑動(dòng)型三維隔震支座,建立其非線性豎向剛度的計(jì)算理論,并通過靜力試驗(yàn)及數(shù)值模擬探討其變化趨勢(shì)。

        1 三維隔震支座構(gòu)造及變形機(jī)理

        斜向滑動(dòng)摩擦三維隔震支座(three-dimensional oblique sliding friction seismic isolator, 3D OSFSI)主要由上部用于水平隔震的鉛芯橡膠支座,中部的滑動(dòng)連接件、滑動(dòng)摩擦塊以及下部用于豎向隔震的鉛芯橡膠支座組合構(gòu)成,該新型三維隔震支座構(gòu)造示意圖如圖1 所示。通過將鉛芯橡膠支座傾斜放置,將結(jié)構(gòu)的豎向變形轉(zhuǎn)化為下部支座的壓剪變形和摩擦滑塊的摩擦滑動(dòng),同時(shí)實(shí)現(xiàn)較大的豎向承載能力和較小的豎向隔震剛度,從而實(shí)現(xiàn)豎向隔震。

        圖1 三維隔震支座構(gòu)造示意圖Fig. 1 Structural sketch of 3D OSFSI

        由于上下連接角塊具有一定的傾斜角度,作用在該裝置上的豎向荷載P分解為軸向分量F1和切向分量F2,傾斜支座發(fā)生壓剪變形,摩擦滑塊與滑動(dòng)連接件發(fā)生摩擦滑動(dòng),如圖2 所示。整個(gè)裝置在豎向由傾斜支座提供豎向剛度,由鉛芯的屈服變形和滑動(dòng)摩擦耗能。在三維隔震支座的變形過程中,由于摩擦的存在,支座的豎向初始剛度和耗能能力將得到顯著的提高。

        圖2 豎向隔震支座組合受力示意圖Fig. 2 Force sketch of vertical isolators

        水平方向上,在滑動(dòng)連接件處設(shè)置水平限位擋板,由于限位擋板限制了傾斜支座的水平變形,豎向隔震支座組合的水平剛度較大。水平地震作用下,水平支座發(fā)生剪切變形,將地震作用經(jīng)隔離后傳遞給上部結(jié)構(gòu),使上部結(jié)構(gòu)水平向處于剛性平動(dòng)狀態(tài)。由此,可實(shí)現(xiàn)結(jié)構(gòu)的三維隔震。

        2 三維隔震支座力學(xué)性能分析

        斜向滑動(dòng)摩擦三維隔震支座的力學(xué)性能包括水平剪切力學(xué)性能和豎向壓縮力學(xué)性能。整個(gè)隔震系統(tǒng)的力學(xué)性能可以看作是上部水平隔震系統(tǒng)和下部豎向隔震系統(tǒng)串聯(lián)組合的力學(xué)性能。

        水平隔震系統(tǒng)是一個(gè)鉛芯橡膠支座,其相關(guān)力學(xué)性能的研究已較為成熟,其剪切剛度[21]與壓縮剛度[22]分別為:

        豎向隔震支座組合可視為2 個(gè)傾斜的鉛芯橡膠支座并聯(lián),并聯(lián)后組合的力學(xué)模型如圖3 所示,Kv,L為豎向屈服前加載剛度,Kv,u,L為豎向屈服前卸載剛度,Kvd,L為豎向屈服后加載剛度,Kvd,u,L為豎向屈服后卸載剛度。豎向屈服后加支座加載剛度與卸載剛度隨豎向位移的增加呈現(xiàn)非線性變化,且加載剛度與卸載剛度呈現(xiàn)不平行、非對(duì)稱特性,與鉛芯橡膠支座的傳統(tǒng)雙線性力學(xué)模型存在顯著差異。

        圖3 三維隔震支座豎向力學(xué)模型Fig. 3 Vertical mechanical model of 3D OSFSI

        圖4 試驗(yàn)?zāi)P椭ё鵉ig. 4 Model bearing for static test

        由式(9)和式(10)發(fā)現(xiàn),三維隔震支座豎向剛度與支座面壓、摩擦系數(shù)以及支座傾斜角度有關(guān),支座豎向剛度隨支座面壓的改變呈現(xiàn)非線性變化趨勢(shì)。

        3 支座力學(xué)性能試驗(yàn)研究

        3.1 試驗(yàn)?zāi)P图霸囼?yàn)工況

        為探究三維隔震支座的豎向非線性力學(xué)性能,組裝了小試支座模型進(jìn)行豎向壓縮試驗(yàn)。試驗(yàn)?zāi)P椭ё褂? 只LRB100 組裝,一只水平放置用作水平隔震,兩只傾斜放置用作豎向隔震,如圖4 所示,LRB100 力學(xué)參數(shù)見表1。連接角塊傾斜角度為12°和15°,摩擦滑塊的材料為聚四氟乙烯(PTFE)。為保證滑動(dòng)連接件僅發(fā)生豎向位移,使兩只傾斜支座協(xié)同變形,本次試驗(yàn)中使用鋼擋板作為限位裝置。

        試驗(yàn)加載采用一臺(tái)電液伺服加載試驗(yàn)機(jī),其最大壓力為2000 kN,作動(dòng)器行程300 mm,最大加載速度為50 mm/min,最大采樣頻率512 Hz,試驗(yàn)加載采用位移控制,具體加載工況如表2 所示。

        3.2 試驗(yàn)結(jié)果

        水平支座與傾斜支座的靜力試驗(yàn)結(jié)果如圖5和圖6 所示,在不同壓應(yīng)力和剪應(yīng)變下支座的豎向和水平向力學(xué)性能均保持穩(wěn)定狀態(tài)。水平工況下各曲線包絡(luò)面積較大,滯回曲線豐滿,具有較強(qiáng)的耗能能力。隨著剪應(yīng)變?cè)黾?,滯回曲線越來越飽滿,耗能面積隨之增加。支座豎向壓縮試驗(yàn)從初始空載狀態(tài)開始加載,循環(huán)加載3 次,豎向工況下各滯回曲線呈扁梭形,循環(huán)加載3 次后曲線斜率幾乎沒有變化,且加載過程中未出現(xiàn)破壞,卸載后支座基本恢復(fù)原狀,說明支座豎向力學(xué)性能保持穩(wěn)定。

        表1 LRB100 力學(xué)參數(shù)設(shè)計(jì)值Table 1 Mechanical parameters of LRB100

        表2 試驗(yàn)加載工況Table 2 Experimental cases

        圖5 水平LRB100 支座靜力試驗(yàn)結(jié)果Fig. 5 Test results of horizontal LRB100

        圖6 傾斜LRB100 支座靜力試驗(yàn)結(jié)果Fig. 6 Test results of oblique LRB100

        靜力試驗(yàn)得到的支座力學(xué)性能參數(shù)及其與理論值的對(duì)比如表3 所示,支座水平及豎向剛度試驗(yàn)參數(shù)值和理論值的誤差在15%以內(nèi)。

        表3 LRB100 支座力學(xué)性能參數(shù) /(kN/mm)Table 3 Mechanical parameters of LRB100

        三維隔震支座模型的靜力試驗(yàn)結(jié)果如圖7 所示,支座加載剛度與卸載剛度不同,滯回曲線呈現(xiàn)出典型的非對(duì)稱特性,不同幅值下滯回曲線形狀相近,力學(xué)性能保持穩(wěn)定。對(duì)比12°與15°三維支座試驗(yàn)曲線,可見相同剪切變形下,支座傾斜角度越大,豎向變形越大,支座耗能能力增強(qiáng)。

        由于不同豎向位移下三維支座的豎向剛度有所不同,故取對(duì)應(yīng)傾斜支座剪應(yīng)變?yōu)?00%±25%時(shí)三維隔震支座豎向變形,由靜力試驗(yàn)得到的支座力學(xué)性能參數(shù)及其與理論值的對(duì)如表4 所示,試驗(yàn)參數(shù)值和理論值較吻合,豎向加載剛度和卸載剛度的誤差均在5%以內(nèi)。

        線性理論與非線性理論得到的豎向剛度理論值與試驗(yàn)值的對(duì)比如圖8 和圖9 所示,支座豎向加載剛度與卸載剛度隨著豎向位移增加而衰減,呈現(xiàn)非線性變化。12°和15°三維隔震支座在100%剪應(yīng)變(豎向變形8 mm 和10 mm)下豎向卸載剛度衰減率分別為16.6%和14.2%,15°三維隔震支座125%剪應(yīng)變(豎向變形13 mm)時(shí)線性理論和非線性理論與試驗(yàn)的誤差分別為17.5%和1.47%。故考慮壓縮應(yīng)力對(duì)水平剛度的影響以及長(zhǎng)期荷載作用下剛度衰減對(duì)三維支座豎向剛度影響的非線性理論更加符合三維支座的力學(xué)性能。

        圖7 三維支座滯回曲線Fig. 7 Hysteresis curves of 3D isolation bearings

        表4 三維支座理論值與試驗(yàn)結(jié)果對(duì)比Table 4 Comparison of vertical stiffness of experimental and theoretical values for 3D OSFSI

        圖8 12°三維支座豎向剛度不同理論值與試驗(yàn)值對(duì)比Fig. 8 Comparison of vertical stiffness of different theoretical values and experimental values for 3D OSFSI of 12°

        圖9 15°三維支座豎向剛度不同理論值與試驗(yàn)值對(duì)比Fig. 9 Comparison of vertical stiffness of different theoretical values and experimental values for 3D OSFSI of 15°

        4 支座力學(xué)性能數(shù)值分析

        4.1 支座有限元模型

        采用通用有限元軟件ABAQUS 建立三維隔震支座有限元分析模型如圖10 所示,鉛芯是理想的彈塑性體,采用雙線性等效強(qiáng)化模型,彈性模量為16.46 GPa,泊松比為0.44。橡膠采用雜交單元C3D8H 模擬,選用Mooney-Rivlin 模型參數(shù)取值為:C10=0.15 MPa,C01=0.039 MPa,D1=0.002 MPa?1。

        圖10 三維隔震支座有限元分析模型Fig. 10 Model of 3D isolation bearing

        4.2 結(jié)果對(duì)比分析

        圖11 為水平支座在剪應(yīng)變?yōu)椤?00%和±200%工況下滯回曲線試驗(yàn)結(jié)果與數(shù)值模擬結(jié)果的對(duì)比,不同剪應(yīng)變下,水平支座模擬的滯回曲線與試驗(yàn)滯回曲線基本吻合。圖12 為12°與15°三維支座數(shù)值模擬和試驗(yàn)滯回曲線的對(duì)比結(jié)果,對(duì)比工況為傾斜支座剪應(yīng)變100%±25%的工況。三維支座模擬效果與試驗(yàn)結(jié)果相互吻合,滯回曲線呈現(xiàn)為飽滿的不平行的四邊形形狀。

        數(shù)值模擬所得三維支座力學(xué)性能參數(shù)與試驗(yàn)結(jié)果的對(duì)比見表5。12°支座在豎向位移6 mm、8 mm 和10 mm 下加載剛度的誤差分別為?1.29%、?7.26%和?6.27%;卸載剛度誤差分別為13.53%、10.25%和9.3%。15°支座在豎向位移8 mm、10 mm和13 mm 下加載剛度的誤差分別為12.72%、?4.06%和?2.48%;卸載剛度誤差分別為?4.37%、?6.35%和?5.9%。表明數(shù)值模擬具有較高精度,可有效地評(píng)估三維隔震支座的力學(xué)性能。建立10°與18°三維支座的有限元模型,對(duì)其動(dòng)力荷載下的力學(xué)性能進(jìn)行數(shù)值模擬,作為試驗(yàn)結(jié)果的補(bǔ)充。

        10°和18°三維支座力學(xué)性能模擬值與理論值對(duì)比如表6 所示,模擬值與理論值吻合較好,不同豎向變形下,三維支座的屈服后加載剛度和卸載剛度誤差在均在10%以內(nèi)。

        圖11 不同剪應(yīng)變下水平支座試驗(yàn)和模擬滯回曲線對(duì)比Fig. 11 Comparison of test and simulation hysteretic curves of horizontal LRB100 in different shear strains

        圖12 三維支座滯回曲線試驗(yàn)值和模擬值對(duì)比Fig. 12 Comparison of hysteretic curves of test and simulation for 3D OSFSI

        表5 12°和15°三維支座數(shù)值模擬力學(xué)性能與試驗(yàn)對(duì)比Table 5 Comparison of vertical stiffness of test and simulation for 3D OSFSI of 12° and 15°

        4.3 支座應(yīng)力分布

        不同層鋼板、橡膠的應(yīng)力分布大致相似,分別取三維支座中傾斜支座的頂層鋼板與頂層橡膠,其應(yīng)力分布如圖13 所示。12°與15°傾斜支座頂層鋼板的Mises 等效應(yīng)力均主要集中在鉛芯圓孔四周,其他部分應(yīng)力分布均勻,大小相近。橡膠整體應(yīng)力較小,鉛芯四周應(yīng)力稍大。

        4.4 參數(shù)分析

        斜向滑動(dòng)摩擦三維隔震裝置的豎向性能主要與傾斜角度、摩擦系數(shù)、剪應(yīng)變相關(guān),基于前文建立的三維支座的非線性理論、靜力試驗(yàn)及數(shù)值模擬結(jié)果,可得到豎向剛度的參數(shù)相關(guān)性能如圖14和圖15 所示。

        表6 10°和18°三維支座數(shù)值模擬力學(xué)性能與理論對(duì)比Table 6 Comparison of vertical stiffness of theory and simulation for 3D OSFSI of 10° and 18°

        圖13 三維支座應(yīng)力分布圖Fig. 13 Stress distribution of 3D isolation bearing

        圖14 100%剪應(yīng)變下豎向剛度相關(guān)性能Fig. 14 Correlation of vertical stiffness at 100% shear strain

        圖15 不同剪應(yīng)變下豎向剛度相關(guān)性能Fig. 15 Correlation of vertical stiffness at different shear strains

        圖14 給出了傾斜支座100%剪應(yīng)變下豎向剛度隨傾斜角度和摩擦系數(shù)的變化趨勢(shì)。隨著傾斜角度的增加,加載剛度與卸載剛度均會(huì)減小,傾斜角度越小,豎向剛度變化幅度越大。同一傾斜角度下,摩擦系數(shù)越大,豎向加載剛度越大,卸載剛度越小。隨著傾斜角度的增加,摩擦系數(shù)對(duì)豎向剛度的影響逐漸減小,當(dāng)傾斜角度達(dá)到20°時(shí),摩擦系數(shù)對(duì)其豎向剛度基本無影響。圖15 給出了傾斜支座不同剪應(yīng)變下豎向剛度隨傾斜角度的變化趨勢(shì)。同一傾斜角度下,剪應(yīng)變?cè)酱?,豎向加載剛度和卸載剛度越小,隨著傾斜角度的增加,剪應(yīng)變對(duì)支座豎向剛度的影響逐漸減小。對(duì)比圖15 線性理論與非線性理論可發(fā)現(xiàn),隨著傾斜角度的增加,壓縮應(yīng)力及剛度衰減等因素對(duì)三維支座豎向剛度的影響逐漸減小,傾斜角度大于20°時(shí),兩種理論結(jié)果趨于一致,偏差僅為0.8%,壓縮應(yīng)力及剛度衰減等因素對(duì)三維支座的豎向剛度的影響可忽略。

        5 結(jié)論

        本文針對(duì)新型斜向滑動(dòng)型三維隔震支座,建立了其非線性力學(xué)性能的計(jì)算理論,并通過靜力試驗(yàn)及數(shù)值模擬探討了其參數(shù)相關(guān)性,得到主要結(jié)論如下:

        (1) 基于鉛芯橡膠支座的剪切/壓縮剛度計(jì)算公式,建立了斜向滑動(dòng)型三維隔震支座的非線性力學(xué)性能計(jì)算理論,在不同的壓應(yīng)力和豎向位移下,三維隔震支座呈現(xiàn)出非線性的力學(xué)性能。

        (2) 進(jìn)行了三維隔震支座縮尺模型的靜力試驗(yàn),支座力學(xué)性能穩(wěn)定,在不同變形和壓力作用下滯回曲線呈現(xiàn)非線性變化,豎向剛度理論值和試驗(yàn)值誤差在5%以內(nèi)。采用ABAQUS 通用有限元軟件進(jìn)行了靜力試驗(yàn)的數(shù)值反演分析,試驗(yàn)值與模擬值的誤差總體在15%以內(nèi)。所提出的計(jì)算理論和數(shù)值模擬方法可以有效模擬三維支座的非線性力學(xué)性能。

        (3) 進(jìn)行了傾角為10°與18°的支座力學(xué)性能數(shù)值模擬,并基于數(shù)值結(jié)果討論了支座的力學(xué)性能相關(guān)性。隨著傾斜角度增大,支座剛度減小,同時(shí)其非線性衰減幅度也逐漸減小,在相同傾斜角度下,支座變形越大,其剛度的衰減越大。

        猜你喜歡
        鉛芯剪應(yīng)變理論值
        改良隔震橡膠支座的試驗(yàn)及模擬分析*
        擴(kuò)招百萬背景下各省區(qū)高職院校新增招生規(guī)模測(cè)度研究
        水泥改良黃土路基動(dòng)力穩(wěn)定性評(píng)價(jià)參數(shù)試驗(yàn)研究
        組合變形實(shí)驗(yàn)中主應(yīng)力方位角理論值的確定
        鄢家橋水庫土壩剪應(yīng)變及穩(wěn)定分析
        ASME規(guī)范與JB/T4730對(duì)接焊縫超聲檢測(cè)的靈敏度差異探討
        利用GPS分析新疆地區(qū)形變特征
        科技視界(2015年27期)2015-10-08 10:43:12
        為什么橡皮擦能擦掉鉛筆道
        保護(hù)煤柱寬度的理論值分析
        河南科技(2014年5期)2014-02-27 14:08:25
        鉛芯橡膠支座在簡(jiǎn)支梁橋減隔震技術(shù)中的應(yīng)用研究
        国产精品一区二区 尿失禁| 91久久国产香蕉视频| 琪琪色原网站在线观看| 麻豆久久久9性大片| 亚洲最新版无码AV| 又爽又猛又大又湿的视频| 中文字幕日韩有码在线| 欧美一性一乱一交一视频| 欧美亚洲综合另类| av黄片免费在线观看| 人妖一区二区三区视频| 久久久久久国产精品免费免费| 中文字幕一区二区三区久久网站| 国产成人精品蜜芽视频| 亚洲国产综合精品中久| 国产免费爽爽视频在线观看| 国产亚洲精品久久久久久久久动漫| 国产午夜精品久久久久九九| 国产成人亚洲精品91专区高清| 五月停停开心中文字幕| 亚洲成人av在线第一页| 亚洲精品久久久久中文字幕一福利| 亚洲国产美女精品久久| 日韩极品免费在线观看| 激情综合五月开心婷婷| 无遮挡又爽又刺激的视频| 亚洲 日韩 在线精品| 网红尤物泛滥白浆正在播放| 日本动漫瀑乳h动漫啪啪免费| 国产精品成人99一区无码| 无码伊人久久大蕉中文无码| 大香蕉国产av一区二区三区| 人妻丰满熟妇av无码区| 国产精品入口牛牛影视| 少妇性l交大片免费快色| 精品国产粉嫩内射白浆内射双马尾| 无套内射蜜桃小视频| 国产成人精品无码一区二区老年人| 91九色视频在线国产| 国产伦理一区二区| 99久久人妻无码精品系列蜜桃|