Next Article in Journal
Control of Laves Precipitation in a FeCrAl-based Alloy Through Severe Thermomechanical Processing
Previous Article in Journal
Electronic Transport Properties of Silicane Determined from First Principles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Simplified Calculation Method of Heat Source Model for Induction Heating

1
School of Naval Architecture and Ocean Engineering, Huazhong University of Science and Technology, Wuhan 430074, China
2
Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration (CISSE), Shanghai 200240, China
3
Hubei Key Laboratory of Naval Architecture and Ocean Engineering Hydrodynamics (HUST), Wuhan 430074, China
4
Shanghai Waigaoqiao Shipbuilding Co., Ltd., Shanghai 200137, China
*
Author to whom correspondence should be addressed.
Materials 2019, 12(18), 2938; https://doi.org/10.3390/ma12182938
Submission received: 4 August 2019 / Revised: 4 September 2019 / Accepted: 9 September 2019 / Published: 11 September 2019

Abstract

:
Line heating is used in forming the complex curve plates of ships, and this process is becoming integrated into automated tools. Induction heating equipment has become commonly used in automatic line heating. When applying automated equipment, it is necessary to calculate the relationship between the heating parameters and the temperature field. Numerical methods are primarily used to accomplish the calculations for induction heating. This computation process requires repeated iterations to obtain a stable heat generation rate. Once the heat generation rate changes significantly, a recalculation takes place. Due to the relative position of the coil and plate changes during heating, the grid needs to be frequently re-divided during computation, which dramatically increases the total computation time. In this paper, through an analysis of the computation process for induction heating, the root node that restricts the computation efficiency in the conventional electromagnetic-thermal computation process was found. A method that uses a Gaussian function to represent the heat flux was proposed to replace the electromagnetic computation. The heat flux is the input for calculating the temperature field, thus avoiding the calculation of the electromagnetic analysis during induction heating. Besides, an equivalence relationship for multi-coil was proposed in this paper. By comparing the results of the experiment and the numerical method, the proposed heat source model’s effectiveness was verified.

1. Introduction

