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

    隨機人群運動荷載作用下大跨度連廊的振動響應

    2019-06-10 00:00:00張瓊南娜娜朱前坤杜永峰
    土木建筑與環(huán)境工程 2019年2期
    關(guān)鍵詞:步速

    張瓊 南娜娜 朱前坤 杜永峰

    摘 要:基于社會力模型,研究了隨機人群運動荷載作用下大跨連廊結(jié)構(gòu)的振動舒適度問題。通過結(jié)合交通學領(lǐng)域、生物力學領(lǐng)域?qū)θ巳哼\動特性及人體運動參數(shù)的統(tǒng)計分析,對隨機人群運動進行仿真,得到不同密度工況下,任意行人在任意時刻的步速、落點位置等參數(shù),建立隨機人群荷載模型,計算不同工況密度下結(jié)構(gòu)的加速度峰值。研究表明,加速度峰值與行人步速關(guān)系很大,隨著行人密度增大,行人速度減小,結(jié)構(gòu)加速度響應則先增大后減小。根據(jù)得到的結(jié)構(gòu)加速度峰值與行人密度關(guān)系曲線,可以找出出現(xiàn)結(jié)構(gòu)振動舒適度問題的行人密度區(qū)間。

    關(guān)鍵詞:社會力模型;隨機人群荷載;步速;峰值加速度;振動舒適度

    中圖分類號:TU393.3

    文獻標志碼:A? 文章編號:2096-6717(2019)02-0099-07

    Abstract:Based on social force model which is combined statistics and analysis of crowd motion characters and body motion parameters in the field of transportation and biodynamic, vibration comfort of long-span corridor under stochastic crowd load is studied and stochastic crowd motion is simulated .The critical parameters such as pace and landing position are obtained in different densities, the stochastic crowd model is built and the peak acceleration is calculated .The study shows, peak acceleration is closely related to the pacing frequency, and pedestrian velocity becomes large along with increase of pedestrian density while the acceleration response of structure increased initially and then decreased. By the relationship curve of structure peak acceleration and pedestrian density the pedestrian density interval in which the structure suffering vibration serviceability problems can be found.

    Keywords:social force model; stochastic crowd load; walking speed; peak acceleration; vibration serviceability

    近年來,隨著新型輕質(zhì)材料的發(fā)展,以及人們對結(jié)構(gòu)新穎、美觀的追求,使得輕質(zhì)、低頻結(jié)構(gòu)的使用越來越廣泛。人群荷載作為連廊、人行橋等結(jié)構(gòu)的主要荷載,是引起結(jié)構(gòu)振動的主要原因,當結(jié)構(gòu)自振頻率較低時,易與行人步頻接近而發(fā)生共振,結(jié)構(gòu)振動響應超過人體對振動舒適度的要求,可能引起行人的恐慌甚至導致結(jié)構(gòu)的安全問題[1-3]。

    人群荷載模型是研究人致結(jié)構(gòu)振動的基礎,目前對單個行人對結(jié)構(gòu)的響應研究較多,但對人群荷載的研究較少,研究表明,人群對結(jié)構(gòu)的作用并不是單個行人的簡單疊加[4]。行人荷載與行人步頻有較大關(guān)系[5],人群平均步速與人群密度有關(guān),行人在不同人群密度下受到其他行人的影響不同,每個行人的移動方式及移動特征的表現(xiàn)也不盡相同,鑒于人群移動的復雜性以及隨機性,目前,人群對結(jié)構(gòu)的作用大都采用簡化的計算方法,主要包括固定多人步行力法[6]、移動多人步行力法、等效同步人數(shù)法[7]等,這些模型都在一定程度上提供了人群對結(jié)構(gòu)動力響應的計算方法。以上假設均沒有考慮行人個體的差異性以及個體間的隨機性。丁國等[8]通過引入步頻變異系數(shù)反映個體內(nèi)的隨機性,引入體重及均值步頻的變化反映個體間的隨機性,研究了行人荷載隨機性對樓蓋結(jié)構(gòu)振動響應的影響,結(jié)果表明,荷載隨機性對結(jié)構(gòu)響應影響顯著,在舒適度評價中不可忽略。楊娜等[9]根據(jù)生物力學以及交通領(lǐng)域內(nèi)對人體質(zhì)量、行人步頻、步速等參數(shù)的統(tǒng)計分析,得到各參數(shù)的分布形式,再利用蒙特卡洛法產(chǎn)生相應的偽隨機數(shù),并結(jié)合人群密度的影響作用生成考慮人群集度的隨機荷載模型。上述方法考慮了行人之間步態(tài)參數(shù)分布的隨機性以及個體在時間域內(nèi)的差異性,但是忽略了行人之間以及行人與環(huán)境的交互作用。因此,得到的行人落點、步速等并不符合行人的實際情況,計算結(jié)果與實際結(jié)構(gòu)響應可能會有較大偏差。

    由Helbing等[10]提出的社會力模型,將人的主觀愿望、行人與行人之間、行人與環(huán)境之間的相互作用用社會力的概念表示,成功模擬了“快即是慢”“瓶頸擺動”等行人現(xiàn)象。筆者結(jié)合交通學以及生物力學領(lǐng)域內(nèi)對行人步速、體重以及行人尺寸等參數(shù)的統(tǒng)計,考慮行人個體差異性以及個體之間的隨機性,利用社會力模型模擬人群的運動過程,得到任意時刻所有的行人落點、步速等;建立基于社會力模型的人群荷載隨機模型,并對不同行人密度下連廊結(jié)構(gòu)動力響應進行計算,得到結(jié)構(gòu)峰值加速度與行人密度的關(guān)系曲線,對隨機人群荷載作用下結(jié)構(gòu)振動舒適度進行評價。

    1 單人步行荷載模型

    人行走時兩條腿交替移動,步行荷載對結(jié)構(gòu)的作用具有周期性,其對連廊等結(jié)構(gòu)產(chǎn)生的激勵主要為豎向激勵,時域內(nèi)步行荷載廣泛采用傅里葉級數(shù)的形式[11]。

    式中:αi是i階諧波的動力荷載因子,定義為αi=Ai/G;φi為第i階諧波的相位角;G為行人重量;fp為行人步頻;n為模型中考慮的階數(shù);相位角φi=0??紤]前3階的影響,動載因子為

    2 隨機人群荷載模型的建立

    2.1 基于社會力模型隨機人群運動仿真

    社會力模型認為行人的運動在社會力的作用下發(fā)生,其中包括自驅(qū)力、行人之間的作用力、行人與周邊環(huán)境之間的作用力。筆者在此基礎上引入了減速避讓機制[12]以及自停止機制[13]。當行人的動態(tài)空間被其他行人占據(jù)時,當前行人會采取減速機制以避免與其他行人發(fā)生碰撞,而當周圍行人對當前行人斥力過大,使行人的速度小于0,則賦與當前行人速度為0。自驅(qū)力

    式中:Aw為行人與障礙物作用力強度;Bw為行人與障礙物作用范圍;ri-diw為行人半徑與行人到障礙物的法向距離差;niw為由障礙物指向行人的單位法向量; tiw為平行于障礙物的單位切向力;vi·tiw為行人速度在障礙物方向上的投影。

    利用MATLAB軟件對社會力模型進行仿真,仿真環(huán)境為一個長22.0 m、寬3.0 m的連廊。假定行人質(zhì)量在50~80 kg的范圍內(nèi)服從均值為65、方差為5的正態(tài)分布,行人期望速度與行人中老人、婦女的所占比重有關(guān),假定在0.7~1.4 m/s的范圍內(nèi)服從均值為1.29、方差為0.29的正態(tài)分布,半徑在0.2~0.3 m的范圍內(nèi)服從均勻分布。其他參數(shù)的取值參考文獻[13],如表1,仿真模擬如圖1所示。模擬結(jié)果表明,行人能夠有效地對其他行人以及障礙物進行避讓,將得到的密度速度曲線與其他人實測的或模擬結(jié)果進行對比,結(jié)果吻合度較高。

    2.2 行人速度密度關(guān)系

    利用社會力模型分別對2~90人范圍內(nèi)選取幾種不同密度工況進行仿真,行人位置分布會影響每個行人的受力,進而影響行人速度以及結(jié)構(gòu)的動力響應,考慮到行人的位置分布、期望速度等的隨機性以及仿真過程的波動性,每種密度工況下產(chǎn)生7組數(shù)據(jù),計算每種工況下速度的平均值,進而得到行人速度密度曲線,見圖2。將均值速度作為行人在某一密度下的移動速度,并與其他學者的結(jié)果進行對比,證明該社會力模型的可靠性。將得到的點進行擬合得到了密度速度關(guān)系式

    式中:k為行人密度,即單位面積的行人數(shù)。由圖2速度密度曲線可知,在行人密度小于0.5時行人速度的變化并不穩(wěn)定,原因是,當密度小于0.5時,行人處于自由運動的狀態(tài),速度具有較大的波動性,與實際情況符合;當密度大于0.5時,處于中密度狀態(tài),行人之間的相互作用開始變得明顯,速度隨著密度的增大呈穩(wěn)定下降趨勢。由文獻[14],速度與行人步頻的關(guān)系可表示為

    2.3 隨機人群荷載

    隨機人群荷載是單個行人荷載的組合。社會力模型中考慮到個體內(nèi)的隨機性,質(zhì)量在50~80 kg的范圍內(nèi)服從均值為65、方差為5的正態(tài)分布,期望速度在0.7~1.4 m/s的范圍內(nèi)服從均值為1.29、方差為0.29的正態(tài)分布。而行人真實速度則由社會力模型中任意時刻行人受力以及前一步速度共同決定。上述過程中,不僅在時間域內(nèi)考慮了行人個體內(nèi)以及個體間的隨機性,并且在空間域內(nèi)考慮了行人運動過程中周圍環(huán)境以及與周圍行人之間的相互作用。根據(jù)社會力模型,行人下一時刻的速度以及位置由當前行人受到的來自周圍行人以及障礙物的力共同決定,不同密度下行人受到其他行人影響不同。計算不同密度下任意行人任意時刻的位置分布、速度等,由式(10)可得到任意時刻行人的步頻,再由式(11)對隨機人群中所有荷載進行疊加則可得到人群隨機荷載,而人體自重引起的結(jié)構(gòu)動力響應相較于簡諧荷載而言可忽略不計[5]。

    式中:m為行人數(shù);n為在計算中考慮的總階數(shù);αi是i階諧波的動力荷載因子,定義為步行荷載傅里葉幅值譜峰值與人體體重之比;φi為第i階諧波的相位角;G為行人體重;fp為行人步頻。

    圖3與圖4中FP為所有行人步行荷載之和,G總為結(jié)構(gòu)上所有行人的體重之和。圖3為單個行人荷載時程曲線圖,從圖3可以明顯看出人行荷載具有周期性。圖4對比了行人同步行走荷載與隨機人群荷載時程曲線,其中人數(shù)為10人,同步行走中將10人排成3列,前2列均為4人,第3列為2人,行、列間距一致均為0.75 m。行走過程中行人步速均為1.3 m/s,步頻均為1.93 Hz。經(jīng)過對比可知,同步人群荷載時程與單個行人荷載時程類似,具有很明顯的周期性,而隨機人群荷載由于步頻并不統(tǒng)一,無明顯周期,行人行走過程中出現(xiàn)同步情況較少,因此,荷載峰值也比同步行走時的小。

    4 算例

    某商貿(mào)城的鋼結(jié)構(gòu)連廊跨度L=21.8 m,寬度為3.0 m,結(jié)構(gòu)的阻尼比ξ為0.01。該鋼結(jié)構(gòu)連廊截面豎向等效抗彎剛度為EI=3.268×109 N·m2,單位長度質(zhì)量m=1 063.5 kg,單位長度阻尼系數(shù)c為950.8 N/(m·s-1),連廊上基于社會力模型的人群隨機荷載分布如圖6所示?!督ㄖ巧w結(jié)構(gòu)振動舒適度技術(shù)規(guī)范(送審稿)》中規(guī)定,連廊的第一階豎向自振頻率不宜小于3 Hz,不封閉連廊豎向峰值加速度限制為0.5 m/s2。該結(jié)構(gòu)豎向基頻為4.72 Hz,滿足頻率限制的要求。當單個行人步頻為結(jié)構(gòu)基頻的1/3通過該連廊,由步頻和速度的關(guān)系式得行人的速度為0.838 m/s結(jié)構(gòu)加速度峰值如圖7所示為0.12 m/s2,小于加速度限制。

    《建筑樓蓋結(jié)構(gòu)振動舒適度技術(shù)規(guī)范(送審稿)》中指出,當行人密度為0.1人/m2是屬于完全自由行走狀態(tài)。當為0.5人/m2時行人較多,行人仍然能夠自由選擇行走,當行人密度等于1.0人/m2時,行人行走有不舒適感,人群擁擠。分別對行人密度為0.167、0.333、0.5、1.0人/m2共4種行人密度工況下的人群荷載引起的結(jié)構(gòu)響應進行計算。得到結(jié)果如圖8所示。從圖8可以看出當行人密度為0.167、0.333、0.5、1.0 人/m2時結(jié)構(gòu)加速度最大響應分別為0.066、0.143、0.255、0.535 m/s2,當行人密度為1.0人/m2時,結(jié)構(gòu)加速度響應大于0.5 m/s2,超過連廊振動舒適度限制。由于不同行人以及期望速度等對行人的運動特性影響較大,而在同一種行人密度工況下,行人位置分布隨機,期望速度取值服從正態(tài)分布??紤]到行人的隨機性以及仿真過程的波動性,同一行人密度工況下引起的結(jié)構(gòu)響應不同,針對10種行人密度工況各產(chǎn)生了8組不同行人模擬。表2為在各行人情況下8組不同行人分布下計算得到的結(jié)構(gòu)加速度峰值,其計算流程如圖5所示,取8組結(jié)構(gòu)加速度峰值的平均值作為對應行人密度下加速度響應的代表值,用以評估結(jié)構(gòu)的振動舒適度。為了得到更加明顯的變化規(guī)律,將得到的峰值加速度與人群密度進行函數(shù)關(guān)系擬和,該擬合只針對該結(jié)構(gòu)。擬合曲線如圖9,R2=0.922擬合度較高,擬合結(jié)果為

    為了驗證曲線的可靠性,分別計算0.25、0.583、0.75、0.917人/m2幾種行人密度工況下結(jié)構(gòu)加速度峰值,與公式擬合的值進行對比,結(jié)果見表3。通過對比實際峰值加速度與擬合公式計算得到的加速度峰值,發(fā)現(xiàn)結(jié)果較為吻合。由圖9密度加速度峰值關(guān)系曲線可知,當密度在1.15~1.30人/m2之間時,速度在0.89~0.84 m/s范圍內(nèi)變化,相應的由擬合公式計算得到的速度平均值為0.841 m/s,與之對應的步頻在1.57~1.63 Hz范圍內(nèi)變化,接近結(jié)構(gòu)基頻的1/3。故而結(jié)構(gòu)在1.15~1.30人/m2的密度范圍內(nèi)取最大值。當行人密度在1.0~1.4人/m2時,結(jié)構(gòu)加速度響應超出了舒適度限制,應該采取一定措施控制人群密度出現(xiàn)在這個范圍的概率,或者采取減振措施。

    5 結(jié)論

    1)以社會力模型為基礎,結(jié)合生物力學以及交通領(lǐng)域內(nèi)對于行人特性的研究,考慮了行人時間和空間內(nèi)的隨機性以及行人之間、行人與環(huán)境之間的交互性,建立了更符合實際的隨機人群荷載模型。

    2)人群隨機荷載引起的結(jié)構(gòu)振動與行人步速有較大關(guān)系,行人步速與行人密度有直接關(guān)系,可以通過行人密度與速度的關(guān)系曲線來預測不同密度下人群引起的結(jié)構(gòu)加速度響應,找出人群荷載引起結(jié)構(gòu)過量振動的行人密度范圍,采取措施進行減振控制。

    3)人群隨機荷載引起的結(jié)構(gòu)動力響應,并不會一直呈現(xiàn)增大的趨勢,當人群密度過大時,人群平均步頻遠小于結(jié)構(gòu)基頻的1/3,此時可將行人荷載看作是靜載,引起的動力響應很微弱。

    參考文獻:

    [1] 樊健生, 李泉, 李全旺,等. 考慮人-結(jié)構(gòu)相互作用的樓蓋振動控制研究[J]. 振動與沖擊, 2010, 29(11): 230-236.

    FAN J S, LI Q, LI Q W, et al. Floor vibration control considering human- structure interaction [J].Journal of Vibration and Shock, 2010, 29(11): 230-236. (in Chinese)

    [2] 簡方梁, 吳定俊, 李奇. 上海虹橋車站人行走廊人致振動分析[J]. 振動與沖擊, 2010, 29(8):136-140.

    JIAN F L, WU D J, LI Q. Human-induced vibration analysis of pedestrian suspension corridor in Shanghai Hongqiao Station [J]. Journal of Vibration and Shock, 2010, 29(8): 136-140. (in Chinese)

    [3] 張家華, 高承勇, 秦衛(wèi)紅,等. 大跨度樓面人群荷載作用響應分析[J]. 建筑結(jié)構(gòu)學報, 2011, 32(11): 177-182.

    ZHANG J H, GAO C Y, QIN W H, et al. Analysis of vibration responses of long-span floor under crowd loads [J]. Journal of Building Structures, 2011, 32(11): 177-182. (in Chinese)

    [4] VENUTI F, RACIC V, CORBETTA A. Pedestrian-structure interaction in the vertical direction: coupled oscillator-force model for vibration serviceability assessment [C] // IX International Conference on Structural Dynamics Eurodyn, Porto, Portugal, 2014.

    [5] 朱前坤, 李宏男, 杜永峰,等. 不同行走步速下人行橋振動舒適度定量化評估[J]. 工程力學, 2016, 33(10):97-104.

    ZHU? Q K, LI H N, DU Y F, et al. Quantitative evaluation of vibration serviceability of pedestrian bridge under different walking speed [J]. Engineering Mechanics, 2016, 33(10): 97-104. (in Chinese)

    [6] 甄偉, 盛平, 張力,等. 財富中心懸吊樓蓋結(jié)構(gòu)舒適度設計及檢測[J]. 建筑結(jié)構(gòu), 2011,41(1):52-54,34.

    ZHEN W, SHENG P, ZHANG L, et al. Vibration comfort design and experiment of the suspended floor of prosper center [J]. Building Structure, 2011,41(1):52-54,34. (in Chinese)

    [7] 潘寧, 杜文博, 倪志偉,等. 計算人群行走樓板響應的隨機振動方法[J]. 建筑科學, 2016, 32(5):39-43.

    PAN N, DU W B, NI Z W, et al. Random vibration method for calculating floor vibration response under crowd walking [J]. Building Science, 2016, 32(5): 39-43. (in Chinese)

    [8] 丁國, 陳雋. 行人荷載隨機性對樓蓋振動響應的影響研究[J]. 振動工程學報, 2016, 29(1): 123-131.

    DING G, CHEN J. Influences of walking load randomness on vibration responses of long-span floors [J]. Journal of Vibration Engineering, 2016, 29(1): 123-131. (in Chinese)

    [9] 楊娜, 邱雪. 人群集度荷載模型及人致振動模擬分析[J]. 北京交通大學學報, 2016, 40(3): 88-96.

    YANG N, QIU X. Crowd density load model and analysis of crowd-induced vibration [J]. Journal of Beijing Jiaotong University, 2016, 40(3): 88-96. (in Chinese)

    [10] HELBING D, BUZNA L, JOHANSSON A, et al. Self-organized pedestrian crowd dynamics: experiments, simulations, and design solutions [J]. Transportation Science, 2005, 39(1): 1-24.

    [11] PFEIL M, AMADOR N, PIMENTEL R, et al. Analytic-numerical model for walking person-footbridge structure interaction [C] // International Conference on Structural Dynamics, Eurodyn, 2014.

    [12] 李珊珊, 錢大琳, 王九州. 考慮行人減速避讓的改進社會力模型[J]. 吉林大學學報(工學版), 2012, 42(3): 623-628.

    LI S S, QIAN D L, WANG J Z. Improved social force model considering pedestrian deceleration to avoid collision [J]. Journal of Jilin University (Engineering and Technology Edition), 2012, 42(3): 623-628. (in Chinese)

    [13] PARISI D R, GILMAN M, MOLDOVAN H. A modification of the Social Force Model can reproduce experimental data of pedestrian flows in normal conditions [J]. Physica A Statistical Mechanics and Its Applications, 2009, 388(17): 3600-3608.

    [14] BRUNO L, VENUTI F. Crowd-structure interaction in footbridges: Modelling, application to a real case-study and sensitivity analyses [J]. Journal of Sound and Vibration, 2009, 323(1): 475-493.

    (編輯 胡英奎)

    猜你喜歡
    步速
    面向典型信號控制場景的行人過街特性探究
    試試4米步速測驗法
    握力和步速對急性加重期慢性阻塞性肺疾病病人BODE指數(shù)的影響研究
    全科護理(2023年20期)2023-07-21 01:06:30
    基于老年人行為特性的行人過街信號優(yōu)化
    杜克活動狀況指數(shù)與日常體力活動問卷的步速問題在中國與英國血液透析患者中的效度檢驗
    步速與住院高齡老年患者全因死亡的相關(guān)性
    老年人步速與跌倒風險的研究進展
    步速與住院老年高血壓患者全因死亡的相關(guān)性
    健康老年人舒適步速測評的Meta分析
    社區(qū)高齡老年人步速、握力與跌倒風險的相關(guān)性研究
    亚洲国产精品国产精品| 99精国产麻豆久久婷婷| 叶爱在线成人免费视频播放| 母亲3免费完整高清在线观看| 精品久久久久久电影网| 人人妻人人添人人爽欧美一区卜| 一级黄片播放器| 欧美av亚洲av综合av国产av| 国产免费福利视频在线观看| 又粗又硬又长又爽又黄的视频| 欧美亚洲 丝袜 人妻 在线| 日韩一卡2卡3卡4卡2021年| 50天的宝宝边吃奶边哭怎么回事| 亚洲五月色婷婷综合| 亚洲成人免费av在线播放| www.av在线官网国产| 十八禁人妻一区二区| 一边亲一边摸免费视频| 老司机亚洲免费影院| 丁香六月天网| 国产男女超爽视频在线观看| 99热网站在线观看| 美国免费a级毛片| 美女视频免费永久观看网站| 亚洲国产精品999| 久久免费观看电影| 在线观看国产h片| 下体分泌物呈黄色| 免费黄频网站在线观看国产| 亚洲欧美精品自产自拍| 日韩大片免费观看网站| 高清欧美精品videossex| 亚洲av在线观看美女高潮| 精品少妇黑人巨大在线播放| 99精品久久久久人妻精品| 亚洲国产精品一区三区| 国产精品国产三级国产专区5o| 两个人免费观看高清视频| 成年美女黄网站色视频大全免费| 一级片'在线观看视频| 亚洲欧美清纯卡通| 巨乳人妻的诱惑在线观看| 大陆偷拍与自拍| 黄色片一级片一级黄色片| 免费女性裸体啪啪无遮挡网站| 免费人妻精品一区二区三区视频| 色视频在线一区二区三区| 亚洲欧美精品综合一区二区三区| 亚洲视频免费观看视频| 久久免费观看电影| 50天的宝宝边吃奶边哭怎么回事| 日韩制服丝袜自拍偷拍| 成在线人永久免费视频| 大码成人一级视频| 曰老女人黄片| 80岁老熟妇乱子伦牲交| 欧美在线一区亚洲| 这个男人来自地球电影免费观看| 亚洲精品日韩在线中文字幕| 国产精品久久久久久精品古装| 久久人人97超碰香蕉20202| 午夜福利视频在线观看免费| 国产日韩欧美在线精品| www.999成人在线观看| 日本黄色日本黄色录像| 性色av乱码一区二区三区2| 亚洲一区二区三区欧美精品| 亚洲av男天堂| 亚洲 欧美一区二区三区| 咕卡用的链子| 欧美 日韩 精品 国产| 黑人巨大精品欧美一区二区蜜桃| 日韩电影二区| 尾随美女入室| 性高湖久久久久久久久免费观看| 国产精品.久久久| av有码第一页| 大香蕉久久网| 免费观看av网站的网址| 日韩av在线免费看完整版不卡| 日韩av在线免费看完整版不卡| 久久天躁狠狠躁夜夜2o2o | 亚洲伊人久久精品综合| √禁漫天堂资源中文www| 丰满迷人的少妇在线观看| 亚洲国产欧美网| 亚洲精品成人av观看孕妇| 久久亚洲精品不卡| 男女边摸边吃奶| 国产片内射在线| 一级黄片播放器| 亚洲成av片中文字幕在线观看| 日本午夜av视频| 午夜免费鲁丝| 人人妻,人人澡人人爽秒播 | 亚洲国产欧美网| 久久久国产欧美日韩av| 伊人亚洲综合成人网| 亚洲欧洲日产国产| 黄色 视频免费看| 成人免费观看视频高清| 操美女的视频在线观看| 97精品久久久久久久久久精品| 免费在线观看完整版高清| 成年人黄色毛片网站| 日本av手机在线免费观看| 午夜视频精品福利| 欧美xxⅹ黑人| 国产日韩一区二区三区精品不卡| 午夜视频精品福利| 大型av网站在线播放| 日本猛色少妇xxxxx猛交久久| 97精品久久久久久久久久精品| 国产91精品成人一区二区三区 | 欧美精品一区二区大全| 国产成人av激情在线播放| 飞空精品影院首页| 成人亚洲欧美一区二区av| 欧美黄色淫秽网站| 精品少妇一区二区三区视频日本电影| 欧美黄色片欧美黄色片| 99久久精品国产亚洲精品| 不卡av一区二区三区| 国产精品香港三级国产av潘金莲 | 国产一级毛片在线| 亚洲av日韩精品久久久久久密 | 老司机靠b影院| 国产在线免费精品| 亚洲综合色网址| 一级a爱视频在线免费观看| 国产91精品成人一区二区三区 | 日韩中文字幕视频在线看片| 久热爱精品视频在线9| 18禁国产床啪视频网站| 少妇人妻久久综合中文| 丝袜脚勾引网站| 高清视频免费观看一区二区| 精品国产乱码久久久久久小说| 99久久综合免费| www.自偷自拍.com| 亚洲av日韩精品久久久久久密 | 久久热在线av| 我的亚洲天堂| 9热在线视频观看99| 精品国产超薄肉色丝袜足j| 欧美 日韩 精品 国产| 18禁国产床啪视频网站| 久久久久久亚洲精品国产蜜桃av| 老司机在亚洲福利影院| 欧美日本中文国产一区发布| 免费久久久久久久精品成人欧美视频| 国产日韩欧美视频二区| 人体艺术视频欧美日本| 亚洲av电影在线观看一区二区三区| 在线观看www视频免费| 人人妻人人澡人人爽人人夜夜| 操美女的视频在线观看| 亚洲欧美一区二区三区国产| 久久精品人人爽人人爽视色| 免费av中文字幕在线| 亚洲精品国产区一区二| 国产精品av久久久久免费| 免费久久久久久久精品成人欧美视频| 中文字幕亚洲精品专区| 国产女主播在线喷水免费视频网站| 9色porny在线观看| 在线观看一区二区三区激情| 制服诱惑二区| 亚洲欧美中文字幕日韩二区| 大香蕉久久网| 国产高清videossex| 蜜桃国产av成人99| 性色av乱码一区二区三区2| 欧美大码av| 三上悠亚av全集在线观看| 狠狠婷婷综合久久久久久88av| 欧美另类一区| 亚洲精品一区蜜桃| 老鸭窝网址在线观看| 日日夜夜操网爽| 新久久久久国产一级毛片| 国产女主播在线喷水免费视频网站| 欧美日韩视频高清一区二区三区二| av线在线观看网站| 亚洲激情五月婷婷啪啪| 99九九在线精品视频| 国产av一区二区精品久久| 欧美人与性动交α欧美精品济南到| 中文字幕精品免费在线观看视频| 中文字幕人妻丝袜一区二区| 亚洲国产欧美在线一区| 韩国精品一区二区三区| 久久久亚洲精品成人影院| 亚洲精品av麻豆狂野| 国产真人三级小视频在线观看| 国产精品国产三级专区第一集| 色播在线永久视频| 热99久久久久精品小说推荐| 一本色道久久久久久精品综合| 精品少妇一区二区三区视频日本电影| 肉色欧美久久久久久久蜜桃| 看免费av毛片| 午夜免费男女啪啪视频观看| 色婷婷av一区二区三区视频| 国产精品久久久久久精品古装| 国产亚洲精品久久久久5区| 国产日韩一区二区三区精品不卡| 美女午夜性视频免费| 亚洲色图 男人天堂 中文字幕| 女人爽到高潮嗷嗷叫在线视频| 亚洲人成网站在线观看播放| 久久久久国产一级毛片高清牌| 欧美亚洲日本最大视频资源| 永久免费av网站大全| 午夜免费鲁丝| 亚洲视频免费观看视频| 国产亚洲欧美精品永久| 欧美成人精品欧美一级黄| 天堂中文最新版在线下载| 欧美久久黑人一区二区| 一区二区日韩欧美中文字幕| e午夜精品久久久久久久| 又紧又爽又黄一区二区| 久久人妻熟女aⅴ| 校园人妻丝袜中文字幕| 一区福利在线观看| 亚洲av男天堂| 免费看av在线观看网站| 久久久久久久久久久久大奶| 午夜日韩欧美国产| 久久鲁丝午夜福利片| 高清av免费在线| 亚洲 国产 在线| 久久精品成人免费网站| 欧美日韩国产mv在线观看视频| 黄片小视频在线播放| 久久久久久亚洲精品国产蜜桃av| 男人舔女人的私密视频| 制服诱惑二区| 午夜两性在线视频| 欧美精品人与动牲交sv欧美| 人妻人人澡人人爽人人| www.av在线官网国产| 欧美国产精品一级二级三级| 日本欧美视频一区| 少妇人妻 视频| 又紧又爽又黄一区二区| 天天操日日干夜夜撸| av在线播放精品| av欧美777| 人人妻人人澡人人看| 狂野欧美激情性bbbbbb| 男女之事视频高清在线观看 | h视频一区二区三区| 国产麻豆69| 一边摸一边做爽爽视频免费| 你懂的网址亚洲精品在线观看| xxxhd国产人妻xxx| 免费在线观看完整版高清| 婷婷丁香在线五月| 极品少妇高潮喷水抽搐| 亚洲av成人精品一二三区| 建设人人有责人人尽责人人享有的| 精品国产乱码久久久久久男人| 精品少妇久久久久久888优播| h视频一区二区三区| 色婷婷av一区二区三区视频| 婷婷成人精品国产| 永久免费av网站大全| netflix在线观看网站| 亚洲国产精品999| av网站在线播放免费| 久久久久久久国产电影| 婷婷色综合大香蕉| 侵犯人妻中文字幕一二三四区| 人体艺术视频欧美日本| 精品国产国语对白av| 一本一本久久a久久精品综合妖精| 成人三级做爰电影| 天天躁夜夜躁狠狠躁躁| 中文字幕人妻熟女乱码| 国产成人欧美| 亚洲五月婷婷丁香| 一级黄色大片毛片| 精品福利永久在线观看| 日本欧美国产在线视频| 亚洲九九香蕉| 国产免费一区二区三区四区乱码| 午夜福利免费观看在线| 久久精品久久精品一区二区三区| 亚洲欧洲国产日韩| 亚洲专区中文字幕在线| 观看av在线不卡| 一级a爱视频在线免费观看| 免费观看人在逋| 啦啦啦中文免费视频观看日本| 午夜福利乱码中文字幕| 日本猛色少妇xxxxx猛交久久| 在线观看人妻少妇| 日本av免费视频播放| 老司机在亚洲福利影院| 国产精品久久久久久人妻精品电影 | 亚洲国产看品久久| 91麻豆av在线| 国产亚洲精品第一综合不卡| 99香蕉大伊视频| 国产亚洲av高清不卡| 亚洲国产精品一区二区三区在线| 欧美人与性动交α欧美软件| 最近中文字幕2019免费版| 久久久久久久精品精品| 精品人妻熟女毛片av久久网站| 中文字幕最新亚洲高清| 成人午夜精彩视频在线观看| 亚洲国产欧美在线一区| 欧美97在线视频| 欧美黄色淫秽网站| 国产在线一区二区三区精| 一区在线观看完整版| 日韩中文字幕欧美一区二区 | 中文字幕制服av| 亚洲精品久久久久久婷婷小说| 老司机在亚洲福利影院| 亚洲国产精品一区二区三区在线| 色婷婷久久久亚洲欧美| 亚洲七黄色美女视频| 9191精品国产免费久久| 99精国产麻豆久久婷婷| 最近手机中文字幕大全| 国产精品av久久久久免费| 国产成人a∨麻豆精品| 久久亚洲国产成人精品v| www.精华液| 男人添女人高潮全过程视频| 一本—道久久a久久精品蜜桃钙片| 好男人电影高清在线观看| 亚洲美女黄色视频免费看| 亚洲熟女毛片儿| 1024视频免费在线观看| 亚洲欧美激情在线| 久久性视频一级片| 天堂俺去俺来也www色官网| a级毛片在线看网站| 黄色视频在线播放观看不卡| av网站免费在线观看视频| 欧美 亚洲 国产 日韩一| 亚洲人成77777在线视频| 中文字幕高清在线视频| 五月开心婷婷网| 啦啦啦在线免费观看视频4| 国产欧美日韩精品亚洲av| 国产高清国产精品国产三级| 久久天堂一区二区三区四区| 亚洲国产精品一区三区| 91精品三级在线观看| 天堂8中文在线网| 国产伦理片在线播放av一区| 亚洲情色 制服丝袜| 宅男免费午夜| 国产一区二区 视频在线| 亚洲人成网站在线观看播放| 丰满迷人的少妇在线观看| 亚洲成人国产一区在线观看 | 国产精品一区二区在线观看99| 久久99热这里只频精品6学生| 久久精品国产亚洲av涩爱| 亚洲av电影在线观看一区二区三区| 18禁国产床啪视频网站| 国产成人av激情在线播放| 亚洲,欧美,日韩| 久久午夜综合久久蜜桃| 亚洲精品美女久久久久99蜜臀 | 人体艺术视频欧美日本| 不卡av一区二区三区| 国产成人精品久久久久久| 日本欧美国产在线视频| 一级片免费观看大全| 午夜激情av网站| 亚洲av欧美aⅴ国产| 亚洲自偷自拍图片 自拍| 天堂中文最新版在线下载| 日本欧美国产在线视频| 国产福利在线免费观看视频| 亚洲欧美中文字幕日韩二区| 这个男人来自地球电影免费观看| 日韩一区二区三区影片| 波多野结衣一区麻豆| 老司机靠b影院| 别揉我奶头~嗯~啊~动态视频 | 亚洲av日韩精品久久久久久密 | 国产成人一区二区在线| 日本av手机在线免费观看| 看免费成人av毛片| 好男人视频免费观看在线| 99国产综合亚洲精品| 亚洲欧美一区二区三区久久| 亚洲七黄色美女视频| 下体分泌物呈黄色| 欧美亚洲 丝袜 人妻 在线| 啦啦啦啦在线视频资源| 国产国语露脸激情在线看| 男人舔女人的私密视频| 欧美黄色淫秽网站| h视频一区二区三区| 久久久久精品国产欧美久久久 | 国产精品.久久久| 亚洲国产欧美网| 丝袜美腿诱惑在线| 亚洲欧美中文字幕日韩二区| 汤姆久久久久久久影院中文字幕| 每晚都被弄得嗷嗷叫到高潮| 性少妇av在线| 国产xxxxx性猛交| 国产精品 欧美亚洲| 久久久欧美国产精品| 青青草视频在线视频观看| 国产一区二区三区综合在线观看| 男男h啪啪无遮挡| 色网站视频免费| 99国产精品一区二区三区| 美女脱内裤让男人舔精品视频| 热99国产精品久久久久久7| 交换朋友夫妻互换小说| 亚洲精品成人av观看孕妇| 欧美成狂野欧美在线观看| 黄色 视频免费看| 一区二区av电影网| 午夜福利视频在线观看免费| 夫妻午夜视频| 一本一本久久a久久精品综合妖精| 日日夜夜操网爽| 成在线人永久免费视频| 9热在线视频观看99| 亚洲人成电影观看| 飞空精品影院首页| 亚洲免费av在线视频| 黄色a级毛片大全视频| h视频一区二区三区| 免费av中文字幕在线| 午夜福利,免费看| 亚洲欧美日韩高清在线视频 | 国产一卡二卡三卡精品| 高清不卡的av网站| 成年美女黄网站色视频大全免费| 亚洲精品第二区| 日本色播在线视频| av在线app专区| 可以免费在线观看a视频的电影网站| 欧美97在线视频| 久久ye,这里只有精品| 国产成人一区二区在线| 亚洲成人手机| 天天操日日干夜夜撸| 日韩中文字幕欧美一区二区 | 久久狼人影院| 久久久久久久久免费视频了| 人人妻人人澡人人爽人人夜夜| 欧美性长视频在线观看| 欧美日韩成人在线一区二区| 国产成人欧美| 亚洲国产av新网站| 在线观看免费视频网站a站| 黑人欧美特级aaaaaa片| 国产精品一区二区在线观看99| 热99国产精品久久久久久7| 中国美女看黄片| 晚上一个人看的免费电影| 精品国产一区二区久久| 午夜影院在线不卡| 精品第一国产精品| videosex国产| 五月天丁香电影| 亚洲人成电影观看| 成人国产一区最新在线观看 | 成年动漫av网址| 青青草视频在线视频观看| 国产精品一区二区在线观看99| 老司机影院成人| 汤姆久久久久久久影院中文字幕| www.自偷自拍.com| 午夜久久久在线观看| 午夜福利,免费看| 好男人视频免费观看在线| 国产女主播在线喷水免费视频网站| 日韩中文字幕欧美一区二区 | 久久精品aⅴ一区二区三区四区| 青春草亚洲视频在线观看| 久久精品人人爽人人爽视色| 下体分泌物呈黄色| 中文字幕人妻熟女乱码| 国产真人三级小视频在线观看| 成人18禁高潮啪啪吃奶动态图| 欧美日本中文国产一区发布| 色播在线永久视频| 日韩一本色道免费dvd| 在线亚洲精品国产二区图片欧美| xxx大片免费视频| kizo精华| 男男h啪啪无遮挡| 纵有疾风起免费观看全集完整版| 亚洲国产欧美日韩在线播放| 欧美黄色淫秽网站| 欧美激情 高清一区二区三区| 99久久综合免费| 久久久久久久久久久久大奶| 午夜精品国产一区二区电影| 自拍欧美九色日韩亚洲蝌蚪91| 91麻豆精品激情在线观看国产 | 丝瓜视频免费看黄片| 国产激情久久老熟女| 日韩,欧美,国产一区二区三区| 无限看片的www在线观看| 十八禁人妻一区二区| 国产在线一区二区三区精| 欧美国产精品一级二级三级| 看十八女毛片水多多多| 国产精品一区二区精品视频观看| 日本一区二区免费在线视频| 国产真人三级小视频在线观看| 天堂8中文在线网| 少妇裸体淫交视频免费看高清 | 久久精品国产a三级三级三级| 欧美中文综合在线视频| 亚洲国产日韩一区二区| 老司机亚洲免费影院| 香蕉国产在线看| av在线app专区| 国产精品一区二区在线不卡| 欧美日韩综合久久久久久| 秋霞在线观看毛片| e午夜精品久久久久久久| 1024香蕉在线观看| 人妻人人澡人人爽人人| 丝袜人妻中文字幕| 精品亚洲成a人片在线观看| 男女高潮啪啪啪动态图| 夫妻性生交免费视频一级片| 国产精品免费视频内射| 菩萨蛮人人尽说江南好唐韦庄| 久久毛片免费看一区二区三区| 丝袜美足系列| 国产黄色免费在线视频| 乱人伦中国视频| 一级毛片电影观看| 欧美日韩综合久久久久久| 麻豆乱淫一区二区| 一级毛片 在线播放| 午夜福利,免费看| 免费在线观看视频国产中文字幕亚洲 | 亚洲熟女精品中文字幕| 男人添女人高潮全过程视频| 久久国产亚洲av麻豆专区| 日韩人妻精品一区2区三区| 国产深夜福利视频在线观看| 搡老乐熟女国产| 精品国产乱码久久久久久男人| 久久久精品区二区三区| 精品欧美一区二区三区在线| 看免费成人av毛片| 九色亚洲精品在线播放| 日本一区二区免费在线视频| 另类亚洲欧美激情| 天天躁日日躁夜夜躁夜夜| 久久性视频一级片| 日韩视频在线欧美| 91精品国产国语对白视频| 麻豆av在线久日| 涩涩av久久男人的天堂| 50天的宝宝边吃奶边哭怎么回事| 国产成人精品久久二区二区91| 天天躁夜夜躁狠狠躁躁| 伦理电影免费视频| 久久国产精品大桥未久av| 久久久精品区二区三区| 啦啦啦在线观看免费高清www| 十八禁网站网址无遮挡| 水蜜桃什么品种好| 亚洲欧美成人综合另类久久久| 午夜免费男女啪啪视频观看| 最近手机中文字幕大全| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久国产欧美日韩av| 啦啦啦中文免费视频观看日本| 母亲3免费完整高清在线观看| 亚洲av成人不卡在线观看播放网 | 精品国产超薄肉色丝袜足j| 91麻豆av在线| 视频区欧美日本亚洲| 一级黄色大片毛片| 亚洲成国产人片在线观看| 看免费成人av毛片| 精品少妇一区二区三区视频日本电影| 国产精品亚洲av一区麻豆| 黄色视频不卡| 亚洲专区国产一区二区| 一本—道久久a久久精品蜜桃钙片| 国产在线视频一区二区| 国产老妇伦熟女老妇高清| 日本wwww免费看| 高清视频免费观看一区二区| 亚洲欧美中文字幕日韩二区| 日本vs欧美在线观看视频| 亚洲精品美女久久久久99蜜臀 | 免费在线观看黄色视频的| 国产欧美日韩一区二区三 | 男人爽女人下面视频在线观看| 亚洲国产av新网站| 18禁国产床啪视频网站| 免费在线观看视频国产中文字幕亚洲 | 亚洲免费av在线视频| 亚洲精品美女久久av网站| 狂野欧美激情性bbbbbb| 爱豆传媒免费全集在线观看| 日韩,欧美,国产一区二区三区| 亚洲欧美日韩另类电影网站|