冲击载荷下煤的动态力学特性与破坏机制研究

冲击载荷下煤的动态力学特性与破坏机制研究一、引言1.1研究背景与意义煤炭作为一种重要的能源资源，在全球能源结构中占据着举足轻重的地位。我国是煤炭生产和消费大国，煤炭在一次能源生产和消费结构中均占较大比重。然而，在煤炭开采过程中，冲击地压等动力灾害严重威胁着煤矿的安全生产和人员生命安全。冲击地压是指煤矿井巷或工作面周围煤（岩）体由于弹性变形能的瞬时释放而产生的突然、剧烈破坏的动力现象，常伴有煤（岩）体瞬间位移、抛出、巨响及气浪等。这种灾害具有突发性、瞬时震动性及破坏性等特点，发生时一般无明显宏观前兆，以突然、急剧、猛烈的形式释放煤岩体变形能，煤岩体被抛出，造成支护体损坏、片帮冒顶、巷道堵塞、伤及人员，并伴随巨大声响和岩体震动，震动时间从几秒到几十秒，冲出的煤岩从几吨到几百吨，甚至上千吨，往往给生产与生命财产造成严重威胁。例如，2023年11月28日，黑龙江龙煤双鸭山矿业有限责任公司双阳煤矿发生冲击地压顶板事故，给企业和社会带来了巨大的损失。煤的动态力学特性是影响冲击地压等灾害发生的关键因素之一。在冲击载荷作用下，煤的力学行为与静态加载条件下有很大的不同，其强度、变形特性、能量耗散等都会发生显著变化。研究煤在冲击载荷下的动态力学特性，对于揭示冲击地压等灾害的发生机理，制定有效的防治措施具有重要的理论意义。同时，通过深入了解煤的动态力学特性，可以为煤矿开采设计、支护结构选型、开采工艺优化等提供科学依据，从而提高煤矿开采的安全性和效率，降低生产成本，具有重要的工程实际意义。1.2国内外研究现状国内外众多学者对煤的动态力学特性开展了大量的研究工作。在实验研究方面，分离式霍普金森压杆（SHPB）试验系统被广泛应用于煤的动态力学性能测试。张文清等人利用SHPB试验装置对型煤与原煤进行了不同应变率条件下的单轴冲击破坏试验，获得了两种煤的动态应力应变曲线和相关动态力学参数，并分析得到两种煤的破坏变形特征及能量耗散规律，发现两种煤在冲击荷载下的变形在初始阶段即表现出非线性，在试验应变率70-280s-1内，动态抗压强度和峰值应变均具有应变率强化的特性，随应变率呈指数增加，且原煤的应变率硬化特性更为显著。付玉凯等人运用直径为50mm的SHPB试验系统，在不同的冲击速率下对无烟煤和砂岩进行冲击加载实验，研究了它们在动态荷载下的破坏机理，包括动态应力应变、破坏形态、破碎耗散能和块度分形等方面的差异，结果显示无烟煤在初始应力应变阶段就表现出强烈的非线性特征，且冲击速率越低，这种非线性表现越明显，随着冲击速率的增加，两者都显示出塑性变形和极限强度的增强，但砂岩的塑性变形程度明显小于无烟煤。在理论分析方面，一些学者建立了煤的动态本构模型来描述其在冲击载荷下的力学行为。穆朝民等构建了煤体的动态本构方程，发现煤体在冲击荷载下呈现快速跃升式和慢速阶梯式两种损伤形式。同时，也有学者从能量理论、强度理论等角度对煤在冲击载荷下的破坏机理进行了探讨。然而，当前研究仍存在一些不足与空白。一方面，研究主要集中在单一煤种或特定条件下的煤样，对于不同煤种、不同地质条件下煤的动态力学特性对比研究较少。另一方面，虽然对煤的动态力学参数进行了大量测试，但对于这些参数的影响因素及内在联系的研究还不够深入，缺乏系统性和全面性。此外，在实验研究中，部分实验设备和方法还存在一定的局限性，导致实验结果的准确性和可靠性有待提高。1.3研究内容与方法本研究将从以下几个方面开展对煤在冲击载荷下动态力学特性的研究：煤样制备：选取不同煤种、不同产地的煤样，通过现场采集和实验室加工，制备出符合实验要求的标准煤样，包括原煤试件和型煤试件。对于原煤试件，将采煤工作面刚暴露的大块煤体现场取出，用密封袋包裹后抽真空，运抵实验室后，利用岩石钻机钻取直径合适的圆柱型煤体，再打磨切割成所需的动力学测试试件；对于型煤试件，从软分层提取松软煤，将其破碎，筛取一定粒径范围的煤粉装入自制的压煤罐，在大型刚性压力机上用特定压力压制成型煤试件。实验研究：采用分离式霍普金森压杆（SHPB）试验系统，对制备好的煤样进行不同应变率下的冲击压缩试验。通过改变子弹的速度来调节冲击加载的应变率，获取煤样在冲击载荷作用下的应力-应变曲线、应变率-应力曲线等动态力学响应数据。同时，利用高速摄像机记录煤样在冲击过程中的破坏形态和变形过程，为后续的分析提供直观依据。理论分析：基于实验结果，结合岩石动力学、材料力学等相关理论，分析煤在冲击载荷下的变形特征、强度特性、能量耗散规律等。探讨应变率、煤种、煤体结构等因素对煤动态力学特性的影响机制，建立煤在冲击载荷下的动态力学模型，揭示其破坏机理。数据分析与处理：运用统计学方法和数据拟合技术，对实验数据进行分析和处理，得到煤的动态力学参数与应变率等因素之间的定量关系。利用Origin、MATLAB等软件绘制图表，直观展示煤的动态力学特性变化规律。本研究采用的实验设备主要为分离式霍普金森压杆（SHPB）试验系统，该系统能够实现对材料在高应变率下的动态加载，具有加载速率可控、测试精度高等优点。在分析方法上，综合运用理论分析、实验研究和数值模拟相结合的手段，全面深入地研究煤在冲击载荷下的动态力学特性，为煤矿安全生产和冲击地压防治提供理论支持和技术指导。二、实验材料与方法2.1煤样采集与制备本次实验选取了[某煤矿名称]的煤样，该煤矿位于[具体地理位置]，开采的煤层为[煤层编号]。在采样过程中，为了确保采集的煤样能够代表该煤层的整体特性，我们严格遵循相关采样标准和规范。采样地点选择在采煤工作面的不同位置，避开了断层、褶皱等地质构造复杂区域以及煤柱破碎带等可能导致煤样性质异常的地方。使用专业的采样钻，其具备足够的硬度和耐磨性，以保证能够顺利采集到煤样。在采样钻进入煤层前，仔细清除钻孔周围的杂物，并确保钻孔的直径和深度符合要求。采样钻进入煤层后，以匀速、稳定的方式进行旋转，避免因过快或过慢造成采样不均匀。采集完成后，将采样桶中的煤样倒入密封的采样袋中，避免煤样受到外界空气、水分及其他杂质的污染。对每个采样点的采样过程进行详细记录，包括采样地点、采样时间、采样深度、采样工具使用情况等信息，以确保采样数据的准确性和可追溯性。将采集到的煤样运输至实验室后，进行进一步的加工处理以制成标准试件。对于原煤试件，首先将大块煤体用岩石切割机切割成大致的圆柱形状，然后利用高精度的磨床对其两端面进行打磨，使其平整度和垂直度满足实验要求。最终加工成直径为[X]mm，高度为[X]mm的标准圆柱试件，根据相关标准，试件两端面的不平行度不超过0.05mm，垂直度不超过0.25°。对于型煤试件，从软分层提取松软煤，将其破碎后筛取粒径在[具体粒径范围]的煤粉，按照一定的配比加入粘结剂（如地层水或标准盐水），充分搅拌均匀。将混合均匀的样品分多次放入自制的煤心模具中，每加入一次样品进行一次压实。例如，采用5MPa的压力进行预压制，每次预压制时间为2min。完成装填和预压制后，在2min内匀速加压至35MPa，并保持10min，以确保型煤试件的成型质量和密度均匀性。压制完成后，打开煤心模具，将制好的型煤试件取出，两端放上导流垫块进行封装。2.2实验设备与原理2.2.1分离式霍普金森压杆（SHPB）系统本实验采用的分离式霍普金森压杆（SHPB）系统主要由动力装置、子弹、入射杆、透射杆、吸收杆、阻尼器、应变片以及数据采集仪等部分组成。动力装置用于发射子弹，为整个实验提供冲击加载的动力，通常采用气体驱动的方式，如高压空气或氮气，通过调节气体压力来控制子弹的发射速度，从而实现不同应变率下的冲击加载。子弹在高速运动下撞击入射杆，在入射杆中产生应力波。入射杆和透射杆一般采用高强度的金属材料制成，如合金钢，其弹性模量和波阻抗等参数经过精心设计和测试，以保证应力波在杆中的稳定传播和准确测量。在入射杆和透射杆的表面粘贴高精度的应变片，用于测量应力波在杆中传播时引起的应变变化。SHPB系统的工作原理基于一维应力波理论。当子弹撞击入射杆时，在入射杆中产生一个入射应力波\varepsilon_{i}(t)，该应力波以弹性波的速度C_{0}在杆中传播。当入射波到达试件与入射杆的界面时，一部分应力波会透射入试件，使试件发生变形，同时另一部分应力波会被反射回入射杆，形成反射应力波\varepsilon_{r}(t)。经过试件变形后的应力波继续传播进入透射杆，形成透射应力波\varepsilon_{t}(t)。根据一维应力波理论和力平衡条件，可以通过测量得到的入射波、反射波和透射波的应变信号，计算出作用在试件上的轴向应力\sigma(t)、轴向应变率\dot{\varepsilon}(t)和应变历史\varepsilon(t)。其计算公式如下（三波法公式）：\sigma(t)=\frac{EA}{A_{s}}(\varepsilon_{i}(t)+\varepsilon_{r}(t)+\varepsilon_{t}(t))\dot{\varepsilon}(t)=\frac{C_{0}}{L_{s}}(\varepsilon_{i}(t)-\varepsilon_{r}(t)-\varepsilon_{t}(t))\varepsilon(t)=\frac{C_{0}}{L_{s}}\int_{0}^{t}(\varepsilon_{i}(t)-\varepsilon_{r}(t)-\varepsilon_{t}(t))dt其中，E为压杆的弹性模量，A为压杆的横截面积，A_{s}为试件的横截面积，L_{s}为试件的长度。SHPB系统用于煤动态力学特性测试具有多方面的优势。首先，它能够实现高应变率加载，应变率范围通常可达到10^{2}-10^{4}s^{-1}，这与煤在实际冲击地压等灾害发生时所承受的应变率条件相近，能够真实模拟煤在冲击载荷下的力学响应。其次，SHPB系统的测试原理基于成熟的一维应力波理论，通过对应力波的测量和分析，可以较为准确地获取煤样在冲击过程中的动态力学参数。再者，该系统结构相对简单，操作方便，并且具有较高的可靠性和重复性，能够保证实验结果的准确性和稳定性，便于对不同煤样和不同实验条件下的煤动态力学特性进行系统研究。2.2.2其他辅助设备在实验过程中，除了SHPB系统外，还使用了多种辅助设备。应变片作为关键的测量元件，被粘贴在入射杆和透射杆的表面，用于测量应力波传播过程中杆表面的应变。选用的应变片具有高精度、高灵敏度和良好的稳定性，其电阻值和灵敏系数经过精确校准，能够准确地将应变信号转换为电信号。数据采集仪用于采集应变片输出的电信号，并将其转换为数字信号进行存储和后续处理。本实验采用的高速数据采集仪具有高采样频率，能够满足在冲击加载过程中快速变化的应力波信号采集需求，确保不会丢失重要的数据信息。例如，采样频率可设置为1MHz以上，以保证能够精确捕捉到应力波的瞬态变化过程。此外，还配备了高精度的游标卡尺和电子天平，用于测量煤样试件的尺寸和质量。游标卡尺的精度可达到0.02mm，能够准确测量煤样的直径和高度，为后续计算煤样的横截面积和体积提供准确的数据。电子天平的精度为0.001g，用于测量煤样的质量，以便根据质量和体积计算煤样的密度，为实验数据分析提供必要的参数。同时，为了记录煤样在冲击过程中的破坏形态和变形过程，使用了高速摄像机。高速摄像机能够以高帧率拍摄煤样的冲击过程，帧率可达到10000帧/秒以上，通过对拍摄的图像序列进行分析，可以直观地观察煤样的破坏模式、裂纹扩展路径等信息，为深入研究煤的动态破坏机理提供重要的依据。2.3实验方案设计2.3.1变量控制本实验主要考虑冲击速度、含水率、层理等因素对煤动态力学特性的影响，并通过控制变量法进行实验研究。在冲击速度方面，通过调节SHPB系统中动力装置的气体压力来改变子弹的发射速度，从而实现不同冲击速度下的加载。设置多个冲击速度水平，如3m/s、5m/s、7m/s等，每个速度水平下进行多次重复实验，以保证实验结果的可靠性和统计学意义。对于含水率变量，将制备好的煤样分为不同的组，一组保持原煤样的自然含水率，其他组通过喷雾加湿或烘干等方法调整煤样的含水率至设定值，如5%、10%、15%等。在调整含水率后，将煤样密封保存一定时间，使水分在煤样内部充分均匀分布，然后再进行冲击实验，以研究含水率对煤动态力学性能的影响。层理是煤体的一个重要结构特征，对其动态力学特性有显著影响。在实验中，选取具有明显层理结构的煤样，通过定向切割的方式制备出层理方向与加载方向不同夹角的试件，如0°（层理平行于加载方向）、45°、90°（层理垂直于加载方向）等。通过对不同层理夹角煤样的冲击实验，分析层理方向对煤动态强度、变形特性和破坏模式的影响规律。在整个实验过程中，除了上述研究的变量外，尽量保持其他实验条件一致，如煤样的制备工艺、试件尺寸、实验环境温度和湿度等，以确保实验结果能够准确反映各变量对煤动态力学特性的影响。2.3.2实验步骤试件准备：将加工好的煤样试件从保存环境中取出，用高精度游标卡尺测量其直径和高度，每个试件在不同位置测量多次并取平均值，以确保尺寸测量的准确性。用电子天平测量煤样试件的质量，记录相关数据。对于需要调整含水率的煤样，按照预定的方法进行加湿或烘干处理，并在处理后再次测量质量以监控含水率的变化，直至达到目标含水率。对于具有特定层理方向要求的试件，进行标记和分类，以便在实验中准确安装。设备调试：检查SHPB系统各部件的连接是否牢固，动力装置、子弹、入射杆、透射杆、吸收杆、阻尼器等是否处于正常工作状态。对应变片进行检查和校准，确保其粘贴牢固且测量准确。将应变片与数据采集仪连接，并对数据采集仪进行参数设置，包括采样频率、量程等，根据预计的应力波幅值和变化频率，合理设置参数以保证能够准确采集到信号。调试高速摄像机，设置合适的拍摄帧率、曝光时间和拍摄角度，确保能够清晰地记录煤样的冲击过程。安装试件：在入射杆和透射杆的接触端涂抹适量的凡士林，以减小摩擦对实验结果的影响。将煤样试件小心地放置在入射杆和透射杆之间，确保试件的轴线与入射杆和透射杆的轴线重合，对于有层理方向要求的试件，按照预定的层理方向进行安装。使用夹具或固定装置将试件夹紧，防止在冲击过程中试件发生移动或脱落。冲击实验：根据预定的实验方案，调节动力装置的气体压力，设定冲击速度。启动动力装置，发射子弹，使子弹撞击入射杆，产生应力波对煤样试件进行冲击加载。在冲击过程中，数据采集仪实时采集应变片输出的信号，记录入射波、反射波和透射波的应变数据。高速摄像机同步拍摄煤样的冲击过程，捕捉煤样从开始变形到最终破坏的整个动态过程。数据记录与整理：冲击实验结束后，停止数据采集仪和高速摄像机的工作。从数据采集仪中导出采集到的应变数据，并进行初步的处理和分析，根据三波法公式计算出煤样在冲击过程中的应力、应变率和应变历史等动态力学参数。从高速摄像机中导出拍摄的图像序列，对煤样的破坏形态和变形过程进行观察和分析，记录煤样的破坏模式、裂纹扩展方向和范围等信息。重复实验：按照上述步骤，在相同的实验条件下对同一组煤样进行多次重复实验，以提高实验结果的可靠性和精度。对于不同变量水平的煤样，如不同冲击速度、含水率或层理方向的煤样，依次进行相应的冲击实验，记录和整理实验数据。在完成所有实验后，对实验数据进行全面的分析和总结，绘制应力-应变曲线、应变率-应力曲线等图表，分析各变量对煤动态力学特性的影响规律。三、实验结果与分析3.1煤样的动态应力-应变关系3.1.1不同冲击速度下的应力-应变曲线通过分离式霍普金森压杆（SHPB）试验，获得了不同冲击速度下煤样的动态应力-应变曲线，如图1所示。从图中可以清晰地看出，随着冲击速度的增加，煤样的应力-应变曲线呈现出明显的变化特征。在低冲击速度下，如3m/s时，应力-应变曲线在初始阶段上升较为平缓，表现出一定的弹性变形特征，随着应变的增加，应力逐渐增大，当达到一定应变值后，应力增长速率加快，进入屈服阶段，随后煤样逐渐发生破坏，应力开始下降。当冲击速度提高到5m/s时，曲线的上升段斜率明显增大，表明煤样在相同应变下承受的应力更大，即动态抗压强度增加。这是因为冲击速度的提高使得煤样在短时间内受到更大的冲击力，煤样内部的颗粒来不及发生相对位移和调整，从而导致其抵抗变形的能力增强，表现为应力的快速上升。在峰值应力之后，曲线的下降段也更为陡峭，说明煤样在高冲击速度下的破坏更为剧烈，脆性特征更加明显。当冲击速度进一步提高到7m/s时，这种变化趋势更加显著。应力-应变曲线几乎呈直线上升，迅速达到峰值应力，随后急剧下降，煤样在极短的时间内发生破坏。这表明冲击速度对煤样的动态力学行为具有显著影响，高冲击速度会使煤样的强度和脆性增加，变形能力减小。[此处插入不同冲击速度下煤样的应力-应变曲线图片，图1：不同冲击速度下煤样的应力-应变曲线]3.1.2不同含水率下的应力-应变曲线不同含水率煤样的应力-应变曲线如图2所示。可以发现，随着含水率的增加，煤样的应力-应变曲线呈现出不同的变化规律。当含水率较低时，如自然含水率状态下，煤样的应力-应变曲线具有一定的斜率，在弹性阶段表现出较好的力学性能，随着应变的增加，应力逐渐上升，达到峰值后开始下降，煤样发生破坏。当含水率增加到5%时，曲线的上升段斜率有所减小，峰值应力降低，说明煤样的动态抗压强度有所下降。这是因为水分的存在削弱了煤颗粒之间的粘结力，使得煤样在受到冲击载荷时更容易发生变形和破坏。随着含水率进一步增加到10%和15%，曲线的变化趋势更为明显，峰值应力进一步降低，且在较低应变时就开始出现应力下降的趋势，表明煤样的脆性增加，变形能力进一步减弱。此外，从曲线的形态可以看出，含水率较高的煤样在破坏前的变形阶段相对较短，说明水分的增加加速了煤样的破坏过程。这是由于水分进入煤体孔隙后，在冲击载荷作用下，孔隙内的水分产生较高的孔隙水压力，进一步削弱了煤体的内部结构，导致煤样更容易发生破坏。[此处插入不同含水率下煤样的应力-应变曲线图片，图2：不同含水率下煤样的应力-应变曲线]3.1.3含层理煤样的应力-应变曲线含层理煤样在不同层理方向与加载方向夹角下的应力-应变曲线如图3所示。当层理平行于加载方向（夹角为0°）时，应力-应变曲线在初始阶段上升较为缓慢，随着应变的增加，应力逐渐增大，但增长速率相对较慢，达到峰值应力后，曲线下降较为平缓，煤样表现出一定的塑性变形特征。这是因为层理平行时，煤样在加载方向上的结构相对较为连续，能够承受一定的变形而不发生突然破坏。当层理与加载方向夹角为45°时，曲线的上升段斜率明显增大，峰值应力高于层理平行时的情况，说明此时煤样的强度有所提高。然而，在峰值应力之后，曲线下降迅速，煤样的破坏较为突然，表现出较强的脆性。这是由于45°夹角时，层理面与加载方向形成一定的剪切作用，在冲击载荷下，层理面之间更容易发生相对滑动和错动，导致煤样的破坏。当层理垂直于加载方向（夹角为90°）时，应力-应变曲线的上升段较为陡峭，峰值应力较高，但煤样在达到峰值应力后几乎瞬间破坏，曲线急剧下降，脆性特征最为明显。这是因为层理垂直时，煤样在加载方向上的结构完整性最差，层理面成为薄弱面，在冲击载荷作用下，煤样极易沿着层理面发生断裂破坏。[此处插入含层理煤样在不同层理方向与加载方向夹角下的应力-应变曲线图片，图3：含层理煤样在不同层理方向与加载方向夹角下的应力-应变曲线]3.2煤样的动态强度特性3.2.1动态抗压强度根据实验数据，计算得到了不同条件下煤样的动态抗压强度，结果如表1所示。从表中可以看出，冲击速度对煤样的动态抗压强度有显著影响。随着冲击速度的增加，煤样的动态抗压强度呈指数增长。当冲击速度从3m/s增加到7m/s时，煤样的动态抗压强度从[X1]MPa增加到[X2]MPa，增长幅度较大。这是因为在高冲击速度下，煤样内部的应变率迅速增大，煤样来不及发生塑性变形，更多地表现出弹性变形特征，从而使得其抵抗外力的能力增强，动态抗压强度提高。含水率对煤样的动态抗压强度也有明显的影响。随着含水率的增加，煤样的动态抗压强度逐渐降低。当含水率从自然含水率增加到15%时，煤样的动态抗压强度从[X3]MPa降低到[X4]MPa。这是由于水分的侵入削弱了煤颗粒之间的粘结力，降低了煤样的内部结构强度，使得煤样在冲击载荷下更容易发生破坏，从而导致动态抗压强度下降。层理方向对煤样的动态抗压强度同样具有重要影响。当层理垂直于加载方向时，煤样的动态抗压强度最低；当层理平行于加载方向时，煤样的动态抗压强度相对较高；而层理与加载方向夹角为45°时，煤样的动态抗压强度介于两者之间。这表明层理的存在改变了煤样的力学各向异性，使得煤样在不同方向上的强度表现出明显差异。[此处插入不同条件下煤样的动态抗压强度数据表格，表1：不同条件下煤样的动态抗压强度（单位：MPa）]3.2.2动态抗拉强度煤样动态抗拉强度的测试采用巴西劈裂法。在SHPB试验中，通过对煤样施加径向冲击载荷，使其在直径方向上产生拉伸应力，当拉伸应力达到煤样的抗拉强度时，煤样发生劈裂破坏。根据实验结果，得到了不同条件下煤样的动态抗拉强度，结果如图4所示。从图中可以看出，冲击速度对煤样的动态抗拉强度有一定的影响。随着冲击速度的增加，煤样的动态抗拉强度呈现出先增加后减小的趋势。在较低冲击速度范围内，冲击速度的增加使得煤样内部的应力分布更加不均匀，局部产生的拉伸应力增大，从而导致动态抗拉强度增加。但当冲击速度超过一定值后，煤样的破坏模式发生改变，更多地表现为剪切破坏，使得动态抗拉强度下降。含水率对煤样的动态抗拉强度影响较为显著。随着含水率的增加，煤样的动态抗拉强度明显降低。这是因为水分的存在削弱了煤颗粒之间的粘结力，使得煤样在拉伸载荷作用下更容易发生断裂，从而降低了动态抗拉强度。层理方向对煤样的动态抗拉强度也有较大影响。当层理平行于拉伸方向时，煤样的动态抗拉强度相对较高；当层理垂直于拉伸方向时，煤样的动态抗拉强度最低。这是由于层理垂直时，煤样在拉伸方向上的结构完整性最差，容易沿着层理面发生断裂，导致动态抗拉强度降低。[此处插入不同条件下煤样的动态抗拉强度变化曲线图片，图4：不同条件下煤样的动态抗拉强度变化曲线]3.3煤样的破坏模式与特征3.3.1宏观破坏模式不同实验条件下煤样的宏观破坏形态如图5所示。在低冲击速度下，煤样的破坏相对较为轻微，主要表现为表面出现少量裂纹，部分煤样可能会有小块煤体脱落，但整体结构基本保持完整。随着冲击速度的增加，煤样的破碎程度逐渐加剧，裂纹数量增多且扩展范围增大，煤样被破碎成多个大块，呈现出明显的脆性破坏特征。对于不同含水率的煤样，含水率较低的煤样在破坏时裂纹相对较细且数量较少，而含水率较高的煤样在破坏时裂纹更宽且数量更多，煤样破碎得更为严重，呈现出更为明显的塑性破坏特征。这是因为水分的增加使得煤样的韧性降低，在冲击载荷下更容易发生塑性变形和破碎。含层理煤样的破坏模式与层理方向密切相关。当层理平行于加载方向时，煤样主要沿着层理面发生滑移和错动，破坏后的煤样呈现出层状剥落的形态；当层理与加载方向夹角为45°时，煤样除了沿着层理面发生破坏外，还会产生与加载方向成一定角度的剪切裂纹，破坏后的煤样呈现出不规则的块状；当层理垂直于加载方向时，煤样主要沿着层理面发生断裂，破坏后的煤样呈现出明显的劈裂形态。[此处插入不同实验条件下煤样的宏观破坏形态图片，图5：不同实验条件下煤样的宏观破坏形态]3.3.2微观破坏特征利用扫描电子显微镜（SEM）对煤样的微观破坏特征进行了观察，结果如图6所示。在微观层面上，煤样的破坏主要表现为孔隙和裂隙的扩展、连通以及煤颗粒的破碎。在冲击载荷作用下，煤样内部原有的孔隙和裂隙首先发生扩展，随着载荷的增加，这些孔隙和裂隙逐渐连通，形成宏观裂纹。同时，煤颗粒之间的粘结力被破坏，煤颗粒发生破碎和分离。对于不同冲击速度下的煤样，高冲击速度下煤样的微观破坏更为严重，孔隙和裂隙的扩展更加明显，煤颗粒的破碎程度更大，且出现了更多的破碎碎屑。这表明高冲击速度使得煤样内部的能量释放更加剧烈，导致微观结构的破坏更为严重。不同含水率的煤样在微观破坏特征上也存在差异。含水率较高的煤样中，由于水分的存在，孔隙和裂隙中填充了水分，在冲击载荷作用下，水分的流动和压力变化进一步加剧了孔隙和裂隙的扩展，使得煤样的微观结构更加松散，煤颗粒之间的粘结力进一步降低。含层理煤样在微观层面上，层理面处的结构相对较为薄弱，在冲击载荷作用下，层理面处的孔隙和裂隙首先发生扩展和连通，导致层理面的破坏。当层理与加载方向夹角不同时，层理面与煤颗粒之间的受力状态不同，从而导致微观破坏特征也有所不同。例如，当层理与加载方向夹角为45°时，层理面与煤颗粒之间会产生较大的剪切应力，使得层理面处的煤颗粒更容易发生破碎和滑移，微观结构的破坏更为复杂。[此处插入煤样微观破坏特征的扫描电子显微镜图片，图6：煤样微观破坏特征的扫描电子显微镜图片]四、煤在冲击载荷下的力学模型与理论分析4.1煤的动态本构模型4.1.1现有本构模型概述岩石类材料本构模型众多，这些模型基于不同的假设和理论，旨在描述岩石在各种载荷条件下的力学行为。其中，弹性本构模型假设岩石材料在受力过程中仅发生弹性变形，应力与应变之间满足线性关系，如胡克定律所描述。这种模型在岩石受力较小、变形处于弹性阶段时具有一定的适用性，但无法解释岩石在塑性变形和破坏阶段的力学行为。弹塑性本构模型则考虑了岩石在受力过程中既有弹性变形又有塑性变形的情况。它将应变增量分为弹性应变增量和塑性应变增量两部分，弹性应变增量通过广义虎克定律计算，而塑性应变增量依据塑性增量理论确定。在经典塑性理论中，常采用屈服准则来判断材料是否进入塑性状态，如常用的莫尔-库仑屈服准则，它考虑了岩石的抗剪强度与正应力之间的关系，适用于描述岩石在复杂应力状态下的屈服和破坏。然而，该准则对于岩石的拉伸破坏等情况的描述存在一定局限性。粘弹性本构模型主要用于描述岩石材料的流变特性，即岩石的应变或应力随时间而变化的性质。此类模型通常将岩石视为由弹性元件和粘性元件组合而成的复合体，通过不同元件的组合和参数设置来模拟岩石的流变行为。例如，Maxwell模型由一个弹簧（代表弹性）和一个阻尼器（代表粘性）串联而成，它能够较好地描述岩石在恒定应力作用下的蠕变现象以及在卸载后的弹性恢复和粘性流动，但对于岩石在复杂加载条件下的力学行为模拟不够全面。对于煤在冲击载荷下的力学行为，现有本构模型存在一定的局限性。煤是一种复杂的非均质材料，内部含有大量的孔隙、裂隙以及矿物质等，其力学性质具有明显的各向异性和应变率相关性。传统的弹性和弹塑性本构模型难以准确描述煤在冲击载荷下的非线性变形、应变率强化效应以及能量耗散等特征。粘弹性本构模型虽然考虑了时间因素对岩石力学行为的影响，但对于煤在高应变率冲击载荷下的瞬态力学响应的模拟效果也不理想。此外，煤的微观结构和组成成分的差异导致其力学性质变化较大，使得通用的岩石本构模型难以全面准确地反映煤在冲击载荷下的力学行为。4.1.2基于实验结果的本构模型改进根据本次实验获得的数据，包括不同冲击速度、含水率、层理条件下煤样的应力-应变曲线、动态强度等，对现有本构模型进行修正。考虑到煤在冲击载荷下的应变率强化特性，在本构模型中引入应变率相关项。例如，将应变率作为一个变量纳入屈服准则或硬化规律中，以描述煤的动态强度随应变率增加而提高的现象。通过对实验数据的拟合和分析，确定应变率相关项的具体形式和参数，使模型能够更准确地反映煤在不同应变率下的力学行为。针对煤的非线性变形特征，对本构模型中的弹性和塑性部分进行改进。在弹性阶段，采用非线性弹性模型来描述煤的应力-应变关系，考虑煤内部孔隙、裂隙等结构对弹性变形的影响。在塑性阶段，引入损伤变量来描述煤在冲击载荷作用下内部结构的损伤演化过程，损伤变量与塑性应变相关联，随着塑性变形的增加，损伤逐渐发展，导致煤的力学性能劣化。通过实验数据确定损伤变量的演化方程和相关参数，从而使本构模型能够更好地模拟煤在冲击载荷下从弹性到塑性再到破坏的全过程。考虑煤的各向异性，尤其是层理结构对其力学性能的影响，在本构模型中引入方向相关的参数或张量。对于含层理煤样，根据层理方向与加载方向的夹角，调整模型中的弹性常数、屈服准则和硬化参数等，以反映煤在不同方向上力学性能的差异。通过对不同层理角度煤样的实验数据进行分析，确定各向异性参数的取值范围和变化规律，使改进后的本构模型能够准确描述含层理煤在冲击载荷下的力学行为。改进后的本构模型与实验结果的对比验证表明，新模型能够更准确地预测煤在冲击载荷下的应力-应变响应、动态强度以及破坏模式等。通过将模型计算结果与不同条件下的实验数据进行对比，发现改进后的模型在拟合应力-应变曲线时，相关系数更高，曲线的趋势和关键特征点（如峰值应力、峰值应变等）与实验结果更为接近。对于不同冲击速度、含水率和层理条件下煤样的动态强度预测，改进后的模型误差明显减小，能够更可靠地为工程实际提供理论支持和参考。4.2煤的断裂力学分析4.2.1裂纹起裂与扩展判据在冲击载荷下，煤中裂纹的起裂和扩展是导致其破坏的关键因素。应力强度因子判据是常用的裂纹起裂判据之一，它基于线弹性断裂力学理论。当煤体中裂纹尖端的应力强度因子K达到材料的断裂韧性K_{IC}时，裂纹开始起裂。对于I型裂纹（张开型裂纹），应力强度因子K_{I}的计算公式为：K_{I}=Y\sigma\sqrt{\pia}其中，\sigma为作用在裂纹面上的应力，a为裂纹长度，Y为与裂纹形状和加载方式有关的几何因子。当K_{I}\geqK_{IC}时，裂纹起裂。除了应力强度因子判据，还有其他一些判据用于判断裂纹的起裂和扩展。如能量释放率判据，该判据从能量的角度出发，认为当裂纹扩展单位面积时系统释放的能量（即能量释放率G）达到材料的临界能量释放率G_{C}时，裂纹起裂。能量释放率G与应力强度因子K之间存在一定的关系，对于线弹性材料，有G=\frac{K^{2}}{E}（平面应力状态，E为弹性模量）或G=\frac{(1-\nu^{2})K^{2}}{E}（平面应变状态，\nu为泊松比）。在冲击载荷下，由于加载速率高，煤体的惯性效应不可忽略，此时裂纹的起裂和扩展还受到加载速率的影响。一些学者提出了基于动态断裂力学的判据，如考虑加载速率影响的动态应力强度因子判据，动态应力强度因子K_{Id}不仅与应力和裂纹几何尺寸有关，还与加载速率相关。当K_{Id}\geqK_{IdC}（动态断裂韧性）时，裂纹起裂。4.2.2断裂过程的能量分析从能量角度分析煤在冲击载荷下的断裂过程，有助于深入理解其破坏机制。在冲击载荷作用下，煤体吸收的能量主要用于两个方面：一是使煤体发生弹性变形，储存为弹性应变能；二是用于裂纹的产生、扩展以及煤体的破碎，这部分能量以不可逆的方式耗散。在裂纹起裂前，煤体主要发生弹性变形，吸收的能量主要以弹性应变能的形式储存。弹性应变能密度U_{e}可表示为：U_{e}=\frac{1}{2}\sigma_{ij}\varepsilon_{ij}其中，\sigma_{ij}和\varepsilon_{ij}分别为应力张量和应变张量。随着冲击载荷的增加，煤体内部的应力逐渐增大，弹性应变能不断积累。当裂纹起裂后，能量开始用于裂纹的扩展和煤体的破碎。裂纹扩展过程中，需要克服裂纹表面的表面能以及煤体内部的摩擦阻力等，这部分能量消耗称为裂纹扩展能量耗散。煤体的破碎过程也伴随着能量的耗散，破碎后的煤块具有一定的动能和表面能，这些能量均来自于冲击载荷输入的能量。通过对实验数据的分析，可以计算出煤在冲击载荷下的能量耗散情况。例如，根据应力-应变曲线下的面积可以计算出煤体吸收的总能量，再通过测量破碎后煤块的动能和估算裂纹扩展的能量耗散，从而得到能量耗散与破坏机制之间的定量关系。研究发现，能量耗散与煤的破坏模式密切相关，脆性破坏时能量主要用于裂纹的快速扩展和煤体的突然破碎，能量耗散相对较大；而塑性破坏时，能量在较长时间内逐渐耗散，用于煤体的塑性变形和缓慢的裂纹扩展。此外，应变率、含水率、层理等因素也会影响能量耗散的大小和方式。高应变率下，能量耗散更快，煤体更容易发生脆性破坏；含水率增加会降低煤体的强度，使能量耗散在较低的应力水平下发生；层理结构会改变煤体的裂纹扩展路径和能量耗散方式，导致不同层理方向的煤体在冲击载荷下的能量耗散特征不同。五、工程应用与展望5.1在煤矿开采中的应用5.1.1冲击地压防治基于对煤在冲击载荷下动态力学特性的研究成果，可提出一系列预防和治理冲击地压的有效措施。在开采顺序优化方面，通过对煤体应力分布和能量积聚规律的深入了解，应优先开采应力集中程度较低、冲击危险性较小的区域。例如，在多煤层开采中，可选择先开采保护层，通过保护层开采引起的岩层移动和应力重新分布，降低被保护层的应力集中程度，释放部分弹性变形能，从而减小冲击地压发生的可能性。以某煤矿为例，该矿在开采过程中，先对上层的薄煤层进行开采作为保护层，在开采过程中，通过微震监测系统对被保护层的应力变化进行实时监测，结果显示，被保护层的应力峰值降低了约30%，有效降低了后续开采的冲击地压风险。开采方法的选择对冲击地压防治也至关重要。对于冲击地压危险区域，应避免采用高应力集中的开采方法，如采用长壁开采时，应合理控制工作面长度和推进速度。根据煤的动态力学特性，当工作面推进速度过快时，煤体在短时间内受到的扰动增加，应变率增大，导致煤体的强度和脆性增加，更容易引发冲击地压。因此，应根据具体的地质条件和煤的力学性质，确定合理的推进速度，使煤体有足够的时间释放能量，避免能量过度积聚。例如，在某冲击地压危险矿井，通过将工作面推进速度从原来的每天8m降低到每天5m，并加强对煤体应力和能量的监测，在后续的开采过程中，冲击地压的发生次数明显减少，从之前每月平均发生3次降低到每月平均1次，有效保障了开采的安全。此外，还可以采用卸压措施来降低煤体的应力集中程度。常见的卸压方法包括钻孔卸压、爆破卸压等。钻孔卸压是通过在煤体中钻进一定深度和直径的钻孔，使煤体中的应力向钻孔周围转移，从而降低钻孔附近煤体的应力集中程度。爆破卸压则是利用炸药在煤体中爆炸产生的冲击波和应力波，使煤体产生裂隙，释放能量，降低应力。在实施卸压措施时，应根据煤的动态力学参数，如动态抗压强度、断裂韧性等，合理确定钻孔的参数（如间距、深度、直径）和爆破的参数（如装药量、炮孔布置），以确保卸压效果。例如，在某煤矿的冲击地压防治实践中，通过对煤样的动态力学测试，确定了合理的钻孔间距为2m，钻孔深度为10m，采用爆破卸压时，根据煤的动态抗压强度和断裂韧性，确定了装药量为每孔1kg，炮孔呈梅花形布置。经过实际应用，该卸压方案有效降低了煤体的应力集中程度，成功预防了多起潜在的冲击地压事故。5.1.2巷道支护设计根据煤的动态力学特性设计合理的巷道支护方案，对于提高支护效果、保障巷道安全具有重要意义。在支护方式选择上，应充分考虑煤体在冲击载荷下的变形和破坏特征。对于冲击地压危险区域的巷道，锚杆支护是一种常用的支护方式，但传统的普通锚杆在抵抗冲击载荷时存在一定的局限性。因此，可采用高强度、高预应力的锚杆，如左旋无纵筋螺纹钢锚杆，其具有较高的屈服强度和抗拉强度，能够在冲击载荷作用下提供更大的锚固力，有效抑制煤体的变形和破坏。同时，配合使用锚索加强支护，锚索能够提供更大的锚固力和支护范围，与锚杆形成联合支护体系，提高巷道的整体稳定性。例如，在某煤矿的深部巷道支护中，采用了直径22mm的左旋无纵筋螺纹钢锚杆，长度为2.5m，间排距为0.8m，同时每隔2m布置一根锚索，锚索直径为17.8mm，长度为6m。通过现场监测，该支护方案在经受多次冲击地压的作用下，巷道变形量得到了有效控制，顶板下沉量控制在50mm以内，两帮移近量控制在80mm以内，保障了巷道的正常使用。支护参数的优化也是提高支护效果的关键。根据煤的动态抗压强度、抗拉强度以及巷道的埋深、围岩应力等因素，合理确定锚杆和锚索的长度、直径、间距等参数。例如，对于埋深较大、煤体强度较低的巷道，应适当增加锚杆和锚索的长度和直径，减小间距，以提高支护强度。同时，还应考虑煤体的层理结构对支护的影响，对于层理发育的煤体，应调整锚杆和锚索的布置方向，使其与层理面形成一定的夹角，增强对层理面的锚固作用。在某煤矿的含层理煤巷支护中，通过数值模拟分析，将锚杆与层理面的夹角从原来的30°调整为45°，并增加了锚索的数量，经过现场应用，巷道的稳定性得到了显著提高，层理面的滑移和错动现象得到了有效控制。此外，还可以采用一些新型的支护材料和技术，如可拉伸锚杆、让压锚索等。可拉伸锚杆在受到冲击载荷时，能够通过自身的拉伸变形吸收能量，减小对煤体的冲击作用；让压锚索则在一定的压力下能够产生让压变形，释放部分能量，避免支护结构因承受过大的压力而破坏。这些新型支护材料和技术的应用，能够更好地适应煤在冲击载荷下的力学行为，提高巷道支护的可靠性和稳定性，为煤矿安全生产提供更有力的保障。5.2研究不足与展望尽管本研究在煤在冲击载荷下的动态力学特性方面取得了一定的成果，但仍存在一些不足之处。在实验研究方面，由于实验设备和技术的限制，难以模拟煤在实际开采过程中复杂的应力状态和加载条件。实际煤矿开采中，煤体不仅受到冲击载荷的作用，还会受到地应力、瓦斯压力、地下水等多种因素的耦合作用，而本实验仅考虑了单一的冲击载荷，无法全面反映煤体的真实力学行为。此外，实验煤样的选取虽然具有一定的代表性，但无法涵盖所有煤种和地质条件，对于一些特殊煤种和复杂地质条件下煤的动态力学特性研究还不够深入。在理论分析方面，虽然对煤的动态本构模型进行了改进，但仍存在一定的局限性。改进后的本构模型虽然能够较好地描述煤在特定条件下的力学行为，但对于一些复杂的力学现象，如煤体在高应变率下的损伤演化过程、不同煤种和地质条件下的力学参数变化规律等，还不能完全准确地进行模拟和预测。同时，对于煤在冲击载荷下的断裂过程和能量耗散机制的研究还需要进一步深入，目前的理论分析还无法全面解释煤体在冲击破坏过程中的各种物理现象。未来煤动态力学特性的研究可以从多场耦合的角度展开，综合考虑地应力、瓦斯压力、地下水等因素对煤动态力学特性的影响，建立更加完善的多场耦合力学模型。例如，研究地应力与冲击载荷耦合作用下煤体的力学响应，分析地应力的大小和方向对煤体在冲击载荷下的强度、变形和破坏模式的影响规律；研究瓦斯压力对煤体动态力学特性的影响，探讨瓦斯在煤体孔隙中的吸附、解吸和渗流过程对煤体力学性能的作用机制；研究地下水对煤体动态力学特性的影响，分析水分对煤体的物理和化学作用，以及水分在煤体中的存在形式和运移规律对煤体力学行为的影响。在实验技术方面，需要进一步改进和创新，研发更加先进的实验设备和测试方法，以更准确地模拟煤在实际开采过程中的复杂加载条件和应力状态。例如，开发能够同时施加多种载荷的实验装置，实现冲击载荷与地应力、瓦斯压力、地下水压力等的联合加载；采用先进的测试技术，如数字图像相关技术（DIC）、声发射技术等，对煤样在冲击过程中的变形、损伤演化等进行实时监测和