CN112380749A - 基于滚压轮优化设计数学模型的制造方法 - Google Patents
基于滚压轮优化设计数学模型的制造方法 Download PDFInfo
- Publication number
- CN112380749A CN112380749A CN202011277550.XA CN202011277550A CN112380749A CN 112380749 A CN112380749 A CN 112380749A CN 202011277550 A CN202011277550 A CN 202011277550A CN 112380749 A CN112380749 A CN 112380749A
- Authority
- CN
- China
- Prior art keywords
- rolling
- rolling wheel
- screw
- wheel
- axis
- 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
Links
- 238000005096 rolling process Methods 0.000 title claims abstract description 123
- 238000013461 design Methods 0.000 title claims abstract description 31
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 17
- 238000005457 optimization Methods 0.000 title claims abstract description 17
- 238000013178 mathematical model Methods 0.000 title claims abstract description 9
- 238000000034 method Methods 0.000 claims abstract description 43
- 238000012360 testing method Methods 0.000 claims abstract description 21
- 238000004088 simulation Methods 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 abstract description 9
- 238000004458 analytical method Methods 0.000 abstract description 5
- 238000013433 optimization analysis Methods 0.000 abstract description 4
- 238000013386 optimize process Methods 0.000 abstract 1
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000013401 experimental design Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
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)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了基于滚压轮优化设计数学模型的制造方法，本发明建立了滚压轮成形段角度、滚压轮转速与滚压结果之间的关系，基于有限元方法获得滚压轮在滚压加工过程中的总做功，利用响应面法进行仿真试验分析，获得优化工艺参数。本发明在滚压成形机理分析的基础上，采用数值仿真和试验优化设计相结合的方法，解决了传统丝杠滚压工艺优化方法不能面向关键工艺变量进行针对性优化分析，且需要依靠大量物理试验的问题，缩短了工艺设计周期，节约制造成本。本发明特别适用于大长径比滚珠丝杠轴向进给加工，也适用于滑动丝杠及行星滚柱丝杠加工制造。
Description
技术领域：
本发明涉及基于滚压轮优化设计数学模型的制造方法。
背景技术：
丝杠副(包括滚珠丝杠副)是机械设备上最常使用的传动部件，其主要功能是将旋转运动转换成线性运动，或将扭矩转换成轴向反复作用力，同时兼具高精度、可逆性和高效率的特点。在工程机械、能源装备、海洋装备等大型设备中使用的滚珠丝杠具有超大长径比，轴向进给滚压成形是其核心制造技术。然而，丝杠滚压成形过程属于大塑性变形过程，局部高应变及应变率导致其成形机理复杂，难以通过有限的物理试验建立滚压工艺参数与滚压结果之间的定量关系。现有技术主要存在的问题：1)工艺设计长期依靠工人经验及试错方法，面对新型丝杠滚压工艺设计缺少实践与理论基础，导致加工成本高，工艺设计效率低，加工质量难以保证；2)多头丝杠滚压成形工艺系统参数繁多，缺少快速高效的方法针对关键工艺参数进行优化。
发明内容：
本发明是为了解决上述现有技术存在的问题而提供了基于滚压轮优化设计数学模型的制造方法。本发明在滚压成形机理分析的基础上，采用数值仿真和试验优化设计相结合的方法，解决了传统丝杠滚压工艺优化方法不能面向关键工艺变量进行针对性优化分析，且需要依靠大量物理试验的问题，缩短了工艺设计周期，节约制造成本。
本发明所采用的技术方案有：
基于滚压轮优化设计数学模型的制造方法，包括：
S1)确定丝杠滚压工艺基本参数，包括丝杠毛坯轴线与滚压轮轴线之间的夹角α，丝杠毛坯轴线与滚压轮轴线之间的距离da；
S2)确定丝杠滚压工艺优化参数，包括滚压轮转速ωr以及滚压轮进料段锥角θ；
S3)计算丝杠轴向进给滚压工艺过程中滚压轮与丝杠毛坯之间相对瞬时速度v：
v＝πdr1ωr/60
其中，dr1为滚压轮中径，单位为：mm；ωr为滚压轮转速，单位为：转/分钟；
S4)计算丝杠轴向进给滚压工艺过程中丝杠毛坯的成形截面每次滚压余量ap：
其中，θ为滚压轮进料段锥角，α为丝杠毛坯轴线与滚压轮轴线之间的夹角，ds1为丝杠中径；
S5)根据丝杠毛坯参数和滚压轮结构设计参数，利用三维建模软件构建丝杠毛坯及滚压轮的参数化几何模型；
S6)将步骤S5中得到的丝杠毛坯及滚压轮的参数化几何模型导入到有限元软件中，然后在有限元软件中选择单元类型、定义材料模型、划分网格、定义接触关系、确定边界及载荷，得丝杠滚压工艺有限元模型；
S7)采用响应面试验设计方法，以滚压瞬时速度v和每次滚压余量ap为变量，以丝杠滚压过程中滚压轮总做功为输出，进行丝杠滚压仿真试验；
S8)根据步骤S7中得到的试验结果，以滚压轮总做功最小为目标，求解滚压参数v和ap，进一步根据步骤S3和S4求解获得优化后的滚压轮转速ωr以及滚压轮进料段锥角θ。
进一步地，步骤S7中的试验设计方法为中心组合设计。
本发明具有如下有益效果：
本发明在滚压成形机理分析的基础上，采用数值仿真和试验优化设计相结合的方法，解决了传统丝杠滚压工艺优化方法不能面向关键工艺变量进行针对性优化分析，且需要依靠大量物理试验的问题，缩短了工艺设计周期，节约制造成本。此方法特别适用于大长径比滚珠丝杠轴向进给加工，也适用于滑动丝杠及行星滚柱丝杠的加工制造。
附图说明：
图1为本发明方法的流程图。
图2和图3为本发明中丝杠滚压工艺基本参数示意图。
图4为本发明中滚压轮示意图。
具体实施方式：
下面结合附图对本发明作进一步的说明。
如图1至图4，本发明基于滚压轮优化设计数学模型的制造方法，本实施设置滚压轮中径为dr1＝200mm，丝杠中径ds1＝17.4mm，具体步骤为：
S1)确定丝杠滚压工艺基本参数，包括丝杠毛坯轴线与滚压轮轴线之间的夹角α＝1.6°，丝杠毛坯轴线与滚压轮轴线之间的距离da＝(dr1+ds1)/2＝108.7mm。
S2)确定丝杠滚压工艺优化参数，包括滚压轮转速ωr以及滚压轮进料段锥角θ；
S3)计算丝杠轴向进给滚压工艺过程中滚压轮与毛坯之间相对瞬时速度v：
ν＝πdr1ωr/60
其中，dr1为滚压轮外径，单位为“mm”；ωr为滚压轮转速，单位为“转/分钟”；
ds1为丝杠中径，单位为“mm”；ωs为丝杠毛坯转速，单位为“转/分钟”，drz和dsz分别为滚压轮和丝杠毛坯之间做纯滚动时滚压轮和丝杠毛坯对应的直径；
S4)计算丝杠轴线进给滚压工艺过程中成形截面每次滚压余量ap：
其中，θ为滚压轮进料段锥角，α为丝杠毛坯轴线与滚压轮轴线之间的夹角；
S5)根据丝杠毛坯参数和滚压轮结构设计参数，利用三维建模软件构建丝杠毛坯及滚压轮的参数化几何模型；
S6)将步骤S5得到的丝杠毛坯及滚压轮的参数化几何模型导入到有限元软件中，然后在有限元软件中选择单元类型、定义材料模型、划分网格、定义接触关系、确定边界即载荷，得丝杠滚压工艺有限元模型；
S7)采用响应面试验设计方法，以滚压瞬时速度v和每次滚压余量ap为变量，以丝杠滚压过程中滚压轮总做功为输出，进行丝杠滚压仿真试验；进行二因子五水平的正交复合(CCD)试验设计，如表1。根据响应曲面法的CCD试验方案，利用最小二乘法拟合响应值与自变量之间的关系方程。对试验数据进行多元回归拟合，可得以滚压轮做功为响应值的多元线性回归模型：
Y＝235.17+77.35A-28.03B+0.63AB-8.41A2+6.81B2+18.65A2B-9.44AB2-44.43A2B2
其中，A为滚压轮的进给量，B为滚压轮的的滚压速度。
S8)根据步骤S7得到的试验结果，以滚压轮总做功最小为目标，考虑滚压效率，以滚压速度为约束，设滚压轮最低转速为5RPM,最高转速为15RPM,则滚压速度v最小值约为52mm/s，最大值约为157mm/s，以此为约束根据响应面模型，以做功最小为优化目标，获得结果：当每次滚压进给量为0.01mm，滚压速度为157时，需要做功最小，为139N.m。即滚压轮转速：ωr＝60v/πdr1＝15RPM，且滚压轮几何参数满足关系式：
表1为本实施例响应面设计及结果
本发明在滚压成形机理分析的基础上，采用数值仿真和试验优化设计相结合的方法，解决了传统丝杠滚压工艺优化方法不能面向关键工艺变量进行针对性优化分析，且需要依靠大量物理试验的问题，缩短了工艺设计周期，节约制造成本。此方法特别适用于大长径比滚珠丝杠轴向进给加工，也适用于滑动丝杠及行星滚柱丝杠的加工制造。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下还可以作出若干改进，这些改进也应视为本发明的保护范围。
Claims (2)
1.基于滚压轮优化设计数学模型的制造方法，其特征在于：包括：
S1)确定丝杠滚压工艺基本参数，包括丝杠毛坯轴线与滚压轮轴线之间的夹角α，丝杠毛坯轴线与滚压轮轴线之间的距离da；
S2)确定丝杠滚压工艺优化参数，包括滚压轮转速ωr以及滚压轮进料段锥角θ；
S3)计算丝杠轴向进给滚压工艺过程中滚压轮与丝杠毛坯之间相对瞬时速度v：
v＝πdr1ωr/60
其中，dr1为滚压轮中径，单位为：mm；ωr为滚压轮转速，单位为：转/分钟；
S4)计算丝杠轴向进给滚压工艺过程中丝杠毛坯的成形截面每次滚压余量ap：
其中，θ为滚压轮进料段锥角，α为丝杠毛坯轴线与滚压轮轴线之间的夹角，ds1为丝杠中径；
S5)根据丝杠毛坯参数和滚压轮结构设计参数，利用三维建模软件构建丝杠毛坯及滚压轮的参数化几何模型；
S6)将步骤S5中得到的丝杠毛坯及滚压轮的参数化几何模型导入到有限元软件中，然后在有限元软件中选择单元类型、定义材料模型、划分网格、定义接触关系、确定边界及载荷，得丝杠滚压工艺有限元模型；
S7)采用响应面试验设计方法，以滚压瞬时速度v和每次滚压余量ap为变量，以丝杠滚压过程中滚压轮总做功为输出，进行丝杠滚压仿真试验；
S8)根据步骤S7中得到的试验结果，以滚压轮总做功最小为目标，求解滚压参数v和ap，进一步根据步骤S3和S4求解获得优化后的滚压轮转速ωr以及滚压轮进料段锥角θ。
2.如权利要求1所述的基于滚压轮优化设计数学模型的制造方法，其特征在于：步骤S7中的试验设计方法为中心组合设计。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011277550.XA CN112380749A (zh) | 2020-11-16 | 2020-11-16 | 基于滚压轮优化设计数学模型的制造方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011277550.XA CN112380749A (zh) | 2020-11-16 | 2020-11-16 | 基于滚压轮优化设计数学模型的制造方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112380749A true CN112380749A (zh) | 2021-02-19 |
Family
ID=74584753
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011277550.XA Pending CN112380749A (zh) | 2020-11-16 | 2020-11-16 | 基于滚压轮优化设计数学模型的制造方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112380749A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113297758A (zh) * | 2021-05-06 | 2021-08-24 | 南昌航空大学 | 一种大型复杂筋板件预锻成形初始毛坯优化设计方法 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106815407A (zh) * | 2016-12-22 | 2017-06-09 | 四川大学 | 一种数控机床整机动态性能优化方法 |
| CN107784172A (zh) * | 2017-10-26 | 2018-03-09 | 南京工程学院 | 波纹换热管制造方法 |
| CN108279121A (zh) * | 2018-02-08 | 2018-07-13 | 湘潭大学 | 一种大切深下滚刀刃底接触力特性测试系统及其测试方法 |
| CN108572072A (zh) * | 2018-04-09 | 2018-09-25 | 上海理工大学 | 基于频响函数法的滚珠丝杠进给系统动态参数辨识方法 |
| CN110362951A (zh) * | 2019-07-24 | 2019-10-22 | 常州工学院 | 滚珠丝杠进给系统结合部虚拟材料模型参数获取方法 |
| CN110569586A (zh) * | 2019-08-28 | 2019-12-13 | 北京自动化控制设备研究所 | 行星滚柱丝杠副精密螺纹切削建模方法 |
| CN110578781A (zh) * | 2019-09-17 | 2019-12-17 | 李敬宇 | 一种小微型行星滚柱丝杠副的丝杠或滚柱的加工方法 |
| CN111222273A (zh) * | 2020-01-03 | 2020-06-02 | 同济大学 | 一种轻质薄板含胶滚合成型耦合建模及其质量预测方法 |
-
2020
- 2020-11-16 CN CN202011277550.XA patent/CN112380749A/zh active Pending
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106815407A (zh) * | 2016-12-22 | 2017-06-09 | 四川大学 | 一种数控机床整机动态性能优化方法 |
| CN107784172A (zh) * | 2017-10-26 | 2018-03-09 | 南京工程学院 | 波纹换热管制造方法 |
| CN108279121A (zh) * | 2018-02-08 | 2018-07-13 | 湘潭大学 | 一种大切深下滚刀刃底接触力特性测试系统及其测试方法 |
| CN108572072A (zh) * | 2018-04-09 | 2018-09-25 | 上海理工大学 | 基于频响函数法的滚珠丝杠进给系统动态参数辨识方法 |
| CN110362951A (zh) * | 2019-07-24 | 2019-10-22 | 常州工学院 | 滚珠丝杠进给系统结合部虚拟材料模型参数获取方法 |
| CN110569586A (zh) * | 2019-08-28 | 2019-12-13 | 北京自动化控制设备研究所 | 行星滚柱丝杠副精密螺纹切削建模方法 |
| CN110578781A (zh) * | 2019-09-17 | 2019-12-17 | 李敬宇 | 一种小微型行星滚柱丝杠副的丝杠或滚柱的加工方法 |
| CN111222273A (zh) * | 2020-01-03 | 2020-06-02 | 同济大学 | 一种轻质薄板含胶滚合成型耦合建模及其质量预测方法 |
Non-Patent Citations (5)
| Title |
|---|
| 付振彪 等: "滚珠丝杠进给系统动力学建模与动态特性分析", 振动与冲击, vol. 38, no. 16, 28 August 2019 (2019-08-28), pages 56 - 63 * |
| 徐光远 等: "基于Pro/E和ANSYS Workbench的滚珠丝杠副造型与有限元分析", 组合机床与自动化加工技术, no. 4, 20 April 2014 (2014-04-20), pages 1 - 5 * |
| 李敬宇 等: "丝杠轴线进给滚压成形工艺参数的优化设计", 世界制造技术与装备市场, no. 6, 15 December 2020 (2020-12-15), pages 65 - 70 * |
| 王生武 等: "滚压强化的残余应力的数值仿真及工艺分析", 计算力学学报, vol. 25, no. 1, 15 December 2008 (2008-12-15), pages 113 - 118 * |
| 钱文海 等: "基于ANSYS的滚珠丝杠副有限元分析", 机械工程师, no. 10, 10 October 2019 (2019-10-10), pages 38 - 39 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113297758A (zh) * | 2021-05-06 | 2021-08-24 | 南昌航空大学 | 一种大型复杂筋板件预锻成形初始毛坯优化设计方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104400214B (zh) | 一种搅拌摩擦焊机床顶锻力及前进抗力测控装置及方法 | |
| CN109189001B (zh) | 拖拉机变速箱用立铣刀端截形图像扫描获取及标定的方法 | |
| CN106825340A (zh) | 一种钛合金c型截面锥度环件轧制的坯料设计方法 | |
| CN112380749A (zh) | 基于滚压轮优化设计数学模型的制造方法 | |
| CN105223814A (zh) | 渐开线齿轮成形砂轮计算方法 | |
| CN110414138B (zh) | 伺服电机的型号确定方法及装置、机器设备 | |
| CN107273612A (zh) | 螺旋锥齿轮凹模齿面设计方法及螺旋锥齿轮加工方法 | |
| CN101733794A (zh) | 竹材去青、去黄疏解机 | |
| CN107248047A (zh) | 基于加工过程状态熵动态计算的加工过程状态评价方法 | |
| Li et al. | Numerical-experimental investigation on the rabbit ear formation mechanism in gear rolling | |
| CN112464526A (zh) | 无芯模旋压成形旋轮加载路径智能优化方法 | |
| CN109976255A (zh) | 一种用于并联主轴头的运动学标定方法 | |
| CN106695023B (zh) | 一种循环球式变比转向器齿轮副齿条齿廓的加工方法 | |
| CN104985244B (zh) | 大型螺杆压缩机转子齿廓面沿螺旋线逐行铣削的方法 | |
| CN102825342B (zh) | 一种人字齿对中误差检验方法及检验量球 | |
| CN108984951B (zh) | 一种优化机床整机动态特性的方法 | |
| Li et al. | Effect of process factors on the rabbit ear based on numerical simulation and experimental study in gear rolling | |
| CN104148397A (zh) | 一种螺旋孔型斜轧轧辊柔性设计方法 | |
| CN110175372A (zh) | 一种基于母面特征参数的包络面表征方法 | |
| CN109407614B (zh) | 一种数控滚齿机滚齿加工工艺参数优化方法 | |
| CN103934727B (zh) | 一种等角螺旋线分区变参数控精密加工方法 | |
| CN103838182B (zh) | 数控机床编程辅助方法 | |
| Khodaee et al. | Finite element simulation as a tool to evaluate gear quality after gear rolling | |
| CN108628250B (zh) | 一种机床主轴与刀柄配合结构尺寸的优化方法 | |
| CN102102670A (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 |