CN106363283B - A kind of tandem double wire hidden arc welding numerical simulation heat source model determination method for parameter - Google Patents
A kind of tandem double wire hidden arc welding numerical simulation heat source model determination method for parameter Download PDFInfo
- Publication number
- CN106363283B CN106363283B CN201611064245.6A CN201611064245A CN106363283B CN 106363283 B CN106363283 B CN 106363283B CN 201611064245 A CN201611064245 A CN 201611064245A CN 106363283 B CN106363283 B CN 106363283B
- Authority
- CN
- China
- Prior art keywords
- parameter
- heat source
- welding
- double
- source model
- 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.)
- Expired - Fee Related
Links
- 238000003466 welding Methods 0.000 title claims abstract description 121
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000004088 simulation Methods 0.000 title claims abstract description 30
- 230000004927 fusion Effects 0.000 claims abstract description 32
- 230000035515 penetration Effects 0.000 claims abstract description 31
- 230000004907 flux Effects 0.000 claims abstract description 23
- 238000012887 quadratic function Methods 0.000 claims abstract description 23
- 238000009826 distribution Methods 0.000 claims abstract description 18
- 230000007246 mechanism Effects 0.000 claims abstract description 18
- 230000003993 interaction Effects 0.000 claims abstract description 7
- 238000005315 distribution function Methods 0.000 claims description 12
- 238000010891 electric arc Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 7
- 238000010206 sensitivity analysis Methods 0.000 claims description 7
- 230000035945 sensitivity Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 abstract description 7
- 238000012360 testing method Methods 0.000 abstract description 7
- 238000011161 development Methods 0.000 abstract description 4
- 238000002474 experimental method Methods 0.000 description 10
- 230000008018 melting Effects 0.000 description 6
- 238000002844 melting Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 229910000831 Steel Inorganic materials 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000010959 steel Substances 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 239000003500 flue dust Substances 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23K—SOLDERING OR UNSOLDERING; WELDING; CLADDING OR PLATING BY SOLDERING OR WELDING; CUTTING BY APPLYING HEAT LOCALLY, e.g. FLAME CUTTING; WORKING BY LASER BEAM
- B23K9/00—Arc welding or cutting
- B23K9/18—Submerged-arc welding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Plasma & Fusion (AREA)
- Mechanical Engineering (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Arc Welding In General (AREA)
Abstract
The invention discloses a kind of tandem double wire hidden arc welding numerical simulation heat source model determination method for parameter, include the following steps：Establish the double stripping mechanism that heat flow density decays along depth direction quadratic function；According to the interaction of forward and backward two heat sources, tandem double wire hidden arc welding heat flux distribution is divided into three regions；Finite element inversion is carried out to heat source model, obtains corresponding heat source model parameter；Heat source model parameter is adjusted with certain amplitude, finite element numerical simulation is carried out respectively to each parameter combination, obtains multigroup molten wide fusion penetration sample values；To determine the sensibility equation of double stripping mechanism parameter and speed of welding to molten wide, fusion penetration；Sensibility equation is fitted by multigroup sample values.The present invention can directly predict the heat source model parameter corresponding to different welding conditions using sensibility equation, the tentative calculation workload of welding analog be substantially reduced, and the precision of simulation is substantially improved, and reduce engineer testing, save development cost.
Description
Technical field
The invention belongs to welding value heat source model technical fields, and in particular to a kind of tandem double wire hidden arc welding numerical value
Simulate heat source model determination method for parameter.
Background technology
Submerged-arc welding is that a kind of electric arc burns the method welded under welding flux layer (containing submerged arc overlay welding and electroslag surfacing etc.).
Due to welding quality stable, welding productivity is high, seldom without arc light and flue dust the advantages that, become pressure vessel, pipeline section system
It makes, the main welding method in the important structure steel fabrication such as box beam column.In recent years, although successively occur many kinds efficiently,
Good new welding method, but the application field of submerged-arc welding is not affected still.Submerged-arc welding can according to welding wire number difference
Be divided into monofilament submerged-arc welding and multiplewire submerged arc welding, opposite monofilament submerged-arc welding, double wire hidden arc welding have can be improved speed of welding 30~
40%, welding quality is good, welding wire deposition rate and the characteristics of current utilization rate higher.
Submerged-arc welding numerical simulation primarily solves the problems, such as the problem of being submerged-arc welding heat source model parameter selection.Mariages buries at present
Arc-welding heat source model parameter determination method is mainly tentative calculation, due to the experience and time restriction of researcher, it is difficult to ensure that heat source
The precision of model, while development cost is increased again.In addition, due to not simple between two electric arcs of tandem double wire hidden arc welding
Superposition, two electric arcs, which exist, to influence each other.And in view of two interactional numerical simulation studies are seldom, therefore cause existing
There is the submerged-arc welding heat source model parameters precision simulated using simple superposition in technology not high.
Invention content
It is an object of the invention to overcome deficiency in the prior art, a kind of tandem double wire hidden arc welding numerical simulation is provided
Heat source model determination method for parameter solves model parameter in the prior art and determines the not high technical problem of of high cost, precision.
In order to solve the above technical problems, the present invention provides a kind of tandem double wire hidden arc welding numerical simulation heat source model parameters
Determination method, characterized in that include the following steps：
Step 1：The double stripping mechanism that heat flow density decays along depth direction quadratic function is established, determines mariages mould
The heat flux distribution function of type；
Attenuation function is quadratic function, and heat flux distribution function is：
Heat flux distribution function is in first half ellipsoid：
Latter half ellipsoid heat flux distribution function is：
F in formulafAnd frRespectively total input power is in the distribution index of molten bath front and rear part, and ff+fr=2, QiIt is i-th
The electric arc effective power of root wire, v are speed of welding；αiFor the welding inclination angle of the i-th root wire, afi、ariIt is corresponding for i-th silk
The length of the forward and backward hemisphere major semiaxis of double stripping mechanism, biFor the corresponding double semiminor axis of ellipsoid length of i-th silk, ciIt is i-th
The corresponding double ellipsoid depth of root silk, in formulaPart is that quadratic function is decayed, and t is welding process
The time of progress, i=1,2；
Tandem double wire hidden arc welding heat flux distribution is divided into three by step 2 according to the interaction of forward and backward two heat sources
A region determines the heat flow density in each region；
Tandem double wire hidden arc welding heat flux distribution can be divided into three areas by two double forward and backward hemisphere lines of demarcation of ellipsoid
Domain：1 heat flow density of region is double ellipsoid first halfs superposition of the double ellipsoid first halfs and rear silk of preceding silk；2 hot-fluid of region is close
Degree is double ellipsoid first halfs superposition of the double ellipsoid latter halfs and rear silk of preceding silk；3 heat flow density of region is the double ellipse of preceding silk
Double ellipsoid latter halfs of ball latter half and rear silk are superimposed；
Step 3 establishes finite element model, take the matched weldingvoltage of one group of submerged-arc welding, welding current, speed of welding and
Inclination angle is welded as simulation basic parameter, applies the heat source load of each region heat flow density description determined above, to the above heat
Source model parameter carries out finite element inversion, obtains corresponding heat source model parameter；
Step 4 is adjusted heat source model parameter obtained as above and speed of welding with the amplitude set, is based on each ginseng
Array closes to obtain corresponding molten wide, fusion penetration sample values；Sensitivity analysis is carried out to each group sample values, it is known that parameter afiWith
ariInfluence to molten wide, fusion penetration is smaller, with double stripping mechanism parameter：b1、b2、c1、c2It is input parameter with speed of welding v,
Using molten wide w, fusion penetration p as output parameter, corresponding sample is substituted into regression equation and obtains sensitivity analysis result：
Wherein, x2w、x3w、x4w、x5w、x6wFor b1、c1、b2、c2, v is to the sensitivity coefficient of molten wide w；x2p、x3p、x4p、x5p、
x6pFor b1、c1、b2、c2, v is to the sensitivity coefficient of molten wide p；x1w、x1pFor with Q1、Q2Related quadratic function f (Q1,Q2)、g(Q1,
Q2)；
Step 5 is fitted and simplifies to gained sensibility equation, and detailed process is：
Take x1w、x1pFor Q1、Q2Related quadratic function f (Q1,Q2)、g(Q1,Q2), the expression formula of quadratic function is as follows：
f(Q1,Q2)=α0+α1Q1+α2Q2+α3Q1Q2+α4Q1 2+α5Q2 2
g(Q1,Q2)=β0+β1Q1+β2Q2+β3Q1Q2+β4Q1 2+β5Q2 2
Wherein, αj, βjFor pre- fitting constant, j=0,1,2,3,4,5；
Multigroup sample values are substituted into sensibility equation to f (Q1,Q2)、g(Q1,Q2) carry out quadratic fit；It obtains and simplifies
Sensibility equation afterwards：；
Step 6 obtains heat source model parameter b by monofilament submerged-arc welding predictor formula1、b2Value, recycle simplification above
Sensibility prediction equation heat source model other parameters.
Further, double stripping mechanism is calculated using pattern search method in step 3, in finite element inversion to join
Number.
Further, in step 4, gained heat source model parameter and speed of welding are adjusted with 10% amplitude.
Further, increase modifying factor Δ x in sensibility equation.
Compared with prior art, the advantageous effect of the invention reached is：The present invention is considering two heat source interactions
Under the premise of, heat source density distribution is divided into three different regions of density function, being used as with the distribution of this density function has
The heat source load of limit member simulation, can improve the precision of numerical simulation.For untested welding condition, pass through prediction
Formula can be extended result, and make result of calculation serialization, so as to obtain arbitrary technique in a certain range
Heat source model parameter corresponding to parameter combination.The heat corresponding to different welding conditions can be directly obtained using the present invention
The tentative calculation workload of welding analog is substantially reduced, and the precision of simulation is substantially improved by source model parameter, reduces engineer testing,
Save development cost.
Description of the drawings
Fig. 1 is the flow diagram of the method for the present invention；
Fig. 2 is two heat source interaction schematic diagrames of tandem double wire hidden arc welding；
Fig. 3 is single, tandem double wire hidden arc welding model molten wide to 1.2 times of susceptibility test results of each parameter；
Fig. 4 is single, tandem double wire hidden arc welding model fusion penetration to 1.2 times of susceptibility test results of each parameter；
Fig. 5 is compared with experimental result using model parameter analog result obtained by predictor formula.
Specific implementation mode
The invention will be further described below in conjunction with the accompanying drawings.Following embodiment is only used for clearly illustrating the present invention
Technical solution, and not intended to limit the protection scope of the present invention.
As shown in Figure 1, a kind of tandem double wire hidden arc welding numerical simulation heat source model determination method for parameter of the present invention, packet
Include following steps：
Step 1：The double stripping mechanism that heat flow density decays along depth direction quadratic function is established, determines mariages mould
The heat flux distribution function of type；
In the prior art, the first half of double stripping mechanism setting welding pool is latter half of as 1/4 ellipsoid
It being allocated as another 1/4 ellipsoid, heat flow density is in Gaussian function normal distribution in semiellipsoid, and central part has maximum value,
Exponentially curve declines from center to edge.Through analysis, Pool and reality obtained by double stripping mechanism in the prior art
Welding pool pattern has larger difference in the depth direction.Therefore, the present invention is on the depth direction of double stripping mechanism
Attenuation function is added, fusion wire shaped is effectively controlled by this attenuation function, achievees the purpose that control weld pool shape, from
And melting pool shape obtained by finite element analysis and practical melting pool shape is made to have preferable coincide.Attenuation function generally includes a letter
Number, exponential function and quadratic function pass through a large amount of attenuation function tentative calculations for differential declines function, find molten obtained by quadratic function
Pond shape is similar to experiment gained melting pool shape.Specific modeling process is referring to application No. is 2016110116037 patents of invention
《A kind of modeling method of tandem double wire hidden arc welding numerical simulation heat source model》, therefore heat flow density is along the decaying letter of depth direction
Number selects the quadratic function, the heat flux distribution function of mariages model to be：
Heat flux distribution function is in first half ellipsoid：
Latter half ellipsoid heat flux distribution function is：
F in formulafAnd frRespectively total input power is in the distribution index of molten bath front and rear part, and ff+fr=2, QiIt is i-th
The electric arc effective power of root wire, wherein Q=η UI；U is weldingvoltage, and I is welding current, and η is electric arc effective thermal efficiency system
Number, value 0.77~0.9；V is speed of welding；αiFor the welding inclination angle of the i-th root wire, afi、ariIt is corresponding double ellipse for i-th silk
The length of the forward and backward hemisphere major semiaxis of ball heat source model, biFor the corresponding double semiminor axis of ellipsoid length of i-th silk, ciFor i-th silk
Corresponding double ellipsoid depth, in formulaPart is attenuation function (quadratic function), and t is to weld
The time of Cheng Jinhang, i=1,2.
It is found that weldingvoltage, welding current, speed of welding and welding inclination angle are known welding condition, determine
Double stripping mechanism shape needs to solve four parameter afi、ari、biAnd ci.And the size of this four parameters will directly influence
The distribution in temperature field in welding process, therefore, it is crucial to solve this four parameters.
Tandem double wire hidden arc welding heat flux distribution is divided into three by step 2 according to the interaction of forward and backward two heat sources
A region determines the heat flow density in each region；
It is simple overlaying relation not to be between the forward and backward silk heat source of tandem double wire hidden arc welding, but there are interactional
Relationship.Two heat sources of tandem double wire hidden arc welding interact schematic diagram as shown in Fig. 2, tandem double wire hidden arc welding heat flux distribution
It can be divided into three regions by two double forward and backward hemisphere lines of demarcation of ellipsoid：1 heat flow density of region is double ellipsoid first halfs of preceding silk
Divide and double ellipsoid first halfs of rear silk are superimposed；2 heat flow density of region be preceding silk double ellipsoid latter halfs and rear silk it is double ellipse
Ball first half is superimposed；3 heat flow density of region is double ellipsoid latter halfs superposition of the double ellipsoid latter halfs and rear silk of preceding silk.
These three regions are only referred to the region of heat flow density effect, because preceding silk and the mutual influence of rear silk, cause
Heat flow density is (can be understood as in these three Regional land surface heat fluxes difference) of different sizes.These three regions are by finite element mould
Intend the subprogram of software to write heat source program, to heat source model subregion is loaded on finite element model anti-to carry out
It drills.
Step 3 establishes finite element model, take the matched weldingvoltage of one group of submerged-arc welding, welding current, speed of welding and
Inclination angle is welded as simulation basic parameter, finite element inversion is carried out to the above heat source model parameter, obtains corresponding heat source mould
Shape parameter；
It being described in detail this process with an embodiment, preferred dimension is the Q345 steel plates of 18.4mm × 200mm × 300mm,
Three-dimensional welding finite element model is established in ABAQUS platforms carries out analog temperature field；Define temperature governing equation and perimeter strip
The weldment material properties parameter arrived involved in part equation, including the density of base material and weld seam, material phase transformation latent heat, convection coefficient,
The coefficient of heat conduction, specific heat capacity, radiation heat transfer coefficient, coefficient of thermal expansion, elasticity modulus and Poisson's ratio etc.；Set absolute zero and glass
The graceful constant of Wurz.
Non-uniform grid division is carried out to finite element model：It is molten to weld that unit size is used at Seam and heat effected zone
/ to eight/10ths of pond width, and use unit size to be arrived for 1/5th of welding pool width in base material periphery
A quarter.
Take the matched weldingvoltage of one group of submerged-arc welding, welding current, speed of welding as simulation basic parameter, before taking welding
Silk electric current is 600A, and preceding silk weldingvoltage is 36V, and rear silk welding current is 500A, and rear silk weldingvoltage is 40V, speed of welding
It for 3.5m/min, tests as a contrast, preceding silk aclinal, it is 15 ° that rear wire bond, which connects inclination angle, and forward and backward distance between weldingwires is 55mm, thermal effect
Rate coefficient η takes 0.9, ffTake 1.1, frTake 0.9.Experiment gained molten wide w=23.96mm, fusion penetration p=8.96mm, this experimental data can
It is directly measured using vernier caliper obtained by weldment.The above parameter is substituted into heat flux distribution formula, three regions are applied respectively
The heat source load for adding each region heat flow density determined above to describe；Inverting is carried out using pattern search method in the prior art,
Mentioned by this pattern search method also has in " accelerating step length inverting multiplewire submerged arc welding double stripping mechanism parameter " text,
Belong to the prior art, details are not described herein again.By inverting, an optimal group model parameter is obtained.To according to pattern search method into
Numerous parametric statistics conclusions obtained from row finite element modelling can obtain, and take afi=bi, ari=2.5bi, obtain model parameter
Respectively：
b1=8.6, c1=5.8, af1=8.6, ar1=21.5
b2=8.1, c2=5.4, af2=8.1, ar2=20.25
The fusion penetration molten wide in simulation molten bath can be obtained based on the above parameter, the fusion penetration simulated gained fusion penetration molten wide and actually measured
The data of molten wide are as shown in table 1, it is known that and simulation gained fusion penetration molten wide is differed with the fusion penetration molten wide data actually measured is less than 10%,
Therefore the above heat source model can be used for simulating the heat flux distribution of biserial mariages heat source.
1 analog result of table is compared with experimental result
Step 4 is adjusted heat source model parameter obtained as above and speed of welding with 10% amplitude, i.e., by afi、
ari、bi、ciIt is adjusted with v to the 0.8 of original parameter, 0.9,1.1,1.2 ... times, carrying out simulation to each parameter combination obtains molten bath
Molten wide, fusion penetration obtain the corresponding sample values of multigroup parameter and molten wide, fusion penetration.Sensitivity analysis is carried out to each group sample values.
Compare the parameter of the respective heat source model of monofilament submerged-arc welding, mariages (forward and backward silk) submerged-arc welding to the shadow of fusion penetration, molten wide
It rings, to obtain the Heat-Source Parameters being affected, carries out sensitivity analysis.The heat source model of monofilament is also double-ellipsoid heat source mould
Type.The sensitivity analysis of monofilament parameter is referring to application No. is 2016110210664 patents of invention《A kind of monofilament submerged-arc welding numerical value
Simulate heat source model determination method for parameter》.Monofilament heat source model parameter and mariages parameter comparison are carried out herein, are in order to for after
The parameter operation of formula provides foundation after the simplification in face.
It is preceding silk or rear silk parameter list when comparing the influence of preceding silk or the individual heat source model parameter of rear silk to fusion penetration molten wide
Solely variation, i.e., before silk parameter adjustment when, rear silk parameter remains unchanged, and when rear silk parameter adjustment, preceding silk parameter remains unchanged.
When each parameter adjustment of model is to 1.2 times, silk model after silk, double wire hidden arc welding before monofilament submerged-arc welding, double wire hidden arc welding
In influence of each parameter to molten wide as shown in figure 3, before monofilament submerged-arc welding, double wire hidden arc welding after silk, double wire hidden arc welding in silk model
Influence of each parameter to fusion penetration is as shown in Figure 4, it is known that parameter afiAnd ariInfluence to molten wide fusion penetration is smaller, double wire hidden arc welding b1、
b2Influence to molten wide or fusion penetration is close.With double stripping mechanism parameter：b1、b2、c1、c2It is input parameter with speed of welding v,
Using molten wide w, fusion penetration p as output parameter, corresponding sample is substituted into regression equation and obtains sensitivity analysis result：
Wherein, x2w、x3w、x4w、x5w、x6wFor b1、c1、b2、c2, v is to the sensitivity coefficient of molten wide w；x2p、x3p、x4p、x5p、
x6pFor b1、c1、b2、c2, v is to the sensitivity coefficient of molten wide p；x1w、x1pFor with Q1、Q2(Q1、Q2Respectively forward and backward silk electric arc is effective
Power) related function, take x1w、x1pFor Q1、Q2(electric arc effective power) related quadratic function f (Q1,Q2)、g(Q1,Q2)。
Step 5 is simplified and is fitted to gained sensibility equation；
Multigroup sample values are substituted into formula (3) (4) to be fitted, sensibility equation (5), (6) can be obtained：
Wherein Q1、Q2Respectively forward and backward silk electric arc effective power, Δ x is modifying factor, due to sensibility equation ignored
Af、arInfluence of the parameter to fusion penetration p, molten wide w, casts out parameter afi、ariIndex coefficient, in order to increase the accurate of formula
Property, part will be cast out and be approximately equal to modifying factor Δ x.
Take x1w、x1pFor Q1、Q2(electric arc effective power) related quadratic function f (Q1,Q2)、g(Q1,Q2), quadratic function
Expression formula is as follows：
f(Q1,Q2)=α0+α1Q1+α2Q2+α3Q1Q2+α4Q1 2+α5Q2 2 (7)
g(Q1,Q2)=β0+β1Q1+β2Q2+β3Q1Q2+β4Q1 2+β5Q2 2 (8)
Wherein, αj, βj(j=0,1,2,3,4,5) it is pre- fitting constant.Q can be effectively indicated in quadratic function1, Q2It is mutual
Influence relationship, and Q in quadratic function formula1, Q2Output parameter is influenced to be peer-to-peer, is tallied with the actual situation.
Multigroup sample values are substituted into sensibility equation to f (Q1,Q2)、g(Q1,Q2) carry out quadratic fit；It obtains and simplifies
Sensibility equation afterwards：
Wherein △ x are fitted to by multigroup sample valuesApproximation obtains 0.98.
The simplification sensibility equation obtained is used for heat source model parameter prediction by step 6, will be pre- using ABAQUS softwares
Parameter is surveyed for heat source model and loading, obtains temperature field result, weld pool shape and corresponding defeated is obtained by temperature field result
Go out parameter molten wide w, fusion penetration p as a result, and experiment show.
Silk electric current is 615A before taking welding, and preceding silk weldingvoltage is 34V, and rear silk welding current is 480A, and rear wire bond connects electricity
Pressure is 40V, and speed of welding 3.5m/min, preceding silk welding inclination angle, it is 15 ° that rear wire bond, which connects inclination angle, and forward and backward distance between weldingwires is
55mm, as confirmatory experiment, experiment measures weldment molten bath molten wide w=23.58mm, fusion penetration p=10.31mm, and according to fig. 3, Fig. 4
It can be seen that silk parameter b before double wire hidden arc welding1, rear silk parameter b2It is similar to the affecting laws of output parameter w, p, b1、b2Value
It approximate can be obtained by monofilament submerged-arc welding predictor formula.The specific solution procedure of monofilament submerged-arc welding predictor formula referring to application No. is
2016110210664 patent of invention《A kind of monofilament submerged-arc welding numerical simulation heat source model determination method for parameter》.First test
The molten wide of silk, fusion penetration numerical value before measuring, by preceding silk parameter b1, welding current, weldingvoltage, speed of welding substitute into the molten of monofilament
In wide, two predictor formulas of fusion penetration, parameter b is calculated1；Then the molten wide of silk, fusion penetration numerical value, the silk by after after experiment is measured
Parameter b2, welding current, weldingvoltage, speed of welding substitute into two molten wide of monofilament, fusion penetration predictor formulas, calculate parameter
b2。
Again by calculated parameter b1And b2It substitutes into predictor formula (9) (10) and calculates parameter c1And c2, take afi=bi,ari
=2.5bi, thus show that model parameter is：
b1=8.5, c1=3.9, af1=8.5, ar1=21.3
b2=8.0, c2=7.9, af2=8.0, ar2=20.
The above parameter is substituted into double stripping mechanism, temperature field and molten bath shape are obtained with the simulation of ABAQUS softwares
Shape.It is as shown in table 2 to simulate obtained molten wide, fusion penetration numerical value and experiment gained molten wide, fusion penetration numerical value, experimental result and predictor formula
The analog result error only 2.5% or so of gained model parameter illustrates that the two goodness of fit is higher.Simulate obtained melting pool shape with
For experiment gained melting pool shape comparison diagram as shown in figure 5, the left side is experiment gained melting pool shape in figure, the right is simulation gained molten bath
Shape graph as can be seen from the figure also has with experimental result on simulation gained Pool and preferably coincide, and predictor formula obtains
The verification of experiment.
Analog result obtained by 2 predictor formula of table is compared with experimental result
The present invention is divided into density function difference under the premise of considering two heat source interactions, by heat source density distribution
Three regions, the heat source load of finite element modelling is used as with the distribution of this density function, the precision of numerical simulation can be improved.
For untested welding condition, result can be extended by predictor formula, and makes result of calculation serialization,
So as to obtain the heat source model parameter corresponding to arbitrary combination of process parameters in a certain range.It can using the present invention
The heat source model parameter corresponding to different welding conditions is directly obtained, the tentative calculation workload of welding analog is substantially reduced,
And the precision of simulation is substantially improved, engineer testing is reduced, development cost is saved.
The above is only a preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For member, without departing from the technical principles of the invention, several improvements and modifications, these improvements and modifications can also be made
Also it should be regarded as protection scope of the present invention.
Claims (4)
1. a kind of tandem double wire hidden arc welding numerical simulation heat source model determination method for parameter, characterized in that include the following steps：
Step 1：The double stripping mechanism that heat flow density decays along depth direction quadratic function is established, determines mariages model
Heat flux distribution function；
Attenuation function is quadratic function, and heat flux distribution function is：
Heat flux distribution function is in first half ellipsoid：
Latter half ellipsoid heat flux distribution function is：
F in formulafiAnd friThe total input power of respectively i-th silk is in the distribution index of molten bath front and rear part, and ffi+fri=2,
QiFor the electric arc effective power of i-th silk, v is speed of welding；αiFor the welding inclination angle of i-th silk, afi、ariIt is right for i-th silk
The length for the forward and backward hemisphere major semiaxis of double stripping mechanism answered, biFor the corresponding double semiminor axis of ellipsoid length of i-th silk, ciFor
The corresponding double ellipsoid depth of i-th silk, in formulaPart is that quadratic function is decayed, and t is to weld
The time of Cheng Jinhang, i=1,2；
Tandem double wire hidden arc welding heat flux distribution is divided into three areas by step 2 according to the interaction of forward and backward two heat sources
Domain determines the heat flow density in each region；
Tandem double wire hidden arc welding heat flux distribution can be divided into three regions by two double forward and backward hemisphere lines of demarcation of ellipsoid：Area
1 heat flow density of domain is double ellipsoid first halfs superposition of the double ellipsoid first halfs and rear silk of preceding silk；Before 2 heat flow density of region is
Double ellipsoid first halfs superposition of the double ellipsoid latter halfs and rear silk of silk；3 heat flow density of region is that double ellipsoids of preceding silk are later half
Double ellipsoid latter halfs of part and rear silk are superimposed；
Step 3 establishes finite element model, takes the matched weldingvoltage of one group of submerged-arc welding, welding current, speed of welding as mould
Quasi- basic parameter applies the heat source load of each region heat flow density description determined above, is carried out to the above heat source model parameter
Finite element inversion obtains corresponding heat source model parameter；
Step 4 is adjusted heat source model parameter obtained as above and speed of welding with the amplitude set, is based on each parameter group
Conjunction obtains corresponding molten wide, fusion penetration sample values；With double stripping mechanism parameter：b1、b2、c1、c2It is defeated with speed of welding v
Enter parameter, using molten wide w, fusion penetration p as output parameter, corresponding sample is substituted into regression equation and obtains sensitivity analysis result：
Wherein, x2w、x3w、x4w、x5w、x6wFor b1、c1、b2、c2, v is to the sensitivity coefficient of molten wide w；x2p、x3p、x4p、x5p、x6pFor
b1、c1、b2、c2, v is to the sensitivity coefficient of molten wide p；x1w、x1pFor with Q1、Q2Related quadratic function f (Q1,Q2)、g(Q1,Q2)；
Step 5 is fitted and simplifies to gained sensibility equation, and detailed process is：
Take x1w、x1pFor Q1、Q2Related quadratic function f (Q1,Q2)、g(Q1,Q2), the expression formula of quadratic function is as follows：
f(Q1,Q2)=α0+α1Q1+α2Q2+α3Q1Q2+α4Q1 2+α5Q2 2
g(Q1,Q2)=β0+β1Q1+β2Q2+β3Q1Q2+β4Q1 2+β5Q2 2
Wherein, αj, βjFor pre- fitting constant, j=0,1,2,3,4,5；
Multigroup sample values are substituted into sensibility equation to f (Q1,Q2)、g(Q1,Q2) carry out quadratic fit；It obtains after simplifying
Sensibility equation；
Step 6 obtains heat source model parameter b by monofilament submerged-arc welding predictor formula1、b2Value, recycle simplification above sensitive
Property prediction equation heat source model other parameters.
2. a kind of tandem double wire hidden arc welding numerical simulation heat source model determination method for parameter according to claim 1,
It is characterized in, in step 3, double stripping mechanism parameter is calculated using pattern search method in finite element inversion.
3. a kind of tandem double wire hidden arc welding numerical simulation heat source model determination method for parameter according to claim 1,
It is characterized in, in step 4, heat source model parameter obtained as above and speed of welding is adjusted with 10% amplitude.
4. a kind of tandem double wire hidden arc welding numerical simulation heat source model determination method for parameter according to claim 1,
It is characterized in, in step 5, increases in sensibility equation and be equal to parameter afi、ariThe modifying factor Δ x of index coefficient.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611064245.6A CN106363283B (en) | 2016-11-28 | 2016-11-28 | A kind of tandem double wire hidden arc welding numerical simulation heat source model determination method for parameter |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611064245.6A CN106363283B (en) | 2016-11-28 | 2016-11-28 | A kind of tandem double wire hidden arc welding numerical simulation heat source model determination method for parameter |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106363283A CN106363283A (en) | 2017-02-01 |
| CN106363283B true CN106363283B (en) | 2018-11-09 |
Family
ID=57892372
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611064245.6A Expired - Fee Related CN106363283B (en) | 2016-11-28 | 2016-11-28 | A kind of tandem double wire hidden arc welding numerical simulation heat source model determination method for parameter |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106363283B (en) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107598370B (en) * | 2017-08-28 | 2018-11-30 | 温州大学 | A kind of technique optimization method of steel/aluminium laser welding |
| CN107649804B (en) * | 2017-10-17 | 2023-06-23 | 华中科技大学鄂州工业技术研究院 | Online detection and control system for penetration of additive manufacturing |
| CN108319740A (en) * | 2017-12-04 | 2018-07-24 | 吉林亚新工程检测有限责任公司 | The vertical bulk heat treatmet Numerical Model of Temperature Field modeling method of pressure vessel internal combustion method |
| CN108681644A (en) * | 2018-05-21 | 2018-10-19 | 河海大学常州校区 | A method of the prediction welding heat affected sector width of double wire hidden arc welding |
| CN109190322B (en) * | 2018-11-07 | 2023-03-24 | 桂林电子科技大学 | Electron beam cladding process parameter optimization method and system based on temperature field |
| CN110705159A (en) * | 2019-09-26 | 2020-01-17 | 华中科技大学 | Heat source model parameter solving method, device, equipment and storage medium |
| CN112276313A (en) * | 2020-10-19 | 2021-01-29 | 上海振华重工(集团)股份有限公司 | Method for predicting hot and cold multi-wire composite submerged arc welding thermal cycle parameters of large steel structural part |
| CN116011221B (en) * | 2022-12-31 | 2023-09-29 | 华中科技大学 | Method and system for rapidly checking welding heat source model parameters based on welding morphology |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0333195A1 (en) * | 1988-03-18 | 1989-09-20 | Hitachi, Ltd. | Method and Apparatus for Automatic Welding Control |
| JP2010221302A (en) * | 2010-06-16 | 2010-10-07 | National Institute For Materials Science | Arc welding method |
| CN102637235A (en) * | 2012-05-02 | 2012-08-15 | 中国石油集团渤海石油装备制造有限公司 | Determination method for heat source model parameters in multiplewire submerged-arc welding by numerical simulation |
| CN103246774A (en) * | 2013-05-13 | 2013-08-14 | 天津大学 | Numerical simulation method for P92 steel tube welding heat affected zone |
| CN105975708A (en) * | 2016-05-16 | 2016-09-28 | 中国计量大学 | Steel tube welding parameter optimization method based on numerical simulation and data analysis |
-
2016
- 2016-11-28 CN CN201611064245.6A patent/CN106363283B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0333195A1 (en) * | 1988-03-18 | 1989-09-20 | Hitachi, Ltd. | Method and Apparatus for Automatic Welding Control |
| JP2010221302A (en) * | 2010-06-16 | 2010-10-07 | National Institute For Materials Science | Arc welding method |
| CN102637235A (en) * | 2012-05-02 | 2012-08-15 | 中国石油集团渤海石油装备制造有限公司 | Determination method for heat source model parameters in multiplewire submerged-arc welding by numerical simulation |
| CN103246774A (en) * | 2013-05-13 | 2013-08-14 | 天津大学 | Numerical simulation method for P92 steel tube welding heat affected zone |
| CN105975708A (en) * | 2016-05-16 | 2016-09-28 | 中国计量大学 | Steel tube welding parameter optimization method based on numerical simulation and data analysis |
Non-Patent Citations (2)
| Title |
|---|
| 基于SYSWELD的多丝埋弧直缝焊管三维热过程数值模拟研究;赵波 等;《焊管》;20120331;第35卷(第3期);41-46 * |
| 窄间隙埋弧焊温度场数值分析;张磊 等;《焊接学报》;20160331;第37卷(第3期);83-89 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106363283A (en) | 2017-02-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106363283B (en) | A kind of tandem double wire hidden arc welding numerical simulation heat source model determination method for parameter | |
| CN106529051B (en) | A kind of monofilament submerged-arc welding numerical simulation heat source model determination method for parameter | |
| Gunaraj et al. | Application of response surface methodology for predicting weld bead quality in submerged arc welding of pipes | |
| Dhas et al. | Optimization of parameters of submerged arc weld using non conventional techniques | |
| Nie et al. | Experimental study and modeling of H13 steel deposition using laser hot-wire additive manufacturing | |
| Murugan et al. | Prediction and control of weld bead geometry and shape relationships in submerged arc welding of pipes | |
| Karaoğlu et al. | Sensitivity analysis of submerged arc welding process parameters | |
| Juang et al. | Process parameter selection for optimizing the weld pool geometry in the tungsten inert gas welding of stainless steel | |
| Moysan et al. | Modelling the grain orientation of austenitic stainless steel multipass welds to improve ultrasonic assessment of structural integrity | |
| CN104298817B (en) | A kind of bielliptic(al) exponential damping body heat source model for being used to simulate high energy beam welding | |
| Vasudevan et al. | Genetic-algorithm-based computational models for optimizing the process parameters of A-TIG welding to achieve target bead geometry in type 304 L (N) and 316 L (N) stainless steels | |
| Little et al. | The effect of thermal properties and weld efficiency on transient temperatures during welding | |
| CN106529047B (en) | A kind of modeling method of tandem double wire hidden arc welding numerical simulation heat source model | |
| Perić et al. | Numerical prediction and experimental validation of temperature and residual stress distributions in buried‐arc welded thick plates | |
| JP2010201474A (en) | Method and system for optimization of welding, and welding method | |
| CA2383569C (en) | Method for providing an analytical solution for a thermal history of a welding process | |
| Sreeraj et al. | Simulation and parameter optimization of GMAW process using neural networks and particle swarm optimization algorithm | |
| CN108681644A (en) | A method of the prediction welding heat affected sector width of double wire hidden arc welding | |
| Sharma et al. | A practical approach towards mathematical modeling of deposition rate during twin-wire submerged arc welding | |
| Vasudevan et al. | Genetic algorithm for optimisation of A-TIG welding process for modified 9Cr–1Mo steel | |
| CN113139314A (en) | Heat source numerical simulation method for laser additive manufacturing process | |
| CN108581165A (en) | A kind of oxygen-free copper welding parameter of electron beam prediction computational methods | |
| CN106529129A (en) | Cladding layer cross section contour curve under broadband laser effect and modeling method | |
| Sreeraj et al. | Optimization of GMAW process parameters using particle swarm optimization | |
| Frank et al. | An investigation of cause-and-effect relationships within a 3D-Printing system and the applicability of optimum printing parameters from experimental models to different printing jobs |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | 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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181109 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |