• <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
    久久久精品大字幕| 看免费成人av毛片| 日日摸夜夜添夜夜爱| 永久免费av网站大全| 精品少妇黑人巨大在线播放 | 国产单亲对白刺激| 亚洲精品,欧美精品| 黄色欧美视频在线观看| 久久精品久久精品一区二区三区| 国产精品福利在线免费观看| 国产精品.久久久| 91精品国产九色| 日日干狠狠操夜夜爽| 色吧在线观看| 亚洲精品aⅴ在线观看| 91精品一卡2卡3卡4卡| 国产片特级美女逼逼视频| 国产精品三级大全| 夜夜看夜夜爽夜夜摸| 国产91av在线免费观看| 少妇人妻一区二区三区视频| 老司机福利观看| 亚洲第一区二区三区不卡| 亚洲国产精品久久男人天堂| 国产在视频线精品| 午夜a级毛片| 国产精品电影一区二区三区| 亚洲人与动物交配视频| a级毛片免费高清观看在线播放| 日韩成人伦理影院| 国语自产精品视频在线第100页| 日韩 亚洲 欧美在线| 天美传媒精品一区二区| 国产伦精品一区二区三区视频9| 一级黄色大片毛片| 日本三级黄在线观看| 久99久视频精品免费| 亚洲精品国产成人久久av| 天天躁夜夜躁狠狠久久av| 精品国内亚洲2022精品成人| 日韩 亚洲 欧美在线| 女人十人毛片免费观看3o分钟| 神马国产精品三级电影在线观看| 天堂影院成人在线观看| 男女那种视频在线观看| 国产高清国产精品国产三级 | 亚洲精品影视一区二区三区av| 哪个播放器可以免费观看大片| 久久久久国产网址| 久久久久久久亚洲中文字幕| 秋霞在线观看毛片| av视频在线观看入口| 亚洲av男天堂| 午夜福利网站1000一区二区三区| 欧美人与善性xxx| 亚洲人成网站高清观看| 免费黄色在线免费观看| 91在线精品国自产拍蜜月| 亚洲真实伦在线观看| 日日干狠狠操夜夜爽| 91午夜精品亚洲一区二区三区| 免费av观看视频| 成人午夜高清在线视频| 网址你懂的国产日韩在线| 国产成人a区在线观看| 亚洲人成网站高清观看| 男人的好看免费观看在线视频| 蜜桃久久精品国产亚洲av| 免费观看人在逋| 国产高清视频在线观看网站| 久久99蜜桃精品久久| 亚洲国产精品合色在线| 国产精品嫩草影院av在线观看| 久久久精品欧美日韩精品| 波野结衣二区三区在线| 久久久久久久亚洲中文字幕| 桃色一区二区三区在线观看| 一级二级三级毛片免费看| h日本视频在线播放| 亚洲最大成人手机在线| 中文字幕av成人在线电影| 亚洲国产色片| 纵有疾风起免费观看全集完整版 | 欧美+日韩+精品| 国产精品美女特级片免费视频播放器| 日本-黄色视频高清免费观看| 日韩欧美国产在线观看| 能在线免费观看的黄片| 日韩大片免费观看网站 | 亚洲经典国产精华液单| 国产成人a区在线观看| 久久99热6这里只有精品| 国产高清国产精品国产三级 | 少妇裸体淫交视频免费看高清| 少妇人妻一区二区三区视频| 亚洲欧美日韩东京热| 日韩亚洲欧美综合| 永久网站在线| 能在线免费观看的黄片| 亚洲欧洲国产日韩| 乱码一卡2卡4卡精品| 亚洲国产欧洲综合997久久,| 伊人久久精品亚洲午夜| 一级黄片播放器| 亚洲欧美日韩高清专用| 日韩欧美精品免费久久| 丰满少妇做爰视频| 欧美性猛交╳xxx乱大交人| 高清日韩中文字幕在线| 中文天堂在线官网| 国产精品1区2区在线观看.| 欧美日韩国产亚洲二区| 啦啦啦观看免费观看视频高清| 欧美性猛交黑人性爽| 日韩av不卡免费在线播放| 亚洲av二区三区四区| 久久精品熟女亚洲av麻豆精品 | 身体一侧抽搐| 亚洲欧美精品综合久久99| av国产久精品久网站免费入址| 亚洲一级一片aⅴ在线观看| 美女大奶头视频| 村上凉子中文字幕在线| 天天躁夜夜躁狠狠久久av| 亚洲av日韩在线播放| 国产成人一区二区在线| 午夜福利成人在线免费观看| 成人av在线播放网站| 性色avwww在线观看| 日韩高清综合在线| 我要搜黄色片| 久久精品夜色国产| 欧美bdsm另类| 2021天堂中文幕一二区在线观| 亚洲国产高清在线一区二区三| 亚洲欧美成人精品一区二区| 国产精品麻豆人妻色哟哟久久 | 丰满少妇做爰视频| 亚洲国产欧美在线一区| 久久婷婷人人爽人人干人人爱| 神马国产精品三级电影在线观看| 大话2 男鬼变身卡| www.av在线官网国产| 日韩一区二区视频免费看| 观看免费一级毛片| 又爽又黄无遮挡网站| 简卡轻食公司| 国产片特级美女逼逼视频| 亚洲高清免费不卡视频| 日韩强制内射视频| 亚洲国产精品合色在线| 大又大粗又爽又黄少妇毛片口| 久久久成人免费电影| 一个人观看的视频www高清免费观看| 纵有疾风起免费观看全集完整版 | 国产精品久久久久久精品电影小说 | 熟妇人妻久久中文字幕3abv| 欧美+日韩+精品| 91狼人影院| 国产成人午夜福利电影在线观看| 欧美日本视频| 三级国产精品片| 少妇人妻精品综合一区二区| 国产中年淑女户外野战色| 观看免费一级毛片| av专区在线播放| 久久久久久国产a免费观看| 爱豆传媒免费全集在线观看| 少妇丰满av| 男人舔女人下体高潮全视频| 久久久久久久久大av| 亚洲欧美日韩卡通动漫| 成人二区视频| 亚洲久久久久久中文字幕| 精品久久久久久成人av| 久久久久久久国产电影| 欧美成人午夜免费资源| 国产精品久久久久久久久免| 成年av动漫网址| 免费黄网站久久成人精品| 伦理电影大哥的女人| 精品国产一区二区三区久久久樱花 | 在线观看av片永久免费下载| 国产一级毛片七仙女欲春2| 不卡视频在线观看欧美| 波野结衣二区三区在线| 亚洲欧美精品自产自拍| 在线观看av片永久免费下载| www日本黄色视频网| 亚洲综合色惰| 国产亚洲精品久久久com| 女人十人毛片免费观看3o分钟| 最近视频中文字幕2019在线8| 国产视频首页在线观看| 天美传媒精品一区二区| 最后的刺客免费高清国语| 不卡视频在线观看欧美| 免费看日本二区| 别揉我奶头 嗯啊视频| 中文在线观看免费www的网站| 国产色爽女视频免费观看| 十八禁国产超污无遮挡网站| 精品久久久久久久久久久久久| 色尼玛亚洲综合影院| 亚洲av成人精品一区久久| 亚洲人成网站高清观看| 99久国产av精品| 国产精品久久电影中文字幕| 尤物成人国产欧美一区二区三区| 国产精品麻豆人妻色哟哟久久 | 日韩 亚洲 欧美在线| 啦啦啦啦在线视频资源| 亚洲图色成人| 久久久精品大字幕| 99久久精品一区二区三区| 日日摸夜夜添夜夜添av毛片| 国产亚洲精品av在线| 伦精品一区二区三区| 欧美激情久久久久久爽电影| 久久久久久久久久久免费av| 可以在线观看毛片的网站| 久久人妻av系列| 成人一区二区视频在线观看| 国产亚洲精品av在线| 伊人久久精品亚洲午夜| 精品熟女少妇av免费看| 日韩三级伦理在线观看| 日本av手机在线免费观看| 六月丁香七月| 精品人妻熟女av久视频| 成年免费大片在线观看| 亚洲久久久久久中文字幕| 成人综合一区亚洲| 亚洲性久久影院| 亚洲成人av在线免费| 身体一侧抽搐| 日本午夜av视频| 亚洲国产欧美在线一区| 视频中文字幕在线观看| 国产亚洲av嫩草精品影院| 久久99热6这里只有精品| 久久精品国产亚洲网站| 最近中文字幕高清免费大全6| 欧美xxxx黑人xx丫x性爽| 亚洲国产色片| 99久国产av精品| 赤兔流量卡办理| 九草在线视频观看| 中文字幕av成人在线电影| 老司机影院毛片| 最近中文字幕2019免费版| 久久久亚洲精品成人影院| 神马国产精品三级电影在线观看| 26uuu在线亚洲综合色| 精品久久国产蜜桃| 少妇猛男粗大的猛烈进出视频 | 亚洲欧美日韩东京热| 三级国产精品欧美在线观看| 五月伊人婷婷丁香| 亚洲欧洲日产国产| 亚洲精品影视一区二区三区av| 简卡轻食公司| 国产乱人偷精品视频| 国产伦在线观看视频一区| 成人毛片60女人毛片免费| 少妇丰满av| 淫秽高清视频在线观看| 亚洲最大成人手机在线| 国产v大片淫在线免费观看| 亚洲在线观看片| 国内精品一区二区在线观看| 国产高潮美女av| 最新中文字幕久久久久| 国产亚洲午夜精品一区二区久久 | 亚洲人成网站在线观看播放| 波野结衣二区三区在线| 只有这里有精品99| 长腿黑丝高跟| 好男人视频免费观看在线| 桃色一区二区三区在线观看| 别揉我奶头 嗯啊视频| 亚洲国产日韩欧美精品在线观看| 国产不卡一卡二| 国产精品熟女久久久久浪| 日韩视频在线欧美| 日本与韩国留学比较| 嫩草影院入口| 秋霞伦理黄片| 26uuu在线亚洲综合色| 久久久久网色| 波多野结衣高清无吗| 99久久精品热视频| 深爱激情五月婷婷| 能在线免费看毛片的网站| 桃色一区二区三区在线观看| 中文字幕久久专区| or卡值多少钱| 插阴视频在线观看视频| 亚洲av电影不卡..在线观看| 中文字幕av在线有码专区| 国产午夜福利久久久久久| 九九久久精品国产亚洲av麻豆| 校园人妻丝袜中文字幕| 国产午夜精品久久久久久一区二区三区| 久久99精品国语久久久| 国产成年人精品一区二区| 卡戴珊不雅视频在线播放| 欧美日韩精品成人综合77777| 国产探花极品一区二区| 国产精品国产高清国产av| 99久久九九国产精品国产免费| 国产精华一区二区三区| 99久国产av精品| 国产又色又爽无遮挡免| 国产午夜福利久久久久久| 国产精品久久久久久精品电影| 美女xxoo啪啪120秒动态图| 亚洲成人久久爱视频| 免费观看精品视频网站| 亚洲中文字幕日韩| 亚洲美女视频黄频| 久久精品久久精品一区二区三区| 韩国高清视频一区二区三区| 秋霞伦理黄片| 国内少妇人妻偷人精品xxx网站| 两个人视频免费观看高清| 秋霞在线观看毛片| 久久久久久久久中文| 亚洲av日韩在线播放| 色吧在线观看| 久99久视频精品免费| 国产亚洲av片在线观看秒播厂 | 中文天堂在线官网| 97人妻精品一区二区三区麻豆| 噜噜噜噜噜久久久久久91| 久久久久九九精品影院| 青春草视频在线免费观看| 国产av在哪里看| 国产精品国产三级国产av玫瑰| av在线天堂中文字幕| 国产精品熟女久久久久浪| 久久人妻av系列| 人体艺术视频欧美日本| 免费播放大片免费观看视频在线观看 | 国内少妇人妻偷人精品xxx网站| 久久亚洲精品不卡| 国产高清有码在线观看视频| 久久精品国产99精品国产亚洲性色| 久久亚洲国产成人精品v| 欧美又色又爽又黄视频| av免费观看日本| 亚洲国产欧美在线一区| 最后的刺客免费高清国语| 欧美丝袜亚洲另类| 麻豆一二三区av精品| 我的老师免费观看完整版| 26uuu在线亚洲综合色| 99久久无色码亚洲精品果冻| 嘟嘟电影网在线观看| 日韩国内少妇激情av| 久久99蜜桃精品久久| 国产精品精品国产色婷婷| 免费黄色在线免费观看| 秋霞伦理黄片| 高清在线视频一区二区三区 | 欧美一区二区精品小视频在线| 纵有疾风起免费观看全集完整版 | 男女啪啪激烈高潮av片| 成人鲁丝片一二三区免费| a级毛片免费高清观看在线播放| 在线a可以看的网站| 又爽又黄a免费视频| 久久精品国产99精品国产亚洲性色| 国产 一区精品| 一级毛片电影观看 | 国产精品久久电影中文字幕| 国产一级毛片在线| 国产成人精品久久久久久| 青青草视频在线视频观看| 久久99热6这里只有精品| 亚洲成av人片在线播放无| 亚洲精品久久久久久婷婷小说 | 久久久精品大字幕| 天堂中文最新版在线下载 | av视频在线观看入口| 超碰97精品在线观看| 亚洲在久久综合| 婷婷色综合大香蕉| 国产高清有码在线观看视频| 老师上课跳d突然被开到最大视频| 精品国产一区二区三区久久久樱花 | 禁无遮挡网站| 亚洲人成网站在线播| 亚洲经典国产精华液单| 免费av毛片视频| 午夜福利视频1000在线观看| 国产又黄又爽又无遮挡在线| 午夜福利在线观看免费完整高清在| 国产男人的电影天堂91| 51国产日韩欧美| 亚洲无线观看免费| 日本色播在线视频| 少妇猛男粗大的猛烈进出视频 | av又黄又爽大尺度在线免费看 | 国产v大片淫在线免费观看| 两性午夜刺激爽爽歪歪视频在线观看| 26uuu在线亚洲综合色| 卡戴珊不雅视频在线播放| 十八禁国产超污无遮挡网站| av在线蜜桃| 九九热线精品视视频播放| 成人鲁丝片一二三区免费| 亚洲最大成人中文| 搡女人真爽免费视频火全软件| 一本久久精品| 午夜精品一区二区三区免费看| 永久免费av网站大全| 色哟哟·www| 日韩一区二区三区影片| 熟妇人妻久久中文字幕3abv| 在线观看一区二区三区| 高清av免费在线| 超碰av人人做人人爽久久| 水蜜桃什么品种好| 天堂网av新在线| 国产精品一区二区三区四区免费观看| 国产老妇女一区| 亚洲精品日韩av片在线观看| 亚洲五月天丁香| 干丝袜人妻中文字幕| 精品久久久噜噜| 波野结衣二区三区在线| 亚洲成人精品中文字幕电影| 国产av不卡久久| 久久午夜福利片| 蜜桃亚洲精品一区二区三区| 国产成年人精品一区二区| 天天躁日日操中文字幕| 久久精品熟女亚洲av麻豆精品 | 少妇的逼水好多| 国产精品久久久久久精品电影| 国产精品美女特级片免费视频播放器| 91在线精品国自产拍蜜月| 亚洲国产精品sss在线观看| 久久精品久久精品一区二区三区| 国产在视频线在精品| 国内少妇人妻偷人精品xxx网站| 激情 狠狠 欧美| 国产久久久一区二区三区| 国产av码专区亚洲av| 亚洲av成人精品一二三区| 午夜a级毛片| 国产成人a∨麻豆精品| 啦啦啦观看免费观看视频高清| 22中文网久久字幕| 偷拍熟女少妇极品色| 国产中年淑女户外野战色| 国产午夜精品论理片| 亚洲av男天堂| 国产美女午夜福利| 国产免费男女视频| 国产伦精品一区二区三区视频9| 亚洲四区av| 精品国内亚洲2022精品成人| 永久免费av网站大全| 免费看美女性在线毛片视频| 国产又黄又爽又无遮挡在线| av在线蜜桃| 搡老妇女老女人老熟妇| 亚洲人成网站高清观看| 午夜福利在线在线| 好男人视频免费观看在线| 国产激情偷乱视频一区二区| 欧美bdsm另类| 欧美色视频一区免费| 欧美日韩国产亚洲二区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产亚洲一区二区精品| 美女脱内裤让男人舔精品视频| 欧美性感艳星| 国产 一区 欧美 日韩| 免费观看性生交大片5| 18禁动态无遮挡网站| 亚洲av日韩在线播放| 久久精品久久精品一区二区三区| 国产欧美日韩精品一区二区| 国产不卡一卡二| 少妇人妻精品综合一区二区| 免费av不卡在线播放| 中文字幕av在线有码专区| 搡女人真爽免费视频火全软件| 欧美xxxx性猛交bbbb| 免费av不卡在线播放| 18禁在线播放成人免费| 国产又黄又爽又无遮挡在线| 久久国产乱子免费精品| 简卡轻食公司| 国产高清不卡午夜福利| 国产精品女同一区二区软件| 夜夜爽夜夜爽视频| 国内揄拍国产精品人妻在线| 三级毛片av免费| 夜夜看夜夜爽夜夜摸| 欧美一区二区亚洲| 青春草视频在线免费观看| 国产免费一级a男人的天堂| 亚洲欧美成人综合另类久久久 | 岛国毛片在线播放| 久久欧美精品欧美久久欧美| 亚洲av熟女| 亚洲av中文字字幕乱码综合| 日韩大片免费观看网站 | 永久免费av网站大全| 97在线视频观看| 最近最新中文字幕大全电影3| 国产午夜精品久久久久久一区二区三区| 久久国内精品自在自线图片| 高清视频免费观看一区二区 | 久久精品夜夜夜夜夜久久蜜豆| 国产黄色视频一区二区在线观看 | 97在线视频观看| 日本三级黄在线观看| 国产午夜精品久久久久久一区二区三区| videos熟女内射| 成人三级黄色视频| 非洲黑人性xxxx精品又粗又长| 国产老妇女一区| 午夜亚洲福利在线播放| 国产精品野战在线观看| 国产一区二区在线av高清观看| 一卡2卡三卡四卡精品乱码亚洲| 亚洲国产精品sss在线观看| 一本一本综合久久| 久久久久久久久久久免费av| 国产精品女同一区二区软件| 99久久精品一区二区三区| 国产国拍精品亚洲av在线观看| 国产成人freesex在线| 国产亚洲最大av| 乱码一卡2卡4卡精品| 午夜a级毛片| 色视频www国产| 中文天堂在线官网| 久久久久九九精品影院| 九草在线视频观看| 我要看日韩黄色一级片| 国产在线男女| 国产成人精品婷婷| 久久久久久久久久久丰满| 亚洲精品日韩av片在线观看| 国产精品国产三级国产专区5o | 午夜福利在线观看吧| 国产免费男女视频| 我要看日韩黄色一级片| 国产淫语在线视频| av又黄又爽大尺度在线免费看 | 免费看日本二区| 男女啪啪激烈高潮av片| 老司机影院成人| av.在线天堂| 欧美成人一区二区免费高清观看| 黄色一级大片看看| 天美传媒精品一区二区| 在线播放无遮挡| 国产精品爽爽va在线观看网站| av专区在线播放| 国产亚洲91精品色在线| 国产成人精品一,二区| 久久久久久久久大av| 日本三级黄在线观看| 久久久久精品久久久久真实原创| 婷婷色综合大香蕉| 国产美女午夜福利| 成人毛片a级毛片在线播放| 国产淫片久久久久久久久| 两个人的视频大全免费| 精品欧美国产一区二区三| 国产免费福利视频在线观看| 国产亚洲精品久久久com| 夜夜看夜夜爽夜夜摸| 国产麻豆成人av免费视频| 久久国产乱子免费精品| 中文字幕亚洲精品专区| 久久精品人妻少妇| 好男人在线观看高清免费视频| 国产成人精品久久久久久| 国产69精品久久久久777片| 婷婷色麻豆天堂久久 | 在线a可以看的网站| 国产在视频线精品| 亚洲第一区二区三区不卡| 男人和女人高潮做爰伦理| 国产免费一级a男人的天堂| 欧美成人一区二区免费高清观看| 国产精品.久久久| 亚洲五月天丁香| 亚洲精品成人久久久久久| 晚上一个人看的免费电影| 国产 一区精品| 午夜精品国产一区二区电影 | 乱系列少妇在线播放| 日韩欧美 国产精品| 1024手机看黄色片| 国内少妇人妻偷人精品xxx网站| 亚洲欧美中文字幕日韩二区| 国产乱人视频| 青春草国产在线视频| 欧美日本亚洲视频在线播放| 免费看av在线观看网站| 精华霜和精华液先用哪个| 啦啦啦啦在线视频资源| 日本午夜av视频| 免费黄网站久久成人精品| 免费av毛片视频| 亚洲人成网站高清观看|