Precision and Flexible Bending Control Strategy Based on Analytical Models and Data Models

Forming of various customized bending parts, small batches, as well as numerous types of materials is a new challenges for Industry 4.0, the current control strategies can not meet the precision and flexibility requirement, expected control strategy of bending processes need to not only resist unknown interferences of process condition and models, but also produce various new parts automatically and efficiently. In this paper, a precision and flexible bending control strategy based on analytical models and data models is proposed to build adaptive bending systems. New analytical prediction models for loading and unloading are established and suitable for various materials, a sequential identification strategy is proposed to search nominal properties using the four sub-optimization models. A data-based feedback model is established to prevent over-bending and eliminate online deviation. Above models are merged into a precision and flexible control strategy. The system firstly uses sub-optimization models to search the nominal point which is near to target point, secondly the system further uses feedback model to eliminate residual error between the nominal point and target point. Compared with four kinds sheet metals, the allowable ranges for variables are determined for a good convergence. The target bending angles were set to 20°, 40°, and 60°. Forty parts were tracked for each kind material, the adaptive bending system converged after one iteration, and exhibited better performances.


Introduction
Owing to flexibility of the air bending process [1], it is widely used in the bending of sheets, pipes, and profiles [2][3][4]. This approach can be used to form multi-shaped parts by controlling the punch stroke [5], and many scholars have proposed automatic control methods [6]. With development of Industry 4.0 and artificial intelligence, process is more closely integrated to production lines [7]. If parts change, the equipment can adjust the process parameters adaptively. To deal with the diversification of materials and part shapes, as well as the interference from production lines [8][9][10], novel flexible predication models are urgently needed.
At present, many scholars propose analytical models, which are based on constraint mechanical equations that can easily reflect the internal mechanism. Stelson et al. [11] proposed a method to determine the related parameters of the sheet using the force-stroke curve, and then used identified material properties to calculate the subsequently process parameters. Moreover, Stelson [12] developed a recognition system with a better real-time performance. Stelson [13] further studied the adaptive control for strain-hardened materials. Antonelli et al. [14] use a simplified model to conduct identification. Ton et al. [15] studied strip metal bending based on his model. Dallinger et al. [16] established an air bending control strategy through a two-step method based on a simplified model. Wang et al. [17] controlled the springback based on a simplified model. Presently, most scholars use simplified model to describe the loading process, which is difficult to precisely describe the behavior of various materials. If the bending deformation is not fully described by the hypothesis model, the control algorithm can not precisely and flexibly find the best working point.
In addition, the research on unloading model is rare. Actually, it is very difficult to thoroughly build bending model based on mechanical principles [18], and some assumptions have to be introduced. These assumptions and other unmodeled residuals cause deviations of prediction results [5]. Thus, many scholars have built models based on online data to predict bending process parameters [19]. Online data can include online interference, some is non-measurable, such as changes in lubrication conditions, tool wear, and equipment temperature rise, these can cause poor product consistency [20]. Many scholars use the feedback method to eliminate the deviation. Peter et al. [21] established a simplified control model for the blanking and bending processes based on the product chain, and realized the feedforward and feedback control of bending springback. Zafer [22] directly used experimental methods to generate data, and then performed springback control. Hazra et al. [23] studied the effects of changes in materials and processes on product formability. Some interference is measurable, such as material property fluctuation, if which can be evaluated in advance, the measurement data can be used to control the bending process. Therefore, many scholars used feedforward control methods to control the process. Havinga et al. [24] used bending force for feedforward control, and realized the consistency control of the product by product fluctuation. Yang et al. [25] used an online database and a fuzzy control model to regulate the precision bending process. Longo et al. [26] used control algorithms to predict bending springback in advance. Lewin et al. [27] used the feedforward control method to identify and control system. Yang et al. [28] developed an intelligent bending process based on a feedforward control system. Groche et al. [29] identified property based on statistical models, and then, controlled the bending process. From above analysis, it is relatively easy to fit a springback feedback model using online data, such as above feedforward and feedback models, But which need to accumulate a certain amount of data. In addition, when a model go beyond data boundary, the extrapolation accuracy of this model is poor [30,31]. If the target angle of control system is repeatedly adjusted within a wide range, a large amount of data collected from history record are required to fit the model, as a result, a lot of discarded parts are produced. So the disadvantages of the fitting model from online data are obvious, a feedback model based on data from the analytical model will be researched. In this paper, the first originality lies that new analytical prediction models for loading and unloading are established for applicability with various materials, and the sequential identification strategy is proposed to search nominal properties of sheet metal automatically, based on this fact, flexible control is realized. The second originality lies that a data-based feedback model is established to prevent over-bending and eliminate online deviation from optimized point, in this way, precision control is realized.
Research on precision and flexible bending control strategy is the focus herein. The advances of this article are as follows: Firstly, in response to the material diversity problem, novel prediction models are established to ensure the applicability, and a sequential identification strategy is proposed to search nominal properties of sheet metal to be manufactured for each of the four suboptimization problems, so flexible control can be done automatically. Secondly, a data-based feedback model is established to eliminate online deviation from optimized point, whose parameter within a certain safety margin is given to ensure that part is not over-bending, so precision control can be achieved steadily. Thirdly, the control strategy of the flexible system is finally determined as two steps.
Step one is that the system uses sub-optimization models to search the nominal point which is near to target point, then step two is that the system further uses feedback model to eliminate residual error between the nominal point and target point. Last but not least, the flexible system convergence effect using this control strategy is verified by experiments.

