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

    強(qiáng)風(fēng)作用下超大型冷卻塔結(jié)構(gòu)失效準(zhǔn)則與強(qiáng)健性分析

    2021-09-08 05:59:26王飛天柯世堂王曉海
    振動(dòng)工程學(xué)報(bào) 2021年4期

    王飛天 柯世堂 王曉海

    摘要: 為分析超大型冷卻塔在強(qiáng)風(fēng)作用下的結(jié)構(gòu)失效準(zhǔn)則與強(qiáng)健性,以超規(guī)范高度限值的世界最高220 m超大型冷卻塔為例,基于ANSYS/LS?DYNA平臺(tái)建立結(jié)構(gòu)三維有限元模型,對(duì)其風(fēng)致倒塌全過程進(jìn)行數(shù)值仿真分析。研究了典型風(fēng)速下位移響應(yīng)與塔筒內(nèi)力分布特性,得出最不利響應(yīng)發(fā)生的位置、數(shù)值及分布規(guī)律,探討了強(qiáng)風(fēng)下塔筒結(jié)構(gòu)“穩(wěn)定?失穩(wěn)?倒塌”全過程演化機(jī)理,并提出了強(qiáng)風(fēng)作用下冷卻塔的結(jié)構(gòu)失效準(zhǔn)則。研究發(fā)現(xiàn):強(qiáng)風(fēng)作用下超大型冷卻塔結(jié)構(gòu)強(qiáng)健性最薄弱部位為迎風(fēng)面塔筒喉部區(qū)域,與上風(fēng)向夾角為-70°和70°子午向出現(xiàn)“褶皺變形”并導(dǎo)致塔頂局部結(jié)構(gòu)失效破壞,進(jìn)而引起塔筒整體結(jié)構(gòu)傾覆倒塌;塔筒應(yīng)變能密度呈現(xiàn)指數(shù)增長(zhǎng)趨勢(shì),以喉部相對(duì)水平位移與喉部直徑百分比S>1%、指數(shù)應(yīng)變能密度和值破壞系數(shù)K≥0.3作為結(jié)構(gòu)失效判斷準(zhǔn)則可較好評(píng)價(jià)強(qiáng)健性指標(biāo)。

    關(guān)鍵詞: 超大型冷卻塔; 強(qiáng)風(fēng)作用; 連續(xù)倒塌; 失效準(zhǔn)則; 強(qiáng)健性分析

    引 言

    結(jié)構(gòu)強(qiáng)健性[1]是由結(jié)構(gòu)連續(xù)倒塌極限狀態(tài)發(fā)展而來,其基本含義是結(jié)構(gòu)在施工和使用階段不發(fā)生與初始破壞原因不相稱的破壞。結(jié)合強(qiáng)風(fēng)破壞屬性和超大型冷卻塔自身特征,冷卻塔結(jié)構(gòu)強(qiáng)健性是指在強(qiáng)風(fēng)環(huán)境下塔筒局部損傷失效不應(yīng)導(dǎo)致整體結(jié)構(gòu)倒塌破壞或屈曲失穩(wěn)。1965年英國(guó)渡橋電廠和1973年蘇格蘭阿德爾曼電廠大型冷卻塔的風(fēng)毀事件[2?3]均是由局部風(fēng)荷載過大使得塔筒局部受拉損傷,進(jìn)而引發(fā)整個(gè)塔筒結(jié)構(gòu)強(qiáng)度破壞和屈曲失穩(wěn)倒塌。事實(shí)上,此類由局部損壞引發(fā)的整體風(fēng)毀破壞均屬于典型的結(jié)構(gòu)強(qiáng)健性破壞形式。然而,目前國(guó)內(nèi)外冷卻塔設(shè)計(jì)規(guī)范和已有研究文獻(xiàn)鮮有引入強(qiáng)健性設(shè)計(jì)理念,尤其是在強(qiáng)臺(tái)風(fēng)、雷暴風(fēng)等極端風(fēng)災(zāi)作用下超大型冷卻塔的破壞極限承載力狀態(tài)設(shè)計(jì)中。

    針對(duì)大型冷卻塔抗風(fēng)穩(wěn)定性研究,目前國(guó)內(nèi)外學(xué)者的研究主要集中在靜風(fēng)失穩(wěn)[4]、整體與局部屈曲失穩(wěn)[5?8]、施工全過程穩(wěn)定性能[9]、幾何缺陷導(dǎo)致的失穩(wěn)[10]、不同氣動(dòng)措施對(duì)冷卻塔穩(wěn)定性影響[11?12]等方面,這些研究成果很好地解決了超大型冷卻塔在施工與運(yùn)營(yíng)期內(nèi)塔筒抗風(fēng)穩(wěn)定性能的問題。針對(duì)結(jié)構(gòu)強(qiáng)健性研究,文獻(xiàn)[13]為提升鋼管混凝土拱橋體系承載能力,從概念設(shè)計(jì)角度提出了拱橋橋面系和吊桿的強(qiáng)健性設(shè)計(jì)方案;文獻(xiàn)[14?16]針對(duì)不同形式的橋梁結(jié)構(gòu)特點(diǎn),基于構(gòu)件拆除計(jì)算方法,確定橋梁的主要構(gòu)件并系統(tǒng)研究了主要構(gòu)件破壞對(duì)其余結(jié)構(gòu)的靜力和動(dòng)力響應(yīng),提出了增強(qiáng)橋梁強(qiáng)健性的措施;文獻(xiàn)[17]制定和實(shí)施了結(jié)構(gòu)系統(tǒng)與基礎(chǔ)設(shè)施網(wǎng)絡(luò)性能標(biāo)準(zhǔn),建立了結(jié)構(gòu)強(qiáng)健性概念及網(wǎng)絡(luò)性能指標(biāo);文獻(xiàn)[18]針對(duì)中、下承式拱橋技術(shù)狀況評(píng)定中的不足,提出考慮懸吊橋面系強(qiáng)健性的技術(shù)狀況評(píng)定方法更能反映橋梁的實(shí)際狀況,適應(yīng)維護(hù)管理需求。盡管結(jié)構(gòu)強(qiáng)健性研究已有開展,但是國(guó)內(nèi)外學(xué)者關(guān)于結(jié)構(gòu)強(qiáng)健性的研究多針對(duì)于橋梁結(jié)構(gòu)方向,而對(duì)于大型冷卻塔結(jié)構(gòu)抗風(fēng)強(qiáng)健性分析處于空白階段。

    鑒于此,本文以已建成的世界最高的220 m山西潞安電廠大型冷卻塔為工程背景,基于ANSYS/LS?DYNA平臺(tái)對(duì)超大型冷卻塔進(jìn)行精細(xì)化建模和風(fēng)致倒塌全過程數(shù)值仿真。在此基礎(chǔ)上,探討了強(qiáng)風(fēng)作用下冷卻塔結(jié)構(gòu)抗風(fēng)穩(wěn)定性能,分析了塔筒在強(qiáng)風(fēng)作用下的變形、位移與內(nèi)力響應(yīng)規(guī)律,揭示塔筒結(jié)構(gòu)“穩(wěn)定?失穩(wěn)?倒塌”全過程的強(qiáng)健性薄弱部位和連續(xù)性倒塌演化機(jī)理,最終提出超大型冷卻塔在強(qiáng)風(fēng)激勵(lì)下的結(jié)構(gòu)失效準(zhǔn)則。

    1 工程概況和有限元建模

    1.1 工程概況

    山西潞安電廠超大型冷卻塔是目前世界已建的最高的冷卻塔,主體結(jié)構(gòu)由地基、環(huán)基、支柱和塔筒構(gòu)成,塔高220.0 m,進(jìn)風(fēng)口標(biāo)高30.75 m,喉部標(biāo)高165.0 m,底部直徑185.0 m,喉部直徑123.0 m,頂部出風(fēng)口直徑128.1 m,塔筒殼體為指數(shù)變厚,最小與最大壁厚為0.34和1.85 m,由64榀X支柱與環(huán)板基礎(chǔ)連接,X型支柱采用矩形截面(尺寸為1.70 m× 1.00 m)鋼筋混凝土結(jié)構(gòu),環(huán)板基礎(chǔ)為寬10.5 m,高2.20 m現(xiàn)澆鋼筋混凝土結(jié)構(gòu)。主要結(jié)構(gòu)尺寸參數(shù)如表1所示。

    1.2 有限元建模

    基于ANSYS/LS?DYNA平臺(tái)建立該超大型冷卻塔的三維有限元足尺模型,其中塔筒及頂部剛性環(huán)采用殼單元(Shell163),環(huán)向和子午向分別劃分為256和128個(gè)單元;塔筒下部64榀X型支柱均采用空間梁?jiǎn)卧˙eam161),X支柱上端與塔筒底部采用節(jié)點(diǎn)剛域耦合方式連接,支柱下端固支作為模型邊界條件,如圖1所示。

    考慮到冷卻塔塔筒環(huán)向與縱向配筋率較大,鋼筋的強(qiáng)大拉接力使塔壁形成整體受力良好的聯(lián)合體,此種聯(lián)合體可將局部荷載“稀釋”,在本構(gòu)關(guān)系上塑性特點(diǎn)突出。因此,基于材料的塑性隨動(dòng)理論,塔壁與支柱均采用塑性隨動(dòng)材料模型(*MAT_PLASTIC_KINEMATIC),在硬化模式上采用完全隨動(dòng)硬化模式,即β=0;接觸類型采用適應(yīng)性較強(qiáng)的自動(dòng)單面接觸, 材料摩擦系數(shù)設(shè)定為0.25。具體參數(shù)如表2所示。

    1.3 動(dòng)力特性分析

    圖2給出ANSYS/LS?DYNA與文獻(xiàn)[19?21]中已有的超大型冷卻塔模型前100階自振頻率隨振型階數(shù)變化對(duì)比曲線。前者采用顯式單元Shell163和Beam161建立冷卻塔模型,后者采用隱式單元Shell63和Beam188建立冷卻塔模型。對(duì)比發(fā)現(xiàn)兩種建模方法下自振頻率趨勢(shì)相同,大致呈線性增長(zhǎng),0?60階范圍內(nèi)兩種計(jì)算結(jié)果基本一致,隨著階數(shù)的增長(zhǎng)二者的差異有所增大,第100階自振頻率僅相差4%;自振頻率分布密集,基頻僅為0.58 Hz,前50階振型頻率集中分布在0.58?1.56 Hz之間,說明ANSYS/LS?DYNA超大型冷卻塔有限元模型具有一定的有效性。

    圖3給出冷卻塔結(jié)構(gòu)前100階典型模態(tài)振型圖,分析發(fā)現(xiàn)塔筒振型在環(huán)向和子午向差異較大,隨著階數(shù)增大,環(huán)向和子午向諧波數(shù)顯著增加;低階振型子午向呈現(xiàn)1?3個(gè)諧波,環(huán)向呈現(xiàn)6?12個(gè)諧波,且多以塔頂振動(dòng)為主。

    2 風(fēng)致響應(yīng)分析

    該超大型冷卻塔平均風(fēng)壓系數(shù)按光滑塔取值,地貌類型為B類,設(shè)計(jì)基本風(fēng)速為23.7 m/s。以-0.5的內(nèi)壓系數(shù)考慮塔筒內(nèi)吸力效應(yīng)[22],以10 m高度處25 m/s的初始風(fēng)速為基礎(chǔ)進(jìn)行逐級(jí)加載,加載風(fēng)速步長(zhǎng)為5 m/s,時(shí)程分析步長(zhǎng)為0.5 s。大型冷卻塔在強(qiáng)風(fēng)作用下,其結(jié)構(gòu)的破壞形式是類似屈曲模態(tài)下的材料破壞,可將風(fēng)荷載視為一種等效靜風(fēng)荷載。為了降低加載過程中的動(dòng)力效應(yīng),防止出現(xiàn)類似沖擊荷載的效應(yīng),將平均風(fēng)壓在20 s內(nèi)線性遞增形成時(shí)程分析輸入所需的風(fēng)荷載,再以點(diǎn)荷載加載到塔筒相應(yīng)的加載點(diǎn)。

    2.1 位移響應(yīng)

    圖4分別給出25 m/s,50 m/s和75 m/s三種典型風(fēng)速下塔筒順風(fēng)向位移三維圖和最大位移時(shí)程曲線。分析可知三種典型風(fēng)速只引起殼體的結(jié)構(gòu)變形,冷卻塔未發(fā)生倒塌;變形主要發(fā)生在與上風(fēng)方向夾角-70°?70°的喉部區(qū)域,塔筒背風(fēng)面變形較小;25和50 m/s風(fēng)速下的塔筒最大位移均出現(xiàn)在迎風(fēng)面喉部,隨著殼體厚度的增加,喉部上、下變形逐漸減小;75 m/s風(fēng)速下最大位移位置有所下降但依然出現(xiàn)在迎風(fēng)面塔筒中上部;冷卻塔塔筒在風(fēng)荷載激勵(lì)的初始階段會(huì)出現(xiàn)外加瞬態(tài)激勵(lì)引起的瞬態(tài)振動(dòng),最大位移?時(shí)程曲線表現(xiàn)出明顯的波動(dòng)狀態(tài)。風(fēng)荷載作用5 s后,隨著系統(tǒng)阻尼的影響,瞬態(tài)振動(dòng)逐漸衰弱,最大位移?時(shí)程曲線趨于平穩(wěn),且最大位移分別為0.08,0.3和2.0 m,表明超大型冷卻塔塔筒在風(fēng)荷載的作用下處于穩(wěn)定狀態(tài)。

    2.2 內(nèi)力響應(yīng)

    圖5和6分別給出了25, 50和75 m/s三種典型風(fēng)速下塔筒子午向應(yīng)力三維等勢(shì)線圖和0°子午向應(yīng)力隨高度的變化曲線,圖中正值表示拉應(yīng)力,分析可知不同風(fēng)速下的子午向應(yīng)力三維等勢(shì)線圖具有對(duì)稱性,分布模式相同,變化幅度不同,塔筒子午向應(yīng)力隨風(fēng)速的增加而增大。0°子午向應(yīng)力主要集中于塔體中下部,且三種風(fēng)速下應(yīng)力均未超過屈服應(yīng)力,塔筒仍處于直立穩(wěn)定狀態(tài);沿子午向高度,塔筒底部至中部應(yīng)力較大,頂部應(yīng)力較小,接近于零值;沿環(huán)向角度,0°?45°子午線應(yīng)力為正,45°?110°應(yīng)力為負(fù),背風(fēng)面110°?180°子午向應(yīng)力為較小的正值。

    3 強(qiáng)健性分析與失效準(zhǔn)則

    3.1 強(qiáng)健性分析

    3.1.1 臨界倒塌風(fēng)速

    首先,給出超大型冷卻塔的結(jié)構(gòu)強(qiáng)健性定義:冷卻塔結(jié)構(gòu)在承受強(qiáng)(臺(tái))風(fēng)、雷暴風(fēng)或龍卷風(fēng)等極端事件時(shí),局部構(gòu)件的損傷失效或破壞不應(yīng)導(dǎo)致塔筒整體出現(xiàn)影響使用的連續(xù)倒塌或屈曲失穩(wěn)。然后,確定冷卻塔結(jié)構(gòu)強(qiáng)健性分析的內(nèi)容,包括強(qiáng)風(fēng)作用下塔筒倒塌機(jī)理、塔筒結(jié)構(gòu)的臨界倒塌風(fēng)速、塔筒連續(xù)倒塌全過程及強(qiáng)健性薄弱位置。

    基于增量動(dòng)力分析(IDA)方法對(duì)塔筒結(jié)構(gòu)進(jìn)行非線性分析,圖7給出不同風(fēng)速下0°子午線塔頂、喉部和塔底部位的順風(fēng)向位移變化曲線。分析可知風(fēng)速低于70 m/s時(shí),三個(gè)不同部位的順風(fēng)向位移隨風(fēng)速基本呈線性關(guān)系,70?80 m/s風(fēng)速范圍內(nèi)順風(fēng)向位移隨風(fēng)速呈非線性關(guān)系,風(fēng)速為85 m/s時(shí)三個(gè)不同部位的順風(fēng)向位移出現(xiàn)4?10倍突增,并且在80 m/s風(fēng)速時(shí)三個(gè)不同部位的順風(fēng)向位移均出現(xiàn)明顯拐點(diǎn)(圖中U點(diǎn))。同時(shí),數(shù)值模擬結(jié)果顯示當(dāng)風(fēng)速大于80 m/s時(shí)超大型冷卻塔塔筒出現(xiàn)失穩(wěn)而整體倒塌,可判定超大型冷卻塔的臨界倒塌風(fēng)速為80 m/s。

    為驗(yàn)證冷卻塔數(shù)值模擬倒塌的有效性,以文獻(xiàn)[4]對(duì)英國(guó)研究者Der的冷卻塔高壓風(fēng)洞穩(wěn)定性試驗(yàn)的介紹和再現(xiàn)計(jì)算進(jìn)行對(duì)比。圖8給出Der試驗(yàn)?zāi)P褪Х€(wěn)形態(tài)與本文數(shù)值模擬失穩(wěn)形態(tài)對(duì)比圖??芍呔哂邢嗤氖Х€(wěn)形態(tài),環(huán)向變形均與塔筒風(fēng)壓分布曲線類似,并且Der風(fēng)洞試驗(yàn)?zāi)P偷脑褪Х€(wěn)風(fēng)速均達(dá)到87 m/s以上甚至達(dá)到百米級(jí)風(fēng)速,由此可知本文的冷卻塔數(shù)值模擬臨界倒塌風(fēng)速為80 m/s具有一定的有效性。

    3.1.2 連續(xù)性倒塌全過程

    圖9給出超大型冷卻塔塔筒連續(xù)性倒塌發(fā)展全過程。在強(qiáng)風(fēng)作用初始,塔筒中部最先產(chǎn)生變形,與上風(fēng)向夾角-70°?70°范圍內(nèi)塔筒呈現(xiàn)凹形,橫風(fēng)向兩側(cè)呈現(xiàn)凸形,環(huán)向變形與風(fēng)壓分布曲線一致;t=5.99 s至t=9.49 s塔筒順風(fēng)向凹陷程度增大,最大變形位置逐漸上升,內(nèi)吸力作用導(dǎo)致環(huán)向70°?110°和250°?290°范圍沿橫風(fēng)向凹陷;t=9.49 s時(shí)最大變形部位由塔筒中部上升至喉部區(qū)域,且與上風(fēng)向夾角-70°和70°子午向的塔筒中部出現(xiàn)“褶皺變形”,塔筒因剪切撕裂出現(xiàn)豎向裂縫,但冷卻塔整體仍保持直立狀態(tài);t=10.49 s至t=12.49 s喉部結(jié)構(gòu)單元失效,塔筒變形凹陷程度不斷增大,最大變形位置逐漸上升,最終引起塔頂局部構(gòu)件破壞,塔筒和支柱出現(xiàn)明顯的幾何大變形,冷卻塔整體進(jìn)入傾覆狀態(tài);t=17.49 s時(shí)塔筒傳力路徑發(fā)生嚴(yán)重破壞,超大型冷卻塔塔筒因無法繼續(xù)承受強(qiáng)風(fēng)作用而發(fā)生連續(xù)性倒塌。

    由此可以得知,在強(qiáng)風(fēng)作用下,塔筒喉部區(qū)域?yàn)槌笮屠鋮s塔的強(qiáng)健性最薄弱部位。喉部結(jié)構(gòu)局部破壞失效使冷卻塔整體承載力降低,導(dǎo)致大型冷卻塔整體從直立穩(wěn)定狀態(tài)迅速進(jìn)入傾覆倒塌狀態(tài)。

    3.2 失效準(zhǔn)則

    進(jìn)行超大型冷卻塔在強(qiáng)風(fēng)作用下的結(jié)構(gòu)強(qiáng)健性分析,必須判斷出在何種條件下塔筒結(jié)構(gòu)出現(xiàn)失效而破壞。目前,主要通過以下幾種分析準(zhǔn)則來判斷結(jié)構(gòu)的失效:變形失效準(zhǔn)則[23]、強(qiáng)度失效準(zhǔn)則、能量失效準(zhǔn)則[24]、變形和能量雙重失效準(zhǔn)則。本文從變形和能量?jī)蓚€(gè)角度提出超大型冷卻塔的結(jié)構(gòu)失效判別準(zhǔn)則。其中,變形失效準(zhǔn)則以塔筒喉部相對(duì)水平位移為目標(biāo),分析不同風(fēng)速下的喉部相對(duì)水平位移與喉部直徑的比例關(guān)系;能量失效準(zhǔn)則以塔筒整體所有單元的應(yīng)變能密度和值為目標(biāo),分析塔筒應(yīng)變能密度和值增長(zhǎng)速度在不同風(fēng)速下的發(fā)展規(guī)律。

    3.2.1 變形失效準(zhǔn)則

    大型冷卻塔的強(qiáng)健性薄弱位置為迎風(fēng)面塔筒喉部區(qū)域,塔筒喉部的相對(duì)變形直接影響冷卻塔的完整性。美國(guó)DOD2010規(guī)范[25]將結(jié)構(gòu)連續(xù)倒塌判定準(zhǔn)則定義為部分構(gòu)件發(fā)生初始破壞后,剩余結(jié)構(gòu)中的結(jié)構(gòu)構(gòu)件不能出現(xiàn)規(guī)定的極限變形。文獻(xiàn)[23]依據(jù)該規(guī)范對(duì)混凝土框架的連續(xù)倒塌失效準(zhǔn)則規(guī)定為梁兩端的相對(duì)豎向位移達(dá)到跨度的1/5,依據(jù)DOD2010規(guī)范及文獻(xiàn)[20]的研究思路分析超大型冷卻塔喉部相對(duì)水平位移與喉部直徑的比例關(guān)系。

    表3給出不同風(fēng)速、不同時(shí)間點(diǎn)塔筒喉部迎風(fēng)面與背風(fēng)面的相對(duì)水平位移(Dr)與喉部中面直徑(R)的百分比(S)。圖10給出不同風(fēng)速、不同時(shí)間點(diǎn)S?V0變化曲線。分析可知不同時(shí)間點(diǎn)的相對(duì)水平位移與風(fēng)速基本呈線性關(guān)系,塔筒處于彈性變形階段,冷卻塔整體處于小變形范圍。風(fēng)速低于75 m/s時(shí),喉部相對(duì)水平位移與喉部直徑的百分比均小于1%。當(dāng)風(fēng)速大于75 m/s時(shí),喉部相對(duì)水平位移與喉部直徑的百分比呈現(xiàn)急劇增加,表現(xiàn)出明顯的拐點(diǎn)(圖中P點(diǎn)),冷卻塔結(jié)構(gòu)的單元開始出現(xiàn)失效,強(qiáng)風(fēng)持續(xù)作用將導(dǎo)致塔筒出現(xiàn)連續(xù)倒塌。

    3.1.1節(jié)的數(shù)值模擬結(jié)果表明風(fēng)速達(dá)到80 m/s時(shí)超大型冷卻塔塔筒出現(xiàn)失穩(wěn)而整體倒塌,在此基礎(chǔ)上,以1%作為臨界點(diǎn)來判斷超大型冷卻塔塔筒是否出現(xiàn)失穩(wěn)而倒塌。即冷卻塔塔筒處于正常工作狀態(tài)時(shí),塔筒喉部迎風(fēng)面、背風(fēng)面的相對(duì)水平位移與喉部直徑的百分比S1%;當(dāng)S>1%時(shí),塔筒進(jìn)入結(jié)構(gòu)單元失穩(wěn)狀態(tài)。

    3.2.2 能量失效準(zhǔn)則

    考慮到能量可以較好表征地震、強(qiáng)風(fēng)作用下結(jié)構(gòu)的真實(shí)彈塑性性能,引入能量法研究強(qiáng)風(fēng)作用下超大型冷卻塔結(jié)構(gòu)倒塌破壞與失效準(zhǔn)則。依據(jù)文獻(xiàn)[26]對(duì)弦支穹頂結(jié)構(gòu)失效準(zhǔn)則的研究,分析超大型冷卻塔塔筒在不同風(fēng)速激勵(lì)下的塔筒應(yīng)變能密度和值Pd?;贚S?Prepost有限元軟件提取不同荷載域、不同時(shí)程下塔筒全部單元的應(yīng)變能及應(yīng)變能密度。

    表4列出數(shù)值模擬過程中t=20 s時(shí)不同風(fēng)速(V0)下冷卻塔模型所有單元的應(yīng)變能密度和值(Pd)。并采用MATLAB對(duì)不同風(fēng)速下所有單元的應(yīng)變能密度和值分布規(guī)律進(jìn)行擬合,得到擬合公式如下式所示

    分析可知數(shù)值模擬與擬合公式(2)的應(yīng)變能密度和值關(guān)系曲線基本吻合,擬合效果較好;超大型冷卻塔在強(qiáng)風(fēng)作用t=20 s時(shí)所有單元的應(yīng)變能密度隨著風(fēng)速的增大而增加;當(dāng)風(fēng)速低于80 m/s時(shí),冷卻塔在不同風(fēng)速作用下的應(yīng)變能密度和值Pd符合指數(shù)律增長(zhǎng)趨勢(shì),應(yīng)變能密度增長(zhǎng)速度較小;當(dāng)風(fēng)速增大至某一值,冷卻塔應(yīng)變能密度和值出現(xiàn)突增現(xiàn)象,表現(xiàn)出明顯的拐點(diǎn)(圖中T點(diǎn)),此時(shí),數(shù)值模擬的Pd?V0關(guān)系曲線斜率從0.30突變?yōu)?.11。

    基于Pd?V0關(guān)系曲線的分析,以Pd?V0關(guān)系曲線的斜率進(jìn)行度量,提出判別超大型冷卻塔結(jié)構(gòu)失效的指數(shù)應(yīng)變能密度和值破壞系數(shù)K,表達(dá)式如下式所示

    圖12給出不同風(fēng)速下的K?V0關(guān)系曲線??梢园l(fā)現(xiàn)在風(fēng)速達(dá)到80 m/s時(shí),破壞系數(shù)K出現(xiàn)明顯的拐點(diǎn)(圖中V點(diǎn)),3.1.1節(jié)的數(shù)值模擬結(jié)果表明風(fēng)速達(dá)到80 m/s時(shí)超大型冷卻塔塔筒出現(xiàn)失穩(wěn)而整體倒塌。因此,以拐點(diǎn)處破壞系數(shù)作為臨界點(diǎn)來判斷超大型冷卻塔塔筒是否出現(xiàn)失穩(wěn)而倒塌。即當(dāng)冷卻塔塔筒處于穩(wěn)定工作狀態(tài)時(shí),K位于0?0.3之間。K ≥0.3時(shí),超大型冷卻塔在強(qiáng)風(fēng)作用下進(jìn)入結(jié)構(gòu)單元失穩(wěn)狀態(tài)。

    4 結(jié) 論

    本文系統(tǒng)開展了強(qiáng)風(fēng)作用下超大型冷卻塔倒塌全過程數(shù)值模擬與破壞機(jī)理研究,內(nèi)容涉及動(dòng)力特性、風(fēng)致響應(yīng)、結(jié)構(gòu)強(qiáng)健性和失效判別準(zhǔn)則。主要結(jié)論如下:

    (1)基于ANSYS/LS?DYNA顯式結(jié)構(gòu)單元建立超大型冷卻塔有限元模型,與已有的超大型冷卻塔模型自振頻率結(jié)果基本一致,發(fā)展趨勢(shì)相同,最大差異僅為4%,驗(yàn)證了本文ANSYS/LS?DYNA冷卻塔模型的有效性。

    (2)最不利順風(fēng)向位移響應(yīng)位置出現(xiàn)在超大型冷卻塔迎風(fēng)面喉部區(qū)域,沿高度方向塔筒子午向應(yīng)力逐漸減小,不同風(fēng)速下0°子午向應(yīng)力集中分布于塔筒中部,塔頂應(yīng)力最小,接近于零值。

    (3)超大型冷卻塔臨界倒塌基本風(fēng)速達(dá)到80 m/s,結(jié)構(gòu)強(qiáng)健性最薄弱部位發(fā)生在喉部區(qū)域;喉部結(jié)構(gòu)單元局部破壞失效引發(fā)與上風(fēng)向夾角為-70°和70°子午向塔筒中部發(fā)生“褶皺變形”,冷卻塔從直立狀態(tài)進(jìn)入傾覆狀態(tài),塔頂局部結(jié)構(gòu)失效破壞激發(fā)超大型冷卻塔出現(xiàn)連續(xù)性倒塌。

    (4)超大型冷卻塔塔筒應(yīng)變能密度隨風(fēng)速增加呈指數(shù)增長(zhǎng)趨勢(shì),喉部相對(duì)水平位移與喉部直徑百分比(S)和塔筒應(yīng)變能密度和值(Pd)在風(fēng)速增長(zhǎng)過程中均出現(xiàn)明顯的轉(zhuǎn)折點(diǎn);基于變形和能量失效準(zhǔn)則,初步提出了喉部相對(duì)水平位移與喉部直徑百分比S>1%、指數(shù)應(yīng)變能密度和值破壞系數(shù)K≥0.3作為超大型冷卻塔在強(qiáng)風(fēng)作用下的結(jié)構(gòu)失效準(zhǔn)則。

    參考文獻(xiàn):

    [1] Gharaibeh E S, Dan M F, Onoufriou T. Reliability-based importance assessment of structural members with applications to complex structures[J]. Computers & Structures, 2002, 80(12): 1113-1131.

    [2] Report of the COI into collapse of cooling towers at Ferrybridge, Monday 1 November 1965[R]. Central Electricity Generating Board, 1966.

    [3] Pope R A. Structural deficiencies of natural draught cooling towers at UK power stations. Part 1: Failures at Ferrybridge and Fiddlers Ferry[J]. Proceedings of the ICE—Structures and Buildings, 1994, 104(1): 1-10.

    [4] 張軍鋒, 葛耀君, 趙林. 基于風(fēng)洞試驗(yàn)的雙曲線冷卻塔靜風(fēng)整體穩(wěn)定研究[J]. 工程力學(xué), 2012, 29(5): 68-77.

    Zhang Junfeng, Ge Yaojun, Zhao Lin. Study on global aerostatic stability of hyperboloidal cooling tower based on the wind tunnel tests[J]. Engineering Mechanics, 2012, 29(5): 68-77.

    [5] Ke Shitang, Wang Hao, Wang Tongguang, et al. Comparison of comprehensive stress performances of super-large cooling tower in different four-tower arrangements under 3D asymmetric wind loads[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2018, 179: 158-172.

    [6] 張軍鋒, 葛耀君, 趙林. 群塔布置對(duì)冷卻塔整體風(fēng)荷載和風(fēng)致響應(yīng)的不同干擾效應(yīng)[J]. 工程力學(xué), 2016, 33(8): 15-23.

    Zhang Junfeng, Ge Yaojun, Zhao Lin. Interference effects on global wind loads and wind induced responses for group hyperboloidal cooling towers[J]. Mechanics in Engineering, 2016, 33(8): 15-23.

    [7] Sabouri-Ghoumi S, Kharrazi M H K, Javidan P. Effect of stiffening rings on buckling stability of R.C. hyperbolic cooling towers[J]. Thin-Walled Structures, 2006, 44(2): 152-158.

    [8] 趙林, 李鵬飛, 葛耀君. 等效靜風(fēng)荷載下超大型冷卻塔受力性能分析[J]. 工程力學(xué), 2008, 25(7):79-86.

    Zhao Lin, Li Pengfei, Ge Yaojun. Numerical investigation on performance of super large cooling towers under equivalent static wind load[J]. Engineering Mechanics, 2008, 25(7): 79-86.

    [9] Ke Shitang, Xu Lu, Ge Yaojun. Study of random characteristics of fluctuating wind loads on ultra-large cooling towers in full construction process[J]. Wind and Structures, 2018, 26(4): 191-204.

    [10] 肖 ?南, 苗永志, 趙文爭(zhēng). 子午向幾何缺陷冷卻塔的結(jié)構(gòu)性能分析[J]. 浙江大學(xué)學(xué)報(bào)(工學(xué)版), 2010, 44(4): 819-825.

    Xiao Nan, Miao Yongzhi, Zhao Wenzheng. Structural performance analysis of a cooling-tower shell with meridional geometric imperfections[J]. Journal of Zhejiang University (Engineering Science), 2010, 44(4): 819-825.

    [11] 鄒云峰, 牛華偉, 陳政清. 基于完全氣動(dòng)彈性模型的冷卻塔風(fēng)致響應(yīng)風(fēng)洞試驗(yàn)研究[J]. 建筑結(jié)構(gòu)學(xué)報(bào), 2013, 34(6): 60-67.

    Zou Yunfeng, Niu Huawei, Chen Zhengqing. Wind tunnel test on wind-induced response of cooling tower based on full aero-elastic model [J]. Journal of Building Structures, 2013, 34(6): 60-67.

    [12] 柯世堂, 杜凌云. 不同氣動(dòng)措施對(duì)特大型冷卻塔風(fēng)致響應(yīng)及穩(wěn)定性能影響分析[J]. 湖南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2016,43(5): 79-89.

    Ke Shitang, Du Lingyun. Impact study on wind-induced response and stability for super large cooling tower with different aerodynamic measures[J]. Journal of Hunan University (Natural Sciences), 2016, 43(5): 79-89.

    [13] 陳寶春, 范冰輝, 余印根, 等. 鋼管混凝土拱橋強(qiáng)健性設(shè)計(jì)[J]. 橋梁建設(shè), 2016, 46(6): 88-93.

    Chen Baochun, Fan Binghui, Yu Yingen, et al. Robustness design of concrete-filled steel tube arch bridges[J]. Bridge Construction, 2016, 46(6): 88-93.

    [14] 鄭小博, 趙 ?煜, 賀拴海, 等. 雙塔鋼桁斜拉橋結(jié)構(gòu)強(qiáng)健性計(jì)算方法[J]. 交通工程運(yùn)輸學(xué)報(bào), 2017, 14(5): 27-38.

    Zheng Xiaobo, Zhao Yu, He Shuanhai, et al. Calculating method of structural robustness of double-tower cable-stayed bridge with steel truss girder[J]. Journal of Traffic and Transportation Engineering, 2017, 14(5): 27-38.

    [15] 鄧海斌, 陳 ?文, 王曉明. 基于AP法的無背索斜拉橋強(qiáng)健性分析[J]. 石家莊鐵道大學(xué)學(xué)報(bào), 2016, 29(3): 13-20.

    Deng Haibin, Chen Wen, Wang Xiaoming. Research on structural robustness for cable-stayed bridge without back-stays [J]. Journal of Shijiazhuang Tiedao University, 2016, 29 (3): 13-20.

    [16] 沈銳利, 房 ?凱, 官 ?快. 單根吊索斷裂時(shí)自錨式懸索橋強(qiáng)健性分析[J]. 橋梁建設(shè), 2014, 44 (6): 35-39.

    Shen Ruili, Fang Kai, Guan Kuai. Robustness analysis of self-anchored suspension bridge with loss of a single sling[J]. Bridge Construction, 2014, 44 (6): 35-39.

    [17] Ghosn M, Due?as-Osorio L, Frangopol D M, et al. Performance indicators for structural systems and infrastructure networks[J]. Journal of Structural Engineering, 2016, 142(9): F4016003.

    [18] 范冰輝, 陳寶春, 吳慶雄. 考慮強(qiáng)健性的中、下承式拱橋技術(shù)狀況評(píng)定[J]. 橋梁建設(shè), 2018, 48(5): 64-68.

    Fan Binghui, Chen Baochun, Wu Qingxiong. Technical condition evaluation of half-through and through arch bridges considering robustness [J]. Journal of Structural Engineering, 2018, 48(5): 64-68.

    [19] 柯世堂, 徐 ?璐, 朱 ?鵬. 基于大渦模擬超大型冷卻塔施工期風(fēng)荷載時(shí)域特性分析[J]. 湖南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2018, 45(11): 62-72.

    Ke Shitiang, Xu Lu, Zhu Peng. Analysis on time domain characteristics of wind loads during construction period of super large cooling towers based on large eddy simulation[J]. Journal of Hunan University (Natural Sciences), 2018, 45(11): 62-72.

    [20] 徐 ?璐, 柯世堂. 加勁環(huán)和子午肋對(duì)特大型冷卻塔靜風(fēng)響應(yīng)的影響研究[J]. 特種結(jié)構(gòu), 2016, 33(6): 23-31.

    Xu Lu, Ke Shitang. The influence study of stiffening ring and meridian rib for the static wind-induced responses on super large cooling towers[J]. Special Structures, 2016, 33(6): 23-31

    [21] 柯世堂, 王 ?浩. 考慮不同四塔形式特大型冷卻塔群風(fēng)荷載極值分布特征與取值探討[J]. 空氣動(dòng)力學(xué)學(xué)報(bào), 2020, 38(1): 48-57.

    Ke Shitang, Wang Hao. Distribution characteristics and value of extreme wind loads of four super-large cooling towers considering different layouts[J]. Acta Aerodynamica Sinica, 2020, 38(1): 48-57.

    [22] 電力規(guī)劃設(shè)計(jì)總院.DL/T 5339-2018, 火力發(fā)電廠水工設(shè)計(jì)規(guī)范[S]. 北京:中國(guó)計(jì)劃出版社, 2006.

    China Eleetric Power Planning & Engineering Institute. DL/T 5339-2018, Code for hydraulic design of fossil fuel power plants[S]. Beijing: China Planning Power Press, 2006.

    [23] 程小衛(wèi), 李 ?易,陸新征, 等. 基于多尺度模型的RC框架撞擊倒塌響應(yīng)數(shù)值分析[J]. 振動(dòng)與沖擊, 2016, 35(5): 82-88.

    Cheng Xiaowei, Li Yi, Lu Xinzheng, et al. Numerical analysis for collapse reponse of a RC frame subjected to impact loading based on multi-scale model [J]. Journal of Vibration and Shock, 2016, 35(5): 82-88.

    [24] 劉紅軍, 李洪江. 基于能量法的錨桿失效模糊判別方法研究[J]. 巖土工程學(xué)報(bào), 2013, 35(8): 1435-1441.

    Liu Hongjun, Li Hongjiang, Fuzzy identification method for failure of anchors based on energy principle[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(8): 1435-1441.

    [25] DoD2010. Design of structures to resist progressive collapse[S]. Washington, D.C: Department of Defense.

    [26] 張 ?明, 田始軒, 劉占輝, 等. 基于指數(shù)應(yīng)變能密度的弦支穹頂結(jié)構(gòu)失效準(zhǔn)則研究[J]. 振動(dòng)與沖擊, 2018, 37(11): 145-152.

    九九久久精品国产亚洲av麻豆| 一本精品99久久精品77| 国产精品一区二区三区四区久久| 欧美又色又爽又黄视频| 伊人久久精品亚洲午夜| 国产亚洲精品av在线| 天堂av国产一区二区熟女人妻| 国产在线精品亚洲第一网站| 免费观看精品视频网站| 色5月婷婷丁香| 搡老熟女国产l中国老女人| 99国产精品一区二区三区| 看免费av毛片| 亚洲精品乱码久久久v下载方式| 天堂影院成人在线观看| 免费在线观看日本一区| 亚洲国产日韩欧美精品在线观看| 国产白丝娇喘喷水9色精品| 女人十人毛片免费观看3o分钟| 久久久成人免费电影| 国产精品自产拍在线观看55亚洲| 国产精品久久视频播放| 日韩高清综合在线| aaaaa片日本免费| 久久精品国产99精品国产亚洲性色| 午夜福利成人在线免费观看| 亚洲精品日韩av片在线观看| 一级黄色大片毛片| 18+在线观看网站| 在线播放国产精品三级| 老司机午夜福利在线观看视频| 亚洲精品影视一区二区三区av| 国产高清视频在线观看网站| 3wmmmm亚洲av在线观看| 国产爱豆传媒在线观看| 我的老师免费观看完整版| 亚洲av中文字字幕乱码综合| 青草久久国产| 波野结衣二区三区在线| 看免费av毛片| 尤物成人国产欧美一区二区三区| 中出人妻视频一区二区| 亚洲人与动物交配视频| 99国产精品一区二区蜜桃av| 精品国产亚洲在线| 麻豆成人午夜福利视频| 毛片女人毛片| 亚洲五月婷婷丁香| 久久久久久国产a免费观看| 日本 av在线| 午夜福利高清视频| av专区在线播放| 欧美xxxx性猛交bbbb| 精品乱码久久久久久99久播| 丰满的人妻完整版| 欧美又色又爽又黄视频| 亚洲精品日韩av片在线观看| 在线观看免费视频日本深夜| 中文字幕高清在线视频| 欧美潮喷喷水| 久久精品国产99精品国产亚洲性色| www.999成人在线观看| 亚洲18禁久久av| 色噜噜av男人的天堂激情| 成人一区二区视频在线观看| 成熟少妇高潮喷水视频| 男人狂女人下面高潮的视频| 老司机午夜十八禁免费视频| 亚洲熟妇中文字幕五十中出| www日本黄色视频网| 国产野战对白在线观看| 嫩草影视91久久| 免费无遮挡裸体视频| 欧美日本视频| 窝窝影院91人妻| 一本精品99久久精品77| 亚洲人成网站在线播| 成人av在线播放网站| 亚洲av电影在线进入| 男女床上黄色一级片免费看| 99久久久亚洲精品蜜臀av| 在线天堂最新版资源| 18+在线观看网站| 怎么达到女性高潮| 成人一区二区视频在线观看| 又黄又爽又刺激的免费视频.| 中文字幕人妻熟人妻熟丝袜美| 国产精品美女特级片免费视频播放器| 国产欧美日韩一区二区精品| 桃红色精品国产亚洲av| 亚洲精华国产精华精| 久久久久久久精品吃奶| 欧美午夜高清在线| 久久热精品热| 精品不卡国产一区二区三区| 成人永久免费在线观看视频| 麻豆一二三区av精品| 国产又黄又爽又无遮挡在线| 亚洲久久久久久中文字幕| 精品人妻熟女av久视频| 两个人视频免费观看高清| 一二三四社区在线视频社区8| 欧美日韩黄片免| 欧美日韩综合久久久久久 | 一边摸一边抽搐一进一小说| 三级毛片av免费| 日韩中字成人| 国产黄a三级三级三级人| 日韩高清综合在线| 精品久久久久久久久久免费视频| 日韩人妻高清精品专区| 日本免费a在线| 国产精品伦人一区二区| 国产亚洲精品久久久久久毛片| 国产精品免费一区二区三区在线| 能在线免费观看的黄片| 午夜免费男女啪啪视频观看 | 亚洲色图av天堂| 草草在线视频免费看| 一边摸一边抽搐一进一小说| 国产在线男女| 美女大奶头视频| 99久久无色码亚洲精品果冻| 窝窝影院91人妻| 婷婷精品国产亚洲av在线| 我要搜黄色片| 色尼玛亚洲综合影院| 动漫黄色视频在线观看| 俄罗斯特黄特色一大片| 日本在线视频免费播放| 精品人妻偷拍中文字幕| 91在线精品国自产拍蜜月| 桃色一区二区三区在线观看| 午夜日韩欧美国产| 国产乱人视频| 国产三级黄色录像| 中亚洲国语对白在线视频| 精品99又大又爽又粗少妇毛片 | 国产精品野战在线观看| 黄色一级大片看看| 国产高清视频在线播放一区| 日本精品一区二区三区蜜桃| 日本一二三区视频观看| 极品教师在线视频| 亚洲av二区三区四区| 色吧在线观看| 中文字幕精品亚洲无线码一区| 男女那种视频在线观看| 黄色日韩在线| 九色成人免费人妻av| 亚洲成a人片在线一区二区| 欧美精品啪啪一区二区三区| 亚洲自偷自拍三级| 中出人妻视频一区二区| 99久久久亚洲精品蜜臀av| 18禁在线播放成人免费| 亚洲第一电影网av| 午夜福利在线在线| 性插视频无遮挡在线免费观看| 黄色丝袜av网址大全| 久久精品91蜜桃| 国产精品日韩av在线免费观看| 国产亚洲欧美98| 国产精品一区二区三区四区免费观看 | 亚洲欧美清纯卡通| 亚洲,欧美,日韩| 国产精品久久久久久精品电影| 久久久色成人| 日韩国内少妇激情av| 国产免费av片在线观看野外av| 最近在线观看免费完整版| 国产精品不卡视频一区二区 | 啦啦啦观看免费观看视频高清| 波多野结衣高清无吗| 欧美激情在线99| 高潮久久久久久久久久久不卡| 国产三级在线视频| 高清毛片免费观看视频网站| 欧美绝顶高潮抽搐喷水| 有码 亚洲区| 久久久国产成人精品二区| 波野结衣二区三区在线| 深爱激情五月婷婷| 成人美女网站在线观看视频| 亚洲精品成人久久久久久| 亚洲五月天丁香| 啦啦啦韩国在线观看视频| 美女大奶头视频| 欧美成人一区二区免费高清观看| 男女那种视频在线观看| 1000部很黄的大片| 亚洲,欧美精品.| 免费看美女性在线毛片视频| 嫩草影院入口| 精品99又大又爽又粗少妇毛片 | 赤兔流量卡办理| 夜夜看夜夜爽夜夜摸| 国产私拍福利视频在线观看| 亚洲精华国产精华精| 欧美极品一区二区三区四区| 国内精品一区二区在线观看| 欧美bdsm另类| 亚洲欧美清纯卡通| 亚洲在线自拍视频| 精品久久久久久久久亚洲 | 淫妇啪啪啪对白视频| 国产三级中文精品| 国产亚洲av嫩草精品影院| 国产精品人妻久久久久久| 中文字幕人成人乱码亚洲影| 51国产日韩欧美| 国产精品国产高清国产av| 亚洲国产欧洲综合997久久,| 亚洲久久久久久中文字幕| 午夜福利高清视频| 最后的刺客免费高清国语| 免费av观看视频| 草草在线视频免费看| 欧美xxxx性猛交bbbb| 婷婷亚洲欧美| 久久精品综合一区二区三区| 真实男女啪啪啪动态图| 国产午夜精品久久久久久一区二区三区 | 欧美日韩黄片免| 午夜亚洲福利在线播放| 赤兔流量卡办理| 日韩欧美一区二区三区在线观看| 欧美丝袜亚洲另类 | 久久婷婷人人爽人人干人人爱| 日韩 亚洲 欧美在线| 老熟妇乱子伦视频在线观看| 中文字幕久久专区| 亚洲精品成人久久久久久| 99精品久久久久人妻精品| 激情在线观看视频在线高清| 我要搜黄色片| 亚洲中文日韩欧美视频| 久久久久久久久大av| 日本撒尿小便嘘嘘汇集6| 午夜两性在线视频| 精品免费久久久久久久清纯| 欧美丝袜亚洲另类 | 波多野结衣高清作品| 一级av片app| 精品一区二区三区视频在线观看免费| 欧美激情在线99| 波多野结衣巨乳人妻| 日本 av在线| 久久久国产成人免费| 国产毛片a区久久久久| 亚洲精品粉嫩美女一区| 亚洲美女搞黄在线观看 | 亚洲精品一卡2卡三卡4卡5卡| 久久人人精品亚洲av| 亚洲成a人片在线一区二区| 老司机午夜福利在线观看视频| 麻豆一二三区av精品| 少妇高潮的动态图| 欧美激情国产日韩精品一区| 亚洲乱码一区二区免费版| 淫秽高清视频在线观看| 日韩 亚洲 欧美在线| 在线国产一区二区在线| 国产午夜精品久久久久久一区二区三区 | 丝袜美腿在线中文| 日本三级黄在线观看| 男人舔奶头视频| 性色avwww在线观看| 欧美+日韩+精品| 天堂√8在线中文| 国产精品久久电影中文字幕| 色在线成人网| 日韩亚洲欧美综合| 国产亚洲精品久久久久久毛片| 嫩草影院新地址| 日日摸夜夜添夜夜添小说| 国产91精品成人一区二区三区| 蜜桃亚洲精品一区二区三区| 神马国产精品三级电影在线观看| 最后的刺客免费高清国语| 韩国av一区二区三区四区| 黄色女人牲交| 一本一本综合久久| 精品熟女少妇八av免费久了| 国内毛片毛片毛片毛片毛片| 亚洲av电影不卡..在线观看| 亚洲欧美日韩东京热| 黄片小视频在线播放| 国产真实乱freesex| 一个人看视频在线观看www免费| 日韩欧美一区二区三区在线观看| 热99re8久久精品国产| 国产精品一区二区三区四区免费观看 | 日韩中字成人| 嫩草影院入口| 又黄又爽又刺激的免费视频.| 12—13女人毛片做爰片一| 亚洲av日韩精品久久久久久密| 国产精品伦人一区二区| 网址你懂的国产日韩在线| 精品一区二区免费观看| 精品人妻熟女av久视频| 国产精品一区二区三区四区免费观看 | 我要搜黄色片| 亚洲精品乱码久久久v下载方式| 久久久久久九九精品二区国产| 变态另类丝袜制服| 精品熟女少妇八av免费久了| 真人做人爱边吃奶动态| 日韩人妻高清精品专区| 丰满乱子伦码专区| 狂野欧美白嫩少妇大欣赏| 狂野欧美白嫩少妇大欣赏| 别揉我奶头 嗯啊视频| 看十八女毛片水多多多| 国产成年人精品一区二区| 午夜精品久久久久久毛片777| 久久人妻av系列| 午夜两性在线视频| 哪里可以看免费的av片| 午夜福利成人在线免费观看| 久久精品人妻少妇| 中文字幕人成人乱码亚洲影| 亚洲欧美清纯卡通| 国产精品不卡视频一区二区 | 51国产日韩欧美| 日韩欧美在线乱码| 亚洲自偷自拍三级| 亚洲avbb在线观看| 不卡一级毛片| 亚洲熟妇中文字幕五十中出| 久99久视频精品免费| 久久精品国产亚洲av涩爱 | 欧美日韩福利视频一区二区| 在线国产一区二区在线| 此物有八面人人有两片| 综合色av麻豆| 男人和女人高潮做爰伦理| 亚洲在线观看片| 精华霜和精华液先用哪个| 69人妻影院| 亚洲天堂国产精品一区在线| 欧美一区二区精品小视频在线| 亚洲av熟女| 看片在线看免费视频| 日韩欧美三级三区| 又爽又黄无遮挡网站| 麻豆久久精品国产亚洲av| 免费大片18禁| 国产又黄又爽又无遮挡在线| 97碰自拍视频| 欧美成人a在线观看| 成人一区二区视频在线观看| 日韩欧美精品v在线| 51国产日韩欧美| 国产精品久久久久久久电影| 日本免费a在线| 人妻丰满熟妇av一区二区三区| 此物有八面人人有两片| 精品久久久久久久久av| 久久精品影院6| 三级毛片av免费| 悠悠久久av| 国产黄色小视频在线观看| 别揉我奶头 嗯啊视频| 亚洲 国产 在线| 69av精品久久久久久| 久久国产精品人妻蜜桃| 久久精品91蜜桃| 一区二区三区四区激情视频 | 成人国产综合亚洲| 久久精品影院6| 欧美极品一区二区三区四区| 悠悠久久av| 一区二区三区免费毛片| 赤兔流量卡办理| 国产欧美日韩精品亚洲av| 男女那种视频在线观看| 十八禁网站免费在线| 国产成人a区在线观看| 我的女老师完整版在线观看| 亚洲人成网站高清观看| 色哟哟·www| 国产精品电影一区二区三区| 久久久精品欧美日韩精品| 日韩欧美在线乱码| 午夜亚洲福利在线播放| 欧美日韩瑟瑟在线播放| 99国产精品一区二区蜜桃av| 国产老妇女一区| 最近最新免费中文字幕在线| 白带黄色成豆腐渣| a在线观看视频网站| 九九在线视频观看精品| a在线观看视频网站| 精品国产亚洲在线| 国产成人aa在线观看| 国产乱人视频| 性欧美人与动物交配| www日本黄色视频网| 99久久精品国产亚洲精品| 无遮挡黄片免费观看| 精华霜和精华液先用哪个| 男人舔女人下体高潮全视频| 亚洲精品一卡2卡三卡4卡5卡| 午夜福利在线在线| 欧美zozozo另类| 老女人水多毛片| 乱码一卡2卡4卡精品| 欧美又色又爽又黄视频| 久久6这里有精品| 国产精品一区二区免费欧美| 国产精品一区二区性色av| 欧美激情在线99| 日本在线视频免费播放| 精品人妻偷拍中文字幕| 深爱激情五月婷婷| 亚洲国产精品久久男人天堂| 不卡一级毛片| 国产国拍精品亚洲av在线观看| 极品教师在线视频| 欧美日韩综合久久久久久 | 午夜久久久久精精品| 国产极品精品免费视频能看的| 国产av不卡久久| 久久久国产成人精品二区| 老司机深夜福利视频在线观看| 91麻豆精品激情在线观看国产| 国产精品伦人一区二区| 成人毛片a级毛片在线播放| 岛国在线免费视频观看| 精品一区二区三区av网在线观看| 成人三级黄色视频| 1024手机看黄色片| 变态另类成人亚洲欧美熟女| 天美传媒精品一区二区| 亚洲中文日韩欧美视频| 亚洲欧美日韩无卡精品| 欧美xxxx黑人xx丫x性爽| 在线播放无遮挡| 国产一区二区在线观看日韩| 免费在线观看成人毛片| 日本黄大片高清| 99久久无色码亚洲精品果冻| 搡老熟女国产l中国老女人| 国产91精品成人一区二区三区| 亚洲人成伊人成综合网2020| 99热精品在线国产| 级片在线观看| .国产精品久久| 国产精品免费一区二区三区在线| 午夜福利视频1000在线观看| 99热精品在线国产| 99国产极品粉嫩在线观看| 免费大片18禁| 成人亚洲精品av一区二区| 亚洲经典国产精华液单 | 99热只有精品国产| 亚洲av免费高清在线观看| xxxwww97欧美| 亚洲成人免费电影在线观看| 国产淫片久久久久久久久 | 亚洲av一区综合| 日日夜夜操网爽| 国产精品一区二区三区四区久久| 国产探花在线观看一区二区| 国产一区二区三区在线臀色熟女| 婷婷亚洲欧美| 男人的好看免费观看在线视频| 91久久精品电影网| 亚洲成a人片在线一区二区| 色噜噜av男人的天堂激情| 亚洲精品乱码久久久v下载方式| 精品乱码久久久久久99久播| 国产精品野战在线观看| 噜噜噜噜噜久久久久久91| 真人做人爱边吃奶动态| 亚洲内射少妇av| 好看av亚洲va欧美ⅴa在| 欧美成人a在线观看| 国产精品亚洲一级av第二区| 色播亚洲综合网| 国产免费男女视频| 免费av毛片视频| 亚洲内射少妇av| 精品午夜福利视频在线观看一区| 老熟妇仑乱视频hdxx| 午夜免费激情av| h日本视频在线播放| 69人妻影院| 国产亚洲精品综合一区在线观看| 精品午夜福利在线看| 日本在线视频免费播放| 亚洲美女视频黄频| 特大巨黑吊av在线直播| 国产三级在线视频| 国产成人av教育| 一进一出抽搐gif免费好疼| 亚洲精品久久国产高清桃花| 757午夜福利合集在线观看| 偷拍熟女少妇极品色| 亚洲欧美日韩高清在线视频| av在线观看视频网站免费| 午夜免费激情av| 在线十欧美十亚洲十日本专区| 丰满乱子伦码专区| 十八禁国产超污无遮挡网站| 久久99热这里只有精品18| 夜夜夜夜夜久久久久| 在线观看av片永久免费下载| 久9热在线精品视频| 成年女人看的毛片在线观看| 岛国在线免费视频观看| 午夜福利视频1000在线观看| 9191精品国产免费久久| 国产三级黄色录像| 国产色婷婷99| 在线免费观看的www视频| 久久久国产成人精品二区| 久久精品人妻少妇| 国产欧美日韩一区二区精品| 亚洲电影在线观看av| 听说在线观看完整版免费高清| 日本三级黄在线观看| 午夜福利18| 成年女人毛片免费观看观看9| 村上凉子中文字幕在线| 在线看三级毛片| 国产视频一区二区在线看| 一二三四社区在线视频社区8| 国产综合懂色| 亚洲av免费在线观看| 久久精品国产99精品国产亚洲性色| 国产精品久久久久久久电影| 性插视频无遮挡在线免费观看| 免费av不卡在线播放| 毛片一级片免费看久久久久 | 久久99热6这里只有精品| 神马国产精品三级电影在线观看| 俺也久久电影网| 国产精品综合久久久久久久免费| 极品教师在线视频| 哪里可以看免费的av片| 亚洲欧美日韩东京热| 99在线视频只有这里精品首页| 免费看a级黄色片| 最近最新中文字幕大全电影3| 欧美性猛交黑人性爽| 激情在线观看视频在线高清| 99精品久久久久人妻精品| 18禁黄网站禁片免费观看直播| 国产色爽女视频免费观看| 欧美性猛交╳xxx乱大交人| 看十八女毛片水多多多| 精品熟女少妇八av免费久了| www日本黄色视频网| 国产色婷婷99| 日韩av在线大香蕉| 欧美区成人在线视频| 中文字幕精品亚洲无线码一区| 亚洲综合色惰| 国产亚洲精品av在线| 18美女黄网站色大片免费观看| 嫩草影视91久久| 一本精品99久久精品77| 听说在线观看完整版免费高清| 夜夜看夜夜爽夜夜摸| 一级av片app| 亚洲18禁久久av| 老司机深夜福利视频在线观看| 最近在线观看免费完整版| 精华霜和精华液先用哪个| 美女被艹到高潮喷水动态| 嫩草影院入口| 国产爱豆传媒在线观看| 欧美高清性xxxxhd video| 琪琪午夜伦伦电影理论片6080| 国语自产精品视频在线第100页| 亚洲成人久久爱视频| av欧美777| 真人做人爱边吃奶动态| 最后的刺客免费高清国语| 亚洲欧美日韩无卡精品| 成人一区二区视频在线观看| 男女视频在线观看网站免费| 嫩草影视91久久| 最近最新中文字幕大全电影3| 搡老妇女老女人老熟妇| 亚洲成人久久性| 99热这里只有精品一区| 欧美日韩中文字幕国产精品一区二区三区| 黄色日韩在线| 亚洲欧美激情综合另类| 在线免费观看的www视频| 757午夜福利合集在线观看| 狂野欧美白嫩少妇大欣赏| 色综合亚洲欧美另类图片| 亚洲国产日韩欧美精品在线观看| 日本一二三区视频观看| 少妇人妻一区二区三区视频| 免费在线观看日本一区| 免费搜索国产男女视频| 波多野结衣高清无吗| 国产又黄又爽又无遮挡在线| 欧美成人性av电影在线观看| 亚洲专区中文字幕在线| 国产 一区 欧美 日韩| 日日干狠狠操夜夜爽| 国内少妇人妻偷人精品xxx网站| 国产在视频线在精品| 制服丝袜大香蕉在线| 两个人的视频大全免费| 高清在线国产一区| 一夜夜www| 国内少妇人妻偷人精品xxx网站|