Intelligent analysis method for change of shaft force of anchor rod under dynamic load

Abstract

This application relates to the field of anchor bolt support engineering technology and discloses an intelligent analysis method for anchor bolt axial force changes under dynamic load. The method includes acquiring measured waveform files of several anchor bolts deployed on the support structure during blasting; filtering the measured waveform files to obtain an initial waveform file; inputting the initial waveform file into a preset anchor bolt axial force model, and combining it with the corresponding anchor bolt static load analysis result file, outputting the dynamic load response characteristics under vertical, horizontal, and bidirectional excitation through dynamic-static coupling numerical calculation; based on the dynamic load response characteristics, determining the direction of the maximum hazardous excitation for each support structure, providing a reference for optimizing support parameters. This application can improve the accuracy of anchor bolt axial force analysis under dynamic load.

Description

Technical Field

[0001] This application relates to the field of anchor bolt support engineering technology, and in particular to an intelligent analysis method for changes in axial force of anchor bolts under dynamic load. Background Technology

[0002] Large-diameter deep-hole mines typically employ a combination of bolts (cables), anchor mesh, and shotcrete for support. When the effectiveness of the bolts decreases, the support structure of the mine's exit roadways is easily damaged. Currently, research on bolt support theory is mainly based on static conditions, but bolt support in actual engineering projects is inevitably affected by dynamic loads.

[0003] Existing anchor bolt support systems typically employ simplified methods such as equivalent static load method or design dynamic load coefficient when considering the effects of dynamic loads. These methods often fail to accurately represent the real dynamic load environment, resulting in significant analytical errors in anchor bolt support systems under complex dynamic loads.

[0004] Regarding the aforementioned technologies, the inventors discovered that existing methods for analyzing the axial force of anchor bolts under dynamic loads suffer from low accuracy. Summary of the Invention

[0005] To improve the analysis accuracy of anchor bolt axial force under dynamic load, this application provides an intelligent analysis method for the change of anchor bolt axial force under dynamic load.

[0006] Firstly, this application provides an intelligent analysis method for the change of axial force in anchor bolts under dynamic load.

[0007] This application is achieved through the following technical solution:

[0008] An intelligent analysis method for changes in axial force of anchor bolts under dynamic load includes the following steps: Obtain the measured waveform files of several anchor bolts deployed on the support structure during blasting; The measured waveform file is filtered to obtain an initial waveform file; The initial waveform file is input into the preset anchor bolt axial force model. Combined with the corresponding anchor bolt static load analysis result file, the dynamic load response characteristics under vertical excitation, horizontal excitation and bidirectional excitation are output through dynamic-static coupling numerical calculation. Based on the dynamic load response characteristics, the direction of the maximum harmful vibration for each support structure is determined, providing a reference for the optimization of support parameters.

[0009] In a preferred embodiment, this application can be further configured such that the steps for constructing the anchor bolt axial force model include: Determine the size of the model and define the range as 3 times the orifice radius from the blast hole as the simulation range; Set the minimum center frequency and critical damping ratio of the model; Set static dynamic boundaries for the model, and set independent dampers in the normal and tangential directions of the model boundaries.

[0010] In a preferred embodiment, this application can be further configured such that the step of constructing the anchor bolt axial force model also includes, The maximum mesh size of the control model is within 8m, and the maximum mesh size is less than or equal to 1 / 8 to 1 / 10 of the wavelength corresponding to the highest frequency of the input waveform.

[0011] In a preferred embodiment, this application can be further configured such that the training steps of the anchor bolt axial force model include, Set up monitoring points for the model's dynamic response; First, a gravity field is applied to the model to achieve static stress equilibrium, forming an initial static stress field; Dynamic loads are applied inside the boundaries or nodes of the model. The types of dynamic loads include acceleration time history, velocity time history, pressure time history, or concentrated force time history. The dynamic loads are respectively excited vertically, horizontally, and bidirectionally. Through dynamic-static coupling calculation, the dynamic load response characteristics, including dynamic displacement, peak axial force, and volume of the plastic zone of the structure, are extracted from the monitoring points.

