• <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)仿真
    老司机靠b影院| 日韩免费高清中文字幕av| 观看美女的网站| 好男人视频免费观看在线| 夜夜骑夜夜射夜夜干| 老司机亚洲免费影院| 精品一区在线观看国产| 久久精品亚洲av国产电影网| 肉色欧美久久久久久久蜜桃| 欧美精品一区二区免费开放| 久久精品久久久久久噜噜老黄| 亚洲欧洲国产日韩| 可以免费在线观看a视频的电影网站 | 最近最新中文字幕大全免费视频 | 91aial.com中文字幕在线观看| 日韩制服骚丝袜av| 中文精品一卡2卡3卡4更新| 久久综合国产亚洲精品| 中国国产av一级| 成人三级做爰电影| 又粗又硬又长又爽又黄的视频| 国产麻豆69| 丝袜人妻中文字幕| 99精品久久久久人妻精品| 少妇的丰满在线观看| 欧美久久黑人一区二区| 亚洲一区中文字幕在线| 久热爱精品视频在线9| 成年av动漫网址| 亚洲国产看品久久| 大香蕉久久网| 日韩av不卡免费在线播放| 亚洲一卡2卡3卡4卡5卡精品中文| 极品人妻少妇av视频| 亚洲欧洲精品一区二区精品久久久 | 国产福利在线免费观看视频| 亚洲人成77777在线视频| 伊人亚洲综合成人网| 成人黄色视频免费在线看| 色婷婷av一区二区三区视频| 亚洲av电影在线观看一区二区三区| 亚洲综合精品二区| 久久久久精品性色| 下体分泌物呈黄色| 国产精品久久久久久精品电影小说| 亚洲三区欧美一区| 超碰97精品在线观看| 一边摸一边做爽爽视频免费| 中文字幕精品免费在线观看视频| 久久精品亚洲熟妇少妇任你| 精品久久蜜臀av无| 亚洲欧洲国产日韩| 成人黄色视频免费在线看| 欧美日韩成人在线一区二区| 国产成人精品久久二区二区91 | 欧美精品一区二区免费开放| 日韩大码丰满熟妇| 久久人妻熟女aⅴ| 国产一区二区激情短视频 | 久久这里只有精品19| www.熟女人妻精品国产| 亚洲色图综合在线观看| 热99国产精品久久久久久7| 99九九在线精品视频| av不卡在线播放| 最近手机中文字幕大全| 亚洲精品日本国产第一区| 国语对白做爰xxxⅹ性视频网站| 黄片播放在线免费| 国产av国产精品国产| 精品久久蜜臀av无| svipshipincom国产片| 国产视频首页在线观看| 久久久久久久久免费视频了| 久久午夜综合久久蜜桃| 久久热在线av| 观看美女的网站| 精品卡一卡二卡四卡免费| 国产免费一区二区三区四区乱码| 一级毛片我不卡| 午夜激情av网站| 在线免费观看不下载黄p国产| 国产老妇伦熟女老妇高清| 一边摸一边抽搐一进一出视频| 亚洲欧美中文字幕日韩二区| 女的被弄到高潮叫床怎么办| 日韩一本色道免费dvd| 国产精品麻豆人妻色哟哟久久| 国产精品一区二区在线不卡| 操出白浆在线播放| 亚洲一区二区三区欧美精品| 男女免费视频国产| 大陆偷拍与自拍| 国产一区二区激情短视频 | 精品午夜福利在线看| 日韩人妻精品一区2区三区| 亚洲精品aⅴ在线观看| 男女免费视频国产| 免费观看av网站的网址| 你懂的网址亚洲精品在线观看| 国产熟女欧美一区二区| 国产亚洲一区二区精品| 国产精品免费视频内射| 97人妻天天添夜夜摸| 青春草国产在线视频| 欧美日韩精品网址| 亚洲成人手机| 热re99久久国产66热| 亚洲国产欧美网| 黄色怎么调成土黄色| 国产高清国产精品国产三级| 啦啦啦在线观看免费高清www| 观看美女的网站| 在线观看免费高清a一片| 91精品国产国语对白视频| 悠悠久久av| 另类亚洲欧美激情| 永久免费av网站大全| 亚洲伊人久久精品综合| 观看美女的网站| 久久人人爽人人片av| 国产一区有黄有色的免费视频| 深夜精品福利| 国产女主播在线喷水免费视频网站| 久久久精品国产亚洲av高清涩受| 男女边吃奶边做爰视频| 日韩中文字幕视频在线看片| 少妇人妻 视频| 久久久亚洲精品成人影院| 中文欧美无线码| 成人国产麻豆网| 精品一区二区三区四区五区乱码 | 久久狼人影院| 亚洲成色77777| 99精国产麻豆久久婷婷| 叶爱在线成人免费视频播放| 午夜福利在线免费观看网站| www.av在线官网国产| 亚洲精品美女久久久久99蜜臀 | av国产精品久久久久影院| 观看av在线不卡| 亚洲欧美成人综合另类久久久| 久久这里只有精品19| 久久热在线av| 老司机深夜福利视频在线观看 | 多毛熟女@视频| 99久国产av精品国产电影| 男人舔女人的私密视频| 免费在线观看完整版高清| 亚洲欧美精品自产自拍| 少妇人妻久久综合中文| 人人妻人人添人人爽欧美一区卜| 亚洲精品第二区| 国产成人欧美| 日本91视频免费播放| 亚洲成av片中文字幕在线观看| 老鸭窝网址在线观看| 十分钟在线观看高清视频www| 国产精品三级大全| 中国国产av一级| 毛片一级片免费看久久久久| 亚洲精品第二区| 久久人人97超碰香蕉20202| 黄色 视频免费看| 欧美精品高潮呻吟av久久| 亚洲成人国产一区在线观看 | 满18在线观看网站| 国产成人91sexporn| 久久韩国三级中文字幕| 国产精品久久久久久久久免| 国产毛片在线视频| 国产精品二区激情视频| 伦理电影大哥的女人| 久久久久久久久久久久大奶| 香蕉丝袜av| 这个男人来自地球电影免费观看 | 久久99热这里只频精品6学生| 99久国产av精品国产电影| 国产无遮挡羞羞视频在线观看| 欧美av亚洲av综合av国产av | av不卡在线播放| 精品人妻一区二区三区麻豆| 美女高潮到喷水免费观看| 中文欧美无线码| 多毛熟女@视频| 色婷婷久久久亚洲欧美| 久久婷婷青草| 久久性视频一级片| 国产精品 欧美亚洲| 亚洲自偷自拍图片 自拍| 国产日韩欧美亚洲二区| 欧美日韩亚洲综合一区二区三区_| 99热网站在线观看| 大片免费播放器 马上看| 狠狠精品人妻久久久久久综合| 热99久久久久精品小说推荐| 大香蕉久久网| 免费女性裸体啪啪无遮挡网站| 母亲3免费完整高清在线观看| 亚洲精品日韩在线中文字幕| 亚洲熟女精品中文字幕| 99久久精品国产亚洲精品| 一边亲一边摸免费视频| 久久综合国产亚洲精品| 中文字幕av电影在线播放| 大片免费播放器 马上看| 少妇 在线观看| 美女午夜性视频免费| 亚洲男人天堂网一区| 国产欧美亚洲国产| 欧美日韩综合久久久久久| 久久精品国产亚洲av高清一级| 少妇人妻久久综合中文| 欧美日韩视频高清一区二区三区二| 国产片特级美女逼逼视频| 丝袜喷水一区| 在线观看免费视频网站a站| 中文字幕人妻丝袜制服| 国产熟女午夜一区二区三区| 香蕉国产在线看| 黑人欧美特级aaaaaa片| a级毛片黄视频| 国产欧美日韩一区二区三区在线| 国产不卡av网站在线观看| 大香蕉久久成人网| 久久青草综合色| 如何舔出高潮| 成人免费观看视频高清| 免费不卡黄色视频| 午夜福利影视在线免费观看| 免费黄频网站在线观看国产| 免费女性裸体啪啪无遮挡网站| 中文字幕色久视频| 亚洲欧美中文字幕日韩二区| 男女床上黄色一级片免费看| 成人手机av| 久久韩国三级中文字幕| 国产成人系列免费观看| 制服丝袜香蕉在线| av免费观看日本| 亚洲国产欧美网| 伊人久久大香线蕉亚洲五| 黑丝袜美女国产一区| 精品午夜福利在线看| 一级黄片播放器| 国产精品国产av在线观看| 大香蕉久久成人网| 亚洲色图综合在线观看| 亚洲精品日本国产第一区| 国产欧美日韩一区二区三区在线| 国产精品久久久人人做人人爽| 在线天堂最新版资源| 观看av在线不卡| 成人三级做爰电影| 国产99久久九九免费精品| 美女福利国产在线| 下体分泌物呈黄色| 九九爱精品视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 欧美激情 高清一区二区三区| 波野结衣二区三区在线| 成人免费观看视频高清| 国产精品一区二区在线观看99| 日韩一本色道免费dvd| 一区二区av电影网| 高清av免费在线| 久久韩国三级中文字幕| 在线免费观看不下载黄p国产| 亚洲美女搞黄在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 欧美黄色片欧美黄色片| 国产男女超爽视频在线观看| 男女边摸边吃奶| 又大又黄又爽视频免费| 午夜激情久久久久久久| 久久综合国产亚洲精品| 国产成人a∨麻豆精品| 色婷婷久久久亚洲欧美| 久久久久久久国产电影| 国产淫语在线视频| 看免费成人av毛片| 欧美日韩一区二区视频在线观看视频在线| 日韩人妻精品一区2区三区| 69精品国产乱码久久久| www.精华液| 国产av国产精品国产| 老司机在亚洲福利影院| 亚洲国产欧美日韩在线播放| 精品人妻在线不人妻| 老汉色av国产亚洲站长工具| 国产精品久久久人人做人人爽| 国产成人精品福利久久| 高清黄色对白视频在线免费看| 人体艺术视频欧美日本| 色网站视频免费| 午夜影院在线不卡| 亚洲中文av在线| 国产精品偷伦视频观看了| 啦啦啦在线观看免费高清www| 午夜福利视频精品| 一二三四在线观看免费中文在| 伊人久久国产一区二区| 国产成人午夜福利电影在线观看| 肉色欧美久久久久久久蜜桃| 如日韩欧美国产精品一区二区三区| 亚洲综合色网址| 卡戴珊不雅视频在线播放| 一二三四在线观看免费中文在| 精品亚洲成国产av| 国产欧美亚洲国产| 丝袜人妻中文字幕| 日韩大片免费观看网站| 久久综合国产亚洲精品| 无限看片的www在线观看| 午夜精品国产一区二区电影| 久久精品aⅴ一区二区三区四区| 色94色欧美一区二区| 国产极品天堂在线| 亚洲av男天堂| 国产视频首页在线观看| 咕卡用的链子| 精品国产超薄肉色丝袜足j| 巨乳人妻的诱惑在线观看| 在线天堂中文资源库| 午夜福利网站1000一区二区三区| 69精品国产乱码久久久| 国产精品欧美亚洲77777| 亚洲婷婷狠狠爱综合网| 亚洲av综合色区一区| 精品少妇黑人巨大在线播放| 精品午夜福利在线看| 久久久久久免费高清国产稀缺| 亚洲精品美女久久av网站| 亚洲伊人久久精品综合| 97精品久久久久久久久久精品| 久久女婷五月综合色啪小说| svipshipincom国产片| 国产精品国产三级国产专区5o| 两个人免费观看高清视频| 永久免费av网站大全| 亚洲国产中文字幕在线视频| 日本欧美国产在线视频| 97在线人人人人妻| 女性生殖器流出的白浆| 91aial.com中文字幕在线观看| 成年女人毛片免费观看观看9 | 久久久精品94久久精品| 最近手机中文字幕大全| 熟妇人妻不卡中文字幕| 久久久欧美国产精品| 在线看a的网站| 国产精品成人在线| 亚洲成国产人片在线观看| 一区二区三区四区激情视频| 欧美日韩av久久| 涩涩av久久男人的天堂| 美女高潮到喷水免费观看| 这个男人来自地球电影免费观看 | 国产乱来视频区| 少妇精品久久久久久久| 久久久久久久久久久免费av| 久久久亚洲精品成人影院| 18禁裸乳无遮挡动漫免费视频| 国产av国产精品国产| 国产日韩欧美视频二区| 久热这里只有精品99| 午夜福利在线免费观看网站| 国产黄频视频在线观看| 少妇的丰满在线观看| 久久精品国产综合久久久| 国产精品久久久久久精品古装| 久久av网站| 亚洲七黄色美女视频| 日日撸夜夜添| 久久精品久久精品一区二区三区| 一区二区三区激情视频| 亚洲精品久久成人aⅴ小说| 成人毛片60女人毛片免费| e午夜精品久久久久久久| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产欧美日韩在线播放| 国产成人一区二区在线| 电影成人av| 国产精品免费大片| 日韩免费高清中文字幕av| 女的被弄到高潮叫床怎么办| 宅男免费午夜| 午夜福利网站1000一区二区三区| 一级毛片 在线播放| av卡一久久| 国产欧美日韩一区二区三区在线| 99re6热这里在线精品视频| 汤姆久久久久久久影院中文字幕| 国产97色在线日韩免费| 中文字幕av电影在线播放| 纯流量卡能插随身wifi吗| 亚洲,欧美,日韩| 中文欧美无线码| 久久久精品94久久精品| 亚洲精品日本国产第一区| 国产野战对白在线观看| 免费高清在线观看日韩| 国产极品天堂在线| 桃花免费在线播放| 在线观看免费高清a一片| 丝瓜视频免费看黄片| 久久天堂一区二区三区四区| 美女大奶头黄色视频| 日韩中文字幕视频在线看片| 我要看黄色一级片免费的| 国产精品麻豆人妻色哟哟久久| av有码第一页| 别揉我奶头~嗯~啊~动态视频 | 99九九在线精品视频| 精品国产超薄肉色丝袜足j| 熟女av电影| 91国产中文字幕| 又粗又硬又长又爽又黄的视频| 下体分泌物呈黄色| 一区二区日韩欧美中文字幕| 在线精品无人区一区二区三| 亚洲美女黄色视频免费看| 亚洲综合色网址| 在线 av 中文字幕| 国产野战对白在线观看| 国产有黄有色有爽视频| 十分钟在线观看高清视频www| 多毛熟女@视频| a级毛片黄视频| avwww免费| 欧美日韩亚洲综合一区二区三区_| 99久久精品国产亚洲精品| 精品久久久精品久久久| 亚洲美女搞黄在线观看| 国产欧美亚洲国产| 99热国产这里只有精品6| 搡老乐熟女国产| 在线免费观看不下载黄p国产| 精品人妻在线不人妻| 免费观看人在逋| 大码成人一级视频| 国产在视频线精品| 国产成人欧美| 精品人妻熟女毛片av久久网站| 波多野结衣一区麻豆| av在线播放精品| 飞空精品影院首页| www.av在线官网国产| 99精品久久久久人妻精品| 国精品久久久久久国模美| 精品国产超薄肉色丝袜足j| 国语对白做爰xxxⅹ性视频网站| 亚洲熟女毛片儿| 国产在线免费精品| 久久久久国产一级毛片高清牌| 一级a爱视频在线免费观看| 欧美日本中文国产一区发布| 中文字幕色久视频| 建设人人有责人人尽责人人享有的| www.av在线官网国产| 亚洲综合色网址| 啦啦啦啦在线视频资源| 看免费成人av毛片| 一本色道久久久久久精品综合| 午夜福利影视在线免费观看| 国产午夜精品一二区理论片| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜日韩欧美国产| 少妇精品久久久久久久| 亚洲精品日韩在线中文字幕| 国产免费一区二区三区四区乱码| 欧美黑人精品巨大| 最黄视频免费看| 夜夜骑夜夜射夜夜干| 欧美黄色片欧美黄色片| 一本色道久久久久久精品综合| 亚洲三区欧美一区| 国产av码专区亚洲av| www.熟女人妻精品国产| www.自偷自拍.com| 久久久久久久久久久久大奶| 亚洲国产精品一区三区| 女人高潮潮喷娇喘18禁视频| 亚洲国产欧美一区二区综合| 色播在线永久视频| 亚洲视频免费观看视频| 七月丁香在线播放| 一区二区三区激情视频| 中文字幕人妻丝袜制服| 日本av手机在线免费观看| 日韩视频在线欧美| 在线观看免费午夜福利视频| 欧美日韩精品网址| 午夜福利免费观看在线| 男女国产视频网站| 精品人妻一区二区三区麻豆| av在线观看视频网站免费| 嫩草影院入口| 国产av一区二区精品久久| 中文字幕人妻熟女乱码| 亚洲欧美精品综合一区二区三区| 伦理电影大哥的女人| 国产片内射在线| 国产爽快片一区二区三区| 1024视频免费在线观看| 日韩免费高清中文字幕av| 中文字幕人妻丝袜一区二区 | 亚洲精品国产av蜜桃| 中文字幕人妻熟女乱码| 热re99久久精品国产66热6| 中文字幕人妻丝袜制服| 精品一区二区三区av网在线观看 | 少妇的丰满在线观看| 国产乱来视频区| 大话2 男鬼变身卡| 精品少妇内射三级| 久久亚洲国产成人精品v| xxx大片免费视频| 丝袜在线中文字幕| 男女之事视频高清在线观看 | 美女脱内裤让男人舔精品视频| 激情五月婷婷亚洲| 国产在线一区二区三区精| 超色免费av| a级毛片黄视频| 国产高清不卡午夜福利| 欧美日韩视频高清一区二区三区二| 男女无遮挡免费网站观看| 人人妻人人澡人人看| 一边摸一边做爽爽视频免费| 成年美女黄网站色视频大全免费| 男女国产视频网站| av在线app专区| 亚洲美女黄色视频免费看| 国产黄色视频一区二区在线观看| 中文字幕人妻丝袜一区二区 | 亚洲图色成人| 只有这里有精品99| 亚洲国产精品成人久久小说| 最近中文字幕2019免费版| 亚洲色图综合在线观看| 成年动漫av网址| 爱豆传媒免费全集在线观看| 青春草亚洲视频在线观看| 九色亚洲精品在线播放| 亚洲精品一二三| 91aial.com中文字幕在线观看| 久久97久久精品| av线在线观看网站| 亚洲第一区二区三区不卡| 2018国产大陆天天弄谢| 不卡av一区二区三区| 久久精品aⅴ一区二区三区四区| 亚洲精品久久久久久婷婷小说| 校园人妻丝袜中文字幕| 伊人久久大香线蕉亚洲五| 老司机影院毛片| 一级毛片 在线播放| 老司机影院成人| 久久精品aⅴ一区二区三区四区| 超色免费av| 久久韩国三级中文字幕| 丝袜喷水一区| 满18在线观看网站| 国产精品一区二区在线不卡| 热re99久久精品国产66热6| 看十八女毛片水多多多| 亚洲一级一片aⅴ在线观看| 色婷婷久久久亚洲欧美| 国产一卡二卡三卡精品 | 99re6热这里在线精品视频| 日韩制服骚丝袜av| av不卡在线播放| 少妇被粗大的猛进出69影院| 又粗又硬又长又爽又黄的视频| 亚洲成人免费av在线播放| 午夜老司机福利片| 18禁动态无遮挡网站| 国产一区二区三区av在线| 国产精品熟女久久久久浪| 欧美人与善性xxx| 亚洲国产精品一区三区| 国产成人午夜福利电影在线观看| 亚洲国产中文字幕在线视频| 欧美97在线视频| 777米奇影视久久| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人91sexporn| 黑人巨大精品欧美一区二区蜜桃| 99香蕉大伊视频| 晚上一个人看的免费电影| 超色免费av| 国产日韩欧美视频二区| 日韩av在线免费看完整版不卡| 国产 精品1| 男的添女的下面高潮视频| 超碰成人久久| 欧美日韩亚洲综合一区二区三区_| 亚洲精品av麻豆狂野| 国产成人一区二区在线| 多毛熟女@视频| 考比视频在线观看| 日韩制服骚丝袜av| 亚洲伊人色综图| 亚洲av男天堂| 丁香六月天网| 国产精品亚洲av一区麻豆 | a级毛片在线看网站| 一本久久精品| 狠狠婷婷综合久久久久久88av| 最近手机中文字幕大全| 黄频高清免费视频| 99久国产av精品国产电影| 成年美女黄网站色视频大全免费| 国语对白做爰xxxⅹ性视频网站|