CN109522610A - A kind of weld metal zone mixing crystalline region finite element modeling method - Google Patents
A kind of weld metal zone mixing crystalline region finite element modeling method Download PDFInfo
- Publication number
- CN109522610A CN109522610A CN201811251317.7A CN201811251317A CN109522610A CN 109522610 A CN109522610 A CN 109522610A CN 201811251317 A CN201811251317 A CN 201811251317A CN 109522610 A CN109522610 A CN 109522610A
- Authority
- CN
- China
- Prior art keywords
- area
- finite element
- crystal
- crystalline region
- coordinate
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention provides the methods that a kind of pair of weld metal zone mixing crystalline region carries out finite element modeling, using improved Thiessen polygon method, Excel is combined to carry out coordinate integration using Matlab programming, the intersecting point coordinate of the mixing crystalline region crystal grain center position coordinates comprising column crystal and different grain size equiaxed grain structures and crystal boundary can be obtained, to set up the geometrical model of mixing crystalline region.The intersecting point coordinate of crystal grain center position coordinates and crystal boundary is exported, text file is generated and is read for finite element software Abaqus；Text file is read using the Python script that Abaqus operation writes, establishes the finite element model that mix-crystal district's groups are knitted.The present invention considers that welding point microstructure is non-homogeneous, can establish the finite element model of weld seam difference microcell, and realize and different crystalline regions are assigned with different material parameters, keep analog result relatively reliable；The method that the modeling of mixing crystalline region is combined using Matlab program with Python script, eliminates cumbersome Abaqus interface operation, can be realized across the limited dimension member simulation of welding point and calculates.
Description
Technical field
The present invention relates to metal welding joints simulations and simulation technical field more particularly to a kind of weld metal zone mixing crystalline region to have
Limit Meta Model method.
Background technique
In order to material mechanical property carry out micro-scale research, consider the microcosmic grain structure characteristic of material across ruler
Degree simulation and emulation are more and more.Therefore in order to improve the feasibility and accuracy of simulation, foundation meets Fine Texture of Material spy
The finite element model of property seems most important.
Common melting welding head weld metal zone microstructure is uneven, generally includes column crystal, thick equiax crystal and thin isometric
The mix-crystals tissue such as crystalline substance.Traditional Thiessen polygon method is only used for establishing the physics mould of the equiaxed grain structure of single grain size
Type is not suitable for establishing the weld metal zone finite element model comprising column crystal and different grain size equiax crystal.Existing electronic diffraction
Figure method (EBSD) combines finite element software, although complicated microcosmic grain structure finite element image can be obtained, it exists as follows not
Foot: it though 1, EBSD technology has been able to realize full-automatic acquisition microcell orientation information, can not achieve from microstructure modeling to splitting
The Whole Process Simulation of line germinating, and subsequent processing is cumbersome.2, the equipment cost is higher, and the difficulty that popularity uses is larger.
Summary of the invention
Technical problem to be solved by the present invention lies in provide a kind of to calculate at low cost, versatile mixing crystalline region modeling
Method.
The present invention is to solve above-mentioned technical problem by the following technical programs:
A kind of weld metal zone mixing crystalline region finite element modeling method, comprising the following steps:
Step 1: choose the actual welds area of certain size, obtain fusion area, heat affected area and the size in base material area and
The mean grain size of column crystal mean grain size, heat affected area and the equiax crystal in base material area in fusion area, and calculated with this
Column crystal quantity x in the fusion area and equiax crystal quantity z in equiax crystal quantity y and base material area in heat affected area；
Step 2: selected in MATLAB dimension scale and the identical region of actual welds area ratio as simulation weld metal zone,
Weld metal zone, which will be simulated, according to actual ratio is transversely divided into simulation fusion area, simulated heat-affected zone and simulation base material area；
Step 3: generating x abscissa spacing at random in the simulation fusion area of step 2 respectively greater than ordinate spacing
Coordinate points generate y, z coordinate points in simulated heat-affected zone and simulation base material area, by the coordinate of all coordinate points at random respectively
Export；
Step 4: reading derived point coordinate in MATLAB, use Thiessen polygon method using coordinate as crystal grain center
Coordinate obtains the physical model of mixing crystalline region, and crystal boundary intersecting point coordinate is exported；
Step 5: reading crystal boundary intersecting point coordinate in Finite Element Simulation Software, obtain the finite element model of mixing crystalline region.
Preferably, in step 1 in specific region equiax crystal quantity P calculation method are as follows:
Wherein, S is the area of corresponding region, and a is the grain size of equiax crystal in corresponding region.
Preferably, the abscissa that adjacent coordinates point in fusion area is simulated in step 3 is about three times of ordinate.
The advantages of weld metal zone mixing provided by the invention crystalline region finite element modeling method, is: can be used in establishing comprising more
The seam organization finite element model of kind crystal, the material mechanical performance knitted for simulating mix-crystal district's groups realize across size simulation
With emulation, and large-scale experiment equipment is not needed, save the cost, using the mixing crystalline region finite element model of resume of the present invention, to column
Shape crystalline substance and equiax crystal impart different material parameters respectively, and analog result is relatively reliable.
Detailed description of the invention
Fig. 1 is the schematic diagram of weld metal zone mixing crystalline region finite element modeling method provided by the embodiment of the present invention；
Fig. 2 is electron beam welded joint metallographic structure figure provided by the embodiment of the present invention；
Fig. 3 is that physical model schematic diagram in crystalline region is mixed provided by the embodiment of the present invention；
Fig. 4 is the schematic diagram of mixing crystalline region finite element model provided by the embodiment of the present invention.
Specific embodiment
To make the objectives, technical solutions, and advantages of the present invention clearer, below in conjunction with specific embodiment, and reference
Attached drawing, the present invention is described in further detail.
As shown in Figure 1, a kind of weld metal zone mixing crystalline region finite element modeling method, comprising the following steps:
Step 1: choose the actual welds area of certain size, obtain fusion area, heat affected area and the size in base material area and
The mean grain size of column crystal mean grain size, heat affected area and the equiax crystal in base material area in fusion area, and calculated with this
Column crystal quantity x in the fusion area and equiax crystal quantity z in equiax crystal quantity y and base material area in heat affected area；
For titanium alloy electron beam welded joint metallographic structure as shown in Figure 2, the present embodiment is classified as symmetrical along y-axis
Two models, and region of the side selection wherein having a size of 0.5mm × 0.5mm is studied, and along weld metal zone, laterally it is molten
Area's range substantially 0~0.2mm is closed, the column crystal catercorner length in fusion area is about 0.2mm；Heat affected area range is substantially
0.2mm~0.4mm, the equiax crystal grain size in heat affected area is about 20 μm；Base material area range substantially 0.4mm~0.5mm, it is female
Equiax crystal grain size in material area is about 40 μm.Then the column crystal quantity in fusion area is about x=6,
The calculation method of equiax crystal quantity P in specific region are as follows:
Wherein, S is the area of corresponding region, and a is the grain size of equiax crystal in corresponding region.
Specific in the present embodiment, the equiax crystal quantity y=250 of heat affected area, the equiax crystal quantity z=31 in base material area.
Step 2: selected in MATLAB dimension scale and the identical region of actual welds area ratio as simulation weld metal zone,
Weld metal zone, which will be simulated, according to actual ratio is transversely divided into simulation fusion area, simulated heat-affected zone and simulation base material area；
Step 3: generating x abscissa spacing at random in the simulation fusion area of step 2 respectively greater than ordinate spacing
Coordinate points, in preferred embodiment, the abscissa of adjacent coordinates point is about three times of ordinate；It is female in simulated heat-affected zone and simulation
It generates y, z coordinate points in material area at random respectively, the coordinate of all coordinate points is exported as into EXCEL file；
Step 4: reading derived point coordinate in MATLAB, use Thiessen polygon method to simulate the point coordinate of generation work
The physical model of mixing crystalline region shown in Fig. 3 is obtained for crystal grain centre coordinate, and is in text file by the export of crystal boundary intersecting point coordinate；
Step 5: reading crystal boundary intersecting point coordinate in Finite Element Simulation Software, obtain the limited of mixing crystalline region shown in Fig. 4
Meta-model.
In preferred embodiment, finite element software selects Abaqus software, and running python script wherein keeps its identification brilliant
The specific procedure of boundary's coordinate and the schematic diagram for generating finite element model, selection is conventional technical means, and details are not described herein again.
Particular embodiments described above has carried out further in detail the purpose of the present invention, technical scheme and beneficial effects
It describes in detail bright, it should be understood that the above is only a specific embodiment of the present invention, is not intended to restrict the invention, not
It is any modification that those of ordinary skill in the art make the present invention, equivalent under the premise of being detached from the spirit and principles in the present invention
Replacement, improvement etc., should all fall within the protection scope that claims of the present invention determines.
Claims (3)
1. a kind of weld metal zone mixing crystalline region finite element modeling method, it is characterised in that: the following steps are included:
Step 1: choosing the actual welds area of certain size, obtain the size and fusion of fusion area, heat affected area and base material area
The mean grain size of column crystal mean grain size, heat affected area and the equiax crystal in base material area in area, and fusion is calculated with this
The equiax crystal quantity z in the equiax crystal quantity y and base material area in column crystal quantity x, heat affected area in area；
Step 2: selected in MATLAB dimension scale and the identical region of actual welds area ratio as simulation weld metal zone, according to
Actual ratio will simulate weld metal zone and transversely be divided into simulation fusion area, simulated heat-affected zone and simulation base material area；
Step 3: generating the coordinate that x abscissa spacing is greater than ordinate spacing at random in the simulation fusion area of step 2 respectively
Point generates y, z coordinate points in simulated heat-affected zone and simulation base material area at random respectively, the coordinate of all coordinate points is led
Out；
Step 4: reading derived point coordinate in MATLAB, use Thiessen polygon method using coordinate as crystal grain centre coordinate
The physical model of mixing crystalline region is obtained, crystal boundary intersecting point coordinate is exported；
Step 5: reading crystal boundary intersecting point coordinate in Finite Element Simulation Software, obtain the finite element model of mixing crystalline region.
2. a kind of weld metal zone mixing crystalline region according to claim 1 finite element modeling method, it is characterised in that: in step 1
The calculation method of equiax crystal quantity P in specific region are as follows:
Wherein, S is the area of corresponding region, and a is the grain size of equiax crystal in corresponding region.
3. a kind of weld metal zone mixing crystalline region according to claim 1 finite element modeling method, it is characterised in that: in step 3
The abscissa of simulation fusion area adjacent coordinates point is about three times of ordinate.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811251317.7A CN109522610A (en) | 2018-10-25 | 2018-10-25 | A kind of weld metal zone mixing crystalline region finite element modeling method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811251317.7A CN109522610A (en) | 2018-10-25 | 2018-10-25 | A kind of weld metal zone mixing crystalline region finite element modeling method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109522610A true CN109522610A (en) | 2019-03-26 |
Family
ID=65774240
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811251317.7A Pending CN109522610A (en) | 2018-10-25 | 2018-10-25 | A kind of weld metal zone mixing crystalline region finite element modeling method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109522610A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110400362A (en) * | 2019-07-12 | 2019-11-01 | 重庆大学 | A kind of ABAQUS two dimension crack modeling method, system and computer readable storage medium based on image |
| CN110516378A (en) * | 2019-08-30 | 2019-11-29 | 桂林电子科技大学 | A kind of Finite Element Method that tinbase binary eutectic mutually separates |
| CN111259585A (en) * | 2020-01-14 | 2020-06-09 | 大连交通大学 | Simulation method and system for weld grain nucleation growth |
| CN111881604A (en) * | 2020-07-24 | 2020-11-03 | 西安建筑科技大学 | Three-dimensional finite element model modeling method for Thiessen polygon subdivision |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103761374A (en) * | 2014-01-10 | 2014-04-30 | 上海工程技术大学 | Method for building laser deep penetration welding heat source model |
| CN104573141A (en) * | 2013-10-10 | 2015-04-29 | 上海汽车集团股份有限公司 | Weld joint connection simulation method based on finite element analysis |
| CN107358005A (en) * | 2017-07-24 | 2017-11-17 | 哈尔滨工业大学 | FEM model modeling method based on true heterogeneous microstructure SEM EBSD images |
-
2018
- 2018-10-25 CN CN201811251317.7A patent/CN109522610A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104573141A (en) * | 2013-10-10 | 2015-04-29 | 上海汽车集团股份有限公司 | Weld joint connection simulation method based on finite element analysis |
| CN103761374A (en) * | 2014-01-10 | 2014-04-30 | 上海工程技术大学 | Method for building laser deep penetration welding heat source model |
| CN107358005A (en) * | 2017-07-24 | 2017-11-17 | 哈尔滨工业大学 | FEM model modeling method based on true heterogeneous microstructure SEM EBSD images |
Non-Patent Citations (3)
| Title |
|---|
| JUNGHO CHO: "Voronoi diagram and microstructure of weldment", 《JOURNAL OF MECHANICAL SCIENCE AND TECHNOLOGY》 * |
| 余明敏: "基于元胞自动机法的TC4-DT钛合金电子束焊接头晶粒组织模拟", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
| 陈小进: "TC4-DT钛合金电子束焊接接头裂纹萌生数值模拟及试验研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110400362A (en) * | 2019-07-12 | 2019-11-01 | 重庆大学 | A kind of ABAQUS two dimension crack modeling method, system and computer readable storage medium based on image |
| CN110400362B (en) * | 2019-07-12 | 2023-02-24 | 重庆大学 | ABAQUS two-dimensional crack modeling method and system based on image and computer readable storage medium |
| CN110516378A (en) * | 2019-08-30 | 2019-11-29 | 桂林电子科技大学 | A kind of Finite Element Method that tinbase binary eutectic mutually separates |
| CN110516378B (en) * | 2019-08-30 | 2022-08-09 | 桂林电子科技大学 | Finite element simulation method for tin-based binary eutectic phase separation |
| CN111259585A (en) * | 2020-01-14 | 2020-06-09 | 大连交通大学 | Simulation method and system for weld grain nucleation growth |
| CN111881604A (en) * | 2020-07-24 | 2020-11-03 | 西安建筑科技大学 | Three-dimensional finite element model modeling method for Thiessen polygon subdivision |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109522610A (en) | A kind of weld metal zone mixing crystalline region finite element modeling method | |
| CN108399307A (en) | A kind of laser 3D printing Finite Element Method | |
| Toyserkani et al. | 3-D finite element modeling of laser cladding by powder injection: effects of laser pulse shaping on the process | |
| CN108022236B (en) | The welding stress deformation prediction technique and system of fitting are extracted based on weld image | |
| Pavlyk et al. | Simulation of weld solidification microstructure and its coupling to the macroscopic heat and fluid flow modelling | |
| CN106057027A (en) | Virtual testing and inspection of a virtual weldment | |
| CN104298817B (en) | A kind of bielliptic(al) exponential damping body heat source model for being used to simulate high energy beam welding | |
| CN110489821B (en) | Nickel-based alloy surfacing molten pool dendritic crystal growth numerical simulation method | |
| CN106363283B (en) | A kind of tandem double wire hidden arc welding numerical simulation heat source model determination method for parameter | |
| CN106529051A (en) | Method for determining heat source model parameters of single wire submerged arc welding numerical simulation | |
| CN109420819A (en) | A kind of robot increasing material manufacturing method based on MIG/MAG welding procedure | |
| CN109359337B (en) | Simulation method for variable light spot and powder particle size coupling effect | |
| CN105138772A (en) | Finite element modeling method for electron beam welding of variable-section component | |
| Huang et al. | The transient behaviours of free surface in a fully penetrated weld pool in gas tungsten arc welding | |
| CN101719333A (en) | Grid realization method for simulation and emulation of plate surfacing forming process | |
| Bi et al. | Anisotropic microstructure, properties and molecular dynamics simulation of CoCrNi medium entropy alloy fabricated by laser directed energy deposition | |
| CN113284570A (en) | Simulation method of microstructure of aluminum alloy welding pool | |
| Seufzer | Additive Manufacturing Modeling and Simulation A Literature Review for Electron Beam Free Form Fabrication | |
| CN111581820B (en) | Novel simulation method for melting additive manufacturing process of laser area array selected area | |
| Bauer et al. | Numerical simulation and calibration of a single seam WAAM process with a commercial and an open source software | |
| Tong et al. | Impacts of torch moving on phase change and fluid flow in weld pool of SMAW | |
| CN109885946B (en) | Method for determining energy distribution of composite heat source and welding simulation method | |
| Yang et al. | Numerical study of heat transfer characteristics in positional wire and arc additive manufacturing | |
| Han et al. | Correlation research on multi-field coupling of multi-track overlapping laser cladding process for nodular cast iron | |
| CN114472924B (en) | Laser path planning method, system, computer and readable storage medium |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190326 |
|
| RJ01 | Rejection of invention patent application after publication |