Method, device and equipment for characterizing dynamic load response characteristics of an anchor rod and storage medium

Abstract

The application provides a kind of anchorage rod dynamic load response characteristic characterization method, device, equipment and storage medium, by obtaining preprocessed test data set;Preprocessed test data set includes the load of the rod body of anchorage rod, rod body cumulative elongation;By moving average algorithm, standard deviation algorithm, confidence interval algorithm, maximum minimum value algorithm or quantile algorithm, target data set is obtained;With rod body cumulative elongation as abscissa, with rod body load as ordinate, draw to obtain the target curve corresponding to target data set;Target curve includes upper envelope line, lower envelope line and median line.The application can obtain the target curve of anchorage rod dynamic load response characteristic, the trend and volatility of test data can be intuitively reflected through the target curve, which helps to detect abnormal values, and can provide key basic data support in the selection of anchorage rod support, optimization of roadway support scheme, prevention and control of roadway rock burst disaster and other aspects of mine roadway.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of anchor rod support engineering, and in particular to a method and device for characterizing dynamic load response characteristics of an anchor rod, an equipment and a storage medium. BACKGROUND

[0002] An anchor rod is an important component in the field of engineering construction, which can withstand tensile force generated by soil pressure, water pressure or wind load, and is used to maintain the stability of structures or reinforce the rock-soil environment. For example, in coal mining, a large number of roadways are often excavated under the mine, and anchor rods are installed in the rock layers around the roadways to form a whole stable rock belt, which uses the anchor rods and surrounding rock to jointly maintain the stability of the roadway. In addition, anchor rods are also applied in the fields of civil tunnel construction, slope reinforcement, dam reinforcement, etc.

[0003] With the increasing depth of engineering construction year by year, the construction environment is also changing. Due to the high ground stress in deep mines, a large amount of strain energy is stored in the rock mass. When the energy of this part exceeds the supporting capacity of the anchor rod support system, rock burst is likely to occur. Therefore, the selection of anchor rods needs to consider not only the safety and quality of construction, but also the cost of engineering construction.

[0004] It can be seen that, in the face of the problem that the current construction environment changes and the adaptability of anchor rods cannot match, it is an important task for the industry to study the response characteristics of anchor rods under impact load in the selection of anchor rod support types. SUMMARY

[0005] To solve the above technical problems, the present application provides a method and device for characterizing dynamic load response characteristics of an anchor rod, an equipment and a storage medium.

[0006] The present application provides a method for characterizing dynamic load response characteristics of an anchor rod, comprising:

[0007] Obtaining a pre-processed test data set, the pre-processed test data set comprising the rod body load and the cumulative elongation of the rod body collected in each impact during the process of the anchor rod being subjected to one or more impacts until breaking;

[0008] Calculating the pre-processed rod body load and the cumulative elongation of the rod body by a preset algorithm to obtain a target data set; wherein the preset algorithm comprises at least one of a moving average algorithm, a standard deviation algorithm, a confidence interval algorithm, a maximum and minimum value algorithm, and a quantile algorithm;

[0009] Taking the cumulative elongation of the rod body as the abscissa and the rod body load as the ordinate, a target curve corresponding to the target data set is drawn; the target curve is a dynamic load response characteristic curve of the anchor rod, and the dynamic load response characteristic curve comprises an upper envelope line, a lower envelope line and a median line.

[0010] According to the method for characterizing the dynamic load response characteristics of an anchor rod provided by the application, the preset algorithm comprises a moving average algorithm; the target data set is obtained by calculating the pretreated rod body load and the rod body cumulative elongation through the preset algorithm, and the method comprises the following steps:

[0011] The pretreated rod body cumulative elongation is sorted in ascending order to obtain a rod body cumulative elongation sequence;

[0012] For each sequence value in the rod body cumulative elongation sequence, the rod body load corresponding to each sequence value in a first preset sliding window corresponding to the sequence value is obtained to obtain a first rod body load set, and the average value of the first rod body load set is calculated;

[0013] The average value is offset by a preset percentage to obtain an upper limit and a lower limit of the average value of the rod body load corresponding to each sequence value;

[0014] All sequence values in the rod body cumulative elongation sequence and the corresponding upper limit and lower limit of the average value of the rod body load are taken as the target data set.

[0015] According to the method for characterizing the dynamic load response characteristics of an anchor rod provided by the application, the preset algorithm comprises a standard deviation algorithm; the target data set is obtained by calculating the pretreated rod body load and the rod body cumulative elongation through the preset algorithm, and the method comprises the following steps:

[0016] The pretreated rod body cumulative elongation is sorted in ascending order to obtain a rod body cumulative elongation sequence;

[0017] For each sequence value in the rod body cumulative elongation sequence, the rod body load corresponding to each sequence value in a second preset sliding window corresponding to the sequence value is obtained to obtain a second rod body load set, and the standard deviation of the second rod body load set is calculated;

[0018] The standard deviation is offset by a second preset percentage to obtain an upper limit and a lower limit of the standard deviation of the rod body load corresponding to each sequence value;

[0019] All sequence values in the rod body cumulative elongation sequence and the corresponding upper limit and lower limit of the standard deviation of the rod body load are taken as the target data set.

[0020] According to the method for characterizing the dynamic load response characteristics of an anchor rod provided by the application, the preset algorithm further comprises a confidence interval algorithm; the target data set is obtained by calculating the pretreated rod body load and the rod body cumulative elongation through the preset algorithm, and the method comprises the following steps:

[0021] obtaining a Z value based on a preset confidence level;

[0022] for each sequence value in the sequence of the cumulative elongation of the rod, calculating a rod load confidence interval upper limit and a rod load confidence interval lower limit corresponding to each sequence value based on the Z value and the standard deviation;

[0023] all sequence values in the sequence of the cumulative elongation of the rod and the corresponding rod load confidence interval upper limit and rod load confidence interval lower limit are taken as the target data set.

[0024] According to the method for characterizing the dynamic load response characteristics of the anchor rod provided by the application, the preset algorithm includes a maximum and minimum value algorithm; the calculation of the pretreated rod load and the cumulative elongation of the rod by the preset algorithm to obtain a target data set includes:

[0025] The pretreated cumulative elongation of the rod is sorted in ascending order to obtain a sequence of the cumulative elongation of the rod;

[0026] for each sequence value in the sequence of the cumulative elongation of the rod, obtaining the rod load corresponding to each sequence value in a third preset sliding window adjacent to the sequence value to obtain a third rod load set, and calculating the maximum value and the minimum value of the rod load in the third rod load set;

[0027] all sequence values in the sequence of the cumulative elongation of the rod and the corresponding maximum value and minimum value of the rod load are taken as the target value set.

[0028] According to the method for characterizing the dynamic load response characteristics of the anchor rod provided by the application, the preset algorithm includes a quantile algorithm; the calculation of the pretreated rod load and the cumulative elongation of the rod by the preset algorithm to obtain a target data set includes:

[0029] The pretreated cumulative elongation of the rod is sorted in ascending order to obtain a sequence of the cumulative elongation of the rod;

[0030] for each sequence value in the sequence of the cumulative elongation of the rod, a fourth preset sliding window is used to sequentially select a sub-window on the sequence of the cumulative elongation of the rod;

[0031] The maximum value and the minimum value of the rod load in each sub-window are calculated;

[0032] all sequence values in the sequence of the cumulative elongation of the rod and the corresponding maximum value and minimum value of the rod load are taken as the target value set.

[0033] The application further provides an anchor rod selection method, which comprises:

[0034] obtaining occurrence characteristic data of a coal mine stratum;

[0035] designing a roadway excavation scheme of a roadway section according to engineering needs, such as a roadway shape and geometric dimensions;

[0036] determining a bolt support scheme of the roadway according to the occurrence characteristic data and the roadway excavation scheme by using a corresponding support design method; the bolt support scheme includes a demand for dynamic parameters of the bolt;

[0037] matching the demand for dynamic parameters of the bolt with a dynamic load response characteristic curve of the bolt obtained by using the characterization method of the bolt dynamic load response characteristic, and determining optimal characteristic parameters of the bolt according to a matching result, so as to complete selection of the bolt.

[0038] The application further provides a characterization device of a bolt dynamic load response characteristic, and the device comprises:

[0039] a data acquisition module, configured to acquire a preprocessed test data set; the preprocessed test data set includes a rod body load and a rod body cumulative elongation of the bolt in each impact in a process in which the bolt is subjected to one or more impacts until fracture;

[0040] a data processing module, configured to calculate the preprocessed rod body load and the rod body cumulative elongation by using a preset algorithm to obtain a target data set; wherein the preset algorithm includes at least one of a moving average algorithm, a standard deviation algorithm, a confidence interval algorithm, a maximum and minimum value algorithm and a quantile algorithm;

[0041] a curve drawing module, configured to draw a target curve corresponding to the target data set by taking the rod body cumulative elongation as an abscissa and the rod body load as an ordinate; the target curve is a dynamic load response characteristic curve of the bolt, and the dynamic load response characteristic curve includes an upper envelope line, a lower envelope line and a median line.

[0042] The application further provides an electronic device, including a memory, a processor and a computer program stored in the memory and executable on the processor; when the processor executes the program, a characterization method of a bolt dynamic load response characteristic as described in any one of the above is realized.

[0043] The application further provides a non-transitory computer readable storage medium, which stores a computer program; when the computer program is executed by a processor, a characterization method of a bolt dynamic load response characteristic as described in any one of the above is realized.

[0044] The application further provides a computer program product, which includes a computer program; when the computer program is executed by a processor, a characterization method of a bolt dynamic load response characteristic as described in any one of the above is realized.

[0045] The application provides a method, device, equipment and storage medium for characterizing an anchor rod dynamic load response feature, which comprises the following steps: obtaining a pretreated test data set; the pretreated test data set comprises a rod body load and a rod body cumulative elongation of the anchor rod collected in each impact in a process in which the anchor rod is subjected to one or more impacts until the anchor rod is broken; calculating the pretreated rod body load and the pretreated rod body cumulative elongation by using a preset algorithm to obtain a target data set; wherein the preset algorithm comprises at least one of a moving average algorithm, a standard deviation algorithm, a confidence interval algorithm, a maximum and minimum value algorithm and a quantile algorithm; taking the rod body cumulative elongation as the abscissa and the rod body load as the ordinate to draw a target curve corresponding to the target data set; and the target curve comprises an upper envelope line, a lower envelope line and a median line. The application can obtain the target curve of the anchor rod dynamic load response feature, the trend and volatility of the test data can be intuitively reflected through the target curve, abnormal values can be detected, and key basic data support can be provided for anchor rod support selection, roadway support scheme optimization and roadway rock burst disaster prevention in a mine roadway. BRIEF DESCRIPTION OF DRAWINGS

[0046] In order to more clearly illustrate the technical solutions in the application or prior art, the following will briefly introduce the drawings needed to be used in the embodiments or prior art description. Obviously, the drawings in the following description are some embodiments of the application, and for those skilled in the art, other drawings can also be obtained without creative labor.

[0047] Figure 1 is one of the flowcharts of the method for characterizing the anchor rod dynamic load response feature provided by the application;

[0048] Figure 2 is another flowchart of the method for characterizing the anchor rod dynamic load response feature provided by the application;

[0049] Figure 3 is a schematic diagram of the anchor rod dynamic load response feature curve provided by the application;

[0050] Figure 4 is a structural schematic diagram of the device for characterizing the anchor rod dynamic load response feature provided by the application;

[0051] Figure 5 is a structural schematic diagram of the anchor rod selection device provided by the application;

[0052] Figure 6 is a structural schematic diagram of the electronic device provided by the application. DETAILED DESCRIPTION

[0053] In order to make the objects, technical solutions and advantages of the present application clearer, the technical solutions will be described clearly and completely below in conjunction with the accompanying drawings in the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative work fall within the protection scope of the present application.

[0054] It should be noted that, in the description of the embodiments of the present application, the terms "comprise", "contain" or any other variants thereof are intended to cover the non-exclusive inclusion, so that the process, method, article or equipment including a series of elements not only includes those elements, but also includes other elements not explicitly listed or inherent to such process, method, article or equipment. Without more limitation, the element defined by the statement "comprising a" does not exclude the presence of another identical element in the process, method, article or equipment comprising the element. The terms "upper", "lower" and the like indicate the orientation or positional relationship shown in the drawings, and are only for the convenience of describing the present application and simplifying the description, and do not indicate or imply that the devices or elements referred to must have a particular orientation, be constructed and operated in a particular orientation, and therefore cannot be understood as a limitation on the present application. Unless otherwise explicitly specified and limited, the terms "mount", "connect", "connect" should be understood broadly, for example, it can be fixedly connected, or it can be detachably connected, or integrally connected; it can be mechanically connected, or it can be electrically connected; it can be directly connected, or it can be indirectly connected through an intermediate medium, or it can be connected inside two elements. For those of ordinary skill in the art, the specific meaning of the above terms in the present application can be understood according to the specific circumstances.

[0055] The terms "first", "second", and the like in the present application are used to distinguish similar objects, and are not used to describe a specific order or sequence. It should be understood that the data thus used can be interchanged under appropriate circumstances, so that the embodiments of the present application can be implemented in an order other than those illustrated or described herein, and the objects distinguished by "first", "second", etc. are generally a class, and do not limit the number of objects, for example, the first object can be one or more. In addition, "and / or" means at least one of the connected objects, and the character " / ", generally means that the front and rear associated objects are in an "or" relationship.

[0056] The specific embodiments of the present application will be described below in conjunction with Figures 1-5 the accompanying drawings.

[0057] In one embodiment, as Figure 1 shown, a method for characterizing the dynamic load response characteristics of an anchor rod is provided, comprising the following steps:

[0058] Step 101, obtaining a pre-processed test data set; the pre-processed test data set includes the rod body load of the anchor rod and the cumulative elongation of the rod body in each impact collected in the process of the anchor rod being subjected to one or more impacts until breaking;

[0059] The rod body load refers to various direct actions applied on the engineering structure to make the engineering structure or component produce effects. The cumulative elongation of the rod body refers to the cumulative algebraic sum of the elongation of the rod body in the previous impact overloads.

[0060] Specifically, a certain form of anchor rod dynamic load impact test bench (such as a drop hammer impact test bench) is used to impact the anchor rod by one or more impacts with a certain energy (higher than the energy threshold of anchor rod impact deformation) until the anchor rod is finally broken. During each impact process, the load of the anchor rod body and the elongation of the rod body and other index data are collected in real time. The elongation of the rod body at the starting point of the subsequent impact is the elongation of the rod body after the end of the previous impact. The elongation increment of this impact is added to the starting point to obtain the dynamic cumulative elongation in this impact process. In this way, the load of the rod body and the cumulative elongation of the rod body in the previous impact process of the anchor rod are obtained. The data pair is used as a sample to obtain the data set of the dynamic load impact response of the anchor rod, which is used as a data sample for drawing the envelope curve in the next step.

[0061] Before the algorithm processing of the test data set, pre-processing is also needed. The pre-processing process includes data cleaning, missing value filling and other data preprocessing operations on the test data to ensure the reliability and integrity of the test data. The data cleaning includes removing interference data. The interference points in the original data are removed, and the effective data is retained to improve the accuracy of the analysis sample. This is because: in the process of anchor rod impact test, there are 2 parts of data that need to be removed in each impact process: one is that the impact load rises sharply from 0 kN to the yield load in an instant, but the value in the rising process is much lower than the load value in the yield and strengthening stage of the anchor rod; the second is that the rod body rebounds at the end of the impact, and the impact load caused by the rebound decreases sharply from more than 200 kN to less than 100 kN, and the value in the load steep drop process is much lower than the load value in the yield and strengthening stage of the anchor rod. If these data are retained, it will affect the accuracy of the analysis result and reduce the reliability of the envelope curve. Therefore, the abnormal load points need to be removed.

[0062] After removing the abnormal data points, the pre-processed test data set is obtained.

[0063] Step 102, calculating the rod body load and the cumulative elongation of the rod body by a preset algorithm to obtain a target data set; the preset algorithm includes at least one of a moving average algorithm, a standard deviation algorithm, a confidence interval algorithm, a maximum and minimum value algorithm, and a quantile algorithm.

[0064] Specifically, for the pre-processed test data, various data analysis algorithms such as moving average method, standard deviation method, maximum and minimum value method, confidence interval method, quantile method, etc. can be used to obtain the data set corresponding to the algorithm, and then the data set of the envelope line is obtained. The data set includes the data pair set required for the upper envelope line, the lower envelope line and the median line.

[0065] Step 103, taking the cumulative elongation of the rod as the horizontal coordinate and the load of the rod as the vertical coordinate, a target curve corresponding to the target data set is drawn; the target curve is the dynamic load response characteristic curve of the anchor rod, and the dynamic load response characteristic curve includes the upper envelope line, the lower envelope line and the median line.

[0066] Specifically, taking the cumulative elongation of the rod as the horizontal coordinate and the load of the rod as the vertical coordinate, a target curve corresponding to the target data set is drawn; the target curve is the dynamic load response characteristic curve of the anchor rod, and the dynamic load response characteristic curve includes the upper envelope line, the lower envelope line or the median line.

[0067] The above embodiment processes the pre-processed test data set by moving average algorithm, standard deviation algorithm, confidence interval algorithm, maximum and minimum value algorithm or quantile algorithm to obtain the target curve of the anchor rod dynamic load response characteristic. The trend and volatility of the test data can be intuitively reflected through the target curve, which helps to detect abnormal values and provides key basic data support for anchor support selection, roadway support scheme optimization and roadway rock burst disaster prevention in mine roadway.

