Remaining Useful Life for Heavy-Duty Railway Cast Steel Knuckles Based on Crack Growth Behavior with Hypothetical Distributions

The current research on the integrity of critical structures of rail vehicles mainly focuses on the design stage, which needs an effective method for assessing the service state. This paper proposes a framework for predicting the remaining useful life (RUL) of in-service structures with and without visible cracks. The hypothetical distribution and delay time models were used to apply the equivalent crack growth life data of heavy-duty railway cast steel knuckles, which revealed the evolution characteristics of the crack length and life scores of the knuckle under different fracture failure modes. The results indicate that the method effectively predicts the RUL of service knuckles in different failure modes based on the cumulative failure probability curves for different locations and surface crack lengths. This study proposes an RUL prediction framework that supports the dynamic overhaul and state maintenance of knuckle fatigue cracks.


Introduction
With the gradual development of railway vehicles for heavy-duty and rapid transport, the service conditions of crucial vehicle components are gradually deteriorating, and higher requirements for their service safety and reliability [1][2][3].The existing maintenance system for railway freight vehicles focuses on daily inspections, regular maintenance, and planned preventive maintenance [4].Although this system can ensure the safe operation of vehicles, it will result in excessive maintenance and unnecessary maintenance costs.The overhaul mode based on condition monitoring data is widely used in energy, aerospace, rail transit, and other industrial fields owing to its forward-looking nature, high overhaul efficiency, safety, and reliability [5,6].Condition monitoring of the critical components of heavy-duty railway vehicles is a passive operation conducted under harsh environments and complex conditions.Owing to the incomplete monitoring methods, the lack of process monitoring for internal damage to system structures has increased the difficulty of realizing fault diagnosis and identification based on the service status.
The heavy-duty cast steel coupler is crucial for connecting and transmitting longitudinal forces between the locomotive head and the vehicle and between two vehicles.Specific results have been reported for the dynamic characteristics of heavy-duty coupler systems and the analysis of the vehicle dynamics performance.The coupler is affected by the bending moment caused by the height difference of the coupler, the lateral force caused by the lateral swing and shaking of the head, and the longitudinal impact caused by the gap of the coupler [7,8].At the same time, the fatigue performance of the coupler structure is also an essential factor affecting life expectancy.The casting process produces random individual defects, such as surface pores and inclusions.The randomness and stress redistribution of the contact relationship caused by the wear of the contact part of the coupler, incomplete non-destructive surface crack testing, and uncertainty regarding the crack size result in a very high dispersion of the service life of the coupler [9,10].Establishing a prediction model for the remaining useful life (RUL) takes much work.To perform a life assessment based on the service state of the structure, a test study based on the equivalent load can effectively reflect the condition damage process of the structure and predict the RUL of the equivalent service state of the vital structure [11,12].In terms of predicting the RUL of the railway vehicle structure, simulation methods are frequently used to investigate the crack growth behavior and RUL under specific damage and crack conditions [13][14][15].In addition, it is difficult to establish a uniform crack growth model due to the service environment and the randomness of the crack state.
The critical components of heavy-duty railway vehicles are robustly designed and regularly maintained.When detected during inspections, cracks are repaired or replaced, and the residual value of damaged structures still needs to be effectively utilized.With the gradual improvement in testing and monitoring methods, more research has been carried out based on the reliability data analysis method for the life prediction of crucial components [16].In the study on life predictions based on coupler failure data, Yin et al. [17] obtained the curve for the relationship between the reliability of the coupler body and the service mileage through a reliability analysis of a small sample of failure data obtained through fatigue bench tests of the coupler.The manuscript only reports the life distribution of one fracture failure mode and does not describe its multi-site cracking.Ren et al. [18] evaluated the RUL of a heavy-duty coupler for a railway freight car with initial defects based on a fracture mechanics approach.The life assessment was limited to 10000 tons of load as a simulation, and the results were not verified experimentally.The RUL prediction of crack growth life in a complex mechanical system needs to be more thorough.In the reliability life assessment of section crack growth behavior, analytical methods and applications based on equivalent service states have not been developed to achieve life prediction under multiple failure modes [19,20], and an RUL prediction framework for the service crack growth process has not been established.
At the same time, the overall flow chart of the methodology used in this paper is necessary, as shown in Figure 1.A hypothetical distributed RUL prediction framework was proposed for visible and invisible cracks in service structures.The novelty is that the RUL distribution of different failure modes of complex structures is predicted based on the cumulative failure probability curve differently, and it is applied to the evolution of nondestructive testing cracks and life prediction under the equivalent service load of the heavy-duty railway knuckles.With new life data, the service mileage, crack location, crack size, and reliability are used as inputs to predict the RUL of the vital structural components based on the service status.

