WO2023162386A1 - 回転機械のフランジ変位量推定方法、この方法を実行するためのプログラム、及び、この方法を実行する装置 - Google Patents

回転機械のフランジ変位量推定方法、この方法を実行するためのプログラム、及び、この方法を実行する装置 Download PDF

Info

Publication number
WO2023162386A1
WO2023162386A1 PCT/JP2022/044204 JP2022044204W WO2023162386A1 WO 2023162386 A1 WO2023162386 A1 WO 2023162386A1 JP 2022044204 W JP2022044204 W JP 2022044204W WO 2023162386 A1 WO2023162386 A1 WO 2023162386A1
Authority
WO
WIPO (PCT)
Prior art keywords
coordinate data
dimensional coordinate
effective
midpoint
positions
Prior art date
Application number
PCT/JP2022/044204
Other languages
English (en)
French (fr)
Inventor
理 熊谷
俊介 水見
光司 石橋
Original Assignee
三菱パワー株式会社
三菱重工業株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 三菱パワー株式会社, 三菱重工業株式会社 filed Critical 三菱パワー株式会社
Priority to JP2024502837A priority Critical patent/JPWO2023162386A1/ja
Priority to KR1020237043921A priority patent/KR20240010039A/ko
Priority to CN202280045297.6A priority patent/CN117561368A/zh
Publication of WO2023162386A1 publication Critical patent/WO2023162386A1/ja

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01DNON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D25/00Component parts, details, or accessories, not provided for in, or of interest apart from, other groups
    • F01D25/24Casings; Casing parts, e.g. diaphragms, casing fastenings
    • F01D25/243Flange connections; Bolting arrangements
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01DNON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D21/00Shutting-down of machines or engines, e.g. in emergency; Regulating, controlling, or safety means not otherwise provided for
    • F01D21/003Arrangements for testing or measuring
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01DNON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D25/00Component parts, details, or accessories, not provided for in, or of interest apart from, other groups
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01DNON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D25/00Component parts, details, or accessories, not provided for in, or of interest apart from, other groups
    • F01D25/24Casings; Casing parts, e.g. diaphragms, casing fastenings
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/02Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
    • G01B21/04Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness by measuring coordinates of points
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05DINDEXING SCHEME FOR ASPECTS RELATING TO NON-POSITIVE-DISPLACEMENT MACHINES OR ENGINES, GAS-TURBINES OR JET-PROPULSION PLANTS
    • F05D2220/00Application
    • F05D2220/30Application in turbines
    • F05D2220/31Application in turbines in steam turbines
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05DINDEXING SCHEME FOR ASPECTS RELATING TO NON-POSITIVE-DISPLACEMENT MACHINES OR ENGINES, GAS-TURBINES OR JET-PROPULSION PLANTS
    • F05D2260/00Function
    • F05D2260/30Retaining components in desired mutual position
    • F05D2260/31Retaining bolts or nuts

