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

    型砂壓實(shí)過程的非線性有限元數(shù)值模擬

    2014-06-06 07:17:54李文珍吳浚郊周克林
    鑄造設(shè)備與工藝 2014年2期
    關(guān)鍵詞:型砂砂箱砂型

    李文珍,吳浚郊,周克林

    (1.清華大學(xué)材料學(xué)院,北京 100084;2.上海汽車集團(tuán),上海 200041)

    在濕型砂鑄造生產(chǎn)中,為了保證鑄件的幾何尺寸和力學(xué)性能符合設(shè)計(jì)和使用要求,高質(zhì)量砂型的生產(chǎn)起決定作用。型砂緊實(shí)過程是生產(chǎn)砂型的關(guān)鍵工藝環(huán)節(jié),利用現(xiàn)代模擬仿真技術(shù)對型砂緊實(shí)過程進(jìn)行模擬,預(yù)測砂型緊實(shí)度分布并進(jìn)而預(yù)測砂型質(zhì)量,是提高造型生產(chǎn)效率獲得高質(zhì)量鑄型的有效手段。但是黏土濕型砂屬于離散材料體系,與一般的連續(xù)材料體系相比具有如下特點(diǎn):

    1)材料非線性即砂型由大量非連續(xù)的外層包覆黏土-水膜的細(xì)小砂粒組成,砂粒間相互接觸;

    2)幾何非線性即型砂在緊實(shí)過程中砂型的變形量很大,可達(dá)30%以上;

    3)接觸非線性即在緊實(shí)過程中型砂與砂箱和模樣之間存在非線性接觸。因此在型砂緊實(shí)過程中同時(shí)存在著材料非線性、幾何非線性和接觸非線性,這類問題一般用非線性有限元法進(jìn)行求解,如圖1所示。目前所能查到的濕砂造型數(shù)值模擬領(lǐng)域最早的研究是Leone和R.L.Lewis[1]使用有限元法對壓實(shí)造型中型砂密實(shí)過程進(jìn)行的數(shù)值模擬。

    1 用ANSYS軟件求解非線性問題方法簡介

    ANSYS軟件是一個(gè)商品化的大型通用有限元分析軟件,目前廣泛應(yīng)用于土木工程、水利工程、機(jī)械工程、航空航天和汽車交通等領(lǐng)域。ANSYS具有一些突出的技術(shù)特點(diǎn):將前處理模型的信息(單元、節(jié)點(diǎn)、材料),邊界條件(載荷、約束)以及后處理信息(求解結(jié)果)集成在統(tǒng)一的數(shù)據(jù)庫中,極大地方便了用戶的使用。ANSYS具有強(qiáng)大的非線性分析功能,可以處理幾何非線性、材料非線性、狀態(tài)非線性及單元非線性問題。ANSYS提供了100余種非線性材料模型如橡膠、泡沫、土壤和巖石等。另外ANSYS軟件能夠靈活地在前處理中進(jìn)行網(wǎng)格劃分,確保了單元的形態(tài)和精度。用ANSYS軟件求解非線性問題需要一系列帶校正的線性近似來實(shí)現(xiàn),如圖2.

    圖1 型砂壓實(shí)過程中非線性問題的求解流程

    圖2 純粹增量近似與牛頓-拉普森近似的關(guān)系[2]

    ANSYS軟件使用牛頓-拉普森平衡迭代進(jìn)行計(jì)算,它迫使在每一個(gè)載荷增量的末端解達(dá)到平衡收斂(在某個(gè)容限范圍內(nèi))。圖2b)描述了在單自由度非線性分析中牛頓-拉普森平衡迭代的使用。在每次求解前,NR方法估算出殘差矢量,這個(gè)矢量是回復(fù)力(對應(yīng)于單元應(yīng)力的載荷)和所加載荷的差值。程序然后使用非平衡載荷進(jìn)行線性求解,且核查收斂性。如果不滿足收斂準(zhǔn)則,重新估算非平衡載荷,修改剛度矩陣,獲得新解。持續(xù)這種迭代過程直到問題收斂[2]。

    ANSYS程序提供了一系列命令來增強(qiáng)問題的收斂性,如自適應(yīng)下降,線性搜索,自動(dòng)載荷步,及二分等,可被激活來加強(qiáng)問題的收斂性,如果不能得到收斂,那么程序依據(jù)輸入的命令或者繼續(xù)計(jì)算下一個(gè)載荷前或者終止[2]。

    2 型砂材料本構(gòu)關(guān)系建立

    型砂本構(gòu)關(guān)系即應(yīng)力應(yīng)變關(guān)系的建立是對砂型緊實(shí)過程進(jìn)行有限元分析的基礎(chǔ)。在土力學(xué)中已提出了很多測定和建立砂土本構(gòu)關(guān)系的方法,主要實(shí)驗(yàn)方法為有側(cè)限壓縮,它是目前最常見的用來測量土的壓縮性能的試驗(yàn)方法[3]。

    2.1 型砂本構(gòu)關(guān)系測試原理[4]

    圖3所示為型砂靜態(tài)本構(gòu)關(guān)系測試原理。具體方法是:將混好的型砂加入試樣筒中,隨后利用液壓式萬能試驗(yàn)機(jī)進(jìn)行慢速加載,通過試樣上的壓板作用在松散的型砂表面。在加載過程中,壓桿上的BLR-1載荷傳感器測量壓板載荷P的變化,位移傳感器測量壓板相應(yīng)的位移值U。這兩路傳感器信號經(jīng)Y6DL-1型應(yīng)變儀放大后,由LZ3-304函數(shù)記錄儀以P為縱坐標(biāo),U為橫坐標(biāo)輸出P-U曲線。經(jīng)下列轉(zhuǎn)換得到在靜態(tài)載荷作用下,有側(cè)限情況下的型砂應(yīng)力-應(yīng)變(本構(gòu)關(guān)系)曲線。

    式中:S-壓板面積;

    h-試樣筒的高度。

    圖3 測量型砂靜態(tài)本構(gòu)關(guān)系實(shí)驗(yàn)原理圖[4]

    2.2 試驗(yàn)結(jié)果及分析

    取型砂緊實(shí)率40%,膨潤土含量為10%,初始密度分別為 0.9g/cm3、1.0g/cm3、1.1 g/cm3、和 1.2g/cm3,測得的型砂應(yīng)力應(yīng)變關(guān)系如圖4所示。

    從圖4可以看出型砂應(yīng)力應(yīng)變關(guān)系為非線性,呈下凹形狀,這點(diǎn)與大多數(shù)連續(xù)材料的特性不同。未緊實(shí)的松散型砂在整個(gè)壓縮過程中具有非常大的變形量。其原因如下:型砂在開始受壓時(shí)的變形主要表現(xiàn)為砂粒重新排列填補(bǔ)空隙,因此應(yīng)力較??;當(dāng)型砂緊實(shí)到一定階段,砂粒之間的接觸變得更加緊密,砂粒運(yùn)動(dòng)不僅需要克服黏土膜之間的黏滯力,還要克服砂粒之間的摩擦力;在緊實(shí)的最后階段,砂粒之間石英內(nèi)核發(fā)生相互作用的正應(yīng)力,因?yàn)槭㈩w粒非常堅(jiān)硬,所以緊實(shí)后期應(yīng)力急劇增大。

    圖4 型砂本構(gòu)關(guān)系曲線

    3 型砂壓實(shí)過程的數(shù)值模擬

    在獲得型砂側(cè)限條件的本構(gòu)關(guān)系后,本文在有限元軟件ANSYS的平臺上,采用APDL進(jìn)行二次開發(fā),編寫了相應(yīng)的計(jì)算程序并耦合到現(xiàn)有ANSYS軟件中,對型砂壓實(shí)過程進(jìn)行了數(shù)值模擬。

    3.1 材料模型的建立

    從上述型砂本構(gòu)關(guān)系試驗(yàn)結(jié)果可以看出,型砂的本構(gòu)關(guān)系是一種非線性應(yīng)力應(yīng)變關(guān)系,其應(yīng)力對應(yīng)變的一階導(dǎo)數(shù)即彈性模量在變形過程中隨著應(yīng)變的增加而遞增,對于這樣的模型,通常的商業(yè)分析軟件無法直接予以處理,必須對實(shí)際本構(gòu)關(guān)系進(jìn)行一定的簡化處理,并采用APDL語言進(jìn)行編程來完成計(jì)算,以實(shí)現(xiàn)對這種復(fù)雜應(yīng)力應(yīng)變關(guān)系的求解。

    本文根據(jù)型砂材料的實(shí)際特點(diǎn),對型砂材料模型作了如下假定:

    1)型砂是均勻、連續(xù)、各向同性的介質(zhì);

    2)由于型砂緊實(shí)過程相對較快,因此,計(jì)算過程中不考慮回彈;

    3)應(yīng)力加載過程為從零單調(diào)增加到最大,卸載過程為無位移(很小)回彈。根據(jù)這些假定,對試驗(yàn)測得的型砂非線性本構(gòu)關(guān)系(圖4)采用多線性本構(gòu)關(guān)系進(jìn)行簡化,如圖5所示。在曲線上,每兩點(diǎn)之間為一直線段,彈性模量ΔE=Δσ/Δε為常數(shù),型砂表現(xiàn)為線彈性材料。本文選取了初始密度為0.9 g/cm3,膨潤土含量為8%,緊實(shí)率為40%的型砂作為計(jì)算對象,線性段階躍點(diǎn)的應(yīng)力應(yīng)變值如表1所示。

    圖5 型砂的多線性本構(gòu)關(guān)系

    表1 初始密度為0.9 g/cm3的型砂各線性段應(yīng)力應(yīng)變

    3.2 型砂幾何模型

    計(jì)算中,型砂幾何模型如圖6所示。設(shè)砂箱和余砂框的總高度為A,壓實(shí)后型砂頂面距離分型面的高度為B,未壓實(shí)前成型壓頭距離模樣頂面的距離為a,壓實(shí)后型砂定面距離模樣頂面的距離為b.為了使型砂壓實(shí)后型砂更加均勻,則成型壓頭尺寸滿足以下等式條件:

    圖6 幾何模型尺寸

    3.3 加載過程和邊界條件的處理

    在型砂壓實(shí)過程中,型砂與砂箱之間存在很大的摩擦力,這對型砂壓實(shí)過程有著不可忽略的影響,因此在計(jì)算過程中,在砂箱和型砂之間采用接觸條件,并參考文獻(xiàn)[3]中的數(shù)據(jù),取平均摩擦系數(shù)為0.5.

    由于型砂和砂箱之間存在摩擦力,型砂采用單面壓實(shí)后,緊實(shí)度分布情況一般為上緊下松,這對于處在分型面處的型腔在充填金屬液體時(shí)是極為不利的。因此,本文分別計(jì)算了兩種加載情況:

    1)采用上加壓單面壓實(shí);

    2)采用上下加壓雙面壓實(shí),且按照上面加載量為下面加載量的3倍進(jìn)行加載,并比較了兩種加載情況的計(jì)算結(jié)果。本文計(jì)算了緊實(shí)率為40%,膨潤土含量為8%,初始密度為0.9 g/cm3的型砂。采用位移加載,且按照壓縮40%的位移量加載;由于型砂的本構(gòu)關(guān)系特殊,將型砂分為幾個(gè)線性段進(jìn)行近似模擬,采用分步加載。在加載過程中,根據(jù)每個(gè)線性段的彈性應(yīng)變量進(jìn)行分步位移加載。當(dāng)前一個(gè)分步加載求解完畢后,提取計(jì)算結(jié)果并改變模型中各個(gè)節(jié)點(diǎn)的位置,重新劃分網(wǎng)格,然后改變彈性模量,在進(jìn)行下一個(gè)分步加載求解時(shí),將前一個(gè)分步計(jì)算的結(jié)果以預(yù)應(yīng)力載荷加在重新劃分網(wǎng)格的模型上,重復(fù)該計(jì)算過程,直到計(jì)算求解到最后一個(gè)線性段結(jié)束為止。

    單元采用8節(jié)點(diǎn)平面等參元,初始網(wǎng)格劃分結(jié)果如圖7所示。由于型砂在壓實(shí)過程中具有流動(dòng)性特征,若以彈性材料進(jìn)行求解,用ANSYS求解后會(huì)在成型壓頭附近出現(xiàn)一個(gè)缺口,為了解決此問題,將成型壓頭下面的節(jié)點(diǎn)在Y方向的位移進(jìn)行耦合求解。此外,成型壓頭兩邊的型砂變形量大,單元若太小,在計(jì)算過程中容易發(fā)生不收斂而導(dǎo)致計(jì)算中止,因此特意將成型壓頭兩邊的型砂單元?jiǎng)澐值妮^大,由于我們感興趣的區(qū)域?yàn)榉中兔娓浇膽?yīng)力分布,因此在型砂頂面將單元?jiǎng)澊髮Y(jié)果不會(huì)有很大的影響[5]。

    圖7 模型網(wǎng)格劃分結(jié)果

    4 數(shù)值模擬結(jié)果及分析

    4.1 單面壓實(shí)和雙面壓實(shí)數(shù)值模擬結(jié)果的比較

    數(shù)值模擬結(jié)果如圖8a)和圖b)分別是單面壓實(shí)和雙面壓實(shí)后的計(jì)算結(jié)果。為了比較兩種加載情況的計(jì)算結(jié)果,每一種加載情況都提取了兩條線上的應(yīng)力分布(如圖6):一條靠近砂箱邊緣,一條則為模型的對稱中心線,兩條線方向均由下向上。兩條線上的應(yīng)力值如表2和表3.將表中的應(yīng)力制成圖,如圖9所示。

    表2 沿砂箱邊緣高度方向的應(yīng)力分布各點(diǎn)的值

    表3 沿砂箱中心線高度方向的應(yīng)力分布各點(diǎn)的值

    圖8 數(shù)值模擬應(yīng)力場的分布

    圖9 沿砂箱邊緣和模型中心兩條線上的應(yīng)力分布

    從圖8和圖9的模擬計(jì)算結(jié)果可以看出,由于型砂砂粒之間摩擦力及型砂與砂箱壁和模型表面之間摩擦力的影響,在靠近砂箱的附近應(yīng)力沿砂箱高度方向增加(如圖9),而在砂箱中心線附近應(yīng)力則是沿砂箱高度降低(如圖9),模擬結(jié)果符合實(shí)際緊實(shí)情況。同時(shí)值得注意的是,在模樣拐角處和型砂狹窄區(qū)域入口處,應(yīng)力場為顏色較深,表明在該處出現(xiàn)一個(gè)低應(yīng)力區(qū)域,應(yīng)力低于其上方區(qū)域的應(yīng)力,模擬出了文獻(xiàn)[4]試驗(yàn)中的“搭橋”現(xiàn)象。

    此外從圖9還可以看出,在高度較低的地方,雙面壓實(shí)砂箱邊緣應(yīng)力曲線的值要比單面壓實(shí)砂箱邊緣應(yīng)力曲線的值要大;而在高度較高的地方,雙面壓實(shí)砂箱邊緣應(yīng)力曲線的值反而要比單面壓實(shí)砂箱邊緣應(yīng)力曲線的值要小。由此說明,雙面壓實(shí)后,可以提高分型面處型砂的緊實(shí)度,且整個(gè)型砂應(yīng)力分布要更加均勻,符合雙面壓實(shí)要優(yōu)于單面壓實(shí)的規(guī)律。

    4.2 數(shù)值模擬結(jié)果和實(shí)驗(yàn)測量值的比較

    由于實(shí)驗(yàn)條件所限,本文參考了文獻(xiàn)[4]中的試驗(yàn)數(shù)據(jù)(如圖10所示,圖中應(yīng)力單位為MPa),并將模擬結(jié)果與之作了數(shù)量級上的比較以驗(yàn)證模擬結(jié)果。

    圖10 參考文獻(xiàn)[4]的實(shí)測應(yīng)力分布圖

    文獻(xiàn)[4]中的試驗(yàn)條件如下:砂箱平面尺寸為390 mm×350 mm,型砂初始密度為0.9 g/cm3,緊實(shí)率40%,膨潤土含量為8%。這些試驗(yàn)條件除了砂箱尺寸和本文計(jì)算模型有一定的差別(本文模型尺寸為400 mm×300 mm),其他條件均相等,因而有一定的參考意義。

    圖9中計(jì)算出的沿砂箱高度方向的應(yīng)力曲線和圖10b)實(shí)測應(yīng)力曲線非常相似,趨勢相同,應(yīng)力均隨著高度的增加而增加,且都在砂箱和模樣之間存在一塊應(yīng)力比較小的區(qū)域,也即“拱橋效應(yīng)”。在應(yīng)力值的數(shù)量級上,在“拱橋效應(yīng)”區(qū)域,計(jì)算出的應(yīng)力值大約在0.34 MPa~0.9 MPa之間,實(shí)測的應(yīng)力值在0.36 MPa~0.7MPa之間;而在高出模樣的區(qū)域,計(jì)算出的應(yīng)力值大約在0.9 MPa~1.9 MPa之間,而試驗(yàn)測得的值也在0.8 MPa以上。此外,從計(jì)算結(jié)果和試驗(yàn)結(jié)果可以看出,在模樣上方都有一塊比較緊實(shí)的區(qū)域。上述結(jié)果表明數(shù)值模擬的結(jié)果和試驗(yàn)結(jié)果比較吻合。

    5 結(jié)論

    利用商品化有限元分析軟件ANSYS,通過APDL二次開發(fā)編程對黏土砂造型的壓實(shí)過程進(jìn)行了數(shù)值模擬。通過對兩種加載方法的計(jì)算結(jié)果比較以及計(jì)算結(jié)果與相應(yīng)的試驗(yàn)數(shù)據(jù)比較,結(jié)果表明ANSYS可以很好地解決型砂緊實(shí)這樣的非線性問題。濕型砂壓實(shí)過程的數(shù)值模擬結(jié)果可以較好地解釋如下現(xiàn)象:

    1)“拱橋效應(yīng)”現(xiàn)象;

    2)“雙面加壓”要優(yōu)于“單面加壓”,從而可以提高型砂分型面的強(qiáng)度;

    3)在砂箱邊緣應(yīng)力分布從下至上逐漸變大。

    通過ANSYS軟件對型砂緊實(shí)過程進(jìn)行數(shù)值模擬取得了較好的效果。但是,由于實(shí)驗(yàn)數(shù)據(jù)來自參考文獻(xiàn),因而以后需要通過進(jìn)一步的實(shí)驗(yàn)來驗(yàn)證模擬結(jié)果的精度。此外,本文模擬仍局限在二維模型上,對于三維模型模擬有待進(jìn)一步的開發(fā)。

    [1]Leone M,Lewis R L.Numerical Modeling of Green Sand During Compaction[J].AFSTransactions,1988,136:763-774.

    [2]王國強(qiáng).實(shí)用工程數(shù)值模擬技術(shù)及其在ANSYS上的實(shí)踐[M].西安:西北工業(yè)大學(xué)出版社,1999.

    [3]張翼飛.射壓造型的數(shù)值模擬和試驗(yàn)研究[D].北京:清華大學(xué)機(jī)械工程系.2001.

    [4]謝濱.壓實(shí)造型和空氣沖擊造型砂型應(yīng)力場研究[D].北京:清華大學(xué)機(jī)械工程系,1993.

    [5]Li Wenzhen,Wu Junjiao.Numerical Simulation of Compacting Process of Green Sand Molding Based on Sand Filling[J],Materials Science Forum,2007,561-565:1879-1882.

    猜你喜歡
    型砂砂箱砂型
    淺談砂箱斷面尺寸對砂箱剛性的影響
    基于Ansys workbench 的單壁砂箱有限元分析
    鋼軌鋁熱焊砂型硬化工藝試驗(yàn)研究和失效分析
    型砂壓力傳遞在復(fù)合材料成型中的應(yīng)用研究
    3DP 砂型打印在鑄造生產(chǎn)中的應(yīng)用研究
    一種用于熔鑄磚的新型水玻璃砂型
    玻璃(2020年2期)2020-12-11 22:20:04
    一種嵌裝式砂型
    玻璃(2020年2期)2020-12-11 22:20:04
    鑄造粉塵回收技術(shù)探討
    砂處理型砂水分調(diào)節(jié)與影響分析
    砂型鑄造模具設(shè)計(jì)及制造
    爱豆传媒免费全集在线观看| 国产男女内射视频| 精品国产一区二区三区久久久樱花 | h视频一区二区三区| 国产视频内射| 啦啦啦中文免费视频观看日本| 国国产精品蜜臀av免费| 亚洲天堂av无毛| 国产免费一级a男人的天堂| 人体艺术视频欧美日本| 国产毛片在线视频| 亚洲图色成人| 最近2019中文字幕mv第一页| 一二三四中文在线观看免费高清| 国产精品国产三级国产av玫瑰| 亚洲精品久久久久久婷婷小说| 国产熟女欧美一区二区| 日本黄色日本黄色录像| 99久久精品一区二区三区| 国产日韩欧美亚洲二区| 少妇熟女欧美另类| 国产片特级美女逼逼视频| 91久久精品国产一区二区成人| 午夜激情久久久久久久| 欧美亚洲 丝袜 人妻 在线| 我的老师免费观看完整版| 欧美最新免费一区二区三区| 人人妻人人爽人人添夜夜欢视频 | 成人漫画全彩无遮挡| 国产亚洲91精品色在线| 自拍欧美九色日韩亚洲蝌蚪91 | 成人国产麻豆网| 免费观看a级毛片全部| 日产精品乱码卡一卡2卡三| 亚洲精品亚洲一区二区| 制服丝袜香蕉在线| 色网站视频免费| 五月开心婷婷网| 国产毛片在线视频| 亚洲精华国产精华液的使用体验| 日韩视频在线欧美| 又爽又黄a免费视频| 精品亚洲成a人片在线观看 | 在线看a的网站| av福利片在线观看| 人人妻人人爽人人添夜夜欢视频 | 欧美精品亚洲一区二区| 亚洲av免费高清在线观看| 一区二区三区免费毛片| 啦啦啦中文免费视频观看日本| 在线观看国产h片| 亚洲精品国产成人久久av| 丰满乱子伦码专区| 色哟哟·www| 欧美国产精品一级二级三级 | 色婷婷av一区二区三区视频| 国产成人免费观看mmmm| 3wmmmm亚洲av在线观看| 少妇丰满av| 亚洲成人中文字幕在线播放| 国产精品无大码| 亚洲精品亚洲一区二区| av在线观看视频网站免费| 国产精品.久久久| 国产伦理片在线播放av一区| 国产伦精品一区二区三区四那| 国产 一区 欧美 日韩| 欧美高清性xxxxhd video| 内地一区二区视频在线| 蜜桃久久精品国产亚洲av| 特大巨黑吊av在线直播| 日本黄色片子视频| 在线免费观看不下载黄p国产| av女优亚洲男人天堂| 五月天丁香电影| 熟女av电影| 久久精品国产亚洲网站| 丝袜脚勾引网站| 性色avwww在线观看| 18禁在线无遮挡免费观看视频| av卡一久久| 国产高清有码在线观看视频| 啦啦啦在线观看免费高清www| 国产精品人妻久久久影院| 18禁裸乳无遮挡动漫免费视频| 18禁裸乳无遮挡动漫免费视频| 我的女老师完整版在线观看| 日本vs欧美在线观看视频 | 亚洲第一区二区三区不卡| 高清午夜精品一区二区三区| 国产伦在线观看视频一区| 亚洲精品456在线播放app| 精品人妻熟女av久视频| av国产久精品久网站免费入址| 亚洲内射少妇av| 男女下面进入的视频免费午夜| 国产女主播在线喷水免费视频网站| 成人毛片a级毛片在线播放| 日韩欧美 国产精品| 欧美精品一区二区大全| 亚洲精品乱码久久久久久按摩| 亚洲av欧美aⅴ国产| 欧美日韩综合久久久久久| 一级毛片黄色毛片免费观看视频| 日本免费在线观看一区| 亚洲成人av在线免费| 国产熟女欧美一区二区| 亚洲国产高清在线一区二区三| 交换朋友夫妻互换小说| 久久久亚洲精品成人影院| 国产黄色免费在线视频| 深夜a级毛片| 蜜臀久久99精品久久宅男| 亚洲激情五月婷婷啪啪| 日本免费在线观看一区| 午夜福利网站1000一区二区三区| 全区人妻精品视频| 久久久久久久久久久免费av| 插逼视频在线观看| 网址你懂的国产日韩在线| 三级经典国产精品| 看免费成人av毛片| 久久鲁丝午夜福利片| 波野结衣二区三区在线| 国产av精品麻豆| 性高湖久久久久久久久免费观看| 波野结衣二区三区在线| 欧美精品人与动牲交sv欧美| 日韩av在线免费看完整版不卡| 大话2 男鬼变身卡| 中文字幕av成人在线电影| 日本免费在线观看一区| 午夜福利网站1000一区二区三区| 日日啪夜夜撸| 这个男人来自地球电影免费观看 | 欧美日韩综合久久久久久| 九草在线视频观看| 国产精品人妻久久久久久| 女人久久www免费人成看片| 精品熟女少妇av免费看| 国产女主播在线喷水免费视频网站| 欧美精品人与动牲交sv欧美| 我的老师免费观看完整版| 精品人妻偷拍中文字幕| 亚洲av国产av综合av卡| 亚洲精品乱码久久久久久按摩| 成人毛片a级毛片在线播放| 成人毛片a级毛片在线播放| 成人毛片60女人毛片免费| av在线蜜桃| 午夜福利在线在线| 哪个播放器可以免费观看大片| 少妇人妻久久综合中文| 久久久久视频综合| 激情 狠狠 欧美| 91精品国产国语对白视频| 欧美成人一区二区免费高清观看| av国产精品久久久久影院| 久久韩国三级中文字幕| 国语对白做爰xxxⅹ性视频网站| 久久久色成人| 99热这里只有是精品在线观看| 国产爽快片一区二区三区| 国产v大片淫在线免费观看| 亚洲精品乱码久久久v下载方式| 国产高清三级在线| 黄色欧美视频在线观看| 日韩欧美一区视频在线观看 | 18禁动态无遮挡网站| 直男gayav资源| 午夜激情久久久久久久| 搡女人真爽免费视频火全软件| 亚洲怡红院男人天堂| 大片电影免费在线观看免费| 各种免费的搞黄视频| 在线观看美女被高潮喷水网站| 午夜老司机福利剧场| 少妇 在线观看| 久久久久久久久久久丰满| 在线播放无遮挡| 精品国产一区二区三区久久久樱花 | 纯流量卡能插随身wifi吗| 国产成人aa在线观看| 免费人妻精品一区二区三区视频| 国产精品免费大片| 国产精品国产三级国产专区5o| 色视频www国产| 精品午夜福利在线看| 免费久久久久久久精品成人欧美视频 | 18+在线观看网站| www.av在线官网国产| av.在线天堂| 国产午夜精品久久久久久一区二区三区| 日韩三级伦理在线观看| 少妇裸体淫交视频免费看高清| 好男人视频免费观看在线| 亚洲精品成人av观看孕妇| 精品亚洲乱码少妇综合久久| 欧美日本视频| 亚洲成人av在线免费| 日韩av不卡免费在线播放| 亚洲成人中文字幕在线播放| 大码成人一级视频| 啦啦啦视频在线资源免费观看| av国产久精品久网站免费入址| 免费在线观看成人毛片| 中文在线观看免费www的网站| 久久精品熟女亚洲av麻豆精品| 亚洲综合色惰| 最黄视频免费看| 久久久久久伊人网av| 成人漫画全彩无遮挡| 美女cb高潮喷水在线观看| 成年av动漫网址| 尾随美女入室| 亚洲经典国产精华液单| 在线看a的网站| 3wmmmm亚洲av在线观看| 免费不卡的大黄色大毛片视频在线观看| 男女国产视频网站| 国产精品久久久久成人av| 蜜桃亚洲精品一区二区三区| 日日摸夜夜添夜夜添av毛片| 久久久久久九九精品二区国产| 中国三级夫妇交换| 国产精品久久久久久av不卡| 天堂俺去俺来也www色官网| 赤兔流量卡办理| 国产高清不卡午夜福利| 亚洲av中文av极速乱| 亚洲欧洲国产日韩| 国产精品精品国产色婷婷| 老司机影院成人| 伊人久久国产一区二区| 亚洲成人av在线免费| 亚洲精品,欧美精品| 久久久久视频综合| 嫩草影院新地址| 国产 精品1| 青春草视频在线免费观看| 国产精品一区二区在线观看99| 边亲边吃奶的免费视频| 亚洲一级一片aⅴ在线观看| 91精品伊人久久大香线蕉| 国产精品伦人一区二区| 成人午夜精彩视频在线观看| 欧美区成人在线视频| 欧美日韩精品成人综合77777| 久久人人爽人人爽人人片va| 久久久久久久精品精品| 国产成人aa在线观看| 精品久久久噜噜| 插逼视频在线观看| 一区二区av电影网| 亚洲色图综合在线观看| 免费黄色在线免费观看| 日韩视频在线欧美| 全区人妻精品视频| 日本与韩国留学比较| 亚洲欧洲日产国产| 亚洲精品日本国产第一区| 高清毛片免费看| 亚洲性久久影院| 亚洲av中文av极速乱| 少妇的逼好多水| 午夜福利影视在线免费观看| 九九爱精品视频在线观看| 欧美+日韩+精品| 国产精品无大码| 99九九线精品视频在线观看视频| 狂野欧美激情性bbbbbb| 国产淫语在线视频| 亚洲第一区二区三区不卡| 丰满人妻一区二区三区视频av| 王馨瑶露胸无遮挡在线观看| 六月丁香七月| 国产中年淑女户外野战色| 成人美女网站在线观看视频| 三级经典国产精品| 久久久久精品久久久久真实原创| 成人影院久久| 51国产日韩欧美| 亚洲三级黄色毛片| 久久人妻熟女aⅴ| 亚洲精品乱久久久久久| 好男人视频免费观看在线| 亚洲国产成人一精品久久久| 久久精品国产亚洲av涩爱| 色综合色国产| 男人舔奶头视频| 亚洲欧美成人综合另类久久久| 七月丁香在线播放| videossex国产| 在线亚洲精品国产二区图片欧美 | 欧美一区二区亚洲| av女优亚洲男人天堂| 联通29元200g的流量卡| 亚洲,欧美,日韩| 老司机影院毛片| 精品人妻一区二区三区麻豆| 免费少妇av软件| 国产有黄有色有爽视频| 日韩av在线免费看完整版不卡| 少妇被粗大猛烈的视频| 大陆偷拍与自拍| 欧美日韩在线观看h| 日韩三级伦理在线观看| 天堂8中文在线网| 麻豆成人午夜福利视频| 亚洲精品456在线播放app| 麻豆成人av视频| 直男gayav资源| 91精品一卡2卡3卡4卡| 我的老师免费观看完整版| 自拍偷自拍亚洲精品老妇| 国产精品国产三级国产av玫瑰| 一级毛片我不卡| 国产精品伦人一区二区| 国产成人精品一,二区| 美女中出高潮动态图| 成年美女黄网站色视频大全免费 | 男人舔奶头视频| 国产亚洲最大av| 午夜福利影视在线免费观看| 欧美成人精品欧美一级黄| 在线播放无遮挡| 久久久久久九九精品二区国产| 亚洲一区二区三区欧美精品| 少妇被粗大猛烈的视频| 欧美zozozo另类| 亚洲欧美清纯卡通| 国产亚洲一区二区精品| 亚洲,欧美,日韩| 亚洲av欧美aⅴ国产| 观看av在线不卡| 美女xxoo啪啪120秒动态图| 一区二区三区免费毛片| 一本一本综合久久| 国产精品一区www在线观看| 欧美变态另类bdsm刘玥| 欧美日韩综合久久久久久| 国产黄频视频在线观看| 中国国产av一级| 亚洲欧美清纯卡通| 毛片女人毛片| 黄色欧美视频在线观看| 色网站视频免费| 国产亚洲5aaaaa淫片| 黄色怎么调成土黄色| 国产 精品1| 尤物成人国产欧美一区二区三区| 99精国产麻豆久久婷婷| 精品人妻视频免费看| 国产亚洲午夜精品一区二区久久| 午夜福利视频精品| 国产爽快片一区二区三区| av视频免费观看在线观看| 欧美xxⅹ黑人| 一级毛片电影观看| 国产精品嫩草影院av在线观看| 如何舔出高潮| 国产av码专区亚洲av| 久久久成人免费电影| 国产高清国产精品国产三级 | av视频免费观看在线观看| 国产成人精品久久久久久| 亚洲精品国产色婷婷电影| 国产老妇伦熟女老妇高清| tube8黄色片| 亚洲精品aⅴ在线观看| 色婷婷av一区二区三区视频| 天美传媒精品一区二区| 欧美极品一区二区三区四区| 熟女av电影| 我的老师免费观看完整版| 只有这里有精品99| 国产国拍精品亚洲av在线观看| 日日摸夜夜添夜夜爱| 中文字幕精品免费在线观看视频 | 国产一区二区三区av在线| 又大又黄又爽视频免费| 国产亚洲午夜精品一区二区久久| 一级毛片我不卡| 天堂俺去俺来也www色官网| 在线观看免费视频网站a站| 婷婷色麻豆天堂久久| 亚洲精品乱久久久久久| 三级国产精品片| 在线免费观看不下载黄p国产| 亚洲精品一二三| 只有这里有精品99| 一级a做视频免费观看| 三级国产精品片| 久久精品熟女亚洲av麻豆精品| 国产成人精品福利久久| 亚洲av.av天堂| 亚洲,一卡二卡三卡| 亚洲精品成人av观看孕妇| 国产探花极品一区二区| 午夜福利高清视频| 亚州av有码| 狂野欧美白嫩少妇大欣赏| 精品亚洲乱码少妇综合久久| 欧美激情极品国产一区二区三区 | 日本av手机在线免费观看| 国产黄色视频一区二区在线观看| 中文字幕人妻熟人妻熟丝袜美| 1000部很黄的大片| 我的女老师完整版在线观看| 国产精品三级大全| 国产人妻一区二区三区在| 亚洲美女搞黄在线观看| 国产成人freesex在线| 亚洲av国产av综合av卡| 久久 成人 亚洲| 嘟嘟电影网在线观看| 欧美老熟妇乱子伦牲交| 26uuu在线亚洲综合色| 国产综合精华液| 亚洲人与动物交配视频| 日本黄色日本黄色录像| 久久国产乱子免费精品| 国产精品欧美亚洲77777| av不卡在线播放| 精品亚洲成a人片在线观看 | 精品一区二区免费观看| av在线播放精品| 国产精品伦人一区二区| 看十八女毛片水多多多| 精品久久久久久久久亚洲| 深爱激情五月婷婷| 2021少妇久久久久久久久久久| 在线观看人妻少妇| 高清午夜精品一区二区三区| 只有这里有精品99| 欧美xxxx黑人xx丫x性爽| 尤物成人国产欧美一区二区三区| 美女脱内裤让男人舔精品视频| 夫妻午夜视频| 亚洲精品成人av观看孕妇| a级一级毛片免费在线观看| 精品久久久久久久久亚洲| 久久国产亚洲av麻豆专区| 女的被弄到高潮叫床怎么办| 久久国产精品大桥未久av | 人人妻人人看人人澡| 精品久久国产蜜桃| 偷拍熟女少妇极品色| 99久久精品国产国产毛片| 成人美女网站在线观看视频| 欧美日韩精品成人综合77777| 91精品伊人久久大香线蕉| 一本久久精品| 99国产精品免费福利视频| 免费播放大片免费观看视频在线观看| 啦啦啦中文免费视频观看日本| 欧美一级a爱片免费观看看| 人妻夜夜爽99麻豆av| 午夜日本视频在线| 一级片'在线观看视频| 简卡轻食公司| 国产在线男女| 成人毛片60女人毛片免费| 少妇猛男粗大的猛烈进出视频| 熟妇人妻不卡中文字幕| 伊人久久精品亚洲午夜| 五月玫瑰六月丁香| 乱系列少妇在线播放| 久久毛片免费看一区二区三区| 午夜日本视频在线| 91久久精品国产一区二区三区| 中文字幕免费在线视频6| av天堂中文字幕网| 国产精品三级大全| 新久久久久国产一级毛片| 亚洲色图av天堂| 搡女人真爽免费视频火全软件| 高清欧美精品videossex| 国产乱人偷精品视频| 国产精品一区二区在线观看99| 欧美性感艳星| 亚洲精品456在线播放app| 一个人看的www免费观看视频| 国产精品爽爽va在线观看网站| 秋霞在线观看毛片| 欧美高清成人免费视频www| 天天躁夜夜躁狠狠久久av| xxx大片免费视频| 亚洲va在线va天堂va国产| 久久久久久久久久成人| 精品久久久久久电影网| 91精品国产国语对白视频| 亚洲av中文字字幕乱码综合| 黄片无遮挡物在线观看| 美女cb高潮喷水在线观看| 夜夜爽夜夜爽视频| 男女啪啪激烈高潮av片| av不卡在线播放| 赤兔流量卡办理| 中文欧美无线码| 久久久久精品性色| 熟女av电影| 少妇人妻一区二区三区视频| 国产 精品1| 国产免费又黄又爽又色| 国产乱人视频| 亚洲精品日韩av片在线观看| av国产久精品久网站免费入址| 免费黄频网站在线观看国产| 五月伊人婷婷丁香| 一本久久精品| av视频免费观看在线观看| 日本wwww免费看| 国产精品熟女久久久久浪| 免费人成在线观看视频色| 搡女人真爽免费视频火全软件| 国产老妇伦熟女老妇高清| 日本免费在线观看一区| a级毛片免费高清观看在线播放| 亚洲av综合色区一区| 日产精品乱码卡一卡2卡三| 欧美成人a在线观看| 在线观看一区二区三区激情| 日本午夜av视频| 精品久久久精品久久久| 蜜桃在线观看..| 精品人妻视频免费看| 黄色怎么调成土黄色| 高清在线视频一区二区三区| av国产免费在线观看| 国内精品宾馆在线| 18+在线观看网站| 啦啦啦在线观看免费高清www| 久久久久性生活片| 伦理电影大哥的女人| 一级毛片久久久久久久久女| 中文字幕制服av| 久久久午夜欧美精品| 欧美zozozo另类| 久久精品国产鲁丝片午夜精品| 波野结衣二区三区在线| 在线观看三级黄色| 色视频www国产| 久久久亚洲精品成人影院| 久久青草综合色| 欧美成人精品欧美一级黄| 欧美精品一区二区大全| 久久99热这里只频精品6学生| a级毛色黄片| 免费观看av网站的网址| 成人18禁高潮啪啪吃奶动态图 | 日本黄色片子视频| 不卡视频在线观看欧美| 日本猛色少妇xxxxx猛交久久| 狂野欧美白嫩少妇大欣赏| 1000部很黄的大片| 性色avwww在线观看| 亚洲三级黄色毛片| 亚洲综合色惰| 精品99又大又爽又粗少妇毛片| 哪个播放器可以免费观看大片| 下体分泌物呈黄色| 身体一侧抽搐| 国内精品宾馆在线| 看十八女毛片水多多多| 欧美成人a在线观看| 亚洲在久久综合| 国产综合精华液| 成人美女网站在线观看视频| 日本免费在线观看一区| 中文字幕精品免费在线观看视频 | 久久久久久人妻| 两个人的视频大全免费| 丝袜脚勾引网站| 精品酒店卫生间| av国产精品久久久久影院| 赤兔流量卡办理| 欧美bdsm另类| 中文字幕免费在线视频6| 男的添女的下面高潮视频| 国产人妻一区二区三区在| 国产精品一区二区三区四区免费观看| 久久韩国三级中文字幕| 国产精品人妻久久久久久| 超碰av人人做人人爽久久| 精品一区二区三卡| 一区在线观看完整版| 亚洲国产成人一精品久久久| 男女啪啪激烈高潮av片| 麻豆乱淫一区二区| 国产免费福利视频在线观看| 中文乱码字字幕精品一区二区三区| 成人免费观看视频高清| 毛片一级片免费看久久久久| 91精品一卡2卡3卡4卡| 国产精品秋霞免费鲁丝片| 直男gayav资源| 91aial.com中文字幕在线观看| 婷婷色综合www| 亚洲性久久影院| 久久ye,这里只有精品| 久久精品国产a三级三级三级| 国产精品偷伦视频观看了| 蜜桃亚洲精品一区二区三区| 亚洲精品视频女| 少妇人妻精品综合一区二区| 国产av精品麻豆| 午夜激情久久久久久久| 欧美变态另类bdsm刘玥| www.色视频.com| 大又大粗又爽又黄少妇毛片口| 日韩大片免费观看网站| 日韩电影二区| 高清日韩中文字幕在线| 亚洲精品中文字幕在线视频 |