Special Issue

## Optimization with Surrogate Models: Flow and Heat Transfer Applications

View this Special Issue

Research Article | Open Access

Volume 2018 |Article ID 8430578 | https://doi.org/10.1155/2018/8430578

Xiurong Fang, Jia Lu, Junfeng Wang, Jinhui Yang, "Parameter Optimization and Prediction Model of Induction Heating for Large-Diameter Pipe", Mathematical Problems in Engineering, vol. 2018, Article ID 8430578, 12 pages, 2018. https://doi.org/10.1155/2018/8430578

# Parameter Optimization and Prediction Model of Induction Heating for Large-Diameter Pipe

Guest Editor: Arshad Afzal
Revised06 Oct 2018
Accepted01 Nov 2018
Published26 Nov 2018

#### Abstract

The parameters of induction heating of large-diameter pipes have a direct effect on the final processing quality of the elbow, and the complexity of multifield coupling of magnetothermal force in induction heating can make it impossible to quantitatively optimize the design parameters of the induction heating device. In this paper, X80 pipeline steel induction heating is taken as the research object, and a corresponding numerical model is established. The influence of induction heating process parameters on the heating temperature of pipeline steel under the skin effect is determined. First, the influence of process parameters on the heating effect of pipeline steel is quantified by orthogonal test. Then, taking the optimum temperature difference between the inner and outer wall of X80 pipeline steel during the induction heating process as a target, the optimal process parameter set of the pipe induction heating is determined by using neural network genetic algorithm. Finally, comparing the relevant test criteria of the regression equation, the optimum mathematical prediction model of the outer wall temperature of the pipe induction heating process is obtained, which provides a theoretical basis for optimization of the process parameters of the pipe-based induction heating device.

#### 1. Introduction

The induction heating of pipe is one of the common techniques for production of elbows thanks to its rapid heating, easy cooling, small thinning rate of wall thickness, and easy adjustment of the bending radius of elbow. In the actual project of oil and gas transmission engineering, the bending of large-diameter pipes in oil and gas transportation pipe network is mostly processed by hydraulic push-type induction heating pipe bending machine. The system diagram is shown in Figure 1 [1], and the specific induction heating diagram is shown in Figure 2 [2]. It is required to realize continuous operation of propulsion→heating→bending→cooling→forming while the induction heating device fulfills the local heating of the pipe when the pipe bending machine is working [3]. At present, many researchers have conducted a lot of studies on the bending process of different materials, the mechanical properties, and the microstructure of the bent pipe [47]. However, the influence of various parameters on the heating effect of the induction heating forming process for large-diameter pipe is not yet clear. The detailed theoretical guidance and basis for the selection and optimization of design parameters on the induction heating device are not yet available. Thus, it is impossible to ensure a reasonable matching between the heating time of the pipe and the pushing speed of the hydraulic cylinder, resulting in the fact that the reliability of the pipe quality cannot be guaranteed [8]. Therefore, the induction heating device of the pipe has become one of the important factors affecting the final processing quality of the elbow. In this paper, the optimization of the process parameters of pipe induction heating is investigated by means of numerical simulation and mathematical methods.

#### 2. Analysis on Induction Heating Parameters of Pipes

The basic principle of induction heating is shown in Figure 3: the coil is connected to the heating power source, and an alternating magnetic field is generated in the pipe when the alternating current in the heating power source flows through the coil. Under the action of the alternating magnetic field, vortex current is generated inside the pipe, and the vortex current and alternating current in the induction coil have opposite directions with the same frequency. When the vortex current passes through the pipe, the pipe generates a large amount of Joule heat under the action of the material resistance, so that the pipe temperature rises rapidly to achieve the purpose of heating [9].

It can be seen that electromagnetic induction is used to generate vortex current heating in the conductor to realize the heating of the work piece. Therefore, by the electromagnetic induction law (1), the Joule-Electric law (2), and technical requirements of the large-diameter pipe induction heating forming device, we can obtain the main influential factors on the heating effect of the pipe, such as the induction frequency, current density, and air gap [10].where is the number of turns of the induction coil; Q is Joule heat (unit is J); i is the effective value of induction current (unit is A); R is conductor resistance (unit is Ω); and is the time of the current flowing through the conductor (unit is s).