Definitions

  • the present disclosure provides a flange displacement estimation method, a program for executing this method, and a method for estimating the displacement on the flange surfaces of the upper half casing and the lower half casing that cover the outer circumference of the rotor in a rotary machine.
  • This application claims priority based on Japanese Patent Application No. 2022-027442 filed in Japan on February 25, 2022, the content of which is incorporated herein.
  • a rotating machine such as a steam turbine includes a rotor rotatable about a horizontally extending axis, a casing that covers the outer circumference of the rotor, stationary parts such as a diaphragm that is placed in the casing and attached to the casing, Prepare.
  • the casing generally has an upper half casing, a lower lower casing half, and a plurality of bolts that fasten the upper half casing and the lower half casing.
  • the upper half casing has an upper flange with an upper flange surface facing downward.
  • the lower half casing has a lower flange formed with a lower flange surface facing upward and vertically opposed to the upper flange surface.
  • the upper half casing When inspecting the rotating machine, the upper half casing is removed from the lower half casing, and the multiple parts that make up the rotating machine are inspected and, if necessary, repaired.
  • a casing in a rotating machine such as a steam turbine may undergo inelastic deformation such as creep deformation under the influence of heat during operation. Therefore, strictly speaking, the lower half casing and the upper half casing in the open state after being operated once are deformed from the time of shipment from the factory.
  • the multiple parts After the inspection is completed, the multiple parts are assembled.
  • This assembling process includes a process of fastening the upper half casing to the lower half casing using a plurality of bolts. In the process of changing the lower half casing and the upper half casing from the open state to the fastened state, the lower half casing and the upper half casing are further deformed.
  • the radial distance between the stationary parts attached to the casing and the rotor must be within the predetermined allowable dimensions. However, when the casing changes from the open state to the closed state and the shapes of the lower half casing and the upper half casing change, the radial distance between the stationary parts attached to the casing and the rotor changes. may fall outside the range of allowable dimensions.
  • the deformation amounts of the lower half casing and the upper half casing when changing from the open state to the fastened state are estimated by the following procedure.
  • a finite element model for the three-dimensional shapes of the lower half casing and the upper half casing is obtained.
  • three-dimensional shape data of the lower half casing and the upper half casing in the open state are acquired by actual measurement.
  • the finite element model is corrected using the measured 3D shape data so that the finite element model matches the measured 3D shape data.
  • the fastened state is simulated using the corrected finite element model indicating the open state to create a finite element model indicating the fastened state.
  • the amount of deformation of the predetermined portion of the lower half casing and the upper half casing is estimated.
  • the predetermined portions of the lower half casing and the upper half casing are the lower flange surface of the lower half casing and the upper flange surface of the upper half casing.
  • the fastening state is simulated using a finite element model that indicates the open state, and the finite element model that indicates the fastening state obtained by this simulation is used to calculate the lower flange surface of the lower half casing. and the amount of displacement of the upper flange surface of the upper half casing.
  • Patent Document 1 uses a finite element model that indicates the open state to simulate the fastening state, so there is a problem that the calculation load for executing this simulation is large. For this reason, the technique described in Patent Document 1 has inherent problems that the preparation period is prolonged and the cost for estimating the displacement amount of the flange surface is increased.
  • the present disclosure can reduce the calculation load, thereby shortening the preparation period for estimating the flange surfaces and reducing the estimation cost.
  • the purpose is to provide technology.
  • This rotary machine includes a rotor rotatable about a horizontally extending axis, a casing covering the outer circumference of the rotor, stationary parts arranged in the casing and attached to the casing, and the casing being mounted downward. and a pedestal supporting from the side.
  • the casing has an upper half casing, a lower lower half casing, and a plurality of bolts that fasten the upper half casing and the lower half casing.
  • the upper half casing has an upper flange formed with an upper flange surface facing downward.
  • the lower half casing includes a lower flange formed with a lower flange surface facing upward and facing the upper flange surface in the vertical direction; a first supported portion and a second supported portion separated from each other in an extending axial direction;
  • the upper flange and the lower flange are formed with bolt holes penetrating in the vertical direction and through which the plurality of bolts can be respectively inserted.
  • the measured three-dimensional coordinate data at a plurality of positions on the upper flange surface and the lower A measured coordinate receiving step of receiving measured three-dimensional coordinate data at a plurality of positions on the flange surface; and a lower first position and a lower second position using the measured three-dimensional coordinate data at the plurality of positions on the lower flange surface.
  • the effective three-dimensional coordinate data at the lower target position and the lower target midpoint position and using the actually measured three-dimensional coordinate data at a plurality of positions in the upper flange surface, the upper first position and the upper an effective coordinate grasping step of grasping effective three-dimensional coordinate data at a second position, an upper target position, and an upper target midpoint position; and an effective three-dimensional coordinate of said lower first position grasped in said effective coordinate grasping step.
  • the data match the effective three-dimensional coordinate data of the first upper position, and the effective three-dimensional coordinate data of the second lower position and the effective three-dimensional coordinate data of the second upper position grasped in the effective coordinate grasping step.
  • the lower first position is a position where the first representative position of the first supported portion coincides with the position in the horizontal direction in the plane continuous with the upper flange surface.
  • the lower second position is a position that coincides with the second representative position of the second supported portion in the horizontal direction in the plane that continues to the lower flange surface.
  • the lower symmetrical position is a position on the lower flange surface where it is desired to obtain an amount of displacement in the vertical direction when changing from the open state to the fastened state.
  • the lower symmetrical midpoint position is the position of the midpoint in the lateral direction, which is the horizontal direction and the direction perpendicular to the axial direction, in the lower flange surface, and the position in the axial direction is the same as the lower symmetrical position. It is the position that matches.
  • the upper first position is a position where the first representative position of the first supported portion coincides with the position in the horizontal direction in the plane continuous with the upper flange surface.
  • the upper second position is a position where the position in the horizontal direction coincides with the second representative position of the second supported portion in the plane continuous with the upper flange surface.
  • the upper symmetrical position is a position in the upper flange surface that coincides with the lower symmetrical position in the horizontal direction.
  • the upper symmetrical midpoint position is a midpoint in the lateral direction in the upper flange surface, and the position in the axial direction coincides with the lower symmetrical position.
  • the vertical position indicated by the effective three-dimensional coordinate data of the lower target position changed in the coordinate changing step and the vertical position indicated by the effective three-dimensional coordinate data of the target contact position.
  • the difference is defined as a vertical displacement amount of the lower target position
  • the vertical position indicated by the effective three-dimensional coordinate data of the upper target position changed in the coordinate changing step and the effective three-dimensional coordinate data of the target contact position are
  • the difference from the indicated position in the vertical direction is defined as the amount of vertical displacement of the upper target position.
  • the target contact position is the midpoint in the vertical direction between the upper symmetrical midpoint position on the upper flange surface and the lower symmetrical midpoint position on the lower flange surface.
  • the displacement amount of the upper symmetrical position is defined as the difference in the vertical direction between the upper symmetrical position and the target contact position where the displacement in the vertical direction is to be obtained on the upper flange surface.
  • the difference in the vertical direction between the target contact position and the lower target position where the vertical displacement is to be obtained on the lower flange surface is used as the displacement of the lower target position.
  • the finite element models of the lower half casing and the upper half casing are used to simulate the deformation of the lower half casing and the upper half casing.
  • a displacement amount can be obtained. Therefore, in this aspect, it is possible to suppress the calculation load when obtaining the displacement amount.
  • the target contact position the middle point in the vertical direction between the upper symmetrical position on the upper flange surface and the lower symmetrical position on the lower flange surface.
  • Deformation of the flange surface includes not only vertical deformation accompanying changes in the axial direction, but also vertical deformation accompanying changes in the lateral direction.
  • the lower symmetrical position and the upper symmetrical position are positions of the inner edge in the flange surface, and the symmetrical contact position is obtained using the upper symmetrical position and the lower symmetrical position as described above.
  • the target is the vertical midpoint between the upper symmetrical midpoint position, which is the midpoint in the lateral direction in the upper flange surface, and the lower symmetrical midpoint position, which is the lateral midpoint in the lower flange surface. in contact position.
  • the vertical deformation accompanying the lateral change in the flange surface is not extremely reflected in the target contact position to be obtained, and the error in the vertical direction of the target contact position can be reduced. , the error in the amount of displacement between the upper target position and the lower target position can be reduced.
  • This rotary machine includes a rotor rotatable about a horizontally extending axis, a casing covering the outer circumference of the rotor, stationary parts arranged in the casing and attached to the casing, and the casing being mounted downward. and a pedestal supporting from the side.
  • the casing has an upper half casing, a lower lower half casing, and a plurality of bolts that fasten the upper half casing and the lower half casing.
  • the upper half casing has an upper flange formed with an upper flange surface facing downward.
  • the lower half casing includes a lower flange formed with a lower flange surface facing upward and facing the upper flange surface in the vertical direction; a first supported portion and a second supported portion separated from each other in an extending axial direction;
  • the upper flange and the lower flange are formed with bolt holes penetrating in the vertical direction and through which the plurality of bolts can be respectively inserted.
  • the program for estimating the flange displacement amount of the rotating machine described above is After the rotating machine is disassembled, the measured three-dimensional coordinate data at a plurality of positions on the upper flange surface and the lower A measured coordinate receiving step of receiving measured three-dimensional coordinate data at a plurality of positions on the flange surface; and a lower first position and a lower second position using the measured three-dimensional coordinate data at the plurality of positions on the lower flange surface.
  • the effective three-dimensional coordinate data at the lower target position and the lower target midpoint position and using the actually measured three-dimensional coordinate data at a plurality of positions in the upper flange surface, the upper first position and the upper an effective coordinate grasping step of grasping effective three-dimensional coordinate data at a second position, an upper target position, and an upper target midpoint position; and an effective three-dimensional coordinate of said lower first position grasped in said effective coordinate grasping step.
  • the data match the effective three-dimensional coordinate data of the first upper position, and the effective three-dimensional coordinate data of the second lower position and the effective three-dimensional coordinate data of the second upper position grasped in the effective coordinate grasping step.
  • the lower first position is a position where the first representative position of the first supported portion coincides with the position in the horizontal direction in the plane continuous with the upper flange surface.
  • the lower second position is a position that coincides with the second representative position of the second supported portion in the horizontal direction in the plane that continues to the lower flange surface.
  • the lower symmetrical position is a position on the lower flange surface where it is desired to obtain an amount of displacement in the vertical direction when changing from the open state to the fastened state.
  • the lower symmetrical midpoint position is the position of the midpoint in the lateral direction, which is the horizontal direction and the direction perpendicular to the axial direction, in the lower flange surface, and the position in the axial direction is the same as the lower symmetrical position. It is the position that matches.
  • the upper first position is a position where the first representative position of the first supported portion coincides with the position in the horizontal direction in the plane continuous with the upper flange surface.
  • the upper second position is a position where the position in the horizontal direction coincides with the second representative position of the second supported portion in the plane continuous with the upper flange surface.
  • the upper symmetrical position is a position in the upper flange surface that coincides with the lower symmetrical position in the horizontal direction.
  • the upper symmetrical midpoint position is a midpoint in the lateral direction in the upper flange surface, and the position in the axial direction coincides with the lower symmetrical position.
  • the vertical position indicated by the effective three-dimensional coordinate data of the lower target position changed in the coordinate changing step and the vertical position indicated by the effective three-dimensional coordinate data of the target contact position.
  • the difference is defined as a vertical displacement amount of the lower target position
  • the vertical position indicated by the effective three-dimensional coordinate data of the upper target position changed in the coordinate changing step and the effective three-dimensional coordinate data of the target contact position are
  • the difference from the indicated position in the vertical direction is defined as the amount of vertical displacement of the upper target position.
  • a flange displacement amount device for a rotating machine as one aspect for achieving the above object is applied to the following rotating machines.
  • This rotary machine includes a rotor rotatable about a horizontally extending axis, a casing covering the outer circumference of the rotor, stationary parts arranged in the casing and attached to the casing, and the casing being mounted downward. and a pedestal supporting from the side.
  • the casing has an upper half casing, a lower lower half casing, and a plurality of bolts that fasten the upper half casing and the lower half casing.
  • the upper half casing has an upper flange formed with an upper flange surface facing downward.
  • the lower half casing includes a lower flange formed with a lower flange surface facing upward and facing the upper flange surface in the vertical direction; a first supported portion and a second supported portion separated from each other in an extending axial direction;
  • the upper flange and the lower flange are formed with bolt holes penetrating in the vertical direction and through which the plurality of bolts can be respectively inserted.
  • the flange displacement amount estimating device for the rotating machine described above is After the rotating machine is disassembled, the measured three-dimensional coordinate data at a plurality of positions on the upper flange surface and the lower A measured coordinate reception unit that receives measured three-dimensional coordinate data at a plurality of positions on the flange surface; and a lower first position and a lower second position using the measured three-dimensional coordinate data at a plurality of positions on the lower flange surface.
  • the effective three-dimensional coordinate data at the lower target position and the lower target midpoint position and using the actually measured three-dimensional coordinate data at a plurality of positions in the upper flange surface, the upper first position and the upper an effective coordinate grasping unit for grasping effective three-dimensional coordinate data at the second position, the upper target position, and the upper target midpoint position; and the effective three-dimensional coordinates of the lower first position grasped by the effective coordinate grasping unit.
  • the data match the effective three-dimensional coordinate data of the first upper position, and the effective three-dimensional coordinate data of the second lower position and the effective three-dimensional coordinate data of the second upper position grasped by the effective coordinate grasping unit are matched.
  • a coordinate changing unit that changes the effective three-dimensional coordinate data grasped by the effective coordinate grasping unit so that the effective cubic coordinates at the lower target midpoint position and the upper target midpoint position changed by the coordinate changing unit a contact position estimator that obtains effective three-dimensional coordinate data of a target contact position, which is an intermediate position in the vertical direction between the lower target midpoint position and the upper target midpoint position, using the original coordinate data; a displacement amount calculation unit that calculates an amount of vertical displacement of the upper symmetrical position and the lower symmetrical position when the upper half casing and the lower half casing are fastened with the plurality of bolts.
  • the lower first position is a position where the first representative position of the first supported portion coincides with the position in the horizontal direction in the plane continuous with the upper flange surface.
  • the lower second position is a position that coincides with the second representative position of the second supported portion in the horizontal direction in the plane that continues to the lower flange surface.
  • the lower symmetrical position is a position on the lower flange surface where it is desired to obtain an amount of displacement in the vertical direction when changing from the open state to the fastened state.
  • the lower symmetrical midpoint position is the position of the midpoint in the lateral direction, which is the horizontal direction and the direction perpendicular to the axial direction, in the lower flange surface, and the position in the axial direction is the same as the lower symmetrical position.
  • the upper first position is a position where the first representative position of the first supported portion coincides with the position in the horizontal direction in the plane continuous with the upper flange surface.
  • the upper second position is a position where the position in the horizontal direction coincides with the second representative position of the second supported portion in the plane continuous with the upper flange surface.
  • the upper symmetrical position is a position in the upper flange surface that coincides with the lower symmetrical position in the horizontal direction.
  • the upper symmetrical midpoint position is a midpoint in the lateral direction in the upper flange surface, and the position in the axial direction coincides with the lower symmetrical position.
  • the displacement amount calculator calculates the difference between the vertical position indicated by the effective three-dimensional coordinate data of the lower target position changed by the coordinate changing unit and the vertical position indicated by the effective three-dimensional coordinate data of the target contact position. is the vertical displacement amount of the lower target position, and the vertical position indicated by the effective three-dimensional coordinate data of the upper target position changed by the coordinate changing unit and the vertical position indicated by the effective three-dimensional coordinate data of the target contact position Let the difference from the position of the direction be the amount of displacement in the vertical direction of the upper symmetrical position.
  • FIG. 1 is a schematic diagram showing a schematic configuration of a steam turbine as a rotating machine in one embodiment according to the present disclosure
  • FIG. 1 is a schematic diagram showing a schematic outline of a steam turbine as a rotating machine in one embodiment according to the present disclosure
  • FIG. FIG. 4 is a plan view of the main part of the upper half casing and the main part of the lower half casing in one embodiment according to the present disclosure
  • 1 is a cross-sectional view of an open casing in one embodiment of the present disclosure
  • FIG. FIG. 4 is a cross-sectional view of a fastened casing in one embodiment of the present disclosure
  • 1 is a functional block diagram of a flange displacement amount estimation device in one embodiment according to the present disclosure
  • FIG. 4 is a flow chart showing procedures of a method for estimating a flange displacement amount according to an embodiment of the present disclosure
  • FIG. 4 is an explanatory diagram showing a position for grasping effective three-dimensional coordinate data on a flange surface in one embodiment according to the present disclosure
  • FIG. 4 is an explanatory diagram showing the processing contents of a contact position estimation step and a displacement amount calculation step in an embodiment according to the present disclosure
  • FIG. 10 is an explanatory diagram showing the processing contents in a coordinate changing step S3 in one embodiment according to the present disclosure
  • FIG. 5 is an explanatory diagram showing positions of actually measured three-dimensional shape data required when executing a first grasping method in an embodiment according to the present disclosure
  • FIG. 4 is an explanatory diagram showing how to obtain effective three-dimensional coordinate data when executing a first grasping method in an embodiment according to the present disclosure
  • FIG. 10 is an explanatory diagram showing another method of obtaining effective three-dimensional coordinate data when executing the first grasping method in one embodiment according to the present disclosure
  • FIG. 5 is an explanatory diagram showing positions of actually measured three-dimensional coordinate data required when executing a second grasping method in an embodiment according to the present disclosure
  • FIG. 5 is an explanatory diagram showing how to obtain effective three-dimensional coordinate data when executing a second grasping method in an embodiment according to the present disclosure
  • FIG. 10 is an explanatory diagram showing positions of actually measured three-dimensional coordinate data required when executing a third grasping method in an embodiment according to the present disclosure; Relative positions between the flange surface indicated by the reference three-dimensional shape data and the points indicated by the actually measured three-dimensional coordinate data at a plurality of positions on the actual flange surface when executing the third grasping method in an embodiment according to the present disclosure It is an image diagram showing a relationship.
  • FIG. 4 is an explanatory diagram for explaining a plurality of polygon data in one embodiment according to the present disclosure; FIG. FIG.
  • FIG. 4 is an explanatory diagram for explaining extraction of a plurality of specific polygon data from a plurality of polygon data in one embodiment according to the present disclosure; Actual measurement at a plurality of positions after polygon data extraction among a plurality of positions of the flange surface indicated by the reference three-dimensional shape data and the actual flange surface when executing the third grasping method in an embodiment according to the present disclosure
  • FIG. 4 is an image diagram showing a relative positional relationship with a point indicated by three-dimensional coordinate data
  • FIG. 11 is an explanatory diagram showing how to obtain a reference position when executing a third comprehension method in an embodiment according to the present disclosure;
  • FIG. 4 is an image diagram showing a relative positional relationship with a point indicated by three-dimensional coordinate data
  • FIG. 11 is an explanatory diagram showing how to obtain a reference position when executing a third comprehension method in an embodiment according to the present disclosure
  • FIG. 11 is an explanatory diagram showing positions of actually measured three-dimensional coordinate data required when executing a fourth grasping method in an embodiment according to the present disclosure; Relative positions between the flange surface indicated by the reference three-dimensional shape data and the points indicated by the measured three-dimensional coordinate data at a plurality of positions of the actual flange surface when executing the fourth grasping method in an embodiment according to the present disclosure It is an image diagram showing a relationship. Actual measurement at a plurality of positions after polygon data extraction among a plurality of positions of the flange surface indicated by the reference three-dimensional shape data and the actual flange surface when executing the fourth grasping method in an embodiment according to the present disclosure
  • FIG. 4 is an image diagram showing a relative positional relationship with a point indicated by three-dimensional coordinate data
  • FIG. 11 is an explanatory diagram showing how to obtain a reference position when executing a fourth grasping method in an embodiment according to the present disclosure;
  • Embodiments of a flange displacement estimation method for a rotating machine according to the present disclosure, a program for executing this method, and an apparatus for executing this method will be described below.
  • "Embodiment of Rotating Machinery" A rotary machine according to the present embodiment will be described with reference to FIGS. 1 to 5.
  • FIG. 1 A rotary machine according to the present embodiment will be described with reference to FIGS. 1 to 5.
  • the rotary machine of this embodiment is a steam turbine 10.
  • the steam turbine 10 includes a rotor 15 that rotates around a horizontally extending axis Ar, a casing 30 that covers the outer peripheral side of the rotor 15, a first bearing device 12a and a second bearing device that rotatably support the rotor 15. 12b, a plurality of diaphragms 20, a first shaft sealing device 13a and a second shaft sealing device 13b for sealing a gap between the casing 30 and the rotor 15, and a base 11 for supporting the casing 30 from below. .
  • the direction in which the axis Ar extends is the axial direction Dy
  • the horizontal direction perpendicular to the axial direction Dy is the lateral direction Dx
  • the circumferential direction with respect to the axis Ar is simply the circumferential direction Dc
  • the radial direction with respect to the axis Ar is simply the radial direction Dr
  • the side closer to the axis Ar is called the radial inner Dri
  • the side farther from the axis Ar is called the radial outer Dro.
  • U used in the symbols in the figure means the upper half
  • L means the lower half.
  • the rotor 15 has a rotor shaft 16 extending in the axial direction Dy, and a plurality of rotor blade rows 17 attached to the rotor shaft 16 aligned in the axial direction Dy.
  • Each of the plurality of rotor blade rows 17 has a plurality of rotor blades aligned in the circumferential direction Dc with respect to the axis Ar.
  • Both ends of the rotor shaft 16 protrude from the casing 30 in the axial direction Dy.
  • One end of the rotor shaft 16 in the axial direction Dy is rotatably supported by a first bearing device 12 a attached to the frame 11 .
  • the other end of the rotor shaft 16 in the axial direction Dy is rotatably supported by a second bearing device 12 b attached to the base 11 .
  • the first shaft sealing device 13a is provided at one end of the casing 30 in the axial direction Dy.
  • the second shaft sealing device 13b is provided at the other end of the casing 30 in the axial direction Dy.
  • Both the first shaft sealing device 13 a and the second shaft sealing device 13 b are devices for sealing the gap between the rotor shaft 16 and the casing 30 .
  • a plurality of diaphragms 20 are arranged in the axial direction Dy within the casing 30 .
  • Each of the plurality of diaphragms 20 has a lower half diaphragm 20L forming a portion below the axis Ar and an upper half diaphragm 20U forming a portion above the axis Ar.
  • Each of the lower half diaphragm 20L and the upper half diaphragm 20U includes a plurality of stator vanes 22 arranged in the circumferential direction Dc, a diaphragm inner ring 23 connecting radially inner Dri portions of the plurality of stator vanes 22, and a plurality of It has a diaphragm outer ring 24 that connects radially outer Dro portions of the stationary blades 22 to each other, and a sealing device 25 that is attached to the radially inner Dri of the diaphragm inner ring 23 .
  • This sealing device 25 is a sealing device that seals the gap between the diaphragm inner ring 23 and the rotor shaft 16 .
  • the first shaft sealing device 13a, the second shaft sealing device 13b, and the plurality of diaphragms 20 described above are all stationary parts that extend in the circumferential direction with respect to the axis Ar and are attached to the casing 30.
  • the casing 30 includes a lower half casing 30L constituting a portion below the axis Ar, an upper half casing 30U constituting a portion above the axis Ar, and a lower half casing 30L. and a plurality of bolts 39 for fastening the upper half casing 30U.
  • the lower half casing 30L includes a lower half casing main body 31L extending in the circumferential direction Dc, a lower flange 32L projecting radially outward Dro from both ends of the lower half casing main body 31L in the circumferential direction Dc, and a base 11 connected to the lower flange 32L. It has a first supported portion 35a and a second supported portion 35b that are supported from below by the .
  • the upper half casing 30U has an upper half casing body 31U extending in the circumferential direction Dc, and upper flanges 32U projecting radially outward Dro from both ends of the upper half casing body 31U in the circumferential direction Dc.
  • the upper flange 32U does not have a portion facing the first supported portion 35a and the second supported portion 35b in the lower flange 32L.
  • the upper flange 32U may be provided with portions of the lower flange 32L facing the first supported portion 35a and the second supported portion 35b.
  • the surface of the lower flange 32L facing upward forms a lower flange surface 33L.
  • the surface of the upper flange 32U facing downward forms an upper flange surface 33U.
  • the lower flange surface 33L and the upper flange surface 33U face each other in the vertical direction Dz.
  • the first supported portion 35a protrudes from one of both sides in the axial direction Dy of the lower flange 32L to one side.
  • the second supported portion 35b protrudes from the other side of both sides in the axial direction Dy of the lower flange 32L to the other side. Therefore, the second supported portion 35b is separated from the first supported portion 35a in the axial direction Dy.
  • the upper surface 35ap of the first supported portion 35a and the upper surface 35bp of the second supported portion 35b are surfaces that continue to the lower flange surface 33L. That is, the upper surface 35ap of the first supported portion 35a and the upper surface 35bp of the second supported portion 35b are continuous with the lower flange surface 33L and have no step with respect to the lower flange surface 33L.
  • the lower flange 32L and the upper flange 32U are formed with bolt holes 34 penetrating in the vertical direction Dz and through which the plurality of bolts 39 can be inserted.
  • the lower half casing 30L and the upper half casing 30U are fastened by bolts 39 inserted through the bolt holes 34 of the lower flange 32L and the upper flange 32U.
  • a plurality of stationary component storage portions 36 are formed on the inner peripheral surface of the lower half casing main body 31L and the inner peripheral surface of the upper half casing 30U to respectively store the above-described plurality of stationary components.
  • Each stationary component storage portion 36 of the lower casing main body 31L is a groove that is recessed radially outward Dro from the inner peripheral surface of the lower casing main body 31L and extends in the circumferential direction Dc.
  • each stationary component storage portion 36 of the upper half casing body 31U is a groove extending in the circumferential direction Dc, recessed radially outward Dro from the inner peripheral surface of the upper half casing body 31U.
  • the diaphragm 20, which is a type of stationary component, is supported by a portion near the flange surface of the stationary component storage portion 36 extending in the circumferential direction Dc.
  • the distance in the radial direction Dr between the stationary parts attached to the casing 30 and the rotor 15 must be within a predetermined allowable dimension.
  • the distance between the rotor shaft 16 and the first shaft sealing device 13a and the second shaft sealing device 13b, which are types of stationary parts, and the distance between the seal device 25 of the diaphragm 20 and the rotor shaft 16 The spacing must be within a predetermined allowable dimension.
  • shape data of the lower half casing 30L in the open state and shape data of the upper half casing 30U the shape of the lower half casing 30L and the upper half casing 30U changes when the casing 30 changes from the open state to the fastened state. If changed, the spacing in the radial direction Dr between the stationary part and the rotor 15 will change, and this spacing may fall outside the range of allowable dimensions.
  • the inventors have found that the change in the radial distance Dr between the stationary part and the rotor 15 due to the deformation of the lower casing 30L and the upper casing 30U due to the change from the open state to the fastened state is It was found to be dominant to the deformation of the flange surface 33U. Therefore, the inventor estimated the amount of displacement of the lower flange surface 33L and the amount of displacement of the upper flange surface 33U due to the change from the open state to the fastened state, and based on these displacement amounts, the static component and the rotor in the fastened state. 15 in the radial direction Dr.
  • a flange displacement amount estimating device and a flange displacement amount estimating method for estimating the displacement amount of the lower flange surface 33L and the displacement amount of the upper flange surface 33U will be described below.
  • the flange displacement estimation device 50 is a computer.
  • This flange displacement amount estimating device 50 includes a CPU (Central Processing Unit) 60 that performs various calculations, a memory 57 that serves as a work area for the CPU 60, an auxiliary storage device 58 such as a hard disk drive, and a keyboard and mouse.
  • a CPU Central Processing Unit
  • memory 57 that serves as a work area for the CPU 60
  • auxiliary storage device 58 such as a hard disk drive
  • keyboard and mouse a keyboard and mouse.
  • the input device (input device) 51 the display device (output device) 52, the input/output interface 53 of the manual input device 51 and the display device 52, and the three-dimensional shape measuring device 69 such as a three-dimensional laser measuring instrument
  • a device interface (input device) 54 for receiving and transmitting data
  • a communication interface (input/output device) 55 for communicating with the outside via the network N
  • a disk type that is a type of non-temporary storage medium.
  • a storage/reproduction device (input/output device) 56 for storing and reproducing data on the storage medium D is provided.
  • the auxiliary storage device 58 pre-stores a flange displacement amount estimation program 58p and reference three-dimensional shape data 58d for each of a plurality of parts constituting the steam turbine 10 .
  • the reference three-dimensional shape data 58d may be three-dimensional design data or, for example, three-dimensional data obtained by actual measurement before shipping the steam turbine 10 from the factory. That is, the reference three-dimensional shape data 58d may be three-dimensional data obtained before the operation before the periodic inspection. From the reference three-dimensional shape data 58d, it is possible to obtain three-dimensional coordinate data at each position of a plurality of parts.
  • the flange displacement amount estimation program 58p is loaded into the auxiliary storage device 58 via the storage/reproduction device 56, for example, from a disk-type storage medium D, which is a type of non-temporary storage medium.
  • the flange displacement amount estimation program 58p may be loaded into the auxiliary storage device 58 from an external device via the communication interface 55.
  • the CPU 60 functionally includes a measured coordinate reception unit 61 , an effective coordinate grasping unit 62 , a coordinate changing unit 63 , a contact position estimation unit 64 and a displacement amount calculation unit 65 .
  • Each of these functional units 61 to 65 functions when the CPU 60 executes a flange displacement amount estimation program 58p stored in the auxiliary storage device 58. FIG. Operations of these functional units 61 to 65 will be described later.
  • the steam turbine 10 is disassembled and assembled each time it is inspected.
  • the upper half casing 30U is removed from the lower half casing 30L as shown in FIG.
  • the casing 30 is in an open state in which the bolts 39 are not fastened between the upper half casing 30U and the lower half casing 30L.
  • the rotor 15 , the plurality of diaphragms 20 , the first shaft sealing device 13 a and the second shaft sealing device 13 b are removed from the casing 30 and arranged outside the casing 30 .
  • the lower casing 30L may be removed from the pedestal 11 when the disassembly of the steam turbine 10 is completed.
  • the operator uses a three-dimensional shape measuring device 69 such as a three-dimensional laser measurement device to measure the three-dimensional shape at a plurality of positions on the upper flange surface 33U.
  • Original coordinate values and three-dimensional coordinate values at a plurality of positions on the lower flange surface 33L are measured.
  • the operator uses the three-dimensional coordinate values at a plurality of positions on the upper flange surface 33U and the three-dimensional coordinate values at a plurality of positions on the lower flange surface 33L as actually measured three-dimensional coordinate data, and measures the flange displacement from the three-dimensional shape measuring device 69. It is transferred to the quantity estimation device 50 .
  • the measured coordinate reception unit 61 of the flange displacement amount estimation device 50 receives the measured three-dimensional coordinate data at a plurality of positions on the upper flange surface 33U and the measured three-dimensional coordinate data at a plurality of positions on the lower flange surface 33L (measured coordinate reception step S1).
  • the three-dimensional coordinate data in this embodiment includes a coordinate value indicating the position in the axial direction Dy extending in the horizontal direction, a coordinate value indicating the position in the vertical direction Dz perpendicular to the axial direction Dy, and a coordinate value indicating the position in the vertical direction Dz perpendicular to the axial direction Dy. and coordinate values indicating the position in the lateral direction Dx extending to the .
  • the effective coordinate grasping unit 62 of the flange displacement amount estimation device 50 uses the plurality of actually measured three-dimensional coordinate data to obtain a plurality of a lower target position 71L, a lower first position 72La, a lower second position 72Lb, a plurality of lower target midpoint positions 75L, a plurality of upper target positions 71U, an upper first position 72Ua, and an upper second position 72Ub and the effective three-dimensional coordinate data at a plurality of upper target midpoint positions 75U are grasped (effective coordinate grasping step S2).
  • the effective three-dimensional coordinate data refers to the cubic coordinates of points on the surfaces of the lower flange surface 33L and the upper flange surface 33U, including imaginary surfaces, calculated based on a plurality of actually-measured three-dimensional coordinate data received. This is the original coordinate data. This data is necessary for estimating the amount of displacement of the lower flange surface 33L and the amount of displacement of the upper flange surface 33U due to the change from the open state to the fastened state. A method for grasping this effective three-dimensional coordinate data will be described later in detail.
  • the lower first position 72La is a position that coincides with the first representative position 74a of the first supported portion 35a in the horizontal direction in the plane that continues to the lower flange surface 33L.
  • the first representative position 74a is the position where the largest load is applied in the first supported portion 35a.
  • the lower second position 72Lb is a position that coincides in the horizontal direction with the second representative position 74b of the second supported portion 35b in the plane that continues to the lower flange surface 33L.
  • the second representative position 74b is the position where the largest load is applied in the second supported portion 35b.
  • the "surface connected to the lower flange surface 33L" may be an actually existing surface or an imaginary surface.
  • the upper surface 35ap of the first supported portion 35a and the upper surface 35bp of the second supported portion 35b are surfaces that continue to the lower flange surface 33L.
  • the plurality of lower target positions 71L are positions where it is desired to obtain a displacement amount in the vertical direction Dz when the casing 30 is changed from the open state to the fastened state in the lower flange surface 33L.
  • the position where the displacement amount in the vertical direction Dz is to be obtained on the lower flange surface 33L is the position where the stationary component storage portion 36 is formed in the axial direction Dy on the lower flange surface 33L. position of the inner edge in 33L. As shown in FIGS.
  • the lower symmetrical midpoint position 75L is the position of the midpoint in the lateral direction Dx on the lower flange surface 33L, and the lower symmetrical position 71L and the position in the axial direction Dy match. position.
  • the upper first position 72Ua is a position that coincides in the horizontal direction with the first representative position 74a of the first supported portion 35a in the plane that continues to the upper flange surface 33U.
  • the upper second position 72Ub is a position that coincides in the horizontal direction with the second representative position 74b of the second supported portion 35b in the plane that continues to the upper flange surface 33U.
  • the "surface connected to the upper flange surface 33U" may be an actually existing surface or a virtual surface.
  • the plurality of upper target positions 71U are positions where it is desired to obtain a displacement amount in the vertical direction Dz when the casing 30 is changed from the open state to the fastened state in the upper flange surface 33U.
  • the position where the displacement amount in the vertical direction Dz is to be obtained on the upper flange surface 33U is the position where the stationary component storage portion 36 is formed in the axial direction Dy on the upper flange surface 33U. position of the inner edge in 33U.
  • the upper symmetrical midpoint position 75U is the position of the midpoint in the lateral direction Dx in the upper flange surface 33U, and the upper symmetrical position and the position in the axial direction Dy match. position.
  • a method for grasping effective three-dimensional coordinate data at each position will be described later in detail.
  • each of the plurality of upper target positions 71U coincides in horizontal position with any one of the plurality of lower target positions 71L. Therefore, the upper symmetrical midpoint position 75U is the position of the midpoint in the lateral direction Dx in the upper flange surface 33U, and is also the position where the lower symmetrical position 71L and the position in the axial direction Dy match.
  • matching the positions in the XX direction includes not only the case in which the positions in the XX direction are completely matching but also the case in which the positions in the XX direction are substantially matching.
  • matching the horizontal position means not only that the coordinate values indicating the position in the axial direction Dy are the same, but also that the coordinate values indicating the position in the lateral direction Dx are the same. and that the coordinate values indicating the positions in the horizontal direction Dx are substantially the same.
  • a change in the radial distance Dr between the stationary part and the rotor 15 due to the deformation of the lower half casing 30L and the upper half casing 30U due to the change from the open state to the fastened state is Deformation of the position of the inner edge of the lower flange surface 33L, which is the position where the stationary component storage portion 36 is formed, and the position of the upper flange surface 33U, where the stationary component storage portion 36 is formed in the axial direction Dy and is dominant to the deformation of the position of the inner edge in the upper flange surface 33U.
  • the lower target position 71L for which the displacement amount in the vertical direction Dz is desired is set to the above-described position
  • the upper target position 71U for which the displacement amount for the vertical direction Dz is desired is set to the above-described position.
  • the lower target position 71L does not have to be the position of the inner edge of the lower flange surface 33L. may be at any position within
  • the upper target position 71U does not have to be the position of the inner edge of the upper flange surface 33U. It can be anywhere within the range.
  • the flange width direction Dw is the direction along the flange surface that connects the outer edge and the inner edge of the flange surface, and is the direction in which the distance from the reference position to the outer edge or the inner edge of the flange surface is the shortest. Note that the reference positions are each of the upper target position 71U and the lower target position 71L.
  • the coordinate changing section 63 of the flange displacement amount estimating device 50 changes the effective three-dimensional coordinate data grasped by the effective coordinate grasping section 62 (coordinate changing step S3).
  • the coordinate changing unit 63 causes the effective three-dimensional coordinate data of the lower first position 72La and the effective three-dimensional coordinate data of the upper first position 72Ua to match, and the lower second position 72Lb Effective three-dimensional coordinate data grasped by the effective coordinate grasping unit 62 by coordinate transformation such as parallel movement and/or rotational movement so that the effective three-dimensional coordinate data of the upper second position 72Ub and the effective three-dimensional coordinate data of the upper second position 72Ub match change the data.
  • the contact position estimating unit 64 of the flange displacement amount estimating device 50 uses the effective three-dimensional coordinate data at the lower object midpoint position 75L and the upper object midpoint position 75U changed by the coordinate changing unit 63 to Effective three-dimensional coordinate data of the target contact position 73, which is the intermediate position in the vertical direction between the midpoint position 75L and the upper target midpoint position 75U, is obtained (contact position estimation step S4).
  • the displacement amount calculator 65 of the flange displacement amount estimating device 50 obtains the amount of displacement in the vertical direction Dz of the upper target position 71U and the lower target position 71L when the casing 30 changes from the open state to the fastened state. These displacement amounts are output in response to a request from (displacement amount calculation step S5). Specifically, as shown in FIG. 9 and the following equations, the displacement amount calculation unit 65 calculates the coordinate value ZL in the vertical direction Dz included in the effective three-dimensional coordinate data of the lower target position 71L after the coordinate change and the coordinate value ZL in the vertical direction Dz.
  • the difference from the coordinate value ZC in the vertical direction Dz included in the effective three-dimensional coordinate data of the target contact position 73 with respect to the target position 71L is defined as the displacement amount ZdL in the vertical direction Dz of the lower target position 71L.
  • the displacement amount calculation unit 65 calculates the coordinate value ZU in the vertical direction Dz included in the effective three-dimensional coordinate data of the upper target position 71U after the coordinate change, and the effective three-dimensional coordinates of the target contact position 73 with respect to the upper target position 71U.
  • the difference from the coordinate value ZC in the vertical direction Dz included in the data is set as the displacement amount ZdU in the vertical direction Dz of the upper target position 71U.
  • ZdL ZL-ZC
  • ZdU ZU-ZC
  • the first grasping method When the first grasping method is executed in the effective coordinate grasping step S2, in the measured coordinate receiving step S1, as shown in FIG. It accepts actually measured three-dimensional coordinate data at a position 71U and a plurality of upper midpoint positions 75Ux.
  • the lower midpoint position 75Lx is the position of the midpoint in the lateral direction Dx in the lower flange surface 33L.
  • the plurality of lower midpoint positions 75Lx are different in position in the axial direction Dy.
  • One of the plurality of lower midpoint positions 75Lx is the lower target midpoint position 75L.
  • the upper midpoint position 75Ux is the position of the midpoint in the lateral direction Dx in the upper flange surface 33U.
  • the plurality of upper midpoint positions 75Ux are different in position in the axial direction Dy.
  • One of the plurality of upper midpoint positions 75Ux is the upper target midpoint position 75U.
  • the effective coordinate grasping unit 62 directly obtains the measured three-dimensional coordinate data at the plurality of lower target positions 71L and the plurality of upper target positions 71U received in the measured coordinate receiving step S1. , a plurality of lower target positions 71L, and a plurality of upper target positions 71U.
  • the valid coordinate grasping unit 62 obtains the three-dimensional coordinates of the plurality of lower target positions 71L and the plurality of upper target positions 71U at the time when this data was created. data can be obtained. Therefore, for example, the effective coordinate grasping unit 62 selects the measured three-dimensional coordinate data of the lower target position 71L from among the measured three-dimensional coordinate data of the plurality of positions received by the measured coordinate receiving unit 61 in the following manner.
  • the effective coordinate grasping unit 62 selects the three-dimensional coordinate data of the one lower target position 71L indicated by the reference three-dimensional shape data 58d from among the measured three-dimensional coordinate data at the plurality of positions received by the measured coordinate receiving unit 61, and the horizontal direction coordinate data. is extracted, and this measured three-dimensional coordinate data is recognized as the measured three-dimensional coordinate data of the lower target position 71L.
  • the effective coordinate grasping unit 62 obtains the effective three-dimensional coordinate data at the lower target middle point position 75L from the change tendency of the actually measured three-dimensional coordinate data at the plurality of lower middle point positions 75Lx.
  • the effective coordinate grasping unit 62 further obtains the effective three-dimensional coordinate data at the upper target middle point position 75U from the change tendency of the effective three-dimensional coordinate data at the plurality of upper middle point positions 75Ux.
  • a high-order function F such as a quadratic function that approximately indicates the surface shape in a region along a plurality of lower midpoint positions 75Lx is obtained.
  • the effective coordinate grasping unit 62 obtains the coordinate value in the vertical direction Dz with respect to the coordinate value in the horizontal direction of the lower target midpoint position 75L indicated by the reference three-dimensional shape data 58d. Then, the effective coordinate grasping unit 62 replaces the coordinate value in the vertical direction Dz among the coordinate values in each direction with respect to the lower target midpoint position 75L indicated by the reference three-dimensional shape data 58d by using the previously obtained coordinate value in the vertical direction Dz. and use this as effective three-dimensional coordinate data for the lower target midpoint position 75L.
  • the measured coordinate receiving unit 61 can receive the measured three-dimensional coordinate data for all of the plurality of lower target midpoint positions 75L and the measured three-dimensional coordinate data for all of the plurality of upper target midpoint positions 75U, these measured The three-dimensional coordinate data may be used as the effective three-dimensional coordinate data.
  • one of the plurality of lower symmetrical middle point positions 75L may be the position of the bolt hole 34 or the like. In this case, the actually measured three-dimensional coordinate data of the lower target midpoint position 75L cannot be obtained. Therefore, here, the effective three-dimensional coordinate data at the lower target midpoint position 75L is obtained from the change tendency of the actually measured three-dimensional coordinate data at the plurality of lower midpoint positions 75Lx.
  • the effective coordinate ascertaining unit 62 uses the above-described high-order function F that indicates the change tendency of the actually measured three-dimensional coordinate data at the plurality of lower middle point positions 75Lx to determine the lower first position 72La and Valid three-dimensional coordinate data at the lower second position 72Lb is obtained by extrapolation.
  • the effective coordinate grasping unit 62 further uses a high-order function that indicates the tendency of change in the effective three-dimensional coordinate data at the plurality of upper midpoint positions 75Ux to determine the effective three-dimensional coordinates at the upper first position 72Ua and the upper second position 72Ub. Data are extrapolated.
  • the surface shapes of the lower flange surface 33L and the upper flange surface 33U are approximated by the high-order function F.
  • the surface shape of part of the lower flange surface 33L and the surface shape of part of the upper flange surface 33U may be approximated by a linear function.
  • the effective coordinate grasping unit 62 calculates the effective three-dimensional coordinates at a plurality of lower middle point positions 75Lx near the lower target middle point position 75L for which the effective three-dimensional coordinate data should be grasped, among the plurality of lower middle point positions 75Lx.
  • the lower target midpoint position 75L for which effective three-dimensional coordinate data should be grasped is approximated by linear functions Fa and Fb.
  • the effective coordinate grasping unit 62 obtains effective three-dimensional coordinate data at a plurality of upper middle point positions 75Ux near the upper target middle point position 75U for which effective three-dimensional coordinate data should be grasped, among the plurality of upper middle point positions 75Ux. is used to approximate the upper target midpoint position 75U for which effective three-dimensional coordinate data should be grasped by a linear function. Then, using this linear function, the coordinate value in the vertical direction Dz of the effective three-dimensional coordinate data at the upper target midpoint position 75U is obtained.
  • the effective coordinate grasping unit 62 obtains the effective three-dimensional coordinate data at the lower first position 72La and the lower second position 72Lb by extrapolation using the linear functions Fa and Fb described above, and uses the linear functions to obtain Valid three-dimensional coordinate data at the upper first position 72Ua and the upper second position 72Ub are obtained by extrapolation.
  • a plurality of lower target positions 71L, a plurality of lower target midpoint positions 75L, a lower first position 72La, a lower second position 72Lb, a plurality of upper target positions 71U, and a plurality of upper target midpoint positions Effective three-dimensional coordinate data at 75U, upper first position 72Ua, and upper second position 72Ub are grasped.
  • the number of three-dimensional coordinate data to be handled can be reduced. can be reduced.
  • the upper half casing 30U may also have a first supported portion and a second supported portion that are connected to the upper flange 32U.
  • the measured coordinate receiving step S1 the measured three-dimensional coordinate data of the upper first position 72Ua on the lower surface of the first supported portion contiguous with the upper flange surface 33U, and the When the actually measured three-dimensional coordinate data of the upper second position 72Ub on the lower surface of the two supported parts is received, each of these actually measured three-dimensional data is directly used as the effective three-dimensional coordinate data of the upper first position 72Ua and the upper second position.
  • 72Ub effective three-dimensional coordinate data may be used.
  • the "YY position" in the "Z virtual latitude line extending to” is the "YY position” indicated by the reference three-dimensional shape data 58d. Further, the number of positions on the virtual line for which the measured three-dimensional coordinate data are received in the measured coordinate receiving step S1 is, for example, 2 or more and less than 10.
  • the effective coordinate grasping unit 62 uses the plurality of measured three-dimensional coordinate data received in the measured coordinate receiving step S1 to obtain a plurality of lower target positions 71L and a plurality of upper target positions. Obtain valid three-dimensional coordinate data in 71U.
  • the effective coordinate grasping unit 62 uses actually measured three-dimensional coordinate data at a plurality of positions 78 on a virtual line 76 passing through the reference position 71 and extending in the flange width direction Dw to obtain coordinates on the virtual line 76.
  • a function F2 that approximately indicates coordinate values in the vertical direction Dz at a plurality of positions 78 is obtained.
  • the effective coordinate grasping unit 62 uses this function F2 to extrapolate the coordinate values in the vertical direction Dz at the reference position 71 .
  • the effective coordinate grasping unit 62 replaces the coordinate values in the vertical direction Dz among the coordinate values in each direction with respect to the reference position 71 indicated by the reference three-dimensional shape data 58d with the coordinate values in the vertical direction Dz obtained previously. , is used as effective three-dimensional coordinate data of the reference position 71 .
  • the actually measured three-dimensional coordinate data at the following positions are received in the measured coordinate receiving step S1.
  • the effective coordinate grasping unit 62 grasps the effective three-dimensional coordinate data at the plurality of lower target midpoint positions 75L and the plurality of upper target midpoint positions 75U as the reference positions described above, in the same manner as described above. .
  • the effective coordinate grasping unit 62 uses the effective three-dimensional coordinate data at the plurality of lower target midpoint positions 75L to calculate the effective three-dimensional coordinates at the lower first position 72La and the lower second position 72Lb. Estimate data. Furthermore, the effective coordinate grasping unit 62 estimates effective three-dimensional coordinate data at the upper first position 72Ua and the upper second position 72Ub using the effective three-dimensional coordinate data at the plurality of upper target midpoint positions 75U. When the effective coordinate grasping unit 62 estimates the effective three-dimensional coordinate data at the lower first position 72La and the lower second position 72Lb, the effective coordinate grasping unit 62 detects a plurality of lower objects in the same manner as the first grasping method.
  • Effective three-dimensional coordinate data at the lower first position 72La and the lower second position 72Lb are obtained from the trend of change in the actually measured three-dimensional coordinate data at the midpoint position 75L. Further, when the effective coordinate grasping unit 62 estimates the effective three-dimensional coordinate data at the upper first position 72Ua and the upper second position 72Ub, the effective coordinate grasping unit 62 uses a plurality of Effective three-dimensional coordinate data at the upper first position 72Ua and the upper second position 72Ub are obtained from the trend of change in the actually measured three-dimensional coordinate data at the upper target midpoint position 75U.
  • a plurality of lower target positions 71L, a plurality of lower target midpoint positions 75L, a lower first position 72La, a lower second position 72Lb, a plurality of upper target positions 71U, and a plurality of upper target midpoint positions Effective three-dimensional coordinate data at 75U, upper first position 72Ua, and upper second position 72Ub are grasped.
  • the actually measured three-dimensional coordinate data of the reference position is used as the effective three-dimensional coordinate data of this reference position.
  • the effective three-dimensional coordinate data of the reference position is susceptible to local shape changes and may contain large measurement errors.
  • the three-dimensional shape measuring device 69 is a three-dimensional laser measuring instrument
  • the three-dimensional shape measured by the three-dimensional laser measuring instrument The original position data will contain an error.
  • the three-dimensional coordinate data of the reference position 71 is estimated from the measured three-dimensional coordinate data at a plurality of positions, and this three-dimensional coordinate data is used as effective three-dimensional coordinate data. Therefore, the second grasping method is less likely to be affected by local shape changes than the first grasping method, and can suppress the possibility of including a large measurement error.
  • the upper half casing 30U may also have the first supported portion and the second supported portion that are continuous with the upper flange 32U.
  • this actually measured three-dimensional coordinate Using the data, effective three-dimensional coordinate data of the upper first position 72Ua and the upper second position 72Ub may be obtained.
  • a. Measured three-dimensional coordinate data at a plurality of positions on a virtual line passing through the upper first position 72Ua and extending in the projecting direction of the first supported portion b.
  • the measured three-dimensional coordinate data at a plurality of positions 78 over the entire lower flange surface 33L are collected.
  • actually measured three-dimensional coordinate data at a plurality of positions over the entire upper flange surface 33U. 17 shows the relative positional relationship between the flange surface 80 and the reference position 81 indicated by the reference three-dimensional shape data 58d and points 85 indicated by the actually measured three-dimensional coordinate data at a plurality of positions over the actual flange surface. It is an image diagram showing.
  • Polygon data is data that defines a plane of a polygon.
  • the effective coordinate grasping unit 62 connects a plurality of points 85 that are close to each other with line segments, and draws a polygonal plane surrounded by these line segments as a polygon. 86.
  • the effective coordinate grasping unit 62 next extracts a plurality of polygon data satisfying certain conditions from among the plurality of polygon data, as shown in FIG. In FIG. 19, the polygon 86a specified by the polygon data to be extracted is patterned, and the polygon 86b specified by the polygon data not to be extracted is not patterned. Also, the XY plane in FIG. 19 is a plane parallel to the flange surface 80 indicated by the reference three-dimensional shape data 58d.
  • the aforementioned condition is that the inclination of the polygon 86 specified by the polygon data with respect to the flange surface 80 indicated by the reference three-dimensional shape data 58d is within a predetermined inclination.
  • the effective coordinate grasping unit 62 first obtains the normal line n of the polygon 86 for each of the plurality of polygons 86 . Next, for each of the plurality of polygons 86, the effective coordinate grasping unit 62 obtains the angle ⁇ between the normal line n to the polygon 86 and the normal p to the flange surface 80 indicated by the reference three-dimensional shape data 58d. Then, the effective coordinate grasping unit 62 selects a plurality of polygon data in which the angle ⁇ between the normal p to the flange surface 80 and the normal n to the polygon 86 is within a predetermined angle (predetermined inclination) from among the plurality of polygon data. Extract.
  • This data extraction process is performed by extracting the measured three-dimensional coordinate data of the plurality of points 85 received in the measured coordinate receiving step S1 from the points in the wall of the edge of the flange surface and the inner peripheral surface of the bolt hole 34 passing through the flange surface. Executed to remove the measured 3D coordinate data at points in the middle. Therefore, the number of points 85 after this extraction process is less than the number of points 85 before it, as shown in FIG. In particular, in the reference shape model indicated by the reference three-dimensional shape data 58d, the number of points 85 after the extraction process is significantly smaller than the number of points 85 before that with respect to the surface 82 that is inclined with respect to the flange surface 80. Become.
  • the effective coordinate grasping unit 62 next divides the virtual three-dimensional space including the flange surface 80 into a plurality of three-dimensional blocks 83, as shown in FIG. Then, the effective coordinate grasping unit 62 determines a representative point 87 in the target three-dimensional block 83 for each of the plurality of three-dimensional blocks 83 . Specifically, the effective coordinate grasping unit 62 selects a plurality of points 85 included in the target three-dimensional block 83 among the plurality of points 85 included in the polygon 86a specified by the plurality of polygon data extracted in the extraction process. A point that is the median value of the points 85 is set as a representative point 87 in the target three-dimensional block 83 .
  • the representative point 87 may be determined by robust estimation or biweight estimation based on the Lorenz distribution of the plurality of points 85 included in the polygon 86a specified by the plurality of polygon data extracted by the extraction process.
  • the effective coordinate grasping unit 62 connects the representative points 87 of each of the plurality of three-dimensional blocks 83 with a plane or curved surface as a complementary surface, and obtains a surface of the complementary surface including the representative points 87 of each of the plurality of three-dimensional blocks 83.
  • Create shape data This surface shape data is represented by a function F3 that indicates the shape of the entire flange surface.
  • the effective coordinate grasping unit 62 obtains effective three-dimensional coordinate data at the aforementioned reference position 71 using the surface shape data of the entire flange surface represented by the function F3.
  • the reference positions 71 in the third grasping method are each of the plurality of lower target positions 71L, the plurality of lower target midpoint positions 75L, the plurality of upper target positions 71U, and the plurality of upper target midpoint positions 75U.
  • the effective coordinate grasping unit 62 uses the effective three-dimensional coordinate data at the plurality of lower target midpoint positions 75L to Effective three-dimensional coordinate data at the lower first position 72La and the lower second position 72Lb are estimated. Furthermore, the effective coordinate grasping unit 62 estimates effective three-dimensional coordinate data at the upper first position 72Ua and the upper second position 72Ub using the effective three-dimensional coordinate data at the plurality of upper target midpoint positions 75U.
  • a plurality of lower target positions 71L, a plurality of lower target midpoint positions 75L, a lower first position 72La, a lower second position 72Lb, a plurality of upper target positions 71U, and a plurality of upper target midpoint positions Effective three-dimensional coordinate data at 75U, upper first position 72Ua, and upper second position 72Ub are grasped.
  • the third grasping method is less susceptible to local shape changes than the second grasping method, and can suppress the possibility of including large measurement errors. Furthermore, in the third grasping method, even when there is a wide range of data loss due to obstacles or the like, effective three-dimensional coordinate data at the reference position can be grasped.
  • the overall surface shape data of the upper surface 35ap of the first supported portion 35a, the upper surface 35bp of the second supported portion 35b, and the lower flange surface 33L are obtained.
  • the effective three-dimensional coordinate data of the lower first position 72La and the lower second position 72Lb are obtained using the surface shape data of the entire surface represented by the function.
  • the upper half casing 30U may also have a first supported portion and a second supported portion that are connected to the upper flange 32U.
  • a plurality of positions over the entire lower surface of the first supported portion connected to the upper flange 32U and a plurality of positions over the entire lower surface of the second supported portion connected to the upper flange 32U the effective three-dimensional coordinate data of the upper first position 72Ua and the upper second position 72Ub may be obtained by the following method.
  • the reference measurement area 79 is also a lower measurement area including the lower target position 71L on the lower flange surface 33L and an upper measurement area including the upper target position 71U on the upper flange surface 33U. Furthermore, this reference measurement area 79 is also a lower midpoint measurement area including the lower target midpoint position 75L on the lower flange surface 33L, and an upper midpoint measurement area including the upper target midpoint position 75U on the upper flange surface 33U. But also.
  • the three-dimensional coordinate data of the reference position 71 here is the three-dimensional coordinate data of the reference position indicated by the standard three-dimensional shape data 58d.
  • the number of measured three-dimensional coordinate data in the reference measurement area 79 received in the measured coordinate receiving step S1 is, for example, 10 or more. Therefore, the number of actually measured three-dimensional coordinate data in the reference measurement area 79 received in the measured coordinate receiving step S1 of the fourth grasping method is the actual measured three-dimensional coordinates of the position on the virtual line received in the measured coordinate receiving step S1 of the second grasping method. More than the number of data.
  • the effective coordinate grasping unit 62 first uses actually measured three-dimensional coordinate data at a plurality of positions 78 to obtain a plurality of polygons in the same manner as the method described in the third grasping method. Data is created, and a plurality of polygon data satisfying certain conditions are extracted from the plurality of polygon data. As a result, the number of points 85 indicated by the actually measured three-dimensional coordinate data after this extraction process is smaller than the number of points 85 before it, as shown in FIG.
  • the effective coordinate grasping unit 62 next divides the virtual three-dimensional space including the flange surface 80 into a plurality of three-dimensional blocks 83, as shown in FIG. Then, the effective coordinate grasping unit 62 determines a representative point 87 in the target three-dimensional block 83 for each of the plurality of three-dimensional blocks 83 .
  • the effective coordinate grasping unit 62 connects the representative points 87 of each of the plurality of three-dimensional blocks 83 with a plane or curved surface as a complementary surface, and obtains a surface of the complementary surface including the representative points 87 of each of the plurality of three-dimensional blocks 83.
  • Create shape data This surface shape data is represented by a function F4 that indicates the shape within the reference measurement area 79 on the flange surface.
  • the effective coordinate grasping unit 62 obtains effective three-dimensional coordinate data at the aforementioned reference position 71 using the surface shape data represented by the function F4.
  • the effective coordinate grasping unit 62 obtains the effective three-dimensional coordinate data at a plurality of lower target midpoint positions 75L. is used to estimate effective three-dimensional coordinate data at the lower first position 72La and the lower second position 72Lb. Furthermore, the effective coordinate grasping unit 62 estimates effective three-dimensional coordinate data at the upper first position 72Ua and the upper second position 72Ub using the effective three-dimensional coordinate data at the plurality of upper target midpoint positions 75U.
  • a plurality of lower target positions 71L, a plurality of lower target midpoint positions 75L, a lower first position 72La, a lower second position 72Lb, a plurality of upper target positions 71U, and a plurality of upper target midpoint positions Effective three-dimensional coordinate data at 75U, upper first position 72Ua, and upper second position 72Ub are grasped.
  • the fourth grasping method is less susceptible to local shape changes than the second grasping method, and can suppress the possibility of including large measurement errors. Furthermore, in the fourth grasping method, even if there is a wide range of data loss due to an obstacle or the like, effective three-dimensional coordinate data at the reference position can be grasped.
  • a plurality of polygon data are created, and a certain condition is selected from among the plurality of polygon data. Extract multiple polygon data that satisfies Then, a representative point is determined from among the plurality of points indicated by the extracted polygon data, and the coordinate value in the vertical direction Dz at this representative point is set as the coordinate value in the vertical direction Dz at the lower first position 72La.
  • a representative point is determined, and the coordinate value of the representative point in the vertical direction Dz is calculated as the lower second position. 72Lb in the vertical direction Dz.
  • the upper half casing 30U may also have a first supported portion and a second supported portion that are connected to the upper flange 32U.
  • representative points on each surface are determined using actually measured three-dimensional coordinate data at a plurality of positions received, and the coordinate values of the representative points on each surface in the vertical direction Dz are calculated as above.
  • the effective three-dimensional coordinate data at the upper target midpoint position 75U and the lower target midpoint position 75L are grasped by the first grasping method, and the effective three-dimensional coordinate data at the upper target position 71U and the lower target position 71L are grasped by the second grasping method. Alternatively, it may be grasped by the fourth grasping method. Also, the effective three-dimensional coordinate data at the upper target midpoint position 75U and the lower target midpoint position 75L are grasped by the second grasping method, and the effective three-dimensional coordinate data at the upper target position 71U and the lower target position 71L are grasped by the first grasping method. Alternatively, it may be grasped by the fourth grasping method.
  • the effective three-dimensional coordinate data at the upper target midpoint position 75U and the lower target midpoint position 75L are grasped by the fourth grasping method, and the effective three-dimensional coordinate data at the upper target position 71U and the lower target position 71L are grasped by the first grasping method. Alternatively, it may be grasped by the second grasping method.
  • the target contact position 73 is the midpoint in the vertical direction Dz between the upper symmetrical midpoint position 75U on the upper flange surface 33U and the lower symmetrical midpoint position 75L on the lower flange surface 33L.
  • the difference in the vertical direction Dz between the target contact position 73 and the upper target position 71U where the displacement in the vertical direction Dz is desired on the upper flange surface 33U is defined as the displacement of the upper target position 71U.
  • the difference in the vertical direction Dz between the target contact position 73 and the lower target position 71L where the displacement in the vertical direction Dz is desired on the lower flange surface 33L is defined as the displacement amount of the lower target position 71L. Therefore, in the present embodiment, the finite element models of the lower half casing 30L and the upper half casing 30U are used to simulate the deformation of the lower half casing 30L and the upper half casing 30U. A displacement amount in the vertical direction Dz of the target position 71L can be obtained. Therefore, in this embodiment, it is possible to reduce the calculation load when obtaining the displacement amount. Therefore, in this embodiment, it is possible to shorten the preparation period for estimating the flange surface and suppress the estimation cost.
  • the target contact position 73 at the middle point in the vertical direction Dz between the upper target position 71U on the upper flange surface 33U and the lower target position 71L on the lower flange surface 33L.
  • the deformation of the flange surface includes not only deformation in the vertical direction Dz accompanying changes in the axial direction Dy, but also deformation in the vertical direction Dz accompanying changes in the lateral direction Dx, as shown in FIG.
  • the lower target position 71L and the upper target position 71U are the inner edge positions in the flange surface, and the target contact position 73 is obtained using the upper target position 71U and the lower target position 71L as described above.
  • the deformation in the vertical direction Dz accompanying the change in the lateral direction Dx on the flange surface is extremely reflected in the target contact position 73 to be obtained, and the error in the target contact position 73 in the vertical direction Dz increases.
  • the error in the amount of displacement between the target position 71U and the lower target position 71L may increase.
  • the target contact position 73 is the middle point in the vertical direction Dz.
  • deformation in the vertical direction Dz accompanying change in the lateral direction Dx on the flange surface is not extremely reflected in the target contact position 73 to be obtained, and an error in the vertical direction Dz of the target contact position 73 is can be reduced, and as a result, the error in the amount of displacement between the upper target position 71U and the lower target position 71L can be reduced.
  • This rotary machine includes a rotor 15 rotatable around a horizontally extending axis Ar, a casing 30 covering the outer circumference of the rotor 15, and stationary parts arranged in the casing 30 and attached to the casing 30. and a frame 11 that supports the casing 30 from below.
  • the casing 30 has an upper half casing 30U, a lower lower half casing 30L, and a plurality of bolts 39 that fasten the upper half casing 30U and the lower half casing 30L.
  • the upper half casing 30U has an upper flange 32U formed with an upper flange surface 33U facing downward.
  • the lower half casing 30L is connected to a lower flange 32L having a lower flange surface 33L facing upward and facing the upper flange surface 33U in the vertical direction Dz, and the lower flange 32L. and a first supported portion 35a and a second supported portion 35b separated from each other in the axial direction Dy in which the axis Ar extends.
  • the upper flange 32U and the lower flange 32L are formed with bolt holes 34 penetrating in the vertical direction Dz and through which the plurality of bolts 39 can be inserted.
  • the lower first position 72La is a position that coincides in the horizontal direction with the first representative position 74a of the first supported portion 35a in the plane continuous with the upper flange surface 33U.
  • the lower second position 72Lb is a position that coincides in the horizontal direction with the second representative position 74b of the second supported portion 35b in the plane that continues to the lower flange surface 33L.
  • the lower target position 71L is a position in the lower flange surface 33L where it is desired to obtain a displacement amount in the vertical direction Dz when changing from the open state to the fastened state.
  • the lower symmetrical midpoint position 75L is the position of the midpoint in the lateral direction Dx, which is the direction perpendicular to the axial direction Dy in the horizontal direction, in the lower flange surface 33L. This is the position where the positions in the direction Dy match.
  • the upper first position 72Ua is a position in the horizontal direction that coincides with the first representative position 74a of the first supported portion 35a in the plane connecting to the upper flange surface 33U.
  • the upper second position 72Ub is a position in the horizontal direction that coincides with the second representative position 74b of the second supported portion 35b in the plane connecting to the upper flange surface 33U.
  • the upper symmetrical position 71U is a position in the upper flange surface 33U that coincides with the lower symmetrical position 71L in the horizontal direction.
  • the upper symmetrical midpoint position 75U is a midpoint in the lateral direction Dx on the upper flange surface 33U, and is a position where the lower symmetrical position 71L and the position in the axial direction Dy match.
  • the position in the vertical direction Dz indicated by the effective three-dimensional coordinate data of the lower target position 71L changed in the coordinate changing step S3 and the vertical position indicated by the effective three-dimensional coordinate data of the target contact position 73 are changed.
  • the difference from the position in the direction Dz is defined as the amount of displacement in the vertical direction Dz of the lower target position 71L, and the position in the vertical direction Dz indicated by the effective three-dimensional coordinate data of the upper target position 71U changed in the coordinate changing step S3.
  • the difference from the position in the vertical direction Dz indicated by the effective three-dimensional coordinate data of the target contact position 73 is defined as the amount of displacement in the vertical direction Dz of the upper target position 71U.
  • the target contact position 73 is the midpoint in the vertical direction Dz between the upper symmetrical midpoint position 75U on the upper flange surface 33U and the lower symmetrical midpoint position 75L on the lower flange surface 33L.
  • the difference in the vertical direction Dz between the target contact position 73 and the upper target position 71U where the displacement in the vertical direction Dz is desired on the upper flange surface 33U is defined as the displacement of the upper target position 71U.
  • the difference in the vertical direction Dz between the target contact position 73 and the lower target position 71L where the displacement amount in the vertical direction Dz is to be obtained on the lower flange surface 33L is the displacement amount of the lower target position 71L.
  • the upper symmetrical position 71U and the lower symmetrical position 71U can be obtained without simulating the deformation of the lower half casing 30L and the upper half casing 30U using the finite element models of the lower half casing 30L and the upper half casing 30U.
  • a displacement amount in the vertical direction Dz of the position 71L can be obtained. Therefore, in this aspect, it is possible to suppress the calculation load when obtaining the displacement amount.
  • the target contact position 73 at the middle point in the vertical direction Dz between the upper target position 71U on the upper flange surface 33U and the lower target position 71L on the lower flange surface 33L.
  • the deformation of the flange surface includes not only deformation in the vertical direction Dz due to change in the axial direction Dy, but also deformation in the vertical direction Dz due to change in the lateral direction Dx.
  • the lower target position 71L and the upper target position 71U are the inner edge positions in the flange surface, and the target contact position 73 is obtained using the upper target position 71U and the lower target position 71L as described above.
  • the deformation in the vertical direction Dz accompanying the change in the lateral direction Dx on the flange surface is extremely reflected in the target contact position 73 to be obtained, and the error in the target contact position 73 in the vertical direction Dz increases.
  • the error in the amount of displacement between the target position 71U and the lower target position 71L may increase.
  • the target contact position 73 is the midpoint of the direction Dz.
  • the deformation in the vertical direction Dz accompanying the change in the lateral direction Dx on the flange surface is not extremely reflected in the target contact position 73 to be obtained, and the error in the vertical direction Dz of the target contact position 73 can be reduced.
  • the error in the amount of displacement between the upper target position 71U and the lower target position 71L can be reduced.
  • the lower target position 71L is a position where the stationary component is arranged in the axial direction Dy and a position of an inner edge in the lower flange surface 33L.
  • the method for estimating the amount of flange displacement of a rotating machine in the third aspect includes: In the method for estimating a flange displacement amount of a rotating machine according to the first aspect or the second aspect, in the measured coordinate receiving step S1, measured three-dimensional coordinate data at the lower target midpoint position 75L and the upper target midpoint position 75U are obtained. accept. In the effective coordinate grasping step S2, the measured three-dimensional coordinate data at the lower target midpoint position 75L is grasped as it is as effective three-dimensional coordinate data at the lower target midpoint position 75L, and the measured coordinate receiving step S1 acquires the The actually measured three-dimensional coordinate data at the upper target midpoint position 75U is grasped as it is as effective three-dimensional coordinate data at the upper target midpoint position 75U.
  • the measured three-dimensional coordinate data at the lower target midpoint position 75L and the upper target midpoint position 75U received in the measured coordinate receiving step S1 are used as they are for the effective three-dimensional coordinate data at the lower target midpoint position 75L and the upper target midpoint position 75U. Since it is grasped as coordinate data, the calculation load can be suppressed.
  • the method for estimating the flange displacement amount of a rotating machine in the fourth aspect includes: In the method for estimating a flange displacement amount of a rotating machine according to the first aspect or the second aspect, in the measured coordinate reception step S1, on a lower midpoint imaginary line passing through the lower target midpoint position 75L and extending in the flange width direction Dw Measured three-dimensional coordinate data at a plurality of positions are accepted, and actually measured three-dimensional coordinate data at a plurality of positions on an upper midpoint virtual line passing through the upper target midpoint position 75U and extending in the flange width direction Dw.
  • effective three-dimensional coordinate data at the lower object middle point position 75L is obtained from actually measured three-dimensional coordinate data at a plurality of positions on the lower middle point virtual line, and a plurality of coordinates on the upper middle point virtual line are obtained.
  • Effective three-dimensional coordinate data at the upper target midpoint position 75U is obtained from the actually measured three-dimensional coordinate data at the position.
  • the effective three-dimensional coordinate data of the lower target midpoint position 75L is obtained from the measured three-dimensional coordinate data at a plurality of positions on the lower midpoint virtual line, and the measured three-dimensional coordinate data at a plurality of positions on the upper midpoint virtual line. Then, the effective three-dimensional coordinate data of the upper target midpoint position 75U is obtained. Therefore, in this aspect, the effective three-dimensional coordinate data at the lower target midpoint position 75L and the upper target midpoint position 75U are less likely to be affected by local shape changes and may include large measurement errors. can be suppressed.
  • a method for estimating a flange displacement amount of a rotating machine in the fifth aspect In the method for estimating a flange displacement amount of a rotating machine according to the first aspect or the second aspect, in the measured coordinate reception step S1, in the lower middle point measurement region including the lower target middle point position 75L in the lower flange surface 33L and receive measured three-dimensional coordinate data at a plurality of positions in the upper middle point measurement area including the upper target middle point position 75U on the upper flange surface 33U.
  • the effective three-dimensional coordinates at the lower target middle point position 75L are calculated using actually measured three-dimensional coordinate data at a plurality of positions in the lower middle point measurement area received in the actually measured coordinate receiving step S1.
  • Data is obtained, and effective three-dimensional coordinate data at the upper target midpoint position 75U is obtained using actually measured three-dimensional coordinate data at a plurality of positions in the upper midpoint measurement area received in the measured coordinate receiving step S1.
  • the effective three-dimensional coordinate data of the lower target midpoint position 75L is obtained from the actually measured three-dimensional coordinate data at a plurality of positions in the lower midpoint measurement area, and the actually measured three-dimensional coordinate data at a plurality of positions in the upper midpoint measurement area.
  • Effective three-dimensional coordinate data of the upper target midpoint position 75U is obtained from the coordinate data. Therefore, in this aspect, the effective three-dimensional coordinate data at the lower target midpoint position 75L and the upper target midpoint position 75U are less likely to be affected by local shape changes and may include large measurement errors. can be suppressed.
  • the method for estimating the flange displacement amount of a rotating machine in the sixth aspect includes: In the method for estimating a flange displacement amount of a rotating machine according to the first aspect or the second aspect, in the measured coordinate receiving step S1, the position of the midpoint in the lateral direction Dx in the lower flange surface 33L, the axis line A plurality of lower midpoint positions 75Lx whose positions are different from each other in the direction Dy, and a plurality of upper midpoints whose positions are different from each other in the axial direction Dy and which are positions of midpoints in the lateral direction Dx on the upper flange surface 33U.
  • the actually measured three-dimensional coordinate data at positions 75Ux and 75Ux are accepted.
  • effective three-dimensional coordinate data at the lower target middle point position 75L is obtained from the change tendency of actually measured three-dimensional coordinate data at the plurality of lower middle point positions 75Lx, and the plurality of upper middle point positions are determined.
  • Effective three-dimensional coordinate data at the upper target midpoint position 75U is obtained from the change tendency of the effective three-dimensional coordinate data at 75Ux.
  • effective three-dimensional coordinate data at the lower target midpoint position 75L is obtained from the trend of change in actually measured three-dimensional coordinate data at the plurality of lower midpoint positions 75Lx, and effective three-dimensional coordinate data at the plurality of upper midpoint positions 75Ux.
  • Effective three-dimensional coordinate data at the upper target midpoint position 75U is obtained from the change tendency of . Therefore, in this embodiment, even if the actually measured three-dimensional coordinate data at the lower target midpoint position 75L and the upper target midpoint position 75U are not received in the measured coordinate receiving step S1, the effective three-dimensional coordinate data at these positions can be grasped. can do.
  • the method for estimating the amount of flange displacement of a rotating machine in the seventh aspect includes: In the method for estimating a flange displacement amount of a rotating machine according to any one of the first to sixth aspects, in the measured coordinate receiving step S1, the lower first position 72La and the lower second position 72Lb are actually measured. Accepts three-dimensional coordinate data. In the effective coordinate grasping step S2, the actually measured three-dimensional coordinate data at the lower first position 72La and the lower second position 72Lb acquired in the actually measured coordinate receiving step S1 are used as they are for the lower first position 72La and the lower second position. It is grasped as effective three-dimensional coordinate data in 72Lb.
  • the measured three-dimensional coordinate data at the lower first position 72La and the lower second position 72Lb received in the measured coordinate receiving step S1 are grasped as they are as effective three-dimensional coordinate data at the lower first position 72La and the lower second position 72Lb. Therefore, the computational load can be reduced.
  • the method for estimating the flange displacement amount of a rotating machine in the eighth aspect includes: In the method for estimating a flange displacement amount of a rotating machine according to any one of the first to sixth aspects, in the measured coordinate reception step S1, at a plurality of positions on the surface of the first supported portion 35a The actually measured three-dimensional coordinate data and the actually measured three-dimensional coordinate data at a plurality of positions on the surface of the second supported portion 35b are accepted.
  • the surface of the first supported portion 35a and the surface of the second supported portion 35b are the surfaces of the upper flange surface 33U and the lower flange surface 33L that are continuous with the flange surface of the one half casing. is.
  • the effective three-dimensional coordinates at the lower first position 72La are determined from the measured three-dimensional coordinate data at a plurality of positions on the surface of the first supported portion 35a acquired in the measured coordinate receiving step S1.
  • Data is obtained, and effective three-dimensional coordinate data of the lower second position 72Lb is obtained from the measured three-dimensional coordinate data at a plurality of positions on the surface of the second supported portion 35b obtained in the measured coordinate receiving step S1.
  • the effective three-dimensional coordinate data of the lower first position 72La is obtained from a plurality of measured three-dimensional coordinate data in the upper surface 35ap of the first supported portion 35a received in the measured coordinate receiving step S1, and the measured coordinate receiving step S1.
  • Effective three-dimensional coordinate data of the lower second position 72Lb is obtained from a plurality of actually measured three-dimensional coordinate data in the upper surface 35bp of the second supported portion 35b received in . Therefore, in this aspect, the effective three-dimensional coordinate data at the lower first position 72La and the lower second position 72Lb are less likely to be affected by local shape changes, and the possibility of including a large measurement error is suppressed. can be done.
  • the method for estimating the flange displacement amount of a rotating machine in the ninth aspect includes: In the method for estimating the flange displacement amount of a rotating machine according to the sixth aspect, in the effective coordinate grasping step S2, from the change tendency of the effective three-dimensional coordinate data at the plurality of lower middle point positions 75Lx, the lower first position 72La and the Effective three-dimensional coordinate data at the lower second position 72Lb is obtained, and from the change tendency of the effective three-dimensional coordinate data at the plurality of upper midpoint positions 75Ux, the effective three-dimensional coordinate data at the upper first position 72Ua and the upper second position 72Ub Find coordinate data.
  • the effective three-dimensional coordinates at the lower first position 72La and the lower second position 72Lb are obtained from the change tendency of the effective three-dimensional coordinate data at the plurality of lower middle point positions 75Lx, and the effective three-dimensional coordinates at the plurality of upper middle point positions 75Ux are obtained.
  • Effective three-dimensional coordinate data at the upper first position 72Ua and the upper second position 72Ub are obtained from the change tendency of the three-dimensional coordinate data. Therefore, in this aspect, even if the measured three-dimensional coordinate data at the lower first position 72La, the lower second position 72Lb, the upper first position 72Ua, and the upper second position 72Ub are not received in the measured coordinate receiving step S1, Valid three-dimensional coordinate data at these positions can be grasped.
  • a method for estimating a flange displacement amount of a rotary machine in the tenth aspect includes: In the method for estimating a flange displacement amount of a rotating machine according to any one of the first to ninth aspects, in the measured coordinate receiving step S1, the measured three-dimensional coordinates at the lower target position 71L and the upper target position 71U are measured. Accepts coordinate data. In the effective coordinate grasping step S2, the actually measured three-dimensional coordinate data at the lower target position 71L is grasped as it is as the effective three-dimensional coordinate data at the lower target position 71L, and the upper target position 71U acquired in the actually measured coordinate receiving step S1. The actually measured three-dimensional coordinate data in .
  • the actually measured three-dimensional coordinate data at the lower target position 71L and the upper target position 71U received in the measured coordinate receiving step S1 are grasped as they are as effective three-dimensional coordinate data at the lower target position 71L and the upper target position 71U. load can be reduced.
  • the method for estimating the flange displacement amount of a rotating machine in the eleventh aspect includes: In the method for estimating a flange displacement amount of a rotating machine according to any one of the first to ninth aspects, in the measured coordinate reception step S1, a lower position passing through the lower target position 71L and extending in the flange width direction Dw Measured three-dimensional coordinate data at a plurality of positions on the virtual line are accepted, and actually measured three-dimensional coordinate data at a plurality of positions on the upper virtual line passing through the upper target position 71U and extending in the flange width direction Dw are accepted.
  • effective three-dimensional coordinate data at the lower target position 71L is obtained from measured three-dimensional coordinate data at a plurality of positions on the lower virtual line, and measured three-dimensional coordinates at a plurality of positions on the upper virtual line.
  • Effective three-dimensional coordinate data at the upper target position 71U is obtained from the data.
  • effective three-dimensional coordinate data at the lower target position 71L is obtained from measured three-dimensional coordinate data at a plurality of positions on the lower virtual line
  • effective three-dimensional coordinate data at the upper target position 71U from the measured three-dimensional coordinate data at a plurality of positions on the upper virtual line
  • the effective three-dimensional coordinate data at the lower target position 71L and the upper target position 71U are less likely to be affected by local shape changes, and the possibility of including a large measurement error can be suppressed. .
  • the method for estimating the flange displacement amount of a rotating machine in the twelfth aspect includes: In the method for estimating a flange displacement amount of a rotating machine according to any one of the first to ninth aspects, in the measured coordinate reception step S1, the lower flange surface 33L includes the lower target position 71L. The measured three-dimensional coordinate data at a plurality of positions in the measurement area are accepted, and the measured three-dimensional coordinate data at a plurality of positions in the upper measurement area including the upper target position 71U on the upper flange surface 33U are accepted.
  • effective three-dimensional coordinate data at the lower target position 71L is obtained by using actually measured three-dimensional coordinate data at a plurality of positions in the lower measurement area received in the measured coordinate receiving step S1
  • Effective three-dimensional coordinate data at the upper target position 71U is obtained using actually measured three-dimensional coordinate data at a plurality of positions in the upper measurement area received in the actually measured coordinate receiving step S1.
  • the effective three-dimensional coordinate data of the lower target position 71L is obtained from the measured three-dimensional coordinate data at a plurality of positions in the lower measurement region
  • the upper target position is obtained from the measured three-dimensional coordinate data at a plurality of positions in the upper measurement region.
  • the method for estimating the flange displacement amount of a rotating machine in the thirteenth aspect includes: In the method for estimating a flange displacement amount of a rotating machine according to the first aspect or the second aspect, in the measured coordinate receiving step S1, the measured three-dimensional coordinate data at a plurality of positions over the entire lower flange surface 33L are received, and the It accepts actually measured three-dimensional coordinate data at a plurality of positions over the entire upper flange surface 33U. In the effective coordinate grasping step S2, the three-dimensional shape of the entire lower flange surface 33L is indicated by using the measured three-dimensional coordinate data at a plurality of positions over the entire lower flange surface 33L received in the measured coordinate receiving step S1.
  • the shape data of the lower flange surface 33L is obtained using the measured three-dimensional coordinate data at a plurality of positions over the entire upper flange surface 33U received in the measured coordinate receiving step S1.
  • Shape data of the upper flange surface 33U indicating the shape is obtained.
  • the shape data of the lower flange surface 33L is used to determine the effective three-dimensional coordinate data of the lower target midpoint position 75L
  • the shape data of the upper flange surface 33U is used to obtain the effective three-dimensional coordinate data.
  • the effective three-dimensional coordinate data of the upper target midpoint position 75U is obtained.
  • the effective three-dimensional coordinate data of the lower target midpoint position 75L is obtained from the measured three-dimensional coordinate data at a plurality of positions over the entire lower flange surface 33L, and the actual measurement at a plurality of positions over the entire upper flange surface 33U.
  • Effective three-dimensional coordinate data of the upper target midpoint position 75U is obtained from the three-dimensional coordinate data. Therefore, in this aspect, the effective three-dimensional coordinate data at the lower target midpoint position 75L and the upper target midpoint position 75U are less likely to be affected by local shape changes and may include large measurement errors. can be suppressed. Furthermore, in this aspect, even if there is a wide range of missing data due to obstacles or the like, effective three-dimensional coordinate data at these positions can be grasped.
  • the method for estimating the flange displacement amount of a rotating machine in the fourteenth aspect includes: In the method for estimating the flange displacement amount of a rotating machine according to the thirteenth aspect, in the effective coordinate grasping step S2, using the shape data of the lower flange surface 33L, Effective three-dimensional coordinate data is obtained, and effective three-dimensional coordinate data at the upper first position 72Ua and the upper second position 72Ub are obtained using the shape data of the upper flange surface 33U.
  • the effective three-dimensional coordinate data of the lower first position 72La and the lower second position 72Lb are obtained from the measured three-dimensional coordinate data at a plurality of positions over the entire lower flange surface 33L, and the effective three-dimensional coordinate data are obtained over the entire upper flange surface 33U.
  • Effective three-dimensional coordinate data of the upper first position 72Ua and the upper second position 72Ub are obtained from the actually measured three-dimensional coordinate data at a plurality of positions. Therefore, in this aspect, even if the measured three-dimensional coordinate data at the lower first position 72La, the lower second position 72Lb, the upper first position 72Ua, and the upper second position 72Ub are not received, effective three-dimensional coordinates at these positions Coordinate data can be grasped.
  • the effective three-dimensional coordinate data at the lower first position 72La, the lower second position 72Lb, the upper first position 72Ua, and the upper second position 72Ub are less likely to be affected by local shape changes. can reduce the possibility of including large measurement errors. Furthermore, in this aspect, even if there is a wide range of missing data due to obstacles or the like, effective three-dimensional coordinate data at these positions can be grasped.
  • the method for estimating the flange displacement amount of a rotating machine in the fifteenth aspect includes: In the method for estimating a flange displacement amount of a rotating machine according to the thirteenth aspect or the fourteenth aspect, in the effective coordinate grasping step S2, using the shape data of the lower flange surface 33L, the effective three-dimensional position at the lower target position 71L is calculated. Along with obtaining the coordinate data, using the shape data of the upper flange surface 33U, effective three-dimensional coordinate data at the upper target position 71U is obtained.
  • the effective three-dimensional coordinate data of the lower target position 71L is obtained from the actually measured three-dimensional coordinate data at a plurality of positions over the entire lower flange surface 33L, and the actually measured three-dimensional coordinate data at a plurality of positions over the entire upper flange surface 33U.
  • Effective three-dimensional coordinate data of the upper target position 71U is obtained from the coordinate data. Therefore, in this aspect, the effective three-dimensional coordinate data at the lower target position 71L and the upper target position 71U are less likely to be affected by local shape changes, and the possibility of including a large measurement error can be suppressed. . Furthermore, in this aspect, even if there is a wide range of missing data due to obstacles or the like, effective three-dimensional coordinate data at these positions can be grasped.
  • a program for estimating the amount of flange displacement of a rotary machine in the above embodiment is grasped, for example, as follows.
  • This rotary machine includes a rotor 15 rotatable around a horizontally extending axis Ar, a casing 30 covering the outer circumference of the rotor 15, and stationary parts arranged in the casing 30 and attached to the casing 30. and a frame 11 that supports the casing 30 from below.
  • the casing 30 has an upper half casing 30U, a lower lower half casing 30L, and a plurality of bolts 39 that fasten the upper half casing 30U and the lower half casing 30L.
  • the upper half casing 30U has an upper flange 32U formed with an upper flange surface 33U facing downward.
  • the lower half casing 30L is connected to a lower flange 32L having a lower flange surface 33L facing upward and facing the upper flange surface 33U in the vertical direction Dz, and the lower flange 32L. and a first supported portion 35a and a second supported portion 35b separated from each other in the axial direction Dy in which the axis Ar extends.
  • the upper flange 32U and the lower flange 32L are formed with bolt holes 34 penetrating in the vertical direction Dz and through which the plurality of bolts 39 can be inserted.
  • the flange displacement amount estimating program 58p of the rotating machine described above is configured to perform the above-described A measured coordinate reception step S1 for receiving measured three-dimensional coordinate data at a plurality of positions on the upper flange surface 33U and measured three-dimensional coordinate data at a plurality of positions on the lower flange surface 33L; Using the actually measured three-dimensional coordinate data, the effective three-dimensional coordinate data at the lower first position 72La, the lower second position 72Lb, the lower target position 71L, and the lower target midpoint position 75L are grasped, and the upper Using the actually measured three-dimensional coordinate data at a plurality of positions on the flange surface 33U, effective three-dimensional data at an upper first position 72Ua, an upper second position 72Ub, an upper target position 71U, and an upper target midpoint position 75U An effective coordinate grasping step S2 of grasping coordinate data, and the effective three-dimensional coordinate data of the lower first position 72La grasped in the effective coordinate grasping step S2 and the effective three-dimensional coordinate
  • the effective coordinate grasping step S2 grasped so that the effective three-dimensional coordinate data of the lower second position 72Lb grasped in the effective coordinate grasping step S2 and the effective three-dimensional coordinate data of the upper second position 72Ub match.
  • a coordinate changing step S3 for changing the effective three-dimensional coordinate data and the effective three-dimensional coordinate data at the lower target midpoint position 75L and the upper target midpoint position 75U changed in the coordinate changing step S3, the a contact position estimation step S4 for obtaining effective three-dimensional coordinate data of a target contact position 73, which is an intermediate position in the vertical direction Dz between the lower target midpoint position 75L and the upper target midpoint position 75U;
  • the lower first position 72La is a position that coincides in the horizontal direction with the first representative position 74a of the first supported portion 35a in the plane continuous with the upper flange surface 33U.
  • the lower second position 72Lb is a position that coincides in the horizontal direction with the second representative position 74b of the second supported portion 35b in the plane that continues to the lower flange surface 33L.
  • the lower target position 71L is a position in the lower flange surface 33L where it is desired to obtain a displacement amount in the vertical direction Dz when changing from the open state to the fastened state.
  • the lower symmetrical midpoint position 75L is the position of the midpoint in the lateral direction Dx, which is the direction perpendicular to the axial direction Dy in the horizontal direction, in the lower flange surface 33L. This is the position where the positions in the direction Dy match.
  • the upper first position 72Ua is a position in the horizontal direction that coincides with the first representative position 74a of the first supported portion 35a in the plane connecting to the upper flange surface 33U.
  • the upper second position 72Ub is a position in the horizontal direction that coincides with the second representative position 74b of the second supported portion 35b in the plane connecting to the upper flange surface 33U.
  • the upper symmetrical position 71U is a position in the upper flange surface 33U that coincides with the lower symmetrical position 71L in the horizontal direction.
  • the upper symmetrical midpoint position 75U is a midpoint in the lateral direction Dx on the upper flange surface 33U, and is a position where the lower symmetrical position 71L and the position in the axial direction Dy match.
  • the position in the vertical direction Dz indicated by the effective three-dimensional coordinate data of the lower target position 71L changed in the coordinate changing step S3 and the vertical position indicated by the effective three-dimensional coordinate data of the target contact position 73 are changed.
  • the difference from the position in the direction Dz is defined as the amount of displacement in the vertical direction Dz of the lower target position 71L, and the position in the vertical direction Dz indicated by the effective three-dimensional coordinate data of the upper target position 71U changed in the coordinate changing step S3.
  • the difference from the position in the vertical direction Dz indicated by the effective three-dimensional coordinate data of the target contact position 73 is defined as the amount of displacement in the vertical direction Dz of the upper target position 71U.
  • the rotating machine flange displacement estimation program in the seventeenth aspect is in the flange displacement amount estimation program 58p of the rotating machine according to the sixteenth aspect, the lower target position 71L is the position where the stationary component is arranged in the axial direction Dy and the position of the inner edge in the lower flange surface 33L. is.
  • the rotating machine flange displacement amount estimation program in the eighteenth aspect includes: In the rotating machine flange displacement amount estimation program 58p according to the sixteenth aspect or the seventeenth aspect, in the measured coordinate reception step S1, the position of the middle point in the lateral direction Dx in the lower flange surface 33L is , a plurality of lower midpoint positions 75Lx having different positions in the axial direction Dy, and a plurality of midpoint positions in the lateral direction Dx on the upper flange surface 33U having different positions in the axial direction Dy. Measured three-dimensional coordinate data at the upper midpoint position 75Ux is accepted.
  • effective three-dimensional coordinate data at the lower target middle point position 75L is obtained from the change tendency of actually measured three-dimensional coordinate data at the plurality of lower middle point positions 75Lx, and the plurality of upper middle point positions are obtained.
  • Effective three-dimensional coordinate data at the upper target midpoint position 75U is obtained from the change tendency of the effective three-dimensional coordinate data at 75Ux.
  • the computer by causing the computer to execute this program, in the measured coordinate receiving step S1, as in the sixth aspect, the actually measured three-dimensional coordinate data at the lower object midpoint position 75L and the upper object midpoint position 75U are accepted. Even if not, valid three-dimensional coordinate data at these positions can be grasped.
  • the rotating machine flange displacement amount estimation program in the nineteenth aspect includes: In the program 58p for estimating the flange displacement amount of the rotating machine according to the eighteenth aspect, in the effective coordinate grasping step S2, the lower first position 72La and effective three-dimensional coordinate data at the lower second position 72Lb, and from the change tendency of the effective three-dimensional coordinate data at the plurality of upper midpoint positions 75Ux, effective at the upper first position 72Ua and the upper second position 72Ub Obtain three-dimensional coordinate data.
  • the apparatus for estimating the amount of flange displacement for a rotating machine in the above embodiment is grasped, for example, as follows.
  • the flange displacement amount estimating device for a rotating machine is applied to the following rotating machines.
  • This rotary machine includes a rotor 15 rotatable around a horizontally extending axis Ar, a casing 30 covering the outer circumference of the rotor 15, and stationary parts arranged in the casing 30 and attached to the casing 30. and a frame 11 that supports the casing 30 from below.
  • the casing 30 has an upper half casing 30U, a lower lower half casing 30L, and a plurality of bolts 39 that fasten the upper half casing 30U and the lower half casing 30L.
  • the upper half casing 30U has an upper flange 32U formed with an upper flange surface 33U facing downward.
  • the lower half casing 30L is connected to a lower flange 32L having a lower flange surface 33L facing upward and facing the upper flange surface 33U in the vertical direction Dz, and the lower flange 32L. and a first supported portion 35a and a second supported portion 35b separated from each other in the axial direction Dy in which the axis Ar extends.
  • the upper flange 32U and the lower flange 32L are formed with bolt holes 34 penetrating in the vertical direction Dz and through which the plurality of bolts 39 can be inserted.
  • the above-mentioned A measured coordinate reception unit 61 that receives measured three-dimensional coordinate data at a plurality of positions on the upper flange surface 33U and measured three-dimensional coordinate data at a plurality of positions on the lower flange surface 33L; Using the actually measured three-dimensional coordinate data, the effective three-dimensional coordinate data at the lower first position 72La, the lower second position 72Lb, the lower target position 71L, and the lower target midpoint position 75L are grasped, and the upper Using the actually measured three-dimensional coordinate data at a plurality of positions on the flange surface 33U, effective three-dimensional data at an upper first position 72Ua, an upper second position 72Ub, an upper target position 71U, and an upper target midpoint position 75U An effective coordinate grasping section 62 grasping coordinate data, and the effective three-dimensional coordinate data of the lower first position 72La grasped by the effective coordinate grasping section
  • the effective coordinate grasping unit 62 grasps such that the effective three-dimensional coordinate data of the lower second position 72Lb grasped by the effective coordinate grasping unit 62 and the effective three-dimensional coordinate data of the upper second position 72Ub match.
  • a coordinate changing section 63 for changing effective three-dimensional coordinate data and the effective three-dimensional coordinate data at the lower target midpoint position 75L and the upper target midpoint position 75U changed by the coordinate changing section 63
  • the lower a contact position estimator 64 for obtaining effective three-dimensional coordinate data of a target contact position 73, which is an intermediate position in the vertical direction Dz between the target midpoint position 75L and the upper target midpoint position 75U
  • a displacement amount calculation unit 65 that calculates the amount of displacement in the vertical direction Dz of the upper symmetrical position 71U and the lower symmetrical position 71L when the 30U and the lower half casing 30L are fastened with the plurality of bolts 39.
  • the lower first position 72La is a position that coincides in the horizontal direction with the first representative position 74a of the first supported portion 35a in the plane continuous with the upper flange surface 33U.
  • the lower second position 72Lb is a position that coincides in the horizontal direction with the second representative position 74b of the second supported portion 35b in the plane that continues to the lower flange surface 33L.
  • the lower target position 71L is a position in the lower flange surface 33L where it is desired to obtain a displacement amount in the vertical direction Dz when changing from the open state to the fastened state.
  • the lower symmetrical midpoint position 75L is the position of the midpoint in the lateral direction Dx, which is the direction perpendicular to the axial direction Dy in the horizontal direction, in the lower flange surface 33L. This is the position where the positions in the direction Dy match.
  • the upper first position 72Ua is a position in the horizontal direction that coincides with the first representative position 74a of the first supported portion 35a in the plane connecting to the upper flange surface 33U.
  • the upper second position 72Ub is a position in the horizontal direction that coincides with the second representative position 74b of the second supported portion 35b in the plane connecting to the upper flange surface 33U.
  • the upper symmetrical position 71U is a position in the upper flange surface 33U that coincides with the lower symmetrical position 71L in the horizontal direction.
  • the upper symmetrical midpoint position 75U is a midpoint in the lateral direction Dx on the upper flange surface 33U, and is a position where the lower symmetrical position 71L and the position in the axial direction Dy match.
  • the displacement amount calculation unit 65 calculates the position in the vertical direction Dz indicated by the effective three-dimensional coordinate data of the lower target position 71L changed by the coordinate changing unit 63 and the vertical direction indicated by the effective three-dimensional coordinate data of the target contact position 73.
  • the difference from the position of Dz is defined as the displacement amount in the vertical direction Dz of the lower target position 71L, and the position in the vertical direction Dz indicated by the effective three-dimensional coordinate data of the upper target position 71U changed by the coordinate changing unit 63 and the target
  • the difference from the position in the vertical direction Dz indicated by the effective three-dimensional coordinate data of the contact position 73 is defined as the amount of displacement in the vertical direction Dz of the upper target position 71U.
  • the lower target position 71L is a position where the stationary component is arranged in the axial direction Dy and a position of an inner edge in the lower flange surface 33L. is.
  • the actual measurement coordinate receiving part 61 is the midpoint position in the lateral direction Dx in the lower flange surface 33L. a plurality of lower midpoint positions 75Lx whose positions in the axial direction Dy are different from each other; and the upper midpoint position 75Ux of the measured three-dimensional coordinate data.
  • the effective coordinate grasping unit 62 obtains the effective three-dimensional coordinate data at the lower target middle point position 75L from the change tendency of the measured three-dimensional coordinate data at the plurality of lower middle point positions 75Lx, and calculates the effective three-dimensional coordinate data at the plurality of upper middle point positions.
  • Effective three-dimensional coordinate data at the upper target midpoint position 75U is obtained from the change tendency of the effective three-dimensional coordinate data at 75Ux.
  • the measured coordinate reception unit 61 does not receive the measured three-dimensional coordinate data at the lower target midpoint position 75L or the upper target midpoint position 75U, effective cubic coordinates at these positions It is possible to grasp the original coordinate data.
  • the effective coordinate grasping unit 62 detects the lower first position from the change tendency of the effective three-dimensional coordinate data at the plurality of lower middle point positions 75Lx. 72La and the second lower position 72Lb are obtained, and from the change tendency of the effective three-dimensional coordinate data at the plurality of upper midpoint positions 75Ux, the first upper position 72Ua and the second upper position 72Ub Obtain valid three-dimensional coordinate data.
  • the measured coordinate reception unit 61 receives measured three-dimensional coordinate data at the lower first position 72La, the lower second position 72Lb, the upper first position 72Ua, and the upper second position 72Ub. Even if not, valid three-dimensional coordinate data at these positions can be grasped.

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

