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

    質(zhì)量約束的三維模型建模方法

    2022-03-07 06:12:37卜寧遠(yuǎn)謝宇洋
    軟件工程 2022年3期

    卜寧遠(yuǎn) 謝宇洋

    摘? 要:在當(dāng)今的智能制造和工業(yè)軟件設(shè)計(jì)領(lǐng)域,CAD和CAE是極為重要的技術(shù)。但現(xiàn)有的CAD和CAE模型表示方法不統(tǒng)一,在數(shù)據(jù)交換上需要耗費(fèi)大量時(shí)間,造成了計(jì)算資源極大的浪費(fèi)。非均勻有理B樣條(NURBS)模型作為一種兼容CAD和CAE的模型表達(dá)方式,以樣條曲線為基函數(shù),無須進(jìn)行交換即可進(jìn)行等幾何分析。本文提出了一種基于質(zhì)量約束的NURBS體建模方法,將模型質(zhì)量作為建模時(shí)的約束,使構(gòu)建的模型符合等幾何分析的要求。本文以帶有復(fù)連通域的模型為例進(jìn)行等幾何分析,經(jīng)過與主流商業(yè)軟件的對(duì)比,最小值相同,最大值誤差在10%以內(nèi),運(yùn)行時(shí)間減少了4.61%,驗(yàn)證了此方法的正確性。

    關(guān)鍵詞:NURBS;建模方法;等幾何分析

    中圖分類號(hào):TP391.41? ? ?文獻(xiàn)標(biāo)識(shí)碼:A

    3D Modeling based on Mass Dependence

    BU Ningyuan, XIE Yuyang

    (School of Mechanical Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China)

    BNY19971997@163.com; 826647688@qq.com

    Abstract: CAD (Computer Aided Design) and CAE (Computer Aided Engineering) are extremely important technologies in today's intelligent manufacturing and industrial software design fields. However, representation methods of the existing CAD and CAE model vary, and it consumes a lot of time to exchange data, resulting in a great waste of computing resources. Non-uniform rational B-spline (NURBS) model, as a model expression compatible with CAD and CAE, takes spline curve as the basis function, and can perform isogeometric analysis without exchange. This paper proposes a NURBS modeling method based on mass dependence, which takes model mass as a constraint during modeling, so that the constructed model meets the requirements of isogeometric analysis. This paper takes a model with a complex connected domain as an example to perform isogeometric analysis. Compared with popular commercial software, the minimum value is the same, the maximum error is within 10%, and the running time decreases by 4.61%, which verifies the correctness of this method.

    Keywords: NURBS; modeling method; isogeometric analysis

    1? ?引言(Introduction)

    隨著第四次工業(yè)革命的到來,多個(gè)國家都提出了工業(yè)數(shù)字化轉(zhuǎn)型的戰(zhàn)略規(guī)劃,在智能制造領(lǐng)域,關(guān)鍵之一就是工業(yè)軟件[1]。對(duì)于產(chǎn)品設(shè)計(jì)而言,專業(yè)的設(shè)計(jì)和仿真軟件是必不可少的。為了降低設(shè)計(jì)成本,減少設(shè)計(jì)時(shí)間,將設(shè)計(jì)(CAD)與分析(CAE)集成一體化已是大勢所趨。但迄今為止,CAD與CAE仍分為兩個(gè)獨(dú)立的領(lǐng)域[2],因?yàn)椴捎玫哪P捅磉_(dá)方式有所不同,導(dǎo)致兩者之間產(chǎn)生了縫隙。當(dāng)今主流CAD模型的表達(dá)多為邊界表達(dá)B-Rep[3]和構(gòu)造實(shí)體CSG[4],而CAE模型為了適用于有限元分析(FEA)[5-6],一般采用網(wǎng)格表示法。因此在整個(gè)設(shè)計(jì)與分析的過程中,模型會(huì)在CAD和CAE兩種表示方式之間反復(fù)切換,拓?fù)浣Y(jié)構(gòu)被破壞,模型精度明顯降低,時(shí)間成本大幅度提高[7-8]。

    而NURBS參數(shù)化模型采用樣條曲線作為基函數(shù),能精確表示邊界和曲面,且能提供模型內(nèi)部的連續(xù)性參數(shù)化表達(dá)式,無須進(jìn)行CAD和CAE之間的反復(fù)轉(zhuǎn)化,產(chǎn)生的模型能直接用來進(jìn)行等幾何分析。等幾何方法是有限元方法的繼承與發(fā)展,該方法采用樣條曲線作為基函數(shù),具有較高的幾何精確性。

    但是,體參數(shù)化模型為了適用于等幾何分析,必須要表示為零虧格的雙變量張量曲面或三變量張量體,且相鄰的體參數(shù)化單元之間的曲面要一致。因此需要對(duì)不滿足等幾何分析要求的體參數(shù)化模型進(jìn)行分割處理,且等幾何分析對(duì)模型的要求較高,分割后建立的模型還要在幾何上滿足等幾何分析的要求。

    2? ?相關(guān)研究(Related works)

    對(duì)于CAD和CAE融合,可以進(jìn)行CAD二次開發(fā),在成熟的商業(yè)平臺(tái)上進(jìn)行融合;也可以使用重建法,在已知CAD模型邊界的基礎(chǔ)上得到內(nèi)部控制點(diǎn)的信息,但需要進(jìn)行大量的運(yùn)算。體參數(shù)化建模不同于以上兩種,具有計(jì)算簡單、創(chuàng)建方便的優(yōu)點(diǎn)。

    在CAD模型剖分方面,相關(guān)學(xué)者已做了大量工作。大致方向可分為兩種:(1)多邊形凸分解,將多邊形剖分為三角形、四邊形及任意凸多邊形集合;(2)網(wǎng)格分解,主要針對(duì)有限元模型,將CAD轉(zhuǎn)換為CAE形式表達(dá),如狄利克雷三角剖分法。由于等幾何分析的特性,模型不需要再被轉(zhuǎn)換成CAE表達(dá),因此體參數(shù)化剖分主要采用第一種思路。

    常見的多邊形凸分解算法有兩種:(1)幾何單元區(qū)域剖分算法,如梯形分割算法[9];(2)考慮是否具有Steiner點(diǎn)的凸分解算法,如編碼劃分算法[10]和加權(quán)剖分算法[11]。但上述剖分算法僅在數(shù)學(xué)方面給出了多邊形分割方案,沒有考慮剖分結(jié)果是否符合分析標(biāo)準(zhǔn),因此十分有必要評(píng)估四邊形質(zhì)量及其對(duì)等幾何分析結(jié)果的影響。

    本文提出了一種以NURBS體參數(shù)化模型質(zhì)量為約束的建模方法,改進(jìn)了適用于等幾何分析的模型質(zhì)量評(píng)估方法,將模型質(zhì)量嵌入建模過程之中。首先遞歸分割模型截面,在分割過程中施加質(zhì)量約束,使分割后的截面滿足等幾何分析的要求;然后使用掃描、拉伸、旋轉(zhuǎn)、放樣等操作進(jìn)行體參數(shù)化映射,得到符合條件的體參數(shù)化模型。

    3? 體參數(shù)化建模與模型質(zhì)量(Volume parametric modeling and model mass)

    對(duì)于機(jī)械零件而言,大部分模型都可以定義截面及路徑,通過掃描、拉伸、旋轉(zhuǎn)、放樣等操作完成模型構(gòu)建。因此,為了使一般體參數(shù)化模型符合等幾何分析的要求,要先對(duì)截面進(jìn)行分割處理,將其分割成四邊形的集合,以滿足等幾何分析中雙變量張量曲面的要求。

    3.1? ?截面剖分算法

    本文使用的是遞歸剖分的方法,將整個(gè)模型控制點(diǎn)看作無向圖,采用深度優(yōu)先遍歷,逐層向下尋找合適的剖分方案。該算法時(shí)間空間復(fù)雜度與傳統(tǒng)算法相同,但在遞歸過程中加上了體參數(shù)建模所要求的約束,從數(shù)學(xué)問題轉(zhuǎn)化為工程問題,因此最終結(jié)果更符合等幾何分析的要求。

    設(shè)算法輸入為多邊形子域集合POL,輸出為四邊形集合Quad,算法流程圖如圖1所示。

    算法流程如下:

    Step 1:若不為空,取作為,否則結(jié)束并返回。

    Step 2:求解可行剖分線集合,設(shè)共有 種剖分連線方式,初始選擇剖分線。

    Step 3:,表示該子域不滿足任何剖分條件;,表示已遍歷所有的剖分連線方式。如結(jié)果不滿足條件,則回溯至上一個(gè)剖分的子域,設(shè)為,。

    Step 4:根據(jù)當(dāng)前選擇的剖分連線方式,創(chuàng)建一條新的NURBS曲線,同時(shí)獲得兩個(gè)新子域、,刪除當(dāng)前被剖分的子域。

    Step 5:為凸四邊形,存入集合;為五邊形及以上,存入集合尾部;若均不符合,則回溯至上一個(gè)剖分的子域,。

    Step 6:遞歸執(zhí)行Step 1至Step 5。

    3.2? ?剖分中的模型質(zhì)量

    在等幾何分析中,計(jì)算域的參數(shù)化對(duì)應(yīng)于有限元分析中的網(wǎng)格生成,對(duì)等幾何分析結(jié)果的準(zhǔn)確性是十分重要的[11],需要將模型從物理域映射到參數(shù)域進(jìn)行分析。使用雅可比矩陣將模型從物理域映射到參數(shù)域,雅可比矩陣如下:

    根據(jù)已有質(zhì)量理論,最標(biāo)準(zhǔn)的參數(shù)化模型應(yīng)滿足以下三個(gè)條件:(1)模型沒有自交;(2)等參單元盡可能正交;(3)等參單元盡可能均勻。這樣會(huì)使雅可比的值接近或等于1,參數(shù)域映射的效果最佳。而扭曲、自交、不均勻的曲面均會(huì)使雅可比的質(zhì)量變差從而破壞參數(shù)域映射,也會(huì)使等幾何分析求解偏微分方程的過程變得十分困難,因此在剖分中我們要盡量避免出現(xiàn)這種曲面。對(duì)于不合格的曲面,需要制定一系列檢測方法。

    對(duì)于六面體網(wǎng)格,主流方法已經(jīng)給出了一系列檢測指標(biāo):對(duì)角線度量、縮放雅可比、扭曲度和寬高比等。對(duì)于NURBS曲面而言,可以在曲面上取若干等參單元看作四邊形進(jìn)行質(zhì)量評(píng)估。我們在此選擇了區(qū)域面積A、寬高比asp和扭曲度dist三個(gè)指標(biāo)作為衡量四邊形的依據(jù),以下是這三個(gè)指標(biāo)的詳細(xì)定義。

    3.2.1? ?區(qū)域面積

    如圖2所示,我們把每一個(gè)等參單元看作一個(gè)小的四邊形,為了衡量它是否均勻,是否接近最好的分析模型,我們首先定義邊向量,然后將其分為四個(gè)區(qū)域,定義坐標(biāo)軸。

    其中,為四邊形的角點(diǎn)坐標(biāo)。為了衡量四個(gè)角點(diǎn)區(qū)域是否在統(tǒng)一平面上,定義每一個(gè)角點(diǎn)處的法向量:

    為了方便計(jì)算,設(shè)它們的單位向量:

    定義一個(gè)中心法向量作為全局的法向量和單位向量,分別為:

    將整個(gè)區(qū)域按照坐標(biāo)的上下左右軸分成四塊,如圖3(a)所示,則可以根據(jù)已知的參數(shù)求得每塊的區(qū)域面積:

    根據(jù)幾何知識(shí)可以知道,若(即該四邊形某一區(qū)域的面積)小于0,則該片一定是退化的,如圖3(b)所示,即為凹四邊形。

    對(duì)于NURBS曲面而言,由于邊界是曲線,因此直接按照角點(diǎn)求法向量是沒有意義的,但每個(gè)面的控制點(diǎn)和等參線可以構(gòu)成若干個(gè)小的等參單元。等參單元最好的情況是四條邊兩兩正交,而退化的等參單元一定不滿足邊邊正交,因此會(huì)拉低整體的圖形質(zhì)量。因此,所有等參單元面積大于0的數(shù)量越多越好。當(dāng)凹四邊形出現(xiàn)時(shí),必然會(huì)有一個(gè)區(qū)域面積等于或小于零。凹多邊形的出現(xiàn)使得等幾何分析中的參數(shù)域映射和求解偏微分方程的效果很差,是一定不能出現(xiàn)在分析模型之中的,因此利用區(qū)域面積法可以很簡單地判別出不符合情況的條件之一,且與普通四邊形剖分中的射線法相比更具有魯棒性。

    3.2.2? ?寬高比

    NURBS等參單元的結(jié)構(gòu)越均勻越好,因此可以定義高寬比來評(píng)價(jià)等參單元是否均勻:

    其中,A為四邊形面積,寬高比值域?yàn)?,?jīng)驗(yàn)表明為可接受域。寬高比排除的是看起來很“窄”的四邊形,這類四邊形的出現(xiàn)會(huì)使控制點(diǎn)的分布極其不均勻,導(dǎo)致等幾何參數(shù)域的映射效果很差,也是一種不能出現(xiàn)在體參數(shù)化模型上的情況之一。寬高比的求解只需要在區(qū)域面積的基礎(chǔ)上再加上一步運(yùn)算即可,消耗的資源很少,產(chǎn)生的效果卻很好。

    3.2.3? ?扭曲度

    為了衡量四邊形是否扭曲,首先要建立一個(gè)參照系,設(shè)一個(gè)位于坐標(biāo)原點(diǎn)的正方形的面積;再根據(jù)式(3)可以算出所有高斯積分點(diǎn)的最小雅可比值,則扭曲度定義為:

    扭曲是衡量從參數(shù)空間到世界坐標(biāo)映射行為的良好程度的一種度量,值域?yàn)?,越接?效果越好,因此將可接受域設(shè)為。

    根據(jù)以上三個(gè)指標(biāo),可以較為準(zhǔn)確地檢查剖分截面的均勻性、正交性和扭曲度,且計(jì)算也十分簡單,不會(huì)占用太多時(shí)間和空間。這幾個(gè)指標(biāo)可以在遞歸時(shí)的判別條件中加入,在連接分割線后,依次檢查這三個(gè)值,若有不在接收域的情況則不合格,否則放入備選庫,作為可能的選項(xiàng)之一。

    3.3? ?體參數(shù)化映射

    得到了剖分好的NURBS四邊形曲面后,可以使用掃描、拉伸、旋轉(zhuǎn)、放樣等操作進(jìn)行體參數(shù)化映射。拉伸體模型可以通過截面控制點(diǎn)沿直線方向的仿射變換得到,在CAD中多用于凸臺(tái)等結(jié)構(gòu)的創(chuàng)建。如果模型的路徑不是直線而是圓弧,則可以進(jìn)行旋轉(zhuǎn)體映射。若拉伸的路徑不是直線也不是規(guī)則圓弧,則可以用掃掠體造型。

    對(duì)于無法用一個(gè)截面生成的模型,可以考慮采用放樣體造型,沿曲線路徑插值放樣體的造型方法和掃掠體相似,需要先統(tǒng)一各個(gè)截面的階數(shù)、節(jié)點(diǎn)矢量及控制點(diǎn)數(shù)量,確保截面與該點(diǎn)曲線切矢垂直,再通過反求控制點(diǎn)即可得到沿曲線路徑的放樣體。而不沿路徑插值的放樣方式則無須保證各截面與曲線切矢垂直。

    對(duì)于機(jī)械零件而言,大部分模型都可以通過這幾種建模方法得到。在得到了體參數(shù)化模型之后便可以直接進(jìn)行等幾何分析,而不需要再花大量時(shí)間進(jìn)行網(wǎng)格劃分。

    4? ?實(shí)例(Examples)

    本部分以一個(gè)帶有通孔的箱體配件為例,該配件經(jīng)常用于減速箱上,頂部受力,底部有固定約束,是典型的機(jī)械部件之一。我們在其上表面施加力約束,下表面固定,進(jìn)行靜力學(xué)分析。由于配件帶有通孔,因此進(jìn)行剖分之前需要在圓上插入若干的點(diǎn),再進(jìn)行剖分,且圓孔周圍的分割線也要保持一定的弧度以免與圓相交,再進(jìn)行拉伸操作,最后得到的體參數(shù)化模型如圖4所示。再使用MATLAB軟件對(duì)模型進(jìn)行等幾何分析,得出的結(jié)果與ANSYS軟件對(duì)比如表1和圖5所示,位移的最小值相同,最大值相差在10%以內(nèi),同時(shí)平均運(yùn)行時(shí)間減少了4.61%,證明了本文方法的有效性。

    5? ?結(jié)論(Conclusion)

    本文為CAD和CAE的無縫集成提供了一種新的思路,即以NURBS樣條曲線模型代替原有的CSG和B-Rep模型,這樣使得進(jìn)行分析的時(shí)候無須數(shù)據(jù)格式和表達(dá)方式的轉(zhuǎn)換,大大降低了時(shí)間復(fù)雜度和內(nèi)存空間的消耗。本文以箱體配件為例,展示了如何對(duì)帶有復(fù)連通域的圖形進(jìn)行剖分和建模,經(jīng)過與商業(yè)軟件的對(duì)比得到了較好的結(jié)果。

    但此方法仍然存在一些問題。等幾何分析對(duì)模型的要求較高,一些不規(guī)則圖形如點(diǎn)陣材料等還無法通過簡單的掃描、拉伸、旋轉(zhuǎn)、放樣等簡單的體參數(shù)化映射方法得到。同樣的,對(duì)于一些細(xì)微的特征如倒角、圓角和小孔等也無法完全保留,這也是我們下一步主要研究的方向。另外,為了提高體參數(shù)化建模保留原有幾何特征的能力,使構(gòu)建的模型更加精確,對(duì)體參數(shù)化模型布爾運(yùn)算也需要進(jìn)一步研究。

    參考文獻(xiàn)(References)

    [1] 寧振波.力筑工業(yè)之基重鑄智造靈魂[J].軟件導(dǎo)刊,2021,? 20(01):7-11,6.

    [2] SMIT M S, BRONSVOORT W F. Integration of design and analysis models[J]. Computer-Aided Design and Applications, 2009, 6(6):795-808.

    [3] KIM B C, MUN D. Feature-based simplification of boundary representation models using sequential iterative volume? ?decomposition[J]. Computers & Graphics, 2014, 38:97-107.

    [4] 羅月童,樊曉菁,俞盛朋,等.基于面殼封閉的B-Rep至CSG轉(zhuǎn)換算法[J].計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào),2014,26(10):1673-1680.

    [5] AVERSA R, PETRESCU F I T, PETRESCU R V V, et al. Biofidel FEA modeling of customized hybrid biological hip joint design part Ⅱ: Flexible stem trabecular prostheses[J]. American Journal of Biochemistry and Biotechnology, 2016, 12(4):23-29.

    [6] KASIK D J, BUXTON W, FERGUSON D R. Ten CAD challenges[J]. IEEE Computer Graphics and Applications, 2005, 25(2):81-92.

    [7] PARK H S, DANG X P. Structural optimization based on CAD-CAE integration and metamodeling techniques[J].Computer-Aided Design, 2010, 42(10):889-902.

    [8] PAN Z, WANG X, TENG R, et al. Computer-aided design-while-engineering technology in top-down modeling of mechanical product[J]. Computers in Industry, 2016, 75: 151-161.

    [9] 王金敏,劉季燁,方沂.基于梯形分解的不規(guī)則多邊形干涉算法[J].工程圖學(xué)學(xué)報(bào),2005(06):52-57.

    [10] 肖忠暉,盧振榮,張謙.簡單多邊形凸單元剖分的編碼算法[J].計(jì)算機(jī)學(xué)報(bào),1996(06):477-481.

    [11] 王博,李笑牛,李華.一種加權(quán)剖分簡單多邊形為三角形和凸四邊形子域的算法[J].中國圖象圖形學(xué)報(bào),2002,7(005):? ?486-490.

    作者簡介:

    卜寧遠(yuǎn)(1997-),男,碩士生.研究領(lǐng)域:CAD/CAE.

    謝宇洋(1997-),男,碩士生.研究領(lǐng)域:CAD/CAE.

    亚洲精品中文字幕一二三四区| 日韩国内少妇激情av| 午夜日韩欧美国产| 国产麻豆69| 亚洲人成电影观看| 久久精品成人免费网站| 大陆偷拍与自拍| 国产三级在线视频| 变态另类成人亚洲欧美熟女 | 亚洲激情在线av| 国产激情欧美一区二区| 久久久久国产一级毛片高清牌| 视频区图区小说| 精品免费久久久久久久清纯| 国产欧美日韩一区二区三| 黄色 视频免费看| 又紧又爽又黄一区二区| 怎么达到女性高潮| 精品久久久久久电影网| 99riav亚洲国产免费| 亚洲国产精品sss在线观看 | 国产精品98久久久久久宅男小说| 国产精品一区二区在线不卡| 少妇的丰满在线观看| av超薄肉色丝袜交足视频| 91成人精品电影| 中文字幕精品免费在线观看视频| 色婷婷av一区二区三区视频| 99国产综合亚洲精品| 国产激情欧美一区二区| 婷婷丁香在线五月| 久久久国产成人免费| 亚洲黑人精品在线| 亚洲国产毛片av蜜桃av| 999精品在线视频| 12—13女人毛片做爰片一| 丝袜美足系列| 久久久久久久久免费视频了| 很黄的视频免费| www.自偷自拍.com| 国产无遮挡羞羞视频在线观看| 在线观看免费视频日本深夜| 亚洲avbb在线观看| 电影成人av| 精品人妻在线不人妻| 国产av一区二区精品久久| 丁香欧美五月| 一级黄色大片毛片| 亚洲国产中文字幕在线视频| 久久午夜亚洲精品久久| 中出人妻视频一区二区| 免费观看人在逋| 亚洲人成电影观看| 欧美一级毛片孕妇| 国产麻豆69| 亚洲精品一卡2卡三卡4卡5卡| 丁香六月欧美| 国产三级在线视频| 国产一区二区激情短视频| 免费看十八禁软件| 国产成人一区二区三区免费视频网站| 午夜a级毛片| 亚洲黑人精品在线| av网站在线播放免费| av超薄肉色丝袜交足视频| 在线国产一区二区在线| 丰满人妻熟妇乱又伦精品不卡| 亚洲情色 制服丝袜| 久热这里只有精品99| 激情视频va一区二区三区| 免费观看精品视频网站| 国产精品久久视频播放| 亚洲av成人av| 国产又色又爽无遮挡免费看| 国产男靠女视频免费网站| 成年版毛片免费区| x7x7x7水蜜桃| 午夜福利在线观看吧| 亚洲久久久国产精品| 在线永久观看黄色视频| 国产99白浆流出| 欧美久久黑人一区二区| www.精华液| 亚洲成国产人片在线观看| 精品久久久久久电影网| 亚洲人成电影免费在线| 中文欧美无线码| 久久精品成人免费网站| 免费女性裸体啪啪无遮挡网站| 成人免费观看视频高清| 午夜a级毛片| 国产精品99久久99久久久不卡| 一本大道久久a久久精品| 无人区码免费观看不卡| 国产成人精品久久二区二区91| 99国产综合亚洲精品| 国产精品日韩av在线免费观看 | 午夜福利影视在线免费观看| 久久青草综合色| 欧美大码av| 天堂影院成人在线观看| 可以免费在线观看a视频的电影网站| 亚洲人成伊人成综合网2020| 桃红色精品国产亚洲av| 琪琪午夜伦伦电影理论片6080| 国产真人三级小视频在线观看| 女人高潮潮喷娇喘18禁视频| 97碰自拍视频| 亚洲成人久久性| 国产亚洲av高清不卡| 午夜激情av网站| 亚洲色图av天堂| 高清在线国产一区| 精品欧美一区二区三区在线| 精品熟女少妇八av免费久了| 一级作爱视频免费观看| 精品福利永久在线观看| 国产精品久久久久久人妻精品电影| 亚洲精品一二三| 男女下面进入的视频免费午夜 | 电影成人av| 国产精品 国内视频| 99久久国产精品久久久| 久久天躁狠狠躁夜夜2o2o| 亚洲黑人精品在线| 国产av精品麻豆| 久久久国产欧美日韩av| 三上悠亚av全集在线观看| 国产伦一二天堂av在线观看| 欧美日韩视频精品一区| 国产成人系列免费观看| 午夜91福利影院| 精品无人区乱码1区二区| 一进一出抽搐动态| 天堂动漫精品| 一级a爱视频在线免费观看| 99久久国产精品久久久| 80岁老熟妇乱子伦牲交| 一区二区日韩欧美中文字幕| 天堂俺去俺来也www色官网| 丰满迷人的少妇在线观看| 国产精品成人在线| 亚洲成av片中文字幕在线观看| 一进一出抽搐gif免费好疼 | 亚洲成人精品中文字幕电影 | 中文字幕精品免费在线观看视频| 成人三级做爰电影| 日韩大码丰满熟妇| 老汉色av国产亚洲站长工具| 国产一区二区三区在线臀色熟女 | 嫩草影院精品99| 久久精品国产亚洲av香蕉五月| 最好的美女福利视频网| 悠悠久久av| 亚洲精品在线美女| 91成人精品电影| 亚洲国产精品999在线| 99国产极品粉嫩在线观看| 久久久久久人人人人人| 欧美日韩亚洲高清精品| 亚洲精品中文字幕在线视频| 亚洲 欧美 日韩 在线 免费| 免费一级毛片在线播放高清视频 | 欧美大码av| 久久久久久久久免费视频了| 夜夜躁狠狠躁天天躁| 中文字幕色久视频| 校园春色视频在线观看| 在线av久久热| 欧美日韩亚洲国产一区二区在线观看| 亚洲第一欧美日韩一区二区三区| 操美女的视频在线观看| 亚洲一区高清亚洲精品| 一个人免费在线观看的高清视频| 亚洲国产欧美一区二区综合| 免费高清在线观看日韩| 又黄又粗又硬又大视频| 亚洲一区二区三区不卡视频| 一边摸一边抽搐一进一出视频| 黄网站色视频无遮挡免费观看| 欧美午夜高清在线| 成在线人永久免费视频| 国产精品久久视频播放| 可以免费在线观看a视频的电影网站| 精品日产1卡2卡| 亚洲成人免费av在线播放| 国产伦人伦偷精品视频| 久久久久久久久久久久大奶| 久久人妻av系列| 中文字幕高清在线视频| 日本免费一区二区三区高清不卡 | 黄网站色视频无遮挡免费观看| 久久狼人影院| 久久亚洲真实| 国产欧美日韩综合在线一区二区| 中文字幕最新亚洲高清| 亚洲国产精品合色在线| 亚洲精品一区av在线观看| 亚洲狠狠婷婷综合久久图片| 免费在线观看亚洲国产| 亚洲精品久久午夜乱码| 伊人久久大香线蕉亚洲五| 久久精品人人爽人人爽视色| 首页视频小说图片口味搜索| 亚洲中文日韩欧美视频| 亚洲自拍偷在线| 国产亚洲欧美在线一区二区| 操出白浆在线播放| 黑人操中国人逼视频| 免费在线观看影片大全网站| 啦啦啦免费观看视频1| 视频区欧美日本亚洲| 真人一进一出gif抽搐免费| 精品午夜福利视频在线观看一区| 国产成人免费无遮挡视频| 一进一出抽搐gif免费好疼 | 国产精品永久免费网站| av在线天堂中文字幕 | 99精国产麻豆久久婷婷| 久久久久九九精品影院| 18禁黄网站禁片午夜丰满| 亚洲免费av在线视频| 亚洲色图综合在线观看| 身体一侧抽搐| 日本欧美视频一区| 大陆偷拍与自拍| 午夜免费激情av| 国产亚洲欧美精品永久| 国产乱人伦免费视频| 亚洲欧美日韩无卡精品| 成人亚洲精品一区在线观看| 免费不卡黄色视频| 国产伦一二天堂av在线观看| 巨乳人妻的诱惑在线观看| 18禁美女被吸乳视频| 精品久久久久久电影网| 国产伦一二天堂av在线观看| 村上凉子中文字幕在线| 成年人免费黄色播放视频| 18禁黄网站禁片午夜丰满| 热re99久久国产66热| 亚洲一区二区三区色噜噜 | 大香蕉久久成人网| av国产精品久久久久影院| 国产亚洲精品第一综合不卡| 久久99一区二区三区| 欧美久久黑人一区二区| 乱人伦中国视频| 露出奶头的视频| 亚洲熟妇熟女久久| 麻豆国产av国片精品| 最近最新中文字幕大全电影3 | 在线观看一区二区三区激情| 91麻豆精品激情在线观看国产 | 精品熟女少妇八av免费久了| 美女福利国产在线| 国产三级在线视频| 午夜福利免费观看在线| 狂野欧美激情性xxxx| 香蕉久久夜色| 悠悠久久av| 天堂√8在线中文| 国产99白浆流出| 日韩人妻精品一区2区三区| 人人妻人人澡人人看| 夫妻午夜视频| 看黄色毛片网站| 欧美日韩av久久| 国产欧美日韩一区二区三区在线| 亚洲精品久久成人aⅴ小说| 免费在线观看日本一区| 亚洲精品av麻豆狂野| 麻豆成人av在线观看| 一区二区三区国产精品乱码| 成人手机av| 国产精品 欧美亚洲| 热re99久久国产66热| 国产精品99久久99久久久不卡| 日韩免费av在线播放| 久久国产精品影院| 精品卡一卡二卡四卡免费| 亚洲色图av天堂| 精品少妇一区二区三区视频日本电影| 神马国产精品三级电影在线观看 | 一级a爱视频在线免费观看| 三上悠亚av全集在线观看| 欧美黄色片欧美黄色片| 国产日韩一区二区三区精品不卡| 欧美性长视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 成人亚洲精品一区在线观看| 亚洲成国产人片在线观看| 一个人免费在线观看的高清视频| tocl精华| 亚洲 欧美一区二区三区| 亚洲美女黄片视频| 国产精品 国内视频| 久久国产亚洲av麻豆专区| 啦啦啦免费观看视频1| 午夜福利在线观看吧| 国产97色在线日韩免费| 手机成人av网站| 国产精品国产av在线观看| 女性生殖器流出的白浆| e午夜精品久久久久久久| 亚洲国产精品合色在线| 国产精品 国内视频| 在线观看www视频免费| 热99国产精品久久久久久7| 精品卡一卡二卡四卡免费| 手机成人av网站| 午夜视频精品福利| 亚洲欧美日韩无卡精品| 在线观看一区二区三区| 免费高清视频大片| 性色av乱码一区二区三区2| 18美女黄网站色大片免费观看| 一区在线观看完整版| 9色porny在线观看| 美女福利国产在线| 亚洲成人免费电影在线观看| 免费久久久久久久精品成人欧美视频| 嫁个100分男人电影在线观看| 久久午夜亚洲精品久久| 精品欧美一区二区三区在线| 亚洲一区二区三区不卡视频| 国产成人欧美| 日本精品一区二区三区蜜桃| 国产精品久久久久久人妻精品电影| 无人区码免费观看不卡| 成人18禁高潮啪啪吃奶动态图| 久热爱精品视频在线9| 男女下面进入的视频免费午夜 | 亚洲黑人精品在线| 91麻豆av在线| 叶爱在线成人免费视频播放| 欧美日韩视频精品一区| 日韩欧美在线二视频| 色老头精品视频在线观看| 另类亚洲欧美激情| 国产精品秋霞免费鲁丝片| 亚洲午夜精品一区,二区,三区| 91精品三级在线观看| 中出人妻视频一区二区| 手机成人av网站| 人人妻人人爽人人添夜夜欢视频| 两人在一起打扑克的视频| 欧美日韩黄片免| 亚洲第一青青草原| 满18在线观看网站| av天堂久久9| 亚洲精品在线美女| videosex国产| 精品福利永久在线观看| 成熟少妇高潮喷水视频| 99精品久久久久人妻精品| 热99国产精品久久久久久7| a级片在线免费高清观看视频| 一区在线观看完整版| 亚洲自偷自拍图片 自拍| 亚洲精品一二三| 国产国语露脸激情在线看| 色老头精品视频在线观看| 超碰97精品在线观看| 免费女性裸体啪啪无遮挡网站| 亚洲欧美精品综合久久99| 一区二区三区激情视频| 久久中文字幕人妻熟女| 日韩大码丰满熟妇| 日日摸夜夜添夜夜添小说| 国产欧美日韩一区二区三区在线| 亚洲人成77777在线视频| 亚洲人成电影免费在线| 国产成人精品久久二区二区免费| 久久久久久免费高清国产稀缺| 90打野战视频偷拍视频| 国产成人欧美在线观看| 国产成+人综合+亚洲专区| 亚洲国产精品合色在线| 亚洲成人免费电影在线观看| 国产黄a三级三级三级人| 97超级碰碰碰精品色视频在线观看| 欧美精品啪啪一区二区三区| 久久精品aⅴ一区二区三区四区| 免费搜索国产男女视频| 9色porny在线观看| av在线播放免费不卡| 日韩大码丰满熟妇| 亚洲专区字幕在线| 男女午夜视频在线观看| 夜夜爽天天搞| 男女床上黄色一级片免费看| 免费在线观看黄色视频的| 一区二区三区激情视频| 90打野战视频偷拍视频| 日韩大尺度精品在线看网址 | 久久国产亚洲av麻豆专区| 亚洲专区中文字幕在线| www.999成人在线观看| 99久久国产精品久久久| 日本 av在线| 欧美一区二区精品小视频在线| 天天影视国产精品| 久久天躁狠狠躁夜夜2o2o| 在线播放国产精品三级| 国产亚洲精品久久久久久毛片| 久久狼人影院| 嫁个100分男人电影在线观看| 最新美女视频免费是黄的| 亚洲成人免费电影在线观看| 亚洲男人天堂网一区| 国产av又大| 91老司机精品| 亚洲全国av大片| 露出奶头的视频| 极品人妻少妇av视频| 真人做人爱边吃奶动态| 伊人久久大香线蕉亚洲五| 日韩欧美免费精品| 天堂影院成人在线观看| 国产精品乱码一区二三区的特点 | 激情视频va一区二区三区| 国产一区二区在线av高清观看| 交换朋友夫妻互换小说| 淫妇啪啪啪对白视频| 亚洲精品中文字幕一二三四区| 女警被强在线播放| av电影中文网址| av欧美777| 欧美日韩一级在线毛片| 欧美 亚洲 国产 日韩一| avwww免费| 欧美日韩视频精品一区| 日韩国内少妇激情av| av在线播放免费不卡| 久久精品亚洲av国产电影网| 精品高清国产在线一区| 精品国产一区二区三区四区第35| 亚洲欧美一区二区三区黑人| 久久草成人影院| 欧美激情久久久久久爽电影 | 免费在线观看影片大全网站| 亚洲男人天堂网一区| 18禁黄网站禁片午夜丰满| 欧美 亚洲 国产 日韩一| 另类亚洲欧美激情| 日本黄色视频三级网站网址| 波多野结衣高清无吗| 岛国视频午夜一区免费看| 久久九九热精品免费| 精品久久久久久成人av| 99国产精品99久久久久| 亚洲自拍偷在线| 日韩人妻精品一区2区三区| 看免费av毛片| 国产一区在线观看成人免费| 亚洲aⅴ乱码一区二区在线播放 | 国产在线精品亚洲第一网站| 国产欧美日韩一区二区三区在线| 99国产精品一区二区蜜桃av| 9色porny在线观看| 波多野结衣高清无吗| 露出奶头的视频| 亚洲一区高清亚洲精品| 激情在线观看视频在线高清| 宅男免费午夜| 色婷婷av一区二区三区视频| 97人妻天天添夜夜摸| 久久久国产精品麻豆| 久久伊人香网站| 精品国产国语对白av| 18禁国产床啪视频网站| 夜夜夜夜夜久久久久| 国产免费男女视频| 精品国产一区二区三区四区第35| 国产色视频综合| 国产免费男女视频| 国产男靠女视频免费网站| 亚洲精品粉嫩美女一区| 久久九九热精品免费| 亚洲自偷自拍图片 自拍| 亚洲av美国av| 91麻豆av在线| 三上悠亚av全集在线观看| 亚洲黑人精品在线| 真人做人爱边吃奶动态| av视频免费观看在线观看| 可以免费在线观看a视频的电影网站| 精品一区二区三区视频在线观看免费 | 国产日韩一区二区三区精品不卡| 天天躁夜夜躁狠狠躁躁| 一级,二级,三级黄色视频| 久久精品国产综合久久久| 午夜福利一区二区在线看| 免费在线观看亚洲国产| 亚洲欧美精品综合久久99| 国产精品一区二区免费欧美| 色尼玛亚洲综合影院| 国产av一区在线观看免费| 人妻丰满熟妇av一区二区三区| 久久精品亚洲熟妇少妇任你| 日本五十路高清| 日韩三级视频一区二区三区| 免费高清视频大片| 十八禁网站免费在线| 女性生殖器流出的白浆| 中文欧美无线码| 我的亚洲天堂| 成人av一区二区三区在线看| 99久久久亚洲精品蜜臀av| 欧美老熟妇乱子伦牲交| 午夜福利在线观看吧| 精品无人区乱码1区二区| 免费观看精品视频网站| 精品一区二区三区视频在线观看免费 | 看片在线看免费视频| 欧美最黄视频在线播放免费 | 少妇粗大呻吟视频| 国产精品九九99| 日本撒尿小便嘘嘘汇集6| 欧美最黄视频在线播放免费 | 五月开心婷婷网| 久久热在线av| 80岁老熟妇乱子伦牲交| 女人被狂操c到高潮| 精品久久久精品久久久| 中文字幕人妻丝袜一区二区| 国产成人啪精品午夜网站| 国产av在哪里看| 亚洲精品粉嫩美女一区| 色婷婷av一区二区三区视频| 丰满饥渴人妻一区二区三| 99国产精品免费福利视频| 日韩欧美在线二视频| 国产亚洲精品综合一区在线观看 | 国产乱人伦免费视频| 日韩欧美一区二区三区在线观看| 99热国产这里只有精品6| 亚洲熟女毛片儿| 中文字幕人妻丝袜制服| 亚洲情色 制服丝袜| 成人亚洲精品av一区二区 | 在线观看午夜福利视频| 精品高清国产在线一区| 久久精品影院6| www.精华液| 亚洲国产精品999在线| 久久亚洲真实| 高潮久久久久久久久久久不卡| 一区二区日韩欧美中文字幕| 久久精品国产99精品国产亚洲性色 | 午夜福利在线观看吧| 人人妻人人添人人爽欧美一区卜| 男女做爰动态图高潮gif福利片 | 91av网站免费观看| 十分钟在线观看高清视频www| 成人国产一区最新在线观看| 中文字幕另类日韩欧美亚洲嫩草| 亚洲中文日韩欧美视频| 国产黄a三级三级三级人| 日日干狠狠操夜夜爽| 在线观看免费日韩欧美大片| 老熟妇乱子伦视频在线观看| 日韩欧美一区视频在线观看| 性欧美人与动物交配| 一级黄色大片毛片| 欧美在线一区亚洲| 岛国在线观看网站| 嫩草影视91久久| 在线观看日韩欧美| 久久久久精品国产欧美久久久| 国产99白浆流出| 九色亚洲精品在线播放| 18禁美女被吸乳视频| 久久人妻av系列| 婷婷精品国产亚洲av在线| 天堂中文最新版在线下载| 中文字幕高清在线视频| 一二三四社区在线视频社区8| 天堂影院成人在线观看| 麻豆久久精品国产亚洲av | 老司机亚洲免费影院| 亚洲七黄色美女视频| 两性午夜刺激爽爽歪歪视频在线观看 | 精品一区二区三区视频在线观看免费 | 中文字幕人妻丝袜制服| 老司机深夜福利视频在线观看| 国产av一区二区精品久久| 高清在线国产一区| 美女大奶头视频| 久久久久久久午夜电影 | 久久亚洲精品不卡| 美女扒开内裤让男人捅视频| 精品一区二区三卡| 久久草成人影院| 国产精品乱码一区二三区的特点 | 日韩国内少妇激情av| 97碰自拍视频| 久久精品亚洲av国产电影网| 亚洲中文字幕日韩| 母亲3免费完整高清在线观看| 看黄色毛片网站| 999久久久国产精品视频| 国产精品久久视频播放| 久久伊人香网站| 一本大道久久a久久精品| 一进一出好大好爽视频| 又大又爽又粗| 国产高清videossex| 久久婷婷成人综合色麻豆| 欧美日韩瑟瑟在线播放| 国产精品影院久久| 日本a在线网址| 黑人欧美特级aaaaaa片| 国产一区在线观看成人免费|