Tuber crop skin breaking rate prediction method based on cumulative damage model

Abstract

The invention provides a tuber crop skin breaking rate prediction method based on an accumulative damage model, and the method comprises the steps: taking the force, speed, time, energy and other parameters of tuber crops in the mechanical harvesting and conveying process as the input; the instantaneous damage rate of the friction and abrasion channel and the instantaneous damage rate of the impact and contusion channel are calculated through the friction and abrasion channel and the impact and contusion channel respectively, meanwhile, the long-term friction and instantaneous impact effects are represented, and the mechanical mechanism of tuber crop skin breaking is completely described; based on the instantaneous damage rate, establishing an accumulated damage rate model, converting the tuber crop epidermis damage process into a quantifiable energy damage evolution model, and calculating an accumulated damage degree; calculating a predicted value of the skin breakage rate based on a probability mapping function; compared with the prior art, the method has the advantages that the experiment workload is obviously reduced, the tedious multi-parameter fitting process is avoided, and meanwhile, by optimizing the friction sensitivity coefficient and the impact sensitivity coefficient, the method can be used for predicting the skin breaking rate of various tuber crops.

Description

Technical Field

[0001] This application relates to the field of damage control and quality detection technology for tuber crops, and in particular to a method for predicting the skin breakage rate of tuber crops based on a cumulative damage model. Background Technology

[0002] Tuber crops, such as potatoes, sweet potatoes, carrots, and sugar beets, often suffer skin damage due to stress during mechanized harvesting, transportation, and sorting. Skin damage not only affects the appearance and commercial value of the crops but also increases disease infection and moisture evaporation, leading to a higher rate of spoilage during storage and severely restricting industry efficiency.

[0003] In recent years, scholars both domestically and internationally have conducted some research on crop bark breakage. Some international research teams have used field experiments and statistical analysis to establish empirical relationships between operating parameters and bark breakage rates, guiding equipment parameter selection. Other studies have analyzed the impact of single factors such as drop and collision on bark breakage based on finite element simulation. Domestic research has largely focused on optimizing the structure of harvester working parts and developing regression models for operating parameters, revealing to some extent the correlation between factors such as normal pressure, friction speed, and drop height and bark breakage rates.

[0004] However, these methods generally have the following limitations: on the one hand, the prediction results are difficult to generalize under different varieties or different equipment conditions; on the other hand, most of the existing prediction methods are empirical formulas or statistical models, and the reliability of the results is limited and they lack universality. Summary of the Invention

[0005] To address the aforementioned technical problems, this application provides a method for predicting the bark breakage rate of tuberous crops based on a cumulative damage model. The prediction method includes the following steps: Step S1: Obtain physical parameter data of tuberous crops during mechanized harvesting and transportation to construct a dataset. The physical parameter data includes normal force. Tangential force Relative slip velocity Duration of action Sliding distance Impact energy and observation of skin breakage rate ; Wherein, the tangential force ; The coefficient of friction; Step S2: Preprocess the data to obtain a standardized dataset. The preprocessing includes denoising, normalization, and dimensionality unification. Step S3: Based on the friction and wear channels and the impact damage channels, calculate the instantaneous damage rate of the friction and wear channels. and the instantaneous damage rate of the impact contusion channel Calculate the instantaneous damage rate; Step S4: Based on the instantaneous damage rate, establish a cumulative damage rate model and calculate the cumulative damage degree; Step S5: Calculate the predicted value of the skin damage rate using a probability mapping function.

[0006] In some embodiments of this application, in step S3, the instantaneous damage rate of the friction and wear channel is... The calculation method is as follows: Where A is the original damage coefficient of the friction and wear channel; Normal critical force; It is the critical tangential force; This is the critical value for velocity; This is the time critical value; This is the critical distance value; These are the power exponents of each parameter.

[0007] In some embodiments of this application, step S3 further includes the following steps: A friction sensitivity coefficient is constructed based on the instantaneous damage rate of the friction wear channel. The friction sensitivity coefficient Configured as follows: ; The friction sensitivity coefficient It is used to characterize the energy absorption and damage sensitivity of tuberous crops under continuous friction, and comprehensively reflects the combined effects of normal force, friction coefficient, relative slip velocity and time on frictional wear damage.

