CN105956254A - 一种用于转子系统动力学建模的刚体单元法 - Google Patents
一种用于转子系统动力学建模的刚体单元法 Download PDFInfo
- Publication number
- CN105956254A CN105956254A CN201610272794.6A CN201610272794A CN105956254A CN 105956254 A CN105956254 A CN 105956254A CN 201610272794 A CN201610272794 A CN 201610272794A CN 105956254 A CN105956254 A CN 105956254A
- Authority
- CN
- China
- Prior art keywords
- rigid unit
- jth
- rotor
- beta
- alpha
- 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
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]
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)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开一种用于转子系统动力学建模的刚体单元法,包括:1)收集转子的几何参数、材料属性和运行状态;2)将转子划分为若干个离散的刚体单元,任意相邻的两个刚体单元之间通过四个假想的弹簧相连接;3)根据转子相邻刚体单元之间相互作用模型计算作用在每个刚体单元上的合力和合力矩,并建立每个刚体单元的动力学微分方程;4)使用变步长的四阶龙格‑库塔‑费尔伯格方法求解各个刚体单元的动力学方程来获得刚体单元的振动响应。利用本发明建立的转子系统动力学模型,不需要任何假设,就能够与Gupta轴承动力学模型实现耦合,建立整个滚动轴承‑转子系统的完全动力学模型。
Description
技术领域
本发明涉及一种用于转子系统动力学建模的方法。
背景技术
在转子动力学研究中,合理地对转子系统进行动力学建模是进行后续分析和计算的基础。从现有文献检索发现,转子动力学建模方法主要有理论法、传递矩阵法和有限元法,其中基于理论力学的建模方法一般只能用于简单的转子系统建模,对于复杂的转子系统,较常用的是传递矩阵法和有限元法。传递矩阵法的主要特点是矩阵的阶数不随系统的自由度数增大而增加,因此编程简单,所需内存小,运算速度快,但是该方法在考虑转子周围结构时分析计算比较困难,而且传递矩阵法对转子系统进行了简化和等效,无法保证模型的完整性和分析的准确度。
有限元法基于转子的几何结构,将模型离散为一系列单元进行分析,其表达式简单规范,能够同时保证模型的完整性和计算准确度,也很容易借助有限元软件实现,很好地兼顾到模型的完整性和计算效率。另外,利用有限元分析软件可以很好地解决转子动力特性分析中的“陀螺效应”问题,可以提高准确度和分析效率。申请号为201510816868.3的中国发明专利申请公开了一种基于双转子简化动力学模型设计的航空发动机仿真试验台,其特点在于首先采用有限元法建立多轮盘的复杂离散动力学模型,然后基于质心集中方法对该模型进行了简化。可以看出,利用有限元方法建立的模型结构复杂,适合用于线性振动特性的定量计算,但对于非线性问题的定量分析计算效率较低,定性分析难以实现;且存在自由度数特别大,耗费计算机时等问题。
此外,在研究轴承-转子系统的动力学行为时,须建立转子和滚动轴承的耦合模型。为了考虑转子的弹性振动和柔性,目前较为理想的方法是采用转子有限元模型与Gupta滚动轴承动力学模型实现耦合,但两者进行耦合时,为了求解有限元模型,需要采用一定的假设,而这些假设的合理性还有待于进一步验证,因此采用转子有限元模型与Gupta轴承模型耦合的方法具有一定的局限性。
发明内容
本发明的目的在于提供一种用于转子系统动力学建模的刚体单元法:通过将转子分隔成由弹簧相连接的离散的刚体单元,根据转子相邻刚体单元之间相互作用模型计算作用在每个刚体单元上的合力和合力矩,并建立每个刚体单元的动力学微分方程,然后使用变步长的四阶龙格-库塔-费尔伯格方法求解各个刚体单元的动力学方程来获得刚体单元的振动响应。本发明提出了一种全新的转子系统动力学建模方法,与传统方法相比较,其显著优点是:考虑了转子的弹性振动和柔性,操作方法稳定可靠、效率较高;并且利用该方法建立的转子系统动力学模型,不需要任何假设,就能够与Gupta轴承动力学模型实现耦合,建立整个滚动轴承-转子系统的完全动力学模型。
为了实现上述目的,本发明采用如下技术方案:
一种用于转子系统动力学建模的刚体单元法,包括以下步骤:
1)收集转子的几何参数、材料属性和运行状态;
2)将转子划分为若干个离散的刚体单元,任意相邻的两个刚体单元之间通过四个假想的弹簧相连接,所述四个假想的弹簧包括一个约束平移运动的拉伸弹簧和三个约束旋转运动的扭转弹簧;
3)根据转子各个相邻刚体单元之间相互作用模型,计算作用在每个转子刚体单元上的合力和合力矩,并建立每个刚体单元的动力学微分方程;
4)使用变步长的四阶龙格-库塔-费尔伯格方法求解各个刚体单元的动力学微分方程来获得每个刚体单元的振动响应。
所述步骤2)中,在将转子划分为若干个离散的刚体单元时,为了使建立的转子动力学模型足够合理,并方便计算作用在刚体单元上的合力和合力矩,需要遵循以下标准:
(1)划分的每一个刚体单元的内径或外径必须一致,同一个刚体单元不能具有不同的内径或外径;
(2)为了保证各个节点振幅相对误差在合理范围内,建议划分后的每一个刚体单元的长径比不超过5;
(3)每个轴承在转子上的支承位置的中心必须与它所对应的刚体单元的节点对齐,而不能与两个相邻刚体单元之间的面对齐。
进一步的,所述步骤3)中,定义惯性坐标系i(xi,yi,zi)和第j个刚体单元的定体坐标系rj(xrj,yrj,zrj),可通过刚体单元的姿态角确定。每个刚体单元有两个端面,左端被命名为α平面,右端被命名为β平面。Mj和Nj分别是第j个刚体单元α平面和β平面的中心。同样,Mj+1和Nj+1分别是第j+1个刚体单元α平面和β平面的中心。第j个刚体单元和第j+1个刚体单元之间的几何作用关系为:当转子在自由状态时,第j个刚体单元的β平面与第j+1个刚体单元的α平面重合,此时它们之间的弹簧处于松弛状态;当转子受到外力时,由于相邻刚体单元运动的不一致,连接相邻刚体单元平面之间的弹簧将发生压缩或者拉伸,从而产生了相互作用力。以第j个刚体单元的β平面和第j+1个刚体单元的α平面的相互作用为例,使用上标i表示惯性坐标系。在惯性坐标系中,Nj和Mj+1分别由向量和确定。由拉伸弹簧产生的作用在第j个刚体单元上的相互作用力和作用在第j+1个刚体单元上的相互作用力表示为:
式中:——第j个刚体单元β面与第j+1个刚体单元α面之间的拉伸弹簧的刚度/N.m-1;E——转子的弹性模量/Pa;lj和lj+1——第j个和第j+1个刚体单元的长度/m;Aj和Aj+1——第j个和第j+1个刚体单元端面的面积/m2;函数min(a,b)计算a和b的最小值。
除了拉伸弹簧的相互作用力,同时还有由三个方向上的扭转弹簧产生的作用在第j个刚体单元的作用力矩和作用在第j+1个刚体单元的作用力矩相邻刚体单元的姿态角为(ηj,ξj,λj)和(ηj+1,ξj+1,λj+1),相互作用力矩和由下式确定:
式中:——第j个刚体单元β面与第j+1个刚体单元α面之间x方向扭转弹簧的刚度/N.m.rad-1,G——转子的剪切模量/Pa;Ipj和Ipj+1——第j个和第j+1个刚体单元的极惯性矩/m4;
力矩分量和通过下式确定:
式中:——第j个刚体单元β面与第j+1个刚体单元α面之间y方向扭转弹簧的刚度/N.m.rad-1;Iyj和Iyj+1——第j个刚体单元和第j+1个刚体单元相对于y轴的惯性矩/m4;
力矩分量和表示为:
式中:——第j个刚体单元β面与第j+1个刚体单元α面之间z方向扭转弹簧的刚度/N.m.rad-1;Izj和Izj+1——第j个刚体单元和第j+1个刚体单元相对于z轴的惯性矩/m4。
进一步的,步骤3)中,将每个刚体单元的运动分解成为质心的平动运动和绕其质心的转动;和分别是由不平衡质量产生的不平衡力和施加在第j个刚体单元上的外力,则第j个刚体单元的平动运动方程为:
式中:mj——第j个刚体单元的质量/kg;(xj,yj,zj)——质心Orj在惯性坐标系中的坐标分量;Gj——刚体单元的重力/N;不平衡力表示为:
式中:mujruj——不平衡质量和其半径的乘积/kg.m;αj——不平衡质量的初始方位角/rad;Trj,i——从刚体单元定体坐标系到惯性坐标系的变换矩阵;
第j个刚体单元的角速度为ωj(ωjx,ωjy,ωjz),则刚体单元的旋转运动方程表示为:
式中施加在第j个刚体单元上的外部力矩和是由和产生的力矩,表示为:
式中:Ti,rj——从惯性坐标系到刚体单元定体坐标系的变换矩阵。
所述步骤4)中,使用变步长四阶龙格-库塔-费尔伯格方法进行积分,可得到下一时刻各刚体单元的速度和位移。重复此步骤可求得每个刚体单元的振动响应,对应于转子不同位置的振动响应。
本发明的进一步改进在于,转子模型的求解步骤为:首先根据转子的几何参数、材料属性、运行状态计算各转子单元初始速度和位移,然后根据转子相邻刚体单元之间相互作用模型计算作用在每个刚体单元上的合力和合力矩,并建立每个刚体单元的动力学微分方程;使用变步长的四阶龙格-库塔-费尔伯格方法求解各个刚体单元的动力学方程来获得刚体单元的振动响应。
与现有技术相比,本发明具有下述优点:
1、提出了一个新的转子系统动力学建模方法(DM),将转子分隔成离散的刚体单元,通过求解各个刚体单元的动力学方程来获得刚体单元的动力学响应;
2、利用本发明的方法建立的转子动力学模型,不需要做任何假设,就能够与Gupta轴承动力学模型实现耦合,建立整个轴承-转子系统的完全动力学模型。
本发明通过将转子分隔成由弹簧相连接的离散的刚体单元,根据转子相邻刚体单元之间相互作用模型计算作用在每个刚体单元上的合力和合力矩,并建立每个刚体单元的动力学微分方程,使用变步长的四阶龙格-库塔-费尔伯格方法求解各个刚体单元的动力学方程来获得刚体单元的振动响应。利用该方法建立的转子系统动力学模型,不需要任何假设,就能够与Gupta轴承动力学模型实现耦合,建立整个滚动轴承-转子系统的完全动力学模型。
附图说明
图1是本发明中转子的任意两个相连接刚体单元之间的相互作用:坐标系分别为惯性坐标系i(xi,yi,zi)和第j个刚体单元的定体坐标系rj(xrj,yrj,zrj);
图2是本发明的流程示意图;
图3是本发明的转子模型,其中图3(a)为转子系统动力学模型,图3(b)为转子系统有限元模型;
图4是本发明中节点4在y方向的时域响应图;
图5是本发明中节点4在z方向的时域响应图;
图6是本发明中转子各节点的位移图;
图7是本发明中转子转速为10000r/min时不同节点的轨迹。
具体实施方式
请参阅图1及图2所示,本发明一种用于转子系统动力学建模的刚体单元法,包括以下步骤:
1)收集转子的几何参数、材料属性和运行状态,为转子系统动力学建模提供数据支持;
2)将转子划分为若干个离散的刚体单元,任意相邻的两个刚体单元之间通过四个假想的弹簧相连接,即一个约束平移运动的拉伸弹簧和三个约束旋转运动的扭转弹簧;
3)根据转子各个相邻刚体单元之间相互作用模型,计算作用在每个转子刚体单元上的合力和合力矩,并建立每个刚体单元的动力学微分方程;
4)使用变步长的四阶龙格-库塔-费尔伯格方法求解各个刚体单元的动力学方程来获得每个刚体单元的振动响应。
步骤1)中,转子的材料特性参数包括转子长度、直径,以及弹性模量、泊松比和密度。
步骤2)中,在将转子划分为若干个离散的刚体单元时,为了使建立的转子动力学模型足够合理,并方便计算作用在刚体单元上的合力和合力矩,需要遵循以下标准:
(1)划分的每一个刚体单元的内径或外径必须一致,同一个刚体单元不能具有不同的内径或外径;
(2)为了保证各个节点振幅相对误差在合理范围内,建议划分后的每一个刚体单元的长径比不超过5;
(3)每个轴承在转子上的支承位置的中心必须与它所对应的刚体单元的节点对齐,而不能与两个相邻刚体单元之间的面对齐。
步骤4)中,使用变步长四阶龙格-库塔-费尔伯格方法进行积分,可得到下一时刻各刚体单元的速度和位移。重复此步骤可求得每个刚体单元的振动响应。
下面结合一个实例对本发明作进一步详细说明,但本实例并不用于限制本发明。其步骤为:
1)收集转子的几何参数、材料属性和运行状态,为转子系统动力学建模提供数据支持,转子的参数如表1所示。
表1转子参数
2)将转子等分成9个刚体单元(其中第2个和第8个刚体单元通过弹簧和阻尼支承),任意相邻的两个刚体单元之间通过四个假想的弹簧相连接,即一个约束平移运动的拉伸弹簧和三个约束旋转运动的扭转弹簧。
3)根据相邻刚体单元之间的相互作用模型,对转子各刚体单元之间的相互作用进行分析,计算作用在每个转子单元上的力和力矩,建立每个刚体单元的运动微分方程,转子动力学模型(DM)如图3(a)所示。具体计算方法如下:
选取任意两个相邻刚体单元(第j个和第j+1个),定义惯性坐标系i(xi,yi,zi)和第j个刚体单元的定体坐标系rj(xrj,yrj,zrj),可通过刚体单元的姿态角确定。每个刚体单元的左端面被命名为α平面,右端面被命名为β平面。Mj和Nj分别是第j个刚体单元α平面和β平面的中心。同样,Mj+1和Nj+1分别是第j+1个刚体单元α平面和β平面的中心。使用上标i表示惯性坐标系,在惯性坐标系中,Nj和Mj+1分别由向量和确定。当转子受到外力时,由拉伸弹簧产生的作用在第j个刚体单元上的相互作用力和作用在第j+1个刚体单元上的相互作用力可根据下式来计算
式中:——第j个刚体单元β面与第j+1个刚体单元α面之间的拉伸弹簧的刚度/N.m-1;E——转子的弹性模量/Pa;lj和lj+1——第j个和第j+1个刚体单元的长度/m;Aj和Aj+1——第j个和第j+1个刚体单元端面的面积/m2。函数min(a,b)计算a和b的最小值。除了拉伸弹簧的相互作用力,同时还有由三个方向上的扭转弹簧产生的作用在第j个刚体单元的作用力矩和作用在第j+1个刚体单元的作用力矩相邻刚体单元的姿态角为(ηj,ξj,λj)和(ηj+1,ξj+1,λj+1),则相互作用力矩和可由下式确定
式中:——第j个刚体单元β面与第j+1个刚体单元α面之间x方向扭转弹簧的刚度/N.m.rad-1,G——转子的剪切模量/Pa;Ipj和Ipj+1——第j个和第j+1个刚体单元的极惯性矩/m4。
力矩分量和可以通过下式确定
式中:——第j个刚体单元β面与第j+1个刚体单元α面之间y方向扭转弹簧的刚度/N.m.rad-1;Iyj和Iyj+1——第j个刚体单元和第j+1个刚体单元相对于y轴的惯性矩/m4。
力矩分量和可表示为
式中:——第j个刚体单元β面与第j+1个刚体单元α面之间z方向扭转弹簧的刚度/N.m.rad-1;Izj和Izj+1——第j个刚体单元和第j+1个刚体单元相对于z轴的惯性矩/m4。
每个刚体单元的运动可以分解成为质心的平动运动和绕其质心的转动。假设和分别是由不平衡质量产生的不平衡力和施加在第j个刚体单元上的外力,则第j个刚体单元的平动运动方程为
式中:mj——第j个刚体单元的质量/kg;(xj,yj,zj)——质心Orj在惯性坐标系中的坐标分量;Gj——刚体单元的重力/N。不平衡力可以表示为
式中:mujruj——不平衡质量和其半径的乘积/kg.m;αj——不平衡质量的初始方位角/rad;Trj,i——从刚体单元定体坐标系到惯性坐标系的变换矩阵。
假设第j个刚体单元的角速度为ωj(ωjx,ωjy,ωjz),则刚体单元的旋转运动方程可以表示为
式中施加在第j个刚体单元上的外部力矩和是由和产生的力矩,可以表示为
式中:Ti,rj——从惯性坐标系到刚体单元定体坐标系的变换矩阵。
4)使用变步长四阶龙格-库塔-费尔伯格方法求解各个刚体单元的动力学方程,最大允许截断误差被设定为1.0E-5,可以得到任意一个或所关注刚体单元的振动响应。
本实例取转子动力学模型的节点4为研究对象,当转子转速为1000r/min时,在节点4的z方向施加一个大小为100N的正弦激振力,通过仿真分析,可以得到它的y方向和z方向的时域响应,分别如图4和图5所示,可以看出节点4在y方向和z方向的时域响应曲线均为正弦曲线。重复以上方法,就可以求解任意一个或所关注刚体单元的振动响应。
另外,图3(b)给出了包含10个单元的转子的有限元模型(FEM),第1个和第9个单元的长度均为20mm,其他单元的长度都为40mm。转子的参数如表1示所列。有限元模型使用Newmarkβ方法求解,计算时间步长设定为2.0E-4s。
1)假设当转子转速为0时,在节点4的z方向施加一个径向力,图6为不同大小的径向力作用下转子各节点在z方向的位移。从图6可以看出当转子受到静载荷时,FEM的计算结果和DM的计算结果匹配较好。
2)假设在节点4处存在不平衡质量,不平衡质量和其半径的乘积是0.001kg.m。当转子旋转时,转子在不平衡力作用下会产生振动响应。图7是当转速为10000r/min时各个节点的轨迹(图中:FEM——有限元模型计算的结果,DM——动力学模型计算的结果),由于支承弹簧的限制,节点2和节点8的振幅要小于其他节点。
Claims (5)
1.一种用于转子系统动力学建模的刚体单元法,其特征在于,包括以下步骤:
1)收集转子的几何参数、材料属性和运行状态;
2)将转子划分为若干个离散的刚体单元,任意相邻的两个刚体单元之间通过四个假想的弹簧相连接,所述四个假想的弹簧包括一个约束平移运动的拉伸弹簧和三个约束旋转运动的扭转弹簧;
3)根据转子各个相邻刚体单元之间相互作用模型,计算作用在每个转子刚体单元上的合力和合力矩,并建立每个刚体单元的动力学微分方程;
4)使用变步长的四阶龙格-库塔-费尔伯格方法求解各个刚体单元的动力学微分方程来获得每个刚体单元的振动响应。
2.根据权利要求1所述的一种用于转子系统动力学建模的刚体单元法,其特征在于,所述步骤3)中,定义惯性坐标系i(xi,yi,zi)和第j个刚体单元的定体坐标系rj(xrj,yrj,zrj),每个刚体单元有两个端面,左端被命名为α平面,右端被命名为β平面,Mj和Nj分别是第j个刚体单元α平面和β平面的中心,同样,Mj+1和Nj+1分别是第j+1个刚体单元α平面和β平面的中心,第j个刚体单元和第j+1个刚体单元之间的几何作用关系为:当转子在自由状态时,第j个刚体单元的β平面与第j+1个刚体单元的α平面重合,此时它们之间的弹簧处于松弛状态;当转子受到外力时,由于相邻刚体单元运动的不一致,连接相邻刚体单元平面之间的弹簧将发生压缩或者拉伸,从而产生了相互作用力;使用上标i表示惯性坐标系,在惯性坐标系中,Nj和Mj+1分别由向量和确定,由拉伸弹簧产生的作用在第j个刚体单元上的相互作用力和作用在第j+1个刚体单元上的相互作用力表示为:
式中:——第j个刚体单元β面与第j+1个刚体单元α面之间的拉伸弹簧的刚度/N.m-1;E——转子的弹性模量/Pa;lj和lj+1——第j个和第j+1个刚体单元的长度/m;Aj和Aj+1——第j个和第j+1个刚体单元端面的面积/m2;函数min(a,b)计算a和b的最小值;
除了拉伸弹簧的相互作用力,同时还有由三个方向上的扭转弹簧产生的作用在第j个刚体单元的作用力矩和作用在第j+1个刚体单元的作用力矩相邻刚体单元的姿态角为(ηj,ξj,λj)和(ηj+1,ξj+1,λj+1),相互作用力矩和由下式确定:
式中:——第j个刚体单元β面与第j+1个刚体单元α面之间x方向扭转弹簧的刚度/N.m.rad-1,G——转子的剪切模量/Pa;Ipj和Ipj+1——第j个和第j+1个刚体单元的极惯性矩/m4;
力矩分量和通过下式确定:
式中:——第j个刚体单元β面与第j+1个刚体单元α面之间y方向扭转弹簧的刚度/N.m.rad-1;Iyj和Iyj+1——第j个刚体单元和第j+1个刚体单元相对于y轴的惯性矩/m4;
力矩分量和表示为:
式中:——第j个刚体单元β面与第j+1个刚体单元α面之间z方向扭转弹簧的刚度/N.m.rad-1;Izj和Izj+1——第j个刚体单元和第j+1个刚体单元相对于z轴的惯性矩/m4。
3.根据权利要求1所述的一种用于转子系统动力学建模的刚体单元法,其特征在于,步骤3)中,将每个刚体单元的运动分解成为质心的平动运动和绕其质心的转动;和分别是由不平衡质量产生的不平衡力和施加在第j个刚体单元上的外力,则第j个刚体单元的平动运动方程为:
式中:mj——第j个刚体单元的质量/kg;(xj,yj,zj)——质心Orj在惯性坐标系中的坐标分量;Gj——刚体单元的重力/N;不平衡力表示为:
式中:mujruj——不平衡质量和其半径的乘积/kg.m;αj——不平衡质量的初始方位角/rad;Trj,i——从刚体单元定体坐标系到惯性坐标系的变换矩阵;
第j个刚体单元的角速度为ωj(ωjx,ωjy,ωjz),则刚体单元的旋转运动方程表示为:
式中施加在第j个刚体单元上的外部力矩和是由和产生的力矩,表示为:
式中:Ti,rj——从惯性坐标系到刚体单元定体坐标系的变换矩阵。
4.根据权利要求1所述的一种用于转子系统动力学建模的刚体单元法,其特征在于,步骤4)中,使用变步长四阶龙格-库塔-费尔伯格方法进行积分,得到下一时刻各刚体单元的速度和位移;重复此步骤求得每个刚体单元的振动响应,对应于转子不同位置的振动响应。
5.根据权利要求1所述的一种用于转子系统动力学建模的刚体单元法,其特征在于,步骤2)中,采用如下标准将转子划分为若干个离散的刚体单元:
(1)同一个刚体单元不能具有不同的内径,同时也不能具有不同的外径;
(2)划分后的每一个刚体单元的长径比不超过5;
(3)每个轴承在转子上的支承位置的中心必须与它所对应的刚体单元的节点对齐,而不能与两个相邻刚体单元之间的面对齐。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610272794.6A CN105956254A (zh) | 2016-04-27 | 2016-04-27 | 一种用于转子系统动力学建模的刚体单元法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610272794.6A CN105956254A (zh) | 2016-04-27 | 2016-04-27 | 一种用于转子系统动力学建模的刚体单元法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105956254A true CN105956254A (zh) | 2016-09-21 |
Family
ID=56916977
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610272794.6A Pending CN105956254A (zh) | 2016-04-27 | 2016-04-27 | 一种用于转子系统动力学建模的刚体单元法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105956254A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106528959A (zh) * | 2016-10-19 | 2017-03-22 | 天津大学 | 一种旋转对称结构固有频率和稳定性的简化分析方法 |
CN106599422A (zh) * | 2016-12-02 | 2017-04-26 | 长沙山水节能研究院有限公司 | 叶片泵转子系统的振动仿真分析方法及其装置 |
CN108959829A (zh) * | 2018-08-16 | 2018-12-07 | 居鹤华 | 基于轴不变量的非理想关节机器人动力学建模与解算方法 |
CN109194050A (zh) * | 2018-09-05 | 2019-01-11 | 西安交通大学 | 一种基于中心复合设计的转子结构优化方法 |
CN109813540A (zh) * | 2019-01-30 | 2019-05-28 | 西北工业大学 | 避免航空发动机盘腔积液失稳转子的设计及监测运行方法 |
CN111266032A (zh) * | 2020-02-12 | 2020-06-12 | 淮南北新建材有限公司 | 一种混合机转子及其构建方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7798720B1 (en) * | 2006-11-16 | 2010-09-21 | Florida Turbine Technologies, Inc. | Squeeze film damper with highly variable support stiffness |
CN102915388A (zh) * | 2012-09-19 | 2013-02-06 | 上海交通大学 | 一种面向对象的转子动力学系统非线性非因果建模仿真方法 |
CN103712785A (zh) * | 2013-09-12 | 2014-04-09 | 西安交通大学 | 一种滚动轴承‐转子系统损伤振动响应定量计算方法 |
CN104331565A (zh) * | 2014-11-10 | 2015-02-04 | 河海大学常州校区 | 轴类磁悬浮刚性转子系统的动力学建模方法及控制方法 |
-
2016
- 2016-04-27 CN CN201610272794.6A patent/CN105956254A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7798720B1 (en) * | 2006-11-16 | 2010-09-21 | Florida Turbine Technologies, Inc. | Squeeze film damper with highly variable support stiffness |
CN102915388A (zh) * | 2012-09-19 | 2013-02-06 | 上海交通大学 | 一种面向对象的转子动力学系统非线性非因果建模仿真方法 |
CN103712785A (zh) * | 2013-09-12 | 2014-04-09 | 西安交通大学 | 一种滚动轴承‐转子系统损伤振动响应定量计算方法 |
CN104331565A (zh) * | 2014-11-10 | 2015-02-04 | 河海大学常州校区 | 轴类磁悬浮刚性转子系统的动力学建模方法及控制方法 |
Non-Patent Citations (3)
Title |
---|
HONGRUI CAO等: "A New Dynamic Model of Ball-Bearing Rotor Systems based on Rigid Body Element", 《JOURNAL OF MANUFACTURING SCIENCE AND ENGINEER》 * |
LI YAMIN: "A General Method for the Dynamic Modeling of Ball Bearing–Rotor Systems", 《JOURNAL OF MANUFACTURING SCIENCE AND ENGINEER》 * |
曹宏瑞等: "高速滚动轴承-转子系统时变轴承刚度及振动响应分析", 《机械工程学报》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106528959A (zh) * | 2016-10-19 | 2017-03-22 | 天津大学 | 一种旋转对称结构固有频率和稳定性的简化分析方法 |
CN106599422A (zh) * | 2016-12-02 | 2017-04-26 | 长沙山水节能研究院有限公司 | 叶片泵转子系统的振动仿真分析方法及其装置 |
CN106599422B (zh) * | 2016-12-02 | 2020-05-26 | 长沙山水节能研究院有限公司 | 叶片泵转子系统的振动仿真分析方法及其装置 |
CN108959829A (zh) * | 2018-08-16 | 2018-12-07 | 居鹤华 | 基于轴不变量的非理想关节机器人动力学建模与解算方法 |
CN109194050A (zh) * | 2018-09-05 | 2019-01-11 | 西安交通大学 | 一种基于中心复合设计的转子结构优化方法 |
CN109813540A (zh) * | 2019-01-30 | 2019-05-28 | 西北工业大学 | 避免航空发动机盘腔积液失稳转子的设计及监测运行方法 |
CN111266032A (zh) * | 2020-02-12 | 2020-06-12 | 淮南北新建材有限公司 | 一种混合机转子及其构建方法 |
CN111266032B (zh) * | 2020-02-12 | 2022-03-25 | 淮南北新建材有限公司 | 一种混合机转子及其构建方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105956254A (zh) | 一种用于转子系统动力学建模的刚体单元法 | |
Gan et al. | Constraint-based limb synthesis and mobility-change-aimed mechanism construction | |
Gosselin et al. | Synthesis and design of reactionless three-degree-of-freedom parallel mechanisms | |
Dakel et al. | Nonlinear dynamics of a support-excited flexible rotor with hydrodynamic journal bearings | |
Hosseini-Hashemi et al. | Accurate natural frequencies and critical speeds of a rotating functionally graded moderately thick cylindrical shell | |
Zhao et al. | Deployment dynamics of a simplified spinning IKAROS solar sail via absolute coordinate based method | |
Cai et al. | The foldability of cylindrical foldable structures based on rigid origami | |
Jianguo et al. | Nonrigidly foldability analysis of Kresling cylindrical origami | |
CN109079784B (zh) | 一种基于轴不变量的多轴机器人系统建模与控制方法 | |
Jiang et al. | Trajectory generation for three-degree-of-freedom cable-suspended parallel robots based on analytical integration of the dynamic equations | |
Liu et al. | Simple formulations of imposing moments and evaluating joint reaction forces for rigid-flexible multibody systems | |
O'Shea et al. | Reference conditions and substructuring techniques in flexible multibody system dynamics | |
Chen et al. | A unified quasi-three-dimensional solution for vibration analysis of rotating pre-twisted laminated composite shell panels | |
Song et al. | Mobility analysis of the threefold-symmetric Bricard linkage and its network | |
Zhao et al. | Type synthesis and analysis of rotational parallel mechanisms with a virtual continuous axis | |
Jiang et al. | Dynamic optimization of reactionless four-bar linkages | |
CN112149245B (zh) | 一种柔性轴-盘系统耦合动力学建模及分析方法 | |
Luo et al. | Rigid-flexible coupling dynamics simulation of 3-RPS parallel robot based on ADAMS and ANSYS | |
Yang et al. | Kinematic modeling and optimization of a clustered tensegrity mobile robot | |
Zahariev | Generalized finite element approach to dynamics modeling of rigid and flexible systems | |
Fan et al. | Discretization methods of a rotating flexible rectangular thin plate | |
Jiang et al. | A four-rigid-body element model and computer simulation for flexible components of wind turbines | |
Weiler | Finite element method based analysis and modeling in rotordynamics | |
Stoykov | Nonlinear vibrations of 3D beams | |
Cornette et al. | Aeroelastic tailoring of helicopter blades |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160921 |