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

    爆炸荷載多尺度分析方法在仿真計算中的應(yīng)用

    2023-05-10 06:40:14陳龍明李述濤陳葉青寶鑫
    北京理工大學(xué)學(xué)報 2023年5期
    關(guān)鍵詞:爆源子結(jié)構(gòu)尺度

    陳龍明,李述濤,陳葉青,寶鑫

    (1.軍事科學(xué)院 國防工程研究院, 北京 100036;2.清華大學(xué) 土木工程系, 北京 100084)

    研究巖土介質(zhì)中爆炸對深埋地下結(jié)構(gòu)的動力作用以及水中爆炸對艦船、大壩等裝備或工程目標(biāo)的破壞作用時,均涉及到介質(zhì)中封閉爆炸問題研究,爆炸荷載產(chǎn)生的強沖擊作用及介質(zhì)動力響應(yīng)涉及多材料、非線性的復(fù)雜耦合作用[1-3].相關(guān)理論研究較為復(fù)雜,目前多采用顯式有限元方法進行數(shù)值模擬研究[4-5].爆炸近區(qū)產(chǎn)生的幾何大變形是高頻爆炸沖擊波直接作用的結(jié)果[6],數(shù)值計算需要設(shè)置滿足截止頻率的精細化網(wǎng)格來捕捉爆炸荷載中的高頻成分,以提高計算精度.而顯式時域逐步積分的最小時間步長受模型中最小單元尺寸限制,因此整體模型的計算效率較低.

    隨著大當(dāng)量深鉆地武器技術(shù)[7]的不斷發(fā)展,越來越多的研究人員開展了爆炸荷載作用下的大范圍工程場地動力反應(yīng)研究[8].張月輝等[9]使用有限元方法對濱海軟土中的爆炸過程進行了動力響應(yīng)分析,總結(jié)了軟土中不同爆炸量和埋深時爆炸沖擊波的傳播規(guī)律.YANG 等[10]通過建立全耦合的大規(guī)模歐拉-拉格朗日模型研究地表爆炸作用下土-結(jié)構(gòu)-水的動力作用,研究發(fā)現(xiàn)地表爆炸將對淺埋引水箱涵的整體安全性造成極大威脅,并且水位會顯著影響箱涵的動態(tài)特性.爆炸荷載作用下的大范圍工程場地數(shù)值模擬研究中,爆源近區(qū)的精細化網(wǎng)格與大尺度整體模型尺度不協(xié)調(diào)[11],存在著復(fù)雜異形網(wǎng)格過渡的難題.另一方面,較小尺寸的單元降低了滿足顯式時域逐步積分穩(wěn)定性條件的臨界時間步長[12],導(dǎo)致計算成本增加,從而使得整體模型的數(shù)值計算難以高效完成.為解決該問題,李述濤等[13]提出了一種爆炸荷載多尺度分析方法,利用波場分解理論[14]及爆源子結(jié)構(gòu)實現(xiàn)等效爆炸荷載的提取和輸入,分步完成爆炸作用下大范圍中遠場地的動力響應(yīng)分析,在保證計算精度的前提下大幅提高了計算效率.該方法需要在近爆源自由場模型計算基礎(chǔ)上,構(gòu)建爆源子結(jié)構(gòu)模型,加載自由波場位移荷載,求解等效爆炸荷載并施加到大尺度中遠場模型中.雖然在有限元軟件中可以手動完成以上操作,但工作量較大且極易出錯,也不便重復(fù)實施,不利于該方法的推廣應(yīng)用.

    本文簡要介紹爆炸荷載多尺度分析方法的基礎(chǔ)理論,研究該方法在仿真計算中的應(yīng)用技術(shù).基于Python 語言環(huán)境自行開發(fā)一種k 文件處理程序,輔助完成建模和后處理工作,最后使用該方法對單爆源和多爆源2 種算例進行比對分析,給出適用性結(jié)論.

    1 爆炸荷載多尺度分析方法

    爆炸荷載多尺度分析方法基于波場分解理論,利用爆源子結(jié)構(gòu)模型,將爆炸近區(qū)的自由場節(jié)點運動轉(zhuǎn)化為節(jié)點反作用力,提取并施加于大尺度目標(biāo)模型中,在完成等效爆炸荷載加載的同時實現(xiàn)網(wǎng)格尺寸過渡[13].操作流程主要分為3 個步驟,具體見圖1.

    圖1 爆炸荷載多尺度分析方法示意圖(1/2 模型)Fig.1 Schematic diagram of multi-scale analysis method of explosion load (1/2 model)

    步驟1 求解近爆源區(qū)域的節(jié)點運動.建立爆源小尺度模型,設(shè)置精細化網(wǎng)格,獲取近爆源區(qū)域節(jié)點位移時程數(shù)據(jù),并在爆源小尺度模型中選取一定間隔的節(jié)點構(gòu)建爆源子結(jié)構(gòu)模型.最后提取子結(jié)構(gòu)模型外2 層節(jié)點所對應(yīng)的近爆源區(qū)域節(jié)點位移數(shù)據(jù).

    步驟2 通過爆源子結(jié)構(gòu)將位移荷載轉(zhuǎn)化為節(jié)點反作用力荷載.將位移數(shù)據(jù)按空間對應(yīng)位置加載到爆源子結(jié)構(gòu)節(jié)點上,進行動力分析后提取內(nèi)側(cè)2層節(jié)點的節(jié)點反力(reaction force).

    步驟3 對大尺度目標(biāo)模型進行動力計算.建立大尺度網(wǎng)格目標(biāo)模型,在空間坐標(biāo)與子結(jié)構(gòu)重合的節(jié)點上施加上一步提取到的節(jié)點反力,即爆源產(chǎn)生的等效爆炸荷載.通過動力計算獲得大尺度模型在爆炸荷載下的動力響應(yīng).

    2 爆炸荷載多尺度分析方法應(yīng)用

    爆炸荷載多尺度分析方法理論清晰、方法簡便,可在仿真計算中應(yīng)用.以巖土中封閉爆炸為例,以下詳細介紹多尺度分析方法在仿真計算中的應(yīng)用方法.

    2.1 黏彈性人工邊界單元的設(shè)置

    應(yīng)用多尺度分析方法研究巖土中的封閉爆炸問題時,首先需獲得爆源附近的自由場節(jié)點位移荷載,合理地解決爆源小尺度模型的邊界反射問題是實現(xiàn)多尺度方法的基礎(chǔ).本節(jié)主要介紹一種無反射邊界-黏彈性人工邊界在仿真計算中的設(shè)置方法.

    在模擬半無限巖土地下封閉爆炸過程中可通過在截斷邊界位置施加人工邊界條件,防止四周外行的地沖擊波在邊界處反射影響計算結(jié)果.常見的人工邊界方法包括邊界元方法、透射邊界[15]、黏性邊界及黏彈性邊界[16-17]等.基于柱面波理論和球面波動理論[18],劉晶波等[19]及谷音等[20]提出的一致黏彈性人工邊界和等效邊界單元提供了一種操作簡便、效果良好的無反射人工邊界方法.利用普通有限元單元在邊界截斷處構(gòu)造黏彈性邊界單元模擬黏彈性邊界,該人工邊界單元具備良好的吸波效果.李述濤等[21]詳細介紹了黏彈性人工邊界單元在ABAQUS中的實現(xiàn)方法,在仿真計算中使用可以參考此文獻.李述濤等[12,22]介紹了黏彈性人工邊界單元顯式時域逐步積分穩(wěn)定性的分析和改善方法,給出了邊界單元厚度和最小時間步長的確定方法.在近爆源小尺度模型中,通常包含炸藥、空氣和巖土介質(zhì).炸藥和空氣域使用歐拉網(wǎng)格建模,空氣域邊界設(shè)置為無反射邊界;巖土介質(zhì)采用拉格朗日網(wǎng)格建模,在截斷邊界處設(shè)置人工邊界單元.建立黏彈性人工邊界時需在模型最外層延伸一層網(wǎng)格,二維和三維的人工邊界單元如圖2 所示.

    圖2 人工邊界單元示意圖Fig.2 Schematic diagram of artificial boundary element

    黏彈性人工邊界單元的具體設(shè)置如下:

    ①建立實體模型,劃分網(wǎng)格后選中模型截斷邊界處的最外層網(wǎng)格單元,將其劃分到單獨的集合.人工邊界單元應(yīng)盡量劃分為規(guī)則的四邊形/六面體單元.

    ②選中人工邊界的最外層節(jié)點,建立節(jié)點集合,對該節(jié)點的所有運動自由度進行約束.

    ③在線彈性材料模型中設(shè)置邊界單元的等效彈性模量、密度以及泊松比.有限元計算中,若邊界單元質(zhì)量為0,會導(dǎo)致顯式算法穩(wěn)定性下降[12,22],因此將邊界單元的材料密度設(shè)置為較小的值,設(shè)置泊松比為0,同時輸入等效彈性模量.在阻尼設(shè)置模塊下輸入等效阻尼,此處采取以時間為單位的瑞利阻尼形式.等效彈性模量和等效阻尼的計算方法可參考文獻[19 - 20].

    2.2 多尺度方法實施步驟

    以LS-DYNA 為例,多尺度分析方法的實施步驟如下:

    ①建立小尺度的近爆源模型,如圖3 所示.模型包含4 個PART,分別為炸藥、空氣、介質(zhì)(土)和人工邊界單元,流體和固體介質(zhì)使用流固耦合算法實現(xiàn)力學(xué)參量的傳遞.

    圖3 近爆源模型Fig.3 Near-source model

    近爆源模型建立完成后,再確定爆源子結(jié)構(gòu)的空間位置:在炸藥四周選擇一個適當(dāng)?shù)恼叫危ǘS)/正方體(三維)區(qū)域,“跳躍”式選取近爆源自由場模型中的節(jié)點,建立一個由2 層單元和3 層節(jié)點組成的爆源子結(jié)構(gòu)模型,見圖4.找到近爆源模型中與子結(jié)構(gòu)外2 層節(jié)點空間坐標(biāo)一致的節(jié)點編號,按k 文件格式寫入關(guān)鍵字*DATABASE_HISTORY_NODE 模塊下,并在代碼中儲存為列表Node1.該步驟旨在建立需要提取位移的節(jié)點集合,便于在計算輸出文件中快速獲取以上節(jié)點位移數(shù)據(jù).

    圖4 爆源子結(jié)構(gòu)模型Fig.4 Substructure model of explosion source

    ②處理近爆源模型輸出的位移數(shù)據(jù).對近爆源模型進行動力流固耦合計算后,將記錄著Node1 中節(jié)點位移數(shù)據(jù)的結(jié)果輸出文件轉(zhuǎn)化輸入到Excel 文件中.分別將各個分量的位移數(shù)據(jù)儲存在不同的工作表(Sheet)中,二維時每個Excel 文件中有x、y2 個工作表,三維時每個Excel 文件中則有x、y、z3 個工作表.每個工作表中第一列數(shù)據(jù)為時間,其余列則為各個節(jié)點與時間相對應(yīng)的位移,節(jié)點數(shù)據(jù)按節(jié)點編號遞增排列.

    ③確定子結(jié)構(gòu)節(jié)點與近爆源模型對應(yīng)節(jié)點的空間對應(yīng)關(guān)系.在近爆源模型k 文件*NODE 模塊下,記錄列表Node1 中節(jié)點的坐標(biāo)值.然后在爆源子結(jié)構(gòu)模型k 文件的*NODE 模塊下,找到與Node1 中節(jié)點坐標(biāo)值相同的節(jié)點,記錄其編號,儲存為列表NodeSet1.此步驟為后續(xù)加載節(jié)點位移數(shù)據(jù)做準(zhǔn)備.

    ④在爆源子結(jié)構(gòu)模型中加載自由波場運動.固定最內(nèi)側(cè)一層節(jié)點,將近爆源自由場中提取的節(jié)點位移數(shù)據(jù)(Node1)按對應(yīng)關(guān)系加載到列表NodeSet1中的節(jié)點上.二維模型每個節(jié)點有x、y2 個方向分量,三維模型則有x、y、z3 個方向分量.使用關(guān)鍵字*BOUNDARY_PRESCRIBED_MOTION_NODE 對指定編號節(jié)點加載位移數(shù)據(jù),通過關(guān)鍵字*DEFINE_CURVE 定義數(shù)據(jù)的時程曲線.首先將Excel 文件中的節(jié)點時程數(shù)據(jù)按規(guī)則轉(zhuǎn)化到*DEFINE_CURVE_TITLE 模塊中.例如寫入第1 889 593 號節(jié)點x和y方向位移時程曲線的語句如下:

    隨后按照步驟③建立的節(jié)點對應(yīng)關(guān)系寫出*BOUNDARY_PRESCRIBED_MOTION_NODE_ID 關(guān)鍵字及其內(nèi)容.例如對第1 889 593 號節(jié)點加載x方向和y方向位移語句如下:

    將2 個模塊輸出為單獨的k 文件,最后在爆源子結(jié)構(gòu)k 文件中利用關(guān)鍵字*INCLUDE 調(diào)用以上2 個模塊即可實現(xiàn)位移加載.

    ⑤設(shè)置子結(jié)構(gòu)模型中的節(jié)點力輸出.選取子結(jié)構(gòu)內(nèi)2 層節(jié)點,定義為新的節(jié)點集*SET_NODE.利用關(guān)鍵字*DATABASE_NODAL_FORCE_GROUP 定義輸出節(jié)點力的集合與坐標(biāo)系.最后在*DATABASE_NODFOR 中定義節(jié)點力的輸出間隔,完成節(jié)點反力的輸出設(shè)置.求解內(nèi)側(cè)2 層節(jié)點的反作用力時,因讀寫的數(shù)據(jù)量較大,應(yīng)設(shè)置合理的數(shù)據(jù)輸出間隔,在保證輸出曲線光滑的同時減少處理的數(shù)據(jù)量.

    ⑥處理子結(jié)構(gòu)模型輸出的節(jié)點反力數(shù)據(jù).對爆源子結(jié)構(gòu)模型開展動力計算后,利用后處理軟件將記錄著節(jié)點集*SET_NODE 模塊中節(jié)點反力的數(shù)據(jù)文件轉(zhuǎn)化輸入到Excel 文件中,將*SET_NODE 中的節(jié)點儲存為列表Node2.分別將各分量的節(jié)點力數(shù)據(jù)儲存在不同的工作表(Sheet)中,二維時每個Excel 文件中有x、y2 個工作表,三維時每個Excel 文件中則有x、y、z3 個工作表.每個工作表中第一列為時間,其余列則為各個節(jié)點的節(jié)點反力.此節(jié)點反力即下一步加載于介質(zhì)-目標(biāo)模型的等效爆炸荷載.

    ⑦確定大尺度模型中加載等效爆炸荷載的節(jié)點與子結(jié)構(gòu)節(jié)點的空間對應(yīng)關(guān)系.根據(jù)子結(jié)構(gòu)單元尺寸建立大尺度介質(zhì)-結(jié)構(gòu)模型,如圖5 所示,并在截斷位置處設(shè)置人工邊界單元.介質(zhì)-結(jié)構(gòu)模型的網(wǎng)格尺寸與子結(jié)構(gòu)網(wǎng)格尺寸需保持一致.在介質(zhì)-結(jié)構(gòu)模型中尋找與子結(jié)構(gòu)模型內(nèi)側(cè)2 層節(jié)點坐標(biāo)值相同的節(jié)點以加載節(jié)點反力.具體流程如下:首先在子結(jié)構(gòu)模型k 文件的*NODE 模塊下,記錄列表Node2 中節(jié)點的坐標(biāo)值.然后在大尺度介質(zhì)-結(jié)構(gòu)模型文件的*NODE 模塊下,找到與Node2 中節(jié)點坐標(biāo)值相同的節(jié)點,記錄其編號并儲存為列表NodeSet2.此步驟是為后續(xù)加載節(jié)點反力數(shù)據(jù)做準(zhǔn)備.⑧在大尺度介質(zhì)-結(jié)構(gòu)模型中加載等效爆炸荷載.通過關(guān)鍵字*LOAD_NODE_POINT 對列表Node-Set2 中的節(jié)點加載力-時程數(shù)據(jù).通過關(guān)鍵字*DEFINE_CURVE 定義數(shù)據(jù)的時程曲線.首先將Excel 中的節(jié)點時程數(shù)據(jù)按k 文件格式寫入*DEFINE_CURVE_TITLE 模塊中.例如寫入第1 888 886 號節(jié)點x和y方向位移時程曲線的語句如下:

    圖5 介質(zhì)-結(jié)構(gòu)目標(biāo)模型Fig.5 Medium-structure model

    隨后按照步驟⑦建立的節(jié)點對應(yīng)關(guān)系寫出*LOAD_NODE_POINT 關(guān)鍵字及其內(nèi)容.例如對第1 888 886 號節(jié)點加載x方向和y方向節(jié)點力語句如下:

    將2 個模塊輸出為單獨的k 文件,最后在大尺度介質(zhì)-結(jié)構(gòu)k 文件中利用關(guān)鍵字*INCLUDE 調(diào)用以上2 個模塊即可實現(xiàn)等效爆炸荷載加載.

    2.3 編程的處理流程

    由于多尺度方法需要比對和讀寫大量空間節(jié)點數(shù)據(jù)信息,如2.2 節(jié)中的步驟③④⑦⑧,特別是三維模型,涉及的數(shù)據(jù)量龐大,若全部由人工手動操作難度極大.Python 語言簡單易讀、功能強大,因此利用該語言自編程序?qū)Ψ抡嬗嬎汴P(guān)鍵字文件進行自動改寫,可簡單快捷實現(xiàn)上述步驟.為便于理解實施流程的邏輯關(guān)系,結(jié)合上一節(jié)的實施步驟,本節(jié)給出了編程的處理流程,如圖6 所示.

    圖6 編程處理流程圖Fig.6 Programming process flow chart

    處理流程大致分為3 部分,各部分的標(biāo)志環(huán)節(jié)分別為近爆源模型、爆源子結(jié)構(gòu)模型和大尺度模型的有限元分析.在處理流程中,工作量較大的環(huán)節(jié)分別為比對節(jié)點坐標(biāo)值、寫入數(shù)據(jù)以及輸出加載位移/節(jié)點反力的關(guān)鍵字.通過自編程序來解決這些問題,例如可通過條件循環(huán)語句來比對節(jié)點坐標(biāo)值,使用循環(huán)語句批量輸出加載位移/節(jié)點反力的關(guān)鍵字模塊.給出的處理流程中,仍有部分環(huán)節(jié)需手動操作,例如設(shè)置指定節(jié)點的集合等,但這些環(huán)節(jié)的工作量和操作難度都較小.圖6 中人工操作難度大的環(huán)節(jié)都通過自編Python 程序進行處理,極大地降低了多尺度分析方法在有限元軟件中應(yīng)用的操作難度.同時也應(yīng)當(dāng)指出,除Python 語言外,也可利用其他程序語言進行編程處理,例如C 語言、Fortran 語言等,本節(jié)僅提供編程處理的思路.

    3 多尺度方法的應(yīng)用算例

    隧道開挖等工程實踐中需分析地下封閉爆炸對深部地下結(jié)構(gòu)和地表建筑的影響,本節(jié)應(yīng)用多尺度方法分別研究地下結(jié)構(gòu)在單個爆源荷載作用下的動態(tài)響應(yīng)和地表橋梁結(jié)構(gòu)在多個爆源荷載作用下的動態(tài)響應(yīng).

    3.1 算例1:單個爆源作用于地下隧道

    應(yīng)用本文的實現(xiàn)方法建立二維半空間單爆源地下封閉爆炸計算模型,計算區(qū)域為180 m×180 m 的平面,如圖7 所示.對模型進行均勻網(wǎng)格劃分,網(wǎng)格尺寸為25 cm,單元總數(shù)為51.6 萬個.同時使用普通建模方法建立了同規(guī)模的對比計算模型,近爆源區(qū)域采取網(wǎng)格過渡,炸藥與介質(zhì)之間采取流固耦合方法實現(xiàn)相互作用.爆源區(qū)域進行網(wǎng)格細化,遠離爆源區(qū)域采用均勻網(wǎng)格劃分,兩者之間采取漸變網(wǎng)格過渡,最小網(wǎng)格尺寸為2 cm,最大網(wǎng)格尺寸為25 cm,單元總數(shù)為74.8 萬個.爆源位于地下66 m 深處,裝藥為TNT,使用JWL 狀態(tài)方程,材料模型選取為*MAT_HIGH_EXPLOSIVE_BURN.介質(zhì)為彈性材料,材料模型選取為*MAT_ELASTIC,并按2.1 節(jié)的方法設(shè)置人工邊界.具體材料參數(shù)參考文獻[13].

    圖7 單爆源地下封閉爆炸二維模型示意圖Fig.7 2D model of a single-source underground closed explosion

    在相同計算條件下完成計算,普通方法計算耗時322 min,多尺度方法耗時44 min.圖8 為多尺度方法算例和普通建模方法算例的位移場演化云圖.在使用多尺度方法的算例中,應(yīng)力波剛開始傳播時,子結(jié)構(gòu)對應(yīng)的內(nèi)部方形區(qū)域的節(jié)點位移為0,隨著應(yīng)力波向外傳播,兩者位移云圖趨于一致.當(dāng)t= 0.130 s時,除爆心區(qū)域外,應(yīng)力波的位移云圖結(jié)構(gòu)基本相同.當(dāng)t> 0.130 s 時,位移云圖的結(jié)構(gòu)差異已不明顯.

    圖8 二維介質(zhì)-結(jié)構(gòu)模型位移云圖Fig.8 Displacement of 2D medium-structure model

    圖9 為測點A和B在2 種建模方法下測到的節(jié)點位移曲線.圖9 中,兩者的位移時程曲線基本吻合,在單爆源的情況下多尺度分析方法可獲得與普通建模方法相同的計算精度.

    圖9 測點A、B 位移時程曲線對比Fig.9 Comparison of displacement curves of A and B

    3.2 算例2:多個爆源作用于地表橋梁

    應(yīng)用本文的實現(xiàn)方法建立二維半空間多爆源地下封閉爆炸計算模型,計算區(qū)域為110 m×340 m 的平面,地表處為橋梁橫截面,如圖10 所示.對模型進行均勻網(wǎng)格劃分,網(wǎng)格尺寸為25 cm,單元總數(shù)為60.6萬個.同時使用普通建模方法建立了同規(guī)模的對比計算模型.爆源區(qū)域進行網(wǎng)格細化,遠離爆源區(qū)域采用均勻網(wǎng)格劃分,兩者之間采取漸變網(wǎng)格過渡,最小網(wǎng)格尺寸為2 cm,最大網(wǎng)格尺寸為25 cm,單元總數(shù)為130.3 萬個.橋梁結(jié)構(gòu)材料模型為*MAT_RHT,其余材料模型及計算條件同算例1.

    圖10 多爆源地下封閉爆炸二維模型示意圖Fig.10 2D model of multiple sources underground closed explosion

    在相同計算條件下完成計算,普通方法計算耗時537 min,多尺度方法耗時31 min.圖11 為分別為多尺度方法算例和普通建模方法算例的位移場演化云圖.在使用多尺度方法的算例中,應(yīng)力波剛開始傳播時,子結(jié)構(gòu)對應(yīng)的內(nèi)部方形區(qū)域的節(jié)點位移為0,隨著應(yīng)力波向外傳播,兩者的位移云圖差異減少.應(yīng)力波作用于橋梁結(jié)構(gòu)時,產(chǎn)生反射、繞射等現(xiàn)象,橋梁結(jié)構(gòu)因爆炸荷載作用而產(chǎn)生不均勻的位移分布.由圖11(c)和11(f)可知,橋梁產(chǎn)生的位移響應(yīng)結(jié)果基本一致,表明在多爆源的算例中多尺度方法與普通方法的計算效果吻合得較好.

    圖11 二維介質(zhì)-結(jié)構(gòu)模型位移云圖Fig.11 Displacement of 2D medium-structure model

    圖12 為測點C和D在2 種建模方法下測到的節(jié)點位移曲線.圖12 中,兩者的位移時程曲線基本吻合,在多爆源的情況下多尺度分析方法可獲得與普通建模方法相同的計算精度.

    圖12 測點C、D 點位移時程曲線對比Fig.12 Comparison of displacement curves of C and D

    多尺度方法在單爆源和多爆源的算例中均能獲得與普通方法基本相同的計算效果,且計算效率遠高于普通建模方法.爆源條件相同時,利用本文提出的操作方法可避免重復(fù)計算近場爆炸荷載,便于獲得大批次相似算例的計算結(jié)果.以上的研究結(jié)果表明,與普通建模方法相比,在有限元軟件中應(yīng)用多尺度方法計算地下封閉爆炸模型具有突出的優(yōu)勢.

    4 結(jié) 論

    針對基于爆源子結(jié)構(gòu)的爆炸荷載多尺度分析方法在仿真計算中存在一定的應(yīng)用困難問題,本文簡要介紹了該方法的基本理論,研究了該方法在仿真計算中的應(yīng)用技術(shù),對單爆源和多爆源2 種算例進行了對比分析.具體結(jié)論如下:

    ①基于爆源子結(jié)構(gòu)的爆炸荷載多尺度分析方法可以較好地應(yīng)用于大當(dāng)量爆炸荷載作用下的大范圍工程場地動力響應(yīng)計算.在保證計算精度的前提下,計算效率顯著優(yōu)于炸藥-介質(zhì)-目標(biāo)整體計算模型.

    ②對小尺度爆源模型進行一次計算,即可形成一個參數(shù)化的標(biāo)準(zhǔn)荷載,快速加載到介質(zhì)中的任意位置.對于爆炸荷載位置變化或多爆源情況下的大范圍工程場地動力計算問題,不必重復(fù)建模,避免對炸藥爆炸過程進行低效率的流固耦合計算,大幅提高計算效率.

    ③建模計算中,需要應(yīng)用基于Python 的自編程序?qū)崿F(xiàn)空間節(jié)點坐標(biāo)和數(shù)據(jù)信息的比對和讀寫.該程序邏輯清晰,可快速讀寫大批量數(shù)據(jù)文件,減少了人工操作可能帶來的數(shù)據(jù)錯誤,極大地提高了模型前處理效率.

    猜你喜歡
    爆源子結(jié)構(gòu)尺度
    爆源和測點深度對水下爆炸沖擊波載荷的影響
    完全對換網(wǎng)絡(luò)的結(jié)構(gòu)連通度和子結(jié)構(gòu)連通度
    財產(chǎn)的五大尺度和五重應(yīng)對
    中型浮動沖擊平臺試驗爆源定位方法研究
    中國測試(2018年10期)2018-11-17 01:58:24
    鋼框架腹板雙角鋼連接梁柱子結(jié)構(gòu)抗倒塌性能分析
    宇宙的尺度
    太空探索(2016年5期)2016-07-12 15:17:55
    不同攻角對實船爆炸試驗沖擊響應(yīng)影響研究
    基于子結(jié)構(gòu)的柴油機曲軸有限元建模方法研究
    某型艦抗沖擊試驗爆源定位方法研究*
    9
    国产成人精品在线电影| 久久久久精品国产欧美久久久| 91在线观看av| 亚洲精品中文字幕在线视频| 久久中文看片网| 国产精品自产拍在线观看55亚洲| 如日韩欧美国产精品一区二区三区| 91麻豆av在线| 国产精品野战在线观看| 午夜视频精品福利| 国产一区二区三区综合在线观看| 又大又爽又粗| 亚洲第一电影网av| 亚洲色图av天堂| 波多野结衣av一区二区av| 久久久水蜜桃国产精品网| 亚洲五月天丁香| 国产成人av教育| 精品久久久久久久久久免费视频| 亚洲三区欧美一区| 男女做爰动态图高潮gif福利片 | 亚洲成a人片在线一区二区| 正在播放国产对白刺激| 久久久久国产精品人妻aⅴ院| 神马国产精品三级电影在线观看 | 天天躁狠狠躁夜夜躁狠狠躁| 大陆偷拍与自拍| 一级作爱视频免费观看| 国内精品久久久久精免费| 欧美日韩精品网址| 天天躁狠狠躁夜夜躁狠狠躁| 淫秽高清视频在线观看| 三级毛片av免费| 高潮久久久久久久久久久不卡| 午夜激情av网站| 亚洲美女黄片视频| 嫩草影院精品99| 一级作爱视频免费观看| 国产精品野战在线观看| 亚洲av成人不卡在线观看播放网| 久久久久久久午夜电影| 无人区码免费观看不卡| 日日夜夜操网爽| 国产成人精品久久二区二区免费| 亚洲三区欧美一区| 久久精品国产清高在天天线| 在线观看午夜福利视频| 久久这里只有精品19| 中出人妻视频一区二区| 亚洲精品av麻豆狂野| 欧美乱码精品一区二区三区| 午夜福利18| videosex国产| 无人区码免费观看不卡| 欧美最黄视频在线播放免费| 亚洲成人免费电影在线观看| 亚洲,欧美精品.| 人妻久久中文字幕网| netflix在线观看网站| 999精品在线视频| 国产欧美日韩一区二区三| 69av精品久久久久久| 97超级碰碰碰精品色视频在线观看| 看片在线看免费视频| 色哟哟哟哟哟哟| 久久精品国产清高在天天线| av超薄肉色丝袜交足视频| 精品欧美一区二区三区在线| 91九色精品人成在线观看| 97碰自拍视频| 91麻豆精品激情在线观看国产| 亚洲av熟女| 欧美中文日本在线观看视频| 超碰成人久久| 日韩精品免费视频一区二区三区| 日本 欧美在线| 自拍欧美九色日韩亚洲蝌蚪91| 美国免费a级毛片| 亚洲伊人色综图| 国产三级黄色录像| 国产精品秋霞免费鲁丝片| 窝窝影院91人妻| 精品国产超薄肉色丝袜足j| 色综合亚洲欧美另类图片| 一个人观看的视频www高清免费观看 | 久久久精品国产亚洲av高清涩受| 搡老熟女国产l中国老女人| 嫩草影视91久久| 国产一区二区在线av高清观看| 一本久久中文字幕| 中文字幕久久专区| 韩国精品一区二区三区| 国内毛片毛片毛片毛片毛片| 欧美丝袜亚洲另类 | 老汉色av国产亚洲站长工具| 人人妻人人澡人人看| 午夜免费鲁丝| 午夜视频精品福利| 亚洲九九香蕉| 欧美另类亚洲清纯唯美| 国产亚洲av高清不卡| 精品国产美女av久久久久小说| 日本 欧美在线| 日本三级黄在线观看| 曰老女人黄片| 1024香蕉在线观看| 一区在线观看完整版| 黄色a级毛片大全视频| 日本黄色视频三级网站网址| 91麻豆精品激情在线观看国产| 在线观看午夜福利视频| 免费在线观看日本一区| 乱人伦中国视频| 成在线人永久免费视频| 欧美日韩亚洲国产一区二区在线观看| 在线观看免费视频日本深夜| 99在线人妻在线中文字幕| 男人操女人黄网站| 国产亚洲精品av在线| 亚洲伊人色综图| 99国产综合亚洲精品| 国产精品亚洲av一区麻豆| 国产伦一二天堂av在线观看| 久热这里只有精品99| 午夜福利在线观看吧| 日本在线视频免费播放| 身体一侧抽搐| 成人亚洲精品一区在线观看| 91麻豆av在线| 成人特级黄色片久久久久久久| 国产精品免费一区二区三区在线| 国产极品粉嫩免费观看在线| 欧美日本视频| 午夜a级毛片| 一级,二级,三级黄色视频| www.熟女人妻精品国产| 12—13女人毛片做爰片一| 亚洲精品国产精品久久久不卡| 在线观看日韩欧美| 亚洲自偷自拍图片 自拍| 啦啦啦免费观看视频1| 男女下面插进去视频免费观看| 老司机午夜福利在线观看视频| av天堂久久9| 精品高清国产在线一区| 精品一区二区三区视频在线观看免费| av在线播放免费不卡| 精品国产美女av久久久久小说| 色综合站精品国产| 亚洲色图 男人天堂 中文字幕| 日韩一卡2卡3卡4卡2021年| www.精华液| 啦啦啦观看免费观看视频高清 | 国产精品av久久久久免费| 黄片大片在线免费观看| svipshipincom国产片| 中文字幕av电影在线播放| 老熟妇乱子伦视频在线观看| 午夜福利18| 国产精品一区二区三区四区久久 | 亚洲视频免费观看视频| 如日韩欧美国产精品一区二区三区| 免费一级毛片在线播放高清视频 | 在线播放国产精品三级| 18禁观看日本| 久久亚洲真实| 成人精品一区二区免费| 一区福利在线观看| 免费在线观看黄色视频的| 在线观看舔阴道视频| 一a级毛片在线观看| 丝袜人妻中文字幕| 99久久99久久久精品蜜桃| 亚洲第一av免费看| tocl精华| 亚洲国产精品合色在线| 亚洲av成人av| 欧美成人免费av一区二区三区| 一级a爱视频在线免费观看| 欧美黑人欧美精品刺激| 日韩一卡2卡3卡4卡2021年| 色老头精品视频在线观看| 99久久久亚洲精品蜜臀av| 变态另类丝袜制服| 美女扒开内裤让男人捅视频| 天天躁夜夜躁狠狠躁躁| 久久精品91无色码中文字幕| 久久精品国产综合久久久| 亚洲自拍偷在线| 夜夜爽天天搞| 亚洲国产日韩欧美精品在线观看 | 高清黄色对白视频在线免费看| 久久影院123| 女人被躁到高潮嗷嗷叫费观| 免费搜索国产男女视频| 色综合婷婷激情| 在线观看免费视频网站a站| 一区二区三区精品91| 啦啦啦韩国在线观看视频| 村上凉子中文字幕在线| 亚洲色图 男人天堂 中文字幕| 91九色精品人成在线观看| 日韩欧美一区视频在线观看| 亚洲成a人片在线一区二区| 亚洲专区国产一区二区| 久久中文字幕一级| 国产高清videossex| 精品电影一区二区在线| 男女下面插进去视频免费观看| 夜夜躁狠狠躁天天躁| 午夜免费成人在线视频| 给我免费播放毛片高清在线观看| 岛国在线观看网站| 午夜福利免费观看在线| 久久国产精品人妻蜜桃| 在线视频色国产色| 国产国语露脸激情在线看| 欧美激情久久久久久爽电影 | 精品国产国语对白av| avwww免费| 久9热在线精品视频| 91九色精品人成在线观看| 精品国产美女av久久久久小说| 精品国产一区二区久久| 日本黄色视频三级网站网址| 99国产精品一区二区蜜桃av| 脱女人内裤的视频| 久久精品国产亚洲av香蕉五月| av福利片在线| 女同久久另类99精品国产91| 男女之事视频高清在线观看| 夜夜躁狠狠躁天天躁| 美女免费视频网站| 一进一出抽搐动态| 国产午夜福利久久久久久| 欧美日韩乱码在线| 精品久久久久久久人妻蜜臀av | 国产亚洲欧美98| 侵犯人妻中文字幕一二三四区| 色综合站精品国产| 久久人人爽av亚洲精品天堂| avwww免费| 国产精品美女特级片免费视频播放器 | 久久精品国产清高在天天线| 婷婷丁香在线五月| 精品人妻1区二区| 久久精品成人免费网站| 亚洲国产精品合色在线| 天天添夜夜摸| 脱女人内裤的视频| 久久中文看片网| 欧美中文日本在线观看视频| 嫩草影视91久久| 久久精品91无色码中文字幕| 亚洲最大成人中文| 露出奶头的视频| 国产日韩一区二区三区精品不卡| 日本 av在线| 美女 人体艺术 gogo| 精品少妇一区二区三区视频日本电影| 一夜夜www| 可以在线观看毛片的网站| 精品人妻在线不人妻| 精品一区二区三区四区五区乱码| 老司机靠b影院| 欧美在线一区亚洲| 欧美绝顶高潮抽搐喷水| 91成年电影在线观看| 级片在线观看| 国产又色又爽无遮挡免费看| 亚洲一区高清亚洲精品| 黑人巨大精品欧美一区二区蜜桃| 香蕉国产在线看| 激情在线观看视频在线高清| 久久中文字幕人妻熟女| 免费一级毛片在线播放高清视频 | 国产精品野战在线观看| 国产精品二区激情视频| 亚洲第一电影网av| 亚洲国产毛片av蜜桃av| 无人区码免费观看不卡| 亚洲中文字幕日韩| 最近最新免费中文字幕在线| www日本在线高清视频| 女生性感内裤真人,穿戴方法视频| 欧美日韩亚洲国产一区二区在线观看| 色在线成人网| 国产极品粉嫩免费观看在线| 老汉色∧v一级毛片| 黄频高清免费视频| 女人爽到高潮嗷嗷叫在线视频| 人人妻人人爽人人添夜夜欢视频| 亚洲精品国产色婷婷电影| 精品久久久久久,| 免费在线观看黄色视频的| 精品不卡国产一区二区三区| 亚洲精品中文字幕一二三四区| 一级毛片精品| 亚洲熟妇熟女久久| 亚洲在线自拍视频| 欧美激情 高清一区二区三区| 精品欧美一区二区三区在线| 国产男靠女视频免费网站| 亚洲色图 男人天堂 中文字幕| 国产精品98久久久久久宅男小说| 69av精品久久久久久| 香蕉久久夜色| 国产成人精品在线电影| 一级毛片高清免费大全| 日韩三级视频一区二区三区| 欧美一级a爱片免费观看看 | 亚洲国产精品999在线| 精品无人区乱码1区二区| 色哟哟哟哟哟哟| videossex国产| 天天一区二区日本电影三级| 国产伦人伦偷精品视频| 欧美日本视频| 欧美性猛交黑人性爽| 日本一本二区三区精品| 国产精品,欧美在线| 亚洲avbb在线观看| 色噜噜av男人的天堂激情| 97超级碰碰碰精品色视频在线观看| 12—13女人毛片做爰片一| 18禁裸乳无遮挡免费网站照片| 免费看光身美女| 嫩草影院精品99| 日韩一区二区视频免费看| 久久久久久久久久黄片| 国产又黄又爽又无遮挡在线| 久久精品国产亚洲av香蕉五月| 成年女人毛片免费观看观看9| 欧美+亚洲+日韩+国产| 1024手机看黄色片| 1000部很黄的大片| 亚洲国产色片| av在线蜜桃| 国产成人影院久久av| 国产精品国产高清国产av| 极品教师在线免费播放| 在线播放无遮挡| 日本精品一区二区三区蜜桃| 99热这里只有是精品50| 久久久久久伊人网av| 欧美另类亚洲清纯唯美| 久久精品夜夜夜夜夜久久蜜豆| 国产69精品久久久久777片| 久久精品人妻少妇| 午夜激情福利司机影院| 久久久午夜欧美精品| 国产成人aa在线观看| 人妻夜夜爽99麻豆av| 男插女下体视频免费在线播放| 波野结衣二区三区在线| 午夜福利在线观看吧| 老熟妇乱子伦视频在线观看| 欧美黑人巨大hd| 日韩,欧美,国产一区二区三区 | 亚洲av五月六月丁香网| 男女视频在线观看网站免费| 黄色日韩在线| 在线观看午夜福利视频| 亚洲中文字幕日韩| 久久国内精品自在自线图片| 国产一级毛片七仙女欲春2| 一区二区三区四区激情视频 | 欧美最黄视频在线播放免费| 少妇的逼好多水| 亚洲欧美日韩高清在线视频| 哪里可以看免费的av片| 欧美一级a爱片免费观看看| 免费不卡的大黄色大毛片视频在线观看 | 久久精品国产亚洲av天美| 少妇被粗大猛烈的视频| 亚洲不卡免费看| 国产高清激情床上av| 色av中文字幕| 亚洲美女视频黄频| 精品久久久久久久久久久久久| 久久精品国产亚洲av香蕉五月| 国产精品久久电影中文字幕| 久久婷婷人人爽人人干人人爱| 精品久久久久久成人av| 久久婷婷人人爽人人干人人爱| 免费黄网站久久成人精品| 欧美日韩精品成人综合77777| 女人被狂操c到高潮| 国产精品福利在线免费观看| 免费看av在线观看网站| 少妇被粗大猛烈的视频| 校园春色视频在线观看| 日韩欧美一区二区三区在线观看| 老师上课跳d突然被开到最大视频| 99热网站在线观看| 成人美女网站在线观看视频| 变态另类丝袜制服| 人人妻人人看人人澡| a级毛片免费高清观看在线播放| 波多野结衣巨乳人妻| 欧美性猛交黑人性爽| 日日摸夜夜添夜夜添小说| 国产一区二区在线av高清观看| 国产欧美日韩精品一区二区| 综合色av麻豆| 特级一级黄色大片| 久久6这里有精品| 成人国产综合亚洲| 97热精品久久久久久| 国产老妇女一区| 桃红色精品国产亚洲av| 12—13女人毛片做爰片一| 亚洲成人久久爱视频| 亚洲专区中文字幕在线| 国产亚洲欧美98| 国产 一区精品| 特大巨黑吊av在线直播| www.色视频.com| 国产 一区 欧美 日韩| 国产精品人妻久久久影院| 又黄又爽又刺激的免费视频.| 亚州av有码| 亚洲人成网站高清观看| 国产精品日韩av在线免费观看| 小蜜桃在线观看免费完整版高清| 俄罗斯特黄特色一大片| 黄色一级大片看看| 99riav亚洲国产免费| 亚洲成人免费电影在线观看| 亚洲国产色片| 亚洲真实伦在线观看| 99热这里只有是精品在线观看| 午夜福利18| 精品午夜福利视频在线观看一区| 国产精品av视频在线免费观看| 久久精品国产亚洲av天美| 久久久久久伊人网av| 男女那种视频在线观看| 国产精品日韩av在线免费观看| 日韩在线高清观看一区二区三区 | 免费观看精品视频网站| 久久久久久久午夜电影| 欧美绝顶高潮抽搐喷水| 久久午夜福利片| 亚洲第一区二区三区不卡| 午夜精品在线福利| 波多野结衣高清无吗| 亚洲一区高清亚洲精品| 久久亚洲精品不卡| 精品久久久久久久久久免费视频| 免费av毛片视频| 久久久久精品国产欧美久久久| 国产精品免费一区二区三区在线| 国产欧美日韩精品一区二区| 免费大片18禁| 亚洲精品一卡2卡三卡4卡5卡| 高清毛片免费观看视频网站| 国产精品综合久久久久久久免费| 麻豆av噜噜一区二区三区| 97人妻精品一区二区三区麻豆| 九九爱精品视频在线观看| 色吧在线观看| 日日夜夜操网爽| 99九九线精品视频在线观看视频| 欧美不卡视频在线免费观看| 高清毛片免费观看视频网站| 午夜福利高清视频| 免费看日本二区| 热99在线观看视频| 日韩大尺度精品在线看网址| 欧美性猛交╳xxx乱大交人| 特级一级黄色大片| 日韩欧美精品免费久久| 亚洲专区中文字幕在线| 如何舔出高潮| 97超视频在线观看视频| 免费一级毛片在线播放高清视频| 日本免费a在线| .国产精品久久| 久久午夜亚洲精品久久| 午夜福利在线在线| 一卡2卡三卡四卡精品乱码亚洲| 99久久中文字幕三级久久日本| 日韩大尺度精品在线看网址| a级毛片a级免费在线| 日韩,欧美,国产一区二区三区 | 日韩一区二区视频免费看| 国内毛片毛片毛片毛片毛片| 又紧又爽又黄一区二区| 可以在线观看的亚洲视频| 久久人人精品亚洲av| 日本三级黄在线观看| 长腿黑丝高跟| 国产午夜精品论理片| 亚洲精品影视一区二区三区av| 999久久久精品免费观看国产| 国内精品久久久久久久电影| 日韩精品有码人妻一区| 人妻丰满熟妇av一区二区三区| 亚洲狠狠婷婷综合久久图片| 国产黄片美女视频| 人妻少妇偷人精品九色| 少妇人妻一区二区三区视频| 观看免费一级毛片| 国产精品自产拍在线观看55亚洲| 国产精品98久久久久久宅男小说| 国产人妻一区二区三区在| 亚洲自拍偷在线| 精品久久久久久久久久久久久| 久久久久久国产a免费观看| 天天一区二区日本电影三级| 成年人黄色毛片网站| 午夜免费成人在线视频| 国产av一区在线观看免费| 国产成人福利小说| 国内精品一区二区在线观看| 亚洲最大成人手机在线| 自拍偷自拍亚洲精品老妇| 国产精品永久免费网站| 精品久久国产蜜桃| 国产毛片a区久久久久| 男人和女人高潮做爰伦理| 亚洲乱码一区二区免费版| 精品久久久久久久久久久久久| 色综合站精品国产| 亚洲经典国产精华液单| 狂野欧美白嫩少妇大欣赏| 久久精品久久久久久噜噜老黄 | 亚洲av一区综合| 国产亚洲精品久久久久久毛片| 非洲黑人性xxxx精品又粗又长| 在线观看一区二区三区| 美女被艹到高潮喷水动态| 岛国在线免费视频观看| 国产女主播在线喷水免费视频网站 | 极品教师在线视频| 亚洲国产精品合色在线| 国产精品嫩草影院av在线观看 | 亚洲成人久久爱视频| 午夜福利欧美成人| 乱码一卡2卡4卡精品| 亚洲狠狠婷婷综合久久图片| 欧美日韩综合久久久久久 | 亚洲国产高清在线一区二区三| 国产精品美女特级片免费视频播放器| 午夜a级毛片| 99久久精品一区二区三区| 亚洲欧美日韩卡通动漫| 日韩欧美 国产精品| 成熟少妇高潮喷水视频| 精品久久国产蜜桃| h日本视频在线播放| 狂野欧美白嫩少妇大欣赏| 婷婷六月久久综合丁香| 午夜日韩欧美国产| 在线观看午夜福利视频| 亚洲精品一区av在线观看| 老师上课跳d突然被开到最大视频| 欧美性猛交╳xxx乱大交人| 中文亚洲av片在线观看爽| 麻豆精品久久久久久蜜桃| 国产精品,欧美在线| 国产精品亚洲美女久久久| 亚洲四区av| 精品人妻视频免费看| 天堂网av新在线| 一个人免费在线观看电影| 亚洲国产高清在线一区二区三| 干丝袜人妻中文字幕| 波多野结衣高清作品| 黄色配什么色好看| 日韩 亚洲 欧美在线| 久久99热这里只有精品18| 国产亚洲精品av在线| 99久久精品国产国产毛片| 成人三级黄色视频| 久久久久久久亚洲中文字幕| 91在线精品国自产拍蜜月| 欧美一区二区精品小视频在线| 日韩欧美在线二视频| 欧美日韩中文字幕国产精品一区二区三区| av福利片在线观看| 桃色一区二区三区在线观看| 亚洲自偷自拍三级| 欧美日韩瑟瑟在线播放| 18禁黄网站禁片午夜丰满| www.色视频.com| 亚洲性久久影院| 成年人黄色毛片网站| 午夜福利欧美成人| 少妇的逼水好多| 嫁个100分男人电影在线观看| 无遮挡黄片免费观看| 欧美最黄视频在线播放免费| 日本 欧美在线| 最好的美女福利视频网| 淫妇啪啪啪对白视频| 我要搜黄色片| 12—13女人毛片做爰片一| 久久精品国产亚洲av天美| 俄罗斯特黄特色一大片| 亚洲一区二区三区色噜噜| 精品99又大又爽又粗少妇毛片 | 看十八女毛片水多多多| 免费在线观看成人毛片| 国产成年人精品一区二区| 日韩精品中文字幕看吧| 在线免费观看不下载黄p国产 | netflix在线观看网站| 亚洲四区av| 热99在线观看视频| av在线亚洲专区| 亚洲国产色片| 久久久久久伊人网av| 日本熟妇午夜|