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

    基于廣義位移法的薄壁結(jié)構(gòu)超單元建模方法

    2018-12-26 07:55:30王志瑾KretovAnatolii
    關(guān)鍵詞:計(jì)算精度錐度薄壁

    原 瀟,王志瑾,夏 津,Kretov Anatolii

    (1. 南京航空航天大學(xué),飛行器先進(jìn)設(shè)計(jì)技術(shù)國防重點(diǎn)學(xué)科實(shí)驗(yàn)室,南京,210016;2. 上海機(jī)電工程研究所,上海,201109)

    0 引 言

    飛行器結(jié)構(gòu)中廣泛采用大長細(xì)比的薄壁梁結(jié)構(gòu)型式,如機(jī)身、機(jī)翼。對類似結(jié)構(gòu)進(jìn)行準(zhǔn)確的強(qiáng)度計(jì)算分析是非常復(fù)雜并且耗時(shí)耗力的事情。在計(jì)算機(jī)技術(shù)出現(xiàn)之前,要對這樣的結(jié)構(gòu)進(jìn)行分析計(jì)算,經(jīng)常采用一系列假設(shè),通過工程模型得到近似解。隨著計(jì)算技術(shù)和數(shù)值計(jì)算方法的快速發(fā)展,解析法和近似法逐漸被有限元方法取代。采用有限元方法對結(jié)構(gòu)進(jìn)行應(yīng)力應(yīng)變狀態(tài)分析前,一般先對結(jié)構(gòu)進(jìn)行必要的簡化,建立等效的有限元數(shù)學(xué)模型[1,2]。然而,有限元模型中包含大量的結(jié)構(gòu)信息數(shù)據(jù),在初始設(shè)計(jì)階段,無法給出完整的結(jié)構(gòu)信息。另外,在結(jié)構(gòu)選型階段,還需要通過初始的計(jì)算分析對結(jié)構(gòu)布置提供參考依據(jù)。因此,建立簡單、高效、可靠并有足夠計(jì)算精度的數(shù)學(xué)模型仍然是非常必要。

    在超單元法中,整體結(jié)構(gòu)被分為若干個(gè)超單元,每個(gè)超單元可以視為一系列有限元素的集合。超單元法通過對系統(tǒng)自由度進(jìn)行減縮,從而大幅減少計(jì)算量,適用于大型復(fù)雜結(jié)構(gòu)的分析[3]。Guyan[4]將結(jié)構(gòu)自由度分為主自由度和副自由度,通過強(qiáng)制忽略副自由度的影響,對剛度矩陣和質(zhì)量矩陣進(jìn)行縮聚,適用于靜力分析和低階模態(tài)分析;陳國平[5,6]將界面自由度和內(nèi)部自由度均進(jìn)行減縮,提出了適用于結(jié)構(gòu)動特性分析的超單元法;在靜力學(xué)領(lǐng)域,Ireneusz Kreja[7]提出一種梁超單元,適用于具有對稱開剖面的薄壁梁結(jié)構(gòu);康龍輝[8]提出了適用于恒定截面薄壁結(jié)構(gòu)靜力分析的超單元法,但只適用于無錐度的薄壁結(jié)構(gòu)。

    本文針對具有錐度的薄壁梁結(jié)構(gòu),基于橫截面內(nèi)剛性假設(shè),將節(jié)點(diǎn)位移縮聚為廣義位移,對于3種簡化模型分別推導(dǎo)出超單元?jiǎng)偠染仃?,通過Matlab編寫相應(yīng)求解程序。通過3個(gè)算例結(jié)果與有限元結(jié)果的對比,驗(yàn)證了該方法的準(zhǔn)確性和有效性。

    1 薄壁梁中的超單元法

    薄壁梁結(jié)構(gòu)一般由板、桿組成,可以有任意橫剖面形狀,其典型外形特征是大長細(xì)比的圓柱或錐形殼體。結(jié)構(gòu)設(shè)計(jì)條件是滿足強(qiáng)度和剛度要求下結(jié)構(gòu)重量最輕。飛機(jī)機(jī)翼、機(jī)身,飛行器機(jī)體都是具有這種特征的薄壁梁結(jié)構(gòu),它們主要由橫向構(gòu)件(翼肋或框)、縱向構(gòu)件(梁緣條、桁條)、腹板及蒙皮組成,如圖1所示。為了建立計(jì)算模型,在橫向構(gòu)件處引入與xy面平行的剖面,共(c+1)個(gè)。這樣,長度為 l的薄壁梁就由順序相連的c個(gè)單元組成,每個(gè)單元的長度為le。每個(gè)橫剖面同 z軸的交點(diǎn)就是廣義節(jié)點(diǎn),兩個(gè)橫剖面之間的一段結(jié)構(gòu)實(shí)際上就是在兩個(gè)剖面上具有廣義位移的超單元。在對超單元進(jìn)行分析時(shí),忽略了橫向構(gòu)件,只考慮桁條、蒙皮、腹板的力學(xué)行為,并將桁條(緣條)等效為n個(gè)桿單元,蒙皮、腹板等效為m個(gè)板單元。

    圖1 典型薄壁梁結(jié)構(gòu)模型Fig.1 Typical Thin-walled Beam Structure Model

    續(xù)圖1

    根據(jù)計(jì)算精度從低到高,假設(shè)每個(gè)橫剖面從絕對剛性到自由變形。假定每個(gè)剖面上的面內(nèi)廣義位移函數(shù)為:?1(z),?2(z),…,?M(z),面外廣義位移函數(shù)為:ψ1(z),ψ2(z),…,ψN(z)。其中 M、N的值根據(jù)計(jì)算精度要求選取,且Mmin=3,Mmax=m;Nmin=3,Nmax=n。剖面上任意節(jié)點(diǎn)k的位移分量可以表示為

    式中 z為與橫剖面垂直的坐標(biāo)軸;s為與橫剖面輪廓相切的方向;函數(shù)?(z),ψ(z)為尋找的廣義位移列向量,僅與坐標(biāo)z相關(guān)。橫向和縱向分布函數(shù)J(s),r(s)為廣義坐標(biāo)s的函數(shù)。

    2 超單元?jiǎng)偠染仃?/h2>

    2.1 無錐度結(jié)構(gòu)

    直翼盒模型如圖2所示。蒙皮和腹板簡化為受剪板,縱向筋條受軸向拉壓。

    圖2 直翼盒模型示意Fig.2 Diagram of Straight-wing Box Model

    模型中,M=3,N=n,即每個(gè)截面具有平面內(nèi)分別沿x軸和y軸的平動自由度?1,?2,繞z軸的轉(zhuǎn)動自由度?3以及 n個(gè)節(jié)點(diǎn)分別沿 z軸的平動自由度ψ1,ψ2,ψ3,…,ψn。相應(yīng)地,?1,?2,?3分別對應(yīng)于截面上的廣義集中力 Qx,Qy,Mz,ψ1,ψ2,ψ3,…,ψn對應(yīng)的廣義力為作用在相應(yīng)節(jié)點(diǎn)的沿z軸方向的集中力Z1,Z2,Z3,…,ZN。

    由于結(jié)構(gòu)沒有錐度,所以每根桿的長度相同,且桿軸向與整體坐標(biāo)系 z軸平行,故沿桿軸向的積分可以轉(zhuǎn)換為沿 z軸的積分,同時(shí)可以交換積分與求和的順序。

    對第i根縱向筋條:

    對第j塊板:

    式中 Sj為第j個(gè)板單元沿著s方向的長度;vje和vjb分別表示板單元末端節(jié)點(diǎn)和起始端節(jié)點(diǎn)在 z軸方向上的位移。

    寫出由n個(gè)桿單元和m個(gè)板單元組成的第e個(gè)超單元的應(yīng)變能和外力功并進(jìn)行變分,以矩陣形式表示如下:

    在一個(gè)超單元范圍內(nèi),認(rèn)為可能的位移及虛位移?(z)、ψ(z)、δ?(z)、δψ(z)符合線性變化規(guī)律,即:

    將式(9)代入式(7),根據(jù)線性位移假設(shè)進(jìn)行積分。根據(jù)能量原理:δAe=δUe,經(jīng)過整理得到:

    對于小錐度結(jié)構(gòu),需要對式(11)中的廣義力向量進(jìn)行修正:

    得到單個(gè)超單元的剛度矩陣后,按照圖3所示的方式進(jìn)行組集,即可以得到整個(gè)結(jié)構(gòu)的剛度矩陣。

    圖3 面內(nèi)剛性假設(shè)下小錐度結(jié)構(gòu)剛度矩陣組裝Fig.3 Assembly of Stiffness Matrix of Small Taper Structure under the Assumption of Cross-section Rigidity

    2.2 梯形平面正應(yīng)力板

    實(shí)際計(jì)算時(shí),當(dāng)無錐度或者結(jié)構(gòu)錐度很小時(shí),符合2.1節(jié)中對應(yīng)的情況,當(dāng)錐度較大時(shí),就必須有另外的方法。

    與此相關(guān)的主要問題就是要求解四邊形板的剛度矩陣。借助有限元中的經(jīng)典方法,用4個(gè)三角形代替四邊形,如圖4所示。

    圖4 大錐度盒段模型示意Fig.4 Diagram of Large Taper Box Segment Model

    在局部坐標(biāo)系下桿元的剛度矩陣[9]為

    在局部坐標(biāo)系下,對 4個(gè)三角形單元?jiǎng)偠染仃嚱M集,就得到四邊形單元的剛度矩陣,它有5個(gè)節(jié)點(diǎn):

    由于中心5號節(jié)點(diǎn)沒有節(jié)點(diǎn)力,根據(jù)此條件,基于式(15),可以得到5號節(jié)點(diǎn)位移與其余4個(gè)節(jié)點(diǎn)位移之間的關(guān)系,進(jìn)一步得到4節(jié)點(diǎn)四邊形平面應(yīng)力板的剛度矩陣:

    對節(jié)點(diǎn)力和位移進(jìn)行坐標(biāo)矩陣變換,就將剛度矩陣從局部坐標(biāo)變換成了總體坐標(biāo)。

    對于桿:

    對于板:

    根據(jù)上一步結(jié)果,基于節(jié)點(diǎn)力平衡原則,將 n根桿和m塊板的剛度矩陣在整體坐標(biāo)系下進(jìn)行組集,得到第e個(gè)超單元以有限元法為基礎(chǔ)的剛度矩陣,即:

    式中 ue,ue+1分別對應(yīng)于 e和 е+1剖面的位移向量;Fe,F(xiàn)e+1為相應(yīng)的節(jié)點(diǎn)力向量。

    將超單元兩端面的真實(shí)位移轉(zhuǎn)換成廣義位移:?e、ψe和?e+1、ψe+1分別為剖面 e和 e+1的廣義位移向量;Re、Ze和Re+1、Ze+1分別為剖面e和e+1的廣義力向量。

    在超單元剖面上,真實(shí)位移和廣義位移之間、廣義力和真實(shí)力之間的關(guān)系分別為

    式中 Т為位移轉(zhuǎn)置矩陣;λ為力轉(zhuǎn)置矩陣。

    且有:

    最終,對每個(gè)超單元,有以下表達(dá)式:

    式中 KF_p為對應(yīng)于節(jié)點(diǎn)力、節(jié)點(diǎn)位移的剛度矩陣;KS_p為對應(yīng)于超單元廣義力、廣義位移的剛度矩陣。

    由于每一個(gè)剖面只有 3個(gè)面內(nèi)自由度,N個(gè)面外自由度,則通過真實(shí)位移向廣義位移轉(zhuǎn)換,超單元的剛度矩陣階數(shù)由 6n×6n降為 2(3+N)×2(3+N)。

    2.3 梯形受剪板

    2.2 節(jié)中超單元?jiǎng)偠染仃嚨挠?jì)算工作量很大,對于需要多次迭代的設(shè)計(jì)和優(yōu)化過程不方便也不可取?;诖?,設(shè)想用更簡單的方法形成板的剛度矩陣,假設(shè)此時(shí)板只受剪切。

    梯形受剪板如圖5所示,其板厚為t,剪切彈性系數(shù)為G。取梯形板4個(gè)角點(diǎn)為節(jié)點(diǎn),分別記為1、2、3、4。圖中表示節(jié)點(diǎn)i的節(jié)點(diǎn)力。

    圖5 受剪梯形板剛度矩陣簡化方式Fig.5 Simplified Method of Stiffness Matrix of Shear Trapezoidal Plate

    文獻(xiàn)[9]中給出了梯形受剪板在整體坐標(biāo)系下的剛度矩陣:

    其中,

    類似于2.2節(jié),每一超單元的剛度矩陣由n根桿和m塊板的剛度矩陣按照節(jié)點(diǎn)力平衡原則進(jìn)行組集,然后通過廣義坐標(biāo)變換得到。

    式中 KF_sp為對應(yīng)于節(jié)點(diǎn)力、節(jié)點(diǎn)位移的剛度矩陣;KS_sp為對應(yīng)于超單元廣義力、廣義位移的剛度矩陣。

    3 不同模型的計(jì)算比較

    3.1 四梁盒段

    盒段的形狀、尺寸、載荷如圖6所示,材料為鋁合金,彈性模量E=72 GPa,泊松比μ=0.33。

    圖6 算例1模型示意Fig.6 Diagram of Example1 Model

    模型1、模型2和模型3是分別對應(yīng)2.1,2.2,2.3節(jié)所提出的建模方法,采用Matlab語言編寫相應(yīng)程序得到的計(jì)算結(jié)果。為了便于比較,也利用通用有限元軟件ABAQUS對該算例進(jìn)行了計(jì)算。式(28)為本算例采用虛力原理得到的解析解[10]。

    式中 第1項(xiàng)為彎曲引起的變形,第2項(xiàng)為剪切引起的變形。

    圖7給出了采用各種計(jì)算模型和方法所需的計(jì)算時(shí)間(以ABAQUS為基準(zhǔn))。采用模型1計(jì)算最快,模型3與ABAQUS相當(dāng),模型2計(jì)算時(shí)間顯著高于其它模型和方法。圖8給出采用模型1、模型2、模型3、錐度角β=0°時(shí)單個(gè)超單元長細(xì)比對計(jì)算誤差的影響。當(dāng)結(jié)構(gòu)的總長度一定時(shí),從圖8可以看出,隨著超單元個(gè)數(shù)的增加,即單個(gè)超單元長細(xì)比的減小,結(jié)構(gòu)最大位移的計(jì)算精度逐漸提高,即可通過控制單個(gè)超單元長細(xì)比(超單元個(gè)數(shù))來獲得不同的計(jì)算精度。同時(shí),模型1與模型3在錐度角β=0時(shí)計(jì)算結(jié)果的精度相當(dāng),且均低于模型2的計(jì)算精度。圖9給出了采用模型2、模型3時(shí),結(jié)構(gòu)錐度角對計(jì)算誤差的影響。從圖9中可以看出,超單元法的計(jì)算精度隨著結(jié)構(gòu)錐度角的增加而下降。

    圖7 各計(jì)算方法計(jì)算時(shí)間比較Fig.7 Comparison of Calculation Time of Each Calculation Method

    圖8 錐度角為0時(shí)計(jì)算誤差與單個(gè)超單元長細(xì)比的關(guān)系Fig.8 Relationship Curve between the Calculation Error and Long Fineness Ratio of a Single Superelement at a Taper Angle of 0°

    圖9 計(jì)算誤差與錐度角的關(guān)系Fig.9 Ralationship Curve between Calculation Error and Taper Angle

    3.2 小展弦比后掠翼

    該算例為帶4根梁和3組桁條的后掠機(jī)翼。結(jié)構(gòu)的外形、尺寸及載荷見圖10,材料為鋁合金,彈性模量E=72 GPa,泊松比μ=0.3;梁緣條面積為173 mm2,外面兩組桁條的面積為50 mm2,中間一組桁條面積為76.4 mm2,蒙皮和腹板厚度為3 mm。在結(jié)構(gòu)端部施加集中力Qy,大小為15 000 N,與y軸平行。

    圖10 后掠翼示意Fig.10 Rear Wing Diagram

    文獻(xiàn)[11]中給出了這個(gè)結(jié)構(gòu)的試驗(yàn)結(jié)果以及計(jì)算結(jié)果;模型1中,計(jì)算板剛度矩陣時(shí)采用2.2節(jié)方法,即將蒙皮和腹板當(dāng)作平面應(yīng)力板處理;模型2中,把板當(dāng)作受剪板處理,此時(shí)必須把板的面積(20~30δ2)等效到兩邊的筋條上。圖11給出了沿結(jié)構(gòu)展長方向各方法得到的位移對比。模型 1的計(jì)算結(jié)果與文獻(xiàn)[11]值大致相當(dāng),而最大位移值較文獻(xiàn)值更接近試驗(yàn)值。與試驗(yàn)值相比,模型2計(jì)算結(jié)果最大誤差為7.7%,計(jì)算精度低于模型1的最大誤差5.8%。

    圖11 沿展長方向y向位移計(jì)算結(jié)果對比Fig.11 Comparison of the Calculation Results of y-direction Displacement Along the Length of the Exhibition

    3.3 小錐度機(jī)身盒段

    該小錐度機(jī)身段長度為8 m,半徑由0.6 m過渡到0.75 m,其示意見圖12。沿機(jī)身軸向分布24根長桁,其橫截面積為50 mm2;蒙皮厚度為1.6 mm,材料為鋁合金,彈性模量 E=72 GPa,泊松比 μ=0.3。從機(jī)身固支端起沿機(jī)身長度方向,每隔0.5 m,在橫截面上施加大小為4000 N的集中力Qy,總載荷大小為64 000 N。將機(jī)身沿長度方向劃分為16個(gè)超單元。

    圖12 小錐度機(jī)身盒段示意Fig.12 Diagram of Small Taper Fuselage Box Segment

    模型1、模型2分別對應(yīng)于應(yīng)力板和受剪板兩種簡化方式。將 Matlab程序求解結(jié)果與有限元軟件ABAQUS分析結(jié)果進(jìn)行比較,沿機(jī)身長度方向橫截面y向位移結(jié)果對比見圖 13。由計(jì)算結(jié)果可知,兩種模型的計(jì)算結(jié)果與有限元分析結(jié)果相接近,模型1橫截面y向最大位移的計(jì)算誤差大于模型2,為3.2%。則超單元法亦適用于機(jī)身類薄壁結(jié)構(gòu)。沿著結(jié)構(gòu)的長度方向,超單元法的計(jì)算誤差逐漸增大。對于大長細(xì)比的薄壁結(jié)構(gòu),可以通過增加超單元的劃分個(gè)數(shù)(即減小單個(gè)超單元的長細(xì)比)來提高超單元法的計(jì)算精度。

    圖13 沿機(jī)身長度方向y向位移計(jì)算結(jié)果對比Fig.13 Comparison of the Calculation Results of y-direction Displacement Along the Length of the Fuselage

    4 結(jié) 論

    本文提出的針對大長細(xì)比薄壁結(jié)構(gòu)的超單元建模方法,單元數(shù)量少,剛度矩陣階數(shù)低,可以大大減少數(shù)據(jù)準(zhǔn)備工作量,減少人力消耗。同時(shí)超單元方法的計(jì)算精度可以根據(jù)不同位移假設(shè)依次提高,適合不同設(shè)計(jì)階段的精度要求。對于3種設(shè)計(jì)模型,2.1節(jié)中的模型和2.3節(jié)中模型計(jì)算效率高,2.2節(jié)中的模型給出的結(jié)果精度最高,但計(jì)算效率相對較低。3種模型都可以有效地用于實(shí)際結(jié)構(gòu)設(shè)計(jì)中。

    猜你喜歡
    計(jì)算精度錐度薄壁
    磁流變液仿生薄壁吸能管及其耐撞性可控度的研究
    汽車工程(2021年12期)2021-03-08 02:34:06
    次級線圈繞線錐度對LVDT靜態(tài)特性的影響
    超小型薄壁線圈架注射模設(shè)計(jì)
    模具制造(2019年7期)2019-09-25 07:29:58
    高速鋼電火花線切割錐度加工試驗(yàn)研究
    等效錐度曲線非線性特性及影響研究
    基于SHIPFLOW軟件的某集裝箱船的阻力計(jì)算分析
    廣東造船(2018年1期)2018-03-19 15:50:50
    無鍵錐度液壓聯(lián)接在大功率高轉(zhuǎn)速偶合器中的應(yīng)用
    電力與能源(2017年6期)2017-05-14 06:19:48
    一種薄壁件的加工工藝介紹
    單元類型和尺寸對拱壩壩體應(yīng)力和計(jì)算精度的影響
    鋼箱計(jì)算失效應(yīng)變的沖擊試驗(yàn)
    丁香六月欧美| 淫秽高清视频在线观看| 青草久久国产| 制服人妻中文乱码| 成人特级黄色片久久久久久久| 99国产极品粉嫩在线观看| 99精品欧美一区二区三区四区| 亚洲一码二码三码区别大吗| 中出人妻视频一区二区| 久久久久国产精品人妻aⅴ院| 老熟妇乱子伦视频在线观看| 午夜精品国产一区二区电影| 一级a爱片免费观看的视频| 男女做爰动态图高潮gif福利片 | 制服诱惑二区| 人妻丰满熟妇av一区二区三区| 亚洲精品成人av观看孕妇| 国产又色又爽无遮挡免费看| 电影成人av| 99久久综合精品五月天人人| 国产免费男女视频| 99re在线观看精品视频| 亚洲av第一区精品v没综合| 免费不卡黄色视频| 欧美性长视频在线观看| 国产又色又爽无遮挡免费看| 午夜福利在线观看吧| 久久久国产一区二区| 不卡av一区二区三区| 五月开心婷婷网| 在线观看午夜福利视频| 亚洲情色 制服丝袜| 国产精品免费视频内射| 欧美黄色片欧美黄色片| 99riav亚洲国产免费| 国内久久婷婷六月综合欲色啪| 又黄又爽又免费观看的视频| 午夜影院日韩av| 天天添夜夜摸| 国产精品亚洲av一区麻豆| 看黄色毛片网站| 精品国产亚洲在线| 久久天堂一区二区三区四区| 丝袜人妻中文字幕| 国产精品一区二区在线不卡| 动漫黄色视频在线观看| 在线天堂中文资源库| 男人的好看免费观看在线视频 | 在线免费观看的www视频| 日韩av在线大香蕉| 欧洲精品卡2卡3卡4卡5卡区| 变态另类成人亚洲欧美熟女 | 国产成人精品无人区| 久久久精品国产亚洲av高清涩受| 久久精品国产综合久久久| 大型黄色视频在线免费观看| 亚洲性夜色夜夜综合| 国产精品美女特级片免费视频播放器 | 国产亚洲欧美98| 精品高清国产在线一区| 一进一出抽搐gif免费好疼 | 亚洲精品一卡2卡三卡4卡5卡| 男人舔女人下体高潮全视频| 国产成人影院久久av| 中文字幕最新亚洲高清| 亚洲精品中文字幕在线视频| 亚洲av日韩精品久久久久久密| 黄色片一级片一级黄色片| 麻豆一二三区av精品| 国产精品综合久久久久久久免费 | 久久久精品欧美日韩精品| 无人区码免费观看不卡| 一本大道久久a久久精品| 亚洲精品在线观看二区| 90打野战视频偷拍视频| 老司机午夜福利在线观看视频| 99久久99久久久精品蜜桃| 国产欧美日韩一区二区三区在线| 国产欧美日韩精品亚洲av| 制服诱惑二区| 热re99久久国产66热| 老熟妇仑乱视频hdxx| 日韩三级视频一区二区三区| 日韩一卡2卡3卡4卡2021年| 一区二区三区国产精品乱码| 涩涩av久久男人的天堂| 国产亚洲精品久久久久5区| 亚洲人成伊人成综合网2020| 一本综合久久免费| 日本欧美视频一区| 中国美女看黄片| 丰满人妻熟妇乱又伦精品不卡| 国产精品1区2区在线观看.| 欧美黑人精品巨大| 中文欧美无线码| 国产在线精品亚洲第一网站| 日韩欧美国产一区二区入口| 免费女性裸体啪啪无遮挡网站| 国产97色在线日韩免费| 午夜日韩欧美国产| 精品国产一区二区三区四区第35| 久久精品91无色码中文字幕| 亚洲精品久久午夜乱码| 久久人人97超碰香蕉20202| 国产亚洲欧美在线一区二区| 色婷婷av一区二区三区视频| 老熟妇乱子伦视频在线观看| 深夜精品福利| 极品人妻少妇av视频| 亚洲午夜理论影院| 国产精品国产av在线观看| 电影成人av| 欧美大码av| 欧美亚洲日本最大视频资源| 国产主播在线观看一区二区| 黄频高清免费视频| 在线观看午夜福利视频| 国产极品粉嫩免费观看在线| 国产区一区二久久| 人妻久久中文字幕网| 久久伊人香网站| 成人亚洲精品av一区二区 | 欧美一区二区精品小视频在线| 欧美老熟妇乱子伦牲交| 满18在线观看网站| 中亚洲国语对白在线视频| 精品无人区乱码1区二区| 国产成人av教育| 制服人妻中文乱码| 国产片内射在线| av天堂久久9| 黄片播放在线免费| 欧美乱码精品一区二区三区| 脱女人内裤的视频| 欧美人与性动交α欧美软件| 少妇粗大呻吟视频| 久久精品国产清高在天天线| 欧美在线黄色| 99re在线观看精品视频| 亚洲一码二码三码区别大吗| 18禁黄网站禁片午夜丰满| 91精品国产国语对白视频| 国产乱人伦免费视频| 日韩欧美三级三区| 一级作爱视频免费观看| 欧美日韩国产mv在线观看视频| 午夜福利影视在线免费观看| 亚洲,欧美精品.| 国产国语露脸激情在线看| 无遮挡黄片免费观看| 中亚洲国语对白在线视频| 我的亚洲天堂| 国产高清视频在线播放一区| 51午夜福利影视在线观看| 丰满人妻熟妇乱又伦精品不卡| 丝袜在线中文字幕| 免费人成视频x8x8入口观看| 夜夜躁狠狠躁天天躁| 老熟妇仑乱视频hdxx| 欧美激情高清一区二区三区| 最新美女视频免费是黄的| 12—13女人毛片做爰片一| 亚洲性夜色夜夜综合| 成人三级做爰电影| 超碰97精品在线观看| 日韩欧美三级三区| 亚洲成国产人片在线观看| 九色亚洲精品在线播放| 日本黄色日本黄色录像| 精品久久蜜臀av无| 日韩免费av在线播放| 中文字幕另类日韩欧美亚洲嫩草| 国产精品一区二区免费欧美| 国产精品久久久av美女十八| 国产欧美日韩一区二区精品| 欧美国产精品va在线观看不卡| 一级a爱片免费观看的视频| 超色免费av| av国产精品久久久久影院| 久久人人精品亚洲av| 久久亚洲真实| 免费在线观看日本一区| 精品人妻在线不人妻| 亚洲国产精品一区二区三区在线| 一级毛片精品| 色在线成人网| 国产欧美日韩一区二区精品| 777久久人妻少妇嫩草av网站| 国产成年人精品一区二区 | 久久精品aⅴ一区二区三区四区| 欧美日韩视频精品一区| 热99国产精品久久久久久7| 亚洲,欧美精品.| aaaaa片日本免费| 脱女人内裤的视频| 丰满的人妻完整版| 亚洲免费av在线视频| 亚洲精品成人av观看孕妇| 午夜福利免费观看在线| 精品久久久久久久久久免费视频 | 午夜a级毛片| 亚洲一卡2卡3卡4卡5卡精品中文| 怎么达到女性高潮| 搡老岳熟女国产| 男人的好看免费观看在线视频 | 日韩人妻精品一区2区三区| 一级毛片女人18水好多| 日韩 欧美 亚洲 中文字幕| 一区在线观看完整版| 久久精品亚洲熟妇少妇任你| 国产亚洲欧美在线一区二区| 免费久久久久久久精品成人欧美视频| 嫩草影视91久久| 黄色视频不卡| 91大片在线观看| 村上凉子中文字幕在线| 国产成人av教育| 亚洲精品成人av观看孕妇| 天天添夜夜摸| 一级毛片精品| 亚洲av美国av| 夫妻午夜视频| 成人国产一区最新在线观看| 在线观看免费日韩欧美大片| 成年版毛片免费区| 欧美激情高清一区二区三区| 亚洲色图综合在线观看| 久久狼人影院| 最新在线观看一区二区三区| 国产无遮挡羞羞视频在线观看| 久久人妻熟女aⅴ| xxxhd国产人妻xxx| 宅男免费午夜| 日本三级黄在线观看| 国产深夜福利视频在线观看| 国产熟女午夜一区二区三区| 在线观看66精品国产| 日韩视频一区二区在线观看| 淫妇啪啪啪对白视频| 美女午夜性视频免费| 欧美精品一区二区免费开放| 麻豆成人av在线观看| 91精品三级在线观看| 老司机福利观看| 少妇粗大呻吟视频| 欧美成狂野欧美在线观看| 老司机深夜福利视频在线观看| 欧美日本中文国产一区发布| 校园春色视频在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲人成电影免费在线| 日韩中文字幕欧美一区二区| 亚洲一码二码三码区别大吗| 久久精品亚洲精品国产色婷小说| 亚洲人成77777在线视频| 天堂影院成人在线观看| 国产精品香港三级国产av潘金莲| 女人精品久久久久毛片| 久久精品91蜜桃| 精品免费久久久久久久清纯| 女人爽到高潮嗷嗷叫在线视频| 亚洲,欧美精品.| 欧美成人免费av一区二区三区| av免费在线观看网站| 大陆偷拍与自拍| 欧美久久黑人一区二区| 欧美+亚洲+日韩+国产| 麻豆av在线久日| 超碰成人久久| 久久中文字幕一级| 亚洲欧美一区二区三区久久| 国产一区二区激情短视频| 国内久久婷婷六月综合欲色啪| 久久久国产欧美日韩av| 18禁黄网站禁片午夜丰满| 亚洲熟女毛片儿| 高清黄色对白视频在线免费看| 丰满人妻熟妇乱又伦精品不卡| 在线视频色国产色| 精品久久久久久电影网| 国产欧美日韩一区二区三区在线| 国产色视频综合| 午夜激情av网站| 极品教师在线免费播放| 精品国产乱子伦一区二区三区| 久久亚洲精品不卡| 亚洲avbb在线观看| 欧美av亚洲av综合av国产av| 日本撒尿小便嘘嘘汇集6| 热re99久久精品国产66热6| 人成视频在线观看免费观看| 国产熟女xx| 久久精品影院6| 黄色丝袜av网址大全| 亚洲自拍偷在线| 国产国语露脸激情在线看| 丰满饥渴人妻一区二区三| 久久久久久久久久久久大奶| 热99国产精品久久久久久7| 成在线人永久免费视频| 亚洲av成人一区二区三| 国产精品日韩av在线免费观看 | 老汉色∧v一级毛片| 一级片免费观看大全| 色综合婷婷激情| 亚洲精品美女久久久久99蜜臀| 成人18禁在线播放| 国产精品久久久av美女十八| 亚洲欧美激情综合另类| 老熟妇仑乱视频hdxx| 丝袜人妻中文字幕| 亚洲人成网站在线播放欧美日韩| 日韩高清综合在线| 中文亚洲av片在线观看爽| 欧美成人性av电影在线观看| www日本在线高清视频| 丝袜美腿诱惑在线| 9191精品国产免费久久| 黄色视频,在线免费观看| 免费观看精品视频网站| 久久久水蜜桃国产精品网| 国产高清videossex| av在线播放免费不卡| 十分钟在线观看高清视频www| ponron亚洲| 久久九九热精品免费| 国产精品98久久久久久宅男小说| 最近最新免费中文字幕在线| 欧美精品亚洲一区二区| a级毛片在线看网站| 亚洲午夜精品一区,二区,三区| 国产伦人伦偷精品视频| 老司机深夜福利视频在线观看| 午夜福利在线免费观看网站| 欧美日韩国产mv在线观看视频| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲成国产人片在线观看| 99久久人妻综合| 女警被强在线播放| 91字幕亚洲| 午夜免费激情av| 国产精品99久久99久久久不卡| 大型av网站在线播放| 日韩免费高清中文字幕av| 日韩视频一区二区在线观看| 欧美乱码精品一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 国产精品久久久久成人av| 淫秽高清视频在线观看| 久久国产亚洲av麻豆专区| 天堂动漫精品| 窝窝影院91人妻| 中文字幕人妻丝袜一区二区| 怎么达到女性高潮| 久久香蕉国产精品| 国产1区2区3区精品| 国产黄a三级三级三级人| 50天的宝宝边吃奶边哭怎么回事| 午夜福利,免费看| 级片在线观看| 99精品欧美一区二区三区四区| 丰满迷人的少妇在线观看| 亚洲视频免费观看视频| 国产高清国产精品国产三级| 深夜精品福利| 亚洲 欧美一区二区三区| 在线天堂中文资源库| 精品国内亚洲2022精品成人| 国产一区二区在线av高清观看| 国产成人av激情在线播放| 91麻豆精品激情在线观看国产 | 色综合站精品国产| 男男h啪啪无遮挡| 日韩三级视频一区二区三区| 色综合婷婷激情| 国产成人精品在线电影| 国产成+人综合+亚洲专区| 久久中文字幕一级| 757午夜福利合集在线观看| 亚洲一码二码三码区别大吗| av免费在线观看网站| 久久精品91无色码中文字幕| 一区二区三区激情视频| 中文字幕人妻熟女乱码| 日韩精品青青久久久久久| 精品无人区乱码1区二区| 99国产极品粉嫩在线观看| 亚洲成a人片在线一区二区| 免费看a级黄色片| 黄色片一级片一级黄色片| www.自偷自拍.com| 国产精品秋霞免费鲁丝片| 97碰自拍视频| 精品无人区乱码1区二区| 亚洲伊人色综图| 亚洲免费av在线视频| 日韩成人在线观看一区二区三区| 国产成人欧美在线观看| 在线观看免费视频网站a站| 午夜成年电影在线免费观看| 午夜免费观看网址| 精品国内亚洲2022精品成人| 免费高清在线观看日韩| 无限看片的www在线观看| 日本免费a在线| 女性被躁到高潮视频| 免费搜索国产男女视频| 天堂俺去俺来也www色官网| 91国产中文字幕| 亚洲一区二区三区色噜噜 | 亚洲 欧美 日韩 在线 免费| 欧美日韩一级在线毛片| 亚洲精品一区av在线观看| 成年女人毛片免费观看观看9| 天天影视国产精品| 亚洲第一欧美日韩一区二区三区| 亚洲片人在线观看| 亚洲在线自拍视频| 国产精品日韩av在线免费观看 | 国产成人欧美| x7x7x7水蜜桃| 在线观看一区二区三区激情| 高清在线国产一区| 免费高清视频大片| 色播在线永久视频| 最近最新中文字幕大全免费视频| 黑人巨大精品欧美一区二区mp4| 高清毛片免费观看视频网站 | 日韩高清综合在线| 国产成人欧美| 日本vs欧美在线观看视频| 日本五十路高清| 久久久久国产一级毛片高清牌| av在线播放免费不卡| 免费av中文字幕在线| 国产高清激情床上av| 国产人伦9x9x在线观看| 亚洲成av片中文字幕在线观看| 久久亚洲精品不卡| 级片在线观看| 亚洲,欧美精品.| 可以在线观看毛片的网站| 国产av精品麻豆| 日韩免费高清中文字幕av| 交换朋友夫妻互换小说| 亚洲精品一卡2卡三卡4卡5卡| 1024视频免费在线观看| 日韩中文字幕欧美一区二区| 免费在线观看黄色视频的| 最好的美女福利视频网| 国产成人啪精品午夜网站| 欧美黄色片欧美黄色片| 在线国产一区二区在线| 满18在线观看网站| 久久精品91蜜桃| 热99re8久久精品国产| 波多野结衣av一区二区av| 美女扒开内裤让男人捅视频| 国产精品二区激情视频| 咕卡用的链子| 夜夜爽天天搞| 色哟哟哟哟哟哟| 性少妇av在线| 精品人妻在线不人妻| 国产精品一区二区在线不卡| 国产精品久久久人人做人人爽| 成年女人毛片免费观看观看9| 男女做爰动态图高潮gif福利片 | 99久久综合精品五月天人人| 亚洲精品一二三| 免费一级毛片在线播放高清视频 | 免费不卡黄色视频| av在线播放免费不卡| 久久国产精品人妻蜜桃| 一个人免费在线观看的高清视频| 一边摸一边抽搐一进一出视频| 91国产中文字幕| 亚洲av成人不卡在线观看播放网| 国产激情久久老熟女| 天堂动漫精品| 国产成+人综合+亚洲专区| 中亚洲国语对白在线视频| 丝袜在线中文字幕| 国产单亲对白刺激| 国产精品免费一区二区三区在线| 中文字幕高清在线视频| 亚洲七黄色美女视频| 国产精品亚洲一级av第二区| bbb黄色大片| 国产成人系列免费观看| 天堂影院成人在线观看| 午夜福利免费观看在线| 国产精品一区二区精品视频观看| 欧美日韩一级在线毛片| 久久国产精品男人的天堂亚洲| 欧美最黄视频在线播放免费 | 久久国产亚洲av麻豆专区| 18禁美女被吸乳视频| 夜夜看夜夜爽夜夜摸 | 亚洲精品国产区一区二| 女人精品久久久久毛片| 18禁国产床啪视频网站| 9191精品国产免费久久| 免费看十八禁软件| 韩国精品一区二区三区| 亚洲国产中文字幕在线视频| 亚洲av成人av| 国产精品乱码一区二三区的特点 | 黄色a级毛片大全视频| 国产成人精品在线电影| av免费在线观看网站| 免费在线观看黄色视频的| 日韩 欧美 亚洲 中文字幕| 亚洲精品国产一区二区精华液| 婷婷六月久久综合丁香| 国内久久婷婷六月综合欲色啪| 亚洲第一av免费看| 日韩欧美一区视频在线观看| 亚洲第一欧美日韩一区二区三区| 波多野结衣高清无吗| av片东京热男人的天堂| 亚洲欧美日韩无卡精品| tocl精华| 色哟哟哟哟哟哟| 国产欧美日韩一区二区精品| 亚洲三区欧美一区| 夜夜爽天天搞| 亚洲国产欧美网| 免费av毛片视频| 欧美大码av| 在线观看一区二区三区激情| 成年女人毛片免费观看观看9| 亚洲成人精品中文字幕电影 | 一进一出抽搐动态| 亚洲精品在线美女| 在线观看免费日韩欧美大片| 精品国产亚洲在线| 国产熟女午夜一区二区三区| 亚洲成人久久性| 母亲3免费完整高清在线观看| 老司机深夜福利视频在线观看| 亚洲精品粉嫩美女一区| 精品卡一卡二卡四卡免费| 久久久国产精品麻豆| 午夜91福利影院| 香蕉丝袜av| 岛国在线观看网站| a级毛片在线看网站| 侵犯人妻中文字幕一二三四区| 久久久久国内视频| 超碰97精品在线观看| 国产极品粉嫩免费观看在线| 日韩欧美国产一区二区入口| 亚洲成国产人片在线观看| 啦啦啦免费观看视频1| 欧美中文综合在线视频| 久久亚洲精品不卡| 欧美久久黑人一区二区| 国产黄a三级三级三级人| 欧美日韩国产mv在线观看视频| 一夜夜www| 国产精品一区二区精品视频观看| 99riav亚洲国产免费| 99久久久亚洲精品蜜臀av| 国产欧美日韩综合在线一区二区| 女人高潮潮喷娇喘18禁视频| 老司机午夜福利在线观看视频| 精品一品国产午夜福利视频| 久久国产精品人妻蜜桃| 欧美日本中文国产一区发布| 91九色精品人成在线观看| 狠狠狠狠99中文字幕| 亚洲人成网站在线播放欧美日韩| 久久香蕉激情| 久久人人精品亚洲av| 美女扒开内裤让男人捅视频| 男女床上黄色一级片免费看| 国产成人精品久久二区二区免费| 色精品久久人妻99蜜桃| 亚洲欧美精品综合一区二区三区| 久99久视频精品免费| 精品国产美女av久久久久小说| 久久国产精品人妻蜜桃| 美女高潮喷水抽搐中文字幕| 国产成人精品久久二区二区免费| 午夜免费激情av| 新久久久久国产一级毛片| 日本一区二区免费在线视频| 久久精品国产清高在天天线| 99国产精品一区二区三区| 欧美日韩福利视频一区二区| 久久天堂一区二区三区四区| 亚洲自偷自拍图片 自拍| 一级作爱视频免费观看| 两个人免费观看高清视频| 亚洲国产欧美网| 制服诱惑二区| 国产主播在线观看一区二区| 亚洲国产欧美网| 制服诱惑二区| 日韩欧美一区视频在线观看| av电影中文网址| 老鸭窝网址在线观看| 久久久久久人人人人人| 亚洲成国产人片在线观看| 18禁国产床啪视频网站| 亚洲精品久久午夜乱码| 国产一区二区三区综合在线观看| 国产熟女xx| 高清黄色对白视频在线免费看| 夜夜看夜夜爽夜夜摸 | 大码成人一级视频| 99re在线观看精品视频| 大码成人一级视频| 国产av在哪里看|