CN117195662A - 一种大型结构件焊接变形预测方法及*** - Google Patents
一种大型结构件焊接变形预测方法及*** Download PDFInfo
- Publication number
- CN117195662A CN117195662A CN202311399013.6A CN202311399013A CN117195662A CN 117195662 A CN117195662 A CN 117195662A CN 202311399013 A CN202311399013 A CN 202311399013A CN 117195662 A CN117195662 A CN 117195662A
- Authority
- CN
- China
- Prior art keywords
- welding
- deformation
- inherent
- measuring point
- joint
- 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.)
- Pending
Links
- 238000003466 welding Methods 0.000 title claims abstract description 437
- 238000000034 method Methods 0.000 title claims abstract description 68
- 238000012360 testing method Methods 0.000 claims abstract description 114
- 238000004088 simulation Methods 0.000 claims abstract description 61
- 238000005259 measurement Methods 0.000 claims description 60
- 210000001503 joint Anatomy 0.000 claims description 40
- 238000005452 bending Methods 0.000 claims description 30
- 238000004364 calculation method Methods 0.000 claims description 20
- 239000000463 material Substances 0.000 claims description 14
- NCGICGYLBXGBGN-UHFFFAOYSA-N 3-morpholin-4-yl-1-oxa-3-azonia-2-azanidacyclopent-3-en-5-imine;hydrochloride Chemical compound Cl.[N-]1OC(=N)C=[N+]1N1CCOCC1 NCGICGYLBXGBGN-UHFFFAOYSA-N 0.000 claims description 12
- 238000006073 displacement reaction Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 6
- 239000000758 substrate Substances 0.000 claims description 5
- 239000000945 filler Substances 0.000 claims description 4
- 238000013461 design Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000000399 orthopedic effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
Landscapes
- Butt Welding And Welding Of Specific Article (AREA)
Abstract
本发明公开了一种大型结构件焊接变形预测方法及***,包括基于仿真固有变形参数Fk对固有变形参数Gk进行更新,重复对焊接接头有限元模型进行焊接变形仿真过程进行迭代,直至仿真固有变形参数Fk与固有变形参数Gk的比值在设定范围Q内,将固有变形参数Gk及相对应焊接试验条件和焊接接头形式作为固有变形数据录入焊接变形数据库中;通过对固有变形数据进行修正,从而保证固有变形数据的采集精度;调取焊接变形数据库内的固有变形数据对结构件有限元模型进行焊接仿真,获取结构件的焊接变形云图,实现对当前大型结构件焊接进行快速和精准的变形预测。
Description
技术领域
本发明属于焊接变形预测技术领域,具体涉及大型结构件焊接变形预测方法及***。
背景技术
工程机械领域常常涉及大型结构件的焊接,这些结构件的焊缝具有数量多、尺寸大、变截面的特点,由于焊接过程的不均匀加热和冷却,不可避免的会存在残余应力与变形。焊接变形的存在会影响成形精度,继而影响结构件的局部强度和稳定性,降低产品使用性能。因此需要在焊后花费大量的人力、物力和时间对结构件进行矫形处理,甚至部分焊接结构件无法矫正到满足生产要求,大大影响产品生产效率和生产成本。如何有效的控制焊接变形成为工程师亟待考虑的问题。
随着数值模拟技术的发展,技术人员将其应用于焊接变形的仿真和预测。目前焊接变形预测方法主要有热弹塑性有限元法和固有应变法。热弹塑性有限元法具有计算结果精确,但是对于大型结构件则耗时过长,计算效率较低,不能满足产品设计周期要求。
固有应变法是基于弹性理论的有限元分析方法,其通过获取焊接接头固有应变,并将其作为初始载荷施加到整条焊缝上进行一次弹性计算,从而得到整个结构的焊接变形及残余应力,因此计算效率非常高。但是存在固有应变参数测量困难、测量精度不高的技术难题。
发明内容
本发明的目的在于提供一种大型结构件焊接变形预测方法及***,解决了当前大型结构件焊接变形预测周期长,固有应变参数测量困难、测量精度不高的技术难题。
为达到上述目的,第一方面本发明所采用的技术方案是:一种大型结构件焊接变形预测方法,包括:
根据待预测的结构件建立结构件有限元模型,调取焊接变形数据库内的固有变形数据对结构件有限元模型进行焊接仿真,获取结构件的焊接变形云图;
所述焊接变形数据库内的固有变形数据的获取过程为:
通过分别改变焊接试验条件和焊接接头形式进行焊接试验,其中,所述焊接试验中由多个焊接试板焊接形成所述焊接接头,记录各焊接试验对应的焊接变形数据;基于焊接变形数据计算焊接试板的固有变形参数Gk;
基于壳单元对焊接试验中焊接接头建立焊接接头有限元模型;根据固有变形参数Gk对焊接接头有限元模型进行焊接变形仿真,获得焊接变形仿真数据;基于焊接变形仿真数据计算焊接试板的仿真固有变形参数Fk;
若仿真固有变形参数Fk与固有变形参数Gk的比值在设定范围Q外,基于仿真固有变形参数Fk对固有变形参数Gk进行更新,重复对焊接接头有限元模型进行焊接变形仿真过程进行迭代,直至仿真固有变形参数Fk与固有变形参数Gk的比值在设定范围Q内,将固有变形参数Gk及相对应焊接试验条件和焊接接头形式作为固有变形数据录入焊接变形数据库中。
优选的,所述焊接试验条件包括焊接试板厚度、焊接方法、焊接材料、填充材料、焊接位置、焊接电流、焊接电压和焊接速度。
优选的,所述焊接接头形式包括对接接头、T形接头和搭接接头。
优选的,当焊接试验中焊接接头形式为对接接头时,基于焊接变形数据计算焊接试板的固有变形参数Gk的方法包括:
当焊接试验中焊接接头形式为焊接板HA1和焊接板HB1组成的对接接头时,所述焊接板HA1和焊接板HB1上均匀设定有呈2×3矩阵分布的6个测量点,记为测量点Pi α,i∈[1,…,12];测量点测量点/>和测量点/>设置为同一行并设置于焊接板HA1远离所述对接接头的一侧;测量点/>测量点/>和测量点/>设置为同一行并设置于焊接板HA1靠近所述对接接头的一侧;测量点/>测量点/>和测量点/>设置为同一行并设置于焊接板HB1靠近所述对接接头的一侧;测量点/>测量点/>和测量点/>设置为同一行并设置于焊接板HB1远离所述对接接头的一侧;所述焊接变形数据包括焊接前后测量点Pi α的位置坐标;
当焊接试验中焊接接头形式为对接接头时,所述固有变形参数Gk包含固有横向收缩δT,固有纵向收缩δL,固有横向弯曲θT和固有纵向弯曲θL,所述固有变形参数Gk计算公式分别为:
公式中,表示为焊接前测量点Pi α的初始横坐标,/>表示为焊接前测量点Pi α的初始纵坐标,/>表示为焊接变形后测量点Pi α的横坐标,/>表示为焊接变形后测量点Pi α的变形后的纵坐标,/>表示为焊接变形后由测量点/>测量点/>和测量点/>拟合形成的曲线对应的弦长,/>表示为焊接变形后由测量点/>测量点/>和测量点/>拟合形成的曲线对应的弦长,/>为焊接变形后任一测量点Pi α在z轴方向上的位移;Bα为对接接头宽度。
优选的,当焊接试验中焊接接头形式为T形接头时,基于焊接变形数据计算焊接试板的固有变形参数Gk的方法包括:
当焊接试验中焊接接头形式为焊接板HA2和焊接板HB2组成的T形接头时,所述焊接板HA2设定为底板,所述焊接板HB2设定为腹板,所述焊接板HB2焊接于焊接板HA2中间位置;所述T形接头上的测量点记为测量点其中,测量点/>至测量点/>以及测量点/>和测量点/>设置于所述焊接板HA2,测量点/>测量点/>测量点/>测量点/>以及测量点/>与测量点/>测量点/>测量点/>测量点/>以及测量点/>对称设置于所述焊接板HB2两侧;所述测量点/>和测量点/>设置于所述焊接板HA2和焊接板HB2焊接处;所述测量点/>至测量点/>设置于所述焊接板HB2上;
当焊接试验中焊接接头形式为T形接头时,所述固有变形参数Gk包含底板固有横向收缩δTP,底板固有纵向收缩δLP,底板固有横向弯曲θTP,底板固有纵向弯曲θLP,腹板固有横向收缩δTW,腹板固有纵向收缩δLW,腹板固有横向弯曲θTW和腹板固有纵向弯曲θLW,其计算公式分别为:
δTP=b-[(b+Δb)/2]/cos(θT1)-[(b+Δb)/2]/cos(θT2)
δLP=ΔL×B/L
θTP=sin-1(e1/d1)+sin-1(e2/d2)
δTW=H-[(H-ΔH)/cos(tan-1(c/H))]
δLW=ΔLw×H/L
θTW=sin-1(c/H)-[sin-1(e1/d1)-sin-1(e2/d2)]/2
θLP=θLW=0
公式中,表示为焊接前测量点/>的初始横坐标,/>表示为焊接前测量点/>的初始纵坐标,/>表示为焊接变形后测量点/>的横坐标,/>表示为焊接后测量点/>的变形后的纵坐标,b表示为测量点/>与/>测量点/>与/>测量点/>与/>Y方向之间距离的平均值,Δb表示为焊接后b值的变化量,L表示为测量点/>与/>测量点/>与/>X方向之间距离的平均值,ΔL表示为焊接后焊接板HA2的L的变化量,B表示为T形接头宽度,e1表示为焊接后焊接板HA2左侧在Z方向的变形量,e2表示为焊接板HA2右侧在Z方向的变形量,d1表示为测量点/>与/>测量点/>与/>Y方向之间距离的平均值,d2表示为测量点/>与/>测量点/>与/>Y方向之间距离的平均值,H表示腹板高度,ΔH表示焊接后腹板高度的收缩量,c表示焊接后腹板顶端横向位移,ΔLw表示为焊接后焊接板HB2的L的变化量。
优选的,当焊接试验中焊接接头形式为搭接接头时,基于焊接变形数据计算焊接试板的固有变形参数Gk的方法包括:
焊接试验中焊接接头形式为焊接板HA3和焊接板HB3组成的搭接接头时,所述焊接板HA3和焊接板HB3上均匀设定有呈2×3矩阵分布的6个测量点,记为测量点测量点/>测量点/>和测量点/>设置为同一行并设置于焊接板HA3远离所述搭接接头的一侧;测量点/>测量点/>和测量点/>设置为同一行并设置于焊接板HA3靠近所述搭接接头的一侧;测量点/>测量点/>和测量点/>设置为同一行并设置于焊接板HB3靠近所述搭接接头的一侧;测量点/>测量点/>和测量点/>设置为同一行并设置于焊接板HB3远离所述搭接接头的一侧;所述焊接变形数据包括焊接前后测量点/>的位置坐标;
当焊接试验中焊接接头形式为搭接接头时,所述固有变形参数Gk包含固有横向收缩δT1,固有纵向收缩δL1,固有横向弯曲θT1,固有纵向弯曲θL1,固有横向收缩δT2,固有纵向收缩δL2,固有横向弯曲θT2和固有纵向弯曲θl2,其计算公式分别为:
δT2=δT1
公式中,表示为焊接前测量点/>的初始横坐标,/>表示为焊接前测量点/>的初始纵坐标,/>表示为焊接后测量点/>的横坐标,/>表示为焊接后测量点/>的变形后的纵坐标,/>为任一测量点/>在z轴方向上的位移;Bγ为搭接接头长度,cγ表示为接头搭接宽度,tγ表示为焊接板HA3和焊接板HB3的厚度。
优选的,基于仿真固有变形参数Fk对固有变形参数Gk进行更新的方法包括:
公式中,Gk+1表示为第k+1次更新迭代后的固有变形参数。
优选的,仿真固有变形参数Fk与固有变形参数Gk比值的设定范围Q为0.9至1.1。
第二方面本发明提供了一种大型结构件焊接变形预测***,包括:
焊接变形预测模块,用于根据待预测的结构件建立结构件有限元模型,调取焊接变形数据库内的固有变形数据对结构件有限元模型进行焊接仿真,获取结构件的焊接变形云图;
试验数据获取模块,用于通过分别改变焊接试验条件和焊接接头形式进行焊接试验,其中,所述焊接试验中由多个焊接试板焊接形成所述焊接接头,记录各焊接试验对应的焊接变形数据;基于焊接变形数据计算焊接试板的固有变形参数Gk;
计算模块,用于基于壳单元对焊接试验中焊接接头建立焊接接头有限元模型;根据固有变形参数Gk对焊接接头有限元模型进行焊接变形仿真,获得焊接变形仿真数据;基于焊接变形仿真数据计算焊接试板的仿真固有变形参数Fk;
迭代模块,用于若仿真固有变形参数Fk与固有变形参数Gk的比值在设定范围Q外,基于仿真固有变形参数Fk对固有变形参数Gk进行更新,重复对焊接接头有限元模型进行焊接变形仿真过程进行迭代,直至仿真固有变形参数Fk与固有变形参数Gk的比值在设定范围Q内,将固有变形参数Gk及相对应焊接试验条件和焊接接头形式作为固有变形数据录入焊接变形数据库中。
优选的,所述焊接试验条件包括焊接试板厚度、焊接方法、焊接材料、填充材料、焊接位置、焊接电流、焊接电压和焊接速度。
优选的,所述焊接接头形式包括对接接头、T形接头和搭接接头。
优选的,基于仿真固有变形参数Fk对固有变形参数Gk进行更新的方法包括:
公式中,Gk+1表示为第k+1次更新迭代后的固有变形参数。
优选的,仿真固有变形参数Fk与固有变形参数Gk比值的设定范围Q为0.9至1.1。
第三方面本发明提供了电子设备包括存储介质和处理器;所述存储介质用于存储指令;所述处理器用于根据所述指令进行操作以执行第一方面所述的方法。
与现有技术相比,本发明所达到的有益效果:
本发明若仿真固有变形参数Fk与固有变形参数Gk的比值在设定范围Q外,基于仿真固有变形参数Fk对固有变形参数Gk进行更新,重复对焊接接头有限元模型进行焊接变形仿真过程进行迭代,直至仿真固有变形参数Fk与固有变形参数Gk的比值在设定范围Q内,将固有变形参数Gk及相对应焊接试验条件和焊接接头形式作为固有变形数据录入焊接变形数据库中,通过对固有变形数据进行修正,从而保证固有变形数据的采集精度;根据待预测的结构件建立结构件有限元模型,调取焊接变形数据库内的固有变形数据对结构件有限元模型进行焊接仿真,获取结构件的焊接变形云图,实现对当前大型结构件焊接进行快速和精准的变形预测。
附图说明
图1是实施例1提供的大型结构件焊接变形预测方法的流程图;
图2是实施例1提供的焊接接头形式的对比示意图;
图3是实施例1提供的对接接头的计算模型示意图;
图4是实施例1提供的T形接头的尺寸计算模型示意图;
图5是实施例1提供的T形接头的测试点分布示意图;
图6是实施例1提供的搭接接头的计算模型示意图;
图7是实施例1提供的固有变形数据修正的流程图;
图8是实施例2提供的大型结构件焊接变形预测***的结构图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
需要说明的是,在本发明的描述中,术语“前”、“后”、“左”、“右”、“上”、“下”、“内”、“外”等指示的方位或位置关系为基于附图中所示的方位或位置关系,仅是为了便于描述本发明而不是要求本发明必须以特定的方位构造和操作,因此不能理解为对本发明的限制。本发明描述中使用的术语“前”、“后”、“左”、“右”、“上”、“下”指的是附图中的方向,术语“内”、“外”分别指的是朝向或远离特定部件几何中心的方向。
实施例1
如图1至图7所示,一种大型结构件焊接变形预测方法,包括
利用Hypermesh软件根据待预测的结构件建立结构件有限元模型,为结构件有限元模型中每一组件分别赋予材料弹性模量、泊松比和厚度属性,基于实际情况为模型施加约束条件,调取焊接变形数据库内的固有变形数据对结构件有限元模型进行焊接仿真,获取结构件的焊接变形云图;
所述焊接变形数据库内的固有变形数据的获取过程为:
通过分别改变焊接试验条件和焊接接头形式进行焊接试验,其中,所述焊接试验中由多个焊接试板焊接形成所述焊接接头,记录各焊接试验对应的焊接变形数据;所述焊接试验条件包括焊接试板厚度、焊接方法、焊接材料、填充材料、焊接位置、焊接电流、焊接电压和焊接速度;如图2所示,所述焊接接头形式包括对接接头、T形接头和搭接接头;当焊接试验中焊接接头形式为对接接头时,基于焊接变形数据计算焊接试板的固有变形参数Gk的方法包括:
如图3所示,当焊接试验中焊接接头形式为焊接板HA1和焊接板HB1组成的对接接头时,所述焊接板HA1和焊接板HB1上均匀设定有呈2×3矩阵分布的6个测量点,记为测量点Pi α,i∈[1,…,12];测量点测量点/>和测量点/>设置为同一行并设置于焊接板HA1远离所述对接接头的一侧;测量点/>测量点/>和测量点/>设置为同一行并设置于焊接板HA1靠近所述对接接头的一侧;测量点/>测量点/>和测量点/>设置为同一行并设置于焊接板HB1靠近所述对接接头的一侧;测量点/>测量点/>和测量点/>设置为同一行并设置于焊接板HB1远离所述对接接头的一侧;所述焊接变形数据包括焊接前后测量点Pi α的位置坐标;如表1所示对接接头板材焊接试验推荐尺寸。
表1对接接头板材焊接试验推荐尺寸表
当焊接试验中焊接接头形式为对接接头时,所述固有变形参数Gk包含固有横向收缩δT,固有纵向收缩δL,固有横向弯曲θT和固有纵向弯曲θL,所述固有变形参数Gk计算公式分别为:
公式中,表示为焊接前测量点Pi α的初始横坐标,/>表示为焊接前测量点Pi α的初始纵坐标,/>表示为焊接变形后测量点Pi α的横坐标,/>表示为焊接变形后测量点Pi α的变形后的纵坐标,/>表示为焊接变形后由测量点/>测量点/>和测量点/>拟合形成的曲线对应的弦长,/>表示为焊接变形后由测量点/>测量点/>和测量点/>拟合形成的曲线对应的弦长,/>为焊接变形后任一测量点Pi α在z轴方向上的位移;Bα为对接接头宽度。
如图4和图5所示,当焊接试验中焊接接头形式为T形接头时,基于焊接变形数据计算焊接试板的固有变形参数Gk的方法包括:
当焊接试验中焊接接头形式为焊接板HA2和焊接板HB2组成的T形接头时,所述焊接板HA2设定为底板,所述焊接板HB2设定为腹板,所述焊接板HB2焊接于焊接板HA2中间位置;所述T形接头上的测量点记为测量点其中,测量点/>至测量点/>以及测量点/>和测量点/>设置于所述焊接板HA2,测量点/>测量点/>测量点/>测量点/>以及测量点/>与测量点/>测量点/>测量点/>测量点/>以及测量点/>对称设置于所述焊接板HB2两侧;所述测量点/>和测量点/>设置于所述焊接板HA2和焊接板HB2焊接处;所述测量点/>至测量点/>设置于所述焊接板HB2上;如表2所示T形接头板材焊接试验相关参数定义及推荐尺寸。
表2T形接头板材焊接试验相关参数定义及推荐尺寸
当焊接试验中焊接接头形式为T形接头时,所述固有变形参数Gk包含底板固有横向收缩δTP,底板固有纵向收缩δLP,底板固有横向弯曲θTP,底板固有纵向弯曲θLP,腹板固有横向收缩δTW,腹板固有纵向收缩δLW,腹板固有横向弯曲θTW和腹板固有纵向弯曲θLW,其计算公式分别为:
δTP=b-[(b+Δb)/2]/cos(θT1)-[(b+Δb)/2]/cos(θT2)
δLP=ΔL×B/L
θTP=sin-1(e1/d1)+sin-1(e2/d2)
δTW=H-[(H-ΔH)/cos(tan-1(c/H))]
δLW=ΔLw×H/L
θTW=sin-1(c/H)-[sin-1(e1/d1)-sin-1(e2/d2)]/2
θLP=θLW=0
公式中,表示为焊接前测量点/>的初始横坐标,/>表示为焊接前测量点/>的初始纵坐标,/>表示为焊接变形后测量点/>的横坐标,/>表示为焊接后测量点/>的变形后的纵坐标,b表示为测量点/>与/>测量点/>与/>测量点/>与/>Y方向之间距离的平均值,Δb表示为焊接后b值的变化量,L表示为测量点/>与/>测量点/>与/>X方向之间距离的平均值,ΔL表示为焊接后焊接板HA2的L的变化量,B表示为T形接头宽度,e1表示为焊接后焊接板HA2左侧在Z方向的变形量,e2表示为焊接板HA2右侧在Z方向的变形量,d1表示为测量点/>与/>测量点/>与/>Y方向之间距离的平均值,d2表示为测量点/>与/>测量点/>与/>Y方向之间距离的平均值,H表示腹板高度,ΔH表示焊接后腹板高度的收缩量,c表示焊接后腹板顶端横向位移,ΔLw表示为焊接后焊接板HB2的L的变化量。
如图6所示,当焊接试验中焊接接头形式为搭接接头时,基于焊接变形数据计算焊接试板的固有变形参数Gk的方法包括:
当焊接试验中焊接接头形式为焊接板HA3和焊接板HB3组成的搭接接头时,所述焊接板HA3和焊接板HB3上均匀设定有呈2×3矩阵分布的6个测量点,记为测量点测量点/>测量点/>和测量点/>设置为同一行并设置于焊接板HA3远离所述搭接接头的一侧;测量点/>测量点/>和测量点/>设置为同一行并设置于焊接板HA3靠近所述搭接接头的一侧;测量点/>测量点/>和测量点/>设置为同一行并设置于焊接板HB3靠近所述搭接接头的一侧;测量点/>测量点/>和测量点/>设置为同一行并设置于焊接板HB3远离所述搭接接头的一侧;所述焊接变形数据包括焊接前后测量点/>的位置坐标;如表3所示搭接接头板材焊接试验推荐尺寸。
表3搭接接头板材焊接试验推荐尺寸
OL | L | B | t | c | |
推荐 | 100t+c | 60t | 50t+c |
当焊接试验中焊接接头形式为搭接接头时,所述固有变形参数Gk包含固有横向收缩δT1,固有纵向收缩δL1,固有横向弯曲θT1,固有纵向弯曲θL1,固有横向收缩δT2,固有纵向收缩δL2,固有横向弯曲θT2和固有纵向弯曲θL2,其计算公式分别为:
δT2=δT1
/>
公式中,表示为焊接前测量点/>的初始横坐标,/>表示为焊接前测量点/>的初始纵坐标,/>表示为焊接后测量点/>的横坐标,/>表示为焊接后测量点/>的变形后的纵坐标,/>为任一测量点/>在z轴方向上的位移;Bγ为搭接接头长度,cγ表示为接头搭接宽度,tγ表示为焊接板HA3和焊接板HB3的厚度;
如图7所示,基于壳单元对焊接试验中焊接接头建立焊接接头有限元模型;根据固有变形参数Gk对焊接接头有限元模型进行焊接变形仿真,获得焊接变形仿真数据;本实施中采用焊接变形预测软件JWRIAN对焊接接头有限元模型进行焊接变形仿真;基于焊接变形仿真数据计算焊接试板的仿真固有变形参数Fk;
若仿真固有变形参数Fk与固有变形参数Gk的比值在设定范围Q外,基于仿真固有变形参数Fk对固有变形参数Gk进行更新的方法包括:
公式中,Gk+1表示为第k+1次更新迭代后的固有变形参数。
重复对焊接接头有限元模型进行焊接变形仿真过程进行迭代,直至仿真固有变形参数Fk与固有变形参数Gk的比值在设定范围Q内,将固有变形参数Gk及相对应焊接试验条件和焊接接头形式作为固有变形数据录入焊接变形数据库中;仿真固有变形参数Fk与固有变形参数Gk比值的设定范围Q为0.9至1.1;根据仿真固有变形参数Fk对固有变形参数Gk进行迭代,以消除试验测量带来的误差,实现固有变形参数Gk的高精度获取。
基于固有应变法利用壳单元对大型结构件进行焊接变形预测,与基于热弹塑性有限元的焊接变形仿真相比,预测时间由天级降到了小时级,计算效率显著提高,且与基于实体单元的固有应变法焊接变形仿真相比,模型前处理工作由复杂变为简单,降低了对仿真人员的要求;输入基于本发明获得的高精度接头固有变形参数,可实现对大型结构件焊接变形的高精高效预测。
通过该方法不断扩充焊接变形数据库,设计人员在产品设计周期可随时便捷的调用数据以获得焊接变形云图,根据焊接变形云图获得变形趋势及变形大小并确定关键焊接位置,在关键位置施加反变形或约束措施,通过焊接仿真结果指导产品设计,可替代大部分的工艺试验,从而有效缩短设计周期和节约产品研发成本。
实施例2
如图8所示,一种大型结构件焊接变形预测***,本实施例提供的***可以应用实施例1所述的方法,控制方法包括:
焊接变形预测模块,用于根据待预测的结构件建立结构件有限元模型,调取焊接变形数据库内的固有变形数据对结构件有限元模型进行焊接仿真,获取结构件的焊接变形云图;
试验数据获取模块,用于通过分别改变焊接试验条件和焊接接头形式进行焊接试验,其中,所述焊接试验中由多个焊接试板焊接形成所述焊接接头,记录各焊接试验对应的焊接变形数据;所述焊接试验条件包括焊接试板厚度、焊接方法、焊接材料、填充材料、焊接位置、焊接电流、焊接电压和焊接速度;所述焊接接头形式包括对接接头、T形接头和搭接接头;基于仿真固有变形参数Fk对固有变形参数Gk进行更新的方法包括:
公式中,Gk+1表示为第k+1次更新迭代后的固有变形参数。
计算模块,用于基于壳单元对焊接试验中焊接接头建立焊接接头有限元模型;根据固有变形参数Gk对焊接接头有限元模型进行焊接变形仿真,获得焊接变形仿真数据;基于焊接变形仿真数据计算焊接试板的仿真固有变形参数Fk;
迭代模块,用于若仿真固有变形参数Fk与固有变形参数Gk的比值在设定范围Q外,基于仿真固有变形参数Fk对固有变形参数Gk进行更新,重复对焊接接头有限元模型进行焊接变形仿真过程进行迭代,直至仿真固有变形参数Fk与固有变形参数Gk的比值在设定范围Q内,将固有变形参数Gk及相对应焊接试验条件和焊接接头形式作为固有变形数据录入焊接变形数据库中;仿真固有变形参数Fk与固有变形参数Gk比值的设定范围Q为0.9至1.1。
实施例3
电子设备包括存储介质和处理器;所述存储介质用于存储指令;所述处理器用于根据所述指令进行操作以执行实施例1所述的方法。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (14)
1.一种大型结构件焊接变形预测方法,其特征在于,包括
根据待预测的结构件建立结构件有限元模型,调取焊接变形数据库内的固有变形数据对结构件有限元模型进行焊接仿真,获取结构件的焊接变形云图;
所述焊接变形数据库内的固有变形数据的获取过程为:
通过分别改变焊接试验条件和焊接接头形式进行焊接试验,其中,所述焊接试验中由多个焊接试板焊接形成所述焊接接头,记录各焊接试验对应的焊接变形数据;基于焊接变形数据计算焊接试板的固有变形参数Gk;
基于壳单元对焊接试验中焊接接头建立焊接接头有限元模型;根据固有变形参数Gk对焊接接头有限元模型进行焊接变形仿真,获得焊接变形仿真数据;基于焊接变形仿真数据计算焊接试板的仿真固有变形参数Fk;
若仿真固有变形参数Fk与固有变形参数Gk的比值在设定范围Q外,基于仿真固有变形参数Fk对固有变形参数Gk进行更新,重复对焊接接头有限元模型进行焊接变形仿真过程进行迭代,直至仿真固有变形参数Fk与固有变形参数Gk的比值在设定范围Q内,将固有变形参数Gk及相对应焊接试验条件和焊接接头形式作为固有变形数据录入焊接变形数据库中。
2.根据权利要求1所述的大型结构件焊接变形预测方法,其特征在于,所述焊接试验条件包括焊接试板厚度、焊接方法、焊接材料、填充材料、焊接位置、焊接电流、焊接电压和焊接速度。
3.根据权利要求1所述的大型结构件焊接变形预测方法,其特征在于,所述焊接接头形式包括对接接头、T形接头和搭接接头。
4.根据权利要求3所述的大型结构件焊接变形预测方法,其特征在于,当焊接试验中焊接接头形式为对接接头时,基于焊接变形数据计算焊接试板的固有变形参数Gk的方法包括:
当焊接试验中焊接接头形式为焊接板HA1和焊接板HB1组成的对接接头时,所述焊接板HA1和焊接板HB1上均匀设定有呈2×3矩阵分布的6个测量点,记为测量点Pi α,i∈[1,…,12];测量点测量点/>和测量点/>设置为同一行并设置于焊接板HA1远离所述对接接头的一侧;测量点/>测量点/>和测量点/>设置为同一行并设置于焊接板HA1靠近所述对接接头的一侧;测量点/>测量点/>和测量点/>设置为同一行并设置于焊接板HB1靠近所述对接接头的一侧;测量点/>测量点/>和测量点/>设置为同一行并设置于焊接板HB1远离所述对接接头的一侧;所述焊接变形数据包括焊接前后测量点Pi α的位置坐标;
当焊接试验中焊接接头形式为对接接头时,所述固有变形参数Gk包含固有横向收缩δT,固有纵向收缩δL,固有横向弯曲θT和固有纵向弯曲θL,所述固有变形参数Gk计算公式分别为:
公式中,表示为焊接前测量点Pi α的初始横坐标,/>表示为焊接前测量点Pi α的初始纵坐标,/>表示为焊接变形后测量点Pi α的横坐标,/>表示为焊接变形后测量点Pi α的变形后的纵坐标,/>表示为焊接变形后由测量点/>测量点/>和测量点/>拟合形成的曲线对应的弦长,/>表示为焊接变形后由测量点/>测量点/>和测量点/>拟合形成的曲线对应的弦长,/>为焊接变形后任一测量点Pi α在z轴方向上的位移;Bα为对接接头宽度。
5.根据权利要求3所述的大型结构件焊接变形预测方法,其特征在于,当焊接试验中焊接接头形式为T形接头时,基于焊接变形数据计算焊接试板的固有变形参数Gk的方法包括:
当焊接试验中焊接接头形式为焊接板HA2和焊接板HB2组成的T形接头时,所述焊接板HA2设定为底板,所述焊接板HB2设定为腹板,所述焊接板HB2焊接于焊接板HA2中间位置;所述T形接头上的测量点记为测量点其中,测量点/>至测量点/>以及测量点/>和测量点/>设置于所述焊接板HA2,测量点/>测量点/>测量点/>测量点/>以及测量点/>与测量点/>测量点/>测量点/>测量点/>以及测量点/>对称设置于所述焊接板HB2两侧;所述测量点/>和测量点/>设置于所述焊接板HA2和焊接板HB2焊接处;所述测量点/>至测量点/>设置于所述焊接板HB2上;
当焊接试验中焊接接头形式为T形接头时,所述固有变形参数Gk包含底板固有横向收缩δTP,底板固有纵向收缩δLP,底板固有横向弯曲θTP,底板固有纵向弯曲θLP,腹板固有横向收缩δTW,腹板固有纵向收缩δLW,腹板固有横向弯曲θTW和腹板固有纵向弯曲θLW,其计算公式分别为:
δTP=b-[(b+Δb)/2]/cos(θT1)-[(b+Δb)/2]/cos(θT2)
δLP=ΔL×B/L
θTP=sin-1(e1/d1)+sin-1(e2/d2)
δTW=H-[(H-ΔH)/cos(tan-1(c/H))]
δLW=ΔLw×H/L
θTW=sin-1(c/H)-[sin-1(e1/d1)-sin-1(e2/d2)]/2
θLP=θLW=0
公式中,表示为焊接前测量点/>的初始横坐标,/>表示为焊接前测量点/>的初始纵坐标,/>表示为焊接变形后测量点/>的横坐标,/>表示为焊接后测量点/>的变形后的纵坐标,b表示为测量点/>与/>测量点/>与/>测量点/>与/>Y方向之间距离的平均值,Δb表示为焊接后b值的变化量,L表示为测量点/>与/>测量点/>与/>X方向之间距离的平均值,ΔL表示为焊接后焊接板HA2的L的变化量,B表示为T形接头宽度,e1表示为焊接后焊接板HA2左侧在Z方向的变形量,e2表示为焊接板HA2右侧在Z方向的变形量,d1表示为测量点/>与/>测量点/>与/>Y方向之间距离的平均值,d2表示为测量点/>与测量点/>与/>Y方向之间距离的平均值,H表示腹板高度,ΔH表示焊接后腹板高度的收缩量,c表示焊接后腹板顶端横向位移,ΔLw表示为焊接后焊接板HB2的L的变化量。
6.根据权利要求3所述的大型结构件焊接变形预测方法,其特征在于,当焊接试验中焊接接头形式为搭接接头时,基于焊接变形数据计算焊接试板的固有变形参数Gk的方法包括:
当焊接试验中焊接接头形式为焊接板HA3和焊接板HB3组成的搭接接头时,所述焊接板HA3和焊接板HB3上均匀设定有呈2×3矩阵分布的6个测量点,记为测量点测量点/>测量点/>和测量点/>设置为同一行并设置于焊接板HA3远离所述搭接接头的一侧;测量点/>测量点/>和测量点/>设置为同一行并设置于焊接板HA3靠近所述搭接接头的一侧;测量点/>测量点/>和测量点/>设置为同一行并设置于焊接板HB3靠近所述搭接接头的一侧;测量点/>测量点/>和测量点/>设置为同一行并设置于焊接板HB3远离所述搭接接头的一侧;所述焊接变形数据包括焊接前后测量点/>的位置坐标;
当焊接试验中焊接接头形式为搭接接头时,所述固有变形参数Gk包含固有横向收缩δT1,固有纵向收缩δL1,固有横向弯曲θT1,固有纵向弯曲θL1,固有横向收缩δT2,固有纵向收缩δL2,固有横向弯曲θT2和固有纵向弯曲θL2,其计算公式分别为:
δT2=δT1
公式中,表示为焊接前测量点/>的初始横坐标,/>表示为焊接前测量点/>的初始纵坐标,/>表示为焊接后测量点/>的横坐标,/>表示为焊接后测量点/>的变形后的纵坐标,/>为任一测量点/>在z轴方向上的位移;Bγ为搭接接头长度,cγ表示为接头搭接宽度,tγ表示为焊接板HA3和焊接板HB3的厚度。
7.根据权利要求1所述的大型结构件焊接变形预测方法,其特征在于,基于仿真固有变形参数Fk对固有变形参数Gk进行更新的方法包括:
公式中,Gk+1表示为第k+1次更新迭代后的固有变形参数。
8.根据权利要求1所述的大型结构件焊接变形预测方法,其特征在于,仿真固有变形参数Fk与固有变形参数Gk比值的设定范围Q为0.9至1.1。
9.一种大型结构件焊接变形预测***,其特征在于,包括:
焊接变形预测模块,用于根据待预测的结构件建立结构件有限元模型,调取焊接变形数据库内的固有变形数据对结构件有限元模型进行焊接仿真,获取结构件的焊接变形云图;
试验数据获取模块,用于通过分别改变焊接试验条件和焊接接头形式进行焊接试验,其中,所述焊接试验中由多个焊接试板焊接形成所述焊接接头,记录各焊接试验对应的焊接变形数据;基于焊接变形数据计算焊接试板的固有变形参数Gk;
计算模块,用于基于壳单元对焊接试验中焊接接头建立焊接接头有限元模型;根据固有变形参数Gk对焊接接头有限元模型进行焊接变形仿真,获得焊接变形仿真数据;基于焊接变形仿真数据计算焊接试板的仿真固有变形参数Fk;
迭代模块,用于若仿真固有变形参数Fk与固有变形参数Gk的比值在设定范围Q外,基于仿真固有变形参数Fk对固有变形参数Gk进行更新,重复对焊接接头有限元模型进行焊接变形仿真过程进行迭代,直至仿真固有变形参数Fk与固有变形参数Gk的比值在设定范围Q内,将固有变形参数Gk及相对应焊接试验条件和焊接接头形式作为固有变形数据录入焊接变形数据库中。
10.根据权利要求9所述的大型结构件焊接变形预测***,其特征在于,所述焊接试验条件包括焊接试板厚度、焊接方法、焊接材料、填充材料、焊接位置、焊接电流、焊接电压和焊接速度。
11.根据权利要求9所述的大型结构件焊接变形预测***,其特征在于,所述焊接接头形式包括对接接头、T形接头和搭接接头。
12.根据权利要求9所述的大型结构件焊接变形预测***,其特征在于,基于仿真固有变形参数Fk对固有变形参数Gk进行更新的方法包括:
公式中,Gk+1表示为第k+1次更新迭代后的固有变形参数。
13.根据权利要求9所述的大型结构件焊接变形预测***,其特征在于,仿真固有变形参数Fk与固有变形参数Gk比值的设定范围Q为0.9至1.1。
14.电子设备包括存储介质和处理器;所述存储介质用于存储指令;所述处理器用于根据所述指令进行操作以执行权利要求1至8任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311399013.6A CN117195662A (zh) | 2023-10-26 | 2023-10-26 | 一种大型结构件焊接变形预测方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311399013.6A CN117195662A (zh) | 2023-10-26 | 2023-10-26 | 一种大型结构件焊接变形预测方法及*** |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117195662A true CN117195662A (zh) | 2023-12-08 |
Family
ID=88988957
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311399013.6A Pending CN117195662A (zh) | 2023-10-26 | 2023-10-26 | 一种大型结构件焊接变形预测方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117195662A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117436321A (zh) * | 2023-12-21 | 2024-01-23 | 中钧科技(深圳)有限公司 | 基于有限元仿真法的工件焊接分析方法及*** |
-
2023
- 2023-10-26 CN CN202311399013.6A patent/CN117195662A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117436321A (zh) * | 2023-12-21 | 2024-01-23 | 中钧科技(深圳)有限公司 | 基于有限元仿真法的工件焊接分析方法及*** |
CN117436321B (zh) * | 2023-12-21 | 2024-04-16 | 中钧科技(深圳)有限公司 | 基于有限元仿真法的工件焊接分析方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106041344B (zh) | 一种用于航空发动机支架的激光喷丸校形质量控制方法 | |
US6560498B1 (en) | Formation method and device for curved plates | |
CN117195662A (zh) | 一种大型结构件焊接变形预测方法及*** | |
CN106529055A (zh) | 一种基于应变模态振型相关性的模型修正方法 | |
KR101284999B1 (ko) | 면내 압축 강도 평가 장치 및 방법 | |
JPWO2005093612A1 (ja) | 溶接変形算出方法、溶接変形算出装置、及びコンピュータプログラム | |
CN109359336A (zh) | 一种基于多目标优化的绑扎桥相似畸变模型构建方法 | |
CN108009327A (zh) | 一种基于钢构件变形分析的虚拟预拼装误差判定方法 | |
JP2020504859A (ja) | しわを抑制する空間薄膜構造治具形状を最適化する設計方法 | |
CN110516365A (zh) | 一种测试螺栓连接刚度的方法 | |
CN114942000B (zh) | 一种直升机尾传动轴检测校准方法 | |
CN113523656A (zh) | 一种结构件焊接反变形设计方法 | |
CN108804725A (zh) | 一种焊接反变形量的计算方法 | |
CN116579207A (zh) | 焊点焊后残余应力预测与优化的变可信度模型确定方法 | |
Zongo et al. | Geometric deviations of laser powder bed–fused AlSi10Mg components: numerical predictions versus experimental measurements | |
CN115544842A (zh) | 一种焊接变形的模拟方法、装置、设备和存储介质 | |
CN107679315A (zh) | 车身板件焊接变形的几何补偿方法及*** | |
CN115795770A (zh) | 一种车身焊接变形量的估算方法及*** | |
CN106323587B (zh) | 机翼风洞试验模型弹性变形的单目视频高精度测量方法 | |
CN108959702A (zh) | 一种基于面偏差传递的火箭贮箱偏差分析方法 | |
CN108229009A (zh) | 设计装配零件模型的生成方法、装置及*** | |
CN115438535A (zh) | 一种利用有限元综合分析结果驱动模具加工数据型面变形的方法 | |
Turlier et al. | FEA shell element model for enhanced structural stress analysis of seam welds | |
CN114880904A (zh) | 一种橡胶衬套的大规模有限元分析方法 | |
KR20130101792A (ko) | 구조물의 좌굴 시뮬레이션 방법 |
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 |