一种水下爆破减震孔布设模拟方法及装置的制造方法

一种水下爆破减震孔布设模拟方法及装置的制造方法
【专利摘要】本发明提供的水下爆破减震孔布设模拟方法及装置，属于水下爆破工程领域。该水下爆破减震孔布设模拟方法及装置在通过炮孔及减震孔的初始布设参数、材料模型及状态方程建立有限元分析模型并计算出相应的模拟结果后，每一次更改减震孔初始布设参数中的一个，获得新的有限元分析模型，通过比较不同的模拟结果，得出所更改的减震孔布设参数的较优值，对提高水下爆破施工效率、降低施工成本、确保紧邻建(构)筑物安全使用具有重要的促进意义，对进一步补充、丰富水下爆破基本理论，具有积极的推进作用。
【专利说明】
-种水下爆破减震孔布设模拟方法及装置
技术领域
[0001] 本发明设及水下爆破工程领域，具体而言，设及一种水下爆破减震孔布设模拟方 法及装置。
【背景技术】
[0002] 在紧邻建(构)筑物环境下进行水下爆破，确定合理的减震方案是控制技术的核屯、 工作，目前国内的研究成果或工程实践均仅停留在水下爆破减震孔的设计方法，而对减震 孔的不同布设参数下的减震机理及减震效果进行数值仿真模拟研究能够有效的完善水下 爆破基本理论，是目前亟待研究的课题。

【发明内容】