[0008] In some embodiments of this application, in step S3, the instantaneous damage rate of the impact injury channel is... The calculation method is as follows: Where B is the original damage coefficient of the impact contusion channel; To reach the critical energy level; For impact event functions; For the power exponent parameter.

[0009] In some embodiments of this application, step S3 further includes the following steps: An impact sensitivity coefficient is constructed based on the instantaneous damage rate of the impact injury channel. The impact sensitivity coefficient Configured as follows: ; The impact sensitivity coefficient Used to characterize the damage resistance of the epidermis of tuberous crops under instantaneous impact or drop, it is equivalent to including impact energy and contact deformation characteristics.

[0010] In some embodiments of this application, the method for calculating the cumulative damage degree in step S4 is as follows: .

[0011] In some embodiments of this application, in step S5, the skin breakage rate is calculated using a probability mapping function. The method for predicting the value is as follows: ; in, The probability mapping shape parameter is used to control the slope and response amplitude of the skin breakage rate output curve.

[0012] In some embodiments of this application, the method further includes the following steps: Based on observed skin breakage rate Fitting friction sensitivity coefficient Impact sensitivity coefficient The power exponents of each parameter , Probability mapping shape parameters The optimal solution.

[0013] In some embodiments of this application, based on the observed skin breakage rate Fitting friction sensitivity coefficient Impact sensitivity coefficient The optimal solution is obtained as follows: Obtain the observed skin breakage rate of the i-th group of experiments. ; Based on observed skin breakage rate Obtain the equivalent damage level for the corresponding group. ; Based on instantaneous damage rate Integrating over time yields the cumulative damage. ; ; in, ;; ; For the channel excitation amount of friction and wear channel, The channel excitation quantity for impact damage to the channel; Configure the equivalent damage level to be equal to the cumulative damage level, that is: ;in, Let be the channel excitation amount of the friction and wear channel in the i-th group of experiments; Let be the channel excitation amount of the impact injury channel in the i-th group of experiments; Based on observed skin breakage rate With predicted skin breakage rate To minimize the error, an objective function is established. ; Based on the objective function, the friction sensitivity coefficient is obtained after iterative convergence. Impact sensitivity coefficient The power exponents of each parameter , Probability mapping shape parameters Find the optimal solution and obtain the parameter set of the optimal solution. .

[0014] In some embodiments of this application, the following steps are also included: The parameter set of the obtained optimal solution Substitute into the instantaneous damage rate equation; The cumulative damage is obtained by integration. ; The predicted skin lesion rate for n groups is calculated using a probability mapping function, where the predicted skin lesion rate for the i-th group is... ; Based on the observed skin damage rate of n groups And predict the rate of skin breakage Calculate the average relative error and coefficient of determination ; in, ; ;in, This represents the average observed skin breakage rate across n groups of experiments. like ,and If the cumulative damage rate model meets the preset accuracy requirements, it means that the model meets the preset accuracy requirements.

[0015] Compared with existing technologies, this invention has the following advantages and beneficial effects: The method for predicting the skin breakage rate of tuberous crops based on a cumulative damage model in this application uses parameters such as force, speed, time, and energy of tuberous crops during mechanized harvesting and transportation as inputs. It calculates the instantaneous damage rate of the friction and wear channel and the instantaneous damage rate of the impact and contusion channel through dual channels, simultaneously characterizing long-term friction and instantaneous impact, thus fully depicting the mechanical mechanism of skin breakage in tuberous crops. Based on the instantaneous damage rate, a cumulative damage rate model is established, transforming the skin damage process of tuberous crops into a quantifiable energy damage evolution model and calculating the cumulative damage degree. The predicted value of the skin breakage rate is calculated based on a probability mapping function. This significantly reduces the experimental workload and avoids the cumbersome multi-parameter fitting process. Furthermore, by optimizing the friction sensitivity coefficient and impact sensitivity coefficient, this method can be used to predict the skin breakage rate of various tuberous crops.

[0016] It should be understood that the above general description and the following detailed description are merely exemplary and explanatory, and do not limit this document. Attached Figure Description

[0017] The accompanying drawings, which form part of this document, are used to provide a further understanding of the document. The illustrative embodiments and descriptions herein are used to explain the document and do not constitute an undue limitation thereof. In the drawings: Figure 1 This is a schematic diagram of the principle of a method for predicting the bark breakage rate of tuberous crops based on a cumulative damage model, provided in an exemplary embodiment of this application. Detailed Implementation

