CN107742030B - Simulation method for intermediate frequency heating and pulse current application of TP2 internal thread copper pipe - Google Patents
Simulation method for intermediate frequency heating and pulse current application of TP2 internal thread copper pipe Download PDFInfo
- Publication number
- CN107742030B CN107742030B CN201710990538.5A CN201710990538A CN107742030B CN 107742030 B CN107742030 B CN 107742030B CN 201710990538 A CN201710990538 A CN 201710990538A CN 107742030 B CN107742030 B CN 107742030B
- Authority
- CN
- China
- Prior art keywords
- simulation
- pulse
- electromagnetic
- pulse current
- copper pipe
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a simulation method for intermediate frequency heating and pulse current application of a TP2 internal thread copper pipe, which comprises the following steps: static simulation of the TP2 internal thread copper tube is simulated, namely a coil and the TP2 internal thread copper tube are always relatively static, and the TP2 internal thread copper tube is subjected to intermediate frequency heat treatment and applied pulse current whole simulation; dynamic simulation of the TP2 internal thread copper tube is simulated, namely, when the coil and the TP2 internal thread copper tube move relatively, the TP2 internal thread copper tube is subjected to intermediate frequency heat treatment and applied with pulse current at the same time; in static simulation and dynamic simulation, the temperature field distribution, electromagnetic field distribution and pulse current distribution conditions in the overall and cross section directions of the copper pipe at each moment are obtained, and the distribution rules of each parameter on the temperature field, the electromagnetic field and the pulse current are summarized. The method can be used for carrying out simulation on the experiment of applying pulse current and medium-frequency heat treatment to the TP2 internal thread copper pipe, and the result is relatively accurate and is very close to the effect of an actual experiment.
Description
Technical Field
The invention relates to a simulation method for improving the texture of a threaded copper pipe in TP2, in particular to a method for simulating pulse current to assist in medium-frequency induction heating treatment of a copper pipe, which can realize numerical simulation of applying pulse current while performing medium-frequency heat treatment on the copper pipe.
Background
The TP2 internal thread copper tube is widely applied to national production and life such as air conditioners, refrigerators and the like, and is widely applied to refrigeration and heat conduction industries due to good heat transfer performance. The method mainly adopted for producing the TP2 internal thread copper pipe in China is a casting and rolling method, but due to the limitation of the prior production technology, certain defects generally exist, and in the using process, after the internal thread tooth part playing an important role in heat conductivity is used for a period of time, the phenomena of cracks, thread tooth defect, tooth breakage and the like caused by fatigue damage are often easy to occur, so that the heat conductivity of the threaded pipe is greatly reduced, and the service life of the threaded pipe is greatly prolonged. In order to improve the performance of the copper pipe, the processing difficulty and the investment of manpower and material resources which need to be overcome in the production process are continuously increased, which both increase the burden and the cost of enterprises.
The existing research results show that when the pulse current acts on the high-temperature metal structure, the structure can be refined, and the mechanical property can be improved; on the other hand, the electromagnetic induction heating treatment has the advantages of fast heating, high efficiency and small pollution, and is generally applied to industrial production, and by combining the two points, the pulse current is applied while the intermediate frequency heat treatment is carried out on the copper pipe, so that the structure of the copper pipe can be refined, and the mechanical property of the copper pipe can be improved. Thereby avoiding the defects of high processing difficulty, high production cost and the like caused by reducing the influence caused by the defects of the production process.
However, under the combined action of induction heating and pulse current application to the TP2 internal thread copper tube, the effect of obtaining the structure performance of the copper tube is generally summarized and verified by a large number of experiments, although the influence and change of the structure performance of the copper tube under two physical actions can be really and reliably mastered through the experiments, the organization can be finished by a large amount of manpower and material resources, and much time is spent, more importantly, the experiment can only be used for checking the structure through metallographic experiments, and in the most important experimental process, the temperature field, the electromagnetic field distribution, the distribution condition and the change condition of the pulse current of each part on the copper tube can not be obtained in the experiment, so that the metallographic experiment can only be performed on a test piece after the experiment to obtain the change of the copper tube structure crystal grains under different parameters, thereby obtaining a certain rule, however, in the experimental process, all physical fields playing an important role are not obtained at all, the defects are too large, and ANSYS software has the coupling of multiple physical fields such as an electromagnetic physical field, a thermoelectric coupling field and the like and a strong numerical simulation function, so that the effect close to the actual effect can be obtained through accurate calculation of an electronic computer, and the ANSYS software is widely used.
Therefore, based on the defects in the experiment and the powerful functions of the computer technology and ANSYS software, the method for simulating the experiment by using the ANSYS software as the simulation platform is an excellent method which saves cost, time, manpower and material resources and is more accurate and close to the actual situation, more importantly, a large amount of preparation work and guiding work are made for the experiment, a corresponding rule is well laid for the experiment summary, an action mechanism of a main physical field in the experiment process is vividly, vividly and more accurately given, and important support and materials are provided for the summary rule.
Disclosure of Invention
The invention aims to provide a simulation calculation method for performing simulation calculation on the whole experimental process of a TP2 internal thread copper pipe based on an ANSYS simulation platform, accurately simulating two physical actions of intermediate frequency heat treatment and pulse current loading for improving the performance of the copper pipe, and being closer to the actual experimental process.
In order to realize the purpose, the invention is realized by the following technical scheme:
a simulation method for intermediate frequency heating and pulse current application of a TP2 internal thread copper pipe comprises static simulation and dynamic simulation;
the static simulation comprises the following steps:
extracting geometric parameters and physical parameters of a copper pipe and an induction coil, and converting the geometric parameters and the physical parameters into a command stream language in ANSYS; establishing a physical model of the copper pipe, the induction coil and the air according to the parameters; carrying out mesh division to obtain a finite element model;
step two, importing a finite element model, and entering a physical environment library to simulate each physical action: the concrete content is as follows:
(1) medium frequency induction heating was simulated by a do-cycle language command:
introducing an electromagnetic physical environment of medium-frequency induction heating, reading in a temperature load and taking the temperature load as an initial condition, and starting electromagnetic analysis and calculation; after the electromagnetic analysis is finished, introducing a medium-frequency induction heating thermal physical environment, and performing thermal calculation by taking the temperature load in the electromagnetic analysis and the heat generation rate in the electromagnetic analysis as initial conditions; after the calculation is finished, under a dot-do circulation command, returning to the beginning of electromagnetic analysis, restarting again, and sequentially repeating the operation according to the set circulation times until the required heat treatment temperature range is reached;
(2) the pulse current was simulated by a do loop language command:
first, pulsed thermoelectric simulation: importing a pulse thermoelectric physical environment, converting pulse current key parameters, namely pulse voltage, pulse frequency and pulse width, into ANSYS language, importing temperature load in induction heating as an initial condition, and calculating; after the specified calculation time is reached, storing the calculation result into a corresponding result file;
secondly, pulse electromagnetic simulation: introducing a pulse electromagnetic physical environment, setting parameters of pulse voltage, pulse frequency and pulse width of pulse current through an APDL (advanced peripheral component description language) and a dot cycle command, taking a temperature load of pulse thermoelectricity as a temperature initial condition, introducing an electromagnetic force load of induction heating as an electromagnetic field initial condition, and calculating; after the specified calculation time is reached, storing the calculation result into a corresponding result file;
step three, clearing away the finite element model introduced in the simulation of the step one and the step two at the last time; extracting geometric parameters and physical parameters of the copper pipe and the newly established induction coil again, and converting the geometric parameters and the physical parameters into a command stream language in ANSYS; establishing a physical model of the copper pipe, the induction coil and the air according to the parameters; carrying out mesh division again to generate a new finite element model;
step four, importing a new finite element model, and entering a physical environment library to simulate each physical action: the concrete content is as follows:
(1) medium frequency induction heating was simulated by a do-cycle language command:
introducing an electromagnetic physical environment of medium-frequency induction heating, reading in the temperature load at the final moment in the last simulation as an initial condition, and starting electromagnetic analysis and calculation; after the electromagnetic analysis is finished, introducing a medium-frequency induction heating thermal physical environment, and performing thermal calculation by taking the temperature load in the electromagnetic analysis and the heat generation rate in the electromagnetic analysis as initial conditions; after the calculation is finished, under a dot-do circulation command, returning to the beginning of electromagnetic analysis, restarting again, and sequentially repeating the operation according to the set circulation times until the required heat treatment temperature range is reached;
(2) the pulse current was simulated by a do loop language command:
first, pulsed thermoelectric simulation: importing a pulse thermoelectric physical environment, converting pulse current key parameters, namely pulse voltage, pulse frequency and pulse width, into ANSYS language, importing temperature load in induction heating as an initial condition, and calculating; when the specified calculation time is up, the calculation result is saved in a corresponding result file;
secondly, pulse electromagnetic simulation: introducing a pulse electromagnetic physical environment, setting parameters of pulse voltage, pulse frequency and pulse width of pulse current through an APDL (advanced peripheral component description language) and a dot cycle command, taking a temperature load of pulse thermoelectricity as a temperature initial condition, introducing an electromagnetic force load of induction heating as an electromagnetic field initial condition, and calculating; when the specified calculation time is up, the calculation result is saved in a corresponding result file;
repeating the steps in sequence until the specified time is reached according to the time setting in the static simulation;
the dynamic simulation has the same basic steps as the static simulation, and only in an ANSYS software platform, in order to realize the relative motion effect between the copper pipe and the induction coil, the relative position between the induction coil and the copper pipe is changed once when the induction coil moves once, and the corresponding finite element model needs to be established again;
the dynamic simulation comprises the following steps:
step one, the same as static simulation step one;
step two, the same static simulation step two;
step three, clearing away the finite element model introduced in the simulation of the step one and the step two at the last time; extracting geometric parameters and physical parameters of the copper pipe and the newly established induction coil, and converting the geometric parameters and the physical parameters into a command stream language in ANSYS; establishing a physical model of the copper pipe, the induction coil and the air according to the parameters; carrying out mesh division again to generate a new finite element model; compared with the induction coil which is newly established in the last simulation, the size of the induction coil is not changed, and only the relative position of the induction coil and the copper pipe is changed;
step four, importing a new finite element model, and entering a physical environment library to simulate each physical action: the concrete content is as follows:
(1) medium frequency induction heating was simulated by a do-cycle language command:
introducing an electromagnetic physical environment of medium-frequency induction heating, reading in the temperature load at the final moment in the last simulation as an initial condition, and starting electromagnetic analysis and calculation; after the electromagnetic analysis is finished, introducing a medium-frequency induction heating thermal physical environment, and performing thermal calculation by taking the temperature load in the electromagnetic analysis and the heat generation rate in the electromagnetic analysis as initial conditions; after the calculation is finished, under a dot-do circulation command, returning to the beginning of electromagnetic analysis, restarting again, and sequentially repeating the operation according to the set circulation times until the required heat treatment temperature range is reached;
(2) the pulse current was simulated by a do loop language command:
first, pulsed thermoelectric simulation: importing a pulse thermoelectric physical environment, converting pulse current key parameters, namely pulse voltage, pulse frequency and pulse width, into ANSYS language, importing temperature load in induction heating as an initial condition, and calculating; when the specified calculation time is up, the calculation result is saved in a corresponding result file;
secondly, pulse electromagnetic simulation: introducing a pulse electromagnetic physical environment, setting parameters of pulse voltage, pulse frequency and pulse width of pulse current through an APDL (advanced peripheral component description language) and a dot cycle command, taking a temperature load of pulse thermoelectricity as a temperature initial condition, introducing an electromagnetic force load of induction heating as an electromagnetic field initial condition, and calculating; when the specified calculation time is up, the calculation result is saved in a corresponding result file;
and repeating the steps in sequence until the specified time is reached according to the time setting in the dynamic simulation.
The static simulation process can be automatically completed by a computer only by inputting a short APDL language on an ANSYS software platform, and the specific process is as follows:
1. importing a finite element model;
2. executing one to a plurality of ANSYS macro libraries;
3. the individual macros in the macro library are executed in order.
According to the dynamic simulation, the specific process only needs to input a short APDL command in ANSYS, so that a computer can automatically calculate to complete the whole dynamic simulation process, and the specific process is as follows:
1. importing a finite element model;
2. executing one to a plurality of ANSYS macro libraries;
3. executing each macro in the macro library in sequence;
4. removing the model, moving the coil and generating a new integral model;
5. commands 1 to 4 are executed again until the prescribed time is reached.
The static simulation and the dynamic simulation are that pulse current is applied to the TP2 internal thread copper pipe while the intermediate frequency heat treatment is carried out by adopting a sequential coupling method; moreover, in order to realize the equivalent effect of the medium-frequency heat treatment and the pulse current simultaneously acting on the copper tube, in each of the static simulation and the dynamic simulation, the time for the medium-frequency heat treatment simulation and the pulse current simulation must be the same, and the equivalent effect that the medium-frequency heat treatment and the pulse current simultaneously act on the copper tube becomes more obvious the shorter the time.
In each simulation of the static simulation and the dynamic simulation, the time for the intermediate frequency heat treatment simulation and the pulse current simulation is not more than 1 s.
Compared with the prior art, the invention has the following advantages:
by adopting the method, the simulation of the experiment of applying the pulse current and the medium-frequency heat treatment effect to the TP2 internal thread copper pipe can be carried out, the result is more accurate and is very close to the effect of the actual experiment; the method can take the data of each parameter group in the experiment to carry out simulation, and can obtain the influence rule and mechanism of each parameter on the copper pipe in the whole experiment, thereby obtaining the optimal parameter collocation; the method provides precious and abundant experience and rules for the actual experiment, so that experimenters can make certain prejudgment on the rules obtained by the actual experiment and lay a cushion for the actual experiment effect; for data of some parameters beyond the actual experiment specification, due to dangerousness or experiments which are difficult to reach by experimental equipment, the simulation can be carried out by the method, and a certain material is provided for the aspect. The method can obtain the temperature distribution, the electromagnetic force distribution, the magnetic induction intensity, the magnetic flux density, the pulse current density distribution cloud chart and the corresponding node data of the TP2 internal thread copper pipe at all times on the whole and the cross section, provides the capability of researching the change of the electromagnetic field, the temperature field and the pulse current distribution of all parts on the copper pipe in real time for the experiment, can accurately and vividly show the change influence mechanism of all parameter groups on the electromagnetic field, the temperature field and the pulse current in the whole experiment, and is hard to reach in the actual experiment. The invention can also carry out simulation on the experiment of various pipes, plates and the like made of other metal materials under the action of induction heating and pulse current, and has universality; experimenters can input simple programming languages to carry out the whole experiment through a computer, and the operation is simple and convenient, and a large amount of manpower and material resources are saved.
Drawings
FIG. 1 is a flow chart of simulation of the whole experiment of the invention in a static state of a copper pipe;
FIG. 2 is a flow chart of the simulation of the whole experiment of the invention during the dynamic state of the copper pipe;
FIG. 3 is a finite element model diagram of a TP2 internal threaded copper tube according to the present invention;
FIG. 4 is a schematic diagram of a cloud of the overall temperature field of the copper tube at a certain time in the simulation of the present invention;
FIG. 5 is a schematic illustration of a cloud of a copper tube cross-sectional temperature field at a certain time in a simulation of the present invention;
FIG. 6 is a schematic view of the overall electromagnetic force cloud of the copper tube at a certain time in the simulation of the present invention;
FIG. 7 is a schematic diagram of a copper tube cross-section electromagnetic force cloud at a certain point in the simulation of the present invention;
FIG. 8 is a schematic diagram of a cloud of the distribution of the overall pulse current of the copper tube at a certain time in the simulation of the present invention;
FIG. 9 is a schematic diagram of a cloud of the distribution of pulse current in a cross section of a copper tube at a certain time in the simulation of the present invention.
Detailed Description
The invention is further described below with reference to the accompanying drawings:
the invention relates to a simulation method for intermediate frequency heating and pulse current application of a TP2 internal thread copper pipe, which comprises static simulation and dynamic simulation;
as shown in fig. 1, the static simulation includes the following steps: extracting geometric parameters and physical parameters of a copper pipe and an induction coil, and converting the geometric parameters and the physical parameters into a command stream language in ANSYS; establishing a physical model of the copper pipe, the coil and the air according to the parameters; carrying out mesh division to obtain a finite element model;
extracting geometric parameters and physical parameters of a copper pipe and an induction coil according to experimental requirements, and setting each group of key parameters of induction heating and pulse current; like the experiment, set up each group's parameter of induction heating and pulse current, also take different numerical values to set up the parameter group in the important parameter in two physical effects in the numerical simulation, wherein the key parameter in the induction heating is: the heating frequency and the alternating current density are two, so that a proper value range is selected based on experiments, and the same number of proper values are respectively taken for the two parameters in the value range, so that two parameter sets are generated; similarly, the three parameters important in the pulse current are: the pulse voltage, the pulse width and the pulse frequency are obtained by respectively selecting appropriate numerical values with the same number for the three parameters based on experiments, and setting three parameter groups; extracting various parameters of the copper pipe and the induction coil, such as geometric parameters of an inner diameter, an outer diameter, a length, a thread angle and the like, and physical parameters of relative magnetic conductivity, resistivity, a heat conductivity coefficient, a specific heat capacity and the like, and respectively storing the physical parameters in a file according to physical attributes after extraction.
Each set of key parameters is converted into a command stream language in ANSYS. Establishing a physical model of the copper pipe, the coil and the air according to the parameters; carrying out grid division to obtain a finite element model, and storing the model into a database;
opening ANSYS software, and performing the following operations in the ANSYS software:
reading each material file, and importing each material characteristic into an ANSYS material library; establishing a physical model of a copper pipe, a coil and air; carrying out mesh division on each entity model, endowing material characteristics and setting unit types to corresponding meshes to obtain a finite element model, and storing the model into a database; establishing a plurality of ANSYS macro libraries which are core parts of the whole simulation, wherein the ANSYS macro libraries are established:
establishing a physical environment macro through APDL language programming: establishing a physical environment for the finite element model, and establishing an electromagnetic environment, a thermal environment, a pulse thermoelectric environment and a pulse electromagnetic environment in sequence according to two physical actions of induction heating and pulse current. The electromagnetic physical environment endows the finite element model with an electromagnetic unit, sets boundary conditions, a load range and a loading mode, an electromagnetic calculation convergence criterion and a solving type; the thermal physical environment is a heat unit assigned for finite element model setting, a thermal convergence criterion, load step setting and solving type; the two physical environments are all physical environments when medium-frequency induction heating is carried out.
Step two, importing a finite element model, and entering a physical environment library to simulate each physical action:
1. medium frequency induction heating was simulated by a do-cycle language command:
and (3) introducing an electromagnetic physical environment, reading a temperature load as an initial condition, starting calculation, introducing a thermal physical environment after electromagnetic analysis is finished, taking the temperature load in the electromagnetic analysis as the initial condition, introducing thermal generation rate data in an electromagnetic analysis result file, and calculating. And after the calculation is finished, under the command of the dot-do cycle, returning to the beginning of the electromagnetic analysis, restarting again, and repeating the steps in sequence according to the set cycle times until the required heat treatment temperature range is reached.
2. The pulsed current is simulated by a do-cycle language command, wherein the pulsed current causes a temperature increase and generates an electromagnetic field, thus requiring two-part simulation thereof, namely, pulsed thermoelectric simulation and pulsed electromagnetic simulation.
First, pulsed thermoelectric simulation: and (3) introducing a pulse thermoelectric physical environment, converting key parameters of pulse current, namely pulse voltage, pulse frequency and pulse width, into scalar parameters of ANSYS, introducing temperature load in induction heating as an initial condition, and starting calculation.
Secondly, pulse electromagnetic simulation: and (2) introducing a pulse electromagnetic physical environment, setting important parameters of pulse voltage, pulse frequency and pulse width of pulse current through an APDL (advanced design description language) and a dot cycle command as same as pulse thermoelectricity, introducing a temperature load of the pulse thermoelectricity as a temperature initial condition and an electromagnetic force load of induction heating as an electromagnetic field initial condition after setting, starting calculation, and realizing the simulation of the whole experiment in pulse electromagnetic simulation.
Step three, clearing away the finite element model introduced in the simulation of the step one and the step two at the last time; extracting geometric parameters and physical parameters of the copper pipe and the newly established induction coil again, and converting the geometric parameters and the physical parameters into a command stream language in ANSYS; establishing a physical model of the copper pipe, the induction coil and the air according to the parameters; carrying out mesh division again to generate a new finite element model;
step four, importing a new finite element model, and entering a physical environment library to simulate each physical action: the concrete content is as follows:
1. medium frequency induction heating was simulated by a do-cycle language command:
introducing an electromagnetic physical environment of medium-frequency induction heating, reading in the temperature load at the final moment in the last simulation as an initial condition, and starting electromagnetic analysis and calculation; after the electromagnetic analysis is finished, introducing a medium-frequency induction heating thermal physical environment, and performing thermal calculation by taking the temperature load in the electromagnetic analysis and the heat generation rate in the electromagnetic analysis as initial conditions; after the calculation is finished, under a dot-do circulation command, returning to the beginning of electromagnetic analysis, restarting again, and sequentially repeating the operation according to the set circulation times until the required heat treatment temperature range is reached;
2. the pulse current was simulated by a do loop language command:
first, pulsed thermoelectric simulation: importing a pulse thermoelectric physical environment, converting pulse current key parameters, namely pulse voltage, pulse frequency and pulse width, into ANSYS language, importing temperature load in induction heating as an initial condition, and calculating; when the specified calculation time is up, the calculation result is saved in a corresponding result file;
secondly, pulse electromagnetic simulation: introducing a pulse electromagnetic physical environment, setting parameters of pulse voltage, pulse frequency and pulse width of pulse current through an APDL (advanced peripheral component description language) and a dot cycle command, taking a temperature load of pulse thermoelectricity as a temperature initial condition, introducing an electromagnetic force load of induction heating as an electromagnetic field initial condition, and calculating; when the specified calculation time is up, the calculation result is saved in a corresponding result file;
the above-mentioned specific procedure for static simulation is that, when the whole experiment is simulated, the computer can automatically calculate and complete the whole simulation only by inputting a short APDL language command on an ANSYS software platform, as follows:
1. importing a finite element model;
2. performing ANSYS macro library (one to several);
3. the individual macros in the macro library are executed in order.
The whole static simulation of the copper pipe can be realized.
Fig. 2 shows a dynamic simulation process of a TP2 internal thread copper pipe, which is basically the same as a static simulation process, except that the dynamic simulation makes the copper pipe achieve an equivalent motion effect, and the specific process is as follows:
firstly, extracting geometric parameters and physical parameters of a copper pipe and an induction coil, and converting the geometric parameters and the physical parameters into a command stream language in ANSYS. Establishing a physical model of the copper pipe, the coil and the air according to the parameters; and carrying out mesh division to obtain a finite element model.
And secondly, importing a finite element model, and entering a physical environment library to simulate each physical action:
the pulse thermoelectric physical environment endows the finite element model with a thermoelectric coupling unit and one end surface of the selected model, and couples voltage freedom degrees to all nodes on the thermoelectric coupling unit and the selected model; the pulse electromagnetic physical environment endows the finite element model with an electromagnetic unit and sets boundary conditions. The two physical environments are all physical environments when the pulse current is applied.
Then, establishing an analog macro through APDL language:
1. medium frequency induction heating was simulated by a do-cycle language command:
introducing scalar parameters into a dot-do circulation command, storing all updated scalar parameters into a text, merging the scalar parameters into a database, introducing an electromagnetic physical environment, and storing a mode and a process of electromagnetic analysis to be performed and an obtained result in a corresponding file form in ANSYS, so that file names are named before the electromagnetic analysis, then temperature loads are introduced and used as initial conditions, calculation is started, after the electromagnetic analysis is completed, a thermal physical environment is introduced, corresponding files in the subsequent thermal analysis are named, the initial temperature conditions are introduced, thermal generation rate data in an electromagnetic analysis result file are introduced, a parameter file is read, and the temperature loads in the electromagnetic analysis and the thermal generation rate in the electromagnetic analysis result are used as initial conditions to perform thermal calculation; and after the thermal analysis is calculated, returning to the beginning electromagnetic analysis under the command of the dot-do cycle, restarting again, and repeating the cycle times set according to the command of the dot-do in sequence until the required thermal treatment temperature range is reached. The above-mentioned steps in the first cycle of the start are the same as the first cycle in each subsequent cycle, except that the first cycle of the start, the temperature initial conditions in each electromagnetic analysis are the temperature results obtained by the thermal analysis in the previous cycle, and the temperature initial conditions in each thermal analysis are also the results obtained by the thermal analysis in the previous cycle.
2. The pulse current was simulated by a do-cycle language command. In which a pulsed current causes a temperature increase and an electromagnetic field is generated, and therefore it is necessary to simulate it in two parts, namely, a pulsed thermoelectric simulation and a pulsed electromagnetic simulation. First, pulsed thermoelectric simulation: and importing a pulse thermoelectric physical environment, naming each file generated in the pulse thermoelectric process in advance, converting key parameters of pulse current, namely pulse voltage, pulse frequency and pulse width, into scalar parameters of ANSYS, setting the scalar parameters through APDL language programming and dot-do circulation commands, importing the temperature load in induction heating as an initial condition after setting, and starting to calculate. It should be noted that, in the present invention, there are many physical fields involved in the intermediate frequency heat treatment and the simulation of the pulse current, so the present invention adopts the sequential coupling method to realize the whole experimental process of applying the pulse current to the intermediate frequency heat treatment of the TP2 internal thread copper pipe at the same time. And in order to realize the effect close to the actual situation, the time for each intermediate frequency heat treatment simulation and the pulse current simulation is the same and short.
Secondly, pulse electromagnetic simulation: and (2) introducing a pulse electromagnetic physical environment, naming each file generated in pulse electromagnetism in advance, setting each important parameter of pulse current by an APDL (advanced programming language) and a dot cycle command as same as pulse thermoelectricity, introducing a temperature load of the pulse thermoelectricity as a temperature initial condition and an electromagnetic force load of induction heating as an electromagnetic field initial condition after setting, starting calculation, and realizing the simulation of the whole experiment in pulse electromagnetic simulation.
And thirdly, after the specified calculation time is reached, storing the calculation result into a corresponding result file, removing the simulated finite element model, reintroducing the copper pipe entity model which is not subjected to simulated calculation, and establishing a new coil, wherein the size of the coil is not changed, only the relative position of the coil and the copper pipe is changed, namely the movement of the copper pipe is equivalent through the stepping movement of the coil. And re-meshing, generating a new integral finite element model, and then sequentially carrying out simulation calculation according to the flow until the specified time is reached.
The above-mentioned specific process for dynamic simulation is the same as static simulation in the whole experiment, and only a short APDL language command needs to be input in the ANSYS software platform to make the computer automatically calculate and complete the whole simulation, as follows:
1. importing a finite element model
2. Executing ANSYS macro library (one to several)
3. Executing each macro in macro library in sequence
4. Removing the model, moving the coil, and creating a new overall model
5. Commands 1 to 4 are executed again until the prescribed time is reached.
FIG. 3 is a finite element model diagram of a TP2 internal threaded copper tube according to the present invention; when ANSYS finite element software is used for carrying out static simulation and dynamic simulation on the TP2 internal thread copper pipe, a three-dimensional model passing through the copper pipe is subjected to grid division to obtain a finite element model, the finite element model is used as an entity, and corresponding load and boundary conditions are applied to the entity according to actual physical conditions, so that simulation is realized.
As shown in fig. 4, the cloud chart is a schematic cloud chart of the temperature distribution of the whole copper pipe at each moment obtained after the static simulation or the dynamic simulation of the copper pipe is completed. Through the acquisition of the cloud picture of the overall temperature of the copper pipe, the change condition of the overall temperature of the copper pipe with the internal thread of TP2 and the influence rule of different parameter changes on the overall temperature change of the copper pipe in the whole static simulation and dynamic simulation processes can be obtained, and the parameter rule obtained from the change condition can be used for prejudging the overall temperature distribution condition of the copper pipe.
As shown in fig. 5, a cloud chart is a schematic diagram of the temperature distribution on the cross section of the copper tube at each moment obtained after the static simulation or the dynamic simulation of the copper tube is completed. By acquiring the temperature cloud pictures on the cross section of the copper pipe, the change condition of the temperature on the cross section of the copper pipe with the internal thread of TP2 and the influence rule of different parameter changes on the whole temperature change of the copper pipe in the whole static simulation and dynamic simulation processes can be obtained, and the parameter rule obtained from the change condition can be used for prejudging the temperature distribution condition on the cross section of the copper pipe.
As shown in fig. 6, the schematic diagram is an electromagnetic force cloud chart of the whole copper tube at each moment obtained after the static simulation or the dynamic simulation of the copper tube is completed. Through obtaining the overall electromagnetic force cloud picture of the copper pipe, the change condition of the overall electromagnetic force of the internal thread copper pipe in TP2 and the influence rule of different parameter changes on the overall electromagnetic force change of the copper pipe in the whole static simulation and dynamic simulation processes can be obtained, and the parameter rule obtained from the change condition can be used for making a prejudgment on the overall electromagnetic force distribution condition of the copper pipe.
As shown in fig. 7, the schematic diagram of the electromagnetic force cloud chart of the cross section of the copper tube at each moment is obtained after the static simulation or the dynamic simulation of the copper tube is completed. By obtaining the electromagnetic force cloud picture of the cross section of the copper pipe, the change condition of the electromagnetic force of the cross section of the copper pipe with the internal thread of TP2 and the influence rule of different parameter changes on the change of the electromagnetic force of the cross section of the copper pipe in the whole static simulation and dynamic simulation processes can be obtained, and the distribution condition of the electromagnetic force on the cross section of the copper pipe can be pre-judged according to the obtained parameter rule.
As shown in fig. 8, the schematic diagram is a cloud chart of the distribution of the whole pulse current of the copper tube at each time obtained after the static simulation or the dynamic simulation of the copper tube is completed. By obtaining the schematic diagram of the copper pipe integral pulse current distribution cloud picture, the change condition of the integral pulse current of the TP2 internal thread copper pipe and the influence rule of different parameter changes on the change of the copper pipe integral pulse current in the whole static simulation and dynamic simulation processes can be obtained, and the obtained parameter rule can make prejudgment on the distribution condition of the copper pipe integral pulse current. And the distribution conditions of the integral temperature cloud picture and the electromagnetic force cloud picture of the copper pipe at corresponding moments are compared, so that the influence rule of the integral temperature distribution and the electromagnetic force distribution of the copper pipe on the integral pulse current distribution of the copper pipe can be obtained.
As shown in fig. 9, the schematic diagram of the distribution cloud chart of the pulse current of the cross section of the copper tube at each time is obtained after the static simulation or the dynamic simulation of the copper tube is completed. By obtaining the schematic diagram of the copper pipe cross section pulse current distribution cloud picture, the change condition of the TP2 internal thread copper pipe cross section pulse current and the influence rule of different parameter changes on the copper pipe cross section pulse current change in the whole static simulation and dynamic simulation processes can be obtained, and the parameter rule obtained from the change condition can be used for prejudging the copper pipe cross section pulse current distribution condition. And the distribution conditions of the copper pipe cross section temperature cloud picture and the electromagnetic force cloud picture at corresponding moments are compared, so that the influence rule of the copper pipe cross section temperature distribution and the electromagnetic force distribution on the copper pipe cross section pulse current distribution can be obtained.
Claims (4)
1. A simulation method for intermediate frequency heating and pulse current application of a TP2 internal thread copper pipe comprises static simulation and dynamic simulation;
the static simulation comprises the following steps:
extracting geometric parameters and physical parameters of a copper pipe and an induction coil, and converting the geometric parameters and the physical parameters into a command stream language in ANSYS; establishing a physical model of the copper pipe, the induction coil and the air according to the parameters; carrying out mesh division to obtain a finite element model;
step two, importing a finite element model, and entering a physical environment library to simulate each physical action: the concrete content is as follows:
(1) medium frequency induction heating was simulated by a do-cycle language command:
introducing an electromagnetic physical environment of medium-frequency induction heating, reading in a temperature load and taking the temperature load as an initial condition, and starting electromagnetic analysis and calculation; after the electromagnetic analysis is finished, introducing a medium-frequency induction heating thermal physical environment, and performing thermal calculation by taking the temperature load in the electromagnetic analysis and the heat generation rate in the electromagnetic analysis as initial conditions; after the calculation is finished, under a dot-do circulation command, returning to the beginning of electromagnetic analysis, restarting again, and sequentially repeating the operation according to the set circulation times until the required heat treatment temperature range is reached;
(2) the pulse current was simulated by a do loop language command:
first, pulsed thermoelectric simulation: importing a pulse thermoelectric physical environment, converting pulse current key parameters, namely pulse voltage, pulse frequency and pulse width, into ANSYS language, importing temperature load in induction heating as an initial condition, and calculating; after the specified calculation time is reached, storing the calculation result into a corresponding result file;
secondly, pulse electromagnetic simulation: introducing a pulse electromagnetic physical environment, setting parameters of pulse voltage, pulse frequency and pulse width of pulse current through an APDL (advanced peripheral component description language) and a dot cycle command, taking a temperature load of pulse thermoelectricity as a temperature initial condition, introducing an electromagnetic force load of induction heating as an electromagnetic field initial condition, and calculating; after the specified calculation time is reached, storing the calculation result into a corresponding result file;
step three, clearing away the finite element model introduced in the simulation of the step one and the step two at the last time; extracting geometric parameters and physical parameters of the copper pipe and the newly established induction coil again, and converting the geometric parameters and the physical parameters into a command stream language in ANSYS; establishing a physical model of the copper pipe, the induction coil and the air according to the parameters; carrying out mesh division again to generate a new finite element model;
step four, importing a new finite element model, and entering a physical environment library to simulate each physical action: the concrete content is as follows:
(1) medium frequency induction heating was simulated by a do-cycle language command:
introducing an electromagnetic physical environment of medium-frequency induction heating, reading in the temperature load at the final moment in the last simulation as an initial condition, and starting electromagnetic analysis and calculation; after the electromagnetic analysis is finished, introducing a medium-frequency induction heating thermal physical environment, and performing thermal calculation by taking the temperature load in the electromagnetic analysis and the heat generation rate in the electromagnetic analysis as initial conditions; after the calculation is finished, under a dot-do circulation command, returning to the beginning of electromagnetic analysis, restarting again, and sequentially repeating the operation according to the set circulation times until the required heat treatment temperature range is reached;
(2) the pulse current was simulated by a do loop language command:
first, pulsed thermoelectric simulation: importing a pulse thermoelectric physical environment, converting pulse current key parameters, namely pulse voltage, pulse frequency and pulse width, into ANSYS language, importing temperature load in induction heating as an initial condition, and calculating; when the specified calculation time is up, the calculation result is saved in a corresponding result file;
secondly, pulse electromagnetic simulation: introducing a pulse electromagnetic physical environment, setting parameters of pulse voltage, pulse frequency and pulse width of pulse current through an APDL (advanced peripheral component description language) and a dot cycle command, taking a temperature load of pulse thermoelectricity as a temperature initial condition, introducing an electromagnetic force load of induction heating as an electromagnetic field initial condition, and calculating; when the specified calculation time is up, the calculation result is saved in a corresponding result file;
repeating the steps in sequence until the specified time is reached according to the time setting in the static simulation;
the dynamic simulation has the same basic steps as the static simulation, and only in an ANSYS software platform, in order to realize the relative motion effect between the copper pipe and the induction coil, the relative position between the induction coil and the copper pipe is changed once when the induction coil moves once, and the corresponding finite element model needs to be established again;
the dynamic simulation comprises the following steps:
step one, the same as static simulation step one;
step two, the same static simulation step two;
step three, clearing away the finite element model introduced in the simulation of the step one and the step two at the last time; extracting geometric parameters and physical parameters of the copper pipe and the newly established induction coil, and converting the geometric parameters and the physical parameters into a command stream language in ANSYS; establishing a physical model of the copper pipe, the induction coil and the air according to the parameters; carrying out mesh division again to generate a new finite element model; compared with the induction coil which is newly established in the last simulation, the size of the induction coil is not changed, and only the relative position of the induction coil and the copper pipe is changed;
step four, importing a new finite element model, and entering a physical environment library to simulate each physical action: the concrete content is as follows:
(1) medium frequency induction heating was simulated by a do-cycle language command:
introducing an electromagnetic physical environment of medium-frequency induction heating, reading in the temperature load at the final moment in the last simulation as an initial condition, and starting electromagnetic analysis and calculation; after the electromagnetic analysis is finished, introducing a medium-frequency induction heating thermal physical environment, and performing thermal calculation by taking the temperature load in the electromagnetic analysis and the heat generation rate in the electromagnetic analysis as initial conditions; after the calculation is finished, under a dot-do circulation command, returning to the beginning of electromagnetic analysis, restarting again, and sequentially repeating the operation according to the set circulation times until the required heat treatment temperature range is reached;
(2) the pulse current was simulated by a do loop language command:
first, pulsed thermoelectric simulation: importing a pulse thermoelectric physical environment, converting pulse current key parameters, namely pulse voltage, pulse frequency and pulse width, into ANSYS language, importing temperature load in induction heating as an initial condition, and calculating; when the specified calculation time is up, the calculation result is saved in a corresponding result file;
secondly, pulse electromagnetic simulation: introducing a pulse electromagnetic physical environment, setting parameters of pulse voltage, pulse frequency and pulse width of pulse current through an APDL (advanced peripheral component description language) and a dot cycle command, taking a temperature load of pulse thermoelectricity as a temperature initial condition, introducing an electromagnetic force load of induction heating as an electromagnetic field initial condition, and calculating; when the specified calculation time is up, the calculation result is saved in a corresponding result file;
repeating the steps in sequence until the specified time is reached according to the time setting in the dynamic simulation;
the static simulation and the dynamic simulation are that pulse current is applied to the TP2 internal thread copper pipe while the intermediate frequency heat treatment is carried out by adopting a sequential coupling method; moreover, in order to realize the equivalent effect of the medium-frequency heat treatment and the pulse current simultaneously acting on the copper tube, in each of the static simulation and the dynamic simulation, the time for the medium-frequency heat treatment simulation and the pulse current simulation must be the same, and the equivalent effect that the medium-frequency heat treatment and the pulse current simultaneously act on the copper tube becomes more obvious the shorter the time.
2. The method for simulating intermediate frequency heating and pulse current application to a TP2 internal thread copper tube according to claim 1, wherein the method comprises the following steps: the static simulation process can be automatically completed by a computer only by inputting a short APDL language on an ANSYS software platform, and the specific process is as follows:
1) importing a finite element model;
2) executing one to a plurality of ANSYS macro libraries;
3) the individual macros in the macro library are executed in order.
3. The method for simulating intermediate frequency heating and pulse current application to a TP2 internal thread copper tube according to claim 1, wherein the method comprises the following steps: according to the dynamic simulation, the specific process only needs to input a short APDL command in ANSYS, so that a computer can automatically calculate to complete the whole dynamic simulation process, and the specific process is as follows:
1) importing a finite element model;
2) executing one to a plurality of ANSYS macro libraries;
3) executing each macro in the macro library in sequence;
4) removing the model, moving the coil and generating a new integral model;
5) commands 1 to 4 are executed again until the prescribed time is reached.
4. The method for simulating intermediate frequency heating and pulse current application to a TP2 internal thread copper tube according to claim 1, wherein the method comprises the following steps: in each simulation of the static simulation and the dynamic simulation, the time for the intermediate frequency heat treatment simulation and the pulse current simulation is not more than 1 s.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710990538.5A CN107742030B (en) | 2017-10-23 | 2017-10-23 | Simulation method for intermediate frequency heating and pulse current application of TP2 internal thread copper pipe |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710990538.5A CN107742030B (en) | 2017-10-23 | 2017-10-23 | Simulation method for intermediate frequency heating and pulse current application of TP2 internal thread copper pipe |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107742030A CN107742030A (en) | 2018-02-27 |
| CN107742030B true CN107742030B (en) | 2020-05-12 |
Family
ID=61237962
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710990538.5A Active CN107742030B (en) | 2017-10-23 | 2017-10-23 | Simulation method for intermediate frequency heating and pulse current application of TP2 internal thread copper pipe |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107742030B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110135073B (en) * | 2019-05-17 | 2023-06-16 | 沈阳大学 | Ultrahigh-strength aluminum alloy pulse current regulation casting and rolling simulation method |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2269758A1 (en) * | 2009-07-03 | 2011-01-05 | EWM Hightec Welding GmbH | DC Arc welding method and welding device |
| CN104141101A (en) * | 2013-05-10 | 2014-11-12 | 贵州大学 | Pulse current assisted aluminum alloy homogenization heat treatment method |
| CN105543468A (en) * | 2015-12-24 | 2016-05-04 | 燕山大学 | Pulse current-assisted weld joint structure refining device for high-frequency welded pipe |
| CN105740536A (en) * | 2016-01-28 | 2016-07-06 | 燕山大学 | Method for dynamically simulating intermediate-frequency heat treatment process of longitudinal welded pipes |
| CN106077092A (en) * | 2016-06-24 | 2016-11-09 | 燕山大学 | The device of the low electric heating Pulse Electric Current with High Density assistant metal operation of rolling |
-
2017
- 2017-10-23 CN CN201710990538.5A patent/CN107742030B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2269758A1 (en) * | 2009-07-03 | 2011-01-05 | EWM Hightec Welding GmbH | DC Arc welding method and welding device |
| CN104141101A (en) * | 2013-05-10 | 2014-11-12 | 贵州大学 | Pulse current assisted aluminum alloy homogenization heat treatment method |
| CN105543468A (en) * | 2015-12-24 | 2016-05-04 | 燕山大学 | Pulse current-assisted weld joint structure refining device for high-frequency welded pipe |
| CN105740536A (en) * | 2016-01-28 | 2016-07-06 | 燕山大学 | Method for dynamically simulating intermediate-frequency heat treatment process of longitudinal welded pipes |
| CN106077092A (en) * | 2016-06-24 | 2016-11-09 | 燕山大学 | The device of the low electric heating Pulse Electric Current with High Density assistant metal operation of rolling |
Non-Patent Citations (5)
| Title |
|---|
| Coupled electro-magneto-thermal model for induction heating process of a moving billet;Kee-Hyeon Cho;《International Journal of Thermal Sciences》;20120608;第60卷;第195-204页 * |
| Induction heat treatment of a ISO C45 steel bar:Experimental and numerical analysis;I. Magnabosco 等;《Computational Materials Science》;20061231;第35卷;第98-106页 * |
| Simulation Analysis of the Temperature Field in an Induction Launcher;Li San-Qun 等;《IEEE TRANSACTIONS ON PLASMA SCIENCE》;20130531;第41卷(第5期);第1055-1060页 * |
| 脉冲电流对金属材料塑性变形和组织结构与性能的影响;姚可夫 等;《机械强度》;20031231;第25卷(第3期);第2340-342页 * |
| 高频直缝焊管动态中频感应加热过程数值分析及实验研究;赵天旭;《万方数据库.学位论文库》;20170725;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107742030A (en) | 2018-02-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105740536B (en) | The analogy method of dynamic simulation is carried out to straight seam welded pipe intermediate frequency heat treatment process | |
| CN107169166B (en) | Large gear ring mobile induction heating calculation method based on ANSYS simulation platform | |
| CN114919181B (en) | Continuous fiber 3D printing process dynamic simulation and printed part buckling deformation prediction method | |
| CN107742030B (en) | Simulation method for intermediate frequency heating and pulse current application of TP2 internal thread copper pipe | |
| Tan et al. | Hierarchical genetic-particle swarm optimization for bistable permanent magnet actuators | |
| Bönig et al. | Explicit dynamics process simulation of linear coil winding for electric drives production | |
| CN104636537A (en) | Pretreatment system and method for umbilical cable finite element model | |
| CN108959793B (en) | Bevel gear double-frequency segmented circulating induction heating numerical simulation method | |
| CN113065272A (en) | Method for rapidly realizing numerical simulation of welding deformation | |
| CN116484668A (en) | Electron beam additive manufacturing process simulation method | |
| CN106529100A (en) | Solenoid superconducting magnet coil magnetic field simulation method based on finite element analysis | |
| CN116050135A (en) | Superconducting coil simulation method, system and computer storage medium | |
| CN117077473B (en) | Two-dimensional axisymmetric electromagnetic thermal multi-field coupling calculation method based on transmission line method | |
| Dong et al. | Electromagneto-mechanical coupling analysis of a test module in J-TEXT Tokamak during plasma disruption | |
| CN102542118A (en) | Method and device for analyzing mechanical properties of geometric model | |
| Li et al. | Development of space-time-controlled multi-stage pulsed magnetic field forming and manufacturing technology at the WHMFC | |
| CN118013665B (en) | Thermal stress simulation method, device and storage medium for service process of casting mold | |
| Si et al. | Development of simulation program for temperature field of mass concrete structures | |
| Xingkui et al. | Modeling and simulating for multi-pass welding process of large welded structures | |
| CN113255186B (en) | Concrete temperature control cooling water pipe separation type automatic modeling method | |
| CN118332952A (en) | Simulation method for full-flow sintering process of full-size sintering furnace | |
| CN116822287A (en) | Method and system for reducing turn-to-turn short circuit multi-coupling field model of converter transformer | |
| Meng et al. | Effect of powder initial distribution on densification of Ti6A14V powder during hot isostatic pressing | |
| JP2011184239A (en) | Simulation method of sintering process of porous body, and program | |
| Zhao et al. | Study on the conversion method of three-dimensional shape prediction model of BGA solder joints |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20201130 Address after: 233000 Room 301, building 6, Shangli University Science Park, No. 8319, Huangshan Avenue, Bengbu City, Anhui Province Patentee after: Anhui huajingdao Technology Service Co.,Ltd. Address before: Hebei Street West Harbor area, 066004 Hebei city of Qinhuangdao province No. 438 Patentee before: Yanshan University |
|
| TR01 | Transfer of patent right |