Line heating is a typical process in the plate forming of ships. In the past, the operation was mostly manual and included the use of a flame heater. With advancements in research and the development of automated line heating equipment, flame heating by oxyacetylene combustion has shown its limitations because of the difficulty with temperature control. Induction heating is used extensively for the forming of plates because it is easily controlled and has high efficiency.
The principle of electromagnetic heating is to generate an eddy current in a workpiece, which is then heated by the Joule heat generated by the workpiece’s resistance. An induction heater includes a power source and an induction coil. The shape of the induction coil varies based on the object to be heated. The induction heating process includes electromagnetic induction and heat conduction. It is tough to obtain the temperature distribution of the workpiece during induction heating through a theoretical method and experiment-based methods are too costly. Therefore, numerical simulation has become an effective solution for obtaining the temperature fields of induction heating. A typical process for induction heating simulation is presented in Figure 1. First, time-harmonic electromagnetic analysis is performed with the designated initial conditions and constraint conditions under a given frequency and power. The Joule heat distribution of the induction coil is calculated when a stable heat generation rate is obtained. The temperature field of the workpiece can be calculated using the Joule heat distribution as the input of thermal analysis. The thermal analysis ends when a stable node temperature is obtained. If the set calculation time has not been reached, the induction coil is then moved to the next position. A new analysis is carried out when establishing the constraint equations about the coil-workpiece-air interface.
In the computation process mentioned above, it is necessary to calculate the electromagnetic field and the temperature field. The Joule heat obtained from the electromagnetic calculation is the input of the heat transfer calculation. It is necessary to frequently re-divide the grid based on the coil position when performing moving induction heating forming simulations, which increase the time of computations and decrease computation efficiency.
Some studies have investigated methods to improve computation efficiency. Monzel [1] used the neural network method to carry out preliminary research on the features of the inductive current distribution during transverse flux induction heating. He gave measures for reducing cable loss by changing each coil’s current to change the heat source distribution. The typical form of heat source density was given in the study, but the corresponding simplified numerical method was not mentioned. Bay [2] postulated that the temperature calculation relative to the electromagnetic calculation was a quasi-static process and proposed mathematical and numerical models for the coupled axisymmetric induction heating process. The model simplified the full-coupling process into a weak coupling process, used the eddy current obtained in the electromagnetic calculation as the heat generation rate in the temperature calculation, and presented a preliminary calculation method for the heat generation rate, but the expression for the heat source model was not given. Luo [3,4,5,6] carried out finite element calculations on the temperature field and the deformation field during high-frequency induction heating. He compared the relevant results with the experimental results, which avoided the electromagnetic coupling calculation. He gave the corresponding heat source model, but the calculation method for the heat source model and the corresponding relationship of various parameters in the model with the actual coil were not discussed. Liu [7] presented a numerical method for the temperature field of induction heating and studied the law of temperature distribution. When handling the electromagnetic-thermal coupling problem, a heat source model based on an empirical formula was used, but the calculation process for the model was not given. Hu [8] proposed the use of an equivalent heat source to replace the electromagnetic coupling calculation in induction heating and carried out experimental verification of the proposed heat source model, but did not give an explicit expression for the heat source model. Zhang [9] analyzed the similarities and differences of induction heating and flame heating in the line heating of ship plates and thought that induction heating was feasible in heat forming. When handling the induction heating process, he provided a heat source model for induction heating and used it as the initial input for the thermal distribution calculation. Compared with the experimental results, the errors of the calculation results satisfied practical engineering requirements. However, the model was simplified and the meaning of various model parameters and the corresponding numerical methods were not explained. Bae [10] proposed a two-dimensional circular heat input model to simulate the induction heating process and obtained satisfactory results, but the expression of the heat source model was not taken into account. Jeong [11] simplified coupled induction heating into the electromagnetic calculations and thermal calculations. He used statistical methods to describe the correlation between deformation and the input parameters, but the expression of the heat source model was not described in the study. Bai [12] used stepwise analysis to carry out coupling analysis within a typical time to obtain the state of induction heat distribution. He used the distribution as the moving heat source for a subsequent calculation, but the expression for the heat source was not given. Kubota [13] used a three-dimensional coupling and heat source model to perform calculations on the induction heating process and found that the results were consistent with measured results. The heat source model was obtained based on the obtained temperature field. Yang [14] proposed a heat source model based on the characteristics of the induction heating eddy current and used the finite element method (FEM) and relevant tests to verify the effectiveness of the model. The proposed heat source model was based on empirical methods and only targeted a particular coil form without giving the expression form of the corresponding heat source for other coil forms. Through theoretical analysis, Li [15] studied the analytical solution for obtaining the heat generation rate in the semi-infinite space during induction heating and discussed the effects of relevant parameters on the heat generation rate. To compare the difference in residual stress in a plate after single heating and double heating, Aung [16] proposed numerical methods based on a surface heat source and a body heat source. The effectiveness of the heat source models was verified through the experiment. Aung provided the expression form for the heat source model but did not discuss its calculation. Riccio [17] used different numerical models to study the bonding of carbon fiber reinforced polymer components with induction heating. ABAQUS was used to carry out electromagnetic analysis to obtain the energy loss caused by the Joule effect. He provided the finite element calculation model but not the expression for the Joule heat of the coil. In other studies related to the induction heating process and induction coil design [18,19,20], the full-coupling numerical method was also used for the induction heating calculations. Although the results were satisfactory, the time cost was high and the applicability was low for coils of different shapes and sizes or with different processing conditions. Zhang [21] replaced the electromagnetic-thermal coupling calculations with a heat source model corresponding to a discrete form and proposed a simplified calculation method for the high-frequency induction heat source. The computation efficiency for the calculation of a certain specific induction heating parameter was greatly improved, but any change in the induction heating parameters necessitated repeating the entire calculation.
The induction heating calculation, especially the direct calculation of the moving induction heating process, results in a long computation time due to the need to continuously change the grid. However, the use of the heat flux to replace the coil can achieve high precision and reduce the computation time. Currently, most studies focus on empirical methods, and there are comparatively few studies that calculate the heat source and investigate the effect of the coil on the model. Based on these, the analytical method was employed in this study to calculate the heat source model for the induction coil. The effects of coil shapes on the heat source were examined. In addition, related experiments were conducted to verify the obtained heat source model for single-coil and multi-coil.

2. Simplification of Heat Source

2.1. Finite Element Calculation