Equivalent Service Fatigue Life Crack Growth Test
Heavy-duty railway wagons have a variety of structural forms, and typical coupler structures include 16-type rotary cast steel couplers and 17-type cast steel couplers, both of which are made of E-cast steel [9,17].The couplers have the advantages of a small connecting gap, high structural strength, good interlocking performance, and high vertical anti-falling performance.They are currently mainly used on 70-ton and 80-ton heavy-duty railway vehicles.The 16H knuckle in the middle cavity connects the two couplers, and the traction pins at both ends of the knuckle are connected to the vehicle under the frame to realize the transmission of traction between the vehicles.The failure rate of coupler cracks continues to increase under the action of a large traction impact load.For the knuckle, the crack rate is as high as 60%-70%, generally occurring on the upper and lower traction platforms and the S traction surface of the knuckle.The composition of the coupler system, crack rates, and location of the knuckle during actual service are shown in Figure 2.
The service load studied in this work is 20000 tons of heavy-duty railway.Among them, the service load spectrum is derived from the coupler force obtained from the actual formation test, including operational factors such as traction braking during the vehicle formation operation.The bench test of the equivalent service life is carried out based on actual service conditions, considering the marshaling and coupler force during the actual service of the coupler.At the same time, the comprehensiveness of the test sample and the consistency of the distribution pattern of fault locations are ensured.The test samples and fracture failure patterns were statistically analyzed, as shown in Figure 3.
Two fracture failures were observed during fatigue bench testing of the knuckle: a fracture of the lower traction platform and a fracture of the S traction surface.The coexistence of the two parts of the cracks during the test  or the S traction surface under fatigue cyclic loading.These cracks expand to the size of the fracture crack, which can be observed from a macroscopic perspective.The morphology of the knuckle fracture is a typical fatigue failure fracture, with a crack source, a "beach strip" fatigue crack growth zone, and an instantaneous fracture zone.The typical crack growth fracture failure locations and shapes are shown in Figure 5.

Definition of the Crack Size Model
Based on the actual service status, it is necessary to conduct a fracture failure study of the fatigue crack growth behavior of the structure through equivalent tests, which can be used to predict the RUL of the structure.Considering the complexity of the knuckle structure of the heavy-duty railway and the defects of the casting process, the crack growth process of the knuckle structure is not a main crack.Under the action of fatigue load, the knuckle structure is accompanied by the initiation of several small cracks in the early stage, which finally merge into a main crack and propagate to failure.In order to facilitate the recording of the crack length and the statistical analysis for crack growth on the weak part of the knuckle (root arc of the upper traction platform arc, root arc of the lower traction platform, S straight line in the middle of the traction surface), the following four crack size models are defined.
Based on the multiple crack coplanarity and merging criterion in the safety assessment criteria for defectcontaining structures, the coplanarity criterion is applied for two neighboring cracks to merge them into one main crack.Multiple cracks in different knuckle regions were defined as four equivalent crack patterns [21,22].Model I: In the case of a single crack, the straight-line distance between the start and end of the crack is defined as the crack length.Model II: The distance between the two cracks is less than 5 mm, and the two cracks are merged into one crack.Model III: Two cracks are parallel (closer distance) and partially overlapped, the overlapped part and parallel distance are ignored.Model IV: If there are n collinear cracks and the distance between the cracks is greater than 5 mm, the crack sizes are recorded separately, and the total length of the crack is defined as the sum of n cracks.The equivalent crack size models at different parts are shown in Figure 6.

