CN113361030A - Method for establishing rotary machine dynamic balance calculation application system based on LabVIEW - Google Patents

Method for establishing rotary machine dynamic balance calculation application system based on LabVIEW Download PDF

Info

Publication number
CN113361030A
CN113361030A CN202110616249.5A CN202110616249A CN113361030A CN 113361030 A CN113361030 A CN 113361030A CN 202110616249 A CN202110616249 A CN 202110616249A CN 113361030 A CN113361030 A CN 113361030A
Authority
CN
China
Prior art keywords
calculation
dynamic balance
vector
module
labview
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
Application number
CN202110616249.5A
Other languages
Chinese (zh)
Inventor
陶传龙
蒋顺军
刘卫军
黄建利
张喜林
刘从敏
何凯
张亚朴
屈兴东
赵爽
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huaneng Shandong Shidaobay Nuclear Power Co Ltd
Original Assignee
Huaneng Shandong Shidaobay Nuclear Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huaneng Shandong Shidaobay Nuclear Power Co Ltd filed Critical Huaneng Shandong Shidaobay Nuclear Power Co Ltd
Priority to CN202110616249.5A priority Critical patent/CN113361030A/en
Publication of CN113361030A publication Critical patent/CN113361030A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M1/00Testing static or dynamic balance of machines or structures
    • G01M1/14Determining imbalance
    • G01M1/16Determining imbalance by oscillating or rotating the body to be tested
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling

Landscapes

  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

A rotary machine dynamic balance calculation application system based on LabVIEW is established, through LabVIEW programming, the common calculation functions of field dynamic balance calculation such as vector calculation, single-sided dynamic balance, double-sided dynamic balance, trial/balance weight quality estimation and the like are integrated, the key calculation and implementation process of the field dynamic balance is simplified and programmed, in the vibration analysis and diagnosis and fault treatment of field equipment, only a plurality of important relevant data are input, the accurate unbalance fault position and balance weight result can be calculated and given, the traditional methods such as manual vector method calculation, drawing method and the like are replaced, convenience, rapidness and accuracy are realized, and the working efficiency is greatly improved; the application program and the installation package file are generated through LabVIEW programming, the installation and use of the application program and the installation package file can be realized on a common conventional computer, the program interface is simple, the operation is clear, the program can be used on other computers without LabVIEW software, and the application range is wide.

Description