[0003] 有鉴于此，本发明的目的在于提供一种水下爆破减震孔布设模拟方法及装置。
[0004] 第一方面，本发明实施例提供的一种水下爆破减震孔布设模拟方法，包括:获得建 立实体模型所需的初始布设参数和所述实体模型包括的材料属性；
[0005] 根据所述材料属性获得对应的材料模型和状态方程;根据所述材料模型和所述状 态方程W及所述初始布设参数建立实体模型;对所述实体模型进行网格划分，生成对应的 有限元分析模型；
[0006] 采用连续介质力学运动描述方法对所述有限元分析模型进行计算，得到水下爆破 减震孔布设的模拟结果，所述模拟结果包括减震保护对象所受到的等效应力模拟数据、水 击波压力模拟数据W及位移距离模拟数据。
[0007] 优选地，所述初始布设参数包括炮孔初始布设参数和减震孔初始布设参数，所述 减震孔初始布设参数包括减震孔距炮孔的距离、减震孔填充介质W及减震孔的直径，所述 方法还包括:更改所述减震孔初始布设参数中的一个，得到新的有限元分析模型；
[000引对新得到的有限元分析模型采用同一连续介质力学运动描述方法进行计算，得到 新的模拟结果;根据新的模拟结果，确定所更改参数的优选值。
[0009] 优选地，所述方法在对所述有限元分析模型进行计算之前还包括：设定所述有限 元分析模型的边界条件为无反射边界。
[0010] 优选地，所述方法在对所述有限元分析模型进行计算之前还包括:对计算所述有 限元分析模型的分析选项进行设置，所述分析选项包括求解时间和时间步，所述求解时间 为对所述有限元分析模型预设的计算时间，所述时间步为每一次对所述实体模型进行初始 化求解使其计算收敛。
[0011] 优选地，所述方法在所述对所述实体模型进行网格划分，生成对应的有限元分析 模型的步骤之后还包括:检测生成的有限元分析模型中网格单元是否满足预设的计算精度 的要求;若否，则调整网格大小，对所述实体模型重新进行网格划分。
[0012] 第二方面，本发明另一实施例提供的一种水下爆破减震孔布设模拟装置，包括:第 一获得模块，用于获得建立实体模型所需的初始布设参数和所述实体模型包括的材料属 性；
[0013] 第二获得模块，用于根据所述材料属性获得对应的材料模型和状态方程；
[0014] 实体模型建立模块，用于根据所述材料模型、所述状态方程W及所述初始布设参 数建立实体模型；
[0015] 网格划分模块，用于对所述实体模型进行网格划分，生成对应的有限元分析模型；
[0016] 计算模块，用于采用连续介质力学运动描述方法对所述有限元分析模型进行计 算，得到水下爆破减震孔布设的模拟结果，所述模拟结果包括减震保护对象所受到的等效 应力模拟数据、水击波压力模拟数据W及位移距离模拟数据。
[0017] 优选地，所述初始布设参数包括炮孔初始布设参数和减震孔初始布设参数，所述 减震孔初始布设参数包括减震孔距炮孔的距离、减震孔填充介质W及减震孔的直径，所述 装置还包括:参数更改模块，用于更改所述减震孔初始布设参数中的一个，得到新的有限元 分析模型;所述计算模块，还用于对新得到的有限元分析模型采用同一连续介质力学运动 描述方法进行计算，得到新的模拟结果;参数优选值确定模块，用于根据新的模拟结果，确 定所更改参数的优选值。
[0018] 优选地，所述装置还包括边界条件设定模块，用于设定所述有限元分析模型的边 界条件为无反射边界。
[0019] 优选地，所述装置还包括:分析选项设置模块，用于对计算所述有限元分析模型的 分析选项进行设置，所述分析选项包括求解时间和时间步，所述求解时间为对所述有限元 分析模型预设的计算时间，所述时间步为每一次对所述实体模型进行初始化求解使其计算 收敛。
[0020] 优选地，所述装置还包括:检测模块，用于检测生成的有限元分析模型中网格单元 是否满足预设的计算精度的要求;所述网格划分模块，还用于在划分的网格单元不满则预 设的计算精度时，对所述实体模型重新进行网格划分。
[0021] 与现有技术相比，本发明提供的水下爆破减震孔布设模拟方法及装置，在通过炮 孔及减震孔的初始布设参数、材料模型及状态方程建立有限元分析模型并计算出相应的模 拟结果后，每一次更改减震孔初始布设参数中的一个，获得新的有限元分析模型，通过比较 不同的模拟结果，得出所更改的减震孔布设参数的较优值，对提高水下爆破施工效率、降低 施工成本、确保紧邻建(构)筑物安全使用具有重要的促进意义，对进一步补充、丰富水下爆 破基本理论，具有积极的推进作用。
[0022] 为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合 所附附图，作详细说明如下。
【附图说明】
[0023] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所 需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施 例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可W根据运些附图获 得其他的附图。通过附图所示，本发明的上述及其它目的、特征和优势将更加清晰。在全部 附图中相同的附图标记指示相同的部分。并未刻意按实际尺寸等比例缩放绘制附图，重点 在于示出本发明的主旨。
[0024] 图I为本发明较佳实施例提供的应用水下爆破减震孔布设模拟装置的本地终端的 结构框图；
[0025] 图2为本发明第一实施例提供的一种水下爆破减震孔布设模拟方法的流程图；
[0026] 图3为本发明第二实施例提供的另一种水下爆破减震孔布设模拟方法的流程图；
[0027] 图4为本发明第=实施例提供的另一种水下爆破减震孔布设模拟方法的流程图；
[0028] 图5为本发明第四实施例提供的一种水下爆破减震孔布设模拟装置的功能模块框 图。
[00巧]主要元件符号说明 [0030]
[0031]

【具体实施方式】
[0032] 下面将结合本发明实施例中附图，对本发明实施例中的技术方案进行清楚、完整 地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。通常在 此处附图中描述和示出的本发明实施例的组件可WW各种不同的配置来布置和设计。因 此，W下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的 范围，而是仅仅表示本发明的选定实施例。基于本发明的实施例，本领域技术人员在没有做 出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0033] 应注意到:相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一 个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。同时，在本发明的 描述中，术语"第一"、"第二"等仅用于区分描述，而不能理解为指示或暗示相对重要性。
[0034] 如图1所示，是本发明较佳实施例提供的本地终端100的方框示意图，所述本地终 端100包括水下爆破减震孔布设模拟装置101、存储器102、存储控制器103W及处理器104。 其中，存储器102、存储控制器103W及处理器104各元件相互之间直接或间接地电性连接， W实现数据的传输或交互。例如，运些元件相互之间可通过一条或多条通讯总线或信号线 实现电性连接。所述水下爆破减震孔布设模拟装置101包括至少一个可W软件或固件 (firmware)的形式存储于所述存储器102中的软件功能模块。所述处理器104用于执行存储 器102中存储的可执行模块，例如水下爆破减震孔布设模拟装置101包括的软件功能模块或 计算机程序。
[0035] 存储器102可包括高速随机存储器，还可包括非易失性存储器，如一个或者多个磁 性存储装置、闪存、或者其他非易失性固态存储器。处理器104W及其他可能的组件对存储 器102的访问可在存储控制器103的控制下进行。
[0036] 处理器104可能是一种集成电路忍片，具有信号的处理能力。上述的处理器可W是 通用处理器，包括中央处理器、网络处理器等;还可W是数字信号处理器、专用集成电路、现 成可编程口阵列或者其他可编程逻辑器件、分立口或者晶体管逻辑器件、分立硬件组件，可 W实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可W是微 处理器或者该处理器104也可W是任何常规的处理器等。
[0037] 请参阅图2,是本发明第一实施例提供的水下爆破减震孔布设模拟方法流程图，所 应说明的是，本发明所述的水下爆破减震孔布设模拟方法并不W图2W及W下所述的具体 顺序为限制。下面将对图2所示的具体流程进行详细阐述。
[0038] 步骤S110,获得建立实体模型所需的初始布设参数和所述实体模型包括的材料属 性。
[0039] 本较佳实施例中，初始布设参数包括炮孔的初始布设参数和减震孔的初始布设参 数，其中炮孔的初始布设参数包括但不限于，炮孔直径、炮孔间距、炸药装药直径及装药高 度等;减震孔的初始布设参数包括但不限于，减震孔设置的组数及组距、每一组减震孔设置 的排数及排距、减震孔距爆源的距离、减震孔填充介质W及减震孔的尺寸等;实体模型所 包括材料属性由具体模拟的工程所包括的材料决定。
[0040] 步骤Slll，根据所述材料属性获得对应的材料模型和状态方程。
[0041] 步骤S112,根据所述材料模型、所述状态方程W及所述初始布设参数建立实体模 型。
[0042] 步骤Sl 13，对所述实体模型进行网格划分，生成对应的有限元分析模型。
[0043] 定义单元类型后，应用该单元类型对建立的实体模型进行网格划分，生成有限元 分析模型。现有的网格划分方法有多种，如自由网格划分、映射网格划分及混合网格划分 等，本实施例中优选地，采用映射网格划分方法对建立的实体模型进行划分，因为映射网格 划分方法划分后生成的有限元分析模型中，网格单元更加的规整，利于获得较优的模拟结 果，提高模拟结果的计算精度，且可W根据操作者的意愿生成确定数目的网格单元，加快后 期的计算速度，另外，映射网格划分方法可W避免产生一些特别崎形的网格单元，运也是其 主要优势之一。
[0044] 步骤S114,采用连续介质力学运动描述方法对所述有限元分析模型进行计算，得 到水下爆破减震孔布设的模拟结果，所述模拟结果包括减震保护对象所受到的等效应力模 拟数据、水击波压力模拟数据W及位移距离模拟数据。
[0045] 连续介质力学的运动描述方法包括=个经典算法:拉格朗日算法、欧拉算法W及 任意拉格朗日欧拉算法，运=个算法均可W对水下爆破减震孔布设模拟的有限元分析模型 进行计算得到相应的模拟结果。下面对=个算法做简要的介绍：
[0046] 拉格朗日算法中，炸药单元与结构单元之间可W通过共用节点方式建立联系，具 体方法是将炸药与结构采用8节点实体单元模型，炸药与结构单元之间有相同的节点。也可 W通过接触定义二者之间的荷载传递，与共节点方法不同之处在于，炸药单元与被爆炸结 构单元之间网格是独立的，通过接触来定义彼此之间的关系，具体的接触方式有侵彻接触、 滑动接触与面面接触。通过定义接触模拟爆炸的主要问题在于，计算时会由于接触的计算 而导致不收敛，计算过程突然中断而不进行任何的错误及警告提示。
[0047] 欧拉算法中，网格被固定在空间是不变形的，物质通过网格流进流出，物质的大变 形不影响时间步长的计算。欧拉算法不需要定义各部分材料之间的复杂接触关系，且计算 稳定性好。其特点是，网格数目巨大、求解时间长、网格中物质流动界面不清晰，但是计算稳 定性好。
[004引对于ALEUrbitrary Lagrange-Euler,任意拉格朗日欧拉）算法，其特点是它采用 的网格是每一步或每隔若干步，根据物质区域的边界构造一个合适的网格，W避免在严重 扭曲的网格上进行计算。该算法提供了输运和禪合算法来保证得到一个准确、稳定、守恒和 收敛的结果。
[0049] 本实施例中，通过比较=种算法的优缺点，优选任意拉格朗日欧拉算法对水下爆 破减震孔布设模拟的有限元分析模型进行计算。
[0050] 下面，本实施例将W具体的科研项目"重庆万州长江公路大桥防撞设施水下基坑 及基槽开挖"为工程依托，采用显式动力分析有限元软件ANSYS/LS-DYNA，对具体的模拟方 法做进一步的阐述，如下所示：
[0051] 计算模型在考虑实际情况的基础上，结合数值模拟技术的特点，将岩体设置在水 下30米处，岩层厚度8米，水上空气层10米;炮孔深度6米，炮孔底部预留2米岩体，炮孔直径 设为100毫米，炸药的装药直径为70毫米，装药高度5米，炮孔与炸药间采用水禪合，炮孔间 距2.5米;减震孔设置两组，组距5.5米，每一组设置两排，排距为1.5米，每一排设置5个减震 孔，孔距300毫米，减震孔距离爆源和长江大桥桥墳的距离分别15米和46米，减震孔采用水 填充，减震孔的直径为120毫米。
[0052] 炸药采用TNT材料模型，并WJWL状态方程描述炸药爆轰产生的压力，压力与体积 变形之间的单"累味
[0化3]
[0054] 式中：P为压力，4、8、扣、1?2和《为1胖1状态方程的五个参数，￥为相对体积，6为单位 体积内能。
[0055] 炸药选用防水性能好、爆能较高、运输使用安全系数较大的2#岩石乳化炸药，该炸 药的材料参数及JWL状态方程参数如表1所示：
[0化6]
[0057]表 1
[005引水对水在冲击波下的压力、密度关系采用Mie-Gruneisen状态方程描述：
[0化9]
[0060]式中：P为压力，E为单位体积的比内能，PO为介质的初始密度，y=l/(l-V)，V为相 对体积，C为介质中声速，丫日为Gruneisen参数，Si、S2、S3为冲击波波速-波后质点粒子速度 曲线斜率的系数，a为Mie-Gruneisen系数修正项。
[0061 ]水的材料参数及状态方程参数如表2所示：
[0063] 表 2
[0064] 其中，C为冲击波波速-波后质点粒子速度曲线的截距。
[0(?日]空气材料状态方程如下所示：
[0066] P = Co+Ciii+C2ii^+C3ii^+ (C4+C5ii+C6ii^) E
[0067] 式中：Co~Cs为状态方程系数，y为比体积，空气材料的参数及状态方程参数如表3 所示：
[006引
[0069] 表 3
[0070] 本实施例中，岩体采用塑性随动模型，目化S-DYNA提供的MAT_PLASTIC_KI肥MATIC 模型，此模型不仅考虑了岩石介质的弹塑性特性，而其能够对材料的强化效应(随动强化与 各向同性强化）和应变率变化效应加 W描述，同时还带有失效应变。应变率用Cowper- Symonds横巧夹考虑，屈服应力Ov与应巧率e之间的关系如下所示：
[0071]
[0072] 式中：O日为岩体的初始屈服应力，e为加载应变率，C、P为Cowper-Symonds应变率参 数，Eo为杨氏模量，Etan为切线模量，峽J各向同性硬化和随动硬化贡献的硬化参数且 1，巧为有效塑性应变。
[0073] 由于水下爆破岩体为中风化砂岩，经室内物理学实验，得到中风化砂岩岩体物理 力学参数如表4所示： 「007/11
[00巧]表4
[0076]工程重点减震保护对象为万州公路长江大桥，其材料属性为钢筋混凝±，材料的 物理力学参数如表5所示：
[0077]
[0078] 表 5
[0079] 对于表5中的参数需要说明的是，由于钢筋混凝±中钢筋含量较少，因此，取混凝 ±的弹性模量代替钢筋混凝±的弹性模量，取混凝±的剪切模量代替钢筋混凝±的剪切模 量。
[0080] 基于上述的初始布设参数、材料模型W及状态方程，建立针对上述工程的减震孔 布设模拟的四孔水下爆破=维实体模型，并利用Solid64单元类型对实体模型进行映射网 格划分，得到相应的有限元分析模型，最后采用ALE算法进行计算，得到减震模拟结果，具体 地，当时间t = 1999.9US时，长江大桥桥墳处的冯?米塞斯应力，即Mises应力达到最大值 2.2144e^MPa，其中MPa表示压强的单位兆帕，当时间达至Ijt= 19000iis后，桥墳处的Mises应 力逐渐变小趋于常数，其中iis表示时间单位微秒；当时间t = 1999.9iis时，桥墳处受到的水 击波压力达到最大值，为0.523MPa，在t = 1999.9iis~4499.4iis时间段内，桥墳处所受的水 击波压力值相对其他时间段有一定程度的波动，在该时间段后，桥墳处所受的水击波压力 值逐渐减小，可推测该时间段是爆破冲击波传播到桥墳处的时间段，在t = 1900化S后，桥墳 处所受水击波压力逐渐减小并趋于常数，表明在该时刻后，桥墳处不再受到因为爆破而产 生的水击波压力；桥墳处的位移模拟数据为，在t = 1499.7iis~1999.9iis时间段，桥墳受到 爆破冲击波的作用并发生微小位移，并在t = 1999.9iis时刻位移达到最大值1.431e^cm，cm 表示距离单位厘米。
[0081] 为了更好的分析减震模拟结果，优选地，还可W在桥墳处自上至下取=个特征单 元，作为当前活动单元，W分析它们的压力和位移随时间变化特征，运里不再寶述。
[0082] 请参阅图3,是本发明第二实施例提供的水下爆破减震孔布设模拟方法流程图，所 应说明的是，本发明所述的水下爆破减震孔布设模拟方法并不W图3W及W下所述的具体 顺序为限制。下面将对图3所示的具体流程进行详细阐述。
[0083] 步骤S210,获得建立实体模型所需的初始布设参数和所述实体模型包括的材料属 性。
[0084] 步骤S211，根据所述材料属性获得对应的材料模型和状态方程。
[0085] 步骤S212,根据所述材料模型、所述状态方程W及所述初始布设参数建立实体模 型。
[0086] 步骤S213,对所述实体模型进行网格划分，生成对应的有限元分析模型。
[0087] 步骤S214,检测生成的有限元分析模型中网格单元是否满足预设的计算精度的要 求，若否，则返回步骤S213,即对所述实体模型重新进行网格划分。
[0088] 本实施例中，采用映射网格划分方法对所建立的实体模型进行网格划分后并得到 相应的有限元分析模型后，需要判断生成的有限元分析模型中，网格单元数目、单元类型等 是否满足预设的计算精度要求，比如，在单个网格单元所预设的求解时间相同的条件下，网 格单元数目越多求解时间越长，计算精度也越高，反之，网格单元数目越少求解时间越短， 但是计算精度也越差。一般，在网格单元数目达到某一数值或者某一范围内后，再继续增加 网格单元的数目，计算精度便不会出现明显提升，在真实的模拟过程中，可多次试验得到较 优的网格单元数目值，W使得在保证计算精度的前提下，尽可能的缩短计算时间。
[0089] 步骤S215,设定所述有限元分析模型的边界条件为无反射边界。
[0090] 无反射边界又称透射边界或是无反应边界，主要应用于无限体或半无限体中，为 了减小研究对象的尺寸而采用的边界条件。无反射边界根据虚功原理将边界上的分布阻尼 转化成等效节点力加到边界上，即列出所有组成无反射边界的单元，在所有无反射边界的 单元上加上粘性正应力和剪应力。
[0091] 步骤S216,对计算所述有限元分析模型的分析选项进行设置，所述分析选项包括 求解时间和时间步，所述求解时间为对所述有限元分析模型预设的计算时间，所述时间步 为每一次对所述实体模型进行初始化求解使其计算收敛。
[0092] 步骤S217,采用连续介质力学运动描述方法对所述有限元分析模型进行计算，得 到水下爆破减震孔布设的模拟结果，所述模拟结果包括减震保护对象所受到的等效应力 模拟数据、水击波压力模拟数据W及位移距离模拟数据。
[0093] 本较佳实施例提供的水下爆破减震孔布设模拟方法同样可W科研项目"重庆万州 长江公路大桥防撞设施水下基坑及基槽开挖"为具体工程依托，采用显式动力分析有限元 软件ANSYS/LS-DYNA，对该工程中减震孔的布设参数做具体的模拟，W得到相应的模拟结 果，本实施例在第一实施例提供的方法的基础上，在对有限元分析模型进行计算之间，还包 括了设定无反射边界条件的步骤、检测网格单元是否满足预设计算精度要求的步骤W及分 析选项设置步骤，本实施例中未提到的具体操作过程可参见第一实施例的描述，运里不再 寶述。
[0094] 请参阅图4,是本发明第=实施例提供的水下爆破减震孔布设模拟方法的流程图， 其与第一实施例和第二实施例不同的是，本实施例提供的水下爆破减震孔布设模拟方法还 包括：
[00M]步骤S310,该步骤在步骤S114或步骤S217后执行，该步骤S310为更改减震孔初始 布设参数中的一个，得到新的有限元分析模型。
[0096] 减震孔初始布设参数包括多个，本实施例中则主要是W减震孔距炮孔的距离、减 震孔填充介质W及减震孔的直径=个参数为例做相应的阐述。在第一实施例或第二实施例 所提供的水下爆破减震孔布设模拟方法的基础上，每一次只更改减震孔距炮孔的距离、减 震孔填充介质或减震孔的直径中的一个参数，而保持其余所有参数不变，建立对应的新的 实体模型，同样地对新的实体模型进行网格划分得到新的有限元分析模型。
[0097] 步骤S311，对新得到的有限元分析模型采用同一连续介质力学运动描述方法进行 计算，得到新的模拟结果。
[0098] W及步骤S312,根据新的模拟结果，确定所更改参数的优选值。
[0099] 通过计算更改上述的减震孔距炮孔的距离、减震孔填充介质或减震孔的直径中的 一个参数所得到的新的有限元分析模型，得到新的模拟结果，对比当前的模拟结果与之前 的有限元分析模型的模拟结果，确定设置哪一个参数能够使减震效果更佳，从而使减震保 护对象避免受到损害。
[0100] 下面将在第一实施例中所述的W科研项目"重庆万州长江公路大桥防撞设施水下 基坑及基槽开挖"为工程依托的【具体实施方式】的基础上，对本实施例提供的水下爆破减震 孔布设模拟方法做详细的阐述。
[0101] 首先，在保持其余初始布设参数不变的条件下，依次建立减震孔距爆源的距离为 30米、40米W及49米的计算模型，并采用ALE算法对新得到的=个新的有限元分析模型进行 计算，得到对应的模拟结果，具体地，结合第一实施例中的结果可得：
[0102]
[010；3]表 6
[0104] 如表6所示，在水下爆破中，减震孔距爆源距离的不同影响了爆击波的传递，从而 改变了桥墳所受到的Mises应力、水击波压力W及位移。当减震孔距爆源15米时，桥墳处所 受到的Mises应力最大值、水击波压力最大值W及位移最大值分别为2.2144e^MPa、 3.315MPa、1.431 e-3cm。当时间达至Ij t = 19000iis后，桥梁处的位移逐渐稳定，表明水下爆破使 得桥墳发生的微小位移；当减震孔位置距离爆源30米时，桥墳处所受到的Mises应力最大 值、水击波压力最大值W及位移最大值分别为4.5234e^MPa、1.421MPa、3.99e^cm，当时间 达到t=1900化S后，桥墳处的Mises应力、水击波压力逐渐减小趋于稳定常数，表明此时刻 后桥墳不再受到爆破产生的水击波压力的影响；当减震孔的位置距离爆源40米时，桥墳处 所受到的Mises应力最大值、水击波压力最大值W及位移最大值分别为3.2743e-4M化、 1.963MPa、5.86e^cm;当减震孔的位置距离爆源49米时，桥墳处所受到的Mises应力最大值、 水击波压力最大值W及位移最大值分别为5.5423e^MPa、0.523MPa、6.65e^cm。
[0105] 通过所上述模拟结果的研究可得出结论:其一，无论减震孔的位置距离爆源的参 数如何布设，桥墳所受到的冲击波压力在其首次达到峰值后，将随着时间的增长而不断衰 减直至归零，而桥墳的位移将随之间的增长而加剧，不断累积直至达到稳定值;其二，当减 震孔距爆源40米时，桥墳处受到的冲击波压力与位移均大于减震孔距爆源30米时的竖直， 运表明水下爆破中减震孔的减震效果不是随着减震孔与爆源的距离增大而持续减弱，中间 可能存在着波动;其=，从整体数据看，随着减震孔与爆源距离的增大，桥墳处受到的Mises 应力与位移逐渐增大，运表明减震孔与爆源距离越大，减震效果越差;最后，综合上述的结 果，在现有建立的计算模型的基础上，可W得出将减震孔与爆源之间的距离设为15米，能够 达到较优的减震效果。
[0106] 其次，在第一实施例中所举出的【具体实施方式】的基础上，保持其余参数不变的情 况下，将减震孔的填充介质改为空气，建立相应的新的计算模型，同样采用ALE算法对新的 计算模型进行求解，得到新的模拟结果。
[0107] 具体地，当减震孔采用水介质填充时，Mises应力最大值、位移最大值W及水击波 压力最大值分别为2.151 X l(T3MPa、3.354cm、5.790 X ICT3MPa;当减震孔采用空气介质填充 时，Mises应力最大值、位移最大值W及水击波压力最大值分别为2.153 X l(T3MPa、5.107mm、 4.796X ICT3MPa,mm表示距离单位毫米。通过上述模拟结果可知，在其他因素条件相同时，减 震孔采用空气介质填充相比采用水介质填充能够更大程度的破坏岩层的整体性，从而在一 定程度上更容易的隔断地震波的传播，减震爆破引起的振动，因此该工程采用空气介质填 充减震孔。
[0108] 最后，同样采用单一控制变量法，对减震孔的几何尺寸进行模拟W得到优选值，依 次建立直径为110毫米、100毫米、90毫米W及80毫米的计算模型，利用ALE算法对新得到的 模型进行计算，得到相应的模拟结果，具体地：
[0109]
[0110] 表7
[0111] 通过上述模拟数值结果，并结合工程实际情况，可得出如下结论:其一，在各模型 的等效应力分析中，减震孔直径为80毫米时的Mises应力值是最大的，对桥墳的剪切破坏影 响最大，减震孔直径大于110毫米时，Mises应力值较小;其二，减震孔直径越大，其可压缩性 越大，减震孔产生位移也越大，由于减震孔中填充的介质具有缓冲作用，总体位移随着减震 孔尺寸增大稍有减小且效果不明显;其=，在爆破产生的冲击波对桥墳影响中，减震孔几何 尺寸过小对桥墳的减震效果达不到要求，减震孔尺寸达到一定数值后，再继续增大减震孔 尺寸，减震效果没有明显提升，且不符合经济要求;最后，综合各方面因素，减震孔的直径优 选值为110毫米。
[0112] 需要说明的是，本实施例给出的减震孔布设参数优选值是在建立的各个计算模型 基础上得到的，若采用不同的参数进行模拟，自然是得到不同的优选值，甚至是在计算模型 足够多时，能够获得最优值，本领域技术人员可根据具体情况进行相应的优选建模。
[0113] 请参阅图5,是本发明第四实施例提供的一种水下爆破减震孔布设模拟装置101。 下面将对图5所示的具体功能模块进行阐述。
[0114] 第一获得模块401，用于获得建立实体模型所需的初始布设参数和所述实体模型 包括的材料属性。具体地，该第一获得模块401可用于执行步骤SllO和步骤S210,具体的操 作方法可参照上述对步骤SllO和步骤S210的详细描述。
[0115] 第二获得模块402,用于根据所述材料属性获得对应的材料模型和状态方程。具体 地，该第二获得模块402可用于执行步骤Slll和步骤S211，具体的操作方法可参照上述对步 骤SI 11和步骤S211的详细描述。
[0116] 实体模型建立模块403,用于根据所述材料模型、所述状态方程中W及所述初始布 设参数建立实体模型。具体地，该实体模型建立模块403可用于执行步骤S112和步骤S212， 具体的操作方法可参照上述对步骤S112和步骤S212的详细描述。
[0117] 网格划分模块404,用于对所述实体模型进行网格划分，生成对应的有限元分析模 型;还用于再划分后得到的网格单元不符合预设计算精度要求时，对实体模型重新进行网 格划分。具体地，该网格划分模块404可用于执行步骤S113和步骤S213,具体的操作方法可 参照上述对步骤S113和步骤S213的详细描述。
[0118] 参数更改模块405,用于更改所述减震孔初始布设参数中的一个，得到新的有限元 分析模型。具体地，该参数更改模块405可用于执行步骤S310,具体的操作方法可参照上述 对步骤S310的详细描述。
[0119] 边界条件设定模块406,用于设定所述有限元分析模型的边界条件为无反射边界。 具体地，该边界条件设定模块406可用于执行步骤S215，具体的操作方法可参照上述对步骤 S215的详细描述。
[0120] 分析选项设置模块407,用于对计算所述有限元分析模型的分析选项进行设置，所 述分析选项包括求解时间和时间步。具体地，该分析选项设置模块407可用于执行步骤 S216,具体的操作方法可参照上述对步骤S216的详细描述。
[0121 ]计算模块408，用于采用连续介质力学运动描述方法对所述有限元分析模型进行 计算，得到水下爆破减震孔布设的模拟结果;还用于更改减震孔布设参数中的一个后，对新 得到的计算模型进行求解。具体地，该计算模块408可用于执行步骤S114、步骤S217W及步 骤S311，具体的操作方法可参照上述对步骤S114、步骤S217W及步骤S311的详细描述。
[0122] 检测模块409,用于检测生成的有限元分析模型中网格单元是否满足预设的计算 精度的要求。具体地，该检测模块409可用于执行步骤S214,具体的操作方法可参照上述对 步骤S214的详细描述。
[0123] 参数优选值确定模块410,用于根据新的模拟结果，确定所更改参数的优选值。具 体地，该参数优选值确定模块410可用于执行步骤S312,具体的操作方法可参照上述对步骤 S312的详细描述。
[0124] 综上所述，本发明提供的水下爆破减震孔布设模拟方法及装置，在通过炮孔及减 震孔的初始布设参数、材料模型及状态方程建立有限元分析模型并计算出相应的模拟结果 后，每一次更改减震孔初始布设参数中的一个，获得新的有限元分析模型，通过比较不同 计算模型所得出的模拟结果，得出所更改的减震孔布设参数的较优值，对提高水下爆破施 工效率、降低施工成本、确保紧邻建（构）筑物安全使用具有重要的促进意义，对进一步补 充、丰富水下爆破基本理论，具有积极的推进作用。
[0125] 在本申请所提供的几个实施例中，应该理解到，所掲露的装置和方法，也可W通过 其它的方式实现。W上所描述的装置实施例仅仅是示意性的，例如，附图中的流程图和框图 显示了根据本发明的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、 功能和操作。在运点上，流程图或框图中的每个方框可W代表一个模块、程序段或代码的一 部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执 行指令。也应当注意，在有些作为替换的实现方式中，方框中所标注的功能也可W W不同于 附图中所标注的顺序发生。例如，两个连续的方框实际上可W基本并行地执行，它们有时也 可W按相反的顺序执行，运依所设及的功能而定。也要注意的是，框图和/或流程图中的每 个方框、W及框图和/或流程图中的方框的组合，可W用执行规定的功能或动作的专用的基 于硬件的系统来实现，或者可W用专用硬件与计算机指令的组合来实现。
[0126] 另外，在本发明各个实施例中的各功能模块可W集成在一起形成一个独立的部 分，也可W是各个模块单独存在，也可W两个或两个W上模块集成形成一个独立的部分。
[0127] 所述功能如果W软件功能模块的形式实现并作为独立的产品销售或使用时，可W 存储在一个计算机可读取存储介质中。基于运样的理解，本发明的技术方案本质上或者说 对现有技术做出贡献的部分或者该技术方案的部分可WW软件产品的形式体现出来，该计 算机软件产品存储在一个存储介质中，包括若干指令用W使得一台计算机设备(可W是个 人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。 而前述的存储介质包括:U盘、移动硬盘、只读存储器、随机存取存储器、磁碟或者光盘等各 种可W存储程序代码的介质。
[0128] 需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实 体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示运些实体或操作之间存 在任何运种实际的关系或者顺序。而且，术语"包括"、"包含"或者其任何其他变体意在涵盖 非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要 素，而且还包括没有明确列出的其他要素，或者是还包括为运种过程、方法、物品或者设备 所固有的要素。在没有更多限制的情况下，由语句"包括一个……"限定的要素，并不排除在 包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
[0129] W上所述，仅为本发明的【具体实施方式】，但本发明的保护范围并不局限于此，任何 熟悉本技术领域的技术人员在本发明掲露的技术范围内，可轻易想到变化或替换，都应涵 盖在本发明的保护范围之内。因此，本发明的保护范围应所述W权利要求的保护范围为准。
【主权项】
1. 一种水下爆破减震孔布设模拟方法，其特征在于，所述方法包括： 获得建立实体模型所需的初始布设参数和所述实体模型包括的材料属性； 根据所述材料属性获得对应的材料模型和状态方程； 根据所述材料模型、所述状态方程以及所述初始布设参数建立实体模型； 对所述实体模型进行网格划分，生成对应的有限元分析模型； 采用连续介质力学运动描述方法对所述有限元分析模型进行计算，得到水下爆破减震 孔布设的模拟结果，所述模拟结果包括减震保护对象所受到的等效应力模拟数据、水击波 压力模拟数据以及位移距离模拟数据。2. 根据权利要求1所述的水下爆破减震孔布设模拟方法，其特征在于，所述初始布设参 数包括炮孔初始布设参数和减震孔初始布设参数，所述减震孔初始布设参数包括减震孔距 炮孔的距离、减震孔填充介质以及减震孔的直径，所述方法还包括： 更改所述减震孔初始布设参数中的一个，得到新的有限元分析模型； 对新得到的有限元分析模型采用同一连续介质力学运动描述方法进行计算，得到新的 模拟结果； 根据新的模拟结果，确定所更改参数的优选值。3. 根据权利要求1所述的水下爆破减震孔布设模拟方法，其特征在于，对所述有限元分 析模型进行计算之前还包括： 设定所述有限元分析模型的边界条件为无反射边界。4. 根据权利要求1所述的水下爆破减震孔布设模拟方法，其特征在于，对所述有限元分 析模型进行计算之前还包括： 对计算所述有限元分析模型的分析选项进行设置，所述分析选项包括求解时间和时间 步，所述求解时间为对所述有限元分析模型预设的计算时间，所述时间步为每一次对所述 实体模型进行初始化求解使其计算收敛。5. 根据权利要求1所述的水下爆破减震孔布设模拟方法，其特征在于，所述对所述实体 模型进行网格划分，生成对应的有限元分析模型的步骤之后还包括： 检测生成的有限元分析模型中网格单元是否满足预设的计算精度的要求； 若否，则调整网格大小，对所述实体模型重新进行网格划分。6. -种水下爆破减震孔布设模拟装置，其特征在于，所述装置包括： 第一获得模块，用于获得建立实体模型所需的初始布设参数和所述实体模型包括的材 料属性； 第二获得模块，用于根据所述材料属性获得对应的材料模型和状态方程； 实体模型建立模块，用于根据所述材料模型和所述状态方程以及所述初始布设参数建 立实体模型； 网格划分模块，用于对所述实体模型进行网格划分，生成对应的有限元分析模型； 计算模块，用于采用连续介质力学运动描述方法对所述有限元分析模型进行计算，得 到水下爆破减震孔布设的模拟结果，所述模拟结果包括减震保护对象所受到的等效应力模 拟数据、水击波压力模拟数据以及位移距离模拟数据。7. 根据权利要求6所述的水下爆破减震孔布设模拟装置，其特征在于，所述初始布设参 数包括炮孔初始布设参数和减震孔初始布设参数，所述减震孔初始布设参数包括减震孔距 炮孔的距离、减震孔填充介质以及减震孔的直径，所述装置还包括： 参数更改模块，用于更改所述减震孔初始布设参数中的一个，得到新的有限元分析模 型； 所述计算模块，还用于对新得到的有限元分析模型采用同一连续介质力学运动描述方 法进行计算，得到新的模拟结果； 参数优选值确定模块，用于根据新的模拟结果，确定所更改参数的优选值。8. 根据权利要求6所述的水下爆破减震孔布设模拟装置，其特征在于，所述装置还包 括： 边界条件设定模块，用于设定所述有限元分析模型的边界条件为无反射边界。9. 根据权利要求6所述的水下爆破减震孔布设模拟装置，其特征在于，所述装置还包 括： 分析选项设置模块，用于对计算所述有限元分析模型的分析选项进行设置，所述分析 选项包括求解时间和时间步，所述求解时间为对所述有限元分析模型预设的计算时间，所 述时间步为每一次对所述实体模型进行初始化求解使其计算收敛。10. 根据权利要求6所述的水下爆破减震孔布设模拟装置，其特征在于，所述装置还包 括： 检测模块，用于检测生成的有限元分析模型中网格单元是否满足预设的计算精度的要 求； 所述网格划分模块，还用于在划分的网格单元不满则预设的计算精度时，对所述实体 模型重新进行网格划分。
【文档编号】G06F17/50GK105956353SQ201610564247
【公开日】2016年9月21日
【申请日】2016年7月18日
【发明人】李春军, 舒利, 吴立, 李波, 李红勇, 彭亚雄, 陈春晖, 胡伟才, 叶小敏, 冉华桥, 袁青, 苏莹, 周瑞锋, 周玉纯, 汪煜烽, 马晨阳, 李平, 杜健如
【申请人】长江重庆航道工程局, 中国地质大学（武汉）