[发明专利]一种爆破块度空间全方位分区耦合数值仿真方法

权利要求书

1.一种爆破块度空间全方位分区耦合数值仿真方法，其特征在于，包括：

步骤1，根据开挖爆破对象的几何尺寸，利用LS-DYNA工具建立爆破对象的连续介质数值计算模型，记为LS-DYNA模型；

步骤2，在LS-DYNA模型中导入小粒径爆破块度预测模型，利用动力有限元方法预测连续介质的小粒径爆破块度的分布曲线，并确定小粒径和中大粒径爆破块度的分布边界，提取分布边界的质点振动速度和应力的时程曲线；

步骤3，利用DDA工具构建中大粒径的爆破块度预测模型，记为DDA模型，以步骤2确定的分布边界作为DDA模型的初始空腔边界，DDA模型的参数直接使用LS-DYNA模型的参数值；

步骤4，将步骤2提取的质点振动速度和应力施加于DDA模型的初始空腔边界，进行非连续计算，获得中大粒径爆破块度的分布特性，从而获得中大粒径爆破块度分布曲线；

步骤5，合并小粒径爆破块度和中大粒径爆破块度的分布曲线，得到爆破块度空间全方位耦合分布曲线图。

2.如权利要求1所述的爆破块度空间全方位分区耦合数值仿真方法，其特征在于：

小粒径和中大粒径爆破块度的分布边界为10mm爆破块度的分布边界，所述小粒径爆破块度指粒径不大于10mm的爆破块度；所述中大粒径的爆破块度指粒径大于10mm的爆破块度。

3.如权利要求1所述的爆破块度空间全方位分区耦合数值仿真方法，其特征在于：

步骤2具体为：

基于LS-DYNA的二次开发模块，将小粒径爆破块度预测模型导入，采用历史变量记录爆破块度值，预测连接介质的小粒径爆破块度的分布图；采用历史变量记录的爆破块度值，确定10mm爆破块度的分布边界，对预测结果云图进行筛选，并提取分布边界上的质点振动速度和应力。

4.如权利要求1所述的爆破块度空间全方位分区耦合数值仿真方法，其特征在于：

步骤3中，所述构建中大粒径的爆破块度预测模型，即将DDA模型中炮孔至10mm爆破块度的分布边界之间的区域删除。