一种基于固定界面法的铣削时变信息更新方法技术

本发明专利技术属于智能制造与加工技术领域，公开一种基于固定界面法的铣削时变信息更新方法，包括以下步骤：根据待加工工件的特征，将整体工件分为恒定工件与移除材料工件两个子结构部分，分别对两个子结构模型使用固定界面法进行分析，并按照界面协调条件，耦合模态坐标空间中的两个子结构模型，获得薄壁构件整体的自由度降阶模型与模态特性

【技术实现步骤摘要】
一种基于固定界面法的铣削时变信息更新方法


[0001]本专利技术属于智能制造与加工
，涉及到铣削过程中薄壁构件的颤振预测与特征映射问题，具体为一种基于固定界面法的铣削时变信息更新方法
。

技术介绍

[0002]大型薄壁结构，如框架，肋条和翼梁等被广泛用于飞机，以提高结构效率和轻量化
。
这些结构主要由多个薄壁构件组成，其壁高和长度与厚度的纵横比通常可分别达到
10∶1
～
50∶1
，甚至达到
50∶1
～
250∶1。
考虑到多个型腔的特点，通常使用铣削工艺来实现最终形状
。
但是，由于结构的刚性低，容易发生颤振，导致表面质量差和加工不准确，进而造成一系列工程问题
。
[0003]从理论上讲，颤振是由于切屑厚度在切削力作用下再生而产生的一种自激振动，它与刀具主轴系统和加工过程中工件动力学系统有着紧密的联系
。
然而，从实验角度出发，对加工过程中工件的实验模态进行重复的离线测量
、
分析是不切实际的
。
有限元方法作为一种高效的现代计算技术，可以很好地表示铣削过程中遇到的工件和夹具布局的复杂几何形状
。
然而，由于材料去除的影响，加工过程中工件的有限元模型必须在每个刀具位置通过重新网格进行更新，并且必须在精加工过程中对全阶有限元模型进行模态重新分析
。
为了规避这一难题，相关领域内学者开发了结构动力学修改技术来计算加工过程中工件动力学，与此同时，根据研究对象不同，相关学者对上述技术不断进行优化
。
通过在加工工件上沿刀具路径以相反的顺序添加去除的材料来计算频域中的
IPW
动态
(
如文献
Budak E
，
Altintas Y.Analytical prediction of chatter stability in milling—Part II
：
Application of the general formulation to common milling systems[J].Journal of Dynamic Systems Measurement and Control
‑
transactions of the Asme
，
1998
，
120
：
31—36.)。
提出了一种时域更新策略，即从初始工件中减去去除的材料
(
如文献
Tuysuz O
，
Altintas Y.Time
‑
domain modeling 0fvarying dynamic characteristics in thin
‑
wall machining using perturbation and reduced
‑
order substructuring methods[J].Journal of Manufacturing Science and Engineering
‑
Transactions of the Asme
，
2018
，
140
，
011015.)。
使用迭代改进的还原系统技术减少去除材料和初始工件的模型阶数
(
如文献
Yang Y,Zhang WH
，
et al.An efficient decomposition
‑
condensation method for chatter prediction in milling large
‑
scale thin
‑
walled structures[J].Mechanical Systems And Signal Processing
，
2019
，
121
：
58
‑
76.)。
他们通过添加与去除的材料具有相反动力学的虚构子结构，消除了去除的材料在频域中对初始工件动力学的影响；后来，他们继续将这种更新策略扩展到使用矩阵扰动技术对薄壁工件进行时域建模
。
上述这些策略的一个共同优点是，只需要对整个初始工件或加工工件的有限元模型进行一次耗时的分析
。
[0004]在材料去除过程中，整个工件只有一小部分区域在变化
。
从结构动力修饰技术的角度来看，它属于局部修饰
。
然而，上述相关研究将这一过程视为一个全局性的修改问题，
在修改过程中涉及初始工件或加工工件的整个领域
。
也就是说，无论模型是否真的减少，这些策略在修改过程中都比局部修改策略使用更多的自由度，即它们需要更多的计算时间和内存
。
为此，亟需研究一种高效
、
稳定的铣削时变信息更新技术，有效解决薄壁零件铣削加工过程中的时变信息运算难题
。

技术实现思路

[0005]为有效解决现有技术存在的问题，本专利技术提供了一种固定界面法的铣削时变信息更新技术，通过机理模型高效获取铣削过程中工件的模态特性，结合数据模型完成其时变信息的特征映射
。
根据待加工工件的特征，将整体工件分为恒定工件与移除材料工件两个子结构部分，分别对两个子结构模型使用固定界面法进行分析，并按照界面协调条件，耦合模态坐标空间中的两个子结构模型，获得薄壁构件整体的自由度降阶模型与模态特性
。
同时引入结构动力学修改方法，揭示铣削过程中材料去除的变化规律，实现对包含移除材料子结构的去除量更新，并结合数据模型，完成铣削时变信息的表征
。
以时变信息的更新精度与速度为综合效率指标，对比传统有限元方法，本专利技术大幅度提升了整体运算过程的效率
。
[0006]本专利技术的技术方案：
[0007]一种基于固定界面法的铣削时变信息更新方法，包括以下步骤：
[0008]步骤1：铣削动力学模型的建立；
[0009]依据铣削过程中刀具与薄壁构件的接触特征，沿刀具轴向变化情况，对两者之间的多点接触区域建立动力学模型，同时引入模态变换方法，对刀具和工件的动力学模型由物理空间转换为模态空间，并将动力学方程表示为矩阵形式，以方便系统方程的解耦与运算
。
[0010]步骤2：子结构模型的划分；
[0011]分析薄壁构件的形态与尺寸，对其建立有限元模型，同时根据待加工工件的特征，将整体工件分为恒定工件与移除材料工件两个子结构部分
。
[0012]步骤3：
α
子结构的分析
[0013]按照工件的内部条件和边界条件，重新本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.
一种基于固定界面法的铣削时变信息更新方法，其特征在于，包括以下步骤：步骤1：铣削动力学模型的建立依据铣削过程中刀具与薄壁构件的接触特征，沿刀具轴向变化情况，对刀具与薄壁构件之间的多点接触区域建立动力学模型；同时引入模态变换方法，对刀具和薄壁构件的动力学模型由物理空间转换为模态空间，并将动力学方程表示为矩阵形式，以方便系统方程的解耦与运算；刀具与薄壁构件动力学模型的矩阵形式如下所示：其中，下标
f
代表
t
或
w
，分别对应着刀具或薄壁构件；
X
f
(t)、
和分别代表物理空间中
3n
f
*1
维度的位移
、
速度以及加速度矢量；
M
f
、C
f
和
K
f
分别代表物理空间中
3n
f
*3n
f
维度的质量
、
阻尼和刚度矩阵；
F
f
(t)
代表
3n
f
*1
维度的刀具或薄壁构件切削力矢量，并且存在关系
F
w
(t)
＝
‑
F
t
(t)
＝
[F1(t),...,F
j
(t),...,F
q
(t)]
T
；引入模态变换方法，对刀具和薄壁构件的动力学模型由物理空间转换为模态空间，并将动力学方程表示为矩阵形式，该铣削系统在状态空间中的动力学方程如下所示：其中，和分别代表模态空间中
r
f
*r
f
维度的质量
、
阻尼和刚度矩阵；和
P
f
(t)
分别代表加速度
、
速度和位移矩阵，为模态矩阵转置；
F
f
(t)
为力矩阵，且满足下列条件：其中，
I
f
代表单位矩阵，
ζ
f
代表阻尼比矩阵，
ω
f
代表固有频率矩阵，均为
r
f
*r
f
维度的对角矩阵；对式
(3)
整理获得模态空间中刀具
‑
主轴系统和薄壁构件的动力学方程，具体如下所示：其中，和
P
t,w
(t)
分别代表加速度
、
速度和位移矩阵，为模态矩阵转置；参数具体表示如下：其中，
U
t,w
是刀具
‑
主轴系统和薄壁构件之间接触区域节点的质量归一化模态振型矩阵，具体如下所示：
U
t,w
＝
[
‑
U
t
,U
w
]
ꢀꢀꢀꢀ
(6)
通过求解上述公式的系统动力学参数：物体固有属性
、
状态变量
、
传递函数响应特性，便求解模态空间内铣削系统的动力学模型；
步骤2：子结构模型的划分分析薄壁构件的形态与尺寸，对其建立有限元模型；同时根据待加工工件的特征，将整体工件分为恒定工件与移除材料工件两个子结构部分，分别为
α
子结构和
β
子结构；步骤3：
α
子结构的分析当薄壁构件被分为
α
子结构与
β
子结构时，按照两物体的接触情况，定义两物体接触面为边界区域，其余部分各位两物体的内部区域；因此，重新排列
α
子结构各个位置的节点；其中，边界节点位于
α
子结构与
β
子结构的连接面处，剩余的节点为
α
子结构的内部节点；具体如下所示：其中，下标
i
与
b
分别代表
α
子结构的内部自由度和边界自由度，即
α
M
ii
代表内部
‑
内部节点处相关自由度参数值；
α
M
ib
代表内部
‑
边界节点处相关自由度参数值；
α
M
bi
代表边界
‑
内部节点处相关自由度参数值；
α
M
bb
代表边界
‑
边界节点处相关自由度参数值；刚度矩阵
α
K
的参数值分布规则与
α
M
相同，此处不再赘述；代表内部节点处相关自由度的加速度，代表边界节点处相关自由度的加速度，
α
X
i
代表内部节点处相关自由度的速度，
α
X
b
代表边界节点处相关自由度的速度，由于铣削过程中的铣削力仅作用于边界节点，因此
{
α
F}
＝
{0 α
F
b
}
T
；同时引入固定界面模态综合法，通过处理正常保留模态和边界约束模态，构建子结构的模态转换矩阵，实现该部分的模型降阶；其中，在位移变换中使用固定边界下的正常保留模态和边界约束模态，获得位移变换矩阵：其中，
[
Φ
ik
]
...

【专利技术属性】
技术研发人员：李宏坤，王朝东，刘淑杰，李智，
申请(专利权)人：大连理工大学，
类型：发明
国别省市：