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

        ?

        跳躍荷載作用下大跨結(jié)構(gòu)動(dòng)力響應(yīng)分析及振動(dòng)控制

        2019-02-28 07:33:32歐陽國(guó)泉呂亞兵
        科學(xué)與財(cái)富 2019年1期
        關(guān)鍵詞:有限元模擬

        歐陽國(guó)泉 呂亞兵

        摘 要:本文系統(tǒng)地介紹了樓蓋結(jié)構(gòu)的跳躍荷載模型,并采用ANSYS有限元方法進(jìn)行樓蓋模擬。收集他人實(shí)測(cè)與模型模擬的進(jìn)展,對(duì)有限元模型的準(zhǔn)確性進(jìn)行研究,總結(jié)出舒適度的評(píng)價(jià)標(biāo)準(zhǔn)。最后,提出大跨樓蓋結(jié)構(gòu)人致振動(dòng)響應(yīng)需要進(jìn)一步研究的問題。

        關(guān)鍵詞:人致振動(dòng);跳躍荷載;大跨樓蓋;有限元模擬;荷載模型

        0 引言

        20世紀(jì)以來,建筑結(jié)構(gòu)形式日益優(yōu)化,大跨樓蓋結(jié)構(gòu)開始偏向于跨度大、自重輕、基頻低等形式,在人的活動(dòng)荷載作用下會(huì)產(chǎn)生較大豎向振動(dòng),甚至發(fā)生共振,對(duì)結(jié)構(gòu)的使用功能造成不利的影響。一些研究已經(jīng)在國(guó)內(nèi)外及有關(guān)舒適性控制標(biāo)準(zhǔn)完成已經(jīng)制定,如表1所示[1]。一般樓板的固有頻率均能保證滿足該規(guī)范,同時(shí),振動(dòng)加速度值應(yīng)不超過0.07m·s-2 。

        0 引言

        20世紀(jì)以來,建筑結(jié)構(gòu)形式日益優(yōu)化,大跨樓蓋結(jié)構(gòu)開始偏向于跨度大、自重輕、基頻低等形式,在人的活動(dòng)荷載作用下會(huì)產(chǎn)生較大豎向振動(dòng),甚至發(fā)生共振,對(duì)結(jié)構(gòu)的使用功能造成不利的影響。一些研究已經(jīng)在國(guó)內(nèi)外及有關(guān)舒適性控制標(biāo)準(zhǔn)完成已經(jīng)制定,如表1所示[1]。一般樓板的固有頻率均能保證滿足該規(guī)范,同時(shí),振動(dòng)加速度值應(yīng)不超過0.07m·s-2 。

        1 跳躍荷載模型

        1.1 確定性荷載模型:?jiǎn)未翁S荷載模型

        單次跳躍荷載模型是將原地跳躍視作單次跳躍的周期性重復(fù)過程,只需要給出單次跳躍的荷載曲線,其核心是單次脈沖曲線形狀,即一次躍的荷載曲線形狀[4]。半正弦模型將單次跳躍視作半正弦[5],可表示為

        式中,Kp=Fmax/G ,為脈沖系數(shù)或峰值因子;Fmax 為跳躍荷載的峰值;G為跳躍著靜體重;tp為腳與地面接觸時(shí)間;Tp 為跳躍周期。

        1.2 確定性荷載模型:傅里葉級(jí)數(shù)模型

        將半正弦模型進(jìn)一步用傅里葉級(jí)數(shù)的形式表示,可以表示為:

        式中,F(xiàn)(t) 為tt時(shí)刻的跳躍力;G為體重;rn 為n階傅里葉系數(shù),亦為動(dòng)載因子(DLFn);φn 為相位角。目前,大量的文獻(xiàn)針對(duì)以上公式中的動(dòng)力荷載因子(DLFs)進(jìn)行了系統(tǒng)性的研究。

        2 跳躍荷載試驗(yàn)

        常見的跳躍荷載試驗(yàn)一般分為直接測(cè)試和間接測(cè)試。直接測(cè)試方法是采用測(cè)力板研究人行走、跳躍、跑步時(shí)產(chǎn)生的作用力,跳躍的固定作用點(diǎn)點(diǎn)比較特殊,所以對(duì)于有剛性的的測(cè)力板可以發(fā)揮它的優(yōu)勢(shì),即測(cè)試者在比較接近自然步行的狀態(tài)下進(jìn)行試驗(yàn),從而可以直接了當(dāng)?shù)墨@得連續(xù)荷載記錄,測(cè)得的數(shù)據(jù)精度也較高。然而由于測(cè)力板位置固定,測(cè)試者往往為了準(zhǔn)確踏上測(cè)力板而主動(dòng)調(diào)整步態(tài),形成所謂“非自然”行走狀態(tài),造成測(cè)試結(jié)果失真。此類實(shí)驗(yàn)同時(shí)缺少對(duì)實(shí)驗(yàn)者的身體運(yùn)動(dòng)信息的測(cè)試,例如人體重心上下波動(dòng),跳躍頻率及節(jié)拍無法得到統(tǒng)一,所以也只能進(jìn)行粗略的估計(jì)。

        3 有限元模擬

        3.1 有限元模型的建立

        在進(jìn)行對(duì)大跨度樓蓋結(jié)構(gòu)的荷載模擬計(jì)算時(shí),首先要建立合適的樓蓋模型。由于實(shí)際情況中可能對(duì)模型的簡(jiǎn)化并不能保證一定的合理性,比如,邊緣條件的等效不合理,對(duì)構(gòu)件的過度簡(jiǎn)化,造成有限元計(jì)算結(jié)果與實(shí)際情況不匹配。在1994年研究中就有前人指出樓板響應(yīng)主要受基波模式的貢獻(xiàn)大小的影響,而在對(duì)稱動(dòng)荷載(如跳躍等)作用下,具有簡(jiǎn)支梁邊界條件的結(jié)構(gòu)是不可忽略的。因此,只需要考慮基本模式的響應(yīng)。即高階級(jí)數(shù)的能量所引起的舒適度也可以暫時(shí)不用考慮。

        3.2 有限元模型優(yōu)化及精確化

        在進(jìn)行土木工程領(lǐng)域的研究中,混凝土開裂及彈性模量、剛度等都是非確定性因素,都存在著不精確的誤差?,F(xiàn)如今,也有許多學(xué)者在進(jìn)行模型設(shè)計(jì)以及計(jì)算時(shí),對(duì)模型也進(jìn)行了許多的優(yōu)化與改進(jìn)。他們結(jié)合了現(xiàn)代計(jì)算機(jī)技術(shù)使得有限元模型的精確度有了很大的提高。此外,國(guó)內(nèi)外學(xué)者對(duì)模型邊界條件的約束也做了大量的研究。比如,當(dāng)邊界條件由簡(jiǎn)單支承改為固定支承時(shí),由于剛度的增加,位移和加速度響應(yīng)大大減小,加速度響應(yīng)大大減小,位移響應(yīng)略有減小。在邊界條件固定的情況下,樓板的振動(dòng)響應(yīng)比簡(jiǎn)單支承的振動(dòng)響應(yīng)下降快得多,樓板結(jié)構(gòu)的振動(dòng)規(guī)律更加穩(wěn)定。邊界條件簡(jiǎn)支時(shí),振動(dòng)具有明顯的三級(jí)特性,這是大跨度現(xiàn)澆后張樓板的共同特征。

        4 模型的現(xiàn)場(chǎng)實(shí)測(cè)動(dòng)力特性

        4.1 游泳館樓蓋實(shí)測(cè)

        朱萍對(duì)某游泳館進(jìn)行了實(shí)地檢測(cè)和有限元模擬分析,現(xiàn)場(chǎng)案列如下,該游泳館為混凝土框架結(jié)構(gòu),總共有3層,整體結(jié)構(gòu)總高16.52m。游泳館1樓層高為3.9m,二樓樓蓋為大跨度預(yù)應(yīng)力次梁樓蓋,二樓結(jié)構(gòu)存在夾層,夾層層高為4.5m,結(jié)構(gòu)層高為3.6m.,,3樓結(jié)構(gòu)層高4.5m。

        通過對(duì)誤差的產(chǎn)生進(jìn)行分析,提出在舒適性評(píng)價(jià)的動(dòng)力計(jì)算的基礎(chǔ)上,一般非承重構(gòu)件(如填充墻)及其與主體結(jié)構(gòu)相連接的構(gòu)件通常不會(huì)斷裂破壞。實(shí)際情況中匹配條件的限制,少量的輕量級(jí)加密墻是建立在大跨度區(qū)域和少量的木模板是分散在大跨度的地板上,導(dǎo)致測(cè)量大跨度樓板的自然振動(dòng)頻率略高于相同頻率的仿真分析。因此,大跨度樓板實(shí)測(cè)固有振動(dòng)頻率略高于有限元模擬分析,相對(duì)誤差約為5%。

        5 人致振動(dòng)舒適度評(píng)價(jià)指標(biāo)

        一般來說,人致振動(dòng)舒適度與地板振動(dòng)的固有頻率、人們生活的環(huán)境、人們自身活動(dòng)的狀態(tài)以及人們的心理反應(yīng)有關(guān)。為保證梁板結(jié)構(gòu)在正常使用條件下的變形要求,國(guó)家規(guī)范一般規(guī)定了彎構(gòu)件截面的撓度跨比和高度跨比的最小值。同時(shí),針對(duì)基頻較低類型的樓板結(jié)構(gòu),須考慮由于人的走、跑等運(yùn)動(dòng)帶來的結(jié)構(gòu)振動(dòng)舒適度問題。

        對(duì)于以上問題,目前有大量學(xué)者已近提出了基于不同理論的舒適度指標(biāo),然而樓蓋在動(dòng)力荷載作用下的加速度峰值指標(biāo)依然在結(jié)構(gòu)設(shè)計(jì)中應(yīng)用得較多且較為成熟。

        在結(jié)構(gòu)自振頻率低于4Hz時(shí),兩種規(guī)范的峰值加速度差距很大,而在4Hz時(shí),兩者的峰值加速度限值基本保持一致。

        通??紤]人行走作用下梁板結(jié)構(gòu)的豎向振動(dòng)響應(yīng),結(jié)構(gòu)豎向振動(dòng)的總體頻率計(jì)權(quán)函數(shù)可以表達(dá)為:

        (3)

        式中f是振動(dòng)頻率。對(duì)于一個(gè)頻率為f,峰因子為_ ,峰值加速度為

        的振動(dòng)信號(hào),其振動(dòng)強(qiáng)度 為:

        (4)

        文獻(xiàn)[5]中建議取為2.5。

        6 結(jié)語

        近些年來,自國(guó)內(nèi)外諸多學(xué)者對(duì)大跨樓蓋的豎向振動(dòng)響應(yīng)及動(dòng)力特性的大量研究之后,本綜述打算提出如下幾條總結(jié)與建議:

        (1)對(duì)于有限元模擬,國(guó)內(nèi)外學(xué)者對(duì)模型的參數(shù)設(shè)置優(yōu)化包括網(wǎng)格細(xì)分化、材料選取精確化非結(jié)構(gòu)模型的精細(xì)化發(fā)展做的越來越詳細(xì)。但在對(duì)于模型計(jì)算時(shí),對(duì)于邊界條件的微分方程的計(jì)算選取還略有不足,仍需要進(jìn)一步的研究。

        (2)關(guān)于跳躍荷載的傅里葉級(jí)數(shù)模型中模型階數(shù)的選取研究較少,大多數(shù)認(rèn)為3階足以計(jì)算結(jié)構(gòu)的振動(dòng)響應(yīng),但往往卻忽略了結(jié)果的精確性。應(yīng)根據(jù)樓板的振動(dòng)響應(yīng)貢獻(xiàn)的大小來確定模型階數(shù)的選擇,做到簡(jiǎn)便化與精確化。

        參考文獻(xiàn):

        [1] 葉正強(qiáng),李愛群,丁幼亮,等.某大跨人行天橋的消能減振設(shè)計(jì)(一)[J].特種結(jié)構(gòu),2003,20(1): 68-70.

        [2] 于敬海,曹建東,李敬明. 跳躍荷載作用下大跨度鋼筋混凝土樓蓋舒適度分析[J]. 工程抗震與加固改造. 2013(06):64-67.

        [3] 朱萍. 跳躍荷載作用下的大跨樓蓋振動(dòng)響應(yīng)與控制研究[D]. 南京:東南大學(xué),2015.

        [4] 陳雋. 人致荷載與人致結(jié)構(gòu)振動(dòng)[M]. 北京:科學(xué)出版社,2016.

        [5] G.Pernica. Dynamic load factors for pedestrian movements and rhythmic exercises[J]. Canadian Acoustics. 2015,18(2):3-18.

        1 跳躍荷載模型

        1.1 確定性荷載模型:?jiǎn)未翁S荷載模型

        單次跳躍荷載模型是將原地跳躍視作單次跳躍的周期性重復(fù)過程,只需要給出單次跳躍的荷載曲線,其核心是單次脈沖曲線形狀,即一次躍的荷載曲線形狀[4]。半正弦模型將單次跳躍視作半正弦[5],可表示為

        式中,Kp=Fmax/G ,為脈沖系數(shù)或峰值因子;Fmax 為跳躍荷載的峰值;G為跳躍著靜體重;tp為腳與地面接觸時(shí)間;Tp 為跳躍周期。

        1.2 確定性荷載模型:傅里葉級(jí)數(shù)模型

        將半正弦模型進(jìn)一步用傅里葉級(jí)數(shù)的形式表示,可以表示為:

        式中,F(xiàn)(t) 為tt時(shí)刻的跳躍力;G為體重;rn 為n階傅里葉系數(shù),亦為動(dòng)載因子(DLFn);φn 為相位角。目前,大量的文獻(xiàn)針對(duì)以上公式中的動(dòng)力荷載因子(DLFs)進(jìn)行了系統(tǒng)性的研究。

        2 跳躍荷載試驗(yàn)

        常見的跳躍荷載試驗(yàn)一般分為直接測(cè)試和間接測(cè)試。直接測(cè)試方法是采用測(cè)力板研究人行走、跳躍、跑步時(shí)產(chǎn)生的作用力,跳躍的固定作用點(diǎn)點(diǎn)比較特殊,所以對(duì)于有剛性的的測(cè)力板可以發(fā)揮它的優(yōu)勢(shì),即測(cè)試者在比較接近自然步行的狀態(tài)下進(jìn)行試驗(yàn),從而可以直接了當(dāng)?shù)墨@得連續(xù)荷載記錄,測(cè)得的數(shù)據(jù)精度也較高。然而由于測(cè)力板位置固定,測(cè)試者往往為了準(zhǔn)確踏上測(cè)力板而主動(dòng)調(diào)整步態(tài),形成所謂“非自然”行走狀態(tài),造成測(cè)試結(jié)果失真。此類實(shí)驗(yàn)同時(shí)缺少對(duì)實(shí)驗(yàn)者的身體運(yùn)動(dòng)信息的測(cè)試,例如人體重心上下波動(dòng),跳躍頻率及節(jié)拍無法得到統(tǒng)一,所以也只能進(jìn)行粗略的估計(jì)。

        3 有限元模擬

        3.1 有限元模型的建立

        在進(jìn)行對(duì)大跨度樓蓋結(jié)構(gòu)的荷載模擬計(jì)算時(shí),首先要建立合適的樓蓋模型。由于實(shí)際情況中可能對(duì)模型的簡(jiǎn)化并不能保證一定的合理性,比如,邊緣條件的等效不合理,對(duì)構(gòu)件的過度簡(jiǎn)化,造成有限元計(jì)算結(jié)果與實(shí)際情況不匹配。在1994年研究中就有前人指出樓板響應(yīng)主要受基波模式的貢獻(xiàn)大小的影響,而在對(duì)稱動(dòng)荷載(如跳躍等)作用下,具有簡(jiǎn)支梁邊界條件的結(jié)構(gòu)是不可忽略的。因此,只需要考慮基本模式的響應(yīng)。即高階級(jí)數(shù)的能量所引起的舒適度也可以暫時(shí)不用考慮。

        3.2 有限元模型優(yōu)化及精確化

        在進(jìn)行土木工程領(lǐng)域的研究中,混凝土開裂及彈性模量、剛度等都是非確定性因素,都存在著不精確的誤差?,F(xiàn)如今,也有許多學(xué)者在進(jìn)行模型設(shè)計(jì)以及計(jì)算時(shí),對(duì)模型也進(jìn)行了許多的優(yōu)化與改進(jìn)。他們結(jié)合了現(xiàn)代計(jì)算機(jī)技術(shù)使得有限元模型的精確度有了很大的提高。此外,國(guó)內(nèi)外學(xué)者對(duì)模型邊界條件的約束也做了大量的研究。比如,當(dāng)邊界條件由簡(jiǎn)單支承改為固定支承時(shí),由于剛度的增加,位移和加速度響應(yīng)大大減小,加速度響應(yīng)大大減小,位移響應(yīng)略有減小。在邊界條件固定的情況下,樓板的振動(dòng)響應(yīng)比簡(jiǎn)單支承的振動(dòng)響應(yīng)下降快得多,樓板結(jié)構(gòu)的振動(dòng)規(guī)律更加穩(wěn)定。邊界條件簡(jiǎn)支時(shí),振動(dòng)具有明顯的三級(jí)特性,這是大跨度現(xiàn)澆后張樓板的共同特征。

        4 模型的現(xiàn)場(chǎng)實(shí)測(cè)動(dòng)力特性

        4.1 游泳館樓蓋實(shí)測(cè)

        朱萍對(duì)某游泳館進(jìn)行了實(shí)地檢測(cè)和有限元模擬分析,現(xiàn)場(chǎng)案列如下,該游泳館為混凝土框架結(jié)構(gòu),總共有3層,整體結(jié)構(gòu)總高16.52m。游泳館1樓層高為3.9m,二樓樓蓋為大跨度預(yù)應(yīng)力次梁樓蓋,二樓結(jié)構(gòu)存在夾層,夾層層高為4.5m,結(jié)構(gòu)層高為3.6m.,,3樓結(jié)構(gòu)層高4.5m。

        通過對(duì)誤差的產(chǎn)生進(jìn)行分析,提出在舒適性評(píng)價(jià)的動(dòng)力計(jì)算的基礎(chǔ)上,一般非承重構(gòu)件(如填充墻)及其與主體結(jié)構(gòu)相連接的構(gòu)件通常不會(huì)斷裂破壞。實(shí)際情況中匹配條件的限制,少量的輕量級(jí)加密墻是建立在大跨度區(qū)域和少量的木模板是分散在大跨度的地板上,導(dǎo)致測(cè)量大跨度樓板的自然振動(dòng)頻率略高于相同頻率的仿真分析。因此,大跨度樓板實(shí)測(cè)固有振動(dòng)頻率略高于有限元模擬分析,相對(duì)誤差約為5%。

        5 人致振動(dòng)舒適度評(píng)價(jià)指標(biāo)

        一般來說,人致振動(dòng)舒適度與地板振動(dòng)的固有頻率、人們生活的環(huán)境、人們自身活動(dòng)的狀態(tài)以及人們的心理反應(yīng)有關(guān)。為保證梁板結(jié)構(gòu)在正常使用條件下的變形要求,國(guó)家規(guī)范一般規(guī)定了彎構(gòu)件截面的撓度跨比和高度跨比的最小值。同時(shí),針對(duì)基頻較低類型的樓板結(jié)構(gòu),須考慮由于人的走、跑等運(yùn)動(dòng)帶來的結(jié)構(gòu)振動(dòng)舒適度問題。

        對(duì)于以上問題,目前有大量學(xué)者已近提出了基于不同理論的舒適度指標(biāo),然而樓蓋在動(dòng)力荷載作用下的加速度峰值指標(biāo)依然在結(jié)構(gòu)設(shè)計(jì)中應(yīng)用得較多且較為成熟。

        在結(jié)構(gòu)自振頻率低于4Hz時(shí),兩種規(guī)范的峰值加速度差距很大,而在4Hz時(shí),兩者的峰值加速度限值基本保持一致。

        通常考慮人行走作用下梁板結(jié)構(gòu)的豎向振動(dòng)響應(yīng),結(jié)構(gòu)豎向振動(dòng)的總體頻率計(jì)權(quán)函數(shù)可以表達(dá)為:

        (3)

        式中f是振動(dòng)頻率。對(duì)于一個(gè)頻率為f,峰因子為_ ,峰值加速度為

        的振動(dòng)信號(hào),其振動(dòng)強(qiáng)度 為:

        (4)

        文獻(xiàn)[5]中建議取為2.5。

        6 結(jié)語

        近些年來,自國(guó)內(nèi)外諸多學(xué)者對(duì)大跨樓蓋的豎向振動(dòng)響應(yīng)及動(dòng)力特性的大量研究之后,本綜述打算提出如下幾條總結(jié)與建議:

        (1)對(duì)于有限元模擬,國(guó)內(nèi)外學(xué)者對(duì)模型的參數(shù)設(shè)置優(yōu)化包括網(wǎng)格細(xì)分化、材料選取精確化非結(jié)構(gòu)模型的精細(xì)化發(fā)展做的越來越詳細(xì)。但在對(duì)于模型計(jì)算時(shí),對(duì)于邊界條件的微分方程的計(jì)算選取還略有不足,仍需要進(jìn)一步的研究。

        (2)關(guān)于跳躍荷載的傅里葉級(jí)數(shù)模型中模型階數(shù)的選取研究較少,大多數(shù)認(rèn)為3階足以計(jì)算結(jié)構(gòu)的振動(dòng)響應(yīng),但往往卻忽略了結(jié)果的精確性。應(yīng)根據(jù)樓板的振動(dòng)響應(yīng)貢獻(xiàn)的大小來確定模型階數(shù)的選擇,做到簡(jiǎn)便化與精確化。

        參考文獻(xiàn):

        [1] 葉正強(qiáng),李愛群,丁幼亮,等.某大跨人行天橋的消能減振設(shè)計(jì)(一)[J].特種結(jié)構(gòu),2003,20(1): 68-70.

        [2] 于敬海,曹建東,李敬明. 跳躍荷載作用下大跨度鋼筋混凝土樓蓋舒適度分析[J]. 工程抗震與加固改造. 2013(06):64-67.

        [3] 朱萍. 跳躍荷載作用下的大跨樓蓋振動(dòng)響應(yīng)與控制研究[D]. 南京:東南大學(xué),2015.

        [4] 陳雋. 人致荷載與人致結(jié)構(gòu)振動(dòng)[M]. 北京:科學(xué)出版社,2016.

        [5] G.Pernica. Dynamic load factors for pedestrian movements and rhythmic exercises[J]. Canadian Acoustics. 2015,18(2):3-18.

        猜你喜歡
        有限元模擬
        CEFR—MOX燃料組件運(yùn)輸容器的實(shí)驗(yàn)分析與驗(yàn)證
        科技視界(2017年1期)2017-04-20 23:02:13
        SMA—疊層橡膠復(fù)合隔震支座耗能能力的影響因素研究
        地震研究(2017年1期)2017-04-10 18:26:56
        房屋結(jié)構(gòu)在地震作用下的有限元分析
        模具圓心角對(duì)AZ31鎂合金劇烈塑性變形的影響
        快速熱循環(huán)注塑模具及工藝關(guān)鍵技術(shù)研究
        深圳平安金融中心施工監(jiān)測(cè)與模擬研究
        加載速率對(duì)SAC系列焊點(diǎn)蠕變性能影響的研究
        科技視界(2016年12期)2016-05-25 00:25:20
        模擬技術(shù)在特形三通半管凈邊精沖成形中的應(yīng)用研究
        新型外包鋼—砼組合梁彎剪扭的非線性有限元模擬
        FRP加固鋼筋混凝土圓柱的有限元分析
        亚洲伊人伊成久久人综合| 亚洲经典三级| 久久精品无码一区二区三区蜜费| 人妻少妇喷水意淫诱惑| 亚洲成人精品久久久国产精品| 日本又色又爽又黄的a片18禁 | 国产极品美女高潮无套| 初女破初的视频| 亚洲精品中文字幕无乱码麻豆| 亚洲免费视频一区二区三区| 中文字幕av熟女中文av| 国产永久免费高清在线| 久久艹影院| 亚洲av套图一区二区| 国产流白浆视频在线观看| 欧美乱妇高清无乱码在线观看| 欧美色精品91av| 一区二区三区精彩视频在线观看| 国产精品第一二三区久久| 亚洲av无码专区首页| 无码午夜剧场| 在线看高清中文字幕一区| 精品亚洲成a人在线观看| av蓝导航精品导航| 91情侣在线精品国产免费| 中文字幕34一区二区| 日本艳妓bbw高潮一19| 精品人无码一区二区三区 | 天天看片视频免费观看| 国产视频不卡在线| 久久精品国产亚洲av久按摩| 色天使综合婷婷国产日韩av| 亚洲国产一区二区三区网| 日日噜噜夜夜久久密挑| 欧洲美熟女乱又伦av影片| 日日av拍夜夜添久久免费| 亚洲成AV人国产毛片| 亚洲综合自拍偷拍一区| 狠狠色噜噜狠狠狠狠7777米奇| 女的把腿张开男的猛戳出浆| 久久精品国产一区老色匹|