CN112556931B - 基于粒子群算法的高速轴承转子***模态动平衡方法 - Google Patents
基于粒子群算法的高速轴承转子***模态动平衡方法 Download PDFInfo
- Publication number
- CN112556931B CN112556931B CN202011319498.XA CN202011319498A CN112556931B CN 112556931 B CN112556931 B CN 112556931B CN 202011319498 A CN202011319498 A CN 202011319498A CN 112556931 B CN112556931 B CN 112556931B
- Authority
- CN
- China
- Prior art keywords
- correction
- rotor system
- bearing rotor
- balance
- particle swarm
- 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
- 239000002245 particle Substances 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000012937 correction Methods 0.000 claims abstract description 55
- 238000005457 optimization Methods 0.000 claims abstract description 12
- 238000004458 analytical method Methods 0.000 claims abstract description 8
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000013016 damping Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 12
- 239000000306 component Substances 0.000 description 5
- 238000012360 testing method Methods 0.000 description 4
- 238000013178 mathematical model Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000004323 axial length Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M1/00—Testing static or dynamic balance of machines or structures
- G01M1/30—Compensating imbalance
- G01M1/32—Compensating imbalance by adding material to the body to be tested, e.g. by correcting-weights
-
- 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
- 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/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Computer Hardware Design (AREA)
- Artificial Intelligence (AREA)
- Mathematical Optimization (AREA)
- Biophysics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pure & Applied Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Computer Graphics (AREA)
- Testing Of Balance (AREA)
Abstract
本发明提供一种基于粒子群算法的高速轴承转子***模态动平衡方法。首先建立轴承转子***的三维模型,然后导入ANSYS软件中进行模态分析,确定出平衡校正面的位置,其次根据模态平衡法的平衡条件建立***的固有不平衡量与平衡校正面上需要添加的校正质量之间的函数关系,最后将需要添加的校正质量作为一个粒子,所有的校正质量作为一个粒子群,采用粒子群算法求解每个校正质量的最优值,本发明可以有效提高转子的平衡精度,利用粒子群优化可以使转子的振动幅值进一步的下降,提高了转子***的平衡精度,使转子***的运行更加平稳。
Description
技术领域
本发明涉及轴承转子***在线动平衡技术领域,具体涉及一种基于粒子群算法的高速轴承转子***模态动平衡方法。
背景技术
高档数控机床的关键部件电主轴单元以及高端航空发动机等作为高端装备制造业的重要装备,其运行都依赖于旋转机械的高速高精密稳定运转,而轴承转子***是直接影响上述高端旋转机械运行的最核心组件。高速运行的轴承转子***由于运行故障造成质量偏心等因素会同时存在不平衡力以及不平衡力偶,从而会产生振动位移。当转子在其速度接近临界转速时转轴会产生挠度变形,甚至会共振引起机械破坏。由于转子横向振动的固有频率是多阶的,其相应的临界转速也是多阶的,会存在多阶振型不平衡量。因此,研究针对存在多阶不平衡的高速轴承转子***的动平衡方法具有重要的理论和应用价值。
近些年国内外学者也对高速转子多阶振型下的模态动平衡调控方法方面开展了诸多研究。李晓丰等在模态平衡理论的基础上,推导出一种高速柔性转子无试重模态动平衡方法,解决了影响系数法和模态平衡法需添加试重的问题,能有效、准确、快速地获取柔性转子的不平衡量大小及方位。宾光富等提出一种基于模态振型分析的大型柔性转子高速动平衡方法,采取逐阶次平衡的方法,有效的平衡了大型柔性转子,减少了起机次数。然而,目前对于如何进一步提高高速电主轴***的动平衡振动抑制水平是需要解决的重要问题。
发明内容
针对现有技术的不足,本发明提出一种基于粒子群算法的高速轴承转子***模态动平衡方法,包括以下步骤:
步骤1:根据轴的结构构建轴的三维模型,并在轴的两端添加轴承支撑,得到轴承转子***的三维模型;
步骤2:对轴承转子***的三维模型进行模态分析,将轴上变形量最大的N个截面作为平衡校正面;
步骤3:根据模态平衡法的平衡条件建立***的固有不平衡量与平衡校正面上需要添加的校正质量之间的函数关系;
步骤4:将需要添加的校正质量作为一个粒子,所有的校正质量作为一个粒子群,采用粒子群算法求解每个校正质量的最优值。
所述步骤3包括:
步骤3.1:利用公式(1)建立轴承转子***的振动微分方程,
式中,M表示轴承转子***的质量矩阵,C表示轴承转子***的阻尼矩阵,K表示轴承转子***的刚度矩阵,x(z)表示轴承转子***的振动量,F(z)表示轴承转子***受到的不平衡力;
步骤3.2:当F(z)=0时,且不考虑轴承转子阻尼的情况下,将公式(1)化简为:
步骤3.3:将F(z)表示为F(z)=q(z)esz并带入公式(1),其中s为实数,则x(z)的通解表示为x(z)=zkq(z)esz,k=0,1,2,将x(z)=zkq(z)esz带入公式(2)化简为:
步骤3.4:求解特征方程|Ms2+K|=0所有的特征值λi以及特征函数ξi(z);
步骤3.5:根据前N阶主振型对应的特征函数ξi(z)建立轴承转子***的固有不平衡量与校正质量之间的函数关系:
式中,ξi(z)表示第i阶主振型对应的特征函数,Qi表示第i个平衡校正面上需要添加的校正质量,Fi(z)表示F(z)在第i个平衡校正面上的分量;
将公式(4)表示为矩阵形式为:
所述步骤4包括:
步骤4.1:设计粒子群算法中的适应度函数为f(z):
式中,Qlim表示校正平面所允许添加的校正质量的最大值;
步骤4.2:设置最大迭代次数J并设置f(z)的阈值ψ,进行迭代计算,当达到最大迭代次数或第j次计算得到的fj(z)≤ψ时停止迭代计算,f(z)取最小值时对应的输出值作为每个校正质量的最优值,j=1,2,…,J。
本发明的有益效果是:
本发明提出了一种基于粒子群算法的高速轴承转子***模态动平衡方法,首先通过模态分析得到平衡校正面的截面位置,通过模态平衡法建立校正质量与***固有不平衡量之间的数学模型,最后通过粒子群算法求解出每个平衡校正面上的最优质量值,基于粒子群算法的高速轴承转子***模态动平衡方法可以有效提高转子的平衡精度,利用粒子群优化可以使转子的振动幅值进一步的下降,提高了转子***的平衡精度,使转子***的运行更加平稳。
附图说明
图1为本发明中基于粒子群算法的高速轴承转子***模态动平衡方法流程图;
图2为本发明中转子***的一阶模态振型;
图3为本发明中转子***的二阶模态振型;
图4为本发明中转子***的优化前后的振幅曲线图,其中图(a)为转子***平衡前的振幅曲线图,图(b)为未采用粒子群算法优化的模态平衡后的振幅曲线图,图(c)为采用粒子群算法优化后的振幅曲线图;
图5为本发明中转子***的坎贝尔图;
图6为本发明中粒子群算法的流程图。
具体实施方式
本发明提出一种基于粒子群优化的高速轴承转子***模态动平衡方法。建立轴承转子***的三维模型,进行模态分析得到转子***的前N阶模态振型和坎贝尔图,将轴上变形量最大的N个截面作为平衡校正面,测试并记录不平衡量。利用模态平衡法建立每个平衡校正面需要添加的校正质量的数学模型,建立以最低残余振动量为目标的优化模型,基于粒子群算法优化计算质量平衡补偿校正量,并进行优化平衡,得出轴承转子***在不平衡量下的校正质量,完成动平衡质量补偿。下面结合附图和具体实施实例对发明做进一步说明。
如图1所示,一种基于粒子群算法的高速轴承转子***模态动平衡方法,以实现对高速轴承转子***多阶不平衡振动有效抑制,包括如下步骤:
步骤1:根据轴的结构构建轴的三维模型,并在轴的两端添加轴承支撑,得到轴承转子***的三维模型;本实施例采用SolidWorks软件构建三维模型,在建模时设置的轴长为0.4m,直径为0.05m,弹性模量为2.1*1011,轴一端的轴承为7012滚子轴承,另一端的轴承为7009滚子轴承。
步骤2:采用ANSYS软件对轴承转子***的三维模型进行模态分析,将轴上变形量最大的N个截面作为平衡校正面;轴承转子***的一阶模态振型图如图2表示,二阶模态振型图如图3所示。
进行模态分析同时可以得到转子***的坎贝尔图,如图5所示,在坎贝尔图中可以的到转子的临界转速,图5中的交点A为一阶临界转速,交点B为二阶临界转速。转子的前二阶临界转速如表1所示,转子会在临界转速处发生强烈的振动,故进行实验时要尽量避开临界转速,选择了接近临界转速的6000r/min以及14000r/min作为平衡转速。故选择2000r/min、4000r/min、6000r/min、8000r/min、10000r/min、12000r/min、14000r/min、16000r/min、18000r/min、20000r/min共十种不同的转速下,对转子的不平衡进行研究。
表1转子的前二阶临界转速
从图2~3中可以看出在轴接近轴两端轴承的截面上变形量最大,故本实施例将两个端面作为不平衡量的平衡校正面。采用在转子的前端加不平衡量的方法来模拟***的固有不平衡量,分别添加上5g·mm、10g·mm和15g·mm三种不平衡量,相位分别是0°,90°和180°,后端不添加不平衡量,后端不平衡量为零。
步骤3:根据模态平衡法的平衡条件建立***的固有不平衡量与平衡校正面上需要添加的校正质量之间的函数关系,包括:
步骤3.1:利用公式(1)建立轴承转子***的振动微分方程,
式中,M表示轴承转子***的质量矩阵,C表示轴承转子***的阻尼矩阵,K表示轴承转子***的刚度矩阵,x(z)表示轴承转子***的振动量,F(z)表示轴承转子***受到的不平衡力;
步骤3.2:当F(z)=0时,且不考虑轴承转子阻尼的情况下,将公式(1)化简为:
步骤3.3:将F(z)表示为F(z)=q(z)esz并带入公式(1),其中s为实数,则x(z)的通解表示为x(z)=zkq(z)esz,k=0,1,2,将x(z)=zkq(z)esz带入公式(2)化简为:
N阶振型相互正交,有:
式中,Nf表示第f阶振型的模态质量,m为轴的质量,将转子***的不平衡量F(z)代入时:
转子***的振动曲线可以看成是各阶主振型按照af的比例叠加生成的,即:
令Af(Ω)为转子运行时的动态方法系数,可以看出,当转子的转速Ω接近于转子的某阶固频率λf时,Af(Ω)趋于无穷大,此时其他阶振型可以忽略。频率超过转子转速的高阶振型的影响很小,可以忽略,在平衡时只需要考虑前N阶振型。
当转子达到平衡时,则轴承的支撑反力为零,可得:
步骤3.5:根据前N阶主振型对应的特征函数ξi(z)建立轴承转子***的固有不平衡量与校正质量之间的函数关系:
式中,ξi(z)表示第i阶主振型对应的特征函数,Qi表示第i个平衡校正面上需要添加的校正质量,Fi(z)表示F(z)在第i个平衡校正面上的分量;
将公式(4)表示为矩阵形式为:
步骤4:将需要添加的校正质量作为一个粒子,所有的校正质量作为一个粒子群,采用粒子群算法求解每个校正质量的最优值,包括:
步骤4.1:设计粒子群算法中的适应度函数为f(z):
式中,Qlim表示校正平面所允许添加的校正质量的最大值;
步骤4.2:设置最大迭代次数J并设置f(z)的阈值ψ,进行迭代计算,当达到最大迭代次数或第j次计算得到的fj(z)≤ψ时停止迭代计算,f(z)取最小值时对应的输出值作为每个校正质量的最优值,j=1,2,…,J。
利用粒子群算法求解最优值的过程如图6所示,其运算过程如下:
(1)随机生成N组配重质量作为初始解,即为粒子群算法的初始种群;
(2)随机生成2个随机数作为配重块质量的飞行速度以及距离;
(3)计算每一组配重质量所对应的适应度函数;
(4)比较计算出的适应度函数值,记录配重块自身的最优质量和全局最优质量。
(5)判断是否达到得到所允许范围内的全局最优配重质量,若是则输出结果,若不是则继续进行下一步;
(6)根据设置的初始飞行速度以及距离更新配重块质量,并改变配重质量的下一次的飞行速度以及距离;返回步骤(3),进行下一轮迭代,直到得到最优配重质量。
通过粒子群算法计算出来的主轴在5g·mm、10g·mm、15g·mm不平衡量的优化校正量的大小,如表1~3所示,
表1 5g·mm不平衡质量的试验数据
表2 10g·mm不平衡质量的试验数据
表3 15g·mm不平衡质量的试验数据
在上面的三个表格中,可以看到在不平衡量5g·mm、10g·mm和15g·mm下,转子***平衡时需要的配重质量。可以得到如图4所示的转子***的振幅。从图4中可以看出,转子***在6000r/min以及14000r/min时的振幅波动较大。利用模态平衡法平衡转子***后的振幅如图4中图(b)所示。平衡效率最大出现在转速为2000r/min,不平衡量为5g·mm时,最大的平衡率为71.31%,平均的平衡率为66.0%。经过优化平衡后转子的振幅如图4中图(c)所示,平衡效率最大出现在转速为6000r/min,不平衡量为10g·mm时,优化后转子的振幅最大的平衡率为75.91%,平均的平衡率为73.13%,相较于模态平衡后振幅的最大下降为34.75%,平均下降20.86%。
Claims (1)
1.一种基于粒子群算法的高速轴承转子***模态动平衡方法,其特征在于,包括如下步骤:
步骤1:根据轴的结构构建轴的三维模型,并在轴的两端添加轴承支撑,得到轴承转子***的三维模型;
步骤2:对轴承转子***的三维模型进行模态分析,将轴上变形量最大的N个截面作为平衡校正面;
步骤3:根据模态平衡法的平衡条件建立***的固有不平衡量与平衡校正面上需要添加的校正质量之间的函数关系;
步骤4:将需要添加的校正质量作为一个粒子,所有的校正质量作为一个粒子群,采用粒子群算法求解每个校正质量的最优值;
所述步骤3包括:
步骤3.1:利用公式(1)建立轴承转子***的振动微分方程,
式中,M表示轴承转子***的质量矩阵,C表示轴承转子***的阻尼矩阵,K表示轴承转子***的刚度矩阵,x(z)表示轴承转子***的振动量,F(z)表示轴承转子***受到的不平衡力;
步骤3.2:当F(z)=0时,且不考虑轴承转子阻尼的情况下,将公式(1)化简为:
步骤3.3:将F(z)表示为F(z)=q(z)esz并带入公式(1),其中s为实数,则x(z)的通解表示为x(z)=zkq(z)esz,k=0,1,2,将x(z)=zkq(z)esz带入公式(2)化简为:
步骤3.4:求解特征方程|Ms2+K|=0所有的特征值λi以及特征函数ξi(z);
步骤3.5:根据前N阶主振型对应的特征函数ξi(z)建立轴承转子***的固有不平衡量与校正质量之间的函数关系:
式中,ξi(z)表示第i阶主振型对应的特征函数,Qi表示第i个平衡校正面上需要添加的校正质量,Fi(z)表示F(z)在第i个平衡校正面上的分量;
将公式(4)表示为矩阵形式为:
所述步骤4包括:
步骤4.1:设计粒子群算法中的适应度函数为f(z):
式中,Qlim表示校正平面所允许添加的校正质量的最大值;
步骤4.2:设置最大迭代次数J并设置f(z)的阈值ψ,进行迭代计算,当达到最大迭代次数或第j次计算得到的fj(z)≤ψ时停止迭代计算,f(z)取最小值时对应的输出值作为每个校正质量的最优值,j=1,2,…,J。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011319498.XA CN112556931B (zh) | 2020-11-23 | 2020-11-23 | 基于粒子群算法的高速轴承转子***模态动平衡方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011319498.XA CN112556931B (zh) | 2020-11-23 | 2020-11-23 | 基于粒子群算法的高速轴承转子***模态动平衡方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112556931A CN112556931A (zh) | 2021-03-26 |
CN112556931B true CN112556931B (zh) | 2022-09-06 |
Family
ID=75044767
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011319498.XA Active CN112556931B (zh) | 2020-11-23 | 2020-11-23 | 基于粒子群算法的高速轴承转子***模态动平衡方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112556931B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114838872B (zh) * | 2022-05-10 | 2023-04-04 | 浙江大学 | 一种超重力离心机动平衡长期监测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101949753A (zh) * | 2010-09-16 | 2011-01-19 | 西北工业大学 | 一种高速柔性转子动平衡方法 |
JP2013003021A (ja) * | 2011-06-17 | 2013-01-07 | Nagahama Seisakusho Ltd | 不釣合い測定および修正装置ならびに不釣合い測定および修正方法 |
GB201310478D0 (en) * | 2013-06-12 | 2013-07-24 | Universal Balancing Ltd | Method and apparatus for balancing a rotor |
CN103257018A (zh) * | 2013-05-14 | 2013-08-21 | 上海工程技术大学 | 刚性转子动平衡的获取方法 |
CN107255549A (zh) * | 2017-06-14 | 2017-10-17 | 西安交通大学 | 一种转子多点不平衡量识别方法 |
-
2020
- 2020-11-23 CN CN202011319498.XA patent/CN112556931B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101949753A (zh) * | 2010-09-16 | 2011-01-19 | 西北工业大学 | 一种高速柔性转子动平衡方法 |
JP2013003021A (ja) * | 2011-06-17 | 2013-01-07 | Nagahama Seisakusho Ltd | 不釣合い測定および修正装置ならびに不釣合い測定および修正方法 |
CN103257018A (zh) * | 2013-05-14 | 2013-08-21 | 上海工程技术大学 | 刚性转子动平衡的获取方法 |
GB201310478D0 (en) * | 2013-06-12 | 2013-07-24 | Universal Balancing Ltd | Method and apparatus for balancing a rotor |
CN107255549A (zh) * | 2017-06-14 | 2017-10-17 | 西安交通大学 | 一种转子多点不平衡量识别方法 |
Non-Patent Citations (2)
Title |
---|
基于改进粒子群算法的最小二乘影响系数法的理论及实验研究;王星星 等;《振动与冲击》;20130428;第32卷(第8期);第1-5节 * |
高速主轴动平衡方法研究现状;王展 等;《现代制造工程》;20190618(第6期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112556931A (zh) | 2021-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103821567B (zh) | 一种航空发动机高压转子结构动力学设计方法 | |
CN111475903A (zh) | 基于多偏置误差同步补偿的大型高速回转装备多级零部件动力学特性逐级测调和分配方法 | |
CN110457740B (zh) | 一种机械结构在基础激励下的响应特性分析方法 | |
Xul et al. | Vibration characteristics of unbalance response for motorized spindle system | |
CN112525180B (zh) | 一种基于混合重复控制器的谐波振动力矩抑制方法 | |
CN109632188B (zh) | 虑及多转速状态的柔性转子动平衡双重优化方法及*** | |
Phadatare et al. | Nonlinear frequencies and unbalanced response analysis of high speed rotor-bearing systems | |
CN112556931B (zh) | 基于粒子群算法的高速轴承转子***模态动平衡方法 | |
CN110145541A (zh) | 一种基于相位稳定的磁悬浮轴承转子不平衡运动控制方法 | |
CN111141518A (zh) | 一种基于模型的非对称转子轴承***不平衡量识别方法 | |
CN114330083A (zh) | 一种轴承参数不确定的离心泵转子动力学稳健预测方法 | |
Zhu et al. | Dynamic analysis of a flexible rotor supported by ball bearings with damping rings based on FEM and lumped mass theory | |
Messager et al. | Discrete optimization of rigid rotor balancing | |
CN108534996B (zh) | 一种基于微振动瀑布图的飞轮刚度参数提取方法 | |
Fegade et al. | Unbalanced response of rotor using ansys parametric design for different bearings | |
Kang et al. | A modified influence coefficient method for balancing unsymmetrical rotor-bearing systems | |
CN110646139B (zh) | 根据弯曲度确定轴弹性的转子的不平衡度的方法 | |
Hundal et al. | Balancing of flexible rotors having arbitrary mass and stiffness distribution | |
CN114577397B (zh) | 一种高速永磁电机转子动平衡方法及*** | |
Genta et al. | DYNROT: a finite element code for rotordynamic analysis based on complex co‐ordinates | |
Hong et al. | Vibration failure analysis of multi-disk high-speed rotor based on rotary inertia load model | |
CN113268908B (zh) | 一种转子***的响应求解方法和装置 | |
CN1103443C (zh) | 非对称转子的全息动平衡方法 | |
Liu et al. | An efficient dynamic modeling technique for a central tie rod rotor | |
Nelson | Steady synchronous response and balancing of rotor systems with residual shaft bow |
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 |