基于数值模拟的开槽法测量混凝土工作应力研究 Study on in - situ stress of concrete measured by slotting method based on numerical simulation 1 first-author 1 corresp 徐赵东，博士，教授，从事防灾减灾与健康监测研究；E-mail：xzdsubmission@163.com 徐赵东，博士，教授，从事防灾减灾与健康监测研究；E-mail：xzdsubmission@163.com xzdsubmission@163.com 2 3 东南大学 中国-巴基斯坦重大基础设施智慧防灾一带一路联合实验室，江苏 南京，210096 东南大学 中国-巴基斯坦重大基础设施智慧防灾一带一路联合实验室，江苏 南京，210096 China-Pakistan Belt and Road Joint Laboratory on Smart Disaster Prevention of Major Infrastructures, Southeast University, Nanjing 210096, China China-Pakistan Belt and Road Joint Laboratory on Smart Disaster Prevention of Major Infrastructures, Southeast University, Nanjing 210096, China 江苏东南特种技术工程有限公司，江苏 南京，210008 江苏东南特种技术工程有限公司，江苏 南京，210008 Southeast Special Engineering Co. Ltd., Nanjing 210008, China Southeast Special Engineering Co. Ltd., Nanjing 210008, China 南京东瑞减震控制工程有限公司，江苏 南京，210033 南京东瑞减震控制工程有限公司，江苏 南京，210033 Nanjing Dongrui Damping Control Technology Co. Ltd., Nanjing 210033, China Nanjing Dongrui Damping Control Technology Co. Ltd., Nanjing 210033, China 杨幼平 提出一种基于应力释放原理测量混凝土工作应力的分步开槽检测方法，针对单向应力状态下的混凝土受压构件，采用数值模拟的方法分析开槽形状、开槽间距、开槽长度以及开槽深度等因素对应力释放程度的影响。研究结果表明：开槽形状为横槽时，开槽间距越小，开槽长度越长，应力释放越快；开槽形状为方槽时，开槽间距越小，应力释放越快，而开槽长度对方槽的应力释放影响较小；当开槽深度为0~ 40 mm时，方槽与相同开槽间距、长宽比 λ > 1的横槽可以相互替换，两种情况中心测点处的应力差值较小。当应力释放时，槽内区域比槽外区域的应力变化更敏感和准确，开槽法解决了圆孔法孔内应变信息无法连续输出的问题。通过测量每步开槽深度下的应力释放值，与数值模拟的结果相结合，利用多步优化结果反推测点处工作应力，减小测量误差。 A step-by-step slotting measurement method based on the stress-release method was proposed to measure in - situ stress of concrete. For concrete compressed members under unidirectional stress, numerical simulation was used to analyze the effects of slotting shape, slotting spacing, slotting length and slotting depth on stress-release degree. The results show that the stress is released faster with the smaller slotting spacing and the longer slotting length when the slotting shape is horizontal. The stress is released faster with the smaller slotting spacing when the slotting shape is square and the slotting length has less influence on the stress release of the square slot. When the slotting depth is 0-40 mm, the square slot can be replaced with the horizontal slot of the same slotting spacing and aspect ratio λ > 1, and the stress difference at the center measurement point of two cases is small. When the stress is released, the stress change in the area inside the slot is more sensitive and accurate than that in the area outside the slot. The slotting method solves the problem that the strain information inside the hole of the circular hole method cannot be continuously output. By measuring the stress release value at each slotting depth step and combining it with the results of numerical simulation, the multi-step optimization results can be used to infer the working stress at the measurement point and reduce the measurement error. 建筑结构 应力释放法 数值模拟 工作应力 开槽法 building structure stress release method numerical simulation in - situ stress slotting method http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44706592&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44706595&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44706594&type= 14.66796875 14.66793823 在混凝土结构的寿命期内，应力衰减、混凝土徐变、施工偏差和不均匀沉降会导致混凝土结构应力发生变化，工作应力的理论计算总是与实际应力相差很大，因此准确评估混凝土结构内的工作应力对于确保结构安全性和可靠性非常重要。测量现有混凝土结构的工作应力一直是一个难题，混凝土材料的不均匀性和随机性、应变测量技术等因素导致应力测量结果存在较大的离散性和误差。 局部破损检测技术即应力释放法，最初由MATHAR [ 1 1 ] 提出，用来对钢结构构件中的残余应力进行测量。对金属构件所测区域进行切割或钻孔，使所测区域的应力得到释放，测量应力释放前后测点处的应变变化量，经计算得到残余应力 [ 2 2 - 3 3 2-3 ] 。 近几十年来，科研人员对应力释放法进行了大量的研究，并逐渐将其应用于混凝土材料的应力检测中。为了克服解析方法仅适用于应力沿厚度方向均匀分布的薄板的局限性，SCHAJER [ 4 4 - 6 6 4-6 ] 提出了以标定系数为核心的积分法与幂级数法来识别沿深度方向的非均匀分布应力场，其适用性、稳定性和精度都得到了显著增强。 TRAUTNER等 [ 7 7 - 8 8 7-8 ] 将应力释放法应用于混凝土结构，使用基于表面位移的影响函数，结果具有了良好的理论精度，但由于释放位移量级极小，实际应用难度较大。影响函数是标定系数的连续形式，目前已发展为广泛使用的方法 [ 9 9 - 12 12 9-12 ] 。国内外学者逐步将影响函数应用到混凝土工作应力的检测之中，对影响函数不断进行优化以提高识别精度 [ 13 13 - 14 14 13-14 ] 。 混凝土应力释放法按钻孔或切割形式的不同，可分为盲孔法、圆孔法和开槽法 [ 15 15 - 19 19 15-19 ] 。开槽法是在圆孔法的基础上进行改进和创新得来的，按开槽形状可以分为直线形开槽法(横槽)和方形开槽法(方槽)。开槽法在实际工程中易于施工，尤其是在只关注单向应力问题时，对混凝土结构产生的损伤较小。因此，近年来部分学者对开槽法进行了深入研究，并将其应用到混凝土的应力检测中 [ 20 20 - 25 25 20-25 ] 。但是实际混凝土结构大多处于复杂受力状态，开槽法的测量结果必定与混凝土结构真实受力状态存在一定的差别，这一缺点限制了开槽法在实际工程中的应用。 本文作者对应力释放法中的开槽法进行研究，针对单向应力状态下的混凝土受压构件，采用数值模拟的分析方法研究开槽形状、开槽间距、开槽长度等因素在不同开槽深度下对应力释放程度的影响，并提出了一种新的分步测量方法。 1 基本原理 1 开槽法是一种局部损伤检测方法，用于评估混凝土结构内的工作应力。开槽过程中，开槽区域内外应力逐渐释放，开槽区域周围的混凝土发生弹性变形。应变传感器根据工作应力方向固定在测点上，随着开槽深度增加，应变传感器开始测量混凝土释放的弹性应变。开槽达到一定深度后，测点处混凝土的工作应力完全释放。根据应变传感器测得的应变变化值，通过混凝土弹性模量计算构件的工作应力。 <math display="block" id="M1"><mi>σ</mi><mo>=</mo><mi>E</mi><mi>ε</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44706442&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44706441&type= 8.89000034 2.62466669 式中： σ 为混凝土的工作应力； E 为混凝土的弹性模量； ε 为释放的弹性应变。 2 数值模拟与分析 1 2.1 　计算模型与参数 2 本文以单向应力状态下的受压构件为研究对象进行数值模拟，使用有限元软件建模分析，选用的模型长×宽×高为3 000 mm×400 mm×400 mm，软件内模型示意图如 图1 图1 所示。模型中的混凝土强度等级为C30，钢筋的强度等级为HRB400。为方便后续对计算结果进行分析，对开槽间距 H 、开槽长度 L 、中心测点以及应力释放区域长宽比 λ 等因素进行规定，如 图2 图2 所示。其中应力释放区域长宽比 λ 的计算公式为 <math display="block" id="M2"><mi>λ</mi><mo>=</mo><mfrac><mrow><mi>L</mi></mrow><mrow><mi>H</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44706465&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44706458&type= 9.22866726 7.28133297 式中： λ 为应力释放区域长宽比； L 为开槽长度即应力释放区域长度； H 为开槽间距即应力释放区域宽度。 (a) 横槽；(b) 方槽 测点应力释放率 α 的计算公式为 <math display="block" id="M3"><mi>α</mi><mo>=</mo><mn mathvariant="normal">1</mn><mo>-</mo><mfrac><mrow><msub><mrow><mi>σ</mi></mrow><mrow><mi>y</mi></mrow></msub></mrow><mrow><msub><mrow><mi>σ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44706479&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44706478&type= 14.98600006 8.63599968 式中： σ y 为不同开槽深度下的测点处的 Y 向应力； σ 0 为测点处的初始 Y 向应力。 2.2 　压应力对应力释放程度的影响 2 针对压应力(即构件工作应力)对应力释放的影响程度进行分析，采用轴压比 u 进行控制，选取轴压比 u 为0.1、0.3、0.5和0.7的4种工况。开槽形状选取方槽 H =100 mm、 L =100 mm、 λ =1，开槽深度取10 mm为一步，从0 mm分步计算至150 mm得到不同开槽深度下开槽区域中心测点处 Y 向应力并计算得到每一步开槽深度下的应力释放率。 Y 向应力的分析结果如 图3 图3 所示，应力释放率的分析结果如 图4 图4 所示。 由 图3 图3 可以看出：压应力只对构件测点处初始 Y 向应力和每一步应力释放值有影响，对应力完全释放即应力释放率为0时的开槽深度没有影响。随着开槽深度增加， Y 向应力从初值逐渐减小，4种工况的 Y 向应力均在开槽深度为36 mm时完全释放。应力释放率到达0时，开槽深度继续增加， Y 向应力随之反向增大，在开槽深度为50 mm时达到反向最大值，其中初始 Y 向应力越大的工况反向应力最大值越大。当开槽深度为50~150 mm时， Y 向应力逐渐减小至0 MPa附近。 由 图4 图4 可以看出：在开槽深度相同时，4种轴压比下测点的应力释放率基本相同。由此可见，压应力对每一步的应力释放率没有影响，即不同开槽深度下的应力释放程度与构件当前的工作应力无关。因此，开槽法可以广泛用于检测混凝土构件的工作应力，不同开槽深度下的应力释放值不受构件工作应力的影响。 2.3 　开槽长度对应力释放程度的影响 2 2.3.1　横槽的开槽长度对应力释放程度的影响 3 针对横槽的开槽长度 L 对应力释放程度的影响进行分析，开槽间距 H 选取50 mm，开槽长度 L 选取30、35、40、45、50、60、70、80、90和100 mm这10种工况进行分析，即长宽比 λ 分别为0.6、0.7、0.8、0.9、1.0、1.2、1.4、1.6、1.8和2.0。从0 mm分步计算至100 mm得到不同开槽深度下开槽区域中心测点处 Y 向应力，结果如 图5 图5 所示。 由 图5 图5 可以看出：开槽长度对横槽应力释放程度的影响可以根据长宽比 λ 分为两种情况进行分析，即长宽比 λ ≤1与 λ > 1。 当长宽比 λ ≤1时，即开槽长度 L 小于等于开槽间距 H 时， Y 向应力随开槽深度的变化趋势一致，但每种工况之间差值较大。随着开槽深度增加， Y 向应力从初值开始减小，在开槽深度为30 mm时应力释放率达到最大值，随后 Y 向应力逐渐稳定到某一负值附近。开槽长度 L 越长，应力释放越快，应力完全释放时的开槽深度也较小。但当 λ =0.6时，不管开槽深度为多大，中心测点的 Y 向应力均未达到完全释放的程度。 当长宽比 λ > 1时，即开槽长度 L 大于开槽间距 H 时， Y 向应力随开槽深度的变化趋势一致，但每种工况之间差值较小。 λ > 1的5种工况中 Y 向应力最大差值出现在开槽深度为30 mm时， Y 向应力最大差值与初始 Y 向应力的比值为7.24%；每步开槽深度下 Y 向应力最大差值与初始 Y 向应力的平均比值为2.74%。由此可见，长宽比 λ > 1的情况下，不同开槽长度之间的 Y 向应力差值较小。因此，在检测切割时只要保证长宽比 λ > 1，检测结果较为精准。 2.3.2　方槽的开槽长度对应力释放程度的影响 3 针对方槽的开槽长度 L 进行分析，与2.3.1中研究横槽开槽长度的工况相同，开槽间距 H 选取50 mm，开槽长度 L 选取30、35、40、45、50、60、70、80、90和100 mm这10种工况进行分析，即长宽比 λ 分别为0.6、0.7、0.8、0.9、1.0、1.2、1.4、1.6、1.8和2.0。从0 mm分步计算至100 mm得到不同开槽深度下开槽区域中心测点处 Y 向应力，结果如 图6 图6 所示。 由 图6 图6 可以看出： Y 向应力随开槽深度的变化趋势一致， Y 向应力从初值开始减小，在开槽深度为30 mm时应力释放率达到最大值，随后 Y 向应力逐渐稳定到0 MPa附近。 Y 向应力最大差值出现在开槽深度为30 mm时， Y 向应力最大差值与 Y 向应力初值的比值为7.26%；每步开槽深度下 Y 向应力最大差值与 Y 向应力初值的平均比值为2.81%。由此可见，当开槽形状为方槽时，开槽长度对不同开槽深度下中心测点的 Y 向应力影响较小。因此，在检测切割时，方槽可以用于结构狭小区域，开槽间距一定时，开槽长度变化所导致的应力释放值差值较小，检测结果较为精准。 2.3.3　横槽与方槽的对比 3 由2.3.2可知，在开槽形状为方槽时，开槽长度 L 对不同开槽深度下中心测点的 Y 向应力影响较小，因此选取方槽长宽比 λ =1的工况为代表值，代入到2.3.1节中横槽不同开槽长度下 Y 向应力随开槽深度 H 的变化中。为便于区分，方槽长宽比用 λ 0 表示，结果如 图7 图7 所示。 由 图7 图7 可见：当0≤ H ≤40 mm时，长宽比 λ 0 =1的方槽 Y 向应力与长宽比 λ > 1的横槽 Y 向应力之间变化趋势一致，差值较小。当开槽深度为30 mm时， λ 0 =1的方槽 Y 向应力与 λ =1.2的横槽 Y 向应力之间的差值与初始 Y 向应力的比值为3.57%；与 λ =2.0的横槽 Y 向应力之间的差值与初始 Y 向应力的比值为3.67%。当40 < H ≤80 mm时，长宽比 λ 0 =1的方槽 Y 向应力逐渐减小至0附近，与长宽比 λ > 1的横槽 Y 向应力之间的差值逐渐变大。当80 < H ≤ 100 mm时，长宽比 λ 0 =1的方槽 Y 向应力逐渐稳定在0附近，与长宽比 λ =1的横槽 Y 向应力之间差值较小。因此，在开槽深度较小时，方槽与相同开槽间距、长宽比 λ > 1的横槽之间可以相互替代，两种开槽形状的中心测点的 Y 向应力之间差值较小，趋势一致。 2.4 　开槽间距对应力释放程度的影响 2 针对开槽间距 H 进行分析，开槽长度 L 选取50、60、70、80、90和100 mm这6种工况进行分析。为了避免开槽长度 L 不同对结果造成影响，控制6种工况的长宽比 λ =1即开槽长度 L 均等于开槽间距 H 。从0 mm分步计算至100 mm得到不同开槽深度下开槽区域中心测点处 Y 向应力，结果如 图8 图8 和 图9 图9 所示。 由 图8 图8 和 图9 图9 可见：横槽与方槽不同开槽间距下的应力释放趋势相同。在开槽过程中，中心测点的 Y 向应力将经历三个阶段。在第一阶段，随着开槽深度增加， Y 向应力逐渐减小，直到应力完全释放。对于应力完全释放的开槽深度称为零应力深度或应力完全释放深度，该深度在第一阶段结束后首次达到。开槽深度继续增加，应力释放进入第二阶段。此时，应力沿相反方向增加，直到达到最大值(远小于初始 Y 向应力)。最后进入第三阶段， Y 向应力从反向最大值再次逐渐接近于 0 MPa。在实际应力测量中，第一阶段的应力释放程度的变化应该是关注的焦点。因此，开槽间距越小，中心测点处的应力释放速率越大，应力完全释放深度越小。 2.5 　数值模拟与试验结果对比分析 2 为验证数值模拟结果的准确性和适用性，选取文献[ 22 22 ]中6个测点的试验结果进行对比分析，如 表1 表1 所示。由 表1 表1 可以看出：在相同的开槽深度下，横槽与方槽之间的应力释放率差值较小，这也验证了前文得到的结论：长宽比 λ= 1的方槽与相同开槽间距 H 、长宽比 λ > 1的横槽之间可以相互替代，两种情况中心测点的应力差值较小。试验应力释放率与数值模拟结果之间的差值均较小，验证了数值模拟方法的准确性和开槽法的适用性。值得注意的是，文献[ 22 22 ]中扣除了水、温度和切割扰动影响后的试验结果具有较高的精度，但实测过程中对测点应变的影响与许多因素有关：包括混凝土含水量、骨料质量、钢筋位置和应变片尺寸等，由此产生扰动对测点应变的影响也会有所不同，这些因素的影响也将是未来研究的重点。 注： 开槽深度均为34 mm，测点1~3为 H =100 mm、 L > 100 mm、 λ > 1的横槽，测点4~6为 H =100 mm、 L= 100 mm、 λ= 1的方槽。 测点 开槽形状 初始应变/10 -6 开槽后应变/10 -6 试验应力释放率/% 数值模拟应力释放率/% 应力释放率差值/% 1 横槽 -155.0 -3 98.1 97.6 0.5 2 横槽 -193.0 -6 96.9 97.6 -0.7 3 横槽 -149.0 0 100.0 97.6 2.4 4 方槽 -147.4 -9 93.9 94.5 -0.6 5 方槽 -144.2 -8 94.4 94.5 -0.1 6 方槽 -178.1 7 103.9 94.5 9.4 3 开槽区域的应力分布规律分析 1 根据应力释放的原理对混凝土进行开槽处理后，开槽区域周围的混凝土将发生应力释放，导致局部应力重新分布。对于测量混凝土工作应力的槽孔，相邻开槽的位置应不受其他开槽区域的局部影响。在实际测量中，应在开槽区域周围合理选取测点，保证能够灵敏、准确地测量开槽引起的应变变化。测点选取不合理、测点处应力释放率小以及应变传感器的灵敏度弱等因素，都可能导致应变测量失败。 提取槽内区域和槽外区域的各测点不同开槽深度时的 Y 向应力，分析开槽区域的应力释放范围。选取 H =50 mm、 L =50 mm、 λ =1的横槽与方槽进行分析，由于横槽与方槽均为对称开槽，因此测点从中心测点开始向一侧选取，测点位置如 图10 图10 所示。开槽深度从0 mm分步计算至100 mm得到不同开槽深度下各测点处 Y 向应力，结果如 图11 图11 和 图12 图12 所示。 数据单位：mm。 横槽与方槽的槽内区域内选取5个测点，5个测点与槽内边缘的距离分别为5、10、15、20和25 mm(中心测点)。由 图11 图11 和 图12 图12 可见：25 mm(中心测点)与20 mm两个测点之间 Y 向应力差值较小，横槽每步开槽深度下的平均差值与初始 Y 向应力的比值为1.43%；方槽的比值为1.45%。开槽深度为0~30 mm时，10、15、20和25 mm四条应力释放曲线变化趋势一致，越靠近凹槽，应力释放速度越快；5 mm工况下的应力释放曲线与中心测点的差值较大，这是由于距离凹槽过近，受到应力集中等因素的影响。考虑到凹槽切割扰动、测量差值较大等因素的影响，离凹槽过近的位置不适合进行应变测量，因此建议选取中心测点以及上下10 mm范围内作为外贴应变片的位置，测量结果稳定，误差在可接受的范围内。 横槽与方槽的槽外区域选取7个测点，7个测点与槽外边缘的距离分别为20、40、60、80、100、140和180 mm。测点越靠近凹槽，应力变化越大，应力释放曲线变化趋势越接近于槽内区域测点的应力释放曲线变化趋势。因此，圆孔法测点通常选取距离切割圆孔外0~15 mm的位置。测点距离凹槽180 mm时，应力释放曲线接近为直线，开槽深度为100 mm时，横槽的应力释放率为3.22%；方槽的应力释放率为5.89%。因此，相邻凹槽至少应保持在200 mm之外，可认为不受相邻凹槽的影响。 对于横槽开槽法和方槽开槽法，与圆孔法相比，其优点在于可以将测点选取在槽内区域，解决了圆孔法孔内应变信息不能连续输出的问题。圆孔法由于钻孔机器的影响，测点只能选取在圆孔外部区域，因此测量结果受测点与圆孔之间的距离影响较大，越靠近圆孔测量结果越精准，但同时会受切割扰动影响。而开槽法的测点可以选取在测量区域中心位置，测量结果精准，受切割扰动等因素影响小。 4 分步测量方法 1 在使用圆孔法或开槽法进行实际测量时，通常的做法是钻孔或开槽至零应力深度或应力完全释放深度，得到测点处的释放应变，从而反推测点处的工作应力。但是这种方法测量得到的工作应力误差较大，分析其原因：第一，零应力深度与开槽形状、开槽间距、开槽长度等因素关系密切，在切割过程中开槽长度等因素容易发生变化，因此，零应力深度也会随之变化，导致误差变大；第二，在通常情况下，由于切割扰动等因素的影响，难以精确控制混凝土切割机器切割深度恰好达到零应力深度，导致误差存在。 为了减小测量误差，提出了一种新的分步测量方法。设置合理的开槽分步距离，通过测量每一步开槽深度下的应力释放值，与数值模拟的结果相结合，利用多步优化结果反推测点处工作应力。分步测量方法与常规测量方法最大的区别在于不需要时刻关注零应力深度的变化，而是可以采用合理的分析步距和分析步数，最大程度地减小误差。 由2.2分析结果可知，测点处现存工作应力对每一步的应力释放率没有影响。设置每一步开槽深度为5 mm，由数值模拟可以得到开槽深度从 0 mm开始每步深度(5 mm)下的应力释放率 α 0 ， α 1 ， α 2 ， α 3 ，…， α n ( α 0 =0)。同时，在构件上采用相同开槽形状从0 mm开始切割，得到每步深度(5 mm)下中心测点处的应变测量值 ε 0 ， ε 1 ， ε 2 ， ε 3 ，…， ε n ( ε 0 =0)。因此，工作应力 σ 可以通过以下公式进行计算： <math display="block" id="M4"><mi>σ</mi><mo>=</mo><mfrac><mrow><mi>E</mi></mrow><mrow><mi>n</mi></mrow></mfrac><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>n</mi></mrow></munderover></mstyle><mrow><mfrac><mrow><msub><mrow><mi>ε</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>-</mo><msub><mrow><mi>ε</mi></mrow><mrow><mi>i</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msub></mrow><mrow><msub><mrow><mi>α</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>-</mo><msub><mrow><mi>α</mi></mrow><mrow><mi>i</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfrac></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44706613&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44706590&type= 26.24666595 8.72066593 式中： n 为测量步数； α i 为数值模拟中第 i 步的应力释放率； ε i 为实际测量中第 i 步的应变测量值。 5 结论 1 1) 提出了一种分步方法测量混凝土的工作应力，无需关注开槽过程中的零应力深度变化，通过测量多步开槽深度下的应力释放值，与数值模拟的结果相结合，减小测量误差。 2) 压应力(工作应力)只对初始应力和每一步应力释放值有影响，对不同开槽深度下的应力释放率没有影响。 3) 开槽长度对长宽比 λ ≤1的横槽影响较大，开槽长度 L 越大，应力释放速率越大，应力完全释放深度越小；对长宽比 λ > 1的横槽影响较小。开槽长度对方槽的影响较小，在开槽深度较小时，方槽与相同开槽间距、长宽比 λ > 1的横槽之间可以相互替代，两种情况中心测点的 Y 向应力之间差值较小。 4) 横槽与方槽在不同开槽间距下的应力释放趋势相同，开槽间距越小，应力释放速率越大，应力完全释放深度越小。 5) 开槽法解决了圆孔法孔内应变信息不能连续输出的问题。应力释放时，槽内区域比槽外区域的应力变化更敏感和准确，选取中心测点以及上下10 mm范围内粘贴应变片，测量结果稳定。 参考文献 MATHAR J . Determination of initial stresses by measuring the deformations around drilled holes [J]. Journal of Fluids Engineering , 1934 , 56 ( 3 ): 249 - 254 . Determination of initial stresses by measuring the deformations around drilled holes RENDLER N J , VIGNESS I . Hole-drilling strain-gage method of measuring residual stresses [J]. Experimental Mechanics , 1966 , 6 ( 12 ): 577 - 586 . Hole-drilling strain-gage method of measuring residual stresses ASTM E837-08e 1. Standard test method for determining residual stresses by the hole-drilling strain-gage method [S]. SCHAJER G S . Application of finite element calculations to residual stress measurements [J]. Journal of Engineering Materials and Technology , 1981 , 103 ( 2 ): 157 - 163 . Application of finite element calculations to residual stress measurements SCHAJER G S . Measurement of non-uniform residual stresses using the hole-drilling method. part I: stress calculation procedures [J]. Journal of Engineering Materials and Technology , 1988 , 110 ( 4 ): 338 - 343 . Measurement of non-uniform residual stresses using the hole-drilling method. part I: stress calculation procedures SCHAJER G S . Measurement of non-uniform residual stresses using the hole-drilling method. part II: practical application of the integral method [J]. Journal of Engineering Materials and Technology , 1988 , 110 ( 4 ): 344 - 349 . Measurement of non-uniform residual stresses using the hole-drilling method. part II: practical application of the integral method TRAUTNER C , MCGINNIS M , PESSIKI S . Analytical and numerical development of the incremental core-drilling method of non-destructive determination of in situ stresses in concrete structures [J]. The Journal of Strain Analysis for Engineering Design , 2010 , 45 ( 8 ): 647 - 658 . Analytical and numerical development of the incremental core-drilling method of non-destructive determination of in situ stresses in concrete structures TRAUTNER C , MCGINNIS M , PESSIKI S . Application of the incremental core-drilling method to determine in-situ stresses in concrete [J]. ACI Materials Journal , 2011 , 108 ( 3 ): 290 - 299 . Application of the incremental core-drilling method to determine in-situ stresses in concrete SCHAJER G S . Hole-drilling residual stress profiling with automated smoothing [J]. Journal of Engineering Materials and Technology , 2007 , 129 ( 3 ): 440 - 445 . Hole-drilling residual stress profiling with automated smoothing BEGHINI M , BERTINI L , MORI L F , et al . Genetic algorithm optimization of the hole-drilling method for non-uniform residual stress fields [J]. The Journal of Strain Analysis for Engineering Design , 2009 , 44 ( 1 ): 105 - 115 . Genetic algorithm optimization of the hole-drilling method for non-uniform residual stress fields BEGHINI M , BERTINI L , MORI L F . Evaluating non-uniform residual stress by the hole-drilling method with concentric and eccentric holes. part I: definition and validation of the influence functions [J]. Strain , 2010 , 46 ( 4 ): 324 - 336 . Evaluating non-uniform residual stress by the hole-drilling method with concentric and eccentric holes. part I: definition and validation of the influence functions BEGHINI M , BERTINI L , MORI L F . Evaluating non-uniform residual stress by the hole-drilling method with concentric and eccentric holes. part II: application of the influence functions to the inverse problem [J]. Strain , 2010 , 46 ( 4 ): 337 - 346 . Evaluating non-uniform residual stress by the hole-drilling method with concentric and eccentric holes. part II: application of the influence functions to the inverse problem 张煜 , 阮欣 , 石雪飞 , 等 . 基于应力释放原理的桥梁混凝土结构单轴原位应力识别方法 [J]. 工程力学 , 2014 , 31 ( 11 ): 154 - 160 . 基于应力释放原理的桥梁混凝土结构单轴原位应力识别方法 ZHANG Yu , RUAN Xin , SHI Xuefei , et al . Identification method for uniaxial in situ stresses in bridge concrete structures [J]. Engineering Mechanics , 2014 , 31 ( 11 ): 154 - 160 . Identification method for uniaxial in situ stresses in bridge concrete structures RUAN Xin , ZHANG Yu . In-situ stress identification of bridge concrete components using core-drilling method [J]. Structure and Infrastructure Engineering , 2015 , 11 ( 2 ): 210 - 222 . In-situ stress identification of bridge concrete components using core-drilling method PARIVALLAL S , RAVISANKAR K , NAGAMANI K , et al . In situ stress assessment in concrete structures under bi-axial condition by trepanning technique [J]. Experimental Techniques , 2015 , 39 ( 4 ): 57 - 63 . In situ stress assessment in concrete structures under bi-axial condition by trepanning technique DENG Nianchun , TANG Pengfei . Research on in situ stress measurements in reinforced concrete beams based on the core-drilling method [J]. Advances in Civil Engineering , 2020 , 2020 : 1 - 11 . Research on in situ stress measurements in reinforced concrete beams based on the core-drilling method 刘永淼 . 环孔法测试混凝土工作应力试验研究 [D]. 杭州 : 浙江大学 , 2006 : 29 - 37 . 环孔法测试混凝土工作应力试验研究 LIU Yongmiao . The experimental research on the working stress measurement of concrete by the ring-hole method [D]. Hangzhou : Zhejiang University , 2006 : 29 - 37 . The experimental research on the working stress measurement of concrete by the ring-hole method 邱兆国 , 张凤鹏 . 混凝土梁内现存应力钻孔释放过程数值模拟 [J]. 哈尔滨工业大学学报 , 2011 , 43 ( S1 ): 232 - 236 . 混凝土梁内现存应力钻孔释放过程数值模拟 QIU Zhaoguo , ZHANG Fengpeng . Numerical simulation of the releasing process of existing stress in concrete beam by drilling hole method [J]. Journal of Harbin Institute of Technology , 2011 , 43 ( S1 ): 232 - 236 . Numerical simulation of the releasing process of existing stress in concrete beam by drilling hole method MCGINNIS M J , PESSIKI S . Experimental study of the core-drilling method for evaluating in situ stresses in concrete structures [J]. Journal of Materials in Civil Engineering , 2016 , 28 ( 2 ): 04015099 . Experimental study of the core-drilling method for evaluating in situ stresses in concrete structures DABLI A G , BAMBOLE A N , BAJORIA K M . Evaluation of in-place stress in concrete by incremental hole drilling [J]. ACI Materials Journal , 2020 , 117 ( 4 ): 27 - 35 . Evaluation of in-place stress in concrete by incremental hole drilling RICKERT T . Residual stress measurement by ESPI hole-drilling [J]. Procedia CIRP , 2016 , 45 : 203 - 206 . Residual stress measurement by ESPI hole-drilling 沈旭凯 . 开槽法测试混凝土工作应力试验与研究 [D]. 杭州 : 浙江大学 , 2007 : 33 - 58 . 开槽法测试混凝土工作应力试验与研究 SHEN Xukai . The testing and research on the working stress measurement of concrete by the grooving method [D]. Hangzhou : Zhejiang University , 2007 : 33 - 58 . The testing and research on the working stress measurement of concrete by the grooving method 王柏生 , 沈旭凯 , 林湘祁 . 开槽法测试混凝土工作应力的试验研究 [J]. 浙江大学学报(工学版) , 2010 , 44 ( 9 ): 1754 - 1759 . 开槽法测试混凝土工作应力的试验研究 WANG Baisheng , SHEN Xukai , LIN Xiangqi . Experimental study on working stress measurement of concrete by grooving method [J]. Journal of Zhejiang University(Engineering Science) , 2010 , 44 ( 9 ): 1754 - 1759 . Experimental study on working stress measurement of concrete by grooving method 李新凯 , 马松林 , 侯相深 . 钻环法测量水泥路面表面残余应力的初步研究 [J]. 工程力学 , 2010 , 27 ( 2 ): 161 - 165 . 钻环法测量水泥路面表面残余应力的初步研究 LI Xinkai , MA Songlin , HOU Xiangshen . A preliminary research on using core ring method to measure residual stress in concrete pavement surface [J]. Engineering Mechanics , 2010 , 27 ( 2 ): 161 - 165 . A preliminary research on using core ring method to measure residual stress in concrete pavement surface 李新凯 , 侯相深 , 马松林 . 钻环测量水泥混凝土板表面残余应力的试验研究 [J]. 工程力学 , 2012 , 29 ( 7 ): 244 - 249 . 钻环测量水泥混凝土板表面残余应力的试验研究 LI Xinkai , HOU Xiangshen , MA Songlin . The research of core ring method to measure residual stress on the surface of cement concrete slab [J]. Engineering Mechanics , 2012 , 29 ( 7 ): 244 - 249 . The research of core ring method to measure residual stress on the surface of cement concrete slab 10.11817/j.issn.1672-7207.2023.05.022 TU375 A 1672-7207（2023）05-1885-09 张一, 徐赵东, 李今保, 等. 基于数值模拟的开槽法测量混凝土工作应力研究[J]. 中南大学学报(自然科学版), 2023, 54(5): 1885-1893. 张一, 徐赵东, 李今保, 等. 基于数值模拟的开槽法测量混凝土工作应力研究[J]. 中南大学学报(自然科学版), 2023, 54(5): 1885-1893. ZHANG Yi, XU Zhaodong, LI Jinbao, et al. Study on in - situ stress of concrete measured by slotting method based on numerical simulation[J]. Journal of Central South University(Science and Technology), 2023, 54(5): 1885-1893. ZHANG Yi, XU Zhaodong, LI Jinbao, et al. Study on in-situ stress of concrete measured by slotting method based on numerical simulation[J]. Journal of Central South University(Science and Technology), 2023, 54(5): 1885-1893. 国家自然科学基金重点资助项目(52130807)；江苏省前沿引领技术基础研究重大项目(BK20222006) Project(52130807) supported by the National Natural Science Foundation of China; Project(BK20222006) supported by the Fundamental Research on Frontier Leading Technology of Jiangsu Province 2022-08-06 2022-10-12 2023-05-26 2023-06-20 38910722 中南大学学报(自然科学版) 5 54