In the induction heating process, the input is alternating current. The alternating current generates an alternating magnetic field in the coil which generates an induced electric field in the workpiece, which in turn generates heat in the workpiece. The electromagnetic field in the induction heating process is given by Maxwell’s equations:
× H = J + D t
× E = B t
· B = 0
· D = ρ
where H is the magnetic field strength, J is the conduction current density, D t is the displacement current density, t represents time, E is the electric field, B is the magnetic field, D is the electric flux density, and ρ is the volume charge density.
The corresponding auxiliary equations are:
D = ε E
B = μ H
J = σ E
where ε represents the dielectric constant, μ is the magnetic permeability, and σ is the electrical conductivity.
The eddy current density distribution in the workpiece can be obtained by solving the aforementioned equations and the distribution of the heat generation rate can be obtained using Joule’s law:
q = | J | 2 σ
where q represents the heat generation rate.
The heat generation rate is used as the input in the thermal conduction equation for calculating the thermal distribution, and the temperature field of the plate is obtained:
ρ c T t = x ( λ T x ) + y ( λ T y ) + z ( λ T z ) + q
where ρ is the density, c is the specific heat, T represents the temperature, and λ represents the thermal conductivity.
The above procedure is the detailed process for the calculation of the induction heating. It is necessary to continuously update the corresponding parameters with temperature changes during the calculation as the workpiece’s electromagnetic and thermal property parameters are related to the temperature, which leads to an extremely time-consuming computation process. When the coil and the workpiece have relative movement, the computation process, as mentioned earlier, must be repeated to obtain the temperature field during movement.

2.2. Simplification of Heat Source Model

The key to temperature calculation is in obtaining the distribution of the heat generation rate   q . Due to the skin effect, the heat generated by the induction is concentrated on the surface of the workpiece so the surface heat flux can be used as the input during induction heating. The heat flux model used in this paper is shown in Figure 2. The cumbersome electromagnetic calculation is unnecessary once the heat flux model is used to represent the coil.
The heat flux along the surface of the workpiece in the radial direction of a single coil was assumed to follow a Gaussian function:
q ( r ) = q m e ( r R 0 r H ) 2
where q ( r ) represents the heat flux at a distance r from the center of the heat source, r is the coordinate along the radial direction of the coil, q m is the maximum heat flux, R 0 represents the radius of the heat source, and r H represents the effective radius of the heat flux.
The total power on the workpiece is equal to the effective power Q of the induction coil, that is:
Q = 0 2 π R 0 d θ + q ( r ) d r = 2 π π q m R 0 r H
which leads to:
q m = Q 2 π π R 0 r H
Thus:
q ( r ) = Q 2 π π R 0 r H e ( r R 0 r H ) 2
Equation (13) is the heat flux corresponding to the coil. The temperature field of a workpiece can be directly calculated after the heat flux is substituted into Equation (9). Since the Joule heat distribution is obtained directly from Equation (13) to represent the induction coil, the electromagnetic calculation in Figure 1 is unnecessary, and the heat transfer is directly calculated. Moreover, grid redivision during coil movement is unnecessary, resulting in fast computation.

3. Verifications of Heat Source Model

3.1. Experimental Verification for Single Coil

The heat flux model of a coil during induction heating has been previously stated, and a static induction heating experiment was conducted on a three-axis motion platform to verify the effectiveness of the model. The platform had degrees of freedom in the three directions of X-Y-Z with a positional precision of 1 mm and a velocity precision of 1 mm/s. The schematic diagram of the induction heating experiment is shown in Figure 3.
The experimental plate was low carbon steel and the material properties were similar to those in the Reference [22], as listed in Table 1 and Table 2. The size of the plate was 1000 mm × 800 mm × 16 mm.
A ring-type coil was used in the experiment, as shown in Figure 3. The induction coil was copper with an inner diameter of 90 mm and an outer diameter of 110 mm. The cross-section of the coil was a rectangle with a size of 10 mm × 8 mm and a wall thickness of 1.5 mm. The schematic diagram of the single coil is illustrated in Figure 4.
The frequency of induction heating during the experiment was 12 kHz and the input power was 7.5 kW. The distance between the coil and the surface of the plate was 5 mm, which remained unchanged during heating. The heating time was 120 s. The induction heating efficiency was assumed to be 80%. An infrared temperature sensor (with a measurement range of 250 °C–1450 °C) was used to measure the temperature history of the plate. The distances between the sensors and the center of the coil were 0 mm, 30 mm, and 70 mm. A schematic diagram of the temperature measurement points is shown in Figure 3b.

3.2. FEM Calculation and Comparison

Following experimental conditions, the heat flux function of the coil during heating can be calculated based on the input parameters from Equation (13):
q ( r ) = 5.39 × 10 5 × e ( r 0.05 0.02 ) 2
The heat flux obtained was used as the input of the heat transfer calculation to obtain the temperature field of the plate. ABAQUS was used to carry out the heat transfer calculation. The plate size was 1000 mm × 800 mm × 16 mm. The element type of the plate is DC3D8 for heat transfer simulation. The finite element model under different grid sizes was established to investigate the effect of the element size. The correlation between the maximum temperature and calculation time and the number of grids in the heat affected zone (HAZ) was compared. The width of the HAZ was 400 mm. The grids of the length direction and the width direction were the same sizes, and the thickness direction was unified into four grids. All calculations were performed on the same computer: Windows 10, 64 bits, four cores with 4.2 GHz and 8 GB of RAM. The relationship between the maximum temperature of the plate and calculation time and the number of elements in HAZ is given in Figure 5.
As seen from Figure 5 that the maximum temperature gradually decreases and reaches a stable value as the number of grids increases to 80 (which means the element size is 5 mm). However, the computing time increases sharply when the number of grids is over 100 (which means the element size is 4 mm). We can conclude that the usage of an element of size equal to 5 mm is appropriate in reasonable computational time without losing appreciable accuracy. A variable grid was used for element division. The grids in the HAZ were 5 mm × 4 mm × 4 mm and the elements far from the heating zone were 50 mm × 20 mm × 8 mm. The total number of elements was 68,600 and the number of grids was 86,399. The finite element result is given in Figure 6.
Comparisons of the experiment results and finite element calculation using the heat source model are shown in Figure 7, Figure 8 and Figure 9, the relative locations of C1, C2, and C3 are shown in Figure 3.
It can be seen from Figure 7, Figure 8 and Figure 9 that the temperature histories obtained using the heat source model for calculation are in agreement with the experiment results. Table 3 summarized the relative error between simulated results and experimental results. It can be seen that the maximum relative error was 5.42% and the average relative error was 3.64%. The comparison indicates that the temperature field of induction heating can be accurately simulated using a heat source model.
The computation efficiency using the heat source model was compared with that using the full-coupling method. COMSOL Multiphysics was used to carry out the coupling calculations. The maximum element size of the model was 1 mm. The total number of elements was 34,274. The coupling model is shown in Figure 10.
The comparisons between the coupling model and the experiment are shown in Figure 11, Figure 12 and Figure 13. The comparison of relative error between coupling model and the experiment are listed in Table 4.
As can be seen from Figure 11, Figure 12 and Figure 13, the temperature histories obtained by the coupling model are in agreement with the experiment results. The relative error between the simulated results and experimental results, listed in Table 4, showed that the average relative error was 7.55% at most, which means the chosen simulation method was correct for calculating induction heating.
Table 5 shows the number of degrees and computing time between the coupling model and the heat source model. The number of degrees of freedom and the computation time during coupling computing was higher than the heat source model. It can be considered that the heat source model has sufficient accuracy and efficiency for induction heating calculation.

3.3. Application to Multi-Coil

It can be seen from Equation (13) that the Gaussian function that represents the heat source has a relationship with the project area of the coil. To further verify the effectiveness of the proposed heat source model, the previously described coil, shown in Figure 3b, was used to conduct a multi-coil experiment.
A static induction heating experiment of multi-coil was conducted on the previously described experimental platform. The schematic diagram of the multi-coil is shown in Figure 14. The heating power was 25 kW and the heating time was 100 s. The bottom of the coil was 5 mm from the surface of the plate. The size of the plate and positions of the temperature measurement points were the same as those in the previously described experiment. The experimental results were compared with the finite element results calculated using the heat source model.
During the experiment, the heat flux was more diffused along the radial direction of the coil and its scope of influence had become larger. Therefore, the action radius of the heat flux was assumed to be   r H = 0.04 , the position of the heat flux center was at R 0 = 0.085 in the middle coil, and the heat flux was given by:
q ( r ) = 5.28 × 10 5 × e ( r 0.085 0.04 ) 2
Comparisons of the simulation results and the experiment results are shown in Figure 15 and Figure 16. The error comparison between simulation results and experimental results are listed in Table 6.
It can be seen from Figure 15 and Figure 16 and Table 6 that the results of the heat source model are in agreement with the results obtained in the experiment in terms of the trend with a maximum relative error of 18.6%. It means that the proposed heat source model is also applicable for use with multi-coil induction heating, which expands the application of the heat source model.

