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

    混合B樣條實體模型的等幾何拓?fù)鋬?yōu)化

    2021-07-12 01:16:28楊佳明郭馬一杜孝孝
    圖學(xué)學(xué)報 2021年3期
    關(guān)鍵詞:控制頂點實體模型樣條

    楊佳明,趙 罡,2,王 偉,3,郭馬一,杜孝孝

    混合B樣條實體模型的等幾何拓?fù)鋬?yōu)化

    楊佳明1,趙 罡1,2,王 偉1,3,郭馬一1,杜孝孝1

    (1. 北京航空航天大學(xué)機械工程及自動化學(xué)院,北京 100191;2. 航空高端裝備智能制造技術(shù)工業(yè)和信息化部重點實驗室,北京 100191;3. 北京市高效綠色數(shù)控加工工藝及裝備工程技術(shù)研究中心,北京 100191)

    等幾何拓?fù)鋬?yōu)化方法將經(jīng)典拓?fù)鋬?yōu)化理論中的有限元分析過程更改為等幾何分析計算,從而提高了拓?fù)鋬?yōu)化的效率與穩(wěn)定性。針對現(xiàn)有的等幾何拓?fù)鋬?yōu)化方法在處理復(fù)雜實體結(jié)構(gòu)優(yōu)化問題時具有一定的局限性,提出一種非結(jié)構(gòu)化樣條實體等幾何拓?fù)鋬?yōu)化方法。基于混合B樣條構(gòu)造技術(shù),在非結(jié)構(gòu)化六面體網(wǎng)格上構(gòu)造具有復(fù)雜結(jié)構(gòu)的樣條實體,并將其作為拓?fù)鋬?yōu)化問題的設(shè)計域。用于描述這一樣條實體的基函數(shù)被直接應(yīng)用于材料密度分布的表達(dá)以及等幾何分析計算。在數(shù)值算例中,該方法表現(xiàn)出應(yīng)用于復(fù)雜結(jié)構(gòu)時的良好穩(wěn)定性和魯棒性。研究成果對等幾何拓?fù)鋬?yōu)化方法應(yīng)用于實際工程問題具有一定的參考意義。

    拓?fù)鋬?yōu)化;等幾何分析;體參數(shù)化;B樣條;非結(jié)構(gòu)化樣條

    拓?fù)鋬?yōu)化是一種重要的工程結(jié)構(gòu)優(yōu)化方法。目前常見的拓?fù)鋬?yōu)化方法包括:均勻化方法[1]、固體各向同性懲罰微結(jié)構(gòu)模型(solid isotropic microstructures with penalization,SIMP)法[2-3]、漸進結(jié)構(gòu)優(yōu)化方法[4]、水平集方法[5-7]、移動組件法[8-9]等。由于這些方法可以在設(shè)計域中按照給定約束條件高效地搜索最優(yōu)的設(shè)計結(jié)構(gòu),因此被廣泛應(yīng)用于各類工程問題。其中,AAGE等[10]將千兆體素級的SIMP法應(yīng)用于全尺寸機翼的結(jié)構(gòu)設(shè)計并將該成果發(fā)表于《Nature》期刊上,體現(xiàn)了拓?fù)鋬?yōu)化技術(shù)的應(yīng)用潛力。

    這些經(jīng)典的拓?fù)鋬?yōu)化方法一般都會通過有限元方法來實現(xiàn)結(jié)構(gòu)響應(yīng)的計算。但傳統(tǒng)的有限元方法存在以下問題和局限性:①低階連續(xù)性的有限元有可能影響優(yōu)化計算的準(zhǔn)確性并導(dǎo)致數(shù)值計算的不穩(wěn)定 性[11];②繁瑣耗時的有限元前處理過程割裂了幾何設(shè)計與力學(xué)分析過程,從而使拓?fù)鋬?yōu)化難以作為一種設(shè)計工具完全融入到現(xiàn)有的計算機輔助設(shè)計系統(tǒng)之中。

    HUGHES等[12-13]推廣的等幾何分析方法作為有限元方法,被認(rèn)為有可能突破上述局限性。其核心思想是直接利用CAD中描述幾何模型的樣條基函數(shù)來表達(dá)力學(xué)分析中的未知場。由于該方法統(tǒng)一了設(shè)計與分析過程中模型的基礎(chǔ)表達(dá)形式,避免了模型轉(zhuǎn)換引起的誤差,因此具有諸多優(yōu)勢。國內(nèi)外很多學(xué)者將其應(yīng)用于拓?fù)鋬?yōu)化,從而形成了等幾何拓?fù)鋬?yōu)化方法。其中,SEO等[14]最先提出在拓?fù)鋬?yōu)化中使用等幾何分析方法及裁剪樣條曲面;隨后,HASSANI等[15]提出了在點密度SIMP法基礎(chǔ)上引入NURBS基函數(shù)的等幾何拓?fù)鋬?yōu)化方法;QIAN[16]將設(shè)計域嵌入到一個B樣條空間中,證明了基于樣條函數(shù)的材料密度分布表示方法不僅具有良好的效率和魯棒性,并且具有固有的過濾性質(zhì)。國內(nèi)WANG和BENSON[17]將等幾何分析引入到水平集方法中,并針對等幾何拓?fù)鋬?yōu)化方法提出了一種高效的計算框架[18];XIE等[19-20]探索了等幾何拓?fù)鋬?yōu)化在移動組件法中的應(yīng)用;GAO等[21-22]對基于NURBS的等幾何拓?fù)鋬?yōu)化方法及相關(guān)應(yīng)用開展了深入的研究。

    可以看到,現(xiàn)有的等幾何拓?fù)鋬?yōu)化研究大多基于經(jīng)典的樣條構(gòu)造方法,如B樣條、NURBS。而這些張量積樣條方法受制于矩形的拓?fù)浣Y(jié)構(gòu)。盡管可以使用裁剪、拼接等造型方法進行處理,但要保證復(fù)雜設(shè)計域模型的水密性依舊是一件極具挑戰(zhàn)的任務(wù)。現(xiàn)有的等幾何拓?fù)鋬?yōu)化方法大多只能在簡單的矩形設(shè)計域上進行驗證,制約了在復(fù)雜問題上的探索。尤其對于三維實體問題,一些在復(fù)雜平面設(shè)計域問題中表現(xiàn)良好的方法難以直接推廣至實體層面,如ZHAO等[23]提出的基于非結(jié)構(gòu)化T樣條曲面的等幾何拓?fù)鋬?yōu)化方法。

    針對上述問題,本文將一種可以在非結(jié)構(gòu)化體網(wǎng)格上構(gòu)造樣條實體的方法引入等幾何拓?fù)鋬?yōu)化。由于允許奇異邊(點)的存在,因此可以對高虧格的復(fù)雜結(jié)構(gòu)進行表達(dá),突破傳統(tǒng)方法的矩形拓?fù)湎拗疲贤負(fù)鋬?yōu)化過程的需求。本文使用WEI等[24]提出的混合B樣條構(gòu)造方法。此方法構(gòu)造的樣條函數(shù)不僅具有優(yōu)良的性質(zhì),如非負(fù)性、規(guī)范性、線性無關(guān)性,并且構(gòu)造過程可靠,具有很好的適應(yīng)性。本文基于此方法構(gòu)造非結(jié)構(gòu)化的樣條實體模型,并給出了在此類模型上實現(xiàn)等幾何拓?fù)鋬?yōu)化的方法,通過算例驗證了其可行性,及對復(fù)雜實體結(jié)構(gòu)的處理能力。該研究成果有望推動等幾何拓?fù)鋬?yōu)化方法在實際工程問題中的應(yīng)用與推廣。

    1 混合B樣條實體構(gòu)造方法

    本節(jié)簡要回顧B樣條及非結(jié)構(gòu)化六面體網(wǎng)格上的混合B樣條構(gòu)造方法,相關(guān)細(xì)節(jié)可以參考文獻(xiàn)[24-25]。

    1.1 B樣條

    雙變量/三變量B樣條基函數(shù)可以通過單變量B樣條基函數(shù)的張量積得到。

    1.2 混合B樣條構(gòu)造方法

    其中,為中的一部分控制頂點。

    由節(jié)點插入算法[27]可知,,和0之間存在線性變換關(guān)系

    其中,可以直接通過節(jié)點插入算法得到。將其代入式(4)可以得到基函數(shù)之間的線性關(guān)系

    圖2 三變量B樣條基函數(shù)的截斷((a)一個C2控制頂點;(b)基函數(shù)的等值面;(c) C2控制點周圍存在一個激活的C0控制頂點;(d)從中截取后得到的截斷基函數(shù)的等值面)

    文獻(xiàn)[24]將上述情況稱為02構(gòu)造。除此之外,根據(jù)激活的0,1控制頂點不同還存在012和02構(gòu)造形式。相比于012和02,本文采用的02方法盡管在奇異點鄰域的連續(xù)性上存在不足,但具有構(gòu)造形式簡單以及基函數(shù)局部線性無關(guān)等優(yōu)勢。

    1.3 C02混合B樣條實體構(gòu)造方法

    一個非結(jié)構(gòu)化六面體網(wǎng)格中包含了頂點、邊、面和體單元4類元素。定義邊(頂點)的價為:包含這條邊(頂點)的體單元的個數(shù)。定義奇異邊為:價不為4的內(nèi)部邊以及價不為1或2的邊界邊。定義奇異點為:奇異邊的端點。

    如果一個體單元的頂點中包含了奇異點,則稱這個體單元為不規(guī)則單元。其余的體單元稱為規(guī)則單元。規(guī)則單元又可分為2類:1-鄰域內(nèi)存在不規(guī)則單元的規(guī)則單元稱為規(guī)則過渡單元,否則稱為規(guī)則非過渡單元。圖3為一個包含三價奇異邊的簡單非結(jié)構(gòu)化體網(wǎng)格,并通過顏色區(qū)分了不同的單元類型。

    圖3 非結(jié)構(gòu)化體網(wǎng)格及3種單元類型

    1.3.1 不規(guī)則單元

    其中,0包含了64個三-三次貝奇爾基函數(shù)。

    1.3.2 規(guī)則非過渡單元

    1.3.3 規(guī)則過渡單元

    圖4 eRT與相鄰非規(guī)則單元之間不同的相對位置關(guān)系

    1.3.4 統(tǒng)一構(gòu)造形式

    綜上,由式(13)、式(15)和式(17),可將非結(jié)構(gòu)化六面體網(wǎng)格上的02混合B樣條實體構(gòu)造總結(jié)為

    可以發(fā)現(xiàn),3類單元均可寫為統(tǒng)一的形式,即

    1.3.5 連續(xù)性

    按照式(19)構(gòu)造的樣條實體模型,在規(guī)則單元間為2連續(xù),在非規(guī)則單元間為0連續(xù),在非規(guī)則單元和規(guī)則單元間為0連續(xù)。02混合B樣條實體的連續(xù)性分析的簡要證明如下:

    為了討論過渡單元與其相鄰規(guī)則單元間的連續(xù)性,引入以下結(jié)論:

    引理[24]1. 對于任意一個固定的單元,其混合B樣條表達(dá)形式(式(12))等價于這個單元的B樣條表達(dá)形式(式(4))。

    為了使模型插值于邊界,本文將邊界單元均按照非規(guī)則單元處理。因此,此時邊界曲面為0連續(xù)曲面。為了提高邊界曲面的連續(xù)性,可使用非結(jié)構(gòu)化四邊形網(wǎng)格的012混合構(gòu)造方法對其進行處理。當(dāng)存在原始的B-rep模型時,也可通過調(diào)整邊界控制頂點使其恢復(fù)原始曲面的光順性。

    2 等幾何拓?fù)鋬?yōu)化

    本節(jié)對混合B樣條實體的等幾何拓?fù)鋬?yōu)化方法進行介紹。為了簡化問題的討論同時又不失一般性,本文著重考慮拓?fù)鋬?yōu)化問題中的最小化柔順度問題。

    2.1 基于混合B樣條基函數(shù)的材料密度分布

    其中,0為材料本身的楊氏模量;min為設(shè)定的最小剛度值,用以避免剛度矩陣奇異的情況;為懲罰因子。

    2.2 等幾何分析

    在等幾何方法中,可以將靜態(tài)載荷下線彈性問題的剛度方程表示為

    其中,為幾何映射(e)的雅可比矩陣。同時,可以通過式(23)對物理域坐標(biāo)的偏導(dǎo)求得;則通過式(22)表示為

    其中,0為楊氏模量等于1時的彈性系數(shù)矩陣。

    使用高斯積分法來完成式(26)的數(shù)值積分運算,則可表示為

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

    2.4 靈敏度分析

    則目標(biāo)函數(shù)的靈敏度,即式(31)可寫為

    2.5 優(yōu)化流程

    圖5 混合B樣條實體等幾何拓?fù)鋬?yōu)化流程圖

    3 數(shù)值算例

    3.1 算例1

    圖6(b)為用于構(gòu)造此設(shè)計域樣條實體模型的非結(jié)構(gòu)化六面體網(wǎng)格,通過彩色點標(biāo)出了模型表面的部分奇異點。可以看到本文所使用的構(gòu)造方法可以支持邊界奇異點,以及奇異點相鄰的情況。圖6(c)和6(d)分別給出了本文構(gòu)造的02混合B樣條實體及其控制頂點分布。此模型的控制頂點數(shù)(同時也為設(shè)計變量數(shù))為17 745,單元數(shù)為540??梢钥吹侥P蛢?nèi)部包含了一些奇異邊(點)。當(dāng)直接使用傳統(tǒng)張量積樣條方法進行構(gòu)造時,需要考慮多塊拼接、裁剪等問題。而本文使用的02混合構(gòu)造方法在處理時則較為簡潔。

    圖6 算例1的設(shè)計域及混合B樣條實體模型((a)設(shè)計域;(b)用于構(gòu)造混合B樣條實體的非結(jié)構(gòu)化六面體網(wǎng)格;(c)構(gòu)造的混合B樣條實體模型及等幾何分析結(jié)果(x向位移);(d)樣條實體模型的控制點)

    等幾何拓?fù)鋬?yōu)化結(jié)果如圖7(a)所示,最終目標(biāo)函數(shù)值為60.283 8,迭代次數(shù)為78。作為對比,圖7(b)中展示了SIMP方法的優(yōu)化結(jié)果。該模型的單元數(shù)為22 680,最終目標(biāo)函數(shù)值為66.488 1。相比于SIMP方法,本文方法以較少的單元數(shù)得到了較好的優(yōu)化結(jié)果(更低的柔順度)。主要是由于在SIMP方法中,材料密度分布通過離散單元的形式來表達(dá),每一個單元對應(yīng)一個單一的密度值。而本文方法以連續(xù)的樣條函數(shù)來表達(dá)材料密度分布,每個單元都具有內(nèi)部的密度分布函數(shù),且單元間密度分布連續(xù),因此可以通過較少的單元表達(dá)復(fù)雜的材料密度分布結(jié)構(gòu)。

    圖7 算例1的拓?fù)鋬?yōu)化結(jié)果((a)混合B樣條實體等幾何拓?fù)鋬?yōu)化結(jié)果;(b) SIMP法優(yōu)化結(jié)果)

    圖8為本例中等幾何拓?fù)鋬?yōu)化的收斂曲線以及部分中間結(jié)果??梢钥吹秸麄€優(yōu)化過程的收斂速度較快,且較為光滑,沒有出現(xiàn)明顯的振蕩。說明本文的方法具有較好的收斂性質(zhì)。

    圖8 算例1的優(yōu)化收斂曲線

    3.2 算例2

    本例為一個復(fù)雜實體結(jié)構(gòu)的拓?fù)鋬?yōu)化問題。設(shè)計域及邊界條件如圖9(a)所示,圖中紅色圓柱面固定,3個角上的藍(lán)色平面施加垂直于平面向外的均布載荷,大小為1。本文將體積分?jǐn)?shù)設(shè)為40%,設(shè)計變量的初始值為0.4。

    圖9(b)為用于構(gòu)造此設(shè)計域樣條實體模型的非結(jié)構(gòu)化六面體網(wǎng)格。在圖9(c)中展示了在樣條實體構(gòu)造過程中識別的不同單元類型。其中,紅色為非規(guī)則單元,白色為規(guī)則單元。規(guī)則單元間具有2連續(xù)性。圖9(d)為最終構(gòu)造的02混合B樣條實體及等幾何分析結(jié)果??梢钥吹轿灰茍鲈谀P蛢?nèi)部具有較高的連續(xù)性。該模型共有162 303個控制頂點,如圖10(a)所示。如果完全采用貝奇爾提取方法來處理此樣條實體模型,則控制頂點數(shù)為197 895,如圖10(b)所示。由于02構(gòu)造方法在模型內(nèi)部具有更高的連續(xù)性,因此使用了更少的控制頂點,其分布更為稀疏,同時減少了設(shè)計變量的規(guī)模。

    圖9 算例2的設(shè)計域及混合B樣條實體模型((a)設(shè)計域;(b)用于構(gòu)造混合B樣條實體的非結(jié)構(gòu)化六面體網(wǎng)格;(c)不同的單元類型:紅色為非規(guī)則單元,白色為規(guī)則單元;(d)構(gòu)造的混合B樣條實體模型及等幾何分析結(jié)果(x向位移))

    圖11為本例中等幾何拓?fù)鋬?yōu)化的結(jié)果,最終目標(biāo)函數(shù)值為112.159,迭代次數(shù)為100。圖12為目標(biāo)函數(shù)收斂曲線、體積分?jǐn)?shù)收斂曲線,以及部分中間結(jié)果。可以看到,同算例1,在本例中優(yōu)化收斂速度較快,無明顯振蕩,說明本文方法在處理復(fù)雜結(jié)構(gòu)時也具有較好的穩(wěn)定性和魯棒性。

    圖10 混合B樣條實體與完全采用貝奇爾提取方法時的控制頂點分布對比,圖中每一個紅色點表示一個控制頂點,紅色點越密集的區(qū)域控制頂點數(shù)目越多((a)混合B樣條實體的控制頂點分布;(b)完全采用貝奇爾提取方法時的控制頂點分布)

    圖11 算例2的等幾何拓?fù)鋬?yōu)化結(jié)果

    圖12 算例2的優(yōu)化收斂曲線

    3.3 算例3

    本例為一個多工況復(fù)雜實體結(jié)構(gòu)的拓?fù)鋬?yōu)化問題。設(shè)計域及2個工況狀態(tài)如圖13(a)和圖13(b)所示。圖中4個紅色圓柱面固定。在工況1時兩側(cè)懸臂的藍(lán)色圓柱面區(qū)域分別受到軸正向和負(fù)向的載荷,大小為1。在工況2時兩側(cè)懸臂的藍(lán)色圓柱面區(qū)域分別受到軸正向和負(fù)向的載荷,大小為1。取體積分?jǐn)?shù)約束為30%。

    按照本文方法所構(gòu)造的樣條實體如圖13(c)所示,圖中同時顯示了工況1下的等幾何分析結(jié)果。此模型的控制頂點分布如圖13(d)所示,總個數(shù)為110 283。可以看到由于此模型包含了多個孔洞結(jié)構(gòu),因此如果直接使用張量積NUBRS實體進行構(gòu)造將會遇到很大困難。

    圖13 算例3的設(shè)計域及混合B樣條實體模型((a)設(shè)計域及工況1;(b)設(shè)計域及工況2;(c)構(gòu)造的混合B樣條實體模型及等幾何分析結(jié)果(z向位移);(d)樣條實體模型的控制點)

    在優(yōu)化過程中本文設(shè)定目標(biāo)函數(shù)為2種工況下柔順度的平均值。等幾何拓?fù)鋬?yōu)化結(jié)果如圖14所示,最終目標(biāo)函數(shù)值為141.173,其中工況1的柔順度為267.901,工況2為14.445,迭代次數(shù)為100次。此算例進一步證明了該方法在復(fù)雜實體結(jié)構(gòu)拓?fù)鋬?yōu)化問題中具有良好的適應(yīng)性。

    圖14 算例3的等幾何拓?fù)鋬?yōu)化結(jié)果

    4 結(jié)束語

    本文針對現(xiàn)有等幾何拓?fù)鋬?yōu)化方法難以處理復(fù)雜實體結(jié)構(gòu)設(shè)計域這一問題,提出了一種基于混合B樣條構(gòu)造的非結(jié)構(gòu)化樣條實體等幾何拓?fù)鋬?yōu)化方法。此方法可在非結(jié)構(gòu)化六面體網(wǎng)格上構(gòu)造具有一定內(nèi)部連續(xù)性的樣條實體。本文將構(gòu)造得到的樣條基函數(shù)直接用于描述拓?fù)鋬?yōu)化問題中的材料密度分布,以及結(jié)構(gòu)響應(yīng)計算中的未知場。給出了等幾何拓?fù)鋬?yōu)化表達(dá)式以及相應(yīng)的靈敏度計算公式。最后,通過3個算例驗證了本文方法在處理復(fù)雜結(jié)構(gòu)優(yōu)化問題時的有效性和穩(wěn)定性。在這一研究成果基礎(chǔ)之上,有望形成面向?qū)嶋H工程中復(fù)雜模型的等幾何拓?fù)鋬?yōu)化設(shè)計工具。在后續(xù)研究中,如何提高奇異點鄰域的連續(xù)性將會是一個關(guān)鍵的研究方向。

    [1] BENDS?E M P, Kikuchi N. Generating optimal topologies in structural design using a homogenization method[J]. Computer Methods in Applied Mechanics and Engineering, 1988, 71(2): 197-224.

    [2] Zhou M, Rozvany G I N. The COC algorithm, part II: topological, geometry and generalized shape optimization[J]. Computer Methods in Applied Mechanics and Engineering, 1991, 89(1-3): 309-336.

    [3] Bends?e M P, Sigmund O. Material interpolation schemes in topology optimization[J]. Archive of Applied Mechanics, 1999, 69(9-10): 635-654.

    [4] Xie Y M, Steven G P. A simple evolutionary procedure for structural optimization[J]. Computers & Structures, 1993, 49(5): 885-896.

    [5] Sethian J A, Wiegmann A. Structural boundary design via level set and immersed interface methods[J]. Journal of Computational Physics, 2000, 163(2): 489-528.

    [6] Wang M Y, Wang X, Guo D. A level set method for structural topology optimization[J]. Computer Methods in Applied Mechanics and Engineering, 2003, 192(1-2): 227-246.

    [7] Allaire G, Jouve F, Toader A M. Structural optimization using sensitivity analysis and a level-set method[J]. Journal of Computational Physics, 2004, 194(1): 363-393.

    [8] Guo X, Zhang W S, Zhong W L. Doing topology optimization explicitly and geometrically - a new moving morphable components based framework[J]. Journal of Applied Mechanics, 2014, 81(8): 081009.

    [9] Zhang W S, Song J F, Zhou J H, et al. Topology optimization with multiple materials via Moving Morphable Component (MMC) method[J]. International Journal for Numerical Methods in Engineering, 2017, 113(11): 1653-1675.

    [10] Aage N, Andreassen E, Lazarov B S, et al. Giga-voxel computational morphogenesis for structural design[J]. Nature, 2017, 550: 84-86.

    [11] Sigmund O, Petersson J. Numerical instabilities in topology optimization: a survey on procedures dealing with checkerboards, mesh-dependencies and local minima[J]. Structural and Multidisciplinary Optimization, 1998, 16(1): 68-75.

    [12] Hughes T J R, Cottrell J A, Bazilevs Y. Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement[J]. Computer Methods in Applied Mechanics and Engineering, 2005, 194(39-41): 4135-4195.

    [13] Cottrell J A, Hughes T J R, Bazilevs Y. Isogeometric analysis: toward Integration of CAD and FEA[M]. Chichester: John Wiley & Sons, 2009: 1-18.

    [14] Seo Y D, Kim H J, Youn S K. Isogeometric topology optimization using trimmed spline surfaces[J]. Computer Methods in Applied Mechanics and Engineering, 2010, 199(49-52): 3270-3296.

    [15] Hassani B, Khanzadi M, Tavakkoli S M. An isogeometrical approach to structural topology optimization by optimality criteria[J]. Structural and Multidisciplinary Optimization, 2012, 45(2): 223-233.

    [16] Qian X P. Topology optimization in B-spline space[J]. Computer Methods in Applied Mechanics and Engineering, 2013, 265(1): 15-35.

    [17] Wang Y J, Benson D J. Isogeometric analysis for parameterized LSM-based structural topology optimization[J]. Computational Mechanics, 2016, 57(1): 19-35.

    [18] Wang Y, Liao Z Y, Ye M, et al. An efficient isogeometric topology optimization using multilevel mesh, MGCG and local-update strategy[J]. Advances in Engineering Software, 2020, 139: 102733.

    [19] Xie X D, WANG S T, Xu M M, et al. A new isogeometric topology optimization using moving morphable components based on R-functions and collocation schemes[J]. Computer Methods in Applied Mechanics and Engineering, 2018, 339(1): 61-90.

    [20] Gai Y D, Zhu X F, Zhang Y J, et al. Explicit isogeometric topology optimization based on moving morphable voids with closed B-spline boundary curves[J]. Structural and Multidisciplinary Optimization, 2019, 61: 963-982.

    [21] Gao J, Luo Z, Xiao M, et al. A NURBS-based Multi-Material Interpolation (N-MMI) for isogeometric topology optimization of structures[J]. Applied Mathematical Modelling, 2020, 81: 818-843.

    [22] Gao J, Xue H P, Gao L, et al. Topology optimization for auxetic metamaterials based on isogeometric analysis[J]. Computer Methods in Applied Mechanics and Engineering, 2019, 352(1): 211-236.

    [23] Zhao G, Yang J M, Wang W, et al. T-splines based isogeometric topology optimization with arbitrarily shaped design domains[J]. Computer Modeling in Engineering and Sciences, 2020, 123(3): 1033-1059.

    [24] Wei X D, Zhang Y J, Toshniwal D, et al. Blended B-spline construction on unstructured quadrilateral and hexahedral meshes with optimal convergence rates in isogeometric analysis[J]. Computer Methods in Applied Mechanics and Engineering, 2018, 341(1): 609-639.

    [25] Wei X D, Zhang Y J, Hughes T J R. Truncated hierarchical tricubic C0 spline construction on unstructured hexahedral meshes for isogeometric analysis applications[J]. Computers & Mathematics with Applications, 2017, 74(9): 2203-2220.

    [26] Giannelli C, Jüttler B, Speleers H. THB-splines: the truncated basis for hierarchical splines[J]. Computer Aided Geometric Design, 2012, 29(7): 485-498.

    [27] Boehm W. Inserting new knots into B-spline curves[J]. Computer Aided Design, 1980, 12(4): 199-201.

    [28] Scott M A, Simpson R N, Evans J A, et al. Isogeometric boundary element analysis using unstructured T-splines[J]. Computer Methods in Applied Mechanics and Engineering, 2013, 254: 197-221.

    [29] Bends?e M P, Sigmund O. Topology optimization: theory, methods, and applications[M]. Heidelberg: Springer, 2004: 9-28.

    [30] Svanberg K. The method of moving asymptotes - a new method for structural optimization[J]. International Journal for Numerical Methods in Engineering, 1987, 24(2): 359-373.

    Isogeometric topology optimization of blended B-spline solid model

    YANG Jia-ming1, ZHAO Gang1,2, WANG Wei1,3, GUO Ma-yi1, DU Xiao-xiao1

    (1. School of Mechanical Engineering & Automation, Beihang University, Beijing 100191, China;2. Key Laboratory of Aeronautics Smart Manufacturing, Ministry of Industry and Information Technology, Beijing 100191, China;3. Beijing Engineering Technological Research Center of High-Efficient & Green CNC Machining Process and Equipment, Beijing 100191, China)

    For isogeometric topology optimization (ITO) methods, isogeometric analysis (IGA) is adopted for topology optimization to address the limitation of the finite element method, which can improve the efficiency and stability of the optimization. However, it is of great challenge for existing ITO methods to manage arbitrarily shaped design domains, especially in three-dimensional solid problems. Therefore, a new ITO method was proposed to handle unstructured solid models. A spline solid with complex structures was obtained from an unstructured hexahedral mesh based on the blended B-spline construction. The basis functions describing the unstructured spline solid were applied to the representation of density distribution and the calculation of IGA. Several examples proved the flexibility and robustness of the proposed method in dealing with complex structures. These results may shed light on the application of ITO in practical engineering problems.

    topology optimization; isogeometric analysis; volume parameterization; B-spline; unstructured splines

    TP 391

    10.11996/JG.j.2095-302X.2021030501

    A

    2095-302X(2021)03-0501-10

    2020-09-17;

    2020-10-15

    17 September,2020;

    15 October,2020

    國家自然科學(xué)基金項目(61972011,61572056)

    National Natural Science Foundation of China (61972011, 61572056)

    楊佳明(1995-),男,北京人,博士研究生。主要研究方向為等幾何拓?fù)鋬?yōu)化、T樣條。E-mail:yangjiaming@buaa.edu.cn

    YANG Jia-ming (1995-), male, PhD candidate. His main research interests cover isogeometric topology optimization and T-splines. E-mail: yangjiaming@buaa.edu.cn

    王 偉(1978-),男,河北衡水人,副教授,博士。主要研究方向為CAD/CAE、智能計算等。E-mail:jrrt@buaa.edu.cn

    WANG Wei (1978-), male, associate professor, Ph.D. His main research interests cover CAD/CAE, intelligence computation, etc. E-mail:jrrt@buaa.edu.cn

    猜你喜歡
    控制頂點實體模型樣條
    帶互異權(quán)值的B樣條曲線的最小二乘漸進迭代逼近
    一元五次B樣條擬插值研究
    三次參數(shù)樣條在機床高速高精加工中的應(yīng)用
    三次樣條和二次刪除相輔助的WASD神經(jīng)網(wǎng)絡(luò)與日本人口預(yù)測
    軟件(2017年6期)2017-09-23 20:56:27
    基于樣條函數(shù)的高精度電子秤設(shè)計
    建筑設(shè)計中的實體模型
    建筑設(shè)計中的實體模型
    有理二次Bézier形式共軛雙曲線段的幾何計算
    面向控制頂點優(yōu)化的自由曲線交互擬合技術(shù)
    曲率對車身A級曲面控制頂點排列的影響*
    又大又黄又爽视频免费| 赤兔流量卡办理| 国产成人免费观看mmmm| 久久鲁丝午夜福利片| 韩国av在线不卡| 男女高潮啪啪啪动态图| 欧美黄色片欧美黄色片| 成人18禁高潮啪啪吃奶动态图| 国产成人精品久久久久久| 老女人水多毛片| av不卡在线播放| 中文精品一卡2卡3卡4更新| 满18在线观看网站| 考比视频在线观看| 999精品在线视频| 精品第一国产精品| 少妇人妻久久综合中文| 美女大奶头黄色视频| 国产av码专区亚洲av| 一二三四中文在线观看免费高清| 中文字幕人妻熟女乱码| 人成视频在线观看免费观看| 久久久久久久久久久免费av| 久久久久国产网址| 亚洲欧美成人精品一区二区| 少妇熟女欧美另类| 午夜老司机福利剧场| 在线观看三级黄色| 如日韩欧美国产精品一区二区三区| 久久精品国产亚洲av高清一级| 精品少妇黑人巨大在线播放| 欧美成人精品欧美一级黄| 国产精品嫩草影院av在线观看| 99久久综合免费| 青青草视频在线视频观看| 亚洲欧洲日产国产| 亚洲一级一片aⅴ在线观看| 精品国产乱码久久久久久小说| 精品福利永久在线观看| 最近中文字幕2019免费版| 国产白丝娇喘喷水9色精品| 国产片特级美女逼逼视频| 黄色毛片三级朝国网站| 麻豆av在线久日| 婷婷成人精品国产| 日韩一区二区三区影片| 在线免费观看不下载黄p国产| 日韩伦理黄色片| 1024香蕉在线观看| 欧美精品亚洲一区二区| 又大又黄又爽视频免费| 亚洲av欧美aⅴ国产| 精品视频人人做人人爽| 丝袜在线中文字幕| 国产精品人妻久久久影院| 亚洲欧美清纯卡通| 国产精品无大码| 国产免费现黄频在线看| 人妻一区二区av| 丝瓜视频免费看黄片| 人体艺术视频欧美日本| 女性生殖器流出的白浆| 久久久久久人妻| 亚洲国产精品成人久久小说| 在线观看美女被高潮喷水网站| 日本wwww免费看| 高清黄色对白视频在线免费看| 十八禁网站网址无遮挡| 不卡视频在线观看欧美| 99re6热这里在线精品视频| www.av在线官网国产| 成年人免费黄色播放视频| 伦精品一区二区三区| 亚洲,欧美,日韩| 老汉色av国产亚洲站长工具| 欧美日本中文国产一区发布| 中文字幕人妻丝袜一区二区 | 精品人妻偷拍中文字幕| 亚洲精品美女久久av网站| 国产精品国产av在线观看| 精品国产一区二区三区四区第35| 国产亚洲av片在线观看秒播厂| 丝袜脚勾引网站| videossex国产| 丝袜美足系列| 美女福利国产在线| 亚洲国产最新在线播放| 日日摸夜夜添夜夜爱| 国产高清国产精品国产三级| 黑人猛操日本美女一级片| 久久影院123| 熟女少妇亚洲综合色aaa.| 日韩av免费高清视频| 青春草亚洲视频在线观看| 日本av手机在线免费观看| 国产综合精华液| 高清不卡的av网站| av又黄又爽大尺度在线免费看| 熟女电影av网| 纯流量卡能插随身wifi吗| 欧美日韩视频高清一区二区三区二| 亚洲天堂av无毛| 国产成人精品福利久久| 国产精品国产三级国产专区5o| 国产精品 欧美亚洲| 成人毛片a级毛片在线播放| 成年人免费黄色播放视频| 在线天堂中文资源库| 十八禁网站网址无遮挡| 亚洲av电影在线进入| 国产欧美亚洲国产| av有码第一页| 久久国产精品大桥未久av| 精品国产一区二区三区久久久樱花| 美女脱内裤让男人舔精品视频| 亚洲成人av在线免费| 高清黄色对白视频在线免费看| 欧美精品一区二区大全| 亚洲欧美精品自产自拍| 欧美精品一区二区大全| 肉色欧美久久久久久久蜜桃| 国产精品av久久久久免费| 亚洲精品美女久久久久99蜜臀 | 另类亚洲欧美激情| 日韩 亚洲 欧美在线| 亚洲综合精品二区| 午夜免费鲁丝| 大片电影免费在线观看免费| 欧美日韩精品成人综合77777| 亚洲精品自拍成人| 精品国产乱码久久久久久男人| 日本黄色日本黄色录像| 80岁老熟妇乱子伦牲交| 久久精品国产鲁丝片午夜精品| 亚洲精品第二区| 日本欧美国产在线视频| 男男h啪啪无遮挡| 亚洲情色 制服丝袜| 丝袜喷水一区| 天堂中文最新版在线下载| 午夜精品国产一区二区电影| 9热在线视频观看99| 久久午夜福利片| 亚洲,一卡二卡三卡| 久久精品国产综合久久久| 国产精品久久久久久精品电影小说| 成人二区视频| 天天躁夜夜躁狠狠躁躁| 欧美 日韩 精品 国产| 老汉色∧v一级毛片| 伊人久久国产一区二区| 成年动漫av网址| 999久久久国产精品视频| 中文字幕亚洲精品专区| 国产精品无大码| 飞空精品影院首页| 综合色丁香网| 免费久久久久久久精品成人欧美视频| 一本色道久久久久久精品综合| 国产一区二区在线观看av| 桃花免费在线播放| 欧美国产精品一级二级三级| 色视频在线一区二区三区| 中文字幕av电影在线播放| 国产精品熟女久久久久浪| 国产精品一二三区在线看| 天天躁夜夜躁狠狠躁躁| 国产淫语在线视频| 熟女av电影| 亚洲经典国产精华液单| 免费高清在线观看日韩| 1024香蕉在线观看| 国产成人精品一,二区| 少妇 在线观看| 免费日韩欧美在线观看| 日韩欧美精品免费久久| 亚洲美女视频黄频| 欧美人与性动交α欧美精品济南到 | 丰满少妇做爰视频| 毛片一级片免费看久久久久| 黄色毛片三级朝国网站| 一区二区三区四区激情视频| av国产精品久久久久影院| 制服人妻中文乱码| 爱豆传媒免费全集在线观看| 日韩av免费高清视频| 国产精品欧美亚洲77777| 免费久久久久久久精品成人欧美视频| 亚洲精品乱久久久久久| 青春草亚洲视频在线观看| 一本大道久久a久久精品| 婷婷成人精品国产| 国产欧美日韩一区二区三区在线| 久久精品国产亚洲av天美| 久久婷婷青草| 国产欧美亚洲国产| 男的添女的下面高潮视频| 99热国产这里只有精品6| 国产精品欧美亚洲77777| 啦啦啦视频在线资源免费观看| 欧美日韩av久久| 999精品在线视频| 午夜影院在线不卡| 国产一区二区三区综合在线观看| 大话2 男鬼变身卡| 午夜福利网站1000一区二区三区| 精品国产乱码久久久久久男人| 边亲边吃奶的免费视频| 99热国产这里只有精品6| 欧美另类一区| 亚洲色图综合在线观看| 国产成人免费无遮挡视频| 日本欧美国产在线视频| 亚洲精品国产av蜜桃| 美女中出高潮动态图| 亚洲中文av在线| 老鸭窝网址在线观看| 久久久精品国产亚洲av高清涩受| 丰满乱子伦码专区| 夜夜骑夜夜射夜夜干| 18禁动态无遮挡网站| 在线观看免费视频网站a站| 亚洲成色77777| 久久99精品国语久久久| xxx大片免费视频| 久久韩国三级中文字幕| 性色av一级| 日本91视频免费播放| a级毛片黄视频| 99九九在线精品视频| 嫩草影院入口| 青青草视频在线视频观看| 有码 亚洲区| 一边亲一边摸免费视频| 国产一区二区在线观看av| 亚洲av在线观看美女高潮| 咕卡用的链子| 免费在线观看视频国产中文字幕亚洲 | 久久久久精品人妻al黑| 天天躁夜夜躁狠狠久久av| 欧美bdsm另类| 日本免费在线观看一区| 热99久久久久精品小说推荐| 波野结衣二区三区在线| av有码第一页| 午夜福利在线免费观看网站| av.在线天堂| 久久99精品国语久久久| 亚洲激情五月婷婷啪啪| 中文欧美无线码| 国产亚洲最大av| 一级a爱视频在线免费观看| 久久人人97超碰香蕉20202| 精品99又大又爽又粗少妇毛片| 成人黄色视频免费在线看| 亚洲成人一二三区av| 国产精品 欧美亚洲| 国产xxxxx性猛交| 99热网站在线观看| 日本av免费视频播放| 国产在线一区二区三区精| xxx大片免费视频| 侵犯人妻中文字幕一二三四区| 亚洲精品日韩在线中文字幕| 香蕉丝袜av| 日韩,欧美,国产一区二区三区| 在线精品无人区一区二区三| 男女免费视频国产| 欧美日韩成人在线一区二区| 久久久久久伊人网av| 男女边吃奶边做爰视频| 国产1区2区3区精品| 久久久久国产精品人妻一区二区| 国产亚洲精品第一综合不卡| 不卡视频在线观看欧美| 一级,二级,三级黄色视频| 在线看a的网站| 哪个播放器可以免费观看大片| 中文精品一卡2卡3卡4更新| 日本爱情动作片www.在线观看| 色播在线永久视频| 久久人人爽人人片av| 国产又爽黄色视频| 久久国产精品大桥未久av| 久久久亚洲精品成人影院| 最新中文字幕久久久久| 男女高潮啪啪啪动态图| 热re99久久国产66热| 成年人午夜在线观看视频| 久久精品国产a三级三级三级| www.自偷自拍.com| 另类精品久久| 伊人久久国产一区二区| 考比视频在线观看| 香蕉丝袜av| 最黄视频免费看| 日本猛色少妇xxxxx猛交久久| 观看美女的网站| 午夜福利影视在线免费观看| 97人妻天天添夜夜摸| 亚洲成人一二三区av| 欧美成人精品欧美一级黄| 99热全是精品| 两性夫妻黄色片| av.在线天堂| 老汉色av国产亚洲站长工具| 精品人妻在线不人妻| 亚洲图色成人| 亚洲成人一二三区av| 下体分泌物呈黄色| 久久午夜福利片| 午夜免费鲁丝| 99久久中文字幕三级久久日本| av卡一久久| 搡老乐熟女国产| 国产日韩一区二区三区精品不卡| 极品人妻少妇av视频| 久久精品亚洲av国产电影网| 国产高清不卡午夜福利| 欧美人与善性xxx| 久久亚洲国产成人精品v| 久久这里有精品视频免费| 亚洲三级黄色毛片| 国产不卡av网站在线观看| 亚洲视频免费观看视频| 成人午夜精彩视频在线观看| 女人高潮潮喷娇喘18禁视频| videos熟女内射| 欧美在线黄色| 国产xxxxx性猛交| 丁香六月天网| 久久精品国产鲁丝片午夜精品| 一边摸一边做爽爽视频免费| 欧美人与性动交α欧美软件| 亚洲国产欧美网| 综合色丁香网| 一区二区三区四区激情视频| 国产成人av激情在线播放| 国产一区二区 视频在线| 在线观看国产h片| www.熟女人妻精品国产| 有码 亚洲区| 免费久久久久久久精品成人欧美视频| 大片电影免费在线观看免费| 久久久国产一区二区| 国产成人91sexporn| 成人国语在线视频| 国产一区二区在线观看av| 成人手机av| 国产极品天堂在线| 精品人妻偷拍中文字幕| 一级片免费观看大全| 纵有疾风起免费观看全集完整版| 满18在线观看网站| 午夜福利在线观看免费完整高清在| 久久免费观看电影| 亚洲av福利一区| 2021少妇久久久久久久久久久| 涩涩av久久男人的天堂| 成人免费观看视频高清| 美女国产视频在线观看| av国产久精品久网站免费入址| 精品人妻一区二区三区麻豆| 日韩熟女老妇一区二区性免费视频| 国产欧美日韩综合在线一区二区| 最近中文字幕高清免费大全6| 另类精品久久| 亚洲成av片中文字幕在线观看 | 99热网站在线观看| 精品国产国语对白av| 久久精品国产鲁丝片午夜精品| 精品国产露脸久久av麻豆| 又粗又硬又长又爽又黄的视频| 少妇人妻精品综合一区二区| 欧美黄色片欧美黄色片| 国产一区二区三区av在线| 久久97久久精品| 天天躁狠狠躁夜夜躁狠狠躁| 日韩av在线免费看完整版不卡| 久久热在线av| 极品少妇高潮喷水抽搐| 成人毛片a级毛片在线播放| 亚洲精品国产色婷婷电影| 国产成人aa在线观看| 国产毛片在线视频| 日韩精品有码人妻一区| 免费在线观看视频国产中文字幕亚洲 | 男女啪啪激烈高潮av片| 国产人伦9x9x在线观看 | a级毛片在线看网站| av国产精品久久久久影院| 乱人伦中国视频| 免费观看在线日韩| 亚洲欧洲日产国产| 黄网站色视频无遮挡免费观看| videosex国产| 日韩大片免费观看网站| 国产片特级美女逼逼视频| 久久毛片免费看一区二区三区| 哪个播放器可以免费观看大片| 亚洲精品国产av成人精品| 久久精品久久精品一区二区三区| 午夜日韩欧美国产| av在线app专区| 中文天堂在线官网| 欧美精品av麻豆av| 免费看av在线观看网站| 亚洲中文av在线| 精品久久久久久电影网| 国产精品熟女久久久久浪| 国产淫语在线视频| 汤姆久久久久久久影院中文字幕| 国产综合精华液| 狠狠精品人妻久久久久久综合| 人体艺术视频欧美日本| 又粗又硬又长又爽又黄的视频| 国产精品国产三级国产专区5o| 人妻一区二区av| 国产熟女午夜一区二区三区| av电影中文网址| av视频免费观看在线观看| 精品亚洲成国产av| 亚洲精品aⅴ在线观看| 国产片内射在线| 国产欧美日韩一区二区三区在线| 成人亚洲欧美一区二区av| 国产男女内射视频| 国产精品一二三区在线看| 国产精品一区二区在线观看99| 欧美精品国产亚洲| 国产在视频线精品| 欧美精品高潮呻吟av久久| 国产精品国产三级专区第一集| 久久久久久久精品精品| 国产亚洲av片在线观看秒播厂| 亚洲激情五月婷婷啪啪| 啦啦啦在线免费观看视频4| 国产精品麻豆人妻色哟哟久久| 电影成人av| 午夜老司机福利剧场| 精品人妻在线不人妻| 极品人妻少妇av视频| 侵犯人妻中文字幕一二三四区| 女人久久www免费人成看片| 男女边吃奶边做爰视频| 免费黄网站久久成人精品| 成人毛片a级毛片在线播放| 一级片'在线观看视频| 久久综合国产亚洲精品| 新久久久久国产一级毛片| 母亲3免费完整高清在线观看 | 一边摸一边做爽爽视频免费| 成人亚洲欧美一区二区av| 成人漫画全彩无遮挡| 老女人水多毛片| 视频在线观看一区二区三区| 王馨瑶露胸无遮挡在线观看| 宅男免费午夜| 久久久久视频综合| 香蕉精品网在线| 男女边摸边吃奶| 国产精品久久久久成人av| 精品一品国产午夜福利视频| 在线观看免费日韩欧美大片| 黄色怎么调成土黄色| 999久久久国产精品视频| 90打野战视频偷拍视频| 美女国产视频在线观看| 婷婷色综合www| 男女无遮挡免费网站观看| 亚洲人成电影观看| 青青草视频在线视频观看| 美女视频免费永久观看网站| 免费大片黄手机在线观看| 亚洲欧洲国产日韩| 午夜日韩欧美国产| 在线观看免费视频网站a站| 亚洲av.av天堂| 久久久久久人妻| 亚洲熟女精品中文字幕| 久久人人爽人人片av| 精品福利永久在线观看| 卡戴珊不雅视频在线播放| 欧美人与善性xxx| 成人影院久久| 一级爰片在线观看| 日韩一卡2卡3卡4卡2021年| 女人被躁到高潮嗷嗷叫费观| 国产极品天堂在线| 丝袜脚勾引网站| 久久久久久久久久久免费av| 搡老乐熟女国产| 在线观看免费视频网站a站| 精品久久久久久电影网| 在线亚洲精品国产二区图片欧美| 欧美日韩成人在线一区二区| 下体分泌物呈黄色| 一级毛片黄色毛片免费观看视频| 国产av国产精品国产| 中文字幕人妻丝袜一区二区 | 国产片内射在线| 国产精品香港三级国产av潘金莲 | 1024香蕉在线观看| 国产亚洲午夜精品一区二区久久| 老女人水多毛片| 国产成人91sexporn| 久久久久久久国产电影| 国产精品偷伦视频观看了| 男人爽女人下面视频在线观看| 国产熟女欧美一区二区| 久久午夜福利片| 久久久国产精品麻豆| 啦啦啦视频在线资源免费观看| 人人妻人人爽人人添夜夜欢视频| 丰满乱子伦码专区| 美女中出高潮动态图| 十分钟在线观看高清视频www| 久久久久久久大尺度免费视频| 国产极品天堂在线| 久久久久国产精品人妻一区二区| 日本-黄色视频高清免费观看| 人妻人人澡人人爽人人| 最近中文字幕2019免费版| 久久精品国产综合久久久| 女性生殖器流出的白浆| 亚洲精品av麻豆狂野| 一边摸一边做爽爽视频免费| 国产精品国产av在线观看| 久久精品国产综合久久久| 久热这里只有精品99| 国产男人的电影天堂91| 国产淫语在线视频| 色哟哟·www| 国产淫语在线视频| 多毛熟女@视频| 国产精品国产av在线观看| 亚洲欧洲精品一区二区精品久久久 | 寂寞人妻少妇视频99o| 一级毛片我不卡| 黄色一级大片看看| 国产免费视频播放在线视频| 成人免费观看视频高清| 99久久综合免费| 99精国产麻豆久久婷婷| 春色校园在线视频观看| 老司机影院毛片| 9191精品国产免费久久| 亚洲成人av在线免费| 大片免费播放器 马上看| 黄色一级大片看看| 99久久综合免费| 国产成人一区二区在线| 亚洲av电影在线观看一区二区三区| 18禁国产床啪视频网站| 免费高清在线观看日韩| 日韩av免费高清视频| 在线观看免费视频网站a站| 国产精品三级大全| 精品视频人人做人人爽| 国产亚洲午夜精品一区二区久久| 久热这里只有精品99| 日韩av免费高清视频| 亚洲国产看品久久| 亚洲国产欧美在线一区| 精品视频人人做人人爽| 97在线视频观看| 国产精品国产av在线观看| 搡老乐熟女国产| 高清视频免费观看一区二区| xxxhd国产人妻xxx| 欧美激情极品国产一区二区三区| 久久99热这里只频精品6学生| 亚洲精品美女久久av网站| av卡一久久| 叶爱在线成人免费视频播放| 亚洲欧美清纯卡通| 国产成人精品无人区| 亚洲欧美一区二区三区国产| 亚洲av国产av综合av卡| 婷婷色av中文字幕| 国产欧美日韩综合在线一区二区| 久久精品国产a三级三级三级| 十分钟在线观看高清视频www| 国产亚洲精品第一综合不卡| 欧美精品国产亚洲| 另类亚洲欧美激情| 成人漫画全彩无遮挡| 亚洲男人天堂网一区| 在线观看免费高清a一片| 老司机影院毛片| 老司机影院成人| 亚洲精品日韩在线中文字幕| 热re99久久精品国产66热6| 亚洲欧美成人综合另类久久久| tube8黄色片| av又黄又爽大尺度在线免费看| 久久婷婷青草| 美女xxoo啪啪120秒动态图| 国产一区二区激情短视频 | 美女xxoo啪啪120秒动态图| 成年人免费黄色播放视频| 国产av国产精品国产| 色94色欧美一区二区| 久久国内精品自在自线图片| 亚洲av欧美aⅴ国产| 日韩 亚洲 欧美在线| 在线观看免费日韩欧美大片| 婷婷成人精品国产| 一区二区av电影网| 国产亚洲av片在线观看秒播厂| 在线观看三级黄色| 亚洲视频免费观看视频| 天美传媒精品一区二区| 波野结衣二区三区在线| 免费日韩欧美在线观看|