Automatic processing method and system for stress-strain data of roots under tension

Abstract

The present invention relates to a method and system for automatically processing stress-strain data of roots under tension. The method uses a computer device to run program code to automatically read stress-strain data of vegetation roots under tension. After automatic calculation, the virtual origin, yield point (yield strength), fracture point (fracture strength), elastic modulus, ramp slope, and yield ratio of the stress-strain data of samples can be obtained in large quantities. The yield point and fracture point positions in the root stress-strain data can be automatically obtained, solving the problem of tedious manual drawing and processing of data and the inability to complete batch processing in existing research, and avoiding the influence of human factors on data calculation.

Description

Technical Field

[0001] The present invention relates to the field of environmental ecological technology, and more particularly to a method for automatically processing stress-strain data of a root system under tension. Background Art

[0002] In the process of plant reinforcement and soil stabilization, such as improving slope stability and strengthening the root-soil bond, the root system plays a primary role. Roots are the plant's nutritional organs that directly contact the external soil, forming a root network with soil particles, thereby firmly anchoring the plant in the soil. The tensile strength of plant roots is one of the most important factors affecting the stability of the root-soil bond and is directly related to the size of the rhizome. The root's ductility before tensioning, and its elastic modulus and yield strength after tensioning also play a crucial role in soil reinforcement.

[0003] The virtual origin, yield point (yield strength), fracture point (fracture strength), elastic modulus, ramp slope, and yield ratio of the root system obtained from the stress-strain data curve are important indicators for evaluating the tensile characteristics of the root system. They are also evaluation indicators for selecting plants for ecological vegetation protection against landslides. However, the relevant parameters are currently obtained by graphical methods. The data obtained by the graphical method not only has large errors, but also requires a high level of data processing personnel and is inefficient. It is only suitable for processing a small amount of data. The graphical method is obviously not capable of processing a wide variety of plant root data. Therefore, the graphical method for processing plant root data is not convenient for large-scale promotion. A more reasonable method for automatically processing stress-strain data under the tensile state of the root system is needed. Summary of the Invention

[0004] The technical problem to be solved by the present invention is to provide a method for automatically processing stress-strain data of a root system under tension, and also provide a system for automatically processing stress-strain data of a root system under tension, in response to the above-mentioned defects of the prior art.

[0005] The technical solution adopted by the present invention to solve its technical problem is:

[0006] A method for automatically processing stress-strain data of a root system under tension is constructed, which includes the following steps:

[0007] Processing the stress-strain data file of the root system under tension state by a file processing computer program to generate a data information file;

[0008] The stress-strain data files are read in sequence by a data processing computer program, and the following operations are performed on each stress-strain data file:

[0009] Determine the strain data as the abscissa and the stress data as the ordinate, obtain each data point of the stress-strain data file, and obtain a stress-strain curve by connecting a plurality of the data points;

[0010] Looping through each data point of the stress-strain data file to obtain the maximum stress value, where the data point corresponding to the maximum stress value is the fracture point;

[0011] Determine a horizontal line equation based on the first data point of each data point in the stress-strain data file, and the straight line determined by the horizontal line equation serves as a virtual horizontal axis;

[0012] The line connecting the coordinates of the breakpoint and the coordinates of the first data point is used as a reference line; the general formula of the straight line equation of the reference line is used to establish a function using the linear programming method, and all data points are substituted into the function at once to obtain the corresponding function value; when the function value is greater than 0, the data point is located on the lower right side of the reference line, and when the function value is less than 0, the data point is located on the upper left side of the reference line; the distance from each data point to the reference line is calculated in sequence using the point-to-straight line distance formula; the data point with a function value less than 0 and the largest straight line distance from the point to the reference line is recorded as point n; the data point with a function value greater than 0 and the largest straight line distance from the point to the reference line between the first data point and point n is recorded as point m; if the If there is no data point with a function value greater than 0 between a data point and point n, the first data point is recorded as point m; point m is the first turning point; a portion of data points between points m and n is intercepted and recorded as array segment R1; array segment R1 is the steep increase segment data; the slope of the steep increase segment of the stress-strain curve is obtained using the stress-strain values ​​of the first data point p and the last data point q of array segment R1. The slope of the steep increase segment of the stress-strain curve represents the elastic modulus of the root system; the regression line equation of array segment R1 is calculated using the least squares method using the first data point p and the last data point q of array segment R1. The regression line corresponding to the regression line equation is the elastic line;

[0013] Finding an intersection point between the elastic line and the virtual horizontal axis, the intersection point is the virtual origin, and the plumb line passing through the virtual origin is the virtual vertical axis;

[0014] Obtaining a dividing point between a steep increase and a slow increase of the stress-strain curve, wherein the dividing point is a yield point;

[0015] Obtain the slow-increasing line of the stress-strain curve. The slope of the slow-increasing line is the plastic modulus, and the ratio of the plastic modulus to the elastic modulus is the yield ratio.

[0016] A graph reflecting the changes in stress-strain data of the root system under tension is obtained.

[0017] The method for automatically processing stress-strain data of a root system under tension according to the present invention, wherein the processing of the stress-strain data file of the root system under tension using a file processing computer program to generate a data information file comprises:

[0018] A working folder is established at a set location in the computer, a data folder is established in the working folder, stress-strain data files of all roots under tension are stored in the data folder, a file processing computer program is run, all stress-strain data files in the data folder are read, the total number of data files, data file name and the total number of data rows of the corresponding data files are recorded, and the records are written into the data information file in the working folder; when writing the records into the data information file in the working folder, the first line of data in the data information file is the total number of data files, the subsequent two lines are the first data file name and the total number of data rows of the first data file, and the subsequent two lines are the second data file name and the total number of data rows of the second data file, until all data file name and the total number of data rows of the corresponding data files are written.

[0019] The method for automatically processing stress-strain data of a root system under tension according to the present invention comprises the following steps: a data processing computer program is used to sequentially read stress-strain data files in a set order; and

[0020] The data processing computer program reads the first line of the data information file to obtain the total number of data files in the data folder, reads the second and third lines of the data information file to obtain the file name of the first data file and the total number of data lines in the data file;

[0021] According to the file name of the obtained data file and the total number of data rows of the data file, the data program reads the content of the corresponding file in the data folder and stores the stress-strain data sequence in the computer memory in the form of an array.

[0022] The method for automatically processing stress-strain data of a root system under tension according to the present invention comprises the following steps: sequentially reading stress-strain data files according to a set sequence using a data processing computer program; and determining whether the data information file is empty. This step adopts:

[0023] Run the data processing computer program to determine whether the data information file is empty. If it is empty, exit the program.

[0024] The method for automatically processing stress-strain data of a root system under tension according to the present invention comprises the following steps: a data processing computer program is used to sequentially read stress-strain data files in a set order; and

[0025] The data processing computer program automatically determines whether the reading of the data information file is completed. If the reading is not completed, the data program continues to automatically read the next two lines of the data information file, obtains the file name of a new data file and the total number of data lines of the data file, and automatically performs a new round of calculations; if the reading is completed, the program exits.

[0026] The method for automatically processing stress-strain data of a root system under tension according to the present invention comprises the following steps: obtaining the dividing point between a steep increase and a slow increase of the stress-strain curve; and

[0027] The data point where the function value is less than 0 and the straight-line distance from the point to the reference line is the maximum is recorded as point a, and point a is the yield point.

[0028] The method for automatically processing stress-strain data of a root system under tension according to the present invention, wherein the method for obtaining the slow-increasing line is as follows:

[0029] With point a as the center and the distances from point p and point q to point a as the radii, determine two concentric circles. All data points between point a and point b and located between these two concentric circles are recorded as array segment R2. The first and last data of array segment R2 are recorded as points u and v, respectively. The least squares method is used to calculate the regression line equation of array segment R2. The straight line corresponding to the regression line equation is a slowly increasing line.

[0030] An automatic processing system for stress-strain data of a root system under tension, comprising a data acquisition module, a data processing module and a data output module;

[0031] The data acquisition module is used to acquire stress-strain data files;

[0032] The data processing module uses a file processing computer program to process the stress-strain data file of the root system under tension to generate a data information file; uses the data processing computer program to read the stress-strain data file in a set order, and performs the following operations on each stress-strain data file: determines the strain data as the horizontal coordinate and the stress data as the vertical coordinate, obtains each data point of the stress-strain data file, and obtains the stress-strain curve by connecting multiple data points; loops through each data point of the stress-strain data file to obtain the maximum stress, and the data point corresponding to the maximum stress is the fracture point; determines a horizontal line equation with the first data point among each data point of the stress-strain data file, and the straight line determined by the horizontal line equation is the virtual horizontal axis; uses the coordinates of the fracture point and the line connecting the coordinates of the first data point as a reference line; uses the general formula of the straight line equation of the reference line to establish a function with a linear programming method, and substitutes all the data points into the function at once to obtain the corresponding function value; when the function value is greater than 0, the data point is located at the lower right side of the reference line, and when the function value is less than 0, the data point is located at the lower right side of the reference line. The upper left side of the test line; use the point-to-straight-line distance formula to calculate the distance of each data point to the reference line in turn; record the data point with a function value less than 0 and the largest straight-line distance from the point to the reference line as point n; record the data point between the first data point and point n with a function value greater than 0 and the largest straight-line distance from the point to the reference line as point m; if there is no data point with a function value greater than 0 between the first data point and point n, then record the first data point as point m; point m is the first turning point; intercept a part of the data points between point m and point n, record it as array segment R1; array segment R1 is the steep increase segment data; using the first data point p and the last data point q of the array segment R1, the least squares method is used to calculate the regression line equation of the array segment R1, and the regression line corresponding to the regression line equation is the elastic line; the intersection of the elastic line and the virtual horizontal axis is obtained, the obtained intersection is the virtual origin, and the plumb line passing through the virtual origin is the virtual vertical axis; the dividing point between the steep increase and the slow increase of the stress-strain curve is obtained, and the dividing point is the yield point; the slow increase line of the stress-strain curve is obtained, the slope of the slow increase line is the plastic modulus, and the ratio of the plastic modulus to the elastic modulus is the yield ratio;

[0033] The data output module is used to output the obtained graph reflecting the stress-strain data changes under the tension state of the root system.

[0034] The beneficial effects of the present invention are as follows: the present invention uses computer programs to process data, which can realize large-scale batch processing and greatly reduce the time cost of data processing. Furthermore, the present invention proposes data processing rules based on mathematical principles, which effectively reduces the calculation errors caused by human factors in data processing; the data processing method proposed by the present invention has a simple calculation principle, is easy to operate on a computer, and requires less computing power of the computer, which shows that the data automation processing method provided by the present invention has strong practicality and is easy to popularize and apply; the data processing method proposed by the present invention is not limited to plant root research, but is also applicable to a large number of other physical stress-strain parameter acquisition studies under tensile conditions, further broadening the scope of application of the method. BRIEF DESCRIPTION OF THE DRAWINGS

[0035] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the present invention will be further described below with reference to the accompanying drawings and embodiments. The drawings described below are only some embodiments of the present invention. Those skilled in the art can also derive other drawings based on these drawings without inventive work.

[0036] Figure 1 This is a flow chart of a method for automatically processing stress-strain data of a root system under tension according to a preferred embodiment of the present invention;

[0037] Figure 2 This is a diagrammatic representation of the changes in stress-strain data of a root system under tension in a preferred embodiment of the present invention;

[0038] Figure 3 This is a file processing diagram of a method for automatically processing stress-strain data of a root system under tension according to a preferred embodiment of the present invention;

[0039] Figure 4 This is a principle block diagram of the automatic processing system for stress-strain data of roots under tension in a preferred embodiment of the present invention. DETAILED DESCRIPTION

[0040] To make the purpose, technical solutions, and advantages of the embodiments of the present invention more clear, the following will be a clear and complete description of the technical solutions in the embodiments of the present invention. Obviously, the embodiments described are some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work shall fall within the scope of protection of the present invention.

[0041] The preferred embodiment of the present invention is a method for automatically processing stress-strain data of a root system under tension, such as Figure 1 See also Figure 2 and Figure 3 , including the following steps:

[0042] Processing the stress-strain data file of the root system under tension state by a file processing computer program to generate a data information file;

[0043] The stress-strain data files are read in sequence by a data processing computer program, and the following operations are performed on each stress-strain data file:

[0044] S01: Determine strain data as the horizontal coordinate and stress data as the vertical coordinate, obtain each data point of the stress-strain data file, and obtain a stress-strain curve by connecting multiple data points;

[0045] S02: looping through each data point of the stress-strain data file to obtain the maximum stress value, where the data point corresponding to the maximum stress value is the fracture point;

[0046] S03: determining a horizontal line equation based on the first data point among the data points in the stress-strain data file, wherein the straight line determined by the horizontal line equation serves as a virtual horizontal axis;

[0047] S04: Obtaining the slope of the steep increase section of the stress-strain curve by a graphical method, wherein the slope of the steep increase section of the stress-strain curve represents the elastic modulus of the root system, and the straight line obtained by the graphical method is the elastic line;

[0048] S05: finding an intersection point between the elastic line and the virtual horizontal axis, the intersection point being the virtual origin, and the plumb line passing through the virtual origin being the virtual vertical axis;

[0049] S06: Obtaining a dividing point between a steep increase and a slow increase of the stress-strain curve, where the dividing point is a yield point;

[0050] S07: Obtain the slow-increasing line of the stress-strain curve. The slope of the slow-increasing line is the plastic modulus, and the ratio of the plastic modulus to the elastic modulus is the yield ratio.

[0051] S08: Obtain a graph reflecting the stress-strain data changes of the root system under tension;

[0052] Optionally, in order to automatically process the stress-strain data, this solution provides a possible implementation method for generating a data information file. Specifically, the method includes the following steps:

[0053] Step a: Create a working folder at a set location on the computer, create a data folder in the working folder, and store all stress-strain data files of the root system under tension in the data folder;

[0054] Step b, running a file processing computer program to read all stress-strain data files in the data folder;

[0055] Step c: record the total number of data files, the data file names, and the total number of data rows in the corresponding data files, and write the records into the data information file in the working folder.

[0056] Among them, for step c, the format of writing to the working folder is: the first line of data in the data information file is the total number of data files, the following two lines are the name of the first data file and the total number of data lines of the first data file, and the following two lines are the name of the second data file and the total number of data lines of the second data file, until all data file names and the total number of data lines of the corresponding data files are written.

[0057] Optionally, with respect to the above-mentioned “reading stress-strain data files sequentially in a set order by a data processing computer program”, an embodiment of the present invention provides a possible implementation method:

[0058] The data processing computer program reads the first line of the data information file to obtain the total number of data files in the data folder; reads the second and third lines of the data information file to obtain the file name of the first data file and the total number of data lines in the data file;

[0059] According to the file name of the obtained data file and the total number of data rows of the data file, the data program reads the content of the corresponding file in the data folder and stores the stress-strain data sequence in the computer memory in the form of an array.

[0060] Optionally, a data processing computer program is run to determine whether the data information file is empty, and if so, the program is exited.

[0061] Regarding “reading stress-strain data files sequentially in a set order by a data processing computer program”, it can be processed based on each reading cycle, and whether to trigger a new reading cycle can be implemented based on the following possible implementation methods:

[0062] The data processing computer program determines whether the reading of the data information file is completed. If the reading is not completed, the data program continues to read the next two lines of the data information file, obtains the file name of a new data file and the total number of data lines of the data file, and performs a new round of calculations;

[0063] If the reading is completed, exit the program;

[0064] The data reading process of the data processing computer program provided by the above strength is explained below with reference to specific examples:

[0065] Combine Figure 2 and Figure 3 , the automated processing method of stress-strain data of roots under tension is described in detail as follows:

[0066] right Figure 2 and Figure 3 Description of content:

[0067] Figure 2 The horizontal axis is strain, which is represented by x in the equation, and the vertical axis is stress, which is represented by y in the equation; point O′ is the virtual origin; point a is the yield point (yield strength), which is also the second turning point; point b is the fracture point (fracture strength); point m is the first turning point; points m, point p, point q, point n, point a, point u, point v and point b are on a curve, which is the stress-strain data curve of the root system; point q is on circle C1 with point a as the center; point p is on circle C2 with point a as the center; point u is the point between point a and point b and closest to circle C1 between circle C1 and circle C2; ​​point v is the point between point a and point b and closest to circle C1 between circle C1 and circle C2 C2 is the nearest point; point m and point n are the points on both sides of line L0 that are farthest away from line L0; line L0 is the line determined by the stress-strain starting point of the root system and point b; line L1 is the regression line obtained by the least squares method in mathematics for the data points between point p and point q (including point p and point q); line L2 is the horizontal line passing through point O′, that is, the imaginary horizontal axis; line L3 is the lower boundary line passing through point O′ and point b; line L4 is the upper boundary line passing through point a and parallel to line L3; line L5 is the plumb line passing through point O′, that is, the imaginary vertical axis; line L6 is the regression line obtained by the least squares method in mathematics for the data points between point u and point v (including point u and point v);

