CN113268908A - 一种转子系统的响应求解方法和装置 - Google Patents
一种转子系统的响应求解方法和装置 Download PDFInfo
- Publication number
- CN113268908A CN113268908A CN202110564291.7A CN202110564291A CN113268908A CN 113268908 A CN113268908 A CN 113268908A CN 202110564291 A CN202110564291 A CN 202110564291A CN 113268908 A CN113268908 A CN 113268908A
- Authority
- CN
- China
- Prior art keywords
- matrix
- dynamic model
- rotor system
- modal
- equivalent
- 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
- 238000000034 method Methods 0.000 title claims abstract description 78
- 230000004044 response Effects 0.000 title claims abstract description 42
- 239000011159 matrix material Substances 0.000 claims abstract description 185
- 238000004458 analytical method Methods 0.000 claims abstract description 18
- 238000006073 displacement reaction Methods 0.000 claims description 20
- 238000013016 damping Methods 0.000 claims description 15
- 239000013598 vector Substances 0.000 claims description 15
- 230000009467 reduction Effects 0.000 claims description 12
- 230000009466 transformation Effects 0.000 claims description 7
- 230000002194 synthesizing effect Effects 0.000 claims description 6
- 238000013519 translation Methods 0.000 claims description 6
- 230000001131 transforming effect Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 abstract description 14
- 238000000926 separation method Methods 0.000 abstract description 3
- 230000008569 process Effects 0.000 description 5
- 238000005452 bending Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000010008 shearing Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000005096 rolling process 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
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Complex Calculations (AREA)
Abstract
本发明提供一种转子系统的响应求解方法和装置,方法包括:根据有限元法对转子系统建立系统动力学模型;对系统动力学模型添加等效虚拟支撑刚度矩阵,获得转子系统的系统等效动力学模型;对系统等效动力学模型进行模态分析,获得模态矩阵;对模态矩阵进行降阶,并基于降阶后的模态矩阵获得降阶的系统等效动力学模型;根据状态空间法对降阶的系统等效动力学模型进行响应求解,获得转子系统的动力学响应。本发明的方法通过在系统动力学模型中添加等效虚拟支撑,从而使确保算法稳定的模态分离,再通过对模态矩阵进行降阶,获得降阶的系统等效动力学模型,以提高转子系统响应的计算效率和计算精度。本发明的装置用于实现上述方法。
Description
技术领域
本发明涉及转子系统技术领域,尤其涉及一种转子系统的响应求解方法和装置。
背景技术
为了提高功率密度和能量转化效率,转子的工作转速往往无法避免地高于弯曲临界转速。转子在转速上升的过程中,需要超越系统的刚性和低阶弯曲临界转速,此时由于转子存在不平衡质量,通常会产生剧烈的共振,因此对转子系统动力特性的分析显得十分重要。
当针对无约束自由支撑的转子系统进行模态分析时,有限元方法由于节点与自由度数量较多导致各结构参数矩阵复杂,并导致矩阵求逆及求解特征值产生数值误差,最终将导致矩阵模态求解不准确等问题。
发明内容
为了解决上述技术问题的至少一个,本发明提供了一种转子系统的响应求解方法和装置。
本发明的技术方案是这样实现的:
第一方面,本发明的实施例提供了一种转子系统的响应求解方法,包括:
根据有限元法对转子系统建立系统动力学模型;
对所述系统动力学模型添加等效虚拟支撑刚度矩阵,获得转子系统的系统等效动力学模型;
对所述系统等效动力学模型进行模态分析,获得模态矩阵;
对所述模态矩阵进行降阶,并基于降阶后的模态矩阵获得降阶的系统等效动力学模型;
根据状态空间法对降阶的系统等效动力学模型进行响应求解,获得转子系统的动力学响应。
作为可选的实施方式,所述模态矩阵为r列特征向量组成的r×r阶矩阵;r为转子系统动力学模型的自由度数量,r>1;
对所述模态矩阵进行降阶,包括:
将所述r×r阶矩阵变换成r×s阶矩阵,获得降阶后的模态矩阵,其中s<r。
作为可选的实施方式,所述基于降阶后的模态矩阵获得降阶的系统等效动力学模型,包括:
对系统等效动力学模型做模态变换,获得系统等效动力学模态方程;
将降阶后的模态矩阵代入所述系统等效动力学模态方程,获得所述降阶的系统等效动力学模型。
作为可选的实施方式,所述根据有限元法对转子系统建立系统动力学模型,包括:
以转子的轴线为z轴,以转子两个相互垂直的径向分别为x轴和y轴,建立坐标系;
基于所述坐标系,采用欧拉梁单元对转子系统进行有限元建模,获得系统动力学模型,所述系统动力学模型具有n个节点和r个自由度,并且r=4n,n≥1。
作为可选的实施方式,所述转子系统包括均质等截面弹性轴段和刚性圆盘单元;
所述根据有限元法对转子系统建立系统动力学模型,包括:
基于所述均质等截面弹性轴段,根据有限元法建立转子系统的数个弹性轴段单元;每个所述弹性轴段单元的两个端面分别为转子系统的两个节点;
基于每个所述弹性轴段单元的节点力与节点位移之间的关系,建立每个弹性轴段单元的基本运动微分方程;
综合全部所述弹性轴端单元的基本运动微分方程,获得转子系统在x轴和y轴方向的系统运动微分方程;
在确定转子系统的其中任意一个或多个节点处存在所述刚性圆盘单元时,将所述刚性圆盘单元的质量矩阵对应添加至所述系统运动微分方程;以及
在确定转子系统的其中任意一个或多个节点处存在外力或外力矩时,在所述系统运动微分方程中添加外力或外力矩,构成转子系统的动力学方程,建立系统动力学模型。
作为可选的实施方式,对所述系统动力学模型添加等效虚拟支撑刚度矩阵,获得转子系统的系统等效动力学模型,包括:
根据系统动力学模型,获取转子系统的质量矩阵、阻尼矩阵和刚度矩阵;
基于所述质量矩阵、阻尼矩阵和刚度矩阵,获得转子系统的动力学方程如下:
在所述动力学方程的左侧添加等效虚拟支撑矩阵,在所述动力学方程的右侧添加对应的载荷,获得转子系统的系统等效动力学模型;
其中,M为转子系统的质量矩阵;D为转子系统的阻尼矩阵;K为转子系统的刚度矩阵;f为外力或外力矩;q为广义位移矢量;
q=[x1,θy1,x2,θy2,...,xn,θyn,y1,-θx1,y2,-θx2,...,yn,-θxn]T;
其中,n为系统动力学模型的节点数量,n>1;xi为系统动力学模型的第i个节点处x轴方向的平动位移;yi为第i个节点处y轴方向的平动位移;θxi为第i个节点处绕x轴的偏转角;θyi为第i个节点处绕y轴的偏转角;1≤i≤n。
作为可选的实施方式,所述系统等效动力学模型的动力学方程为:
其中,Kev为r×r阶的等效虚拟支撑刚度矩阵,r为转子系统动力学模型的自由度数量。
第二方面,本发明的实施例提供一种转子系统的响应求解装置,包括:
动力学模型建立模块,用于根据有限元法对转子系统建立系统动力学模型;
等效动力学模型获取模块,用于对所述系统动力学模型添加等效虚拟支撑刚度矩阵,获得转子系统的系统等效动力学模型;
模态分析模块,用于对所述系统等效动力学模型进行模态分析,获得模态矩阵;
等效动力学模型降阶模块,用于对所述模态矩阵进行降阶,并基于降阶后的模态矩阵获得降阶的系统等效动力学模型;
响应求解模块,用于根据状态空间法对降阶的系统等效动力学模型进行响应求解,获得转子系统的动力学响应。
作为可选的实施方式,所述模态矩阵为r列特征向量组成的r×r阶矩阵;r为转子系统动力学模型的自由度数量,r>1;
对所述模态矩阵进行降阶,包括:
将所述r×r阶矩阵变换成r×s阶矩阵,获得降阶后的模态矩阵,其中s<r。
作为可选的实施方式,所述等效动力学模型降阶模块,具体用于:
对所述模态矩阵进行降阶;
对系统等效动力学模型做模态变换,获得系统等效动力学模态方程;
将降阶后的模态矩阵代入所述系统等效动力学模态方程,获得降阶的系统等效动力学模型。
作为可选的实施方式,所述转子系统包括均质等截面弹性轴段和刚性圆盘单元;
所述动力学模型建立模块,具体用于:
基于所述均质等截面弹性轴段,根据有限元法建立转子系统的数个弹性轴段单元;每个所述弹性轴段单元的两个端面分别为转子系统的两个节点;
基于每个所述弹性轴段单元的节点力与节点位移之间的关系,建立每个弹性轴段单元的基本运动微分方程;
综合全部所述弹性轴端单元的基本运动微分方程,获得转子系统在x轴和y轴方向的系统运动微分方程;
在确定转子系统的其中任意一个或多个节点处存在所述刚性圆盘单元时,将所述刚性圆盘单元的质量矩阵对应添加至所述系统运动微分方程;以及
在确定转子系统的其中任意一个或多个节点处存在外力或外力矩时,在所述系统运动微分方程中添加外力或外力矩,构成转子系统的动力学方程,建立系统动力学模型。
本发明的实施例与现有技术相比,至少具有如下优点:
本发明通过有限元法转子系统的有限元系统动力学模型,然后通过在系统动力学模型中添加等效虚拟支撑,获得等效虚拟支撑的等效动力学模型,从而使确保算法稳定的模态分离,再通过对模态矩阵进行降阶,获得降阶的系统等效动力学模型,以提高转子系统响应的计算效率和计算精度。
附图说明
附图示出了本发明的示例性实施方式,并与其说明一起用于解释本发明的原理,其中包括了这些附图以提供对本发明的进一步理解,并且附图包括在本说明书中并构成本说明书的一部分。
图1为本发明的转子系统的响应求解方法的流程示意图;
图2为本发明的建立系统动力学模型的流程示意图;
图3为本发明的弹性轴段的质量矩阵组装示意图;
图4为本发明的弹性轴段的刚度矩阵组装示意图;
图5为本发明考虑刚性圆盘单元的转子系统的质量矩阵组装示意图;
图6为本发明的转子系统的响应求解装置的结构示意图。
具体实施方式
下面结合附图和实施方式对本发明作进一步的详细说明。可以理解的是,此处所描述的具体实施方式仅用于解释相关内容,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分。
需要说明的是,在不冲突的情况下,本发明中的实施方式及实施方式中的特征可以相互组合。下面将参考附图并结合实施方式来详细说明本发明。
需要说明的是,文中的步骤编号,仅为了方便具体实施例的解释,不作为限定步骤执行先后顺序的作用。本实施例提供的方法可以由相关的服务器执行,且下文均以服务器或计算机等电子设备作为执行主体为例进行说明。
实施例一
目前,柔性转子系统动力特性的分析方法主要有传递矩阵法和有限元法。在传递矩阵法中,传递矩阵的阶数不随转子系统自由度增大而增大,因而编程简单,运算速度快。因此,在很长一段时间里,传递矩阵法在柔性转子系统动力学的研究中占主导地位,但这种方法用于求解高速大型转子系统的动力学问题时,有可能出现数值不稳定现象,而且传递矩阵法不易解决转子系统的非线性问题。
有限元法分析转子动力学问题,起初只用于研究转子的弯曲振动问题,随着研究的逐步深入,转子的有限元模型也不断得到完善,在模型中逐渐包括了转动惯量、陀螺力矩、轴向载荷、外内阻以及剪切变形等因素的影响。有限元法比传递矩阵法更为精确的建模方法,但是,其计算矩阵的阶数按转子节点数量的增加而剧烈增长,计算复杂度较大。随着计算机技术的迅猛发展,有限元法越来越体现出其优势。在近几十年来,用有限元分析在研究转子的临界转速、不平衡响应以及稳定性等方面取得了很大的进步,已经成为柔性转子系统动力学常用的分析工具。
由于有限元模型可以用很多的自由度来描述复杂结构,其计算复杂度较大,当针对无约束自由支撑的转子系统进行模态分析时,有限元方法由于节点与自由度数量较多导致各结构参数矩阵复杂,并导致矩阵求逆及求解特征值产生数值误差,该误差将导致矩阵模态求解不准确,从而无法有效降阶。本方法通过采用等效虚拟支撑作用于转子系统,实现转子系统的高精度降阶,提升转子系统响应的数值计算效率。
参照图1,本发明的实施例提供了一种转子系统的响应求解方法,包括:
S1,根据有限元法对转子系统建立系统动力学模型;
S2,对系统动力学模型添加等效虚拟支撑刚度矩阵,获得转子系统的系统等效动力学模型;
S3,对系统等效动力学模型进行模态分析,获得模态矩阵;
S4,对模态矩阵进行降阶,并基于降阶后的模态矩阵获得降阶的系统等效动力学模型;
S5,根据状态空间法对降阶的系统等效动力学模型进行响应求解,获得转子系统的动力学响应。
在本实施例中,转子系统优选为柔性转子系统,本发明通过有限元法建立转子系统的有限元系统动力学模型,然后通过在系统动力学模型中添加等效虚拟支撑,获得等效虚拟支撑的等效动力学模型,从而使确保算法稳定的模态分离,再通过对模态矩阵进行降阶,获得降阶的系统等效动力学模型,以提高转子系统响应的计算效率和计算精度。
作为本实施例的优选实施方式,根据有限元法对转子系统建立系统动力学模型,包括:
以转子的轴线为z轴,以转子两个相互垂直的径向分别为x轴和y轴,建立oxyz坐标系;基于坐标系,采用欧拉梁(Euler)单元对转子系统进行有限元建模,获得系统动力学模型,由于转子上任一截面的位移可由截面轴心沿x轴和y轴的平移坐标绕x轴和y轴的截面的偏转角θy及θx,因此每个节点处具有4个自由度,假设根据有限元法建立的系统动力学模型共有n个节点,r个自由度,那么r=4n,n≥1。本实施例采用欧拉梁单元对转子系统进行有限元建模,建模时可不考虑转子的陀螺效应。
组成转子系统的典型单元包括均质等截面弹性轴段、刚性圆盘单元、支撑等;
参照图2,在步骤S1中,根据有限元法对转子系统建立系统动力学模型,具体包括:
S11,建立弹性轴段单元:
基于均质等截面弹性轴段,根据有限元法建立转子系统的数个弹性轴段单元;每个弹性轴段单元的两个端面分别为转子系统的两个节点;基于每个弹性轴段单元的节点力与节点位移之间的关系,建立每个弹性轴段单元的基本运动微分方程;
本实施例方式中,按从左至右顺序排列,其中弹性轴段单元的左端面为第i节点,右端面为第i+1节点。弹性轴段单元的广义坐标决定于两个端面的位移。每个端面的位置可由其端面轴心的坐标x及y和偏转角θy及θx来表示,即:
在不考虑轴段剪切和扭转变形影响的条件下,一个弹性轴段单元在x和y两个方向上的运动微分方程分别为
其中,fx与fy分别为弹性轴段两个端面节点处在x轴及y轴两个方向上所受的作用力和力矩的广义外力矢量,包括由不平衡引起的外力。ux与uy为弹性轴段两个端面在x轴和y轴两个方向上响应位移的广义坐标。mx、dx、kx和my、dy、ky分别为转子在x轴和y轴两个方向上的质量矩阵、阻尼矩阵和刚度矩阵。由于转子在径向一般是对称的,因此以上矩阵在x轴和y轴两个方向上相同,可用m、d、k表示。对于单个轴段在一个坐标方向内的运动方程而言,其m、d、k矩阵都是4×4的矩阵,具体表达式为:
其中,μ为材料单位长质量(线密度),l为轴段长度,EI为等截面轴段的弯曲刚度,阻尼矩阵d为零矩阵。
S12,将弹性轴段单元的基本运动微分方程组成转子系统的系统运动微分方程:
综合全部弹性轴端单元的基本运动微分方程,获得转子系统的弹性轴段在x轴和y轴方向的系统运动微分方程;
对于具有n个节点,n-1个轴段构成的转子系统,综合各轴段的运动方程,就可以得到整个转子系统在x及y方向上的运动微分方程分别为
其中,M0为弹性轴段单元的整体质量矩阵,D0为弹性轴段单元的整体阻尼矩阵,K0为弹性轴段单元的整体刚度矩阵,F0为各个弹性轴段端面节点处在x轴及y轴两个方向上所受的作用力和力矩的广义外力矢量。
图3为本发明的弹性轴段的质量矩阵组装示意图,其中,M0为弹性轴段单元的整体质量矩阵;图3中,小方块表示相应单个轴段的4×4的矩阵m;每个轴段的矩阵m按图3所示以对角线排列,并且在相邻矩阵的2×2部分重叠相加;
D0为弹性轴段单元的整体阻尼矩阵,对于本实施例,该矩阵为零矩阵。
图4为本发明的弹性轴段的刚度矩阵组装示意图,其中,K0为弹性轴段单元的整体刚度矩阵,图4中,小方块表示相应单个轴段的4×4的矩阵k。每个轴段的矩阵k按图4所示以对角线排列,并且在相邻矩阵的2×2部分重叠相加。
S13,建立刚性圆盘单元:
设刚性圆盘单元的质量为m、过轴心的直径转动惯量和极转动惯量分别为Jd和Jp。由于刚性圆盘不考虑刚度特性,因此刚性圆盘单元的运动微分方程无刚度项,刚性圆盘单元的质量矩阵为:
在确定转子系统的其中任意一个或多个节点处存在刚性圆盘单元时,将刚性圆盘单元的质量矩阵对应添加至系统运动微分方程;例如:如果转子在节点a处存在一个刚性圆盘,则需要在原来未考虑刚性圆盘时的整个转子系统运动方程中的质量矩阵M0的基础上作如图5所示的变化,在相应位置叠加2×2矩阵。其中,质量矩阵M0为弹性轴段的质量矩阵,即未计圆盘时整个转子的质量矩阵,i=2a-1,j=2n+2a-1。若转子有多个圆盘,矩阵M就在各自相应位置多次叠加。
S14,添加外力或外力矩:
外力或外力矩包括支撑力与不平衡力等。其中支撑力施加在转子系统中对应的节点位置,通常是轴承、滑动轴承或电磁轴承支撑等形式。当考虑不平衡激励对转子系统动力特性的影响时,转子的不平衡会产生与转速同频的离心力,导致转子同频振动,称之为不平衡力干扰fu(t)。根据力学原理,在恒定转速下,某节点位置不平衡力干扰可表示为
式中,ω为转子旋转角速度,m为不平衡质量,θ0为初始不平衡相角,e为实际偏心距,即转子当前的旋转轴心与质心的距离。
在确定转子系统的其中任意一个或多个节点处存在外力或外力矩时,在系统运动微分方程中添加外力或外力矩,构成转子系统的动力学方程,建立系统动力学模型,然后获取转子系统的质量矩阵、阻尼矩阵和刚度矩阵;并基于质量矩阵、阻尼矩阵和刚度矩阵,获得转子系统的动力学方程如下:
其中,M为转子系统的质量矩阵;D为转子系统的阻尼矩阵,可采用比例阻尼表示D=αM+βK,;K为转子系统的刚度矩阵;f为外力或外力矩;q为广义位移矢量;
q=[x1,θy1,x2,θy2,...,xn,θyn,y1,-θx1,y2,-θx2,...,yn,-θxn]T;
其中,n为系统动力学模型的节点数量,n>1;xi为系统动力学模型的第i个节点处x轴方向的平动位移;yi为第i个节点处y轴方向的平动位移;θxi为第i个节点处绕x轴的偏转角;θyi为第i个节点处绕y轴的偏转角;1≤i≤n。
作为本实施例的优选技术方案,对系统动力学模型添加等效虚拟支撑刚度矩阵,获得转子系统的系统等效动力学模型,包括:
在动力学方程的左侧添加等效虚拟支撑矩阵,
在动力学方程的右侧添加对应的载荷,获得转子系统的系统等效动力学模型;同时为了保证系统的一致性,系统等效动力学模型的动力学方程为:
其中,Kev为r×r阶的等效虚拟支撑刚度矩阵,依据广义位移矢量q的节点分布规律,确认实际支撑位置对应的节点序号,例如对应的第i、j、c、d个节点位置施加了实际的外部约束,则在虚拟支撑刚度矩阵Kev中对应元素kii、kjj、kcc、kdd为设置的等效虚拟刚度,其他元素数值为0。其中,效虚拟刚度数值不宜过大或过小,数值过大对系统临界转速影响较大,数值过小则其作用不明显,使得模态求解过程中数值误差仍然较大,等效虚拟刚度可以为正值,也可以为负值。
本发明通过在模型中添加等效虚拟支撑,在此基础上对转子系统进行模态分析与降阶,对于添加等效虚拟支撑的系统,其模态降阶相比于无支撑约束系统更加准确,对于无约束转子系统,通常其零频对应着刚体模态,此时系统矩阵特征值理论上应该为零,但是由于矩阵求逆、求特征值过程中的数值误差,导致系统矩阵特征值及特征向量计算存在误差,从而导致模态矩阵求解不准确。而本发明的实施例通过采用等效虚拟支撑方法后,所计算的特征值并不为零,此时产生的数值误差相对较小,能够较好地提升后续模态分析的准确性,为后续模态降阶与响应求解奠定了良好基础。
本实例中,对系统等效动力学模型进行模态分析,获得模态矩阵,包括:
求解模态矩阵,设系统矩阵S=M-1(K+Kev),则系统矩阵S的特征值λi的开方即为对应的各阶临界转速,特征向量Φi组成的矩阵即为模态矩阵Tm,且Tm=[Φ1,Φ2,......,Φr,]。
利用模态矩阵对等效系统作模态变换
q=Tmqm;
其中,qm为模态坐标,则系统在模态坐标中可表示为:
其中,Mm、Dm、Kem、和Kvm分别为模态坐标下的模态质量矩阵、模态阻尼矩阵、模态刚度矩阵和模态等效虚拟支撑刚度矩阵,fm为模态坐标下转子系统所受的作用力和力矩的广义外力矢量。各模态矩阵和向量求解方法为:
本实施方式中,全阶的模态矩阵为r列特征向量组成的r×r阶矩阵;通常转子系统中高阶模态对系统影响较小,可以忽略,因此忽略模态矩阵Tm中对应的高阶模态,因此,对模态矩阵进行降阶,包括:
将r×r阶矩阵变换成r×s阶矩阵,获得降阶后的模态矩阵,其中s<r。
作为优选的,基于降阶后的模态矩阵获得降阶的系统等效动力学模型,包括:
对系统等效动力学模型做模态变换,获得系统等效动力学模态方程;将降阶后的模态矩阵代入系统等效动力学模态方程,获得降阶的系统等效动力学模型。
作为本实施例的优选实施方式,根据状态空间法对降阶的系统等效动力学模型进行响应求解,获得转子系统的动力学响应,包括:
定义状态变量:
利用状态空间方法,转子系统状态方程为:
其中,Am为系统矩阵,Bm1、Bm2为输入矩阵,y为系统位移输出,Cm为输出矩阵。各矩阵表达如下:
Cm=[Ts Tm O];
其中,O为零矩阵,I为单位矩阵,Ts为定义传感器位置矩阵。通过求解该状态方程可以确定转子各节点对应的响应。
综上,本发明的实施例的方法针对滚动轴承、滑动轴承、电磁轴承等形式支撑的转子系统,首先建立转子系统有限元动力学模型,然后在转子系统动力学模型中添加等效虚拟支撑,进而对等效虚拟支撑的转子系统进行模态分析与降阶,最后完成降阶后的等效虚拟支撑转子系统的响应求解。本发明能够在保证计算精度的同时,实现转子系统阶次的合理简化,显著提升计算效率。
实施例二
参照图6,本发明的实施例还提供一种转子系统的响应求解装置,包括:
动力学模型建立模块,用于根据有限元法对转子系统建立系统动力学模型;
等效动力学模型获取模块,用于对系统动力学模型添加等效虚拟支撑刚度矩阵,获得转子系统的系统等效动力学模型;
模态分析模块,用于对系统等效动力学模型进行模态分析,获得模态矩阵;
等效动力学模型降阶模块,用于对模态矩阵进行降阶,并基于降阶后的模态矩阵获得降阶的系统等效动力学模型;
响应求解模块,用于根据状态空间法对降阶的系统等效动力学模型进行响应求解,获得转子系统的动力学响应。
作为可选的实施方式,模态矩阵为r列特征向量组成的r×r阶矩阵;r为转子系统动力学模型的自由度数量,r>1;
对模态矩阵进行降阶,包括:
将r×r阶矩阵变换成r×s阶矩阵,获得降阶后的模态矩阵,其中s<r。
作为可选的实施方式,等效动力学模型降阶模块,具体用于:
对模态矩阵进行降阶;
对系统等效动力学模型做模态变换,获得系统等效动力学模态方程;
将降阶后的模态矩阵代入系统等效动力学模态方程,获得降阶的系统等效动力学模型。
作为可选的实施方式,转子系统包括均质等截面弹性轴段和刚性圆盘单元;
动力学模型建立模块,具体用于:
基于均质等截面弹性轴段,根据有限元法建立转子系统的数个弹性轴段单元;每个弹性轴段单元的两个端面分别为转子系统的两个节点;
基于每个弹性轴段单元的节点力与节点位移之间的关系,建立每个弹性轴段单元的基本运动微分方程;
综合全部弹性轴端单元的基本运动微分方程,获得转子系统在x轴和y轴方向的系统运动微分方程;
在确定转子系统的其中任意一个或多个节点处存在刚性圆盘单元时,将刚性圆盘单元的质量矩阵对应添加至系统运动微分方程;以及
在确定转子系统的其中任意一个或多个节点处存在外力或外力矩时,在系统运动微分方程中添加外力或外力矩,构成转子系统的动力学方程,建立系统动力学模型。
本实施例的装置中的各模块的原理和功能与实施例一中的一致,本实施例不再重复描述。
应理解的是,本发明的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中,多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。上述实施例方法的全部或部分步骤是可以通过程序来指令相关的硬件完成,该程序可以存储于一种计算机可读存储介质中,该程序在执行时,包括方法实施例的步骤之一或其组合。
此外,在本发明各个实施例中的各功能单元可以集成在一个处理模块中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现,也可以采用软件功能模块的形式实现。上述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时,也可以存储在一个计算机可读存储介质中。该存储介质可以是只读存储器,磁盘或光盘等。
本发明的流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为,表示包括一个或更多个用于实现特定逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分。并且本发明的优选实施方式的范围包括另外的实现,其中可以不按所示出或讨论的顺序,包括根据所涉及的功能按基本同时的方式或按相反的顺序,来执行功能。
本发明在流程图中表示或在此以其他方式描述的逻辑和/或步骤,例如,可以被认为是用于实现逻辑功能的可执行指令的定序列表,可以具体实现在任何计算机可读介质中,以供指令执行系统、装置或设备(如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用,或结合这些指令执行系统、装置或设备而使用。
本领域的技术人员应当理解,上述实施方式仅仅是为了清楚地说明本发明,而并非是对本发明的范围进行限定。对于所属领域的技术人员而言,在上述公开的基础上还可以做出其它变化或变型,并且这些变化或变型仍处于本发明的范围内。
Claims (10)
1.一种转子系统的响应求解方法,其特征在于,包括:
根据有限元法对转子系统建立系统动力学模型;
对所述系统动力学模型添加等效虚拟支撑刚度矩阵,获得转子系统的系统等效动力学模型;
对所述系统等效动力学模型进行模态分析,获得模态矩阵;
对所述模态矩阵进行降阶,并基于降阶后的模态矩阵获得降阶的系统等效动力学模型;
根据状态空间法对降阶的系统等效动力学模型进行响应求解,获得转子系统的动力学响应。
2.如权利要求1所述的方法,其特征在于,
所述模态矩阵为r列特征向量组成的r×r阶矩阵;r为转子系统动力学模型的自由度数量,r>1;
对所述模态矩阵进行降阶,包括:
将所述r×r阶矩阵变换成r×s阶矩阵,获得降阶后的模态矩阵,其中s<r。
3.如权利要求1或2所述的方法,其特征在于,
所述基于降阶后的模态矩阵获得降阶的系统等效动力学模型,包括:
对系统等效动力学模型做模态变换,获得系统等效动力学模态方程;
将降阶后的模态矩阵代入所述系统等效动力学模态方程,获得所述降阶的系统等效动力学模型。
4.如权利要求1所述的方法,其特征在于,
所述转子系统包括均质等截面弹性轴段和刚性圆盘单元;
所述根据有限元法对转子系统建立系统动力学模型,包括:
基于所述均质等截面弹性轴段,根据有限元法建立转子系统的数个弹性轴段单元;每个所述弹性轴段单元的两个端面分别为转子系统的两个节点;
基于每个所述弹性轴段单元的节点力与节点位移之间的关系,建立每个弹性轴段单元的基本运动微分方程;
综合全部所述弹性轴端单元的基本运动微分方程,获得转子系统在x轴和y轴方向的系统运动微分方程;
在确定转子系统的其中任意一个或多个节点处存在所述刚性圆盘单元时,将所述刚性圆盘单元的质量矩阵对应添加至所述系统运动微分方程;以及
在确定转子系统的其中任意一个或多个节点处存在外力或外力矩时,在所述系统运动微分方程中添加外力或外力矩,构成转子系统的动力学方程,建立系统动力学模型。
5.如权利要求1所述的方法,其特征在于,对所述系统动力学模型添加等效虚拟支撑刚度矩阵,获得转子系统的系统等效动力学模型,包括:
根据系统动力学模型,获取转子系统的质量矩阵、阻尼矩阵和刚度矩阵;
基于所述质量矩阵、阻尼矩阵和刚度矩阵,获得转子系统的动力学方程如下:
在所述动力学方程的左侧添加等效虚拟支撑矩阵,在所述动力学方程的右侧添加对应的载荷,获得转子系统的系统等效动力学模型;
其中,M为转子系统的质量矩阵;D为转子系统的阻尼矩阵;K为转子系统的刚度矩阵;f为外力或外力矩;q为广义位移矢量;
q=[x1,θy1,x2,θy2,...,xn,θyn,y1,-θx1,y2,-θx2,...,yn,-θxn]T;
其中,n为系统动力学模型的节点数量,n>1;xi为系统动力学模型的第i个节点处x轴方向的平动位移;yi为第i个节点处y轴方向的平动位移;θxi为第i个节点处绕x轴的偏转角;θyi为第i个节点处绕y轴的偏转角;1≤i≤n。
7.一种转子系统的响应求解装置,其特征在于,包括:
动力学模型建立模块,用于根据有限元法对转子系统建立系统动力学模型;
等效动力学模型获取模块,用于对所述系统动力学模型添加等效虚拟支撑刚度矩阵,获得转子系统的系统等效动力学模型;
模态分析模块,用于对所述系统等效动力学模型进行模态分析,获得模态矩阵;
等效动力学模型降阶模块,用于对所述模态矩阵进行降阶,并基于降阶后的模态矩阵获得降阶的系统等效动力学模型;
响应求解模块,用于根据状态空间法对降阶的系统等效动力学模型进行响应求解,获得转子系统的动力学响应。
8.如权利要求7所述的装置,其特征在于,
所述模态矩阵为r列特征向量组成的r×r阶矩阵;r为转子系统动力学模型的自由度数量,r>1;
对所述模态矩阵进行降阶,包括:
将所述r×r阶矩阵变换成r×s阶矩阵,获得降阶后的模态矩阵,其中s<r。
9.如权利要求7或8所述的装置,其特征在于,所述等效动力学模型降阶模块,具体用于:
对所述模态矩阵进行降阶;
对系统等效动力学模型做模态变换,获得系统等效动力学模态方程;
将降阶后的模态矩阵代入所述系统等效动力学模态方程,获得所述降阶的系统等效动力学模型。
10.如权利要求7所述的装置,其特征在于,
所述转子系统包括均质等截面弹性轴段和刚性圆盘单元;
所述动力学模型建立模块,具体用于:
基于所述均质等截面弹性轴段,根据有限元法建立转子系统的数个弹性轴段单元;每个所述弹性轴段单元的两个端面分别为转子系统的两个节点;
基于每个所述弹性轴段单元的节点力与节点位移之间的关系,建立每个弹性轴段单元的基本运动微分方程;
综合全部所述弹性轴端单元的基本运动微分方程,获得转子系统在x轴和y轴方向的系统运动微分方程;
在确定转子系统的其中任意一个或多个节点处存在所述刚性圆盘单元时,将所述刚性圆盘单元的质量矩阵对应添加至所述系统运动微分方程;以及
在确定转子系统的其中任意一个或多个节点处存在外力或外力矩时,在所述系统运动微分方程中添加外力或外力矩,构成转子系统的动力学方程,建立系统动力学模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110564291.7A CN113268908B (zh) | 2021-05-24 | 2021-05-24 | 一种转子系统的响应求解方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110564291.7A CN113268908B (zh) | 2021-05-24 | 2021-05-24 | 一种转子系统的响应求解方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113268908A true CN113268908A (zh) | 2021-08-17 |
CN113268908B CN113268908B (zh) | 2024-03-01 |
Family
ID=77232357
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110564291.7A Active CN113268908B (zh) | 2021-05-24 | 2021-05-24 | 一种转子系统的响应求解方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113268908B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114611358A (zh) * | 2022-03-15 | 2022-06-10 | 清华大学 | 利用等效梁单元建模的阶梯转子动力学分析方法 |
CN118378388A (zh) * | 2024-06-24 | 2024-07-23 | 中国航空工业集团公司金城南京机电液压工程研究中心 | 一种多跨转子系统临界转速优化方法和装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20140065692A (ko) * | 2012-11-20 | 2014-05-30 | 국방과학연구소 | 등가 유한요소 모델링을 이용한 감쇠평판구조물 해석방법 |
CN111209639A (zh) * | 2020-02-17 | 2020-05-29 | 合肥工业大学 | 一种叶轮-轴承-转子系统的高效定量建模方法 |
CN112528411A (zh) * | 2020-12-10 | 2021-03-19 | 中国运载火箭技术研究院 | 一种基于模态减缩的几何非线性结构噪声振动响应计算方法 |
EP3809305A1 (de) * | 2019-10-15 | 2021-04-21 | Siemens Aktiengesellschaft | Verfahren und vorrichtung zum erzeugen einer digitalen repräsentation eines technischen gebildes und entsprechendes computerprogrammprodukt |
-
2021
- 2021-05-24 CN CN202110564291.7A patent/CN113268908B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20140065692A (ko) * | 2012-11-20 | 2014-05-30 | 국방과학연구소 | 등가 유한요소 모델링을 이용한 감쇠평판구조물 해석방법 |
EP3809305A1 (de) * | 2019-10-15 | 2021-04-21 | Siemens Aktiengesellschaft | Verfahren und vorrichtung zum erzeugen einer digitalen repräsentation eines technischen gebildes und entsprechendes computerprogrammprodukt |
CN111209639A (zh) * | 2020-02-17 | 2020-05-29 | 合肥工业大学 | 一种叶轮-轴承-转子系统的高效定量建模方法 |
CN112528411A (zh) * | 2020-12-10 | 2021-03-19 | 中国运载火箭技术研究院 | 一种基于模态减缩的几何非线性结构噪声振动响应计算方法 |
Non-Patent Citations (2)
Title |
---|
彭航;秦剑波;: "基于有限元法的高速测扭器动力学分析", 航空科学技术, no. 03 * |
骞朋波;钱林方;尹晓春;: "基于降阶模型的弹塑性瞬态响应求解", 机械工程学报, no. 05 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114611358A (zh) * | 2022-03-15 | 2022-06-10 | 清华大学 | 利用等效梁单元建模的阶梯转子动力学分析方法 |
CN114611358B (zh) * | 2022-03-15 | 2023-06-16 | 清华大学 | 利用等效梁单元建模的阶梯转子动力学分析方法 |
CN118378388A (zh) * | 2024-06-24 | 2024-07-23 | 中国航空工业集团公司金城南京机电液压工程研究中心 | 一种多跨转子系统临界转速优化方法和装置 |
CN118378388B (zh) * | 2024-06-24 | 2024-09-06 | 中国航空工业集团公司金城南京机电液压工程研究中心 | 一种多跨转子系统临界转速优化方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN113268908B (zh) | 2024-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zou et al. | Coupled longitudinal-transverse dynamics of a marine propulsion shafting under superharmonic resonances | |
Banerjee | Dynamic stiffness formulation and free vibration analysis of centrifugally stiffened Timoshenko beams | |
Genta | Whirling of unsymmetrical rotors: a finite element approach based on complex co-ordinates | |
CN113268908B (zh) | 一种转子系统的响应求解方法和装置 | |
Hosseini-Hashemi et al. | Accurate natural frequencies and critical speeds of a rotating functionally graded moderately thick cylindrical shell | |
CN108776734B (zh) | 一种螺栓连接鼓筒转子结构的响应特性分析方法 | |
Cai et al. | Model study and active control of a rotating flexible cantilever beam | |
CN113190930B (zh) | 柔性轴-盘-壳连接转子系统的耦合动力学建模与分析方法 | |
CN112149245B (zh) | 一种柔性轴-盘系统耦合动力学建模及分析方法 | |
Kaya et al. | Energy expressions and free vibration analysis of a rotating uniform Timoshenko beam featuring bending—torsion coupling | |
CN111382503A (zh) | 弹性支承下旋转的柔性圆环的振动分析方法及系统 | |
Pan et al. | Coupled dynamic modeling and analysis of the single gimbal control moment gyroscope driven by ultrasonic motor | |
Patil et al. | A scalable time-parallel solution of periodic rotor dynamics in X3D | |
Zhou et al. | Three-dimensional dynamics analysis of rotating functionally gradient beams based on Timoshenko beam theory | |
Klanner et al. | Quasi-analytical solutions for the whirling motion of multi-stepped rotors with arbitrarily distributed mass unbalance running in anisotropic linear bearings | |
Chao et al. | Dynamic analysis of the optical disk drives equipped with an automatic ball balancer with consideration of torsional motions | |
Chen et al. | New analytic method for free torsional vibration analysis of a shaft with multiple disks and elastic supports | |
CN114577397B (zh) | 一种高速永磁电机转子动平衡方法及系统 | |
CN112556931B (zh) | 基于粒子群算法的高速轴承转子系统模态动平衡方法 | |
CN114722673A (zh) | 一种电磁轴承转子系统的转子位移分析方法及装置 | |
KOPMAZ et al. | On the eigenfrequencies of a flexible arm driven by a flexible shaft | |
Chen et al. | Dynamic analysis of a helical geared rotor-bearing system with composite rotating shafts | |
Lin et al. | Prediction of vibration of rotating damped beams with arbitrary pretwist | |
Aleyaasin et al. | Multivariable hybrid models for rotor-bearing systems | |
Cruz et al. | Experimental design and analysis of a gyroelastic beam |
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 |