Method for establishing rotary machine dynamic balance calculation application system based on LabVIEW
Technical Field
The invention belongs to the technical field of vibration testing, and particularly relates to a method for establishing a rotary machine dynamic balance calculation application system based on LabVIEW.
Background
Vibration monitoring is one of the most effective methods in condition monitoring related work, and vibration faults are also a type of faults which often occur in a production field and are difficult to handle. By carrying out vibration monitoring and analysis on the equipment, the running condition of the equipment can be evaluated, a basis is provided for troubleshooting of fault reasons, and guidance is provided for improving the vibration state. The unbalance of the rotating machinery is the most common fault in the industrial production field, the mass unbalance of the rotor system is the main reason of the vibration fault, and the practice shows that the dynamic balance implemented in the field is one of the effective ways for solving the unbalance vibration fault.
At present, the conventional methods for calculating and implementing the on-site dynamic balance of the rotary machine comprise a mapping method or a graphical method, a vector calculation method, an influence coefficient method and the like, the methods are relatively complex to calculate and implement, if the vector calculation method and the influence coefficient method are used for manual calculation, the workload is large, errors are easy to make mistakes, and errors of the graphical method are relatively large; some advanced vibration analyzers integrate dynamic balance calculation programs, but the dynamic balance calculation programs can only calculate vibration data measured by the vibration analyzer, but cannot calculate and apply measurement data of other equipment, and in addition, the vibration analyzer with dynamic balance calculation is high in price; some units outsource the service work needing field dynamic balance to professional vibration service companies, but the price is relatively high, and meanwhile, because the units are not resident on the field, the units cannot guarantee that the treatment is finished in the shortest time when emergency treatment is needed.
Because on-site mechanical equipment structure restriction still often need carry out some other calculations when carrying out dynamic balance, for example the balance calculates and can't install the balancing weight (structural restriction or occupation) in the position that obtains to add the counter weight, consequently need carry out vector calculation and decompose it, be convenient for install near available position in order to reach the same balanced effect. The problems can be solved by manual calculation, but the manual calculation is time-consuming, labor-consuming and error-prone, inconvenience is brought to field vibration problem treatment, and work efficiency is reduced.
Disclosure of Invention
In order to overcome the defects in the prior art, the invention aims to provide the method for establishing the rotary machine dynamic balance calculation application system based on the LabVIEW, which is used for reducing the complexity of calculation such as dynamic balance and the like on site, improving the fault processing efficiency and conveniently and quickly eliminating the problem of unbalance of the rotary machine by vibration processing personnel; the method has the advantages of simple interface, convenience in operation, reliability in calculation and wide applicability.
In order to achieve the purpose, the invention adopts the technical scheme that:
a method for establishing a rotary machine dynamic balance calculation application system based on LabVIEW specifically comprises the following steps:
step one, manufacturing a vector calculation basic module
1) Vector calculation
The simple harmonic vibration can be represented by the projection of a rotation vector on a rectangular coordinate axis, the size of the rotation vector is the amplitude of the simple harmonic vibration, the angular velocity is the circular frequency of the simple harmonic vibration, and the included angle between the angular velocity and a horizontal axis is the phase angle of the simple harmonic vibration; and (3) performing vector addition and subtraction calculation by using a trigonometric function calculation method: assuming that the two vectors have a magnitude of A, B and phases of α and β, respectively, when the two vectors are added, the magnitude C is:
Figure BDA0003097727030000021
the angle θ is:
Figure BDA0003097727030000022
when these two vectors are subtracted, the magnitude D is:
Figure BDA0003097727030000023
angle of rotation
Figure BDA0003097727030000024
Comprises the following steps:
Figure BDA0003097727030000025
similarly, by using the trigonometric function calculation method, the vector can be decomposed into two specified directions, and assuming that the magnitude of a certain vector is E and the phase angle is gamma, and the vector needs to be decomposed into two directions of delta and epsilon, the magnitude F of the vector in the delta direction is
Figure BDA0003097727030000026
A size G in the epsilon direction of
Figure BDA0003097727030000027
2) Establishment of vector calculation basic module
a making vector addition calculation basic module
Placing a numerical value input control and a numerical value display control on a LabVIEW front panel, renaming the corresponding controls for data input and output display, and displaying the controls created on the front panel in a corresponding program block diagram; editing the vector addition calculation formula (1) into a program diagram by using LabVIEW built-in sine, cosine and arc tangent calculation and conventional addition, subtraction, multiplication, division, square root and other calculators so as to realize vector addition operation; the module is packaged into a basic module, so that calling in subsequent programming is facilitated;
b making a vector subtraction calculation base module
Placing a numerical value input control and a numerical value display control on a LabVIEW front panel, renaming the corresponding controls for data input and output display, and displaying the controls created on the front panel in a corresponding program block diagram; editing the vector subtraction calculation formula (2) into a program block diagram by using LabVIEW built-in sine, cosine and arc tangent calculation and conventional addition, subtraction, multiplication, division, square root and other calculators so as to realize vector subtraction operation; the module is packaged into a basic module, so that calling in subsequent programming is facilitated;
c making angle conversion basic module
After vector calculation, for the angles which are obtained in the step 1) and exceed the range of 0-360 degrees, the angles do not conform to the daily use habit, the calculated angles are converted into the range of 0-360 degrees by manufacturing an angle conversion basic module, and the method is conveniently applied to conversion of output results of all angles in the subsequent calculation process, and specifically comprises the following steps of: nesting two layers of condition structures in the FOR loop structure: the first layer of condition structure setting judgment form: if the input data is larger than 360, outputting the data after subtracting 360, and if the input data is smaller than or equal to 360, entering a second-layer condition structure, wherein the second-layer condition structure is set with a judgment form: if the input data is less than 0, the input data is added with 360 degrees and then output, and if the input data is more than or equal to 0, the input data is directly output;
step two, programming design of vector calculation function module is carried out
Designing multi-vector addition calculation and separate vector subtraction and decomposition calculation functions, and simultaneously, according to a centrifugal force formula F ═ m omega2r, wherein F represents centrifugal force, m represents mass of the balancing weight, omega represents rotation angular velocity, r represents radius of the balancing weight, and the radius r is used as an adjusting coefficient for realizing vector calculation of different radii; the specific method comprises the following steps: placing a plurality of numerical value input controls on a LabVIEW front panel andthe numerical display controls are named according to corresponding physical meanings, the controls are connected with corresponding ports of the vector addition calculation basic module, the vector subtraction calculation basic module and the angle conversion basic module which are manufactured in the step one in a corresponding program block diagram, the physical meanings of the control names and the corresponding ports of the calculation basic modules are ensured to be consistent, the controls are placed in a While circulating structure, the waiting time is set according to the actual situation, the continuous operation of the program is facilitated, and the vector calculation function is realized;
step three, establishing a single-sided and double-sided dynamic balance calculation module, and carrying out detailed design on the single-sided and double-sided dynamic balance calculation module;
step four, establishing a trial/counterweight estimation module;
step five, main interface arrangement
Placing a vector calculation Boolean switch, a single-sided dynamic balance Boolean switch, a double-sided dynamic balance Boolean switch, a trial/balance weight estimation Boolean switch and an end Boolean switch on a LabVIEW front panel; placing an event structure in a corresponding program block diagram, placing corresponding edited function modules of vector calculation, single-sided dynamic balance, double-sided dynamic balance and trial/balance weight mass estimation under corresponding event branches of the event structure, setting the event as that the corresponding Boolean switch value is changed, executing the event under the corresponding branch, namely entering the corresponding function module, and placing the function module in a while loop; a sequence structure is placed in the program block diagram, and the sequence structure is designed to firstly carry out Boolean switch to restore the default value and then execute while circulation corresponding content during each operation, so that the program can be continuously executed conveniently; the corresponding function module can be entered or the program function can be exited by clicking the corresponding icon on the front panel of the main interface;
step six, generating an installation program
And generating a modularized application program and an installation package file by using a LabVIEW self-contained application program generator.
And step one, 2) nested two condition structures are circularly calculated FOR multiple times in the FOR loop so as to ensure that the output result is in the range of 0-360 after data passes through the module.
In the third step, a single-sided and double-sided dynamic balance calculation module is established, and specifically:
1) single-sided dynamic balance calculation module establishment
Under the working rotating speed or the balance rotating speed of the equipment, the initial vibration at the bearing is measured to
Figure BDA0003097727030000031
Applying a test weight to a correction plane
Figure BDA0003097727030000032
After starting to the same rotation speed, measuring the vibration
Figure BDA0003097727030000033
The device is tested
Figure BDA0003097727030000034
The resulting exacerbation effect can be expressed as
Figure BDA0003097727030000035
Coefficient of influence
Figure BDA0003097727030000041
Take off the test weight
Figure BDA0003097727030000042
Applying a balancing weight
Figure BDA0003097727030000043
Just can balance out initial vibration
Figure BDA0003097727030000044
Then
Figure BDA0003097727030000045
Namely, it is
Figure BDA0003097727030000046
If actually aggravated
Figure BDA0003097727030000047
Residual vibration
Figure BDA0003097727030000048
Similarly, according to the known influence coefficient, the size and the position of the balance weight can be directly deduced without adding a test weight;
placing a plurality of numerical value input controls and numerical value display controls on a LabVIEW front panel, connecting corresponding controls with a vector addition calculation basic module, a vector subtraction calculation basic module and a trigonometric function calculation module of the LabVIEW in a corresponding program block diagram, editing the calculation process of the single-sided dynamic balance calculation formula, converting the obtained angle into a range of 0-360 degrees through the angle conversion basic module manufactured in the step one, placing the angle in a While circulating structure, and facilitating continuous operation of a program according to the actual setting of waiting time, namely realizing the single-sided dynamic balance calculation;
2) double-sided dynamic balance calculation module establishment
When the rotor is weighted at I, II two calibration planes, the A, B bearings near I, II planes are affected differently and the initial vibration at the bearings is measured at balanced speed
Figure BDA0003097727030000049
Figure BDA00030977270300000410
Will try out
Figure BDA00030977270300000411
After being applied to the plane I, the second start-up, vibration after test weight measurement
Figure BDA00030977270300000412
Take down the test weight
Figure BDA00030977270300000413
Will try out
Figure BDA00030977270300000414
Applying to plane II, third starting to measure vibration after test weight
Figure BDA00030977270300000415
Figure BDA00030977270300000416
From this it can be calculated: the influence coefficients of the weighting on the A, B bearing on the I surface are respectively
Figure BDA00030977270300000417
The influence coefficients of the weighting on the A, B bearing on the II surface are respectively
Figure BDA00030977270300000418
Suppose weighting on I, II sides
Figure BDA00030977270300000419
The residual vibration of the two bearings after balance can be 0
Figure BDA00030977270300000420
Figure BDA00030977270300000421
Placing a plurality of numerical value input controls and numerical value display controls on a LabVIEW front panel, connecting corresponding controls with a vector addition basic module, a vector subtraction basic module and a trigonometric function calculation module of the LabVIEW in a corresponding program block diagram, editing the calculation process of the double-sided dynamic balance calculation formula into a program, converting the obtained angle into a range of 0-360 degrees through the angle conversion basic module manufactured in the step I, placing the converted angle into a While circulating structure, and facilitating continuous operation of the program according to actual set waiting time, namely realizing double-sided dynamic balance calculation.
The step three is a detailed design of the single and double-sided dynamic balance calculation module, which specifically comprises the following steps:
1) refinement design of single-side dynamic balance calculation module
a no influence coefficient calculation function
When the field dynamic balance work is carried out, a trial-weight test is required, namely an influence coefficient is calculated through an original vibration value and a vibration value after the trial weight is increased, the test needs to start the shutdown equipment for at least 1 time, and the trial-weight test is not required under the condition of known influence coefficient, so that the startup and shutdown times can be reduced; according to the basic principle of the dynamic balance work, the known influence coefficients are used as input parameters to be edited into a LabVIEW single-sided dynamic balance calculation program, so that the function of calculating directly through inputting the influence coefficients is realized;
b balance data comparison function
When dynamic balance calculation is carried out, any data of a certain measuring point in the horizontal or vertical direction is used for carrying out balance calculation, and calculation needs to be carried out for a plurality of times in order to compare the calculation results of data in two directions, so that no influence coefficient and influence coefficient calculation program compilation needs to be carried out repeatedly on a single-sided dynamic balance calculation program, and after an edited program diagram is copied, pasted and the names of all controls are modified, two parallel mutually independent calculation sub-modules can be formed, and the calculation results obtained by using different input data are visually compared;
c vector diagram display function
The LabVIEW self-contained two-dimensional compass diagram is used for displaying the vector diagram in the program block diagram, process data generated by the first row of calculation in the non-influence coefficient dynamic balance calculation in the single-side dynamic balance calculation module is connected to the two-dimensional compass diagram, the first row of calculation data in the non-influence coefficient dynamic balance calculation is displayed as the vector diagram, and the comparison with the conventional vector diagram solution dynamic balance is realized;
d data screenshot saving function
The dynamic balance calculation data needs to be stored as equipment data, and a dynamic balance data storage function is designed: a plurality of numerical value input controls are arranged on a front panel of the program and are used for basic parameters of balanced equipment, such as equipment names, equipment KKS codes, names of dynamic balance operators and other information needing to be recorded; LabVIEW is used to display the operation time with 'format date/time character string'; placing a screenshot Boolean switch on a front panel, using an event structure in a program block diagram, and executing an event in an event structure branch when an editing event branch is that a screenshot button value is changed; using LabVIEW with 'create report' and 'application Front Panel Image' modules in the event structure branch, and capturing a Front Panel when executing the event;
2) the double-sided dynamic balance calculation module is designed in a detailed mode, only the function of storing the data screenshots is added to double-sided dynamic balance calculation, and the design method of the double-sided dynamic balance calculation module is the same as that of the single-sided dynamic balance calculation module.
Establishing a trial/counterweight estimation module in the fourth step, which specifically comprises:
1) establishing a trial weight quality estimation module: reference formula
Figure BDA0003097727030000051
M represents the test weight mass, M represents the rotor mass, g represents the gravity acceleration, f represents the rotation frequency, e represents the weighted radius, and the test weight mass is estimated according to the basic parameters of the balanced equipment;
placing a plurality of numerical value input controls and numerical value display controls on a LabVIEW front panel, and enabling the corresponding controls and the self-contained mathematical operators of the LabVIEW in the corresponding program block diagram according to a reference formula
Figure BDA0003097727030000052
Connecting, placing the connection in a While circulating structure, and setting waiting time so as to facilitate continuous operation of a program and realize test weight quality estimation;
2) establishing a counterweight mass estimation module: estimating a counterweight mass according to the geometric dimensions and material densities of different counterweights, wherein m represents the counterweight mass, ρ represents the counterweight material density, and V represents the counterweight material volume; the common counterweight materials in the production field comprise a hexagon head bolt, a hexagon nut, a round flat gasket and a square gasket, the common materials of the materials comprise carbon steel, stainless steel, brass and aluminum, the external dimension of the material can be measured by using a ruler, and the volume V of the material can be calculated according to the external dimension;
placing a plurality of numerical value input controls and numerical value display controls on a LabVIEW front panel, connecting corresponding controls in a corresponding program diagram by using mathematical operators carried by the LabVIEW, calculating the volume V of the controls according to different external dimensions, placing density numerical values of common materials in a condition format for selection during calculation, multiplying the calculated volume by the density value of the corresponding material to obtain an output result, and placing the output result in the same weight mass estimation While circulating structure to realize weight mass estimation.
The invention has the following beneficial effects:
1) the invention integrates the common calculation functions of vector calculation, single-sided dynamic balance, double-sided dynamic balance and trial/balance weight mass estimation field dynamic balance calculation, and compared with the traditional methods such as manual vector method calculation, drawing method and the like, the method is convenient, quick and accurate to carry out field dynamic balance by utilizing the invention, thereby greatly improving the working efficiency;
2) and the LabVIEW programming is adopted to simplify and program the key calculation and implementation process of the field dynamic balance. In the vibration analysis and diagnosis and fault treatment of the field equipment, accurate unbalance fault position and counterweight result can be calculated and given as long as a plurality of important relevant data are input;
3) all complex functions of calculation, correction, storage and the like of field implementation of dynamic balance can be completed on a computer, a balance calculation method based on an influence coefficient is provided in single-sided dynamic balance calculation, the starting and stopping times can be reduced under the condition of complete historical data, and the time cost can be greatly saved.
4) In the most common single-side dynamic balance calculation, parallel multi-column data simultaneous calculation is provided, and simultaneously, a vector diagram is provided, so that comparison with a conventional vector diagram solution is facilitated, and the establishment of a balance scheme is greatly facilitated.
5) When the dynamic balance analysis method is used for dynamic balance calculation, a simple vibration analyzer with a phase measurement function is required to be matched, the vibration analyzer does not need to have the balance calculation related function, and the purchase cost of related analysis equipment is reduced.
6) The invention can be installed and used in a common conventional computer by programming the LabVIEW, generating the application program and the installation package file, has simple program interface and clear operation, can be used in other computers without LabVIEW software, and improves the working efficiency.
Drawings
FIG. 1 is a main interface diagram of the present invention.
FIG. 2 is a block diagram of a main interface process of the present invention.
FIG. 3 is a block diagram of a two-vector addition calculation process and terminal diagrams according to the present invention.
FIG. 4 is a block diagram of a two vector subtraction calculation process and terminal icons according to the present invention.
FIG. 5 is a block diagram of angle collection calculation and terminal icons according to the present invention.
FIG. 6 is a block diagram of a vector computing sub-functional module program of the present invention.
FIG. 7 is a vector computing sub-function block operating interface of the present invention.
FIG. 8 is a block diagram of a single dynamic balance function module according to the present invention.
FIG. 9 is a single-sided dynamic balance function module operation interface of the present invention.
FIG. 10 is a block diagram of the dual dynamic balance function module of the present invention.
FIG. 11 is a diagram of a double-sided dynamic balance function module operation interface according to the present invention.
Fig. 12 is a block diagram of a trial/counterweight mass estimation function according to the present invention.
Fig. 13 is an interface for operation of the trial/counterweight mass estimation function module of the present invention.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings.
A method for establishing a rotary machine dynamic balance calculation application system based on LabVIEW specifically comprises the following steps:
step one, manufacturing a vector calculation basic module
1) Vector calculation
The simple harmonic vibration can be represented by the projection of a rotation vector on a rectangular coordinate axis, the size of the rotation vector is the amplitude of the simple harmonic vibration, the angular velocity is the circular frequency of the simple harmonic vibration, and the included angle between the angular velocity and a horizontal axis is the phase angle of the simple harmonic vibration; and (3) performing vector addition and subtraction calculation by using a trigonometric function calculation method: assuming that the two vectors have a magnitude of A, B and phases of α and β, respectively, when the two vectors are added, the magnitude C is:
Figure BDA0003097727030000071
the angle θ is:
Figure BDA0003097727030000072
when these two vectors are subtracted, the magnitude D is:
Figure BDA0003097727030000073
angle of rotation
Figure BDA0003097727030000074
Comprises the following steps:
Figure BDA0003097727030000075
similarly, by using the trigonometric function calculation method, the vector can be decomposed into two specified directions, and assuming that the magnitude of a certain vector is E and the phase angle is gamma, and the vector needs to be decomposed into two directions of delta and epsilon, the magnitude F of the vector in the delta direction is
Figure BDA0003097727030000076
A size G in the epsilon direction of
Figure BDA0003097727030000077
2) Establishment of vector calculation basic module
a making vector addition calculation basic module
Placing a numerical value input control and a numerical value display control on a LabVIEW front panel, renaming the corresponding controls for data input and output display, and displaying the controls created on the front panel in a corresponding program block diagram; editing the vector addition calculation formula (1) into a program diagram by using LabVIEW built-in sine, cosine and arc tangent calculation and conventional addition, subtraction, multiplication, division, square root and other calculators so as to realize vector addition operation; the module is packaged into a basic module, so that calling in subsequent programming is facilitated;
b making a vector subtraction calculation base module
Placing a numerical value input control and a numerical value display control on a LabVIEW front panel, renaming the corresponding controls for data input and output display, and displaying the controls created on the front panel in a corresponding program block diagram; editing the vector subtraction calculation formula (2) into a program block diagram by using LabVIEW built-in sine, cosine and arc tangent calculation and conventional addition, subtraction, multiplication, division, square root and other calculators so as to realize vector subtraction operation; the module is packaged into a basic module, so that calling in subsequent programming is facilitated;
c making angle conversion basic module
After vector calculation, for the angles which are obtained in the step c) of the step 2) and exceed the range of 0-360 degrees, the angles do not conform to the daily use habit, and the calculated angles are converted into the range of 0-360 degrees by manufacturing an angle conversion basic module, so that the method is conveniently applied to conversion of output results of all angles in the subsequent calculation process, and specifically comprises the following steps of: nesting two layers of condition structures in the FOR loop structure: the first layer of condition structure setting judgment form: if the input data is larger than 360, outputting the data after subtracting 360, and if the input data is smaller than or equal to 360, entering a second-layer condition structure, wherein the second-layer condition structure is set with a judgment form: if the input data is less than 0, the input data is added with 360 degrees and then output, and if the input data is more than or equal to 0, the input data is directly output;
step two, programming design of vector calculation function module is carried out
Combining with the actual needs of an industrial field, in the process of carrying out dynamic balance or balance optimization of equipment, vector calculation is needed for facilitating work, and the vector calculation comprises the steps of carrying out multi-vector addition calculation, two-vector subtraction and single-vector angle decomposition calculation to carry out integration or splitting adjustment on balancing weights, and carrying out adjustment calculation on the balancing weights on different balance radiuses;
designing 4 internal vector addition calculation and separate vector subtraction and decomposition calculation functions, and simultaneously, according to a centrifugal force formula F ═ m omega2r, wherein F represents centrifugal force, m represents mass of the balancing weight, omega represents rotation angular velocity, r represents radius of the balancing weight, and the radius r is used as an adjusting coefficient for realizing vector calculation of different radii; the specific method comprises the following steps: placing a plurality of numerical value input controls and a plurality of numerical value display controls on a LabVIEW front panel, naming the numerical value input controls and the numerical value display controls according to corresponding physical meanings, connecting each control with corresponding ports of the vector addition calculation basic module, the vector subtraction calculation basic module and the angle conversion basic module which are manufactured in the first step in a corresponding program block diagram, ensuring that the names of the controls and the physical meanings of the corresponding ports of each calculation basic module are consistent, placing the controls in a While circulating structure, setting waiting time according to reality, facilitating continuous operation of a program and realizing a vector calculation function;
step three, establishing a single-sided and double-sided dynamic balance calculation module, and carrying out detailed design on the single-sided and double-sided dynamic balance calculation module;
establishing a single-sided and double-sided dynamic balance calculation module:
1) establishment of single-side dynamic balance calculation module
Under the working rotating speed or the balance rotating speed of the equipment, the initial vibration at the bearing is measured to
Figure BDA0003097727030000081
Applying a test weight to a correction plane
Figure BDA0003097727030000082
After starting to the same rotation speed, measuring the vibration
Figure BDA0003097727030000083
The device is tested
Figure BDA0003097727030000084
The resulting exacerbation effect can be expressed as
Figure BDA0003097727030000085
Coefficient of influence
Figure BDA0003097727030000091
Take off the test weight
Figure BDA0003097727030000092
Applying a balancing weight
Figure BDA0003097727030000093
Just can balance out initial vibration
Figure BDA0003097727030000094
Then
Figure BDA0003097727030000095
Namely, it is
Figure BDA0003097727030000096
If actually aggravated
Figure BDA0003097727030000097
Residual vibration
Figure BDA0003097727030000098
Similarly, according to the known influence coefficient, the size and the position of the balance weight can be directly deduced without adding a test weight;
placing a plurality of numerical value input controls and numerical value display controls on a LabVIEW front panel, connecting corresponding controls with a vector addition calculation and vector subtraction calculation basic module manufactured in the step one and a trigonometric function calculation module of the LabVIEW in a corresponding program block diagram, editing the calculation process of the single-sided dynamic balance calculation formula, converting the obtained angle into a range of 0-360 degrees through the angle conversion basic module manufactured in the step one, placing the angle in a While circulating structure, and facilitating continuous operation of a program according to the actual set waiting time, namely realizing the single-sided dynamic balance calculation;
2) establishment of double-sided dynamic balance calculation module
I, II two at two ends of the rotorWhen the weight is added on the correction plane, A, B two bearings close to I, II surfaces are affected differently, and the initial vibration at the two bearings is measured at the balanced rotating speed
Figure BDA0003097727030000099
Figure BDA00030977270300000910
Will try out
Figure BDA00030977270300000911
After being applied to the plane I, the second start-up, vibration after test weight measurement
Figure BDA00030977270300000912
Take down the test weight
Figure BDA00030977270300000913
Will try out
Figure BDA00030977270300000914
Applying to plane II, third starting to measure vibration after test weight
Figure BDA00030977270300000915
Figure BDA00030977270300000916
From this it can be calculated: the influence coefficients of the weighting on the A, B bearing on the I surface are respectively
Figure BDA00030977270300000917
The influence coefficients of the weighting on the A, B bearing on the II surface are respectively
Figure BDA00030977270300000918
Suppose weighting on I, II sides
Figure BDA00030977270300000919
The residual vibration of the two bearings after balance can be 0
Figure BDA00030977270300000920
Figure BDA00030977270300000921
Placing a plurality of numerical value input controls and numerical value display controls on a LabVIEW front panel, connecting corresponding controls with a vector addition and vector subtraction basic module and a trigonometric function calculation module carried by the LabVIEW in a corresponding program block diagram, editing the calculation process of the double-sided dynamic balance calculation formula into a program, converting the obtained angle into a range of 0-360 degrees through the angle conversion basic module manufactured in the step I, placing the converted angle into a While circulating structure, and facilitating the continuous operation of the program according to the actual set waiting time, namely realizing double-sided dynamic balance calculation;
refined design for single-sided and double-sided dynamic balance calculation module
1) Refinement design of single-side dynamic balance calculation module
a no influence coefficient calculation function
When the field dynamic balance work is carried out, a trial-weight test is required, namely an influence coefficient is calculated through an original vibration value and a vibration value after the trial weight is increased, the test needs to start the shutdown equipment for at least 1 time, and the trial-weight test is not required under the condition of known influence coefficient, so that the startup and shutdown times can be reduced; according to the basic principle of the dynamic balance work, the known influence coefficients are used as input parameters to be edited into a LabVIEW single-sided dynamic balance calculation program, so that the function of calculating directly through inputting the influence coefficients is realized;
b balance data comparison function
When dynamic balance calculation is carried out, any data of a certain measuring point in the horizontal or vertical direction is used for carrying out balance calculation, and calculation needs to be carried out for a plurality of times in order to compare the calculation results of data in two directions, so that no influence coefficient exists in a single-side dynamic balance calculation program, the calculation program of the influence coefficient is compiled, and for the convenience of work, after an edited program diagram is copied, pasted and the names of all controls are modified, two parallel calculation sub-modules which are independent of each other can be formed, and the calculation results obtained by using different input data are visually compared;
c vector diagram display function
In order to compare with the conventional vector graphic solution dynamic balance, a LabVIEW self-contained two-dimensional compass diagram is used for displaying the vector diagram in a program diagram, process data generated by the first column in the non-influence coefficient dynamic balance calculation in a single-side dynamic balance calculation module is connected to the two-dimensional compass diagram, and the first column of calculation data in the non-influence coefficient dynamic balance calculation is displayed as the vector diagram;
d data screenshot saving function
The dynamic balance calculation data needs to be stored as equipment data, and a dynamic balance data storage function is designed: a plurality of numerical value input controls are arranged on a front panel of the program and are used for basic parameters of balanced equipment, such as equipment names, equipment KKS codes, names of dynamic balance operators and other information needing to be recorded; LabVIEW is used to display the operation time with 'format date/time character string'; placing a screenshot Boolean switch on a front panel, using an event structure in a program block diagram, and executing an event in an event structure branch when an editing event branch is that a screenshot button value is changed; using LabVIEW with 'create report' and 'application Front Panel Image' modules in the event structure branch, and capturing a Front Panel when executing the event;
2) the double-sided dynamic balance calculation module is designed in a detailed mode, only the function of storing the data screenshots is added to double-sided dynamic balance calculation, and the design method of the double-sided dynamic balance calculation module is the same as that of the single-sided dynamic balance calculation module.
Step four, establishing a trial/balance weight estimation module
1) The trial weight quality estimation module is established as follows: reference formula
Figure BDA0003097727030000101
M represents the test weight mass, M represents the rotor mass, g represents the gravity acceleration, f represents the rotation frequency, e represents the weighted radius, and the test weight mass is estimated according to the basic parameters of the balanced equipment;
placing a plurality of numerical value input controls and numerical value displays on a LabVIEW front panelShowing a control element, and enabling the corresponding control element and a mathematic operator carried by LabVIEW to be in accordance with a reference formula in a corresponding program block diagram
Figure BDA0003097727030000102
Connecting, placing the connection in a While circulating structure, and setting waiting time so as to facilitate continuous operation of a program and realize test weight quality estimation;
2) and establishing a counterweight mass estimation module: estimating a counterweight mass according to the geometric dimensions and material densities of different counterweights, wherein m represents the counterweight mass, ρ represents the counterweight material density, and V represents the counterweight material volume; the common counterweight materials in the production field comprise a hexagon head bolt, a hexagon nut, a round flat gasket and a square gasket, the common materials of the materials comprise carbon steel, stainless steel, brass and aluminum, the external dimension of the material can be measured by using a ruler, and the volume V of the material can be calculated according to the external dimension;
placing a plurality of numerical value input controls and numerical value display controls on a LabVIEW front panel, connecting corresponding controls in a corresponding program diagram by using mathematical operators carried by the LabVIEW, calculating the volume V of the controls according to different external dimensions, placing density numerical values of common materials in a condition format for selection during calculation, multiplying the calculated volume by the density value of the corresponding material to obtain an output result, and placing the output result in the same weight mass estimation While circulating structure to realize weight mass estimation;
step five, main interface arrangement
Placing a vector calculation Boolean switch, a single-sided dynamic balance Boolean switch, a double-sided dynamic balance Boolean switch, a trial/balance weight estimation Boolean switch and an end Boolean switch on a LabVIEW front panel; placing an event structure in a corresponding program block diagram, placing corresponding edited function modules of vector calculation, single-sided dynamic balance, double-sided dynamic balance and trial/balance weight mass estimation under corresponding event branches of the event structure, setting the event as that the corresponding Boolean switch value is changed, executing the event under the corresponding branch, namely entering the corresponding function module, and placing the function module in a while loop; a sequence structure is placed in the program block diagram, and the sequence structure is designed to firstly carry out Boolean switch to restore the default value and then execute while circulation corresponding content during each operation, so that the program can be continuously executed conveniently; the corresponding function module can be entered or the program function can be exited by clicking the corresponding icon on the front panel of the main interface;
step six, generating an installation program
The LabVIEW self-contained application program generator is used for generating the modularized application program and the installation package file, so that the LabVIEW self-contained application program generator can be conveniently used on other computers without LabVIEW software.
And step one, 2) nested two condition structures are circularly calculated FOR multiple times in the FOR loop so as to ensure that the output result is in the range of 0-360 after data passes through the module.
Referring to fig. 1 and 2, fig. 1 is a main interface of an application program of the present invention, and fig. 2 is a program block diagram corresponding to the main interface, when the present invention is used to perform related calculations, a corresponding calculation function module can be accessed by clicking a button of the corresponding function module on the interface.
Referring to fig. 3, 4 and 5, fig. 3 is a block diagram of a calculation procedure for adding two vectors and a terminal icon, fig. 4 is a block diagram of a calculation procedure for subtracting two vectors and a terminal icon, and fig. 5 is a block diagram of a calculation procedure for collecting angles and a terminal icon. The three figures show the design method of three basic modules of vector addition, vector subtraction and angle collection, namely, the related mathematical calculation formula is converted into a program readable language, and related input and output terminals are designed, so that the method is convenient to use when a functional module is edited subsequently.
Referring to fig. 6 and 7, a program block diagram of a vector computing subfunction module of fig. 6, and an operation interface of a vector computing subfunction module of fig. 7. Fig. 6 shows a method for implementing the vector calculation function module, that is, mathematical relationships such as common vector addition, vector subtraction, vector decomposition, etc. (including the basic module manufactured as described above) are converted into program readable languages, and the calculation function module can be pushed out to return to the main interface when the "home" button value is changed by the condition structure. Fig. 7 shows a vector calculation sub-function module operation interface, which can be specifically used for performing the following operations as required: 1. vector addition calculation: the size and the angle of the vector are respectively input at the value 1 and the value 2 (the default radius is 1, and the default radius can be modified according to the actual situation), so that the result of the addition of the two vectors can be obtained at the vector sum; in the same way, the addition calculation of 4 vectors can be simultaneously carried out, and the result is obtained once. 2. And (3) vector subtraction calculation: similar to vector addition, the corresponding vector is input at the position of vector subtraction calculation, and the operation result can be obtained. 3. And (3) vector decomposition calculation: and inputting the vector to be decomposed and the angle to be decomposed, and obtaining the size of the decomposed vector on the corresponding angle. After the use is finished, clicking a 'homepage' button, exiting the calculation module and returning to the program main interface. In order to facilitate data use, when the 'homepage' button is used for exiting the calculation module, the result calculated by using the module is not cleared, the 'vector calculation' button is clicked again on the main interface, and when the calculation module is entered again, the related calculation result is still reserved and can be continuously used.
Referring to fig. 8 and 9, fig. 8 is a single-sided dynamic balance function module program block diagram, and fig. 9 is an operation interface of the single-sided dynamic balance function module. Fig. 8 shows a specific implementation method of the single-sided dynamic balance function module, that is, a mathematical formula for single-sided dynamic balance calculation is converted into a programming language, and meanwhile, according to the use requirement, practical functions such as balance data comparison, vector diagram display, no-influence coefficient calculation, data storage and the like are added. Fig. 9 is a single-side operation interface of the dynamic balance function module, and a specific using method of the calculation module is briefly described as follows:
for new equipment or equipment without balance record, the balance calculation can be carried out by using the function of the part of the single-sided dynamic balance sub-interface without reference influence coefficient. Inputting the measured initial 1 frequency multiplication vibration amplitude value of the rotating equipment at the position of initial vibration 1, and inputting the initial 1 frequency multiplication vibration phase value of the rotating equipment at the position of initial phase 1; inputting a test weight at the position of the test weight 1, and inputting the angle of the added test weight on the rotor at the test weight angle; after the apparatus is restarted, the vibration amplitude and phase after the trial weight are measured, and the values are input at "vibration after trial weight 1" and "phase after trial weight 1". At the moment, a dynamic balance scheme can be obtained, such as the parameters of the weight and the angle which need to be added after the test weight is removed, the weight and the angle which need to be added after the test weight is reserved, the influence coefficient of the equipment and the like, and a vector diagram for changing the dynamic balance work is drawn on a vector diagram.
The parameters can also be input in the initial vibration 2 part by adopting vibration data in different directions, and compared with a balance scheme.
For the equipment with the historical data of successful trial weighing and the influence coefficient of the same kind of equipment, partial functions of 'having reference influence coefficient' in a single-side dynamic balance sub-interface can be used for balance calculation.
The measured initial 1-frequency multiplication vibration amplitude value and phase of the rotating equipment are input at the initial vibration 3 and the initial phase 3, and known influence coefficients are respectively input at the reference influence coefficient 3 and the reference coefficient phase 3, so that a balance scheme can be directly obtained. In order to compare the influences of different influence coefficients, corresponding numerical values can be input at the reference influence coefficient 4 and the reference coefficient phase 4, and the calculation results can be compared, so that the balance scheme is optimized conveniently.
After the balance is finished, the parameters of the equipment and the data required to be stored in the test process can be input at the lower part of the interface of FIG. 9, and the operator can click the screenshot button to store the screenshot of the dynamic balance data.
The single-side dynamic balance calculation function module is used for solving balance weight and position by a vector analysis method and obtaining an influence coefficient; the satisfactory effect can be achieved only by one-time counterweight adjustment, the operation is convenient, the interface is visual, and the calculation is accurate; residual vibration calculation is provided, and a balance scheme is conveniently formulated according to the actual situation on site; the vector diagram is provided, comparison with the traditional vector diagram method is convenient, and meanwhile, the vector diagram can improve the data analysis efficiency of vibration analysts; two completely independent dynamic balance calculation input columns are provided, so that vibration analysis personnel can conveniently compare vibration data; meanwhile, the data storage function is realized; generally, historical data and influence coefficients are successfully tested by similar equipment, and balancing efficiency and accuracy are greatly improved.
Referring to fig. 10 and 11, fig. 10 is a block diagram of a double-sided dynamic balance function module, and fig. 11 is an operation interface of the double-sided dynamic balance function module. The related calculations and operations are similar to the single sided dynamic balance calculations and are not described in detail herein.
Referring to fig. 12 and 13, fig. 12 is a block diagram of a trial/counterweight mass estimation function block, and fig. 13 is an operation interface of the trial/counterweight mass estimation function block. Fig. 12 shows a method for implementing trial/weight mass estimation calculation, i.e., converting mathematical formulas into programming language, in the weight mass estimation calculation, a conditional format structure is used to distinguish common weight materials, so that a user can directly select the weight materials by using a pull-down list when using the weight mass estimation calculation. Fig. 13 is an operation interface of the computing module, and the specific use method is as follows: 1. the trial weight estimation is based on an empirical formula, the force generated according to the trial weight is estimated to be one tenth of the static weight of the rotor, and when the test weight is used, the related parameters of the balanced rotor are respectively input into input frames of 'the weight of the rotor', 'the radius of the counterweight' and 'the balanced rotating speed', so that the trial weight can be obtained; 2. the counterweight mass estimation uses a ruler to measure the sizes of common counterweight blocks such as bolts, nuts, circular gaskets, rectangular gaskets and the like, inputs the sizes into a corresponding input frame, selects materials of the common counterweight blocks, and inputs the using quantity of each counterweight block to obtain the total weight of the counterweight blocks.
Since no chamfer, thread, etc. is taken into account and the density is not well defined according to the material of the specific composition, the estimated quality results will deviate from reality. For example, a plain carbon steel bolt of M12 x 50 has an actual mass of about 60.6g, and the mass estimated using the present invention is 60.9g, and such mass deviations are negligible at actual equilibrium.
In the application process of the invention, the dynamic balance work of field complex conditions can be conveniently completed by matching and using all the functional modules, the following processes of implementing a dynamic single-sided balance process by equipment without balance record, implementing a single-sided dynamic balance process by equipment with reference influence coefficients and implementing a double-sided dynamic balance are respectively used, and the matching and using method of all the modules is briefly described:
1) the equipment without balance record implements dynamic single-side balance process (corresponding data is input into the input box with name and tail number 1 without special description, hereinafter referred to as data line 1)
s1. when the rotary machine is stopped, a reflective sheet is adhered to the exposed part of the rotary shaft of the device, such as the position of the coupling, to serve as the phase zero point on the rotary shaft, to provide pulse signals for the laser phase sensor, and the phase angle is increased in the reverse rotation direction, that is, if the rotary shaft rotates clockwise, the cursor is used as the zero point, and the counterclockwise rotation direction is used as the angle increase;
s2, aligning the laser phase sensor with a reflector adhered to the rotating shaft by using a dial indicator frame, connecting the laser phase sensor with a vibration analyzer, slightly shaking the laser phase sensor to enable a laser spot to sweep the reflector back and forth, confirming that a rotating speed display exists on the vibration analyzer, and ensuring that a phase measuring system is available;
s3. installing vibration sensor on the bearing shell of the tested device, using magnetic seat to adsorb to collect the vibration signal of the bearing seat, calculating in the direction of larger vibration amplitude during balancing;
s4., confirming the safety of the field environment, and starting the tested equipment to the working rotation speed;
s5. collecting vibration signal of the tested device at working speed by using a vibration analyzer, checking a frequency doubling vibration amplitude and phase in a spectrogram, namely an initial vibration amplitude and an initial vibration phase, and then stopping the device;
s6., analyzing the vibration signal vibration data of the tested equipment collected at s5 at the working speed, and determining that the vibration problem of the equipment can accord with the unbalance characteristic;
s7. opening the program interface of the invention, clicking the 'trial/balance weight estimation' button on the main interface, entering the trial/balance weight mass estimation operator program interface, obtaining the trial weight mass by inputting the rotor mass, the balance radius, namely the distance between the position of the balance weight and the rotation center, and the balance rotation speed, namely the current equipment rotation speed, wherein the mass is estimated by an empirical formula for reference, or determining the trial weight mass according to experience, and clicking the 'homepage' button to return to the main interface after memorizing the mass;
s8. clicking a single-sided dynamic balance button on the main interface, entering a single-sided dynamic balance subprogram interface, inputting initial vibration amplitude nuclear initial vibration phase values acquired in the step s5 in the input boxes of initial vibration 1 and initial phase 1, inputting test weight mass obtained by calculation in the step s7 in the input box of test weight 1, and inputting corresponding angles in the input box of test weight angle 1 according to analysis of vibration data and combination of experience;
s9. applying the test weight in step s8 on the corresponding position of the balancing disk or the shaft coupling of the tested device;
s10, starting the rotating driven equipment to a working rotating speed again, collecting vibration data of each bearing vibration measuring point after test weight is increased, and finishing the balancing process if the vibration is reduced to a qualified level; because the sensitivity of each rotating device to the unbalanced response is inconsistent, the vibration result is still not ideal, and the change situation of the 1x amplitude and the phase of the rotating device needs to be judged: when the phase change is less than 30 degrees and the 1x amplitude change is less than 20 percent, the trial weight quality or the angle selection is not matched, and the trial weight needs to be carried out again at the moment; when the phase change is more than 30 degrees and the 1x amplitude change is more than 20 percent, the test weight is properly selected, the operation of the tested equipment is stopped, and the measured amplitude and phase are respectively filled into the input boxes of 'vibration after test weight 1' and 'phase after test weight 1';
and s11, after data is input, the rotary mechanical dynamic balance calculation application system based on LabVIEW automatically calculates a subsequent balance scheme: if the weight and the angle of the test weight are removed and the weight and the angle of the test weight are reserved, an influence coefficient and an aggravation effect are given;
s12, when the calculated size and the calculated position of the counterweight are not convenient to match or apply on site, inputting approximate mass in the 'correction mass 1' and inputting an approximate angle in the 'correction angle', automatically calculating residual vibration by using a rotary mechanical dynamic balance calculation application system based on LabVIEW, and estimating a balance effect by using the function; when the calculated result is obtained and the counter weight can not be applied to the counter weight position, the 'homepage' button can be clicked to return to the main interface, the data of the single-side dynamic balance interface can not be lost after entering the single-side dynamic balance interface again, and the 'vector calculation' button is clicked to enter the vector calculation sub-interface, such as: inputting 130 weights needing to be added with counter weights in a vector sum input box in a vector decomposition function area, inputting 68 angles needing to be added with counter weights in an angle input box, inputting 60 angles of the vacant counter weight holes in the balance disc in a decomposition angle input box, inputting 90 in a decomposition angle 2 input box, automatically calculating the weights needing to be added with the counter weights in two directions of the angles of the vacant counter weight holes in the balance disc to be 88.9g and 45g respectively by a rotary mechanical dynamic balance calculation application system based on LabVIEW, and returning to a main interface by clicking a home button;
s13, applying the balance weight obtained in the step s11-12 on the balance disc, and restarting the tested device after safety is confirmed;
s14, testing the vibration condition after balancing by using a spectrum analyzer, inputting the weight and the test result into the data line 2 so as to facilitate the comparative analysis of vibration analysis personnel, making a next processing scheme, repeating the steps and finally completing the balancing;
s15, synchronously providing a vector diagram drawn by using the 1 st row of data on a single-side dynamic balance module sub-interface of the rotary mechanical dynamic balance calculation application system based on LabVIEW, and facilitating comparison with a conventional vector graphic method;
s16, after the vibration of the tested equipment is confirmed to be qualified, the equipment is closed or the operation of the equipment is continuously maintained according to the field production requirement;
s17, after the name, the code and the name of a dynamic balance operator of the tested equipment are input and other information needing to be explained are input below a sub-interface of the single-side dynamic balance module, a screenshot button is clicked, a rotary mechanical dynamic balance calculation application system based on LabVIEW can screenshot the current interface, all data are stored for calling, and a program interface can be closed after the data are stored.
2) Single-sided dynamic balancing process with reference influence coefficient equipment
For equipment with balance record or equipment with approximate structure size, carrying out balance calculation by using a region with reference influence coefficient;
s 1-s 5 corresponds to 1) the general case;
s6. clicking a single-sided dynamic balance button on a main interface, entering a single-sided dynamic balance module subprogram interface, inputting the acquired initial vibration amplitude kernel initial vibration phase value in an initial vibration 3 and initial phase 3 input box, inputting the equipment pre-stored influence coefficient or similar equipment reference influence coefficient in a reference influence coefficient 3 and reference coefficient phase 3, and automatically calculating the mass and position of the counterweight to be applied by a rotary mechanical dynamic balance calculation application system based on LabVIEW;
s 7-the subsequent steps of s17 are identical to those in the case of 1) in general, and are not described in detail here.
3) Double-sided dynamic balance implementation process
s1. when the rotary machine is stopped, a reflective sheet is adhered to the exposed part of the rotary shaft of the device to serve as the phase zero point of the rotary shaft, to provide pulse signals for the laser phase sensor, and the phase angle is increased in the reverse rotation direction, i.e. if the rotary shaft rotates clockwise, the cursor is used as the zero point, and the counterclockwise rotation direction is increased;
s2, aligning the laser phase sensor with a reflector adhered to the rotating shaft by using a dial indicator frame, connecting the sensor with a vibration analyzer, slightly shaking the laser phase sensor to enable a laser spot to sweep the reflector back and forth, confirming that a rotating speed display exists on the vibration analyzer, and ensuring that a phase measuring system is available;
s3. vibration sensors are respectively mounted on the bearing housings at the two ends of the tested device, usually by magnetic seat adsorption, for collecting the vibration signal of the bearing seat (calculation is carried out in the direction of larger vibration amplitude during balancing);
s4., confirming the safety of the field environment, and starting the tested equipment to the working rotation speed;
s5. collecting vibration signal of the tested device at working speed by using a vibration analyzer, checking a frequency doubling vibration amplitude and phase in a spectrogram, namely an initial vibration amplitude and an initial vibration phase, and then stopping the device;
s6. analyzing the vibration data to determine that the equipment vibration problem can meet the imbalance characteristic;
s7. clicking the button of double-sided dynamic balance on the main interface, entering the interface of the double-sided dynamic balance subprogram, and inputting the initial vibration and initial phase measured by two planes under the item of initial vibration;
s8. applying a test weight 1 to the corresponding position of the balance disc on the A surface of the tested device;
s9., the rotating passive equipment is started again to work speed, vibration data of each bearing vibration measuring point after test weight is added is collected, and the measured amplitude and phase are filled in the corresponding input frames under the item of 'A surface test weight vibration';
s10, stopping the machine, taking down the trial weight of the surface A, and applying a trial weight 2 to the surface B (generally, the weight of the trial weight is consistent with that of the trial weight 1);
s11, starting the machine again, testing the vibration conditions of the two bearing seats, inputting the vibration conditions into a condition of 'test weight vibration on the B surface', and stopping the machine from running;
s12, after data are input, the rotary mechanical dynamic balance calculation application system based on LabVIEW can automatically calculate a subsequent balance scheme and give an influence coefficient and an aggravation effect, and meanwhile, a residual vibration result is given for a vibration analyst to refer to if a certain aggravation is given;
s13, after the calculation balance weights are applied to the two balance discs, restarting the tested equipment after safety is confirmed;
and S14, testing the vibration condition after balancing by using a spectrum analyzer, and repeating the steps to finally complete the balancing.