Taking the elbow forming of X80 pipeline steel [11] with diameter Φ1219mm and wall thickness δ 27 mm as an example, the numerical modeling of induction heating is carried out.

##### 2.1. Mathematical Model of Induction Heating Electromagnetic Field

In the numerical calculation of finite element, the mathematical model of the induction heating vortex current field for the heated workpiece area can be obtained by the control equation (Maxwell equations) and the magnetic vector-scalar potential () theory of the induction heating electromagnetic field [1214]:where is Hamilton operator; μ is magnetic conductivity (H/m); σ is material conductivity (S/m); ε is capacitivity (F/m); is magnetic vector potential; Φ is scalar potential; and is induction vortex current density vector (A/).

In order to ensure the continuity condition of the field quantity at the interface of different media, consider the boundary condition as a ferromagnetic boundary with surface current density . The expression is shown as

##### 2.2. Mathematical Model of Induction Heating Temperature Field

It is assumed that the Joule heat generated by the vortex current in the alternating magnetic field is completely converted into its own heat; that is, the induced current Joule heat obtained by the electromagnetic field is used as an internal heat source to heat the workpiece, and the heating process is an unsteady heat conduction process [15, 16]. Based on Fourier’s law and the first law of thermodynamics [17], the differential form of the transient temperature field governing equation in induction heating can be derived:where is thermal conductivity (W/(m·°C)); T is temperature field distribution function; is intensity of the heat source of induction vortex current in the pipeline (W/m3); ρ is material density (Kg/); c is the specific heat capacity of material (J/(Kg· K)).

Use the Newton convection boundary as a temperature boundary condition:where refers to the normal direction on the boundary; is thermal conductivity (W/m·°C); h is convective heat transfer coefficient (W/m2·K); T is work surface temperature (K); and indicates ambient temperature (K).

##### 2.3. Coupling of Induction Heating Electromagnetic Field-Temperature Field

The induction heating is involved in the interaction of two physical fields of electromagnetic field and temperature field. The complexity of solving the magnetothermal coupling problem makes it difficult to obtain accurate solutions for many engineering issues [18]. Hopefully, the large-scale finite element numerical simulation software can be used for engineers to solve complex problems on the basis of theoretical analysis. This paper is aimed at the characteristics of induction heating magneto-thermal coupling of large-diameter pipes. We adopt ABAQUS software to solve the magneto-thermal coupling issue by sequential coupling [19]. First, the vortex current is obtained in the electromagnetic field; then the Joule heat of the workpiece under the action of vortex current field is got. Eventually, the Joule heat is used as the temperature field input condition to obtain the temperature distribution on the workpiece. The analysis process of magnetothermal coupling during the induction heating is shown in Figure 4.

##### 2.4. Building of Numerical Model

The geometric model of pipe induction heating consists of three parts: X80 pipe, induction coil, and air. Since the pipe is only involved in the solution of temperature field, it is assumed that, in the solution of the electromagnetic field, the physical properties of the material corresponding to the induction coil and air are not affected by the temperature change. Since the X80 pipe, induction coil and air together constitute the electromagnetic path in the induction heating process, the thermal property parameters of the relative magnetoconductivity, resistivity, specific heat capacity, thermal conductivity, density, etc. corresponding to the three parts should be considered in solving the electromagnetic field (as shown in Table 1). The material of induction coil is T3 copper, its resistivity is  Ω•m, and the relative magneto-conductivity is 1. As air is not affected by the current, the relative magnetoconductivity should be set as 1 when performing numerical simulation.

 Temperature(°C) Relative Permeability μ Resistivity ρ (10−6Ω·m) Specific Heat (J/Kg·°C) Thermal Conductivity (W/m·°C) Density ρ(Kg/m3) Emissivity (J/m2sK4) 0 200 0.195 450 47.4 7896 1 25 198 0.202 450 46.2 7896 0.998 100 194 0.251 479 43.7 7870 0.996 200 186 0.336 520 40.3 7849 0.992 300 179 0.434 564 38 7824 0.987 400 167 0.539 612 35.9 7796 0.983 500 156 0.659 668 33.8 7769 0.978 600 139 0.792 785 31.7 7744 0.978 700 93 0.947 936 28.4 7716 0.879 800 1 1.073 848 25.9 7691 0/648 900 1 1.164 655 23.5 7665 0.594 1000 1 1.199 620 23.9 7639 0.535

It is one of the important influential factors on the quality of the pipe to reasonably match the heating time of the pipe and the pushing speed of the hydraulic cylinder during the bending process of the large-diameter pipe (see Figure 1). Thus, in order to further study the problem of matching the induction heating time with the pushing speed, the finite element model is built in Cartesian coordinate system. At the same time, considering that the large-diameter pipe and its applied load are axisymmetric, 1/4 of the pipe is used to build a finite element model. The working loads are applied to the model according to the characteristics of each physical field; the corresponding analysis steps, solution variables, and boundary conditions are setup; different grid types are selected to mesh the model (the unit type of EMC3D8 is selected in the electromagnetic field and the unit type of DC3D8 is selected in the temperature field). In order to ensure the accuracy and convergence of the calculation, the pipe mesh is refined, as shown in Figure 5.

#### 3. Results Analysis

##### 3.1. Analysis of Numerical Simulation Results

In order to visually display the dynamic changes of the electromagnetic field and temperature field inside the pipe during the heating process, the main parameters such as current frequency (ƒ=1000Hz), current density (J=6 × 107A/m2), and air gap (d=30mm) are determined according to the engineering practice, and the distribution clouds of the Joule heat and temperature of the pipe under induction heating are obtained through simulation (see Figures 6 and 7), and the curves of temperature for the inner and outer wall of the pipe with the times are shown in Figure 8. It can be seen that the distribution of Joule heat and temperature on the pipe is that the outer wall’s is much greater than that of the inner wall due to the effects of skin effect, proximity effect, and ring effect in the induction heating process [20], which ultimately leads to the greater temperatures difference of the inner and outer walls of the pipe, causing a great inconvenience to the pushing speed of the elbow, which in turn affects the quality of the pipe bending.

It can be seen from Figure 8 that although there is unevenness in the temperature rise of the inner and outer walls in the pipe, the difference in temperatures between the inner and outer walls tends to be stable as the increase of heating time. Therefore, by parameter optimization, the difference in temperature gradient between the inner and outer walls of the pipe is reduced as much as possible.

In order to further explore the relationship between the current frequency, current density, air gap, and pipe heating effect, the control variable method is used to change the three process parameters for simulation calculation. The results are shown in Figure 9.

Figure 9 shows the effect of the current frequency, current density, and the air gap between coil and pipe on the temperature rise of the inner and outer wall of the pipe, respectively. It can be seen from the figure although that relationship between the variation of parameters and the temperature of the inner and outer wall of the pipe is basically linear; the influence of three parameters on the heating effect of pipe cannot be quantified. Moreover, the temperature difference between the inner and outer wall of the pipe is the key quantity; thus it is necessary to analyze the variation rule of the temperature difference between the inner and outer wall of the pipe to find the optimal parameter group.

##### 3.2. Parameter Optimization Based on Orthogonal Test Design
###### 3.2.1. Orthogonal Test Design for Induction Heating of Pipe

The three process parameters of the current frequency, current density, and air gap are the factors of the orthogonal test [21], and each factor is set at 5 levels, as shown in Table 2.

 Level Factor Current Frequency (Hz) Current Density (A/m2) Air Gap (mm) 1 600 5×107 10 2 800 5.5×107 20 3 1000 6×107 30 4 1200 6.5×107 40 5 1500 7×107 50

It can be seen from Table 2 that there is a total of 125 experimental parameter combinations, and the orthogonal test is mainly used to analyze the primary and secondary relationship of the influence of various process parameters on the heating effect. Therefore, according to the orthogonality and the engineering practice, some parameter groups with the characteristics of “uniform dispersion, neatness, and comparability” are selected as the representative ones for the experiment.