4. Conclusions

In this paper, the heat flux of the coil during induction heating was studied. The expression of the heat flux was obtained and the accuracy was verified through calculation and experiment. The following conclusions were drawn:
1. The heat flux during induction heating can be simplified into a Gaussian function, and the values of the Gaussian function are related to the size of the induction coil and the heating conditions.
2. The heat source model can be used to calculate the temperature distribution of induction heating. The proposed model can improve computation efficiency while ensuring precision. Fast computation of moving induction heating can be achieved using this model.
3. The proposed heat source can also be used in multi-coil induction heating, as the goodness of the finite element results and the experimental result is within an acceptable range.

Author Contributions

Y.Z. contributed to the discussions and reviewed the paper; H.D. carried out the experiments, performed the numerical simulations and wrote the paper; H.Y. assisted in the experiment; X.H. guided the experiment; Z.Y. reviewed the paper.

Funding

This work was supported by the Ministry of Science and Technology of the People’s Republic of China (Grant No.: 2012DFR80390) and Ministry of Industry and Information Technology of the People’s Republic of China (Grant No.:2016545).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Monzel, C.; Henneberger, G. Optimizing heat source distribution for transverse flux inductive heating devices for thin strips with genetic algorithms. In Proceedings of the Power Electronics and Variable Speed Drives, London, UK, 18–19 September 2000; pp. 426–430. [Google Scholar]
  2. Bay, F.; Labbe, V.; Favennec, Y.; Chenot, J.L. A numerical model for induction heating processes coupling electromagnetism and thermomechanics. Int. J. Numer. Methods Eng. 2003, 58, 839–867. [Google Scholar] [CrossRef]
  3. Shuai, K.; Luo, Y.; Sha, W.; Xie, L. FEM Analysis of the Temperature Field of High Frequency Induction Heating in Plate Bend Molding. Boil. Technol. 2004, 35, 52–54. [Google Scholar]
  4. Luo, Y.; Ishiyama, M.; Murakawa, H. Study of Temperature Field and Inherent Strain Produced by High Frequency Induction Heating on Flat Plate. Trans. JWRI 2004, 33, 59–63. [Google Scholar]
  5. Zhou, H.; Jiang, Z.; Luo, Y.; Luo, P. Research on Hull Plate Bending Based on the Technology of High-Frequency Induction. Shipbuild. China 2014, 55, 128–138. [Google Scholar]
  6. Zhou, H.; Jiang, Z.; Luo, Y.; Luo, P. Influence of Thermal Properties of Plate on Its High-Frequency Induction Bending in Forming Ship Hulls. Shipbuild. China 2015, 56, 101–108. [Google Scholar]
  7. Liu, F.; Shi, Y.; Lei, X. Numerical investigation of the temperature field of a metal plate during high-frequency induction heat forming. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2008, 223, 979–986. [Google Scholar] [CrossRef]
  8. Hu, X. A Study on Heating Source Model and Deformation Law for Line Heating by High Frequency Induction Based on Numerical Experiment. Master’s Thesis, Dalian University of Technology, Dalian, China, November 2008. [Google Scholar]
  9. Zhang, X.; Yang, Y.; Liu, Y. Feasibility research on application of a high frequency induction heat to line heating technology. J. Mar. Sci. Appl. 2011, 10, 456–464. [Google Scholar] [CrossRef]
  10. Bae, K.Y.; Yang, Y.S.; Hyun, C.M. Analysis of Triangle Heating Technique using High Frequency Induction Heating in Forming Process of Steel Plate. Int. J. Precis. Eng. Manuf. 2012, 13, 539–545. [Google Scholar] [CrossRef]
  11. Jeong, C.M.; Yang, Y.S.; Bae, K.Y.; Hyun, C.M. Prediction of deformation of steel plate with forced displacement and initial curvature in a forming process with high frequency induction heating. Int. J. Precis. Eng. Manuf. 2013, 14, 785–790. [Google Scholar] [CrossRef]
  12. Bai, X.; Zhang, H.; Wang, G. Modeling of the moving induction heating used as secondary heat source in weld-based additive manufacturing. Int. J. Adv. Manuf. Technol. 2014, 77, 717–727. [Google Scholar] [CrossRef]
  13. Kubota, H.; Tomizawa, A.; Yamamoto, K.; Okada, N.; Hama, T.; Takuda, H. Finite Element Analysis of Three-Dimensional Hot Bending and Direct Quench Process Considering Phase Transformation and Temperature Distribution by Induction Heating. ISIJ Int. 2014, 54, 1856–1865. [Google Scholar] [CrossRef] [Green Version]
  14. Yang, S. Research on Heat Source Model and Numerical Simulation of High Frequency Induction Heating Forming. Dev. Appl. Mater. 2015, 30, 20–23. [Google Scholar]
  15. Li, J. Numerical Studies on the Induction Quenching Process of Crankshaft. Doctor’s Thesis, Beijing Institute of Technology, Beijing, China, January 2015. [Google Scholar]
  16. Aung, M.P.; Nakamura, M.; Hirohata, M. Characteristics of Residual Stresses Generated by Induction Heating on Steel Plates. Metals 2018, 8, 25. [Google Scholar] [CrossRef]
  17. Riccio, A.; Russo, A.; Raimondo, A.; Cirillo, P.; Caraviello, A. A numerical/experimental study on the induction heating of adhesives for composite materials bonding. Mater. Today Commun. 2018, 15, 203–213. [Google Scholar] [CrossRef]
  18. Bao, L.; Qi, X.W.; Mei, R.B.; Zhang, X.; Li, G.L. Investigation and modeling of work roll temperature in induction heating by finite element method. J. S. Afr. Inst. Min. Met. 2018, 118, 735–743. [Google Scholar]
  19. Li, F.; Li, X.K.; Qin, X.F.; Rong, Y. Study on the plane induction heating process strengthened by magnetic flux concentrator based on response surface methodology. J. Mech. Sci. Technol. 2018, 32, 2347–2356. [Google Scholar] [CrossRef]
  20. Sun, R.; Shi, Y.J.; Pei, Z.F.; Li, Q.; Wang, R.H. Heat transfer and temperature distribution during high-frequency induction cladding of 45 steel plate. Appl. Therm. Eng. 2018, 139, 1–10. [Google Scholar] [CrossRef]
  21. Zhang, Z.; Zhao, Y.; Hu, X.; Yang, Z. Simplified calculation method for high frequency induction heating source of ship plate. Chin. J. Ship Res. 2018, 5, 18–24. [Google Scholar]
  22. Li, Y. Steel Structure Welding Residual Stress Analysis. Master’s Thesis, Wuhan University of Technology, Wuhan, China, April 2007. [Google Scholar]
