JP5622177B2 - How to obtain influence coefficient - Google Patents

Description

  The present invention relates to an influence coefficient acquisition method for acquiring an influence coefficient of a rotating machine.

  The rotating machine is, for example, a fluid machine in which rotating blades that exert force on a fluid are provided on a rotating body. There are a prime mover and a driven machine in the fluid machine. The prime mover converts the energy of the fluid into rotational kinetic energy by rotationally driving the rotating body by the pressure that the fluid acts on the rotor blades. As a prime mover, for example, there is a gas turbine (axial turbine, radial turbine). The driven machine applies rotational kinetic energy to the fluid by rotating the rotor blades that are rotationally driven to apply pressure to the fluid. Examples of the driven machine include a compressor (an axial flow compressor, a mixed flow compressor, a cross flow compressor, and a pump provided in a centrifugal compressor, an aircraft engine, or the like). In addition, some fluid machines have a turbocharger that functions as both a prime mover and a driven machine.

  The influence coefficient indicates the vibration change of the rotating body with respect to the balance change of the rotating body provided in the rotating machine. The influence coefficient is used for correcting the balance of the rotating body. Therefore, the influence coefficient of the rotating machine is acquired before correcting the balance of the rotating body.

  The influence coefficient can be obtained by using a test weight as follows, for example. First, without using the trial weight, the rotating body is rotated, and the vibration of the supporting body that supports the rotating body is measured. Next, the trial weight is attached to the rotating body, the rotating body is rotated, and the vibration of the support that supports the rotating body is measured. Then, the influence coefficient is calculated from the vibration when the trial weight is not used, the vibration when the trial weight is attached, and the mass and attachment position of the trial weight. The vibration used for calculating the influence coefficient is preferably a vibration having the same frequency component as the rotation speed of the rotating body (that is, the number of rotations per second).

  Using the influence coefficient obtained in this way, the balance of the rotating body is corrected as follows. The imbalance data of the rotating body is calculated using the influence coefficient. At the correction position indicated by the unbalance data, the rotating body is cut by the mass indicated by the unbalance data. Thereby, the balance of the rotating body is corrected. Such a balance correction method is called an influence coefficient method, and is described, for example, in Patent Document 1 below.

JP 2008-102049 A

  When the above-described balance correction is performed on a plurality of rotating machines of the same model, unbalance data of the rotating body of each rotating machine is calculated using the same influence coefficient.

  In this case, the accuracy of balance correction may be lowered. Since a plurality of rotating machines of the same model have processing errors and assembly errors, there are individual differences in the influence coefficients of these rotating machines. Accordingly, the unbalanced data calculated using the influence coefficient having individual differences also includes an error. Therefore, the accuracy of balance correction may be lowered.

  As a result, there is a possibility that the number of times that the balance correction is performed again for one rotating machine increases, or that the rotating machine cannot be shipped as a product without completing the balance correction.

  Accordingly, an object of the present invention is to efficiently obtain an influence coefficient that enables the balance correction of each rotating machine to be performed with high accuracy when the balance correction is performed on a plurality of rotating machines of the same model. Is to provide a way

In order to achieve the above object, according to the present invention, an influence coefficient acquisition method for acquiring an influence coefficient indicating a vibration change of a rotating body with respect to a balance change of the rotating body provided in a rotating machine,
(A) The rotating body of the rotating machine is rotated, and in this state, the vibration of the support body that supports the rotating body is measured, the rotation angle of the rotating body is measured, and the vibration is based on the measured vibration and the rotation angle. Generate data,
(B) Calculate unbalance data based on the vibration data generated in (A) and the basic influence coefficient,
(C) cutting the rotating body based on the unbalance data;
(D) Thereafter, the rotating body of the rotating machine is rotated, and in this state, the vibration of the support that supports the rotating body is measured, the rotation angle of the rotating body is measured, and the measured vibration and rotation angle are used. Vibration data,
(E) Vibration data generated in (A) and (D), data of cutting performed in (C), and influence on the rotating machine performed in (A), (C), and (D) Get the relational expression with the coefficient,
(F) The above-mentioned (A) to (E) are performed on a plurality of rotating machines of the same model, and a new influence on the model is obtained on the basis of the plurality of relational expressions obtained in (E) above for these rotating machines. An influence coefficient acquisition method characterized by obtaining a coefficient is provided.

  According to a preferred embodiment of the present invention, after obtaining the new influence coefficient, the above-mentioned (A) to (E) are again performed for the same new rotating machine as the model using the new influence coefficient as the basic influence coefficient. Do.

Preferably, a plurality of unbalance amount ranges indicated by the unbalance data are set in advance,
For each rotating machine, determine the range including the unbalance amount based on the unbalance data calculated in (B) for the rotating machine, and assume that the rotating machine belongs to the range.
In (F), a new influence coefficient is obtained for each of the ranges based on the plurality of relational expressions generated in (E) for the rotating machines belonging to the range.

According to a preferred embodiment of the present invention, in (A),
(A-1) The rotating body of the rotating machine is rotated at the first rotation speed, and in this state, the vibration of the supporting body that supports the rotating body is measured, and the rotation angle of the rotating body is measured and measured. Generate vibration data based on vibration and rotation angle,
(A-2) Next, the rotating body of the rotating machine is rotated at the second rotation speed, and in this state, the vibration of the supporting body that supports the rotating body is measured, and the rotation angle of the rotating body is measured. Generate vibration data based on the measured vibration and rotation angle,
In (B), first and second unbalance data are calculated based on the vibration data generated in (A-1) and (A-2) and the basic influence coefficient,
In (C) above,
(C-1) Based on the first unbalance data, the first cutting target portion in the rotating body is cut,
(C-2) Based on the second unbalance data, the second cutting target portion in the rotating body is cut,
In (D) above,
(D-1) The rotating body of the rotating machine is rotated at the first rotation speed, and in this state, the vibration of the supporting body that supports the rotating body is measured, and the rotation angle of the rotating body is measured and measured. Generate vibration data based on vibration and rotation angle,
(D-2) Next, the rotating body of the rotating machine is rotated at the second rotation speed, and in this state, the vibration of the supporting body that supports the rotating body is measured, and the rotation angle of the rotating body is measured. Generate vibration data based on the measured vibration and rotation angle,
In (E), the vibration data generated in (A-1), (A-2), (D-1), and (D-2), and (C-1) and (C-2) A relational expression between the data of the performed cutting and the influence coefficient regarding the rotating machine that performed the (A), (C), and (D) is generated.

According to another embodiment of the present invention, the first and second vibration sensors are respectively attached to different portions of the support,
In (A), the rotating body of the rotating machine is rotated, and in this state, the vibration of the supporting body that supports the rotating body is measured by the first and second vibration sensors, and the rotation angle of the rotating body is set to an angle. Vibration data is generated based on the vibration measured by the sensor, the vibration measured by the first vibration sensor and the rotation angle measured by the angle sensor, and based on the vibration measured by the second vibration sensor and the rotation angle measured by the angle sensor. Vibration data,
In (B), vibration data generated based on the vibration measured by the first vibration sensor in (A) and vibration data generated based on the vibration measured by the second vibration sensor in (A). And the first and second imbalance data based on the basic influence coefficient,
In (C) above,
(C-1) Based on the first unbalance data, the first cutting target portion in the rotating body is cut,
(C-2) Based on the second unbalance data, the rotary body is cut at the second cutting target part,
In (D), the rotating body of the rotating machine is rotated, and in this state, the vibration of the supporting body that supports the rotating body is measured by the first and second vibration sensors, and the rotation angle of the rotating body is set to the angle. Vibration data is generated based on the vibration measured by the sensor, the vibration measured by the first vibration sensor and the rotation angle measured by the angle sensor, and based on the vibration measured by the second vibration sensor and the rotation angle measured by the angle sensor. Vibration data,
In (E), vibration data generated based on the vibration measured by the first vibration sensor in (A) and vibration data generated based on the vibration measured by the second vibration sensor in (A). Vibration data generated based on the vibration measured by the first vibration sensor in (D), vibration data generated based on the vibration measured by the second vibration sensor in (D), A relational expression with the influence coefficient regarding the rotating machine on which A), (C), and (D) have been performed is generated.

Preferably, after the (D), the rotating machine that has performed the (A) to (D),
(A) calculating unbalance data based on the vibration data generated in (D) and the basic influence coefficient;
(B) When the unbalance amount indicated by the unbalance data is larger than the threshold value, the rotating body of the rotating machine is cut based on the unbalance data,
(C) Thereafter, for the rotating machine, until the unbalance amount becomes equal to or less than the threshold value, generation of vibration data, calculation of unbalance data, and cutting of the rotating body based on the unbalance data are repeated. Then, after cutting the rotating body a set number of times, based on the data of the cutting performed on the rotating machine and the vibration data generated on the rotating machine, calculate the influence coefficient for replacement of the rotating machine,
(D) calculating unbalanced data based on the influence coefficient for replacement and the vibration data generated last for the rotating machine;
(E) The rotating body of the rotating machine is cut based on the unbalance data.

According to the present invention described above, in the balance correction of each rotating machine, vibration data is acquired as in (A), then unbalance data is calculated as in (B), and then (C). Then, the rotating body is cut as described above, and then vibration data is obtained (in order to obtain the unbalance amount after cutting) as in (D) above. Thereafter, as in (E), a relational expression between the vibration data, the cutting data, and the influence coefficient obtained in the process of balance correction is acquired for each rotating machine, and based on these relational expressions, Calculate the impact coefficient. Therefore, a new influence coefficient can be calculated efficiently.
Further, since the new influence coefficient reflects a plurality of influence coefficients specific to a plurality of rotating machines of the same model, the new influence coefficient is a relatively high-precision influence coefficient relating to the rotating machine of the model. Therefore, by using the new influence coefficient for the balance correction of the next rotating machine, the balance can be corrected with relatively high accuracy.

The influence coefficient acquisition apparatus which can be used for the influence coefficient acquisition method of this invention is shown. It is explanatory drawing of vibration data. It is a flowchart which shows the influence coefficient acquisition method by 1st Embodiment or 2nd Embodiment of this invention. (A) is a figure explaining the range of unbalance amount, (B) is a graph which shows the relationship between unbalance amount and the absolute value of an influence coefficient. It is a flowchart which shows the influence coefficient acquisition method by 3rd Embodiment of this invention. The influence coefficient acquisition apparatus in case a rotary machine is a supercharger is shown. The influence coefficient acquisition apparatus in the case of using two vibration sensors is shown.

  An embodiment for carrying out the present invention will be described with reference to the drawings. In addition, the same code | symbol is attached | subjected to the common part in each figure, and the overlapping description is abbreviate | omitted.

  FIG. 1A shows an influence coefficient acquisition apparatus 10 that can be used in the influence coefficient acquisition method of the present invention. FIG. 1B is a BB arrow view of FIG. The influence coefficient acquisition device 10 includes a support 3, a vibration sensor 5, an angle sensor 7, a calculator 9, and a cutting device 11.

  The support 3 supports the rotating body 13 of the rotating machine. The rotating body 13 is rotatable about the central axis C of the rotating body 13 while being supported by the support body 3. In addition, a part of the support body 3 may be comprised by the stationary side member of a rotary machine.

  The vibration sensor 5 is attached to the support 3. The vibration sensor 5 measures the vibration (that is, acceleration, speed, displacement, or load) of the support 3 while the rotating body 13 is rotating, and outputs a vibration signal indicating the vibration to the calculator 9. . As the vibration sensor 5, a known sensor that can be used for obtaining the influence coefficient can be used.

  The angle sensor 7 measures the rotation angle of the rotating body 13 and outputs a rotation angle signal indicating the rotation angle to the calculator 9. The rotation angle changes from zero degrees to 360 degrees when the rotating body 13 rotates once.

The computing unit 9 generates vibration data based on the vibration signal (measured vibration) and the rotation angle signal (measured rotation angle). The vibration data is composed of vibration amplitude and phase θ. FIG. 2A shows the amplitude and phase θ of the vibration. In FIG. 2A, the horizontal axis indicates the rotation angle of the rotating body 13 measured by the angle sensor 7, and the vertical axis indicates the intensity of the primary vibration among the vibrations detected by the vibration sensor 5. The primary vibration is a vibration having the same frequency component as the rotational speed of the rotating body 13. That is, the primary vibration amplitude is a component having the same frequency [Hz] as the rotational speed of the rotating body 13 (the number of rotations per second) at the time of vibration measurement by the vibration sensor 5 from the vibration signal output from the vibration sensor 5. It is the extracted vibration. In FIG. 2A, the phase θ represents the deviation of the primary vibration with respect to the reference rotation angle (zero degree in the example of FIG. 2A). That is, the phase θ represents a deviation of the rotation angle that is the starting point of the period of the primary vibration with respect to the reference rotation angle.
The vibration data is represented by complex numbers. FIG. 2B shows vibration data represented by complex numbers. As shown in FIG. 2B, vibration data is represented by complex numbers with R as the magnitude (absolute value) of the amplitude of the primary vibration and the above-mentioned phase θ as a declination (the same applies hereinafter). The computing unit 9 generates complex vibration data from the vibration signal from the vibration sensor 5 and the rotation angle signal from the angle sensor 7 (hereinafter the same).

  The cutting device 11 includes a cutting tool 11a (for example, an end mill) that cuts the cutting target portion 13a of the rotating body 13, and a three-dimensional (for example, X axis direction in FIG. A drive mechanism 11b that moves in the Y-axis direction and the Z-axis direction) and a position control unit 11c that controls the position of the cutting tool 11a by controlling the operation of the drive mechanism 11b. The position control unit 11c cuts the cutting target portion 13a according to the input unbalance data. The cutting target portion 13a is, for example, a nut or a columnar member attached to the tip of the rotating shaft that constitutes the rotating body 13.

[First Embodiment]
FIG. 3 is a flowchart illustrating an influence coefficient acquisition method according to the first embodiment of the present invention. This influence coefficient acquisition method is also an unbalance correction method for the rotating body 13 and acquires a new influence coefficient during the execution of this unbalance correction method (in the second or third embodiment described later). The same applies to the case).

  In step S1, the rotating body 13 of the m-th rotating machine is rotated, and in this state, the vibration sensor 5 detects the vibration (that is, acceleration, speed, displacement, or load) of the support body 3 that supports the rotating body 13. While measuring, the angle sensor 7 measures the rotation angle. Based on the vibration and rotation angle measured in this way, the calculator 9 generates vibration data a. Note that with respect to the number m of the rotating machine, when performing step S1 for the first rotating machine, m = 1, and thereafter the value of m increases by one each time the process returns to step S1.

In step S2, the unbalance data U is calculated based on the vibration data a generated in step S1 and the basic influence coefficient α. The unbalanced data U is calculated by the arithmetic unit 9 according to the following equation (1).

U = a / α (1)

If step S9, which will be described later, has not been performed yet at the time of performing step S2, the basic influence coefficient α used in step S2 of this time is the same model as the rotating machine currently performing the processing of FIG. It may be an influence coefficient (hereinafter referred to as an initial influence coefficient) acquired in advance by the following equation (2) for the rotating machine.

α = (a _i1 −a _i0 ) / ΔM _T (2)

Here, a _i0 is vibration data generated as described above in a state where the trial weight is not attached to the rotating body 13, and a _i1 is as described above in a state where the trial weight is attached to the rotating body 13. And ΔM _T is expressed by the following equation (3).

ΔM _T = A _T (cos θ _T + j sin θ _T ) (3)

Here, j is an imaginary unit, _AT is the product of the mass of the trial weight and the distance between the position where the trial weight is attached to the rotating body 13 and the center axis C of the rotating body 13, and θ _T is This is a phase indicating the circumferential position where the trial weight is attached to the rotating body 13. The circumferential direction is a direction around the central axis C, the radial direction is a direction with respect to the central axis C, and the axial direction is a direction parallel to the central axis C (hereinafter the same).

When step S9 described later is already performed at the time of performing step S2, the basic influence coefficient used in step S2 is the latest influence coefficient updated in the previous step S9.
When a new influence coefficient is calculated for each range of the unbalance amount in step S9 to be described later when step S2 is performed this time, the following is performed. First, in step S2, the reference unbalance data U0 = a / α is calculated based on the initial influence coefficient α and the vibration data a generated in the immediately preceding step S1. The range including the reference unbalance data U0 is determined. If a new influence coefficient has already been calculated in the previous step S9 for this determined range, in this step S2, the new influence coefficient α and the vibration data a generated in the immediately preceding step S1 are converted. Based on this, unbalanced data U = a / α is calculated. Otherwise, the basic influence coefficient α used in the current step S2 is an initial influence coefficient.

In step S3, the cutting target portion 13a of the rotating body 13 is cut based on the cutting position and the cutting amount indicated by the unbalance data U calculated in step S2. The unbalance data U is input from the computing unit 9 to the position control unit 11c. The position control unit 11c moves the cutting tool 11a by controlling the drive mechanism 11b based on the unbalance data U, whereby the cutting target unit 13a is cut according to the unbalance data U.
The unbalanced data U is expressed by the following equation (4), where A is an absolute value, the declination is θ, and j is an imaginary unit.

U = A (cos θ + jsin θ) (4)

Here, A indicates an unbalance amount, and A is a product of x and y. x is a radial distance from the central axis C of the rotating body 13 to the cutting position in step S3, and y is a mass to be cut in step S3. On the other hand, θ represents the circumferential position of the cutting position. Accordingly, in step S3, the cutting tool 11a moves in the direction of the central axis C until the cutting mass becomes y while being positioned at the radial position and the circumferential position indicated by x and θ. At this time, the position control unit 11c calculates a distance for moving the cutting tool 11a in the direction of the central axis C based on the density of the cutting target part 13a and, if necessary, the shape of the cutting target part 13a, The drive mechanism 11b is controlled based on the distance.
In step S3, the position controller 11c may position the cutting tool 11a at the circumferential position with respect to the direction around the central axis C based on the rotation angle of the angle sensor 7. In step S3, for example, the axial end of the rotating body 13 opposite to the cutting target portion 13a may be gripped by an appropriate means so that the rotating body 13 does not rotate.

  In step S4, the rotating body 13 of the m-th rotating machine is rotated, and in this state, the angle sensor 7 measures the rotation angle while the vibration sensor 5 measures the vibration of the support body 3 that supports the rotating body 13. To do. Based on the vibration and rotation angle measured in this way, the calculator 9 generates vibration data a.

  In step S5, the unbalance data U is calculated by the same method as in step S2, based on the vibration data a and the basic influence coefficient α generated in step S4.

If step S9, which will be described later, has not yet been performed when step S5 is performed, the basic influence coefficient α used in step S5 is the initial influence coefficient.
If step S9, which will be described later, has already been performed at the time of performing step S5, the basic influence coefficient used in step S5 is the latest influence coefficient updated in the previous step S9.
When a new influence coefficient is calculated for each range of the unbalance amount in step S9 described later at the time of performing step S5 this time, the following is performed. First, in this step S5, reference unbalance data U0 = a / α is calculated based on the initial influence coefficient and the vibration data a generated in the immediately preceding step S4. The range including the reference unbalance data U0 is determined. If a new influence coefficient has already been calculated in the previous step S9 for this determined range, in this step S5, the new influence coefficient α and the vibration data a generated in the immediately preceding step S4 are converted. Based on this, unbalanced data U = a / α is calculated. Otherwise, the basic influence coefficient α used in this step S5 is an initial influence coefficient.

  In step S6, it is determined whether the unbalance amount indicated by the unbalance data U calculated in step S5 is equal to or less than a threshold value. If this determination is YES, the rotary machine currently being processed is shipped, and the process proceeds to step S7.

  On the other hand, if the determination in step S6 is no, the process proceeds to step S61. In step S61, it is determined whether or not step S3 has already been performed for the set number of times (two or more times) for the rotating machine currently performing the process of FIG. If the determination in step S61 is yes, it is determined that the balance of the rotating machine cannot be corrected, and the process proceeds to step S7. If the determination in step S61 is no, the process returns to step S3. In the returned step S3, the cutting of the cutting target portion 13a is performed based on the unbalance data U calculated in the immediately preceding step S5, but the other points are the same as described above.

In step S7, the vibration data a _{m, n} , a _{m, n-1} generated in the immediately preceding step S1 or S4, the data of the cutting performed in the immediately preceding step S3, and the influence coefficient α regarding the mth rotating machine, Generate a relational expression.
This relational expression is expressed by the following expression (5).

a _{m, n} −a _{m, n−1} = α × (−ΔM _{m, n} ) (5)

Here, the subscript n or n−1 of the vibration data a indicates the number of times step S3 has already been performed at the time of generating the vibration data, and ΔM _{m, n} represents the vibration data a _m in step S1 or S4. _{, N-1} and before the vibration data am _{, n} are generated, the cutting data in step S3 is shown. That is, ΔM _{m, n} = A (cos θ + jsin θ). Here, A indicates an unbalance amount removed by cutting, and A is represented by a product of x and y. In this case, x is from the central axis C of the rotating body 13 to the cutting position in step S3. It is a radial distance, and y is the mass cut in step S3. On the other hand, θ represents the circumferential position of the cutting position.
Moreover, the subscript m of the vibration data a and the cutting data is the number of the rotating machine that is the target of the current step S7.

With respect to the subscript n of the vibration data a, for example, if the current rotating machine has never returned from step S6 to step S3, n = 1, n−1 = 0, and a _{m, 1} Is the vibration data generated in step S4, am _{, 0} is the vibration data generated in step S1, and ΔM _{m, 1} is the data of the cutting performed in step S3.
Further, when the process returns from step S6 to step S3 k times, (k + 1) number of the relational expressions are generated. For example, when returning from step S6 to step S3 only twice _{, am, 1−am , 0} = α × (−ΔM _{m, 1} ), am _{, 2} −am _{, 1} = α × (−ΔM _{m, 2} ), am _{, 3} −am _{, 2} = α × (−ΔM _{m, 3} ) is generated, and when the process returns from step S6 to step S3 only three times _{, 1−am , 0} = α × (−ΔMm _{, 1} ), am _{, 2−am , 1} = α × (−ΔMm _{, 2} ), am _{, 3−am , 2} = α × ( −ΔM _{m, 3} ), am _{, 4} −am _{, 3} = α × (−ΔM _{m, 4} ).

  In step S8, it is determined whether or not the number of rotating machines that have performed steps S1 to S7 is equal to or greater than a set number (for example, 2). That is, it is determined whether the above m is equal to or greater than the set number. If this determination is YES, the process proceeds to step S9; otherwise, the value of m is incremented by 1, and a rotating body 13 of a new rotating machine of the same model is installed on the support 3; Return to S1. Thereafter, steps S1 to S8 are performed again for the new rotating machine.

  In step S9, a new number of rotating machines of the same model (that is, m rotating machines) equal to or greater than the set number are newly created based on the plurality of relational expressions obtained by performing steps S1 to S7. Calculate the influence coefficient. The basic influence coefficient is updated to the new influence coefficient. That is, the next time the step S2 and step S5 are performed, the new influence coefficient is used as the basic influence coefficient. In addition, the number of the relational expressions used in step S9 increases as the number of rotating machines that performed steps S1 to S7 (that is, the value of m) increases.

  In step S9, preferably, a new influence coefficient is calculated from the plurality of relational expressions by the least square method with the influence coefficient in each relational expression as an unknown. Specifically, a new influence coefficient α is calculated by solving the following [Equation 1].

In [Expression 1], a _{m, n} , a _{m, n−1} , α, ΔM _{m, n} are respectively a _{m, n} , a _{m, n−1} , α, ΔM _m in the above equation (5). _{, N.}
In [Equation 1], Σ _{m, n} represents the sum of all combinations of m and n with respect to the obtained relational expression. For this, when the value of m at the present time and m _c, m is an integer from 1 to m _c, for each m, if the back by k times from step S6 to step S3, n is, It is an integer from 1 to (k + 1).
In [Expression 1], α _R and α _I are a real part and an imaginary part of the influence coefficient α, respectively. The real part of the new influence coefficient α is obtained by solving the upper expression of the above expression (6), and the imaginary part of the new influence coefficient α is obtained by solving the lower expression of the above expression (6). Is obtained.

  When a new influence coefficient is calculated in step S9, the value of m is incremented by one and the process returns to step S1, and the above-described processing is performed again for a new rotating machine of the same model. However, in the next step S2 and step S5, unbalanced data is calculated using the new influence coefficient calculated in the immediately preceding step S9 as the basic influence coefficient.

Preferably, step S9 is performed as follows.
First, a plurality of unbalance amount ranges indicated by the unbalance data are set in advance. For example, three ranges R1, R2, and R3 shown in FIG. 4A are set in advance for the unbalance amount. The range R1 is a range where the unbalance amount is 0 to X1, the range R2 is a range where the unbalance amount is X1 to X2, and the range R3 is a range where the unbalance amount is X2 or more.
Then, for each of the rotating machines, the range including the unbalance amount indicated by the unbalance data calculated in step S2 for the rotating machine is determined, and the rotating machine belongs to the range. In step S9, for each of the ranges, a new influence coefficient is obtained as described above based on the plurality of relational expressions generated in step S7 for the rotating machines belonging to the range.
As shown in FIG. 4B, since the magnitude (absolute value) of the influence coefficient tends to change as the unbalance quantity increases, a new influence coefficient is obtained for each range of the unbalance quantity. Unbalanced data can be calculated with high accuracy using the influence coefficient.

[Second Embodiment]
Next, an influence coefficient acquisition method according to the second embodiment of the present invention using the above-described influence coefficient acquisition apparatus 10 will be described. An influence coefficient acquisition method according to the second embodiment of the present invention will also be described with reference to FIG.

  In the second embodiment, as described below, in step S3, the first and second cutting target portions 13a and 13b in the rotating body 13 are partially cut. The first and second cutting target portions 13a and 13b have different radial positions with respect to the central axis C of the rotating body 13 as in the example of FIG. 1, or have different axial positions. In the example of FIG. 1, the first cutting target portion 13 a is an outer peripheral portion of the tip end portion in the axial direction of the rotating body 13, and the second cutting target portion 13 b is a disk provided at an intermediate portion in the axial direction of the rotating body 13. It is the outer peripheral part of a shaped member.

In step S1, the rotating body 13 of the m-th rotating machine is rotated at the first rotating speed, and in this state, the vibration sensor 5 vibrates the supporting body 3 that supports the rotating body 13 (that is, acceleration, speed, While measuring the displacement or load), the angle sensor 7 measures the rotation angle. Based on the vibration and the rotation angle thus measured, the computing unit 9 generates vibration data a _{1, m, 0} . Note that the first rotation speed described here is the same as each first rotation speed described below.
Furthermore, in step S1, the rotating body 13 of the m-th rotating machine is rotated at a second rotational speed different from the first rotational speed, and in this state, the vibration sensor 5 supports the supporting body 3 that supports the rotating body 13. The angle sensor 7 measures the rotation angle while measuring the vibration. Based on the vibration and the rotation angle thus measured, the computing unit 9 generates vibration data a2 _{, m, 0} . Note that the second rotation speed described here is the same as each second rotation speed described below.
Here, the subscript 1 of the vibration data a indicates the first rotation speed, the subscript 2 of the vibration data a indicates the second rotation speed, and the subscript m of the vibration data a is the m-th rotation. A machine is shown (the same applies hereinafter). The subscript 0 of the vibration data a indicates that step S3 has not been performed for the m-th rotating machine.

In step S2, the vibration data _{a 1} generated in step _{S1, m, 0, a 2} , m, 0 and based on the basic influence coefficient, the calculator 9, first and second unbalanced data _U b, U _d is calculated. The first and second unbalanced data U _b and U _d are expressed by the following determinant of [Equation 2].

In [Equation 2], α _{1, b} , α _{1, d} , α _{2, b} , α _{2, d} are basic influence coefficients, and −1 indicates an inverse matrix.

If step S9, which will be described later, has not yet been performed at the time of performing step S2, the basic influence coefficient used in step S2 of this time is the rotation of the same model as the rotating machine currently performing the processing of FIG. For the machine, the initial influence coefficients α _{1, b} , α _{1, d} , α _{2, b} , α _{2, d} obtained in advance by the following [Equation 3] may be used.

Here, a ₁₀ is not attached trial weight to the rotating body 13, as described above, the vibration generated on the basis of the rotation angle and the vibration measured while rotating the rotary body 13 at the first rotation speed a data, a ₂₀ is not attached trial weight to the rotating body 13, as described above, were generated based on the rotation angle and the vibration measured while rotating the rotary body 13 at a second rotational speed A _1b is vibration data, and a trial weight is attached to the first cutting target portion 13a and, as described above, based on vibration and rotation angle measured in a state where the rotating body 13 is rotated at the first rotation speed. A _2b is vibration and rotation measured with a trial weight attached to the first cutting target portion 13a and the rotating body 13 rotated at the second rotational speed as described above. a vibration data that is generated based on the corner, a _1d The vibration data generated based on the vibration and the rotation angle measured in a state where the trial weight is attached to the second cutting target portion 13b and the rotating body 13 is rotated at the first rotational speed as described above. The vibration data a _2d is generated based on the vibration and rotation angle measured with the trial weight attached to the second cutting target portion 13b and the rotating body 13 rotated at the second rotation speed as described above. Vibration data.
ΔM _{Tb in} [Expression 3] is expressed by the following equation (9).

ΔM _Tb = A _Tb (cos θ _Tb + j sin θ _Tb ) (9)

Here, j is an imaginary unit, and _ATb is the mass of the trial weight attached to the first cutting target portion 13a, the position where the trial weight is attached to the first cutting target portion 13a, and the rotating body 13. It is a product of the distance from the central axis C, and θ _Tb is a phase indicating the circumferential position of the rotating body 13 where the trial weight is attached.
Similarly, ΔM _Td in [Equation 3] is expressed by the following equation (10).

ΔM _Td = A _Td (cos θ _Td + j sin θ _Td ) (10)

