一种气泡帷幕对水下爆破水击波阻隔分析数值建模方法技术

本发明专利技术公开了一种气泡帷幕对水下爆破水击波阻隔分析数值建模方法，基于程序和内嵌于有限元软件

【技术实现步骤摘要】
一种气泡帷幕对水下爆破水击波阻隔分析数值建模方法


[0001]本专利技术涉及爆破工程
，具体涉及一种气泡帷幕对水下爆破水击波阻隔分析数值建模方法
。

技术介绍

[0002]水下爆破方式广泛应用于水利水电工程，如水下岩坎
、
坞门
、
围堰等的爆破拆除；为保护水下爆破作业邻近区域的水生物
、
水下设备
、
水下构筑物等，气泡帷幕得到广泛应用，其意义与作用也得到了重视
。
为评估气泡帷幕的防护效果，物理模型试验
、
现场试验和数值仿真是使用最广泛的方法
。
[0003]然而，物理模型试验和现场试验耗资巨大
、
效率低下，因而关于水下气泡帷幕的数值仿真方法得到了广大相关从业人员的青睐
。
[0004]水下气泡帷幕由放置于水体底部的发生装置产生，气泡由底部向水面移动，在这一过程中，气泡逐渐增大，呈现出“由下至上，大尺寸气泡逐渐增多”这一特征；此外，气泡的分布范围亦表现出“距离水面越近，范围越大”特征，大体呈倒梯形
。
然而，当前相关领域的现有模拟技术尚有以下不足：一
.
商业有限元软件，如
ABAQUS
，仅能建立单个气泡
、
需要重复多次建立气泡帷幕模型，不能直接构建气泡帷幕模型，需大量的手动
、
反复操作才能建立能够完全反映上述气泡帷幕空间分布特征的数值模型；二
.
当前已发表的气泡帷幕数值模拟方法未考虑“气泡距离水面越近，尺寸越大”、“气泡帷幕大体呈倒梯形”的特征
。

技术实现思路

