CN112487575B - 双转子系统高压转子及低压转子连接节点的检查方法 - Google Patents

双转子系统高压转子及低压转子连接节点的检查方法 Download PDF

Info

Publication number
CN112487575B
CN112487575B CN202011335291.1A CN202011335291A CN112487575B CN 112487575 B CN112487575 B CN 112487575B CN 202011335291 A CN202011335291 A CN 202011335291A CN 112487575 B CN112487575 B CN 112487575B
Authority
CN
China
Prior art keywords
rotor
shaft
pressure rotor
node
low
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
Application number
CN202011335291.1A
Other languages
English (en)
Other versions
CN112487575A (zh
Inventor
陈茉莉
罗贵火
陈丽佳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN202011335291.1A priority Critical patent/CN112487575B/zh
Publication of CN112487575A publication Critical patent/CN112487575A/zh
Application granted granted Critical
Publication of CN112487575B publication Critical patent/CN112487575B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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
    • 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)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Turbine Rotor Nozzle Sealing (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明公开了一种双转子系统高压转子及低压转子连接节点的检查方法。主要步骤是构建有轴段部件、节点信息的轴段矩阵；检查所有轴段的连接关系，分解出高压和低压转子几何模型组件；利用低压转子长度大于高压转子长度特点，区分高压转子及低压转子轴组件；结合高压转子及低压转子轴组件区分结果，利用中介轴承单元两个连接节点分别处于高压和低压轴组件特点，实现了连接节点的几何检查。本发明对不含静子部件、转子轴向方向平行于任一笛卡尔坐标系方向、中介轴承连接的双转子模型，在高压转子及低压转子连接节点被设定后，通过几何检查，可自动分解高压转子及低压转子轴组件，判断中介轴承连接节点的合理性。

Description

双转子系统高压转子及低压转子连接节点的检查方法
技术领域
本发明涉及转子动力学特性有限元建模领域，具体涉及以中介轴承为连接部件的高低压双转子系统的连接节点建模及几何检查方法。
背景技术
中介轴承是一种广泛使用的双转子航空发动机连接部件，大多数涡扇、涡喷发动机都采用中介轴承连接高压转子及低压转子，以达到缩短支点距离、降低发动机重量目的。
在有限元建模中，中介轴承是独立的连接单元模型，每个中介轴承连接单元由高压转子及低压转子的连接节点成对联接而成。一个双转子系统模型可以有多个中介轴承单元，一个中介轴承可以有多个连接节点对，但有限元模型不会对中介轴承的个数进行区分，它们在模型中都表现为连接节点对。
中介轴承单元模型可以是弹簧单元模型、弹簧-阻尼单元模型以及接触单元模型等。无论哪种单元模型，附属于中介轴承单元的高压转子及低压转子的连接节点都必须手动设置。这导致连接节点经常因人为疏忽而产生错误。
在手动设置中介轴承单元节点时，中介轴承单元节点分别属于高压转子和低压转子轴段部件，单个或多个联接的轴段部件，可形成高压转子或低于转子轴组件。两个转子轴组件中的中介轴承连接节点对应的转速、部件单元质量、刚度及陀螺矩阵均不同，又因中介轴承单元模型由连接节点对组成，这要求高压转子及低压转子连接节点的几何建模具备防错能力。
同时，中介轴承是高压转子及低压转子非线性振动的影响因素之一，连接节点几何参数设置不当可能会引起非线性微分方程无法收敛、结果发散、计算结果不可靠等问题。对复杂的双转子系统，连接节点的几何属性如所属部件、几何坐标、单元信息等应仔细检查，避免连接节点几何参数设置不当引起的有限元模型缺陷。
目前的转子动力学有限元软件如Ansys、Samcef、HyperMesh等，均具有一定的几何检查功能，但几何检查主要针对单元网格，目的是优化梁、面及体这类部件单元的尺寸、拓扑结构等，使单元形状匹配单元插值形函数，获得更为精确且收敛的计算结果。而转子中介轴承单元为连接单元，不存在网格划分及拓扑结构，因此多数商业软件无专门针对中介轴承的几何检查功能。
作为连接部件，中介轴承几何检查前，需对双转子结构的高、低压转子组件分解。目前的商业软件通过设定高、低压转子单元组件，将高、低压转子在模型上区分，从而便于不同转速载荷的施加。但是，组件的构成需手动设定，一个部件可隶属于多个组件，也可以施加多种转速载荷，导致高、低压转子单元组件易因人为疏忽分解错误，而目前商业软件对含中介轴承的双转子结构高、低压转子组件的分解合理与否无相应的几何检查。
因此，对这类广泛使用的、含中介轴承的航空发动机双转子系统，在高压转子及低压转子及中介轴承有限元单元模型建立后，有必要寻找一种几何检查方法，自动分解高、低压转子组件，并基于高压转子及低压转子组件判断中介轴承连接节点的合理性，但目前未有较好的相关几何检查方法的技术方案所公开。
发明内容
发明目的：本发明要解决的问题是提供一种双转子结构高压转子及低压转子连接节点的几何检查方法，解决如何能够自动分解高、低压转子轴组件，判断中介轴承连接节点的合理性的技术问题。
技术方案：为实现上述目的，本发明可采用以下技术方案。
一种双转子系统高压转子及低压转子连接节点的检查方法，检查对象为不含静子部件、转子轴向方向平行于任一笛卡尔坐标系方向的双转子模型的高压转子及低压转子间连接节点；该检查方法包括以下步骤：
(1)导出双转子结构的网格模型数据，其中含有高压转子及低压转子的轴段部件、单元、节点及中介轴承数据；
轴段部件数据包括轴段部件编号，各轴段部件上的所有单元编号以及每个单元上的所有节点编号；
节点数据包括节点编号及其六自由度坐标值，六自由度坐标值依次为节点在笛卡尔坐标系OXYZ中X、Y、Z三个方向的坐标值以及节点绕X、Y、Z三个方向的角度坐标值；
中介轴承数据包括中介轴承编号及各中介轴承上两个连接节点编号；
(2)从轴段部件数据中，依次读取各个轴段部件编号及所含单元，然后根据单元数据读取该轴段的节点数据，提取出轴段部件编号、轴段部件包含的所有单元编号、单元包含的所有节点编号，所有节点的坐标值，构造轴段矩阵P，轴段矩阵P表达式为
Figure BDA0002797021430000031
q为属于高压转子及低压转子轴组件的轴段部件总个数；
pi为轴段部件编号，i＝1,2,…,q；
Pi为轴段部件pi的矩阵；
I为元素全部为1的列向量；
Li为轴段部件pi所包含的单元总数；
eij为轴段部件pi的第j个单元编号，j＝1,2,…,Li；
mi为轴段部件pi单元节点数，特别说明，同一轴段部件的单元节点数相同。
nijk为单元eij的第k个节点编号，k＝1,2,…,mi；
Nijk为节点nijk的坐标值，它是1×6的行向量，
Figure BDA0002797021430000032
xijk、yijk、zijk依次为节点nijk在笛卡尔坐标系OXYZ中X、Y、Z三个方向的坐标值；θijk、/>
Figure BDA0002797021430000033
ψijk依次为节点绕X、Y、Z三个方向的角度坐标值，在建模未加载时，有/>
Figure BDA0002797021430000034
(3)令i＝1，读取矩阵P中的子矩阵P1中的节点编号，形成节点数据集合A0＝{n1jk|j＝1,2,…,L1，k＝1,2,…,m1}及轴段部件编号集合B0＝{p1}；
(4)令B＝B0；构造另一轴段部件编号集合C0，令
Figure BDA0002797021430000035
Figure BDA0002797021430000036
表示空序列，i＝2；
(5)读取子矩阵Pi；
(6)令
Figure BDA0002797021430000041
(7)若轴段Pi存在某一节点nijk，有nijk∈A0，则构造集合
Figure BDA0002797021430000042
Li，k＝1,2,…,mi}；令/>
Figure BDA0002797021430000043
否则，/>
Figure BDA0002797021430000044
(8)令i＝i+1，A0＝A，B0＝B，C0＝C；重复步骤(5)和(7)，直至完成子矩阵Pq的读取；
(9)令i＝1，构造集合D和E，令D＝C，E＝B，
Figure BDA0002797021430000045
(10)若D中的第i元素所代表的轴段编号为di，对轴段di对应的子矩阵
Figure BDA0002797021430000046
重复步骤(5)～(8)，直至i>card(D)，完成D中最后一个元素对应编号的轴段读取，其中card(D)表示集合D的长度；
(11)判断E＝B是否成立，若不成立，转入步骤(12)，否则转入步骤(13)；
(12)重复步骤(9)～(11)；
(13)对集合B中的轴段，在所有轴段中查找X、Y、Z三个方向最大和最小坐标值，分别记为xBmax、xBmin、yBmax、yBmin、zBmax、zBmin，计算集合B中轴段在三个方向的跨度，分别为xB0＝xBmax-xBmin、yB0＝yBmax-yBmin、zB0＝zBmax-zBmin，求取三个方向的最大跨度szB＝max{xB0,yB0,zB0}；
(14)对集合D中的轴段，在所有轴段中查找X、Y、Z三个方向最大和最小坐标值，分别记为xDmax、xDmin、yDmax、yDmin、zDmax、zDmin，计算集合D中轴段在三个方向的跨度，分别为xD0＝xDmax-xDmin，yD0＝yDmax-yDmin，zD0＝zDmax-zDmin，求取三个方向的最大跨度szD＝max{xD0,yD0,zD0}；
(15)若szB>szD，则B为低压转子轴段集合，D为高压转子轴段集合；反之则B为高压转子轴段集合，D为低压转子轴段集合；输出轴段编号集合B和D，并标识出高压转子轴段及低压转子轴段；
(16)依次导出各中介轴承的两个连接节点，形成二维连接节点对矩阵K，表示为
Figure BDA0002797021430000047
/>
c为中介轴承的个数；hv、lv分别表示编号为v的中介轴承上的两个连接节点编号，称它们为第v个连接节点对，v＝1,2,…,c；
(17)令v＝1，构造不合理的连接节点对编号集合R和R0，令
Figure BDA0002797021430000048
(18)对中介轴承的连接节点hv和lv，组成集合{hv，lv}，若{hv，lv}∈B或{hv，lv}∈D，则令
Figure BDA0002797021430000049
(19)令v＝v+1，R0＝R；
(20)若v>c，判断集合
Figure BDA0002797021430000051
是否为真，若是，则连接节点合理，反之输出R中的元素，R中的元素是不合理连接节点对编号；若v≤c，重复步骤(18)和步骤(19)。
有益效果：对不含静子部件、转子轴向方向平行于任一笛卡尔坐标系方向的双转子模型，提供一种双转子结构高压转子及低压转子连接节点的几何检查方法。该方法专用于中介轴承连接的双转子结构建模，在连接节点被设定后，通过几何检查，可自动分解高、低压转子轴组件，判断中介轴承连接节点的合理性。
进一步的，中介轴承连接高压转子及低压转子，高压转子及低压转子的轴段部件无公共节点，利用了这一特性对高压转子及低压转子连接节点的合理性进行检查。
进一步的，以公共节点判断轴段部件的连接关系，利用各转子(高压或低压转子)包含的轴段部件相互联接、具备公共节点这一特点，通过检查所有轴段模型中单元的连接关系，将整个双转子系统模型分解为独立的高压和低压转子几何组件模型，便于各转子转速等载荷的施加。
进一步的，考虑了含中介轴承的航空发动机双转子结构特点，利用转子轴组件结构的轴向尺寸大于径向尺寸、低压转子长度大于高压转子长度这两个特性，实现了高、低压转子轴组件的区分，获得了各转子的轴段编号集合。
进一步的，结合高、低压转子轴组件区分结果，利用中介轴承单元两个连接节点分别处于高压和低压轴组件这一特性，实现了连接节点的几何检查，判断了连接节点的合理性，并输出了不合理连接节点对(单元)编号。
附图说明
图1为高低压转子与中介轴承连接关系；
图2为单个转子轴组件与轴段部件关系图，其中A1、A2为轴段部件，整体为轴组件；
图3为轴段部件与单元的关系；
图4为单元与节点的关系；
图5为几何检查简图；
图6为几何检查详细流程图；
图7为双转子网格模型数据示意图。
具体实施方式
本实施例为高压转子及低压转子间连接节点的几何检查方法。
请参考图1所示建立高低压双转子及中介轴承连接模型。转子轴向方向平行于X方向、低压转子含2个轴段部件、高压转子含2个轴段部件、仅含1个中介轴承。
轴段部件为单个转子轴组件上的各个轴段。例如，图2中低压转子轴组件包含两个轴段部件A1和A3，其中字母“A”表征轴段部件，其后的数字表示轴段部件编号。
每个轴段部件包含若干个单元。例如，图3中轴段部件A1包含E1、E2两个单元，轴段部件A3包含E3、E4、E5、E6四个单元，其中字母“E”表征单元，其后的数字表示单元编号。
每个单元包含若干个节点。例如，图4所示的单元E2包含N2和N3两个节点，其中字母“N”表征节点，其后的数字表示节点编号。
高压转子及低压转子间连接节点的几何检查方法流程简图及详细图示分别如图5和图6所示。具体步骤如下：
(1)请参考图7所示，导出双转子结构的网格模型数据，其中含有轴段部件、单元、节点及中介轴承连接节点数据信息。
轴段部件数据包括轴段部件编号，各轴段部件上的所有单元编号以及每个单元上的所有节点编号；
节点数据包括节点编号及其六自由度坐标值，六自由度坐标值依次为节点在笛卡尔坐标系OXYZ中X、Y、Z三个方向的坐标值以及节点绕X、Y、Z三个方向的角度坐标值；
中介轴承数据包括中介轴承编号及各中介轴承上两个连接节点编号。
(2)从轴段部件数据中，依次读取各个轴段部件编号及所含单元，然后根据单元信息读取该轴段的节点信息。提取出轴段部件编号、轴段部件包含的所有单元编号、单元包含的所有节点编号，所有节点的坐标值，构造轴段矩阵P，有
Figure BDA0002797021430000061
Figure BDA0002797021430000062
Figure BDA0002797021430000071
/>
Figure BDA0002797021430000072
Figure BDA0002797021430000073
其中，对应于发明内容中的步骤(2)，有：
q＝4；
p1＝1，p2＝2，p3＝3，p4＝4；
L1＝2，L2＝5，L3＝4，L4＝3；
e11＝1，e12＝2；
e21＝7，e22＝8，e23＝9，e24＝10，e25＝11；
e31＝3，e32＝4，e33＝5，e34＝6；
e41＝12，e42＝13，e43＝14；
m1＝m2＝m3＝m4＝2；
nijk为单元eij的第k个节点，k＝1,2，如n311＝3；
Nijk为节点nijk的坐标值，如N311＝{0.65，0，0，0，0，0}。
(3)令i＝1，读取矩阵P中的子矩阵P1中的节点编号，形成节点数据集合A0＝{1,2,3}及轴段部件编号集合B0＝{1}。
(4)令B＝B0；构造另一轴段部件编号集合C0，令
Figure BDA0002797021430000081
Figure BDA0002797021430000082
表示空序列，i＝2。
(5)读取子矩阵Pi。
(6)令
Figure BDA0002797021430000083
(7)若轴段Pi存在某一节点nijk，有nijk∈A0，则构造集合
Figure BDA0002797021430000084
Figure BDA0002797021430000085
令/>
Figure BDA0002797021430000086
否则，/>
Figure BDA0002797021430000087
(8)令i＝i+1，A0＝A，B0＝B，C0＝C；重复步骤(5)和(7)，直至完成子矩阵P4的读取。
说明：完成P4读取后，有A0＝{1,2,3，4,5,6,7}，B0＝{1,3}，C0＝{2,4}
(9)令i＝1，构造轴段部件编号集合D和E，令D＝C，E＝B，
Figure BDA0002797021430000088
(10)若D中的第i元素所代表的轴段编号为di，对di对应的轴段
Figure BDA0002797021430000089
重复步骤(5)～(8)，直至i>2。
(11)判断E＝B是否成立。因成立，转入如下步骤(12)。
(12)对集合B中的轴段，在所有轴段中查找X、Y、Z三个方向最大和最小坐标值，得
xBmax＝1.5、xBmin＝0、yBmax＝yBmin＝0、zBmax＝zBmin＝0
计算集合B中轴段在三个方向的跨度，得xB0＝1.5、zB0＝yB0＝0。因此三个方向的最大跨度szB＝1.5。
(13)对集合D中的轴段，在所有轴段中查找X、Y、Z三个方向最大和最小坐标值，得
xDmax＝0.85、xDmin＝0.45、yDmax＝yDmin＝0、zDmax＝zDmin＝0
计算集合D中轴段在三个方向的跨度，分别为xD0＝0.4，yD0＝zD0＝0。因此三个方向的最大跨度szD＝0.4。
(14)因szB>szD，故B为低压转子轴段，D为高压转子轴段。输出轴段编号集合B和D，并标识出高压转子轴段及低压转子轴段。
(15)请参考图7，导出中介轴承的连接节点对，构造二维矩阵K，表示为
K＝[3 12]
对应于发明内容中的步骤(16)，有c＝1，h1＝3，l1＝12。
(16)令v＝1，构造不合理的连接节点对编号集合R和R0，令
Figure BDA0002797021430000091
(17)对中介轴承的连接节点hv和lv，组成集合{hv，lv}，若{hv，lv}∈B或{hv，lv}∈D，则令
Figure BDA0002797021430000092
(18)令v＝v+1，R0＝R。
(19)因v>1，集合
Figure BDA0002797021430000093
为真，连接节点合理。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (5)

