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

    厚薄通用四邊形平板殼元在薄壁結(jié)構(gòu)加筋布局優(yōu)化中的應(yīng)用

    2023-07-12 02:47:26丁曉紅
    空天防御 2023年2期
    關(guān)鍵詞:舵面加強(qiáng)筋四邊形

    王 謙,丁曉紅,張 橫

    (上海理工大學(xué) 機(jī)械工程學(xué)院,上海 200093)

    0 引 言

    舵面作為彈體的重要部件,是飛行器主要承力部件之一,具有輕質(zhì)、高承載的設(shè)計(jì)要求[1]。舵面結(jié)構(gòu)由薄壁板殼結(jié)構(gòu)組成,單純薄壁板殼結(jié)構(gòu)的剛度、強(qiáng)度、振動(dòng)等力學(xué)性能不佳,在實(shí)際應(yīng)用中一般采用加筋結(jié)構(gòu)形式。布置加強(qiáng)筋能夠以較小的質(zhì)量增加為代價(jià),大幅提高板殼結(jié)構(gòu)的靜動(dòng)態(tài)性能。因此,如何布局加強(qiáng)筋成為了加筋設(shè)計(jì)中的重點(diǎn)問(wèn)題。

    近年來(lái),研究人員將結(jié)構(gòu)拓?fù)鋬?yōu)化引入舵面結(jié)構(gòu)加強(qiáng)筋分布的設(shè)計(jì)中,提升了舵面結(jié)構(gòu)的各項(xiàng)性能[2-4]。但是基于變密度法(density method)尋求最優(yōu)材料布局進(jìn)而間接確定加強(qiáng)筋位置的方法需要特殊的后處理手段,才能判定加強(qiáng)筋的實(shí)際位置、走向、厚度等,但后處理過(guò)程存在一定的主觀性,會(huì)影響優(yōu)化設(shè)計(jì)結(jié)果的最優(yōu)性。丁曉紅等[5-6]提出了一種啟發(fā)式幾何顯式加筋拓?fù)鋬?yōu)化方法“自適應(yīng)成長(zhǎng)法(adaptive growth method)”,該方法基于自然界分枝系統(tǒng)成長(zhǎng)、分歧與退化機(jī)理發(fā)展而來(lái)。相較變密度法等幾何隱式方法,自適應(yīng)成長(zhǎng)法直接以加強(qiáng)筋截面幾何參數(shù)作優(yōu)化設(shè)計(jì)變量,具有優(yōu)化效率高、優(yōu)化結(jié)果清晰、可制造性強(qiáng)等優(yōu)點(diǎn),近年來(lái)被廣泛運(yùn)用于各類板殼結(jié)構(gòu)加筋靜動(dòng)態(tài)優(yōu)化設(shè)計(jì)中。

    自適應(yīng)成長(zhǎng)方法通過(guò)有限元方法進(jìn)行結(jié)構(gòu)分析,在迭代中以筋板截面參數(shù)變化表示“成長(zhǎng)”和“退化”。季學(xué)榮[7]利用自適應(yīng)成長(zhǎng)法對(duì)汽車發(fā)動(dòng)機(jī)罩板加筋布局進(jìn)行了設(shè)計(jì)。侯劍云[8]同樣受自然界生物自適應(yīng)成長(zhǎng)機(jī)理啟發(fā),提出基于連續(xù)階躍參考應(yīng)力的軟殺法(soft kill optimization,SKO)。這兩項(xiàng)研究的結(jié)構(gòu)分析部分利用商業(yè)有限元分析軟件,采用差分法計(jì)算單元靈敏度,因此,計(jì)算精度受到限制且計(jì)算效率較低;Shen 等[9-10]基于自適應(yīng)成長(zhǎng)法提出了3D 箱型結(jié)構(gòu)的剛度、阻尼協(xié)同優(yōu)化方法,結(jié)構(gòu)分析部分利用Long 等[11]提出的基于廣義協(xié)調(diào)理論的平板型矩形殼元GCR24 處理的,該單元由平面薄膜元GR12 和薄板彎曲元GPL-R12 單元耦合構(gòu)建。GCR24 單元用于幾何顯式板殼加筋拓?fù)鋬?yōu)化時(shí)具有列式簡(jiǎn)單、計(jì)算效率高的優(yōu)勢(shì),同時(shí)在筋板厚度變化過(guò)程中避免了剪切閉鎖問(wèn)題,但GCR24 單元是矩形單元,不能滿足舵面等復(fù)雜結(jié)構(gòu)的有限元分析對(duì)網(wǎng)格靈活性的需求。

    綜上所述,利用自適應(yīng)成長(zhǎng)法對(duì)舵面結(jié)構(gòu)的加筋布局進(jìn)行拓?fù)鋬?yōu)化,構(gòu)建一種厚薄通用、高準(zhǔn)確度的四邊形平板殼元十分必要。本文通過(guò)耦合GQ12膜元和MITC4 板元構(gòu)造了四邊形平板殼元并構(gòu)建基結(jié)構(gòu),隨后引入自適應(yīng)成長(zhǎng)法對(duì)舵面結(jié)構(gòu)進(jìn)行加筋布局拓?fù)鋬?yōu)化,以此實(shí)現(xiàn)結(jié)構(gòu)的輕量化設(shè)計(jì)。

    1 GQ12 膜元與MITC4 板元耦合的四邊形平板殼元

    1.1 具有旋轉(zhuǎn)自由度的四邊形膜元GQ12

    為克服低階位移型單元計(jì)算精度差的問(wèn)題,Long等[12-13]采用廣義協(xié)調(diào)條件構(gòu)造了具有旋轉(zhuǎn)自由度的四邊形膜元GQ12,該單元列式簡(jiǎn)單,能通過(guò)任意四邊形分片檢驗(yàn)收斂單元,計(jì)算精度較高。GQ12引入剛體轉(zhuǎn)角作為單元節(jié)點(diǎn)的旋轉(zhuǎn)自由度,每個(gè)節(jié)點(diǎn)有兩個(gè)線自由度和一個(gè)平面內(nèi)轉(zhuǎn)角自由度,每個(gè)節(jié)點(diǎn)的自由度為

    式中:ui、vi為節(jié)點(diǎn)線位移分量;θi為節(jié)點(diǎn)附加平面內(nèi)轉(zhuǎn)角位移分量。

    GQ12膜元的位移場(chǎng)由兩個(gè)部分組成:

    式中:u0=[u0v0]T,為雙線性協(xié)調(diào)位移場(chǎng);uθ=[uθvθ]T,為由單元節(jié)點(diǎn)角位移分量θi引起的附加位移場(chǎng)。附加位移場(chǎng)為

    式中:(ξ,η)為單元的自然坐標(biāo)系。

    將附加剛體轉(zhuǎn)角位移而產(chǎn)生的單元邊界位移沿單元邊界設(shè)為三次分布,并對(duì)每一邊引入廣義協(xié)調(diào)條件。待定參數(shù)αi、βi(i=1,2,3,4)的詳細(xì)推導(dǎo)過(guò)程參見文獻(xiàn)[12],代入式(3)可得GQ12 單元的總位移場(chǎng),表示為

    式中:N0i、Nuθi、Nvθi為單元對(duì)不同位移場(chǎng)的形函數(shù)。

    式中:ξ、η和ξi、ηi分別為單元和單元節(jié)點(diǎn)的自然坐標(biāo)系坐標(biāo)。

    GQ12的應(yīng)變場(chǎng)為

    單元?jiǎng)偠染仃嚍?/p>

    式中:B為單元應(yīng)力-應(yīng)變矩陣;D為彈性系數(shù)矩陣;|J|為雅可比矩陣行列式;t為單元厚度。

    1.2 剪應(yīng)變混合插值板元MITC4

    基于一階剪切變形理論的Mindlin 板元對(duì)位移場(chǎng)和轉(zhuǎn)角場(chǎng)各自獨(dú)立插值,由于單元邊界只要求C0類連續(xù),因此滿足網(wǎng)格靈活性的需求,以及舵面等復(fù)雜外形結(jié)構(gòu)的分析需要。一階剪切變形理論放棄了Kirchhoff 理論的直法線假設(shè),因此對(duì)于橫向剪切變形明顯的中厚板的分析,Mindlin 板元的預(yù)測(cè)效果相比Kirchhoff板元更好。但是,Mindlin板用于薄板時(shí)存在剪切閉鎖和應(yīng)力振蕩等問(wèn)題,不能滿足幾何顯式加筋拓?fù)鋬?yōu)化對(duì)有限元分析的要求。

    縮減積分和選擇積分方法列式簡(jiǎn)單,可解決剪切閉鎖的問(wèn)題,但根據(jù)約束類型等的不同,會(huì)帶來(lái)零能模式問(wèn)題?;旌喜逯凳且环N比較理想的構(gòu)造厚薄通用板元的方法,通過(guò)對(duì)單元面內(nèi)應(yīng)變和橫向剪切應(yīng)變使用不同的插值函數(shù),避免剪切閉鎖和零能模式問(wèn)題。Bathe等[14-15]基于Mindlin板理論提出剪應(yīng)變混合插值板元(mixed interpolated tensorial component,MITC4),計(jì)算結(jié)果較好,滿足復(fù)雜結(jié)構(gòu)的工程分析需要,如圖1所示。

    圖1 板單元位移場(chǎng)Fig.1 Displacement field of plate element

    由圖1可得,板單元的橫向剪切應(yīng)變?yōu)?/p>

    式中:θx、θy為x軸和y軸轉(zhuǎn)角位移分量;ω為橫向位移。

    Mindlin板理論中,對(duì)于系統(tǒng)總位能公式中的彎曲應(yīng)變能項(xiàng),MITC4 單元使用與基于位移方法的Mindlin板單元相同的插值方法,廣義位移通過(guò)相同的形函數(shù)獨(dú)立插值為

    式中:ωi、θxi、θyi分別為單元各節(jié)點(diǎn)處橫向位移分量、x軸和y軸轉(zhuǎn)角位移分量。

    但是,對(duì)于橫向剪切應(yīng)變能項(xiàng),如果采用相同的插值形式,會(huì)在薄板情況下出現(xiàn)剪切閉鎖問(wèn)題。Bathe 等[14-15]提出的MITC4 單元方案中,構(gòu)造了一種在單元邊和單元內(nèi)滿足剪應(yīng)變一致的插值函數(shù),插值形式為

    式中:γAξz、γBηz、γCξz、γDηz分別為MITC4 單元在ξ-η-z坐標(biāo)系下A、B、C、D四邊中點(diǎn)的橫向剪應(yīng)變。

    如圖2 所示,對(duì)于x-y平面內(nèi)任意凸四邊形單元I-J-K-L,γξz、γηz的計(jì)算式為

    圖2 任意四邊形單元Fig.2 Quadrilateral element

    式中:detJ為雅可比矩陣行列式,

    則有

    同時(shí),由于剪應(yīng)變表達(dá)式定義在ξ-η-z坐標(biāo)系下,需要轉(zhuǎn)化到x-y-z坐標(biāo)系下,轉(zhuǎn)換形式為

    式中:α為ξ軸與x軸夾角;β為η軸與y軸夾角。

    1.3 驗(yàn)證算例

    為驗(yàn)證GQ12 膜元與MITC4 板元耦合平板殼元的分析精度,設(shè)計(jì)一個(gè)四邊形折板,其幾何尺寸如圖3所示。折板一側(cè)約束,另一側(cè)一點(diǎn)受集中力載荷F=-1 000 N 的作用,取結(jié)構(gòu)E=210 GPa,v=0.3,分別對(duì)5 mm厚度及10 mm厚度折板算例求解結(jié)構(gòu)應(yīng)變。

    圖3 四邊形折板算例Fig.3 Example of quadrilateral folded-plate

    5 mm 厚度折板GQ12 膜元與MITC4 板元耦合平板殼元自編Matlab 程序的分析結(jié)果如圖4(a)所示,采用商業(yè)有限元分析軟件Ansys的分析結(jié)果如圖4(b)所示。不同厚度折板加載節(jié)點(diǎn)和最大位移節(jié)點(diǎn)Z軸方向位移數(shù)值對(duì)照見表1和表2。

    表1 5 mm厚四邊形折板Z軸方向位移Tab.1 Z-axis displacement of the quadrilateral foldedplate with 5 mm thickness

    表2 10 mm厚四邊形折板Z軸方向位移Tab.2 Z-axis displacement of the quadrilateral foldedplate with 10 mm thickness

    圖4 四邊形折板分析結(jié)果對(duì)比Fig.4 Analysis results comparison of the quadrilateral folded-plate

    對(duì)以上算例結(jié)果進(jìn)行對(duì)比,考察加載節(jié)點(diǎn)與最大位移節(jié)點(diǎn)Z軸方向位移值,相同約束和載荷條件下,5 mm 厚度折板與10 mm 厚度折板自編程序GQ12+MITC4 耦合平板殼元與商業(yè)有限元軟件分析結(jié)果誤差均<4.5%,說(shuō)明GQ12+MITC4耦合平板殼元具備厚薄通用性,計(jì)算精度良好。

    2 自適應(yīng)成長(zhǎng)法

    在自適應(yīng)成長(zhǎng)法結(jié)構(gòu)加筋拓?fù)鋬?yōu)化中,采用基結(jié)構(gòu)法建立優(yōu)化幾何模型,其優(yōu)點(diǎn)在于避免優(yōu)化迭代過(guò)程中網(wǎng)格的重劃分。參考承載和約束條件確定“種子點(diǎn)”的位置,如圖5(a)所示。從“種子點(diǎn)”生長(zhǎng)出來(lái)的加強(qiáng)筋根據(jù)一定規(guī)則成長(zhǎng)或退化;生長(zhǎng)到一定尺度可以“分歧”,再?gòu)摹胺制琰c(diǎn)”繼續(xù)生長(zhǎng),在下步迭代中參與成長(zhǎng)與退化;退化到一定尺度則“消失”,如圖5(b)—圖5(c)所示。這個(gè)過(guò)程在優(yōu)化迭代中反復(fù)進(jìn)行,直至滿足收斂條件。

    圖5 加強(qiáng)筋拓?fù)鋬?yōu)化的自適應(yīng)成長(zhǎng)過(guò)程Fig.5 Adaptive growth process of stiffener topology optimization

    2.1 數(shù)學(xué)模型及靈敏度分析

    靜載荷下,薄壁結(jié)構(gòu)自適應(yīng)成長(zhǎng)法加筋拓?fù)鋬?yōu)化問(wèn)題的數(shù)學(xué)模型可以表述為:在給定的體積約束下,選擇加強(qiáng)筋厚度作為設(shè)計(jì)變量,以結(jié)構(gòu)在載荷條件下的總應(yīng)變能為目標(biāo)函數(shù),逐步尋求最優(yōu)解。其統(tǒng)一的數(shù)學(xué)模型可以表達(dá)為

    式中:X為設(shè)計(jì)變量;xi為第i個(gè)活動(dòng)筋板單元的厚度;n為加強(qiáng)筋單元的總數(shù);Φ(X)為目標(biāo)函數(shù),即總應(yīng)變能;V、V0分別為整體結(jié)構(gòu)總體積和結(jié)構(gòu)初始體積;η為體積約束因子;xmin、xmax為設(shè)計(jì)變量xi的取值下限和上限;g(X)為約束函數(shù)。

    薄壁結(jié)構(gòu)的柔度關(guān)于加強(qiáng)筋厚度的靈敏度計(jì)算公式為

    式中:ui為第i加強(qiáng)筋單元的位移矢量;ki為第i加強(qiáng)筋單元的剛度矩陣。

    2.2 優(yōu)化設(shè)計(jì)流程

    基于優(yōu)化數(shù)學(xué)模型和自適應(yīng)成長(zhǎng)法基本原理,以移動(dòng)漸近線法(method of moving asymptotes,MMA)作為迭代更新準(zhǔn)則對(duì)結(jié)構(gòu)分布進(jìn)行尋優(yōu),具體流程如圖6所示。

    圖6 自適應(yīng)成長(zhǎng)法優(yōu)化設(shè)計(jì)流程Fig.6 Flow chart of optimization design of adaptive growth method

    1) 建立基結(jié)構(gòu),設(shè)置設(shè)計(jì)域。

    2) 選取“種子點(diǎn)”,設(shè)置優(yōu)化參數(shù)。參考結(jié)構(gòu)的載荷和支撐情況,選取基結(jié)構(gòu)上若干節(jié)點(diǎn)作為“種子點(diǎn)”,給定加強(qiáng)筋的初始厚度x0、分歧臨界值xb、體積約束因子η、收斂容差ε以及最大迭代次數(shù)N。

    3) 靈敏度分析。對(duì)結(jié)構(gòu)進(jìn)行有限元分析,通過(guò)解析法計(jì)算活動(dòng)筋板單元的靈敏度。

    4) 成長(zhǎng)和退化。判斷收斂狀態(tài),如果滿足筋板成長(zhǎng)條件,通過(guò)MMA 算法更新xi;如果筋板厚度達(dá)到分歧臨界值xb,該單元進(jìn)行分歧;如果滿足退化條件,該單元退化。

    3 舵面結(jié)構(gòu)加筋布局優(yōu)化

    3.1 舵面基結(jié)構(gòu)構(gòu)建

    根據(jù)舵面幾何形狀,在優(yōu)化設(shè)計(jì)前構(gòu)建拓?fù)鋬?yōu)化基結(jié)構(gòu),舵面外層為蒙皮結(jié)構(gòu),定義厚度為0.5 mm。選擇舵面內(nèi)部空間區(qū)域作為設(shè)計(jì)域,定義基結(jié)構(gòu)初始厚度為0.5 mm,如圖7 所示。由于舵面的幾何特征,構(gòu)成基結(jié)構(gòu)的單元為任意凸四邊形單元,GQ12 膜元和MITC4 板元耦合構(gòu)造的四邊形等參平板殼元適用于舵面基結(jié)構(gòu)的有限元分析。

    圖7 舵面優(yōu)化初始基結(jié)構(gòu)Fig.7 Initial ground structure of rudder structure for optimization

    3.2 基于自適應(yīng)成長(zhǎng)法的加筋布局

    根據(jù)舵面的實(shí)際工作狀態(tài),對(duì)其約束施加于舵軸,受單側(cè)橫向均布載荷。因此,本文以舵面舵軸處約束點(diǎn)為自適應(yīng)成長(zhǎng)法的“種子點(diǎn)”,以整體結(jié)構(gòu)應(yīng)變能最小為設(shè)計(jì)目標(biāo),骨架結(jié)構(gòu)體積為約束條件,對(duì)舵面進(jìn)行加筋拓?fù)鋬?yōu)化,尋找最優(yōu)的骨架結(jié)構(gòu)分布。優(yōu)化的數(shù)學(xué)模型為

    式中:T為設(shè)計(jì)變量,為每個(gè)活動(dòng)筋板的厚度;U(T)為目標(biāo)函數(shù),即整體結(jié)構(gòu)應(yīng)變能;V、V0分別為加筋舵面結(jié)構(gòu)總體積和初始實(shí)體舵面結(jié)構(gòu)體積;η=20%為體積約束因子。

    經(jīng)過(guò)拓?fù)鋬?yōu)化設(shè)計(jì)得到的舵面內(nèi)部最優(yōu)加筋布局如圖8 所示。圖中呈現(xiàn)出清晰的分枝結(jié)構(gòu),拓?fù)湫螒B(tài)合理且具有較優(yōu)的可制造性。

    圖8 拓?fù)鋬?yōu)化后的舵面加筋布局Fig.8 Rudder stiffener layout obtained by topology optimization

    優(yōu)化過(guò)程迭代曲線如圖9 所示。優(yōu)化迭代過(guò)程中,應(yīng)變能與初始應(yīng)變能比值在前10個(gè)迭代步中下降較快,隨后迭代過(guò)程逐漸平緩直至收斂,優(yōu)化效果理想,驗(yàn)證了GQ12膜元和MITC4板元耦合四邊形等參平板殼元分析的準(zhǔn)確性,自適應(yīng)成長(zhǎng)法的有效性及優(yōu)化結(jié)構(gòu)的合理性。

    圖9 舵面結(jié)構(gòu)加筋布局優(yōu)化迭代曲線Fig.9 Iteration history curve of rudder stiffener layout

    3.3 設(shè)計(jì)結(jié)果分析

    為了評(píng)估加筋拓?fù)鋬?yōu)化效果,利用耦合四邊形平板殼元建立舵面的傳統(tǒng)蒙皮-網(wǎng)格式骨架設(shè)計(jì)模型,如圖10所示。對(duì)文獻(xiàn)[1]中介紹的傳統(tǒng)骨架設(shè)計(jì)進(jìn)行建模并分析,計(jì)算得到拓?fù)鋬?yōu)化前后舵面的整體結(jié)構(gòu)應(yīng)變能,如圖11所示。

    圖10 舵面?zhèn)鹘y(tǒng)骨架設(shè)計(jì)Fig.10 Traditional rudder skeleton design

    圖11 舵面結(jié)構(gòu)優(yōu)化前后整體結(jié)構(gòu)應(yīng)變Fig.11 Strain of overall structure before and after rudder structure optimization

    對(duì)比傳統(tǒng)網(wǎng)格式骨架設(shè)計(jì)與加筋拓?fù)鋬?yōu)化后舵面整體結(jié)構(gòu)的應(yīng)變能數(shù)值以及結(jié)構(gòu)質(zhì)量,結(jié)果見表3。由表3可見,優(yōu)化后結(jié)構(gòu)的應(yīng)變能為17.28 J,較傳統(tǒng)設(shè)計(jì)25.80 J降低了33%,結(jié)構(gòu)剛度得到增強(qiáng);傳統(tǒng)骨架設(shè)計(jì)的質(zhì)量為1.19 kg,優(yōu)化后結(jié)構(gòu)的質(zhì)量降低了29%。

    表3 舵面結(jié)構(gòu)優(yōu)化前后性能比較Tab.3 Performance comparison before and after rudder structure optimization

    4 結(jié)束語(yǔ)

    本文利用GQ12 膜元與MITC4 板元耦合構(gòu)造了四邊形等參平板殼元,通過(guò)有限元分析實(shí)現(xiàn)了基于自適應(yīng)成長(zhǎng)法的不規(guī)則薄壁結(jié)構(gòu)的加筋設(shè)計(jì)。隨后以輕質(zhì)、高承載的舵面結(jié)構(gòu)為例,獲得了靜載荷狀態(tài)下舵面內(nèi)部的最優(yōu)骨架布局。經(jīng)過(guò)對(duì)優(yōu)化前后舵面結(jié)構(gòu)的性能對(duì)比,驗(yàn)證了GQ12膜元與MITC4板元耦合平板殼元分析工具的準(zhǔn)確性和實(shí)用性,同時(shí)驗(yàn)證了自適應(yīng)成長(zhǎng)法對(duì)舵面優(yōu)化設(shè)計(jì)的有效性。

    猜你喜歡
    舵面加強(qiáng)筋四邊形
    基于舵面振蕩激勵(lì)試飛的飛機(jī)動(dòng)特性識(shí)別與修正
    基于拔銷器鎖定的飛行器氣動(dòng)控制面解鎖控制方法
    圓錐曲線內(nèi)接四邊形的一個(gè)性質(zhì)
    汽車零件加強(qiáng)筋結(jié)構(gòu)設(shè)計(jì)
    汽車工程師(2019年7期)2019-08-12 07:25:04
    越野車側(cè)壁加強(qiáng)筋板工藝改進(jìn)
    模具制造(2019年3期)2019-06-06 02:10:52
    一種折疊舵面設(shè)計(jì)及分析
    四邊形逆襲記
    4.4 多邊形和特殊四邊形
    加強(qiáng)筋的截面尺寸對(duì)注塑制品翹曲變形的影響
    塑料制造(2016年5期)2016-06-15 20:27:39
    多舵面飛機(jī)電力作動(dòng)系統(tǒng)協(xié)調(diào)控制策略研究
    kizo精华| 国产精品 国内视频| 波野结衣二区三区在线| 成人手机av| 国产97色在线日韩免费| 男男h啪啪无遮挡| 少妇 在线观看| 国产在线一区二区三区精| 另类亚洲欧美激情| 超色免费av| 又粗又硬又长又爽又黄的视频| 国产av国产精品国产| 亚洲三区欧美一区| 精品一品国产午夜福利视频| 伊人亚洲综合成人网| 可以免费在线观看a视频的电影网站 | 亚洲精品aⅴ在线观看| 亚洲综合色惰| 一本大道久久a久久精品| 熟妇人妻不卡中文字幕| 看十八女毛片水多多多| 日韩制服骚丝袜av| 亚洲综合精品二区| 日日啪夜夜爽| 国产精品.久久久| 国产精品成人在线| 久久青草综合色| 韩国高清视频一区二区三区| 久久精品国产亚洲av涩爱| 最新的欧美精品一区二区| 国产 精品1| tube8黄色片| 亚洲人成77777在线视频| 一区在线观看完整版| 国产一区二区 视频在线| 少妇人妻精品综合一区二区| 两个人免费观看高清视频| 满18在线观看网站| 一级毛片电影观看| 热re99久久精品国产66热6| 欧美老熟妇乱子伦牲交| 肉色欧美久久久久久久蜜桃| 日韩大片免费观看网站| 成年人午夜在线观看视频| 午夜影院在线不卡| 午夜福利网站1000一区二区三区| √禁漫天堂资源中文www| 老司机影院成人| 国产一区有黄有色的免费视频| 欧美最新免费一区二区三区| 亚洲精品乱久久久久久| 免费在线观看视频国产中文字幕亚洲 | 青草久久国产| 免费观看性生交大片5| 999精品在线视频| 99re6热这里在线精品视频| 国产精品偷伦视频观看了| 亚洲精品国产色婷婷电影| 国产精品免费视频内射| 国产毛片在线视频| 国产欧美日韩综合在线一区二区| 日日爽夜夜爽网站| 97在线人人人人妻| 中文天堂在线官网| 制服人妻中文乱码| 丝袜在线中文字幕| 欧美最新免费一区二区三区| 欧美日韩亚洲国产一区二区在线观看 | 精品福利永久在线观看| 丰满迷人的少妇在线观看| 久久人人爽人人片av| 亚洲欧美成人精品一区二区| 日本av免费视频播放| 色视频在线一区二区三区| 啦啦啦啦在线视频资源| 成人免费观看视频高清| 少妇的丰满在线观看| 国产精品久久久久成人av| 精品亚洲成国产av| 精品福利永久在线观看| 欧美日韩视频高清一区二区三区二| 性色av一级| 国产成人一区二区在线| 亚洲伊人久久精品综合| 99九九在线精品视频| 亚洲精品一区蜜桃| 国产成人精品无人区| 大片电影免费在线观看免费| 九色亚洲精品在线播放| 午夜免费男女啪啪视频观看| 中文字幕人妻丝袜制服| 最黄视频免费看| 亚洲国产av影院在线观看| 一区二区av电影网| 日韩av免费高清视频| 毛片一级片免费看久久久久| 97人妻天天添夜夜摸| 一级a爱视频在线免费观看| 美女视频免费永久观看网站| 国产一区二区三区综合在线观看| 少妇人妻 视频| 女人久久www免费人成看片| 少妇被粗大猛烈的视频| 日韩免费高清中文字幕av| 欧美av亚洲av综合av国产av | 成人国产av品久久久| 国产成人91sexporn| 精品久久久精品久久久| kizo精华| 视频区图区小说| 丝袜美足系列| 久久狼人影院| 成人黄色视频免费在线看| 亚洲 欧美一区二区三区| 久久久久久久国产电影| 午夜福利在线观看免费完整高清在| 又粗又硬又长又爽又黄的视频| 狂野欧美激情性bbbbbb| 高清视频免费观看一区二区| av免费观看日本| 国产成人a∨麻豆精品| av在线app专区| 亚洲第一区二区三区不卡| 国产男人的电影天堂91| 亚洲欧美一区二区三区国产| 岛国毛片在线播放| 大片电影免费在线观看免费| 国产精品蜜桃在线观看| 午夜福利在线观看免费完整高清在| 日韩av不卡免费在线播放| av免费在线看不卡| 亚洲精品第二区| 男人舔女人的私密视频| 亚洲一码二码三码区别大吗| 美女xxoo啪啪120秒动态图| 飞空精品影院首页| 亚洲欧美一区二区三区黑人 | 在线天堂中文资源库| 成人国产麻豆网| 一本大道久久a久久精品| 亚洲精品久久成人aⅴ小说| 亚洲人成电影观看| 咕卡用的链子| 熟女av电影| 久久久久精品久久久久真实原创| 成人毛片60女人毛片免费| av又黄又爽大尺度在线免费看| 国产精品人妻久久久影院| a 毛片基地| 免费少妇av软件| 国产白丝娇喘喷水9色精品| 亚洲成人av在线免费| 在线观看免费日韩欧美大片| 一区在线观看完整版| 人人妻人人澡人人看| 爱豆传媒免费全集在线观看| 国产精品一二三区在线看| 熟女av电影| 日产精品乱码卡一卡2卡三| 国产精品无大码| 日本黄色日本黄色录像| 欧美激情高清一区二区三区 | 18禁动态无遮挡网站| 黑人巨大精品欧美一区二区蜜桃| 精品第一国产精品| 美女xxoo啪啪120秒动态图| 99久久综合免费| 国产无遮挡羞羞视频在线观看| 国产精品久久久av美女十八| kizo精华| 国产福利在线免费观看视频| 一区二区av电影网| 国产精品一区二区在线不卡| 人妻人人澡人人爽人人| 久热久热在线精品观看| 欧美97在线视频| 亚洲人成网站在线观看播放| 国产一区二区激情短视频 | 欧美精品人与动牲交sv欧美| 亚洲欧洲国产日韩| 欧美bdsm另类| 久久久久视频综合| 少妇人妻 视频| 国产成人免费无遮挡视频| 老司机亚洲免费影院| 纯流量卡能插随身wifi吗| 日韩电影二区| 99精国产麻豆久久婷婷| 精品午夜福利在线看| 麻豆乱淫一区二区| 亚洲第一区二区三区不卡| 91久久精品国产一区二区三区| 999精品在线视频| 国产熟女欧美一区二区| 国产成人精品婷婷| www.自偷自拍.com| 天天操日日干夜夜撸| 日本爱情动作片www.在线观看| 午夜福利视频在线观看免费| 免费观看a级毛片全部| 99久久综合免费| 80岁老熟妇乱子伦牲交| 亚洲一级一片aⅴ在线观看| 欧美成人午夜精品| 在线看a的网站| 啦啦啦中文免费视频观看日本| 国产一区有黄有色的免费视频| 只有这里有精品99| 丰满少妇做爰视频| 人人妻人人爽人人添夜夜欢视频| www.精华液| 男人舔女人的私密视频| 一级片免费观看大全| 国精品久久久久久国模美| 男女无遮挡免费网站观看| 69精品国产乱码久久久| 亚洲欧美一区二区三区国产| 中文字幕人妻丝袜制服| 日韩成人av中文字幕在线观看| 国产亚洲av片在线观看秒播厂| 午夜日本视频在线| 久久午夜福利片| 午夜免费男女啪啪视频观看| 国产男人的电影天堂91| 国产乱人偷精品视频| 国产极品天堂在线| av不卡在线播放| 亚洲,欧美,日韩| 一本色道久久久久久精品综合| 青春草国产在线视频| 99热国产这里只有精品6| 伦理电影免费视频| 久久国产精品男人的天堂亚洲| 91aial.com中文字幕在线观看| 久久国内精品自在自线图片| 日本黄色日本黄色录像| 国产精品三级大全| 99久久人妻综合| 国产黄色免费在线视频| 丝袜脚勾引网站| 国产白丝娇喘喷水9色精品| 男女国产视频网站| 成年女人在线观看亚洲视频| 中国三级夫妇交换| 久热久热在线精品观看| av在线观看视频网站免费| 午夜福利影视在线免费观看| 国产精品.久久久| 国产视频首页在线观看| 精品一区二区三卡| 一本色道久久久久久精品综合| 黄色 视频免费看| 母亲3免费完整高清在线观看 | 亚洲在久久综合| 日韩av免费高清视频| 精品少妇一区二区三区视频日本电影 | 精品午夜福利在线看| 国产精品熟女久久久久浪| 天堂俺去俺来也www色官网| 一级爰片在线观看| 男女无遮挡免费网站观看| 国产午夜精品一二区理论片| 精品人妻偷拍中文字幕| 亚洲欧洲国产日韩| 9191精品国产免费久久| 丝袜在线中文字幕| 熟女av电影| 亚洲av日韩在线播放| 成年女人毛片免费观看观看9 | 亚洲国产精品成人久久小说| 宅男免费午夜| 色哟哟·www| 亚洲av在线观看美女高潮| 亚洲美女黄色视频免费看| 色婷婷久久久亚洲欧美| 狠狠婷婷综合久久久久久88av| 最近中文字幕高清免费大全6| 人妻 亚洲 视频| 中文字幕人妻熟女乱码| av国产久精品久网站免费入址| 亚洲国产毛片av蜜桃av| 在线观看三级黄色| 在线观看免费高清a一片| 天天躁夜夜躁狠狠躁躁| 国产免费又黄又爽又色| 满18在线观看网站| 国产精品秋霞免费鲁丝片| av电影中文网址| 香蕉精品网在线| 亚洲国产精品999| 男女边摸边吃奶| 国产成人91sexporn| 99九九在线精品视频| 久久狼人影院| 午夜av观看不卡| 大片免费播放器 马上看| 国产野战对白在线观看| 久久久国产精品麻豆| 国产 一区精品| 精品少妇久久久久久888优播| 亚洲精品中文字幕在线视频| 欧美黄色片欧美黄色片| 1024视频免费在线观看| 欧美少妇被猛烈插入视频| 日韩一区二区三区影片| 秋霞伦理黄片| 国产成人精品久久久久久| 亚洲av福利一区| 热re99久久国产66热| 欧美97在线视频| 亚洲av.av天堂| 久久久久久久国产电影| 一区在线观看完整版| 国产精品国产三级国产专区5o| 在线看a的网站| 99久久综合免费| 在线天堂中文资源库| 国产一区二区 视频在线| av视频免费观看在线观看| av又黄又爽大尺度在线免费看| 国产乱人偷精品视频| 精品国产乱码久久久久久小说| 国产精品一二三区在线看| 久久久久久久久久久久大奶| 国产老妇伦熟女老妇高清| 99精国产麻豆久久婷婷| 美国免费a级毛片| 成人免费观看视频高清| 国产精品久久久久久av不卡| 美女高潮到喷水免费观看| 国产综合精华液| 18禁国产床啪视频网站| 丰满迷人的少妇在线观看| 午夜免费男女啪啪视频观看| 伦理电影免费视频| 宅男免费午夜| 国产精品女同一区二区软件| 欧美+日韩+精品| 国产日韩一区二区三区精品不卡| 尾随美女入室| 青青草视频在线视频观看| 欧美老熟妇乱子伦牲交| 人人澡人人妻人| 国产免费一区二区三区四区乱码| 国产在线免费精品| 美女福利国产在线| 成人亚洲欧美一区二区av| 中文字幕av电影在线播放| 久久久久精品人妻al黑| 亚洲色图综合在线观看| 成人免费观看视频高清| 午夜91福利影院| 99re6热这里在线精品视频| 亚洲精品乱久久久久久| 国产乱来视频区| 免费高清在线观看日韩| 美女国产视频在线观看| 国产在视频线精品| 日韩精品免费视频一区二区三区| 伊人亚洲综合成人网| 久久久久国产网址| 麻豆av在线久日| 欧美国产精品一级二级三级| 国产精品一二三区在线看| 久久久久久久久久久免费av| 美女xxoo啪啪120秒动态图| 欧美老熟妇乱子伦牲交| 观看美女的网站| 91精品三级在线观看| 天美传媒精品一区二区| 日韩,欧美,国产一区二区三区| 国产综合精华液| 99re6热这里在线精品视频| 综合色丁香网| 五月开心婷婷网| 久久韩国三级中文字幕| 欧美日韩综合久久久久久| 男男h啪啪无遮挡| 欧美最新免费一区二区三区| 一本大道久久a久久精品| 免费观看在线日韩| 国产 精品1| 一级,二级,三级黄色视频| 亚洲欧美精品综合一区二区三区 | 午夜日韩欧美国产| 欧美日韩综合久久久久久| 国语对白做爰xxxⅹ性视频网站| 黑丝袜美女国产一区| 精品国产国语对白av| 90打野战视频偷拍视频| 亚洲精品国产av成人精品| 午夜福利在线观看免费完整高清在| 不卡av一区二区三区| 午夜福利网站1000一区二区三区| 亚洲欧美日韩另类电影网站| 国产精品偷伦视频观看了| 日韩中文字幕欧美一区二区 | 国产福利在线免费观看视频| 一二三四在线观看免费中文在| 亚洲一区二区三区欧美精品| 亚洲av综合色区一区| 欧美日韩精品成人综合77777| 飞空精品影院首页| 日韩一卡2卡3卡4卡2021年| 国产野战对白在线观看| 午夜激情久久久久久久| 日韩在线高清观看一区二区三区| 久久青草综合色| 女性生殖器流出的白浆| 高清欧美精品videossex| 久久人妻熟女aⅴ| 欧美国产精品va在线观看不卡| 纵有疾风起免费观看全集完整版| 一级a爱视频在线免费观看| 久久亚洲国产成人精品v| 我的亚洲天堂| 波野结衣二区三区在线| 国产精品熟女久久久久浪| 伊人亚洲综合成人网| 亚洲五月色婷婷综合| 18在线观看网站| 超碰97精品在线观看| av福利片在线| 在线观看人妻少妇| 一级片'在线观看视频| 国产午夜精品一二区理论片| 最近中文字幕高清免费大全6| 成人免费观看视频高清| 国语对白做爰xxxⅹ性视频网站| 亚洲欧美一区二区三区黑人 | 26uuu在线亚洲综合色| 97在线人人人人妻| 国产精品免费大片| 亚洲欧美成人精品一区二区| 黄色怎么调成土黄色| 岛国毛片在线播放| 日韩欧美精品免费久久| 中文字幕人妻丝袜制服| 最近的中文字幕免费完整| 欧美日韩成人在线一区二区| 熟妇人妻不卡中文字幕| 爱豆传媒免费全集在线观看| 午夜免费男女啪啪视频观看| 啦啦啦在线免费观看视频4| 欧美最新免费一区二区三区| 男人爽女人下面视频在线观看| 久久久亚洲精品成人影院| 热re99久久精品国产66热6| 免费在线观看黄色视频的| 丁香六月天网| 欧美中文综合在线视频| 色播在线永久视频| 少妇被粗大的猛进出69影院| 精品久久蜜臀av无| 伊人久久大香线蕉亚洲五| 一区二区三区精品91| a 毛片基地| 国产免费福利视频在线观看| 99热全是精品| 免费在线观看黄色视频的| 最近手机中文字幕大全| 成人亚洲欧美一区二区av| 中文乱码字字幕精品一区二区三区| 亚洲激情五月婷婷啪啪| 三上悠亚av全集在线观看| 日韩大片免费观看网站| 国产成人午夜福利电影在线观看| 啦啦啦在线免费观看视频4| 香蕉丝袜av| 下体分泌物呈黄色| 青草久久国产| av在线app专区| 亚洲欧洲日产国产| 成年人免费黄色播放视频| 国语对白做爰xxxⅹ性视频网站| 美女大奶头黄色视频| 亚洲精品在线美女| 亚洲av男天堂| av视频免费观看在线观看| 久久精品久久精品一区二区三区| 国产男女超爽视频在线观看| 美女国产高潮福利片在线看| 国产精品一国产av| av线在线观看网站| 欧美国产精品一级二级三级| 99re6热这里在线精品视频| 欧美日韩国产mv在线观看视频| 亚洲av欧美aⅴ国产| 婷婷成人精品国产| 少妇人妻久久综合中文| 一区二区日韩欧美中文字幕| 黄色怎么调成土黄色| 国产有黄有色有爽视频| 日韩一区二区三区影片| 精品久久久久久电影网| 赤兔流量卡办理| 黄色毛片三级朝国网站| 国产成人精品无人区| 超碰97精品在线观看| 七月丁香在线播放| 国产又色又爽无遮挡免| 色播在线永久视频| 日韩av在线免费看完整版不卡| 国产精品久久久久久久久免| 成人国产av品久久久| 亚洲中文av在线| 春色校园在线视频观看| 18在线观看网站| 欧美成人午夜免费资源| 97在线视频观看| 亚洲中文av在线| 国产精品嫩草影院av在线观看| 国产女主播在线喷水免费视频网站| 天堂俺去俺来也www色官网| 日韩制服骚丝袜av| 久久久欧美国产精品| 99久久综合免费| 久久免费观看电影| 韩国精品一区二区三区| 午夜久久久在线观看| 欧美日本中文国产一区发布| 久久国内精品自在自线图片| 亚洲色图 男人天堂 中文字幕| 波多野结衣av一区二区av| 国产亚洲最大av| 亚洲av电影在线进入| 在线 av 中文字幕| 久久久久久久久久久久大奶| 一二三四在线观看免费中文在| 日产精品乱码卡一卡2卡三| 国产乱来视频区| 大香蕉久久成人网| 久久久久久久国产电影| 日韩制服骚丝袜av| 国产亚洲av片在线观看秒播厂| 777米奇影视久久| 美女中出高潮动态图| 国产片内射在线| 日韩熟女老妇一区二区性免费视频| 亚洲国产精品成人久久小说| av在线观看视频网站免费| 婷婷成人精品国产| 又大又黄又爽视频免费| 日日啪夜夜爽| 亚洲av福利一区| 久久精品aⅴ一区二区三区四区 | 飞空精品影院首页| 精品久久久久久电影网| 国产精品久久久av美女十八| 国产精品免费大片| 80岁老熟妇乱子伦牲交| 国产精品一区二区在线不卡| 高清视频免费观看一区二区| 精品少妇久久久久久888优播| 免费大片黄手机在线观看| 丝袜美足系列| 免费女性裸体啪啪无遮挡网站| 热99国产精品久久久久久7| 交换朋友夫妻互换小说| 亚洲欧美精品综合一区二区三区 | 91在线精品国自产拍蜜月| 国产免费现黄频在线看| 亚洲精品国产av成人精品| 男女下面插进去视频免费观看| 国产精品免费视频内射| 两性夫妻黄色片| 久久久久久久大尺度免费视频| 久久久久精品久久久久真实原创| 国产野战对白在线观看| 欧美日韩亚洲国产一区二区在线观看 | 欧美日韩一级在线毛片| 国产色婷婷99| 观看美女的网站| 永久免费av网站大全| 国产欧美日韩一区二区三区在线| 国产精品欧美亚洲77777| 欧美最新免费一区二区三区| 日韩中文字幕欧美一区二区 | 欧美日本中文国产一区发布| 曰老女人黄片| 视频区图区小说| 赤兔流量卡办理| 亚洲精品国产av成人精品| 欧美人与性动交α欧美精品济南到 | 亚洲一码二码三码区别大吗| 久久久久久人妻| 制服人妻中文乱码| 国产成人精品久久久久久| 在线亚洲精品国产二区图片欧美| 激情五月婷婷亚洲| 精品国产露脸久久av麻豆| 啦啦啦啦在线视频资源| 熟女少妇亚洲综合色aaa.| 女人久久www免费人成看片| 日韩一区二区三区影片| 亚洲国产欧美在线一区| 久久精品人人爽人人爽视色| 韩国高清视频一区二区三区| 大片电影免费在线观看免费| 人人妻人人澡人人爽人人夜夜| 岛国毛片在线播放| 欧美精品亚洲一区二区| 99久久中文字幕三级久久日本| 欧美+日韩+精品| 水蜜桃什么品种好| 最近最新中文字幕免费大全7| 久久久久国产网址| 在线免费观看不下载黄p国产| 成年女人毛片免费观看观看9 | 亚洲,一卡二卡三卡| 中文字幕另类日韩欧美亚洲嫩草| 黄色 视频免费看| www.熟女人妻精品国产| 欧美日韩一区二区视频在线观看视频在线| 国产精品久久久久久久久免|