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

    工程結(jié)構(gòu)優(yōu)化設(shè)計(jì)研究進(jìn)展

    2011-04-09 12:11:38李洪煊武穎利
    關(guān)鍵詞:模擬退火遺傳算法變量

    蔡 新,李洪煊,武穎利,朱 杰

    (1.河海大學(xué)力學(xué)與材料學(xué)院,江蘇南京 210098;2.河海大學(xué)水利水電學(xué)院,江蘇南京 210098; 3.南京市市政設(shè)計(jì)研究院,江蘇南京 210008)

    最優(yōu)化就是追求最好的結(jié)果或最優(yōu)的目標(biāo).最優(yōu)化設(shè)計(jì)是從可能設(shè)計(jì)中選擇最合理的設(shè)計(jì)以達(dá)到最優(yōu)目標(biāo).搜尋最優(yōu)設(shè)計(jì)的方法就是最優(yōu)化設(shè)計(jì)法,這種方法的數(shù)學(xué)理論就是最優(yōu)化設(shè)計(jì)理論[1].工程結(jié)構(gòu)優(yōu)化設(shè)計(jì)是把力學(xué)概念和優(yōu)化技術(shù)有機(jī)地結(jié)合,根據(jù)設(shè)計(jì)要求,使部分參與計(jì)算的量以變量出現(xiàn),形成全部可能的結(jié)構(gòu)設(shè)計(jì)方案域,再利用數(shù)學(xué)手段在域中搜索出滿足預(yù)定要求的不僅可行而且最好的設(shè)計(jì)方案.實(shí)踐證明,工程結(jié)構(gòu)優(yōu)化設(shè)計(jì)能縮短設(shè)計(jì)周期,提高設(shè)計(jì)質(zhì)量和水平,較之原設(shè)計(jì)方案,優(yōu)化設(shè)計(jì)方案一般可以降低工程造價(jià)5%~30%[2].

    工程結(jié)構(gòu)優(yōu)化設(shè)計(jì)早期采用的是基于直覺(jué)的準(zhǔn)則法,如滿應(yīng)力準(zhǔn)則法,滿應(yīng)變準(zhǔn)則法等.準(zhǔn)則法的顯著特點(diǎn)是單步設(shè)計(jì)變量修改幅度大,因而收斂速度快,且與結(jié)構(gòu)大小及復(fù)雜程度無(wú)關(guān).

    1960年,Schmit[3]首先將數(shù)學(xué)規(guī)劃法引入到結(jié)構(gòu)優(yōu)化設(shè)計(jì)中.通過(guò)數(shù)學(xué)規(guī)劃法對(duì)一些典型問(wèn)題的研究表明,準(zhǔn)則法求得的解不一定是真實(shí)的最優(yōu)解,其缺乏嚴(yán)格的理論依據(jù).數(shù)學(xué)規(guī)劃法的特點(diǎn)是理論依據(jù)嚴(yán)謹(jǐn),算法的穩(wěn)定性和可靠性較好;求解問(wèn)題的規(guī)模受到限制且要求能顯式表達(dá),計(jì)算效率低.

    實(shí)際的工程結(jié)構(gòu)優(yōu)化設(shè)計(jì)一般都是有約束、非線性和隱式的優(yōu)化問(wèn)題.傳統(tǒng)尋優(yōu)算法中準(zhǔn)則法和數(shù)學(xué)規(guī)劃法由于存在各自的缺陷,不能很好地解決上述復(fù)雜問(wèn)題,因而新的現(xiàn)代尋優(yōu)算法應(yīng)運(yùn)而生.20世紀(jì)70年代,人們?cè)诜治鲅芯孔匀唤缟矬w的進(jìn)化及其進(jìn)化的過(guò)程中,發(fā)現(xiàn)了新的仿生學(xué)方法,主要有基因遺傳算法(genetic algorithm,GA)、神經(jīng)網(wǎng)絡(luò)算法(artificial neural network,ANN)、免疫算法(immunity algorithm,IA)、蟻群算法(ant colony algorithm,ACA)等.在研究退火對(duì)金屬晶體的發(fā)育完善過(guò)程中,產(chǎn)生了模擬自然界物質(zhì)生成的模擬退火算法(simulated annealing,SA).這些新的尋優(yōu)方法具有能夠搜索到全局最優(yōu)解、程序?qū)崿F(xiàn)方便等優(yōu)點(diǎn).

    隨著計(jì)算機(jī)軟、硬件的發(fā)展,基于數(shù)學(xué)力學(xué)方法的結(jié)構(gòu)優(yōu)化設(shè)計(jì)的理論也有較快發(fā)展,由尺寸優(yōu)化,發(fā)展到形狀優(yōu)化、拓?fù)鋬?yōu)化、布局優(yōu)化和類(lèi)型優(yōu)化;由單目標(biāo)優(yōu)化發(fā)展到多目標(biāo)優(yōu)化;由確定性優(yōu)化發(fā)展到不確定性優(yōu)化;由工程結(jié)構(gòu)靜力優(yōu)化發(fā)展到靜動(dòng)力優(yōu)化,優(yōu)化設(shè)計(jì)的層次在不斷提高.

    根據(jù)實(shí)際工程特點(diǎn),對(duì)現(xiàn)有的優(yōu)化算法進(jìn)行改進(jìn)、重組和綜合,從而產(chǎn)生滿足實(shí)際工程結(jié)構(gòu)要求的行之有效的算法將是今后優(yōu)化算法的主要發(fā)展方向.借助于高效率的尋優(yōu)算法、高性能并行服務(wù)器和高效率的結(jié)構(gòu)重分析程序,形狀優(yōu)化、拓?fù)鋬?yōu)化、多目標(biāo)優(yōu)化、不確定性優(yōu)化和動(dòng)力優(yōu)化也將是未來(lái)發(fā)展的主要方向.本文主要介紹近年來(lái)工程結(jié)構(gòu)優(yōu)化設(shè)計(jì)研究的進(jìn)展情況,并指出工程結(jié)構(gòu)優(yōu)化設(shè)計(jì)的發(fā)展方向.

    1 現(xiàn)代工程結(jié)構(gòu)優(yōu)化設(shè)計(jì)

    1.1 工程結(jié)構(gòu)的不確定性優(yōu)化設(shè)計(jì)

    工程結(jié)構(gòu)設(shè)計(jì)中變量存在不同程度的不確定性,基于不確定性理論的工程結(jié)構(gòu)優(yōu)化設(shè)計(jì)主要考慮變量的隨機(jī)性和模糊性.考慮事物的隨機(jī)性,出現(xiàn)了基于可靠度理論的結(jié)構(gòu)優(yōu)化設(shè)計(jì);研究事物的模糊性,提出了工程結(jié)構(gòu)模糊優(yōu)化設(shè)計(jì).

    1.1.1 基于可靠度理論的工程結(jié)構(gòu)優(yōu)化設(shè)計(jì)

    從可靠度的觀點(diǎn)來(lái)看,傳統(tǒng)的確定性的結(jié)構(gòu)優(yōu)化設(shè)計(jì)不免存在一些缺陷:(a)傳統(tǒng)結(jié)構(gòu)優(yōu)化設(shè)計(jì)不能給出結(jié)構(gòu)可靠性的清晰概念和定量描述,經(jīng)傳統(tǒng)結(jié)構(gòu)優(yōu)化設(shè)計(jì)得出的最優(yōu)結(jié)構(gòu)并不能保證結(jié)構(gòu)有適當(dāng)?shù)目煽啃运?(b)傳統(tǒng)結(jié)構(gòu)優(yōu)化設(shè)計(jì)忽略了構(gòu)件尺寸及材料性能參數(shù)的隨機(jī)性,其結(jié)果不能反映其參數(shù)不確定性這一特點(diǎn).因此,設(shè)計(jì)者應(yīng)把結(jié)構(gòu)的可靠度要求加入到工程結(jié)構(gòu)優(yōu)化設(shè)計(jì)數(shù)學(xué)模型中,按照工程結(jié)構(gòu)功能的要求,根據(jù)所擁有的資源,在結(jié)構(gòu)可靠度分析基礎(chǔ)上進(jìn)行結(jié)構(gòu)優(yōu)化設(shè)計(jì),尋求合理、滿意、安全的設(shè)計(jì)方案.

    目前許多隨機(jī)變量的概率統(tǒng)計(jì)特征還不明了,積極地去收集、統(tǒng)計(jì)、分析、研究隨機(jī)變量的概率特征是基于可靠度理論的工程結(jié)構(gòu)優(yōu)化設(shè)計(jì)的先行條件.

    求解基于可靠度的工程結(jié)構(gòu)優(yōu)化設(shè)計(jì)的直觀解法是把可靠度和優(yōu)化的各自算法組合成一兩層嵌套迭代,該法計(jì)算效率較低.為了改善其收斂性,提高其計(jì)算效率,人們提出單層次算法、功能測(cè)度法、半無(wú)限規(guī)劃法等,蔡迎建等[4]給出了離散變量框架結(jié)構(gòu)可靠性優(yōu)化的相對(duì)差商方法及迭代格式,程耿東等[5]提出了基于可靠度的工程結(jié)構(gòu)優(yōu)化的序列近似規(guī)劃算法,該算法系對(duì)經(jīng)典序列近似規(guī)劃法進(jìn)行改造和擴(kuò)展.

    Dimou等[6]將粒子群優(yōu)化方法引入到基于可靠度的桁架結(jié)構(gòu)優(yōu)化設(shè)計(jì)中,可靠指標(biāo)采用Ditlevsen界限法求解.高亞強(qiáng)等[7]建立基于可靠度,以剛架重量為目標(biāo)函數(shù)的平面剛架優(yōu)化模型,可靠指標(biāo)采用JC法求解,優(yōu)化算法采用約束變尺度法,結(jié)果表明該方法是有效的.王福來(lái)等[8-9]闡述了基于可靠度的橋梁結(jié)構(gòu)優(yōu)化設(shè)計(jì)一般方法及其基本思想,討論了基于可靠度的橋梁結(jié)構(gòu)優(yōu)化設(shè)計(jì)模型的建立,提出了該模型重點(diǎn)研究的內(nèi)容和需要解決的問(wèn)題.

    1.1.2 工程結(jié)構(gòu)模糊優(yōu)化設(shè)計(jì)

    模糊性是指客觀事物差異的中間過(guò)渡中的“不分明性”[1],是由于問(wèn)題過(guò)于復(fù)雜,不能給事物以明確的定義或評(píng)價(jià)標(biāo)準(zhǔn)而形成的不確定性.具有模糊性信息的系統(tǒng)稱為模糊系統(tǒng),也稱軟系統(tǒng).

    與工程結(jié)構(gòu)設(shè)計(jì)有關(guān)的一些因素如:設(shè)計(jì)優(yōu)劣標(biāo)準(zhǔn)、結(jié)構(gòu)響應(yīng)的允許范圍、地震烈度的劃分、建筑場(chǎng)地的分類(lèi)、結(jié)構(gòu)經(jīng)濟(jì)性的衡量、計(jì)算模型和設(shè)計(jì)參數(shù)等都具有模糊性.

    20世紀(jì)70年代,Bellman等[10]給出了模糊決策的定義,并構(gòu)建了模糊環(huán)境下的決策模型[11].隨后,模糊線性規(guī)劃、模糊非線性規(guī)劃、模糊整數(shù)規(guī)劃、模糊多目標(biāo)規(guī)劃及模糊動(dòng)態(tài)規(guī)劃等模型得到積極的研究.王光遠(yuǎn)在文獻(xiàn)[11]中詳細(xì)論述了工程結(jié)構(gòu)的模糊優(yōu)化設(shè)計(jì),并將其應(yīng)用于抗震結(jié)構(gòu)的模糊優(yōu)化中.

    葉萬(wàn)軍等[12]根據(jù)影響黃土路塹邊坡斷面設(shè)計(jì)的經(jīng)濟(jì)性、安全可靠性、抗沖刷性、施工技術(shù)的成熟性及環(huán)保性等5個(gè)指標(biāo),建立了邊坡斷面設(shè)計(jì)各指標(biāo)的模糊判斷矩陣,將定性指標(biāo)轉(zhuǎn)換為模糊值,有效組合主、客觀權(quán)重,建立了基于主客觀權(quán)重綜合模糊識(shí)別的邊坡斷面優(yōu)化決策模型.劉劍等[13]運(yùn)用模糊集合論理論,提出了重力式擋土墻穩(wěn)定性分析的模糊優(yōu)化設(shè)計(jì)方法,從而使擋土墻的穩(wěn)定性分析更加接近客觀實(shí)際.

    1.2 形狀優(yōu)化與拓?fù)鋬?yōu)化

    1.2.1 形狀優(yōu)化

    形狀優(yōu)化是指通過(guò)調(diào)整工程結(jié)構(gòu)內(nèi)外邊界形狀來(lái)改善結(jié)構(gòu)性能和降低工程結(jié)構(gòu)造價(jià),其主要用來(lái)發(fā)掘工程系統(tǒng)構(gòu)件的合理內(nèi)外邊界形狀.形狀優(yōu)化主要包括桁架、剛架類(lèi)的離散變量和塊體、板、殼類(lèi)的連續(xù)變量.

    a.離散變量形狀優(yōu)化.離散變量形狀優(yōu)化一般以節(jié)點(diǎn)坐標(biāo)等幾何變量為設(shè)計(jì)變量,且通常需要同時(shí)考慮尺寸優(yōu)化和形狀優(yōu)化.其主要方法有兩種:一是將兩類(lèi)變量統(tǒng)一處理,采用無(wú)量綱化,構(gòu)造近似問(wèn)題求解.該法的優(yōu)點(diǎn)是可以同時(shí)考慮兩類(lèi)變量的耦合效應(yīng),缺點(diǎn)是計(jì)算工作量大.另一種方法是把尺寸優(yōu)化與形狀優(yōu)化分成兩個(gè)優(yōu)化層次,分別對(duì)兩類(lèi)優(yōu)化交替進(jìn)行,即每步固定一類(lèi)優(yōu)化變量只對(duì)另一類(lèi)優(yōu)化變量進(jìn)行尋優(yōu),兩類(lèi)優(yōu)化層次間通過(guò)迭代協(xié)調(diào).該法的優(yōu)點(diǎn)是求解問(wèn)題規(guī)模較大,缺點(diǎn)是兩類(lèi)優(yōu)化的耦合性較弱,一般僅能求出局部最優(yōu)解.

    b.連續(xù)變量形狀優(yōu)化.連續(xù)型結(jié)構(gòu)的邊界形狀常采用參數(shù)化的曲線或曲面方程等來(lái)描述.其優(yōu)化設(shè)計(jì)方法可歸納為兩類(lèi):數(shù)值法和解析法.數(shù)值法可以應(yīng)用較成熟的約束非線性規(guī)劃法,如GRG法(general reduced gradient)、SQP法(sequential quadratic programming)等[14-15].解析法通過(guò)泛函分析,得到優(yōu)化設(shè)計(jì)函數(shù)的變分形式,進(jìn)而導(dǎo)出形狀函數(shù)對(duì)設(shè)計(jì)變量靈敏度的解析式及最優(yōu)解滿足的條件.亦可先用解析法將結(jié)構(gòu)優(yōu)化設(shè)計(jì)函數(shù)簡(jiǎn)單化,再用數(shù)值法尋優(yōu).

    與尺寸優(yōu)化相比,形狀優(yōu)化發(fā)展緩慢.其主要原因是:(a)自編結(jié)構(gòu)優(yōu)化程序網(wǎng)格自適應(yīng)劃分存在一定難度;(b)形狀優(yōu)化的敏度分析計(jì)算量大,程序收斂性差.綜合解決上述問(wèn)題需要數(shù)學(xué)、力學(xué)、計(jì)算機(jī)以及工程人員的共同努力.

    劉齊茂等[16]提出一種求解桁架形狀優(yōu)化設(shè)計(jì)問(wèn)題的改進(jìn)模擬退火算法,算例表明該方法能獲得較高質(zhì)量的解.易少華等[17]利用遺傳算法(GA)對(duì)桁架結(jié)構(gòu)分別進(jìn)行桿件截面和節(jié)點(diǎn)坐標(biāo)的優(yōu)化設(shè)計(jì),并通過(guò)25桿輸電塔架算例的優(yōu)化設(shè)計(jì),驗(yàn)證了此方法的實(shí)際應(yīng)用是可行的.曹天捷[18]針對(duì)集中載荷、均布載荷以及填料重量載荷同時(shí)作用下的三鉸拱結(jié)構(gòu),以高次拋物線為拱軸線方程,提出了合理拱軸線的方法.王棟等[19]對(duì)空間桁架結(jié)構(gòu)的形狀進(jìn)行優(yōu)化設(shè)計(jì),使結(jié)構(gòu)重量達(dá)到最小.

    1.2.2 拓?fù)鋬?yōu)化

    拓?fù)鋬?yōu)化是指尋求結(jié)構(gòu)剛度在設(shè)計(jì)空間最佳的分布形式或結(jié)構(gòu)最佳的傳力形式,其主要用于求解工程結(jié)構(gòu)受力、位移等的最合理形式,有助于發(fā)掘目前還未知的最優(yōu)結(jié)構(gòu)拓?fù)?

    與尺寸優(yōu)化和形狀優(yōu)化相比,拓?fù)鋬?yōu)化能夠在工程結(jié)構(gòu)設(shè)計(jì)的初始階段搜索出結(jié)構(gòu)的最優(yōu)布局方案,進(jìn)而產(chǎn)生可觀的經(jīng)濟(jì)效益,且容易被工程設(shè)計(jì)人員所接受.拓?fù)鋬?yōu)化變量亦包括離散型和連續(xù)型.

    a.離散變量拓?fù)鋬?yōu)化.1964年,Dorn等[20]以結(jié)構(gòu)節(jié)點(diǎn)、支座點(diǎn)及荷載作用點(diǎn)為節(jié)點(diǎn)集合,集合中所有節(jié)點(diǎn)之間采用桿件單元連接的基結(jié)構(gòu),并以內(nèi)力為設(shè)計(jì)變量,以應(yīng)力為約束函數(shù),建立了單工況線性規(guī)劃優(yōu)化設(shè)計(jì)模型.該法計(jì)算效率較高,但不能應(yīng)用于多工況和有位移約束的優(yōu)化設(shè)計(jì)問(wèn)題上.Dobbs等[21]以截面面積為設(shè)計(jì)變量,采用最速下降法(steepest descent method,SDM)成功地解決了多工況應(yīng)力約束下桁架結(jié)構(gòu)的拓?fù)鋬?yōu)化問(wèn)題.Kirsch等[22]提出了兩階段算法,第一階段以桿件截面積和贅余內(nèi)力為設(shè)計(jì)變量,不考慮位移約束和變形協(xié)調(diào)條件,將離散變量拓?fù)鋬?yōu)化轉(zhuǎn)化為線性規(guī)劃優(yōu)化設(shè)計(jì);第二階段考慮所有約束,在已有的拓?fù)浣Y(jié)構(gòu)上,將離散變量拓?fù)鋬?yōu)化轉(zhuǎn)化為非線性規(guī)劃優(yōu)化設(shè)計(jì).Lipson等[23]建議在多工況下以桿件內(nèi)力為準(zhǔn)則來(lái)判斷應(yīng)刪除的桿件.

    b.連續(xù)變量拓?fù)鋬?yōu)化.連續(xù)變量拓?fù)鋬?yōu)化設(shè)計(jì)是一種0-1離散變量的組合優(yōu)化問(wèn)題.其基本思想是將設(shè)計(jì)區(qū)域離散為有限網(wǎng)格,根據(jù)相應(yīng)的準(zhǔn)則,刪除某些網(wǎng)格.其主要方法有:均勻化法、變密度法和變厚法.

    均勻化法以微結(jié)構(gòu)的單胞尺寸為設(shè)計(jì)變量,以單胞尺寸的增減實(shí)現(xiàn)微結(jié)構(gòu)的增刪和復(fù)合.其特點(diǎn)是:數(shù)學(xué)理論推導(dǎo)嚴(yán)密,可獲得宏觀的彈性常數(shù)和局部應(yīng)力應(yīng)變,容易收斂到局部最優(yōu)解,計(jì)算量大,求解的問(wèn)題類(lèi)型有限,容易引起棋盤(pán)效應(yīng).

    變密度法假設(shè)優(yōu)化設(shè)計(jì)對(duì)象的材料密度是可變的,且假定材料物理參數(shù)與密度間存在某種數(shù)學(xué)關(guān)系,以材料密度為優(yōu)化設(shè)計(jì)變量,以材料的最優(yōu)分布為目標(biāo)函數(shù).通常單元密度與彈性模量之間的關(guān)系采用人為給出的冪函數(shù)規(guī)律.目前主要采用SIMP(solid isotropic microstructure with penalization)法.其特點(diǎn)是:與采用尺寸變量相比,更能反映拓?fù)鋬?yōu)化的本質(zhì)特征,程序?qū)崿F(xiàn)簡(jiǎn)單,計(jì)算效率高,但解的精度較差.

    變厚度法以基結(jié)構(gòu)中單元厚度為優(yōu)化設(shè)計(jì)變量,以結(jié)構(gòu)中厚度的最優(yōu)分布為目標(biāo)函數(shù).其特點(diǎn)是:方法簡(jiǎn)單,應(yīng)用范圍狹窄.

    姜冬菊等[24-25]將相對(duì)差商法和混沌優(yōu)化結(jié)合起來(lái),導(dǎo)出求解離散變量桁架結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)的混合算法,該算法具有差商法的離散變量快速尋優(yōu)以及混沌法的全局遍歷性特點(diǎn).穆春燕等[26]利用均勻化拓?fù)鋬?yōu)化方法,根據(jù)拱壩特點(diǎn),以體積最小為目標(biāo),對(duì)拱壩的體形優(yōu)化進(jìn)行了探討分析,得到了較為滿意的結(jié)果.張莉等[27]以各向異性連續(xù)體為基結(jié)構(gòu),采用類(lèi)桁架連續(xù)體材料模型進(jìn)行結(jié)構(gòu)拓?fù)鋬?yōu)化,以材料在節(jié)點(diǎn)位置的密度和方向作為優(yōu)化設(shè)計(jì)變量,使材料在設(shè)計(jì)域內(nèi)得到最優(yōu)分布.

    1.3 多目標(biāo)優(yōu)化

    多目標(biāo)優(yōu)化是指優(yōu)化過(guò)程中所考慮的優(yōu)化目標(biāo)不是單一的.一般情況下各目標(biāo)函數(shù)之間往往相互矛盾,比如要取得好的安全性,就要求結(jié)構(gòu)的截面面積要大,而為了取得最少重量,又要求截面面積較小.因此不存在使所有目標(biāo)都達(dá)到最優(yōu)的“絕對(duì)最優(yōu)解”,只能求得“滿意解集”,由決策者最終選定某一個(gè)滿意解作為最后定解.

    實(shí)際工程中,多目標(biāo)優(yōu)化一般用于工程系統(tǒng)決策,即在工程決策方面先采用多目標(biāo)優(yōu)化進(jìn)行方案確定,再優(yōu)化各個(gè)分目標(biāo).

    不同的優(yōu)化設(shè)計(jì)數(shù)學(xué)模型有不同的求解方法.主要有以下幾種方法[14-15]:(a)約束法.在多個(gè)分目標(biāo)中選擇一個(gè)為主目標(biāo),對(duì)其余分目標(biāo)給出希望值,進(jìn)而轉(zhuǎn)化為單目標(biāo)優(yōu)化問(wèn)題求解.(b)功效系數(shù)法.將各分目標(biāo)的“壞”價(jià)值用統(tǒng)一的功效系數(shù)表達(dá),而后采用幾何平均構(gòu)成評(píng)價(jià)函數(shù),進(jìn)而轉(zhuǎn)化成單目標(biāo)優(yōu)化問(wèn)題求解. (c)評(píng)價(jià)函數(shù)法.采用線性加權(quán)、平方和加權(quán)等方法將分目標(biāo)函數(shù)綜合成一個(gè)總函數(shù)進(jìn)而轉(zhuǎn)化為單目標(biāo)優(yōu)化問(wèn)題求解.(d)目的規(guī)劃法.希望值與真實(shí)值之間的差值稱為約束偏差,以約束偏差和目標(biāo)偏差的某種組合作為總函數(shù)進(jìn)而轉(zhuǎn)化為單目標(biāo)優(yōu)化問(wèn)題求解.(e)多屬性效用函數(shù)法.實(shí)際多目標(biāo)優(yōu)化時(shí)往往得到的不是某一個(gè)最優(yōu)解,而是最優(yōu)解的一個(gè)集合,再在這個(gè)集合中選出需要的最優(yōu)解.為此可應(yīng)用效用理論建立決策者的效用函數(shù)(曲線),按此曲線從有限解集中選出最終的合適方案.

    謝能剛等[28-29]提出了基于博弈理論、無(wú)私合作博弈理論和模糊理論,考慮壩體體積、最大主拉應(yīng)力、整體應(yīng)變能、高應(yīng)力區(qū)范圍為目標(biāo)函數(shù)的拱壩體型多目標(biāo)優(yōu)化設(shè)計(jì).曹洪濤等[30]以Matlab多目標(biāo)優(yōu)化函數(shù)為工具,研究了不同截面鋼與混凝土組合梁有效發(fā)揮抗彎承載力和降低截面尺寸的多目標(biāo)優(yōu)化問(wèn)題.黃煒等[31]研究了密肋復(fù)合墻體多目標(biāo)優(yōu)化設(shè)計(jì)方法,提出了綜合考慮墻體的安全性(墻體斜截面抗剪承載力)與經(jīng)濟(jì)性(墻體造價(jià))的多目標(biāo)優(yōu)化設(shè)計(jì)數(shù)學(xué)模型,并用序列二次規(guī)劃法(sequential quadratic programming,SQP)求解.

    1.4 工程系統(tǒng)優(yōu)化

    實(shí)際工程項(xiàng)目大多是由若干個(gè)建筑物或構(gòu)筑物組成,它們構(gòu)成一個(gè)結(jié)構(gòu)系列.把結(jié)構(gòu)系列中那些主要的、投資較大的結(jié)構(gòu)所構(gòu)成的集合稱為“工程結(jié)構(gòu)系統(tǒng)”或簡(jiǎn)稱“工程系統(tǒng)”[32].對(duì)“工程系統(tǒng)”進(jìn)行全局優(yōu)化稱為系統(tǒng)優(yōu)化.

    工程系統(tǒng)可以分為串聯(lián)、并聯(lián)、復(fù)合的并串聯(lián)系統(tǒng),土木、建筑、水利工程項(xiàng)目主要簡(jiǎn)化為串聯(lián)工程系統(tǒng).針對(duì)不同的工程系統(tǒng)有不同的全局優(yōu)化方法.王光遠(yuǎn)在文獻(xiàn)[32]中,以工程系統(tǒng)的總造價(jià)和總損失期望為目標(biāo)函數(shù),詳細(xì)闡述了簡(jiǎn)單工程系統(tǒng)全局優(yōu)化的直接解法、串聯(lián)和并聯(lián)工程系統(tǒng)的的全局優(yōu)化、簡(jiǎn)單復(fù)合工程系統(tǒng)的全局優(yōu)化、表決系統(tǒng)的全局優(yōu)化以及網(wǎng)絡(luò)工程系統(tǒng)的全局優(yōu)化.仿照文獻(xiàn)[32]中相關(guān)理論和方法將基于最優(yōu)設(shè)防水平下的工程系統(tǒng)優(yōu)化推廣到基于最大可靠指標(biāo)下的工程系統(tǒng)優(yōu)化,將是未來(lái)發(fā)展的一個(gè)主要方向.

    朱靖華等在文獻(xiàn)[33-36]中闡述了多級(jí)遞階工程系統(tǒng)的整體優(yōu)化和全局優(yōu)化,提出了全局優(yōu)化計(jì)算的混合法,該方法具有耦合法與分解法的優(yōu)點(diǎn),能靈活地處理系統(tǒng)結(jié)構(gòu)邊界條件,使優(yōu)化空間達(dá)到高度降維,對(duì)于大規(guī)模、分布復(fù)雜的多級(jí)遞階系統(tǒng)具有很高的優(yōu)化效率.

    2 尋優(yōu)算法的新思路

    優(yōu)化設(shè)計(jì)方法大致可以分為三類(lèi):最優(yōu)準(zhǔn)則法、數(shù)學(xué)規(guī)劃法和現(xiàn)代優(yōu)化算法[14-15].最優(yōu)準(zhǔn)則法考慮的問(wèn)題比較局限,對(duì)不同性質(zhì)的約束要采用不同的準(zhǔn)則,所得的解不一定是最優(yōu)解;數(shù)學(xué)規(guī)劃法計(jì)算量大,收斂慢,因而新的現(xiàn)代尋優(yōu)算法應(yīng)運(yùn)而生.隨著人們對(duì)自然界認(rèn)識(shí)的加深,將自然界的過(guò)程和結(jié)構(gòu)抽象出來(lái),形成許多仿生學(xué)算法,如:遺傳算法、模擬退火算法、神經(jīng)網(wǎng)絡(luò)算法.

    2.1 遺傳算法(GA)

    遺傳算法是模擬Darwin和Mendel的遺傳選擇和自然淘汰的生物進(jìn)化過(guò)程的計(jì)算方法,是一種通過(guò)模擬自然進(jìn)化過(guò)程搜索最優(yōu)解的方法.

    1975年,Holland教授用二進(jìn)制位串來(lái)表示染色體,根據(jù)遺傳法則,由父代二進(jìn)制位串經(jīng)過(guò)再生、變異和交叉等操作后產(chǎn)生下代二進(jìn)制位串來(lái)模擬生物群體的進(jìn)化歷程[14-15],并出版了頗有影響的專(zhuān)著Adaptation in Natural and Artificial Systems[37].1998年,Michalewicz出版了 Genetic Algorithms+Data Structures=Evolution Programs[38],推動(dòng)了遺傳算法應(yīng)用于最優(yōu)化問(wèn)題.

    遺傳算法的主要優(yōu)點(diǎn)是:解題能力較強(qiáng),適應(yīng)性較廣,平行性良好,全局優(yōu)化性能較強(qiáng)[14-15].其缺點(diǎn)是:算法操作隨機(jī)性較大,計(jì)算工作量大.通過(guò)引入模式理論、積木塊假設(shè),良好的編碼方式、并行遺傳算法及算法控制參數(shù)(群體大小、變異概率、交叉概率)的智能調(diào)整等,遺傳算法的尋優(yōu)效率得到較大提高.

    工程結(jié)構(gòu)中,遺傳算法主要用于桁架結(jié)構(gòu)、框架結(jié)構(gòu)、網(wǎng)格結(jié)構(gòu)以及振動(dòng)控制結(jié)構(gòu)的優(yōu)化.柳春光等[39]將遺傳算法引入橋梁的震害預(yù)測(cè)中,并與傳統(tǒng)的單獨(dú)BP神經(jīng)網(wǎng)絡(luò)相比較,結(jié)果表明該方法能夠有效、準(zhǔn)確地對(duì)橋梁結(jié)構(gòu)進(jìn)行震害預(yù)測(cè).汪大洋等[40]以結(jié)構(gòu)優(yōu)化中的遺傳算法理論為基礎(chǔ),建立6種不同的遺傳算法目標(biāo)函數(shù),研究表明遺傳算法優(yōu)化方法能依據(jù)結(jié)構(gòu)的具體動(dòng)力響應(yīng)在結(jié)構(gòu)層間合理分配阻尼系數(shù).陳文英等[41]采用免疫遺傳算法解決同時(shí)含有離散和連續(xù)變量的優(yōu)化設(shè)計(jì)問(wèn)題.楊海霞等[42]提出了一種可以任意控制離散度的改進(jìn)實(shí)數(shù)編碼遺傳算法,該算法的優(yōu)點(diǎn)是可以任意選擇變量的精度,適合復(fù)雜結(jié)構(gòu)的優(yōu)化設(shè)計(jì).趙同彬等[43]等將遺傳算法引入到擋土墻庫(kù)侖土壓力求解分析中,計(jì)算結(jié)果表明遺傳算法在計(jì)算擋土墻庫(kù)侖主動(dòng)土壓力的過(guò)程中,收斂速度快,具有較高的計(jì)算精度.

    2.2 模擬退火算法(SA)

    退火是將固體升溫到一定溫度,再讓其慢慢冷卻,升溫時(shí),固體內(nèi)部粒子隨溫升變?yōu)闊o(wú)序狀,內(nèi)能增大,而慢慢冷卻時(shí)固體內(nèi)部粒子漸漸趨于有序狀.固體內(nèi)部粒子在每個(gè)溫度都達(dá)到平衡態(tài),最后在常溫時(shí)達(dá)到基態(tài),內(nèi)能減為最小.

    1953年,Metropolis等提出了模擬退火算法.1983年,Kirkpatriek等將此法加以改進(jìn)后用于大規(guī)模集成電路優(yōu)化問(wèn)題.1987年,Laarhovent等出版了Simulated Annealing:Theory and Applications[44]一書(shū),對(duì)模擬退火算法作了比較系統(tǒng)的總結(jié),促進(jìn)了模擬退火算法的理論研究和實(shí)際應(yīng)用的發(fā)展.

    模擬退火算法的主要優(yōu)點(diǎn):適用于離散型、連續(xù)型及混合型變量;魯棒性、全局收斂性、隱含并行性較強(qiáng),適應(yīng)性較廣.其主要缺點(diǎn):效率不高;參數(shù)、試探數(shù)的控制上限選取準(zhǔn)則性較弱.通過(guò)與遺傳算法、混沌理論、粒子群算法和直接搜索法進(jìn)行組合,可以大大提高模擬退火算法的尋優(yōu)效率.

    陳昌富等[45]將復(fù)合形模擬退火算法引入到水泥土墻支護(hù)結(jié)構(gòu)參數(shù)優(yōu)化設(shè)計(jì)中,結(jié)果表明了該方法的可行性和可靠性.王永德等[46]首次將模擬退火算法應(yīng)用于電力網(wǎng)線路優(yōu)化設(shè)計(jì)中,同時(shí)提出網(wǎng)架優(yōu)化的模擬退火算法.高占遠(yuǎn)等[47]將模擬退火算法引入到網(wǎng)架結(jié)構(gòu)優(yōu)化設(shè)計(jì)中,并對(duì)模擬退火算法中參數(shù)的選擇進(jìn)行了探討.項(xiàng)寶衛(wèi)等[48]針對(duì)桁架形狀優(yōu)化問(wèn)題,提出了一種改進(jìn)的模擬退火算法,計(jì)算結(jié)果表明該法計(jì)算精度高、通用性強(qiáng).顧元憲等[49]將模擬退火算法應(yīng)用于桁架結(jié)構(gòu)截面尺寸優(yōu)化設(shè)計(jì),結(jié)果表明,改進(jìn)模擬退火算法具有很高的優(yōu)化求解精度,計(jì)算效率顯著提高且優(yōu)于遺傳算法.

    2.3 人工神經(jīng)網(wǎng)絡(luò)算法(ANN)

    人工神經(jīng)網(wǎng)絡(luò)是由大量簡(jiǎn)單的神經(jīng)元按某種方式連接形成的智能仿生動(dòng)態(tài)網(wǎng)絡(luò),它以簡(jiǎn)單非線性神經(jīng)元作為處理單元,通過(guò)廣泛連接構(gòu)成大規(guī)模分布式并行處理非線性動(dòng)力學(xué)系統(tǒng).主要由輸入層、隱含層和輸出層組成.

    1943年,法國(guó)心理學(xué)家W.S.McCuloch和W.Pitts在分析、綜合神經(jīng)元基本特征的基礎(chǔ)上提出了第一個(gè)神經(jīng)元數(shù)學(xué)模型(M-P模型),開(kāi)創(chuàng)了人工神經(jīng)網(wǎng)絡(luò)算法的研究.1982年Hopfield[50]提出了人工神經(jīng)網(wǎng)絡(luò)模型.Hopfield開(kāi)始主要針對(duì)線性規(guī)劃和0-1離散規(guī)劃問(wèn)題,隨后,將其推廣到非線性優(yōu)化問(wèn)題上.1985年Rumelhart等提出了多層網(wǎng)絡(luò)的反向傳播學(xué)系算法(back propagation,BP),該算法能夠較大幅度提高神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)與知識(shí)表達(dá)能力.

    人工神經(jīng)網(wǎng)絡(luò)算法主要優(yōu)點(diǎn):運(yùn)算能力、自學(xué)習(xí)、自適應(yīng)能力及非線性映射能力較強(qiáng),容錯(cuò)性較好.其主要缺點(diǎn):容易陷入局部最優(yōu)解,計(jì)算工作量大.通過(guò)提高網(wǎng)絡(luò)的泛化性能、與遺傳算法和粒子群算法相結(jié)合,可以大大提高人工神經(jīng)網(wǎng)絡(luò)算法的尋優(yōu)效率.

    姚宏等[51]利用進(jìn)化人工神經(jīng)網(wǎng)絡(luò)算法預(yù)測(cè)結(jié)構(gòu)抗震可靠度,克服了BP網(wǎng)絡(luò)輸入層節(jié)點(diǎn)數(shù)、隱含層數(shù)和隱含層節(jié)點(diǎn)數(shù)確定的盲目性,其適應(yīng)性更強(qiáng),精度更高.任遠(yuǎn)等[52]以提高網(wǎng)絡(luò)的泛化性能為主旨,提出了GA-BP貝葉斯算法,并應(yīng)用于方形扁平封裝器中,該算法克服了現(xiàn)有網(wǎng)絡(luò)訓(xùn)練算法對(duì)初始權(quán)值敏感、建模精度不高的缺點(diǎn).

    2.4 其他新尋優(yōu)算法

    新的現(xiàn)代尋優(yōu)算法還有1989年F.Glover提出的禁忌搜索算法(tabu search)[53];Toyoo Fukuda等提出的一種基于信息熵概念的免疫算法,該法全局收斂性良好;蟻群算法(ant colony optimization,ACO)是通過(guò)模擬自然界中真實(shí)蟻群的覓食行為而產(chǎn)生的;粒子群優(yōu)化算法(partical swarm optimization,PSO)是近年來(lái)發(fā)展起來(lái)的一種新的進(jìn)化算法,源于對(duì)鳥(niǎo)群捕食的行為研究.粒子群優(yōu)化算法和遺傳算法相似,它均從隨機(jī)解出發(fā),通過(guò)迭代尋找最優(yōu)解,采用適應(yīng)度來(lái)評(píng)價(jià)解的品質(zhì),但是比遺傳算法規(guī)則更為簡(jiǎn)單,沒(méi)有遺傳算法的“交叉”和“變異”操作.

    3 結(jié) 語(yǔ)

    a.工程結(jié)構(gòu)優(yōu)化設(shè)計(jì)從尺寸優(yōu)化到形狀優(yōu)化、拓?fù)鋬?yōu)化,從單目標(biāo)優(yōu)化到多目標(biāo)優(yōu)化,從確定性優(yōu)化到不確定性優(yōu)化,從靜力優(yōu)化到靜動(dòng)力優(yōu)化,其優(yōu)化設(shè)計(jì)的層次向更高方向發(fā)展.優(yōu)化算法從傳統(tǒng)的準(zhǔn)則法、數(shù)學(xué)規(guī)劃法,發(fā)展到仿生學(xué)算法.優(yōu)化的理論和方法都得到長(zhǎng)足發(fā)展.

    b.尋求目標(biāo)函數(shù)和約束函數(shù)的高精度近似顯式解析式仍然是結(jié)構(gòu)重分析算法發(fā)展的一個(gè)重要方向.

    c.最優(yōu)準(zhǔn)則法、數(shù)學(xué)規(guī)劃法及現(xiàn)代的仿生學(xué)算法都有其自身局限性.針對(duì)實(shí)際工程,對(duì)上述各尋優(yōu)算法進(jìn)行改進(jìn)、組合,產(chǎn)生新的優(yōu)化算法,是工程結(jié)構(gòu)優(yōu)化算法發(fā)展的主要方向.

    d.工程結(jié)構(gòu)優(yōu)化設(shè)計(jì)中尺寸優(yōu)化已相對(duì)完善,形狀優(yōu)化逐漸走向成熟,拓?fù)鋬?yōu)化及布局優(yōu)化呈現(xiàn)加速發(fā)展的態(tài)勢(shì).

    [1]蔡新,郭興文,張旭明.工程結(jié)構(gòu)優(yōu)化設(shè)計(jì)[M].北京:中國(guó)水利水電出版社,2003.

    [2]張炳華,侯昶.土建結(jié)構(gòu)優(yōu)化設(shè)計(jì)[M].2版.上海:同濟(jì)大學(xué)出版社,1998.

    [3]SCHMIT L A.Structural design by systematic synthesis[C]//Proceeding of the second Conference on Electronic Computation.New York: ASCE,1960:105-132.

    [4]蔡迎建,孫煥純.框架結(jié)構(gòu)的可靠性優(yōu)化[J].計(jì)算力學(xué)學(xué)報(bào),1999,16(3):288-294.(CAI Ying-jian,SUN Huan-chu.Reliability based on optimization for frame structur[J].Chinese Journal of Computational Mechanics,1999,16(3):288-294.(in Chinese))

    [5]程耿東,許林.基于可靠度的結(jié)構(gòu)優(yōu)化的序列近似規(guī)劃算法[J].計(jì)算力學(xué)學(xué)報(bào),2006,23(6):641-646.(CHENG Geng-dong,XU Lin.Sequential approximate programming approach to reliability based structural optimization[J].Chinese Journal of Computational Mechanics,2006,23(6):641-646.(in Chinese))

    [6]DIMOU C K,KOUMOUSIS V K.Reliability-based optimal design of truss structures using particle swarm optimization[J].Journal of Computing in Civil Engineering,2009,23(2):100-109.

    [7]高亞強(qiáng),王興國(guó).基于可靠性的結(jié)構(gòu)優(yōu)化設(shè)計(jì)[J].河北理工大學(xué)學(xué)報(bào):自然科學(xué)版,2008,30(1):120-123.(GAO Ya-qiang,WANG Xing-guo.A Study on reliability-based structural optimization and design[J].Journal of Hebei Polytechnic University:Natural Science Edition,2008,30(1):120-123.(in Chinese))

    [8]王福來(lái),葉少有.基于可靠度理論的橋梁結(jié)構(gòu)優(yōu)化設(shè)計(jì)[J].山西建筑,2008,34(2):89-90.(WANG Fu-lai,YE Shao-you. Optimal design of bridge structures based on reliability[J].Shanxi Architecture,2008,34(2):89-90.(in Chinese))

    [9]張鵬飛.淺談可靠度在橋梁結(jié)構(gòu)優(yōu)化設(shè)計(jì)中的應(yīng)用[J].山西建筑,2007,33(20):78-79.(ZHANG Peng-fei.Application of the reliability in bridge structural optimum design[J].Shanxi Architecture,2007,33(20):78-79.(in Chinese))

    [10]BELLMAN R E,ZADEH L A.Decision-making in a fuzzy environment[R].Berkeley:University of California Berkeley,1970.

    [11]王光遠(yuǎn).工程軟設(shè)計(jì)理論[M].北京:科學(xué)出版社,1992.

    [12]葉萬(wàn)軍,彭小云,折學(xué)森.綜合模糊識(shí)別下的黃土邊坡斷面優(yōu)化[J].西安科技大學(xué)學(xué)報(bào),2007,27(4):584-589.(YE Wan-jun,PENG Xiao-yun,SHE Xue-sen.Loess slope section optimization based on comprehensive fuzzy recognition[J].Journal of Xi'an University of Science and Technoidgy,2007,27(4):584-589.(in Chinese))

    [13]劉劍,汪志勇,肖盛燮.擋土墻穩(wěn)定性分析中的模糊優(yōu)化設(shè)計(jì)[J].中南公路工程,2006,31(3):102-104.(LIU Jian,WANG Zhiyong,XIAO Sheng-xie.Fuzzy optimization designin the stability of retaining wall analysed[J].Central SouthHighway Engineering,2006,31(3):102-104.(in Chinese))

    [14]汪樹(shù)玉,劉國(guó)華,包志仁.結(jié)構(gòu)優(yōu)化設(shè)計(jì)的現(xiàn)狀與進(jìn)展(續(xù))[J].基建優(yōu)化,1999,20(5):3-7.(WANG Shu-yu,LIU Guo-hua,BAO Zhi-ren.Present situation and development of structuraloptimizationmethod(continuation)[J].Optimizqtionof Capital Construction,1999,20(5):3-7.(in Chinese))

    [15]汪樹(shù)玉,劉國(guó)華,包志仁.結(jié)構(gòu)優(yōu)化設(shè)計(jì)的現(xiàn)狀與進(jìn)展[J].基建優(yōu)化,1999,20(4):3-14.(WANG Shu-yu,LIU Guo-hua,BAO Zhi-ren.Present situation and development of structural optimization method[J].Optimizqtion of Capital Construction,1999,20(4):3-14. (in Chinese))

    [16]劉齊茂,燕柳斌,鄧?yán)誓?桁架形狀優(yōu)化的一種改進(jìn)模擬退火算法研究[J].計(jì)算機(jī)工程與應(yīng)用,2007,43(23):218-221.(LIU Qi-mao,YAN Liu-bin,DENG Lang-ni.Improved simulated annealing algorithm study on truss shape optimal design[J].Computer Engineering and Applications,2007,43(23):218-221.(in Chinese))

    [17]易少華,劉樹(shù)堂.BP遺傳算法在結(jié)構(gòu)形狀優(yōu)化設(shè)計(jì)中的應(yīng)用[J].廣州大學(xué)學(xué)報(bào):自然科學(xué)版,2007,6(3):82-86.(YI Shaohua,LIU Shu-tang.Application of BP genetic algorithms on shape optimization design of truss structure[J].Journal of Guangzhou University:Natural Science Edition,2007,6(3):82-86.(in Chinese))

    [18]曹天捷.幾種載荷同時(shí)作用下三鉸拱軸線的形狀優(yōu)化[J].中國(guó)民航大學(xué)學(xué)報(bào),2007,25(2):33-35.(CAO Tian-jie.Optimization of axial curve of three-hinged arch under action of few kind of loads[J].Journal of Civil Aviation University of China,2007,25(2):33-35.(in Chinese))

    [19]王棟,李晶.空間桁架結(jié)構(gòu)動(dòng)力學(xué)形狀優(yōu)化設(shè)計(jì)[J].工程力學(xué),2007,24(4):129-134.(WANG Dong,LI Jing.Shape optimization of space trusses subject to frequency constrains[J].Engineering Mechanics,2007,24(4):129-134.(in Chinese))

    [20]DOR N W S,GEMORY R E,GREENBERG H J.Automatic design of optimal structures[J].J Mech,1964,3(1):25-52.

    [21]DOBBS M W,F(xiàn)ELTON L P.Optimization of truss geometry[J].J Struct Div Proc ASCE,1969,95(ST10):2105-2118.

    [22]KIRSCH U,TOPPING B H V.Minimum weight design of structural topologies[J].J Struc Eng,1992,118(7):1770-1785.

    [23]LIPSON S L,GWIN L B.The complex method applied to optimal truss configuration[J].Computers&Structures,1977,7(3):461-468.

    [24]姜冬菊,王德信.離散變量桁架結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)的混合算法[J].工程力學(xué),2007,24(1):112-116.(JIANG Dong-ju,WANG De-xin.A hybrid algorithm for topology optimization of truss structures with discrete variables[J].Engineering Mechanics,2007,24(1): 112-116.(in Chinese))

    [25]姜冬菊.結(jié)構(gòu)拓?fù)浜筒季謨?yōu)化及工程應(yīng)用研究[D].南京:河海大學(xué),2008.

    [26]穆春燕,蘇超,叢賽飛.拱壩二維水平拱圈的拓?fù)鋬?yōu)化研究[J].黑龍江水專(zhuān)學(xué)報(bào),2007,34(4):1-3.(MU Chun-yan,SU Chao,CONG Sai-fei.Study on the topology optimization of 2-d horizontal arch ring[J].Journal of Heilongjiang Hydraulic Engineering College,2007,34(4):1-3.(in Chinese))

    [27]張莉,周克民.基于類(lèi)桁架連續(xù)體的結(jié)構(gòu)拓?fù)鋬?yōu)化方法與應(yīng)用[J].應(yīng)用力學(xué)學(xué)報(bào),2007,24(3):412-415.(ZHANG Li,ZHOU Ke-min.Topological optimization based on truss-like continuum with applications[J].Chinese Journal of Applied Mechanics,2007,24 (3):412-415.(in Chinese))

    [28]謝能剛,孫林松,包家漢,等.基于無(wú)私合作博弈模型的拱壩體型多目標(biāo)優(yōu)化設(shè)計(jì)[J].水利學(xué)報(bào),2007,38(9):1109-1114. (XIENeng-gang,SUN Lin-song,BAO Jia-han,et al.Multi-objective optimization design of arch dam shape based on unselfish cooperation game model[J].Journal of Hydraulic Engineering,2007,38(9):1109-1114.(in Chinese))

    [29]孫林松,張偉華,謝能剛,等.基于博弈理論的拱壩體形多目標(biāo)優(yōu)化設(shè)計(jì)[J].河海大學(xué)學(xué)報(bào):自然科學(xué)版,2006,34(4):392-396.(SUN Lin-song,ZHANG Wei-hua,XIE Neng-gang,et al.Multi-objective optimization for shape design of arch dams based on game theory[J].Journal of Hohai University:Natural Sciences,2006,34(4):392-396.(in Chinese))

    [30]曹洪濤,董軍,居雅楠.基于MATLAB的鋼-混凝土組合梁截面多目標(biāo)優(yōu)化設(shè)計(jì)[J].四川建筑科學(xué)研究,2008,34(1):5-8. (CAO Hong-tao,DONG Jun,JU Ya-nan.Multi-objective optimal section design of steel-concrete composite beams basedon theMATLAB [J].Sichuan Building Science,2008,34(1):5-8.(in Chinese))

    [31]黃煒,侯莉娜,姚謙峰,等.密肋復(fù)合墻體多目標(biāo)優(yōu)化設(shè)計(jì)方法研究[J].西安建筑科技大學(xué)學(xué)報(bào):自然科學(xué)版,2007,39(5): 641-645.(HUANG Wei,HOU Li-na,YAO Qian-feng,et al.Research on multi-objective optimization design method of multi-ribbed composite wall[J].Journal of Xi'an University of Architecture&Technology:Natural Science Edition,2007,39(5):641-645.(in Chinese))

    [32]王光遠(yuǎn).工程結(jié)構(gòu)與系統(tǒng)抗震優(yōu)化設(shè)計(jì)的實(shí)用方法[M].北京:中國(guó)建筑工業(yè)出版社,1999.

    [33]朱靖華,王光遠(yuǎn).多級(jí)遞階工程結(jié)構(gòu)系統(tǒng)抗地震設(shè)計(jì)的全局優(yōu)化[J].工程力學(xué),2002,19(6):171-180.(ZHU Jing-hua,WANG Guang-yuan.Global optimization of anti-seismic design for multilevel hierarchical engineering structure systems[J].Engineering Mechanics,2002,19(6):171-180.(in Chinese))

    [34]朱靖華,王光遠(yuǎn).多級(jí)遞階工程系統(tǒng)全局優(yōu)化抗震設(shè)計(jì)的混合法[J].土木工程學(xué)報(bào),2005,38(4):18-24,52.(ZHU Jing-hua,WANG Guang-yuan.Global optimization of anti-seismic design for multilevel hierarchical engineering systems by using a mixed method [J].China Civil Engineering Journal,2005,38(4):18-24,52.(in Chinese))

    [35]朱靖華,王光遠(yuǎn).多級(jí)遞階工程系統(tǒng)抗地震設(shè)計(jì)全局優(yōu)化的對(duì)偶分解法[J].建筑結(jié)構(gòu)學(xué)報(bào),2002,23(6):33-40.(ZHU Jinghua,WANG Guang-yuan.Global optimization of anti-seismic design for multilevel hierarchical engineering systems by a dualdecomposition method[J].Journal of Building Structures,2002,23(6):33-40.(in Chinese))

    [36]朱靖華,王光遠(yuǎn).多級(jí)遞階工程系統(tǒng)抗地震設(shè)計(jì)整體優(yōu)化的計(jì)算實(shí)現(xiàn)[J].計(jì)算力學(xué)學(xué)報(bào),2002,19(4):384-392.(ZHU Jinghua,WANG Guang-yuan.Computational implementation of global optimization of anti-seismic design for multilevel hierarchical engineering systems[J].Chinese Journal of Computational Mechanics,2002,19(4):384-392.(in Chinese))

    [37]HOLLAND J H.Adaptation in natural and artificial systems[M].Cambridge:MIT Press,1992.

    [38]MICHALEWICZ Z.Genetic algorithms+data structures=evolution programs[M].Berlin:Springer,1998.

    [39]柳春光,張利華.遺傳優(yōu)化神經(jīng)網(wǎng)絡(luò)方法在橋梁震害預(yù)測(cè)中的應(yīng)用[J].地震工程與工程振動(dòng),2008,28(1):139-145.(LIU Chun-guang,ZHANG Li-hua.The application of genetic algorithms and BP neural network to seismic damage prediction of bridges[J]. Journal of Earthquake Engineering and Engineering Vibration,2008,28(1):139-145.(in Chinese))

    [40]汪大洋,周云.基于GA優(yōu)化理論的高層結(jié)構(gòu)風(fēng)振控制研究[J].土木工程學(xué)報(bào),2008,41(1):42-47.(WANG Da-yang,ZHOU Yun.Control of wind-induced vibration of high-rise bulidings based onGA optimum theory[J].China Civil Engineering Journal,2008,41 (1):42-47.(in Chinese))

    [41]陳文英,閻紹澤,褚福磊.免疫遺傳算法在智能桁架結(jié)構(gòu)振動(dòng)主動(dòng)控制系統(tǒng)優(yōu)化設(shè)計(jì)中的應(yīng)用[J].機(jī)械工程學(xué)報(bào),2008,44 (2):196-200.(CHEN Wen-ying,YAN Shao-ze,CHU Fu-lei.Immune genetic algorithm used to integrated optimal design of active vibration control system for piezoelectric intelligent truss structures[J].Chinese Journal of Mechanical Engineering,2008,44(2):196-200.(in Chinese))

    [42]楊海霞,謝軼鵬.一種改進(jìn)的遺傳算法及其在結(jié)構(gòu)優(yōu)化設(shè)計(jì)中的應(yīng)用[J].河海大學(xué)學(xué)報(bào):自然科學(xué)版,2007,35(6):659-663.(YANG Hai-xia,XIE Yi-peng.An improved genetic algorithm and its application to structural design optimization[J].Journal of Hohai University:Natural Sciences,2007,35(6):659-663.(in Chinese))

    [43]趙同彬,譚云亮,王虹,等.擋土墻庫(kù)侖土壓力的遺傳算法求解分析[J].巖土力學(xué),2009,30(4):1170-1174.(ZHAO Tong-bin,TAN Yun-liang,WANG Hong,et al.Analytical solution of Coulomb earth pressure on retaining wall by genetic algorithm[J].Rock and Soil Mechanics,2009,30(4):1170-1174.(in Chinese))

    [44]VAN LAAR HOVEN P J M,AARTS E H L.Simulated annealing:theory and applications[M].Berlin:Springer,1987.

    [45]陳昌富,吳子儒,龔曉南.復(fù)合形模擬退火算法及其在水泥土墻優(yōu)化設(shè)計(jì)中的應(yīng)用[J].巖土力學(xué),2007,28(12):2543-2548. (CHEN Chang-fu,WU Zi-ru,Gong Xiao-nan.Complex simulated annealing and its application to optimization design of cement-soil retaining walls[J].Rock and Soil Mechanics,2007,28(12):2543-2548.(in Chinese))

    [46]王永德,李蘭,趙宏才,等.電力網(wǎng)網(wǎng)架結(jié)構(gòu)優(yōu)化設(shè)計(jì)的模擬退火算法[J].青島理工大學(xué)學(xué)報(bào),2007,28(3):108-110,117. (WANG Yong-de,LI Lan,ZHAO Hong-cai,et al.Simulated annealing algorithm for the optimal design of the electric power grid structure [J].Journal of Qingdao Technological University,2007,28(3):108-110,117.(in Chinese))

    [47]高占遠(yuǎn),鹿曉陽(yáng),陳世英,等.基于單純形模擬退火算法的網(wǎng)架結(jié)構(gòu)優(yōu)化設(shè)計(jì)[J].山東建筑大學(xué)學(xué)報(bào),2007,22(1):21-23,40.(GAO Zhan-yuan,LU Xiao-yang,CHEN Shi-ying,et al.Optimal design of lattice grid structure based on simulated simplex-annealing method[J].Journal of Shandong Jianzhu University,2007,22(1):21-23.(in Chinese))

    [48]項(xiàng)寶衛(wèi),胡永良,凌塑勇.桁架形狀優(yōu)化設(shè)計(jì)的模擬退火算法[J].力學(xué)與實(shí)踐,2006,28(6):46-50.(XIANG Bao-wei,HU Yong-liang,LING Su-yong.Shape optimization of trusses using simulated annealing[J].Mechanics in Engineering,2006,28(6):46-50. (in Chinese))

    [49]顧元憲,項(xiàng)寶衛(wèi),趙國(guó)忠.桁架結(jié)構(gòu)截面優(yōu)化設(shè)計(jì)的改進(jìn)模擬退火算法[J].計(jì)算力學(xué)學(xué)報(bào),2006,23(5):546-552.(GU Yuanxian,XIANG Bao-wei,ZHAO Guo-zhong.Improved simulated annealing algorithm for cross-section design optimization of truss tructures [J].Chinese Journal of Computational Mechanics,2006,23(5):546-552.(in Chinese))

    [50]HOPFIELD J J.Neural networks and physical systems with emergent collective computational abilities[J].Proceedings of the National Academy of Sciences of the United States of America,1982,79(8):2554-2558.

    [51]姚宏,劉小文.基于神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)抗震可靠度預(yù)測(cè)及優(yōu)化[J].南昌工程學(xué)院學(xué)報(bào),2007,26(6):22-25.(YAO Hong,LIU Xiao-wen.Prediction and optimization of reliability of seismic structures based on neural network[J].Journal of Nanchang Institute of Technology,2007,26(6):22-25.(in Chinese))

    [52]任遠(yuǎn),白廣忱.新型NN訓(xùn)練算法及其在優(yōu)化設(shè)計(jì)中的應(yīng)用[J].計(jì)算機(jī)工程與應(yīng)用,2009,45(8):5-8.(REN Yuan,BAI Guang-chen.New NN training algorithm and its application in optimization design[J].Computer Engineering and Applications,2009,45 (8):5-8.(in Chinese))

    [53]GLOVER F.Tabu search:Part I[J].ORSA Journal on Computing,1989,1(3):190-206.

    猜你喜歡
    模擬退火遺傳算法變量
    抓住不變量解題
    也談分離變量
    模擬退火遺傳算法在機(jī)械臂路徑規(guī)劃中的應(yīng)用
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    一種基于遺傳算法的聚類(lèi)分析方法在DNA序列比較中的應(yīng)用
    基于遺傳算法和LS-SVM的財(cái)務(wù)危機(jī)預(yù)測(cè)
    基于模糊自適應(yīng)模擬退火遺傳算法的配電網(wǎng)故障定位
    基于改進(jìn)的遺傳算法的模糊聚類(lèi)算法
    SL(3,3n)和SU(3,3n)的第一Cartan不變量
    SOA結(jié)合模擬退火算法優(yōu)化電容器配置研究
    成人无遮挡网站| 美女cb高潮喷水在线观看 | 12—13女人毛片做爰片一| 亚洲午夜精品一区,二区,三区| 国产免费男女视频| 国产激情偷乱视频一区二区| 成人三级做爰电影| 日本撒尿小便嘘嘘汇集6| 成年人黄色毛片网站| 麻豆av在线久日| 免费观看人在逋| 日韩欧美免费精品| 国产成人福利小说| 黄片小视频在线播放| 国产伦精品一区二区三区视频9 | 怎么达到女性高潮| 欧美国产日韩亚洲一区| 首页视频小说图片口味搜索| 午夜精品久久久久久毛片777| 一边摸一边抽搐一进一小说| 日本 欧美在线| 欧美激情在线99| 性色av乱码一区二区三区2| 视频区欧美日本亚洲| АⅤ资源中文在线天堂| 午夜免费观看网址| 欧美在线黄色| 欧美日本亚洲视频在线播放| 2021天堂中文幕一二区在线观| 久久精品91蜜桃| 久久精品91蜜桃| а√天堂www在线а√下载| 亚洲av美国av| 国产欧美日韩一区二区三| 热99在线观看视频| 人人妻人人澡欧美一区二区| 午夜福利高清视频| 国产精品久久久久久久电影 | 国产激情欧美一区二区| 久久伊人香网站| 夜夜爽天天搞| 亚洲激情在线av| 欧美又色又爽又黄视频| 国产爱豆传媒在线观看| 亚洲成人久久性| 亚洲人成网站在线播放欧美日韩| 一本综合久久免费| 国产亚洲精品一区二区www| 少妇的逼水好多| 两个人看的免费小视频| 一本精品99久久精品77| 级片在线观看| 国产av一区在线观看免费| 亚洲国产中文字幕在线视频| 性欧美人与动物交配| 女人被狂操c到高潮| 国产日本99.免费观看| 久久这里只有精品19| 中文字幕精品亚洲无线码一区| 国产淫片久久久久久久久 | 亚洲人与动物交配视频| 欧美xxxx黑人xx丫x性爽| 91久久精品国产一区二区成人 | 一级a爱片免费观看的视频| 2021天堂中文幕一二区在线观| 午夜福利18| av天堂在线播放| 久久国产乱子伦精品免费另类| 国产高清视频在线播放一区| 成人av一区二区三区在线看| 国产高潮美女av| 伊人久久大香线蕉亚洲五| 男女午夜视频在线观看| 嫩草影视91久久| 12—13女人毛片做爰片一| 男人舔女人的私密视频| ponron亚洲| 欧美日韩综合久久久久久 | 亚洲欧美精品综合久久99| 亚洲欧美日韩无卡精品| 免费在线观看亚洲国产| 亚洲 国产 在线| 日韩精品中文字幕看吧| 一进一出抽搐动态| 中文字幕精品亚洲无线码一区| 国产单亲对白刺激| 精品日产1卡2卡| 精品99又大又爽又粗少妇毛片 | 啦啦啦观看免费观看视频高清| 91老司机精品| 亚洲在线观看片| 99国产精品99久久久久| 老司机深夜福利视频在线观看| 一个人看视频在线观看www免费 | 久久精品国产综合久久久| 国产精品亚洲美女久久久| 久久久精品大字幕| 可以在线观看毛片的网站| 天堂av国产一区二区熟女人妻| 欧美乱码精品一区二区三区| 91在线观看av| 国产精品乱码一区二三区的特点| 亚洲色图av天堂| 老司机深夜福利视频在线观看| 婷婷亚洲欧美| 日韩中文字幕欧美一区二区| 国产极品精品免费视频能看的| 全区人妻精品视频| 午夜免费成人在线视频| 日本五十路高清| 午夜免费成人在线视频| 九九久久精品国产亚洲av麻豆 | 色播亚洲综合网| 毛片女人毛片| 很黄的视频免费| 精华霜和精华液先用哪个| 国产精品久久久人人做人人爽| 国产三级黄色录像| 少妇裸体淫交视频免费看高清| 亚洲黑人精品在线| 久久久色成人| bbb黄色大片| 三级毛片av免费| 女人被狂操c到高潮| 国产久久久一区二区三区| 免费观看精品视频网站| xxxwww97欧美| 国产亚洲欧美98| 久久精品91无色码中文字幕| 国产亚洲精品久久久com| 国产成人av教育| 精品无人区乱码1区二区| 成人性生交大片免费视频hd| 久久久久性生活片| 国产单亲对白刺激| 制服人妻中文乱码| 99精品在免费线老司机午夜| 国产午夜精品论理片| 免费看光身美女| 亚洲一区二区三区色噜噜| 国产三级中文精品| 日韩 欧美 亚洲 中文字幕| 欧美在线黄色| 国产私拍福利视频在线观看| 国产三级黄色录像| 日本免费a在线| 一本久久中文字幕| 亚洲avbb在线观看| 国产97色在线日韩免费| 成人18禁在线播放| 欧美一级a爱片免费观看看| 国产成人精品久久二区二区免费| 岛国在线观看网站| 色哟哟哟哟哟哟| 18禁黄网站禁片免费观看直播| 天堂√8在线中文| 一夜夜www| 免费搜索国产男女视频| av女优亚洲男人天堂 | 在线观看免费视频日本深夜| 成人国产一区最新在线观看| 国产精品av视频在线免费观看| 免费高清视频大片| 国产高清有码在线观看视频| 一本久久中文字幕| 国产日本99.免费观看| 成人18禁在线播放| 色综合婷婷激情| 欧美日韩瑟瑟在线播放| 亚洲国产高清在线一区二区三| 成年女人看的毛片在线观看| 亚洲国产精品999在线| 亚洲av成人av| 午夜亚洲福利在线播放| 这个男人来自地球电影免费观看| 麻豆国产av国片精品| 人妻久久中文字幕网| 午夜精品一区二区三区免费看| 欧美绝顶高潮抽搐喷水| 三级国产精品欧美在线观看 | 真人一进一出gif抽搐免费| 国产精品亚洲av一区麻豆| 国产69精品久久久久777片 | 欧美又色又爽又黄视频| 久久国产乱子伦精品免费另类| 草草在线视频免费看| 男女床上黄色一级片免费看| 色综合站精品国产| 亚洲性夜色夜夜综合| 久久久色成人| 久久精品国产99精品国产亚洲性色| 级片在线观看| 欧美日韩国产亚洲二区| 午夜福利在线观看免费完整高清在 | 国产黄色小视频在线观看| 亚洲精品久久国产高清桃花| 亚洲精品456在线播放app | 操出白浆在线播放| 在线观看日韩欧美| 精品熟女少妇八av免费久了| 国产久久久一区二区三区| 午夜日韩欧美国产| 757午夜福利合集在线观看| 精品不卡国产一区二区三区| 又粗又爽又猛毛片免费看| 精品国产三级普通话版| 噜噜噜噜噜久久久久久91| 嫁个100分男人电影在线观看| 国内精品久久久久久久电影| 国产色爽女视频免费观看| 日本熟妇午夜| 神马国产精品三级电影在线观看| 中文精品一卡2卡3卡4更新| 久久韩国三级中文字幕| 干丝袜人妻中文字幕| 波多野结衣巨乳人妻| 国产精品99久久久久久久久| 日韩国内少妇激情av| 久久久亚洲精品成人影院| 国产欧美另类精品又又久久亚洲欧美| 丝袜美腿在线中文| 麻豆精品久久久久久蜜桃| 夜夜爽夜夜爽视频| 精品国产一区二区三区久久久樱花 | 久99久视频精品免费| 国产久久久一区二区三区| 亚洲性久久影院| 欧美一区二区精品小视频在线| 国产高清视频在线观看网站| 性色avwww在线观看| 亚洲天堂国产精品一区在线| 久久久久国产网址| 亚洲乱码一区二区免费版| 男女边吃奶边做爰视频| 毛片女人毛片| 乱码一卡2卡4卡精品| 99视频精品全部免费 在线| 欧美日韩精品成人综合77777| or卡值多少钱| 国产淫片久久久久久久久| 国产成人精品一,二区| 亚洲电影在线观看av| 日日撸夜夜添| 色播亚洲综合网| 亚洲av电影不卡..在线观看| 伦精品一区二区三区| 亚洲激情五月婷婷啪啪| 长腿黑丝高跟| 午夜精品在线福利| 在线播放国产精品三级| 国产毛片a区久久久久| 日本午夜av视频| 国产高清国产精品国产三级 | 国产成人a∨麻豆精品| 免费看日本二区| 人人妻人人看人人澡| 国产精品一区二区性色av| 七月丁香在线播放| 久久午夜福利片| 欧美成人午夜免费资源| 成人高潮视频无遮挡免费网站| 22中文网久久字幕| 国产成人91sexporn| 欧美激情在线99| 一级av片app| 搞女人的毛片| 日本五十路高清| 噜噜噜噜噜久久久久久91| 午夜激情欧美在线| 一级爰片在线观看| 色网站视频免费| 精品久久久久久电影网 | 亚洲最大成人中文| 国产精品福利在线免费观看| 中文字幕熟女人妻在线| 国产精品一区www在线观看| 国产综合懂色| 欧美区成人在线视频| 亚洲av.av天堂| 午夜福利在线观看吧| 国产在线一区二区三区精 | 亚洲成人中文字幕在线播放| 中国国产av一级| 久久人人爽人人片av| 日韩欧美在线乱码| 亚洲av中文av极速乱| 91久久精品国产一区二区三区| 国产精品美女特级片免费视频播放器| 精品人妻一区二区三区麻豆| 小说图片视频综合网站| 免费电影在线观看免费观看| videossex国产| 丝袜喷水一区| 插逼视频在线观看| 婷婷六月久久综合丁香| 天天躁日日操中文字幕| 男人狂女人下面高潮的视频| 午夜视频国产福利| 久久精品人妻少妇| 欧美一区二区国产精品久久精品| 男人和女人高潮做爰伦理| 国产精品麻豆人妻色哟哟久久 | 美女脱内裤让男人舔精品视频| 久久久久久伊人网av| 午夜福利视频1000在线观看| 日本一本二区三区精品| 久久精品91蜜桃| 长腿黑丝高跟| 国产精品熟女久久久久浪| 最近最新中文字幕大全电影3| 国产熟女欧美一区二区| 亚洲一级一片aⅴ在线观看| 69av精品久久久久久| 亚洲国产精品久久男人天堂| 日本熟妇午夜| 日本三级黄在线观看| 99视频精品全部免费 在线| 亚洲综合色惰| 国产在视频线精品| 午夜久久久久精精品| 日本免费在线观看一区| 色吧在线观看| 亚洲欧美日韩高清专用| 日韩一区二区视频免费看| 亚洲精品亚洲一区二区| 听说在线观看完整版免费高清| 桃色一区二区三区在线观看| 日韩国内少妇激情av| 婷婷六月久久综合丁香| 禁无遮挡网站| 秋霞伦理黄片| 免费看a级黄色片| 国产黄色小视频在线观看| 中文在线观看免费www的网站| 永久网站在线| 亚洲精品久久久久久婷婷小说 | 日本欧美国产在线视频| 国产成人a区在线观看| 国内精品美女久久久久久| 听说在线观看完整版免费高清| 久久久国产成人精品二区| 黄色配什么色好看| 久久久久久伊人网av| 99久国产av精品国产电影| 少妇的逼水好多| 亚洲精品日韩av片在线观看| 一夜夜www| 免费黄网站久久成人精品| 亚洲av日韩在线播放| 一个人免费在线观看电影| 久久人妻av系列| 日本与韩国留学比较| 日本一二三区视频观看| 亚洲欧美精品自产自拍| 99在线视频只有这里精品首页| 中文亚洲av片在线观看爽| 一区二区三区免费毛片| 人体艺术视频欧美日本| 国产亚洲精品久久久com| 亚洲中文字幕日韩| 亚洲欧美日韩高清专用| 精品人妻偷拍中文字幕| 免费在线观看成人毛片| 91aial.com中文字幕在线观看| 日本黄色片子视频| 亚洲成人中文字幕在线播放| 免费黄网站久久成人精品| 午夜a级毛片| 精品久久久久久成人av| 久久久成人免费电影| 国产在线一区二区三区精 | 91精品国产九色| 日韩欧美精品v在线| 国产单亲对白刺激| 久久久久精品久久久久真实原创| 国产av一区在线观看免费| 色综合亚洲欧美另类图片| 国产久久久一区二区三区| 国产伦精品一区二区三区视频9| 亚洲精品乱码久久久v下载方式| 婷婷六月久久综合丁香| 国产大屁股一区二区在线视频| 国产私拍福利视频在线观看| 夜夜爽夜夜爽视频| 欧美一区二区国产精品久久精品| 国产伦理片在线播放av一区| 日本爱情动作片www.在线观看| 国产淫语在线视频| 一边亲一边摸免费视频| 日韩欧美在线乱码| 国产黄色小视频在线观看| 亚洲最大成人中文| 18禁在线无遮挡免费观看视频| 免费av不卡在线播放| 日韩欧美精品免费久久| av在线蜜桃| 亚洲欧美日韩高清专用| 国产三级中文精品| 成人二区视频| 久久精品夜夜夜夜夜久久蜜豆| 精品99又大又爽又粗少妇毛片| 亚洲av成人精品一二三区| 99久久成人亚洲精品观看| 久久久久久伊人网av| 国产高清三级在线| 国产亚洲午夜精品一区二区久久 | 国产综合懂色| 变态另类丝袜制服| 日日摸夜夜添夜夜添av毛片| 亚洲欧美日韩卡通动漫| 三级男女做爰猛烈吃奶摸视频| 丰满少妇做爰视频| 永久免费av网站大全| 三级毛片av免费| 2021天堂中文幕一二区在线观| 日本黄色片子视频| 免费播放大片免费观看视频在线观看 | 国产成人午夜福利电影在线观看| 国产黄色小视频在线观看| 国产精品国产高清国产av| 六月丁香七月| 国内精品美女久久久久久| 亚洲成人久久爱视频| 国产av在哪里看| 蜜桃亚洲精品一区二区三区| 婷婷六月久久综合丁香| 精品久久久噜噜| 久久久精品94久久精品| 不卡视频在线观看欧美| 欧美xxxx黑人xx丫x性爽| 免费观看在线日韩| 亚洲欧美一区二区三区国产| 久久国产乱子免费精品| av黄色大香蕉| 一边亲一边摸免费视频| av在线亚洲专区| 两个人视频免费观看高清| 国产精品.久久久| 亚洲综合精品二区| 日韩一本色道免费dvd| 欧美丝袜亚洲另类| 美女大奶头视频| 高清毛片免费看| 黑人高潮一二区| 看片在线看免费视频| 一级二级三级毛片免费看| www.av在线官网国产| 一级毛片我不卡| 日韩欧美 国产精品| 国产精品一二三区在线看| 少妇被粗大猛烈的视频| 男女那种视频在线观看| 色综合站精品国产| 久久鲁丝午夜福利片| 精品久久久久久久末码| 一区二区三区高清视频在线| 成人鲁丝片一二三区免费| 亚洲美女视频黄频| 国产探花在线观看一区二区| 欧美一区二区亚洲| 成人亚洲欧美一区二区av| 精品一区二区三区视频在线| 中文欧美无线码| av线在线观看网站| 久久精品影院6| 中文资源天堂在线| 色综合色国产| 日韩高清综合在线| 日韩av在线免费看完整版不卡| 超碰97精品在线观看| 九九久久精品国产亚洲av麻豆| 天堂网av新在线| 最近视频中文字幕2019在线8| 黄色一级大片看看| 国内精品宾馆在线| 国产av在哪里看| 蜜桃久久精品国产亚洲av| 国产精品人妻久久久影院| 亚洲精华国产精华液的使用体验| 丝袜美腿在线中文| 麻豆久久精品国产亚洲av| 久久99热6这里只有精品| 久久这里只有精品中国| 一级毛片aaaaaa免费看小| 免费观看人在逋| 亚洲成人av在线免费| 熟女人妻精品中文字幕| 婷婷色综合大香蕉| 国产成人a区在线观看| 日韩强制内射视频| 亚洲经典国产精华液单| 国产淫片久久久久久久久| 国产午夜精品久久久久久一区二区三区| av线在线观看网站| 久久久久久久午夜电影| 伦理电影大哥的女人| 精品欧美国产一区二区三| 国产三级在线视频| 老师上课跳d突然被开到最大视频| 国产精品日韩av在线免费观看| 久久鲁丝午夜福利片| 毛片女人毛片| 国产精品无大码| 乱码一卡2卡4卡精品| 男女国产视频网站| 一个人看的www免费观看视频| 在线免费十八禁| 免费观看a级毛片全部| 国产精品麻豆人妻色哟哟久久 | 欧美丝袜亚洲另类| 七月丁香在线播放| 久久亚洲国产成人精品v| 日韩欧美精品免费久久| 2022亚洲国产成人精品| 国产乱人偷精品视频| 一级毛片电影观看 | 18禁在线播放成人免费| 欧美区成人在线视频| 精品无人区乱码1区二区| 亚洲av成人精品一二三区| 秋霞伦理黄片| 免费搜索国产男女视频| 日本免费一区二区三区高清不卡| 欧美成人免费av一区二区三区| 亚洲国产欧美在线一区| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲av免费高清在线观看| 日韩欧美精品v在线| 国产色爽女视频免费观看| 国产精品女同一区二区软件| 日本黄色视频三级网站网址| 成人鲁丝片一二三区免费| 国产精品,欧美在线| 可以在线观看毛片的网站| 91在线精品国自产拍蜜月| 国产成人精品久久久久久| 久久久国产成人精品二区| 爱豆传媒免费全集在线观看| 国产免费男女视频| 久久久色成人| 日韩中字成人| av专区在线播放| 久久韩国三级中文字幕| 又粗又爽又猛毛片免费看| av线在线观看网站| 激情 狠狠 欧美| 九九在线视频观看精品| 美女高潮的动态| 黑人高潮一二区| 两性午夜刺激爽爽歪歪视频在线观看| 日韩亚洲欧美综合| 看非洲黑人一级黄片| 日韩一区二区视频免费看| 少妇猛男粗大的猛烈进出视频 | 成年女人看的毛片在线观看| 亚洲精品456在线播放app| 亚洲av成人av| 女人久久www免费人成看片 | 亚洲欧美成人综合另类久久久 | 嫩草影院入口| 亚洲精品久久久久久婷婷小说 | 欧美bdsm另类| 成人午夜精彩视频在线观看| 欧美成人一区二区免费高清观看| 免费看a级黄色片| 99热这里只有是精品50| 亚洲av熟女| 男女啪啪激烈高潮av片| 级片在线观看| 在线观看美女被高潮喷水网站| 国内精品美女久久久久久| 欧美极品一区二区三区四区| 欧美一区二区国产精品久久精品| 欧美日韩综合久久久久久| 国产69精品久久久久777片| 五月玫瑰六月丁香| 少妇熟女aⅴ在线视频| 水蜜桃什么品种好| 国产精品久久久久久精品电影| 日韩在线高清观看一区二区三区| 国产精华一区二区三区| 日韩中字成人| 久久久久久久久中文| 午夜免费激情av| 亚洲成人中文字幕在线播放| a级一级毛片免费在线观看| 午夜老司机福利剧场| 少妇人妻一区二区三区视频| 亚洲av成人精品一区久久| 日日摸夜夜添夜夜爱| 波多野结衣高清无吗| 十八禁国产超污无遮挡网站| av视频在线观看入口| 日韩精品青青久久久久久| 真实男女啪啪啪动态图| 亚洲欧美日韩卡通动漫| 波多野结衣高清无吗| 日本一本二区三区精品| 丝袜美腿在线中文| 日产精品乱码卡一卡2卡三| 亚洲在线观看片| 男女视频在线观看网站免费| 亚洲精品,欧美精品| 亚洲在线自拍视频| 免费搜索国产男女视频| 人人妻人人澡欧美一区二区| 国产探花在线观看一区二区| 午夜福利在线观看吧| 久久国产乱子免费精品| 波多野结衣巨乳人妻| 99久久精品国产国产毛片| 午夜a级毛片| 国产伦一二天堂av在线观看| 午夜爱爱视频在线播放| 国产成人免费观看mmmm|