1.一种双转子系统高压转子及低压转子连接节点的检查方法，其特征在于，检查对象为不含静子部件、转子轴向方向平行于任一笛卡尔坐标系方向的双转子模型的高压转子及低压转子间连接节点；该检查方法包括以下步骤：
(1)导出双转子结构的网格模型数据，其中含有高压转子及低压转子的轴段部件、单元、节点及中介轴承数据；
轴段部件数据包括轴段部件编号，各轴段部件上的所有单元编号以及每个单元上的所有节点编号；
节点数据包括节点编号及其六自由度坐标值，六自由度坐标值依次为节点在笛卡尔坐标系OXYZ中X、Y、Z三个方向的坐标值以及节点绕X、Y、Z三个方向的角度坐标值；
中介轴承数据包括中介轴承编号及各中介轴承上两个连接节点编号；
(2)从轴段部件数据中，依次读取各个轴段部件编号及所含单元，然后根据单元数据读取该轴段的节点数据，提取出轴段部件编号、轴段部件包含的所有单元编号、单元包含的所有节点编号，所有节点的坐标值，构造轴段矩阵P，轴段矩阵P表达式为
Figure FDA0002797021420000011
q为属于高压转子及低压转子轴组件的轴段部件总个数；
pi为轴段部件编号，i＝1,2,…,q；
Pi为轴段部件pi的矩阵；
I为元素全部为1的列向量；
Li为轴段部件pi所包含的单元总数；
eij为轴段部件pi的第j个单元编号，j＝1,2,…,Li；
mi为轴段部件pi单元节点数，同一轴段部件的单元节点数相同；
nijk为单元eij的第k个节点编号，k＝1,2,…,mi；
Nijk为节点nijk的坐标值，它是1×6的行向量，
Figure FDA0002797021420000021
xijk、yijk、zijk依次为节点nijk在笛卡尔坐标系OXYZ中X、Y、Z三个方向的坐标值；θijk、/>
Figure FDA0002797021420000022
ψijk依次为节点绕X、Y、Z三个方向的角度坐标值，在建模未加载时，有/>
Figure FDA0002797021420000023
/>
(3)令i＝1，读取矩阵P中的子矩阵P1中的节点编号，形成节点数据集合A0＝{n1jk|j＝1,2,…,L1，k＝1,2,…,m1}及轴段部件编号集合B0＝{p1}；
(4)令B＝B0；构造另一轴段部件编号集合C0，令
Figure FDA0002797021420000024
Figure FDA0002797021420000025
表示空序列，i＝2；
(5)读取子矩阵Pi；
(6)令
Figure FDA0002797021420000026
(7)若轴段Pi存在某一节点nijk，有nijk∈A0，则构造集合
Figure FDA0002797021420000027
Li，k＝1,2,…,mi}；令/>
Figure FDA0002797021420000028
否则，/>
Figure FDA0002797021420000029
(8)令i＝i+1，A0＝A，B0＝B，C0＝C；重复步骤(5)和(7)，直至完成子矩阵Pq的读取；
(9)令i＝1，构造集合D和E，令D＝C，E＝B，
Figure FDA00027970214200000210
(10)若D中的第i元素所代表的轴段编号为di，对轴段di对应的子矩阵
Figure FDA00027970214200000211
重复步骤(5)～(8)，直至i>card(D)，完成D中最后一个元素对应编号的轴段读取，其中card(D)表示集合D的长度；
(11)判断E＝B是否成立，若不成立，转入步骤(12)，否则转入步骤(13)；
(12)重复步骤(9)～(11)；
(13)对集合B中的轴段，在所有轴段中查找X、Y、Z三个方向最大和最小坐标值，分别记为xBmax、xBmin、yBmax、yBmin、zBmax、zBmin，计算集合B中轴段在三个方向的跨度，分别为xB0＝xBmax-xBmin、yB0＝yBmax-yBmin、zB0＝zBmax-zBmin，求取三个方向的最大跨度szB＝max{xB0,yB0,zB0}；
(14)对集合D中的轴段，在所有轴段中查找X、Y、Z三个方向最大和最小坐标值，分别记为xDmax、xDmin、yDmax、yDmin、zDmax、zDmin，计算集合D中轴段在三个方向的跨度，分别为xD0＝xDmax-xDmin，yD0＝yDmax-yDmin，zD0＝zDmax-zDmin，求取三个方向的最大跨度szD＝max{xD0,yD0,zD0}；
(15)若szB>szD，则B为低压转子轴段集合，D为高压转子轴段集合；反之则B为高压转子轴段集合，D为低压转子轴段集合；输出轴段编号集合B和D，并标识出高压转子轴段及低压转子轴段；
(16)依次导出各中介轴承的两个连接节点，形成二维连接节点对矩阵K，表示为
Figure FDA0002797021420000031
c为中介轴承的个数；hv、lv分别表示编号为v的中介轴承上的两个连接节点编号，称它们为第v个连接节点对，v＝1,2,…,c；
(17)令v＝1，构造不合理的连接节点对编号集合R和R0，令
Figure FDA0002797021420000032
(18)对中介轴承的连接节点hv和lv，组成集合{hv，lv}，若，则令
Figure FDA0002797021420000033
Figure FDA0002797021420000034
(19)令v＝v+1，R0＝R；
(20)若v>c，判断集合
Figure FDA0002797021420000035
是否为真，若是，则连接节点合理，反之输出R中的元素，R中的元素是不合理连接节点对编号；若v≤c，重复步骤(18)和步骤(19)。
2.根据权利要求1所述的检查方法，其特征在于：中介轴承连接高压转子及低压转子，高压转子及低压转子的轴段部件无公共节点。
3.根据权利要求1或2所述的检查方法，其特征在于：以公共节点判断轴段部件的连接关系，高压或低压转子包含的轴段部件相互联接并具备公共节点。
4.根据权利要求1所述的检查方法，其特征在于：转子轴组件结构的轴向尺寸大于径向尺寸、低压转子长度大于高压转子长度。
5.根据权利要求4所述的检查方法，其特征在于：中介轴承单元两个连接节点分别处于高压和低压轴组件。
CN202011335291.1A 2020-11-25 2020-11-25 双转子系统高压转子及低压转子连接节点的检查方法 Active CN112487575B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011335291.1A CN112487575B (zh) 2020-11-25 2020-11-25 双转子系统高压转子及低压转子连接节点的检查方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011335291.1A CN112487575B (zh) 2020-11-25 2020-11-25 双转子系统高压转子及低压转子连接节点的检查方法

