Arc finite element modeling method based on reabsorption process

Abstract

A finite element modeling method for an arc based on a reabsorption process, wherein current passes through a gas to generate an arc, and temperature distribution data of the arc is monitored in real time; finite element modeling of the arc is performed, and a two-dimensional longitudinal plane segmentation diagram is drawn with the arc center direction as the x-axis and the y-axis perpendicular to the arc center direction, and the interior of the arc is divided into multiple triangular mesh units, wherein unit attributes include the coordinates of three nodes, the geometric center coordinates of the unit, and the average temperature of the unit; radial units at predetermined positions of the arc are sorted; based on the sorting and temperature distribution data, the unit is determined to be an absorbing state region of a radiation absorbing part or a radiating state region of a radiation emitting part; the absorbed heat is set in the absorbing state region and the radiated heat is set in the radiating state region to correct the arc simulation results.

Description

Technical Field

[0001] The present invention relates to the technical field of arc simulation, and in particular to an arc finite element modeling method based on a reabsorption process. Background Art

[0002] Arc simulation uses computer simulation technology to study arc phenomena and their physical properties. Through numerical simulation, engineers can analyze various arc physical processes, including heat transfer, fluid dynamics, and electromagnetic fields, to optimize applications such as welding and cutting. Arc simulation provides a deep understanding of arc behavior, thereby improving process efficiency and product quality while reducing experimental costs and time.

[0003] In arc simulation, the reabsorption phenomenon is an important factor that affects the radial distribution of arc temperature. The arc reabsorption phenomenon refers to the process in which the radiant energy generated in the arc is reabsorbed by the gas near the arc, which involves the energy transfer inside the arc. When the maximum arc temperature Tmax is greater than 16000K, considering the arc reabsorption process is important for accurately simulating the characteristics of the arc. Current studies have shown that arcs with monotonic radial temperature distribution radiate heat from the highest point of the arc center with temperature Tmax to the part where the temperature drops to 0.83Tmax, while the part of the arc where the temperature drops from 0.83Tmax to 5000K absorbs the radiated heat. However, in arc simulation, the temperature is not monotonically distributed radially at every moment. If a non-monotonic situation occurs at a certain moment, it is easy to cause distortion of subsequent simulation results, and then cause the final simulation results to diverge.

[0004] The above information disclosed in this Background section is only for enhancement of understanding of the background of the invention and therefore it may contain information that does not form the prior art that is already known to a person of ordinary skill in the art. Summary of the Invention

[0005] The present invention provides an arc finite element modeling method based on the reabsorption process. When a non-monotonic distribution of radial temperature occurs, it is determined which part of the arc is in a radiation state and which part is in an absorption state. Then, the radiation domain and the absorption domain are reset, so that the simulation result at the next moment regresses to the monotonic distribution of radial temperature, and finally the simulation regression converges.

[0006] An arc finite element modeling method based on the reabsorption process includes:

[0007] Step 1: Current passes through the gas to generate an arc, and the temperature distribution data of the arc is monitored in real time;

[0008] Step 2: Perform finite element modeling on the arc. Draw a two-dimensional longitudinal plane segmentation diagram with the arc center direction as the x-axis and the y-axis perpendicular to the arc center direction. Divide the interior of the arc into multiple triangular mesh units. The unit attributes include the coordinates of the three nodes, the geometric center coordinates of the unit, and the average temperature of the unit.

[0009] Step 3: Sort the radial units of the arc at the predetermined position, and judge each unit. If the predetermined position is between the maximum horizontal coordinate x of the three nodes of a unit, max and the minimum value of the horizontal coordinate x min If the intersection line passes through the node, write the unit object into the matrix Slab[ ]. If the intersection line passes through the node, write all the units associated with the node. Sort all the units in the matrix Slab[ ] from small to large according to the vertical coordinate value of the geometric center.

[0010] Step 4, determining whether the unit is an absorbing state region of the radiation absorbing part or a radiating state region of the radiation emitting part based on the sorting and temperature distribution data;

[0011] Step 5: Set the absorbed heat in the absorption state region and the radiated heat in the radiation state region to correct the arc simulation results.

[0012] In the arc finite element modeling method based on the reabsorption process, step 4 includes:

[0013] Starting from the first one, calculate the average temperature of every five cells in the Slab[ ] matrix, and store them in the matrix Temp[ ] in sequence. At the same time, calculate the average vertical coordinates of the geometric centers of these five cells, and store them in the corresponding Yaver[ ] in sequence. The first value in the Slab[ ] matrix is ​​the maximum value in the matrix.

[0014] Take the Temp[ ] matrix, and compare each element Temp[i] with the two elements after it, Temp[i+1] and Temp[i+2]. If it is found that Temp[i] and Temp[i+2] are both smaller than Temp[i+1], record the serial number i+1 of Temp[i+1] and write it into the matrix Index[ ].

[0015] In the arc finite element modeling method based on the reabsorption process, if Index[ ] is empty, the arc is judged to present a spatial temperature at a predetermined position that monotonically decreases with increasing distance from the arc center, the temperature distribution space with a temperature not lower than 83% of the arc's maximum temperature Tmax is a radiation state region of the radiation portion during the reabsorption process, the temperature distribution space with a temperature higher than 5000K and lower than 83% of the arc's maximum temperature Tmax is an absorption state region of the radiation portion during the reabsorption process, and the temperature distribution space with a temperature lower than 5000K is determined not to be part of the arc.

[0016] In the arc finite element modeling method based on the reabsorption process, if Index[ ] has 1 element, the initial values ​​of counter 1 and counter 2 are set to 0, and the temperature values ​​of every two units in the Slab[ ] matrix are taken starting from the first unit. If the 83% arc maximum temperature Tmax is between these two values, counter 1 is increased by 1; if 13000K is between these two values, counter 2 is increased by 1. Each time the counter changes, the two values ​​in Yaver[ ] corresponding to the two Temp[ ] are averaged and stored in the matrix Y[ ]. After the traversal is completed, the values ​​of the counters are observed.

[0017] If counter 1 of A = 3 and counter 2 = 1, the arc is judged to have a peak at a predetermined position and not at the arc center. The temperature at the arc center is higher than 83% of the arc's maximum temperature Tmax. When the temperatures before and after the peak are monotonically distributed, the temperature distribution space from the arc center to where the temperature begins to drop and drops to 83% of the arc's maximum temperature Tmax is the radiation state region of the radiation portion emitted during the reabsorption process. The temperature distribution space from the first temperature drop to 83% of the arc's maximum temperature Tmax to the end of the temperature drop to 5000K is the absorption state region of the radiation portion absorbed during the reabsorption process. The remaining temperature distribution space is judged not to be part of the arc.

[0018] If counter 1 = 3 and counter 2 = 3, then the arc is judged to have two peaks at the predetermined position, one of which is at the arc center, and there is a valley position between the two peaks. The temperature at the arc center is higher than 83% of the maximum arc temperature Tmax, the temperature peak outside the arc center is higher than 83% of the maximum arc temperature Tmax, the valley temperature is greater than 5000K and less than 13000K. When the temperatures between the extreme values ​​are monotonically distributed, the temperature distribution space from the arc center to the first time the temperature rises to 13000K is the absorption state region of the radiation part absorbed during the reabsorption process. The temperature distribution space from the first time the temperature rises to 13000K to the temperature drops to 83% of the maximum arc temperature Tmax again is the radiation state region of the radiation part emitted during the reabsorption process. The remaining temperature distribution space is judged not to be part of the arc.

[0019] If counter 1 = 1 and counter 2 = 3, it is judged that the arc has two peaks at the predetermined position, one of which is at the arc center, and there is a valley position between the two peaks. The temperature at the arc center is higher than 83% of the maximum arc temperature Tmax, the non-arc center temperature peak is lower than 83% of the maximum arc temperature Tmax and is greater than 13000K, and the valley temperature is greater than 5000K and less than 13000K. When the temperatures between the extreme values ​​are monotonically distributed, the temperature distribution space from the arc center to the first time the temperature rises to 13000K is the absorption state area of ​​the radiation part absorbed during the reabsorption process, and the temperature distribution space from the first time the temperature rises to 13000K to the temperature drops to 13000K again is the radiation state area of ​​the radiation part emitted during the reabsorption process. The remaining temperature distribution space is determined not to be part of the arc; if none of the above situations occurs.

[0020] In the arc finite element modeling method based on the reabsorption process, if Index has two or more values, it is judged that the arc has three peaks at a predetermined position and one is at the arc center, and there are two valley positions between the three peaks. The temperature at the arc center is 83% of the arc maximum temperature Tmax, the inner non-arc center temperature peak is less than 83% of the arc maximum temperature Tmax and greater than 13000K, the outer non-arc center temperature peak is greater than 5000K and less than 13000K, the inner valley temperature is greater than 13000K and less than 83% of the arc maximum temperature Tmax, and the outer valley temperature is greater than 5000K and less than 13000K. When the temperatures between the extreme values ​​are monotonically distributed, the temperature distribution space from the arc center to the second temperature peak is the absorption state area of ​​the radiation absorbing part during the reabsorption process, and the temperature distribution space outside the second temperature peak is the radiation state area of ​​the radiation emitting part during the reabsorption process.

[0021] In the arc finite element modeling method based on the reabsorption process, step 4 includes:

[0022] Starting from the first one, calculate the average temperature of every five cells in the Slab[ ] matrix, and store them in the matrix Temp[ ] in sequence. At the same time, calculate the average vertical coordinates of the geometric centers of these five cells, and store them in the corresponding Yaver[ ] in sequence. The first value in the Slab[ ] matrix is ​​not the maximum value in the matrix.

[0023] Set the initial value of counter 1, counter 2 and counter 3 to 0, and take the temperature value of every two cells in the Slab[ ] matrix starting from the first cell. If the 83% arc maximum temperature Tmax is between these two values, counter 1 is incremented by 1; if 13000K is between these two values, counter 2 is incremented by 1; if 5000K is between these two values, counter 3 is incremented by 1. Each time the counter changes, the two values ​​in Yaver[ ] corresponding to the two Temp[ ] are averaged and stored in the matrix Y[ ]. After the traversal is completed, observe the value of the counter:

[0024] If counter 1 = 2, counter 2 = 2, and counter 3 = 2, then the arc is judged to have a peak at a predetermined position and not at the arc center, the temperature at the arc center is lower than 5000K, the non-arc center temperature peak is greater than 83% of the arc's maximum temperature Tmax, and the temperatures before and after the peak are monotonically distributed, then the temperature distribution space from the arc center to 13000K is the absorption state region of the radiation absorption part during the reabsorption process, and the temperature distribution space from the first time the temperature exceeds 13000K to the end of the temperature dropping back to 83% of the arc's maximum temperature Tmax is the radiation state region of the radiation emission part during the reabsorption process, and the remaining temperature distribution space is judged not to be part of the arc;

[0025] If counter 1 = 2, counter 2 = 2, and counter 3 = 1, then the arc is judged to have a peak at a predetermined position and not at the arc center, the temperature at the arc center is lower than 13000K but greater than 5000K, the non-arc center temperature peak is greater than 83% of the arc's maximum temperature Tmax, and the temperatures before and after the peak are monotonically distributed, then the temperature distribution space from the arc center to 13000K is the absorption state region of the radiation absorption part during the reabsorption process, and the temperature distribution space from the first time the temperature exceeds 13000K to the end of the temperature dropping back to 83% of the arc's maximum temperature Tmax is the radiation state region of the radiation emission part during the reabsorption process, and the remaining temperature distribution space is judged not to be part of the arc;

[0026] If counter 1 = 2, counter 2 = 1, and counter 3 = 1, then the arc is judged to have a peak at a predetermined position and is not at the arc center, the temperature at the arc center is higher than 13000K but lower than 83% of the arc's maximum temperature Tmax, the non-arc center temperature peak is greater than 83% of the arc's maximum temperature Tmax, and when the temperatures before and after the peak are monotonically distributed, the temperature distribution space from the arc center to where the temperature begins to drop and drops to 83% of the arc's maximum temperature Tmax is the radiation state region of the radiation portion emitted during the reabsorption process, and the temperature distribution space from the first temperature drop to 83% of the arc's maximum temperature Tmax to the end of the temperature drop to 5000K is the absorption state region of the radiation portion absorbed during the reabsorption process, and the remaining temperature distribution space is judged not to be part of the arc;

[0027] If counter 1 = 3, counter 2 = 1, and counter 3 = 1, then the arc is judged to have a peak at a predetermined position and not at the arc center, the temperature at the arc center is higher than 83% of the arc's maximum temperature Tmax, and when the temperatures before and after the peak are monotonically distributed, the temperature distribution space from the arc center to where the temperature begins to drop and drops to 83% of the arc's maximum temperature Tmax is the radiation state region of the radiation portion emitted during the reabsorption process, and the temperature distribution space from the first temperature drop to 83% of the arc's maximum temperature Tmax to the end of the temperature drop to 5000K is the absorption state region of the radiation portion absorbed during the reabsorption process, and the remaining temperature distribution space is judged not to be part of the arc;

[0028] If counter 1 = 1, counter 2 = 1, and counter 3 = 1, then it is judged that the arc has two peaks at the predetermined position, one of which is at the arc center, and there is a valley position between the two peaks. The temperature at the arc center is higher than 83% of the arc maximum temperature Tmax, the temperature peak outside the arc center is higher than 83% of the arc maximum temperature Tmax, the valley temperature is greater than 13000K and lower than 83% of the arc maximum temperature Tmax, and when the temperatures between the extreme values ​​are monotonically distributed, the temperature distribution space from the arc center to the first time the temperature rises to 83% of the arc maximum temperature Tmax is the absorption state region of the radiation part absorbed during the reabsorption process, and the temperature distribution space from the first time the temperature rises to 83% of the arc maximum temperature Tmax to the time the temperature drops to 83% of the arc maximum temperature Tmax again is the radiation state region of the radiation part emitted during the reabsorption process. The remaining temperature distribution space is determined not to be part of the arc.

[0029] If counter 1 = 3, counter 2 = 3, and counter 3 = 1, then it is judged that the arc has two peaks at the predetermined position, one of which is at the arc center, and there is a valley position between the two peaks. The temperature at the arc center is higher than 83% of the maximum arc temperature Tmax, the non-arc center temperature peak is higher than 83% of the maximum arc temperature Tmax, the valley temperature is greater than 5000K and less than 13000K. When the temperatures between the extreme values ​​are monotonically distributed, the temperature distribution space from the arc center to the first time the temperature rises to 13000K is the absorption state region of the radiation part absorbed during the reabsorption process, and the temperature distribution space from the first time the temperature rises to 13000K to the temperature drops to 83% of the maximum arc temperature Tmax again is the radiation state region of the radiation part emitted during the reabsorption process. The remaining temperature distribution space is determined not to be part of the arc.