[0018] To make the objectives, technical solutions, and advantages of the embodiments of this application clearer, the technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application. It should be noted that, unless otherwise specified, the embodiments and features in the embodiments of this application can be arbitrarily combined with each other.

[0019] Tuber crops, such as potatoes, sweet potatoes, carrots, and sugar beets, often suffer skin damage due to stress during mechanized harvesting, transportation, and sorting. Skin damage not only affects the appearance and commercial value of the crops but also increases disease infection and moisture evaporation, leading to a higher rate of spoilage during storage and severely restricting industry efficiency.

[0020] In recent years, scholars both domestically and internationally have conducted some research on crop bark breakage. Some international research teams have used field experiments and statistical analysis to establish empirical relationships between operating parameters and bark breakage rates, guiding equipment parameter selection. Other studies have analyzed the impact of single factors such as drop and collision on bark breakage based on finite element simulation. Domestic research has largely focused on optimizing the structure of harvester working parts and developing regression models for operating parameters, revealing to some extent the correlation between factors such as normal pressure, friction speed, and drop height and bark breakage rates.

[0021] However, these methods generally have the following limitations: on the one hand, the prediction results are difficult to generalize under different varieties or different equipment conditions; on the other hand, most of the existing prediction methods are empirical formulas or statistical models, and the reliability of the results is limited and they lack universality.

[0022] Based on this, an exemplary embodiment of this application provides a method for predicting the skin breakage rate of tuberous crops based on a cumulative damage model. Using parameters such as force, velocity, time, and energy during the mechanized harvesting and transportation of tuberous crops as input, the method calculates the instantaneous damage rate of the friction and wear channel and the instantaneous damage rate of the impact and contusion channel through dual channels, respectively. This simultaneously characterizes long-term friction and instantaneous impact, comprehensively depicting the mechanical mechanism of skin breakage in tuberous crops. Based on the instantaneous damage rate, a cumulative damage rate model is established, transforming the skin damage process of tuberous crops into a quantifiable energy damage evolution model, and calculating the cumulative damage degree. The predicted skin breakage rate is calculated based on a probability mapping function. This significantly reduces the experimental workload and avoids the cumbersome multi-parameter fitting process. Furthermore, by optimizing the friction sensitivity coefficient and impact sensitivity coefficient, this method can be used to predict the skin breakage rate of various tuberous crops.

[0023] During mechanized harvesting and transportation, tuberous crops are prone to bark breakage due to external forces. Bark breakage is defined as the irreversible separation or rupture of the epidermal cell layer from the underlying tissue caused by external forces. Bark breakage can be caused by friction, shear abrasion, or impact damage.

[0024] An exemplary embodiment of this application provides a method for predicting the bark breakage rate of tuberous crops based on a cumulative damage model, such as... Figure 1 As shown, the prediction method includes the following steps: Step S1: Obtain physical parameter data of tuberous crops during mechanized harvesting and transportation to construct a dataset. The physical parameter data includes normal force. Tangential force Relative slip velocity Duration of action Sliding distance Impact energy and observation of skin breakage rate ; Among them, tangential force ; The coefficient of friction; during mechanized harvesting and transportation, tuber crops are mainly subjected to normal force. and tangential force The functions of the force are as follows: normal force primarily compresses cells, enhances frictional contact, and covers the area from light to heavy pressure; tangential force primarily tears the epidermis. Relative slip velocity. The friction power and strain rate are determined, covering both low-speed and high-speed conveying conditions; the duration of action... Sliding distance The cumulative effect of wear can be controlled by the conveyor belt length or test time; impact energy The potential for epidermal fracture caused by instantaneous impact can be quantified by varying drop heights or impact velocities; the rate of epidermal perforation can be observed. Record whether the sample epidermis is broken and the proportion of broken area, as reference data for model calibration, validation, and optimization. Completely cover multi-dimensional physical parameters related to broken skin to ensure the comprehensiveness of model input parameters and avoid prediction bias caused by single-factor analysis.

[0025] Data is collected using high-precision equipment such as mechanical sensors and image recognition systems, and processed synchronously using a unified timestamp and sampling frequency, which significantly improves the consistency and comparability of the data.