フランジ変位量推定方法では、第一被支持部の下フランジ面に連なる面中における下第一位置と、下フランジ面中の下対象中点位置と、上フランジ面に連なる面中で前記下第一位置と水平方向における位置が一致している上第一位置と、上フランジ面中の上対象中点位置と、における有効三次元座標データを把握する。前記下第一位置の有効三次元座標データと前記上第一位置の有効三次元座標データとが一致するよう、各位置における有効三次元座標データを変更する。座標変更後の下対象中点位置と上対象中点位置との上下方向における中間位置を対象接触位置とする。上対象位置と対象接触位置との上下方向の差を上対象位置の変位量とし、下対象位置と対象接触位置との上下方向の差を下対象位置の変位量とする。

Description

回転機械のフランジ変位量推定方法、この方法を実行するためのプログラム、及び、この方法を実行する装置
 本開示は、回転機械で、ロータの外周を覆う上半ケーシング及び下半ケーシングのフランジ面における変位量を推定する、フランジ変位量推定方法、この方法を実行するためのプログラム、及び、この方法を実行する装置に関する。
 本願は、2022年2月25日に、日本国に出願された特願2022-027442号に基づき優先権を主張し、この内容をここに援用する。
 蒸気タービン等の回転機械は、水平方向に延びる軸線を中心として回転可能なロータと、ロータの外周を覆うケーシングと、ケーシング内に配置され、このケーシングに取り付けられているダイヤフラム等の静止部品と、を備える。ケーシングは、一般的に、上側の上半ケーシングと、下側の下半ケーシングと、上半ケーシングと下半ケーシングとを締結する複数のボルトと、を有する。上半ケーシングは、下側を向く上フランジ面が形成されている上フランジを有する。下半ケーシングは、上側を向き、上フランジ面と上下方向で対向する下フランジ面が形成されている下フランジを有する。
 回転機械の点検時には、下半ケーシングから上半ケーシングを外した開放状態にして、回転機械を構成する複数の部品を点検、必要に応じて修理する。蒸気タービン等の回転機械におけるケーシングは、運転中の熱等の影響でクリープ変形等の非弾性変形が発生することがある。このため、一旦運転した後における開放状態での下半ケーシング及び上半ケーシングは、工場出荷時から厳密には変形している。点検が終了すると、複数の部品を組み立てる。この組立の工程には、複数のボルトを用いて、下半ケーシングに上半ケーシングを締結して締結状態にする工程が含まれる。下半ケーシング及び上半ケーシングを開放状態から締結状態にする過程で、下半ケーシング及び上半ケーシングはさらに変形する。
 ケーシングに取り付けられる静止部品とロータとの間の径方向の間隔は、予め定められた許容寸法の範囲内に収める必要がある。しかしながら、ケーシングが開放状態から締結状態になって、下半ケーシング及び上半ケーシングの形状が変わってしまうと、ケーシングに取り付けられた静止部品とロータとの間の径方向の間隔が変わり、この間隔が許容寸法の範囲から外れることがある。
 そこで、以下の特許文献1に記載の技術では、以下の手順で、開放状態から締結状態になった際の下半ケーシング及び上半ケーシングの変形量を推定している。まず、下半ケーシング及び上半ケーシングの三次元形状に関する有限要素モデルを取得する。続いて、開放状態における下半ケーシング及び上半ケーシングの三次元形状データを実測により取得する。次に、有限要素モデルが実測三次元形状データに合うよう、実測三次元形状データを用いて有限要素モデルを補正する。次に、開放状態を示す補正後の有限要素モデルを用いて、締結状態をシミュレーションして、締結状態を示す有限要素モデルを作成する。そして、開放状態を示す有限要素モデルと締結状態を示す有限要素モデルとの差から、下半ケーシング及び上半ケーシングにおける所定部位の変形量を推定する。なお、下半ケーシング及び上半ケーシングにおける所定部位とは、下半ケーシングの下フランジ面及び上半ケーシングの上フランジ面である。
 すなわち、特許文献1に記載の技術では、開放状態を示す有限要素モデルを用いて、締結状態をシミュレーションし、このシミュレーションで得られた締結状態を示す有限要素モデルから、下半ケーシングの下フランジ面及び上半ケーシングの上フランジ面の変位量を推定している。