Here, j is an imaginary unit, and _ATd is the mass of the trial weight attached to the second cutting target portion 13b, the position where the trial weight is attached to the second cutting target portion 13b, and the rotating body 13. It is a product of the distance from the central axis C, and θ _Td is a phase indicating a circumferential position where the trial weight is attached to the rotating body 13.

  When step S9 described later is already performed at the time of performing step S2, the basic influence coefficient used in step S2 is the latest influence coefficient updated in the previous step S9.

  When a new influence coefficient is calculated for each range of the unbalance amount in step S9 to be described later when step S2 is performed this time, the following is performed.

First, in step S2, based on the initial influence coefficient and vibration data _{a 1} generated in the previous step _{S1, m, 0, a 2} , m, 0, follows by Equation 4], the first and second Calculate reference unbalance data.

In [Equation 4], α _{1, b} , α _{1, d} , α _{2, b} , α _{2, d} are initial influence coefficients, and a _{1, m, 0} and a _{2, m, 0} are the immediately preceding coefficients. The vibration data generated in step S1 in the above, U _b is the first reference unbalance data, and U _d is the second reference unbalance data.

Next, the range including the unbalance amount based on the first and second reference unbalance data U _b and U _d is determined.
When following (C1) described later in step S9, the sum or average of the unbalance amount indicated by the first reference unbalance data U _b and the unbalance amount indicated by the second reference unbalance data U _d It is determined whether the value belongs to one of the preset ranges, and it is assumed that the currently targeted rotating machine belongs to the range to which the sum or average value belongs.
If according to later in step S9 (C2), and unbalance amount indicated by the first reference for unbalanced data U _b, the larger of the unbalanced amount indicated by the second reference for unbalanced data U _d It is determined which of the plurality of ranges set in advance falls within the range, and the currently targeted rotating machine belongs to the range to which the larger value belongs.
If according to the described later in step S9 (C3), for the first reference for unbalanced data U _b, amount of unbalance the first reference for unbalanced data U _b is indicated, the first unbalance data It is determined which of the plurality of ranges set in advance with respect to U _b belongs, and the rotating machine currently targeted for the range to which the unbalance amount indicated by the first reference unbalance data U _b belongs is determined. and it belongs, for the second reference for unbalanced data U _d, unbalance amount which the second reference for unbalanced data U _d is indicated, a plurality of the preset range for the second unbalance data U _d of determines fall within any of the ranges, the current in a range in which the unbalance amount indicated by the second reference for unbalanced data U _d belongs And rotating machinery belongs you are elephants.

When following (C1) or (C2) described later in step S9, if a new influence coefficient has already been calculated in the previous step S9 for the range determined to belong to the rotating machine, this step is performed. In S2, using the new influence coefficient as a basic influence coefficient, the basic influence coefficient and the vibration data a _{1, m, 0} , a _{2, m, 0} generated in the immediately preceding step S1 are expressed in [Equation 2]. By applying, the first and second unbalanced data U _b and U _d are calculated. Otherwise, the basic influence coefficient used in the current step S2 is an initial influence coefficient.

On the other hand, when according to later in step S9 (C3), the first reference for unbalanced data U _b, with respect to the range where the rotating machine is determined to belong, as described above, the first unbalance data U _{b 1} alpha new influence coefficients corresponding _{to, b,} when alpha _{2, b} is already calculated in the previous step S9 is that in this step S2, the new influence coefficients alpha _{1, b,} alpha _{2, b} is a basic influence coefficient α _{1, b} , α _{2, b} corresponding to the first unbalance data U _b , and the second reference unbalance data U _d is determined to belong to the rotating machine as described above. with respect to said ranges, the second unbalanced data U new influence coefficients corresponding to the _d alpha _{1, d,} when the alpha _{2, d} is already calculated in the previous step S9 is that in this step S2, The new influence coefficients alpha _{1, d,} alpha _2, the basic influence coefficient corresponding to _d in the second unbalance data _{U d α 1, d, α} 2, and _d, the basic influence coefficient α _{1, b,} α _{2, b 1} , α _{1, d} , α _{2, d} and the vibration data a _{1, m, n} , a _{2, m, n} generated in the immediately preceding step S ₁ are applied to [Equation 2], First and second unbalanced data U _b and U _d are calculated.
Incidentally, with respect to the range where the rotating machine is determined to belong for the first reference for unbalanced data U _b, the basic influence coefficient alpha _{1, b} corresponding to the first unbalance data U _b, alpha _{2, b} is step S9 If the basic influence coefficients α _{1, b} and α _{2, b} used in step S2 are the initial influence coefficients α _{1, b} and α _{2, b} , the second reference is used. The basic influence coefficients α _{1, d} and α _{2, d} corresponding to the second unbalance data U _d have not yet been calculated in step S9 with respect to the range in which the rotating machine is determined to belong to the unbalanced data U _d for use. In this case, the basic influence coefficients α _{1, d} and α _{2, d} used in the current step S2 are the initial influence coefficients α _{1, d} and α _{2, d} .

In step S3, the first cutting target portion 13a of the rotating body 13 is cut based on the cutting position and the cutting amount indicated by the first unbalance data U _b calculated in step S2. This cutting is performed similarly to step S3 of the first embodiment. Note that the cutting position here is on the outer peripheral portion 13a of the tip end portion in the axial direction of the rotating body 13, and is on the alternate long and short dash line L1 in FIG. In FIG. 1B, a circle drawn with a broken line overlapping on the alternate long and short dash line L1 indicates the cutting tool 11a positioned on the alternate long and short dash line L1.
Further, in step S3, the second cutting target portion 13b of the rotating body 13 is cut based on the cutting position and the cutting amount indicated by the second unbalance data U _d calculated in step S2. This cutting is performed similarly to step S3 of the first embodiment. The cutting position here is on the outer peripheral portion 13b of the disk-like member provided on the rotating body 13, and is on the alternate long and short dash line L2 in FIG. In FIG. 1B, a circle drawn with a broken line overlapping on the alternate long and short dash line L2 indicates the cutting tool 11a positioned on the alternate long and short dash line L2.

In step S4, the rotating body 13 of the rotating machine is rotated at the first rotation speed. In this state, the vibration sensor 5 measures the vibration of the support 3 that supports the rotating body 13, while the angle sensor 7 Measure the rotation angle. Thus, based on the vibration and rotation angle measured at the first rotation speed, the calculator 9 generates vibration data a1 _{, m, n} .
Further, in step S4, the rotary body 13 of the rotary machine is rotated at the second rotational speed, and in this state, the vibration sensor 5 measures the vibration of the support body 3 that supports the rotary body 13, while the angle sensor 7 measures the rotation angle. Thus, based on the vibration and rotation angle measured at the second rotation speed, the computing unit 9 generates vibration data a2 _{, m, n} .
Here, the subscript n of the vibration data a indicates the number of times step S3 has already been performed in the flowchart of FIG. 3 when the current step S4 is performed for the m-th rotating machine (hereinafter the same), and vibration data. Subscripts 1, 2, and m of a are the same as described above.

In step S5, the first and second unbalanced data U _b and U _d are obtained based on the vibration data a _{1, m, n} , a _{2, m, n} generated in the immediately preceding step S4 and the basic influence coefficient. calculate. The first and second unbalanced data U _b and U _d are expressed by the following determinant of [Equation 5].

If step S9 described later has not yet been performed at the time of performing step S5, the basic influence coefficient used in step S5 is the initial influence coefficient.
If step S9, which will be described later, has already been performed at the time of performing step S5, the basic influence coefficient used in step S5 is the latest influence coefficient updated in the previous step S9.

  Further, when a new influence coefficient is calculated for each range of the unbalance amount in step S9 described later at the time of performing step S5 this time, the following is performed.

First, in step S5, based on the initial influence coefficient and the vibration data a1 _{, m, n} , a2 _{, m, n} generated in the immediately preceding step S4, the first and second _values are obtained by the following [Equation 6]. Calculate reference unbalance data.

In this [Equation 6], α _{1, b} , α _{1, d} , α _{2, b} , α _{2, d} are initial influence coefficients, and a _{1, m, n} , a _{2, m, n} are the immediately preceding coefficients. The vibration data generated in step S4 in the above, U _b is the first reference unbalance data, and U _d is the second reference unbalance data.

Next, the range including the unbalance amount based on the first and second reference unbalance data U _b and U _d is determined.
When following (C1) described later in step S9, the sum or average of the unbalance amount indicated by the first reference unbalance data U _b and the unbalance amount indicated by the second reference unbalance data U _d It is determined whether the value belongs to one of the preset ranges, and it is assumed that the currently targeted rotating machine belongs to the range to which the sum or average value belongs.
If according to later in step S9 (C2), and unbalance amount indicated by the first reference for unbalanced data U _b, the larger of the unbalanced amount indicated by the second reference for unbalanced data U _d It is determined which of the plurality of ranges set in advance falls within the range, and the currently targeted rotating machine belongs to the range to which the larger value belongs.
If according to the described later in step S9 (C3), for the first reference for unbalanced data U _b, amount of unbalance the first reference for unbalanced data U _b is indicated, the first unbalance data It is determined which of the plurality of ranges set in advance with respect to U _b belongs, and the rotating machine currently targeted for the range to which the unbalance amount indicated by the first reference unbalance data U _b belongs is determined. As for the second reference unbalance data U _b , the unbalance amount indicated by the second reference unbalance data U _d is a plurality of the ranges set in advance for the second unbalance data U _d. of determines fall within any of the ranges, the current in a range in which the unbalance amount indicated by the second reference for unbalanced data U _d belongs And rotating machinery belongs you are elephants.

When following (C1) or (C2) described later in step S9, if a new influence coefficient has already been calculated in the previous step S9 for the range determined to belong to the rotating machine, this step is performed. In S5, using this new influence coefficient as a basic influence coefficient, the basic influence coefficient and the vibration data a _{1, m, n} , a _{2, m, n} generated in the immediately preceding step S4 are _expressed in [Equation 5]. By applying, the first and second unbalanced data U _b and U _d are calculated. Otherwise, the basic influence coefficient used in step S5 of this time is the initial influence coefficient.

