本发明涉及振动时效技术领域，特指一种基于应变振型的振动时效激振频率的确定方法。
技术背景
消除构件加工制造过程中的残余应力是机械制造工业领域中的一项重要研究课题。振动时效技术因为具有设备简单、处理时间短、处理效果好等特点而成为备受关注的残余应力消除技术。国内外专家学者已经对振动时效工艺和理论开展了深入的研究，并取得了丰富的理论成果和应用实例。通过国内外专家学者的研究发现，影响振动时效效果的关键工艺参数主要包括：激振频率、激振振幅和激振时间。目前对于激振频率的选择主要采用扫频法，即对构件进行扫频振动，寻找构件振动幅值较大的共振频率，然后在其共振频率或者亚共振频率下进行振动时效处理。然而通过扫频法确定振动时效的激振频率时，通常并不考虑材料的残余应力分布状态以及振型的分布规律，极大地限制了振动时效消除残余应力的效果。针对扫频法确定振动时效激振频率存在的不足，本发明提出一种基于应变振型的振动时效激振频率的确定方法。
技术实现要素：
针对扫频法确定振动时效激振频率存在的不足，本发明提出一种基于应变振型的振动时效激振频率的确定方法。
基于应变振型的振动时效激振频率的确定方法，包括以下步骤：
(1)、获取构件的表面残余应力分布状态；
(2)、推导应变振型与位移振型的关系，建立应变振型的控制方程；
(3)、采用有限元软件建立构件的有限元模型，对构件进行数值模态分析；
(4)、获取构件的各阶应变振型和固有频率；
(5)、在步骤(1)和(4)的基础上，确定振动时效的激振频率f0；
(6)、以(f0-f0×20％)hz为初始激振频率，对构件进行扫频振动，获取振动幅值最大时的频率作为振动时效的激振频率f1。
步骤(2)中建立应变振型的控制方程过程为：以位移振型为基础，推导应变振型与位移振型的关系，建立应变振型的控制方程，为后续的数值模态分析奠定基础。一个典型的有限单元是由节点j，k，m和相应的连接直线来定义的，则该单元内任意一点的位移ui可以表示为
ui＝niai(1)
式中：ni是基于节点位置事先给定的函数，称为形状函数矩阵，ai表示第i个单元所有节点位移组成的节点位移向量。
当已知单元内任意位置的位移，则相应的可以确定出单元内任意位置的应变εi，其表达式为
εi＝siui(2)
式中：si为微分算子矩阵。若系统被划分为n个单元，则共有n个式(2)这样的方程，可以表示为
可以简写为
ε＝su(4)
式中：ε为n个单元的应变向量，u为n个单元的节点位移向量。式(4)建立的应变与位移的关系是在单元坐标系中建立的，考虑到各单元在连接节点上位移是连续变化的，则单元节点位移向量可以通过式(5)的坐标变换转换到总体坐标系中
u＝tug(5)
式中：t为坐标变换矩阵，ug为在总体坐标系中的节点位移向量。将式(5)代入到(4)中，可以得到总体坐标系中的应变方程为
采用有限元模型，在总体坐标系中，系统的无阻尼运动方程为
式中：üg为总体坐标系中的节点加速度向量，fg为载荷向量。令fg＝fmaxejωt，则节点的位移向量为ug＝umaxejωt，代入方程(7)可以得到
(kg-ω2m)umax＝fmax(8)
采用模态叠加法，节点位移的响应可以表示为
式中：
将方程(9)代入到方程(6)中，可以得到
令
步骤(5)中，在步骤(1)和(4)的基础上，确定振动时效的激振频率f0包括以下步骤：
(5.1)、在步骤(1)的基础上，确定构件较大残余应力所在的区域；
(5.2)、在步骤(4)的基础上，确定应变振型较大应变所在的区域；
(5.3)、当应变振型较大应变所在的区域与构件较大残余应力所在的区域一致时，记录下该应变振型所对应的固有频率，并记为f0。
步骤(6)中，以(f0-f0×20％)hz为初始激振频率，对构件进行扫频振动，获取振动幅值最大时的频率作为振动时效的激振频率f1包括以下步骤：
(6.1)、扫频振动的初始激振频率f2设置为(f0-f0×20％)hz，然后以10hz为步长逐步增加扫频振动的频率；分别记录下每一激振频率时作用在构件上的振动幅值，得到振动幅值最大时的频率，并记为f21；
(6.2)、扫频振动的初始激振频率设置为(f21-10)hz，以1hz为步长逐步增加扫频振动的频率；然后重复步骤(6.1)的过程，得到振动幅值最大时的频率，并记为f22，即为振动时效的激振频率f1。
进一步，所述的有限元软件为ansys有限元软件。
本发明的技术构思是：通过本发明提出的方法确定的激振频率是以构件表面的残余应力分布为依据，共振应变函数ψr为准则的，有利于获得较为理想的振动时效效果。
本发明的有益效果如下：
1、传统的扫频振动法确定振动时效的激振频率是以0hz为初始激振频率，而本发明通过采用有限元数值模拟的方式初步确定振动时效的激振频率，然后以此频率为基础，对构件进行扫频振动确定振动时效的激振频率，能够减少扫频振动的时间，提高振动时效的工作效率。
2、传统的扫频振动法确定振动时效的激振频率时并不考虑构件的残余应力分布状态，而本发明提出的基于应变振型的振动时效激振频率的确定方法，是以构件表面的残余应力分布状态为基础，并结合构件的应变振型，这样确定的激振频率能够提高振动时效的效果。
3、构件的应变振型相比于位移振型对构件的局部变化更加敏感，能够反映构件局部位置的能量变化情况，即反映构件局部位置的动应力分布情况，因此以应变振型为基础确定的构件振动时效激振频率能够提高振动时效的效果。
附图说明
图1基于应变振型的振动时效激振频率的确定方法的流程示意图。
图2模态分析构件示意图。
图3激光表面处理工艺的原理图。
图4残余应力测点示意图。
图5振动时效前后构件1#残余应力测试结果:(a)x轴向残余应力；(b)y轴向残余应力。
图6振动时效前后构件2#残余应力测试结果:(a)x轴向残余应力；(b)y轴向残余应力。
图7无孔和有孔构件的低阶和高阶模态位移。
图8无孔和有孔构件的低阶和高阶模态应变。
具体实施方式
参照附图，进一步说明本发明：
基于应变振型的振动时效激振频率的确定方法，包括以下步骤：
(1)、获取构件的表面残余应力分布状态；
具体实施细节为：
实验中选用45#优质碳素结构钢来制备实验构件。实验构件(即数值模态分析的构件)的尺寸为230mm×210mm×6mm，如图2所示。为了使得构件的表面产生较大的拉伸残余应力，实验中采用激光表面处理工艺来产生残余应力，如图3所示。这种工艺是通过激光与构件相互作用区域快速的升温来改变构件的微观组织形态，构件的表面产生热损伤而产生初始拉伸残余应力。激光表面处理工艺在slc-2030d龙门式co2多功能激光加工机上展开，激光器的型号为rofinslabdc030。
构件残余应力测点以及应变花粘贴方向如图4所示。实验中对试样的两个区域进行激光表面处理，每个激光热损伤区域的中心距离试样的边缘均为55mm。实验中选取14个测点用于评价振动前后试样的残余应力，其中测点1#～7#被分为一组，用于测试振动前试样的残余应力，测点1'#～7'#被分为另一组，用于测试振动后试样的残余应力。每一组测点对称的分布在垂直于垂直中心线的直线上，距离垂直中心线的距离分别为0mm、40mm、60mm、80mm，每组测点距离试样边缘的距离均为80mm。以测点4为例，其应变花的粘贴方向如图4所示，其余测点处应变花的粘贴方向与测点4处应变花的粘贴方向一致。通过应变仪测试试样钻孔处的释放应变ε1、ε2、ε3，并将释放应变代入如下的方程就可以求解出试样钻孔处x轴向和y轴向的残余应力σx和σy
式中:e为弹性模量、ν为泊松比；
振动时效处理前构件1#和2#残余应力测试结果如图5和图6所示。构件1#和构件2#振动时效前残余应力的测试结果表明，构件较大的残余应力区域为激光表面烧伤处理中心区域，随着远离激光表面烧伤中心区域距离的增加，残余应力则逐渐降低。
(2)、推导应变振型与位移振型的关系，建立应变振型的控制方程；
(3)、采用有限元软件建立构件的有限元模型，对构件进行数值模态分析；
具体实施细节为：
建立物理模型：数值模态分析的构件就是实验中采用的45#钢板状构件，长230mm，宽210mm、厚6mm，如图2所示。为了研究构件局部结构的变化对其位移振型和应变振型的影响，对两种结构的板状构件进行数值模态分析，其中板状构件1#的中心无通孔，称为无孔构件，另外一板状构件2#的中心有φ6mm的通孔，称为有孔构件。由于数值模态分析的构件三维几何模型相对简单，此处选择直接在ansys经典gui界面中建立构件的三维几何模型。
材料设置：数值模态分析时需要用到45#钢的密度ρ、弹性模量e、泊松比ν，具体数值分别为ρ＝7850kg/m3、e＝200gpa、ν＝0.28。
有限元模型：采用20节点的solid95单元对构件划分网格，建立其三维有限元模型。为了提高求解的精度，在试件的长度、宽度以及厚度三个方向设置单元划分的数目分别为55、45以及3，如图2所示。然后采用扫掠(volumesweep)的方式对构件划分网格，可以得到规则的网格，也可以提高求解的精度。另外，solid95单元具有中间节点，这将使得求解的结果更为精确。
求解方法：采用ansys模态分析模块中的blocklanczos方法(分块兰索斯法)求解构件的位移模态振型和应变模态振型。blocklanczos法特征值求解器是ansys模态分析模块中的默认求解器，采用的是lanczos算法。这种算法不仅可以保证求解的精确性，更能提高求解的速度。
数值模态分析的结果为：无孔构件低阶模态振型(模态频率为572hz)；有孔构件低阶模态振型(模态频率为571hz)。无孔构件高阶模态振型(模态频率为3231hz)；有孔构件高阶模态振型(模态频率为3229hz)。无孔和有孔构件的低阶和高阶模态位移沿着图2所示的映射路径的分布规律见图7所示；无孔和有孔构件的低阶和高阶模态应变沿着图2所示的映射路径的分布规律见图8所示。从图7可见，无孔和有孔构件的低阶和高阶模态位移值基本是重叠在一起的，即使在钻孔处无孔和有孔构件的低阶和高阶模态位移值也没有明显的变化，表明位移振型对结构的局部变化并不敏感。从图8可见，无孔和有孔构件的低阶和高阶模态应变值在钻孔以外的地方基本是重叠在一起的，但是在钻孔处无孔和有孔构件的低阶和高阶模态应变值却发生了明显的突变，表明应变振型对结构的局部变化是敏感的，并且高阶模态应变值在钻孔处的突变更加明显，这也说明了为什么从应变振型的角度确定振动时效的激振频率。
(4)、获取构件的各阶应变振型和固有频率；
(5)、在步骤(1)和(4)的基础上，确定振动时效的激振频率f0；
(6)、以(f0-f0×20％)hz为初始激振频率，对构件进行扫频振动，获取振动幅值最大时的频率作为振动时效的激振频率f1。
步骤(2)中建立应变振型的控制方程过程为：以位移振型为基础，推导应变振型与位移振型的关系，建立应变振型的控制方程，为后续的数值模态分析奠定基础。一个典型的有限单元是由节点j，k，m和相应的连接直线来定义的，则该单元内任意一点的位移ui可以表示为
ui＝niai(1)
式中：ni是基于节点位置事先给定的函数，称为形状函数矩阵，ai表示第i个单元所有节点位移组成的节点位移向量。
当已知单元内任意位置的位移，则相应的可以确定出单元内任意位置的应变εi，其表达式为
εi＝siui(2)
式中：si为微分算子矩阵。若系统被划分为n个单元，则共有n个式(2)这样的方程，可以表示为
可以简写为
ε＝su(4)
式中：ε为n个单元的应变向量，u为n个单元的节点位移向量。式(4)建立的应变与位移的关系是在单元坐标系中建立的，考虑到各单元在连接节点上位移是连续变化的，则单元节点位移向量可以通过式(5)的坐标变换转换到总体坐标系中
u＝tug(5)
式中：t为坐标变换矩阵，ug为在总体坐标系中的节点位移向量。将式(5)代入到(4)中，可以得到总体坐标系中的应变方程为
采用有限元模型，在总体坐标系中，系统的无阻尼运动方程为
式中：üg为总体坐标系中的节点加速度向量，fg为载荷向量。令fg＝fmaxejωt，则节点的位移向量为ug＝umaxejωt，代入方程(7)可以得到
(kg-ω2m)umax＝fmax(8)
采用模态叠加法，节点位移的响应可以表示为
式中：
将方程(9)代入到方程(6)中，可以得到
令
步骤(5)中，在步骤(1)和(4)的基础上，确定振动时效的激振频率f0包括以下步骤：
(5.1)、在步骤(1)的基础上，确定构件较大残余应力所在的区域；
(5.2)、在步骤(4)的基础上，确定应变振型较大应变所在的区域；
(5.3)、当应变振型较大应变所在的区域与构件较大残余应力所在的区域一致时，记录下该应变振型所对应的固有频率，并记为f0。
具体实施细节为：
综合图5和图6振动时效处理前的构件残余应力测试结果以及图8的模态应变分析结果可知，高阶应变振型较大应变所在的区域能够覆盖构件较大残余应力所在的区域，因此振动时效的激振频率为3231hz，即构件2#通过数值模态分析的方法确定的振动时效的激振频率为3231hz。
为了验证本发明提出的方法能够提高振动时效效果的有效性，同时选取振动时效的激振频率为572hz，即构件1#通过数值模态分析的方法确定的振动时效的激振频率为572hz。
步骤(6)中，以(f0-f0×20％)hz为初始激振频率，对构件进行扫频振动，获取振动幅值最大时的频率作为振动时效的激振频率f1包括以下步骤：
(6.1)、扫频振动的初始激振频率f2设置为(f0-f0×20％)hz，然后以10hz为步长逐步增加扫频振动的频率；分别记录下每一激振频率时作用在构件上的振动幅值，得到振动幅值最大时的频率，并记为f21；
(6.2)、扫频振动的初始激振频率设置为(f21-10)hz，以1hz为步长逐步增加扫频振动的频率；然后重复步骤(6.1)的过程，得到振动幅值最大时的频率，并记为f22，即为振动时效的激振频率f1。
具体实施细节：
对于构件1#设置初始激振频率为(572-572×20％)hz，然后按照步骤(6.1)和步骤(6.2)的方法寻找到振动幅值最大的频率为717hz。
对于构件2#设置初始激振频率为(3231-3231×20％)，然后按照步骤(6.1)和步骤(6.2)的方法寻找到振动幅值最大的频率为3311hz。
通过本发明提出的方法确定的振动时效激振频率对构件1#和构件2#进行振动时效处理。构件1#和构件2#振动时效处理后的残余应力测试结果如图5和图6所示。构件1#和构件2#振动时效处理后残余应力峰值消除率均较大，而对于残余应力较小的地方(远离激光表面烧伤中心区域)，振动时效处理后残余应力的消除率相对较小。由图5可见，构件1#经过振动时效处理后其x轴向和y轴向峰值残余应力的释放量分别为106.8mpa和131.8mpa；由图6可见，构件2#经过振动时效处理后其x轴向和y轴向峰值残余应力的释放量分别为183.6mpa和226.1mpa，表明根据残余应力分布状态和应变振型分布规律，合理的选择振动时效的激振频率能够获得较为理想的振动时效消除残余应力效果。由图5和图6可见，构件2#具有更好的振后残余应力均化分布效果，表明根据残余应力分布状态和应变振型分布规律，合理的选择振动时效的激振频率，可以获得较为理想的振动时效均化残余应力效果。综上所述，对于构件进行振动时效处理时，需要根据构件表面的残余应力分布状态和应变振型分布规律，合理选取振动时效的激振频率，能够获得较为理想的振动时效消除和均化残余应力效果。
进一步，所述的有限元软件为ansys有限元软件。
本说明书实施例所述的内容仅仅是对发明构思的实现形式的列举，本发明的保护范围不应当被视为仅限于实施例所陈述的具体形式，本发明的保护范围也及于本领域技术人员根据本发明构思所能够想到的等同技术手段。