CN109325292B - 一种基于ansys的复杂焊接轨迹热源加载局部坐标系设定方法 - Google Patents

一种基于ansys的复杂焊接轨迹热源加载局部坐标系设定方法 Download PDF

Info

Publication number
CN109325292B
CN109325292B CN201811114131.7A CN201811114131A CN109325292B CN 109325292 B CN109325292 B CN 109325292B CN 201811114131 A CN201811114131 A CN 201811114131A CN 109325292 B CN109325292 B CN 109325292B
Authority
CN
China
Prior art keywords
welding
heat source
line
local coordinate
coordinate system
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.)
Active
Application number
CN201811114131.7A
Other languages
English (en)
Other versions
CN109325292A (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.)
Liaoning Shihua University
Original Assignee
Liaoning Shihua 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 Liaoning Shihua University filed Critical Liaoning Shihua University
Priority to CN201811114131.7A priority Critical patent/CN109325292B/zh
Publication of CN109325292A publication Critical patent/CN109325292A/zh
Application granted granted Critical
Publication of CN109325292B publication Critical patent/CN109325292B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Arc Welding In General (AREA)
  • Numerical Control (AREA)

Abstract

本发明公开了一种基于ANSYS的复杂焊接轨迹热源加载局部坐标系设定方法。该方法的主要内容包括，在给定焊接轨迹线、焊枪轨迹线、焊接起点和热源移动步长的前提下，利用ANSYS的内部命令，自动计算出热源移动过程中以热源中心为原点的动态坐标系各坐标轴的实时方向，为热源加载提供准确的局部坐标系。在复杂焊接轨迹的情况下，热源加载方向是实时变动的，本发明通过焊枪轨迹线的辅助作用，严格确定局部坐标系的坐标轴方向，从而定位热源加载方向。利用本发明公开的局部坐标设定方法，使焊接模拟能够很方便地进行复杂轨迹线的模拟，增大了焊接模拟的适用范围，增强了焊接模拟的实际应用性。

Description

一种基于ANSYS的复杂焊接轨迹热源加载局部坐标系设定 方法
技术领域
本发明涉及有限元仿真技术领域，尤其涉及一种基于ANSYS的复杂焊接轨迹热源加载局部坐标系设定方法。
背景技术
在焊接模拟过程中，为了实现焊接热源移动的模拟，需要将整个焊接轨迹线按照焊接模拟所设定的空间步长分段，在每段上建立热源加载所需的局部坐标系。通过文献查阅，焊接模拟的研究通常局限于直线或圆弧等规则轨迹或规则轨迹的组合，随着热源中心的移动，局部坐标系原点做相应的改变，坐标轴方向在移动过程中保持不变或做相应的有规律的转向。而实际的焊接过程中焊接轨迹可能是很复杂的曲线，甚至是空间不规则曲线。本发明所公开的一种基于ANSYS的复杂焊接轨迹热源加载局部坐标系设定方法，正是解决上述技术问题，从而可有效地解决复杂焊接轨迹焊接模拟中热源加载的难题，提高焊接模拟的适用范围，兼容实践中复杂情况下的焊接。
发明内容
本发明所要解决的技术问题，就是提供一种基于ANSYS的复杂焊接轨迹热源加载局部坐标系设定方法，本方法使焊接模拟能够很方便地进行复杂轨迹线的模拟，增大了焊接模拟的适用范围，增强了焊接模拟的实际应用性。
本发明采用的技术解决方案如下:
（1）准备工作
本发明需要在以下工作完成的情况下才能实施：
已进行符合焊接模拟要求的合理的有限元网格的剖分，通常为了平衡计算精度和计算时间，焊缝区的网格划分得密一些，远离焊缝区疏松一些，以及其他的与焊接模拟有关的设定，如散热面的指定，装夹点的设定，各部分材料属性的设定，环境温度、预热温度的设定等。
绘制焊接轨迹线和焊枪轨迹线。焊枪轨迹线是指焊枪上的一固定点在整个焊接过程中随着热源移动所形成的轨迹线。
在实际的焊接中，焊枪轨迹线可能不平行焊接轨迹线，本发明也适用。通常情况下，可基于焊接轨迹线采用曲线偏移的命令获得焊枪轨迹线。
再有一个必备的参数是焊接热源沿着焊接轨迹线移动的步长需要给定。
（2）焊接轨迹热源加载局部坐标系设定方法
a. 起始点的指定
指定焊接轨迹线上的焊接起始点，包括焊枪轨迹线上的起始点，比如一个完整的圆，必须要有个开始点。可通过拾取的方式，或通过对话框的方式填写开始点的编号，来指明焊接轨迹和焊枪轨迹的起始点。
b. 焊接轨迹线和焊枪轨迹线的指定和长度的计算
因为Ansys软件曲线实体的特点，复杂的曲线是由多段线的组合而成，例如，圆就是由四段圆弧线组成，因此，对于复杂曲线的焊接轨迹，要按照焊接顺序依次点选，同样，焊接轨迹的选择也要按照焊接顺序依次点选。或者通过对话框的方式填写曲线组合段的数量和编号。
如下面通过循环命令获取焊接轨迹线的分线段的编号和线段的长度：
*do,i,1,WNUM ! WNUM为组成焊接轨迹的分线段的数量
LSEL,S,P !通过拾取选择焊接线
*GET,WL(i),LINE,,NUM,MAX !通过拾取的方式获得线的编号
*GET,wleng(i),LINE,wl(i),LENG !提取线的长度
……
*enddo
同样可通过类似的命令获取焊枪轨迹线的分线段的编号和线段的长度。
将循环过程获得的分线段求和，计算焊接线和焊枪轨迹线的总长。
c. 计算局部坐标系的原点坐标
计算分段加载法焊接热源的中心坐标
根据设定的移动的空间步长，计算焊接线的分段数k。然后利用循环过程及下列组合命令计算焊接线各空间步长分段起点处的坐标。
frac=(k-1)*step/WLENG(i) !从焊接起点到上第k段步长起点处分数
wxc(k)=LX(WL(i),frac) !焊接线编号为WL(i)段的长度百分数为frac处的x坐标值
d. 计算局部坐标轴相对于全局坐标轴的转角
根据坐标差的比值计算局部坐标轴相对于全局坐标轴的转角。
如局部坐标轴的x轴向全局坐标轴的y轴方向的转角：当wyc(k+1)-wyc(k)>0时，THXY(k) = -ATAN ((wxc(k+1) - wxc(k)) / (wyc(k+1) - wyc(k)))
同理，可计算局部坐标轴的y轴向全局坐标轴的z轴方向的转角，局部坐标轴的z轴向全局坐标轴的x轴方向的转角。
e. 局部坐标系的设定
LOCAL,c,0,wxc(k),wyc(k),wzc(k),THXY(k),THYZ(k),THZX(k),1,1,
由此确定了焊接热源沿着焊接轨迹及焊枪运动轨迹局部坐标的原点及坐标轴的方向。
此方法可处理复杂的焊接轨迹线。
上述办法就是本发明公开的一种基于ANSYS的复杂焊接轨迹热源加载局部坐标系设定方法。后续的热源沿着焊接轨迹的加载，可利用ANSYS的循环命令，基于上述局部坐标系进行循环加载和卸载与焊接方法匹配的热源模型（如高斯热源，椭圆热源，双椭圆模型等，或更为复杂的组合热源），实现焊接热源移动的模拟。
与现有技术相比，本发明的有益效果为：
利用本发明公开的局部坐标设定方法，使焊接模拟能够很方便地进行复杂轨迹线的模拟，增大了焊接模拟的适用范围，增强了焊接模拟的实际应用性。
附图说明
图1 是本发明实施例中的焊接轨迹线，即焊接轨迹线为三段相接的圆弧。
图2 是本发明实施例中焊接轨迹热源加载的局部坐标系，及不同时刻焊接热源的运动位置及温度场分布图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细说明，但本发明的实施方式不限于此，在不脱离本发明上述思想的情况下，根据本领域普通技术知识和惯用手段，做出各种替换和变更，均应包括在本发明的范围内。
实施例：
本实例公开了一种基于ANSYS的复杂焊接轨迹热源加载局部坐标系设定方法。本实例的焊接轨迹线如图1所示，为平板表面三段弧线连接而成。
本实例拟进行焊接模拟的试板为5mm厚的平板，已进行符合焊接模拟要求的合理的有限元网格的剖分，即为了平衡计算精度和计算时间，焊接轨迹线附近的网格划分得密一些，远离焊接轨迹线的区域划分得疏松一些，本实例为了简化采用统一的网格尺寸1mm，散热条件为对流散热，焊接热源采用双椭球热源模型，材料采用普通碳素结构钢。
焊枪轨迹线位在试板的上方，与焊接轨迹线平行，两线相距5mm。
本实例首先要指定焊接轨迹线和焊枪轨迹线上的焊接起始点，本实例通过点选的方法获取焊接轨迹线和焊枪轨迹线上的焊接起始点的编号：
NSEL,S,P !拾取焊接起点
*GET,Ws,KP,,NUM,MAX !通过拾取的方式获得节点的编号
……
本实例的焊接轨迹线和焊枪轨迹线均由三段圆弧线连接而成，APDL程序需要知道组成线段的顺序和编号，本实例通过循环命令获取焊接轨迹线和焊枪轨迹线的分线段的编号和线段的长度：
*do,i,1,WNUM ! WNUM为组成焊接轨迹的分线段的数量
LSEL,S,P !通过拾取选择焊接线
*GET,WL(i),LINE,,NUM,MAX !通过拾取的方式获得线的编号
*GET,wleng(i),LINE,wl(i),LENG !提取线的长度
……
*enddo
将循环过程获得的分线段求和，计算焊接轨迹线和焊枪轨迹线的总长。
根据设定的空间步长，计算焊接轨迹线和焊枪轨迹上为了加载焊接热源设定的分段数k。然后利用循环过程及下列组合命令计算焊接线各空间步长分段起点处的坐标。
frac=(k-1)*step/WLENG(i) !从焊接起点到上第k段步长起点处分数
wxc(k)=LX(WL(i),frac) !焊接线编号为WL(i)段的长度百分数为frac处的x坐标值
为了保证热源运行趋势与焊接轨迹线相切及焊枪瞬时方位，本实例根据步长起始点的坐标差的比值计算局部坐标轴相对于全局坐标轴的转角。
如局部坐标轴的x轴向全局坐标轴的y轴方向的转角：当wyc(k+1)-wyc(k)>0时，THXY(k) = -ATAN ((wxc(k+1) - wxc(k)) / (wyc(k+1) - wyc(k)))
同理，可计算局部坐标轴的y轴向全局坐标轴的z轴方向的转角，局部坐标轴的z轴向全局坐标轴的x轴方向的转角。
根据上面计算获得第k段步长起点处的坐标值和转角值，设定热源在该点处加载的局部坐标系：
LOCAL,c,0,wxc(k),wyc(k),wzc(k),THXY(k),THYZ(k),THZX(k),1,1,
由此确定了焊接热源在某一瞬时在焊接轨迹线上加载位置及方位，如图二所示。图2展示了焊接热源在焊接轨迹线上加载的局部坐标系，及不同时刻焊接热源的运动位置及温度场分布。

Claims (3)

1.一种基于ANSYS的复杂焊接轨迹热源加载局部坐标系设定方法，其特征在于该方法包括下列主要步骤：
a.首先指定焊接轨迹线和焊枪轨迹线上的焊接起始点，可通过拾取的方式，或通过对话框的方式填写开始点的编号；
b.焊接轨迹线和焊枪轨迹线的指定和长度的计算：
因为在Ansys中复杂的曲线是由多段线的组合而成，因此，对于复杂曲线的焊接轨迹，要按照焊接顺序依次点选，或者通过对话框的方式填写曲线组合段的数量和编号，下面的代码段是通过循环命令获取焊接轨迹线的分线段的编号和线段的长度：
*do,i,1,WNUM！WNUM为组成焊接轨迹的分线段的数量
LSEL,S,P！通过拾取选择焊接线
*GET,WL(i),LINE,,NUM,MAX！通过拾取的方式获得线的编号WL(i)
*GET,wleng(i),LINE,wl(i),LENG！提取线编号为wl(i)的长度wleng(i)
*enddo
同理，获取焊枪轨迹线的分线段的编号和线段的长度，将循环过程获得的分线段求和，计算焊接线和焊枪轨迹线的总长；
c.计算局部坐标系的原点坐标：
计算分段加载法焊接热源的中心坐标，根据设定的移动的空间步长step，计算焊接线的分段数k，然后利用循环过程及下列组合命令计算焊接线各空间步长分段起点处的坐标：
frac＝(k-1)*step/WLENG(i)！从焊接起点到上第k段步长起点处分数frac
wxc(k)＝LX(WL(i),frac)！焊接线编号为WL(i)段的长度百分数为frac处的x坐标值wxc(k)
d.计算局部坐标轴相对于全局坐标轴的转角：
根据坐标差的比值计算局部坐标轴相对于全局坐标轴的转角，如局部坐标轴的x轴向全局坐标轴的y轴方向的转角THXY(k)：当wyc(k+1)-wyc(k)>0时，THXY(k)＝-ATAN((wxc(k+1)-wxc(k))/(wyc(k+1)-wyc(k)))，同理，可计算局部坐标轴的y轴向全局坐标轴的z轴方向的转角THYZ(k)，局部坐标轴的z轴向全局坐标轴的x轴方向的转角THZX(k)，
e.局部坐标系的设定：
LOCAL,c,0,wxc(k),wyc(k),wzc(k),THXY(k),THYZ(k),THZX(k),1,1,！c为局部坐标系编号
由此确定了焊接热源沿着焊接轨迹及焊枪运动轨迹局部坐标的原点及坐标轴的方向。
2.权利要求1所述的基于ANSYS的复杂焊接轨迹热源加载局部坐标系设定方法，复杂焊接轨迹包括空间不规则曲线的焊接轨迹。
3.权利要求1所述的基于ANSYS的复杂焊接轨迹热源加载局部坐标系设定方法，复杂焊接轨迹包括与焊接轨迹线非平行的焊枪轨迹线的焊接情况。
CN201811114131.7A 2018-09-25 2018-09-25 一种基于ansys的复杂焊接轨迹热源加载局部坐标系设定方法 Active CN109325292B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811114131.7A CN109325292B (zh) 2018-09-25 2018-09-25 一种基于ansys的复杂焊接轨迹热源加载局部坐标系设定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811114131.7A CN109325292B (zh) 2018-09-25 2018-09-25 一种基于ansys的复杂焊接轨迹热源加载局部坐标系设定方法