Publications (2)

Publication Number Publication Date
CN112487575A CN112487575A (zh) 2021-03-12
CN112487575B true CN112487575B (zh) 2023-06-06

Family

ID=74934325

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011335291.1A Active CN112487575B (zh) 2020-11-25 2020-11-25 双转子系统高压转子及低压转子连接节点的检查方法

Country Status (1)

Country Link
CN (1) CN112487575B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105278349A (zh) * 2015-11-23 2016-01-27 哈尔滨工业大学 基于双转子简化动力学模型设计的航空发动机仿真试验台
CN105469316A (zh) * 2014-08-29 2016-04-06 华北电力大学 配电网任意两点间理论线损计算方法及系统
CN107133377A (zh) * 2017-03-28 2017-09-05 上海数设科技有限公司 螺栓连接的建模方法和装置
CN108197360A (zh) * 2017-12-22 2018-06-22 哈尔滨汽轮机厂有限责任公司 汽轮机转子网格自动划分系统和方法
CN111209639A (zh) * 2020-02-17 2020-05-29 合肥工业大学 一种叶轮-轴承-转子系统的高效定量建模方法
CN111462838A (zh) * 2020-03-25 2020-07-28 华南理工大学 一种直接将图像像素转换成有限元单元的方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105469316A (zh) * 2014-08-29 2016-04-06 华北电力大学 配电网任意两点间理论线损计算方法及系统
CN105278349A (zh) * 2015-11-23 2016-01-27 哈尔滨工业大学 基于双转子简化动力学模型设计的航空发动机仿真试验台
CN107133377A (zh) * 2017-03-28 2017-09-05 上海数设科技有限公司 螺栓连接的建模方法和装置
CN108197360A (zh) * 2017-12-22 2018-06-22 哈尔滨汽轮机厂有限责任公司 汽轮机转子网格自动划分系统和方法
CN111209639A (zh) * 2020-02-17 2020-05-29 合肥工业大学 一种叶轮-轴承-转子系统的高效定量建模方法
CN111462838A (zh) * 2020-03-25 2020-07-28 华南理工大学 一种直接将图像像素转换成有限元单元的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
刘玉 ; 王建军 ; 陈立强 ; 王志 ; .低压转子初始弯曲双转子系统跨坐标系建模方法.振动与冲击.2020,(第03期),全文. *
罗贵火 ; 周海仑 ; 王飞 ; 杨喜关 ; .含滚动轴承的同向和反向旋转双转子系统动力学响应.航空动力学报.2012,(第08期),全文. *
胡绚 ; 罗贵火 ; 高德平 ; .反向旋转双转子系统动力特性分析.现代机械.2007,(第04期),全文. *

Also Published As

Publication number Publication date
CN112487575A (zh) 2021-03-12

Similar Documents

Publication Publication Date Title
US6606612B1 (en) Method for constructing composite response surfaces by combining neural networks with other interpolation or estimation techniques
CN105808829A (zh) 一种基于cpu+gpu异构并行计算的透平机械叶片固有频率特性分析方法
CN111695193B (zh) 一种全局相关三维气动力数学模型的建模方法及系统
CN115841548B (zh) 一种叶片模型的计算机辅助生成方法及系统
EP2407902A2 (en) System and method for configuring a simulation model utilizing a tool for automatic input/output assignment
CN114372318A (zh) 设计参数确定方法、装置、电子设备及可读存储介质
CN112149262B (zh) 一种集束线及集束线三维设计模型构建方法
CN112231844A (zh) 零件模型的曲面离散方法、装置、设备及存储介质
CN112487575B (zh) 双转子系统高压转子及低压转子连接节点的检查方法
CN112287484B (zh) 一种基于矢量代理模型的复杂工程系统可靠性设计方法
CN113536709A (zh) 一种基于湍流建模机器学习的湍流模型修正方法
Reitenbach et al. Parametric datamodel for collaborative preliminary aircraft engine design
Henry et al. A fleet risk prediction methodology for mistuned ibrs using geometric mistuning models
CN116956619A (zh) 一种轻量级环境下的土木工程模型虚拟装配方法
CN114417508B (zh) 一种大气数据传感器仿真软件的验证方法及装置
US20060253269A9 (en) Method and apparatus for machine vector loop extraction
CN112329279B (zh) 坎贝尔图模态排序的方法
CN112861285B (zh) 压气机几何建模方法及可读存储介质
CN111177370B (zh) 一种自然语言处理的算法
CN107563041A (zh) 一种风电机组大部件静强度快速评估方法
Kim et al. Self-similarity in fractal and non-fractal networks
Liu et al. Study on burning surface regression algorithm under erosive burning based on CT images of solid rocket motor grain
CN116204974A (zh) 航空发动机叶片零件cad模型几何一致性评价方法
CN117390812B (zh) 一种cad图纸暖通风管结构化信息提取方法及系统
CN115048730B (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