CN106353033B - A kind of aero-engine centroid computing method - Google Patents

A kind of aero-engine centroid computing method Download PDF

Info

Publication number
CN106353033B
CN106353033B CN201611007353.XA CN201611007353A CN106353033B CN 106353033 B CN106353033 B CN 106353033B CN 201611007353 A CN201611007353 A CN 201611007353A CN 106353033 B CN106353033 B CN 106353033B
Authority
CN
China
Prior art keywords
coordinate system
relative
master
component
transformation matrix
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
CN201611007353.XA
Other languages
Chinese (zh)
Other versions
CN106353033A (en
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.)
Shenyang Jianzhu University
Original Assignee
Shenyang Jianzhu University
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 Shenyang Jianzhu University filed Critical Shenyang Jianzhu University
Priority to CN201611007353.XA priority Critical patent/CN106353033B/en
Publication of CN106353033A publication Critical patent/CN106353033A/en
Application granted granted Critical
Publication of CN106353033B publication Critical patent/CN106353033B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M1/00Testing static or dynamic balance of machines or structures
    • G01M1/12Static balancing; Determining position of centre of gravity
    • G01M1/122Determining position of centre of gravity

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Manufacture Of Motors, Generators (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The present invention discloses a kind of aero-engine centroid computing method;This method has many advantages, such as that computational accuracy is high, universality is strong, convenient and efficient;The technical scheme is that:A kind of aero-engine centroid computing method, mainly includes the following steps that:A. transformation matrix of the component local coordinate system relative to relative coordinate system is established;B. transformation matrix of the component local coordinate system relative to master coordinate system is established;C. Motor Mass Centre is sought;Purposes of the present invention:One, the enough barycenter parameters for accurately calculating engine complete machine of the present invention, facilitate lifting and the work such as the design of transport frock and the assembly of complete machine;Two, the present invention has very strong universality, can calculate the barycenter parameter of any one component;Three, the present invention equally has reference and reference for other rotation classes machinery.

Description

A kind of aero-engine centroid computing method
Technical field
The present invention relates to aero-engine field more particularly to a kind of computational methods of aero-engine barycenter.
Background technology
Aero-engine barycenter is an important feature parameter of engine, is that influence Engine Hoisting and transporting equipment are set One important technology index of the work such as meter.Design feature and installation of the calculating of Motor Mass Centre at present generally according to parts Angular position classify:Fan, Middle casing, compressor, combustion chamber, turbine etc. had into axial symmetry feature and rotating shaft The host part overlapped with engine axis is classified as one kind;Oil tank, combustion oil cooler, switching gear-box etc. had into non-axis Symmetrical feature and the unfixed attachment of angular position are classified as another kind of.It, can be by three-dimensionals such as Pro/E, UG for host part Modeling software directly reads its quality characteristics data, and computational accuracy is higher;But for attachment, then difficulty is larger, due to attachment number Amount is more, complicated, and it includes local coordinate system that the manufacturer of attachment, which generally only provides a simplified accessory model and one, Mass property file, therefore can not read by three-dimensional software merely.For problems, the measure taken at present is This class formation is carried out first to simplify processing, each attachment is regarded as a particle and is analyzed, and the selection of particle position, Randomness is stronger, and human error is larger.Although the quality and volume of most attachmentes are all smaller, since engine includes tens of A attachment, by accessory tapes Lai influence of the cumulative errors to engine entirety it is very important, still lack a kind of accurately hair at present Motivation barycenter calculation method of parameters.
Invention content
The purpose of the present invention is to propose to a kind of aero-engine centroid computing method, this method has computational accuracy high, general The advantages that adaptive is strong, convenient and efficient.
In order to realize the purpose of foregoing invention, the technical scheme is that:
A kind of aero-engine centroid computing method, is calculated based on homogeneous coordinate transformation method, includes mainly following Step:
B. transformation matrix of the component local coordinate system relative to relative coordinate system is established
A11. start owner's installation section to set, establish engine master coordinate system;
A12. according to part quality characteristic file, component local coordinate system is determined;
A13. local coordinate system origin is crossed to establish and the equidirectional relative coordinate system of each axis of master coordinate system;
A14. according to local coordinate system and relative coordinate system, change of the component local coordinate system relative to relative coordinate system is established Change matrix.
B. transformation matrix of the component local coordinate system relative to master coordinate system is established
B11. coordinate of the relative coordinate system origin in master coordinate system is sought;
B12. transformation matrix of the component relative coordinate system relative to master coordinate system is established;
B13. transformation matrix of the component local coordinate system relative to master coordinate system is established.
C. Motor Mass Centre is sought
C11. coordinate of the component barycenter in local coordinate system is sought;
C12. the transformation matrix according to component local coordinate system relative to master coordinate system seeks component barycenter in master coordinate system In coordinate;
C13. according to part quality and resultant force theorem of moments, Motor Mass Centre is sought.
It is an advantage of the invention that:
1, aero-engine centroid computing method of the invention can accurately calculate the barycenter parameter of engine complete machine, side Just lifting and the work such as the design of transport frock and the assembly of complete machine.
2, aero-engine centroid computing method of the invention is applicable not only to attachment, is also applied for host part, has very Strong universality can calculate the barycenter parameter of any one component.
3, aero-engine centroid computing method of the invention is applicable not only to aero-engine field, for other rotations The same meaning of class machinery.
Description of the drawings
Fig. 1 is the aero-engine centroid computing method schematic diagram of the present invention.
Fig. 2 is the aero-engine centroid computing method implementing procedure figure of the present invention.
Fig. 3 is the parts coordinate transform geometrical model of the aero-engine centroid computing method of the present invention.
Specific implementation mode
Illustrate present embodiment, its theoretical base of the aero-engine centroid computing method of present embodiment in conjunction with attached drawing 1-3 Plinth is homogeneous coordinate transformation method and resultant force theorem of moments, is calculated according to the following steps:
A. component local coordinate system σ is established2Relative to relative coordinate system σ3Transformation matrix
A11. start owner's installation section to set, establish engine master coordinate system σ1=[O1;x1,y1,z1];
A12. according to part quality characteristic file, component local coordinate system σ is determined2=[O2;x2,y2,z2];
A13. local coordinate system σ is crossed2Origin O2It establishes and master coordinate system σ1Each equidirectional relative coordinate system σ of axis3=[O3; x3,y3,z3], coordinate system σ2With σ3The angle of corresponding each axis is respectively α, β, γ;
A14. according to local coordinate system σ2With relative coordinate system σ3, establish component local coordinate system σ2Relative to relative coordinate It is σ3Transformation matrixWherein
B. component local coordinate system σ is established2Relative to master coordinate system σ1Transformation matrix
B11. relative coordinate system σ is sought3Origin is in master coordinate system σ1In coordinate
B12. component relative coordinate system σ is established3Relative to master coordinate system σ1Transformation matrix
B13. component local coordinate system σ is established2Relative to master coordinate system σ1Transformation matrix
C. Motor Mass Centre (X, Y, Z) is sought
C11. coordinate P of the component barycenter in local coordinate system is soughti=[x 'i y′i z′i 1]T, wherein i=1,2 ..., n;
C12. according to component local coordinate system σ2Relative to master coordinate system σ1Transformation matrixComponent barycenter is sought in master Coordinate in coordinate system(α, β, γ) in formula,(x 'i y′i z′i) three groups of parameters can easily be obtained by mass property file and three-dimensional software;
C13. last, according to part quality miWith resultant force theorem of moments, Motor Mass Centre (X, Y, Z) is sought.

Claims (4)

1. a kind of aero-engine centroid computing method, is calculated based on homogeneous coordinate transformation method, which is characterized in that
It mainly includes the following steps that:
A. transformation matrix of the component local coordinate system relative to relative coordinate system is established;
A11. start owner's installation section to set, establish engine master coordinate system;
A12. according to part quality characteristic file, component local coordinate system is determined;
A13. it is established and the equidirectional relative coordinate system of each axis of master coordinate system by component local coordinate system origin;
A14. according to component local coordinate system and relative coordinate system, change of the component local coordinate system relative to relative coordinate system is established Change matrix;
B. transformation matrix of the component local coordinate system relative to master coordinate system is established;
B11. coordinate of the relative coordinate system origin in master coordinate system is sought;
B12. transformation matrix of the component relative coordinate system relative to master coordinate system is established;
B13. transformation matrix of the component local coordinate system relative to master coordinate system is established;
C. Motor Mass Centre is sought;
C11. coordinate of the component barycenter in component local coordinate system is sought;
C12. the transformation matrix according to component local coordinate system relative to master coordinate system seeks component barycenter in master coordinate system Coordinate;
C13. according to part quality and resultant force theorem of moments, Motor Mass Centre is sought.
2. a kind of aero-engine centroid computing method according to claim 1, it is characterised in that:
A. component local coordinate system σ is established2Relative to relative coordinate system σ3Transformation matrix
A11. start owner's installation section to set, establish engine master coordinate system σ1=[O1;x1,y1,z1];
A12. according to part quality characteristic file, component local coordinate system σ is determined2=[O2;x2,y2,z2];
A13. pass through component local coordinate system σ2Origin O2It establishes and master coordinate system σ1Each equidirectional relative coordinate system σ of axis3= [O3;x3,y3,z3], coordinate system σ2With σ3The angle of corresponding each axis is respectively α, β, γ;
A14. according to component local coordinate system σ2With relative coordinate system σ3, establish component local coordinate system σ2Relative to relative coordinate It is σ3Transformation matrixWherein
3. a kind of aero-engine centroid computing method according to claim 1, it is characterised in that:
B. component local coordinate system σ is established2Relative to master coordinate system σ1Transformation matrix
B11. relative coordinate system σ is sought3Origin is in master coordinate system σ1In coordinate
B12. relative coordinate system σ is established3Relative to master coordinate system σ1Transformation matrix;
B13. component local coordinate system σ is established2Relative to master coordinate system σ1Transformation matrix
4. a kind of aero-engine centroid computing method according to claim 1, it is characterised in that:
C. Motor Mass Centre (X, Y, Z) is sought;
C11. coordinate P of the component barycenter in component local coordinate system is soughti=[x 'i y′i z′i 1]T, wherein i=1,2 ..., n;
C12. according to component local coordinate system σ2Relative to master coordinate system σ1Transformation matrixComponent barycenter is sought in principal coordinate Coordinate in system(α, β, γ) in formula,(x 'i y′i z′i ) Three groups of parameters can easily be obtained by part quality characteristic file and three-dimensional software;
C13. last, according to part quality mi and resultant force theorem of moments,Seek Motor Mass Centre (X, Y, Z)。
CN201611007353.XA 2016-11-15 2016-11-15 A kind of aero-engine centroid computing method Active CN106353033B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611007353.XA CN106353033B (en) 2016-11-15 2016-11-15 A kind of aero-engine centroid computing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611007353.XA CN106353033B (en) 2016-11-15 2016-11-15 A kind of aero-engine centroid computing method

Publications (2)

Publication Number Publication Date
CN106353033A CN106353033A (en) 2017-01-25
CN106353033B true CN106353033B (en) 2018-11-02

Family

ID=57861925

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611007353.XA Active CN106353033B (en) 2016-11-15 2016-11-15 A kind of aero-engine centroid computing method

Country Status (1)

Country Link
CN (1) CN106353033B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106979843B (en) * 2016-01-15 2019-08-27 中国航发常州兰翔机械有限责任公司 A kind of aero-engine center of gravity detection method
CN107238468B (en) * 2017-08-01 2019-06-18 中科新松有限公司 Discrimination method, device, equipment and the storage medium of mechanical arm tail end tool mass center
CN108267266B (en) * 2017-12-21 2019-11-22 北汽福田汽车股份有限公司 The measurement method and device of vehicle centroid
CN111307370A (en) * 2020-03-19 2020-06-19 青岛航空技术研究院(中国科学院工程热物理研究所青岛研究中心) Method for measuring rotational inertia of unmanned aerial vehicle

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105243172A (en) * 2015-07-22 2016-01-13 湖南大学 Engine suspension system optimization method capable of considering dynamic rigidity and damping
CN105913464A (en) * 2016-04-05 2016-08-31 中国科学院自动化研究所 Multi-body target online measurement method based on videos

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014124319A (en) * 2012-12-26 2014-07-07 Tokyo Univ Of Agriculture & Technology Ultrasonic calibration system, and ultrasonic calibration method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105243172A (en) * 2015-07-22 2016-01-13 湖南大学 Engine suspension system optimization method capable of considering dynamic rigidity and damping
CN105913464A (en) * 2016-04-05 2016-08-31 中国科学院自动化研究所 Multi-body target online measurement method based on videos

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于齐次坐标变换的制造误差建模研究;徐旭松等;《浙江大学学报(工学版)》;20080630;第42卷(第6期);全文 *
航空发动机质量、质心及转动惯量的一种计算方法及应用;喻双喜;《中国航空学会重量工程专业学术会》;19981231;全文 *
齐次坐标变换的测量机误差修正模型;陈宝刚等;《黑龙江科技学院学报》;20100331;第20卷(第2期);第143-145页第3-4节 *

Also Published As

Publication number Publication date
CN106353033A (en) 2017-01-25

Similar Documents

Publication Publication Date Title
CN106353033B (en) A kind of aero-engine centroid computing method
EP3314570B1 (en) Real-time, model-based object detection and pose estimation
CN106919944B (en) ORB algorithm-based large-view-angle image rapid identification method
CN103886593B (en) A kind of based on three-dimensional point cloud curved surface circular hole detection method
Korn et al. Color supported generalized-ICP
US9868212B1 (en) Methods and apparatus for determining the pose of an object based on point cloud data
US11833692B2 (en) Method and device for controlling arm of robot
CN110879400A (en) Method, equipment and storage medium for fusion positioning of laser radar and IMU
JP5706576B2 (en) Offset estimation apparatus, offset estimation method, offset estimation program, and information processing apparatus
CN109465822A (en) Based on 3D vision hand and eye calibrating method
CN106991705B (en) Position parameter estimation method based on P3P algorithm
JP2013205175A (en) Device, method and program for recognizing three-dimensional target surface
CN110595479B (en) SLAM track evaluation method based on ICP algorithm
Dan et al. An algorithm of NURBS surface fitting for reverse engineering
CN107063237A (en) A kind of method and apparatus for measuring gestures of object angle
CN105574884B (en) A kind of medical robot DH parameter calibration methods
CN108108555B (en) Assembly size chain calculation method, device and equipment
CN107066726A (en) A kind of Digit Control Machine Tool rotary shaft error of perpendicularity modeling method
CN109591051A (en) The compensation method of composite machine people's homework precision, system and storage medium
CN104573144A (en) System and method for simulating offline point cloud of measuring equipment
CN112591132A (en) Method for transforming rotational inertia coordinates of control surface of airplane
JP5462662B2 (en) Position / orientation measurement apparatus, object identification apparatus, position / orientation measurement method, and program
KR101604319B1 (en) Geometrically exact isogeometric shape sensitivity analysis method in curvilinear coordinate system of shell structure
Wang et al. Positionit: an image-based remote target localization system on smartphones
CN106373183A (en) Model establishing method based on three-dimensional data processing

Legal Events

Date Code Title Description
C06 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