Publications (2)

Publication Number Publication Date
CN109325292A CN109325292A (zh) 2019-02-12
CN109325292B true CN109325292B (zh) 2023-01-10

Family

ID=65266418

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811114131.7A Active CN109325292B (zh) 2018-09-25 2018-09-25 一种基于ansys的复杂焊接轨迹热源加载局部坐标系设定方法

Country Status (1)

Country Link
CN (1) CN109325292B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113673124A (zh) * 2021-07-06 2021-11-19 华南理工大学 三通相贯线焊接温度场数值模拟预测方法、系统及介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013154451A1 (en) * 2012-04-12 2013-10-17 Siemens Aktiengesellschaft Method for a welding process control of nickel based superalloy products
CN103366044A (zh) * 2013-06-04 2013-10-23 广东电网公司电力科学研究院 一种基于ansys预测cmt焊缝形状尺寸的方法
CN104400279A (zh) * 2014-10-11 2015-03-11 南京航空航天大学 基于ccd的管道空间焊缝自动识别与轨迹规划的方法及系统
CN106649956A (zh) * 2016-09-26 2017-05-10 辽宁石油化工大学 一种基于轴测图的管道三维重建方法
CN106875783A (zh) * 2017-02-28 2017-06-20 深圳市为汉科技有限公司 虚拟焊接系统及方法
CN108563878A (zh) * 2018-04-19 2018-09-21 大连五州轨道科技有限公司 多功能的焊接结构疲劳寿命计算与焊接质量等级确认一体化的集成系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9851332B2 (en) * 2014-09-19 2017-12-26 King Fahd University Of Petroleum And Minerals Process for determining weld quality using flexural characteristics

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013154451A1 (en) * 2012-04-12 2013-10-17 Siemens Aktiengesellschaft Method for a welding process control of nickel based superalloy products
CN103366044A (zh) * 2013-06-04 2013-10-23 广东电网公司电力科学研究院 一种基于ansys预测cmt焊缝形状尺寸的方法
CN104400279A (zh) * 2014-10-11 2015-03-11 南京航空航天大学 基于ccd的管道空间焊缝自动识别与轨迹规划的方法及系统
CN106649956A (zh) * 2016-09-26 2017-05-10 辽宁石油化工大学 一种基于轴测图的管道三维重建方法
CN106875783A (zh) * 2017-02-28 2017-06-20 深圳市为汉科技有限公司 虚拟焊接系统及方法
CN108563878A (zh) * 2018-04-19 2018-09-21 大连五州轨道科技有限公司 多功能的焊接结构疲劳寿命计算与焊接质量等级确认一体化的集成系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于ANSYS的焊接温度场和应力的数值模拟研究;姬庆玲;《时代农机》;20161228(第12期);全文 *
基于ANSYS的输电杆塔塔脚主筋板焊接过程的温度场仿真;陈华等;《机械工程师》;20180410(第04期);全文 *

