CN111695256B - Modeling method of laser arc composite heat source based on energy distribution coefficient - Google Patents

Modeling method of laser arc composite heat source based on energy distribution coefficient Download PDF

Info

Publication number
CN111695256B
CN111695256B CN202010520886.8A CN202010520886A CN111695256B CN 111695256 B CN111695256 B CN 111695256B CN 202010520886 A CN202010520886 A CN 202010520886A CN 111695256 B CN111695256 B CN 111695256B
Authority
CN
China
Prior art keywords
heat source
laser
model
double
ellipsoid
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010520886.8A
Other languages
Chinese (zh)
Other versions
CN111695256A (en
Inventor
严春妍
张�浩
朱子江
易思
黄晏程
侯佳兵
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changzhou Campus of Hohai University
Original Assignee
Changzhou Campus of Hohai University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Changzhou Campus of Hohai University filed Critical Changzhou Campus of Hohai University
Priority to CN202010520886.8A priority Critical patent/CN111695256B/en
Publication of CN111695256A publication Critical patent/CN111695256A/en
Application granted granted Critical
Publication of CN111695256B publication Critical patent/CN111695256B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling

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)
  • Laser Beam Processing (AREA)

Abstract

The invention discloses a modeling method of a laser-arc composite heat source based on an energy distribution coefficient in the technical field of laser welding, and aims to solve the technical problems that the heat source check period is long and the coincidence degree of a simulated molten pool and an actual molten pool shape is not high when welding simulation is performed in the prior art. The method comprises the following steps: acquiring actual molten pool parameters; the method comprises the steps of constructing a laser arc composite heat source model based on an energy distribution coefficient model between any two pre-constructed component heat sources, obtaining a simulation molten pool parameter based on the laser arc composite heat source model, comparing the simulation molten pool parameter with an actual molten pool parameter to obtain the goodness of fit of the simulation molten pool relative to the actual molten pool, judging whether the goodness of fit reaches a preset threshold value, adjusting the heat source parameter, and taking the laser arc composite heat source model when the goodness of fit reaches the preset threshold value as a final model of the laser arc composite heat source.

Description

