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

    地面式鋼筋混凝土圓形水池液固耦合模態(tài)分析★

    2017-03-02 10:06:19王明振
    山西建筑 2017年2期
    關(guān)鍵詞:液固儲液水池

    高 霖 王明振,2

    (1.重慶文理學(xué)院建筑工程學(xué)院,重慶 402160; 2.中國地震局工程力學(xué)研究所 中國地震局地震工程與工程振動重點實驗室,黑龍江 哈爾濱 150080)

    地面式鋼筋混凝土圓形水池液固耦合模態(tài)分析★

    高 霖1王明振1,2

    (1.重慶文理學(xué)院建筑工程學(xué)院,重慶 402160; 2.中國地震局工程力學(xué)研究所 中國地震局地震工程與工程振動重點實驗室,黑龍江 哈爾濱 150080)

    對地面式鋼筋混凝土圓形水池進行了考慮液固耦合效應(yīng)的模態(tài)分析,并使用ADINA有限元軟件,分析了儲液狀態(tài)、儲液高度、水池半徑以及液體表面重力波對模態(tài)的影響,對比了儲液晃動頻率的理論和有限元計算結(jié)果,推導(dǎo)得到了水池儲液晃動基本周期計算公式。

    水池,液固耦合,ADINA,模態(tài)分析,晃動周期

    0 引言

    給水系統(tǒng)是典型的城市生命線系統(tǒng),是城市健康運營的物質(zhì)基礎(chǔ)。生命線系統(tǒng)的功能以整個系統(tǒng)的運營狀態(tài)為體現(xiàn),其中每一個環(huán)節(jié)對整個系統(tǒng)的功能均有不可忽略的作用。給水系統(tǒng)功能有取水、儲水、凈水和輸配水。給水系統(tǒng)中,完成儲水和凈水功能角色的則是水池,水池是給水系統(tǒng)中的重要構(gòu)筑物,對整個系統(tǒng)的功能有重要影響。歷次地震震害表明:水池結(jié)構(gòu)在強烈地震作用下易發(fā)生不同程度的破壞。結(jié)構(gòu)的模態(tài)特性對其在地震作用下的響應(yīng)有重要影響。通過模態(tài)分析可以得到結(jié)構(gòu)的頻率、振型及模態(tài)參與質(zhì)量等動力特性,由動力特性根據(jù)反應(yīng)譜可確定結(jié)構(gòu)的地震影響系數(shù)。1957年,Housner教授對剛性儲液結(jié)構(gòu)的模態(tài)和液動壓力進行了分析[1],后續(xù)較多專家學(xué)者也對儲液結(jié)構(gòu)模態(tài)進行了分析[2-5],但大多數(shù)結(jié)果是針對剛性儲液結(jié)構(gòu)給出的。我國GB 50032—2003室外給水排水和燃?xì)鉄崃こ炭拐鹪O(shè)計規(guī)范中尚未給出儲液水池的模態(tài)計算方法,不便于儲液水池的地震響應(yīng)計算。同時,由于水池這類結(jié)構(gòu)在正常使用時包括固體結(jié)構(gòu)和液體兩相,內(nèi)部液體的存在會影響水池結(jié)構(gòu)的模態(tài)特性,因此在研究儲液結(jié)構(gòu)模態(tài)特征時必須考慮兩相間液固耦合效應(yīng)的影響。

    綜上所述,為了得到合理的儲液水池模態(tài)分布特征,本文使用ADINA有限元軟件和理論公式法對地面式鋼筋混凝土圓形水池進行考慮液固耦合效應(yīng)的模態(tài)分析,研究模態(tài)分布的規(guī)律,為計算儲液水池的地震響應(yīng)奠定基礎(chǔ)。

    1 液固耦合水池有限元模型

    水池結(jié)構(gòu)尺寸的確定主要參照相關(guān)規(guī)范、圖集和實際工程案例,所分析水池的結(jié)構(gòu)形式均為地面式、無頂蓋的圓形水池。將水池的結(jié)構(gòu)形式、名稱、容量、底厚、壁厚、儲液半徑、最大儲液高度及配筋等基本參數(shù)列于表1中。為分析不同儲液量對水池振動模態(tài)的影響,同時考慮水池的日常工作狀態(tài),選取從無水至70%儲液量之間共8個儲液工況(如表2所示)對水池A,B,C分別進行模態(tài)分析。為便于分析結(jié)果的表述,使用編號“C-70%”表示容量2 000 m3水池、儲水高度為2.45 m的工況,其他工況以此類推。ADINA中液體材料的密度為1 000 kg/m3,體積模量2.3×109Pa,阻尼比0.16%[6]。

    表1 圓形水池基本參數(shù)

    表2 儲水量與儲水高度對應(yīng)表

    使用ADINA Parasolid幾何建模方式建立水池模型。液體采用3D-Fuild勢流體單元,可模擬與結(jié)構(gòu)單元的耦合問題。池體結(jié)構(gòu)采用3D-Solid單元,為混凝土材料,使用Truss單元中的Rebar選項設(shè)置鋼筋[7]。模型網(wǎng)格劃分時兼顧計算精度和速度。根據(jù)研究對象的特點,選用行列式搜索法(Determinant-Search Method)進行模態(tài)分析,同時得到模態(tài)參與因子,用以表征某階模態(tài)對結(jié)構(gòu)振動的貢獻程度。

    2 模態(tài)計算結(jié)果分析

    進行液固耦合系統(tǒng)模態(tài)分析時,首先對模型施加重力荷載進行靜力分析,在模態(tài)分析時選擇靜力計算結(jié)果進行重啟動,進行模態(tài)和參與因子計算。3種不同尺寸水池、8種儲液狀態(tài)、是否考慮表面重力波影響2種情況,共計算48個工況。對比分析各工況的計算結(jié)果,分析不同因素對模態(tài)的影響程度和趨勢。

    2.1 儲液狀態(tài)對模態(tài)的影響

    計算無液和儲液狀態(tài)時水池的模態(tài),分析儲液狀態(tài)對模態(tài)的影響。無液時,水池的振型主要為池壁的振動,且振動頻率較高。當(dāng)水池內(nèi)儲有液體時,無論儲液量的多少,液固耦合系統(tǒng)包括剛體模態(tài)、儲液低階晃動模態(tài)和結(jié)構(gòu)高階振動模態(tài)三種。模態(tài)計算時第一階模態(tài)為剛體模態(tài),其頻率非常小,無限接近于零,主要表現(xiàn)為流體所有節(jié)點上的壓力值相同,固體結(jié)構(gòu)在均布內(nèi)壓下的變形,剛體模態(tài)的參與質(zhì)量百分比幾乎為零。以A-50%工況為例,圖1為其剛體模態(tài)。從第二階開始直至結(jié)構(gòu)高階振動模態(tài)之間全部為儲液晃動模態(tài),儲液晃動模態(tài)的頻率相對較小,主要表現(xiàn)為池內(nèi)液體的晃動,水池結(jié)構(gòu)雖在晃動壓力下變形,但相對液體晃動而言,池壁的振動響應(yīng)是微小的。圖2為A-50%工況的儲液晃動一階模態(tài)。結(jié)構(gòu)振動模態(tài)頻率較高,不易被激發(fā),主要表現(xiàn)為池壁的振動,池內(nèi)液體為脈沖模態(tài),結(jié)構(gòu)一階振動模態(tài)如圖3所示。當(dāng)儲液量不同時,出現(xiàn)結(jié)構(gòu)振動模態(tài)的階數(shù)略有差異。

    對比分析各工況模態(tài)計算結(jié)果可知:液體的存在改變了結(jié)構(gòu)的振動模態(tài),由單一的結(jié)構(gòu)振動轉(zhuǎn)變?yōu)榘▌傮w模態(tài)、儲液低階晃動模態(tài)和結(jié)構(gòu)高階振動模態(tài)三種。不同儲液量時,相同階數(shù)三種模態(tài)的振型幾乎完全相同;儲液時的結(jié)構(gòu)高階振型與無液時的結(jié)構(gòu)振型幾乎完全相同,且出現(xiàn)的順序基本相同。儲液水池的儲液晃動模態(tài)很容易被激發(fā),且其振型參與質(zhì)量較大,對結(jié)構(gòu)的振動響應(yīng)影響較大。由于結(jié)構(gòu)的旋轉(zhuǎn)對稱性,振動模態(tài)大多成對出現(xiàn),僅振型方向不同。

    2.2 儲液高度對模態(tài)的影響

    將A水池不同儲液高度下的剛體模態(tài)頻率、前三階儲液晃動頻率和結(jié)構(gòu)高階振動頻率列于表3中。

    表3 A水池不同儲液高度時若干階模態(tài)頻率

    Hz

    從表3中可以看出:儲液時的第一階模態(tài)均為剛體模態(tài),其頻率無限小接近于零;隨著儲液高度的增大,同一階次的儲液晃動頻率逐漸增大,且儲液高度越大,數(shù)值增大越緩慢;對比結(jié)構(gòu)振動頻率的變化趨勢可以分析液體對結(jié)構(gòu)振動的影響,同一階次結(jié)構(gòu)振動頻率隨著儲液高度的增大而不斷減小,說明儲液高度越大,液體對結(jié)構(gòu)振動的影響程度越大,使得液固耦合效應(yīng)越顯著。

    2.3 液體表面重力波對模態(tài)的影響

    表面重力波存在于液體和氣體兩種流體的分界面上,以表面波的形式存在,并以重力為恢復(fù)力[8]。前述計算水池動力特性時均考慮了液體表面重力波,即液面為自由液面,其約束為Free Surface。當(dāng)約束住液面的Z向位移,即約束了液面的自由晃動,也即不考慮液體表面重力波,通過ADINA計算可直接得到結(jié)構(gòu)振動模態(tài)。將水池A儲水高度為10%,30%,50%和70%的考慮和不考慮液體表面重力波時的結(jié)構(gòu)振動頻率列于表4中。

    表4 考慮和不考慮液體表面重力波時水池A的結(jié)構(gòu)振動頻率 Hz

    從表4中可以看出:儲液高度相同時,對于同一階結(jié)構(gòu)振動頻率,考慮液體表面重力波比不考慮時的略小,但相差甚微,說明表面重力波對結(jié)構(gòu)振動模態(tài)有一定的影響,但影響程度較小。同時發(fā)現(xiàn)考慮與不考慮表面重力波時的模態(tài)振型基本相同,這與文獻[9]的研究結(jié)果相同。因此,求解結(jié)構(gòu)模態(tài)時可以不考慮液體表面重力波的影響。由于約束了液面的Z向位移,不考慮液體表面重力波時得不到儲液晃動模態(tài)。

    2.4 水池半徑對模態(tài)的影響

    將儲液高度為70%的A,B和C水池前三階儲液晃動頻率和結(jié)構(gòu)振動頻率列于表5中。計算時均考慮液體表面重力波。

    對于儲液晃動模態(tài),從表5中可以看出,儲液晃動頻率數(shù)值隨水池半徑的增大而減小;對比圖2,圖4和圖5,水池半徑不同時同一階儲液晃動振型基本相同。對于結(jié)構(gòu)振動模態(tài),頻率數(shù)值也隨水池半徑的增大而減小,水池半徑越大,液固耦合效應(yīng)越顯著;水池半徑不同時,結(jié)構(gòu)振動振型被激發(fā)的順序不同。水池半徑對儲液晃動和結(jié)構(gòu)振動頻率均有影響,水池半徑對儲液晃動振型影響較小,對結(jié)構(gòu)振動振型影響較大。

    表5 儲液高度70%的水池A,B和C的模態(tài)頻率計算結(jié)果 Hz

    3 儲液晃動頻率有限元法和理論法計算結(jié)果對比

    GB 50032—2003室外給水排水和燃?xì)鉄崃こ炭拐鹪O(shè)計規(guī)范中尚未給出儲液水池晃動基本周期計算公式。文獻[3]中給出了基于如下假定的儲液結(jié)構(gòu)內(nèi)液體振動頻率計算公式,見式(1),認(rèn)為儲液結(jié)構(gòu)為絕對剛體,液體為無粘、無旋、不可壓縮的理想液體,且做小振幅運動。

    (1)

    由式(1)計算所分析水池的儲液晃動頻率,列于表6中,并與ADINA計算結(jié)果進行對比。

    表6 儲液晃動頻率理論與有限元計算結(jié)果對比 Hz

    從表6中可以看出:不同儲液高度、不同水池半徑下,儲液晃動前三階頻率的理論和ADINA計算結(jié)果相差很小,也說明了所建立的有限元模型合理可靠。由于有限元計算時儲液結(jié)構(gòu)具有一

    定的彈性模量,而并非剛性,液面為自由液面且考慮表面重力波的影響,使得兩種方法計算的頻率值有一定的差異?;谏鲜龇治?,推導(dǎo)得到圓形水池儲液晃動基本周期計算公式,見式(2)。儲液晃動基本周期是計算儲液結(jié)構(gòu)地震響應(yīng)的重要參數(shù)。

    (2)

    其中,T1為儲液晃動基本周期,s;R為儲液半徑,m;HW為儲液高度,m。

    4 結(jié)語

    本文對地面式鋼筋混凝土圓形水池進行了考慮液固耦合效應(yīng)的模態(tài)分析,主要得到以下5條結(jié)論:1)儲液結(jié)構(gòu)由于內(nèi)部液體的存在,改變了結(jié)構(gòu)的動力特性。儲液結(jié)構(gòu)模態(tài)包括剛體模態(tài)、以液體晃動為主的儲液低階晃動模態(tài)和以結(jié)構(gòu)振動為主的結(jié)構(gòu)高階振動模態(tài)三種。2)水池半徑相同時,儲液高度越大,同一階次儲液晃動頻率越大,同一階次結(jié)構(gòu)振動頻率越小。儲液高度越大,液固耦合效應(yīng)越顯著。3)液體表面重力波對儲液結(jié)構(gòu)模態(tài)影響較小。不考慮液體表面重力波時,ADINA軟件可直接求得結(jié)構(gòu)振動模態(tài)。4)儲液高度相同時,水池半徑越大,同一階次儲液晃動頻率和結(jié)構(gòu)振動頻率越小,儲液晃動振型基本相同,但結(jié)構(gòu)振動振型被激發(fā)的順序不同。5)有限元法和理論法計算儲液晃動頻率相差很小,驗證了所建立有限元模型的可靠性。推導(dǎo)得到了水池儲液晃動基本周期計算公式,可用于水池結(jié)構(gòu)地震響應(yīng)計算。

    [1] Housner GW.Dynamic pressures on accelerated fluid containers[J].Bulletin of the Seismological Society of America,1957,47(1):15-35.

    [2] Haroun MA.Vibration studies and tests of liquid storage tanks[J].Earthquake Engineering and Structural Dynamics,1983,11(2):179-206.

    [3] 居榮初,曾心傳.彈性結(jié)構(gòu)與液體的耦聯(lián)振動理論[M].北京:地震出版社,1983.

    [4] 項忠權(quán),李清林.立式儲罐抗震[M].北京:地震出版社,1990.

    [5] 陳志濤.卵形消化池流—固耦合系統(tǒng)的ADINA有限元分析[D].天津:天津大學(xué),2007.

    [6] 沈世杰.貯水池動水壓力計算探討[J].特種結(jié)構(gòu),1995,12(4):4-7.

    [7] ADINA北京代表處.ADINA土木工程分析實例手冊[M].北京:ADINA代表處,2005.

    [8] 高 霖.地面式鋼筋混凝土水池自愈、滲漏試驗及地震響應(yīng)分析[D].哈爾濱:中國地震局工程力學(xué)研究所,2015.

    [9] 舒亞俐.考慮液固耦合振動的貯液結(jié)構(gòu)動力特性和地震響應(yīng)研究[J].特種結(jié)構(gòu),2013,30(1):69-75.

    The modal analyses considering liquid-structure coupling effect of ground reinforced concrete circular tanks★

    Gao Lin1Wang Mingzhen1,2

    (1.CollegeofArchitectureEngineering,ChongqingUniversityofArtsandSciences,Chongqing402160,China; 2.InstituteofEngineeringMechanics,ChinaEarthquakeAdministration,KeyLaboratoryofEarthquakeEngineeringandEngineeringVibration,ChinaEarthquakeAdministration,Harbin150080,China)

    The modal of ground reinforced concrete tanks are analyzed considering liquid-structure coupling effects. The influence of liquid state, liquid height, tank radius and surface gravity waves on modal are analyzed by using ADINA software. And the theory and finite element calculation results of sloshing frequency are compared. And the sloshing foundation period formula of circular reinforced concrete tank is derived.

    tank, liquid-structure coupling, ADINA, modal analysis, sloshing period

    1009-6825(2017)02-0047-03

    2016-11-08★:重慶文理學(xué)院校級科技項目(項目編號:R2015JJ06);國家自然科學(xué)基金資助項目(項目編號:51678544);國家科技支撐計劃課題(課題編號:2015BAK17B05)

    高 霖(1988- ),女,博士,講師

    P315.8

    A

    猜你喜歡
    液固儲液水池
    實現(xiàn)光參與液固催化界面局域溫度精準(zhǔn)測量
    一種橡膠制品的硫化裝置
    小區(qū)的水池
    適合海上60~90 ℃儲層的液固相變支撐劑相變性能評價
    廣州化工(2020年4期)2020-03-12 08:56:06
    把住醫(yī)?;鹚亻l門
    找水池
    基于ADINA的矩形儲液池流固耦合動力響應(yīng)分析
    山西建筑(2016年25期)2016-09-28 09:59:00
    機載儲液罐研制方法研究
    機械工程師(2014年3期)2014-11-22 02:11:58
    基于CFD的液固兩相流沖刷腐蝕預(yù)測研究
    壓縮機儲液罐大小對空調(diào)整機系統(tǒng)的影響
    家電科技(2013年6期)2013-11-15 07:20:56
    免费一级毛片在线播放高清视频| 国产91精品成人一区二区三区| 精品福利观看| 免费无遮挡裸体视频| 亚洲国产精品合色在线| 成人国产一区最新在线观看| 高清日韩中文字幕在线| 桃红色精品国产亚洲av| 亚洲av免费高清在线观看| 久久久国产成人免费| 一级黄片播放器| 久久草成人影院| 欧美成人a在线观看| 成人性生交大片免费视频hd| 国产激情偷乱视频一区二区| 色哟哟哟哟哟哟| 色噜噜av男人的天堂激情| 无遮挡黄片免费观看| 亚洲av成人av| 成人av在线播放网站| 久久精品国产亚洲av涩爱 | 网址你懂的国产日韩在线| 女人高潮潮喷娇喘18禁视频| 乱人视频在线观看| 欧美丝袜亚洲另类 | 婷婷丁香在线五月| 网址你懂的国产日韩在线| eeuss影院久久| 国产黄色小视频在线观看| 久久午夜亚洲精品久久| 亚洲国产中文字幕在线视频| 波野结衣二区三区在线 | 在线观看日韩欧美| 99久久精品一区二区三区| 亚洲一区二区三区色噜噜| 91麻豆精品激情在线观看国产| 丝袜美腿在线中文| 国产精品久久电影中文字幕| 国产毛片a区久久久久| 18禁黄网站禁片免费观看直播| 国产中年淑女户外野战色| 中文字幕人妻熟人妻熟丝袜美 | 黄片大片在线免费观看| 久久久久久久亚洲中文字幕 | 国产欧美日韩一区二区三| 亚洲国产精品成人综合色| 亚洲人成伊人成综合网2020| 久久久久久久久中文| 中文亚洲av片在线观看爽| 久久这里只有精品中国| 亚洲精华国产精华精| 亚洲av不卡在线观看| 国产亚洲精品一区二区www| 欧美一区二区精品小视频在线| 深夜精品福利| 国模一区二区三区四区视频| av视频在线观看入口| 美女 人体艺术 gogo| 脱女人内裤的视频| 成人国产一区最新在线观看| 超碰av人人做人人爽久久 | 国产探花极品一区二区| 亚洲成人中文字幕在线播放| 欧美黑人巨大hd| 18+在线观看网站| 一本综合久久免费| 色老头精品视频在线观看| 亚洲国产精品久久男人天堂| 精品国产三级普通话版| 国产精品影院久久| 2021天堂中文幕一二区在线观| 看黄色毛片网站| 国产高清videossex| 天堂网av新在线| 欧美日本视频| 亚洲精品一卡2卡三卡4卡5卡| 中文亚洲av片在线观看爽| 日韩欧美国产一区二区入口| 一a级毛片在线观看| 波多野结衣高清作品| 在线观看美女被高潮喷水网站 | 亚洲一区二区三区色噜噜| 丰满人妻一区二区三区视频av | 精品久久久久久久人妻蜜臀av| 精品一区二区三区人妻视频| 亚洲国产欧美人成| 成人av在线播放网站| 人妻丰满熟妇av一区二区三区| 国产一级毛片七仙女欲春2| 岛国视频午夜一区免费看| 日韩欧美精品免费久久 | 又黄又粗又硬又大视频| 精品国产三级普通话版| 男女视频在线观看网站免费| 成人鲁丝片一二三区免费| 成人永久免费在线观看视频| 男人舔女人下体高潮全视频| 亚洲精华国产精华精| 69av精品久久久久久| 成人午夜高清在线视频| 一个人观看的视频www高清免费观看| 男女视频在线观看网站免费| 给我免费播放毛片高清在线观看| 首页视频小说图片口味搜索| 精品久久久久久久久久免费视频| 国产精品亚洲美女久久久| 日本黄色视频三级网站网址| 国产av不卡久久| 母亲3免费完整高清在线观看| 亚洲五月婷婷丁香| 日本黄色片子视频| 亚洲人成网站高清观看| 岛国在线观看网站| 亚洲美女黄片视频| 麻豆国产97在线/欧美| 国产精品女同一区二区软件 | 偷拍熟女少妇极品色| 无人区码免费观看不卡| 亚洲中文日韩欧美视频| 国产伦精品一区二区三区视频9 | 一级黄色大片毛片| 国产精品三级大全| 成人亚洲精品av一区二区| 亚洲精品一卡2卡三卡4卡5卡| 精品乱码久久久久久99久播| 精品无人区乱码1区二区| 国产伦在线观看视频一区| 中文字幕人成人乱码亚洲影| 日本 欧美在线| 九色国产91popny在线| 一区二区三区高清视频在线| 日本精品一区二区三区蜜桃| 99在线视频只有这里精品首页| 精品久久久久久久人妻蜜臀av| 色在线成人网| 欧美黑人巨大hd| xxx96com| 亚洲 国产 在线| 国产精品精品国产色婷婷| 中出人妻视频一区二区| 欧美区成人在线视频| 99精品欧美一区二区三区四区| 成人av在线播放网站| 国产伦在线观看视频一区| 国产主播在线观看一区二区| 久99久视频精品免费| 国产亚洲精品一区二区www| 欧美又色又爽又黄视频| 最新中文字幕久久久久| 好男人电影高清在线观看| 一区二区三区高清视频在线| 亚洲av成人不卡在线观看播放网| 亚洲,欧美精品.| 51午夜福利影视在线观看| 无人区码免费观看不卡| 精品久久久久久久久久免费视频| 老司机在亚洲福利影院| 亚洲精品日韩av片在线观看 | 国产麻豆成人av免费视频| 97超级碰碰碰精品色视频在线观看| 免费人成在线观看视频色| 欧美区成人在线视频| 亚洲精品成人久久久久久| 久久久国产精品麻豆| 女警被强在线播放| 丁香六月欧美| eeuss影院久久| 18禁裸乳无遮挡免费网站照片| 国产av麻豆久久久久久久| 久久国产乱子伦精品免费另类| 欧美xxxx黑人xx丫x性爽| 窝窝影院91人妻| 亚洲成人中文字幕在线播放| 别揉我奶头~嗯~啊~动态视频| 韩国av一区二区三区四区| 精品国产三级普通话版| 97超级碰碰碰精品色视频在线观看| 欧美日韩精品网址| 九九热线精品视视频播放| 色av中文字幕| 国产爱豆传媒在线观看| 男女之事视频高清在线观看| 亚洲精品在线观看二区| 九九热线精品视视频播放| 亚洲avbb在线观看| 好男人电影高清在线观看| 丝袜美腿在线中文| 一级毛片女人18水好多| 草草在线视频免费看| 欧美成狂野欧美在线观看| 亚洲一区二区三区不卡视频| 免费高清视频大片| 亚洲久久久久久中文字幕| 1000部很黄的大片| 欧美成人一区二区免费高清观看| 国产激情偷乱视频一区二区| 在线视频色国产色| 18禁裸乳无遮挡免费网站照片| 9191精品国产免费久久| 久久久久久久久久黄片| 午夜福利成人在线免费观看| 亚洲精品在线观看二区| 成年人黄色毛片网站| 久久久久免费精品人妻一区二区| 精品国产美女av久久久久小说| 狂野欧美白嫩少妇大欣赏| 欧美日韩瑟瑟在线播放| 久久中文看片网| 97超级碰碰碰精品色视频在线观看| 特级一级黄色大片| 日韩欧美三级三区| 内地一区二区视频在线| 午夜免费成人在线视频| 免费观看的影片在线观看| 午夜日韩欧美国产| 免费高清视频大片| 国产精品久久电影中文字幕| 村上凉子中文字幕在线| 一本一本综合久久| 在线a可以看的网站| 国产亚洲欧美在线一区二区| 在线播放国产精品三级| 搡老岳熟女国产| 国模一区二区三区四区视频| 少妇熟女aⅴ在线视频| 99精品欧美一区二区三区四区| 免费看日本二区| 国产一区二区三区视频了| 制服丝袜大香蕉在线| 日本三级黄在线观看| 成年女人看的毛片在线观看| 日本免费a在线| 最近最新中文字幕大全电影3| 欧美又色又爽又黄视频| 国产淫片久久久久久久久 | 成人特级黄色片久久久久久久| 最近视频中文字幕2019在线8| www.熟女人妻精品国产| 一级a爱片免费观看的视频| tocl精华| 亚洲va日本ⅴa欧美va伊人久久| 天堂√8在线中文| svipshipincom国产片| 国内少妇人妻偷人精品xxx网站| 欧美黑人巨大hd| 亚洲午夜理论影院| 国产欧美日韩精品亚洲av| 97超级碰碰碰精品色视频在线观看| 亚洲中文字幕日韩| 精品久久久久久久末码| 亚洲va日本ⅴa欧美va伊人久久| 九色成人免费人妻av| 欧洲精品卡2卡3卡4卡5卡区| 中文字幕av成人在线电影| 中文字幕熟女人妻在线| 国产日本99.免费观看| 国产精品99久久久久久久久| 国产一区二区亚洲精品在线观看| 99热这里只有是精品50| 免费av毛片视频| 最新在线观看一区二区三区| 亚洲av成人不卡在线观看播放网| 成人av在线播放网站| 又紧又爽又黄一区二区| 日日干狠狠操夜夜爽| 国产av麻豆久久久久久久| 国产精品久久久久久精品电影| 熟女人妻精品中文字幕| 成人特级黄色片久久久久久久| 岛国在线免费视频观看| 亚洲中文日韩欧美视频| 久久精品影院6| 亚洲第一电影网av| 男人舔女人下体高潮全视频| 免费看日本二区| 好男人电影高清在线观看| 黄片大片在线免费观看| 女同久久另类99精品国产91| 欧美在线一区亚洲| 国内精品一区二区在线观看| 亚洲av一区综合| 国产精品电影一区二区三区| 国产高清视频在线观看网站| 日本一本二区三区精品| 久久精品国产亚洲av香蕉五月| 免费在线观看成人毛片| 亚洲 国产 在线| 禁无遮挡网站| 小说图片视频综合网站| 久久精品91无色码中文字幕| 亚洲av二区三区四区| 香蕉丝袜av| 成人精品一区二区免费| 九色成人免费人妻av| 一本精品99久久精品77| 国产主播在线观看一区二区| 熟女少妇亚洲综合色aaa.| 亚洲 国产 在线| 精品欧美国产一区二区三| 手机成人av网站| 免费观看人在逋| 久久久久久久久中文| 国产欧美日韩精品亚洲av| 天堂影院成人在线观看| 97人妻精品一区二区三区麻豆| 熟女人妻精品中文字幕| 久久久久久人人人人人| 成人av在线播放网站| 成年女人看的毛片在线观看| 青草久久国产| 日韩中文字幕欧美一区二区| 欧美黑人巨大hd| 熟妇人妻久久中文字幕3abv| 露出奶头的视频| 欧美一级a爱片免费观看看| 国产欧美日韩一区二区三| 色av中文字幕| 久久国产乱子伦精品免费另类| 午夜激情福利司机影院| 波多野结衣高清无吗| 国产欧美日韩一区二区精品| 国内精品美女久久久久久| a在线观看视频网站| 99国产精品一区二区三区| 女生性感内裤真人,穿戴方法视频| 久久人妻av系列| 久久亚洲精品不卡| 日本三级黄在线观看| 午夜福利欧美成人| 特大巨黑吊av在线直播| 神马国产精品三级电影在线观看| bbb黄色大片| 国产视频一区二区在线看| 亚洲欧美日韩卡通动漫| 好看av亚洲va欧美ⅴa在| 19禁男女啪啪无遮挡网站| 亚洲在线观看片| 精品乱码久久久久久99久播| 日韩欧美国产在线观看| 亚洲成人中文字幕在线播放| 婷婷丁香在线五月| 最新中文字幕久久久久| 亚洲精品乱码久久久v下载方式 | 久久国产精品人妻蜜桃| 国产成年人精品一区二区| 国产aⅴ精品一区二区三区波| 久久99热这里只有精品18| 成熟少妇高潮喷水视频| 一a级毛片在线观看| 一本一本综合久久| 欧美精品啪啪一区二区三区| 一级毛片高清免费大全| 亚洲国产精品合色在线| 九色成人免费人妻av| 国产日本99.免费观看| 黄色成人免费大全| 少妇的逼好多水| 美女高潮的动态| 日本一二三区视频观看| 欧美成人免费av一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 99热只有精品国产| 色视频www国产| 亚洲电影在线观看av| 制服丝袜大香蕉在线| 日本三级黄在线观看| 国产精品久久久久久亚洲av鲁大| 91久久精品国产一区二区成人 | 国产伦在线观看视频一区| 精品乱码久久久久久99久播| 亚洲成a人片在线一区二区| 又紧又爽又黄一区二区| 亚洲狠狠婷婷综合久久图片| 成人18禁在线播放| 又粗又爽又猛毛片免费看| 日本 欧美在线| a级一级毛片免费在线观看| 18禁黄网站禁片免费观看直播| 18禁裸乳无遮挡免费网站照片| 一级毛片高清免费大全| 欧美午夜高清在线| 狂野欧美激情性xxxx| 午夜福利在线观看免费完整高清在 | 人妻久久中文字幕网| 日本免费一区二区三区高清不卡| 午夜激情福利司机影院| 久久人人精品亚洲av| 国产精品亚洲av一区麻豆| 国产三级在线视频| e午夜精品久久久久久久| 婷婷精品国产亚洲av在线| av片东京热男人的天堂| 欧美在线黄色| 国产成人欧美在线观看| 欧美在线黄色| 成人18禁在线播放| 五月玫瑰六月丁香| 免费观看精品视频网站| 综合色av麻豆| 99国产精品一区二区三区| 欧美激情在线99| 激情在线观看视频在线高清| 18美女黄网站色大片免费观看| 女同久久另类99精品国产91| 十八禁人妻一区二区| 又黄又粗又硬又大视频| 舔av片在线| 午夜视频国产福利| 国产精品一区二区三区四区久久| 欧美日韩亚洲国产一区二区在线观看| 久久精品夜夜夜夜夜久久蜜豆| 亚洲男人的天堂狠狠| 色综合亚洲欧美另类图片| 国产精品永久免费网站| 一级作爱视频免费观看| 午夜福利在线在线| www国产在线视频色| 国产不卡一卡二| 在线播放国产精品三级| 天堂√8在线中文| 中文字幕高清在线视频| 哪里可以看免费的av片| 99久久精品一区二区三区| 久久久国产精品麻豆| 好看av亚洲va欧美ⅴa在| 亚洲七黄色美女视频| 亚洲av二区三区四区| 黄片大片在线免费观看| 91字幕亚洲| xxxwww97欧美| 国产真人三级小视频在线观看| 亚洲成人久久爱视频| 成熟少妇高潮喷水视频| 国产亚洲精品综合一区在线观看| 日本精品一区二区三区蜜桃| 午夜激情福利司机影院| 国产午夜福利久久久久久| 欧美在线一区亚洲| 韩国av一区二区三区四区| 亚洲片人在线观看| 国产成人av教育| 成人特级黄色片久久久久久久| 69av精品久久久久久| 特级一级黄色大片| av视频在线观看入口| 精品不卡国产一区二区三区| 亚洲欧美一区二区三区国产| 国产真实伦视频高清在线观看| 国产成人aa在线观看| 街头女战士在线观看网站| 亚洲人成网站在线播| av免费在线看不卡| 99热这里只有是精品50| 最近视频中文字幕2019在线8| 久久久久国产网址| 毛片女人毛片| 久久精品人妻少妇| 国产精品蜜桃在线观看| 七月丁香在线播放| 欧美变态另类bdsm刘玥| 女人十人毛片免费观看3o分钟| av福利片在线观看| 老司机影院成人| 神马国产精品三级电影在线观看| 狠狠精品人妻久久久久久综合| 精品久久久久久电影网| 大又大粗又爽又黄少妇毛片口| 久久久精品欧美日韩精品| 日韩国内少妇激情av| 能在线免费看毛片的网站| 日韩欧美 国产精品| 麻豆av噜噜一区二区三区| 少妇被粗大猛烈的视频| 精品久久久噜噜| 国产视频首页在线观看| 国产白丝娇喘喷水9色精品| 大香蕉97超碰在线| 自拍偷自拍亚洲精品老妇| 欧美成人午夜免费资源| 亚洲内射少妇av| 国产精品99久久久久久久久| 亚洲va在线va天堂va国产| 天天躁夜夜躁狠狠久久av| 亚洲精品456在线播放app| 亚洲国产精品专区欧美| 黄色欧美视频在线观看| 国产综合精华液| 免费av不卡在线播放| 五月天丁香电影| 国产精品福利在线免费观看| 日韩欧美国产在线观看| 日韩强制内射视频| 国产精品1区2区在线观看.| 在线免费十八禁| 91av网一区二区| 丝袜喷水一区| 99九九线精品视频在线观看视频| 成人亚洲精品一区在线观看 | 又大又黄又爽视频免费| 精品久久久噜噜| 一级毛片久久久久久久久女| 国产有黄有色有爽视频| 麻豆精品久久久久久蜜桃| 99久久精品热视频| 亚洲人成网站在线观看播放| 毛片一级片免费看久久久久| 午夜福利网站1000一区二区三区| 人妻系列 视频| 99久久精品热视频| 中文字幕av在线有码专区| av在线播放精品| 欧美xxⅹ黑人| 亚洲精品456在线播放app| 97超视频在线观看视频| 高清视频免费观看一区二区 | 亚洲精品日韩在线中文字幕| 欧美成人午夜免费资源| 一级毛片电影观看| 国产精品一区二区在线观看99 | 全区人妻精品视频| 青春草国产在线视频| 亚洲综合色惰| 在线播放无遮挡| 欧美成人a在线观看| 久久热精品热| 高清毛片免费看| 97在线视频观看| 亚洲经典国产精华液单| 欧美xxxx性猛交bbbb| 热99在线观看视频| 国产欧美日韩精品一区二区| 91精品伊人久久大香线蕉| 亚洲精品aⅴ在线观看| 欧美丝袜亚洲另类| 精品一区二区免费观看| 国产激情偷乱视频一区二区| 大片免费播放器 马上看| 国产永久视频网站| 色尼玛亚洲综合影院| av一本久久久久| 国产精品一区二区在线观看99 | 亚洲欧美精品自产自拍| 国产一级毛片在线| 亚洲18禁久久av| 国产一区二区三区综合在线观看 | 国精品久久久久久国模美| 日本-黄色视频高清免费观看| 女人被狂操c到高潮| xxx大片免费视频| 久久鲁丝午夜福利片| 久久久久精品性色| 搡女人真爽免费视频火全软件| 久久这里有精品视频免费| 国产精品久久久久久av不卡| 婷婷色麻豆天堂久久| 精品人妻一区二区三区麻豆| 狂野欧美激情性xxxx在线观看| 网址你懂的国产日韩在线| 免费黄色在线免费观看| 99久久精品热视频| 中文资源天堂在线| 99久国产av精品国产电影| 亚洲精品成人久久久久久| av专区在线播放| 午夜福利视频精品| 2022亚洲国产成人精品| 国产高清三级在线| 国产高清有码在线观看视频| ponron亚洲| 久久精品国产鲁丝片午夜精品| 最近视频中文字幕2019在线8| 亚洲国产高清在线一区二区三| 亚洲精品456在线播放app| 高清av免费在线| 一个人看的www免费观看视频| 国产精品久久久久久精品电影小说 | 尾随美女入室| 在现免费观看毛片| 精品国产三级普通话版| 在线播放无遮挡| 国内精品宾馆在线| 精华霜和精华液先用哪个| 国产激情偷乱视频一区二区| 亚洲国产最新在线播放| 中文欧美无线码| 青春草视频在线免费观看| 少妇猛男粗大的猛烈进出视频 | av一本久久久久| 国产精品久久久久久久电影| 成人二区视频| 亚洲国产欧美人成| 精品国产露脸久久av麻豆 | av网站免费在线观看视频 | av免费观看日本| 久久久精品免费免费高清| 一区二区三区免费毛片| 亚洲欧美日韩无卡精品| 日韩一区二区视频免费看| 最后的刺客免费高清国语| 91av网一区二区| 欧美变态另类bdsm刘玥| 春色校园在线视频观看| 男人舔女人下体高潮全视频| kizo精华| 老司机影院成人| 亚洲国产av新网站| 国产亚洲精品av在线| 亚洲最大成人手机在线| 3wmmmm亚洲av在线观看| 18+在线观看网站| 精品欧美国产一区二区三| 亚洲国产av新网站| 亚洲在久久综合| 简卡轻食公司| 色吧在线观看|