In the orthogonal experimental analysis, the Taguchi design method is used to neglect the interaction between different process parameters, the L25 (56) orthogonal table is selected, and the empty column is adopted as the error column. According to the 25 sets of parameters selected in Table 3, the numerical simulation analysis of the induction heating of X80 pipe (Φ1219mm, wall thickness 27 mm) is carried out to obtain the temperatures of the inner and outer walls of the pipe. The orthogonal experimental scheme is designed with the temperatures of the inner and outer walls of the pipe as the test indexes (see Table 3).

 Test No. Factor Test Index Current Frequency (Hz) Current Density (A/m2) Interval (mm) Outer Wall Temperature (°C) Inner Wall Temperature (°C) 1 1(600) 1(5×107) 1(10) 649 127 2 1(600) 2(5.5×107) 2(20) 618 125 3 1(600) 3(6×107) 3(30) 574 119 4 1(600) 4(6.5×107) 4(40) 528 114 5 1(600) 5(7×107) 5(50) 474 105 6 2(800) 1(5×107) 2(20) 587 121 7 2(800) 2(5.5×107) 3(30) 569 120 8 2(800) 3(6×107) 4(40) 546 117 9 2(800) 4(6.5×107) 5(50) 508 113 10 2(800) 5(7×107) 1(10) 1526 215 11 3(1000) 1(5×107) 3(30) 532 114 12 3(1000) 2(5.5×107) 4(40) 528 115 13 3(1000) 3(6×107) 5(50) 508 113 14 3(1000) 4(6.5×107) 1(10) 1420 208 15 3(1000) 5(7×107) 2(20) 1349 205 16 4(1200) 1(5×107) 4(40) 509 112 17 4(1200) 2(5.5×107) 5(50) 508 114 18 4(1200) 3(6×107) 1(10) 1319 201 19 4(1200) 4(6.5×107) 2(20) 1297 202 20 4(1200) 5(7×107) 3(30) 1248 200 21 5(1500) 1(5×107) 5(50) 510 115 22 5(1500) 2(5.5×107) 1(10) 989 173 23 5(1500) 3(6×107) 2(20) 1255 200 24 5(1500) 4(6.5×107) 3(30) 1253 202 25 5(1500) 5(7×107) 4(40) 1239 204
###### 3.2.2. Analysis of the Orthogonal Test Results of Induction Heating of Pipe

The variance method is used to determine the significant degree and contribution rate of the influence of various process parameters on the heating effect of medium-frequency induction heating of large-diameter pipes, and the primary and secondary relationship of the influence of various process parameters on the heating effect is obtained. Then we use the range analysis to verify the primary and secondary relationships again to provide a guidance for engineering process design.

(1) Variance Analysis. The analysis results of the influence of each parameter on temperatures of the inner and outer walls of the pipe are shown in Tables 4 and 5.

 Sources of Variance Degree of Freedom Sum of Squares Mean Square F Contribution Rate (%) Significance Current Frequency 4 0.996 0.25 31.63 21.28 Current Density 4 1.49 0.372 47.32 31.84 Air Gap 4 2.1 0.525 66.68 44.86 Error 12 0.095 0.008 -- 2.02 Total 24 4.671 1.155 -- 100
Note: indicates that the difference is significant; indicates that the difference is extremely significant.
 Sources of Variance Degree of Freedom Sum of Squares Mean Square F Contribution Rate (%) Significance Current Frequency 4 0.478 0.12 37.77 25.78 Current Density 4 0.589 0.147 46.47 31.72 Air Gap 4 0.751 0.188 59.26 40.45 Error 12 0.038 0.003 -- 2.05 Total 24 1.856 0.458 -- 100
Note: indicates that the difference is significant; indicates that the difference is extremely significant.

(2) Range Analysis. First, the comprehensive mean and range of each factor are obtained, and then the primary and secondary relationships of the influence of the process parameters on the heating effect are got by comparing the extremes [22]. The range analysis of the inner and outer wall temperatures of the pipe under different levels in the orthogonal test is shown in Table 6:

 Outer Wall Temperature (°C) Inner Wall Temperature (°C) Current Frequency Current Density Air Gap Current Frequency Current Density Air Gap K1 568.6 557.4 1180.6 118 117.6 184.8 K2 747.2 642.4 1021.2 137.2 129.4 170.6 K3 867.4 840.4 835.2 151 150 151 K4 976.2 1001.2 670 165.8 167.8 132.4 K5 1049.2 1167.2 501.6 178.8 185.8 112 R 480.6 609.8 679 60.8 68 72.8
