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

    平紋編織面板泡沫夾芯結(jié)構(gòu)修補(bǔ)后側(cè)向壓縮有限元漸進(jìn)失效分析

    2022-05-29 08:05:41王軒冮慶庸張世秋
    航空科學(xué)技術(shù) 2022年4期
    關(guān)鍵詞:修補(bǔ)復(fù)合材料

    王軒 冮慶庸 張世秋

    摘要:為研究修補(bǔ)對平紋編織面板泡沫夾芯結(jié)構(gòu)側(cè)向壓縮性能的影響,通過USDFLD材料子程序,設(shè)置Tsai-Wu準(zhǔn)則面板材料的失效判據(jù),將線性屈曲模態(tài)考慮為初始缺陷,建立平紋編織面板泡沫夾芯結(jié)構(gòu)修補(bǔ)后側(cè)向壓縮漸進(jìn)損傷失效分析模型。結(jié)果表明,所建立模型的有效性通過試驗得到了驗證;面板挖補(bǔ)傾角在1∶10~1∶22范圍內(nèi)修補(bǔ)質(zhì)量較好,且側(cè)向壓縮強(qiáng)度隨面板挖補(bǔ)傾角比例的減小而增大;附加層重疊寬度在5~15mm范圍內(nèi)修補(bǔ)質(zhì)量較好,且側(cè)向壓縮強(qiáng)度隨附加層重疊寬度的增大而減小,但變化幅度不大;對完好泡沫夾芯結(jié)構(gòu)或者母板損傷小的修補(bǔ)件,純屈曲分析方法對其最大側(cè)向壓縮載荷的預(yù)測更加準(zhǔn)確;夾芯結(jié)構(gòu)側(cè)向壓縮的失效機(jī)理為芯材首先發(fā)生失效,然后面板發(fā)生局部屈曲,當(dāng)載荷達(dá)到試件側(cè)向壓縮強(qiáng)度后,面板最終失效。

    關(guān)鍵詞:平紋編織;復(fù)合材料;修補(bǔ);泡沫夾芯;側(cè)向壓縮;漸進(jìn)失效分析

    中圖分類號:TB332文獻(xiàn)標(biāo)識碼:ADOI:10.19452/j.issn1007-5453.2022.04.011

    基金項目:航空科學(xué)基金(2018ZF67011)

    民用飛機(jī)機(jī)載雷達(dá)罩常以玻璃纖維平紋編織材料為面板,以聚甲基丙烯酰亞胺泡沫(PMI泡沫)為芯材,與傳統(tǒng)的蜂窩夾芯結(jié)構(gòu)相比,泡沫夾芯結(jié)構(gòu)具有密封性好、整體成形工藝易實現(xiàn)等優(yōu)點[1]。由于復(fù)合材料面板較薄,夾芯結(jié)構(gòu)在使用過程中易受到面外沖擊而損壞,需要修補(bǔ)來恢復(fù)其性能。由于雷達(dá)罩結(jié)構(gòu)常處于側(cè)向受壓狀態(tài)[2],故研究平紋編織面板泡沫夾芯結(jié)構(gòu)修補(bǔ)后側(cè)向壓縮性能具有重要的工程價值。

    國內(nèi)外很多學(xué)者都對復(fù)合材料修補(bǔ)結(jié)構(gòu)進(jìn)行了仿真分析研究。張鐵純等[3]總結(jié)了挖補(bǔ)修理復(fù)合材料夾芯結(jié)構(gòu)側(cè)向壓縮性能的研究現(xiàn)狀,分析了材料、工藝、構(gòu)型及環(huán)境等影響因素,指出最新的修補(bǔ)建模方法和力學(xué)試驗結(jié)果分析思路。李劍峰等[4]在對蜂窩夾芯結(jié)構(gòu)單面貼補(bǔ)彎曲性能研究中,建立了修補(bǔ)結(jié)構(gòu)的三維有限元模型,分析結(jié)果表明,隨著補(bǔ)片的直徑和厚度的增加,修補(bǔ)結(jié)構(gòu)的承載能力呈現(xiàn)出先增加后減小的趨勢。劉國春等[1]建立了泡沫夾芯修補(bǔ)結(jié)構(gòu)有限元分析模型,結(jié)合夾芯結(jié)構(gòu)的穩(wěn)定性理論,認(rèn)為復(fù)合材料夾芯修補(bǔ)結(jié)構(gòu)的主要側(cè)壓破壞模式為面板一階與二階屈曲失效。T. Zhang等[5]對蜂窩復(fù)合材料的開孔損傷與修補(bǔ)后側(cè)向壓縮進(jìn)行了準(zhǔn)靜態(tài)模擬,考慮了層內(nèi)損傷、層間分層和蜂窩狀損傷,研究了不同挖補(bǔ)傾角比例、不同鋪層順序和不同附加層數(shù)量對挖補(bǔ)夾芯板的極限強(qiáng)度和應(yīng)力分布的影響。王躍全等[6]建立了復(fù)合材料層合板貼補(bǔ)結(jié)構(gòu)的漸進(jìn)損傷三維有限元模型,發(fā)現(xiàn)相同的補(bǔ)片搭接寬度,損傷面積越大的層合板修補(bǔ)強(qiáng)度越低。R.D.S.G.Campilho等[7]利用ABAQUS中的cohesive單元模擬膠層的力學(xué)行為建立了有限元模型,對斜接修補(bǔ)的膠結(jié)接頭進(jìn)行純壓縮屈曲驗證研究,試驗結(jié)果證明模型是有效的。鄧健等[8]基于非線性有限元方法,研究了雙面貼補(bǔ)復(fù)合材料層合板的壓縮屈曲強(qiáng)度和損傷演化過程,發(fā)現(xiàn)補(bǔ)片尺寸較大時,膠結(jié)修補(bǔ)區(qū)域的外部產(chǎn)生局部屈曲現(xiàn)象。白瑞祥等[9]利用有限元軟件ANSYS建立有限元模型,對受載裂紋損傷結(jié)構(gòu)的復(fù)合材料膠結(jié)修補(bǔ)效果進(jìn)行了斷裂力學(xué)分析,發(fā)現(xiàn)增加補(bǔ)片的剛度和增加膠層的彈性模量均會降低裂紋尖端的應(yīng)力強(qiáng)度因子。程起有等[10]開展了復(fù)合材料層合板膠結(jié)修補(bǔ)參數(shù)優(yōu)化研究,采用改進(jìn)的BP神經(jīng)網(wǎng)絡(luò)建立參數(shù)優(yōu)化預(yù)測模型,試驗結(jié)果表明模型有較高準(zhǔn)確性。呂勝利等[11]對復(fù)合材料單面膠結(jié)修補(bǔ)結(jié)構(gòu)建立了分析模型,利用該模型研究補(bǔ)片參數(shù)對膠結(jié)修補(bǔ)效果的影響,得到優(yōu)化后的補(bǔ)片參數(shù)。

    在夾芯結(jié)構(gòu)的壓縮性能仿真研究方面,王寶芹等[12]基于改進(jìn)通用單胞模型預(yù)測蜂窩夾芯結(jié)構(gòu)的壓縮失穩(wěn)載荷和損傷破壞載荷,并發(fā)現(xiàn)蜂窩夾芯結(jié)構(gòu)破壞模式為先發(fā)生屈曲失穩(wěn),然后迅速破壞。劉智良等[13]對含有分層缺陷的蜂窩夾芯結(jié)構(gòu)建立了子層局部屈曲模型,結(jié)果表明當(dāng)分層形狀及其大小不同時,夾芯板表現(xiàn)出不同的破壞機(jī)制。萬玉敏等[14]通過有限元方法對蜂窩夾層復(fù)合材料結(jié)構(gòu)的面內(nèi)壓縮整體屈曲開展研究,分析了典型薄壁夾層結(jié)構(gòu)的承載能力和破壞模式。陳悅等[15]基于非線性屈曲理論,利用三維內(nèi)聚力單元模擬面芯分層并引入初始預(yù)變形及材料損傷判據(jù),發(fā)現(xiàn)結(jié)構(gòu)首先出現(xiàn)的破壞模式為面芯脫膠分層破壞,與試驗現(xiàn)象吻合,軸壓極限承載的試驗和仿真結(jié)果最大誤差為9.04%,證明了該方法預(yù)報極限承載的可靠性。L. Fagerberg等[16]在對夾芯結(jié)構(gòu)的壓縮研究中,發(fā)現(xiàn)由于由初始缺陷引起的變形會產(chǎn)生較大的局部應(yīng)變,引發(fā)材料失效,這是強(qiáng)度問題,而不是穩(wěn)定性問題,所以建議在夾芯結(jié)構(gòu)的壓縮強(qiáng)度計算中考慮初始缺陷。曹景斌等[17]研究了濕熱條件對于層壓結(jié)構(gòu)壓縮強(qiáng)度和蜂窩夾層結(jié)構(gòu)側(cè)壓強(qiáng)度的影響,發(fā)現(xiàn)濕熱環(huán)境是造成結(jié)構(gòu)承載能力下降的重要因素。綜合上述文獻(xiàn)發(fā)現(xiàn)復(fù)合材料層合板修補(bǔ)的仿真研究方面已經(jīng)有了許多成果,蜂窩夾芯結(jié)構(gòu)較早用于航空航天領(lǐng)域,其力學(xué)性能研究成果也有不少。但關(guān)于泡沫夾芯結(jié)構(gòu)的壓縮性能研究尚少,因此本文通過建立修補(bǔ)后的平紋編織面板泡沫夾芯結(jié)構(gòu)的漸進(jìn)損傷失效分析模型,對其側(cè)向壓縮性能進(jìn)行研究。

    首先建立了基于Tsai-Wu失效準(zhǔn)則的平紋編織面板泡沫夾芯結(jié)構(gòu)漸進(jìn)損傷有限元分析模型,該模型通過芯材與面板共節(jié)點建模的方法,實現(xiàn)結(jié)構(gòu)的連續(xù)損傷,通過屈曲分析引入模型線性屈曲模態(tài)的位移,乘以一個比例因子(夾芯板厚度1%)作為模型的初始缺陷,再利用ABAQUS隱式求解器中的靜態(tài)Riks方法進(jìn)行側(cè)向壓縮強(qiáng)度分析,研究修補(bǔ)參數(shù)對平紋編織面板泡沫夾芯結(jié)構(gòu)側(cè)向壓縮強(qiáng)度和破壞模式的影響,并與試驗結(jié)果對比,驗證模型有效性。

    1試驗

    側(cè)向壓縮試驗件面板的材料是平紋編織玻璃纖維3218/SW280A,材料參數(shù)見表1,表中的E11、E22和E33為彈性模量,G12、G13和G23為切變模量,ν12、ν13和ν23為泊松比,S12、S13和S23為剪切強(qiáng)度,XT為1方向拉伸強(qiáng)度,XC為1方向壓縮強(qiáng)度,YT為2方向拉伸強(qiáng)度,YC為2方向壓縮強(qiáng)度,ZT為3方向拉伸強(qiáng)度,ZC為3方向壓縮強(qiáng)度。面板為三層平紋編織預(yù)浸料固化制成,厚度0.75mm。芯材的材料為PMI泡沫Cascell?WH,厚度為5mm,廠家提供的材料參數(shù)見表2。表2中,ρ為泡沫密度,E為泡沫彈性模量,XPC為泡沫壓縮強(qiáng)度。平紋編織面板泡沫夾芯完好試驗件的長寬高分別為150mm、100mm和6.5mm,面板由三層單板組成,每層0.25mm,尺寸如圖1所示,單側(cè)面板修補(bǔ)后試驗件的尺寸如圖2所示,附加層厚度為0.25mm,圖中的?d1,?d2,?d3分別為第1~3層圓形補(bǔ)片的直徑,?d0為圓形附加層的直徑。

    為了清楚地顯示修補(bǔ)參數(shù)和修補(bǔ)結(jié)構(gòu),以單側(cè)面板加芯材修補(bǔ)試驗件為例,在中部橫截面處繪制修補(bǔ)示意圖,如圖3所示。圖中,d為附加層重疊寬度,h為泡沫夾芯結(jié)構(gòu)單側(cè)面板厚度,也是挖補(bǔ)的厚度,α為面板挖補(bǔ)傾角,本文中的挖補(bǔ)傾角比例的值為h/L。雙側(cè)面板加芯材修補(bǔ)試驗件中部橫截面修補(bǔ)示意圖如圖4所示。各組試驗件編號對應(yīng)的修補(bǔ)方式及修補(bǔ)參數(shù)見表3,每個組有6個試驗件。

    根據(jù)美國材料與試驗協(xié)會(ASTM)D7137標(biāo)準(zhǔn)進(jìn)行側(cè)向壓縮試驗,壓縮試驗在Instron 5892型萬能材料試驗機(jī)上進(jìn)行,試驗機(jī)壓頭加載速度為1.25mm/min,具體試驗條件如圖5所示。萬能材料試驗機(jī)內(nèi)置計算記錄軟件,試驗結(jié)束后可以獲得壓頭在壓縮過程中所受力與位移曲線、壓縮強(qiáng)度和應(yīng)力應(yīng)變曲線等相關(guān)記錄。

    2漸進(jìn)損傷失效模型

    對試驗現(xiàn)象觀察可發(fā)現(xiàn),在泡沫夾芯結(jié)構(gòu)修補(bǔ)后側(cè)向壓縮過程中,結(jié)構(gòu)會出現(xiàn)面板損傷和芯材損傷。故選取連續(xù)介質(zhì)損傷力學(xué)模型和可壓碎泡沫(crushable foam)模型分別模擬復(fù)合材料面板和芯材的損傷。

    2.1面板損傷模型

    林國偉等[18]在復(fù)合材料后屈曲分析失效準(zhǔn)則適應(yīng)性研究中,發(fā)現(xiàn)不區(qū)分失效模式的Tsai-Wu失效準(zhǔn)則在計算結(jié)果精度上比區(qū)分失效模式的Hashin準(zhǔn)則、Puck準(zhǔn)則和ChangChang準(zhǔn)則要稍高,側(cè)向壓縮仿真的主要目的是預(yù)測較為精確的側(cè)向壓縮強(qiáng)度,故采用Tsai-Wu失效準(zhǔn)則作為面板材料的失效準(zhǔn)則,Tsai-Wu失效準(zhǔn)則具體形式如式(1)和式(2)所示[19]。復(fù)合材料面板性能退化具體參數(shù)見表4[20]。

    3有限元建模

    對試驗件邊界條件與加載方式進(jìn)行簡化處理,模型右側(cè)面約束所有6個自由度方向,設(shè)置為固支,對模型上下兩個側(cè)面約束Z方向的面外位移,在模型左側(cè)面選取中部一個節(jié)點作為參考點,將參考點與左側(cè)面上的所有節(jié)點的相互作用設(shè)置為梁類型的多點約束,便于在參考點施加集中力或位移載荷,約束參考點除X方向外所有5個自由度方向,有限元的模型邊界條件示意圖如圖6所示。

    根據(jù)前期已發(fā)表論文試驗結(jié)果發(fā)現(xiàn)[23],不論是完好試件還是修補(bǔ)試件,在側(cè)向壓縮試驗過程中均沒有發(fā)現(xiàn)補(bǔ)片脫黏、面芯脫膠、面板分層等膠層破壞的現(xiàn)象,說明膠層在此模型中并沒有被破壞。徐緋等[24]在蜂窩夾芯結(jié)構(gòu)側(cè)向壓縮行為研究中,發(fā)現(xiàn)膠膜的強(qiáng)度、剛度改變后,不改變蜂窩夾層結(jié)構(gòu)的破壞位置和破壞形式,結(jié)構(gòu)承載能力的變化最大不超過5%,認(rèn)為膠層參數(shù)對結(jié)構(gòu)壓縮載荷的承載能力影響很小。因此,選用拋棄膠層Cohesive單元建模的方法,將膠層對各組分的良好黏結(jié)作用,以各實體單元層之間的共節(jié)點方法來實現(xiàn),這既實現(xiàn)了膠層的黏結(jié)作用,保證了泡沫夾芯結(jié)構(gòu)損傷的連續(xù)性,也免去了引入膠層Cohesive單元帶來的模型復(fù)雜化。

    將面板和PMI芯材的單元類型設(shè)置為實體,三維8節(jié)點縮減積分單元(C3D8R)。采用這種單元時,單元形狀對縮減積分單元的計算精度影響不大,在積分點上的應(yīng)力結(jié)果相對準(zhǔn)確,這對于維修區(qū)域不規(guī)則單元形狀的有限元計算十分有利。完好試驗件的有限元模型如圖7所示。根據(jù)實際修補(bǔ)方式和修補(bǔ)參數(shù)建好的典型修補(bǔ)試件的有限元模型如圖8所示。

    復(fù)合材料層合板的強(qiáng)度計算常采用單純靜強(qiáng)度漸進(jìn)損傷分析方法,但泡沫夾芯復(fù)合材料有其獨(dú)特結(jié)構(gòu),若采用單純靜強(qiáng)度漸進(jìn)損傷分析會與真實情況差距很大。為了解決這一問題,本文通過引入線性屈曲模態(tài)來考慮初始缺陷,即利用弧長法,將夾芯結(jié)構(gòu)的初始缺陷考慮為線性屈曲模態(tài),用后屈曲強(qiáng)度分析Riks算法的結(jié)果來代表模型的側(cè)向壓縮強(qiáng)度。

    根據(jù)上述模型,利用ABAQUS有限元商用軟件的腳本接口,編寫關(guān)于面板的材料參數(shù)剛度退化腳本程序,以實現(xiàn)面板材料的剛度退化。用ABAQUS軟件USDFLD子程序格式,將含有漸進(jìn)損傷分析的Tsai-Wu失效準(zhǔn)則用FORTRAN語言編寫出來,以實現(xiàn)復(fù)合材料面板在Tsai-Wu失效準(zhǔn)則下的漸進(jìn)損傷過程。

    有限元漸進(jìn)失效分析流程如圖9所示。建立有限元模型并進(jìn)行線性屈曲分析,得到泡沫夾芯結(jié)構(gòu)的各階屈曲模態(tài),引入泡沫夾芯結(jié)構(gòu)的前5階屈曲模態(tài)作為結(jié)構(gòu)的初始狀態(tài),取初始缺陷破壞因子為夾芯厚度1%,即0.065mm,采用線性逐步加載的方式對結(jié)構(gòu)進(jìn)行強(qiáng)度分析,求解每個增量步中結(jié)構(gòu)的應(yīng)力應(yīng)變。由Tsai-Wu損傷準(zhǔn)則判斷面板材料是否失效,并對失效后單元的材料參數(shù)進(jìn)行剛度退化處理,在下一增量步中對剛度退化后的結(jié)構(gòu)再進(jìn)行應(yīng)力應(yīng)變分析,并對單元的材料失效進(jìn)行判斷和剛度退化處理,循環(huán)往復(fù),直至結(jié)構(gòu)完全失效無法迭代計算并停止加載。

    4結(jié)果分析與討論

    4.1側(cè)向壓縮強(qiáng)度

    4.1.1考慮初始缺陷

    在正式開始分析前,先對比引入線性屈曲模態(tài)和不引入屈曲模態(tài)方法對強(qiáng)度計算的不同,并分析其與實際試驗結(jié)果的吻合度。以完好試件組為例,分別利用這兩種計算方法,得出計算結(jié)果和與試驗時的誤差見表5??梢娨刖€性屈曲模態(tài)的計算方法在夾芯結(jié)構(gòu)的側(cè)向壓縮強(qiáng)度計算上與實際試驗結(jié)果更加符合。

    分別對試驗件建立有限元模型并進(jìn)行分析,將側(cè)向壓縮試驗的試驗結(jié)果與有限元模型的仿真結(jié)果進(jìn)行對比。完好及各修補(bǔ)方式和修補(bǔ)參數(shù)下有限元模型與對應(yīng)的試驗件側(cè)向壓縮強(qiáng)度數(shù)值見表6。

    從表6中可以看出,除xb13組和xb33組外,有限元仿真計算結(jié)果在各組上的計算是相當(dāng)準(zhǔn)確的,其他修補(bǔ)模型有限元仿真值與試驗均值誤差絕對值不超過9%,可見有限元仿真模型的有效性。然而xb13組,也就是雙側(cè)面板加芯材修補(bǔ)試驗件的仿真計算結(jié)果與試驗均值相差27.10%。研究認(rèn)為,該組試驗件有兩個面板要修補(bǔ),采用了二次固化的方式,且修補(bǔ)為手工操作,修補(bǔ)涉及到易碎裂的低密度PMI泡沫芯材的替換,存在較多的內(nèi)部缺陷引入可能,因此導(dǎo)致其該組試驗件均值偏低,與仿真值偏離較遠(yuǎn)。xb33組模型的仿真值與試驗值相比誤差較大,分析后認(rèn)為附加層重疊寬度為25mm時,附加層的直徑已經(jīng)達(dá)到98mm,尺寸接近試件的寬度,附加層容易受到壓縮夾具的影響。而在仿真計算條件下,對模型邊界條件做了簡化處理,不會影響附加層的受力狀況,附加層直徑尺寸過大也增加了引入膠結(jié)缺陷的可能性,而在仿真分析中,用共節(jié)點方法設(shè)置膠結(jié)的質(zhì)量良好,不會因為膠結(jié)缺陷導(dǎo)致試驗件側(cè)向壓縮強(qiáng)度降低,正是這些原因使得附加層重疊寬度為25mm試驗件(xb33)側(cè)向壓縮強(qiáng)度與仿真值偏差較大。

    從表6的數(shù)據(jù)中也可以看到,無論是仿真計算數(shù)值還是試驗結(jié)果數(shù)值,均顯示單側(cè)面板修補(bǔ)泡沫夾芯結(jié)構(gòu)的側(cè)向壓縮強(qiáng)度高于單側(cè)面板加芯材修補(bǔ)的形式,隨面板挖補(bǔ)傾角的比例減小,修補(bǔ)試驗件側(cè)向壓縮強(qiáng)度也減小,當(dāng)附加層重疊寬度過大時,修補(bǔ)試驗件側(cè)向壓縮強(qiáng)度減小。

    4.1.2載荷位移曲線

    通過對試驗件建立仿真模型,利用模型左側(cè)面上參考點得到夾芯結(jié)構(gòu)壓縮方向位移載荷曲線。由于試驗件修補(bǔ)前后采用的泡沫夾芯結(jié)構(gòu)的材料與組成完全相同,因此修補(bǔ)后試驗件側(cè)向壓縮過程的位移載荷曲線形式基本相同,這里取修補(bǔ)類xb11中的一個典型試驗件與有限元模型側(cè)向壓縮位移載荷曲線進(jìn)行對比,如圖10所示。

    有限元模型與試驗件的位移載荷曲線形式基本相同,在未達(dá)到側(cè)向壓縮強(qiáng)度前,兩者的位移與載荷基本上呈線性變化,達(dá)到側(cè)向壓縮強(qiáng)度后載荷值突然下降,試件失效。在位移加載初始階段,試驗件的側(cè)向壓縮載荷相比仿真模型增加較慢,這是由于試驗件與夾具之間存在微小縫隙,隨位移增加到0.32mm附近時,試驗件與夾具間的縫隙完全消失,試驗件側(cè)向壓縮載荷隨位移的增加也近似呈線性,與仿真模型的載荷變化形式相同,且最終達(dá)到基本相等的側(cè)向壓縮強(qiáng)度。為此,將仿真的位移載荷曲線向右沿位移坐標(biāo)軸平移0.32mm,如圖10中所示,平移后兩者在側(cè)向壓縮載荷和失效位移上數(shù)值基本相同,這驗證了上文的說法。

    圖10中試驗件側(cè)向壓縮載荷突變后穩(wěn)定在5kN左右,而有限元模型側(cè)向壓縮載荷突變后穩(wěn)定在10kN左右,這是因為試驗中PMI泡沫在達(dá)到壓縮強(qiáng)度后就會碎裂,幾乎不再給面板提供任何支撐,而仿真模型中泡沫在達(dá)到強(qiáng)度極限后,會保持它達(dá)到側(cè)向壓縮強(qiáng)度時能提供的支撐力,這個特性是由ABAQUS中可壓碎泡沫模型決定的。

    4.1.3面板挖補(bǔ)傾角比例的影響

    在上文對有限元模型驗證時已經(jīng)有三個不同面板挖補(bǔ)傾角1∶10、1∶20、1∶30,因此僅在有限元模型中增加了1∶5、1∶15、1∶22、1∶25這4種面板挖補(bǔ)傾角比例模型進(jìn)行分析。因為1∶22這種挖補(bǔ)傾角比例在本次挖補(bǔ)修補(bǔ)試驗中使用,故也將其納入為有限元分析模型的一個參數(shù)值。

    分析時又加入對各階屈曲模態(tài)下側(cè)向壓縮強(qiáng)度的研究,因為在仿真研究中以模型各階屈曲模態(tài)節(jié)點的位移,乘以一個為試驗件厚度1/100的缺陷因子來作為模型的初始狀態(tài),通過這個初始狀態(tài)來模擬試驗件存在的各種微小形變和側(cè)向壓縮試驗初期試驗件的形變,每一階屈曲模態(tài)的位移變形都代表了試驗件可能發(fā)生的一種微小形變?nèi)毕?。通過比較同一模型各階屈曲模態(tài)下側(cè)向壓縮強(qiáng)度差異,就能知道修補(bǔ)后的結(jié)構(gòu)在稍有不同的初始狀態(tài)下側(cè)向壓縮強(qiáng)度的變化范圍,對修補(bǔ)后強(qiáng)度變化范圍大小進(jìn)行分析,變化范圍越小,說明結(jié)構(gòu)在各種微小缺陷的情況下都處于穩(wěn)定的狀態(tài),修補(bǔ)質(zhì)量較好。

    不同面板挖補(bǔ)傾角模型各階屈曲模態(tài)下側(cè)向壓縮強(qiáng)度的數(shù)據(jù)如圖11所示,圖中的多階仿真組所取初始缺陷破壞因子為夾芯厚度1%,即0.065mm,因為和前三階屈曲模態(tài)下側(cè)向壓縮強(qiáng)度進(jìn)行對比,混合模態(tài)模型節(jié)點位移變形僅取前三階。

    從圖11中對仿真(多階)組側(cè)向壓縮強(qiáng)度載荷的變化看,模型側(cè)向壓縮強(qiáng)度隨面板挖補(bǔ)傾角比例的減小,呈現(xiàn)出先增大后減小,而后再增大并最終減小的趨勢。面板挖補(bǔ)傾角比例為1∶10的模型側(cè)向壓縮強(qiáng)度較大,為先出現(xiàn)的峰值位置,面板挖補(bǔ)傾角為1∶25的模型側(cè)向壓縮強(qiáng)度也較大,為出現(xiàn)的第二個峰值位置。

    在圖11中,從模型各階屈曲模態(tài)下的側(cè)向壓縮強(qiáng)度變化范圍大小來看,面板挖補(bǔ)傾角比例大于1:22時,模型的各階模態(tài)下側(cè)向壓縮強(qiáng)度變化范圍較小,修補(bǔ)質(zhì)量較好;而面板挖補(bǔ)傾角比例小于1∶22,在圖11中即為1∶25和1∶30時,模型的各階模量下側(cè)向壓縮強(qiáng)度變化范圍較大,修補(bǔ)后質(zhì)量不佳。

    綜合試驗結(jié)果與上述有限元仿真分析計算認(rèn)為,泡沫夾芯結(jié)構(gòu)在單側(cè)面板修補(bǔ)時,試驗件面板挖補(bǔ)傾角比例為1∶10~1∶22時,修補(bǔ)質(zhì)量較好,推薦使用。在這個范圍內(nèi)修補(bǔ)后試驗件側(cè)向壓縮強(qiáng)度隨面板挖補(bǔ)傾角比例的減小而減小。

    4.1.4附加層重疊寬度的影響

    劉國春等[25]認(rèn)為在修補(bǔ)平紋編織復(fù)合材料拉伸試件時,存在使得試驗件強(qiáng)度取得最大值的附加層重疊寬度。Beylergil等[26]在對階梯狀接頭的軸壓研究中,發(fā)現(xiàn)隨附加層重疊寬度增加,試驗件破壞載荷會變大。為進(jìn)一步分析附加層重疊寬度對側(cè)向壓縮強(qiáng)度的影響,將附加層重疊寬度修補(bǔ)參數(shù)的范圍由原模型中的15mm、20mm、25mm這三個進(jìn)行擴(kuò)充,增加附加層重疊寬度參數(shù)為0mm、2.5mm、5mm、7.5mm、10mm、12.5mm、17.5mm、22.5mm的模型,對各修補(bǔ)模型進(jìn)行側(cè)向壓縮有限元仿真,計算結(jié)果如圖12所示,仿真多階組模型的混合模態(tài)節(jié)點位移僅取前三階。

    從圖12中仿真(多階)有限元模型側(cè)向壓縮強(qiáng)度的變化可以看出,隨著附加層重疊寬度的增加,修補(bǔ)模型的側(cè)向壓縮強(qiáng)度呈現(xiàn)出波浪形的變化趨勢,附加層重疊寬度為2.5mm處為一個波形底部,附加層重疊寬度為17.5mm處為一個波形頂部,在它們之間附加層重疊寬度為5~15mm的區(qū)域內(nèi)側(cè)向壓縮強(qiáng)度值變化較平緩,而附加層重疊寬度為17.5mm模型的側(cè)向壓縮強(qiáng)度為0~25mm附加層重疊范圍中出現(xiàn)的最大峰值。

    在圖12中,從模型各階屈曲模態(tài)下的側(cè)向壓縮強(qiáng)度變化范圍大小來看,單側(cè)面板修補(bǔ)中,附加層重疊寬度為2.5~ 15mm時,各階屈曲模態(tài)下壓縮強(qiáng)度差異較小,修補(bǔ)質(zhì)量較好。附加層重疊寬度為0mm在實際修補(bǔ)中無意義,僅作對照。附加層重疊寬度為2.5mm模型的各階模態(tài)下側(cè)向壓縮強(qiáng)度變化很小,但其側(cè)向壓縮強(qiáng)度在所有模型中最低,不推薦采用。附加層重疊寬度為17.5~22.5mm時,模型各階模態(tài)下的側(cè)向壓縮強(qiáng)度差異很大,修補(bǔ)質(zhì)量不佳。附加層重疊寬度為25mm時,模型各階模態(tài)下側(cè)向壓縮強(qiáng)度變化范圍小,其側(cè)向壓縮強(qiáng)度同附加層重疊寬度為5~15mm時基本相同,卻增加修補(bǔ)的材料和修補(bǔ)質(zhì)量,這得不償失。

    綜合試驗結(jié)果和上述有限元仿真分析計算,認(rèn)為在泡沫夾芯結(jié)構(gòu)單側(cè)面板修補(bǔ)時,附加層重疊寬度為5~15mm時修補(bǔ)質(zhì)量較好,推薦使用。在附加層重疊寬度為5~ 10mm和12.5~15mm范圍中,修補(bǔ)后試驗件側(cè)向壓縮強(qiáng)度分別隨附加層重疊寬度的增大而減小,但變化幅度不大。雖然模型附加層重疊寬度在5~15mm范圍內(nèi)5mm處的側(cè)向壓縮強(qiáng)度取得最大值,但考慮到實際修補(bǔ)時可能會有誤差,從能夠獲得可靠的修補(bǔ)效果考慮,推薦盡可能采用附加層重疊寬度為12.5mm的修補(bǔ)方案。

    4.2損傷模式

    試驗后各組面板的典型損傷形式如圖13所示,試驗件面板損傷處用矩形框進(jìn)行示意。修補(bǔ)件的損傷主要在附加層區(qū)域的外圍邊界,如圖13所示,完好試件的損傷在面板的左右兩側(cè)和中部均有出現(xiàn)。圖13中展示了完好試件類別(xb00)中的一個損傷在中部偏左側(cè)的試件與其他修補(bǔ)件進(jìn)行對比,修補(bǔ)件的中部區(qū)域因附加層的加強(qiáng)作用,沒有出現(xiàn)中部損傷的情況,這說明修補(bǔ)的有效性。有限元模型側(cè)向壓縮仿真模擬后的面板的損傷情況如圖14所示,圖中失效單元通過顏色標(biāo)記高亮展示,仿真結(jié)果中面板的失效位置與試驗中基本一致。

    4.3失效機(jī)理分析

    為了分析試驗件達(dá)到側(cè)向壓縮強(qiáng)度時的失效原因,注意到在位移載荷曲線圖10上有轉(zhuǎn)折點A、B兩點,A點是有限元模型發(fā)生結(jié)構(gòu)失效的起點,B點是有限元模型的受力達(dá)到了新的平衡的起點,弄清楚了從A點到B點有限元模型各部分應(yīng)力應(yīng)變是如何變化,有哪些部分發(fā)生了損傷或損傷擴(kuò)展,也就知道了泡沫夾芯結(jié)構(gòu)的失效機(jī)理。

    在第一個轉(zhuǎn)折點A處,通過對有限元模型整體變形應(yīng)力云圖觀察,發(fā)現(xiàn)此時平紋編織復(fù)合材料面板遠(yuǎn)遠(yuǎn)未達(dá)到其應(yīng)力強(qiáng)度極限,面板未損傷,應(yīng)力云圖如圖15所示。但通過觀察PMI泡沫芯材有限元模型變形應(yīng)力云圖,發(fā)現(xiàn)芯材模型一部分的應(yīng)力已經(jīng)達(dá)到材料壓縮強(qiáng)度,芯材發(fā)生了損傷,芯材模型的變形應(yīng)力云圖如圖16所示。因此,認(rèn)為泡沫夾芯結(jié)構(gòu)達(dá)到其側(cè)向壓縮強(qiáng)度,不是由于面板的破壞,而是強(qiáng)度較弱芯材破壞導(dǎo)致的。

    繼續(xù)對圖10曲線中轉(zhuǎn)折點B處模型的變形應(yīng)力狀態(tài)進(jìn)行研究,觀察此時有限元模型的整體變形應(yīng)力云圖,如圖17所示。從圖10中知夾芯結(jié)構(gòu)由A狀態(tài)到B狀態(tài)壓縮方向位移幾乎沒有增加,但芯材破壞處面板的應(yīng)力急劇增加,從泡沫芯材模型的變形應(yīng)力云圖18與圖16的對比中能看出,芯材應(yīng)力達(dá)到其壓縮強(qiáng)度的部分變多,芯材損傷面積擴(kuò)大,這一點則說明了面板的應(yīng)力突變與芯材損傷擴(kuò)展有關(guān)。PMI泡沫芯材的彈性模量與平紋編織面板相比很小,結(jié)構(gòu)的承載主要靠面板材料,而從此時的圖17來看,面板并未達(dá)到材料強(qiáng)度極限,也就是說此時面板材料并沒有失效,但其承載載荷的能力卻隨著芯材失效面積的擴(kuò)大而急劇降低,這印證了上文所說的泡沫夾芯結(jié)構(gòu)達(dá)到側(cè)向壓縮強(qiáng)度極限是泡沫芯材的破壞導(dǎo)致的。

    側(cè)向壓縮試驗結(jié)束后,試驗件面板確實出現(xiàn)了損傷,而在有限元模型側(cè)向壓縮過程中,直至側(cè)向壓縮載荷突降后達(dá)到B點,有限元模型再次受力平衡,模型面板也未出現(xiàn)損傷。推測模型面板的損傷發(fā)生在后續(xù)的側(cè)向壓縮過程,注意到圖10中仿真的位移載荷曲線B點后還有一個轉(zhuǎn)折點C,觀察該處有限元模型面板變形應(yīng)力云圖情況,發(fā)現(xiàn)有限元模型面板確實在此時出現(xiàn)損傷,將損傷處面板的失效單元標(biāo)記出來,并在面板上下兩側(cè)邊將破壞處框選放大,如圖19所示。

    在側(cè)向壓縮過程的最后,面板發(fā)生了失效,根據(jù)上文知泡沫夾芯結(jié)構(gòu)的側(cè)向壓縮強(qiáng)度的降低與芯材的損傷有關(guān),然而芯材的損傷如何影響夾芯結(jié)構(gòu)的面板,進(jìn)而使整個泡沫夾芯結(jié)構(gòu)的側(cè)向承載能力下降,并最終使得面板損傷,這一點可以通過觀察模型位移載荷曲線A點、B點和C點處模型變形的側(cè)向視圖來研究,有限元模型在各點處的變形側(cè)向視圖如圖20所示,從左到右依次為A點、B點和C點處模型變形的側(cè)向視圖。

    A點處面板未出現(xiàn)明顯面外位移變形,但泡沫芯材開始出現(xiàn)損傷,B點處面板面外位移變形明顯,泡沫芯材損傷面積擴(kuò)大,同時從圖17中可知面板應(yīng)力相對于A點有了突增,C點處面板的面外位移變形加大,此時面板材料出現(xiàn)損傷??梢钥闯觯静膿p傷導(dǎo)致面板局部無支撐,在側(cè)向壓縮載荷的作用下面板發(fā)生了局部屈曲,并使得泡沫夾芯結(jié)構(gòu)失效,整個結(jié)構(gòu)的側(cè)向壓縮承載能力突降,隨著位移加載的繼續(xù),面板變形過大,最終發(fā)生材料的損傷。

    4.4純屈曲分析

    從表7中看到試驗件側(cè)向壓縮時能夠承受的最大載荷值均小于仿真模型一階屈曲載荷,這說明試驗件未達(dá)到屈曲狀態(tài)就已經(jīng)失效,通過上文分析可知,試驗件未達(dá)到屈曲狀態(tài)的原因是泡沫芯材先達(dá)到了壓縮強(qiáng)度,發(fā)生了失效。根據(jù)試驗件應(yīng)變計的記錄結(jié)果顯示,試驗件在側(cè)向壓縮過程中沒有發(fā)生整體屈曲破壞,但觀察到了面板在芯材失效后,發(fā)生局部屈曲的試驗現(xiàn)象。因此認(rèn)為試驗件失效機(jī)理為芯材失效后,面板局部無支撐,發(fā)生了局部屈曲。

    表7中多數(shù)模型的一階屈曲載荷值與試驗中真實的最大側(cè)向壓縮載荷有很大差距,但也觀察到xb00組和xb21組的試驗承載能力與一階屈曲載荷值比較接近,xb00組是完好試件,xb11組是修補(bǔ)時對母板去除最小的試件,這都說明一階屈曲載荷對完好試驗件或接近完好試驗件的最大側(cè)向壓縮載荷能提供一個預(yù)測值,這個預(yù)測值稍微偏大,但誤差分別僅為3.14%和6.15%,準(zhǔn)確率很高。

    通過上文分析可知,試驗件到達(dá)屈曲前泡沫芯材就已經(jīng)失效,而屈曲特征值分析得到的屈曲載荷值只考慮材料的彈性模量和邊界條件,不考慮材料的失效,因此與試驗得到的最大側(cè)向壓縮載荷存在誤差。若芯材的強(qiáng)度變大,足以支撐試驗件達(dá)到屈曲狀態(tài),用仿真模型進(jìn)行一階屈曲載荷值計算,來預(yù)測試驗件的側(cè)向壓縮強(qiáng)度會更加準(zhǔn)確,這一點有助于優(yōu)化泡沫夾芯結(jié)構(gòu)的強(qiáng)度設(shè)計。A.G.Mamalis等[27]也認(rèn)為泡沫夾芯結(jié)構(gòu)中芯材的性能和強(qiáng)度非常重要,是決定泡沫夾芯結(jié)構(gòu)側(cè)向壓縮失效模式的重要因素。

    5結(jié)論

    通過對平紋編織面板泡沫夾芯結(jié)構(gòu)修補(bǔ)試驗件建立漸進(jìn)損傷的有限元模型,進(jìn)行側(cè)向壓縮數(shù)值仿真,與試驗結(jié)果進(jìn)行對比分析。

    (1)從泡沫夾芯結(jié)構(gòu)修補(bǔ)后的側(cè)向壓縮強(qiáng)度、失效過程、面板損傷位置、失效機(jī)理等多方面驗證了模型的有效性。

    (2)根據(jù)有限元仿真結(jié)合試驗結(jié)果分析,認(rèn)為泡沫夾芯結(jié)構(gòu)的側(cè)向壓縮失效機(jī)理為芯材損傷導(dǎo)致面板局部無支撐,在側(cè)向壓縮載荷的作用下面板發(fā)生了局部屈曲,并使得泡沫夾芯結(jié)構(gòu)失效,達(dá)到側(cè)向壓縮強(qiáng)度。

    (3)分析面板挖補(bǔ)傾角比例對模型側(cè)向壓縮強(qiáng)度的影響,得出面板挖補(bǔ)傾角為1∶10到1∶22時,修補(bǔ)質(zhì)量較好,在這個范圍內(nèi)修補(bǔ)后,試驗件側(cè)向壓縮強(qiáng)度隨面板挖補(bǔ)傾角比例的減小而增大。分析面板附加層重疊寬度對模型側(cè)向壓縮強(qiáng)度的影響,為了獲得可靠的修補(bǔ)效果,推薦盡可能采用附加層重疊寬度為12.5mm的修補(bǔ)方案。

    (4)通過純屈曲分析與試驗結(jié)果對比,認(rèn)為對完好泡沫夾芯結(jié)構(gòu)或者母板損傷小的修補(bǔ)件,純屈曲分析預(yù)測其最大側(cè)向壓縮載荷準(zhǔn)確率高。

    參考文獻(xiàn)

    [1]劉國春,魏桂明,楊文鋒,等.復(fù)合材料泡沫夾芯板膠結(jié)修理的壓縮性能[J].宇航材料工藝, 2020, 50(2): 16-21. LiuGuochun,WeiGuiming,YangWenfeng,etal. Compression performance of adhesively repaired foam core sandwich composite[J]. Aerospace Materials & Technology, 2020, 50(2): 16-21. (in Chinese)

    [2]周春蘋,劉武帥,王軒,等.含表面劃傷平紋編織面板蜂窩夾芯結(jié)構(gòu)側(cè)向壓縮漸進(jìn)失效分析[J].玻璃鋼/復(fù)合材料, 2019(7): 5-10. Zhou Chunping, Liu Wushuai, Wang Xuan, et al. Study on edgewise compression performance of honeycomb sandwich structure with scratch damage[J]. Fiber Reinforced Plastics/ Composites, 2019(7): 5-10. (in Chinese)

    [3]張鐵純,張世秋,王軒,等.挖補(bǔ)修理復(fù)合材料夾芯結(jié)構(gòu)側(cè)向壓縮性能研究[J].航空科學(xué)技術(shù),2021,32(8):1-11. Zhang Tiechun, Zhang Shiqiu, Wang Xuan, et al. Research on lateral compression performance of scarf repaired composite sandwich structure[J]. Aeronautical Science & Technology, 2021,32(8):1-11. (in Chinese)

    [4]李劍峰,燕瑛,廖寶華,等.復(fù)合材料蜂窩夾芯結(jié)構(gòu)單面貼補(bǔ)彎曲性能的分析模型與試驗研究[J].航空學(xué)報, 2013, 34(8): 1884-1891. Li Jianfeng, Yan Ying, Liao Baohua, et al. Analysis model and experimentalstudyofbendingbehaviorofcomposite honeycomb sandwich structures with one side bonded repair[J]. Acta Aeronautica et Astronautica Sinica, 2013, 34(8): 1884-1891. (in Chinese)

    [5]Zhang T,Yan Y,Jin C. Experimental and uumerical investigations of honeycomb sandwich composite panels with open-hole damage and scarf repair subjected to compressive loads[J]. The Journal ofAdhesion,2016,92(5):380-401.

    [6]王躍全,童明波,朱書華.復(fù)合材料層合板膠結(jié)貼補(bǔ)修理漸進(jìn)損傷分析[J].復(fù)合材料學(xué)報, 2011, 3(28): 197-202. Wang Yuequan, Tong Mingbo, Zhu Shuhua. Progressive damage analysis on adhesively bonding patch repair of composite laminates[J]. Acta Materiae Compositae Sinica, 2011, 3(28): 197-202. (in Chinese)

    [7]Campilho R D S G,De Moura M F S F,Ramantani D A,et al. Buckling behaviour of carbon epoxy adhesively-bonded scarf repairs[J]. Journal of Adhesion Science and Technology,2009,23(10-11):1493-1513.

    [8]鄧健,周光明,尹喬之,等.雙面貼補(bǔ)層合板壓縮屈曲漸進(jìn)損傷分析[J].航空學(xué)報, 2016, 37(5): 1526-1535. Deng Jian, Zhou Guangming, Yin Qiaozhi, et al. Progressive damage analysis of double strap bonding repaired laminates underbucklingcompression[J].ActaAeronauticaet Astronautica Sinica, 2016, 37(5): 1526-1535. (in Chinese)

    [9]白瑞祥,王秀秀,王蔓.碳纖維復(fù)合材料修補(bǔ)結(jié)構(gòu)的ANSYS參數(shù)化建模技術(shù)與修補(bǔ)參數(shù)討論[C]//第19屆全國結(jié)構(gòu)工程學(xué)術(shù)會議.濟(jì)南,2010:321-324. Bai Ruixiang, Wang Xiuxiu, Wang Man. ANSYS parametric modeling technology and repair parameters discussion of carbon fiber composite repair structure[C]//The 19th National Academic Conference on Structural Engineering. Jinan, 2010: 321-324. (in Chinese)

    [10]程起有,童小燕,姚磊江,等.基于神經(jīng)網(wǎng)絡(luò)的復(fù)合材料膠結(jié)修補(bǔ)參數(shù)優(yōu)化[J].計算機(jī)仿真, 2008, 25(12): 205-207. Cheng Qiyou, Tong Xiaoyan, Yao Leijiang, et al. Parameter optimization for bonding repair of composite laminates based on neural network[J]. Computer Simulation, 2008, 25(12): 205-207. (in Chinese)

    [11]呂勝利,程起有,姚磊江,等.損傷復(fù)合材料層板膠結(jié)修理的優(yōu)化設(shè)計[J].機(jī)械強(qiáng)度, 2007(4): 598-600. Lyu Shengli, Cheng Qiyou, Yao Leijiang, et al. Optimization design of adhesive bonding repair of damaged composite laminates[J]. Journal of Mechanical Strength, 2007(4): 598-600. (in Chinese)

    [12]王寶芹,王沫楠,劉長喜.基于多尺度方法的蜂窩夾層復(fù)合材料結(jié)構(gòu)軸向壓縮穩(wěn)定性[J].復(fù)合材料學(xué)報, 2020, 37(3): 601-608. Wang Baoqin, Wang Monan, Liu Changxi. Stability of honeycombsandwichcompositestructureunderanxial compression based on multi-scale method[J]. Acta Materiae Compositae Sinica, 2020, 37(3): 601-608.(in Chinese)

    [13]劉智良,程小全,胡仁偉,等.含不同形狀分層缺陷蜂窩夾層板的壓縮性能[J].復(fù)合材料學(xué)報, 2009, 26(5): 153-158. Liu Zhiliang, Cheng Xiaoquan, Hu Renwei, et al. Compression properties of honeycomb sandwich plates with different shape delamination defects[J]. Acta Materiae Compositae Sinica, 2009, 26(5): 153-158. (in Chinese)

    [14]萬玉敏,張發(fā),劉長喜,等.飛機(jī)典型薄壁復(fù)合材料夾層結(jié)構(gòu)整體屈曲[J].復(fù)合材料學(xué)報, 2018, 35(8): 2235-2245. Wan Yumin, Zhang Fa, Liu Changxi, et al. Overall buckling of typical thin-wall sandwich composites applied on the aircraft[J]. Acta Materiae Compositae Sinica, 2018, 35(8): 2235-2245.(in Chinese)

    [15]陳悅,朱錫,李華東,等.復(fù)合材料夾芯梁屈曲破壞模式及極限承載[J].復(fù)合材料學(xué)報, 2016, 33(5): 991-997. Chen Yue, Zhu Xi, Li Huadong, et al. Buckling failure mode and ultimate load of composite sandwich beam[J]. Acta Materiae Compositae Sinica, 2016, 33(5): 991-997.(in Chinese)

    [16]Fagerberg L,Zenkert D. Imperfection-induced wrinkling material failure in sandwich panels[J]. Journal of Sandwich Structures & Materials,2016,7(3):195-219.

    [17]曹景斌,王松,章強(qiáng).樹脂基復(fù)合材料濕熱壓縮性能研究[J].航空科學(xué)技術(shù),2020,31(3):47-52. Cao Jingbin, Wang Song, Zhang Qiang. Research on resin composite hydrothermal compressive properties[J]. Aeronautical Science & Technology, 2020, 31(3):47-52. (in Chinese)

    [18]林國偉,李新祥.復(fù)合材料壁板后屈曲分析失效準(zhǔn)則的適用性研究[J].機(jī)械強(qiáng)度, 2020, 42(3): 565-572. Lin Guowei, Li Xinxiang. Applicability investigation of failure criterion for post-bucking analysis of stiffened composite panels[J]. Journal of Mechanical Strength, 2020, 42(3): 565-572. (in Chinese)

    [19]Tsai S W,Wu E M. A general theory of strength for anisotropic materials[J]. Journal of Composite Materials,1971,5(1):58-80.

    [20]孔祥宏,王志瑾.基于Abaqus的有限元漸進(jìn)損傷強(qiáng)度分析[J].計算機(jī)應(yīng)用與軟件, 2012, 29(10): 236-240. Kong Xianghong, Wang Zhijin. Finite element incremental damage strength analysis based on Abaqus[J]. Computer Applications and Software, 2012, 29(10): 236-240. (in Chinese)

    [21]Wang J,Chen B,Wang H,et al. Experimental study on the compression-after-impact behavior of foam-core sandwich panels[J]. Journal of Sandwich Structures & Materials,2015,17(4):446-465.

    [22]譚年富.復(fù)合材料泡沫夾層結(jié)構(gòu)低速沖擊后的面內(nèi)力學(xué)性能研究[D].上海:上海交通大學(xué), 2012. TanNianfu.Mechanicalpropertiesoffoamsandwich composites with low-velocity impact damage[D]. Shanghai: Shanghai Jiao Tong University, 2012. (in Chinese)

    [23]張鐵純,張世秋,王軒,等.平紋編織面板泡沫夾芯結(jié)構(gòu)修補(bǔ)后側(cè)向壓縮性能[J].復(fù)合材料科學(xué)與工程,2021(3):51-59. Zhang Tiechun, Zhang Shiqiu, Wang Xuan, et al. Research on edgewise compression properties of plain weave panel foam sandwich composite material after repair[J]. Composites Science and Engineering ,2021(3):51-59. (in Chinese)

    [24]徐緋,劉亞各,閆慧敏.蜂窩夾芯結(jié)構(gòu)沖擊損傷后的壓縮行為研究[J].應(yīng)用力學(xué)學(xué)報, 2013, 30(5): 726-730. Xu Fei, Liu Yage, Yan Huimin. Honeycomb sandwich compositestructurewithpre-damageunderuniaxial compression[J]. Chinese Journal of Applied Mechanics, 2013, 30(5): 726-730. (in Chinese)

    [25]劉國春,秦文峰,楊文鋒,等.平紋編織復(fù)合材料挖補(bǔ)修理附加層優(yōu)化研究[J].合肥工業(yè)大學(xué)學(xué)報(自然科學(xué)版), 2015, 38(8): 1045-1048. Liu Guochun, Qin Wenfeng, Yang Wenfeng, et al. Parameters optimization of covering layer of plain-woven composite laminates with flush repair[J]. Journal of Hefei University of Technology, 2015, 38(8): 1045-1048. (in Chinese)

    [26]BeylergilB,AktasA,CunediogluY.Bucklingand compressive failure of stepped-lap joints repaired with composite patches[J]. Journal of Composite Materials,2012,46(26):3213-3230.

    [27]Mamalis A G,Manolakos D E,Ioannidis M B,et al. On the crushing response of composite sandwich panels subjected to edgewise compression:experimental[J]. Composite Structures,2005,71(2):246-257.

    Progressive Failure Analysis of Repaired Foam Sandwich Structure with Plain Weave Faceplate Under Edgewise Compression

    Wang Xuan,Gang Qingyong,Zhang Shiqiu

    Civil Aviation University of China,Tianjin 300300,China

    Abstract: In order to study the influence of the scarf repair on the edgewise compression performance of the foam sandwich structure with plain weave faceplate, by the use of USDFLD material subroutine, Tsai-Wu criterion was set up for material failure of the faceplates. The linear buckling mode is considered as the initial defect, and the progressive failure analysis model for the repaired foam sandwich structure with plain woven faceplates is established. The analysis results show that the validity of the model is verified by experiments. The repair quality of the sandwich structure is better when the slope ratio of scarf repair is 1:10 to 1:22, and the edgewise compression strength increases with the decrease of the slope ratio of scarf repair. When the overlap width of the additional layer is in the range of 5mm to 15mm, the repair quality is better. The edgewise compression strength decreases with the increase of the overlap width of the additional layer, but the change range is not great. The analytical methodology involving the pure buckling analysis is more accurate to predict the maximum edgewise compression load for the intact test pieces and the repaired test pieces with little damage. The failure mechanism of the sandwich structure under edgewise compression is that after the failure of the core material, the faceplate buckles locally, and then the faceplate fails when the test piece reaches its edgewise compression strength.

    Key Words: plain weave; composite; repair; foam sandwich; edgewise compression; progressive failure analysis

    猜你喜歡
    修補(bǔ)復(fù)合材料
    金屬復(fù)合材料在機(jī)械制造中的應(yīng)用研究
    纖維素基多孔相變復(fù)合材料研究
    民機(jī)復(fù)合材料的適航鑒定
    復(fù)合材料無損檢測探討
    電子測試(2017年11期)2017-12-15 08:57:13
    淺析橋梁工程混凝土裂縫控制及修補(bǔ)技術(shù)
    科技視界(2016年18期)2016-11-03 23:35:38
    淺談建筑混凝土結(jié)構(gòu)裂縫的修補(bǔ)設(shè)計與方法
    試議混凝土施工技術(shù)在道橋施工中的應(yīng)用
    科技視界(2016年16期)2016-06-29 07:46:19
    淺談瀝青路面常見病害及維護(hù)措施
    科技視界(2016年14期)2016-06-08 16:40:41
    瀝青路面坑槽成因及修補(bǔ)措施分析
    科技視界(2016年11期)2016-05-23 20:11:09
    G30線永山高速公路瀝青路面修補(bǔ)常見問題原因分析及質(zhì)量控制措施
    科技視界(2016年8期)2016-04-05 08:50:00
    成人毛片60女人毛片免费| 黄片无遮挡物在线观看| 久久人人爽人人爽人人片va| 嫩草影院入口| 热re99久久精品国产66热6| 国产亚洲精品第一综合不卡 | 久久人人爽av亚洲精品天堂| 国产精品一二三区在线看| 午夜影院在线不卡| 欧美xxⅹ黑人| 亚洲一级一片aⅴ在线观看| av国产精品久久久久影院| 中文乱码字字幕精品一区二区三区| 熟妇人妻不卡中文字幕| 国产白丝娇喘喷水9色精品| a级毛片在线看网站| 久久久久视频综合| 热99国产精品久久久久久7| 免费黄网站久久成人精品| 尾随美女入室| 18在线观看网站| 2018国产大陆天天弄谢| 性色avwww在线观看| 久久久久久人人人人人| 日韩,欧美,国产一区二区三区| 91午夜精品亚洲一区二区三区| 久久久久视频综合| 最后的刺客免费高清国语| 晚上一个人看的免费电影| 亚洲精品一二三| 国产精品久久久久久久电影| 国产欧美另类精品又又久久亚洲欧美| 亚洲成人手机| 伊人久久国产一区二区| 精品国产露脸久久av麻豆| 日韩伦理黄色片| 丰满迷人的少妇在线观看| a 毛片基地| 亚洲欧美色中文字幕在线| 国产成人精品在线电影| 女人精品久久久久毛片| 最近手机中文字幕大全| 一个人免费看片子| 大话2 男鬼变身卡| 国产男人的电影天堂91| 成年人午夜在线观看视频| av在线老鸭窝| 你懂的网址亚洲精品在线观看| 日韩制服骚丝袜av| 亚洲婷婷狠狠爱综合网| 老女人水多毛片| videosex国产| 女人精品久久久久毛片| 在线 av 中文字幕| 80岁老熟妇乱子伦牲交| 午夜福利网站1000一区二区三区| 久久久久视频综合| 久久久久精品性色| 国产成人午夜福利电影在线观看| 国产精品 国内视频| 少妇的逼好多水| 国内精品宾馆在线| 日韩精品免费视频一区二区三区 | av国产精品久久久久影院| 久久精品久久久久久久性| 97超碰精品成人国产| 一级片免费观看大全| 中文精品一卡2卡3卡4更新| 一级毛片我不卡| 久久精品国产亚洲av天美| 国产成人精品福利久久| www.色视频.com| 久久国产亚洲av麻豆专区| 亚洲精品456在线播放app| 亚洲性久久影院| 中文欧美无线码| 成人18禁高潮啪啪吃奶动态图| 最近2019中文字幕mv第一页| 久久亚洲国产成人精品v| 国产精品 国内视频| 日本猛色少妇xxxxx猛交久久| 中文乱码字字幕精品一区二区三区| 国产又爽黄色视频| 日韩一本色道免费dvd| 一边亲一边摸免费视频| a级毛色黄片| 欧美激情 高清一区二区三区| 2022亚洲国产成人精品| 成人漫画全彩无遮挡| 欧美成人午夜精品| 精品人妻熟女毛片av久久网站| 久久99蜜桃精品久久| 欧美性感艳星| 妹子高潮喷水视频| 国产高清不卡午夜福利| 高清av免费在线| 成年人午夜在线观看视频| 成人综合一区亚洲| 国产精品 国内视频| 精品一区在线观看国产| 飞空精品影院首页| 色网站视频免费| 我的女老师完整版在线观看| 水蜜桃什么品种好| 国产精品不卡视频一区二区| 国产一区二区在线观看日韩| 一区二区三区乱码不卡18| 午夜福利视频在线观看免费| 欧美精品人与动牲交sv欧美| 十八禁高潮呻吟视频| 日韩,欧美,国产一区二区三区| 99香蕉大伊视频| 国产成人a∨麻豆精品| 免费黄网站久久成人精品| 免费av不卡在线播放| 国产成人91sexporn| 性色avwww在线观看| 国产一区有黄有色的免费视频| 你懂的网址亚洲精品在线观看| 精品久久蜜臀av无| 51国产日韩欧美| 免费看不卡的av| 欧美日韩成人在线一区二区| 国产成人免费无遮挡视频| 欧美日韩亚洲高清精品| 日本-黄色视频高清免费观看| av又黄又爽大尺度在线免费看| 久久久久网色| 久久久精品免费免费高清| 亚洲国产精品一区二区三区在线| √禁漫天堂资源中文www| 亚洲精品第二区| 一区二区日韩欧美中文字幕 | 亚洲av免费高清在线观看| 日韩视频在线欧美| 国产黄频视频在线观看| 咕卡用的链子| 成人午夜精彩视频在线观看| 日韩人妻精品一区2区三区| 久久ye,这里只有精品| 一边摸一边做爽爽视频免费| 国产探花极品一区二区| 国产女主播在线喷水免费视频网站| √禁漫天堂资源中文www| 青春草国产在线视频| 国产探花极品一区二区| 日韩制服骚丝袜av| 在线 av 中文字幕| 久久久国产欧美日韩av| 亚洲国产看品久久| 深夜精品福利| 精品国产露脸久久av麻豆| 波多野结衣一区麻豆| 七月丁香在线播放| 高清视频免费观看一区二区| 在线精品无人区一区二区三| 黑丝袜美女国产一区| 巨乳人妻的诱惑在线观看| 亚洲国产日韩一区二区| 免费高清在线观看视频在线观看| 国产精品人妻久久久影院| 一级毛片 在线播放| 黄色怎么调成土黄色| 国产乱来视频区| 亚洲第一av免费看| 国产精品人妻久久久久久| 国产熟女欧美一区二区| 在线观看免费日韩欧美大片| 日韩成人av中文字幕在线观看| 久久人人爽人人片av| 在线观看免费视频网站a站| 母亲3免费完整高清在线观看 | 在线天堂最新版资源| 少妇被粗大猛烈的视频| 我的女老师完整版在线观看| 久久99热这里只频精品6学生| 我要看黄色一级片免费的| 黄色 视频免费看| 99国产精品免费福利视频| 午夜福利在线观看免费完整高清在| 制服诱惑二区| 男女国产视频网站| 国产精品一区www在线观看| 天天影视国产精品| 美女脱内裤让男人舔精品视频| 街头女战士在线观看网站| 爱豆传媒免费全集在线观看| 校园人妻丝袜中文字幕| a 毛片基地| 亚洲图色成人| 亚洲精品一二三| 制服丝袜香蕉在线| 中文字幕人妻熟女乱码| 2018国产大陆天天弄谢| 国产1区2区3区精品| 9热在线视频观看99| 亚洲国产欧美日韩在线播放| 久久久久久久亚洲中文字幕| 国产成人一区二区在线| 亚洲激情五月婷婷啪啪| av一本久久久久| 亚洲精品成人av观看孕妇| 欧美精品高潮呻吟av久久| 精品国产一区二区三区久久久樱花| 国产精品一国产av| 欧美性感艳星| 激情视频va一区二区三区| 久久精品aⅴ一区二区三区四区 | 成年av动漫网址| 国产精品不卡视频一区二区| 国产女主播在线喷水免费视频网站| 亚洲伊人久久精品综合| 桃花免费在线播放| 91精品伊人久久大香线蕉| 国产日韩欧美视频二区| 在线观看一区二区三区激情| 国产精品国产av在线观看| 亚洲一码二码三码区别大吗| 亚洲av中文av极速乱| 久久人人爽av亚洲精品天堂| 亚洲国产欧美日韩在线播放| 免费av不卡在线播放| 十分钟在线观看高清视频www| 丰满迷人的少妇在线观看| 韩国高清视频一区二区三区| 亚洲欧美成人精品一区二区| 亚洲av日韩在线播放| 精品久久蜜臀av无| 黄色毛片三级朝国网站| 在线免费观看不下载黄p国产| 亚洲第一av免费看| 国产精品国产三级国产av玫瑰| 综合色丁香网| 亚洲欧美一区二区三区国产| 免费人成在线观看视频色| 中文字幕另类日韩欧美亚洲嫩草| 免费在线观看黄色视频的| 国产日韩一区二区三区精品不卡| 欧美精品亚洲一区二区| 日本欧美视频一区| 国产熟女欧美一区二区| 国产一区亚洲一区在线观看| 亚洲国产精品成人久久小说| 成年动漫av网址| 欧美精品一区二区大全| 午夜久久久在线观看| av女优亚洲男人天堂| 欧美精品高潮呻吟av久久| 午夜影院在线不卡| 成年人午夜在线观看视频| 国产亚洲精品第一综合不卡 | 欧美日韩成人在线一区二区| 亚洲av电影在线观看一区二区三区| 日韩av免费高清视频| 欧美 亚洲 国产 日韩一| 亚洲精品第二区| 国产日韩欧美亚洲二区| 五月开心婷婷网| 精品国产乱码久久久久久小说| 校园人妻丝袜中文字幕| 丝瓜视频免费看黄片| 老司机影院毛片| 爱豆传媒免费全集在线观看| 亚洲国产色片| 中国美白少妇内射xxxbb| 久久99蜜桃精品久久| 在线观看美女被高潮喷水网站| 亚洲国产欧美日韩在线播放| 黄色一级大片看看| 制服诱惑二区| 精品少妇内射三级| 91成人精品电影| 全区人妻精品视频| 免费不卡的大黄色大毛片视频在线观看| 亚洲高清免费不卡视频| 亚洲精品一区蜜桃| 亚洲欧美精品自产自拍| 国产男女内射视频| 青春草国产在线视频| 高清视频免费观看一区二区| a级毛片在线看网站| 一级片'在线观看视频| 五月天丁香电影| 日本欧美视频一区| 亚洲精品aⅴ在线观看| 国产精品三级大全| 秋霞在线观看毛片| 赤兔流量卡办理| 精品人妻偷拍中文字幕| 国产激情久久老熟女| 亚洲av免费高清在线观看| 精品一品国产午夜福利视频| av网站免费在线观看视频| 黑人猛操日本美女一级片| 黄网站色视频无遮挡免费观看| 亚洲性久久影院| 久久精品国产自在天天线| 丝袜美足系列| 美女中出高潮动态图| 国产精品 国内视频| 搡老乐熟女国产| 在线亚洲精品国产二区图片欧美| 国产在线免费精品| 爱豆传媒免费全集在线观看| 久久久久精品性色| 中文字幕另类日韩欧美亚洲嫩草| av国产久精品久网站免费入址| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 视频区图区小说| 三级国产精品片| 边亲边吃奶的免费视频| 老司机影院成人| 国产男女超爽视频在线观看| 丝袜美足系列| 成人亚洲精品一区在线观看| 啦啦啦视频在线资源免费观看| 亚洲精品一二三| 少妇的逼好多水| 日本91视频免费播放| 一区在线观看完整版| 午夜老司机福利剧场| 久久婷婷青草| 一区在线观看完整版| 男男h啪啪无遮挡| 又黄又爽又刺激的免费视频.| 伊人亚洲综合成人网| 久久综合国产亚洲精品| 99热全是精品| 18+在线观看网站| av免费观看日本| 美国免费a级毛片| 国产午夜精品一二区理论片| 亚洲天堂av无毛| 久久 成人 亚洲| 乱码一卡2卡4卡精品| 成人毛片a级毛片在线播放| 狂野欧美激情性xxxx在线观看| 蜜桃国产av成人99| 欧美日韩一区二区视频在线观看视频在线| 一级片'在线观看视频| 久久99热6这里只有精品| 欧美人与善性xxx| 亚洲美女视频黄频| 国内精品宾馆在线| 亚洲精华国产精华液的使用体验| 国产精品成人在线| 精品少妇久久久久久888优播| 日韩制服丝袜自拍偷拍| 国产免费福利视频在线观看| 丰满乱子伦码专区| 久久精品久久久久久久性| 亚洲av电影在线进入| 久久人人爽人人爽人人片va| 建设人人有责人人尽责人人享有的| 人人妻人人爽人人添夜夜欢视频| 亚洲美女视频黄频| 国产精品99久久99久久久不卡 | 2021少妇久久久久久久久久久| 高清视频免费观看一区二区| 一本—道久久a久久精品蜜桃钙片| 麻豆精品久久久久久蜜桃| 成人亚洲欧美一区二区av| 五月伊人婷婷丁香| 又黄又粗又硬又大视频| 国语对白做爰xxxⅹ性视频网站| 久久国产精品男人的天堂亚洲 | 在线观看www视频免费| 免费大片黄手机在线观看| 中文精品一卡2卡3卡4更新| 人妻系列 视频| 免费av不卡在线播放| 精品久久国产蜜桃| 性色av一级| 久久久久久久久久人人人人人人| 在线观看三级黄色| 久久99热6这里只有精品| 成人亚洲欧美一区二区av| 国国产精品蜜臀av免费| 纯流量卡能插随身wifi吗| 亚洲第一区二区三区不卡| 免费大片18禁| 免费黄频网站在线观看国产| 丝袜人妻中文字幕| 好男人视频免费观看在线| 赤兔流量卡办理| 免费播放大片免费观看视频在线观看| 国内精品宾馆在线| 亚洲国产日韩一区二区| 亚洲欧洲日产国产| 丝瓜视频免费看黄片| 免费人妻精品一区二区三区视频| 午夜老司机福利剧场| 免费在线观看完整版高清| 最新的欧美精品一区二区| 免费黄频网站在线观看国产| 国产黄频视频在线观看| 欧美人与性动交α欧美软件 | 22中文网久久字幕| 日本av手机在线免费观看| 精品一品国产午夜福利视频| 女人精品久久久久毛片| 欧美3d第一页| 天天操日日干夜夜撸| 日韩成人av中文字幕在线观看| 精品亚洲乱码少妇综合久久| 汤姆久久久久久久影院中文字幕| 丝袜美足系列| 日韩在线高清观看一区二区三区| 99热国产这里只有精品6| 亚洲欧美中文字幕日韩二区| 在线天堂中文资源库| 少妇的逼好多水| 高清毛片免费看| 国产永久视频网站| av不卡在线播放| 亚洲精华国产精华液的使用体验| 日日撸夜夜添| 在线观看www视频免费| 肉色欧美久久久久久久蜜桃| 人人澡人人妻人| 欧美97在线视频| 又大又黄又爽视频免费| 国产精品嫩草影院av在线观看| 国产一区二区三区综合在线观看 | 在线观看一区二区三区激情| 国产 一区精品| 男人爽女人下面视频在线观看| 宅男免费午夜| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品中文字幕在线视频| 一级片免费观看大全| 久久精品国产综合久久久 | 欧美国产精品一级二级三级| 久久精品国产a三级三级三级| 只有这里有精品99| 我要看黄色一级片免费的| 国产一区二区三区综合在线观看 | 国产av一区二区精品久久| 亚洲一级一片aⅴ在线观看| 亚洲色图 男人天堂 中文字幕 | 亚洲久久久国产精品| 人人澡人人妻人| 看免费成人av毛片| 精品久久久精品久久久| 纵有疾风起免费观看全集完整版| 欧美成人午夜精品| 精品国产乱码久久久久久小说| 久久精品国产a三级三级三级| 我的女老师完整版在线观看| 在线观看免费视频网站a站| 久久ye,这里只有精品| 在线观看免费视频网站a站| av片东京热男人的天堂| 天天躁夜夜躁狠狠躁躁| 成人黄色视频免费在线看| 国产精品99久久99久久久不卡 | 青春草国产在线视频| 国产福利在线免费观看视频| 美国免费a级毛片| 亚洲中文av在线| 青春草亚洲视频在线观看| 久久狼人影院| 少妇的丰满在线观看| 丝瓜视频免费看黄片| 成人黄色视频免费在线看| 少妇人妻精品综合一区二区| 亚洲av欧美aⅴ国产| 亚洲欧美清纯卡通| 亚洲人成77777在线视频| 久久久久久久精品精品| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美xxⅹ黑人| 美国免费a级毛片| 成人免费观看视频高清| 国产av国产精品国产| 中国美白少妇内射xxxbb| 精品午夜福利在线看| 国产日韩欧美视频二区| 色94色欧美一区二区| 久久人人爽人人爽人人片va| 国产精品99久久99久久久不卡 | av网站免费在线观看视频| 制服人妻中文乱码| 中文字幕制服av| 亚洲欧美一区二区三区黑人 | 国产男人的电影天堂91| 亚洲国产精品国产精品| 久久精品久久久久久噜噜老黄| 亚洲av国产av综合av卡| 国产精品久久久久成人av| 国产高清不卡午夜福利| 国产免费视频播放在线视频| 国产在线视频一区二区| 国产综合精华液| 亚洲性久久影院| 黑人欧美特级aaaaaa片| 国产又色又爽无遮挡免| 免费久久久久久久精品成人欧美视频 | 91久久精品国产一区二区三区| 国产精品麻豆人妻色哟哟久久| 最黄视频免费看| 美女视频免费永久观看网站| 汤姆久久久久久久影院中文字幕| 下体分泌物呈黄色| 寂寞人妻少妇视频99o| 亚洲欧美清纯卡通| 大陆偷拍与自拍| 国产av国产精品国产| 久久久久久久久久成人| av一本久久久久| 亚洲人成网站在线观看播放| 欧美精品国产亚洲| 国产在线一区二区三区精| 日韩伦理黄色片| 国产精品99久久99久久久不卡 | 精品卡一卡二卡四卡免费| 国产精品一区二区在线观看99| 国产av国产精品国产| 午夜福利,免费看| 激情五月婷婷亚洲| 午夜福利乱码中文字幕| 国产成人精品婷婷| 男女无遮挡免费网站观看| 一级片'在线观看视频| 天天影视国产精品| 一本—道久久a久久精品蜜桃钙片| 涩涩av久久男人的天堂| 色婷婷av一区二区三区视频| 一级爰片在线观看| 亚洲情色 制服丝袜| 亚洲经典国产精华液单| 九色亚洲精品在线播放| 久久精品国产自在天天线| 巨乳人妻的诱惑在线观看| 久久综合国产亚洲精品| 国产有黄有色有爽视频| 黄网站色视频无遮挡免费观看| 精品久久久久久电影网| 午夜福利乱码中文字幕| 免费观看无遮挡的男女| 天天操日日干夜夜撸| 女人被躁到高潮嗷嗷叫费观| 五月玫瑰六月丁香| 亚洲成av片中文字幕在线观看 | 国产 精品1| 国产熟女午夜一区二区三区| 久久国产精品男人的天堂亚洲 | 99视频精品全部免费 在线| 日韩 亚洲 欧美在线| 最黄视频免费看| 精品人妻一区二区三区麻豆| 桃花免费在线播放| 精品国产一区二区三区四区第35| 极品人妻少妇av视频| 涩涩av久久男人的天堂| 亚洲一级一片aⅴ在线观看| 国产色爽女视频免费观看| 免费不卡的大黄色大毛片视频在线观看| 丝袜喷水一区| 少妇人妻久久综合中文| 精品人妻一区二区三区麻豆| 99久久中文字幕三级久久日本| 国产高清三级在线| 日韩中文字幕视频在线看片| 巨乳人妻的诱惑在线观看| 国产精品免费大片| 日韩精品免费视频一区二区三区 | 丝袜脚勾引网站| 久久久亚洲精品成人影院| 秋霞在线观看毛片| 91aial.com中文字幕在线观看| 香蕉精品网在线| 久久久久精品久久久久真实原创| 青春草视频在线免费观看| 日日爽夜夜爽网站| av网站免费在线观看视频| 熟妇人妻不卡中文字幕| 欧美另类一区| 午夜福利,免费看| 免费观看无遮挡的男女| 全区人妻精品视频| 下体分泌物呈黄色| 国产黄色免费在线视频| a级毛色黄片| 69精品国产乱码久久久| 丝瓜视频免费看黄片| 免费不卡的大黄色大毛片视频在线观看| 欧美97在线视频| 欧美精品av麻豆av| 免费看光身美女| 久久精品国产综合久久久 | 精品久久久精品久久久| 2022亚洲国产成人精品| 国产亚洲精品第一综合不卡 | 99国产综合亚洲精品| 妹子高潮喷水视频| 久久人人爽人人爽人人片va| 日韩熟女老妇一区二区性免费视频| 免费黄色在线免费观看| 日韩大片免费观看网站| 不卡视频在线观看欧美| 人妻少妇偷人精品九色| 亚洲国产av新网站| 人妻一区二区av| 亚洲精品日本国产第一区| 国产国语露脸激情在线看| 国产午夜精品一二区理论片| 国产在视频线精品| 男人操女人黄网站| 一区二区日韩欧美中文字幕 | 亚洲,欧美,日韩| 男人爽女人下面视频在线观看| 亚洲欧美成人综合另类久久久| 免费看不卡的av| 欧美 亚洲 国产 日韩一|