[0026] Step S2: Preprocess the data to obtain a standardized dataset. Preprocessing includes noise reduction, normalization, and dimensionality unification. Introducing preprocessing steps such as noise removal, normalization, and dimensionality unification effectively eliminates measurement errors and system biases, providing a high-quality standardized dataset for subsequent operations. Through real-time data acquisition and preprocessing, dynamic responses to undefined operating conditions can be implemented, enabling the prediction and control of skin breakage rates.

[0027] Step S3: Based on the friction and wear channels and the impact damage channels, calculate the instantaneous damage rate of the friction and wear channels. and the instantaneous damage rate of the impact contusion channel Calculate the instantaneous damage rate; To reflect the "threshold effect," that is, when the external force is below a certain critical value, skin rupture will not occur, and when the force exceeds this value, the risk of skin rupture increases, this invention introduces a critical parameter. , , , , , And dimensionless processing is achieved through normalization: , , , , , ; Based on the aforementioned mechanisms and critical parameters, to reflect the dual mechanisms of the skin-breaking process, this invention divides the force behavior of crop epidermis into two typical channels: continuous friction and instantaneous impact. In the frictional wear channel, factors such as normal force, friction coefficient, slip velocity, contact time, and slip distance occur simultaneously and are coupled with each other. Their interaction is expressed in a product form, reflecting the energy accumulation effect of multi-parameter parallel driving. In the impact damage channel, the skin-breaking behavior is controlled by a single energy input and uses an impulse function. The event trigger is characterized, and the response intensity is described in the form of an energy power. The two channels have an energy superposition relationship, indicating that frictional wear and impact damage can accumulate independently and jointly determine the epidermal damage rate. By decomposing the crop bark breakage process into two independent damage channels, frictional wear and impact damage, and by establishing instantaneous damage rate calculation formulas for the two channels, the damage mechanism of different types of external forces is fully characterized, which can improve the accuracy of prediction.

[0028] In step S3, the instantaneous damage rate of the friction wear channel The calculation method is as follows: Where A is the original damage coefficient of the friction and wear channel; Normal critical force; It is the critical tangential force; This is the critical value for velocity; This is the time critical value; This is the critical distance value; These are the power exponents of each parameter, used to describe the nonlinear interaction between the parameters.

[0029] In step S3, the instantaneous damage rate of the impact contusion channel. The calculation method is as follows: Where B is the original damage coefficient of the impact contusion channel; To reach the critical energy level; For impact event functions; For the power exponent parameter.

[0030] Thus, the instantaneous damage rate of the present invention is as follows: ; In this application, the mechanism of skin breakage in tuberous crops is divided into a frictional wear channel and an impact injury channel. The frictional wear channel adopts a continuous variable integral form, while the impact injury channel adopts a discrete event accumulation form, so as to form a unified instantaneous damage rate expression.

[0031] When the external force is less than the corresponding critical value, the damage rate tends to zero; when the external force exceeds the critical value, the risk of skin breakage increases significantly in a power-law manner, which is consistent with the actual mechanism of the skin breakage process. This ensures dimensional consistency at the mathematical level and reflects the complete evolutionary process of tuberous crop epidermis from continuous abrasion to localized rupture at the physical level, laying the foundation for subsequent cumulative damage integrals and probability mapping.

[0032] Considering that critical parameters are determined by the material properties of the crop itself, and that they cannot be measured individually in practical applications due to the interference of various factors, this invention proposes a friction sensitivity coefficient based on these critical parameters and the original damage coefficient. and impact sensitivity coefficient This is to avoid the problems associated with multi-parameter calibration.

[0033] Step S3 also includes the following steps: A friction sensitivity coefficient is constructed based on the instantaneous damage rate of the friction wear channel. Friction sensitivity coefficient Configured as follows: ; Wiping sensitivity coefficient It is used to characterize the energy absorption and damage sensitivity of tuberous crops under continuous friction, and comprehensively reflects the combined effects of normal force, friction coefficient, relative slip velocity and time on frictional wear damage.

[0034] An impact sensitivity coefficient is constructed based on the instantaneous damage rate of the impact injury channel. Impact sensitivity coefficient Configured as follows: ; Impact sensitivity coefficient Used to characterize the damage resistance of the epidermis of tuberous crops under instantaneous impact or drop, it is equivalent to including impact energy and contact deformation characteristics.

