CN109614724A - 一种关于机床关重件残余应力水平的评估方法 - Google Patents

一种关于机床关重件残余应力水平的评估方法 Download PDF

Info

Publication number
CN109614724A
CN109614724A CN201811539939.XA CN201811539939A CN109614724A CN 109614724 A CN109614724 A CN 109614724A CN 201811539939 A CN201811539939 A CN 201811539939A CN 109614724 A CN109614724 A CN 109614724A
Authority
CN
China
Prior art keywords
creep
important parts
lathe
lathe important
residual stress
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
Application number
CN201811539939.XA
Other languages
English (en)
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.)
Chongqing University
Original Assignee
Chongqing 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 Chongqing University filed Critical Chongqing University
Priority to CN201811539939.XA priority Critical patent/CN109614724A/zh
Publication of CN109614724A publication Critical patent/CN109614724A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种关于机床关重件残余应力水平的评估方法,包括:步骤一,通过蠕变实验,获取关重件材料在常温25℃下的蠕变仿真参数;步骤二,将蠕变参数代入已知的机床关重件abaqus模型,运用有限元分析软件进行机床关重件的蠕变仿真分析;步骤三,对蠕变仿真进行修正;步骤四,根据机床关重件几何精度要求,判断出机床关重件蠕变后的几何精度是否要求;步骤五,根据机床关重件蠕变后几何精度是否符合要求,判断出机床关重件蠕变之前的残余内应力是否可靠。本发明能对机床关重件的残余应力水平的评估客观准确,解决了现有技术中缺少相关的机床零件残余应力的评判方法的问题,对指导企业生产工艺的改进、零件内应力合格检验具有很大的帮助。

Description

一种关于机床关重件残余应力水平的评估方法
技术领域
本发明涉及机械制造技术领域,特别涉及一种基于abaqus蠕变仿真的机床关重件残余内应力的风险评估方法。
背景技术
机床零件在铸造、热处理、焊接、锻压、切削加工等制造过程中,都会在零件内部或表面出现不同程度的残余应力。残余应力是影响机床零件特别是关重件尺寸稳定性、结构强度及疲劳寿命等指标的主要因素之一,同时也是引起变形、开裂等问题的罪魁祸首,对机床产品的可靠性、性能稳定性和精度寿命都有重要影响。然而机床零件的残余内应力往往是不可预测的,同种工序下所产生的残余应力都是不一样的,现在有许多测量零件局部的残余内应力的方法,但是对于机床零件整体的残余应力的描述,一般采用的是有限元仿真的方法。
机床零件的残余应力是充满不确定性的,如何去评估一个机床零件的残余应力是否合格,始终没有一个评估的标准,因为不同的零件结构,材料属性,使用场合,都会对评估结果产生影响。
发明内容
有鉴于此,本发明的目的是提供一种关于机床关重件残余应力水平的评估方法,其通过蠕变试验获取材料蠕变参数,结合abaqus软件,建立机床关重件蠕变有限元模型,提出蠕变模型修正方法,计算得出包括机床关重件蠕应力、蠕应变的仿真结果,整理仿真得出的机床关重件蠕变仿真几何误差,对比机床关重件几何误差要求标准,判断机床关重件几何精度是否符合标准,由此反向证明、判断机床关重件的残余内应力是否可靠。
本发明关于机床关重件残余应力水平的评估方法,关于机床关重件残余应力水平的评估方法,包括以下步骤:
步骤一,通过蠕变实验,获取关重件材料在常温25℃下的蠕变仿真参数;
步骤二,将蠕变参数代入已知的机床关重件abaqus模型,运用abaqus 有限元分析软件进行机床关重件的蠕变仿真分析;
步骤三,对机床关重件进行相关蠕变测试实验,对abaqus蠕变仿真进行修正;
步骤四,通过修正后的蠕变仿真分析得到由于机床关重件残余应力松弛、机床关重件产生蠕变所造成的几何误差,再将蠕变仿真产生的几何误差和机床关重件的初始几何误差进行叠加,再将叠加后的误差与机床关重件几何精度要求进行对比,判断出机床关重件蠕变后的几何精度是否符合要求;
步骤五,根据机床关重件蠕变后几何精度是否符合要求,判断出机床关重件蠕变之前的残余内应力是否可靠。
进一步,在步骤三中对abaqus蠕变仿真进行修正的步骤如下:
1)将机床关重件静置一段时间,静置时间为仿真设置时间的一半;
2)测出机床关重件在这一段时间内测点的几何精度变化;
3)利用abaqus软件对机床关重件在这一段时间的蠕变进行分析;
4)对实验测得的数据和仿真得到的相应的数据进行拟合修正,得到相应的修正公式或者修正系数;
5)将得到的修正系数或者公式用于对步骤二中的仿真结果进行修正。
本发明的有益效果:
本发明关于机床关重件残余应力水平的评估方法,通过蠕变试验获取蠕变仿真所需参数,再将数据代入abaqus软件进行机床关重件蠕变仿真,并通过修正方案提高仿真可信度,然后将机床关重件蠕变仿真结果的几何误差与机床关重件的初始误差叠加后与机床关重件几何精度标准进行对比,根据对比结果反向判断出机床关重件蠕变前残余应力是否可靠,其对机床关重件的残余应力水平的评估客观准确,解决了现有技术中无相关的机床零件残余应力的评判方法的问题,对指导企业生产工艺的改进、零件内应力合格检验具有很大的帮助。
附图说明
图1是实施例中关于机床关重件残余应力水平的评估方法的技术路线图。
图2是实施例中机床横梁初始模型残余应力分布云图。
图3是实施例中机床横梁蠕变仿真计算之后的应力分布云图。
图4是实施例中横梁蠕变仿真修正方法的流程图。
图5是实施例中横梁蠕变仿真上导轨面法向位移拟合曲线。
图6是实施例中横梁蠕变仿真下导轨面法向位移拟合曲线。
具体实施方式
下面结合实施例和附图,对本发明技术方案进行清楚、完整的描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本实施例例关于机床关重件残余应力水平的评估方法,包括以下步骤:
步骤一,通过蠕变实验,获取关重件材料在常温25℃下的蠕变仿真参数;本实施例中,关重件为机床横梁,机床横梁的材料为HT300。待获取的蠕变仿真参数,包括abaqus蠕变时效本构模型中采用时间硬化率(time hardening) 的幂次法则模式(Power-law model),其蠕变方程如下:
蠕应变
等效应力
t:蠕变时间
通过蠕变曲线,利用蠕变方程进行拟合得到三个蠕变仿真参数n、m和A;本实施例中n=2.794,m=-0.698,A=4.38×10-12
步骤二,将蠕变参数代入已知的机床横梁abaqus模型,运用abaqus有限元分析软件进行机床横梁的蠕变仿真分析。机床横梁的初始abaqus模型为已有的机床横梁残余内应力分布云图,如图2所示,利用abaqus软件自带的时效蠕变本构模型,经过abaqus蠕变仿真得到机床横梁蠕变效应之后的蠕变应力云图如图3。在本步骤中,利用abaqus软件中自带的Visco分析步和实验得出的creep材料属性,对机床横梁进行蠕变分析,得出机床横梁的蠕应力和蠕应变。
步骤三,对机床关重件进行相关蠕变测试实验,对abaqus蠕变仿真进行修正,步骤如下:
1)将机床关重件静置一段时间,静置时间为仿真设置时间的一半;
2)测出机床关重件在这一段时间内测点的几何精度变化;
3)利用abaqus软件对机床关重件在这一段时间的蠕变进行分析;
4)对实验测得的数据和仿真得到的相应的数据进行拟合修正,得到相应的修正公式或者修正系数;
5)将得到的修正系数或者公式用于对步骤二中的仿真结果进行修正。
对abaqus蠕变仿真进行修正具体到本实施例中为:将步骤2)中所测得的横梁导轨面直线度实验数据和步骤3)横梁蠕变仿真导轨面法向位移数据进行整理,利用实验数据对仿真数据进行拟合修正,得到横梁导轨面法向位移数据修正系数或者修正公式。以测试点的实验数据拟合修正仿真数据,再推广到整个模型,将得到的修正系数或者公式用于步骤二中的横梁蠕变仿真导轨面法向位移结果数据的修正。
步骤四,通过修正后的蠕变仿真分析得到由于机床关重件残余应力松弛、机床关重件产生蠕变所造成的几何误差,再将蠕变仿真产生的几何误差和机床关重件的初始几何误差进行叠加,再将叠加后的误差与机床关重件几何精度要求进行对比,判断出机床关重件蠕变后的几何精度是否符合要求。
具体到本实施例中,修正后的蠕变仿真可以得到由于机床横梁残余应力松弛、机床横梁产生蠕变所产生的导轨面直线度误差,再将该蠕变产生的几何误差与横梁初始状态导轨面的直线度误差进行叠加,并将叠加结果与机床横梁导轨面直线度要求进行对比,可以判断出机床横梁蠕变后的导轨面直线度是否要求。
因为初始仿真模型经过处理只存在残余内应力场,所以蠕变前机床横梁本身存在的加工误差需要实验测得。机床横梁精度要求为横梁导轨的直线度误差,经过实际测试,横梁上下导轨面的直线度误差都为0.01mm,所以对蠕变仿真之后的导轨面法线方向上的位移进行离散点的提取。然后将离散点的数据进行曲线拟合如图5,可以得到导轨面的直线度误差,考虑最坏的情况,即蠕变之前的直线度和蠕变之后的直线度在最大误差处叠加,即将两个直线度误差直接相加得到最坏情况出现的最大直线度误差,通过机床横梁导轨面直线度的要求标准0-0.02mm,可以对比判断出机床横梁蠕变之后是否符合几何精度要求,仿真数据结果判断如下表1。
表1仿真数据判断
步骤五,根据机床横梁蠕变后几何精度是否符合要求,判断出机床横梁蠕变之前的残余内应力是否可靠。
如果步骤四中所有指标判断结果都为1,那么认为横梁蠕变后直线度符合标准要求。得出机床横梁直线度符合要求,那么就可以认为机床横梁蠕变之前的残余应力是可靠的。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,但不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (2)

