CN110008554A - 一种基于数值模拟和深度学习的搅拌摩擦焊缝成形预测优化方法 - Google Patents
一种基于数值模拟和深度学习的搅拌摩擦焊缝成形预测优化方法 Download PDFInfo
- Publication number
- CN110008554A CN110008554A CN201910238706.4A CN201910238706A CN110008554A CN 110008554 A CN110008554 A CN 110008554A CN 201910238706 A CN201910238706 A CN 201910238706A CN 110008554 A CN110008554 A CN 110008554A
- Authority
- CN
- China
- Prior art keywords
- soldering set
- network
- weld
- input
- parameter
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提出了一种基于数值模拟和深度学习的搅拌摩擦焊缝成形预测优化方法,属于搅拌摩擦焊技术领域。所述方法包括:步骤一、设定三次模拟试验作为数据测试集;步骤二、计算焊接过程中材料流动场和温度场分布情况;步骤三、计算不同参数下搅拌摩擦焊具断裂失效情况,并计算不同参数下焊缝成形质量及缺陷分布情况;步骤四、利用生成对抗网络深度学习模型遍历所有工艺参数和焊具结构的焊缝成形结果,获得在保证焊具可靠工作前提下焊缝成形最优化结果。本发明所述方法目的在于为生产中最优化搅拌摩擦焊工艺提供有效的普适性预测方法,具有降低时间消耗、材料成本以及预测精确度高等优点。
Description
技术领域
本发明涉及一种基于数值模拟和深度学习的搅拌摩擦焊缝成形预测优化方法,属于搅拌摩擦焊技术领域。
背景技术
搅拌摩擦焊(Friction Stir Welding)作为自激光焊接问世以来最引人注目的焊接方法,在克服传统熔焊轻质合金过程中易产生的裂纹和气孔等焊接缺陷问题上具有显著优势,使难于熔焊的部分材料实现了可靠的优质连接,其工作基本原理是用一个具有轴肩(或带有搅拌针)的焊具进行高速旋转,工具与工件之间产生大量的摩擦热和塑性变形热,使工件材料局部塑化,当焊具沿着焊接界面向前移动时,塑化材料在焊具的作用下用前部移向后部,并在焊具的挤压作用下形成致密的固相接头。
然而,尚不成熟的搅拌摩擦焊存在着如下缺点:(1)焊具设计尚无普适性理论。在搅拌摩擦焊过程中,被焊材料主要受焊具驱动流动并形成致密焊缝结构,同时焊具会受到流动的被焊材料的强阻力作用,不可靠的焊具设计会导致焊缝成形差、焊具折断等问题的发生;(2)目前关于搅拌摩擦焊工艺参数优化的研究主要以实验为主,然而由于搅拌摩擦热力耦合系统本身的复杂性、焊缝成形设计的因素过多、正交实验设计次数较多、工作量巨大且重复性较弱等,优化预测的实验研究和预测成本极高,实际工程应用中参数优化过程通常较为简陋。因此,工程上对经济、便捷且普适的基于数值模拟和深度学习的搅拌摩擦焊缝成形预测优化方法有强烈需求。
发明内容
本发明针对上述搅拌摩擦焊工艺参数和焊具结构优化存在的不足,提供一种基于数值模拟和深度学习的搅拌摩擦焊缝成形预测优化方法,目的在于为生产中最优化搅拌摩擦焊工艺提供有效的普适性预测方法。由于数值模拟和深度学习只需依托于一般的计算机硬件和相应的软件平台,所以其可以经济地实现对不同工艺参数和焊具结构的快捷测试对比,实现焊缝成形最优化;同时,结合充分的实验数据作为深度学习测试集,整个预测优化系统的可靠性得到了有力保障,为达到上述目的,本发明所采用的技术方法如下:
一种基于数值模拟和深度学习的搅拌摩擦焊缝成形预测优化方法,所述搅拌摩擦焊缝成形预测优化方法包括:
步骤一、可选地,根据实际生产需求,初步设置至少三次针对搅拌摩擦焊缝成形的模拟试验,作为数值模拟模型的数据基础和生成对抗网络(Generative AdversarialNetwork)深度学习模型的测试集,然后生成对抗网络深度学习模型;
步骤二、根据给定的实际工况及材料物性参数建立用于数值模拟模型,计算焊接过程中材料流动场和温度场分布情况;
步骤三、根据步骤二获得的材料流动场和温度场分布情况,结合单向流固耦合模型,计算不同参数下搅拌摩擦焊具断裂失效情况;然后,结合离散相模型,计算不同参数下焊缝成形质量及缺陷分布情况;
步骤四、根据上述计算所得焊缝成形质量及缺陷分布情况和焊具断裂失效情况,利用生成对抗网络深度学习模型遍历所有工艺参数和焊具结构的焊缝成形结果,从而获得在保证焊具可靠工作前提下焊缝成形最优化结果。
进一步地,所述模拟实验所采集的参数包括:焊具结构,焊接参数,焊具断裂失效状况,焊缝成形质量结果,且其中应包括至少一组用以完成焊接的参数。
进一步地,所述搅拌摩擦焊具断裂失效情况的确定是采用流体动力学方法,将被焊材料视为非牛顿、不可压缩的连续流体,并采用有限体积法进行计算,计算软件可选地为Fluent、STAR-CCM、OpenFOAM等流体动力学仿真软件;所述搅拌摩擦焊具断裂失效情况的数值计算的步骤如下:
第一步、建立被焊材料流域模型:建立符合实际生产需求结构的被焊材料流域模型,并利用滑移网格将所述焊接材料流域分割为圆柱状的滑移网格部分和除滑移网格部分外的固定部分;所述滑移网格部分旋转速率等同于焊具旋转速率;
第二步、设置被焊材料物性参数,所述被焊材料物性参数包括被焊材料的密度、流动应力、导热率和定压比热容;
第三步、确定近焊具壁面材料流动准则:将所述近焊具壁面材料流动准则分解为两部分:一部分是竖直速度vvertical,一部分是水平速度vhorizon;然后将所述水平速度vhorizon分为剪切速度vshear和法向速度vnormal;
第四步、设定产热模型:将所述产热模型分解为两部分,所述两部分包括:摩擦热和变形热,其中,摩擦热表达粘着摩擦过程中伴随的能量转换和耗散,变形热为材料层间剪切产生;
第五步、利用牛顿换热原理和蒂芬-玻尔兹曼定律确定所述被焊材料流域的流域边界换热条件,本步骤中,流域边界换热条件可由牛顿换热公式和斯蒂芬-玻尔兹曼定律给出,速度边界条件可视为完全自由流动或完全约束流动;
第六步、从所述流域的入口添加离散相模型,所述离散相模型的输入入口同于流域速度入口,离散相为无质量粒子;
第七步、利用单向流固耦合模型计算焊具断裂失效状况,
第八步、提交运算并输出结果,数值计算完成后,即可获取不同焊接参数和焊具结构下的焊缝成形图像和焊具断裂失效状况。
进一步地,第二步所述流动应力由Sellars-Tegart材料本构模型计算,所述Sellars-Tegart材料本构模型为:
其中,σe为材料流动应力,T为材料当前温度,Tmelt为材料熔点(固相线),α、A和n为材料特征常数,σmelt为熔化态材料流动应力,Z为Zener-Hollomon系数,公式为:
其中,为材料当前应变速率,Q为材料活化能,R为气体常数,e为自然指数。
进一步地,第三步所述竖直速度vvertical、水平速度vhorizon、剪切速度vshear和法向速度vnormal的计算公式为:
其中,ijk代指焊具上各微元体;ω为焊具旋转速度;r表示微元体距焊具轴线距离;ξ表示轴向方位角;θ表示切向方位角;η为滑移系数,即壁面材料流动速率与焊具表面线速度的比值,计算方式如下:
其中,ijk代指焊具上各微元体;η0为滑移归一化常数;e表示自然指数;ω表示角速度;r表示微元体距焊具轴线距离;ω0表示角速度归一化常数;rs表示参考半径。
进一步地,第四步所述产热模型表示为:
其中,ijk代指焊具上各微元体;F1、F2分别表示为摩擦热输入和形变热输入;β1和β2为经验参数,取值为0.1~1.0;r表示微元体距焊具轴线距离;ω表示角速度;η表示滑移系数;σe表示微元体等效应力;μ表示动摩擦因数;P表示焊接下压力;ξ表示轴向方位角;θ表示切向方位角;表示等效应变速率。
进一步地,所述生成对抗网络深度学习模型的过程如下:
步骤1、建立图像标签库:将数值模拟结果所得焊缝成形图像与焊接断裂失效状况进行分别和标注,每张图像具有唯一与之对应的焊接参数和焊具结构标签k;
步骤2、构建生成器网络G:在所述生成器网络G中,将参数标签k与随机噪声组合为新的向量输入,将所述生成器网络G生成的数据作为判别器网络D的输入,生成器网络G包括第一输入层、第一全连接层、第一反卷积层、第一批归一化层和第一激活层,以随机的噪声信号和类别标签组合而成的向量作为所示第一输入层的输入,所述输入经过所述第一全连接层、第一反卷积层、第一批归一化层和第一激活层生成数据,其中反卷积次数可选地为4~6次,以生成不低于三通道128*128分辨率的焊缝成形特征图,并将其作为判别器网络D的输入;
步骤3、构建判别器网络D:所述判别器网络D的损失函数由用于判别图像真假的第一损失函数和用于判别图像类别的第二损失函数组成,叛变期网络D包括第二输入层、第二反卷积层、第二批归一化层、第二激活层和第二全连接层,分别以步骤2中所示生成器网络K生成的数据和数值模拟所得焊缝成形图像作为所示第二输入层的输入,所述第二全连接层由用于判断图像真假的全连接层和判断图像类别的全连接层组成,所述第二输入层的输入经过所述第二反卷积层、第二批归一化层、第二激活层和第二全连接层后,得到指定类别的图像,其中第一损失函数如下:
第二损失函数如下:
其中,DT、DC分别表示判别器网络的图像真假部分和图像类别部分;
步骤4、训练生成对抗网络,训练步骤包括:更新判别网络D的参数,更新生成网络G的参数;
步骤5、生成最优化结果:采用遍历方式,将所有焊接参数和焊具结构输入到训练好的生成器网络G中,得到指定类别的图像;
步骤6、可选地,将模拟试验参数输入生成器网络G中,输出结果与实际结果对照以验证生成器网络G的可靠性。
本发明有益效果:
本发明提出的一种基于数值模拟和深度学习的搅拌摩擦焊缝成形预测优化方法是一种经济可靠的普适性方法,有效避免了试验手段优化中所消耗的大量时间,有效降低了材料消耗和时间成本,并为焊具设计提供了理论基础;同时,本发明提出的一种基于数值模拟和深度学习的搅拌摩擦焊缝成形预测优化方法直接获取最优化工艺参数和焊具结构,相比试验手段显著提高了预测的精确度。
附图说明
图1是本发明所提供的搅拌摩擦焊缝成形预测优化方法流程图;
图2是本发明中的数值模拟流域模型图;
图3是本发明中的生成对抗网络架构图。
具体实施方式
下面结合具体实施例对本发明做进一步说明,但本发明不受实施例的限制。
实施例1:
一种基于数值模拟和深度学习的搅拌摩擦焊缝成形预测优化方法,如图1所示,所述搅拌摩擦焊缝成形预测优化方法包括:
步骤一、可选地,根据实际生产需求,初步设置至少三次针对搅拌摩擦焊缝成形的模拟试验,作为数值模拟模型的数据基础和生成对抗网络(Generative AdversarialNetwork)深度学习模型的测试集,然后生成对抗网络深度学习模型;
步骤二、根据给定的实际工况及材料物性参数建立用于数值模拟模型,计算焊接过程中材料流动场和温度场分布情况;
步骤三、根据步骤二获得的材料流动场和温度场分布情况,结合单向流固耦合模型,计算不同参数下搅拌摩擦焊具断裂失效情况;然后,结合离散相模型,计算不同参数下焊缝成形质量及缺陷分布情况;
步骤四、根据上述计算所得焊缝成形质量及缺陷分布情况和焊具断裂失效情况,利用生成对抗网络深度学习模型遍历所有工艺参数和焊具结构的焊缝成形结果,从而获得在保证焊具可靠工作前提下焊缝成形最优化结果。
其中,所述模拟实验所采集的参数包括:焊具结构,焊接参数,焊具断裂失效状况,焊缝成形质量结果,且其中应包括至少一组用以完成焊接的参数。
本实施例所述搅拌摩擦焊具断裂失效情况的确定是采用流体动力学方法,如图2所示,将被焊材料视为非牛顿、不可压缩的连续流体,并采用有限体积法进行计算,计算软件可选地为Fluent、STAR-CCM、OpenFOAM等流体动力学仿真软件;所述搅拌摩擦焊具断裂失效情况的数值计算的步骤如下:
第一步、建立被焊材料流域模型:建立符合实际生产需求结构的被焊材料流域模型,并利用滑移网格将所述焊接材料流域分割为圆柱状的滑移网格部分和除滑移网格部分外的固定部分;所述滑移网格部分旋转速率等同于焊具旋转速率;
第二步、设置被焊材料物性参数,所述被焊材料物性参数包括被焊材料的密度、流动应力、导热率和定压比热容;所述流动应力由Sellars-Tegart材料本构模型计算,所述Sellars-Tegart材料本构模型为:
其中,σe为材料流动应力,T为材料当前温度,Tmelt为材料熔点(固相线),α、A和n为材料特征常数,σmelt为熔化态材料流动应力,Z为Zener-Hollomon系数,公式为:
其中,为材料当前应变速率,Q为材料活化能,R为气体常数,e为自然指数。
第三步、确定近焊具壁面材料流动准则:将所述近焊具壁面材料流动准则分解为两部分:一部分是竖直速度vvertical,一部分是水平速度vhorizon;然后将所述水平速度vhorizon分为剪切速度vshear和法向速度vnormal;所述竖直速度vvertical、水平速度vhorizon、剪切速度vshear和法向速度vnormal的计算公式为:
其中,ijk代指焊具上各微元体;ω为焊具旋转速度;r表示微元体距焊具轴线距离;ξ表示轴向方位角;θ表示切向方位角;η为滑移系数,即壁面材料流动速率与焊具表面线速度的比值,计算方式如下:
其中,ijk代指焊具上各微元体;η0为滑移归一化常数;e表示自然指数;ω表示角速度;r表示微元体距焊具轴线距离;ω0表示角速度归一化常数;rs表示参考半径。
第四步、设定产热模型:将所述产热模型分解为两部分,所述两部分包括:摩擦热和变形热,其中,摩擦热表达粘着摩擦过程中伴随的能量转换和耗散,变形热为材料层间剪切产生;第四步所述产热模型表示为:
其中,ijk代指焊具上各微元体;F1、F2分别表示为摩擦热输入和形变热输入;β1和β2为经验参数,取值为0.1~1.0;r表示微元体距焊具轴线距离;ω表示角速度;η表示滑移系数;σe表示微元体等效应力;μ表示动摩擦因数;P表示焊接下压力;ξ表示轴向方位角;θ表示切向方位角;表示等效应变速率。
第五步、利用牛顿换热原理和蒂芬-玻尔兹曼定律确定所述被焊材料流域的流域边界换热条件,本步骤中,流域边界换热条件可由牛顿换热公式和斯蒂芬-玻尔兹曼定律给出,速度边界条件可视为完全自由流动或完全约束流动;
第六步、从所述流域的入口添加离散相模型,所述离散相模型的输入入口同于流域速度入口,离散相为无质量粒子;
第七步、利用单向流固耦合模型计算焊具断裂失效状况,本实施方式所示单向流固耦合模型基于流体动力学软件和有限元分析软件间的单向压力和温度传输,传输界面间采用相同网格结构,有限元软件可选地为ABAQUS或ANSYS等;
第八步、提交运算并输出结果,数值计算完成后,即可获取不同焊接参数和焊具结构下的焊缝成形图像和焊具断裂失效状况。
所述生成对抗网络深度学习模型的过程如下:
步骤1、建立图像标签库:将数值模拟结果所得焊缝成形图像与焊接断裂失效状况进行分别和标注,每张图像具有唯一与之对应的焊接参数和焊具结构标签k;
步骤2、构建生成器网络G:在所述生成器网络G中,将参数标签k与随机噪声组合为新的向量输入,将所述生成器网络G生成的数据作为判别器网络D的输入,生成器网络G包括第一输入层、第一全连接层、第一反卷积层、第一批归一化层和第一激活层,以随机的噪声信号和类别标签组合而成的向量作为所示第一输入层的输入,所述输入经过所述第一全连接层、第一反卷积层、第一批归一化层和第一激活层生成数据,其中反卷积次数可选地为4~6次,以生成不低于三通道128*128分辨率的焊缝成形特征图,并将其作为判别器网络D的输入;
步骤3、构建判别器网络D:所述判别器网络D的损失函数由用于判别图像真假的第一损失函数和用于判别图像类别的第二损失函数组成,叛变期网络D包括第二输入层、第二反卷积层、第二批归一化层、第二激活层和第二全连接层,分别以步骤2中所示生成器网络K生成的数据和数值模拟所得焊缝成形图像作为所示第二输入层的输入,所述第二全连接层由用于判断图像真假的全连接层和判断图像类别的全连接层组成,所述第二输入层的输入经过所述第二反卷积层、第二批归一化层、第二激活层和第二全连接层后,得到指定类别的图像,其中第一损失函数如下:
第二损失函数如下:
其中,DT、DC分别表示判别器网络的图像真假部分和图像类别部分;
步骤4、训练生成对抗网络,训练步骤包括:更新判别网络D的参数,更新生成网络G的参数;
步骤5、生成最优化结果:采用遍历方式,将所有焊接参数和焊具结构输入到训练好的生成器网络G中,得到指定类别的图像;
步骤6、可选地,将模拟试验参数输入生成器网络G中,输出结果与实际结果对照以验证生成器网络G的可靠性。
本发明提出的一种基于数值模拟和深度学习的搅拌摩擦焊缝成形预测优化方法是一种经济可靠的普适性方法,有效避免了试验手段优化中所消耗的大量时间,有效降低了材料消耗和时间成本,并为焊具设计提供了理论基础;同时,本发明提出的一种基于数值模拟和深度学习的搅拌摩擦焊缝成形预测优化方法直接获取最优化工艺参数和焊具结构,相比试验手段显著提高了预测的精确度。
虽然本发明已以较佳的实施例公开如上,但其并非用以限定本发明,任何熟悉此技术的人,在不脱离本发明的精神和范围内,都可以做各种改动和修饰,因此本发明的保护范围应该以权利要求书所界定的为准。
Claims (7)
1.一种基于数值模拟和深度学习的搅拌摩擦焊缝成形预测优化方法,其特征在于,所述搅拌摩擦焊缝成形预测优化方法包括:
步骤一、进行至少三次针对搅拌摩擦焊缝成形的模拟试验,作为数值模拟模型的数据基础和生成对抗网络深度学习模型的测试集,然后生成对抗网络深度学习模型;
步骤二、根据给定的实际工况及材料物性参数建立用于数值模拟模型,计算焊接过程中材料流动场和温度场分布情况;
步骤三、根据步骤二获得的材料流动场和温度场分布情况,结合单向流固耦合模型,计算不同参数下搅拌摩擦焊具断裂失效情况;然后,结合离散相模型,计算不同参数下焊缝成形质量及缺陷分布情况;
步骤四、根据上述计算所得焊缝成形质量及缺陷分布情况和焊具断裂失效情况,利用生成对抗网络深度学习模型遍历所有工艺参数和焊具结构的焊缝成形结果,从而获得在保证焊具可靠工作前提下焊缝成形最优化结果。
2.根据权利要求1所述搅拌摩擦焊缝成形预测优化方法,其特征在于,所述模拟实验所采集的参数包括:焊具结构,焊接参数,焊具断裂失效状况,焊缝成形质量结果,且其中应包括至少一组用以完成焊接的参数。
3.根据权利要求1所述搅拌摩擦焊缝成形预测优化方法,其特征在于,所述搅拌摩擦焊具断裂失效情况的确定是采用流体动力学方法,将被焊材料视为非牛顿、不可压缩的连续流体,并采用有限体积法进行计算;所述搅拌摩擦焊具断裂失效情况的数值计算的步骤如下:
第一步、建立被焊材料流域模型:建立被焊材料流域模型,并利用滑移网格技术将所述焊接材料流域分割为圆柱状的滑移网格部分和除滑移网格部分外的固定部分;所述滑移网格部分旋转速率等同于焊具旋转速率;
第二步、设置被焊材料物性参数,所述被焊材料物性参数包括被焊材料的密度、流动应力、导热率和定压比热容;
第三步、确定近焊具壁面材料流动准则:将所述近焊具壁面材料流动准则分解为两部分:一部分是竖直速度vvertical,一部分是水平速度vhorizon;然后将所述水平速度vhorizon分为剪切速度vshear和法向速度vnormal;
第四步、设定产热模型:将所述产热模型分解为两部分,所述两部分包括:摩擦热和变形热;
第五步、利用牛顿换热原理和蒂芬-玻尔兹曼定律确定所述被焊材料流域的流域边界换热条件,
第六步、从所述流域的入口添加离散相模型,所述离散相模型的输入入口同于流域速度入口,离散相为无质量粒子;
第七步、利用单向流固耦合模型计算焊具断裂失效状况
第八步、提交运算并输出结果,数值计算完成后,即可获取不同焊接参数和焊具结构下的焊缝成形图像和焊具断裂失效状况。
4.根据权利要求3所述搅拌摩擦焊缝成形预测优化方法,其特征在于,第二步所述流动应力由Sellars-Tegart材料本构模型计算,所述Sellars-Tegart材料本构模型为:
其中,σe为材料流动应力,T为材料当前温度,Tmelt为材料熔点(固相线),α、A和n为材料特征常数,σmelt为熔化态材料流动应力,Z为Zener-Hollomon系数,公式为:
其中,为材料当前应变速率,Q为材料活化能,R为气体常数,e为自然指数。
5.根据权利要求3所述搅拌摩擦焊缝成形预测优化方法,其特征在于,第三步所述竖直速度vvertical、水平速度vhorizon、剪切速度vshear和法向速度vnormal的计算公式为:
其中,ijk代指焊具上各微元体;ω为焊具旋转速度;r表示微元体距焊具轴线距离;ξ表示轴向方位角;θ表示切向方位角;η为滑移系数,即壁面材料流动速率与焊具表面线速度的比值,计算方式如下:
其中,ijk代指焊具上各微元体;η0为滑移归一化常数;e表示自然指数;ω表示角速度;r表示微元体距焊具轴线距离;ω0表示角速度归一化常数;rs表示参考半径。
6.根据权利要求3所述搅拌摩擦焊缝成形预测优化方法,其特征在于,第四步所述产热模型表示为:
其中,ijk代指焊具上各微元体;F1、F2分别表示为摩擦热输入和形变热输入;β1和β2为经验参数,取值为0.1~1.0;r表示微元体距焊具轴线距离;ω表示角速度;η表示滑移系数;σe表示微元体等效应力;μ表示动摩擦因数;P表示焊接下压力;ξ表示轴向方位角;θ表示切向方位角;表示等效应变速率。
7.根据权利要求1所述搅拌摩擦焊缝成形预测优化方法,其特征在于,所述生成对抗网络深度学习模型的过程如下:
步骤1、建立图像标签库:将数值模拟结果所得焊缝成形图像与焊接断裂失效状况进行分别和标注,每张图像具有唯一与之对应的焊接参数和焊具结构标签k;
步骤2、构建生成器网络G:在所述生成器网络G中,将参数标签k与随机噪声组合为新的向量输入,将所述生成器网络G生成的数据作为判别器网络D的输入,生成器网络G包括第一输入层、第一全连接层、第一反卷积层、第一批归一化层和第一激活层,以随机的噪声信号和类别标签组合而成的向量作为所示第一输入层的输入,所述输入经过所述第一全连接层、第一反卷积层、第一批归一化层和第一激活层生成数据,其中反卷积次数可选地为4~6次,以生成不低于三通道128*128分辨率的焊缝成形特征图,并将其作为判别器网络D的输入;
步骤3、构建判别器网络D:所述判别器网络D的损失函数由用于判别图像真假的第一损失函数和用于判别图像类别的第二损失函数组成,判别器网络D包括第二输入层、第二反卷积层、第二批归一化层、第二激活层和第二全连接层,分别以步骤2中所示生成器网络K生成的数据和数值模拟所得焊缝成形图像作为所示第二输入层的输入,所述第二全连接层由用于判断图像真假的全连接层和判断图像类别的全连接层组成,所述第二输入层的输入经过所述第二反卷积层、第二批归一化层、第二激活层和第二全连接层后,得到指定类别的图像,其中第一损失函数如下:
第二损失函数如下:
其中,DT、DC分别表示判别器网络的图像真假部分和图像类别部分;
步骤4、训练生成对抗网络,训练步骤包括:更新判别网络D的参数,更新生成网络G的参数;
步骤5、生成最优化结果:采用遍历方式,将所有焊接参数和焊具结构输入到训练好的生成器网络G中,得到指定类别的图像;
步骤6、将模拟试验参数输入生成器网络G中,输出结果与实际结果对照以验证生成器网络G的可靠性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910238706.4A CN110008554B (zh) | 2019-03-27 | 2019-03-27 | 一种基于数值模拟和深度学习的搅拌摩擦焊缝成形预测的工艺参数、焊具结构优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910238706.4A CN110008554B (zh) | 2019-03-27 | 2019-03-27 | 一种基于数值模拟和深度学习的搅拌摩擦焊缝成形预测的工艺参数、焊具结构优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110008554A true CN110008554A (zh) | 2019-07-12 |
CN110008554B CN110008554B (zh) | 2022-10-18 |
Family
ID=67168420
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910238706.4A Active CN110008554B (zh) | 2019-03-27 | 2019-03-27 | 一种基于数值模拟和深度学习的搅拌摩擦焊缝成形预测的工艺参数、焊具结构优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110008554B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110405388A (zh) * | 2019-08-05 | 2019-11-05 | 蕴硕物联技术(上海)有限公司 | 预测焊接质量的方法、装置及电子设备 |
CN111177976A (zh) * | 2019-12-25 | 2020-05-19 | 广东省焊接技术研究所(广东省中乌研究院) | 一种基于深度学习的弧焊焊缝成形精准预测方法 |
CN111767666A (zh) * | 2020-06-24 | 2020-10-13 | 中国第一汽车股份有限公司 | 一种汽车零部件激光焊连接的cae仿真模拟方法 |
CN112101432A (zh) * | 2020-09-04 | 2020-12-18 | 西北工业大学 | 一种基于深度学习的材料显微图像与性能双向预测方法 |
CN112548308A (zh) * | 2020-11-12 | 2021-03-26 | 广州瑞松智能科技股份有限公司 | 一种铝与钢搅拌摩擦焊焊接变形预测方法、电子设备、存储介质 |
CN113579462A (zh) * | 2021-08-26 | 2021-11-02 | 天津大学 | 一种双侧同时搅拌摩擦焊工艺参数优化方法 |
CN114367734A (zh) * | 2022-02-11 | 2022-04-19 | 黄山学院 | 一种自适应板材表层搅拌摩擦加工改性的优化工艺方法 |
CN114519312A (zh) * | 2022-02-24 | 2022-05-20 | 中车青岛四方机车车辆股份有限公司 | 焊缝气孔产生倾向的判别方法、装置、设备及产品 |
CN116140786A (zh) * | 2023-03-06 | 2023-05-23 | 四川艾庞机械科技有限公司 | 搅拌摩擦焊接方法及其系统 |
CN116680944A (zh) * | 2023-05-05 | 2023-09-01 | 哈尔滨工业大学 | 基于数值模拟及深度学习的电弧增材构件工艺预测方法 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101559511A (zh) * | 2009-05-22 | 2009-10-21 | 清华大学 | 一种以温度为控制变量的焊接数值模拟计算方法 |
CN101966621A (zh) * | 2010-11-04 | 2011-02-09 | 哈尔滨工业大学 | 一种电弧预热搅拌摩擦塞补焊的方法 |
US20110222321A1 (en) * | 2010-03-10 | 2011-09-15 | Kabushiki Kaisha Toshiba | Power converter and method for controlling the same |
US20150130464A1 (en) * | 2012-03-12 | 2015-05-14 | Koninklijke Philips N.V. | Power converter for powering an mri gradient coil and method of operating a power converter |
CN104794266A (zh) * | 2015-04-03 | 2015-07-22 | 上海交通大学 | 基于摩擦功和数值模拟的成形界面损伤预测方法 |
CN106202803A (zh) * | 2016-07-25 | 2016-12-07 | 四川大学 | 摩擦焊工艺热力流微观组织多物理场数值计算方法 |
CN107220600A (zh) * | 2017-05-17 | 2017-09-29 | 清华大学深圳研究生院 | 一种基于深度学习的图片生成方法及生成对抗网络 |
CN107442964A (zh) * | 2017-09-21 | 2017-12-08 | 中国科学院金属研究所 | 一种预测铝合金搅拌摩擦焊接中未焊透与弱连接缺陷的方法 |
CN108027911A (zh) * | 2015-07-29 | 2018-05-11 | 伊利诺斯工具制品有限公司 | 促进焊接即服务软件的系统及方法 |
CN108920792A (zh) * | 2018-06-21 | 2018-11-30 | 北京工业大学 | 一种基于小裂纹扩展的搅拌摩擦焊构件疲劳寿命预测方法 |
CN109345507A (zh) * | 2018-08-24 | 2019-02-15 | 河海大学 | 一种基于迁移学习的大坝图像裂缝检测方法 |
CN109522857A (zh) * | 2018-11-26 | 2019-03-26 | 山东大学 | 一种基于生成式对抗网络模型的人数估计方法 |
-
2019
- 2019-03-27 CN CN201910238706.4A patent/CN110008554B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101559511A (zh) * | 2009-05-22 | 2009-10-21 | 清华大学 | 一种以温度为控制变量的焊接数值模拟计算方法 |
US20110222321A1 (en) * | 2010-03-10 | 2011-09-15 | Kabushiki Kaisha Toshiba | Power converter and method for controlling the same |
CN101966621A (zh) * | 2010-11-04 | 2011-02-09 | 哈尔滨工业大学 | 一种电弧预热搅拌摩擦塞补焊的方法 |
US20150130464A1 (en) * | 2012-03-12 | 2015-05-14 | Koninklijke Philips N.V. | Power converter for powering an mri gradient coil and method of operating a power converter |
CN104794266A (zh) * | 2015-04-03 | 2015-07-22 | 上海交通大学 | 基于摩擦功和数值模拟的成形界面损伤预测方法 |
CN108027911A (zh) * | 2015-07-29 | 2018-05-11 | 伊利诺斯工具制品有限公司 | 促进焊接即服务软件的系统及方法 |
CN106202803A (zh) * | 2016-07-25 | 2016-12-07 | 四川大学 | 摩擦焊工艺热力流微观组织多物理场数值计算方法 |
CN107220600A (zh) * | 2017-05-17 | 2017-09-29 | 清华大学深圳研究生院 | 一种基于深度学习的图片生成方法及生成对抗网络 |
CN107442964A (zh) * | 2017-09-21 | 2017-12-08 | 中国科学院金属研究所 | 一种预测铝合金搅拌摩擦焊接中未焊透与弱连接缺陷的方法 |
CN108920792A (zh) * | 2018-06-21 | 2018-11-30 | 北京工业大学 | 一种基于小裂纹扩展的搅拌摩擦焊构件疲劳寿命预测方法 |
CN109345507A (zh) * | 2018-08-24 | 2019-02-15 | 河海大学 | 一种基于迁移学习的大坝图像裂缝检测方法 |
CN109522857A (zh) * | 2018-11-26 | 2019-03-26 | 山东大学 | 一种基于生成式对抗网络模型的人数估计方法 |
Non-Patent Citations (5)
Title |
---|
HAITAO LUO 等: "Static Simulation Analysis on Friction Stir Welding (FSW) Robot", 《THE 6TH ANNUAL IEEE INTERNATIONAL CONFERENCE ON CYBER TECHNOLOGY IN AUTOMATION, CONTROL AND INTELLIGENT SYSTEMS》 * |
YONGXIAN HUANG 等: "Numerical design of high depth-to-width ratio friction stir welding", 《JOURNAL OF MATERIALS PROCESSING TECH》 * |
史雄华 等: "船体结构焊接变形的预测与控制研究进展", 《造船技术》 * |
张利国 等: "基于支持向量机的搅拌摩擦焊接头力学性能预测", 《沈阳航空航天大学学报》 * |
谢聿铭: "大深宽比搅拌摩擦焊数字化设计与过程仿真", 《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅰ辑》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110405388B (zh) * | 2019-08-05 | 2022-02-11 | 蕴硕物联技术(上海)有限公司 | 预测焊接质量的方法、装置及电子设备 |
CN110405388A (zh) * | 2019-08-05 | 2019-11-05 | 蕴硕物联技术(上海)有限公司 | 预测焊接质量的方法、装置及电子设备 |
CN111177976B (zh) * | 2019-12-25 | 2024-04-26 | 广东省科学院中乌焊接研究所 | 一种基于深度学习的弧焊焊缝成形精准预测方法 |
CN111177976A (zh) * | 2019-12-25 | 2020-05-19 | 广东省焊接技术研究所(广东省中乌研究院) | 一种基于深度学习的弧焊焊缝成形精准预测方法 |
CN111767666A (zh) * | 2020-06-24 | 2020-10-13 | 中国第一汽车股份有限公司 | 一种汽车零部件激光焊连接的cae仿真模拟方法 |
CN112101432A (zh) * | 2020-09-04 | 2020-12-18 | 西北工业大学 | 一种基于深度学习的材料显微图像与性能双向预测方法 |
CN112101432B (zh) * | 2020-09-04 | 2022-06-07 | 西北工业大学 | 一种基于深度学习的材料显微图像与性能双向预测方法 |
CN112548308A (zh) * | 2020-11-12 | 2021-03-26 | 广州瑞松智能科技股份有限公司 | 一种铝与钢搅拌摩擦焊焊接变形预测方法、电子设备、存储介质 |
CN113579462A (zh) * | 2021-08-26 | 2021-11-02 | 天津大学 | 一种双侧同时搅拌摩擦焊工艺参数优化方法 |
CN114367734A (zh) * | 2022-02-11 | 2022-04-19 | 黄山学院 | 一种自适应板材表层搅拌摩擦加工改性的优化工艺方法 |
CN114367734B (zh) * | 2022-02-11 | 2022-10-28 | 黄山学院 | 一种自适应板材表层搅拌摩擦加工改性的优化工艺方法 |
CN114519312A (zh) * | 2022-02-24 | 2022-05-20 | 中车青岛四方机车车辆股份有限公司 | 焊缝气孔产生倾向的判别方法、装置、设备及产品 |
CN116140786B (zh) * | 2023-03-06 | 2023-07-14 | 四川艾庞机械科技有限公司 | 搅拌摩擦焊接方法及其系统 |
CN116140786A (zh) * | 2023-03-06 | 2023-05-23 | 四川艾庞机械科技有限公司 | 搅拌摩擦焊接方法及其系统 |
CN116680944A (zh) * | 2023-05-05 | 2023-09-01 | 哈尔滨工业大学 | 基于数值模拟及深度学习的电弧增材构件工艺预测方法 |
CN116680944B (zh) * | 2023-05-05 | 2023-12-22 | 哈尔滨工业大学 | 基于数值模拟及深度学习的电弧增材构件工艺预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110008554B (zh) | 2022-10-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110008554A (zh) | 一种基于数值模拟和深度学习的搅拌摩擦焊缝成形预测优化方法 | |
Yan et al. | An integrated process–structure–property modeling framework for additive manufacturing | |
Stubblefield et al. | A meshfree computational framework for the numerical simulation of the solid-state additive manufacturing process, additive friction stir-deposition (AFS-D) | |
Dialami et al. | A fast and accurate two-stage strategy to evaluate the effect of the pin tool profile on metal flow, torque and forces in friction stir welding | |
Li et al. | Modeling of continuous drive friction welding of mild steel | |
Zhang et al. | A deformable surface model for real-time water drop animation | |
Dialami et al. | Material flow visualization in friction stir welding via particle tracing | |
CN106021659A (zh) | 一种冲蚀-二氧化碳腐蚀耦合作用下的天然气注采井管柱腐蚀速率的确定方法 | |
Wu et al. | CFD simulations of GMA welding of horizontal fillet joints based on coordinate rotation of arc models | |
Li et al. | Metal flow of weld pool and keyhole evolution in gas focusing plasma arc welding | |
Shams et al. | Large eddy simulation of a randomly stacked nuclear pebble bed | |
Dialami et al. | Numerical simulation and visualization of material flow in friction stir welding via particle tracing | |
Bai et al. | Numerical simulation of the Richtmyer-Meshkov instability in initially nonuniform flows and mixing with reshock | |
Li et al. | A level set based immersed boundary method for simulation of non-isothermal viscoelastic melt filling process | |
Kumar et al. | Model designed to acquire an optimized performance implementing l27 orthogonal array for the prandtl fluid flow maneuvering grey relational theory | |
Jia et al. | Numerical and experimental study on the thermal process, material flow and welding defects during high-speed friction stir welding | |
Li et al. | Quantitative analysis and experimental study of the influence of process parameters on the evolution of laser cladding | |
Li et al. | Study on parameter optimization of laser cladding Fe60 based on GA-BP neural network | |
Lang et al. | Stability mechanism of the molten pool in variable polarity plasma arc welding of medium thickness aluminum alloy | |
Song et al. | Numerical and experimental study on the properties of laser cladding of 6061Al alloy | |
CN109215100A (zh) | 一种混合流体相变动画生成方法及装置 | |
Lu et al. | Optimization of process parameters in friction stir welding of 2219 aluminum alloy thick plate | |
CN117330470A (zh) | 一种基于单方向渗透率确定三维裂隙岩体渗透率张量的方法 | |
Chen et al. | Investigation of thermal behavior of powder stream and molten pool during laser-based directed energy deposition | |
Zhou et al. | Research and application on 6D integrated system in metro construction based on BIM |
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 |