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

    毫米級球狀沙塵顆粒連續(xù)沖擊損傷的力學(xué)模型與數(shù)值模擬研究

    2023-02-12 08:19:14游文俊王計(jì)真楊強(qiáng)詹志新胡偉平
    航空科學(xué)技術(shù) 2023年1期
    關(guān)鍵詞:沙粒凹坑沖擊

    游文俊,王計(jì)真,楊強(qiáng),詹志新,胡偉平

    1.北京航空航天大學(xué),北京 100191

    2.中國飛機(jī)強(qiáng)度研究所 結(jié)構(gòu)沖擊動(dòng)力學(xué)航空科技重點(diǎn)實(shí)驗(yàn)室,陜西 西安 710065

    鈦合金是航空、航天等領(lǐng)域重要的結(jié)構(gòu)材料,航空發(fā)動(dòng)機(jī)葉片等多采用鈦合金材料。在沙漠環(huán)境中,航空發(fā)動(dòng)機(jī)葉片容易受到空氣中沙塵微粒的連續(xù)沖擊作用,產(chǎn)生凹陷、撕裂、微裂紋等損傷,從而影響其高周疲勞性能。因此,有必要開展沙塵顆粒連續(xù)沖擊損傷的力學(xué)模型與數(shù)值模擬研究。

    目前,國內(nèi)外學(xué)者大多采用基于試驗(yàn)的分析方法和數(shù)值模擬方法研究外物對結(jié)構(gòu)的沖擊作用,并取得了一定的成果。外物沖擊的試驗(yàn)分析方法主要有機(jī)械加工缺口法、低速沖擊法、準(zhǔn)靜態(tài)擠壓法和高速彈道沖擊法。低速沖擊法包括擺錘法和螺線槍法,高速沖擊法包括槍彈法和空氣炮法。J. Dunham等[1]利用機(jī)械加工缺口的方法分析了三種金屬材料對缺口的敏感程度。Hu Dianyin 等[2]對鈦合金葉片邊緣采用落錘沖擊的方法模擬外物沖擊,并通過電鏡掃描分析沖擊損傷部位,發(fā)現(xiàn)沖擊切口具有大塑性變形、材料堆積、材料褶皺、階梯解理、液滴結(jié)構(gòu)和微裂紋等典型特征。Mall 等[3]用靜態(tài)擠壓法在鈦合金平板上模擬外物沖擊,研究結(jié)果表明,準(zhǔn)靜態(tài)擠壓方法能模擬出明顯的塑性變形,損傷深度越大,越容易出現(xiàn)剪切裂紋。B. Lin等[4]利用高速彈道沖擊法模擬了立方體沖擊物沖擊激光噴丸強(qiáng)化后的翼型葉片前緣,并分析了沖擊殘余應(yīng)力的分布規(guī)律。M.W.Crowell 等[5]利用氣槍進(jìn)行外物沖擊葉片涂層試驗(yàn),研究外物沖擊對涂層的破壞作用。Nesredin 等[6]利用改進(jìn)的氣槍進(jìn)行外物沖擊有涂層窄葉片前沿的試驗(yàn),分析了沖擊對葉片造成的損傷。吳松波等[7]在新型沖擊磨損試驗(yàn)機(jī)上開展了硬質(zhì)沙粒沖擊TC4鈦合金的試驗(yàn),研究了沙粒粒徑、沖擊次數(shù)對TC4鈦合金沖擊磨損行為的影響,研究結(jié)果表明,粒徑和沖擊次數(shù)能顯著影響TC4 鈦合金的沖擊磨損行為。Yan Cheng 等[8]利用帶拉瓦爾噴嘴的空氣炮裝置研究了沙粒高速沖擊Ti-6Al-4V葉片,揭示了葉片破壞與沙粒沖擊速度和沖擊角度的關(guān)系。Wheeler 等[9]利用高速沖擊裝置研究了沙粒沖擊金剛石涂層的損傷機(jī)制,通過電鏡掃描損傷部位,發(fā)現(xiàn)周向裂紋和微空洞是主要破壞特征。張宇等[10]基于二級空氣炮系統(tǒng)對V形鋁合金靶板進(jìn)行了高速沖擊試驗(yàn),發(fā)現(xiàn)其破壞模式為局部剪切沖塞,射彈角和射彈間距對塑性區(qū)大小有顯著影響。王計(jì)真[11]采用一種初始應(yīng)力復(fù)合材料高速沖擊試驗(yàn)方法,對復(fù)合材料層合板開展了高速沖擊試驗(yàn),研究發(fā)現(xiàn)面內(nèi)初始應(yīng)力對復(fù)合材料板高速沖擊行為影響顯著。

    數(shù)值模擬方法的相關(guān)研究主要是基于有限元方法的仿真計(jì)算。韓璟圳[12]在ABAQUS平臺利用J-C本構(gòu)模型和損傷演化模型,進(jìn)行外物沖擊輕型點(diǎn)陣結(jié)構(gòu)數(shù)值模擬,并通過分析損傷面積,得出不同角度對點(diǎn)陣結(jié)構(gòu)沖擊損傷的影響。凡金金[13]進(jìn)行了沙粒碰撞平板數(shù)值模擬,并研究了沙粒碰撞后的反彈特性。唐飛燕等[14]采用計(jì)算流體力學(xué)(CFD)技術(shù)模擬風(fēng)沙場得到沙粒速度,基于動(dòng)量定理和沖量等效原則研究了沙粒對建筑物的沖擊作用。Chen[15]借助ABAQUS 軟件的顯示動(dòng)力分析模塊對鈦合金的外物沖擊損傷過程進(jìn)行了模擬,采用率相關(guān)本構(gòu)模型計(jì)算了沖擊變形過程。P. Duo等[16]利用LSDYNA軟件動(dòng)力分析模塊開展外物沖擊數(shù)值模擬,研究了外物損傷的殘余應(yīng)力分布,并且與試驗(yàn)結(jié)果進(jìn)行了對比,發(fā)現(xiàn)該數(shù)值模擬方法與試驗(yàn)結(jié)果吻合較好。D. Williams 等[17]采用 DYNA3D 軟件分析了FOD 的動(dòng)態(tài)沖擊過程,結(jié)果表明,低速沖擊條件下的數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果吻合很好,沖擊過程中只發(fā)生了材料屈服,沒有失效。P. H. Tranter 等[18]基于修正的 Zerilli-Armstrong模型對FOD進(jìn)行了模擬,研究表明Goldthorpe路徑相關(guān)失效模型能夠很好地模擬多種沖擊損傷。

    本文開展了毫米級球狀沙塵顆粒連續(xù)沖擊損傷的力學(xué)模型與數(shù)值模擬研究。首先,推導(dǎo)了損傷耦合的J-C本構(gòu)模型和連續(xù)沖擊損傷模型。其次,編寫Vufield 子程序和Python 腳本,實(shí)現(xiàn)了連續(xù)沖擊損傷計(jì)算的數(shù)值計(jì)算。通過將計(jì)算結(jié)果與文獻(xiàn)中的試驗(yàn)數(shù)據(jù)進(jìn)行對比,驗(yàn)證了該方法的有效性。最后,分別進(jìn)行單個(gè)和多個(gè)沙粒連續(xù)沖擊葉片的數(shù)值模擬,分析了沖擊變形、殘余應(yīng)力和沖擊損傷的變化規(guī)律。

    1 理論模型

    1.1 損傷耦合的J-C本構(gòu)模型

    高應(yīng)變率情況下,材料的應(yīng)力應(yīng)變關(guān)系通常采用率相關(guān)本構(gòu)方程來描述,其中Johnson-Cook 方程就是一種很好的描述大多數(shù)金屬動(dòng)態(tài)塑性行為的本構(gòu)模型,其方程形式如式(1)所示

    式中,σeq為等效應(yīng)力;ε為等效塑性應(yīng)變;ε?為應(yīng)變率;T為工件當(dāng)前溫度;ε?0為基準(zhǔn)應(yīng)變率;Tr和Tm分別為室溫和材料熔化溫度;A、B、C、m、n是材料常數(shù)。

    沖擊會造成損傷,造成材料劣化,沖擊引起的材料劣化可以用連續(xù)損傷力學(xué)中的損傷變量來表示。對于各向同性損傷,損傷變量D可以表示為代表性體積單元(RVE)剛度的劣化率。

    式中:E為無損傷的彈性模量;ED為含損傷RVE 的彈性模量。

    則耦合損傷的J-C本構(gòu)模型表示如下

    1.2 連續(xù)沖擊損傷模型

    基于上述損傷耦合的J-C本構(gòu)模型所計(jì)算得到的沖擊應(yīng)力應(yīng)變,可根據(jù)式(5)計(jì)算由塑性變形引起的損傷

    若沙塵顆粒連續(xù)沖擊i次,設(shè)第i次沖擊結(jié)束后,損傷為Di,在進(jìn)行第i次沖擊時(shí),將第i-1 次沖擊結(jié)束后的場變量(包括應(yīng)力應(yīng)變)作為第i次沖擊的初始場,并將第i-1次的損傷Di-1作為第i次沖擊的初始損傷代入式(4)中,可計(jì)算第i次沖擊的等效應(yīng)力和等效塑性應(yīng)變,同時(shí)將等效應(yīng)力和等效塑性應(yīng)變以及Di-1代入式(7)可計(jì)算損傷增量ΔDi,然后由式(8)得到第i次沙塵顆粒沖擊后的損傷Di。進(jìn)而得到第i次沖擊后的殘余應(yīng)力、殘余應(yīng)變及沖擊損傷

    2 材料參數(shù)標(biāo)定

    有兩類材料參數(shù)需要標(biāo)定,即J-C本構(gòu)模型參數(shù)和連續(xù)沖擊損傷模型參數(shù)。本文葉片的材料為Ti-6Al-4V鈦合金,其基本力學(xué)性能參數(shù)見表1[2],該材料的J-C 本構(gòu)模型的參數(shù)取自參考文獻(xiàn)[2],見表2。根據(jù)Ti-6Al-4V鈦合金試驗(yàn)件的低周疲勞試驗(yàn),采用最小二乘法可確定沖擊損傷模型中的參數(shù)S和m,見表3。

    表1 Ti-6Al-4V靜力性能參數(shù)[2]Table 1 Static mechanical properties of Ti-6Al-4V[2]

    表2 J-C本構(gòu)模型參數(shù)的材料參數(shù)[2]Table 2 Material parameters of the J-C constitutive model[2]

    表3 Ti-6Al-4V沖擊損傷計(jì)算模型參數(shù)[19]Table 3 Parameters of Ti-6Al-4V impact damage calculation model[19]

    3 連續(xù)沖擊損傷的數(shù)值模擬

    3.1 數(shù)值計(jì)算方法

    基于ABAQUS 平臺,實(shí)現(xiàn)損傷耦合的J-C 本構(gòu)模型和沖擊損傷模型的有限元數(shù)值計(jì)算。計(jì)算流程簡述如下:(1)初始化所有參數(shù)。(2)在ABAQUS 平臺建立有限元模型,包括劃分網(wǎng)格、材料屬性設(shè)置、葉片邊界條件設(shè)置、接觸屬性設(shè)置(包括接觸算法設(shè)置和摩擦系數(shù)設(shè)置)、設(shè)置分析步為顯式動(dòng)力學(xué)。此外,在設(shè)置分析步時(shí),分析步時(shí)間的設(shè)置應(yīng)滿足沖擊物與葉片接觸后完全分離。(3)將損傷耦合的J-C本構(gòu)模型和沖擊損傷模型編寫為Vufield 子程序并嵌入ABAQUS Explicit Dynamic 模塊計(jì)算應(yīng)力應(yīng)變以及沖擊損傷。(4)前一次沖擊結(jié)束后,利用Python 程序?qū)⒃摯萎a(chǎn)生的損傷、殘余應(yīng)力應(yīng)變和變形導(dǎo)入下一次沖擊計(jì)算模型作為初始場。利用Python程序設(shè)置工況,包括沙粒尺寸、沙粒個(gè)數(shù)、沖擊速度、沖擊角度、連續(xù)沖擊次數(shù), 以實(shí)現(xiàn)不同工況的連續(xù)沖擊。

    3)1985-2000年,巢湖流域的8個(gè)評價(jià)單元全部為負(fù)值,表明土地利用程度衰退是該時(shí)間段內(nèi)巢湖流域的主要趨勢,人類活動(dòng)對土地利用程度的影響較小。

    3.2 數(shù)值模擬方法的有效性驗(yàn)證

    本文采用一次沖擊的試驗(yàn)數(shù)據(jù)對上述數(shù)值模擬方法進(jìn)行驗(yàn)證。雖然沙粒和鋼球的特性差異很大,但是二者的硬度均遠(yuǎn)大于鈦合金,本文采用了鋼球撞擊試驗(yàn)來驗(yàn)證所提的分析方法。 J. O. Peters等[20]利用直徑為3.2mm的鋼球在高速(200m/s、250m/s、300m/s)下垂直沖擊Ti-6Al-4V 試件,得到的沖擊凹坑的結(jié)果見表4。該試驗(yàn)中,鋼球的密度為7850kg/m3,試件的沖擊示意圖如圖1所示。

    表4 沖擊凹坑深度的數(shù)值模擬結(jié)果Table 4 Numerical simulation results of pit depth

    圖1 沖擊示意圖[19]Fig. 1 Ⅰmpact diagram[19]

    首先,建立有限元模型,采用網(wǎng)格局部細(xì)化法對模型劃分網(wǎng)格,如圖2 所示,模型共有120600 個(gè)單元、127402 個(gè)節(jié)點(diǎn),單元最小尺寸為0.1mm。邊界條件為兩端固定。材料屬性方面,鋼珠設(shè)置為剛體,鈦合金葉片的材料參數(shù)見表1 和表2。其次,設(shè)置接觸屬性,采用通用接觸,摩擦系數(shù)設(shè)置為0.15。然后,設(shè)置分析步時(shí)間為1e-5s,以保證鋼珠和葉片接觸后可完全分離。最后,編寫Vufield 子程序并嵌入ABAQUS Explicit Dynamic模塊,實(shí)現(xiàn)沖擊模擬,得到?jīng)_擊響應(yīng)。本文施加的沖擊物的沖擊速度分別為200m/s、250m/s和300m/s,沖擊凹坑的云圖如圖3所示,提取凹坑最大深度的結(jié)果見表4,可以發(fā)現(xiàn),數(shù)值模擬與試驗(yàn)結(jié)果相近,誤差在允許范圍內(nèi),驗(yàn)證了所提數(shù)值模擬方法的適用性。此外,對于300m/s 的沖擊速度,計(jì)算得到如圖4 所示的初始沖擊損傷分布云圖,將其與參考文獻(xiàn)[20]中的實(shí)驗(yàn)觀測的電鏡圖(見圖5)進(jìn)行對比,發(fā)現(xiàn)電鏡圖中的微裂紋分布規(guī)律與初始損傷云圖的分布規(guī)律有一定的相似性。

    圖2 有限元模型圖Fig.2 Finite element model

    圖3 沖擊變形圖Fig.3 Ⅰmpact deformation diagram

    圖4 300m/s時(shí)的沖擊損傷云圖Fig.4 Ⅰmpact damage diagram with 300m/s

    圖5 沖擊凹坑電鏡圖[20]Fig.5 Electron microscope of impact pit [20]

    4 單個(gè)沙粒連續(xù)沖擊葉片數(shù)值模擬

    4.1 有限元模型

    本節(jié)利用上述方法開展單個(gè)沙粒連續(xù)沖擊葉片的數(shù)值模擬研究。沙粒半徑在0.3~1mm[21],沙粒的基本力學(xué)性能見表5。考慮到飛機(jī)發(fā)動(dòng)機(jī)葉片容易受到沙粒不同角度的高速沖擊,因此,為了研究沙粒尺寸、沖擊速度、沖擊角度,以及沖擊次數(shù)對葉片的沖擊損傷影響,本節(jié)開展以下工況的數(shù)值模擬研究:沙粒半徑分別取為0.5mm、0.7mm、0.9mm,沖擊速度分別取為100m/s、150m/s、200m/s,沖擊角度分別取為15°、30°、45°、60°、75°、90°,設(shè)定沖擊次數(shù)為10次,見表6。沖擊示意圖如圖6所示。葉片形狀參考了參考文獻(xiàn)[22]中葉片的基本尺寸,如圖7 所示,長65mm,寬25mm,厚2mm。

    圖6 單個(gè)沙粒沖擊示意圖Fig.6 Ⅰmpact schematic diagram of single sand particle

    表5 沙粒力學(xué)性能參數(shù)[23]Table 5 Mechanical property parameters of sand particles[23]

    表6 單個(gè)沙粒沖擊工況Table 6 Single sand impact condition

    基于ABAQUS平臺,建立有限元模型。在葉片沖擊接觸部位進(jìn)行網(wǎng)格細(xì)化,該部分網(wǎng)格采用C3D8R 單元,最小網(wǎng)格尺寸為0.05mm。為提高計(jì)算速度,葉片非接觸區(qū)域采用C3D4單元開展網(wǎng)格劃分。考慮到沙粒的硬度遠(yuǎn)大于鈦合金的硬度,在數(shù)值模擬中將沙粒視為剛體,故未考慮沙粒的率相關(guān)的特性。沙粒采用C3D8R 單元?jiǎng)澐志W(wǎng)格。通過有限元網(wǎng)格的收斂性驗(yàn)證分析,最終建立的有限元模型如圖8所示,包含223498個(gè)單元、280486個(gè)節(jié)點(diǎn)。

    圖8 單個(gè)沙粒沖擊有限元模型Fig.8 Finite element model of single sand impact

    4.2 單個(gè)沙粒連續(xù)沖擊數(shù)值計(jì)算結(jié)果

    4.2.1 沖擊變形分析

    為了研究單個(gè)沙粒沖擊后鈦合金板的變形情況,這里取沙粒半徑為0.7mm,沖擊速度為200m/s,沖擊角度分別為30°、60°、90°的變形云圖進(jìn)行分析,如圖9所示??梢钥吹剑瑳_擊角度為30°、60°時(shí)沖擊局部會出現(xiàn)材料堆積,這是由于葉片受到了切應(yīng)力,從而產(chǎn)生滑移,最終導(dǎo)致材料堆積,如圖9(a)中A處和圖9(b)中C處所示,同時(shí)會產(chǎn)生材料凹陷即沖擊凹坑,如圖9(a)中B處和圖9(b)中D處所示;沖擊角度為90°時(shí)幾乎不會產(chǎn)生材料堆積,僅產(chǎn)生材料凹陷即沖擊凹坑,如圖9(c)所示。綜上,當(dāng)沖擊角度小于90°時(shí),沖擊局部區(qū)域產(chǎn)生材料堆積和凹陷,而垂直沖擊時(shí)僅產(chǎn)生材料凹陷。然后,分析不同角度沖擊后在路徑1(見圖9(a)Z向)上的法向位移(U2方向)分布規(guī)律,如圖9(d)所示,其中,U2>0即為材料堆積區(qū),U2<0即為材料凹陷區(qū)??梢园l(fā)現(xiàn):沖擊角度為30°時(shí),U2 的正值最大,即材料堆積高度最大;沖擊角度為90°時(shí),U2 的負(fù)值最小,即凹陷最深。本文統(tǒng)一把U2<0 的情況統(tǒng)稱為凹坑,U2 的最小值即為凹坑的深度。圖10給出了沖擊凹坑深度隨沖擊次數(shù)的變化規(guī)律,可以看到:沖擊次數(shù)越多,沖擊凹坑深度越大;隨著沖擊次數(shù)的增加,曲線越來越平緩,表明沖擊凹坑深度的變化率越來越小。圖11 給出了單個(gè)沙粒一次沖擊后凹坑深度隨沖擊速度的變化規(guī)律,可以看到,沖擊速度越大,沖擊凹坑深度越大。沖擊凹坑深度隨沖擊角度的變化情況如圖12所示,隨著沖擊角度的增加,沖擊凹坑越來越深。圖13給出了沖擊凹坑深度隨沙粒半徑的變化趨勢,沙粒半徑越大,凹坑越深。這是由于半徑越大,其質(zhì)量越大,沖擊動(dòng)能越大,從而在沙粒與葉片接觸時(shí),葉片會吸收更多的沖擊動(dòng)能轉(zhuǎn)化為塑性應(yīng)變能,意味著更大的不可恢復(fù)的變形,在宏觀上表現(xiàn)為更大尺寸的凹坑。

    圖9 不同沖擊角度葉片變形情況分析Fig.9 Analysis on blade deformation at different impact angles

    圖10 沖擊凹坑隨沖擊次數(shù)的變化Fig.10 Variation of impact pit with the number of impact

    圖11 沖擊凹坑隨沖擊速度的變化Fig.11 Variation of impact pit with impact velocity

    圖12 沖擊凹坑隨沖擊角度的變化Fig.12 Variation of impact pit with impact angle

    圖13 沖擊凹坑隨沙粒尺寸的變化Fig.13 Variation of impact pit with sand particle size

    4.2.2 殘余應(yīng)力分析

    葉片受到外物沖擊后,一部分沖擊能量以塑性變形能的形式耗散,產(chǎn)生塑性變形,在受沖擊區(qū)域產(chǎn)生較大的殘余應(yīng)力。本文所述的殘余應(yīng)力是指在移除沙粒沖擊作用后,仍留在葉片內(nèi)的自相平衡的內(nèi)應(yīng)力。本節(jié)主要分析Z方向的殘余應(yīng)力(S33)。圖14 展示了沙粒半徑為0.5mm,沖擊速度為150m/s,沖擊角度分別為30°、60°、90°的Z向殘余應(yīng)力分布云圖,可以看到,沖擊后產(chǎn)生了殘余壓應(yīng)力和殘余拉應(yīng)力,且殘余壓應(yīng)力主要出現(xiàn)在次表層。沖擊角度為30°時(shí),最大殘余壓應(yīng)力值為-1308MPa,如圖14(a)所示;沖擊角度為60°時(shí),最大殘余壓應(yīng)力值為-1682MPa,如圖14(b)所示;沖擊角度為90°時(shí),最大殘余壓應(yīng)力值為-1806MPa,如圖14(c)所示。對比三種情況的最大殘余壓應(yīng)力,發(fā)現(xiàn)沖擊角度為90°時(shí),殘余壓應(yīng)力最大,出現(xiàn)在次表層。

    圖15展示了圖14中路徑H上的殘余應(yīng)力分布規(guī)律,可以發(fā)現(xiàn),沖擊角度為90°時(shí),殘余應(yīng)力在路徑H 上呈現(xiàn)對稱分布,且殘余應(yīng)力影響的跨度最大,而沖擊角度為30°時(shí),殘余應(yīng)力影響的跨度最小。圖16給出了沙粒半徑為0.7mm,沖擊速度為150m/s,沖擊角度為90°時(shí),不同沖擊次數(shù)后路徑F上的殘余應(yīng)力分布規(guī)律,可以看到,在路徑F上殘余壓應(yīng)力先變大后變小,最后變成拉伸應(yīng)力再減小趨于0。此外,一次沖擊后,殘余應(yīng)力分布曲線基本重合,這說明一次沖擊后殘余應(yīng)力變化率變小,但隨著沖擊次數(shù)的增加,殘余應(yīng)力影響的跨度變大。

    圖14 不同沖擊角度葉片Z向殘余應(yīng)力分布Fig.14 Z-direction residual stress distribution of blades ith different impact angles

    圖15 不同角度沖擊下路徑H上的軸向殘余應(yīng)力分布Fig.15 Axial residual stress distribution on path H with different impact angles

    圖16 沖擊角度為90°時(shí),不同沖擊次數(shù)下路徑F上的軸向殘余應(yīng)力分布Fig.16 Axial residual stress distribution on path F with different number of impacts under the impact angle of 90°

    為了研究殘余應(yīng)力與沖擊速度的關(guān)系,圖17給出了沖擊角度為90°,沙粒半徑為0.5mm,一次沖擊速度分別為100m/s、150m/s、200m/s 時(shí)的路徑F 上的殘余應(yīng)力分布曲線??梢钥吹?,沖擊速度越大,曲線越靠下且曲線跨度越大,即接觸部位產(chǎn)生的整體殘余壓應(yīng)力越大且殘余應(yīng)力影響跨度越大。圖18 給出了不同沙粒半徑對殘余應(yīng)力的影響,發(fā)現(xiàn)沙粒半徑越大,殘余壓應(yīng)力越大且殘余應(yīng)力影響跨度越大。

    圖17 不同沖擊速度下路徑F的軸向殘余應(yīng)力分布Fig. 17 Axial residual stress distribution on path F with different impact speeds

    圖18 不同沙粒半徑下路徑F的軸向殘余應(yīng)力分布Fig. 18 Axial residual stress distribution on path F with different sand radius

    4.2.3 沖擊損傷分析

    本節(jié)所述的沖擊損傷是由損傷度D 表征,基于前述的沖擊損傷模型,可計(jì)算沙塵顆粒一次沖擊所產(chǎn)生的損傷增量,可以得到單個(gè)沙粒在不同沖擊角度(30°、60°及90°)下的沖擊損傷分布云圖。如圖19 所示,可以看到,當(dāng)沖擊角度為30°時(shí),最大損傷出現(xiàn)在葉片表面,這是由于當(dāng)沖擊角度為30°時(shí),會產(chǎn)生切削作用。圖20 展示了沙粒半徑為0.7mm,沖擊速度為200m/s,且同一次沖擊后,沖擊損傷最大值隨著沖擊次數(shù)的變化規(guī)律,可以發(fā)現(xiàn),沖擊次數(shù)越多,沖擊損傷最大值越大,但隨著沖擊次數(shù)的增加,沖擊損傷最大值的變化率越來越小,這是由于沖擊在產(chǎn)生損傷的同時(shí)也會使材料得到硬化,這種硬化減緩了沖擊損傷的增加。圖21給出了沖擊損傷與沖擊速度的關(guān)系,從圖21中可以看出,沖擊速度越大,沖擊損傷越大。圖22 給出了沙粒半徑為0.9mm、沖擊速度為150m/s情況下,不同沖擊角度下沖擊損傷的變化規(guī)律,可以看出沖擊損傷在30°時(shí)最大。圖23展示了沖擊速度為200m/s、沖擊角度為90°時(shí),不同沙粒尺寸造成的沖擊損傷隨沖擊次數(shù)的變化規(guī)律,可以發(fā)現(xiàn),沙粒半徑越大,產(chǎn)生的損傷越大。

    圖19 不同沖擊角度損傷分布Fig.19 Damage distribution at different impact angles

    圖20 不同沖擊次數(shù)下最大損傷變化規(guī)律Fig.20 Variation of maximum damage versus different number of impacts

    圖21 不同沖擊速度下最大損傷變化規(guī)律Fig. 21 Variation of maximum damage versus different impact speeds

    圖22 不同沖擊角度下最大損傷變化規(guī)律Fig.22 Variation of maximum damage versus different impact angles

    圖23 不同沙粒尺寸下最大損傷變化規(guī)律Fig.23 Variation of maximum damage versus different sand particle size

    5 多個(gè)沙粒連續(xù)沖擊葉片數(shù)值模擬分析

    5.1 有限元模型

    在沙漠環(huán)境中,葉片會同時(shí)受到多個(gè)沙粒的沖擊作用,本節(jié)開展多個(gè)沙粒連續(xù)沖擊葉片的數(shù)值模擬分析,研究沖擊損傷與沙粒個(gè)數(shù)及沖擊位置的關(guān)系。為了更加真實(shí)有效地模擬多個(gè)沙粒的沖擊作用,這里考慮沙粒沖擊位置的隨機(jī)分布。設(shè)置沙粒個(gè)數(shù)分別為10個(gè)、20個(gè)、30個(gè)、40個(gè),同時(shí)考慮沙粒沖擊時(shí)沖擊位置的隨機(jī)性,分兩種情況,即第一次沖擊位置隨機(jī)后不變和每次沖擊位置都隨機(jī)。沙粒沖擊速度為150m/s,沙粒半徑為0.5mm,沖擊角度為90°,計(jì)算模擬工況見表7。多個(gè)沙粒沖擊示意圖如圖24所示。

    表7 多個(gè)沙粒沖擊工況說明Table 7 Impact conditions of multiple sand

    圖24 多個(gè)沙粒沖擊示意圖Fig.24 Schematic diagram of multiple sand impact

    建立多個(gè)沙粒沖擊的有限元模型如圖25 所示,采用C3D8R 單元,在主要接觸部位細(xì)化網(wǎng)格,通過有限元網(wǎng)格的收斂性驗(yàn)證分析,最終,該模型共包含280097 個(gè)單元、387838個(gè)節(jié)點(diǎn),最小單元尺寸為0.05mm。模型的邊界條件為一端固支,沙粒與鈦合金板之間設(shè)為通用接觸,摩擦系數(shù)為0.15。為保證沙粒與葉片沖擊后完全分離,本文設(shè)置的時(shí)間間隔為10-5s。

    圖25 多個(gè)沙粒沖擊的有限元模型圖Fig.25 Finite element models of multiple sand impact

    5.2 多個(gè)沙粒連續(xù)沖擊數(shù)值計(jì)算結(jié)果

    5.2.1 沖擊變形分析

    對于兩種不同的隨機(jī)沖擊方式,30個(gè)沙粒連續(xù)沖擊10次后的變形云圖如圖26 所示。圖27 給出了兩種不同隨機(jī)沖擊方式?jīng)_擊下,最大沖擊凹坑深度與沙粒個(gè)數(shù)、沖擊次數(shù)的變化規(guī)律,紅線代表每次沖擊位置隨機(jī)的情況,藍(lán)線代表第一次沖擊位置隨機(jī)后不變的情況。可以發(fā)現(xiàn),沙粒個(gè)數(shù)越多,最大沖擊凹坑深度越深,且隨著沖擊次數(shù)增加,沖擊凹坑深度越大,與此同時(shí),最大凹坑深度的變化率也增大。對比這兩種不同的沖擊位置隨機(jī)分布情況,可以看到,紅線均位于藍(lán)線下方,且紅線比藍(lán)線更加陡峭,這說明每次沖擊位置隨機(jī)的沖擊方式產(chǎn)生的最大沖擊凹坑深度更深,可以推斷隨機(jī)沖擊會減弱沖擊后的硬化作用。

    圖26 30個(gè)沙粒連續(xù)沖擊10次后變形云圖Fig.26 Deformation diagram after 10 consecutive impacts of 30 sand particles

    圖27 最大凹坑深度與沙粒個(gè)數(shù)和沖擊次數(shù)的關(guān)系Fig.27 Relationship between the maximum pit depth and the number of sand particles and impact times

    5.2.2 沖擊損傷分析

    對于兩種不同的隨機(jī)沖擊方式,30個(gè)沙粒連續(xù)沖擊10次后所造成的沖擊損傷云圖如圖28 所示,可以發(fā)現(xiàn),沖擊損傷主要發(fā)生在沙粒與葉片接觸部位。圖29 給出了兩種不同隨機(jī)沖擊方式?jīng)_擊下,最大沖擊損傷與沙粒個(gè)數(shù)、沖擊次數(shù)的變化規(guī)律。紅線代表每次沖擊位置隨機(jī)的情況,藍(lán)線代表第一次沖擊位置隨機(jī)后不變的情況??梢钥吹?,對于藍(lán)線,沙粒個(gè)數(shù)越多,最大沖擊損傷越大,且第二次沖擊之后,不同沙粒個(gè)數(shù)的最大沖擊損傷有明顯差異,與此同時(shí),這種差異隨著沖擊次數(shù)的增加而增加;對于紅線,沖擊損傷隨沖擊次數(shù)變化呈階段性突增,且沙粒個(gè)數(shù)為20 個(gè)、30 個(gè)、40 個(gè)時(shí),每次沖擊位置隨機(jī)的沖擊方式造成的損傷大于第一次沖擊位置隨機(jī)后不變的沖擊方式,這說明隨機(jī)沖擊會減弱強(qiáng)化作用,更容易造成沖擊損傷。

    圖28 30個(gè)沙粒連續(xù)10次沖擊后沖擊損傷云圖Fig.28 Ⅰmpact damage diagram after 10 consecutive impacts of 30 sand particles

    圖29 最大沖擊損傷與沙粒個(gè)數(shù)和沖擊次數(shù)的關(guān)系Fig.29 Relationship between the maximum impact damage and the number of sand particles and impact times

    6 結(jié)論

    本文建立了毫米級球狀沙塵顆粒連續(xù)沖擊損傷的力學(xué)模型與數(shù)值模擬方法,并以鈦合金葉片為例,進(jìn)行了單個(gè)及多個(gè)沙粒連續(xù)沖擊的變形、殘余應(yīng)力及損傷分析。主要結(jié)論如下:

    (1)基于ABAQUS 平臺,通過編寫Vufield 子程序和Python 腳本,實(shí)現(xiàn)了連續(xù)沖擊損傷計(jì)算的數(shù)值模擬。通過將計(jì)算結(jié)果與參考文獻(xiàn)中的試驗(yàn)數(shù)據(jù)進(jìn)行對比,驗(yàn)證了該方法的有效性。

    (2)基于建立的連續(xù)沖擊模擬方法,開展了單個(gè)沙粒連續(xù)沖擊葉片數(shù)值模擬分析。沖擊凹坑深度、沖擊殘余應(yīng)力的影響跨度及沖擊損傷最大值與沙粒尺寸、沖擊速度、沖擊次數(shù)成正相關(guān)。沖擊角度為30°時(shí),沖擊損傷最大。此外,隨著沖擊次數(shù)的增加,沖擊凹坑深度、殘余應(yīng)力跨度和沖擊損傷的變化率減小,主要是由于前一次沖擊會對材料局部產(chǎn)生一定的硬化作用,進(jìn)而減緩沖擊損傷增加的速率。

    (3)針對兩種不同的隨機(jī)沖擊方式,即每次沖擊位置隨機(jī)情況和第一次沖擊位置隨機(jī)后不變的情況,開展了多個(gè)沙粒連續(xù)沖擊葉片的數(shù)值模擬分析。沙粒個(gè)數(shù)越多,最大沖擊凹坑深度越深,且隨著沖擊次數(shù)增加,沖擊凹坑深度越大,與此同時(shí),最大凹坑深度的變化率也增大。此外,沖擊損傷主要發(fā)生在沙粒與葉片接觸部位。沖擊損傷隨沖擊次數(shù)變化呈階段性突增,每次沖擊位置隨機(jī)的沖擊方式造成的損傷大于第一次沖擊位置隨機(jī)后不變的沖擊方式。

    猜你喜歡
    沙粒凹坑沖擊
    沙粒和水珠
    腐蝕凹坑干涉效應(yīng)與疲勞壽命計(jì)算
    含有不同間距凹坑缺陷的發(fā)酵罐應(yīng)力分析與計(jì)算
    想看山的小沙粒
    想看山的小沙粒
    含凹坑缺陷結(jié)構(gòu)疲勞壽命的快速估算方法
    核技術(shù)(2016年4期)2016-08-22 09:05:32
    奧迪Q5換擋沖擊
    奧迪A8L換擋沖擊
    一汽奔騰CA7165AT4尊貴型車換擋沖擊
    巴菲特給我沖擊最大
    女人久久www免费人成看片| 免费观看a级毛片全部| 亚洲精品国产av蜜桃| 啦啦啦视频在线资源免费观看| 又大又黄又爽视频免费| 国产欧美亚洲国产| 久久精品久久久久久久性| 久久人人爽av亚洲精品天堂| 国产熟女午夜一区二区三区| 久久久久国产网址| 另类精品久久| 综合色丁香网| 制服诱惑二区| videosex国产| 精品亚洲乱码少妇综合久久| av卡一久久| 亚洲伊人久久精品综合| 亚洲人成77777在线视频| 青春草国产在线视频| 天堂8中文在线网| 免费播放大片免费观看视频在线观看| 久久久久精品人妻al黑| 不卡av一区二区三区| 亚洲av福利一区| 午夜日本视频在线| 人妻一区二区av| 韩国高清视频一区二区三区| 日韩成人av中文字幕在线观看| 国产欧美日韩综合在线一区二区| av有码第一页| 熟女电影av网| 中文字幕亚洲精品专区| 久久久久久久精品精品| 爱豆传媒免费全集在线观看| 国产极品天堂在线| 成年动漫av网址| 免费黄色在线免费观看| 久久久久久久久免费视频了| 成人国产麻豆网| 永久网站在线| 亚洲欧美成人精品一区二区| 18禁观看日本| 嫩草影院入口| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品自拍成人| 成年人免费黄色播放视频| 女人高潮潮喷娇喘18禁视频| 欧美日韩av久久| 久久99热这里只频精品6学生| 精品一区二区三区四区五区乱码 | 黄色一级大片看看| 国产亚洲精品第一综合不卡| 97人妻天天添夜夜摸| 综合色丁香网| 国产精品亚洲av一区麻豆 | 国产精品av久久久久免费| 日本欧美视频一区| 日韩 亚洲 欧美在线| 久久综合国产亚洲精品| 亚洲伊人久久精品综合| 男女边摸边吃奶| 久久久久国产网址| 亚洲国产毛片av蜜桃av| 一区二区三区激情视频| 午夜91福利影院| 成人影院久久| 精品久久久久久电影网| 日韩人妻精品一区2区三区| 久久久久久久久久久免费av| 天天影视国产精品| 999精品在线视频| 日韩免费高清中文字幕av| 搡老乐熟女国产| 日韩一卡2卡3卡4卡2021年| 一级,二级,三级黄色视频| 国产在线免费精品| 久久精品aⅴ一区二区三区四区 | 一区福利在线观看| 狠狠精品人妻久久久久久综合| 97人妻天天添夜夜摸| 嫩草影院入口| 制服人妻中文乱码| 成人国语在线视频| 成人手机av| 欧美激情高清一区二区三区 | 午夜福利视频在线观看免费| 国产亚洲一区二区精品| 精品国产国语对白av| 国产精品一区二区在线不卡| 亚洲av.av天堂| 国产一区二区三区av在线| 波野结衣二区三区在线| 亚洲精品久久成人aⅴ小说| 亚洲国产毛片av蜜桃av| 建设人人有责人人尽责人人享有的| 天天躁日日躁夜夜躁夜夜| 中文字幕精品免费在线观看视频| 欧美日韩国产mv在线观看视频| 精品少妇内射三级| 不卡av一区二区三区| 老司机亚洲免费影院| 热99国产精品久久久久久7| 亚洲一区中文字幕在线| 亚洲国产精品成人久久小说| 狂野欧美激情性bbbbbb| 国产精品嫩草影院av在线观看| 亚洲精品日本国产第一区| 男人爽女人下面视频在线观看| 91精品国产国语对白视频| 啦啦啦在线免费观看视频4| 久久 成人 亚洲| 9色porny在线观看| av天堂久久9| 亚洲精品av麻豆狂野| 国产1区2区3区精品| 18禁裸乳无遮挡动漫免费视频| 97在线视频观看| 欧美日韩亚洲国产一区二区在线观看 | 国产精品亚洲av一区麻豆 | 午夜日本视频在线| av网站免费在线观看视频| 成人毛片a级毛片在线播放| 亚洲国产精品一区二区三区在线| 丝袜美腿诱惑在线| 永久免费av网站大全| 日韩一本色道免费dvd| 亚洲av日韩在线播放| 你懂的网址亚洲精品在线观看| 亚洲av电影在线进入| 九草在线视频观看| 人妻一区二区av| 久久亚洲国产成人精品v| 精品人妻在线不人妻| 日本91视频免费播放| 亚洲三级黄色毛片| 两个人免费观看高清视频| 日韩在线高清观看一区二区三区| 巨乳人妻的诱惑在线观看| 精品人妻在线不人妻| 成年女人在线观看亚洲视频| 丁香六月天网| 女人高潮潮喷娇喘18禁视频| 18禁动态无遮挡网站| 国产精品久久久久久久久免| 美女视频免费永久观看网站| 午夜激情av网站| 国产老妇伦熟女老妇高清| 亚洲国产欧美网| xxxhd国产人妻xxx| 美女xxoo啪啪120秒动态图| 国产精品女同一区二区软件| 多毛熟女@视频| av免费在线看不卡| 精品午夜福利在线看| 亚洲av福利一区| 国产麻豆69| 99香蕉大伊视频| 欧美精品一区二区免费开放| 精品福利永久在线观看| 熟女电影av网| 五月开心婷婷网| 亚洲成人av在线免费| 久久久久久久久久人人人人人人| 久久女婷五月综合色啪小说| 丝袜美腿诱惑在线| 亚洲av电影在线进入| 岛国毛片在线播放| 午夜福利,免费看| 男女无遮挡免费网站观看| 国产欧美日韩一区二区三区在线| 精品国产超薄肉色丝袜足j| 欧美成人午夜免费资源| 亚洲少妇的诱惑av| 中文字幕亚洲精品专区| 国产精品蜜桃在线观看| 丝袜在线中文字幕| 国产精品久久久久久久久免| 在线亚洲精品国产二区图片欧美| 香蕉精品网在线| 国产精品久久久久久精品古装| 久久久久精品性色| av不卡在线播放| 春色校园在线视频观看| 午夜福利影视在线免费观看| 2021少妇久久久久久久久久久| 男人爽女人下面视频在线观看| 一级黄片播放器| 香蕉精品网在线| videos熟女内射| 免费黄频网站在线观看国产| 欧美变态另类bdsm刘玥| 在线天堂中文资源库| 色94色欧美一区二区| 免费观看av网站的网址| 黑丝袜美女国产一区| 黄频高清免费视频| 国产免费福利视频在线观看| 欧美xxⅹ黑人| 亚洲一区二区三区欧美精品| 国产精品国产三级专区第一集| 日韩人妻精品一区2区三区| 视频在线观看一区二区三区| 亚洲熟女精品中文字幕| 久久精品久久久久久噜噜老黄| 日韩伦理黄色片| 国产熟女午夜一区二区三区| 久久免费观看电影| 天天操日日干夜夜撸| 天天操日日干夜夜撸| 久久久久人妻精品一区果冻| 最近最新中文字幕大全免费视频 | 亚洲色图综合在线观看| 日本av手机在线免费观看| 啦啦啦视频在线资源免费观看| 色视频在线一区二区三区| 天美传媒精品一区二区| 亚洲国产最新在线播放| 欧美日韩国产mv在线观看视频| 老司机影院成人| 丰满饥渴人妻一区二区三| 美女视频免费永久观看网站| 久久久久久人人人人人| 黑人巨大精品欧美一区二区蜜桃| 蜜桃国产av成人99| 老女人水多毛片| 国产在线一区二区三区精| 午夜激情av网站| 麻豆乱淫一区二区| 国产爽快片一区二区三区| 18禁动态无遮挡网站| 婷婷精品国产亚洲av在线| 婷婷精品国产亚洲av在线| 亚洲激情在线av| 日本精品一区二区三区蜜桃| 亚洲性夜色夜夜综合| 国产成人精品无人区| 亚洲黑人精品在线| 手机成人av网站| 免费高清在线观看日韩| 日韩欧美免费精品| 桃红色精品国产亚洲av| 中文字幕人妻熟女乱码| 久久久久久亚洲精品国产蜜桃av| 久久草成人影院| 人人妻人人添人人爽欧美一区卜| 中文字幕av电影在线播放| 国产亚洲精品一区二区www| 亚洲精品一卡2卡三卡4卡5卡| 悠悠久久av| 国产又爽黄色视频| 免费搜索国产男女视频| 午夜免费观看网址| 免费在线观看黄色视频的| 午夜福利影视在线免费观看| 日韩精品免费视频一区二区三区| 美女 人体艺术 gogo| 国产精品一区二区精品视频观看| 日韩视频一区二区在线观看| 黄色怎么调成土黄色| 国产精品久久久av美女十八| 人人妻人人澡人人看| 九色亚洲精品在线播放| 欧美日韩乱码在线| 99国产精品99久久久久| 国产av又大| av网站在线播放免费| 久久人妻熟女aⅴ| 18禁美女被吸乳视频| tocl精华| 欧洲精品卡2卡3卡4卡5卡区| 久久午夜综合久久蜜桃| 黄色丝袜av网址大全| 人人澡人人妻人| 亚洲自偷自拍图片 自拍| 麻豆久久精品国产亚洲av | 三级毛片av免费| 免费久久久久久久精品成人欧美视频| 午夜免费鲁丝| 久热爱精品视频在线9| 在线视频色国产色| 日韩三级视频一区二区三区| 99在线视频只有这里精品首页| 亚洲国产精品合色在线| 99精国产麻豆久久婷婷| 久久国产亚洲av麻豆专区| 一区二区日韩欧美中文字幕| 亚洲第一欧美日韩一区二区三区| 日本五十路高清| 久久九九热精品免费| 老司机亚洲免费影院| 美女高潮到喷水免费观看| 丰满人妻熟妇乱又伦精品不卡| 久久欧美精品欧美久久欧美| 欧美色视频一区免费| 97人妻天天添夜夜摸| 精品久久久久久久毛片微露脸| xxx96com| 欧美午夜高清在线| 亚洲国产中文字幕在线视频| 在线观看一区二区三区激情| 亚洲成a人片在线一区二区| 热re99久久精品国产66热6| 一夜夜www| 窝窝影院91人妻| 久久精品国产清高在天天线| 日韩欧美在线二视频| 国产三级在线视频| 夜夜夜夜夜久久久久| 12—13女人毛片做爰片一| 老鸭窝网址在线观看| 一本综合久久免费| 国产有黄有色有爽视频| 99久久精品国产亚洲精品| 久久婷婷成人综合色麻豆| 日韩人妻精品一区2区三区| 日韩免费高清中文字幕av| 欧美激情 高清一区二区三区| 亚洲男人的天堂狠狠| 午夜影院日韩av| 日韩有码中文字幕| 亚洲国产欧美日韩在线播放| 久久精品亚洲av国产电影网| 又紧又爽又黄一区二区| 国产黄色免费在线视频| 欧美日韩精品网址| а√天堂www在线а√下载| 亚洲av片天天在线观看| 色综合站精品国产| 亚洲精品久久午夜乱码| a级片在线免费高清观看视频| 老司机福利观看| 免费看a级黄色片| 在线观看日韩欧美| 不卡一级毛片| 91国产中文字幕| 97人妻天天添夜夜摸| 免费女性裸体啪啪无遮挡网站| 91精品国产国语对白视频| 99精品欧美一区二区三区四区| 色尼玛亚洲综合影院| 亚洲国产精品sss在线观看 | 欧美亚洲日本最大视频资源| 中文欧美无线码| 97人妻天天添夜夜摸| 精品免费久久久久久久清纯| 中文字幕精品免费在线观看视频| 日韩欧美免费精品| 午夜日韩欧美国产| 一个人免费在线观看的高清视频| 精品高清国产在线一区| 琪琪午夜伦伦电影理论片6080| 免费在线观看日本一区| 国产精品二区激情视频| 夜夜躁狠狠躁天天躁| 国产91精品成人一区二区三区| 亚洲国产中文字幕在线视频| av欧美777| 亚洲第一欧美日韩一区二区三区| 国产在线观看jvid| 精品欧美一区二区三区在线| 啦啦啦免费观看视频1| 国产97色在线日韩免费| 欧美av亚洲av综合av国产av| 日本vs欧美在线观看视频| 一进一出抽搐gif免费好疼 | 91成人精品电影| 精品无人区乱码1区二区| 日本免费一区二区三区高清不卡 | 久久久久久久久免费视频了| 欧美黑人欧美精品刺激| 高潮久久久久久久久久久不卡| 欧美激情 高清一区二区三区| 男女高潮啪啪啪动态图| 80岁老熟妇乱子伦牲交| 中文字幕最新亚洲高清| 亚洲精品一二三| 国产亚洲欧美精品永久| 国产精品香港三级国产av潘金莲| 亚洲国产中文字幕在线视频| 国产精品国产av在线观看| √禁漫天堂资源中文www| 美女午夜性视频免费| 中出人妻视频一区二区| 久久狼人影院| 欧美另类亚洲清纯唯美| 在线观看一区二区三区| 色尼玛亚洲综合影院| 久久亚洲精品不卡| 久久久久精品国产欧美久久久| 他把我摸到了高潮在线观看| 亚洲欧洲精品一区二区精品久久久| 精品一品国产午夜福利视频| 男女高潮啪啪啪动态图| 日韩欧美一区二区三区在线观看| 久久午夜综合久久蜜桃| a级片在线免费高清观看视频| 黑人猛操日本美女一级片| 在线观看一区二区三区| 免费久久久久久久精品成人欧美视频| 日韩精品免费视频一区二区三区| svipshipincom国产片| 久久精品国产亚洲av香蕉五月| 欧美中文综合在线视频| 久久精品国产清高在天天线| 成人影院久久| 99国产精品免费福利视频| 国产高清国产精品国产三级| 99精品欧美一区二区三区四区| 黄色a级毛片大全视频| 天堂动漫精品| 色尼玛亚洲综合影院| 国产免费av片在线观看野外av| 老司机福利观看| 国产精品香港三级国产av潘金莲| 在线播放国产精品三级| 男女床上黄色一级片免费看| 91精品国产国语对白视频| 久久久精品欧美日韩精品| 精品久久久精品久久久| 久久人妻熟女aⅴ| 女性被躁到高潮视频| 久久亚洲真实| 午夜两性在线视频| 一级作爱视频免费观看| 国产av在哪里看| 国产单亲对白刺激| 国产主播在线观看一区二区| 久久人妻av系列| www.999成人在线观看| 日日爽夜夜爽网站| 99国产综合亚洲精品| 国产免费男女视频| 亚洲国产欧美一区二区综合| 在线观看www视频免费| 久久香蕉激情| 国产野战对白在线观看| 成人三级做爰电影| 久久天堂一区二区三区四区| 久久青草综合色| 欧美激情极品国产一区二区三区| 91九色精品人成在线观看| 丝袜美腿诱惑在线| 欧美不卡视频在线免费观看 | ponron亚洲| 一进一出抽搐动态| 日韩欧美国产一区二区入口| 欧美成人免费av一区二区三区| 亚洲国产毛片av蜜桃av| 午夜久久久在线观看| 亚洲精品一二三| 免费观看精品视频网站| 成人亚洲精品一区在线观看| 亚洲欧美一区二区三区久久| 精品久久久久久电影网| 国产精品影院久久| 一级片'在线观看视频| 精品人妻在线不人妻| 我的亚洲天堂| 亚洲人成77777在线视频| 亚洲国产欧美网| 久久中文字幕人妻熟女| 欧美国产精品va在线观看不卡| 亚洲一码二码三码区别大吗| 又紧又爽又黄一区二区| 中文字幕人妻熟女乱码| 久久国产乱子伦精品免费另类| 成年女人毛片免费观看观看9| 一区二区日韩欧美中文字幕| 50天的宝宝边吃奶边哭怎么回事| 久久久精品欧美日韩精品| 久久精品亚洲熟妇少妇任你| 18禁裸乳无遮挡免费网站照片 | 亚洲 欧美一区二区三区| 午夜福利在线免费观看网站| 一边摸一边抽搐一进一出视频| 在线观看午夜福利视频| 国产欧美日韩一区二区三区在线| 国产单亲对白刺激| 亚洲激情在线av| 亚洲色图综合在线观看| 国产亚洲精品久久久久5区| 水蜜桃什么品种好| 麻豆国产av国片精品| 国产深夜福利视频在线观看| 久久久国产一区二区| 黄片小视频在线播放| 久久 成人 亚洲| 性色av乱码一区二区三区2| 久久中文字幕人妻熟女| 国产伦一二天堂av在线观看| 午夜福利在线观看吧| 黄色 视频免费看| 操美女的视频在线观看| 国产av一区二区精品久久| 午夜免费鲁丝| 黑丝袜美女国产一区| 女性被躁到高潮视频| 桃红色精品国产亚洲av| 亚洲伊人色综图| 91麻豆av在线| 日本 av在线| 黄色a级毛片大全视频| 嫩草影院精品99| 国产精品免费视频内射| 国产成人一区二区三区免费视频网站| 久久亚洲精品不卡| 久久 成人 亚洲| 大型av网站在线播放| 国产99白浆流出| 欧美丝袜亚洲另类 | 精品国产乱码久久久久久男人| 亚洲一区中文字幕在线| 制服诱惑二区| 波多野结衣av一区二区av| 亚洲精品国产色婷婷电影| 日本黄色日本黄色录像| 99国产精品一区二区三区| 十八禁人妻一区二区| www.自偷自拍.com| 人妻久久中文字幕网| 深夜精品福利| 亚洲专区字幕在线| 嫁个100分男人电影在线观看| 三上悠亚av全集在线观看| 久久久久国产一级毛片高清牌| 久久中文看片网| 在线观看舔阴道视频| 久久天躁狠狠躁夜夜2o2o| 午夜精品国产一区二区电影| 婷婷丁香在线五月| 91精品三级在线观看| 欧美乱色亚洲激情| 超色免费av| av中文乱码字幕在线| 黄色丝袜av网址大全| 久久久国产精品麻豆| 又黄又爽又免费观看的视频| 69精品国产乱码久久久| 美女午夜性视频免费| 日韩大尺度精品在线看网址 | 在线av久久热| 午夜免费观看网址| 亚洲精品av麻豆狂野| 最近最新中文字幕大全免费视频| 少妇粗大呻吟视频| 国产精品久久电影中文字幕| 亚洲熟妇中文字幕五十中出 | 一级a爱片免费观看的视频| 黄色女人牲交| 成人国语在线视频| 日本免费一区二区三区高清不卡 | 不卡av一区二区三区| 美女高潮到喷水免费观看| 国产成人av激情在线播放| 一区二区日韩欧美中文字幕| 国产97色在线日韩免费| 亚洲专区中文字幕在线| 久久久国产成人免费| 精品一区二区三区视频在线观看免费 | 国产一卡二卡三卡精品| 狂野欧美激情性xxxx| 神马国产精品三级电影在线观看 | 免费看十八禁软件| 制服人妻中文乱码| 成人三级黄色视频| 国产精华一区二区三区| www日本在线高清视频| 国产精品av久久久久免费| 国产成人精品久久二区二区免费| 免费高清视频大片| 精品一区二区三区av网在线观看| 亚洲激情在线av| 在线播放国产精品三级| 日韩 欧美 亚洲 中文字幕| 成人国语在线视频| 亚洲精品在线观看二区| 99热国产这里只有精品6| 亚洲av五月六月丁香网| 亚洲国产精品一区二区三区在线| 韩国av一区二区三区四区| 日韩欧美一区二区三区在线观看| 午夜久久久在线观看| 日韩精品免费视频一区二区三区| 狂野欧美激情性xxxx| 久久精品91蜜桃| 这个男人来自地球电影免费观看| 久久影院123| 日韩欧美国产一区二区入口| 久99久视频精品免费| 搡老乐熟女国产| 国产一卡二卡三卡精品| 法律面前人人平等表现在哪些方面| 国产一区二区三区视频了| 亚洲人成电影观看| 激情在线观看视频在线高清| 成人免费观看视频高清| 人妻丰满熟妇av一区二区三区| 天堂俺去俺来也www色官网| 国产精品久久视频播放| 欧美日韩av久久| 久久精品影院6| 欧美最黄视频在线播放免费 | 男女做爰动态图高潮gif福利片 | 国产又色又爽无遮挡免费看| 久久久国产成人免费| 亚洲专区中文字幕在线| 一级作爱视频免费观看| 色婷婷久久久亚洲欧美| 在线看a的网站| √禁漫天堂资源中文www| 日韩欧美一区二区三区在线观看| 夫妻午夜视频| 成人国产一区最新在线观看| 免费在线观看黄色视频的| 每晚都被弄得嗷嗷叫到高潮|