Modeling method of laser arc composite heat source based on energy distribution coefficient
Technical Field
The invention relates to a modeling method of a laser-arc composite heat source based on an energy distribution coefficient, and belongs to the technical field of laser welding.
Background
With the rapid development of the calculation simulation technology, the welding simulation technology gradually draws attention, and more engineers study the welding process and phenomenon through the simulation technology, so as to optimize the welding process and effectively control the welding deformation and the welding defects. At present, the simulation technology related to welding is not mature enough, and the technical problems that the heat source check period is long, the coincidence degree of the shape of the simulated molten pool and the actual molten pool is not high and the like exist.
Disclosure of Invention
Aiming at the defects of the prior art, the invention aims to provide a modeling method of a laser-arc composite heat source based on an energy distribution coefficient, so as to solve the technical problems that the heat source check period is long and the fit degree between a simulated molten pool and an actual molten pool is not high when welding simulation is performed in the prior art.
In order to solve the technical problems, the technical scheme adopted by the invention is as follows:
a modeling method of a laser arc composite heat source based on an energy distribution coefficient comprises the following steps:
acquiring actual molten pool parameters;
constructing a laser arc composite heat source model based on an energy distribution coefficient model between any two pre-constructed component heat sources, wherein the component heat source is one of the heat sources forming the laser arc composite heat source, and variables of the energy distribution coefficient model comprise heat source parameters;
acquiring a simulated molten pool parameter based on the laser-arc composite heat source model, comparing the simulated molten pool parameter with the actual molten pool parameter to acquire the goodness of fit of the simulated molten pool relative to the actual molten pool, and judging whether the goodness of fit reaches a preset threshold value or not;
responding to the goodness of fit not reaching a preset threshold value, adjusting heat source parameters, then repeatedly simulating a molten pool and an obtaining process of the goodness of fit relative to an actual molten pool, and judging whether the goodness of fit reaches the preset threshold value or not;
and in response to the coincidence degree reaching a preset threshold value, taking the laser arc composite heat source model under the corresponding heat source parameter as a final model of the laser arc composite heat source.
Furthermore, the heat source comprises a double-ellipsoid heat source and two cylindrical heat sources, wherein the double-ellipsoid heat source is used for representing electric arcs and the two cylindrical heat sources are used for representing laser, the double-ellipsoid heat source acts on the upper half part of the action area of the laser and electric arc composite heat source, and the cylindrical heat source acts on the lower half part of the action area of the laser and electric arc composite heat source.
Further, the heat source parameter includes at least any one of a shape parameter of the double-ellipsoid heat source, a radius of the cylindrical heat source, an effective operation depth of the double-ellipsoid heat source, and an effective operation depth of the cylindrical heat source.
Further, the expression of the double-ellipsoid heat source model is as follows:
Figure BDA0002532008100000021
in the formula, q 1 (x, y, z) is a heat flow density function of the double-ellipsoid heat source model, x is a coordinate value of the double-ellipsoid heat source model in the x direction of a space coordinate system, y is a coordinate value of the double-ellipsoid heat source model in the y direction of the space coordinate system, z is a coordinate value of the double-ellipsoid heat source model in the z direction of the space coordinate system, A 1 Is the energy coefficient of a double ellipsoid heat source, f 1 (x, y, z) is a shape function of the double-ellipsoid heat source, a, b and c are shape parameters of the double-ellipsoid heat source, h 1 Effective depth of action, Q, of a double ellipsoidal heat source 1 Effective power of double ellipsoid heat source, Q 1 =η 1 P 1 ,η 1 Is the power effective coefficient, P, of a double ellipsoid heat source 1 The actual power of the double-ellipsoid heat source is shown, beta is the included angle between the arc main shaft and the x direction, gamma is the included angle between the arc main shaft and the y direction, and theta is the included angle between the arc main shaft and the z direction.
Further, the expression of the cylinder heat source model is as follows:
Figure BDA0002532008100000022
in the formula, q i+1 (x, y, z) is a heat flow density function of the ith cylinder heat source model, x is a coordinate value of the cylinder heat source model in the x direction of a space coordinate system, y is a coordinate value of the cylinder heat source model in the y direction of the space coordinate system, z is a coordinate value of the cylinder heat source model in the z direction of the space coordinate system, A i+1 Is the energy coefficient of the ith cylinder heat source; f. of i+1 (x, y, z) is a shape function of the ith cylinder heat source; h is i+1 Is the effective depth of action, r, of the ith cylinder heat source i Radius of the ith cylinder heat source, Q i+1 Effective power of the i-th cylinder heat source, Q i+1 =η i+1 P i+1 ,η i+1 Is the power efficiency coefficient, P, of the ith cylinder heat source i+1 R (z) is a heat flow distribution function of the two cylindrical heat sources.
Further, the laser arc composite heat source model has the following expression:
when z is more than or equal to 0 and less than or equal to h 1 When the temperature of the water is higher than the set temperature,
Figure BDA0002532008100000031
when h is generated 1 ≤z≤h 1 +h 2 When the temperature of the water is higher than the set temperature,
Figure BDA0002532008100000032
when h is generated 1 +h 2 ≤z≤h 1 +h 2 +h 3 When the temperature of the water is higher than the set temperature,
Figure BDA0002532008100000033
wherein q (x, y, z) is a heat flow density function of the laser arc composite heat source model, x is a coordinate value of the laser arc composite heat source model in the x direction of a space coordinate system, y is a coordinate value of the laser arc composite heat source model in the y direction of the space coordinate system, z is a coordinate value of the laser arc composite heat source model in the z direction of the space coordinate system, h 1 Is the effective action depth h of the double ellipsoid heat source 2 Is the energy coefficient of the 1 st cylinder heat source, h 3 Is the energy coefficient, Q, of the 2 nd cylinder heat source 1 Effective power of double ellipsoid heat source, Q 2 Effective power of the 1 st cylinder heat source, Q 3 Is the effective power of the 2 nd cylinder heat source, beta is the included angle between the main axis of the electric arc and the x direction, gamma is the included angle between the main axis of the electric arc and the y direction, theta is the included angle between the main axis of the electric arc and the z direction, a, b and c are the shape parameters of the double ellipsoid heat source, r 1 Radius of the 1 st cylinder heat source, r 2 Radius of the 2 nd cylinder heat source, R (z) is a heat flow distribution function of the two cylinder heat sources, f 1 Energy distribution coefficient f of double-ellipsoid heat source constructed according to interaction relation of electric arc and laser 2 Constructed according to the interaction of laser and electric arcEnergy distribution coefficient of first cylinder heat source, f 3 The energy distribution coefficient of the second cylinder heat source constructed according to the interaction relationship of the laser and the laser.
Further, the energy distribution coefficient model is expressed as follows:
Figure BDA0002532008100000041
Figure BDA0002532008100000042
Figure BDA0002532008100000043
further, the laser-arc hybrid heat source comprises a laser-MIG hybrid heat source.
Compared with the prior art, the invention has the following beneficial effects: according to an interaction mechanism of electric arcs and laser, a laser electric arc composite heat source model is established based on an energy distribution coefficient model, a molten pool parameter is simulated by the laser electric arc composite heat source model, the simulated molten pool parameter is compared with an actual molten pool parameter to obtain goodness of fit, the goodness of fit can be adjusted by adjusting a heat source parameter in the energy distribution coefficient model, and the laser electric arc composite heat source model with the goodness of fit reaching a preset threshold value is extracted as a final model. Because the laser arc composite heat source model is constructed based on the energy distribution coefficient model, when heat source checking is carried out, heat source parameters including the size of the heat source model and the like are used as initial values to be input into the energy distribution coefficient model, then the energy distribution coefficient model is led into the laser arc composite heat source model, the shape of a simulated molten pool and the shape of an actual workpiece molten pool obtained through numerical calculation are compared, basic matching can be achieved only by finely adjusting the energy input values and other related heat source parameters, repeated adjustment is not needed, the checking period of the composite heat source is greatly shortened, and the time consumption of the whole welding simulation is further shortened.
Drawings
FIG. 1 is a schematic flow diagram of an embodiment of the method of the present invention;
FIG. 2 is a schematic diagram showing the composition of a laser-arc hybrid heat source model in an embodiment of the method of the present invention;
FIG. 3 is a schematic view of a composite weld joint grid model in an embodiment of the method of the present invention;
FIG. 4 is a schematic cross-sectional view comparing the cross-sectional shape of an experimental weld with the simulated weld obtained by the method of the present invention;
FIG. 5 is a schematic cross-sectional view comparing the experimental weld seam in the second embodiment of the method of the present invention with the simulated weld seam obtained by the method of the present invention.
Detailed Description
The invention is further described below with reference to the accompanying drawings. The following examples are only for illustrating the technical solutions of the present invention more clearly, and the protection scope of the present invention is not limited thereby.
The specific embodiment of the invention provides a modeling method of a laser arc composite heat source based on an energy distribution coefficient, as shown in fig. 1, a flow diagram of an embodiment of the method of the invention is as follows: firstly, acquiring actual molten pool parameters; then, according to the action area of the laser-arc composite heat source in metal, representing the arc by using a double-ellipsoid heat source model, representing the laser by using two cylindrical heat source models, and obtaining an energy distribution coefficient model of the double-ellipsoid heat source and the two cylindrical heat sources based on the interaction mechanism of the arc and the laser, wherein variables of the energy distribution coefficient model are provided with a plurality of heat source parameters; then, constructing an initial model of the laser-arc composite heat source based on the energy distribution coefficient model, and obtaining a molten pool in a simulated manner based on the initial model of the laser-arc composite heat source so as to obtain parameters of the simulated molten pool; and then, comparing the obtained simulated molten pool parameters with the actual molten pool parameters to obtain the goodness of fit of the simulated molten pool relative to the actual molten pool. And when the goodness of fit of the simulated molten pool relative to the actual molten pool reaches a preset threshold, extracting a laser arc composite heat source model under the corresponding heat source parameter to serve as a determination model of the laser arc composite heat source. Based on the technical principle, when the initial model of the laser-arc composite heat source is adjusted, only the heat source parameters in the energy distribution coefficient model need to be adjusted, and a complicated adjusting program is not needed, so that the checking period of the composite heat source is greatly shortened; the energy distribution coefficient model obtained by deduction based on the method has great correlation with the shape of the molten pool, and the error of distributing heat source energy according to experience in the past is abandoned, so that the shape of the molten pool obtained by simulation is more consistent with the actual shape of the molten pool. The method specifically comprises the following steps:
(1) analyzing the characteristics of the action area of the laser-arc composite heat source in the metal
In the laser-arc composite heat source, a double-ellipsoid heat source and two cylindrical heat sources are combined, the double-ellipsoid heat source acts on the upper half part of the action area of the laser-arc composite heat source, and the two cylindrical heat sources act on the lower half part of the action area of the laser-arc composite heat source in the effective action depth of the laser-arc composite heat source. In this embodiment, the double-ellipsoid heat source and the cylindrical heat source may be set as the constituent heat sources.
(2) Establishing a double-ellipsoid heat source model
Figure BDA0002532008100000051
This formula can be deformed into:
Figure BDA0002532008100000061
in the formula, q 1 (x, y, z) is a heat flow density function of the double-ellipsoid heat source model, x is a coordinate value of the double-ellipsoid heat source model in the x direction of a space coordinate system, y is a coordinate value of the double-ellipsoid heat source model in the y direction of the space coordinate system, z is a coordinate value of the double-ellipsoid heat source model in the z direction of the space coordinate system, A 1 Is the energy coefficient of a double ellipsoid heat source, f 1 (x, y, z) is a double ellipsoid heat sourceThe shape function a, b and c are the shape parameters of the double-ellipsoid heat source h 1 Effective depth of action, Q, of a double ellipsoidal heat source 1 Effective power of double ellipsoid heat source, Q 1 =η 1 P 1 ,η 1 Is the power effective coefficient, P, of a double ellipsoid heat source 1 The actual power of the double-ellipsoid heat source is shown, beta is the included angle between the arc main shaft and the x direction, gamma is the included angle between the arc main shaft and the y direction, and theta is the included angle between the arc main shaft and the z direction.
(3) Establishing two cylinder heat source models
Figure BDA0002532008100000062
This formula can be deformed into:
Figure BDA0002532008100000063
in the formula, q i+1 (x, y, z) is a heat flow density function of the ith cylinder heat source model, x is a coordinate value of the cylinder heat source model in the x direction of a space coordinate system, y is a coordinate value of the cylinder heat source model in the y direction of the space coordinate system, z is a coordinate value of the cylinder heat source model in the z direction of the space coordinate system, A i+1 Is the energy coefficient of the ith cylinder heat source; f. of i+1 (x, y, z) is a shape function of the ith cylinder heat source; h is i+1 Is the effective depth of action, r, of the ith cylinder heat source i Radius of the ith cylinder heat source, Q i+1 Effective power of the i-th cylinder heat source, Q i+1 =η i+1 P i+1 ,η i+1 Is the power efficiency coefficient, P, of the ith cylinder heat source i+1 R (z) is a heat flow distribution function of the two cylindrical heat sources.
(4) Compounding two heat source models to establish a compound heat source model consisting of three heat sources
The laser arc composite heat source model is divided into an upper part and a lower part, the upper part is a double-ellipsoid heat source, the lower part is two cylindrical heat sources, and a laser arc composite heat source model q (x, y, z) is formed by combination, wherein the expression is as follows:
when z is more than or equal to 0 and less than or equal to h 1 When the temperature of the water is higher than the set temperature,
Figure BDA0002532008100000071
when h is generated 1 ≤z≤h 1 +h 2 When the temperature of the water is higher than the set temperature,
Figure BDA0002532008100000072
when h is generated 1 +h 2 ≤z≤h 1 +h 2 +h 3 When the temperature of the water is higher than the set temperature,
Figure BDA0002532008100000073
wherein q (x, y, z) is a heat flow density function of the laser arc composite heat source model, x is a coordinate value of the laser arc composite heat source model in the x direction of a space coordinate system, y is a coordinate value of the laser arc composite heat source model in the y direction of the space coordinate system, z is a coordinate value of the laser arc composite heat source model in the z direction of the space coordinate system, h 1 Is the effective action depth h of the double ellipsoid heat source 2 Is the energy coefficient of the 1 st cylinder heat source, h 3 Is the energy coefficient, Q, of the 2 nd cylinder heat source 1 Effective power of double ellipsoid heat source, Q 2 Effective power of the 1 st cylinder heat source, Q 3 Is the effective power of the 2 nd cylinder heat source, beta is the included angle between the main axis of the electric arc and the x direction, gamma is the included angle between the main axis of the electric arc and the y direction, theta is the included angle between the main axis of the electric arc and the z direction, a, b and c are the shape parameters of the double ellipsoid heat source, r 1 Radius of the 1 st cylinder heat source, r 2 Radius of the 2 nd cylinder heat source, R (z) is a heat flow distribution function of the two cylinder heat sources, f 1 Energy distribution coefficient f of double-ellipsoid heat source constructed according to interaction relation of electric arc and laser 2 Is root ofEnergy distribution coefficient, f, of the first cylindrical heat source constructed from the interaction of laser and arc 3 The energy distribution coefficient of the second cylinder heat source constructed according to the interaction relationship of the laser and the laser.
(5) Establishing an energy distribution proportion model
f 1 、f 2 And f 3 The relationship of the energy distribution coefficient as a composite heat source model is as follows: f. of 1 +f 2 +f 3 1. And the energy distribution coefficient has great correlation with the size of the heat source model. Thus according to the above A 1 A and A 3 The formula of (c) can derive the energy distribution coefficient as:
Figure BDA0002532008100000081
further, it is possible to obtain:
Figure BDA0002532008100000082
Figure BDA0002532008100000083
Figure BDA0002532008100000084
the heat source parameters comprise shape parameters of the double-ellipsoid heat source, radius of the cylindrical heat source, effective acting depth of the double-ellipsoid heat source, and effective acting depth of the cylindrical heat source, i.e. a, b, c, r in the above formula 1 、r 2 、h 1 、h 2 、h 3 And the like.
(6) Simulation of temperature field
Establishing a finite element model, taking the welding voltage, the welding current, the laser power, the welding speed and the welding inclination angle of two groups of medium plates compounding multiple channels as known parameters, adjusting the action area of a laser-arc compound heat source through the shape of a fusion line to obtain the corresponding heat source model size, inputting the heat source model size as an initial value, and then importing an energy distribution coefficient model into the laser-arc compound heat source model to obtain the corresponding simulated molten pool shape. And (3) determining whether the constructed composite heat source energy distribution coefficient model is accurate or not by using the coincidence degree of the simulated molten pool shape and the actual molten pool shape to obtain an optimal numerical simulation heat source model. In the embodiment of the method, the laser-arc composite heat source is a laser-MIG composite heat source.
This process is described in detail below with two examples.
The first embodiment is as follows:
an alloy steel plate of 150mm x 80mm x 18.4mm is subjected to cutting processing, a V-shaped groove of 30 degrees is formed, and three composite heat source multi-pass welding are carried out. The X100 high strength pipeline steel may have chips and oil stains on its surface after cutting and surface treatment. Therefore, the test specimens were sanded with sandpaper before welding, then cleaned with acetone, and finally wiped with alcohol. The welding was carried out with the process parameters as shown in table 1.
Table 1: example one adopted welding Process parameters
Figure BDA0002532008100000091
And obtaining the macroscopic morphology of the laser-arc hybrid welding seam through a hybrid welding experiment. Specifically, as shown in fig. 3, a schematic diagram of a mesh model of a composite welding joint in the embodiment of the method of the present invention is shown. Then, various sizes of the laser-arc composite welding heat source model are obtained according to the appearance of the welding seam, and then a temperature field is obtained according to the composite heat source model formula provided by the invention, namely the laser-arc composite heat source model based on the energy distribution coefficient, so that the simulated welding pool shape is obtained.
Specifically, as shown in FIG. 4, it is a schematic diagram comparing the cross-sectional shapes of the experimental weld and the simulated weld obtained by the method of the present invention in the first embodiment of the method of the present invention, the left side is the shape of the actual molten pool obtained by the experiment, and the right side is the laser-arc recombination obtained by the method of the present inventionThe shape of a simulated molten pool obtained by a heat source model is shown, wherein the part higher than 1460 degrees is a welding seam, A 1 、A 2 、A 3 For penetration, B 1 、B 2 、B 3 The melt width is measured. The shape of the molten pool in the temperature field obtained by the laser-arc composite heat source model obtained by the method is compared with the shape of the molten pool in the experimental result, so that the similarity between the shape of the molten pool in the temperature field and the shape of the molten pool in the experimental result is higher.
The second embodiment:
an alloy steel plate with the thickness of 150mm multiplied by 80mm multiplied by 18.4mm is subjected to cutting processing, a V-shaped groove with the angle of 30 degrees is formed, and three composite heat source multi-pass welding are carried out, specifically as shown in fig. 2, the composite heat source model is a schematic composition diagram of a laser arc composite heat source model in the embodiment of the method. The X100 high strength pipeline steel may have chips and oil stains on its surface after cutting and surface treatment. Therefore, the test specimens were sanded with sandpaper before welding, then cleaned with acetone, and finally wiped with alcohol. The welding was performed using the process parameters as shown in table 2.
Table 2: example two welding Process parameters
Figure BDA0002532008100000092
And obtaining the macroscopic morphology of the laser-arc hybrid welding seam through a hybrid welding experiment. Then, various sizes of the laser-arc composite welding heat source model are obtained according to the appearance of the welding seam, and then a temperature field is obtained according to the composite heat source model formula provided by the invention, namely the laser-arc composite heat source model based on the energy distribution coefficient, so that the simulated welding pool shape is obtained.
Specifically, as shown in fig. 5, a schematic diagram of a comparison between the cross-sectional shapes of the experimental weld and the simulated weld obtained by the method in the second embodiment of the method of the present invention is shown, where the left side is the shape of the actual molten pool obtained by the experiment, and the right side is the shape of the simulated molten pool obtained by the laser-arc composite heat source model obtained by the method of the present invention. In the figure, the part higher than 1460 degrees is a welding seam, C 1 、C 2 、C 3 For penetration, D 1 、D 2 、D 3 The melt width is measured. Comparing the left side shape and the right side shape, the similarity between the shape of the molten pool of the temperature field obtained by the laser arc composite heat source model obtained by the method and the shape of the molten pool of the experimental result is higher.
In summary, the method of the present invention has the following advantages: because the laser arc composite heat source model is constructed based on the energy distribution coefficient model, when heat source checking is carried out, heat source parameters including the size of the heat source model and the like are used as initial values to be input into the energy distribution coefficient model, then the energy distribution coefficient model is led into the laser arc composite heat source model, the shape of a simulated molten pool and the shape of an actual workpiece molten pool obtained through numerical calculation are compared, basic matching can be achieved only by finely adjusting the energy input values and other related heat source parameters, repeated adjustment is not needed, the checking period of the composite heat source is greatly shortened, and the time consumption of the whole welding simulation is further shortened.
The above description is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, several modifications and variations can be made without departing from the technical principle of the present invention, and these modifications and variations should also be regarded as the protection scope of the present invention.

Claims (6)

1. A modeling method of a laser arc composite heat source based on an energy distribution coefficient is characterized by comprising the following steps:
acquiring actual molten pool parameters;
constructing a laser arc composite heat source model based on an energy distribution coefficient model between any two pre-constructed component heat sources, wherein the component heat source is one of the heat sources forming the laser arc composite heat source, and variables of the energy distribution coefficient model comprise heat source parameters;
acquiring a simulated molten pool parameter based on the laser-arc composite heat source model, comparing the simulated molten pool parameter with the actual molten pool parameter to acquire the goodness of fit of the simulated molten pool relative to the actual molten pool, and judging whether the goodness of fit reaches a preset threshold value or not;
responding to the goodness of fit not reaching a preset threshold value, adjusting heat source parameters, then repeatedly simulating a molten pool and an obtaining process of the goodness of fit relative to an actual molten pool, and judging whether the goodness of fit reaches the preset threshold value or not;
in response to that the goodness of fit reaches a preset threshold, taking a laser arc composite heat source model under corresponding heat source parameters as a final model of the laser arc composite heat source;
the heat source comprises a double-ellipsoid heat source for representing electric arc and two cylindrical heat sources for representing laser, wherein the double-ellipsoid heat source acts on the upper half part of the action area of the laser-electric arc composite heat source, and the cylindrical heat source acts on the lower half part of the action area of the laser-electric arc composite heat source;
the laser arc composite heat source model has the following expression:
when z is more than or equal to 0 and less than or equal to h 1 When the temperature of the water is higher than the set temperature,
Figure FDA0003696962990000011
when h is generated 1 ≤z≤h 1 +h 2 When the temperature of the water is higher than the set temperature,
Figure FDA0003696962990000012
when h is generated 1 +h 2 ≤z≤h 1 +h 2 +h 3 When the temperature of the water is higher than the set temperature,
Figure FDA0003696962990000013
wherein q (x, y, z) is a heat flow density function of the laser arc composite heat source model, x is a coordinate value of the laser arc composite heat source model in the x direction of a space coordinate system, y is a coordinate value of the laser arc composite heat source model in the y direction of the space coordinate system, z is a coordinate value of the laser arc composite heat source model in the z direction of the space coordinate system, h 1 Is the effective action depth h of the double ellipsoid heat source 2 As a 1 st cylinder heat sourceEnergy coefficient of (h) 3 Is the energy coefficient, Q, of the 2 nd cylinder heat source 1 Effective power of double ellipsoid heat source, Q 2 Effective power of the 1 st cylinder heat source, Q 3 Is the effective power of the 2 nd cylinder heat source, beta is the included angle between the main axis of the electric arc and the x direction, gamma is the included angle between the main axis of the electric arc and the y direction, theta is the included angle between the main axis of the electric arc and the z direction, a, b and c are the shape parameters of the double ellipsoid heat source, r 1 Radius of the 1 st cylinder heat source, r 2 Radius of the 2 nd cylinder heat source, R (z) is a heat flow distribution function of the two cylinder heat sources, f 1 Energy distribution coefficient f of double-ellipsoid heat source constructed according to interaction relation of electric arc and laser 2 Energy distribution coefficient f of the first cylindrical heat source constructed from the interaction of laser and arc 3 The energy distribution coefficient of the second cylinder heat source constructed according to the interaction relationship of the laser and the laser.
2. The modeling method of a laser arc composite heat source based on energy distribution coefficient as claimed in claim 1, wherein the heat source parameter comprises at least any one of a shape parameter of a double ellipsoid heat source, a radius of a cylinder heat source, an effective depth of action of a double ellipsoid heat source, and an effective depth of action of a cylinder heat source.
3. The modeling method of the laser arc composite heat source based on the energy distribution coefficient as claimed in claim 1, wherein the expression of the double-ellipsoid heat source model is as follows:
Figure FDA0003696962990000021
in the formula, q 1 (x, y, z) is a heat flow density function of the double-ellipsoid heat source model, x is a coordinate value of the double-ellipsoid heat source model in the x direction of a space coordinate system, y is a coordinate value of the double-ellipsoid heat source model in the y direction of the space coordinate system, z is a coordinate value of the double-ellipsoid heat source model in the z direction of the space coordinate system, A 1 Is a double ellipseEnergy coefficient of spherical heat source, f 1 (x, y, z) is a shape function of the double-ellipsoid heat source, a, b and c are shape parameters of the double-ellipsoid heat source, h 1 Effective depth of action, Q, of a double ellipsoidal heat source 1 Effective power of double ellipsoid heat source, Q 1 =η 1 P 1 ,η 1 Is the power effective coefficient, P, of a double ellipsoid heat source 1 The actual power of the double-ellipsoid heat source is shown, beta is the included angle between the arc main shaft and the x direction, gamma is the included angle between the arc main shaft and the y direction, and theta is the included angle between the arc main shaft and the z direction.
4. The modeling method of the energy distribution coefficient-based laser arc composite heat source as claimed in claim 1, wherein the expression of the cylindrical heat source model is as follows:
Figure FDA0003696962990000031
in the formula, q i+1 (x, y, z) is a heat flow density function of the ith cylinder heat source model, x is a coordinate value of the cylinder heat source model in the x direction of a space coordinate system, y is a coordinate value of the cylinder heat source model in the y direction of the space coordinate system, z is a coordinate value of the cylinder heat source model in the z direction of the space coordinate system, A i+1 Is the energy coefficient of the ith cylinder heat source; f. of i+1 (x, y, z) is a shape function of the ith cylinder heat source; h is i+1 Is the effective depth of action, r, of the ith cylinder heat source i Radius of the ith cylinder heat source, Q i+1 Effective power of the ith cylinder heat source, Q i+1 =η i+1 P i+1 ,η i+1 Is the power efficiency coefficient, P, of the ith cylinder heat source i+1 R (z) is a heat flow distribution function of the two cylindrical heat sources.
5. The modeling method of the energy distribution coefficient-based laser arc composite heat source according to claim 1, wherein the energy distribution coefficient model is expressed as follows:
Figure FDA0003696962990000032
Figure FDA0003696962990000033
Figure FDA0003696962990000034
6. the modeling method for an energy distribution coefficient based laser arc composite heat source of claim 1 wherein the laser arc composite heat source comprises a laser-MIG composite heat source.
CN202010520886.8A 2020-06-10 2020-06-10 Modeling method of laser arc composite heat source based on energy distribution coefficient Active CN111695256B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010520886.8A CN111695256B (en) 2020-06-10 2020-06-10 Modeling method of laser arc composite heat source based on energy distribution coefficient

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010520886.8A CN111695256B (en) 2020-06-10 2020-06-10 Modeling method of laser arc composite heat source based on energy distribution coefficient

Publications (2)

Publication Number Publication Date
CN111695256A CN111695256A (en) 2020-09-22
CN111695256B true CN111695256B (en) 2022-08-16

Family

ID=72480020

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010520886.8A Active CN111695256B (en) 2020-06-10 2020-06-10 Modeling method of laser arc composite heat source based on energy distribution coefficient

Country Status (1)

Country Link
CN (1) CN111695256B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113139314B (en) * 2021-04-29 2022-09-27 四川大学 Heat source numerical simulation method for laser additive manufacturing process
CN113779793B (en) * 2021-09-10 2023-09-26 电子科技大学 Heat source modeling method for selective laser melting based on ray tracing

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106529051A (en) * 2016-11-15 2017-03-22 河海大学常州校区 Method for determining heat source model parameters of single wire submerged arc welding numerical simulation
CN106670666A (en) * 2016-12-14 2017-05-17 南京航空航天大学 Construction method of energy distribution coefficient model of laser-electric-arc combined machining based on precise energy distribution

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106529051A (en) * 2016-11-15 2017-03-22 河海大学常州校区 Method for determining heat source model parameters of single wire submerged arc welding numerical simulation
CN106670666A (en) * 2016-12-14 2017-05-17 南京航空航天大学 Construction method of energy distribution coefficient model of laser-electric-arc combined machining based on precise energy distribution

Also Published As

Publication number Publication date
CN111695256A (en) 2020-09-22

Similar Documents

Publication Publication Date Title
CN111695256B (en) Modeling method of laser arc composite heat source based on energy distribution coefficient
CN105117547B (en) The emulation mode of micro- milling nickel base superalloy prediction residue stress
CN112149330B (en) Welding residual stress prediction and welding process optimization method for wind power tower oil seal platform
CN105385839A (en) System and method for automatic control over laser shock peening
CN106077901A (en) Electric arc based on hot-work die failure site increases material manufacture method
CN112380752B (en) Method for improving welding process of metal sheet by predicting welding heat treatment value of metal sheet
CN109933007A (en) A kind of method for building up of Laser-MIG Composite Welding heat source model
CN105279291A (en) Built-up welding technological parameter optimization method for batch processing
CN106670666A (en) Construction method of energy distribution coefficient model of laser-electric-arc combined machining based on precise energy distribution
CN112528535A (en) Mortise broaching process simulation analysis method based on heat-force-flow multi-field coupling
CN104657565B (en) The hot-work die design method of plesiotype face water channel
CN113618245A (en) X100 pipeline steel laser-MIG composite welding parameter optimization method based on numerical simulation
WO2024087524A1 (en) Oscillating laser welding equivalent heat source modeling method and oscillating laser welding simulation method
CN109885946B (en) Method for determining energy distribution of composite heat source and welding simulation method
Cai et al. Optimization of shielding gas composition in narrow gap GMA welding based on response surface methodology
CN116228766A (en) Intelligent regulation and control method and system for plasma processing equipment
CN111390299B (en) Method for predicting machining deformation of floating support friction plate
Wang et al. Research on shape change of multimaterial electrode for EDM
Chen et al. Optimization of processing parameters for beam blank continuous casting using MOGA combined with FEM
CN112163361A (en) Method for realizing dynamic refinement of simulation local grid in thin-wall casing machining process
CN110899522A (en) Process method for improving mold closing efficiency of automobile mold
CN112214863A (en) Method for acquiring coarse grain region characteristic heat cycle curve of austenitic stainless steel welding heat affected zone
CN116468187B (en) Green welding process decision method for large-diameter thick plate
CN113033039B (en) Method for establishing welding heat source model of single-wire electro-gas welding
Long et al. Parameter optimization using multi-objective taguchi method and response surface methodology for the laser welding of dissimilar materials

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant