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

    基于Chebyshev零點多項式區(qū)間不確定的可靠性拓?fù)鋬?yōu)化設(shè)計

    2020-06-17 05:30:38蘇海亮蘭鳳崇賀裕雁陳吉清
    關(guān)鍵詞:橢球不確定性區(qū)間

    蘇海亮 蘭鳳崇 賀裕雁 陳吉清?

    (1.華南理工大學(xué) 機(jī)械與汽車工程學(xué)院,廣東 廣州 510640;2.華南理工大學(xué) 廣東省汽車工程重點實驗室,廣東 廣州 510640;3.華南理工大學(xué) 工商管理學(xué)院,廣東 廣州 510640)

    拓?fù)鋬?yōu)化方法為探索理想和優(yōu)化的結(jié)構(gòu)提供了可能,并在諸多領(lǐng)域得到了廣泛的應(yīng)用。對于結(jié)構(gòu)拓?fù)鋬?yōu)化而言,可以將結(jié)構(gòu)視為離散體或連續(xù)體結(jié)構(gòu),其本質(zhì)是通過迭代計算,在給定的約束條件下以最優(yōu)化性能的方式計算設(shè)計域內(nèi)材料的最佳布局[1]。在過去的幾十年里,關(guān)于連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化的研究大多建立在所有參數(shù)都是確定性的這一假設(shè)的基礎(chǔ)上,而在設(shè)計的整個生命周期中,材料參數(shù)(楊氏模量、許用應(yīng)力等)、外載荷、制造誤差、邊界條件和幾何尺寸等設(shè)計變量都具有結(jié)構(gòu)固有的本質(zhì)不確定性[2]。確定性假設(shè)可能導(dǎo)致設(shè)計不能滿足預(yù)期的性能目標(biāo),最終導(dǎo)致整體結(jié)構(gòu)設(shè)計不可行。因此,為了對工程系統(tǒng)進(jìn)行真實的可靠性評估,在優(yōu)化過程的數(shù)值建模中,定量考慮不確定性是很有必要的。

    國內(nèi)外學(xué)者對于不確定參數(shù)建模及可靠性評估進(jìn)行了大量的研究,這種尋找不確定性的優(yōu)化設(shè)計被稱為基于可靠性的優(yōu)化設(shè)計(RBDO)[3]?;诳煽啃缘脑O(shè)計優(yōu)化理念是一種在考慮不確定性影響的情況下確定最優(yōu)設(shè)計方案的工程設(shè)計方法,它將不確定性作為一個量化的失效概率方式,優(yōu)化結(jié)構(gòu)性能指標(biāo),保持各種設(shè)計約束滿足預(yù)期的可靠性水平。在傳統(tǒng)的可靠性分析方法中,不確定參數(shù)通常被設(shè)定為隨機(jī)變量處理,其概率分布是基于完整數(shù)據(jù)信息定義的。然而,由于工程問題存在的多樣復(fù)雜性,確定某種精確的概率分布形式非常費時,甚至是一項不可能完成的任務(wù)。在這樣的情況下,可以使用基于非概率概念的分析方法來量化和處理這種不確定性,例如目前研究較多的凸模型、區(qū)間模型、模糊集理論等[4- 6]。區(qū)間分析作為一種特殊的凸集模型已成功地應(yīng)用于一系列具有不確定但有界參數(shù)的工程設(shè)計問題中,近年來在學(xué)術(shù)界和工程界引起了越來越多的關(guān)注[7- 9]?;趨^(qū)間分析的非概率可靠度方法在測量、建模思想、評估方式等方面與傳統(tǒng)的結(jié)構(gòu)概率可靠度方法類似,已有相關(guān)學(xué)者對此進(jìn)行了多方面的比較分析。Ben-Haim等[6]于1990年提出了基于凸模型的結(jié)構(gòu)非概率可靠度概念,此后,基于凸模型的非概率可靠性設(shè)計優(yōu)化作為可選擇的RBDO方法得到了廣泛應(yīng)用。Kang等[10]提出了一種基于性能的方法,這種方法避免了基于多橢球凸模型的方法在非概率可靠性指標(biāo)求解過程中所需的冗長計算時間。Jiang等[11]基于凸模型提出一種不確定性的優(yōu)化方法,這種方法將約束以顯式形式來表達(dá)復(fù)雜的結(jié)構(gòu)問題,能夠有效地計算約束區(qū)間。Au等[12]利用凸模型提出一種結(jié)構(gòu)魯棒可靠性方法,引入不滿意度函數(shù)將各種不確定變量映射為無量綱量來實現(xiàn)基準(zhǔn)測試,通過最小化最差的不滿意度來求解目標(biāo)函數(shù)的最優(yōu)值。Meng等[13]提出一種基于單橢球凸模型的混合非概率方法,研究了單橢球凸模型的搜索過程,并利用粒子群優(yōu)化算法來保證全局優(yōu)化設(shè)計。

    針對拓?fù)鋬?yōu)化問題涉及不確定性參數(shù)的情形,基于數(shù)值理論的可靠性拓?fù)鋬?yōu)化設(shè)計方法(RTOD)被提出。Eigel等[14]考慮隨機(jī)材料屬性和隨機(jī)力,建立了一種不確定性情況下的結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計方法,采用相場方法使得材料分布可以在迭代優(yōu)化期間有任意拓?fù)渥兓?。Jalalpour等[15]基于隨機(jī)擾動方式估計不確定性參數(shù)結(jié)構(gòu)系統(tǒng)的失效模式,提出一種基于可靠性的拓?fù)鋬?yōu)化算法。Liu等[16]考慮了載荷方向不確定性的概率性和模糊性質(zhì),并應(yīng)用云模型描述這種不確定性,提出了體積約束的不確定拓?fù)鋬?yōu)化方法。Keshavarzzadeh等[17]將非侵入式多項式混沌擴(kuò)展與可靠性相結(jié)合,提出一種在不確定性下進(jìn)行拓?fù)鋬?yōu)化的系統(tǒng)方法,從而準(zhǔn)確且有效地估計統(tǒng)計矩、失效概率及其靈敏度。Luo等[18- 20]提出了一種基于非概率多橢球凸模型的拓?fù)鋬?yōu)化方法,該模型表示有界參數(shù)的不確定性。McWilliam[21]采用區(qū)間有理系列展開策略來評估區(qū)間位移和應(yīng)力的下界和上界,從而進(jìn)行拓?fù)鋬?yōu)化設(shè)計。Wu等[22]將Chebyshev區(qū)間方法應(yīng)用于具有區(qū)間不確定性的車輛動力學(xué)優(yōu)化,并通過案例證明該方法可進(jìn)行有效估計,并避免優(yōu)化中出現(xiàn)嵌套雙環(huán)。雖然RTOD是一個具有學(xué)術(shù)研究價值且迅速發(fā)展的活躍領(lǐng)域,但具有區(qū)間參數(shù)的約束函數(shù)的拓?fù)鋬?yōu)化集成仍是當(dāng)前值得探索且極具挑戰(zhàn)性的問題。目前,RTOD在工程上的應(yīng)用相對比較少,其主要限制之一是優(yōu)化過程中難以保證計算精度與計算效率的平衡。RTOD的實質(zhì)是一種內(nèi)外層嵌套優(yōu)化問題,可靠性分析嵌入優(yōu)化循環(huán)中,優(yōu)化過程中需要多次計算約束函數(shù)的可靠度,致使整個優(yōu)化過程的計算效率較低,導(dǎo)致無法在實際工程中應(yīng)用。

    鑒于此,文中在傳統(tǒng)的拓?fù)鋬?yōu)化設(shè)計基礎(chǔ)上,基于變密度法提出了一種有效的計算方法來解決結(jié)構(gòu)可靠性拓?fù)鋬?yōu)化中的重要問題:找到使剛度和質(zhì)量達(dá)到最佳平衡的材料分布,這可以表示為受應(yīng)變能約束的最小化體積(最輕量化)。針對實際應(yīng)用中難以獲得足夠的不確定參數(shù)信息的情況,應(yīng)用區(qū)間變量來表示文中所描述的參數(shù)不確定性,采用非概率凸模型建模方法研究具有區(qū)間特性的設(shè)計參數(shù)與結(jié)構(gòu)響應(yīng)之間的映射關(guān)系??紤]到Chebyshev多項式在區(qū)間范圍內(nèi)具有最佳一致逼近的性質(zhì),利用Chebyshev級數(shù)逼近多項式函數(shù)求解函數(shù)關(guān)系,建立Chebyshev多項式模型。通過引進(jìn)區(qū)間數(shù)定義進(jìn)一步計算結(jié)構(gòu)的可靠性指標(biāo),采用分解技術(shù)來提高計算效率。通過內(nèi)循環(huán)是單獨進(jìn)行可靠性分析、外循環(huán)是考慮可靠性的拓?fù)鋬?yōu)化設(shè)計這種雙層分離尋優(yōu)策略實現(xiàn)對考慮不確定性的結(jié)構(gòu)的拓?fù)鋬?yōu)化研究。為了驗證所提方法的可行性和適用性,文中還通過二維和三維模型算例進(jìn)行了探討。

    1 基于凸模型的非概率可靠性指標(biāo)定義

    (1)

    (2)

    式中,Xi為第i個不確定變量X的歸一化變量,w表示具不確定性的參數(shù)。

    超橢球(凸)模型思想是利用超橢球定義不確定參數(shù)的波動區(qū)域,因其不需要精確的信息就能定義不確定性而被廣泛使用。文獻(xiàn)[24- 25]中總結(jié)了超橢球面的測定方法,同時也說明了區(qū)間變量與凸集模型之間的聯(lián)系。二維超橢球模型如圖1所示。根據(jù)文獻(xiàn)[19- 20]提出的橢球凸模型非概率可靠度指標(biāo)的定義,所有不確定參數(shù)的可能值都可以假設(shè)為一個多維(超)橢球。因此,不確定參數(shù)向量w的界限范圍用一個超橢球集合來界定,表示為

    w∈(Q,ε)={w:wTQw≤ε2}

    (3)

    式中,Q為超橢球面的特征矩陣,ε為定義參數(shù)變異性大小的實數(shù)。

    圖1 二維超橢球模型Fig.1 Two-dimension super ellipsoid model

    假設(shè)不確定參數(shù)采用k個橢球體集合表示,則

    (4)

    對于每組參數(shù)wi(i=1,2,…,k),分別由一組橢球無量綱形式有界表示為

    wi∈

    (5)

    通過引入向量定義

    (6)

    (7)

    式中,E為標(biāo)準(zhǔn)空間下的單位超球體(如圖2所示),q為不確定參數(shù)在標(biāo)準(zhǔn)空間的標(biāo)準(zhǔn)化向量。經(jīng)過式(7)所表示的歸一化后,極限狀態(tài)函數(shù)g(X)同樣可以轉(zhuǎn)變?yōu)間(q)。此時將q空間劃分為一個安全區(qū)域(對應(yīng)g(q)> 0)和一個失效區(qū)域(對應(yīng)g(q)<0),則極限狀態(tài)函數(shù)的失效概率Pf同樣可以表示為

    Pf=Pr{g(q)≤0}

    (8)

    圖2 單橢球模型非概率可靠性指標(biāo)的定義

    Fig.2 Definition of non-probabilistic reliability index of single ellipsoid model

    根據(jù)可靠性指標(biāo)的意義,定義結(jié)構(gòu)單橢球模型的可靠性指標(biāo)γ為

    (9)

    式中,sgn(υ)是實數(shù)υ的sgn函數(shù),定義為

    (10)

    之所以存在sgn函數(shù),是因為當(dāng)不確定性變量對應(yīng)的名義值處所求極限狀態(tài)函數(shù)為負(fù)時,產(chǎn)生的可靠性指標(biāo)應(yīng)為負(fù)值。圖3展示了如何用單橢球模型來描述參數(shù)的不確定性。在變換后的標(biāo)準(zhǔn)空間中,不確定性參數(shù)以一個單位球為界來呈現(xiàn)。圖中單位圓能夠很好地反映這種參數(shù)不確定結(jié)構(gòu)的可靠度。根據(jù)圖3所示,從原點到極限狀態(tài)曲面上任意點的向量的最小無窮大范數(shù)能夠度量結(jié)構(gòu)的非概率可靠性,因此,非概率可靠性指標(biāo)γ定義為

    γ=sgn(g(0))·min(max(|q1|,|q2|))

    (11)

    s.t.g(q)=0。

    式中,qi(i=1,2)為第i組不確定變量wi的歸一化向量。根據(jù)標(biāo)準(zhǔn)多橢球集合的邊界形狀,將其推廣到多橢球凸模型,即存在多橢球體模型的k組不確定性。由此,將γ推廣至多維空間下的變量可靠性指標(biāo),表示為

    (12)

    s.t.gj(q)=0。

    式中的γ受制于對應(yīng)的極限狀態(tài)函數(shù)。式(12)表明,可靠性指標(biāo)γ=1時,極限狀態(tài)曲面與單位圓凸模型存在著相交邊界,這種相交意味著凸模型與設(shè)定結(jié)構(gòu)波動性的極限狀態(tài)函數(shù)處于一種臨界狀態(tài)。γ>1時,所有存在的不確定性變量都處于可靠的結(jié)構(gòu)域和相應(yīng)的安全裕度,此時結(jié)構(gòu)被認(rèn)為是可靠的。γ值越大表明結(jié)構(gòu)將會被允許有更大程度的參數(shù)波動。

    圖3 兩個區(qū)間參數(shù)的不確定可靠性量化標(biāo)準(zhǔn)

    Fig.3 Reliability quantification criteria for two interval uncertain parameters

    2 區(qū)間Chebyshev多項式逼近函數(shù)模型

    極限狀態(tài)函數(shù)的定義在可靠性評估中非常重要,由于結(jié)構(gòu)拓?fù)鋬?yōu)化分析中的目標(biāo)函數(shù)或約束函數(shù)不能明確表達(dá),導(dǎo)致進(jìn)行RTOD分析時會產(chǎn)生高昂的費用甚至無法實施這樣的分析。因此,可考慮在優(yōu)化過程中使用精確的多項式函數(shù)來顯性表達(dá)性能與變量之間的關(guān)系。變量在某一范圍內(nèi)波動的目標(biāo)函數(shù)的表達(dá)同樣是一個難點?;趨^(qū)間概念并考慮計算效率,文中使用Chebyshev多項式擬合變量與性能之間的隸屬關(guān)系來實現(xiàn)精確的區(qū)間函數(shù)擬合。文獻(xiàn)[26]表明,Chebyshev多項式在區(qū)間函數(shù)逼近方面比泰勒展開式具有更高的精度。

    考慮到Chebyshev多項式在區(qū)間范圍內(nèi)具有最佳一致逼近的性質(zhì),假設(shè)存在一維函數(shù),該函數(shù)可近似通過一維k階Chebyshev多項式逼近:

    (13)

    式中:pk為k階Chebyshev多項式;f表示多項式的系數(shù),其計算公式為

    (14)

    m為插值點數(shù)目;θj=arccosxj,x為變量;Cj(x)為構(gòu)造Chebyshev多項式的向量,定義為

    Cj(x)=cos(jθ)

    (15)

    (16)

    式中,xL和xR分別表示變量x的上界和下界,x∈[-1,1],設(shè)xj為Chebyshev插值點,其可以通過下式計算得到:

    (17)

    處理多維不確定參數(shù)問題時,引入張量積可以直接構(gòu)造多維Chebyshev插值點。如一維插值點節(jié)點基為X={x1,x2,…,xm},為了平衡效率與精度的關(guān)系,插值點數(shù)目m=ki+1,基于張量積構(gòu)造的L維插值多項式節(jié)點集合為X=X1?X2?…?Xl,從而可知,L維問題的插值點總數(shù)Ntotal應(yīng)為

    Ntotal=(k1+1)×(k2+1)×…×(kL+1)

    (18)

    因此,Chebyshev級數(shù)略去高階的k階多項式可近似表示為

    (19)

    式中,r表示下標(biāo)等于0的總數(shù),Ci1…ik為k維Chebyshev多項式向量,fi1…ik為k維系數(shù)。

    對每一維積分采用插值積分公式計算,則Chebyshev多項式的每個系數(shù)都可通過相應(yīng)的積分公式計算得到:

    (20)

    式中,k為數(shù)值積分公式的階數(shù),每個維的插值點都是k+1次Chebyshev多項式的零點。為了計算的便利性,同樣可以用最小二乘法來計算Chebyshev多項式的系數(shù)[22],將式(19)轉(zhuǎn)化為

    (21)

    β=[β1…βs]T=[f0,…,0…fii,…,ik]T

    (22)

    δ=[δ1…δs]T=[C0,…,0…Cii,…,ik]T

    (23)

    式中,0≤i1+…+ik≤n,s表示系數(shù)向量中系數(shù)的數(shù)量,也是多項式基向量δ中Chebyshev多項式的數(shù)量。拓?fù)鋬?yōu)化設(shè)計由于其域的幾何形狀變化復(fù)雜,概率約束往往是隱式的,計算起來也比較困難,而采用逼近函數(shù)方法進(jìn)行近似,計算起來就容易得多,成本也低。通過Chebyshev代理模型替換有限元模型進(jìn)行計算,僅需采用有限的樣本數(shù)據(jù)構(gòu)造Chebyshev多項式,就能夠針對RTOD進(jìn)行可靠性分析。

    3 基于非概率可靠性的拓?fù)鋬?yōu)化設(shè)計

    3.1 確定性拓?fù)鋬?yōu)化

    連續(xù)體拓?fù)鋬?yōu)化問題通常是通過有限元方法求解,其中元素體積分?jǐn)?shù)用于指示每一個元素域內(nèi)的材料分布。在一個典型的拓?fù)鋬?yōu)化問題中,目標(biāo)通常是通過確定材料的分布來獲得性能最佳的結(jié)構(gòu)。因此,文中所建立的拓?fù)鋬?yōu)化設(shè)計模型表達(dá)為在剛度約束下最小化體積V,以應(yīng)變能作為約束的形式,其表達(dá)式為

    (24)

    s.t.K(ρ)d(ρ)=F;

    gi(ρ)=CT-Ci(ρ)≥0,i=1,2,…,N;

    0≤ρmin≤ρe≤1。

    上述公式中,ρ為材料密度,V為結(jié)構(gòu)材料體積,ve為單元體積,ρe為單元密度,K(ρ)、d(ρ)、F為整體剛度矩陣、位移矩陣以及有限元離散得到的力向量,N為位移約束的個數(shù),CT為ρ在第i個應(yīng)變能約束Ci下的上限值,ρmin為單元密度的下限。使用常用的固體各向同性材料懲罰(SIMP)方法[27]懲罰單元設(shè)計域中體積分?jǐn)?shù)和驅(qū)動單元0-1分布進(jìn)行求解。在SIMP[28]插值下,單元剛度矩陣為Ke表示如下:

    (25)

    3.2 非概率可靠性拓?fù)鋬?yōu)化設(shè)計

    利用區(qū)間變量不確定性的概率定義,基于凸模型非概率可靠性指標(biāo),將可靠性約束作為結(jié)構(gòu)拓?fù)鋬?yōu)化問題的約束條件,以應(yīng)變能作為約束問題的一般RTOD公式的形式如下:

    (26)

    s.t.K(ρ)d(ρ)=F;

    0≤ρmin≤ρe≤1。

    γ=sgn(gj(0))·min{max(|q1|,|q2|,…,|qk|)}

    (27)

    s.t.gj(q)=0。

    γ越大意味著結(jié)構(gòu)允許的參數(shù)變化越大。

    3.3 基于Chebyshev多項式的可靠性拓?fù)鋬?yōu)化

    本節(jié)基于Chebyshev代理模型來解決上述區(qū)間不確定情形下的可靠性拓?fù)鋬?yōu)化問題。在優(yōu)化過程中,將Chebyshev零點作為不確定設(shè)計向量的初始樣本,然后利用Chebyshev多項式模型建立不確定參數(shù)約束條件的近似模型。因此,式(26)可轉(zhuǎn)變?yōu)橄率剑?/p>

    (28)

    (29)

    (30)

    s.t.wTw≤(γ*)2。

    當(dāng)區(qū)間極限狀態(tài)函數(shù)引起的目標(biāo)可靠性指標(biāo)存在區(qū)間范圍時,式(30)同樣可以描述為

    (31)

    圖4 區(qū)間約束可能的取值范圍Fig.4 Possible value range of interval constraint

    采用Karush-Kuhn-Tucker(KKT)最優(yōu)性條件求解約束優(yōu)化問題,于是有

    (32)

    (33)

    s.t.Zk(u*(ρ),x*(ρ))≥0,k=1,2,…,m。

    式中:Zk為目標(biāo)性能值,對應(yīng)于規(guī)定的目標(biāo)可靠性指標(biāo)γ;k代表第k種失效模式。因此,約束條件關(guān)于最佳設(shè)計點的問題可以通過下面這個最小化問題求解:

    (34)

    對于這種最小化問題,可通過可靠性優(yōu)化設(shè)計中比較常見的性能度量方法求解。

    在求解RTOD問題的過程中,需構(gòu)建Chebyshev多項式模型,包括兩個主要步驟:獲取采樣所需的數(shù)據(jù),計算模型的系數(shù)。采樣點是Chebyshev多項式的零點,其數(shù)量主要由區(qū)間變量的維數(shù)和所構(gòu)建多項式的階數(shù)來決定。在得到足夠采樣數(shù)據(jù)的基礎(chǔ)上,可利用最小二乘法求解系數(shù),之后構(gòu)建函數(shù)的替代模型。Chebyshev代理模型能夠與外環(huán)優(yōu)化相結(jié)合,實現(xiàn)不確定優(yōu)化的求解(式(34)的解析)?;趨^(qū)間Chebyshev模型函數(shù)的主要優(yōu)點是所建立的約束函數(shù)的近似模型精度高,優(yōu)化結(jié)果更加可靠。

    3.4 迭代機(jī)制流程

    可靠性分析方法是指用于計算系統(tǒng)響應(yīng)值在參數(shù)波動影響下滿足某種約束條件的可靠度的方法。文中的拓?fù)鋬?yōu)化設(shè)計是以約束條件滿足給定的可靠性指標(biāo)為前提的搜索尋優(yōu)問題,在優(yōu)化方程(33)中,將不確定參數(shù)確定為一種目標(biāo)非概率可靠性的確定性形式,所得結(jié)果用作拓?fù)鋬?yōu)化的輸入變量進(jìn)行新的確定性拓?fù)鋬?yōu)化設(shè)計,這樣設(shè)計的結(jié)果可滿足結(jié)構(gòu)的可靠度設(shè)計要求。其主要步驟如下:

    步驟1 針對結(jié)構(gòu)建立有限元模型,定義加載工況,構(gòu)建非概率凸模型;

    步驟2 進(jìn)行單元敏感度分析,基于Chebyshev零點原理構(gòu)造初始樣本點;

    步驟3 根據(jù)樣本點進(jìn)行有限元分析,建立Chebyshev多項式模型;

    步驟4 利用Chebyshev多項式模型結(jié)合非概率凸模型進(jìn)行可靠性分析,求解式(34)使用式(32)的迭代機(jī)制,得出給定目標(biāo)可靠性指標(biāo)下的性能最佳設(shè)計參數(shù);

    步驟5 將步驟4所得參數(shù)值返回拓?fù)鋬?yōu)化設(shè)計模型,建立新的拓?fù)鋬?yōu)化模型來尋找最佳的材料分布,直到滿足拓?fù)鋬?yōu)化收斂條件為止,最終獲得RBDO拓?fù)鋯栴}的最優(yōu)解。

    基于非概率可靠性的優(yōu)化設(shè)計是一種嵌套的雙優(yōu)化問題,其中可靠性分析和優(yōu)化是嵌套的。文中使用分離解耦方式,使得評估約束的可靠性與最小化目標(biāo)函數(shù)單獨進(jìn)行,從而精確、有效地得到不確定性結(jié)構(gòu)的最優(yōu)設(shè)計方案。

    4 算例分析

    在算例分析中,都是假定材料的楊氏模量在整個域中是不確定的。設(shè)計目標(biāo)是使受應(yīng)變能約束的結(jié)構(gòu)質(zhì)量最小化,其中包含不同可靠性指標(biāo)下的材料體積最小值。文中分別提供了二維和三維的數(shù)值算例來說明所提方法的有效性和適用性,以及在結(jié)構(gòu)拓?fù)鋬?yōu)化中考慮不確定性的重要性。

    4.1 算例1

    圖5 L形結(jié)構(gòu)的設(shè)計區(qū)域Fig.5 Design field of L-shaped structure

    設(shè)計目標(biāo)是在應(yīng)變能約束下找到使材料總體積最小的最優(yōu)拓?fù)湫螤睢2捎萌AChebyshev多項式擬合變量與應(yīng)變能之間的關(guān)系。按照Chebyshev多項式零點抽樣法則,每個變量抽取4個樣本點,故在變量區(qū)間內(nèi)共抽取64(43=64)個樣本點,即采用有限元模型計算64組樣本數(shù)據(jù),并基于這些樣本數(shù)據(jù)構(gòu)建約束函數(shù)的Chebyshev多項式模型。對考慮確定性(DTO)的拓?fù)鋬?yōu)化設(shè)計和文中所提出的RTOD進(jìn)行了結(jié)果分析與對比。從圖6可以看出,與確定性拓?fù)鋬?yōu)化所得模型相比,RTOD的解具有不同的結(jié)構(gòu)拓?fù)湫螤?。毫無疑問,這種優(yōu)化設(shè)計的模型結(jié)構(gòu)是采用確定性變量進(jìn)行拓?fù)鋬?yōu)化設(shè)計所不能實現(xiàn)的。

    圖6 結(jié)構(gòu)確定性與不確定性拓?fù)鋬?yōu)化結(jié)果

    Fig.6 Results of structural deterministic and uncertain topology optimization

    表1 算例1中不確定參數(shù)在標(biāo)準(zhǔn)化空間的最優(yōu)值

    Table 1 Optimal values of uncertain parameters in the norma-lized space in example 1

    設(shè)計方法E?F?xF?yRBTO,γ=1.0-0.5497-0.5589-0.6208RBTO,γ=1.2-0.6692-0.6663-0.7404RBTO,γ=1.4-0.7917-0.7722-0.8584RBTO,γ=1.6-0.9171-0.8766-0.9749

    圖7顯示了不同可靠性指標(biāo)的迭代過程,由圖可見,RTOD所需的拓?fù)鋬?yōu)化迭代次數(shù)比基于確定性的拓?fù)鋬?yōu)化迭代次數(shù)更少。從數(shù)值算例中可以看出,該方法在求解整體拓?fù)浣Y(jié)構(gòu)優(yōu)化問題時將拓?fù)渑c可靠性分析解耦分離,通過將變量與目標(biāo)函數(shù)映射成顯式呈現(xiàn),再經(jīng)過可靠性分析,即可求得指定可靠性指標(biāo)所對應(yīng)的可靠性能值。

    表2 非概率可靠性拓?fù)鋬?yōu)化設(shè)計結(jié)果

    Table 2 Results of non-probabilistic reliability topology optimization design

    設(shè)計方法材料體積/mm3總迭代次數(shù)可靠性分析評估次數(shù)拓?fù)鋬?yōu)化迭代次數(shù)DTO5.666×102——34RTOD,γ=1.07.636×102401228RTOD,γ=1.27.842×10235728RTOD,γ=1.48.083×10235728RTOD,γ=1.68.557×10230723

    圖7 拓?fù)鋬?yōu)化迭代進(jìn)程Fig.7 Topology optimization iterative process

    4.2 算例2

    將應(yīng)變能約束下的可靠性應(yīng)用于如圖8所示的受變動載荷影響的三維框架結(jié)構(gòu)(算例2)。結(jié)構(gòu)左側(cè)端面固定,隨機(jī)載荷作用于右側(cè)端面的中心。采用單元大小為4 mm的有限元方法,設(shè)計域包含4 961個節(jié)點和4 000個單元。所施加的載荷為均勻分布的隨機(jī)變量,平均載荷F1為200 N,F(xiàn)2為

    圖8 三維框架設(shè)計區(qū)域Fig.8 3D frame design domain

    表3 算例2中不確定參數(shù)在標(biāo)準(zhǔn)化空間的最優(yōu)值

    Table 3 Example 2 optimal values of uncertain parameters in the normalized space

    設(shè)計方法E?F?1F?2F?3RBTO,γ=1.0-0.50830.84470.10650.1292RBTO,γ=1.2-0.63171.00000.12960.1553RBTO,γ=1.4-0.93451.00000.18750.2271

    表4 三維框架結(jié)構(gòu)優(yōu)化結(jié)果Table 4 3D frame structure optimization results

    圖9 結(jié)構(gòu)確定性與不確定性拓?fù)鋬?yōu)化結(jié)果分析

    Fig.9 Results analysis of structural deterministic and uncertain topology optimization

    圖10 原始密度灰度圖Fig.10 Original density grayscales

    在算例2中,內(nèi)優(yōu)化搜索在每個內(nèi)部循環(huán)中的迭代次數(shù)獨立于外部迭代次數(shù),需調(diào)用函數(shù)的次數(shù)非常小,表明了這種性能函數(shù)評估方法所固有的穩(wěn)健特性。圖11、12所示為整體結(jié)構(gòu)設(shè)計域優(yōu)化的迭代過程。整個優(yōu)化過程中,目標(biāo)函數(shù)不斷下降,經(jīng)過25次迭代后逐漸達(dá)到穩(wěn)定,結(jié)構(gòu)的優(yōu)化迭代趨勢基本一致。

    圖12 拓?fù)鋬?yōu)化應(yīng)變能迭代進(jìn)程Fig.12 Topology optimization iterative process of strain energy

    5 結(jié)語

    文中利用基于超橢球凸模型的可靠性優(yōu)化設(shè)計方法,結(jié)合區(qū)間Chebyshev多項式提出一種基于非概率模型的可靠性優(yōu)化設(shè)計方法。該方法將RTOD問題分離為獨立的可靠性分析與等價形式的確定性拓?fù)鋬?yōu)化,運算中的迭代次數(shù)相比確定性拓?fù)鋬?yōu)化的迭代次數(shù)更少。本研究為參數(shù)具有隨機(jī)性但有界的結(jié)構(gòu)體的可靠性拓?fù)湓O(shè)計提供了一種有效的方法,并驗證了材料和荷載的不確定性在應(yīng)變能約束問題的拓?fù)湓O(shè)計中起著重要的作用。文中算例分析表明,基于可靠性的拓?fù)鋬?yōu)化設(shè)計由于定量考慮了影響結(jié)構(gòu)性能的各種不確定因素,因此能更加有效地衡量結(jié)構(gòu)的安全性和經(jīng)濟(jì)性,使設(shè)計結(jié)果更加合理。

    在拓?fù)鋬?yōu)化過程中,載荷、邊界條件等因素可能會影響拓?fù)鋬?yōu)化的最終拓?fù)湫螤?,因此,后續(xù)研究中將深入探討拓?fù)鋬?yōu)化設(shè)計的最優(yōu)算法,并對同時存在隨機(jī)變量與區(qū)間變量的情形進(jìn)行討論。

    猜你喜歡
    橢球不確定性區(qū)間
    解兩類含參數(shù)的復(fù)合不等式有解與恒成立問題
    你學(xué)會“區(qū)間測速”了嗎
    法律的兩種不確定性
    法律方法(2022年2期)2022-10-20 06:41:56
    獨立坐標(biāo)系橢球變換與坐標(biāo)換算
    橢球槽宏程序編制及其Vericut仿真
    智能制造(2021年4期)2021-11-04 08:54:44
    英鎊或繼續(xù)面臨不確定性風(fēng)險
    中國外匯(2019年7期)2019-07-13 05:45:04
    橢球精加工軌跡及程序設(shè)計
    基于外定界橢球集員估計的純方位目標(biāo)跟蹤
    具有不可測動態(tài)不確定性非線性系統(tǒng)的控制
    區(qū)間對象族的可鎮(zhèn)定性分析
    色播在线永久视频| 丝袜美腿诱惑在线| 一级片免费观看大全| 老司机午夜十八禁免费视频| 免费观看精品视频网站| 亚洲午夜精品一区,二区,三区| 久久久久亚洲av毛片大全| www.999成人在线观看| 欧美成人性av电影在线观看| 无人区码免费观看不卡| 国产av一区在线观看免费| 精品高清国产在线一区| 亚洲专区国产一区二区| netflix在线观看网站| 国产一区二区三区综合在线观看| 人人妻人人爽人人添夜夜欢视频| а√天堂www在线а√下载| 真人一进一出gif抽搐免费| 欧美成人免费av一区二区三区| 亚洲第一欧美日韩一区二区三区| 成人三级做爰电影| 亚洲欧美日韩高清在线视频| 亚洲一区中文字幕在线| 欧美亚洲日本最大视频资源| 日本 av在线| 在线永久观看黄色视频| 久久天堂一区二区三区四区| 午夜亚洲福利在线播放| 高清毛片免费观看视频网站| 国产区一区二久久| 免费在线观看日本一区| 9热在线视频观看99| 亚洲电影在线观看av| 成人国产一区最新在线观看| av视频免费观看在线观看| 午夜福利,免费看| 淫秽高清视频在线观看| 欧美乱妇无乱码| 欧美色视频一区免费| 国产麻豆成人av免费视频| 国产av又大| av超薄肉色丝袜交足视频| 99riav亚洲国产免费| www.www免费av| 欧美乱妇无乱码| 热re99久久国产66热| 欧美精品亚洲一区二区| 91字幕亚洲| 身体一侧抽搐| 涩涩av久久男人的天堂| 91精品三级在线观看| 两个人视频免费观看高清| 在线观看免费日韩欧美大片| 大香蕉久久成人网| 后天国语完整版免费观看| 国产色视频综合| 首页视频小说图片口味搜索| 最新美女视频免费是黄的| 亚洲一区中文字幕在线| 12—13女人毛片做爰片一| 国产亚洲精品第一综合不卡| 日本欧美视频一区| 精品国产美女av久久久久小说| 中文字幕人妻丝袜一区二区| 少妇被粗大的猛进出69影院| 一进一出好大好爽视频| 欧美午夜高清在线| 热re99久久国产66热| 曰老女人黄片| 香蕉久久夜色| 国产成人av激情在线播放| 91精品三级在线观看| 亚洲熟妇中文字幕五十中出| 天堂动漫精品| 在线观看舔阴道视频| 免费在线观看黄色视频的| 这个男人来自地球电影免费观看| 69av精品久久久久久| 午夜精品久久久久久毛片777| 国产国语露脸激情在线看| 欧美国产日韩亚洲一区| 国产99久久九九免费精品| 亚洲自偷自拍图片 自拍| 老汉色∧v一级毛片| 老汉色∧v一级毛片| 波多野结衣av一区二区av| 99精品欧美一区二区三区四区| 午夜福利欧美成人| 国产亚洲精品第一综合不卡| 一区福利在线观看| 精品久久久久久久毛片微露脸| 久久狼人影院| 成人亚洲精品一区在线观看| 亚洲五月婷婷丁香| 高清毛片免费观看视频网站| 在线永久观看黄色视频| 国产成人精品久久二区二区免费| 69av精品久久久久久| 国产精品亚洲美女久久久| 一个人观看的视频www高清免费观看 | 免费高清在线观看日韩| 两个人视频免费观看高清| 香蕉国产在线看| 日韩欧美国产在线观看| 日韩精品中文字幕看吧| 婷婷精品国产亚洲av在线| 女性生殖器流出的白浆| 在线观看免费日韩欧美大片| 久久久久久久久免费视频了| 国产午夜福利久久久久久| 男女做爰动态图高潮gif福利片 | 看黄色毛片网站| 精品电影一区二区在线| 长腿黑丝高跟| 一进一出抽搐动态| 色综合婷婷激情| 国产高清视频在线播放一区| 国产三级黄色录像| 亚洲一卡2卡3卡4卡5卡精品中文| 丁香六月欧美| 乱人伦中国视频| 欧美成狂野欧美在线观看| 亚洲精品国产色婷婷电影| 搞女人的毛片| 脱女人内裤的视频| 国产亚洲精品av在线| 热99re8久久精品国产| 99国产精品免费福利视频| 男人舔女人下体高潮全视频| 极品教师在线免费播放| 国产欧美日韩一区二区三区在线| 少妇粗大呻吟视频| 亚洲中文字幕日韩| 亚洲第一欧美日韩一区二区三区| 亚洲成人精品中文字幕电影| 久久久国产欧美日韩av| 亚洲欧美一区二区三区黑人| 亚洲avbb在线观看| 亚洲人成网站在线播放欧美日韩| 亚洲久久久国产精品| 一区在线观看完整版| 亚洲精品av麻豆狂野| 一进一出抽搐gif免费好疼| 日韩三级视频一区二区三区| 极品人妻少妇av视频| 91九色精品人成在线观看| 人人妻,人人澡人人爽秒播| 欧美日韩乱码在线| 欧洲精品卡2卡3卡4卡5卡区| 日韩中文字幕欧美一区二区| 变态另类丝袜制服| 久久这里只有精品19| 咕卡用的链子| 啦啦啦免费观看视频1| av网站免费在线观看视频| 高清在线国产一区| 国产高清视频在线播放一区| 一夜夜www| 亚洲国产精品久久男人天堂| 999精品在线视频| 国产成人精品在线电影| 老司机福利观看| 国产视频一区二区在线看| 91精品三级在线观看| 国产成人av教育| 国产精品电影一区二区三区| 亚洲成人精品中文字幕电影| 高清毛片免费观看视频网站| 国产亚洲精品综合一区在线观看 | 一本综合久久免费| 亚洲激情在线av| 久9热在线精品视频| 久久影院123| www.熟女人妻精品国产| 色综合站精品国产| 国产欧美日韩一区二区三| 久久久久久久午夜电影| 国产单亲对白刺激| 一边摸一边做爽爽视频免费| 午夜福利免费观看在线| 精品日产1卡2卡| 精品不卡国产一区二区三区| АⅤ资源中文在线天堂| 桃色一区二区三区在线观看| 国产精品一区二区精品视频观看| 天堂√8在线中文| 免费在线观看亚洲国产| 亚洲五月色婷婷综合| av有码第一页| 欧美另类亚洲清纯唯美| 纯流量卡能插随身wifi吗| 免费观看人在逋| 色尼玛亚洲综合影院| 久久婷婷成人综合色麻豆| 国产一区二区三区在线臀色熟女| 亚洲国产精品久久男人天堂| 自线自在国产av| 欧美成人性av电影在线观看| 我的亚洲天堂| 亚洲国产精品sss在线观看| 免费看十八禁软件| 91九色精品人成在线观看| 少妇被粗大的猛进出69影院| 无人区码免费观看不卡| 手机成人av网站| 深夜精品福利| 91麻豆av在线| 亚洲国产欧美日韩在线播放| 日本a在线网址| 两个人看的免费小视频| 91成人精品电影| 午夜福利在线观看吧| 中文字幕av电影在线播放| 日韩成人在线观看一区二区三区| 亚洲五月婷婷丁香| 久久青草综合色| 亚洲熟女毛片儿| 午夜福利免费观看在线| 多毛熟女@视频| 悠悠久久av| 精品乱码久久久久久99久播| 十八禁人妻一区二区| 亚洲色图综合在线观看| 精品国产超薄肉色丝袜足j| АⅤ资源中文在线天堂| 日韩欧美在线二视频| 每晚都被弄得嗷嗷叫到高潮| 欧美在线一区亚洲| 一区二区三区国产精品乱码| av视频免费观看在线观看| av福利片在线| 变态另类成人亚洲欧美熟女 | 人人澡人人妻人| 日韩成人在线观看一区二区三区| 国产熟女午夜一区二区三区| 又黄又爽又免费观看的视频| 久久精品人人爽人人爽视色| 天堂影院成人在线观看| 亚洲自拍偷在线| 大香蕉久久成人网| 超碰成人久久| 日韩三级视频一区二区三区| 国产精品久久电影中文字幕| 此物有八面人人有两片| 亚洲成av人片免费观看| 亚洲午夜精品一区,二区,三区| 国产精品二区激情视频| 午夜精品久久久久久毛片777| 精品国内亚洲2022精品成人| www.自偷自拍.com| 老司机在亚洲福利影院| 久久婷婷成人综合色麻豆| 啦啦啦 在线观看视频| 国产欧美日韩综合在线一区二区| 国产精品1区2区在线观看.| 亚洲无线在线观看| 首页视频小说图片口味搜索| 在线十欧美十亚洲十日本专区| x7x7x7水蜜桃| 男人的好看免费观看在线视频 | 久久久久久人人人人人| 99精品欧美一区二区三区四区| 亚洲成a人片在线一区二区| 亚洲三区欧美一区| 啪啪无遮挡十八禁网站| 国产成人精品久久二区二区免费| 精品乱码久久久久久99久播| 国产激情久久老熟女| 成人国产一区最新在线观看| 黄片大片在线免费观看| 777久久人妻少妇嫩草av网站| 久久久精品欧美日韩精品| 国产精品一区二区在线不卡| 国产又爽黄色视频| 国内精品久久久久久久电影| 身体一侧抽搐| 国产精品久久视频播放| 欧美日韩瑟瑟在线播放| 亚洲九九香蕉| 看片在线看免费视频| 少妇被粗大的猛进出69影院| 老司机福利观看| 亚洲精品在线观看二区| 久9热在线精品视频| 久久国产精品影院| 国产成人精品久久二区二区91| 免费一级毛片在线播放高清视频 | 真人做人爱边吃奶动态| 亚洲男人天堂网一区| 欧美成人免费av一区二区三区| 97碰自拍视频| 一区在线观看完整版| 一级黄色大片毛片| 满18在线观看网站| 好看av亚洲va欧美ⅴa在| www.999成人在线观看| 国产激情久久老熟女| 久久人人精品亚洲av| 国产av又大| 欧美精品亚洲一区二区| 精品免费久久久久久久清纯| 9191精品国产免费久久| 国产一区二区在线av高清观看| 国产片内射在线| 美女免费视频网站| 欧美色欧美亚洲另类二区 | 久久精品国产99精品国产亚洲性色 | 久久婷婷人人爽人人干人人爱 | 亚洲第一青青草原| 午夜成年电影在线免费观看| 一a级毛片在线观看| 欧美黄色淫秽网站| 国产精品久久久久久亚洲av鲁大| 一区在线观看完整版| 亚洲欧美激情在线| 一级,二级,三级黄色视频| 欧美乱色亚洲激情| 成在线人永久免费视频| 国产亚洲精品综合一区在线观看 | 精品国产美女av久久久久小说| 日韩一卡2卡3卡4卡2021年| 日韩国内少妇激情av| 每晚都被弄得嗷嗷叫到高潮| 两性午夜刺激爽爽歪歪视频在线观看 | 村上凉子中文字幕在线| 日本免费a在线| 欧美黄色淫秽网站| 免费看a级黄色片| 日本一区二区免费在线视频| 丝袜在线中文字幕| 高清毛片免费观看视频网站| 国产亚洲精品久久久久久毛片| 亚洲av熟女| 91九色精品人成在线观看| 亚洲中文av在线| 我的亚洲天堂| 亚洲av成人av| 怎么达到女性高潮| 亚洲久久久国产精品| 制服人妻中文乱码| 午夜福利,免费看| 久久伊人香网站| av在线播放免费不卡| 欧美丝袜亚洲另类 | 国产av又大| 日本黄色视频三级网站网址| 成人手机av| 麻豆成人av在线观看| 日本 欧美在线| 午夜福利影视在线免费观看| 黄色视频,在线免费观看| 国产成人精品久久二区二区免费| 一个人免费在线观看的高清视频| 麻豆一二三区av精品| 啦啦啦 在线观看视频| 日韩大尺度精品在线看网址 | 看黄色毛片网站| 国产三级黄色录像| 精品久久久久久久人妻蜜臀av | 亚洲第一电影网av| 久久婷婷人人爽人人干人人爱 | 亚洲精品粉嫩美女一区| 欧美最黄视频在线播放免费| 久久久久久国产a免费观看| 高潮久久久久久久久久久不卡| 国产91精品成人一区二区三区| svipshipincom国产片| 国产又色又爽无遮挡免费看| 亚洲avbb在线观看| 丰满的人妻完整版| 亚洲情色 制服丝袜| 国产精品久久久久久精品电影 | 黄色成人免费大全| 免费搜索国产男女视频| 亚洲五月色婷婷综合| www.www免费av| 一边摸一边做爽爽视频免费| 欧美乱妇无乱码| 亚洲中文字幕一区二区三区有码在线看 | 国产精品亚洲美女久久久| 俄罗斯特黄特色一大片| 搡老熟女国产l中国老女人| 国产av精品麻豆| 欧美乱码精品一区二区三区| 免费看美女性在线毛片视频| 精品久久久精品久久久| 国产伦人伦偷精品视频| 精品福利观看| 又大又爽又粗| 麻豆久久精品国产亚洲av| 国产亚洲av嫩草精品影院| 不卡一级毛片| 在线观看免费午夜福利视频| 国产av在哪里看| 成人国语在线视频| 国产成人系列免费观看| 亚洲av电影不卡..在线观看| 亚洲视频免费观看视频| 一级毛片女人18水好多| 国产精品爽爽va在线观看网站 | 手机成人av网站| 电影成人av| 国产精品一区二区精品视频观看| 男人的好看免费观看在线视频 | 久久人妻福利社区极品人妻图片| 女性被躁到高潮视频| av视频在线观看入口| 色老头精品视频在线观看| 电影成人av| 在线观看免费视频日本深夜| 亚洲成av人片免费观看| 自线自在国产av| 18美女黄网站色大片免费观看| 亚洲国产日韩欧美精品在线观看 | 欧美日本中文国产一区发布| 中文字幕人妻熟女乱码| 无限看片的www在线观看| 啦啦啦 在线观看视频| 男男h啪啪无遮挡| 亚洲无线在线观看| 不卡av一区二区三区| 啪啪无遮挡十八禁网站| 动漫黄色视频在线观看| 久99久视频精品免费| 欧美+亚洲+日韩+国产| av片东京热男人的天堂| 亚洲第一av免费看| 天堂影院成人在线观看| 日韩三级视频一区二区三区| 香蕉久久夜色| cao死你这个sao货| 女同久久另类99精品国产91| 后天国语完整版免费观看| 午夜视频精品福利| 欧美丝袜亚洲另类 | 国产精品 欧美亚洲| 色综合亚洲欧美另类图片| 国产精品一区二区免费欧美| 午夜福利免费观看在线| www日本在线高清视频| 色精品久久人妻99蜜桃| 国产一级毛片七仙女欲春2 | 91九色精品人成在线观看| 亚洲欧美一区二区三区黑人| 亚洲五月色婷婷综合| 国产av又大| 神马国产精品三级电影在线观看 | 人人妻人人澡人人看| 午夜免费成人在线视频| 精品日产1卡2卡| 一进一出好大好爽视频| 丁香六月欧美| 日韩欧美国产一区二区入口| av片东京热男人的天堂| 久久精品aⅴ一区二区三区四区| av天堂久久9| 一二三四社区在线视频社区8| 在线国产一区二区在线| 两个人视频免费观看高清| 亚洲五月色婷婷综合| 久9热在线精品视频| 免费看a级黄色片| 亚洲男人的天堂狠狠| 人人妻人人爽人人添夜夜欢视频| 成人精品一区二区免费| 51午夜福利影视在线观看| 女同久久另类99精品国产91| 99久久国产精品久久久| 搡老岳熟女国产| 丰满人妻熟妇乱又伦精品不卡| 亚洲男人天堂网一区| 丝袜美足系列| 国产成人精品久久二区二区91| 中文亚洲av片在线观看爽| 操出白浆在线播放| 国产真人三级小视频在线观看| 国产片内射在线| 女生性感内裤真人,穿戴方法视频| 中文字幕精品免费在线观看视频| 日韩大码丰满熟妇| 亚洲 国产 在线| 国产单亲对白刺激| 欧美黄色片欧美黄色片| 在线观看舔阴道视频| 69av精品久久久久久| 国产午夜福利久久久久久| 老司机午夜十八禁免费视频| 久久婷婷人人爽人人干人人爱 | 久久人妻熟女aⅴ| 黄色成人免费大全| 国产麻豆成人av免费视频| 国产真人三级小视频在线观看| 久久人妻av系列| 免费高清在线观看日韩| 在线观看66精品国产| 在线观看免费午夜福利视频| 99re在线观看精品视频| 中文字幕久久专区| 91av网站免费观看| 欧美成人午夜精品| 操美女的视频在线观看| 两个人免费观看高清视频| 窝窝影院91人妻| 免费观看人在逋| 亚洲成人免费电影在线观看| 欧美色欧美亚洲另类二区 | e午夜精品久久久久久久| 欧美成人午夜精品| 妹子高潮喷水视频| 免费看十八禁软件| 欧美精品亚洲一区二区| 久久久久国内视频| 又紧又爽又黄一区二区| 欧美激情极品国产一区二区三区| 久久青草综合色| 精品久久久久久久毛片微露脸| 国产熟女xx| 色播亚洲综合网| 久久国产精品影院| 91国产中文字幕| 亚洲精品国产区一区二| 久久精品成人免费网站| 亚洲人成网站在线播放欧美日韩| 亚洲精华国产精华精| 精品福利观看| 色精品久久人妻99蜜桃| 久久 成人 亚洲| 亚洲男人天堂网一区| 人妻丰满熟妇av一区二区三区| 亚洲片人在线观看| 欧美黄色淫秽网站| 亚洲国产精品999在线| 性少妇av在线| 国产av在哪里看| 很黄的视频免费| 男女午夜视频在线观看| 国产精品永久免费网站| 天天添夜夜摸| 一区二区三区激情视频| av片东京热男人的天堂| 精品乱码久久久久久99久播| 黄片播放在线免费| av欧美777| 午夜精品在线福利| 欧美日本视频| 美女高潮喷水抽搐中文字幕| 欧美丝袜亚洲另类 | av天堂久久9| 99国产极品粉嫩在线观看| 51午夜福利影视在线观看| 淫妇啪啪啪对白视频| 国产亚洲精品第一综合不卡| 大香蕉久久成人网| 中文字幕精品免费在线观看视频| 美女高潮喷水抽搐中文字幕| 国产成人欧美在线观看| 老汉色∧v一级毛片| 最近最新中文字幕大全免费视频| 一级毛片精品| 国产不卡一卡二| 精品一区二区三区av网在线观看| 欧美av亚洲av综合av国产av| 国产欧美日韩综合在线一区二区| 免费看a级黄色片| 欧美久久黑人一区二区| 精品免费久久久久久久清纯| 成人av一区二区三区在线看| 成熟少妇高潮喷水视频| 91麻豆精品激情在线观看国产| 国产日韩一区二区三区精品不卡| 激情视频va一区二区三区| 丰满的人妻完整版| 国产片内射在线| 91成人精品电影| www国产在线视频色| 中文字幕人成人乱码亚洲影| 美女扒开内裤让男人捅视频| av有码第一页| av免费在线观看网站| 精品国产一区二区久久| 精品国产美女av久久久久小说| 国产伦一二天堂av在线观看| 国产熟女午夜一区二区三区| 精品人妻在线不人妻| 丝袜美足系列| 久久性视频一级片| 国产在线观看jvid| 99精品在免费线老司机午夜| 国产精品久久久久久人妻精品电影| 亚洲av成人不卡在线观看播放网| 国产日韩一区二区三区精品不卡| 久久精品国产亚洲av高清一级| 天天躁夜夜躁狠狠躁躁| 欧美大码av| 欧美成人性av电影在线观看| 纯流量卡能插随身wifi吗| www.熟女人妻精品国产| 黄色视频,在线免费观看| 亚洲精品国产区一区二| 欧美中文综合在线视频| 99国产综合亚洲精品| 麻豆久久精品国产亚洲av| svipshipincom国产片| 香蕉丝袜av| ponron亚洲| 国产av一区在线观看免费| 丝袜美腿诱惑在线| 久久久久久人人人人人| 天堂影院成人在线观看| 亚洲七黄色美女视频| 亚洲av成人不卡在线观看播放网| 国产在线精品亚洲第一网站| 亚洲黑人精品在线| 涩涩av久久男人的天堂| 高潮久久久久久久久久久不卡|