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

    液體火箭發(fā)動(dòng)機(jī)結(jié)構(gòu)動(dòng)力學(xué)設(shè)計(jì)關(guān)鍵技術(shù)綜述

    2023-06-27 11:35:00杜大華李斌
    航空學(xué)報(bào) 2023年10期
    關(guān)鍵詞:壽命裂紋動(dòng)力學(xué)

    杜大華,李斌

    1.液體火箭發(fā)動(dòng)機(jī)技術(shù)重點(diǎn)實(shí)驗(yàn)室,西安 710100

    2.航天推進(jìn)技術(shù)研究院,西安 710100

    大推力液體火箭發(fā)動(dòng)機(jī)作為大型火箭的主動(dòng)力,是火箭的“心臟”,其性能和可靠性直接關(guān)系到航天事業(yè)的發(fā)展[1]。液體火箭主發(fā)動(dòng)機(jī)結(jié)構(gòu)復(fù)雜、服役環(huán)境極端,極致輕量化,復(fù)雜燃燒、流動(dòng)與機(jī)械振動(dòng)等誘發(fā)的結(jié)構(gòu)動(dòng)力學(xué)問題異常突出,發(fā)動(dòng)機(jī)是導(dǎo)彈武器、航天運(yùn)載器故障率最高的子系統(tǒng),可高達(dá)51%,且以渦輪泵、推力室、管路連接故障最為常見。據(jù)NASA統(tǒng)計(jì),在飛行器所發(fā)生的重大事故中,40%均與振動(dòng)有關(guān)。除去人為、超常環(huán)境因素外,在重復(fù)載荷作用下疲勞破壞是影響結(jié)構(gòu)安全性的最主要因素[2]。

    在國外,發(fā)動(dòng)機(jī)結(jié)構(gòu)振動(dòng)問題早已引起從事航空航天等科研部門的充分重視,他們對發(fā)動(dòng)機(jī)結(jié)構(gòu)進(jìn)行了系統(tǒng)性研究。在標(biāo)準(zhǔn)規(guī)范及設(shè)計(jì)準(zhǔn)則中提出了各種設(shè)計(jì)、分析和試驗(yàn)考核要求,如美軍標(biāo)《Engine structural integrity program》[3]指出,發(fā)動(dòng)機(jī)必須滿足結(jié)構(gòu)完整性的要求。建立了《Spacecraft mechanical loads analysis handbook》[4]、《Dynamic environmental criteria》[5]、NASA/ESA模型質(zhì)量評價(jià)標(biāo)準(zhǔn)[6-7]、發(fā)動(dòng)機(jī)強(qiáng)度與壽命評定條件[8]、液體火箭發(fā)動(dòng)機(jī)評估與試驗(yàn)要求[9]等,可對發(fā)動(dòng)機(jī)結(jié)構(gòu)強(qiáng)度裕度與工作可靠性進(jìn)行系統(tǒng)性分析。

    錢令希院士指出“傳統(tǒng)設(shè)計(jì)在某種程度上可以說是一種藝術(shù),而優(yōu)化設(shè)計(jì)才是科學(xué)”。航天器設(shè)計(jì)主要走以試驗(yàn)驗(yàn)證為主的路線,結(jié)構(gòu)設(shè)計(jì)已從靜態(tài)設(shè)計(jì)轉(zhuǎn)為靜態(tài)、動(dòng)態(tài)設(shè)計(jì)[10]。傳統(tǒng)的“靜力/疲勞強(qiáng)度設(shè)計(jì)、動(dòng)強(qiáng)度校核”以及“非主承力結(jié)構(gòu)等出了動(dòng)強(qiáng)度問題再說”的設(shè)計(jì)習(xí)慣會(huì)造成新研制發(fā)動(dòng)機(jī)動(dòng)力學(xué)品質(zhì)的先天不足[11]。

    發(fā)動(dòng)機(jī)結(jié)構(gòu)動(dòng)力學(xué)設(shè)計(jì)技術(shù)的發(fā)展直接關(guān)系到航天器與運(yùn)載火箭的設(shè)計(jì)水平,甚至決定著航天任務(wù)的成?。?2]。由于中國的航天工業(yè)起步相對較晚,在發(fā)動(dòng)機(jī)結(jié)構(gòu)動(dòng)力學(xué)設(shè)計(jì)、評估及標(biāo)準(zhǔn)規(guī)范等方面,與美俄航天強(qiáng)國相比仍有一定的差距。另外,隨著航天技術(shù)的快速發(fā)展,世界各航天大國逐步開展了航天運(yùn)輸系統(tǒng)的升級換代,商業(yè)化驅(qū)動(dòng)下的低成本、重復(fù)使用和航班化航天運(yùn)輸已成為世界各國關(guān)注的焦點(diǎn)[13],開展液體火箭發(fā)動(dòng)機(jī)結(jié)構(gòu)動(dòng)力學(xué)設(shè)計(jì)是可重復(fù)使用發(fā)動(dòng)機(jī)研制必須解決的關(guān)鍵技術(shù)之一。

    因此,本文針對發(fā)動(dòng)機(jī)結(jié)構(gòu)關(guān)鍵動(dòng)力學(xué)問題與動(dòng)力學(xué)設(shè)計(jì)實(shí)際工程需要,詳細(xì)介紹發(fā)動(dòng)機(jī)載荷預(yù)計(jì)、動(dòng)力學(xué)建模及模型修正、動(dòng)強(qiáng)度評估與壽命評定、結(jié)構(gòu)動(dòng)力學(xué)優(yōu)化及抗疲勞設(shè)計(jì)等關(guān)鍵技術(shù),以期為液體火箭發(fā)動(dòng)機(jī)結(jié)構(gòu)動(dòng)力學(xué)設(shè)計(jì)技術(shù)的進(jìn)一步發(fā)展提供支撐。

    1 載荷預(yù)計(jì)

    將發(fā)動(dòng)機(jī)服役的載荷環(huán)境分為靜態(tài)與動(dòng)態(tài)2大類。靜態(tài)載荷包括安裝載荷、預(yù)冷載荷、壓力、溫度、離心力和火箭加速慣性力等。動(dòng)態(tài)載荷有:燃燒組件燃燒產(chǎn)生的隨機(jī)載荷、渦輪泵流體流動(dòng)產(chǎn)生流動(dòng)“噪聲”、通道分流諧波與轉(zhuǎn)子動(dòng)不平衡、軸承退化/摩擦產(chǎn)生的旋轉(zhuǎn)諧波載荷、管道/閥門內(nèi)流體脈動(dòng)與沖擊、發(fā)動(dòng)機(jī)及火箭飛行中的氣動(dòng)與噪聲、級間分離沖擊載荷、地面運(yùn)輸載荷,二級、上面級發(fā)動(dòng)機(jī)工作前還將經(jīng)歷一級或下面級發(fā)動(dòng)機(jī)工作時(shí)產(chǎn)生的聲振環(huán)境。

    Harry[14]認(rèn)為,航天飛機(jī)主發(fā)動(dòng)機(jī)(SSME)故障主要是由于設(shè)計(jì)不足,其根源是缺乏對動(dòng)態(tài)載荷的充分了解。MIL-STD-810強(qiáng)調(diào)了載荷預(yù)示、響應(yīng)分析和數(shù)據(jù)統(tǒng)計(jì)方法的重要性,需將強(qiáng)度與載荷環(huán)境一并考慮。目前,對靜態(tài)載荷的研究比較充分,而針對動(dòng)態(tài)載荷的產(chǎn)生機(jī)理、獲取途徑、分析方法等方面還有待進(jìn)一步研究。

    1.1 動(dòng)力學(xué)試驗(yàn)

    NASA研究人員Jack Bunting說:“不在現(xiàn)在試驗(yàn),就在將來試驗(yàn)”。發(fā)動(dòng)機(jī)典型零部組件及整機(jī)動(dòng)力學(xué)試驗(yàn)在其研制中是一項(xiàng)非常重要的工作。近年來,隨著試驗(yàn)仿真技術(shù)的日漸成熟,“計(jì)算試驗(yàn)”已成為動(dòng)力學(xué)試驗(yàn)的有力補(bǔ)充。

    長期以來,航天結(jié)構(gòu)的靜/動(dòng)強(qiáng)度主要依賴于力學(xué)環(huán)境試驗(yàn)考核[15]。然而,由于試驗(yàn)?zāi)芰Φ南拗坪驮囼?yàn)成本的約束,在地面上對預(yù)期使用環(huán)境的真實(shí)模擬有一定難度。因此,即便是對力學(xué)環(huán)境考核試驗(yàn),仍需解決以下2個(gè)問題:一是實(shí)驗(yàn)室環(huán)境試驗(yàn)條件如何代表真實(shí)的工作載荷環(huán)境;二是環(huán)境試驗(yàn)條件與真實(shí)的工作載荷環(huán)境不一致所造成的影響如何評估(如天地一致性問題)。其本質(zhì)是地面試驗(yàn)與工作環(huán)境下結(jié)構(gòu)損傷的等效問題,而問題的解決依賴于準(zhǔn)確的動(dòng)力學(xué)分析和試驗(yàn),以及動(dòng)強(qiáng)度與壽命定量評估方法。

    1.2 載荷數(shù)據(jù)測試

    目前,發(fā)動(dòng)機(jī)載荷數(shù)據(jù)的主要獲取途徑有數(shù)值計(jì)算、地面試驗(yàn)測試與飛行遙測。利用多場計(jì)算軟件及高性能仿真分析平臺,可得到渦輪泵、燃燒組件、管路及自動(dòng)器的內(nèi)部載荷。在振動(dòng)試驗(yàn)、水力試驗(yàn)及熱試車等的基礎(chǔ)上,采用先進(jìn)傳感器技術(shù)與測試手段(如光纖應(yīng)變測量、數(shù)字圖像相關(guān)(DIC)技術(shù)等)對信息智能感知以獲得環(huán)境載荷。另外,也可借助飛行遙測數(shù)據(jù),采取包絡(luò)、統(tǒng)計(jì)分析的方法得到載荷數(shù)據(jù)。

    在對測試數(shù)據(jù)的處理中,目前已制定了相應(yīng)的規(guī)范,如《液體火箭發(fā)動(dòng)機(jī)試驗(yàn)數(shù)據(jù)處理方法》(QJ1492A-2005)等。通過多源信息融合、大數(shù)據(jù)快速挖掘及數(shù)據(jù)特征準(zhǔn)確提取,可掌握發(fā)動(dòng)機(jī)結(jié)構(gòu)的動(dòng)力特征與動(dòng)態(tài)響應(yīng)規(guī)律,并可依此制定力學(xué)環(huán)境試驗(yàn)條件等。

    1.3 載荷識別

    發(fā)動(dòng)機(jī)3大振源載荷呈現(xiàn)空間分布場的形式,主要集中在中、高頻,振動(dòng)能量分布的隨機(jī)性較大。發(fā)動(dòng)機(jī)試車及飛行中,各主部件和界面處所測得的振動(dòng)加速度是結(jié)構(gòu)響應(yīng)而并非激勵(lì),將響應(yīng)數(shù)據(jù)作為振動(dòng)激勵(lì)的輸入并不合適。因此,需要將激勵(lì)源附近測得的加速度轉(zhuǎn)換成主振源相應(yīng)位置上的力,即多源復(fù)雜載荷識別與等效,重構(gòu)結(jié)構(gòu)所服役的載荷環(huán)境是進(jìn)行發(fā)動(dòng)機(jī)動(dòng)力學(xué)設(shè)計(jì)面臨的技術(shù)挑戰(zhàn)之一。

    載荷識別屬于振動(dòng)的第2類反問題[4]。數(shù)學(xué)上大部分反問題都是不適定的,即反問題至少不滿足解的存在性、唯一性和穩(wěn)定性3個(gè)條件中的1個(gè)。動(dòng)態(tài)載荷識別是利用已知系統(tǒng)的動(dòng)態(tài)特性和響應(yīng),反演作用于結(jié)構(gòu)上的外部載荷。發(fā)動(dòng)機(jī)多源載荷識別是對發(fā)動(dòng)機(jī)部件界面激勵(lì)載荷的辨識[16-17],即動(dòng)態(tài)激勵(lì)力的直接測量或間接計(jì)算。

    Christensen等[18]給出了幾種再現(xiàn)發(fā)動(dòng)機(jī)振動(dòng)環(huán)境的方法。Belelloch[19]提出了一種進(jìn)行發(fā)動(dòng)機(jī)振動(dòng)響應(yīng)匹配的方法。通過經(jīng)驗(yàn)貝葉斯(EB) 方法估計(jì)有限脈沖響應(yīng) (FIR) 函數(shù),利用FIR函數(shù) 、加速度響應(yīng)和貝葉斯公式,并運(yùn)用哈密頓蒙特卡羅 (HMC) 采樣器開展了發(fā)動(dòng)機(jī)的多源載荷識別[20]。目前,關(guān)于外部動(dòng)載荷(特別是集中載荷)識別方面的研究較多,而關(guān)于內(nèi)部動(dòng)載荷、分布隨機(jī)動(dòng)載荷辨識方面的研究則較少。另外,研究工作主要集中在針對確定性結(jié)構(gòu)動(dòng)載荷的識別方法,而對非確定性結(jié)構(gòu)動(dòng)載荷準(zhǔn)確識別卻有一定的難度[21]。

    對于液體火箭發(fā)動(dòng)機(jī),可采用基于響應(yīng)一致的力學(xué)環(huán)境復(fù)現(xiàn)方法,以進(jìn)行發(fā)動(dòng)機(jī)多源振動(dòng)載荷傳遞路徑分析(TPA)、識別、等效與重構(gòu)。路廣霖等[22]采用了一種基于Tikhonov加權(quán)正則化的TPA技術(shù),通過多源振動(dòng)載荷傳遞路徑試驗(yàn)系統(tǒng),對3大振源同時(shí)激勵(lì)以獲取激振力與結(jié)構(gòu)響應(yīng)的相關(guān)特征,并分析不同激勵(lì)源在不同頻率下對結(jié)構(gòu)響應(yīng)的貢獻(xiàn)量,如圖1所示。

    圖1 多源振動(dòng)載荷傳遞路徑試驗(yàn)Fig.1 Transfer path test of multi-source vibration loads

    此外,在多源振動(dòng)載荷等效、傳遞機(jī)理及力學(xué)環(huán)境重構(gòu)中,西安航天動(dòng)力研究所研究了發(fā)動(dòng)機(jī)多源振動(dòng)載荷時(shí)/頻域等效方法,多源多路徑振動(dòng)耦合傳遞規(guī)律及多源振動(dòng)對響應(yīng)的貢獻(xiàn)量,提出了基于逆矩陣法的界面載荷識別技術(shù)及基于響應(yīng)一致原則的力學(xué)環(huán)境復(fù)現(xiàn)方法[23],以此為基礎(chǔ)并結(jié)合試驗(yàn)構(gòu)建了液體火箭發(fā)動(dòng)機(jī)試驗(yàn)全景數(shù)據(jù)庫(圖2)。

    圖2 發(fā)動(dòng)機(jī)試驗(yàn)全景數(shù)據(jù)庫Fig.2 Test panorama database of engines

    1.4 振動(dòng)疲勞載荷譜編制

    發(fā)動(dòng)機(jī)振動(dòng)疲勞載荷譜編制是進(jìn)行結(jié)構(gòu)壽命預(yù)測的重要工作??紤]到發(fā)動(dòng)機(jī)振動(dòng)載荷的復(fù)雜性,進(jìn)行準(zhǔn)確的載荷譜編制有一定難度。因此,在考慮綜合環(huán)境效應(yīng)下,通過開展編譜方法研究,提出適用于火箭發(fā)動(dòng)機(jī)振動(dòng)疲勞載荷譜的編制技術(shù)。

    對載荷譜簡化與加重,必須保證與原始譜的一致性,遵循損傷等效、有限目標(biāo)、主要損傷部位和損傷模式一致等原則。對載荷譜加速編輯處理[24]分時(shí)域、頻域及時(shí)頻域方法。時(shí)域載荷譜只保留原始數(shù)據(jù)的幅值、均值、載荷頻次信息。頻域譜能有效反映數(shù)據(jù)的頻率及能量信息,而現(xiàn)有頻域疲勞加速方法研究主要集中在單軸方面,相關(guān)多軸同步加速的方法較少;在加速耐久性試驗(yàn)中對多軸載荷譜進(jìn)行編制,可采用小波變換法[25]。時(shí)頻耦合疲勞載荷譜編制[26]是基于雨流計(jì)數(shù)法,可獲得時(shí)序下各載荷點(diǎn)的幅值和均值,以及相應(yīng)發(fā)生的頻次;在數(shù)據(jù)處理中,保留時(shí)間軸信息,采取改進(jìn)的包含時(shí)間信息的雨流計(jì)數(shù)法(RFC),實(shí)現(xiàn)對載荷序列的載荷、頻次、頻率的耦合統(tǒng)計(jì),但此方法多針對共振模式的振動(dòng)疲勞問題。

    2 動(dòng)力學(xué)建模及模型修正

    2.1 結(jié)構(gòu)動(dòng)力學(xué)建模

    有限元法(FEM)、子結(jié)構(gòu)法已成為解決結(jié)構(gòu)動(dòng)力學(xué)問題的主要方法。整體FEM在處理大型復(fù)雜結(jié)構(gòu)動(dòng)力學(xué)問題時(shí)效率偏低,有時(shí)可能會(huì)失效;單級子結(jié)構(gòu)法在獲取各子結(jié)構(gòu)間組集關(guān)系及模態(tài)信息時(shí)有很多不足,導(dǎo)致求解復(fù)雜結(jié)構(gòu)動(dòng)力學(xué)問題的精度與效率得不到保證。然而,多重多級動(dòng)態(tài)子結(jié)構(gòu)法[27-28]除可彌補(bǔ)上述不足之外,對分析大型復(fù)雜結(jié)構(gòu)動(dòng)力學(xué)問題優(yōu)勢更加明顯。

    對于低頻動(dòng)力學(xué)問題,通常采取“積木式”由零部組件到整機(jī)的分層級建模思想[29-30]和子結(jié)構(gòu)試驗(yàn)建模綜合技術(shù)。先進(jìn)行建模技術(shù)研究,建立各部件有限元模型;其次,在保證部組件準(zhǔn)確建模的基礎(chǔ)上,組集得到整機(jī)三維動(dòng)力學(xué)模型;最后,依據(jù)試驗(yàn)數(shù)據(jù)對數(shù)學(xué)模型進(jìn)行修正和確認(rèn),最終得到準(zhǔn)確的發(fā)動(dòng)機(jī)動(dòng)力學(xué)模型。在對原始復(fù)雜的物理模型進(jìn)行合理等效時(shí),需考慮對整機(jī)動(dòng)力學(xué)影響較大的重要組件、主要連接等,而忽略對低頻特性影響較小的部件和細(xì)節(jié)等次要因素,使模型既簡單又能真實(shí)地反映結(jié)構(gòu)的低頻特性。因此,總的建模原則是準(zhǔn)確模擬發(fā)動(dòng)機(jī)的剛度和質(zhì)量分布,使建立的模型能準(zhǔn)確地反映所關(guān)心模態(tài)的變形能和動(dòng)能。建模流程如圖3所示。

    圖3 發(fā)動(dòng)機(jī)結(jié)構(gòu)低頻動(dòng)力學(xué)建模流程Fig.3 Modeling process of low frequency dynamics for engine structures

    目前,中國已建立健全了在役、預(yù)研火箭主發(fā)動(dòng)機(jī)的結(jié)構(gòu)動(dòng)力學(xué)模型(圖4)。在此基礎(chǔ)上,開展了某型120 t級液氧煤油發(fā)動(dòng)機(jī)低頻動(dòng)力學(xué)優(yōu)化設(shè)計(jì),并成功應(yīng)用于新一代三型火箭總體控制及彈簧蹺振(POGO)振動(dòng)設(shè)計(jì)[31]。

    圖4 發(fā)動(dòng)機(jī)結(jié)構(gòu)動(dòng)力學(xué)模型庫Fig.4 Dynamics model library of engine structures

    對于中頻動(dòng)力學(xué)問題,一般采用FEM與邊界元法(BEM)[32]、BEM與統(tǒng)計(jì)能量法(SEA)組合的混合元法[33],而分析高頻動(dòng)力學(xué)問題則直接采用SEA[34]。對于FEM,在低頻易于實(shí)現(xiàn)且精度較高,但隨著頻率增高,模態(tài)數(shù)急劇增多,計(jì)算量快速增大,且在處理不確定因素較多的問題時(shí)能力不足。SEA把研究對象從總體中抽取出來,忽略某些細(xì)節(jié),而主要關(guān)注結(jié)構(gòu)響應(yīng)在頻域、空間上的統(tǒng)計(jì)特征,該方法適用于分析含有高頻、密模態(tài)復(fù)雜系統(tǒng)的耦合動(dòng)力學(xué)問題。

    2.2 模型修正技術(shù)

    由于發(fā)動(dòng)機(jī)結(jié)構(gòu)復(fù)雜、連接與約束形式多樣,材料、工藝、安裝誤差及非線性,無法準(zhǔn)確定量的動(dòng)態(tài)阻尼,各種理論假設(shè)的引入,對邊界條件的近似,大量的簡化及不確定因素等,使得理論有限元模型難以真實(shí)地反映實(shí)際結(jié)構(gòu)的動(dòng)力特征與行為,需要采用模型修正技術(shù)對建立的初始有限元模型進(jìn)行修正[35]。

    通常,模型修正最終都要轉(zhuǎn)化為一個(gè)優(yōu)化問題。主要的模型修正方法有:① 基于模態(tài)參數(shù)的模型修正[36],由于模態(tài)阻尼受環(huán)境及測量條件限制的影響較大,目前基于模態(tài)參數(shù)的模型修正研究,主要集中在以模態(tài)頻率及模態(tài)振型為目標(biāo)函數(shù)的修正方法[37]。② 基于頻響函數(shù)(FRF)的模型修正[38],模態(tài)參數(shù)只表示結(jié)構(gòu)在共振頻率附近的動(dòng)力學(xué)特征,而FRF則表示結(jié)構(gòu)在一個(gè)較寬頻帶范圍內(nèi)的動(dòng)力學(xué)特性,故直接使用FRF進(jìn)行模型修正可以利用到更為豐富的信息,還避免了模態(tài)參數(shù)識別帶來的誤差影響。③ 基于時(shí)域動(dòng)力學(xué)響應(yīng)的模型修正[39],其實(shí)驗(yàn)數(shù)據(jù)獲取簡單,對于結(jié)構(gòu)非線性、運(yùn)行狀態(tài)結(jié)構(gòu)等難以利用模態(tài)參數(shù)或FRF數(shù)據(jù)進(jìn)行模型修正的情況更加適用。

    目前,NASA[6]、ESA[7]已建立較為完善的模型質(zhì)量評價(jià)標(biāo)準(zhǔn),如SMC-S-004、Chuang[40]、Link和Friswell[41]提出了評價(jià)模型修正方法及模型修正效果的準(zhǔn)則。此外,西安航天動(dòng)力研究所也提出了液體火箭發(fā)動(dòng)機(jī)結(jié)構(gòu)動(dòng)力學(xué)建模、模型修正與評估的具體要求[12]。

    3 動(dòng)強(qiáng)度評估與壽命評定技術(shù)

    動(dòng)強(qiáng)度評估、壽命評定是進(jìn)行發(fā)動(dòng)機(jī)結(jié)構(gòu)動(dòng)力學(xué)設(shè)計(jì)的重要工作[42]。發(fā)動(dòng)機(jī)結(jié)構(gòu)壽命是通過分析和試驗(yàn)確定結(jié)構(gòu)在給定使用環(huán)境載荷下的可靠性壽命。Schijve[43]從哲學(xué)的角度提出材料與結(jié)構(gòu)的疲勞應(yīng)該研究什么。疲勞問題涉及到的研究內(nèi)容廣泛豐富,是一個(gè)多學(xué)科問題,見圖5。

    圖5 結(jié)構(gòu)疲勞多學(xué)科問題Fig.5 Multidisciplinary problems of structural fatigue

    確定疲勞壽命有試驗(yàn)分析法和試驗(yàn)法。研究疲勞分析方法的目標(biāo)之一,是降低疲勞分析對于大量試驗(yàn)的依賴性,減少分析方法中的經(jīng)驗(yàn)成分,目前已形成了多種疲勞壽命分析方法。

    MIL-HDBK-1530指出,發(fā)動(dòng)機(jī)必須滿足耐久性、損傷容限設(shè)計(jì)等要求。然而,液體火箭發(fā)動(dòng)機(jī)因結(jié)構(gòu)復(fù)雜、載荷環(huán)境極端、失效模式多樣、源數(shù)據(jù)少、子樣小、地面試驗(yàn)難度大等諸多因素影響,目前雖然對疲勞現(xiàn)象已有了較好的定性了解,但是壽命定量的準(zhǔn)確度仍然有限,針對發(fā)動(dòng)機(jī)具體結(jié)構(gòu)/問題的壽命預(yù)測方法還處于發(fā)展階段,因此很好地認(rèn)識并理解所有的相關(guān)問題非常重要。

    3.1 耐久性評定方法

    耐久性[44]是發(fā)動(dòng)機(jī)在一定時(shí)間內(nèi)抵抗開裂、腐蝕、材料退化、熱降格、涂鍍層熱失配、磨損及內(nèi)/外部物體損傷的能力。根據(jù)應(yīng)力應(yīng)變響應(yīng)表征方式的不同,可分為頻域法與時(shí)域法。壽命預(yù)測流程如圖6所示。

    圖6 疲勞壽命預(yù)測流程Fig.6 Process of fatigue life prediction

    細(xì)節(jié)疲勞額定值法(DFR)[45],是在名義應(yīng)力法的基礎(chǔ)上提出的一種有效耐久性分析方法,以結(jié)構(gòu)細(xì)節(jié)疲勞額定值作為固有疲勞特性的度量,將結(jié)構(gòu)的疲勞/耐久性分析簡化為類似靜強(qiáng)度校核的形式,目前已廣泛應(yīng)用于結(jié)構(gòu)的疲勞壽命初步設(shè)計(jì)及耐久性分析。當(dāng)量原始疲勞質(zhì)量法[46-47],是通過特定載荷譜下的疲勞試驗(yàn)直接獲得結(jié)構(gòu)細(xì)節(jié)處的裂紋形成和小裂紋擴(kuò)展特性, 進(jìn)而估算細(xì)節(jié)在類似載荷譜作用下的疲勞壽命。

    對于如渦輪盤、燃燒室等典型結(jié)構(gòu),在起動(dòng)/關(guān)機(jī)瞬變、穩(wěn)定工作極端載荷作用下,將引起結(jié)構(gòu)初始損傷至疲勞、蠕變與棘輪等多種損傷交互與時(shí)序效應(yīng)[48],需研究多態(tài)交互場景下的結(jié)構(gòu)動(dòng)力行為、演變規(guī)律與失效機(jī)理,清楚壽命主控因素,并提出相應(yīng)的延壽技術(shù)。

    3.2 損傷容限評估技術(shù)

    損傷容限[49]是結(jié)構(gòu)在規(guī)定的不維修使用周期內(nèi),結(jié)構(gòu)抵抗由缺陷、裂紋或其他損傷引起破壞的能力,損傷容限評估是對含損傷結(jié)構(gòu)的裂紋擴(kuò)展壽命和剩余強(qiáng)度進(jìn)行分析。

    基于斷裂力學(xué)的裂紋擴(kuò)展分析方法是疲勞壽命預(yù)測的另一條技術(shù)途徑,通常適用于宏觀裂紋擴(kuò)展階段(圖7)。依據(jù)經(jīng)驗(yàn),液體火箭發(fā)動(dòng)機(jī)管路等結(jié)構(gòu)的疲勞總壽命受裂紋萌生階段的短裂紋主導(dǎo),短裂紋階段的擴(kuò)展行為與同種材料的長裂紋擴(kuò)展行為有明顯差異。在相同的應(yīng)力強(qiáng)度因子作用下,短裂紋的擴(kuò)展速率高于長裂紋,并且在低于長裂紋門檻值的情況下短裂紋仍能擴(kuò)展,即所謂的“短裂紋效應(yīng)”。

    圖7 裂紋擴(kuò)展過程Fig.7 Process of crack growth

    依據(jù)Paris裂紋擴(kuò)展速率定義,考慮殘余塑性變形,采用Newman塑性閉合模型,引入裂紋前緣三維約束效應(yīng)對長裂紋進(jìn)行修正。拓展應(yīng)用到自然萌生的短裂紋壽命分析,可采用NASGRO裂紋擴(kuò)展速率公式

    式中:a為裂紋長度;N為循環(huán)次數(shù);C、n、p和q為與材料有關(guān)的參數(shù);f為Newman裂紋閉合公式;R為應(yīng)力比;ΔK和Kmax分別為應(yīng)力強(qiáng)度因子范圍和最大值;ΔKth是裂紋擴(kuò)展閾值;Kc是材料斷裂韌性。

    為滿足安全壽命的要求,小于無損檢測(NDI)閾值的裂紋在結(jié)構(gòu)部件的服役壽命周期內(nèi)不得擴(kuò)展到非穩(wěn)定狀態(tài)[9]。劉士杰和梁國柱[50]指出,盡管裂紋擴(kuò)展理論研究早已開展,但目前在國內(nèi)實(shí)際工程應(yīng)用中對渦輪葉片壽命關(guān)鍵件一般不采取損傷容限設(shè)計(jì)。

    3.3 多/跨尺度疲勞分析技術(shù)

    疲勞發(fā)生在微觀尺度上,但最終結(jié)果則是宏觀尺度上的疲勞斷裂。疲勞損傷是一個(gè)從微觀到宏觀的跨尺度行為,一般都會(huì)經(jīng)歷微觀小裂紋、物理小裂紋、長裂紋等階段(見圖8),小裂紋擴(kuò)展占總壽命的70%~80%。微觀疲勞模型無法直接用于宏觀結(jié)構(gòu)的疲勞壽命預(yù)測,因此需利用多尺度技術(shù)建立微觀、宏觀模型之間的聯(lián)系,將微觀尺度模型的計(jì)算結(jié)果拓展到宏觀模型上[51-52]。多尺度疲勞分析法為裂紋形成階段與擴(kuò)展階段進(jìn)行統(tǒng)一分析奠定了基礎(chǔ)[53]。通過對跨尺度損傷過程進(jìn)行數(shù)值模擬,弄清結(jié)構(gòu)損傷的演化過程及機(jī)理,對結(jié)構(gòu)的剩余壽命預(yù)估從細(xì)觀的維度提供可行、高效的數(shù)值研究手段。Schijve[42]指出,研究工作的第一步應(yīng)該聚焦于對微觀尺度的觀察,采用金屬物理的理論來描述該現(xiàn)象。

    圖8 疲勞壽命的不同階段及影響因素Fig.8 Different stages and influencing factors of fatigue life

    液體火箭發(fā)動(dòng)機(jī)是由材料-結(jié)構(gòu)-發(fā)動(dòng)機(jī)組成的復(fù)雜系統(tǒng),極端環(huán)境下結(jié)構(gòu)失效是由不同尺度和層級下相關(guān)聯(lián)的事件所導(dǎo)致。因此,需研究極端環(huán)境下材料-結(jié)構(gòu)-發(fā)動(dòng)機(jī)一體化的跨尺度評估體系,進(jìn)而合理評價(jià)發(fā)動(dòng)機(jī)系統(tǒng)的概率風(fēng)險(xiǎn)。

    目前,多尺度壽命評定方法正在迅速發(fā)展。然而,至今所開展的研究多是在材料層面探討疲勞損傷如何在微/細(xì)觀尺度上產(chǎn)生、發(fā)展,并進(jìn)而引起宏觀尺度上的破壞;研究對象主要集中在實(shí)驗(yàn)室試驗(yàn)件,缺乏工程應(yīng)用實(shí)踐。其次,為深入研究損傷機(jī)制,需在模型中植入較多微/細(xì)觀細(xì)節(jié),這會(huì)導(dǎo)致模型規(guī)模龐大,從而造成高計(jì)算成本甚至無法求解,且現(xiàn)有多尺度疲勞分析方法中存在計(jì)算精度與成本之間的矛盾,以現(xiàn)有技術(shù)條件這種矛盾很難在短期內(nèi)得到解決。因此,如何建立多尺度模型以最低的計(jì)算代價(jià)從微、細(xì)觀維度研究結(jié)構(gòu)的跨尺度損傷機(jī)理及演變過程,以進(jìn)行宏觀層面上剩余壽命預(yù)測是一項(xiàng)極具挑戰(zhàn)性的工作[2]。

    3.4 疲勞壽命可靠性

    分散系數(shù)是結(jié)構(gòu)壽命可靠性的重要指標(biāo),為中值壽命N50與指定可靠度指標(biāo)下使用壽命Np的比值。因加工與使用中材料、尺寸、環(huán)境載荷、初始/邊界條件等的隨機(jī)性,疲勞壽命通常存在較大的分散性。Schijve[42]認(rèn)為,影響因素多導(dǎo)致壽命分散性大,對疲勞壽命的準(zhǔn)確定量是困難的,但在數(shù)量級上具有指導(dǎo)意義,故壽命準(zhǔn)確定量是不科學(xué)的,只能是可靠估計(jì)。確定性壽命預(yù)估很難保證壽命的可靠度,傳統(tǒng)的安全系數(shù)方法,為了保守設(shè)計(jì)通常取較大的安全系數(shù),這容易造成使用壽命的浪費(fèi),且用概率壽命進(jìn)行評估將更加科學(xué)。由于發(fā)動(dòng)機(jī)系統(tǒng)的高功效及結(jié)構(gòu)的高效率,系統(tǒng)對于設(shè)計(jì)和性能參數(shù)不確定性的敏感度就越高[54]。用概率設(shè)計(jì)裕度代替?zhèn)鹘y(tǒng)的安全系數(shù)或確定性裕度更能準(zhǔn)確地表征結(jié)構(gòu)響應(yīng)的變化。

    文獻(xiàn)[3]提到采用概率設(shè)計(jì)的思想,要求疲勞設(shè)計(jì)應(yīng)當(dāng)包含“頻率概率設(shè)計(jì)裕度”與“概率振動(dòng)應(yīng)力裕度”2個(gè)層次。文獻(xiàn)[9]對使用壽命系數(shù)(SLF)提出要求,在考慮疲勞分析系數(shù)(FAF)下,低周疲勞(LCF)取4.0,高周疲勞(HCF)為10.0。姜金朋等[55]在考慮材料屬性、幾何參數(shù)、工作載荷等變量的隨機(jī)性下,通過建立可靠性分析模型來預(yù)測液體火箭發(fā)動(dòng)機(jī)渦輪葉片概率疲勞壽命,見圖9[55]。近年來,基于數(shù)字孿生[56]、代理模型[57]等技術(shù)研究概率疲勞壽命的分析方法也得到了發(fā)展。

    圖9 概率壽命分析流程[55]Fig.9 Analysis process of probabilistic life[55]

    4 結(jié)構(gòu)動(dòng)力學(xué)優(yōu)化及抗疲勞設(shè)計(jì)

    對發(fā)動(dòng)機(jī)結(jié)構(gòu)按照動(dòng)力學(xué)指標(biāo)要求進(jìn)行設(shè)計(jì),以滿足對振動(dòng)特性、響應(yīng)及穩(wěn)定性的要求。通常情況,結(jié)構(gòu)動(dòng)力學(xué)設(shè)計(jì)實(shí)際上是動(dòng)力學(xué)優(yōu)化設(shè)計(jì)。結(jié)構(gòu)優(yōu)化[58]涉及研制全過程,需考慮發(fā)動(dòng)機(jī)結(jié)構(gòu)各設(shè)計(jì)準(zhǔn)則要求,采用優(yōu)化設(shè)計(jì)技術(shù)以實(shí)現(xiàn)結(jié)構(gòu)的預(yù)計(jì)功能。近年來,基于可靠性[59]、魯棒性[60]等優(yōu)化設(shè)計(jì)方法已成為研究的重點(diǎn)。Schijve[42]建議,承受動(dòng)態(tài)載荷的結(jié)構(gòu)應(yīng)進(jìn)行抗疲勞設(shè)計(jì),如何進(jìn)行抗疲勞設(shè)計(jì)顯然與所要實(shí)現(xiàn)的目標(biāo)相關(guān)(圖10)。此外,抗疲勞設(shè)計(jì)需要理解、想象和經(jīng)驗(yàn),是一項(xiàng)極具挑戰(zhàn)性的工作。

    圖10 與抗疲勞設(shè)計(jì)有關(guān)的課題Fig.10 Topics related to anti-fatigue design

    4.1 渦輪泵系統(tǒng)抗復(fù)合疲勞設(shè)計(jì)

    渦輪泵是泵壓式液體推進(jìn)劑供應(yīng)系統(tǒng)與控制系統(tǒng)的重要部件,是液體火箭發(fā)動(dòng)機(jī)故障率最高、危害性最大的部件,也是可重復(fù)使用發(fā)動(dòng)機(jī)中使用壽命最短的主部件[61]。幾乎所有大型液體火箭發(fā)動(dòng)機(jī)渦輪泵在研制或服役過程中均發(fā)生過疲勞導(dǎo)致的故障,涉及到渦輪泵所有零部組件,如掉葉、斷軸、軸承故障、密封失效及靜子葉柵裂紋等。渦輪泵疲勞是可重復(fù)使用發(fā)動(dòng)機(jī)必須重點(diǎn)解決的問題之一。因此,研究渦輪泵復(fù)雜載荷預(yù)示與耦合作用機(jī)理、振動(dòng)抑制及抗疲勞設(shè)計(jì)方法對于大推力發(fā)動(dòng)機(jī)研制有著重要的意義,如圖11所示。

    圖11 渦輪泵流體激振模式Fig.11 Modes of fluid excitation in turbine pumps

    中國為解決某型高壓補(bǔ)燃液氧煤油火箭發(fā)動(dòng)機(jī)渦輪泵振動(dòng)大的問題,開展了泵流體激振機(jī)理及振動(dòng)控制技術(shù)研究。通過建立渦輪泵流體激振分析模型、流固耦合模型、阻尼密封和轉(zhuǎn)子動(dòng)力學(xué)模型,揭示泵流體激振(空化[62]、動(dòng)靜干涉[63]、間隙密封流體激振等)的特征、機(jī)理及影響規(guī)律,形成系統(tǒng)有效的渦輪泵流體激振分析、評估技術(shù),突破了補(bǔ)燃發(fā)動(dòng)機(jī)渦輪泵流體激振控制技術(shù)(額定工況下振動(dòng)量級降低約62%)。

    渦輪泵動(dòng)力學(xué)設(shè)計(jì)主要面臨2大類問題:一是渦輪泵軸系及其支撐結(jié)構(gòu)的轉(zhuǎn)子動(dòng)力學(xué)問題,涉及轉(zhuǎn)子臨界轉(zhuǎn)速、動(dòng)平衡、穩(wěn)定性、次/超同步振動(dòng)、動(dòng)力學(xué)優(yōu)化、關(guān)鍵裝配參數(shù)影響、相關(guān)的轉(zhuǎn)子/密封動(dòng)力特性試驗(yàn)等;二是流致振動(dòng)問題,主要包括誘導(dǎo)輪空化激振[64-65]、動(dòng)靜干涉流體激振[63]、間隙密封流動(dòng)激振[66]及渦輪流體誘發(fā)振動(dòng)[67]等。對于一次性使用發(fā)動(dòng)機(jī),應(yīng)重點(diǎn)關(guān)注誘導(dǎo)輪空化激振、轉(zhuǎn)子系統(tǒng)振動(dòng)水平過高、流固耦合共振、輪盤行波振動(dòng)、應(yīng)力集中與加工缺陷等。對于可重復(fù)使用發(fā)動(dòng)機(jī),應(yīng)進(jìn)行軸低/高周振動(dòng)疲勞、泵高低周復(fù)合疲勞/超高周疲勞與渦輪葉盤低周熱機(jī)疲勞、蠕變加速損傷下高周振動(dòng)疲勞等研究。渦輪泵抗疲勞設(shè)計(jì)具體內(nèi)容包括:

    1) 渦輪泵流體激振控制的首要任務(wù)是控制誘發(fā)渦輪泵振動(dòng)的流體激振力,即進(jìn)行流場優(yōu)化,控制空化、流固耦合以及湍流等各種誘發(fā)振動(dòng)的因素。如對渦輪靜子葉片數(shù)及其不對稱非諧設(shè)計(jì),可降低渦輪腔內(nèi)壓力脈動(dòng);改變動(dòng)、靜組件葉片數(shù)以提高葉片的耦合階次,增加轉(zhuǎn)動(dòng)件、支撐結(jié)構(gòu)的剛度與緊固件的擰緊力矩。

    2) 動(dòng)特性分析與設(shè)計(jì),利用Campbell圖、SAFE公式、“三重點(diǎn)”法進(jìn)行振動(dòng)安全評估與抗共振設(shè)計(jì)。

    3) 動(dòng)響應(yīng)分析,考慮流體黏性阻尼與結(jié)構(gòu)阻尼,開展基于動(dòng)響應(yīng)測試反演系統(tǒng)阻尼參數(shù)的識別技術(shù)研究,如“whirligig”試驗(yàn)。

    4) 開展振動(dòng)傳遞路徑優(yōu)化與振動(dòng)響應(yīng)控制,如阻尼密封(圖12)可降低振動(dòng)傳遞率及提高轉(zhuǎn)子穩(wěn)定性,將渦輪盤-軸連接方式由整體式改為分體式(圖13)、盤上采用隔熱阻尼層可增加結(jié)構(gòu)阻尼,從而減小振動(dòng)響應(yīng)。

    圖12 阻尼密封Fig.12 Damper seals

    圖13 渦輪盤-軸連接方式Fig.13 Connection methods between turbine disk and shaf

    5) 對于渦輪內(nèi)聲共振問題,聲共振的能量源自分離流引起的脫落渦,目前控制聲共振常見的方法主要是對尾跡進(jìn)行干預(yù)[68]。

    6) 對于渦輪葉盤,需考慮起動(dòng)/關(guān)機(jī)、變工況過渡過程力熱沖擊效應(yīng)、時(shí)序效應(yīng)與交互耦合損傷的影響,研究復(fù)合疲勞分析方法。

    7) 疲勞小裂紋擴(kuò)展技術(shù),研究燃料泵帶小裂紋結(jié)構(gòu)在振動(dòng)環(huán)境下的壽命評估技術(shù),給出裂紋量化檢測標(biāo)準(zhǔn)。

    8) 對于超高周疲勞問題,需研究材料的超高周疲勞性能及壽命評估、試驗(yàn)方法。

    9) 其他還需解決寬范圍軸向力平衡、高速重載軸承、高速密封、轉(zhuǎn)子穩(wěn)定性以及增材制造結(jié)構(gòu)強(qiáng)度評估等一系列關(guān)鍵技術(shù)。

    4.2 推力室結(jié)構(gòu)熱機(jī)疲勞控制

    燃燒部件包括推力室與燃?xì)獍l(fā)生器,發(fā)生器的主要失效模式及研究方法與推力室類似,在此不做專門介紹。推力室是發(fā)動(dòng)機(jī)的主體,其故障率位于整個(gè)發(fā)動(dòng)機(jī)故障的第2位,故障主要涉及流熱聲振耦合[69]、熱機(jī)疲勞(TMF)[70]、噴注耦合振動(dòng)[71]和噴管側(cè)向力問題[72]等。研究表明,推力室的失效模式受多種因素影響,主要是塑性應(yīng)變及低周疲勞、高溫蠕變、棘輪應(yīng)變效應(yīng)等;同時(shí)這些因素以相互耦合的形式作用到結(jié)構(gòu)上,使對結(jié)構(gòu)失效機(jī)理的分析變得更加復(fù)雜,熱機(jī)疲勞是燃燒部件熱結(jié)構(gòu)最主要的失效模式之一。

    SSME、Vuclain和RS-68等高壓大熱流密度火箭發(fā)動(dòng)機(jī)在多次熱試車后,在推力室喉部上游收斂段均出現(xiàn)不同程度的裂紋,并產(chǎn)生“Dog house”現(xiàn)象[73]。RD-170、RD-180推力室在多次熱試車后,在推力室喉部附近的涂鍍層上產(chǎn)生了明顯的裂紋。中國某型氫氧火箭發(fā)動(dòng)機(jī)推力室內(nèi)壁在經(jīng)過4次試車后也出現(xiàn)梳狀裂紋,內(nèi)壁向燃燒室內(nèi)鼓起。

    目前,中國已進(jìn)行了再生冷卻推力室[70]及渦輪盤[74]等典型熱端結(jié)構(gòu)低周疲勞壽命預(yù)測及延壽技術(shù)研究。通過分析起動(dòng)/關(guān)機(jī)瞬變、穩(wěn)定工作過程復(fù)雜多場載荷耦合作用機(jī)制,研究了主要失效模式、交互損傷機(jī)理及壽命預(yù)測方法,并進(jìn)行了抗疲勞設(shè)計(jì),研究成果已為發(fā)動(dòng)機(jī)典型熱應(yīng)力結(jié)構(gòu)可靠性提升、可重復(fù)使用發(fā)動(dòng)機(jī)壽命型可靠性設(shè)計(jì)等提供關(guān)鍵技術(shù)支撐。

    對于再生冷卻推力室熱機(jī)疲勞控制與重復(fù)使用技術(shù),主要從以下幾方面考慮:

    1) 開展發(fā)動(dòng)機(jī)熱力組件高效穩(wěn)定燃燒及裕度評估技術(shù)研究,提高穩(wěn)定燃燒質(zhì)量,降低振(聲)源強(qiáng)度。

    2) 采用長壽命的內(nèi)壁材料(高激活能、高強(qiáng)度與蠕變強(qiáng)度,合理熱導(dǎo)率與延展性)。

    3) 采用可靠冷卻技術(shù),如內(nèi)壁采用多條液氣膜內(nèi)冷卻環(huán)帶[75]、人為粗糙度強(qiáng)化換熱[76]、內(nèi)壁銑槽結(jié)構(gòu)和熱障涂層(TBC)[77]等。

    4) 在保證極限承載能力的前提下,采用允許內(nèi)壁變形的低剛度外壁,可降低推力室內(nèi)壁應(yīng)力比R,但這與外壁具備良好的承載能力相矛盾。

    5) 優(yōu)化技術(shù),如新工藝、冷卻方式、冷卻通道結(jié)構(gòu)形式與尺寸參數(shù)[78]、工作時(shí)序(特別針對力熱沖擊問題)、壓力/溫度參數(shù)等。

    4.3 管路與連接結(jié)構(gòu)抗疲勞設(shè)計(jì)

    管路系統(tǒng)是發(fā)動(dòng)機(jī)介質(zhì)傳輸與能量傳遞的重要通道,是發(fā)動(dòng)機(jī)故障率較高的組件之一。故障模式多為低周大應(yīng)力疲勞與密封失效,主要是由強(qiáng)振動(dòng)、異源/多源機(jī)械載荷激勵(lì)與非協(xié)調(diào)振動(dòng)所引起,且多發(fā)生在焊接結(jié)構(gòu)、尺寸/缺口效應(yīng)顯著等位置。

    管路與連接的故障模式多樣,隨機(jī)影響因素眾多,多源/異源/多場耦合作用機(jī)理復(fù)雜,目前尚缺乏有效的評估方法,試驗(yàn)驗(yàn)證難度大,且對管路結(jié)構(gòu)進(jìn)行抗疲勞設(shè)計(jì)存在一定的困難。疲勞是連接接頭的藝術(shù)[42],采用完全理性的方法來定量考慮對疲勞的影響因素是困難的。

    至今,中國已開展了大量的發(fā)動(dòng)機(jī)管路連接結(jié)構(gòu)振動(dòng)疲勞壽命預(yù)測與驗(yàn)證技術(shù)研究工作,見圖14。建立了液體火箭發(fā)動(dòng)機(jī)力學(xué)環(huán)境載荷分析、疲勞載荷譜編制、基于損傷累積的管路連接結(jié)構(gòu)疲勞壽命預(yù)估、基于損傷等效與壽命一致性原則的壽命驗(yàn)證試驗(yàn)的統(tǒng)一分析方法[79],形成了發(fā)動(dòng)機(jī)管路結(jié)構(gòu)強(qiáng)振動(dòng)下壽命定量分析及判斷準(zhǔn)則。另外,進(jìn)行了泵后擺發(fā)動(dòng)機(jī)高溫高壓燃?xì)鈸u擺軟管關(guān)鍵結(jié)構(gòu)設(shè)計(jì)技術(shù)研究;解決了搖擺軟管構(gòu)型及其數(shù)學(xué)描述問題,完成了結(jié)構(gòu)安全裕度評價(jià)方法和設(shè)計(jì)準(zhǔn)則研究,開展了燃?xì)夂腿剂蠐u擺軟管的結(jié)構(gòu)優(yōu)化設(shè)計(jì),并在高壓搖擺特性、壽命、穩(wěn)定性及熱態(tài)搖擺等試驗(yàn)技術(shù)方面取得了重大的技術(shù)突破。

    圖14 管路及連接結(jié)構(gòu)壽命預(yù)測流程Fig.14 Life prediction process of pipeline and connection structures

    對液體火箭發(fā)動(dòng)機(jī)管路系統(tǒng)振動(dòng)進(jìn)行主被動(dòng)控制,開展管路與連接結(jié)構(gòu)的抗疲勞設(shè)計(jì),主要從以下幾方面考慮:

    1) 合理選取結(jié)構(gòu)材料。

    2) 一體化設(shè)計(jì),降低缺口效應(yīng)、尺寸效應(yīng)、剛度不匹配的影響。

    3) 減小加工、焊接、裝配等對疲勞性能的影響。

    4) 進(jìn)行管路頻率管理設(shè)計(jì)與動(dòng)增益控制設(shè)計(jì),如俄羅斯中央機(jī)械研究院要求“所有管路振動(dòng)放大系數(shù)不能超過7”。

    5) 降低結(jié)構(gòu)形式、敷設(shè)、約束對動(dòng)強(qiáng)度的影響,提高接頭疲勞強(qiáng)度;優(yōu)化管線空間走向,利用約束阻尼層;優(yōu)化卡箍結(jié)構(gòu)、數(shù)量及約束位置,采用高阻尼支撐卡箍,并提高固定基礎(chǔ)剛度[80]。

    6) 避免異源激勵(lì)、基礎(chǔ)運(yùn)動(dòng)等非協(xié)調(diào)振動(dòng)問題。

    7) 對于柔性管路流體引發(fā)的振動(dòng)問題,可按文獻(xiàn)[81]進(jìn)行分析評估。

    8) 模擬真實(shí)工作力熱環(huán)境下,開展管路振動(dòng)疲勞可靠性等效試驗(yàn)驗(yàn)證。

    焊接結(jié)構(gòu)的疲勞性能受多種因素影響,疲勞開裂中裂紋萌生對壽命的貢獻(xiàn)是可以忽略的,不能將焊接結(jié)構(gòu)的疲勞失效問題僅歸結(jié)為焊接質(zhì)量問題,需要考慮殘余應(yīng)力的影響;另外既要重視設(shè)計(jì),也要重視驗(yàn)證[82]。因此,需要加強(qiáng)焊接結(jié)構(gòu)的抗疲勞設(shè)計(jì)理論與方法研究。

    4.4 自動(dòng)器結(jié)構(gòu)動(dòng)力學(xué)設(shè)計(jì)

    液體火箭發(fā)動(dòng)機(jī)自動(dòng)器包括閥門、流量調(diào)節(jié)器等,其控制、調(diào)節(jié)過程涉及機(jī)械運(yùn)動(dòng)、瞬變流、力熱環(huán)境等多場耦合,工作條件復(fù)雜。在自動(dòng)器液流試驗(yàn)及發(fā)動(dòng)機(jī)熱試車中,曾多次發(fā)生閥芯-流體自激振蕩[83]、開/關(guān)瞬態(tài)過程水擊效應(yīng)[84]、結(jié)構(gòu)力熱環(huán)境不適應(yīng)等,導(dǎo)致閥芯磨損卡滯、密封失效、與自動(dòng)器連接結(jié)構(gòu)疲勞破壞等故障,控制系統(tǒng)組件的可靠性是決定發(fā)動(dòng)機(jī)安全工作的重要因素之一。因此,自動(dòng)器除具有優(yōu)良的靜動(dòng)態(tài)特性外,還應(yīng)具備較高的工作可靠性。在流體高壓、高速流動(dòng)、高低溫、振動(dòng)、力熱沖擊、多場耦合等復(fù)雜環(huán)境下,對自動(dòng)器進(jìn)行動(dòng)力學(xué)設(shè)計(jì)主要考慮有:

    1) 應(yīng)具備高精度調(diào)節(jié)和控制性能。

    2) 結(jié)構(gòu)應(yīng)能可靠適應(yīng)各類復(fù)雜的載荷環(huán)境。

    3) 自動(dòng)器的動(dòng)力學(xué)問題多與其接入的系統(tǒng)、工作狀態(tài)參數(shù)等相關(guān),需要研究自動(dòng)器與接入系統(tǒng)的匹配性,使所組成新系統(tǒng)的性能、穩(wěn)定性等滿足要求。

    4) 采用先進(jìn)材料、3D打印與智能優(yōu)化等新技術(shù),進(jìn)行結(jié)構(gòu)優(yōu)化設(shè)計(jì)并提高工作安全性。

    5)研究嚴(yán)苛環(huán)境下密封耐磨性、抗疲勞及長壽命等高可靠性設(shè)計(jì)技術(shù)。

    4.5 發(fā)動(dòng)機(jī)整機(jī)振動(dòng)控制

    考慮發(fā)動(dòng)機(jī)結(jié)構(gòu)的復(fù)雜性和工作環(huán)境的極端性及某些特殊要求,開展振動(dòng)精細(xì)化控制與抗強(qiáng)振設(shè)計(jì)技術(shù)研究,以增強(qiáng)結(jié)構(gòu)的力熱環(huán)境適應(yīng)性。對于火箭結(jié)構(gòu)動(dòng)強(qiáng)度,結(jié)構(gòu)頻率需滿足總體頻率要求,響應(yīng)需滿足結(jié)構(gòu)設(shè)計(jì)要求,動(dòng)態(tài)位移需滿足結(jié)構(gòu)總體干涉和運(yùn)動(dòng)包絡(luò)要求。同時(shí),針對火箭發(fā)動(dòng)機(jī)采取的主要控制措施包括:

    1) 振動(dòng)傳遞路徑優(yōu)化。進(jìn)行振動(dòng)傳遞路徑分析(TPA)[85],分析振動(dòng)的主要傳遞路線及貢獻(xiàn)量,采取相應(yīng)的隔振/減振/抑振控制措施,以降低結(jié)構(gòu)振動(dòng)的輸入載荷及響應(yīng)水平,從而提高疲勞壽命,如進(jìn)行傳力路徑優(yōu)化[86]與采取復(fù)合材料機(jī)架、阻尼機(jī)架、隔/減振墊設(shè)計(jì)等。

    2) 結(jié)構(gòu)耦合共振檢測與頻率管理設(shè)計(jì)。振動(dòng)環(huán)境下疲勞失效的主要特征是共振疲勞及疲勞損傷累積。頻率管理的目的是避開結(jié)構(gòu)共振,通過對發(fā)動(dòng)機(jī)激勵(lì)源頻率特征提取、動(dòng)力學(xué)傳遞特性獲取、振動(dòng)頻率管理原則、改進(jìn)措施研究,以進(jìn)行頻率耦合檢查、頻率管理設(shè)計(jì)[87]及相關(guān)標(biāo)準(zhǔn)制定。此外,應(yīng)避開“有害”模態(tài)、共振或動(dòng)增益較大的振動(dòng)。

    3) 振動(dòng)響應(yīng)優(yōu)化。針對發(fā)動(dòng)機(jī)典型結(jié)構(gòu),通過增加卡箍、附加阻尼、改進(jìn)結(jié)構(gòu)設(shè)計(jì)等措施,以降低結(jié)構(gòu)的振動(dòng)響應(yīng)水平,從而提高結(jié)構(gòu)的疲勞抗力。

    5 總結(jié)與展望

    液體火箭發(fā)動(dòng)機(jī)結(jié)構(gòu)動(dòng)力學(xué)設(shè)計(jì)以極端環(huán)境下復(fù)雜結(jié)構(gòu)可能涉及到的各種動(dòng)力學(xué)問題為出發(fā)點(diǎn),通過研究諸多力學(xué)基礎(chǔ)與前沿技術(shù),在實(shí)現(xiàn)發(fā)動(dòng)機(jī)系統(tǒng)設(shè)計(jì)預(yù)定功能的前提下,使結(jié)構(gòu)設(shè)計(jì)合理、多功能一體化、極致輕量化、極限承載、長壽命與結(jié)構(gòu)完整性,從而確保發(fā)動(dòng)機(jī)安全可靠地工作。目前,發(fā)動(dòng)機(jī)結(jié)構(gòu)動(dòng)力學(xué)設(shè)計(jì)技術(shù)已成為學(xué)術(shù)界與工程界研究的熱點(diǎn)之一。迄今為止,國內(nèi)外學(xué)者在該領(lǐng)域已開展了大量的研究工作,進(jìn)行了材料/結(jié)構(gòu)的力學(xué)行為及本構(gòu)理論、載荷預(yù)計(jì)、動(dòng)力學(xué)建模分析、耦合損傷與動(dòng)態(tài)失效機(jī)制、強(qiáng)度評估及壽命預(yù)測、抗疲勞與延壽設(shè)計(jì)、動(dòng)力學(xué)試驗(yàn)與試驗(yàn)仿真等技術(shù)研究,取得了豐碩的成果,促進(jìn)了液體火箭發(fā)動(dòng)機(jī)結(jié)構(gòu)動(dòng)力學(xué)設(shè)計(jì)技術(shù)的發(fā)展。

    然而,由于液體火箭發(fā)動(dòng)機(jī)結(jié)構(gòu)動(dòng)力學(xué)是一個(gè)多學(xué)科交叉、龐大的系統(tǒng)性工程,動(dòng)力學(xué)問題非常復(fù)雜,其研究與解決的難度巨大。因此,在構(gòu)建涵蓋全壽命周期的液體火箭發(fā)動(dòng)機(jī)結(jié)構(gòu)動(dòng)力學(xué)設(shè)計(jì)技術(shù)體系時(shí),該研究方向的后續(xù)重點(diǎn)工作包括:

    1) 頂層策劃。從設(shè)計(jì)理念、基礎(chǔ)理論、技術(shù)方法及標(biāo)準(zhǔn)等方面開展研究,建立液體火箭發(fā)動(dòng)機(jī)動(dòng)強(qiáng)度設(shè)計(jì)、分析評估、試驗(yàn)及抗疲勞設(shè)計(jì)技術(shù)體系。

    2) 材料力學(xué)性能與疲勞性能數(shù)據(jù)??紤]材料在實(shí)際批次散差、熱處理、加工工藝及使用循環(huán)載荷影響下,建設(shè)發(fā)動(dòng)機(jī)材料基礎(chǔ)力學(xué)性能、疲勞性能數(shù)據(jù)庫,構(gòu)建材料力學(xué)性能與疲勞性能試驗(yàn)新方法與表征。

    3) 先進(jìn)測試分析技術(shù)。開展基于視覺觀測的全場非接觸先進(jìn)測試、智能感知技術(shù)研究,開展基于信息融合技術(shù)的大數(shù)據(jù)快速挖掘與數(shù)據(jù)特征準(zhǔn)確提取技術(shù)研究,發(fā)展復(fù)雜載荷準(zhǔn)確高效預(yù)計(jì)技術(shù)。

    4) 復(fù)雜結(jié)構(gòu)動(dòng)力學(xué)建模分析方法。建立大型復(fù)雜結(jié)構(gòu)多層多級、精準(zhǔn)、高分辨率模型,不斷提升新的面向真實(shí)服役場景的多場耦合分析能力,推進(jìn)高精度、高效及穩(wěn)定的數(shù)值計(jì)算方法。

    5) 動(dòng)強(qiáng)度評估與壽命評定技術(shù)。研究極端環(huán)境下結(jié)構(gòu)動(dòng)態(tài)力學(xué)行為與多態(tài)交互耦合損傷機(jī)理、動(dòng)態(tài)失效判據(jù)與評估方法,開展結(jié)構(gòu)疲勞壽命可靠性與概率壽命研究,建立發(fā)動(dòng)機(jī)結(jié)構(gòu)動(dòng)強(qiáng)度評估與壽命評定技術(shù)。

    6) 抗疲勞設(shè)計(jì)與延壽技術(shù)。從安全壽命設(shè)計(jì)到經(jīng)濟(jì)壽命設(shè)計(jì),開展新材料與智能制造、振動(dòng)主/被動(dòng)控制、結(jié)構(gòu)動(dòng)力學(xué)優(yōu)化設(shè)計(jì)、減損控制(DMC)與延壽技術(shù)研究,提高結(jié)構(gòu)的健壯性、動(dòng)強(qiáng)度可靠性與疲勞抗力。

    7) 靜動(dòng)態(tài)復(fù)合加載試驗(yàn)技術(shù)。考慮發(fā)動(dòng)機(jī)結(jié)構(gòu)實(shí)際工作環(huán)境載荷,采用載荷等效、結(jié)構(gòu)模擬等方法,推動(dòng)靜動(dòng)態(tài)載荷聯(lián)合加載試驗(yàn)技術(shù),對結(jié)構(gòu)動(dòng)強(qiáng)度與壽命進(jìn)行真實(shí)、準(zhǔn)確地等效考核與驗(yàn)證。

    8) 低成本、可重復(fù)使用液體火箭發(fā)動(dòng)機(jī)設(shè)計(jì)技術(shù)。隨著航天技術(shù)的飛速發(fā)展,各航天大國陸續(xù)開展了航天運(yùn)輸系統(tǒng)的升級換代,商業(yè)化驅(qū)動(dòng)下的低成本、可重復(fù)使用和航班化航天運(yùn)輸是世界各國關(guān)注的焦點(diǎn)及未來航天技術(shù)發(fā)展的主方向。

    猜你喜歡
    壽命裂紋動(dòng)力學(xué)
    《空氣動(dòng)力學(xué)學(xué)報(bào)》征稿簡則
    人類壽命極限應(yīng)在120~150歲之間
    中老年保健(2021年8期)2021-12-02 23:55:49
    倉鼠的壽命知多少
    Epidermal growth factor receptor rs17337023 polymorphism in hypertensive gestational diabetic women: A pilot study
    馬烈光養(yǎng)生之悟 自靜其心延壽命
    人類正常壽命為175歲
    奧秘(2017年12期)2017-07-04 11:37:14
    微裂紋區(qū)對主裂紋擴(kuò)展的影響
    基于隨機(jī)-動(dòng)力學(xué)模型的非均勻推移質(zhì)擴(kuò)散
    TNAE的合成和熱分解動(dòng)力學(xué)
    預(yù)裂紋混凝土拉壓疲勞荷載下裂紋擴(kuò)展速率
    女人高潮潮喷娇喘18禁视频| 成年女人毛片免费观看观看9| 精品国内亚洲2022精品成人| 国产无遮挡羞羞视频在线观看| 久久精品aⅴ一区二区三区四区| 嫩草影院精品99| 国产欧美日韩一区二区精品| 淫妇啪啪啪对白视频| 欧美日韩精品网址| av视频免费观看在线观看| 久久中文字幕一级| aaaaa片日本免费| 亚洲色图综合在线观看| 丝袜在线中文字幕| 19禁男女啪啪无遮挡网站| 国产精品成人在线| 久久国产精品人妻蜜桃| 身体一侧抽搐| 国产91精品成人一区二区三区| 黄频高清免费视频| 自拍欧美九色日韩亚洲蝌蚪91| 99久久99久久久精品蜜桃| 纯流量卡能插随身wifi吗| 精品乱码久久久久久99久播| 丝袜人妻中文字幕| 久久久久久久久免费视频了| 黄色 视频免费看| 国产91精品成人一区二区三区| 一级,二级,三级黄色视频| 成熟少妇高潮喷水视频| av免费在线观看网站| 国产真人三级小视频在线观看| 国产av一区在线观看免费| 欧美 亚洲 国产 日韩一| 日本五十路高清| 国产精品日韩av在线免费观看 | 日韩欧美在线二视频| 超碰成人久久| 在线观看一区二区三区| 免费看十八禁软件| 成人18禁在线播放| 亚洲中文字幕日韩| 美女高潮喷水抽搐中文字幕| 69精品国产乱码久久久| 丝袜美足系列| 精品久久久久久,| 丝袜美腿诱惑在线| 久久精品国产99精品国产亚洲性色 | 日本vs欧美在线观看视频| 免费一级毛片在线播放高清视频 | 99在线人妻在线中文字幕| 国产色视频综合| 最新美女视频免费是黄的| 变态另类成人亚洲欧美熟女 | 麻豆av在线久日| 多毛熟女@视频| 9热在线视频观看99| 99re在线观看精品视频| 亚洲一区二区三区欧美精品| www.熟女人妻精品国产| 麻豆久久精品国产亚洲av | 在线观看舔阴道视频| 神马国产精品三级电影在线观看 | 9热在线视频观看99| 日本 av在线| 免费人成视频x8x8入口观看| 制服人妻中文乱码| 悠悠久久av| 国产一区二区三区综合在线观看| 久久香蕉激情| 91九色精品人成在线观看| 人人妻,人人澡人人爽秒播| 久久精品人人爽人人爽视色| 中文字幕色久视频| 精品福利观看| 国产成年人精品一区二区 | 国产野战对白在线观看| 69av精品久久久久久| 精品午夜福利视频在线观看一区| 在线观看免费视频日本深夜| 在线天堂中文资源库| 在线观看免费午夜福利视频| 久久国产亚洲av麻豆专区| 三上悠亚av全集在线观看| 琪琪午夜伦伦电影理论片6080| 热99re8久久精品国产| 成人精品一区二区免费| cao死你这个sao货| 搡老乐熟女国产| 亚洲三区欧美一区| 90打野战视频偷拍视频| 大香蕉久久成人网| 亚洲欧美精品综合久久99| 久久精品国产综合久久久| 99国产极品粉嫩在线观看| 极品教师在线免费播放| 正在播放国产对白刺激| 嫩草影视91久久| 丝袜美腿诱惑在线| 久久午夜综合久久蜜桃| 久久人人97超碰香蕉20202| 国产精品爽爽va在线观看网站 | av在线天堂中文字幕 | 亚洲中文字幕日韩| 成在线人永久免费视频| 法律面前人人平等表现在哪些方面| av免费在线观看网站| 国产精品免费一区二区三区在线| 国产成人av教育| netflix在线观看网站| 国产精品久久久av美女十八| 一区二区三区精品91| 淫妇啪啪啪对白视频| 久久久久久人人人人人| 男女之事视频高清在线观看| av天堂久久9| 热re99久久精品国产66热6| 欧美日本中文国产一区发布| 国产单亲对白刺激| 亚洲 欧美 日韩 在线 免费| 国产精品自产拍在线观看55亚洲| 老司机午夜福利在线观看视频| 女性被躁到高潮视频| 男女做爰动态图高潮gif福利片 | 老司机福利观看| 水蜜桃什么品种好| 欧美日本中文国产一区发布| 人妻久久中文字幕网| 国产精品1区2区在线观看.| 亚洲激情在线av| 国产av又大| 看片在线看免费视频| 在线视频色国产色| 91国产中文字幕| 一级毛片高清免费大全| 男人舔女人下体高潮全视频| 欧美日韩视频精品一区| 咕卡用的链子| 久久天堂一区二区三区四区| 极品人妻少妇av视频| 超碰成人久久| 日韩三级视频一区二区三区| 人妻久久中文字幕网| 中文亚洲av片在线观看爽| 久久久水蜜桃国产精品网| 久久久精品国产亚洲av高清涩受| 国产精品99久久99久久久不卡| 国产精品免费一区二区三区在线| 亚洲一区高清亚洲精品| 国产一区二区三区视频了| 色综合站精品国产| 国产片内射在线| 色哟哟哟哟哟哟| 久久久久久大精品| 久久伊人香网站| 亚洲精品中文字幕在线视频| 狠狠狠狠99中文字幕| 国产高清国产精品国产三级| 亚洲精品在线观看二区| 18禁裸乳无遮挡免费网站照片 | 18禁黄网站禁片午夜丰满| 高清欧美精品videossex| 日韩欧美一区视频在线观看| 在线观看日韩欧美| 变态另类成人亚洲欧美熟女 | 国产成人av教育| 欧美日韩瑟瑟在线播放| 成年人免费黄色播放视频| 桃色一区二区三区在线观看| 亚洲成人免费电影在线观看| 香蕉丝袜av| 正在播放国产对白刺激| 自线自在国产av| 欧美成人性av电影在线观看| 亚洲精品av麻豆狂野| 精品久久久久久,| 99riav亚洲国产免费| 中文字幕高清在线视频| 精品国产美女av久久久久小说| 一区在线观看完整版| 午夜福利影视在线免费观看| 一二三四社区在线视频社区8| aaaaa片日本免费| 亚洲欧美精品综合一区二区三区| 一个人免费在线观看的高清视频| 多毛熟女@视频| 午夜精品国产一区二区电影| 中亚洲国语对白在线视频| 免费在线观看日本一区| 日本 av在线| 在线观看日韩欧美| 国产精品亚洲一级av第二区| 香蕉丝袜av| 精品第一国产精品| 久久狼人影院| 精品国产一区二区久久| 久久香蕉精品热| 黄色视频,在线免费观看| 欧美中文综合在线视频| 大型黄色视频在线免费观看| 久久久久九九精品影院| 搡老熟女国产l中国老女人| 丰满的人妻完整版| 女人被狂操c到高潮| 亚洲 欧美 日韩 在线 免费| 亚洲黑人精品在线| 黑人操中国人逼视频| 欧美黑人精品巨大| 女人被狂操c到高潮| 丰满饥渴人妻一区二区三| 天天添夜夜摸| 人妻丰满熟妇av一区二区三区| 亚洲自拍偷在线| 超色免费av| 女人爽到高潮嗷嗷叫在线视频| 欧美午夜高清在线| 亚洲欧美精品综合久久99| 亚洲国产欧美日韩在线播放| 国产有黄有色有爽视频| 久久国产亚洲av麻豆专区| 精品人妻在线不人妻| 精品国产一区二区久久| 亚洲精品中文字幕在线视频| 在线观看午夜福利视频| 岛国在线观看网站| 极品教师在线免费播放| 丁香六月欧美| 91麻豆精品激情在线观看国产 | 另类亚洲欧美激情| 少妇的丰满在线观看| 999久久久精品免费观看国产| 久久中文字幕一级| 在线观看免费高清a一片| 亚洲九九香蕉| 亚洲精品久久午夜乱码| 另类亚洲欧美激情| 高清av免费在线| 欧美日韩亚洲高清精品| 99riav亚洲国产免费| 高清欧美精品videossex| 长腿黑丝高跟| 国产aⅴ精品一区二区三区波| 亚洲精品一区av在线观看| 国产成人影院久久av| 一本综合久久免费| 一个人免费在线观看的高清视频| 最近最新中文字幕大全电影3 | 热99国产精品久久久久久7| 久久精品影院6| 久久久国产一区二区| 91字幕亚洲| 亚洲一区二区三区欧美精品| 手机成人av网站| 欧美日韩中文字幕国产精品一区二区三区 | 精品熟女少妇八av免费久了| 老熟妇乱子伦视频在线观看| 国产成人精品在线电影| 黄片播放在线免费| 18禁国产床啪视频网站| 亚洲男人的天堂狠狠| 国产深夜福利视频在线观看| 中文字幕人妻熟女乱码| 身体一侧抽搐| 成人手机av| 老司机亚洲免费影院| 长腿黑丝高跟| 亚洲在线自拍视频| 国产欧美日韩一区二区三| 久久草成人影院| 老熟妇仑乱视频hdxx| 18美女黄网站色大片免费观看| 黑人操中国人逼视频| 免费在线观看完整版高清| 成人三级做爰电影| 搡老岳熟女国产| 日本一区二区免费在线视频| 精品少妇一区二区三区视频日本电影| 久久精品国产亚洲av高清一级| 黄色毛片三级朝国网站| 成年人免费黄色播放视频| 亚洲欧洲精品一区二区精品久久久| 午夜影院日韩av| 一区二区三区精品91| 欧美不卡视频在线免费观看 | 欧美性长视频在线观看| 国产有黄有色有爽视频| 91九色精品人成在线观看| 99在线人妻在线中文字幕| 国产高清国产精品国产三级| 国产精品国产av在线观看| 国产单亲对白刺激| 久久国产乱子伦精品免费另类| 国产极品粉嫩免费观看在线| 老汉色∧v一级毛片| 电影成人av| 一本综合久久免费| 国产有黄有色有爽视频| 精品高清国产在线一区| 久久久国产精品麻豆| 精品国产乱子伦一区二区三区| 丰满的人妻完整版| 国产av一区在线观看免费| 美女午夜性视频免费| 亚洲成a人片在线一区二区| 少妇粗大呻吟视频| 久久午夜亚洲精品久久| 变态另类成人亚洲欧美熟女 | 欧美另类亚洲清纯唯美| 一级a爱片免费观看的视频| 精品无人区乱码1区二区| 亚洲一区二区三区色噜噜 | 国产精品久久久久成人av| 国产精品日韩av在线免费观看 | 精品少妇一区二区三区视频日本电影| 51午夜福利影视在线观看| 久久久久久久久中文| 免费av毛片视频| 可以免费在线观看a视频的电影网站| 99久久人妻综合| 男人操女人黄网站| 欧美日韩黄片免| 国产黄色免费在线视频| 水蜜桃什么品种好| 一个人观看的视频www高清免费观看 | 女人被狂操c到高潮| 嫩草影视91久久| av福利片在线| 久久久久久人人人人人| 日本一区二区免费在线视频| 天堂中文最新版在线下载| xxxhd国产人妻xxx| 国产精品一区二区三区四区久久 | 大香蕉久久成人网| 日本黄色视频三级网站网址| 国产精品偷伦视频观看了| 国产乱人伦免费视频| 午夜亚洲福利在线播放| 久久久久久久久中文| 美女高潮到喷水免费观看| 国产日韩一区二区三区精品不卡| 18禁国产床啪视频网站| 91大片在线观看| 亚洲国产看品久久| 精品久久蜜臀av无| av在线播放免费不卡| 精品午夜福利视频在线观看一区| 国产精品香港三级国产av潘金莲| 人妻久久中文字幕网| 巨乳人妻的诱惑在线观看| 中文字幕色久视频| 天天影视国产精品| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲伊人色综图| 黄色片一级片一级黄色片| 老司机靠b影院| 久久九九热精品免费| 久久精品91蜜桃| 大香蕉久久成人网| 欧美日韩乱码在线| 亚洲中文日韩欧美视频| 国产免费现黄频在线看| 国产熟女xx| 叶爱在线成人免费视频播放| 亚洲va日本ⅴa欧美va伊人久久| 极品教师在线免费播放| 两性午夜刺激爽爽歪歪视频在线观看 | 久久人人爽av亚洲精品天堂| 深夜精品福利| 人妻丰满熟妇av一区二区三区| 日韩欧美三级三区| 夜夜躁狠狠躁天天躁| 十八禁网站免费在线| 亚洲精品久久午夜乱码| 国产精品日韩av在线免费观看 | 午夜视频精品福利| 日韩有码中文字幕| 亚洲精品一区av在线观看| 一区二区三区精品91| 亚洲专区字幕在线| 男女之事视频高清在线观看| 高清欧美精品videossex| 亚洲熟妇中文字幕五十中出 | 国产一区二区在线av高清观看| 一级毛片女人18水好多| 日韩三级视频一区二区三区| 亚洲av成人一区二区三| 中文字幕另类日韩欧美亚洲嫩草| 757午夜福利合集在线观看| 大型av网站在线播放| 岛国在线观看网站| 又大又爽又粗| 亚洲激情在线av| 级片在线观看| 免费高清视频大片| 国产成人免费无遮挡视频| 国产一区二区在线av高清观看| 国产精品野战在线观看 | 9热在线视频观看99| 最近最新免费中文字幕在线| 激情视频va一区二区三区| 男人操女人黄网站| 99久久人妻综合| 亚洲免费av在线视频| 美女大奶头视频| 熟女少妇亚洲综合色aaa.| 在线观看免费视频网站a站| 国产成人啪精品午夜网站| 一区二区日韩欧美中文字幕| 国产人伦9x9x在线观看| 99在线人妻在线中文字幕| 一区二区三区国产精品乱码| 丝袜美腿诱惑在线| 老司机在亚洲福利影院| 黄片大片在线免费观看| 动漫黄色视频在线观看| 久久久国产成人免费| 欧美日韩亚洲高清精品| 国产亚洲欧美98| 在线视频色国产色| 男女高潮啪啪啪动态图| 深夜精品福利| 看黄色毛片网站| 亚洲性夜色夜夜综合| 久久精品亚洲熟妇少妇任你| 在线观看午夜福利视频| 欧美日韩亚洲高清精品| 十分钟在线观看高清视频www| 日韩成人在线观看一区二区三区| 亚洲精品国产一区二区精华液| 欧美黑人精品巨大| 50天的宝宝边吃奶边哭怎么回事| 视频区欧美日本亚洲| 一本综合久久免费| 最近最新免费中文字幕在线| 999精品在线视频| 久99久视频精品免费| 最好的美女福利视频网| 人妻久久中文字幕网| 他把我摸到了高潮在线观看| 黄色怎么调成土黄色| 91成年电影在线观看| 亚洲专区字幕在线| 18美女黄网站色大片免费观看| 国产无遮挡羞羞视频在线观看| 日本一区二区免费在线视频| 搡老岳熟女国产| 国产高清videossex| 88av欧美| 精品国产国语对白av| 免费女性裸体啪啪无遮挡网站| 九色亚洲精品在线播放| 久久香蕉激情| 欧美av亚洲av综合av国产av| 亚洲欧洲精品一区二区精品久久久| 一二三四社区在线视频社区8| 午夜日韩欧美国产| 亚洲av熟女| 又紧又爽又黄一区二区| 在线看a的网站| 久热爱精品视频在线9| 亚洲色图 男人天堂 中文字幕| 国产深夜福利视频在线观看| 美女扒开内裤让男人捅视频| 欧美黑人欧美精品刺激| 久久久久久久久中文| 91av网站免费观看| 亚洲精品成人av观看孕妇| 一级作爱视频免费观看| 亚洲成人免费av在线播放| 多毛熟女@视频| 国产欧美日韩精品亚洲av| 欧美日韩国产mv在线观看视频| 久久精品国产亚洲av高清一级| 精品少妇一区二区三区视频日本电影| 国产精品美女特级片免费视频播放器 | av视频免费观看在线观看| 人妻丰满熟妇av一区二区三区| 久久久精品欧美日韩精品| 老司机亚洲免费影院| 桃红色精品国产亚洲av| 成人国产一区最新在线观看| 亚洲精品国产一区二区精华液| 岛国视频午夜一区免费看| 在线国产一区二区在线| 亚洲国产欧美一区二区综合| 国产精品综合久久久久久久免费 | 国产精品一区二区免费欧美| 午夜福利免费观看在线| 水蜜桃什么品种好| 一级作爱视频免费观看| 在线观看一区二区三区激情| 久久久久久久久免费视频了| 欧美日韩亚洲综合一区二区三区_| 欧美不卡视频在线免费观看 | 亚洲国产精品sss在线观看 | 久久精品国产亚洲av香蕉五月| 人妻丰满熟妇av一区二区三区| 久久中文字幕一级| 亚洲中文字幕日韩| 桃色一区二区三区在线观看| 午夜免费观看网址| 国产亚洲精品综合一区在线观看 | 天天躁夜夜躁狠狠躁躁| 国产亚洲欧美精品永久| 欧美日本中文国产一区发布| 欧美精品啪啪一区二区三区| 成人黄色视频免费在线看| 亚洲成人精品中文字幕电影 | 一a级毛片在线观看| 午夜免费鲁丝| 亚洲欧美一区二区三区黑人| 亚洲精品中文字幕在线视频| 国产单亲对白刺激| 亚洲成av片中文字幕在线观看| 9191精品国产免费久久| 国产精品 国内视频| 天天躁夜夜躁狠狠躁躁| 午夜久久久在线观看| 亚洲精品美女久久久久99蜜臀| 91九色精品人成在线观看| 在线观看午夜福利视频| av中文乱码字幕在线| 中文字幕最新亚洲高清| 99国产综合亚洲精品| 亚洲精品一卡2卡三卡4卡5卡| 99热国产这里只有精品6| 国产精品香港三级国产av潘金莲| 亚洲精品av麻豆狂野| 高清黄色对白视频在线免费看| 色哟哟哟哟哟哟| 精品欧美一区二区三区在线| 久久精品影院6| 欧美日韩瑟瑟在线播放| av福利片在线| 亚洲片人在线观看| 看黄色毛片网站| 亚洲一区二区三区欧美精品| 神马国产精品三级电影在线观看 | 成人亚洲精品一区在线观看| 国产精品一区二区精品视频观看| 国产成人一区二区三区免费视频网站| av视频免费观看在线观看| 男男h啪啪无遮挡| 久久亚洲真实| 神马国产精品三级电影在线观看 | 90打野战视频偷拍视频| 亚洲欧美日韩另类电影网站| 久久这里只有精品19| 两个人免费观看高清视频| 黄色毛片三级朝国网站| a级毛片在线看网站| 天天躁夜夜躁狠狠躁躁| 国产片内射在线| 90打野战视频偷拍视频| 在线观看66精品国产| 日本vs欧美在线观看视频| 国产精品国产av在线观看| 欧美乱码精品一区二区三区| 高清毛片免费观看视频网站 | 十八禁人妻一区二区| 亚洲成人免费av在线播放| 黄色a级毛片大全视频| 日本a在线网址| 国产欧美日韩一区二区三区在线| svipshipincom国产片| 国产一区二区三区在线臀色熟女 | 久久久精品国产亚洲av高清涩受| 91麻豆精品激情在线观看国产 | 久久精品成人免费网站| av天堂在线播放| 久久人妻av系列| 亚洲九九香蕉| 老司机靠b影院| 夫妻午夜视频| 一个人免费在线观看的高清视频| 日韩中文字幕欧美一区二区| 99久久综合精品五月天人人| 亚洲精品美女久久av网站| a级毛片黄视频| 精品无人区乱码1区二区| 黄片小视频在线播放| 欧美国产精品va在线观看不卡| 一区在线观看完整版| 精品人妻在线不人妻| 搡老熟女国产l中国老女人| 制服诱惑二区| 日韩欧美一区视频在线观看| 亚洲欧美一区二区三区久久| 午夜精品在线福利| 在线观看免费视频网站a站| 黑人操中国人逼视频| tocl精华| 精品人妻在线不人妻| av电影中文网址| 中文字幕人妻熟女乱码| 国产欧美日韩一区二区三| 两性午夜刺激爽爽歪歪视频在线观看 | 少妇 在线观看| 操美女的视频在线观看| 十八禁网站免费在线| 国产精品98久久久久久宅男小说| 18禁美女被吸乳视频| 久久人人97超碰香蕉20202| 国产aⅴ精品一区二区三区波| 狠狠狠狠99中文字幕| 99久久人妻综合| 人人妻人人添人人爽欧美一区卜| 久久亚洲真实| 波多野结衣av一区二区av| av国产精品久久久久影院| 久久精品成人免费网站| 天天影视国产精品| 午夜日韩欧美国产| 后天国语完整版免费观看| 亚洲精品久久成人aⅴ小说|