CN110877751B - 基于矢量投影的大型高速回转装备转动惯量堆叠方法 - Google Patents
基于矢量投影的大型高速回转装备转动惯量堆叠方法 Download PDFInfo
- Publication number
- CN110877751B CN110877751B CN201911227462.6A CN201911227462A CN110877751B CN 110877751 B CN110877751 B CN 110877751B CN 201911227462 A CN201911227462 A CN 201911227462A CN 110877751 B CN110877751 B CN 110877751B
- Authority
- CN
- China
- Prior art keywords
- rotor
- matrix
- stage
- inertia
- global
- 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.)
- Active
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
- B64F5/00—Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
- B64F5/10—Manufacturing or assembling aircraft, e.g. jigs therefor
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
- B64F5/00—Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Manufacturing & Machinery (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Aviation & Aerospace Engineering (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Transportation (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Computing Systems (AREA)
- Turbine Rotor Nozzle Sealing (AREA)
Abstract
本发明提供了一种基于矢量投影的大型高速回转装备转动惯量堆叠方法,包括首先定义装配体全局坐标系、确定n级转子装配过程中,第i级转子的加工误差矩阵、安装相位矩阵;确定n级转子装配后各级转子由于坐标系平移变换引起的惯性张量传递矩阵、各级转子由于坐标系旋转变换引起的惯性张量传递矩阵,其次确定经过平移与旋转综合变换后第k级转子的全局惯性张量矩阵;然后确定装配体的全局惯性张量矩阵、确定n级转子装配绕转轴的转动惯量堆叠优化模型,最后,根据优化模型利用遗传算法寻优计算各级转子安装相位。本发明可以指导航空发动机转子多级盘装配,实现整体转动惯量最优,使航空发动机转子具有良好的启停特性和精确的姿态控制。
Description
技术领域
本发明属于机械装配技术领域,尤其是涉及一种基于矢量投影的大型高速回转装备转动惯量堆叠方法。
背景技术
以航空发动机、舰船燃机、电站燃机、天然气/石油输送增压燃机为代表的大型精密回转装备,由于技术复杂,研制难度大,再加上国外技术封锁,已成为我国高端装备领域中的“卡脖子”问题。
以航空发动机为例,其核心机作为航空发动机的核心部件,其装配质量将直接影响航空发动机的性能,为保证航空发动机核心机能够达到设计目标,必须对每一环节进行严格把控。其中就包含转子质量特性参数的调整装配。航空发动机转子质量特性参数主要是指发动机转子的重心和惯性张量。
其中惯性张量是衡量物体做定点转动时惯性大小的量,其作为航空航天器运动建模的重要参数,是进行动力学分析、姿态控制特性分析和气动弹性分析不可或缺的数据。作为高速旋转的部件,惯性张量同样是航空发动机性能分析的重要参数,其中航空发动机多级转子相对回转轴线的转动惯量尤为重要,其会影响发动机的启停特性。同时对于急剧加速和减速旋转过程也有较大影响,这关系到战机的急速攀升、俯冲等飞行性能。通过合理的方法来有效控制航空发动机多级转子相对回转轴线转动惯量,是改善航空发动机启停、加减速性能的一个重要途径。
航空发动机作为精密机械设备,在生产制造过程中除了单级盘提高加工测量精度外,提高其装配精度也是重中之重。在实际生产中,往往单件加工精度符合要求,但装配后的整机验证不合格。因此,对于自研自制高性能的航空发动机,提高装配精度是当务之急。提高装配精度对于多级转子装配就是实现堆叠优化,而优化的目的是直接改善转子质量特性,减小转动惯量,提高发动机的振动和启停特性,同时也保证一次装配合格率,提高装配效率
国内航空发动机装配没有将转动惯量的调整装配优化纳入到考虑范围之内。无法保证装配后转动惯量最优,非常影响安装的效率和一次成功率。
发明内容
有鉴于此,本发明旨在提出一种基于矢量投影的大型高速回转装备转动惯量堆叠方法,可以指导航空发动机转子多级盘装配,实现整体转动惯量最优,使航空发动机转子具有良好的启停特性和精确的姿态控制。
为达到上述目的,本发明的技术方案是这样实现的:
一种基于矢量投影的大型高速回转装备转动惯量堆叠方法,具体包括以下步骤:
步骤一、定义装配体全局坐标系,在装配过程中,选择第一级转子的下端面圆心作为坐标原点,通过坐标原点且垂直于第一级转子的下端面的轴为Z轴;通过坐标原点,垂直于Z轴且通过第一级转子初始安装位置的轴为X轴,建立装配体全局坐标;
步骤二、确定n级转子装配过程中,第i级转子的加工误差矩阵,所述加工误差矩阵为:
步骤三、确定n级转子装配过程中,第i级转子的安装相位矩阵,其中i=1,2,…,n,所述安装相位矩阵为:
步骤四、确定n级转子装配后各级转子由于坐标系平移变换引起的惯性张量传递矩阵,所述惯性张量传递矩阵为:
其中mk为第k级转子的质量,为第k级转子重心在全局坐标系下的位置坐标的位置向量,为n级转子装配后第k级转子重心在全局坐标系下X轴坐标;为n级转子装配后第k级转子重心在全局坐标系下Y轴坐标;为n级转子装配后第k级转子重心在全局坐标系下Z轴坐标;
步骤五、确定n级转子装配后各级转子由于坐标系旋转变换引起的惯性张量传递矩阵,所述惯性张量传递矩阵为:
其中Ik为第k级转子的重心坐标系下的3×3惯性张量矩阵;
步骤六、将步骤四中由于坐标系平移变换引起的惯性张量传递矩阵与步骤五中由于坐标系旋转变换引起的惯性张量传递矩阵相叠加,确定n级转子装配后,经过平移与旋转综合变换后第k级转子的全局惯性张量矩阵所述全局惯性张量矩阵为:
步骤七、将各级转子的全局惯性张量矩阵逐级累加,确定n级转子装配后,装配体的全局惯性张量矩阵,所述装配体的全局惯性张量矩阵为:
步骤八、根据确定n级转子装配后,装配体的全局惯性张量矩阵,确定n级转子装配绕转轴的转动惯量堆叠优化模型:
相对于现有技术,本发明所述的基于矢量投影的大型高速回转装备转动惯量堆叠方法具有以下优势:
本发明所述的基于矢量投影的大型高速回转装备转动惯量堆叠方法,利用该优化模型指导转子装配,有效降低装配后转子绕其回转轴线的转动惯量,以此指导航空发动机转子多级盘装配,实现整体转动惯量最优,使航空发动机转子具有良好的启停特性和精确的姿态控制,并且本方法还实现了通过各级转子惯量参数,预测装配后的装配体惯量参数,对于转子设计和加工制造以及航空发动机配套零部件的选型提供参考转动惯量数据。
具体实施方式
需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。
下面以实施例来详细说明本发明。
基于矢量投影的大型高速回转装备转动惯量堆叠方法,具体包括以下步骤:
步骤一、定义装配体全局坐标系,在装配过程中,选择第一级转子的下端面圆心作为坐标原点,通过坐标原点且垂直于第一级转子的下端面的轴为Z轴;通过坐标原点,垂直于Z轴且通过第一级转子初始安装位置的轴为X轴,建立装配体全局坐标;
步骤二、确定n级转子装配过程中,第i级转子的加工误差矩阵,其中i=1,2,…,n,所述加工误差矩阵为:
步骤三、确定n级转子装配过程中,第i级转子的安装相位矩阵,i=1,2,…,n,所述安装相位矩阵为:
步骤四、确定n级转子装配后各级转子由于坐标系平移变换引起的惯性张量传递矩阵,所述惯性张量传递矩阵为:
其中mk为第k级转子的质量,为已知条件,为第k级转子重心在全局坐标系下的位置坐标的位置向量,为n级转子装配后第k级转子重心在全局坐标系下X轴坐标;为n级转子装配后第k级转子重心在全局坐标系下Y轴坐标;为n级转子装配后第k级转子重心在全局坐标系下Z轴坐标;
步骤五、确定n级转子装配后各级转子由于坐标系旋转变换引起的惯性张量传递矩阵,所述惯性张量传递矩阵为:
步骤六、将步骤四中由于坐标系平移变换引起的惯性张量传递矩阵与步骤五中由于坐标系旋转变换引起的惯性张量传递矩阵相叠加,确定n级转子装配后,经过平移与旋转综合变换后第k级转子的全局惯性张量矩阵所述全局惯性张量矩阵为:
步骤七、将各级转子的全局惯性张量矩阵逐级累加,确定n级转子装配后,装配体的全局惯性张量矩阵,所述装配体的全局惯性张量矩阵为:
步骤八、根据确定n级转子装配后,装配体的全局惯性张量矩阵,确定n级转子装配绕转轴的转动惯量堆叠优化模型:
本发明结合各级转子自身的惯性张量测量数据,以n级转子装配绕转轴的转动惯量堆叠优化公式为优化目标进行寻优计算,即可得到在特定的各级转子惯性张量数据与转子级数下,每级转子的最佳安装相位,以此指导航空发动机多级盘的装配,实现装配后整体转动惯量最优。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种基于矢量投影的大型高速回转装备转动惯量堆叠方法,其特征在于:具体包括以下步骤:
步骤一、定义装配体全局坐标系,在装配过程中,选择第一级转子的下端面圆心作为坐标原点,通过坐标原点且垂直于第一级转子的下端面的轴为Z轴;通过坐标原点,垂直于Z轴且通过第一级转子初始安装位置的轴为X轴,建立装配体全局坐标;
步骤二、确定n级转子装配过程中,第i级转子的加工误差矩阵,所述加工误差矩阵为:
步骤三、确定n级转子装配过程中,第i级转子的安装相位矩阵,其中i=1,2,…,n,所述安装相位矩阵为:
步骤四、确定n级转子装配后各级转子由于坐标系平移变换引起的惯性张量传递矩阵,所述惯性张量传递矩阵为:
其中mk为第k级转子的质量,为第k级转子重心在全局坐标系下的位置坐标的位置向量,为n级转子装配后第k级转子重心在全局坐标系下X轴坐标;为n级转子装配后第k级转子重心在全局坐标系下Y轴坐标;为n级转子装配后第k级转子重心在全局坐标系下Z轴坐标;
步骤五、确定n级转子装配后各级转子由于坐标系旋转变换引起的惯性张量传递矩阵,所述惯性张量传递矩阵为:
其中Ik为第k级转子的重心坐标系下的3×3惯性张量矩阵;
步骤六、将步骤四中由于坐标系平移变换引起的惯性张量传递矩阵与步骤五中由于坐标系旋转变换引起的惯性张量传递矩阵相叠加,确定n级转子装配后,经过平移与旋转综合变换后第k级转子的全局惯性张量矩阵所述全局惯性张量矩阵为:
步骤七、将各级转子的全局惯性张量矩阵逐级累加,确定n级转子装配后,装配体的全局惯性张量矩阵,所述装配体的全局惯性张量矩阵为:
步骤八、根据确定n级转子装配后,装配体的全局惯性张量矩阵,确定n级转子装配绕转轴的转动惯量堆叠优化模型:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911227462.6A CN110877751B (zh) | 2019-12-04 | 2019-12-04 | 基于矢量投影的大型高速回转装备转动惯量堆叠方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911227462.6A CN110877751B (zh) | 2019-12-04 | 2019-12-04 | 基于矢量投影的大型高速回转装备转动惯量堆叠方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110877751A CN110877751A (zh) | 2020-03-13 |
CN110877751B true CN110877751B (zh) | 2021-04-30 |
Family
ID=69730065
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911227462.6A Active CN110877751B (zh) | 2019-12-04 | 2019-12-04 | 基于矢量投影的大型高速回转装备转动惯量堆叠方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110877751B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111666639B (zh) * | 2020-05-09 | 2022-09-20 | 大连理工大学 | 一种消除曲轴型转子装配的优化方法 |
CN114013685B (zh) * | 2021-11-15 | 2023-09-08 | 中国航发沈阳发动机研究所 | 一种航空发动机核心机试验加温加压程序设计方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103912315A (zh) * | 2014-04-14 | 2014-07-09 | 西北工业大学 | 一种航空发动机转子结构动力学的设计方法 |
CN105426566A (zh) * | 2015-10-15 | 2016-03-23 | 哈尔滨工业大学 | 基于蒙特卡罗技术的大型高速回转装备形位公差分配方法 |
CN107255549A (zh) * | 2017-06-14 | 2017-10-17 | 西安交通大学 | 一种转子多点不平衡量识别方法 |
CN107600173A (zh) * | 2017-09-20 | 2018-01-19 | 南京航空航天大学 | 一种汽车液压变传动比转向系统及其多目标优化方法 |
CN109871947A (zh) * | 2019-03-21 | 2019-06-11 | 哈尔滨工业大学 | 基于卷积神经网络的大型高速回转装备多级零部件初始不平衡量逐级堆叠方法和装置 |
CN109960869A (zh) * | 2019-03-21 | 2019-07-02 | 哈尔滨工业大学 | 基于功能性滤波处理的大型高速回转装备零部件形心和惯性中心数据处理方法 |
CN109960867A (zh) * | 2019-03-21 | 2019-07-02 | 哈尔滨工业大学 | 一种基于惯性回转中心调控最优化及智能学习的大型高速回转装备多级零部件选配方法 |
CN109974749A (zh) * | 2019-04-09 | 2019-07-05 | 哈尔滨工业大学 | 一种三轴转台综合指向误差的评定方法 |
CN110153664A (zh) * | 2019-03-21 | 2019-08-23 | 哈尔滨工业大学 | 基于bp神经网络的大型高速回转装备多级零部件装配方法和装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017020097A2 (ru) * | 2015-02-23 | 2017-02-09 | Национальная Академия Авиации | Способ и устройство повышения стабилизации и маневренности беспилотных летательных аппаратов (бла) с применением гироскопического эффекта |
US10369702B2 (en) * | 2016-10-17 | 2019-08-06 | Raytheon Company | Automated work piece moment of inertia (MOI) identification system and method for same |
-
2019
- 2019-12-04 CN CN201911227462.6A patent/CN110877751B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103912315A (zh) * | 2014-04-14 | 2014-07-09 | 西北工业大学 | 一种航空发动机转子结构动力学的设计方法 |
CN105426566A (zh) * | 2015-10-15 | 2016-03-23 | 哈尔滨工业大学 | 基于蒙特卡罗技术的大型高速回转装备形位公差分配方法 |
CN107255549A (zh) * | 2017-06-14 | 2017-10-17 | 西安交通大学 | 一种转子多点不平衡量识别方法 |
CN107600173A (zh) * | 2017-09-20 | 2018-01-19 | 南京航空航天大学 | 一种汽车液压变传动比转向系统及其多目标优化方法 |
CN109871947A (zh) * | 2019-03-21 | 2019-06-11 | 哈尔滨工业大学 | 基于卷积神经网络的大型高速回转装备多级零部件初始不平衡量逐级堆叠方法和装置 |
CN109960869A (zh) * | 2019-03-21 | 2019-07-02 | 哈尔滨工业大学 | 基于功能性滤波处理的大型高速回转装备零部件形心和惯性中心数据处理方法 |
CN109960867A (zh) * | 2019-03-21 | 2019-07-02 | 哈尔滨工业大学 | 一种基于惯性回转中心调控最优化及智能学习的大型高速回转装备多级零部件选配方法 |
CN110153664A (zh) * | 2019-03-21 | 2019-08-23 | 哈尔滨工业大学 | 基于bp神经网络的大型高速回转装备多级零部件装配方法和装置 |
CN109974749A (zh) * | 2019-04-09 | 2019-07-05 | 哈尔滨工业大学 | 一种三轴转台综合指向误差的评定方法 |
Non-Patent Citations (1)
Title |
---|
基于矢量投影的多级转子同轴度测量方法研究;孙传智;《中国博士学位论文全文数据库工程科技Ⅱ辑》;20180115(第01期);C031-28 * |
Also Published As
Publication number | Publication date |
---|---|
CN110877751A (zh) | 2020-03-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110877751B (zh) | 基于矢量投影的大型高速回转装备转动惯量堆叠方法 | |
CN110595690B (zh) | 基于形心质心重心惯性中心矢量极小化的大型高速回转装备测量与智能学习装配方法和装置 | |
CN108844524B (zh) | 一种基于振动修正的直升机旋翼桨叶运动参数测量方法 | |
CN109960867B (zh) | 一种基于惯性回转中心调控最优化及智能学习的大型高速回转装备多级零部件选配方法 | |
CN105426566B (zh) | 基于蒙特卡罗技术的大型高速回转装备形位公差分配方法 | |
CN109871947B (zh) | 基于卷积神经网络的大型高速回转装备多级零部件初始不平衡量逐级堆叠方法 | |
CN111413030B (zh) | 基于刚度矢量空间投影极大化的大型高速回转装备测量与神经网络学习调控方法及其装置 | |
CN109597357B (zh) | 一种面向叶片旋铣工艺的数控编程方法及装置 | |
CN107480356B (zh) | 基于catia和激光跟踪仪的部件设计检验一体化方法 | |
CN109190248B (zh) | 一种用于滑翔飞行器的滑翔射程解析方法及解析系统 | |
CN105426565A (zh) | 基于概率密度技术的大型高速回转装备形位公差分配方法 | |
CN110119553A (zh) | 一种止口连接的航空发动机转子零件选配优化方法 | |
CN107977491A (zh) | 一种非稳态情况下飞行器空气舵缝隙的气动热评估方法 | |
CN109255149A (zh) | 一种反击式转轮最高效率数值计算方法和系统 | |
CN113311869B (zh) | 一种多约束条件下无人机路径规划方法 | |
CN113280800A (zh) | 一种磁悬浮控制敏感陀螺角动量包络框架等效分析方法 | |
CN109960870B (zh) | 基于接触面积极大化测调的大型高速回转装备多级零部件刚度预测方法 | |
CN104915500B (zh) | 粉末注射成形涡轮及其结构优化设计方法 | |
CN111076867B (zh) | 基于质心和惯性中心同步测调的大型高速回转装备多级零部件不平衡量分配方法 | |
CN111666639B (zh) | 一种消除曲轴型转子装配的优化方法 | |
CN115979687A (zh) | 一种基于三轴转台的低成本四轴转台改装方法 | |
CN114519232B (zh) | 一种箭体运动方程系数计算方法及系统 | |
CN113110630A (zh) | 一种换电机器人举升系统节能性集成参数优化方法 | |
CN115145301A (zh) | 一种基于改进粒子群算法的多无人机任务分配方法 | |
CN106342019B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |