為求有限元仿真自然破片的平均迎風(fēng)面積,基于蒙特卡洛剖分投影法建立了平均迎風(fēng)面積仿真計(jì)算模型,得到了破片投影點(diǎn)集圍成圖形的邊界�；诓娉朔ň帉懥擞邢拊匀黄破挠L(fēng)面積仿真計(jì)算程序,并針對(duì)所處不同象限的凹、凸多邊形頂點(diǎn)進(jìn)行了驗(yàn)證。結(jié)果表明,仿真程序?qū)Π�、凸多邊形面積求解均適用,同一多邊形所處不同象限的結(jié)果完全相同,該程序適用于求解任意不規(guī)則形狀仿真破片的平均迎風(fēng)面積。 

【文章來源】：電腦編程技巧與維護(hù). 2020,(03)

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

【部分圖文】：

爆震彈有限元破片不同時(shí)刻的破碎情況

基于蒙特卡洛剖分投影法的平均迎風(fēng)面積仿真算法流程圖如圖3所示。其中，Pe為破片單元數(shù)據(jù)庫，P為破片匯總數(shù)據(jù)庫，i2、i分為Pe庫行搜索標(biāo)記的首端和末端，q為破片序號(hào)，mtkl為蒙特卡洛次數(shù)變量，N為蒙特卡洛總次數(shù)，設(shè)置為50000次，dt為用于存放三角剖分?jǐn)?shù)據(jù)的結(jié)構(gòu)體，area＿mtkl為蒙特卡洛后的投影面積矩陣。3.2 叉乘法求多邊形面積

根據(jù)仿真程序，可得到按順時(shí)針或逆時(shí)針排列的點(diǎn)集邊界索引矩陣Ploygon。假設(shè)其點(diǎn)集坐標(biāo)按順時(shí)針方向排列，分別為P1(x1,y1）、P2(x2,y2）、…、Pn(xn,yn），由圖4可知，該多邊形面積A等于頂點(diǎn)為順時(shí)針排列的三角形面積之和減去頂點(diǎn)為逆時(shí)針排列的三角形面積之和，即圖5 叉乘法求凹、凸多邊形面積Matlab代碼

【參考文獻(xiàn)】：
期刊論文
[1]不同類型預(yù)制破片速度衰減規(guī)律及迎風(fēng)面積的修正[J]. 全嘉林,苗潤源,梁爭峰.  爆破器材. 2019(06)
[2]Delaunay三角網(wǎng)約束的Harris-SURF圖像匹配算法[J]. 張貝貝,舒紅,江萬壽.  地理與地理信息科學(xué). 2019(06)
[3]導(dǎo)氣式步槍單發(fā)精度仿真分析方法[J]. 苑大威,沙金龍,王雪皎.  彈道學(xué)報(bào). 2019(03)
[4]非金屬殼體低附帶戰(zhàn)斗部實(shí)驗(yàn)與破片飛散分析[J]. 楊世全,孫傳杰,錢立新,衛(wèi)劍鋒.  高壓物理學(xué)報(bào). 2018(04)
[5]破片式戰(zhàn)斗部的破片初始運(yùn)動(dòng)參數(shù)初始化方法研究[J]. 吳景超.  電子測(cè)試. 2017(06)
[6]應(yīng)用陰影圖的航天器迎風(fēng)面積計(jì)算方法[J]. 楊成,唐歌實(shí),李勰,陳光明.  計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào). 2015(11)
[7]閃光爆震彈沖擊波效應(yīng)模擬[J]. 秦華楊,郭三學(xué).  遼寧工程技術(shù)大學(xué)學(xué)報(bào)(自然科學(xué)版). 2014(11)



本文編號(hào)：3061748