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

    基于Info—Gap決策的結(jié)構(gòu)抗震穩(wěn)健性優(yōu)化設(shè)計(jì)

    2015-06-16 15:00唐和生范德偉李大偉薛松濤
    關(guān)鍵詞:穩(wěn)健性抗震設(shè)計(jì)

    唐和生+范德偉+李大偉+薛松濤

    摘要:針對(duì)結(jié)構(gòu)抗震設(shè)計(jì)中存在嚴(yán)重不確定性問(wèn)題,基于Info-Gap理論建立一種考慮地震設(shè)計(jì)譜參數(shù)不確定的結(jié)構(gòu)抗震穩(wěn)健性優(yōu)化設(shè)計(jì)方法.該設(shè)計(jì)方法采用Info-Gap模型來(lái)描述地震設(shè)計(jì)譜中反映地面運(yùn)動(dòng)強(qiáng)度的水平地震影響系數(shù)最大值αmax和場(chǎng)地特征周期Tg的不確定,通過(guò)嵌套優(yōu)化使結(jié)構(gòu)設(shè)計(jì)滿足結(jié)構(gòu)的臨界性能要求的同時(shí)實(shí)現(xiàn)最大化不確定的穩(wěn)健性.通過(guò)對(duì)一個(gè)6層3跨的鋼框架的抗震穩(wěn)健性優(yōu)化設(shè)計(jì)驗(yàn)證分析表明:這種基于滿足性能的結(jié)構(gòu)抗震設(shè)計(jì)需在滿足結(jié)構(gòu)性能需求和提高不確定穩(wěn)健性之間進(jìn)行權(quán)衡取舍;同時(shí)也證實(shí)該方法為在不易得到不確定性因素足夠信息情況下的基于性能結(jié)構(gòu)抗震可靠性設(shè)計(jì)提供了一條新思路.

    關(guān)鍵詞:Info-Gap理論;不確定;穩(wěn)健性;抗震設(shè)計(jì)

    中圖分類號(hào):TU318;TU323.5 文獻(xiàn)標(biāo)識(shí)碼:A

    Info-Gap Decision for the Robust Seismic

    Design Optimization of Structures

    TANG He-sheng1, 2, FAN De-wei2, LI Da-wei2, XUE Song-tao2,3

    (1.State Key Laboratory of Disaster Prevention in Civil Engineering, Tongji Univ, Shanghai200092, China;

    2. Research Institute of Structural Engineering and Disaster Reduction, Tongji Univ, Shanghai200092, China;

    3. Dept of Architecture, Tohoku Institute of Technology, Sendai982-8577, Japan)

    Abstract:Seismic design for buildings is usually subject to various uncertainties, often severe, which have the potential to undermine engineering decisions. It is crucial that these uncertainties be accounted for in seismic design. We formulated a performance-based seismic design model that takes into account uncertainty in the seismic design spectrum of the αmax and Tg. We used ?info-gap theory for satisfying the critical performance requirements, while at the same time maximized the robustness to uncertainty through nested optimization. The design implications of this robust-satisfying approach were demonstrated with a three-span six-floor steel frame design example. It is shown that design preferences depend upon the performance requirements considering the trade-off between robustness to uncertainty. Also, the result reveals that the proposed method provides a novel tool for the performance-based seismic reliability design under the lack of knowledge.

    Key words:Info-Gap theory; uncertainty; robust; seismic design

    結(jié)構(gòu)抗震設(shè)計(jì)中涉及到很多不確定的因素,如中國(guó)現(xiàn)行規(guī)范中的設(shè)計(jì)反應(yīng)譜便是通過(guò)對(duì)大量實(shí)際地震記錄的反應(yīng)譜進(jìn)行統(tǒng)計(jì)分析,并結(jié)合工程經(jīng)驗(yàn)判斷后得到的平均譜,以地震影響系數(shù)的形式給出.建筑結(jié)構(gòu)的地震影響系數(shù)應(yīng)根據(jù)烈度、場(chǎng)地類別、設(shè)計(jì)地震動(dòng)分組和結(jié)構(gòu)的自振周期以及阻尼比等確定.而對(duì)于給定的結(jié)構(gòu),其自振周期和阻尼比通常是確定的,進(jìn)而可直接由抗震設(shè)防烈度和場(chǎng)地特征周期來(lái)確定地震影響系數(shù)的取值[1].由此可知,規(guī)范中的反應(yīng)譜主要取決于反映地面運(yùn)動(dòng)強(qiáng)度的地震影響系數(shù)最大值αmax和場(chǎng)地的特征周期Tg.而地震影響系數(shù)最大值αmax和場(chǎng)地特征周期Tg會(huì)因震源機(jī)制、傳播途徑特性及局部地質(zhì)土質(zhì)條件等不同而存在著一定幅度的變化,并不是一個(gè)確定的值.再者由于規(guī)范反應(yīng)譜的統(tǒng)計(jì)值會(huì)受到統(tǒng)計(jì)分析所用資料的完備性、計(jì)算回歸過(guò)程中數(shù)據(jù)的離散性以及各種人為的主觀性的影響,使得αmax和Tg的取值都是不確定的.

    目前不確定性的結(jié)構(gòu)優(yōu)化設(shè)計(jì)主要采用的概率和模糊方法都依賴于大量的統(tǒng)計(jì)數(shù)據(jù),并需要對(duì)概率分布或隸屬度進(jìn)行主觀假定.研究表明概率模型參數(shù)的小偏差可導(dǎo)致結(jié)構(gòu)可靠性計(jì)算出現(xiàn)較大誤差[2-4].然而,很多結(jié)構(gòu)工程統(tǒng)計(jì)樣本的采集相對(duì)比較困難,難以得到不確定參量的足夠信息來(lái)確定其概率密度分布,只能得到其幅度或界限.因此,以概率方法為代表的不確定優(yōu)化方法并不完全適用于結(jié)構(gòu)工程領(lǐng)域[5].

    20世紀(jì)90年代,Ben-Haim[2]和Elishakoff等人[6]提出并倡導(dǎo)應(yīng)用凸集合模型描述區(qū)間不確定性.凸集合模型僅需要少量不確定性參數(shù)的統(tǒng)計(jì)信息,特別適用于嚴(yán)重不確定的優(yōu)化設(shè)計(jì).Ben-Haim在文獻(xiàn)[7]中基于凸集合理論,提出了非概率可靠性的思想:若系統(tǒng)能容許不確定參量在一定范圍內(nèi)的波動(dòng),則認(rèn)為系統(tǒng)是可靠的,并指出結(jié)構(gòu)可靠性可以用結(jié)構(gòu)在失效前能夠承受的不確定性總量來(lái)度量.在結(jié)構(gòu)設(shè)計(jì)方面,文獻(xiàn)[8-9]基于凸模型將非概率可靠性理論應(yīng)用于結(jié)構(gòu)設(shè)計(jì).蘇永華等[10]基于區(qū)間數(shù)學(xué)理論,建立了深部工程圍巖力學(xué)參數(shù)及其對(duì)開(kāi)挖響應(yīng)的區(qū)間非概率分析模型.但是當(dāng)不確定性參數(shù)采用區(qū)間數(shù)或凸集合描述時(shí),由于復(fù)雜結(jié)構(gòu)情形的可靠性問(wèn)題通常所呈現(xiàn)的非凸特性給此類方法帶來(lái)了分析計(jì)算上的困難,而且區(qū)間算法可能引起擴(kuò)張,從而使該方法的應(yīng)用受到較大限制.

    在凸集合理論研究的基礎(chǔ)上,Ben-Haim提出了“嚴(yán)重不確定性”的概念,即僅能獲知不確定性參數(shù)假定的名義值,但無(wú)法獲得不確定性參數(shù)其他信息的情況;并針對(duì)這類問(wèn)題創(chuàng)立了Info-Gap決策理論(Info-Gap Decision Theory, IGDT)[11].IGDT作為一種考慮嚴(yán)重不確定性的非概率優(yōu)化方法,已經(jīng)成功地應(yīng)用于各個(gè)領(lǐng)域,包括神經(jīng)網(wǎng)絡(luò)[12]、生物保護(hù)[13- 14]、金融經(jīng)濟(jì)[15]等.一些學(xué)者也嘗試將IGDT應(yīng)用于結(jié)構(gòu)工程領(lǐng)域中[5, 16],但這些研究都僅限于結(jié)構(gòu)的穩(wěn)健性函數(shù)解析求解簡(jiǎn)單問(wèn)題,對(duì)復(fù)雜的嵌套優(yōu)化問(wèn)題則鮮有涉及.鑒于此,本文提出了一種將IGDT與自然計(jì)算方法相結(jié)合的優(yōu)化方法,應(yīng)用人工智能算法來(lái)處理IGDT的嵌套優(yōu)化問(wèn)題,并將該方法應(yīng)用于基于性能的結(jié)構(gòu)抗震優(yōu)化中.最后通過(guò)一例鋼框架結(jié)構(gòu)的抗震性能優(yōu)化設(shè)計(jì),驗(yàn)證該方法處理基于性能的抗震穩(wěn)健優(yōu)化設(shè)計(jì)的有效性.

    1 Info-Gap不確定穩(wěn)健性優(yōu)化模型

    基于非概率Info-Gap不確定模型的結(jié)構(gòu)穩(wěn)健性設(shè)計(jì)是一種滿足性能的優(yōu)化設(shè)計(jì),其目的是設(shè)計(jì)出滿足所有約束條件及臨界性能的穩(wěn)健性最優(yōu)的結(jié)構(gòu).其基本思想為:首先選擇一個(gè)合適的臨界性能作為約束,然后在保證結(jié)構(gòu)系統(tǒng)正常工作的前提下,最大化結(jié)構(gòu)系統(tǒng)抵抗不確定變量波動(dòng)的能力,該模型主要由以下幾部分組成.

    1)設(shè)計(jì)變量:在結(jié)構(gòu)優(yōu)化設(shè)計(jì)中,首先是確定設(shè)計(jì)變量以達(dá)到理想設(shè)計(jì)的結(jié)構(gòu)性能或其他目標(biāo).可以表示為:

    q = [ q1, q2, …, qn ]T , qmin ≤q≤qmax . (1)

    式中:n為設(shè)計(jì)變量的個(gè)數(shù);qmin, qmax分別為設(shè)計(jì)變量的下限與上限.在結(jié)構(gòu)優(yōu)化設(shè)計(jì)中,設(shè)計(jì)變量包括幾何尺寸、形狀或結(jié)構(gòu)拓?fù)涞冉Y(jié)構(gòu)設(shè)計(jì)參數(shù).

    2)不確定變量:在結(jié)構(gòu)設(shè)計(jì)中,除設(shè)計(jì)變量外,還存在一些不確定變量.這些不確定變量u來(lái)源于許多方面,如外荷載、邊界條件、幾何參數(shù)以及材料參數(shù)等.圖1為一簡(jiǎn)單的Info-Gap模型,假定不確定變量的名義值為,而其真實(shí)值u未知,這兩者之間存在的差距以α來(lái)表示,且這種差距的程度是未知的.Info-Gap不確定模型將α稱為不確定性水平.

    3)約束條件:通常結(jié)構(gòu)設(shè)計(jì)問(wèn)題還必須考慮約束條件,以滿足設(shè)計(jì)變量空間的可行域要求.約束條件的表達(dá)式為:

    gi(q)≤0,i=1,2,…,I. ?(3)

    式中:I為約束條件的個(gè)數(shù).結(jié)構(gòu)設(shè)計(jì)中,約束條件通常需要考慮應(yīng)力約束和位移約束等.

    4)性能函數(shù):性能函數(shù)R(q, u) 是設(shè)計(jì)變量q和不確定變量u的函數(shù),代表結(jié)構(gòu)性能指標(biāo)或響應(yīng).性能函數(shù)可以是結(jié)構(gòu)造價(jià)、材料體積或結(jié)構(gòu)質(zhì)量,也可以是結(jié)構(gòu)柔度、節(jié)點(diǎn)位移、應(yīng)力、自振頻率或屈曲荷載等結(jié)構(gòu)性能.臨界性能rc為性能函數(shù)R(q, u)必須滿足的前提條件.這兩者構(gòu)成了非概率Info-Gap穩(wěn)健性設(shè)計(jì)過(guò)程時(shí)必須滿足的臨界性能條件(假設(shè)性能函數(shù)必須小于臨界性能rc):

    R(q,u)≤rc. ?(4)

    通常情況下,性能函數(shù)是設(shè)計(jì)變量與不確定變量的隱函數(shù),需要通過(guò)結(jié)構(gòu)分析(如有限元分析)來(lái)計(jì)算.

    5)穩(wěn)健性函數(shù):傳統(tǒng)的優(yōu)化設(shè)計(jì)是將性能函數(shù)R(q, u)作為目標(biāo)函數(shù),以實(shí)現(xiàn)性能最優(yōu)化的結(jié)構(gòu)設(shè)計(jì).而基于Info-Gap的結(jié)構(gòu)穩(wěn)健性優(yōu)化設(shè)計(jì)的目的是在滿足臨界性能條件下最大化不確定性水平αi.數(shù)學(xué)表達(dá)式如下:

    i(q,rc)=max {αi:(max R(q,u)≤rc)},

    i=1,2,…,I.(5)

    式中:I為設(shè)計(jì)中所需考慮的不確定性水平的個(gè)數(shù),當(dāng)存在多個(gè)不確定性水平時(shí),基于Info-Gap模型的結(jié)構(gòu)穩(wěn)健性設(shè)計(jì)轉(zhuǎn)變?yōu)槎嗄繕?biāo)設(shè)計(jì).穩(wěn)健性i是第i個(gè)不確定性水平αi的最大值,即結(jié)構(gòu)在滿足約束條件及臨界性能條件下所能承受的不確定變量實(shí)際值偏離其名義值的最大極限.

    2基于Info-Gap理論的結(jié)構(gòu)抗震性能穩(wěn)健

    性優(yōu)化設(shè)計(jì)

    傳統(tǒng)的性能優(yōu)化設(shè)計(jì)將結(jié)構(gòu)的性能指標(biāo)作為優(yōu)化目標(biāo),例如得到造價(jià)最小的結(jié)構(gòu)設(shè)計(jì).在不確定性結(jié)構(gòu)設(shè)計(jì)中,以概率方法為代表的性能優(yōu)化設(shè)計(jì)需要足夠充分的統(tǒng)計(jì)樣本,或?qū)Σ淮_定性參數(shù)進(jìn)行經(jīng)驗(yàn)假設(shè),難以建立不確定性參數(shù)信息嚴(yán)重缺失的結(jié)構(gòu)設(shè)計(jì)模型.為解決此類問(wèn)題,本文基于Info-Gap模型[11]以穩(wěn)健性函數(shù)來(lái)表示結(jié)構(gòu)系統(tǒng)的最大非概率不確定性水平,建立最大化穩(wěn)健性函數(shù)的結(jié)構(gòu)穩(wěn)健性設(shè)計(jì)策略及優(yōu)化方法,將穩(wěn)健性函數(shù)的最大化問(wèn)題轉(zhuǎn)化為考慮約束條件的嵌套優(yōu)化問(wèn)題,并采用人工智能算法處理穩(wěn)健性函數(shù)的嵌套優(yōu)化問(wèn)題.

    如上所述,基于Info-Gap理論的結(jié)構(gòu)抗震性能穩(wěn)健性優(yōu)化設(shè)計(jì)的過(guò)程就是通過(guò)公式(5)的嵌套優(yōu)化來(lái)求解穩(wěn)健性函數(shù).該優(yōu)化過(guò)程包括外優(yōu)化和內(nèi)優(yōu)化兩個(gè)部分,外優(yōu)化的目標(biāo)是求解穩(wěn)健性函數(shù),內(nèi)優(yōu)化的目標(biāo)是求解滿足臨界性能條件的最不利性能函數(shù).嵌套優(yōu)化的復(fù)雜程度遠(yuǎn)大于單一的優(yōu)化過(guò)程,這使得精確計(jì)算結(jié)構(gòu)穩(wěn)健性函數(shù)變得非常困難.因此,很多學(xué)者轉(zhuǎn)而研究可行的數(shù)值優(yōu)化問(wèn)題.在Info-Gap模型的發(fā)展過(guò)程中,也是將不確定性問(wèn)題假定為滿足凸集合條件下的不確定性問(wèn)題,并針對(duì)穩(wěn)健性函數(shù)和其他部分提出了相應(yīng)的計(jì)算方法.但是在結(jié)構(gòu)優(yōu)化中,結(jié)構(gòu)不確定性參數(shù)采用區(qū)間數(shù)或凸集合描述時(shí),結(jié)構(gòu)響應(yīng)值域通常是非凸的,且具有不止一個(gè)局部極值點(diǎn)[17].因此,實(shí)際結(jié)構(gòu)的不確定性優(yōu)化問(wèn)題相對(duì)于單一的優(yōu)化問(wèn)題來(lái)說(shuō)更需要采用全局優(yōu)化方法.

    遺傳算法(GA)[18]是模擬生物自然進(jìn)化機(jī)制的隨機(jī)優(yōu)化方法,它借用了仿真生物遺傳學(xué)和自然選擇機(jī)理,通過(guò)自然選擇、遺傳、變異等作用機(jī)制,實(shí)現(xiàn)個(gè)體的適應(yīng)性的提高,末代適應(yīng)性最高的個(gè)體作為近似最優(yōu)解.GA對(duì)設(shè)計(jì)變量特殊編碼,其尋優(yōu)搜索由適應(yīng)度值為引導(dǎo),對(duì)不可微以至不連續(xù)的函數(shù)優(yōu)化,GA能以較大概率求得全局解;具有較強(qiáng)的魯棒性、全局收斂性、隱含并行性及廣泛的適應(yīng)性;并且能處理不同類型的優(yōu)化變量(離散的、連續(xù)的和混合型的);直接搜索法不需要任何輔助信息,對(duì)目標(biāo)函數(shù)和約束函數(shù)也沒(méi)有任何要求.由于上述優(yōu)點(diǎn),GA在眾多學(xué)者的努力下已經(jīng)發(fā)展成為一種自組織、自適應(yīng)啟發(fā)式概率性搜索優(yōu)化技術(shù),能有效地求解屬于NPC(不可計(jì)算的)類型的組合優(yōu)化問(wèn)題及非線性多模型、多目標(biāo)的函數(shù)優(yōu)化問(wèn)題,從而得到了多學(xué)科的廣泛重視.GA開(kāi)始滲透到自適應(yīng)行為、人工智能、神經(jīng)網(wǎng)絡(luò)[19]、機(jī)器人、運(yùn)籌學(xué)等學(xué)科[20-21],并逐漸擴(kuò)展到工程技術(shù)、計(jì)算機(jī)科學(xué)、管理科學(xué)和社會(huì)科學(xué)等領(lǐng)域[22].基于上述優(yōu)點(diǎn),本文采用GA來(lái)實(shí)現(xiàn)基于Info-Gap的結(jié)構(gòu)抗震性能穩(wěn)健性優(yōu)化的計(jì)算.

    2.1 Info-Gap結(jié)構(gòu)抗震性能穩(wěn)健性設(shè)計(jì)步驟

    基于Info-Gap結(jié)構(gòu)抗震性能穩(wěn)健性設(shè)計(jì)的GA算法主要步驟如下.

    1) 開(kāi)始外優(yōu)化:設(shè)置GA算法參數(shù),確定臨界性能rc,在可行空間內(nèi)初始化設(shè)計(jì)變量q;

    2) 開(kāi)始內(nèi)優(yōu)化:初始化不確定變量u;

    3) 以設(shè)計(jì)變量q和不確定性變量u進(jìn)行結(jié)構(gòu)分析,求解性能函數(shù)R(q, u);

    4) 更新不確定性變量u;

    5) 返回3),直至達(dá)到內(nèi)優(yōu)化最大迭代代數(shù);

    6) 輸出當(dāng)前滿足臨界性能條件的最不利性能R(q, u)及相應(yīng)的不確定變量u;

    7) 更新設(shè)計(jì)變量q;

    8) 返回2),直至達(dá)到外優(yōu)化最大迭代代數(shù);

    9) 輸出當(dāng)前穩(wěn)健性函數(shù)及相應(yīng)最優(yōu)設(shè)計(jì)變量.

    結(jié)構(gòu)穩(wěn)健性設(shè)計(jì)過(guò)程中,嵌套內(nèi)優(yōu)化和外優(yōu)化分別需要計(jì)算性能函數(shù)R(q, u)與不確定性水平αi的極值,極值的求解通過(guò)GA來(lái)實(shí)現(xiàn).在本文中定義不確定性水平u為抗震設(shè)計(jì)反應(yīng)譜中的水平地震影響系數(shù)最大值αmax和場(chǎng)地特征周期Tg.定義內(nèi)優(yōu)化性能函數(shù)為結(jié)構(gòu)的最大層間位移角θi,max(q, u).

    2.2 抗震性能優(yōu)化設(shè)計(jì)中性能函數(shù)的求解

    本文采用基于文獻(xiàn)[23]中μ-ξy-T關(guān)系的延性需求譜法而開(kāi)發(fā)求解鋼框架地震性能指標(biāo)的程序,求解結(jié)構(gòu)性能點(diǎn)的位移和層間位移角,求解步驟如下.

    1) 求解結(jié)構(gòu)的動(dòng)力特性:周期T以及振型Φ,其中振型以頂點(diǎn)幅值為1進(jìn)行標(biāo)準(zhǔn)化.

    2) 對(duì)第一階模態(tài),建立結(jié)構(gòu)的pushover曲線.將pushover曲線等效為二折線型,并求得等效單自由度(ESDOF)的周期Te與屈服位移Δy.

    3) 根據(jù)彈性加速度與位移之間關(guān)系求解等效單自由度體系的需求位移Sde并與Δy相對(duì)比:①若Sde< Δy則說(shuō)明單自由度體系尚處于彈性狀態(tài),即Sde就是等效單自由度體系的最大響應(yīng)位移Δu;②若Sde>Δy則說(shuō)明單自由度體系在該水準(zhǔn)下已處于彈塑性狀態(tài),利用文獻(xiàn)[23]中μ-ξy-T關(guān)系求解延性系數(shù)μ,進(jìn)而求解單自由度體系最大響應(yīng)位移Δu.

    4) 利用Δu=φTMφφTMΔu求解結(jié)構(gòu)頂部最大位移,再根據(jù)標(biāo)準(zhǔn)化振型Φ求解各樓層最大位移Δi,max,計(jì)算樓層最大位移角θi,max.

    通過(guò)上述介紹,基于Info-Gap理論的結(jié)構(gòu)抗震性能穩(wěn)健性優(yōu)化設(shè)計(jì)流程如圖2所示.

    3算例分析

    3.1模型及算法描述

    本文算例為一個(gè)6層3跨鋼框架結(jié)構(gòu)的抗震穩(wěn)健性優(yōu)化設(shè)計(jì).設(shè)計(jì)資料如下:跨度7.5 m,首層層高4.5 m,其余層高3.9 m,抗震設(shè)防烈度7度,場(chǎng)地類別一組三類場(chǎng)地.在本文優(yōu)化中首先對(duì)梁柱進(jìn)行編組,柱子分組為C1~C12,梁分組為B1~B6,單榀框架模型如圖3所示.

    采用本文算法在滿足最大層間位移角小于2.0%前提下,實(shí)現(xiàn)鋼框架總質(zhì)量最小并且抵抗的地震影響系數(shù)最大.本文優(yōu)化部分采用MATLAB嵌套DRAIN編程進(jìn)行結(jié)構(gòu)抗震性能優(yōu)化設(shè)計(jì),在設(shè)計(jì)中將3D3S的H型鋼庫(kù)截面尺寸定義為設(shè)計(jì)變量并對(duì)其進(jìn)行編碼以便在優(yōu)化設(shè)計(jì)中調(diào)用.GA算法參數(shù)選取如下:群體規(guī)模pop=50;最大迭代進(jìn)化代數(shù)gen=150;因?yàn)樗阉骺臻g是離散變量,所以采用正整數(shù)編碼,模擬二進(jìn)制編碼;選擇方法采用輪盤賭法;交叉方法采用多點(diǎn)交叉,交叉概率pc=50%;變異方法采用非均勻變異,變異概率pm=10%.約束采用限定搜索空間與罰函數(shù)相結(jié)合的方法,首先限定解搜索空間為柱選用寬翼緣H型鋼,梁選用窄翼緣H型鋼.分析程序的參數(shù)如下:以位移控制,每次位移增加push step=2.5 cm;材料參數(shù):屈服后的強(qiáng)度比SHR=0.03;屈服強(qiáng)度Fy=398 MPa;考慮P-Δ效應(yīng).

    3.2考慮αmax不確定的穩(wěn)健性抗震優(yōu)化設(shè)計(jì)

    將αmax作為不確定變量,性能函數(shù)R(X,u)為結(jié)構(gòu)層間最大位移角θi,max.將水平地震影響系數(shù)最大值αmax對(duì)于名義值max 所增大的百分比定義為模型的穩(wěn)健性,數(shù)學(xué)表達(dá)式如下:

    Πα,max =αmax :αmax -max max ≤α,α≥0.(6)

    對(duì)應(yīng)地將式(6)中α定義為穩(wěn)健性指標(biāo),框架總質(zhì)量隨αmax穩(wěn)健性指標(biāo)變化,如圖4所示.

    將圖4對(duì)應(yīng)①~⑥不同穩(wěn)健性指標(biāo)的反映頂點(diǎn)位移角和側(cè)向力關(guān)系的pushover曲線列于圖5(a)~5(f) 中,并將對(duì)應(yīng)的框架整體特性列于表1中(其中T為結(jié)構(gòu)的自振周期,M為結(jié)構(gòu)的質(zhì)量,θi,max為pushover分析中結(jié)構(gòu)最大層間位移值).

    由表1可知,反映αmax穩(wěn)健性指標(biāo)越高所需要鋼材質(zhì)量越大,穩(wěn)健性指標(biāo)為30%時(shí),需要的鋼材質(zhì)量要增加17.68%,穩(wěn)健性指標(biāo)為80%時(shí),需要增加的鋼材質(zhì)量為40.05%;結(jié)構(gòu)自振周期隨αmax穩(wěn)健性的增高而減小,這也說(shuō)明隨著αmax的增大可通過(guò)增大結(jié)構(gòu)的剛度抵抗地震作用.

    3.3考慮Tg不確定的穩(wěn)健性抗震優(yōu)化設(shè)計(jì)

    將Tg作為不確定變量,性能函數(shù)定義與3.2節(jié)中相同.將Tg對(duì)于名義值g波動(dòng)幅度值(百分比)定義為穩(wěn)健性指標(biāo)α,數(shù)學(xué)表達(dá)式為:

    Πα,g=Tg:Tg-gg≤α,α≥0. (7)

    框架總質(zhì)量隨Tg穩(wěn)健性指標(biāo)變化,如圖6所示.

    將圖6中對(duì)應(yīng)①~⑥不同穩(wěn)性指標(biāo)的反映頂點(diǎn)位移角和側(cè)向力關(guān)系的pushover曲線列于圖7(a)~7(f) 中,并將對(duì)應(yīng)的框架整體特性列于表2中(其中T為結(jié)構(gòu)的自振周期,M為結(jié)構(gòu)的質(zhì)量,θi,max為pushover分析中結(jié)構(gòu)最大層間位移值).

    由表2可知,場(chǎng)地周期不確定性的穩(wěn)健性指標(biāo)越大,需要的結(jié)構(gòu)質(zhì)量越大,說(shuō)明要抵抗場(chǎng)地周期不確定需要更大的結(jié)構(gòu)質(zhì)量,如穩(wěn)健性指標(biāo)為40%時(shí),需要增大6.29%的結(jié)構(gòu)質(zhì)量,而當(dāng)穩(wěn)健性指標(biāo)為80%時(shí),則需要增大12.49%的結(jié)構(gòu)質(zhì)量.

    3.4兩種不確定性穩(wěn)健性優(yōu)化對(duì)比分析

    圖8為兩種穩(wěn)健性優(yōu)化結(jié)果對(duì)比圖.由圖8可知,同等地震作用條件下穩(wěn)健性指標(biāo)為0時(shí),兩點(diǎn)基本重合,隨穩(wěn)健性指標(biāo)的增大兩條曲線逐漸分離開(kāi)來(lái),其中地震影響系數(shù)的穩(wěn)健性指標(biāo)增大引起優(yōu)化質(zhì)量增加得更快;從圖8中還可以看出,地震影響系數(shù)不確定對(duì)結(jié)構(gòu)質(zhì)量的影響相對(duì)于場(chǎng)地周期不確定性對(duì)結(jié)構(gòu)質(zhì)量的影響要大得多.

    4結(jié)論

    本文將Info-Gap 決策理論引入基于性能的結(jié)構(gòu)抗震優(yōu)化設(shè)計(jì)中,建立了基于Info-Gap 決策的結(jié)構(gòu)穩(wěn)健設(shè)計(jì)優(yōu)化模型,采用GA實(shí)現(xiàn)該模型的嵌套優(yōu)化問(wèn)題.分別就考慮水平地震影響系數(shù)最大值αmax和場(chǎng)地周期Tg不確定性進(jìn)行了鋼框架穩(wěn)健優(yōu)化設(shè)計(jì)分析,通過(guò)對(duì)6層3跨框架分析表明:反映水平地震影響系數(shù)最大值αmax和場(chǎng)地周期Tg不確定性的穩(wěn)健性指標(biāo)越大所需要的鋼材質(zhì)量越大.結(jié)構(gòu)設(shè)計(jì)師和業(yè)主可根據(jù)自己的要求選擇合適的構(gòu)件組合,付出的代價(jià)就是要通過(guò)增加結(jié)構(gòu)的質(zhì)量抵抗αmax不確定和Tg不確定性.通過(guò)兩種穩(wěn)健性設(shè)計(jì)的對(duì)比可以看出,αmax不確定對(duì)結(jié)構(gòu)質(zhì)量的影響相對(duì)于Tg不確定性對(duì)結(jié)構(gòu)質(zhì)量的影響要大得多.

    參考文獻(xiàn)

    [1]GB 50011—2010建筑抗震設(shè)計(jì)規(guī)范[S]. 北京:中國(guó)建筑工業(yè)出版社,2010:33-35.

    GB 50011—2010Code for seismic design of building[S]. Beijing: China Architecture & Building Press, 2010:33-35. (In Chinese)

    [2]BEN-HAIM Y, ELISHAKOFF I. Convex models of uncertainty in applied mechanics [M]. Amsterdam: Elsevier, 1990:325-330.

    [3]BEN-HAIM ?Y. Convex models of uncertainty: applications and implications[J]. Erkenntnis, 1994, 41(2): 139-156.

    [4]MOENS D, VANDEPITTE D. A survey of non-probabilistic uncertainty treatment in finite element analysis[J]. Computer Methods in Applied Mechanics and Engineering, 2005, 194(12): 1527-1555.

    [5]TAKEWAKI I, BEN-HAIM ?Y. Info-Gap robust design with load and model uncertainties[J]. Journal of Sound and Vibration, 2005, 288(3): 551-570.

    [6]ELISHAKOFF I, ELISSEEFF P, GLEGG S A. Nonprobabilistic, convex-theoretic modeling of scatter in material properties[J]. AIAA Journal, 1994, 32(4): 843-849.

    [7]BEN-HAIM Y. A non-probabilistic concept of reliability[J]. Structural Safety, 1994, 14(4):227-245.

    [8]亢戰(zhàn),羅陽(yáng)軍. 基于凸模型的結(jié)構(gòu)非概率可靠性優(yōu)化[J]. 力學(xué)學(xué)報(bào), 2006, 38(6): 807-815.

    KANG Zhan, LUO Yang-jun. On structural optimization for non-probabilistic reliability based on convex models [J]. Acta Mechanica Sinica:Chinese Edition, 2006, 38(6): 807-815. (In Chinese)

    [9]羅陽(yáng)軍,亢戰(zhàn). 超橢球模型下結(jié)構(gòu)非概率可靠性指標(biāo)的迭代算法[J]. 計(jì)算力學(xué)學(xué)報(bào), 2009, 25(6): 747-752.

    LUO Yang-jun, KANG Zhan. An iteration approach for structural non-probabilistic reliability analysis based on hyper-ellipsoidal models[J]. Chinese Journal of Computational Mechanics, 2009,25(6):747-752. (In Chinese)

    [10]蘇永華, 常偉濤, 趙明華. 深部巷道圍巖穩(wěn)定的區(qū)間非概率指標(biāo)分析[J]. 湖南大學(xué)學(xué)報(bào):自然科學(xué)版, 2007, 34(7): 17-21.

    SU Yong-hua, CHANG Wei-tao, ZHAO Ming-hua. Non-probabilistic index analysis method of the stability of deep rock around roadway based on interval variable[J]. Journal of Hunan University:Natural Sciences, 2007, 34(7): 17-21. (In Chinese)

    [11]BEN-HAIM ?Y. Info-Gap decision theory: decisions under severe uncertainty[M]. 2nd ed. London: Academic Press, 2006: 260.

    [12]PIERCE S G, WORDEN K, MANSON G. A novel information-Gap technique to assess reliability of neural network-based damage detection[J]. Journal of Sound and Vibration, 2006, 293(1): 96-111.

    [13]MOILANMEN A, WINTE B A. Uncertainty analysis favours selection of spatially aggregated reserve networks[J]. Biological Conservation, 2006, 129(3): 427-434.

    [14]DUNCAN S J, BRAS B, PAREDIS C J. An approach to robust decision making under severe uncertainty in life cycle design[J]. International Journal of Sustainable Design, 2008, 1(1): 45-59.

    [15]BEN-HAIM ?Y. Value-at-risk with info-Gap uncertainty[J]. Journal of Risk Finance, 2005, 6(5): 388-403.

    [16]BEN-HAIM ?Y. Uncertainty, probability and information-Gaps[J]. Reliability Engineering & System Safety, 2004, 85(1): 249-266.

    [17]RAO S S, BERKE L. Analysis of uncertain structural systems using interval analysis[J]. AIAA Journal, 1997, 35(4): 727-735.

    [18]DAVIS L. Handbook of genetic algorithms[M].New York:Van Nostrand Reinhold, 1991: 325-330.

    [19]劉思思,趙明華,楊明輝, 等. 基于自組織神經(jīng)網(wǎng)絡(luò)與遺傳算法的邊坡穩(wěn)定性分析方法[J]. 湖南大學(xué)學(xué)報(bào):自然科學(xué)版, 2008, 35 (12): 7-12.

    LIU Si-si, ZHAO Ming-hua, YANG Ming-hui,et al. Slope stability analysis method based on self-organizing neural network and genetic algorithm[J]. Journal of Hunan University:Natural Sciences, 2008, 35(12): 7-12. (In Chinese)

    [20]LEITE J P, TOPPING B H. Improved genetic operators for structural engineering optimization[J]. Advances in Engineering Software, 1998, 29(7): 529-562.

    [21]AABEK O, DION C M, YEDDER A B H. Evolutionary algorithms for the optimal laser control of molecular orientation[J]. Journal of Physics B: Atomic, Molecular and Optical Physics, 2003, 36(23): 4667.

    [22]BRAMLTTE M F, CUSIC R. A comparative evaluation of search methods applied to parametric design of aircraft[C]//Proceedings of 3th International Conference on Genetic Algorithms. Burlington,Massachusetts: Morgan Kaufmann Publishers Inc,1989: 213-218.

    [23]呂西林,周定松. 考慮場(chǎng)地類別與設(shè)計(jì)分組的延性需求譜和彈塑性位移反應(yīng)譜[J]. 地震工程與工程振動(dòng), 2004, 24(1): 39-48.

    LV Xi-lin, ZHOU Ding-song. Ductility demand spectra and inelastic displacement spectra considering soil conditions and design characteristic periods[J]. Earthquake Engineering and Engineering Vibration, 2004, 24(1): 39-48. (In Chinese)

    猜你喜歡
    穩(wěn)健性抗震設(shè)計(jì)
    會(huì)計(jì)穩(wěn)健性的定義和計(jì)量
    會(huì)計(jì)穩(wěn)健性的研究評(píng)述
    會(huì)計(jì)穩(wěn)健性的文獻(xiàn)綜述
    不確定性、會(huì)計(jì)穩(wěn)健性與投資效率
    會(huì)計(jì)穩(wěn)健性影響因素與經(jīng)濟(jì)結(jié)果分析
    會(huì)計(jì)穩(wěn)健性的經(jīng)濟(jì)后果研究述評(píng)
    高層建筑結(jié)構(gòu)抗震設(shè)計(jì)要點(diǎn)探析
    關(guān)于建筑結(jié)構(gòu)設(shè)計(jì)中的抗震設(shè)計(jì)探析
    高層建筑結(jié)構(gòu)抗震設(shè)計(jì)要點(diǎn)探析
    樓梯間的抗震設(shè)計(jì)問(wèn)題
    久久亚洲精品不卡| 午夜福利视频1000在线观看 | 国产亚洲精品一区二区www| 侵犯人妻中文字幕一二三四区| 天天躁夜夜躁狠狠躁躁| 少妇粗大呻吟视频| 老汉色∧v一级毛片| 欧美 亚洲 国产 日韩一| 91老司机精品| 亚洲精品国产色婷婷电影| 9色porny在线观看| 欧美成人午夜精品| 男女床上黄色一级片免费看| 给我免费播放毛片高清在线观看| 国产在线精品亚洲第一网站| 一级a爱视频在线免费观看| 一级a爱片免费观看的视频| 成人亚洲精品一区在线观看| 国产极品粉嫩免费观看在线| 亚洲av片天天在线观看| 亚洲少妇的诱惑av| 国产亚洲精品久久久久5区| 精品无人区乱码1区二区| 国内精品久久久久精免费| 侵犯人妻中文字幕一二三四区| 精品一品国产午夜福利视频| 午夜福利在线观看吧| 欧美另类亚洲清纯唯美| 一级黄色大片毛片| 精品人妻在线不人妻| 日本a在线网址| 免费在线观看亚洲国产| 国产免费男女视频| 涩涩av久久男人的天堂| а√天堂www在线а√下载| 欧美色欧美亚洲另类二区 | 久久久精品欧美日韩精品| 99精品久久久久人妻精品| 久久国产精品人妻蜜桃| 少妇粗大呻吟视频| 俄罗斯特黄特色一大片| 色综合欧美亚洲国产小说| 成人永久免费在线观看视频| 美女免费视频网站| 麻豆av在线久日| 日本在线视频免费播放| 亚洲 欧美 日韩 在线 免费| 19禁男女啪啪无遮挡网站| 精品欧美国产一区二区三| 黑人欧美特级aaaaaa片| 99久久久亚洲精品蜜臀av| 人人妻人人爽人人添夜夜欢视频| 人人妻人人澡人人看| 国产精品野战在线观看| 91成人精品电影| 国产精品自产拍在线观看55亚洲| 成年女人毛片免费观看观看9| 美国免费a级毛片| 国产男靠女视频免费网站| 久久久久久大精品| 97超级碰碰碰精品色视频在线观看| 中文字幕高清在线视频| 波多野结衣av一区二区av| 欧美日韩乱码在线| 欧美中文综合在线视频| 精品一区二区三区视频在线观看免费| 亚洲熟妇熟女久久| 免费不卡黄色视频| 曰老女人黄片| 国产精品,欧美在线| 十八禁人妻一区二区| 正在播放国产对白刺激| 欧美亚洲日本最大视频资源| 久久久国产精品麻豆| 无人区码免费观看不卡| 午夜日韩欧美国产| 亚洲伊人色综图| 国产成人精品在线电影| 亚洲第一欧美日韩一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 亚洲视频免费观看视频| 99香蕉大伊视频| 国产精品电影一区二区三区| 亚洲美女黄片视频| 啪啪无遮挡十八禁网站| 久久久久精品国产欧美久久久| 少妇 在线观看| 国产精品久久久久久精品电影 | 高清在线国产一区| 一级黄色大片毛片| 国产精品日韩av在线免费观看 | 日韩精品免费视频一区二区三区| 这个男人来自地球电影免费观看| 看免费av毛片| 欧美av亚洲av综合av国产av| 久久久久国内视频| 正在播放国产对白刺激| 日本免费一区二区三区高清不卡 | 9色porny在线观看| 精品国产超薄肉色丝袜足j| 久久婷婷人人爽人人干人人爱 | 欧美不卡视频在线免费观看 | 亚洲自拍偷在线| 亚洲欧美一区二区三区黑人| 亚洲天堂国产精品一区在线| 老熟妇仑乱视频hdxx| 国产一区二区激情短视频| 国产成人影院久久av| 午夜老司机福利片| 久久久久国产一级毛片高清牌| 国产激情久久老熟女| 一卡2卡三卡四卡精品乱码亚洲| 国内久久婷婷六月综合欲色啪| 成人18禁高潮啪啪吃奶动态图| 久久 成人 亚洲| 久久精品国产综合久久久| 日韩精品中文字幕看吧| 一区二区三区激情视频| 久久影院123| 50天的宝宝边吃奶边哭怎么回事| 国语自产精品视频在线第100页| 久久香蕉精品热| 午夜影院日韩av| 日本a在线网址| 久久久久精品国产欧美久久久| 国产99白浆流出| 亚洲精品国产精品久久久不卡| 美女 人体艺术 gogo| 中文字幕精品免费在线观看视频| 亚洲视频免费观看视频| 久99久视频精品免费| 午夜日韩欧美国产| 丰满的人妻完整版| 久久久国产成人免费| 9191精品国产免费久久| 婷婷六月久久综合丁香| 精品熟女少妇八av免费久了| 少妇熟女aⅴ在线视频| 一进一出好大好爽视频| 97超级碰碰碰精品色视频在线观看| 亚洲国产精品合色在线| 亚洲av成人av| 很黄的视频免费| 少妇熟女aⅴ在线视频| 国产精品久久久人人做人人爽| 男人操女人黄网站| 搞女人的毛片| 亚洲av成人av| 美女扒开内裤让男人捅视频| 少妇粗大呻吟视频| av网站免费在线观看视频| 欧美久久黑人一区二区| 无遮挡黄片免费观看| 一级片免费观看大全| 91老司机精品| 中文字幕人成人乱码亚洲影| 69精品国产乱码久久久| 岛国在线观看网站| 亚洲精品美女久久av网站| 精品欧美国产一区二区三| 国产亚洲精品av在线| 国产又爽黄色视频| 亚洲人成电影免费在线| 曰老女人黄片| or卡值多少钱| 国产aⅴ精品一区二区三区波| 伊人久久大香线蕉亚洲五| 亚洲精品中文字幕一二三四区| www.熟女人妻精品国产| 一进一出抽搐gif免费好疼| 国产精品精品国产色婷婷| 精品久久久久久成人av| 欧美一区二区精品小视频在线| 国产伦人伦偷精品视频| 夜夜躁狠狠躁天天躁| 人人妻,人人澡人人爽秒播| 国产高清有码在线观看视频 | 色精品久久人妻99蜜桃| 麻豆一二三区av精品| 搡老岳熟女国产| 亚洲五月婷婷丁香| 日韩国内少妇激情av| 午夜免费观看网址| 女人被狂操c到高潮| 国产视频一区二区在线看| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲情色 制服丝袜| 亚洲人成网站在线播放欧美日韩| 国产欧美日韩一区二区三| 高清黄色对白视频在线免费看| 69av精品久久久久久| 在线观看午夜福利视频| 国产成人欧美在线观看| 久久热在线av| 久久精品亚洲熟妇少妇任你| 亚洲成人免费电影在线观看| 97人妻精品一区二区三区麻豆 | 久久久久久人人人人人| avwww免费| 国内毛片毛片毛片毛片毛片| 99riav亚洲国产免费| 国产成人欧美在线观看| 午夜日韩欧美国产| 老司机靠b影院| 男女床上黄色一级片免费看| 夜夜爽天天搞| 夜夜躁狠狠躁天天躁| 成人三级黄色视频| 欧美日韩福利视频一区二区| 一级a爱视频在线免费观看| 成人免费观看视频高清| 久久影院123| 国产亚洲精品第一综合不卡| 国产av又大| 在线十欧美十亚洲十日本专区| 88av欧美| 国产精品久久电影中文字幕| 久久久国产欧美日韩av| av福利片在线| 精品国产亚洲在线| 日本 欧美在线| 国产1区2区3区精品| 男女下面进入的视频免费午夜 | 成人亚洲精品av一区二区| 免费在线观看影片大全网站| 国产精品久久视频播放| 国产熟女xx| 亚洲av美国av| 久久久久精品国产欧美久久久| 人人妻,人人澡人人爽秒播| a在线观看视频网站| 国产精品亚洲美女久久久| 亚洲第一青青草原| 亚洲avbb在线观看| 午夜福利高清视频| a在线观看视频网站| 无限看片的www在线观看| 一二三四社区在线视频社区8| 成人亚洲精品av一区二区| 亚洲专区中文字幕在线| 禁无遮挡网站| 少妇 在线观看| 欧美日韩亚洲综合一区二区三区_| 久久影院123| 免费少妇av软件| 午夜精品久久久久久毛片777| 亚洲一区中文字幕在线| 久久中文字幕一级| 午夜日韩欧美国产| 免费在线观看完整版高清| 性欧美人与动物交配| 精品一区二区三区视频在线观看免费| 1024视频免费在线观看| 女性被躁到高潮视频| 久久久水蜜桃国产精品网| 国产激情久久老熟女| 亚洲精品在线观看二区| 午夜精品久久久久久毛片777| 黄色成人免费大全| 欧美激情高清一区二区三区| 一级毛片高清免费大全| 天天一区二区日本电影三级 | 9191精品国产免费久久| 国产亚洲欧美在线一区二区| 精品国内亚洲2022精品成人| 久久精品国产清高在天天线| 亚洲国产中文字幕在线视频| 一级a爱片免费观看的视频| 国产精品美女特级片免费视频播放器 | 亚洲无线在线观看| 久久国产精品人妻蜜桃| 国产欧美日韩一区二区三| 久久久久国产一级毛片高清牌| 男人舔女人的私密视频| 真人一进一出gif抽搐免费| 国产精品国产高清国产av| 亚洲av电影在线进入| xxx96com| 国产91精品成人一区二区三区| 成人国语在线视频| 神马国产精品三级电影在线观看 | 老鸭窝网址在线观看| 精品久久久精品久久久| 久久精品人人爽人人爽视色| 免费观看精品视频网站| 黄片小视频在线播放| 精品午夜福利视频在线观看一区| 老司机靠b影院| 丝袜人妻中文字幕| 国产精品一区二区精品视频观看| 色播亚洲综合网| 成人永久免费在线观看视频| 露出奶头的视频| 久久热在线av| 午夜福利18| av超薄肉色丝袜交足视频| 9色porny在线观看| 成人av一区二区三区在线看| 亚洲,欧美精品.| 久久午夜亚洲精品久久| 脱女人内裤的视频| 亚洲专区中文字幕在线| 国产av精品麻豆| 三级毛片av免费| 国产99久久九九免费精品| 亚洲欧美日韩另类电影网站| 亚洲国产毛片av蜜桃av| 国产成人精品久久二区二区91| 熟女少妇亚洲综合色aaa.| 99国产精品一区二区三区| 美女免费视频网站| 亚洲黑人精品在线| 成人亚洲精品av一区二区| 777久久人妻少妇嫩草av网站| 美女免费视频网站| 男人的好看免费观看在线视频 | 久久久久久久精品吃奶| 99久久99久久久精品蜜桃| 在线天堂中文资源库| 国产aⅴ精品一区二区三区波| 欧美黑人精品巨大| 免费久久久久久久精品成人欧美视频| 国产真人三级小视频在线观看| 国产亚洲欧美98| 亚洲 国产 在线| 99精品久久久久人妻精品| 妹子高潮喷水视频| 国产欧美日韩一区二区精品| 国产一级毛片七仙女欲春2 | 可以在线观看毛片的网站| 久久中文看片网| svipshipincom国产片| 国产精品久久久av美女十八| 国产成人欧美在线观看| 丁香欧美五月| 成人亚洲精品av一区二区| 97人妻精品一区二区三区麻豆 | 制服人妻中文乱码| 欧美黄色片欧美黄色片| 欧美一级a爱片免费观看看 | 在线永久观看黄色视频| av天堂久久9| 国产精品久久久人人做人人爽| 国产精品久久电影中文字幕| 国产av一区二区精品久久| 日韩欧美在线二视频| 国产精品香港三级国产av潘金莲| 国产精品 欧美亚洲| 亚洲片人在线观看| 日韩欧美在线二视频| 国产精品香港三级国产av潘金莲| √禁漫天堂资源中文www| 美女午夜性视频免费| 性欧美人与动物交配| 男女午夜视频在线观看| 欧美黄色淫秽网站| 国产真人三级小视频在线观看| 亚洲第一欧美日韩一区二区三区| 午夜视频精品福利| 国产精品久久久久久人妻精品电影| 91精品国产国语对白视频| 国产xxxxx性猛交| 国产av又大| 超碰成人久久| 88av欧美| 欧美日韩亚洲综合一区二区三区_| 国产高清videossex| 亚洲伊人色综图| cao死你这个sao货| 成年人黄色毛片网站| 欧美绝顶高潮抽搐喷水| 亚洲成人国产一区在线观看| 性少妇av在线| 亚洲人成网站在线播放欧美日韩| 99精品欧美一区二区三区四区| 18禁黄网站禁片午夜丰满| 一本大道久久a久久精品| 巨乳人妻的诱惑在线观看| 午夜老司机福利片| 国产亚洲av高清不卡| 午夜福利免费观看在线| 亚洲精品国产一区二区精华液| 黄色 视频免费看| 亚洲性夜色夜夜综合| 国产精品一区二区在线不卡| 亚洲欧美激情在线| 国产精品香港三级国产av潘金莲| 精品欧美国产一区二区三| 国产精品一区二区免费欧美| 亚洲成人免费电影在线观看| 精品国产乱码久久久久久男人| 久久热在线av| 国产在线观看jvid| 成人三级做爰电影| 精品人妻在线不人妻| 亚洲精品国产色婷婷电影| 1024香蕉在线观看| 欧美最黄视频在线播放免费| 琪琪午夜伦伦电影理论片6080| 国产1区2区3区精品| 露出奶头的视频| 美国免费a级毛片| 一夜夜www| 色婷婷久久久亚洲欧美| 麻豆av在线久日| 黄网站色视频无遮挡免费观看| 欧美成狂野欧美在线观看| 老司机靠b影院| 人人妻人人澡人人看| av天堂在线播放| 国产成人欧美在线观看| 久久精品aⅴ一区二区三区四区| 亚洲精品国产色婷婷电影| 99热只有精品国产| 亚洲色图综合在线观看| 亚洲人成77777在线视频| 国产精品亚洲美女久久久| 日日干狠狠操夜夜爽| 国产乱人伦免费视频| a在线观看视频网站| 身体一侧抽搐| 热99re8久久精品国产| 69精品国产乱码久久久| www.自偷自拍.com| 精品久久蜜臀av无| 欧美午夜高清在线| 国产激情欧美一区二区| 久久久精品欧美日韩精品| 欧美绝顶高潮抽搐喷水| 麻豆av在线久日| 999久久久精品免费观看国产| 好男人在线观看高清免费视频 | 19禁男女啪啪无遮挡网站| 亚洲av片天天在线观看| 国产高清有码在线观看视频 | 级片在线观看| 亚洲av成人av| 精品久久久精品久久久| 一级片免费观看大全| 夜夜看夜夜爽夜夜摸| 波多野结衣一区麻豆| 黑人巨大精品欧美一区二区蜜桃| 精品午夜福利视频在线观看一区| av视频免费观看在线观看| 两个人看的免费小视频| av有码第一页| 欧美乱色亚洲激情| 电影成人av| 国产午夜福利久久久久久| bbb黄色大片| 男女午夜视频在线观看| 欧美国产日韩亚洲一区| 亚洲久久久国产精品| 91精品国产国语对白视频| 欧美成人免费av一区二区三区| 精品人妻1区二区| 精品午夜福利视频在线观看一区| 久久中文看片网| 国产1区2区3区精品| 中文字幕另类日韩欧美亚洲嫩草| 9热在线视频观看99| 久久久久九九精品影院| 老汉色∧v一级毛片| 极品人妻少妇av视频| 色在线成人网| 中文字幕最新亚洲高清| videosex国产| www日本在线高清视频| 十八禁人妻一区二区| 午夜精品国产一区二区电影| 久久久久国产一级毛片高清牌| 亚洲视频免费观看视频| 精品高清国产在线一区| 亚洲国产精品成人综合色| 少妇熟女aⅴ在线视频| 亚洲人成电影观看| 久久久久亚洲av毛片大全| 亚洲成国产人片在线观看| 午夜福利一区二区在线看| 高潮久久久久久久久久久不卡| 黑丝袜美女国产一区| 少妇熟女aⅴ在线视频| 黄片大片在线免费观看| 韩国精品一区二区三区| 国产成人啪精品午夜网站| 亚洲成av片中文字幕在线观看| 黑人操中国人逼视频| 久久亚洲真实| 亚洲成人久久性| 久99久视频精品免费| 好男人电影高清在线观看| 国产国语露脸激情在线看| 少妇粗大呻吟视频| 99精品久久久久人妻精品| 精品国产超薄肉色丝袜足j| 1024视频免费在线观看| 亚洲精品一卡2卡三卡4卡5卡| 一夜夜www| 国产伦一二天堂av在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 91老司机精品| 亚洲,欧美精品.| 国产日韩一区二区三区精品不卡| 久久精品国产亚洲av香蕉五月| 大码成人一级视频| 免费看a级黄色片| 桃色一区二区三区在线观看| 最好的美女福利视频网| 国内精品久久久久久久电影| www.自偷自拍.com| 国产精品秋霞免费鲁丝片| 成人av一区二区三区在线看| 麻豆国产av国片精品| 在线观看www视频免费| 日韩高清综合在线| 一边摸一边做爽爽视频免费| 久久久久九九精品影院| 国产亚洲av高清不卡| 高清毛片免费观看视频网站| 麻豆成人av在线观看| 黄片大片在线免费观看| 黄网站色视频无遮挡免费观看| 免费在线观看黄色视频的| 免费少妇av软件| 亚洲欧美日韩高清在线视频| 亚洲七黄色美女视频| 亚洲电影在线观看av| 97超级碰碰碰精品色视频在线观看| 他把我摸到了高潮在线观看| 成人精品一区二区免费| 色精品久久人妻99蜜桃| 亚洲伊人色综图| 国产高清videossex| 露出奶头的视频| 少妇 在线观看| 很黄的视频免费| 国产欧美日韩一区二区三| 色播在线永久视频| 99riav亚洲国产免费| 亚洲精品一区av在线观看| 成人18禁高潮啪啪吃奶动态图| 国产av精品麻豆| 中文亚洲av片在线观看爽| 免费在线观看亚洲国产| 国产精品,欧美在线| 亚洲一区二区三区不卡视频| 久久九九热精品免费| 精品电影一区二区在线| 国内久久婷婷六月综合欲色啪| 国产精品爽爽va在线观看网站 | 999久久久国产精品视频| www.精华液| 午夜福利免费观看在线| 在线观看免费视频日本深夜| 夜夜夜夜夜久久久久| 美女高潮喷水抽搐中文字幕| 亚洲专区字幕在线| 欧美色视频一区免费| 男人舔女人下体高潮全视频| 波多野结衣巨乳人妻| 麻豆av在线久日| 人人妻人人澡欧美一区二区 | 亚洲狠狠婷婷综合久久图片| 亚洲成人免费电影在线观看| 国产成人欧美在线观看| 在线天堂中文资源库| 国产精品98久久久久久宅男小说| 日韩一卡2卡3卡4卡2021年| 欧美激情 高清一区二区三区| 熟妇人妻久久中文字幕3abv| 69精品国产乱码久久久| 婷婷六月久久综合丁香| 91九色精品人成在线观看| 久久九九热精品免费| 精品国产一区二区三区四区第35| 国产1区2区3区精品| 一边摸一边抽搐一进一小说| 女人高潮潮喷娇喘18禁视频| 国产成人av激情在线播放| 亚洲aⅴ乱码一区二区在线播放 | 午夜福利成人在线免费观看| 午夜福利欧美成人| 国产亚洲av高清不卡| 在线观看日韩欧美| 日韩精品青青久久久久久| 性少妇av在线| 亚洲国产看品久久| 久久人妻av系列| 91成年电影在线观看| 欧美乱码精品一区二区三区| 国产成人精品无人区| 国语自产精品视频在线第100页| www.www免费av| 少妇粗大呻吟视频| 成人国语在线视频| 欧美黑人精品巨大| 久久青草综合色| 欧美日韩亚洲综合一区二区三区_| 国产精品一区二区精品视频观看| 久久久国产成人免费| 日本五十路高清| 精品日产1卡2卡| 亚洲,欧美精品.| 国产成人免费无遮挡视频| 国产成人精品在线电影| 久久精品亚洲精品国产色婷小说| 一个人观看的视频www高清免费观看 | 午夜成年电影在线免费观看| 国产精品,欧美在线| 99国产精品99久久久久| 最近最新中文字幕大全免费视频| 亚洲欧美激情在线| 美女 人体艺术 gogo| 国产精品香港三级国产av潘金莲| 我的亚洲天堂|