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
trajectory
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.)
Trelon Fushun Intelligent Technology Co ltd
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)
  • Numerical Control (AREA)
  • Arc Welding In General (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
CN102785246B (zh) 一种可实现自动轨迹修正的机器人标定方法
CN103885390B (zh) 基于逆向工程技术的免示教激光三维测量方法和设备
CN110411338B (zh) 机器人电弧增材修复的焊枪工具参数三维扫描标定方法
CN108332652A (zh) 一种农田作业面积计量方法及装置
Geng et al. A novel welding path planning method based on point cloud for robotic welding of impeller blades
CN109658524A (zh) 一种三维模型的编辑方法、系统及相关装置
CN109325292B (zh) 一种基于ansys的复杂焊接轨迹热源加载局部坐标系设定方法
WO2024149402A1 (zh) 多层多道焊接方法、装置、处理器和焊接系统
CN111761406B (zh) 基于在机测量的工件各自由度分步找正的快速装夹方法
CN106839987A (zh) 一种基于树干曲面模型的树干材积获取方法
CN106777561A (zh) 一种涡轮叶片精铸蜡型陶芯夹紧元件的布局方法
CN117876358A (zh) 焊接路径规划方法、装置、计算机设备和存储介质
CN106020194A (zh) 一种基于传热分析的工业机器人路径优化方法及系统
CN114998446A (zh) 相机与机械手之间的标定方法、控制器、设备及存储介质
CN110457860A (zh) 网格的自动生成、装置、设备和存储介质
CN108595886A (zh) 一种回弹补偿网格模型构建方法
CN109313819A (zh) 线路模型实现方法、装置和计算机可读存储介质
CN109530860B (zh) 一种可用于大尺寸零件制造的电弧增材方法
CN117428773A (zh) 一种焊接机器人的编程方法及系统
CN116833606A (zh) 一种焊接方法、装置、设备及存储介质
CN109812260A (zh) 用于海上勘探的定向井轨迹设计方法和系统、电子设备
CN112560276A (zh) 一种载荷谱的生成方法
CN115167278A (zh) 一种密集型点位加工的最优运动轨迹规划方法
CN113119123B (zh) 一种办公转椅生产的运动控制方法
CN104729440B (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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20241128

Address after: No. 102, Unit 5, Building 18, Jianshe Street, Wanghua District, Fushun City, Liaoning Province 113001

Patentee after: Trelon (Fushun) Intelligent Technology Co.,Ltd.

Country or region after: China

Address before: Wanghua District West Dandong Road 113001 Liaoning city of Fushun province No. 1

Patentee before: Liaoming Petrochemical University

Country or region before: China