Flexible Hardening Model for Sheet Metal
To improve the adaptability of the flexible system, the material hardening model which have a wider description capability will be selected preferentially. From references, one type is constructed in the form of a power function, while the other is in the form of an exponential function. Among the power function type are the Hollomon and Ludwik models, whereas the exponential function models include the Voce, Voce+Voce, Voce++, and Hockett-Sherby models. where σ is true stress, ε is true strain, E is elastic modulus, ε s is yield strain, σ s is yield strength, K is strength coefficient, and n is hardening index. Referring to the Voce constitutive equation, the second recommended elementary function is determined as follows: where A and B are Voce model coefficients.
In order to test the hardening models, three representative curves are selected, as shown in Figure 1. High-strength steel DP980 have a higher yield stress, lower uniform elongation, and no evident yield turning point. The aluminum alloy 6061 sheets exhibits an apparent elastic-plastic transition. With the continuation of plastic deformation, the material hardening rate decreases slowly in a small range. For the 304 stainless ε ≥ ε s , steel material, there is a visible elastic-plastic transition curve, the hardening rate of the material, which remains larger than those of the other two, changes a little. Other test materials with similar behaviors are not listed there individually.
After fitting the stress and strain data of the material by the above two elementary functions, the fitting results are shown in the Figure 1(b)-(d). DP980 sheet metal is better fitted by Voce model than Ludwik model. The correlations of Voce model and Ludwik model with aluminum alloy 6061 sheet are similar. For the 304 stainless steel, it is same as aluminum alloy 6061 sheet.
When the elementary functions are used to extrapolate the high strain data, the Voce hardening rate changes little and tends to be saturated earlier, while the Ludwik model can continue to harden and be consistent with the cold working hardening behavior of the material. Therefore, when the material hardens slowly, the Ludwik and Voce models can be linearly combined to build a flexible hardening model: When α = 1 and α = 0, the combination function degenerates to elementary Eqs. (1) and (2), respectively. When α∈(0,1), the combination function can describe hardening behavior of diverse materials precisely. The parameter identification performance of different hardening models for various materials will be analyzed later.

Flexible Prediction Model in Loading Process
The combination constitutive model has been determined above. It is necessary to give the corresponding analytical model for the hardening model. In the forming process, as the bending stroke of the punch increases, the bending angle of the sheet gradually increases. The mechanical model is shown in Figure 2.
According to the mechanical model established in Figure 2, the static equilibrium equation and the geometric relationship in y direction are obtained as Among them, μ is the friction coefficient, θ is the rotation angle at the die fillet, P is the bending force, r d is die radius, r p is punch radius, N is the direct reaction force at the fulcrum, M is bending moment, l is the distance from the origin point in x direction to the tangent point of the sheet and die, L is the distance from the origin point to the center of the round corner on the die. The different mass points on curved section are in different deformation stages. ρ x is the curvature radius on the neutral layer at x, ρ lim is the limited curvature radius of elastic. When ρ x is larger than ρ lim , particle which is on the section of the sheet is in the state of complete elastic deformation.
When ρ x is less than ρ lim , particle which is on the section of the sheet is in the state of complete elastic deformation or plastic deformation. The two state calculation methods are as follows: 1) In the process of elastic deformation, the limited curvature radius of elastic ρ lim is obtained by generalized Hooke's law which is under the principal stress state: where ν is Poisson's ratio, y̅ is the ordinate value under the elastic limit. 2) In the process of elastic-plastic deformation, the deformation zone can be simplified to plane strain state, which means ε z =0. According to the assumptions, the stress in thickness direction during the air bending process is neglected, which means σ y =0.
And the deformation zone conforms to the small deformation assumption, which means the strain distribution on the section is linear. The bending moment of the sheet is calculated by the following formula: where M e is elastic bending moment, M p is plastic bending moment. According to the first elementary function 1 and the second elementary function 2, the combined bending moment calculation model is established as where y is the ordinate value of sheet cross section, b is the width of cross section , t is the thickness of the sheet metal, when α=1, the combined model degenerates to Eq. (1); when α=0, the combined model degenerates to Eq. (2); when α∈(0,1), the combined model can describe diverse relationships between bending moment and curvature. As the bending stroke of the punch increases, the radius of curvature under the punch gradually decreases. When the radius of curvature under the punch is ρ min >r p +t/2, the area closing to the punch of the sheet did not appear, as shown in Figure 3(a); when the radius of curvature under the punch is ρ min =r p +t/2, as shown in Figure 3(b), the area closing to the punch of the sheet appeared.
According to the situations whether the die is attached or not, the curvature of each point on the sheet is calculated by the following formula: Next, these control equations are used to establish the relationship between bending force and bending stroke during loading for each deformation state. Through a trial analysis, when the die span is larger, there is no monotonous relationship between the bending force and the stroke, in other words, when the stroke is increasing, the bending force which is corresponding to the next stroke increment may be larger or smaller than the current force. In order to simplify solution, a return back algorithm is used for calculation. At first, it is assumed that the bending force is larger than the current force, then the algorithm will find the corresponding stroke. If the stroke can not be found in the forward search direction, the search direction will be adjusted (9) backward. Until the convergence condition is meet, the algorithm exits the search program. This calculation process is shown in Figure 4.

Prediction Model in Unloading Process
The accurate descriptions of the bending loading and unloading process are important for the springback  prediction. In the elastic theory, it is assumed that the elastic unloading deformation is linear, so some highorder terms are ignored. But after certain plastic deformation occurs, these high-order terms will bring a high nonlinear effect on the unloading deformation of some materials, such as high strength steel. So in the unloading stage, the stress-strain nonlinear relationship of sheet metal is necessary to incorporate, and its typical unloading stress-strain curve is shown in Figure 5.
According to the reference, the elastic modulus is exponentially described as a function of the total strain, as shown in the following equation: where E u is elastic modulus during unloading, E 0 is initial elastic modulus, E a is saturated modulus.
It can be concluded that under the bending moment, the relationship between the springback curvature and the bending moment is shown as follows: where ρ u x is the radius of curvature at a point x after unloading the bending moment M. The amount of variation in curvature after springback can be further expressed as

