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

    水反應(yīng)金屬燃料發(fā)動(dòng)機(jī)藥柱結(jié)構(gòu)完整性分析

    2022-01-15 08:18:16王志龍韓新波尹韶平
    關(guān)鍵詞:擋塊藥柱本構(gòu)

    王志龍, 韓新波, 喬 宏, 尹韶平, 韓 艾, 雷 鳴

    水反應(yīng)金屬燃料發(fā)動(dòng)機(jī)藥柱結(jié)構(gòu)完整性分析

    王志龍1, 韓新波1, 喬 宏1, 尹韶平1, 韓 艾2, 雷 鳴2

    (1. 中國(guó)船舶集團(tuán)有限公司 第705研究所, 陜西西安, 710077; 2. 西北工業(yè)大學(xué) 航天學(xué)院, 陜西西安, 710072)

    為研究魚雷水反應(yīng)金屬燃料發(fā)動(dòng)機(jī)高夾雜比燃料藥柱在局部限位條件下的結(jié)構(gòu)完整性問題, 開展了燃料藥柱試樣力學(xué)性能試驗(yàn), 獲得了藥柱力學(xué)性能參數(shù)。在有限元軟件中建立了燃料藥柱試樣的數(shù)值分析模型, 進(jìn)行了藥柱試樣的力學(xué)性能仿真, 通過仿真結(jié)果和試驗(yàn)結(jié)果的對(duì)比, 驗(yàn)證了所建立的藥柱力學(xué)本構(gòu)模型的合理性。利用該力學(xué)本構(gòu)模型對(duì)水反應(yīng)金屬燃料發(fā)動(dòng)機(jī)中燃料藥柱的結(jié)構(gòu)完整性進(jìn)行了仿真分析,得到了燃料藥柱的力-位移曲線、應(yīng)力分布云圖以及位移云圖,仿真結(jié)果表明, 燃料藥柱在典型工況下未發(fā)生破壞,徑向和軸向的變形程度符合發(fā)動(dòng)機(jī)整體結(jié)構(gòu)要求。

    魚雷; 水反應(yīng)金屬燃料發(fā)動(dòng)機(jī); 燃料藥柱結(jié)構(gòu)完整性; 有限元分析

    0 引言

    動(dòng)力系統(tǒng)是魚雷的主要組成部分, 它對(duì)魚雷的航速、航程、航深、可靠性及成本起著決定性作用[1]。由于海水具有高密度特性, 魚雷受到的阻力為空中飛行器的800多倍, 使得航速和航程成為常規(guī)動(dòng)力魚雷技術(shù)發(fā)展的“軟肋”[2]。水反應(yīng)金屬燃料發(fā)動(dòng)機(jī)是超高速魚雷的主發(fā)動(dòng)機(jī), 具有比沖高、推力大的特點(diǎn)[3], 可以有效提高魚雷的航速[4]。水反應(yīng)金屬燃料發(fā)動(dòng)機(jī)所使用的燃料不同于固體火箭發(fā)動(dòng)機(jī)的推進(jìn)劑[5], 該燃料金屬含量極高, 氧化劑和粘合劑含量很少, 粘性特性不明顯[6], 屬于高夾雜比貧氧燃料, 成型方法與傳統(tǒng)的復(fù)合推進(jìn)劑以及雙基推進(jìn)劑均不相同, 力學(xué)性能目前尚未有深入研究。

    在水反應(yīng)金屬燃料發(fā)動(dòng)機(jī)中, 燃料藥柱受力條件惡劣, 不僅要承受點(diǎn)火時(shí)的壓強(qiáng)載荷, 還要承受活塞和局部限位擋塊的擠壓作用。設(shè)計(jì)時(shí)需要對(duì)燃料藥柱進(jìn)行結(jié)構(gòu)完整性分析[7], 計(jì)算燃料藥柱在工作載荷下的力學(xué)響應(yīng), 獲取藥柱發(fā)生破壞時(shí)的載荷極限以及徑向位移量, 評(píng)價(jià)其工作的可靠性, 保證其內(nèi)部或表面不產(chǎn)生裂紋破壞、不出現(xiàn)藥柱與發(fā)動(dòng)機(jī)殼體接觸等現(xiàn)象。若出現(xiàn)裂紋破壞, 燃料藥柱的燃面會(huì)突然增大, 導(dǎo)致發(fā)動(dòng)機(jī)的內(nèi)壓急劇上升, 造成推力特性畸形甚至爆炸等嚴(yán)重后果[8]。

    燃料藥柱力學(xué)性能和本構(gòu)模型的研究是發(fā)動(dòng)機(jī)藥柱結(jié)構(gòu)完整性評(píng)估的基礎(chǔ)。目前, 國(guó)內(nèi)外學(xué)者對(duì)傳統(tǒng)固體火箭發(fā)動(dòng)機(jī)藥柱結(jié)構(gòu)完整性分析進(jìn)行了大量的研究, 此類藥柱由于粘合劑組分含量高, 呈現(xiàn)出明顯的粘彈塑特性[9]。申志彬[10]、李磊[11]等采用三維粘彈性有限元法分析了傘盤結(jié)構(gòu)位置對(duì)燃料藥柱結(jié)構(gòu)完整性的影響; 檀葉[12]、溫瑞珩[13]等基于粘彈性本構(gòu)關(guān)系分析得到了套管裝藥發(fā)動(dòng)機(jī)工作時(shí)最危險(xiǎn)部位的應(yīng)力分布; Chyuan[14-15]采用粘彈性模型分析了均布內(nèi)壓載荷作用下, 材料參數(shù)對(duì)燃料藥柱結(jié)構(gòu)完整性的影響。而對(duì)于高夾雜比貧氧燃料藥柱結(jié)構(gòu)完整性分析目前卻少有公開報(bào)道。

    文中針對(duì)高夾雜比貧氧燃料藥柱結(jié)構(gòu)完整性問題, 開展力學(xué)性能試驗(yàn), 采用數(shù)字圖像相關(guān)法(digital image correlation, DIC)技術(shù)進(jìn)行了燃料藥柱試樣損傷狀態(tài)監(jiān)測(cè), 利用商用有限元軟件對(duì)發(fā)動(dòng)機(jī)中的燃料藥柱結(jié)構(gòu)完整性進(jìn)行了仿真分析, 得到典型工況下燃料藥柱的力-位移曲線、應(yīng)力云圖以及位移云圖, 對(duì)燃料藥柱的結(jié)構(gòu)完整性進(jìn)行了評(píng)估。

    1 推進(jìn)劑材料力學(xué)本構(gòu)模型

    1.1 彈-塑-脆性力學(xué)本構(gòu)模型

    圖1 彈-塑-脆性本構(gòu)模型

    1.2 剪切損傷理論

    其中

    2 抗壓試驗(yàn)與驗(yàn)證

    2.1 試驗(yàn)方案與結(jié)果分析

    參照GJB770B-2005《火藥試驗(yàn)方法》中的抗壓強(qiáng)度壓縮法, 采用萬(wàn)能試驗(yàn)機(jī)以3 mm/min的恒定加載速率對(duì)30 mm×30 mm的燃料藥柱試樣施加載荷, 在加載過程中為了保證燃料藥柱試樣所受壓力方向與其幾何中心線同軸, 需要確保萬(wàn)能試驗(yàn)機(jī)壓頭、燃料藥柱試樣以及放置試樣的圓盤同心。靜態(tài)壓縮試驗(yàn)裝置與燃料藥柱試樣如圖2所示。

    圖2 燃料藥柱壓縮試樣及試驗(yàn)裝置

    通過計(jì)算機(jī)采集并記錄燃料藥柱試樣在壓縮過程中所受壓力以及軸向位移的大小, 獲得燃料藥柱試樣一維軸向的力-位移曲線, 如圖3所示。

    圖3 試驗(yàn)測(cè)量的力-位移曲線

    從圖3可以看出, 在燃料藥柱試樣加載的初始階段, 材料尚處于彈性變形階段, 力-位移曲線基本呈線性變化趨勢(shì); 此后, 隨著載荷量的增大, 材料發(fā)生了塑性變形, 力-位移曲線也呈現(xiàn)非線性變化的趨勢(shì), 一直持續(xù)到該燃料藥柱試樣所能承受的最大壓力43 kN, 此時(shí)萬(wàn)能試驗(yàn)機(jī)壓縮裝置的位移量約為1.1 mm。若繼續(xù)對(duì)燃料藥柱試樣施加位移載荷, 力-位移曲線就呈下降趨勢(shì), 這表明燃料藥柱試樣此時(shí)已經(jīng)開始發(fā)生損傷破壞, 其內(nèi)部或表面已產(chǎn)生裂紋。

    2.2 材料參數(shù)

    燃料藥柱試樣為各向同性的均質(zhì)材料, 通過力學(xué)性能試驗(yàn)和DIC技術(shù), 獲得的藥柱材料參數(shù)如表1所示。

    表1 材料參數(shù)列表

    2.3 力學(xué)本構(gòu)模型驗(yàn)證

    為了保證燃料藥柱力學(xué)本構(gòu)模型選取的合理性, 首先需要采用該模型對(duì)燃料藥柱試樣進(jìn)行仿真計(jì)算, 并比較試樣的仿真結(jié)果和試驗(yàn)結(jié)果, 驗(yàn)證力學(xué)本構(gòu)模型的合理性。在有限元分析軟件中建立燃料藥柱試樣受壓的物理模型, 并對(duì)其進(jìn)行網(wǎng)格劃分, 如圖4所示。將萬(wàn)能試驗(yàn)機(jī)上下壓頭視為剛體部件。燃料藥柱試樣共劃分為71 036個(gè)六面體單元。

    圖4 燃料藥柱試樣三維有限元模型

    2.4 結(jié)果對(duì)比

    將燃料藥柱試樣試驗(yàn)得到的力-位移曲線和基于本構(gòu)模型仿真得到的力-位移曲線進(jìn)行對(duì)比, 如圖5所示??梢钥闯? 在壓縮速率恒定的條件下, 仿真和試驗(yàn)結(jié)果在燃料藥柱試樣發(fā)生損傷前的彈塑性階段比較吻合, 均在點(diǎn)達(dá)到了最大承載壓力, 約43 kN; 在點(diǎn)以后, 由于燃料藥柱試樣開始發(fā)生損傷, 力-位移曲線呈現(xiàn)明顯的下降趨勢(shì)。為了保證水反應(yīng)金屬燃料發(fā)動(dòng)機(jī)工作的可靠性, 需要確保燃料藥柱所承受的最大壓力不大于點(diǎn)處的壓力值, 下降階段并不代表研究中燃料藥柱的主要力學(xué)特性, 所以可不考慮點(diǎn)以后燃料藥柱的受力情況。從整體趨勢(shì)來看, 仿真和試驗(yàn)結(jié)果基本吻合, 說明了本構(gòu)模型的適用性。

    圖5 試驗(yàn)和本構(gòu)模型仿真結(jié)果對(duì)比

    圖6為DIC檢測(cè)得到的燃料藥柱試樣破壞模式, 可以看出, 燃料藥柱試樣中間部位的應(yīng)力值明顯低于其他部位, 剪切帶方向與水平方向間的夾角約為45°。

    圖7展示了燃料藥柱試樣在壓縮過程中4個(gè)典型時(shí)刻的應(yīng)力分布云圖。通過對(duì)比圖7(a)和(b)可看出, 在以3 mm/min的恒定加載速率對(duì)燃料藥柱進(jìn)行加載時(shí), 當(dāng)加載時(shí)間達(dá)到32.42 s, 燃料藥柱試樣中間部位的應(yīng)力值出現(xiàn)明顯下降, 表明該部位已發(fā)生損傷。根據(jù)損傷帶的形狀可判定燃料藥柱試樣的破壞形式為剪切損傷; 若繼續(xù)施加載荷, 當(dāng)加載時(shí)間為32.47 s時(shí), 由圖7(c)可看出, 此時(shí)燃料藥柱試樣中間一部分單元已被刪除, 出現(xiàn)了明顯的宏觀裂紋, 該位置的應(yīng)力值也因此急劇下降, 裂紋方向與水平方向間的夾角約45°; 為了研究損傷的發(fā)展, 進(jìn)一步加載, 由圖7(d)可看出, 圖7(c)中出現(xiàn)的裂紋沿著45°方向繼續(xù)擴(kuò)展, 直至整個(gè)燃料藥柱試樣完全破壞, 這與圖6的試驗(yàn)結(jié)果一致。綜上可以證明燃料藥柱力學(xué)本構(gòu)模型選擇是合理的。

    圖6 燃料藥柱試樣破壞模式

    圖7 加載過程中不同時(shí)刻應(yīng)力云圖

    3 發(fā)動(dòng)機(jī)內(nèi)燃料藥柱結(jié)構(gòu)完整性分析

    3.1 幾何模型與網(wǎng)格劃分

    采用上述力學(xué)本構(gòu)模型對(duì)210 mm口徑水反應(yīng)金屬燃料發(fā)動(dòng)機(jī)中的燃料藥柱結(jié)構(gòu)完整性進(jìn)行仿真分析。燃料藥柱置于活塞和限位擋塊之間, 一端受到活塞的推動(dòng), 另一端局部受到限位擋塊的擠壓, 限位擋塊為8個(gè)大小相同的圓柱體結(jié)構(gòu), 其余部分受到燃?xì)鈮簭?qiáng)作用。由于燃料藥柱是軸對(duì)稱結(jié)構(gòu), 為了提高網(wǎng)格劃分質(zhì)量, 加快仿真計(jì)算速度, 根據(jù)對(duì)稱性準(zhǔn)則, 把燃料藥柱均勻地分成完全對(duì)稱的8等份, 每一份的夾角為45°, 取其中的一份建立三維有限元計(jì)算模型, 并在模型的一側(cè)施加對(duì)稱邊界條件。對(duì)整個(gè)模型進(jìn)行網(wǎng)格劃分, 考慮到限位擋塊與燃料藥柱的接觸, 對(duì)該區(qū)域附近劃分了更細(xì)的網(wǎng)格。建立了燃料藥柱三維局部有限元模型, 燃料藥柱劃分成453 219個(gè)六面體單元, 如圖8所示。

    圖8 燃料藥柱三維局部有限元模型

    3.2 載荷與邊界條件

    由于限位擋塊與發(fā)動(dòng)機(jī)殼體相連, 且擋塊的剛度遠(yuǎn)遠(yuǎn)大于燃料藥柱的剛度, 在仿真中給擋塊施加剛體約束和完全固定約束, 限制其6個(gè)方向的自由度。發(fā)動(dòng)機(jī)工作時(shí), 在流體動(dòng)壓作用下, 活塞一端會(huì)受到水的壓強(qiáng)作用, 達(dá)到額定工況后, 活塞對(duì)燃料藥柱的作用力為101.8 kN, 同時(shí)燃料藥柱其他方向受到3 MPa的燃?xì)鈮簭?qiáng)。針對(duì)此工況, 對(duì)活塞一側(cè)的燃料藥柱端面施加大小為101.8 kN的壓力, 對(duì)擋塊一側(cè)的端面和燃料藥柱周向施加大小為3 MPa的均布?jí)簭?qiáng)。

    3.3 計(jì)算結(jié)果分析

    通過仿真計(jì)算, 得到了燃料藥柱的力學(xué)響應(yīng)特性。圖9反映了燃料藥柱在整個(gè)加載過程中的力-位移變化情況, 可以看出該曲線尚未出現(xiàn)下降階段, 說明燃料藥柱在該典型工況下只發(fā)生了彈塑性變形, 其內(nèi)部或表面沒有產(chǎn)生裂紋破壞, 滿足結(jié)構(gòu)完整性要求。

    圖9 燃料藥柱力-位移曲線

    圖10反映了燃料藥柱的應(yīng)力分布情況, 從圖中可以明顯看出, 在燃料藥柱與擋塊接觸區(qū)域會(huì)產(chǎn)生應(yīng)力集中現(xiàn)象, 且燃料藥柱與擋塊接觸區(qū)域內(nèi)側(cè)邊緣處的應(yīng)力值大于外側(cè)邊緣處的應(yīng)力值。為減小此處的應(yīng)力值, 在發(fā)動(dòng)機(jī)結(jié)構(gòu)允許的條件下, 可增加擋塊與燃料藥柱接觸面積。在軸向上, 活塞附近區(qū)域的燃料藥柱和擋塊附近區(qū)域的燃料藥柱所受到的應(yīng)力均比燃料藥柱中間部位所受到的應(yīng)力大, 活塞附近區(qū)域的燃料藥柱應(yīng)力值大約為1 MPa, 擋塊附近區(qū)域的燃料藥柱應(yīng)力值大約為2~4 MPa, 這主要是由于擋塊附近區(qū)域的燃料藥柱存在應(yīng)力集中。

    圖11和圖12分別反映了燃料藥柱徑向和軸向的變形情況, 由圖11可知燃料藥柱由于周向受到了燃?xì)獾膲簭?qiáng)作用, 整體的徑向尺寸變小, 在遠(yuǎn)離擋塊附近區(qū)域的燃料藥柱均處于收縮狀態(tài), 最大徑向收縮量為0.078 3 mm。與燃料藥柱所受最大應(yīng)力位置相同, 燃料藥柱的徑向最大擴(kuò)張量出現(xiàn)在了其與擋塊接觸區(qū)域內(nèi)側(cè)邊緣處, 大小為0.0271 mm, 遠(yuǎn)小于燃料藥柱和發(fā)動(dòng)機(jī)殼體內(nèi)表面的間隙, 兩者不會(huì)發(fā)生接觸; 由圖12可以看出燃料藥柱由于受到活塞和擋塊的擠壓作用, 在軸向上整體處于壓縮狀態(tài), 在燃料藥柱與擋塊接觸區(qū)域, 燃料藥柱軸向位移量達(dá)到最大值0.023 mm, 相對(duì)于燃料藥柱實(shí)際的運(yùn)動(dòng)路程可忽略不計(jì)。綜上所述, 在此典型工況下, 該燃料藥柱滿足結(jié)構(gòu)完整性要求。

    圖10 燃料藥柱應(yīng)力分布云圖

    圖11 燃料藥柱徑向位移分布云圖

    圖12 燃料藥柱軸向位移分布云圖

    結(jié)束語(yǔ)

    文中對(duì)燃料藥柱試樣進(jìn)行了壓縮試驗(yàn)和DIC損傷監(jiān)測(cè), 獲得了材料的基本參數(shù)以及破壞模式; 根據(jù)所選的力學(xué)本構(gòu)模型對(duì)燃料藥柱試樣的壓縮過程進(jìn)行了仿真計(jì)算, 通過試驗(yàn)結(jié)果和仿真結(jié)果的對(duì)比, 驗(yàn)證了燃料藥柱力學(xué)本構(gòu)模型選擇的合理性, 為水反應(yīng)金屬燃料發(fā)動(dòng)機(jī)內(nèi)燃料藥柱結(jié)構(gòu)完整性分析提供了理論依據(jù)。

    對(duì)典型工況下燃料藥柱結(jié)構(gòu)完整性進(jìn)行了分析, 仿真結(jié)果表明燃料藥柱在該典型工況下不會(huì)發(fā)生破壞, 燃料藥柱徑向和軸向的變形程度符合發(fā)動(dòng)機(jī)整體結(jié)構(gòu)要求。由于應(yīng)力集中現(xiàn)象的存在, 燃料藥柱最大應(yīng)力和最大徑向位移均出現(xiàn)在了燃料藥柱與擋塊接觸區(qū)域內(nèi)側(cè)邊緣處。

    [1] 查志武, 史小鋒, 錢志博. 魚雷熱動(dòng)力技術(shù)[M]. 北京: 國(guó)防工業(yè)出版社, 2006.

    [2] 張亞俊, 王祎, 李吉禎, 等. 高速魚雷水沖壓發(fā)動(dòng)機(jī)用金屬/水反應(yīng)燃料研究進(jìn)展[J]. 四川兵工學(xué)報(bào), 2013, 34 (5): 123-126.

    Zhang Ya-jun, Wang Yi, Li Ji-zhen, et al. Review on the Hydroreactive Metal Fuel Used in Water Ramjet Engine of High Speed Torpedo[J]. Journal of Sichuan Ordnance, 2013, 34(5): 123-126.

    [3] 晁侃, 牛楠, 陸賀建. 一次水燃比對(duì)高金屬含量鎂基推進(jìn)劑水沖壓發(fā)動(dòng)機(jī)比沖性能影響分析[J]. 水下無(wú)人系統(tǒng)學(xué)報(bào), 2017, 25(1): 52-56.

    Chao Kan, Niu Nan, Lu He-jian. Effect of Primary Water-to-fuel Ratio on Specific Impulse Performance of Water Ramjet with High Metal Content Propellant[J]. Journal of Unmanned Undersea Systems, 2017, 25(1): 52-56.

    [4] 黃海濤. 高能鎂基水反應(yīng)金屬燃料推進(jìn)劑研究[D]. 北京: 北京理工大學(xué), 2015.

    [5] Ozerov E S, Yurinov A A. Combustion of Particlesof Aluminum-magnesium Alloys in Water Vapor[J]. Com- bustion Explosion and Shock Waves, 1977, 13(6): 778- 780.

    [6] 趙衛(wèi)兵, 史小鋒, 伊寅, 等. 水反應(yīng)金屬燃料在超高速魚雷推進(jìn)系統(tǒng)中的應(yīng)用[J]. 火炸藥學(xué)報(bào), 2006(5): 53-56.

    Zhao Wei-bing, Shi Xiao-feng, Yi Yin, et al. Application of Hydroreactive Metal Fuel in Super-Cavitation Torpedo Propulsion System[J]. Chinese Journal of Explosives & Propellants, 2006(5): 53-56.

    [7] 鮑福廷, 侯曉. 固體火箭發(fā)動(dòng)機(jī)設(shè)計(jì)[M]. 第1版. 北京: 中國(guó)宇航出版社, 2016.

    [8] 宋仕雄. 低溫點(diǎn)火狀態(tài)下固體發(fā)動(dòng)機(jī)藥柱結(jié)構(gòu)完整性分析[D]. 西安: 航天動(dòng)力技術(shù)研究院, 2018.

    [9] 王哲君, 強(qiáng)洪夫, 王廣, 等. 固體推進(jìn)劑力學(xué)性能和本構(gòu)模型的研究進(jìn)展[J]. 含能材料, 2016, 24(4): 403-416.

    Wang Zhe-jun, Qiang Hong-fu, Wang Guang, et al. Review on the Mechanical Properties and Constitutive Models of Solid Propellants[J]. Chinese Journal of Energetic Materials,2016, 24(4): 403-416.

    [10] 申志彬, 姜人偉, 職世君. 傘盤結(jié)構(gòu)對(duì)大長(zhǎng)徑比發(fā)動(dòng)機(jī)藥柱結(jié)構(gòu)完整性的影響[J]. 國(guó)防科技大學(xué)學(xué)報(bào), 2017, 39(1): 148-152.

    Shen Zhi-bin, Jiang Ren-wei, Zhi Shi-jun. Effects of Umbrella Slot on the Structure Integrity of Solid Rocket Motor Grain with Large Aspect Ratio[J]. Journal of National University of Defense Technology,2017, 39(1): 148-152.

    [11] 李磊, 唐國(guó)金, 雷勇軍, 等. 固體火箭發(fā)動(dòng)機(jī)藥柱傘盤結(jié)構(gòu)應(yīng)力應(yīng)變分析[J]. 推進(jìn)技術(shù), 2008(4): 477-480, 507.

    Li Lei, Tang Guo-jin, Lei Yong-jun, et al. Stress and Strain Analysis for the Umbrella Slot Configuration of Solid Rocketmotor Grain[J]. Tactical Missile Technology, 2008(4): 477-480, 507.

    [12] 檀葉, 陳科. 套管裝藥發(fā)動(dòng)機(jī)藥柱結(jié)構(gòu)完整性分析[J].戰(zhàn)術(shù)導(dǎo)彈技術(shù), 2018(5): 90-94.

    Tan Ye, Chen Ke. Analysis of Structural Integrity for the Canular Solid Rocket Motor Grainl[J]. Tactical Missile Technology, 2018(5): 90-94.

    [13] 溫瑞珩, 李健. 固體火箭發(fā)動(dòng)機(jī)結(jié)構(gòu)破壞分析[J]. 戰(zhàn)術(shù)導(dǎo)彈技術(shù), 2019(3): 88-93.

    Wen Rui-heng, Li Jian. Structural Failure Analysis of Solid Rocket Engine[J]. Tactical Missile Technology, 2019(3): 88-93.

    [14] Chyuan S W. Dynamic Analysis of Solid Propellant Grains Subjected to Ignition Pressurization Loading[J]. Journal of Sound & Vibration, 2003, 268(3): 465-483.

    [15] Chyuan S W. Studies of Poisson’s Ratio Variation for Solid Propellant Grains under Ignition Pressure Loading[J]. International Journal of Pressure Vessels & Piping, 2003, 80(12): 871-877.

    [16] 朱浩, 朱亮, 呂先鋒, 等. 鋁合金薄壁擠壓管在準(zhǔn)靜態(tài)和動(dòng)態(tài)載荷下的有限元模擬[J]. 蘭州理工大學(xué)學(xué)報(bào), 2007(2): 34-38.

    Zhu Hao, Zhu Liang, Lü Xian-feng, et al. FEM Simulation of Thin-wall Aluminum Extruded Tubes under Quasi-static and Dynamic Loads[J]. Journal of Lanzhou University of Technology, 2007(2): 34-38.

    [17] 張學(xué)廣, 劉純國(guó), 鄭愿, 等. 基于延性損傷和剪切損傷的鋁合金成形極限預(yù)測(cè)[J]. 吉林大學(xué)學(xué)報(bào)(工學(xué)版), 2016, 46(5): 1558-1566.

    Zhang Xue-guang, Liu Chun-guo, Zheng Yuan, et al. Forming Limit Prediction of Aluminum Alloy Based on Ductile Damage and Shear Damage[J]. Journal of Jilin University(Engineering and Technology Edition), 2016, 46(5): 1558-1566.

    Analysis of Structural Integrity for the Hydro-reactive Metal Fuel Motor Grainl

    WANG Zhi-long1, HAN Xin-bo1, QIAO Hong1, YIN Shao-ping1, HAN Ai2, LEI Ming2

    (1.The 705 Research Institute, China State Shipbuilding Corporation Limited, Xi’an 710077, China,2. School of Astronautics, Northwestern Polytechnical University, Xi’an 710072, China)

    To study the structural integrity of a high-inclusion ratio grainl for a torpedo hydro-reactive metal fuel motor under local limit conditions, the mechanical properties of the grainl sample are tested, and the mechanical parameters are obtained. The numerical analysis model for the grainl sample is established using finite element software, and the simulation results are compared with the test results, which verified the rationality of the established constitutive model of grainl. The structural integrity of the hydro-reactive metal fuel motor grainl is analyzed using the constitutive model, and the force displacement curve, stress distribution, and displacement contour of the grainl are obtained. The simulation results indicate that the grainl is not damaged under typical working conditions, and the radial and axial deformation of the grainl can meet the requirements of the structure for a hydro-reactive motor.

    torpedo; hydro-reactive metal fuel motor; structural integrity for grainl; finite element analysis

    TJ630.32; TQ56

    A

    2096-3920(2021)06-0695-07

    10.11993/j.issn.2096-3920.2021.06.008

    王志龍, 韓新波, 喬宏, 等. 水反應(yīng)金屬燃料發(fā)動(dòng)機(jī)藥柱結(jié)構(gòu)完整性分析[J]. 水下無(wú)人系統(tǒng)學(xué)報(bào), 2021, 29(6): 695-701.

    2021-05-28;

    2021-07-20.

    國(guó)家自然科學(xué)基金重大項(xiàng)目(52090051).

    王志龍(1996-), 男, 在讀碩士, 主要研究方向?yàn)槟茉磩?dòng)力推進(jìn)技術(shù).

    (責(zé)任編輯: 楊力軍)

    猜你喜歡
    擋塊藥柱本構(gòu)
    高聚物黏結(jié)炸藥沖擊波感度試驗(yàn)方法
    四川化工(2022年6期)2023-01-15 10:54:54
    探討民用飛機(jī)艙門擋塊設(shè)計(jì)及公差
    一款彈簧式導(dǎo)向柱塞擋塊替換方案設(shè)計(jì)
    新型雙重?fù)鯄K構(gòu)造及地震反應(yīng)分析
    四川建筑(2020年1期)2020-07-21 07:26:06
    離心SC柱混凝土本構(gòu)模型比較研究
    鋸齒形結(jié)構(gòu)面剪切流變及非線性本構(gòu)模型分析
    橋梁抗震擋塊研究進(jìn)展
    四川水泥(2017年9期)2017-11-03 06:33:35
    一種新型超固結(jié)土三維本構(gòu)模型
    更 正
    含能材料(2017年6期)2017-03-07 06:32:07
    減面燃燒規(guī)律的藥柱選用準(zhǔn)則*
    久久草成人影院| 国产精品乱码一区二三区的特点| 欧洲精品卡2卡3卡4卡5卡区| 熟女人妻精品中文字幕| 最近最新中文字幕大全免费视频| 久久这里只有精品中国| 三级国产精品欧美在线观看 | 在线视频色国产色| 国产精品99久久99久久久不卡| 国产亚洲精品综合一区在线观看| 欧美绝顶高潮抽搐喷水| 亚洲国产精品合色在线| 国内精品久久久久精免费| 成人鲁丝片一二三区免费| 午夜激情欧美在线| 69av精品久久久久久| 国产精品自产拍在线观看55亚洲| 国产精品亚洲一级av第二区| 白带黄色成豆腐渣| 久久午夜亚洲精品久久| 国产欧美日韩一区二区三| 每晚都被弄得嗷嗷叫到高潮| 精品午夜福利视频在线观看一区| 国产精品一区二区三区四区久久| 最好的美女福利视频网| 女同久久另类99精品国产91| 亚洲av片天天在线观看| 国产高清有码在线观看视频| 国产精品久久久人人做人人爽| 曰老女人黄片| xxxwww97欧美| 99热只有精品国产| 精品国产亚洲在线| 国内精品久久久久精免费| 久久久久免费精品人妻一区二区| 88av欧美| 欧美成狂野欧美在线观看| 日韩大尺度精品在线看网址| 美女高潮喷水抽搐中文字幕| 黄色丝袜av网址大全| 麻豆国产97在线/欧美| 久久香蕉国产精品| 欧美激情在线99| 亚洲国产色片| 日本在线视频免费播放| 中文字幕av在线有码专区| 小说图片视频综合网站| 亚洲中文字幕一区二区三区有码在线看 | 午夜久久久久精精品| 在线十欧美十亚洲十日本专区| 999久久久国产精品视频| 婷婷六月久久综合丁香| 亚洲精品国产精品久久久不卡| 亚洲国产欧美一区二区综合| 亚洲欧美激情综合另类| 香蕉国产在线看| 夜夜夜夜夜久久久久| 中文在线观看免费www的网站| 国产成人啪精品午夜网站| 性欧美人与动物交配| 一边摸一边抽搐一进一小说| 亚洲精品粉嫩美女一区| 日本黄色视频三级网站网址| 夜夜躁狠狠躁天天躁| 给我免费播放毛片高清在线观看| 人人妻人人澡欧美一区二区| 国产精品香港三级国产av潘金莲| 一二三四社区在线视频社区8| 午夜亚洲福利在线播放| 精品欧美国产一区二区三| 亚洲第一欧美日韩一区二区三区| 99在线视频只有这里精品首页| 日本 欧美在线| 久久性视频一级片| 国产69精品久久久久777片 | 国产伦一二天堂av在线观看| 国产精品亚洲一级av第二区| 一区福利在线观看| а√天堂www在线а√下载| 在线观看免费午夜福利视频| 亚洲在线自拍视频| 国产精品国产高清国产av| 久久久久久九九精品二区国产| 久久99热这里只有精品18| 久久精品aⅴ一区二区三区四区| 亚洲av成人不卡在线观看播放网| 国产男靠女视频免费网站| 色综合站精品国产| 亚洲人成电影免费在线| 亚洲在线观看片| 亚洲成人免费电影在线观看| 国内揄拍国产精品人妻在线| 亚洲av熟女| 色视频www国产| 高潮久久久久久久久久久不卡| 国产 一区 欧美 日韩| 精品午夜福利视频在线观看一区| 午夜福利免费观看在线| 激情在线观看视频在线高清| 日韩欧美国产在线观看| 久久久色成人| 丰满的人妻完整版| 欧美日韩国产亚洲二区| 午夜福利在线观看吧| 一区二区三区激情视频| 国产美女午夜福利| 精品国产乱码久久久久久男人| 欧美性猛交╳xxx乱大交人| 免费看十八禁软件| 97超级碰碰碰精品色视频在线观看| 亚洲av成人精品一区久久| 欧美高清成人免费视频www| 最新在线观看一区二区三区| 国产日本99.免费观看| 日本a在线网址| 久久九九热精品免费| 久久精品91无色码中文字幕| 青草久久国产| 观看美女的网站| 国模一区二区三区四区视频 | 三级国产精品欧美在线观看 | 亚洲国产精品成人综合色| 国产精品99久久99久久久不卡| 国产激情偷乱视频一区二区| 国产欧美日韩一区二区精品| 亚洲av日韩精品久久久久久密| 一级a爱片免费观看的视频| 亚洲欧美精品综合一区二区三区| 人妻丰满熟妇av一区二区三区| 两性夫妻黄色片| 国产av不卡久久| 男女视频在线观看网站免费| 国产日本99.免费观看| 日日夜夜操网爽| 国产亚洲精品久久久久久毛片| 欧美在线黄色| 久久中文字幕人妻熟女| 丰满人妻一区二区三区视频av | 色播亚洲综合网| 亚洲成人久久性| 国产欧美日韩精品一区二区| 男人的好看免费观看在线视频| 久久久国产精品麻豆| 国产精品九九99| 午夜福利在线观看吧| 欧美成人免费av一区二区三区| 日本在线视频免费播放| 亚洲欧美日韩高清在线视频| 老司机深夜福利视频在线观看| 成人18禁在线播放| 黄片大片在线免费观看| 一个人免费在线观看的高清视频| 18禁黄网站禁片免费观看直播| 老熟妇仑乱视频hdxx| 亚洲va日本ⅴa欧美va伊人久久| 欧美一区二区精品小视频在线| 亚洲av美国av| 欧美xxxx黑人xx丫x性爽| 校园春色视频在线观看| 韩国av一区二区三区四区| 欧美一区二区精品小视频在线| 精品久久久久久,| АⅤ资源中文在线天堂| 亚洲美女视频黄频| 999久久久国产精品视频| 欧美绝顶高潮抽搐喷水| 男女那种视频在线观看| 99精品久久久久人妻精品| a级毛片a级免费在线| 一a级毛片在线观看| 两个人的视频大全免费| 两个人的视频大全免费| 又粗又爽又猛毛片免费看| 国产高清视频在线观看网站| 毛片女人毛片| 色综合欧美亚洲国产小说| 亚洲国产欧美网| 日日夜夜操网爽| 亚洲午夜理论影院| 小蜜桃在线观看免费完整版高清| 一本久久中文字幕| 国产高清有码在线观看视频| 淫秽高清视频在线观看| 国产高清激情床上av| 日本 av在线| 欧美激情在线99| 成熟少妇高潮喷水视频| 亚洲国产精品合色在线| 国产真人三级小视频在线观看| 国产黄色小视频在线观看| 国产97色在线日韩免费| www国产在线视频色| 欧美一级a爱片免费观看看| 国产人伦9x9x在线观看| 欧美一级毛片孕妇| 日韩有码中文字幕| 1024手机看黄色片| 少妇的逼水好多| 一夜夜www| 69av精品久久久久久| 99在线视频只有这里精品首页| 欧美3d第一页| 亚洲中文日韩欧美视频| 亚洲人成网站在线播放欧美日韩| 999久久久精品免费观看国产| 天堂av国产一区二区熟女人妻| 三级毛片av免费| 香蕉国产在线看| 十八禁人妻一区二区| 五月玫瑰六月丁香| 欧美在线一区亚洲| 亚洲精品乱码久久久v下载方式 | 99热这里只有精品一区 | 女人高潮潮喷娇喘18禁视频| 琪琪午夜伦伦电影理论片6080| 日本一二三区视频观看| 精品国产乱子伦一区二区三区| 欧美+亚洲+日韩+国产| 国产主播在线观看一区二区| 久久久久久久午夜电影| 亚洲国产精品成人综合色| 国产69精品久久久久777片 | 国产精品免费一区二区三区在线| 脱女人内裤的视频| 757午夜福利合集在线观看| 日本 av在线| tocl精华| 九色国产91popny在线| 精品一区二区三区视频在线观看免费| 亚洲真实伦在线观看| 国产主播在线观看一区二区| 宅男免费午夜| 亚洲性夜色夜夜综合| av片东京热男人的天堂| 成人高潮视频无遮挡免费网站| 久久精品国产清高在天天线| 成人性生交大片免费视频hd| 不卡av一区二区三区| 日韩中文字幕欧美一区二区| 久久久久亚洲av毛片大全| 草草在线视频免费看| 久久久久久久午夜电影| 一个人观看的视频www高清免费观看 | 国产精品久久久人人做人人爽| 婷婷精品国产亚洲av在线| 女生性感内裤真人,穿戴方法视频| 后天国语完整版免费观看| 桃红色精品国产亚洲av| 亚洲男人的天堂狠狠| 欧美最黄视频在线播放免费| 日韩欧美 国产精品| 中文字幕熟女人妻在线| 国产97色在线日韩免费| 国产av一区在线观看免费| 国产精品乱码一区二三区的特点| 美女午夜性视频免费| 色播亚洲综合网| 久久午夜综合久久蜜桃| 国产一区二区三区视频了| 欧美黑人欧美精品刺激| 午夜免费激情av| 免费在线观看日本一区| 亚洲av成人不卡在线观看播放网| 在线观看美女被高潮喷水网站 | 亚洲专区国产一区二区| 国产精品,欧美在线| 99在线视频只有这里精品首页| 人人妻人人看人人澡| 亚洲自拍偷在线| 亚洲国产欧美一区二区综合| 久久久久久久久免费视频了| 一区二区三区高清视频在线| 听说在线观看完整版免费高清| 色综合欧美亚洲国产小说| 国产亚洲av高清不卡| 亚洲国产欧洲综合997久久,| 热99re8久久精品国产| 日韩大尺度精品在线看网址| 99久久精品热视频| 美女扒开内裤让男人捅视频| 天天躁狠狠躁夜夜躁狠狠躁| 老司机午夜十八禁免费视频| 午夜福利高清视频| 真人一进一出gif抽搐免费| 中出人妻视频一区二区| 老汉色∧v一级毛片| 中国美女看黄片| 亚洲最大成人中文| 性色av乱码一区二区三区2| 一区二区三区国产精品乱码| 亚洲精品国产精品久久久不卡| 久久午夜综合久久蜜桃| 一区福利在线观看| 精品久久久久久成人av| 18禁黄网站禁片免费观看直播| 亚洲人成伊人成综合网2020| 亚洲 国产 在线| 午夜福利视频1000在线观看| 亚洲中文字幕日韩| 久久亚洲真实| 白带黄色成豆腐渣| 日本免费a在线| 波多野结衣巨乳人妻| 免费av毛片视频| 成人av在线播放网站| 国产精品电影一区二区三区| 好看av亚洲va欧美ⅴa在| 在线观看一区二区三区| 精品欧美国产一区二区三| 国产亚洲精品av在线| 午夜a级毛片| 别揉我奶头~嗯~啊~动态视频| 少妇熟女aⅴ在线视频| 精品国产美女av久久久久小说| 又大又爽又粗| 最近最新中文字幕大全电影3| 国产精品亚洲美女久久久| 日韩欧美精品v在线| 中国美女看黄片| 美女被艹到高潮喷水动态| 欧美日本亚洲视频在线播放| 欧美xxxx黑人xx丫x性爽| 一级毛片精品| 搡老岳熟女国产| 久久久久九九精品影院| 美女被艹到高潮喷水动态| 日韩欧美精品v在线| 国产成人影院久久av| 国产又黄又爽又无遮挡在线| 特级一级黄色大片| 国产精品免费一区二区三区在线| 亚洲精华国产精华精| 久久久久精品国产欧美久久久| 免费在线观看影片大全网站| 性色avwww在线观看| 国产精品国产高清国产av| 性色av乱码一区二区三区2| 91九色精品人成在线观看| 久久香蕉国产精品| 亚洲 欧美一区二区三区| 亚洲av免费在线观看| 高清在线国产一区| 午夜免费观看网址| 亚洲国产精品久久男人天堂| 91字幕亚洲| 一个人免费在线观看的高清视频| 国产av一区在线观看免费| 1000部很黄的大片| 国产精品久久视频播放| 国产久久久一区二区三区| 欧美日韩黄片免| 欧美日韩乱码在线| 亚洲在线观看片| 丁香欧美五月| 欧美日韩一级在线毛片| 日韩欧美三级三区| 亚洲人成电影免费在线| 国产精品亚洲美女久久久| 天堂√8在线中文| 熟女少妇亚洲综合色aaa.| 日本黄大片高清| 一个人免费在线观看电影 | 色在线成人网| 精品免费久久久久久久清纯| 久久人妻av系列| 久久欧美精品欧美久久欧美| 99久久无色码亚洲精品果冻| 亚洲 欧美一区二区三区| 国产不卡一卡二| 黄色片一级片一级黄色片| 麻豆成人av在线观看| 免费搜索国产男女视频| 三级男女做爰猛烈吃奶摸视频| 亚洲五月天丁香| 亚洲色图av天堂| 欧美zozozo另类| 曰老女人黄片| www.www免费av| 小说图片视频综合网站| 两性夫妻黄色片| 亚洲精品在线美女| 可以在线观看的亚洲视频| 亚洲 国产 在线| 久久这里只有精品19| 男人舔奶头视频| 久久久久性生活片| 天天躁狠狠躁夜夜躁狠狠躁| 校园春色视频在线观看| 黄色女人牲交| 亚洲人成网站高清观看| 国产一区二区在线av高清观看| 窝窝影院91人妻| 亚洲熟妇熟女久久| 看黄色毛片网站| 极品教师在线免费播放| 国内精品一区二区在线观看| 18美女黄网站色大片免费观看| 亚洲中文日韩欧美视频| 久久精品国产亚洲av香蕉五月| 久久性视频一级片| 国产成+人综合+亚洲专区| 五月玫瑰六月丁香| 老司机福利观看| 嫩草影视91久久| 禁无遮挡网站| 757午夜福利合集在线观看| 露出奶头的视频| 久久人妻av系列| 美女午夜性视频免费| 成人性生交大片免费视频hd| 亚洲国产精品合色在线| 九九热线精品视视频播放| 国产在线精品亚洲第一网站| www.自偷自拍.com| 给我免费播放毛片高清在线观看| 757午夜福利合集在线观看| 久久久久久大精品| 91老司机精品| 在线观看一区二区三区| av片东京热男人的天堂| 无遮挡黄片免费观看| 精品日产1卡2卡| 亚洲片人在线观看| 色视频www国产| 伊人久久大香线蕉亚洲五| 一本久久中文字幕| 搞女人的毛片| 久久久久久久午夜电影| 日本一本二区三区精品| 一区福利在线观看| 国产日本99.免费观看| 午夜视频精品福利| 男人的好看免费观看在线视频| 激情在线观看视频在线高清| 久久久精品欧美日韩精品| 久久热在线av| 97超视频在线观看视频| 黑人欧美特级aaaaaa片| 亚洲最大成人中文| 久久婷婷人人爽人人干人人爱| 又爽又黄无遮挡网站| 日韩欧美精品v在线| 成年免费大片在线观看| 精品电影一区二区在线| 禁无遮挡网站| 在线a可以看的网站| 久久中文字幕人妻熟女| 脱女人内裤的视频| 宅男免费午夜| 国产精品 欧美亚洲| 97人妻精品一区二区三区麻豆| 国产精品美女特级片免费视频播放器 | av中文乱码字幕在线| 婷婷丁香在线五月| 免费电影在线观看免费观看| 看片在线看免费视频| 俄罗斯特黄特色一大片| 亚洲欧美精品综合一区二区三区| 1024手机看黄色片| 久久久国产欧美日韩av| av黄色大香蕉| www.999成人在线观看| 真实男女啪啪啪动态图| 久久天堂一区二区三区四区| 午夜成年电影在线免费观看| 国产一区在线观看成人免费| 成人三级做爰电影| 男人舔奶头视频| 成人性生交大片免费视频hd| 国内精品美女久久久久久| 18禁裸乳无遮挡免费网站照片| 国产精品一区二区精品视频观看| 性色av乱码一区二区三区2| 午夜日韩欧美国产| 综合色av麻豆| 欧美高清成人免费视频www| 91在线精品国自产拍蜜月 | 在线观看一区二区三区| xxx96com| 淫秽高清视频在线观看| www日本在线高清视频| 成人无遮挡网站| 每晚都被弄得嗷嗷叫到高潮| 久久久国产成人精品二区| 色综合亚洲欧美另类图片| 日日夜夜操网爽| 一个人免费在线观看电影 | 一进一出抽搐动态| 美女大奶头视频| 中文在线观看免费www的网站| 国产一区二区激情短视频| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久九九精品影院| 日本五十路高清| 久久伊人香网站| 亚洲精品美女久久久久99蜜臀| 韩国av一区二区三区四区| 国产精品 国内视频| 日韩欧美 国产精品| 狂野欧美激情性xxxx| 黄片小视频在线播放| 一区二区三区国产精品乱码| 久久亚洲真实| 老熟妇仑乱视频hdxx| 成年女人毛片免费观看观看9| 国产欧美日韩一区二区三| 成人精品一区二区免费| 色播亚洲综合网| 悠悠久久av| 久久精品国产99精品国产亚洲性色| 午夜免费观看网址| 在线观看舔阴道视频| 嫩草影院入口| ponron亚洲| 一卡2卡三卡四卡精品乱码亚洲| 岛国视频午夜一区免费看| 久久久久免费精品人妻一区二区| 丝袜人妻中文字幕| 婷婷精品国产亚洲av在线| 丰满人妻熟妇乱又伦精品不卡| 成年女人永久免费观看视频| 亚洲av五月六月丁香网| 操出白浆在线播放| 久久精品aⅴ一区二区三区四区| 脱女人内裤的视频| 欧美zozozo另类| 亚洲成人久久性| 精品午夜福利视频在线观看一区| 一二三四社区在线视频社区8| 国产免费男女视频| 成年人黄色毛片网站| 在线视频色国产色| 精品国产乱子伦一区二区三区| 精品免费久久久久久久清纯| 老汉色∧v一级毛片| 久久久久久久午夜电影| 欧美日韩精品网址| 午夜激情欧美在线| 免费av毛片视频| 黑人操中国人逼视频| 十八禁人妻一区二区| 亚洲乱码一区二区免费版| xxx96com| 免费观看人在逋| 夜夜看夜夜爽夜夜摸| 国产日本99.免费观看| 久久天堂一区二区三区四区| 可以在线观看毛片的网站| 久久久久久久精品吃奶| 亚洲专区中文字幕在线| 亚洲男人的天堂狠狠| 国产亚洲欧美在线一区二区| 国产黄片美女视频| 99热6这里只有精品| 草草在线视频免费看| 人人妻人人看人人澡| 亚洲精品久久国产高清桃花| 日本与韩国留学比较| 亚洲人成网站高清观看| 嫩草影院精品99| 成人特级黄色片久久久久久久| 久久久国产成人精品二区| 国产高清视频在线观看网站| 欧美日韩瑟瑟在线播放| 操出白浆在线播放| 国产精品久久久久久精品电影| 91九色精品人成在线观看| 禁无遮挡网站| 欧美在线黄色| 精品久久久久久久人妻蜜臀av| 嫩草影院入口| 人人妻人人澡欧美一区二区| 久久久久久久久免费视频了| 欧美成人性av电影在线观看| 午夜视频精品福利| 丰满人妻熟妇乱又伦精品不卡| 听说在线观看完整版免费高清| 欧美三级亚洲精品| 老司机福利观看| 色综合站精品国产| 国产亚洲av高清不卡| 久久久国产成人精品二区| 天堂√8在线中文| 国语自产精品视频在线第100页| 99热6这里只有精品| 丰满人妻一区二区三区视频av | 亚洲午夜精品一区,二区,三区| 麻豆成人午夜福利视频| 黄色女人牲交| 欧美乱色亚洲激情| 在线观看日韩欧美| 一进一出好大好爽视频| 亚洲一区二区三区不卡视频| 午夜精品在线福利| 国产主播在线观看一区二区| 在线观看舔阴道视频| 成年人黄色毛片网站| av欧美777| 黑人操中国人逼视频| 男人和女人高潮做爰伦理| 久久精品夜夜夜夜夜久久蜜豆| 人人妻人人看人人澡| 青草久久国产| 中文字幕人成人乱码亚洲影| 在线国产一区二区在线| 久久九九热精品免费| av欧美777| 好看av亚洲va欧美ⅴa在| 久久久国产成人免费| 国产精品久久久久久精品电影| 老司机深夜福利视频在线观看| 国产精品美女特级片免费视频播放器 | 国产成人精品无人区| 成年女人看的毛片在线观看| 无人区码免费观看不卡|