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

Authors: Xiurong Fang; Jia Lu; Junfeng Wang; Jinhui Yang
Journal: Mathematical Problems in Engineering
 (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 ofX80 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.


---

## Body

## 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 Figure1 [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 [4–7]. 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.Figure 1
Structure of the bend of induction heating pipe bending machine.Figure 2
Specific induction heating.


## 2. Analysis on Induction Heating Parameters of Pipes

The basic principle of induction heating is shown in Figure3: 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].Figure 3
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 [10].(1)e=-NdΦdt(2)Q=0.24i2Rtwhere N is the number of turns of the induction coil;Q is Joule heat (unit isJ);i is the effective value of induction current (unit isA);Ris conductor resistance (unit isΩ); and t is the time of the current flowing through the conductor (unit iss).Taking the elbow forming ofX80 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 (A→-Φ) theory of the induction heating electromagnetic field [12–14]:(3)∇×1μ∇×A→-∇1μ∇·A→+σ∂A→∂t+σ∇Φ=0(4)∇2Φ+∂∂t∇·A→=0(5)Je→=-σ∂A→∂t-σ∇Φ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 densityJ→. 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 [17], the differential form of the transient temperature field governing equation in induction heating can be derived:(7)ρc∂T∂t-∂∂xλx∂T∂x-∂∂yλy∂T∂y-∂∂zλz∂T∂z-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);cis the specific heat capacity of material (J/(Kg· K)).Use the Newton convection boundary as a temperature boundary condition:(8)-λ∂T∂n=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 [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.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 theX80 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.Table 1
The thermal properties ofX80.
Temperature(°C) Relative Permeabilityμ Resistivityρ (10−6Ω·m) Specific HeatCp (J/Kg·°C) Thermal ConductivityK (W/m·°C) Densityρ(Kg/m3) EmissivityCs(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.535It 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 Figure1). 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.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.
(a) (b) (c)


## 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 [12–14]:(3)∇×1μ∇×A→-∇1μ∇·A→+σ∂A→∂t+σ∇Φ=0(4)∇2Φ+∂∂t∇·A→=0(5)Je→=-σ∂A→∂t-σ∇Φ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 densityJ→. 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 [17], the differential form of the transient temperature field governing equation in induction heating can be derived:(7)ρc∂T∂t-∂∂xλx∂T∂x-∂∂yλy∂T∂y-∂∂zλz∂T∂z-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);cis the specific heat capacity of material (J/(Kg· K)).Use the Newton convection boundary as a temperature boundary condition:(8)-λ∂T∂n=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 [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.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 theX80 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.Table 1
The thermal properties ofX80.
Temperature(°C) Relative Permeabilityμ Resistivityρ (10−6Ω·m) Specific HeatCp (J/Kg·°C) Thermal ConductivityK (W/m·°C) Densityρ(Kg/m3) EmissivityCs(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.535It 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 Figure1). 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.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.
(a) (b) (c)


## 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.Figure 6
Distribution cloud of Joule heat in pipe (Unit:J).Figure 7
Distribution cloud of temperature in pipe (Unit:°C).Figure 8
Temperature curves of the inner and outer wall of pipe with the times.It can be seen from Figure8 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 Figure9.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.
(a) (b) (c)Figure9 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.Table 2
Test parameters and levels.
Level Factor Current Frequencyf(Hz) Current DensityJ(A/m2) Air Gapd(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 50It can be seen from Table2 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 ofX80 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).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.Table 4
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.Table 5
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 [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: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 Tables4 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 temperaturesT1 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.Figure 10
Structure of pipe induction heating BP neural network.The data (see Table3) 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.Figure 11
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 Table7.Table 7
Error between numerical simulation and BP neural network training.
Current Frequencyf(Hz) Current DensityJ(A/m2) Air Gapd(mm) Outer Wall TemperatureT1(°C) Inner Wall TemperatureT2(°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 Table7, 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.Figure 12
Genetic algorithm optimization process.It can be seen from the Figure12 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 densityJ(A/m2), and air gapd(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 [30].Regression Model I:T=b0ƒ+b1J+b2d+CRegression Model II:T=Cƒb0Jb1db2Assume that the above two regression models meet the following conditions:(1)εi~N(0,σ2), Cov(εi,εj)=0, i≠j(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 variableT and the transformed independent variable in the regression model IIThe two regression models are analyzed and solved by a multiple regression calculation program, and the results are shown in Table8.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−19The regression equations obtained by stepwise regression calculation are as follows:Equation of Model I:T=0.525f+3.092×10-5J-17.184d-1049.604Equation of Model II:T=e-33.103ƒ0.642J2.061d-0.481It can be seen from Figure13 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.Table 9
Errors between the numerical simulation and regression model prediction.
Current Frequencyf(Hz) Current DensityJ(A/m2) Air Gapd(mm) Outer Wall TemperatureT (°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%Figure 13
Regression standardized residual histogram.
(a)
Model I (b)
Model IIIt can be seen from Table9 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.


## 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.Figure 6
Distribution cloud of Joule heat in pipe (Unit:J).Figure 7
Distribution cloud of temperature in pipe (Unit:°C).Figure 8
Temperature curves of the inner and outer wall of pipe with the times.It can be seen from Figure8 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 Figure9.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.
(a) (b) (c)Figure9 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.Table 2
Test parameters and levels.
Level Factor Current Frequencyf(Hz) Current DensityJ(A/m2) Air Gapd(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 50It can be seen from Table2 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 ofX80 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).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.Table 4
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.Table 5
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 [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: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 Tables4 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.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.Table 2
Test parameters and levels.
Level Factor Current Frequencyf(Hz) Current DensityJ(A/m2) Air Gapd(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 50It can be seen from Table2 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 ofX80 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).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.Table 4
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.Table 5
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 [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: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 Tables4 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 temperaturesT1 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.Figure 10
Structure of pipe induction heating BP neural network.The data (see Table3) 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.Figure 11
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 Table7.Table 7
Error between numerical simulation and BP neural network training.
Current Frequencyf(Hz) Current DensityJ(A/m2) Air Gapd(mm) Outer Wall TemperatureT1(°C) Inner Wall TemperatureT2(°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 Table7, 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.Figure 12
Genetic algorithm optimization process.It can be seen from the Figure12 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.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 temperaturesT1 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.Figure 10
Structure of pipe induction heating BP neural network.The data (see Table3) 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.Figure 11
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 Table7.Table 7
Error between numerical simulation and BP neural network training.
Current Frequencyf(Hz) Current DensityJ(A/m2) Air Gapd(mm) Outer Wall TemperatureT1(°C) Inner Wall TemperatureT2(°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 Table7, 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.Figure 12
Genetic algorithm optimization process.It can be seen from the Figure12 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 densityJ(A/m2), and air gapd(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 [30].Regression Model I:T=b0ƒ+b1J+b2d+CRegression Model II:T=Cƒb0Jb1db2Assume that the above two regression models meet the following conditions:(1)εi~N(0,σ2), Cov(εi,εj)=0, i≠j(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 variableT and the transformed independent variable in the regression model IIThe two regression models are analyzed and solved by a multiple regression calculation program, and the results are shown in Table8.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−19The regression equations obtained by stepwise regression calculation are as follows:Equation of Model I:T=0.525f+3.092×10-5J-17.184d-1049.604Equation of Model II:T=e-33.103ƒ0.642J2.061d-0.481It can be seen from Figure13 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.Table 9
Errors between the numerical simulation and regression model prediction.
Current Frequencyf(Hz) Current DensityJ(A/m2) Air Gapd(mm) Outer Wall TemperatureT (°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%Figure 13
Regression standardized residual histogram.
(a)
Model I (b)
Model IIIt can be seen from Table9 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-diameterX80 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 ofX80 pipe, the numerical simulation process parameters sets are somewhat limited. Hopefully, our mathematical analysis method is of a referential significance for practical engineering.



---

*Source: 8430578-2018-11-26.xml*