[0068] In an embodiment, the above-mentioned preset algorithm includes a moving average algorithm; the above-mentioned step 102 includes: sorting the pre-processed cumulative elongation of the rod in ascending order to obtain a rod cumulative elongation sequence; for each sequence value in the rod cumulative elongation sequence, obtaining the rod load corresponding to each sequence value in the preset sliding window corresponding to the sequence value to obtain a first rod load set, and calculating the average value of the first rod load set; offsetting the average value by a preset percentage to obtain the upper limit and the lower limit of the average value of the rod load corresponding to each sequence value; taking all sequence values in the rod load sequence and their corresponding upper limit and lower limit of the average value of the rod load as the target data set.

[0069] The moving average algorithm is used to calculate the average value of the data set. In this algorithm, for any data point, a certain number of data points in its neighborhood window are taken to perform the average operation. Each data point in the data set is traversed to obtain a sequence of average values.

[0070] Specifically, the moving average algorithm in this embodiment mainly includes:

[0071] ① Local description result of scatter points: take the cumulative elongation of the rod as the horizontal coordinate, take the load of the rod as the vertical coordinate, and select a fixed number of points W for the sliding window points The pre-processed test data set is sorted in ascending order of the cumulative elongation of the rod (i.e. the horizontal coordinate). For the ith point, the vertical coordinate is described by the mean value of the vertical coordinates of the records from the (i-W)th to the ith record, and the description result of the vertical coordinate of any scatter point is obtained as follows:

[0072]

[0073] In the formula, is the mean value of the vertical coordinates of the neighborhood of the ith point in the sequence of the cumulative elongation of the rod; W points is the number of points in the sliding window; y i is the vertical coordinate value (i.e. the load of the rod) of the ith point in the original data sequence, and the corresponding horizontal coordinate is x i .

[0074] ② Deviation setting: in order to obtain the upper and lower envelope lines, the local description result of the scatter points needs to be offset. The moving average method adopts a percentage form, so the deviation percentage Deviation needs to be set.

[0075]

[0076]

[0077] In the formula, is the mean value of the vertical coordinates of the neighborhood of the ith point in the sequence of the cumulative elongation of the rod, u i is the upper limit of the average value of the load of the rod of the ith point, l i is the lower limit of the average value of the load of the rod of the ith point, and Deviation is the deviation percentage of the upper and lower limits.

[0078] In the above embodiment, the upper limit and the lower limit of the average value of the load of the rod corresponding to each cumulative elongation of the rod are calculated by the moving average algorithm, which is beneficial to provide a data basis for subsequent detection of rod abnormalities.

[0079] ​​In one embodiment, the preset algorithm includes a standard deviation algorithm, and step 102 includes: sorting the preprocessed cumulative elongation of the rod in ascending order to obtain a sequence of cumulative elongation of the rod; for each sequence value in the sequence of cumulative elongation of the rod, obtaining the rod load corresponding to each sequence value within a preset sliding window to obtain a second set of rod loads, and calculating the standard deviation of the second set of rod loads; deviating the standard deviation from the second preset percentage above and below to obtain the upper bound and lower bound of the standard deviation of the rod load corresponding to each sequence value; and using all sequence values ​​in the rod load sequence and their corresponding upper and lower bounds of the standard deviation of the rod load as the target dataset.

[0080] The standard deviation algorithm is a statistical algorithm that measures the dispersion of a data series and is used to analyze the dispersion and volatility of data. This algorithm calculates the mean of the dataset and the sum of the squares of the differences between each data point and the mean. Finally, it divides the sum of squares by the number of data points to obtain the standard deviation. A larger standard deviation indicates a more dispersed distribution of the dataset, and vice versa.

[0081] Specifically, the envelope of a data series can be obtained using the concept of standard deviation. The envelope of a data series can be calculated by appropriately adding or subtracting the mean and standard deviation of the data sequence. Specifically, the upper envelope of a data series can be obtained by adding a certain multiple of the standard deviation to the mean; this multiple needs to be determined based on the volatility of the data. Similarly, the lower envelope of a data series can be obtained by subtracting a certain multiple of the standard deviation from the mean. The mathematical description of the algorithm is as follows:

[0082] ① Scattered local description results: A sliding window with a fixed width is used to sort the original dataset by the magnitude of the x-axis values. For the k-th window, its y-axis uses the x-axis range of... arrive The mean of the ordinates of the scatter points within the window is used to describe the scatter points, thus obtaining the description of any scatter point.

[0083]

[0084]

[0085] In the formula, and Let mean(y) and standard deviation be the calculated mean and standard deviation of all scatter points of the calculated sequence within the k-th window, respectively. i ) represents the mean function, std(y i ) represents the standard deviation function.

[0086] ② Setting the deviation: In order to obtain the upper and lower bound envelopes, the local description results of the scatter points need to be offset. This algorithm uses a multiple of the standard deviation α to set the deviation.

[0087]

[0088]

[0089] wherein, is the vertical coordinate of the kth window of the upper bound sequence of the data set, is the vertical coordinate of the kth window of the lower bound sequence of the data set, and both of them correspond to the horizontal coordinate k·W width , and α is the multiple of the standard deviation of the upper and lower bounds.

[0090] In the above embodiment, the standard deviation algorithm is used to calculate the standard deviation upper bound and the standard deviation lower bound of the rod load corresponding to the cumulative elongation of each rod, which is beneficial to provide a data basis for subsequent detection of rod abnormalities.

[0091] In an embodiment, the preset algorithm further includes a confidence interval algorithm; and the step 102 further includes:

[0092] obtaining a Z value based on the preset confidence level; for each sequence value in the sequence of cumulative elongation of the rod, calculating the upper bound and the lower bound of the confidence interval of the rod load corresponding to the sequence value based on the Z value and the standard deviation; and taking all sequence values in the sequence of rod load and the upper bound and the lower bound of the confidence interval of the rod load corresponding to the sequence values as the target data set.

[0093] The confidence interval method is an interval estimation method in statistics, which uses an interval to estimate an unknown population parameter, such as population mean, population proportion, etc. The meaning of the confidence interval is that, in the case of repeated sampling, there is a certain probability (confidence level) to ensure that the interval contains the true population parameter, i.e. the probability that the interval statistical parameter can represent the whole statistical parameter.

[0094] Similar to the idea of obtaining the envelope line of the data series by the standard deviation method, the description of the scatter points in the interval includes two parts of mean and variance. The upper and lower envelope lines of the data series are replaced by the upper and lower bounds of the confidence interval of the scatter points in the window range W width The mathematical description of the algorithm is as follows:

[0095]

[0096]

[0097] wherein, is the vertical coordinate of the kth window of the upper bound sequence of the data set, is the vertical coordinate of the kth window of the lower bound sequence of the data set, and Z(β) represents the statistical Z value corresponding to the confidence level β.

[0098] The above embodiment, by confidence interval algorithm, the upper limit and the lower limit of the confidence interval of the rod load corresponding to the cumulative elongation of each rod body are calculated, which is beneficial to provide data basis for subsequent detection of abnormal anchor rod.

[0099] In an embodiment, the preset algorithm includes a maximum and minimum value algorithm; and the step 102 includes:

[0100] The pretreated cumulative elongation of the rod body is sorted in ascending order to obtain a rod body cumulative elongation sequence; for each sequence value in the rod body cumulative elongation sequence, the rod load corresponding to each sequence value in a third preset sliding window adjacent to the sequence value is obtained to obtain a third rod load set, and the maximum value and the minimum value of the rod load in the third rod load set are calculated; and all sequence values in the rod load sequence and the corresponding maximum value and minimum value of the rod load are taken as a target value set.

[0101] Specifically, in order to ensure that all sample points are included in the envelope line, the maximum and minimum value method respectively uses the maximum value and the minimum value of all sampling points in the specified sliding window point number W points or the window size W width as the upper envelope line and the lower envelope line, which has the advantages of simple calculation and ensuring that the envelope includes all scattered points. Correspondingly, since the sampling points are not statistically analyzed, the representative ability of the data is insufficient. If the sampling points are few, the envelope line obtained is insufficient in representativeness and accuracy, and the smoothness of the curve is insufficient. The mathematical description of the algorithm is as follows:

[0102] ①Sliding window point number W points Get: Select a fixed number of points for the sliding window, sort the original data according to the horizontal coordinates, and for the ith point, select to scatter points, and the upper and lower boundary description results of any scatter point can be obtained.

[0103]

[0104]

[0105] In the formula, u i is the vertical coordinate of the ith point of the upper boundary sequence of the calculation data set, and l i is the vertical coordinate of the ith point of the lower boundary sequence of the calculation data set, and max(…), min(…) represent the maximum value and minimum value functions, respectively.

[0106] ②Fixed window size W width Get: Respectively describe the horizontal coordinates as W width , 2W width , until nW widthupper boundary ordinate of the corresponding position and lower boundary ordinate

[0107]

[0108]

[0109] wherein, is the ordinate of the kth window of the upper boundary sequence of the data set, is the ordinate of the kth window of the lower boundary sequence of the data set.

[0110] In the above embodiment, the maximum and minimum values of the rod load corresponding to the cumulative elongation of each rod are calculated by the maximum and minimum value algorithm, which is beneficial to provide a data basis for subsequent detection of rod abnormalities.

[0111] In an embodiment, the preset algorithm includes a quantile algorithm; and the step 102 includes:

[0112] The pretreated cumulative elongation of the rod is sorted in ascending order to obtain a cumulative elongation sequence of the rod; for each sequence value in the cumulative elongation sequence of the rod, a fourth preset sliding window is used to sequentially select a sub-window on the cumulative elongation sequence of the rod; the maximum and minimum values of the rod load in each sub-window are calculated; and all sequence values in the cumulative elongation sequence of the rod and the corresponding maximum and minimum values of the rod load are taken as a target value set.

[0113] Specifically, the algorithm uses quantile as an index to describe the original data set, which divides the data into one-fourth, median, three-fourths, etc., to reflect the distribution of the data. According to the set upper and lower quantiles, it can be determined where a sample data is in the entire data set, so as to obtain the envelope line of the data. This algorithm is very useful when processing dynamic data, and can quickly understand the trend and change of the data. The mathematical description of the algorithm is as follows:

[0114] Finding method: sliding window with fixed width, sorting the original data by horizontal coordinates, for the kth window, get the scatter sequence with horizontal coordinates ranging from to and sort them by vertical coordinates from small to large, the horizontal and vertical coordinates of the new sequence are respectively denoted as There are M scatter points in total. The upper and lower boundaries of the envelope line are determined according to the following formula:

[0115]

[0116]

[0117] wherein, to calculate the ordinate of the kth window of the upper bound sequence of the data set, to calculate the ordinate of the kth window of the lower bound sequence of the data set, W width P is the window size, u P is the window size, l are the upper and lower quantiles respectively, M is the number of sampling points in the interval, k u are the upper and lower quantiles respectively, M is the number of sampling points in the interval, k l respectively represent the ranking index positions of the upper and lower quantiles.

[0118] The above embodiment is beneficial to provide a data basis for subsequent detection of anchor rod abnormalities by calculating the maximum value and the minimum value of the rod load corresponding to the cumulative elongation of each rod body through the quantile algorithm.

[0119] The application also provides an anchor rod selection method, comprising: obtaining occurrence characteristic data of a coal mine stratum; designing a roadway excavation scheme of a roadway section according to engineering needs; determining an anchor rod support scheme of the roadway and proposing a demand for dynamic parameters of the anchor rod by using a corresponding support design method according to the occurrence characteristic data of the coal mine stratum and the roadway excavation scheme; and matching the demand for dynamic parameters of the anchor rod with a dynamic load response characteristic curve of the anchor rod obtained by using the characterization method of the dynamic load response characteristic of the anchor rod as recorded in the above embodiment, and determining characteristic parameters of the anchor rod according to a matching result.

[0120] The occurrence characteristic data of the coal mine stratum refers to stratum conditions of a construction site, including rock strength, fissure development degree, ground stress, etc., which can be flexibly selected according to actual conditions, and is not limited in the present document.

[0121] Specifically, when construction work needs to be performed, for example, when an underground roadway needs to be constructed, the stratum conditions of the local area need to be detected first, and the occurrence characteristic data of the coal mine stratum is obtained through different measuring devices; a roadway excavation scheme of a roadway section is designed according to engineering needs, and the scheme includes information such as roadway shape and roadway geometric size; an anchor rod support scheme of the roadway is determined according to the occurrence characteristic data and the roadway excavation scheme, and a demand for dynamic parameters of the anchor rod is proposed; the demand for dynamic parameters of the anchor rod is compared and matched with the dynamic load response characteristic curve of the anchor rod, and finally the characteristic parameters of the required anchor rod are obtained, including anchor rod thickness, length, spacing, yield strength, etc.

[0122] The above embodiment obtains the dynamic load response characteristic curve of the anchor rod through prior tests, wherein the dynamic load response characteristic curves of different types of anchor rods are unique, and when engineering construction is needed, the stratum condition obtained through prior measurement is compared and matched with the dynamic load response characteristic curve of the anchor rod to finally obtain the optimal anchor rod characteristic parameter. Compared with the traditional method of determining the anchor rod characteristic parameter through actual tests, the method can improve the anchor rod selection efficiency, save the test cost of repeated tests, and provide an effective data basis for subsequent anchor rod maintenance and replacement.

[0123] The characterization device of the anchor rod dynamic load response characteristic provided by the present application is described below, and the characterization device of the anchor rod dynamic load response characteristic described below can be correspondingly referred to the characterization method of the anchor rod dynamic load response characteristic described above.

[0124] In one embodiment, as shown in Figure 4 A characterization device of an anchor rod dynamic load response characteristic is provided, comprising a data acquisition module 401, a data processing module 402, and a curve drawing module 403, wherein

[0125] The data acquisition module 401 is configured to acquire a preprocessed test data set, wherein the preprocessed test data set comprises the rod body load and the rod body cumulative elongation of the anchor rod in each impact during the process of the anchor rod being subjected to one or more impacts until the anchor rod is broken.

[0126] The data processing module 402 is configured to calculate the preprocessed rod body load and the rod body cumulative elongation by a preset algorithm to obtain a target data set, wherein the preset algorithm comprises at least one of a moving average algorithm, a standard deviation algorithm, a confidence interval algorithm, a maximum and minimum value algorithm, and a quantile algorithm.

[0127] The curve drawing module 403 draws a target curve corresponding to the target data set with the rod body cumulative elongation as the horizontal coordinate and the rod body load as the vertical coordinate, wherein the target curve comprises an upper envelope line, a lower envelope line, and a median line.

[0128] In one embodiment, the preset algorithm comprises a moving average algorithm, and the data processing module 402 is further configured to:

[0129] The pretreated rod cumulative elongation is sorted in ascending order to obtain a rod cumulative elongation sequence; for each sequence value in the rod cumulative elongation sequence, the rod load corresponding to each sequence value in a first preset sliding window corresponding to the sequence value is obtained to obtain a first rod load set, and the average value of the first rod load set is calculated; the average value is offset by a preset percentage to obtain a rod load average upper limit and a rod load average lower limit corresponding to each sequence value; and all sequence values in the rod cumulative elongation sequence and the corresponding rod load average upper limit and rod load average lower limit are taken as the target data set.

[0130] In one of the embodiments, the preset algorithm includes a standard deviation algorithm; and the data processing module 402 is further configured to:

[0131] The pretreated rod cumulative elongation is sorted in ascending order to obtain a rod cumulative elongation sequence; for each sequence value in the rod cumulative elongation sequence, the rod load corresponding to each sequence value in a second preset sliding window corresponding to the sequence value is obtained to obtain a second rod load set, and the standard deviation of the second rod load set is calculated; the standard deviation is offset by a second preset percentage to obtain a rod load standard deviation upper limit and a rod load standard deviation lower limit corresponding to each sequence value; and all sequence values in the rod cumulative elongation sequence and the corresponding rod load standard deviation upper limit and rod load standard deviation lower limit are taken as the target data set.

[0132] In one of the embodiments, the preset algorithm further includes a confidence interval algorithm; and the data processing module 402 is further configured to:

[0133] A Z value is obtained based on a preset confidence level; for each sequence value in the rod cumulative elongation sequence, a rod load confidence interval upper limit and a rod load confidence interval lower limit corresponding to each sequence value are calculated based on the Z value and the standard deviation; and all sequence values in the rod cumulative elongation sequence and the corresponding rod load confidence interval upper limit and rod load confidence interval lower limit are taken as the target data set.

[0134] In one of the embodiments, the preset algorithm includes a maximum and minimum value algorithm; and the data processing module 402 is further configured to:

[0135] The pretreated rod cumulative elongation is sorted in ascending order to obtain a rod cumulative elongation sequence; for each sequence value in the rod cumulative elongation sequence, the rod load corresponding to each sequence value in a third preset sliding window adjacent to the sequence value is obtained to obtain a third rod load set, and the maximum rod load and the minimum rod load in the third rod load set are calculated; and the entire sequence value in the rod cumulative elongation sequence and the corresponding maximum rod load and minimum rod load are taken as the target value set.

[0136] In one of the embodiments, the preset algorithm comprises a quantile algorithm; and the data processing module 402 is further configured to:

[0137] The pretreated rod cumulative elongation is sorted in ascending order to obtain a rod cumulative elongation sequence; for each sequence value in the rod cumulative elongation sequence, a sub-window is selected on the rod cumulative elongation sequence in turn by using a fourth preset sliding window; the maximum rod load and the minimum rod load in each sub-window are calculated; and the entire sequence value in the rod cumulative elongation sequence and the corresponding maximum rod load and minimum rod load are taken as the target value set.

[0138] In one of the embodiments, as shown in Figure 5 An anchor rod selection device is provided, comprising: a stratum condition data acquisition module 501, a roadway excavation scheme design module 502, and an anchor rod characteristic parameter determination module 503, wherein,

[0139] The stratum condition data acquisition module 501 is configured to acquire occurrence characteristic data of a coal mine stratum.

[0140] The roadway excavation scheme design module 502 is configured to determine an anchor rod support scheme of a roadway by using a corresponding support design method according to the occurrence characteristic data and the roadway excavation scheme; the anchor rod support scheme comprises a requirement for dynamic parameters of the anchor rod.

[0141] The anchor rod characteristic parameter determination module 503 is configured to match the requirement for dynamic parameters of the anchor rod with a dynamic load response characteristic curve of the anchor rod obtained by using the above-mentioned representation method of the dynamic load response characteristic of the anchor rod, and determine the best characteristic parameters of the anchor rod according to a matching result, thereby completing selection of the anchor rod.

[0142] Figure 6 An example of a schematic diagram of a physical structure of an electronic device is shown in Figure 6As shown, the electronic device can include a processor 610, a communications interface 620, a memory 630, and a communications bus 640, wherein the processor 610, the communications interface 620, and the memory 630 complete mutual communication through the communications bus 640. The processor 610 can invoke a logic instruction in the memory 630 to execute the method for characterizing the dynamic load response characteristics of the anchor rod, the method including: obtaining a preprocessed test data set; the preprocessed test data set including the rod body load and the rod body cumulative elongation of the anchor rod in each impact collected in the process of the anchor rod being subjected to one or more impacts until fracture; calculating the preprocessed rod body load and the rod body cumulative elongation by a preset algorithm to obtain a target data set; wherein the preset algorithm includes at least one of a moving average algorithm, a standard deviation algorithm, a confidence interval algorithm, a maximum and minimum value algorithm, and a quantile algorithm; drawing a target curve corresponding to the target data set with the rod body cumulative elongation as the abscissa and the rod body load as the ordinate; the target curve including an upper envelope line, a lower envelope line, and a median line.

[0143] In addition, the logic instruction in the memory 630 described above can be implemented in the form of a software functional unit and sold or used as an independent product, and can be stored in a computer readable storage medium. Based on such understanding, the technical solutions of the present application essentially or the part that contributes to the prior art or part of the technical solutions can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes several instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in various embodiments of the present application. The aforementioned storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM, Read-Only Memory), a random access memory (RAM, Random Access Memory), a magnetic disk or an optical disk, and various program code storage media.

[0144] In another aspect, the present application also provides a computer program product, which comprises a computer program, the computer program being stored on a non-transitory computer-readable storage medium and executable by a processor to enable a computer to perform the anchoring rod dynamic load response characteristic representation method or the anchoring rod selection method provided by the above-mentioned methods. The anchoring rod dynamic load response characteristic representation method comprises: obtaining a pre-processed test data set; the pre-processed test data set comprises the anchoring rod body load and the anchoring rod body cumulative elongation collected in each impact in the process of the anchoring rod being subjected to one or more impacts until the anchoring rod is broken; performing calculation on the pre-processed anchoring rod body load and the anchoring rod body cumulative elongation by using a preset algorithm to obtain a target data set; the preset algorithm comprises at least one of a moving average algorithm, a standard deviation algorithm, a confidence interval algorithm, a maximum and minimum value algorithm and a quantile algorithm; a target curve corresponding to the target data set is drawn by taking the anchoring rod body cumulative elongation as the abscissa and the anchoring rod body load as the ordinate; the target curve is an anchoring rod dynamic load response characteristic curve, and the anchoring rod dynamic load response characteristic curve comprises an upper envelope line, a lower envelope line and a median line. The anchoring rod selection method comprises: obtaining stratum condition data; matching the stratum condition data with the anchoring rod dynamic load response characteristic curve obtained by using the anchoring rod dynamic load response characteristic representation method, and determining the characteristic parameters of the anchoring rod according to the matching result.

[0145] In another aspect, the present application also provides a non-transitory computer-readable storage medium, which stores a computer program executable by a processor to perform the anchoring rod dynamic load response characteristic representation method or the anchoring rod selection method provided by the above-mentioned methods. The anchoring rod dynamic load response characteristic representation method comprises: obtaining a pre-processed test data set; the pre-processed test data set comprises the anchoring rod body load and the anchoring rod body cumulative elongation collected in each impact in the process of the anchoring rod being subjected to one or more impacts until the anchoring rod is broken; performing calculation on the pre-processed anchoring rod body load and the anchoring rod body cumulative elongation by using a preset algorithm to obtain a target data set; the preset algorithm comprises at least one of a moving average algorithm, a standard deviation algorithm, a confidence interval algorithm, a maximum and minimum value algorithm and a quantile algorithm; a target curve corresponding to the target data set is drawn by taking the anchoring rod body cumulative elongation as the abscissa and the anchoring rod body load as the ordinate; the target curve is an anchoring rod dynamic load response characteristic curve, and the anchoring rod dynamic load response characteristic curve comprises an upper envelope line, a lower envelope line and a median line. The anchoring rod selection method comprises: obtaining stratum condition data; matching the stratum condition data with the anchoring rod dynamic load response characteristic curve obtained by using the anchoring rod dynamic load response characteristic representation method, and determining the characteristic parameters of the anchoring rod according to the matching result.

[0146] The device embodiments described above are merely illustrative, wherein the units described as separate components can or can not be physically separate, and the components displayed as units can or can not be physical units, i.e., can be located in one place, or can be distributed to multiple network units. Part or all of the modules can be selected to achieve the purposes of the embodiments according to actual needs. Those skilled in the art can understand and implement without creative labor.

[0147] Through the description of the above embodiments, those skilled in the art can clearly understand that the embodiments can be realized by means of software and the necessary general hardware platform, and of course can also be realized by hardware. Based on such understanding, the above technical solutions can be embodied in the form of a software product, which can be stored in a computer readable storage medium, such as a ROM / RAM, a magnetic disk, an optical disk, etc., and includes a number of instructions to make a computer device (which can be a personal computer, a server, or a network device, etc.) execute the methods described in each embodiment or some parts of the embodiments.

[0148] Finally, it should be noted that: the above embodiments are only used to illustrate the technical solutions of the present application, and not to limit them; although the present application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that: it can still modify the technical solutions recorded in the foregoing embodiments, or make equivalent replacement to part of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims

1. A method for selecting an anchor rod, characterized in that, The method comprises the following steps: obtaining occurrence characteristic data of a coal mine stratum; designing a roadway excavation scheme of a roadway section according to engineering needs; determining an anchor rod support scheme of the roadway by using a corresponding support design method according to the occurrence characteristic data and the roadway excavation scheme; the anchor rod support scheme comprises dynamic parameter requirements of the anchor rod; matching the dynamic parameter requirements of the anchor rod with dynamic load response characteristic curves of the anchor rod obtained by a characterization method of the dynamic load response characteristic of the anchor rod, and determining optimal characteristic parameters of the anchor rod according to a matching result, so as to complete selection of the anchor rod; wherein the characterization method of the dynamic load response characteristic of the anchor rod comprises the following steps: obtaining a preprocessed test data set; the preprocessed test data set comprises a rod body load and a rod body cumulative elongation of the anchor rod in each impact collected in a process in which the anchor rod is subjected to multiple impacts until the anchor rod is broken; calculating the preprocessed rod body load and the rod body cumulative elongation by using a preset algorithm to obtain a target data set; drawing a target curve corresponding to the target data set by taking the rod body cumulative elongation as the horizontal coordinate and the rod body load as the vertical coordinate; the target curve is a dynamic load response characteristic curve of the anchor rod, and the dynamic load response characteristic curve comprises an upper envelope line, a lower envelope line and a median line; the preset algorithm is a moving average algorithm; at this time, the calculation of the preprocessed rod body load and the rod body cumulative elongation by using the preset algorithm to obtain the target data set comprises the following steps: sorting the preprocessed rod body cumulative elongation in ascending order to obtain a rod body cumulative elongation sequence; for each sequence value in the rod body cumulative elongation sequence, a first preset sliding window is used to sequentially select a sub-window on the rod body cumulative elongation sequence; calculating an average value of the rod body load corresponding to the sequence value in each sub-window; deviating the average value by a first preset percentage up and down respectively to obtain an upper limit and a lower limit of the average value of the rod body load in each sub-window; taking all sequence values in the rod body cumulative elongation sequence and the corresponding upper limit and lower limit of the average value of the rod body load as the target data set.

2. A method of anchoring rod selection according to claim 1, characterised in that, Alternatively, the preset algorithm is a moving standard deviation algorithm; at this time, the calculation of the preprocessed rod body load and the rod body cumulative elongation by using the preset algorithm to obtain the target data set comprises the following steps: sorting the preprocessed rod body cumulative elongation in ascending order to obtain a rod body cumulative elongation sequence; for each sequence value in the rod body cumulative elongation sequence, a second preset sliding window is used to sequentially select a sub-window on the rod body cumulative elongation sequence; calculating a first standard deviation of the rod body load corresponding to the sequence value in each sub-window; deviating the first standard deviation by a second preset percentage up and down respectively to obtain an upper limit and a lower limit of the standard deviation of the rod body load in each sub-window; taking all sequence values in the rod body cumulative elongation sequence and the corresponding upper limit and lower limit of the standard deviation of the rod body load as the target data set.

3. A method of anchoring rod selection according to claim 1, c h a r a c t e r i s e d i n that Alternatively, the preset algorithm is a moving confidence interval algorithm; At this time, the target data set is obtained by calculating the pretreated rod body load and the pretreated rod body cumulative elongation through a preset algorithm, including: The pretreated rod body cumulative elongation is sorted in ascending order to obtain a rod body cumulative elongation sequence; For each sequence value in the rod body cumulative elongation sequence, a fifth preset sliding window is used to sequentially select a sub-window on the rod body cumulative elongation sequence; The second standard deviation of the rod body load corresponding to the sequence value in each sub-window is calculated; The Z value is obtained based on the preset confidence level; The upper and lower confidence intervals of the rod body load in each sub-window are calculated based on the Z value and the second standard deviation; All sequence values in the rod body cumulative elongation sequence and their corresponding upper and lower confidence intervals of the rod body load are taken as the target data set.

4. The method according to claim 1, c h a r a c t e r i z e d in that, Alternatively, the preset algorithm is a moving maximum and minimum value algorithm; At this time, the target data set is obtained by calculating the pretreated rod body load and the pretreated rod body cumulative elongation through a preset algorithm, including: The pretreated rod body cumulative elongation is sorted in ascending order to obtain a rod body cumulative elongation sequence; For each sequence value in the rod body cumulative elongation sequence, a third preset sliding window is used to sequentially select a sub-window on the rod body cumulative elongation sequence; The maximum value and the minimum value of the rod body load corresponding to the sequence value in each sub-window are calculated; All sequence values in the rod body cumulative elongation sequence and their corresponding maximum and minimum values of the rod body load are taken as the target data set.

5. The method according to claim 1, c h a r a c t e r i z e d in that, Alternatively, the preset algorithm is a moving quantile algorithm; At this time, the target data set is obtained by calculating the pretreated rod body load and the pretreated rod body cumulative elongation through a preset algorithm, including: The pretreated rod body cumulative elongation is sorted in ascending order to obtain a rod body cumulative elongation sequence; For each sequence value in the rod body cumulative elongation sequence, a fourth preset sliding window is used to sequentially select a sub-window on the rod body cumulative elongation sequence; The preset upper quantile value and the preset lower quantile value of the rod body load corresponding to the sequence value in each sub-window are calculated; All sequence values in the rod body cumulative elongation sequence and their corresponding preset upper quantile values and preset lower quantile values of the rod body load are taken as the target data set.

6. An electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, The processor executes the program to implement the anchor rod selection method of any one of claims 1-5.

7. A non-transitory computer-readable storage medium having stored thereon a computer program, characterized in that, The computer program is executed by the processor to implement the anchor rod selection method of any one of claims 1-5.

Images