• <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)則*
    亚洲国产精品专区欧美| 春色校园在线视频观看| 亚洲最大成人av| 99在线人妻在线中文字幕| 午夜精品一区二区三区免费看| av在线播放精品| 不卡视频在线观看欧美| 日本一二三区视频观看| 超碰av人人做人人爽久久| 免费看a级黄色片| 两个人的视频大全免费| 麻豆一二三区av精品| 国产亚洲av片在线观看秒播厂 | 超碰97精品在线观看| 日本wwww免费看| 国产精品一区二区性色av| 日本与韩国留学比较| 久久久久免费精品人妻一区二区| 中文精品一卡2卡3卡4更新| 三级男女做爰猛烈吃奶摸视频| 国产久久久一区二区三区| 成人鲁丝片一二三区免费| 在线免费观看的www视频| 男人舔女人下体高潮全视频| 欧美精品一区二区大全| 久久这里有精品视频免费| 性插视频无遮挡在线免费观看| 欧美区成人在线视频| 亚洲最大成人中文| 欧美一级a爱片免费观看看| 日韩一区二区三区影片| 国产精品爽爽va在线观看网站| 在线观看一区二区三区| 2021天堂中文幕一二区在线观| av.在线天堂| 岛国毛片在线播放| 精品久久久久久久人妻蜜臀av| 精品久久久久久久末码| 成人三级黄色视频| 亚洲国产高清在线一区二区三| 国产av一区在线观看免费| 日韩一区二区三区影片| 一级二级三级毛片免费看| 波野结衣二区三区在线| 亚洲国产成人一精品久久久| 国产精华一区二区三区| 日韩制服骚丝袜av| 午夜福利在线在线| 亚洲精品,欧美精品| 婷婷色麻豆天堂久久 | 少妇人妻精品综合一区二区| 久久久久久伊人网av| 村上凉子中文字幕在线| 亚洲精品乱码久久久久久按摩| 免费人成在线观看视频色| 美女cb高潮喷水在线观看| 国产高清不卡午夜福利| 国产在线男女| 51国产日韩欧美| 韩国高清视频一区二区三区| 亚洲av成人精品一区久久| 免费黄色在线免费观看| 亚洲精品一区蜜桃| 免费人成在线观看视频色| 国产精品.久久久| 国产黄片美女视频| 看非洲黑人一级黄片| 变态另类丝袜制服| 国产精品精品国产色婷婷| 69人妻影院| 97热精品久久久久久| 最近2019中文字幕mv第一页| 在线播放国产精品三级| 人体艺术视频欧美日本| 亚洲国产精品sss在线观看| 在线观看美女被高潮喷水网站| 99在线视频只有这里精品首页| 乱人视频在线观看| 亚洲欧美清纯卡通| 久热久热在线精品观看| 欧美变态另类bdsm刘玥| 免费看a级黄色片| 我的女老师完整版在线观看| 亚洲国产精品成人综合色| 国产精品美女特级片免费视频播放器| 波多野结衣巨乳人妻| 插逼视频在线观看| 国产一区二区在线观看日韩| 日本猛色少妇xxxxx猛交久久| 亚洲va在线va天堂va国产| 日本黄色片子视频| 免费一级毛片在线播放高清视频| 欧美不卡视频在线免费观看| 免费播放大片免费观看视频在线观看 | 麻豆成人av视频| 日本-黄色视频高清免费观看| 日韩精品青青久久久久久| 97人妻精品一区二区三区麻豆| 久久精品夜夜夜夜夜久久蜜豆| 丰满乱子伦码专区| 99久久精品国产国产毛片| 久久久国产成人精品二区| 日本一本二区三区精品| 精品少妇黑人巨大在线播放 | 69人妻影院| 中文字幕亚洲精品专区| 亚洲在久久综合| 免费看光身美女| 99久久精品热视频| 亚洲av一区综合| 国产老妇伦熟女老妇高清| av专区在线播放| 亚洲国产高清在线一区二区三| 精华霜和精华液先用哪个| 噜噜噜噜噜久久久久久91| 尤物成人国产欧美一区二区三区| 国产在线男女| 精品一区二区三区人妻视频| 国产精品一及| 性色avwww在线观看| 又爽又黄无遮挡网站| 可以在线观看毛片的网站| 你懂的网址亚洲精品在线观看 | 亚洲四区av| 成人亚洲精品av一区二区| 国产精品人妻久久久久久| av.在线天堂| 国产一区二区在线观看日韩| 搡老妇女老女人老熟妇| 国产精品一二三区在线看| 中文字幕久久专区| 欧美三级亚洲精品| 99久久精品热视频| 欧美成人一区二区免费高清观看| 午夜精品国产一区二区电影 | 亚洲熟妇中文字幕五十中出| 好男人视频免费观看在线| 啦啦啦啦在线视频资源| 色视频www国产| 国产精品女同一区二区软件| 人妻制服诱惑在线中文字幕| 国产精品蜜桃在线观看| 国产视频内射| 欧美成人免费av一区二区三区| 午夜激情福利司机影院| 日韩 亚洲 欧美在线| 欧美3d第一页| 99久国产av精品| 中文字幕免费在线视频6| 亚洲精品乱码久久久v下载方式| 99久国产av精品| 中文字幕人妻熟人妻熟丝袜美| 三级毛片av免费| 亚洲欧美清纯卡通| 国产精品久久久久久精品电影小说 | 乱系列少妇在线播放| 毛片一级片免费看久久久久| 欧美xxxx黑人xx丫x性爽| 国产精品爽爽va在线观看网站| 色哟哟·www| 五月玫瑰六月丁香| 亚洲精品乱久久久久久| 亚洲国产精品专区欧美| 一级黄色大片毛片| 日本黄大片高清| 日本-黄色视频高清免费观看| 免费av不卡在线播放| 精品不卡国产一区二区三区| 夜夜看夜夜爽夜夜摸| 国产一区有黄有色的免费视频 | 三级毛片av免费| 精品一区二区三区视频在线| 2021少妇久久久久久久久久久| 国内精品宾馆在线| 国产精品女同一区二区软件| 天美传媒精品一区二区| av在线播放精品| 校园人妻丝袜中文字幕| 久久久久久大精品| 亚洲欧美精品专区久久| 美女大奶头视频| 欧美日韩精品成人综合77777| 亚洲自拍偷在线| av.在线天堂| 看免费成人av毛片| 免费看光身美女| 午夜视频国产福利| 高清视频免费观看一区二区 | 精品国产三级普通话版| 日韩大片免费观看网站 | 白带黄色成豆腐渣| 亚洲欧美日韩东京热| 啦啦啦观看免费观看视频高清| 国产精品一区二区在线观看99 | 久久婷婷人人爽人人干人人爱| 在线免费观看的www视频| 精品人妻视频免费看| 少妇熟女aⅴ在线视频| 国产精品国产三级国产专区5o | 嫩草影院入口| 午夜福利在线观看吧| 男女国产视频网站| 亚洲婷婷狠狠爱综合网| 亚洲最大成人手机在线| 欧美色视频一区免费| 最近2019中文字幕mv第一页| 亚洲成人久久爱视频| 青春草国产在线视频| 国产老妇伦熟女老妇高清| 久久人人爽人人片av| 啦啦啦韩国在线观看视频| 日本黄色视频三级网站网址| 久久精品影院6| 午夜精品在线福利| 国产探花在线观看一区二区| 波多野结衣高清无吗| 成人特级av手机在线观看| 中文字幕精品亚洲无线码一区| 欧美极品一区二区三区四区| 伊人久久精品亚洲午夜| 亚洲国产精品sss在线观看| 精品免费久久久久久久清纯| 亚洲久久久久久中文字幕| 久久久久久久久中文| 国产高潮美女av| 国产国拍精品亚洲av在线观看| 最近中文字幕高清免费大全6| 国产精品一二三区在线看| 日日摸夜夜添夜夜添av毛片| 日本免费一区二区三区高清不卡| 亚洲成人精品中文字幕电影| 日韩av在线大香蕉| 欧美一级a爱片免费观看看| 欧美区成人在线视频| 国产精品一及| 能在线免费看毛片的网站| 国产又黄又爽又无遮挡在线| 91久久精品电影网| 国产av不卡久久| 少妇的逼好多水| 欧美成人a在线观看| av在线老鸭窝| 神马国产精品三级电影在线观看| 精品久久久久久久末码| 成人综合一区亚洲| 亚洲精品亚洲一区二区| 日韩欧美精品免费久久| 成人综合一区亚洲| 亚洲欧美日韩无卡精品| 青青草视频在线视频观看| 黄色欧美视频在线观看| av播播在线观看一区| 国产视频首页在线观看| 亚洲丝袜综合中文字幕| 黄片无遮挡物在线观看| 97热精品久久久久久| 高清视频免费观看一区二区 | 综合色丁香网| 国产精品久久电影中文字幕| 精品国产一区二区三区久久久樱花 | 日韩av在线大香蕉| 午夜激情福利司机影院| 国产精品乱码一区二三区的特点| 亚洲美女视频黄频| 久久韩国三级中文字幕| 欧美性感艳星| 日韩欧美精品免费久久| 一级毛片aaaaaa免费看小| 又粗又爽又猛毛片免费看| www日本黄色视频网| 亚洲欧美精品自产自拍| 国产免费一级a男人的天堂| 亚洲怡红院男人天堂| 国产精品三级大全| 少妇裸体淫交视频免费看高清| 色播亚洲综合网| av女优亚洲男人天堂| 国产精品一二三区在线看| 午夜精品国产一区二区电影 | 嫩草影院精品99| 国产黄a三级三级三级人| 国产精品一区二区三区四区久久| 嫩草影院精品99| 天天躁日日操中文字幕| 午夜爱爱视频在线播放| 成人二区视频| 麻豆乱淫一区二区| 小蜜桃在线观看免费完整版高清| 日韩一本色道免费dvd| 日本一二三区视频观看| 午夜视频国产福利| 美女国产视频在线观看| 欧美xxxx性猛交bbbb| 国产激情偷乱视频一区二区| 熟女电影av网| 婷婷色麻豆天堂久久 | 中文字幕免费在线视频6| 久久精品国产亚洲网站| 久久精品国产鲁丝片午夜精品| 伦精品一区二区三区| 国产成人精品久久久久久| www日本黄色视频网| 国产亚洲最大av| 久久久久九九精品影院| 国产69精品久久久久777片| 精品久久久噜噜| 免费av不卡在线播放| av在线观看视频网站免费| 国产成人精品婷婷| 久久午夜福利片| 亚洲最大成人手机在线| 国产一级毛片在线| 小蜜桃在线观看免费完整版高清| 身体一侧抽搐| 日韩制服骚丝袜av| 91在线精品国自产拍蜜月| 深爱激情五月婷婷| 亚洲精品影视一区二区三区av| 亚洲av不卡在线观看| 亚洲成av人片在线播放无| 1000部很黄的大片| 最近最新中文字幕免费大全7| 一二三四中文在线观看免费高清| 精品人妻熟女av久视频| 久久精品国产亚洲av天美| 国产午夜精品论理片| 亚洲av日韩在线播放| 久久久久久伊人网av| 免费观看人在逋| 三级国产精品片| 又爽又黄a免费视频| 看免费成人av毛片| 少妇的逼水好多| 一个人看视频在线观看www免费| 欧美人与善性xxx| 亚洲电影在线观看av| 在线观看66精品国产| 高清午夜精品一区二区三区| 日本午夜av视频| 日本猛色少妇xxxxx猛交久久| 亚洲五月天丁香| 久久国内精品自在自线图片| 一卡2卡三卡四卡精品乱码亚洲| 久久这里有精品视频免费| 精品一区二区三区人妻视频| 国产一区二区三区av在线| 在线观看美女被高潮喷水网站| 最近2019中文字幕mv第一页| 国产成人a区在线观看| 欧美三级亚洲精品| 免费av毛片视频| 乱系列少妇在线播放| 国产大屁股一区二区在线视频| 国产精品一区二区三区四区免费观看| av在线观看视频网站免费| 国产精品爽爽va在线观看网站| 最近2019中文字幕mv第一页| 亚洲图色成人| 国产淫语在线视频| 久久人妻av系列| 亚洲精华国产精华液的使用体验| www.av在线官网国产| 男人狂女人下面高潮的视频| 大又大粗又爽又黄少妇毛片口| 你懂的网址亚洲精品在线观看 | 天天躁夜夜躁狠狠久久av| 亚洲熟妇中文字幕五十中出| 日韩av在线大香蕉| 成年av动漫网址| 天堂影院成人在线观看| 国产精品国产高清国产av| 久久99蜜桃精品久久| 男女国产视频网站| 亚洲欧美日韩高清专用| 久久草成人影院| 午夜免费男女啪啪视频观看| 国产三级中文精品| 成人漫画全彩无遮挡| 啦啦啦啦在线视频资源| 成年版毛片免费区| 亚洲美女视频黄频| 插阴视频在线观看视频| 高清在线视频一区二区三区 | 亚洲欧美日韩东京热| 亚洲一级一片aⅴ在线观看| 韩国高清视频一区二区三区| 国产免费一级a男人的天堂| 久久精品久久精品一区二区三区| 成人欧美大片| eeuss影院久久| 国产黄色小视频在线观看| 免费看av在线观看网站| 男人舔奶头视频| 国产麻豆成人av免费视频| 男人舔奶头视频| 建设人人有责人人尽责人人享有的 | 亚洲综合色惰| 欧美一级a爱片免费观看看| 国产国拍精品亚洲av在线观看| 成人一区二区视频在线观看| or卡值多少钱| 国产成人freesex在线| 国产黄片视频在线免费观看| 搞女人的毛片| 天堂√8在线中文| 欧美成人精品欧美一级黄| 成年女人看的毛片在线观看| 久久综合国产亚洲精品| 欧美潮喷喷水| 精品国产三级普通话版| 国产精品一二三区在线看| 欧美丝袜亚洲另类| 夜夜爽夜夜爽视频| 一级av片app| 两个人视频免费观看高清| 九色成人免费人妻av| 99国产精品一区二区蜜桃av| 2021少妇久久久久久久久久久| 久久久a久久爽久久v久久| 国产精品99久久久久久久久| 亚洲成人av在线免费| 亚洲怡红院男人天堂| 毛片一级片免费看久久久久| 69av精品久久久久久| 色视频www国产| 高清视频免费观看一区二区 | 在线天堂最新版资源| 亚洲18禁久久av| 99国产精品一区二区蜜桃av| 99久久成人亚洲精品观看| videos熟女内射| 婷婷色麻豆天堂久久 | 亚洲国产成人一精品久久久| 只有这里有精品99| 偷拍熟女少妇极品色| 国产在视频线在精品| 尾随美女入室| 少妇丰满av| 精品无人区乱码1区二区| 中文字幕制服av| 婷婷色麻豆天堂久久 | 亚洲最大成人手机在线| 欧美激情久久久久久爽电影| 床上黄色一级片| 毛片女人毛片| 国产精品人妻久久久久久| 国产亚洲91精品色在线| 久久久a久久爽久久v久久| 久久久久久久久久黄片| 99久久成人亚洲精品观看| 国产伦理片在线播放av一区| 亚洲欧洲国产日韩| 97在线视频观看| 欧美性猛交黑人性爽| 久久99热6这里只有精品| 午夜精品一区二区三区免费看| 啦啦啦韩国在线观看视频| 亚洲av二区三区四区| 亚洲伊人久久精品综合 | 精品熟女少妇av免费看| 国产一区有黄有色的免费视频 | 草草在线视频免费看| 国产亚洲精品久久久com| 99热精品在线国产| 久久99蜜桃精品久久| 少妇的逼好多水| 国产亚洲一区二区精品| 人人妻人人看人人澡| 2021天堂中文幕一二区在线观| 亚洲精品乱久久久久久| 成人性生交大片免费视频hd| 熟妇人妻久久中文字幕3abv| 你懂的网址亚洲精品在线观看 | 1000部很黄的大片| 日韩欧美精品免费久久| 亚洲国产精品合色在线| 成人漫画全彩无遮挡| 小蜜桃在线观看免费完整版高清| 夫妻性生交免费视频一级片| 亚洲综合色惰| 亚洲欧美中文字幕日韩二区| 啦啦啦观看免费观看视频高清| 免费看a级黄色片| 日本爱情动作片www.在线观看| 精品久久久久久久久久久久久| 欧美三级亚洲精品| 久久久久久国产a免费观看| 禁无遮挡网站| 特级一级黄色大片| 国产亚洲5aaaaa淫片| 成人特级av手机在线观看| 久热久热在线精品观看| 网址你懂的国产日韩在线| 九草在线视频观看| 少妇被粗大猛烈的视频| 国产黄片视频在线免费观看| 日韩欧美国产在线观看| 中文天堂在线官网| 欧美一区二区国产精品久久精品| 成人二区视频| 成年女人看的毛片在线观看| 亚洲国产精品成人综合色| 伊人久久精品亚洲午夜| 亚洲美女视频黄频| 老师上课跳d突然被开到最大视频| 只有这里有精品99| 国产精品国产三级国产专区5o | 性色avwww在线观看| 成人亚洲欧美一区二区av| 视频中文字幕在线观看| 99久久精品一区二区三区| 国产精品一区二区三区四区久久| 精品久久久噜噜| 变态另类丝袜制服| 亚洲人成网站高清观看| 免费看美女性在线毛片视频| 久久久久久久久久成人| 深爱激情五月婷婷| 久久久久免费精品人妻一区二区| 欧美一区二区国产精品久久精品| 久久精品国产鲁丝片午夜精品| 少妇的逼好多水| 午夜视频国产福利| 成人特级av手机在线观看| 波多野结衣巨乳人妻| 国产男人的电影天堂91| 久久久a久久爽久久v久久| 国产视频首页在线观看| 亚洲国产日韩欧美精品在线观看| 久久久精品欧美日韩精品| 在线播放国产精品三级| 深夜a级毛片| 国产视频内射| 欧美日本视频| 免费看a级黄色片| 午夜日本视频在线| 免费大片18禁| 天堂影院成人在线观看| av.在线天堂| 久久久久久国产a免费观看| 国产精品蜜桃在线观看| av天堂中文字幕网| av在线天堂中文字幕| 色5月婷婷丁香| 天堂网av新在线| 免费无遮挡裸体视频| 精品一区二区三区视频在线| 精品国产露脸久久av麻豆 | 亚洲激情五月婷婷啪啪| 1024手机看黄色片| 日本免费一区二区三区高清不卡| 三级毛片av免费| 亚洲成人精品中文字幕电影| 国产色爽女视频免费观看| 综合色丁香网| 亚洲av日韩在线播放| 女的被弄到高潮叫床怎么办| 国产精品一区二区三区四区久久| 黄色配什么色好看| 狂野欧美白嫩少妇大欣赏| 亚洲国产精品专区欧美| 欧美高清成人免费视频www| 国语对白做爰xxxⅹ性视频网站| 自拍偷自拍亚洲精品老妇| 久久久久久久亚洲中文字幕| 国产高清有码在线观看视频| 日本免费一区二区三区高清不卡| av线在线观看网站| 国产视频首页在线观看| 国产精品久久久久久久电影| 亚洲精品456在线播放app| 高清av免费在线| 一个人看的www免费观看视频| 日韩 亚洲 欧美在线| 国产成人91sexporn| 国产伦一二天堂av在线观看| 国产亚洲av嫩草精品影院| 97人妻精品一区二区三区麻豆| 免费人成在线观看视频色| 男女视频在线观看网站免费| 日本五十路高清| 亚洲在线自拍视频| 国产免费视频播放在线视频 | 一级毛片电影观看 | 91久久精品国产一区二区成人| 日韩一本色道免费dvd| 日韩欧美国产在线观看| 日韩成人av中文字幕在线观看| 精品国产一区二区三区久久久樱花 | 高清av免费在线| 亚洲av二区三区四区| 99热6这里只有精品| 97超视频在线观看视频| 不卡视频在线观看欧美| 国产美女午夜福利| 婷婷色综合大香蕉| 国产在线一区二区三区精 | 亚洲欧美日韩高清专用| 亚洲真实伦在线观看| 成年av动漫网址| 黑人高潮一二区| 免费av观看视频| 国产黄色小视频在线观看| 人体艺术视频欧美日本| 精品人妻熟女av久视频| 国产一区二区三区av在线| 亚洲不卡免费看| 国产高清视频在线观看网站| 国产亚洲精品av在线| 日韩av在线大香蕉| 国产熟女欧美一区二区| 成人性生交大片免费视频hd| 久久久午夜欧美精品| 99热这里只有精品一区| 中文精品一卡2卡3卡4更新|