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

    冰雹動(dòng)態(tài)本構(gòu)建模與驗(yàn)證

    2023-11-01 10:17:33王計(jì)真
    航空科學(xué)技術(shù) 2023年8期

    王計(jì)真

    摘 要:飛機(jī)飛行過程中,機(jī)身外表容易遭受冰雹撞擊,造成機(jī)體破壞,危害乘客安全。為研究結(jié)構(gòu)抗冰雹沖擊性能,本文基于拉格朗日(Lagrange)方法、耦合歐拉-拉格朗日方法(CEL)和光滑粒子動(dòng)力學(xué)(SPH)方法,建立冰雹動(dòng)態(tài)本構(gòu)模型,并根據(jù)冰撞剛性靶試驗(yàn)數(shù)據(jù),修正冰雹模型參數(shù),驗(yàn)證三種方法的有效性。數(shù)值計(jì)算結(jié)果表明,三種方法均能較好地描述冰雹動(dòng)態(tài)破碎行為,且在中低速冰撞分析時(shí),基于應(yīng)變率強(qiáng)化效應(yīng)和張力失效準(zhǔn)則的Lagrange方法與試驗(yàn)有更好的一致性。

    關(guān)鍵詞:冰雹撞擊; 動(dòng)態(tài)本構(gòu); 冰撞剛性靶; 應(yīng)變率強(qiáng)化; 張力失效

    中圖分類號(hào):TG146.23 文獻(xiàn)標(biāo)識(shí)碼:A DOI:10.19452/j.issn1007-5453.2023.08.007

    基金項(xiàng)目: 航空科學(xué)基金(2016ZA23005)

    民機(jī)在飛行或起降過程中,容易遭受冰雹的沖擊威脅。冰雹沖擊過程屬于小質(zhì)量高速?zèng)_擊,對(duì)飛機(jī)結(jié)構(gòu)尤其是復(fù)合材料結(jié)構(gòu)安全產(chǎn)生重大影響,嚴(yán)重威脅乘員的生命安全,帶來嚴(yán)重的經(jīng)濟(jì)損失。飛機(jī)結(jié)構(gòu)抗冰雹撞擊研究具有重要理論與工程意義。

    冰有20余種結(jié)構(gòu)形式(晶體結(jié)構(gòu)或非晶體狀態(tài)),因此,構(gòu)建一種普適材料模型描述其動(dòng)態(tài)行為極為困難,且冰雹撞擊分析僅在航空、軌道交通和風(fēng)電等少數(shù)行業(yè)存在需求,因而針對(duì)冰雹動(dòng)態(tài)本構(gòu)研究有限。20世紀(jì)70年代,Haynes[1]首次實(shí)現(xiàn)了冰的靜態(tài)及準(zhǔn)靜態(tài)壓縮性能測(cè)試。隨后,Schulson等[2-3]和Dempsey等[4]研究了不同溫度和晶體結(jié)構(gòu)下冰的靜態(tài)力學(xué)性能和裂紋擴(kuò)展屬性。研究結(jié)果表明,冰雹的拉伸強(qiáng)度遠(yuǎn)小于壓縮強(qiáng)度,且靜態(tài)壓縮過程表現(xiàn)出韌性破壞特性。近年來,冰雹動(dòng)態(tài)本構(gòu)的研究逐步開展,表明冰為典型的應(yīng)變率敏感材料。Botto[5]開展了低應(yīng)變率下(10?8 ~10?3s?1)冰結(jié)構(gòu)的動(dòng)態(tài)力學(xué)性能測(cè)試,觀察到隨應(yīng)變率增加材料特性由韌性向脆性轉(zhuǎn)變。Kim等[6]和Shazly等[7]采用霍普金森桿測(cè)試高應(yīng)變率動(dòng)態(tài)壓縮力學(xué)性能,結(jié)果表明冰的壓縮強(qiáng)度隨應(yīng)變率的提高而增大。宋振華等[8-9]修正前人試驗(yàn)中的不合理因素,獲取更高應(yīng)變率下壓縮強(qiáng)度。相關(guān)文獻(xiàn)中很少有關(guān)于冰雹動(dòng)態(tài)本構(gòu)建模的研究,彈塑性本構(gòu)和彈塑性水動(dòng)力本構(gòu)[10]都曾被用于冰雹撞擊分析,但由于未考慮高速?zèng)_擊過程的應(yīng)變率效應(yīng),模擬精度不高。Tippmann等[11]也采用彈塑性本構(gòu),但考慮了屈服強(qiáng)度的應(yīng)變率強(qiáng)化效應(yīng)和張力失效特性,建立了一種冰雹模型,該模型得到的撞擊載荷峰值較為準(zhǔn)確,但峰值出現(xiàn)時(shí)間不夠準(zhǔn)確,且網(wǎng)格依賴性太高,求解效率有待提高。最近Sain等[12]較為深入地研究了冰雹撞擊破壞機(jī)理,并在總結(jié)前人研究成果的基礎(chǔ)上,建立了一種唯象本構(gòu)模型,該模型將彈塑性模型和狀態(tài)方程進(jìn)行結(jié)合,考慮流變應(yīng)力的效應(yīng)以更真實(shí)模擬材料失效后的承載特性。張笑宇等[13]針對(duì)碳纖維復(fù)合材料蜂窩夾芯板,進(jìn)行了冰雹沖擊及多次沖擊數(shù)值研究,針對(duì)多次沖擊工況,分析了不同沖擊能量、沖擊角度以及沖擊位置的多次沖擊損傷,以及在冰雹多次沖擊下的損傷累積和疊加。

    近十年來,美國聯(lián)邦航空條例(FAR)、歐洲聯(lián)合航空要求(JAR)和中國民用航空規(guī)章(CCAR)都對(duì)近代飛機(jī)結(jié)構(gòu)提出專門要求,以保證結(jié)構(gòu)受到大量冰雹撞擊后仍能保持其安全飛行和著陸的功能。1984年美國聯(lián)邦航空局(FAA)在其發(fā)布的具有指導(dǎo)作用的咨詢通報(bào)《復(fù)合材料飛機(jī)結(jié)構(gòu)》(AC 20-107A)中,對(duì)包括冰雹在內(nèi)的離散源撞擊飛機(jī)結(jié)構(gòu)的沖擊安全問題有專門規(guī)定,所有即將投入商業(yè)運(yùn)行的飛機(jī)的設(shè)計(jì)必須嚴(yán)格滿足相應(yīng)的離散源撞擊指標(biāo)。中國民用航空規(guī)章第25部《運(yùn)輸類飛機(jī)適航標(biāo)準(zhǔn)》(CCAR-25-R4)也規(guī)定了結(jié)構(gòu)元件不能因包括冰雹在內(nèi)的離散源損傷而失效。

    本文在冰雹撞擊剛性靶試驗(yàn)的基礎(chǔ)上,基于有限元軟件ABAQUS,采用拉格朗日(Lagrange)、CEL和SPH三種方法建立了冰雹本構(gòu)模型,并對(duì)比研究了三種模型的適用性。

    1 冰雹撞擊剛性靶試驗(yàn)

    1.1 冰彈制備

    用注射器將去離子純凈水注入模具;將模具放置在冰箱中,調(diào)節(jié)冰箱為-18℃±3℃溫度環(huán)境,保持溫差±2℃以內(nèi)冰凍6h以上;取出冰彈及模具,浸入25℃±5℃的水中直至冰彈可在模具中滑動(dòng),將冰彈取出放置于-10℃±3℃的環(huán)境中待用。本次試驗(yàn)根據(jù)ASTM F320[14]要求選擇最大尺寸(直徑50mm)的球形冰彈,質(zhì)量約為60g。

    1.2 試驗(yàn)方法

    剛性靶板冰撞試驗(yàn)系統(tǒng)如圖1所示,包括D80氣炮系統(tǒng)、控制系統(tǒng)、激光瞄準(zhǔn)系統(tǒng)、數(shù)據(jù)采集系統(tǒng)、力傳感器、剛性靶板和高速攝像機(jī)測(cè)速系統(tǒng)等。試驗(yàn)過程中通過力傳感器測(cè)量冰撞載荷數(shù)據(jù),高速攝像機(jī)對(duì)冰撞過程進(jìn)行全程拍攝。力傳感器安裝如圖2所示,單個(gè)力傳感器垂向量程為150kN,根據(jù)載荷大小選擇單載荷傳感器或雙載荷傳感器并聯(lián)。

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

    自然環(huán)境中的冰雹凝結(jié)溫度在一定范圍內(nèi)變化,為研究?jī)鼋Y(jié)溫度對(duì)冰雹力學(xué)性能的影響,選擇溫度范圍-5~-30℃制作冰彈,以200m/s速度開展剛性靶板冰撞試驗(yàn),獲取載荷峰值與凍結(jié)溫度的關(guān)系,見表1??梢钥闯?,冰撞載荷峰值變化小于1.25%,表明冰彈力學(xué)性能對(duì)凍結(jié)溫度不敏感,后續(xù)冰彈制作溫度均選為-10℃。

    為研究冰雹撞擊過程,在50~300m/s速度范圍內(nèi)開展純凈冰彈剛性靶板沖擊試驗(yàn)。圖3是沖擊速度為50m/s時(shí)冰彈沖擊剛性靶板過程。由圖3可以看出,在t=0時(shí),冰彈開始與靶板垂直接觸,在0~0.3ms時(shí),冰彈開始碰撞剛性靶板至完全破碎,冰彈碎屑開始沿靶板平面橫向飛散。同時(shí)通過觀測(cè)高速攝像可以發(fā)現(xiàn),冰雹沖擊破壞形式主要以脆性斷裂為主,斷裂后,冰雹開始沿靶板平面橫向擴(kuò)散,呈現(xiàn)明顯的塑性流動(dòng)特性。

    圖4給出了50m/s、150m/s和300m/s三組不同速度下,冰彈的沖擊載荷曲線,圖5給出了冰彈沖擊載荷峰值及達(dá)到載荷峰值時(shí)刻隨沖擊速度的變化。由圖4、圖5可以看出,隨著沖擊速度的增加,沖擊載荷峰值逐漸增大,擬合載荷峰值與沖擊速度關(guān)系曲線,得到?jīng)_擊載荷峰值與沖擊速度的平方近似呈線性關(guān)系;隨著沖擊速度的提高,達(dá)到載荷峰值的時(shí)刻提前。

    2 冰雹動(dòng)態(tài)本構(gòu)模型

    2.1 材料模型

    冰從力學(xué)屬性上可分為三類:?jiǎn)尉П?、多晶冰和弱冰(如雪)。本文中冰從制作過程及其微觀檢測(cè)結(jié)果均表明,在結(jié)構(gòu)上和性能上屬于多晶冰范疇。

    結(jié)合前人分析,總結(jié)冰雹的材料屬性包括:冰雹與應(yīng)變率有很強(qiáng)的依賴關(guān)系,當(dāng)從低應(yīng)變率增大到高應(yīng)變率時(shí),冰的力學(xué)行為從韌性轉(zhuǎn)變?yōu)榇嘈?;冰雹的壓縮強(qiáng)度隨應(yīng)變率的增大而增大,拉伸強(qiáng)度比壓縮強(qiáng)度要小得多;冰雹高速撞擊時(shí),在相當(dāng)大的變形下完全破裂,之后表現(xiàn)出流體特性。結(jié)合實(shí)驗(yàn)室制備的冰球模型和前人的研究,給出冰球基本力學(xué)參數(shù),見表2。

    2.2 Lagrange方法

    基于Lagrange方法的有限元模型,采用彈塑性本構(gòu)關(guān)系、應(yīng)變率強(qiáng)化效應(yīng)和張力失效準(zhǔn)則定義冰的本構(gòu)模型。在冰雹的高速撞擊條件下,碰撞物體會(huì)發(fā)生硬化現(xiàn)象,高應(yīng)變率的硬化作用改變了原定的冰雹力學(xué)過程,屈服強(qiáng)度發(fā)生動(dòng)態(tài)上升,因?yàn)榍?qiáng)度并不唯一,所以本研究對(duì)塑性階段屈服強(qiáng)度值的設(shè)定,采用了基于應(yīng)變率的強(qiáng)度設(shè)置,本研究的應(yīng)變率強(qiáng)化效應(yīng)見表3。

    ABAQUS 的軟件動(dòng)態(tài)失效模型,其中常見的有考慮應(yīng)變率效應(yīng)的屈服準(zhǔn)則、Johnson-Cook失效準(zhǔn)則、剪切失效準(zhǔn)則、張力失效準(zhǔn)則等,其中剪切失效準(zhǔn)則與張力失效準(zhǔn)則均適用于高應(yīng)變率的動(dòng)態(tài)問題。張力失效準(zhǔn)則在宏觀上描述了物質(zhì)因?yàn)閺埩虞d作用而被破壞失效的過程,在微觀層面上,它用靜水壓強(qiáng)(壓力)應(yīng)力來衡量模型的動(dòng)態(tài)碎裂和切斷力,同時(shí)該準(zhǔn)則也可以作為中間準(zhǔn)則與其他準(zhǔn)則相連接,共同定義破壞失效?;趶埩κУ膭?dòng)態(tài)失效模型,可較好地模擬冰雹沖擊破壞過程。張力失效準(zhǔn)則受預(yù)設(shè)場(chǎng)變量值的作用,當(dāng)材料單元節(jié)點(diǎn)靜水壓力應(yīng)力達(dá)到準(zhǔn)則設(shè)定值時(shí),該點(diǎn)就發(fā)生失效。

    冰雹材料進(jìn)行的設(shè)置是將偏應(yīng)力失效類型設(shè)置為脆性,壓應(yīng)力失效類型設(shè)置為韌性。在這種設(shè)置下,一旦壓應(yīng)力達(dá)到失效臨界值,偏應(yīng)力將會(huì)被設(shè)置為零,并長(zhǎng)期保持為零,根據(jù)準(zhǔn)則此時(shí)在單元上的應(yīng)力受到切斷應(yīng)力臨界值的限制,當(dāng)靜水張力或靜水壓力超過了切斷應(yīng)力,應(yīng)力就只表現(xiàn)靜水壓力和靜水張力,失效單元表現(xiàn)出類似于流體的行為。當(dāng)應(yīng)力狀態(tài)將要超過截?cái)鄩毫?,?yīng)力狀態(tài)進(jìn)入后失效狀態(tài)時(shí),切應(yīng)力減小為零,單元將不再承受切應(yīng)力只承受靜水壓應(yīng)力。通過與試驗(yàn)對(duì)比,此模型可較好地模擬冰雹的動(dòng)態(tài)碎裂過程。

    2.3 SPH方法和CEL方法

    該模型可正確地代表冰雹在撞擊早期階段的冰雹行為,開始時(shí)是用大剛度來表征其特點(diǎn),而在隨后階段,即撞擊后冰雹開裂,其行為特性表現(xiàn)為流體特性。材料模型是用拉伸張力失效判據(jù)來表征的,當(dāng)達(dá)到拉伸失效應(yīng)力時(shí),偏應(yīng)力分量置零,并且材料只能保持壓縮應(yīng)力。在計(jì)算中,只采用彈塑性水動(dòng)力材料模型的失效判據(jù),而不用塑性失效應(yīng)變判據(jù),設(shè)置張力失效準(zhǔn)則,同Lagrange方法一致。

    3 冰雹撞擊傳感器的仿真分析

    3.1 有限元模型

    圖6所示為冰雹沖擊傳感器的數(shù)值仿真模型和試驗(yàn)?zāi)P汀楹?jiǎn)化分析,認(rèn)為傳感器為剛體,且底部固定。為提高計(jì)算效率,在有限元仿真時(shí)采用1/4模型。

    3.2 仿真結(jié)果

    圖7為通過Lagrange方法、CEL方法和SPH方法進(jìn)行數(shù)值仿真給出的冰雹撞擊變形過程。從變形來看,三者冰雹變形形式較為相似,均能體現(xiàn)出冰雹碎裂之后的塑性流動(dòng)過程,和試驗(yàn)現(xiàn)象較為一致。

    圖8為當(dāng)沖擊速度為49.75m/s時(shí),采用不同方法的直徑55mm冰彈撞擊力-時(shí)間歷程的比較,仿真中的撞擊力采用的冰雹與傳感器的法向撞擊力。從圖8中可以看出,三種分析方法均能模擬出冰雹的材料硬化性質(zhì)和脆性斷裂后的塑性流動(dòng)過程,在中低速時(shí)SPH方法收斂性較差,不能較好地模擬出載荷歷程,而CEL方法的載荷峰值與試驗(yàn)結(jié)果誤差偏小,Lagrange方法的峰值較為接近。但仿真中的峰值時(shí)間較試驗(yàn)值提前,主要原因是本文未考慮測(cè)力傳感器的剛度,仿真過程中對(duì)接觸過程的處理也較為理想化,還需要進(jìn)一步完善。

    4 結(jié)論

    本研究根據(jù)冰雹粒子的物理特性,分析冰雹沖擊動(dòng)態(tài)本構(gòu),研究復(fù)合材料動(dòng)態(tài)本構(gòu)和損傷模型,并基于ABAQUS軟件對(duì)冰雹撞擊進(jìn)行數(shù)值模擬??傻玫饺缦陆Y(jié)論:

    (1) 開展50~300m/s的沖擊試驗(yàn),從試驗(yàn)數(shù)據(jù)得出沖擊載荷峰值與沖擊速度的平方近似呈線性關(guān)系,隨著沖擊速度的提高載荷峰值時(shí)刻提前。

    (2) 冰雹沖擊破壞形式主要以脆性斷裂為主,脆性斷裂后,呈現(xiàn)明顯的塑性流動(dòng)特性,破壞前無明顯的屈服階段,沖擊壓縮載荷作用下的冰雹呈現(xiàn)明顯應(yīng)變率強(qiáng)化效應(yīng)。

    (3) CEL方法、SPH方法和Lagrange方法均能描述冰雹動(dòng)態(tài)破壞行為,中低速時(shí)基于應(yīng)變率強(qiáng)化效應(yīng)和張力失效準(zhǔn)則的Largrange方法和試驗(yàn)一致性更好。

    參考文獻(xiàn)

    [1]Haynes F D. Effect of temperature on the strength of snow-ice[D]. Hanover: Cold Regions Research and Engineering Lab, 1978.

    [2]Currier J H, Schulson E M. The tensile strength of ice as a func‐tion of grain size[J]. Acta Metallurgica, 1982, 30(8): 1511-1514.

    [3]Nixon W A, Schulson E M. A micromechanical view of the fracture toughness of ice[J]. Le Journal de Physique Colloques, 1987, 48(C1): 313-319.

    [4]Dempsey J P, DeFranco S J, Adamson R M, et al. Scale effects on the in-situ tensile strength and fracture of ice. Part I: Large grained freshwater ice at Spray Lakes Reservoir, Alberta[M]. Fracture Scaling: Springer Netherlands, 1999.

    [5]Batto R A, Schulson E M. On the ductile-to-brittle transition in ice under compression[J]. Acta Metallurgica et Materialia, 1993, 41(7): 2219-2225.

    [6]Kim H, Keune J N. Compressive strength of ice at impact strain rates[J]. Journal of Materials Science, 2007, 42(8): 2802-2806.

    [7]Shazly M, Prakash V, Lerch B A. High strain-rate behavior of ice under uniaxial compression[J]. International Journal of Solids and Structures, 2009, 46(6): 1499-1515.

    [8]宋振華. 冰載荷作用下碳纖維復(fù)合材料桁條加筋曲面板的沖擊動(dòng)力響應(yīng)研究[D]. 廣州:暨南大學(xué), 2014. Song Zhenhua. The dynamic response of stringer stiffened curved composite panels under the hail ice impact[D]. Guang‐zhou: Ji’nan University,2014.(in Chinese)

    [9]Song Z H, Wang Z, Kim H, et al. Pulse shaper design and dynamic mechanical property tests on ice using the split hopkinson pressure bar[J]. International Journal of Impact Engineering, 2016, 12(2):1-5.

    [10]Anghileri M, Castelletti L M L, Invernizzi F, et al. A survey of numerical models for hail impact analysis using explicit finite element codes[J]. International Journal of Impact Engineering, 2005, 31(8): 929-944.

    [11]Tippmann J D, Kim H, Rhymer J D. Experimentally validated strain rate dependent material model for spherical ice impact simulation[J]. International Journal of Impact Engineering, 2013, 57:43-54.

    [12]Sain T, Narasimhan R. Constitutive modeling of ice in the high strain rate regime[J]. International Journal of Solids and Structures, 2011, 48(5): 817-827.

    [13]張笑宇,徐緋,張玉林,等.復(fù)合材料蜂窩夾芯板冰雹沖擊及多次沖擊數(shù)值研究[J].航空科學(xué)技術(shù),2021,32(12):74-83. Zhang Xiaoyu, Xu Fei, Zhang Yulin, et al. Numerical study on single and repeated impact of hail on composite honeycomb sandwich panels[J].Aeronautical Science & Technology, 2021, 32(12):74-83.(in Chinese)

    [14]ASTM International. Standard test method for hail impact resistance of aerospace transparent enclosures[R]. Report No. F320-10, 2010.

    Modeling and Verification of the Dynamic Constitutive of the Hailstone

    Wang Jizhen

    National Key Laboratory of Strength and Structural Integrity, Aircraft Strength Research Institute of China, Xi’an 710065, China

    Abstract: During the flight, the fuselage may be hit by the hailstone, causing structures damage and endangering passenger safety. In order to study the structural hail impact resistance, a hailstone dynamic constitutive model based on Lagrange method, Coupled Euler-Lagrange (CEL) method and Smooth Particle Dynamics(SPH ) method were established. Then the model was verificated by comparing with test data of ice impact rigid panel. The numerical results show that the three methods can describe the dynamic breaking behavior of hailstone. During the low and medium velocity impact, the Lagrange method, based on strain rate enhancement and tension failure criterion, has better consistency with test data.

    Key Words: hailstone impact; dynamic constitutive; ice impact rigid panel; strain rate enhancement; tension failure

    少妇的丰满在线观看| 欧美日韩综合久久久久久| 97精品久久久久久久久久精品| 欧美国产精品一级二级三级| 亚洲成国产人片在线观看| 久久这里只有精品19| 丝袜人妻中文字幕| 人人妻人人澡人人爽人人夜夜| 国产高清三级在线| 水蜜桃什么品种好| 91精品国产国语对白视频| 大片电影免费在线观看免费| 午夜日本视频在线| 伊人久久国产一区二区| 曰老女人黄片| 亚洲,一卡二卡三卡| 国产一区二区在线观看av| 插逼视频在线观看| 国产av一区二区精品久久| 午夜精品国产一区二区电影| 日韩免费高清中文字幕av| 9色porny在线观看| 有码 亚洲区| 欧美日韩精品成人综合77777| 免费女性裸体啪啪无遮挡网站| 欧美变态另类bdsm刘玥| 高清毛片免费看| 国产精品一区www在线观看| 日本爱情动作片www.在线观看| 日韩精品有码人妻一区| 午夜av观看不卡| 你懂的网址亚洲精品在线观看| 亚洲性久久影院| 高清在线视频一区二区三区| 97人妻天天添夜夜摸| 精品国产国语对白av| 国产激情久久老熟女| 欧美bdsm另类| 久久99蜜桃精品久久| 亚洲精品自拍成人| 中文字幕人妻丝袜制服| 欧美激情极品国产一区二区三区 | 午夜影院在线不卡| 一二三四在线观看免费中文在 | 纯流量卡能插随身wifi吗| 熟女人妻精品中文字幕| 日韩免费高清中文字幕av| 免费观看无遮挡的男女| 丰满迷人的少妇在线观看| 最近最新中文字幕免费大全7| 免费在线观看黄色视频的| 欧美精品人与动牲交sv欧美| 一级毛片黄色毛片免费观看视频| 女性生殖器流出的白浆| 国产成人一区二区在线| 久久这里有精品视频免费| tube8黄色片| 高清毛片免费看| 国产乱人偷精品视频| 三上悠亚av全集在线观看| 久热久热在线精品观看| 九色亚洲精品在线播放| 欧美国产精品va在线观看不卡| 嫩草影院入口| xxxhd国产人妻xxx| 天堂8中文在线网| 欧美另类一区| 国国产精品蜜臀av免费| 成人国语在线视频| 最新中文字幕久久久久| 欧美亚洲 丝袜 人妻 在线| 午夜福利网站1000一区二区三区| 国产高清不卡午夜福利| 精品熟女少妇av免费看| 99国产综合亚洲精品| 国产白丝娇喘喷水9色精品| 免费大片黄手机在线观看| 久久精品国产亚洲av涩爱| 国产午夜精品一二区理论片| 久久久欧美国产精品| 热99久久久久精品小说推荐| 午夜免费观看性视频| 精品国产一区二区三区四区第35| 国产色爽女视频免费观看| 国产色爽女视频免费观看| 美女中出高潮动态图| 美女中出高潮动态图| 久久精品aⅴ一区二区三区四区 | 丰满乱子伦码专区| 亚洲综合色网址| 又粗又硬又长又爽又黄的视频| 日韩不卡一区二区三区视频在线| 中国三级夫妇交换| 中文字幕最新亚洲高清| 美女xxoo啪啪120秒动态图| 有码 亚洲区| 亚洲av电影在线观看一区二区三区| 狂野欧美激情性xxxx在线观看| 中文精品一卡2卡3卡4更新| 亚洲国产av影院在线观看| 又大又黄又爽视频免费| 日本与韩国留学比较| 亚洲欧美中文字幕日韩二区| 亚洲综合色惰| 久久这里只有精品19| 女人精品久久久久毛片| 日本vs欧美在线观看视频| 日韩伦理黄色片| 草草在线视频免费看| 成人国语在线视频| 一区二区av电影网| 欧美bdsm另类| 国产欧美日韩综合在线一区二区| av免费观看日本| 国产熟女午夜一区二区三区| 久久精品国产自在天天线| 国产欧美日韩综合在线一区二区| 国产欧美日韩综合在线一区二区| 黄片无遮挡物在线观看| 一级,二级,三级黄色视频| 亚洲av国产av综合av卡| 欧美亚洲 丝袜 人妻 在线| 一区二区三区精品91| freevideosex欧美| 免费人妻精品一区二区三区视频| 最近最新中文字幕大全免费视频 | 成人二区视频| 亚洲欧美清纯卡通| 欧美成人午夜精品| 亚洲,欧美精品.| 少妇 在线观看| 亚洲第一区二区三区不卡| 中文字幕最新亚洲高清| 热re99久久国产66热| 亚洲av综合色区一区| 国产午夜精品一二区理论片| a级毛片在线看网站| 大片电影免费在线观看免费| 又黄又粗又硬又大视频| videos熟女内射| 天堂中文最新版在线下载| 中文精品一卡2卡3卡4更新| 黑人高潮一二区| 精品国产一区二区三区久久久樱花| 久久久亚洲精品成人影院| 激情五月婷婷亚洲| 久久99热6这里只有精品| 亚洲,欧美,日韩| 如日韩欧美国产精品一区二区三区| 在线观看美女被高潮喷水网站| 香蕉精品网在线| 女人久久www免费人成看片| 又大又黄又爽视频免费| 精品国产国语对白av| 日日啪夜夜爽| 香蕉丝袜av| 久久人人爽人人爽人人片va| 欧美xxxx性猛交bbbb| 欧美激情 高清一区二区三区| 制服诱惑二区| 国产av国产精品国产| 在线天堂最新版资源| 国产精品国产三级国产av玫瑰| 热99国产精品久久久久久7| 一区二区日韩欧美中文字幕 | 亚洲欧美清纯卡通| 天天影视国产精品| 美女大奶头黄色视频| 91久久精品国产一区二区三区| 国产精品一二三区在线看| 亚洲精品美女久久av网站| 日本wwww免费看| 人妻人人澡人人爽人人| 久久狼人影院| 2022亚洲国产成人精品| 亚洲精品久久午夜乱码| 乱人伦中国视频| 日本黄大片高清| 国产视频首页在线观看| 亚洲精品456在线播放app| 90打野战视频偷拍视频| 少妇的丰满在线观看| 最近2019中文字幕mv第一页| 欧美+日韩+精品| 久久这里只有精品19| 99久久精品国产国产毛片| 亚洲精品久久成人aⅴ小说| 日本与韩国留学比较| 91精品国产国语对白视频| 精品亚洲成国产av| 天天操日日干夜夜撸| 制服人妻中文乱码| 欧美日韩亚洲高清精品| 久久久久久久久久人人人人人人| 久久久国产一区二区| 久久久久国产精品人妻一区二区| 欧美激情国产日韩精品一区| 桃花免费在线播放| 国产成人aa在线观看| 亚洲国产日韩一区二区| 一级毛片 在线播放| 丝袜脚勾引网站| 99re6热这里在线精品视频| 精品久久久久久电影网| 亚洲精品日本国产第一区| 欧美老熟妇乱子伦牲交| 精品国产乱码久久久久久小说| 高清欧美精品videossex| 亚洲丝袜综合中文字幕| 一区在线观看完整版| 亚洲精品国产av成人精品| 久久毛片免费看一区二区三区| 观看av在线不卡| 亚洲精品国产av成人精品| 国产不卡av网站在线观看| 免费播放大片免费观看视频在线观看| 女性被躁到高潮视频| 国产成人一区二区在线| 99热网站在线观看| 少妇人妻精品综合一区二区| 大香蕉久久网| 日韩av免费高清视频| 成人二区视频| 亚洲国产精品一区二区三区在线| 9色porny在线观看| 97在线人人人人妻| 视频在线观看一区二区三区| 啦啦啦视频在线资源免费观看| 国产免费又黄又爽又色| av播播在线观看一区| 高清欧美精品videossex| 久久av网站| av免费在线看不卡| 麻豆乱淫一区二区| 精品人妻熟女毛片av久久网站| 亚洲第一区二区三区不卡| 高清黄色对白视频在线免费看| 国产男女超爽视频在线观看| 夫妻性生交免费视频一级片| 黄色配什么色好看| 国产成人精品久久久久久| 美女主播在线视频| 亚洲精品国产色婷婷电影| 99九九在线精品视频| 国产69精品久久久久777片| 国产一区有黄有色的免费视频| 如何舔出高潮| 在线观看国产h片| 久久韩国三级中文字幕| 2021少妇久久久久久久久久久| 热re99久久精品国产66热6| 日本91视频免费播放| 欧美人与善性xxx| 久久久久久久久久久久大奶| 亚洲欧洲国产日韩| 一级a做视频免费观看| 亚洲成人手机| 久久99热6这里只有精品| 日韩一区二区视频免费看| 久久久久久久久久久免费av| 最新的欧美精品一区二区| 久久久久视频综合| 国产成人a∨麻豆精品| 亚洲五月色婷婷综合| 男女午夜视频在线观看 | 国产麻豆69| xxxhd国产人妻xxx| 97人妻天天添夜夜摸| 九九爱精品视频在线观看| 热99久久久久精品小说推荐| 欧美日韩一区二区视频在线观看视频在线| 男人操女人黄网站| 欧美国产精品一级二级三级| 国产在视频线精品| 久久久久久人妻| 夜夜爽夜夜爽视频| 久久久久精品人妻al黑| 亚洲欧美精品自产自拍| www.色视频.com| 久久99热这里只频精品6学生| 国产精品99久久99久久久不卡 | 久久久久视频综合| 一区二区三区乱码不卡18| 热re99久久精品国产66热6| 日韩一区二区视频免费看| 国产av码专区亚洲av| 亚洲天堂av无毛| 又黄又粗又硬又大视频| 婷婷成人精品国产| 赤兔流量卡办理| 日本黄大片高清| 精品一区二区免费观看| 免费观看性生交大片5| 欧美成人午夜精品| 在线观看免费视频网站a站| 九九在线视频观看精品| av福利片在线| 久久久久久久精品精品| 欧美 日韩 精品 国产| 国产成人一区二区在线| 日韩av不卡免费在线播放| 一本—道久久a久久精品蜜桃钙片| 亚洲精品av麻豆狂野| 久热久热在线精品观看| 亚洲精华国产精华液的使用体验| 在线看a的网站| 桃花免费在线播放| 亚洲婷婷狠狠爱综合网| 精品亚洲乱码少妇综合久久| 成人亚洲精品一区在线观看| 日本午夜av视频| 久久ye,这里只有精品| 9热在线视频观看99| 国产色爽女视频免费观看| 波野结衣二区三区在线| 国产极品粉嫩免费观看在线| 观看美女的网站| 国产乱来视频区| 亚洲精品久久成人aⅴ小说| 亚洲精品aⅴ在线观看| 日本vs欧美在线观看视频| 日韩成人av中文字幕在线观看| 下体分泌物呈黄色| av在线老鸭窝| 九九在线视频观看精品| 免费播放大片免费观看视频在线观看| 视频中文字幕在线观看| 黄网站色视频无遮挡免费观看| 精品一品国产午夜福利视频| 欧美xxⅹ黑人| 777米奇影视久久| 一区二区三区乱码不卡18| 精品亚洲乱码少妇综合久久| 欧美最新免费一区二区三区| 观看av在线不卡| 亚洲,一卡二卡三卡| 女人被躁到高潮嗷嗷叫费观| 深夜精品福利| 中文欧美无线码| 久久久久人妻精品一区果冻| 在线免费观看不下载黄p国产| 国产av一区二区精品久久| 亚洲av欧美aⅴ国产| 大片电影免费在线观看免费| 成年人午夜在线观看视频| 国产精品麻豆人妻色哟哟久久| 久久国产精品大桥未久av| 日韩制服骚丝袜av| 美女内射精品一级片tv| 另类精品久久| 9191精品国产免费久久| 美女内射精品一级片tv| 久热久热在线精品观看| 一级爰片在线观看| 亚洲精华国产精华液的使用体验| a 毛片基地| 国产精品久久久久久久久免| 美女中出高潮动态图| 男女免费视频国产| 国产国拍精品亚洲av在线观看| 最近的中文字幕免费完整| 欧美激情国产日韩精品一区| 26uuu在线亚洲综合色| 国产1区2区3区精品| 乱人伦中国视频| 永久网站在线| 国产精品免费大片| 在线免费观看不下载黄p国产| 人妻一区二区av| 亚洲欧洲国产日韩| 午夜老司机福利剧场| 夜夜爽夜夜爽视频| 久久99精品国语久久久| 中国国产av一级| 在线 av 中文字幕| 高清毛片免费看| 人人妻人人澡人人爽人人夜夜| 成人国语在线视频| 曰老女人黄片| 国产不卡av网站在线观看| 成人手机av| 日日爽夜夜爽网站| 国产亚洲av片在线观看秒播厂| 国产成人精品无人区| 国产精品女同一区二区软件| 久久久久久人人人人人| 亚洲成人av在线免费| 精品午夜福利在线看| 国产精品久久久久久精品电影小说| 插逼视频在线观看| 亚洲国产最新在线播放| 亚洲成人av在线免费| 最后的刺客免费高清国语| 色网站视频免费| 国精品久久久久久国模美| av有码第一页| 在线观看免费高清a一片| 成人免费观看视频高清| 男男h啪啪无遮挡| 日韩熟女老妇一区二区性免费视频| 卡戴珊不雅视频在线播放| 黄色一级大片看看| 在线观看免费高清a一片| 成人亚洲精品一区在线观看| 99九九在线精品视频| 精品国产露脸久久av麻豆| 又黄又粗又硬又大视频| 亚洲国产精品一区二区三区在线| 亚洲欧美中文字幕日韩二区| 一级片免费观看大全| 97人妻天天添夜夜摸| 亚洲av国产av综合av卡| a级毛色黄片| 如日韩欧美国产精品一区二区三区| 欧美亚洲日本最大视频资源| 亚洲美女黄色视频免费看| 婷婷成人精品国产| av在线观看视频网站免费| 春色校园在线视频观看| 最近中文字幕2019免费版| 成人无遮挡网站| 久久av网站| 免费久久久久久久精品成人欧美视频 | 九色亚洲精品在线播放| 免费高清在线观看视频在线观看| 日韩一本色道免费dvd| 水蜜桃什么品种好| 大香蕉久久成人网| 人妻少妇偷人精品九色| 99视频精品全部免费 在线| 成人18禁高潮啪啪吃奶动态图| 97在线人人人人妻| 香蕉国产在线看| 午夜老司机福利剧场| 亚洲,欧美精品.| 十分钟在线观看高清视频www| 久热久热在线精品观看| 亚洲精品自拍成人| 亚洲欧美清纯卡通| 女人久久www免费人成看片| 成人18禁高潮啪啪吃奶动态图| 男人添女人高潮全过程视频| www.色视频.com| 飞空精品影院首页| 蜜臀久久99精品久久宅男| 美女国产高潮福利片在线看| 成人亚洲欧美一区二区av| 啦啦啦中文免费视频观看日本| 免费看不卡的av| 久久人人爽人人爽人人片va| 精品99又大又爽又粗少妇毛片| 日韩大片免费观看网站| 99热全是精品| 日本黄大片高清| 咕卡用的链子| 黄色一级大片看看| 日本av手机在线免费观看| 两个人免费观看高清视频| 大码成人一级视频| 91成人精品电影| 国产麻豆69| 中文字幕制服av| 最黄视频免费看| 欧美丝袜亚洲另类| 亚洲一区二区三区欧美精品| 五月天丁香电影| 最近最新中文字幕大全免费视频 | 人妻系列 视频| 在线看a的网站| 亚洲国产最新在线播放| 久久久久国产网址| 久久综合国产亚洲精品| 人人妻人人爽人人添夜夜欢视频| 女人精品久久久久毛片| 高清毛片免费看| 人人妻人人澡人人看| 亚洲精品中文字幕在线视频| 国产极品天堂在线| 国产精品熟女久久久久浪| 亚洲成色77777| 久久国产亚洲av麻豆专区| 成年人免费黄色播放视频| 成人影院久久| 国产精品人妻久久久久久| 母亲3免费完整高清在线观看 | 亚洲精华国产精华液的使用体验| 免费久久久久久久精品成人欧美视频 | 97在线人人人人妻| 国产成人91sexporn| 女人被躁到高潮嗷嗷叫费观| 街头女战士在线观看网站| 欧美xxxx性猛交bbbb| 交换朋友夫妻互换小说| 性色avwww在线观看| 国产成人精品福利久久| 91久久精品国产一区二区三区| 在线观看免费视频网站a站| 国产成人免费观看mmmm| 中文字幕制服av| 国产极品粉嫩免费观看在线| 如何舔出高潮| 婷婷色综合大香蕉| 一二三四在线观看免费中文在 | 大片电影免费在线观看免费| 国产xxxxx性猛交| 国产国语露脸激情在线看| 在线观看人妻少妇| 亚洲国产欧美在线一区| 国产一区有黄有色的免费视频| 黑人高潮一二区| 一级爰片在线观看| 精品午夜福利在线看| 久久婷婷青草| 黑人猛操日本美女一级片| 日韩一区二区视频免费看| av一本久久久久| 国产成人一区二区在线| 国产日韩一区二区三区精品不卡| 午夜免费男女啪啪视频观看| 亚洲综合色网址| 男人添女人高潮全过程视频| 亚洲精品色激情综合| 久久国产精品男人的天堂亚洲 | 美女中出高潮动态图| 亚洲欧洲日产国产| 另类亚洲欧美激情| 男女高潮啪啪啪动态图| 人妻一区二区av| 青青草视频在线视频观看| 大香蕉久久网| 9191精品国产免费久久| 美女脱内裤让男人舔精品视频| 美女内射精品一级片tv| 黄色一级大片看看| 亚洲精品久久午夜乱码| 日韩,欧美,国产一区二区三区| 久久精品久久精品一区二区三区| 性高湖久久久久久久久免费观看| 一二三四中文在线观看免费高清| 亚洲久久久国产精品| 国产精品女同一区二区软件| 精品人妻熟女毛片av久久网站| 久久综合国产亚洲精品| 91精品三级在线观看| 啦啦啦视频在线资源免费观看| 伦理电影大哥的女人| 啦啦啦视频在线资源免费观看| 最近中文字幕2019免费版| 女性被躁到高潮视频| 1024视频免费在线观看| 飞空精品影院首页| 一区二区av电影网| 一个人免费看片子| 少妇的逼好多水| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久这里只有精品19| 久久久久久人妻| 亚洲av男天堂| videos熟女内射| 免费久久久久久久精品成人欧美视频 | 国产午夜精品一二区理论片| 搡老乐熟女国产| 日韩,欧美,国产一区二区三区| 午夜福利视频精品| 波野结衣二区三区在线| 美国免费a级毛片| 欧美精品一区二区大全| 少妇熟女欧美另类| 中文欧美无线码| 日日啪夜夜爽| 亚洲内射少妇av| 免费在线观看完整版高清| 校园人妻丝袜中文字幕| 最近中文字幕2019免费版| 亚洲国产精品成人久久小说| 国产极品天堂在线| 国产精品蜜桃在线观看| 免费在线观看黄色视频的| 国产乱人偷精品视频| 精品午夜福利在线看| 97超碰精品成人国产| 国精品久久久久久国模美| 欧美日韩一区二区视频在线观看视频在线| 亚洲国产色片| 亚洲精品乱码久久久久久按摩| 丝袜在线中文字幕| 十分钟在线观看高清视频www| 日韩欧美一区视频在线观看| 两个人看的免费小视频| 国产在线视频一区二区| 日韩免费高清中文字幕av| 永久免费av网站大全| 热re99久久精品国产66热6| 少妇人妻 视频| 巨乳人妻的诱惑在线观看| 丰满乱子伦码专区| 成年人免费黄色播放视频| 欧美人与性动交α欧美软件 | 中文字幕另类日韩欧美亚洲嫩草| 免费黄色在线免费观看| 69精品国产乱码久久久| 又黄又爽又刺激的免费视频.| 免费av中文字幕在线| 欧美日本中文国产一区发布| 王馨瑶露胸无遮挡在线观看| 男女免费视频国产| 99热这里只有是精品在线观看| 看非洲黑人一级黄片| 丝袜喷水一区| 国产高清国产精品国产三级| 伦理电影大哥的女人| 女性被躁到高潮视频| 黄片无遮挡物在线观看| 新久久久久国产一级毛片| 成人影院久久|