特開2019-070334号公報
 特許文献1に記載の技術では、開放状態を示す有限要素モデルを用いて、締結状態をシミュレーションするため、このシミュレーションを実行するための計算負荷が大きい、という問題点がある。このため、特許文献1に記載の技術では、準備期間が長期化する上に、フランジ面の変位量の推定コストが嵩む、という問題点も内在している。
 そこで、本開示は、上半ケーシング及び下半ケーシングのフランジ面の変位量を推定するにあたり、計算負荷を抑えることで、フランジ面の推定の準備期間を短縮し且つその推定コストを抑えることができる技術を提供することを目的とする。
 前記目的を達成するための一態様としての回転機械のフランジ変位量推定方法は、以下の回転機械に適用される。
 この回転機械は、水平方向に延びる軸線を中心として回転可能なロータと、前記ロータの外周を覆うケーシングと、前記ケーシング内に配置され、前記ケーシングに取り付けられている静止部品と、前記ケーシングを下側から支える架台と、を備える。前記ケーシングは、上側の上半ケーシングと、下側の下半ケーシングと、前記上半ケーシングと前記下半ケーシングとを締結する複数のボルトと、を有する。前記上半ケーシングは、下側を向く上フランジ面が形成されている上フランジを有する。前記下半ケーシングは、上側を向き、前記上フランジ面と上下方向で対向する下フランジ面が形成されている下フランジと、前記下フランジに連なり、前記架台により下側から支えられ、前記軸線が延びる軸線方向で互いに離れた第一被支持部及び第二被支持部と、を有する。前記上フランジ及び前記下フランジには、上下方向に貫通して、前記複数のボルトのそれぞれが挿通可能なボルト孔が形成されている。
 以上の回転機械のフランジ変位量推定方法では、
 前記回転機械を分解した後で、前記上半ケーシングと前記下半ケーシングとが前記複数のボルトで締結されていない開放状態における、前記上フランジ面中の複数位置における実測三次元座標データ及び前記下フランジ面中の複数位置における実測三次元座標データを受け付ける実測座標受付工程と、前記下フランジ面中の複数位置における前記実測三次元座標データを用いて、下第一位置と、下第二位置と、下対象位置と、下対象中点位置と、における有効三次元座標データを把握すると共に、前記上フランジ面中の複数位置における前記実測三次元座標データを用いて、上第一位置と、上第二位置と、上対象位置と、上対象中点位置と、における有効三次元座標データを把握する有効座標把握工程と、前記有効座標把握工程で把握した前記下第一位置の有効三次元座標データと前記上第一位置の有効三次元座標データとが一致し、前記有効座標把握工程で把握した前記下第二位置の有効三次元座標データと前記上第二位置の有効三次元座標データとが一致するよう、前記有効座標把握工程で把握した有効三次元座標データを変更する座標変更工程と、前記座標変更工程で変更された、前記下対象中点位置及び前記上対象中点位置における有効三次元座標データを用いて、前記下対象中点位置と前記上対象中点位置との上下方向における中間位置である対象接触位置の有効三次元座標データを求める接触位置推定工程と、前記開放状態から前記上半ケーシングと前記下半ケーシングとが前記複数のボルトで締結された締結状態になったときの前記上対象位置及び前記下対象位置の上下方向の変位量を求める変位量演算工程と、を実行する。前記下第一位置は、前記上フランジ面に連なる面中で前記第一被支持部の第一代表位置と水平方向における位置が一致している位置である。前記下第二位置は、前記下フランジ面に連なる面中で前記第二被支持部の第二代表位置と水平方向における位置が一致している位置である。前記下対象位置は、前記下フランジ面中で、前記開放状態から前記締結状態になったときの上下方向の変位量を得たい位置である。前記下対象中点位置は、前記下フランジ面中で、水平方向で前記軸線方向に垂直な方向である横方向における中点の位置であって、前記下対象位置と前記軸線方向の位置が一致している位置である。前記上第一位置は、前記上フランジ面に連なる面中で前記第一被支持部の前記第一代表位置と水平方向における位置が一致している位置である。前記上第二位置は、前記上フランジ面に連なる面中で前記第二被支持部の前記第二代表位置と水平方向における位置が一致している位置である。前記上対象位置は、前記上フランジ面中で前記下対象位置と水平方向の位置が一致している位置である。前記上対象中点位置は、前記上フランジ面中で前記横方向における中点であって前記下対象位置と前記軸線方向の位置が一致している位置である。前記変位量演算工程では、前記座標変更工程で変更された前記下対象位置の有効三次元座標データが示す上下方向の位置と前記対象接触位置の有効三次元座標データが示す上下方向の位置との差を前記下対象位置の上下方向の変位量とし、前記座標変更工程で変更された前記上対象位置の有効三次元座標データが示す上下方向の位置と前記対象接触位置の有効三次元座標データが示す上下方向の位置との差を前記上対象位置の上下方向の変位量とする。
 本態様では、上フランジ面中の上対象中点位置と、下フランジ面中の下対象中点位置との上下方向の中点を対象接触位置とする。そして、本態様では、上フランジ面中で上下方向の変位量を得たい上対象位置と対象接触位置との上下方向の差を上対象位置の変位量とする。さらに、本態様では、下フランジ面中で上下方向の変位量を得たい下対象位置と対象接触位置との上下方向の差を下対象位置の変位量とする。このため、本態様では、下半ケーシング及び上半ケーシングの有限要素モデルを用いて、下半ケーシング及び上半ケーシングの変形をシミュレートしなくても、上対象位置及び下対象位置の上下方向の変位量を求めることができる。よって、本態様では、変位量を求める際の計算負荷を抑えることができる。
 ところで、上フランジ面中の上対象位置と、下フランジ面中の下対象位置との上下方向の中点を対象接触位置とすることも可能である。フランジ面の変形には、軸線方向の変化に伴う上下方向の変形のみならず、横方向の変化に伴う上下方向の変形も含まれる。仮に、下対象位置や上対象位置がフランジ面中の内側縁の位置であるとし、前述したように、上対象位置及び下対象位置を用いて対象接触位置を求めるとする。この場合、フランジ面における横方向の変化に伴う上下方向の変形が、求める対象接触位置に極端に反映されて、対象接触位置の上下方向の誤差が大きくなり、結果として、上対象位置及び下対象位置の変位量の誤差が大きくなる場合がある。一方、本態様では、上フランジ面中で横方向における中点である上対象中点位置と下フランジ面中で横方向における中点である下対象中点位置との上下方向の中点を対象接触位置にしている。このため、本態様では、フランジ面における横方向の変化に伴う上下方向の変形が、求める対象接触位置に極端に反映されることがなく、対象接触位置の上下方向の誤差を小さくでき、結果として、上対象位置及び下対象位置の変位量の誤差を小さくすることができる。
 前記目的を達成するための一態様としての回転機械のフランジ変位量推定プログラムは、以下の回転機械に適用される。
 この回転機械は、水平方向に延びる軸線を中心として回転可能なロータと、前記ロータの外周を覆うケーシングと、前記ケーシング内に配置され、前記ケーシングに取り付けられている静止部品と、前記ケーシングを下側から支える架台と、を備える。前記ケーシングは、上側の上半ケーシングと、下側の下半ケーシングと、前記上半ケーシングと前記下半ケーシングとを締結する複数のボルトと、を有する。前記上半ケーシングは、下側を向く上フランジ面が形成されている上フランジを有する。前記下半ケーシングは、上側を向き、前記上フランジ面と上下方向で対向する下フランジ面が形成されている下フランジと、前記下フランジに連なり、前記架台により下側から支えられ、前記軸線が延びる軸線方向で互いに離れた第一被支持部及び第二被支持部と、を有する。前記上フランジ及び前記下フランジには、上下方向に貫通して、前記複数のボルトのそれぞれが挿通可能なボルト孔が形成されている。
 以上の回転機械のフランジ変位量推定プログラムは、
 前記回転機械を分解した後で、前記上半ケーシングと前記下半ケーシングとが前記複数のボルトで締結されていない開放状態における、前記上フランジ面中の複数位置における実測三次元座標データ及び前記下フランジ面中の複数位置における実測三次元座標データを受け付ける実測座標受付工程と、前記下フランジ面中の複数位置における前記実測三次元座標データを用いて、下第一位置と、下第二位置と、下対象位置と、下対象中点位置と、における有効三次元座標データを把握すると共に、前記上フランジ面中の複数位置における前記実測三次元座標データを用いて、上第一位置と、上第二位置と、上対象位置と、上対象中点位置と、における有効三次元座標データを把握する有効座標把握工程と、前記有効座標把握工程で把握した前記下第一位置の有効三次元座標データと前記上第一位置の有効三次元座標データとが一致し、前記有効座標把握工程で把握した前記下第二位置の有効三次元座標データと前記上第二位置の有効三次元座標データとが一致するよう、前記有効座標把握工程で把握した有効三次元座標データを変更する座標変更工程と、前記座標変更工程で変更された、前記下対象中点位置及び前記上対象中点位置における有効三次元座標データを用いて、前記下対象中点位置と前記上対象中点位置との上下方向における中間位置である対象接触位置の有効三次元座標データを求める接触位置推定工程と、前記開放状態から前記上半ケーシングと前記下半ケーシングとが前記複数のボルトで締結された締結状態になったときの前記上対象位置及び前記下対象位置の上下方向の変位量を求める変位量演算工程と、をコンピュータに実行させる。前記下第一位置は、前記上フランジ面に連なる面中で前記第一被支持部の第一代表位置と水平方向における位置が一致している位置である。前記下第二位置は、前記下フランジ面に連なる面中で前記第二被支持部の第二代表位置と水平方向における位置が一致している位置である。前記下対象位置は、前記下フランジ面中で、前記開放状態から前記締結状態になったときの上下方向の変位量を得たい位置である。前記下対象中点位置は、前記下フランジ面中で、水平方向で前記軸線方向に垂直な方向である横方向における中点の位置であって、前記下対象位置と前記軸線方向の位置が一致している位置である。前記上第一位置は、前記上フランジ面に連なる面中で前記第一被支持部の前記第一代表位置と水平方向における位置が一致している位置である。前記上第二位置は、前記上フランジ面に連なる面中で前記第二被支持部の前記第二代表位置と水平方向における位置が一致している位置である。前記上対象位置は、前記上フランジ面中で前記下対象位置と水平方向の位置が一致している位置である。前記上対象中点位置は、前記上フランジ面中で前記横方向における中点であって前記下対象位置と前記軸線方向の位置が一致している位置である。前記変位量演算工程では、前記座標変更工程で変更された前記下対象位置の有効三次元座標データが示す上下方向の位置と前記対象接触位置の有効三次元座標データが示す上下方向の位置との差を前記下対象位置の上下方向の変位量とし、前記座標変更工程で変更された前記上対象位置の有効三次元座標データが示す上下方向の位置と前記対象接触位置の有効三次元座標データが示す上下方向の位置との差を前記上対象位置の上下方向の変位量とする。
 本態様では、このプログラムをコンピュータに実行させることにより、前記一態様におけるフランジ変位量推定方法と同様に、変位量を求める際の計算負荷を抑えることができる。
 前記目的を達成するための一態様としての回転機械のフランジ変位量装置は、以下の回転機械に適用される。
 この回転機械は、水平方向に延びる軸線を中心として回転可能なロータと、前記ロータの外周を覆うケーシングと、前記ケーシング内に配置され、前記ケーシングに取り付けられている静止部品と、前記ケーシングを下側から支える架台と、を備える。前記ケーシングは、上側の上半ケーシングと、下側の下半ケーシングと、前記上半ケーシングと前記下半ケーシングとを締結する複数のボルトと、を有する。前記上半ケーシングは、下側を向く上フランジ面が形成されている上フランジを有する。前記下半ケーシングは、上側を向き、前記上フランジ面と上下方向で対向する下フランジ面が形成されている下フランジと、前記下フランジに連なり、前記架台により下側から支えられ、前記軸線が延びる軸線方向で互いに離れた第一被支持部及び第二被支持部と、を有する。前記上フランジ及び前記下フランジには、上下方向に貫通して、前記複数のボルトのそれぞれが挿通可能なボルト孔が形成されている。
 以上の回転機械のフランジ変位量推定装置は、
 前記回転機械を分解した後で、前記上半ケーシングと前記下半ケーシングとが前記複数のボルトで締結されていない開放状態における、前記上フランジ面中の複数位置における実測三次元座標データ及び前記下フランジ面中の複数位置における実測三次元座標データを受け付ける実測座標受付部と、前記下フランジ面中の複数位置における前記実測三次元座標データを用いて、下第一位置と、下第二位置と、下対象位置と、下対象中点位置と、における有効三次元座標データを把握すると共に、前記上フランジ面中の複数位置における前記実測三次元座標データを用いて、上第一位置と、上第二位置と、上対象位置と、上対象中点位置と、における有効三次元座標データを把握する有効座標把握部と、前記有効座標把握部が把握した前記下第一位置の有効三次元座標データと前記上第一位置の有効三次元座標データとが一致し、前記有効座標把握部が把握した前記下第二位置の有効三次元座標データと前記上第二位置の有効三次元座標データとが一致するよう、前記有効座標把握部が把握した有効三次元座標データを変更する座標変更部と、前記座標変更部が変更した、前記下対象中点位置及び前記上対象中点位置における有効三次元座標データを用いて、前記下対象中点位置と前記上対象中点位置との上下方向における中間位置である対象接触位置の有効三次元座標データを求める接触位置推定部と、前記開放状態から前記上半ケーシングと前記下半ケーシングとが前記複数のボルトで締結された締結状態になったときの前記上対象位置及び前記下対象位置の上下方向の変位量を求める変位量演算部と、を備える。前記下第一位置は、前記上フランジ面に連なる面中で前記第一被支持部の第一代表位置と水平方向における位置が一致している位置である。前記下第二位置は、前記下フランジ面に連なる面中で前記第二被支持部の第二代表位置と水平方向における位置が一致している位置である。前記下対象位置は、前記下フランジ面中で、前記開放状態から前記締結状態になったときの上下方向の変位量を得たい位置である。前記下対象中点位置は、前記下フランジ面中で、水平方向で前記軸線方向に垂直な方向である横方向における中点の位置であって、前記下対象位置と前記軸線方向の位置が一致している位置である。前記上第一位置は、前記上フランジ面に連なる面中で前記第一被支持部の前記第一代表位置と水平方向における位置が一致している位置である。前記上第二位置は、前記上フランジ面に連なる面中で前記第二被支持部の前記第二代表位置と水平方向における位置が一致している位置である。前記上対象位置は、前記上フランジ面中で前記下対象位置と水平方向の位置が一致している位置である。前記上対象中点位置は、前記上フランジ面中で前記横方向における中点であって前記下対象位置と前記軸線方向の位置が一致している位置である。前記変位量演算部は、前記座標変更部が変更した前記下対象位置の有効三次元座標データが示す上下方向の位置と前記対象接触位置の有効三次元座標データが示す上下方向の位置との差を前記下対象位置の上下方向の変位量とし、前記座標変更部が変更した前記上対象位置の有効三次元座標データが示す上下方向の位置と前記対象接触位置の有効三次元座標データが示す上下方向の位置との差を前記上対象位置の上下方向の変位量とする。
 本態様では、前記一態様におけるフランジ変位量推定方法と同様に、変位量を求める際の計算負荷を抑えることができる。
 本開示の一態様では、上半ケーシング及び下半ケーシングのフランジ面の変位量を推定するにあたり、計算負荷を抑えることで、フランジ面の推定の準備期間を短縮し且つその推定コストを抑えることができる。
本開示に係る一実施形態における回転機械としての蒸気タービンの概略構成を示す模式図である。 本開示に係る一実施形態における回転機械としての蒸気タービンの概略外形を示す模式図である。 本開示に係る一実施形態における上半ケーシングの要部、及び下半ケーシングの要部の平面図である。 本開示に係る一実施形態における開放状態のケーシングの断面図である。 本開示に係る一実施形態における締結状態のケーシングの断面図である。 本開示に係る一実施形態におけるフランジ変位量推定装置の機能ブロック図である。 本開示に係る一実施形態におけるフランジ変位量推定方法の手順を示すフローチャートである。 本開示に係る一実施形態におけるフランジ面中で、有効三次元座標データを把握する位置を示す説明図である。 本開示に係る一実施形態における接触位置推定工程及び変位量演算工程での処理内容を示す説明図である。 本開示に係る一実施形態における座標変更工程S3での処理内容を示す説明図である。 本開示に係る一実施形態における第一把握方法を実行する場合に必要な実測三次元形状データの位置を示す説明図である。 本開示に係る一実施形態における第一把握方法を実行する場合の有効三次元座標データの求め方を示す説明図である。 本開示に係る一実施形態における第一把握方法を実行する場合の有効三次元座標データの他の求め方を示す説明図である。 本開示に係る一実施形態における第二把握方法を実行する場合に必要な実測三次元座標データの位置を示す説明図である。 本開示に係る一実施形態における第二把握方法を実行する場合の有効三次元座標データの求め方を示す説明図である。 本開示に係る一実施形態における第三把握方法を実行する場合に必要な実測三次元座標データの位置を示す説明図である。 本開示に係る一実施形態における第三把握方法を実行する場合の、基準三次元形状データが示すフランジ面と、実際のフランジ面の複数の位置における実測三次元座標データが示す点との相対位置関係を示すイメージ図である。 本開示に係る一実施形態における複数のポリゴンデータを説明するための説明図である。 本開示に係る一実施形態における複数のポリゴンデータから特定の複数のポリゴンデータの抽出を説明するための説明図である。 本開示に係る一実施形態における第三把握方法を実行する場合の、基準三次元形状データが示すフランジ面と、実際のフランジ面の複数の位置のうち、ポリゴンデータ抽出後の複数の位置における実測三次元座標データが示す点との相対位置関係を示すイメージ図である。 本開示に係る一実施形態における第三把握方法を実行する場合の、参照位置の求め方を示す説明図である。 本開示に係る一実施形態における第四把握方法を実行する場合に必要な実測三次元座標データの位置を示す説明図である。 本開示に係る一実施形態における第四把握方法を実行する場合の、基準三次元形状データが示すフランジ面と、実際のフランジ面の複数の位置における実測三次元座標データが示す点との相対位置関係を示すイメージ図である。 本開示に係る一実施形態における第四把握方法を実行する場合の、基準三次元形状データが示すフランジ面と、実際のフランジ面の複数の位置のうち、ポリゴンデータ抽出後の複数の位置における実測三次元座標データが示す点との相対位置関係を示すイメージ図である。 本開示に係る一実施形態における第四把握方法を実行する場合の、参照位置の求め方を示す説明図である。
 以下、本開示に係る回転機械のフランジ変位量推定方法、この方法を実行するためのプログラム、及び、この方法を実行する装置の実施形態について説明する。
 「回転機械の実施形態」
 本実施形態における回転機械について、図1~図5を参照して説明する。
 図1~図3に示すように、本実施形態の回転機械は、蒸気タービン10である。この蒸気タービン10は、水平方向に延びる軸線Arを中心として回転するロータ15と、ロータ15の外周側を覆うケーシング30と、ロータ15を回転可能に支持する第一軸受装置12a及び第二軸受装置12bと、複数のダイヤフラム20と、ケーシング30とロータ15との隙間を封止する第一軸封装置13a及び第二軸封装置13bと、ケーシング30を下側から支持する架台11と、を備える。
 ここで、軸線Arが延びる方向を軸線方向Dy、水平方向で軸線方向Dyに垂直な方向を横方向Dx、軸線Arに対する周方向を単に周方向Dc、軸線Arに対する径方向を単に径方向Dr、とする。また、この径方向Drで、軸線Arに近づく側を径方向内側Dri、軸線Arから遠ざかる側を径方向外側Droとする。また、図中の符号で用いているUは上半を意味し、Lは下半を意味する。
 ロータ15は、軸線方向Dyに延びるロータ軸16と、軸線方向Dyに並んでロータ軸16に取り付けられている複数の動翼列17と、を有する。複数の動翼列17は、いずれも、軸線Arに対する周方向Dcに並ぶ複数の動翼を有する。ロータ軸16の両端部は、ケーシング30から軸線方向Dyに突出している。ロータ軸16で、軸線方向Dyにおける一方の端部は、架台11に取り付けられている第一軸受装置12aにより、回転可能に支持されている。ロータ軸16で、軸線方向Dyにおける他方の端部は、架台11に取り付けられている第二軸受装置12bにより、回転可能に支持されている。
 第一軸封装置13aは、ケーシング30の軸線方向Dyにおける一方の端部に設けられている。第二軸封装置13bは、ケーシング30の軸線方向Dyにおける他方の端部に設けられている。第一軸封装置13a及び第二軸封装置13bは、いずれも、ロータ軸16とケーシング30との隙間を封止する装置である。
 複数のダイヤフラム20は、ケーシング30内で軸線方向Dyに並んでいる。複数のダイヤフラム20は、いずれも、軸線Arよりも下側の部分を構成する下半ダイヤフラム20Lと、軸線Arよりも上側の部分を構成する上半ダイヤフラム20Uと、を有する。下半ダイヤフラム20L及び上半ダイヤフラム20Uは、いずれも、周方向Dcに並ぶ複数の静翼22と、複数の静翼22の径方向内側Driの部分を相互に連結するダイヤフラム内輪23と、複数の静翼22の径方向外側Droの部分を相互に連結するダイヤフラム外輪24と、ダイヤフラム内輪23の径方向内側Driに取り付けられているシール装置25と、を有する。このシール装置25は、ダイヤフラム内輪23とロータ軸16との間の隙間をシールするシール装置である。
 以上で説明した第一軸封装置13a及び第二軸封装置13b、さらに、複数のダイヤフラム20は、いずれも、軸線Arに対する周方向に延びて、ケーシング30に取り付けられている静止部品である。
 ケーシング30は、図2に示すように、軸線Arよりも下側の部分を構成する下半ケーシング30Lと、軸線Arよりも上側の部分を構成する上半ケーシング30Uと、下半ケーシング30Lに対して上半ケーシング30Uを締結するための複数のボルト39と、を有する。下半ケーシング30Lは、周方向Dcに延びる下半ケーシング本体31Lと、下半ケーシング本体31Lの周方向Dcの両端部から径方向外側Droに突出する下フランジ32Lと、下フランジ32Lに連なり架台11により下側から支えられる第一被支持部35a及び第二被支持部35bと、を有する。また、上半ケーシング30Uは、周方向Dcに延びる上半ケーシング本体31Uと、上半ケーシング本体31Uの周方向Dcの両端部から径方向外側Droに突出する上フランジ32Uと、を有する。なお、上フランジ32Uには、下フランジ32Lにおける第一被支持部35a及び第二被支持部35bに対向する部分が設けられていない。但し、上フランジ32Uに、下フランジ32Lにおける第一被支持部35a及び第二被支持部35bに対向する部分が設けられていてもよい。
 図2~図5に示すように、下フランジ32Lで上側を向く面が下フランジ面33Lを成す。また、上フランジ32Uで下側を向く面が上フランジ面33Uを成す。下フランジ面33Lと上フランジ面33Uとは、上下方向Dzで互いに対向している。
 第一被支持部35aは、下フランジ32Lの軸線方向Dyにおける両側のうちの一方側から一方側に突出している。第二被支持部35bは、下フランジ32Lの軸線方向Dyにおける両側のうちの他方側から他方側に突出している。よって、第一被支持部35aに対して、第二被支持部35bは、軸線方向Dyに離れている。本実施形態において、第一被支持部35aの上面35ap及び第二被支持部35bの上面35bpは、下フランジ面33Lに連なる面である。すなわち、第一被支持部35aの上面35ap及び第二被支持部35bの上面35bpは、下フランジ面33Lと連続しており、下フランジ面33Lに対して段差がない。
 下フランジ32L及び上フランジ32Uには、上下方向Dzに貫通して、複数のボルト39のそれぞれが挿通可能なボルト孔34が形成されている。下半ケーシング30Lと上半ケーシング30Uとは、下フランジ32Lのボルト孔34及び上フランジ32Uのボルト孔34に挿通されたボルト39により締結される。
 下半ケーシング本体31Lの内周面、及び上半ケーシング30Uの内周面には、前述した複数の静止部品がそれぞれ格納される複数の静止部品格納部36が形成されている。下半ケーシング本体31Lの各静止部品格納部36は、下半ケーシング本体31Lの内周面から径方向外側Droに凹み、周方向Dcに延びる溝である。また、上半ケーシング本体31Uの各静止部品格納部36は、上半ケーシング本体31Uの内周面から径方向外側Droに凹み、周方向Dcに延びる溝である。なお、静止部品の一種であるダイヤフラム20は、周方向Dcに延びる静止部品格納部36のうち、フランジ面近傍の部分で支持されている。
 ケーシング30の内周面は、蒸気タービン10の運転により、高温の蒸気に晒される。このため、ケーシング30は、蒸気タービン10の運転により、クリープ変形等の非弾性変形が発生することがある。この変形の結果、下半ケーシング30Lに対して上半ケーシング30Uが締結されていない開放状態では、図4に示すように、下フランジ面33L及び上フランジ面33Uの上下方向Dzの位置が、軸線方向Dyの位置に応じて変わる。
 以上のように変形した下半ケーシング30Lに、以上のように変形した上半ケーシング30Uを締結して、ケーシング30を締結状態にすると、図5に示すように、下フランジ面33L及び上フランジ面33Uの上下方向Dzの位置が、軸線方向Dyの位置に応じてさらに変わる。
 ケーシング30に取り付けられる静止部品とロータ15との間の径方向Drの間隔は、予め定められた許容寸法の範囲内に収める必要がある。具体的に、例えば、静止部品の一種である第一軸封装置13a及び第二軸封装置13bとロータ軸16との間の間隔や、ダイヤフラム20のシール装置25とロータ軸16との間の間隔は、予め定められた許容寸法の範囲内に収める必要がある。しかしながら、開放状態の下半ケーシング30Lの形状データ及び上半ケーシング30Uの形状データがあったとしても、ケーシング30が開放状態から締結状態になって、下半ケーシング30L及び上半ケーシング30Uの形状が変わってしまうと、静止部品とロータ15との間の径方向Drの間隔が変わり、この間隔が許容寸法の範囲から外れることがある。
 発明者は、開放状態から締結状態にしたことによる下半ケーシング30L及び上半ケーシング30Uの変形に伴う静止部品とロータ15との間の径方向Drの間隔の変化は、下フランジ面33L及び上フランジ面33Uの変形に対して支配的であることを見出した。そこで、発明者は、開放状態から締結状態にしたことによる下フランジ面33Lの変位量及び上フランジ面33Uの変位量の推定し、これらの変位量に基づき、締結状態のときの静止部品とロータ15との間の径方向Drの間隔を把握するようにした。
 以下、下フランジ面33Lの変位量及び上フランジ面33Uの変位量を推定するフランジ変位量推定装置、及びフランジ変位量推定方法について説明する。
 「フランジ変位量推定装置の実施形態」
 本実施形態におけるフランジ変位量推定装置について、図6を参照して説明する。
 フランジ変位量推定装置50は、コンピュータである。このフランジ変位量推定装置50は、各種演算を行うCPU(Central Processing Unit)60と、CPU60のワークエリア等になるメモリ57と、ハードディスクドライブ装置等の補助記憶装置58と、キーボードやマウス等の手入力装置(入力装置)51と、表示装置(出力装置)52と、手入力装置51及び表示装置52の入出力インタフェース53と、三次元レーザ計測器等の三次元形状測定装置69との間でデータの受送信を行うための装置インタフェース(入力装置)54と、ネットワークNを介して外部と通信するための通信インタフェース(入出力装置)55と、非一時的な記憶媒体の一種であるディスク型記憶媒体Dに対してデータの記憶処理や再生処理を行う記憶・再生装置(入出力装置)56と、を備えている。
 補助記憶装置58には、フランジ変位量推定プログラム58pや、蒸気タービン10を構成する複数の部品毎の基準三次元形状データ58dが予め格納されている。この基準三次元形状データ58dは、三次元設計データであってもよいし、例えば、蒸気タービン10を工場から出荷する前に実測で得た三次元データであってもよい。すなわち、この基準三次元形状データ58dは、定期検査前の運転よりも前に得られた三次元データであればよい。基準三次元形状データ58dからは、複数の部品毎の各位置における三次元座標データを得ることができる。フランジ変位量推定プログラム58pは、例えば、記憶・再生装置56を介して、非一時的記憶媒体の一種であるディスク型記憶媒体Dから補助記憶装置58に取り込まれる。なお、このフランジ変位量推定プログラム58pは、通信インタフェース55を介して外部の装置から補助記憶装置58に取り込まれてもよい。
 CPU60は、機能的に、実測座標受付部61と、有効座標把握部62と、座標変更部63と、接触位置推定部64と、変位量演算部65と、を有する。これらの各機能部61~65は、いずれも、CPU60が補助記憶装置58に格納されているフランジ変位量推定プログラム58pを実行することで機能する。これらの各機能部61~65における動作については、後述する。
 「フランジ変位量推定方法の実施形態」
 本実施形態におけるフランジ変位量推定方法について、図7に示すフローチャートに従って説明する。なお、このフランジ変位量推定方法は、前述したフランジ変位量推定装置により実行される。
 蒸気タービン10は、点検等を行う毎に、分解、組立が行われる。蒸気タービン10は、分解が完了した時点では、図4に示すように、上半ケーシング30Uが下半ケーシング30Lから外される。この結果、ケーシング30は、上半ケーシング30Uと下半ケーシング30Lとがボルト39により締結されていない開放状態になる。さらに、ロータ15、複数のダイヤフラム20、第一軸封装置13a及び第二軸封装置13bは、ケーシング30から外され、このケーシング30外に配置される。なお、蒸気タービン10の分解が完了した時点で、下半ケーシング30Lが架台11から外されていてもよいが、ここでは、下半ケーシング30Lが架台11に支持されているとする。
 作業者は、以上のように蒸気タービン10を分解し、ケーシング30が開放状態になると、三次元レーザ計測器等の三次元形状測定装置69を用いて、上フランジ面33U中の複数位置における三次元座標値及び下フランジ面33L中の複数位置における三次元座標値を測定する。そして、作業者は、上フランジ面33U中の複数位置における三次元座標値及び下フランジ面33L中の複数位置における三次元座標値を実測三次元座標データとして、三次元形状測定装置69からフランジ変位量推定装置50に転送させる。フランジ変位量推定装置50の実測座標受付部61は、上フランジ面33U中の複数位置における実測三次元座標データ及び下フランジ面33L中の複数位置における実測三次元座標データを受け付ける(実測座標受付工程S1)。
 本実施形態における三次元座標データは、水平方向に延びる軸線方向Dyの位置を示す座標値と、軸線方向Dyに垂直な上下方向Dzの位置を示す座標値と、軸線方向Dyに垂直で水平方向に延びる横方向Dxの位置を示す座標値と、を含む。
 実測座標受付部61が複数の実測三次元座標データを受け付けると、フランジ変位量推定装置50の有効座標把握部62は、複数の実測三次元座標データを用いて、図8に示すように、複数の下対象位置71Lと、下第一位置72Laと、下第二位置72Lbと、複数の下対象中点位置75Lと、複数の上対象位置71Uと、上第一位置72Uaと、上第二位置72Ubと、複数の上対象中点位置75Uと、における有効三次元座標データを把握する(有効座標把握工程S2)。ここで、有効三次元座標データとは、受け付けた複数の実測三次元座標データを元に算出された、仮想上の面も含めた下フランジ面33L及び上フランジ面33Uの面上の点の三次元座標データである。このデータは開放状態から締結状態にしたことによる下フランジ面33Lの変位量及び上フランジ面33Uの変位量を推定するために必要なデータである。この有効三次元座標データの把握方法については、後ほど、詳細に説明する。
 ここで、下第一位置72Laは、下フランジ面33Lに連なる面中で第一被支持部35aの第一代表位置74aと水平方向における位置が一致している位置である。第一代表位置74aは、第一被支持部35a中で最も大きな荷重がかかる位置である。下第二位置72Lbは、下フランジ面33Lに連なる面中で第二被支持部35bの第二代表位置74bと水平方向における位置が一致している位置である。第二代表位置74bは、第二被支持部35b中で最も大きな荷重がかかる位置である。なお、「下フランジ面33Lに連なる面」は、実際に存在する面でも、仮想上の面であってもよい。本実施形態では、第一被支持部35aの上面35ap及び第二被支持部35bの上面35bpは、下フランジ面33Lに連なる面である。複数の下対象位置71Lは、下フランジ面33L中でケーシング30が開放状態から締結状態になったときの上下方向Dzの変位量を得たい位置である。ここで、下フランジ面33L中で上下方向Dzの変位量を得たい位置とは、下フランジ面33L中で軸線方向Dyにおける静止部品格納部36が形成されている位置であって、下フランジ面33L中の内側縁の位置である。下対象中点位置75Lは、図8及び図9に示すように、下フランジ面33L中で横方向Dxにおける中点の位置であって、下対象位置71Lと軸線方向Dyの位置が一致している位置である。上第一位置72Uaは、上フランジ面33Uに連なる面中で第一被支持部35aの第一代表位置74aと水平方向における位置が一致している位置である。上第二位置72Ubは、上フランジ面33Uに連なる面中で第二被支持部35bの第二代表位置74bと水平方向における位置が一致している位置である。なお、「上フランジ面33Uに連なる面」は、実際に存在する面でも、仮想上の面であってもよい。複数の上対象位置71Uは、上フランジ面33U中でケーシング30が開放状態から締結状態になったときの上下方向Dzの変位量を得たい位置である。ここで、上フランジ面33U中で上下方向Dzの変位量を得たい位置とは、上フランジ面33U中で軸線方向Dyにおける静止部品格納部36が形成されている位置であって、上フランジ面33U中の内側縁の位置である。上対象中点位置75Uは、図8及び図9に示すように、上フランジ面33U中で横方向Dxにおける中点の位置であって、上対象位置と軸線方向Dyの位置が一致している位置である。なお、各位置における有効三次元座標データの把握方法については、後ほど、詳細に説明する。
 複数の上対象位置71Uは、いずれも、複数の下対象位置71Lのうちのいずれか一の下対象位置71Lと水平方向の位置が一致している。従って、上対象中点位置75Uは、上フランジ面33U中で横方向Dxにおける中点の位置であって、下対象位置71Lと軸線方向Dyの位置が一致している位置でもある。ここで、XX方向の位置が一致しているとは、XX方向の位置が完全に一致している場合のみならず、XX方向の位置が実質的に一致している場合も含まれる。例えば、水平方向の位置が一致しているとは、軸線方向Dyの位置を示す座標値が同じで且つ横方向Dxの位置を示す座標値も同じであるという意味でのみならず、軸線方向Dyの位置を示す座標値が実質的に同じで且つ横方向Dxの位置を示す座標値も実質的に同じであるという意味も含まれる。
 開放状態から締結状態にしたことによる下半ケーシング30L及び上半ケーシング30Uの変形に伴う静止部品とロータ15との間の径方向Drの間隔の変化は、下フランジ面33L中で軸線方向Dyにおける静止部品格納部36が形成されている位置であって下フランジ面33L中の内側縁の位置の変形、及び、上フランジ面33U中で軸線方向Dyにおける静止部品格納部36が形成されている位置であって上フランジ面33U中の内側縁の位置の変形に対して支配的である。このため、上下方向Dzの変位量を得たい下対象位置71Lを前述の位置にし、上下方向Dzの変位量を得たい上対象位置71Uを前述の位置にしている。
 なお、下対象位置71Lは、下フランジ面33Lの内側縁の位置でなくてもよく、例えば、フランジ幅方向Dwで、下フランジ面33Lの内側縁からフランジ幅の1/3の位置までの範囲内の、いずれかの位置であってもよい。同様に、上対象位置71Uは、上フランジ面33Uの内側縁の位置でなくてもよく、例えば、フランジ幅方向Dwで、上フランジ面33Uの内側縁からフランジ幅の1/3の位置までの範囲内の、いずれかの位置であってもよい。ここで、フランジ幅方向Dwとは、フランジ面に沿って、フランジ面の外縁と内縁とを結ぶ方向であって、参照位置からフランジ面の外縁又は内縁までの距離が最も短くなる方向である。なお、参照位置とは、上対象位置71U、及び下対象位置71Lのそれぞれである。
 次に、フランジ変位量推定装置50の座標変更部63が、有効座標把握部62で把握された有効三次元座標データを変更する(座標変更工程S3)。具体的に、座標変更部63は、図10に示すように、下第一位置72Laの有効三次元座標データと上第一位置72Uaの有効三次元座標データとが一致し、下第二位置72Lbの有効三次元座標データと上第二位置72Ubの有効三次元座標データとが一致するよう、平行移動及び/又は回転移動などの座標変換により、有効座標把握部62で把握された有効三次元座標データを変更する。
 次に、フランジ変位量推定装置50の接触位置推定部64が、座標変更部63により変更された下対象中点位置75L及び上対象中点位置75Uにおける有効三次元座標データを用いて、下対象中点位置75Lと上対象中点位置75Uとの上下方向における中間位置である対象接触位置73の有効三次元座標データを求める(接触位置推定工程S4)。
 次に、フランジ変位量推定装置50の変位量演算部65が、ケーシング30が開放状態から締結状態になったときの上対象位置71U及び下対象位置71Lの上下方向Dzの変位量を求め、外部からの要求に応じて、これらの変位量を出力する(変位量演算工程S5)。具体的に、変位量演算部65は、図9及び以下の式に示すように、座標変更後の下対象位置71Lの有効三次元座標データに含まれる上下方向Dzの座標値ZLと、この下対象位置71Lに対する対象接触位置73の有効三次元座標データに含まれる上下方向Dzの座標値ZCと差を下対象位置71Lの上下方向Dzの変位量ZdLとする。さらに、変位量演算部65は、座標変更後の上対象位置71Uの有効三次元座標データに含まれる上下方向Dzの座標値ZUと、この上対象位置71Uに対する対象接触位置73の有効三次元座標データに含まれる上下方向Dzの座標値ZCと差を上対象位置71Uの上下方向Dzの変位量ZdUとする。
  ZdL=ZL-ZC
  ZdU=ZU-ZC
 以上で、フランジ変位量推定装置50による、下フランジ32Lの下対象位置71L、及び上対象位置71Uの下対象位置71Lにおける上下方向Dzの変位量の推定が終了する。
 次に、有効座標把握部62における有効三次元座標データの複数種類の把握方法について説明する。
 「第一把握方法」
 有効座標把握工程S2で第一把握方法を実行する場合、実測座標受付工程S1では、図11に示すように、複数の下対象位置71Lと、複数の下中点位置75Lxと、複数の上対象位置71Uと、複数の上中点位置75Uxと、における実測三次元座標データを受け付ける。ここで、下中点位置75Lxは、下フランジ面33L中で横方向Dxにおける中点の位置である。複数の下中点位置75Lxは、軸線方向Dyにおける位置が互いに異なっている。複数の下中点位置75Lxのうち、いずれかは、下対象中点位置75Lである。上中点位置75Uxは、上フランジ面33U中で横方向Dxにおける中点の位置である。複数の上中点位置75Uxは、軸線方向Dyにおける位置が互いに異なっている。複数の上中点位置75Uxのうち、いずれかは、上対象中点位置75Uである。
 第一把握方法における有効座標把握工程S2では、有効座標把握部62が、実測座標受付工程S1で受け付けた複数の下対象位置71L、及び複数の上対象位置71Uにおける実測三次元座標データを、そのまま、複数の下対象位置71L、及び複数の上対象位置71Uにおける有効三次元座標データとする。
 有効座標把握部62は、補助記憶装置58に記憶されている基準三次元形状データ58dから、このデータが作成された時点における複数の下対象位置71L、及び複数の上対象位置71Uの三次元座標データを得ることができる。そこで、有効座標把握部62は、例えば、以下のようにして、実測座標受付部61が受け付けた複数の位置における実測三次元座標データのうちから、一の下対象位置71Lの実測三次元座標データを認識する。有効座標把握部62は、実測座標受付部61が受け付けた複数の位置における実測三次元座標データのうちから、基準三次元形状データ58dが示す一の下対象位置71Lの三次元座標データと水平方向の座標値が一致している実測三次元座標データを抽出し、この実測三次元座標データを一の下対象位置71Lの実測三次元座標データを認識する。
 有効座標把握工程S2では、有効座標把握部62が、複数の下中点位置75Lxにおける実測三次元座標データの変化傾向から、下対象中点位置75Lにおける有効三次元座標データを求める。有効座標把握部62は、さらに、複数の上中点位置75Uxにおける有効三次元座標データの変化傾向から、上対象中点位置75Uにおける有効三次元座標データを求める。有効座標把握部62が下対象中点位置75Lにおける有効三次元座標データを求める場合、図12に示すように、複数の下中点位置75Lxにおける有効三次元座標データを用いて、下フランジ面33L中で複数の下中点位置75Lxに沿った領域における面形状を近似的に示す二次関数等の高次関数Fを求める。有効座標把握部62は、この高次関数Fを用いて、基準三次元形状データ58dが示す下対象中点位置75Lの水平方向の座標値に対する、上下方向Dzの座標値を求める。そして、有効座標把握部62は、基準三次元形状データ58dが示す下対象中点位置75Lに関する各方向の座標値のうち、上下方向Dzの座標値を、先に求めた上下方向Dzの座標値に置き換えて、これを下対象中点位置75Lの有効三次元座標データとする。
 実測座標受付部61が、複数の下対象中点位置75Lの全てにおける実測三次元座標データ、及び複数の上対象中点位置75Uの全てにおける実測三次元座標データを受け付けることができれば、これらの実測三次元座標データをそのまま有効三次元座標データにしてもよい。しかしながら、現実には、例えば、複数の下対象中点位置75Lのうち、いずれかの下対象中点位置75Lが、ボルト孔34等の位置である場合がある。この場合には、この下対象中点位置75Lの実測三次元座標データを得ることができない。このため、ここでは、複数の下中点位置75Lxにおける実測三次元座標データの変化傾向から、下対象中点位置75Lにおける有効三次元座標データを求めるようにしている。
 さらに、有効座標把握工程S2では、有効座標把握部62が、複数の下中点位置75Lxにおける実測三次元座標データの変化傾向を示す前述の高次関数Fを用いて、下第一位置72La及び下第二位置72Lbにおける有効三次元座標データを外挿により求める。有効座標把握部62が、さらに、複数の上中点位置75Uxにおける有効三次元座標データの変化傾向を示す高次関数を用いて、上第一位置72Ua及び上第二位置72Ubにおける有効三次元座標データを外挿により求める。
 以上では、高次関数Fで下フランジ面33L及び上フランジ面33Uの面形状を近似する。しかしながら、図13に示すように、一次関数で、下フランジ面33Lの一部の面形状や、上フランジ面33Uの一部の面形状を近似してもよい。この場合、有効座標把握部62は、複数の下中点位置75Lxのうちで、有効三次元座標データを把握すべき下対象中点位置75Lに近い複数の下中点位置75Lxにおける有効三次元座標データを用いて、有効三次元座標データを把握すべき下対象中点位置75Lを一次関数Fa,Fbで近似する。そして、この一次関数Fa,Fbを用いて、この下対象中点位置75Lにおける有効三次元座標データの上下方向Dzの座標値を求める。また、有効座標把握部62は、複数の上中点位置75Uxのうちで、有効三次元座標データを把握すべき上対象中点位置75Uに近い複数の上中点位置75Uxにおける有効三次元座標データを用いて、有効三次元座標データを把握すべき上対象中点位置75Uを一次関数で近似する。そして、この一次関数を用いて、この上対象中点位置75Uにおける有効三次元座標データの上下方向Dzの座標値を求める。さらに、有効座標把握部62は、前述の一次関数Fa,Fbを用いて、下第一位置72La及び下第二位置72Lbにおける有効三次元座標データを外挿により求めると共に、一次関数を用いて、上第一位置72Ua及び上第二位置72Ubにおける有効三次元座標データを外挿により求める。
 以上で、複数の下対象位置71Lと、複数の下対象中点位置75Lと、下第一位置72Laと、下第二位置72Lbと、複数の上対象位置71Uと、複数の上対象中点位置75Uと、上第一位置72Uaと、上第二位置72Ubと、における有効三次元座標データが把握される。
 以上のように、第一把握方法では、取り扱う三次元座標データの数を少なくすることができるため、作業者が三次元座標値を測定する手間を軽減することができる上に、コンピュータによる計算負荷を軽減することができる。
 以上では、下第一位置72La及び下第二位置72Lbの上下方向Dzの座標値を推定する。しかしながら、実測座標受付工程S1で、下フランジ面33Lと連なっている第一被支持部35aの上面35apにおける下第一位置72Laの実測三次元座標データ、さらに、下フランジ面33Lと連なっている第二被支持部35bの上面35bpにおける下第二位置72Lbの実測三次元座標データを受け付けた場合には、これらの実測三次元データのそれぞれを、そのまま、下第一位置72Laの有効三次元座標データ、下第二位置72Lbの有効三次元座標データとしてもよい。
 前述したように、上半ケーシング30Uも、上フランジ32Uに連なる第一被支持部及び第二被支持部を有する場合がある。この場合に、実測座標受付工程S1で、上フランジ面33Uと連なっている第一被支持部の下面における上第一位置72Uaの実測三次元座標データ、さらに、上フランジ面33Uと連なっている第二被支持部の下面における上第二位置72Ubの実測三次元座標データを受け付けると、これらの実測三次元データのそれぞれを、そのまま、上第一位置72Uaの有効三次元座標データ、上第二位置72Ubの有効三次元座標データとしてもよい。
 「第二把握方法」
 有効座標把握工程S2で第二把握方法を実行する場合、実測座標受付工程S1では、図14に示す以下の位置における実測三次元座標データを受け付ける。
a.複数の下対象位置71L毎に、下対象位置71Lを通り且つフランジ幅方向Dwに延びる下仮想線76L上の複数の位置78における実測三次元座標データ
b.複数の上対象位置71U毎に、上対象位置71Uを通り且つフランジ幅方向Dwに延びる上仮想線76U上の複数の位置78における実測三次元座標データ
 ここで、「YY位置を通り且つフランジ幅方向に延びるZ仮想緯線」における「YY位置」は、基準三次元形状データ58dが示す「YY位置」である。また、実測座標受付工程S1で実測三次元座標データを受け付ける仮想線上の位置の数は、例えば、2以上で10未満である。
 第二把握方法における有効座標把握工程S2では、実測座標受付工程S1で受け付けた複数の実測三次元座標データを用いて、有効座標把握部62が、複数の下対象位置71L及び複数の上対象位置71Uにおける有効三次元座標データを把握する。
 有効座標把握部62は、図15に示すように、参照位置71を通り且つフランジ幅方向Dwに延びる仮想線76上の複数の位置78における実測三次元座標データを用いて、仮想線76上の複数の位置78における上下方向Dzの座標値を近似的に示す関数F2を求める。有効座標把握部62は、この関数F2を用いて、参照位置71における上下方向Dzの座標値を外挿により求める。そして、有効座標把握部62は、基準三次元形状データ58dが示す参照位置71に関する各方向の座標値のうち、上下方向Dzの座標値を、先に求めた上下方向Dzの座標値に置き換えて、これを参照位置71の有効三次元座標データとする。
 さらに、有効座標把握工程S2で第二把握方法を実行する場合、実測座標受付工程S1では、以下の位置における実測三次元座標データを受け付ける。