[0030] If counter 1 = 1, counter 2 = 3, and counter 3 = 1, then it is determined that the arc has two peaks at the predetermined position, one of which is at the arc center, and there is a valley position between the two peaks. The temperature at the arc center is higher than 83% of the maximum arc temperature Tmax, the non-arc center temperature peak is lower than 83% of the maximum arc temperature Tmax and is greater than 13000K, and the valley temperature is greater than 5000K and less than 13000K. When the temperatures between the extreme values ​​are monotonically distributed, the temperature distribution space from the arc center to the first time the temperature rises to 13000K is the absorption state region of the radiation absorbing part during the reabsorption process, and the temperature distribution space from the first time the temperature rises to 13000K to the temperature drops to 13000K again is the radiation state region of the radiation emitting part during the reabsorption process. The remaining temperature distribution space is determined not to be part of the arc.

[0031] If counter 1 = 2, counter 2 = 3, and counter 3 = 1, then the arc is judged to have two peaks at the predetermined position, one of which is at the arc center, and there is a valley position between the two peaks. The temperature at the arc center is lower than 83% of the arc maximum temperature Tmax and higher than 13000K. The non-arc center temperature peak is greater than 83% of the arc maximum temperature Tmax, and the valley temperature is greater than 5000K and less than 13000K. When the temperatures between the extreme values ​​are monotonically distributed, the temperature distribution space from the arc center to the first time the temperature rises to 13000K is the absorption state region of the radiation part absorbed during the reabsorption process. The temperature distribution space from the first time the temperature rises to 13000K to the temperature drops to 83% of the arc maximum temperature Tmax again is the radiation state region of the radiation part emitted during the reabsorption process. The remaining temperature distribution space is judged not to be part of the arc. If none of the above situations occurs, an error is reported.

[0032] In the arc finite element modeling method based on the reabsorption process, integration with Fluent software is performed to avoid divergence of simulation temperature results.

[0033] In the arc finite element modeling method based on the reabsorption process, in step 2, a plane is divided using a grid of 0.5 mm size.

[0034] In the arc finite element modeling method based on the reabsorption process, in step 2, the predetermined position is a position 1.0 mm away from the arc center.

[0035] In the arc finite element modeling method based on the reabsorption process, in step 2, each cross section of the longitudinal plane segmentation diagram has only one layer of grid.

[0036] Compared with the prior art, the present invention has the following advantages: In arc simulation, many schemes ignore the radiation reabsorption process. Although this treatment method simplifies the calculation, it deviates greatly from the actual situation. This scheme is designed for arc simulation that takes into account the radiation reabsorption process, and is advanced and accurate. The scheme covers almost all distortion situations that may occur in arc simulation through ten types of distorted temperature distributions caused by the reabsorption phenomenon in arc simulation, and is comprehensive and universal. In addition to 0.83Tmax and 5000K, the third key threshold for judging the radiation domain and absorption domain of gas in the arc is 13000K, and all situations are divided into radiation domain / absorption domain based on these three thresholds. The method of taking the average value within the range and then comparing before and after eliminates the influence of extreme values ​​caused by small-scale fluctuations on the overall judgment. The ten radial temperature distribution curves are illustrated in a visual way, which is helpful for reading and understanding, and is easy to understand. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] Various other advantages and benefits of the present invention will become apparent to those skilled in the art upon reading the detailed description of the preferred embodiments below. The accompanying drawings are intended only to illustrate preferred embodiments and are not to be construed as limiting the present invention. It should be understood that the drawings described below are merely examples of the present invention, and that those skilled in the art will be able to derive other drawings from these drawings without inventive effort. Throughout the drawings, identical reference numerals are used to denote identical components.

[0038] In the attached figure:

[0039] Figure 1 This is a schematic diagram of finite element unit division of an arc longitudinal cross section according to an arc finite element modeling method based on a reabsorption process provided by an embodiment of the present disclosure;

[0040] Figure 2 This is a schematic diagram of determining the radiation state / absorption state based on the arc radial temperature distribution in an arc finite element modeling method based on the reabsorption process provided by one embodiment of the present disclosure;

[0041] Figure 3 This is a flow chart of an arc finite element modeling method based on a reabsorption process provided by an embodiment of the present disclosure;

[0042] Figure 4 This is a schematic diagram of an application example provided by an embodiment of the present disclosure.

[0043] The present invention will be further explained below with reference to the accompanying drawings and embodiments. DETAILED DESCRIPTION

[0044] Specific embodiments of the present invention will be described in more detail below with reference to the accompanying drawings. Although specific embodiments of the present invention are shown in the accompanying drawings, it should be understood that the present invention can be implemented in various forms and should not be limited by the embodiments set forth herein. Rather, these embodiments are provided to enable a more thorough understanding of the present invention and to fully convey the scope of the present invention to those skilled in the art.

[0045] It should be noted that certain words are used in the specification and claims to refer to specific components. Those skilled in the art should understand that technicians may use different nouns to refer to the same component. This specification and claims do not use the difference in nouns as a way to distinguish components, but use the difference in the functions of the components as the criterion for distinction. As mentioned throughout the specification and claims, "including" or "comprising" is an open term, so it should be interpreted as "including but not limited to". The subsequent description of the specification is a preferred embodiment of the present invention, but the description is based on the general principles of the specification and is not intended to limit the scope of the invention. The scope of protection of the present invention shall be as defined in the attached claims.

[0046] To facilitate understanding of the embodiments of the present invention, further explanation will be given below using specific embodiments as examples in conjunction with the accompanying drawings, and the accompanying drawings do not constitute a limitation on the embodiments of the present invention.

[0047] like Figures 1 to 4 As shown in FIG, the arc finite element modeling method based on the reabsorption process includes the following steps:

[0048] Step 1: Current passes through the gas to generate an arc, and the temperature distribution data of the arc is monitored in real time;

[0049] Step 2: Perform finite element modeling on the arc. Draw a two-dimensional longitudinal plane segmentation diagram with the arc center direction as the x-axis and the y-axis perpendicular to the arc center direction. Divide the interior of the arc into multiple triangular mesh units. The unit attributes include the coordinates of the three nodes, the geometric center coordinates of the unit, and the average temperature of the unit.

[0050] Step 3: Sort the radial units of the arc at the predetermined position, and judge each unit. If the predetermined position is between the maximum horizontal coordinate x of the three nodes of a unit, max and the minimum value of the horizontal coordinate x min If the intersection line passes through the node, write the unit object into the matrix Slab[ ]. If the intersection line passes through the node, write all the units associated with the node. Sort all the units in the matrix Slab[ ] from small to large according to the vertical coordinate value of the geometric center.

[0051] Step 4, determining whether the unit is an absorbing state region of the radiation absorbing part or a radiating state region of the radiation emitting part based on the sorting and temperature distribution data;

[0052] Step 5: Set the absorbed heat in the absorption state region and the radiated heat in the radiation state region to correct the arc simulation results.