Figure 1. Electromagnetic-thermal computation process.
Figure 1. Electromagnetic-thermal computation process.
Materials 12 02938 g001
Figure 2. Diagrams of computation model and heat source distribution: (a) Computation model and (b) Gaussian distribution of the heat source.
Figure 2. Diagrams of computation model and heat source distribution: (a) Computation model and (b) Gaussian distribution of the heat source.
Materials 12 02938 g002
Figure 3. Schematic diagram of the induction heating experiment: (a) Experimental platform for induction heating and (b) Positions of the temperature measurement points.
Figure 3. Schematic diagram of the induction heating experiment: (a) Experimental platform for induction heating and (b) Positions of the temperature measurement points.
Materials 12 02938 g003
Figure 4. Schematic diagram of the single coil.
Figure 4. Schematic diagram of the single coil.
Materials 12 02938 g004
Figure 5. Relationship between maximum temperature and computational time and number of elements in the heat affected zone (HAZ).
Figure 5. Relationship between maximum temperature and computational time and number of elements in the heat affected zone (HAZ).
Materials 12 02938 g005
Figure 6. The finite element result of the heat transfer calculation.
Figure 6. The finite element result of the heat transfer calculation.
Materials 12 02938 g006
Figure 7. Comparison of the temperature histories of C1.
Figure 7. Comparison of the temperature histories of C1.
Materials 12 02938 g007
Figure 8. Comparison of the temperature histories of C2.
Figure 8. Comparison of the temperature histories of C2.
Materials 12 02938 g008
Figure 9. Comparison of the temperature histories of C3.
Figure 9. Comparison of the temperature histories of C3.
Materials 12 02938 g009
Figure 10. Coupling calculation model.
Figure 10. Coupling calculation model.
Materials 12 02938 g010
Figure 11. Comparison of the temperature histories of C1 between the coupling model and the experiment.
Figure 11. Comparison of the temperature histories of C1 between the coupling model and the experiment.
Materials 12 02938 g011
Figure 12. Comparison of the temperature histories of C2 between the coupling model and the experiment.
Figure 12. Comparison of the temperature histories of C2 between the coupling model and the experiment.
Materials 12 02938 g012
Figure 13. Comparison of the temperature histories of C3 between the coupling model and the experiment.
Figure 13. Comparison of the temperature histories of C3 between the coupling model and the experiment.
Materials 12 02938 g013
Figure 14. Schematic diagram of sizes of multiple coils.
Figure 14. Schematic diagram of sizes of multiple coils.
Materials 12 02938 g014
Figure 15. The temperature histories of C3 during the multi-coil experiment and simulation.
Figure 15. The temperature histories of C3 during the multi-coil experiment and simulation.
Materials 12 02938 g015
Figure 16. The temperature histories of C2 during the multi-coil experiment and simulation.
Figure 16. The temperature histories of C2 during the multi-coil experiment and simulation.
Materials 12 02938 g016
Table 1. Thermal parameters of low carbon steel.
Table 1. Thermal parameters of low carbon steel.
Temperature (°C)Density (kg/m3)Specific Heat (J/(kg·°C))Heat Conductivity Coefficient (W/(m·°C))
07842450.3666.97
50-464.665.21
2007822498.157.38
250-502.2654.91
300-514.8253
4007802537.4247.92
450-623.6445.83
500-707.3543.53
600778281239.3
650-904.0736.37
700-967.6934.74
80077611026.3231.02
Table 2. Mechanical parameters of low carbon steel.
Table 2. Mechanical parameters of low carbon steel.
Temperature (°C)Young‘s Modulus (GPa)Poisson RatioHeat Expansion Coefficient (1/°C)Yield Strength (MPa)
02060.2671.20 × 10−5235
501960.291.25 × 10−5-
2001960.3221.40 × 10−5163
2501860.2961.43 × 10−5-
3001860.2621.47 × 10−5-
4001660.241.54 × 10−5130
4501570.2291.57 × 10−5-
5001570.2231.59 × 10−5-
6001350.2231.64 × 10−5119
6501170.2231.66 × 10−5-
7001120.2231.67 × 10−5-
8001130.2231.69 × 10−5109
Table 3. Comparison between simulated results and experimental results.
Table 3. Comparison between simulated results and experimental results.
SensorMaximum Relative ErrorAverage Relative ErrorMinimum Relative Error
C13.55%2.43%1.15%
C23.84%2.42%1.16%
C35.42%3.64%2.46%
Table 4. Comparison of coupling model and experiment.
Table 4. Comparison of coupling model and experiment.
SensorMaximum Relative ErrorAverage Relative ErrorMinimum Relative Error
C18.65%7.55%6.38%
C22.71%1.68%0.50%
C31.61%1.17%0.27%
Table 5. Comparison of computation efficiency.
Table 5. Comparison of computation efficiency.
Number of Degrees of FreedomComputation Time/Min
Full-coupling model230,9176.7
Heat source model86,3993.7
Table 6. Comparison of simulation and experiment.
Table 6. Comparison of simulation and experiment.
SensorMaximum Relative ErrorAverage Relative ErrorMinimum Relative Error
C218.63%12.72%4.82%
C315.42%7.65%0.15%

Share and Cite

MDPI and ACS Style

Dong, H.; Zhao, Y.; Yuan, H.; Hu, X.; Yang, Z. A Simplified Calculation Method of Heat Source Model for Induction Heating. Materials 2019, 12, 2938. https://doi.org/10.3390/ma12182938

AMA Style

Dong H, Zhao Y, Yuan H, Hu X, Yang Z. A Simplified Calculation Method of Heat Source Model for Induction Heating. Materials. 2019; 12(18):2938. https://doi.org/10.3390/ma12182938

Chicago/Turabian Style

Dong, Hongbao, Yao Zhao, Hua Yuan, Xiaocai Hu, and Zhen Yang. 2019. "A Simplified Calculation Method of Heat Source Model for Induction Heating" Materials 12, no. 18: 2938. https://doi.org/10.3390/ma12182938

APA Style

Dong, H., Zhao, Y., Yuan, H., Hu, X., & Yang, Z. (2019). A Simplified Calculation Method of Heat Source Model for Induction Heating. Materials, 12(18), 2938. https://doi.org/10.3390/ma12182938

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop