一种用于转子系统动力学建模的刚体单元法
【专利摘要】本发明公开一种用于转子系统动力学建模的刚体单元法，包括：1)收集转子的几何参数、材料属性和运行状态；2)将转子划分为若干个离散的刚体单元，任意相邻的两个刚体单元之间通过四个假想的弹簧相连接；3)根据转子相邻刚体单元之间相互作用模型计算作用在每个刚体单元上的合力和合力矩，并建立每个刚体单元的动力学微分方程；4)使用变步长的四阶龙格?库塔?费尔伯格方法求解各个刚体单元的动力学方程来获得刚体单元的振动响应。利用本发明建立的转子系统动力学模型，不需要任何假设，就能够与Gupta轴承动力学模型实现耦合，建立整个滚动轴承?转子系统的完全动力学模型。
【专利说明】
-种用于转子系统动力学建模的刚体单元法
技术领域
[0001] 本发明设及一种用于转子系统动力学建模的方法。
【背景技术】
[0002] 在转子动力学研究中，合理地对转子系统进行动力学建模是进行后续分析和计算 的基础。从现有文献检索发现，转子动力学建模方法主要有理论法、传递矩阵法和有限元 法，其中基于理论力学的建模方法一般只能用于简单的转子系统建模，对于复杂的转子系 统，较常用的是传递矩阵法和有限元法。传递矩阵法的主要特点是矩阵的阶数不随系统的 自由度数增大而增加，因此编程简单，所需内存小，运算速度快，但是该方法在考虑转子周 围结构时分析计算比较困难，而且传递矩阵法对转子系统进行了简化和等效，无法保证模 型的完整性和分析的准确度。
[0003] 有限元法基于转子的几何结构，将模型离散为一系列单元进行分析，其表达式简 单规范，能够同时保证模型的完整性和计算准确度，也很容易借助有限元软件实现，很好地 兼顾到模型的完整性和计算效率。另外，利用有限元分析软件可W很好地解决转子动力特 性分析中的叩它螺效应"问题，可W提高准确度和分析效率。申请号为201510816868.3的中 国发明专利申请公开了一种基于双转子简化动力学模型设计的航空发动机仿真试验台，其 特点在于首先采用有限元法建立多轮盘的复杂离散动力学模型，然后基于质屯、集中方法对 该模型进行了简化。可W看出，利用有限元方法建立的模型结构复杂，适合用于线性振动特 性的定量计算，但对于非线性问题的定量分析计算效率较低，定性分析难W实现;且存在自 由度数特别大，耗费计算机时等问题。
[0004] 此外，在研究轴承-转子系统的动力学行为时，须建立转子和滚动轴承的禪合模 型。为了考虑转子的弹性振动和柔性，目前较为理想的方法是采用转子有限元模型与Gupta 滚动轴承动力学模型实现禪合，但两者进行禪合时，为了求解有限元模型，需要采用一定的 假设，而运些假设的合理性还有待于进一步验证，因此采用转子有限元模型与Gupta轴承模 型禪合的方法具有一定的局限性。
【发明内容】
[0005] 本发明的目的在于提供一种用于转子系统动力学建模的刚体单元法:通过将转子 分隔成由弹黃相连接的离散的刚体单元，根据转子相邻刚体单元之间相互作用模型计算作 用在每个刚体单元上的合力和合力矩，并建立每个刚体单元的动力学微分方程，然后使用 变步长的四阶龙格-库塔-费尔伯格方法求解各个刚体单元的动力学方程来获得刚体单元 的振动响应。本发明提出了一种全新的转子系统动力学建模方法，与传统方法相比较，其显 著优点是:考虑了转子的弹性振动和柔性，操作方法稳定可靠、效率较高;并且利用该方法 建立的转子系统动力学模型，不需要任何假设，就能够与Gupta轴承动力学模型实现禪合， 建立整个滚动轴承-转子系统的完全动力学模型。
[0006] 为了实现上述目的，本发明采用如下技术方案：
[0007] -种用于转子系统动力学建模的刚体单元法，包括W下步骤：
[000引1)收集转子的几何参数、材料属性和运行状态；
[0009] 2)将转子划分为若干个离散的刚体单元，任意相邻的两个刚体单元之间通过四个 假想的弹黃相连接，所述四个假想的弹黃包括一个约束平移运动的拉伸弹黃和=个约束旋 转运动的扭转弹黃；
[0010] 3)根据转子各个相邻刚体单元之间相互作用模型，计算作用在每个转子刚体单元 上的合力和合力矩，并建立每个刚体单元的动力学微分方程；
[0011] 4)使用变步长的四阶龙格-库塔-费尔伯格方法求解各个刚体单元的动力学微分 方程来获得每个刚体单元的振动响应。
[0012] 所述步骤2)中，在将转子划分为若干个离散的刚体单元时，为了使建立的转子动 力学模型足够合理，并方便计算作用在刚体单元上的合力和合力矩，需要遵循W下标准：
[0013] (1)划分的每一个刚体单元的内径或外径必须一致，同一个刚体单元不能具有不 同的内径或外径；
[0014] (2)为了保证各个节点振幅相对误差在合理范围内，建议划分后的每一个刚体单 元的长径比不超过5;
[0015] (3)每个轴承在转子上的支承位置的中屯、必须与它所对应的刚体单元的节点对 齐，而不能与两个相邻刚体单元之间的面对齐。
[0016] 进一步的，所述步骤3)中，定义惯性坐标系i(xi，yi，zi)和第j个刚体单元的定体坐 标系rj(xu，yu，zu)，可通过刚体单元的姿态角确定。每个刚体单元有两个端面，左端被命 名为a平面，右端被命名为e平面。M神日的分别是第j个刚体单元a平面和0平面的中屯、。同样， Mw和Nw分别是第j+1个刚体单元a平面和e平面的中屯、。第j个刚体单元和第j+1个刚体单 元之间的几何作用关系为：当转子在自由状态时，第j个刚体单元的e平面与第j+1个刚体单 元的a平面重合，此时它们之间的弹黃处于松弛状态；当转子受到外力时，由于相邻刚体单 元运动的不一致，连接相邻刚体单元平面之间的弹黃将发生压缩或者拉伸，从而产生了相 互作用力。W第j个刚体单元的e平面和第j+1个刚体单元的a平面的相互作用为例，使用上 标i表示惯性坐标系。在惯性坐标系中，神PMw分别由向量i和喊确定。由拉伸弹黃产生 的作用在第^'个刚体单元上的相互作用力^;(^^'..'^',1|'，厂：)和作用在第^'+1个刚体单元上的 相互作用力把表示为：
[0017]
[001引式中--第j个刚体单元e面与第j + 1个刚体单元a面之间的拉伸弹黃的刚 度/N.nfi;E--转子的弹性模量/Pa; Ij和Ij+i--第j个和第j+1个刚体单元的长度/m;Aj和 Aj+1--第j个和第j+1个刚体单元端面的面积/m2;函数min(a，b)计算a和b的最小值。
[0019] 除了拉伸弹黃的相互作用力，同时还有由S个方向上的扭转弹黃产生的作用在第 j个刚体单元的作用力矩M；^'(Mf ,Mf ,Mf)和作用在第_]'"个刚体单元的作用力矩 品，M品相邻刚体单元的姿态角为扣点，、）和咕+1点+i，Vi)，相互作用 力矩Mf和尽!由下式确定：
[0020]
[0021] 式中：--第j个刚体单元e面与第j+1个刚体单元a面之间X方向扭转弹黃的 巧IJ度/N.m.racfU-转子的剪切模量/Pa; Ip麻Ipw-第j个和第j+1个刚体单元的极惯 性矩/V;
[0022] 力矩分量Mfy和M:,通过下式确定：
[0023]
[0024] AT --呆_rr刚14^平兀口凹*^呆_] 十1/|、刚14^平兀U面之间y方向扭转弹黃的 刚度/^.111^曰^1;1、1和^^--第j个刚体单元和第j+1个刚体单元相对于y轴的惯性矩/m4; [00巧]力矩分量Mf和M品表示为：
[0026]
[0027] 式中：/、' 一一第j个刚体单元0面与第j+1个刚体单元a面之间Z方向扭转弹黃的 刚度/N.m.rad-1; Izj和IzW-一第j个刚体单元和第j + 1个刚体单元相对于z轴的惯性矩/ m*。
[0028] 进一步的，步骤3)中，将每个刚体单元的运动分解成为质屯、的平动运动和绕其质 屯、的转动;巧(巧?巧，巧)和巧(巧，巧，巧)分别是由不平衡质量产生的不平衡力和施加 在第j个刚体单元上的外力，则第j个刚体单元的平动运动方程为：
[0029]
[0030] 式中吨一-第j个刚体单元的质量Ag; (?，化Zj)--质屯、Ou在惯性坐标系中的 坐标分量;Gj-一刚体单元的重力/N;不平衡力巧表示为：
[0031]
[0032] 式中:mujruj-一不平衡质量和其半径的乘积Ag.m;aj--不平衡质量的初始方位 角/rad;Tu,i-一从刚体单元定体坐标系到惯性坐标系的变换矩阵；
[0033] 第j个刚体单元的角速度为则刚体单元的旋转运动方程表示 为：
[0034]
[0035] 式中施加在第j个刚体单元上的外部力矩)和 如贫，M貧',M穿)是由F;;和与,产生的力矩，表示为：
[0036]
[0037] 式中：Tl,U-一从惯性坐标系到刚体单元定体坐标系的变换矩阵。
[0038] 所述步骤4)中，使用变步长四阶龙格-库塔-费尔伯格方法进行积分，可得到下一 时刻各刚体单元的速度和位移。重复此步骤可求得每个刚体单元的振动响应，对应于转子 不同位置的振动响应。
[0039] 本发明的进一步改进在于，转子模型的求解步骤为:首先根据转子的几何参数、材 料属性、运行状态计算各转子单元初始速度和位移，然后根据转子相邻刚体单元之间相互 作用模型计算作用在每个刚体单元上的合力和合力矩，并建立每个刚体单元的动力学微分 方程;使用变步长的四阶龙格-库塔-费尔伯格方法求解各个刚体单元的动力学方程来获得 刚体单元的振动响应。
[0040] 与现有技术相比，本发明具有下述优点：
[0041] 1、提出了一个新的转子系统动力学建模方法(DM)，将转子分隔成离散的刚体单 元，通过求解各个刚体单元的动力学方程来获得刚体单元的动力学响应；
[0042] 2、利用本发明的方法建立的转子动力学模型，不需要做任何假设，就能够与Gupta 轴承动力学模型实现禪合，建立整个轴承-转子系统的完全动力学模型。
[0043] 本发明通过将转子分隔成由弹黃相连接的离散的刚体单元，根据转子相邻刚体单 元之间相互作用模型计算作用在每个刚体单元上的合力和合力矩，并建立每个刚体单元的 动力学微分方程，使用变步长的四阶龙格-库塔-费尔伯格方法求解各个刚体单元的动力学 方程来获得刚体单元的振动响应。利用该方法建立的转子系统动力学模型，不需要任何假 设，就能够与Gupta轴承动力学模型实现禪合，建立整个滚动轴承-转子系统的完全动力学 模型。
【附图说明】
[0044] 图1是本发明中转子的任意两个相连接刚体单元之间的相互作用：坐标系分别为 惯性坐标系i(xl，yl，zl)和第j个刚体单元的定体坐标系;rj(xU，yU，zU);
[0045] 图2是本发明的流程示意图；
[0046] 图3是本发明的转子模型，其中图3(a)为转子系统动力学模型，图3(b)为转子系统 有限元模型；
[0047] 图4是本发明中节点4在y方向的时域响应图；
[0048] 图5是本发明中节点4在Z方向的时域响应图；
[0049] 图6是本发明中转子各节点的位移图；
[0050] 图7是本发明中转子转速为100(K)r/min时不同节点的轨迹。
【具体实施方式】
[0051] 请参阅图1及图2所示，本发明一种用于转子系统动力学建模的刚体单元法，包括 W下步骤：
[0052] 1)收集转子的几何参数、材料属性和运行状态，为转子系统动力学建模提供数据 支持；
[0053] 2)将转子划分为若干个离散的刚体单元，任意相邻的两个刚体单元之间通过四个 假想的弹黃相连接，即一个约束平移运动的拉伸弹黃和=个约束旋转运动的扭转弹黃；
[0054] 3)根据转子各个相邻刚体单元之间相互作用模型，计算作用在每个转子刚体单元 上的合力和合力矩，并建立每个刚体单元的动力学微分方程；
[0055] 4)使用变步长的四阶龙格-库塔-费尔伯格方法求解各个刚体单元的动力学方程 来获得每个刚体单元的振动响应。
[0056] 步骤1)中，转子的材料特性参数包括转子长度、直径，W及弹性模量、泊松比和密 度。
[0057] 步骤2)中，在将转子划分为若干个离散的刚体单元时，为了使建立的转子动力学 模型足够合理，并方便计算作用在刚体单元上的合力和合力矩，需要遵循W下标准：
[0058] (1)划分的每一个刚体单元的内径或外径必须一致，同一个刚体单元不能具有不 同的内径或外径；
[0059] (2)为了保证各个节点振幅相对误差在合理范围内，建议划分后的每一个刚体单 元的长径比不超过5;
[0060] (3)每个轴承在转子上的支承位置的中屯、必须与它所对应的刚体单元的节点对 齐，而不能与两个相邻刚体单元之间的面对齐。
[0061] 步骤4)中，使用变步长四阶龙格-库塔-费尔伯格方法进行积分，可得到下一时刻 各刚体单元的速度和位移。重复此步骤可求得每个刚体单元的振动响应。
[0062] 下面结合一个实例对本发明作进一步详细说明，但本实例并不用于限制本发明。 其步骤为：
[0063] 1)收集转子的几何参数、材料属性和运行状态，为转子系统动力学建模提供数据 支持，转子的参数如表1所示。
[0064] 表1转子参数
[00 化]
[0066] 2)将转子等分成9个刚体单元（其中第2个和第8个刚体单元通过弹黃和阻尼支 承），任意相邻的两个刚体单元之间通过四个假想的弹黃相连接，即一个约束平移运动的拉 伸弹黃和=个约束旋转运动的扭转弹黃。
[0067] 3)根据相邻刚体单元之间的相互作用模型，对转子各刚体单元之间的相互作用进 行分析，计算作用在每个转子单元上的力和力矩，建立每个刚体单元的运动微分方程，转子 动力学模型(DM)如图3(a)所示。具体计算方法如下：
[006引选取任意两个相邻刚体单元(第j个和第j+1个），定义惯性坐标系i(xi，yi，zi)和第 j个刚体单元的定体坐标系rj(x",y"，z")，可通过刚体单元的姿态角确定。每个刚体单元 的左端面被命名为a平面，右端面被命名为e平面。心和的分别是第j个刚体单元a平面和0平 面的中屯、。同样，Mj+i和Nj+i分别是第j+1个刚体单元a平面和0平面的中屯、。使用上标i表示惯 性坐标系，在惯性坐标系中，的和Mw分别由向量^和^^,?确定。当转子受到外力时，由拉伸 弹黃产生的作用在第j个刚体单元上的相互作用力厂'，厂-和作用在第j + 1个刚体 单元上的相互作用力，厂I、F-~ )可根据下式来计算 財\ 卽1 a州. a江I /
[0069]
[0070] 式中：--第j个刚体单元e面与第j + 1个刚体单元a面之间的拉伸弹黃的刚 度/N.nfi;E--转子的弹性模量/Pa; Ij和Ij+i--第j个和第j+1个刚体单元的长度/m;Aj和 Ajn一一第j个和第j+1个刚体单元端面的面积/V。函数min(a，b)计算a和b的最小值。除了 拉伸弹黃的相互作用力，同时还有由=个方向上的扭转弹黃产生的作用在第j个刚体 单元的作用力矩，Mf)和作用在第j + 1个刚体单元的作用力矩 前三品，W品，始品相邻刚体单元的姿态角为咕点，、)和咕+1点+1，、+1)，则相互作 用力矩Mf和M品可由下式确定
[0071:
[0072] 式中：一一第j个刚体单元e面与第j+1个刚体单元a面之间X方向扭转弹黃的 刚度/N.m.rad i，G--转子的剪切模量/Pa; Ipj和Ipj+i--第j个和第j+1个刚体单兀的极惯 性矩/V。
[0073] 力矩分量敵尸和姐|：；可^通过下式确定
[0074]
[00巧]式中：--第j个刚体单元e面与第j+1个刚体单元a面之间y方向扭转弹黃的 刚度/N.m.racri;^神日-一第j个刚体单元和第j+1个刚体单元相对于y轴的惯性矩/m4。
[0076] 力矩分量Mf和M品可表示为
[0077]
[0078] 式中：-一第j个刚体单元e面与第j+1个刚体单元a面之间Z方向扭转弹黃的 刚度/N.m.racTi;Iz神日IzW-一第j个刚体单元和第j+1个刚体单元相对于z轴的惯性矩/m4。
[0079] 每个刚体单元的运动可W分解成为质屯、的平动运动和绕其质屯、的转动。假设 巧(巧，巧，巧)和巧(巧，巧，巧)分别是由不平衡质量产生的不平衡力和施加在第j个刚体 单元上的外力，则第j个刚体单元的平动运动方程为
[0080]
[0081 ]式中:mj--第j个刚体单元的质量Ag; (Xj，yj，zj)--质屯、Orj在惯性坐标系中的 坐标分量;Gj-一刚体单元的重力/N。不平衡力巧可W表示为
[0082
[0083] 式中:mujruj-一不平衡质量和其半径的乘积Ag.m;aj--不平衡质量的初始方位 角/rad;Tu,i-一从刚体单元定体坐标系到惯性坐标系的变换矩阵。
[0084] 假设第j个刚体单元的角速度为COjx, 则刚体单元的旋转运动方程 可W表示为
[0085]
[0086] 式中施加在第j个刚体单元上的外部力矩M克晋，M备，M篇;)和 貧，M貧，M貧)是由捉.和巧,产生的力矩，可W表示为
[0087]
[0088] 式中：Tl,U-一从惯性坐标系到刚体单元定体坐标系的变换矩阵。
[0089] 4)使用变步长四阶龙格-库塔-费尔伯格方法求解各个刚体单元的动力学方程，最 大允许截断误差被设定为1.0E-5,可W得到任意一个或所关注刚体单元的振动响应。
[0090] 本实例取转子动力学模型的节点4为研究对象，当转子转速为lOOOr/min时，在节 点4的Z方向施加一个大小为IOON的正弦激振力，通过仿真分析，可W得到它的y方向和Z方 向的时域响应，分别如图4和图5所示，可W看出节点4在y方向和Z方向的时域响应曲线均为 正弦曲线。重复W上方法，就可W求解任意一个或所关注刚体单元的振动响应。
[0091] 另外，图3(b)给出了包含10个单元的转子的有限元模型(FEM)，第1个和第9个单元 的长度均为20mm，其他单元的长度都为40mm。转子的参数如表1示所列。有限元模型使用 化wmarkP方法求解，计算时间步长设定为2.0E-4s。
[0092] 1)假设当转子转速为0时，在节点4的Z方向施加一个径向力，图6为不同大小的径 向力作用下转子各节点在Z方向的位移。从图6可W看出当转子受到静载荷时，FEM的计算结 果和DM的计算结果匹配较好。
[0093] 2)假设在节点4处存在不平衡质量，不平衡质量和其半径的乘积是0.001kg.m。当 转子旋转时，转子在不平衡力作用下会产生振动响应。图7是当转速为100(K)r/min时各个节 点的轨迹（图中：FEM-一有限元模型计算的结果，DM-一动力学模型计算的结果），由于支 承弹黃的限制，节点2和节点8的振幅要小于其他节点。
【主权项】
1. 一种用于转子系统动力学建模的刚体单元法，其特征在于，包括以下步骤： 1) 收集转子的几何参数、材料属性和运行状态； 2) 将转子划分为若干个离散的刚体单元，任意相邻的两个刚体单元之间通过四个假想 的弹簧相连接，所述四个假想的弹簧包括一个约束平移运动的拉伸弹簧和三个约束旋转运 动的扭转弹簧； 3) 根据转子各个相邻刚体单元之间相互作用模型，计算作用在每个转子刚体单元上的 合力和合力矩，并建立每个刚体单元的动力学微分方程； 4) 使用变步长的四阶龙格-库塔-费尔伯格方法求解各个刚体单元的动力学微分方程 来获得每个刚体单元的振动响应。2. 根据权利要求1所述的一种用于转子系统动力学建模的刚体单元法，其特征在于，所 述步骤3)中，定义惯性坐标系i (X1，y1，z1)和第j个刚体单元的定体坐标系r j (X1··3，yU，z1··3)， 每个刚体单元有两个端面，左端被命名为α平面，右端被命名为β平面，MjPNj分别是第j个刚 体单兀α平面和β平面的中心，同样，Mj+i和Nj+i分别是第j+Ι个刚体单兀α平面和β平面的中 心，第j个刚体单元和第j+Ι个刚体单元之间的几何作用关系为：当转子在自由状态时，第j 个刚体单元的β平面与第j+1个刚体单元的α平面重合，此时它们之间的弹簧处于松弛状态； 当转子受到外力时，由于相邻刚体单元运动的不一致，连接相邻刚体单元平面之间的弹簧 将发生压缩或者拉伸，从而产生了相互作用力；使用上标i表示惯性坐标系，在惯性坐标系 中，N#PMj+1分别由向量<和4/+|确定，由拉伸弹簧产生的作用在第j个刚体单元上的相互作 用力4(5,5,/^和作用在第_]_+1个刚体单元上的相互作用力€,+1|< 1，<1,/^)表示 为：式中-一第j个刚体单元β面与第j + Ι个刚体单元α面之间的拉伸弹簧的刚度/ N.nfSE--转子的弹性模量/Pa;lj和lj+1--第j个和第j + Ι个刚体单元的长度/m;Aj和 Aj+i--第j个和第j+Ι个刚体单元端面的面积/m2;函数min(a,b)计算a和b的最小值； 除了拉伸弹簧的相互作用力，同时还有由三个方向上的扭转弹簧产生的作用在第j个 刚体单元的作用力矩和作用在第j + Ι个刚体单元的作用力矩 相邻刚体单元的姿态角为和(rb+i，L+i，Vi)，相互作用 力矩Mf和由下式确定：式中：/(^ν,--第j个刚体单元β面与第j+1个刚体单元α面之间χ方向扭转弹簧的刚度/ N.m.racf^G-一转子的剪切模量/Pa; IPjPIpj+1-一第j个和第j+1个刚体单元的极惯性矩/ m4; 力矩分量Tkff和通过下式确定：式中：巧U+1--第j个刚体单元β面与第j+1个刚体单元α面之间y方向扭转弹簧的刚度/ N.m.racTSlyj和Iyj+i--第j个刚体单元和第j+1个刚体单元相对于y轴的惯性矩/m4; 力矩分量Mf和表示为：式中：--第j个刚体单元β面与第j+1个刚体单元α面之间z方向扭转弹簧的刚度/ N.m.rad-1 ;Izj和Izj+i--第j个刚体单元和第j+1个刚体单元相对于z轴的惯性矩/m4。3.根据权利要求1所述的一种用于转子系统动力学建模的刚体单元法，其特征在于，步 骤3)中，将每个刚体单元的运动分解成为质心的平动运动和绕其质心的转动； /^/^/^/^和^^/^/^/^分别是由不平衡质量产生的不平衡力和施加在第^冲刚 体单元上的外力，则第j个刚体单元的平动运动方程为：式中:mj--第j个刚体单元的质量/kg; (xj，yj，zj)--质心Orj在惯性坐标系中的坐标 分量一一刚体单元的重力/N;不平衡力i?表示为：式中：mUjrUj--不平衡质量和其半径的乘积/kg. m;aj--不平衡质量的初始方位角/ rad;Trj>1一一从刚体单元定体坐标系到惯性坐标系的变换矩阵； 第j个刚体单元的角速度为ω」（ω jx, ω jy, ω jz)，则刚体单元的旋转运动方程表示为： ./々 j- \ .ja jy / jy jj -1 j i .f 式中施加在第j个刚体单元上的外部力矩和是由产生的力矩，表示为：式中:T1>rj一一从惯性坐标系到刚体单元定体坐标系的变换矩阵。4. 根据权利要求1所述的一种用于转子系统动力学建模的刚体单元法，其特征在于，步 骤4)中，使用变步长四阶龙格-库塔-费尔伯格方法进行积分，得到下一时刻各刚体单元的 速度和位移;重复此步骤求得每个刚体单元的振动响应，对应于转子不同位置的振动响应。5. 根据权利要求1所述的一种用于转子系统动力学建模的刚体单元法，其特征在于，步 骤2)中，采用如下标准将转子划分为若干个离散的刚体单元： (1) 同一个刚体单元不能具有不同的内径，同时也不能具有不同的外径； (2) 划分后的每一个刚体单元的长径比不超过5; (3) 每个轴承在转子上的支承位置的中心必须与它所对应的刚体单元的节点对齐，而 不能与两个相邻刚体单元之间的面对齐。
【文档编号】G06F17/50GK105956254SQ201610272794
【公开日】2016年9月21日
【申请日】2016年4月27日
【发明人】曹宏瑞, 魏江, 陈雪峰, 张兴武, 李亚敏
【申请人】西安交通大学