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

    典型裝藥水下爆炸的殉爆規(guī)律研究

    2014-02-27 06:30:08魯忠寶胡宏偉
    關(guān)鍵詞:模型

    魯忠寶, 胡宏偉, 劉 銳, 楊 帆

    ?

    典型裝藥水下爆炸的殉爆規(guī)律研究

    魯忠寶1, 胡宏偉2, 劉 銳1, 楊 帆1

    (1. 中國(guó)船舶重工集團(tuán)公司 第705研究所, 陜西 西安, 710075; 2. 西安近代化學(xué)研究所, 陜西 西安, 710065)

    針對(duì)典型的殼裝炸藥殉爆更接近炸藥實(shí)際使用狀態(tài), 采用ANSYS/LS_DYNA軟件建立了典型裝藥水下殉爆的有限元仿真模型, 通過(guò)計(jì)算得到了殉爆距離與安全距離, 基于此加工了試驗(yàn)樣彈, 并進(jìn)行了相應(yīng)的水下殉爆試驗(yàn)。試驗(yàn)結(jié)果與數(shù)值仿真結(jié)果較為吻合, 表明了本文的仿真模型能夠有效描述帶殼裝藥的水下殉爆情況, 水下試驗(yàn)中由沖擊波壓力和氣泡周期來(lái)判斷被發(fā)裝藥是否殉爆是可行的。最后, 在水下殉爆理論分析基礎(chǔ)上, 根據(jù)試驗(yàn)結(jié)果預(yù)計(jì)了大藥量水中兵器戰(zhàn)斗部的安全性。本文的研究可為水中兵器戰(zhàn)場(chǎng)環(huán)境的安全性設(shè)計(jì)提供依據(jù), 并可為水中兵器戰(zhàn)斗部殉爆毀傷能力評(píng)估提供參考。

    水中兵器; 戰(zhàn)斗部; 水下爆炸; 殉爆; 有限元仿真模型; 安全性

    0 引言

    炸藥作為一種含能物質(zhì), 在受到足夠強(qiáng)的外界刺激如力、熱、光和電等作用時(shí)會(huì)發(fā)生爆炸, 并可能造成臨近炸藥的殉爆[1]。炸藥殉爆主要有裸裝炸藥殉爆和殼裝炸藥殉爆2類, 其中殼裝炸藥殉爆更接近炸藥實(shí)際使用狀態(tài)。根據(jù)所處的環(huán)境, 主要有空氣中殉爆與水下殉爆2類。目前, 國(guó)內(nèi)外對(duì)炸藥空氣中的殉爆研究較多, 如Howe等人采用歐拉程序2DE對(duì)殼裝炸藥殉爆進(jìn)行了2D數(shù)值仿真[2]; Lu 等人對(duì)裸裝PBXN2109炸藥殉爆進(jìn)行了數(shù)值仿真[3]; Fisher等人對(duì)殼裝PBXN29 炸藥殉爆進(jìn)行了數(shù)值仿真[4]; 陳朗等人計(jì)算了不同距離下裸裝炸藥的殉爆距離[5]; 文獻(xiàn)[6]的表1中, 給出了當(dāng)主發(fā)藥和被發(fā)藥均為柱狀TNT散裝單質(zhì)炸藥時(shí)的殉爆試驗(yàn)數(shù)據(jù); 王晨等人進(jìn)行殼裝固黑鋁(GHL)炸藥殉爆試驗(yàn)[7]; 周保順等對(duì)沖擊波作用下炸藥殉爆過(guò)程進(jìn)行了數(shù)值仿真并進(jìn)行了試驗(yàn)驗(yàn)證[8]。不過(guò), 對(duì)炸藥水中殉爆研究尚未見(jiàn)詳細(xì)報(bào)道。

    對(duì)于水中兵器, 戰(zhàn)斗部水下爆炸可能激起一定距離上的另一水中兵器戰(zhàn)斗部發(fā)生爆炸, 即產(chǎn)生殉爆, 殉爆也是摧毀敵方水中兵器的作用方式之一。研究典型裝藥水下爆炸的殉爆機(jī)理, 分析水中兵器戰(zhàn)斗部對(duì)來(lái)襲水中兵器戰(zhàn)斗部的安全距離與殉爆距離, 可為水中兵器戰(zhàn)場(chǎng)環(huán)境的安全性設(shè)計(jì)提供依據(jù), 同時(shí)也可為打擊水中兵器用的水中兵器戰(zhàn)斗部的毀傷威力評(píng)定提供一種新方法。

    1 水中殉爆的理論分析

    1.1 影響水中殉爆的因素

    水中殉爆規(guī)律很復(fù)雜, 引起殉爆的原因可能有主發(fā)裝藥爆炸產(chǎn)生的沖擊波、氣泡、二次壓力波及爆轟產(chǎn)物直接作用等。影響水中殉爆的因素很多, 主要有以下幾個(gè)方面。

    1.1.1主發(fā)裝藥的影響

    不同種類主發(fā)炸藥其爆炸能量釋放與輸出特性存在差異, 其殉爆距離也不同。相同種類的炸藥, 主發(fā)裝藥的藥量越大, 能量越高, 被發(fā)裝藥越容易被起爆。主發(fā)裝藥帶外殼時(shí), 爆轟產(chǎn)物受到約束, 增大了主發(fā)裝藥的爆轟性能, 使殉爆距離增大。外殼強(qiáng)度越大, 這種效應(yīng)越顯著。

    1.1.2被發(fā)裝藥的影響

    被發(fā)裝藥對(duì)外界刺激越敏感, 殉爆距離越大。裝藥密度較低時(shí), 也容易殉爆。另外, 外殼對(duì)爆轟產(chǎn)物、沖擊波和破片有衰減、阻擋作用, 能導(dǎo)致殉爆距離減小。

    1.1.3裝藥的布局

    主發(fā)裝藥與被發(fā)裝藥按同軸線的擺放方式比按軸線垂直擺放方式更容易殉爆, 起爆方式對(duì)殉爆距離也有影響。

    當(dāng)主發(fā)裝藥和被發(fā)裝藥的種類確定時(shí), 影響殉爆的主要因素為殼體厚度和材料, 殼體厚度和材料都會(huì)對(duì)沖擊波大小產(chǎn)生影響。有研究表明, 帶殼裝藥水下爆炸的沖擊波峰值壓力隨殼厚與裝藥半徑比的增加先增大后減小, 當(dāng)殼厚裝藥半徑比為0.07~0.08時(shí), 峰值壓力隨殼厚變化的趨勢(shì)發(fā)生改變, 對(duì)主發(fā)裝藥適當(dāng)?shù)臍んw厚度可以增大殉爆距離。但對(duì)被發(fā)裝藥, 可能會(huì)導(dǎo)致殉爆距離減小。當(dāng)主、被發(fā)裝藥都帶殼體時(shí), 難以判斷出殼體對(duì)殉爆距離的影響。

    1.2 水中爆炸的起爆判據(jù)和殉爆閾值

    1.2.1水中爆炸的起爆判據(jù)

    起爆判據(jù)是判斷主發(fā)裝藥爆炸后, 被發(fā)裝藥是否爆炸的依據(jù), 空氣中一般用爆炸產(chǎn)生的特征參量(沖擊波超壓)或試驗(yàn)現(xiàn)象(殘藥或見(jiàn)證板的破壞情況)來(lái)判別。對(duì)于水中爆炸, 被發(fā)裝藥的殘藥或殘留外殼、見(jiàn)證板難于回收, 可采用水中沖擊波壓力和氣泡周期來(lái)判斷被發(fā)裝藥是否發(fā)生爆轟。

    1.2.2水中爆炸的起爆閾值

    沖擊起爆是炸藥對(duì)壓力脈沖的一種動(dòng)態(tài)響應(yīng)形式。對(duì)于均質(zhì)炸藥, 早期認(rèn)為只要沖擊波壓力超過(guò)臨界起爆壓力, 炸藥就會(huì)被引爆。對(duì)于非均質(zhì)炸藥, 其起爆閾值既與沖擊波壓力有關(guān), 又與其持續(xù)時(shí)間有關(guān)。因此, 水中爆炸殉爆的臨界起爆判據(jù)需要綜合考慮臨界起爆壓力或臨界起爆能量, 殉爆臨界起爆閾值應(yīng)符合:>,≥(為沖擊波能,為臨界起爆能)。

    1.3 帶殼裝藥水中爆炸的殉爆相似律

    函數(shù)的具體形式可通過(guò)試驗(yàn)確定, 當(dāng)以沖擊波超壓作為判別殉爆的主要參數(shù)時(shí), 水中殉爆距離與主發(fā)藥量、殼體厚度和裝藥半徑(包括主發(fā)和被發(fā)裝藥)的比值有關(guān)。

    2 數(shù)值仿真

    有限元分析軟件ANSYS/LS_DYNA強(qiáng)大的自適應(yīng)線譜增強(qiáng)算法(adaptive line enhancement, ALE)和Euler算法及炸藥的材料與狀態(tài)方程廣泛應(yīng)用于各種水下爆炸分析中, 相關(guān)資料表明[9], ANSYS/LS-DYNA在水中爆炸分析中有很多成功的案例。對(duì)炸藥水下殉爆試驗(yàn)的數(shù)值仿真可以描述炸藥在殉爆中的反應(yīng)規(guī)律, 同時(shí)能夠獲得殉爆過(guò)程中的細(xì)節(jié)規(guī)律, 在很大程度上可以減少試驗(yàn)數(shù)量。

    2.1 計(jì)算模型

    2.1.1 物理模型及算法

    炸藥裝藥水下爆炸殉爆模型見(jiàn)圖1,為2個(gè)裝藥的中心距,為2個(gè)裝藥的邊距。主發(fā)裝藥與被發(fā)裝藥外形尺寸均為f90 mm×90 mm的圓柱體, 裝藥類型均為GUHL-1, 主發(fā)裝藥包裹有 1.5 mm的鋁殼, 被發(fā)裝藥包裹有1 mm的鋁殼。

    圖1 裝藥水下殉爆示意圖

    計(jì)算模型由主發(fā)裝藥、被發(fā)裝藥、主發(fā)裝藥外殼、被發(fā)裝藥外殼和水介質(zhì)5部分組成, 其中主發(fā)裝藥、水介質(zhì)采用歐拉網(wǎng)格建模, 單元使用多物質(zhì)ALE算法, 被發(fā)裝藥、主發(fā)裝藥外殼、被發(fā)裝藥外殼采用LAGRANGE算法, 在被發(fā)裝藥、主發(fā)裝藥外殼、被發(fā)裝藥外殼與水域間采用關(guān)鍵字*CONSTRAINED_LAGRANGE_IN_SOLID實(shí)現(xiàn)耦合算法[10]。主發(fā)裝藥采用圓柱體中心起爆, 根據(jù)模型對(duì)稱性特點(diǎn), 主發(fā)裝藥及其外殼采用1/2有限元建模。選擇SOLID164單元和相應(yīng)的材料屬性對(duì)幾何模型劃分網(wǎng)格, 創(chuàng)建PART, 施加對(duì)稱約束和透射邊界條件, 調(diào)整主發(fā)裝藥與被發(fā)裝藥間的距離, 設(shè)置求解時(shí)間等操作。形成有限元模型如圖2所示。

    圖2 裝藥水下殉爆有限元模型

    2.1.2材料與狀態(tài)方程

    主發(fā)裝藥及被發(fā)裝藥的鋁殼選用各項(xiàng)同性隨動(dòng)強(qiáng)化的彈塑性金屬材料模型MAT_PLASTIC_ KINEMATIC, 密度=2.7×103kg/m3, 彈性模量=70 GPa,泊松比=0.35, 屈服應(yīng)力=0.3 GPa; 水選用空白材料模型NULL, 狀態(tài)方程為EOS_GRUNEISEN, 密度為水=1.025×103kg/m3; 主發(fā)裝藥選用高能炸藥燃燒與增長(zhǎng)模型MAT_ HIGH_EXPLOSIVE_BURN, EOS_JWL狀態(tài)方程, 其密度藥=1.8×103kg/m3, 爆速=7020 m/s, 爆壓P=21.5 GPa, 炸藥及水的具體參數(shù)可參考文獻(xiàn)[11]。被發(fā)裝藥采用流體彈塑性材料模型*MAT _ELASTIC_PLASTIC_HYDRO[12], 采用點(diǎn)火增長(zhǎng)反應(yīng)速率方程*EOS_IGNITION_AND_ GROWTH_ OF_REACTION_IN_HE[12], 方程如下

    GUHL-1炸藥點(diǎn)火增長(zhǎng)模型參數(shù)選取如表1所示。

    2.2 仿真計(jì)算結(jié)果

    選取幾種不同作用距離進(jìn)行仿真模型的計(jì)算。當(dāng)=3 cm時(shí), 被發(fā)裝藥未被殉爆, 沖擊波傳播、被發(fā)裝藥變形、被發(fā)裝藥內(nèi)部#120154單元的壓力時(shí)程曲線(峰值壓力P= 2.31714GPa)見(jiàn)圖3。當(dāng)=1 cm時(shí), 被發(fā)裝藥被殉爆, 沖擊波傳播、被發(fā)裝藥殉爆、被發(fā)裝藥內(nèi)部#120154單元的壓力時(shí)程曲線(峰值壓力P=23.224 5 GPa)見(jiàn)圖4。

    表1 GUHL-1炸藥點(diǎn)火增長(zhǎng)模型反應(yīng)速率方程參數(shù)

    圖3 作用距離R=3 cm時(shí), 沖擊波傳播、被發(fā)裝藥變形以及內(nèi)部單元壓力時(shí)程曲線圖

    圖4 R=1 cm時(shí), 沖擊波傳播、被發(fā)裝藥殉爆以及內(nèi)部單元壓力時(shí)程曲線圖

    3 殉爆試驗(yàn)

    按照2.1中的仿真研究模型, 加工了試驗(yàn)樣彈, 尺寸都為f90mm×90mm, 裝藥品種都為GUHL-1, 密度大于1.8×103kg/m3, 主發(fā)裝藥鋁殼厚1.5 mm, 傳爆藥柱為20 g JH-14, 采用8#銅殼電雷管端面中心起爆。被發(fā)裝藥鋁殼厚1 mm。

    試驗(yàn)水池為圓柱形結(jié)構(gòu), 試驗(yàn)時(shí)將主發(fā)裝藥、被發(fā)裝藥懸吊在水池中心, 入水深度5 m, 并處于同一水平面上, 主發(fā)裝藥和被發(fā)裝藥的距離通過(guò)鐵絲編制的架子來(lái)控制。同時(shí)在距爆心2.5 m處布放3個(gè)水下傳感器, 傳感器與兩裝藥處于同一水平面。測(cè)定炸藥爆炸時(shí)水中壓力隨時(shí)間而變化的-曲線。

    測(cè)試系統(tǒng)由高速采集儀(見(jiàn)圖5)、適配器、水下沖擊波壓力傳感器(見(jiàn)圖6)、低噪聲電纜和專用軟件組成。傳感器采用美國(guó)PCB公司生產(chǎn)的PCB138A型大量程電氣石壓力傳感器, 測(cè)量系統(tǒng)有專用軟件, 能自動(dòng)找出氣泡脈動(dòng)周期和沖擊波峰值壓力等參數(shù)。

    首先測(cè)量了單個(gè)1 kg的GUHL-1炸藥(主發(fā)裝藥, 見(jiàn)圖7)水中爆炸的沖擊波壓力和氣泡周期, 測(cè)距為2.5 m, 3點(diǎn)的峰值壓力平均值為16.27 MPa, 氣泡周期平均值為322.9 ms。GUHL-1和GUHL-1的殉爆試驗(yàn)裝配見(jiàn)圖8。殉爆的試驗(yàn)結(jié)果見(jiàn)表2。

    圖5 水下數(shù)據(jù)采集系統(tǒng)

    圖6 傳感器

    圖7 單個(gè)1 kg的GUHL-1炸藥

    圖8 GUHL-1和GUHL-1的殉爆試驗(yàn)裝配圖

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

    針對(duì)相同的主發(fā)裝藥、被發(fā)裝藥的水下殉爆研究模型, 采用數(shù)值仿真以及試驗(yàn)測(cè)試的方式展開(kāi)研究。

    由圖4可看出, 在=1 cm時(shí), 被發(fā)裝藥發(fā)生膨脹, 其內(nèi)部一單元的峰值壓力達(dá)23.224 5 GPa,大于該炸藥的爆壓P, 被發(fā)裝藥被殉爆; 由圖3可看出, 在= 3 cm時(shí), 被發(fā)裝藥發(fā)生壓縮變形, 內(nèi)部同一單元的峰值壓力僅為2.31714 GPa, 遠(yuǎn)遠(yuǎn)小于該炸藥的爆壓P, 被發(fā)裝藥未被殉爆。同樣可得在= 6 cm和= 8 cm工況下, 被發(fā)裝藥均未被殉爆。采用數(shù)值仿真方法, 可以明顯地觀測(cè)到被發(fā)裝藥是否被殉爆。數(shù)值仿真結(jié)果與試驗(yàn)結(jié)果的吻合表明, 采用數(shù)值仿真方法能較好地模擬水下殉爆, 本文2.1節(jié)中的材料與狀態(tài)方程參數(shù)是合理可行的。

    表2 GUHL-1和GUHL-1的殉爆試驗(yàn)結(jié)果

    當(dāng)殉爆距離較近時(shí), 難以準(zhǔn)確測(cè)量沖擊波峰值壓力, 通過(guò)數(shù)值仿真計(jì)算可知, 被發(fā)裝藥GUHL-1的殉爆沖擊波壓力閾值約為7.49 GPa。

    由表2可知, 2個(gè)裝藥表面距離為3 cm, 6 cm和8 cm時(shí), 測(cè)點(diǎn)處的沖擊波峰值壓力和氣泡周期與單個(gè)1 kg的GUHL-1裝藥一致, 被發(fā)裝藥未殉爆。2個(gè)裝藥表面距離為1 cm時(shí), 測(cè)點(diǎn)處的沖擊波峰值壓力和氣泡周期都大于單個(gè)1 kg的GUHL-1裝藥的沖擊波峰值壓力和氣泡周期, 被發(fā)裝藥殉爆, 與數(shù)值仿真結(jié)果吻合。研究表明, 在水下殉爆試驗(yàn)中, 由沖擊波壓力和氣泡周期來(lái)判斷被發(fā)裝藥是否殉爆是可行的。

    當(dāng)主發(fā)裝藥和被發(fā)裝藥都為帶殼的1 kg的GUHL-1時(shí), 水中爆炸的殉爆中心距離約為10 cm, 安全中心距離約為12 cm, 距離很近。根據(jù)水中爆炸的殉爆相似律可知, 大藥量的水中兵器戰(zhàn)斗部在水下戰(zhàn)場(chǎng)環(huán)境下的安全性也很高。

    5 結(jié)束語(yǔ)

    本文選取了水下典型的主發(fā)裝藥與被發(fā)裝藥模型, 采用有限元軟件ANSYS/LS-DYNA對(duì)典型帶殼裝藥的水下殉爆進(jìn)行了數(shù)值仿真, 并加工了試驗(yàn)樣彈, 開(kāi)展了水下殉爆的試驗(yàn)研究, 得到了典型帶殼裝藥的水下殉爆距離與安全距離。數(shù)值仿真結(jié)果與試驗(yàn)結(jié)果較為吻合, 表明數(shù)值仿真模型能夠有效描述帶殼裝藥的水下殉爆情況, 試驗(yàn)中由沖擊波壓力和氣泡周期來(lái)判斷被發(fā)裝藥是否殉爆是可行的。由水下殉爆的理論分析, 可以預(yù)計(jì)大藥量水中兵器戰(zhàn)斗部的殉爆安全性。

    本文僅選取了典型的水下裝藥模型, 有關(guān)不同裝藥品種、不同裝藥形式與起爆方式、不同裝藥外殼的水下殉爆現(xiàn)象與規(guī)律的研究尚沒(méi)有涉及, 有待后續(xù)深入研究。

    [1] 章冠人, 陳大年. 凝聚炸藥起爆動(dòng)力學(xué)[M]. 北京: 國(guó)防工業(yè)出版社, 1991.

    [2] Howe P M, Huang Y K, Arbuckle A L. A Numerical Study of Detonation Propagation Between Munitions[C]//Pro- cee-dings of the 7th Symposium(International) on Detonation, 1982.

    [3] Lu J P, Lochert I J, Kennedy D L, et al. Simulation of Sy- mpathetic Reaction Rests for PBXN2109[C]//Proceedings of 13th Symposium(International) on Detonation, 2006.

    [4] Fisher S, Baker E L, Wells L, et al. XM982 Excalibur Sympathetic Detonation Modeling and Experimentation

    [C]//Insensitive Munitions & Energetic Materials Te- chnology Symposium, 2006.

    [5] 陳朗, 王晨, 魯建英, 等. 炸藥殉爆實(shí)驗(yàn)和數(shù)值模擬[C] //第四屆爆炸計(jì)算力學(xué)會(huì)議論文集, 2008.

    [6] 李錚, 項(xiàng)續(xù)章. 各種炸藥的殉爆距離[J]. 爆炸與沖擊, 1994, 14(3): 231-241.

    [7] 王晨. 殼裝炸藥殉爆實(shí)驗(yàn)和數(shù)值模擬[J]. 爆炸與沖擊, 2010, 30(2): 152-158.Wang Chen. Experiments and Numerical Simulations of Sympathetic Detonation of Explosives in Shell[J]. Explosion and Shock Waves, 2010, 30(2): 152-158.

    [8] 周保順. 非均質(zhì)炸藥殉爆試驗(yàn)數(shù)值模擬[J]. 彈箭與制導(dǎo)學(xué)報(bào), 2009, 29(5): 145-148.Zhou Bao-shun. Numerical Simulation of Sympathetic Detonation of Heterogeneous Condensed Explosives[J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2009, 29(5): 145-148.

    [9] 張勝民. 基于有限元軟件ANSYS7.0的結(jié)構(gòu)分析[M]. 北京: 清華大學(xué)出版社, 2003.

    [10] ANSYS股份有限公司北京辦事處. ANSYS/LS-DYNA算法基礎(chǔ)和使用方法[M], 2000.

    [11] 劉銳. 不同條件下水中爆炸沖擊壓力場(chǎng)分布特性研究[D]. 北京: 北京理工大學(xué), 2010.

    [12] 美國(guó)ANSYS公司. ANSYS/LS-DYNA User′s Manual[M], 2003.

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

    Research on Law of Sympathetic Detonation of Typical Charge Subjected to Underwater Explosion

    LU Zhong-bao,HU Hong-wei,LIU Riu,YANG Fan

    (1. The 705 Research Institute, China Shipbuilding Industry Corporation, Xi′an 710075, China; 2. Xi′an Modern Chemistry Research Institute, Xi′an 710065, China)

    Aiming at typical charge in shell, a finite element analysis (FEA) model of sympathetic detonation of typical charge subjected to underwater explosion is established using the simulation software ANSYS/LS_DYNA. The distance of sympathetic detonation and the distance of safety are achieved by calculation, and a prototype ammunition is manufactured for underwater sympathetic detonation experiment. The experimental result indicates that the FEA model can describe the underwater sympathetic detonation of charge in shell effectively, and that explosive pressure and bubble period can be used to estimate the occurrence of sympathetic detonation of acceptor charge. Based on the theoretical analysis of underwater sympathetic detonation, the warhead security of typical underwater weapon is estimated according to the experiment result. This study may give a reference for the safety design and the destroy ability evaluation of underwater weapon warhead.

    underwater weapon; warhead; underwater explosion; sympathetic detonation; finite element analysis(FEA) model; security

    TJ630.1; TP391.9

    A

    1673-1948(2014)03-0230-06

    2013-12-16;

    2014-01-27.

    魯忠寶(1978-), 男, 碩士, 高級(jí)工程師, 主要從事水中兵器戰(zhàn)斗部的研究與設(shè)計(jì).

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機(jī)模型
    提煉模型 突破難點(diǎn)
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    在线观看一区二区三区激情| 免费高清在线观看视频在线观看| 69精品国产乱码久久久| 中文字幕色久视频| 久久热在线av| 性色avwww在线观看| 丝袜脚勾引网站| 卡戴珊不雅视频在线播放| 丰满迷人的少妇在线观看| 在线观看免费高清a一片| 亚洲国产色片| 午夜免费观看性视频| 人人妻人人澡人人爽人人夜夜| 亚洲av欧美aⅴ国产| 亚洲av日韩在线播放| 亚洲男人天堂网一区| 国产成人精品一,二区| 91精品三级在线观看| 老汉色av国产亚洲站长工具| 麻豆av在线久日| 亚洲欧美精品自产自拍| 免费高清在线观看日韩| 在线观看一区二区三区激情| 妹子高潮喷水视频| 亚洲国产精品成人久久小说| 亚洲婷婷狠狠爱综合网| 亚洲av电影在线观看一区二区三区| 十分钟在线观看高清视频www| 午夜福利网站1000一区二区三区| 亚洲精品美女久久久久99蜜臀 | 亚洲少妇的诱惑av| av免费在线看不卡| 亚洲综合色惰| 少妇的丰满在线观看| 精品一区二区免费观看| 国产欧美日韩综合在线一区二区| 欧美日韩视频精品一区| 亚洲精品自拍成人| 国产福利在线免费观看视频| 永久网站在线| 精品人妻偷拍中文字幕| 这个男人来自地球电影免费观看 | 黄频高清免费视频| 日韩人妻精品一区2区三区| 日韩不卡一区二区三区视频在线| 多毛熟女@视频| 久久韩国三级中文字幕| 男人爽女人下面视频在线观看| 少妇熟女欧美另类| 天天躁夜夜躁狠狠躁躁| 少妇精品久久久久久久| 少妇精品久久久久久久| 国产精品香港三级国产av潘金莲 | 精品国产国语对白av| 人人澡人人妻人| 国产综合精华液| 亚洲欧洲国产日韩| 日韩精品有码人妻一区| 日韩精品有码人妻一区| 男的添女的下面高潮视频| 最近最新中文字幕大全免费视频 | 国产日韩欧美亚洲二区| 日本欧美国产在线视频| 我的亚洲天堂| 日本爱情动作片www.在线观看| 免费av中文字幕在线| 一边亲一边摸免费视频| 七月丁香在线播放| 国产成人精品婷婷| 亚洲综合色惰| 午夜免费鲁丝| 啦啦啦在线观看免费高清www| 91国产中文字幕| 久久国产精品大桥未久av| 亚洲少妇的诱惑av| 成人二区视频| 亚洲国产色片| 日韩免费高清中文字幕av| 99久久综合免费| 国精品久久久久久国模美| 黄色毛片三级朝国网站| 大片电影免费在线观看免费| 可以免费在线观看a视频的电影网站 | 天天操日日干夜夜撸| 国产成人精品久久二区二区91 | 岛国毛片在线播放| 麻豆av在线久日| 亚洲一级一片aⅴ在线观看| av国产精品久久久久影院| 精品久久久精品久久久| 国产在线一区二区三区精| 久久精品夜色国产| 国产精品免费视频内射| 久久精品国产a三级三级三级| 欧美日韩亚洲高清精品| 少妇的逼水好多| 国产精品不卡视频一区二区| videosex国产| 成人国产麻豆网| 久久综合国产亚洲精品| 777久久人妻少妇嫩草av网站| 九九爱精品视频在线观看| 欧美+日韩+精品| 寂寞人妻少妇视频99o| 国产国语露脸激情在线看| 中国国产av一级| 国产精品嫩草影院av在线观看| 欧美成人精品欧美一级黄| 久久国产精品男人的天堂亚洲| 伦精品一区二区三区| 中文天堂在线官网| 啦啦啦在线观看免费高清www| 人人妻人人添人人爽欧美一区卜| 啦啦啦中文免费视频观看日本| 午夜福利视频精品| 男人舔女人的私密视频| 精品人妻熟女毛片av久久网站| 1024视频免费在线观看| av网站在线播放免费| 9191精品国产免费久久| 久久久久久久大尺度免费视频| 亚洲国产日韩一区二区| 午夜福利在线观看免费完整高清在| 高清黄色对白视频在线免费看| 不卡视频在线观看欧美| 国产av一区二区精品久久| 久热这里只有精品99| 免费少妇av软件| av在线老鸭窝| 欧美精品高潮呻吟av久久| 亚洲av综合色区一区| 日韩中文字幕欧美一区二区 | 色网站视频免费| 大话2 男鬼变身卡| 青草久久国产| 美女主播在线视频| 一边亲一边摸免费视频| 国产免费一区二区三区四区乱码| 日韩中文字幕欧美一区二区 | 精品第一国产精品| 99久国产av精品国产电影| 久久久久久人人人人人| 亚洲国产毛片av蜜桃av| 不卡视频在线观看欧美| 欧美精品人与动牲交sv欧美| 99久久中文字幕三级久久日本| 在线观看美女被高潮喷水网站| 最近2019中文字幕mv第一页| av线在线观看网站| 国产精品二区激情视频| 国产成人精品久久久久久| 少妇人妻 视频| 黑人猛操日本美女一级片| 精品国产国语对白av| 国产老妇伦熟女老妇高清| 亚洲五月色婷婷综合| 亚洲欧美色中文字幕在线| 美女高潮到喷水免费观看| av在线观看视频网站免费| 777米奇影视久久| 日日摸夜夜添夜夜爱| 日韩中文字幕欧美一区二区 | 国产精品久久久久久av不卡| 纯流量卡能插随身wifi吗| 精品视频人人做人人爽| 久久久久久久久久久免费av| 久久免费观看电影| 日本vs欧美在线观看视频| 日本av手机在线免费观看| 久久久久网色| 国产亚洲精品第一综合不卡| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲成国产人片在线观看| av女优亚洲男人天堂| 99国产精品免费福利视频| 精品一区二区三区四区五区乱码 | www.精华液| 久热久热在线精品观看| av.在线天堂| 欧美精品av麻豆av| 久久久国产欧美日韩av| 婷婷色综合大香蕉| 中文字幕精品免费在线观看视频| 日韩,欧美,国产一区二区三区| av片东京热男人的天堂| 婷婷成人精品国产| 亚洲经典国产精华液单| 国产一区二区在线观看av| 少妇猛男粗大的猛烈进出视频| 婷婷成人精品国产| 久久亚洲国产成人精品v| 最近最新中文字幕大全免费视频 | 亚洲成色77777| 街头女战士在线观看网站| 欧美另类一区| 亚洲精品日本国产第一区| 国产乱人偷精品视频| 五月开心婷婷网| 男人操女人黄网站| 成人18禁高潮啪啪吃奶动态图| 老熟女久久久| 亚洲美女黄色视频免费看| 国产在视频线精品| 大香蕉久久网| 中国国产av一级| 青春草国产在线视频| 欧美97在线视频| 天天操日日干夜夜撸| 免费av中文字幕在线| 国产精品免费视频内射| 国产精品麻豆人妻色哟哟久久| 女性生殖器流出的白浆| 日本色播在线视频| 午夜免费鲁丝| 秋霞在线观看毛片| 亚洲精品久久成人aⅴ小说| 最近最新中文字幕大全免费视频 | 高清视频免费观看一区二区| 欧美精品人与动牲交sv欧美| 免费久久久久久久精品成人欧美视频| 美国免费a级毛片| 成年av动漫网址| 亚洲国产精品999| 亚洲av欧美aⅴ国产| 丰满乱子伦码专区| 午夜激情av网站| 在线观看一区二区三区激情| 欧美激情高清一区二区三区 | 菩萨蛮人人尽说江南好唐韦庄| 9色porny在线观看| 亚洲欧美一区二区三区国产| 午夜福利乱码中文字幕| 亚洲第一区二区三区不卡| 搡女人真爽免费视频火全软件| 国产成人精品婷婷| 婷婷色综合www| 99久久人妻综合| 亚洲经典国产精华液单| 视频区图区小说| 在线观看美女被高潮喷水网站| 欧美日韩精品网址| 男女午夜视频在线观看| 成年美女黄网站色视频大全免费| 在线观看三级黄色| 午夜激情av网站| 午夜福利在线免费观看网站| 免费观看性生交大片5| 日韩视频在线欧美| av片东京热男人的天堂| 免费黄网站久久成人精品| 男女边摸边吃奶| 欧美日韩国产mv在线观看视频| 熟女av电影| 久久久久国产一级毛片高清牌| 飞空精品影院首页| freevideosex欧美| 欧美日韩精品网址| 高清黄色对白视频在线免费看| 男男h啪啪无遮挡| 赤兔流量卡办理| 色吧在线观看| 欧美成人午夜免费资源| 中文字幕人妻熟女乱码| 一本大道久久a久久精品| 国产在线视频一区二区| 欧美少妇被猛烈插入视频| 极品少妇高潮喷水抽搐| 国产国语露脸激情在线看| 三级国产精品片| 成年美女黄网站色视频大全免费| 天天躁日日躁夜夜躁夜夜| 亚洲,欧美,日韩| 日韩av免费高清视频| 黑人巨大精品欧美一区二区蜜桃| 日本色播在线视频| 国产无遮挡羞羞视频在线观看| 最近最新中文字幕大全免费视频 | 天堂8中文在线网| 亚洲精品久久成人aⅴ小说| 亚洲国产av影院在线观看| 热99国产精品久久久久久7| 美女午夜性视频免费| 天天躁夜夜躁狠狠久久av| 亚洲熟女精品中文字幕| 啦啦啦中文免费视频观看日本| 欧美亚洲日本最大视频资源| 女性生殖器流出的白浆| 亚洲精品国产一区二区精华液| 一级a爱视频在线免费观看| 国产精品国产av在线观看| 亚洲美女黄色视频免费看| 国产综合精华液| 亚洲精品国产av蜜桃| 成人亚洲欧美一区二区av| 搡女人真爽免费视频火全软件| 久久国产亚洲av麻豆专区| 国产亚洲av片在线观看秒播厂| 亚洲第一区二区三区不卡| 亚洲av福利一区| 黄色视频在线播放观看不卡| 亚洲视频免费观看视频| 日本av免费视频播放| 日日爽夜夜爽网站| 亚洲精品一区蜜桃| 中文字幕另类日韩欧美亚洲嫩草| 777米奇影视久久| 午夜激情av网站| 国产精品二区激情视频| 菩萨蛮人人尽说江南好唐韦庄| 日本欧美国产在线视频| 永久网站在线| 色婷婷久久久亚洲欧美| 亚洲国产精品一区二区三区在线| 午夜福利影视在线免费观看| 精品国产乱码久久久久久男人| 精品少妇内射三级| 亚洲av电影在线进入| 亚洲精品久久午夜乱码| 十八禁网站网址无遮挡| 久久青草综合色| 人妻一区二区av| www.熟女人妻精品国产| 18+在线观看网站| 国产精品免费大片| 在线观看人妻少妇| 免费不卡的大黄色大毛片视频在线观看| 亚洲图色成人| 精品久久久久久电影网| 亚洲伊人色综图| 国产欧美日韩综合在线一区二区| 婷婷成人精品国产| 香蕉国产在线看| 亚洲三区欧美一区| 欧美日韩成人在线一区二区| 久久久精品区二区三区| 黄片无遮挡物在线观看| 美女视频免费永久观看网站| av免费观看日本| 久久精品久久久久久久性| freevideosex欧美| 国产精品一二三区在线看| 一级,二级,三级黄色视频| 欧美人与性动交α欧美软件| 国产免费福利视频在线观看| 我要看黄色一级片免费的| 18+在线观看网站| 一本—道久久a久久精品蜜桃钙片| www.自偷自拍.com| 超色免费av| 丰满少妇做爰视频| 免费观看在线日韩| 国产老妇伦熟女老妇高清| 制服丝袜香蕉在线| 亚洲欧美中文字幕日韩二区| 老鸭窝网址在线观看| 曰老女人黄片| 亚洲av中文av极速乱| 欧美 亚洲 国产 日韩一| 成年人免费黄色播放视频| 午夜久久久在线观看| www.自偷自拍.com| 交换朋友夫妻互换小说| 啦啦啦中文免费视频观看日本| 久久久a久久爽久久v久久| 欧美日韩成人在线一区二区| av国产久精品久网站免费入址| 在线观看免费视频网站a站| 中文字幕人妻熟女乱码| 久久精品国产综合久久久| 久久精品aⅴ一区二区三区四区 | 97在线视频观看| 天天躁夜夜躁狠狠躁躁| 永久免费av网站大全| 一区二区三区四区激情视频| 国产一区有黄有色的免费视频| 制服诱惑二区| 国产片特级美女逼逼视频| 捣出白浆h1v1| 极品人妻少妇av视频| 在线免费观看不下载黄p国产| 在线看a的网站| 卡戴珊不雅视频在线播放| 亚洲视频免费观看视频| 9色porny在线观看| 九草在线视频观看| 久久精品国产综合久久久| 熟女电影av网| av女优亚洲男人天堂| 电影成人av| 少妇 在线观看| 日产精品乱码卡一卡2卡三| 建设人人有责人人尽责人人享有的| 亚洲人成电影观看| 女人被躁到高潮嗷嗷叫费观| 少妇人妻 视频| 伊人久久国产一区二区| av免费在线看不卡| 久久精品熟女亚洲av麻豆精品| 国产精品免费大片| 黄片小视频在线播放| 人人妻人人添人人爽欧美一区卜| 五月开心婷婷网| 丝袜美足系列| 日本免费在线观看一区| 下体分泌物呈黄色| 久久精品亚洲av国产电影网| 黄色 视频免费看| 国产极品粉嫩免费观看在线| 新久久久久国产一级毛片| 中文字幕色久视频| 欧美变态另类bdsm刘玥| 黄片播放在线免费| 岛国毛片在线播放| 日韩免费高清中文字幕av| 久久精品人人爽人人爽视色| 亚洲国产最新在线播放| 永久免费av网站大全| 久久精品夜色国产| 最新的欧美精品一区二区| 热re99久久精品国产66热6| 日韩中字成人| 老女人水多毛片| 一区二区av电影网| 国产免费又黄又爽又色| 中文字幕最新亚洲高清| 岛国毛片在线播放| 我要看黄色一级片免费的| 成人毛片a级毛片在线播放| 精品少妇一区二区三区视频日本电影 | 日韩精品有码人妻一区| 久久久a久久爽久久v久久| 大片电影免费在线观看免费| 亚洲国产精品999| 高清不卡的av网站| av.在线天堂| 777久久人妻少妇嫩草av网站| 免费高清在线观看视频在线观看| 18禁观看日本| 最近最新中文字幕免费大全7| 免费少妇av软件| 男女国产视频网站| 精品少妇黑人巨大在线播放| 精品卡一卡二卡四卡免费| 狠狠婷婷综合久久久久久88av| 亚洲国产色片| 午夜av观看不卡| 老熟女久久久| 熟女少妇亚洲综合色aaa.| 国产精品不卡视频一区二区| 久久韩国三级中文字幕| 高清黄色对白视频在线免费看| 精品少妇一区二区三区视频日本电影 | 欧美xxⅹ黑人| 国产日韩欧美视频二区| 亚洲欧洲国产日韩| 青春草视频在线免费观看| 女的被弄到高潮叫床怎么办| 欧美日韩视频精品一区| 新久久久久国产一级毛片| 久久久国产一区二区| 午夜免费男女啪啪视频观看| 校园人妻丝袜中文字幕| 国产熟女欧美一区二区| 在线免费观看不下载黄p国产| 国产熟女午夜一区二区三区| 日日啪夜夜爽| 免费av中文字幕在线| 亚洲精品美女久久久久99蜜臀 | 天美传媒精品一区二区| 精品国产露脸久久av麻豆| 下体分泌物呈黄色| 哪个播放器可以免费观看大片| 国产在线一区二区三区精| 国产亚洲欧美精品永久| 久久av网站| 91精品国产国语对白视频| 一级毛片黄色毛片免费观看视频| 国产精品99久久99久久久不卡 | 免费观看无遮挡的男女| 寂寞人妻少妇视频99o| 麻豆乱淫一区二区| 菩萨蛮人人尽说江南好唐韦庄| 青春草视频在线免费观看| 国产av码专区亚洲av| 欧美日韩国产mv在线观看视频| 精品福利永久在线观看| 飞空精品影院首页| 大码成人一级视频| 人妻一区二区av| 久久精品国产综合久久久| 久久久久久久久久久久大奶| 中文字幕制服av| 免费不卡的大黄色大毛片视频在线观看| 国产午夜精品一二区理论片| 在线观看免费高清a一片| 亚洲一码二码三码区别大吗| 免费看av在线观看网站| 亚洲成国产人片在线观看| 亚洲精华国产精华液的使用体验| 精品卡一卡二卡四卡免费| 老鸭窝网址在线观看| 中文字幕亚洲精品专区| 2018国产大陆天天弄谢| 午夜福利乱码中文字幕| av在线播放精品| 久久精品国产亚洲av高清一级| 丰满饥渴人妻一区二区三| 亚洲国产av新网站| 两个人免费观看高清视频| 中文字幕av电影在线播放| 制服丝袜香蕉在线| 国产成人欧美| 精品久久蜜臀av无| 国产探花极品一区二区| 春色校园在线视频观看| 两性夫妻黄色片| 丝袜喷水一区| 欧美日韩视频高清一区二区三区二| 亚洲成人av在线免费| 少妇 在线观看| 日本wwww免费看| 国产精品香港三级国产av潘金莲 | 欧美精品av麻豆av| 久久久久久久久久人人人人人人| 久久人人爽人人片av| 亚洲伊人久久精品综合| 18禁国产床啪视频网站| 美国免费a级毛片| 免费人妻精品一区二区三区视频| 亚洲少妇的诱惑av| 成人18禁高潮啪啪吃奶动态图| 最近中文字幕2019免费版| 欧美成人午夜精品| 成人国产av品久久久| 国产免费现黄频在线看| 最近的中文字幕免费完整| 精品国产一区二区久久| 免费看不卡的av| 好男人视频免费观看在线| 亚洲av中文av极速乱| 日日啪夜夜爽| 9色porny在线观看| 午夜福利在线观看免费完整高清在| 超碰97精品在线观看| 日韩大片免费观看网站| 亚洲在久久综合| 女性被躁到高潮视频| 亚洲熟女精品中文字幕| 国产淫语在线视频| 男男h啪啪无遮挡| 亚洲av国产av综合av卡| 涩涩av久久男人的天堂| 天堂8中文在线网| 寂寞人妻少妇视频99o| 欧美 日韩 精品 国产| 国产男女内射视频| 精品酒店卫生间| 最近2019中文字幕mv第一页| 亚洲综合色惰| 在线精品无人区一区二区三| 色婷婷久久久亚洲欧美| 国产欧美日韩综合在线一区二区| 欧美变态另类bdsm刘玥| 精品久久久精品久久久| 国产精品一国产av| 成年动漫av网址| 99久久综合免费| av在线观看视频网站免费| 天天躁狠狠躁夜夜躁狠狠躁| 黄片播放在线免费| 国产女主播在线喷水免费视频网站| 女的被弄到高潮叫床怎么办| 免费av中文字幕在线| 国产 一区精品| 男女午夜视频在线观看| 亚洲精品久久午夜乱码| 老鸭窝网址在线观看| 高清在线视频一区二区三区| 国产精品女同一区二区软件| 亚洲一码二码三码区别大吗| 蜜桃国产av成人99| 免费不卡的大黄色大毛片视频在线观看| 捣出白浆h1v1| 一级a爱视频在线免费观看| 亚洲av男天堂| 十八禁网站网址无遮挡| 欧美精品亚洲一区二区| 中文字幕人妻熟女乱码| 国产成人午夜福利电影在线观看| 日韩欧美精品免费久久| 一区在线观看完整版| 成人午夜精彩视频在线观看| 毛片一级片免费看久久久久| 高清黄色对白视频在线免费看| 亚洲美女搞黄在线观看| 中国国产av一级| 满18在线观看网站| 99久久人妻综合| 欧美日韩视频高清一区二区三区二| 人成视频在线观看免费观看| 精品亚洲乱码少妇综合久久| 天美传媒精品一区二区| 日韩av不卡免费在线播放| 日本91视频免费播放| 高清不卡的av网站| 亚洲人成电影观看| 精品一区二区三卡| 国产免费一区二区三区四区乱码| 香蕉国产在线看| 成人午夜精彩视频在线观看| 夜夜骑夜夜射夜夜干| 国产极品粉嫩免费观看在线| 久久久欧美国产精品| 欧美国产精品va在线观看不卡| 亚洲美女视频黄频|