[0068] Figure 3 It means that the working folder project folder contains the data folder data folder, the file processing computer program calc.sh and the data program calc.exe; the data folder data folder contains root stress-strain data files such as s1.txt and s2.txt; running calc.sh can generate the data information file info.txt; running calc.exe can generate the result data file out.txt.

[0069] Reference Figure 2 and Figure 3 The automated processing method for stress-strain data of roots under tension includes the following steps:

[0070] Step 1. Create a working folder named project in any location on the computer that does not endanger the safety of computer operation; create a data folder named data in the project folder; copy all stress-strain data files under the tension state of the root system (such as s1.txt, s2.txt, etc.) to the data folder; run the file processing computer program calc.sh (written in a computer shell script language), automatically read all data files in the data folder data, record the total number of data files, data file name and the total number of data lines of the corresponding data file, and write this record to the data information file info.txt in the project folder; the first line of info.txt saves the total number of data files in the data folder, and then two lines are grouped to store the data file name and the total number of data lines of the corresponding data file.

[0071] Step 2. Run the data program calc.exe. First, it automatically determines whether the info.txt file is empty. If it is empty, it exits the program. If it is not empty, it automatically reads the first three lines of info.txt to obtain the total number of data files in the data folder, the file name of the first data file, and the total number of data lines in the data file.

[0072] Step 2-1: Based on the file name of the obtained data file and the total number of data rows in the data file, the data program calc.exe reads the content of the corresponding file in the data folder and stores the stress-strain data sequence in the computer memory in the form of an array.

[0073] For the above-mentioned method of processing files, see Figure 3 It should be noted that this method is only for the purpose of more conveniently illustrating the specific implementation of the present invention, which is conducive to demonstrating the principles and methods of processing data, and does not make any restrictive settings on the corresponding file or folder names.

[0074] The data program calc.exe automatically loops through the array to find the maximum stress value and records the corresponding strain data at the same time. This stress-strain data is the coordinate of point b, which is the fracture point (fracture strength). The stress-strain data is recorded in the form of point coordinates (the coordinates are composed of two data, the first data is the strain at the point, and the second data is the stress at the point. The following coordinates are marked in this way). The coordinate of point b can be recorded as (x b ,y b ), the coordinates of the first data point are denoted as (x0, y0). Point b and the first data point can form a straight line, which is the reference line.

[0075] y=k0x+b0 (4)

[0076] Where: y is stress, x is strain, k0 and b0 are unknown coefficients.

[0077]

[0078] b0=y b -k0x b (6)

[0079] Step 2-3: The data program calc.exe automatically uses the general formula of the reference line equation to establish a function using the mathematical linear programming method, as shown in formula (7).

[0080] f(x,y)=k0x-y+b0 (7)

[0081] The data program calc.exe automatically brings the data points in the array into function (7) for calculation. When f(x,y)>0, the data point is located on the lower right side of the straight line equation (4). When f(x,y)<0, the data point is located on the upper left side of the straight line equation (4). At the same time, the distance d from the point to the straight line equation (4) is calculated using the point-to-straight line distance formula (8). The coordinates of any point in the array are marked as (x s ,y s ), then the distance formula (8) from a point to a straight line is

[0082]

[0083] Step 2-4: The data program calc.exe automatically finds the data point where f(x,y)<0 and the distance d is the maximum, and records it as point n(x n ,y n ); Find the data point between the first data point in the array and point n when f(x,y)>0 and the distance d is the largest, and record it as point m(x m ,y m ); if there is no data point with f(x,y)>0 between the first data point in the array and point n, then the value of the first data point in the array is recorded as point m, and the coordinates of point m are (x0,y0); point m is the first turning point.

[0084] Step 2-5, the data program calc.exe automatically calculates the number of data points between point m and point n (including point m and point n), finds the data of the middle section between point m and point n (generally the data of the middle 1 / 3 section is taken, but it can also be selected as needed). This part of the data is also a series of data points, recorded as array segment R1; array segment R1 is the steep increase segment data; the first and last data of array segment R1 are recorded as point p(x p ,y p ) and point q(x q ,y q), the least squares method in mathematics is used to calculate the regression line equation of the array segment R1. This equation is the equation of the straight line L1, that is, the elastic line, as shown in formula (9).

[0085] y=k1x+b1 (9)

[0086] In the formula: k1 and b1 are unknown coefficients; for all data between point p and point q (including point p and point q), k1 and b1 can be determined by using the least squares method in mathematics to find the regression line method; k1 is the elastic modulus of the root system after tensioning. The larger its value, the less likely the root system is to deform after being tensioned.

[0087] Step 2-6: The data program calc.exe automatically determines the equation of a horizontal line through the first data point in the array. This equation is the equation of the straight line L2, that is, the imaginary horizontal axis, as shown in formula (10).

[0088] y=y0 (10)

[0089] Where y0 is a constant, which is the y coordinate of the first data point in the array, that is, the stress value of the first data point in the array.

[0090] Step 2-7: The data program calc.exe automatically solves equations (9) and (10) to obtain the coordinates of point O′.

[0091]

[0092] Point O′ is the virtual origin. The magnitude of the abscissa of point O′ represents the ductility of the root system before tensioning; a larger value indicates greater ductility. Use the coordinates of point O′ and point b to determine the equation of the line L3, the lower boundary, as shown in Equation (11).

[0093] y=k3x+b3 (11)

[0094] In the formula: k3, b3 are unknown coefficients.

[0095]

[0096] b3=y b -k3x b (13)

[0097] Step 2-8, the data program calc.exe automatically puts the data points in the array into function (7) in turn, and finds the data point where f(x,y)<0 and the distance d is the maximum, which is recorded as point a(x a ,y a ), this point is the second turning point of the stress-strain data curve, that is, the yield point (yield strength). The equation of the straight line passing through point a and parallel to the straight line L3 is L4, that is, the upper limit, as shown in formula (14).

[0098] y=k4x+b4 (14)

[0099] In the formula: k4, b4 are unknown coefficients.

[0100] k4=k3 (15)

[0101] b4=y a -k4x a (16)

[0102] Step 2-9: The data program calc.exe automatically obtains the plumb line passing through the imaginary origin O′ as the straight line L5, that is, the imaginary vertical axis, as shown in formula (17).

[0103]

[0104] Step 2-10: The data program calc.exe automatically determines two concentric circles (circle C1 and circle C2) with point a as the center and the distances from point p and point q to point a as the radius. All data points between point a and point b and between these two concentric circles are recorded as array segment R2. The first and last data points of array segment R2 are recorded as point u(x u ,y u ) and point v(x v ,y v ), the least squares method in mathematics is used to calculate the regression line equation of the array segment R2. This equation is the equation of the straight line L6, that is, the slowly increasing line, as shown in formula (18).

[0105] y=k6x+b6 (18)

[0106] In the formula: k6 and b6 are unknown coefficients; k6 and b6 can be determined by using the least squares method to find the regression line method for all data between point u and point v, and all data include point u and point v; k6 is the slowly increasing slope.

[0107] Step 2-11: The data program calc.exe automatically determines the yield ratio z by the ratio between the ramp slope k6 and the elastic modulus k1, as shown in formula (19).

[0108]

[0109] Step 2-12. The data program calc.exe automatically outputs the data file name, virtual origin coordinates, yield point coordinates, fracture point coordinates, yield ratio z, reference line slope k0, reference line intercept b0, elastic line slope k1, elastic line intercept b1, lower boundary slope k3, lower boundary intercept b3, upper boundary slope k4, upper boundary intercept b4, ramp line slope k6, and ramp line intercept b6 to the result data file out.txt.

[0110] In step 2-13, the data program calc.exe automatically determines whether reading info.txt has been completed. If so, it proceeds to step 3. If not, the data program calc.exe automatically reads the next two lines of info.txt, obtains the file name and total number of lines of data in the new data file, and automatically jumps to step 2-1, executing steps 2-1 through 2-13 again.

[0111] Step 3: Analyze the data in the result data file out.txt and select the appropriate root system.

[0112] It is worth noting that the naming contents such as file names and program names in the specific implementation methods of the present invention are only for the purpose of more conveniently explaining the specific implementation methods of the present invention, which is conducive to demonstrating the principles and methods of processing data, and do not set any limitative settings for the names; the file processing computer program calc.sh and the data processing computer program cals.exe described in the present invention can be compiled in any programming language that can meet the contents of this specification.

[0113] 1) The present invention introduces a calculation method for reference lines, thereby proposing a mathematical solution method for elastic lines, providing a theoretical basis for computer programming to solve elastic lines. The advantages of this method are that it solves the tediousness of manual drawing and data processing, reduces the level requirements of data processing personnel, facilitates large-scale promotion, and greatly improves the efficiency of data processing, providing convenience for the optimal vegetation selection for ecological vegetation protection slopes.

[0114] 2) The present invention can process data using computer programs, enabling large-scale batch processing and significantly reducing the time cost of data processing. Furthermore, the present invention proposes data processing rules based on mathematical principles, effectively reducing computational errors caused by human factors during data processing.

[0115] 3) The data processing method proposed in the present invention has a simple calculation principle, is easy to operate on a computer, and requires relatively low computer computing power. This shows that the data automation processing method provided by the present invention has strong practicality and is easy to popularize and apply.

[0116] 4) The data processing method proposed in the present invention is not limited to plant root research, but is also applicable to the acquisition of stress-strain parameters in a large number of other physical tensile states, such as the testing of tensile stress-strain properties of vulcanized rubber or thermoplastic rubber, further broadening the scope of application of this method.

[0117] Related instructions:

[0118] In the process of plant-based soil reinforcement and stabilization, such as improving slope stability and strengthening the root-soil bond, the root system plays a primary role. Roots are the plant's nutritional organs that directly contact the external soil and form a root network with soil particles, thereby firmly anchoring the plant in the soil. The tensile strength of plant roots is one of the most important factors affecting the stability of the root-soil bond and is directly related to the size of the rhizome. The root's ductility before tensioning, elastic modulus after tensioning, and yield strength also play a crucial role in soil reinforcement. Existing research has largely focused on the relationship between root tensile strength and rhizome diameter. However, there are few clear reports on the tensile properties of roots (which primarily focus on ductility before tensioning, elastic modulus after tensioning, yield strength, and fracture strength). Therefore, further research on different root tensile properties is necessary to broaden the scope of applications for vegetation-based slope protection and stabilization, which has certain practical implications.

[0119] The tensile properties of plant roots need to be measured using equipment that can simultaneously measure tension and elongation. The present invention uses a universal material testing machine with the advantages of high measurement accuracy, wide range, and high test control as an example for measurement. The present invention does not limit the model of equipment that can simultaneously measure tension and elongation, and any equipment that can meet the test requirements can be used. The average diameter of the plant root system is measured by equipment such as a vernier caliper and a micrometer. Multiple measurements of different parts are taken to obtain the average value. There are also no restrictions on the measuring equipment and instruments. The effective initial length of the plant root system refers to the distance between the clamps at both ends of the root system after the plant root system is installed on the universal material testing machine. This length is measured using a ruler that can meet the accuracy requirements. There are also no restrictions on the measuring equipment and instruments. After measurement, the average diameter of the plant root system, the effective initial length, the tensile elongation and tensile force relationship data table, and the maximum tensile force can be obtained. These data can be calculated to obtain the stress-strain data sequence of the root system and the tensile strength of the root system.

[0120] The calculation formula is as follows:

[0121]

[0122]

[0123]

[0124] Where:

[0125] σ——Stress (MPa);

[0126] F——root tensile force (N);

[0127] D——average root diameter (mm);

[0128] π——pi;

[0129] ΔL——root tensile elongation (mm);

[0130] L——effective initial length of the root system, i.e. gauge length (mm);

[0131] ε——strain, i.e. elongation, dimensionless;

[0132] σ max ——Ultimate (maximum) stress, i.e. root tensile strength (MPa);

[0133] F max ——Maximum tensile force of the root system (N).

[0134] In the above data, the root tensile force F and the root tensile elongation ΔL are a series of corresponding data, which are obtained by a device that can simultaneously measure the root tensile force and elongation; F max is the maximum value of the root tensile force F; the average root diameter and the effective initial length of the root are directly obtained by the corresponding measuring equipment and instruments; the stress and strain data series with corresponding relationships can be obtained through equations (1) and (2); the root tensile strength can be obtained through equation (3).

[0135] First, the plant root to be measured is installed on a universal material testing machine that can measure tension and elongation simultaneously, and the average root diameter and the effective initial length of the root are measured using appropriate instruments (including but not limited to vernier calipers, micrometer screws, rulers, etc.) (after the root system is installed, clamps are fixed at both ends of the root system, and the length of the root system in the middle of the clamps). Secondly, the average root diameter and the effective initial length of the root system are input into the control computer of the universal material testing machine, and then the universal material testing machine is started to automatically collect the root tensile elongation and root tensile force, and automatically calculate the stress-strain data sequence according to formulas (1) and (2), and save it in the computer. Finally, different root systems are replaced and the operation is repeated to obtain a series of data files, each of which stores stress and strain data. There are two data cases in the data file, one column is strain data and the other column is stress data.

[0136] Using strain as the abscissa and stress as the ordinate, a root stress-strain curve is plotted in the first quadrant of a rectangular coordinate system. This curve generally exhibits a trend: a slow increase, a steep increase after the first turning point, a slow increase again after the second turning point, and finally a continuous increase to the end point. The following describes the eight parameters of the stress-strain curve and their representative meanings.

[0137] 1. Tensile strength: The end point of the stress-strain curve, also known as the breaking point, represents the maximum stress and strain that the root system can withstand. Its value is relatively easy to obtain and can be used to determine the breaking strength of the root system. The larger the strength value, the greater the tensile stress or tensile strain that the root system can withstand.

[0138] 2. Young's modulus: The slope of the steep increase section of the stress-strain curve represents the elastic modulus of the root system. Its value is generally obtained by graphical method. This straight line is called the elastic line and can be used to determine whether the root system is easily deformed under tension. The larger this value, the less likely the root system is to deform under tension.

[0139] 3. Crimp point: The intersection of the trend line of the steep increase section of the stress-strain curve and the horizontal line passing through the starting point is the virtual origin. Its value is obtained by graphical construction. The horizontal line passing through this point is called the virtual horizontal axis, and the vertical axis passing through the crimp point is called the virtual vertical axis. It can be used to determine whether the root system is in a tensioned state. The larger this value is, the stronger the ductility of the root system before tensioning.

[0140] 4. Lower boundary: The line connecting the imaginary origin and the fracture point represents the lower boundary of the stress-strain curve, which is called the lower boundary. The stress-strain curve is between this line and the vertical axis of the coordinate axis.

[0141] 5. Yield point: The point on the stress-strain curve farthest from the lower limit is the second turning point, also known as the yield point. It is the dividing point between the steep increase and the slow increase of the stress-strain curve, representing the yield stress and yield strain of the root system. Its value is generally obtained by graphical method and can be used to judge whether the root system is easy to yield. The larger this value is, the less likely the root system is to yield after being pulled.

[0142] 6. Upper boundary: The straight line passing through the yield point and parallel to the lower boundary is called the upper boundary, and the stress-strain curve is between this line and the lower boundary.

[0143] 7. Plastic modulus: As the tensile force increases beyond the yield point, the sample enters the plastic deformation stage. The slowly increasing section of the stress-strain curve is called the slow-increasing line, and its slope is called the plastic modulus or slow-increasing slope. The ratio of the plastic modulus to the elastic modulus is called the yield ratio. The yield ratio can be used to determine the degree of damage to the root system after yielding. The larger the value, the greater the damage to the root system after yielding.

[0144] 8. Toughness: The integral area of ​​the stress-strain curve is the toughness (J / m3), which describes the overall toughness of the sample.

[0145] An automatic processing system for stress-strain data of roots under tension, such as Figure 4 , including a data acquisition module 100, a data processing module 101 and a data output module 102;

[0146] The data acquisition module 100 is used to acquire stress-strain data files;

[0147] The data processing module 101 processes the stress-strain data file of the root system under tension using a file processing computer program to generate a data information file; uses the data processing computer program to sequentially read the stress-strain data file in a set order, and performs the following operations on each stress-strain data file: determines the strain data as the horizontal coordinate and the stress data as the vertical coordinate, obtains each data point of the stress-strain data file, and obtains the stress-strain curve by connecting multiple data points; loops through each data point of the stress-strain data file to obtain the maximum stress value, and the data point corresponding to the maximum stress value is the fracture point; and uses the stress-strain curve as the vertical coordinate. The first data point of each data point in the variable data file determines a horizontal line equation, and the straight line determined by the horizontal line equation is the virtual horizontal axis; the slope of the steep increase section of the stress-strain curve is obtained by the drawing method, and the slope of the steep increase section of the stress-strain curve represents the elastic modulus of the root system, and the straight line obtained by the drawing method is the elastic line; the intersection point of the elastic line and the virtual horizontal axis is obtained, and the obtained intersection point is the virtual origin, and the plumb line passing through the virtual origin is the virtual vertical axis; the dividing point between the steep increase and the slow increase of the stress-strain curve is obtained, and the dividing point is the yield point; the slow increase line of the stress-strain curve is obtained, and the slope of the slow increase line is the plastic modulus, and the ratio of the plastic modulus to the elastic modulus is the yield ratio;

[0148] The data output module 102 is used to output the obtained graph reflecting the stress-strain data changes of the root system under tension.

[0149] By applying the system of the present application, the yield point and fracture point positions in the root stress-strain data can be automatically obtained, which solves the tediousness of manual drawing and processing of data and the problem of inability to complete batch quantity processing in existing research, and avoids the influence of human factors on data calculation.

[0150] Example

[0151] The working folder is the project folder, which contains the data folder, the file processing computer program calc.sh, and the data processing computer program cals.exe. The data folder contains two data files: data file s1.txt (see Table 1) and data file s2.txt (see Table 2). Running the file processing computer program calc.sh outputs the data information file info.txt (see Table 3); then running the data processing computer program cals.exe outputs the result data file out.txt (see Table 4).

[0152] It should be noted that the first column of the s1.txt and s2.txt files is strain data, and the second column is stress data; the first line of the info.txt file indicates that there are only two data files in the data folder data folder, the second and third lines indicate that there are 199 lines of data in the data file s1.txt in the data folder data folder, and the next two lines (the fourth and fifth lines) indicate that there are 199 lines of data in the data file s2.txt in the data folder data folder; there are 3 lines of output in the result data file out.txt, the first line indicates the meaning of each column of data, and each subsequent line indicates the corresponding output result, where fname, o1x, o1y, ax, ay, bx, by, z, k0, b0, k1, b1, k3, b3, k4, b4, k6 and b6 in the first line correspond to the data file name, virtual origin, and so on. The horizontal coordinate, the vertical coordinate of the imaginary origin, the horizontal coordinate of the yield point, the total coordinate of the yield point, the horizontal coordinate of the fracture point, the vertical coordinate of the fracture point, the yield ratio z, the slope of the reference line k0, the intercept of the reference line b0, the slope of the elastic line k1, the intercept of the elastic line b1, the slope of the lower boundary k3, the intercept of the lower boundary b3, the slope of the upper boundary k4, the intercept of the upper boundary b4, the slope of the ramp line k6 and the intercept of the ramp line b6, the second and third rows represent the output results obtained after the data in the corresponding files are processed; the table lines of Table 1, Table 2, Table 3 and Table 4 in this embodiment are only for more convenient explanation of the contents of the embodiments of the present invention, and do not limit whether the data must be placed in the table lines; the values ​​in the contents of the s1.txt, s2.txt, info.txt and out.txt files are as shown in the corresponding tables. Such a data storage method is only for more convenient explanation of the contents of the present invention, and does not limit the data storage and reading format of the present invention.

[0153] Table 1 Data file s1.txt content:

[0154]

[0155]

[0156]

[0157]

[0158]

[0159]

[0160] Table 2 Data file s2.txt content:

[0161]

[0162]

[0163]

[0164]

[0165]

[0166] Table 3 Data information file info.txt content:

[0167] 2 s1.txt 199 s2.txt 199

[0168] Table 4 Results data file out.txt content:

[0169] fname o1x o1y ax ay bx by z k0 b0 k1 b1 k3 b3 k4 b4 k6 b6 s1.txt 0.00 0.00 0.11 75.59 0.42 132.46 0.25 287.45 0.00 860.22 -67.87 314.23 -24.70 314.23 17.47 212.52 37.94 s2.txt 0.00 0.00 0.11 42.94 0.37 72.58 0.22 186.07 0.00 464.11 -18.63 196.09 -7.82 196.09 13.70 103.48 29.02

[0170] It should be understood that those skilled in the art can make improvements or changes based on the above description, and all such improvements and changes should fall within the scope of protection of the appended claims of the present invention.

Claims

1. A method for automatically processing stress-strain data of a root system under tension, characterized in that: The following steps are involved: Processing the stress-strain data file of the root system under tension state by using a file processing computer program to generate a data information file; The stress-strain data files are read in sequence using a data processing computer program, and the following calculations are performed on each stress-strain data file: Determine the strain data as the abscissa and the stress data as the ordinate, obtain each data point of the stress-strain data file, and obtain a stress-strain curve by connecting a plurality of the data points; Looping through each data point of the stress-strain data file to obtain the maximum stress value, where the data point corresponding to the maximum stress value is the fracture point; Determine a horizontal line equation based on the first data point of each data point in the stress-strain data file, and the straight line determined by the horizontal line equation serves as a virtual horizontal axis; Using the line connecting the coordinates of the breakpoint and the coordinates of the first data point as a reference line; Using the general formula of the straight line equation of the reference line to establish a function using a linear programming method, and substituting all the data points into the function at once to obtain the corresponding function value; When the function value is greater than 0, the data point is located on the lower right side of the reference line; when the function value is less than 0, the data point is located on the upper left side of the reference line. Use the point-to-line distance formula to calculate the distance from each data point to the reference line in turn. The data point with a function value less than 0 and the largest straight-line distance from the point to the reference line is recorded as point n. The data point between the first data point and point n with a function value greater than 0 and the largest straight-line distance from the point to the reference line is recorded as point m. If there is no data point with a function value greater than 0 between the first data point and point n, then the first data point is recorded as point m; point m is the first turning point; a portion of data points between point m and point n is intercepted and recorded as array segment R1; array segment R1 is the steep increase segment data; the slope of the steep increase segment of the stress-strain curve is obtained using the stress-strain values ​​of the first data point p and the last data point q of array segment R1, and the slope of the steep increase segment of the stress-strain curve represents the elastic modulus of the root system; the regression line equation of array segment R1 is calculated using the least squares method using the first data point p and the last data point q of array segment R1, and the regression line corresponding to the regression line equation is the elastic line; Finding an intersection point between the elastic line and the virtual horizontal axis, the intersection point is the virtual origin, and the plumb line passing through the virtual origin is the virtual vertical axis; Obtaining a dividing point between a steep increase and a slow increase of the stress-strain curve, wherein the dividing point is a yield point; Obtain the slow-increasing line of the stress-strain curve. The slope of the slow-increasing line is the plastic modulus, and the ratio of the plastic modulus to the elastic modulus is the yield ratio. A graph reflecting the changes in stress-strain data of the root system under tension is obtained.

2. The method for automatically processing stress-strain data of a root system under tension according to claim 1, characterized in that: The process of using a file processing computer program to process the stress-strain data file of the root system under tension to generate a data information file includes: Establish a working folder at a set location in the computer, establish a data folder in the working folder, store all stress-strain data files of the root system under tension in the data folder, run the file processing computer program, read all stress-strain data files in the data folder, record the total number of data files, data file name and the total number of data rows of the corresponding data files, and write the records into the data information file in the working folder in the following manner: the first line of data in the data information file is the total number of data files, the subsequent two lines are the first data file name and the total number of data rows of the first data file, and the subsequent two lines are the second data file name and the total number of data rows of the second data file, until all data file name and the total number of data rows of the corresponding data files are written.

3. The method for automatically processing stress-strain data of a root system under tension according to claim 2, characterized in that: Use a data processing computer program to read the stress-strain data files in the set order: The data processing computer program reads the first line of the data information file to obtain the total number of data files in the data folder, reads the second and third lines of the data information file to obtain the file name of the first data file and the total number of data lines in the data file; According to the file name of the obtained data file and the total number of data rows of the data file, the data program reads the content of the corresponding file in the data folder and stores the stress-strain data sequence in the computer memory in the form of an array.

4. The method for automatically processing stress-strain data of a root system under tension according to claim 3, characterized in that: Reading the stress-strain data files in sequence according to a set order using a data processing computer program also includes a step of determining whether the data information file is empty, which step adopts: Run the data processing computer program to determine whether the data information file is empty. If it is empty, exit the program.

5. The method for automatically processing stress-strain data of a root system under tension according to claim 3, characterized in that: Use a data processing computer program to read the stress-strain data files in the set order: The data processing computer program determines whether the reading of the data information file is completed. If the reading is not completed, the data program continues to read the next two lines of the data information file, obtains the file name of a new data file and the total number of data lines of the data file, and performs a new round of calculations; if the reading is completed, the program exits.

6. The method for automatically processing stress-strain data of a root system under tension according to claim 1, characterized in that: The method for obtaining the dividing point between the steep increase and the slow increase of the stress-strain curve is: The data point where the function value is less than 0 and the straight-line distance from the point to the reference line is the maximum is recorded as point a, and point a is the yield point.

7. The method for automatically processing stress-strain data of a root system under tension according to claim 6, characterized in that: The method for obtaining the slowly increasing line is as follows: With point a as the center and the distances from point p and point q to point a as the radii, determine two concentric circles. All data points between point a and point b and located between these two concentric circles are recorded as array segment R2. The first and last data of array segment R2 are recorded as points u and v, respectively. The least squares method is used to calculate the regression line equation of array segment R2. The straight line corresponding to the regression line equation is a slowly increasing line.

8. An automated processing system for stress-strain data of roots under tension, characterized in that: It includes a data acquisition module, a data processing module and a data output module; The data acquisition module is used to acquire stress-strain data files; The data processing module processes the stress-strain data file of the root system under tension using a file processing computer program to generate a data information file; uses the data processing computer program to sequentially read the stress-strain data files according to a set order, and performs the following operations on each stress-strain data file: determines the strain data as the horizontal coordinate and the stress data as the vertical coordinate, obtains each data point of the stress-strain data file, and obtains a stress-strain curve by connecting multiple data points; loops through each data point of the stress-strain data file to obtain the maximum stress value, and the data point corresponding to the maximum stress value is the fracture point; Determine a horizontal line equation using the first data point among the data points in the stress-strain data file, with the straight line determined by the horizontal line equation serving as a virtual horizontal axis; use the line connecting the coordinates of the fracture point and the coordinates of the first data point as a reference line; establish a function using a linear programming method using the general formula of the straight line equation of the reference line, and substitute all the data points into the function at once to obtain corresponding function values; When the function value is greater than 0, the data point is located on the lower right side of the reference line; when the function value is less than 0, the data point is located on the upper left side of the reference line. Use the point-to-line distance formula to calculate the distance from each data point to the reference line in turn. The data point with a function value less than 0 and the largest straight-line distance from the point to the reference line is recorded as point n. The data point between the first data point and point n with a function value greater than 0 and the largest straight-line distance from the point to the reference line is recorded as point m. If there is no data point with a function value greater than 0 between the first data point and point n, then the first data point is recorded as point m; point m is the first turning point; a portion of data points between points m and n is intercepted and recorded as array segment R1; array segment R1 is the steep increase segment data; using the first data point p and the last data point q of array segment R1, the least squares method is used to calculate the regression line equation of array segment R1, and the regression line corresponding to the regression line equation is the elastic line; The intersection of the elastic line and the imaginary horizontal axis is obtained, and the obtained intersection is the imaginary origin. The plumb line passing through the imaginary origin is the imaginary vertical axis; the boundary point of the steep increase and the slow increase of the stress-strain curve is obtained, and the boundary point is the yield point; Obtain the slow-increasing line of the stress-strain curve. The slope of the slow-increasing line is the plastic modulus, and the ratio of the plastic modulus to the elastic modulus is the yield ratio. The data output module is used to output the obtained graph reflecting the stress-strain data changes under the tension state of the root system.

Images