CN111597742A - 一种基于焊接构件联动效应的农机结构残余变形预测方法 - Google Patents
一种基于焊接构件联动效应的农机结构残余变形预测方法 Download PDFInfo
- Publication number
- CN111597742A CN111597742A CN202010309792.6A CN202010309792A CN111597742A CN 111597742 A CN111597742 A CN 111597742A CN 202010309792 A CN202010309792 A CN 202010309792A CN 111597742 A CN111597742 A CN 111597742A
- Authority
- CN
- China
- Prior art keywords
- linkage
- agricultural machinery
- damaged
- machinery structure
- value
- 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
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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/28—Fuselage, exterior or interior
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种基于焊接构件联动效应的农机结构残余变形预测方法,包括如下步骤:S1、农机结构中焊接部分危险区域的确定;S2、谐响应冲击下损伤区域之间联动因子的确定;S3、焊接完毕的损伤区域残余应力的测量;S4、不同损伤区域的联动递进函数的确定;S5、基于联动效应的农机结构残余变形的预测。该方法预测精度高,可准确、高效实现农机结构残余变形的预测。
Description
技术领域
本发明涉及农机结构残余变形预测方法,特别涉及一种基于焊接构件联动效应的农机结构残余变形预测方法。
背景技术
随着农业经济的迅猛发展,越来越多同种类的农机装备出现在农耕作业中。农机装备的结构大多是通过不同零部件焊接组成的,焊接工艺在农机结构的总制造工序中占据了较大的比例。但是由于焊接工艺的加入,多个结构件的关键位置势必会产生一定的残余应力,从而造成结构出现一定程度的残余变形,进而影响到农耕作业的质量。目前对农耕机械结构残余变形的获取大多是通过实测方式获得的,且往往是单一零件变形量的检测,而没有考虑不同结构之间焊接联动关系对整体农机结构残余变形的影响。因此,从结构中所有焊接零部件之间的焊接关系出发来预测农机结构的残余变形量是必要的。
发明内容
发明目的:本发明目的是提供一种基于焊接构件联动效应的农机结构残余变形预测方法。
技术方案:本发明提供一种基于焊接构件联动效应的农机结构残余变形预测方法,包括如下步骤:
S1、农机结构中焊接部分危险区域的确定;
S2、谐响应冲击下损伤区域之间联动因子的确定;
S3、焊接完毕的损伤区域残余应力的测量;
S4、不同损伤区域的联动递进函数的确定;
S5、基于联动效应的农机结构残余变形的预测。
进一步地,所述S1农机结构中焊接部分危险区域的确定方法:在有限元分析软件ANSYS中导入农机结构的三维模型,采用自动网格模式对农机结构三维模型进行网格化划分,在前处理模块中对不同零部件之间的连接焊缝运用焊缝模块进行定义,然后定义工况的载荷条件与约束条件,启用后处理模块对农机结构中焊接零部件中应力值及其分布情况进行分析,确定农机整体结构在工况条件下发生损坏的区域。
进一步地,所述S2谐响应冲击下损伤区域之间联动因子的确定方法:在有限元分析软件中导入农机结构的三维模型,选用谐响应分析模块,采用自动网格模式对整体结构的网格进行划分,然后采用加倍细化的处理方法对S1中确定的发生损坏的区域进行网格的细化,并对这些发生损坏的区域进行编号i,i=1…e,e为损伤区域的个数,在设置谐响应冲击载荷中定义载荷的施加位置为S1中分析的结果中应力最小的区域,然后分析每个损坏区域的谐响应频率值记为fi,同时,在模态分析模块中同理分析农机结构的前三阶固有频率值标记为W1,W2,W3,计算每个损伤区域的频率共享系数Xi,
其中,Xi为每个损伤区域的频率共享系数;fi为每个损伤区域的谐响应频率值;W1为农机结构第1阶固有频率;W2为农机结构第2阶固有频率;W3为农机结构第3阶固有频率,进而计算损伤区域之间的联动因子Li,
其中,Li为农机结构中损伤区域之间的联动因子;Xi为每个损伤区域的频率共享系数;e为损伤区域的个数;t为工况过程的载荷时间点;M为一个工况周期的总时间值;i为发生损坏的区域。
进一步地,所述S3焊接完毕的损伤区域残余应力的测量方法:采用X射线残余应力测量装置对S1中确定的三维残余应力值进行测量,标记结果为σi,然后计算不同损伤区域间的残余应力联动系数Ci,
其中,Ci为农机结构中不同损伤区域间的残余应力联动系数;σi为损伤区域的三维残余应力值;M为一个工况周期的总时间值;t为工况过程的载荷时间点;i为发生损坏的区域;e为损伤区域的个数。
进一步地,所述S4不同损伤区域的联动递进函数的确定方法:在S2与S3的基础上,计算综合联动因子Zi的值,
其中,Zi为综合联动因子;Li为农机结构中损伤区域之间的联动因子;β为三维晶体残余应力方向适用值,农机结构取值为1.2;Ci为农机结构中不同损伤区域间的残余应力联动系数,然后在MATLAB软件中对每个损伤区域的Zi按照编号进行函数关系的拟合,获得5次条件下的联动递进函数方程F(i),
F(i)=α1i5+α2i4+α3i3+α4i2+α5i+α6
其中,F(i)为5次条件下的联动递进函数方程;i为发生损坏的区域;α1为联动递进函数中五次项的方程拟合系数;α2为联动递进函数中四次项的方程拟合系数;α3为联动递进函数中三次项的方程拟合系数;α4为联动递进函数中二次项的方程拟合系数;α5为联动递进函数中一次项的方程拟合系数;α6为常数项。
进一步地,所述S5基于联动效应的农机结构残余变形的预测方法:将S3中测量得到的残余应力值σi与联动递进函数一起代入下式对某个损伤位置i处的残余变形值εi进行计算,
其中,εi为损伤位置i处的残余变形值;σi为损伤区域的三维残余应力值;γ为弹性模量,取值为2.06×1011Pa;e为损伤区域的个数;F(i)为5次条件下的联动递进函数方程;i为发生损坏的区域。
有益效果:本发明可以实现农机结构残余变形的实时预测,通过零部件焊接联动因子来表达残余变形在不同零件之间的传递,建立反映残余变形集中递进的函数,并进一步通过不同递进函数的联动集聚效应来确定农机结构的残余变形,精确、高效,从而更有效的确保农机结构的作业精准度与使用安全性。
附图说明
图1为本发明方法的工艺流程图。
具体实施方式
如图1所示,本实施例基于焊接构件联动效应的农机结构残余变形预测方法,包括以下步骤:
S1.农机结构中焊接部分危险区域的确定:
在有限元分析软件ANSYS中导入农机结构的三维模型,采用自动网格模式对农机结构三维模型进行网格化划分,在前处理模块中对不同零部件之间的连接焊缝运用焊缝模块进行定义,然后定义工况的载荷条件与约束条件,启用后处理模块对农机结构中焊接零部件中应力值及其分布情况进行分析,确定农机整体结构在工况条件下最发生损坏的区域。
S2.谐响应冲击下损伤区域之间联动因子的确定:
在有限元分析软件中导入农机结构的三维模型,选用谐响应分析模块,采用自动网格模式对整体结构的网格进行粗略划分,然后采用加倍细化的处理方法对S1中确定的最发生损坏的区域进行网格的细化。并对这些发生损坏的区域进行编号i(i=1…e,e为损伤区域的个数)。在设置谐响应冲击载荷中定义载荷的施加位置为S1中分析的结果中应力最小的区域,载荷的大小为100N,然后分析每个损坏区域的谐响应频率值记为fi,同时,在模态分析模块中同理分析农机结构的前三阶固有频率值标记为W1,W2,W3,计算每个损伤区域的频率共享系数Xi,
其中,Xi为每个损伤区域的频率共享系数;fi为每个损伤区域的谐响应频率值;W1为农机结构第1阶固有频率;W2为农机结构第2阶固有频率;W3为农机结构第3阶固有频率,进而计算损伤区域之间的联动因子Li,
其中,Li为农机结构中损伤区域之间的联动因子;Xi为每个损伤区域的频率共享系数;e为损伤区域的个数;t为工况过程的载荷时间点;M为一个工况周期的总时间值;i为发生损坏的区域。
S3.焊接完毕的损伤区域残余应力的测量:
采用X射线残余应力测量装置对S1中确定的三维残余应力值进行测量,标记结果为σi,然后计算不同损伤区域间的残余应力联动系数Ci,
其中,Ci为农机结构中不同损伤区域间的残余应力联动系数;σi为损伤区域的三维残余应力值;M为一个工况周期的总时间值;t为工况过程的载荷时间点;i为发生损坏的区域;e为损伤区域的个数。
S4.不同损伤区域的联动递进函数的确定:
在S2与S3的基础上,计算综合联动因子Zi的值,
其中,Zi为综合联动因子;Li为农机结构中损伤区域之间的联动因子;β为三维晶体残余应力方向适用值,农机结构取值为1.2;Ci为农机结构中不同损伤区域间的残余应力联动系数。
然后在MArLAB软件中对每个损伤区域的Zi按照编号进行函数关系的拟合,获得5次条件下的联动递进函数方程F(i),
F(i)=α1i5+α2i4+α3i3+α4i2+α5i+α6 (5)
其中,F(i)为5次条件下的联动递进函数方程;i为发生损坏的区域;α1为联动递进函数中五次项的方程拟合系数;α2为联动递进函数中四次项的方程拟合系数;α3为联动递进函数中三次项的方程拟合系数;α4为联动递进函数中二次项的方程拟合系数;α5为联动递进函数中一次项的方程拟合系数;α6为常数项。
S5.基于联动效应的农机结构残余变形的预测:
将S3中测量得到的残余应力值σi与联动递进函数一起代入下式对某个损伤位置i处的残余变形值εi进行计算,
其中,εi为损伤位置i处的残余变形值;σi为损伤区域的三维残余应力值;γ为弹性模量,取值为2.06×1011Pa;e为损伤区域的个数;F(i)为5次条件下的联动递进函数方程;i为发生损坏的区域。
Claims (6)
1.一种基于焊接构件联动效应的农机结构残余变形预测方法,其特征在于:包括如下步骤:
S1、农机结构中焊接部分危险区域的确定;
S2、谐响应冲击下损伤区域之间联动因子的确定;
S3、焊接完毕的损伤区域残余应力的测量;
S4、不同损伤区域的联动递进函数的确定;
S5、基于联动效应的农机结构残余变形的预测。
2.根据权利要求1所述的基于焊接构件联动效应的农机结构残余变形预测方法,其特征在于:所述S1农机结构中焊接部分危险区域的确定方法:在有限元分析软件ANSYS中导入农机结构的三维模型,采用自动网格模式对农机结构三维模型进行网格化划分,在前处理模块中对不同零部件之间的连接焊缝运用焊缝模块进行定义,然后定义工况的载荷条件与约束条件,启用后处理模块对农机结构中焊接零部件中应力值及其分布情况进行分析,确定农机整体结构在工况条件下发生损坏的区域。
3.根据权利要求2所述的基于焊接构件联动效应的农机结构残余变形预测方法,其特征在于:所述S2谐响应冲击下损伤区域之间联动因子的确定方法:在有限元分析软件中导入农机结构的三维模型,选用谐响应分析模块,采用自动网格模式对整体结构的网格进行划分,然后采用加倍细化的处理方法对S1中确定的发生损坏的区域进行网格的细化,并对这些发生损坏的区域进行编号i,i=1...e,e为损伤区域的个数,在设置谐响应冲击载荷中定义载荷的施加位置为S1中分析的结果中应力最小的区域,然后分析每个损坏区域的谐响应频率值记为fi,同时,在模态分析模块中同理分析农机结构的前三阶固有频率值标记为W1,W2,W3,计算每个损伤区域的频率共享系数Xi,
其中,Xi为每个损伤区域的频率共享系数;fi为每个损伤区域的谐响应频率值;W1为农机结构第1阶固有频率;W2为农机结构第2阶固有频率;W3为农机结构第3阶固有频率,进而计算损伤区域之间的联动因子Li,
其中,Li为农机结构中损伤区域之间的联动因子;Xi为每个损伤区域的频率共享系数;e为损伤区域的个数;t为工况过程的载荷时间点;M为一个工况周期的总时间值;i为发生损坏的区域。
5.根据权利要求4所述的基于焊接构件联动效应的农机结构残余变形预测方法,其特征在于:所述S4不同损伤区域的联动递进函数的确定方法:在S2与S3的基础上,计算综合联动因子Zi的值,
其中,Zi为综合联动因子;Li为农机结构中损伤区域之间的联动因子;β为三维晶体残余应力方向适用值,农机结构取值为1.2;Ci为农机结构中不同损伤区域间的残余应力联动系数,然后在MATLAB软件中对每个损伤区域的Zi按照编号进行函数关系的拟合,获得5次条件下的联动递进函数方程F(i),
F(i)=α1i5+α2i4+α3i3+α4i2+α5i+α6
其中,F(i)为5次条件下的联动递进函数方程;i为发生损坏的区域;α1为联动递进函数中五次项的方程拟合系数;α2为联动递进函数中四次项的方程拟合系数;α3为联动递进函数中三次项的方程拟合系数;α4为联动递进函数中二次项的方程拟合系数;α5为联动递进函数中一次项的方程拟合系数;α6为常数项。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010309792.6A CN111597742B (zh) | 2020-04-17 | 2020-04-17 | 一种基于焊接构件联动效应的农机结构残余变形预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010309792.6A CN111597742B (zh) | 2020-04-17 | 2020-04-17 | 一种基于焊接构件联动效应的农机结构残余变形预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111597742A true CN111597742A (zh) | 2020-08-28 |
CN111597742B CN111597742B (zh) | 2023-05-23 |
Family
ID=72189039
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010309792.6A Active CN111597742B (zh) | 2020-04-17 | 2020-04-17 | 一种基于焊接构件联动效应的农机结构残余变形预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111597742B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113177337A (zh) * | 2021-04-20 | 2021-07-27 | 扬州大学 | 基于关联因素特征值波动区间的芦苇收割机安全评估方法 |
CN114002332A (zh) * | 2021-09-29 | 2022-02-01 | 西安交通大学 | 一种结构损伤监测预警方法及结构完整性数字孪生系统 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU897441A1 (ru) * | 1979-12-12 | 1982-01-15 | Всесоюзный Проектно-Конструкторский Институт Сварочного Производства | Способ снижени остаточных напр жений и деформаций в сварных конструкци х |
JP2004122222A (ja) * | 2002-10-07 | 2004-04-22 | Toshiba Corp | 溶接最適化システムおよび溶接最適化方法 |
CN104625403A (zh) * | 2015-01-27 | 2015-05-20 | 湖南大学 | 一种抑制万瓦级激光焊接厚板上表面缺陷的方法 |
CN106312299A (zh) * | 2016-08-25 | 2017-01-11 | 广东工业大学 | 航空发动机支架激光喷丸校形形状精度在线控制的方法与装置 |
CN106485029A (zh) * | 2016-10-27 | 2017-03-08 | 北京市市政工程研究院 | 基于残余应变的钢筋砼梁桥损伤后承载能力评估方法 |
CN107220410A (zh) * | 2017-05-03 | 2017-09-29 | 中车青岛四方机车车辆股份有限公司 | 参量对焊接残余应力及变形影响敏感度获取方法 |
CN109926767A (zh) * | 2019-02-26 | 2019-06-25 | 中国石油大学(华东) | 超大结构的焊接残余应力及变形的集成计算方法 |
CN110472283A (zh) * | 2019-07-12 | 2019-11-19 | 中国核动力研究设计院 | 一种用于估算环形焊接结构体积残余应力的无损方法 |
CN110598357A (zh) * | 2019-09-25 | 2019-12-20 | 华中科技大学 | 一种焊接接头应力变形仿真方法、装置、设备及存储介质 |
-
2020
- 2020-04-17 CN CN202010309792.6A patent/CN111597742B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU897441A1 (ru) * | 1979-12-12 | 1982-01-15 | Всесоюзный Проектно-Конструкторский Институт Сварочного Производства | Способ снижени остаточных напр жений и деформаций в сварных конструкци х |
JP2004122222A (ja) * | 2002-10-07 | 2004-04-22 | Toshiba Corp | 溶接最適化システムおよび溶接最適化方法 |
CN104625403A (zh) * | 2015-01-27 | 2015-05-20 | 湖南大学 | 一种抑制万瓦级激光焊接厚板上表面缺陷的方法 |
CN106312299A (zh) * | 2016-08-25 | 2017-01-11 | 广东工业大学 | 航空发动机支架激光喷丸校形形状精度在线控制的方法与装置 |
CN106485029A (zh) * | 2016-10-27 | 2017-03-08 | 北京市市政工程研究院 | 基于残余应变的钢筋砼梁桥损伤后承载能力评估方法 |
CN107220410A (zh) * | 2017-05-03 | 2017-09-29 | 中车青岛四方机车车辆股份有限公司 | 参量对焊接残余应力及变形影响敏感度获取方法 |
CN109926767A (zh) * | 2019-02-26 | 2019-06-25 | 中国石油大学(华东) | 超大结构的焊接残余应力及变形的集成计算方法 |
CN110472283A (zh) * | 2019-07-12 | 2019-11-19 | 中国核动力研究设计院 | 一种用于估算环形焊接结构体积残余应力的无损方法 |
CN110598357A (zh) * | 2019-09-25 | 2019-12-20 | 华中科技大学 | 一种焊接接头应力变形仿真方法、装置、设备及存储介质 |
Non-Patent Citations (2)
Title |
---|
BEZERRA A C 等: "On the influence of welding residual stresses on the dynamic behavior of structure", 《SHOCK & VIBRATION》 * |
朱林 等: "基于贝叶斯理论的结构件健康状态评估方法研究", 《振动与冲击》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113177337A (zh) * | 2021-04-20 | 2021-07-27 | 扬州大学 | 基于关联因素特征值波动区间的芦苇收割机安全评估方法 |
CN113177337B (zh) * | 2021-04-20 | 2023-05-26 | 扬州大学 | 基于关联因素特征值波动区间的芦苇收割机安全评估方法 |
CN114002332A (zh) * | 2021-09-29 | 2022-02-01 | 西安交通大学 | 一种结构损伤监测预警方法及结构完整性数字孪生系统 |
CN114002332B (zh) * | 2021-09-29 | 2023-07-25 | 西安交通大学 | 一种结构损伤监测预警方法及结构完整性数字孪生系统 |
Also Published As
Publication number | Publication date |
---|---|
CN111597742B (zh) | 2023-05-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3267167B1 (en) | Residual stress estimation method and residual stress estimation device | |
CN107766254B (zh) | 一种基于层次分析的软件质量评估方法及系统 | |
US10422733B2 (en) | Method and device for testing wettability of dense oil reservoir | |
US20120259593A1 (en) | Method for the prediction of fatigue life for welded structures | |
Dupuis et al. | Robust estimation of the Birnbaum-Saunders distribution | |
CN111597742A (zh) | 一种基于焊接构件联动效应的农机结构残余变形预测方法 | |
US8756028B2 (en) | Fault detection method of semiconductor manufacturing processes and system architecture thereof | |
JP3094191B2 (ja) | プラントの自己学習診断、予測方法及び装置 | |
US8346584B2 (en) | Method and apparatus for determining design modification effects on a computerized engineering model | |
Sonmez et al. | Using accurately measured production amounts to obtain calibration curve corrections of production line speed and stoppage duration consisting of measurement errors | |
Cook et al. | Examining performance differences in workload execution phases | |
Yang et al. | Crack identification driven by the fusion of mechanism and data for the variable-cross-section cantilever beam | |
EP2738701A2 (en) | Method and apparatus for statistical electronic circuit simulation | |
CN112329108B (zh) | 一种地铁车站优化抗浮验算方法及系统 | |
Abollado et al. | Uncertainty assessment for measurement processes in the aerospace manufacturing industry | |
CN110823443B (zh) | 一种静力试验应变数据处理方法 | |
CN109858699B (zh) | 水质定量模拟方法、装置、电子设备及存储介质 | |
CN110008560B (zh) | 叶片动应力测点确定方法 | |
Giagopoulos et al. | Dynamic response estimation and fatigue prediction in a linear substructure of a complex mechanical assembly | |
Sinha | Efficient wall thickness analysis methods for optimal design of casting parts | |
KR20230006308A (ko) | 금형수명 예측 시스템 및 방법 | |
Xue et al. | Unified Elastic–Plastic Analytical Method for Estimating Notch Local Strains in Real Time under Multiaxial Irregular Loading | |
Veltri | FEM techniques for high stress detection in accelerated fatigue simulation | |
CN111339699A (zh) | 基于宏观应力与残余应力拓度的箱型结构重变形量预测方法 | |
JP6665023B2 (ja) | 躯体の変位推定モニタリングシステム。 |
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 |