• <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)仿真
    人人妻人人澡人人看| 国产精品综合久久久久久久免费 | 久久精品亚洲熟妇少妇任你| 精品人妻1区二区| av欧美777| 国产精品自产拍在线观看55亚洲| 欧美精品一区二区免费开放| 最新美女视频免费是黄的| 两性午夜刺激爽爽歪歪视频在线观看 | 女人高潮潮喷娇喘18禁视频| 国产亚洲欧美在线一区二区| 久久香蕉激情| 可以在线观看毛片的网站| 9热在线视频观看99| 日韩免费av在线播放| 黄色女人牲交| 国产精品 国内视频| 在线av久久热| 国产精品久久久久成人av| 18禁观看日本| 色综合站精品国产| 热99re8久久精品国产| 国产人伦9x9x在线观看| 美国免费a级毛片| 一进一出抽搐gif免费好疼 | 国产乱人伦免费视频| 一进一出抽搐gif免费好疼 | 老熟妇乱子伦视频在线观看| 亚洲国产看品久久| 亚洲黑人精品在线| 欧美中文综合在线视频| 欧美国产精品va在线观看不卡| 国产欧美日韩一区二区三区在线| 欧美激情极品国产一区二区三区| 成人黄色视频免费在线看| 免费av中文字幕在线| 不卡av一区二区三区| www.自偷自拍.com| 黑丝袜美女国产一区| 午夜精品久久久久久毛片777| 成人18禁高潮啪啪吃奶动态图| 国产亚洲精品一区二区www| 成人av一区二区三区在线看| 亚洲成av片中文字幕在线观看| 最近最新免费中文字幕在线| 久久午夜综合久久蜜桃| 国产精品98久久久久久宅男小说| 久久精品国产亚洲av香蕉五月| 精品免费久久久久久久清纯| 国产激情久久老熟女| 久久精品亚洲熟妇少妇任你| 国产高清激情床上av| 亚洲av第一区精品v没综合| 操美女的视频在线观看| 欧美人与性动交α欧美软件| 欧美在线黄色| 一级,二级,三级黄色视频| 亚洲欧美精品综合久久99| 99国产精品99久久久久| 婷婷精品国产亚洲av在线| 琪琪午夜伦伦电影理论片6080| 欧美色视频一区免费| 久久久久久久久免费视频了| 国产高清视频在线播放一区| 色综合站精品国产| 欧美中文日本在线观看视频| 曰老女人黄片| 18禁国产床啪视频网站| 99在线视频只有这里精品首页| 久久影院123| 水蜜桃什么品种好| 久久中文看片网| 日韩高清综合在线| 脱女人内裤的视频| 神马国产精品三级电影在线观看 | 亚洲色图av天堂| 美女午夜性视频免费| 日本a在线网址| 精品欧美一区二区三区在线| 成年版毛片免费区| 在线免费观看的www视频| 又大又爽又粗| 亚洲色图综合在线观看| 波多野结衣高清无吗| 亚洲av日韩精品久久久久久密| 中文字幕人妻熟女乱码| 日韩免费高清中文字幕av| 三上悠亚av全集在线观看| 免费在线观看完整版高清| 桃红色精品国产亚洲av| 亚洲 欧美一区二区三区| 亚洲精品中文字幕在线视频| 女性生殖器流出的白浆| 久久草成人影院| 好看av亚洲va欧美ⅴa在| 亚洲美女黄片视频| 黑人巨大精品欧美一区二区蜜桃| 99精品欧美一区二区三区四区| www.www免费av| 欧美大码av| 一边摸一边抽搐一进一出视频| 一级a爱片免费观看的视频| 免费观看人在逋| 精品午夜福利视频在线观看一区| 黄色成人免费大全| 91精品三级在线观看| 91国产中文字幕| 9色porny在线观看| 97碰自拍视频| 久久99一区二区三区| 91麻豆av在线| 亚洲人成网站在线播放欧美日韩| 成年人免费黄色播放视频| 精品熟女少妇八av免费久了| 亚洲自偷自拍图片 自拍| 国产精品久久久av美女十八| a在线观看视频网站| 淫妇啪啪啪对白视频| 国产精品综合久久久久久久免费 | 一夜夜www| 国产97色在线日韩免费| 亚洲精品成人av观看孕妇| 免费av毛片视频| 久久精品亚洲av国产电影网| 久久久久久久久中文| 国产av在哪里看| 黑丝袜美女国产一区| 中文字幕色久视频| 在线观看66精品国产| 欧美黄色片欧美黄色片| 免费在线观看黄色视频的| 精品国内亚洲2022精品成人| 新久久久久国产一级毛片| 91av网站免费观看| 国产av一区在线观看免费| 国产视频一区二区在线看| 精品乱码久久久久久99久播| 免费在线观看亚洲国产| 色老头精品视频在线观看| 亚洲国产看品久久| 91国产中文字幕| av国产精品久久久久影院| 丰满饥渴人妻一区二区三| 99精品欧美一区二区三区四区| 嫁个100分男人电影在线观看| 国产精品一区二区在线不卡| 伊人久久大香线蕉亚洲五| 国产精品爽爽va在线观看网站 | 欧美在线一区亚洲| 国产在线精品亚洲第一网站| av中文乱码字幕在线| 91av网站免费观看| 午夜久久久在线观看| 9热在线视频观看99| 欧美激情 高清一区二区三区| 久9热在线精品视频| 91在线观看av| 国产av一区二区精品久久| 久久天躁狠狠躁夜夜2o2o| 国产人伦9x9x在线观看| 午夜免费成人在线视频| 电影成人av| 久久国产精品人妻蜜桃| 亚洲av成人不卡在线观看播放网| 日日干狠狠操夜夜爽| 伦理电影免费视频| 国产精品98久久久久久宅男小说| 1024香蕉在线观看| 麻豆国产av国片精品| 色综合欧美亚洲国产小说| 手机成人av网站| 深夜精品福利| 自拍欧美九色日韩亚洲蝌蚪91| 侵犯人妻中文字幕一二三四区| 欧美日本中文国产一区发布| 婷婷丁香在线五月| ponron亚洲| 亚洲九九香蕉| 精品一区二区三区av网在线观看| 国产精品影院久久| 国产成人影院久久av| 一级毛片高清免费大全| avwww免费| www.www免费av| 亚洲av片天天在线观看| 在线免费观看的www视频| 色婷婷av一区二区三区视频| 超碰97精品在线观看| 嫩草影院精品99| 丰满人妻熟妇乱又伦精品不卡| 国产精品综合久久久久久久免费 | 制服诱惑二区| 欧美乱色亚洲激情| 91精品三级在线观看| www国产在线视频色| 欧美日韩av久久| 欧美日本中文国产一区发布| 高清在线国产一区| 日韩成人在线观看一区二区三区| 后天国语完整版免费观看| 欧美乱码精品一区二区三区| 日本黄色日本黄色录像| 人妻久久中文字幕网| 热re99久久精品国产66热6| av在线天堂中文字幕 | 无遮挡黄片免费观看| 国产精品久久久久久人妻精品电影| 人人妻人人添人人爽欧美一区卜| 久热这里只有精品99| 人人妻人人澡人人看| 九色亚洲精品在线播放| 久久久精品欧美日韩精品| 一进一出抽搐动态| 侵犯人妻中文字幕一二三四区| 19禁男女啪啪无遮挡网站| 欧美日韩福利视频一区二区| 亚洲国产中文字幕在线视频| 亚洲五月婷婷丁香| 国产极品粉嫩免费观看在线| 国产97色在线日韩免费| 中文字幕av电影在线播放| 亚洲成人免费电影在线观看| 久热这里只有精品99| 老鸭窝网址在线观看| 欧美成人免费av一区二区三区| 美国免费a级毛片| 国产在线观看jvid| 高清黄色对白视频在线免费看| 女人被狂操c到高潮| 超碰成人久久| 成人三级黄色视频| 男女下面插进去视频免费观看| 亚洲欧美精品综合久久99| 国产精品国产高清国产av| 国产真人三级小视频在线观看| av网站免费在线观看视频| 久久久久精品国产欧美久久久| 亚洲人成网站在线播放欧美日韩| 在线免费观看的www视频| 性欧美人与动物交配| 亚洲国产欧美日韩在线播放| 免费人成视频x8x8入口观看| 搡老熟女国产l中国老女人| 国产成年人精品一区二区 | 少妇粗大呻吟视频| 成人手机av| 国产精品久久久久成人av| 久久久水蜜桃国产精品网| 欧美日韩亚洲国产一区二区在线观看| 99国产极品粉嫩在线观看| 成人手机av| 免费日韩欧美在线观看| 久久人妻福利社区极品人妻图片| 黄色怎么调成土黄色| 久久狼人影院| 亚洲免费av在线视频| 999精品在线视频| 亚洲精品中文字幕一二三四区| 极品人妻少妇av视频| 久久香蕉精品热| 大陆偷拍与自拍| 搡老乐熟女国产| 村上凉子中文字幕在线| 99久久国产精品久久久| 欧美亚洲日本最大视频资源| 久久天躁狠狠躁夜夜2o2o| 看片在线看免费视频| 国产高清激情床上av| 国产区一区二久久| 亚洲成人久久性| 亚洲国产欧美一区二区综合| 欧美国产精品va在线观看不卡| 色综合站精品国产| 丁香六月欧美| 午夜福利,免费看| 国产单亲对白刺激| 欧美+亚洲+日韩+国产| 国产精品久久视频播放| 午夜福利影视在线免费观看| 久久精品国产清高在天天线| 成人亚洲精品av一区二区 | 老司机福利观看| 波多野结衣高清无吗| 久久精品亚洲av国产电影网| 变态另类成人亚洲欧美熟女 | 欧美丝袜亚洲另类 | 90打野战视频偷拍视频| 国产成人免费无遮挡视频| 亚洲精品久久成人aⅴ小说| 在线av久久热| 欧美黄色淫秽网站| 亚洲国产欧美日韩在线播放| 大香蕉久久成人网| 黄色成人免费大全| 电影成人av| 不卡av一区二区三区| 午夜福利在线免费观看网站| 黄色怎么调成土黄色| 亚洲激情在线av| 大码成人一级视频| 天天影视国产精品| 在线av久久热| 天天躁狠狠躁夜夜躁狠狠躁| 性欧美人与动物交配| 波多野结衣高清无吗| 久久久国产成人免费| 日韩中文字幕欧美一区二区| av超薄肉色丝袜交足视频| 亚洲自偷自拍图片 自拍| 国产一区二区三区在线臀色熟女 | 99国产综合亚洲精品| 免费在线观看日本一区| 欧美午夜高清在线| 久久这里只有精品19| 老熟妇乱子伦视频在线观看| 99在线人妻在线中文字幕| 成熟少妇高潮喷水视频| 自线自在国产av| 在线观看免费日韩欧美大片| 男男h啪啪无遮挡| 国产高清激情床上av| 黄色片一级片一级黄色片| 色播在线永久视频| 男女做爰动态图高潮gif福利片 | 精品国产乱码久久久久久男人| 日本三级黄在线观看| 色哟哟哟哟哟哟| 51午夜福利影视在线观看| 国产aⅴ精品一区二区三区波| 神马国产精品三级电影在线观看 | 91大片在线观看| 午夜影院日韩av| 久久天躁狠狠躁夜夜2o2o| 看片在线看免费视频| 两个人免费观看高清视频| 日韩有码中文字幕| 亚洲中文字幕日韩| 国产成人精品无人区| 一级毛片女人18水好多| 在线天堂中文资源库| 黄片播放在线免费| 亚洲熟妇中文字幕五十中出 | 午夜福利免费观看在线| 丁香六月欧美| 午夜福利在线免费观看网站| 国产亚洲精品一区二区www| 嫩草影院精品99| 亚洲午夜理论影院| 国产精品爽爽va在线观看网站 | 操美女的视频在线观看| 亚洲自拍偷在线| x7x7x7水蜜桃| 99国产极品粉嫩在线观看| 国产精品秋霞免费鲁丝片| 大型黄色视频在线免费观看| 亚洲欧美激情综合另类| 久久国产乱子伦精品免费另类| 精品无人区乱码1区二区| 纯流量卡能插随身wifi吗| 如日韩欧美国产精品一区二区三区| 午夜两性在线视频| 成人手机av| 超碰97精品在线观看| 香蕉国产在线看| 午夜两性在线视频| 真人一进一出gif抽搐免费| 精品国产亚洲在线| 国产熟女午夜一区二区三区| 12—13女人毛片做爰片一| 午夜精品在线福利| 丝袜在线中文字幕| 99香蕉大伊视频| 日韩免费av在线播放| bbb黄色大片| 国产精品综合久久久久久久免费 | 自拍欧美九色日韩亚洲蝌蚪91| 国产xxxxx性猛交| 自拍欧美九色日韩亚洲蝌蚪91| 国产真人三级小视频在线观看| av网站在线播放免费| 男女下面进入的视频免费午夜 | a在线观看视频网站| 亚洲成人免费电影在线观看| 久久亚洲真实| 国产视频一区二区在线看| 男人的好看免费观看在线视频 | 国产精品秋霞免费鲁丝片| 亚洲精品在线观看二区| www国产在线视频色| 国产xxxxx性猛交| 高清av免费在线| 婷婷六月久久综合丁香| 正在播放国产对白刺激| 十八禁网站免费在线| 精品国内亚洲2022精品成人| 免费高清在线观看日韩| 欧美国产精品va在线观看不卡| 国产成人精品久久二区二区91| 亚洲av五月六月丁香网| 男人的好看免费观看在线视频 | 久久久久国产一级毛片高清牌| 成人三级黄色视频| 免费av中文字幕在线| 亚洲精品一二三| 欧美久久黑人一区二区| 久久久久久久午夜电影 | 午夜福利在线观看吧| 中文亚洲av片在线观看爽| 成人手机av| 女人被躁到高潮嗷嗷叫费观| aaaaa片日本免费| 亚洲美女黄片视频| 法律面前人人平等表现在哪些方面| 在线永久观看黄色视频| 国产日韩一区二区三区精品不卡| 啦啦啦 在线观看视频| 一级毛片高清免费大全| 一夜夜www| 午夜激情av网站| 丝袜人妻中文字幕| 成人av一区二区三区在线看| av片东京热男人的天堂| 日韩精品免费视频一区二区三区| 两个人看的免费小视频| 欧美激情极品国产一区二区三区| 91九色精品人成在线观看| 精品免费久久久久久久清纯| 少妇 在线观看| 香蕉国产在线看| 91字幕亚洲| 国产97色在线日韩免费| 精品人妻在线不人妻| 亚洲欧洲精品一区二区精品久久久| 国产成年人精品一区二区 | 欧美在线一区亚洲| 久久精品91蜜桃| 国产精品美女特级片免费视频播放器 | 女人高潮潮喷娇喘18禁视频| 国产亚洲欧美精品永久| 91麻豆精品激情在线观看国产 | 欧美日韩亚洲综合一区二区三区_| 亚洲自偷自拍图片 自拍| av超薄肉色丝袜交足视频| 精品一区二区三区视频在线观看免费 | 9色porny在线观看| 国产成人精品在线电影| 久久精品亚洲熟妇少妇任你| 午夜免费激情av| 日韩免费av在线播放| 欧美丝袜亚洲另类 | 大香蕉久久成人网| 日韩欧美国产一区二区入口| 欧美日韩亚洲国产一区二区在线观看| 免费在线观看完整版高清| 女人爽到高潮嗷嗷叫在线视频| 国产黄色免费在线视频| 成人特级黄色片久久久久久久| 中文字幕人妻熟女乱码| 成人三级黄色视频| 久久影院123| 亚洲中文av在线| 欧美人与性动交α欧美软件| 在线观看午夜福利视频| 91字幕亚洲| 中出人妻视频一区二区| 在线免费观看的www视频| 美女福利国产在线| 欧美日韩乱码在线| 97人妻天天添夜夜摸| 青草久久国产| 久久欧美精品欧美久久欧美| 久久久国产一区二区| 99精品欧美一区二区三区四区| 在线国产一区二区在线| 啪啪无遮挡十八禁网站| 国产无遮挡羞羞视频在线观看| 美女 人体艺术 gogo| 亚洲欧美日韩无卡精品| 国产aⅴ精品一区二区三区波| 国产精品日韩av在线免费观看 | 日韩精品免费视频一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 一个人免费在线观看的高清视频| 日韩精品中文字幕看吧| 波多野结衣一区麻豆| 最好的美女福利视频网| 国产高清videossex| 亚洲熟妇熟女久久| 欧美人与性动交α欧美软件| 一区二区三区精品91| 神马国产精品三级电影在线观看 | 美女高潮到喷水免费观看| 欧美人与性动交α欧美软件| 大香蕉久久成人网| 国产亚洲精品第一综合不卡| 国产精品一区二区在线不卡| 嫩草影院精品99| av免费在线观看网站| 亚洲人成电影免费在线| 99re在线观看精品视频| 亚洲欧美日韩高清在线视频| 欧美日韩乱码在线| 国产99白浆流出| 女人精品久久久久毛片| 国产精品免费一区二区三区在线| 精品卡一卡二卡四卡免费| av超薄肉色丝袜交足视频| 午夜精品国产一区二区电影| 91麻豆av在线| 99在线人妻在线中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 丝袜人妻中文字幕| 国产av在哪里看| 国产高清videossex| 国产免费现黄频在线看| 午夜激情av网站| 国产极品粉嫩免费观看在线| 国产主播在线观看一区二区| 成人精品一区二区免费| 亚洲aⅴ乱码一区二区在线播放 | 91麻豆av在线| 日本黄色日本黄色录像| 色老头精品视频在线观看| 69精品国产乱码久久久| 国产一区在线观看成人免费| 99精品在免费线老司机午夜| 欧美在线一区亚洲| 天天躁夜夜躁狠狠躁躁| √禁漫天堂资源中文www| 国产精品日韩av在线免费观看 | 另类亚洲欧美激情| 亚洲中文日韩欧美视频| 久久精品成人免费网站| 精品人妻1区二区| 欧美大码av| 国产精品av久久久久免费| 黄网站色视频无遮挡免费观看| 国产亚洲欧美98| 精品国产美女av久久久久小说| 久久久久久久久中文| 中国美女看黄片| 国产av在哪里看| 最新美女视频免费是黄的| 久久人妻福利社区极品人妻图片| 精品国产超薄肉色丝袜足j| 成人国语在线视频| 中文字幕人妻丝袜一区二区| 中文字幕精品免费在线观看视频| 在线观看免费日韩欧美大片| 美女大奶头视频| 精品国内亚洲2022精品成人| av网站免费在线观看视频| 成人国产一区最新在线观看| 多毛熟女@视频| 精品久久久久久久毛片微露脸| 国产亚洲av高清不卡| 久久久久久亚洲精品国产蜜桃av| 人人妻人人爽人人添夜夜欢视频| 欧美乱色亚洲激情| 欧美人与性动交α欧美精品济南到| 丁香欧美五月| 国产激情欧美一区二区| 国产精品久久视频播放| 91成人精品电影| 久久人妻av系列| 叶爱在线成人免费视频播放| 国产一区二区三区视频了| 一级片'在线观看视频| 琪琪午夜伦伦电影理论片6080| 一级黄色大片毛片| 国产精品亚洲av一区麻豆| 深夜精品福利| 一边摸一边抽搐一进一出视频| 国产成人av教育| 在线永久观看黄色视频| 欧美乱色亚洲激情| 99香蕉大伊视频| 国产免费av片在线观看野外av| 搡老熟女国产l中国老女人| 如日韩欧美国产精品一区二区三区| 日韩 欧美 亚洲 中文字幕| 美女福利国产在线| 黄片大片在线免费观看| 久久久国产欧美日韩av| 亚洲色图综合在线观看| 1024视频免费在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 可以在线观看毛片的网站| 国产精品 国内视频| 人人妻,人人澡人人爽秒播| 黄片播放在线免费| 人成视频在线观看免费观看| 午夜免费成人在线视频| 精品国产乱子伦一区二区三区| av网站免费在线观看视频| 女人高潮潮喷娇喘18禁视频| 怎么达到女性高潮| 色精品久久人妻99蜜桃| 新久久久久国产一级毛片| 看免费av毛片| 中亚洲国语对白在线视频| 日韩成人在线观看一区二区三区| 亚洲色图综合在线观看| 日韩高清综合在线| 国产一区在线观看成人免费| 色播在线永久视频| tocl精华| 一区二区日韩欧美中文字幕| 91在线观看av| 久久香蕉激情| 国产一卡二卡三卡精品| 久久久久国产精品人妻aⅴ院| 久久久精品国产亚洲av高清涩受| 99精品欧美一区二区三区四区| 国产精品一区二区精品视频观看|