Joint roughness quantification method and system considering damage range and damage degree

Abstract

The present application relates to the technical field of data processing, and provides a joint roughness quantification method and system that takes into account the damage range and damage degree. By digitizing the standard contour line, obtaining the coordinates of each sampling point, calculating the geometric inclination, protrusion height and length of each micro-segment that constitutes the contour line, and then establishing the quantitative relationship between the inclination and the potential damage area, and the height and the potential damage degree, respectively, and proposing two independent morphological parameters based on the quantitative relationship, namely the effective angle and the effective height. Finally, the quantitative relationship between the effective angle, the effective height and the joint roughness is fitted to obtain a joint roughness quantification model. The model fully considers the geometric and mechanical characteristics of the joint, can reflect the joint influence of the two independent morphological parameters of the effective angle and the effective height on the joint roughness, can reasonably explain the relationship between the roughness and the angle, height and other morphologies, and the model structure is simple, which is easier to promote and apply.

Description

Technical Field

[0001] The present application relates to the technical field of digital data processing, and in particular to a joint roughness quantification method and system taking into account damage range and damage degree. Background Art

[0002] Rock mass is composed of rock blocks and cracks. The existence of cracks significantly reduces the strength of the rock mass. Joints and cracks become the weak parts inside the rock mass and are the key factors controlling the mechanical properties of the rock mass. The morphological characteristics of the joint surface have a significant impact on its shear mechanical properties. Domestic and foreign scholars have conducted a lot of research on the quantitative methods of the morphological characteristics of the joint surface and have achieved many research results.

[0003] However, the current technology has the following problems: (1) The morphological characterization parameters are all based on the geometric characteristics of the contour line, and fail to consider the compression and shear failure process of the contour line, fail to reflect the damage range and degree of the contour line, and fail to combine the geometric characteristics with the mechanical characteristics; (2) The morphological parameters usually only consider the inclination of the contour line, and a small number of parameters consider the influence of the contour line height. There is a lack of indicators that consider the combined effect of the two at the same time; (3) The expression of some morphological characteristics is very cumbersome, and the combination of parameters lacks a reasonable explanation.

[0004] Therefore, it is necessary to provide an improved technical solution to address the above-mentioned deficiencies in the prior art. Summary of the invention

[0005] The purpose of the present application is to provide a joint roughness quantification method and system that takes into account the damage range and damage degree, so as to solve or alleviate the problems existing in the above-mentioned prior art.

[0006] In order to achieve the above objectives, this application provides the following technical solutions:

[0007] This application provides a joint roughness quantification method that takes into account the damage range and damage degree, including:

[0008] Step S101, digitizing the standard contour line and obtaining the coordinates of each sampling point on the standard contour line;

[0009] Step S102, adjusting any standard contour line to a horizontal position, and calculating the total length of the shear side of the contour line and the geometric inclination and convex height of each micro segment according to the coordinates of each sampling point on the contour line; wherein a micro segment refers to a segment obtained by connecting any two adjacent sampling points on the contour line;

[0010] Step S103: fitting the relationship between the critical angle and the potential damage area based on the total length of the shear side and the geometric inclination of each micro-segment, and drawing a first relationship curve between the critical angle and the potential damage area;

[0011] Step S104: fitting the relationship between the critical height and the potential damage degree based on the total length of the shear side and the convex height of each micro-segment, and drawing a second relationship curve between the critical height and the potential damage degree;

[0012] Step S105, defining the area of ​​the geometric space enclosed by the first relationship curve and the X-axis and the Y-axis as the effective angle of the contour line; defining the area of ​​the geometric space enclosed by the second relationship curve and the X-axis and the Y-axis as the effective height of the contour line;

[0013] Step S106: Fit the relationship between the effective angle, the effective height and the joint roughness by the least square method to obtain a quantitative model of the joint roughness that takes into account the potential damage range and damage degree; the expression of the model is:

[0014] JRC=0.83·θ′ 1.14 ·h′ 0.31 ,

[0015] Where JRC is the joint roughness, θ′ is the effective angle, and h′ is the effective height.

[0016] Preferably, after step S102 and before step S103, the method further includes:

[0017] According to the geometric inclinations of each micro-segment, the maximum geometric inclination is obtained;

[0018] The value range from 0 to the maximum geometric inclination angle is sampled at a specified angle interval to obtain multiple angle sampling values, and each angle sampling value is used as a critical angle.

[0019] Preferably, step S103 is specifically:

[0020] Based on the total length of the shear side and the geometric inclination of each micro-segment, a database of the relationship between the critical angle and the potential damage area is established;

[0021] The database of the relationship between the critical angle and the potential damage area is statistically analyzed, and the relationship between the critical angle and the potential damage area is fitted according to the result of the statistical analysis to obtain a fitting model of the potential damage area; the expression of the fitting model of the potential damage area is as follows:

[0022]

[0023] In the formula, represents the potential contact length caused by the angle, θ 临界 represents the critical angle, and n represent fitting parameters, exp represents an exponential function with the natural constant e as the base;

[0024] Using the fitting model of the potential damage area, a first relationship curve with critical angle values ​​ranging from 0 to the maximum geometric inclination angle is drawn.

[0025] Preferably, any critical angle is taken as the current critical angle, and the following steps are performed to establish a relationship between the current critical angle and the potential damage area:

[0026] Search all micro-segments whose geometric inclination angles are greater than the current critical angle to obtain a first critical micro-segment set;

[0027] Calculate the sum of the lengths of all micro-segments in the first critical micro-segment set to obtain a first contact length;

[0028] Calculate the ratio of the first contact length to the total length of the micro-segments with inclination angles greater than 0 in all the micro-segments, use the ratio as the quantitative index value of the potential damage area corresponding to the current critical angle, and record the corresponding relationship between the current critical angle and the quantitative index value of the potential damage area in a database;

[0029] Traverse all critical angles and execute the above steps until all critical angles in the database have corresponding quantitative index values ​​of potential damage areas, and end the loop to obtain a database of the relationship between critical angles and potential damage areas.

[0030] Preferably, after step S102 and before step S104, the method further includes:

[0031] According to the convex heights of each micro-segment, the maximum convex height is obtained;

[0032] The value range from 0 to the maximum bulge height is sampled at specified height intervals to obtain multiple height sampling values, and each height sampling value is used as a critical height.

[0033] Preferably, step S104 is specifically:

[0034] Based on the total length of the shear side and the convex height of each micro-segment, a database of the relationship between the critical height and the potential damage degree is established;

[0035] The database of the relationship between the critical height and the potential damage degree is statistically analyzed, and the relationship between the critical height and the potential damage degree is fitted according to the results of the statistical analysis to obtain a fitting model of the potential damage degree. The expression of the fitting model of the potential damage degree is as follows:

[0036]

[0037] In the formula, represents the potential contact length caused by the protrusion height, h临界 represents the critical height, and m represent fitting parameters; exp represents an exponential function with the natural constant e as the base;

[0038] Using the fitting model of the potential damage degree, a second relationship curve is drawn with the critical height value ranging from 0 to the maximum protrusion height.

[0039] Preferably, any critical height is taken as the current critical height, and the following steps are performed to establish a relationship between the current critical height and the potential damage degree:

[0040] Search all micro-segments whose protrusion height is greater than the current critical height to obtain a second critical micro-segment set;

[0041] Calculate the sum of the lengths of all micro-segments in the second critical micro-segment set to obtain a second contact length;

[0042] Calculate the ratio of the second contact length to the total length of the micro-segments whose protrusion height is greater than 0 in all the micro-segments, use the ratio as the quantitative index value of the potential damage degree corresponding to the current critical height, and record the corresponding relationship between the current critical height and the quantitative index value of the potential damage degree in the database;

[0043] Traverse all critical heights and execute the above steps until all critical heights in the database have corresponding quantitative index values ​​of potential damage levels, and end the loop to obtain a database of the relationship between critical heights and potential damage levels.

[0044] Preferably, the angular interval is 0.1 degrees.

[0045] Preferably, the height interval is 0.1 mm.

[0046] The embodiment of the present application also provides a joint roughness quantification system that takes into account damage range and damage degree, including:

[0047] A digitizing unit configured to digitize the standard contour line and obtain the coordinates of each sampling point on the standard contour line;

[0048] The first calculation unit is configured to adjust any standard contour line to a horizontal position, and calculate the total length of the shear side of the contour line and the geometric inclination angle and the convex height of each micro segment according to the coordinates of each sampling point on the contour line; wherein the micro segment refers to a segment obtained by connecting any two adjacent sampling points on the contour line; the total length of the shear side is the sum of the lengths of the micro segments with geometric inclination angles greater than 0 among all the micro segments;

[0049] A first fitting unit is configured to fit the relationship between the critical angle and the potential damage area based on the total length of the micro-segments with inclination angles greater than 0 in all the micro-segments and the geometric inclination angles of each micro-segment, and draw a first relationship curve between the critical angle and the potential damage area;

[0050] a second fitting unit configured to fit the relationship between the critical height and the potential damage degree based on the total length of the micro-segments with convex heights greater than 0 in all the micro-segments and the convex heights of the respective micro-segments, and draw a second relationship curve between the critical height and the potential damage degree;

[0051] The second calculation unit is configured to define the area of ​​the geometric space enclosed by the first relationship curve and the X-axis and the Y-axis as the effective angle of the contour line; and define the area of ​​the geometric space enclosed by the second relationship curve and the X-axis and the Y-axis as the effective height of the contour line;

[0052] The third fitting unit is configured to fit the relationship between the effective angle, the effective height and the joint roughness by the least square method, and obtain a quantitative model of the joint roughness that takes into account the potential damage range and the damage degree; the expression of the model is:

[0053] JRC=0.83·θ′ 1.14 ·h′ 0.31 ,

[0054] Where JRC is the joint roughness, θ′ is the effective angle, and h′ is the effective height.

[0055] The technical solution provided by the embodiment of the present application has the following advantages:

[0056] This technical solution fully considers the influence of the geometric inclination and uplift height of the micro-segment on the damage range and degree during the contour line compression and shear failure process, and establishes the quantitative relationship between the geometric inclination and uplift height and the potential damage area and potential damage degree respectively. On this basis, two independent morphological parameters are proposed, namely the effective height h′ and the effective angle θ′. These two morphological parameters jointly reflect the roughness of the joint surface, laying the foundation for subsequent research to use the joint surface damage area as a bridge to establish the functional relationship between the joint surface morphological characteristics and the shear strength. BRIEF DESCRIPTION OF THE DRAWINGS

[0057] The drawings constituting part of the present application are used to provide a further understanding of the present application. The exemplary embodiments and descriptions of the present application are used to explain the present application and do not constitute an improper limitation on the present application. Among them:

[0058] Figure 1 A schematic flow chart of a joint roughness quantification method taking into account damage range and damage degree according to some embodiments of the present application.

[0059] Figure 2 An example diagram of a relationship curve between a critical angle and a potential damage area provided according to some embodiments of the present application.

[0060] Figure 3 Another example diagram of a relationship curve between a critical angle and a potential damage area provided according to some embodiments of the present application.

[0061] Figure 4 This is another example diagram of a relationship curve between a critical angle and a potential damage area according to some embodiments of the present application.

[0062] Figure 5 An example diagram of a relationship curve between critical height and potential damage degree provided according to some embodiments of the present application.

[0063] Figure 6 Another example diagram of a relationship curve between critical height and potential damage degree provided according to some embodiments of the present application.

[0064] Figure 7 This is another example diagram of a relationship curve between critical height and potential damage degree provided according to some embodiments of the present application.

[0065] Figure 8 This is a schematic diagram of the comparison results between the calculated value and the true value of the joint roughness quantification model JRC provided according to some embodiments of the present application.

[0066] Fig. 9 A schematic diagram of the structure of a joint roughness quantification system taking into account damage range and damage degree according to some embodiments of the present application. DETAILED DESCRIPTION

[0067] The present application will be described in detail below with reference to the accompanying drawings and in conjunction with embodiments. Each example is provided by way of explanation of the present application and does not limit the present application. In fact, it will be clear to those skilled in the art that modifications and variations may be made in the present application without departing from the scope or spirit of the present application. For example, a feature shown or described as a part of an embodiment may be used in another embodiment to produce yet another embodiment. Therefore, it is desired that the present application includes such modifications and variations within the scope of the appended claims and their equivalents.

[0068] In the following description, the terms "first / second / third" involved are merely used to distinguish similar objects and do not represent a specific ordering of the objects. It can be understood that "first / second / third" can be interchanged with a specific order or sequence where permitted, so that the embodiments of the present application described herein can be implemented in an order other than that illustrated or described herein.

[0069] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as those commonly understood by those skilled in the art to which the present disclosure belongs. The terms used herein are only for the purpose of describing the embodiments of the present disclosure and are not intended to limit the present disclosure.

[0070] Embodiment 1:

[0071] The present application embodiment provides a joint roughness quantification method that takes into account the damage range and damage degree, such as Figures 1 to 9 As shown, the method includes:

[0072] Step S101: digitize the standard contour line and obtain the coordinates of each sampling point on the standard contour line.

[0073] It should be noted that the standard contour line described in this embodiment refers to the JRC (Joint Roughness Coefficient) standard contour line, also known as the JRC standard profile line, which is a standardized representation method for describing the surface morphology characteristics of rock joints. The JRC standard was proposed by Barton and Choubey, which is based on the observation of the joint surface contour line and divides the contour of the joint surface into a variety of standardized morphologies, each of which represents the surface morphology of the joint with different degrees of roughness and irregularity. Standardized representation of the joint surface morphology characteristics can quantify the roughness and irregularity of the joint surface, which has important applications in fields such as rock mechanics, geological engineering, and underground engineering.

[0074] In this embodiment, digital processing refers to converting the standard contour line into digital data for more in-depth analysis and quantitative description. There are many ways to implement specific digital processing. For example, a laser scanner, a surface profiler, etc. can be used to measure the standardized joint surface to obtain the morphological data of the joint surface. Exemplarily, this embodiment uses a document scanner to obtain a high-quality picture of the JRC standard section line, and then uses GetData Graph Digitizer, MATLAB and other software to digitize the standard contour line, and then obtain the digital expression result of the JRC standard section line.

[0075] On the basis of obtaining the digital expression result of the JRC standard profile, the digital standard contour is sampled and the coordinates of each sampling point are extracted. It should be noted that there are many sampling methods, such as uniform interval sampling method, feature point sampling method, etc. This embodiment does not limit the specific sampling method.

[0076] Exemplarily, in order to obtain enough sampling points to effectively represent the shape and characteristics of the curve, this embodiment uses a uniformly spaced sampling method with a sampling interval of 0.5 mm to divide the standard contour line into several equal parts, and records the coordinates of each intersection point and stores them in a file as the initial file for calculation.

[0077] Step S102: adjust any standard contour line to a horizontal position, and calculate the total length of the shear side of the contour line and the geometric inclination angle and convex height of each micro-segment according to the coordinates of each sampling point on the contour line.

[0078] Based on step S101, it can be known that after digital processing, any contour line contains a series of endpoints (i.e., sampling points). Among them, a micro segment refers to a segment obtained by connecting any two adjacent sampling points on the contour line. All micro segments are combined together to form a complete contour line.

[0079] In this embodiment, first, for any contour line, the contour line is adjusted to a horizontal position, and the fitting plane of the contour line is Y = 0. Then, according to the coordinates of each sampling point obtained in step S101, the geometric inclination angle, protrusion height and length of each micro-segment in all micro-segments on the contour line are calculated according to formulas (1), (2) and (3):

[0080]

[0081]

[0082]

[0083] In the formula, θ i 、h i , L i are the geometric inclination, protrusion height and length of the i-th micro-line segment, i is a positive integer, (x i ,y i )、(x i+1 ,y i+1 ) are the starting point coordinates and the ending point coordinates of the i-th micro-segment respectively.

[0084] Then, search for all micro-segments with geometric inclination angles greater than 0, sum up the lengths of all micro-segments in the micro-segments, and obtain the total length of the shear side. 迎剪侧 In other words, the total length of the shear side is the sum of the lengths of the micro-segments with geometric inclination angles greater than 0 in all micro-segments.

[0085] Studies have shown that the geometric inclination angle and the height of the protrusions are both key factors affecting the roughness of the joints. The geometric inclination angle mainly affects the damage area of ​​the damaged area, while the height of the protrusion mainly affects the degree of damage. In this embodiment, quantitative relationships between the two and the potential damage area and the degree of damage are established respectively. On this basis, two independent morphological parameters, namely the effective height h′ and the effective angle θ′, are proposed to jointly reflect the roughness of the joint surface. The following is a detailed introduction to the process of establishing the quantitative relationship between the geometric inclination angle, the protrusion height, the potential damage area, and the degree of damage.

[0086] First, the process of establishing the quantitative relationship between the angle and the potential damage area is described, which is specifically implemented through step S103.

[0087] Step S103 is as follows: based on the total length of the shear side and the geometric inclination of each micro-segment, the relationship between the critical angle and the potential damage area is fitted, and a first relationship curve between the critical angle and the potential damage area is drawn.

[0088] In this embodiment, considering that the larger the geometric inclination of the micro-segment is, the easier it is to be sheared during the shearing process, the relationship between the critical angle and the potential damage area is fitted according to the geometric inclination of each micro-segment. Among them, the critical angle refers to the critical value of the angle, exceeding this angle may cause the destruction or instability of the joint, and the potential damage area refers to the area where reaching or exceeding the critical angle may cause damage, destruction or deformation of the rock. It should be understood that for different joint morphological characteristics, different critical angles may cause different damage areas. The geometric inclination of each micro-segment collected on the joint surface can be compared with the critical angle, and the micro-segment can be judged whether it will be damaged based on the comparison result. For any specified rock, its critical angle can be set based on prior knowledge, or it can be obtained through interval sampling.

[0089] Specifically, the critical angle can be determined based on the geometric inclination of each micro segment through an interval sampling method. That is, after the geometric inclination of each micro segment is calculated in step S102 and before step S103, the method further includes: obtaining the maximum geometric inclination according to the geometric inclination of each micro segment; performing interval sampling on the value range from 0 to the maximum geometric inclination at a specified angle interval to obtain multiple angle sampling values, and taking each angle sampling value as a critical angle θ 临界 .

[0090] Preferably, the sampling interval (ie, the angular interval) is 0.1 degrees.

[0091] For example, the maximum geometric inclination angle (θ max Then specify the angle interval as 0.1°, and use θ to represent the angle from 0 to the maximum geometric inclination.max Sampling is performed at intervals of 0.1 degrees to obtain multiple angle sampling values, and the critical angle θ 临界 For example, it can be expressed as: 临界 =0, 0.1, 0.2...θ max .

[0092] Since there is a corresponding relationship between the critical angle and the damaged area, it is possible to determine whether the micro-segment has the possibility of being damaged based on the size relationship between the geometric inclination of each micro-segment and the critical angle. Then, based on the total length of the micro-segments with inclination angles greater than 0 in all micro-segments and the length of the micro-segments that may be damaged, the potential damaged area is quantified, thereby collecting the corresponding relationship between the critical angle and the potential damaged area. Finally, the relationship between the critical angle and the potential damaged area is fitted based on the corresponding relationship.

[0093] It should be noted that when fitting the relationship between the critical angle and the potential damage area, there are many specific fitting methods, such as linear fitting, polynomial fitting, nonlinear fitting, etc. This embodiment does not limit this.

[0094] Exemplarily, the process may specifically include the following steps: establishing a relationship database between critical angles and potential damage areas based on the total length of all micro-segments and the geometric inclination of each micro-segment; performing statistical analysis on the relationship database between critical angles and potential damage areas, and fitting the relationship between critical angles and potential damage areas according to the results of the statistical analysis to obtain a fitting model of the potential damage area; the expression of the fitting model of the potential damage area is as follows:

[0095]

[0096] In the formula, represents the potential contact length caused by the angle, θ 临界 represents the critical angle, and n represent fitting parameters, and exp represents an exponential function with the natural constant e as the base.

[0097] Using the fitting model of the potential damage area, a first relationship curve with critical angle values ​​ranging from 0 to the maximum geometric inclination angle is drawn.

[0098] It should be noted that the potential contact length caused by the angle is used to quantitatively characterize the potential damage area.

[0099] Figure 2 to Figure 4 Some examples of relationship curves (first relationship curves) between the critical angle and the potential contact length caused by the angle are shown, where the X-axis is the critical angle and the Y-axis is the potential contact length caused by the angle. It can be seen from the figure that there is a negative correlation between the potential contact length caused by the angle and the critical angle.

[0100] Furthermore, establishing a relationship database between critical angles and potential damage areas can be broken down into the following steps: taking any critical angle as the current critical angle, and executing the following steps to establish a relationship between the current critical angle and the potential damage area: searching all micro-segments for micro-segments whose geometric inclination angles are greater than the current critical angle to obtain a first critical micro-segment set; calculating the sum of the lengths of all micro-segments in the first critical micro-segment set to obtain a first contact length; calculating the ratio of the first contact length to the total length of the shear side, and using the ratio as a quantitative index value of the potential damage area corresponding to the current critical angle, and recording the corresponding relationship between the current critical angle and the quantitative index value of the potential damage area in the database; traversing all critical angles and executing the above steps until all critical angles in the database have corresponding quantitative index values ​​of the potential damage area, and ending the loop to obtain a relationship database between critical angles and potential damage areas.

[0101] The above steps are illustrated below with examples.

[0102] Continuing from the above, when setting the critical angle θ 临界 =0, 0.1, 0.2...θ max After that, for any critical angle θ 临界 , search for a geometric inclination angle greater than the critical angle θ 临界 All micro-segments of θ i >θ 临界 , get the first critical micro-segment set, calculate the total length of this part of micro-segments, and get the first contact length, recorded as L 接触 , and then use L 接触 With L 迎剪侧 Ratio To quantify the potential contact / damage area, and finally establish its contour Database, namely the relationship database between critical angles and potential damage areas.

[0103] Then, for the contour line The database is statistically analyzed, and based on the results of the statistical analysis, an exponential function is proposed to fit the potential damage area, that is, to establish the θ 临界 and The experiment shows that the model has a good fitting effect. The specific fitting results are as follows: Figure 2 to Figure 4 As shown. When the potential contact length is Equal to 1 / e.

[0104] After completing the fitting of the relationship between the angle and the damage area, the following describes the process of establishing the quantitative relationship between the protrusion height and the potential damage degree, which is mainly implemented by step S104.

[0105] Step S104: fitting the relationship between the critical height and the potential damage degree based on the total length of the shear side and the protrusion height of each micro-segment, and drawing a second relationship curve between the critical height and the potential damage degree.

[0106] Considering that the greater the convex height of the micro-segment, the more obvious the damage degree of the joint surface during the shearing process, in step S104, the relationship between the critical height and the potential damage degree is fitted according to the convex height of each micro-segment.

[0107] Among them, the critical height of the joint surface refers to the minimum height at which the joint surface is damaged or slides when the rock or rock body is subjected to a sufficiently large shear stress under stress. When the height of the joint surface exceeds the critical height, it may cause slippage or damage. Since different critical heights cause different degrees of damage, and the critical height is related to the joint characteristics of the rock, such as joint gaps, friction coefficients between joint surfaces, etc., for any given rock joint surface, its critical height can be set based on prior knowledge, or it can be obtained by interval sampling based on the results of digital processing.

[0108] Specifically, the critical height can be determined based on the convex height of each micro-segment through an interval sampling method. That is, after the convex height of each micro-segment is calculated in step S102 and before step S104, it also includes: obtaining the maximum convex height according to the convex height of each micro-segment; performing interval sampling on the value range from 0 to the maximum convex height according to the specified height interval to obtain multiple height sampling values, and taking each height sampling value as a critical height.

[0109] Preferably, the height interval may be 0.1 mm.

[0110] For example, the maximum convex height (expressed as h) can be found by counting the convex height values ​​of each micro-segment. max Indicates), and then specify the height interval as 0.1mm, from 0 to the maximum convex height h max The sample is taken at intervals of 0.1 mm to obtain multiple height sampling values ​​and set the critical height h 临界 =0, 0.1, 0.2...h max .

[0111] After setting the critical height, the relationship between the protrusion height and the degree of damage can be fitted based on the set critical height. The fitting process can include the following specific steps: based on the total length of the micro-segments with protrusion heights greater than 0 in all micro-segments and the protrusion heights of each micro-segment, a database of the relationship between the critical height and the potential degree of damage is established; statistical analysis is performed on the database of the relationship between the critical height and the potential degree of damage, and the relationship between the critical height and the potential degree of damage is fitted based on the results of the statistical analysis to obtain a fitting model for the potential degree of damage, the expression of the fitting model for the potential degree of damage is as follows:

[0112]

[0113] In the formula, represents the potential contact length caused by the protrusion height, h 临界 represents the critical height, and m represent fitting parameters; exp represents the exponential function with the natural constant e as the base.

[0114] Using the fitting model of the potential damage degree, a second relationship curve is drawn with the critical height value ranging from 0 to the maximum protrusion height.

[0115] It should be noted that the potential contact length caused by the protrusion height is used to quantitatively characterize the potential damage degree.

[0116] Figure 5 to Figure 7 Some examples of the relationship curve between the critical height and the potential contact length caused by the protrusion height (the second relationship curve) are shown, wherein the X-axis is the critical height and the Y-axis is the potential contact length caused by the protrusion height. It can be seen from the figure that there is a negative correlation between the potential contact length caused by the protrusion height and the critical height.

[0117] Furthermore, establishing a relationship database between critical height and potential damage degree can be refined into the following steps: taking any critical height as the current critical height, and executing the following steps to establish a relationship between the current critical height and the potential damage degree: searching all micro-segments for micro-segments whose protrusion height is greater than the current critical height to obtain a second critical micro-segment set; calculating the sum of the lengths of all micro-segments in the second critical micro-segment set to obtain a second contact length; calculating the ratio of the second contact length to the total length of the micro-segments whose protrusion height is greater than 0 in all micro-segments, using the ratio as a quantitative index value of the potential damage degree corresponding to the current critical height, and recording the corresponding relationship between the current critical height and the quantitative index value of the potential damage degree in the database; traversing all critical heights and executing the above steps until all critical heights in the database have corresponding quantitative index values ​​of the potential damage degree, ending the loop, and obtaining a relationship database between critical heights and potential damage degrees.

[0118] The following is an example to illustrate the above steps.

[0119] Continuing from the above, when setting the critical height h 临界 =0, 0.1, 0.2...h max Afterwards, for any critical height h 临界 , search for a bump height exceeding the critical value (i.e. h i >h 临界 ) to obtain the second critical micro-segment set, calculate the total length of this part of micro-segments, and obtain the second contact length, which is recorded as L 损伤 , using L 损伤 With L 迎剪侧 Ratio Quantify the potential damage area. Database, namely, database of relationship between critical height and potential damage degree.

[0120] Then, by The database was statistically analyzed and fitted in the form of exponential function to obtain the fitting model of potential damage degree and establish h 临界 and The experiment shows that the model has a good fitting effect. The specific fitting effect examples are as follows Figure 5 to Figure 7 As shown. When the potential contact length is Equal to 1 / e.

[0121] In addition, since the fitting model of the potential damage area reflecting the relationship between the angle and the damage area has the same structure as the fitting model of the potential damage degree reflecting the relationship between the height and the damage degree, the consistent structural expression makes the model easier to read and understand, which is conducive to making a reasonable explanation of the relationship between the two, and it is easier to understand the overall connection between them and the joint roughness quantification model. In addition, the consistent structure simplifies the maintenance and updating of the model. When the parameters need to be modified or updated, the unified structure can reduce the complexity of operation and adjustment.

[0122] Step S105, defining the area of ​​the geometric space enclosed by the first relationship curve and the X-axis and the Y-axis as the effective angle θ′ of the contour line; defining the area of ​​the geometric space enclosed by the second relationship curve and the X-axis and the Y-axis as the effective height h′ of the contour line.

[0123] The effective angle θ′ is used to measure the effect of the micro-segment geometric inclination on the roughness of the contour line, and the effective height h′ is used to measure the effect of the micro-segment protrusion height on the roughness of the contour line. Solving the area of ​​the geometric space enclosed by the first relationship curve, the second relationship curve, the X-axis, and the Y-axis can be achieved by, for example, an integral method, which is not limited in this embodiment.

[0124] Step S106: Fit the relationship between the effective angle θ′, the effective height h′ and the joint roughness by the least square method to obtain a quantitative model of joint roughness that takes into account the potential damage range and damage degree; the expression of the model is:

[0125] JRC=0.83·θ′ 1.14 ·h′ 0.31 (6)

[0126] Where JRC is the joint roughness, θ′ is the effective angle, and h′ is the effective height.

[0127] In the traditional method, the influence of joint height and inclination on the roughness can be expressed by characteristic angle, characteristic height and two fitting parameters, which requires four indicators in total, which is very cumbersome. The joint roughness quantification model provided in this embodiment only needs two indicators, effective angle and effective height, to achieve a reasonable expression of joint roughness. Moreover, the effective angle and effective height parameters used in the model have clear physical meanings and can jointly express the influence of micro-convex height and angle on joint roughness.

[0128] In addition, this example also verifies the accuracy of the joint roughness quantification model by making a series of mortar joint replicas and conducting three-dimensional laser scanning and direct shear tests on the joint specimens. Figure 8 shown. Figure 8 In the figure, the X-axis is the true value of the joint roughness (true value of JRC), and the Y-axis is the calculated value of JRC obtained by the model calculation. Figure 8 It can be seen that the JRC calculated value obtained by the model provided in this embodiment is basically consistent with the true value, with high accuracy, and can meet the needs of subsequent research on rock shear strength.