Equivalent Test Verification
The overhaul period of a knuckle is approximately 400000 km in service.First, the service load of the coupler was obtained according to the actual service state of the line, and then the bench test load spectrum was obtained.Based on the equivalent service load, the fatigue bench test of the new knuckle was carried out, and the fatigue cracks of the structure at different positions under the action of the fatigue load were monitored.The corresponding relationship between different crack lengths and the number of fatigue load cycles was recorded from crack initiation to final fracture failure.Finally, the life data of the knuckle of a limited sample were obtained.The mileage and bench test load spectra are equivalently converted to obtain the bench test knuckle after approximately 40000 fatigue cycles.Therefore, according to the crack sizes of different parts corresponding to the equivalent service mileage of the bench test, the crack information for 16 test knuckles in this cycle is obtained, as summarized in Figure 7.
According to the bench test load history of the knuckle in the equivalent service period corresponding to one repair period in Figure 7, the crack size of the lower traction platform is fitted to a normal distribution, and the data for the knuckle cracks in the actual service are obtained during a repair period.The fitting results are compared, and the mean, standard deviation, coefficient of variation, and correlation coefficient of the two fitted datasets are listed in Table 1.
The correlation coefficients of the two fitting models are more significant than 0.8, indicating that the normal distribution can satisfactorily fit the crack data of the test equivalent knuckle in service for a period of service.The coefficient of variation is relatively large, indicating that the data dispersion is relatively large.Based on the mean values, the test knuckle was used for a period equivalent to the maintenance period, and the length of the cracks on the lower traction platform was approximately 19 mm.In the survey data, the average crack length on the lower traction platform was approximately 18 mm, and the average value in the actual service was slightly lower.The main reason is that the load spectrum used in the test is equivalent to a relatively harsh knuckle position.In contrast, the knuckle position in the field investigation was random.In summary, the fatigue test of the coupler is conservative compared with the actual service condition of the line and can adequately meet the actual safety requirements of the project.Moreover, the test results can reflect the cracking condition of the service state of the knuckle to a certain extent.Therefore, the subsequent analysis is based on the crack growth test data of the knuckle to develop a prediction method for the RUL.

Hypothetical Distribution Evaluation Method
The hypothetical distribution is a statistical hypothesis of random variables, which can effectively describe the random distribution of system failure data and predict component life by fitting failure data under different failure modes [23].Based on the assumed distribution and improved model, many research fields have applied the parameter estimation method for failure samples under complex conditions [24,25].In mechanical system structure, the reliability index refers to the probability that the structure will complete a predetermined function under specified conditions and time.Commonly used reliability indicators include reliability, failure rate, and average working time.Assuming that the probability density function is f(t), the reliability R(t) is defined as the probability that the service life is greater than a certain specified value t, as shown in Eq. ( 1): The cumulative failure probability is defined as shown in Eq. ( 2): The key to calculating the structural reliability is obtaining the structural failure probability density function f(t), for large-sample problems (capacity of more than 50 or 70), as the probability density function is generally determined directly based on the data.For smallsample problems, the overall fitting effect, fatigue failure mechanism, and algorithm safety can be comprehensively considered.Commonly used statistical distribution types include the three-parameter Weibull distribution (3PWD), two-parameter Weibull distribution (2PWD), normal distribution (ND), lognormal distribution (LND), ( 1) (2)  a) The data were sorted from smallest to largest and the empirical failure probability was calculated based on the median P em (x i ): b) Based on linear regression, the common statistical distributions are used to fit the data (x i , P em (x i )) , and the statistical distribution parameters and fitting correlation coefficients R XY is obtained.c) For high reliability, the safety of the tail prediction can be evaluated by d F 1 and d F 2 .
For the left tail, the equations are as follows: For the right tail, the equations are as follows: (3) Considering the fitting correlation coefficient, fatigue failure mechanism, and safety of the tail prediction comprehensively [29].The above hypothetical distributions are used as the hypothetical distribution of crack length under a specific life score.The goodness of fit test is carried out on the crack size based on the hypothesized distribution.The linear correlation coefficients under each distribution function are compared to judge whether the hypothesized distribution is reasonable.The cumulative failure probability curve of the structure is obtained to realize the RUL prediction based on the overall characteristics.

