(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202210610507.3 (22)申请日 2022.05.31 (71)申请人 西安近代化学研究所 地址 710065 陕西省西安市雁塔区丈 八东 路168号 (72)发明人 韩璐　李芝绒　蒋海燕　张玉磊　 林钦栋　 (74)专利代理 机构 西安恒泰知识产权代理事务 所 61216 专利代理师 王孝明 (51)Int.Cl. G06F 30/15(2020.01) G06F 30/20(2020.01) G06F 119/14(2020.01) (54)发明名称 一种考虑舱壁遮挡的舰船舱室内爆炸冲击 波载荷快速预测方法 (57)摘要 本发明提供了一种考虑舱壁遮挡的舰船舱 室内爆炸冲击波载荷快速预测方法， 该方法包括 以下步骤： 步骤S1， 三维舰船舱室结构模型离散 化及等效药量转换； 步骤S2， 划分计算步长； 步骤 S3， 计算各舱壁压力与破坏状态； 步骤S4， 剔除实 际未破坏的舱壁； 步骤S5， 计算破坏舱壁的压力 载荷折减； 本发 明通过按步长求解的方式进行舱 壁载荷和毁伤状态计算， 对于大尺寸舰船目标， 可大大提高其计算效率， 减少计算时长。 本发明 通过构建炸点与舱壁几何中心射线的方式， 基于 循环扫描的方法判断在被破坏舱壁前是否存在 未破坏壁面的遮挡情况， 从而实现了将实际不会 发生破坏的舱 壁剔除， 使计算结果更符合实际。 权利要求书3页 说明书7页 附图3页 CN 115081110 A 2022.09.20 CN 115081110 A 1.一种考虑舱壁遮挡的舰船舱室内爆炸冲击波载荷快速预测方法， 其特征在于， 该方 法包括以下步骤： 步骤S1， 三维舰船舱室结构模型离 散化及等效药量 转换： 首先， 将三维舰船舱室实体结构进行退化， 采用离散化薄壳单元对各舱壁结构进行模 型替代， 三维舰船舱室结构模型 的尺寸与舰船目标 的实际尺寸保持一致， 对实际结构 中各 舱壁相连接的位置采用共线处 理， 以保证舱室结构的封闭； 其次， 炸药的等效TNT装药质量 We可按式Ⅰ进行换算： 式中： W为实际炸药的装药质量； Qe为实际炸药的爆热； Q为TNT炸药的爆热； 步骤S2， 划分计算 步长： 统计所有舱壁到炸点(x0,y0,z0)的距离， 依据各舱壁所包含的所有面元节点坐标， 按式 Ⅱ进行其几何中心(xci,yci,zci)的求解计算； 每个舱壁的几何中心至 炸点的距离dci按式Ⅲ进行计算： 依据所有的舱壁与炸点的距离dci， 统计出其中最小值dci_min和最大值dci_max， 依据 (dci_max‑dci_min‑d0)确定计算 步骤数量 n， 并以此确定计算 步长： Δ(dci_max‑dci_min‑d0)/n； 式中： xi， yi， zi分别为该舱 壁包含的所有离 散化细分面元节点的三维坐标； i为离散化细分面元节点编号； n为步骤数量； dc为两点之间的距离； d0为所选择的初始计算距离； 步骤S3， 计算各舱 壁压力与破坏状态： 按d0+Δ(dci_max‑dci_min‑d0)/n×i(i＝0,1,2 …n)依次计算步长范 围之内所有舱室壁面 各点的反射压力， 按式 Ⅳ进行计算： Δpref＝Δpref·cos2αV·cos2αH  式Ⅳ； 式中： αV和 αH分别为“炸点‑测点”连线与舱 壁两向垂面的夹角； Δpref为炸点在该舱 壁平面投影处的冲击波正反射压力； Δpref按式Ⅴ进行计算；权　利　要　求　书 1/3 页 2 CN 115081110 A 2式中： We为由式Ⅰ计算的炸药的等效TNT装药质量， 单位 为： kg； rcic为炸点距离舱 壁最近点的距离， 单位 为： m， 由空间点到平面 最短距离计算得到； p0为环境压力， 取0.101MPa； 按d0+Δ(dci_max‑dci_min‑d0)/n×i(i＝0,1,2 …n)依次计算步长范 围之内所有舱室壁面 的破坏情况， 判断依据为rcic小于造成舱壁结构破坏的临界爆距dcic， dcic可由式Ⅵ计算获 得： 式中： h为舱壁的等效厚度， 单位 为： m； ρ 为舱壁结构材 料密度， 单位 为： kg/m3； σd为舱壁结构材 料的动态屈服强度， 单位 为： Pa； εf为舱壁结构材 料的断裂 应变； 步骤S4， 剔除实际未破坏的舱 壁： 统计所有被破坏的舱壁结构， 依次构建由炸点(x0,y0,z0)到每个被破坏舱壁几何中心 [(xd1,yd1,zd1),(xd2,yd3,zd4),…,(xdi,ydi,zdi)]的空间射线： 式中： x， y， z分别为空间构造直线方程的变量； 将该射线与目标所有舱室壁面进行求交， 并统计该射线路径上所有穿过的舱壁， 统计 方法为首 先计算空间直线与舱 壁平面的交点， 其次确认该交点在舱 壁平面内； 空间一点P不在凸多边形平面A(n1,n2,...,nk)内的充要条件为式 Ⅷ： 其中： 为三角形； n1,n2,...,nk为凸多边形平面A的第1,2,. ..,k个顶点； i为顶点的序号； 统计射线与该舱壁平面交汇前的穿透路径， 判断路径上共有多少个其它舱壁， 同时判 断是否存在未毁伤的壁面， 如果有则说明当前壁面将受到前述未发生破坏壁面的遮挡作 用， 从而将该舱 壁从破坏壁 面集合中去除； 步骤S5， 计算破坏舱 壁的压力载荷折减： 由于实际破坏舱室壁面将损耗一部分爆炸冲击波能量， 因此在统计完所有受到毁伤的权　利　要　求　书 2/3 页 3 CN 115081110 A 3