1. Introduction
As an important part of support and transmission, a bar workpiece is widely used in mechanical products. The straightness of the bar workpiece not only affects the accuracy of assembly but also affects the performance of the equipment. Thus, improving the straightness of the bar workpiece has been widely studied [
1]. At present, straightening is an effective method to improve the straightness of rod parts. Straightening methods are divided into the following categories: roller straightening, stretching straightening, and bending-straightening, among others. Bending-straightening is widely used due to its high precision and suitability for various sections of the workpiece [
2].
The basic principle of bending-straightening is to make reverse the bending of the workpiece and produce elastic-plastic deformation. After the force is removed, the workpiece bounces back, finally having the straightness of the workpiece meet certain requirements [
3]. When straightening begins, the shaft is placed on two support blocks, usually at the midpoint of the two blocks where the deflection of the shaft is greatest, and the position of the straightening load is also at this midpoint [
4,
5]. When the initial deflection
at the midpoint is measured, the straightening system calculates the straightening stroke
. When straightening is complete, the instantaneous springback of the midpoint of the shaft is
. If Equation (1) is satisfied, the straightness of the shaft meets the requirements.
Research on bending-straightening has been developing for many years. Li et al. [
6] established a load-deflection bending-straightening model, studied the relationship between the load and deflection, and finally verified the validity of the model through finite element simulation and experiments. Zhao et al. [
7] developed a one-time leveling control strategy for multi-point bending and proposed a method to obtain the corresponding straightening parameters by discretization and linearization of theoretical torque curves. They also introduced load-correction coefficients. Finally, the validity of the proposed method was verified by finite element analysis and experiments. Mao et al. [
8] proposed a bending-straightening algorithm for seamless welded rail. The bending-straightening model was established by simplifying the section. In order to reduce the error caused by section simplification, the shape compensation coefficient was introduced, and the effectiveness of the proposed algorithm was verified through experiments.
Although there have been many studies on bending–straightening, in the process of developing a straightening machine, we found that, for slender shafts, there is a time-dependent springback phenomenon seen in the straightening process. Specifically, when the straightening actuator is separated from the shaft, the shaft will continue to spring back slowly; this has not been reported in the field of straightening. In
Figure 1,
is the time-dependent springback.
To date, some scholars have carried out research on time-dependent springback in other fields. Lim H et al. [
9] conducted draw–bend tests on different types of high-strength steel, and the time-dependent springback angle reached saturation after about one day. Sun and E. [
10] conducted a uniaxial tensile test on AC170PX aluminum alloy and studied the time-dependent springback of the sample under different tensile and strain rates. They found that the time-dependent springback was related to both rates.
The finite element method (FEM) is an effective tool, often used for various types of research, and the development of this approach is becoming increasingly refined. For example, Pinnola F P et al. [
11] studied the bending behavior of straight elastic beam systems with different scales by using the stress-driven, non-local, continuum mechanics method, and developed an effective calculation method based on a non-local, two-node, finite element method. At present, the finite element method is mainly used in the prediction of time-dependent springback. Liu and E [
12] used this method to predict the time-dependent springback of 1Cr18Ni9Ti stainless steel tubes. They adopted the one-dimensional idealization of the two-layer viscoelastic model [
13], fitted the parameters of the model by uniaxial tensile testing, and then used finite element software to predict the time-dependent springback. The predicted results are close to the experimental values. Fu and E [
14] adopted creep theory and then used finite element software for prediction. The prediction result was small, and the adjustment coefficient was introduced to make the prediction result close to the experimental value. Li et al. [
15] used creep theory and a viscoelastic model to predict the time-dependent springback of high-strength titanium tubular materials.
However, straightening is different from other processing methods, and the processing parameters are not the same. Secondly, although the finite element method generates results close to experimental values, they still have a large error due to the idealization of material parameters. Moreover, the finite element method is time-consuming and cannot be directly used in engineering applications.
There are few studies on the mathematical model of time-dependent springback. Using the theory of viscoelasticity, Sun [
16] predicted the time-dependent springback of DP600 under uniaxial tension. The predicted result was close to the experimental value and could well reflect the springback law. However, elastic-plastic deformation is involved in the straightening process, the mechanism of time-dependent springback is very complicated, and it is difficult to establish a mathematical model.
Although less attention is paid to time-dependent springback in the field of straightening, we found that the time-dependent springback had a great impact on the straightness evaluation of the slender shaft after straightening when we developed the straightening equipment. Considering the efficiency of straightening, a key procedure after straightening—the evaluation of the straightness of the shaft—will be carried out immediately after straightening. We connect the ends of the shaft as a reference. As shown in
Figure 2, we define the deflection value of the shaft as negative when the position of maximum deflection is lower than the reference, and positive when it is higher than the reference. Range ① is the evaluation standard without considering time-dependent springback. When the maximum deflection of the shaft is in range ①, the straightness error of the shaft can be judged to meet the requirements. However, after straightening, there will be time-dependent springback in the shaft, which makes the evaluation of straightness after straightening risky. After straightening, for example, when the straightness error of the shaft is in range ②, it is qualified without considering time-dependent springback, but time-dependent springback ultimately causes the straightness error of the shaft to exceed range ①. Thus, the straightness error does not meet the straightness requirements, and there will be
risk. If the straightness error of the shaft is in the range ③, it means the shaft is unqualified without considering time-dependent springback. However, due to the existence of time-dependent springback, the maximum deflection of the shaft is eventually within range ①, and the straightness error of the shaft finally meets the requirements. So, there is
risk. Therefore, the motivation of this paper is to propose a dynamic evaluation method of straightness considering time-dependent springback so as to reduce the two types of risk.
Therefore, understanding how to reduce the risk of straightness evaluation is critical. The key problem in this regard is determining the dynamic acceptability limit of straightness, and the key problem in determining this acceptability limit is predicting the time-dependent springback.
According to the literature review, there is no good method to predict the time-dependent springback at present, and it cannot be directly applied in engineering projects. However, machine learning [
17,
18,
19] methods are a better choice for predicting time-dependent springback. In particular, the BP neural network is widely used in this regard. Liu et al. [
20] used a BP neural network to predict the farthest distance and maximum velocity of flame propagation and achieved good prediction results. Wang et al. [
21] established a short-term solar irradiance prediction model based on the BP neural network and time series and verified its validity through simulation results. Liu et al. [
22] established a dynamic prediction model for the operating parameters of a nuclear power plant using a back-propagation neural network, and the verification results show that the model has high accuracy. However, the BP neural network has disadvantages such as slow learning speed and ease of falling into local minimum [
23,
24]. The genetic algorithm (GA) [
25,
26] has the ability of global searching and can make up for the shortcomings of the BP neural network. The GA-BP algorithm is also widely used in prediction and has achieved good prediction results. Yao et al. [
27] used a BP neural network optimized by the genetic algorithm to predict the foreclosure property market trend during the pandemic, providing an effective method for predicting housing prices. Wang et al. [
28] used a wind-speed prediction model based on improved empirical mode decomposition and the GA-BP neural network method and verified the accuracy of the model with MATLAB. The results were relatively accurate. He et al. [
29] established a prediction model for mold breakout prediction in slab continuous casting based on the GA-BP neural network and logic rules, and the prediction accuracy reached 100%. According to our literature survey, there is no relevant research that applies the GA-BP algorithm to the prediction of time-dependent springback in the straightening process. In this paper, the GA-BP algorithm is used to predict the time-dependent springback, offering a new idea for the research of time-dependent springback in straightening. However, as the model input is a key factor affecting the accuracy of the algorithm, no one has explored the influencing factors in the straightening process, which is also an important contribution of our paper.
Based on the above analysis, the purpose of this paper is to identify the factors affecting time-dependent springback in the process of straightening on the basis of mechanism analysis. We also propose a dynamic straightness-evaluation method considering time-dependent springback. The experimental results show that the genetic algorithm has high prediction accuracy and that the proposed method can effectively reduce the risk of misjudgment.
The innovative work of this paper is as follows: (1) Since there is no relevant research on the phenomenon of time-dependent springback in the process of straightening, there has been no systematic analysis of its influencing factors. So, based on bending-straightening and viscoelastic mechanics, influencing factors affecting time-dependent springback are analyzed on the basis of making certain assumptions. (2) Due to the phenomenon of time-dependent springback, there is a risk of misjudgment in the evaluation of straightness after straightening, so a dynamic evaluation method of straightness is proposed. Based on the GA-BP algorithm, a model to predict the time-dependent springback in straightening is established and compared with the linear regression model. Finally, the dynamic evaluation method of straightness is verified experimentally. (3) In order to explore the contribution of different influencing factors to time-dependent springback, the importance and sensitivity of said factors are analyzed.
The chapters of this paper are arranged as follows:
Section 2 analyzes the influencing factors of time-dependent springback in straightening;
Section 3 proposes a dynamic evaluation method of straightness based on
Section 2; In
Section 4, the dynamic evaluation method is verified experimentally;
Section 5 analyzes the sensitivity and importance of the influencing factors;
Section 6 presents the conclusion.
2. Analysis of Influencing Factors of Time-Dependent Springback Based on Mechanism Analysis
It is worth mentioning that the constitutive model of time-dependent springback is difficult to establish. The purpose of this section is to identify the parameters affecting time-dependent springback as far as possible through mechanism analysis on the premise of making certain assumptions. Thus, we aim to establish a model to quickly predict time-dependent springback based on GA-BP.
Figure 3 shows the Kelvin model, which is used in this paper to analyze the influencing factors of time-dependent springback.
E is the elastic modulus of the material,
is the viscosity parameter of the material.
is the stress of the model, and
is the stress of the elastic element, and
is the stress of the viscous element.
Assuming that constant stress
is loaded at
t = 0, where
is the time calculated from the loading start time,
is unit step function
The following relationship exists, where
is the strain of the model and
is the strain rate of the model.
The creep model of the Kelvin model can be obtained.
where,
The creep compliance
of the Kelvin model can then be obtained.
Figure 4 is the schematic diagram of shaft straightening loading.
When straightening, the load is usually carried in the middle of the two fulcrums. From material mechanics, we know that in the elastic bending stage, the relationship between deflection
and force
in the middle of the two fulcrums can be expressed as follows:
where
I is the inertia moment of the workpiece and
is half the fulcrum distance. According to the corresponding principle of elasticity and viscoelasticity [
30], when the shaft is subjected to a step force
, the deflection of the shaft can be expressed as follows:
Figure 5 shows the changes in force and deflection during the straightening process.
is the straightening stroke and
is instantaneous springback. In the loading stage, the shaft will first undergo elastic bending and then undergo elastic–plastic bending, and then the unloading process is elastic deformation. Phase 1 (loading procedure) can be expressed as follows:
Phase 2 (unloading procedure) can be expressed as follows:
The loading history (
Figure 6) of the shaft can be obtained from Phase 1 can be expressed as follows:
Phase 2 can be expressed as follows:
If the loading speed of the straightening actuator is
,
- (1)
Straightening loading stage
Equation (6) is only suitable for elastic bending, but not for elastic–plastic bending. Due to the complex mechanism of viscous behavior in elastic–plastic bending, in order to analyze the influencing factors, we make the following assumptions:
As shown in
Figure 6,
is the end time of straightening loading, and
is the end time of straightening unloading. The force in the loading stage is divided into countless
, and then during the
period, we think of the shaft as viscoelastic bending, in which the elasticity modulus
and the viscosity coefficient
vary in equal proportion, that is,
is constant.
According to the Boltzmann accumulate principle [
26], we can obtain the following:
where
, and
Substituting Equation (15) into Equation (14), we can obtain the following:
where
is also the time calculated from the loading start time.
- (2)
Straightening unloading stage
We can think of the unloading phase as reverse loading [
31], with the force increasing linearly from 0 to
. Then, the deflection of straightening unloading stage can be expressed as follows:
Therefore, the deflection change throughout the whole straightening process can be expressed as follows:
Thus, time-dependent springback can be expressed as follows:
However, there is a problem. Since the time-dependent springback analytical model is derived according to the creep model, during the straightening process of loading, we control the displacement of the straightening actuator precisely, which will inevitably generate errors. At the same time, the derivation process of the model also made some assumptions. Considering the effect of error compensation, we use Equation (21) to express time-dependent springback.
When evaluating straightness, we focus on the time-dependent springback after saturation, not the springback process itself. When
,
is the saturation time-dependent springback.
It can be seen from Equation (22) that, on the premise that the loading history of the straightening is determined, time-dependent springback is related to the straightening actuator speed , straightening stroke , unloading stroke , and viscoelastic parameter . The straightening actuator speed is related to the straightening time (time of the shaft-bending deformation process), . Therefore, we can preliminarily screen out the three influencing factors , , and .
According to the bending-straightening theoretical model [
2,
32], material parameters and fulcrum distance
L will affect the loading history of straightening. As can be seen between Equations (7) and (12),
also affects the straightening loading history. As shown in
Figure 6, the straightening loading history is closely related to force and
.
Since material parameters cannot be obtained in real time, but related physical parameters affected by material parameters can be, after comprehensive evaluation, in order to establish the fast prediction model of time-dependent springback, the following are selected as the influencing factors: , , , L, and .
3. Dynamic Method for Evaluation of Straightness Considering Time-Dependent Springback
3.1. Principle of Method for Dynamic Evaluation of Straightness
Based on the above analysis, a method for the dynamic evaluation of straightness is proposed.
Figure 7 shows the structure of the straightening actuator of the self-made straightening machine. The pressure head can detect separation and contact with the shaft, the force sensor can detect the straightening force in real-time, and the grating ruler can detect the displacement of the pressure head. The grating displacement sensor collects data on time-dependent springback, and it does not exist in the normal straightening process.
As shown in
Figure 8a, during straightening, the pressure head first detects contact with the shaft. At this time, the system records the position information
of the grating ruler and time
. When the straightening load is complete, the system records the position information
of the grating ruler and time
, as well as the force value
. Then, straightening and unloading begin, and when the pressure head is determined to be separated from the shaft, the system records the position information
of the grating ruler and time
.
These parameters can be obtained in real-time. At this point, as shown in
Figure 8b, the influencing factors analyzed in
Section 2 are taken as input in real-time, and saturation time-dependent springback
is predicted through the prediction model of time-dependent springback in real-time. Then, the acceptability limit of straightness can be directly obtained as follows:
The control system can then evaluate the straightening effect according to the acceptability limit in real-time. Therefore, the most critical problem is establishing the prediction model of time-dependent springback.
3.2. Time-Dependent-Springback Prediction Model Based on GA-BP Neural Network
3.2.1. BP Neural Network Prediction Model
The BP neural network is a typical feedforward neural network, which can approximate the nonlinear continuous functions with arbitrary precision. In this paper, a three-layer (including input layer) neural network is selected to establish the prediction model. According to the previous analysis, there are five main parameters that affect time-dependent springback. Therefore, there are five nodes in the input layer, namely, straightening stroke , straightening force , straightening time , fulcrum distance L, and instantaneous springback stroke . The output is saturated with time-dependent springback .
The transfer function between the input layer and the hidden layer adopts a sigmoid function. This is expressed as follows:
The transfer function between the output layer and the hidden layer adopts a purelin function. This expressed as follows:
Therefore, the prediction model can be expressed as follows:
where
is the weight matrix from the input layer to the hidden layer,
is the weight matrix from the hidden layer to the output layer, and
and
are threshold matrices of the hidden layer and output layer.
The number of nodes in the hidden layer has an important influence on the performance of the prediction model, and there is no relevant theory regarding the selection of the number of nodes in the hidden layer. The selection can be based on empirical formulas.
where
h,
m and
n are the number of nodes in the hidden layer, the number of nodes in the input layer, and the number of nodes in the output layer, respectively, and
a is a constant from 1 to 10. Thus, the number of nodes of the hidden layer is 3–13.
In this paper, the performance evaluation function of the BP neural network is the root mean square error (RMSE). The maximum number of training points is 1000, and the learning rate is 0.01, setting the training accuracy at 0.001. The Levenberg–Marquardt algorithm is used for network training.
However, the BP neural network also has disadvantages, including the fact that the weight and deviation of the network easily fall into the local minimum value, and there is no possibility for global searching. Thus, this paper uses the genetic algorithm (GA) to optimize the BP neural network.
3.2.2. BP Neural Network was Optimized by GA
The genetic algorithm (GA) simulates the genetic evolution process of organisms in nature and has the advantage of global searching for optimal values, which can make up for the disadvantage of the BP neural network in terms of easily falling into a local minimum value. The process of BP neural network optimization by GA is shown in
Figure 9.
Firstly, the structure of the BP neural network is determined;
The population is initialized, and then the weights and thresholds of the BP neural network are coded;
Populations are selected, crossed, and mutated;
The fitness values,
f are calculated to determine the fitness of individuals in a population. In this paper, the root mean square error (RMSE) of the predicted output and the target output of the BP neural network is selected as the fitness function, and the calculation formula is as follows:
is target output,
is prediction output;
We judge whether the stop condition of the GA algorithm is met. If so, the weight and threshold values after optimization are passed to the BP neural network. If the stop condition is not met, steps 3–5 is repeated.
3.3. Experimental Data Collection
In this paper, the shaft, which is commonly used by relevant enterprises, is a regular hexagonal section shaft; the maximum inner circle diameter is 12 mm; the material is 304 stainless steel.
As shown in
Figure 7, in order to collect data on time-dependent springback, a grating displacement sensor is placed at the midpoint of two fulcrum points below the shaft, and the displacement information on the grating displacement sensor is collected as soon as the shaft is separated from the pressure head. The resolution of the displacement sensor is 0.001 mm, and the uncertainty is ±0.3 μm.
In the process of bending forming, the greater the bending angle of the titanium tube, the greater the time-dependent springback angle, and the longer the saturation time of the time-dependent springback [
15]. Thus, in the straightening process, the greater the straightening stroke
, the greater the time-dependent springback and the longer the saturation time. Therefore, in order to determine the saturation time, the straightening speed was 20 mm/s, the fulcrum distance L = 665 mm, and the straightening of the stroke was 30 mm. The experimental results are shown in
Figure 10.
It can be seen that after 10 days of data collection, the time-dependent springback essentially reaches saturation, and the small springback amount after that has no great significance to engineering applications. Therefore, we chose 10 days as the saturation time for time-dependent springback.
In order to save time during data collection, we performed curve fitting on the data of the first 4 h from
Figure 10 to predict the saturation data after 10 days.
Figure 11 shows the predicted results.
The fitting curve was compared with the experimental data of 10 days, and the root mean square error (RMSE) was calculated as follows:
where
is the experimental value and
is the predicted value.
After calculation, the RMSE of the curve fitted with the 4-h experimental data and the 10-day experimental data is 0.003 mm. The saturation data predicted by the curve fitting after 10 days is 0.317 mm, which is only 0.002 mm different from the experimental value of 0.319 mm. Therefore, in this paper, in order to save on the cost of data collection, only 4 h of data were collected, and then MATLAB was used for curve fitting to predict saturation data.
In our straightening equipment, the fulcrum distances are fixed at 475 mm, 570 mm, 665 mm, and 760 mm. Therefore, under the condition that the straightening speed is 20 mm/s, the saturated time-dependent springback data (
Table 1,
Table 2,
Table 3 and
Table 4) are obtained by changing the straightening stroke
as follows:
3.4. The Establishment of Time-Dependent Springback Prediction Model
A total of 32 groups of experimental data were collected, among which 4 groups of data (
Table 5) were used as test sets to verify the generalization ability of the model.
Before model training, data needed to be normalized. The normalization process was carried out according to the following formula:
where
is the normalized value,
are the experimental data,
is the minimum value of the experimental data, and
is the maximum value of the experimental data.
After extensive testing, the number of nodes in the hidden layer is finally determined to be 8, so the structure of the BP neural network in this paper is 5-8-1.
The individual coding length of the genetic algorithm is , the population number is 20, the crossover rate is 0.4, the mutation rate is 0.1, and the number of iterations is 200.
Figure 12 shows the iterative process of the genetic algorithm. In around the 20th generation, the optimal interval between the weight and threshold value of the BP neural network is found. It can be seen that the genetic algorithm is able to continuously optimize the BP neural network, and it can compensate for the BP neural network’s disadvantage of easily falling into a local minimum value. Then, the weights and thresholds are transferred to the BP neural network for model training.
Figure 13 shows the training results of the BP neural network optimized by the genetic algorithm, and the regression value, R, represents the correlation between the output results and target values.
The regression values, R, of the training, verification, and test sets were 0.99837, 0.99952, and 0.99959, respectively, and the regression value of the whole set was 0.99867. This indicates that the correlation between the output value and the target value is high and that the model is good.
Figure 14 shows the predicted and target values of the test set. They almost coincide.
Figure 15 shows the error between the target value and the prediction value, and
Figure 16 shows the percentage error against the target value. As can be seen from the figure, among the four groups of test data, the maximum error is 0.0038 mm, and the maximum error percentage is 0.035%. This shows that the BP neural network optimized by the genetic algorithm has good generalization ability and high prediction accuracy.
3.5. Prediction Model Comparison
In order to further compare the performances of the prediction models, a multiple linear regression model was used to calculate the predicted values of the test sets. Using data from the training set and validation set for fitting, the regression model is shown in Equation (30). The degree of fit
was 0.9759.
The regression model was used to predict the test set data, and the prediction error and GA-BP model prediction error are plotted in
Figure 17. As can be seen from the figure, the prediction accuracy of the GA-BP model is much higher than that of the regression model. The maximum prediction error of the regression model is 0.014 mm, with a maximum error percentage of 0.127%, much higher than that of the GA-BP model. The RMSE and degree of fit
are used to further evaluate the accuracy of the prediction models. The formula for
R2 is as follows:
where
is the target value,
is the prediction value,
is the average of the target value. Two evaluation indexes were calculated using test set data.
As shown in
Table 6, the RMSE of the GA-BP model is 0.0042, and that of the regression model is 0.0098. The
of the GA-BP model is 0.9969 and that of the regression model is 0.9416. This shows that the relationship between input variables and output variables is not a simple linear relationship but a more complex nonlinear one. Furthermore, the GA-BP neural network has a good nonlinear fitting ability; it can learn the complex relationship between variables and is more suitable for the prediction of time-dependent springback in straightening.
At the same time, it also shows that the influencing factors screened in Chapter 2 can better predict the time-dependent springback. The model is successful.