<正>為了求得有限元結(jié)果中自然破片的平均迎風面積,根據(jù)蒙特卡洛剖分投影法建立了平均迎風面積仿真計算模型,獲得了破片投影點集圍成圖形的邊界�；谔菪畏▽τ邢拊匀黄破队斑吔绲亩噙呅蚊娣e進行了求解,并針對凹、凸多邊形的適用性和不同象限的頂點進行了驗證,獲得了有限元結(jié)果中自然破片的平均迎風面積。結(jié)果表明,基于梯形法的平均迎風面積仿真程序,對凹、凸多邊形面積求解均適用;同一多邊形所處不同象限的結(jié)果完全相同;該程序適用于求解任意不規(guī)則形狀仿真破片的平均迎風面積、質(zhì)量、數(shù)量、初速、初始位置等。在規(guī)則形狀破片平均迎風面積的計算方面,相關學者很早就開始了理論與試驗研究,并基于試驗數(shù)據(jù)形成了圓柱形、球形、長方

【文章頁數(shù)】：2 頁

【部分圖文】：

圖1 爆震彈有限元破片不同時刻的破碎情況

本文根據(jù)爆震彈的有限元數(shù)值模擬結(jié)果，獲得了1/4封裝殼體在40us時刻下的破碎情況，如圖1所示。通過LSPP后處理程序可將該時刻的單元節(jié)點數(shù)據(jù)導出，并根據(jù)Matlab編寫的同一破片搜索排序算法，可得到同一破片上的單元編號、坐標、速度、體積等，如圖2所示。圖2同一破片搜索排序后的....

圖2 同一破片搜索排序后的單元數(shù)據(jù)

圖1爆震彈有限元破片不同時刻的破碎情況2基于梯形法的平均迎風面積仿真程序

圖3 基于MCSPS方法的平均迎風面積算法流程圖

基于MCSPS法的平均迎風面積仿真算法流程圖如圖3所示。其中，Pe為破片單元數(shù)據(jù)庫，P為破片匯總數(shù)據(jù)庫，i2、i分為Pe庫行搜索標記的首端和末端，q為破片序號，mtkl為蒙特卡洛次數(shù)變量，N為蒙特卡洛總次數(shù)，本文設置為50000次，dt為用于存放三角剖分數(shù)據(jù)的結(jié)構(gòu)體，area＿m....

圖6 四頂點的凹、凸多邊形

為驗證梯形法對凹、凸多邊形的適用性，選取便于幾何計算的凹、凸四邊形進行討論，如圖6所示。根據(jù)梯形法計算，得到凸多邊形面積為6，凹多邊形面積為2.25。而通過幾何計算，凸多邊形面積為左三角形面積3與右三角形面積3之和，凹多邊形面積為兩三角形之差，與程序結(jié)果完全吻合。由此說明，梯形法....

本文編號：4037687