CN104656562A - 一种柔性结构件铣削加工的工艺参数优化方法 - Google Patents
一种柔性结构件铣削加工的工艺参数优化方法 Download PDFInfo
- Publication number
- CN104656562A CN104656562A CN201410848655.4A CN201410848655A CN104656562A CN 104656562 A CN104656562 A CN 104656562A CN 201410848655 A CN201410848655 A CN 201410848655A CN 104656562 A CN104656562 A CN 104656562A
- Authority
- CN
- China
- Prior art keywords
- flexible
- workpiece
- cutting
- displacement
- flexible workpiece
- 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
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
- G05B19/4097—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by using design data to control NC machines, e.g. CAD/CAM
Landscapes
- Engineering & Computer Science (AREA)
- Human Computer Interaction (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Automatic Control Of Machine Tools (AREA)
- Numerical Control (AREA)
Abstract
本发明提供了一种柔性结构件铣削加工的工艺参数优化方法，其能在保证柔性结构件加工稳定性的同时保证加工精度。将工件安装于柔性工件加工实验装置形成柔性工件系统，然后按以下步骤进行工艺参数优化：（1）获取柔性工件系统的模态参数；（2）根据柔性工件系统的模态参数，进行切削动力学建模；（3）利用时域数值方法获取铣刀切削柔性工件系统的稳定域边界；（4）采用柔性工件加工实验装置进行切削试验对切削动力学建模进行验证；（5）以机床主轴切削材料去除率最大化和工件振动最小化为优化目标，以稳定域边界为约束条件，建立柔性结构件铣削参数优化模型，通过非线性优化算法来获得优化后的机床主轴转速和切深。
Description
技术领域
本发明涉及计算机数控加工技术领域，尤其是涉及数控铣削加工工艺参数优化领域，具体为一种柔性结构件铣削加工的工艺参数优化方法。
背景技术
在航空、航天以及能源等领域，柔性薄壁结构零件的应用非常广泛，零件的加工精度和加工效率很大程度上依赖于加工参数的合理选择。为了达到柔性结构零件高效精密的加工要求，不仅需要确保加工过程稳定性，即不发生颤振，还需要保证零件的加工精度和低成本。
现有的柔性结构件加工工艺参数优化方法主要基于常数边界的加工动力学模型（Altintas Y, Budak E. Analytical Prediction of Stability Lobes in Milling. CIRP Annals-Manufacturing Technology. 1995;44(1):357-62.），在动力学建模过程中将铣刀刀刃切入切出角当做恒定的常数，即不考虑切削振动对其影响，利用基于摄动假设将切厚动态变化量导致的切削力瞬态分量看做是额定激励的高阶小量，从而在稳定性分析中忽略额定激励导致的切削力稳态量，以切削力瞬态分量的稳定性来判断切削过程是否稳定。以该建模方法获取的稳定性以及基于此稳定性进行的参数优化在刚性大的工件切削过程有较好的预报结果，但是在柔性结构加工过程中往往误差较大，究其原因是由于柔性结构零件由于其结构刚度较低，在切削过程中振动较大，振动位移大大改变实际刀刃的切入切出角，从而无法对真实的工况进行准确地建模，因此获得的稳定性以及工艺参数不是真实的优化值，可能导致高幅值振动乃至颤振，难以实现柔性结构零件的高效、精密加工。
发明内容
针对上述问题，本发明提供了一种柔性结构件铣削加工的工艺参数优化方法，其能在保证柔性结构件加工稳定性的同时保证加工精度，从而实现柔性结构件的精密、高效加工。
一种柔性结构件铣削加工的工艺参数优化方法，其特征在于：首先搭建一套柔性工件加工实验装置，将工件安装于所述柔性工件加工实验装置从而形成一个柔性工件系统，然后按以下步骤进行工艺参数优化：
（1）获取所述柔性工件系统的模态参数，所述模态参数包括模态质量、模态阻尼和模态频率；
（2）根据所述柔性工件系统的模态参数，进行切削动力学建模，建立刀具几何参数、所述柔性工件系统的模态参数、机床主轴转速、工件的振动位移对瞬时切厚、切入切出角以及切削力的函数关系；
（3）利用时域数值方法获取铣刀所述切削柔性工件系统的稳定域边界：根据所述步骤（2）确定的函数关系，基于时域数值仿真和数值积分法来综合求解，并以稳态振动位移与设定值比较为判稳条件，求解铣削稳定域的边界；
（4）采用所述柔性工件加工实验装置进行切削试验对所述切削动力学建模进行实验验证；
（5）以机床主轴切削材料去除率最大化和工件振动最小化为优化目标，以所述稳定域边界为约束条件，建立柔性结构件铣削参数优化模型，通过非线性优化算法来获得优化后的机床主轴转速和切深。
其进一步特征在于：
所述柔性工件加工实验装置包括一双板支承的柔性铰链机构以及与信号采集工控机相连的非接触加工位移传感系统；
所述获取所述柔性工件系统的模态参数的具体过程为，以脉冲实验力作为激励输入信号，以粘贴在工件上的加速度信号为振动输出，通过传统函数计算获取所述柔性工件的所述模态参数；
所述步骤（3）的具体过程为，采用时域数值方法对所述柔性工件系统响应进行输出仿真，首先将整个输出时间段进行离散化，在单个时间步内进行计算和位移、切入切出边界的更新，将输出的稳态位移和设定值作对比，当仿真输出值大于设定值时认为所述柔性工件系统不稳定，当仿真输出值小于设定值时认为所述柔性工件系统稳定，以此来获得整个转速-切深参数平面内的稳定域；
所述单个时间步进行计算和位移、切入切出边界的更新流程为：依次进行切厚计算、切削力计算、振动位移计算、切入切出边界计算以及开关函数计算，所述开关函数计算结果返回至所述切削计算进行循环迭代计算，所述每一次振动位移计算输出后对位移、切入切出边界进行更新；
所述步骤（4）的具体过程为，利用激光位移传感器对所述柔性工件的位移信息和主轴转速进行同步采集，获取真实转速下的工件振动位移；对所述工件位移采用相图法进行判稳分析，当位移相图为点时系统判为稳定，当位移相图为极限环时系统判为不稳定。
由于本发明中考虑了柔性结构零件切削振动位移对切入切出角变化的影响，更加符合真实的工况，提高了柔性结构件加工稳定性预报的准确性；本发明以工件的振动位移幅值的大小来判稳，符合实际工况可操作性条件，对于切削颤振实时监测具有很大的便捷性。利用本发明所提的方法，不仅可保证柔性结构零件加工过程的稳定性，还能保证零件加工过程中具有较小的绝对振幅，即具有较高的加工精度。在同时满足加工稳定性和加工精度的参数范围内对主轴转速和工件振动最小化进行优化选取，可保证柔性结构零件的无颤振高效精密加工。
附图说明
图1为本发明中柔性工件加工实验装置示意图；
图2-1为本发明实施例中柔性工件系统的频响函数图；
图2-2为本发明实施例中刀具系统的频响函数图；
图3为本发明实施例中切入角随工件位移振动的变化示意图；
图4为本发明实施例的单个时间步的时域计算输出仿真流程示意图；
图5为本发明实施例的柔性结构件铣削稳定域与实测对比图；
图6-1为本发明实施例中失稳点的位移和相图；
图6-2为本发明实施例中稳定点参数点处的位移和相图。
具体实施方式
以下结合附图和实施例对本发明的接受方案作进一步的详细描述，但本实施例并不用于限制本发明。
本发明实施例采用的柔性工件加工实验装置如图1所示，将工件1安装在双板支承的铰链柔性实验装置2顶部，然后将整个柔性工件-支承系统安装在机床工作平台4上，其中柔性铰链支承板的结构参数为：板厚度5mm，高、宽分别为100mm。选取超细晶粒硬质合金两齿平底立铣刀5为加工刀具，铣刀直径8mm，螺旋角30°，铣刀悬长30mm。采用激光位移传感器对柔性工件系统和转速进行同步采集，获取真实转速下的工件振动位移。
按照本发明的方法进行柔性结构工件系统切削加工参数优化，其步骤为：
1. 将加速度传感器粘贴在工件上，用脉冲锤激励工件系统，将激励力信号和加速度信号进行同步采集，获取柔性工件系统的传递函数；采用相同的实验模态测试过程对刀具系统进行频响函数的采集和刀尖模态参数的拟合。工件系统的频响函数以及刀具系统的频响函数如图2-1、图2-2所示。由图2-1、2-2可知，刀具系统的刚性比工件系统的刚性高两个数量级，因此刀具被认为是刚性的，采用有理分式多项式拟合算法分别获取工件模态参数为：工件想模态质量为8kg，模态阻尼为3060N·s/m，模态刚度为165820000 N/m；
2. 建立柔性工件铣削动力学方程，即刀具几何参数、工件系统的模态参数、机床主轴转速、工件的振动位移对瞬时切厚、切入切出角以及切削力的函数关系：
在上述模型中，y为振动位移，f z 为名义每齿进给，τ = 60 / (N·n)为单个刀齿的切削周期，n为转速，N为铣刀齿数。g j 为开关函数，由切入切出角来决定其取值：
可以看出，切削过程中工件的大幅振动，将影响切入切出角值，如图3所示，从而影响开关函数和切削力，使得上述系统不再具有周期性，因此对系统的输出响应将采用时域数值法进行求解。
3. 采用数值方法对系统响应进行输出仿真：首先将整个输出时间段进行离散化，在每个小时间段内进行计算和位移、切入切出边界的更新，单个时间步的计算流程如图4所示，将输出稳态位移和设定值做对比，当仿真输出值大于设定值，认为系统不稳定；当仿真输出值小于设定值时，认为系统稳定，以此来获得整个转速-切深参数平面内的稳定域。
4. 采用柔性铰链支承系统实验装置进行切削试验对上述模型进行实验验证。采用激光位移传感器3对柔性工件系统和转速进行同步采集，获取真实转速下的工件振动位移。对工件位移采用相图法进行判稳分析，当位移相图为点时，系统是稳定；当位移相图为极限环时，系统是不稳定的。本实施例中对转速为2000-7800转/分钟，切深为0-2mm内的参数平面域内的预报稳定边界进行部分参数点测试，并与传统模型预报的稳定域边界进行对比，如图5所示，图5中的不可行区域包含了颤振失稳参数点和高振幅参数点；可以看出，本发明中所提新模型预报的稳定域边界与实测数据吻合地很好，相对传统模型预报的稳定域边界来说，不仅排除了不稳定的点，而且能够保证切削振幅保持在一定的范围内，即保证了柔性结构零件的加工精度。
5. 优化材料去除率和加工精度，在动力学方程中，可由
求解出振动位移y和某转速下的切深极限a p，振动位移直接影响加工精度，建立材料去除率和加工精度的优化模型
min{-a p·n·a r ·N·f z }
s.t. y≤y 0, l(m,c,k,N,f z ,Kt,Kr, y 0)≤ 1.05.
在上述优化模型中，优化目标为机床切削用量最大化，以稳定域边界和工件设定振幅为约束条件，该条件给出了机床主轴转速和切深极限的隐式函数关系。通过非线性优化算法可求解出优化后的主轴转速和切削深度，优化模型求解得到的结果如图6所示，转速和切深参数分别为7766rpm和1.5mm，利用该参数点加工测得的其振幅及相图曲线是稳定的。选取传统稳定域分析稳定域最高的区域内的参数点4205rpm，1.5mm进行对比分析，由实验结果表面该参数点的振幅很大，相图分析表明其发生了不稳定极限环。
采用本方面获得的加工参数，在薄壁件加工过程中不仅能保证加工稳定，还能保证加工精度。
Claims (6)
1.一种柔性结构件铣削加工的工艺参数优化方法，其特征在于：首先搭建一套柔性工件加工实验装置，将工件安装于所述柔性工件加工实验装置从而形成一个柔性工件系统，然后按以下步骤进行工艺参数优化：
（1）获取所述柔性工件系统的模态参数，所述模态参数包括模态质量、模态阻尼和模态频率；
（2）根据所述柔性工件系统的模态参数，进行切削动力学建模，建立刀具几何参数、所述柔性工件系统的模态参数、机床主轴转速、工件的振动位移对刀刃瞬时切厚、切入切出角以及切削力的函数关系；
（3）利用时域数值方法获取铣刀所述切削柔性工件系统的稳定域边界：根据所述步骤（2）确定的函数关系，基于时域数值仿真和数值积分法来综合求解，并以稳态振动位移与设定值比较为判稳条件，求解铣削稳定域的边界；
（4）采用所述柔性工件加工实验装置进行切削试验对所述切削动力学建模进行实验验证；
（5）以机床主轴切削材料去除率最大化和工件振动最小化为优化目标，以所述稳定域边界为约束条件，建立柔性结构件铣削参数优化模型，通过非线性优化算法来获得优化后的机床主轴转速和切深。
2.根据权利要求1所述的一种柔性结构件铣削加工的工艺参数优化方法，其特征在于：所述柔性工件加工实验装置包括一双板支承的柔性铰链机构以及与信号采集工控机相连的非接触加工位移传感系统。
3.根据权利要求2所述的一种柔性结构件铣削加工的工艺参数优化方法，其特征在于：所述获取所述柔性工件系统的模态参数的具体过程为，以脉冲实验力作为激励输入信号，以粘贴在工件上的加速度信号为振动输出，通过传统函数计算获取所述柔性工件的所述模态参数。
4.根据权利要求3所述的一种柔性结构件铣削加工的工艺参数优化方法，其特征在于：所述步骤（3）的具体过程为，采用时域数值方法对所述柔性工件系统响应进行输出仿真，首先将整个输出时间段进行离散化，在单个时间步内进行计算和位移、切入切出边界的更新，将输出的稳态位移和设定值作对比，当仿真输出值大于设定值时认为所述柔性工件系统不稳定，当仿真输出值小于设定值时认为所述柔性工件系统稳定，以此来获得整个转速-切深参数平面内的稳定域。
5.根据权利要求4所述的一种柔性结构件铣削加工的工艺参数优化方法，其特征在于：所述单个时间步进行计算和位移、切入切出边界的更新流程为：依次进行切厚计算、切削力计算、振动位移计算、切入切出边界计算以及开关函数计算，所述开关函数计算结果返回至所述切削计算进行循环迭代计算，所述每一次振动位移计算输出后对位移、切入切出边界进行更新。
6.根据权利要求3所述的一种柔性结构件铣削加工的工艺参数优化方法，其特征在于：所述步骤（4）的具体过程为，利用激光位移传感器对所述柔性工件系统的位移信息和主轴转速进行同步采集，获取真实转速下的工件振动位移；对所述工件位移采用相图法进行判稳分析，当位移相图为点时系统判为稳定，当位移相图为极限环时系统判为不稳定。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410848655.4A CN104656562B (zh) | 2014-12-31 | 2014-12-31 | 一种柔性结构件铣削加工的工艺参数优化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410848655.4A CN104656562B (zh) | 2014-12-31 | 2014-12-31 | 一种柔性结构件铣削加工的工艺参数优化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104656562A true CN104656562A (zh) | 2015-05-27 |
| CN104656562B CN104656562B (zh) | 2017-07-07 |
Family
ID=53247826
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410848655.4A Active CN104656562B (zh) | 2014-12-31 | 2014-12-31 | 一种柔性结构件铣削加工的工艺参数优化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104656562B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110007644A (zh) * | 2019-03-07 | 2019-07-12 | 中南大学 | 一种加工综合误差建模方法 |
| CN110377940A (zh) * | 2019-06-10 | 2019-10-25 | 大连理工大学 | 一种获取cfrp工件铣削振动幅度的方法 |
| CN111032281A (zh) * | 2017-08-07 | 2020-04-17 | 弗兰茨·翰默机械制造两合公司 | 在机加工中心创建数字孪生 |
| CN112783138A (zh) * | 2020-12-30 | 2021-05-11 | 上海交通大学 | 生产线设备加工稳定性智能监控和异常诊断方法及装置 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103753353A (zh) * | 2013-12-31 | 2014-04-30 | 无锡透平叶片有限公司 | 一种快速测定铣刀偏心的非接触式激光测量方法 |
-
2014
- 2014-12-31 CN CN201410848655.4A patent/CN104656562B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103753353A (zh) * | 2013-12-31 | 2014-04-30 | 无锡透平叶片有限公司 | 一种快速测定铣刀偏心的非接触式激光测量方法 |
Non-Patent Citations (4)
| Title |
|---|
| P. PALANISAMY等: "Optimization of machining parameters using genetic algorithm", 《INT J ADV MANUF TECHNOL》 * |
| YUSUF ALTINTAS等: "Analytical Prediction of Stability Lobes in Milling", 《CIRP ANNALS一MANUFACTURING TECHNOLOGY》 * |
| 张小俭: "柔性结构铣削时滞工艺系统的稳定性理论与实验研究", 《中国博士学位论文全文数据库工程科技Ⅰ辑》 * |
| 张小明: "五轴数控铣削加工几何-动力学建模与工艺参数优化", 《中国博士学位论文全文数据库工程科技Ⅰ辑》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111032281A (zh) * | 2017-08-07 | 2020-04-17 | 弗兰茨·翰默机械制造两合公司 | 在机加工中心创建数字孪生 |
| CN111032281B (zh) * | 2017-08-07 | 2022-11-25 | 弗兰茨·翰默机械制造两合公司 | 在机加工中心创建数字孪生 |
| CN110007644A (zh) * | 2019-03-07 | 2019-07-12 | 中南大学 | 一种加工综合误差建模方法 |
| CN110007644B (zh) * | 2019-03-07 | 2022-02-11 | 中南大学 | 一种加工综合误差建模方法 |
| CN110377940A (zh) * | 2019-06-10 | 2019-10-25 | 大连理工大学 | 一种获取cfrp工件铣削振动幅度的方法 |
| CN110377940B (zh) * | 2019-06-10 | 2020-11-13 | 大连理工大学 | 一种获取cfrp工件铣削振动幅度的方法 |
| CN112783138A (zh) * | 2020-12-30 | 2021-05-11 | 上海交通大学 | 生产线设备加工稳定性智能监控和异常诊断方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104656562B (zh) | 2017-07-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Shi et al. | An improved cutting power-based model for evaluating total energy consumption in general end milling process | |
| CN110270705B (zh) | 考虑刀具跳动影响的微细铣削加工切削力仿真预测方法 | |
| Shi et al. | A novel energy consumption model for milling process considering tool wear progression | |
| Huang et al. | Mechanistic modeling of process damping in peripheral milling | |
| CN103567815B (zh) | 基于铣削小孔的数控机床切削热误差测试和评价方法 | |
| Wan et al. | Numerical prediction of static form errors in peripheral milling of thin-walled workpieces with irregular meshes | |
| CN102929210B (zh) | 基于特征的数控加工过程控制和优化系统及方法 | |
| CN107457609B (zh) | 基于刚度变化的铣削颤振抑制方法及铣削颤振优化系统 | |
| Costes et al. | Surface roughness prediction in milling based on tool displacements | |
| CN104656562A (zh) | 一种柔性结构件铣削加工的工艺参数优化方法 | |
| CN102601434B (zh) | 一种整体叶轮开槽插铣加工优化方法 | |
| CN101905340B (zh) | 一种高速铣削稳定性快速判定方法 | |
| CN106021796B (zh) | 一种铬钢叶片型面加工用球头铣刀的剩余寿命预测方法 | |
| CN103676782B (zh) | 数控铣床加工过程中能量效率在线检测方法 | |
| CN105414616A (zh) | 螺旋铣孔过程中切削力预报和稳定性判别方法 | |
| Sun et al. | Chatter stability of orthogonal turn-milling analyzed by complete discretization method | |
| Tang et al. | Stability prediction based effect analysis of tool orientation on machining efficiency for five-axis bull-nose end milling | |
| CN103268430A (zh) | 基于机床刀具动刚度测量的铣削工艺参数优化方法 | |
| Reith et al. | Optimal detuning of a parallel turning system—theory and experiments | |
| Cai et al. | An improved cutting force prediction model in the milling process with a multi-blade face milling cutter based on FEM and NURBS | |
| CN106096146A (zh) | 切削过程中薄壁件动力学参数的预测方法 | |
| CN103394972A (zh) | 基于声发射信号的铣削加工表面粗糙度在线预测方法 | |
| CN108647413B (zh) | 一种微细表面位置误差与稳定性综合预测方法 | |
| Nithyanandam et al. | A Six Sigma approach for precision machining in milling | |
| CN109446721B (zh) | 基于标识符软件线程执行顺序排列的机床工艺交互算法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |