MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi 10.1155/2018/8430578 8430578 Research Article Parameter Optimization and Prediction Model of Induction Heating for Large-Diameter Pipe http://orcid.org/0000-0003-1096-2791 Fang Xiurong 1 Lu Jia 1 Wang Junfeng 1 Yang Jinhui 1 Afzal Arshad School of Mechanical Engineering Xi’an University of Science and Technology Xi’an 710054 China xust.edu.cn 2018 26112018 2018 26 07 2018 06 10 2018 01 11 2018 26112018 2018 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.

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.

National Natural Science Foundation of China 51775427
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 , and the specific induction heating diagram is shown in Figure 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 . 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 . 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 . 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.

Structure of the bend of induction heating pipe bending machine.

Specific induction heating.

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 .

Schematic diagram of pipe induction heating.

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 .(1)e=-NdΦdt(2)Q=0.24i2Rtwhere N 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 t is the time of the current flowing through the conductor (unit is s).

Taking the elbow forming of X80 pipeline steel  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 (A-Φ) theory of the induction heating electromagnetic field :(3)×1μ×A-1μ·A+σAt+σΦ=0(4)2Φ+t·A=0(5)Je=-σAt-σΦwhere is Hamilton operator; μ is magnetic conductivity (H/m); σ is material conductivity (S/m); ε is capacitivity (F/m); A is magnetic vector potential; Φ is scalar potential; and Je is induction vortex current density vector (A/m2).

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 J. The expression is shown as(6)n×1μ×A=-J

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 , the differential form of the transient temperature field governing equation in induction heating can be derived:(7)ρcTt-xλxTx-yλyTy-zλzTz-qv=0where λ is thermal conductivity (W/(m·°C)); T is temperature field distribution function; qv is intensity of the heat source of induction vortex current in the pipeline (W/m3); ρ is material density (Kg/m3); c is the specific heat capacity of material (J/(Kg· K)).

Use the Newton convection boundary as a temperature boundary condition:(8)-λTn=hT-T0where n 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 T0 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 . 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 . 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.

Analysis process of magnetothermal coupling during the induction heating.

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 1.75×10-8 Ω•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.

The thermal properties of X80.

Temperature(°C) Relative Permeability μ Resistivity ρ (10−6Ω·m) Specific Heat Cp (J/Kg·°C) Thermal Conductivity K (W/m·°C) Density ρ(Kg/m3) Emissivity Cs(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.

Numerical model of induction heating for thin-walled pipe. (a) Geometric model. (b) Finite element model. (c) Finite element model of the pipe.

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 , 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.

Distribution cloud of Joule heat in pipe (Unit: J).

Distribution cloud of temperature in pipe (Unit:°C).

Temperature curves of the inner and outer wall of pipe with the times.

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.

Relationship between three parameters and the inner and outer wall temperature of pipe. (a) Current frequency. (b) Current density. (c) Air gap between pipe and induction coil.

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 , and each factor is set at 5 levels, as shown in Table 2.

Test parameters and levels.

Level Factor
Current Frequency f(Hz) Current Density J(A/m2) Air Gap d(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).

Orthogonal test design scheme and numerical calculation results.

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.

Variance analysis results of the orthogonal test of temperatures for pipe’s outer wall.

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.

Variance analysis results of orthogonal test of temperatures for pipe’s inner wall.

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 . 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:

Range analysis results of orthogonal test for pipe’s outer wall temperature.

Outer Wall Temperature (°C) Inner Wall Temperature (°C)
Current Frequency Current Density Air Gap Current Frequency Current Density Air Gap
K 1 568.6 557.4 1180.6 118 117.6 184.8
K 2 747.2 642.4 1021.2 137.2 129.4 170.6
K 3 867.4 840.4 835.2 151 150 151
K 4 976.2 1001.2 670 165.8 167.8 132.4
K 5 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  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 T1 and T2 of the pipe are taken as output to build a three-layer BP neural network model. The network structure is shown in Figure 10.

Structure of pipe induction heating BP neural network.

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.

BP neural network training error curve.

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.

Error between numerical simulation and BP neural network training.

Current Frequency f(Hz) Current Density J(A/m2) Air Gap d(mm) Outer Wall Temperature T1(°C) Inner Wall Temperature T2(°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.

Genetic algorithm optimization process.

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 T1 of the outer wall and the temperature T2 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 T 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 .

Regression Model I: T=b0ƒ+b1J+b2d+C

Regression Model II: T=Cƒb0Jb1db2

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

(1) εi  ~N(0,σ2), Cov(εi,εj)=0, ij

(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 T 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.

Results of regression analysis.

Model I: T=b0f+b1J+b2d+C Model II: T=Cfb0Jb1db2
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: T=0.525f+3.092  ×  10-5J-17.184d-1049.604

Equation of Model II: T=e-33.103ƒ0.642J2.061d-0.481

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.

Errors between the numerical simulation and regression model prediction.

Current Frequency f(Hz) Current Density J(A/m2) Air Gap d(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%

Regression standardized residual histogram.

Model I

Model II

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)

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