Distribution Fitting RUL Results
During the knuckle equivalent service load fatigue bench tests, two crack growth behaviors occur under the same fracture failure mode: crack growth on the lower traction platform and the S traction surface.The RUL is defined as the period from the fatigue load cycle life of (6)   the structure with a certain crack size to the final failure.
Based on bench test data, the life score is defined as the ratio of fatigue loading cycles at the structure's current crack length to the cycles at fracture failure.The life score is distinguished for crack growth on these two parts under the fracture failure modes.The following four crack growth behavior modes are defined: Mode I: Lower traction platform cracks under a lower traction platform fracture; Mode II: Lower traction platform cracks under an S traction surface fracture; Mode III: S traction surface cracks under a lower traction platform fracture; Mode IV: S traction surface cracks under an S traction surface fracture.In order to facilitate the analysis of the fracture behavior, the surface crack lengths with identical life scores under the four modes are calculated through linear interpolation, as summarized in Table 2.
Taking the Mode I life score data (samples #1-#8) as an example, the calculated curve in the average crack length and the variation coefficient are shown in Figure 8.
According to the crack growth behavior of the lower traction platform in Mode I, the mean value increases continuously with the increase in the life score.However, the coefficient of variation steadily declines, moving from a strong to a weak variation.This demonstrates that the crack length exhibits considerable dispersion at the initial crack growth stage.However, the dispersion steadily decreases as the crack length increases.Combining the structure and measurement error of the crack definition and length during the growth process indicates that in the initial stage of crack growth, multiple small cracks are generated at the stress concentration area at the root of the traction platform, and the randomness of the locations and sizes of these tiny cracks leads to high dispersion.Under cyclic loading, the position of the fracture becomes relatively fixed as the minor cracks progressively merge to form the significant fatigue crack, and the   A comparative analysis using the following procedures should determine an appropriate hypothetical distribution: checking the total fit effects of possible assumed distributions on the crack length under a specific life score data.The approach allows consistency with the fatigue physics and checks tail fit effects.The six parameter distributions listed above are considered hypothetical distributions of the crack length under a specific life score and the linear correlation coefficients of each distribution function.Finally, the R XY values for each set of data are compared.The fitted linear correlation coefficients for the hypothetical distributions are listed in Table 3.
As shown in Table 4, the |R XY | value closer to 1 indi- cates a higher degree of fit.No data are available for life scores below, it cannot be assumed that their distributions are 3PWD, 2PWD, or LND.From the perspective of an |R XY | close to 1, EMVD2 provides the best fit for data with life scores of 0.1, 0.2, 0.3, 0.4, and 1.The ND function fitted best for data with a life value of 0.5, while EMVD1 fitted best for data with life values of 0.6, 0.7, 0.8, and 0.9.The data were further analyzed in terms of the safety of the tail prediction in order to obtain more reliable results.The deviations between the empirical value and the predicted value of the failure probability at the right tail of the six hypothetical distributions ( d F1 and d F2 ) are summarized in Table 4.
The results in Table 4 indicate that d F1 < d F2 means that as the reliability increases, there is a tendency to be more conservative.For d F1 < 0 , x > x n can provide a conservative estimate.For data that satisfy d F1 < d F2 , 3PWD has one group, 2PWD has one group, ND has four groups, LND has one group, EMVD1 has five groups, and EMVD2 has four groups.For data that satisfy d F1 < 0 , 3PWD has one group, 2PWD has one group, ND has five  groups, LND has one group, EMVD1 has five groups, and EMVD2 has four groups [29].
Considering the degree of fit and the safety of the tail prediction, EMVD2 can be used as a good hypothetical distribution for the Mode I failure crack data.Based on the above analysis method, the parameters of the cumulative distribution functions of the four crack growth behaviors and fracture failure modes under different life scores can be obtained, as summarized in Table 5.The corresponding cumulative distribution curves are presented in Figure 9.
Considering the degree of fit and the safety of the tail prediction, different distributions can be used as good statistical distributions under different failure modes.Based on the given crack surface length and cumulative failure probability of the different failure modes in Figure 9, the RUL of the knuckle can be predicted.Taking failure Mode I as an example, the specific method is as follows: assuming that a knuckle has been in service for k kilometers, the root crack length of the lower traction platform is determined through flaw detection to be 10 mm (corresponding to the transverse coordinates in Figure 9a).Then, for a cumulative failure probability of 0.1 (corresponding to the ordinate in Figure 9a), an evaluation point can be identified in Figure 9a, and this point falls in the middle of the curve with a life score of 0.6.This indicates that the service life is approximately 60% of the entire life, which is expected to be 1.67 k km, and the knuckle can thus still be used for 0.67 k km (k is a dimensionless unit).This conclusion has 90% reliability.The method for using the RUL cumulative failure probability curves under the other failure modes is the same.When high reliability is required (above 99.9%), a knuckle that has been in service for a repair period is tested for flaws.Once a crack (greater than 5 mm) is detected, the knuckle can be expected to survive the next repair period.

Validation of the Crack Growth Analysis
The crack growth analysis of a service-damaged knuckle selects two damaged knuckles under each of four failure modes.The crack growth curve of the new test knuckle  For the knuckle with crack growth on the lower traction platform, the knuckle in the lower traction platform failure Mode I is shown in Figure 10(a).The crack growth curve of the damaged knuckle tends to be flat in the first half of the life scores.The crack length ratio to the life scores rises gradually from the life fraction of 0.6 onwards, and the overall trend is consistent with the new knuckle in the equivalent bench test.Similarly, the failure Mode II in Figure 10(b) is the same as that for Mode I.For the third failure mode, as shown in Figure 10(c), the crack growth occurs on the S surface of the knuckle and the lower traction platform at the same time, the appearance of cracks in the lower traction platform leads to a reduction in the overall stiffness of the knuckle, and the crack growth rate of the S surface is relatively flat.For the cracks in the S surface of the knuckle and finally breaks in this mode, the crack growth location is more centralized in a plane, and the crack growth curve of the damaged knuckle is similar to that of the new knuckle.As shown in Figure 10(d), it can be analyzed that the crack growth curve of the damaged knuckle agrees well with the equivalent bench test data of the new knuckle.Overall, the crack growth curve of the damaged knuckle is almost entirely within the dispersion zone of the new knuckle test data, and it is consistent with the trend of the crack growth curve of the new knuckle.Therefore, these results verify that the fatigue bench test of the knuckle can reflect the actual service status of the knuckle to a certain extent, and the results of the crack growth analysis of the knuckle based on the bench test can guide the implementation of the condition detection and maintenance.

Remaining Useful Life Prediction Invisible Cracks
Section 3 realizes the prediction of RUL for different failure modes based on the knuckle non-destructive monitoring of crack growth life data.When there is an invisible fatigue crack in the knuckle, it is considered that the structure is in the stage before crack initiation, and damage has been generated relative to the new knuckle.However, existing flaw detection means that internal damage cannot be quantified.The crack rate of a knuckle in service for 400000 km is approximately 60%-70%.Once cracks are found, they will be scrapped.In order to evaluate the condition of the knuckle, a guideline for extending the maintenance period is proposed.Determining the RUL of a knuckle that has been in service for 30%-40% invisible cracks are also essential.Based on the relationship between the crack growth life data and service mileage of different knuckle samples under the equivalent service load above, the reliability curve for a knuckle with 400000 km of driving range and invisible cracks can be calculated as follows: considering the detection accuracy and dispersion of the surface defects of the structure, the whole process of structural service failure is divided into two stages.The surface crack length of 5 mm is defined as the crack initiation size, and the extension from 5 mm to critical failure is regarded as the growth life.The service mileage of the line corresponding to a surface crack length of 5 mm is obtained through linear interpolation, and the results are listed in Figure 11.The Grubbs test is used to ensure that there is no abnormal data in the life data for the knuckle crack initiation and growth life data [30,31].A two-parameter Weibull distribution is used to fit the life data.The fitting parameters of the life data are listed in Table 6.
The probability density function for the crack initiation life N i is The cumulative distribution function as: The probability density function for the crack growth life of the knuckle can be calculated by: (8)

The cumulative distribution function is
Based on the cumulative distribution function of crack initiation and propagation, the RUL prediction model based on delay time is constructed [32,33].It is assumed that the failure mileage of the knuckle is t × 10 4 km (t > 40).Event A is defined as follows: the knuckle does not detect cracks at 400000 km, and the crack initiation life of the knuckle is µ ∈ (40, t) .Event B is defined as fol- lows: the knuckle fails before time t, and the crack growth life of the knuckle is ν ∈ (0, t−µ) .P(B|A) represents the reliability of the knuckle that has been in service for (11) G N p = 1 − exp − N p 128.0482 3.9694 .400000 km invisible cracks are detected [34,35].Here, µ and ν are independent random variables.
According to the above equations, the variation in the reliability of the knuckle without detected cracks after service of 400000 km with the service mileage is shown in Figure 12.The service mileage parameters under different reliability levels are listed in Table 7.
It can be concluded from Figure 12 that the RUL of the knuckles that have been in service for 400000 km invisible cracks are detected under different reliability levels is less than the full life of the new knuckles.The RUL prediction model based on the invisible cracks can (12) reflect the difference in service performance between the old and the new knuckles.In addition, under different reliability levels, the full life of the knuckles that have been in service for 400000 km invisible cracks detected is greater than that of the new knuckles, indicating that the service performance of this knuckle is at the middle and upper level among the new knuckles.It can be seen from Table 7 that when the required reliability is greater than 99.9%, the remaining service mileage of the knuckle is less than 400000 km.It is recommended that a life assessment be carried out based on the crack sizes of the different parts of the steering knuckle detected during maintenance to improved knuckle utilization.

Discussions and Verification
Based on the framework proposed in this paper for predicting the RUL of in-service structures with visible and invisible cracks, the RUL of visible cracks at different degrees of reliability and the RUL prediction of the In Section 4, the reliability RUL prediction curves based on the delay time model are given for the invisible cracked condition of the knuckle.The reliable life of the knuckles that have served 400000 km of invisible cracks is worth the results shown in Figure 12.When high reliability (99.9% or more) is required, the knuckles that have been in service for 400000 km in one section of the repair period are flaw detected.As long as cracks (5 mm) are detected, they cannot meet the service requirements for the next repair cycle.In the case of the four fracture failure modes of the knuckle, Figure 13 compares the crack size states and the reliability RUL curves of the knuckle after 400000 km of service, respectively.The actual service site found that the knuckles in service 400000 km had a crack rate of about 60%~70%, for the remaining 30%~40% of the knuckles did not detect cracks.The RUL of the invisible crack knuckles at 400000 km in service is higher than the knuckle visible cracks, verifying that both methods can guide the life prediction of the actual service structure.Based on the RUL prediction framework proposed in this paper, a life assessment for the knuckle structure in service and utilizing a 400000 km cycle can be carried out.Combined with the crack size found by the actual flaw detection, the RUL prediction curve of visible crack state in different failure areas as shown in Figure 9, is adopted for the knuckle with crack state, and the RUL prediction curve of invisible crack state as shown in Figure 12 is adopted for the knuckle invisible crack state in one cycle of service.Combining the above two states can be selected in Figure 13 under the 99.9% reliability of different crack size curves to realize the dynamic RUL prediction.

Conclusions
A hypothetical distributed RUL prediction framework was proposed for visible and invisible cracks in service structures.Considering the structural component of the cast steel knuckles in heavy-duty railway vehicles as the research object, the failure behavior caused by fatigue crack growth is investigated based on the equivalent service load.Crack growth life data under different failure modes of the structure are obtained.The hypothetical distribution framework is employed to predict the statistical evolution characteristics of the crack length and life scores, and the RUL and maintenance suggestions are obtained for different areas, crack sizes, and reliability levels.The main conclusions are as follows: (1) Based on the crack growth behavior of the knuckle in service, an equivalent service bench test was carried out, and the crack growth life data and distribution rules of the S traction surface and lower traction platform of the knuckle were obtained under different failure modes.(2) A small amount of nondestructive crack growth life data obtained from equivalent bench tests were synthesized, which considered the overall fitting effect of the data, fatigue failure mechanism and algorithm safety.The cumulative failure probability curves for the crack length and life scores of the knuckle under different modes were obtained.These could be used to evaluate the RUL of a knuckle for a given crack size and reliability.(3) During the existing state maintenance cycle of the knuckle of 400000 km, once a 5 mm crack is detected, the RUL will not be able to satisfy the next maintenance cycle with 99.9% or more excellent reliability.It is recommended to combine the crack sizes monitored in the operation and maintenance of the actual service knuckle, and realize dynamic RUL prediction by selecting 99.9% reliability curves under different crack sizes according to the RUL prediction methods of visible and invisible cracks.
It should be noted that for the RUL prediction of the crack growth behavior at multiple areas in a complex structure, the crack growth process monitoring data is highly dispersed, and only the typical characteristics of the crack growth failure modes in different regions are considered in this study.Therefore, the current structural condition, mutual influence of multi-site crack growth, and condition process of individual differences should have been considered.The RUL prediction results from the universal application of the equivalent service state.In the future, data-driven algorithms such as neural networks and support vector machines can be used to expand the sample data, build a data-driven condition model based on physical characteristics, and further improve the accuracy of RUL predictions.

Figure 1 Figure 2
Figure 1 The overall flow chart of the methodology

Figure 3 Figure 4
Figure 3 Equivalent service load and test knuckle samples: (a) Heavy vehicle marshalling model (b) Heavy duty casting coupler structure, (c) Coupler force load spectrum, (d) Equivalent service fatigue bench test load spectra, (e) Non-destructive testing for cracks in service knuckles, (f) Knuckle test sample date distribution, (g) Knuckle test sample, (h) Fracture failure ratio of knuckle at different areas

Figure 5 Figure 6
Figure 5 Typical fracture failure modes of fatigue crack growth in the knuckle: (a) Fracture of the lower traction platform under the knuckle (b) S traction surface under the knuckle Figure 6 Equivalent crack patterns in different regions of the knuckle: (a) Equivalent crack Model I, (b) Equivalent crack Model II, (c) Equivalent crack Model III, (d) Equivalent crack Model IV

Figure 7
Figure 7 Length of crack on the surface of the lower traction platform under the test knuckle 150.0 140.0 141.0 130.0 140.0 140.0 110.0

Figure 8
Figure 8 Relationship between crack length and variation

Figure 9
Figure 9 Cumulative distribution curves for the crack length with different life scores

2 CrackFigure 10
Figure 10Comparative of life scores for damage knuckles

Figure 11
Figure 11 Equivalent life of knuckle crack growth results

Figure 13
Figure 13 Reliability curve for a knuckle in service for 400000 km: (a) Mode I (b) Mode II (c) Mode III (d) Mode IV

Table 1
Fitting parameters of crack length with normal distribution

Table 2
Evolution of the crack length on the lower traction platform with the life score

Table 3
Fitting correlation coefficients for the crack length under each hypothetical distribution

Table 4
The failure probability predicted differences of six hypothetical distributions

Table 5
Parameters of the cumulative distribution functions for four failure modes

Table 6
Two-parameter Weibull distribution fitting parameters for the knuckle crack life data

Table 7
12maining mileages with different reliability levelsFigure12Reliability curve for a knuckle service for 400000 km invisible crack state based on the delay time model is obtained for the knuckle equivalent bench-test fatigue life data in Sections 3 and 4, respectively.Based on the statistical distribution curves of the RUL of visible cracks in Section 3 above, the relationship between the cracks in different regions and the remaining service miles for different failure modes that have been in service for 400000 km can be derived.The RUL for the four methods with varying crack lengths and reliability is calculated based on the cumulative distribution function.The cumulative distribution curves of crack lengths for different life fractions of Figure9are taken to obtain the predicted values of RUL for crack sizes with arbitrary reliability and different failure modes.