CN116205107A - High-speed thin-wall revolution structure unbalance amount calculating method and related device - Google Patents
High-speed thin-wall revolution structure unbalance amount calculating method and related device Download PDFInfo
- Publication number
- CN116205107A CN116205107A CN202310180490.7A CN202310180490A CN116205107A CN 116205107 A CN116205107 A CN 116205107A CN 202310180490 A CN202310180490 A CN 202310180490A CN 116205107 A CN116205107 A CN 116205107A
- Authority
- CN
- China
- Prior art keywords
- assembly
- pose
- unbalance
- space pose
- model
- 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
- 238000000034 method Methods 0.000 title claims abstract description 47
- 230000008859 change Effects 0.000 claims abstract description 18
- 238000004458 analytical method Methods 0.000 claims abstract description 17
- 238000005259 measurement Methods 0.000 claims abstract description 16
- 230000008569 process Effects 0.000 claims abstract description 15
- 238000012546 transfer Methods 0.000 claims abstract description 14
- 238000012937 correction Methods 0.000 claims abstract description 11
- 238000004590 computer program Methods 0.000 claims description 20
- 238000009434 installation Methods 0.000 claims description 13
- 238000005520 cutting process Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 7
- 238000010586 diagram Methods 0.000 claims description 6
- 239000007787 solid Substances 0.000 claims description 6
- 238000000354 decomposition reaction Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 238000003860 storage Methods 0.000 claims description 4
- 230000013011 mating Effects 0.000 claims description 2
- 238000004088 simulation Methods 0.000 claims description 2
- 238000005457 optimization Methods 0.000 abstract description 3
- 238000013461 design Methods 0.000 abstract description 2
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000012876 topography Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 2
- 230000009191 jumping Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
Images
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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/04—Manufacturing
-
- 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)
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Geometry (AREA)
- Manufacturing & Machinery (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Economics (AREA)
- General Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Testing Of Balance (AREA)
Abstract
The invention discloses a high-speed thin-wall revolution structure unbalance amount calculating method and a related device, comprising the following steps: measuring and dividing the area of the single-stage rotary structure; establishing a mounting edge analysis model of the rotary structure by combining the measurement result and the region division result; establishing a geometric error transfer model of a steering engine rotating structure; the unbalance amount generated due to the rotation axis variation is calculated on the two correction surfaces after the assembly of each stage of rotary structure. According to the invention, the actually measured geometric shape data are fused with the finite element model, a high-speed thin-wall rotary structure spigot connection model with an actual surface shape is established, the influence of a rotor assembly process on a rotor rotation axis is considered, the influence of the change of the assembly axis on the initial unbalance amount and the mass center is considered, the accurate prediction of the unbalance amount of the rotor after assembly can be realized, the foundation is laid for the assembly optimization design of the high-speed thin-wall rotary structure, and the assembly one-time success rate and the assembly quality are improved.
Description
Technical Field
The invention belongs to the technical field of intelligent manufacturing, and relates to a high-speed thin-wall revolution structure unbalance amount calculating method and a related device.
Background
The high-speed thin-wall rotary structure is widely applied to important machinery such as aeroengines, gas turbines and the like, has typical high-speed, high-pressure and high-temperature operation characteristics, and has strict requirements on reliability and service life. The presence of unbalance will create large forces and moments during high speed rotation, severely affecting the vibration characteristics of the mechanical parts, which will also directly affect its lifetime and reliability. Therefore, the unbalance amount is an important index in manufacturing and assembling the high-speed thin-wall rotary structure.
The high-speed thin-wall rotary structure is generally assembled by a plurality of parts, and in the manufacturing process of the parts, no matter how precise machinery is used, the shape error and self-unbalance of the surfaces of the parts can be generated. Because the flange at the assembling seam allowance is small in thickness and has morphology errors, large deformation usually occurs in the primary screwing process of the bolt, and the deformation of the flange has large influence on the mass center position and the initial unbalance amount of the rotary structure. Therefore, the unbalance distribution condition of the rotary machine under different assembly processes is accurately predicted, and the method is particularly important to optimizing the overall unbalance of the rotor of the aeroengine and improving the service life and the reliability of the rotor.
The existing unbalance prediction method for the high-speed thin-wall rotary structure has some problems, mainly comprising the following steps: the influence of the assembly process on the unbalance of the multi-stage high-speed thin-wall rotary structure is not considered, and the initial unbalance change of the high-speed thin-wall rotary structure caused by the change of the central axis of the high-speed thin-wall rotary structure in the assembly process is not considered, so that the unbalance change trend prediction of the rotor is inaccurate, the rotor is inaccurate in optimization of the unbalance, and the assembly process is not provided with a clear guiding method.
Disclosure of Invention
The invention aims to solve the problems in the prior art and provides a high-speed thin-wall revolution structure unbalance amount calculating method and a related device.
In order to achieve the above purpose, the invention is realized by adopting the following technical scheme:
in a first aspect, the present invention provides a method for testing unbalance of a high-speed thin-wall rotating structure, comprising the steps of:
measuring and dividing the area of the single-stage rotary structure;
establishing a mounting edge analysis model of the rotary structure by combining the measurement result and the region division result;
establishing a geometric error transfer model of the multi-stage rotary structure by utilizing a homogeneous coordinate change theory, and calculating by combining with a mounting edge analysis model to obtain an assembled integral space pose;
according to the assembled integral space pose, calculating to obtain the space pose of the rotating structure after assembling of unbalanced mass and the space pose of the rotating structure after assembling of mass center, wherein the space pose is obtained after unbalanced decomposition;
and calculating the unbalance amount generated after the rotation axis is changed and the unbalance amount of each level of rotation structure on two correction surfaces after the assembly according to the space pose of the rotation structure after the assembly and the space pose of the centroid of the rotation structure after the assembly.
In a second aspect, the present invention provides a high-speed thin-wall revolution mechanic unbalance test system, comprising:
the measuring module is used for measuring and dividing the area of the single-stage rotary structure;
the first model building module is used for building a mounting edge analysis model of the rotary structure by combining the measurement result and the region division result;
the second model building module is used for building a geometric error transfer model of the steering engine rotating structure by utilizing the homogeneous coordinate change theory and calculating by combining with the installation edge analysis model to obtain an assembled integral space pose;
the first calculation module is used for calculating the space pose of the rotating structure after assembly and the space pose of the center of mass of the rotating structure after assembly according to the assembled whole space pose;
and the second calculation module is used for calculating the unbalance amount generated after the rotation axis is changed and the unbalance amount of each level of rotation structure on two correction surfaces after the assembly according to the space pose of the rotation structure after the assembly and the space pose of the centroid after the assembly of the unbalance mass.
In a third aspect, the present invention provides a computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, the processor implementing the steps of the method as described above when executing the computer program.
In a fourth aspect, the present invention provides a computer readable storage medium storing a computer program which when executed by a processor performs the steps of a method as described above.
Compared with the prior art, the invention has the following beneficial effects:
according to the invention, the actually measured geometric shape data are fused with the finite element model, a high-speed thin-wall rotary structure spigot connection model with an actual surface shape is established, the influence of a rotor assembly process on a rotor rotation axis is considered, the influence of the change of the assembly axis on the initial unbalance amount and the mass center is considered, the accurate prediction of the unbalance amount of the rotor after assembly can be realized, the foundation is laid for the assembly optimization design of the high-speed thin-wall rotary structure, and the assembly one-time success rate and the assembly quality are improved.
Drawings
For a clearer description of the technical solutions of the embodiments of the present invention, the drawings that are needed in the embodiments will be briefly described below, it being understood that the following drawings only illustrate some embodiments of the present invention and should not be considered as limiting the scope, and other related drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of the method of the present invention.
Fig. 2 is a flow chart of a method for calculating unbalance of a high-speed thin-wall rotary structure in consideration of an assembly process according to an embodiment of the present invention.
Fig. 3 is a rigid-flexible area division diagram of a single-stage slewing structure in an embodiment of the invention.
FIG. 4 is a schematic diagram of the center of mass offset caused by axis variation after assembly of a multi-stage rotor according to an embodiment of the present invention.
FIG. 5 is a schematic illustration of the initial imbalance variation caused by axis variation after assembly of a multi-stage rotor according to an embodiment of the present invention.
FIG. 6 is a schematic diagram of a change in centroid shift caused by axis change after assembly of a multi-stage rotor in accordance with an embodiment of the present invention.
Fig. 7 is a schematic diagram of the system of the present invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments of the present invention. The components of the embodiments of the present invention generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations.
Thus, the following detailed description of the embodiments of the invention, as presented in the figures, is not intended to limit the scope of the invention, as claimed, but is merely representative of selected embodiments of the invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
It should be noted that: like reference numerals and letters denote like items in the following figures, and thus once an item is defined in one figure, no further definition or explanation thereof is necessary in the following figures.
In the description of the embodiments of the present invention, it should be noted that, if the terms "upper," "lower," "horizontal," "inner," and the like indicate an azimuth or a positional relationship based on the azimuth or the positional relationship shown in the drawings, or the azimuth or the positional relationship in which the inventive product is conventionally put in use, it is merely for convenience of describing the present invention and simplifying the description, and does not indicate or imply that the apparatus or element to be referred to must have a specific azimuth, be configured and operated in a specific azimuth, and thus should not be construed as limiting the present invention. Furthermore, the terms "first," "second," and the like, are used merely to distinguish between descriptions and should not be construed as indicating or implying relative importance.
Furthermore, the term "horizontal" if present does not mean that the component is required to be absolutely horizontal, but may be slightly inclined. As "horizontal" merely means that its direction is more horizontal than "vertical", and does not mean that the structure must be perfectly horizontal, but may be slightly inclined.
In the description of the embodiments of the present invention, it should also be noted that, unless explicitly specified and limited otherwise, the terms "disposed," "mounted," "connected," and "connected" should be construed broadly, and may be, for example, fixedly connected, detachably connected, or integrally connected; can be mechanically or electrically connected; can be directly connected or indirectly connected through an intermediate medium, and can be communication between two elements. The specific meaning of the above terms in the present invention can be understood by those of ordinary skill in the art according to the specific circumstances.
The invention is described in further detail below with reference to the attached drawing figures:
referring to fig. 1, the embodiment of the invention discloses a method for testing unbalance of a high-speed thin-wall rotary structure, which comprises the following steps:
s1, measuring and dividing areas of a single-stage rotary structure;
s2, combining the measurement result and the region division result to establish a mounting edge analysis model of the rotary structure;
s3, establishing a geometric error transfer model of the steering engine rotating structure by using a homogeneous coordinate change theory, and calculating by combining with a mounting edge analysis model to obtain an assembled integral space pose;
s4, calculating to obtain the space pose of the rotating structure after assembly and the space pose of the centroid of the rotating structure after assembly according to the overall space pose after assembly, wherein the space pose of the rotating structure after assembly is obtained by decomposing the unbalance;
and S5, calculating the unbalance amount generated after the rotation axis is changed and the unbalance amount of each level of rotation structure after assembly on two correction surfaces according to the space pose of the rotation structure after assembly and the space pose of the centroid of the rotation structure after assembly.
Referring to fig. 2, the embodiment of the invention discloses a method for calculating unbalance of a high-speed thin-wall rotary structure by considering an assembly process, which comprises the following steps:
s101, measuring and obtaining part self data.
Measuring and obtaining real surface morphology error data of front and rear spigot of each part, and processing the data to obtain the spatial pose of the spigot position of each part; the initial unbalance of each component is measured by a horizontal hard support balancing machine.
S102, rigid-flexible area division of the single-stage slewing structure is carried out.
Stress and deformation mainly occur near the joint surface in the assembly process, so that the single-stage rotary structure can be cut into a rigid body part and a flexible body part, wherein the flexible body part is a region which is affected by the assembly deformation near the spigot, and the rigid part is a region which is basically not affected by the assembly;
s103, based on the soft body part of the single-stage rotary structure, combining with the real surface topography error data to establish a mounting edge analysis model of the rotary structure considering the assembly process.
Utilizing ANSYS APDL programming language to construct a solid model according to the high-speed thin-wall revolution structural part diagram: establishing a solid model of a flexible body part near the installation edge of each rotary structure in consideration of an assembly process, wherein the solid model comprises bolt contact, spigot end surface contact and spigot interference contact;
after grid division is carried out on the entity model of the soft body part, decoupling the entity model from the finite element model; and selecting nodes at corresponding positions on the grid by using ANSYS APDL programming language, and adding the real surface topography error data by driving the movement of coordinates of the nodes corresponding to the positions of the end surface and the spigot mating surface. In order to simplify the input of jumping data, the jumping error is expanded, so that the diameter jumps at the same angle on each section of the cylindrical surface are consistent, and the end jumps at the same angle on the end surface are consistent;
and (3) applying load and constraint to the finite element model, and constraining the soft body part of the single-stage rotary structure according to actual conditions. And determining a starting point and a tightening sequence of bolt tightening in the installation edge model, sequentially applying pre-tightening force to each bolt, writing in a load step, and simulating a real assembly condition through the load step. If the assembly batch is required, applying corresponding bolt pretightening force to each batch; and after corresponding contact, interference, constraint and bolt tightening are applied, simulation calculation is carried out on the installation edge model.
S104, calculating the space pose of the rigid body part through a similar triangle principle, and obtaining the space pose of the flexible body part by using a mounting edge analysis model
S105, establishing a multistage rotary structure geometric error transfer model considering an assembly process by using a homogeneous coordinate change theory.
Extracting coordinate data calculated by the installation edge model through ANSYA APDL command, and fitting the extracted coordinate points by using a least square method to obtain the space pose of the installation edge model;
calculating the spatial pose of the rigid body part and the cutting surface: as shown in fig. 3, the least square method fits the real topography data, and the spatial pose at the cutting surface is calculated through the eccentric position relationship and can be expressed as:
the end face at the cutting surface is not inclined because the cutting is performed parallel to the axis of the turntable, and A and B are all 0.Xc and Yc can be calculated by the eccentricity of the upper end face and the lower end face of the rotor.
And then calculating to obtain the space pose of the rigid body part, and calculating to obtain the assembled integral space pose according to the space pose of the flexible body part and the space pose of the rigid body part by a homogeneous coordinate change theory, so as to obtain the geometrical error transfer model of the multi-stage rotary structure.
S106 decomposes the initial unbalance amount.
According to the data measured by the horizontal hard support balancing machine, the initial unbalance amount is decomposed into unbalance mass and radius at the measuring surface, and the coordinates of the unbalance mass of each level of rotary structure relative to the front spigot error measuring surface can be obtained according to the position of the unbalance measuring surface and the decomposed radius.
S107, acquiring space positions of all levels of rotary structures.
Determining the space pose of each level of revolution structure after assembly through a geometric error transfer model; according to the space pose coordinates of all levels of parts after assembly, the space pose of the unbalanced mass after decomposition is obtained: the original decomposed unbalanced mass position is relative to the error measuring surface of the front spigot of the single-stage part, and when the space pose of each stage of parts is changed, the unbalanced mass position is also changed;
and obtaining the space pose of the decomposed unbalanced mass after the assembly of the rotary structure according to the space pose of the rotary structure at each stage and the coordinates of the unbalanced mass relative to the front spigot error measurement surface.
And according to the space pose of each level of rotary structure and the coordinates of the centroid position relative to the front spigot error measurement surface, obtaining the space pose of the assembled rotary structure of the centroid.
S108, calculating the unbalance amount of the multi-stage rotary structure considering the assembly process.
And according to the space positions of all levels of assembled parts of the multi-level rotary structure, taking the front end of the first part and the rear end of the last part, and calculating the assembled rotation axis according to the eccentric value.
Calculating the unbalance amount of each level of parts generated by the mass center off-axis: as shown in fig. 4, after the multi-stage rotating structure member is assembled, the rotation axis of the rotating structure is deviated due to the eccentricity, the centroid is deviated from the axis, the distance between the centroid and the rotation axis is calculated according to the centroid position coordinates, and the unbalance amount generated due to the deviation of the centroid position from the axis is calculated.
Calculating initial unbalance amount change of each level of parts caused by axis change: as shown in fig. 5, since the initial unbalance amount of each rotating structure member is measured under the axis of the rotating structure itself, the initial unbalance amount of the rotating structure is also changed by the axis variation after the assembly, the distance from the unbalanced mass obtained by decomposing the initial unbalance amount to the rotation axis is calculated from the unbalanced mass position obtained by decomposing the initial unbalance amount, and the unbalance amount due to the axis deviation is calculated.
Calculating unbalance amount of each level of parts: as shown in fig. 6, the unbalance amount due to the centroid shift caused by the axis shift and the initial unbalance amount due to the axis shift are vector-added and distributed to the two spigot error measurement surfaces, thereby obtaining the unbalance amount on each stage of parts. Vector addition is carried out on the unbalance amount on each level of parts after assembly at the same spigot error measurement surface to obtain the unbalance amount on each spigot error measurement surface, and then the unbalance amount on each spigot error measurement surface is projected to two correction surfaces of the revolution structure after assembly to obtain the unbalance amount on the two correction surfaces after assembly of each level of revolution structure.
As shown in fig. 7, an embodiment of the present invention discloses a high-speed thin-wall revolution structure unbalance amount test system, which includes:
the measuring module is used for measuring and dividing the area of the single-stage rotary structure;
the first model building module is used for building a mounting edge analysis model of the rotary structure by combining the measurement result and the region division result;
the second model building module is used for building a geometric error transfer model of the steering engine rotating structure by utilizing the homogeneous coordinate change theory and calculating by combining with the installation edge analysis model to obtain an assembled integral space pose;
the first calculation module is used for calculating the space pose of the rotating structure after assembly and the space pose of the center of mass of the rotating structure after assembly according to the assembled whole space pose;
and the second calculation module is used for calculating the unbalance amount generated after the rotation axis is changed and the unbalance amount of each level of rotation structure on two correction surfaces after the assembly according to the space pose of the rotation structure after the assembly and the space pose of the centroid after the assembly of the unbalance mass.
The embodiment of the invention provides computer equipment. The computer device of this embodiment includes: a processor, a memory, and a computer program stored in the memory and executable on the processor. The steps of the various method embodiments described above are implemented when the processor executes the computer program. Alternatively, the processor may implement the functions of the modules/units in the above-described device embodiments when executing the computer program.
The computer program may be divided into one or more modules/units, which are stored in the memory and executed by the processor to accomplish the present invention.
The computer equipment can be a desktop computer, a notebook computer, a palm computer, a cloud server and other computing equipment. The computer device may include, but is not limited to, a processor, a memory.
The processor may be a central processing unit (Central Processing Unit, CPU), but may also be other general purpose processors, digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), off-the-shelf programmable gate arrays (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, or the like.
The memory may be used to store the computer program and/or modules, and the processor may implement various functions of the computer device by running or executing the computer program and/or modules stored in the memory, and invoking data stored in the memory.
The modules/units integrated with the computer device may be stored in a computer readable storage medium if implemented in the form of software functional units and sold or used as stand alone products. Based on such understanding, the present invention may implement all or part of the flow of the method of the above embodiment, or may be implemented by a computer program to instruct related hardware, where the computer program may be stored in a computer readable storage medium, and when the computer program is executed by a processor, the computer program may implement the steps of each of the method embodiments described above. Wherein the computer program comprises computer program code which may be in source code form, object code form, executable file or some intermediate form etc. The computer readable medium may include: any entity or device capable of carrying the computer program code, a recording medium, a U disk, a removable hard disk, a magnetic disk, an optical disk, a computer Memory, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), an electrical carrier signal, a telecommunications signal, a software distribution medium, and so forth. It should be noted that the computer readable medium contains content that can be appropriately scaled according to the requirements of jurisdictions in which such content is subject to legislation and patent practice, such as in certain jurisdictions in which such content is subject to legislation and patent practice, the computer readable medium does not include electrical carrier signals and telecommunication signals.
The above is only a preferred embodiment of the present invention, and is not intended to limit the present invention, but various modifications and variations can be made to the present invention by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (9)
1. The method for testing the unbalance of the high-speed thin-wall rotary structure is characterized by comprising the following steps of:
measuring and dividing the area of the single-stage rotary structure;
establishing a mounting edge analysis model of the rotary structure by combining the measurement result and the region division result;
establishing a geometric error transfer model of the multi-stage rotary structure by utilizing a homogeneous coordinate change theory, and calculating by combining with a mounting edge analysis model to obtain an assembled integral space pose;
according to the assembled integral space pose, calculating to obtain the space pose of the rotating structure after assembling of unbalanced mass and the space pose of the rotating structure after assembling of mass center, wherein the space pose is obtained after unbalanced decomposition;
and calculating the unbalance amount generated after the rotation axis is changed and the unbalance amount of each level of rotation structure on two correction surfaces after the assembly according to the space pose of the rotation structure after the assembly and the space pose of the centroid of the rotation structure after the assembly.
2. The method for measuring unbalance of a high-speed thin-wall revolution structure according to claim 1, wherein the measuring and the area division of the single-stage revolution structure comprise:
the single-stage revolution structure is measured as follows:
measuring and obtaining real surface morphology error data of front and rear spigot of each part, and processing the data to obtain the spatial pose of the spigot position of each part; measuring by a horizontal hard support balancing machine to obtain initial unbalance of each component;
the single-stage rotary structure is divided into areas, and the method is as follows:
the single-stage rotary structure is cut into a rigid body part and a flexible body part, wherein the flexible body part is a region which is affected by assembly deformation near the spigot, and the rigid part is a region which is not affected by assembly basically.
3. The method for measuring unbalance of a high-speed thin-wall revolution structure according to claim 2, wherein the combining the measurement result and the area division result creates a mounting-edge analysis model of the revolution structure, comprising:
utilizing ANSYS APDL programming language to construct a solid model according to the high-speed thin-wall revolution structural part diagram: establishing a solid model of a flexible body part near the installation edge of each rotary structure in consideration of an assembly process, wherein the solid model comprises bolt contact, spigot end surface contact and spigot interference contact;
after grid division is carried out on the entity model of the soft body part, decoupling the entity model from the finite element model; selecting nodes at corresponding positions on the grid by using ANSYS APDL programming language, and adding real surface morphology error data by driving the movement of node coordinates corresponding to the positions of the end surface and the spigot mating surface; expanding the runout error to ensure that the diameter runout at the same angle on each section of the cylindrical surface is consistent, and the end runout at the same angle on the end surface is consistent;
load and constraint are applied to the finite element model, and the flexible body part of the rotary structure is constrained according to actual conditions; determining a starting point and a tightening sequence of bolt tightening in the installation edge model, sequentially applying pre-tightening force to each bolt, writing in a load step, and simulating a real assembly condition through the load step; if the assembly batch is required, applying corresponding bolt pretightening force to each batch; and after the corresponding contact, interference, constraint and bolt tightening are applied, performing simulation calculation on the installation edge model.
4. The method for testing unbalance of a high-speed thin-wall rotating structure according to claim 3, wherein the establishing a geometrical error transfer model of the multi-stage rotating structure by using the homogeneous coordinate variation theory and calculating by combining with the installation edge analysis model to obtain the assembled integral space pose comprises the following steps:
extracting coordinate data calculated by the installation edge model through ANSYA APDL command, and fitting the extracted coordinate points by using a least square method to obtain the space pose of the soft body part;
calculating the spatial pose of the rigid body part and the cutting surface: fitting real morphology data by using a least square method, and calculating to obtain the spatial pose of the cutting surface through an eccentric position relation:
wherein T represents the spatial pose of the mounting edge, (a, B, 1) is the normal vector obtained By ax+by+c=1 fitting; (X) c ，Y c ) The eccentricity of the measuring surface of the spigot of the rotor is calculated by the eccentricity of the upper spigot and the lower spigot, and h is the relative height of the measuring surface of the spigot; the end face at the cutting surface is not inclined because the cutting is performed parallel to the axis of the turntable, and both A and B are 0; the method comprises the steps of carrying out a first treatment on the surface of the X is X c And Y c The eccentric calculation of the upper end face and the lower end face of the rotor is adopted to obtain the rotor;
h is the total height of the rotor, H 1 And h 2 Is the cutting distance at both faces of the rotor, (eccX 1 ,eccY 1 ) And (eccX) 2 ,eccY 2 ) Respectively representing the eccentric coordinates of the lower spigot of the rotor and the eccentric coordinates of the upper spigot of the rotor; (X) c1 ，Y c1 ) And (X) c2 ，Y c2 ) Respectively representing eccentric coordinates of a lower cutting surface and an upper cutting surface of the rotor;
calculating according to the spatial pose of the cutting surface and the spatial pose of the upper seam allowance and the lower seam allowance to obtain the spatial pose of the rigid body part; and establishing a geometric error transfer model of the multi-stage rotary structure by using a homogeneous coordinate change theory, and calculating to obtain the assembled integral space pose according to the space pose of the rigid body part and the space pose of the soft body part.
5. The method for measuring unbalance of a high-speed thin-wall rotating structure according to claim 4, wherein the calculating the spatial pose of the unbalanced mass of the rotating structure after assembly and the spatial pose of the centroid of the rotating structure after assembly after the decomposition of the unbalance amount comprises:
decomposing the initial unbalance amount into unbalance mass and radius at the measuring surface according to the initial unbalance amount measured by the horizontal hard support balancing machine, and obtaining the coordinates of the unbalance mass of each level of rotary structure relative to the front spigot error measuring surface according to the position of the unbalance measuring surface and the decomposed radius;
determining the space pose of each level of revolution structure after assembly through a geometric error transfer model, and obtaining the space pose of the revolution structure after assembly of the decomposed unbalanced mass according to the space pose of each level of revolution structure and the coordinates of the unbalanced mass relative to the front spigot error measurement surface;
and determining the space pose of the assembled rotary structures at all levels through a geometric error transfer model, and obtaining the space pose of the centroid after the assembly according to the space pose of the assembled rotary structures at all levels and the coordinates of the centroid position relative to the error measurement surface of the front spigot.
6. The method for measuring unbalance of a high-speed thin-wall revolution structure according to claim 5, wherein the calculating of the unbalance amount due to the variation of the rotation axis and the unbalance amount of each revolution structure after the assembly in two correction surfaces comprises:
calculating according to the assembled integral space pose to obtain an assembled rotation axis, and according to the space pose of the rotation structure after assembly of the rotation axis and the unbalanced mass, obtaining the unbalanced quantity generated by the initial unbalanced quantity after the rotation axis is changed;
according to the space pose of the rotating structure after the rotation axis and the mass center are assembled, the unbalance amount generated by mass center offset caused by the rotation axis variation is obtained;
and calculating the unbalance amount of each level of rotary structure after assembly on two correction surfaces according to the unbalance amount generated after the rotation axis is changed and the unbalance amount generated by the mass center offset caused by the rotation axis is changed.
7. A high-speed thin-wall revolution mechanic unbalance amount test system, comprising:
the measuring module is used for measuring and dividing the area of the single-stage rotary structure;
the first model building module is used for building a mounting edge analysis model of the rotary structure by combining the measurement result and the region division result;
the second model building module is used for building a geometric error transfer model of the steering engine rotating structure by utilizing the homogeneous coordinate change theory and calculating by combining with the installation edge analysis model to obtain an assembled integral space pose;
the first calculation module is used for calculating the space pose of the rotating structure after assembly and the space pose of the center of mass of the rotating structure after assembly according to the assembled whole space pose;
and the second calculation module is used for calculating the unbalance amount generated after the rotation axis is changed and the unbalance amount of each level of rotation structure on two correction surfaces after the assembly according to the space pose of the rotation structure after the assembly and the space pose of the centroid after the assembly of the unbalance mass.
8. A computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that the processor implements the steps of the method according to any of claims 1-6 when the computer program is executed.
9. A computer readable storage medium storing a computer program, characterized in that the computer program when executed by a processor implements the steps of the method according to any of claims 1-6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310180490.7A CN116205107A (en) | 2023-02-28 | 2023-02-28 | High-speed thin-wall revolution structure unbalance amount calculating method and related device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310180490.7A CN116205107A (en) | 2023-02-28 | 2023-02-28 | High-speed thin-wall revolution structure unbalance amount calculating method and related device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116205107A true CN116205107A (en) | 2023-06-02 |
Family
ID=86514449
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310180490.7A Pending CN116205107A (en) | 2023-02-28 | 2023-02-28 | High-speed thin-wall revolution structure unbalance amount calculating method and related device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116205107A (en) |
-
2023
- 2023-02-28 CN CN202310180490.7A patent/CN116205107A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112699504B (en) | Assembly physical digital twin modeling method and device, electronic equipment and medium | |
| US7979233B2 (en) | Rotor assembly system and method | |
| US20050114059A1 (en) | Point cloud measuring system and method | |
| US7194387B1 (en) | Apparatus and method for determining the minimum zone for an array of features with a true position tolerance | |
| KR20220109434A (en) | Turbine assembly method, turbine assembly support program, and turbine assembly support device | |
| CN112985867B (en) | Steering engine testing method, device, equipment and storage medium | |
| JPH10267993A (en) | Fault analyzer | |
| CN116205107A (en) | High-speed thin-wall revolution structure unbalance amount calculating method and related device | |
| CN112434453B (en) | Bearing finite element model simplified equivalent method, system, medium, equipment and terminal | |
| CN111027194B (en) | Machine tool debugging method and device | |
| CN117171914A (en) | High-precision spindle performance test method and system | |
| CN116341038A (en) | Method, system, equipment and storage medium for machining crankshaft center hole | |
| CN116754134A (en) | Rotor imbalance state identification method based on test and simulation data fusion | |
| CN100582659C (en) | Complex position degree computing system and method | |
| CN114357653A (en) | End tooth connecting structure assembling and mounting angle optimization method considering tooth surface morphology | |
| CN114936483A (en) | Finite element equivalent modeling method for statics simulation analysis of large bearing | |
| CN114486033A (en) | Method, device and equipment for determining support axial force and storage medium | |
| CN113536482A (en) | Concentricity prediction method, system and device for high-speed rotating machinery | |
| CN111985122A (en) | Part tolerance optimization design method | |
| CN117824922B (en) | Dynamic balance performance test method and system for ultra-precise speed reducer | |
| CN114858123A (en) | Method and device for determining the concentricity of a toothed disc core relative to a rotation axis | |
| CN114763990B (en) | Method for measuring axial dimension of inner cavity of rotor and stator | |
| CN117561368A (en) | Method for estimating flange displacement of rotary machine, program for executing the method, and device for executing the method | |
| CN117990312A (en) | Method for testing joint rigidity of cooperative robot | |
| CN117830393A (en) | Method, system, equipment and medium for calculating space pose of high-voltage rotor |
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 |