Claims (5)

1. A method for establishing a rotary machine dynamic balance calculation application system based on LabVIEW is characterized in that: the method specifically comprises the following steps:
step one, manufacturing a vector calculation basic module
1) Vector calculation
The simple harmonic vibration can be represented by the projection of a rotation vector on a rectangular coordinate axis, the size of the rotation vector is the amplitude of the simple harmonic vibration, the angular velocity is the circular frequency of the simple harmonic vibration, and the included angle between the angular velocity and a horizontal axis is the phase angle of the simple harmonic vibration; and (3) performing vector addition and subtraction calculation by using a trigonometric function calculation method: assuming that the two vectors have a magnitude of A, B and phases of α and β, respectively, when the two vectors are added, the magnitude C is:
Figure FDA0003097727020000011
the angle θ is:
Figure FDA0003097727020000012
when these two vectors are subtracted, the magnitude D is:
Figure FDA0003097727020000013
angle of rotation
Figure FDA0003097727020000014
Comprises the following steps:
Figure FDA0003097727020000015
similarly, by using the trigonometric function calculation method, the vector can be decomposed into two specified directions, and assuming that the magnitude of a certain vector is E and the phase angle is gamma, and the vector needs to be decomposed into two directions of delta and epsilon, the magnitude F of the vector in the delta direction is
Figure FDA0003097727020000016
A size G in the epsilon direction of
Figure FDA0003097727020000017
2) Establishment of vector calculation basic module
a making vector addition calculation basic module
Placing a numerical value input control and a numerical value display control on a LabVIEW front panel, renaming the corresponding controls for data input and output display, and displaying the controls created on the front panel in a corresponding program block diagram; editing the vector addition calculation formula (1) into a program diagram by using LabVIEW built-in sine, cosine and arc tangent calculation and conventional addition, subtraction, multiplication, division, square root and other calculators so as to realize vector addition operation; the module is packaged into a basic module, so that calling in subsequent programming is facilitated;
b making a vector subtraction calculation base module
Placing a numerical value input control and a numerical value display control on a LabVIEW front panel, renaming the corresponding controls for data input and output display, and displaying the controls created on the front panel in a corresponding program block diagram; editing the vector subtraction calculation formula (2) into a program block diagram by using LabVIEW built-in sine, cosine and arc tangent calculation and conventional addition, subtraction, multiplication, division, square root and other calculators so as to realize vector subtraction operation; the module is packaged into a basic module, so that calling in subsequent programming is facilitated;
c making angle conversion basic module
After vector calculation, for the angles which are obtained in the step 1) and exceed the range of 0-360 degrees, the angles do not conform to the daily use habit, the calculated angles are converted into the range of 0-360 degrees by manufacturing an angle conversion basic module, and the method is conveniently applied to conversion of output results of all angles in the subsequent calculation process, and specifically comprises the following steps of: nesting two layers of condition structures in the FOR loop structure: the first layer of condition structure setting judgment form: if the input data is larger than 360, outputting the data after subtracting 360, and if the input data is smaller than or equal to 360, entering a second-layer condition structure, wherein the second-layer condition structure is set with a judgment form: if the input data is less than 0, the input data is added with 360 degrees and then output, and if the input data is more than or equal to 0, the input data is directly output;
step two, programming design of vector calculation function module is carried out
Designing multi-vector addition calculation and separate vector subtraction and decomposition calculation functions, and simultaneously, according to a centrifugal force formula F ═ m omega2r, wherein F represents centrifugal force, m represents mass of the balancing weight, omega represents rotation angular velocity, r represents radius of the balancing weight, and the radius r is used as an adjusting coefficient for realizing vector calculation of different radii; the specific method comprises the following steps: placing a plurality of numerical value input controls and a plurality of numerical value display controls on a LabVIEW front panel, naming the numerical value input controls and the numerical value display controls according to corresponding physical meanings, and calculating the basic modules of each control and the vector addition prepared in the first step in the corresponding program block diagramThe corresponding ports of the vector subtraction calculation basic module and the angle conversion basic module are connected, so that the physical meanings of the names of the controls and the corresponding ports of the calculation basic modules are ensured to be consistent, the controls and the corresponding ports of the calculation basic modules are placed in a While circulating structure, the programs can run continuously conveniently according to the actual setting of waiting time, and the vector calculation function is realized;
step three, establishing a single-sided and double-sided dynamic balance calculation module, and carrying out detailed design on the single-sided and double-sided dynamic balance calculation module;
step four, establishing a trial/counterweight estimation module;
step five, main interface arrangement
Placing a vector calculation Boolean switch, a single-sided dynamic balance Boolean switch, a double-sided dynamic balance Boolean switch, a trial/balance weight estimation Boolean switch and an end Boolean switch on a LabVIEW front panel; placing an event structure in a corresponding program block diagram, placing corresponding edited function modules of vector calculation, single-sided dynamic balance, double-sided dynamic balance and trial/balance weight mass estimation under corresponding event branches of the event structure, setting the event as that the corresponding Boolean switch value is changed, executing the event under the corresponding branch, namely entering the corresponding function module, and placing the function module in a while loop; a sequence structure is placed in the program block diagram, and the sequence structure is designed to firstly carry out Boolean switch to restore the default value and then execute while circulation corresponding content during each operation, so that the program can be continuously executed conveniently; the corresponding function module can be entered or the program function can be exited by clicking the corresponding icon on the front panel of the main interface;
step six, generating an installation program
And generating a modularized application program and an installation package file by using a LabVIEW self-contained application program generator.
2. The method for building a LabVIEW-based rotating mechanical dynamic balance computing application system according to claim 1, wherein: and step one, 2) nested two condition structures are circularly calculated FOR multiple times in the FOR loop so as to ensure that the output result is in the range of 0-360 after data passes through the module.
3. The method for building a LabVIEW-based rotating mechanical dynamic balance computing application system according to claim 1, wherein: in the third step, a single-sided and double-sided dynamic balance calculation module is established, and specifically:
1) single-sided dynamic balance calculation module establishment
Under the working rotating speed or the balance rotating speed of the equipment, the initial vibration at the bearing is measured to
Figure FDA0003097727020000031
Applying a test weight to a correction plane
Figure FDA0003097727020000032
After starting to the same rotation speed, measuring the vibration
Figure FDA0003097727020000033
The device is tested
Figure FDA0003097727020000034
The resulting exacerbation effect can be expressed as Δ
Figure FDA0003097727020000035
Coefficient of influence
Figure FDA0003097727020000036
Take off the test weight
Figure FDA0003097727020000037
Applying a balancing weight
Figure FDA0003097727020000038
Just can balance out initial vibration
Figure FDA0003097727020000039
Then
Figure FDA00030977270200000310
Namely, it is
Figure FDA00030977270200000311
If actually aggravated
Figure FDA00030977270200000312
Residual vibration
Figure FDA00030977270200000313
Similarly, according to the known influence coefficient, the size and the position of the balance weight can be directly deduced without adding a test weight;
placing a plurality of numerical value input controls and numerical value display controls on a LabVIEW front panel, connecting corresponding controls with a vector addition calculation basic module, a vector subtraction calculation basic module and a trigonometric function calculation module of the LabVIEW in a corresponding program block diagram, editing the calculation process of the single-sided dynamic balance calculation formula, converting the obtained angle into a range of 0-360 degrees through the angle conversion basic module manufactured in the step one, placing the angle in a While circulating structure, and facilitating continuous operation of a program according to the actual setting of waiting time, namely realizing the single-sided dynamic balance calculation;
2) double-sided dynamic balance calculation module establishment
When the rotor is weighted at I, II two calibration planes, the A, B bearings near I, II planes are affected differently and the initial vibration at the bearings is measured at balanced speed
Figure FDA00030977270200000314
Figure FDA00030977270200000315
Will try out
Figure FDA00030977270200000316
After being applied to the plane I, the second start-up, vibration after test weight measurement
Figure FDA00030977270200000317
Take down the test weight
Figure FDA00030977270200000318
Will try out
Figure FDA00030977270200000319
Applying to plane II, third starting to measure vibration after test weight
Figure FDA00030977270200000320
Figure FDA00030977270200000321
From this it can be calculated: the influence coefficients of the weighting on the A, B bearing on the I surface are respectively
Figure FDA00030977270200000322
The influence coefficients of the weighting on the A, B bearing on the II surface are respectively
Figure FDA00030977270200000323
Suppose weighting on I, II sides
Figure FDA00030977270200000324
The residual vibration of the two bearings after balance can be 0
Figure FDA00030977270200000325
Figure FDA00030977270200000326
Placing a plurality of numerical value input controls and numerical value display controls on a LabVIEW front panel, connecting corresponding controls with a vector addition basic module, a vector subtraction basic module and a trigonometric function calculation module of the LabVIEW in a corresponding program block diagram, editing the calculation process of the double-sided dynamic balance calculation formula into a program, converting the obtained angle into a range of 0-360 degrees through the angle conversion basic module manufactured in the step I, placing the converted angle into a While circulating structure, and facilitating continuous operation of the program according to actual set waiting time, namely realizing double-sided dynamic balance calculation.
4. The method for building a LabVIEW-based rotating mechanical dynamic balance computing application system according to claim 1, wherein: the step three is a detailed design of the single and double-sided dynamic balance calculation module, which specifically comprises the following steps:
1) refinement design of single-side dynamic balance calculation module
a no influence coefficient calculation function
When the field dynamic balance work is carried out, a trial-weight test is required, namely an influence coefficient is calculated through an original vibration value and a vibration value after the trial weight is increased, the test needs to start the shutdown equipment for at least 1 time, and the trial-weight test is not required under the condition of known influence coefficient, so that the startup and shutdown times can be reduced; according to the basic principle of the dynamic balance work, the known influence coefficients are used as input parameters to be edited into a LabVIEW single-sided dynamic balance calculation program, so that the function of calculating directly through inputting the influence coefficients is realized;
b balance data comparison function
When dynamic balance calculation is carried out, any data of a certain measuring point in the horizontal or vertical direction is used for carrying out balance calculation, and calculation needs to be carried out for a plurality of times in order to compare the calculation results of data in two directions, so that no influence coefficient and influence coefficient calculation program compilation needs to be carried out repeatedly on a single-sided dynamic balance calculation program, and after an edited program diagram is copied, pasted and the names of all controls are modified, two parallel mutually independent calculation sub-modules can be formed, and the calculation results obtained by using different input data are visually compared;
c vector diagram display function
The LabVIEW self-contained two-dimensional compass diagram is used for displaying the vector diagram in the program block diagram, process data generated by the first row of calculation in the non-influence coefficient dynamic balance calculation in the single-side dynamic balance calculation module is connected to the two-dimensional compass diagram, the first row of calculation data in the non-influence coefficient dynamic balance calculation is displayed as the vector diagram, and the comparison with the conventional vector diagram solution dynamic balance is realized;
d data screenshot saving function
The dynamic balance calculation data needs to be stored as equipment data, and a dynamic balance data storage function is designed: a plurality of numerical value input controls are arranged on a front panel of the program and are used for basic parameters of balanced equipment, such as equipment names, equipment KKS codes, names of dynamic balance operators and other information needing to be recorded; LabVIEW is used to display the operation time with 'format date/time character string'; placing a screenshot Boolean switch on a front panel, using an event structure in a program block diagram, and executing an event in an event structure branch when an editing event branch is that a screenshot button value is changed; using LabVIEW with 'create report' and 'application Front Panel Image' modules in the event structure branch, and capturing a Front Panel when executing the event;
2) the double-sided dynamic balance calculation module is designed in a detailed mode, only the function of storing the data screenshots is added to double-sided dynamic balance calculation, and the design method of the double-sided dynamic balance calculation module is the same as that of the single-sided dynamic balance calculation module.
5. The method for building a LabVIEW-based rotating mechanical dynamic balance computing application system according to claim 1, wherein: establishing a trial/counterweight estimation module in the fourth step, which specifically comprises:
1) establishing a trial weight quality estimation module: reference formula
Figure FDA0003097727020000051
M represents the test weight mass, M represents the rotor mass, g represents the gravity acceleration, f represents the rotation frequency, e represents the weighted radius, and the test weight mass is estimated according to the basic parameters of the balanced equipment;
placing a plurality of numerical value input controls and numerical value display controls on a LabVIEW front panel, and enabling the corresponding controls and the self-contained mathematical operators of the LabVIEW in the corresponding program block diagram according to a reference formula
Figure FDA0003097727020000052
Connecting, placing the connection in a While circulating structure, and setting waiting time so as to facilitate continuous operation of a program and realize test weight quality estimation;
2) establishing a counterweight mass estimation module: estimating a counterweight mass according to the geometric dimensions and material densities of different counterweights, wherein m represents the counterweight mass, ρ represents the counterweight material density, and V represents the counterweight material volume; the common counterweight materials in the production field comprise a hexagon head bolt, a hexagon nut, a round flat gasket and a square gasket, the common materials of the materials comprise carbon steel, stainless steel, brass and aluminum, the external dimension of the material can be measured by using a ruler, and the volume V of the material can be calculated according to the external dimension;
placing a plurality of numerical value input controls and numerical value display controls on a LabVIEW front panel, connecting corresponding controls in a corresponding program diagram by using mathematical operators carried by the LabVIEW, calculating the volume V of the controls according to different external dimensions, placing density numerical values of common materials in a condition format for selection during calculation, multiplying the calculated volume by the density value of the corresponding material to obtain an output result, and placing the output result in the same weight mass estimation While circulating structure to realize weight mass estimation.
CN202110616249.5A 2021-06-02 2021-06-02 Method for establishing rotary machine dynamic balance calculation application system based on LabVIEW Pending CN113361030A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110616249.5A CN113361030A (en) 2021-06-02 2021-06-02 Method for establishing rotary machine dynamic balance calculation application system based on LabVIEW

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110616249.5A CN113361030A (en) 2021-06-02 2021-06-02 Method for establishing rotary machine dynamic balance calculation application system based on LabVIEW

Publications (1)

Publication Number Publication Date
CN113361030A true CN113361030A (en) 2021-09-07

Family

ID=77531370

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110616249.5A Pending CN113361030A (en) 2021-06-02 2021-06-02 Method for establishing rotary machine dynamic balance calculation application system based on LabVIEW

Country Status (1)

Country Link
CN (1) CN113361030A (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6102965A (en) * 1996-09-23 2000-08-15 National Instruments Corporation System and method for providing client/server access to graphical programs
CN101059386A (en) * 2007-06-04 2007-10-24 西安交通大学 Precession vector based holographic on-spot dynamic balance method
CN104568313A (en) * 2015-01-20 2015-04-29 湖南科技大学 Influence coefficient dynamic balance method on rotating machine with multiple plane, multiple-points and multiple revolving speed shafting
CN105890843A (en) * 2016-04-18 2016-08-24 神华集团有限责任公司 Dynamic balance method and dynamic balance device
CN110579312A (en) * 2019-10-17 2019-12-17 江苏方天电力技术有限公司 dynamic balance fault detection method for multi-wheel-disc shafting of non-trial-weight rotating machinery

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6102965A (en) * 1996-09-23 2000-08-15 National Instruments Corporation System and method for providing client/server access to graphical programs
CN101059386A (en) * 2007-06-04 2007-10-24 西安交通大学 Precession vector based holographic on-spot dynamic balance method
CN104568313A (en) * 2015-01-20 2015-04-29 湖南科技大学 Influence coefficient dynamic balance method on rotating machine with multiple plane, multiple-points and multiple revolving speed shafting
CN105890843A (en) * 2016-04-18 2016-08-24 神华集团有限责任公司 Dynamic balance method and dynamic balance device
CN110579312A (en) * 2019-10-17 2019-12-17 江苏方天电力技术有限公司 dynamic balance fault detection method for multi-wheel-disc shafting of non-trial-weight rotating machinery

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
张帝: "矢量法现场动平衡研究", 《中国设备工程》 *
操文盛: "基于LabVIEW的电机转子动平衡***设计", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 *
殷祥超: "振动理论与测试技术", 《振动理论与测试技术 *
王展 等: "高速主轴在线动平衡控制策略与实验研究", 《机床与液压》 *

Similar Documents

Publication Publication Date Title
US10802679B2 (en) Methods and apparatus for configuring a data analyzer
JPS6317167B2 (en)
CN110646141B (en) Office software-based rotary machine on-site dynamic balance calculation and adjustment method
CN106484623A (en) A kind of method of software test, apparatus and system
US10718689B2 (en) Modeling and visualization of vibration mechanics in residual space
CN109115408B (en) Dynamic balance test method of large hydroelectric generating set based on centrifugal force equation
CN113361030A (en) Method for establishing rotary machine dynamic balance calculation application system based on LabVIEW
Kniel et al. Detecting 6 DoF geometrical errors of rotary tables
Ivanova et al. The role of information technology in rotor balancing
CN101192242A (en) Crankshaft parametric computer modeling system
Zhu et al. Compilation of program-loading spectrum for milling of a motorized spindle based on cutting force model
CN102494846A (en) Single-side dynamic balance tool
US4345472A (en) Method and apparatus for digitally analyzing dynamic unbalance of a rotating body
US6622105B2 (en) Dynamic correlation extension for a self-balancing rotatable apparatus
EP2939124A1 (en) Methods and apparatus for defining a probe configuration using a probe configuration tool
CN104697484B (en) A kind of gear rotation angle measurement method and device
CN111038514B (en) Vehicle speed control method and related device
Al Najjar et al. Science federation emulation testbed: Demonstration of vfsie functionalities
CN207147459U (en) Gears tooth pitch measurement apparatus
CN114329674A (en) Cloud high-throughput material automatic screening method and system based on Quantum Espresso
Zhao et al. A Dynamic‐Balancing Testing System Designed for Flexible Rotor
Kyi et al. Development of a user-friendly web-based rainfall-runoff model
CN1264035A (en) Holographic dynamic balance method for asymmetric rotor
Postlethwaite et al. Geometric error analysis software for CNC machine tools
Manzato et al. Order Based Modal Analysis Versus Standard Techniques to Extract Modal Parameters of Operational Wind Turbine Gearboxes

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
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20210907

WD01 Invention patent application deemed withdrawn after publication