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

    基于精確幾何模型梁?jiǎn)卧穆菪龔椈蓜偠确治?/h1>
    2020-02-10 09:56:12齊朝暉卓英鵬國(guó)樹東
    工程力學(xué) 2020年2期
    關(guān)鍵詞:曲梁螺旋彈簧曲率

    張 健,齊朝暉,卓英鵬,國(guó)樹東

    (大連理工大學(xué)工業(yè)裝備結(jié)構(gòu)分析國(guó)家重點(diǎn)實(shí)驗(yàn)室,大連 116024)

    螺旋彈簧在機(jī)械系統(tǒng)中有著重要影響,尤其在車用發(fā)動(dòng)機(jī)配氣機(jī)構(gòu)和懸架系統(tǒng)等復(fù)雜機(jī)構(gòu)上起著重要作用[1-2]。螺旋彈簧的剛度特性是設(shè)計(jì)彈簧的重要評(píng)價(jià)因素。

    在復(fù)雜機(jī)械系統(tǒng)建模中,將螺旋彈簧視為單自由度振動(dòng)模型或是采用將彈簧按照圈數(shù)分割為集中質(zhì)量,質(zhì)量與質(zhì)量之間通過連接具有當(dāng)量剛度的無質(zhì)量彈簧[3],該方法忽略了彈簧的非線性特性,不能準(zhǔn)確地反映彈簧的非線性剛度。采用有限元分析能獲得比較精確的彈簧剛度特性,因而,有限元分析廣泛應(yīng)用于對(duì)彈簧特性的評(píng)價(jià)。早期建立彈簧有限元模型是利用傳統(tǒng)梁?jiǎn)卧?,將彈簧分割成很多部分,每部分之間通過直梁?jiǎn)卧B接,該方法忽略了彈簧本身是空間彎曲結(jié)構(gòu),故而彈簧模型離散的疏密程度和梁?jiǎn)卧趶椈芍械谋壤紩?huì)影響計(jì)算的精度和計(jì)算速度[4]。文獻(xiàn)[5]應(yīng)用直梁?jiǎn)卧慕K枷虢⑶簡(jiǎn)卧C明對(duì)于螺旋彈簧結(jié)構(gòu)中使用曲梁?jiǎn)卧木纫黠@高于直梁?jiǎn)卧?,并且同等精度的情況下使用更少的單元數(shù)量。文獻(xiàn)[6]提出了節(jié)點(diǎn)自由度為 6的 2節(jié)點(diǎn)三維曲梁有限單元模型,在一定程度上是一種高效彈簧建模方法。同時(shí),文獻(xiàn)[7]作為彈簧設(shè)計(jì)指導(dǎo)手冊(cè),對(duì)彈簧設(shè)計(jì)做出了比較系統(tǒng)的分析,提出了經(jīng)典理論剛度公式,在其有限元分析案例中采用三維實(shí)體單元,這種分析手段方便成熟,適合于對(duì)彈簧特性的評(píng)價(jià),但并未考慮彈簧彎曲變形對(duì)彈簧剛度的影響。文獻(xiàn)[8]中的彈簧模型考慮了彎曲變形,但所建立的曲梁?jiǎn)卧贿m合于小變形情況。

    螺旋彈簧模型應(yīng)選擇一種幾何非線性分析方法。對(duì)于細(xì)長(zhǎng)結(jié)構(gòu),近幾年發(fā)展起來的幾何非線性分析方法有絕對(duì)節(jié)點(diǎn)坐標(biāo)法和“精確幾何模型梁?jiǎn)卧薄?/p>

    絕對(duì)節(jié)點(diǎn)坐標(biāo)法是放棄了傳統(tǒng)梁理論中的剛性截面假設(shè),將梁看成一種特殊實(shí)體,可精確描述單元?jiǎng)傮w位移的變形場(chǎng),它是由 Shabana[9]首次提出,由于該方法避免了對(duì)截面有限轉(zhuǎn)動(dòng)插值引起的各種困難,是柔體動(dòng)力學(xué)研究的一個(gè)重要進(jìn)展,受到了國(guó)內(nèi)外學(xué)者的廣泛關(guān)注[10-11]。

    “精確幾何模型梁?jiǎn)卧眲t是一種摒棄了小位移小轉(zhuǎn)動(dòng)假設(shè)的幾何非線性分析方法。選取描述單元端部形心位置的位移變量和描述單元端部橫截面有限轉(zhuǎn)動(dòng)的角度變量作為獨(dú)立的描述變量,在單元域內(nèi)進(jìn)行插值離散。而傳統(tǒng)“精確幾何模型梁?jiǎn)卧敝挥袆傂越孛婕僭O(shè),分析細(xì)長(zhǎng)結(jié)構(gòu)的時(shí)候會(huì)出現(xiàn)“剪切閉鎖”,結(jié)合 Bernoulli梁理論和細(xì)長(zhǎng)梁的幾何關(guān)系,可避免對(duì)截面轉(zhuǎn)動(dòng)矢量進(jìn)行獨(dú)立插值,從而避免傳統(tǒng)方法當(dāng)中的很多問題[12-15]。

    本文基于“精確幾何模型梁?jiǎn)卧碧岢隽艘环N新的螺旋彈簧剛度分析方法。結(jié)合螺旋彈簧的設(shè)計(jì)構(gòu)型特點(diǎn)和變形規(guī)律,尋找螺旋彈簧的描述參數(shù);基于傳統(tǒng)“精確幾何模型梁?jiǎn)卧睅缀畏蔷€性分析的方法,同時(shí)增加形心線與剛截面垂直假設(shè)條件,建立彈簧系統(tǒng)虛功率方程,并通過數(shù)值算例進(jìn)行剛度特性分析,對(duì)比與經(jīng)典理論剛度計(jì)算結(jié)果、傳統(tǒng)有限元計(jì)算結(jié)果的差異。

    1 螺旋彈簧系統(tǒng)描述

    對(duì)于任意形狀的圓柱螺旋彈簧,其基本參數(shù)包含螺旋線圓柱半徑ρ、高度z和極角θ,以描述彈簧截面形心位置。首先,以右手螺旋彈簧為例,在彈簧底圈中心處,建立總體基為如圖1所示。

    圖1 總體基和截面坐標(biāo)系Fig.1 Global coordinate system and section coordinate system

    螺旋彈簧螺旋線的形心矢徑為:

    考慮圓柱螺旋彈簧只受到軸向載荷作用時(shí),彈簧產(chǎn)生軸向變形,在變形過程中形心線的位置仍保持螺旋形狀[7],只是基本參數(shù)螺旋線圓柱半徑ρ、高度z和極角θ會(huì)發(fā)生變化,故高度z、半徑ρ和極角θ是彈簧變形的狀態(tài)變量。

    對(duì)于已有彈簧,往往已知彈簧自由狀態(tài)下高度和半徑隨極角變化的構(gòu)型曲線。只需要通過選擇適當(dāng)?shù)牟逯捣椒▽?gòu)型曲線離散,便可得到螺旋線上任意一點(diǎn)的形心矢量。

    以彈簧的每一圈為一小單元,采用Hermite插值方法,螺旋線圓柱半徑插值函數(shù)為:

    同時(shí),保證2個(gè)樣條單元之間節(jié)點(diǎn)處的二階導(dǎo)數(shù)連續(xù)性。依據(jù)這個(gè)連續(xù)性條件,單元節(jié)點(diǎn)處的半徑對(duì)初始極角的一階導(dǎo)數(shù)表示為:

    式中:ρi=[ρ0ρ1...ρn];S0和S1均為常數(shù)矩陣。

    從而可將式(2)改寫為:

    式中:

    同理,彈簧螺旋線高度為:

    彈簧螺旋線極角為:

    假設(shè)簧絲截面為剛性截面,彈簧在壓縮或者拉伸過程中,不僅彈簧螺旋線的形心矢徑發(fā)生改變,而且簧絲截面也會(huì)發(fā)生相應(yīng)的扭轉(zhuǎn)。除形心矢徑的描述變量外,還需增加一個(gè)簧絲截面的扭轉(zhuǎn)角。

    插值方法同上,有:

    故,螺旋彈簧的完備描述參數(shù)為:

    2 簧絲截面坐標(biāo)系的建立

    基于剛性截面假設(shè),簧絲截面所做的運(yùn)動(dòng)為剛體運(yùn)動(dòng),因而可以用固結(jié)于橫截面的截面坐標(biāo)系描述剛體截面的方位,如圖1所示。截面坐標(biāo)系的原點(diǎn)在彈簧形心線上,es軸是形心線的切線方向,另外et、eb兩根軸分別為彈簧截面的形心主軸。

    螺旋線形心線切矢量為:

    卡爾丹描述的截面法向矢量為:

    考慮彈簧細(xì)長(zhǎng)結(jié)構(gòu)特點(diǎn),采用 Bernoulli梁理論,變形過程中梁截面始終保持與形心線切向垂直,截面法向矢量與螺旋線的形心線切線方向一致。通過式(10)和式(11),求得卡爾丹角為:

    由式(12)和式(13)可知,α、β角完全由彈簧形心線的形狀所決定。因而只需要用第3次定軸轉(zhuǎn)動(dòng)的卡爾丹角就可以描述簧絲截面的扭轉(zhuǎn)。

    簧絲截面的其他2個(gè)矢量為截面形心主軸,分別為:

    3 彈簧單元的慣性力虛功率

    剛性簧絲截面的慣性力虛功率包含平動(dòng)慣性力虛功率和轉(zhuǎn)動(dòng)慣性力虛功率兩部分,在彈簧單元域(0,2nπ)內(nèi)積分可得到彈簧單元慣性力虛功率。

    彈簧單元平動(dòng)慣性力虛功率為:

    角速度具有明確的物理意義[16],有:

    式中,截面角速度ω在連體基中的分解為:

    考慮到方便對(duì)比推導(dǎo)曲率矢量,故將截面角速度在截面坐標(biāo)系下表達(dá)。

    采用卡爾丹描述,在總體基下的3個(gè)轉(zhuǎn)軸為:

    通過式(11)、式(14)和式(15)將其轉(zhuǎn)化到截面坐標(biāo)系下,則3個(gè)轉(zhuǎn)軸為:

    根據(jù)相繼定軸轉(zhuǎn)動(dòng),彈簧截面的角速度為:

    角速度矢量在截面坐標(biāo)系中的分量為:

    彈簧單元轉(zhuǎn)動(dòng)慣性力虛功率為:

    式中,J為剛截面的慣性矩張量。

    4 螺旋彈簧的曲率矢量與軸向應(yīng)變

    螺旋彈簧是空間曲梁結(jié)構(gòu),可以用曲率矢量表征形心線的彎曲扭轉(zhuǎn)程度,與角速度相同,曲率也具有明確的物理意義,曲率為:

    式中:s0為螺旋彈簧形心線初始弧長(zhǎng)坐標(biāo);曲率矢量κ在連體基中的分解為:

    類比角速度的疊加原理,彈簧單元截面的曲率矢量為:

    曲率矢量改寫為:

    曲率矢量在截面坐標(biāo)系中的分量為:

    曲率矢量之所以須在截面坐標(biāo)系下表達(dá),這是由曲梁的本構(gòu)關(guān)系所決定。

    剛截面作剛體運(yùn)動(dòng),可將曲梁微元體簡(jiǎn)化成一條形心線,梁的應(yīng)變只是變形前弧長(zhǎng)坐標(biāo)s0和時(shí)間的函數(shù)。變形后弧長(zhǎng)也是這兩個(gè)坐標(biāo)的函數(shù),則曲梁的軸向應(yīng)變?yōu)椋?/p>

    5 彈簧的變形虛功率

    在彈簧中取一段長(zhǎng)為ds0微元體,該微元體是具有初始曲率的曲梁微元體。

    微元體可當(dāng)成一段剛體,剛體的動(dòng)力學(xué)方程為:

    式中:m、f分別為曲梁微元體左端面合力矩和合力;

    根據(jù)虛功率原理,可得微元體的外力虛功率為:

    微元體的慣性力虛功率為:

    根據(jù)虛功率原理,曲梁微元體的變形虛功率為:

    將其沿極角積分可以得到彈簧變形虛功率:

    根據(jù)截面法向與截面形心線方向相同的假設(shè)條件,由式(10)和式(31)可得:

    式(37)兩端對(duì)時(shí)間求一階導(dǎo)數(shù),得:

    角速度ω對(duì)求一階導(dǎo)數(shù):

    將式(40)代入式(39),有:

    曲率矢量κ的一階時(shí)間導(dǎo)數(shù)為:

    將式(42)代入式(41),則有:

    將式(38)和式(43)代入式(36),可得彈簧的變形虛功率為:

    構(gòu)造曲梁?jiǎn)卧谋緲?gòu)關(guān)系可以描述為:

    整理后,彈簧的變形虛功率為:

    6 彈簧剛度的計(jì)算

    6.1 模型降噪方法

    系統(tǒng)動(dòng)力學(xué)方程具有很強(qiáng)的剛性,而高頻彈性振動(dòng)主要源于應(yīng)力的快速變化,模型降噪的方法是將應(yīng)力在一個(gè)時(shí)間區(qū)間(t,t+h)內(nèi)平均[17],用平均應(yīng)力替換應(yīng)力σ,以此來計(jì)算柔體的變形虛功率。

    平均應(yīng)力近似表達(dá)為:

    則變形虛功率為:

    將式(47)代入式(48),變形虛功率為:

    與修改前的變形虛功率相比,增加了慣性項(xiàng)和阻尼項(xiàng),從而降低了系統(tǒng)的固有頻率,可以獲得更高的計(jì)算效率。

    6.2 系統(tǒng)虛功率方程

    彈簧剛度是通過在彈簧軸向方向上施加測(cè)試外載荷獲得,在彈簧系統(tǒng)中外力虛功率包含重力虛功率和測(cè)試載荷虛功率。

    重力虛功率為:

    測(cè)試載荷虛功率為:

    式中:fn為測(cè)試載荷;

    系統(tǒng)虛功率方程為:

    7 數(shù)值算例

    算例1.圓柱壓縮螺旋彈簧基本參數(shù):中徑30.4 mm,有效圈數(shù) 4,簧絲直徑 4.1 mm,節(jié)距為11.94 mm,材料彈性模量2.1×105MPa,泊松比0.3,密度7800 kg/m3。

    彈簧模型如圖2所示,彈簧的兩端圓并緊并磨平,支撐圈數(shù)為1.5。

    圖2 彈簧模型/mmFig.2 The model of spring

    彈簧模型的曲率如圖3所示。由圖3可知,彈簧曲率不僅在小單元內(nèi)部連續(xù),而且在每圈樣條單元的節(jié)點(diǎn)拼接處也連續(xù),該種彈簧單元具有很好的光順性。

    當(dāng)彈簧支撐磨平面承受承載力的時(shí)候,忽略與支撐圈與活動(dòng)圈的壓緊,將作用載荷等效到有效圈的端部節(jié)點(diǎn)。

    圖3 彈簧初始曲率Fig.3 Spring curvature in initial state

    彈簧單元底部施加固定約束,頂部平面緩慢施加過程壓力載荷,最大載荷為300 N,載荷方程為:

    彈簧各圈高度的變化情況如下圖4所示。

    頂部施加壓力載荷 270 N~540 N,在彈簧受壓過程中,彈簧剛度變化情況,如圖5所示。

    由剛度變化曲線可知,彈簧在壓縮過程中,彈簧的剛度是漸減的,程序計(jì)算剛度值與文獻(xiàn)[7]評(píng)價(jià)一致。

    算例2.圓柱拉伸螺旋彈簧基本參數(shù):中徑18.5 mm,有效圈數(shù)10,簧絲直徑2.5 mm,節(jié)距為3 mm,材料彈性模量2.1×105MPa;泊松比0.3;密度7800 kg/m3。彈簧模型如圖6所示。

    圖4 壓縮過程中各圈高度Fig.4 Each spring height in compression process

    圖5 圓柱壓縮螺旋彈簧剛度變化曲線Fig.5 Cylindrical compression coil spring stiffness curve

    圖6 彈簧模型/mmFig.6 The model of spring

    底部固定約束,頂部施加拉力 100 N~500 N。剛度變化如圖7所示。

    由圖7剛度變化曲線可知,拉伸彈簧的剛度是漸增型,程序計(jì)算剛度值與文獻(xiàn)[7]評(píng)價(jià)一致。進(jìn)一步發(fā)現(xiàn),剛度并不是線性增加,而是非線性的。

    算例3.圓柱螺旋彈簧基本參數(shù):中徑157 mm;有效圈數(shù) 6;簧絲直徑 13 mm;材料彈性模量2.1×105MPa;泊松比0.3。頂部施加壓力1000 N。將剛度仿真結(jié)果與理論值比較。

    圖7 圓柱拉伸螺旋彈簧剛度變化曲線Fig.7 Cylindrical tensile coil spring stiffness curve

    圓柱螺旋彈簧的經(jīng)典理論剛度公式為:

    式中,G為彈簧材料的切變模量。

    計(jì)算結(jié)果如表1所示。

    表1 彈簧剛度的仿真結(jié)果與理論值比較Table 1 Comparison of spring stiffness simulation results with theoretical values

    由表1對(duì)比可知,有限元Solid95計(jì)算剛度結(jié)果與線性經(jīng)典理論計(jì)算結(jié)果接近。對(duì)比節(jié)距為40 mm的 Solid95(小變形、單元數(shù)目為 156456和 7932)與 Beam189(大變形、單元數(shù)目為 50),其計(jì)算結(jié)果基本一致,可以證明該有限元方法的正確性。同時(shí),本文方法與Solid95計(jì)算結(jié)果基本一致,可以證明本文方法的正確性。

    8 結(jié)論

    本文將彈簧半徑、高度、極角和截面扭轉(zhuǎn)角作為螺旋彈簧的描述變量,并根據(jù)設(shè)計(jì)曲線進(jìn)行插值離散?;诟倪M(jìn)后的“精確幾何梁?jiǎn)卧苯7椒?,通過形心曲線切矢量和簧絲截面扭轉(zhuǎn)角建立截面坐標(biāo)系。應(yīng)用高頻彈性振動(dòng)降噪方法,得到彈簧系統(tǒng)虛功率方程,最后求解得到彈簧剛度。通過對(duì)比傳統(tǒng)有限元模型和經(jīng)典理論剛度公式的計(jì)算結(jié)果,驗(yàn)證了該種螺旋彈簧剛度分析方法的合理性和正確性,并為變剛度彈簧設(shè)計(jì)做好了理論準(zhǔn)備。

    猜你喜歡
    曲梁螺旋彈簧曲率
    變曲率均質(zhì)梁結(jié)構(gòu)的振動(dòng)特性研究
    大曲率沉管安裝關(guān)鍵技術(shù)研究
    關(guān)于薄壁曲梁與直梁解析解的進(jìn)一步討論1)
    一類雙曲平均曲率流的對(duì)稱與整體解
    基于有限元法的波紋曲梁結(jié)構(gòu)振動(dòng)特性研究
    某懸架螺旋彈簧斷裂問題分析及改進(jìn)
    Soechting and Her Wild Animals
    半正迷向曲率的四維Shrinking Gradient Ricci Solitons
    一種收口螺旋彈簧設(shè)計(jì)方法
    汽車零部件(2015年5期)2015-01-03 08:00:48
    某車型后螺旋彈簧斷裂問題分析
    汽車零部件(2015年5期)2015-01-03 08:00:42

    人妻久久中文字幕网| 欧美激情 高清一区二区三区| 在线观看www视频免费| 精品福利观看| 大片免费播放器 马上看| 搡老乐熟女国产| 精品国产乱码久久久久久小说| 免费人妻精品一区二区三区视频| 精品免费久久久久久久清纯 | 国产又爽黄色视频| 老司机靠b影院| 侵犯人妻中文字幕一二三四区| 黑人操中国人逼视频| 淫妇啪啪啪对白视频 | 咕卡用的链子| 国产免费一区二区三区四区乱码| 国产1区2区3区精品| 久久青草综合色| 精品第一国产精品| 国产精品 欧美亚洲| 日韩大片免费观看网站| 热99re8久久精品国产| 在线观看免费午夜福利视频| 久久九九热精品免费| 视频在线观看一区二区三区| 久久人妻福利社区极品人妻图片| 十八禁网站网址无遮挡| h视频一区二区三区| 久久狼人影院| 丝袜喷水一区| 巨乳人妻的诱惑在线观看| 少妇精品久久久久久久| 久久久久久人人人人人| 久久久久久人人人人人| 老司机靠b影院| 蜜桃国产av成人99| 一边摸一边做爽爽视频免费| 777久久人妻少妇嫩草av网站| 国产在视频线精品| 搡老熟女国产l中国老女人| 精品卡一卡二卡四卡免费| 水蜜桃什么品种好| 9191精品国产免费久久| 欧美国产精品一级二级三级| 欧美精品av麻豆av| 国产高清videossex| www.精华液| √禁漫天堂资源中文www| 2018国产大陆天天弄谢| 日本一区二区免费在线视频| 日本av免费视频播放| 国产色视频综合| 午夜精品久久久久久毛片777| 99精品久久久久人妻精品| 日本vs欧美在线观看视频| 日韩大片免费观看网站| 一级片'在线观看视频| 亚洲欧洲日产国产| 久久国产精品人妻蜜桃| av一本久久久久| www.熟女人妻精品国产| 人人澡人人妻人| 在线观看免费高清a一片| 美女福利国产在线| 午夜91福利影院| 午夜91福利影院| 欧美日韩精品网址| 超碰97精品在线观看| 18禁观看日本| 欧美av亚洲av综合av国产av| 久热这里只有精品99| 国产一卡二卡三卡精品| 国产在线视频一区二区| 色94色欧美一区二区| 成人亚洲精品一区在线观看| 日本a在线网址| a在线观看视频网站| 亚洲一卡2卡3卡4卡5卡精品中文| 久久香蕉激情| 久久精品成人免费网站| 日本黄色日本黄色录像| 国产精品久久久av美女十八| 久久精品国产a三级三级三级| 免费黄频网站在线观看国产| 午夜影院在线不卡| 少妇被粗大的猛进出69影院| 国产精品二区激情视频| 高清视频免费观看一区二区| 一二三四社区在线视频社区8| 国产欧美亚洲国产| 老司机影院毛片| 国产伦人伦偷精品视频| av网站免费在线观看视频| a级毛片在线看网站| 欧美日韩视频精品一区| 国产欧美日韩一区二区精品| bbb黄色大片| 女人精品久久久久毛片| 多毛熟女@视频| 亚洲 欧美一区二区三区| 亚洲av成人一区二区三| 欧美国产精品va在线观看不卡| 精品少妇一区二区三区视频日本电影| 国产亚洲一区二区精品| 男人舔女人的私密视频| av天堂在线播放| 国产男人的电影天堂91| 青春草亚洲视频在线观看| 精品国产国语对白av| 久热爱精品视频在线9| 亚洲av成人不卡在线观看播放网 | 中亚洲国语对白在线视频| 亚洲色图 男人天堂 中文字幕| 黄色片一级片一级黄色片| 夜夜夜夜夜久久久久| 久久午夜综合久久蜜桃| 国产高清国产精品国产三级| 天天躁日日躁夜夜躁夜夜| 在线观看www视频免费| 成年av动漫网址| 国产精品国产av在线观看| 亚洲av电影在线进入| 99香蕉大伊视频| 丁香六月欧美| 国产精品香港三级国产av潘金莲| 国产成人影院久久av| 人人妻人人添人人爽欧美一区卜| 在线观看人妻少妇| 国产精品偷伦视频观看了| 精品少妇久久久久久888优播| 91字幕亚洲| 香蕉丝袜av| 免费久久久久久久精品成人欧美视频| 啦啦啦啦在线视频资源| 精品国产乱码久久久久久男人| 免费久久久久久久精品成人欧美视频| 欧美日韩亚洲高清精品| 日韩欧美免费精品| 国产精品免费大片| 每晚都被弄得嗷嗷叫到高潮| 精品免费久久久久久久清纯 | 少妇的丰满在线观看| 亚洲国产精品成人久久小说| 日韩三级视频一区二区三区| av网站在线播放免费| 亚洲欧美日韩另类电影网站| 国产精品自产拍在线观看55亚洲 | 夜夜夜夜夜久久久久| 亚洲国产精品成人久久小说| 在线永久观看黄色视频| 制服诱惑二区| 亚洲精品国产av成人精品| 久热这里只有精品99| 亚洲五月婷婷丁香| 不卡av一区二区三区| 考比视频在线观看| 精品一区二区三区四区五区乱码| 丝袜脚勾引网站| 男女床上黄色一级片免费看| 久久久久久亚洲精品国产蜜桃av| 日韩中文字幕欧美一区二区| 午夜福利在线免费观看网站| 亚洲国产中文字幕在线视频| 欧美+亚洲+日韩+国产| 最近最新中文字幕大全免费视频| 中国国产av一级| 韩国精品一区二区三区| 亚洲人成电影观看| 91精品国产国语对白视频| 国产又爽黄色视频| 国产国语露脸激情在线看| 高清欧美精品videossex| av视频免费观看在线观看| av网站免费在线观看视频| 法律面前人人平等表现在哪些方面 | 午夜日韩欧美国产| 91字幕亚洲| 久久人妻福利社区极品人妻图片| 久久久精品国产亚洲av高清涩受| 久久久国产一区二区| 国产成人精品久久二区二区免费| 侵犯人妻中文字幕一二三四区| 国产淫语在线视频| 国产精品久久久久久人妻精品电影 | 一本一本久久a久久精品综合妖精| 纵有疾风起免费观看全集完整版| av又黄又爽大尺度在线免费看| 热re99久久国产66热| 好男人电影高清在线观看| 成在线人永久免费视频| 日本撒尿小便嘘嘘汇集6| 人妻人人澡人人爽人人| 十八禁高潮呻吟视频| 欧美日本中文国产一区发布| 制服诱惑二区| 亚洲欧洲日产国产| 无限看片的www在线观看| 午夜两性在线视频| 国产av一区二区精品久久| 国产免费现黄频在线看| 久久性视频一级片| 国产精品 欧美亚洲| 亚洲自偷自拍图片 自拍| 手机成人av网站| 97在线人人人人妻| 成人亚洲精品一区在线观看| 成人手机av| 成人影院久久| 欧美激情久久久久久爽电影 | 9191精品国产免费久久| 亚洲av成人不卡在线观看播放网 | 亚洲精品久久午夜乱码| 日韩中文字幕欧美一区二区| 国产三级黄色录像| 欧美日韩福利视频一区二区| 男女高潮啪啪啪动态图| 精品久久久精品久久久| 久久人人爽人人片av| 亚洲少妇的诱惑av| 久久精品人人爽人人爽视色| 欧美人与性动交α欧美精品济南到| 多毛熟女@视频| 色婷婷久久久亚洲欧美| 无遮挡黄片免费观看| 老司机影院毛片| 日本撒尿小便嘘嘘汇集6| 在线观看www视频免费| 日韩视频在线欧美| 国产色视频综合| 在线十欧美十亚洲十日本专区| 十八禁网站免费在线| 777米奇影视久久| 悠悠久久av| 母亲3免费完整高清在线观看| 国产成人a∨麻豆精品| 80岁老熟妇乱子伦牲交| 亚洲综合色网址| 亚洲全国av大片| 亚洲国产欧美一区二区综合| 免费日韩欧美在线观看| 国产精品av久久久久免费| 丰满迷人的少妇在线观看| 午夜影院在线不卡| 国产97色在线日韩免费| 黑人巨大精品欧美一区二区蜜桃| 午夜精品久久久久久毛片777| 色老头精品视频在线观看| 国产三级黄色录像| 成人18禁高潮啪啪吃奶动态图| 在线观看一区二区三区激情| 最黄视频免费看| 久久久久国产精品人妻一区二区| 久久人妻熟女aⅴ| av在线播放精品| 国产主播在线观看一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 纯流量卡能插随身wifi吗| 国产精品久久久久久人妻精品电影 | 日日摸夜夜添夜夜添小说| 日韩欧美一区视频在线观看| 久久久久视频综合| 999久久久国产精品视频| 最新的欧美精品一区二区| 我要看黄色一级片免费的| 国产精品九九99| 久久午夜综合久久蜜桃| 另类亚洲欧美激情| 国产精品免费大片| 美女中出高潮动态图| 少妇人妻久久综合中文| 亚洲 国产 在线| 91国产中文字幕| 国产老妇伦熟女老妇高清| 亚洲人成77777在线视频| 国产一区二区 视频在线| 久久国产精品男人的天堂亚洲| 国产麻豆69| 国产一级毛片在线| 久久精品亚洲av国产电影网| 97人妻天天添夜夜摸| 成在线人永久免费视频| 亚洲成人手机| 亚洲精品美女久久av网站| 久久久久精品国产欧美久久久 | 国产国语露脸激情在线看| 午夜两性在线视频| 汤姆久久久久久久影院中文字幕| 久久国产精品影院| 伊人亚洲综合成人网| 精品久久久精品久久久| 成在线人永久免费视频| netflix在线观看网站| 十八禁网站网址无遮挡| 亚洲欧美日韩高清在线视频 | 久久精品亚洲av国产电影网| 亚洲中文字幕日韩| 十八禁人妻一区二区| 老司机影院成人| 久久中文字幕一级| 天天躁日日躁夜夜躁夜夜| 亚洲精品粉嫩美女一区| 日日爽夜夜爽网站| 日韩,欧美,国产一区二区三区| av天堂久久9| 国产主播在线观看一区二区| 亚洲国产中文字幕在线视频| 国产主播在线观看一区二区| 欧美中文综合在线视频| 波多野结衣av一区二区av| 国产麻豆69| 国产欧美日韩一区二区三区在线| 别揉我奶头~嗯~啊~动态视频 | 久久久精品国产亚洲av高清涩受| 成人国产一区最新在线观看| 久久久国产成人免费| 国产视频一区二区在线看| 高清黄色对白视频在线免费看| 午夜福利影视在线免费观看| 中文字幕另类日韩欧美亚洲嫩草| 日韩免费高清中文字幕av| 性高湖久久久久久久久免费观看| 看免费av毛片| 久久久久精品人妻al黑| 日日夜夜操网爽| 在线精品无人区一区二区三| 久久99热这里只频精品6学生| 久久天堂一区二区三区四区| 亚洲av电影在线观看一区二区三区| av视频免费观看在线观看| 2018国产大陆天天弄谢| av欧美777| 可以免费在线观看a视频的电影网站| av天堂在线播放| bbb黄色大片| 12—13女人毛片做爰片一| 国产色视频综合| 亚洲欧美日韩另类电影网站| 超色免费av| 国产精品1区2区在线观看. | 人妻一区二区av| av天堂久久9| 他把我摸到了高潮在线观看 | 亚洲精品成人av观看孕妇| av超薄肉色丝袜交足视频| 老熟妇乱子伦视频在线观看 | 咕卡用的链子| 精品国产乱码久久久久久男人| 午夜91福利影院| 免费高清在线观看日韩| 一本综合久久免费| 肉色欧美久久久久久久蜜桃| 女性被躁到高潮视频| 欧美激情高清一区二区三区| 欧美国产精品va在线观看不卡| 天堂中文最新版在线下载| 一本久久精品| 久久99一区二区三区| 色婷婷久久久亚洲欧美| 老鸭窝网址在线观看| 精品福利观看| 亚洲精品日韩在线中文字幕| 激情视频va一区二区三区| 国产黄频视频在线观看| 婷婷色av中文字幕| 欧美精品一区二区免费开放| 国产高清视频在线播放一区 | 搡老岳熟女国产| 久久久久国产一级毛片高清牌| 日本91视频免费播放| 首页视频小说图片口味搜索| 欧美av亚洲av综合av国产av| 亚洲avbb在线观看| 亚洲视频免费观看视频| 国产人伦9x9x在线观看| 99国产精品免费福利视频| 久久久久久久精品精品| 女人久久www免费人成看片| 午夜福利视频精品| 久久人妻福利社区极品人妻图片| 免费日韩欧美在线观看| 中亚洲国语对白在线视频| 一级黄色大片毛片| 中文欧美无线码| 丝袜在线中文字幕| 国产国语露脸激情在线看| 黑人巨大精品欧美一区二区蜜桃| 成人国语在线视频| 男女无遮挡免费网站观看| 高清av免费在线| 精品人妻在线不人妻| 国产片内射在线| 一区福利在线观看| 久久热在线av| 日韩三级视频一区二区三区| 精品国产国语对白av| 久久中文字幕一级| 下体分泌物呈黄色| a在线观看视频网站| 热re99久久精品国产66热6| 久9热在线精品视频| 男人舔女人的私密视频| 久久久久久亚洲精品国产蜜桃av| 亚洲国产看品久久| 搡老乐熟女国产| 91精品三级在线观看| 欧美av亚洲av综合av国产av| 美女扒开内裤让男人捅视频| 亚洲,欧美精品.| 黄色a级毛片大全视频| 在线天堂中文资源库| 99久久综合免费| 最近中文字幕2019免费版| 免费在线观看视频国产中文字幕亚洲 | 一级毛片精品| 黑丝袜美女国产一区| 亚洲自偷自拍图片 自拍| 人成视频在线观看免费观看| 建设人人有责人人尽责人人享有的| 亚洲熟女精品中文字幕| av天堂在线播放| av又黄又爽大尺度在线免费看| h视频一区二区三区| 搡老乐熟女国产| 十八禁网站网址无遮挡| 丰满少妇做爰视频| 天天操日日干夜夜撸| 99热全是精品| 超色免费av| 少妇 在线观看| 香蕉丝袜av| 少妇猛男粗大的猛烈进出视频| 日韩欧美一区视频在线观看| 黄片大片在线免费观看| 大片免费播放器 马上看| 亚洲精品美女久久久久99蜜臀| 自拍欧美九色日韩亚洲蝌蚪91| 黄色 视频免费看| 午夜福利在线观看吧| 欧美精品啪啪一区二区三区 | 18禁黄网站禁片午夜丰满| 在线观看www视频免费| 久久天躁狠狠躁夜夜2o2o| 五月天丁香电影| 精品国产一区二区三区久久久樱花| 日日爽夜夜爽网站| 人人妻人人爽人人添夜夜欢视频| 精品一区在线观看国产| 亚洲国产看品久久| 爱豆传媒免费全集在线观看| 好男人电影高清在线观看| 国产欧美日韩一区二区精品| 亚洲中文字幕日韩| 国产成人影院久久av| 欧美久久黑人一区二区| 亚洲综合色网址| 国产在线一区二区三区精| 亚洲国产av影院在线观看| 999久久久国产精品视频| 亚洲美女黄色视频免费看| 婷婷成人精品国产| 在线观看免费午夜福利视频| 国产深夜福利视频在线观看| 国产精品秋霞免费鲁丝片| 极品少妇高潮喷水抽搐| 夜夜夜夜夜久久久久| 成年av动漫网址| 极品人妻少妇av视频| 亚洲成人免费av在线播放| 操美女的视频在线观看| 丁香六月天网| 可以免费在线观看a视频的电影网站| 国产日韩欧美亚洲二区| 精品国产一区二区久久| 亚洲国产成人一精品久久久| 亚洲av国产av综合av卡| 精品福利观看| 大型av网站在线播放| 成年人午夜在线观看视频| 成年美女黄网站色视频大全免费| 国产在视频线精品| 最近中文字幕2019免费版| 午夜日韩欧美国产| 成人亚洲精品一区在线观看| 黄色毛片三级朝国网站| 一二三四社区在线视频社区8| 我的亚洲天堂| 美女高潮喷水抽搐中文字幕| 精品一区在线观看国产| 两个人免费观看高清视频| 一级片'在线观看视频| 国产精品麻豆人妻色哟哟久久| 亚洲伊人久久精品综合| 成人国产一区最新在线观看| 三级毛片av免费| 成人国产一区最新在线观看| 国产国语露脸激情在线看| 男人舔女人的私密视频| 国产国语露脸激情在线看| 好男人电影高清在线观看| 久久精品久久久久久噜噜老黄| 999久久久国产精品视频| 在线观看免费午夜福利视频| 亚洲国产欧美日韩在线播放| 两性夫妻黄色片| 婷婷成人精品国产| 亚洲天堂av无毛| 麻豆乱淫一区二区| 亚洲精品日韩在线中文字幕| 亚洲精品一卡2卡三卡4卡5卡 | 欧美精品人与动牲交sv欧美| 9热在线视频观看99| 国产成人精品无人区| 亚洲va日本ⅴa欧美va伊人久久 | 大码成人一级视频| 日本五十路高清| 最新在线观看一区二区三区| 久久久国产欧美日韩av| 欧美黄色淫秽网站| 欧美av亚洲av综合av国产av| 91麻豆av在线| 欧美黑人精品巨大| 伦理电影免费视频| 黄片小视频在线播放| 久久久久久久久免费视频了| 亚洲精品美女久久久久99蜜臀| 国产99久久九九免费精品| 19禁男女啪啪无遮挡网站| 亚洲精品国产av成人精品| 免费观看av网站的网址| 色婷婷久久久亚洲欧美| 国产精品偷伦视频观看了| 丁香六月天网| www.熟女人妻精品国产| 国产伦理片在线播放av一区| 国产99久久九九免费精品| 男女免费视频国产| 国产成人欧美| 免费久久久久久久精品成人欧美视频| 五月开心婷婷网| 国产精品1区2区在线观看. | 久久人妻熟女aⅴ| 超碰97精品在线观看| 亚洲第一欧美日韩一区二区三区 | 国产免费福利视频在线观看| 蜜桃国产av成人99| 国产亚洲欧美精品永久| 啦啦啦 在线观看视频| 老熟妇乱子伦视频在线观看 | 亚洲精品中文字幕在线视频| 老司机福利观看| 黄色视频不卡| 狂野欧美激情性bbbbbb| 免费不卡黄色视频| 亚洲av美国av| 在线 av 中文字幕| 久久精品久久久久久噜噜老黄| 80岁老熟妇乱子伦牲交| 91九色精品人成在线观看| 久久精品亚洲av国产电影网| 免费少妇av软件| 一级毛片精品| 亚洲一区二区三区欧美精品| 91成人精品电影| 亚洲欧美日韩高清在线视频 | av一本久久久久| 亚洲九九香蕉| 丝袜在线中文字幕| 多毛熟女@视频| 最近中文字幕2019免费版| 热re99久久精品国产66热6| 国产男人的电影天堂91| 国产一区二区激情短视频 | 免费在线观看黄色视频的| av视频免费观看在线观看| 天天躁日日躁夜夜躁夜夜| 制服诱惑二区| 成年动漫av网址| 黑人操中国人逼视频| 欧美精品一区二区免费开放| 国产淫语在线视频| 9191精品国产免费久久| 日本av免费视频播放| av欧美777| 欧美日韩国产mv在线观看视频| 久久九九热精品免费| 日韩视频在线欧美| 热re99久久国产66热| 在线观看免费日韩欧美大片| 日韩视频在线欧美| 丝袜人妻中文字幕| 国产日韩欧美在线精品| 91精品国产国语对白视频| 免费高清在线观看视频在线观看| 国产日韩欧美在线精品| 国产成人a∨麻豆精品| 国产精品偷伦视频观看了| 日韩电影二区| 欧美日韩国产mv在线观看视频| 亚洲色图 男人天堂 中文字幕| 午夜福利乱码中文字幕| 黑人操中国人逼视频| 成人黄色视频免费在线看| 日韩 欧美 亚洲 中文字幕| 王馨瑶露胸无遮挡在线观看| 制服诱惑二区| 午夜福利一区二区在线看| 人妻人人澡人人爽人人| 欧美人与性动交α欧美精品济南到| 国产免费av片在线观看野外av| 亚洲自偷自拍图片 自拍| 人人妻人人澡人人看| 欧美老熟妇乱子伦牲交| a级毛片黄视频| 亚洲免费av在线视频| 国产精品九九99|