Sequential Identification Strategy for Nominal Parameters
The prediction model of loading process which is based on combination elementary function and the prediction model of unloading process have been given. From above analysis, loading model and unloading model contain many parameters. Identification the nominal parameters using global optimization model usually fall into local optima, because the global optimization problem is non convex. So sub-optimization problems need to be decoupled from the global model, then the sequential identification strategy for nominal parameters needs to be determined. It is easy to obtain the full force stroke curve of the loading and unloading process from bending equipment, but different sub-optimization problems using different segments from full load curve. So the first task is to segment the load curve for each sub-optimization problem.
The force and stroke sensors collect data from the bending machine synchronously, The original data are recorded as The curves of the data are shown in Figure 6. The interpolation algorithm can be used to resample the experimental data at equal intervals to solve the problem that the experimental data and the predicated data from the mechanical prediction model are unequal on the abscissa. The interpolated experimental data points are recorded as (h e i , P e i ), i=1, 2,…, m int . The sampling interval is Δh, and the calculation formula is presented below. The force stroke data points given by the analytical model are recorded as (h a i , P a i ), i=1, 2,…, m int , and m int is the number of experimental data set after resampling.
Based on the above-described model, these experimental data points can be divided into four segments according to different purposes. The data points before yield bending are used to optimize the elastic modulus, and the transition points between elastic and plastic bending are used to optimize the yield strength, while the data points after yield bending are used to optimize the hardening parameters, and the unloading data points are used to optimize the variable modulus parameters. The window method [32] is used to automatically identify the (14) each segment boundary points of the experimental and predicated force-stroke curves. Finally, automatic curve segmentation is completed. The boundary points for each segment on the curve can be stably extracted by moving the window. The thresholds of the abscissa and ordinate are set to Δh c and ΔF c , respectively. In order to keep the threshold universal, the curve data are normalized, the number of window points is determined by the following formula: stroke are set as the ordinate, and the curve is shown in Figure 7 and Figure 8. It is evident from the curve that the loading and unloading of stroke curve is an isosceles triangle. It is shown from angle curve at the ath loading and unloading, there are three peak points, the points 1, 2, and 3 are shown in Figure 8. For the bth loading and unloading, point 4 is a valley point, while 5, 6, and 7 are peak points. For the cth and subsequent loading and unloading, the four boundary points are generated too.
The program determines the boundary points for each cycle (a-i). Finally, the segmented force curve is shown in Figure 9. The objective functions of the elastic, yield, plastic, and elastic unloading segments are defined as follows. When the predicated and experimental data are sufficiently close, in other words, each of the objective functions takes the minimum value, the selected nominal parameters are reasonable. Now four sub-optimization problems are built. They are optimized one by one sequentially, this is a sequential identification strategy for nominal parameters.
The nominal parameter ranges for each sub-optimization problem are set as where l e , l v , l p , l u are the number of nominal parameters that need to be optimized, and variables X e k , X y k , X p k , X u k are the degree of freedom. For each sub-optimization problem, a quadratic approximation function is used for each objective function. At the kth iteration, there is an approximate function shown as follows: where m∈[(l s +2), (l s +1)(l s +2)/2], s∈{e, y, p, u}. The optimization problems of above four objective functions are transformed into the minimizing problems of approximate functions. A optimization problem is shown as follows: where d is the vector step of each iteration and Δ k is the confidence radius at the kth iteration.

Feedback Model Using Online Data
The deformation control accuracy is affected by some assumptions introduced, mold misaligned and wear, so online closed-loop feedback mechanism is necessary. In the closed-loop bending control model, the target bending angle after springback is defined as α obj , the bending angle after springback in the last forming is defined as α i-1 , then the deviation angle of the previous forming is e i−1 , and it is expressed as There is a assumed nonlinear relationship between the bending stroke and bending angle after springback over a wide range. The spingback angle can be approximated by a linear relationship in a small neighbourhood of the nominal reference point h 0 , and the first-order Taylor is represented as According to Eq. (26), the relationship between the current forming deviation and the stroke increment is where K c is the control coefficient, the plastic bending stroke of the ith bending is The total plastic bending stroke of the ith bending is The bending angle after springback of the ith forming is where K b is the linear coefficient of bending model. The difference equations of the discrete control system are Eqs. (25)- (29), z-transformation is performed on the above-stated difference equation to obtain the transfer function, and the stability of the system is analyzed using these transfer functions. The transformed equations are as follows: (25)  The relationship between bending stroke and the stroke increment obtained from an integral element, as well as the transfer function of Eq. (30) are described as Eqs. (31)-(35) can be used to get the transfer function of the closed-loop feedback control system as The characteristic equation of the system is According to Eq. (37), the characteristic root of the system can be obtained as In the z-plane, the system satisfies the stability condition as The transfer function of control variable is shown as follows: Inverse z-transformation of Eq. (40), the difference equation used for the controller is According to the stability condition, when 0≤K c ≤2/K b , the closed-loop system converges steadily, different values of K c are selected, the output results of the system are shown in Figure 10. When the target bending angle of the springback control system is altered from 0° to 20°, that is increased by 20°, and K c =2/K b (PT 20-2/K b ), the output of the system is an oscillation curve with the same amplitude. When K c =1.5/K b (PT 20-1.5/K b ), some overshoots occur in the early stages of the system. After five adjustments, the output of the system meets the control requirement. When K c =0.9/K b (PT 20-0.9/K b ), the system converges to the steady-state value after one time. When K c =0.5/K b (PT 20-0.5/K b ), the system output is a rising asymptotic curve, and after five adjustments, the output of the system reaches the control requirement. When the target bending angle is altered from 40° to 20°, that is reduced by 20°, the output curve of the system with an angle increased by 20° is symmetrical to that with an angle reduced by 20°. Figure 11 shows the output of target bending angle 60°. It is evident that when the target bending angle increases, so do the oscillation and overshoot amplitude.
According to the above-stated analysis, when α i uo >α i-1 uo , a positive step is produced, and when α i uo <α i-1 uo , a negative step is generated. For the positive step signal, when the output of the system overshoots, the system will produce an over-bending phenomenon, which is irreparable and generates unqualified parts. When the output of the system is a gradual line, then multiple forming can be used for correction, however, too many corrections will reduce production efficiency. Negative step signals are contrary to positive ones, and the maximum overshooting value results in under-bending. Meanwhile, considering that the material properties of different batches and directions exhibit obvious fluctuations, system parameters with a certain safety margin are given to ensure that they are always under-bending during the Part number PT20-0.5/K b PT20-0.9/K b PT20-2.0/K b PT20-1.5/K b NT20-1.5/K b NT20-0.9/K b NT20-2.0/K b NT20-1.

