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

    固體火箭發(fā)動機粘接界面蠕變損傷研究進展

    2022-08-13 09:04:24李康佳強洪夫王哲君王學仁王稼祥
    含能材料 2022年8期
    關(guān)鍵詞:襯層內(nèi)聚力推進劑

    李康佳,強洪夫,王哲君,王學仁,王稼祥

    (火箭軍工程大學導彈工程學院,陜西 西安 710025)

    0 引言

    由于可靠性高、結(jié)構(gòu)簡單、易于長期貯存等優(yōu)點,固體火箭發(fā)動機(Solid Rocket Motor,SRM)作為導彈武器系統(tǒng)的動力裝置,得到了廣泛應用。在全壽命周期內(nèi),SRM 經(jīng)歷固化降溫、長期貯存、點火建壓以及軸向過載等序貫載荷的作用,其裝藥結(jié)構(gòu)完整性會受到影響,一定程度上會影響導彈武器的可靠性及壽命[1-2]。在SRM 服役期間,貯存占據(jù)了絕大部分時間,因此貯存期間裝藥結(jié)構(gòu)受損傷的概率大大增加,特別是長期貯存的大型SRM,其推進劑藥柱和粘接界面是典型失效部位。

    針對SRM 粘接界面的性能,國內(nèi)外研究者已經(jīng)開展了一定研究。例如,Palazotto 等[3]綜述了粘接接頭中粘接斷裂和內(nèi)聚斷裂的理論研究現(xiàn)狀,并認為SRM粘接界面斷裂問題必須考慮接頭材料的本構(gòu)模型。周清春等[4]從粘接機理、組分遷移、損傷測定和斷裂分析等方面,對SRM 粘接界面性能的研究現(xiàn)狀進行了回顧。龐愛民等[5]重點從化學組成分布特征與細觀物理結(jié)構(gòu)方面闡述了高能固體推進劑/襯層界面的粘接機理和失效機理。寗暢等[6]從性能失效、組分遷移、表征方法和在線監(jiān)測4 個方面,分析了SRM 的推進劑/襯層界面貯存性能評估的重要意義。上述研究對深入認識SRM 的粘接界面性能起到了重要作用,但由于粘接界面性能影響因素的復雜性以及失效模式的多樣性,已開展的相關(guān)研究與實際貯存條件下的SRM 結(jié)構(gòu)完整性評估還有一定差距,尤其是針對長期立式貯存條件下SRM 粘接界面性能的研究還比較匱乏。王鑫等[7]研究表明長期立式貯存下SRM 藥柱蠕變損傷較小,應重點關(guān)注粘接界面的強度、脫粘等問題。因此,本文從蠕變條件下界面損傷影響因素、界面損傷測試試驗和界面損傷數(shù)值模擬3 個方面對長期立式貯存條件下SRM 粘接界面損傷研究進展進行綜述,并就研究中存在的不足展開討論和分析。

    1 蠕變條件下粘接界面損傷的影響因素

    1.1 力學因素的影響

    1.1.1 推進劑自重因素的影響

    推進劑自重是影響界面蠕變的最大因素,它使界面長期承受拉剪應力。目前固體火箭發(fā)動機采用的推進劑多為復合固體推進劑,主要由高分子黏合劑、固化劑、氧化劑、燃料添加劑(鋁粉等金屬顆粒)和增塑劑等組成,存在由粘合劑母體構(gòu)成的連續(xù)相和由無機氧化劑和金屬添加劑等固體顆粒組成的分散相2 種相態(tài),因此,推進劑具有粘彈性和高填充固體復合材料的特性[8]。推進劑的靜態(tài)粘彈性質(zhì)之一便是蠕變,它是材料失效的重要影響因素。對于聚合物來說,當溫度超過其玻璃化溫度時,蠕變因速率變大導致材料變形增大;低應力水平和常溫下,蠕變因聚合物組成分子鏈的性質(zhì)和其流動性也會存在。為使發(fā)動機在立式貯存過程中免受環(huán)境因素的影響,工業(yè)部門通常將貯存條件設(shè)置為常溫和恒定濕度。即便如此,推進劑還是會因自重產(chǎn)生蠕變效應。藥柱連續(xù)的下沉變形傳遞至粘接界面處,使其承受載荷,長期載荷導致界面結(jié)構(gòu)也產(chǎn)生蠕變效應。王鑫等[9]設(shè)置了發(fā)動機立式貯存半年的工況條件,數(shù)值模擬得出長時間立式貯存后粘接結(jié)構(gòu)最大應力位于前人工脫粘層根部,達到0.6 MPa。由此可見,不可忽視粘接界面所受的長期應力。

    大多數(shù)情況下,可忽略固體發(fā)動機在1 g 的重力加速度載荷下產(chǎn)生的應力應變,但大直徑固體火箭發(fā)動機裝藥量大,長期經(jīng)歷1 g 的加速度載荷,裝藥粘接界面和人工脫粘根部受到很大的拉伸和剪切應力,藥柱端部界面容易發(fā)生脫粘[10]。因此,必須重視推進劑重力長時間作用導致的粘接界面結(jié)構(gòu)蠕變變形。界面蠕變改變粘合界面處的應力分布,應力隨時間重新分布會降低界面結(jié)構(gòu)的強度,并且傾向于引發(fā)過早的失效。此外,粘接界面蠕變引起斷裂能退化,這就導致界面結(jié)構(gòu)的延性變差,加速斷裂發(fā)生。綜上所述,長時間的拉剪應力作用于界面使其逐漸產(chǎn)生損傷,損傷時間累積,將使粘接界面性能退化和承載能力降低,產(chǎn)生裂紋以至于界面脫粘。

    1.1.2 固化降溫后熱力學因素的影響

    固化降溫過程中SRM 承受非均布溫度載荷,固化過程是指藥漿(以黏合劑為主體,添加固化劑、交聯(lián)劑、氧化劑和金屬粉末等組分的液態(tài)推進劑)在一定的溫度和壓力下,經(jīng)過一定時間使推進劑中的粘合劑系統(tǒng)完成化學交聯(lián),形成空間網(wǎng)狀大分子的過程[11]。推進劑藥漿固化完成后形成藥柱,藥柱成形后黏結(jié)在發(fā)動機殼體內(nèi)壁;降溫階段通常是指藥柱從58 ℃降至20 ℃的過程。溫度的變化和非均勻性會影響藥柱和粘接界面的性能,由于殼體、絕熱層、襯層和藥柱的熱膨脹系數(shù)不同,導致不同部位收縮率不同,而它們相互之間又黏結(jié)在一起,藥柱的變形受到殼體、絕熱層的約束而產(chǎn)生殘余熱應力和熱應變。由固化降溫載荷下燃燒室結(jié)構(gòu)的特點以及國內(nèi)外公開的文獻[12,15]可知,在固化降溫工況下,粘接界面會因為不均勻溫度場作用產(chǎn)生殘余應力,界面材料熱膨脹差異引起的熱應力和熱應變是后續(xù)產(chǎn)生損傷的重要原因。

    國內(nèi)外針對固體火箭發(fā)動機固化降溫過程的研究,主要是基于推進劑粘彈性本構(gòu)方程和熱傳導微分方程,建立計算模型,利用有限元分析軟件對固化降溫過程進行數(shù)值模擬,分析裝藥在固化降溫過程中溫度場和應力應變場的分布規(guī)律,最終得到裝藥經(jīng)歷固化降溫后的危險點。王錕等[13]研究發(fā)現(xiàn)固化降溫載荷后裝藥應變較大區(qū)域為人工脫粘根部和翼槽表面。徐新琦等[16]模擬了五星型藥柱固化降溫過程中的應力應變瞬時響應,發(fā)現(xiàn)在推進劑藥柱前、后部位應力應變較大,尤其是在推進劑頭部有包覆層與星角過渡處。周紅梅等[17]用ABAQUS 有限元軟件數(shù)值模擬了翼柱型HTPB 固體推進劑藥柱在固化降溫時應力應變的變化規(guī)律,發(fā)現(xiàn)前、后封頭降溫速率明顯高于其他部位的降溫速率,而前翼槽的應力明顯高于其他部位。這些研究表明,在固化降溫的過程中,粘接界面特別是人工脫粘根部受殘余應力應變較大,人工脫粘根部出現(xiàn)較大應力集中,產(chǎn)生損傷的概率很大。如果只經(jīng)歷固化降溫產(chǎn)生的載荷,界面失效的可能性較小。然而固化降溫產(chǎn)生的殘余應力貫穿發(fā)動機全壽命周期,殘余應力的存在不會使界面立即表現(xiàn)出缺陷,但界面后續(xù)還要經(jīng)受其它工作載荷,由此產(chǎn)生的疊加效應會使界面所受總應力超過強度極限,最終粘接界面附近推進劑會產(chǎn)生裂紋以致于出現(xiàn)界面脫粘。

    由此可見,發(fā)動機長期貯存時,加劇界面蠕變和失效的力學因素主要為推進劑自重導致界面產(chǎn)生的長期應力和藥柱固化降溫后的殘余熱應力。界面所受的殘余應力和長期應力疊加在一起會在很大程度上損害發(fā)動機結(jié)構(gòu)完整性。

    1.2 化學因素的影響

    粘接界面結(jié)構(gòu)的化學反應種類繁多,其中推進劑/襯層界面的反應最為復雜。固體推進劑的組成以燃燒劑和氧化劑為主,兼有改善推進劑性能的多種添加成分,如固化劑、增塑劑和化學安定劑等。貼壁澆注的固體火箭發(fā)動機推進劑必須先潤濕半固化的襯層基體,讓兩者發(fā)生足夠的化學交聯(lián)反應同時固化成型,這樣推進劑/襯層界面才能形成良好的粘合。推進劑和襯層均為高分子化合物材料,一些化學組分活性高,因此推進劑/襯層界面容易發(fā)生組分遷移和其他化學反應。

    推進劑和襯層中的化學組分都會發(fā)生遷移,大量化學成分富集在界面附近,有研究表明[18-20],組分遷移的主要成分為推進劑中的增塑劑,如鄰苯二甲酸二丁酯和硝酸酯(NG/BTTN)。增塑劑是一種化學上未結(jié)合的液體,這些低分子液體組分在長期貯存期間往往會逐漸通過推進劑-襯層界面擴散。Byrd 等[21]在對HTPB 推進劑/襯層界面性能研究中發(fā)現(xiàn),吸收增塑劑的襯層會發(fā)生膨脹,引起局部應力集中,增塑劑遷移后的推進劑硬化收縮,抵抗變形能力大幅下降,導致推進劑從襯層剝離。黃志萍等[22]研究發(fā)現(xiàn)增塑劑和穩(wěn)定劑的遷移導致NEPE 推進劑/HTPB 襯墊/EPDM 絕緣材料的界面出現(xiàn)分解。龐愛民等[5]研究NEPE 推進劑/襯層界面時,發(fā)現(xiàn)硝酸酯分解遷移導致界面區(qū)域粘合劑斷鏈降解和性能弱化。組分遷移導致界面力學性能下降,如圖1 所示,高溫環(huán)境使增塑劑更易從推進劑遷移至襯層中,兩者的粘結(jié)力下降,導致粘接界面的強度降低。NEPE 推進劑活性組分多,推進劑與襯層的粘合劑體系大多不同,界面化學物理過程復雜,裝藥界面粘接問題更加突出,因此以NEPE 推進劑為代表的高能推進劑粘接界面是今后的研究重點。綜上所述,組分遷移不僅與界面材料自身的性質(zhì)有關(guān),還受到溫度、濕度和時間等因素的顯著影響,從本質(zhì)上說,組分遷移就是界面老化的表現(xiàn)形式。組分遷移會削弱界面附近推進劑的物理性能,會使其在應力集中點開裂[24],影響界面結(jié)合強度,導致發(fā)動機的性能出現(xiàn)不利。

    固體火箭發(fā)動機在長期貯存時受到推進劑自重、固化降溫以及組分遷移和老化等復雜載荷和因素的影響(圖2),導致難以準確分析粘接界面結(jié)構(gòu)的失效原因,雖然界面最終的失效模式還是強度弱化、斷裂和脫粘。因此,必須認識到各影響因素產(chǎn)生的危害,采取有效措施改善界面性能;且最重要的是,為準確掌握固體火箭發(fā)動機粘接界面結(jié)構(gòu)在長期貯存時的性能變化,須獲得載荷大小、持荷時間、老化和殘余應力大小等因素對界面損傷和破壞的影響規(guī)律,建立粘接界面的強度失效和斷裂失效準則,得到精確的界面失效判據(jù)和評估方法模型。

    圖2 貯存期間界面性能劣化因素Fig.2 Deterioration factors of interface performance during storage

    2 粘接界面蠕變損傷的試驗研究

    與許巍等[25]的定義相似,固體火箭發(fā)動機粘接界面并不是指狹義的推進劑/襯層/絕熱層/殼體之間的界面,而是指被粘接的界面層及其所有起連接作用的界面相,故對于發(fā)動機粘接界面來說,界面為功能梯度不連續(xù)的三相材料連接體。在進行試驗分析的時候,考慮的是整個界面結(jié)構(gòu)的失效破壞,這其中尤以薄弱區(qū)域推進劑/襯層界面為主。針對固體火箭發(fā)動機裝藥粘接界面的試驗表征,國內(nèi)外研究者已經(jīng)開展了界面不同部位,不同加載條件以及不同尺度下的試驗研究。本文將從宏觀和細觀2 個方面綜述粘接界面蠕變損傷的試驗研究。

    2.1 粘接界面蠕變損傷的宏觀試驗研究

    2.1.1 恒載試驗

    從工程角度看,往往是載荷而不是應力在工程應用中保持恒定。在SRM 長期貯存的過程中,界面受的是恒載而非恒應力,因此開展恒定載荷下的粘接界面結(jié)構(gòu)件蠕變試驗比恒定應力下的更具有實際參考價值。目前,對粘接界面的試驗研究多采用恒定拉伸速率下的位移加載[26-29],雖然可以深入了解不同應變率下的界面性能,但不能有效表征長期貯存時界面性能的劣化,恒載試驗比恒定應變率下進行的試驗更能代表固體火箭發(fā)動機貯存的實際情況。Gustavson 等[30]較早進行了恒載下粘接系統(tǒng)的拉伸和剝離試驗,如圖3 所示,采用低載荷、長時間的恒載試驗裝置來測定粘接界面的強度,分析表明:①低荷載下失效時間減少的原因為表觀粘結(jié)強度的降低;②界面失效模式總是內(nèi)聚破壞的,發(fā)生在界面附近的推進劑內(nèi);③推進劑水分含量以及硬度和表觀粘結(jié)強度之間具有相關(guān)性。Bills 等[31]對推進劑-粘接系統(tǒng)雙板試件進行了恒定載荷拉伸破壞試驗,發(fā)現(xiàn)了失效時間隨平均拉伸應力的增大而減少,通過失效時間的長短來評估發(fā)動機在儲存期間的壽命。雖然以上試驗可以模擬粘接界面貯存的失效過程,但是試驗的時間短,不能預示長期載荷的作用。

    2.1.2 形變測量試驗

    評判界面性能主要還是選取強度、剛度和斷裂能等參數(shù),在給定時間和載荷水平下,界面形變也可以被認為是強度、性能的損失。在粘接界面蠕變的形變測量方法上,也有研究者提出新的方法。周獻剛等[32]采用如圖4 所示的自制小型推進劑/襯層粘接界面試件進行界面蠕變試驗,通過CCD 相機測量劃定觀測區(qū)的變形,得出劃定界面區(qū)域的蠕變-位移曲線,可以發(fā)現(xiàn)界面附近的推進劑變形占據(jù)了大部分位移。但其試驗的有效性還有待驗證,主要問題在于小型界面試件的破壞能否代表真實的發(fā)動機界面失效,以及界面位移測量區(qū)域的劃定范圍如何確定。崔輝如等[33]設(shè)計了帶有初始裂紋的DCSB 試樣,如圖5 所示,在試樣上安裝軟位移計測量蠕變位移。但其試件是選用工業(yè)粘合劑、推進劑和絕緣材料自行制成,得出的試驗結(jié)果能夠在一定程度上反應粘彈性界面的規(guī)律,但不能反應真實的發(fā)動機界面性能變化。

    圖4 自制界面蠕變變形測量區(qū)域[32]Fig.4 Self-made interface creep deformation measurement area

    圖5 DCSB 試樣蠕變試驗的位移測量[33]Fig.5 Displacement measurement of DCSB specimen in creep test

    其它粘接界面的蠕變試驗也可以給我們帶來一些思考和借鑒,膠接接頭在工業(yè)中應用廣泛,由于膠粘劑通常情況下都為高聚物材料,膠接接頭的蠕變不可忽視。除了失效時間外,還可以采用殘余斷裂能等分析界面的損傷程度。Neto 等[34]利用端缺口彎曲(ENF)試樣使其經(jīng)受不同的蠕變載荷和時間,然后進行準靜態(tài)試驗,研究了蠕變載荷和時間對粘合劑的殘余II 型斷裂能的影響。在界面蠕變位移測量方面,Han 等[35]為了研究蠕變應力對接頭老化響應的影響,將單搭接接頭約束在彈簧加載夾具中,蠕變引起的接頭延伸通過定期測量預壓縮彈簧的位移來確定。Khabazaghdam 等[36]研究了石墨烯增強粘合劑鋁-鋁粘接接頭的蠕變行為,通過使用連接至樣品的引伸計測量單搭接接頭的蠕變位移。以上試驗雖然可以為固體發(fā)動機粘接界面蠕變試驗提供幫助,但是要認識到它們之間的區(qū)別,以膠粘劑為黏合材料的界面,它的被粘物往往是金屬等剛度較大的材料,在蠕變試驗中形變總發(fā)生在膠粘劑上,便于獲得界面蠕變位移。而固體火箭發(fā)動機粘接界面的材料都為具有粘彈性以及粘超彈性的軟材料,粘接界面中不同材料的形變難以從試驗中區(qū)分,應將推進劑/襯層/絕熱層界面作為一個整體結(jié)構(gòu)進行位移測量,建立一個合適的力學模型來表征界面結(jié)構(gòu)蠕變特性。

    從以上總結(jié)來看,要真實模擬發(fā)動機粘接界面在長期貯存時的性能變化,必須設(shè)計合理的試驗和選取有效的試驗參量來表征粘接界面在長期載荷下的蠕變變形和強度損失。合理的試驗應滿足的條件:(1)試件應該和整機界面的制造工藝保持一致;(2)加載方式滿足長期貯存時的特點;(3)試件受載荷后的應力應變響應要和真實發(fā)動機界面保持一致。在試驗分析方面,界面蠕變失效時間的長短、蠕變位移的大小和殘余斷裂能的變化等都可以表征蠕變條件下界面的損傷。

    2.2 粘接界面蠕變損傷的細觀試驗研究

    隨著測試技術(shù)的發(fā)展,目前用于粘接界面結(jié)構(gòu)細觀觀測的工具和技術(shù)主要有光學顯微鏡(Optical Microscope,OM)、掃描電鏡(Scanning Electron Micoscope,SEM)、X 光 電 子 能 譜(X-Ray Photoelectron Spectrometer,XPS)、計算機斷層掃描(Computed Tomography,CT)和同步輻射光源(Synchrotron Radiation Source,SRS)等,這些高精度、高性能的實驗設(shè)備和技術(shù)從細觀尺度揭示了固體發(fā)動機粘接界面的損傷和破壞過程,為分析宏觀力學性能提供了極大的幫助。

    邱欣等[37]采用CCD 光學顯微鏡對小型粘接試件界面破壞進行了觀測,發(fā)現(xiàn)AP 顆粒與襯層基體之間出現(xiàn)的脫濕損傷,是引起界面損傷的主要原因,同時表明顆粒與襯層之間的界面能小于顆粒與基體之間的界面能。李高春等[38-39]對丁羥推進劑/襯層粘接界面試件在拉伸過程的變形破壞過程進行了觀察,并采用DICM 對拉伸過程圖像進行數(shù)字處理,獲得了粘結(jié)界面的形變場,分析了其變形演化過程和破壞機理。結(jié)果表明DICM 可以有效地分析HTPB 推進劑/襯層粘結(jié)界面的形變過程。Hojo 等[40]采用掃描電鏡觀察了在原位加載下單纖維/樹脂界面的細觀損傷過程,并分析了不同加載方式下界面斷裂韌性的變化,揭示了粘結(jié)界面在不同破壞模式下的破壞機理。國內(nèi)一些研究者[41-43]也對基于SEM 原位拉伸的推進劑/襯層粘接界面破壞過程進行了分析,發(fā)現(xiàn)隨著載荷增大,推進劑內(nèi)部顆粒與基體出現(xiàn)脫濕,微孔洞逐漸增大。隨著載荷進一步增加,推進劑內(nèi)部微裂紋出現(xiàn)匯合,導致宏觀裂紋的形成,界面細觀損傷表現(xiàn)如圖6 所示。郭翔等[42]采用微CT 技術(shù),開展了NEPE 推進劑/襯層/絕熱層界面細觀結(jié)構(gòu)研究,發(fā)現(xiàn)圖像可明顯區(qū)分界面各相以及各相的基體與填充物。綜上所述,各種觀測手段和設(shè)備都有其優(yōu)點和不足,具有明顯優(yōu)勢的是SEM 和微CT,它們的分辨率可以達到納米級,特別是其與原位加載裝置結(jié)合,可以對粘接界面受載過程中的細觀結(jié)構(gòu)變化進行實時原位觀測,直觀地看出界面結(jié)構(gòu)內(nèi)部細觀尺度的變化,對損傷進行定性定量分析。但SEM與微CT 實現(xiàn)界面原位觀測的加載方式只有準靜態(tài),針對長期定載加載方式下的原位觀測還沒有開展,長期觀測的經(jīng)濟性可能是制約其開展的重要因素。

    圖6 35%應變下粘接界面細觀損傷圖[42]Fig.6 Meso damage diagram of bonding interface specimen at 35% strain[42]

    利用各種觀測設(shè)備對粘接界面的細觀試驗分為2種,一種是單純觀測界面結(jié)構(gòu),獲取界面在損傷或斷裂后的界面信息;另一種即是設(shè)計小型試件進行原位拉伸觀測試驗,觀察細觀結(jié)構(gòu)的變化及損傷演化規(guī)律。定性描述界面損傷已經(jīng)很成熟,但定量描述界面損傷的方法還有待提高,且針對粘接界面蠕變作用下的細觀損傷觀測開展的也較少,主要難點在于:①蠕變時間長,靜態(tài)加載下的原位連續(xù)觀測不容易實現(xiàn);②找到合適的參量準確表征損傷演化比較困難。細觀試驗的作用就是從不同的尺度解釋說明界面性能變化的原因,并且與宏觀試驗的參數(shù)和結(jié)論構(gòu)效,因此通過細觀試驗準確獲得微細觀力學參數(shù)至關(guān)重要。已有研究者[44]采用微CT 獲取了推進劑在單軸拉伸不同應變下推進劑的孔隙率變化規(guī)律,以孔隙率作為損傷參量分析了推進劑細觀損傷對宏觀力學性能的影響。如圖7所示,本課題組利用微CT 開展了準靜態(tài)拉伸下的界面原位觀測試驗,也采用孔隙率表征了界面的細觀損傷。

    3 粘接界面蠕變損傷的數(shù)值模擬研究

    隨著計算機性能的提升和數(shù)值算法的進步,數(shù)值模擬方法被應用的越來越廣泛,更多表征界面力學行為的方法出現(xiàn)。傳統(tǒng)的有限元模擬不能很好地表示界面在拉伸過程中的損傷情況,基于損傷力學的界面力學模型是研究界面性能退化失效的最常用方法,這其中又以內(nèi)聚力模型(Cohesive Zone Model,CZM)為代表。CZM 能夠很好地預測整個粘接界面的損傷演化 以 及 斷 裂 失 效 過 程。Dugdale[45]和Barenblatt[46]首先提出了內(nèi)聚力概念,他們把內(nèi)聚力表示為裂紋尖端張開位移的函數(shù),避免了線彈性力學中裂紋尖端的奇異性。后來的研究者在這些概念和框架下,對內(nèi)聚力模型進行了不斷地擴充和發(fā)展。

    內(nèi)聚力模型方法認為在待擴展的界面裂尖前沿區(qū)域產(chǎn)生內(nèi)聚力損傷區(qū)域,核心是建立損傷區(qū)域界面張開位移與牽引力的關(guān)系。具體方法是將界面粘結(jié)層視為一排并列的單元,通過牽引-分離法則定義每個單元的力學響應,得到每個單元損傷起始和演化的形式,通過多個單元并列進而模擬出整個界面的損傷起始和演化過程。內(nèi)聚力模型中通常包含2 個準則:損傷起始準則和損傷演化準則[47-48]。研究粘接界面蠕變的損傷特性,損傷起始判據(jù)和損傷演化規(guī)律都是值得重視的。損傷起始準則,即界面性能開始退化時應力或應變滿足的某種條件。關(guān)于內(nèi)聚力模型,ABAQUS 提供四種損傷起始判據(jù),分別為最大應力準則、最大應變準則、二次名義應力準則和二次名義應變準則。如表1所示,已有研究者應用以下準則作為內(nèi)聚力模型中固體發(fā)動機粘接界面損傷的起始判據(jù)。

    此外還有最大應變準則和二次名義應力準則作為界面損傷起始判據(jù),表達式分別為:

    采用有限元軟件進行數(shù)值分析時,現(xiàn)有的損傷起始判據(jù)選取比較隨意,忽略了不同加載方式忽略了不同加載方式下界面損傷起始點不同的事實,損傷起始準則和損傷演化準則不獨立,在Oliveira 等[54]的研究中已經(jīng)得到了證實,因此需要針對相應試驗的加載條件選取力學性能參數(shù)。由于SRM 粘接界面組成材料具有粘彈性,其力學性能具有溫度和率相關(guān)性,在粘接界面內(nèi)聚力牽引分離法則的損傷起始判據(jù)中,對于一般的拉伸剪切載荷來說,4 個主要的參數(shù)即界面的法向最大應力和最大應變、切向最大應力和最大應變,可通過單軸或多軸試驗獲取。在界面蠕變條件下,應力水平低于瞬態(tài)拉伸損傷起始應力,但損傷起始應力隨著蠕變應變增加而非線性減小。卞躍威等[55]在Drucker-Prager屈服準則的基礎(chǔ)上考慮蠕變應變的影響,定義了巖石損傷起始準則,該準則表示隨著蠕變應變的增加損傷起始面逐漸縮小,其包含的材料參數(shù)可以通過常規(guī)三軸試驗和單軸蠕變試驗確定。以上研究為我們提供借鑒,可以建立應力和蠕變應變的界面損傷起始判據(jù)函數(shù),然后編入有限元軟件用戶定義材料子程序中,得到界面蠕變條件下合理的損傷起始準則。損傷演化準則指的是斷裂過程中材料的損傷演變規(guī)律,它由具體的界面本構(gòu)關(guān)系決定。不同形式的內(nèi)聚力模型對應著不同的牽引-分離法則,即損傷演化準則。Volokh[56]利用剝離試驗對雙線性、拋物線、正弦和指數(shù)內(nèi)聚力模型進行了定性和定量比較,發(fā)現(xiàn)不僅內(nèi)聚強度和分離功是重要的,而且不同形狀的內(nèi)聚力模型在斷裂模擬中表現(xiàn)也不盡相同。Alfano[57]也利用了如圖8 所示的不同內(nèi)聚力本構(gòu)形式對界面斷裂行為的計算準確性進行了對比。對于界面具體的應力狀態(tài)和破壞模式,還需要特定形式的內(nèi)聚力模型進行仿真。這就需要根據(jù)加載方式和材料特性選擇合適的內(nèi)聚力界面本構(gòu)模型。在發(fā)動機結(jié)構(gòu)完整性分析方面,運用內(nèi)聚力模型最早和最多的是推進劑固體顆粒/粘合劑基體界面“脫濕”的研究,采用改進的內(nèi)聚力模型研究粘接界面的工作也開展了很多。陳雄等[26]將不同拉伸速率下的界面參數(shù)擬合經(jīng)驗公式,代入雙線性內(nèi)聚力模型函數(shù),構(gòu)建了率相關(guān)的HTPB 推進劑/襯層界面Ⅱ型內(nèi)聚力模型,但其不能表現(xiàn)界面的粘彈特性。周清春等[27]建立了一種新的基于內(nèi)聚力框架的雙線性界面力學模型,模型中考慮了界面分離速率的影響。馬曉琳[28]采用考慮切向和法向斷裂能不同的PPR 內(nèi)聚力模型對粘接界面Ⅰ型和Ⅱ型脫粘進行了數(shù)值仿真,認為這種模型比雙線性內(nèi)聚力模型更加適用于混合脫粘。伍鵬等[49]采用混合模式下雙線型內(nèi)聚力模型對矩形粘接界面試件拉伸過程進行了數(shù)值模擬,獲取了推進劑/襯層/絕熱層界面的界面參數(shù)和損傷破壞模式。以上研究中,粘接界面的加載速度相較蠕變還是很大,改進的內(nèi)聚力模型并不能準確用于表征極慢加載速率下界面的損傷演化,但可以為我們提供借鑒。不同的受力狀態(tài)導致發(fā)動機粘接界面的強度變化以及損傷情況不一,針對粘接界面在長期貯存過程中產(chǎn)生的界面蠕變,采用合適的界面損傷本構(gòu)十分重要。

    圖8 幾種典型內(nèi)聚力模型的比較[58]Fig.8 Comparison of several typical cohesion models

    Time Researcher Research contents 2021 Wu Peng,et al[49]Damage initiation criterion Maximum nominal stress criterion max■■■■■σn |τs |σn max,τs max■■■■■≥1 The uniaxial tensile test of standard rectangular specimens with a tensile rate of 2 mm·min-1 was carried out 2018 2016 Kuhlmann,et al[50]Niu Ranming,et al[51]Revised maximum principal stress criterion σ1 - σth =(σ0 - σth)(tm/aTap)m Maximum nominal stress criterion|Ts ||Tn |Ts max≥1 或Tn max≥1 The tensile tests under multi angle conditions were carried out with rectangular test pieces Shear test with tensile rate of 1mm/min was carried out on HTPB single lap interface specimen Quadratic nominal strain criterion 2 2014 Jiang Aimin,et al[52]■■■■■2 εn εs εn max■■■■■+■■■■■εs max■■■■■≥1 The HTPB propellant / liner rectangular bonded specimens were tested under uniaxial tension at a rate of 5 mm /min Improved quadratic nominal strain criterion 2005 Liu Fu[53]umax 2 2 nδn( )umax sδs H(un)+( )≥1 Considering the mutual coupling of tangential and normal damage,the unit step function H is used to distinguish the different properties of the interface under normal tension and compression

    蠕變條件下發(fā)動機粘接界面內(nèi)聚力模型的構(gòu)建可以參考橡膠、膠粘劑等粘彈性材料界面模型。Zehsaz等[58]通過試驗得到膠粘劑的蠕變本構(gòu)方程,利用有限元模擬了單搭接膠接接頭的蠕變行為。這種常規(guī)的有限元方法只能得出界面處應力隨時間的變化,不能直觀體現(xiàn)界面結(jié)構(gòu)蠕變的損傷演化。Musto 和Alfano[59-61]基于熱力學框架提出了一種結(jié)合損傷和粘彈性的新的速率相關(guān)內(nèi)聚帶模型,引入了一個與能量有關(guān)的損傷變量后,CZM 能很好地模擬極限慢速率下DCB 試件的界面開裂。崔輝如等[33]通過在PPR 模型中引入臨界位移、粘結(jié)強度和初始坡度參數(shù)等,建立了與時間和溫度相關(guān)的粘彈性CZM 并且用DCSB 試件蠕變試驗驗證了模型的合理性。由于推進劑/襯層/絕熱層結(jié)構(gòu)體中材料行為的粘彈特性,界面性能也具有很強的時間依賴性,要準確模擬其蠕變損傷破壞過程,必須將速率效應考慮進分離模型中。有3 種方法可以解決此類問題,一種是將界面參數(shù)與應變率相關(guān)聯(lián),得到率相關(guān)的內(nèi)聚力本構(gòu)模型,這種模型形式簡單,但不能表現(xiàn)出界面的粘彈特性;另一種是將粘彈性本構(gòu)與率無關(guān)內(nèi)聚力模型相結(jié)合,將牽引分離定律直接納入粘彈性本構(gòu)。Xu 等[62-63]的工作可以為我們提供借鑒。第三種是建立含損傷的界面蠕變本構(gòu)方程,Truong等[64]建立了基于K-R 蠕變損傷的內(nèi)聚力本構(gòu)模型,模擬了錫-硅膜界面的蠕變裂紋擴展過程。陳興[65]基于L-M 損傷本構(gòu)提出了可用于不銹鋼焊接接頭蠕變裂紋擴展分析的內(nèi)聚力模型。以上工作是將考慮應力、時間的損傷演化函數(shù)代入Norton 蠕變本構(gòu)方程,進而表征出蠕變損傷的時間效應。無論采用哪種方法定義的內(nèi)聚力本構(gòu)模型,都要通過試驗獲取內(nèi)聚力本構(gòu)方程參數(shù),并且在有限元軟件中進行二次開發(fā),利用內(nèi)聚力模型有效準確地模擬SRM 粘接界面蠕變損傷行為的方法還需要我們更深入地探討。

    4 總結(jié)與展望

    盡管國內(nèi)外已經(jīng)對SRM 粘接界面開展了包括分析模型、試驗研究以及數(shù)值模擬等方面的工作,但對SRM 的失效預測仍然面臨挑戰(zhàn)。首先,我們對界面失效機理和過程尚不能完全清楚,缺乏有效的理論解釋界面性能變化的機理。其次,在SRM 長期貯存中,界面經(jīng)受的載荷復雜多樣,累積損傷難以準確評估。提高發(fā)動機工作的可靠性必須解決這些難點問題。本文從蠕變條件下粘接界面損傷的影響因素、粘接界面蠕變損傷的試驗研究以及粘接界面蠕變損傷的數(shù)值模擬研究3 個方面,對SRM 粘接界面蠕變損傷研究進展以及當前研究中存在的不足等進行了較為詳細的論述和分析,今后還應該投入更多的研究來全面認識和提高界面長期貯存的性能,需要重點研究的內(nèi)容可以概括為:

    (1)累積損傷研究。貯存期間各種因素對界面產(chǎn)生的危害是真實存在的,除了采取有效措施改善界面性能,還需獲得載荷水平、老化時間和殘余應力大小等因素對界面的影響規(guī)律。因此要建立貯存期間界面累積損傷評估模型,綜合固化降溫、老化和組分遷移以及推進劑長期自重等載荷歷史,建立一種累積損傷模型來表征發(fā)動機界面在使用期間所承受的應力-時間-歷史的變化,得到界面損傷程度與失效時間的關(guān)系。

    (2)試驗方法和分析研究。應設(shè)計合理的試驗有效表征真實的界面貯存應力狀態(tài),關(guān)鍵在于小型試件的受力狀態(tài)要近似準確地表征真實SRM 的應力狀態(tài),破壞機制要一致。從不同尺度試驗展開分析,都應選取合適的損傷參量表示界面結(jié)構(gòu)蠕變的損傷程度,才能將試驗得出的結(jié)論應用到解決實際工程問題中。

    (3)界面本構(gòu)模型的構(gòu)建。界面本構(gòu)模型主要考慮內(nèi)聚力模型方法,當前內(nèi)聚力模型已廣泛應用于發(fā)動機粘接界面的數(shù)值模擬,但現(xiàn)有模型表征的都是較高拉伸速率下的界面力學響應,針對界面與時間相關(guān)的損傷破壞機制,還需建立極慢速率加載以及蠕變條件下的粘彈性界面內(nèi)聚力模型,提高數(shù)值模擬的可靠性。

    固體火箭發(fā)動機的廣泛應用帶來了一系列待研究解決的問題,最主要的便是壽命評估和延壽工作。發(fā)動機長期貯存時的性能退化往往發(fā)生在裝藥結(jié)構(gòu)上,而粘接界面脫粘和失效又是一個核心問題。界面設(shè)計生產(chǎn)和應用涉及許多必須同時處理的復雜問題,將問題解耦通常是困難的,不過這些問題正通過經(jīng)驗技術(shù)得到解決。然而,這些解決方案不利于界面的可靠性評估或使用壽命估計,并且對新的發(fā)動機研制和粘接界面設(shè)計幫助有限,因此必須建立一套評估粘接界面性能的科學方法和體系,著力解決工程實際問題。隨著對固體火箭發(fā)動機結(jié)構(gòu)完整性分析的重視,以及對界面力學行為的研究深入發(fā)展,相信困擾粘接界面發(fā)展的難題都會逐步得到解決。

    猜你喜歡
    襯層內(nèi)聚力推進劑
    CRTS Ⅱ型軌道板/CA 砂漿界面內(nèi)聚力模型研究
    工程力學(2022年9期)2022-09-03 03:56:04
    基于內(nèi)聚力模型的輪盤破裂轉(zhuǎn)速預測方法研究
    超磁致伸縮傳感器背襯層參數(shù)優(yōu)化的實驗研究
    電子科技(2021年2期)2021-01-08 02:26:02
    大學英語教學中影響閱讀教學的因素淺析
    青年時代(2018年11期)2018-07-21 20:02:08
    KNSB推進劑最佳配比研究
    固體火箭發(fā)動機襯層固化狀態(tài)超聲波實時監(jiān)測方法
    兵工學報(2014年5期)2014-06-27 05:41:35
    含LLM-105無煙CMDB推進劑的燃燒性能
    火炸藥學報(2014年5期)2014-03-20 13:17:53
    無鋁低燃速NEPE推進劑的燃燒性能
    火炸藥學報(2014年5期)2014-03-20 13:17:53
    DNTF-CMDB推進劑的燃燒機理
    火炸藥學報(2014年1期)2014-03-20 13:17:27
    固體火箭發(fā)動機光固化襯層研究①
    搡老熟女国产l中国老女人| 亚洲色图av天堂| 久久久久久亚洲精品国产蜜桃av| 色婷婷久久久亚洲欧美| 亚洲成av片中文字幕在线观看| 久久久国产成人免费| 亚洲欧美激情在线| 天天添夜夜摸| 久久精品国产99精品国产亚洲性色 | 国产精品 国内视频| 国产亚洲精品久久久久5区| 丝袜人妻中文字幕| 精品人妻1区二区| 别揉我奶头~嗯~啊~动态视频| 黑丝袜美女国产一区| tocl精华| 后天国语完整版免费观看| 丁香六月天网| 久久九九热精品免费| 一个人免费看片子| 天堂8中文在线网| 亚洲九九香蕉| 大片免费播放器 马上看| 巨乳人妻的诱惑在线观看| 在线观看免费午夜福利视频| 国产淫语在线视频| 精品国内亚洲2022精品成人 | 亚洲av日韩精品久久久久久密| 黑丝袜美女国产一区| 国产亚洲精品一区二区www | 老鸭窝网址在线观看| 在线观看66精品国产| 久久午夜亚洲精品久久| 久久久久久久国产电影| 老司机午夜福利在线观看视频 | 少妇的丰满在线观看| 少妇猛男粗大的猛烈进出视频| 99riav亚洲国产免费| 18禁观看日本| 国产男女内射视频| 精品乱码久久久久久99久播| 亚洲av成人一区二区三| 黑丝袜美女国产一区| 女人爽到高潮嗷嗷叫在线视频| 少妇裸体淫交视频免费看高清 | 国产无遮挡羞羞视频在线观看| 男女床上黄色一级片免费看| 丰满人妻熟妇乱又伦精品不卡| 2018国产大陆天天弄谢| 日本av免费视频播放| 老司机福利观看| 精品国内亚洲2022精品成人 | 免费看a级黄色片| 美女视频免费永久观看网站| 精品国产一区二区久久| 国产日韩欧美在线精品| 我的亚洲天堂| 久久99一区二区三区| 国产成人欧美| 国产精品久久电影中文字幕 | 法律面前人人平等表现在哪些方面| 18禁国产床啪视频网站| 成人精品一区二区免费| 国产一区二区三区综合在线观看| 老司机在亚洲福利影院| 亚洲国产av新网站| 肉色欧美久久久久久久蜜桃| 99re在线观看精品视频| 色婷婷久久久亚洲欧美| 操出白浆在线播放| 人成视频在线观看免费观看| 亚洲精品美女久久久久99蜜臀| 男女午夜视频在线观看| 久久人妻av系列| 天天操日日干夜夜撸| 一区福利在线观看| 色播在线永久视频| 成人国语在线视频| 午夜福利在线免费观看网站| 侵犯人妻中文字幕一二三四区| 亚洲成国产人片在线观看| 亚洲成人免费电影在线观看| 免费人妻精品一区二区三区视频| 免费少妇av软件| 色播在线永久视频| 久久中文字幕人妻熟女| 成人特级黄色片久久久久久久 | 动漫黄色视频在线观看| 香蕉丝袜av| 自线自在国产av| 国产成人欧美在线观看 | 99热网站在线观看| 精品一品国产午夜福利视频| 精品卡一卡二卡四卡免费| 最近最新免费中文字幕在线| 国产一区有黄有色的免费视频| 免费不卡黄色视频| 国产有黄有色有爽视频| 十八禁高潮呻吟视频| 久久ye,这里只有精品| 99热国产这里只有精品6| 国产欧美日韩一区二区三区在线| 久久久久久久精品吃奶| 国产成人精品久久二区二区91| 亚洲精品中文字幕在线视频| 国产一区二区 视频在线| 亚洲国产欧美一区二区综合| 欧美日韩精品网址| 男人操女人黄网站| 久久亚洲精品不卡| 亚洲专区字幕在线| 自拍欧美九色日韩亚洲蝌蚪91| 黄色视频,在线免费观看| 黑人巨大精品欧美一区二区蜜桃| 丰满少妇做爰视频| 午夜老司机福利片| 成人亚洲精品一区在线观看| 妹子高潮喷水视频| av一本久久久久| 18禁观看日本| 在线观看免费视频网站a站| 亚洲av国产av综合av卡| 99riav亚洲国产免费| 天天影视国产精品| 久久国产亚洲av麻豆专区| 亚洲一码二码三码区别大吗| 91成人精品电影| 久久精品亚洲熟妇少妇任你| 国产单亲对白刺激| 人人妻人人澡人人爽人人夜夜| 一区二区三区激情视频| 国产成人av激情在线播放| 久久性视频一级片| 午夜福利免费观看在线| 大型av网站在线播放| 国产亚洲午夜精品一区二区久久| 欧美激情极品国产一区二区三区| 日韩免费av在线播放| 亚洲av片天天在线观看| 80岁老熟妇乱子伦牲交| 国产高清视频在线播放一区| 搡老乐熟女国产| 国产精品 国内视频| 欧美+亚洲+日韩+国产| 国产激情久久老熟女| 99在线人妻在线中文字幕 | 一本色道久久久久久精品综合| 欧美老熟妇乱子伦牲交| 国产人伦9x9x在线观看| 国产黄频视频在线观看| 欧美黑人精品巨大| 国产精品免费视频内射| 成人国语在线视频| av有码第一页| 91成人精品电影| 亚洲色图av天堂| 国产aⅴ精品一区二区三区波| 大香蕉久久成人网| 露出奶头的视频| 免费在线观看影片大全网站| 亚洲午夜理论影院| 老司机福利观看| 亚洲精华国产精华精| kizo精华| 黄色视频在线播放观看不卡| 男女无遮挡免费网站观看| 亚洲欧美一区二区三区久久| 国产成人一区二区三区免费视频网站| 免费黄频网站在线观看国产| 久久精品国产99精品国产亚洲性色 | avwww免费| 涩涩av久久男人的天堂| 国产深夜福利视频在线观看| 乱人伦中国视频| 国产高清视频在线播放一区| 天天躁夜夜躁狠狠躁躁| 99国产极品粉嫩在线观看| 亚洲专区中文字幕在线| 精品人妻1区二区| 丝瓜视频免费看黄片| 免费看a级黄色片| 无遮挡黄片免费观看| 亚洲第一av免费看| 天天影视国产精品| 日韩视频在线欧美| 国产精品香港三级国产av潘金莲| 视频在线观看一区二区三区| 女性被躁到高潮视频| 国产高清国产精品国产三级| 国产成人av教育| 国产麻豆69| 成人18禁在线播放| 亚洲欧美色中文字幕在线| 18禁观看日本| 亚洲少妇的诱惑av| 国产成人系列免费观看| 精品一品国产午夜福利视频| 国产不卡一卡二| 精品国产一区二区三区四区第35| h视频一区二区三区| 欧美在线黄色| 欧美黑人精品巨大| 一本一本久久a久久精品综合妖精| 宅男免费午夜| 考比视频在线观看| 考比视频在线观看| 精品乱码久久久久久99久播| 欧美变态另类bdsm刘玥| www.999成人在线观看| 亚洲久久久国产精品| 国产欧美亚洲国产| 大陆偷拍与自拍| 亚洲九九香蕉| 高清黄色对白视频在线免费看| 电影成人av| 桃红色精品国产亚洲av| 日本av免费视频播放| 国产精品美女特级片免费视频播放器 | 久久久国产精品麻豆| 看免费av毛片| 精品一区二区三区av网在线观看 | 十八禁高潮呻吟视频| 一级毛片女人18水好多| 18在线观看网站| videos熟女内射| 国产在线精品亚洲第一网站| 成人黄色视频免费在线看| 亚洲男人天堂网一区| 国产亚洲一区二区精品| 搡老乐熟女国产| 亚洲av国产av综合av卡| 久久久久精品人妻al黑| 日本五十路高清| 丰满少妇做爰视频| 精品久久蜜臀av无| 99国产精品一区二区蜜桃av | 成人三级做爰电影| 国产精品 国内视频| 老汉色av国产亚洲站长工具| 国产精品免费一区二区三区在线 | 色尼玛亚洲综合影院| 亚洲综合色网址| 久久国产精品人妻蜜桃| 丁香六月欧美| 一二三四在线观看免费中文在| 不卡一级毛片| 12—13女人毛片做爰片一| 新久久久久国产一级毛片| 巨乳人妻的诱惑在线观看| 一区二区日韩欧美中文字幕| 亚洲精品美女久久久久99蜜臀| 免费高清在线观看日韩| 大陆偷拍与自拍| 日韩欧美免费精品| 国产成人影院久久av| 国产欧美日韩一区二区精品| 亚洲一区二区三区欧美精品| 高清毛片免费观看视频网站 | 一级片'在线观看视频| 夫妻午夜视频| 久久国产精品人妻蜜桃| 欧美日韩国产mv在线观看视频| 国产精品国产高清国产av | 免费在线观看视频国产中文字幕亚洲| 亚洲av日韩精品久久久久久密| 国产在线免费精品| 欧美激情高清一区二区三区| 一级,二级,三级黄色视频| 欧美精品一区二区大全| 欧美精品av麻豆av| 老司机在亚洲福利影院| 日韩欧美三级三区| 久久久久国内视频| 亚洲av日韩精品久久久久久密| 亚洲,欧美精品.| 人人妻人人爽人人添夜夜欢视频| 国产精品一区二区在线观看99| 国产精品 国内视频| 操出白浆在线播放| 大码成人一级视频| 色视频在线一区二区三区| 中文字幕高清在线视频| 五月天丁香电影| 精品一区二区三区视频在线观看免费 | 亚洲综合色网址| 欧美大码av| 搡老熟女国产l中国老女人| 国产精品98久久久久久宅男小说| 亚洲视频免费观看视频| 9色porny在线观看| 亚洲成人手机| 女性被躁到高潮视频| 国产亚洲av高清不卡| 日韩免费av在线播放| 天天影视国产精品| videosex国产| 精品一区二区三区四区五区乱码| 日日夜夜操网爽| 在线播放国产精品三级| 老司机午夜十八禁免费视频| 男人操女人黄网站| 亚洲av欧美aⅴ国产| 亚洲性夜色夜夜综合| 99国产精品免费福利视频| 国产一区有黄有色的免费视频| 少妇粗大呻吟视频| 9热在线视频观看99| 在线观看免费高清a一片| 精品一区二区三区av网在线观看 | 午夜激情久久久久久久| 脱女人内裤的视频| 亚洲欧美日韩另类电影网站| 久久久久久人人人人人| 亚洲精品久久成人aⅴ小说| 人人妻人人澡人人爽人人夜夜| 欧美精品人与动牲交sv欧美| 日本一区二区免费在线视频| www日本在线高清视频| 又大又爽又粗| 国产欧美日韩一区二区精品| 国产精品亚洲av一区麻豆| 久久国产精品人妻蜜桃| 高清欧美精品videossex| 国产精品免费视频内射| 色在线成人网| 91精品三级在线观看| 天天操日日干夜夜撸| 国产97色在线日韩免费| 午夜福利视频精品| 99国产精品一区二区三区| av免费在线观看网站| 飞空精品影院首页| 午夜福利一区二区在线看| 欧美日韩亚洲综合一区二区三区_| 成人国语在线视频| 激情在线观看视频在线高清 | 国产日韩一区二区三区精品不卡| 日日爽夜夜爽网站| 一区二区三区激情视频| 丝袜喷水一区| 在线天堂中文资源库| 交换朋友夫妻互换小说| 国产欧美日韩综合在线一区二区| 中文字幕最新亚洲高清| 热99国产精品久久久久久7| 国产欧美日韩一区二区三区在线| 日韩欧美一区视频在线观看| 叶爱在线成人免费视频播放| 欧美成狂野欧美在线观看| 最新在线观看一区二区三区| 黄色a级毛片大全视频| 一边摸一边做爽爽视频免费| 久热这里只有精品99| 男女下面插进去视频免费观看| 久久这里只有精品19| 一区二区三区国产精品乱码| 黄网站色视频无遮挡免费观看| 女人被躁到高潮嗷嗷叫费观| 不卡一级毛片| 男人舔女人的私密视频| 精品亚洲乱码少妇综合久久| 99在线人妻在线中文字幕 | 午夜福利视频在线观看免费| 亚洲av美国av| 丝袜美腿诱惑在线| 在线观看免费午夜福利视频| 久久影院123| www日本在线高清视频| 97在线人人人人妻| 国产免费现黄频在线看| 亚洲天堂av无毛| 久久青草综合色| 国产精品98久久久久久宅男小说| 十八禁高潮呻吟视频| 色94色欧美一区二区| 精品国产乱码久久久久久小说| 怎么达到女性高潮| 999久久久精品免费观看国产| kizo精华| 国产黄色免费在线视频| 亚洲精品国产色婷婷电影| 国产精品99久久99久久久不卡| 精品视频人人做人人爽| 久久久久久久久久久久大奶| 91麻豆av在线| 欧美国产精品va在线观看不卡| 一个人免费看片子| 免费女性裸体啪啪无遮挡网站| 老司机福利观看| 国产精品久久久久久精品电影小说| 啦啦啦免费观看视频1| 狠狠狠狠99中文字幕| 高清av免费在线| 国产精品久久久久久精品电影小说| 久久精品aⅴ一区二区三区四区| 中国美女看黄片| 久久国产精品人妻蜜桃| 变态另类成人亚洲欧美熟女 | 熟女少妇亚洲综合色aaa.| 在线 av 中文字幕| 亚洲精品美女久久av网站| 国产精品美女特级片免费视频播放器 | 多毛熟女@视频| 岛国在线观看网站| 亚洲国产毛片av蜜桃av| 男女下面插进去视频免费观看| 亚洲五月婷婷丁香| 成人18禁在线播放| 成人特级黄色片久久久久久久 | 久久久久精品人妻al黑| 日本一区二区免费在线视频| 一本大道久久a久久精品| 欧美在线黄色| 日韩大码丰满熟妇| 国产99久久九九免费精品| 91字幕亚洲| 久久人妻熟女aⅴ| 久久久久精品国产欧美久久久| 757午夜福利合集在线观看| 精品人妻在线不人妻| 国产福利在线免费观看视频| 伦理电影免费视频| 国产精品久久久久久精品古装| 女性生殖器流出的白浆| 国产精品亚洲av一区麻豆| 欧美激情久久久久久爽电影 | 国产视频一区二区在线看| 欧美av亚洲av综合av国产av| av网站免费在线观看视频| 久久久精品94久久精品| 一边摸一边做爽爽视频免费| 国产欧美日韩一区二区三| 国产av又大| 女人爽到高潮嗷嗷叫在线视频| 精品一品国产午夜福利视频| 桃红色精品国产亚洲av| 窝窝影院91人妻| 在线观看一区二区三区激情| 人成视频在线观看免费观看| 亚洲情色 制服丝袜| 亚洲全国av大片| 国产精品av久久久久免费| 色综合欧美亚洲国产小说| 久久婷婷成人综合色麻豆| 亚洲一码二码三码区别大吗| 成人18禁在线播放| 黄色 视频免费看| 免费观看人在逋| 最新美女视频免费是黄的| 黄色 视频免费看| 狠狠婷婷综合久久久久久88av| 可以免费在线观看a视频的电影网站| 9191精品国产免费久久| 在线观看一区二区三区激情| a在线观看视频网站| netflix在线观看网站| 免费少妇av软件| 桃红色精品国产亚洲av| 一级毛片女人18水好多| 精品亚洲乱码少妇综合久久| 亚洲黑人精品在线| 日日爽夜夜爽网站| 亚洲一区中文字幕在线| 国产av精品麻豆| 在线观看免费高清a一片| 1024香蕉在线观看| 精品一区二区三区四区五区乱码| 成年人午夜在线观看视频| 亚洲成人国产一区在线观看| 久久这里只有精品19| 女同久久另类99精品国产91| 精品久久久久久久毛片微露脸| 国产精品免费大片| 飞空精品影院首页| 国产成人精品无人区| 怎么达到女性高潮| 免费少妇av软件| 久久久国产成人免费| 免费日韩欧美在线观看| 深夜精品福利| 久久精品国产99精品国产亚洲性色 | 精品福利永久在线观看| 91精品国产国语对白视频| 久久精品aⅴ一区二区三区四区| 国产极品粉嫩免费观看在线| 亚洲欧美日韩另类电影网站| 国产av精品麻豆| 色视频在线一区二区三区| 99re在线观看精品视频| 亚洲天堂av无毛| 久热爱精品视频在线9| 日韩视频一区二区在线观看| 最新美女视频免费是黄的| 丁香欧美五月| 欧美日韩一级在线毛片| 亚洲精品一卡2卡三卡4卡5卡| 亚洲中文av在线| 高清在线国产一区| 亚洲美女黄片视频| 成年人黄色毛片网站| 欧美日韩av久久| 99国产精品免费福利视频| 精品久久蜜臀av无| 999久久久国产精品视频| 这个男人来自地球电影免费观看| 夜夜骑夜夜射夜夜干| 午夜福利,免费看| 757午夜福利合集在线观看| 看免费av毛片| 午夜福利欧美成人| 别揉我奶头~嗯~啊~动态视频| 色婷婷久久久亚洲欧美| 欧美激情极品国产一区二区三区| 夜夜夜夜夜久久久久| 中文字幕制服av| a在线观看视频网站| 欧美日韩国产mv在线观看视频| 国产免费现黄频在线看| 一级毛片女人18水好多| 国产精品二区激情视频| 一本色道久久久久久精品综合| 麻豆成人av在线观看| 国产精品国产av在线观看| 久久精品亚洲av国产电影网| 久久久久精品国产欧美久久久| 一区二区三区激情视频| 19禁男女啪啪无遮挡网站| 精品熟女少妇八av免费久了| 午夜91福利影院| 午夜免费成人在线视频| 成年人黄色毛片网站| 亚洲午夜精品一区,二区,三区| 国产高清激情床上av| 三级毛片av免费| 国产一区二区三区视频了| 久久免费观看电影| 欧美日韩国产mv在线观看视频| 巨乳人妻的诱惑在线观看| 操出白浆在线播放| 激情视频va一区二区三区| 人成视频在线观看免费观看| 咕卡用的链子| 日韩熟女老妇一区二区性免费视频| 国产日韩一区二区三区精品不卡| 精品国产一区二区三区四区第35| 亚洲国产av新网站| 麻豆乱淫一区二区| 男女之事视频高清在线观看| 丰满迷人的少妇在线观看| 国产淫语在线视频| 老司机午夜十八禁免费视频| 中文字幕制服av| 超色免费av| 不卡av一区二区三区| 黑丝袜美女国产一区| 成人精品一区二区免费| 国产欧美日韩一区二区三| 老司机福利观看| 十八禁高潮呻吟视频| 女警被强在线播放| 999精品在线视频| 欧美成狂野欧美在线观看| 国产精品亚洲一级av第二区| 香蕉丝袜av| av电影中文网址| 天天添夜夜摸| 亚洲成人手机| 亚洲精品国产精品久久久不卡| 国产成人啪精品午夜网站| 国产精品久久久av美女十八| 久久这里只有精品19| 悠悠久久av| 最新的欧美精品一区二区| 亚洲精品一卡2卡三卡4卡5卡| 热99re8久久精品国产| 国产国语露脸激情在线看| 在线观看免费视频日本深夜| 亚洲精华国产精华精| 国产精品久久久久成人av| 国产精品自产拍在线观看55亚洲 | av网站在线播放免费| 免费少妇av软件| 自拍欧美九色日韩亚洲蝌蚪91| 国产成人系列免费观看| 一本色道久久久久久精品综合| 亚洲精品在线观看二区| 黄色片一级片一级黄色片| 亚洲专区中文字幕在线| 日本一区二区免费在线视频| 老熟妇仑乱视频hdxx| 欧美日韩福利视频一区二区| 亚洲精品中文字幕在线视频| tocl精华| 亚洲精品在线美女| 中文欧美无线码| 热re99久久精品国产66热6| 成年人黄色毛片网站| 老汉色∧v一级毛片| 国产色视频综合| 欧美亚洲日本最大视频资源| 欧美成狂野欧美在线观看| 999精品在线视频| 欧美乱妇无乱码| 黄色 视频免费看| 一进一出抽搐动态| 日韩免费av在线播放| 国产成人av教育| 女人被躁到高潮嗷嗷叫费观| 美女主播在线视频| 久久久久久久大尺度免费视频| bbb黄色大片| 免费在线观看影片大全网站| 国产精品美女特级片免费视频播放器 | 亚洲精品久久成人aⅴ小说| 国产亚洲午夜精品一区二区久久|