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

    基于有限元的天平校準系統(tǒng)基座結(jié)構(gòu)優(yōu)化設(shè)計

    2017-07-15 21:53:34趙永豐王均濤
    有色金屬材料與工程 2017年3期
    關(guān)鍵詞:拓撲優(yōu)化模態(tài)分析基座

    趙永豐+王均濤

    摘要: 在Solidworks環(huán)境下建立三維模型,運用有限元軟件對基座結(jié)構(gòu)進行動、靜態(tài)性能分析,得到初始模型的性能參數(shù).以最小質(zhì)量為目標函數(shù),以剛度、應力和頻率為約束條件,建立數(shù)學模型.提出一種基于漸進結(jié)構(gòu)優(yōu)化法的拓撲優(yōu)化方法,引入過濾函數(shù)來解決拓撲優(yōu)化過程中容易出現(xiàn)的數(shù)值不穩(wěn)定、數(shù)值分布呈棋盤模式的現(xiàn)象,保證計算結(jié)果的準確性.對拓撲優(yōu)化后的模型重新分析計算后發(fā)現(xiàn),其性能依然優(yōu)異,基頻提升26.7%,質(zhì)量減少5.8%.

    關(guān)鍵詞: 基座; 靜力分析; 模態(tài)分析; 拓撲優(yōu)化; 漸進結(jié)構(gòu)優(yōu)化

    中圖分類號: TH 715文獻標志碼: A

    Optimization Design of Pedestal Structure of Balance

    Calibration System Based on Finite Element Method

    ZHAO Yongfeng1,WANG Juntao2

    (1.School of Mechanical Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China;

    2.Shanghai Machine Tool Works Ltd., Shanghai 200093, China)

    Abstract: The threedimensional model is established in Solidworks,and the static and dynamic performance of the pedestal structure is analyzed by finite element software,which can help to get the performance parameters of the initial model.Based on the minimummass as the objective function,the mathematical model for topology optimization is established with the stiffness,stress and frequency as the constraint condition.For ensuring the accuracy of the calculation results,a method based on evolution structure optimization for topology optimization is presented and a filtering function is introduced in order to solve the numerical instability and checkerboard phenomenon in the topology optimization process.After analyzing the calculation of the topological optimized model,the structure performance is still excellent.As a result,the fundamental frequency increased by 27.6%,and the mass reduced by 5.8%.

    Keywords: pedestal; static analysis; modal analysis; topology optimization; evolution structure optimization

    風洞天平校準系統(tǒng)是目前世界各國在各類型氣動力測量試驗中廣泛采用的測量裝置,其性能指標直接影響天平校準系統(tǒng)的試驗測量數(shù)據(jù)精度.風洞天平校準裝置主要由龍門架、基座和四自由度機械臺體組成,其系統(tǒng)結(jié)構(gòu)如圖1所示.基座支撐龍門架、加載機構(gòu)和四自由度機械臺等部件,在整個系統(tǒng)中占總質(zhì)量的44.6%.因此,在保證剛度和強度的情況下,合理設(shè)計基座結(jié)構(gòu),進行拓撲優(yōu)化,以使基座質(zhì)量減小,同時降低工藝難度與生產(chǎn)成本,獲得良好的經(jīng)濟效益.基座內(nèi)部的筋板結(jié)構(gòu)和布局在很大程度上決定了基座的動、靜態(tài)性能[1],而基座靜態(tài)特性與動態(tài)特性又影響著天平校準系統(tǒng)的測量性能與精度穩(wěn)定性.低階模態(tài)特性基本決定了產(chǎn)品的動態(tài)性能,因此,在基座結(jié)構(gòu)優(yōu)化過程中也要考慮頻率方面的約束[2].對天平校準系統(tǒng)的基座進行靜力分析與模態(tài)分析,利用漸進結(jié)構(gòu)優(yōu)化法,以最小質(zhì)量為目標函數(shù),以剛度、應力和頻率為約束條件,引入過濾函數(shù)來解決拓撲優(yōu)化過程中容易出現(xiàn)的數(shù)值不穩(wěn)定現(xiàn)象.最后對基座進行拓撲分析,從滿足約束的結(jié)構(gòu)中刪除多余的材料,以此來獲得最小的基座質(zhì)量.

    1.1模型建立

    基座由灰鑄鐵鑄造而成,結(jié)構(gòu)復雜多變,主要由厚度不一、形狀各異的筋板,以及鑄造岀砂孔構(gòu)成,根據(jù)實際使用情況需要,整體結(jié)構(gòu)形狀設(shè)計為T型.Solidworks環(huán)境下建立的三維模型是有限元分析和結(jié)構(gòu)優(yōu)化的前提和基礎(chǔ).為了便于有限元分析計算,必須對三維模型進行結(jié)構(gòu)簡化,省去倒角、圓角、螺紋孔等細微特征.簡化后的模型及有限元模型如圖2、圖3所示.基座長3 100 mm,兩翼寬2 990 mm,中間寬1 260 mm,高520 mm,設(shè)定中間筋板厚度25 mm,兩翼筋板厚度20 mm,壁厚25 mm.

    1.2靜力分析

    靜力分析用于計算在固定不變的載荷作用下結(jié)構(gòu)的效應,它不考慮慣性和阻尼的影響,如結(jié)構(gòu)隨時間變化的載荷的情況.其靜力學方程為:

    [K]{u}={F}(1)

    式中:[K]為剛度矩陣;{u}為各節(jié)點位移;{F}為靜力載荷.

    在靜態(tài)工作環(huán)境下,基座主要受到3個相同的側(cè)力加載機構(gòu)的重力,其值均為9 000 N,四自由度機械臺重力約為30 000 N,龍門架及其上的測量架與升力加載機構(gòu)的重力約為20 000 N.在ANSYS軟件中設(shè)置材料、劃分網(wǎng)格、施加約束和作用力,得到變形與應力分布圖,如圖4、圖5所示.

    由圖4可知,最大變形為0.015 mm,最大應力為1.47 MPa.最大變形發(fā)生于基座兩翼部分,此處正是側(cè)力加載機構(gòu)所處的位置,而這個系統(tǒng)最關(guān)鍵的部分就是此處的側(cè)力加載機構(gòu),其精度直接影響測量的準確度.應力范圍為0.16~1.47 MPa,遠遠小于鑄鐵的極限強度,基座可以進一步減輕質(zhì)量.

    1.3模態(tài)分析

    模態(tài)是結(jié)構(gòu)的一種固有振動特性,是一種研究結(jié)構(gòu)動力特性的方法.每個結(jié)構(gòu)都具有各自的模態(tài)參數(shù),而模態(tài)參數(shù)的識別大致分為兩種方法:一是通過計算機的有限元軟件進行計算獲得;另一種是模態(tài)試驗的方法,通過試驗設(shè)備采集到輸出信號,經(jīng)過參數(shù)識別從而獲得模態(tài)參數(shù).本文采用有限元軟件進行計算獲取模態(tài)參數(shù).

    自由度的振動系統(tǒng)可表示為[3]:

    [M]{x··}+[ξ]{x·}+[K]{x}={f}(2)

    式中:[M],[K]分別為質(zhì)量和剛度矩陣;{f}為外激勵矩陣;{x··}、{x}分別為物理坐標下的位移和加速度矩陣.

    當{f}=0時,阻尼ξ=0,式(2)為自由振動方程,其特征方程為:

    [K]-ω2[M]=0(3)

    解方程式(3)得ω的m個互異正根ω0i(i=1,2,…,m),并按升序排序,0<ω01<ω02<…<ω0m,第i階固有圓頻率

    ω0i=kimi或固有頻率f0i=12πkimi(4)

    各階固有頻率f0i均與kimi成正比,f0i越大,說明單位質(zhì)量的剛度越高,可作為結(jié)構(gòu)動態(tài)設(shè)計的一個優(yōu)化約束條件.在進行動態(tài)結(jié)構(gòu)優(yōu)化時,盡量提高各階固有頻率,避免固有頻率與外界激勵頻率一致引起共振.利用ANSYS軟件對基座進行模態(tài)分析,分析時采用靜力分析模型,對基座底面進行固定約束,不施加任何載荷,由于低階模態(tài)對振動系統(tǒng)影響較大,求得基座的前四階頻率,如表1所示.各階振型如圖5所示.

    通過對基座進行自由模態(tài)分析,得到了基座的固有頻率和固有振型,真實地反映出基座的動態(tài)特性.表1和圖6表明基座具有良好的振動性能,同時也為后續(xù)基座結(jié)構(gòu)優(yōu)化分析設(shè)計做了準備.

    2基座的拓撲優(yōu)化

    靜力學方面的拓撲優(yōu)化問題一般可以按兩種方式建立優(yōu)化模型:一是在體積或質(zhì)量約束下求最小柔度,即最大剛度;二是在剛度(或其他條件)約束下求最小體積或質(zhì)量[4].本文是在頻率、剛度及應力約束下采用漸進結(jié)構(gòu)優(yōu)化法(evolutionary structural optimization,簡稱ESO法[5-6])求解最輕的質(zhì)量.

    所謂ESO法就是預先選取結(jié)構(gòu)設(shè)計域和初始域,給定結(jié)構(gòu)所需滿足的目標函數(shù)和約束條件的貢獻度.根據(jù)其貢獻度大小刪除無效的材料,循環(huán)該過程,達到終止條件,最終得到最佳拓撲結(jié)構(gòu).在優(yōu)化迭代過程中,采用固定的有限元網(wǎng)格,對存在的材料單元,其材料數(shù)編號為非零的數(shù),而對不存在的材料單元,其材料數(shù)編號為零,通過這種零和非零模式實現(xiàn)結(jié)構(gòu)拓撲優(yōu)化.該方法具有良好的通用性,能夠采用已有的有限元分析軟件,通過迭代,在計算機上實現(xiàn).

    2.1靜力優(yōu)化準則

    2.1.1應力準則

    在某一結(jié)構(gòu)中,如果某一部分的材料沒有得到充分的利用,應將其刪除或改進,如此反復進行,直至得到最優(yōu)應力分布[7].基于此原理,通過對基座進行應力分析,確定單元的Von Mises應力σvme與整個結(jié)構(gòu)的最大Von Mises應力σvmmax,然后刪除不滿足條件的單元,減小基座質(zhì)量.對于各項同性的材料,其Von Mises應力定義為:

    σvm=σ2x+σ2y-σxσy+3τ2xy(5)

    式中:σx和σy分別為x和y方向的正應力;τxy為剪應力.

    將單元的Von Mises應力σvme與整個結(jié)構(gòu)的最大Von Mises應力σvmmax進行比較,便可確定每個單元的應力水平.在每一次有限元分析后,從模型中刪除滿足下列條件的所有單元:

    σvmeσvmmax

    式中,RRi為當前的拒絕率.

    采用相同的拒絕率,反復迭代有限元分析和單元刪除循環(huán),直至達到穩(wěn)態(tài),然后引進一個進化率ER,并附加到拒絕率中,即:

    RRi+1=RRi+ER, i=0,1,2,…(7)

    用這個增加的拒絕率再一次進行有限元分析和單元刪除,直至達到一個新的穩(wěn)態(tài).

    2.1.2剛度準則

    對基座而言,在改變其結(jié)構(gòu)時,必須保證其剛度,但是部分材料的剛度不會對結(jié)構(gòu)整體產(chǎn)生很大的影響,因此類比于應力準則將沒有得到充分利用的材料刪除,以此減小基座的質(zhì)量.根據(jù)式(1),結(jié)構(gòu)的平均柔順度,即剛度的對立量:

    C=12[F]T[K]-1[F](8)

    在載荷向量不變的情況下,整體剛度的最大化也就是結(jié)構(gòu)柔順度的最小化.從結(jié)構(gòu)中刪除第i號單元時,式(1)就變?yōu)椋?/p>

    ([K]-[Ki1])·([u]+Δ[u])=[F](9)

    式中:[Ki1]為對第i號單元擴維了的剛度矩陣;Δ[u]為位移向量的變化.

    用式(1)減去式(9),并忽略式(9)中的高階量,得到位移的變化量:

    Δ[u]=[K]-1[Ki1][u](10)

    將式(10)帶入式(8)得:

    ΔC=12[F]TΔ[u]=12[F]T[K]-1[Ki1][u]=

    12[u]T[Ki1][u]=12[ui]T[Ki][ui](11)

    式中:[Ki]為第i號單元的剛度;[ui]為第i號單元的位移向量.

    由此得到第i號單元的靈敏度:

    αi=[ui]T[Ki][ui](12)

    式中:αi表示刪除第i號單元時的結(jié)構(gòu)柔順度的變化情況.

    優(yōu)化目標就是尋找剛度約束的最輕結(jié)構(gòu),給出如下的約束:

    C≤C*(13)

    式中:C*為C的上限.

    為了達到優(yōu)化結(jié)果,在優(yōu)化過程中刪除αi值最小的單元,以使得C的增加量最小.

    2.2建立數(shù)學模型

    根據(jù)上述分析,基座具有良好的動態(tài)性能,其一階頻率遠高于基座上電機的頻率(25 Hz),因此限定某一固定頻率,以頻率、剛度、應力作為約束,最小質(zhì)量作為目標函數(shù)建立如下數(shù)學模型:

    findMinM

    s.t.f=fd

    C≤C*

    σvmeσvmmax≤1%(14)

    式中:M為基座的質(zhì)量;fd為基座的某一固有頻率.

    在優(yōu)化拓撲過程中容易出現(xiàn)數(shù)值計算不穩(wěn)定的情況,產(chǎn)生類似棋盤格的模式.可以通過引入過濾函數(shù)[8-11]來抑制棋盤格,對單元的靈敏度進行松弛調(diào)整,避免相鄰單元靈敏度相差懸殊.將編號為1到j(luò)-1的單元靈敏度通過權(quán)系數(shù)μij加到第j個單元上,即:

    u′j=∑i

    定義加權(quán)系數(shù):

    μij=d(i,j)∑Nj=1d(i,j)

    定義單元靈敏度:

    u′i=∑ji=1d(i,j)ui∑ji=1d(i,j)(15)

    式中:d(i,j)=R-D(Ni-Nj);D(Ni-Nj)為單元i到單元j的距離;R為過濾半徑,其值為2倍的網(wǎng)格單元邊長.

    2.3優(yōu)化結(jié)果分析

    依據(jù)上述數(shù)學模型,設(shè)定基座質(zhì)量下降20%,對基座進行拓撲優(yōu)化,優(yōu)化后的拓撲圖如圖7所示.圖7中深黑色材料部分表示應該切除,但是出于實際制造、加工和使用等方面考慮,不是所有的部分都應該完全切除,也可以采取減小厚度的辦法,降低基座的質(zhì)量.另外,圖7中認為要保留的部分,出于實際考慮有些部分也應該刪除.例如,在鑄造時,需要設(shè)置鑄造出砂孔,就應在內(nèi)部筋板上合理布置一些岀砂孔,優(yōu)化后的基座結(jié)構(gòu)如圖8所示.

    對優(yōu)化后的基座結(jié)構(gòu)進行有限元分析,優(yōu)化前后性能參數(shù)如表2所示.基座最大應力由原來的1.47 MPa變?yōu)?.16 MPa,但仍滿足強度要求,最大變形基本無變化,說明整體剛性良好,基頻增加了47.8 Hz,質(zhì)量由原來的5 970 kg減小到5 620 kg,減小了5.8%.

    3結(jié)論

    在正確建立模型結(jié)構(gòu)的基礎(chǔ)上,對模型做了靜力、模態(tài)分析,得到了初始模型的靜、動態(tài)性能參數(shù),分析了各參數(shù)對基座性能的影響.運用ESO法對基座進行拓撲優(yōu)化,減小了基座質(zhì)量,且整體剛性良好,基頻得到了提升,達到了預期的目標.

    參考文獻:

    [1]陳葉林,丁曉紅,郭春星,等.機床床身結(jié)構(gòu)優(yōu)化設(shè)計方法[J].機械設(shè)計,2010,27(8):65-68.

    [2]倪曉宇,易紅,湯文成,等.機床床身結(jié)構(gòu)的有限元分析與優(yōu)化[J].制造技術(shù)與機床,2005(2):47-50.

    [3]張學玲,徐燕申,鐘偉泓.基于有限元分析的數(shù)控機床床身結(jié)構(gòu)動態(tài)優(yōu)化設(shè)計方法研究[J].機械強度,2005,27(3):353-357.

    [4]焦洪宇,周奇才,李文軍,等.基于變密度法的周期性拓撲優(yōu)化[J].機械工程學報,2013,49(13):132-138.

    [5]林丹益,李芳.基于ANSYS的結(jié)構(gòu)拓撲優(yōu)化[J].機電工程,2012,29(8):898-901.

    [6]榮見華,葛森,鄧果,等.基于位移和應力靈敏度的結(jié)構(gòu)拓撲優(yōu)化設(shè)計[J].力學學報,2009,41(4):518-529.

    [7]榮見華.漸進結(jié)構(gòu)優(yōu)化方法及其應用研究[D].長沙:國防科學技術(shù)大學,2006.

    [8]HUANG X D,XIE Y M.Optimal design of periodic structures using evolutionary topology optimization[J].Structural and Multidisciplinary Optimization,2008,36(6):597-606.

    [9]葉紅玲,沈靜嫻,隋允康.過濾函數(shù)對應力約束連續(xù)體結(jié)構(gòu)拓撲優(yōu)化的影響分析[J].北京工業(yè)大學學報,2013,39(3):321-330.

    [10]葉紅玲,李耀明,張顏明,等.基于對數(shù)型Heaviside近似函數(shù)作為過濾函數(shù)的動力響應結(jié)構(gòu)拓撲優(yōu)化ICM方法[J].工程力學,2014,31(6):13-20.

    [11]尚珍,隋允康.結(jié)構(gòu)拓撲優(yōu)化中不同過濾函數(shù)間關(guān)系的探討[J].力學與實踐,2011,33(2):52-57.有色金屬材料與工程第38卷第3期NONFERROUS METAL MATERIALS AND ENGINEERINGVol.38No.32017

    文章編號:2096-2983(2017)03-0172-05DOI:10.13258/j.cnki.nmme.2017.03.009

    猜你喜歡
    拓撲優(yōu)化模態(tài)分析基座
    基于NXnastran的異步電動機基座有限元強度分析
    防爆電機(2021年2期)2021-06-09 08:14:48
    心臟固定器基座注射模設(shè)計
    模具制造(2019年7期)2019-09-25 07:30:00
    超大型FPSO火炬塔及船體基座設(shè)計
    基于負載均衡的電力通信流量分析和預測方法研究
    基于ANSYS workbench六片斜葉圓盤渦輪攪拌器的模態(tài)分析
    基于Ansys的礦用局部通風機葉輪模態(tài)分析
    電推衛(wèi)星方圓過渡結(jié)構(gòu)輕量化設(shè)計
    某調(diào)速型液力偶合器泵輪的模態(tài)分析
    基于OptiStruct的襟翼搖臂的拓撲優(yōu)化
    科技視界(2016年22期)2016-10-18 15:37:54
    民用飛機沖壓渦輪機的動剛度分析
    科技視界(2015年25期)2015-09-01 16:34:55
    久久久午夜欧美精品| 国产精品三级大全| 精品酒店卫生间| 国产高清三级在线| 麻豆成人午夜福利视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 99久久九九国产精品国产免费| 欧美另类亚洲清纯唯美| 国产成人福利小说| 日韩av不卡免费在线播放| 亚洲性久久影院| 欧美成人精品欧美一级黄| 国产成人福利小说| 亚洲欧美日韩无卡精品| 国产老妇女一区| 国产91av在线免费观看| 亚洲自偷自拍三级| 亚洲色图av天堂| 国产乱来视频区| 日韩欧美三级三区| 1000部很黄的大片| 91精品一卡2卡3卡4卡| 免费看av在线观看网站| 美女内射精品一级片tv| 草草在线视频免费看| 亚洲国产日韩欧美精品在线观看| 国产麻豆成人av免费视频| 免费观看人在逋| 成人国产麻豆网| 岛国毛片在线播放| 国产单亲对白刺激| 成人午夜高清在线视频| 日本黄色视频三级网站网址| 国产麻豆成人av免费视频| 国产一区二区在线观看日韩| 美女高潮的动态| 嫩草影院精品99| 午夜福利高清视频| 中文字幕av在线有码专区| 中文字幕亚洲精品专区| 丝袜美腿在线中文| eeuss影院久久| 国产成人一区二区在线| 国产精品嫩草影院av在线观看| 我的女老师完整版在线观看| 久久精品国产亚洲av涩爱| 久久精品久久久久久久性| 亚洲国产欧美在线一区| 1000部很黄的大片| 亚洲丝袜综合中文字幕| 又爽又黄无遮挡网站| 亚洲美女视频黄频| 国产精品国产高清国产av| 国产av一区在线观看免费| 狠狠狠狠99中文字幕| 精品人妻视频免费看| 插逼视频在线观看| 一个人免费在线观看电影| 免费黄色在线免费观看| 国产 一区精品| 22中文网久久字幕| 丝袜喷水一区| 久久精品国产99精品国产亚洲性色| 亚洲av一区综合| 免费不卡的大黄色大毛片视频在线观看 | 中文资源天堂在线| 男人狂女人下面高潮的视频| videossex国产| 深爱激情五月婷婷| 亚洲欧美精品自产自拍| 久久久久免费精品人妻一区二区| 欧美高清性xxxxhd video| 热99re8久久精品国产| 欧美一区二区亚洲| 国产精品av视频在线免费观看| 亚洲国产高清在线一区二区三| 亚洲在线观看片| 美女高潮的动态| 久久久国产成人免费| 免费观看a级毛片全部| 亚洲av电影在线观看一区二区三区 | 欧美成人免费av一区二区三区| 嫩草影院入口| 99国产精品一区二区蜜桃av| 久久精品久久精品一区二区三区| 最近视频中文字幕2019在线8| 亚洲精品国产成人久久av| 国产高清国产精品国产三级 | 国产精品女同一区二区软件| 最近中文字幕高清免费大全6| 亚洲aⅴ乱码一区二区在线播放| 国模一区二区三区四区视频| 春色校园在线视频观看| 一本久久精品| 国产亚洲精品久久久com| 男插女下体视频免费在线播放| 免费观看的影片在线观看| 久久精品国产亚洲网站| 亚洲人与动物交配视频| 能在线免费观看的黄片| 午夜激情欧美在线| 老司机福利观看| 高清在线视频一区二区三区 | 国产白丝娇喘喷水9色精品| 成人高潮视频无遮挡免费网站| 国语自产精品视频在线第100页| av天堂中文字幕网| 中文亚洲av片在线观看爽| 亚洲中文字幕日韩| 亚洲欧美日韩东京热| 免费无遮挡裸体视频| 免费无遮挡裸体视频| 亚洲国产高清在线一区二区三| 天堂中文最新版在线下载 | 久久久午夜欧美精品| 日韩人妻高清精品专区| 1000部很黄的大片| 91久久精品电影网| 天天躁夜夜躁狠狠久久av| 2022亚洲国产成人精品| 欧美又色又爽又黄视频| 日本与韩国留学比较| 亚洲欧美精品专区久久| 只有这里有精品99| 亚洲精品日韩av片在线观看| 深夜a级毛片| 日日撸夜夜添| av又黄又爽大尺度在线免费看 | 国产淫片久久久久久久久| av女优亚洲男人天堂| 一二三四中文在线观看免费高清| 日本免费a在线| 99国产精品一区二区蜜桃av| 97热精品久久久久久| 男人和女人高潮做爰伦理| 三级毛片av免费| 在线播放国产精品三级| 大又大粗又爽又黄少妇毛片口| 在线播放无遮挡| 日韩在线高清观看一区二区三区| 国产亚洲精品av在线| 亚洲国产精品合色在线| 在线天堂最新版资源| 午夜福利网站1000一区二区三区| 国产午夜福利久久久久久| 两个人的视频大全免费| 天天躁夜夜躁狠狠久久av| 国产精品福利在线免费观看| 国产成人精品婷婷| 免费播放大片免费观看视频在线观看 | 波多野结衣巨乳人妻| 久久这里有精品视频免费| 欧美变态另类bdsm刘玥| 欧美日本视频| 久久鲁丝午夜福利片| 亚洲色图av天堂| 久久热精品热| 男人舔女人下体高潮全视频| 国产男人的电影天堂91| 九九热线精品视视频播放| 国产免费视频播放在线视频 | 男女边吃奶边做爰视频| 91久久精品国产一区二区三区| av视频在线观看入口| 少妇猛男粗大的猛烈进出视频 | 欧美成人精品欧美一级黄| 精品99又大又爽又粗少妇毛片| 级片在线观看| 最近最新中文字幕大全电影3| a级一级毛片免费在线观看| 国产精品乱码一区二三区的特点| 欧美最新免费一区二区三区| 亚洲av二区三区四区| 九九在线视频观看精品| 日韩人妻高清精品专区| 日韩精品有码人妻一区| 亚洲av.av天堂| 床上黄色一级片| 欧美一区二区亚洲| 18禁裸乳无遮挡免费网站照片| 久久久色成人| 超碰97精品在线观看| 中文字幕亚洲精品专区| 十八禁国产超污无遮挡网站| 亚洲,欧美,日韩| 夜夜看夜夜爽夜夜摸| 日本五十路高清| 深夜a级毛片| 亚洲国产欧美在线一区| av视频在线观看入口| 91精品伊人久久大香线蕉| 男插女下体视频免费在线播放| 建设人人有责人人尽责人人享有的 | h日本视频在线播放| 寂寞人妻少妇视频99o| av在线观看视频网站免费| 色吧在线观看| 日韩强制内射视频| 日韩欧美精品免费久久| 午夜爱爱视频在线播放| 欧美成人午夜免费资源| 亚洲成av人片在线播放无| h日本视频在线播放| 中文字幕av成人在线电影| 亚洲国产精品sss在线观看| 国产成人91sexporn| 成人漫画全彩无遮挡| 日韩高清综合在线| 99热这里只有是精品在线观看| 国产在线一区二区三区精 | 国产成人一区二区在线| 夜夜看夜夜爽夜夜摸| 亚洲av中文av极速乱| 亚洲欧洲国产日韩| a级毛片免费高清观看在线播放| 99久久无色码亚洲精品果冻| 久久久久网色| 成年av动漫网址| 国产精品人妻久久久久久| 欧美三级亚洲精品| 免费播放大片免费观看视频在线观看 | 3wmmmm亚洲av在线观看| 床上黄色一级片| 性插视频无遮挡在线免费观看| 久久精品国产99精品国产亚洲性色| 日韩国内少妇激情av| 波多野结衣高清无吗| 亚洲18禁久久av| 性插视频无遮挡在线免费观看| 美女国产视频在线观看| 在线观看美女被高潮喷水网站| 能在线免费观看的黄片| 亚洲人成网站高清观看| 亚洲av二区三区四区| 汤姆久久久久久久影院中文字幕 | 欧美精品一区二区大全| av免费在线看不卡| 成人午夜高清在线视频| 99热6这里只有精品| 精品国内亚洲2022精品成人| 美女脱内裤让男人舔精品视频| 五月伊人婷婷丁香| 激情 狠狠 欧美| 久久精品夜色国产| 国产高清三级在线| 搞女人的毛片| 三级国产精品欧美在线观看| 久久这里只有精品中国| 免费观看人在逋| 午夜福利成人在线免费观看| 精品久久久久久久人妻蜜臀av| 国产高清不卡午夜福利| 在线观看av片永久免费下载| 精品午夜福利在线看| 国产亚洲精品久久久com| 国产乱人偷精品视频| 国产伦一二天堂av在线观看| 亚洲精品亚洲一区二区| 五月玫瑰六月丁香| 精品久久久久久久末码| 国产久久久一区二区三区| 国产精品乱码一区二三区的特点| 国产午夜精品久久久久久一区二区三区| 精品午夜福利在线看| 午夜精品在线福利| 国产精品久久久久久精品电影小说 | 三级国产精品欧美在线观看| 日本-黄色视频高清免费观看| 人妻系列 视频| 欧美日韩国产亚洲二区| 精品久久久久久久人妻蜜臀av| 九九久久精品国产亚洲av麻豆| 久久久成人免费电影| 丝袜美腿在线中文| 91精品国产九色| 久久精品国产亚洲av涩爱| 亚洲激情五月婷婷啪啪| 变态另类丝袜制服| 国产精品嫩草影院av在线观看| 中文在线观看免费www的网站| 色综合色国产| 国产一区二区三区av在线| 最近中文字幕高清免费大全6| 欧美丝袜亚洲另类| a级一级毛片免费在线观看| 国产成人精品久久久久久| 亚洲av.av天堂| 国语对白做爰xxxⅹ性视频网站| 午夜精品国产一区二区电影 | 哪个播放器可以免费观看大片| 免费大片18禁| 久久99热6这里只有精品| 国产精品不卡视频一区二区| 身体一侧抽搐| 亚洲18禁久久av| 2021少妇久久久久久久久久久| 中文字幕人妻熟人妻熟丝袜美| 99视频精品全部免费 在线| 美女国产视频在线观看| 18禁在线无遮挡免费观看视频| 午夜福利成人在线免费观看| www.av在线官网国产| 色噜噜av男人的天堂激情| 日本av手机在线免费观看| 麻豆乱淫一区二区| 久久久欧美国产精品| 国产女主播在线喷水免费视频网站 | 波多野结衣巨乳人妻| 永久免费av网站大全| 国产精品一区二区三区四区久久| 一区二区三区乱码不卡18| 久久草成人影院| 一个人看视频在线观看www免费| 欧美日韩在线观看h| 偷拍熟女少妇极品色| 亚洲一区高清亚洲精品| 欧美人与善性xxx| 欧美日本亚洲视频在线播放| 欧美一区二区亚洲| 少妇熟女欧美另类| 国产午夜精品论理片| 男人舔奶头视频| 国国产精品蜜臀av免费| 国内精品宾馆在线| 黄色欧美视频在线观看| 大话2 男鬼变身卡| 天堂√8在线中文| 国产黄色视频一区二区在线观看 | 亚洲国产欧美在线一区| 久久久久网色| a级一级毛片免费在线观看| 欧美又色又爽又黄视频| 99九九线精品视频在线观看视频| 人妻系列 视频| 国产伦一二天堂av在线观看| 97人妻精品一区二区三区麻豆| 日日干狠狠操夜夜爽| 国产一级毛片在线| 人妻制服诱惑在线中文字幕| 五月伊人婷婷丁香| 国产免费视频播放在线视频 | 性色avwww在线观看| 男人狂女人下面高潮的视频| 久久精品国产亚洲av涩爱| 久久久久精品久久久久真实原创| 欧美成人a在线观看| 亚洲人成网站高清观看| 午夜福利高清视频| av在线天堂中文字幕| 又爽又黄无遮挡网站| 黄色配什么色好看| 热99在线观看视频| 日本一二三区视频观看| 青春草视频在线免费观看| 国产av在哪里看| 久久精品久久精品一区二区三区| 高清毛片免费看| 国产免费福利视频在线观看| 日韩强制内射视频| 欧美另类亚洲清纯唯美| 精品不卡国产一区二区三区| 日韩欧美三级三区| 日韩 亚洲 欧美在线| 中文资源天堂在线| 国产 一区精品| 国产亚洲精品久久久com| 免费人成在线观看视频色| 久久精品人妻少妇| 九九热线精品视视频播放| 国产精品伦人一区二区| 成年版毛片免费区| 国产精品人妻久久久影院| 久99久视频精品免费| 国产又色又爽无遮挡免| 日本午夜av视频| 能在线免费看毛片的网站| 啦啦啦观看免费观看视频高清| 久久欧美精品欧美久久欧美| 日韩成人av中文字幕在线观看| 久久久久久久久久久丰满| 亚洲av日韩在线播放| 欧美性猛交黑人性爽| 又爽又黄无遮挡网站| 搞女人的毛片| 国产伦精品一区二区三区视频9| 日本一二三区视频观看| 赤兔流量卡办理| 男女下面进入的视频免费午夜| 日韩高清综合在线| 美女被艹到高潮喷水动态| 成人三级黄色视频| 亚洲综合精品二区| 天堂av国产一区二区熟女人妻| 午夜福利在线在线| 欧美成人a在线观看| 天堂网av新在线| 黄片无遮挡物在线观看| 毛片一级片免费看久久久久| 国产色婷婷99| 欧美日韩在线观看h| 能在线免费看毛片的网站| 久久久久久久久中文| 免费播放大片免费观看视频在线观看 | 久久精品综合一区二区三区| 男人狂女人下面高潮的视频| 99视频精品全部免费 在线| 久久久精品94久久精品| 午夜福利成人在线免费观看| 全区人妻精品视频| 激情 狠狠 欧美| 亚洲av电影不卡..在线观看| 夜夜看夜夜爽夜夜摸| 久久久精品欧美日韩精品| 男人狂女人下面高潮的视频| 九九在线视频观看精品| 最近视频中文字幕2019在线8| 晚上一个人看的免费电影| 国产免费福利视频在线观看| 欧美一区二区国产精品久久精品| 婷婷色av中文字幕| 你懂的网址亚洲精品在线观看 | 99热这里只有是精品50| 午夜福利在线在线| 性色avwww在线观看| 中文资源天堂在线| 亚洲av.av天堂| 熟女人妻精品中文字幕| 亚洲欧美日韩东京热| 亚洲18禁久久av| 久久久国产成人精品二区| 欧美三级亚洲精品| 深爱激情五月婷婷| 国产一区有黄有色的免费视频 | 97超碰精品成人国产| 亚洲怡红院男人天堂| 国产亚洲一区二区精品| 日本免费在线观看一区| 一级毛片久久久久久久久女| 一级爰片在线观看| 久久99精品国语久久久| 男女视频在线观看网站免费| 日韩一区二区视频免费看| eeuss影院久久| 熟女人妻精品中文字幕| 黄色配什么色好看| 少妇的逼好多水| av线在线观看网站| 蜜臀久久99精品久久宅男| 国产亚洲av嫩草精品影院| 乱人视频在线观看| 国产黄色视频一区二区在线观看 | 成人午夜高清在线视频| 美女被艹到高潮喷水动态| 日本wwww免费看| 午夜激情福利司机影院| 欧美日韩综合久久久久久| 美女xxoo啪啪120秒动态图| 男女那种视频在线观看| av国产免费在线观看| 在现免费观看毛片| 少妇高潮的动态图| 日本猛色少妇xxxxx猛交久久| 蜜桃亚洲精品一区二区三区| 精品人妻偷拍中文字幕| 成人无遮挡网站| 国产探花极品一区二区| 少妇的逼好多水| 内地一区二区视频在线| 日韩av不卡免费在线播放| 禁无遮挡网站| 国产高清三级在线| 我的老师免费观看完整版| 岛国在线免费视频观看| 久久精品熟女亚洲av麻豆精品 | 精品不卡国产一区二区三区| 国产黄色小视频在线观看| 国产激情偷乱视频一区二区| 久久久久久伊人网av| 建设人人有责人人尽责人人享有的 | 欧美色视频一区免费| av福利片在线观看| 少妇裸体淫交视频免费看高清| 精品久久久久久久久久久久久| 人人妻人人澡欧美一区二区| 久久久国产成人免费| 精品国产露脸久久av麻豆 | 色综合亚洲欧美另类图片| 色视频www国产| 日韩欧美 国产精品| 日韩欧美精品免费久久| 内射极品少妇av片p| 97在线视频观看| 中文字幕久久专区| 国产精品久久久久久精品电影| 欧美成人精品欧美一级黄| 亚洲自拍偷在线| 国产精品熟女久久久久浪| 亚洲av二区三区四区| 偷拍熟女少妇极品色| 小蜜桃在线观看免费完整版高清| 亚洲18禁久久av| 乱码一卡2卡4卡精品| 色噜噜av男人的天堂激情| 色视频www国产| 国产精品国产三级专区第一集| 在现免费观看毛片| 国产av码专区亚洲av| 嫩草影院入口| 看免费成人av毛片| 成人特级av手机在线观看| 老女人水多毛片| 日本wwww免费看| 老司机福利观看| 日本午夜av视频| 婷婷色av中文字幕| 国产精品无大码| 韩国av在线不卡| 边亲边吃奶的免费视频| 男插女下体视频免费在线播放| 亚洲欧美日韩东京热| 久久久久久久午夜电影| 精品国内亚洲2022精品成人| 热99在线观看视频| 免费看日本二区| 久久99精品国语久久久| 人妻系列 视频| 亚洲精品日韩在线中文字幕| 在现免费观看毛片| 日本免费在线观看一区| av国产久精品久网站免费入址| 国产精品一区二区三区四区免费观看| 综合色av麻豆| 免费观看精品视频网站| 最近最新中文字幕免费大全7| 中文字幕人妻熟人妻熟丝袜美| 日本色播在线视频| 国产v大片淫在线免费观看| 亚洲欧美精品专区久久| 国产精品一二三区在线看| 国产精品国产三级专区第一集| 精品无人区乱码1区二区| 国产乱人视频| 亚洲欧美日韩高清专用| 你懂的网址亚洲精品在线观看 | 精品酒店卫生间| 国产免费视频播放在线视频 | 国产精品1区2区在线观看.| 女人十人毛片免费观看3o分钟| 熟妇人妻久久中文字幕3abv| 午夜亚洲福利在线播放| 中国美白少妇内射xxxbb| 欧美性感艳星| 日韩在线高清观看一区二区三区| a级毛片免费高清观看在线播放| 国产欧美日韩精品一区二区| 高清日韩中文字幕在线| 精品久久久久久久久久久久久| 草草在线视频免费看| 久久精品国产亚洲av天美| 久久99热这里只有精品18| 直男gayav资源| 久久久久精品久久久久真实原创| 免费黄网站久久成人精品| 热99在线观看视频| 欧美成人午夜免费资源| 精品一区二区免费观看| 中文在线观看免费www的网站| 亚洲av男天堂| 国产乱人偷精品视频| 色吧在线观看| 日韩视频在线欧美| АⅤ资源中文在线天堂| 免费一级毛片在线播放高清视频| 亚洲精品国产成人久久av| 91狼人影院| 精品免费久久久久久久清纯| 欧美日本视频| or卡值多少钱| 亚洲无线观看免费| 日韩成人av中文字幕在线观看| 亚洲最大成人中文| 91久久精品国产一区二区成人| 两性午夜刺激爽爽歪歪视频在线观看| 69人妻影院| .国产精品久久| 午夜福利视频1000在线观看| 大话2 男鬼变身卡| 成人性生交大片免费视频hd| 欧美精品一区二区大全| 免费黄色在线免费观看| 人人妻人人澡欧美一区二区| 国产v大片淫在线免费观看| 在线观看一区二区三区| 亚洲一区高清亚洲精品| 日韩中字成人| 日本三级黄在线观看| 丰满人妻一区二区三区视频av| 国产成人精品一,二区| 国产黄a三级三级三级人| 国产午夜精品一二区理论片| 亚洲av电影在线观看一区二区三区 | 97人妻精品一区二区三区麻豆| 秋霞在线观看毛片| 国产精品国产三级专区第一集| 国产亚洲一区二区精品| 五月伊人婷婷丁香| 久久久久久久国产电影| 99久久精品热视频| 欧美97在线视频| 欧美性感艳星| 精品酒店卫生间| 日韩欧美精品免费久久| 91久久精品电影网| 久久精品久久久久久久性| 哪个播放器可以免费观看大片|