• <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.

    最近手机中文字幕大全| 日本欧美国产在线视频| 午夜福利高清视频| 日本欧美国产在线视频| 国产精品福利在线免费观看| 久久久久精品国产欧美久久久| 亚洲成人精品中文字幕电影| 联通29元200g的流量卡| 亚洲精品一区av在线观看| 精品无人区乱码1区二区| 91久久精品电影网| 亚洲欧美精品综合久久99| 少妇丰满av| 国产成人a区在线观看| 精品乱码久久久久久99久播| 国产毛片a区久久久久| 成人鲁丝片一二三区免费| 一级毛片电影观看 | 久久草成人影院| 国产精品一及| 欧美区成人在线视频| 麻豆精品久久久久久蜜桃| 亚洲精品粉嫩美女一区| 精品熟女少妇av免费看| 亚洲第一电影网av| 男人和女人高潮做爰伦理| 日韩三级伦理在线观看| 成人特级黄色片久久久久久久| 国产成人a区在线观看| 麻豆国产av国片精品| 乱人视频在线观看| 国产免费一级a男人的天堂| 啦啦啦韩国在线观看视频| 国产真实乱freesex| 哪里可以看免费的av片| 午夜福利在线在线| 欧美激情国产日韩精品一区| 欧美成人精品欧美一级黄| 一夜夜www| 国产伦精品一区二区三区视频9| 1000部很黄的大片| 老师上课跳d突然被开到最大视频| 最后的刺客免费高清国语| 91麻豆精品激情在线观看国产| 国产精品乱码一区二三区的特点| 天堂影院成人在线观看| 老司机影院成人| 91久久精品国产一区二区三区| 亚洲自拍偷在线| 亚洲精品一卡2卡三卡4卡5卡| 国产伦一二天堂av在线观看| 久久国内精品自在自线图片| www日本黄色视频网| 波多野结衣巨乳人妻| 丰满乱子伦码专区| 欧美又色又爽又黄视频| 两性午夜刺激爽爽歪歪视频在线观看| 在线国产一区二区在线| 午夜福利视频1000在线观看| 国产黄a三级三级三级人| 亚洲精品一区av在线观看| 欧美激情久久久久久爽电影| 一卡2卡三卡四卡精品乱码亚洲| 欧美成人一区二区免费高清观看| 国产精品人妻久久久影院| 国产av在哪里看| 久久婷婷人人爽人人干人人爱| 三级毛片av免费| 亚洲国产精品久久男人天堂| 国产日本99.免费观看| 久久亚洲精品不卡| 成人欧美大片| a级毛色黄片| 亚洲精品日韩在线中文字幕 | 男女那种视频在线观看| 蜜桃亚洲精品一区二区三区| 51国产日韩欧美| 高清毛片免费观看视频网站| 男女之事视频高清在线观看| 我要搜黄色片| 免费搜索国产男女视频| 69人妻影院| 亚洲国产精品成人综合色| 一级黄片播放器| 国产成人a区在线观看| 波野结衣二区三区在线| 久久99热这里只有精品18| 男女之事视频高清在线观看| 内地一区二区视频在线| 国产欧美日韩精品一区二区| 精品一区二区三区视频在线观看免费| 十八禁网站免费在线| 在线观看午夜福利视频| 亚洲人成网站高清观看| av.在线天堂| 久久久成人免费电影| 亚洲五月天丁香| 欧美高清性xxxxhd video| 亚洲无线观看免费| 婷婷色综合大香蕉| 亚洲aⅴ乱码一区二区在线播放| 男人狂女人下面高潮的视频| 少妇高潮的动态图| 国产精品国产三级国产av玫瑰| 国产精品久久久久久av不卡| 99热这里只有是精品50| 精品人妻偷拍中文字幕| 久久久精品94久久精品| 欧美丝袜亚洲另类| 国产午夜精品论理片| av在线观看视频网站免费| 中国国产av一级| 久久精品国产亚洲网站| 熟妇人妻久久中文字幕3abv| 日本 av在线| 天天躁日日操中文字幕| 国产又黄又爽又无遮挡在线| 12—13女人毛片做爰片一| 国产成人福利小说| 最新在线观看一区二区三区| 国产黄片美女视频| 丰满人妻一区二区三区视频av| 高清毛片免费观看视频网站| 看片在线看免费视频| 又粗又爽又猛毛片免费看| 国产精品日韩av在线免费观看| 国产美女午夜福利| 小蜜桃在线观看免费完整版高清| 欧美又色又爽又黄视频| 深爱激情五月婷婷| 直男gayav资源| 熟女电影av网| 噜噜噜噜噜久久久久久91| 51国产日韩欧美| 人人妻人人看人人澡| 波多野结衣高清无吗| 噜噜噜噜噜久久久久久91| 久久精品国产亚洲网站| 黄色配什么色好看| 久久99热这里只有精品18| 在线观看午夜福利视频| 可以在线观看毛片的网站| 亚洲精品亚洲一区二区| 日本爱情动作片www.在线观看 | 国产蜜桃级精品一区二区三区| 狠狠狠狠99中文字幕| 中出人妻视频一区二区| 亚洲在线自拍视频| 精品一区二区三区人妻视频| 国产一区二区在线av高清观看| 简卡轻食公司| 亚洲av成人精品一区久久| 国产精品亚洲美女久久久| 国产av在哪里看| 午夜a级毛片| 免费观看人在逋| 国产日本99.免费观看| 欧美激情国产日韩精品一区| 自拍偷自拍亚洲精品老妇| 激情 狠狠 欧美| 村上凉子中文字幕在线| 国产精品人妻久久久影院| 插阴视频在线观看视频| 精品久久久久久久久久久久久| 国语自产精品视频在线第100页| avwww免费| 听说在线观看完整版免费高清| 51国产日韩欧美| 免费av观看视频| 此物有八面人人有两片| 亚洲欧美日韩卡通动漫| 中文字幕久久专区| 午夜福利18| а√天堂www在线а√下载| 成熟少妇高潮喷水视频| 成年女人看的毛片在线观看| 国产精品精品国产色婷婷| 久久久国产成人精品二区| 少妇熟女aⅴ在线视频| 亚洲欧美成人综合另类久久久 | 露出奶头的视频| 欧美3d第一页| 你懂的网址亚洲精品在线观看 | 婷婷六月久久综合丁香| 日韩 亚洲 欧美在线| 日日摸夜夜添夜夜添av毛片| 日日撸夜夜添| 一级黄片播放器| 久久久久久久久久成人| 97超视频在线观看视频| 色视频www国产| 色综合色国产| 97在线视频观看| 99久国产av精品| 国产黄片美女视频| 久久久精品94久久精品| 成年av动漫网址| 在线天堂最新版资源| 亚洲内射少妇av| 国产亚洲欧美98| 国语自产精品视频在线第100页| 亚洲av美国av| 国产69精品久久久久777片| 色综合色国产| 九九热线精品视视频播放| 国产v大片淫在线免费观看| 亚洲av二区三区四区| 99热只有精品国产| 午夜爱爱视频在线播放| 大香蕉久久网| 色哟哟·www| 美女高潮的动态| 国产精品1区2区在线观看.| 十八禁网站免费在线| 日本撒尿小便嘘嘘汇集6| 看免费成人av毛片| 国产成人福利小说| 久久精品国产亚洲av天美| 日本黄大片高清| 亚洲最大成人av| 欧美另类亚洲清纯唯美| 欧美一区二区精品小视频在线| 听说在线观看完整版免费高清| 99热这里只有精品一区| 一区二区三区高清视频在线| 精品久久久久久久末码| 亚洲丝袜综合中文字幕| 国产精品久久久久久精品电影| 可以在线观看的亚洲视频| 大香蕉久久网| av在线天堂中文字幕| 亚洲图色成人| 深夜a级毛片| 国产亚洲av嫩草精品影院| 亚洲精华国产精华液的使用体验 | 69av精品久久久久久| 老女人水多毛片| av天堂在线播放| 亚洲婷婷狠狠爱综合网| 成人永久免费在线观看视频| 国产国拍精品亚洲av在线观看| 免费看日本二区| 在线观看免费视频日本深夜| 在线观看美女被高潮喷水网站| 国产乱人偷精品视频| 国产爱豆传媒在线观看| 夜夜爽天天搞| 久久久久免费精品人妻一区二区| 亚洲18禁久久av| 夜夜夜夜夜久久久久| 在线观看66精品国产| 五月玫瑰六月丁香| 国产91av在线免费观看| 日韩欧美 国产精品| 成人午夜高清在线视频| 久久久欧美国产精品| 国产伦精品一区二区三区视频9| 别揉我奶头~嗯~啊~动态视频| 看免费成人av毛片| 国国产精品蜜臀av免费| 观看免费一级毛片| 精华霜和精华液先用哪个| 国产精华一区二区三区| 国产成人freesex在线 | 免费看光身美女| 亚洲av免费在线观看| 啦啦啦啦在线视频资源| 国产乱人视频| 在现免费观看毛片| 婷婷精品国产亚洲av| 91久久精品国产一区二区三区| 精品久久久久久成人av| 日韩欧美精品v在线| 男插女下体视频免费在线播放| 久久久久免费精品人妻一区二区| 国产亚洲91精品色在线| 欧美高清性xxxxhd video| 日韩av不卡免费在线播放| 一边摸一边抽搐一进一小说| 亚洲国产精品合色在线| 欧美一区二区亚洲| 赤兔流量卡办理| 亚洲美女视频黄频| a级毛色黄片| 亚洲中文字幕日韩| 最近中文字幕高清免费大全6| 久久久a久久爽久久v久久| 国产精品福利在线免费观看| 99在线视频只有这里精品首页| 久久久久国产网址| 成人漫画全彩无遮挡| 天堂网av新在线| 少妇丰满av| 欧美xxxx性猛交bbbb| 国产精品乱码一区二三区的特点| 国产精品爽爽va在线观看网站| 久久精品国产鲁丝片午夜精品| 2021天堂中文幕一二区在线观| 99久国产av精品国产电影| 国产成人freesex在线 | 亚洲成人av在线免费| 免费观看人在逋| 日韩一区二区视频免费看| 久久99热这里只有精品18| 97超级碰碰碰精品色视频在线观看| 久久精品91蜜桃| 狠狠狠狠99中文字幕| 少妇的逼水好多| 日韩欧美三级三区| 在线播放国产精品三级| 欧美日韩在线观看h| 久久久久精品国产欧美久久久| 国产真实伦视频高清在线观看| 午夜爱爱视频在线播放| 1000部很黄的大片| 一本一本综合久久| 精品不卡国产一区二区三区| 激情 狠狠 欧美| 黄色欧美视频在线观看| 日日撸夜夜添| 国产黄色视频一区二区在线观看 | 国产综合懂色| 亚洲在线观看片| 精品一区二区三区av网在线观看| 午夜福利成人在线免费观看| 男女那种视频在线观看| 免费黄网站久久成人精品| 99riav亚洲国产免费| 欧美xxxx黑人xx丫x性爽| 观看美女的网站| 男女视频在线观看网站免费| 日韩av在线大香蕉| 床上黄色一级片| 国产黄色小视频在线观看| 国国产精品蜜臀av免费| 露出奶头的视频| 深爱激情五月婷婷| 99九九线精品视频在线观看视频| 久久久久久久午夜电影| 亚洲精品国产成人久久av| 国产男人的电影天堂91| 欧美激情在线99| 成年女人永久免费观看视频| 搡女人真爽免费视频火全软件 | 精品久久久久久久久久久久久| 亚洲国产色片| 久久久色成人| 最后的刺客免费高清国语| 国产精品一区www在线观看| 一本久久中文字幕| 成人特级av手机在线观看| 99久久精品热视频| 免费电影在线观看免费观看| 久久久久久九九精品二区国产| 伊人久久精品亚洲午夜| 蜜桃亚洲精品一区二区三区| 午夜精品在线福利| av在线播放精品| a级毛片a级免费在线| 亚洲精品国产成人久久av| 免费不卡的大黄色大毛片视频在线观看 | 毛片女人毛片| 久久综合国产亚洲精品| 中文字幕熟女人妻在线| 嫩草影院入口| 日韩成人av中文字幕在线观看 | 国内揄拍国产精品人妻在线| 男女啪啪激烈高潮av片| 欧美一区二区精品小视频在线| ponron亚洲| 少妇人妻一区二区三区视频| 99久久九九国产精品国产免费| 亚洲中文字幕日韩| 国产一区二区在线观看日韩| 精华霜和精华液先用哪个| 精品久久国产蜜桃| 少妇猛男粗大的猛烈进出视频 | 欧美激情国产日韩精品一区| 成人一区二区视频在线观看| 性欧美人与动物交配| 中文字幕免费在线视频6| 美女cb高潮喷水在线观看| 激情 狠狠 欧美| 日韩三级伦理在线观看| 国产黄色小视频在线观看| 久久韩国三级中文字幕| 欧美一区二区亚洲| 午夜爱爱视频在线播放| 性欧美人与动物交配| 国产乱人偷精品视频| 亚洲,欧美,日韩| 亚洲色图av天堂| 91精品国产九色| 日韩高清综合在线| 精品一区二区三区视频在线| 亚洲在线观看片| 国产精品99久久久久久久久| 日韩精品中文字幕看吧| 3wmmmm亚洲av在线观看| 色综合站精品国产| 最近的中文字幕免费完整| 久久99热6这里只有精品| 亚洲成人久久爱视频| 国产精品国产高清国产av| 亚洲av五月六月丁香网| 国产精品国产高清国产av| 亚洲欧美日韩高清专用| 最近2019中文字幕mv第一页| 三级毛片av免费| 成人一区二区视频在线观看| 国产成人一区二区在线| 日日啪夜夜撸| 美女大奶头视频| 亚洲一级一片aⅴ在线观看| 欧美xxxx黑人xx丫x性爽| 午夜激情欧美在线| 人人妻人人看人人澡| 亚洲经典国产精华液单| 看非洲黑人一级黄片| a级毛片a级免费在线| 国产69精品久久久久777片| 亚洲丝袜综合中文字幕| 国产三级中文精品| 婷婷六月久久综合丁香| 蜜臀久久99精品久久宅男| 亚洲国产欧美人成| 高清毛片免费观看视频网站| 日本一本二区三区精品| 亚洲av五月六月丁香网| 嫩草影视91久久| 久久精品综合一区二区三区| 亚洲精品久久国产高清桃花| 91精品国产九色| 中文字幕久久专区| 国内精品美女久久久久久| 久久久久久久久大av| 插阴视频在线观看视频| 热99re8久久精品国产| 蜜桃久久精品国产亚洲av| 色综合色国产| 一级a爱片免费观看的视频| 久久久久久久久久黄片| 国产亚洲精品久久久com| 成人鲁丝片一二三区免费| 成人永久免费在线观看视频| 国产亚洲精品久久久久久毛片| 日本五十路高清| 狠狠狠狠99中文字幕| 最近的中文字幕免费完整| 欧美3d第一页| 免费无遮挡裸体视频| 18禁黄网站禁片免费观看直播| 一夜夜www| 99热全是精品| 高清毛片免费看| 永久网站在线| 国产精品国产三级国产av玫瑰| 在线免费十八禁| 观看免费一级毛片| 内射极品少妇av片p| 成人国产麻豆网| 日韩欧美免费精品| 免费黄网站久久成人精品| 国产探花极品一区二区| 免费av毛片视频| 精品免费久久久久久久清纯| 1024手机看黄色片| 天堂av国产一区二区熟女人妻| h日本视频在线播放| 长腿黑丝高跟| 久久午夜亚洲精品久久| 搡老熟女国产l中国老女人| 国产一级毛片七仙女欲春2| 国产精品免费一区二区三区在线| 一本一本综合久久| 国产又黄又爽又无遮挡在线| 在线观看av片永久免费下载| 成人一区二区视频在线观看| 插逼视频在线观看| 久久人人精品亚洲av| videossex国产| 一卡2卡三卡四卡精品乱码亚洲| 色哟哟哟哟哟哟| 国产欧美日韩精品亚洲av| 老司机午夜福利在线观看视频| 国产精品99久久久久久久久| 2021天堂中文幕一二区在线观| 久久久久国内视频| 一卡2卡三卡四卡精品乱码亚洲| 亚洲欧美精品自产自拍| 国产成人91sexporn| 国产日本99.免费观看| 久久久精品大字幕| av福利片在线观看| 久久精品综合一区二区三区| 99riav亚洲国产免费| 亚洲性夜色夜夜综合| 久久久久久久久大av| 国产女主播在线喷水免费视频网站 | 国产熟女欧美一区二区| 精品人妻视频免费看| 久久午夜亚洲精品久久| 直男gayav资源| 国产免费一级a男人的天堂| 亚洲自偷自拍三级| 免费看日本二区| 日本黄大片高清| 国产私拍福利视频在线观看| АⅤ资源中文在线天堂| 国产片特级美女逼逼视频| 国产亚洲精品综合一区在线观看| 一区二区三区高清视频在线| 日韩强制内射视频| 久久人妻av系列| 老熟妇仑乱视频hdxx| 一个人观看的视频www高清免费观看| 日韩欧美国产在线观看| 日日啪夜夜撸| 美女xxoo啪啪120秒动态图| 精品久久久久久成人av| 精品少妇黑人巨大在线播放 | 搡女人真爽免费视频火全软件 | 成人美女网站在线观看视频| 五月伊人婷婷丁香| 国产人妻一区二区三区在| 免费av不卡在线播放| 国产高清激情床上av| 97超级碰碰碰精品色视频在线观看| 亚洲精华国产精华液的使用体验 | 精品福利观看| 久久人人精品亚洲av| 简卡轻食公司| 国产欧美日韩精品一区二区| 久久精品夜夜夜夜夜久久蜜豆| 在线观看午夜福利视频| 国产免费男女视频| 亚洲精品日韩av片在线观看| 国产黄色视频一区二区在线观看 | 久久久久免费精品人妻一区二区| 久久热精品热| 亚洲欧美清纯卡通| 久久久精品大字幕| avwww免费| 毛片女人毛片| 国产毛片a区久久久久| 久久99热这里只有精品18| 干丝袜人妻中文字幕| 搞女人的毛片| 色综合站精品国产| 欧美极品一区二区三区四区| 国产男靠女视频免费网站| 色播亚洲综合网| 国产午夜精品论理片| 免费不卡的大黄色大毛片视频在线观看 | 久久久色成人| 人人妻人人澡人人爽人人夜夜 | 精品午夜福利视频在线观看一区| 午夜视频国产福利| 最好的美女福利视频网| 国产私拍福利视频在线观看| 久久国内精品自在自线图片| 看十八女毛片水多多多| 美女xxoo啪啪120秒动态图| 啦啦啦啦在线视频资源| 三级国产精品欧美在线观看| 俄罗斯特黄特色一大片| 国产伦精品一区二区三区四那| 变态另类丝袜制服| 老女人水多毛片| 草草在线视频免费看| 一个人观看的视频www高清免费观看| 18+在线观看网站| 婷婷色综合大香蕉| 成年免费大片在线观看| 精品福利观看| 日日干狠狠操夜夜爽| 中文字幕久久专区| 国产精品免费一区二区三区在线| 非洲黑人性xxxx精品又粗又长| 麻豆久久精品国产亚洲av| 淫秽高清视频在线观看| 国产色爽女视频免费观看| 国产一级毛片七仙女欲春2| 欧美成人一区二区免费高清观看| 啦啦啦韩国在线观看视频| 日韩欧美三级三区| 97超级碰碰碰精品色视频在线观看| 麻豆乱淫一区二区| 在线观看av片永久免费下载| 欧美色视频一区免费| 亚洲国产日韩欧美精品在线观看| 色哟哟·www| 99久久精品一区二区三区| 日本黄大片高清| 国产黄a三级三级三级人| 晚上一个人看的免费电影| 日韩中字成人| 欧美成人一区二区免费高清观看| 久久久久性生活片| 我的老师免费观看完整版| 欧美一区二区亚洲| 97超碰精品成人国产| 国产一区二区三区av在线 | 国产伦在线观看视频一区| 亚洲欧美成人精品一区二区| 我的女老师完整版在线观看| av在线老鸭窝| 婷婷六月久久综合丁香| 在线免费观看不下载黄p国产| 午夜福利高清视频| av卡一久久| 不卡视频在线观看欧美| 色综合站精品国产| 亚洲中文日韩欧美视频| 成人国产麻豆网|