[0005]本专利技术需要解决的技术问题是提供一种气泡帷幕对水下爆破水击波阻隔分析数值建模方法，以解决当前广泛使用的商业有限元软件不能直接建立反映随机气泡大小
、
气泡尺寸大致由下至上逐渐增大以及气泡分布范围由下至上逐渐扩大的问题
。
[0006]为解决上述技术问题，本专利技术所采取的技术方案如下
。
[0007]一种气泡帷幕对水下爆破水击波阻隔分析数值建模方法，基于程序和内嵌于有限元软件
ABAQUS
的可视化人机交互界面执行，包括以下步骤：
[0008]S1.
用户输入建模所需的所有信息；
[0009]S2.
依据用户输入的模型尺寸，创建整体模型并自动划分出炸药药包
、
水面空气层和水体；
[0010]S3.
根据用户输入的信息，遵循“随机气泡大小”、“由下至上，大尺寸气泡逐渐增多”和“气泡帷幕呈倒梯形”原则，通过逐个循环的方式创建气泡帷幕；
[0011]S4.
依据用户输入的材料参数，自动为炸药
、
空气
、
水体赋予材料属性；
[0012]S5.
依据用户输入的网格信息，自动划分网格
。
[0013]优选的，所述步骤
S3
具体包括以下步骤：
[0014]S31.
依据用户输入的帷幕倒梯形边界尺寸，随机为一个气泡指定其竖向坐标
y
；
[0015]S32.
根据竖向坐标
y
，计算气泡半径取值范围并随机生成其半径尺寸；
[0016]S33.
根据竖向坐标
y
及用户输入的缩放系数
SP
，计算气泡水平坐标
x
取值范围并随机生成其水平坐标
x
；
[0017]S34.
依据上述数据，创建一个气泡部件并放置于水平坐标
x
和竖向坐标
y
位置；
[0018]S35.
循环执行步骤
S31
至
S34
，直至创建的气泡数量达到用户指定的数量
。
[0019]优选的，所述步骤
S34
中单个气泡部件放置于水平坐标
x
和竖向坐标
y
位置与整体模型融合并保留二者的相交轮廓线
。
[0020]优选的，所述步骤
S31
中根据用户输入的帷幕倒梯形边界尺寸确定帷幕倒梯形的步骤为：
[0021]A.
根据用户指定帷幕倒梯形长边的左
、
右端点距整体模型左边界的距离，分别由
LB
和
RB
表示，从而确定倒梯形的最大宽度；
[0022]B.
根据用户指定帷幕倒梯形的长
、
短边距整体模型下边界的距离，分别由
TB
和
BB
表示，从而确定倒梯形的高度；
[0023]C.
帷幕倒梯形的短边尺寸由其长边尺寸根据用户指定的缩放系数
SP
通过比例缩放获得，短边的长度范围为
[LB+SP*(1.0
‑
BB/TB)*(RB
‑
LB),RB
‑
SP*(1.0
‑
BB/TB)*(RB
‑
LB)]。
[0024]优选的，每个所述气泡的的空间位置由水平坐标
x
和竖向坐标
y
控制，水平坐标
x
和竖向坐标
y
的取值范围为帷幕倒梯形范围；每个气泡的
z
向坐标相同，为一固定值；每个所述气泡的竖向坐标
y
服从
[TB
，
BB]上的均匀分布，每个气泡的水平坐标
x
也服从均匀分布，并根据其竖向坐标
y、
缩放系数
SP
以及用户指定的
LB
和
RB
随机产生，随机取值范围为
[LB+SP*(1.0
‑
y/TB)*(RB
‑
LB),RB
‑
SP*(1.0
‑
y/TB)*(RB
‑
LB)]。
[0025]优选的，所述气泡为圆形且其半径服从均匀分布，气泡半径的上
、
下限，分别由
LR
和
HR
表示，气泡半径的随机取值范围为
[LR+LR*(1.0
‑
y/TB),HR
‑
HR*(1.0
‑
y/TB)]，该取值范围表示竖向坐标
y
越大，气泡半径的取值范围越大
。
[0026]优选的，所述水体
、
炸药和空气均使用
Euler
网格并属同一个
Euler
域且可在该
Euler
域内流动而网格不移动
。
[00本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.
一种气泡帷幕对水下爆破水击波阻隔分析数值建模方法，其特征在于：基于程序和内嵌于有限元软件
ABAQUS
的可视化人机交互界面执行，包括以下步骤：
S1.
用户输入建模所需的所有信息；
S2.
依据用户输入的模型尺寸，创建整体模型并自动划分出炸药药包
、
水面空气层和水体；
S3.
根据用户输入的信息，遵循“随机气泡大小”、“由下至上，大尺寸气泡逐渐增多”和“气泡帷幕呈倒梯形”原则，通过逐个循环的方式创建气泡帷幕；
S4.
依据用户输入的材料参数，自动为炸药
、
空气
、
水体赋予材料属性；
S5.
依据用户输入的网格信息，自动划分网格
。2.
根据权利要求1所述的一种气泡帷幕对水下爆破水击波阻隔分析数值建模方法，其特征在于：所述步骤
S3
具体包括以下步骤：
S31.
依据用户输入的帷幕倒梯形边界尺寸，随机为一个气泡指定其竖向坐标
y
；
S32.
根据竖向坐标
y
，计算气泡半径取值范围并随机生成其半径尺寸；
S33.
根据竖向坐标
y
及用户输入的缩放系数
SP
，计算气泡水平坐标
x
取值范围并随机生成其水平坐标
x
；
S34.
依据上述数据，创建一个气泡部件并放置于水平坐标
x
和竖向坐标
y
位置；
S35.
循环执行步骤
S31
至
S34
，直至创建的气泡数量达到用户指定的数量
。3.
根据权利要求2所述的一种气泡帷幕对水下爆破水击波阻隔分析数值建模方法，其特征在于：所述步骤
S34
中单个气泡部件放置于水平坐标
x
和竖向坐标
y
位置与整体模型融合并保留二者的相交轮廓线
。4.
根据权利要求2所述的一种气泡帷幕对水下爆破水击波阻隔分析数值建模方法，其特征在于：所述步骤
S31
中根据用户输入的帷幕倒梯形边界尺寸确定帷幕倒梯形的步骤为：
A.
根据用户指定帷幕倒梯形长边的左
、
右端点距整体模型左边界的距离，分别由
LB
和
RB
表示，从而确定倒梯形的最大宽度；
B.
根据用户指定帷幕倒梯形的长
、
短边距整体模型下边界的距离，分别由
TB
和
BB

【专利技术属性】
技术研发人员：姚二雷，朱莅，黎卫超，饶宇，杨招伟，马晨阳，吴新霞，赵根，刘志芳，王秀杰，袁俊，张继楷，
申请(专利权)人：长江水利委员会长江科学院，
类型：发明
国别省市：