Note: K1~K5 represent the average temperature of each factor at each level, and R is range.

According to the analysis results of Tables 4 and 5, when the temperature of the inner and outer walls of the pipe is taken as an index, the influence of the current frequency and current density on the index is significant, and the influence of the air gap on the index is extremely significant. And in the medium-frequency induction heating of the pipe, the influential degree of each process parameter on the heating effect of the pipe is air gap>current density>current frequency. From the range R corresponding to each process parameter in Table 6, it can be concluded that the influential degree of the process parameters on the heating effect of the pipe is consistent with the variance analysis result, which further verifies the accuracy of the obtained results. Compared with the range analysis, the contribution rate of each process parameter to the survey indicators is defined more clearly by the variance analysis.

It can be known from the engineering that when the induction heating temperatures of the outer wall of the pipe are 900~1000°C [23] and the difference between the inner and outer wall temperature is small, the final processing quality of the elbow is better. Combined with the above analysis, a better process parameter group can be initially determined as current frequency 1500 Hz, current density 5.5 × 107 A/m2, and air gap 10 mm.

##### 3.3. Process Parameter Optimization Based on Neural Network and Genetic Algorithm

The orthogonal experiment can determine the influence degree of the process parameters on the temperature of the inner and outer walls of the pipe and judge the better parameter group from the existing process parameter group. However, it is difficult to obtain the global optimal process parameters for the multiparameter group problems which affect the actual engineering. We use BP neural network method to establish the mapping relationship between the induction heating process parameters and the temperatures of the inner and outer walls of the pipe [24, 25]. The data is trained to predict the heating effect of all combinations of process parameters within the variation range, and the process parameters are optimized by genetic algorithm [26, 27]. Thus, the global optimal process parameter group is obtained.

###### 3.3.1. Prediction Model Based on BP Neural Network

The current frequency, current density, and air gap are selected as the network model input, and the outer and inner wall temperatures and of the pipe are taken as output to build a three-layer BP neural network model. The network structure is shown in Figure 10.

The data (see Table 3) obtained by numerical simulation is used as a training sample, and the sample is trained by running a BP neural network program. The training target is selected as 0.000001 and the learning rate is 0.1. The training of the network model is shown in Figure 11. When the iteration is 274 times, the training is finished with the system error of 9.75 × 10−7; thus its error is within the expected range, and the network output is more accurate.

In order to verify the accuracy of the training prediction results, the five new sets of process parameters are reselected as test samples for training prediction and numerical simulation. The training prediction results and numerical simulation results are shown in Table 7.

 Current Frequency (Hz) Current Density (A/m2) Air Gap (mm) Outer Wall Temperature (°C) Inner Wall Temperature (°C) Simulation Prediction Error Simulation Prediction Error 1200 6×107 30 883 862.5 2.3% 164 163.29 0.43% 1500 6×107 30 1043 1036.5 0.62% 183 185.35 1.3% 1000 5.5×107 30 640 638.2 0.28% 132 130.67 1% 1000 6.5×107 30 900 894.5 0.61% 167 163.40 2.1% 1000 7×107 30 1070 1079.3 0.87% 184 180.71 1.8%

According to the above results in Table 7, the error between the prediction data and the simulation data is within 3%, which indicates that the prediction accuracy of the BP neural network for the heating effect is higher, and the prediction result is close to the actual situation.

###### 3.3.2. Optimization of Induction Heating Parameters Based on Genetic Algorithm

The population with a size of 50 is selected and the binary code is adopted. The crossover probability and mutation probability are 0.75 and 0.01, respectively. The BP neural network prediction result is used as the fitness function to control the temperature of the inner and outer walls of the pipe [28, 29]. The genetic algorithm with 100 times of iteration is designed to optimize the process parameters of pipe induction heating. The software is used to iteratively calculate the BP neural network and genetic algorithm joint optimization program. After 100 iterations, the target fitness is shown in Figure 12.

It can be seen from the Figure 12 that, after 30 iterations, the fitness value of the temperature difference between the inner and outer wall under the processable condition is maintained at 751; meanwhile the temperature of the outer wall and the temperature of the inner wall is 923°C and 172°C, respectively. The process parameter group that is superior to the orthogonal test in the induction heating process of the pipe is obtained with the current frequency, current density, and air gap which are 1460 Hz, 5.65 × 107A/m2, and 15 mm, respectively.