Control Strategy of Flexible Control System
Flexible prediction model and sequential identification strategy are built separately. The control flow chart is shown in Figure 12. When the flexible system is used to bend a kind sheet metal for the first time, the system identifies the nominal parameters of the prediction model using sequential identification strategy automatically. The prediction model completes adjustment of the system working point. After the system arrives at the vicinity of the working point, the feedback correction model participates in compensation, and thereafter, interference and deviation are eliminated near the working point. The control strategy of the flexible system is determined.
When the system is adjusted in a large range, the nominal working point of the system is found by optimal model. When the nominal working point fluctuates in a small range, the system control accuracy is compensated by the data model. These promote flexibility and efficiency of the control algorithm for the air bending process.

Angle Measurement after Springback
To improve the accuracy and efficiency of angle measurement, a camera is used for online measurement. The grids shown in Figure 13 are created on the left-hand and right-hand sides of the image. The corresponding pixel coordinates are extracted at the abrupt pixel edge. The pixel coordinates extracted from the left-hand and righthand grid areas are linearly fitted by the least squares method to obtain the line equations, and the bending angle of the part is determined. Integrating the image processing algorithm into the computer control system program, the result is shown in Figure 13.
In order to verify the measurement accuracy of the online image measurement method used in this paper, the 3000iTM series of flexible three-coordinate measurement system produced by CimCore is used. Using PowerINSPECTTM measurement software and threecoordinate measuring tools to perform point-contact detection on the bending angle after unloading, the deviation is 0.01%. The image method used in this paper can meet the needs of experiments and engineering.

Flexible Control Experiment System and Softwares
In this study, a free V-bending flexible control experimental system is developed. The execution unit of the experimental system utilizes a servo electric cylinder, and a Delta 1.5 kW servo motor is used as the power source. The reduction ratio of the reducer is 1:12, and the screw lead is 10 mm. A grating ruler with an accuracy of 1 µm is installed on the movable beam of the machine, which form a fully closed-loop system with a servo drive. A control software is developed in the industrial control computer IPC-610H. The pulse control cards, pulse acquisition counting cards, and voltage acquisition cards are all connected through the PCI bus. The pulse control card is used for pulse generation to control the servo motor, while the pulse acquisition card is used for counting the pulse signal of the grating ruler, and the voltage acquisition card is used for the force sensor. The control software is written in the G programming language of LabVIEW. The software and hardware structures are shown in Figure 14. The mechanical system for air bending experiment is shown in Figure 15.
The control process of this study has been built into the computer program. When the geometric parameters of the material are adjusted in the interface, the program will identify the parameters automatically. The control accuracy of the experimental system controlled by computer is affected by various interference. In order to eliminate the above interference, after the identification is completed, the system enables the feedback control algorithm to correct the deviation.

