• <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)題
    飞空精品影院首页| 日韩 欧美 亚洲 中文字幕| 精品人妻一区二区三区麻豆| 女人高潮潮喷娇喘18禁视频| 国产xxxxx性猛交| 成年动漫av网址| 考比视频在线观看| 宅男免费午夜| 五月开心婷婷网| 免费女性裸体啪啪无遮挡网站| 一级片免费观看大全| 日韩一本色道免费dvd| 亚洲精品久久久久久婷婷小说| 国产精品国产三级专区第一集| 99精品久久久久人妻精品| 中国三级夫妇交换| 男女国产视频网站| 久久97久久精品| 日本欧美国产在线视频| 日韩一区二区三区影片| 熟妇人妻不卡中文字幕| 一级爰片在线观看| a级片在线免费高清观看视频| 大香蕉久久网| 99国产综合亚洲精品| 一级片免费观看大全| av福利片在线| 大码成人一级视频| 精品第一国产精品| 欧美乱码精品一区二区三区| 亚洲欧美精品综合一区二区三区| 国产精品久久久久久久久免| 精品亚洲成a人片在线观看| 午夜免费观看性视频| 老汉色av国产亚洲站长工具| 国产熟女欧美一区二区| 亚洲精品自拍成人| 久久久久精品久久久久真实原创| av网站在线播放免费| 亚洲精品久久午夜乱码| 夜夜骑夜夜射夜夜干| 亚洲五月色婷婷综合| 成年av动漫网址| 久久狼人影院| 中文欧美无线码| 日韩伦理黄色片| 亚洲精品一二三| 国产成人精品福利久久| 亚洲美女黄色视频免费看| 精品久久久久久电影网| 欧美少妇被猛烈插入视频| 午夜福利视频精品| 9热在线视频观看99| 18禁观看日本| 午夜免费观看性视频| av有码第一页| 久久久久久久久久久久大奶| 99香蕉大伊视频| 精品久久蜜臀av无| 欧美国产精品一级二级三级| 久久青草综合色| 99久久99久久久精品蜜桃| 欧美精品亚洲一区二区| 老熟女久久久| 国产又色又爽无遮挡免| 色94色欧美一区二区| 交换朋友夫妻互换小说| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日韩亚洲国产一区二区在线观看 | 亚洲成国产人片在线观看| 久久ye,这里只有精品| 嫩草影院入口| 国产精品一二三区在线看| 国产日韩欧美在线精品| 十八禁网站网址无遮挡| 久久99一区二区三区| 国产国语露脸激情在线看| 国产精品三级大全| 午夜日韩欧美国产| 秋霞在线观看毛片| 91成人精品电影| 国产毛片在线视频| 国产成人精品无人区| 久久久精品国产亚洲av高清涩受| 伦理电影免费视频| 夜夜骑夜夜射夜夜干| 黄色 视频免费看| 男女无遮挡免费网站观看| 黄色视频不卡| 国产黄色免费在线视频| 成人国语在线视频| 黄色毛片三级朝国网站| 国产在视频线精品| 亚洲熟女毛片儿| 国产一区有黄有色的免费视频| 午夜老司机福利片| 爱豆传媒免费全集在线观看| av在线播放精品| 欧美日韩av久久| 亚洲av在线观看美女高潮| 午夜免费男女啪啪视频观看| 国产高清不卡午夜福利| 18在线观看网站| 18禁裸乳无遮挡动漫免费视频| 午夜免费鲁丝| 亚洲欧美中文字幕日韩二区| 亚洲少妇的诱惑av| 国产探花极品一区二区| 国产精品 国内视频| 国产成人91sexporn| 搡老乐熟女国产| 2021少妇久久久久久久久久久| 亚洲专区中文字幕在线 | 男女国产视频网站| 国产精品久久久久成人av| 亚洲国产欧美日韩在线播放| 一区二区日韩欧美中文字幕| 成人亚洲精品一区在线观看| 国产精品成人在线| 亚洲欧洲日产国产| 国产在视频线精品| 自线自在国产av| 亚洲,欧美,日韩| 美女大奶头黄色视频| 一级毛片 在线播放| 天美传媒精品一区二区| 日韩 亚洲 欧美在线| 在线观看一区二区三区激情| 五月开心婷婷网| 日本av免费视频播放| 国产精品 欧美亚洲| 日韩免费高清中文字幕av| 天天躁夜夜躁狠狠躁躁| 高清视频免费观看一区二区| 欧美在线一区亚洲| 七月丁香在线播放| 国产1区2区3区精品| 天堂俺去俺来也www色官网| 精品一品国产午夜福利视频| 免费观看人在逋| 日韩不卡一区二区三区视频在线| 黑人欧美特级aaaaaa片| 九色亚洲精品在线播放| tube8黄色片| 午夜福利免费观看在线| 亚洲天堂av无毛| 成人影院久久| 久久久久精品性色| 巨乳人妻的诱惑在线观看| 国产精品一区二区在线观看99| 国产精品欧美亚洲77777| 亚洲成人手机| 亚洲熟女毛片儿| 日韩中文字幕欧美一区二区 | 婷婷成人精品国产| 日韩精品免费视频一区二区三区| 香蕉丝袜av| 精品久久久久久电影网| 国精品久久久久久国模美| 精品国产超薄肉色丝袜足j| 国产av一区二区精品久久| 欧美xxⅹ黑人| 免费不卡黄色视频| 热re99久久国产66热| 飞空精品影院首页| 亚洲,欧美,日韩| 国产极品天堂在线| 黑人猛操日本美女一级片| 免费在线观看视频国产中文字幕亚洲 | 欧美日韩成人在线一区二区| 多毛熟女@视频| 老司机深夜福利视频在线观看 | 久久久久视频综合| 99久久99久久久精品蜜桃| 天天躁夜夜躁狠狠躁躁| 国产在线视频一区二区| 两性夫妻黄色片| 一级毛片电影观看| 777久久人妻少妇嫩草av网站| 天天躁狠狠躁夜夜躁狠狠躁| 午夜福利在线免费观看网站| 欧美日韩一级在线毛片| 国产99久久九九免费精品| 久久久精品区二区三区| 1024香蕉在线观看| 亚洲一区二区三区欧美精品| 精品午夜福利在线看| 精品卡一卡二卡四卡免费| 纯流量卡能插随身wifi吗| 精品少妇黑人巨大在线播放| 久久99精品国语久久久| 国产高清国产精品国产三级| 国产精品一区二区精品视频观看| 日韩电影二区| 欧美在线一区亚洲| 免费久久久久久久精品成人欧美视频| 老司机亚洲免费影院| 无遮挡黄片免费观看| 在线天堂最新版资源| 少妇人妻久久综合中文| 一二三四中文在线观看免费高清| 女的被弄到高潮叫床怎么办| 精品少妇黑人巨大在线播放| 国产精品一国产av| 久久久久人妻精品一区果冻| 亚洲欧美一区二区三区国产| 精品一区二区三卡| 国产女主播在线喷水免费视频网站| 少妇的丰满在线观看| 日韩免费高清中文字幕av| 91成人精品电影| 免费观看人在逋| 成年女人毛片免费观看观看9 | 两性夫妻黄色片| 亚洲国产精品一区三区| 麻豆精品久久久久久蜜桃| 欧美黄色片欧美黄色片| 1024视频免费在线观看| 亚洲免费av在线视频| 日本爱情动作片www.在线观看| 三上悠亚av全集在线观看| 老汉色av国产亚洲站长工具| 波多野结衣一区麻豆| 一区二区三区四区激情视频| 国产无遮挡羞羞视频在线观看| 欧美激情 高清一区二区三区| 免费高清在线观看视频在线观看| 欧美国产精品一级二级三级| 亚洲av电影在线观看一区二区三区| 亚洲成av片中文字幕在线观看| 你懂的网址亚洲精品在线观看| 国产一区二区三区av在线| www.精华液| 在线 av 中文字幕| 在线观看免费视频网站a站| 国产成人精品无人区| 欧美日韩亚洲国产一区二区在线观看 | 亚洲欧美清纯卡通| 国产亚洲av片在线观看秒播厂| 99re6热这里在线精品视频| 国产精品一区二区精品视频观看| 日日啪夜夜爽| 大陆偷拍与自拍| 操出白浆在线播放| 国产日韩一区二区三区精品不卡| 日韩中文字幕欧美一区二区 | 少妇精品久久久久久久| 精品国产露脸久久av麻豆| 国产精品亚洲av一区麻豆 | 美女脱内裤让男人舔精品视频| 中国三级夫妇交换| 欧美国产精品va在线观看不卡| 男人添女人高潮全过程视频| 伊人亚洲综合成人网| 伦理电影免费视频| 晚上一个人看的免费电影| 亚洲成色77777| 可以免费在线观看a视频的电影网站 | 丁香六月天网| 最近中文字幕2019免费版| 亚洲欧美激情在线| 国产成人精品福利久久| 国产精品一区二区在线不卡| av在线老鸭窝| 免费不卡黄色视频| 亚洲精品国产av蜜桃| 久久久亚洲精品成人影院| 一区二区三区乱码不卡18| 久久精品国产亚洲av涩爱| 亚洲国产欧美在线一区| 九九爱精品视频在线观看| 亚洲av国产av综合av卡| 丝袜脚勾引网站| 久久国产精品大桥未久av| 青青草视频在线视频观看| 久久婷婷青草| 成人三级做爰电影| 午夜福利视频在线观看免费| 99热网站在线观看| 18禁动态无遮挡网站| 国产精品 国内视频| 日韩免费高清中文字幕av| 啦啦啦 在线观看视频| svipshipincom国产片| 无遮挡黄片免费观看| 日本一区二区免费在线视频| 悠悠久久av| 一级,二级,三级黄色视频| 一本大道久久a久久精品| 精品酒店卫生间| 国产亚洲午夜精品一区二区久久| 国产日韩欧美亚洲二区| 亚洲美女视频黄频| 九九爱精品视频在线观看| 中文欧美无线码| 九色亚洲精品在线播放| 久久久久久久久久久久大奶| 9热在线视频观看99| 精品国产乱码久久久久久男人| 啦啦啦在线免费观看视频4| a级片在线免费高清观看视频| 一边亲一边摸免费视频| 午夜影院在线不卡| 女人精品久久久久毛片| 日日啪夜夜爽| 操美女的视频在线观看| 亚洲图色成人| 亚洲天堂av无毛| 国产日韩欧美视频二区| 亚洲激情五月婷婷啪啪| 一级a爱视频在线免费观看| 中文字幕人妻丝袜制服| 亚洲综合精品二区| 色网站视频免费| av国产精品久久久久影院| 亚洲精华国产精华液的使用体验| 亚洲国产中文字幕在线视频| 亚洲国产看品久久| 可以免费在线观看a视频的电影网站 | 亚洲精品自拍成人| 777久久人妻少妇嫩草av网站| 男女之事视频高清在线观看 | 色播在线永久视频| 久久国产精品男人的天堂亚洲| 美女大奶头黄色视频| 多毛熟女@视频| 亚洲国产精品成人久久小说| 亚洲精品美女久久久久99蜜臀 | 国产亚洲最大av| 韩国精品一区二区三区| 母亲3免费完整高清在线观看| 色视频在线一区二区三区| 十八禁高潮呻吟视频| 男女边摸边吃奶| 国产免费视频播放在线视频| 欧美在线黄色| 69精品国产乱码久久久| 高清欧美精品videossex| 亚洲精品一二三| 亚洲国产av新网站| 777久久人妻少妇嫩草av网站| 国产成人精品福利久久| 国产男女内射视频| 女人被躁到高潮嗷嗷叫费观| 欧美在线一区亚洲| 卡戴珊不雅视频在线播放| 在线观看免费视频网站a站| 国产极品粉嫩免费观看在线| 欧美日韩av久久| 亚洲av日韩精品久久久久久密 | 99久久精品国产亚洲精品| 国产精品国产三级国产专区5o| 欧美日韩综合久久久久久| 成人18禁高潮啪啪吃奶动态图| 在线观看人妻少妇| 国产精品偷伦视频观看了| 欧美 亚洲 国产 日韩一| 精品人妻熟女毛片av久久网站| 可以免费在线观看a视频的电影网站 | 99国产精品免费福利视频| 国产成人欧美在线观看 | 成人国产av品久久久| 激情视频va一区二区三区| 黄频高清免费视频| 色视频在线一区二区三区| 美女中出高潮动态图| av福利片在线| 侵犯人妻中文字幕一二三四区| 亚洲av电影在线进入| 国产亚洲最大av| 午夜福利视频在线观看免费| 男人操女人黄网站| 亚洲精品日本国产第一区| 国产精品.久久久| 国产片内射在线| 少妇精品久久久久久久| av福利片在线| 一边摸一边做爽爽视频免费| 欧美日韩亚洲综合一区二区三区_| 丰满迷人的少妇在线观看| 欧美精品高潮呻吟av久久| 少妇被粗大猛烈的视频| 一级a爱视频在线免费观看| 国产精品av久久久久免费| 国产精品女同一区二区软件| 亚洲国产看品久久| 精品久久久久久电影网| 亚洲av综合色区一区| 美女大奶头黄色视频| 女人爽到高潮嗷嗷叫在线视频| 欧美黄色片欧美黄色片| 欧美日韩亚洲国产一区二区在线观看 | av福利片在线| 久久人人爽人人片av| 成年美女黄网站色视频大全免费| 久久 成人 亚洲| 国产成人午夜福利电影在线观看| 各种免费的搞黄视频| 亚洲精品国产av蜜桃| 亚洲av中文av极速乱| 亚洲综合精品二区| 国产一区二区激情短视频 | 国产精品 国内视频| 国语对白做爰xxxⅹ性视频网站| 日本黄色日本黄色录像| 精品一区在线观看国产| 最新在线观看一区二区三区 | 欧美日韩视频高清一区二区三区二| 国产淫语在线视频| 日日撸夜夜添| 久久久久久久国产电影| 国产精品国产三级专区第一集| 精品国产国语对白av| 国产成人系列免费观看| a级毛片在线看网站| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产日韩一区二区| 久久影院123| 秋霞伦理黄片| 日韩中文字幕视频在线看片| 久久久久久久国产电影| 国产av码专区亚洲av| 欧美中文综合在线视频| videos熟女内射| 欧美人与性动交α欧美精品济南到| 午夜激情av网站| 久久女婷五月综合色啪小说| 可以免费在线观看a视频的电影网站 | 久久精品aⅴ一区二区三区四区| 亚洲精品久久久久久婷婷小说| 婷婷色综合www| 亚洲成色77777| av天堂久久9| 亚洲国产毛片av蜜桃av| 亚洲精品久久久久久婷婷小说| 超碰成人久久| 免费在线观看完整版高清| 国产男人的电影天堂91| 久久久久久人人人人人| 这个男人来自地球电影免费观看 | 一级毛片黄色毛片免费观看视频| 看免费av毛片| 1024香蕉在线观看| 男女国产视频网站| 搡老岳熟女国产| 国产麻豆69| 亚洲国产欧美日韩在线播放| 亚洲av国产av综合av卡| 亚洲自偷自拍图片 自拍| 岛国毛片在线播放| 久久久久久久久免费视频了| 考比视频在线观看| 又粗又硬又长又爽又黄的视频| 男女午夜视频在线观看| 午夜福利免费观看在线| av电影中文网址| 高清黄色对白视频在线免费看| 777久久人妻少妇嫩草av网站| 国产精品偷伦视频观看了| 久久久久久人人人人人| 国产成人欧美在线观看 | 欧美精品一区二区大全| 妹子高潮喷水视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美色中文字幕在线| 日本猛色少妇xxxxx猛交久久| 黄片小视频在线播放| 高清不卡的av网站| 五月开心婷婷网| xxxhd国产人妻xxx| 国产精品国产av在线观看| 欧美日韩一级在线毛片| 一级黄片播放器| 最近最新中文字幕大全免费视频 | 男男h啪啪无遮挡| 精品少妇黑人巨大在线播放| 亚洲欧洲精品一区二区精品久久久 | 亚洲欧美色中文字幕在线| 色婷婷av一区二区三区视频| 三上悠亚av全集在线观看| 亚洲熟女精品中文字幕| 久久久精品区二区三区| 18在线观看网站| 久热爱精品视频在线9| 国产福利在线免费观看视频| 中文字幕精品免费在线观看视频| videosex国产| 久久国产亚洲av麻豆专区| 欧美黑人欧美精品刺激| 国产一区二区 视频在线| 另类亚洲欧美激情| 男女边吃奶边做爰视频| 日韩一区二区视频免费看| 亚洲av福利一区| 国产一区有黄有色的免费视频| 最新在线观看一区二区三区 | 国产亚洲精品第一综合不卡| 激情五月婷婷亚洲| www.av在线官网国产| 卡戴珊不雅视频在线播放| 国产一区二区激情短视频 | 欧美在线黄色| av片东京热男人的天堂| 在线亚洲精品国产二区图片欧美| 最黄视频免费看| 天天躁夜夜躁狠狠躁躁| 成人国产麻豆网| 男人舔女人的私密视频| 日韩精品有码人妻一区| 嫩草影视91久久| 菩萨蛮人人尽说江南好唐韦庄| av在线app专区| 狂野欧美激情性xxxx| 天天影视国产精品| 欧美日韩亚洲国产一区二区在线观看 | 七月丁香在线播放| 免费黄频网站在线观看国产| 精品午夜福利在线看| 在线观看免费视频网站a站| 免费av中文字幕在线| 国产不卡av网站在线观看| 黑人欧美特级aaaaaa片| 欧美日本中文国产一区发布| 亚洲 欧美一区二区三区| 女人久久www免费人成看片| netflix在线观看网站| 叶爱在线成人免费视频播放| 天天躁日日躁夜夜躁夜夜| 精品国产乱码久久久久久男人| 我要看黄色一级片免费的| 女人久久www免费人成看片| 国产日韩一区二区三区精品不卡| 宅男免费午夜| 日本爱情动作片www.在线观看| 美女高潮到喷水免费观看| 热re99久久国产66热| 丝袜美腿诱惑在线| 女性被躁到高潮视频| 成人毛片60女人毛片免费| 国产精品免费视频内射| 在线观看免费日韩欧美大片| 欧美 日韩 精品 国产| 视频区图区小说| 免费黄频网站在线观看国产| 大片电影免费在线观看免费| 少妇被粗大的猛进出69影院| 国产成人av激情在线播放| 在线天堂中文资源库| 国产免费一区二区三区四区乱码| 狠狠婷婷综合久久久久久88av| 国产片特级美女逼逼视频| 日韩一本色道免费dvd| 精品国产超薄肉色丝袜足j| 少妇被粗大猛烈的视频| 新久久久久国产一级毛片| 午夜激情久久久久久久| 亚洲精品一二三| 黄色毛片三级朝国网站| 国产熟女欧美一区二区| 在现免费观看毛片| 亚洲av欧美aⅴ国产| 激情视频va一区二区三区| 看非洲黑人一级黄片| 91成人精品电影| 成人三级做爰电影| 精品久久蜜臀av无| 国产精品蜜桃在线观看| 十八禁人妻一区二区| 黄色视频在线播放观看不卡| 亚洲欧美一区二区三区黑人| 99久久综合免费| 国产一区二区激情短视频 | 国产精品香港三级国产av潘金莲 | 69精品国产乱码久久久| 激情视频va一区二区三区| 十八禁人妻一区二区| 侵犯人妻中文字幕一二三四区| 国产精品久久久久久人妻精品电影 | 亚洲人成电影观看| 侵犯人妻中文字幕一二三四区| 婷婷色麻豆天堂久久| 999久久久国产精品视频| 日韩中文字幕视频在线看片| 亚洲精品一二三| 丁香六月欧美| 精品酒店卫生间| av在线播放精品| 国产av码专区亚洲av| 免费黄色在线免费观看| 观看av在线不卡| 欧美精品高潮呻吟av久久| 久久久久人妻精品一区果冻| 日韩熟女老妇一区二区性免费视频| 99久国产av精品国产电影| 大香蕉久久成人网| 国产乱人偷精品视频| 黑人猛操日本美女一级片| 岛国毛片在线播放| 老司机亚洲免费影院| 两性夫妻黄色片| 国产亚洲av片在线观看秒播厂| 激情五月婷婷亚洲| 久久精品人人爽人人爽视色| 久久久久精品性色| 日本午夜av视频| 国产一区有黄有色的免费视频| 国产日韩一区二区三区精品不卡| 午夜福利在线免费观看网站| 97人妻天天添夜夜摸| 欧美日韩福利视频一区二区| 日韩视频在线欧美| 综合色丁香网| 黄频高清免费视频| 久久久久久免费高清国产稀缺| 十分钟在线观看高清视频www|