CN105426565A - 基于概率密度技术的大型高速回转装备形位公差分配方法 - Google Patents
基于概率密度技术的大型高速回转装备形位公差分配方法 Download PDFInfo
- Publication number
- CN105426565A CN105426565A CN201510664083.9A CN201510664083A CN105426565A CN 105426565 A CN105426565 A CN 105426565A CN 201510664083 A CN201510664083 A CN 201510664083A CN 105426565 A CN105426565 A CN 105426565A
- Authority
- CN
- China
- Prior art keywords
- equipment
- large high
- level
- tolerance
- speed revolution
- 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
Landscapes
- Manufacture Of Motors, Generators (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
基于概率密度技术的大型高速回转装备形位公差分配方法属于机械公差分配技术;分析大型高速回转装备的径向和轴向测量面的定位及定向公差在装配中的传递过程,确定<i>n</i>级装备装配后的圆心坐标的传递关系,得到装配后装备偏心与各级装备定位、定向公差和旋转角度之间的关系;依据同轴度公差的目标函数,得到<i>n</i>级装备同轴度公差的概率密度,最终得到各级大型高速回转装备的径向偏心及轴向垂直度公差与最终多级装备同轴度公差的概率关系,实现大型高速回转装备公差的分配。
Description
技术领域
本发明属于机械公差分配技术,主要涉及一种基于概率密度技术的大型高速回转装备形位公差分配方法。
背景技术
大型高速回转装备为机械领域中重要的特征部件,在各种领域中均有广泛应用,如随着航空事业的不断发展,各种军用、民用飞机研发技术水平不断提升,对涡轮发动机性能的要求也不断提高。发动机在追求更高推重比的同时,也需尽量减小自身振动对发动机性能和飞行安全的影响。
发动机振动是影响飞机安全的一个重要因素,也是反映发动机性能的一项重要指标。发动机涡轮部件转速高、质量大,是发动机的一个主要振源。为了降低此影响,除了保证发动机在装配过程中的质量,还要通过在设计航空发动机转静子公差分配上进行合理赋值,由装配不合理,主要是装配后多级转子的同轴度误差较大,导致振动在高速运转时会放大100至1000倍,装配时消除偏摆可以很大程度上减小动平衡的压力。所以,作为提升航空发动机性能的关键技术,航空发动机转静子公差分配技术越来越受到重视,并且成为研究的热点。
波音公司提出一种统计公差确定方法(统计公差确定方法。公开号:CN1549069A),选择一优选的零件装配顺序并对其有效性进行确认。选择待加工的配合件的定位件、数量和尺寸,借此使各部件彼此相对定位或紧固在一起而形成装配件。将各个单独零部件的统计尺寸公差确定为零部件的制造要求。利用此公差可经济地进行加工,并使所组装的装配件满足装配件的尺寸公差。
北京航空航天大学一种基于公差等级和遗传算法的公差优化分配方法(一种基于公差等级和遗传算法的公差优化分配方法。公开号:CN102622495A),该方法包括以下步骤:(一)建立公差优化的数学模型;(二)建立标准公差数据库;(三)基于公差等级进行编码;(四)确定遗传算法的适应度函数;(五)设定遗传算子;(六)遗传进化;(七)解码得到对应的公差等级,并查询相应的公差值。
北京理工大学一种基于正交试验原理的尺寸链公差分配方法(基于正交试验原理的尺寸链公差分配方法。公开号:CN102622480A),该方法的步骤为:根据装配体或零件设计图获取尺寸链组成环的基本信息、封闭环的基本值以及封闭环的表达式;根据组成环个数,选取相应正交表;根据封闭环表达式,求解每个组成环的灵敏度;根据组成环的灵敏度,构造对应于封闭环最大极限值情况下的组成环偏差集合LU以及对应于封闭环最小极限值情况下的组成环偏差集合LD;分别优化分配组成环偏差集合LU和LD;输出结果。
以上方法均存在的问题在于:提出一种普遍的公差分配方法,并未考虑大型高速回转装备装配过程中的公差传递过程,依据其传递过程的特点进而设计一种应用于大型高速回转装备的公差分配方法。
发明内容
针对上述现有技术存在的不足,提出一种基于概率密度技术的大型高速回转装备形位公差分配方法,以解决大型高速回转装备装配同轴度低、装配质量差,振动大的问题,改善大型高速回转装备的性能。
本发明的目的是这样实现的:
一种基于概率密度技术的大型高速回转装备形位公差分配方法,
基于概率密度技术的大型高速回转装备形位公差分配方法,其特征在于:
定义4×4矩阵,其中R为3×3旋转矩阵,表示绕X、Y和Z轴的角度,p为1×3矢量,表示在X、Y和Z轴方向位置公差;
大型高速回转装备的径向测量面的定位公差在装配中的传递矩阵为T ecc ,则偏心矩阵,其中e为大型高速回转装备的径向测量面偏心,θ e 为大型高速回转装备的径向测量面偏心角,H为大型高速回转装备高度;
大型高速回转装备的轴向测量面的定向公差在装配中的传递矩阵为T tilt ,则倾斜矩阵,其中θ t 为大型高速回转装备的轴向测量面倾斜角,θ l 为大型高速回转装备的轴向测量面的拟合平面最低点所在角度;
则n个装备的装配后圆心坐标的传递关系T n 为,(n=2,3,4,…),其中T ri ,为i级装备的旋转矩阵,T ecci 为i级装备的偏心矩阵和T perpi 为i级装备的倾斜矩阵(i=1,2,…,n),其中1级装备的旋转矩阵T r1=E,n级装备的倾斜矩阵T tiltn =E,旋转矩阵,其中θ ri 为i级大型高速回转装备的旋转角度;则装配后各级装备在X,Y轴方向的分量为(n=2,3,4,…),其中e i 为i级大型高速回转装备的径向测量面偏心,θ ei 为i级大型高速回转装备的径向测量面偏心角,H i 为i级大型高速回转装备的高度,θ ri 为i级大型高速回转装备的旋转角度,θ li 为大型高速回转装备的轴向测量面的拟合平面最低点所在角度,r i 为i级大型高速回转装备的轴向测量面半径,h i 为i级大型高速回转装备的轴向测量面垂直度,(i=1,2,…,n);
根据公式可以得到n级装备的同轴度公差,当级数n(n=2,3,4,…)确定时,可以得到同轴度公差的目标函数,进而得到n级装备同轴度公差的概率密度,最终得到各级大型高速回转装备的径向偏心及轴向垂直度公差与最终多级装备同轴度公差的概率关系,实现大型高速回转装备公差的分配。
与现有技术相比,本发明的特点是:
本发明分析大型高速回转装备的径向和轴向测量面的定位及定向公差在装配中的传递过程,确定n级装备装配后的圆心坐标的传递关系,得到装配后装备偏心与各级装备定位、定向公差和旋转角度之间的关系;依据各级大型高速回转装备的径向偏心及轴向垂直度之间的关系,绘制出分布函数求出概率密度函数,进而得到各级大型高速回转装备的径向偏心及轴向垂直度公差与最终多级装备同轴度公差的概率关系,实现大型高速回转装备公差的分配。
附图说明
图1是大型高速回转装备装配示意图。
具体实施方式
下面结合附图对本发明实施方案作进一步详细描述:
一种基于概率密度技术的大型高速回转装备形位公差分配方法,所述方法是:定义4×4矩阵,其中R为3×3旋转矩阵,表示绕X、Y和Z轴的角度,p为1×3矢量,表示在X、Y和Z轴方向位置公差;
O为部件1的圆心,O 1A为部件1装配面的圆心,O 1B为部件2的圆心,O 2为部件2装配面的圆心;大型高速回转装备的径向测量面的定位公差在装配中的传递矩阵为T ecc ,可以描述部件1中O到O 1A圆心和部件2中O 1B到O 2圆心坐标间的关系,则偏心矩阵,其中e为大型高速回转装备的径向测量面偏心,θ e 为大型高速回转装备的径向测量面偏心角,H为大型高速回转装备高度;
大型高速回转装备的轴向测量面的定向公差在装配中的传递矩阵为T tilt ,可以描述部件1和部件2装配过程中,由于装配面的垂直度参数所引起的O 1A圆心到O 1B圆心坐标间关系,则倾斜矩阵,其中θ t 为大型高速回转装备的轴向测量面倾斜角,θ l 为大型高速回转装备的轴向测量面的拟合平面最低点所在角度;
则n个装备的装配后圆心坐标的传递关系T n 为,(n=2,3,4,…),其中T ri ,为i级装备的旋转矩阵,T ecci 为i级装备的偏心矩阵和T perpi 为i级装备的倾斜矩阵(i=1,2,…,n),其中1级装备的旋转矩阵T r1=E,n级装备的倾斜矩阵T tiltn =E,旋转矩阵,其中θ ri 为i级大型高速回转装备的旋转角度;则装配后各级装备在X,Y轴方向的分量为(n=2,3,4,…),其中e i 为i级大型高速回转装备的径向测量面偏心,θ ei 为i级大型高速回转装备的径向测量面偏心角,H i 为i级大型高速回转装备的高度,θ ri 为i级大型高速回转装备的旋转角度,θ li 为大型高速回转装备的轴向测量面的拟合平面最低点所在角度,r i 为i级大型高速回转装备的轴向测量面半径,h i 为i级大型高速回转装备的轴向测量面垂直度,(i=1,2,…,n);
根据公式可以得到n级装备的同轴度公差,当级数n(n=2,3,4,…)确定时,可以得到同轴度公差的目标函数,进而得到n级装备同轴度公差的概率密度,最终得到各级大型高速回转装备的径向偏心及轴向垂直度公差与最终多级装备同轴度公差的概率关系,实现大型高速回转装备公差的分配。
Claims (1)
1.基于概率密度技术的大型高速回转装备形位公差分配方法,其特征在于:
定义4×4矩阵,其中R为3×3旋转矩阵,表示绕X、Y和Z轴的角度,p为1×3矢量,表示在X、Y和Z轴方向位置公差;
大型高速回转装备的径向测量面的定位公差在装配中的传递矩阵为T ecc ,则偏心矩阵,其中e为大型高速回转装备的径向测量面偏心,θ e 为大型高速回转装备的径向测量面偏心角,H为大型高速回转装备高度;
大型高速回转装备的轴向测量面的定向公差在装配中的传递矩阵为T tilt ,则倾斜矩阵,其中θ t 为大型高速回转装备的轴向测量面倾斜角,θ l 为大型高速回转装备的轴向测量面的拟合平面最低点所在角度;
则n个装备的装配后圆心坐标的传递关系T n 为,(n=2,3,4,…),其中T ri ,为i级装备的旋转矩阵,T ecci 为i级装备的偏心矩阵和T perpi 为i级装备的倾斜矩阵(i=1,2,…,n),其中1级装备的旋转矩阵T r1=E,n级装备的倾斜矩阵T tiltn =E,旋转矩阵,其中θ ri 为i级大型高速回转装备的旋转角度;则装配后各级装备在X,Y轴方向的分量为(n=2,3,4,…),其中e i 为i级大型高速回转装备的径向测量面偏心,θ ei 为i级大型高速回转装备的径向测量面偏心角,H i 为i级大型高速回转装备的高度,θ ri 为i级大型高速回转装备的旋转角度,θ li 为大型高速回转装备的轴向测量面的拟合平面最低点所在角度,r i 为i级大型高速回转装备的轴向测量面半径,h i 为i级大型高速回转装备的轴向测量面垂直度,(i=1,2,…,n);
根据公式可以得到n级装备的同轴度公差,当级数n(n=2,3,4,…)确定时,可以得到同轴度公差的目标函数,进而得到n级装备同轴度公差的概率密度,最终得到各级大型高速回转装备的径向偏心及轴向垂直度公差与最终多级装备同轴度公差的概率关系,实现大型高速回转装备公差的分配。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510664083.9A CN105426565B (zh) | 2015-10-15 | 2015-10-15 | 基于概率密度技术的大型高速回转装备形位公差分配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510664083.9A CN105426565B (zh) | 2015-10-15 | 2015-10-15 | 基于概率密度技术的大型高速回转装备形位公差分配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105426565A true CN105426565A (zh) | 2016-03-23 |
CN105426565B CN105426565B (zh) | 2018-02-16 |
Family
ID=55504775
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510664083.9A Expired - Fee Related CN105426565B (zh) | 2015-10-15 | 2015-10-15 | 基于概率密度技术的大型高速回转装备形位公差分配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105426565B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109960870A (zh) * | 2019-03-21 | 2019-07-02 | 哈尔滨工业大学 | 基于接触面积极大化测调的大型高速回转装备多级零部件刚度预测方法 |
CN110008433A (zh) * | 2019-04-10 | 2019-07-12 | 重庆大学 | 一种基于状态空间方程的元动作单元零件公差分配方法 |
CN110598229A (zh) * | 2019-01-07 | 2019-12-20 | 哈尔滨工业大学 | 基于蒙特卡罗偏差评估的大型高速回转装备多级零部件优化方法 |
CN110595690A (zh) * | 2019-01-07 | 2019-12-20 | 哈尔滨工业大学 | 基于形心质心重心惯性中心矢量极小化的大型高速回转装备测量与智能学习装配方法和装置 |
CN111076867A (zh) * | 2018-10-22 | 2020-04-28 | 哈尔滨工业大学 | 基于质心和惯性中心同步测调的大型高速回转装备多级零部件不平衡量分配方法 |
CN111076866A (zh) * | 2018-10-22 | 2020-04-28 | 哈尔滨工业大学 | 基于形心矢量和极小化的大型高速回转装备多级零部件不平衡量堆叠装配方法和装置 |
CN111475890A (zh) * | 2019-01-07 | 2020-07-31 | 哈尔滨工业大学 | 基于偏心矢量跟随测调的大型高速回转装备零部件相对间隙参数分配方法 |
CN111475903A (zh) * | 2019-01-07 | 2020-07-31 | 哈尔滨工业大学 | 基于多偏置误差同步补偿的大型高速回转装备多级零部件动力学特性逐级测调和分配方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2222324C (en) * | 1995-06-28 | 2005-10-04 | The Boeing Company | Statistical tolerancing |
CN102426615A (zh) * | 2011-09-01 | 2012-04-25 | 北京理工大学 | 一种用于精密机械系统误差传递建模的配合误差计算方法 |
CN102622480A (zh) * | 2012-03-02 | 2012-08-01 | 北京理工大学 | 基于正交试验原理的尺寸链公差分配方法 |
CN102622495A (zh) * | 2012-04-13 | 2012-08-01 | 北京航空航天大学 | 一种基于公差等级和遗传算法的公差优化分配方法 |
CN103902759A (zh) * | 2013-12-27 | 2014-07-02 | 西京学院 | 一种基于遗传算法的装配公差优化设计方法 |
CN104850698A (zh) * | 2015-05-15 | 2015-08-19 | 西安交通大学 | 一种考虑测量与调整工艺的精密机床公差设计方法 |
-
2015
- 2015-10-15 CN CN201510664083.9A patent/CN105426565B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2222324C (en) * | 1995-06-28 | 2005-10-04 | The Boeing Company | Statistical tolerancing |
CN102426615A (zh) * | 2011-09-01 | 2012-04-25 | 北京理工大学 | 一种用于精密机械系统误差传递建模的配合误差计算方法 |
CN102622480A (zh) * | 2012-03-02 | 2012-08-01 | 北京理工大学 | 基于正交试验原理的尺寸链公差分配方法 |
CN102622495A (zh) * | 2012-04-13 | 2012-08-01 | 北京航空航天大学 | 一种基于公差等级和遗传算法的公差优化分配方法 |
CN103902759A (zh) * | 2013-12-27 | 2014-07-02 | 西京学院 | 一种基于遗传算法的装配公差优化设计方法 |
CN104850698A (zh) * | 2015-05-15 | 2015-08-19 | 西安交通大学 | 一种考虑测量与调整工艺的精密机床公差设计方法 |
Non-Patent Citations (3)
Title |
---|
DESROCHERS A, RIVIÈRE A.: "A Matrix Approach to the Representation of Tolerance Zones and Clearances", 《THE INTERNATIONAL JOURNAL OF ADVANCED MANUFACTURING TECHNOLOGY》 * |
WANG L, TAN J B, LIU Y T.: "Research on giant magnetostrictive micro-displacement actuator with self-adaptive control algorithm", 《JOURNAL OF PHYSICS: CONFERENCE SERIES.》 * |
ZHAOBO W L.: "Circular trajectory motion control of an inspection spherical robot", 《PROC. OF SPIE VOL. 2013》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111076867A (zh) * | 2018-10-22 | 2020-04-28 | 哈尔滨工业大学 | 基于质心和惯性中心同步测调的大型高速回转装备多级零部件不平衡量分配方法 |
CN111076866A (zh) * | 2018-10-22 | 2020-04-28 | 哈尔滨工业大学 | 基于形心矢量和极小化的大型高速回转装备多级零部件不平衡量堆叠装配方法和装置 |
CN111076866B (zh) * | 2018-10-22 | 2021-03-30 | 哈尔滨工业大学 | 基于形心矢量和极小化的大型高速回转装备多级零部件不平衡量堆叠装配方法和装置 |
CN110598229A (zh) * | 2019-01-07 | 2019-12-20 | 哈尔滨工业大学 | 基于蒙特卡罗偏差评估的大型高速回转装备多级零部件优化方法 |
CN110595690A (zh) * | 2019-01-07 | 2019-12-20 | 哈尔滨工业大学 | 基于形心质心重心惯性中心矢量极小化的大型高速回转装备测量与智能学习装配方法和装置 |
CN111475890A (zh) * | 2019-01-07 | 2020-07-31 | 哈尔滨工业大学 | 基于偏心矢量跟随测调的大型高速回转装备零部件相对间隙参数分配方法 |
CN111475903A (zh) * | 2019-01-07 | 2020-07-31 | 哈尔滨工业大学 | 基于多偏置误差同步补偿的大型高速回转装备多级零部件动力学特性逐级测调和分配方法 |
CN110598229B (zh) * | 2019-01-07 | 2020-10-16 | 哈尔滨工业大学 | 基于蒙特卡罗偏差评估的大型高速回转装备多级零部件优化方法 |
CN109960870A (zh) * | 2019-03-21 | 2019-07-02 | 哈尔滨工业大学 | 基于接触面积极大化测调的大型高速回转装备多级零部件刚度预测方法 |
CN109960870B (zh) * | 2019-03-21 | 2022-06-14 | 哈尔滨工业大学 | 基于接触面积极大化测调的大型高速回转装备多级零部件刚度预测方法 |
CN110008433A (zh) * | 2019-04-10 | 2019-07-12 | 重庆大学 | 一种基于状态空间方程的元动作单元零件公差分配方法 |
CN110008433B (zh) * | 2019-04-10 | 2023-02-07 | 重庆大学 | 一种基于状态空间方程的元动作单元零件公差分配方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105426565B (zh) | 2018-02-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105426565A (zh) | 基于概率密度技术的大型高速回转装备形位公差分配方法 | |
CN105426566A (zh) | 基于蒙特卡罗技术的大型高速回转装备形位公差分配方法 | |
US11480490B2 (en) | Large-scale high-speed rotary equipment measuring and intelligent learning assembly method and device based on vector minimization of geometry center, mass center, center of gravity and center of inertia | |
CN110598229B (zh) | 基于蒙特卡罗偏差评估的大型高速回转装备多级零部件优化方法 | |
CN109960868B (zh) | 基于多参数调控的大型高速回转装备单级零部件偏心数据处理及跳动公差分配方法 | |
CN111460677B (zh) | 一种基于几何代数理论建立转子堆叠精度预测模型的方法 | |
CN111413030B (zh) | 基于刚度矢量空间投影极大化的大型高速回转装备测量与神经网络学习调控方法及其装置 | |
CN110119553B (zh) | 一种止口连接的航空发动机转子零件选配优化方法 | |
CN109960867A (zh) | 一种基于惯性回转中心调控最优化及智能学习的大型高速回转装备多级零部件选配方法 | |
Wang et al. | Improvement of location and orientation tolerances propagation control in cylindrical components assembly using stack-build assembly technique | |
CN103308021B (zh) | 一种测量工件圆度误差的方法 | |
CN108829970A (zh) | 基于涡动力学的轴流风扇叶片优化设计方法及优化系统 | |
CN108331899A (zh) | 一种rv减速器摆线齿轮齿廓修形设计方法及系统 | |
CN110193708A (zh) | 超大直径多段式齿圈加工方法 | |
CN110877751A (zh) | 基于矢量投影的大型高速回转装备转动惯量堆叠方法 | |
CN111982402B (zh) | 考虑初始不平衡量的航空发动机转子零件选配优化方法 | |
CN111076867B (zh) | 基于质心和惯性中心同步测调的大型高速回转装备多级零部件不平衡量分配方法 | |
CN110717293A (zh) | 一种转子止口螺栓连接结合面变形规律拟合方法 | |
CN104680016A (zh) | 基于几何优化逼近的抛物线轮廓最小区域拟合方法 | |
CN104915500B (zh) | 粉末注射成形涡轮及其结构优化设计方法 | |
CN104359622A (zh) | 一种整体叶轮去除不平衡量方法 | |
US20200217223A1 (en) | Method for Distributing Relative Gap Parameters of Large-Scale High-Speed Rotary Equipment Components Based on Eccentricity Vector Following Measurement and Adjustment | |
CN104180791A (zh) | 一种飞机水平测量方法 | |
CN111666639B (zh) | 一种消除曲轴型转子装配的优化方法 | |
CN104583604A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180216 Termination date: 20201015 |