基于叉乘法的有限元自然破片平均迎风面积编程与实现

发布时间：2021-03-03 18:45

　　为求有限元仿真自然破片的平均迎风面积,基于蒙特卡洛剖分投影法建立了平均迎风面积仿真计算模型,得到了破片投影点集围成图形的边界。基于叉乘法编写了有限元自然破片的迎风面积仿真计算程序,并针对所处不同象限的凹、凸多边形顶点进行了验证。结果表明,仿真程序对凹、凸多边形面积求解均适用,同一多边形所处不同象限的结果完全相同,该程序适用于求解任意不规则形状仿真破片的平均迎风面积。 

【文章来源】：电脑编程技巧与维护. 2020,(03)

【文章页数】：3 页

【部分图文】：

爆震弹有限元破片不同时刻的破碎情况

基于蒙特卡洛剖分投影法的平均迎风面积仿真算法流程图如图3所示。其中，Pe为破片单元数据库，P为破片汇总数据库，i2、i分为Pe库行搜索标记的首端和末端，q为破片序号，mtkl为蒙特卡洛次数变量，N为蒙特卡洛总次数，设置为50000次，dt为用于存放三角剖分数据的结构体，area＿mtkl为蒙特卡洛后的投影面积矩阵。3.2 叉乘法求多边形面积

根据仿真程序，可得到按顺时针或逆时针排列的点集边界索引矩阵Ploygon。假设其点集坐标按顺时针方向排列，分别为P1(x1,y1）、P2(x2,y2）、…、Pn(xn,yn），由图4可知，该多边形面积A等于顶点为顺时针排列的三角形面积之和减去顶点为逆时针排列的三角形面积之和，即图5 叉乘法求凹、凸多边形面积Matlab代码

【参考文献】：
期刊论文
[1]不同类型预制破片速度衰减规律及迎风面积的修正[J]. 全嘉林,苗润源,梁争峰.  爆破器材. 2019(06)
[2]Delaunay三角网约束的Harris-SURF图像匹配算法[J]. 张贝贝,舒红,江万寿.  地理与地理信息科学. 2019(06)
[3]导气式步枪单发精度仿真分析方法[J]. 苑大威,沙金龙,王雪皎.  弹道学报. 2019(03)
[4]非金属壳体低附带战斗部实验与破片飞散分析[J]. 杨世全,孙传杰,钱立新,卫剑锋.  高压物理学报. 2018(04)
[5]破片式战斗部的破片初始运动参数初始化方法研究[J]. 吴景超.  电子测试. 2017(06)
[6]应用阴影图的航天器迎风面积计算方法[J]. 杨成,唐歌实,李勰,陈光明.  计算机辅助设计与图形学学报. 2015(11)
[7]闪光爆震弹冲击波效应模拟[J]. 秦华杨,郭三学.  辽宁工程技术大学学报(自然科学版). 2014(11)



本文编号：3061748