[0129] In summary, the present embodiment provides a joint roughness quantification method that takes into account the damage range and damage degree. By digitizing the standard contour line, the coordinates of each sampling point are obtained. On this basis, the geometric inclination angle, protrusion height and length of each micro-segment that constitutes the contour line are calculated. Then, the quantitative relationship between the inclination angle and the potential damage area, and the height and the potential damage degree are established respectively. Based on the quantitative relationship, two independent morphological parameters, namely the effective angle and the effective height, are proposed. The effective angle is used to characterize the influence of the micro-segment inclination angle on the roughness, and the effective height is used to characterize the influence of the micro-segment height on the roughness. Finally, the effective angle and the effective height are combined, and the quantitative relationship between the two and the joint roughness is fitted. The obtained joint roughness quantification model can jointly reflect the influence of the two independent morphological parameters of the effective angle and the effective height on the joint roughness, which can provide a theoretical basis for subsequent research to use the joint surface damage area as a bridge and establish a functional relationship between the joint surface morphological characteristics and the shear strength.

[0130] Embodiment 2:

[0131] The present application embodiment provides a joint roughness quantification system that takes into account the damage range and damage degree, such as Fig. 9 As shown, the system includes: a digitization unit 501, a first calculation unit 502, a first fitting unit 503, a second fitting unit 504, a second calculation unit 505 and a third fitting unit 506. Specifically:

[0132] The digitizing unit 501 is configured to digitize the standard contour line and obtain the coordinates of each sampling point on the standard contour line.

[0133] The first calculation unit 502 is configured to adjust any standard contour line to a horizontal position, and calculate the total length of all micro-segments on the contour line and the geometric inclination angle and convex height of each micro-segment according to the coordinates of each sampling point on the contour line.

[0134] Among them, a micro-segment refers to a segment obtained by connecting any two adjacent sampling points on the contour line.

[0135] The first fitting unit 503 is configured to fit the relationship between the critical angle and the potential damage area based on the total length of the micro-segments with inclination angles greater than 0 in all the micro-segments and the geometric inclination angles of each micro-segment, and draw a first relationship curve between the critical angle and the potential damage area.

[0136] The second fitting unit 504 is configured to fit the relationship between the critical height and the potential damage degree based on the total length of the micro-segments with protrusion height greater than 0 in all micro-segments and the protrusion height of each micro-segment, and draw a second relationship curve between the critical height and the potential damage degree.

[0137] The second calculation unit 505 is configured to define the area of ​​the geometric space enclosed by the first relationship curve and the X-axis and the Y-axis as the effective angle of the contour line; and define the area of ​​the geometric space enclosed by the second relationship curve and the X-axis and the Y-axis as the effective height of the contour line.

[0138] The third fitting unit 506 is configured to fit the relationship between the effective angle, the effective height and the joint roughness by the least square method to obtain a joint roughness quantification model that takes into account the potential damage range and damage degree; the expression of the model is:

[0139] JRC=0.83·θ′ 1.14 ·h′ 0.31 ,

[0140] Where JRC is the joint roughness, θ′ is the effective angle, and h′ is the effective height.

[0141] The joint roughness quantification system considering the damage range and damage degree provided in the embodiment of the present application can implement the steps and processes of the joint roughness quantification method considering the damage range and damage degree provided in any of the above embodiments, and achieve the same technical effects, which will not be described one by one here.

[0142] The above description is only a preferred embodiment of the present application and is not intended to limit the present application. For those skilled in the art, the present application may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims

1. A joint roughness quantification method considering damage range and damage degree, characterized in that: include: Step S101, digitizing the standard contour line and obtaining the coordinates of each sampling point on the standard contour line; Step S102, adjusting any standard contour line to a horizontal position, and calculating the total length of the shear side of the contour line and the geometric inclination and convex height of each micro segment according to the coordinates of each sampling point on the contour line; wherein a micro segment refers to a segment obtained by connecting any two adjacent sampling points on the contour line; Step S103: fitting the relationship between the critical angle and the potential damage area based on the total length of the shear side and the geometric inclination of each micro-segment, and drawing a first relationship curve between the critical angle and the potential damage area; Step S104: fitting the relationship between the critical height and the potential damage degree based on the total length of the shear side and the convex height of each micro-segment, and drawing a second relationship curve between the critical height and the potential damage degree; Step S105, defining the area of ​​the geometric space enclosed by the first relationship curve and the X-axis and the Y-axis as the effective angle of the contour line; defining the area of ​​the geometric space enclosed by the second relationship curve and the X-axis and the Y-axis as the effective height of the contour line; Step S106: Fit the relationship between the effective angle, the effective height and the joint roughness by the least square method to obtain a quantitative model of the joint roughness that takes into account the potential damage range and damage degree; the expression of the model is: JRC=0.83·θ′ 1.14 ·h′ 0.31 , Where JRC is the joint roughness, θ′ is the effective angle, and h′ is the effective height.

2. The method according to claim 1, characterized in that After step S102 and before step S103, the method further includes: According to the geometric inclinations of each micro-segment, the maximum geometric inclination is obtained; The value range from 0 to the maximum geometric inclination angle is sampled at a specified angle interval to obtain multiple angle sampling values, and each angle sampling value is used as a critical angle.

3. The method according to claim 2, characterized in that Step S103 is specifically as follows: Based on the total length of the shear side and the geometric inclination of each micro-segment, a database of the relationship between the critical angle and the potential damage area is established; The database of the relationship between the critical angle and the potential damage area is statistically analyzed, and the relationship between the critical angle and the potential damage area is fitted according to the result of the statistical analysis to obtain a fitting model of the potential damage area; the expression of the fitting model of the potential damage area is as follows: In the formula, represents the potential contact length caused by the angle, θ 临界 represents the critical angle, and n represent fitting parameters, exp represents an exponential function with the natural constant e as the base; Using the fitting model of the potential damage area, a first relationship curve with critical angle values ​​ranging from 0 to the maximum geometric inclination angle is drawn.

4. The method according to claim 3, characterized in that Take any critical angle as the current critical angle and perform the following steps to establish the relationship between the current critical angle and the potential damage area: Search all micro-segments whose geometric inclination angles are greater than the current critical angle to obtain a first critical micro-segment set; Calculate the sum of the lengths of all micro-segments in the first critical micro-segment set to obtain a first contact length; Calculate the ratio of the first contact length to the total length of the shear side, use the ratio as the quantitative index value of the potential damage area corresponding to the current critical angle, and record the corresponding relationship between the current critical angle and the quantitative index value of the potential damage area in the database; Traverse all critical angles and execute the above steps until all critical angles in the database have corresponding quantitative index values ​​of potential damage areas, and end the loop to obtain a database of the relationship between critical angles and potential damage areas.

5. The method according to claim 1, characterized in that After step S102 and before step S104, the method further includes: According to the convex heights of each micro-segment, the maximum convex height is obtained; The value range from 0 to the maximum bulge height is sampled at specified height intervals to obtain multiple height sampling values, and each height sampling value is used as a critical height.

6. The method according to claim 5, characterized in that Step S104 is specifically as follows: Based on the total length of the shear side and the convex height of each micro-segment, a database of the relationship between the critical height and the potential damage degree is established; The database of the relationship between the critical height and the potential damage degree is statistically analyzed, and the relationship between the critical height and the potential damage degree is fitted according to the results of the statistical analysis to obtain a fitting model of the potential damage degree. The expression of the fitting model of the potential damage degree is as follows: In the formula, represents the potential contact length caused by the protrusion height, h 临界 represents the critical height, and m represent fitting parameters; exp represents an exponential function with the natural constant e as the base; Using the fitting model of the potential damage degree, a second relationship curve is drawn with the critical height value ranging from 0 to the maximum protrusion height.

7. The method according to claim 6, characterized in that Take any critical height as the current critical height and perform the following steps to establish the relationship between the current critical height and the potential damage level: Search all micro-segments whose protrusion height is greater than the current critical height to obtain a second critical micro-segment set; Calculate the sum of the lengths of all micro-segments in the second critical micro-segment set to obtain a second contact length; Calculate the ratio of the second contact length to the total length of the shear side, use the ratio as the quantitative index value of the potential damage degree corresponding to the current critical height, and record the corresponding relationship between the current critical height and the quantitative index value of the potential damage degree in the database; Traverse all critical heights and execute the above steps until all critical heights in the database have corresponding quantitative index values ​​of potential damage levels, and end the loop to obtain a database of the relationship between critical heights and potential damage levels.

8. The method according to claim 2, characterized in that: The angular interval is 0.1 degree.

9. The method according to claim 5, characterized in that The height interval is 0.1 mm.

10. A joint roughness quantification system considering damage range and damage degree, characterized in that: include: A digitizing unit configured to digitize the standard contour line and obtain the coordinates of each sampling point on the standard contour line; The first calculation unit is configured to adjust any standard contour line to a horizontal position, and calculate the total length of the shear side of the contour line and the geometric inclination angle and convex height of each micro segment according to the coordinates of each sampling point on the contour line; wherein the micro segment refers to a segment obtained by connecting any two adjacent sampling points on the contour line; A first fitting unit is configured to fit the relationship between the critical angle and the potential damage area based on the total length of the shear side and the geometric inclination of each micro-segment, and draw a first relationship curve between the critical angle and the potential damage area; a second fitting unit configured to fit the relationship between the critical height and the potential damage degree based on the total length of the shear side and the convex height of each micro-segment, and draw a second relationship curve between the critical height and the potential damage degree; The second calculation unit is configured to define the area of ​​the geometric space enclosed by the first relationship curve and the X-axis and the Y-axis as the effective angle of the contour line; and define the area of ​​the geometric space enclosed by the second relationship curve and the X-axis and the Y-axis as the effective height of the contour line; The third fitting unit is configured to fit the relationship between the effective angle, the effective height and the joint roughness by the least square method, and obtain a quantitative model of the joint roughness that takes into account the potential damage range and the damage degree; the expression of the model is: JRC=0.83·θ′ 1.14 ·h′ 0.31 , Where JRC is the joint roughness, θ′ is the effective angle, and h′ is the effective height.

Images