[0012] In a preferred embodiment, this application can be further configured as follows: the step of inputting the initial waveform file into a preset anchor bolt axial force model, combining it with the corresponding anchor bolt static load analysis result file, and outputting the dynamic load response characteristics under vertical excitation, horizontal excitation, and bidirectional excitation through dynamic-static coupling numerical calculation includes... The static load file command is used to retrieve the static load analysis result file of the monitoring points; The waveform file is read into the measured waveform file collected by the monitoring point using the waveform file read command; The waveform file is loaded onto the static load analysis result file using the waveform file loading command; The dynamic load response results are obtained by solving the problem.

[0013] In a preferred embodiment, this application can be further configured such that the step of filtering the measured waveform file to obtain an initial waveform file includes: Determine the transverse wave velocity of the measured waveform file; Based on the shear wave velocity and the preset filtering frequency, the filtering wavelength corresponding to the filtering frequency is calculated. Calculate the grid size range based on the filtered wavelength; The maximum value of the grid size range is compared with a preset maximum grid size threshold; When the maximum value of the grid size range is greater than the maximum grid size threshold, the filtering frequency is used as the upper limit value; The initial waveform file is obtained by filtering out waveforms with frequencies greater than the upper limit from the measured waveform file.

[0014] Secondly, this application provides an intelligent analysis device for the change of axial force of anchor bolts under dynamic load.

[0015] This application is achieved through the following technical solution: An intelligent analysis device for analyzing the axial force variation of anchor bolts under dynamic load, comprising, The data module is used to acquire the measured waveform files generated during blasting of several anchor bolts deployed on the support structure. The filtering module is used to filter the measured waveform file to obtain an initial waveform file. The dynamic load response module is used to input the initial waveform file into the preset anchor bolt axial force model, combine it with the corresponding anchor bolt static load analysis result file, and output the dynamic load response characteristics under vertical excitation, horizontal excitation and bidirectional excitation through dynamic-static coupling numerical calculation. The analysis module is used to determine the direction of maximum harmful vibration for each support structure based on the dynamic load response characteristics, providing a reference for the optimization of support parameters.

[0016] Thirdly, this application provides a computer device.

[0017] This application is achieved through the following technical solution: A computer device includes a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements the steps of any of the above-described intelligent analysis methods for the change of axial force of anchor bolts under dynamic load.

[0018] Fourthly, this application provides a computer-readable storage medium.

[0019] This application is achieved through the following technical solution: A computer-readable storage medium storing a computer program, which, when executed by a processor, implements the steps of any of the above-described intelligent analysis methods for changes in anchor axial force under dynamic load.

[0020] Fifthly, this application provides a computer program product.

[0021] This application is achieved through the following technical solution: A computer program product includes a computer program that, when executed by a processor, implements the steps of any of the above-mentioned intelligent analysis methods for changes in axial force of anchor bolts under dynamic load.

[0022] In summary, compared with the prior art, the beneficial effects of the technical solution provided in this application include at least the following: Measured waveform files of several anchor bolts deployed on the support structure during blasting were obtained to realistically simulate the response of the support structure under dynamic blasting loads, improving the effectiveness of the initial data and making the subsequent anchor bolt axial force analysis results more reliable. The measured waveform files were filtered to obtain initial waveform files, improving waveform quality while reducing the highest frequency and lengthening the corresponding minimum wavelength, thus increasing the maximum mesh size of the model. This reduced the total number of mesh elements while keeping the original model dimensions unchanged, which helps shorten the model's computation time. The initial waveform files were input into a pre-defined anchor bolt axial force model, combined with the corresponding anchor bolt static load analysis results, and through dynamic-static coupling numerical calculations, the vertical excitation and horizontal... The dynamic load response characteristics under vibration and bidirectional vibration enable coupled simulation of the actual static initial state and complex dynamic load path. This reveals the true mechanical behavior of the anchor bolt evolving from static equilibrium to dynamic axial force response under dynamic load, resulting in more physically realistic output results. Furthermore, the model extracts dynamic load response characteristics under vertical, horizontal, and bidirectional vibration, which helps to clarify the differentiated influence of dynamic load directionality on the anchor bolt force mechanism. The model's output results are more accurate and comprehensive, improving the analysis accuracy of anchor bolt axial force under dynamic load. Based on the dynamic load response characteristics, the maximum hazardous vibration direction of each support structure can be determined. The analysis results of anchor bolt axial force changes under dynamic load can be used to predict and locate risk points in the support system in advance, providing a reference for the optimization of support parameters. Attached Figure Description

[0023] Figure 1 This is a schematic diagram of the overall process of an intelligent analysis method for the axial force variation of anchor bolts under dynamic load, provided as an exemplary embodiment of this application.

[0024] Figure 2 A schematic diagram of a mining area for an exemplary embodiment of this application, illustrating an intelligent analysis method for the axial force variation of anchor bolts under dynamic load.

[0025] Figure 3 The waveform diagram before and after filtering is provided for an exemplary embodiment of this application of an intelligent analysis method for dynamic load axial force variation of anchor bolts.

[0026] Figure 4 This application provides an exemplary embodiment of an intelligent analysis method for the change of axial force of anchor bolts under dynamic load, showing the distribution of the plastic zone of the surrounding rock after different excitation effects.

[0027] Figure 5 This application provides an exemplary embodiment of an intelligent analysis method for dynamic load-bearing anchor axial force variation, showing a distribution map of measuring points in the surrounding rock of a mining area.

[0028] Figure 6 The present application provides a structural block diagram of an intelligent analysis device for dynamic load axial force variation of anchor bolts, which is an exemplary embodiment of this application. Detailed Implementation

[0029] This specific embodiment is merely an explanation of this application and is not intended to limit it. After reading this specification, those skilled in the art can make modifications to this embodiment without contributing any inventive step, but such modifications are protected by patent law as long as they fall within the scope of the claims of this application.

[0030] To make the objectives, technical solutions, and advantages of the embodiments of this application clearer, the technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0031] The embodiments of this application will now be described in further detail with reference to the accompanying drawings.

[0032] Reference Figure 1 This application provides an intelligent analysis method for the change of axial force of anchor bolts under dynamic load. The main steps of the method are described below.

[0033] S1: Obtain the measured waveform files of several anchor bolts deployed on the support structure during blasting; S2: Filter the measured waveform file to obtain an initial waveform file; S3: Input the initial waveform file into the preset anchor bolt axial force model, combine it with the corresponding anchor bolt static load analysis result file, and output the dynamic load response characteristics under vertical excitation, horizontal excitation and bidirectional excitation through dynamic-static coupling numerical calculation; S4: Based on the dynamic load response characteristics, determine the maximum hazardous vibration direction of each support structure to provide a reference for the optimization of support parameters.

[0034] Specifically, to obtain the response characteristics of the support structure under actual dynamic blasting load, the model adopts the method of inputting measured waveform files. These measured waveform files can be velocity time history curve files. The testing instrument is the EXP3850 blasting vibration recorder and analysis system manufactured by Zhongke Dynamic Instruments Co., Ltd. By collecting the dynamic wave waveforms generated by the blasting, the system obtains measured waveform files of several anchor bolts deployed on the support structure during blasting, improving the validity of the initial data and making the subsequent analysis results of the anchor bolt axial force more reliable.

[0035] Next, the measured waveform file is filtered to obtain the initial waveform file. By setting an upper limit value for the frequency, waveforms that do not meet the conditions are filtered out. This improves the quality of the waveform file, reduces the highest frequency in the signal, increases the corresponding minimum wavelength, and thus increases the maximum size of the model's mesh. Under the condition that the original scale of the model remains unchanged, the total number of mesh units is reduced, which helps to shorten the computation time of the model.

[0036] The initial waveform file is input into a preset anchor bolt axial force model. Combined with the corresponding anchor bolt static load analysis result file, dynamic load response characteristics under vertical, horizontal, and bidirectional excitation are output through dynamic-static coupling numerical calculation. In practice, existing software such as FLAC3D can be used to construct the anchor bolt axial force model. The pre-processed waveform file is input into the anchor bolt axial force model, and combined with the corresponding anchor bolt static load analysis result file, dynamic load response characteristics under vertical, horizontal, and bidirectional excitation are output through dynamic-static coupling numerical calculation. The anchor bolt static load analysis result file is obtained by performing static load analysis on the model.

[0037] Preferably, the numerical calculation of dynamic-static coupling specifically includes: 1. Call the static load initial state The static load file call command is used to read and restore the static stress balance result file of the anchor monitoring point to form the initial static stress field, which serves as the initial mechanical state for dynamic load calculation.

[0038] 2. Read in the measured dynamic carrier waveform The filtered measured dynamic carrier waveform of the blasting was read into the model using the waveform file reading command to form an applicable dynamic load time history curve.

[0039] 3. Apply triaxial excitation load in stages Dynamic loads are applied at the model boundary or internal nodes in three directions: vertical, horizontal, and bidirectional superposition. The dynamic load type is one of the following: velocity time history, acceleration time history, pressure time history, or concentrated force time history.

[0040] 4. Dynamic-static coupling iterative solution Based on the initial static stress field, the dynamic carrier curve is dynamically superimposed on the static stress state, and dynamic time history iterative calculation is performed to make the anchor bolt axial force continuously evolve from the static load equilibrium state to the dynamic load response state.

[0041] 5. Extract dynamic load response features Four dynamic load response characteristics are extracted and output from the preset monitoring points: dynamic stress, dynamic displacement, peak value of anchor bolt axial force, and volume of plastic zone of surrounding rock.

[0042] 6. Output the comparison results of triaxial excitation. The dynamic load response results under vertical excitation, horizontal excitation, and bidirectional excitation are output respectively to form a three-dimensional comparison dataset for subsequent determination of the excitation direction of the maximum hazard.

[0043] Finally, based on the dynamic load response characteristics, the analysis results of the change in anchor axial force under dynamic load are used to predict and locate the risk points of the support system in advance, determine the maximum harmful vibration direction of each support structure, and provide a reference for the optimization of support parameters.

[0044] Preferably, the direction of maximum hazard excitation includes: based on the dynamic load response characteristics obtained from the dynamic-static coupling numerical calculation, three types of dynamic load response indicators are extracted from the anchor bolts under vertical excitation, horizontal excitation, and bidirectional excitation, namely, the peak axial force, dynamic displacement amplitude, and volume of the plastic zone of the surrounding rock. The above three types of indicators are comprehensively quantified and scored according to preset weights. Based on the scoring results, the degree of hazard of each support structure under different excitation directions is automatically determined, thereby locating the direction of maximum hazard excitation and providing a quantitative basis for the optimization of support parameters.

[0045] In one embodiment, the steps for constructing the anchor bolt axial force model include: First, determine the overall dimensions of the anchor bolt axial force model, and based on the principle of stress concentration at the borehole opening, set the area three times the borehole opening radius from the blasting hole as the core simulation area of ​​the model in order to accurately capture the real stress response of the surrounding rock and anchor bolt near the borehole opening under the dynamic load of blasting. Set the minimum center frequency and critical damping ratio of the model; In terms of model boundary treatment, in order to reduce the reflection interference of boundary waves during dynamic calculation, a static dynamic boundary is set for the model, and independent dampers are set in the normal and tangential directions of the model boundary, so that the normal damping and tangential damping are independent and act separately. This can efficiently absorb incident waves at different incident angles, significantly reduce the impact of boundary reflection waves on calculation accuracy, and ensure the authenticity and stability of the anchor axial force response under dynamic load.

[0046] Reference Figure 2 For example, in a mining area, the ore body height is 50m, and the total width of the three stops and three pillars is 108m. A 400m×200m model is established with twice the width of this as the center. The shotcrete and anchor support in roadway 2 is affected by the redistribution of stress fields after the mining outage of stops 1 and 3. In the analysis, if roadway 2 is taken as the research object, the radial distance of the anchor bolts is 1m, the model thickness is set to 1m, and the model thickness direction is taken to be consistent with the axial direction of the roadway, simplifying the model into a plane strain problem.

[0047] According to the principle of stress concentration at the orifice in elasticity, creating a hole in an elastic medium will cause the stress near the orifice to be much greater than the stress without a hole, and the stress concentration zone is localized. In this embodiment, the simulation range of the anchor bolt axial force model is set within a range of 3 times the orifice radius from the blast hole.

[0048] Rayleigh damping is set for the anchor bolt axial force model, including the minimum center frequency and critical damping ratio. For geotechnical materials, when using an elasto-plastic model for dynamic analysis, a significant amount of energy is dissipated during the plastic flow stage. Therefore, in large-strain dynamic analysis, only a small critical damping ratio needs to be set. After reaching the plastic zone, as the stress-strain hysteresis loop expands, energy dissipation becomes increasingly apparent. In this embodiment, the critical damping ratio can be set to 0.02, and the minimum center frequency is 250Hz, to determine the mechanical damping that meets the actual application scenario for the model.

[0049] To reduce computational complexity and minimize energy reflection at the model boundaries, a static dynamic boundary is applied to the anchor bolt axial force model. Independent dampers, attached independently to the normal and tangential directions of the model boundaries, effectively absorb body waves with incident angles greater than 30°. Verification shows that, in specific scenarios involving underground deposits with deep mining tunnels and stopes / pillars, applying a static dynamic boundary to the anchor bolt axial force model effectively absorbs reflected waves from all six faces of the model, achieving good dynamic load simulation results.

[0050] In addition, structural damping must be set for anchor bolt elements; otherwise, the damping of the anchor bolt elements and the soil and rock elements will be inconsistent, affecting the calculation accuracy of the model.

[0051] In one embodiment, the step of constructing the anchor bolt axial force model further includes, The maximum mesh size of the control model is within 8m, and the maximum mesh size is less than or equal to 1 / 8 to 1 / 10 of the wavelength corresponding to the highest frequency of the input waveform.

[0052] Since there is an upper limit frequency for energy propagation, the calculation results of the model can only match the actual situation when the input load is less than the upper limit frequency. Therefore, in the process of meshing, in order to ensure the uniformity of the mesh and control the maximum size of the mesh, the maximum size of the mesh is controlled within 8m to improve the calculation accuracy of the model.

[0053] Furthermore, since the frequency of the input waveform directly affects the propagation of waves in the soil and rock mass, the highest frequency of the input load is related to the maximum size of the model element. That is, the frequency of the input waveform affects the accuracy of the model. Therefore, the maximum mesh size is set to be less than or equal to 1 / 8 to 1 / 10 of the wavelength corresponding to the highest frequency of the input waveform.

[0054] In one embodiment, the training steps of the anchor bolt axial force model include, Set up monitoring points for the model's dynamic response; First, a gravity field is applied to the model to achieve static stress equilibrium, forming an initial static stress field; Dynamic loads are applied inside the boundaries or nodes of the model. The types of dynamic loads include acceleration time history, velocity time history, pressure time history, or concentrated force time history. The dynamic loads are respectively excited vertically, horizontally, and bidirectionally. Through dynamic-static coupling calculation, the dynamic load response characteristics, including dynamic displacement, peak axial force, and volume of the plastic zone of the structure, are extracted from the monitoring points.

[0055] Specifically, when training the anchor bolt axial force model, monitoring points for the model's dynamic response are set at intervals within the pre-defined simulation range of the anchor bolt axial force model. First, a gravity field is applied to the model to achieve static stress equilibrium, forming an initial static stress field; then, dynamic loads are applied inside the model's boundaries or nodes. In this embodiment, when analyzing the change in anchor bolt axial force under dynamic load, taking roadway 2 as the research object, a blasting dynamic load is applied to the anchor bolt axial force model to analyze the change in anchor bolt axial force in roadway 2 after the mining of stopes 1 and 3. The dynamic load is a velocity-time history, and excitation is performed vertically, horizontally, and in a superimposed manner in both vertical and horizontal directions.

[0056] In one embodiment, the step of inputting the initial waveform file into a preset anchor bolt axial force model, combining it with the corresponding anchor bolt static load analysis result file, and outputting the dynamic load response characteristics under vertical excitation, horizontal excitation, and bidirectional excitation through dynamic-static coupling numerical calculation includes: The static load file command is used to retrieve the static load analysis result file of the monitoring points; The waveform file is read into the measured waveform file collected by the monitoring point using the waveform file read command; The waveform file is loaded onto the static load analysis result file using the waveform file loading command; The dynamic load response results are obtained by solving the problem.

[0057] Specifically, when performing dynamic-static coupling calculations on the anchor bolt axial force model, the static load analysis result file of the monitoring points is called using the `restore` command, the waveforms collected by the monitoring points are read using the `table` command, and the dynamic load waveforms are loaded using the `apply` command. After solving, the dynamic load response result is obtained and saved as a result file with the .sav extension. The dynamic load response results extract dynamic load response characteristics including dynamic displacement with dynamic stress, peak axial force, and the volume of the structure's plastic zone.

[0058] In one embodiment, the step of filtering the measured waveform file to obtain an initial waveform file includes: Determine the transverse wave velocity; Based on the shear wave velocity and the preset filtering frequency, the filtering wavelength corresponding to the filtering frequency is calculated. Calculate the grid size range based on the filtered wavelength; The maximum value of the grid size range is compared with a preset maximum grid size threshold; When the maximum value of the grid size range is greater than the maximum grid size threshold, the filtering frequency is used as the upper limit value; The initial waveform file is obtained by filtering out waveforms with frequencies greater than the upper limit from the measured waveform file.

[0059] Specifically, to ensure that the mesh size in the numerical calculation meets the requirements of the dynamic calculation, for example, if the upper limit of the frequency is taken as 120Hz, it is verified whether the wavelength corresponding to the highest frequency waveform is greater than the maximum size of the model mesh.

[0060] Given the longitudinal wave velocity ; transverse wave velocity ; In the formula, For longitudinal wave velocity; The transverse wave velocity; R is Lamé constant; G is shear modulus; Density; is Poisson's ratio; E is the elastic modulus.

[0061] Since the wavelength of the transverse wave is shorter than that of the longitudinal wave, after calculating the transverse wave velocity, the maximum frequency is taken as 120Hz, and the wavelength corresponding to the highest frequency is approximately 86m. The maximum mesh size of the model is 8m, meaning that the upper limit of the filtering frequency is 120Hz, which is a reasonable value.

[0062] In this embodiment, the upper limit of the filtering frequency is selected as 120Hz to filter the input dynamic carrier waveform and filter out waveform files with a frequency greater than 120Hz.

[0063] Reference Figure 3 For example, (a) is the horizontal moving carrier pattern before filtering, (b) is the horizontal moving carrier pattern after filtering, (c) is the vertical moving carrier pattern before filtering, and (d) is the vertical moving carrier pattern after filtering.

[0064] This application improves the quality of the measured waveform file by filtering it, thereby reducing the highest frequency in the signal and increasing the corresponding minimum wavelength. As a result, the maximum size of the model's mesh increases, and the total number of mesh cells is reduced while keeping the original scale of the model unchanged, which helps to shorten the computation time of the model.

[0065] By analyzing the results file of the dynamic load response, the variation law of the anchor axial force under dynamic load can be obtained, thereby determining the direction of the maximum harmful vibration of each support structure and guiding the optimization of support parameters.

[0066] For example, refer to Figure 4This study illustrates the distribution of the plastic zone in the surrounding rock under static load, vertical excitation, horizontal excitation, and bidirectional excitation. Comparing the plastic zones under these three conditions, the plastic zone is smallest under horizontal excitation, similar to the distribution under static load; the plastic zone increases under vertical excitation, and is largest under bidirectional excitation. The changes in the plastic zone under these three conditions are mainly observed in the tunnel floor and the left side of the arch. The largest plastic zone under bidirectional excitation indicates that there is coupling and superposition of dynamic loads, leading to accelerated damage to the surrounding rock.

[0067] For example, in a mining area, the main monitoring focus is on the surrounding rock in the middle of the vertical wall of the roadway and the arch. Five monitoring points (1-5) are set up, and the distribution of the monitoring points is as follows: Figure 5 As shown.

[0068] Under vertical vibration, the displacement range of each measuring point in the Y direction is -0.27mm to 0.6mm. The displacement of the measuring points on the left and top of the roadway is positive, while the displacement of the measuring points on the right side of the roadway is negative, indicating that the surrounding rock is moving towards the center of the roadway. The displacement range of each measuring point in the Z direction is -1.02mm to 1.84mm. The maximum displacement in the Z direction occurs at measuring point 4, and the minimum displacement occurs at measuring point 5. The right side of the roadway has been mined out, which releases a certain amount of surrounding rock. Under dynamic load, it is easier to generate a larger displacement. Therefore, the maximum displacement occurs at measuring points 4 and 5 on the right side of the roadway.

[0069] Under horizontal excitation, the displacement range of each measuring point in the Y direction is -0.39mm to 0.48mm. The displacement range of each measuring point in the Z direction is -0.1mm to 0.1mm. The overall displacement is very small, with only the displacement of measuring point 4 being positive, while the displacements of all other measuring points are negative.

[0070] Under bidirectional excitation, the peak negative displacement in the Y direction at each measuring point is -0.40 mm, located at measuring point 1; the peak positive displacement is 0.41 mm, located at measuring point 4. The peak negative displacement in the Z direction at each measuring point is -1.02 mm, located at measuring point 4; the peak positive displacement is 1.86 mm, located at measuring point 5.

[0071] For example, in the analysis results, if the peak value of the axial force of the dynamic stress in the observation area exceeds the upper limit, the observation area can be focused on for observation, thereby enabling early prediction of the risk points of the positioning support system and providing a reference for the optimization of support parameters.

[0072] In summary, an intelligent analysis method for dynamic anchor axial force variation acquires measured waveform files of several anchors deployed on the support structure during blasting to realistically simulate the response law of the support structure under dynamic blasting load, improving the effectiveness of the initial data and making the subsequent anchor axial force analysis results more reliable. The measured waveform files are filtered to obtain the initial waveform file, which improves the quality of the waveform file while reducing the highest frequency in the signal and increasing the corresponding minimum wavelength, thereby increasing the maximum mesh size of the model. While keeping the original model scale unchanged, the total number of mesh elements is reduced, which helps to shorten the model's computation time. The initial waveform file is input into a preset anchor axial force model and combined with the corresponding anchor static load analysis result file, through dynamic-static coupling numerical analysis... The model calculates and outputs the dynamic load response characteristics under vertical, horizontal, and bidirectional excitation, achieving coupled simulation of the actual static initial state and complex dynamic load path. It reveals the true mechanical behavior of the anchor bolt evolving from static equilibrium to dynamic axial force response under dynamic load, resulting in more physically realistic outputs. Furthermore, the model extracts the dynamic load response characteristics under vertical, horizontal, and bidirectional excitations, which helps clarify the differentiated influence of dynamic load directionality on the anchor bolt's force mechanism. The model's outputs are more accurate and comprehensive, improving the analysis accuracy of anchor bolt axial force under dynamic load. Based on the dynamic load response characteristics, the maximum hazardous excitation direction for each support structure is determined. The analysis results of anchor bolt axial force changes under dynamic load can be used to predict and locate risk points in the support system in advance, providing a reference for optimizing support parameters.

[0073] It should be understood that the sequence number of each step in the above embodiments does not imply the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of this application.

[0074] Reference Figure 6 This application also provides an intelligent analysis device for dynamic load anchor bolt axial force variation, which corresponds one-to-one with the intelligent analysis method for dynamic load anchor bolt axial force variation in the above embodiments. The intelligent analysis device for dynamic load anchor bolt axial force variation includes... The data module is used to acquire the measured waveform files generated during blasting of several anchor bolts deployed on the support structure. The filtering module is used to filter the measured waveform file to obtain an initial waveform file. The dynamic load response module is used to input the initial waveform file into the preset anchor bolt axial force model, combine it with the corresponding anchor bolt static load analysis result file, and output the dynamic load response characteristics under vertical excitation, horizontal excitation and bidirectional excitation through dynamic-static coupling numerical calculation. The analysis module is used to determine the direction of maximum harmful vibration for each support structure based on the dynamic load response characteristics, providing a reference for the optimization of support parameters.

[0075] For specific limitations regarding an intelligent analysis device for changes in axial force of anchor bolts under dynamic load, please refer to the limitations of an intelligent analysis method for changes in axial force of anchor bolts under dynamic load mentioned above, which will not be repeated here.

[0076] The various modules in the aforementioned intelligent analysis device for anchor bolt axial force changes under dynamic load can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in or independent of the processor in a computer device, or stored in the memory of a computer device as software, so that the processor can call and execute the corresponding operations of each module.

[0077] In one embodiment, a computer device is provided, which may be a server. The computer device includes a processor, memory, a network interface, and a database connected via a system bus. The processor provides computing and control capabilities. The memory includes a non-volatile storage medium and internal memory. The non-volatile storage medium stores an operating system, computer programs, and a database. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage medium. The network interface is used to communicate with external terminals via a network connection. When the computer program is executed by the processor, it implements any of the above-described intelligent analysis methods for anchor bolt axial force changes under dynamic load.

[0078] In one embodiment, a computer-readable storage medium is provided, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements any of the above-described intelligent analysis methods for changes in anchor axial force under dynamic load.

[0079] In one embodiment, a computer program product is provided, which includes a computer program that, when executed by a processor, implements any of the above-described intelligent analysis methods for changes in anchor axial force under dynamic load.

[0080] Those skilled in the art will understand that all or part of the processes in the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this application. Any references to memory, storage, database, or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory. Non-volatile memory may include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory may include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), dual data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), RAMbus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

[0081] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the above-described division of functional units and modules is used as an example. In practical applications, the above functions can be assigned to different functional units and modules as needed, that is, the internal structure of the system can be divided into different functional units or modules to complete all or part of the functions described above.

Claims

1. A method for intelligent analysis of changes in shaft force of a jacked anchor rod under dynamic load, characterized in that, Includes the following steps, Obtain the measured waveform files of several anchor bolts deployed on the support structure during blasting; The measured waveform file is filtered to obtain an initial waveform file; The initial waveform file is input into the preset anchor bolt axial force model. Combined with the corresponding anchor bolt static load analysis result file, the dynamic load response characteristics under vertical excitation, horizontal excitation and bidirectional excitation are output through dynamic-static coupling numerical calculation. Based on the dynamic load response characteristics, the direction of the maximum harmful vibration for each support structure is determined, providing a reference for the optimization of support parameters.

2. The method according to claim 1, wherein, The steps for constructing the anchor bolt axial force model include: Determine the size of the model and define the range as 3 times the orifice radius from the blast hole as the simulation range; Set the minimum center frequency and critical damping ratio of the model; Set static dynamic boundaries for the model, and set independent dampers in the normal and tangential directions of the model boundaries.

3. The method according to claim 2, wherein, The steps for constructing the anchor bolt axial force model also include, The maximum mesh size of the control model is within 8m, and the maximum mesh size is less than or equal to 1 / 8 to 1 / 10 of the wavelength corresponding to the highest frequency of the input waveform.

4. The method of claim 1, wherein, The training steps for the anchor bolt axial force model include: Set up monitoring points for the model's dynamic response; First, a gravity field is applied to the model to achieve static stress equilibrium, forming an initial static stress field; Dynamic loads are applied inside the boundaries or nodes of the model. The types of dynamic loads include acceleration time history, velocity time history, pressure time history, or concentrated force time history. The dynamic loads are respectively excited vertically, horizontally, and bidirectionally. Through dynamic-static coupling calculation, the dynamic load response characteristics, including dynamic displacement, peak axial force, and volume of the plastic zone of the structure, are extracted from the monitoring points.

5. The method of claim 1, wherein the method further comprises: The steps of inputting the initial waveform file into a preset anchor bolt axial force model, combining it with the corresponding anchor bolt static load analysis result file, and outputting the dynamic load response characteristics under vertical excitation, horizontal excitation, and bidirectional excitation through dynamic-static coupling numerical calculation include: The static load file command is used to retrieve the static load analysis result file of the monitoring points; The waveform file is read into the measured waveform file collected by the monitoring point using the waveform file read command; The waveform file is loaded onto the static load analysis result file using the waveform file loading command; The dynamic load response results are obtained by solving the problem.

6. The method according to any one of claims 1-5, wherein, The step of filtering the measured waveform file to obtain the initial waveform file includes: Determine the transverse wave velocity of the measured waveform file; Based on the shear wave velocity and the preset filtering frequency, the filtering wavelength corresponding to the filtering frequency is calculated. Calculate the grid size range based on the filtered wavelength; The maximum value of the grid size range is compared with a preset maximum grid size threshold; When the maximum value of the grid size range is greater than the maximum grid size threshold, the filtering frequency is used as the upper limit value; The initial waveform file is obtained by filtering out waveforms with frequencies greater than the upper limit from the measured waveform file.

7. A device for intelligent analysis of changes in the axial force of a rock bolt under dynamic load, characterized by include, The data module is used to acquire the measured waveform files generated during blasting of several anchor bolts deployed on the support structure. The filtering module is used to filter the measured waveform file to obtain an initial waveform file; The dynamic load response module is used to input the initial waveform file into the preset anchor bolt axial force model, combine it with the corresponding anchor bolt static load analysis result file, and output the dynamic load response characteristics under vertical excitation, horizontal excitation and bidirectional excitation through dynamic-static coupling numerical calculation. The analysis module is used to determine the direction of maximum harmful vibration for each support structure based on the dynamic load response characteristics, providing a reference for the optimization of support parameters.

8. A computer device, comprising: The method includes a memory, a processor, and a computer program stored in the memory, wherein the processor executes the computer program to implement the steps of the method according to any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, implements the steps of the method according to any one of claims 1 to 6.

10. A computer program product, characterised in that, It includes a computer program that, when executed by a processor, implements the steps of the method according to any one of claims 1 to 6.

Images