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

    基于跨尺度等效彈性的航空裝備損修性能評估

    2025-03-20 00:00:00馮蘊雯宋祉岑路成陳曉宇
    關(guān)鍵詞:系統(tǒng)仿真復(fù)合材料

    摘 要:為解決裝備系統(tǒng)在復(fù)雜損修模式下的性能評估問題,提出一種跨尺度等效彈性仿真方法,在保證計算精度的前提下,通過跨尺度映射減少損修及其幾何特征帶來的建模負(fù)擔(dān)。首先,基于等效彈性理論,在微尺度結(jié)構(gòu)內(nèi),用彈性參數(shù)的等效模擬和取代修理因素對本構(gòu)行為帶來的影響;其次,開展跨尺度建模,將微尺度結(jié)構(gòu)的力學(xué)行為映射到宏觀裝備系統(tǒng)中;最后,以浸水和穿孔修復(fù)后的航空雷達(dá)罩裝備為案例驗證對象進(jìn)行驗證實驗。驗證實驗的結(jié)果表明,分析結(jié)果精度均值為6.3%,且極大地減小了試驗和仿真工作量,為復(fù)雜航空裝備系統(tǒng)的損修性能評估提供思路。

    關(guān)鍵詞: 航空裝備; 系統(tǒng)仿真; 等效彈性; 復(fù)合材料

    中圖分類號: V 229

    文獻(xiàn)標(biāo)志碼: ADOI:10.12305/j.issn.1001 506X.2025.02.17

    Cross scale equivalent elasticity based damage and repair performance

    assessment of aviation equipment

    FENG Yunwen1, SONG Zhicen1,*, LU Cheng1, CHEN Xiaoyu2

    (1. School of Aeronautics, Northwestern Polytechnical University, Xi’an 710072, China; 2. Research Institute for

    Special Structures of Aeronautical Composites, Aviation Industry Corporation of China, LTD., Jinan 250031, China)

    Abstract:To solve the performance evaluation of equipment systems in complex damage and repair mode, a cross scale equivalent elasticity simulation method is proposed, which ensures the calculation accuracy and reduces the modeling burden caused by damage and repair and geometric features by cross scale mapping. Firstly, the equivalent elasticity theory is introduced, and the equivalent elastic parameters are used to simulate and replace the influence of repair factors on the constitutive behavior in the micro scale structure. Secondly, cross scale modeling is carried out to map the mechanical behavior of the micro scale structure to the macroscopic equipment system. Finally, experiment is carried out with the verification case of aviation radar radome after water immersion or perforation repairing and the result of verification experiment shows that the average accuracy of the analysis results is 6.3%. It also greatly reduces the workload of testing and simulation and provides an idea for the damage and repair performance evaluation of complex aviation equipment systems.

    Keywords:aviation equipment; system simulation; equivalent elasticity; composite material

    0 引 言

    要保障航空器在其壽命周期內(nèi)持續(xù)安全、可靠,必須對其結(jié)構(gòu)裝備,如雷達(dá)罩、機身、機翼等,開展修復(fù)后性能和安全性評估。然而,飛機運行過程中的損傷形式復(fù)雜多樣、損修后結(jié)構(gòu)強度測試成本高昂,這些因素都給上述工作帶來阻礙。以雷達(dá)罩裝備為例,其由纖維增強復(fù)合材料夾芯結(jié)構(gòu)制成,作為飛機最前端的裝備,其損傷類型復(fù)雜多樣。使用殼單元建模對雷達(dá)罩這類薄壁裝備的工程計算與理論研究帶來極大的便利,但同時也導(dǎo)致在分析中需要考慮損傷或修復(fù)因素時,對浸水、分層、脫粘或凹坑等復(fù)雜因素在宏觀裝備系統(tǒng)中很難模擬出毫米級別的真實損修模式。對民用飛機而言,開展裝備或結(jié)構(gòu)的修后性能分析是維修保障的重要組成部分。

    針對這類復(fù)雜裝備,尤其是具有曲率的復(fù)合材料蜂窩夾芯裝備或結(jié)構(gòu),判斷其是否在維修后保持和恢復(fù)良好的技術(shù)狀態(tài)已成為不少相關(guān)學(xué)者的研究重點。對復(fù)合材料夾芯結(jié)構(gòu)而言,在拉壓1-2、面載3、彎曲4、沖擊5等載荷下,其修復(fù)后性能依賴于修補方式、黏合界面、斜挖角度等修理參數(shù);連接處也是裝備的薄弱環(huán)節(jié),連接接頭6-7或筋條、緣條等8通過選擇不同的修理界面參數(shù)可以改善應(yīng)力集中等現(xiàn)象。當(dāng)裝備的幾何外形帶有曲率,甚至是不規(guī)則流線形,如機身壁板9、雷達(dá)罩10、機翼壁板11、發(fā)動機整流罩12、渦輪葉盤13等,此時平面結(jié)構(gòu)已經(jīng)不能滿足實際裝備的性能分析需求。針對這類問題,許多學(xué)者開展修補曲面建模與分析,如邱求元等14分析單曲率層壓殼結(jié)構(gòu)穿孔修復(fù)性能表現(xiàn),在曲面結(jié)構(gòu)上,楔形補片相比于膠補能夠弱化應(yīng)力集中;Citil等15評估具有曲面搭接接頭的鋼管對黏合界面的敏感性,曲率對破壞壓力有著較大的影響;Urbanek等16通過定向能量沉積對曲面薄壁結(jié)構(gòu)開展修理,并預(yù)測該工藝制造的薄壁結(jié)構(gòu)在靜態(tài)載荷下的行為。上述研究對象通常選擇修理特征位于中心的小尺寸結(jié)構(gòu)件,以便于分析修理參數(shù)對性能的影響,然而在真實裝備中,損傷或修復(fù)的位置具有不確定性。近年來獲得發(fā)展的尺度建模方法,既具備層板結(jié)構(gòu)的宏觀尺度,又具備組分材料的細(xì)觀尺度,以及孔隙微裂紋等微觀尺度,許多研究同時兼顧微觀模型、細(xì)觀模型,甚至宏觀模型,將材料性能、力學(xué)行為在不同尺度間映射、轉(zhuǎn)換。Aboudi等17提出的單元法,通過將細(xì)觀單胞模型離散為子胞,結(jié)合子胞間的連續(xù)條件,求解細(xì)觀應(yīng)力應(yīng)變場;Babuska等18提出漸進(jìn)均勻化理論,對宏觀應(yīng)力和位移場進(jìn)行漸進(jìn)展開,進(jìn)行宏觀等效參數(shù)的求解;Fu等19利用多尺度建模描述短碳纖維聚合物材料空隙缺陷,從制造、生成到演化的全過程,考慮微觀尺度的數(shù)值結(jié)果,性能偏差值更小;Wineman等20在各向異性等效彈性理論內(nèi),總結(jié)多力學(xué)響應(yīng),如非線性粘彈性、彈熱性的本構(gòu)模型,展示纖維或基體性能如何在數(shù)學(xué)模型中與整體性能相關(guān)聯(lián)。尺度策略的發(fā)展為這種具有損傷或修復(fù)模式的復(fù)合材料夾芯曲面裝備或結(jié)構(gòu)的性能分析提供思路,但上述研究也存在如下兩方面的問題:一方面,將仿真或試驗驗證停留在結(jié)構(gòu)級,曲率結(jié)構(gòu)也多是單曲率或規(guī)則圓柱體,與實際裝備仍存在差距;另一方面,微觀與宏觀的尺度轉(zhuǎn)換多考慮微觀材料間的性能差異,難以映射損傷、修復(fù)等因素給結(jié)構(gòu)特征帶來的影響。

    針對上述研究存在的問題,本文將等效彈性理論引入到尺度建模中,提出一種跨尺度等效彈性方法,通過等效彈性映射修復(fù)因素及其幾何特征對本構(gòu)帶來的影響,結(jié)合跨尺度建模方法,將修復(fù)因素的力學(xué)特性在細(xì)觀與宏觀尺度間傳遞,解決損傷或修復(fù)模式的復(fù)雜曲率裝備性能評估問題,為裝備維修后的性能分析和保障提供理論支撐。以雷達(dá)罩裝備為主要應(yīng)用對象,對方法和案例進(jìn)行詳盡的說明。相較于傳統(tǒng)方法,通過跨尺度等效彈性方法可以在保證計算精度的前提下,極大地減少建模負(fù)擔(dān)和降低計算規(guī)模。

    1 跨尺度等效彈性方法

    1.1 等效彈性理論

    1.1.1 等效彈性模量

    彈性參數(shù)來自于細(xì)觀力學(xué)理論,彈性參數(shù)是人為定義的、材料的二維或三維力學(xué)行為的表征。將由多種材料制成的新材料(如纖維增強復(fù)合材料單向板、合金21、巖石22-23、土介質(zhì)24等)或多維結(jié)構(gòu)(如蜂窩芯子25、層合板26、接頭27、多相復(fù)合材料28等)視作一個整體,用一組參數(shù)來表征其整體力學(xué)行為。等效彈性模量的方法是一種類似于代理模型的數(shù)值解析方法,在一定尺度內(nèi),將具有某種特征及其附近的區(qū)域視作一個整體,用力學(xué)性能參數(shù)的修正等效內(nèi)部或外部結(jié)構(gòu)帶來的影響。在線彈性范圍內(nèi),以拉伸載荷引起的應(yīng)力、應(yīng)變?yōu)槔?,?yīng)力和應(yīng)變之間的關(guān)系如下所示:

    PA=σ=E·ε=E·ΔlL

    Δl=LE·PA(1)

    式中:P為結(jié)構(gòu)受載;σ為應(yīng)力;ε為正應(yīng)變;E為模量;A為橫截面積;L為結(jié)構(gòu)長度;Δl為結(jié)構(gòu)在相應(yīng)載荷下的位移。

    根據(jù)等效彈性理論,損傷或修復(fù)的結(jié)構(gòu)細(xì)節(jié)能夠通過本構(gòu)關(guān)系的改變映射在一個較小的區(qū)域內(nèi),該尺度被稱為一個等效微結(jié)構(gòu)。相對于宏觀結(jié)構(gòu),等效彈性模量Eeq或等效剪切模量Geq能夠在該區(qū)域內(nèi)表達(dá)出相近的本構(gòu)關(guān)系。令等效前后在區(qū)域邊界引起的邊界條件一致,結(jié)構(gòu)與等效結(jié)構(gòu)的力學(xué)表征如下所示:

    PA=P′A′=σ-=Eeq·ε-=Eeq·ΔltotalL

    PA=P′A′=τ-=Geq·γ-=Geq·ΔltotalH

    LE·PA=LEeq·P′A′

    (2)

    式中:P′為修理后微結(jié)構(gòu)邊界受載;σ-為等效應(yīng)力;ε-為等效正應(yīng)變;τ-為等效切應(yīng)變;γ-為等效切應(yīng)力;Δltotal為總位移;A′為修理后微結(jié)構(gòu)橫截面大小,H為結(jié)構(gòu)厚度。上式不僅適用于線彈性階段,對屈服階段也同樣適用,只是在驅(qū)服階段模量會發(fā)生退化。在拉伸或剪切載荷下,相應(yīng)位移的取法如圖1所示。圖1(a)為三維視圖,圖1(b)為俯視圖。

    均勻化和等效化后,等效彈性模量是與材料、結(jié)構(gòu)特征等相關(guān)的參數(shù):

    E-=f(Ea,va12,Eb,vb12,k,g,…)

    G-=f(Gaij,va12,Gbij,vb12,k,g,…)(3)

    式中:E或G為等效微結(jié)構(gòu)內(nèi)組成材料本身的模量,如蜂窩夾芯結(jié)構(gòu)由面板和芯子共同構(gòu)成;Ea為面板彈性模量;va12為面板泊松比;Gaij為面板剪切模量;Eb為芯子彈性模量;vb12為芯子泊松比;Gbij為芯子剪切模量;k為結(jié)構(gòu)細(xì)節(jié)參數(shù),g為與尺度相關(guān)的度量;f(·)為輸入與輸出間的函數(shù)關(guān)系。

    1.1.2 等效彈性模量仿真計算方法

    等效彈性模量仿真的思路是在微結(jié)構(gòu)尺度內(nèi)盡可能詳細(xì)地構(gòu)建和描述損傷和修復(fù)所能包含的幾何特征細(xì)節(jié)。在該微結(jié)構(gòu)內(nèi),由于包含損傷和修復(fù)信息,其模量、泊松比等彈性參數(shù)必然發(fā)生變化。在無損模型上引入等效模量,則能夠減少仿真計算的規(guī)模,提高仿真效率。求解的等效彈性參數(shù)集包括S={E1,E2,v12,G12,G13,G23}。其中:E1為0°拉伸模量,E2為90°拉伸模量,v12為泊松比;G12為面內(nèi)剪切模量;G13和G23為面外剪切模量。

    對周期性結(jié)構(gòu)的復(fù)合材料彈性體Ω,在給定位移邊界條件時,其宏觀結(jié)構(gòu)可看作是微觀單元在空間中周期性的重復(fù)排列,這種將力學(xué)行為在不同尺度間映射的漸進(jìn)式分析方法被稱為多尺度建模(multi scale modeling, MSM)[29-30。單胞的尺度相對于宏觀幾何尺度為小量,如圖2所示。常用的宏觀尺度單位為m,微觀尺度為nm。單胞模型在nm級別內(nèi)很難將修復(fù)模式全面地涵蓋在內(nèi)。因此,本文將修復(fù)結(jié)構(gòu)細(xì)節(jié)映射在一個較小的mm級別區(qū)域內(nèi),將該尺度稱為一個等效彈性的微結(jié)構(gòu)。微結(jié)構(gòu)的尺寸根據(jù)修理的尺寸來設(shè)定,原則上等效彈性模型邊緣應(yīng)不再引起應(yīng)力集中,且等效彈性模型應(yīng)遠(yuǎn)小于裝備。

    穿孔模型的應(yīng)力集中最為明顯,以一個帶有直徑D為25 mm的圓孔的微結(jié)構(gòu)為例,當(dāng)微結(jié)構(gòu)的外部尺寸為修補直徑的1.5倍時,其在拉伸載荷下的應(yīng)力云圖如圖3所示。過多的無損單元會消磨掉帶有特別結(jié)構(gòu)細(xì)節(jié)的單元的模量改變量,且應(yīng)當(dāng)覆蓋本構(gòu)因獨有的幾何特征而發(fā)生改變的區(qū)域范圍。結(jié)合圖3(a)的應(yīng)力云圖也可以觀察到,在寬度W=1.5D、長度L=1.5D的范圍內(nèi),應(yīng)力、應(yīng)變、位移的改變量相對明顯,且又包括應(yīng)力集中的區(qū)域。因此,第一種尺度模型為外部尺寸為1.5D mm×1.5D mm的平板等效微結(jié)構(gòu)模型。

    損傷或修復(fù)蜂窩夾芯結(jié)構(gòu)等效彈性模量求解在拉伸和剪切兩種載荷下開展,各類載荷的加載條件如圖4所示,其中U1為1方向位移;U2為2方向位移;U3為3方向位移,取有限元分析的結(jié)果計算等效彈性模量:

    E1=RF1/(L·H)ΔU1/L

    E2=RF2/(L·H)ΔU2/L

    v12=ΔU1ΔU2 G12=RF1/(L·H)ΔU1/L

    G13=RF1/(L·L)ΔU1/H

    G23=RF2/(L·L)ΔU2/H (4)

    式中:RF1為1方向作用力;ΔU1為1方向位移變化量;ΔU2為2方向位移變化量。

    1.2 跨尺度建模仿真

    雷達(dá)天線罩的所用芯材與面板與其他裝備相比都相對較薄,計算中通常使用殼元。這樣的處理一方面可以減少計算量,另一方面降低分析不收斂的風(fēng)險,但損傷及損傷修復(fù)的參數(shù)和尺寸卻難以體現(xiàn)在殼元的建模和計算中。穿孔和穿孔修補可以通過簡單的靜態(tài)建模來解決,然而對復(fù)雜修復(fù)特征,如分層和脫粘,在只有單層的殼元模型內(nèi)無法區(qū)分夾芯結(jié)構(gòu)不同區(qū)域的分離現(xiàn)象。如浸水修復(fù),在只有單層的殼元模型內(nèi)無法在上蒙皮或下蒙皮均勻打孔并保留芯材的完整性。另外,小孔周圍的區(qū)域在劃分網(wǎng)格時易發(fā)生高度扭曲。因此,在采用跨尺度等效彈性方法時,要進(jìn)行篩選和鑒別,如穿孔和雙層蒙皮夾芯材料損傷修復(fù)這類可以通過簡單的靜態(tài)建模來解決的問題,不必使用等效方法,同時能夠避免均勻化帶來的計算誤差;如分層和脫粘、浸水修復(fù)這類無法在殼元上通過靜態(tài)建模來解決的問題,則應(yīng)當(dāng)基于多尺度的等效彈性方法首先計算出等效彈性參模量,將其代入宏觀模型中進(jìn)行分析和計算。在宏觀模型中,應(yīng)當(dāng)切分出大小盡可能近似1.5D mm×1.5D mm的區(qū)域。為盡可能使曲率、鋪層角度規(guī)則化,則左右應(yīng)當(dāng)平行于雷達(dá)罩母線,上下應(yīng)當(dāng)平行于雷達(dá)罩底面。切分方法如圖5所示。在平行于雷達(dá)罩底面的平面上以修復(fù)區(qū)域中心點坐標(biāo)(x,z)為中心點,畫出距離中心點為1.5D mm的4個點,連接(0,0)和左右兩點,這兩條線在雷達(dá)罩上的投影與其母線應(yīng)當(dāng)重合;過上下兩點做兩條平行線,且兩條線均垂直于(0,0)和修復(fù)區(qū)域中心點的連線。以(x,z)為中心,將平面旋轉(zhuǎn)到平行于曲面的切面方向,由這4條邊所圍住的區(qū)域向雷達(dá)罩投影所形成的區(qū)域,即為損傷等效微結(jié)構(gòu)所在的區(qū)域。因此,第二種尺度模型為經(jīng)過切分的宏觀模型,如圖5所示。

    等效模型以外的應(yīng)變測量點,其計算結(jié)果應(yīng)當(dāng)滿足精度要求。對等效模型上的應(yīng)變測量點,應(yīng)當(dāng)結(jié)合等效模型的載荷和約束條件,進(jìn)行復(fù)合材料蜂窩夾芯結(jié)構(gòu)曲板件的細(xì)觀建模,分析因均勻化而被模糊的真實的應(yīng)變場和應(yīng)力場。因此,第3種尺度模型為帶有曲率和修復(fù)模式細(xì)節(jié)建模的曲板細(xì)觀模型。總體而言,跨尺度損傷修復(fù)等效建模和分析的過程包括:①損修平板精細(xì)化模型將等效模量代入整體裝備;②從整體裝備模型中獲取損修局部區(qū)域的邊界條件;③將邊界條件賦給具有真實曲率和損修細(xì)節(jié)的曲板損修精細(xì)化模型;④對曲板損修精細(xì)化模型開展仿真分析,獲得損修區(qū)域真實的應(yīng)力、應(yīng)變、分布狀態(tài)信息。多尺度跨越建模過程如圖6所示。

    1.3 基于跨尺度等效彈性的性能評估框架

    跨尺度等效彈性方法考慮復(fù)雜修復(fù)模式下的力學(xué)性能改變量,用等效彈性模量來模擬一定尺寸內(nèi)具有復(fù)雜修復(fù)模式結(jié)構(gòu)特征微結(jié)構(gòu)的本構(gòu)行為,再結(jié)合多尺度策略,將微結(jié)構(gòu)的力學(xué)行為投射到宏觀模型中,并結(jié)合細(xì)觀模型判斷真實應(yīng)力場變化。等效模量能夠在整體模型中,在等效區(qū)域內(nèi)模擬出具有復(fù)雜模式的性能狀態(tài)。因此,基于等效彈性理論和多尺度策略的力學(xué)性能分析框架如圖7所示。

    2 雷達(dá)罩裝備損修性能評估

    2.1 損傷修復(fù)雷達(dá)罩的物理信息

    飛機雷達(dá)罩的作用是保護(hù)天線,其所處的復(fù)雜特殊的工作環(huán)境,要求其具有高的比強度、高的比剛度和良好的電磁波透過性能,因此雷達(dá)罩在材料使用方面首選復(fù)合材料蜂窩夾芯結(jié)構(gòu)。民機雷達(dá)罩如圖8所示,根據(jù)材料和結(jié)構(gòu)設(shè)置的不同,可將罩體分為4個部分。其中:①和②區(qū)域為蜂窩夾芯結(jié)構(gòu),蒙皮為3層單層厚度為0.267 mm的7781織物,芯子是兩種芳綸紙六邊形蜂窩芯子;③為復(fù)合材料層合結(jié)構(gòu),為11層單層厚度為0.267 mm的7781織物;④為金屬鋁結(jié)構(gòu)。本文采用雷達(dá)罩整體坐標(biāo)系。雷達(dá)罩頂點坐標(biāo)為(0,0,0),坐標(biāo)系的方向如圖1所示。

    該雷達(dá)罩產(chǎn)生損傷共計7處,并進(jìn)行修理,其中雙層蒙皮加芯子修補4處,浸水修補3處。穿孔的損傷為直徑為100 mm的圓形孔,浸水區(qū)域為直徑為100 mm的圓形區(qū)域。將試驗件通過鉸鏈、連接鎖及5個錐銷安裝于支持夾具上,將支持夾具固定于承力地坪上,如圖8所示。損傷修理在雷達(dá)罩中的相對位置如圖8所示,坐標(biāo)信息被展示在表1中。

    試驗件采用靜力試驗加載木塊固定結(jié)構(gòu),通過作動筒和木塊將集中力轉(zhuǎn)換為面載,集中力的大小根據(jù)實際氣動載荷轉(zhuǎn)換而來,如圖9所示。

    2.2 浸水損修雷達(dá)罩跨尺度建模

    浸水的尺度模型包含3個維度,分別是用來求解等效彈性參數(shù)的平板夾芯微結(jié)構(gòu)、細(xì)觀曲板浸水修復(fù)模型和宏觀蜂窩夾芯雷達(dá)罩。雷達(dá)天線罩的基本材料屬性如表2所示,穿孔修補時蒙皮采用原材料,芯子所用材料不同。為等效區(qū)域內(nèi)的材料屬性賦等效彈性參數(shù),其具體結(jié)果在開展等效彈性模量分析后確定。

    根據(jù)維護(hù)手冊,水浸的修復(fù)方法是根據(jù)水浸周圍的芯材網(wǎng)格孔的位置,在蒙皮上根據(jù)蜂窩尺寸,在相對蜂窩格孔中心的位置開直徑為1.5~1.6 mm的小孔,以排出水,在蒙皮設(shè)置附加的鋪層是為保持結(jié)構(gòu)的承重,其細(xì)節(jié)如圖10所示。按照第1.3節(jié)介紹的多尺度建模方法,在宏觀模型中應(yīng)當(dāng)切分出與微結(jié)構(gòu)尺寸盡可能相似的區(qū)域,對雷達(dá)罩3處浸水修復(fù)和4處穿孔損傷修復(fù)切分。雷達(dá)罩在試驗中共設(shè)置15處應(yīng)變,測量每一處0°、45°和90°的應(yīng)變,并設(shè)置2處位移測量,距離各損傷修復(fù)中心點距離最近的應(yīng)變信息如表3所示。以2處位移和大于2×10-4的應(yīng)變值作為仿真和試驗的比對指標(biāo)。恰有一處應(yīng)變距離損傷修復(fù)中心位置的距離為55 mm,該應(yīng)變在等效的區(qū)域內(nèi)。在該區(qū)域內(nèi)的力學(xué)行為被均勻化,此時應(yīng)當(dāng)通過細(xì)觀模型來觀察真實的應(yīng)力、應(yīng)變場分布。因此,本文的等效微結(jié)構(gòu)、細(xì)觀曲板模型、宏觀雷達(dá)罩模型如圖11和表3所示。

    2.3 浸水損修等效彈性模量仿真計算

    現(xiàn)得到復(fù)合材料蜂窩夾芯結(jié)構(gòu)平板在浸水修復(fù)下共計2類模型(浸水修復(fù)夾芯結(jié)構(gòu)平板150 mm×150 mm(夾芯為蜂窩芯子1)、浸水修復(fù)夾芯結(jié)構(gòu)平板150 mm×150 mm(夾芯為蜂窩芯子2),針對每類模型得到6類等效彈性參數(shù),共計12類。結(jié)合第1.2節(jié)中的分析方法,浸水修復(fù)在尺寸為150 mm×150 mm以內(nèi)的等效彈性模量如表4所示。

    以7781織物和蜂窩芯子2組成的蜂窩夾芯結(jié)構(gòu)為例,其應(yīng)變云圖如圖12所示。在微結(jié)構(gòu)內(nèi),可以通過無損建模并在源修復(fù)區(qū)域設(shè)置等效彈性模量的方法初步判斷這種等效的彈性模量能否模擬出真實的損傷修復(fù)模型的本構(gòu)關(guān)系或力學(xué)行為。因此,構(gòu)建一個150 mm×150 mm的平板,在等效彈性模量的材料屬性下,其應(yīng)變云圖如圖13所示。二者間的誤差如表5所示。分析結(jié)果較清楚地表明,等效后的模型的應(yīng)變場趨于均勻化,其大小接近原微結(jié)構(gòu)內(nèi)的中位數(shù)或平均值;從整個微結(jié)構(gòu)而言,等效模型能夠很好地表達(dá)E1和E2和泊松比,對剪切行為、尤其是層間剪切行為的表達(dá)能力稍弱,但還在可接受的范圍內(nèi);蜂窩夾芯結(jié)構(gòu)的E1、E2和G12受蒙皮的影響更大,G13和G23受芯子材料屬性的影響更大。

    3 性能評估與討論

    在等效區(qū)域內(nèi)賦等效彈性模量,在穿孔修補區(qū)域內(nèi)賦修補材料屬性,在其他區(qū)域賦雷達(dá)罩裝備的材料屬性。為分析在100%等效的實際氣動載荷下,蜂窩夾芯雷達(dá)罩的力學(xué)行為,結(jié)合上文所述的約束條件設(shè)置、載荷條件設(shè)置、載荷加載分區(qū)和等效區(qū)域劃分,完成有限元模型的輸入設(shè)置,并進(jìn)行分析,損傷修復(fù)雷達(dá)罩在100%試驗載荷下的有限元分析結(jié)果如圖14所示。

    對大于2×10-4的應(yīng)變值和位移進(jìn)行精度的對比,除31號應(yīng)變片,其他精度均小于10%。從損傷修復(fù)雷達(dá)罩的仿真結(jié)果可以看出,在等效區(qū)域內(nèi)外,應(yīng)力、應(yīng)變和位移是連續(xù)變化的,不會因帶入等效彈性參數(shù)而發(fā)生突變,雷達(dá)罩的力學(xué)響應(yīng)主要取決于裝備的約束、載荷條件。在等效區(qū)域內(nèi)使用等效彈性模量就意味著在該區(qū)域內(nèi)設(shè)置均勻化的假設(shè),將損修幾何結(jié)構(gòu)在該區(qū)域內(nèi)擴散開。因此,對由預(yù)測得到的等效區(qū)域內(nèi)的真實應(yīng)力-應(yīng)變場,需要進(jìn)一步結(jié)合細(xì)觀模型構(gòu)建及仿真分析以驗證。

    應(yīng)變點31處的誤差為23.61%,這是因為應(yīng)變在等效彈性模型以內(nèi),此處的應(yīng)變值因均勻化的原因被離散開,應(yīng)當(dāng)通過曲板的細(xì)觀建模和仿真分析,判斷該等效彈性模型內(nèi)的應(yīng)變場變化。曲板模型的4邊約束條件取損傷修復(fù)罩在點位3處的等效彈性模型4邊的位移,等效彈性模型網(wǎng)格尺寸為10 mm,4邊約有16個節(jié)點,取16個節(jié)點Ux、Uy、Uz位移的平均值作為曲板的載荷和約束條件。復(fù)合材料蜂窩夾芯結(jié)構(gòu)曲板損傷修復(fù)有限元分析結(jié)果如圖15所示。對應(yīng)的應(yīng)變測量處的應(yīng)變?yōu)椋害?1=-2.91×10-4、ε22=-4.79×10-5、ε12=7.78×10-4。結(jié)合細(xì)觀建模,該點處的絕對誤差由23.61%降低至2.54%,整體的平均誤差由8.93%降至6.3%。分析結(jié)果和誤差對比分別被展示在圖15和圖16中。

    4 結(jié) 論

    本文提出一種結(jié)合多尺度策略的等效彈性模量方法,為在考慮修理因素的同時減少復(fù)雜曲率裝備的建模和計算負(fù)擔(dān),通過等效彈性模量映射修理特征在不同尺度間的性能變化,針對復(fù)雜曲率裝備,如蜂窩夾芯雷達(dá)罩,在復(fù)雜修理模式下的力學(xué)性能評估問題提供一種思路,形成的主要結(jié)論如下:

    (1) 復(fù)雜的損傷或修復(fù)模式,如浸水、分層、脫粘等,在殼單元中無法通過靜態(tài)建模反映真實狀態(tài)。在模型復(fù)雜、易發(fā)生網(wǎng)格扭曲或尖銳時,通過等效彈性模量在不同尺度間映射力學(xué)行為的方式能夠極大地減少建模的負(fù)擔(dān),且其分析精度在工程應(yīng)用可接受的范圍內(nèi);

    (2) 當(dāng)需要分析等效區(qū)域內(nèi)的應(yīng)變場分布時,需要通過細(xì)觀模型來求解。細(xì)觀模型能夠很好地表達(dá)出該區(qū)域內(nèi)的真實狀態(tài)。一方面,其結(jié)合多尺度策略,可以實現(xiàn)真實應(yīng)變和等效應(yīng)變的轉(zhuǎn)換;另一方面,等效彈性模量能夠在整體模型中發(fā)揮應(yīng)有的本構(gòu)行為,實際和等效后的邊界位移是相近的,可以作為細(xì)觀模型的載荷條件;

    (3) 本文采用的等效彈性模量分析方法在E1、E2和v12的表達(dá)能力較好,但在剪切行為、尤其是層間剪切中描述較弱。若想提高這類等效彈性模量方法的分析精度,可以從模量計算精度方面入手,尤其是對剪切試驗行為的模仿,本文將在今后的研究中關(guān)注這一方向。

    參考文獻(xiàn)

    [1]JEBADURAI D S, KUMAR A A J. Influence of different core configurations on the skin core bonding of sandwich composites[J]. Journal of the Brazilian Society of Mechanical Sciences and Engineering, 2023, 45(8): 391.

    [2]CHEN J W, CHEN W J, WANG M Y, et al. Mechanical behaviors and elastic parameters of laminated fabric URETEK3216LV subjected to uniaxial and biaxial loading[J]. Applied Composite Materials, 2017, 24: 1107-1136.

    [3]ZHANG J, CHENG X Q, ZHANG J K, et al. Effect of curing condition on bonding quality of scarf repaired composite laminates[J]. Chinese Journal of Aeronautics, 2020, 33(8): 163-173.

    [4]SONG S J, XIONG C, YIN J H, et al. Flexural behavior, fai lure analysis, and optimization design of a hybrid composite Kagome honeycomb sandwich structure[J]. Thin Walled Structures, 2023, 187: 110743.

    [5]ZORER O O N. An experimental study on low velocity impact behavior of thermoplastic composites repaired by composite patches[J]. Journal of Composite Materials, 2020, 54(28): 4515-4524.

    [6]喻健, 張騰, 何宇廷, 等. 膠鉚混合修補鋁合金板的疲勞性能研究[J]. 北京航空航天大學(xué)學(xué)報, 2021, 47(11): 2399-2406.

    YU J, ZHANG T, HE Y T, et al. Fatigue performance of aluminum alloy plate with hybrid repair by glue rivet[J]. Journal of Beijing University of Aeronautics and Astronautics, 2021, 47(11): 2399-2406.

    [7]RASLAVIIUS L, BAZARAS I, LUKOEVIIUS V, et al. Statistical investigation of the weld joint efficiencies in the repaired WWER pressure vessel[J]. International Journal of Pressure Vessels and Piping, 2020, 189: 1042711.

    [8]TURKI A Y, AL FARTTOOSI M H. Flexural strength of damaged RC beams repaired with carbon fiber reinforced polymer (CFRP) using different techniques[J]. Fibers, 2023, 11(7): 61.

    [9]AABID A, HRAIRI M, ALI J S M. A review on reductions in the stress intensity factor of cracked plates using bonded composite patches[J]. Materials, 2023, 15(9): 3086.

    [10]王樂, 武岳, 薛鵬, 等. 考慮蒙皮效應(yīng)的雷達(dá)罩結(jié)構(gòu)受力性能試驗[J]. 哈爾濱工業(yè)大學(xué)學(xué)報, 2020, 52(8): 176-183.

    WANG L, WU Y, XUE P, et al. Structural stress perfor mance test of radome considering skin effect[J]. Journal of Harbin Institute of Technology, 2020, 52(8): 176-183.

    [11]HIREMATH P, VISWAMURTHY S R, SHETTAR M. Damage tolerance of a stiffened composite panel with an access cutout under fatigue loading and validation using FEM analysis and digital image correlation[J]. Fibers, 2022, 10(12): 105.

    [12]FEI C W, LIU H T, LI S L, et al. Dynamic parametric modeling based model updating strategy of aeroengine casings[J]. Chinese Journal of Aeronautics, 2021, 34(12): 145-157.

    [13]LU C, TENG D, KESHTEGAR B. Extremum hybrid intelligent inspired models for accurate predicting mechanical performances of turbine blisk[J]. Mechanical Systems and Signal Processing, 2023, 190: 110136.

    [14]邱求元, 范海蓉, 林振能. 損傷單曲層壓殼膠接修復(fù)結(jié)構(gòu)力學(xué)性能測試與分析[J]. 玻璃鋼/復(fù)合材料, 2015(4): 5-10.

    QIU Q Y, FAN H R, LIN Z N. Testing and analyzing the mechanical properties of damaged single curved laminated shell glued repair structures[J]. Fiber Reinforced Plastics/Compo sites, 2015(4): 5-10.

    [15]CITIL S, BOZKURT I, AYDIN M D. Experimental and 3D non linear stress analysis of adhesively bonded pipes with curved surface lap joints[J]. The Journal of Adhesion, 2019, 95(5/7): 515-528.

    [16]URBANEK M, HODEK J, MELZER D, et al. Prediction of behaviour of thin walled DED processed structure: experimental numerical approach[J]. Materials, 2022, 15(3): 806.

    [17]ABOUDI J, GILAT R. The effect of local and random fiber waviness on the microbuckling of composite materials[J]. International Journal of Solids and Structures, 2021, 254: 111862.

    [18]BABUSKA I, MOTAMED M. A fuzzy stochastic multiscale model for fiber composites: a one dimensional study[J]. Computer Methods in Applied Mechanics amp; Engineering, 2022, 15(9): 3086.

    [19]FU Y T, LI J, LI Y Q, et al. Full process multi scale morphological and mechanical analyses of 3D printed short carbon fiber reinforced polyetheretherketone composites[J]. Composites Science and Technology, 2023, 236: 109999.

    [20]WINEMAN A, PENCE T J. Fiber reinforced composites: nonlinear elasticity and beyond[J]. Journal of Engineering Mathematics, 2022, 127(1): 30.

    [21]LEHOCKA D, BOTKO F, KLICH J, et al. Effect of pulsating water jet disintegration on hardness and elasticity modulus of austenitic stainless steel AISI 304L[J]. The International Journal of Advanced Manufacturing Technology, 2020, 107(5/6): 2719-2730.

    [22]XU N Z, LIU C Q, SU C, et al. Study on the size and occurrence effects of equivalent elasticity modulus of orthogonal random jointed rock masses[J]. Frontiers in Earth Science, 2022, 10: 888551.

    [23]PANTELIDIS L. The equivalent modulus of elasticity of soil mediums for designing shallow foundations[J]. Geotechnical and Geological Engineering, 2021, 39: 3863-3873.

    [24]PANTELIDIS L. The equivalent modulus of elasticity of la yered soil mediums for designing shallow foundations with the Winkler spring hypothesis: a critical review[J]. Engineering Structures, 2019, 201: 109452.

    [25]HE D, FENG J Y. An analytical model for predicting equivalent elastic moduli of micro/nano honeycombs with nonlocal effects[J]. Semantic Scholar, 2023, 120: 420-435.

    [26]SORIANO J, FIORELLI J, JUNIOR W E L, et al. Numerical modeling for adjustment of the equivalent moduli of elasticity of OSB layers estimated from experimental flexural rigidity[J]. Journal of Materials Research and Technology, 2021, 14(3):1630-1643.

    [27]MARINELLI F, WEISS B A, BERLI M, et al. Equivalent elastic modulus for prediction of deformations in joints[J]. Ingenius, Revista de Ciencia Tecnologia, 2018, 20: 70-82.

    [28]TIAN L, ZHAO H T, WANG G N, et al. Elasticity based locally exact homogenization theory for three phase composites considering the morphological effect of carbon fibers[J]. Composite Structures, 2023, 304(2): 116428.

    [29]鄭曉霞, 鄭錫濤, 緱林虎. 多尺度方法在復(fù)合材料力學(xué)分析中的研究進(jìn)展[J]. 力學(xué)進(jìn)展, 2010, 40(1): 41-56.

    ZHENG X X, ZHENG X T, GOU L H. Research progress of multiscale methods in mechanical analysis of composite mate rials[J]. Advances in Mechanics, 2010, 40(1): 41-56.

    [30]孫李剛, 凌超, 陳浩, 等. 結(jié)構(gòu)完整性分析中的多尺度力學(xué)方法[J]. 機械工程學(xué)報, 2021, 57(16): 106-121.

    SUN L G, LING C, CHEN H, et al. A multiscale mechanics approach to structural integrity analysis[J]. Journal of Mechanical Engineering, 2021, 57(16): 106-121.

    作者簡介

    馮蘊雯(1968—),女,教授,博士,主要研究方向為可靠性分析、修理工程分析、運行支持。

    宋祉岑(1998—),女,博士研究生,主要研究方向為飛行器設(shè)計、損修性能分析、可靠性分析。

    路 成(1989—),男,博士,主要研究方向為維修工程、損修性能分析。

    陳曉宇(1996—),女,助理工程師,碩士,主要研究方向為結(jié)構(gòu)強度分析、損修性能分析。

    猜你喜歡
    系統(tǒng)仿真復(fù)合材料
    金屬復(fù)合材料在機械制造中的應(yīng)用研究
    纖維素基多孔相變復(fù)合材料研究
    民機復(fù)合材料的適航鑒定
    復(fù)合材料無損檢測探討
    電子測試(2017年11期)2017-12-15 08:57:13
    油氣儲運專業(yè)仿真系統(tǒng)的開發(fā)和應(yīng)用
    重大交通設(shè)施項目風(fēng)險復(fù)雜動態(tài)交互演化機理與仿真分析
    預(yù)測(2016年3期)2016-12-29 18:57:28
    基于前饋補償?shù)闹鄙龣C控制律設(shè)計與仿真
    航空兵器(2016年4期)2016-11-28 21:34:37
    基于計算機系統(tǒng)仿真的電力機車模擬駕駛裝置研究
    大滯后系統(tǒng)中單神經(jīng)元PID控制器設(shè)計
    雙饋感應(yīng)發(fā)電機不脫網(wǎng)運行的系統(tǒng)仿真
    亚洲三级黄色毛片| 少妇的逼水好多| 亚洲国产精品成人综合色| 欧美日韩乱码在线| 午夜激情欧美在线| 国产在线精品亚洲第一网站| 三级男女做爰猛烈吃奶摸视频| 特级一级黄色大片| 婷婷精品国产亚洲av在线| 男女做爰动态图高潮gif福利片| 亚州av有码| 嫩草影院新地址| 亚洲av免费高清在线观看| 国产69精品久久久久777片| 国产v大片淫在线免费观看| 中文字幕久久专区| 国产熟女欧美一区二区| 国产蜜桃级精品一区二区三区| 国产成人91sexporn| 嫩草影院精品99| 男插女下体视频免费在线播放| 别揉我奶头~嗯~啊~动态视频| 乱系列少妇在线播放| 级片在线观看| 黄色日韩在线| 国产一区二区在线观看日韩| 女人十人毛片免费观看3o分钟| 国产亚洲精品久久久com| 成人鲁丝片一二三区免费| 一级毛片aaaaaa免费看小| 在线a可以看的网站| 亚洲成人精品中文字幕电影| 免费看美女性在线毛片视频| 欧美色欧美亚洲另类二区| 成人综合一区亚洲| av在线蜜桃| 男女做爰动态图高潮gif福利片| 少妇的逼好多水| 免费高清视频大片| 亚洲高清免费不卡视频| 欧洲精品卡2卡3卡4卡5卡区| 日韩,欧美,国产一区二区三区 | 午夜亚洲福利在线播放| 久久久久久九九精品二区国产| 免费高清视频大片| 99久久中文字幕三级久久日本| 国产一区二区三区av在线 | 日本撒尿小便嘘嘘汇集6| 久久精品夜夜夜夜夜久久蜜豆| 精品福利观看| 亚洲av成人精品一区久久| 精品一区二区三区av网在线观看| 国产精品国产高清国产av| 日日啪夜夜撸| 一级毛片久久久久久久久女| aaaaa片日本免费| 青春草视频在线免费观看| 久久精品夜夜夜夜夜久久蜜豆| 国国产精品蜜臀av免费| 国产成人aa在线观看| 日本一本二区三区精品| 国产成人精品久久久久久| 日韩制服骚丝袜av| 成人无遮挡网站| 91精品国产九色| 看黄色毛片网站| 久久精品91蜜桃| 99久久精品热视频| 一本久久中文字幕| 男女那种视频在线观看| 全区人妻精品视频| 欧美激情久久久久久爽电影| 男人的好看免费观看在线视频| 成人高潮视频无遮挡免费网站| 老司机影院成人| 久久精品人妻少妇| 亚洲无线观看免费| 国产在视频线在精品| 亚洲国产精品国产精品| 色尼玛亚洲综合影院| 我的女老师完整版在线观看| 国产伦精品一区二区三区四那| 日韩成人av中文字幕在线观看 | 久久韩国三级中文字幕| 亚洲人成网站在线播| 在线国产一区二区在线| 国产伦一二天堂av在线观看| 麻豆av噜噜一区二区三区| 一个人看的www免费观看视频| 别揉我奶头~嗯~啊~动态视频| 国产精品一及| 3wmmmm亚洲av在线观看| 国产乱人视频| 成人特级黄色片久久久久久久| 一进一出抽搐gif免费好疼| 草草在线视频免费看| 狠狠狠狠99中文字幕| 最新在线观看一区二区三区| 日本熟妇午夜| 久久久成人免费电影| 18禁在线播放成人免费| 免费黄网站久久成人精品| 免费看光身美女| 搡老岳熟女国产| 日日摸夜夜添夜夜爱| 久久精品国产99精品国产亚洲性色| 国产精品永久免费网站| 免费观看在线日韩| 日韩强制内射视频| 午夜福利成人在线免费观看| 在线天堂最新版资源| 久久亚洲精品不卡| 免费高清视频大片| or卡值多少钱| 一夜夜www| 午夜福利视频1000在线观看| 国产精品久久久久久久久免| 亚洲人成网站高清观看| 日韩成人伦理影院| 看十八女毛片水多多多| 国产精品久久久久久久电影| 精品一区二区三区人妻视频| 国产免费一级a男人的天堂| 在线天堂最新版资源| 日本五十路高清| av在线天堂中文字幕| 亚洲一区高清亚洲精品| 老熟妇乱子伦视频在线观看| 久久久久久久久久黄片| 少妇人妻精品综合一区二区 | aaaaa片日本免费| 欧美高清成人免费视频www| 国产精品1区2区在线观看.| 精品人妻偷拍中文字幕| 观看美女的网站| 人妻久久中文字幕网| 欧美不卡视频在线免费观看| 97超视频在线观看视频| 一级黄片播放器| 国产成人一区二区在线| 亚洲精品一区av在线观看| 日韩成人av中文字幕在线观看 | avwww免费| 欧美xxxx黑人xx丫x性爽| 国产在视频线在精品| 免费人成在线观看视频色| 欧美又色又爽又黄视频| 男女视频在线观看网站免费| 日韩大尺度精品在线看网址| 欧美日韩在线观看h| 两个人视频免费观看高清| 亚洲图色成人| 国产成人精品久久久久久| 色哟哟·www| 欧美xxxx性猛交bbbb| 免费av观看视频| 无遮挡黄片免费观看| 亚洲欧美日韩卡通动漫| 精品国内亚洲2022精品成人| 99国产极品粉嫩在线观看| 免费观看人在逋| 国产精品永久免费网站| 99久久久亚洲精品蜜臀av| 久久中文看片网| 毛片女人毛片| 国产在线男女| 亚洲精品成人久久久久久| а√天堂www在线а√下载| 中文字幕免费在线视频6| 伦理电影大哥的女人| 99久久精品热视频| 国产精品美女特级片免费视频播放器| 欧美一级a爱片免费观看看| 你懂的网址亚洲精品在线观看 | 国产 一区 欧美 日韩| 亚洲va在线va天堂va国产| 性色avwww在线观看| 国产中年淑女户外野战色| 免费在线观看成人毛片| 中文在线观看免费www的网站| 成熟少妇高潮喷水视频| 亚洲丝袜综合中文字幕| 精品久久久久久久久亚洲| 午夜福利在线观看吧| 免费看av在线观看网站| 久久亚洲国产成人精品v| 国产不卡一卡二| 男插女下体视频免费在线播放| 在线播放国产精品三级| 97碰自拍视频| 91在线精品国自产拍蜜月| 久久久久精品国产欧美久久久| 九九在线视频观看精品| 日韩强制内射视频| 午夜激情福利司机影院| 特大巨黑吊av在线直播| 18禁裸乳无遮挡免费网站照片| 男人的好看免费观看在线视频| 免费在线观看成人毛片| 波多野结衣巨乳人妻| 亚洲一区二区三区色噜噜| 国产精品日韩av在线免费观看| 亚洲av免费高清在线观看| 国产在线男女| 亚洲国产精品sss在线观看| 校园人妻丝袜中文字幕| 秋霞在线观看毛片| 免费看光身美女| 亚洲精品日韩在线中文字幕 | 婷婷精品国产亚洲av在线| 丰满乱子伦码专区| 精品人妻熟女av久视频| 精品久久久久久成人av| 欧美中文日本在线观看视频| 日韩欧美在线乱码| 日韩国内少妇激情av| 人人妻人人澡欧美一区二区| 亚洲内射少妇av| 午夜精品国产一区二区电影 | 国产乱人视频| 亚洲精品久久国产高清桃花| 免费看光身美女| 欧美精品国产亚洲| 如何舔出高潮| 97碰自拍视频| 在线观看午夜福利视频| 搡老熟女国产l中国老女人| 美女黄网站色视频| 九九爱精品视频在线观看| 久久精品人妻少妇| 好男人在线观看高清免费视频| 一级毛片aaaaaa免费看小| 一区福利在线观看| 桃色一区二区三区在线观看| 中文亚洲av片在线观看爽| 综合色丁香网| 真人做人爱边吃奶动态| 97超视频在线观看视频| 色尼玛亚洲综合影院| 日本五十路高清| 国产精华一区二区三区| 国产探花极品一区二区| a级毛片免费高清观看在线播放| 女的被弄到高潮叫床怎么办| 国内少妇人妻偷人精品xxx网站| 亚洲七黄色美女视频| 中文字幕免费在线视频6| 99久久成人亚洲精品观看| 午夜视频国产福利| 午夜免费男女啪啪视频观看 | 日韩中字成人| 中国国产av一级| 一级av片app| 少妇猛男粗大的猛烈进出视频 | 日韩精品青青久久久久久| 国产精品日韩av在线免费观看| 51国产日韩欧美| av女优亚洲男人天堂| 日韩成人伦理影院| 波多野结衣高清作品| 一级毛片我不卡| 美女被艹到高潮喷水动态| 免费av不卡在线播放| 色综合站精品国产| 成年免费大片在线观看| 成人美女网站在线观看视频| 97碰自拍视频| 插阴视频在线观看视频| 亚洲丝袜综合中文字幕| 国产精品av视频在线免费观看| 97热精品久久久久久| 日韩欧美国产在线观看| av天堂在线播放| 中文字幕熟女人妻在线| 色尼玛亚洲综合影院| 少妇熟女aⅴ在线视频| 大型黄色视频在线免费观看| 久久久久国内视频| 免费人成在线观看视频色| 一a级毛片在线观看| 亚洲aⅴ乱码一区二区在线播放| 99热6这里只有精品| 久久亚洲国产成人精品v| 精品一区二区免费观看| 简卡轻食公司| 看十八女毛片水多多多| 久久精品综合一区二区三区| 最近手机中文字幕大全| 网址你懂的国产日韩在线| 菩萨蛮人人尽说江南好唐韦庄 | 一边摸一边抽搐一进一小说| 特大巨黑吊av在线直播| 少妇人妻精品综合一区二区 | 亚洲欧美清纯卡通| 成熟少妇高潮喷水视频| 国产精品人妻久久久久久| 欧美绝顶高潮抽搐喷水| 国产精品1区2区在线观看.| 婷婷亚洲欧美| 我的老师免费观看完整版| 菩萨蛮人人尽说江南好唐韦庄 | 高清毛片免费观看视频网站| 天堂动漫精品| av在线天堂中文字幕| 亚洲国产色片| 欧美国产日韩亚洲一区| 夜夜看夜夜爽夜夜摸| 乱系列少妇在线播放| 精品日产1卡2卡| 尾随美女入室| 国产一区二区在线观看日韩| 亚洲精品国产成人久久av| 亚洲最大成人手机在线| 国产激情偷乱视频一区二区| 白带黄色成豆腐渣| 人人妻人人澡欧美一区二区| 97在线视频观看| aaaaa片日本免费| 日韩一区二区视频免费看| 亚洲专区国产一区二区| 国产午夜精品久久久久久一区二区三区 | 啦啦啦啦在线视频资源| 免费看日本二区| 精品久久久久久久久久久久久| 国产久久久一区二区三区| 搞女人的毛片| 午夜福利在线在线| 91麻豆精品激情在线观看国产| 丰满的人妻完整版| 精品久久国产蜜桃| 欧美zozozo另类| 国产女主播在线喷水免费视频网站 | 国产精品1区2区在线观看.| 悠悠久久av| 亚州av有码| 亚洲精品影视一区二区三区av| 日韩av在线大香蕉| 波多野结衣高清作品| av卡一久久| 男人和女人高潮做爰伦理| 国产人妻一区二区三区在| 亚洲欧美精品综合久久99| 亚洲人成网站在线播| 欧美激情国产日韩精品一区| 三级经典国产精品| 少妇被粗大猛烈的视频| 亚洲av成人av| 国内精品宾馆在线| 国产69精品久久久久777片| 男女边吃奶边做爰视频| 自拍偷自拍亚洲精品老妇| 国产伦精品一区二区三区视频9| 久久精品国产清高在天天线| 国产成人91sexporn| 深爱激情五月婷婷| 日韩 亚洲 欧美在线| 成年免费大片在线观看| 欧美丝袜亚洲另类| 国产 一区 欧美 日韩| 久久久久久久久久成人| av福利片在线观看| 国产成人91sexporn| 久久精品国产自在天天线| 丝袜美腿在线中文| 真实男女啪啪啪动态图| 女同久久另类99精品国产91| 亚洲电影在线观看av| 国国产精品蜜臀av免费| 中文字幕免费在线视频6| 又粗又爽又猛毛片免费看| 国产精品女同一区二区软件| aaaaa片日本免费| 看片在线看免费视频| 国产伦一二天堂av在线观看| 国产成人精品久久久久久| 国产av不卡久久| 久久精品国产亚洲网站| 在线国产一区二区在线| 人妻夜夜爽99麻豆av| 99国产极品粉嫩在线观看| 俄罗斯特黄特色一大片| 99国产精品一区二区蜜桃av| 国产精品一二三区在线看| 欧美另类亚洲清纯唯美| 天天一区二区日本电影三级| 99国产极品粉嫩在线观看| 亚洲精品色激情综合| 欧美性感艳星| 国产又黄又爽又无遮挡在线| 搡老岳熟女国产| 一本久久中文字幕| 99热这里只有是精品50| 美女免费视频网站| 国内精品宾馆在线| 波野结衣二区三区在线| 狂野欧美激情性xxxx在线观看| 国产精品久久电影中文字幕| 日日啪夜夜撸| 日韩亚洲欧美综合| 熟妇人妻久久中文字幕3abv| 岛国在线免费视频观看| 亚洲在线观看片| 亚洲精品日韩av片在线观看| 少妇裸体淫交视频免费看高清| 51国产日韩欧美| 精品一区二区三区av网在线观看| 国产精品电影一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 国产一区二区三区av在线 | 亚洲自偷自拍三级| 欧美色视频一区免费| 亚洲av中文字字幕乱码综合| 国产在线男女| 又粗又爽又猛毛片免费看| 一边摸一边抽搐一进一小说| 精品一区二区三区av网在线观看| 看黄色毛片网站| 男人舔女人下体高潮全视频| 亚洲av熟女| 欧美一级a爱片免费观看看| 一进一出抽搐gif免费好疼| 看非洲黑人一级黄片| 久久韩国三级中文字幕| 久久人人爽人人片av| av专区在线播放| 日本与韩国留学比较| 亚洲一级一片aⅴ在线观看| 日韩强制内射视频| 国产综合懂色| 亚洲,欧美,日韩| 少妇人妻一区二区三区视频| 亚洲va在线va天堂va国产| 精品国内亚洲2022精品成人| 国产一区二区在线观看日韩| 女生性感内裤真人,穿戴方法视频| 亚洲综合色惰| 免费人成视频x8x8入口观看| 日本在线视频免费播放| 成年女人毛片免费观看观看9| 国产欧美日韩精品一区二区| 免费av不卡在线播放| 成人国产麻豆网| 狂野欧美激情性xxxx在线观看| 1024手机看黄色片| 非洲黑人性xxxx精品又粗又长| 国产又黄又爽又无遮挡在线| 亚洲中文字幕一区二区三区有码在线看| 一级a爱片免费观看的视频| 最近视频中文字幕2019在线8| 亚洲成人久久爱视频| 中国美白少妇内射xxxbb| 亚洲自拍偷在线| 给我免费播放毛片高清在线观看| 美女 人体艺术 gogo| 国产av麻豆久久久久久久| av在线观看视频网站免费| 国产av不卡久久| 国产真实乱freesex| 中文资源天堂在线| av国产免费在线观看| 99国产精品一区二区蜜桃av| 尤物成人国产欧美一区二区三区| 三级经典国产精品| 亚洲av.av天堂| 丝袜美腿在线中文| 美女xxoo啪啪120秒动态图| 免费观看的影片在线观看| 男女之事视频高清在线观看| 又爽又黄无遮挡网站| 婷婷精品国产亚洲av| 男女边吃奶边做爰视频| 女人被狂操c到高潮| 不卡一级毛片| 欧美一区二区亚洲| 天堂影院成人在线观看| www日本黄色视频网| 亚洲精品久久国产高清桃花| 精品久久久久久久久久免费视频| 天堂av国产一区二区熟女人妻| 久久国产乱子免费精品| 日日干狠狠操夜夜爽| 欧美最新免费一区二区三区| 99国产精品一区二区蜜桃av| 搞女人的毛片| 噜噜噜噜噜久久久久久91| 国产v大片淫在线免费观看| 精品福利观看| 免费av毛片视频| 桃色一区二区三区在线观看| 国产午夜福利久久久久久| 久久午夜福利片| 亚洲av中文字字幕乱码综合| 色噜噜av男人的天堂激情| 国产美女午夜福利| 天堂网av新在线| 俺也久久电影网| 成人性生交大片免费视频hd| 大又大粗又爽又黄少妇毛片口| 天堂影院成人在线观看| 99热只有精品国产| 最后的刺客免费高清国语| 欧美成人免费av一区二区三区| 美女高潮的动态| 国产探花在线观看一区二区| 国产视频一区二区在线看| 露出奶头的视频| 秋霞在线观看毛片| 99精品在免费线老司机午夜| 免费观看在线日韩| 嫩草影院新地址| 国产精品久久久久久久久免| 少妇熟女aⅴ在线视频| 欧美又色又爽又黄视频| 亚洲精品粉嫩美女一区| 亚洲国产精品国产精品| 99视频精品全部免费 在线| 丰满人妻一区二区三区视频av| 亚洲高清免费不卡视频| 亚洲精品一卡2卡三卡4卡5卡| 国产成人91sexporn| 国产淫片久久久久久久久| 亚洲色图av天堂| 成年女人永久免费观看视频| 成人二区视频| 18禁黄网站禁片免费观看直播| 日日干狠狠操夜夜爽| 亚洲欧美日韩高清在线视频| 少妇被粗大猛烈的视频| 亚洲欧美中文字幕日韩二区| 少妇熟女aⅴ在线视频| 亚州av有码| 夜夜夜夜夜久久久久| 少妇高潮的动态图| АⅤ资源中文在线天堂| 国产精品伦人一区二区| 国内精品一区二区在线观看| 久久人人爽人人爽人人片va| 女人被狂操c到高潮| 欧美潮喷喷水| ponron亚洲| videossex国产| 18+在线观看网站| 国产精品精品国产色婷婷| 一本精品99久久精品77| 日韩精品有码人妻一区| 成人av在线播放网站| 老熟妇乱子伦视频在线观看| 少妇猛男粗大的猛烈进出视频 | 天美传媒精品一区二区| 啦啦啦观看免费观看视频高清| 12—13女人毛片做爰片一| 欧美日韩乱码在线| 男人的好看免费观看在线视频| 国内精品久久久久精免费| 日本黄大片高清| 哪里可以看免费的av片| 国内精品久久久久精免费| 日本 av在线| 精品一区二区三区视频在线| 永久网站在线| 久久6这里有精品| 久久精品国产鲁丝片午夜精品| 夜夜看夜夜爽夜夜摸| 一级毛片电影观看 | 国产精品,欧美在线| 精品一区二区三区视频在线| 大香蕉久久网| 色吧在线观看| 在现免费观看毛片| a级一级毛片免费在线观看| 亚洲精品日韩av片在线观看| 夜夜爽天天搞| 日本三级黄在线观看| 亚洲成人av在线免费| 丰满人妻一区二区三区视频av| 久久韩国三级中文字幕| 亚洲久久久久久中文字幕| 国语自产精品视频在线第100页| 亚洲美女黄片视频| 中国美女看黄片| 亚洲国产日韩欧美精品在线观看| avwww免费| 国产日本99.免费观看| 国产精品1区2区在线观看.| 噜噜噜噜噜久久久久久91| 偷拍熟女少妇极品色| 日本一本二区三区精品| 日本色播在线视频| 一本久久中文字幕| 久久亚洲国产成人精品v| 亚洲av免费在线观看| 国产真实乱freesex| 波野结衣二区三区在线| 一级毛片电影观看 | 岛国在线免费视频观看| 亚洲在线观看片| 有码 亚洲区| 国内久久婷婷六月综合欲色啪| 内射极品少妇av片p| 69av精品久久久久久| 91在线观看av| 最近中文字幕高清免费大全6| 啦啦啦观看免费观看视频高清| 国产精品久久久久久精品电影| 免费观看的影片在线观看| 网址你懂的国产日韩在线| 国产成人影院久久av| 国产伦在线观看视频一区| 亚洲18禁久久av| 国产激情偷乱视频一区二区| 国产精品免费一区二区三区在线| 久久久久九九精品影院| 久久鲁丝午夜福利片| 九色成人免费人妻av| 狂野欧美白嫩少妇大欣赏|