Results of Online Flexible Control System
To confirm the accuracy of the above-mentioned method, high-strength steel sheets DP980, 304 stainless steel, ST12, and 6061 aluminum alloy were used in the verification experiment. The thickness of the DP980 sheet, 304 stainless sheet, ST12 sheet, and 6061 sheet were 2.2 mm, 0.8 mm, 1.0 mm, and 2.0 mm, respectively. A set of moulds was used to conduct an experimental study on the springback control of air bending. In this study, the material parameters self-identification optimization algorithm needs to explicitly provide the coefficient feasible ranges of the two constitutive models. The parameter ranges of these two models are determinable by fitting, where the parameter ranges of the   2000] and B∈ [5,200]. In the Ludwik model, the elastic modulus and the yield strength range are same, and the ranges of the parameters K and n are K∈ [10,5000] and n∈[0.001, 1], respectively. To verify the convergence of the algorithm, each parameter was optimized with three different initial values within the feasible range. The selection principle of the three initial values are that a point nears the upper boundary and another nears the lower boundary, and a random point is picked within the feasible range. The identification results of the elastic modulus of SUS304, ST12, AL6061, and DP980 are shown in Figure 16. It is evident from the figure that the elastic modulus is relatively stable during the iterative process, and the optimal approximate surface is constructed by several function value calculations in the trust region that rapidly converges to the target value through approximately 30 iterations. Although the elastic modulus of the four materials is different, the algorithm can always converge to a fixed value. The elastic modulus parameters, the length of the elastic phase, and the number of experimental points collected at a fixed acquisition frequency are different, which affect the calculation of the target function, and consequently influence the convergence efficiency of the elastic modulus, but not affect the convergence results. The yield strength parameter identification process is shown in Figure 17. Figure 17 apparently reveals that the SUS304, ST12, and Al6061 materials enter into the fast convergence stage after ten iterations. Furthermore, after approximately 60 iterations, their yield strength enter a stable convergence stage. The yield strength of AL6061 has a certain fluctuation, which is smaller than the yield strength of DP980. The total number of iterations of AL6061 materials is less than that of the other three materials, which have equivalent total iteration numbers. By comparing the four force-stroke curves, the DP980 transition from linear to nonlinear is longer, and the identification efficiency is relatively lower. The yield strength of the four materials converged to a fixed value after approximately 130 iterations. The identification processes of the SUS304 and ST12 materials are shown in Figure 18.  In the Ludwik model, there are two parameters that comprehensively control the hardening behavior of the material. Therefore, when choosing the initial value, K and n can be selected near the upper boundary, lower boundary or the middle. K and n fluctuates lightly in the early stage, and the subsequent learning process is stable.
By comparing the optimization processes of the two materials, for the initial value from the upper boundary, the SUS304 and ST12 materials converge approximately 240 and 160 iterations, respectively. In the case of the initial value from the lower boundary, both of the SUS304 and ST12 materials converge approximately 400 iterations, so the stop iterations of the initial value from the upper boundary are less than those from the lower boundary. For the Ludwik model, the initial value of the flexible system learning algorithm is recommended near the upper boundary for improving learning efficiency. Compared with the Ludwik model, the Voce model has a better saturation capacity and rapid hardening capacity in a shorter strain range. When the material hardens violently, the Voce model is preferred to describe the constitutive behavior of materials. The processes for identifying AL6061 and DP980 are shown in Figure 19.
There are two parameters in the Voce model that comprehensively control the hardening behavior of the material. It is evident from the formula that when the value of A is larger, the material hardens faster, whereas when the value of B is higher, the material saturates more rapidly, the material hardening rate is too high for larger A and B, but it is inconsistent with conventional materials. Therefore, when identifying the parameters of the Voce model, the value of parameter A is taken from the upper boundary, while that of parameter B is obtained from the lower boundary. in Figure 19, the case that parameters A and B are obtained from the upper and lower boundaries is compared with the other case that parameters B and A are got from the upper and lower boundaries, the total number of iterations of the former (300) is exceeds that of the latter (200). For the Voce model, the flexible system learning algorithm selects the initial value of parameters A and B near the upper and lower boundaries, respectively, to improve the learning efficiency. The above-stated recognition results are summarized in Tables 1, 2.
After the optimization of elasticity, yield, plasticity, and repeated unloading, the force-stroke curves of SUS304, AL6061, ST12, and DP980 are shown in Figure 20. It is evident from the figure that the bending load slightly decreases when the punch stops, that is, the sheet bending process exhibits a certain stress relaxation phenomenon. Among them, SUS304 and AL6061 sheets show a similar decline, while the DP980 sheet exhibit the greatest decline. During reverse unloading, a slight increase in the load. The phenomenon of stress relaxation affects the force-stroke curve for unloading bending. The starting point on the forcestroke curve for unloading bending coincides with the end point of the stress relaxation. The force-stroke curve for unloading bending has a different nonlinearity degree, among which DP980 has the largest nonlinearity degree. Figure 20(a), (b), and (c) reveal that the amount of springback deformation increases with the deformation, and the established variable modulus model has a higher prediction accuracy for springback. However, the nonlinearity of the springback of the DP980 sheet increases highly compared with the other three sheets, so the deviation of the springback prediction via the variable modulus model also increases. Established prediction models have unmodeled deviations, which need to be eliminated by the online feedback algorithm.
The optimization model can quickly find the nominal work point of the flexible system, but it cannot completely eliminate the effects of online interference. Therefore, it is necessary to introduce online methods to eliminate deviations. In order to compare the closedloop control effect, the open-loop control is also studied through one by one forming. Fluctuations were considered by cutting test specimens in the 0° and 90° directions of the sheet. To fully verify the correctness and the effect  of the algorithm, three kinds of bending parts are to be formed. The target bending angles are 20°, 40°, and 60°.
To determine the adaptive capacity of the system in the experiment, ten specimens were tracked.
For considering fluctuation of sheet property, specimens 1-3 and specimens 7-10 used blanks along the 0° direction, whereas specimens 4-6 utilized blanks along the 90° direction, and then, the adaptability of the system was evaluated. The parts were labeled after the experiment, the experimental data obtained are shown in Figure 21(a), bending parts for material SUS304 are shown in Figure 22.
Four groups experiments were conducted for each material: negative step feedback control experiment (G1-60-20), negative step no feedback control experiment   (G2-60-20), positive step feedback control experiment (G3-40-60), and positive step no feedback control experiment (G4-40-60). For each kind of four materials, it is apparent that if the flexible system is controlled with no feedback, there is always a deviation. For feedback control experiment, the bending angle of part number 1 after springback exceeds tolerance band, the bending angle of part number 2 corrected, the deviation is reduced. It can be concluded that the system has the ability to feedback out of tolerance. The results obtained through part number 1 and 2 experiments can converge quickly. The fluctuation also persists for blanks along the 0° direction from part number 1-3. Compare with part number 3 which is cut from 0° direction, part number 4 which is cut from 90° direction has a step interference caused by difference roll direction, the system can further resist this step interference.
To analyze the respond of the system to target angle step, step tests were conducted on the forming experiments of the four kinds sheets. For group 1, target angle of part number 6 is 60°, while target angle of part number 7 is 20°, so a 60°-20° negative step was occurred. In the same way, for group 3, a 40°-60° positive step was occurred. There were out-of-tolerance phenomena at the steps. To prevent the over-bending, the control parameter safety margin β was set to 10%. It can be seen from the curve that the feedback control model has a better control effect.
Furthermore, it is noticeable from the forming experiments of the four materials, materials ST12, AL6061, and SUS304 have approximately the same deviation after stabilization. However, DP980 has high-strength and large springback values, and the initial deviation of the analytical control model is larger, which implies that using online feedback methods is a better choice to eliminate deviation.