c.複数の下対象中点位置75L毎に、下対象中点位置75Lを通り且つフランジ幅方向Dwに延びる下中点仮想線上の複数の位置における実測三次元座標データ
d.複数の上対象中点位置75U毎に、上対象中点位置75Uを通り且つフランジ幅方向Dwに延びる上中点仮想線上の複数の位置における実測三次元座標データ
 有効座標把握工程S2では、実測座標受付工程S1で受け付けた、下中点仮想線上の複数の位置における実測三次元座標データ及び上中点仮想線上の複数の位置における実測三次元座標データを用いて、有効座標把握部62が、前述の方法と同様の方法で、前述の参照位置としての、複数の下対象中点位置75L及び複数の上対象中点位置75Uにおける有効三次元座標データを把握する。
 さらに、有効座標把握工程S2では、有効座標把握部62が、複数の下対象中点位置75Lにおける有効三次元座標データを用いて、下第一位置72La及び下第二位置72Lbにおける有効三次元座標データを推定する。さらに、有効座標把握部62は、複数の上対象中点位置75Uにおける有効三次元座標データを用いて、上第一位置72Ua及び上第二位置72Ubにおける有効三次元座標データを推定する。有効座標把握部62が下第一位置72La及び下第二位置72Lbにおける有効三次元座標データを推定する際、有効座標把握部62は、第一把握方法で説明した方法と同様、複数の下対象中点位置75Lにおける実測三次元座標データの変化傾向から、下第一位置72La及び下第二位置72Lbにおける有効三次元座標データを求める。また、有効座標把握部62が上第一位置72Ua及び上第二位置72Ubにおける有効三次元座標データを推定する際、有効座標把握部62は、第一把握方法で説明した方法と同様、複数の上対象中点位置75Uにおける実測三次元座標データの変化傾向から、上第一位置72Ua及び上第二位置72Ubにおける有効三次元座標データを求める。
 以上で、複数の下対象位置71Lと、複数の下対象中点位置75Lと、下第一位置72Laと、下第二位置72Lbと、複数の上対象位置71Uと、複数の上対象中点位置75Uと、上第一位置72Uaと、上第二位置72Ubと、における有効三次元座標データが把握される。
 第一把握方法では、参照位置の実測三次元座標データをそのままこの参照位置の有効三次元座標データにしている。このため、参照位置の有効三次元座標データには、局所的な形状変化の影響を受け易い上に、大きな計測誤差が含まれる可能性がある。例えば、三次元形状測定装置69が三次元レーザ計測器である場合、計測対象と三次元レーザ計測器との間に、微小な浮遊物が存在すると、この三次元レーザ計測器で計測された三次元位置データには誤差が含まれることになる。一方、第二把握方法では、複数の位置における実測三次元座標データから参照位置71の三次元座標データを推定し、この三次元座標データを有効三次元座標データにしている。このため、第二把握方法では、第一把握方法よりも、局所的な形状変化の影響を受けにくくなる上に、大きな計測誤差が含まれる可能性を抑えることができる。
 以上では、下第一位置72La及び下第二位置72Lbの上下方向Dzの座標値を推定する。しかしながら、実測座標受付工程S1で、以下の実測三次元座標データを受け付けた場合には、前述した参照位置71の有効三次元座標データを求めた方法と同様の方法で、この実測三次元座標データを用いて、下第一位置72La及び下第二位置72Lbの有効三次元座標データを求めてもよい。
a.下第一位置72Laを通り且つ第一被支持部35aの突出方向に延びる仮想線上の複数の位置における実測三次元座標データ
b.下第二位置72Lbを通り且つ第二被支持部35bの突出方向に延びる仮想線上の複数の位置における実測三次元座標データ
 ここで、被支持部35a,35bの突出方向とは、被支持部35a,35bの上面35ap,35bpに沿って、フランジから被支持部35a,35bが突出する方向である。
 前述したように、上半ケーシング30Uも、上フランジ32Uに連なる第一被支持部及び第二被支持部を有する場合がある。この場合に、実測座標受付工程S1で、以下の実測三次元座標データを受け付けた場合には、前述した参照位置の有効三次元座標データを求めた方法と同様の方法で、この実測三次元座標データを用いて、上第一位置72Ua及び上第二位置72Ubの有効三次元座標データを求めてもよい。
a.上第一位置72Uaを通り且つ第一被支持部の突出方向に延びる仮想線上の複数の位置における実測三次元座標データ
b.上第二位置72Ubを通り且つ第二被支持部の突出方向に延びる仮想線上の複数の位置における実測三次元座標データ
 「第三把握方法」
 有効座標把握工程S2で第三把握方法を実行する場合、図16及び図17に示すように、実測座標受付工程S1では、下フランジ面33Lの全体に渡る複数の位置78における実測三次元座標データ、及び上フランジ面33Uの全体に渡る複数の位置における実測三次元座標データを受け付ける。なお、図17は、基準三次元形状データ58dが示すフランジ面80及び参照位置81と、実際のフランジ面の全体に渡る複数の位置における実測三次元座標データが示す点85との相対位置関係を示すイメージ図である。
 第三把握方法における有効座標把握工程S2では、有効座標把握部62が、まず、図18に示すように、フランジ面の全体に渡る複数の位置における実測三次元座標データを用いて、複数のポリゴンデータを作成する。ポリゴンデータとは、多角形の平面を規定するデータである。有効座標把握部62は、複数の位置における実測三次元座標データが示す点85のうち、互に近接する複数の点85を線分で結び、これらの線分で囲まれた多角形平面をポリゴン86とする。
 有効座標把握部62は、次に、複数のポリゴンデータのうちから、図19に示すように、ある条件を満たす複数のポリゴンデータを抽出する。なお、図19では、抽出するポリゴンデータで特定されるポリゴン86aに模様を施し、抽出しないポリゴンデータで特定されるポリゴン86bには模様を施していない。また、図19中のXY平面は、基準三次元形状データ58dが示すフランジ面80に平行な面である。ここで、前述の条件とは、基準三次元形状データ58dが示すフランジ面80に対する、ポリゴンデータで特定されるポリゴン86の傾きが所定の傾き以内である、である。有効座標把握部62は、まず、複数のポリゴン86毎に、ポリゴン86の法線nを求める。次に、有効座標把握部62は、複数のポリゴン86毎に、基準三次元形状データ58dが示すフランジ面80に対する垂線pとポリゴン86の法線nとの角度αを求める。そして、有効座標把握部62は、複数のポリゴンデータのうちから、フランジ面80に対する垂線pとポリゴン86の法線nとの角度αが所定の角度(所定の傾き)以内の複数のポリゴンデータを抽出する。
 このデータの抽出処理は、実測座標受付工程S1で受け付けた複数の点85における実測三次元座標データから、フランジ面の縁の壁中の点や、フランジ面を貫通するボルト孔34の内周面中の点における実測三次元座標データを除くために実行される。このため、この抽出処理後の点85の数は、図20に示すように、その前の点85の数より少なくなる。特に、基準三次元形状データ58dが示す基準形状モデル中で、フランジ面80に対して傾斜している面82に関して、抽出処理後の点85の数は、その前の点85の数より著しく少なくなる。
 有効座標把握部62は、次に、図21に示すように、フランジ面80を含む仮想三次元空間を複数の三次元ブロック83に分割する。そして、有効座標把握部62は、複数の三次元ブロック83毎に、対象とする三次元ブロック83中の代表点87を定める。具体的に、有効座標把握部62は、抽出処理で抽出された複数のポリゴンデータで特定されるポリゴン86aに含まれる複数の点85のうち、対象とする三次元ブロック83中に含まれる複数の点85の中央値となる点を、対象とする三次元ブロック83中の代表点87とする。
 なお、代表点87は、抽出処理で抽出された複数のポリゴンデータで特定されるポリゴン86aに含まれる複数の点85のローレンツ分布に基づくロバスト推定やバイウェイト推定により定めてもよい。
 有効座標把握部62は、複数の三次元ブロック83毎の代表点87を相互に補完面としての平面又は曲面で接続して、複数の三次元ブロック83毎の代表点87を含む補完面の面形状データを作成する。この面形状データは、フランジ面全体の形状を示す関数F3で表される。有効座標把握部62は、関数F3で表されるフランジ面全体の面形状データを用いて、前述の参照位置71における有効三次元座標データを求める。なお、第三把握方法における参照位置71は、複数の下対象位置71L、複数の下対象中点位置75L、複数の上対象位置71U、複数の上対象中点位置75Uのそれぞれである。
 第三把握方法における有効座標把握工程S2でも、第一把握方法及び第二把握方法と同様に、有効座標把握部62が、複数の下対象中点位置75Lにおける有効三次元座標データを用いて、下第一位置72La及び下第二位置72Lbにおける有効三次元座標データを推定する。さらに、有効座標把握部62は、複数の上対象中点位置75Uにおける有効三次元座標データを用いて、上第一位置72Ua及び上第二位置72Ubにおける有効三次元座標データを推定する。
 以上で、複数の下対象位置71Lと、複数の下対象中点位置75Lと、下第一位置72Laと、下第二位置72Lbと、複数の上対象位置71Uと、複数の上対象中点位置75Uと、上第一位置72Uaと、上第二位置72Ubと、における有効三次元座標データが把握される。
 第三把握方法では、第二把握方法よりも、局所的な形状変化の影響を受けにくくなる上に、大きな計測誤差が含まれる可能性を抑えることができる。さらに、第三把握方法では、障害物等に起因した広範囲なデータ欠損がある場合でも、参照位置における有効三次元座標データを把握できる。
 以上では、下第一位置72La及び下第二位置72Lbの上下方向Dzの座標値を推定する。しかしながら、実測座標受付工程S1で、第一被支持部35aの上面35apの全体に渡る複数の位置及び第二被支持部35bの上面35bpの全体に渡る複数の位置における実測三次元座標データを受け付けた場合には、以下の方法で、下第一位置72La及び下第二位置72Lbの有効三次元座標データを求めてもよい。具体的に、まず、第一被支持部35aの上面35apの全体に渡る複数の位置及び第二被支持部35bの上面35bpの全体に渡る複数の位置における実測三次元座標データを含めて、第一被支持部35aの上面35ap、第二被支持部35bの上面35bp、及び下フランジ面33Lの全体の面形状データを求める。そして、関数で表される面全体の面形状データを用いて、下第一位置72La及び下第二位置72Lbの有効三次元座標データを求める。
 前述したように、上半ケーシング30Uも、上フランジ32Uに連なる第一被支持部及び第二被支持部を有する場合がある。この場合に、実測座標受付工程S1で、上フランジ32Uに連なる第一被支持部の下面の全体に渡る複数の位置及び上フランジ32Uに連なる第二被支持部の下面の全体に渡る複数の位置における実測三次元座標データを受け付け、以下の方法で、上第一位置72Ua及び上第二位置72Ubの有効三次元座標データを求めてもよい。具体的に、まず、第一被支持部の下面の全体に渡る複数の位置及び第二被支持部の下面の全体に渡る複数の位置における実測三次元座標データを含めて、第一被支持部の下面、第二被支持部の下面、及び上フランジ面33Uの全体の面形状データを求める。そして、関数で表される面全体の面形状データを用いて、上第一位置72Ua及び上第二位置72Ubの有効三次元座標データを求める。
 「第四把握方法」
 有効座標把握工程S2で第四把握方法を実行する場合、図22及び図23に示すように、実測座標受付工程S1では、フランジ面中で、前述の参照位置71を含む参照計測領域79中の複数の位置78における実測三次元座標データを受け付ける。なお、図23は、基準三次元形状データ58dが示すフランジ面80と、実際のフランジ面中で参照計測領域79中の複数の位置における実測三次元座標データが示す点85との相対位置関係を示すイメージ図である。ここで、参照計測領域79とは、図22に示すように、例えば、参照位置71を起点として、この参照位置71におけるフランジ幅の1/20~1/2の距離の範囲内の領域である。よって、この参照計測領域79は、下フランジ面33L中で下対象位置71Lを含む下計測領域でもあり、上フランジ面33U中で上対象位置71Uを含む上計測領域でもある。さらに、この参照計測領域79は、下フランジ面33L中で下対象中点位置75Lを含む下中点計測領域でもあり、上フランジ面33U中で上対象中点位置75Uを含む上中点計測領域でもある。なお、ここでの参照位置71の三次元座標データは、基準三次元形状データ58dが示す参照位置の三次元座標データである。また、実測座標受付工程S1で受け付ける参照計測領域79内の実測三次元座標データの数は、例えば、10以上である。よって、第四把握方法の実測座標受付工程S1で受け付ける参照計測領域79内の実測三次元座標データの数は、第二把握方法の実測座標受付工程S1で受け付ける仮想線上の位置の実測三次元座標データの数より多い。
 第四把握方法における有効座標把握工程S2では、有効座標把握部62が、まず、第三把握方法で説明した方法と同様に、複数の位置78における実測三次元座標データを用いて、複数のポリゴンデータを作成し、複数のポリゴンデータのうちから、ある条件を満たす複数のポリゴンデータを抽出する。この結果、この抽出処理後の実測三次元座標データが示す点85の数は、図24に示すように、その前の点85の数より少なくなる。
 有効座標把握部62は、次に、第三把握方法で説明した方法と同様に、図25に示すように、フランジ面80を含む仮想三次元空間を複数の三次元ブロック83に分割する。そして、有効座標把握部62が、複数の三次元ブロック83毎に、対象とする三次元ブロック83中の代表点87を定める。
 有効座標把握部62は、複数の三次元ブロック83毎の代表点87を相互に補完面としての平面又は曲面で接続して、複数の三次元ブロック83毎の代表点87を含む補完面の面形状データを作成する。この面形状データは、フランジ面中の参照計測領域79内の形状を示す関数F4で表される。有効座標把握部62は、関数F4で表される面形状データを用いて、前述の参照位置71における有効三次元座標データを求める。
 第四把握方法における二次処理工程でも、第一把握方法及び第二把握方法における二次処理工程と同様に、有効座標把握部62が、複数の下対象中点位置75Lにおける有効三次元座標データを用いて、下第一位置72La及び下第二位置72Lbにおける有効三次元座標データを推定する。さらに、有効座標把握部62は、複数の上対象中点位置75Uにおける有効三次元座標データを用いて、上第一位置72Ua及び上第二位置72Ubにおける有効三次元座標データを推定する。
 以上で、複数の下対象位置71Lと、複数の下対象中点位置75Lと、下第一位置72Laと、下第二位置72Lbと、複数の上対象位置71Uと、複数の上対象中点位置75Uと、上第一位置72Uaと、上第二位置72Ubと、における有効三次元座標データが把握される。
 第四把握方法では、第二把握方法よりも、局所的な形状変化の影響を受けにくくなる上に、大きな計測誤差が含まれる可能性を抑えることができる。さらに、第四把握方法では、障害物等に起因した広範囲なデータ欠損がある場合でも、参照位置における有効三次元座標データを把握できる。
 以上では、フランジ面中の参照計測領域79内の面形状データを用いて、前述の参照位置71における有効三次元座標データを求める。しかしながら、面形状データを作成せず、複数の三次元ブロックのうち、参照位置71を含む三次元ブロック83の代表点87における上下方向Dzの座標値を、参照位置71における上下方向Dzの座標値としてもよい。
 以上では、下第一位置72La及び下第二位置72Lbの上下方向Dzの座標値を推定する。しかしながら、実測座標受付工程S1で、第一被支持部35aの上面35apの全体に渡る複数の位置及び第二被支持部35bの上面35bpの全体に渡る複数の位置における実測三次元座標データを受け付けた場合には、以下の方法で、下第一位置72La及び下第二位置72Lbの有効三次元座標データを求めてもよい。具体的に、まず、第一被支持部35aの上面35apの全体に渡る複数の位置における実測三次元座標データを用いて、複数のポリゴンデータを作成し、複数のポリゴンデータのうちから、ある条件を満たす複数のポリゴンデータを抽出する。そして、抽出したポリゴンデータが示す複数の点のうちから代表点を定め、この代表点における上下方向Dzの座標値を、下第一位置72Laにおける上下方向Dzの座標値とする。同様に、第二被支持部35bの上面35apの全体に渡る複数の位置における実測三次元座標データを用いて、代表点を定め、この代表点における上下方向Dzの座標値を、下第二位置72Lbにおける上下方向Dzの座標値とする。
 前述したように、上半ケーシング30Uも、上フランジ32Uに連なる第一被支持部及び第二被支持部を有する場合がある。この場合に、実測座標受付工程S1で、上フランジ32Uに連なる第一被支持部の下面の全体に渡る複数の位置及び上フランジ32Uに連なる第二被支持部35bの下面の全体に渡る複数の位置における実測三次元座標データを受け付ける。そして、以上と同様に、受け付けた複数の位置における実測三次元座標データを用いて、各面上の代表点を定め、各面上の代表点における上下方向Dzの座標値を、それぞれ、上第一位置72Uaにおける上下方向Dzの座標値、上第二位置72Ubにおける上下方向Dzの座標値とする。
 なお、上対象中点位置75U及び下対象中点位置75Lにおける有効三次元座標データを第一把握方法で把握し、上対象位置71U及び下対象位置71Lにおける有効三次元座標データを第二把握方法又は第四把握方法で把握してもよい。また、上対象中点位置75U及び下対象中点位置75Lにおける有効三次元座標データを第二把握方法で把握し、上対象位置71U及び下対象位置71Lにおける有効三次元座標データを第一把握方法又は第四把握方法で把握してもよい。さらに、上対象中点位置75U及び下対象中点位置75Lにおける有効三次元座標データを第四把握方法で把握し、上対象位置71U及び下対象位置71Lにおける有効三次元座標データを第一把握方法又は第二把握方法で把握してもよい。
 以上のように、本実施形態では、上フランジ面33U中の上対象中点位置75Uと下フランジ面33L中の下対象中点位置75Lとの上下方向Dzの中点を対象接触位置73とする。そして、本実施形態では、上フランジ面33U中で上下方向Dzの変位量を得たい上対象位置71Uと対象接触位置73との上下方向Dzの差を上対象位置71Uの変位量とする。さらに、本実施形態では、下フランジ面33L中で上下方向Dzの変位量を得たい下対象位置71Lと対象接触位置73との上下方向Dzの差を下対象位置71Lの変位量とする。このため、本実施形態では、下半ケーシング30L及び上半ケーシング30Uの有限要素モデルを用いて、下半ケーシング30L及び上半ケーシング30Uの変形をシミュレートしなくても、上対象位置71U及び下対象位置71Lの上下方向Dzの変位量を求めることができる。よって、本実施形態では、変位量を求める際の計算負荷を抑えることができる。このため、本実施形態では、フランジ面の推定の準備期間を短縮し且つその推定コストを抑えることができる。
 ところで、上フランジ面33U中の上対象位置71Uと、下フランジ面33L中の下対象位置71Lとの上下方向Dzの中点を対象接触位置73とすることも可能である。フランジ面の変形には、軸線方向Dyの変化に伴う上下方向Dzの変形のみならず、図9に示すように、横方向Dxの変化に伴う上下方向Dzの変形も含まれる。仮に、下対象位置71Lや上対象位置71Uがフランジ面中の内側縁の位置であるとし、前述したように、上対象位置71U及び下対象位置71Lを用いて対象接触位置73を求めるとする。この場合、フランジ面における横方向Dxの変化に伴う上下方向Dzの変形が、求める対象接触位置73に極端に反映されて、対象接触位置73の上下方向Dzの誤差が大きくなり、結果として、上対象位置71U及び下対象位置71Lの変位量の誤差が大きくなる場合がある。一方、本実施形態では、上フランジ面33U中で横方向Dxにおける中点である上対象中点位置75Uと下フランジ面33L中で横方向Dxにおける中点である下対象中点位置75Lとの上下方向Dzの中点を対象接触位置73にしている。このため、本実施形態では、フランジ面における横方向Dxの変化に伴う上下方向Dzの変形が、求める対象接触位置73に極端に反映されることがなく、対象接触位置73の上下方向Dzの誤差を小さくでき、結果として、上対象位置71U及び下対象位置71Lの変位量の誤差を小さくすることができる。
 以上、本開示の実施形態について詳述したが、本開示は上記実施形態に限定されるものではない。特許請求の範囲に規定された内容及びその均等物から導き出される本発明の概念的な思想と趣旨を逸脱しない範囲において、種々の追加、変更、置き換え、部分的削除等が可能である。
「付記」
 以上の実施形態における回転機械のフランジ変位量推定方法は、例えば、以下のように把握される。
(1)第一態様における回転機械のフランジ変位量推定方法は、以下の回転機械に適用される。
 この回転機械は、水平方向に延びる軸線Arを中心として回転可能なロータ15と、前記ロータ15の外周を覆うケーシング30と、前記ケーシング30内に配置され、前記ケーシング30に取り付けられている静止部品と、前記ケーシング30を下側から支える架台11と、を備える。前記ケーシング30は、上側の上半ケーシング30Uと、下側の下半ケーシング30Lと、前記上半ケーシング30Uと前記下半ケーシング30Lとを締結する複数のボルト39と、を有する。前記上半ケーシング30Uは、下側を向く上フランジ面33Uが形成されている上フランジ32Uを有する。前記下半ケーシング30Lは、上側を向き、前記上フランジ面33Uと上下方向Dzで対向する下フランジ面33Lが形成されている下フランジ32Lと、前記下フランジ32Lに連なり、前記架台11により下側から支えられ、前記軸線Arが延びる軸線方向Dyで互いに離れた第一被支持部35a及び第二被支持部35bと、を有する。前記上フランジ32U及び前記下フランジ32Lには、上下方向Dzに貫通して、前記複数のボルト39のそれぞれが挿通可能なボルト孔34が形成されている。
 以上の回転機械のフランジ変位量推定方法では、
 前記回転機械を分解した後で、前記上半ケーシング30Uと前記下半ケーシング30Lとが前記複数のボルト39で締結されていない開放状態における、前記上フランジ面33U中の複数位置における実測三次元座標データ及び前記下フランジ面33L中の複数位置における実測三次元座標データを受け付ける実測座標受付工程S1と、前記下フランジ面33L中の複数位置における前記実測三次元座標データを用いて、下第一位置72Laと、下第二位置72Lbと、下対象位置71Lと、下対象中点位置75Lと、における有効三次元座標データを把握すると共に、前記上フランジ面33U中の複数位置における前記実測三次元座標データを用いて、上第一位置72Uaと、上第二位置72Ubと、上対象位置71Uと、上対象中点位置75Uと、における有効三次元座標データを把握する有効座標把握工程S2と、前記有効座標把握工程S2で把握した前記下第一位置72Laの有効三次元座標データと前記上第一位置72Uaの有効三次元座標データとが一致し、前記有効座標把握工程S2で把握した前記下第二位置72Lbの有効三次元座標データと前記上第二位置72Ubの有効三次元座標データとが一致するよう、前記有効座標把握工程S2で把握した有効三次元座標データを変更する座標変更工程S3と、前記座標変更工程S3で変更された、前記下対象中点位置75L及び前記上対象中点位置75Uにおける有効三次元座標データを用いて、前記下対象中点位置75Lと前記上対象中点位置75Uとの上下方向Dzにおける中間位置である対象接触位置73の有効三次元座標データを求める接触位置推定工程S4と、前記開放状態から前記上半ケーシング30Uと前記下半ケーシング30Lとが前記複数のボルト39で締結された締結状態になったときの前記上対象位置71U及び前記下対象位置71Lの上下方向Dzの変位量を求める変位量演算工程S5と、を実行する。前記下第一位置72Laは、前記上フランジ面33Uに連なる面中で前記第一被支持部35aの第一代表位置74aと水平方向における位置が一致している位置である。前記下第二位置72Lbは、前記下フランジ面33Lに連なる面中で前記第二被支持部35bの第二代表位置74bと水平方向における位置が一致している位置である。前記下対象位置71Lは、前記下フランジ面33L中で、前記開放状態から前記締結状態になったときの上下方向Dzの変位量を得たい位置である。前記下対象中点位置75Lは、前記下フランジ面33L中で、水平方向で前記軸線方向Dyに垂直な方向である横方向Dxにおける中点の位置であって、前記下対象位置71Lと前記軸線方向Dyの位置が一致している位置である。前記上第一位置72Uaは、前記上フランジ面33Uに連なる面中で前記第一被支持部35aの前記第一代表位置74aと水平方向における位置が一致している位置である。前記上第二位置72Ubは、前記上フランジ面33Uに連なる面中で前記第二被支持部35bの前記第二代表位置74bと水平方向における位置が一致している位置である。前記上対象位置71Uは、前記上フランジ面33U中で前記下対象位置71Lと水平方向の位置が一致している位置である。前記上対象中点位置75Uは、前記上フランジ面33U中で前記横方向Dxにおける中点であって前記下対象位置71Lと前記軸線方向Dyの位置が一致している位置である。前記変位量演算工程S5では、前記座標変更工程S3で変更された前記下対象位置71Lの有効三次元座標データが示す上下方向Dzの位置と前記対象接触位置73の有効三次元座標データが示す上下方向Dzの位置との差を前記下対象位置71Lの上下方向Dzの変位量とし、前記座標変更工程S3で変更された前記上対象位置71Uの有効三次元座標データが示す上下方向Dzの位置と前記対象接触位置73の有効三次元座標データが示す上下方向Dzの位置との差を前記上対象位置71Uの上下方向Dzの変位量とする。
 本態様では、上フランジ面33U中の上対象中点位置75Uと、下フランジ面33L中の下対象中点位置75Lとの上下方向Dzの中点を対象接触位置73とする。そして、本態様では、上フランジ面33U中で上下方向Dzの変位量を得たい上対象位置71Uと対象接触位置73との上下方向Dzの差を上対象位置71Uの変位量とする。さらに、本態様では、下フランジ面33L中で上下方向Dzの変位量を得たい下対象位置71Lと対象接触位置73との上下方向Dzの差を下対象位置71Lの変位量とする。このため、本態様では、下半ケーシング30L及び上半ケーシング30Uの有限要素モデルを用いて、下半ケーシング30L及び上半ケーシング30Uの変形をシミュレートしなくても、上対象位置71U及び下対象位置71Lの上下方向Dzの変位量を求めることができる。よって、本態様では、変位量を求める際の計算負荷を抑えることができる。
 ところで、上フランジ面33U中の上対象位置71Uと、下フランジ面33L中の下対象位置71Lとの上下方向Dzの中点を対象接触位置73とすることも可能である。フランジ面の変形には、軸線方向Dyの変化に伴う上下方向Dzの変形のみならず、横方向Dxの変化に伴う上下方向Dzの変形も含まれる。仮に、下対象位置71Lや上対象位置71Uがフランジ面中の内側縁の位置であるとし、前述したように、上対象位置71U及び下対象位置71Lを用いて対象接触位置73を求めるとする。この場合、フランジ面における横方向Dxの変化に伴う上下方向Dzの変形が、求める対象接触位置73に極端に反映されて、対象接触位置73の上下方向Dzの誤差が大きくなり、結果として、上対象位置71U及び下対象位置71Lの変位量の誤差が大きくなる場合がある。一方、本態様では、上フランジ面33U中で横方向Dxにおける中点である上対象中点位置75Uと下フランジ面33L中で横方向Dxにおける中点である下対象中点位置75Lとの上下方向Dzの中点を対象接触位置73にしている。このため、本態様では、フランジ面における横方向Dxの変化に伴う上下方向Dzの変形が、求める対象接触位置73に極端に反映されることがなく、対象接触位置73の上下方向Dzの誤差を小さくでき、結果として、上対象位置71U及び下対象位置71Lの変位量の誤差を小さくすることができる。
(2)第二態様における回転機械のフランジ変位量推定方法は、
 前記第一態様における回転機械のフランジ変位量推定方法において、前記下対象位置71Lは、前記軸線方向Dyで前記静止部品が配置される位置で且つ前記下フランジ面33L中の内側縁の位置である。
 回転機械の性能等の観点から、静止部品とロータ15との間の径方向Drの間隔を管理する必要がある。発明者は、ケーシング30が開放状態から締結状態になったことによる下半ケーシング30L及び上半ケーシング30Uの変形に伴う静止部品とロータ15との間の径方向Drの間隔の変化は、下フランジ面33L中で軸線方向Dyにおける静止部品格納部36が形成されている位置であって下フランジ面33L中の内側縁の位置の変形、及び、上フランジ面33U中で軸線方向Dyにおける静止部品格納部36が形成されている位置であって上フランジ面33U中の内側縁の位置の変形に対して支配的である、ことを見出した。よって、本態様では、ケーシング30が開放状態から締結状態になったときの静止部品とロータ15との間の径方向Drの間隔を高精度で管理することができる。
(3)第三態様における回転機械のフランジ変位量推定方法は、
 前記第一態様又は前記第二態様における回転機械のフランジ変位量推定方法において、前記実測座標受付工程S1では、前記下対象中点位置75L及び前記上対象中点位置75Uにおける実測三次元座標データを受け付ける。前記有効座標把握工程S2では、前記下対象中点位置75Lにおける実測三次元座標データをそのまま前記下対象中点位置75Lにおける有効三次元座標データとして把握し、前記実測座標受付工程S1で取得した前記上対象中点位置75Uにおける実測三次元座標データをそのまま前記上対象中点位置75Uにおける有効三次元座標データとして把握する。
 本態様では、実測座標受付工程S1で受け付けた下対象中点位置75L及び上対象中点位置75Uにおける実測三次元座標データをそのまま下対象中点位置75L及び上対象中点位置75Uにおける有効三次元座標データとして把握するので、計算負荷を抑えることができる。
(4)第四態様における回転機械のフランジ変位量推定方法は、
 前記第一態様又は前記第二態様における回転機械のフランジ変位量推定方法において、前記実測座標受付工程S1では、前記下対象中点位置75Lを通り且つフランジ幅方向Dwに延びる下中点仮想線上の複数の位置における実測三次元座標データを受け付けると共に、前記上対象中点位置75Uを通り且つフランジ幅方向Dwに延びる上中点仮想線上の複数の位置における実測三次元座標データを受け付ける。前記有効座標把握工程S2では、前記下中点仮想線上の複数の位置における実測三次元座標データから前記下対象中点位置75Lにおける有効三次元座標データを求め、前記上中点仮想線上の複数の位置における実測三次元座標データから前記上対象中点位置75Uにおける有効三次元座標データを求める。
 本態様では、下中点仮想線上の複数の位置における実測三次元座標データから下対象中点位置75Lの有効三次元座標データを求め、上中点仮想線の複数の位置における実測三次元座標データから上対象中点位置75Uの有効三次元座標データを求める。よって、本態様では、下対象中点位置75L及び上対象中点位置75Uにおける有効三次元座標データに関して、局所的な形状変化の影響を受けにくくなる上に、大きな計測誤差が含まれる可能性を抑えることができる。