1.一种关于机床关重件残余应力水平的评估方法,其特征在于:包括以下步骤:
步骤一,通过蠕变实验,获取关重件材料在常温25℃下的蠕变仿真参数;
步骤二,将蠕变参数代入已知的机床关重件abaqus模型,运用abaqus有限元分析软件进行机床关重件的蠕变仿真分析;
步骤三,对机床关重件进行相关蠕变测试实验,对abaqus蠕变仿真进行修正;
步骤四,通过修正后的蠕变仿真分析得到由于机床关重件残余应力松弛、机床关重件产生蠕变所造成的几何误差,再将蠕蠕变仿真产生的几何误差和机床关重件的初始几何误差进行叠加。再将叠加后的误差与机床关重件几何精度要求进行对比,判断出机床关重件蠕变后的几何精度是否符合要求;
步骤五,根据机床关重件蠕变后几何精度是否符合要求,判断出机床关重件蠕变之前的残余内应力是否可靠。
2.根据权利要求1所述的一种关于机床关重件残余应力水平的评估方法,其特征在于:在步骤三中对abaqus蠕变仿真进行修正的步骤如下:
1)将机床关重件静置一段时间,静置时间为仿真设置时间的一半;
2)测出机床关重件在这一段时间内测点的几何精度变化;
3)利用abaqus软件对机床关重件在这一段时间的蠕变进行分析;
4)对实验测得的数据和仿真得到的相应的数据进行拟合修正,得到相应的修正公式或者修正系数;
5)将得到的修正系数或者公式用于对步骤二中的仿真结果进行修正。
CN201811539939.XA 2018-12-17 2018-12-17 一种关于机床关重件残余应力水平的评估方法 Pending CN109614724A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811539939.XA CN109614724A (zh) 2018-12-17 2018-12-17 一种关于机床关重件残余应力水平的评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811539939.XA CN109614724A (zh) 2018-12-17 2018-12-17 一种关于机床关重件残余应力水平的评估方法

Publications (1)

Publication Number Publication Date
CN109614724A true CN109614724A (zh) 2019-04-12

Family

ID=66010290

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811539939.XA Pending CN109614724A (zh) 2018-12-17 2018-12-17 一种关于机床关重件残余应力水平的评估方法

Country Status (1)

Country Link
CN (1) CN109614724A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112378686A (zh) * 2020-11-13 2021-02-19 大自然科技股份有限公司 一种床垫支撑稳定性指标的评判方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20010095458A (ko) * 2000-03-30 2001-11-07 홍상복 리드 프레임 소재의 잔류응력 측정장치
CN103745114A (zh) * 2014-01-17 2014-04-23 北京印刷学院 一种钛合金应力松弛数值及回弹计算方法
CN107430637A (zh) * 2015-03-05 2017-12-01 株式会社神户制钢所 残余应力推定方法及残余应力推定装置
CN108304657A (zh) * 2018-02-02 2018-07-20 重庆大学 基于有限元的机床关重件残余应力多工序连续建模仿真方法
CN108536970A (zh) * 2018-04-13 2018-09-14 清华大学深圳研究生院 一种铝合金预拉伸板内部残余应力的测试方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20010095458A (ko) * 2000-03-30 2001-11-07 홍상복 리드 프레임 소재의 잔류응력 측정장치
CN103745114A (zh) * 2014-01-17 2014-04-23 北京印刷学院 一种钛合金应力松弛数值及回弹计算方法
CN107430637A (zh) * 2015-03-05 2017-12-01 株式会社神户制钢所 残余应力推定方法及残余应力推定装置
CN108304657A (zh) * 2018-02-02 2018-07-20 重庆大学 基于有限元的机床关重件残余应力多工序连续建模仿真方法
CN108536970A (zh) * 2018-04-13 2018-09-14 清华大学深圳研究生院 一种铝合金预拉伸板内部残余应力的测试方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘凤岐、汤心颐, 高等教育出版社 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112378686A (zh) * 2020-11-13 2021-02-19 大自然科技股份有限公司 一种床垫支撑稳定性指标的评判方法
CN112378686B (zh) * 2020-11-13 2024-03-12 大自然科技股份有限公司 一种床垫支撑稳定性指标的评判方法

Similar Documents

Publication Publication Date Title
KR100948035B1 (ko) 인장시험과 유한요소법을 이용한 고 변형률에 대한 진변형률-진응력 곡선의 획득 방법 및 이를 이용한 인장 시험기
US20120053858A1 (en) Fatigue Life Estimation Method and System
Ainsworth R5 procedures for assessing structural integrity of components under creep and creep–fatigue conditions
JP2014071053A (ja) 高温部材のクリープ損傷評価方法および損傷評価システム
CN109146402A (zh) 一种软件开发供应商交付质量的评估方法
CN109388878A (zh) 一种综合考虑喷丸强化效果的疲劳寿命预测方法
CN111428361A (zh) 一种适用于多种循环载荷的寿命预测方法
Livieri et al. Fatigue limit evaluation of notches, small cracks and defects: an engineering approach
CN109614724A (zh) 一种关于机床关重件残余应力水平的评估方法
Kawabata et al. Investigation on η and m factors for J integral in SE (B) specimens
Calvín et al. Distribution of the through-thickness effective stress intensity factor range and its influence on fatigue crack growth rate curves
Méndez-González et al. The alpha power Weibull transformation distribution applied to describe the behavior of electronic devices under voltage stress profile
CN108182503B (zh) 一种工业用离心泵的质量评价方法及其设备
Susmikanti et al. Analysis of 3d semi-elliptical crack on reactor pressure vessel wall with load stress and crack ratio
Bomba et al. Geometric Measurements on a CNC Machining Device as an Element of Closed Door Technology
Toteva et al. Methods for selection of measuring instruments
Meggiolaro et al. Automation of the fatigue design under variable amplitude loading using the ViDa software
JP2019082985A (ja) 非線形応力ひずみ解析装置、非線形応力ひずみ解析方法、及び非線形応力ひずみ解析プログラム
US20220382247A1 (en) Estimation method and estimation system
Aibara et al. Analysis of fatigue crack configuration influence on fatigue life
Domozhirov Accumulation of corrosion-fatigue damage by hydro turbine equipment elements
Mirzaee-Sisan et al. Standardisation on Measurement and Interpretation of Residual Stress Data
Priyadumkol et al. A Study of Mechanical Property of Eggshell with Inverse Method and Finite Element Analysis
KR102403748B1 (ko) 경도 측정기 검증 시스템
CN117744932A (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
RJ01 Rejection of invention patent application after publication

Application publication date: 20190412

RJ01 Rejection of invention patent application after publication