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

    水工混凝土結構日照輻射模擬方法

    2019-11-14 05:04:36頡志強呂興棟
    關鍵詞:日照射線溫度場

    頡志強,付 志,呂興棟

    (1. 長江水利委員會 長江科學院,湖北 武漢430010; 2. 水利部水工程安全與病害防治工程技術研究中心,湖北 武漢430010)

    0 引 言

    日照輻射對混凝土結構溫度、應力時空分布的影響不容忽視,對西藏高海拔地區(qū)的水工混凝土結構尤其如此。研究表明[1-3],日照輻射能大幅加劇結構溫度場時空分布不均勻性,而后者引起的變形與結構的約束共同導致了混凝土的開裂。因此,準確掌握強日照輻射對混凝土結構的溫度、應力的影響是該地區(qū)水工混凝土結構溫控防裂的前提。當前的水工混凝土溫控研究中,日照輻射影響分析絕大多數采用朱伯芳[1]引入的等效算法。該算法將日照輻射影響等效為氣溫增量,均勻施加于結構邊界面,能在平均意義上反映日照輻射引起的“溫升”效應,因此得到了廣泛的應用并取得了大量的研究成果[4-7]。雖然等效算法能夠近似反映日照輻射對結構溫度場時空分布的影響,但仍有明顯不足,主要表現在其無法準確反映周邊環(huán)境(山體或其他建筑物)、結構的方位、朝向等因素所導致的結構受照不均勻性。也就是說,等效算法在實質上弱化了結構受照的不均勻性。為了進一步提高分析精度,陳拯[8]將ASHRAE日照輻射模型及遮蔽算法引入到拱壩運行期溫度場計算中,取得了良好的效果;蘇衛(wèi)強等[9]對主要的日照輻射模型的適應性進行了探討,為仿真過程中日照輻射模型的選取提供了有意義的建議。通過對現有研究成果的總結,筆者認為水工混凝土結構日照輻射模擬主要涉及3個方面的內容:①日照參數計算,主要是對不同時間、經緯度、海拔高程的日照輻射強度、日照時間、入射角進行計算;②結構受照過程模擬,主要是指對特定地區(qū)、特定形式、方位的結構在日照過程中,表面各部位受日照強度、入射角、遮蔽情況進行模擬;③日照輻射影響下溫度應力仿真,主要將日照輻射影響引入混凝土結構溫度場、應力場仿真計算中,評價日照輻射對結構溫度、應力時空分布特性的影響。為此,筆者系統(tǒng)地梳理了相關計算方法,針對傳統(tǒng)遮蔽判斷算法繁瑣,計算效率低的問題,引入相交快速算法提高了邊界面遮蔽判斷效率,簡化了算法流程。開發(fā)了日照輻射精細化仿真分析程序。通過高海拔地區(qū)某電站開展的立桿實驗及墩墻溫度應力分析,驗證了算法的合理性和準確性。

    1 天文參數計算

    太陽光穿過大氣層后,除了散射和吸收的部分外,最終達到地面的輻射為直接輻射,大氣層散射或再發(fā)射的輻射為散射輻射,從其它表面反射的輻射為反射輻射。對于混凝土結構溫度場而言,太陽輻射的影響主要與結構受照表面的輻射強度和太陽入射角有關。

    1.1 輻射強度

    到達地球的太陽輻射以短波輻射為主,對于地表任意受照平面而言,太陽輻射主要包括:直接輻射、天空散射、地面反射3個部分。

    1)直接輻射:是指太陽光經過大氣直接到達地球表面的輻射強度,對于地面太陽入射角為φ的結構表面,直接輻射強度可按式(1)計算:

    IDφ=IDcosφ

    (1)

    式中:ID為與太陽光垂直平面的輻射強度,根據KEHLBECK模型:

    ID=I00.9tukam

    (2)

    式中:I0為太陽常數;tu為林克氏混濁度系數,隨大氣變化而變化;ka為不同海拔高度的相對大氣壓,m為大氣光學質量,按式(3)~式(5)計算或取值。

    (3)

    (4)

    (5)

    式中:Atu、Btu為與計算區(qū)域相關的參數對處于山區(qū)的水利工程而言,一般Atu=2.2,Btu=0.5。

    表1 林克氏混濁度系數參數取值[10]Table 1 Parameters of Rink’s turbidity coefficient

    2)天空散射:是指達到地球表面的被大氣層散射的太陽輻射,天空散射與受照表面的方位角無關,只與表面傾角有關,任意表面天空散射強度為[11]

    (6)

    式中:βn為受照平面法向與地平面的夾角;IdH為水平表面所受到的天空散射強度,按式(7)計算:

    IdH=(0.271I0-0.294ID)sinβs

    (7)

    3)地面反射:是指直接輻射和天空散射投射到地面以后被地表反射的太陽輻射,地面反射強度根據式(8)計算[10]:

    (8)

    式中:re為地面短波反射率,一般地面取值0.2,積雪地面取值0.7。

    由式(1)、式(6)、式(8)可知,對于混凝土結構任意受照表面而言,總的太陽輻射為

    (9)

    1.2 幾何參數

    對于地面上任意表面在空間的傾向可以由表面外法線n的方位角αn以及n與地平面的夾角βn兩者完全確定。而相對于地面上的觀察者,太陽在天空半球內的位置可以由太陽方位角αs和高度角βs兩者所完全確定。

    建立圖1所示的空間直角坐標系OXYZ。坐標原點置于受照面上,OXY位于地平面內,X軸指向正南(S),Y軸指向正東(E),Z軸垂直地平面指向天頂。設n和s分別表示表面外法線n方向和太陽射線方向(從原點指向太陽)的單位向量。由圖1,受照面的入射角φ即為向量n與向量s的夾角。根據幾何關系可知,兩向量可分別表示為

    n=icosβncosαn+jcosβnsinαn+ksinβn

    (10)

    s=icosβscosαs+jcosβssinαs+ksinβs

    (11)

    因此,n和s之間的夾角即太陽入射角為

    cosφ=n·s=cosβncosβscos(αn-αs)+ sinβnsinβs

    (12)

    圖1 向量關系Fig. 1 Relationship of vectors

    采用式(12)計算入射角時,還需確定如下參數:

    1)太陽高度角βs,是指自地面觀察者至太陽的連線與觀察者所在的地平面的夾角,按式(13)計算:

    sinβs=cosφcosδcosτ+ sinφsinδ

    (13)

    2)太陽方位角αs:太陽方位角αs為自地面觀察者至太陽的連線在地平面上的投影與正南方向的夾角,一般規(guī)定太陽方位角以由南轉向東為正,向西為負,即上午為正,下午為負,任意時刻的方位角有:

    (14)

    (15)

    3)太陽傾角δ:亦稱太陽赤緯,是指由地心指向日心的連線與地球赤道平面的夾角。太陽傾角在春分和秋分時刻等于0,夏至和冬至時刻出現極值,分別為+23.45°和-23.45°,其在一年中的變化可按經驗公式(16)計算。

    (16)

    式中:N為日序數。

    實際計算中,首先根據太陽日計算太陽傾角δ,再根據δ、計算時間τ及受照面地理緯度φ,依次計算太陽高度角βs和方位角αs,通過簡單推導有[10]:

    cosφ=(sinφsinβn- cosφcosαncosβn)sinδ+

    (cosφcosβn- sinφcosαncosβn)cosδcosτ+

    sinαncosβncosδsinτ

    (17)

    根據式(17)可知,在已知結構方位、角度、地理位置(經度、緯度)、日期、時間情況下,可以方便地計算出結構任意表面(非遮蔽)日照輻射入射角。

    1.3 受照時段

    在日照輻射精細化仿真分析中,需準確考慮日出日落時間,根據公式(13)令高度角βs=0,則特定經緯度地區(qū)日出日落時間為

    (18)

    基于此,可求得日出日落真太陽時tt和ts分別為

    (19)

    2 結構受照過程模擬算法

    水利工程一般地處高山峽谷地區(qū),在受照過程中不可避免地受到周邊山體或自身不同部位的遮蔽。遮蔽效應進一步加劇了日照輻射的空間不均勻性。因此,在結構受照過程模擬中,必須對有限元模型進行實時遮蔽判斷分析。

    2.1 邊界面遮蔽分析方法

    天文參數計算完成后,即可對于有限元模型進行遮蔽判斷分析,算法思路如圖2,對于某待分析面,從該面沿入射光線的反向引出一條射線,判斷該射線是否與模型其他面相交,若相交則該面被遮蔽,若不相交則該面未被遮蔽。例如,圖2中i號面引出的一條與入射光線平行的射線,該射線與岸坡相交,說明i號面被遮蔽,而j面未被遮蔽因此可以被光線直接照射。需要注意的是,實際分析中亦可對邊界節(jié)點進行逐個分析,邊界節(jié)點所在引出的射線方向為其所在所有面的入射向量的均值。

    圖2 遮蔽效應Fig. 2 Shadow effect

    結構遮蔽分析核心是線、面相交判斷,計算量巨大,以圖2中的小拱壩為例,邊界面總數為6 086個,溫度場計算中,每個時步都需要分別對6 086個臨空面引出的射線與其他邊界面(共6 085個)進行至少3.6億次以上的相交判斷,才能夠準確分析表面受照情況。顯然,在原本耗時巨大的仿真分析中,加入巨量的相交分析,既不經濟也不合理。為此,文獻[8]引入了柵格加速算法減少判斷次數。主要是建立如圖3的虛擬柵格,特定射線只需與所經柵格內的臨空面進行相交判斷。這種處理方式能夠有效減少射線與臨空面的相交判斷次數,但仍然需要大量的相交計算。

    圖3 柵格分區(qū)Fig. 3 Grid partition

    2.2 射線空間面相交快速判斷

    三維有限元分析絕大多數為六面體或四面體單元,因此應用中可以統(tǒng)一簡化為射線與空間三角形的相交判斷。常規(guī)做法是直接判斷射線與三角形面相交,再判斷交點是否在三角形內,則需要二次判斷,計算過程繁瑣,且效率不高?,F引入一種新的射線與空間三角形相交快速判斷方法,進一步提高遮蔽判斷效率。

    設任意射線的方程為

    (20)

    式中:(x0,y0,z0)為射線的起點,向量(m,n,p)為射線的方向向量。令:

    則,射線方程可表示為

    X=Dt+O

    (21)

    對于空間三角形,其內部任意點可表示為

    (22)

    式中:(xi,yi,zi)(i=1,3)為空間三角形的3個角點,令:

    則,空間三角形可表示為

    X=(1-u-v)V1+uV2+vV3

    (23)

    若射線與空間三角形相交,則有:

    Dt+O=(1-u-v)V1+uV2+vV3

    (24)

    即:

    Dt+u(V1-V2)+v(V1-V3)=V1-O

    (25)

    令,P0=V1-O,P1=V1-V2,P2=V1-V3,則有:

    Dt+uP1+vP2=P0

    (26)

    顯然,式(26)為簡單的三階線性方程組,易知式(26)的解為

    (27)

    利用式(27),可以快速地判斷射線是否與三角形相交。

    基于上述算法,筆者采用Fortran開發(fā)了大體積混凝土結構日照輻射過程精細化仿真模擬計算程序。采用PC機(AMD Athlon(tm) II X64 2.8GHz)對圖2所示拱壩進行遮蔽分析,單個時間步天文參數計算及遮蔽判斷總耗時僅0.01 s左右。

    2.3 試驗驗證

    為進一步驗證上述算法和程序的正確性,依托高海拔地區(qū)在建工程營地現場,開展了現場立桿試驗。試驗裝置如圖4,立桿高度2.0 m,立桿底部經緯度為北緯29°15′28″,東經92°25′12″,海拔高程3 377 m,立桿高度為2.0 m,試驗時間為2017年11月3日。根據立桿尺寸,建立如圖5的有限元網格,該網格節(jié)點總數為5 242,單元總數2 510,x軸正方向為正南方向。本試驗主要通過模擬立桿陰影的變化驗證天文參數、遮蔽算法和程序的準確性。

    圖4 陰影試驗裝置Fig. 4 Shadow test device

    圖6 試驗與計算結果比較Fig. 6 Comparison of test and calculation results

    根據計算,試驗地點當天的理論日出時間為06:03,日落時間為17:23,由于兩側山體遮蔽,測點在上午09:00左右才受到日照作用,下午16:00后太陽被山體遮蔽,也無法測得陰影數據。試驗測得的立桿陰影與計算得到的陰影對比見圖6,其中紅色短線為實測陰影,黑色區(qū)域為計算得到的陰影。對比可知,除了上午09:00計算陰影與實測陰影在長度上略有差異外,其他各整點計算結果與實測結果無論陰影角度還是陰影長度均符合良好。09:00誤差主要是由陰影測量誤差(上午陰影不夠清晰)所致,總體而言,所述算法及開發(fā)的日照輻射模擬程序能夠有效地對日照輻射過程及受照情況進行精細化仿真模擬。

    3 考慮日照輻射影響的混凝土結構溫度場算法

    通過天文參數計算、受照過程模擬分析,確定結構的受照情況(即表面日照輻射熱流密度等)后,可進一步在溫度場計算中考慮日照輻射的影響。大體積混凝土溫度場、應力場仿真相關理論在文獻[2]中有詳細的介紹,不再贅述,此處僅對溫度場計算中日照輻射的精細化模擬方法進行介紹。根據熱傳導理論,對于結構任意邊界面,熱流密度可根據式(28)計算:

    (28)

    式中:λ為混凝土導熱系數;T為混凝土溫度;qs、qc、qr分別為通過日照短波輻射、表面對流換熱、物體間長波輻射進入混凝土內的熱流密度,對于水工混凝土結構而言,物體間的長波輻射對混凝土溫度場影響很小,一般不予考慮。而對流換熱引起的熱流量為

    qc=β(T-Ta)

    (29)

    考慮表面日照輻射吸收率為αs,則由于日照輻射引起的表面熱流量為

    qs=-αs×(IDφ+Idβ+Irβ)

    (30)

    進一步,在推導有限元計算格式時,對于受照邊界面,利用式(31)即可。

    (31)

    基于此,筆者將開發(fā)的日照輻射模塊引入長江科學院現有大體積混凝土溫控仿真計算程序中,開發(fā)了能夠充分考慮日照輻射影響的大體積混凝土溫度場、應力場仿真計算程序。

    4 高海拔地區(qū)混凝土結構日照響應分析

    4.1 工程概況

    以前述工程中某南北走向墩墻為研究對象,計算2017年7月20日全天墩墻溫度場、應力場變化過程。該墻高6.0 m、墻體厚度為3.0 m、長10.0 m,墻體及地基有限元模型如圖7,模型單元總數為11 590個、節(jié)點總數為13 560個?;炷翉椥阅A繛?7.5 GPa,線膨脹系數為8.0×10-6/℃,基巖彈性模量為30 GPa。

    圖7 墩墻有限元模型Fig. 7 Finite element model of pier wall

    為便于比較,分別計算考慮和不考慮日照輻射影響情況下,墻體一天內的溫度、應力時空分布。在墻體中截面上,選取T1~T5 5個特征點,其中T1為墩墻東側表面點、T2距東側表面0.4 m,T3為西側表面點,T4距西側表面0.4 m,T5為墩墻中心點,見圖8。

    圖8 特征點Fig. 8 Typical points

    4.2 結果分析

    圖9~圖11為考慮和不考慮日照輻射影響下,墩墻中部特征點T1~T5溫度、應力歷程。

    圖9 陽面特征點、溫度應力Fig. 9 Typical points and temperature stresses in sunward side

    溫度方面,根據計算,考慮日照輻射以后,上午10點東側表面點T1溫度由于日照輻射引起的溫差達到了16.6 ℃,下午17點左右,墩墻西側表面點T3由于日照輻射影響溫度同樣達到了16.7 ℃。考慮日照影響后,兩側表面附近(T2、T4)均明顯高于無日照工況,T2點最大溫差2.0 ℃左右,但與表面不同,溫差最大值出現的時間在20點左右,有明顯的延遲??紤]日照輻射影響后,由于墻體較厚(3.0 m)墻體中心點溫度受日照輻射影響不大。

    應力方面,在墻體兩側(即T1、T3),由于日照輻射導致的溫度升高,表面產生了一定程度的壓應力。由于日照輻射導致內外溫差增大,表面附近(即T2、T4)拉應力大幅增加,以東側T2點為例,不考慮日照影響時T2點第一主應力最大值為0.25 MPa左右,考慮日照輻射影響后,第一主應力達到了0.9 MPa左右。說明對于實際工程而言,日照輻射能夠大幅增加表面混凝土開裂風險。同樣,對于墩墻中心點T5,由于日照輻射影響,第一主應力最大值由原來的0.2 MPa增大至0.9 MPa以上。

    計算成果表明,日照輻射明顯改變了結構溫度場時空分布特性,相應地,對結構溫度應力時空分布也有顯著影響。日照輻射的存在大幅加速了結構表面的開裂風險,對結構工作性態(tài)影響明顯,因此進行日照輻射精細化模擬是必要的。

    圖10 陰面特征點溫度、應力歷程Fig. 10 Temperature and stress history of typical points in shaded side

    圖11 T5特征點溫度、應力歷程Fig. 11 Temperature and stress history of T5 typical points

    5 結 論

    1)準確分析了日照輻射對結構的影響,對于推進仿真分析精細化和掌握結構真實工作性態(tài)至關重要,也是當前智能建造的迫切需求。筆者系統(tǒng)地梳理了水工混凝土結構日照輻射仿真模擬方法,給出了相應的計算公式,并開發(fā)了程序,實現了對水工混凝土結構日照響應的精細化模擬。

    2)針對傳統(tǒng)遮蔽判斷算法繁瑣、效率低下的問題,引入相交判斷快速算法,避免了二次判斷,有效提高了遮蔽分析效率,滿足了大型工程實時日照輻射響應分析對計算效率的要求。

    3)算例表明,對高海拔地區(qū)水工混凝土結構而言,日照輻射會大幅加劇溫度場、應力場時空分布的不均勻性,增加混凝土結構開裂風險。因此在高海拔地區(qū)水工混凝土結構溫控分析中,對日照輻射過程進行精細化仿真模擬是必要的。

    猜你喜歡
    日照射線溫度場
    心之向往,便是日照
    金橋(2021年9期)2021-11-02 05:54:52
    滄海桑田話日照
    金橋(2021年9期)2021-11-02 05:54:30
    鋁合金加筋板焊接溫度場和殘余應力數值模擬
    “直線、射線、線段”檢測題
    『直線、射線、線段』檢測題
    基于紋影法的溫度場分布測量方法
    測控技術(2018年4期)2018-11-25 09:47:10
    MJS工法與凍結法結合加固區(qū)溫度場研究
    建筑科技(2018年6期)2018-08-30 03:41:08
    在日照(節(jié)選)
    赤石脂X-射線衍射指紋圖譜
    中成藥(2017年3期)2017-05-17 06:09:16
    日照行
    小布老虎(2016年18期)2016-12-01 05:47:41
    麻豆精品久久久久久蜜桃| 97超级碰碰碰精品色视频在线观看| 国语自产精品视频在线第100页| 免费观看精品视频网站| 国产成人91sexporn| 亚洲精品一卡2卡三卡4卡5卡| 丝袜美腿在线中文| 欧美三级亚洲精品| 啦啦啦啦在线视频资源| 亚洲真实伦在线观看| 在线观看免费视频日本深夜| 最近2019中文字幕mv第一页| 18禁在线播放成人免费| 在线观看66精品国产| 真人做人爱边吃奶动态| 免费大片18禁| 看非洲黑人一级黄片| 婷婷精品国产亚洲av| 成人一区二区视频在线观看| 亚洲一区高清亚洲精品| 18禁在线无遮挡免费观看视频 | 91久久精品电影网| 久久久欧美国产精品| 十八禁国产超污无遮挡网站| 久久99热6这里只有精品| 别揉我奶头 嗯啊视频| 最好的美女福利视频网| 欧美人与善性xxx| 免费黄网站久久成人精品| 日韩欧美精品免费久久| av中文乱码字幕在线| 国产极品精品免费视频能看的| 久久久久国产精品人妻aⅴ院| 精品人妻一区二区三区麻豆 | 久久九九热精品免费| 精品午夜福利在线看| 一个人看的www免费观看视频| 亚洲欧美精品综合久久99| 最近在线观看免费完整版| 我要看日韩黄色一级片| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲精品一区av在线观看| 日韩精品中文字幕看吧| 免费一级毛片在线播放高清视频| 亚洲欧美精品综合久久99| 久久午夜亚洲精品久久| 校园人妻丝袜中文字幕| 亚洲av美国av| 精华霜和精华液先用哪个| 日韩欧美 国产精品| 中文资源天堂在线| 高清毛片免费看| 精品欧美国产一区二区三| 深爱激情五月婷婷| 丰满人妻一区二区三区视频av| 91久久精品电影网| 熟女电影av网| 亚洲人成网站在线播放欧美日韩| 免费看光身美女| 欧美一级a爱片免费观看看| 国产视频内射| 国产亚洲av嫩草精品影院| 我要搜黄色片| 亚洲18禁久久av| 精品国产三级普通话版| 亚洲精品在线观看二区| 成人特级黄色片久久久久久久| 内射极品少妇av片p| 十八禁国产超污无遮挡网站| 国产精品一二三区在线看| 日韩高清综合在线| 成人漫画全彩无遮挡| 亚洲精品色激情综合| 色噜噜av男人的天堂激情| 欧美+亚洲+日韩+国产| 看十八女毛片水多多多| 久久久久久久久久成人| 中文字幕人妻熟人妻熟丝袜美| 亚洲欧美精品综合久久99| 亚洲国产色片| 老司机影院成人| 日日摸夜夜添夜夜爱| 黑人高潮一二区| 国产蜜桃级精品一区二区三区| 99国产极品粉嫩在线观看| 成人永久免费在线观看视频| 国内精品久久久久精免费| 亚洲av电影不卡..在线观看| 在线看三级毛片| 久久精品国产亚洲网站| 国产人妻一区二区三区在| 中国美女看黄片| 偷拍熟女少妇极品色| 国产精品精品国产色婷婷| avwww免费| 在线观看一区二区三区| 日韩精品青青久久久久久| 国产视频一区二区在线看| 在线天堂最新版资源| 三级经典国产精品| 熟妇人妻久久中文字幕3abv| 国内精品久久久久精免费| 亚洲精品456在线播放app| 国产片特级美女逼逼视频| 国产精华一区二区三区| 少妇猛男粗大的猛烈进出视频 | 少妇熟女欧美另类| 欧美日韩精品成人综合77777| 老司机午夜福利在线观看视频| 久久韩国三级中文字幕| 舔av片在线| 亚洲精品成人久久久久久| 日本成人三级电影网站| 亚洲精品乱码久久久v下载方式| 久久精品夜色国产| 性插视频无遮挡在线免费观看| 日韩成人av中文字幕在线观看 | 国产69精品久久久久777片| 中文字幕人妻熟人妻熟丝袜美| 岛国在线免费视频观看| 亚洲精品粉嫩美女一区| 一进一出好大好爽视频| 久久人人精品亚洲av| 此物有八面人人有两片| 12—13女人毛片做爰片一| 变态另类成人亚洲欧美熟女| 亚洲av五月六月丁香网| 国产免费男女视频| 成人一区二区视频在线观看| 精品久久久久久久久久免费视频| 精品久久久噜噜| 亚洲天堂国产精品一区在线| 99热精品在线国产| 可以在线观看毛片的网站| 国产乱人偷精品视频| 午夜激情福利司机影院| 男人和女人高潮做爰伦理| 亚洲国产精品成人久久小说 | 国产一区二区在线av高清观看| 国产成人福利小说| 少妇的逼好多水| 黑人高潮一二区| 国产精品一及| 色综合亚洲欧美另类图片| 男人舔女人下体高潮全视频| 菩萨蛮人人尽说江南好唐韦庄 | 久久久精品94久久精品| 日韩强制内射视频| 舔av片在线| av视频在线观看入口| 国产精品免费一区二区三区在线| 国产精品,欧美在线| 久久久成人免费电影| 最新中文字幕久久久久| 色在线成人网| 亚洲av中文av极速乱| 婷婷精品国产亚洲av在线| 神马国产精品三级电影在线观看| 国产 一区精品| 国产精品久久久久久亚洲av鲁大| eeuss影院久久| 精品久久国产蜜桃| 亚洲最大成人av| 一级毛片电影观看 | 一夜夜www| 亚洲五月天丁香| 欧美xxxx黑人xx丫x性爽| 嫩草影院入口| 最后的刺客免费高清国语| 91久久精品国产一区二区成人| 一区二区三区免费毛片| 观看美女的网站| 欧美zozozo另类| 高清毛片免费观看视频网站| 国产成人aa在线观看| 一级av片app| 日本在线视频免费播放| 亚洲国产精品国产精品| 观看免费一级毛片| 乱码一卡2卡4卡精品| 看片在线看免费视频| 特级一级黄色大片| 日韩精品中文字幕看吧| 色综合色国产| 欧美最新免费一区二区三区| 大型黄色视频在线免费观看| 亚洲乱码一区二区免费版| 久久婷婷人人爽人人干人人爱| 久久久久久伊人网av| 久久精品影院6| 亚洲最大成人av| 国产成人a∨麻豆精品| 日本欧美国产在线视频| 成人性生交大片免费视频hd| 亚洲久久久久久中文字幕| 老熟妇仑乱视频hdxx| 欧美日韩精品成人综合77777| or卡值多少钱| 成人鲁丝片一二三区免费| 中文字幕精品亚洲无线码一区| 亚洲成a人片在线一区二区| 中出人妻视频一区二区| 国内精品久久久久精免费| 99久久精品一区二区三区| 内射极品少妇av片p| .国产精品久久| av在线蜜桃| 一进一出抽搐gif免费好疼| 欧美日本亚洲视频在线播放| 午夜a级毛片| 国产成年人精品一区二区| 丰满乱子伦码专区| 丝袜美腿在线中文| 狂野欧美白嫩少妇大欣赏| 亚洲第一电影网av| 久久久欧美国产精品| 国产精品99久久久久久久久| 欧美日韩综合久久久久久| 免费观看在线日韩| 成人三级黄色视频| 女生性感内裤真人,穿戴方法视频| 国产黄色小视频在线观看| 国产黄色小视频在线观看| 六月丁香七月| 国产精华一区二区三区| 国产日本99.免费观看| 国产精品久久电影中文字幕| 99久久精品一区二区三区| 国内精品一区二区在线观看| 国产白丝娇喘喷水9色精品| 九九在线视频观看精品| 蜜臀久久99精品久久宅男| 男女做爰动态图高潮gif福利片| 99热只有精品国产| 久久久精品94久久精品| 岛国在线免费视频观看| 乱系列少妇在线播放| 校园春色视频在线观看| a级一级毛片免费在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 九九在线视频观看精品| 少妇人妻一区二区三区视频| 精品午夜福利在线看| 日韩强制内射视频| 久久6这里有精品| 午夜精品在线福利| 国产精品久久久久久久久免| 国产精品久久电影中文字幕| 大又大粗又爽又黄少妇毛片口| 波多野结衣巨乳人妻| 亚洲国产精品合色在线| 国内精品美女久久久久久| 97人妻精品一区二区三区麻豆| 天堂√8在线中文| 日本爱情动作片www.在线观看 | 男人狂女人下面高潮的视频| 观看免费一级毛片| 国产成人freesex在线 | 亚洲性夜色夜夜综合| 久久午夜亚洲精品久久| 一级黄片播放器| 级片在线观看| 中国美白少妇内射xxxbb| 在线观看66精品国产| 最好的美女福利视频网| 色吧在线观看| 欧美激情在线99| 热99re8久久精品国产| 内地一区二区视频在线| 免费一级毛片在线播放高清视频| 免费无遮挡裸体视频| 久久久久久九九精品二区国产| 最新中文字幕久久久久| 亚洲性夜色夜夜综合| 国产精品一区www在线观看| 欧洲精品卡2卡3卡4卡5卡区| 中文字幕久久专区| 看十八女毛片水多多多| 久久人人精品亚洲av| 五月玫瑰六月丁香| 日韩成人av中文字幕在线观看 | 又黄又爽又免费观看的视频| 欧美一区二区精品小视频在线| 久久精品国产自在天天线| 久久久精品欧美日韩精品| 日韩高清综合在线| 99久久成人亚洲精品观看| 亚洲,欧美,日韩| 12—13女人毛片做爰片一| 搡女人真爽免费视频火全软件 | 最近手机中文字幕大全| 久久久久性生活片| 久久久a久久爽久久v久久| 国产成人freesex在线 | av卡一久久| av女优亚洲男人天堂| 国产高清三级在线| 日韩国内少妇激情av| 12—13女人毛片做爰片一| 亚洲av.av天堂| 一级毛片久久久久久久久女| 亚洲成人久久爱视频| 午夜精品一区二区三区免费看| 内射极品少妇av片p| 99国产极品粉嫩在线观看| 国产精品一区二区三区四区久久| 两个人的视频大全免费| 一a级毛片在线观看| 欧美一区二区亚洲| 亚洲国产色片| 大又大粗又爽又黄少妇毛片口| 女同久久另类99精品国产91| 久久精品夜夜夜夜夜久久蜜豆| 精品久久国产蜜桃| 国产一区二区激情短视频| 男女边吃奶边做爰视频| 亚洲av成人av| 久久精品综合一区二区三区| 国产av不卡久久| 亚洲熟妇中文字幕五十中出| 免费不卡的大黄色大毛片视频在线观看 | 国产精品永久免费网站| 亚洲一区高清亚洲精品| 亚洲精品456在线播放app| 桃色一区二区三区在线观看| 国产黄色视频一区二区在线观看 | 日韩高清综合在线| 深爱激情五月婷婷| 亚洲七黄色美女视频| 亚洲av第一区精品v没综合| 成人午夜高清在线视频| 夜夜爽天天搞| av卡一久久| 国产中年淑女户外野战色| 国产精品野战在线观看| 精品熟女少妇av免费看| 91麻豆精品激情在线观看国产| 久久韩国三级中文字幕| 麻豆av噜噜一区二区三区| 国产精品电影一区二区三区| 在线观看av片永久免费下载| 亚洲av免费在线观看| 两个人视频免费观看高清| 成人鲁丝片一二三区免费| 久久久久久久久中文| 久久精品人妻少妇| 久久人人爽人人爽人人片va| 日韩欧美精品免费久久| 色视频www国产| 少妇的逼水好多| 不卡一级毛片| 丰满的人妻完整版| 99久久无色码亚洲精品果冻| 美女高潮的动态| 成人精品一区二区免费| 又粗又爽又猛毛片免费看| 天堂动漫精品| 久久久久九九精品影院| 日韩大尺度精品在线看网址| 久久久久久九九精品二区国产| 国产精品三级大全| 免费观看的影片在线观看| 国内精品久久久久精免费| 国产亚洲精品久久久久久毛片| 色综合亚洲欧美另类图片| 一本久久中文字幕| 啦啦啦韩国在线观看视频| 三级国产精品欧美在线观看| 香蕉av资源在线| 一区二区三区免费毛片| 精品99又大又爽又粗少妇毛片| 看免费成人av毛片| 亚洲精品一区av在线观看| www.色视频.com| 美女内射精品一级片tv| 99热这里只有精品一区| 你懂的网址亚洲精品在线观看 | 亚洲成人久久爱视频| 成人综合一区亚洲| 1000部很黄的大片| 欧美一区二区国产精品久久精品| 桃色一区二区三区在线观看| 不卡视频在线观看欧美| 免费高清视频大片| 最新在线观看一区二区三区| 99久久九九国产精品国产免费| 日本 av在线| 午夜免费激情av| 黄色欧美视频在线观看| 色哟哟哟哟哟哟| 丰满乱子伦码专区| 国产视频内射| 伊人久久精品亚洲午夜| 国产不卡一卡二| 岛国在线免费视频观看| 国产v大片淫在线免费观看| 大香蕉久久网| 色综合亚洲欧美另类图片| 1024手机看黄色片| 少妇的逼水好多| 亚洲精品一区av在线观看| 黄色视频,在线免费观看| 日韩 亚洲 欧美在线| 久久久欧美国产精品| 久久人人精品亚洲av| 亚洲欧美日韩卡通动漫| 别揉我奶头 嗯啊视频| 亚洲欧美精品自产自拍| 免费高清视频大片| 好男人在线观看高清免费视频| 亚洲人成网站在线观看播放| 国产一级毛片七仙女欲春2| 十八禁网站免费在线| 亚洲人与动物交配视频| 99精品在免费线老司机午夜| 午夜福利在线在线| 精品免费久久久久久久清纯| 成年av动漫网址| 插逼视频在线观看| 国产av不卡久久| 午夜激情福利司机影院| 91狼人影院| 欧美日韩精品成人综合77777| 婷婷精品国产亚洲av在线| 亚洲av成人av| 亚洲一区高清亚洲精品| 三级男女做爰猛烈吃奶摸视频| 国产伦精品一区二区三区视频9| 午夜福利视频1000在线观看| 啦啦啦韩国在线观看视频| 2021天堂中文幕一二区在线观| 久久综合国产亚洲精品| 免费黄网站久久成人精品| 亚洲第一区二区三区不卡| 长腿黑丝高跟| 久久人妻av系列| 男女啪啪激烈高潮av片| 精品免费久久久久久久清纯| 亚洲欧美精品综合久久99| 国产精品久久久久久久电影| 成人特级av手机在线观看| 非洲黑人性xxxx精品又粗又长| 少妇熟女aⅴ在线视频| 久久婷婷人人爽人人干人人爱| 在线免费十八禁| 亚洲性久久影院| 精品一区二区免费观看| 成人性生交大片免费视频hd| 中文资源天堂在线| 99热网站在线观看| 好男人在线观看高清免费视频| 日日摸夜夜添夜夜添av毛片| 青春草视频在线免费观看| 人妻制服诱惑在线中文字幕| 久久久久久久久久黄片| 99热全是精品| 99视频精品全部免费 在线| 色综合站精品国产| 中文字幕免费在线视频6| 一区二区三区高清视频在线| 少妇裸体淫交视频免费看高清| 中文字幕免费在线视频6| 麻豆久久精品国产亚洲av| 观看美女的网站| 亚洲在线自拍视频| 国产一区二区三区av在线 | 深爱激情五月婷婷| 国产 一区 欧美 日韩| 亚洲真实伦在线观看| 成人无遮挡网站| 久久久欧美国产精品| 男女做爰动态图高潮gif福利片| 亚洲人与动物交配视频| 看十八女毛片水多多多| 国产国拍精品亚洲av在线观看| 特大巨黑吊av在线直播| 三级国产精品欧美在线观看| 日韩欧美免费精品| 99热这里只有是精品50| 色在线成人网| av女优亚洲男人天堂| 三级经典国产精品| 亚洲精品国产成人久久av| 国内久久婷婷六月综合欲色啪| 国产高潮美女av| av卡一久久| 在线播放无遮挡| 我的女老师完整版在线观看| 久久久久性生活片| 国产午夜福利久久久久久| 亚洲av一区综合| 久久久久久久亚洲中文字幕| 黄片wwwwww| 午夜激情欧美在线| 亚洲专区国产一区二区| 免费看美女性在线毛片视频| 久久久精品大字幕| 看非洲黑人一级黄片| 别揉我奶头~嗯~啊~动态视频| h日本视频在线播放| 两个人视频免费观看高清| 亚洲av中文字字幕乱码综合| 高清日韩中文字幕在线| 久久久国产成人精品二区| av在线蜜桃| 亚洲性夜色夜夜综合| 啦啦啦啦在线视频资源| 女同久久另类99精品国产91| 亚洲精品粉嫩美女一区| 午夜爱爱视频在线播放| 亚洲色图av天堂| 天美传媒精品一区二区| 一级毛片aaaaaa免费看小| 久久99热这里只有精品18| 免费看日本二区| 男人舔奶头视频| 午夜福利视频1000在线观看| 变态另类成人亚洲欧美熟女| 在线观看美女被高潮喷水网站| 国产激情偷乱视频一区二区| 久久久久久伊人网av| 欧美极品一区二区三区四区| 国产爱豆传媒在线观看| 免费观看的影片在线观看| av天堂中文字幕网| 亚洲av熟女| 日本熟妇午夜| 狂野欧美激情性xxxx在线观看| 精品一区二区三区av网在线观看| 精品久久久久久久久亚洲| 国产中年淑女户外野战色| 精品一区二区三区人妻视频| 一级毛片我不卡| 婷婷亚洲欧美| 欧美性感艳星| 99久久精品一区二区三区| 深夜精品福利| 亚洲国产精品国产精品| 精品久久久久久久久亚洲| 成年免费大片在线观看| 1000部很黄的大片| 久久九九热精品免费| 人人妻人人澡人人爽人人夜夜 | 欧美日本亚洲视频在线播放| 狂野欧美白嫩少妇大欣赏| 免费在线观看影片大全网站| 男人和女人高潮做爰伦理| 成人特级黄色片久久久久久久| 色播亚洲综合网| 国内精品一区二区在线观看| 黄色视频,在线免费观看| 免费一级毛片在线播放高清视频| 亚洲av成人av| 老师上课跳d突然被开到最大视频| 日韩中字成人| 欧美成人精品欧美一级黄| 亚洲成人精品中文字幕电影| 三级毛片av免费| 天天躁夜夜躁狠狠久久av| 亚洲欧美日韩高清专用| 精品不卡国产一区二区三区| 欧美在线一区亚洲| 欧美激情国产日韩精品一区| 亚洲人成网站在线播放欧美日韩| 在线观看午夜福利视频| 女同久久另类99精品国产91| 国产精品久久久久久av不卡| 大香蕉久久网| 成人午夜高清在线视频| 国产伦在线观看视频一区| 超碰av人人做人人爽久久| 久久久欧美国产精品| 久久久久久国产a免费观看| 国产精品电影一区二区三区| 色综合站精品国产| 麻豆一二三区av精品| 99在线人妻在线中文字幕| 可以在线观看毛片的网站| 日本精品一区二区三区蜜桃| 十八禁国产超污无遮挡网站| 丰满的人妻完整版| 精品99又大又爽又粗少妇毛片| 欧美日韩精品成人综合77777| 成年女人毛片免费观看观看9| 国产一区二区三区av在线 | 国产精品,欧美在线| 中文字幕免费在线视频6| 精品久久久久久久久亚洲| 欧美成人精品欧美一级黄| 色哟哟哟哟哟哟| 欧美在线一区亚洲| 大香蕉久久网| 免费大片18禁| 成人美女网站在线观看视频| 男人舔女人下体高潮全视频| 校园春色视频在线观看| www日本黄色视频网| 如何舔出高潮| 97人妻精品一区二区三区麻豆| 欧美中文日本在线观看视频| 久久国产乱子免费精品| 免费看日本二区| 成人三级黄色视频| 人妻久久中文字幕网| 黄片wwwwww| 亚洲国产精品久久男人天堂| 99热网站在线观看| 国产aⅴ精品一区二区三区波| 美女黄网站色视频| 久久久a久久爽久久v久久| 深夜精品福利| 免费观看精品视频网站| 国产一区二区激情短视频| 国产激情偷乱视频一区二区| 国产精品日韩av在线免费观看|