CN106326574B - It is a kind of meter and welding deformation and residual stress ultimate strength numerical computation method - Google Patents
It is a kind of meter and welding deformation and residual stress ultimate strength numerical computation method Download PDFInfo
- Publication number
- CN106326574B CN106326574B CN201610741183.1A CN201610741183A CN106326574B CN 106326574 B CN106326574 B CN 106326574B CN 201610741183 A CN201610741183 A CN 201610741183A CN 106326574 B CN106326574 B CN 106326574B
- Authority
- CN
- China
- Prior art keywords
- welding
- analysis
- residual stress
- ultimate strength
- welding deformation
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention discloses a kind of meter and the ultimate strength numerical computation methods of welding deformation and residual stress, it is related to the method for numerical simulation technical field of thin-wall construction Ultimate Strength Calculation；Its analogy method are as follows: tested by material mechanical performance, determine thin-wall construction material parameter；The welding condition of query structure writes welding subprogram；Temperature-displacement coupling analysis model is established, welding analog analysis is carried out, obtains welding deformation and residual stress；Definition restarts analysis on the basis of temperature-displacement coupling analysis model, carries out the analysis of feature buckling mode；Static arc-length methods analysis model is established, will be deformed comprising the low order feature buckling mode including welding deformation and welding residual stress is introduced into progress structural limits Strength co-mputation in arc-length methods analysis model.It not only can efficiently but also accurately calculate the ultimate strength of the thin-wall construction of meter and welding deformation and residual stress, all be of great importance to scientific research and engineering design.
Description
Technical field
The present invention relates to a kind of meter and the ultimate strength numerical computation method of welding deformation and residual stress, particularly one
Method for numerical simulation that welding deformation and residual stress introduced to thin-wall construction Ultimate Strength Calculation of the kind based on ABAQUS.
Background technique
Unstable phenomenon can occur under certain load for thin-wall construction, such as by the pressurized spherical shell of uniform external, by axial in face
The stiffened panel structure etc. of compressive load effect.When with the stability of Analysis on Numerical Simulation Method thin-wall construction, need to consider to tie
Initial deformation and residual stress in structure, currently, the initial deformation of structure is mainly introduced by two methods: first is that using measuring
Equipment measures the practical distortion in structure, and then measurement gained deformation is applied on finite element model node；Second is that passing through spy
Buckling mode analysis is levied, low order buckling mode is introduced into nonlinear buckling analysis model by specified deformation amplitude.First
Kind method really reflects the practical distortion of structure, but process is excessively complicated, and easily malfunctions, second method process
It is relatively easy, but introduced deformation is not the true strain of structure, therefore can not accurately calculate the true buckling of structure
And ultimate strength.The introducing of residual stress is typically all that Utopian residual stress distribution form or empirical formula is direct
It is applied in calculative finite element model, although convenience of calculation, can not really reflect the residual of physical presence in structure
Residue stress.
Summary of the invention
In view of the defects and deficiencies of the prior art, the present invention intends to provide one kind can really consider welding deformation
With the ultimate strength numerical computation method of residual stress.
In order to solve the problems existing in background technology, the ultimate strength of meter of the invention and welding deformation and residual stress
Numerical computation method, its analogy method are as follows:
A, it is tested by material mechanical performance, determines thin-wall construction material parameter；
B, the welding condition of query structure writes welding subprogram；
C, temperature-displacement coupling analysis model is established, welding analog analysis is carried out, obtains welding deformation and residual stress；
D, definition restarts analysis on the basis of temperature-displacement coupling analysis model, carries out the analysis of feature buckling mode；
E, establish static arc-length methods analysis model, will comprising including welding deformation the deformation of low order feature buckling mode and
Welding residual stress is introduced into arc-length methods analysis model, carries out structural limits Strength co-mputation.
The invention has the following beneficial effects: it not only can efficiently but also accurately calculate the thin of meter and welding deformation and residual stress
The ultimate strength of wall construction is all of great importance to scientific research and engineering design.
Detailed description of the invention
Fig. 1 is that the present invention is based on the Finite Element Method flow charts of ABAQUS.
Specific embodiment
With reference to the accompanying drawings and examples, the present invention is further illustrated.
In order to make the objectives, technical solutions, and advantages of the present invention clearer, below in conjunction with attached drawing and specific implementation
Mode, the present invention will be described in further detail.It should be appreciated that the specific embodiments described herein are only to explain this
Invention, is not intended to limit the present invention.
For example, to internal diameter 1.0m, ultimate strength of the titanium alloy TC 4 pressurized spherical shell of wall thickness 30mm under uniform external effect
Calculating analysis is carried out, which is that two identical hemisphere are welded, and is counted using the present invention and welding deformation and remnants
The ultimate strength numerical value of stress calculates, its calculating step are as follows:
A, it is tested by material mechanical performance, determines thin-wall construction material parameter.
1. elasticity modulus: 110GPa；Poisson's ratio: 0.3；Yield strength: 860MPa；
B, the welding condition of query structure writes welding subprogram.
1. welding process uses moving heat source method, the FORTRAN subprogram that establishment is calculated for ABAQUS, the mobile speed of heat source
Degree is 10mm/s；
C, temperature-displacement coupling analysis model is established, welding analog analysis is carried out, obtains welding deformation and residual stress.
1. establishing the three-dimensional entity model of pressurized spherical shell, material and section attribute are assigned, boundary condition is arranged in grid division
And predefined temperature field；
2. two temperature-displacement coupling analysis step is arranged: STEP-1 is welding process, and STEP-2 is cooling procedure, and is set
That sets the two analysis steps restarts request；
3. the interaction of spherical shell and environment is arranged are as follows: surface heat exchange and surface emissivity；Calorifics load is set are as follows: body heat
Flux；
4. FORTRAN subprogram file is added in JOB module and then submits operation；
D, definition restarts analysis on the basis of temperature-displacement coupling analysis model, carries out the analysis of feature buckling mode.
1. duplication temperature-displacement coupling model and edit model restarts attribute；
2. increasing single-sweep polarograpy buckling analysis after STEP-2 walks STEP-3, preceding 5 rank feature buckling mode is sought；
3. editor's keyword enables to export destination file .fil file；
4. submitting operation；
E, establish static arc-length methods analysis model, will comprising including welding deformation the deformation of low order feature buckling mode and
Welding residual stress is introduced into arc-length methods analysis model, carries out structural limits Strength co-mputation.
1. establishing the three-dimensional entity model of pressurized spherical shell, material and section attribute are assigned, perimeter strip is arranged in grid division
Part；
2. creating predefined stress field .odb file obtained in step C is added in predefined field, it is residual to introduce welding
Residue stress；
It will include the feature buckling mode change of welding deformation obtained in step D 3. editing keyword * Imperfection
Shape is introduced into Ultimate Strength Calculation；
4. FORTRAN subprogram file is added in JOB module and then submits operation；
5. being post-processed, ultimate strength value is obtained.
The invention has the following beneficial effects: it not only can efficiently but also accurately calculate the thin of meter and welding deformation and residual stress
The ultimate strength of wall construction is all of great importance to scientific research and engineering design.
The above is only used to illustrate the technical scheme of the present invention and not to limit it, and those of ordinary skill in the art are to this hair
The other modifications or equivalent replacement that bright technical solution is made, as long as it does not depart from the spirit and scope of the technical scheme of the present invention,
It is intended to be within the scope of the claims of the invention.
Claims (1)
1. the ultimate strength numerical computation method of a kind of meter and welding deformation and residual stress, it is characterised in that according to the following steps
It carries out:
A, it is tested by material mechanical performance, determines thin-wall construction material parameter；
1. elasticity modulus: 110GPa；Poisson's ratio: 0.3；Yield strength: 860MPa；
B, the welding condition of query structure writes welding subprogram；
1. welding process uses moving heat source method, the FORTRAN subprogram that establishment is calculated for ABAQUS, heat source movement speed is
10mm/s；
C, temperature-displacement coupling analysis model is established, welding analog analysis is carried out, obtains welding deformation and residual stress；
1. establishing the three-dimensional entity model of pressurized spherical shell, material and section attribute are assigned, boundary condition and pre- is arranged in grid division
Define temperature field；
2. two temperature-displacement coupling analysis step is arranged: STEP-1 is welding process, and STEP-2 is cooling procedure, and this is arranged
Two analysis steps restart request；
3. the interaction of spherical shell and environment is arranged are as follows: surface heat exchange and surface emissivity；Calorifics load is set are as follows: body heat is logical
Amount；
4. FORTRAN subprogram file is added in JOB module and then submits operation；
D, definition restarts analysis on the basis of temperature-displacement coupling analysis model, carries out the analysis of feature buckling mode；
1. duplication temperature-displacement coupling model and edit model restarts attribute；
2. increasing single-sweep polarograpy buckling analysis after STEP-2 walks STEP-3, preceding 5 rank feature buckling mode is sought；
3. editor's keyword enables to export destination file .fil file；
4. submitting operation；
E, static arc-length methods analysis model is established, will deform and weld comprising the low order feature buckling mode including welding deformation
Residual stress is introduced into arc-length methods analysis model, carries out structural limits Strength co-mputation；
1. establishing the three-dimensional entity model of pressurized spherical shell, material and section attribute are assigned, boundary condition is arranged in grid division；
2. creating predefined stress field .odb file obtained in step C is added in predefined field, welding remnants is introduced and answers
Power；
3. editing keyword * Imperfection, the feature buckling mode deformation obtained in step D comprising welding deformation is drawn
Enter into Ultimate Strength Calculation；
4. FORTRAN subprogram file is added in JOB module and then submits operation；
5. being post-processed, ultimate strength value is obtained.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610741183.1A CN106326574B (en) | 2016-08-27 | 2016-08-27 | It is a kind of meter and welding deformation and residual stress ultimate strength numerical computation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610741183.1A CN106326574B (en) | 2016-08-27 | 2016-08-27 | It is a kind of meter and welding deformation and residual stress ultimate strength numerical computation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106326574A CN106326574A (en) | 2017-01-11 |
| CN106326574B true CN106326574B (en) | 2019-10-15 |
Family
ID=57791134
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610741183.1A Expired - Fee Related CN106326574B (en) | 2016-08-27 | 2016-08-27 | It is a kind of meter and welding deformation and residual stress ultimate strength numerical computation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106326574B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109949879B (en) * | 2019-03-18 | 2023-05-30 | 西北工业大学 | Self-reinforcing method for pressure-resistant shell of underwater vehicle |
| CN110737997B (en) * | 2019-08-05 | 2024-01-26 | 中国船舶科学研究中心（中国船舶重工集团公司第七0二研究所） | Analysis method for ultimate strength of ultra-large floating body structure under complex load |
| CN113849923B (en) * | 2021-08-19 | 2022-05-24 | 北京市机械施工集团有限公司 | ABAQUS-based numerical simulation method for overall welding deformation of steel structure |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103246772A (en) * | 2013-05-11 | 2013-08-14 | 天津大学 | ABAQUS-based finite element simulation method of correcting welding deformation through ultrasonic shot-peening |
| CN103714207A (en) * | 2013-12-17 | 2014-04-09 | 天津大学 | Subsea pipeline three-dimensional dynamic numerical simulation method for introducing initial defects based on modes |
| CN105224750A (en) * | 2015-10-10 | 2016-01-06 | 北京工业大学 | A kind of new spatial based on response surface can open up single reed structure optimization method in hinge |
-
2016
- 2016-08-27 CN CN201610741183.1A patent/CN106326574B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103246772A (en) * | 2013-05-11 | 2013-08-14 | 天津大学 | ABAQUS-based finite element simulation method of correcting welding deformation through ultrasonic shot-peening |
| CN103714207A (en) * | 2013-12-17 | 2014-04-09 | 天津大学 | Subsea pipeline three-dimensional dynamic numerical simulation method for introducing initial defects based on modes |
| CN105224750A (en) * | 2015-10-10 | 2016-01-06 | 北京工业大学 | A kind of new spatial based on response surface can open up single reed structure optimization method in hinge |
Non-Patent Citations (4)
| Title |
|---|
| 《单层球面网壳结构的稳定性研究》;刘旭凡等;《四川建筑》;20061231;第26卷(第6期);第116-117页 * |
| 《基于仿真的铝合金型材焊接残余应力研究》;李灿;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20140615(第6期);摘要，正文第14-15、25-42页 * |
| 《结构极限强度分析的三种有限元解法研究》;彭大炜等;《中国海洋平台》;20100430;第25卷(第2期);第1-5页 * |
| 《船体梁极限强度的非线性有限元分析》;白宝强;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20111215(第S1期);摘要，正文第42-88页第4-5章 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106326574A (en) | 2017-01-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106326574B (en) | It is a kind of meter and welding deformation and residual stress ultimate strength numerical computation method | |
| CN103049623B (en) | A kind of method for building up of laser bonding heat source model | |
| CN103246772A (en) | ABAQUS-based finite element simulation method of correcting welding deformation through ultrasonic shot-peening | |
| Xiao et al. | Forming limit in thermal cruciform biaxial tensile testing of titanium alloy | |
| CN102592019B (en) | Simulation method for welding residual stress of quartz flexible accelerometer | |
| CN105138772B (en) | The Finite Element Method of non-uniform member electron beam welding | |
| CN104372285B (en) | A kind of manufacturing method of high-temperature speckle based on plasma spraying | |
| CN106021644B (en) | The method for determining mixed dimensional modelling interface constraint equation coefficient | |
| CN105548005A (en) | Method for evaluating fatigue crack propagation characteristics of CFRP reinforced steel structure based on SIF values | |
| Luo et al. | Effects of process parameters on deformation and temperature uniformity of forged Ti-6Al-4V turbine blade | |
| CN107180131B (en) | Method for determining deformation curvature radius of multipoint laser shock peening thin-walled part | |
| CN103994747B (en) | The measuring method of cloth without basis point of beam-shaped structure stretch bending composite deformation field | |
| Guo et al. | A method to realize high-precision and large laser thermal bending angle | |
| Chen et al. | A criterion for void closure in the porous model during the forging of steel ingot and its application | |
| Cheng et al. | Finite element modeling and experimental investigation on manufacturing TA18 alloy pipes via hot free bending forming technology: Forming characteristics and process optimization | |
| Suárez et al. | Fracture of eutectoid steel bars under tensile loading: Experimental results and numerical simulation | |
| Lu et al. | Numerical simulation of welding distortion using shrinkage force approach and application | |
| CN108563800A (en) | A kind of method for numerical simulation improving residual stress distribution process based on ultrasonic impact | |
| Zhao et al. | Springback theory of plane bending and the progress of study on its engineering application | |
| CN205620080U (en) | Vacuum sheet glass vacuum on -line measuring system | |
| CN105760574B (en) | A kind of method of thermo parameters method in simulation solid ferment process | |
| CN105912746A (en) | Mechanical model building method based on mould temperature rise rate under condition of continuous plastic processing and production and application thereof | |
| Cai et al. | Line Gauss heat source model: an efficient approach for numerical welding simulation | |
| Wang et al. | An effective thermal-mechanical coupling method for simulating friction stir-assisted incremental aluminum alloy sheet forming | |
| CN205898548U (en) | Concrete temperature stress testing machine with temperature deformation is from compensation function |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191015 Termination date: 20210827 |