[0053] In a preferred embodiment of the arc finite element modeling method based on the reabsorption process, step 4 includes:

[0054] Starting from the first one, calculate the average temperature of every five cells in the Slab[ ] matrix, and store them in the matrix Temp[ ] in sequence. At the same time, calculate the average vertical coordinates of the geometric centers of these five cells, and store them in the corresponding Yaver[ ] in sequence. The first value in the Slab[ ] matrix is ​​the maximum value in the matrix.

[0055] Take the Temp[ ] matrix, and compare each element Temp[i] with the two elements after it, Temp[i+1] and Temp[i+2]. If it is found that Temp[i] and Temp[i+2] are both less than Temp[i+1], record the sequence number i+1 of Temp[i+1] and write it into the matrix Index[ ]. The Slab matrix is ​​a row or column vector composed of all the related units that the arc cross section passes through, which can be called the arc cross section unit matrix.

[0056] In a preferred embodiment of the arc finite element modeling method based on the reabsorption process, if Index[ ] is empty, the arc is judged to present a spatial temperature at a predetermined position that decreases monotonically with increasing distance from the arc center, the temperature distribution space with a temperature not lower than 83% of the maximum arc temperature Tmax is the radiation state region of the radiation part during the reabsorption process, the temperature distribution space with a temperature higher than 5000K and lower than 83% of the maximum arc temperature Tmax is the absorption state region of the radiation part during the reabsorption process, and the temperature distribution space with a temperature lower than 5000K is determined not to be part of the arc.

[0057] In a preferred embodiment of the arc finite element modeling method based on the reabsorption process, if Index[ ] has 1 element, the initial values ​​of counter 1 and counter 2 are set to 0, and the temperature values ​​of every two units in the Slab[ ] matrix are taken starting from the first unit. If the 83% arc maximum temperature Tmax is between these two values, counter 1 is incremented by 1; if 13000K is between these two values, counter 2 is incremented by 1. Each time the counter changes, the two values ​​in Yaver[ ] corresponding to the two Temp[ ] are averaged and stored in the matrix Y[ ]. After the traversal is completed, the values ​​of the counters are observed.

[0058] If counter 1 of A = 3 and counter 2 = 1, the arc is judged to have a peak at a predetermined position and not at the arc center. The temperature at the arc center is higher than 83% of the arc's maximum temperature Tmax. When the temperatures before and after the peak are monotonically distributed, the temperature distribution space from the arc center to where the temperature begins to drop and drops to 83% of the arc's maximum temperature Tmax is the radiation state region of the radiation portion emitted during the reabsorption process. The temperature distribution space from the first temperature drop to 83% of the arc's maximum temperature Tmax to the end of the temperature drop to 5000K is the absorption state region of the radiation portion absorbed during the reabsorption process. The remaining temperature distribution space is judged not to be part of the arc.

[0059] If counter 1 = 3 and counter 2 = 3, then the arc is judged to have two peaks at the predetermined position, one of which is at the arc center, and there is a valley position between the two peaks. The temperature at the arc center is higher than 83% of the maximum arc temperature Tmax, the temperature peak outside the arc center is higher than 83% of the maximum arc temperature Tmax, the valley temperature is greater than 5000K and less than 13000K. When the temperatures between the extreme values ​​are monotonically distributed, the temperature distribution space from the arc center to the first time the temperature rises to 13000K is the absorption state region of the radiation part absorbed during the reabsorption process. The temperature distribution space from the first time the temperature rises to 13000K to the temperature drops to 83% of the maximum arc temperature Tmax again is the radiation state region of the radiation part emitted during the reabsorption process. The remaining temperature distribution space is judged not to be part of the arc.

[0060] If counter 1 = 1 and counter 2 = 3, it is judged that the arc has two peaks at the predetermined position, one of which is at the arc center, and there is a valley position between the two peaks. The temperature at the arc center is higher than 83% of the maximum arc temperature Tmax, the non-arc center temperature peak is lower than 83% of the maximum arc temperature Tmax and is greater than 13000K, and the valley temperature is greater than 5000K and less than 13000K. When the temperatures between the extreme values ​​are monotonically distributed, the temperature distribution space from the arc center to the first time the temperature rises to 13000K is the absorption state area of ​​the radiation part absorbed during the reabsorption process, and the temperature distribution space from the first time the temperature rises to 13000K to the temperature drops to 13000K again is the radiation state area of ​​the radiation part emitted during the reabsorption process. The remaining temperature distribution space is determined not to be part of the arc; if none of the above situations occurs.

[0061] In a preferred embodiment of the arc finite element modeling method based on the reabsorption process, if Index has two or more values, it is judged that the arc has three peaks at a predetermined position and one is at the arc center, there are two valley positions between the three peaks, the temperature at the arc center is 83% of the arc maximum temperature Tmax, the inner non-arc center temperature peak is less than 83% of the arc maximum temperature Tmax and greater than 13000K, the outer non-arc center temperature peak is greater than 5000K and less than 13000K, the inner valley temperature is greater than 13000K and less than 83% of the arc maximum temperature Tmax, the outer valley temperature is greater than 5000K and less than 13000K, and when the temperatures between the extreme values ​​are monotonically distributed, the temperature distribution space from the arc center to the second temperature peak is the absorption state area of ​​the radiation absorbing part during the reabsorption process, and the temperature distribution space outside the second temperature peak is the radiation state area of ​​the radiation emitting part during the reabsorption process.

[0062] In a preferred embodiment of the arc finite element modeling method based on the reabsorption process, step 4 includes:

[0063] Starting from the first one, calculate the average temperature of every five cells in the Slab[ ] matrix, and store them in the matrix Temp[ ] in sequence. At the same time, calculate the average vertical coordinates of the geometric centers of these five cells, and store them in the corresponding Yaver[ ] in sequence. The first value in the Slab[ ] matrix is ​​not the maximum value in the matrix.

[0064] Set the initial value of counter 1, counter 2 and counter 3 to 0, and take the temperature value of every two cells in the Slab[ ] matrix starting from the first cell. If the 83% arc maximum temperature Tmax is between these two values, counter 1 is incremented by 1; if 13000K is between these two values, counter 2 is incremented by 1; if 5000K is between these two values, counter 3 is incremented by 1. Each time the counter changes, the two values ​​in Yaver[ ] corresponding to the two Temp[ ] are averaged and stored in the matrix Y[ ]. After the traversal is completed, observe the value of the counter:

[0065] If counter 1 = 2, counter 2 = 2, and counter 3 = 2, then the arc is judged to have a peak at a predetermined position and not at the arc center, the temperature at the arc center is lower than 5000K, the non-arc center temperature peak is greater than 83% of the arc's maximum temperature Tmax, and the temperatures before and after the peak are monotonically distributed, then the temperature distribution space from the arc center to 13000K is the absorption state region of the radiation absorption part during the reabsorption process, and the temperature distribution space from the first time the temperature exceeds 13000K to the end of the temperature dropping back to 83% of the arc's maximum temperature Tmax is the radiation state region of the radiation emission part during the reabsorption process, and the remaining temperature distribution space is judged not to be part of the arc;

[0066] If counter 1 = 2, counter 2 = 2, and counter 3 = 1, then the arc is judged to have a peak at a predetermined position and not at the arc center, the temperature at the arc center is lower than 13000K but greater than 5000K, the non-arc center temperature peak is greater than 83% of the arc's maximum temperature Tmax, and the temperatures before and after the peak are monotonically distributed, then the temperature distribution space from the arc center to 13000K is the absorption state region of the radiation absorption part during the reabsorption process, and the temperature distribution space from the first time the temperature exceeds 13000K to the end of the temperature dropping back to 83% of the arc's maximum temperature Tmax is the radiation state region of the radiation emission part during the reabsorption process, and the remaining temperature distribution space is judged not to be part of the arc;

[0067] If counter 1 = 2, counter 2 = 1, and counter 3 = 1, then the arc is judged to have a peak at a predetermined position and is not at the arc center, the temperature at the arc center is higher than 13000K but lower than 83% of the arc's maximum temperature Tmax, the non-arc center temperature peak is greater than 83% of the arc's maximum temperature Tmax, and when the temperatures before and after the peak are monotonically distributed, the temperature distribution space from the arc center to where the temperature begins to drop and drops to 83% of the arc's maximum temperature Tmax is the radiation state region of the radiation portion emitted during the reabsorption process, and the temperature distribution space from the first temperature drop to 83% of the arc's maximum temperature Tmax to the end of the temperature drop to 5000K is the absorption state region of the radiation portion absorbed during the reabsorption process, and the remaining temperature distribution space is judged not to be part of the arc;

[0068] If counter 1 = 3, counter 2 = 1, and counter 3 = 1, then the arc is judged to have a peak at a predetermined position and not at the arc center, the temperature at the arc center is higher than 83% of the arc's maximum temperature Tmax, and when the temperatures before and after the peak are monotonically distributed, the temperature distribution space from the arc center to where the temperature begins to drop and drops to 83% of the arc's maximum temperature Tmax is the radiation state region of the radiation portion emitted during the reabsorption process, and the temperature distribution space from the first temperature drop to 83% of the arc's maximum temperature Tmax to the end of the temperature drop to 5000K is the absorption state region of the radiation portion absorbed during the reabsorption process, and the remaining temperature distribution space is judged not to be part of the arc;

[0069] If counter 1 = 1, counter 2 = 1, and counter 3 = 1, then it is judged that the arc has two peaks at the predetermined position, one of which is at the arc center, and there is a valley position between the two peaks. The temperature at the arc center is higher than 83% of the arc maximum temperature Tmax, the temperature peak outside the arc center is higher than 83% of the arc maximum temperature Tmax, the valley temperature is greater than 13000K and lower than 83% of the arc maximum temperature Tmax, and when the temperatures between the extreme values ​​are monotonically distributed, the temperature distribution space from the arc center to the first time the temperature rises to 83% of the arc maximum temperature Tmax is the absorption state region of the radiation part absorbed during the reabsorption process, and the temperature distribution space from the first time the temperature rises to 83% of the arc maximum temperature Tmax to the time the temperature drops to 83% of the arc maximum temperature Tmax again is the radiation state region of the radiation part emitted during the reabsorption process. The remaining temperature distribution space is determined not to be part of the arc.

[0070] If counter 1 = 3, counter 2 = 3, and counter 3 = 1, then it is judged that the arc has two peaks at the predetermined position, one of which is at the arc center, and there is a valley position between the two peaks. The temperature at the arc center is higher than 83% of the maximum arc temperature Tmax, the non-arc center temperature peak is higher than 83% of the maximum arc temperature Tmax, the valley temperature is greater than 5000K and less than 13000K. When the temperatures between the extreme values ​​are monotonically distributed, the temperature distribution space from the arc center to the first time the temperature rises to 13000K is the absorption state region of the radiation part absorbed during the reabsorption process, and the temperature distribution space from the first time the temperature rises to 13000K to the temperature drops to 83% of the maximum arc temperature Tmax again is the radiation state region of the radiation part emitted during the reabsorption process. The remaining temperature distribution space is determined not to be part of the arc.

[0071] If counter 1 = 1, counter 2 = 3, and counter 3 = 1, then it is determined that the arc has two peaks at the predetermined position, one of which is at the arc center, and there is a valley position between the two peaks. The temperature at the arc center is higher than 83% of the maximum arc temperature Tmax, the non-arc center temperature peak is lower than 83% of the maximum arc temperature Tmax and is greater than 13000K, and the valley temperature is greater than 5000K and less than 13000K. When the temperatures between the extreme values ​​are monotonically distributed, the temperature distribution space from the arc center to the first time the temperature rises to 13000K is the absorption state region of the radiation absorbing part during the reabsorption process, and the temperature distribution space from the first time the temperature rises to 13000K to the temperature drops to 13000K again is the radiation state region of the radiation emitting part during the reabsorption process. The remaining temperature distribution space is determined not to be part of the arc.

[0072] If counter 1 = 2, counter 2 = 3, and counter 3 = 1, then the arc is judged to have two peaks at the predetermined position, one of which is at the arc center, and there is a valley position between the two peaks. The temperature at the arc center is lower than 83% of the arc maximum temperature Tmax and higher than 13000K. The non-arc center temperature peak is greater than 83% of the arc maximum temperature Tmax, and the valley temperature is greater than 5000K and less than 13000K. When the temperatures between the extreme values ​​are monotonically distributed, the temperature distribution space from the arc center to the first time the temperature rises to 13000K is the absorption state region of the radiation part absorbed during the reabsorption process. The temperature distribution space from the first time the temperature rises to 13000K to the temperature drops to 83% of the arc maximum temperature Tmax again is the radiation state region of the radiation part emitted during the reabsorption process. The remaining temperature distribution space is judged not to be part of the arc. If none of the above situations occurs, an error is reported.

[0073] In a preferred embodiment of the arc finite element modeling method based on the reabsorption process, integration with Fluent software is performed to avoid divergence of the simulated temperature results.

[0074] In a preferred embodiment of the arc finite element modeling method based on the reabsorption process, in step 2, a 0.5 mm-sized grid is used for plane segmentation.

[0075] In a preferred embodiment of the arc finite element modeling method based on the reabsorption process, in step 2, the predetermined position is a position 1.0 mm away from the arc center.

[0076] In a preferred embodiment of the arc finite element modeling method based on the reabsorption process, in step 2, each cross section of the longitudinal plane sectioning diagram has only one layer of grid.

[0077] In one embodiment, the arc can be viewed microscopically as a three-dimensional symmetrical cylinder with countless cross sections. In finite element analysis, a two-dimensional longitudinal plane segmentation diagram is usually drawn with the arc center as the x-axis and the direction perpendicular to the arc center as the y-axis, and the interior of the arc is divided into a finite number of grids. The grid segmentation accuracy should also determine the density of cross-section selection. For example, if a 0.5mm grid is used for plane segmentation, a cross section should be selected every 0.5mm for analysis. This can effectively ensure that each cross section has only about one layer of grid, and can show the temperature distribution on the cross section in a more gradient manner. The steps after docking with the finite element software can be divided into three steps: sorting of cross-section units, judging the temperature distribution and determining the "critical dividing point", and dividing the radiation domain / absorption domain.

[0078] 1. Sectional element sorting

[0079] To study the radial temperature distribution, we should select any arc cross section for investigation. The following figure shows how to sort the radial distribution grid on a cross section.

[0080] Figure 1 This is a schematic diagram of the finite element division of the arc longitudinal section. The unit used for division is a triangular unit, and the x-axis is the axis where the arc center is located. Take the ordering of the radial units at 1.0mm as an example:

[0081] Create an object: a cell. It has five properties, including the coordinates of three nodes, the coordinates of the geometric center of the cell, and the average temperature of the cell.

[0082] All units in the area are numbered from 1 to 6;

[0083] For each element, if 1mm is between the maximum xmax and minimum xmin of the three nodes of a certain element, the element object is written into the matrix Slab[ ]. If the section line passes through a node, all elements associated with that node should be written. The purpose of this step is to find all elements that are located on the 1mm section, which in this case should be [element 2, element 3, element 4, element 5].

[0084] Sort all cells in the matrix Slab[ ] from smallest to largest according to the ordinate value of their geometric centers. In this example, the order should be [cell 4, cell 2, cell 3, cell 5]. This completes the ordering of the cells from the inside out.

[0085] 2. Determination of temperature distribution and determination of “critical dividing points”

[0086] As mentioned earlier, in an arc with a monotonic radial temperature distribution, heat is radiated from the arc center (Tmax) to the arc where the temperature drops to 0.83Tmax, while the arc portion where the temperature drops from 0.83Tmax to 5000K absorbs the radiated heat. However, when the radial temperature distribution is non-monotonic, there is no standard to determine which portion is in the radiating or absorbing state.

[0087] To this end, after a lot of simulation practice, we have summarized a set of judgment methods that are applicable to almost all temperature distribution situations. Figure 2 The threshold of 13000K is obtained through a lot of practice. Together with the original 0.83Tmax and 5000K as thresholds, we can see that there are three important thresholds.

[0088] With the method of judging the radiation state / absorption state according to different temperature distribution conditions, the remaining key is to judge which of the above ten temperature distribution conditions the current cross-section is in at the current moment based on the finite element simulation results. The key among the keys is to find the intersection of the temperature distribution curve and the threshold line, and thus judge which situation it belongs to and find the "critical dividing point" of the radiation domain and the absorption domain.

[0089] For the Slab[ ] matrix processed in step 3, do the following:

[0090] 1) Starting from the first one, calculate the average temperature of every five cells in the Slab[ ] matrix and store them in the Temp[ ] matrix in order. At the same time, calculate the average of the vertical coordinates of the geometric centers of these five cells and store them in the corresponding Yaver[ ] in order. (If there are less than five cells at the end, use all the remaining values ​​to calculate the average.) The values ​​in Temp[ ] and Yaver[ ] are one-to-one corresponding.

[0091] 2) Determine whether the first value in the Slab[ ] matrix is ​​the maximum value in the matrix. If so, jump to 3). If not, jump to 5).

[0092] 3) Take the Temp[ ] matrix, and compare each element Temp[i] with the two elements after it, Temp[i+1] and Temp[i+2]. If it is found that Temp[i] and Temp[i+2] are both less than Temp[i+1], record the sequence number i+1 of Temp[i+1] and write it into the matrix Index[ ]. After traversing the matrix, if Index[ ] is empty, it is the first case (indicating that the temperature distribution decreases monotonically in the radial direction), and no correction is required. Continue the simulation according to the traditional method; if Index[ ] has 1 element, jump to 4); if Index has two or more values, it is the tenth case (indicating that the radial temperature distribution has multiple peaks), and record the position of the first non-arc center peak, that is, let Y[0] = Yaver[i+1] be the "critical dividing point".

[0093] 4) Set the initial values ​​of Counter 1 and Counter 2 to 0. Starting from the first cell, take the temperature value of every two cells in the Slab[ ] matrix. If 0.83Tmax is between these two values, increment Counter 1 by 1; if 13000K is between these two values, increment Counter 2 by 1. Each time the counter changes, average the two values ​​in Yaver[ ] corresponding to the two Temp[ ] values ​​and store them in the matrix Y[ ]. After the traversal is complete, observe the counter values: If Counter 1 = 3 and Counter 2 = 1, then it is case 5; if Counter 1 = 3 and Counter 2 = 3, then it is case 7; if Counter 1 = 1 and Counter 2 = 3, then it is case 8 (this method determines which case it is by checking the number of intersections with the threshold, which will be used later). If it is not any of the above cases, report an error. The judgment ends.

[0094]

[0095] 5) Set the initial values ​​of counter 1, counter 2, and counter 3 to 0. Starting from the first cell, take the temperature values ​​of every two cells in the Slab[ ] matrix. If 0.83Tmax is between these two values, counter 1 increases by 1; if 13000K is between these two values, counter 2 increases by 1; if 5000K is between these two values, counter 3 increases by 1. Each time the counter changes, the two values ​​in Yaver[ ] corresponding to the two Temp[ ] values ​​must be averaged and stored in the matrix Y[ ]. After the traversal is complete, observe the counter values: If Counter 1 = 2, Counter 2 = 2, and Counter 3 = 2, then it is case 2; if Counter 1 = 2, Counter 2 = 2, and Counter 3 = 1, then it is case 3; if Counter 1 = 2, Counter 2 = 1, and Counter 3 = 1, then it is case 4; if Counter 1 = 3, Counter 2 = 1, and Counter 3 = 1, then it is case 5; if Counter 1 = 1, Counter 2 = 1, and Counter 3 = 1, then it is case 6; if Counter 1 = 3, Counter 2 = 3, and Counter 3 = 1, then it is case 7; if Counter 1 = 1, Counter 2 = 3, and Counter 3 = 1, then it is case 8; if Counter 1 = 2, Counter 2 = 3, and Counter 3 = 1, then it is case 9. If none of the above cases apply, report an error. The judgment ends.

[0096]

[0097] 3. Radiation domain / absorption domain division

[0098] Observe according to the obtained situation number Figure 2In the corresponding situation diagram, the key dividing points can be found in the Y[ ] matrix. For example, in the second case, the "key dividing points" are the second and fourth intersections of the curve and the threshold line. The vertical coordinates of these two intersections have been saved in Y[1] and Y[3] (the first value in C language is Y[0]). In this way, in the simulation, the correction effect can be achieved by setting all units with geometric center vertical coordinates between 0 and Y[1] to the absorption state and all units with geometric center vertical coordinates between Y[1] and Y[3] to the radiation state.

[0099] In one embodiment, 0.83Tmax is determined based on the optical thickness of the arc, 5000K is determined based on the conductive boundary of the arc, and 13000K is determined based on experience.

[0100] In one embodiment, from Figure 4 It can be clearly observed that: Figure 4 At the moment shown in the upper half of the figure, the arc's radial temperature distribution is non-monotonous. This analysis reveals that it belongs to Case 5. This distribution has a peak that is not at the arc center, where the temperature is above 83% of Tmax. The temperature distribution before and after the peak is monotonous. At this point, we determine that the arc portion from the arc center to where the temperature begins to drop and reaches 83% of Tmax is in the radiating state, while the portion from the first temperature drop to 83% of Tmax to the end of the temperature drop to 5000K is in the absorbing state. The remaining portion is not considered part of the arc.

[0101] Although the embodiments of the present invention have been described above with reference to the accompanying drawings, the present invention is not limited to the above-mentioned specific embodiments and application fields. The above-mentioned specific embodiments are merely illustrative and instructive, and are not restrictive. A person skilled in the art, guided by this specification and without departing from the scope of protection of the claims of the present invention, may also devise various forms, all of which fall within the scope of protection of the present invention.

Claims

1. An arc finite element modeling method based on the reabsorption process, characterized in that: The steps include: Step 1: Current passes through the gas to generate an arc, and the temperature distribution data of the arc is monitored in real time; Step 2: Perform finite element modeling on the arc. Draw a two-dimensional longitudinal plane segmentation diagram with the arc center direction as the x-axis and the y-axis perpendicular to the arc center direction. Divide the interior of the arc into multiple triangular mesh units. The unit attributes include the coordinates of the three nodes, the geometric center coordinates of the unit, and the average temperature of the unit. Step 3: Sort the radial units of the arc at the predetermined position, and judge each unit. If the predetermined position is between the maximum horizontal coordinate x of the three nodes of a unit, max and the minimum value of the horizontal axis x min If the intersection line passes through the node, write the unit object into the matrix Slab[ ]. If the intersection line passes through the node, write all the units associated with the node. Sort all the units in the matrix Slab[ ] from small to large according to the vertical coordinate value of the geometric center. Step 4, determining whether the unit is an absorbing state region of the radiation absorbing part or a radiating state region of the radiation emitting part based on the sorting and temperature distribution data; Step 5: Set the absorbed heat in the absorption state region and the radiated heat in the radiation state region to correct the arc simulation results.

2. The arc finite element modeling method based on the reabsorption process according to claim 1, characterized in that: Step 4 includes, Starting from the first one, calculate the average temperature of every five cells in the Slab[ ] matrix, and store them in the matrix Temp[ ] in sequence. At the same time, calculate the average vertical coordinates of the geometric centers of these five cells, and store them in the corresponding Yaver[ ] in sequence. The first value in the Slab[ ] matrix is ​​the maximum value in the matrix. Take the Temp[ ] matrix, and compare each element Temp[i] with the two elements after it, Temp[i+1] and Temp[i+2]. If it is found that Temp[i] and Temp[i+2] are both smaller than Temp[i+1], record the serial number i+1 of Temp[i+1] and write it into the matrix Index[ ].

3. The arc finite element modeling method based on the reabsorption process according to claim 2, characterized in that: If Index[ ] is empty, the arc is judged to have a spatial temperature at the predetermined position that decreases monotonically with increasing distance from the arc center. The temperature distribution space with a temperature not lower than 83% of the arc's maximum temperature Tmax is the radiation state region of the radiation portion emitted during the reabsorption process. The temperature distribution space with a temperature higher than 5000K and lower than 83% of the arc's maximum temperature Tmax is the absorption state region of the radiation portion absorbed during the reabsorption process. The temperature distribution space with a temperature lower than 5000K is judged not to be part of the arc.

4. The arc finite element modeling method based on the reabsorption process according to claim 2, characterized in that: If Index[ ] has 1 element, set the initial value of counter 1 and counter 2 to 0, and take the temperature value of every two cells in the Slab[ ] matrix starting from the first cell. If the 83% arc maximum temperature Tmax is between these two values, counter 1 is incremented by 1; if 13000K is between these two values, counter 2 is incremented by 1. Each time the counter changes, the two values ​​in Yaver[ ] corresponding to the two Temp[ ] are averaged and stored in the matrix Y[ ]. After the traversal is completed, observe the value of the counter. If counter 1 = 3 and counter 2 = 1, the arc is judged to have a peak at a predetermined position and not at the arc center. The temperature at the arc center is higher than 83% of the arc's maximum temperature Tmax. When the temperatures before and after the peak are monotonically distributed, the temperature distribution space from the arc center to where the temperature begins to drop and drops to 83% of the arc's maximum temperature Tmax is the radiation state region of the radiation portion emitted during the reabsorption process. The temperature distribution space from the first temperature drop to 83% of the arc's maximum temperature Tmax to the end of the temperature drop to 5000K is the absorption state region of the radiation portion absorbed during the reabsorption process. The remaining temperature distribution space is judged not to be part of the arc. If counter 1 = 3 and counter 2 = 3, then the arc is judged to have two peaks at the predetermined position, one of which is at the arc center, and there is a valley position between the two peaks. The temperature at the arc center is higher than 83% of the maximum arc temperature Tmax, the temperature peak outside the arc center is higher than 83% of the maximum arc temperature Tmax, the valley temperature is greater than 5000K and less than 13000K. When the temperatures between the extreme values ​​are monotonically distributed, the temperature distribution space from the arc center to the first time the temperature rises to 13000K is the absorption state region of the radiation part absorbed during the reabsorption process. The temperature distribution space from the first time the temperature rises to 13000K to the temperature drops to 83% of the maximum arc temperature Tmax again is the radiation state region of the radiation part emitted during the reabsorption process. The remaining temperature distribution space is judged not to be part of the arc. If counter 1 = 1 and counter 2 = 3, the arc is judged to have two peaks at the predetermined position, one of which is at the arc center, and there is a valley position between the two peaks. The temperature at the arc center is higher than 83% of the arc maximum temperature Tmax, the non-arc center temperature peak is lower than 83% of the arc maximum temperature Tmax and greater than 13000K, and the valley temperature is greater than 5000K and less than 13000K. When the temperatures between the extreme values ​​are monotonically distributed, the temperature distribution space from the arc center to the first time the temperature rises to 13000K is the absorption state region of the radiation part absorbed during the reabsorption process, and the temperature distribution space from the first time the temperature rises to 13000K to the temperature drops to 13000K again is the radiation state region of the radiation part emitted during the reabsorption process. The remaining temperature distribution space is determined not to be part of the arc. If none of the above situations occurs, an error is reported.

5. The arc finite element modeling method based on the reabsorption process according to claim 2, characterized in that: If Index has two or more values, it is judged that the arc has three peaks at the predetermined position and one is at the arc center, there are two valley positions between the three peaks, the temperature at the arc center is 83% of the arc maximum temperature Tmax, the inner non-arc center temperature peak is less than 83% of the arc maximum temperature Tmax and greater than 13000K, the outer non-arc center temperature peak is greater than 5000K and less than 13000K, the inner valley temperature is greater than 13000K and less than 83% of the arc maximum temperature Tmax, the outer valley temperature is greater than 5000K and less than 13000K, and when the temperatures between the extreme values ​​are monotonically distributed, the temperature distribution space from the arc center to the second temperature peak is the absorption state region of the radiation part absorbed during the reabsorption process, and the temperature distribution space outside the second temperature peak is the radiation state region of the radiation part emitted during the reabsorption process.

6. The arc finite element modeling method based on the reabsorption process according to claim 1, characterized in that: Step 4 includes, Starting from the first one, calculate the average temperature of every five cells in the Slab[ ] matrix, and store them in the matrix Temp[ ] in sequence. At the same time, calculate the average vertical coordinates of the geometric centers of these five cells, and store them in the corresponding Yaver[ ] in sequence. The first value in the Slab[ ] matrix is ​​not the maximum value in the matrix. Set the initial value of counter 1, counter 2 and counter 3 to 0, and take the temperature value of every two cells in the Slab[ ] matrix starting from the first cell. If the 83% arc maximum temperature Tmax is between these two values, counter 1 is incremented by 1; if 13000K is between these two values, counter 2 is incremented by 1; if 5000K is between these two values, counter 3 is incremented by 1. Each time the counter changes, the two values ​​in Yaver[ ] corresponding to the two Temp[ ] are averaged and stored in the matrix Y[ ]. After the traversal is completed, observe the value of the counter: If counter 1 = 2, counter 2 = 2, and counter 3 = 2, then the arc is judged to have a peak at a predetermined position and not at the arc center, the temperature at the arc center is lower than 5000K, the non-arc center temperature peak is greater than 83% of the arc's maximum temperature Tmax, and the temperatures before and after the peak are monotonically distributed, then the temperature distribution space from the arc center to 13000K is the absorption state region of the radiation absorption part during the reabsorption process, and the temperature distribution space from the first time the temperature exceeds 13000K to the end of the temperature dropping back to 83% of the arc's maximum temperature Tmax is the radiation state region of the radiation emission part during the reabsorption process, and the remaining temperature distribution space is judged not to be part of the arc; If counter 1 = 2, counter 2 = 2, and counter 3 = 1, then the arc is judged to have a peak at a predetermined position and not at the arc center, the temperature at the arc center is lower than 13000K but greater than 5000K, the non-arc center temperature peak is greater than 83% of the arc's maximum temperature Tmax, and the temperatures before and after the peak are monotonically distributed, then the temperature distribution space from the arc center to 13000K is the absorption state region of the radiation absorption part during the reabsorption process, and the temperature distribution space from the first time the temperature exceeds 13000K to the end of the temperature dropping back to 83% of the arc's maximum temperature Tmax is the radiation state region of the radiation emission part during the reabsorption process, and the remaining temperature distribution space is judged not to be part of the arc; If counter 1 = 2, counter 2 = 1, and counter 3 = 1, then the arc is judged to have a peak at a predetermined position and is not at the arc center, the temperature at the arc center is higher than 13000K but lower than 83% of the arc's maximum temperature Tmax, the non-arc center temperature peak is greater than 83% of the arc's maximum temperature Tmax, and when the temperatures before and after the peak are monotonically distributed, the temperature distribution space from the arc center to where the temperature begins to drop and drops to 83% of the arc's maximum temperature Tmax is the radiation state region of the radiation portion emitted during the reabsorption process, and the temperature distribution space from the first temperature drop to 83% of the arc's maximum temperature Tmax to the end of the temperature drop to 5000K is the absorption state region of the radiation portion absorbed during the reabsorption process, and the remaining temperature distribution space is judged not to be part of the arc; If counter 1 = 3, counter 2 = 1, and counter 3 = 1, then the arc is judged to have a peak at a predetermined position and not at the arc center, the temperature at the arc center is higher than 83% of the arc's maximum temperature Tmax, and when the temperatures before and after the peak are monotonically distributed, the temperature distribution space from the arc center to where the temperature begins to drop and drops to 83% of the arc's maximum temperature Tmax is the radiation state region of the radiation portion emitted during the reabsorption process, and the temperature distribution space from the first temperature drop to 83% of the arc's maximum temperature Tmax to the end of the temperature drop to 5000K is the absorption state region of the radiation portion absorbed during the reabsorption process, and the remaining temperature distribution space is judged not to be part of the arc; If counter 1 = 1, counter 2 = 1, and counter 3 = 1, then it is judged that the arc has two peaks at the predetermined position, one of which is at the arc center, and there is a valley position between the two peaks. The temperature at the arc center is higher than 83% of the arc maximum temperature Tmax, the temperature peak outside the arc center is higher than 83% of the arc maximum temperature Tmax, the valley temperature is greater than 13000K and lower than 83% of the arc maximum temperature Tmax, and when the temperatures between the extreme values ​​are monotonically distributed, the temperature distribution space from the arc center to the first time the temperature rises to 83% of the arc maximum temperature Tmax is the absorption state region of the radiation part absorbed during the reabsorption process, and the temperature distribution space from the first time the temperature rises to 83% of the arc maximum temperature Tmax to the time the temperature drops to 83% of the arc maximum temperature Tmax again is the radiation state region of the radiation part emitted during the reabsorption process. The remaining temperature distribution space is determined not to be part of the arc. If counter 1 = 3, counter 2 = 3, and counter 3 = 1, then it is judged that the arc has two peaks at the predetermined position, one of which is at the arc center, and there is a valley position between the two peaks. The temperature at the arc center is higher than 83% of the maximum arc temperature Tmax, the non-arc center temperature peak is higher than 83% of the maximum arc temperature Tmax, the valley temperature is greater than 5000K and less than 13000K. When the temperatures between the extreme values ​​are monotonically distributed, the temperature distribution space from the arc center to the first time the temperature rises to 13000K is the absorption state region of the radiation part absorbed during the reabsorption process, and the temperature distribution space from the first time the temperature rises to 13000K to the temperature drops to 83% of the maximum arc temperature Tmax again is the radiation state region of the radiation part emitted during the reabsorption process. The remaining temperature distribution space is determined not to be part of the arc. If counter 1 = 1, counter 2 = 3, and counter 3 = 1, then it is determined that the arc has two peaks at the predetermined position, one of which is at the arc center, and there is a valley position between the two peaks. The temperature at the arc center is higher than 83% of the maximum arc temperature Tmax, the non-arc center temperature peak is lower than 83% of the maximum arc temperature Tmax and is greater than 13000K, and the valley temperature is greater than 5000K and less than 13000K. When the temperatures between the extreme values ​​are monotonically distributed, the temperature distribution space from the arc center to the first time the temperature rises to 13000K is the absorption state region of the radiation absorbing part during the reabsorption process, and the temperature distribution space from the first time the temperature rises to 13000K to the temperature drops to 13000K again is the radiation state region of the radiation emitting part during the reabsorption process. The remaining temperature distribution space is determined not to be part of the arc. If counter 1 = 2, counter 2 = 3, and counter 3 = 1, then the arc is judged to have two peaks at the predetermined position, one of which is at the arc center, and there is a valley position between the two peaks. The temperature at the arc center is lower than 83% of the arc maximum temperature Tmax and higher than 13000K. The non-arc center temperature peak is greater than 83% of the arc maximum temperature Tmax, and the valley temperature is greater than 5000K and less than 13000K. When the temperatures between the extreme values ​​are monotonically distributed, the temperature distribution space from the arc center to the first time the temperature rises to 13000K is the absorption state region of the radiation part absorbed during the reabsorption process. The temperature distribution space from the first time the temperature rises to 13000K to the temperature drops to 83% of the arc maximum temperature Tmax again is the radiation state region of the radiation part emitted during the reabsorption process. The remaining temperature distribution space is judged not to be part of the arc. If none of the above situations occurs, an error is reported.

7. The arc finite element modeling method based on the reabsorption process according to claim 1, characterized in that: Integration with Fluent software avoids divergence of simulation temperature results.

8. The arc finite element modeling method based on the reabsorption process according to claim 1, characterized in that: In step 2, a 0.5 mm mesh is used for plane segmentation.

9. The arc finite element modeling method based on the reabsorption process according to claim 8, characterized in that: In step 2, the predetermined position is a position 1.0 mm away from the arc center.

10. The arc finite element modeling method based on the reabsorption process according to claim 1, characterized in that: In step 2, each cross section of the longitudinal plane segmentation diagram has only one layer of mesh.

Images