Conclusions
In order to precisely and efficiently form bending parts with various shapes and types, as well as small batches, a new control strategy of the flexible system is obtained. (1) In this new control system, novel prediction models are establish to ensure the applicability in diverse materials. The springback analytical model is built and incorporated the variable modulus function, a sequential identification strategy is proposed to search nominal material properties using the four sub-optimization models. A data-based feedback model is established to prevent over-bending and eliminate online deviation. Above models are merged into the new precision and flexible control strategy. (2) In feedback model, the coefficient K c is determined by Eq. (41) to prevent over-compensation. when α i uo >α i-1 uo , the input signal of the system is a positive step signal, and when α i uo <α i-1 uo , the input signal of the system is a negative step signal. Feedback model parameters with a certain safety margin β are given to ensure that part is not over-bending for positive step signal and negative step signal respectively. β is set to 10% that can meet the online requirement.
(3) In optimization models, compared with four kinds sheet metals, the allowable ranges for variables are determined for a good convergence. When choosing the initial value in the Ludwik model for optimization, K and n can be selected from the upper and lower boundaries at the same time or randomly selected from the middle. Moreover, in selecting the initial value in the Voce model for optimization. When parameters A and B are initialized from the upper and lower boundaries, respectively, the models converge easily. (4) In the experiment of control strategy, when the desired bending angles were separately set to 20°, 40°, and 60°. For each angle, the system firstly uses sub-optimization models to search nominal point which is near to target point, secondly the system further uses feedback model to eliminate residual error between nominal point and target point. Forty parts were tracked for each kind material. The adaptive bending system converged after one iteration, and exhibited better performances.