(5)第五態様における回転機械のフランジ変位量推定方法は、
 前記第一態様又は前記第二態様における回転機械のフランジ変位量推定方法において、前記実測座標受付工程S1では、前記下フランジ面33L中で前記下対象中点位置75Lを含む下中点計測領域中の複数の位置における実測三次元座標データを受け付けると共に、前記上フランジ面33Uで前記上対象中点位置75Uを含む上中点計測領域中の複数の位置における実測三次元座標データを受け付ける。前記有効座標把握工程S2では、前記実測座標受付工程S1で受け付けた前記下中点計測領域中の複数の位置における実測三次元座標データを用いて、前記下対象中点位置75Lにおける有効三次元座標データを求め、前記実測座標受付工程S1で受け付けた前記上中点計測領域中の複数の位置における実測三次元座標データを用いて、前記上対象中点位置75Uにおける有効三次元座標データを求める。
 本態様では、下中点計測領域中の複数の位置における実測三次元座標データから下対象中点位置75Lの有効三次元座標データを求め、上中点計測領域中の複数の位置における実測三次元座標データから上対象中点位置75Uの有効三次元座標データを求める。よって、本態様では、下対象中点位置75L及び上対象中点位置75Uにおける有効三次元座標データに関して、局所的な形状変化の影響を受けにくくなる上に、大きな計測誤差が含まれる可能性を抑えることができる。
(6)第六態様における回転機械のフランジ変位量推定方法は、
 前記第一態様又は前記第二態様における回転機械のフランジ変位量推定方法において、前記実測座標受付工程S1では、前記下フランジ面33L中で前記横方向Dxにおける中点の位置であって、前記軸線方向Dyにおける位置が互いに異なる複数の下中点位置75Lxと、前記上フランジ面33U中で前記横方向Dxにおける中点の位置であって、前記軸線方向Dyにおける位置が互いに異なる複数の上中点位置75Uxと、における実測三次元座標データを受け付ける。前記有効座標把握工程S2では、前記複数の下中点位置75Lxにおける実測三次元座標データの変化傾向から、前記下対象中点位置75Lにおける有効三次元座標データを求め、前記複数の上中点位置75Uxにおける有効三次元座標データの変化傾向から、前記上対象中点位置75Uにおける有効三次元座標データを求める。
 本態様で、複数の下中点位置75Lxにおける実測三次元座標データの変化傾向から、下対象中点位置75Lにおける有効三次元座標データを求め、複数の上中点位置75Uxにおける有効三次元座標データの変化傾向から、前記上対象中点位置75Uにおける有効三次元座標データを求める。よって、本態様では、実測座標受付工程S1で、下対象中点位置75Lや上対象中点位置75Uにおける実測三次元座標データを受け付けていない場合でも、これらの位置における有効三次元座標データを把握することができる。
(7)第七態様における回転機械のフランジ変位量推定方法は、
 前記第一態様から前記第六態様のうちのいずれか一態様における回転機械のフランジ変位量推定方法において、前記実測座標受付工程S1では、前記下第一位置72La及び前記下第二位置72Lbの実測三次元座標データを受け付ける。前記有効座標把握工程S2では、前記実測座標受付工程S1で取得した前記下第一位置72La及び前記下第二位置72Lbにおける実測三次元座標データをそのまま前記下第一位置72La及び前記下第二位置72Lbにおける有効三次元座標データとして把握する。
 本態様では、実測座標受付工程S1で受け付けた下第一位置72La及び下第二位置72Lbにおける実測三次元座標データをそのまま下第一位置72La及び下第二位置72Lbにおける有効三次元座標データとして把握するので、計算負荷を抑えることができる。
(8)第八態様における回転機械のフランジ変位量推定方法は、
 前記第一態様から前記第六態様のうちのいずれか一態様における回転機械のフランジ変位量推定方法において、前記実測座標受付工程S1では、前記第一被支持部35aの表面中の複数の位置における実測三次元座標データ、及び前記第二被支持部35bの表面中の複数の位置における実測三次元座標データを受け付ける。前記第一被支持部35aの前記表面及び前記第二被支持部35bの前記表面は、前記上フランジ面33Uと前記下フランジ面33Lとのうちで前記一方の半ケーシングが有するフランジ面に連なる面である。前記有効座標把握工程S2では、前記実測座標受付工程S1で取得した前記第一被支持部35aの前記表面中の複数の位置における実測三次元座標データから前記下第一位置72Laにおける有効三次元座標データを求め、前記実測座標受付工程S1で取得した前記第二被支持部35bの前記表面中の複数の位置における実測三次元座標データから前記下第二位置72Lbの有効三次元座標データを求める。
 本態様では、実測座標受付工程S1で受け付けた第一被支持部35aにおける上面35ap中の複数の実測三次元座標データから下第一位置72Laの有効三次元座標データを求め、実測座標受付工程S1で受け付けた第二被支持部35bにおける上面35bp中の複数の実測三次元座標データから下第二位置72Lbの有効三次元座標データを求める。よって、本態様では、下第一位置72La及び下第二位置72Lbにおける有効三次元座標データに関して、局所的な形状変化の影響を受けにくくなる上に、大きな計測誤差が含まれる可能性を抑えることができる。
(9)第九態様における回転機械のフランジ変位量推定方法は、
 前記第六態様における回転機械のフランジ変位量推定方法において、前記有効座標把握工程S2では、前記複数の下中点位置75Lxにおける有効三次元座標データの変化傾向から、前記下第一位置72La及び前記下第二位置72Lbにおける有効三次元座標データを求め、前記複数の上中点位置75Uxにおける有効三次元座標データの変化傾向から、前記上第一位置72Ua及び前記上第二位置72Ubにおける有効三次元座標データを求める。
 本態様で、複数の下中点位置75Lxにおける有効三次元座標データの変化傾向から、下第一位置72La及び下第二位置72Lbにおける有効三次元座標を求め、複数の上中点位置75Uxにおける有効三次元座標データの変化傾向から、上第一位置72Ua及び上第二位置72Ubにおける有効三次元座標データを求める。よって、本態様では、実測座標受付工程S1で、下第一位置72La、下第二位置72Lb、上第一位置72Ua、及び上第二位置72Ubにおける実測三次元座標データを受け付けていない場合でも、これらの位置における有効三次元座標データを把握することができる。
(10)第十態様における回転機械のフランジ変位量推定方法は、
 前記第一態様から前記第九態様のうちのいずれか一態様における回転機械のフランジ変位量推定方法において、前記実測座標受付工程S1では、前記下対象位置71L及び前記上対象位置71Uにおける実測三次元座標データを受け付ける。前記有効座標把握工程S2では、前記下対象位置71Lにおける実測三次元座標データをそのまま前記下対象位置71Lにおける有効三次元座標データとして把握し、前記実測座標受付工程S1で取得した前記上対象位置71Uにおける実測三次元座標データをそのまま前記上対象位置71Uにおける有効三次元座標データとして把握する。
 本態様では、実測座標受付工程S1で受け付けた下対象位置71L及び上対象位置71Uにおける実測三次元座標データをそのまま下対象位置71L及び上対象位置71Uにおける有効三次元座標データとして把握するので、計算負荷を抑えることができる。
(11)第十一態様における回転機械のフランジ変位量推定方法は、
 前記第一態様から前記第九態様のうちのいずれか一態様における回転機械のフランジ変位量推定方法において、前記実測座標受付工程S1では、前記下対象位置71Lを通り且つフランジ幅方向Dwに延びる下仮想線上の複数の位置における実測三次元座標データを受け付けると共に、前記上対象位置71Uを通り且つフランジ幅方向Dwに延びる上仮想線上の複数の位置における実測三次元座標データを受け付ける。前記有効座標把握工程S2では、前記下仮想線上の複数の位置における実測三次元座標データから前記下対象位置71Lにおける有効三次元座標データを求め、前記上仮想線上の複数の位置における実測三次元座標データから前記上対象位置71Uにおける有効三次元座標データを求める。
 本態様では、下仮想線上の複数の位置における実測三次元座標データから下対象位置71Lにおける有効三次元座標データを求め、上仮想線上の複数の位置における実測三次元座標データから上対象位置71Uにおける有効三次元座標データを求める。よって、本態様では、下対象位置71L及び上対象位置71Uにおける有効三次元座標データに関して、局所的な形状変化の影響を受けにくくなる上に、大きな計測誤差が含まれる可能性を抑えることができる。
(12)第十二態様における回転機械のフランジ変位量推定方法は、
 前記第一態様から前記第九態様のうちのいずれか一態様における回転機械のフランジ変位量推定方法において、前記実測座標受付工程S1では、前記下フランジ面33L中で前記下対象位置71Lを含む下計測領域中の複数の位置における実測三次元座標データを受け付けると共に、前記上フランジ面33Uで前記上対象位置71Uを含む上計測領域中の複数の位置における実測三次元座標データを受け付ける。前記有効座標把握工程S2では、前記実測座標受付工程S1で受け付けた前記下計測領域中の複数の位置における実測三次元座標データを用いて、前記下対象位置71Lにおける有効三次元座標データを求め、前記実測座標受付工程S1で受け付けた前記上計測領域中の複数の位置における実測三次元座標データを用いて、前記上対象位置71Uにおける有効三次元座標データを求める。
 本態様では、下計測領域中の複数の位置における実測三次元座標データから下対象位置71Lの有効三次元座標データを求め、上計測領域中の複数の位置における実測三次元座標データから上対象位置71Uの有効三次元座標データを求める。よって、本態様では、下対象位置71L及び上対象位置71Uにおける有効三次元座標データに関して、局所的な形状変化の影響を受けにくくなる上に、大きな計測誤差が含まれる可能性を抑えることができる。
(13)第十三態様における回転機械のフランジ変位量推定方法は、
 前記第一態様又は前記第二態様における回転機械のフランジ変位量推定方法において、前記実測座標受付工程S1では、前記下フランジ面33Lの全体渡る複数の位置における実測三次元座標データを受け付けると共に、前記上フランジ面33Uの全体に渡る複数の位置における実測三次元座標データを受け付ける。前記有効座標把握工程S2では、前記実測座標受付工程S1で受け付けた前記下フランジ面33Lの全体渡る複数の位置における実測三次元座標データを用いて、前記下フランジ面33L全体の三次元形状を示す下フランジ面33Lの形状データを求めると共に、前記実測座標受付工程S1で受け付けた前記上フランジ面33Uの全体渡る複数の位置における実測三次元座標データを用いて、前記上フランジ面33U全体の三次元形状を示す上フランジ面33Uの形状データを求める。前記有効座標把握工程S2では、さらに、前記下フランジ面33Lの形状データを用いて、前記下対象中点位置75Lの有効三次元座標データを求めると共に、前記上フランジ面33Uの形状データを用いて、前記上対象中点位置75Uの有効三次元座標データを求める。
 本態様では、下フランジ面33Lの全体に渡る複数の位置における実測三次元座標データから下対象中点位置75Lの有効三次元座標データを求め、上フランジ面33Uの全体に渡る複数の位置における実測三次元座標データから上対象中点位置75Uの有効三次元座標データを求める。よって、本態様では、下対象中点位置75L及び上対象中点位置75Uにおける有効三次元座標データに関して、局所的な形状変化の影響を受けにくくなる上に、大きな計測誤差が含まれる可能性を抑えることができる。さらに、本態様では、さらに、障害物等に起因した広範囲なデータ欠損がある場合でも、これらの位置における有効三次元座標データを把握できる。
(14)第十四態様における回転機械のフランジ変位量推定方法は、
 前記第十三態様における回転機械のフランジ変位量推定方法において、前記有効座標把握工程S2では、前記下フランジ面33Lの形状データを用いて、前記下第一位置72La及び前記下第二位置72Lbにおける有効三次元座標データを求めると共に、前記上フランジ面33Uの形状データを用いて、前記上第一位置72Ua及び前記上第二位置72Ubにおける有効三次元座標データを求める。
 本態様では、下フランジ面33Lの全体に渡る複数の位置における実測三次元座標データから下第一位置72La及び下第二位置72Lbの有効三次元座標データを求め、上フランジ面33Uの全体に渡る複数の位置における実測三次元座標データから上第一位置72Ua及び上第二位置72Ubの有効三次元座標データを求める。よって、本態様では、下第一位置72La、下第二位置72Lb、上第一位置72Ua、及び上第二位置72Ubにおける実測三次元座標データを受け付けていない場合でも、これらの位置における有効三次元座標データを把握することができる。また、本態様では、下第一位置72La、下第二位置72Lb、上第一位置72Ua、及び上第二位置72Ubにおける有効三次元座標データに関して、局所的な形状変化の影響を受けにくくなる上に、大きな計測誤差が含まれる可能性を抑えることができる。さらに、本態様では、さらに、障害物等に起因した広範囲なデータ欠損がある場合でも、これらの位置における有効三次元座標データを把握できる。
(15)第十五態様における回転機械のフランジ変位量推定方法は、
 前記十三態様又は前記第十四態様における回転機械のフランジ変位量推定方法において、前記有効座標把握工程S2では、前記下フランジ面33Lの形状データを用いて、前記下対象位置71Lにおける有効三次元座標データを求めると共に、前記上フランジ面33Uの形状データを用いて、前記上対象位置71Uにおける有効三次元座標データを求める。
 本態様では、下フランジ面33Lの全体に渡る複数の位置における実測三次元座標データから下対象位置71Lの有効三次元座標データを求め、上フランジ面33Uの全体に渡る複数の位置における実測三次元座標データから上対象位置71Uの有効三次元座標データを求める。よって、本態様では、下対象位置71L及び上対象位置71Uにおける有効三次元座標データに関して、局所的な形状変化の影響を受けにくくなる上に、大きな計測誤差が含まれる可能性を抑えることができる。さらに、本態様では、さらに、障害物等に起因した広範囲なデータ欠損がある場合でも、これらの位置における有効三次元座標データを把握できる。
 以上の実施形態における回転機械のフランジ変位量推定プログラムは、例えば、以下のように把握される。
(16)第十六態様における回転機械のフランジ変位量推定プログラムは、以下の回転機械に適用される。
 この回転機械は、水平方向に延びる軸線Arを中心として回転可能なロータ15と、前記ロータ15の外周を覆うケーシング30と、前記ケーシング30内に配置され、前記ケーシング30に取り付けられている静止部品と、前記ケーシング30を下側から支える架台11と、を備える。前記ケーシング30は、上側の上半ケーシング30Uと、下側の下半ケーシング30Lと、前記上半ケーシング30Uと前記下半ケーシング30Lとを締結する複数のボルト39と、を有する。前記上半ケーシング30Uは、下側を向く上フランジ面33Uが形成されている上フランジ32Uを有する。前記下半ケーシング30Lは、上側を向き、前記上フランジ面33Uと上下方向Dzで対向する下フランジ面33Lが形成されている下フランジ32Lと、前記下フランジ32Lに連なり、前記架台11により下側から支えられ、前記軸線Arが延びる軸線方向Dyで互いに離れた第一被支持部35a及び第二被支持部35bと、を有する。前記上フランジ32U及び前記下フランジ32Lには、上下方向Dzに貫通して、前記複数のボルト39のそれぞれが挿通可能なボルト孔34が形成されている。
 以上の回転機械のフランジ変位量推定プログラム58pは、前記回転機械を分解した後で、前記上半ケーシング30Uと前記下半ケーシング30Lとが前記複数のボルト39で締結されていない開放状態における、前記上フランジ面33U中の複数位置における実測三次元座標データ及び前記下フランジ面33L中の複数位置における実測三次元座標データを受け付ける実測座標受付工程S1と、前記下フランジ面33L中の複数位置における前記実測三次元座標データを用いて、下第一位置72Laと、下第二位置72Lbと、下対象位置71Lと、下対象中点位置75Lと、における有効三次元座標データを把握すると共に、前記上フランジ面33U中の複数位置における前記実測三次元座標データを用いて、上第一位置72Uaと、上第二位置72Ubと、上対象位置71Uと、上対象中点位置75Uと、における有効三次元座標データを把握する有効座標把握工程S2と、前記有効座標把握工程S2で把握した前記下第一位置72Laの有効三次元座標データと前記上第一位置72Uaの有効三次元座標データとが一致し、前記有効座標把握工程S2で把握した前記下第二位置72Lbの有効三次元座標データと前記上第二位置72Ubの有効三次元座標データとが一致するよう、前記有効座標把握工程S2で把握した有効三次元座標データを変更する座標変更工程S3と、前記座標変更工程S3で変更された、前記下対象中点位置75L及び前記上対象中点位置75Uにおける有効三次元座標データを用いて、前記下対象中点位置75Lと前記上対象中点位置75Uとの上下方向Dzにおける中間位置である対象接触位置73の有効三次元座標データを求める接触位置推定工程S4と、前記開放状態から前記上半ケーシング30Uと前記下半ケーシング30Lとが前記複数のボルト39で締結された締結状態になったときの前記上対象位置71U及び前記下対象位置71Lの上下方向Dzの変位量を求める変位量演算工程S5と、をコンピュータに実行させる。前記下第一位置72Laは、前記上フランジ面33Uに連なる面中で前記第一被支持部35aの第一代表位置74aと水平方向における位置が一致している位置である。前記下第二位置72Lbは、前記下フランジ面33Lに連なる面中で前記第二被支持部35bの第二代表位置74bと水平方向における位置が一致している位置である。前記下対象位置71Lは、前記下フランジ面33L中で、前記開放状態から前記締結状態になったときの上下方向Dzの変位量を得たい位置である。前記下対象中点位置75Lは、前記下フランジ面33L中で、水平方向で前記軸線方向Dyに垂直な方向である横方向Dxにおける中点の位置であって、前記下対象位置71Lと前記軸線方向Dyの位置が一致している位置である。前記上第一位置72Uaは、前記上フランジ面33Uに連なる面中で前記第一被支持部35aの前記第一代表位置74aと水平方向における位置が一致している位置である。前記上第二位置72Ubは、前記上フランジ面33Uに連なる面中で前記第二被支持部35bの前記第二代表位置74bと水平方向における位置が一致している位置である。前記上対象位置71Uは、前記上フランジ面33U中で前記下対象位置71Lと水平方向の位置が一致している位置である。前記上対象中点位置75Uは、前記上フランジ面33U中で前記横方向Dxにおける中点であって前記下対象位置71Lと前記軸線方向Dyの位置が一致している位置である。前記変位量演算工程S5では、前記座標変更工程S3で変更された前記下対象位置71Lの有効三次元座標データが示す上下方向Dzの位置と前記対象接触位置73の有効三次元座標データが示す上下方向Dzの位置との差を前記下対象位置71Lの上下方向Dzの変位量とし、前記座標変更工程S3で変更された前記上対象位置71Uの有効三次元座標データが示す上下方向Dzの位置と前記対象接触位置73の有効三次元座標データが示す上下方向Dzの位置との差を前記上対象位置71Uの上下方向Dzの変位量とする。
 本態様では、このプログラムをコンピュータに実行させることにより、第一態様と同様に、変位量を求める際の計算負荷を抑えることができる。
(17)第十七態様における回転機械のフランジ変位量推定プログラムは、
 前記第十六態様における回転機械のフランジ変位量推定プログラム58pにおいて、前記下対象位置71Lは、前記軸線方向Dyで前記静止部品が配置される位置で且つ前記下フランジ面33L中の内側縁の位置である。
(18)第十八態様における回転機械のフランジ変位量推定プログラムは、
 前記第十六態様又は前記第十七態様における回転機械のフランジ変位量推定プログラム58pにおいて、前記実測座標受付工程S1では、前記下フランジ面33L中で前記横方向Dxにおける中点の位置であって、前記軸線方向Dyにおける位置が互いに異なる複数の下中点位置75Lxと、前記上フランジ面33U中で前記横方向Dxにおける中点の位置であって、前記軸線方向Dyにおける位置が互いに異なる複数の上中点位置75Uxと、における実測三次元座標データを受け付ける。前記有効座標把握工程S2では、前記複数の下中点位置75Lxにおける実測三次元座標データの変化傾向から、前記下対象中点位置75Lにおける有効三次元座標データを求め、前記複数の上中点位置75Uxにおける有効三次元座標データの変化傾向から、前記上対象中点位置75Uにおける有効三次元座標データを求める。
 本態様で、このプログラムをコンピュータに実行させることにより、第六態様と同様に、実測座標受付工程S1で、下対象中点位置75Lや上対象中点位置75Uにおける実測三次元座標データを受け付けていない場合でも、これらの位置における有効三次元座標データを把握することができる。
(19)第十九態様における回転機械のフランジ変位量推定プログラムは、
 前記第十八態様における回転機械のフランジ変位量推定プログラム58pにおいて、前記有効座標把握工程S2では、前記複数の下中点位置75Lxにおける有効三次元座標データの変化傾向から、前記下第一位置72La及び前記下第二位置72Lbにおける有効三次元座標データを求め、前記複数の上中点位置75Uxにおける有効三次元座標データの変化傾向から、前記上第一位置72Ua及び前記上第二位置72Ubにおける有効三次元座標データを求める。
 本態様で、このプログラムをコンピュータに実行させることにより、第九態様と同様に、実測座標受付工程S1で、下第一位置72La、下第二位置72Lb、上第一位置72Ua、及び上第二位置72Ubにおける実測三次元座標データを受け付けていない場合でも、これらの位置における有効三次元座標データを把握することができる。
 以上の実施形態における回転機械のフランジ変位量推定装置は、例えば、以下のように把握される。
(20)第二十態様における回転機械のフランジ変位量推定装置は、以下の回転機械に適用される。
 この回転機械は、水平方向に延びる軸線Arを中心として回転可能なロータ15と、前記ロータ15の外周を覆うケーシング30と、前記ケーシング30内に配置され、前記ケーシング30に取り付けられている静止部品と、前記ケーシング30を下側から支える架台11と、を備える。前記ケーシング30は、上側の上半ケーシング30Uと、下側の下半ケーシング30Lと、前記上半ケーシング30Uと前記下半ケーシング30Lとを締結する複数のボルト39と、を有する。前記上半ケーシング30Uは、下側を向く上フランジ面33Uが形成されている上フランジ32Uを有する。前記下半ケーシング30Lは、上側を向き、前記上フランジ面33Uと上下方向Dzで対向する下フランジ面33Lが形成されている下フランジ32Lと、前記下フランジ32Lに連なり、前記架台11により下側から支えられ、前記軸線Arが延びる軸線方向Dyで互いに離れた第一被支持部35a及び第二被支持部35bと、を有する。前記上フランジ32U及び前記下フランジ32Lには、上下方向Dzに貫通して、前記複数のボルト39のそれぞれが挿通可能なボルト孔34が形成されている。
 以上の回転機械のフランジ変位量推定装置50において、前記回転機械を分解した後で、前記上半ケーシング30Uと前記下半ケーシング30Lとが前記複数のボルト39で締結されていない開放状態における、前記上フランジ面33U中の複数位置における実測三次元座標データ及び前記下フランジ面33L中の複数位置における実測三次元座標データを受け付ける実測座標受付部61と、前記下フランジ面33L中の複数位置における前記実測三次元座標データを用いて、下第一位置72Laと、下第二位置72Lbと、下対象位置71Lと、下対象中点位置75Lと、における有効三次元座標データを把握すると共に、前記上フランジ面33U中の複数位置における前記実測三次元座標データを用いて、上第一位置72Uaと、上第二位置72Ubと、上対象位置71Uと、上対象中点位置75Uと、における有効三次元座標データを把握する有効座標把握部62と、前記有効座標把握部62が把握した前記下第一位置72Laの有効三次元座標データと前記上第一位置72Uaの有効三次元座標データとが一致し、前記有効座標把握部62が把握した前記下第二位置72Lbの有効三次元座標データと前記上第二位置72Ubの有効三次元座標データとが一致するよう、前記有効座標把握部62が把握した有効三次元座標データを変更する座標変更部63と、前記座標変更部63が変更した、前記下対象中点位置75L及び前記上対象中点位置75Uにおける有効三次元座標データを用いて、前記下対象中点位置75Lと前記上対象中点位置75Uとの上下方向Dzにおける中間位置である対象接触位置73の有効三次元座標データを求める接触位置推定部64と、前記開放状態から前記上半ケーシング30Uと前記下半ケーシング30Lとが前記複数のボルト39で締結された締結状態になったときの前記上対象位置71U及び前記下対象位置71Lの上下方向Dzの変位量を求める変位量演算部65と、を備える。前記下第一位置72Laは、前記上フランジ面33Uに連なる面中で前記第一被支持部35aの第一代表位置74aと水平方向における位置が一致している位置である。前記下第二位置72Lbは、前記下フランジ面33Lに連なる面中で前記第二被支持部35bの第二代表位置74bと水平方向における位置が一致している位置である。前記下対象位置71Lは、前記下フランジ面33L中で、前記開放状態から前記締結状態になったときの上下方向Dzの変位量を得たい位置である。前記下対象中点位置75Lは、前記下フランジ面33L中で、水平方向で前記軸線方向Dyに垂直な方向である横方向Dxにおける中点の位置であって、前記下対象位置71Lと前記軸線方向Dyの位置が一致している位置である。前記上第一位置72Uaは、前記上フランジ面33Uに連なる面中で前記第一被支持部35aの前記第一代表位置74aと水平方向における位置が一致している位置である。前記上第二位置72Ubは、前記上フランジ面33Uに連なる面中で前記第二被支持部35bの前記第二代表位置74bと水平方向における位置が一致している位置である。前記上対象位置71Uは、前記上フランジ面33U中で前記下対象位置71Lと水平方向の位置が一致している位置である。前記上対象中点位置75Uは、前記上フランジ面33U中で前記横方向Dxにおける中点であって前記下対象位置71Lと前記軸線方向Dyの位置が一致している位置である。前記変位量演算部65は、前記座標変更部63が変更した前記下対象位置71Lの有効三次元座標データが示す上下方向Dzの位置と前記対象接触位置73の有効三次元座標データが示す上下方向Dzの位置との差を前記下対象位置71Lの上下方向Dzの変位量とし、前記座標変更部63が変更した前記上対象位置71Uの有効三次元座標データが示す上下方向Dzの位置と前記対象接触位置73の有効三次元座標データが示す上下方向Dzの位置との差を前記上対象位置71Uの上下方向Dzの変位量とする。
 本態様では、第一態様と同様に、変位量を求める際の計算負荷を抑えることができる。
(21)第二十一態様における回転機械のフランジ変位量推定装置は、
 前記第二十態様における回転機械のフランジ変位量推定装置50において、前記下対象位置71Lは、前記軸線方向Dyで前記静止部品が配置される位置で且つ前記下フランジ面33L中の内側縁の位置である。
(22)第二十二態様における回転機械のフランジ変位量推定装置は、
 前記第二十態様又は前記第二十一態様における回転機械のフランジ変位量推定装置50において、前記実測座標受付部61は、前記下フランジ面33L中で前記横方向Dxにおける中点の位置であって、前記軸線方向Dyにおける位置が互いに異なる複数の下中点位置75Lxと、前記上フランジ面33U中で前記横方向Dxにおける中点の位置であって、前記軸線方向Dyにおける位置が互いに異なる複数の上中点位置75Uxと、における実測三次元座標データを受け付ける。前記有効座標把握部62は、前記複数の下中点位置75Lxにおける実測三次元座標データの変化傾向から、前記下対象中点位置75Lにおける有効三次元座標データを求め、前記複数の上中点位置75Uxにおける有効三次元座標データの変化傾向から、前記上対象中点位置75Uにおける有効三次元座標データを求める。
 本態様で、第六態様と同様に、実測座標受付部61が、下対象中点位置75Lや上対象中点位置75Uにおける実測三次元座標データを受け付けていない場合でも、これらの位置における有効三次元座標データを把握することができる。
(23)第二十三態様における回転機械のフランジ変位量推定装置は、
 前記第二十二態様における回転機械のフランジ変位量推定装置50において、前記有効座標把握部62は、前記複数の下中点位置75Lxにおける有効三次元座標データの変化傾向から、前記下第一位置72La及び前記下第二位置72Lbにおける有効三次元座標データを求め、前記複数の上中点位置75Uxにおける有効三次元座標データの変化傾向から、前記上第一位置72Ua及び前記上第二位置72Ubにおける有効三次元座標データを求める。
 本態様で、第九態様と同様に、実測座標受付部61が、下第一位置72La、下第二位置72Lb、上第一位置72Ua、及び上第二位置72Ubにおける実測三次元座標データを受け付けていない場合でも、これらの位置における有効三次元座標データを把握することができる。
 本開示の一態様では、上半ケーシング及び下半ケーシングのフランジ面の変位量を推定するにあたり、計算負荷を抑えることで、フランジ面の推定の準備期間を短縮し且つその推定コストを抑えることができる。
10:蒸気タービン(回転機械)
11:架台
12a:第一軸受装置
12b:第二軸受装置
13a:第一軸封装置(静止部品)
13b:第二軸封装置(静止部品)
15:ロータ
16:ロータ軸
17:動翼列
20:ダイヤフラム(静止部品)
20L:下半ダイヤフラム
20U:上半ダイヤフラム
22:静翼
23:ダイヤフラム内輪
24:ダイヤフラム外輪
25:シール装置
30:ケーシング
30L:下半ケーシング
30U:上半ケーシング
31L:下半ケーシング本体
31U:上半ケーシング本体
32L:下フランジ
32U:上フランジ
33L:下フランジ面
33U:上フランジ面
34:ボルト孔
35a:第一被支持部
35ap:上面
35b:第二被支持部
35bp:上面
36:静止部品格納部
39:ボルト
50:フランジ変位量推定装置
51:手入力装置
52:表示装置
53:入出力インタフェース
54:装置インタフェース
55:通信インタフェース
56:記憶・再生装置
57:メモリ
58:補助記憶装置
58d:基準三次元形状データ
58p:フランジ変位量推定プログラム
60:CPU
61:実測座標受付部
62:有効座標把握部
63:座標変更部
64:接触位置推定部
65:変位量演算部
69:三次元形状測定装置
71:参照位置
71L:下対象位置
71U:上対象位置
72La:下第一位置
72Ua:上第一位置
72Lb:下第二位置
72Ub:上第二位置
73:対象接触位置
74a:第一代表位置
74b:第二代表位置
75L:下対象中点位置
75U:上対象中点位置
75Lx:下中点位置
75Ux:上中点位置
76:仮想線
76L:下仮想線
76U:上仮想線
77La:下第一仮想線
77Ua:上第一仮想線
77Lb:下第二仮想線
77Ub:上第二仮想線
79:参照計測領域
80:基準三次元形状データが示すフランジ面
81:基準三次元形状データが示す参照位置
82:基準三次元形状データが示す、フランジ面に対して傾斜している面
83:三次元ブロック
85:点
86,86a,86b:ポリゴン(多角形平面)
87:代表点
Ar:軸線
Dc:周方向
Dr:径方向
Dri:径方向内側
Dro:径方向外側
Dx:横方向
Dy:軸線方向
Dz:上下方向
Dw:フランジ幅方向

