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

    網(wǎng)架設(shè)計(jì)軟件STADS的計(jì)算風(fēng)工程開發(fā)
    ——風(fēng)振分析

    2022-03-29 03:06:46白澤升王孟鴻
    關(guān)鍵詞:大渦風(fēng)振網(wǎng)架

    白澤升, 王孟鴻

    (北京建筑大學(xué) 土木與交通工程學(xué)院, 北京 100044)

    實(shí)際風(fēng)荷載通常包含長周期部分和短周期部分,其中長周期部分通常以平均風(fēng)速描述,短周期部分以脈動(dòng)風(fēng)來描述。對(duì)結(jié)構(gòu)物來說,由平均風(fēng)速引起的作用力相當(dāng)于靜力作用,不同部位通過體形系數(shù)來計(jì)算;脈動(dòng)風(fēng)速引起的作用力是動(dòng)力性質(zhì)的,將引起結(jié)構(gòu)振動(dòng),通過風(fēng)振系數(shù)來計(jì)算。荷載規(guī)范[1]規(guī)定:對(duì)于基本自振周期T1大于0.25 s的工程結(jié)構(gòu),如房屋、屋蓋及各種高聳結(jié)構(gòu),以及對(duì)于高度大于30 m且高寬比大于1.5的高柔房屋,均應(yīng)考慮風(fēng)壓脈動(dòng)對(duì)結(jié)構(gòu)發(fā)生順風(fēng)向風(fēng)振的影響。規(guī)范對(duì)于復(fù)雜的空間網(wǎng)架結(jié)構(gòu)的風(fēng)振系數(shù)沒有給出相應(yīng)的數(shù)據(jù)和計(jì)算方法。因此,復(fù)雜網(wǎng)架結(jié)構(gòu)的風(fēng)振計(jì)算通常需要借助相關(guān)理論,通過計(jì)算分析確定[2]。

    目前對(duì)復(fù)雜網(wǎng)架結(jié)構(gòu)進(jìn)行風(fēng)振響應(yīng)分析,主要采取對(duì)結(jié)構(gòu)的周圍風(fēng)場進(jìn)行模擬,以獲取結(jié)構(gòu)表面風(fēng)壓數(shù)據(jù),以進(jìn)行后續(xù)的風(fēng)振效應(yīng)分析;而風(fēng)洞試驗(yàn)一般作為剛性模型,可對(duì)結(jié)構(gòu)表面風(fēng)壓分布進(jìn)行測試,直接測試結(jié)構(gòu)的風(fēng)振效應(yīng)較為困難。對(duì)結(jié)構(gòu)周圍風(fēng)場的模擬方法主要有風(fēng)洞試驗(yàn)法和數(shù)值風(fēng)洞法。風(fēng)洞試驗(yàn)法周期長、成本高,試驗(yàn)數(shù)據(jù)需要經(jīng)過復(fù)雜的數(shù)據(jù)處理后導(dǎo)入分析軟件進(jìn)行分析計(jì)算。數(shù)值風(fēng)洞法成本低、周期短、效率高,便于參數(shù)化控制,但需要復(fù)雜的前處理建模工作,而網(wǎng)架結(jié)構(gòu)形式日益多樣化、異型化,無疑增加了建模工作量,同時(shí)計(jì)算流體動(dòng)力學(xué)(CFD)軟件與結(jié)構(gòu)分析軟件脫節(jié),網(wǎng)架結(jié)構(gòu)桿件與節(jié)點(diǎn)眾多,需要大量的荷載施加工作。

    本接口軟件可將復(fù)雜網(wǎng)架模型自動(dòng)導(dǎo)入前處理軟件ICEM,節(jié)省了大量的建模工作,而后處理功能基于FLUENT計(jì)算風(fēng)工程功能,自動(dòng)完成數(shù)據(jù)的讀取以及荷載施加與計(jì)算過程,減少了大量手動(dòng)施加工作。

    1 風(fēng)振分析方法

    由于風(fēng)荷載的隨機(jī)特性,目前大跨空間網(wǎng)架結(jié)構(gòu)的風(fēng)振分析方法主要有頻域法和時(shí)域法。頻域法是通過Fourier變換將風(fēng)壓時(shí)程轉(zhuǎn)化為風(fēng)壓功率譜,再將風(fēng)壓功率譜通過動(dòng)力系數(shù)傳遞轉(zhuǎn)化為動(dòng)力反應(yīng)譜,并根據(jù)隨機(jī)振動(dòng)理論通過對(duì)動(dòng)力反應(yīng)譜進(jìn)行積分求解得到結(jié)構(gòu)的動(dòng)力響應(yīng)。因頻域法計(jì)算效率較高,是目前大多數(shù)學(xué)者常采用的方法[3-6],其進(jìn)行的是線性范圍內(nèi)的分析。而大跨網(wǎng)架結(jié)構(gòu)具有頻率密集性,需通過振型分解解耦處理,考慮振型模態(tài)對(duì)系統(tǒng)應(yīng)變能的貢獻(xiàn)[8],選取主要貢獻(xiàn)模態(tài),可能會(huì)損失一部分精度。時(shí)域法則是將結(jié)構(gòu)風(fēng)壓荷載時(shí)程作用于大跨網(wǎng)架結(jié)構(gòu)進(jìn)行風(fēng)振響應(yīng)動(dòng)力時(shí)程分析,獲得動(dòng)力時(shí)程響應(yīng)的全過程。由于時(shí)域法能夠反映風(fēng)荷載的時(shí)空特性,并能夠用于結(jié)構(gòu)的非線性分析。因此,本文采用了時(shí)域法。

    進(jìn)行時(shí)域內(nèi)的風(fēng)振分析,需要確定作用在單元節(jié)點(diǎn)上的風(fēng)荷載時(shí)程,通常有2種方法。其一采用風(fēng)洞試驗(yàn)的風(fēng)壓時(shí)程,也是本文所采用的方法:從FLUENT模擬中提取風(fēng)荷載時(shí)程,并將其等效為節(jié)點(diǎn)動(dòng)力時(shí)程作用在相應(yīng)的網(wǎng)架節(jié)點(diǎn)上。其二是對(duì)脈動(dòng)風(fēng)進(jìn)行模擬,使模擬出的風(fēng)速盡可能接近和滿足自然風(fēng)特性如平均值、與高度有關(guān)的自功率譜和互功率譜以及相位角關(guān)系等,常用的方法有諧波疊加法和線性濾波器法。通過模擬得到的風(fēng)速時(shí)程在準(zhǔn)定常假定的基礎(chǔ)上轉(zhuǎn)換成風(fēng)荷載的時(shí)程。如何模擬合成輸入風(fēng)速以獲得合理的風(fēng)壓時(shí)程文件需要做進(jìn)一步的研究工作。

    2 STADS風(fēng)振分析流程

    2.1 湍流數(shù)值模擬方法

    CFD數(shù)值模擬方法可分為直接數(shù)值模擬法(DNS)、雷諾平均法(RANS)和大渦模擬法(LES)。直接數(shù)值模擬法(DNS)是直接采用瞬時(shí)的N-S方程對(duì)各尺度的渦流進(jìn)行求解,需要耗費(fèi)較高的計(jì)算資源,目前無法應(yīng)用于真正意義上的工程計(jì)算。雷諾平均法(RANS)將湍流運(yùn)動(dòng)時(shí)均化處理,引入各湍流模型對(duì)流體控制方程進(jìn)行封閉,廣泛應(yīng)用于穩(wěn)態(tài)風(fēng)場求解中。大渦模擬法(LES)的原理是將湍流中的渦旋分成大渦和小渦,用瞬時(shí)N-S方程對(duì)大渦直接求解,采用亞格子模型考慮小渦的能量輸運(yùn)與耗散,實(shí)質(zhì)是對(duì)湍流運(yùn)動(dòng)進(jìn)行空間平均化。由于大渦模擬法(LES)對(duì)高雷諾數(shù)的復(fù)雜湍流運(yùn)動(dòng)具有較高的分辨率,可以捕捉到復(fù)雜的非定常流動(dòng)現(xiàn)象及其各物理量,準(zhǔn)確度較高,是目前求解非定常流較為理想的湍流數(shù)值模擬方法。故本文采用了大渦模擬法(LES)。

    2.2 瞬態(tài)場風(fēng)壓計(jì)算

    基于大渦模擬法(LES)的瞬態(tài)場風(fēng)壓計(jì)算程序:

    1) 模型轉(zhuǎn)換,將STADS的模型轉(zhuǎn)換到ANSYS ICEM CFD[8-9]。

    2) 瞬態(tài)風(fēng)壓計(jì)算,通過FLUENT計(jì)算并輸出控制點(diǎn)的風(fēng)壓時(shí)程[8,10]。

    3) 計(jì)算節(jié)點(diǎn)的風(fēng)壓力時(shí)程,通過節(jié)點(diǎn)的控制面積和風(fēng)壓計(jì)算。

    4) 求解振型模態(tài)[11]計(jì)算結(jié)構(gòu)自振周期,利用前兩階振型計(jì)算阻尼系數(shù)。

    5) 利用Newmark法進(jìn)行時(shí)程計(jì)算,獲得各節(jié)點(diǎn)位移時(shí)程。

    6) 根據(jù)上一步的位移時(shí)程,計(jì)算各個(gè)節(jié)點(diǎn)的風(fēng)振系數(shù)。

    3 算例分析

    本文以跨度93 m、長度94 m的大型扁筒殼為例展示STADS風(fēng)振計(jì)算的全部過程。在此需要說明:由于監(jiān)控點(diǎn)太多,造成計(jì)算時(shí)間過長。因此,除進(jìn)行模型分區(qū)處理以外,程序還增加了計(jì)算輸出監(jiān)控點(diǎn)時(shí)程的歸并處理。通過定義數(shù)個(gè)相鄰分區(qū)共用一個(gè)區(qū)域的風(fēng)壓時(shí)程[12],減少數(shù)據(jù)輸出,從而加快程序運(yùn)行。

    計(jì)算模型如圖1所示:跨度93 m,長度94 m,節(jié)點(diǎn)數(shù)1 472,桿件5 724,上弦塊數(shù)為712;經(jīng)分區(qū),上弦劃分為142塊區(qū)域如圖1所示。

    圖1 計(jì)算模型Fig.1 Computational model

    3.1 模型轉(zhuǎn)換與瞬態(tài)風(fēng)壓計(jì)算

    通過選擇風(fēng)振控制項(xiàng)自動(dòng)選擇大渦模擬法(LES)進(jìn)行計(jì)算,根據(jù)ICEM的腳本語言和函數(shù)控制格式生成模型文件,根據(jù)FLUENT二次開發(fā)方式,生成UDF和腳本控制文件。相應(yīng)參數(shù)控制,如圖2所示。

    圖2 參數(shù)控制Fig.2 Parameter control

    其中,大渦模擬亞格子模型采用壁面適應(yīng)局部渦黏模型(WALE),對(duì)流場的整體流動(dòng)以及分離、附著、轉(zhuǎn)捩等局部流動(dòng)模擬較為準(zhǔn)確。壓力與速度耦合采用SIMPLEC算法,對(duì)流項(xiàng)采用有限中心差分格式離散,時(shí)間項(xiàng)離散采用有限二階精度格式,壓力方程采用二階精度離散格式。

    瞬態(tài)分析中,時(shí)間步長的設(shè)置對(duì)大渦模擬的結(jié)果影響很大,決定了模擬結(jié)果的時(shí)間分辨度,是大渦模擬收斂情況的一個(gè)重要參數(shù),同時(shí)也要保證數(shù)值模擬能在一個(gè)合理時(shí)間內(nèi)完成。針對(duì)以上標(biāo)準(zhǔn)進(jìn)行試算,計(jì)算步長控制為0.06 s,滿足模擬需要。

    在ICEM CFD中讀取模型文件自動(dòng)完成模型的轉(zhuǎn)換和網(wǎng)格劃分等前處理功能,程序通過UDF定義邊界條件,運(yùn)行FLUENT腳本控制文件后,自動(dòng)完成模型的選擇、監(jiān)控點(diǎn)的設(shè)置、求解參數(shù)設(shè)置等操作,并輸出控制點(diǎn)的風(fēng)壓時(shí)程數(shù)據(jù)。控制點(diǎn)風(fēng)壓時(shí)程如圖3所示。100、137號(hào)為山墻端部節(jié)點(diǎn),30、34、70號(hào)分別是筒殼迎風(fēng)面、頂部和背風(fēng)面節(jié)點(diǎn)。

    圖3 風(fēng)壓- 時(shí)間曲線Fig.3 Wind pressure-time curve

    3.2 節(jié)點(diǎn)風(fēng)壓力時(shí)程的生成

    時(shí)程分析前需要將FLUENT計(jì)算得到的控制點(diǎn)的風(fēng)壓時(shí)程轉(zhuǎn)換成網(wǎng)架節(jié)點(diǎn)的風(fēng)壓力時(shí)程,該過程涉及復(fù)雜每個(gè)網(wǎng)格承載面積的計(jì)算,通過STADS-2加載程序自動(dòng)完成該功能,從而形成每一個(gè)節(jié)點(diǎn)的風(fēng)壓時(shí)程??缰?52、769、770號(hào)節(jié)點(diǎn)風(fēng)壓時(shí)程如圖4所示。

    圖4 節(jié)點(diǎn)風(fēng)壓- 時(shí)間曲線Fig.4 Node wind pressure-time curve

    3.3 自振周期和阻尼系數(shù)計(jì)算

    利用子空間迭代求得前兩階振型周期:T1=0.641 s,T2=0.529 s;何艷麗等基于各基本參數(shù)對(duì)某筒殼風(fēng)振響應(yīng)進(jìn)行了研究,隨著阻尼比增大,各節(jié)點(diǎn)風(fēng)振系數(shù)逐步遞減[13],參考其阻尼比對(duì)風(fēng)振系數(shù)的影響和鋼結(jié)構(gòu)阻尼比取值經(jīng)驗(yàn),阻尼比取值ξ=0.025;阻尼矩陣采用Rayleigh阻尼矩陣:[C]=α*[M]+β*[K],其中,[C]、[M]、[K]分別為阻尼、質(zhì)量和剛度矩陣,α和β為比例系數(shù)。程序據(jù)依前兩階角頻率算得:α=0.268 2,β=0.002 3,如圖5所示。

    圖5 計(jì)算系數(shù)求解設(shè)置Fig.5 Calculation coefficient solution setting

    3.4 時(shí)程分析計(jì)算

    時(shí)程分析計(jì)算采用高級(jí)分析軟件STADS-3進(jìn)行。程序讀取節(jié)點(diǎn)風(fēng)荷載時(shí)程,采用Newmark法進(jìn)行時(shí)程分析,計(jì)算得到各節(jié)點(diǎn)位移時(shí)程,以用于下一步計(jì)算各個(gè)節(jié)點(diǎn)的風(fēng)振系數(shù)。計(jì)算所得的跨中727和778號(hào)節(jié)點(diǎn)z向位移時(shí)程如圖6所示。

    圖6 節(jié)點(diǎn)位移- 時(shí)間曲線Fig.6 Node displacement-time curve

    3.5 風(fēng)振系數(shù)的計(jì)算

    根據(jù)荷載規(guī)范[1]定義,在工程應(yīng)用中將風(fēng)荷載的動(dòng)力效應(yīng)以風(fēng)振系數(shù)β的形式等效為靜力荷載,即風(fēng)振系數(shù)為:

    (1)

    采用峰值因子法,由平衡方程,荷載風(fēng)振系數(shù)通過剛度矩陣和位移表示為:

    (2)

    對(duì)于節(jié)點(diǎn)i,荷載風(fēng)振系數(shù)βi為:

    (3)

    式中:j為i節(jié)點(diǎn)的相鄰節(jié)點(diǎn),kij為總剛矩陣第i行的非對(duì)角元。

    采用峰值因子法,可將節(jié)點(diǎn)i處的位移風(fēng)振系數(shù)表示為:

    (4)

    上述計(jì)算過程通過菜單引導(dǎo)執(zhí)行,將各節(jié)點(diǎn)風(fēng)振系數(shù)記錄在相應(yīng)文件中。用戶可以選擇執(zhí)行荷載風(fēng)振系數(shù)式(3)或者位移風(fēng)振系數(shù)式(4),其中在計(jì)算荷載風(fēng)振系數(shù)時(shí)需要調(diào)用整體剛度矩陣。筒殼表面和山墻的位移風(fēng)振系數(shù)如圖7所示。

    3.6 風(fēng)振系數(shù)的利用

    上述無論荷載風(fēng)振系數(shù)還是位移風(fēng)振系數(shù),都僅作為設(shè)計(jì)參考,需要設(shè)計(jì)人員做出合理判斷使用??紤]到其加載過程手工施加的工作量較大,在STADS-2的加載過程中可以自動(dòng)讀入各節(jié)點(diǎn)風(fēng)振系數(shù)數(shù)據(jù),可避免手工輸入的煩瑣與可能產(chǎn)生的數(shù)據(jù)錯(cuò)誤,減小了相應(yīng)工作量。

    4 結(jié) 論

    1) 在STADS網(wǎng)架設(shè)計(jì)軟件基礎(chǔ)上,借助FLUENT軟件的風(fēng)工程計(jì)算功能,編制了軟件模型轉(zhuǎn)換和風(fēng)壓讀取接口,并在高級(jí)分析軟件STADS-3的動(dòng)力時(shí)程分析基礎(chǔ)上,完成了復(fù)雜體型網(wǎng)架結(jié)構(gòu)的風(fēng)振分析。本接口軟件的開發(fā)可為此類鋼結(jié)構(gòu)軟件的風(fēng)振分析開發(fā)提供理論依據(jù)與思路參考。

    2) 接口程序通過菜單引導(dǎo)控制相應(yīng)參數(shù)自動(dòng)完成模型轉(zhuǎn)換、FLUENT風(fēng)場數(shù)值模擬參數(shù)設(shè)置和大量控制點(diǎn)數(shù)據(jù)的輸入與計(jì)算結(jié)果的提取,以及荷載的自動(dòng)讀取與施加,操作簡便,減少了大量手動(dòng)輸入的煩瑣工作,節(jié)約了數(shù)值風(fēng)洞環(huán)節(jié)的時(shí)間和經(jīng)濟(jì)成本。將FLUENT與STADS建立接口,有效彌補(bǔ)了計(jì)算流體動(dòng)力學(xué)軟件(CFD)與結(jié)構(gòu)分析軟件的脫節(jié),加快了程序的運(yùn)行,提高了工作效率。

    3) 接口程序以某實(shí)際工程的筒殼模型為算例,計(jì)算得到網(wǎng)架結(jié)構(gòu)各節(jié)點(diǎn)位移風(fēng)振系數(shù)或荷載風(fēng)振系數(shù),在實(shí)際工程設(shè)計(jì)中可通過靜力加載方式計(jì)算得到結(jié)構(gòu)響應(yīng),簡化了工程設(shè)計(jì)工作,為復(fù)雜體型網(wǎng)架結(jié)構(gòu)抗風(fēng)設(shè)計(jì)提供參考。同時(shí)可適應(yīng)各類異型空間網(wǎng)架結(jié)構(gòu),具有良好的實(shí)用前景。

    猜你喜歡
    大渦風(fēng)振網(wǎng)架
    屋頂大型鋼網(wǎng)架球節(jié)點(diǎn)測量應(yīng)用研究
    基于壁面射流的下?lián)舯┝鞣欠€(wěn)態(tài)風(fēng)場大渦模擬
    張力對(duì)接觸網(wǎng)系統(tǒng)風(fēng)振響應(yīng)影響的風(fēng)洞試驗(yàn)研究
    軸流風(fēng)機(jī)葉尖泄漏流動(dòng)的大渦模擬
    高塔設(shè)備風(fēng)振失效原因分析及改善措施
    網(wǎng)架結(jié)構(gòu)設(shè)計(jì)軟件對(duì)比分析
    某網(wǎng)架桿件彎曲的原因分析及處理
    基于大渦模擬的旋風(fēng)分離器錐體結(jié)構(gòu)影響研究
    大跨屋蓋結(jié)構(gòu)MTMD風(fēng)振控制最優(yōu)性能研究
    復(fù)合屋面板鋼構(gòu)體系風(fēng)振特性試驗(yàn)
    在线精品无人区一区二区三 | 在线观看美女被高潮喷水网站| 丰满迷人的少妇在线观看| 国产视频首页在线观看| 国产精品福利在线免费观看| 少妇人妻久久综合中文| 精品午夜福利在线看| 免费av中文字幕在线| 免费av中文字幕在线| 国产av国产精品国产| 亚洲综合精品二区| 99热全是精品| 蜜桃久久精品国产亚洲av| 亚洲人成网站高清观看| 国产成人午夜福利电影在线观看| 美女脱内裤让男人舔精品视频| 日韩欧美 国产精品| 亚洲激情五月婷婷啪啪| 亚洲av男天堂| 青春草亚洲视频在线观看| 男男h啪啪无遮挡| 99久国产av精品国产电影| 国产精品嫩草影院av在线观看| 一级黄片播放器| 久热久热在线精品观看| 赤兔流量卡办理| 九草在线视频观看| 男女边摸边吃奶| 精品人妻视频免费看| 97超视频在线观看视频| 少妇裸体淫交视频免费看高清| 2021少妇久久久久久久久久久| 国产精品一二三区在线看| 免费黄网站久久成人精品| 久久久a久久爽久久v久久| 亚洲精华国产精华液的使用体验| av福利片在线观看| 熟女电影av网| 欧美亚洲 丝袜 人妻 在线| 国产伦精品一区二区三区视频9| 高清av免费在线| 久久久久久久精品精品| 男人和女人高潮做爰伦理| 18禁动态无遮挡网站| 新久久久久国产一级毛片| 国产av一区二区精品久久 | 国产永久视频网站| 全区人妻精品视频| 日韩国内少妇激情av| 人妻一区二区av| 老熟女久久久| av免费在线看不卡| 日韩欧美 国产精品| 亚洲av二区三区四区| 国产精品久久久久久精品电影小说 | 少妇人妻 视频| 午夜福利在线在线| 国产男女超爽视频在线观看| 久久毛片免费看一区二区三区| 99视频精品全部免费 在线| 免费黄频网站在线观看国产| 男人和女人高潮做爰伦理| 国产黄色视频一区二区在线观看| 国产伦精品一区二区三区四那| 岛国毛片在线播放| 插逼视频在线观看| 街头女战士在线观看网站| 亚洲伊人久久精品综合| 一二三四中文在线观看免费高清| 99久久精品国产国产毛片| 蜜桃亚洲精品一区二区三区| 久久精品熟女亚洲av麻豆精品| 国产探花极品一区二区| 男女边摸边吃奶| 一本色道久久久久久精品综合| 成人18禁高潮啪啪吃奶动态图 | 国产精品99久久99久久久不卡 | 99久久精品热视频| 青青草视频在线视频观看| 一本一本综合久久| 精品亚洲成a人片在线观看 | 91久久精品国产一区二区三区| 久久国产亚洲av麻豆专区| 国产老妇伦熟女老妇高清| 亚洲av不卡在线观看| 久久婷婷青草| 免费人妻精品一区二区三区视频| 高清av免费在线| 精品久久久久久久末码| 国产乱来视频区| 精品人妻一区二区三区麻豆| 不卡视频在线观看欧美| 国产亚洲5aaaaa淫片| 99久久精品热视频| 久久影院123| 亚洲婷婷狠狠爱综合网| a级毛色黄片| 夜夜爽夜夜爽视频| av天堂中文字幕网| 国产伦精品一区二区三区视频9| 青春草亚洲视频在线观看| 久久久久精品性色| 日本黄色片子视频| 人妻 亚洲 视频| 国产精品一区二区在线观看99| 2022亚洲国产成人精品| 国产亚洲午夜精品一区二区久久| 亚洲国产最新在线播放| 美女中出高潮动态图| 看十八女毛片水多多多| 欧美日韩精品成人综合77777| 男女国产视频网站| 亚洲综合精品二区| 老熟女久久久| 成人亚洲精品一区在线观看 | 亚洲欧洲国产日韩| 亚洲精品国产色婷婷电影| 国语对白做爰xxxⅹ性视频网站| 国产69精品久久久久777片| 日韩一区二区视频免费看| 国产有黄有色有爽视频| 丰满人妻一区二区三区视频av| 国产亚洲91精品色在线| 黄色视频在线播放观看不卡| 赤兔流量卡办理| 国产精品国产三级专区第一集| 99热网站在线观看| 欧美成人a在线观看| 在线免费观看不下载黄p国产| 免费黄频网站在线观看国产| 看十八女毛片水多多多| 777米奇影视久久| 一本久久精品| 蜜桃在线观看..| 99视频精品全部免费 在线| 精品久久久久久久久亚洲| 青春草亚洲视频在线观看| 欧美xxxx黑人xx丫x性爽| 国产精品麻豆人妻色哟哟久久| 九草在线视频观看| 成人亚洲欧美一区二区av| xxx大片免费视频| av在线播放精品| 亚洲四区av| 久久这里有精品视频免费| 国产有黄有色有爽视频| 丰满人妻一区二区三区视频av| 水蜜桃什么品种好| 午夜视频国产福利| 男女啪啪激烈高潮av片| 国国产精品蜜臀av免费| 日韩 亚洲 欧美在线| 又黄又爽又刺激的免费视频.| 网址你懂的国产日韩在线| 免费高清在线观看视频在线观看| 舔av片在线| 狂野欧美白嫩少妇大欣赏| 99热这里只有是精品50| 久久人妻熟女aⅴ| 日韩中字成人| 综合色丁香网| 简卡轻食公司| av在线老鸭窝| 亚洲电影在线观看av| 亚洲欧美成人综合另类久久久| 嫩草影院新地址| 99热全是精品| av在线app专区| 亚洲人成网站在线播| 乱系列少妇在线播放| 久久精品夜色国产| 美女xxoo啪啪120秒动态图| 黑人高潮一二区| 观看免费一级毛片| 直男gayav资源| 妹子高潮喷水视频| 欧美丝袜亚洲另类| 精品人妻熟女av久视频| 亚洲欧美中文字幕日韩二区| 热99国产精品久久久久久7| 中文资源天堂在线| 99视频精品全部免费 在线| 久久久亚洲精品成人影院| 2022亚洲国产成人精品| 偷拍熟女少妇极品色| 麻豆乱淫一区二区| videos熟女内射| 亚洲av日韩在线播放| 中国三级夫妇交换| 国产免费一区二区三区四区乱码| 国产 一区 欧美 日韩| 免费av不卡在线播放| 亚洲成人中文字幕在线播放| 少妇裸体淫交视频免费看高清| 啦啦啦视频在线资源免费观看| 免费在线观看成人毛片| 最近最新中文字幕大全电影3| 国产在线视频一区二区| 一区在线观看完整版| 亚洲内射少妇av| 观看美女的网站| 在线免费观看不下载黄p国产| 亚洲av综合色区一区| 成年女人在线观看亚洲视频| 国产精品国产三级专区第一集| 免费久久久久久久精品成人欧美视频 | 免费人妻精品一区二区三区视频| 国产精品一区二区在线不卡| 国产精品99久久99久久久不卡 | 大话2 男鬼变身卡| 国产亚洲午夜精品一区二区久久| 青春草国产在线视频| 91狼人影院| 国产精品福利在线免费观看| 国产成人免费无遮挡视频| 蜜桃久久精品国产亚洲av| 九草在线视频观看| 欧美三级亚洲精品| 男女边吃奶边做爰视频| 精品久久久久久久末码| 久久这里有精品视频免费| 日韩人妻高清精品专区| 男女边吃奶边做爰视频| 97精品久久久久久久久久精品| 久久久欧美国产精品| 交换朋友夫妻互换小说| 青青草视频在线视频观看| 美女视频免费永久观看网站| 岛国毛片在线播放| 久久久久视频综合| 夜夜骑夜夜射夜夜干| 老女人水多毛片| 免费观看在线日韩| 26uuu在线亚洲综合色| 18禁在线播放成人免费| 美女内射精品一级片tv| 免费观看无遮挡的男女| av在线播放精品| 午夜激情久久久久久久| 亚洲欧美日韩卡通动漫| 久久午夜福利片| 日韩强制内射视频| 欧美精品一区二区免费开放| 插逼视频在线观看| 国产免费福利视频在线观看| 亚洲精品自拍成人| 啦啦啦在线观看免费高清www| 99久国产av精品国产电影| 男人添女人高潮全过程视频| 午夜福利在线观看免费完整高清在| 国产精品免费大片| 秋霞在线观看毛片| 日韩,欧美,国产一区二区三区| 国内揄拍国产精品人妻在线| 日本av手机在线免费观看| 老熟女久久久| 美女脱内裤让男人舔精品视频| 777米奇影视久久| 国产精品99久久99久久久不卡 | 久久久亚洲精品成人影院| 99九九线精品视频在线观看视频| 久久人人爽人人片av| 一个人免费看片子| 夜夜爽夜夜爽视频| 老师上课跳d突然被开到最大视频| 99re6热这里在线精品视频| 天堂中文最新版在线下载| 又大又黄又爽视频免费| 免费不卡的大黄色大毛片视频在线观看| 亚洲国产欧美人成| 欧美成人a在线观看| 狂野欧美白嫩少妇大欣赏| 久久精品国产亚洲av天美| 蜜桃亚洲精品一区二区三区| 精品视频人人做人人爽| 欧美激情极品国产一区二区三区 | av女优亚洲男人天堂| 超碰97精品在线观看| 久久精品国产亚洲网站| 国产乱人偷精品视频| 欧美变态另类bdsm刘玥| 女的被弄到高潮叫床怎么办| 国产在线一区二区三区精| 人人妻人人看人人澡| 国产精品一区二区三区四区免费观看| 成人毛片60女人毛片免费| 26uuu在线亚洲综合色| 日韩三级伦理在线观看| 日日摸夜夜添夜夜爱| av女优亚洲男人天堂| 久久久久久人妻| 欧美极品一区二区三区四区| 国产av一区二区精品久久 | 精品一品国产午夜福利视频| 夜夜爽夜夜爽视频| 大片免费播放器 马上看| 亚洲精品日本国产第一区| 伊人久久国产一区二区| 色5月婷婷丁香| 久久女婷五月综合色啪小说| 人妻 亚洲 视频| 免费观看在线日韩| 日韩av不卡免费在线播放| 精品国产三级普通话版| 3wmmmm亚洲av在线观看| 国产精品一区二区在线不卡| 青青草视频在线视频观看| 亚洲四区av| 在线亚洲精品国产二区图片欧美 | 秋霞在线观看毛片| 欧美极品一区二区三区四区| 日韩成人伦理影院| 日本黄大片高清| 精华霜和精华液先用哪个| 黑人猛操日本美女一级片| 美女cb高潮喷水在线观看| 狠狠精品人妻久久久久久综合| 看免费成人av毛片| 在线观看免费日韩欧美大片 | 日韩中字成人| 美女国产视频在线观看| 国产欧美另类精品又又久久亚洲欧美| 亚洲,一卡二卡三卡| 亚洲精品第二区| 国产高清三级在线| 久久久久久久久久久免费av| 欧美三级亚洲精品| 亚洲伊人久久精品综合| 精品少妇久久久久久888优播| 国产乱来视频区| 久热这里只有精品99| 国产淫语在线视频| 欧美三级亚洲精品| 久久久久精品久久久久真实原创| 欧美国产精品一级二级三级 | 精华霜和精华液先用哪个| av免费在线看不卡| 国模一区二区三区四区视频| 国产精品国产三级国产专区5o| 亚洲丝袜综合中文字幕| 成人免费观看视频高清| 熟女电影av网| 一区二区av电影网| 久久6这里有精品| 51国产日韩欧美| 秋霞伦理黄片| 国产免费一级a男人的天堂| 亚洲欧美成人综合另类久久久| 一级片'在线观看视频| 性色avwww在线观看| 亚洲国产精品国产精品| av福利片在线观看| 五月伊人婷婷丁香| 美女xxoo啪啪120秒动态图| 女性被躁到高潮视频| 少妇精品久久久久久久| 国产高清不卡午夜福利| 看非洲黑人一级黄片| 久久99热这里只频精品6学生| 国产精品一区www在线观看| 中文字幕久久专区| 在线看a的网站| 日日撸夜夜添| 少妇裸体淫交视频免费看高清| 狂野欧美白嫩少妇大欣赏| 亚洲av成人精品一区久久| 国产探花极品一区二区| 蜜臀久久99精品久久宅男| 啦啦啦在线观看免费高清www| 18禁裸乳无遮挡免费网站照片| 久久精品国产a三级三级三级| 最近最新中文字幕大全电影3| 99re6热这里在线精品视频| av在线蜜桃| 在线观看三级黄色| 亚洲av免费高清在线观看| 日日啪夜夜爽| 亚洲精品乱码久久久久久按摩| 91狼人影院| 韩国av在线不卡| 夫妻性生交免费视频一级片| 97超视频在线观看视频| a级毛色黄片| 久久人妻熟女aⅴ| 精品一区在线观看国产| 欧美高清性xxxxhd video| 日韩人妻高清精品专区| 大片免费播放器 马上看| 97超碰精品成人国产| 九九久久精品国产亚洲av麻豆| 人妻系列 视频| 熟女电影av网| 久久精品久久精品一区二区三区| 国产 精品1| 51国产日韩欧美| 不卡视频在线观看欧美| 国产黄频视频在线观看| 亚洲四区av| av.在线天堂| 国产在视频线精品| 女性生殖器流出的白浆| 97超视频在线观看视频| 欧美高清性xxxxhd video| 男女边摸边吃奶| 精品人妻偷拍中文字幕| 男女国产视频网站| 国产在线一区二区三区精| 国产精品福利在线免费观看| 国产亚洲精品久久久com| 欧美激情国产日韩精品一区| 欧美精品国产亚洲| 欧美丝袜亚洲另类| 免费av不卡在线播放| 亚洲av电影在线观看一区二区三区| 国产日韩欧美在线精品| 深夜a级毛片| 狂野欧美白嫩少妇大欣赏| 大陆偷拍与自拍| 亚洲精品国产色婷婷电影| 国产成人a区在线观看| 久久久久精品久久久久真实原创| 亚洲精品成人av观看孕妇| 日本-黄色视频高清免费观看| 亚洲,一卡二卡三卡| 欧美高清成人免费视频www| 一本—道久久a久久精品蜜桃钙片| 亚洲精品国产av蜜桃| 这个男人来自地球电影免费观看 | 色视频www国产| 插逼视频在线观看| 少妇猛男粗大的猛烈进出视频| 国产伦在线观看视频一区| 一边亲一边摸免费视频| 久久精品国产a三级三级三级| freevideosex欧美| 日本欧美国产在线视频| 成人高潮视频无遮挡免费网站| 97在线视频观看| 一个人免费看片子| 视频区图区小说| 色视频www国产| av网站免费在线观看视频| 国产黄片视频在线免费观看| 蜜桃在线观看..| 又粗又硬又长又爽又黄的视频| 亚洲图色成人| 日本与韩国留学比较| 国产乱人偷精品视频| 成年女人在线观看亚洲视频| 亚洲伊人久久精品综合| 婷婷色av中文字幕| 一区二区三区精品91| 国产精品久久久久久精品电影小说 | av播播在线观看一区| 亚洲欧美中文字幕日韩二区| 九九在线视频观看精品| av在线播放精品| 色5月婷婷丁香| 国国产精品蜜臀av免费| 成人综合一区亚洲| 最近最新中文字幕免费大全7| 中文字幕免费在线视频6| 亚洲欧美日韩无卡精品| 熟女人妻精品中文字幕| 一级毛片黄色毛片免费观看视频| 亚洲色图综合在线观看| 欧美老熟妇乱子伦牲交| 欧美国产精品一级二级三级 | 国产毛片在线视频| 2018国产大陆天天弄谢| 国产大屁股一区二区在线视频| 少妇人妻 视频| 国产一级毛片在线| 国产精品成人在线| 国产乱来视频区| 国内少妇人妻偷人精品xxx网站| 一区二区三区乱码不卡18| 国产白丝娇喘喷水9色精品| 欧美一区二区亚洲| 欧美日本视频| 国产精品三级大全| 一级爰片在线观看| 色吧在线观看| 尤物成人国产欧美一区二区三区| 深夜a级毛片| 亚洲精品第二区| 男的添女的下面高潮视频| 六月丁香七月| 国产成人精品福利久久| 91在线精品国自产拍蜜月| 午夜免费男女啪啪视频观看| 国产欧美另类精品又又久久亚洲欧美| 极品教师在线视频| 国产极品天堂在线| 日韩 亚洲 欧美在线| 交换朋友夫妻互换小说| 亚洲av免费高清在线观看| 欧美最新免费一区二区三区| 激情 狠狠 欧美| 国产伦理片在线播放av一区| a级毛片免费高清观看在线播放| 久热这里只有精品99| 国产精品三级大全| 另类亚洲欧美激情| 免费av中文字幕在线| 亚洲丝袜综合中文字幕| 又粗又硬又长又爽又黄的视频| 久久久久久久亚洲中文字幕| 美女主播在线视频| 啦啦啦啦在线视频资源| 中文欧美无线码| 国产无遮挡羞羞视频在线观看| 久久99蜜桃精品久久| 韩国高清视频一区二区三区| 成人毛片60女人毛片免费| 在线观看免费高清a一片| 免费看光身美女| 深爱激情五月婷婷| 久久久精品免费免费高清| 欧美国产精品一级二级三级 | 欧美97在线视频| 日本猛色少妇xxxxx猛交久久| videos熟女内射| av不卡在线播放| 国产白丝娇喘喷水9色精品| 欧美zozozo另类| 人妻制服诱惑在线中文字幕| 观看美女的网站| 伊人久久精品亚洲午夜| 91在线精品国自产拍蜜月| 日本黄色日本黄色录像| 精品一区在线观看国产| 成人综合一区亚洲| 亚洲内射少妇av| 内地一区二区视频在线| 欧美精品一区二区大全| 国产亚洲精品久久久com| 性高湖久久久久久久久免费观看| 国产精品国产三级国产av玫瑰| 亚洲精品456在线播放app| 亚洲精品,欧美精品| 黄色视频在线播放观看不卡| 欧美xxⅹ黑人| 热99国产精品久久久久久7| 大话2 男鬼变身卡| 下体分泌物呈黄色| 久久久亚洲精品成人影院| 最近的中文字幕免费完整| 亚洲av国产av综合av卡| 亚洲精品第二区| 日产精品乱码卡一卡2卡三| 黄色视频在线播放观看不卡| 国内少妇人妻偷人精品xxx网站| 国产精品99久久99久久久不卡 | 色哟哟·www| 精品一品国产午夜福利视频| 成人特级av手机在线观看| 另类亚洲欧美激情| 欧美一区二区亚洲| 亚洲久久久国产精品| 亚洲色图综合在线观看| 人人妻人人澡人人爽人人夜夜| 成人无遮挡网站| 插逼视频在线观看| 麻豆成人午夜福利视频| 边亲边吃奶的免费视频| 最近2019中文字幕mv第一页| 日产精品乱码卡一卡2卡三| 国产无遮挡羞羞视频在线观看| 观看av在线不卡| 丝瓜视频免费看黄片| 亚洲伊人久久精品综合| 最近最新中文字幕大全电影3| 国产精品免费大片| 在线播放无遮挡| 欧美一区二区亚洲| 女人久久www免费人成看片| 精品久久久精品久久久| 久久久久久久久久人人人人人人| 亚洲精品日本国产第一区| 欧美亚洲 丝袜 人妻 在线| 久久久久久久大尺度免费视频| 搡女人真爽免费视频火全软件| 内地一区二区视频在线| 亚洲av.av天堂| 国产日韩欧美在线精品| av一本久久久久| 国产探花极品一区二区| 国产成人免费无遮挡视频| 美女脱内裤让男人舔精品视频| 久久久色成人| 国产精品国产三级专区第一集| 国产精品久久久久久久久免| 国产精品成人在线| 久久久久久久国产电影| 97超视频在线观看视频| 久久婷婷青草| 精品午夜福利在线看| 国产又色又爽无遮挡免| 国产黄片视频在线免费观看| 国产高清不卡午夜福利| 亚洲国产精品成人久久小说| 亚洲精品自拍成人| 欧美成人精品欧美一级黄| 九九久久精品国产亚洲av麻豆| 又黄又爽又刺激的免费视频.| 激情五月婷婷亚洲| 看免费成人av毛片| 国产亚洲精品久久久com| 亚洲精品自拍成人| 三级国产精品片| 婷婷色综合大香蕉| av在线app专区| 大片电影免费在线观看免费| 欧美老熟妇乱子伦牲交|