• <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級曲面控制頂點排列的影響*
    久久99精品国语久久久| 日韩欧美精品免费久久| 亚洲国产欧美在线一区| 国产 一区精品| 男女免费视频国产| 午夜福利在线观看免费完整高清在| 边亲边吃奶的免费视频| 伊人亚洲综合成人网| 黄色一级大片看看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲国产成人一精品久久久| 免费在线观看视频国产中文字幕亚洲 | 伦精品一区二区三区| 一本久久精品| 丝袜美足系列| 男女国产视频网站| 婷婷色综合www| 十分钟在线观看高清视频www| 国产精品免费视频内射| 丝袜人妻中文字幕| 精品亚洲乱码少妇综合久久| 日韩免费高清中文字幕av| 亚洲人成77777在线视频| 久久午夜福利片| 欧美黄色片欧美黄色片| 国产视频首页在线观看| 丰满饥渴人妻一区二区三| 天堂8中文在线网| 五月开心婷婷网| 午夜影院在线不卡| 亚洲精品日本国产第一区| 国产成人a∨麻豆精品| 国产1区2区3区精品| 午夜福利视频在线观看免费| 高清视频免费观看一区二区| 日本vs欧美在线观看视频| 国产欧美日韩一区二区三区在线| 日产精品乱码卡一卡2卡三| 欧美亚洲日本最大视频资源| 18在线观看网站| 久久久国产一区二区| 日本午夜av视频| 亚洲精品久久成人aⅴ小说| 亚洲国产欧美网| 视频在线观看一区二区三区| 久久久久久人妻| 亚洲 欧美一区二区三区| 国产成人午夜福利电影在线观看| 国产精品偷伦视频观看了| www日本在线高清视频| 日韩一本色道免费dvd| 欧美成人精品欧美一级黄| 七月丁香在线播放| 99re6热这里在线精品视频| 搡老乐熟女国产| 国产在线免费精品| 啦啦啦中文免费视频观看日本| 纯流量卡能插随身wifi吗| 波多野结衣一区麻豆| 亚洲av综合色区一区| 欧美日韩精品成人综合77777| 亚洲av在线观看美女高潮| 在线观看国产h片| 国产福利在线免费观看视频| 中文字幕制服av| 男人操女人黄网站| 母亲3免费完整高清在线观看 | 久久久久精品久久久久真实原创| 晚上一个人看的免费电影| 又大又黄又爽视频免费| 一级爰片在线观看| 99九九在线精品视频| 成人黄色视频免费在线看| 一区福利在线观看| 欧美日韩国产mv在线观看视频| 精品第一国产精品| 777久久人妻少妇嫩草av网站| 新久久久久国产一级毛片| 国产一区二区三区综合在线观看| a级毛片在线看网站| 可以免费在线观看a视频的电影网站 | 亚洲色图 男人天堂 中文字幕| 97在线视频观看| 99久国产av精品国产电影| 国产乱人偷精品视频| 久久久久人妻精品一区果冻| 欧美精品亚洲一区二区| 美女午夜性视频免费| 国产熟女欧美一区二区| 亚洲av男天堂| 大片免费播放器 马上看| 国产一级毛片在线| 国产伦理片在线播放av一区| 视频在线观看一区二区三区| 黑人猛操日本美女一级片| 亚洲av免费高清在线观看| 日日啪夜夜爽| videosex国产| 国产极品粉嫩免费观看在线| 91在线精品国自产拍蜜月| 麻豆精品久久久久久蜜桃| 国产精品香港三级国产av潘金莲 | 不卡视频在线观看欧美| 国产精品偷伦视频观看了| 国产人伦9x9x在线观看 | av免费观看日本| 中国国产av一级| 亚洲欧美精品综合一区二区三区 | 欧美日韩成人在线一区二区| 男女高潮啪啪啪动态图| 水蜜桃什么品种好| 建设人人有责人人尽责人人享有的| 精品亚洲成国产av| 国产高清不卡午夜福利| 国产欧美日韩综合在线一区二区| 搡女人真爽免费视频火全软件| 青春草亚洲视频在线观看| 青青草视频在线视频观看| 777久久人妻少妇嫩草av网站| 精品一区二区三卡| 黄频高清免费视频| 成人国产av品久久久| 国精品久久久久久国模美| 啦啦啦啦在线视频资源| 午夜福利,免费看| 麻豆av在线久日| 黑人欧美特级aaaaaa片| 三级国产精品片| 亚洲人成网站在线观看播放| 国产黄色免费在线视频| 国产人伦9x9x在线观看 | 国产成人91sexporn| 一级毛片我不卡| 日本欧美视频一区| 亚洲欧洲国产日韩| 日本wwww免费看| 国产成人精品无人区| 精品国产乱码久久久久久男人| 不卡视频在线观看欧美| 成年人免费黄色播放视频| 亚洲人成77777在线视频| 欧美精品人与动牲交sv欧美| 大香蕉久久成人网| 七月丁香在线播放| 9191精品国产免费久久| 黄色 视频免费看| 久久久国产欧美日韩av| 久久av网站| 国产高清不卡午夜福利| 汤姆久久久久久久影院中文字幕| 久久久久精品久久久久真实原创| 欧美国产精品一级二级三级| 青青草视频在线视频观看| 春色校园在线视频观看| 久久精品久久久久久久性| 精品视频人人做人人爽| 国产av精品麻豆| 日日啪夜夜爽| 女性被躁到高潮视频| 久久久久精品性色| 老汉色∧v一级毛片| 日韩电影二区| 夫妻性生交免费视频一级片| 2022亚洲国产成人精品| 亚洲国产av影院在线观看| 搡老乐熟女国产| 日韩免费高清中文字幕av| 国产精品久久久久久精品古装| 亚洲国产毛片av蜜桃av| 国产一区二区三区综合在线观看| 久久99一区二区三区| 观看av在线不卡| www.精华液| 一区二区三区乱码不卡18| 另类精品久久| 久久国产精品大桥未久av| 午夜激情久久久久久久| 国产精品久久久久久精品古装| 亚洲成色77777| 中文精品一卡2卡3卡4更新| 一二三四在线观看免费中文在| 欧美黄色片欧美黄色片| 一级毛片 在线播放| 久久久久精品性色| 免费大片黄手机在线观看| 丝袜喷水一区| 老鸭窝网址在线观看| 久久热在线av| 久久精品国产鲁丝片午夜精品| 国产成人精品久久久久久| 欧美少妇被猛烈插入视频| 另类精品久久| 久久人妻熟女aⅴ| 男女高潮啪啪啪动态图| 韩国精品一区二区三区| 美女大奶头黄色视频| 欧美日韩视频高清一区二区三区二| 熟妇人妻不卡中文字幕| 国产女主播在线喷水免费视频网站| 黑人欧美特级aaaaaa片| 丝袜喷水一区| 国产一区二区三区av在线| 咕卡用的链子| 国产av一区二区精品久久| 亚洲熟女精品中文字幕| 日韩免费高清中文字幕av| 免费看不卡的av| 国产精品国产三级专区第一集| 免费观看无遮挡的男女| 男人添女人高潮全过程视频| 欧美bdsm另类| 久久久久精品性色| 国产免费现黄频在线看| 国产精品三级大全| 熟妇人妻不卡中文字幕| 一级,二级,三级黄色视频| 9色porny在线观看| 国产亚洲av片在线观看秒播厂| a级毛片在线看网站| 久久 成人 亚洲| 天天影视国产精品| 最黄视频免费看| 久久国内精品自在自线图片| 老汉色av国产亚洲站长工具| 国产男女内射视频| 国产精品一区二区在线观看99| 中文精品一卡2卡3卡4更新| 成人二区视频| 成人毛片a级毛片在线播放| 日本av手机在线免费观看| 男女午夜视频在线观看| 制服诱惑二区| 女性被躁到高潮视频| 深夜精品福利| 爱豆传媒免费全集在线观看| 超碰97精品在线观看| 免费观看av网站的网址| 一区二区三区四区激情视频| 久久精品熟女亚洲av麻豆精品| 日本色播在线视频| 水蜜桃什么品种好| 精品国产露脸久久av麻豆| 久久影院123| 十八禁高潮呻吟视频| 国产极品天堂在线| 不卡av一区二区三区| 久久久久精品人妻al黑| 欧美少妇被猛烈插入视频| 亚洲伊人色综图| 久久久亚洲精品成人影院| 免费观看性生交大片5| 国产成人精品一,二区| 日本av免费视频播放| 我要看黄色一级片免费的| 亚洲成av片中文字幕在线观看 | 久久久精品94久久精品| 亚洲国产精品国产精品| 免费高清在线观看视频在线观看| 国语对白做爰xxxⅹ性视频网站| 国产毛片在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品蜜桃在线观看| 午夜日本视频在线| 日韩免费高清中文字幕av| 欧美+日韩+精品| xxx大片免费视频| 999久久久国产精品视频| 中文天堂在线官网| 欧美日韩国产mv在线观看视频| 蜜桃在线观看..| 亚洲av男天堂| 久久久久久久国产电影| 18在线观看网站| 99久久精品国产国产毛片| 欧美人与善性xxx| 国产精品人妻久久久影院| 纯流量卡能插随身wifi吗| 91午夜精品亚洲一区二区三区| 亚洲婷婷狠狠爱综合网| 久久久久精品人妻al黑| 国产高清国产精品国产三级| 国产日韩欧美在线精品| 亚洲欧美成人综合另类久久久| 一级a爱视频在线免费观看| 丝袜美腿诱惑在线| 一区二区三区四区激情视频| 黄片小视频在线播放| 啦啦啦在线观看免费高清www| 国产乱来视频区| 在线精品无人区一区二区三| 婷婷色麻豆天堂久久| 人人妻人人澡人人爽人人夜夜| 日本午夜av视频| 久久 成人 亚洲| 韩国精品一区二区三区| 日韩熟女老妇一区二区性免费视频| 国产一区二区在线观看av| 欧美日韩国产mv在线观看视频| 国产成人精品一,二区| 边亲边吃奶的免费视频| 97在线视频观看| 天天操日日干夜夜撸| 各种免费的搞黄视频| 国产精品嫩草影院av在线观看| 日韩精品免费视频一区二区三区| 亚洲精品久久久久久婷婷小说| 黄片播放在线免费| 亚洲美女搞黄在线观看| 国产精品av久久久久免费| 国产男女内射视频| 性色av一级| 日韩熟女老妇一区二区性免费视频| 色吧在线观看| 水蜜桃什么品种好| 不卡视频在线观看欧美| 菩萨蛮人人尽说江南好唐韦庄| 久久国内精品自在自线图片| 少妇的逼水好多| 欧美日韩一级在线毛片| 一本大道久久a久久精品| kizo精华| 免费在线观看黄色视频的| 亚洲精品日本国产第一区| 国产不卡av网站在线观看| 看非洲黑人一级黄片| 午夜福利影视在线免费观看| 人妻一区二区av| 精品久久久久久电影网| 最近2019中文字幕mv第一页| 亚洲成色77777| 丰满乱子伦码专区| 久久精品国产a三级三级三级| 成年美女黄网站色视频大全免费| 黄片播放在线免费| 少妇精品久久久久久久| 大香蕉久久成人网| av国产久精品久网站免费入址| 国产精品女同一区二区软件| 老汉色∧v一级毛片| 久久国产精品大桥未久av| 亚洲内射少妇av| 性色av一级| 久久久国产精品麻豆| 少妇的丰满在线观看| 国产成人aa在线观看| 日韩在线高清观看一区二区三区| 国产极品粉嫩免费观看在线| 欧美国产精品va在线观看不卡| 18禁国产床啪视频网站| 日韩制服骚丝袜av| 亚洲国产色片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 香蕉国产在线看| 一本久久精品| 热99久久久久精品小说推荐| 国产欧美日韩综合在线一区二区| 曰老女人黄片| 亚洲男人天堂网一区| 国产一区亚洲一区在线观看| 精品亚洲成a人片在线观看| 精品国产一区二区久久| 黑人巨大精品欧美一区二区蜜桃| 国产国语露脸激情在线看| 亚洲av综合色区一区| 一区二区三区精品91| 日韩一区二区三区影片| 亚洲经典国产精华液单| 在线观看三级黄色| 久久久久久免费高清国产稀缺| 亚洲三级黄色毛片| 中文字幕人妻丝袜制服| 中文字幕色久视频| 亚洲av电影在线观看一区二区三区| 少妇精品久久久久久久| 亚洲色图 男人天堂 中文字幕| 2021少妇久久久久久久久久久| 在线天堂中文资源库| 夫妻午夜视频| 国产欧美日韩一区二区三区在线| 日韩人妻精品一区2区三区| 日本爱情动作片www.在线观看| 日本91视频免费播放| 亚洲精品自拍成人| 国产成人免费无遮挡视频| 在线观看免费视频网站a站| 欧美人与性动交α欧美精品济南到 | 免费av中文字幕在线| 18禁观看日本| 人成视频在线观看免费观看| 黄色毛片三级朝国网站| 日韩电影二区| 亚洲欧美成人综合另类久久久| 在线免费观看不下载黄p国产| 老司机影院毛片| 成人18禁高潮啪啪吃奶动态图| 精品亚洲成国产av| 观看av在线不卡| 亚洲欧美精品综合一区二区三区 | 亚洲精品国产一区二区精华液| 在线亚洲精品国产二区图片欧美| 精品人妻偷拍中文字幕| 午夜91福利影院| 欧美激情 高清一区二区三区| 久久久久国产一级毛片高清牌| 王馨瑶露胸无遮挡在线观看| 麻豆精品久久久久久蜜桃| 精品国产一区二区三区久久久樱花| 两个人看的免费小视频| 婷婷色综合www| 人妻一区二区av| 国产男女内射视频| 水蜜桃什么品种好| 免费女性裸体啪啪无遮挡网站| 最近中文字幕高清免费大全6| 男女无遮挡免费网站观看| 日韩一卡2卡3卡4卡2021年| 少妇 在线观看| 天天操日日干夜夜撸| 一区二区三区乱码不卡18| 黄色配什么色好看| 91久久精品国产一区二区三区| 亚洲欧美日韩另类电影网站| 国产高清不卡午夜福利| 精品久久蜜臀av无| 亚洲美女搞黄在线观看| av福利片在线| 亚洲av综合色区一区| 伊人亚洲综合成人网| 两个人看的免费小视频| 亚洲欧美色中文字幕在线| 久久精品久久久久久噜噜老黄| 一二三四在线观看免费中文在| 国产高清不卡午夜福利| 天美传媒精品一区二区| 久久这里有精品视频免费| 一区在线观看完整版| 成人漫画全彩无遮挡| 成年动漫av网址| 亚洲伊人久久精品综合| 一级片免费观看大全| 国产精品麻豆人妻色哟哟久久| 黄片无遮挡物在线观看| 不卡av一区二区三区| 伦理电影大哥的女人| 中国国产av一级| h视频一区二区三区| 美女国产视频在线观看| 免费看av在线观看网站| 亚洲激情五月婷婷啪啪| 欧美人与性动交α欧美软件| 亚洲国产日韩一区二区| 人人澡人人妻人| 亚洲精华国产精华液的使用体验| 午夜激情久久久久久久| 久久久久久久亚洲中文字幕| 久久久精品国产亚洲av高清涩受| 国产一级毛片在线| 国产精品偷伦视频观看了| 国产精品99久久99久久久不卡 | 欧美97在线视频| 亚洲美女视频黄频| 久久久精品区二区三区| 一级片免费观看大全| 99久久精品国产国产毛片| 久久精品人人爽人人爽视色| 久久 成人 亚洲| 久久久国产精品麻豆| av国产精品久久久久影院| 亚洲精品一二三| 亚洲国产欧美在线一区| 精品人妻一区二区三区麻豆| 看非洲黑人一级黄片| 蜜桃在线观看..| 精品久久蜜臀av无| 午夜福利视频在线观看免费| av片东京热男人的天堂| 狠狠精品人妻久久久久久综合| 国产精品国产三级专区第一集| 水蜜桃什么品种好| 成人毛片60女人毛片免费| 9191精品国产免费久久| 黄色配什么色好看| 国产综合精华液| 国产一区二区激情短视频 | 国产女主播在线喷水免费视频网站| 日韩一区二区三区影片| 国产精品秋霞免费鲁丝片| av在线老鸭窝| 激情五月婷婷亚洲| av有码第一页| 中文乱码字字幕精品一区二区三区| 午夜日本视频在线| 狠狠精品人妻久久久久久综合| 少妇的丰满在线观看| 久久毛片免费看一区二区三区| 啦啦啦中文免费视频观看日本| 国产有黄有色有爽视频| 侵犯人妻中文字幕一二三四区| 69精品国产乱码久久久| av在线老鸭窝| 精品第一国产精品| 99国产综合亚洲精品| 美女国产高潮福利片在线看| 日本色播在线视频| 国产野战对白在线观看| 中文字幕最新亚洲高清| 亚洲av电影在线进入| a 毛片基地| 观看av在线不卡| 亚洲国产毛片av蜜桃av| 久久综合国产亚洲精品| 精品国产露脸久久av麻豆| 少妇人妻精品综合一区二区| 久久久久精品久久久久真实原创| 男女无遮挡免费网站观看| 日韩成人av中文字幕在线观看| 国产精品.久久久| 美女国产高潮福利片在线看| 男人操女人黄网站| 女的被弄到高潮叫床怎么办| 国产 一区精品| 99久久人妻综合| 精品午夜福利在线看| 深夜精品福利| 我要看黄色一级片免费的| freevideosex欧美| 久久人人爽人人片av| 18在线观看网站| 如何舔出高潮| 中文精品一卡2卡3卡4更新| 日韩精品免费视频一区二区三区| av网站在线播放免费| 少妇的逼水好多| 在线观看免费高清a一片| 制服丝袜香蕉在线| 在线观看www视频免费| 久久这里只有精品19| 久久精品aⅴ一区二区三区四区 | 亚洲精华国产精华液的使用体验| 王馨瑶露胸无遮挡在线观看| 亚洲精品视频女| 久久精品国产亚洲av天美| 亚洲精品国产一区二区精华液| 国产一区二区三区综合在线观看| 中文字幕精品免费在线观看视频| 久久精品亚洲av国产电影网| 精品亚洲成国产av| 久久综合国产亚洲精品| 狠狠婷婷综合久久久久久88av| 侵犯人妻中文字幕一二三四区| 一区二区av电影网| 午夜福利乱码中文字幕| 在线观看免费高清a一片| 母亲3免费完整高清在线观看 | 欧美日韩成人在线一区二区| 精品亚洲成a人片在线观看| 在线观看国产h片| 国产av码专区亚洲av| 激情视频va一区二区三区| 亚洲成国产人片在线观看| 亚洲成av片中文字幕在线观看 | 女性生殖器流出的白浆| 丰满少妇做爰视频| 国产男女超爽视频在线观看| 午夜福利影视在线免费观看| 日本欧美国产在线视频| 精品国产乱码久久久久久男人| 色播在线永久视频| 亚洲第一区二区三区不卡| 国产精品嫩草影院av在线观看| 男女下面插进去视频免费观看| www.av在线官网国产| av有码第一页| 又黄又粗又硬又大视频| 黄频高清免费视频| 亚洲国产色片| 亚洲精品日韩在线中文字幕| 久久久久国产精品人妻一区二区| 亚洲av国产av综合av卡| 国产乱来视频区| 日本色播在线视频| 久久精品国产亚洲av天美| 国产97色在线日韩免费| 精品一区二区三卡| 国产亚洲最大av| 亚洲精品美女久久久久99蜜臀 | 欧美人与善性xxx| 免费观看无遮挡的男女| 丝袜美腿诱惑在线| 亚洲欧美中文字幕日韩二区| 视频在线观看一区二区三区| 成人二区视频| 日韩不卡一区二区三区视频在线| 高清在线视频一区二区三区| 亚洲国产欧美在线一区| 久久毛片免费看一区二区三区| 人人妻人人添人人爽欧美一区卜| 波野结衣二区三区在线| www日本在线高清视频| 欧美亚洲 丝袜 人妻 在线| 晚上一个人看的免费电影| av福利片在线| 久久国产精品男人的天堂亚洲| 看免费av毛片| 成人毛片a级毛片在线播放| 国产在线一区二区三区精| 国产日韩一区二区三区精品不卡| 制服丝袜香蕉在线| 免费黄网站久久成人精品| 国产片内射在线| 天堂中文最新版在线下载| 晚上一个人看的免费电影| av一本久久久久| av.在线天堂| 成年av动漫网址|