CN103008827A - 一种应用于水环境中焊接物理场耦合的系统分析方法 - Google Patents

一种应用于水环境中焊接物理场耦合的系统分析方法 Download PDF

Info

Publication number
CN103008827A
CN103008827A CN2012105246051A CN201210524605A CN103008827A CN 103008827 A CN103008827 A CN 103008827A CN 2012105246051 A CN2012105246051 A CN 2012105246051A CN 201210524605 A CN201210524605 A CN 201210524605A CN 103008827 A CN103008827 A CN 103008827A
Authority
CN
China
Prior art keywords
welding
curve
water
boiling
condition
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.)
Granted
Application number
CN2012105246051A
Other languages
English (en)
Other versions
CN103008827B (zh
Inventor
李航
王立君
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN201210524605.1A priority Critical patent/CN103008827B/zh
Publication of CN103008827A publication Critical patent/CN103008827A/zh
Application granted granted Critical
Publication of CN103008827B publication Critical patent/CN103008827B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Arc Welding In General (AREA)

Abstract

本发明公开了一种应用于水环境中焊接物理场耦合的系统分析方法，其步骤是：参考水的大容器饱和沸腾曲线，根据公式h=q/Δt计算出上述曲线中沸腾状态变化点的对流散热系数h，并拟合得到对流散热系数随焊件表面温度tw变化的曲线，将此曲线对应的分段函数作为水下湿法焊接有限元模型的边界散热条件；在确定水下湿法焊接条件的基础上，利用Visual Mesh建立焊板几何模型和进行有限元网格划分，结合上边得到的边界散热条件，利用SYSWELD软件计算了水下湿法焊接温度场的分布情况。本发明可以合理预测水下湿法焊接温度场的分布情况，在实际生产中有重要的科学和工程意义，可以减少进行水下焊接工艺试验的工作量，降低工艺试验成本。

Description

一种应用于水环境中焊接物理场耦合的系统分析方法
技术领域
本发明属于材料加工工程学科中的焊接领域，尤其涉及一种应用于水环境中焊接物理场耦合的系统分析方法。
背景技术
海洋蕴含着丰富的自然资源。海洋的开发需要大量海洋工程的建设，而这必定离不开水下焊接技术的应用。水下焊接技术从电弧环境上主要区分为：湿法、干法和局部干法。湿法焊接是指直接在水环境中进行焊接，水与电弧之间没有任何隔离设施，电弧主要靠焊条药皮或药芯燃烧生成的气罩进行保护。干法焊接是指用大型隔离舱将焊接件与水环境隔开，形成类似于陆上焊接的工况环境。局部干法是指将局部焊接区域的水排开，保证施焊位置局部干燥。
焊接过程是个复杂的传热传质问题。在计算焊接温度场时，电弧、焊件与周围环境之间热对流系数的取值决定着焊接温度场分析的准确与否。以往计算大气条件下的焊接温度场时往往设定焊件与空气环境之间的对流换热系数hc为某个常数，取值范围为1~100W/(m2·K)。与此对应，在为数不多的几篇有关水下湿法焊接温度场模拟的研究文献中，把焊件与水环境之间的对流换热系数hc也常设定为某个常数，取值范围为200~1000W/(m2·K)。但对于水下湿法焊接而言，水环境的存在使得该焊接过程的热传递变得更为复杂。随着焊件表面温度的上升，与其接触的水介质会发生沸腾和汽化，水介质的相变必然会导致焊件与水环境之间对流传热系数的变化，此时再将其取为常数显然不能恰当描述水下湿法焊接的物理场耦合情况。
发明内容
针对上述现有技术，本发明提供一种应用于水环境中焊接物理场耦合的系统分析方法。本发明结合水的饱和沸腾传热原理，得到水下湿法焊接过程中，焊件与水环境之间的对流传热系数随焊件温度变化的曲线，并以此为基础给出应用于水下湿法焊接的温度场分析流程。本发明的目的在于改进水下湿法焊接物理场的分析方法，使分析结果更符合实际工况。
为了解决上述技术问题，本发明一种应用于水环境中焊接物理场耦合的系统分析方法，步骤如下：
步骤一：参考水的大容器饱和沸腾曲线，根据焊件表面过热度Δt增加过程中水沸腾传热的不同机理，将该曲线分为自然对流、核状沸腾、过度沸腾和稳定膜态沸腾四个区；其中：Δt=tw-ts，tw表示焊件表面温度，ts表示水环境温度，Δt表示焊件表面过热度；
步骤二：根据公式h=q/Δt分别计算出上述水的大容器饱和沸腾曲线中沸腾状态变化特征点A、B、C、D、E、F的对流散热系数h，并拟合得到对流散热系数随焊件表面温度tw变化的曲线，该曲线的分段函数如下：
h = 205 , t w ≤ 100 ; h = 250 + 353 ( t - 100 ) , t w ≤ 130 ; h = 10840 - 109 ( t - 130 ) , t w ≤ 220 ; h = 1030 , t w ≤ 2000 ; - - - ( 1 )
公式（1）中，tw的单位为℃，h的单位为W/(m2·K)；
步骤三：确定水下湿法焊接条件，至少包括：环境温度、焊接工艺参数、焊板几何尺寸，焊板材质以及坡口形式；
步骤四：依据上述焊板几何尺寸用网格划分软件对焊板建立几何建模和进行有限元网格划分；
步骤五：根据实际焊接条件，初步设定水下湿法焊接模型的热输入条件、边界散热条件和边界约束条件，其中，将公式（1）作为水下湿法堆焊有限元模型的边界散热条件；
步骤六：利用SYSWELD软件计算水下湿法焊接的温度场，提取焊接热循环参数。
与现有技术相比，本发明的有益效果是：
本发明结合水的饱和沸腾传热原理，得到水下湿法焊接过程中固液耦合界面对流传热系数随焊件温度变化的曲线，并以此为基础提出了分析水下湿法焊接温度场的改进方法。本发明提供的技术思想可以合理预测水下湿法焊接温度场的分布情况，在实际生产中有重要的科学和工程意义，可以减少进行水下焊接工艺试验的工作量，降低工艺试验成本。
附图说明
图1是水的大容器饱和沸腾曲线；
图2是对流散热系数随焊件表面温度变化的曲线；
图3水下焊接的网格模型；
图4计算所得的焊缝熔池形状同实际熔池形状的对比；
图5计算所得的焊缝熔合线处的温度循环曲线；
图6计算所得的焊缝热影响区处的温度循环曲线。
具体实施方式
下面对本发明的实施例作详细说明：本实施例以本发明技术方案为前提进行实施，但本发明的保护范围不限于下述的实施例。
本发明一种应用于水环境中焊接物理场耦合的系统分析方法，步骤如下：
步骤一：参考水的大容器饱和沸腾曲线，如图1所示，根据加热面过热度Δt增加过程中水沸腾传热的不同机理，将该曲线分为自然对流、核状沸腾、过度沸腾和稳定膜态沸腾四个区；其中：Δt＝tw-ts，tw表示焊件表面温度，ts表示水环境温度，Δt表示加热面过热度；
步骤二：根据公式h=q/Δt分别计算出上述水的大容器饱和沸腾曲线中沸腾状态变化特征点A、B、C、D、E、F的对流散热系数h，（其中，AB段属于孤立汽泡区，BC属于汽气块区，CD段属于大汽块区，DE属于过渡沸腾区，EF段属于稳定膜态沸腾区）并拟合得到对流散热系数随焊件表面温度tw变化的曲线，如图2所示，该曲线的分段函数如下：
h = 205 , t w ≤ 100 ; h = 250 + 353 ( t - 100 ) , t w ≤ 130 ; h = 10840 - 109 ( t - 130 ) , t w ≤ 220 ; h = 1030 , t w ≤ 2000 ; - - - ( 1 )
公式（1）中，tw的单位为℃，h的单位为W/(m2·K)；
步骤三：确定水下湿法焊接条件，至少包括：环境温度、焊接工艺参数、焊板几何尺寸，焊板材质以及坡口形式；
本实施例采用水下湿法堆焊，焊接条件为：环境温度设为20℃，平板试样，板厚7.5mm，长度100mm宽度60mm；板材为Q345；焊条采用国产的TS202水下专用焊条，堆焊规范如表1所示：
表1水下湿法堆焊焊接规范
焊接电流 焊接电压 焊接速度 焊条直径 焊条倾角
160~170A 35~38V 3.1mm/s 4.0mm 90°±2°
步骤四：依据上述焊板几何尺寸用网格划分软件，诸如Visual Mesh、Hyper Mesh，对焊板建立几何建模和进行有限元网格划分，单元类型为三维六面体单元，在焊缝及其周围5mm处进行网格细化，划分好的网格模型如图3所示。
步骤五：根据实际焊接条件，初步设定水下湿法堆焊模型的热输入条件、边界散热条件和边界约束条件，具体过程如下：
（1）热输入条件：用有限元分析软件SYSWELD对热源模型的热源参数进行校核，热输入系数取0.7，热源形式为双椭球与平面高斯复合热源，能量比为3:1；
（2）将公式（1）作为水下湿法堆焊有限元模型的边界散热条件，即：依据公式（1）编写子程序如下：
Figure BDA00002532172900032
Figure BDA00002532172900041
Figure BDA00002532172900051
将其导入有限元分析软件SWSWELD中的函数库，从而将其设为水下堆焊有限元模型的边界散热条件；
（3）设定水下湿法堆焊模型的边界约束条件：在焊板模型的端面选取不在一条直线上三点，分别记为点1、点2、点3，其中，对点1施加X，Y，Z方向约束，对点2施加X，Y方向约束；对点3施加X方向约束。
步骤六：利用SYSWELD软件计算水下湿法堆焊模型的温度场，提取其熔池截面形状并同实际的熔池形状作对比，如图4所示,通过观察可以发现计算所得的熔池形状与实际熔池形状吻合情况良好，这说明该方法的计算准确度较高，计算结果能够很好地反映真实的焊接温度场。最后提取熔合线及热影响区处的焊接热循环参数分别如图5和6所示，至此完成该水下堆焊温度场的分析模拟。
尽管上面结合图对本发明进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨的情况下，还可以作出很多变形，这些均属于本发明的保护之内。

Claims (1)

1.一种应用于水环境中焊接物理场耦合的系统分析方法，其特征在于，步骤如下：
步骤一：参考水的大容器饱和沸腾曲线，根据焊件表面过热度Δt增加过程中水沸腾传热的不同机理，将该曲线分为自然对流、核状沸腾、过度沸腾和稳定膜态沸腾四个区；其中：Δt＝tw-ts，tw表示焊件表面温度，ts表示水环境温度，Δt表示焊件表面过热度；
步骤二：根据公式h=q/Δt分别计算出上述水的大容器饱和沸腾曲线中沸腾状态变化特征点A、B、C、D、E、F的对流散热系数h，并拟合得到对流散热系数随焊件表面温度tw变化的曲线，该曲线的分段函数如下：
h = 205 , t w ≤ 100 ; h = 250 + 353 ( t - 100 ) , t w ≤ 130 ; h = 10840 - 109 ( t - 130 ) , t w ≤ 220 ; h = 1030 , t w ≤ 2000 ; - - - ( 1 )
公式（1）中，tw的单位为℃，h的单位为W/(m2·K)；
步骤三：确定水下湿法焊接条件，至少包括：环境温度、焊接工艺参数、焊板几何尺寸，焊板材质以及坡口形式；
步骤四：依据上述焊板几何尺寸用网格划分软件对焊板建立几何建模和进行有限元网格划分；
步骤五：根据实际焊接条件，初步设定水下湿法焊接模型的热输入条件、边界散热条件和边界约束条件，其中，将公式（1）作为水下湿法堆焊有限元模型的边界散热条件；
步骤六：利用SYSWELD软件计算水下湿法焊接的温度场，提取焊接热循环参数。
CN201210524605.1A 2012-12-05 2012-12-05 一种应用于水环境中焊接物理场耦合的系统分析方法 Expired - Fee Related CN103008827B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210524605.1A CN103008827B (zh) 2012-12-05 2012-12-05 一种应用于水环境中焊接物理场耦合的系统分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210524605.1A CN103008827B (zh) 2012-12-05 2012-12-05 一种应用于水环境中焊接物理场耦合的系统分析方法

Publications (2)

Publication Number Publication Date
CN103008827A true CN103008827A (zh) 2013-04-03
CN103008827B CN103008827B (zh) 2015-06-17

Family

ID=47958234

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210524605.1A Expired - Fee Related CN103008827B (zh) 2012-12-05 2012-12-05 一种应用于水环境中焊接物理场耦合的系统分析方法

Country Status (1)

Country Link
CN (1) CN103008827B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109657269A (zh) * 2018-11-05 2019-04-19 华中科技大学 一种不同水深压力下局部干法水下激光焊接数值模拟方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001030532A1 (en) * 1999-10-27 2001-05-03 Caterpillar Inc. Method and apparatus for providing a simulation of a welding process using integrated models
CN101559511A (zh) * 2009-05-22 2009-10-21 清华大学 一种以温度为控制变量的焊接数值模拟计算方法
KR20120050794A (ko) * 2010-11-11 2012-05-21 삼성중공업 주식회사 용접 변형 해석 방법
CN102693336A (zh) * 2012-05-09 2012-09-26 天津大学 一种预测大型管道焊接热循环参数的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001030532A1 (en) * 1999-10-27 2001-05-03 Caterpillar Inc. Method and apparatus for providing a simulation of a welding process using integrated models
CN101559511A (zh) * 2009-05-22 2009-10-21 清华大学 一种以温度为控制变量的焊接数值模拟计算方法
KR20120050794A (ko) * 2010-11-11 2012-05-21 삼성중공업 주식회사 용접 변형 해석 방법
CN102693336A (zh) * 2012-05-09 2012-09-26 天津大学 一种预测大型管道焊接热循环参数的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
J.RONDA等: "Thermomechanical simulation of underwater welding processes", 《ARCHIVE OF APPLIED MECHANICS》 *
贾剑平: "基于ANSYS的水下局部干法焊接过程温度场分析", 《电焊机》 *
赵博等: "水下湿法焊接温度场的数值模拟", 《第十六次全国焊接学术会议论文摘要集》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109657269A (zh) * 2018-11-05 2019-04-19 华中科技大学 一种不同水深压力下局部干法水下激光焊接数值模拟方法
CN109657269B (zh) * 2018-11-05 2021-06-08 华中科技大学 一种不同水深压力下局部干法水下激光焊接数值模拟方法

Also Published As

Publication number Publication date
CN103008827B (zh) 2015-06-17

Similar Documents

Publication Publication Date Title
Wang et al. Thermal stress analysis of eccentric tube receiver using concentrated solar radiation
Zhang et al. Numerical simulation investigation on thermal performance of heat pipe flat-plate solar collector
Cui et al. Heat transfer analysis of pile geothermal heat exchangers with spiral coils
Jiandong et al. Numerical simulation for structural parameters of flat-plate solar collector
Liu et al. Flow field optimization for the solar parabolic trough receivers in direct steam generation systems by the variational principle
Zhang et al. Study on heat transfer of pile foundation ground heat exchanger with three-dimensional groundwater seepage
Liang et al. Performance analysis of a new-design filled-type solar collector with double U-tubes
Das et al. Role of multiple solar heaters along the walls for the thermal management during natural convection in square and triangular cavities
CN103279630B (zh) 激光深熔点焊匙孔动态成形数值模拟方法
CN103008827A (zh) 一种应用于水环境中焊接物理场耦合的系统分析方法
Abbas et al. Performance analysis of seasonal soil heat storage system based on numerical simulation and experimental investigation
Miao et al. Industrial experiment and numerical analysis of the effect of pulsing water flow in thermoelectric generators
Xie et al. Heat transfer enhancement for the coil zone of closed wet cooling towers through field synergy analysis
Li et al. Thermal analysis and stress analysis of the heat-exchange pipe based on ANSYS
Hussen et al. Design, optimization and CFD simulation of solar air heater with jet impingement on V-corrugated plate
CN104573206A (zh) 一种基于有限元热力耦合的型钢断面热形状尺寸设计方法
Pengfei et al. RETRACTED: Establishment and solution of the model for loop pipeline network with multiple heat sources
Guang-Pu et al. Review on the Numerical Simulation Methods of Cooling Pipe in Mass Concrete
Li et al. Factors affecting the performance of molten salt heat receiver
Zhao et al. Analysis of thermal performance of solar collector in solar water heating system
Xue et al. Heat transfer performance of wet porous solar collectors under periodic conditions
Zuo et al. The influence of heat collector structural parameters on the solar chimneys integrated with seawater desalination
Zhang et al. 3D Numerical Simulation of Heat Transfer for the Vertical U Type Berried Pipe of the Ground Source Heat Pump
Liu et al. The research progress of earth to air heat exchanger: a review
Yi et al. Numerical Simulation of Straight Rib Solar Air Collector

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150617

Termination date: 20201205

CF01 Termination of patent right due to non-payment of annual fee