Claims (23)

  1.  水平方向に延びる軸線を中心として回転可能なロータと、
     前記ロータの外周を覆うケーシングと、
     前記ケーシング内に配置され、前記ケーシングに取り付けられている静止部品と、
     前記ケーシングを下側から支える架台と、
     を備え、
     前記ケーシングは、上側の上半ケーシングと、下側の下半ケーシングと、前記上半ケーシングと前記下半ケーシングとを締結する複数のボルトと、を有し、
     前記上半ケーシングは、下側を向く上フランジ面が形成されている上フランジを有し、
     前記下半ケーシングは、上側を向き、前記上フランジ面と上下方向で対向する下フランジ面が形成されている下フランジと、前記下フランジに連なり、前記架台により下側から支えられ、前記軸線が延びる軸線方向で互いに離れた第一被支持部及び第二被支持部と、を有し、
     前記上フランジ及び前記下フランジには、上下方向に貫通して、前記複数のボルトのそれぞれが挿通可能なボルト孔が形成されている、
     回転機械のフランジ変位量推定方法において、
     前記回転機械を分解した後で、前記上半ケーシングと前記下半ケーシングとが前記複数のボルトで締結されていない開放状態における、前記上フランジ面中の複数位置における実測三次元座標データ及び前記下フランジ面中の複数位置における実測三次元座標データを受け付ける実測座標受付工程と、
     前記下フランジ面中の複数位置における前記実測三次元座標データを用いて、下第一位置と、下第二位置と、下対象位置と、下対象中点位置と、における有効三次元座標データを把握すると共に、前記上フランジ面中の複数位置における前記実測三次元座標データを用いて、上第一位置と、上第二位置と、上対象位置と、上対象中点位置と、における有効三次元座標データを把握する有効座標把握工程と、
     前記有効座標把握工程で把握した前記下第一位置の有効三次元座標データと前記上第一位置の有効三次元座標データとが一致し、前記有効座標把握工程で把握した前記下第二位置の有効三次元座標データと前記上第二位置の有効三次元座標データとが一致するよう、前記有効座標把握工程で把握した有効三次元座標データを変更する座標変更工程と、
     前記座標変更工程で変更された、前記下対象中点位置及び前記上対象中点位置における有効三次元座標データを用いて、前記下対象中点位置と前記上対象中点位置との上下方向における中間位置である対象接触位置の有効三次元座標データを求める接触位置推定工程と、
     前記開放状態から前記上半ケーシングと前記下半ケーシングとが前記複数のボルトで締結された締結状態になったときの前記上対象位置及び前記下対象位置の上下方向の変位量を求める変位量演算工程と、
     を実行し、
     前記下第一位置は、前記上フランジ面に連なる面中で前記第一被支持部の第一代表位置と水平方向における位置が一致している位置であり、
     前記下第二位置は、前記下フランジ面に連なる面中で前記第二被支持部の第二代表位置と水平方向における位置が一致している位置であり、
     前記下対象位置は、前記下フランジ面中で、前記開放状態から前記締結状態になったときの上下方向の変位量を得たい位置であり、
     前記下対象中点位置は、前記下フランジ面中で、水平方向で前記軸線方向に垂直な方向である横方向における中点の位置であって、前記下対象位置と前記軸線方向の位置が一致している位置であり、
     前記上第一位置は、前記上フランジ面に連なる面中で前記第一被支持部の前記第一代表位置と水平方向における位置が一致している位置であり、
     前記上第二位置は、前記上フランジ面に連なる面中で前記第二被支持部の前記第二代表位置と水平方向における位置が一致している位置であり、
     前記上対象位置は、前記上フランジ面中で前記下対象位置と水平方向の位置が一致している位置であり、
     前記上対象中点位置は、前記上フランジ面中で前記横方向における中点であって前記下対象位置と前記軸線方向の位置が一致している位置であり、
     前記変位量演算工程では、
     前記座標変更工程で変更された前記下対象位置の有効三次元座標データが示す上下方向の位置と前記対象接触位置の有効三次元座標データが示す上下方向の位置との差を前記下対象位置の上下方向の変位量とし、前記座標変更工程で変更された前記上対象位置の有効三次元座標データが示す上下方向の位置と前記対象接触位置の有効三次元座標データが示す上下方向の位置との差を前記上対象位置の上下方向の変位量とする、
     回転機械のフランジ変位量推定方法。
  2.  請求項1に記載の回転機械のフランジ変位量推定方法において、
     前記下対象位置は、前記軸線方向で前記静止部品が配置される位置で且つ前記下フランジ面中の内側縁の位置である、
     回転機械のフランジ変位量推定方法。
  3.  請求項1に記載の回転機械のフランジ変位量推定方法において、
     前記実測座標受付工程では、前記下対象中点位置及び前記上対象中点位置における実測三次元座標データを受け付け、
     前記有効座標把握工程では、前記下対象中点位置における実測三次元座標データをそのまま前記下対象中点位置における有効三次元座標データとして把握し、前記実測座標受付工程で取得した前記上対象中点位置における実測三次元座標データをそのまま前記上対象中点位置における有効三次元座標データとして把握する、
     回転機械のフランジ変位量推定方法。
  4.  請求項1に記載の回転機械のフランジ変位量推定方法において、
     前記実測座標受付工程では、前記下対象中点位置を通り且つフランジ幅方向に延びる下中点仮想線上の複数の位置における実測三次元座標データを受け付けると共に、前記上対象中点位置を通り且つフランジ幅方向に延びる上中点仮想線上の複数の位置における実測三次元座標データを受け付け、
     前記有効座標把握工程では、前記下中点仮想線上の複数の位置における実測三次元座標データから前記下対象中点位置における有効三次元座標データを求め、前記上中点仮想線上の複数の位置における実測三次元座標データから前記上対象中点位置における有効三次元座標データを求める、
     回転機械のフランジ変位量推定方法。
  5.  請求項1に記載の回転機械のフランジ変位量推定方法において、
     前記実測座標受付工程では、前記下フランジ面中で前記下対象中点位置を含む下中点計測領域中の複数の位置における実測三次元座標データを受け付けると共に、前記上フランジ面で前記上対象中点位置を含む上中点計測領域中の複数の位置における実測三次元座標データを受け付け、
     前記有効座標把握工程では、前記実測座標受付工程で受け付けた前記下中点計測領域中の複数の位置における実測三次元座標データを用いて、前記下対象中点位置における有効三次元座標データを求め、前記実測座標受付工程で受け付けた前記上中点計測領域中の複数の位置における実測三次元座標データを用いて、前記上対象中点位置における有効三次元座標データを求める、
     回転機械のフランジ変位量推定方法。
  6.  請求項1に記載の回転機械のフランジ変位量推定方法において、
     前記実測座標受付工程では、前記下フランジ面中で前記横方向における中点の位置であって、前記軸線方向における位置が互いに異なる複数の下中点位置と、前記上フランジ面中で前記横方向における中点の位置であって、前記軸線方向における位置が互いに異なる複数の上中点位置と、における実測三次元座標データを受け付け、
     前記有効座標把握工程では、前記複数の下中点位置における実測三次元座標データの変化傾向から、前記下対象中点位置における有効三次元座標データを求め、前記複数の上中点位置における有効三次元座標データの変化傾向から、前記上対象中点位置における有効三次元座標データを求める、
     回転機械のフランジ変位量推定方法。
  7.  請求項1から6のいずれか一項に記載の回転機械のフランジ変位量推定方法において、
     前記実測座標受付工程では、前記下第一位置及び前記下第二位置の実測三次元座標データを受け付け、
     前記有効座標把握工程では、前記実測座標受付工程で取得した前記下第一位置及び前記下第二位置における実測三次元座標データをそのまま前記下第一位置及び前記下第二位置における有効三次元座標データとして把握する、
     回転機械のフランジ変位量推定方法。
  8.  請求項1から6のいずれか一項に記載の回転機械のフランジ変位量推定方法において、
     前記実測座標受付工程では、前記第一被支持部の上面中の複数の位置における実測三次元座標データ、及び前記第二被支持部の上面中の複数の位置における実測三次元座標データを受け付け、
     前記有効座標把握工程では、前記実測座標受付工程で取得した前記第一被支持部の上面中の複数の位置における実測三次元座標データから前記下第一位置における有効三次元座標データを求め、前記実測座標受付工程で取得した前記第二被支持部の上面中の複数の位置における実測三次元座標データから前記下第二位置の有効三次元座標データを求める、
     回転機械のフランジ変位量推定方法。
  9.  請求項6に記載の回転機械のフランジ変位量推定方法において、
     前記有効座標把握工程では、前記複数の下中点位置における有効三次元座標データの変化傾向から、前記下第一位置及び前記下第二位置における有効三次元座標データを求め、前記複数の上中点位置における有効三次元座標データの変化傾向から、前記上第一位置及び前記上第二位置における有効三次元座標データを求める、
     回転機械のフランジ変位量推定方法。
  10.  請求項1から6のいずれか一項に記載の回転機械のフランジ変位量推定方法において、
     前記実測座標受付工程では、前記下対象位置及び前記上対象位置における実測三次元座標データを受け付け、
     前記有効座標把握工程では、前記下対象位置における実測三次元座標データをそのまま前記下対象位置における有効三次元座標データとして把握し、前記実測座標受付工程で取得した前記上対象位置における実測三次元座標データをそのまま前記上対象位置における有効三次元座標データとして把握する、
     回転機械のフランジ変位量推定方法。
  11.  請求項1から6のいずれか一項に記載の回転機械のフランジ変位量推定方法において、
     前記実測座標受付工程では、前記下対象位置を通り且つフランジ幅方向に延びる下仮想線上の複数の位置における実測三次元座標データを受け付けると共に、前記上対象位置を通り且つフランジ幅方向に延びる上仮想線上の複数の位置における実測三次元座標データを受け付け、
     前記有効座標把握工程では、前記下仮想線上の複数の位置における実測三次元座標データから前記下対象位置における有効三次元座標データを求め、前記上仮想線上の複数の位置における実測三次元座標データから前記上対象位置における有効三次元座標データを求める、
     回転機械のフランジ変位量推定方法。
  12.  請求項1から6のいずれか一項に記載の回転機械のフランジ変位量推定方法において、
     前記実測座標受付工程では、前記下フランジ面中で前記下対象位置を含む下計測領域中の複数の位置における実測三次元座標データを受け付けると共に、前記上フランジ面で前記上対象位置を含む上計測領域中の複数の位置における実測三次元座標データを受け付け、
     前記有効座標把握工程では、前記実測座標受付工程で受け付けた前記下計測領域中の複数の位置における実測三次元座標データを用いて、前記下対象位置における有効三次元座標データを求め、前記実測座標受付工程で受け付けた前記上計測領域中の複数の位置における実測三次元座標データを用いて、前記上対象位置における有効三次元座標データを求める、
     回転機械のフランジ変位量推定方法。
  13.  請求項1又は2に記載の回転機械のフランジ変位量推定方法において、
     前記実測座標受付工程では、前記下フランジ面の全体渡る複数の位置における実測三次元座標データを受け付けると共に、前記上フランジ面の全体に渡る複数の位置における実測三次元座標データを受け付け、
     前記有効座標把握工程では、
     前記実測座標受付工程で受け付けた前記下フランジ面の全体渡る複数の位置における実測三次元座標データを用いて、前記下フランジ面全体の三次元形状を示す下フランジ面の形状データを求めると共に、前記実測座標受付工程で受け付けた前記上フランジ面の全体渡る複数の位置における実測三次元座標データを用いて、前記上フランジ面全体の三次元形状を示す上フランジ面の形状データを求め、
     前記下フランジ面の形状データを用いて、前記下対象中点位置の有効三次元座標データを求めると共に、前記上フランジ面の形状データを用いて、前記上対象中点位置の有効三次元座標データを求める、
     回転機械のフランジ変位量推定方法。
  14.  請求項13に記載の回転機械のフランジ変位量推定方法において、
     前記有効座標把握工程では、前記下フランジ面の形状データを用いて、前記下第一位置及び前記下第二位置における有効三次元座標データを求めると共に、前記上フランジ面の形状データを用いて、前記上第一位置及び前記上第二位置における有効三次元座標データを求める、
     回転機械のフランジ変位量推定方法。
  15.  請求項13に記載の回転機械のフランジ変位量推定方法において、
     前記有効座標把握工程では、前記下フランジ面の形状データを用いて、前記下対象位置における有効三次元座標データを求めると共に、前記上フランジ面の形状データを用いて、前記上対象位置における有効三次元座標データを求める、
     回転機械のフランジ変位量推定方法。
  16.  水平方向に延びる軸線を中心として回転可能なロータと、
     前記ロータの外周を覆うケーシングと、
     前記ケーシング内に配置され、前記ケーシングに取り付けられている静止部品と、
     前記ケーシングを下側から支える架台と、
     を備え、
     前記ケーシングは、上側の上半ケーシングと、下側の下半ケーシングと、前記上半ケーシングと前記下半ケーシングとを締結する複数のボルトと、を有し、
     前記上半ケーシングは、下側を向く上フランジ面が形成されている上フランジを有し、
     前記下半ケーシングは、上側を向き、前記上フランジ面と上下方向で対向する下フランジ面が形成されている下フランジと、前記下フランジに連なり、前記架台により下側から支えられ、前記軸線が延びる軸線方向で互いに離れた第一被支持部及び第二被支持部と、を有し、
     前記上フランジ及び前記下フランジには、上下方向に貫通して、前記複数のボルトのそれぞれが挿通可能なボルト孔が形成されている、
     回転機械のフランジ変位量推定プログラムにおいて、
     前記回転機械を分解した後で、前記上半ケーシングと前記下半ケーシングとが前記複数のボルトで締結されていない開放状態における、前記上フランジ面中の複数位置における実測三次元座標データ及び前記下フランジ面中の複数位置における実測三次元座標データを受け付ける実測座標受付工程と、
     前記下フランジ面中の複数位置における前記実測三次元座標データを用いて、下第一位置と、下第二位置と、下対象位置と、下対象中点位置と、における有効三次元座標データを把握すると共に、前記上フランジ面中の複数位置における前記実測三次元座標データを用いて、上第一位置と、上第二位置と、上対象位置と、上対象中点位置と、における有効三次元座標データを把握する有効座標把握工程と、
     前記有効座標把握工程で把握した前記下第一位置の有効三次元座標データと前記上第一位置の有効三次元座標データとが一致し、前記有効座標把握工程で把握した前記下第二位置の有効三次元座標データと前記上第二位置の有効三次元座標データとが一致するよう、前記有効座標把握工程で把握した有効三次元座標データを変更する座標変更工程と、
     前記座標変更工程で変更された、前記下対象中点位置及び前記上対象中点位置における有効三次元座標データを用いて、前記下対象中点位置と前記上対象中点位置との上下方向における中間位置である対象接触位置の有効三次元座標データを求める接触位置推定工程と、
     前記開放状態から前記上半ケーシングと前記下半ケーシングとが前記複数のボルトで締結された締結状態になったときの前記上対象位置及び前記下対象位置の上下方向の変位量を求める変位量演算工程と、
     をコンピュータに実行させ、
     前記下第一位置は、前記上フランジ面に連なる面中で前記第一被支持部の第一代表位置と水平方向における位置が一致している位置であり、
     前記下第二位置は、前記下フランジ面に連なる面中で前記第二被支持部の第二代表位置と水平方向における位置が一致している位置であり、
     前記下対象位置は、前記下フランジ面中で、前記開放状態から前記締結状態になったときの上下方向の変位量を得たい位置であり、
     前記下対象中点位置は、前記下フランジ面中で、水平方向で前記軸線方向に垂直な方向である横方向における中点の位置であって、前記下対象位置と前記軸線方向の位置が一致している位置であり、
     前記上第一位置は、前記上フランジ面に連なる面中で前記第一被支持部の前記第一代表位置と水平方向における位置が一致している位置であり、
     前記上第二位置は、前記上フランジ面に連なる面中で前記第二被支持部の前記第二代表位置と水平方向における位置が一致している位置であり、
     前記上対象位置は、前記上フランジ面中で前記下対象位置と水平方向の位置が一致している位置であり、
     前記上対象中点位置は、前記上フランジ面中で前記横方向における中点であって前記下対象位置と前記軸線方向の位置が一致している位置であり、
     前記変位量演算工程では、
     前記座標変更工程で変更された前記下対象位置の有効三次元座標データが示す上下方向の位置と前記対象接触位置の有効三次元座標データが示す上下方向の位置との差を前記下対象位置の上下方向の変位量とし、前記座標変更工程で変更された前記上対象位置の有効三次元座標データが示す上下方向の位置と前記対象接触位置の有効三次元座標データが示す上下方向の位置との差を前記上対象位置の上下方向の変位量とする、
     回転機械のフランジ変位量推定プログラム。
  17.  請求項16に記載の回転機械のフランジ変位量推定プログラムにおいて、
     前記下対象位置は、前記軸線方向で前記静止部品が配置される位置で且つ前記下フランジ面中の内側縁の位置である、
     回転機械のフランジ変位量推定プログラム。
  18.  請求項16又は17に記載の回転機械のフランジ変位量推定プログラムにおいて、
     前記実測座標受付工程では、前記下フランジ面中で前記横方向における中点の位置であって、前記軸線方向における位置が互いに異なる複数の下中点位置と、前記上フランジ面中で前記横方向における中点の位置であって、前記軸線方向における位置が互いに異なる複数の上中点位置と、における実測三次元座標データを受け付け、
     前記有効座標把握工程では、前記複数の下中点位置における実測三次元座標データの変化傾向から、前記下対象中点位置における有効三次元座標データを求め、前記複数の上中点位置における有効三次元座標データの変化傾向から、前記上対象中点位置における有効三次元座標データを求める、
     回転機械のフランジ変位量推定プログラム。
  19.  請求項18に記載の回転機械のフランジ変位量推定プログラムにおいて、
     前記有効座標把握工程では、前記複数の下中点位置における有効三次元座標データの変化傾向から、前記下第一位置及び前記下第二位置における有効三次元座標データを求め、前記複数の上中点位置における有効三次元座標データの変化傾向から、前記上第一位置及び前記上第二位置における有効三次元座標データを求める、
     回転機械のフランジ変位量推定プログラム。
  20.  水平方向に延びる軸線を中心として回転可能なロータと、
     前記ロータの外周を覆うケーシングと、
     前記ケーシング内に配置され、前記ケーシングに取り付けられている静止部品と、
     前記ケーシングを下側から支える架台と、
     を備え、
     前記ケーシングは、上側の上半ケーシングと、下側の下半ケーシングと、前記上半ケーシングと前記下半ケーシングとを締結する複数のボルトと、を有し、
     前記上半ケーシングは、下側を向く上フランジ面が形成されている上フランジを有し、
     前記下半ケーシングは、上側を向き、前記上フランジ面と上下方向で対向する下フランジ面が形成されている下フランジと、前記下フランジに連なり、前記架台により下側から支えられ、前記軸線が延びる軸線方向で互いに離れた第一被支持部及び第二被支持部と、を有し、
     前記上フランジ及び前記下フランジには、上下方向に貫通して、前記複数のボルトのそれぞれが挿通可能なボルト孔が形成されている、
     回転機械のフランジ変位量推定装置において、
     前記回転機械を分解した後で、前記上半ケーシングと前記下半ケーシングとが前記複数のボルトで締結されていない開放状態における、前記上フランジ面中の複数位置における実測三次元座標データ及び前記下フランジ面中の複数位置における実測三次元座標データを受け付ける実測座標受付部と、
     前記下フランジ面中の複数位置における前記実測三次元座標データを用いて、下第一位置と、下第二位置と、下対象位置と、下対象中点位置と、における有効三次元座標データを把握すると共に、前記上フランジ面中の複数位置における前記実測三次元座標データを用いて、上第一位置と、上第二位置と、上対象位置と、上対象中点位置と、における有効三次元座標データを把握する有効座標把握部と、
     前記有効座標把握部が把握した前記下第一位置の有効三次元座標データと前記上第一位置の有効三次元座標データとが一致し、前記有効座標把握部が把握した前記下第二位置の有効三次元座標データと前記上第二位置の有効三次元座標データとが一致するよう、前記有効座標把握部が把握した有効三次元座標データを変更する座標変更部と、
     前記座標変更部が変更した、前記下対象中点位置及び前記上対象中点位置における有効三次元座標データを用いて、前記下対象中点位置と前記上対象中点位置との上下方向における中間位置である対象接触位置の有効三次元座標データを求める接触位置推定部と、
     前記開放状態から前記上半ケーシングと前記下半ケーシングとが前記複数のボルトで締結された締結状態になったときの前記上対象位置及び前記下対象位置の上下方向の変位量を求める変位量演算部と、
     を備え、
     前記下第一位置は、前記上フランジ面に連なる面中で前記第一被支持部の第一代表位置と水平方向における位置が一致している位置であり、
     前記下第二位置は、前記下フランジ面に連なる面中で前記第二被支持部の第二代表位置と水平方向における位置が一致している位置であり、
     前記下対象位置は、前記下フランジ面中で、前記開放状態から前記締結状態になったときの上下方向の変位量を得たい位置であり、
     前記下対象中点位置は、前記下フランジ面中で、水平方向で前記軸線方向に垂直な方向である横方向における中点の位置であって、前記下対象位置と前記軸線方向の位置が一致している位置であり、
     前記上第一位置は、前記上フランジ面に連なる面中で前記第一被支持部の前記第一代表位置と水平方向における位置が一致している位置であり、
     前記上第二位置は、前記上フランジ面に連なる面中で前記第二被支持部の前記第二代表位置と水平方向における位置が一致している位置であり、
     前記上対象位置は、前記上フランジ面中で前記下対象位置と水平方向の位置が一致している位置であり、
     前記上対象中点位置は、前記上フランジ面中で前記横方向における中点であって前記下対象位置と前記軸線方向の位置が一致している位置であり、
     前記変位量演算部は、前記座標変更部が変更した前記下対象位置の有効三次元座標データが示す上下方向の位置と前記対象接触位置の有効三次元座標データが示す上下方向の位置との差を前記下対象位置の上下方向の変位量とし、前記座標変更部が変更した前記上対象位置の有効三次元座標データが示す上下方向の位置と前記対象接触位置の有効三次元座標データが示す上下方向の位置との差を前記上対象位置の上下方向の変位量とする、
     回転機械のフランジ変位量推定装置。
  21.  請求項20に記載の回転機械のフランジ変位量推定装置において、
     前記下対象位置は、前記軸線方向で前記静止部品が配置される位置で且つ前記下フランジ面中の内側縁の位置である、
     回転機械のフランジ変位量推定装置。
  22.  請求項20又は21に記載の回転機械のフランジ変位量推定装置において、
     前記実測座標受付部は、前記下フランジ面中で前記横方向における中点の位置であって、前記軸線方向における位置が互いに異なる複数の下中点位置と、前記上フランジ面中で前記横方向における中点の位置であって、前記軸線方向における位置が互いに異なる複数の上中点位置と、における実測三次元座標データを受け付け、
     前記有効座標把握部は、前記複数の下中点位置における実測三次元座標データの変化傾向から、前記下対象中点位置における有効三次元座標データを求め、前記複数の上中点位置における有効三次元座標データの変化傾向から、前記上対象中点位置における有効三次元座標データを求める、
     回転機械のフランジ変位量推定装置。
  23.  請求項22に記載の回転機械のフランジ変位量推定装置において、
     前記有効座標把握部は、前記複数の下中点位置における有効三次元座標データの変化傾向から、前記下第一位置及び前記下第二位置における有効三次元座標データを求め、前記複数の上中点位置における有効三次元座標データの変化傾向から、前記上第一位置及び前記上第二位置における有効三次元座標データを求める、
     回転機械のフランジ変位量推定装置。
PCT/JP2022/044204 2022-02-25 2022-11-30 回転機械のフランジ変位量推定方法、この方法を実行するためのプログラム、及び、この方法を実行する装置 WO2023162386A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2024502837A JPWO2023162386A1 (ja) 2022-02-25 2022-11-30
KR1020237043921A KR20240010039A (ko) 2022-02-25 2022-11-30 회전 기계의 플랜지 변위량 추정 방법, 이 방법을 실행하기 위한 프로그램 및, 이 방법을 실행하는 장치
CN202280045297.6A CN117561368A (zh) 2022-02-25 2022-11-30 旋转机械的凸缘位移量推定方法、用于执行该方法的程序以及执行该方法的装置

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2022027442 2022-02-25
JP2022-027442 2022-02-25

Publications (1)

Publication Number Publication Date
WO2023162386A1 true WO2023162386A1 (ja) 2023-08-31

Family

ID=87765499

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2022/044204 WO2023162386A1 (ja) 2022-02-25 2022-11-30 回転機械のフランジ変位量推定方法、この方法を実行するためのプログラム、及び、この方法を実行する装置

Country Status (4)

Country Link
JP (1) JPWO2023162386A1 (ja)
KR (1) KR20240010039A (ja)
CN (1) CN117561368A (ja)
WO (1) WO2023162386A1 (ja)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03206303A (ja) * 1989-07-28 1991-09-09 Westinghouse Electric Corp <We> 接合部の製作方法
JP2004516415A (ja) * 2000-12-21 2004-06-03 ゼネラル・エレクトリック・カンパニイ 蒸気タービンのハーフシェルを位置合わせする方法
JP2018084169A (ja) * 2016-11-22 2018-05-31 三菱日立パワーシステムズ株式会社 タービンの組立方法
WO2021086208A1 (en) * 2019-10-28 2021-05-06 General Electric Company Method and system for component alignment in turbine casing and related turbine casing
JP2022037334A (ja) * 2020-08-25 2022-03-09 三菱重工業株式会社 タービンの組立方法、タービンの組立支援プログラム、及びタービンの組立支援装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6864596B2 (ja) 2017-10-06 2021-04-28 三菱パワー株式会社 タービン組立支援プログラム、タービン組立支援システム及びタービンの組立方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03206303A (ja) * 1989-07-28 1991-09-09 Westinghouse Electric Corp <We> 接合部の製作方法
JP2004516415A (ja) * 2000-12-21 2004-06-03 ゼネラル・エレクトリック・カンパニイ 蒸気タービンのハーフシェルを位置合わせする方法
JP2018084169A (ja) * 2016-11-22 2018-05-31 三菱日立パワーシステムズ株式会社 タービンの組立方法
WO2021086208A1 (en) * 2019-10-28 2021-05-06 General Electric Company Method and system for component alignment in turbine casing and related turbine casing
JP2022037334A (ja) * 2020-08-25 2022-03-09 三菱重工業株式会社 タービンの組立方法、タービンの組立支援プログラム、及びタービンの組立支援装置

Also Published As

Publication number Publication date
CN117561368A (zh) 2024-02-13
JPWO2023162386A1 (ja) 2023-08-31
KR20240010039A (ko) 2024-01-23

Similar Documents

Publication Publication Date Title
KR101965507B1 (ko) 가상 블레이드 검사
JP6470954B2 (ja) ターボ機械
EP3130887B1 (en) Virtual turbomachine blade contact gap inspection
JP7416674B2 (ja) タービンの組立方法、タービンの組立支援プログラム、及びタービンの組立支援装置
Rybczyński The possibility of evaluating turbo-set bearing misalignment defects on the basis of bearing trajectory features
KR101741885B1 (ko) 회전 로터의 현재 편심을 결정하기 위한 방법 및 회전 로터의 편심 진단 방법
JP4772594B2 (ja) 回転機器の低速バランス法及び低速バランス実施装置
Millecamps et al. Snecma’s viewpoint on the numerical and experimental simulation of blade-tip/casing unilateral contacts
WO2023162386A1 (ja) 回転機械のフランジ変位量推定方法、この方法を実行するためのプログラム、及び、この方法を実行する装置
WO2023162384A1 (ja) 回転機械のフランジ変位量推定方法、この方法を実行するためのプログラム、及び、この方法を実行する装置
WO2023162388A1 (ja) 回転機械のフランジ変位量推定方法、この方法を実行するためのプログラム、及び、この方法を実行する装置
Wang et al. Reduced-order modeling for rotating rotor-bearing systems with cracked impellers using three-dimensional finite element models
CN111695206A (zh) 一种旋转态叶盘的失谐识别和模型更新方法
WO2023162387A1 (ja) 回転機械におけるフランジ面圧分布の推定方法、フランジ面間からの流体のリーク評価方法、これらの方法を実行するためのプログラム及び装置
KR102680961B1 (ko) 터빈의 조립 방법, 터빈의 조립 지원 프로그램, 및 터빈의 조립 지원 장치
Racic et al. Practical approach for solving vibrations of large turbine and generator rotors-Reconciling the discord between theory and practice
CN110231161A (zh) 基于三维坐标系的单跨度转子故障位置诊断方法及系统
CN114867932B (zh) 涡轮机的组装方法、涡轮机的组装辅助程序以及涡轮机的组装辅助装置
Maturkanič et al. Analytical method “Anmena” for blade tip timing measurement systems–Principle of the method
Katinić et al. Unbalance Response Analysis of a Rotor Kit with Two Identical Discs Located Between Bearings
Aguirre Reduced order models for nonlinear dynamic analysis of rotating structures: Application to turbomachinery blades| Theses. fr
TR2023019793T2 (tr) Döner makinede flanş yer değiştirme miktarının tahminine yönelik yöntem, yöntemin yürütülmesi için program ve yöntemin gerçekleştirilmesine yönelik cihaz.
Okabe et al. Rotor-blade coupled torsional vibration analysis using modal parameters based on fem analyses and experiments
CN116205107A (zh) 一种高速薄壁回转结构不平衡量计算方法及相关装置
CN114935430A (zh) 一种仅采用手持测振仪进行动平衡加重点位置的计算方法

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 22928903

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2024502837

Country of ref document: JP

ENP Entry into the national phase

Ref document number: 20237043921

Country of ref document: KR

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 1020237043921

Country of ref document: KR

WWE Wipo information: entry into national phase

Ref document number: 202280045297.6

Country of ref document: CN

WWE Wipo information: entry into national phase

Ref document number: MX/A/2024/000596

Country of ref document: MX

WWE Wipo information: entry into national phase

Ref document number: 2024/002451

Country of ref document: TR