CN106326574B - 一种计及焊接变形和残余应力的极限强度数值计算方法 - Google Patents
一种计及焊接变形和残余应力的极限强度数值计算方法 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
本发明公开了一种计及焊接变形和残余应力的极限强度数值计算方法,它涉及薄壁结构极限强度计算的数值模拟方法技术领域;它的模拟方法为:通过材料力学性能试验,确定薄壁结构材料参数;查询结构的焊接工艺参数,编写焊接子程序;建立温度‑位移耦合分析模型,进行焊接模拟分析,得到焊接变形和残余应力;在温度‑位移耦合分析模型基础上定义重启动分析,进行特征屈曲模态分析;建立静态弧长法分析模型,将包含焊接变形在内的低阶特征屈曲模态变形以及焊接残余应力引入到弧长法分析模型中进行结构极限强度计算。它能够既高效又准确的计算出计及焊接变形和残余应力的薄壁结构的极限强度,对科学研究和工程设计都有重大意义。
Description
技术领域
本发明涉及一种计及焊接变形和残余应力的极限强度数值计算方法,具体的是一种基于ABAQUS的将焊接变形和残余应力引入薄壁结构极限强度计算的数值模拟方法。
背景技术
薄壁结构在某些载荷下会发生失稳现象,如受均匀外压的耐压球壳,受面内轴向压缩载荷作用的加筋板结构等。在用数值模拟方法分析薄壁结构的稳定性时,需要考虑结构中的初始变形和残余应力,目前,结构的初始变形主要通过两种方法引入:一是利用测量设备测出结构中的实际变形,然后将测量所得变形施加到有限元模型节点上;二是通过特征屈曲模态分析,将低阶屈曲模态按指定的变形幅值引入到非线性屈曲分析模型中。第一种方法真实的反映了结构的实际变形,但是过程过于复杂,并且极易出错,第二种方法过程相对简单,但是所引入的变形并不是结构的真实变形,因此无法准确计算结构真实的屈曲和极限强度。残余应力的引入一般都是将理想化的残余应力分布形式或经验计算公式直接施加到需要计算的有限元模型中,虽然计算方便,但却无法真实反映结构中实际存在的残余应力。
发明内容
本发明的目的在于针对现有技术的缺陷和不足,提供一种能够真实考虑焊接变形和残余应力的极限强度数值计算方法。
为了解决背景技术所存在的问题,本发明的计及焊接变形和残余应力的极限强度数值计算方法,它的模拟方法为:
A、通过材料力学性能试验,确定薄壁结构材料参数;
B、查询结构的焊接工艺参数,编写焊接子程序;
C、建立温度-位移耦合分析模型,进行焊接模拟分析,得到焊接变形和残余应力;
D、在温度-位移耦合分析模型基础上定义重启动分析,进行特征屈曲模态分析;
E、建立静态弧长法分析模型,将包含焊接变形在内的低阶特征屈曲模态变形以及焊接残余应力引入到弧长法分析模型中,进行结构极限强度计算。。
本发明有益效果为:它能够既高效又准确的计算出计及焊接变形和残余应力的薄壁结构的极限强度,对科学研究和工程设计都有重大意义。
附图说明
图1为本发明基于ABAQUS的有限元模拟方法流程图。
具体实施方式
下面结合附图和实施例,对本发明作进一步的说明。
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及具体实施方式,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施方式仅用以解释本发明,并不用于限定本发明。
例如,对内径1.0m,壁厚30mm的钛合金TC4耐压球壳在均匀外压作用下的极限强度进行计算分析,该球壳为两个相同的半球焊接而成,利用本发明进行计及焊接变形和残余应力的极限强度数值计算,它的计算步骤为:
A、通过材料力学性能试验,确定薄壁结构材料参数。
①弹性模量:110GPa;泊松比:0.3;屈服强度:860MPa;
B、查询结构的焊接工艺参数,编写焊接子程序。
①焊接过程采用移动热源法,编制用于ABAQUS计算的FORTRAN子程序,热源移动速度为10mm/s;
C、建立温度-位移耦合分析模型,进行焊接模拟分析,得到焊接变形和残余应力。
①建立耐压球壳的三维实体模型,赋予材料和截面属性,划分网格,设置边界条件及预定义温度场;
②设置两个温度-位移耦合分析步:STEP-1为焊接过程,STEP-2为冷却过程,并设置这两个分析步的重启动请求;
③设置球壳与环境的相互作用为:表面热交换和表面辐射;设置热学载荷为:体热通量;
④将FORTRAN子程序文件添加到JOB模块中然后提交运算;
D、在温度-位移耦合分析模型基础上定义重启动分析,进行特征屈曲模态分析。
①复制温度-位移耦合模型并编辑模型的重启动属性;
②在STEP-2之后增加线性摄动屈曲分析步STEP-3,求取前5阶特征屈曲模态;
③编辑关键字使之能够输出结果文件.fil文件;
④提交运算;
E、建立静态弧长法分析模型,将包含焊接变形在内的低阶特征屈曲模态变形以及焊接残余应力引入到弧长法分析模型中,进行结构极限强度计算。
①建立耐压球壳的三维实体模型,赋予材料和截面属性,划分网格,设置边界条件;
②创建预定义应力场,将步骤C中得到的.odb文件添加到预定义场中,引入焊接残余应力;
③编辑关键字*Imperfection,将步骤D中得到的包含焊接变形的特征屈曲模态变形引入到极限强度计算中;
④将FORTRAN子程序文件添加到JOB模块中然后提交运算;
⑤进行后处理,得到极限强度值。
本发明有益效果为:它能够既高效又准确的计算出计及焊接变形和残余应力的薄壁结构的极限强度,对科学研究和工程设计都有重大意义。
以上所述,仅用以说明本发明的技术方案而非限制,本领域普通技术人员对本发明的技术方案所做的其它修改或者等同替换,只要不脱离本发明技术方案的精神和范围,均应涵盖在本发明的权利要求范围当中。
Claims (1)
1.一种计及焊接变形和残余应力的极限强度数值计算方法,其特征在于按照以下步骤进行:
A、通过材料力学性能试验,确定薄壁结构材料参数;
①弹性模量:110GPa;泊松比:0.3;屈服强度:860MPa;
B、查询结构的焊接工艺参数,编写焊接子程序;
①焊接过程采用移动热源法,编制用于ABAQUS计算的FORTRAN子程序,热源移动速度为10mm/s;
C、建立温度-位移耦合分析模型,进行焊接模拟分析,得到焊接变形和残余应力;
①建立耐压球壳的三维实体模型,赋予材料和截面属性,划分网格,设置边界条件及预定义温度场;
②设置两个温度-位移耦合分析步:STEP-1为焊接过程,STEP-2为冷却过程,并设置这两个分析步的重启动请求;
③设置球壳与环境的相互作用为:表面热交换和表面辐射;设置热学载荷为:体热通量;
④将FORTRAN子程序文件添加到JOB模块中然后提交运算;
D、在温度-位移耦合分析模型基础上定义重启动分析,进行特征屈曲模态分析;
①复制温度-位移耦合模型并编辑模型的重启动属性;
②在STEP-2之后增加线性摄动屈曲分析步STEP-3,求取前5阶特征屈曲模态;
③编辑关键字使之能够输出结果文件.fil文件;
④提交运算;
E、建立静态弧长法分析模型,将包含焊接变形在内的低阶特征屈曲模态变形以及焊接残余应力引入到弧长法分析模型中,进行结构极限强度计算;
①建立耐压球壳的三维实体模型,赋予材料和截面属性,划分网格,设置边界条件;
②创建预定义应力场,将步骤C中得到的.odb文件添加到预定义场中,引入焊接残余应力;
③编辑关键字*Imperfection,将步骤D中得到的包含焊接变形的特征屈曲模态变形引入到极限强度计算中;
④将FORTRAN子程序文件添加到JOB模块中然后提交运算;
⑤进行后处理,得到极限强度值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610741183.1A CN106326574B (zh) | 2016-08-27 | 2016-08-27 | 一种计及焊接变形和残余应力的极限强度数值计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610741183.1A CN106326574B (zh) | 2016-08-27 | 2016-08-27 | 一种计及焊接变形和残余应力的极限强度数值计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106326574A CN106326574A (zh) | 2017-01-11 |
CN106326574B true CN106326574B (zh) | 2019-10-15 |
Family
ID=57791134
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610741183.1A Expired - Fee Related CN106326574B (zh) | 2016-08-27 | 2016-08-27 | 一种计及焊接变形和残余应力的极限强度数值计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106326574B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109949879B (zh) * | 2019-03-18 | 2023-05-30 | 西北工业大学 | 一种水下航行器耐压壳体自增强方法 |
CN110737997B (zh) * | 2019-08-05 | 2024-01-26 | 中国船舶科学研究中心(中国船舶重工集团公司第七0二研究所) | 复杂载荷作用下超大型浮体结构极限强度的分析方法 |
CN113849923B (zh) * | 2021-08-19 | 2022-05-24 | 北京市机械施工集团有限公司 | 一种基于abaqus的钢结构整体焊接变形的数值模拟方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103246772A (zh) * | 2013-05-11 | 2013-08-14 | 天津大学 | 一种基于abaqus的超声喷丸矫正焊接变形的有限元模拟方法 |
CN103714207A (zh) * | 2013-12-17 | 2014-04-09 | 天津大学 | 基于模态引入初始缺陷的海底管线三维动力数值模拟方法 |
CN105224750A (zh) * | 2015-10-10 | 2016-01-06 | 北京工业大学 | 一种基于响应面的新型空间可展铰链中单簧片结构优化方法 |
-
2016
- 2016-08-27 CN CN201610741183.1A patent/CN106326574B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103246772A (zh) * | 2013-05-11 | 2013-08-14 | 天津大学 | 一种基于abaqus的超声喷丸矫正焊接变形的有限元模拟方法 |
CN103714207A (zh) * | 2013-12-17 | 2014-04-09 | 天津大学 | 基于模态引入初始缺陷的海底管线三维动力数值模拟方法 |
CN105224750A (zh) * | 2015-10-10 | 2016-01-06 | 北京工业大学 | 一种基于响应面的新型空间可展铰链中单簧片结构优化方法 |
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 (zh) | 2017-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104298818B (zh) | 一种端铣加工表面误差预测及仿真方法 | |
CN106326574B (zh) | 一种计及焊接变形和残余应力的极限强度数值计算方法 | |
CN103049623B (zh) | 一种激光焊接热源模型的建立方法 | |
CN103246772A (zh) | 一种基于abaqus的超声喷丸矫正焊接变形的有限元模拟方法 | |
Xiao et al. | Forming limit in thermal cruciform biaxial tensile testing of titanium alloy | |
CN105138772B (zh) | 变截面构件电子束焊的有限元模拟方法 | |
CN104372285B (zh) | 一种基于等离子喷涂的高温散斑制作方法 | |
CN109977617B (zh) | 一种基于apdl语言的钢板火焰切割热-力耦合仿真分析方法 | |
Luo et al. | Effects of process parameters on deformation and temperature uniformity of forged Ti-6Al-4V turbine blade | |
CN107180131B (zh) | 用于确定多点激光冲击强化薄壁件变形曲率半径的方法 | |
Guo et al. | A method to realize high-precision and large laser thermal bending angle | |
Song et al. | Numerical and experimental study on the properties of laser cladding of 6061Al alloy | |
CN108629110A (zh) | 快速获取修正金属粉末压制成形dpc本构模型参数的方法 | |
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 | |
CN108563800A (zh) | 一种基于超声冲击改善残余应力分布过程的数值模拟方法 | |
Zhao et al. | Springback theory of plane bending and the progress of study on its engineering application | |
CN205620080U (zh) | 一种真空平板玻璃真空度在线检测*** | |
Lorenz et al. | Recent advances and new developments in hot forming simulation with LS-DYNA | |
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 | |
CN103745068A (zh) | 一种套筒铆接结构的有限元建模方法 | |
CN104765928A (zh) | 一种塑性成形摩擦特性测算方法 | |
Novella et al. | Numerical modelling of AlSi7 tubular components flowformed at elevated temperature | |
CN107052212A (zh) | 多腔类构件多向加载成形加载路径的确定方法 |
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 |