On the other hand, when according to later in step S9 (C3), the first reference for unbalanced data U _b, with respect to the range where the rotating machine is determined to belong, as described above, the first unbalance data U _{b 1} alpha new influence coefficients corresponding _{to, b,} when alpha _{2, b} is already calculated in the previous step S9 is that in this step S5, the new influence coefficients alpha _{1, b,} alpha _{2, b} is a basic influence coefficient α _{1, b} , α _{2, b} corresponding to the first unbalance data U _b , and the second reference unbalance data U _d is determined to belong to the rotating machine as described above. with respect to said ranges, the second unbalanced data U new influence coefficients corresponding to the _d alpha _{1, d,} when the alpha _{2, d} is already calculated in the previous step S9 is that in this step S5, The new influence coefficients alpha _{1, d,} alpha _2, the basic influence coefficient corresponding to _d in the second unbalance data _{U d α 1, d, α} 2, and _d, the basic influence coefficient α _{1, b,} α _{2, b} , α _{1, d} , α _{2, d} and the vibration data a _{1, m, n} , a _{2, m, n} generated in the immediately preceding step S4 are applied to [Equation 5], First and second unbalanced data U _b and U _d are calculated.
Incidentally, with respect to the range where the rotating machine is determined to belong for the first reference for unbalanced data U _b, the basic influence coefficient alpha _{1, b} corresponding to the first unbalance data U _b, alpha _{2, b} is step S9 In step S5, the basic influence coefficients α _{1, b} and α _{2, b} used in the current step S5 are initial influence coefficients α _{1, b} and α _{2, b,} which are the second reference. The basic influence coefficients α _{1, d} and α _{2, d} corresponding to the second unbalance data U _d have not yet been calculated in step S9 with respect to the range in which the rotating machine is determined to belong to the unbalanced data U _d for use. In this case, the basic influence coefficients α _{1, d} and α _{2, d} used in the current step S5 are the initial influence coefficients α _{1, d} and α _{2, d} .

In step S6, it is determined whether the unbalance amount between the first unbalance data U _b calculated in step S5 and the second unbalance data U _d calculated in step S5 is equal to or less than a threshold value.
For example, the sum or average value of the unbalance amount indicated by the first unbalance data U _b calculated in step S5 and the unbalance amount indicated by the second unbalance data U _d calculated in step S5 is a threshold value. Determine whether:
Alternatively, the larger value of the unbalance amount indicated by the first unbalance data U _b calculated in step S5 and the unbalance amount indicated by the second unbalance data U _d calculated in step S5 is the threshold value. Determine if it is less than or equal to the value.
Alternatively, the unbalance amount indicated by the first unbalance data U _b calculated in step S5 is less than or equal to the first unbalance data threshold value, and the second unbalance data calculated in step S5. It is determined whether the unbalance amount indicated by U _d is equal to or less than the threshold value for the second unbalance data. In this case, preferably, the threshold value for the first unbalanced data is different from the threshold value for the second unbalanced data.

If the determination in step S6 is YES, the rotary machine currently being processed is shipped, and the process proceeds to step S7.
On the other hand, if the determination in step S6 is no, the process proceeds to step S61. In step S61, it is determined whether or not step S3 has already been performed for the set number of times (two or more times) for the rotating machine currently performing the process of FIG. If the determination in step S61 is yes, it is determined that the balance of the rotating machine cannot be corrected, and the process proceeds to step S7. If the determination in step S61 is no, the process returns to step S3. In the returned step S3, the cutting of the first and second cutting target portions 13a and 13b is performed based on the first and second unbalance data U _b and U _d calculated in the immediately preceding step S5, respectively. The other points are the same as described above.

In step S7, the vibration data a _{1, m, 0} , a _{2, m, 0} generated in step S1, the vibration data a _{1, m, n} , a _{2, m, n} generated in step S4, and step S3 A relational expression between the data of the cutting performed in step 1 and the influence coefficient regarding the m-th rotating machine is generated.
This relational expression is expressed by the following determinant of [Equation 7].

