CN113579462A - 一种双侧同时搅拌摩擦焊工艺参数优化方法 - Google Patents
一种双侧同时搅拌摩擦焊工艺参数优化方法 Download PDFInfo
- Publication number
- CN113579462A CN113579462A CN202110986364.1A CN202110986364A CN113579462A CN 113579462 A CN113579462 A CN 113579462A CN 202110986364 A CN202110986364 A CN 202110986364A CN 113579462 A CN113579462 A CN 113579462A
- Authority
- CN
- China
- Prior art keywords
- workpiece
- friction stir
- stirring
- heat
- stir welding
- 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
Links
Images
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
- B23K20/00—Non-electric welding by applying impact or other pressure, with or without the application of heat, e.g. cladding or plating
- B23K20/12—Non-electric welding by applying impact or other pressure, with or without the application of heat, e.g. cladding or plating the heat being generated by friction; Friction welding
- B23K20/122—Non-electric welding by applying impact or other pressure, with or without the application of heat, e.g. cladding or plating the heat being generated by friction; Friction welding using a non-consumable tool, e.g. friction stir welding
-
- 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
- B23K20/00—Non-electric welding by applying impact or other pressure, with or without the application of heat, e.g. cladding or plating
- B23K20/26—Auxiliary equipment
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Pressure Welding/Diffusion-Bonding (AREA)
Abstract
本发明公开了一种双侧同时搅拌摩擦焊工艺参数优化方法,基于搅拌摩擦焊接过程产热机理、傅里叶传热定律和能量守恒定律,建立双侧同时搅拌摩擦焊的产热模型,确定待优化工艺参数,以焊接单位长度输入能量最小为目标函数,根据焊接工艺设置约束条件,建立优化模型;采用序列二次规划算法求解优化模型,得到最优工艺参数值。本发明基于焊接过程的产热机理,与SQP优化算法相耦合,得到了全局最优的输入能量和最佳的工艺参数组合,优化后的参数提高了焊接的效率和改善了焊接接头的质量。
Description
技术领域
本发明涉及一种搅拌摩擦焊工艺参数的优化方法,特别涉及一种双侧同时搅拌摩擦焊工艺参数优化方法。
背景技术
目前,双侧同时搅拌摩擦焊相对于常规搅拌摩擦焊具有产生对称的焊缝和热输入、提高生产效率、降低夹具上反力矩的优点。在搅拌摩擦焊焊接过程中,热量是决定最终焊接质量的关键因素,它通常提供两个主要功能:(i)使金属材料足够软,以便搅拌头搅拌形成固态焊缝,以及(ii)导致微观结构演变,从而影响接头性能。工艺参数是影响产热的主要因素,开展工艺参数优化方法的研究,有利于改善双侧同时搅拌摩擦焊焊接接头的质量和扩大该焊接方法的应用场景。到目前为止,还没有关于双侧同时搅拌摩擦焊工艺参数的报道,但人们对常规搅拌摩擦焊工艺参数优化的方法开展了大量的研究工作,以增加常规搅拌摩擦焊在工业中的应用。
常规搅拌摩擦焊工艺参数优化的方法主要有实验法和数值法。在采用实验法进行搅拌摩擦焊工艺参数优化时,需要做大量的实验和误差纠正工作,才能得到满意的结果,这种方法需要投入大量的精力、人力、财力且缺乏搅拌摩擦焊焊接机理理论指导,在一定程度上限制了该方法的应用。而数值方法由于其计算效率高、直观性强、可重复的优点逐渐成为预测搅拌摩擦焊参数的有利工具,但在工艺参数优化过程中需要进行大量的数值模拟且不能得到全局最优值,从而得不到满意的结果。因此,如何有效的避免数值法陷入局部最优并将其应用于双侧同时搅拌摩擦焊的工艺参数优化成为主要的问题。
发明内容
本发明的目的在于解决背景技术中针对搅拌摩擦焊工艺参数优化的实验法理论支撑不足、过程复杂和数值法的局部最优问题,提出了一种双侧同时搅拌摩擦焊工艺参数优化方法。
本发明为解决公知技术中存在的技术问题所采取的技术方案是:一种双侧同时搅拌摩擦焊工艺参数优化方法,基于搅拌摩擦焊接过程产热机理、傅里叶传热定律和能量守恒定律,建立双侧同时搅拌摩擦焊的产热模型,确定待优化工艺参数,以焊接单位长度输入能量最小为目标函数,根据焊接工艺设置约束条件,建立优化模型;采用序列二次规划算法求解优化模型,得到最优工艺参数值。
进一步地,待优化工艺参数包括:搅拌头的转速、搅拌头的进给焊接速度以及两侧搅拌针端面之间的距离。
进一步地,建立双侧同时搅拌摩擦焊的产热模型的方法包括:
由于两侧搅拌头关于工件厚度方向几何对称且两侧的搅拌头的焊接工艺参数相同,因此两侧搅拌头的产热模型相同,只建立其中一侧搅拌头的产热模型;
根据搅拌头与工件的接触面不同,将搅拌摩擦焊接过程产生的热功率分成以下四种:分别是搅拌头轴肩端面与工件摩擦产生的热功率Qsb、搅拌头轴肩侧面与工件摩擦产生的热功率Qss、搅拌针端面与工件摩擦产生的热功率Qpb、搅拌针侧面与工件摩擦产生的热功率Qps:
设:R2为搅拌头轴肩半径,R1为搅拌针端部半径,hp为搅拌针长度,h为轴肩下压量,hg为两侧搅拌针端面之间的距离,H为工件厚度,ω为搅拌头的转速,σyield为工件材料的屈服应力,基于搅拌摩擦焊接过程产热机理,分别按下式计算Qsb、Qss、Qpb、Qps:
采用热流通量指标表示产热模型,热流通量包括面热流通量和体热流通量;
设搅拌针和工件摩擦产生的热能被均匀地施加在搅拌针的体积范围内,由Qpb、Qps得到体热流通量qvolume为:
设搅拌头轴肩与工件摩擦产生的热能被线性的热流通量施加在工件表面,设r是工件表面上的点到搅拌头轴线的距离,且R1≤r≤R2;由H、h和hg的关系,以及Qsb、Qss,得到面热流通量qsurface(r)并化简为:
进一步地,设x为一组待优化工艺参数,约束条件包括:
2)设工件与搅拌头轴肩边缘接触区域的最高温度的约束条件为g2(x),g2(x)=Tmax(x)-Tm≤0;式中,Tmax(x)为工件与搅拌针轴肩边缘接触区域的最高温度;Tm为工件材料的熔点;
3)设工件中间厚度的最高温度约束条件为g3(x),g3(x)=Tmd(x)-Tw≥0;式中,Tmd(x)为工件中间厚度区域的最高温度;Tw为焊接温度。
进一步地,设v为搅拌头进给焊接速度,设S为工件与搅拌头轴肩接触面积,设f(x)为双侧同时搅拌摩擦焊焊接单位长度的输入能量,则f(x)表示如下:
进一步地,由约束条件和目标函数,建立如下优化模型:
min f(x);
s.t.xl≤x≤xu;
g2(x)=Tmax(x)-Tm≤0;
g3(x)=Tmd(x)-Tw≥0;
其中,xl和xu为参数x的上下边界。
进一步地,采用序列二次规划算法求解优化模型的方法包括如下步骤:
步骤1,输入初始工艺参数x0,计算初始能量f0(x),并令i=0;
步骤2,令i=i+1,采用序列二次规划算法确定变量的参数值xi;
步骤3,计算当前迭代步的面热流通量qsurface(r)和体热流通量qvolume;
步骤4,由傅里叶传热定律和能量守恒定律,计算得到双侧同时搅拌摩擦焊焊接过程中工件的温度分布;
步骤5,由焊接过程中工件的温度分布,判断是否满足约束条件;不满足约束条件,令i=i-1,返回步骤2,满足约束条件进行步骤6;
步骤6,计算目标函数fi(x);
步骤7,判断是否满足迭代终止条件:i>imax或Δf=fi+1(x)-fi(x)<103;若满足,迭代过程终止,输出当前第i代的x为最优工艺参数;否则,转至步骤2。
进一步地,步骤4中包括如下具体步骤:
由傅里叶传热定律和能量守恒定律得到如下控制方程:
其中其中ρ(T)、k(T)和cp(T)分别为与温度T相关的工件材料密度、工件导热率和工件比热容;q(x,y,z,t)是变量在工件单位体积内产生的与坐标(x,y,z)和时间t相关的热量函数;
将体热流通量qvolume、面热流通量qsurface(r)、环境温度T0、系数h∞和系数作为控制方程的边界条件加载到有限元模型中,采用数值法求解控制方程,得到双侧同时搅拌摩擦焊焊接过程中工件的温度分布。
本发明具有的优点和积极效果是:本发明所提出的半解析产热模型比较全面的考虑了双侧同时搅拌摩擦焊焊接过程中的产热因素,基于该产热过程建立了自适应的热源模型,能够准确的预测双侧同时搅拌摩擦焊焊接过程中的温度场。
本发明提出了一种基于半解析产热模型的双侧同时搅拌摩擦焊工艺参数优化方法,该方法基于焊接过程的产热机理,与SQP优化算法相耦合,得到了全局最优的输入能量和最佳的工艺参数组合,优化后的参数提高了焊接的效率和改善了焊接接头的质量。该方法整体过程操作简单,可方便高效的针对其它不同材料和尺寸的工件进行双侧同时搅拌摩擦焊工艺参数优化,为提高双侧同时搅拌摩擦焊在工程中的更近一步应用提供理论指导。
附图说明
图1为双侧同时搅拌摩擦焊接的工作原理示意图。
图2为双侧同时搅拌摩擦焊接过程中搅拌头产热热量分布图。
图3为本发明的一种工件表面热影响区的宽度示意图。
图4为采用本发明的方法计算得到的热循环与实验中采用热电偶测得的热循环结果对比图。
图中:1、待焊接的铝合金板之一;2、待焊接的铝合金板之二;3、搅拌头;4、工装夹具;5、焊缝;6、搅拌头轴肩;7、搅拌针;8、热影响区;9、热机影响区和焊核区。
Qsb为搅拌头轴肩端面与工件摩擦产生的热功率,Qss为搅拌头轴肩侧面与工件摩擦产生的热功率,Qpb为搅拌针端面与工件摩擦产生的热功率,Qps为搅拌针侧面与工件摩擦产生的热功率,R2为搅拌头轴肩半径,R1为搅拌针半径,hp为搅拌针长度,h为搅拌头轴肩下压量,hg为两侧搅拌针端面之间的距离,H为工件整体厚度,L为工件整体长度,W为工件整体宽度,ω为搅拌头的转速,σyield为工件材料屈服应力,τcontact为搅拌头与工件接触点产生的剪切应力;P1(L/2,yd,0)为实验验证的参考点,P2(xi,R1,0)、P3(xi,R1,-H/2),j=1,2,3为优化过程计算约束函数的两个参考点,LHAZ为板件表面热影响区的宽度。
具体实施方式
为能进一步了解本发明的发明内容、特点及功效,兹列举以下实施例,并配合附图详细说明如下:
请参见图1至图4,一种双侧同时搅拌摩擦焊工艺参数优化方法,基于搅拌摩擦焊接过程产热机理、傅里叶传热定律和能量守恒定律,建立双侧同时搅拌摩擦焊的产热模型,确定待优化工艺参数,以焊接单位长度输入能量最小为目标函数,根据焊接工艺设置约束条件,建立优化模型;采用序列二次规划算法求解优化模型,得到最优工艺参数值。
优化地,待优化工艺参数可包括:搅拌头3的转速、搅拌头3的进给焊接速度以及两侧搅拌针7端面之间的距离。
优化地,建立双侧同时搅拌摩擦焊的产热模型的方法可包括:
由于两侧搅拌头3关于工件厚度方向几何对称且两侧的搅拌头3的焊接工艺参数相同,因此两侧搅拌头3的产热模型相同,可只建立其中一侧搅拌头3的产热模型;
根据搅拌头3与工件的接触面不同,可将搅拌摩擦焊接过程产生的热功率分成以下四种:分别是搅拌头轴肩6端面与工件摩擦产生的热功率Qsb、搅拌头轴肩6侧面与工件摩擦产生的热功率Qss、搅拌针7端面与工件摩擦产生的热功率Qpb、搅拌针7侧面与工件摩擦产生的热功率Qps:
可设:R2为搅拌头轴肩6半径,R1为搅拌针7端部半径,hp为搅拌针7长度,h为轴肩下压量,hg为两侧搅拌针7端面之间的距离,H为工件厚度,ω为搅拌头3的转速,σyield为工件材料的屈服应力,基于搅拌摩擦焊接过程产热机理,可分别按下式计算Qsb、Qss、Qpb、Qps:
可采用热流通量指标表示产热模型,热流通量包括面热流通量和体热流通量;
可设搅拌针7和工件摩擦产生的热能被均匀地施加在搅拌针7的体积范围内,由Qpb、Qps得到体热流通量qvolume为:
可设搅拌头轴肩6与工件摩擦产生的热能被线性的热流通量施加在工件表面,设r是工件表面上的点到搅拌头3轴线的距离,且R1≤r≤R2;由H、h和hg的关系,以及Qsb、Qss,得到面热流通量qsurface(r)并化简为:
优化地,可设x为一组待优化工艺参数,约束条件可包括:
2)设工件与搅拌头轴肩6边缘接触区域的最高温度的约束条件为g2(x),g2(x)=Tmax(x)-Tm≤0;式中,Tmax(x)为工件与搅拌针7轴肩边缘接触区域的最高温度;Tm为工件材料的熔点;
3)设工件中间厚度的最高温度约束条件为g3(x),g3(x)=Tmd(x)-Tw≥0;式中,Tmd(x)为工件中间厚度区域的最高温度;Tw为焊接温度。
优化地,可设v为搅拌头3进给焊接速度,设S为工件与搅拌头轴肩6接触面积,设f(x)为双侧同时搅拌摩擦焊焊接单位长度的输入能量,则f(x)可表示如下:
优化地,由约束条件和目标函数,可建立如下优化模型:
min f(x);
s.t.xl≤x≤xu;
g2(x)=Tmax(x)-Tm≤0;
g3(x)=Tmd(x)-Tw≥0;
其中,xl和xu为参数x的上下边界。
优化地,采用序列二次规划算法求解优化模型的方法可包括如下步骤:
步骤1,输入初始工艺参数x0,计算初始能量f0(x),并令i=0;
步骤2,令i=i+1,采用序列二次规划算法确定变量的参数值xi;
步骤3,计算当前迭代步的面热流通量qsurface(r)和体热流通量qvolume;
步骤4,由傅里叶传热定律和能量守恒定律,计算得到双侧同时搅拌摩擦焊焊接过程中工件的温度分布;
步骤5,由焊接过程中工件的温度分布,判断是否满足约束条件;不满足约束条件,令i=i-1,返回步骤2,满足约束条件进行步骤6;
步骤6,计算目标函数fi(x);
步骤7,判断是否满足迭代终止条件:i>imax或Δf=fi+1(x)-fi(x)<103;若满足,迭代过程终止,输出当前第i代的x为最优工艺参数;否则,转至步骤2。imax为设定的最大迭代次数。
优化地,步骤4中可包括如下具体步骤:
由傅里叶传热定律和能量守恒定律得到如下控制方程:
其中其中ρ(T)、k(T)和cp(T)分别为与温度T相关的工件材料密度、工件导热率和工件比热容;q(x,y,z,t)是变量在工件单位体积内产生的与坐标(x,y,z)和时间t相关的热量函数;
可将体热流通量qvolume、面热流通量qsurface(r)、环境温度T0、系数h∞和系数作为控制方程的边界条件加载到有限元模型中,采用数值法求解控制方程,得到双侧同时搅拌摩擦焊焊接过程中工件的温度分布。
下面以本发明的一个优选实施例来进一步说明本发明的工作流程及工作原理:
一种双侧同时搅拌摩擦焊工艺参数优化方法,该方法包括如下步骤:
步骤1:建立双侧同时搅拌摩擦焊的半解析产热模型
步骤1.1:产热模型假设
研究的模型采用双侧同时搅拌摩擦焊焊接方式,由材料相同的两块铝合金板通过双侧同时搅拌摩擦焊工艺焊接而成,如图1所示,将待焊接的铝合金板之一1和待焊接的铝合金板之二2,拼贴在一起,在铝合金板两侧面同时采用搅拌头3进行搅拌焊接,焊接后整体工件尺寸大小为L×W×H,两个搅拌头3正对放置,hg为两侧搅拌针7端面之间的距离,旋转方向相反,以相同的旋转速度ω和进给速度v沿着焊缝5线向x轴正方向移动。为了降低求解热模型的复杂性,进行了一些必要的假设;整个焊接过程中工件的热传导行为是三维瞬态的、不考虑搅拌头3的倾斜角度、仅对焊接过程的准稳态阶段进行建模、不考虑焊接过程的热辐射。
步骤1.2:控制方程和边界条件
由傅里叶传热定律和能量守恒定律得出热传导的基本方程为:
其中ρ(T)、k(T)和cp(T)是工件材料与温度T相关的密度、导热率和比热容。q(x,y,z,t)是单位体积内产生的热量(是热源在工件中的坐标(x,y,z)和时间t的函数)。为了求解控制方程,需给出焊接工件的边界传热条件,由于焊接工件关于焊缝5几何对称且两块板的材料相同,只对其中一块工件建模,焊接工件表面、工装夹具4与周围环境发生热对流的系数分别用h∞和表示,环境温度为T0,搅拌头3和工件的产热用热流通量q表示,该热流通量由面热流通量和体热流通量组成。
步骤1.3:建立双侧同时搅拌摩擦焊的解析产热模型
焊接过程中的热量分布图如图2所示,图中,R2是搅拌头轴肩6半径,R1和hp是搅拌针7半径和长度,h是搅拌头轴肩6下压量。由于两侧搅拌头3关于工件厚度方向几何对称,以一侧搅拌头3的产热为例进行建模。根据搅拌头3与工件的接触面不同,产热主要分为以下四类热源:搅拌头轴肩6端面与工件接触产生的热功率Qsb,搅拌头轴肩6侧面与工件接触产生的热功率Qss,搅拌针7端面与工件接触产生的热功率Qpb,搅拌针7侧面与工件接触产生的热功率Qps,(在沿焊缝5移动过程中,只有前进侧与工件接触)。焊接过程中工具和工件接触点产生的剪切应力τcontact为:
τcontact=(1-δ)τfriction+δτyield
其中,δ是接触状态变量,τfriction是材料的摩擦剪切应力,τyield是材料的屈服剪切应力,根据von Mises屈服准则,材料的屈服剪切应力τyield与屈服应力σyield存在如下关系:
当焊接过程处于准稳态状态下,摩擦剪切应力和屈服剪切应力是一对平衡力,即τfriction=τyield,因此工具和工件接触点产生的剪切应力τcontact化简为:
通过对搅拌头3和工件接触面的剪切应力τcontact和接触面积进行积分,即可得到各个热源的解析表达式并化简为:
搅拌针7和工件产热的热源被均匀地施加在搅拌针7的体积范围内,由以上产热方程得出体热流通量qvolume为:
搅拌头轴肩6与工件接触面产生的热源被线性的热流通量施加在工件表面,由以上产热方程和轴肩下压量h和搅拌针7之间距离hg的关系h=(H-hg-2hp)/2,得出面热流通量qsurface(r)并化简为:
qsurface(r)=kr (3)
式中,k是比例系数,设r是工件表面上的点到搅拌头3轴线的距离,且R1≤r≤R2;S是搅拌头轴肩6端面从半径R1到半径R2的面积。
将以上求出的体热源qvolume、面热源qsurface(r)、环境温度T0、对流边h∞和作为控制方程(1)的边界条件加载到有限元模型中(另一侧搅拌头3采用相同的建立热源的方法即可得到该侧搅拌头3施加在工件中的面热源和体热源),采用数值法求解控制方程(1),得到双侧同时搅拌摩擦焊焊接过程中工件的温度分布。
步骤2:建立基于半解析产热模型的双侧同时搅拌摩擦焊工艺参数的优化过程步骤2.1建立优化变量和目标
在双侧同时搅拌摩擦焊焊接过程中,搅拌头3的尺寸R1、R2、hp和工件的尺寸L×W×H为提前设定好的常数值,而转速ω、焊接速度v和相对两侧搅拌头之间的距离hg是可调节的变量,且这些变量是影响焊接过程中热输入量的主要因素,因此待优化的工艺参数是ω、v和hg,以向量x表示为:
x=(ω,v,hg)
目标需要被定义为描述焊接过程中能量消耗的函数,由焊接过程中产生的总热量除以焊接速度v即可得到优化过程中的目标函数为:
步骤2.2:建立约束条件
搅拌摩擦焊焊后工件的显微硬度在工件上整体呈“W“形状,临近热机影响区和焊核区9,为热影响区8,热影响区8的硬度最低,其在拉伸过程中最容易发生断裂,所以,热影响区8为危险区域,优化过程应使该区域最小化。根据半解析产热模型的热循环求解结果,在迭代优化过程中,定义该区域的约束条件为g1(x),其表达式为:
焊接过程中,搅拌头轴肩6边缘与工件接触区域的最高温度Tmax(x)必须低于材料的熔点Tm,定义该状态的约束条件为g2(x),其表达式为:
g2(x)=Tmax(x)-Tm≤0
为了确保搅拌头3沿着焊缝5顺利移动,在搅拌针7端部周围区域的材料应足够糊状,该工件中间厚度位置的最高温度Tmd(x)必须高于或等于焊接温度Tw,该状态的约束条件为g3(x),其表达式为:
g3(x)=Tmd(x)-Tw≥0
因此,根据以上约束条件,得到所要求解的优化问题如下:
min f(x)
s.t.xl≤x≤xu
g2(x)=Tmax(x)-Tm≤0
g3(x)=Tmd(x)-Tw≥0
其中,xl和xu为参数x的上下边界。
步骤2.3:优化算法流程
基于以上建立的约束条件和目标函数,耦合序列二次规划算法(SQP)来求解双侧同时搅拌摩擦焊中的最优工艺参数值,当算法的迭代次数i超过最大迭代次数imax=50或者相邻两迭代步的能量函数差值Δf<103时,迭代过程终止。得出优化算法的流程如下:
①输入初始工艺参数x0,计算初始能量f0(x),并令i=0。
②令i=i+1,采用SQP算法确定当前迭代步i中变量的参数值为xi。
③采用方程(2)和(3)计算当前迭代步i中的面热流通量qsurface(r)和体热流通量qvolume。
采用同样的建模方法得到的另一侧搅拌头3的面热流通量qsurface(r)和体热流通量qvolume。
④将步骤③得到的面热流通量qsurface(r)和体热流通量qvolume,作为控制方程(1)中的边界条件,通过有限元软件ABAQUS的DFLUX子程序,加载双侧搅拌摩擦焊的热流密度到产热模型中求解控制方程(1),得到焊接过程中的温度分布。
⑤根据步骤④得到的结果计算约束函数g1(x)、g2(x)、g3(x)和目标函数fi(x)。
⑥判断是否满足迭代终止条件:i>imax或Δf=fi+1(x)-fi(x)<103。若满足,迭代过程终止,输出当前迭代步i的变量x为最优工艺参数。否则,转至步骤②。
步骤3:基于半解析产热模型的双侧同时搅拌摩擦焊工艺参数优化方法分析
步骤3.1:采用实验方法验证半解析产热模型的正确性
待焊接工件的整体几何尺寸为250mm*150mm*4mm,搅拌头3的尺寸为;搅拌头轴肩6的半径R2=6.5mm,搅拌针7半径R1=2.5mm,搅拌针7长度hp=1.8mm,材料为铝合金Al2024,对流边界h∞=10W.m-2.K-1、初始焊接的参数为;搅拌头3旋转速度ω=1100rpm,进给速度v=21mm/min,相对两侧搅拌头3的间隙为hg=0mm。实验中采用热电偶测量工件的温度循环结果,热电偶的位置为P1(125,20,0)。
为了便于对优化过程的结果进行分析,设置了三个参考点P1(L/2,yd,0)、P2(xi,R1,0)和P3(xi,R1,-H/2),j=1,2,3。如图3所示,P1点代表了实验中焊接工件热电偶的位置,将实验得到的结果与半解析产热模型的热循环结果进行对比分析。P2点和P3点分别代表搅拌头轴肩6端部和工件中间厚度的位置,目的是为了便于计算优化问题中的约束条件g1(x)、g2(x)、g3(x),其中x1=L/4,x2=L/2,x3=3L/4。
步骤3.2:在建立半解析产热模型的基础上,通过与序列二次规划算法(SQP)相耦合,采用步骤2.3的工艺参数优化流程即可得到双侧同时搅拌摩擦焊最佳的工艺参数。
采用本发明的优化方法进行优化同时双侧搅拌摩擦焊最佳的工艺参数,热电偶测出的温度结果和本发明的半解析产热模型的热循环结果对比如图4所示,从图中看出,半解析产热模型的温度循环结果和热电偶的实验结果峰值温度相近且温度上升和下降具有较好的一致性,从而验证了该半解析产热模型的正确性,可精确的用于温度分布和冷却速率的预测。工艺参数优化后的结果如下表1所示,从表中看出,优化后的焊接热输入能量低于初始输入能量,TP1表示峰值温度,LHAZ表示热影响区8的长度,热影响区8的长度小于初始参数对应的热影响区8长度,且优化后的焊接焊接速度相对于初始速度有较大的提高。
表1优化结果
以上所述的实施例仅用于说明本发明的技术思想及特点,其目的在于使本领域内的技术人员能够理解本发明的内容并据以实施,不能仅以本实施例来限定本发明的专利范围,即凡本发明所揭示的精神所作的同等变化或修饰,仍落在本发明的专利范围内。
Claims (8)
1.一种双侧同时搅拌摩擦焊工艺参数优化方法,其特征在于,基于搅拌摩擦焊接过程产热机理、傅里叶传热定律和能量守恒定律,建立双侧同时搅拌摩擦焊的产热模型,确定待优化工艺参数,以焊接单位长度输入能量最小为目标函数,根据焊接工艺设置约束条件,建立优化模型;采用序列二次规划算法求解优化模型,得到最优工艺参数值。
2.根据权利要求1所述的双侧同时搅拌摩擦焊工艺参数优化方法,其特征在于,待优化工艺参数包括:搅拌头的转速、搅拌头的进给焊接速度以及两侧搅拌针端面之间的距离。
3.根据权利要求1所述的双侧同时搅拌摩擦焊工艺参数优化方法,其特征在于,建立双侧同时搅拌摩擦焊的产热模型的方法包括:
由于两侧搅拌头关于工件厚度方向几何对称且两侧的搅拌头的焊接工艺参数相同,因此两侧搅拌头的产热模型相同,只建立其中一侧搅拌头的产热模型;
根据搅拌头与工件的接触面不同,将搅拌摩擦焊接过程产生的热功率分成以下四种:分别是搅拌头轴肩端面与工件摩擦产生的热功率Qsb、搅拌头轴肩侧面与工件摩擦产生的热功率Qss、搅拌针端面与工件摩擦产生的热功率Qpb、搅拌针侧面与工件摩擦产生的热功率Qps:
设:R2为搅拌头轴肩半径,R1为搅拌针端部半径,hp为搅拌针长度,h为轴肩下压量,hg为两侧搅拌针端面之间的距离,H为工件厚度,ω为搅拌头的转速,σyield为工件材料的屈服应力,基于搅拌摩擦焊接过程产热机理,分别按下式计算Qsb、Qss、Qpb、Qps:
采用热流通量指标表示产热模型,热流通量包括面热流通量和体热流通量;
设搅拌针和工件摩擦产生的热能被均匀地施加在搅拌针的体积范围内,由Qpb、Qps得到体热流通量qvolume为:
设搅拌头轴肩与工件摩擦产生的热能被线性的热流通量施加在工件表面,设r是工件表面上的点到搅拌头轴线的距离,且R1≤r≤R2;由H、h和hg的关系,以及Qsb、Qss,得到面热流通量qsurface(r)并化简为:
7.根据权利要求6所述的双侧同时搅拌摩擦焊工艺参数优化方法,其特征在于,采用序列二次规划算法求解优化模型的方法包括如下步骤:
步骤1,输入初始工艺参数x0,计算初始能量f0(x),并令i=0;
步骤2,令i=i+1,采用序列二次规划算法确定变量的参数值xi;
步骤3,计算当前迭代步的面热流通量qsurface(r)和体热流通量qvolume;
步骤4,由傅里叶传热定律和能量守恒定律,计算得到双侧同时搅拌摩擦焊焊接过程中工件的温度分布;
步骤5,由焊接过程中工件的温度分布,判断是否满足约束条件;不满足约束条件,令i=i-1,返回步骤2,满足约束条件进行步骤6;
步骤6,计算目标函数fi(x);
步骤7,判断是否满足迭代终止条件:i>imax或Δf=fi+1(x)-fi(x)<103;若满足,迭代过程终止,输出当前第i代的x为最优工艺参数;否则,转至步骤2。
8.根据权利要求7所述的双侧同时搅拌摩擦焊工艺参数优化方法,其特征在于,步骤4中包括如下具体步骤:
由傅里叶传热定律和能量守恒定律得到如下控制方程:
其中其中ρ(T)、k(T)和cp(T)分别为与温度T相关的工件材料密度、工件导热率和工件比热容;q(x,y,z,t)是变量在工件单位体积内产生的与坐标(x,y,z)和时间t相关的热量函数;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110986364.1A CN113579462B (zh) | 2021-08-26 | 2021-08-26 | 一种双侧同时搅拌摩擦焊工艺参数优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110986364.1A CN113579462B (zh) | 2021-08-26 | 2021-08-26 | 一种双侧同时搅拌摩擦焊工艺参数优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113579462A true CN113579462A (zh) | 2021-11-02 |
CN113579462B CN113579462B (zh) | 2022-11-04 |
Family
ID=78239365
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110986364.1A Active CN113579462B (zh) | 2021-08-26 | 2021-08-26 | 一种双侧同时搅拌摩擦焊工艺参数优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113579462B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004216435A (ja) * | 2003-01-16 | 2004-08-05 | Mitsubishi Heavy Ind Ltd | 摩擦攪拌接合装置とその接合方法 |
CN101722360A (zh) * | 2009-12-29 | 2010-06-09 | 湖北工业大学 | 一种搅拌摩擦点焊方法 |
CN102601516A (zh) * | 2012-03-29 | 2012-07-25 | 哈尔滨工业大学 | 双面对称搅拌摩擦焊接方法 |
CN104203481A (zh) * | 2012-04-06 | 2014-12-10 | 杰富意钢铁株式会社 | 钢板的摩擦搅拌接合方法 |
CN109783921A (zh) * | 2019-01-07 | 2019-05-21 | 中石化石油工程技术服务有限公司 | 油气输送管道的热影响区评估方法、装置及计算机设备 |
CN110008554A (zh) * | 2019-03-27 | 2019-07-12 | 哈尔滨工业大学 | 一种基于数值模拟和深度学习的搅拌摩擦焊缝成形预测优化方法 |
CN112548308A (zh) * | 2020-11-12 | 2021-03-26 | 广州瑞松智能科技股份有限公司 | 一种铝与钢搅拌摩擦焊焊接变形预测方法、电子设备、存储介质 |
-
2021
- 2021-08-26 CN CN202110986364.1A patent/CN113579462B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004216435A (ja) * | 2003-01-16 | 2004-08-05 | Mitsubishi Heavy Ind Ltd | 摩擦攪拌接合装置とその接合方法 |
CN101722360A (zh) * | 2009-12-29 | 2010-06-09 | 湖北工业大学 | 一种搅拌摩擦点焊方法 |
CN102601516A (zh) * | 2012-03-29 | 2012-07-25 | 哈尔滨工业大学 | 双面对称搅拌摩擦焊接方法 |
CN104203481A (zh) * | 2012-04-06 | 2014-12-10 | 杰富意钢铁株式会社 | 钢板的摩擦搅拌接合方法 |
CN109783921A (zh) * | 2019-01-07 | 2019-05-21 | 中石化石油工程技术服务有限公司 | 油气输送管道的热影响区评估方法、装置及计算机设备 |
CN110008554A (zh) * | 2019-03-27 | 2019-07-12 | 哈尔滨工业大学 | 一种基于数值模拟和深度学习的搅拌摩擦焊缝成形预测优化方法 |
CN112548308A (zh) * | 2020-11-12 | 2021-03-26 | 广州瑞松智能科技股份有限公司 | 一种铝与钢搅拌摩擦焊焊接变形预测方法、电子设备、存储介质 |
Non-Patent Citations (2)
Title |
---|
刘向前: "铝合金型材高速搅拌摩擦焊接头结构设计与组织性能研究", 《工程科技Ⅰ辑》 * |
肖毅华等: "6061-T6铝合金搅拌摩擦焊温度场的数值模型和参数影响分析", 《机械科学与技术》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113579462B (zh) | 2022-11-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Buffa et al. | A continuum based fem model for friction stir welding—model development | |
Sun et al. | A comparative study on welding temperature fields, residual stress distributions and deformations induced by laser beam welding and CO2 gas arc welding | |
Ulysse | Three-dimensional modeling of the friction stir-welding process | |
Li et al. | Numerical Simulation of Friction Welding Processes Based on ABAQUS Environment. | |
Boukraa et al. | Friction stir welding process improvement through coupling an optimization procedure and three-dimensional transient heat transfer numerical analysis | |
CN110530541B (zh) | 一种可精确模拟大型压力容器焊后热处理温度场计算方法 | |
Darvazi et al. | Thermal modeling of friction stir welding of stainless steel 304L | |
Israr et al. | Simulation of different path strategies for wire-arc additive manufacturing with Lagrangian finite element methods | |
Nguyen et al. | Analysis and compensation of shrinkage and distortion in wire-arc additive manufacturing of thin-walled curved hollow sections | |
Li et al. | Using coupled Eulerian Lagrangian formulation for accurate modeling of the friction stir welding process | |
Liu et al. | A modified analytical heat source model for numerical simulation of temperature field in friction stir welding | |
Klobčar et al. | Finite element modeling of GTA weld surfacing applied to hot-work tooling | |
CN113579462B (zh) | 一种双侧同时搅拌摩擦焊工艺参数优化方法 | |
Xie et al. | In-depth understanding of rotating toolhead-induced heat generation and material flow behavior in friction-rolling additive manufacturing | |
Ahmad et al. | Numerical modelling and experimental analysis on angular strain induced by bead-on-plate SS316L GMAW using inherent strain and thermomechanical methods | |
Abhilash et al. | Simulation of friction stir welding using thermo-mechanical coupled finite element method | |
Leon et al. | Numerical modeling of thermal field during friction stir welding using tool with polygonal pin profile | |
Li et al. | Simplified thermo-mechanical modeling of friction stir welding with a sequential FE method | |
Das et al. | Prediction of surface profile in friction stir welding using coupled eulerian and lagrangian method | |
Chaudhary | Modeling of laser-additive manufacturing processes | |
Meng et al. | A method for process parameter optimization of simultaneous double-sided friction stir welding using a heat transfer model | |
Arif et al. | Finite element modelling for validation of maximum temperature in friction stir welding of aluminium alloy | |
Prajadhiana et al. | Investigation on Welded T-Joint Distortion Using Virtual Manufacturing Tools with Simplified Procedure | |
Deo et al. | Experimental verification of distortion analysis of welded stiffeners | |
Hasan | Modelling of Tool Wear and Metal Flow Behaviour in Friction Stir Welding (FSW). |
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 |