##### 3.4. Pipeline Induction Heating Prediction Model Based on Regression Analysis

Through the aforementioned analysis and calculation, the selection of optimal parameter group for the induction heating parameter design can be realized, but In order to further guide the engineering practice, a prediction model of pipe induction heating is necessary. In this paper, the regression analysis of the data obtained from numerical simulation is carried out to determine a reasonable mathematical model for the prediction of the outer wall temperature, which provides a certain theoretical guidance for the actual processing. The data obtained from the numerical simulation of medium-frequency induction heating of large-diameter pipe are taken as samples, and the process parameters that influence the induction heating effect of pipe materials, such as current frequency ƒ (Hz), current density J (A/m2), and air gap d (mm), are seen as independent variables, while the outer wall temperature of the pipe is seen as the dependent variable. Subsequently, two regression models are established for prediction of the pipe outer wall temperature on the basis of combining the feature of sample data. It can be seen from Figure 9 that the process parameters and the temperature of the inner and outer walls show a linear change law; thereby a regression Model I is established. However, considering the magnitude difference of the selected values of each parameter, the logarithmic conversion of the sample data is performed to build a regression Model II [30].Regression Model I: Regression Model II:

Assume that the above two regression models meet the following conditions:

(1) , ,

(2) The selected variables are not random and have no multiple linear relationship with each other

(3) There is a linear relationship between the dependent variable and the transformed independent variable in the regression model II

The two regression models are analyzed and solved by a multiple regression calculation program, and the results are shown in Table 8.

 Model I: Model II: Regression Coefficients (0.525, 3.092×10−5, -17.184, -1049.604) (-33.103, 0.642, 2.061, -0.481) Sum of Residuals 181575.2825 0.295 Return Square Sum 3712511.687 4.915 Correlation Coefficient r2 0.953 0.943 F-test Value 129.515 120.585 P 1.029×10−20 2.23×10−19

The regression equations obtained by stepwise regression calculation are as follows:Equation of Model I: Equation of Model II:

It can be seen from Figure 13 that although the normalized residuals of Model I and Model II are basically normal distribution, comprehensive analysis of the calculation results of Table 8, including the correlation coefficient, sum of squared residuals, and F-test value. It can been seen that Model I has an optimal fitting degree; thus the regression equation (1) is chosen as the mathematical model for predicting the outer wall temperature of large-diameter pipe during the medium-frequency induction heating. In order to verify the accuracy of the regression models, five sets of process parameters are taken as test samples for the numerical simulation and regression model prediction. The results are shown in Table 9.

 Current Frequency (Hz) Current Density (A/m2) Air Gap (mm) Outer Wall Temperature T (°C) Simulation Return Error 1200 6×107 30 883 920.1 4.2% 1500 6×107 30 1043 1077.6 3.2% 1000 5.5×107 30 640 660.5 3.1% 800 6×107 30 683 710.1 4.0% 1000 6×107 40 626 643.2 2.7%

It can be seen from Table 9 that errors between the simulated data and predicted data are within 5%. The error value needs to be determined according to industry standards. However, it is worth noting that the mathematical model of the outer wall temperature prediction can accurately reflect the mapping relationship between the process parameters of the induction heating and the outer wall temperatures of the pipe. The mathematical model of the outer wall temperature prediction for large-diameter pipe in the induction heating process has certain significance for guiding engineering practice.

#### 4. Conclusions

(1) Based on the basic theory of induction heating and analysis of induction heating forming process for large-diameter X80 pipe, the magnetic-thermal coupling issue in the heating process is simulated by sequential coupling method. The influential law of each group’s induction heating process parameters (current frequency, current density, and air gap) on the temperature of the inner and outer walls of the pipe has been obtained.

(2) Taking the induction heating process parameters as the test factors, we carry out the orthogonal test design of the process parameters with the difference of the inner and outer wall temperature of the pipe as the evaluation index, the influential degree of each process parameter on the heating effect is determined as air gap>current density> current frequency, and the optimized process parameter sets for the heating bending of the pipe have been obtained using the neural network genetic algorithm, with the current frequency, current density, and air gap which are 1460 Hz, 5.65 × 107A/m2, and 15 mm, respectively.

(3) The linear regression models of two types of pipe temperature have been built based on the simulation results of different heating parameters. By comparing and analyzing the relevant test standards of each regression equation, we have obtained the mathematical model of the outer wall temperature prediction during the induction heating process of the pipe.

(4) From the bend forming practice of X80 pipe, the numerical simulation process parameters sets are somewhat limited. Hopefully, our mathematical analysis method is of a referential significance for practical engineering.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This work was supported by National Natural Science Foundation of China (NSFC) (no. 51775427)

#### References

1. R. X. Fang and W. S. Cao, “The research of hydraulic technology for large caliber medium frequency heating pipe hydraulic bending machine,” Chinese Hydraulics & Pneumatics, 2010. View at: Google Scholar
2. X. Fang, J. Wang, and W. Cao, “Induction heating of large-diameter thin-walled elbow based on skin effect,” Oil & Gas Storage and Transportation, vol. 36, no. 8, pp. 958–963, 2017. View at: Google Scholar
3. X. Fang, X. Meng, and H. Fan, “Research on Ultra-low and Constant Speed Hydraulic Control System of Large Push-bending Style Pipe Bender,” Machine Tool & Hydraulics, vol. 19, pp. 61–64, 2017. View at: Google Scholar
4. C. Pandey, M. M. Mahapatra, P. Kumar, N. Saini, and J. G. Thakre, “Nano-size Particle Evolution During Heat Treatment of P91 Steel and Their Effect on Micro Hardness,” Transactions of the Indian Institute of Metals, pp. 1–8, 2017. View at: Publisher Site | Google Scholar
5. C. Pandey and M. M. Mahapatra, “Effect of Heat Treatment on Microstructure and Hot Impact Toughness of Various Zones of P91 Welded Pipes,” Journal of Materials Engineering and Performance, vol. 25, no. 6, pp. 2195–2210, 2016. View at: Publisher Site | Google Scholar
6. C. Zheng, H. Liu, X. Ding, and Q. Fu, “Horizontal Vibration of a Large-Diameter Pipe Pile in Viscoelastic Soil,” Mathematical Problems in Engineering, vol. 2013, Article ID 269493, 13 pages, 2013. View at: Publisher Site | Google Scholar | MathSciNet
7. C. Pandey, M. Mahapatra, P. Kumar et al., “Some studies on P91 steel and their weldments,” Journal of Alloys & Compounds, vol. 743, pp. 332–364, 2018. View at: Google Scholar
8. S. Q. Pan, X. U. Zhe, and H. E. Bi-Shi, “A Study on Temperature Control for the Billet Electromagnetic Induction Heating Based on Simulation,” Industrial Heating, vol. 41, no. 6, pp. 25–28, 2012. View at: Google Scholar
9. H. Du, J. Li, and Y. Qu, “Mathematical Modeling of Eddy-Current Loss for a New Induction Heating Device,” Mathematical Problems in Engineering, vol. 2014, Article ID 923745, 7 pages, 2014. View at: Publisher Site | Google Scholar
10. Y. Wang, C. Liu, and X. Wang, “Influence of induction coil parameter on temperature field of steel plate and choosing of coil size,” Jinshu Rechuli/Heat Treatment of Metals, vol. 41, no. 4, pp. 178–182, 2016. View at: Google Scholar
11. T. D. S. Antonino, P. B. Guimarães, R. D. A. Alécio et al., “Measurements of the Thermophysical Properties of the API 5L X80,” Materials Sciences & Applications, vol. 5, no. 8, pp. 617–627, 2014. View at: Google Scholar
12. Y. H. Zhang, Numerical Simulation of Induction Heating Temperature Field, Jiangnan University, 2008.
13. Z. QianZhe, Numerical Simulation and Experimental Study on Induction Heating Process of Steel Bar Production Line, North China Electric Power University, 2013.
14. Z. Tianxu, Numerical analysis and experimental study on dynamic medium frequency induction heating process of high frequency straight seam welded pipe, Yanshan University, 2016.
15. Y. Chenguang, Numerical Simulation of Variable Power Induction Heating of 42CrMo Steel Shaft Parts, Yanshan University, 2010.
16. S. Wenjie, Numerical simulation and experimental study on electromagnetic field and temperature field of large-diameter molten silicon single crystal growth equipment, Zhejiang University, 2013.
17. J. L. Fleming and J. Moser, “A Fourier Approach to Model Electromagnetic Fields Scattered by a Buried Rectangular Cavity,” Mathematical Problems in Engineering, vol. 2009, Article ID 937327, 22 pages, 2009. View at: Publisher Site | Google Scholar
18. X. Yi, L. Li, J. Yue et al., “Simulation of coupling of electromagnetic and thermal of intermediate frequency induction heating of pipeline,” Heat Treatment of Metals, vol. 41, no. 7, pp. 154–158, 2016. View at: Google Scholar
19. H. U. Jun, J. Chen, and L. I. Ming, “Research on Drawing Technology of High-strength Steel Rectangular Box Based on ABAQUS,” Hot Working Technology, 2014. View at: Google Scholar
20. S. Liu, H. Shi, and L. Feng, “Transformer harmonic loss model considering skin effect and proximity effect,” Dianli Zidonghua Shebei/Electric Power Automation Equipment, vol. 35, no. 3, pp. 133–139, 2015. View at: Google Scholar
21. Z. Zhao and S.-Y. Zhang, “Self-adaptive step sequential-interaction method and its application in transient-fields coupling problems,” Jisuanji Jicheng Zhizao Xitong/Computer Integrated Manufacturing Systems, CIMS, vol. 17, no. 7, pp. 1404–1414, 2011. View at: Google Scholar
22. Z. Wang, W. Yang, and Y. Song, “The R/S analysis and application in the oil-gas gathering and transportation process,” Petroleum Exploration & Development, vol. 37, no. 5, pp. 618–622, 2010. View at: Google Scholar
23. C. Pandey, A. Giri, and M. M. Mahapatra, “Evolution of phases in P91 steel in various heat treatment conditions and their effect on microstructure stability and mechanical properties,” Materials Science & Engineering A, vol. 664, pp. 58–74, 2016. View at: Google Scholar
24. S. Wang, N. Zhang, L. Wu, and Y. Wang, “Wind speed forecasting based on the hybrid ensemble empirical mode decomposition and GA-BP neural network method,” Journal of Renewable Energy, vol. 94, pp. 629–636, 2016. View at: Publisher Site | Google Scholar
25. Z. Sun, H. Yang, and Z. Tang, “Microstructural evolution model of TA15 titanium alloy based on BP neural network method and application in isothermal deformation,” Computational Materials Science, vol. 50, no. 2, pp. 308–318, 2010. View at: Publisher Site | Google Scholar
26. C. Bharathi, D. Rekha, and V. Vijayakumar, “Genetic Algorithm Based Demand Side Management for Smart Grid,” Wireless Personal Communications, vol. 93, no. 2, pp. 481–502, 2017. View at: Publisher Site | Google Scholar
27. H. S. Madraswala and A. S. Deshpande, “Genetic algorithm solution to unit commitment problem,” in Proceedings of the 2017 International Conference on Nascent Technologies in Engineering, ICNTE 2017, India, January 2017. View at: Google Scholar
28. S. Pal, S. K. Pal, and A. K. Samantaray, “Artificial neural network modeling of weld joint strength prediction of a pulsed metal inert gas welding process using arc signals,” Journal of Materials Processing Technology, vol. 202, no. 1-3, pp. 464–474, 2008. View at: Publisher Site | Google Scholar
29. J. Tang, “Study on temperature control effect with pipe-cooling of thick inverted T-shaped concrete pier-wall during construction,” Water Resources & Hydropower Engineering, 2017. View at: Google Scholar
30. C. Pandey, A. Giri, and M. M. Mahapatra, “On the prediction of effect of direction of welding on bead geometry and residual deformation of double-sided fillet welds,” International Journal of Steel Structures, vol. 16, no. 2, pp. 333–345, 2016. View at: Publisher Site | Google Scholar

Copyright © 2018 Xiurong Fang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Related articles

No related content is available yet for this article.