The definition of each symbol in [Equation 7] is as follows.
α _{1, b} , α _{1, d} , α _{2, b} , α _{2, d} are influence coefficients related to the m-th rotating machine.
ΔM _{b, m, n} is data of cutting of the first cutting target portion 13a performed in step S3. That is, ΔM _{b, m, n} = A _b (cos θ _b + j sin θ _b ). Here, A _b represents an unbalance amount, and A _b is a product of x and y. x is the radial distance from the central axis C of the rotating body 13 to the cutting position in step S3, and y is the mass cut in step S3. On the other hand, theta _b indicates the circumferential position of the cutting position.
Similarly, ΔM _{d, m, n} is data of cutting of the second cutting target portion 13b performed in step S3. That is, ΔM _{d, m, n} = A _d (cos θ _d + j sin θ _d ). Here, A _d indicates an unbalance amount, and A _d is a product of x and y. x is the radial distance from the central axis C of the rotating body 13 to the cutting position in step S3, and y is the mass cut in step S3. On the other hand, θ _d indicates the circumferential position of the cutting position.
In addition, the subscript b of ΔM in [Expression 7] indicates that the ΔM is cutting data of the first cutting target portion 13a, and the subscript d of ΔM indicates that the ΔM is the second cutting target portion. 13b indicates that the data is cutting data, and the suffix m of ΔM indicates that the ΔM is data relating to the mth rotating machine, and the suffix n of ΔM indicates that the ΔM is the nth time in step S3. (It is the same in the expressions representing e1 _{, m, n} and e2 _{, m, n} described later).
The subscript n of the vibration data a indicates the number of times of step S3 already performed at the time when the vibration data is generated for the m-th rotating machine.
If the m-th rotating machine has not returned from step S6 to step S3, subscripts n and n-1 of vibration data a are 1 and 0, respectively, and a set of relational expressions ( That is, a set of [Equation 7]) with n = 1 is obtained.
Further, in the case of returning from step S6 to step S3 k times, (k + 1) sets of the relational expressions (that is, [Equation 7] with n = 1 and [Equation 7] with n = 2, ... [Expression 7] k relational expressions where n = k, and (k + 1) relational expression [Expression 7] where n = k + 1.

  In step S8, it is determined whether or not the number of rotating machines that have performed steps S1 to S7 is equal to or greater than a set number (for example, a set number of 3 or more). That is, it is determined whether the above m is equal to or greater than the set number. If this determination is YES, the process proceeds to step S9; otherwise, the value of m is incremented by one and the process returns to step S1, and steps S1 to S8 are performed again for a new rotating machine of the same model. For each rotating machine with the number m equal to or greater than the set number, the process proceeds from step S8 to step S9.

  In step S9, a new influence coefficient is calculated based on the plurality of relational expressions obtained by performing steps S1 to S7 for the same number of rotating machines of the number m equal to or greater than the set number. The basic influence coefficient is updated to the new influence coefficient. That is, the next time the step S2 and step S5 are performed, the new influence coefficient is used as the basic influence coefficient. In addition, the number of the relational expressions used in step S9 increases as the number of rotating machines that performed steps S1 to S7 (that is, the value of m) increases.

In step S9, preferably, a new influence coefficient is calculated from a plurality of relational expressions by the least square method with the influence coefficient in each relational expression as an unknown. Specifically, new influence coefficients α _{1, b} , α _{1, d} , α _{2, b} , α _{2, d} are calculated by solving the following [Equation 8].

In [Equation 8], e1 _{, m, n} and e2 _{, m, n} are expressed by the following equations (16) and (17).

_{e 1, m, n = a} 1, m, n -a 1, m, n-1 + α 1, b × ΔM b, m, n + α 1, d × ΔM d, m, n
... (16)

_{e 2, m, n = a} 2, m, n -a 2, m, n-1 + α 2, b × ΔM b, m, n + α 2, d × ΔM d, m, n
... (17)

In [Equation 8], α _{1, b, R} and α _{1, b, I} are a real part and an imaginary part of α _{1, b} , respectively, and α _{1, d, R} and α _{1, d, I} is the real part and imaginary part of α _{1, d} , respectively, α _{2, b, R} and α _{2, b, I} are the real part and imaginary part of α _{2, b} , respectively, and α _{2 , D, R} and α _{2, d, I} are the real and imaginary parts of α _{2, d} , respectively. That is, the following expressions (18) to (21) are established with j as an imaginary unit.

α _{1, b} = α _{1, b, R} + jα _{1, b, I} (18)

α _{1, d} = α _{1, d, R} + jα _{1, d, I} (19)

α _{2, b} = α _{2, b, R} + jα _{2, b, I} (20)

α _{2, d} = α _{2, d, R} + jα _{2, d, I} (21)

In [Equation 8], Σ _{m, n} represents the sum of all combinations of m and n. For this, when the value of m at the present time and m _c, m is an integer from 1 to m _c, for each m, if the back by k times from step S6 to step S3, n is, It is an integer from 1 to (k + 1).

When a new influence coefficient is calculated in step S9, the value of m is incremented by one and the process returns to step S1, and the above-described processing is performed again for a new rotating machine of the same model. However, in the next step S2 and step S5, the first and second unbalanced data U _b and U _d are calculated using the new influence coefficient calculated in the immediately preceding step S9 as the basic influence coefficient.

Preferably, step S9 is performed as follows.
First, a plurality of unbalance amount ranges indicated by the unbalance data are set in advance. For example, three ranges R1, R2, and R3 shown in FIG. 4A are set in advance for the unbalance amount. The range R1 is a range where the unbalance amount is 0 to X1, the range R2 is a range where the unbalance amount is X1 to X2, and the range R3 is a range where the unbalance amount is X2 or more.
Then, for each of the rotating machines, the range of the unbalance amount based on the first and second unbalance data U _b and U _d calculated for the rotating machine in step S2 is determined, and the rotation is included in the range. Assume that the machine belongs. This determination method includes the following methods (C1), (C2), and (C3).

(C1) The sum or average value of the unbalance amount indicated by the first unbalance data U _b calculated in step S2 and the unbalance amount indicated by the second unbalance data U _d calculated in step S2 is calculated in advance. It is determined which of the plurality of set ranges belongs, and the rotating machine belongs to a range to which the sum or average value belongs.
(C2) The larger value of the unbalance amount indicated by the first unbalance data U _b calculated in step S2 and the unbalance amount indicated by the second unbalance data U _d calculated in step S2 is: It is determined which of the plurality of ranges set in advance belongs, and the rotating machine belongs to the range to which the larger value belongs.
(C3) For the first unbalanced data U _b , a plurality of unbalance amount ranges are set in advance, and the unbalance amount indicated by the first unbalance data U _b calculated in step S2 is within the plurality of ranges. It is determined in which range the rotation machine belongs to the range to which the unbalance amount belongs, and a plurality of unbalance amount ranges are set in advance for the second unbalance data U _d. It is determined which of the plurality of ranges the unbalance amount indicated by the calculated second unbalance data U _d belongs, and the rotating machine belongs to the range to which the unbalance amount belongs.

When the method of (C1) or (C2) described above is adopted, in step S9, for each of the ranges, based on the plurality of relational expressions generated in step S7 for each of the rotating machines belonging to the range. Thus, a new influence coefficient is obtained in the same manner as described above.
When the method (C3) described above is adopted, in step S9, for each of the ranges of the first unbalanced data U _b , a plurality of pieces generated in step S7 for each of the rotating machines belonging to the range. As described above, new influence coefficients α _{1, b} , α _{2, b} corresponding to the first unbalanced data U _b are obtained based on the above relational expression, and the second unbalanced data U _d For each of the ranges, a new influence coefficient corresponding to the second unbalanced data U _d based on the plurality of relational expressions generated in step S7 for each of the rotating machines belonging to the range. α _{1, d} and α _{2, d} are obtained.

  As shown in FIG. 4B, since the magnitude (absolute value) of the influence coefficient tends to change as the unbalance amount increases, a plurality of rotations are performed for each range of the unbalance amount as described above. A new influence coefficient reflecting the characteristics of the machine is obtained, and the imbalance data can be calculated with high accuracy using this influence coefficient. In the case of (C3) described above, a new influence coefficient can be obtained in consideration of the balance change of the first and second cutting target portions 13a and 13b.

[Third Embodiment]
FIG. 5 is a flowchart illustrating an influence coefficient acquisition method according to the third embodiment of the present invention. In the influence coefficient acquisition method according to the third embodiment, the following points are changed from the second embodiment described above, and the other points are the same as those of the second embodiment.
According to the third embodiment, as shown in FIG. 5, after step S3 is performed, step S11 is performed before proceeding to step S4.
In step S11, when step S3 has already been performed for the same rotating machine (that is, the rotating machine currently performing the process of FIG. 5) for the set number of times (two or more times), the process proceeds to step S12. If not, the process proceeds to step S4. When the process proceeds to step S4, the subsequent processing is the same as described above except that step S61 described above is omitted. That is, if the determination in step S6 is NO, the process returns to step S3 without performing step S61 after step S6.

In step S12, the rotating body 13 of the rotating machine is rotated at the first rotation speed. In this state, the angle sensor 7 rotates while the vibration sensor 5 measures the vibration of the support 3 that supports the rotating body 13. Measure the corner. Based on the vibration and rotation angle measured in this way, the computing unit 9 generates vibration data a1 _{, m, n} .
Furthermore, in step S12, the rotating body 13 of the rotating machine is rotated at a second rotational speed different from the first rotational speed. In this state, the vibration sensor 5 vibrates the support 3 that supports the rotating body 13. While measuring, the angle sensor 7 measures the rotation angle. Based on the vibration and rotation angle measured in this way, the computing unit 9 generates vibration data a2 _{, m, n} .

In step S13, the first and second unbalanced data U _b and U _d are calculated based on the vibration data a _{1, m, n} , a _{2, m, n} generated in step S12 and the basic influence coefficient. . The first and second unbalanced data U _b and U _d are expressed by the following determinant of [Equation 9].

  If the above-described step S9 has not been performed yet at the time of performing the present step S13, the basic influence coefficient used in the present step S13 is an initial influence coefficient.

  If the above-described step S9 has already been performed at the time of performing step S13 this time, the basic influence coefficient used in this step S13 is the latest influence coefficient updated in the previous step S9.

  Further, when a new influence coefficient is calculated for each range of the unbalance amount at the above-described step S9 at the time of performing step S13 this time, the following is performed.

First, in step S13, based on the initial influence coefficient and the vibration data a1 _{, m, n} , a2 _{, m, n} generated in the immediately preceding step S12, the first and second _values are obtained by the following [Equation 10]. Calculate reference unbalance data.

In this [Equation 10], α _{1, b} , α _{1, d} , α _{2, b} , α _{2, d} are initial influence coefficients, and a _{1, m, n} , a _{2, m, n} are the immediately preceding coefficients. The vibration data generated in step S12 in the above, U _b is the first reference unbalance data, and U _d is the second reference unbalance data.

Next, the range including the unbalance amount based on the first and second reference unbalance data U _b and U _d is determined.
When (C1) is followed in step S9, the sum or average of the unbalance amount indicated by the first reference unbalance data U _b and the unbalance amount indicated by the second reference unbalance data U _d It is determined whether the value belongs to one of the preset ranges, and it is assumed that the currently targeted rotating machine belongs to the range to which the sum or average value belongs.
When (C2) described above is followed in step S9, the larger of the unbalance amount indicated by the first reference unbalance data U _b and the unbalance amount indicated by the second reference unbalance data U _d It is determined which of the plurality of ranges set in advance falls within the range, and the currently targeted rotating machine belongs to the range to which the larger value belongs.
When (C3) described above in step S9 is followed, for the first reference unbalance data U _b , the unbalance amount indicated by the first reference unbalance data U _b is equal to the first unbalance data. It is determined which of the plurality of ranges set in advance with respect to U _b belongs, and the rotating machine currently targeted for the range to which the unbalance amount indicated by the first reference unbalance data U _b belongs is determined. and it belongs, for the second reference for unbalanced data U _d, unbalance amount which the second reference for unbalanced data U _d is indicated, a plurality of the preset range for the second unbalance data U _d of determines fall within any of the ranges, the current in a range in which the unbalance amount indicated by the second reference for unbalanced data U _d belongs And rotating machinery belongs you are elephants.

When following (C1) or (C2) described later in step S9, if a new influence coefficient has already been calculated in the previous step S9 for the range determined to belong to the rotating machine, this step is performed. In S13, using this new influence coefficient as a basic influence coefficient, the basic influence coefficient and the vibration data a _{1, m, n} , a _{2, m, n} generated in the immediately preceding step S12 are applied to [Equation 9]. As a result, the first and second unbalanced data U _b and U _d are calculated. Otherwise, the basic influence coefficient used in step S13 this time is the initial influence coefficient.

On the other hand, when according to the aforementioned (C3) in step S9, the first reference for unbalanced data U _b, with respect to the range where the rotating machine is determined to belong, as described above, the first unbalance data U _{b 1} alpha new influence coefficients corresponding _{to, b,} when alpha _{2, b} is already calculated in the previous step S9 is that in this step S13, the new influence coefficients alpha _{1, b,} alpha _{2, b} is a basic influence coefficient α _{1, b} , α _{2, b} corresponding to the first unbalance data U _b , and the second reference unbalance data U _d is determined to belong to the rotating machine as described above. with respect to said ranges, the second unbalanced data U new influence coefficients corresponding to the _d alpha _{1, d,} when the alpha _{2, d} is already calculated in the previous step S9 is in this step S13 , The new influence coefficients alpha _{1, d,} alpha _2, the basic influence coefficient corresponding to _d in the second unbalance data _{U d} alpha _{1, d,} alpha _2, and _d, the basic influence coefficient alpha _1, b, By applying α _{2, b} , α _{1, d} , α _{2, d} and the vibration data a _{1, m, n} , a _{2, m, n} generated in the immediately preceding step S12 to [Equation 9] First and second unbalanced data U _b and U _d are calculated.
Incidentally, with respect to the range where the rotating machine is determined to belong for the first reference for unbalanced data U _b, the basic influence coefficient alpha _{1, b} corresponding to the first unbalance data U _b, alpha _{2, b} is step S9 In step S13, the basic influence coefficients α _{1, b} and α _{2, b} used in the current step S13 are initial influence coefficients α _{1, b} and α _{2, b.} The basic influence coefficients α _{1, d} and α _{2, d} corresponding to the second unbalance data U _d have not yet been calculated in step S9 with respect to the range in which the rotating machine is determined to belong to the unbalanced data U _d for use. In this case, the basic influence coefficients α _{1, d} and α _{2, d} used in step S13 are initial influence coefficients α _{1, d} and α _{2, d} .

In step S14, it is determined whether or not the amount of unbalance between the first unbalance data U _b calculated in step S13 and the second unbalance data U _d calculated in step S13 is equal to or less than a threshold value.
For example, the sum or average value of the unbalance amount indicated by the first unbalance data U _b calculated in step S13 and the unbalance amount indicated by the second unbalance data U _d calculated in step S13 is the threshold value. Determine whether:
Alternatively, the larger value of the unbalance amount indicated by the first unbalance data U _b calculated in step S13 and the unbalance amount indicated by the second unbalance data U _d calculated in step S13 is the threshold value. Determine if it is less than or equal to the value.
Alternatively, the unbalance amount indicated by the first unbalance data U _b calculated in step S13 is less than or equal to the first unbalance data threshold value, and the second unbalance data calculated in step S13. It is determined whether the unbalance amount indicated by U _d is equal to or less than the threshold value for the second unbalance data. In this case, preferably, the threshold value for the first unbalanced data is different from the threshold value for the second unbalanced data.
If the determination in step S14 is YES, the rotary machine is shipped, and the method of the present embodiment is performed from step S1 on the next rotary machine. That is, the value of m indicating the number of the rotating machine is increased by 1, and step S1 is performed.
On the other hand, if the determination in step S14 is no, the process proceeds to step S15.

In step S15, it is determined whether or not step S14 is performed more than the allowable number of times for the current rotating machine. That is, it is determined whether the number of times determined as NO in step S14 exceeds the allowable number. If the determination in step S15 is YES, assuming that the balance of the rotating machine cannot be corrected, the process of FIG. 5 is terminated for the rotating machine, and the method of the present embodiment from step S1 for the next rotating machine. To implement. That is, the value of m indicating the number of the rotating machine is increased by 1, and step S1 is performed.
On the other hand, if the determination in step S15 is no, the process proceeds to step S16.

In step S < _b > 16, the computing unit 9 calculates the replacement influence coefficients α _{1, b} , α _{1, d} , α _{2, b} , α _{2, d} for the rotating body 13 of the m-th rotating machine by the following equation (24). Calculated by (27). That is, the influence coefficient for replacement is calculated by solving the following equations (24) to (27) simultaneously.

a1 _{, m, p-} a1 _{, m, q} = [alpha] _{1, b} * (-[Delta] _Mbpq ) + [alpha] _{1, d} * (-[Delta] _Mdpq ) (24)

_{a 1, m, r -a 1} , m, s = α 1, b × (-ΔM brs) + α 1, d × (-ΔM drs) ··· (25)

a _{2, m, p} −a _{2, m, q} = α _{2, b} × (−ΔM _bpq ) + α _{2, d} × (−ΔM _dpq ) (26)

_{a 2, m, r -a 2} , m, s = α 2, b × (-ΔM brs) + α 2, d × (-ΔM drs) ··· (27)

In the above formulas (24) to (27), the subscript p, q, r or s of the vibration data a corresponds to the subscript n described above, and the vibration data is the p, q, r or s step. The vibration data generated immediately after S3 (in step S12 or S4) is indicated. That is, the subscript p, q, r, or s of the vibration data a is the number of times step S3 has already been performed when the vibration data is generated. The subscript m of the vibration data a is the same as described above.
In the above equations (24) to (27), p, q, r, and s are values equal to or less than the set number used in step S11. p is greater than q and r is greater than s. Further, at least one of p ≠ r and q ≠ s holds. If the subscript q or s of the vibration data a is 0, the vibration data is the vibration data generated in step S1. Preferably, the set number of times is 2, p is 2, r is 1, and q and s are 0.

ΔM _bpq , ΔM _dpq Δ, M _brs , and ΔM _drs in the above equations (24) to (27) are the data of the cutting performed in step S3 described above, and are generated by, for example, the position controller 11c and calculated by the calculator 9. Is input. ΔM _bpq , ΔM _dpq Δ, M _brs , and ΔM _drs are specifically as follows.

In the above equations (24) and (26), ΔM _bpq is represented by the following equation (28).

ΔM _bpq = ΣA _bi (cos θ _bi + jsin θ _bi ) (28)

Here, j is an imaginary unit, b indicates cutting data in the first cutting target portion 13a, Σ indicates a sum related to i, and i is an integer greater than q and less than or equal to p. It is. Further, in the equation (28), A _bi is the mass obtained by cutting the first cutting target portion 13a in the i-th step S3, the position where the cutting is performed in the first cutting target portion 13a, and the rotational body 13 It is a product of the distance from the central axis C, and θ _bi is a phase indicating the circumferential position where the first cutting target portion 13a is cut in the i-th step S3. However, since i = 0 indicates that step S3 has not been performed, when i = 0, A _bi = 0.

In the above equations (24) and (26), ΔM _dpq is represented by the following equation (29).

ΔM _dpq = ΣA _di (cos θ _di + j sin θ _di ) (29)

Here, j is an imaginary unit, d indicates cutting data in the second cutting target portion 13b, Σ indicates a sum related to i, and i is an integer greater than q and less than or equal to p. It is. In Formula (29), A _di is the mass of cutting the second cutting target portion 13b in the i-th step S3, the position where the cutting is performed in the second cutting target portion 13b, and the position of the rotating body 13. It is a product of the distance from the central axis C, and θ _di is a phase indicating a circumferential position where the second cutting target portion 13b is cut in the i-th step S3. However, since i = 0 indicates that step S3 has not been performed, when i = 0, A _di = 0.

In the above equations (25) and (27), ΔM _brs is represented by the following equation (30).

ΔM _brs = ΣA _bi (cos θ _bi + jsin θ _bi ) (30)

Here, j is an imaginary unit, b indicates cutting data in the first cutting target portion 13a, Σ indicates a sum related to i, and i is an integer greater than s and equal to or less than r. It is. Further, in the equation (30), A _bi is the mass of cutting the first cutting target portion 13a in the i-th step S3, the position where the cutting is performed in the first cutting target portion 13a, and the rotational body 13. It is a product of the distance from the central axis C, and θ _bi is a phase indicating the circumferential position where the first cutting target portion 13a is cut in the i-th step S3. However, since i = 0 indicates that step S3 has not been performed, when i = 0, A _bi = 0.

In the above equations (25) and (27), ΔM _drs is represented by the following equation (31).

ΔM _drs = ΣA _di (cos θ _di + j sin θ _di ) (31)

Here, j is an imaginary unit, d indicates cutting data in the first cutting target portion 13a, Σ indicates a sum relating to i, and i is an integer greater than s and less than or equal to r. It is. In Formula (31), A _di is the mass of cutting the second cutting target portion 13b in the i-th step S3, the position where the cutting is performed in the second cutting target portion 13b, and the position of the rotating body 13. It is a product of the distance from the central axis C, and θ _di is a phase indicating a circumferential position where the second cutting target portion 13b is cut in the i-th step S3. However, since i = 0 indicates that step S3 has not been performed, when i = 0, A _di = 0.

In step S17, the vibration data a _{1, m, n} , a _{2, m, n} generated in step S12 and the influence coefficients α _{1, b} , α _{1, d} , α _{2, b} for replacement calculated in step S16 are displayed. , Α _{2, d} , the first and second unbalanced data U _b and U _d are calculated by the following [Equation 11].

Step S17 returns to step S3 After making, in step S3 of returning, based on the cutting position and cutting amount indicated by the first unbalanced data U _b calculated in step S17, the first cutting target portion of the rotating body 13 13a with cutting based on the cutting amount and the cutting position indicated by the second unbalance data U _d calculated in step S17, cutting the second cut target portion 13b of the rotating body 13. The cutting of each of the cutting target portions 13a and 13b here is performed in the same manner as Step S3 of the first embodiment described above.
The procedure after cutting in the returned step S3 is the same as described above.

[Example]
FIG. 6A shows a case where the influence coefficient acquisition device 10 according to the above-described embodiment is applied to a supercharger 20 as a rotating machine.

  As shown in FIG. 6A, the rotating body 13 of the supercharger 20 supplies the compressed air to the engine by rotating integrally with the turbine blade 15 that is rotationally driven by the exhaust gas of the engine and the turbine blade 15. And a rotating shaft 19 having a turbine blade 15 coupled to one end and a compressor blade 17 coupled to the other end. Moreover, the supercharger 20 has the stationary side member 21 which supports the rotary body 13 rotatably. In the example of FIG. 6A, the stationary side member 21 is a bearing housing in which bearings 23a and 23b that rotatably support the rotating body 13 (rotating shaft 19) are incorporated. Further, the supercharger 20 includes a turbine housing 25 that houses the turbine blades 15 therein, and a compressor housing (removed in FIG. 6A) that houses the compressor blades 17 inside. The turbine housing 25 is formed with a flow path (scroll) through which a fluid for rotationally driving the turbine blades 15 flows. The turbine housing 25 is attached to the inside of the support body 3. The support 3 is configured so that the fluid that drives the turbine blade 15 can be supplied to the flow path of the turbine housing 25, and the fluid that has driven the turbine blade 15 can be discharged to the outside of the support 3.

  Further, the support 3 supports the bearing housing 21 via the turbine housing 25 or directly. For example, the support body 3 may directly support the bearing housing 21 by attaching the bearing housing 21 to the support body 3 by an appropriate means. Thereby, the support body 3 supports the rotary body 13 through the bearing housing 21 so that rotation is possible.

  In the example of FIG. 6 (A), the turbine housing 25 is different from the turbine housing of the finished supercharger 20 and is dedicated for carrying out the influence coefficient acquisition method according to the first, second, or third embodiment described above. It is a housing. That is, if YES is obtained in step S6 or step S14, a dedicated turbine housing and a compressor housing are attached to the finished turbocharger 20 in the bearing housing 21 of the rotating body 13, whereby the supercharger 20 is completed. To do.

  FIG. 6B is a view taken along the line BB in FIG. 6A, but the illustration of the bearing housing 21 and the like is omitted. The first cutting target portion 13a is an outer peripheral portion of the nut 27 screwed into the axial end portion on the compressor blade 17 side of the rotating body 13 (that is, an annular portion drawn by a one-dot chain line L1 in FIG. 6B). is there. The compressor impeller 12 in which a plurality of compressor blades 17 are formed by the screwing of the nut 27 is fastened to the rotary shaft 19. The compressor impeller 12 has a blade coupling portion 24 extending in the circumferential direction so as to couple the plurality of compressor blades 17 in the circumferential direction. The second cutting target portion 13b is an outer peripheral portion of the blade coupling portion 24 (that is, an annular portion drawn by a broken line L2 in FIG. 6B). The plurality of compressor blades 17 are arranged at intervals in the circumferential direction.

In the case of cutting the second cut target portion 13b on the basis of the second unbalance data U _d, a second unbalance data U _d indicate circumferential position (corresponding to the argument of the complex number U _d) is a compressor When it is at the position of the wing 17, the following is performed. The unbalance data U _d, the outer peripheral portion of the blade coupling part 24 in the (second cutting target portion 13b), a plurality of locations within (hatched portion in FIG. 6 (B)) ranges compressor blade 17 is not positioned To an unbalance amount (hereinafter referred to as a decomposition unbalance amount). Then, at each of the plurality of positions, the second cutting target portion 13b is cut by an amount corresponding to the corresponding disassembly imbalance amount, as described above.
Since the unbalanced data U _d can be regarded as a vector composed of the circumferential position where the unbalance exists (corresponding to the declination of the complex number U _d ) and the unbalanced amount (corresponding to the absolute value of the complex number U _d ), the unbalanced data U _d U _d can be decomposed into a plurality of vector components each consisting of a circumferential position (vector direction) where an imbalance exists and an unbalance amount (vector magnitude).

  In FIG. 6A, the angle sensor 7 and the cutting device 11 are both provided on the compressor blade 17 side. However, when the angle sensor 7 is used, the cutting tool 11a of the cutting device 11 does not interfere with the angle sensor 7. When retracting to a position and using the cutting device 11, the angle sensor 7 may be retracted to a position where the angle sensor 7 does not interfere with the cutting device 11.

  In FIG. 6, other configurations and operations of the influence coefficient acquisition apparatus 10 may be the same as those in the above-described embodiment. Moreover, you may implement the influence coefficient acquisition method by the above-mentioned embodiment with respect to the supercharger 20 using the influence coefficient acquisition apparatus 10 shown in FIG.

  The present invention is not limited to the above-described embodiment, and various changes can be made without departing from the scope of the present invention. For example, the following modification examples 1, 2, and 3 may be arbitrarily combined and employed, or any one of modification examples 1, 2, and 3 may be employed. In this case, other points may be the same as described above, or may be changed as appropriate.

(Modification 1)
For example, in the second embodiment and the third embodiment described above, vibration and rotation angle are measured in a state where the rotating body 13 is rotated at the first rotation speed, and vibration data is generated. Vibration and rotation angle were measured in a state where the rotating body 13 was rotated at the speed, and vibration data was generated. Instead, as shown in FIG. Two vibration sensors 5a and 5b may be attached. In this case, it is performed as follows.
In step S1, the rotating body 13 of the rotating machine is rotated, and in this state, the vibration of the support body that supports the rotating body 13 is measured by the first and second vibration sensors 5a and 5b, and the angle sensor 7 is used. The rotation angle of the rotator 13 is measured, and based on the vibration measured by the first vibration sensor 5a and the rotation angle measured by the angle sensor 7, the calculator 9 generates vibration data a1 _{, m, n} (however, in step S1) In this case, n = 0, and the vibration data a _{2, m, n} (however, the step 9) is calculated by the calculator 9 based on the vibration measured by the second vibration sensor 5b and the rotation angle measured by the angle sensor 7. In S1, n = 0) is generated. The same applies to step S4 and step S12. In this case, the subscript 1 of the vibration data a indicates that the vibration data is generated based on the vibration measured by the first vibration sensor 5a, and the subscript 2 of the vibration data a indicates that the vibration data is the first one. 2 is generated based on the vibration measured by the second vibration sensor 5b, and other suffixes of the vibration data a are the same as described above. Vibration data _{a 1,} m of the modification 1, _{n, a 2, m, n,} respectively, the vibration data _{a 1} in the second embodiment or the third embodiment described _above, m, _{n, a 2, m,} is replaced by _n .

In the first modification, when step S9 is not yet performed when step S2 is performed, the basic influence coefficient used in step S2 at this time is the rotating machine that is currently performing the processing of FIG. 3 or FIG. May be the initial influence coefficients α _{1, b} , α _{1, d} , α _{2, b} , α _{2, d} acquired in advance by the above [Equation 3].
However, in this case, in Equation 3 above, a ₁₀ is not attached trial weight to the rotating body 13, as described above, the first vibration sensor 5a is measured while rotating the rotary member 13 The vibration data generated based on the vibration is a ₂₀ measured by the second vibration sensor 5b with the rotating body 13 rotated as described above without attaching the test weight to the rotating body 13. The vibration data generated based on the vibration, a _1b is measured by the first vibration sensor 5a in a state where the trial weight is attached to the first cutting target portion 13a and the rotating body 13 is rotated as described above. A _1d is a vibration data generated on the basis of the vibration, and the first vibration sensor 5a is attached to the second cutting target portion 13b and the rotating body 13 is rotated as described above. Generated based on measured vibration This is vibration data, and a _2b is generated based on the vibration measured by the second vibration sensor 5b with the trial weight attached to the first cutting target portion 13a and the rotating body 13 rotated as described above. A _2d is based on the vibration measured by the second vibration sensor 5b in a state where the trial weight is attached to the second cutting target portion 13b and the rotating body 13 is rotated as described above. It is the generated vibration data.

(Modification 2)
In the above description, the cutting of the portion to be cut is performed in the axial direction, but may be performed in the radial direction.

(Modification 3)
In the above description, the initial influence coefficient is obtained by giving a balance change to the rotating body 13 using a trial weight. However, by cutting the rotating body 13, the initial influence coefficient is obtained by giving a balance change to the rotating body 13. You may get it.

3 Support body, 5 Vibration sensor, 7 Angle sensor, 9 Calculator, 11 Cutting device

Claims (6)

An influence coefficient acquisition method for acquiring an influence coefficient indicating a vibration change of a rotating body with respect to a balance change of the rotating body provided in the rotating machine,
(A) The rotating body of the rotating machine is rotated, and in this state, the vibration of the support body that supports the rotating body is measured, the rotation angle of the rotating body is measured, and the vibration is based on the measured vibration and the rotation angle. Generate data,
(B) Calculate unbalance data based on the vibration data generated in (A) and the basic influence coefficient,
(C) cutting the rotating body based on the unbalance data;
(D) Thereafter, the rotating body of the rotating machine is rotated, and in this state, the vibration of the support that supports the rotating body is measured, the rotation angle of the rotating body is measured, and the measured vibration and rotation angle are used. Vibration data,
(E) Vibration data generated in (A) and (D), data of cutting performed in (C), and influence on the rotating machine performed in (A), (C), and (D) Get the relational expression with the coefficient,
(F) The above-mentioned (A) to (E) are performed on a plurality of rotating machines of the same model, and a new influence on the model is obtained on the basis of the plurality of relational expressions obtained in (E) above for these rotating machines. An influence coefficient acquisition method characterized by obtaining a coefficient.   The said (A)-(E) is again performed about the new rotary machine same as the said model after calculating | requiring the said new influence coefficient as the said new influence coefficient as the said basic influence coefficient. 2. The influence coefficient acquisition method according to 1. A plurality of unbalance amount ranges indicated by the unbalance data are set in advance,
For each rotating machine, determine the range including the unbalance amount based on the unbalance data calculated in (B) for the rotating machine, and assume that the rotating machine belongs to the range.
In (F), a new influence coefficient is obtained for each of the ranges based on the plurality of relational expressions generated in (E) for the rotating machines belonging to the range. Item 3. An influence coefficient acquisition method according to Item 1 or 2. In (A) above,
(A-1) The rotating body of the rotating machine is rotated at the first rotation speed, and in this state, the vibration of the supporting body that supports the rotating body is measured, and the rotation angle of the rotating body is measured and measured. Generate vibration data based on vibration and rotation angle,
(A-2) Next, the rotating body of the rotating machine is rotated at the second rotation speed, and in this state, the vibration of the supporting body that supports the rotating body is measured, and the rotation angle of the rotating body is measured. Generate vibration data based on the measured vibration and rotation angle,
In (B), first and second unbalance data are calculated based on the vibration data generated in (A-1) and (A-2) and the basic influence coefficient,
In (C) above,
(C-1) Based on the first unbalance data, the first cutting target portion in the rotating body is cut,
(C-2) Based on the second unbalance data, the second cutting target portion in the rotating body is cut,
In (D) above,
(D-1) The rotating body of the rotating machine is rotated at the first rotation speed, and in this state, the vibration of the supporting body that supports the rotating body is measured, and the rotation angle of the rotating body is measured and measured. Generate vibration data based on vibration and rotation angle,
(D-2) Next, the rotating body of the rotating machine is rotated at the second rotation speed, and in this state, the vibration of the supporting body that supports the rotating body is measured, and the rotation angle of the rotating body is measured. Generate vibration data based on the measured vibration and rotation angle,
In (E), the vibration data generated in (A-1), (A-2), (D-1), and (D-2), and (C-1) and (C-2) The relational expression between the data of the performed cutting and the influence coefficient regarding the rotating machine that performed the (A), (C), and (D) is generated. How to get the influence coefficient. The first and second vibration sensors are respectively attached to different places on the support,
In (A), the rotating body of the rotating machine is rotated, and in this state, the vibration of the supporting body that supports the rotating body is measured by the first and second vibration sensors, and the rotation angle of the rotating body is set to an angle. Vibration data is generated based on the vibration measured by the sensor, the vibration measured by the first vibration sensor and the rotation angle measured by the angle sensor, and based on the vibration measured by the second vibration sensor and the rotation angle measured by the angle sensor. Vibration data,
In (B), vibration data generated based on the vibration measured by the first vibration sensor in (A) and vibration data generated based on the vibration measured by the second vibration sensor in (A). And the first and second imbalance data based on the basic influence coefficient,
In (C) above,
(C-1) Based on the first unbalance data, the first cutting target portion in the rotating body is cut,
(C-2) Based on the second unbalance data, the rotary body is cut at the second cutting target part,
In (D), the rotating body of the rotating machine is rotated, and in this state, the vibration of the supporting body that supports the rotating body is measured by the first and second vibration sensors, and the rotation angle of the rotating body is set to the angle. Vibration data is generated based on the vibration measured by the sensor, the vibration measured by the first vibration sensor and the rotation angle measured by the angle sensor, and based on the vibration measured by the second vibration sensor and the rotation angle measured by the angle sensor. Vibration data,
In (E), vibration data generated based on the vibration measured by the first vibration sensor in (A) and vibration data generated based on the vibration measured by the second vibration sensor in (A). Vibration data generated based on the vibration measured by the first vibration sensor in (D), vibration data generated based on the vibration measured by the second vibration sensor in (D), 4. The influence coefficient acquisition method according to claim 1, wherein a relational expression with an influence coefficient relating to the rotating machine on which A), (C), and (D) have been performed is generated. About the rotating machine which performed said (A)-(D), after said (D),
(A) calculating unbalance data based on the vibration data generated in (D) and the basic influence coefficient;
(B) When the unbalance amount indicated by the unbalance data is larger than the threshold value, the rotating body of the rotating machine is cut based on the unbalance data,
(C) Thereafter, for the rotating machine, until the unbalance amount becomes equal to or less than the threshold value, generation of vibration data, calculation of unbalance data, and cutting of the rotating body based on the unbalance data are repeated. Then, after cutting the rotating body a set number of times, based on the data of the cutting performed on the rotating machine and the vibration data generated on the rotating machine, calculate the influence coefficient for replacement of the rotating machine,
(D) calculating unbalanced data based on the influence coefficient for replacement and the vibration data generated last for the rotating machine;
(E) The influence coefficient acquisition method according to claim 4 or 5, wherein the rotating body of the rotating machine is cut based on the unbalanced data.