[0035] Thus, the instantaneous damage rate of the present invention is as follows: ; To ensure direct comparison of parameters under different experimental conditions, this invention adopts the International System of Units (SI) and normalizes each input parameter, transforming variables with different dimensions into dimensionless ratios, thereby establishing a stable computational basis. After processing, the dataset is as follows: Through the above steps, the various energy inputs and material responses during the skin breaking process are integrated into a unified equivalent parameter system, providing standardized dimensions and calculation basis for the subsequent calculation of cumulative damage.

[0036] Friction sensitivity coefficient and impact sensitivity coefficient As equivalent damage parameters, they can comprehensively reflect crop resistance and environmental conditions, facilitating cross-variety and cross-operating condition applications. This invention reconstructs the main force and energy parameters of the bark-breaking process into equivalent parameters to establish a unified modeling framework. Through mathematical induction and dimensional analysis, various force factors and structural response characteristics are integrated into a small number of equivalent parameters in the simplest form, thereby achieving comparable calculations across crops and operating conditions while maintaining physical meaning.

[0037] Each power exponent parameter can be obtained through large-sample statistics and fixed as a universal constant, thus eliminating the need for repeated measurements in every application. This is achieved through the friction sensitivity coefficient. and impact sensitivity coefficient The construction achieves parameter dimensionality reduction, which not only reduces model complexity but also significantly reduces the amount of experimental data required for model calibration. Both factors together form the overall response coefficient of crop materials to energy input. This invention not only achieves accurate prediction of the skin breakage rate of tuber crops such as potatoes but also has the ability to be extended across crops, and can be applied to predict the skin breakage rate of various tuber crops such as sweet potatoes and taro. Friction sensitivity coefficient and impact sensitivity coefficient This comprehensively reflects the characteristics of crop epidermal material properties, moisture content, and tissue structure. When predicting the epidermal damage rate of other crop types, only the friction sensitivity coefficient needs to be recalibrated. and impact sensitivity coefficient Using only two parameters significantly reduces the workload and avoids the tedious process of multi-parameter fitting.

[0038] Step S4: Based on the instantaneous damage rate, establish a cumulative damage rate model and calculate the cumulative damage degree; Based on the instantaneous damage rate mentioned above, the cumulative damage at any time can be obtained through integration. A cumulative damage model is then established, and the calculation method for the cumulative damage is as follows: .

[0039] By performing continuous integration calculations over time, the superimposed effect of damage from both the friction and wear channel and the impact and contusion channel is precisely quantified. Instantaneous fluctuation interference is eliminated through cumulative damage integration, improving the accuracy of crop damage rate prediction. The cumulative damage degree is constructed by integrating the instantaneous damage rate over time, dynamically reflecting the synergistic damage effect of the friction and wear and impact and contusion channels.

[0040] Step S5: Calculate the predicted value of the skin damage rate using a probability mapping function.

[0041] In step S5, the skin abrasion rate is calculated using a probability mapping function. The method for predicting the value is as follows: ; in, These are probabilistic mapping shape parameters used to control the slope and response amplitude of the skin abrasion rate output curve, adjusting the steepness of the skin abrasion rate growth curve with damage accumulation and the rate at which it approaches the saturation value, ensuring... The expression is monotonically increasing in the interval [0,1) and eventually approaches 1. The dataset is obtained as follows: .

[0042] By employing a probability mapping relationship, a quantitative correlation is established between the predicted value of the skin abrasion rate and the degree of damage, and the predicted value of the skin abrasion rate is calculated.

[0043] Based on observed skin breakage rate Fitting friction sensitivity coefficient Impact sensitivity coefficient The power exponents of each parameter , Probability mapping shape parameters The optimal solution.

[0044] In this invention, the friction sensitivity coefficient and impact sensitivity coefficient The optimal solution for the equivalent damage parameters is determined as follows: the working condition parameters obtained through experiments are substituted into the calculation formula for the cumulative damage degree to obtain the cumulative damage degree. Then, the theoretical skin breakage rate is determined by the calculation formula for the predicted skin breakage rate and compared with the measured skin breakage rate. The optimal solution for the parameters is obtained by maximum likelihood estimation or least squares method.

[0045] Based on observed skin breakage rate Fitting friction sensitivity coefficient Impact sensitivity coefficient The optimal solution is obtained as follows: Obtain the observed skin breakage rate of the i-th group of experiments. ; Based on observed skin breakage rate Obtain the equivalent damage level for the corresponding group. ; Based on instantaneous damage rate Integrating over time yields the cumulative damage. ; ; in, ; ; For the channel excitation amount of friction and wear channel, The channel excitation quantity for impact damage to the channel; Configure the equivalent damage level to be equal to the cumulative damage level, that is: ;in, Let be the channel excitation amount of the friction and wear channel in the i-th group of experiments; The channel excitation amount is the impact damage channel of the i-th experimental group; the excitation amounts of the friction and wear channel and the impact channel are included in the equivalent damage calculation to realize the quantitative characterization of the multi-physics coupling effect and expand the applicability of the model under complex working conditions.

[0046] Based on observed skin breakage rate With predicted skin breakage rate To minimize the error, an objective function is established. ; Based on the objective function, the friction sensitivity coefficient is obtained after iterative convergence. Impact sensitivity coefficient The power exponents of each parameter , Probability mapping shape parameters Find the optimal solution and obtain the parameter set of the optimal solution. By optimizing multiple parameters collaboratively, the bias problems caused by optimizing a single parameter can be avoided, thus improving the accuracy of predictions.

[0047] It also includes the following steps: The parameter set of the obtained optimal solution Substitute into the instantaneous damage rate equation; ; The cumulative damage is obtained by integration. ; The predicted skin lesion rate for n groups is calculated using a probability mapping function, where the predicted skin lesion rate for the i-th group is... ; Based on the observed skin damage rate of n groups And predict the rate of skin breakage Calculate the average relative error and coefficient of determination ; in, ; ;in, This represents the average observed skin breakage rate across n groups of experiments. like ,and This indicates that the cumulative damage rate model meets the preset accuracy requirements, demonstrating that the model can stably predict the skin breakage rate variation under different loads, speeds, and other conditions, indicating that the model structure is reasonable and the parameter solution is stable. A value below 0.95 but above 0.70 can still be considered as indicating that the model has basic interpretability of the data; when When the value is below 0.70, the parameters can be refitted.

[0048] The prediction method and model of this invention can be embedded in agricultural machinery operation systems. By collecting signals such as operating load, sliding speed, and contact time in real time, the current risk of crop damage can be calculated. Based on the predicted damage rate output by the model, the system can dynamically adjust operating parameters such as clamping force and roller linear speed to achieve adaptive control of damage risk. In addition, the model of this invention can be used as a simulation module in the design process, combined with 3D modeling and finite element analysis, to evaluate the impact of different structural parameters on crop damage, providing quantitative basis for optimizing operating components.

[0049] In this application, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that an article or device that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such an article or device. Without further limitation, an element defined by the phrase "comprising..." does not exclude the presence of other identical elements in the article or device that includes said element.

[0050] Although preferred embodiments of this application have been described, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments as well as all changes and modifications falling within the scope of this application.

[0051] Obviously, those skilled in the art can make various modifications and variations to this application without departing from the spirit and scope of this application. Therefore, if these modifications and variations fall within the scope of the claims of this application and their equivalents, the intent of this application also includes these modifications and variations.

Claims

1. A method for predicting the bark breakage rate of tuberous crops based on a cumulative damage model, characterized in that, The prediction method includes the following steps: Step S1: Obtain physical parameter data of tuberous crops during mechanized harvesting and transportation to construct a dataset. The physical parameter data includes normal force. Tangential force Relative slip velocity Duration of action Sliding distance Impact energy and observation of skin breakage rate ; Wherein, the tangential force ; The coefficient of friction; Step S2: Preprocess the data to obtain a standardized dataset. The preprocessing includes denoising, normalization, and dimensionality unification. Step S3: Based on the friction and wear channels and the impact damage channels, calculate the instantaneous damage rate of the friction and wear channels. and the instantaneous damage rate of the impact contusion channel Calculate the instantaneous damage rate; Step S4: Based on the instantaneous damage rate, establish a cumulative damage rate model and calculate the cumulative damage degree; Step S5: Calculate the predicted value of the skin damage rate using a probability mapping function.