Also Published As

Publication number Publication date
CN109325292A (zh) 2019-02-12

Similar Documents

Publication Publication Date Title
US20130060369A1 (en) Method and system for generating instructions for an automated machine
CN109584288B (zh) 一种五轴系统中三维模型的重构方法及系统
CN111596614B (zh) 基于云边协同的运动控制误差补偿系统及方法
CN113634958A (zh) 一种基于三维视觉的大型结构件自动化焊接系统及方法
CN109325292B (zh) 一种基于ansys的复杂焊接轨迹热源加载局部坐标系设定方法
CN108490871A (zh) 四轴数控铣床加工方法、装置、计算机设备和存储介质
CN108858187B (zh) 一种带电接搭引流线作业现场的快速重构方法
Geng et al. A novel welding path planning method based on point cloud for robotic welding of impeller blades
CN102990178A (zh) 一种使用Gerber文件提高焊锡机器人编程速度和精度的方法
CN110171000A (zh) 坡口切割方法、装置及控制设备
CN104715109B (zh) 一种基于球干涉的汽车白车身焊点自动划分方法
CN109483545B (zh) 一种焊缝重构方法、智能机器人焊接方法及系统
CN110153582A (zh) 焊接方案生成方法、装置以及焊接系统
CN109581964A (zh) 控制装置、控制装置的控制方法以及记录媒体
CN107066654B (zh) 面向点云模型的包容盒面的轨迹控制点拾取方法和装置
CN104741791B (zh) 一种曲面工件表面阵列微结构图形布局方法
CN109530860B (zh) 一种可用于大尺寸零件制造的电弧增材方法
JP3601303B2 (ja) 有限要素生成方法
CN109812260A (zh) 用于海上勘探的定向井轨迹设计方法和系统、电子设备
CN112560276A (zh) 一种载荷谱的生成方法
CN113523651A (zh) 焊接系统及焊接参数优化方法
US7069102B2 (en) System and method to customize bond programs compensating integrated circuit bonder variability
JP2004106049A (ja) 三次元レーザ加工機による加工方法並びに三次元レーザ加工用のｎｃプログラムの作成方法
CN110722554B (zh) 一种基于激光点云数据的机械手轨迹编辑及纠正方法
CN113119123B (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