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

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

Info

Publication number
CN109325292A
CN109325292A CN201811114131.7A CN201811114131A CN109325292A CN 109325292 A CN109325292 A CN 109325292A CN 201811114131 A CN201811114131 A CN 201811114131A CN 109325292 A CN109325292 A CN 109325292A
Authority
CN
China
Prior art keywords
welding
line
local coordinate
heat source
axis
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
Application number
CN201811114131.7A
Other languages
English (en)
Other versions
CN109325292B (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

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 !通过拾取的方式获得线的编号
*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,
由此确定了焊接热源沿着焊接轨迹及焊枪运动轨迹局部坐标的原点及坐标轴的方向。
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 true CN109325292A (zh) 2019-02-12
CN109325292B 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)

Cited By (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 (7)

* 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的管道空间焊缝自动识别与轨迹规划的方法及系统
US20160084802A1 (en) * 2014-09-19 2016-03-24 King Fahd University Of Petroleum And Minerals Process for determining weld quality using flexural characteristics
CN106649956A (zh) * 2016-09-26 2017-05-10 辽宁石油化工大学 一种基于轴测图的管道三维重建方法
CN106875783A (zh) * 2017-02-28 2017-06-20 深圳市为汉科技有限公司 虚拟焊接系统及方法
CN108563878A (zh) * 2018-04-19 2018-09-21 大连五州轨道科技有限公司 多功能的焊接结构疲劳寿命计算与焊接质量等级确认一体化的集成系统

Patent Citations (7)

* 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焊缝形状尺寸的方法
US20160084802A1 (en) * 2014-09-19 2016-03-24 King Fahd University Of Petroleum And Minerals Process for determining weld quality using flexural characteristics
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的焊接温度场和应力的数值模拟研究", 《时代农机》 *
陈华等: "基于ANSYS的输电杆塔塔脚主筋板焊接过程的温度场仿真", 《机械工程师》 *

Cited By (1)

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

Also Published As

Publication number Publication date
CN109325292B (zh) 2023-01-10

Similar Documents

Publication Publication Date Title
US9933255B2 (en) Process mapping of melt pool geometry
CN106041075B (zh) 一种金属零件悬空结构的高能束增材制造方法
CN102345390B (zh) 一种混凝土泵车臂架变形量的补偿方法
CN106799745A (zh) 一种基于协同克里金的工业机械臂精度标定方法
CN103983963B (zh) 一种多站地基激光雷达数据的自动配准方法
CN105528503B (zh) 一种基于结构分解的大型构件动态优化设计方法
CN110990913B (zh) 一种采用水平面内合力控制的边跨找形方法
CN107876938A (zh) 一种基于机械臂摆动的电弧增材制造道间搭接方法
CN101559511B (zh) 一种以温度为控制变量的焊接数值模拟计算方法
CN107330585A (zh) 一种船舶建造数据集成系统及辅助建造船舶的方法
CN108629453A (zh) 一种基于实测数据的飞机现场装配精度预测方法
CN109325292A (zh) 一种基于ansys的复杂焊接轨迹热源加载局部坐标系设定方法
CN109766597A (zh) 一种考虑几何误差的装配位姿高精度线性化求解方法
Wang et al. Design and implementation of five-axis transformation function in CNC system
CN106021650B (zh) 一种图元结构剪取与贴合的基准融合方法
CN111551178A (zh) 一种基于最短路径的分段轨迹时间规划方法
CN110807276A (zh) 导入残余应力分析铸件结构强度的方法及系统
CN110703689A (zh) 一种数控机床空间误差补偿方法及系统
CN110221313A (zh) Agv小车的激光雷达安装位置修正方法以及agv小车
CN110153582A (zh) 焊接方案生成方法、装置以及焊接系统
CN108376182A (zh) 一种加强增材制造微单元和蒙皮结构结合性能控制方法
CN110348072A (zh) 一种提高电弧增材构件热力分析有限元数值计算效率方法
Elayeb et al. Printing orientation selection based on the dimensional errors modelling in additive manufacturing
CN109886879A (zh) 图形镜像处理方法、电子设备及存储介质
CN105157656A (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