CN113673124A - Numerical simulation prediction method, system and medium for three-way intersection line welding temperature field - Google Patents
Numerical simulation prediction method, system and medium for three-way intersection line welding temperature field Download PDFInfo
- Publication number
- CN113673124A CN113673124A CN202110761134.5A CN202110761134A CN113673124A CN 113673124 A CN113673124 A CN 113673124A CN 202110761134 A CN202110761134 A CN 202110761134A CN 113673124 A CN113673124 A CN 113673124A
- Authority
- CN
- China
- Prior art keywords
- welding
- heat source
- model
- temperature field
- coordinate system
- 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.)
- Pending
Links
- 238000003466 welding Methods 0.000 title claims abstract description 181
- 238000000034 method Methods 0.000 title claims abstract description 71
- 238000004088 simulation Methods 0.000 title claims abstract description 49
- 230000008569 process Effects 0.000 claims abstract description 28
- 239000011324 bead Substances 0.000 claims abstract description 20
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 239000013598 vector Substances 0.000 claims description 20
- 230000008859 change Effects 0.000 claims description 14
- 239000000463 material Substances 0.000 claims description 14
- 238000012546 transfer Methods 0.000 claims description 13
- 238000005315 distribution function Methods 0.000 claims description 10
- 239000007787 solid Substances 0.000 claims description 9
- 230000006870 function Effects 0.000 claims description 8
- 238000010408 sweeping Methods 0.000 claims description 7
- 238000009792 diffusion process Methods 0.000 claims description 5
- 230000036314 physical performance Effects 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 230000001788 irregular Effects 0.000 abstract description 3
- 230000007704 transition Effects 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 4
- 238000001816 cooling Methods 0.000 description 4
- 230000000704 physical effect Effects 0.000 description 4
- 229910000831 Steel Inorganic materials 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 239000010959 steel Substances 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 239000011572 manganese Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- PWHULOQIROXLJO-UHFFFAOYSA-N Manganese Chemical compound [Mn] PWHULOQIROXLJO-UHFFFAOYSA-N 0.000 description 1
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 description 1
- 239000000956 alloy Substances 0.000 description 1
- 229910045601 alloy Inorganic materials 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 229910052748 manganese Inorganic materials 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 229910052698 phosphorus Inorganic materials 0.000 description 1
- 239000011574 phosphorus Substances 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 229910052717 sulfur Inorganic materials 0.000 description 1
- 239000011593 sulfur Substances 0.000 description 1
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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
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)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
The invention discloses a numerical simulation prediction method, a numerical simulation prediction system and a numerical simulation prediction medium for a three-way intersection line welding temperature field, wherein the method comprises the following steps of: establishing a three-dimensional model, and performing grid division on the three-dimensional model, wherein the three-dimensional model comprises a three-way pipe fitting and a welding seam; establishing a welded heat source model, and acquiring a moving track curve of the heat source model according to the grid node coordinates of a welding bead in the three-dimensional model; and (4) performing simulation analysis by combining the heat source model, the track curve and the three-dimensional model to obtain temperature parameters on the welding component in the welding process. The method can realize the analog calculation of the temperature field in the three-way intersection line welding process, and in addition, the welding track curve is fitted by capturing the coordinates of the nodes of the welding bead grid, so that the method is suitable for any irregular space welding curve, and the problem that the welding track is difficult to accurately describe through the traditional geometric parameter equation is solved. The invention can be widely applied to the field of welding numerical simulation.
Description
Technical Field
The invention relates to the field of welding numerical simulation, in particular to a numerical simulation prediction method, a numerical simulation prediction system and a numerical simulation prediction medium for a three-way intersection line welding temperature field.
Background
In the field of petrochemical industry, a pressure pipeline is an important medium for heat transfer and mass transfer, is widely applied, is mostly a high-temperature, inflammable and explosive medium for transportation, and can cause serious casualties and huge economic loss once being damaged. The pressure pipeline often needs the tee bend pipe fitting to change the way or shunt in laying process, has a large amount of pipe looks transversal weld joint welding problems. The intersecting line welding seam belongs to complex space curve welding, the process is complex, the welding seam is difficult to position, and the difficulty coefficient of ensuring the quality is large. In addition, the welding process has the characteristics of high temperature, instantaneity, dynamics and the like, and the change condition of parameters such as the temperature on a welding component in the welding process is difficult to accurately obtain by using the traditional testing method.
Disclosure of Invention
In order to solve at least one of the technical problems in the prior art to a certain extent, the invention aims to provide a numerical simulation prediction method, a numerical simulation prediction system and a numerical simulation prediction medium for a three-way intersection line welding temperature field.
The technical scheme adopted by the invention is as follows:
a numerical simulation prediction method for a three-way intersection line welding temperature field comprises the following steps:
establishing a three-dimensional model, and performing grid division on the three-dimensional model, wherein the three-dimensional model comprises a three-way pipe fitting and a welding seam;
establishing a welded heat source model, and acquiring a moving track curve of the heat source model according to the grid node coordinates of a welding bead in the three-dimensional model;
and (4) performing simulation analysis by combining the heat source model, the track curve and the three-dimensional model to obtain temperature parameters on the welding component in the welding process.
Further, the establishing of the three-dimensional model and the meshing of the three-dimensional model include:
establishing three-dimensional solid models of the three-way pipe fitting and the welding seam, and combining the three-dimensional solid models of the three-way pipe fitting and the welding seam into an integral geometric model;
meshing the geometric model according to the sequence of diffusion from the welding seam to the periphery;
predefining physical performance parameters of the pipe fitting and the welding material along with temperature change, and respectively endowing material properties to corresponding structural parts.
Further, the meshing the geometric model in the order of diffusion from the weld to the surroundings includes:
segmenting an integral model of the geometric model into a plurality of entities;
2D meshing the surface of the entity, and sweeping the generated volumetric mesh.
Further, the establishing of the welded heat source model and the obtaining of the moving trajectory curve of the heat source model according to the grid node coordinates of the weld bead in the three-dimensional model include:
acquiring grid node coordinates of the heat source model moving in each welding pass, and fitting the grid node coordinates to obtain a path of the heat source model moving along with time;
based on a double-ellipsoid heat source model, aiming at the characteristic that the welding of the three-way intersection line has a complex spatial path, establishing a local coordinate system for the heat source model;
and establishing a heat source model for three-way intersecting line welding according to a local coordinate system.
Further, the establishing a local coordinate system for the heat source model includes:
setting the rectangular coordinate system of the three-dimensional model as the global coordinate system x, y and z and the local coordinate system x of the heat source modeli、yi、zi(ii) a Local coordinate system xiCoincidence of tangent vectors of axis and weld curve, ziThe axis is the normal vector of the discretization grid surface unit of the welding bead, yiAxis xiAxial coordinate vector cross-product ziAn axial coordinate vector;
the transformation relation between the global coordinate system and the local coordinate system is as follows:
wherein n isx1(ti) Corresponding coordinate axes for the heat source model under a local coordinate systemComponent of the vector in the global coordinate system, Fx(ti)、Fy(ti)、Fz(ti) Is tiAnd the central position of the heat source at the moment is a coordinate value under the global coordinate system.
Further, the establishing of the integral heat source model for the three-way intersecting line welding according to the local coordinate system comprises:
writing heat source subprogram by using Fortran programming language, t during weldingiThe central position of the heat source is (x)i,yi,zi) Establishing a heat source model for three-way intersecting line welding;
the heat source model is divided into the following two parts:
ellipsoidal internal heat source distribution function along the first half of the welding direction:
ellipsoidal internal heat source distribution function in the second half of the welding direction:
wherein q is1、q2As a function of heat source distribution; a isf、arB and c are shape parameters of a double-ellipsoid heat source; eta is welding thermal efficiency; u is welding voltage; i is welding current; f. offThe heat distribution coefficient of the front half part of the heat source model; f. ofrThe heat distribution coefficient of the latter half of the heat source model.
Further, in the process of simulation analysis, heat transfer calculation analysis is carried out by adopting a differential equation based on heat transfer control:
wherein ρ is density; c is the specific heat capacity of the material; t is a temperature distribution function; t is the heat transfer time;
before the simulation analysis, the method further comprises the following steps:
the initial temperature of the welded component and the boundary conditions of the temperature field are set.
The other technical scheme adopted by the invention is as follows:
a three-way intersection line welding temperature field numerical simulation prediction system comprises:
the welding model building module is used for building a three-dimensional model and carrying out meshing division on the three-dimensional model, wherein the three-dimensional model comprises a three-way pipe fitting and a welding seam;
the heat source model building module is used for building a welded heat source model and obtaining a moving track curve of the heat source model according to the grid node coordinates of a welding bead in the three-dimensional model;
and the temperature field prediction module is used for carrying out simulation analysis by combining the heat source model, the track curve and the three-dimensional model to obtain temperature parameters on the welding component in the welding process.
The other technical scheme adopted by the invention is as follows:
a three-way intersection line welding temperature field numerical simulation prediction system comprises:
at least one processor;
at least one memory for storing at least one program;
when executed by the at least one processor, cause the at least one processor to implement the method described above.
The other technical scheme adopted by the invention is as follows:
a storage medium having stored therein a processor-executable program for performing the method as described above when executed by a processor.
The invention has the beneficial effects that: the method can realize the analog calculation of the temperature field in the three-way intersection line welding process, and in addition, the welding track curve is fitted by capturing the coordinates of the nodes of the welding bead grid, so that the method is suitable for any irregular space welding curve, and the problem that the welding track is difficult to accurately describe through the traditional geometric parameter equation is solved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following description is made on the drawings of the embodiments of the present invention or the related technical solutions in the prior art, and it should be understood that the drawings in the following description are only for convenience and clarity of describing some embodiments in the technical solutions of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to these drawings without creative efforts.
FIG. 1 is a flow chart illustrating the steps of a method for predicting the numerical simulation of a three-way intersection welding temperature field according to an embodiment of the present invention;
FIG. 2 is a schematic structural diagram of a tee in an embodiment of the present invention;
FIG. 3 is a grid-divided view of a tee in an embodiment of the present invention;
FIG. 4 is a schematic diagram of the temperature distribution of the 1 st, 3 rd, 5 th and 7 th welds of the intersecting line of the tee in the embodiment of the invention;
FIG. 5 is a graph of the thermal cycle at each of the nine welds in the weld of the present invention.
Detailed Description
Reference will now be made in detail to embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to the same or similar elements or elements having the same or similar function throughout. The embodiments described below with reference to the accompanying drawings are illustrative only for the purpose of explaining the present invention, and are not to be construed as limiting the present invention. The step numbers in the following embodiments are provided only for convenience of illustration, the order between the steps is not limited at all, and the execution order of each step in the embodiments can be adapted according to the understanding of those skilled in the art.
In the description of the present invention, it should be understood that the orientation or positional relationship referred to in the description of the orientation, such as the upper, lower, front, rear, left, right, etc., is based on the orientation or positional relationship shown in the drawings, and is only for convenience of description and simplification of description, and does not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and thus, should not be construed as limiting the present invention.
In the description of the present invention, the meaning of a plurality of means is one or more, the meaning of a plurality of means is two or more, and larger, smaller, larger, etc. are understood as excluding the number, and larger, smaller, inner, etc. are understood as including the number. If the first and second are described for the purpose of distinguishing technical features, they are not to be understood as indicating or implying relative importance or implicitly indicating the number of technical features indicated or implicitly indicating the precedence of the technical features indicated.
In the description of the present invention, unless otherwise explicitly limited, terms such as arrangement, installation, connection and the like should be understood in a broad sense, and those skilled in the art can reasonably determine the specific meanings of the above terms in the present invention in combination with the specific contents of the technical solutions.
As shown in fig. 1, the present embodiment provides a three-way intersection line welding temperature field numerical simulation prediction method, which includes the following steps:
and S1, establishing a finite element model of the tee intersection line welding process.
Step S1 includes steps S11-S13:
s11, analyzing the simulation object, establishing a three-dimensional geometric model by utilizing SolidWorks three-dimensional modeling software, wherein the three-dimensional geometric model comprises a main pipe, branch pipes and welding line parts, combining and assembling the parts of the model into an integral geometric model, exporting the integral geometric model in an SLDPRT file format, and importing the model into Hypermesh meshing software.
S12, firstly, dividing the integral model into different solid, then carrying out 2D grid division on the welding seam section and part of surface, then sweeping the generated body grid, and realizing grid transition division by adjusting the along bias style; the grid size is adjusted according to the actual situation, in order to guarantee the calculation precision and improve the calculation efficiency, the grids of the welding line and the nearby area are finely divided, the grid size is increased in the area far away from the welding line, if the grid division fails, whether the grid size is reasonable or not and whether the sweeping direction is correct or not can be checked, after the grid division of the whole model is completed, the grid quality is checked, and when the grid quality is not wrong and the warning value is smaller than 10%, the grid is considered to be qualified.
S13, predefining physical performance parameters of the parent metal and the welding material along with temperature change by using Hypermesh, respectively giving material attributes to corresponding structural parts, and modifying specific parameters after finite element software is introduced.
And S2, setting initial conditions and boundary conditions of the welding simulation.
In an initial step, the initial temperature of the weldment is set and the boundary conditions of the temperature field (including thermal convection and thermal radiation) are set to obtain a relatively accurate temperature field.
And S3, establishing a welding heat source model.
Step S3 includes steps S31-S33:
s31, compiling a script for extracting the coordinates of the heat source mobile unit node of each welding bead by utilizing a Python programming language, and fitting the coordinate values to obtain a path of the heat source model moving along with time.
The Python script is operated to obtain the coordinates of the nodes of the heat source mobile unit of each welding bead, then Excel is used for splitting the coordinates to respectively obtain the coordinate values of the nodes along the x direction, the y direction and the z direction along with the time change, and then polynomial curve fitting is carried out to obtain a heat source model movement track equation.
S32, based on a common double-ellipsoid heat source model, aiming at the characteristic that welding of three-way intersecting lines has a complex spatial path, a local coordinate system is established for the heat source.
Setting the rectangular coordinate system of the geometric model as the global coordinate system x, y, z and the local coordinate system x of the heat sourcei、yi、zi. Local coordinate system xiCoincidence of tangent vectors of axis and weld curve, ziThe axis is the normal vector of the discretization grid surface unit of the welding bead, yiAxis xiAxial coordinate vector cross-product ziAn axial coordinate vector. The coordinate transformation relation is as follows:
wherein n isx1(ti) Wait for 9 valuesIs the component of the coordinate axis vector in the global coordinate system under the local coordinate system of the heat source, Fx(ti)、Fy(ti)、Fz(ti) Is tiAnd the central position of the heat source at the moment is a coordinate value under the global coordinate system.
S33, writing a heat source subprogram by using Fortran programming language, and t during weldingiThe central position of the heat source is (x)i,yi,zi) Establishing an integral heat source model for three-way intersecting line welding, wherein the heat source model is divided into the following two parts:
ellipsoidal internal heat source distribution function along the first half of the welding direction:
ellipsoidal internal heat source distribution function in the second half of the welding direction:
wherein q is1、q2As a function of heat source distribution; a isf、arB and c are shape parameters of a double-ellipsoid heat source; eta is welding thermal efficiency; u is welding voltage; i is welding current; f. offThe heat distribution coefficient of the front half part of the heat source model; f. ofrThe heat distribution coefficient of the latter half of the heat source model.
And S4, creating a welding and cooling process analysis step.
And compiling analysis steps of multi-pass welding and each post-welding cooling process by using a Python programming language, and then applying a moving heat source which changes along with time on the whole weldment by selecting a user-defined heat source subprogram.
And S5, submitting the task to calculate and solve and post-processing.
Step S5 includes steps S51-S52:
s51, submitting a solving task to solve the temperature field, and performing heat transfer calculation analysis by finite element software by adopting a differential equation based on heat transfer control:
when the temperature field is solved, the welding process is simulated through the live and dead units, all units of the welding line need to be killed before solving, and the units are gradually activated in the solving process, so that the whole-process simulation of the welding transient temperature field is realized.
And S52, entering a post processor, reading a calculation result, acquiring temperature parameters on a welding component in the welding process, and selecting any point in the weldment to analyze the temperature change condition.
The above method is explained in detail with reference to specific examples below.
A numerical simulation prediction method for a three-way intersection welding temperature field is described by taking the calculation of the prefabricated three-way (Q345R) intersection welding process in the opening plugging of an urban gas pipeline as an example, wherein the chemical components (wt%) of Q345R are shown in Table 1:
TABLE 1
Alloy element | Carbon (C) | Silicon (Si) | Manganese (Mn) | Sulfur (S) | Phosphorus (P) | Aluminum (Al) |
wt% | 0.13 | 0.35 | 1.39 | 0.006 | 0.011 | 0.037 |
The method is implemented according to the following steps:
step 1, establishing a finite element model of a tee intersection line welding process;
step 1.1, selecting a quarter part of a tee joint (namely taking a quarter structure according to a symmetric axis) because the model of the tee joint is symmetrical and focusing on the change condition of a temperature field near the intersecting line of a tee joint pipe fitting, establishing a three-dimensional solid model for the tee joint pipe and a welding line by utilizing SolidWorks three-dimensional modeling software, and combining the three-dimensional solid model into an integral geometric model, as shown in figure 2, wherein the outer diameter of a main pipe is 385.5mm, the length of the main pipe is 600mm, the outer diameter of a branch pipe is 355mm, the length of the branch pipe is 541mm, the wall thickness of the branch pipe is 25mm, and the angle of the welding line is 45 degrees; exporting the model in SLDPRT file format, and importing the model into Hypermesh software for meshing.
And step 1.2, carrying out mesh division on the geometric model obtained in the step 1.1 by using Hypermesh software according to the sequence of diffusion from the welding seam to the periphery.
Step 1.2 specifically, the model is divided into 5 Solid bodies, then 2D meshing is carried out on the surface of the finite element model by using a transition unit, and then the whole welding bead is subjected to meshing in a sweeping mode; the grid size is adjusted according to the actual situation, in order to ensure the calculation precision and improve the calculation efficiency, the grids of the welding line and the nearby area are finely divided, and the grid size is increased in the area far away from the welding line, in the example, the size of the welding line and the nearby grids is set to be 2mm, and the size of the rest grids is larger the farther away from the welding line, as shown in fig. 3; if the grid division fails, whether the grid size is reasonable and the sweeping direction is correct can be checked, after the grid division of the whole model is completed, the grid quality is checked, and when the grid quality has no error and the warning value is less than 10%, the grid is considered to be qualified.
In this embodiment, the 5 solids are respectively a pipe body 1 (such as a main pipe), a pipe body 2 (such as a branch pipe), a transition region 1, a transition region 2 and a weld joint region, 2D meshing is performed on the weld joint region and the sections of the two transition regions, and a body mesh is generated by sweeping along an intersecting line; and then 2D meshing is carried out on the outer walls of the main pipe and the branch pipes, and the generated body meshes are swept along the radial direction.
Step 1.3, predefining physical performance parameters of the material, including density, thermal conductivity coefficient, specific heat capacity, elastic modulus, Poisson's ratio, thermal expansion coefficient and the like, of the material along with temperature change by using Hypermesh software, respectively endowing material properties to corresponding structures, and modifying specific parameters after a model is introduced into finite element software. The physical properties of the E5015 welding material used in the example are similar to those of the Q345R, so that the physical property parameters of the welding seam are replaced by the physical property parameters of the Q345R steel, and the physical property parameters of the Q345R steel along with the change of temperature are shown in the following table 2:
TABLE 2
Step 2 specifically, setting the initial temperature of the weldment and the boundary conditions of the temperature field;
the initial temperature of the weldment and the boundary conditions of the temperature field, including thermal convection and thermal radiation, are set in an initial step to obtain a relatively accurate temperature field. In this embodiment, the initial temperature of the weldment is set to 20 ℃, the heat radiation and the heat convection are considered comprehensively, and the heat transfer coefficient is set to 10 mW/(mm)2.C)。
The step 3 specifically comprises the following steps:
3.1, compiling a script for extracting the coordinates of the heat source mobile unit node of each welding bead by using a Python programming language, and fitting the coordinate values to obtain a path of the heat source model moving along with time;
the Python script is operated to obtain the coordinates of the nodes of the heat source mobile unit of each welding bead, then Excel is used for splitting the coordinates to respectively obtain the coordinate values of the nodes along the x direction, the y direction and the z direction along with the time change, and then polynomial curve fitting is carried out to obtain a heat source model movement track parameter equation.
3.2, based on a common double-ellipsoid heat source model, aiming at the characteristic that the welding of the three-way intersection line has a complex spatial path, establishing a local coordinate system for the heat source;
setting the rectangular coordinate system of the geometric model as the global coordinate system x, y, z and the local coordinate system x of the heat sourcei、yi、zi. Local coordinate system xiCoincidence of tangent vectors of axis and weld curve, ziThe axis is the normal vector of the discretization grid surface unit of the welding bead, yiAxis xiAxial coordinate vector cross-product ziAn axial coordinate vector. The coordinate transformation relation is as follows:
wherein n isx1(ti) The 9 values are components of the coordinate axis vector corresponding to the heat source in the local coordinate system in the global coordinate system, Fx(ti)、Fy(ti)、Fz(ti) Is tiAnd the central position of the heat source at the moment is a coordinate value under the global coordinate system.
Step 3.3, writing a heat source subprogram by using Fortran programming language, and t during weldingiThe central position of the heat source is (x)i,yi,zi) Establishing an integral heat source model for three-way intersecting line welding, wherein the heat source model is divided into the following two parts:
ellipsoidal internal heat source distribution function along the first half of the welding direction:
ellipsoidal internal heat source distribution function in the second half of the welding direction:
wherein q is1、q2As a function of heat source distribution; a isf、arB and c are shape parameters of a double-ellipsoid heat source; eta is welding thermal efficiency; u is welding voltage; i is welding current; f. offThe heat distribution coefficient of the front half part of the heat source model; f. ofrThe heat distribution coefficient of the latter half of the heat source model. In this embodiment, af=15、ar=20、b=10、c=10、η=0.85、ff=0.6、fr=1.4、U=25V、I=100A;
And 4, compiling analysis steps of multi-pass welding and each post-welding cooling process by using a Python programming language, and then applying a moving heat source which changes along with time on the whole weldment by selecting a user-defined heat source subprogram.
In this embodiment, each weld is 225s, and the next weld is performed after cooling for 100s until the end.
The step 5 specifically comprises the following steps:
step 5.1, submitting a solving task to solve the temperature field, and performing heat transfer calculation analysis by finite element software by adopting a differential equation based on heat transfer control:
when the temperature field is solved, the welding process is simulated through the live and dead units, all units of the welding line need to be killed before solving, and the units are gradually activated in the solving process, so that the whole-process simulation of the welding transient temperature field is realized;
and 5.2, entering a post processor, reading a calculation result, acquiring temperature parameters on a welding component in the welding process, and selecting any point in the weldment to analyze the temperature change condition.
FIG. 4 shows the temperature distribution during the 1 st, 3 rd, 5 th and 7 th welding in the multi-pass welding of the intersecting line of the tee pipe fitting, the central temperature of the heat source is about 1783 ℃ when the first welding is started, the central temperature of the heat source is about 1691 ℃ when the third welding is started, the central temperature of the heat source is 1792 ℃ when the fifth welding is started, the central temperature of the heat source is 1504 ℃ when the seventh welding is started, the melting point of Q345R steel is 1500 ℃, and the gray area in the figure is a welding molten pool area which is similar to the shape of an actual welding molten pool.
The change situation of the temperature of a certain point of the weldment along with the time is visually displayed, in the example, a mode of drawing a heat cycle curve of a certain point of a welding seam area is adopted, and as shown in fig. 5, the heat cycle curve of each point in each welding is shown.
In summary, compared with the prior art, the method of the embodiment has the following beneficial effects:
(1) the numerical simulation prediction method for the three-way intersection line welding temperature field can be used for realizing the simulation calculation of the three-way intersection line welding process temperature field.
(2) In the embodiment, a double-ellipsoid heat source model which best accords with an actual welding heat source is adopted, and a welding path is a complex space curve, so that a local coordinate system is established to accurately describe the shape of the double-ellipsoid heat source, the accurate simulation of an intersecting line welding temperature field can be realized, the temperature change condition of any point in a weldment is analyzed, and a large amount of process tests are not required.
(3) The embodiment fits the welding track curve by capturing the coordinates of the nodes of the welding bead grid, can be suitable for any irregular space welding curve, and avoids the problem that the welding track is difficult to accurately describe through the traditional geometric parameter equation.
This embodiment still provides a tee bend intersection line welding temperature field numerical simulation prediction system, includes:
the welding model building module is used for building a three-dimensional model and carrying out meshing division on the three-dimensional model, wherein the three-dimensional model comprises a three-way pipe fitting and a welding seam;
the heat source model building module is used for building a welded heat source model and obtaining a moving track curve of the heat source model according to the grid node coordinates of a welding bead in the three-dimensional model;
and the temperature field prediction module is used for carrying out simulation analysis by combining the heat source model, the track curve and the three-dimensional model to obtain temperature parameters on the welding component in the welding process.
The numerical simulation prediction system for the three-way intersection line welding temperature field can execute the numerical simulation prediction method for the three-way intersection line welding temperature field provided by the embodiment of the method, can execute any combination implementation steps of the embodiment of the method, and has corresponding functions and beneficial effects of the method.
The embodiment also provides a storage medium, which stores instructions or programs capable of executing the numerical simulation prediction method for the three-way intersection line welding temperature field provided by the embodiment of the method of the invention, and when the instructions or the programs are run, the steps can be executed in any combination of the embodiment of the method, so that the method has corresponding functions and beneficial effects.
More specific examples (a non-exhaustive list) of the computer-readable medium would include the following: an electrical connection (electronic device) having one or more wires, a portable computer diskette (magnetic device), a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber device, and a portable compact disc read-only memory (CDROM). Additionally, the computer-readable medium could even be paper or another suitable medium upon which the program is printed, as the program can be electronically captured, via for instance optical scanning of the paper or other medium, then compiled, interpreted or otherwise processed in a suitable manner if necessary, and then stored in a computer memory.
It should be understood that portions of the present invention may be implemented in hardware, software, firmware, or a combination thereof. In the above embodiments, the various steps or methods may be implemented in software or firmware stored in memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, any one or combination of the following techniques, which are known in the art, may be used: a discrete logic circuit having a logic gate circuit for implementing a logic function on a data signal, an application specific integrated circuit having an appropriate combinational logic gate circuit, a Programmable Gate Array (PGA), a Field Programmable Gate Array (FPGA), or the like.
In the foregoing description of the specification, reference to the description of "one embodiment/example," "another embodiment/example," or "certain embodiments/examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
While embodiments of the present invention have been shown and described, it will be understood by those of ordinary skill in the art that: various changes, modifications, substitutions and alterations can be made to the embodiments without departing from the principles and spirit of the invention, the scope of which is defined by the claims and their equivalents.
While the preferred embodiments of the present invention have been illustrated and described, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention as defined by the appended claims.
Claims (10)
1. A numerical simulation prediction method for a three-way intersection line welding temperature field is characterized by comprising the following steps:
establishing a three-dimensional model, and performing grid division on the three-dimensional model, wherein the three-dimensional model comprises a three-way pipe fitting and a welding seam;
establishing a welded heat source model, and acquiring a moving track curve of the heat source model according to the grid node coordinates of a welding bead in the three-dimensional model;
and (4) performing simulation analysis by combining the heat source model, the track curve and the three-dimensional model to obtain temperature parameters on the welding component in the welding process.
2. The method for predicting the numerical simulation of the three-way intersection line welding temperature field according to claim 1, wherein the establishing of the three-dimensional model and the meshing of the three-dimensional model comprise:
establishing three-dimensional solid models of the three-way pipe fitting and the welding seam, and combining the three-dimensional solid models of the three-way pipe fitting and the welding seam into an integral geometric model;
meshing the geometric model according to the sequence of diffusion from the welding seam to the periphery;
predefining physical performance parameters of the pipe fitting and the welding material along with temperature change, and respectively endowing material properties to corresponding structural parts.
3. The method for predicting the numerical simulation of the three-way intersection line welding temperature field according to claim 2, wherein the meshing the geometric model according to the sequence of diffusion from the welding seam to the periphery comprises the following steps:
segmenting an integral model of the geometric model into a plurality of entities;
2D meshing the surface of the entity, and sweeping the generated volumetric mesh.
4. The three-way intersection line welding temperature field numerical simulation prediction method of claim 1, wherein the establishing of the welding heat source model and the obtaining of the moving trajectory curve of the heat source model according to the grid node coordinates of the welding bead in the three-dimensional model comprise:
acquiring grid node coordinates of the heat source model moving in each welding pass, and fitting the grid node coordinates to obtain a path of the heat source model moving along with time;
based on a double-ellipsoid heat source model, aiming at the characteristic that the welding of the three-way intersection line has a spatial path, establishing a local coordinate system for the heat source model;
and establishing a heat source model for three-way intersecting line welding according to a local coordinate system.
5. The method for predicting the numerical simulation of the three-way intersection line welding temperature field according to claim 4, wherein the establishing a local coordinate system for the heat source model comprises the following steps:
setting the rectangular coordinate system of the three-dimensional model as the global coordinate system x, y and z and the local coordinate system x of the heat source modeli、yi、zi(ii) a Local coordinate system xiCoincidence of tangent vectors of axis and weld curve, ziThe axis is the normal vector of the discretization grid surface unit of the welding bead, yiAxis xiAxial coordinate vector cross-product ziAn axial coordinate vector;
the transformation relation between the global coordinate system and the local coordinate system is as follows:
wherein n isx1(ti) Is the component of the coordinate axis vector of the heat source model in the local coordinate system in the global coordinate system, Fx(ti)、Fy(ti)、Fz(ti) Is tiAnd the central position of the heat source at the moment is a coordinate value under the global coordinate system.
6. The numerical simulation prediction method for the three-way intersection line welding temperature field according to claim 5, wherein the establishing of the integral heat source model for the three-way intersection line welding according to the local coordinate system comprises the following steps:
writing heat source subprogram by using Fortran programming language, t during weldingiThe central position of the heat source is (x)i,yi,zi) Establishing a heat source model for three-way intersecting line welding;
the heat source model is divided into the following two parts:
ellipsoidal internal heat source distribution function along the first half of the welding direction:
ellipsoidal internal heat source distribution function in the second half of the welding direction:
wherein q is1、q2As a function of heat source distribution; a isf、arB and c are shape parameters of a double-ellipsoid heat source; eta is welding thermal efficiency; u is welding voltage; i is welding current; f. offThe heat distribution coefficient of the front half part of the heat source model; f. ofrThe heat distribution coefficient of the latter half of the heat source model.
7. The numerical simulation prediction method for the three-way intersection line welding temperature field according to claim 1, characterized in that in the simulation analysis process, a differential equation based on heat transfer control is adopted for heat transfer calculation analysis:
wherein ρ is density; c is the specific heat capacity of the material; t is a temperature distribution function; t is the heat transfer time;
before the simulation analysis, the method further comprises the following steps:
the initial temperature of the welded component and the boundary conditions of the temperature field are set.
8. A tee intersection line welding temperature field numerical simulation prediction system is characterized by comprising:
the welding model building module is used for building a three-dimensional model and carrying out meshing division on the three-dimensional model, wherein the three-dimensional model comprises a three-way pipe fitting and a welding seam;
the heat source model building module is used for building a welded heat source model and obtaining a moving track curve of the heat source model according to the grid node coordinates of a welding bead in the three-dimensional model;
and the temperature field prediction module is used for carrying out simulation analysis by combining the heat source model, the track curve and the three-dimensional model to obtain temperature parameters on the welding component in the welding process.
9. A tee intersection line welding temperature field numerical simulation prediction system is characterized by comprising:
at least one processor;
at least one memory for storing at least one program;
when executed by the at least one processor, cause the at least one processor to implement the method of any one of claims 1-7.
10. A storage medium having stored therein a program executable by a processor, wherein the program executable by the processor is adapted to perform the method of any one of claims 1-7 when executed by the processor.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110761134.5A CN113673124A (en) | 2021-07-06 | 2021-07-06 | Numerical simulation prediction method, system and medium for three-way intersection line welding temperature field |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110761134.5A CN113673124A (en) | 2021-07-06 | 2021-07-06 | Numerical simulation prediction method, system and medium for three-way intersection line welding temperature field |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113673124A true CN113673124A (en) | 2021-11-19 |
Family
ID=78538901
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110761134.5A Pending CN113673124A (en) | 2021-07-06 | 2021-07-06 | Numerical simulation prediction method, system and medium for three-way intersection line welding temperature field |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113673124A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115358115A (en) * | 2022-08-04 | 2022-11-18 | 北京建工集团有限责任公司 | Temperature field analysis method based on actually measured welding temperature field and combined with finite element |
CN116822098A (en) * | 2023-08-18 | 2023-09-29 | 广州誉鑫精密部件有限公司 | Intelligent monitoring control method and system for welding |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101628174B1 (en) * | 2014-12-09 | 2016-06-08 | 국방과학연구소 | Method for generating synthetic infrared image with mutual heat exchange effect between a moving object embedded with heat and a surroundings |
US20170337307A1 (en) * | 2016-05-20 | 2017-11-23 | Dassault Systemes Simulia Corp. | Scalable finite element simulation of additive manufacturing |
CN107967376A (en) * | 2017-10-30 | 2018-04-27 | 钢铁研究总院 | The T-shaped bilateral interval numerical simulation of welding processes computational methods of steel aluminium transit joint |
CN109325292A (en) * | 2018-09-25 | 2019-02-12 | 辽宁石油化工大学 | A kind of complicated welding track heat source load local coordinate system setting method based on ANSYS |
CN111428416A (en) * | 2020-04-20 | 2020-07-17 | 苏州科技大学 | Finite element modeling method for simulating high-energy beam welding |
CN111597739A (en) * | 2020-04-09 | 2020-08-28 | 西安理工大学 | Numerical simulation method for welding tube plate structure by using dissimilar materials |
CN111950182A (en) * | 2020-08-20 | 2020-11-17 | 南京中车浦镇城轨车辆有限责任公司 | Model simplification method for improving finite element calculation efficiency of large complex component |
CN112182921A (en) * | 2020-08-31 | 2021-01-05 | 北京化工大学 | Prediction method for selective laser melting thermal-mechanical coupling behavior of high-performance alloy steel |
-
2021
- 2021-07-06 CN CN202110761134.5A patent/CN113673124A/en active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101628174B1 (en) * | 2014-12-09 | 2016-06-08 | 국방과학연구소 | Method for generating synthetic infrared image with mutual heat exchange effect between a moving object embedded with heat and a surroundings |
US20170337307A1 (en) * | 2016-05-20 | 2017-11-23 | Dassault Systemes Simulia Corp. | Scalable finite element simulation of additive manufacturing |
CN107967376A (en) * | 2017-10-30 | 2018-04-27 | 钢铁研究总院 | The T-shaped bilateral interval numerical simulation of welding processes computational methods of steel aluminium transit joint |
CN109325292A (en) * | 2018-09-25 | 2019-02-12 | 辽宁石油化工大学 | A kind of complicated welding track heat source load local coordinate system setting method based on ANSYS |
CN111597739A (en) * | 2020-04-09 | 2020-08-28 | 西安理工大学 | Numerical simulation method for welding tube plate structure by using dissimilar materials |
CN111428416A (en) * | 2020-04-20 | 2020-07-17 | 苏州科技大学 | Finite element modeling method for simulating high-energy beam welding |
CN111950182A (en) * | 2020-08-20 | 2020-11-17 | 南京中车浦镇城轨车辆有限责任公司 | Model simplification method for improving finite element calculation efficiency of large complex component |
CN112182921A (en) * | 2020-08-31 | 2021-01-05 | 北京化工大学 | Prediction method for selective laser melting thermal-mechanical coupling behavior of high-performance alloy steel |
Non-Patent Citations (3)
Title |
---|
宋芳芳;金亮;恩云飞;李斌;: "异种金属激光焊接的温度场分析方法及应用", 焊接技术, no. 12, 28 December 2016 (2016-12-28), pages 9 - 13 * |
马小明;卢冠鹏;: "基于Simufact的CRDM管座焊接残余应力数值模拟", 焊接技术, no. 05, 28 May 2017 (2017-05-28), pages 51 - 54 * |
马小明;欧清扬;: "盲孔法测曲面残余应力时释放系数的数值模拟", 华南理工大学学报(自然科学版), no. 12, 15 December 2019 (2019-12-15), pages 25 - 31 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115358115A (en) * | 2022-08-04 | 2022-11-18 | 北京建工集团有限责任公司 | Temperature field analysis method based on actually measured welding temperature field and combined with finite element |
CN116822098A (en) * | 2023-08-18 | 2023-09-29 | 广州誉鑫精密部件有限公司 | Intelligent monitoring control method and system for welding |
CN116822098B (en) * | 2023-08-18 | 2024-01-16 | 广州誉鑫精密部件有限公司 | Intelligent monitoring control method and system for welding |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113673124A (en) | Numerical simulation prediction method, system and medium for three-way intersection line welding temperature field | |
CN109226937B (en) | Industrial robot spatial intersecting curve welding off-line programming method | |
CN109926767B (en) | Integrated computing method for welding residual stress and deformation of super-large structure | |
Hu et al. | Simulation and analysis of assembly processes considering compliant, non-ideal parts and tooling variations | |
CN107563072B (en) | Finite element modeling method for weld nut based on automobile fatigue simulation | |
CN103246774B (en) | The method of numerical simulation P92 steel-pipe welding heat-affected zone width | |
CN102637235A (en) | Determination method for heat source model parameters in multiplewire submerged-arc welding by numerical simulation | |
CN111597739A (en) | Numerical simulation method for welding tube plate structure by using dissimilar materials | |
CN112926213B (en) | Method, system, medium, equipment and terminal for measuring thermal damage boundary element | |
CN105069183A (en) | Weld simulation method and device | |
Feng et al. | Energy metrics for product assembly equipment and processes | |
CN102300668B (en) | Automatic estimation of weldgun size using section geometry | |
CN112597603B (en) | PDMS pipeline automatic modeling method based on key points and computer terminal | |
CN105528480A (en) | Welding spot analysis method and apparatus | |
Shahi et al. | Quality-driven optimization of assembly line configuration for multi-station assembly systems with compliant non-ideal sheet metal parts | |
CN115293054A (en) | Neural network-based dented pipeline failure evaluation method | |
CN113158380B (en) | Numerical simulation calculation method for multi-layer multi-pass welding process of header tube seat | |
Lorin et al. | Simulation of non-nominal welds by resolving the melted zone and its implication to variation simulation | |
JP2009128085A (en) | Residual stress analysis method by finite element method | |
CN113624840A (en) | Method, system, equipment and medium for detecting welding seam of BOSS head of nuclear power station | |
CN112395782A (en) | CAE-based method and device for calculating coordinates of welding points in batch | |
CN204790978U (en) | Welded analogue means | |
Gkatzogiannis et al. | Simulation of Welding Residual Stresses in Large Structures–First Results of Thermal Analyses | |
Khandagale et al. | Numerical Approach to Welding Process and its Integration in Assessment of Fatigue life of Component | |
Schoenborn | Fatigue analysis of a welded assembly using ANSYS Workbench environment |
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 |