CN110245425B - Air conditioner compressor excitation identification method and computer device - Google Patents
Air conditioner compressor excitation identification method and computer device Download PDFInfo
- Publication number
- CN110245425B CN110245425B CN201910521858.5A CN201910521858A CN110245425B CN 110245425 B CN110245425 B CN 110245425B CN 201910521858 A CN201910521858 A CN 201910521858A CN 110245425 B CN110245425 B CN 110245425B
- Authority
- CN
- China
- Prior art keywords
- matrix
- compressor
- vibration acceleration
- obtaining
- hammering
- 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
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
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Control Of Positive-Displacement Pumps (AREA)
Abstract
The invention provides an air conditioner compressor excitation identification method and a computer device, wherein the method comprises the following steps: when the compressor runs at a preset frequency, acquiring vibration acceleration matrixes of all acceleration sensors on the surface of the compressor; obtaining a mass matrix of the compressor according to the centroid coordinate and the inertia parameter matrix; acquiring a rigidity matrix of a rubber foot pad system of the compressor; obtaining a mass center vibration acceleration matrix of the compressor according to the vibration acceleration matrixes of all the acceleration sensors; and obtaining the mass center equivalent load of the compressor according to the mass matrix, the rigidity matrix and the mass center vibration acceleration matrix. The computer device comprises a controller, and the controller realizes the excitation identification method of the air conditioner compressor when executing the computer program stored in the memory. The invention can solve the problem that the load of the compressor cannot be equivalent to the mass center.
Description
Technical Field
The invention relates to the technical field of air conditioner compressors, in particular to an air conditioner compressor excitation identification method and a computer device applying the air conditioner compressor excitation identification method.
Background
At present, the exciting force of the compressor is usually calculated by a compressor dynamic model, and mainly comprises a pneumatic torque load, an electromagnetic torque load of a motor, an eccentric load of a rotor in dynamic balance and the like. The main disadvantage of this method is the need to obtain accurate values for various parameters, some of which are difficult to obtain, such as gas torque load, and in some cases the compressor is supplied by a third party, and no internal models and parameters are available. In the early stage of development of a compressor piping system, the compressor piping system needs to be subjected to simulation analysis to predict vibration and stress of a pipeline, but the excitation of the compressor needs to be known. The existing excitation identification cannot be equivalent to the centroid of the compressor, so that the existing excitation identification is not beneficial to being close to the reality and is not beneficial to the simulation analysis of a compressor piping system.
Disclosure of Invention
The first purpose of the invention is to provide an air conditioner compressor excitation identification method for solving the problem that the compressor load cannot be equivalent to the mass center.
A second object of the invention is to provide a computer arrangement that solves the problem of the compressor load not being equivalent to the centroid.
In order to achieve the first object, the present invention provides an excitation recognition method for an air conditioner compressor, including: when the compressor runs at a preset frequency, acquiring vibration acceleration matrixes of all acceleration sensors on the surface of the compressor; obtaining test result data of hammering a single compressor body, obtaining a mass center coordinate and an inertia parameter matrix of the compressor according to the test result data, and obtaining a mass matrix of the compressor according to the mass center coordinate and the inertia parameter matrix; establishing a simplified model of the compressor, and acquiring a rigidity matrix of a rubber foot pad system of the compressor; obtaining a mass center vibration acceleration matrix of the compressor according to the vibration acceleration matrixes of all the acceleration sensors; and obtaining the mass center equivalent load of the compressor according to the mass matrix, the rigidity matrix and the mass center vibration acceleration matrix.
According to the scheme, the vibration acceleration of the surface of the compressor is obtained under the preset frequency, the compressor single body is hammered, the frequency response function of the compressor single body is obtained through the acceleration sensor, the inertial parameters of the compressor are obtained through matrix calculation, and then the rigidity matrix and the mass center vibration acceleration of the foot pad of the compressor are obtained through finite element simulation, so that the equivalent load of the mass center of the compressor is obtained. Therefore, the problem that the load of the compressor cannot be equivalent to the mass center is solved through the excitation identification method of the air conditioner compressor, and the finite element modeling of the piping system of the compressor is closer to the reality.
In a further aspect, the step of obtaining a vibration acceleration matrix of all acceleration sensors on the surface of the compressor includes: and acquiring vibration acceleration time domain data collected by all the acceleration sensors, and performing frequency domain conversion on the vibration acceleration time domain data to obtain a vibration acceleration matrix.
Therefore, the vibration acceleration time domain data is subjected to frequency domain conversion, so that the vibration acceleration can be conveniently analyzed, and the detection efficiency is improved.
In a further aspect, the step of obtaining a vibration acceleration matrix of all acceleration sensors on the surface of the compressor further includes: and judging whether the detection position needs to be replaced or not according to the vibration acceleration time domain data corresponding to each acceleration sensor, and if so, generating prompt information for replacing the detection position.
Therefore, in order to ensure the reliability of the data acquired by the acceleration sensor, the vibration acceleration time domain data needs to be judged, so that whether position acquisition needs to be replaced or not is confirmed, and the acceleration sensor is ensured to be at the optimal acquisition position.
In a further scheme, the step of obtaining test result data of hammering a single compressor body and obtaining a mass center coordinate and an inertia parameter matrix of the compressor according to the test result data comprises the following steps: when hammering is carried out on a compressor by preset points, a hammering vibration acceleration matrix and a hammering excitation matrix of all acceleration sensors on the surface of the compressor are obtained during each hammering; obtaining a hammering vibration acceleration matrix of a coordinate origin of a reference coordinate system of the compressor three-dimensional model according to the hammering vibration acceleration matrix and the hammering excitation matrix; obtaining a mass center coordinate of the compressor according to the hammering vibration acceleration matrix of the coordinate origin and the coordinate of the coordinate origin; and obtaining an inertia parameter matrix according to the mass, the mass center coordinate, the hammering excitation matrix, the hammering vibration acceleration matrix and the hammering point coordinate of the compressor.
Therefore, when the mass center coordinate and the inertia parameter matrix are carried out, the hammering vibration acceleration matrix of the coordinate origin is calculated by acquiring the hammering vibration acceleration matrix and the hammering excitation matrix of the acceleration sensor during hammering, and the mass center coordinate of the compressor is acquired through the hammering vibration acceleration matrix of the coordinate origin, so that the inertia parameter matrix is acquired.
In a further aspect, the step of obtaining a stiffness matrix of the compressor rubber foot pad system comprises: acquiring a single stiffness matrix of a single rubber foot pad in the three main shaft directions; constructing a transformation matrix of intersection point coordinates of the centroid coordinates and the single rubber foot pad in three main shaft directions; and obtaining the rigidity matrix of the compressor rubber foot pad system according to the monomer rigidity matrix and the conversion matrix.
According to the scheme, in order to better simulate the actual working environment of the compressor, the rubber foot pad system of the compressor is added to a simplified model for building the compressor. When the rigidity matrix of the compressor rubber foot pad system is obtained, the actual compressor rubber foot pad system can be simulated by obtaining the rigidity matrix of the compressor rubber foot pad system through the monomer rigidity matrix and the conversion matrix of the single rubber foot pad.
In a further aspect, the step of obtaining the centroid vibration acceleration matrix of the compressor according to the vibration acceleration matrices of all the acceleration sensors includes: constructing a transfer matrix of the vibration acceleration of each acceleration sensor and the vibration acceleration of the mass center of the compressor when the compressor runs; and obtaining a mass center vibration acceleration matrix according to the vibration acceleration matrixes of all the acceleration sensors and the transfer matrix.
Therefore, when the mass center vibration acceleration matrix is obtained, the mass center vibration acceleration matrix can be conveniently obtained by constructing a transfer matrix of the vibration acceleration of the acceleration sensor and the vibration acceleration of the mass center of the compressor.
In a further aspect, the centroid equivalent load is obtained by the following equation:
wherein f is a preset frequency, K is a stiffness matrix, M is a mass matrix,
is a centroid vibration acceleration matrix.
Therefore, the working condition of the actual mass center of the compressor can be more accurately obtained through the formula of the equivalent load of the mass center.
In a further scheme, after the step of obtaining the centroid equivalent load of the compressor according to the mass matrix, the rigidity matrix and the centroid vibration acceleration matrix, the method further comprises the following steps: and loading the centroid equivalent load to a finite element model of the compressor, and correcting the centroid equivalent load.
Therefore, after the centroid equivalent load is obtained, the centroid equivalent load is loaded to a finite element model of the compressor for verification, and the centroid equivalent load is corrected, so that the centroid equivalent load is more accurate.
In a further aspect, the step of correcting the centroid equivalent load comprises: acquiring a corrected mass center vibration acceleration matrix corresponding to a finite element model of the compressor under a preset frequency; and correcting the centroid equivalent load according to the centroid vibration acceleration matrix and the corrected centroid vibration acceleration matrix.
Therefore, the centroid equivalent load is loaded to the finite element model of the compressor, the corrected centroid vibration acceleration matrix is obtained, the centroid equivalent load is corrected through the corrected centroid vibration acceleration matrix and the corrected centroid vibration acceleration, and the accuracy of the centroid equivalent load can be improved.
In order to achieve the second objective of the present invention, the present invention provides a computer device including a processor and a memory, wherein the memory stores a computer program, and the computer program is executed by the processor to implement the steps of the above-mentioned method for identifying the excitation of the air conditioner compressor.
Drawings
FIG. 1 is a flowchart illustrating an embodiment of an excitation recognition method for an air conditioner compressor according to the present invention.
Fig. 2 is a schematic view of an installation position of an acceleration sensor in an embodiment of an excitation identification method for an air conditioner compressor according to the present invention.
Fig. 3 is a schematic view of another view of the installation position of the acceleration sensor in the embodiment of the method for identifying the excitation of the air conditioner compressor according to the invention.
FIG. 4 is a schematic diagram of a simplified model of a compressor in an embodiment of the method for identifying excitation of an air conditioner compressor according to the present invention.
The invention is further explained with reference to the drawings and the embodiments.
Detailed Description
The invention relates to an excitation identification method of an air conditioner compressor, which is a computer program applied to excitation identification equipment of the air conditioner compressor and is used for realizing the excitation identification of the air conditioner compressor. The invention also provides a computer device which comprises a controller, wherein the controller is used for realizing the steps of the excitation identification method of the air conditioner compressor when executing the computer program stored in the memory. The present invention also provides a computer readable storage medium having a computer program stored thereon, the computer program, when executed by a controller, implementing the steps of the above-described excitation identification method for an air conditioner compressor.
The embodiment of the method for identifying the excitation of the air conditioner compressor comprises the following steps:
as shown in fig. 1, when performing the excitation recognition of the air conditioner compressor, the method of the present invention first performs step S1, and obtains the vibration acceleration matrix of all acceleration sensors on the surface of the compressor when the compressor runs at the preset frequency. In this embodiment, the preset frequency may be set according to the working frequency corresponding to the compressor during the nominal cooling, the nominal heating, the overload cooling, and the overload heating of the air conditioner external unit according to the vibration acceleration of the air conditioner external unit under the four working conditions of the nominal cooling, the nominal heating, the overload cooling, and the overload heating.
In this embodiment, the step of obtaining the vibration acceleration matrix of all the acceleration sensors on the surface of the compressor includes: and acquiring vibration acceleration time domain data collected by all the acceleration sensors, and performing frequency domain conversion on the vibration acceleration time domain data to obtain a vibration acceleration matrix.
Before acquiring the vibration acceleration matrix of all the acceleration sensors on the surface of the compressor, the positions and the number of the acceleration sensors can be set as required, and preferably, the number of the acceleration sensors is more than six. Referring to fig. 2 and 3, the surface of the compressor 1 is provided with eight acceleration sensors 2, wherein four acceleration sensors 2 are located on the side near the top of the compressor 1 and are arranged around the circumference of the compressor 1, and the other four acceleration sensors 2 are located on the side near the bottom of the compressor 1 and are arranged around the circumference of the compressor 1.
When the vibration acceleration matrix of all acceleration sensors on the surface of the compressor is obtained, vibration acceleration time domain data x corresponding to each acceleration sensor is obtained by obtaining vibration acceleration time domain data collected by all acceleration sensorsi(t), wherein i is the ith sensor, i is 1,2,3.. n, and n is the total number of sensors. For vibration acceleration time domain data xiAnd (t) performing windowing and Fourier transform, and performing frequency domain conversion on the vibration acceleration time domain data to obtain vibration acceleration frequency domain data. Each vibration acceleration frequency domain data can be divided into components of three-dimensional coordinate axes:
recording all vibration acceleration frequency domain data as an acceleration matrix:
in this embodiment, the step of obtaining the vibration acceleration matrix of all the acceleration sensors on the surface of the compressor further includes: and judging whether the detection position needs to be replaced or not according to the vibration acceleration time domain data corresponding to each acceleration sensor, and if so, generating prompt information for replacing the detection position. When the vibration acceleration time domain number is judged, data collection is carried out after the position of the acceleration sensor needs to be replaced when any one of the conditions is judged to be met, and therefore reliability of the data is guaranteed. And prompting the user to replace the detection position in time by generating prompt information for replacing the detection position.
After the vibration acceleration matrix is obtained, step S2 is executed to obtain test result data of hammering the compressor single body, obtain a centroid coordinate and an inertia parameter matrix of the compressor according to the test result data, and obtain a mass matrix of the compressor according to the centroid coordinate and the inertia parameter matrix. When the compressor monomer is subjected to hammering test, a special tool for installing the compressor is required, and the special tool is well known to those skilled in the art and is not described herein again.
In this embodiment, the step of obtaining test result data of hammering the compressor monomer and obtaining the barycenter coordinate and the inertia parameter matrix of the compressor according to the test result data includes: when hammering is carried out on a compressor by preset points, a hammering vibration acceleration matrix and a hammering excitation matrix of all acceleration sensors on the surface of the compressor are obtained during each hammering; obtaining a hammering vibration acceleration matrix of a coordinate origin of a reference coordinate system of the compressor three-dimensional model according to the hammering vibration acceleration matrix and the hammering excitation matrix; obtaining a mass center coordinate of the compressor according to the hammering vibration acceleration matrix of the coordinate origin and the coordinate of the coordinate origin; and obtaining an inertia parameter matrix according to the mass, the mass center coordinate, the hammering excitation matrix, the hammering vibration acceleration matrix and the hammering point coordinate of the compressor.
When the compressor is subjected to preset point hammering, the preset point can be set as required, and is preferably larger than six. And when a plurality of points are taken on the surface of the compressor and the compressor is unidirectionally excited by adopting a hammering method, the direction of each hammering is consistent with the direction of a reference coordinate system of the three-dimensional model. The excitation direction comprises X, Y, Z three directions, and the hammering vibration acceleration matrix of all acceleration sensors on the surface of the compressor is obtained by:
wherein the content of the first and second substances,
and the frequency domain data of the hammering vibration acceleration are components of three-dimensional coordinate axes. And processing the force acquired by the force resistance head when the compressor is hammered by the force hammer each time to obtain a hammering excitation matrix, and recording the hammering excitation matrix as:
wherein, Fqx,Fqy,FqzRespectively, the force components in the direction X, Y, Z when the hammer blows. And after obtaining the hammering vibration acceleration matrix and the hammering excitation matrix, obtaining the hammering vibration acceleration matrix of the coordinate origin of the reference coordinate system of the compressor three-dimensional model according to the hammering vibration acceleration matrix and the hammering excitation matrix. When acquiring the hammering vibration acceleration matrix, according to the equation system
(wherein, Xi,Yi,ZiIs the coordinate position of the acceleration sensor, and alpha, beta and gamma are the angular displacement of the mass center), and the hammering vibration acceleration matrix of the coordinate origin of the reference coordinate system of the three-dimensional model is obtained through conversion
And after the hammering vibration acceleration matrix of the coordinate origin is obtained, obtaining the mass center coordinate of the compressor according to the hammering vibration acceleration matrix of the coordinate origin and the coordinate of the coordinate origin. When a mass center coordinate is obtained, a relational expression of an excitation matrix, the mass center matrix and a coordinate origin is hammered simultaneously
Obtaining, wherein xc,yc,zcIs the coordinate of the mass center,
being the centre of massAcceleration of vibration. And after the mass center coordinate is obtained, obtaining an inertia parameter matrix according to the mass of the compressor, the mass center coordinate, the hammering excitation matrix, the hammering vibration acceleration matrix and the hammering point coordinate. When obtaining the inertia parameter matrix, the inertia parameter matrix can be obtained by simultaneous torque matrix
Where m is the mass of the compressor, xq,yq,zqObtaining an inertia parameter matrix I ═ I of the compressor by solving for the position coordinates of the excitation pointsxxIyyIzzIxyIyzIxz]. After obtaining the inertial parameter matrix, a mass matrix for the compressor may be established
Wherein, I3×3Is an identity matrix, thereby obtaining a mass matrix of the compressor.
And after obtaining the mass matrix of the compressor, executing the step S3, establishing a simplified model of the compressor, and obtaining the rigidity matrix of the rubber foot pad system of the compressor. Wherein, the step of obtaining the rigidity matrix of the compressor rubber foot pad system comprises the following steps: acquiring a single stiffness matrix of a single rubber foot pad in the three main shaft directions; constructing a transformation matrix of intersection point coordinates of the centroid coordinates and the single rubber foot pad in three main shaft directions; and obtaining the rigidity matrix of the compressor rubber foot pad system according to the monomer rigidity matrix and the conversion matrix.
When a simplified model of the compressor is established, as shown in fig. 4, the compressor 1 is simplified into a mass point, the rubber foot pad is simplified into a spring 3, the spring 3 is rigidly connected with the ground, and the mass point is rigidly connected with the spring 3. The stress of each rubber foot pad can be solved through a mechanical equation, the stress is used as a preload and is applied to the rubber foot pads, the finite element modeling mode of the foot pads is that the foot pads and the machine feet are in contact through different types, friction is set, the contact is close to reality, and therefore the rigidity matrix of the single rubber foot pad in the three main shaft directions is obtained through solving
This is a technique well known to those skilled in the art and will not be described in detail herein.
After a single body rigidity matrix of a single rubber foot pad in three main shaft directions is obtained, a transformation matrix of a centroid coordinate and an intersection point coordinate of the single rubber foot pad in the three main shaft directions is constructed
Then, obtaining a rigidity matrix of the compressor rubber foot pad system according to the single rigidity matrix and the conversion matrix, wherein the rigidity matrix of the compressor rubber foot pad system is obtained by the following formula:
wherein T ═ T1… Ti]。
After the rigidity matrix of the rubber foot pad system of the compressor is obtained, step S4 is executed, and the centroid vibration acceleration matrix of the compressor is obtained according to the vibration acceleration matrices of all the acceleration sensors. When the mass center vibration acceleration matrix of the compressor is obtained according to the vibration acceleration matrixes of all the acceleration sensors, a transfer matrix of the vibration acceleration of each acceleration sensor and the vibration acceleration of the mass center of the compressor is constructed when the compressor runs
Obtaining a mass center vibration acceleration matrix according to the vibration acceleration matrixes of all the acceleration sensors and the transfer matrix, wherein the mass center vibration acceleration matrix is obtained by the following formula:
wherein E ═ E1… Ei]。
And after the mass center vibration acceleration matrix is obtained, executing the step S5, and obtaining the mass center equivalent load of the compressor according to the mass matrix, the rigidity matrix and the mass center vibration acceleration matrix. The centroid equivalent load is obtained by the following formula:
wherein f is a preset frequency, K is a stiffness matrix, M is a mass matrix,
is a centroid vibration acceleration matrix.
And after the centroid equivalent load is obtained, executing step S6, loading the centroid equivalent load to the finite element model of the compressor, and correcting the centroid equivalent load. The step of correcting the centroid equivalent load comprises the following steps: acquiring a corrected mass center vibration acceleration matrix corresponding to a finite element model of the compressor under a preset frequency; and correcting the centroid equivalent load according to the centroid vibration acceleration matrix and the corrected centroid vibration acceleration matrix.
Loading the obtained centroid equivalent load F (f) to a finite element model of the compressor, analyzing and solving natural frequency harmonic response of the finite element model, and obtaining a corresponding corrected centroid vibration acceleration matrix of the finite element model of the compressor under a preset frequency
Comparing the centroid vibration acceleration matrix with the corrected centroid vibration acceleration matrix to obtain an acceleration error value
According to the error value
And (4) solving the load difference value, and accumulating the load difference value on the basis of the original centroid equivalent load so as to obtain the corrected centroid equivalent load. The obtained corrected centroid equivalent load can be applied to a finite element model of a compressor pipeline system, and responses such as vibration, stress and the like of a pipeline can be further solved.
The embodiment of the computer device comprises:
the computer device of the embodiment comprises a controller, and the steps of the method for identifying the excitation of the air conditioner compressor are realized when the controller executes a computer program.
For example, a computer program may be partitioned into one or more modules, which are stored in a memory and executed by a controller to implement the present invention. One or more of the modules may be a sequence of computer program instruction segments for describing the execution of a computer program in a computer device that is capable of performing certain functions.
The computer device may include, but is not limited to, a controller, a memory. Those skilled in the art will appreciate that the computer apparatus may include more or fewer components, or combine certain components, or different components, e.g., the computer apparatus may also include input-output devices, network access devices, buses, etc.
For example, the controller may be a Central Processing Unit (CPU), other general purpose controller, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), an off-the-shelf programmable Gate Array (FPGA) or other programmable logic device, discrete Gate or transistor logic, discrete hardware components, and so on. The general controller may be a microcontroller or the controller may be any conventional controller or the like. The controller is the control center of the computer device and connects the various parts of the entire computer device using various interfaces and lines.
The memory may be used to store computer programs and/or modules, and the controller may implement various functions of the computer apparatus by executing or otherwise executing the computer programs and/or modules stored in the memory and invoking data stored in the memory. For example, the memory may mainly include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required for at least one function (e.g., a sound receiving function, a sound-to-text function, etc.), and the like; the storage data area may store data (e.g., audio data, text data, etc.) created according to the use of the cellular phone, etc. In addition, the memory may include high speed random access memory, and may also include non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a flash memory Card (FlashCard), at least one magnetic disk storage device, a flash memory device, or other volatile solid state storage device.
Computer-readable storage medium embodiments:
the modules integrated by the computer apparatus of the above embodiments, if implemented in the form of software functional units and sold or used as independent products, may be stored in a computer-readable storage medium. Based on such understanding, all or part of the processes of the embodiments of the method for identifying excitation of an air conditioner compressor may also be implemented by a computer program instructing related hardware, where the computer program may be stored in a computer readable storage medium, and when the computer program is executed by a controller, the steps of the embodiments of the method for identifying excitation of an air conditioner compressor may be implemented. Wherein the computer program comprises computer program code, which may be in the form of source code, object code, an executable file or some intermediate form, etc. The storage medium may include: any entity or device capable of carrying computer program code, recording medium, U.S. disk, removable hard disk, magnetic disk, optical disk, computer Memory, Read-Only Memory (ROM), Random Access Memory (RAM), electrical carrier wave signals, telecommunications signals, software distribution media, and the like. It should be noted that the computer readable medium may contain other components which may be suitably increased or decreased as required by legislation and patent practice in jurisdictions, for example, in some jurisdictions, in accordance with legislation and patent practice, the computer readable medium does not include electrical carrier signals and telecommunications signals.
According to the method for identifying the excitation of the air conditioner compressor, the vibration acceleration of the surface of the compressor is obtained under the preset frequency, the compressor single body is hammered, the frequency response function of the compressor single body is obtained through the acceleration sensor, the inertial parameters of the compressor are obtained through matrix calculation, and then the rigidity matrix and the mass center vibration acceleration of the foot pad of the compressor are obtained through finite element simulation, so that the equivalent load of the mass center of the compressor is obtained. Therefore, the problem that the load of the compressor cannot be equivalent to the mass center is solved through the excitation identification method of the air conditioner compressor, and the finite element modeling of the piping system of the compressor is closer to the reality.
It should be noted that the above is only a preferred embodiment of the present invention, but the design concept of the present invention is not limited thereto, and any insubstantial modifications made by using the design concept also fall within the protection scope of the present invention.
Claims (9)
1. An air conditioner compressor excitation identification method is characterized in that: the method comprises the following steps:
when a compressor runs at a preset frequency, acquiring vibration acceleration matrixes of all acceleration sensors on the surface of the compressor;
obtaining test result data of hammering a single compressor, obtaining a mass center coordinate and an inertia parameter matrix of the compressor according to the test result data, and obtaining a mass matrix of the compressor according to the mass center coordinate and the inertia parameter matrix;
establishing a simplified model of the compressor, and acquiring a rigidity matrix of a rubber foot pad system of the compressor;
obtaining a mass center vibration acceleration matrix of the compressor according to the vibration acceleration matrixes of all the acceleration sensors;
obtaining a mass center equivalent load of the compressor according to the mass matrix, the rigidity matrix and the mass center vibration acceleration matrix:
wherein f is a preset frequency, K is the stiffness matrix, M is the mass matrix,
and the mass center vibration acceleration matrix is obtained.
2. The air conditioner compressor excitation identification method according to claim 1, wherein:
the step of obtaining the vibration acceleration matrix of all the acceleration sensors on the surface of the compressor comprises the following steps:
and acquiring vibration acceleration time domain data collected by all the acceleration sensors, and performing frequency domain conversion on the vibration acceleration time domain data to obtain the vibration acceleration matrix.
3. The air conditioner compressor excitation identification method according to claim 2, wherein:
the step of obtaining the vibration acceleration matrix of all the acceleration sensors on the surface of the compressor further comprises:
and judging whether the detection position needs to be replaced or not according to the vibration acceleration time domain data corresponding to each acceleration sensor, and if so, generating prompt information of replacing the detection position.
4. The air conditioner compressor excitation identification method according to any one of claims 1 to 3, wherein:
the step of obtaining test result data of hammering a single compressor and obtaining the barycenter coordinate and the inertia parameter matrix of the compressor according to the test result data comprises the following steps:
when hammering is carried out on the compressor by preset points, acquiring a hammering vibration acceleration matrix and a hammering excitation matrix of all acceleration sensors on the surface of the compressor during each hammering;
obtaining a hammering vibration acceleration matrix of a coordinate origin of a reference coordinate system of a compressor three-dimensional model according to the hammering vibration acceleration matrix and the hammering excitation matrix;
obtaining a mass center coordinate of the compressor according to the hammering vibration acceleration matrix of the coordinate origin and the coordinate of the coordinate origin;
and obtaining the inertia parameter matrix according to the mass of the compressor, the mass center coordinate, the hammering excitation matrix, the hammering vibration acceleration matrix and the hammering point coordinate.
5. The air conditioner compressor excitation identification method according to any one of claims 1 to 3, wherein:
the step of obtaining the stiffness matrix of the compressor rubber foot pad system comprises the following steps:
acquiring a single stiffness matrix of a single rubber foot pad in the three main shaft directions;
constructing a conversion matrix of the barycenter coordinates and the intersection point coordinates of the single rubber foot pad in the three main shaft directions;
and obtaining the rigidity matrix of the compressor rubber foot pad system according to the monomer rigidity matrix and the conversion matrix.
6. The air conditioner compressor excitation identification method according to any one of claims 1 to 3, wherein:
the step of obtaining the vibration acceleration matrix of the mass center of the compressor according to the vibration acceleration matrix of all the acceleration sensors comprises the following steps:
constructing a transfer matrix of the vibration acceleration of each acceleration sensor and the vibration acceleration of the mass center of the compressor when the compressor runs;
and obtaining the centroid vibration acceleration matrix according to the vibration acceleration matrixes of all the acceleration sensors and the transfer matrix.
7. The air conditioner compressor excitation identification method according to any one of claims 1 to 3, wherein:
after the step of obtaining the centroid equivalent load of the compressor according to the mass matrix, the stiffness matrix and the centroid vibration acceleration matrix, the method further comprises:
and loading the centroid equivalent load to a finite element model of the compressor, and correcting the centroid equivalent load.
8. The air conditioner compressor excitation identification method according to claim 7, wherein:
the step of correcting the centroid equivalent load comprises:
acquiring a corrected mass center vibration acceleration matrix corresponding to the finite element model of the compressor under the preset frequency;
and correcting the centroid equivalent load according to the centroid vibration acceleration matrix and the corrected centroid vibration acceleration matrix.
9. A computer arrangement comprising a processor and a memory, said memory storing a computer program which, when executed by said processor, carries out the steps of the air conditioner compressor excitation recognition method according to any one of claims 1 to 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910521858.5A CN110245425B (en) | 2019-06-17 | 2019-06-17 | Air conditioner compressor excitation identification method and computer device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910521858.5A CN110245425B (en) | 2019-06-17 | 2019-06-17 | Air conditioner compressor excitation identification method and computer device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110245425A CN110245425A (en) | 2019-09-17 |
| CN110245425B true CN110245425B (en) | 2020-10-23 |
Family
ID=67887638
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910521858.5A Active CN110245425B (en) | 2019-06-17 | 2019-06-17 | Air conditioner compressor excitation identification method and computer device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110245425B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110673585B (en) * | 2019-10-15 | 2021-01-22 | 新誉轨道交通科技有限公司 | Test method, test device and test system for train air conditioning system |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62170864A (en) * | 1986-01-24 | 1987-07-27 | Mitsubishi Electric Corp | Pulse compression radar equipment |
| JP2008128742A (en) * | 2006-11-17 | 2008-06-05 | Toyota Central R&D Labs Inc | Apparatus and method for analyzing vibration in rotating structure |
| CN103884482A (en) * | 2012-12-21 | 2014-06-25 | 珠海格力电器股份有限公司 | Compressor-based vibration test method and system |
| CN104791234A (en) * | 2015-05-04 | 2015-07-22 | 合肥工业大学 | Load exciting test analytical method for rotor compressor of refrigeration equipment under starting working condition |
| CN106815446A (en) * | 2017-01-24 | 2017-06-09 | 合肥工业大学 | When variable working condition under refrigeration plant rotor compressor load excitation discrimination method |
| CN108920804A (en) * | 2018-06-25 | 2018-11-30 | 四川长虹空调有限公司 | Refrigeration equipment frequency-changeable compressor excitation load emulated computation method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10830230B2 (en) * | 2017-01-04 | 2020-11-10 | Haier Us Appliance Solutions, Inc. | Method for operating a linear compressor |
-
2019
- 2019-06-17 CN CN201910521858.5A patent/CN110245425B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62170864A (en) * | 1986-01-24 | 1987-07-27 | Mitsubishi Electric Corp | Pulse compression radar equipment |
| JP2008128742A (en) * | 2006-11-17 | 2008-06-05 | Toyota Central R&D Labs Inc | Apparatus and method for analyzing vibration in rotating structure |
| CN103884482A (en) * | 2012-12-21 | 2014-06-25 | 珠海格力电器股份有限公司 | Compressor-based vibration test method and system |
| CN104791234A (en) * | 2015-05-04 | 2015-07-22 | 合肥工业大学 | Load exciting test analytical method for rotor compressor of refrigeration equipment under starting working condition |
| CN106815446A (en) * | 2017-01-24 | 2017-06-09 | 合肥工业大学 | When variable working condition under refrigeration plant rotor compressor load excitation discrimination method |
| CN108920804A (en) * | 2018-06-25 | 2018-11-30 | 四川长虹空调有限公司 | Refrigeration equipment frequency-changeable compressor excitation load emulated computation method |
Non-Patent Citations (2)
| Title |
|---|
| Active Surge Control for Magnetically Suspended Centrifugal Compressors Using a Variable Equilibrium Point Approach;Xin Ma 等;《IEEE Transactions on Industrial Electronics ( Volume: 66 , Issue: 12 , Dec. 2019 )》;20190114;第9383-9393页 * |
| 基于虚拟样机的转子式压缩机瞬态工况激励载荷辨识方法;曹兴枫 等;《合肥工业大学学报(自然科学版)》;20180728;第41卷(第7期);第891-894页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110245425A (en) | 2019-09-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6813749B2 (en) | Method, system and computer program product for multidisciplinary design analysis of structural components | |
| EP1131768B1 (en) | Generating a nonlinear model and generating drive signals for simulation testing using the same | |
| JPH10267800A (en) | Durability evaluation method for wheel drum | |
| Giagopoulos et al. | Dynamic response and identification of critical points in the superstructure of a vehicle using a combination of numerical and experimental methods | |
| CN103983412A (en) | Avionic device operating modal measuring method for vibration finite element model correction | |
| CN110245425B (en) | Air conditioner compressor excitation identification method and computer device | |
| US20150073730A1 (en) | Mechanical strain gauge simulation | |
| CN110516334B (en) | Convolution calculation simulation test method and device based on hardware environment and related equipment | |
| US9098673B2 (en) | Structural optimization for vehicle crashworthiness | |
| CN112665810B (en) | Method and system for determining vibration shedding of chip, storage medium and electronic equipment | |
| JP7036209B2 (en) | Diagnostic equipment, diagnostic methods, and programs | |
| CN113138059B (en) | Resonance endurance assessment method, device, equipment, system and storage medium | |
| US20210334435A1 (en) | Method and device for simulating a technical system | |
| JP2004045294A (en) | Determination system and program for risk of damaging structure | |
| JP2008286545A (en) | Method and apparatus for simulating rubber part | |
| CN107678967B (en) | Unit test coverage rate generation method and device, readable storage medium and equipment | |
| CN108312179B (en) | Elastic part testing method and device based on mechanical arm and mechanical arm | |
| CN112035781B (en) | Method, device, equipment and storage medium for predicting endurance life of rubber vibration isolator | |
| JP6552938B2 (en) | Method, apparatus and program for calculating hysteresis loss of filler filled rubber | |
| CN114545500B (en) | Method and device for determining wave velocity distribution information | |
| JP6346526B2 (en) | Deformation simulation method for tire and spring assembly | |
| US20240020439A1 (en) | Device and method for analyzing durability of drive line for vehicle | |
| CN113609590A (en) | Method, device and equipment for predicting noise in vehicle and computer storage medium | |
| CN114297896A (en) | Method, device, equipment and medium for determining landing collision stiffness | |
| Sahin et al. | MATLAB for all steps of dynamic vibration test of structures |
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 |