2. The method for predicting the bark breakage rate of tuberous crops based on a cumulative damage model according to claim 1, characterized in that, In step S3, the instantaneous damage rate of the friction and wear channel The calculation method is as follows: Where A is the original damage coefficient of the friction and wear channel; Normal critical force; It is the critical tangential force; This is the critical value for velocity; This is the time critical value; This is the critical distance value; These are the power exponents of each parameter.

3. The method for predicting the bark breakage rate of tuberous crops based on a cumulative damage model according to claim 2, characterized in that, Step S3 further includes the following steps: A friction sensitivity coefficient is constructed based on the instantaneous damage rate of the friction wear channel. The friction sensitivity coefficient Configured as follows: ； The friction sensitivity coefficient It is used to characterize the energy absorption and damage sensitivity of tuberous crops under continuous friction, and comprehensively reflects the combined effects of normal force, friction coefficient, relative slip velocity and time on frictional wear damage.

4. The method for predicting the bark breakage rate of tuberous crops based on a cumulative damage model according to claim 3, characterized in that, In step S3, the instantaneous damage rate of the impact injury channel The calculation method is as follows: Where B is the original damage coefficient of the impact contusion channel; To reach the critical energy level; For impact event functions; For the power exponent parameter.

5. The method for predicting the bark breakage rate of tuberous crops based on a cumulative damage model according to claim 4, characterized in that, Step S3 further includes the following steps: An impact sensitivity coefficient is constructed based on the instantaneous damage rate of the impact injury channel. The impact sensitivity coefficient Configured as follows: ； The impact sensitivity coefficient Used to characterize the damage resistance of the epidermis of tuberous crops under instantaneous impact or drop, it is equivalent to including impact energy and contact deformation characteristics.

6. The method for predicting the bark breakage rate of tuberous crops based on a cumulative damage model according to claim 1, characterized in that, In step S4, the cumulative damage is calculated as follows: 。 7. The method for predicting the bark breakage rate of tuberous crops based on a cumulative damage model according to claim 6, characterized in that, In step S5, the skin breakage rate is calculated using a probability mapping function. The method for predicting values ​​is as follows: ； in, The probability mapping shape parameter is used to control the slope and response amplitude of the skin breakage rate output curve.

8. The method for predicting the bark breakage rate of tuberous crops based on a cumulative damage model according to claim 5, characterized in that, The method further includes the following steps: Based on observed skin breakage rate Fitting friction sensitivity coefficient Impact sensitivity coefficient The power exponents of each parameter Probability mapping shape parameters The optimal solution.

9. The method for predicting the bark breakage rate of tuberous crops based on a cumulative damage model according to claim 8, characterized in that, Based on observed skin breakage rate Fitting friction sensitivity coefficient Impact sensitivity coefficient The optimal solution is obtained as follows: Obtain the observed skin breakage rate of the i-th group of experiments. ; Based on observed skin breakage rate Obtain the equivalent damage level for the corresponding group. ; Based on instantaneous damage rate Integrating over time yields the cumulative damage. ; ； in, ; ; For the channel excitation amount of friction and wear channel, The channel excitation quantity for impact damage to the channel; Configure the equivalent damage level to be equal to the cumulative damage level, that is: ;in, Let be the channel excitation amount of the friction and wear channel in the i-th group of experiments; Let be the channel excitation amount of the impact injury channel in the i-th group of experiments; Based on observed skin breakage rate With predicted skin breakage rate To minimize the error, an objective function is established. ； Based on the objective function, the friction sensitivity coefficient is obtained after iterative convergence. Impact sensitivity coefficient The power exponents of each parameter , Probability mapping shape parameters Find the optimal solution and obtain the parameter set of the optimal solution. .

10. The method for predicting the bark breakage rate of tuberous crops based on a cumulative damage model according to claim 9, characterized in that, It also includes the following steps: The parameter set of the obtained optimal solution Substitute into the instantaneous damage rate equation; The cumulative damage is obtained by integration. ; The predicted skin lesion rate for n groups is calculated using a probability mapping function, where the predicted skin lesion rate for the i-th group is... ; Based on the observed skin damage rate of n groups And predict the rate of skin breakage Calculate the average relative error and coefficient of determination ; in, ; ;in, This represents the average observed skin breakage rate across n groups of experiments; like ,and If the cumulative damage rate model meets the preset accuracy requirements, it means that the model meets the preset accuracy requirements.

Images