地应力影响下岩石爆破损伤模型构建与数值试验研究

地应力影响下岩石爆破损伤模型构建与数值试验研究一、引言1.1研究背景与意义在现代工程建设中，岩石爆破是一项极为关键且应用广泛的技术手段，对推动工程进展起着不可或缺的作用。在矿山开采领域，岩石爆破是实现矿石高效开采的核心环节，直接关系到矿产资源的获取效率与成本。例如，在大型露天矿山中，通过精确的爆破作业，能够将大量的矿石从岩体中分离出来，为后续的选矿和加工提供原料，满足工业生产对矿产资源的需求。在交通基础设施建设方面，无论是公路、铁路的隧道挖掘，还是桥梁基础的建设，岩石爆破技术都发挥着重要作用。比如，在山区修建公路时，常常需要通过爆破来开辟道路，打通山体，为公路的顺利铺设创造条件。在水利水电工程中，岩石爆破对于大坝基础的开挖、地下厂房的建设等都具有重要意义。以三峡大坝的建设为例，在基础开挖过程中，运用了大规模的岩石爆破技术，确保了基础的稳定性和工程的顺利进行。地应力作为岩石所处地质环境中的固有应力，对岩石爆破损伤有着极为重要的影响，其作用不可忽视。地应力的存在使得岩石内部处于一种复杂的应力状态，这种应力状态会改变岩石的力学性质，进而对爆破过程中岩石的损伤和破坏机制产生显著影响。在高地应力区域进行爆破作业时，地应力可能会导致岩石在爆破前就已经存在一定程度的潜在损伤，使得岩石的强度和稳定性降低。在爆破过程中，地应力与爆炸应力波相互作用，会改变裂纹的扩展方向和速度，影响岩石的破碎效果。若地应力方向与爆破方向不一致，可能会使裂纹偏向地应力方向扩展，导致爆破效果不理想，甚至可能引发安全事故。研究地应力影响下的岩石爆破损伤模型，对于工程实践具有多方面的重要指导意义。准确掌握地应力对岩石爆破损伤的影响规律，能够为爆破参数的优化设计提供科学依据。通过合理调整爆破参数，如炸药的用量、炮孔的布置、起爆顺序等，可以使爆破能量更有效地作用于岩石，提高爆破效率，减少炸药的浪费和对周围环境的影响。在实际工程中，根据不同的地应力条件和岩石特性，精确计算和调整爆破参数，能够使岩石破碎更加均匀，降低大块率，提高工程质量和进度。对岩石爆破损伤模型的深入研究有助于提高工程的安全性。了解地应力作用下岩石的损伤演化过程，可以提前预测爆破可能引发的安全隐患，如岩石的坍塌、飞石的飞溅等，从而采取相应的防护措施，保障施工人员和周边环境的安全。在隧道爆破施工中，通过对岩石爆破损伤模型的分析，合理设计支护结构和施工工艺，能够有效预防隧道坍塌等事故的发生。1.2国内外研究现状在理论研究方面，众多学者围绕地应力对岩石爆破损伤的影响机制展开了深入探索。早期，一些学者基于经典的弹性力学和断裂力学理论，尝试分析地应力与爆破应力波的相互作用原理。他们认为，地应力的存在会改变岩石内部的初始应力状态，使得爆破应力波在传播过程中受到干扰，从而影响岩石的裂纹起裂和扩展。随着研究的不断深入，损伤力学理论逐渐被引入到该领域，为建立岩石爆破损伤模型提供了新的思路。学者们开始从微观层面考虑岩石内部的损伤演化过程，通过定义损伤变量来描述岩石在爆破和地应力共同作用下的力学性能劣化。有学者提出了基于应变等效原理的损伤模型，将地应力作为影响损伤演化的重要因素，分析了不同地应力条件下岩石的损伤阈值和损伤扩展规律。然而，现有的理论研究虽然取得了一定成果，但在考虑地应力的复杂性和岩石材料的多样性方面仍存在不足，部分理论模型的假设条件与实际工程情况存在一定偏差，导致其在实际应用中的准确性和可靠性受到限制。在试验研究领域，国内外学者开展了大量的物理试验，旨在获取地应力影响下岩石爆破损伤的直接数据和现象。早期的试验主要集中在研究不同地应力水平下岩石的基本力学性质变化，如抗压强度、抗拉强度等。通过对岩石试件施加不同方向和大小的地应力，然后进行爆破试验，观察岩石的破坏形态和损伤特征。随着试验技术的不断进步，一些先进的测试手段，如声发射技术、数字图像相关技术等被广泛应用于岩石爆破损伤试验中。声发射技术可以实时监测岩石在爆破过程中的裂纹产生和扩展，获取岩石内部损伤的动态信息；数字图像相关技术则能够精确测量岩石表面的变形和位移，为分析岩石的损伤演化提供直观的数据支持。有研究利用声发射技术监测了高地应力条件下岩石爆破过程中的能量释放规律，发现地应力会显著影响岩石的声发射特征，从而反映出岩石内部损伤的差异。然而，物理试验往往受到试验条件和设备的限制，难以完全模拟实际工程中的复杂地质条件和爆破工况，且试验成本较高、周期较长，这在一定程度上限制了试验研究的全面性和深入性。数值模拟作为一种重要的研究手段，在分析地应力对岩石爆破损伤的影响方面发挥了重要作用。随着计算机技术的飞速发展，各种数值模拟软件不断涌现，如ANSYS、LS-DYNA、FLAC3D等，为岩石爆破损伤的数值模拟提供了强大的工具。学者们通过建立岩石的数值模型，考虑地应力的作用，模拟爆破过程中岩石的应力、应变分布以及裂纹的扩展和贯通，从而预测岩石的爆破损伤范围和程度。在一些研究中，利用有限元软件建立了三维岩石模型，考虑了地应力的各向异性和非线性特性，模拟了不同地应力条件下的爆破过程，分析了地应力对爆破效果的影响规律。尽管数值模拟在研究地应力影响下岩石爆破损伤方面具有诸多优势，但数值模型的建立往往依赖于对岩石材料参数的准确获取，而实际岩石材料的参数具有较大的离散性和不确定性，这会导致数值模拟结果与实际情况存在一定的误差。此外，现有的数值模拟方法在处理复杂的岩石结构和多物理场耦合问题时还存在一定的困难，需要进一步改进和完善。综合来看，目前对于地应力影响下岩石爆破损伤的研究虽然已经取得了一定的进展，但仍存在一些不足之处。现有研究在理论模型的普适性、试验研究的全面性以及数值模拟的准确性等方面还需要进一步提高。在理论模型方面，需要建立更加完善的考虑地应力复杂作用机制和岩石材料非线性特性的损伤模型；在试验研究方面，需要开发更加先进的试验技术和设备，以更好地模拟实际工程条件，获取更准确的试验数据；在数值模拟方面，需要进一步优化数值算法和模型参数，提高模拟结果的可靠性和精度。基于此，本文将针对现有研究的不足，深入开展地应力影响下岩石爆破损伤模型及其数值试验的研究，旨在为岩石爆破工程提供更加科学、准确的理论支持和技术指导。1.3研究内容与方法本文主要研究地应力影响下的岩石爆破损伤模型及其数值试验，具体研究内容如下：岩石爆破损伤理论分析：深入剖析地应力对岩石力学性质的影响机制，包括地应力如何改变岩石的强度、弹性模量、泊松比等参数。探究地应力与爆破应力波的相互作用原理，分析这种相互作用对岩石裂纹起裂、扩展和贯通的影响，为后续建立岩石爆破损伤模型奠定坚实的理论基础。岩石爆破损伤模型构建：基于损伤力学理论，充分考虑地应力的复杂作用以及岩石材料的非线性特性，引入合适的损伤变量来描述岩石在爆破过程中的损伤演化。通过理论推导和数学建模，建立能够准确反映地应力影响下岩石爆破损伤过程的模型，明确模型中各参数的物理意义和取值方法。数值试验方案设计：运用数值模拟软件，构建合理的岩石爆破数值模型，设定不同的地应力条件、爆破参数和岩石材料参数，如地应力的大小、方向、侧压系数，炸药的类型、装药量、起爆方式，岩石的抗压强度、抗拉强度、密度等。设计全面且有针对性的数值试验方案，以系统研究这些因素对岩石爆破损伤的影响规律。数值试验结果分析：对数值试验结果进行深入分析，获取岩石在爆破过程中的应力、应变分布云图，直观展示地应力和爆破作用下岩石内部的力学响应。分析裂纹的扩展路径和损伤区域的分布特征，研究地应力如何影响裂纹的起裂位置、扩展方向和速度，以及损伤区域的范围和形状。通过数据统计和对比分析，总结地应力与爆破参数、岩石材料参数之间的定量关系，为工程实践提供数据支持。模型验证与工程应用：收集实际工程中的岩石爆破案例数据，将建立的岩石爆破损伤模型应用于实际工程案例的模拟分析，对比模拟结果与实际观测数据，验证模型的准确性和可靠性。根据模型验证结果，对模型进行优化和改进，使其能够更好地适应不同的工程地质条件和爆破工况。基于验证后的模型，为实际工程中的爆破参数优化设计提供科学依据，提出合理的爆破方案建议，提高工程爆破的效率和安全性。本文采用理论分析、数值模拟和案例验证相结合的研究方法：理论分析：通过查阅大量国内外相关文献资料，深入研究岩石力学、爆炸力学、损伤力学等领域的基础理论，梳理地应力影响下岩石爆破损伤的研究现状和发展趋势。运用弹性力学、断裂力学等理论知识，分析地应力与爆破应力波的相互作用机制，推导岩石爆破损伤模型的基本方程，为研究提供坚实的理论支撑。数值模拟：选用专业的数值模拟软件，如ANSYS、LS-DYNA、FLAC3D等，根据岩石的物理力学性质和爆破工程的实际情况，建立三维岩石爆破数值模型。在模型中准确施加地应力和爆破荷载，设置合理的材料参数和边界条件，模拟不同工况下的岩石爆破过程。利用数值模拟软件的后处理功能，对模拟结果进行可视化处理和数据分析，获取岩石爆破损伤的相关信息。案例验证：收集矿山开采、隧道挖掘、水利水电工程等实际工程中的岩石爆破案例，详细记录工程的地质条件、地应力测量数据、爆破参数以及爆破后的岩石损伤情况等信息。将建立的岩石爆破损伤模型应用于这些实际案例的模拟分析，将模拟结果与实际工程数据进行对比验证，评估模型的准确性和适用性。根据案例验证结果，对模型进行修正和完善，使其能够更好地应用于实际工程。二、地应力与岩石爆破损伤相关理论基础2.1地应力概述2.1.1地应力的概念与组成地应力是存在于地壳中的未受工程扰动的天然应力，也被称作岩体初始应力、绝对应力或原岩应力，从广义角度来讲，其还涵盖地球体内的应力。它是岩石在漫长地质历史时期中，历经多种复杂地质作用而形成的应力状态，是多种起源应力的综合体现。地应力主要由自重应力和构造应力组成，同时还包括孔隙压力、热应力和残余应力等，其中，重力应力（即自重应力）和构造应力是地应力的主要来源。自重应力是由岩石自身重量产生的应力，其形成机制相对直观。在地球引力的作用下，上覆岩石的重量对下部岩石产生压力，从而形成自重应力。以深度为z处的某一点为例，其垂直方向的自重应力\sigma_{z}可通过公式\sigma_{z}=\gammaz计算得出，其中\gamma为岩石的重度。在水平方向上，由于岩石的泊松效应，会产生一定的水平自重应力，其大小通常与垂直自重应力存在一定的比例关系，可表示为\sigma_{x}=\sigma_{y}=\frac{\mu}{1-\mu}\sigma_{z}，其中\mu为岩石的泊松比。构造应力则是由于地壳运动、板块碰撞、褶皱、断裂等地质构造运动所产生的应力。这种应力的形成机制较为复杂，涉及到地球内部的动力学过程和岩石圈的变形。板块碰撞时，巨大的挤压力会使岩石发生变形和破裂，从而在岩石内部产生构造应力。构造应力在空间分布上具有明显的不均匀性，其大小和方向会受到地质构造的控制。在褶皱构造区域，岩石会受到弯曲和挤压，导致构造应力集中在褶皱的轴部和翼部；在断裂构造附近，由于岩石的错动和变形，构造应力也会呈现出复杂的分布状态。构造应力的方向往往与地质构造的走向相关，在板块碰撞带，构造应力的主方向通常与板块运动方向一致。孔隙压力是指岩石孔隙中流体（如水、石油、天然气等）所产生的压力。当岩石孔隙中充满流体时，流体压力会对岩石骨架产生作用，从而影响地应力的分布。在石油开采过程中，随着油井的开采，孔隙中的油被抽出，孔隙压力降低，会导致岩石骨架所承受的有效应力增加，进而可能引发岩石的变形和破坏。热应力是由于岩石温度变化而产生的应力。当地下岩石的温度发生变化时，岩石会发生热胀冷缩，若这种变形受到约束，就会在岩石内部产生热应力。在地下热水开采区域，由于热水的抽出导致岩石温度降低，可能会引发热应力，对周围岩石的稳定性产生影响。残余应力是指岩石在形成过程中或经历地质构造运动后，由于岩石内部的不均匀变形等原因而残留下来的应力。这些残余应力在一定程度上也会影响岩石的力学性质和地应力场的分布。2.1.2地应力的测量方法地应力的准确测量对于深入了解地质构造、保障工程安全等方面具有重要意义。目前，常用的地应力测量方法包括水压致裂法、应力解除法等，每种方法都有其独特的原理、操作流程和优缺点。水压致裂法是一种应用较为广泛的地应力测量方法，其原理基于岩石的水力破裂现象。在钻孔中，利用封隔器将一段钻孔封闭，然后向封隔段内注入高压水，随着水压的逐渐升高，当水压达到一定程度时，岩石会在钻孔壁上产生破裂。通过测量水压致裂过程中的各个参数，如破裂压力、关闭压力等，再依据相应的理论公式，即可计算出地应力的大小和方向。具体操作流程如下：首先，在预定测量地点进行钻孔，钻孔深度需根据实际测量需求确定；接着，将封隔器下入钻孔中，使其位于待测位置，并将封隔段密封；随后，通过高压泵向封隔段内注入水，同时监测水压的变化以及钻孔壁的变形情况；当观察到岩石发生破裂时，记录此时的破裂压力P_f；之后，停止注水，压力会逐渐下降，当压力稳定时，记录关闭压力P_s。根据水压致裂理论，水平最小主应力\sigma_{h}可通过公式\sigma_{h}=P_s计算得出，而水平最大主应力\sigma_{H}则可通过公式\sigma_{H}=3\sigma_{h}-P_f+T计算，其中T为岩石的抗拉强度。水压致裂法的优点在于其测量深度较大，可达到数千米，适用于深部地应力测量；且操作相对简便，不需要复杂的设备和技术。然而，该方法也存在一定的局限性，它假设岩石为各向同性的弹性体，这在实际地质条件下往往难以完全满足，可能导致测量结果存在一定误差；同时，该方法无法直接测量垂直应力，需要结合其他方法进行综合确定。应力解除法是目前国内外应用最为普遍的一种地应力测量方法，其原理基于弹性力学的应力-应变关系。该方法通过在岩石中钻出一个小孔，然后将小孔周围的岩石进行套芯解除，使小孔周围的岩石应力得到释放，产生弹性恢复变形。利用安装在小孔内的应变计测量这种变形，再根据岩石的弹性参数，通过计算即可得到地应力的大小和方向。以空心包体应变计应力解除法为例，其操作流程如下：首先，在岩石表面钻一个大孔，然后在大孔底部钻出一个同心的小孔；将空心包体应变计安装在小孔中，应变计由多个应变片组成，可测量不同方向的应变；接着，采用套芯技术，将包含应变计的岩石芯体从周围岩体中取出，使岩石芯体的应力得到解除；通过测量应变计的应变变化，利用弹性力学公式进行计算，可得到地应力的各个分量。应力解除法的优点是测量精度较高，能够较为准确地测量地应力的大小和方向；并且可以测量三维地应力状态。但该方法的缺点是测量过程较为复杂，对设备和技术要求较高；测量深度相对有限，一般适用于浅部地应力测量；而且在测量过程中，岩石的扰动可能会对测量结果产生一定影响。2.2岩石爆破损伤理论2.2.1岩石爆破损伤的机理岩石爆破损伤是一个极为复杂的物理力学过程，涉及爆炸应力波、爆生气体等多种因素的综合作用，这些因素相互交织，共同导致了岩石内部结构的破坏和力学性能的劣化。爆炸应力波在岩石爆破损伤过程中扮演着重要角色。当炸药在岩石内部爆炸时，瞬间释放出巨大的能量，形成高温高压的爆轰产物，这些爆轰产物迅速膨胀，产生强烈的冲击载荷，进而在岩石中激发出爆炸应力波。爆炸应力波以极高的速度在岩石中传播，其传播过程伴随着能量的传递和衰减。在传播初期，应力波的峰值应力极高，会使岩石内部产生强烈的压缩和拉伸变形。当岩石所承受的应力超过其动态抗压强度时，在炮孔周围会形成压碎区，岩石结构被严重破坏，呈现出粉碎状。随着应力波的进一步传播，其能量逐渐衰减，当应力波传播到自由面时，会发生反射，反射波与入射波相互叠加，在岩石内部形成复杂的应力状态。由于岩石的抗拉强度远低于抗压强度，在反射拉伸波的作用下，岩石内部会产生大量的拉伸裂纹，这些裂纹不断扩展、连接，逐渐形成裂隙区，导致岩石的完整性受到严重破坏。爆生气体同样对岩石爆破损伤有着重要影响。炸药爆炸后，除了产生爆炸应力波外，还会生成大量的高温高压爆生气体。这些爆生气体在炮孔内迅速膨胀，对炮孔壁产生持续的高压作用。在爆生气体的高压作用下，炮孔周围的岩石进一步被压实和破碎，压碎区的范围可能会进一步扩大。爆生气体还会沿着岩石内部已有的裂纹和孔隙渗透，在裂纹尖端产生楔入作用，使裂纹进一步扩展。爆生气体的膨胀作用还会导致岩石内部的应力重新分布，促进裂纹的相互连接和贯通，形成更大范围的损伤区域。在一些情况下，爆生气体的膨胀作用甚至可以使岩石产生抛掷和破碎，形成爆破漏斗。爆炸应力波和爆生气体在岩石爆破损伤过程中并非孤立作用，而是相互协同、相互影响。爆炸应力波首先在岩石中传播，使岩石产生初始的损伤和裂纹，为爆生气体的侵入和作用提供了通道。爆生气体随后沿着这些裂纹和孔隙渗透，进一步加剧岩石的损伤和破坏，其膨胀作用还会改变岩石内部的应力状态，影响后续爆炸应力波的传播和作用效果。这种相互作用的过程使得岩石爆破损伤呈现出复杂的演化特征，不同阶段的主导因素也有所不同。在爆破初期，爆炸应力波的作用较为突出，主要负责岩石的初始破碎和裂纹的萌生；随着时间的推移，爆生气体的作用逐渐增强，对岩石的后续损伤和破碎起到关键作用。2.2.2岩石爆破损伤的影响因素岩石爆破损伤受到多种因素的综合影响，这些因素涵盖了岩石自身性质、炸药特性、爆破参数以及地应力等多个方面，它们相互作用，共同决定了岩石爆破损伤的程度和特征。岩石性质是影响爆破损伤的内在因素，不同岩石因其矿物成分、结构构造和力学性能的差异，在爆破过程中的损伤表现各不相同。从矿物成分来看，含有较多脆性矿物的岩石，如石英含量高的砂岩，在爆破时更容易产生裂纹和破碎，因为脆性矿物在受到应力作用时难以发生塑性变形，容易达到其强度极限而破裂。而含有较多韧性矿物的岩石，如黏土矿物含量较高的页岩，相对较难破碎，爆破损伤程度相对较小，韧性矿物能够吸收一部分能量，延缓裂纹的扩展。岩石的结构构造也对爆破损伤有重要影响，具有层理、节理等结构面的岩石，爆破时应力波会在这些结构面处发生反射、折射和绕射，导致应力集中，使得岩石更容易沿着结构面发生破裂和损伤。节理间距较小、节理组数较多的岩石，在爆破后往往会形成更多的碎块，损伤范围也更大。岩石的力学性能，如抗压强度、抗拉强度、弹性模量等，直接决定了岩石抵抗爆破应力的能力。抗压强度和抗拉强度较低的岩石，在相同爆破条件下更容易受到损伤，而弹性模量较大的岩石，能够更快地传播应力波，但也可能导致应力集中更加明显，从而增加岩石的损伤程度。炸药特性对岩石爆破损伤起着关键作用，不同类型的炸药，其爆炸性能和能量释放特性存在差异，进而影响爆破效果和岩石损伤情况。炸药的爆速是一个重要参数，爆速较高的炸药能够产生更强烈的爆炸应力波，对岩石的冲击作用更大，更容易使岩石破碎。高爆速炸药在爆炸瞬间释放出的能量更集中，能够在短时间内使岩石达到更高的应力状态，促使岩石产生更多的裂纹和破碎。炸药的爆轰压力也直接影响岩石的爆破损伤，爆轰压力越大，对炮孔壁的冲击作用越强，岩石在高压作用下更容易被压碎和破坏。不同的炸药配方和装药结构也会影响炸药的能量利用率和爆炸作用方式。采用合理的装药结构，如不耦合装药，可以降低爆炸应力波对炮孔壁的直接冲击，减少岩石的过度破碎，同时使爆炸能量更有效地向远处传播，扩大岩石的损伤范围。爆破参数是人为可控的因素，对岩石爆破损伤有着直接且显著的影响。炮孔间距的大小决定了相邻炮孔爆炸应力波的相互作用程度。炮孔间距过小，会导致爆炸应力波相互叠加过度，使岩石过度破碎，增加炸药的浪费和工程成本；炮孔间距过大，则可能使岩石破碎不均匀，出现大块岩石，影响爆破效果。合理的炮孔间距应根据岩石性质、炸药特性等因素进行优化确定。装药结构也对爆破损伤有着重要影响，除了前面提到的不耦合装药结构外，连续装药和间隔装药等不同的装药方式也会影响爆炸能量的分布和作用效果。连续装药能够使爆炸能量在炮孔内连续释放，对炮孔周围岩石的作用较为集中；间隔装药则可以使爆炸能量在不同位置分段释放，有利于岩石的均匀破碎，减少大块率。起爆顺序的设计对于控制岩石的破裂过程和损伤范围至关重要。采用合理的起爆顺序，如微差起爆，可以使相邻炮孔的爆炸应力波在时间和空间上相互叠加，形成有利的应力场，促进岩石的破碎和裂纹的贯通，提高爆破效率。地应力作为岩石所处地质环境中的固有应力，对岩石爆破损伤有着不可忽视的影响。地应力的存在使得岩石内部处于一种初始应力状态，这种应力状态会改变岩石在爆破过程中的力学响应。在高地应力区域，地应力会增加岩石的抗压强度，使得岩石在爆破时更难破碎。地应力还会影响爆炸应力波的传播和裂纹的扩展方向。当爆炸应力波传播方向与地应力方向一致时，地应力会对爆炸应力波起到增强作用，使岩石更容易沿着该方向产生裂纹和破碎；当爆炸应力波传播方向与地应力方向垂直时，地应力会对裂纹的扩展产生阻碍作用，改变裂纹的扩展路径，可能导致爆破效果不理想。地应力的大小和方向还会影响岩石在爆破后的稳定性，高地应力条件下爆破后的岩石更容易发生垮塌等安全事故。2.3数值试验方法基础2.3.1有限元法原理与应用有限元法作为一种强大的数值计算方法，在工程领域中应用广泛，其基本原理基于变分原理和离散化思想。在岩石爆破数值模拟中，有限元法能够有效地分析岩石在爆破过程中的力学响应，为研究岩石爆破损伤提供了重要的手段。有限元法的基本原理是将连续的求解区域离散为有限个相互连接的单元，这些单元通过节点相互连接。在每个单元内，选择合适的插值函数来近似表示单元内的未知场变量，如位移、应力等。通过将实际问题的控制方程（如平衡方程、几何方程和物理方程）在每个单元上进行离散化处理，利用虚功原理或变分原理，建立起单元的刚度方程。对于岩石爆破问题，其控制方程主要包括动力学方程、连续方程和状态方程。动力学方程描述了岩石在爆破应力作用下的运动状态，如\rho\frac{\partial^2u_i}{\partialt^2}=\frac{\partial\sigma_{ij}}{\partialx_j}+f_i，其中\rho为岩石密度，u_i为位移分量，\sigma_{ij}为应力分量，f_i为体积力分量；连续方程反映了岩石的质量守恒，\frac{\partial\rho}{\partialt}+\frac{\partial(\rhov_i)}{\partialx_i}=0，v_i为速度分量；状态方程则描述了岩石的应力-应变关系，如\sigma_{ij}=C_{ijkl}\varepsilon_{kl}，C_{ijkl}为弹性常数，\varepsilon_{kl}为应变分量。通过对这些方程进行离散化，得到单元的刚度方程K^e\delta^e=F^e，K^e为单元刚度矩阵，\delta^e为单元节点位移向量，F^e为单元节点力向量。将所有单元的刚度方程进行组装，形成整个结构的平衡方程组K\delta=F，K为总体刚度矩阵，\delta为总体节点位移向量，F为总体节点力向量。通过求解这个方程组，就可以得到整个结构在给定荷载作用下的节点位移，进而计算出应力、应变等物理量。在岩石爆破数值模拟中，有限元法具有诸多优势。它能够精确地模拟复杂的岩石几何形状和边界条件，对于具有不规则形状的岩石结构体或存在多个炮孔的爆破模型，有限元法可以通过灵活的网格划分技术，准确地描述其几何特征，从而更真实地反映爆破过程中的力学行为。有限元法还可以方便地考虑多种物理场的耦合作用，如在岩石爆破中，爆炸应力波的传播会引起岩石的温度变化，同时岩石的损伤和破坏也会影响其热传导性能，有限元法能够将力学场和温度场等进行耦合分析，全面地研究爆破过程中的多物理现象。通过有限元模拟，可以直观地得到岩石在爆破过程中的应力、应变分布云图，清晰地展示出爆破应力波的传播路径和岩石内部的应力集中区域，以及岩石的变形和破坏过程，为分析爆破效果和优化爆破参数提供了直观的数据支持。然而，有限元法在岩石爆破数值模拟中也存在一些局限性。在处理岩石的大变形和破坏问题时，由于有限元法基于小变形假设，当岩石发生较大变形或破裂时，其计算精度会受到一定影响。岩石在爆破过程中会产生大量的裂纹和破碎，传统的有限元方法难以准确地模拟裂纹的扩展和贯通，需要结合其他方法，如扩展有限元法等进行改进。有限元模拟的计算精度依赖于网格的划分质量，在复杂的岩石结构中，为了保证计算精度，往往需要划分大量的细小网格，这会导致计算量急剧增加，计算效率降低，对计算机的硬件性能提出了较高的要求。2.3.2离散元法原理与应用离散元法是一种适用于模拟离散介质力学行为的数值方法，在分析岩石颗粒间相互作用和裂纹扩展方面具有独特的优势，为研究地应力影响下的岩石爆破损伤提供了新的视角。离散元法的基本原理是将所研究的对象视为由一系列离散的颗粒组成，这些颗粒之间通过接触力相互作用。每个颗粒都具有独立的质量、速度、位移等物理量，通过牛顿第二定律来描述颗粒的运动状态。在离散元模型中，颗粒之间的接触力模型是关键，常用的接触力模型有线性弹簧-阻尼模型、赫兹接触模型等。以线性弹簧-阻尼模型为例，当两个颗粒发生接触时，接触力由弹簧力和阻尼力组成，弹簧力用于模拟颗粒之间的弹性相互作用，F_s=k\Delta\delta，k为弹簧刚度，\Delta\delta为颗粒间的相对位移；阻尼力用于模拟颗粒之间的能量耗散，F_d=c\Deltav，c为阻尼系数，\Deltav为颗粒间的相对速度。通过不断更新颗粒的位置、速度和接触力，来模拟整个离散系统的动态演化过程。在模拟岩石颗粒间相互作用方面，离散元法能够直观地展现颗粒之间的接触状态和力的传递过程。在岩石爆破过程中，爆炸应力波会使岩石颗粒产生运动和相互碰撞，离散元法可以精确地模拟这些颗粒的运动轨迹和碰撞力，分析颗粒间的能量传递和耗散机制。通过离散元模拟，可以观察到岩石颗粒在爆破作用下的排列变化和孔隙结构的演化，从而深入了解岩石的力学性质变化。在模拟裂纹扩展方面，离散元法具有独特的优势。当岩石内部的应力超过颗粒间的接触强度时，颗粒之间的连接会被破坏，从而形成裂纹。离散元法能够自然地模拟裂纹的萌生和扩展过程，通过追踪颗粒间连接的断裂情况，可以清晰地描绘出裂纹的扩展路径和形态。在高地应力条件下，离散元法可以考虑地应力对颗粒间接触力的影响，分析地应力如何改变裂纹的扩展方向和速度，以及对岩石整体破碎模式的影响。离散元法在岩石爆破研究中得到了广泛的应用。在研究岩石的爆破破碎机理方面，离散元模拟可以帮助研究人员深入了解爆炸应力波与岩石颗粒相互作用的微观过程，揭示岩石破碎的内在机制。通过对比不同地应力条件下的离散元模拟结果，可以分析地应力对岩石爆破破碎效果的影响规律，为优化爆破参数提供理论依据。在工程应用中，离散元法可以用于预测爆破后的岩石块度分布，根据模拟结果调整爆破参数，以达到预期的岩石破碎效果，提高工程效率和经济效益。离散元法也存在一些不足之处，由于需要对大量的颗粒进行计算，其计算量较大，计算时间较长，在处理大规模岩石爆破问题时，计算效率有待提高。离散元模型中颗粒的物理参数和接触参数的确定较为困难，这些参数的准确性对模拟结果的可靠性有较大影响。三、地应力影响下岩石爆破损伤模型构建3.1模型假设与基本思路3.1.1模型假设条件为了构建地应力影响下的岩石爆破损伤模型，首先需要明确一系列假设条件，这些假设条件是简化复杂实际问题、建立理论模型的基础，同时也界定了模型的适用范围。假设岩石为连续介质，忽略岩石内部微观结构的离散性和非均匀性。在实际岩石中，存在着矿物颗粒、孔隙、微裂纹等微观结构，这些微观结构的存在会导致岩石力学性质的局部变化和非均匀性。但在本模型中，为了便于分析和计算，假设岩石在宏观上是连续的，即认为岩石内部的物理量（如应力、应变等）在空间上是连续分布的。这一假设在一定程度上简化了问题的复杂性，使得基于连续介质力学的理论和方法能够应用于岩石爆破损伤的研究中。当岩石内部的微观结构尺寸远小于所研究问题的特征尺寸时，这一假设是合理的，例如在研究大规模的岩石爆破工程时，岩石内部微观结构的局部影响可以被平均化，从而可以将岩石视为连续介质进行分析。假设岩石为各向同性材料，即岩石在各个方向上的力学性质相同。实际上，许多岩石由于其内部结构的定向性，如层理、节理等，呈现出各向异性的力学性质。然而，在一些情况下，当岩石的各向异性特征不显著，或者为了初步研究地应力对岩石爆破损伤的一般性影响时，假设岩石为各向同性可以简化模型的建立和分析过程。对于一些较为均匀的块状岩石，在不考虑其内部细微结构影响的前提下，将其视为各向同性材料能够满足一定的工程精度要求。但需要注意的是，当岩石的各向异性对爆破损伤有显著影响时，这一假设可能会导致模型结果与实际情况存在偏差，此时需要进一步考虑岩石的各向异性特性，采用更为复杂的各向异性模型进行研究。假设爆破过程为瞬间完成，不考虑爆破过程中的时间效应。在实际爆破中，炸药爆炸、应力波传播以及岩石损伤破坏是一个动态的时间过程，涉及到复杂的动力学响应。但为了简化模型，本假设将爆破过程视为在极短时间内瞬间完成，忽略了爆破过程中的时间因素对岩石损伤的影响。这一假设适用于一些对爆破过程时间细节要求不高，主要关注爆破最终损伤结果的研究。在分析爆破后岩石的整体破碎形态和损伤范围时，忽略时间效应可以在一定程度上简化计算，提高分析效率。然而，对于一些需要精确研究爆破动态过程的问题，如应力波的传播规律、裂纹的动态扩展过程等，这一假设则不再适用，需要采用考虑时间效应的动力学模型进行深入研究。本模型主要适用于脆性岩石，对于塑性变形较大的岩石，模型的适用性需要进一步验证和修正。脆性岩石在爆破过程中主要表现为裂纹的产生和扩展，最终导致岩石的破碎，其力学行为相对较为简单，基于上述假设建立的模型能够较好地描述其爆破损伤过程。而塑性变形较大的岩石，在受到爆破荷载作用时，除了裂纹扩展外，还会发生显著的塑性流动和变形，其力学行为更为复杂，需要考虑更多的因素，如岩石的塑性本构关系、加工硬化等，本模型在处理这类岩石时可能存在局限性，需要进行针对性的改进和完善。3.1.2模型构建的基本思路构建地应力影响下的岩石爆破损伤模型的总体思路是综合考虑地应力、岩石性质和爆破荷载等多方面因素，从损伤力学的角度出发，建立能够准确描述岩石在爆破过程中损伤演化的数学模型。地应力是岩石所处地质环境中的固有应力，其大小和方向对岩石的力学性质和爆破损伤过程有着重要影响。在模型构建过程中，需要充分考虑地应力的作用。通过地应力测量数据，获取岩石所处位置的初始地应力状态，包括三个主应力的大小和方向。将地应力作为岩石的初始应力条件施加到模型中，分析地应力如何改变岩石的力学性能，如强度、弹性模量等。在高地应力条件下，岩石的抗压强度会增加，这是因为地应力的存在使得岩石内部的颗粒之间的接触更加紧密，抵抗外力的能力增强；同时，地应力还会影响岩石的泊松比，改变岩石在受力时的横向变形特性。地应力对爆破应力波的传播和裂纹扩展也有着重要影响，它会改变应力波的传播速度和方向，使得裂纹更容易沿着地应力的方向扩展。岩石性质是影响爆破损伤的重要内在因素。不同类型的岩石，其矿物成分、结构构造和力学性能存在差异，这些差异会导致岩石在爆破过程中的损伤表现各不相同。在模型中，需要准确描述岩石的性质。通过岩石力学试验，获取岩石的基本力学参数，如抗压强度、抗拉强度、弹性模量、泊松比等。这些参数是建立岩石本构关系的基础，本构关系描述了岩石在受力过程中应力与应变之间的关系，它反映了岩石的力学特性和变形规律。对于脆性岩石，通常采用线弹性本构关系来描述其在弹性阶段的力学行为；而对于具有一定塑性的岩石，则需要采用弹塑性本构关系来考虑其塑性变形。岩石的结构构造，如节理、裂隙等，也会对爆破损伤产生影响，在模型中可以通过引入损伤变量来考虑这些结构面对岩石力学性能的劣化作用。爆破荷载是导致岩石损伤的直接原因，炸药爆炸产生的高温高压气体和应力波会对岩石施加巨大的作用力，使岩石发生损伤和破坏。在模型中，需要合理模拟爆破荷载的作用。根据炸药的类型、装药量和装药结构等参数，计算炸药爆炸产生的能量和应力波的传播特性。采用爆炸力学理论，将炸药爆炸产生的能量等效为作用在岩石上的荷载，通过数值模拟方法，如有限元法、离散元法等，将爆破荷载施加到岩石模型上，分析岩石在爆破荷载作用下的应力、应变分布以及损伤演化过程。在有限元模拟中，可以通过在炮孔位置施加压力荷载来模拟炸药爆炸产生的应力波作用，通过调整压力荷载的大小、作用时间和分布方式，来模拟不同的爆破工况。基于损伤力学理论，引入合适的损伤变量来描述岩石在爆破过程中的损伤演化。损伤变量是一个能够反映岩石内部损伤程度的物理量，它可以是标量、矢量或张量。常见的损伤变量定义方式有基于裂纹密度、基于弹性模量劣化、基于能量耗散等。在本模型中，选择基于弹性模量劣化的损伤变量定义方式，即通过岩石弹性模量的变化来衡量岩石的损伤程度。假设岩石的初始弹性模量为E_0，随着爆破过程的进行，岩石内部产生损伤，弹性模量逐渐降低为E，则损伤变量D可以定义为D=1-\frac{E}{E_0}。通过建立损伤演化方程，描述损伤变量随应力、应变等因素的变化规律，从而实现对岩石爆破损伤过程的定量描述。损伤演化方程可以根据岩石的试验数据和理论分析来确定，它反映了岩石在爆破荷载和地应力共同作用下的损伤发展机制。3.2模型的建立过程3.2.1考虑地应力的岩石本构关系地应力的存在会显著改变岩石的力学性质，为了准确描述这种影响，需要对岩石的弹性模量、泊松比等参数进行修正，从而建立考虑地应力的岩石本构关系。在传统的岩石力学中，岩石的弹性模量E和泊松比\mu通常被视为常数，但在实际工程中，地应力的作用会使这些参数发生变化。研究表明，地应力会导致岩石内部微观结构的改变，如微裂纹的闭合或扩展，从而影响岩石的宏观力学性质。当岩石受到地应力作用时，微裂纹会在应力作用下发生闭合，使得岩石的整体刚度增加，弹性模量增大。根据相关的试验研究和理论分析，考虑地应力影响的弹性模量E'可以通过以下经验公式进行修正：E'=E(1+\alpha\sigma_{max})其中，\alpha为与岩石性质相关的修正系数，可通过试验确定；\sigma_{max}为地应力的最大值。该公式表明，弹性模量随着地应力最大值的增加而增大，且增大的幅度与修正系数\alpha有关。对于花岗岩，其\alpha值可能在0.01-0.05之间，当地应力最大值为10MPa时，若初始弹性模量E=50GPa，则修正后的弹性模量E'=50\times(1+0.03\times10)=65GPa。泊松比\mu也会受到地应力的影响。地应力会改变岩石在受力时的横向变形特性，使得泊松比发生变化。考虑地应力影响的泊松比\mu'可以通过以下公式进行修正：\mu'=\mu+\beta(\sigma_{1}-\sigma_{3})其中，\beta为泊松比修正系数，同样可通过试验确定；\sigma_{1}和\sigma_{3}分别为地应力的最大主应力和最小主应力。这意味着泊松比的变化与地应力的主应力差值有关，当主应力差值增大时，泊松比会相应地改变。在某一岩石样本中，若初始泊松比\mu=0.25，\beta=0.005，\sigma_{1}=15MPa，\sigma_{3}=5MPa，则修正后的泊松比\mu'=0.25+0.005\times(15-5)=0.3。基于修正后的弹性模量和泊松比，采用广义胡克定律建立考虑地应力的岩石本构关系。广义胡克定律描述了线性弹性材料在多轴应力状态下应力与应变之间的关系，其表达式为：\begin{cases}\varepsilon_{x}=\frac{1}{E'}[\sigma_{x}-\mu'(\sigma_{y}+\sigma_{z})]+\alpha_{T}\DeltaT\\\varepsilon_{y}=\frac{1}{E'}[\sigma_{y}-\mu'(\sigma_{x}+\sigma_{z})]+\alpha_{T}\DeltaT\\\varepsilon_{z}=\frac{1}{E'}[\sigma_{z}-\mu'(\sigma_{x}+\sigma_{y})]+\alpha_{T}\DeltaT\end{cases}其中，\varepsilon_{x}、\varepsilon_{y}、\varepsilon_{z}分别为x、y、z方向的应变；\sigma_{x}、\sigma_{y}、\sigma_{z}分别为x、y、z方向的应力；\alpha_{T}为热膨胀系数，考虑到爆破过程中可能产生的温度变化对岩石变形的影响；\DeltaT为温度变化量。在实际应用中，根据岩石所处的地应力状态和受力情况，将修正后的弹性模量E'和泊松比\mu'代入广义胡克定律，即可得到考虑地应力的岩石本构关系，从而准确描述岩石在复杂应力条件下的力学行为。3.2.2爆破荷载的施加方式在岩石爆破损伤模型中，准确施加爆破荷载是模拟爆破过程的关键环节，需要合理考虑爆炸应力波和爆生气体压力的加载方式及参数确定方法。爆炸应力波是炸药爆炸瞬间产生的高强度应力脉冲，其传播速度极快，对岩石的初始破坏起着重要作用。在数值模拟中，通常采用压力时程曲线来模拟爆炸应力波的加载过程。压力时程曲线描述了爆炸应力波随时间的变化规律，其形状和参数与炸药的类型、装药量、装药结构以及岩石的性质等因素密切相关。对于常用的乳化炸药，其爆炸应力波的压力时程曲线可以通过实验测量或理论计算获得。根据相关研究，爆炸应力波的峰值压力P_{max}可以通过以下公式估算：P_{max}=\frac{\rho_{0}D^{2}}{4}其中，\rho_{0}为炸药的初始密度，D为炸药的爆速。在实际应用中，根据炸药的具体参数，如乳化炸药的密度通常为1.1-1.3g/cm^{3}，爆速在3000-5000m/s之间，代入公式即可计算出爆炸应力波的峰值压力。假设乳化炸药的密度\rho_{0}=1.2g/cm^{3}，爆速D=4000m/s，则P_{max}=\frac{1.2\times10^{3}\times4000^{2}}{4}=4.8\times10^{9}Pa。爆炸应力波的作用时间t_{0}也是一个重要参数，它反映了爆炸应力波对岩石作用的持续时间。作用时间通常与炸药的性质和装药结构有关，一般在微秒到毫秒量级。在数值模拟中，可以根据实际情况选择合适的作用时间，如对于小型爆破试验，作用时间可能在几十微秒左右；对于大型工程爆破，作用时间可能在几百微秒到几毫秒之间。通过确定爆炸应力波的峰值压力和作用时间，结合压力时程曲线的形状（如三角形、指数衰减型等），可以在数值模型中准确施加爆炸应力波荷载。爆生气体压力是炸药爆炸后产生的高温高压气体对炮孔壁施加的持续压力，它对岩石的后续破碎和损伤扩展起着关键作用。在数值模拟中，通常采用准静态压力加载方式来模拟爆生气体压力的作用。准静态压力加载假设爆生气体在炮孔内的膨胀过程是缓慢的，其压力变化可以用静态压力来近似。爆生气体压力P_{g}可以通过以下公式计算：P_{g}=\frac{\gamma_{g}nQ}{V_{0}}其中，\gamma_{g}为爆生气体的绝热指数，一般取值在1.2-1.4之间；n为炸药的爆炸反应不完全系数，取值范围为0.9-1.0；Q为炸药的爆热，不同类型的炸药爆热不同，如乳化炸药的爆热一般在3000-4000kJ/kg之间；V_{0}为炮孔的初始容积。在实际工程中，根据炸药和炮孔的参数，代入公式即可计算出爆生气体压力。假设炮孔直径为50mm，深度为2m，则炮孔初始容积V_{0}=\pi\times(0.025)^{2}\times2=0.003927m^{3}，若使用乳化炸药，爆热Q=3500kJ/kg，爆炸反应不完全系数n=0.95，绝热指数\gamma_{g}=1.3，则爆生气体压力P_{g}=\frac{1.3\times0.95\times3500\times10^{3}}{0.003927}\approx1.1\times10^{9}Pa。在数值模拟中，将计算得到的爆生气体压力按照一定的时间步长逐渐施加到炮孔壁上，模拟爆生气体对岩石的持续作用过程。通过合理确定爆炸应力波和爆生气体压力的加载方式及参数，可以在数值模型中真实地模拟爆破荷载对岩石的作用，为研究岩石的爆破损伤提供准确的荷载条件。3.2.3损伤变量的定义与演化方程为了定量描述岩石在爆破过程中的损伤程度及其演化过程，基于岩石微裂纹扩展的原理定义损伤变量，并建立其在爆破过程中的演化方程。岩石在爆破荷载和地应力的共同作用下，内部会产生大量的微裂纹，这些微裂纹的萌生、扩展和贯通是导致岩石损伤和破坏的主要原因。基于此，定义损伤变量D来描述岩石内部微裂纹的发展程度。采用基于裂纹密度的损伤变量定义方法，即损伤变量D与岩石内部单位体积内的裂纹面积A_{c}成正比。假设岩石的初始状态为无损伤状态，此时损伤变量D=0；当岩石内部产生微裂纹后，随着裂纹的扩展和增多，损伤变量D逐渐增大，当岩石完全破坏时，损伤变量D=1。损伤变量D的具体表达式为：D=\frac{A_{c}}{A_{0}}其中，A_{0}为岩石单位体积的初始面积。在实际计算中，通过对岩石内部微裂纹的观测和统计，获取裂纹面积A_{c}，进而计算出损伤变量D。利用显微镜对岩石试件进行观察，统计单位体积内的裂纹数量和尺寸，从而计算出裂纹面积A_{c}，再结合岩石的初始面积A_{0}，得到损伤变量D的值。为了描述损伤变量D在爆破过程中的演化规律，建立损伤演化方程。损伤演化方程考虑了爆破应力、地应力以及岩石自身性质等因素对损伤发展的影响。根据相关的损伤力学理论和试验研究，建立如下的损伤演化方程：\frac{dD}{dt}=f(\sigma,\varepsilon,D)其中，\frac{dD}{dt}表示损伤变量D随时间的变化率，即损伤演化速率；\sigma为岩石所受的应力，包括爆破应力和地应力；\varepsilon为岩石的应变；f(\sigma,\varepsilon,D)为损伤演化函数，它反映了损伤演化速率与应力、应变和损伤变量之间的关系。损伤演化函数f(\sigma,\varepsilon,D)可以通过理论推导和试验数据拟合得到。根据试验结果，发现当岩石所受的应力超过其损伤阈值\sigma_{th}时，损伤演化速率会显著增大。损伤演化函数可以表示为：f(\sigma,\varepsilon,D)=\begin{cases}0,&\sigma\leq\sigma_{th}\\k_{1}(\sigma-\sigma_{th})^{\alpha}(1-D)^{\beta}\varepsilon^{\gamma},&\sigma\gt\sigma_{th}\end{cases}其中，k_{1}为损伤演化系数，与岩石性质和爆破条件有关；\alpha、\beta、\gamma为与损伤演化相关的参数，可通过试验确定。在某一岩石爆破试验中，通过对不同应力水平下岩石损伤演化的观测和分析，确定k_{1}=0.01，\alpha=1.5，\beta=0.5，\gamma=1。当岩石所受应力\sigma=20MPa，损伤阈值\sigma_{th}=15MPa，损伤变量D=0.3，应变\varepsilon=0.001时，损伤演化速率\frac{dD}{dt}=0.01\times(20-15)^{1.5}\times(1-0.3)^{0.5}\times0.001^{1}\approx1.3\times10^{-5}。通过上述损伤变量的定义和演化方程，能够定量地描述岩石在爆破过程中的损伤发展过程，为分析岩石的爆破损伤机制和预测爆破效果提供了重要的理论依据。在数值模拟中，将损伤演化方程与岩石的本构关系相结合，通过迭代计算，可以实时更新损伤变量D的值，从而模拟岩石在爆破过程中损伤的动态演化过程。3.3模型验证与分析3.3.1与已有理论模型对比为了验证所建立的地应力影响下岩石爆破损伤模型的合理性，将本文模型的计算结果与其他经典的岩石爆破损伤模型进行对比分析。选择了两种在岩石爆破损伤研究领域具有代表性的经典模型，分别是Grady模型和Kuszmaul模型。Grady模型是基于断裂力学理论建立的，它假设岩石在爆破过程中的损伤主要是由于裂纹的扩展和贯通导致的。该模型通过引入裂纹扩展速度和断裂能等参数来描述岩石的损伤演化过程。在Grady模型中，裂纹扩展速度与爆炸应力波的强度和岩石的断裂韧性有关，通过求解裂纹扩展的动力学方程来确定裂纹的扩展路径和长度，进而计算岩石的损伤程度。Kuszmaul模型则是从能量的角度出发，认为岩石爆破损伤过程是能量的输入、耗散和转化的过程。该模型考虑了爆炸能量、岩石的弹性应变能以及损伤耗散能之间的关系，通过建立能量平衡方程来描述岩石的损伤演化。在Kuszmaul模型中，损伤耗散能与岩石内部微裂纹的产生和扩展密切相关，通过计算损伤耗散能的变化来确定岩石的损伤变量。在对比分析时，设定了相同的岩石材料参数、爆破参数以及地应力条件。岩石材料选用花岗岩，其基本力学参数为：抗压强度150MPa，抗拉强度10MPa，弹性模量60GPa，泊松比0.25。爆破参数为：炸药采用乳化炸药，装药量为1kg，炮孔直径50mm，炮孔深度2m。地应力条件设定为：最大主应力\sigma_{1}=15MPa，最小主应力\sigma_{3}=5MPa。分别利用本文模型、Grady模型和Kuszmaul模型对岩石爆破损伤进行计算，得到了岩石在爆破后的损伤区域分布和损伤程度。从损伤区域分布来看，本文模型计算得到的损伤区域形状和范围与其他两种模型存在一定差异。本文模型考虑了地应力对岩石力学性质的影响以及地应力与爆破应力波的相互作用，使得损伤区域在主应力方向上有一定的扩展，且损伤区域的边界更加不规则。而Grady模型和Kuszmaul模型由于没有充分考虑地应力的复杂影响，损伤区域相对较为规则，且在主应力方向上的扩展不明显。在损伤程度方面，通过对比三种模型计算得到的损伤变量值，发现本文模型计算得到的损伤变量在炮孔周围区域相对较大，随着距离炮孔距离的增加，损伤变量逐渐减小，且减小的速率与其他两种模型不同。这是因为本文模型中损伤演化方程考虑了地应力和爆破应力的联合作用，以及岩石内部微裂纹在不同应力状态下的扩展特性，使得损伤程度的分布更加符合实际情况。而Grady模型主要侧重于裂纹扩展速度的计算，对其他因素的考虑相对较少；Kuszmaul模型虽然从能量角度出发，但在处理地应力对能量转化的影响时不够全面，导致损伤程度的计算结果与本文模型存在差异。综合对比结果表明，本文建立的地应力影响下岩石爆破损伤模型在考虑地应力的复杂性和岩石材料的非线性特性方面具有一定优势，计算结果更能反映实际岩石爆破损伤的特征，验证了模型的合理性和有效性。然而，由于岩石爆破损伤过程的复杂性，任何模型都存在一定的局限性，未来还需要进一步结合实际工程数据和试验研究，对模型进行优化和完善。3.3.2敏感性分析为了深入了解地应力大小、方向等因素对损伤变量的影响，确定关键参数，开展敏感性分析。敏感性分析是通过改变模型中某个参数的值，观察模型输出结果的变化情况，从而评估该参数对模型结果的影响程度。首先分析地应力大小对损伤变量的影响。保持其他参数不变，仅改变地应力的大小。当地应力较小时，损伤变量的增长相对较为缓慢，损伤区域主要集中在炮孔周围较小的范围内。这是因为较小的地应力对岩石的初始力学性质影响较小，爆破应力波在传播过程中受到的干扰也较小，岩石主要在爆破荷载的直接作用下产生损伤。随着地应力的逐渐增大，损伤变量的增长速度明显加快，损伤区域也显著扩大。这是因为地应力的增大改变了岩石的力学性质，使岩石的抗压强度和弹性模量发生变化，同时地应力与爆破应力波相互作用，增强了岩石内部的应力集中效应，导致裂纹更容易起裂和扩展，从而使损伤程度加剧。当最大主应力从10MPa增加到20MPa时，损伤变量在炮孔周围一定范围内增加了约30\%，损伤区域的半径也增大了约20\%。接着分析地应力方向对损伤变量的影响。设定不同的地应力方向，观察损伤变量的变化。当爆破方向与地应力主方向一致时，损伤变量在该方向上的增长最为明显，损伤区域呈现出沿地应力主方向的长条状扩展。这是因为在这种情况下，地应力对爆破应力波起到了增强作用，使得岩石更容易沿着该方向产生裂纹和破碎。当爆破方向与地应力主方向垂直时，损伤变量的增长相对较慢，损伤区域在垂直方向上的扩展受到一定限制，呈现出相对较窄的分布。这是因为地应力对垂直方向上的裂纹扩展产生了阻碍作用，改变了裂纹的扩展路径，使得损伤程度相对较小。在某一数值模拟中，当爆破方向与地应力主方向夹角为0^{\circ}时，损伤区域在该方向上的长度是夹角为90^{\circ}时的1.5倍。通过敏感性分析，确定地应力的大小和方向是影响岩石爆破损伤的关键参数。在实际工程中，应充分考虑地应力的这些因素，根据具体的地应力条件优化爆破参数，以达到更好的爆破效果，减少岩石的过度损伤和浪费，同时提高工程的安全性和经济性。在高地应力区域进行爆破作业时，应适当调整炮孔的布置和起爆顺序，使爆破方向尽量与地应力主方向相适应，以充分利用地应力的作用，提高爆破效率；在低地应力区域，则可以适当降低炸药的用量，避免岩石的过度破碎。四、基于构建模型的数值试验设计与实施4.1数值试验方案设计4.1.1试验模型的建立利用专业的数值模拟软件ANSYS/LS-DYNA建立岩石爆破模型。在建立模型时，充分考虑实际工程中的各种因素，确保模型能够真实地反映岩石爆破的实际情况。首先确定模型的几何尺寸。根据常见的岩石爆破工程场景，将模型设定为长方体形状，其长、宽、高分别为5m、3m、2m。这样的尺寸既能满足模拟岩石爆破过程中应力波传播和岩石损伤的需求，又能在一定程度上减少计算量，提高计算效率。在实际工程中，如隧道爆破，隧道的截面尺寸和长度会影响爆破应力波的传播和反射，合理的模型尺寸可以更准确地模拟这种影响。对于一些小型的岩石爆破试验，模型尺寸可能会相对较小，但也需要根据试验目的和实际情况进行合理选择。接着对模型进行网格划分，这是数值模拟中至关重要的一步，网格的质量和密度直接影响模拟结果的准确性和计算效率。采用六面体单元对模型进行网格划分，在炮孔周围区域，由于爆破应力波的作用较为集中，岩石的损伤和破坏较为严重，因此对该区域进行加密处理，将单元尺寸设置为0.05m，以更精确地捕捉岩石在爆破过程中的力学响应。在远离炮孔的区域，岩石的损伤相对较小，为了减少计算量，将单元尺寸逐渐增大至0.2m。通过这种渐变的网格划分方式，既能保证在关键区域的计算精度，又能有效地控制计算成本。在进行网格划分时，还需要注意单元的形状和质量，避免出现畸形单元，以确保计算结果的可靠性。明确模型的材料参数，这是描述岩石和炸药物理力学性质的关键。对于岩石材料，通过大量的岩石力学试验和相关文献调研，获取其基本力学参数。其密度设定为2500kg/m³，这是常见岩石密度的典型值，反映了岩石的质量分布特性。弹性模量为40GPa，弹性模量是衡量岩石抵抗弹性变形能力的重要指标，该值表示岩石在弹性阶段受力时的变形难易程度。泊松比为0.2，泊松比描述了岩石在受力时横向应变与纵向应变的比值，反映了岩石的横向变形特性。抗压强度为120MPa，抗拉强度为8MPa，抗压强度和抗拉强度分别表示岩石抵抗压缩和拉伸破坏的能力，这两个参数对于分析岩石在爆破过程中的破坏模式至关重要。对于炸药材料，选用常用的乳化炸药，其密度为1100kg/m³，爆速为4000m/s，爆速是炸药爆炸时爆轰波传播的速度，它直接影响爆炸应力波的强度和传播特性。通过准确设定这些材料参数，能够使模型更真实地模拟岩石爆破的物理过程。4.1.2试验参数的选取与控制为了全面研究地应力和爆破参数对岩石爆破损伤的影响规律，精心选择不同的地应力水平、爆破参数等，并设置多组对比试验。在不同地应力水平设置方面，考虑到实际工程中地应力的变化范围，设置了三种典型的地应力工况。工况一为低地应力水平，最大主应力\sigma_{1}设置为5MPa，最小主应力\sigma_{3}设置为3MPa，侧压系数\lambda（即最小主应力与最大主应力的比值）为0.6。在一些浅埋的岩石工程中，地应力水平相对较低，这种工况可以模拟此类工程场景下的岩石爆破情况。工况二为中地应力水平，\sigma_{1}设置为10MPa，\sigma_{3}设置为6MPa，侧压系数\lambda为0.6，这是一种较为常见的地应力工况，适用于许多一般地质条件下的岩石爆破工程。工况三为高地应力水平，\sigma_{1}设置为15MPa，\sigma_{3}设置为9MPa，侧压系数\lambda同样为0.6，在深部岩石工程中，如深部矿山开采、深埋隧道挖掘等，地应力水平较高，这种工况可以模拟这些工程中的高地应力环境。通过设置这三种不同的地应力水平，能够系统地研究地应力大小对岩石爆破损伤的影响。在爆破参数选择上，装药量的大小直接影响爆破能量的输入，进而影响岩石的损伤程度。设置了三种装药量，分别为0.5kg、1.0kg和1.5kg。较小的装药量可以模拟对岩石破碎程度要求较低的爆破工况，如在一些对周边环境要求较高的岩石开挖工程中，需要控制装药量以减少对周围岩体的影响；较大的装药量则可以模拟需要大量破碎岩石的工程场景，如大型矿山的开采。炮孔间距是影响爆破效果的重要参数之一，它决定了相邻炮孔爆炸应力波的相互作用程度。设置了三种炮孔间距，分别为0.8m、1.0m和1.2m。较小的炮孔间距会使爆炸应力波相互叠加作用增强，可能导致岩石过度破碎；较大的炮孔间距则可能使岩石破碎不均匀，出现大块岩石。通过改变炮孔间距，可以研究其对岩石破碎效果和损伤范围的影响。起爆顺序的设计对于控制岩石的破裂过程和损伤范围至关重要。设置了三种起爆顺序，分别为顺序起爆、中心起爆和对角起爆。顺序起爆是按照一定的顺序依次起爆各个炮孔，这种起爆顺序可以使爆破应力波依次传播，有利于控制岩石的破裂方向；中心起爆是先起爆中心炮孔，然后依次起爆周围炮孔，这种起爆方式可以使爆炸能量集中在中心区域，然后向四周扩散；对角起爆是对角线上的炮孔先起爆，然后依次起爆其他炮孔，这种起爆顺序可以使爆破应力波在岩石中形成复杂的叠加效果，从而影响岩石的破碎模式。通过对比不同起爆顺序下的岩石爆破损伤情况，可以优化起爆顺序，提高爆破效率。4.2数值试验的实施过程4.2.1模拟软件的选择与操作选用ANSYS/LS-DYNA软件进行数值模拟。ANSYS/LS-DYNA是一款功能强大的非线性动力分析软件，在岩石爆破模拟领域具有广泛的应用。它能够精确模拟岩石在爆破过程中的大变形、材料非线性以及接触碰撞等复杂力学行为，为研究地应力影响下的岩石爆破损伤提供了有力的工具。在操作过程中，首先将建立好的岩石爆破模型导入ANSYS/LS-DYNA软件中。在导入模型时，确保模型的几何形状、网格划分以及材料参数等信息的准确性，避免因导入错误而影响模拟结果。对模型进行边界条件的设置，根据实际工程情况，将模型的底部边界设置为固定约束，限制其在三个方向上的位移，以模拟岩石与基础的固定连接；将模型的四周边界设置为无反射边界条件，确保爆破应力波在传播到边界时能够自由传播，避免应力波的反射对模拟结果产生干扰。这种边界条件的设置能够更真实地模拟岩石在实际爆破过程中的受力状态。在软件中设置地应力和爆破荷载。对于地应力的施加，根据预先设定的地应力水平，在软件中通过定义初始应力场的方式，将不同大小和方向的地应力准确地施加到岩石模型上。在设置爆破荷载时，按照前面确定的爆炸应力波和爆生气体压力的加载方式及参数，在炮孔位置施加相应的压力荷载。在施加爆炸应力波荷载时，根据压力时程曲线，设置荷载的峰值、作用时间和加载方式，确保能够准确模拟爆炸应力波的瞬间冲击作用；在施加爆生气体压力荷载时，按照准静态压力加载方式，将计算得到的爆生气体压力按照一定的时间步长逐渐施加到炮孔壁上，模拟爆生气体对岩石的持续作用过程。在施加荷载时，仔细检查参数的设置，确保荷载的施加符合实际爆破情况。设置求解控制参数，如时间步长、求解算法等。时间步长的选择至关重要，它直接影响计算的精度和效率。时间步长过小，会导致计算量急剧增加，计算时间过长；时间步长过大，则可能会导致计算结果不准确，无法捕捉到爆破过程中的一些关键信息。根据模型的特点和计算要求，通过多次试验和验证，选择合适的时间步长，一般在微秒量级。在选择求解算法时，根据岩石爆破问题的特点，选择显式动力学算法，该算法适用于求解高速冲击和大变形等动力学问题，能够有效地模拟爆破过程中岩石的动态响应。在设置求解控制参数后，进行试算，观察计算结果的收敛性和合理性，根据试算结果对参数进行调整和优化，确保模拟计算的准确性和稳定性。4.2.2模拟结果的记录与初步分析在模拟计算完成后，利用ANSYS/LS-DYNA软件的后处理功能，详细记录岩石在爆破过程中的应力、应变、损伤分布等模拟结果。通过软件的云图显示功能，生成不同时刻岩石的应力云图、应变云图和损伤云图，直观地展示岩石内部的力学响应和损伤演化过程。从应力云图中可以清晰地看到，在爆破瞬间，炮孔周围岩石的应力急剧升高，形成明显的应力集中区域。随着爆炸应力波的传播，应力逐渐向周围扩散，应力值也逐渐减小。在高地应力条件下，应力分布呈现出明显的各向异性，地应力方向上的应力值相对较高，这是因为地应力的存在改变了岩石内部的初始应力状态，使得爆破应力波在传播过程中受到地应力的影响，从而导致应力分布的不均匀。在最大主应力方向上，应力集中区域的范围更大，应力值也更高，这表明地应力对该方向上的应力分布影响更为显著。通过对应力云图的分析，可以了解爆破应力波在岩石中的传播规律以及地应力对其的影响，为进一步分析岩石的损伤机制提供依据。应变云图展示了岩石在爆破过程中的变形情况。在炮孔周围，岩石的应变较大，这是由于爆炸应力波的强烈作用导致岩石发生了较大的变形。随着距离炮孔距离的增加，应变逐渐减小。在不同地应力水平下，岩石的应变分布也存在差异。在低地应力水平下，岩石的应变相对较小，变形主要集中在炮孔周围较小的范围内；而在高地应力水平下，岩石的应变范围明显扩大，且在主应力方向上的应变更为显著。这是因为高地应力使得岩石的力学性质发生改变，岩石在受到爆破荷载作用时更容易发生变形，且变形在主应力方向上得到了进一步的增强。通过分析应变云图，可以了解岩石在爆破过程中的变形特征以及地应力对变形的影响，为评估岩石的稳定性提供参考。损伤云图则直观地显示了岩石在爆破后的损伤区域分布。在炮孔周围，损伤程度最为严重，损伤区域呈现出以炮孔为中心的近似圆形分布。随着距离炮孔距离的增加，损伤程度逐渐减小。在不同地应力条件下，损伤区域的形状和范围也有所不同。在高地应力条件下，损伤区域在主应力方向上有明显的扩展，形状更加不规则，这是由于地应力与爆破应力波相互作用，使得裂纹更容易沿着主应力方向扩展，从而导致损伤区域的扩大和形状的改变。通过对损伤云图的分析，可以直观地了解岩石的损伤范围和程度，以及地应力对损伤分布的影响，为评估爆破效果和岩石的剩余强度提供依据。对记录的模拟结果进行初步的可视化分析，对比不同地应力水平和爆破参数下的模拟结果，观察岩石的应力、应变和损伤分布的变化规律。通过这种初步分析，可以发现地应力和爆破参数对岩石爆破损伤有着显著的影响，为后续更深入的定量分析奠定基础。在对比不同装药量下的模拟结果时，可以发现随着装药量的增加，炮孔周围的应力、应变和损伤程度都明显增大，损伤区域的范围也随之扩大。在不同地应力水平下，相同爆破参数的模拟结果也存在明显差异，高地应力条件下岩石的损伤程度和范围都大于低地应力条件下的情况，这进一步验证了地应力对岩石爆破损伤的重要影响。4.3数值试验结果分析4.3.1地应力对岩石爆破损伤范围的影响通过对不同地应力水平下的数值试验结果进行深入分析，清晰地揭示了地应力对岩石爆破损伤范围的显著影响规律。在低地应力水平下，岩石爆破后的损伤范围相对较小，损伤区域主要集中在炮孔周围。这是因为低地应力对岩石的初始力学性质改变较小，岩石在爆破应力波和爆生气体的作用下，其抵抗破坏的能力相对较强，裂纹的扩展受到一定限制，难以向远处延伸。在低地应力工况下，损伤区域的半径大约在1.0-1.2m之间。随着地应力水平的升高，岩石爆破损伤范围呈现出明显的增大趋势。在中地应力水平下，损伤区域的范围明显扩大，不仅在炮孔周围的损伤程度加剧，而且损伤区域向远处的扩展也更为明显。这是由于地应力的增大改变了岩石的力学性质，使得岩石的弹性模量和抗压强度发生变化，同时地应力与爆破应力波相互作用，增强了岩石内部的应力集中效应，促使裂纹更容易起裂和扩展，从而导致损伤范围的扩大。在中地应力工况下，损伤区域的半径增大到1.5-1.8m左右。在高地应力水平下，岩石爆破损伤范围进一步扩大，损伤区域的边界更加不规则，且在主应力方向上的扩展尤为显著。高地应力使得岩石内部的应力状态更加复杂，爆破应力波在传播过程中受到地应力的影响更大，裂纹更容易沿着主应力方向扩展，形成更长、更宽的裂纹，从而导致损伤范围在主应力方向上的延伸。高地应力还会使岩石的脆性增加，更容易发生破碎，进一步扩大了损伤范围。在高地应力工况下，损伤区域在主应力方向上的半径可达2.0-2.5m，且损伤区域的形状呈现出明显的各向异性。地应力对岩石爆破损伤范围的影响主要是通过改变岩石的力学性质以及与爆破应力波的相互作用来实现的。地应力的存在使得岩石内部的微观结构发生变化，如微裂纹的闭合或扩展，从而改变了岩石的宏观力学性能。地应力还会影响爆破应力波的传播速度、方向和衰减特性，进而影响裂纹的起裂和扩展，最终导致岩石爆破损伤范围的改变。在高地应力条件下，地应力会使爆破应力波在传播过程中发生折射和反射，使得应力波的能量分布更加不均匀，从而加剧了岩石的损伤和破坏，扩大了损伤范围。4.3.2地应力对岩石爆破裂纹扩展的影响地应力对岩石爆破裂纹扩展的影响体现在多个方面，包括裂纹扩展方向、长度和数量等。在不同地应力条件下，裂纹扩展呈现出明显不同的特征。在低地应力条件下，裂纹扩展方向相对较为随机，没有明显的优势方向。这是因为低地应力对岩石内部应力场的影响较小，爆破应力波在传播过程中受到的干扰较弱，裂纹主要在爆破应力波和爆生气体的直接作用下产生和扩展，其扩展方向主要取决于炮孔周围的应力分布和岩石的初始缺陷。在低地应力工况下，裂纹扩展方向呈现出以炮孔为中心的近似放射状分布，各方向上的裂纹扩展长度相对较为均匀。随着地应力水平的增加，裂纹扩展方向逐渐呈现出一定的规律性。在中地应力条件下，裂纹扩展方向开始受到地应力的影响，部分裂纹倾向于沿着最大主应力方向扩展。这是因为地应力的存在改变了岩石内部的应力状态，使得沿着最大主应力方向的岩石更容易达到其破坏强度，从而促进了裂纹在该方向上的扩展。在中地应力工况下，沿着最大主应力方向的裂纹扩展长度明显大于其他方向，裂纹扩展路径也更加规则，呈现出一定的方向性。在高地应力条件下，裂纹扩展方向主要沿着最大主应力方向，形成明显的优势扩展方向。高地应力使得岩石内部的应力差增大，沿着最大主应力方向的岩石受到的拉应力作用更强，更容易发生破裂，从而导致裂纹主要沿着该方向扩展。高地应力还会使岩石的各向异性增强，进一步促进裂纹在最大主应力方向上的扩展。在高地应力工况下，沿着最大主应力方向的裂纹扩展长度远远大于其他方向，裂纹形成了较为集中的破裂带，且破裂带的宽度也较大。地应力对裂纹长度和数量也有显著影响。随着地应力水平的升高，裂纹长度明显增加。这是因为地应力的增大为裂纹扩展提供了更大的驱动力，使得裂纹能够克服岩石的阻力向远处延伸。高地应力还会使岩石的脆性增加，裂纹更容易扩展，从而导致裂纹长度的增大。在低地应力工况下，裂纹的平均长度较短，大约在0.2-0.3m之间；而在高地应力工况下，裂纹的平均长度可增加到0.5-0.8m。地应力水平的升高还会导致裂纹数量的增加。这是因为地应力的作用使得岩石内部的应力分布更加不均匀，产生了更多的应力集中点，从而促进了更多裂纹的起裂和扩展。在高地应力条件下，岩石内部的微裂纹更容易被激活和扩展，导致裂纹数量明显增多。4.3.3地应力与其他因素的耦合作用分析地应力与炸药类型、装药量等因素之间存在着复杂的耦合作用，这些耦合作用对爆破效果产生了显著的影响。不同炸药类型具有不同的爆炸性能，如爆速、爆热、爆轰压力等，这些性能参数会与地应力相互作用，从而影响爆破效果。高爆速炸药能够产生更强烈的爆炸应力波，在低地应力条件下，高爆速炸药的爆炸应力波能够迅速在岩石中传播，使岩石产生大量的裂纹和破碎，爆破效果较为明显。但在高地应力条件下，高爆