JP6578250B2 - Turbine component strain estimation method and apparatus, and turbine component evaluation method - Google Patents

Turbine component strain estimation method and apparatus, and turbine component evaluation method Download PDF

Info

Publication number
JP6578250B2
JP6578250B2 JP2016122021A JP2016122021A JP6578250B2 JP 6578250 B2 JP6578250 B2 JP 6578250B2 JP 2016122021 A JP2016122021 A JP 2016122021A JP 2016122021 A JP2016122021 A JP 2016122021A JP 6578250 B2 JP6578250 B2 JP 6578250B2
Authority
JP
Japan
Prior art keywords
turbine component
deformation
shape
strain
turbine
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2016122021A
Other languages
Japanese (ja)
Other versions
JP2017227137A (en
Inventor
浩平 羽田野
浩平 羽田野
貴洋 宮本
貴洋 宮本
豪 李
豪 李
橋本 幸弘
幸弘 橋本
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Mitsubishi Heavy Industries Ltd
Original Assignee
Mitsubishi Heavy Industries Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Mitsubishi Heavy Industries Ltd filed Critical Mitsubishi Heavy Industries Ltd
Priority to JP2016122021A priority Critical patent/JP6578250B2/en
Publication of JP2017227137A publication Critical patent/JP2017227137A/en
Application granted granted Critical
Publication of JP6578250B2 publication Critical patent/JP6578250B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Turbine Rotor Nozzle Sealing (AREA)

Description

本発明は、蒸気タービンやガスタービンなどにタービンコンポーネントとして使用される部材における非弾性ひずみを推定するタービンコンポーネントのひずみ推定方法及び装置、タービンコンポーネントの評価方法に関するものである。   The present invention relates to a turbine component strain estimation method and apparatus for estimating inelastic strain in a member used as a turbine component in a steam turbine, a gas turbine, or the like, and a turbine component evaluation method.

例えば、ガスタービンは、圧縮機と燃焼器とタービンにより構成されている。そして、空気取入口から取り込まれた空気が圧縮機によって圧縮されることで高温・高圧の圧縮空気となり、燃焼器にて、この圧縮空気に対して燃料を供給して燃焼させることで高温・高圧の燃焼ガス(作動流体)を得て、この燃焼ガスによりタービンを駆動し、このタービンに連結された発電機を駆動する。このガスタービンにて、圧縮機やタービンに使用される動翼は、回転時に作用する遠心力や熱の影響により損傷を受ける。そのため、使用された動翼に対して損傷の予測解析を実施し、損傷の予測やき裂の有無などの点検調査などを実施することで健全性の維持に努めている。   For example, a gas turbine includes a compressor, a combustor, and a turbine. The air taken in from the air intake port is compressed by the compressor to become high-temperature and high-pressure compressed air. In the combustor, the fuel is supplied to the compressed air and burned, so that the high-temperature and high-pressure is burned The combustion gas (working fluid) is obtained, the turbine is driven by the combustion gas, and the generator connected to the turbine is driven. In this gas turbine, the rotor blades used in the compressor and turbine are damaged by the influence of centrifugal force and heat acting during rotation. For this reason, we are striving to maintain soundness by conducting predictive analysis of damage to used rotor blades and conducting inspections such as predicting damage and the presence or absence of cracks.

動翼に対するクリープ変形量で損傷を評価する手法として、例えば、下記特許文献1に記載されたものがある。この特許文献1に記載されたタービンコンポーネントのクリープ能力を決定するための方法は、複数のタービンコンポーネントに引張応力、遠心応力、熱応力を測定可能な量のクリープが得られるまで加え、タービンコンポーネントで使用される材料に対する知られているクリープ特性に対して比較し、定義済みの量より大きいクリープを示すコンポーネントを分離して運用から外すものである。   As a method for evaluating damage by the amount of creep deformation with respect to a moving blade, for example, there is one described in Patent Document 1 below. The method for determining the creep capability of a turbine component described in Patent Document 1 is applied to a plurality of turbine components until a measurable amount of tensile stress, centrifugal stress, and thermal stress is obtained. Compared to known creep properties for the materials used, components that exhibit a creep greater than a defined amount are separated and removed from service.

特開2013−253599号公報JP2013-253599A

上述した従来のタービンコンポーネントのクリープ能力を決定するための方法では、試験装置で動翼に対して引張応力、遠心応力、熱応力を与えることでクリープ変形を誘発させ、このときに得られたひずみ量またはその蓄積速度を、クリープ特性を有する材料から定義された公差値と比較して動翼を評価している。ところが、試験装置により各応力が与えられた動翼と、実際のプラントでの運転中に各種の応力が作用した動翼とでは、ひずみ量に差があり、動翼のクリープ損傷を高精度に予測することは困難である。   In the conventional method for determining the creep capability of the turbine component, creep deformation is induced by applying tensile stress, centrifugal stress, and thermal stress to the rotor blade with a test device. The blade is evaluated by comparing the quantity or its accumulation rate with a tolerance value defined from a material having creep properties. However, there is a difference in the amount of strain between the moving blades that have been subjected to various stresses by the test equipment and the moving blades that have been subjected to various stresses during actual plant operation, and the creep damage of the moving blades is highly accurate. It is difficult to predict.

本発明は上述した課題を解決するものであり、タービンコンポーネントに発生した非弾性ひずみを精度良く推定することができるタービンコンポーネントのひずみ推定方法及び装置、タービンコンポーネントの評価方法を提供することを目的とする。   SUMMARY OF THE INVENTION An object of the present invention is to solve the above-described problems, and to provide a turbine component strain estimation method and apparatus and a turbine component evaluation method capable of accurately estimating inelastic strain generated in a turbine component. To do.

上記の目的を達成するための本発明のタービンコンポーネントのひずみ推定方法は、変形前のタービンコンポーネントを3次元計測するステップと、前記タービンコンポーネントの解析モデルを作成するステップと、前記変形前のタービンコンポーネントが分割された複数のメッシュを所定の大きさの単位領域に区画するステップと、前記解析モデルを用いて前記単位領域の前記複数のメッシュに対して一様な非弾性ひずみを入力して前記タービンコンポーネントの解析モデル形状を形成するステップと、変形後のタービンコンポーネントを3次元計測するステップと、前記タービンコンポーネントの解析モデル形状と変形後の前記タービンコンポーネントの形状とを比較するステップと、前記タービンコンポーネントの解析モデル形状と変形後の前記タービンコンポーネントの形状との比較結果が収束条件を満たさないときに前記解析モデルに入力する非弾性ひずみを変更するステップと、前記タービンコンポーネントの解析モデル形状と変形後の前記タービンコンポーネントの形状との比較結果が収束条件を満たすときに前記解析モデルに入力した非弾性ひずみを確定するステップと、を有することを特徴とするものである。   In order to achieve the above object, a turbine component strain estimation method according to the present invention includes a step of three-dimensionally measuring a turbine component before deformation, a step of creating an analysis model of the turbine component, and the turbine component before deformation. Dividing the plurality of meshes divided into unit regions of a predetermined size, and inputting the uniform inelastic strain to the plurality of meshes in the unit regions using the analysis model Forming an analytical model shape of the component; three-dimensionally measuring the turbine component after deformation; comparing the analytical model shape of the turbine component with the shape of the turbine component after deformation; and the turbine component Analytical model shape and deformation A step of changing an inelastic strain input to the analysis model when a result of comparison with the shape of the subsequent turbine component does not satisfy a convergence condition; and an analysis model shape of the turbine component and a shape of the turbine component after deformation And determining the inelastic strain input to the analysis model when the comparison result satisfies the convergence condition.

従って、変形前のタービンコンポーネントが分割された複数のメッシュを所定の大きさの単位領域に区画し、解析モデルを用いて単位領域の各メッシュに一様な非弾性ひずみを入力してタービンコンポーネントの解析モデル形状を形成し、変形後のタービンコンポーネントと比較する。ここで、比較結果が収束条件を満たさないときは、解析モデルに入力する非弾性ひずみを繰り返し変更することで、比較結果が収束条件を満たすようにし、収束条件を満たした非弾性ひずみを確定して出力する。そのため、単位領域に一様な非弾性ひずみを入力することで、逆解析時における変数(パラメータ)の数が減少して発散を防止することができると共に、自由に非弾性ひずみを入力して解析を行うことでタービンコンポーネントに発生した非弾性ひずみを精度良く推定することができる。   Therefore, a plurality of meshes obtained by dividing the turbine component before deformation are divided into unit regions of a predetermined size, and uniform inelastic strain is input to each mesh in the unit region using an analysis model. An analytical model shape is formed and compared with the turbine component after deformation. If the comparison result does not satisfy the convergence condition, the inelastic strain input to the analysis model is repeatedly changed so that the comparison result satisfies the convergence condition, and the inelastic strain that satisfies the convergence condition is determined. Output. Therefore, by inputting uniform inelastic strain to the unit area, the number of variables (parameters) at the time of inverse analysis can be reduced and divergence can be prevented, and analysis can be performed by freely inputting inelastic strain. This makes it possible to accurately estimate the inelastic strain generated in the turbine component.

本発明のタービンコンポーネントのひずみ推定方法では、前記単位領域は、前記タービンコンポーネントに対して複数設定し、複数の前記単位領域ごとに異なる非弾性ひずみを入力することを特徴としている。   The turbine component strain estimation method of the present invention is characterized in that a plurality of unit regions are set for the turbine component, and different inelastic strains are input to the plurality of unit regions.

従って、設定された複数の単位領域に対して異なる大きさの非弾性ひずみを入力することで、タービンコンポーネントの変形に応じた適正な解析モデル形状を形成することができる。   Accordingly, by inputting inelastic strains having different sizes to the plurality of set unit regions, it is possible to form an appropriate analysis model shape corresponding to the deformation of the turbine component.

本発明のタービンコンポーネントのひずみ推定方法では、前記解析モデルを用いて入力する一様な非弾性ひずみとは、同一の非弾性ひずみであることを特徴としている。   The turbine component strain estimation method of the present invention is characterized in that the uniform inelastic strain input using the analysis model is the same inelastic strain.

従って、1つの単位領域における複数のメッシュに同一の大きさの非弾性ひずみを入力することで、変数(パラメータ)の数が減少して発散を防止することができる。   Therefore, by inputting the same inelastic strain to a plurality of meshes in one unit region, the number of variables (parameters) can be reduced and divergence can be prevented.

本発明のタービンコンポーネントのひずみ推定方法では、前記解析モデルを用いて入力する一様な非弾性ひずみとは、線形に変化する非弾性ひずみであることを特徴としている。   The turbine component strain estimation method of the present invention is characterized in that the uniform inelastic strain input using the analysis model is a linearly varying inelastic strain.

従って、1つの単位領域における複数のメッシュに線形に変化する非弾性ひずみを入力することで、変数(パラメータ)の数が減少して発散を防止することができる。   Therefore, by inputting linearly changing inelastic strain to a plurality of meshes in one unit region, the number of variables (parameters) can be reduced and divergence can be prevented.

本発明のタービンコンポーネントのひずみ推定方法では、前記収束条件は、前記タービンコンポーネントの解析モデル形状と3次元計測した変形後の前記タービンコンポーネントの形状との変形差が予め設定された規定値以内にあることを特徴としている。   In the turbine component strain estimation method of the present invention, the convergence condition is such that the deformation difference between the analysis model shape of the turbine component and the shape of the turbine component after deformation measured three-dimensionally is within a preset specified value. It is characterized by that.

従って、解析モデル形状と変形後のタービンコンポーネントの形状との変形差が規定値以内にあるときに、収束条件を満たして解析モデルに入力した非弾性ひずみを確定することで、タービンコンポーネントにマークを付けたり、加工を施す必要がなく、容易に変形差を比較することができる。   Therefore, when the deformation difference between the shape of the analytical model and the shape of the turbine component after deformation is within the specified value, the inelastic strain that satisfies the convergence condition and is input to the analytical model is determined, thereby marking the turbine component. There is no need to attach or process, and the difference in deformation can be easily compared.

本発明のタービンコンポーネントのひずみ推定方法では、前記変形差とは、前記タービンコンポーネントの解析モデルの体積形状と3次元計測した変形後の前記タービンコンポーネントの体積形状とのずれ体積量であることを特徴としている。   In the turbine component strain estimation method of the present invention, the deformation difference is a displacement volume amount between a volume shape of an analysis model of the turbine component and a volume shape of the turbine component after deformation measured three-dimensionally. It is said.

従って、変形差を解析モデルの体積形状と変形後の前記タービンコンポーネントの体積形状とのずれ体積量とすることで、高精度に変形差を求めることができる。   Therefore, the deformation difference can be obtained with high accuracy by setting the deformation difference as the displacement volume amount between the volume shape of the analysis model and the volume shape of the turbine component after the deformation.

また、本発明のタービンコンポーネントのひずみ推定方法は、変形前のタービンコンポーネントを3次元計測するステップと、前記タービンコンポーネントの解析モデルを作成するステップと、前記解析モデルに複数種類の変形条件を入力して前記タービンコンポーネントの複数の解析モデル形状を形成するステップと、変形後のタービンコンポーネントを3次元計測するステップと、前記タービンコンポーネントの複数の解析モデル形状と変形後の前記タービンコンポーネントの形状とを比較するステップと、前記タービンコンポーネントの複数の解析モデル形状と変形後の前記タービンコンポーネントの形状との変形差が最小となる前記解析モデル形状に入力した非弾性ひずみを確定するステップと、を有することを特徴とするものである。   In the turbine component strain estimation method of the present invention, a step of three-dimensionally measuring a turbine component before deformation, a step of creating an analysis model of the turbine component, and inputting a plurality of types of deformation conditions into the analysis model. Forming a plurality of analytical model shapes of the turbine component, measuring a three-dimensional measurement of the turbine component after deformation, and comparing the plurality of analytical model shapes of the turbine component and the shape of the turbine component after deformation And determining an inelastic strain input to the analysis model shape that minimizes a deformation difference between the plurality of analysis model shapes of the turbine component and the shape of the turbine component after deformation. Features .

従って、解析モデルに複数種類の変形条件を入力してタービンコンポーネントの複数の解析モデル形状を形成し、変形後のタービンコンポーネントの形状と比較する。ここで、タービンコンポーネントの複数の解析モデル形状と変形後のタービンコンポーネントの形状との変形差が最小となる非弾性ひずみを確定して出力する。そのため、逆解析時における変数(パラメータ)の数が減少して発散を防止することができると共に、自由に非弾性ひずみを入力して解析を行うことでタービンコンポーネントに発生した非弾性ひずみを精度良く推定することができる。   Accordingly, a plurality of types of deformation conditions are input to the analysis model to form a plurality of analysis model shapes of the turbine component and compared with the shape of the turbine component after the deformation. Here, the inelastic strain that minimizes the deformation difference between the plurality of analysis model shapes of the turbine component and the shape of the turbine component after the deformation is determined and output. Therefore, the number of variables (parameters) at the time of inverse analysis can be reduced to prevent divergence, and the inelastic strain generated in the turbine component can be accurately analyzed by inputting inelastic strain freely and performing analysis. Can be estimated.

また、本発明のタービンコンポーネントのひずみ推定装置は、タービンコンポーネントを3次元計測する計測部と、前記計測部が計測した変形前の前記タービンコンポーネントが分割された複数のメッシュを所定の大きさの単位領域に区画し、解析モデルを用いて前記単位領域の前記複数のメッシュに対して一様な非弾性ひずみを入力して前記タービンコンポーネントの解析モデル形状を求める演算部と、前記タービンコンポーネントの解析モデル形状と前記計測部が計測した変形後の前記タービンコンポーネントの形状とを比較する比較部と、前記比較部の比較結果が収束条件を満たさないときに前記演算部の処理を繰り返す一方、前記比較部の比較結果が収束条件を満たすときに解析モデルに入力する非弾性ひずみを確定する判定部と、を備えることを特徴とするものである。   In addition, the turbine component strain estimation apparatus of the present invention includes a measurement unit that three-dimensionally measures a turbine component, and a plurality of meshes obtained by dividing the turbine component before deformation measured by the measurement unit in units of a predetermined size. A calculation unit that divides into regions and inputs uniform inelastic strain to the plurality of meshes of the unit region using an analysis model to obtain an analysis model shape of the turbine component; and an analysis model of the turbine component A comparison unit that compares the shape and the shape of the turbine component after deformation measured by the measurement unit, and the comparison unit repeats the processing of the calculation unit when a comparison result of the comparison unit does not satisfy a convergence condition. A determination unit that determines the inelastic strain to be input to the analysis model when the comparison result of It is characterized in further comprising.

従って、変形前のタービンコンポーネントが分割された複数のメッシュを所定の大きさの単位領域に区画し、解析モデルを用いて単位領域の各メッシュに一様な非弾性ひずみを入力してタービンコンポーネントの解析モデル形状を形成し、変形後のタービンコンポーネントと比較する。ここで、比較結果が収束条件を満たさないときは、解析モデルに入力する非弾性ひずみを繰り返し変更することで、比較結果が収束条件を満たすようにし、収束条件を満たした非弾性ひずみを確定して出力する。そのため、単位領域に一様な非弾性ひずみを入力することで、逆解析時における変数(パラメータ)の数が減少して発散を防止することができると共に、自由に非弾性ひずみを入力して解析を行うことでタービンコンポーネントに発生した非弾性ひずみを精度良く推定することができる。   Therefore, a plurality of meshes obtained by dividing the turbine component before deformation are divided into unit regions of a predetermined size, and uniform inelastic strain is input to each mesh in the unit region using an analysis model. An analytical model shape is formed and compared with the turbine component after deformation. If the comparison result does not satisfy the convergence condition, the inelastic strain input to the analysis model is repeatedly changed so that the comparison result satisfies the convergence condition, and the inelastic strain that satisfies the convergence condition is determined. Output. Therefore, by inputting uniform inelastic strain to the unit area, the number of variables (parameters) at the time of inverse analysis can be reduced and divergence can be prevented, and analysis can be performed by freely inputting inelastic strain. This makes it possible to accurately estimate the inelastic strain generated in the turbine component.

また、本発明のタービンコンポーネントの評価方法は、前記タービンコンポーネントのひずみ推定方法により推定した非弾性ひずみを用いて前記タービンコンポーネントを評価することを特徴とするものである。   The turbine component evaluation method of the present invention is characterized in that the turbine component is evaluated using inelastic strain estimated by the turbine component strain estimation method.

従って、逆解析時における変数(パラメータ)の数が減少して発散を防止することができると共に、自由に非弾性ひずみを入力して解析を行うことでタービンコンポーネントに発生した非弾性ひずみを精度良く推定することができ、その結果、タービンコンポーネントの評価を精度良く実施することができる。   Therefore, the number of variables (parameters) at the time of inverse analysis can be reduced and divergence can be prevented, and the inelastic strain generated in the turbine component can be accurately analyzed by inputting inelastic strain freely and performing analysis. As a result, the turbine component can be evaluated with high accuracy.

本発明のタービンコンポーネントのひずみ推定方法及び装置、タービンコンポーネントの評価方法によれば、逆解析時における変数(パラメータ)の数が減少して発散を防止することができると共に、自由に非弾性ひずみを入力して解析を行うことでタービンコンポーネントに発生した非弾性ひずみを精度良く推定することができる。   According to the turbine component strain estimation method and apparatus and turbine component evaluation method of the present invention, the number of variables (parameters) at the time of inverse analysis can be reduced to prevent divergence, and inelastic strain can be freely reduced. By inputting and analyzing, inelastic strain generated in the turbine component can be accurately estimated.

図1は、第1実施形態のタービン動翼のひずみ推定装置を表すブロック構成図である。FIG. 1 is a block configuration diagram illustrating a turbine rotor blade strain estimation apparatus according to a first embodiment. 図2は、使用後のタービン動翼の評価方法を表すフローチャートである。FIG. 2 is a flowchart showing a method for evaluating a turbine rotor blade after use. 図3は、変形処理後のタービン動翼の評価方法を表すフローチャートである。FIG. 3 is a flowchart showing a method for evaluating the turbine rotor blade after the deformation process. 図4は、第1実施形態のタービン動翼のひずみ推定方法を表すフローチャートである。FIG. 4 is a flowchart showing the turbine blade distortion estimation method of the first embodiment. 図5は、タービン動翼のひずみ推定方法を説明するための概略図である。FIG. 5 is a schematic diagram for explaining a method for estimating the distortion of the turbine rotor blade. 図6は、第2実施形態のタービン動翼のひずみ推定方法を表すフローチャートである。FIG. 6 is a flowchart showing a turbine blade distortion estimation method according to the second embodiment.

以下に添付図面を参照して、本発明に係るタービンコンポーネントのひずみ推定方法の好適な実施形態を詳細に説明する。なお、この実施形態により本発明が限定されるものではなく、また、実施形態が複数ある場合には、各実施形態を組み合わせて構成するものも含むものである。   Exemplary embodiments of a turbine component strain estimation method according to the present invention will be described below in detail with reference to the accompanying drawings. In addition, this invention is not limited by this embodiment, and when there are two or more embodiments, what comprises combining each embodiment is also included.

[第1実施形態]
第1実施形態では、タービンコンポーネントとしてタービン動翼を用いて説明する。タービン動翼のひずみ推定方法を用いたタービン動翼の評価方法は、変形後のタービン動翼に対して3次元(3D)計測を実施し、3次元計測結果を用いた逆解析(FEM解析)を行うことで、タービン動翼における非弾性ひずみの発生量を予測し、タービン動翼を評価するものである。
[First Embodiment]
In the first embodiment, a turbine blade is used as a turbine component. The turbine blade evaluation method using the turbine blade distortion estimation method is to perform three-dimensional (3D) measurement on the deformed turbine blade and perform inverse analysis (FEM analysis) using the three-dimensional measurement result. By performing the above, the amount of inelastic strain generated in the turbine blade is predicted, and the turbine blade is evaluated.

タービン動翼が変形する要因としては、以下の3つの要因がある。第1の要因は、タービン動翼を実際のガスタービンに装着し、このガスタービンの運転中にタービン動翼に対して作用する各種の応力であり、この応力によりタービン動翼が変形して非弾性ひずみが発生する。第2の要因は、タービン動翼の製造過程で、形状修正のために行う曲げ加工(ひずみ修正)であり、このときにタービン動翼に塑性変形を与えるため残留応力が発生する。第3の要因は、タービン動翼の製造過程で、塑性変形による残留応力の付与を図るために行うショットピーニングであり、タービン動翼の全面に均一な残留応力を与えても翼が撓むが、このときのタービン動翼極表面に塑性変形を与えるため残留応力が発生する。   There are the following three factors that cause deformation of the turbine blade. The first factor is various stresses acting on the turbine blades during operation of the gas turbine when the turbine blades are mounted on an actual gas turbine. Elastic strain occurs. The second factor is a bending process (strain correction) performed for shape correction in the manufacturing process of the turbine blade, and residual stress is generated at this time because plastic deformation is applied to the turbine blade. The third factor is shot peening performed in the manufacturing process of the turbine rotor blade in order to apply residual stress due to plastic deformation, and the blade may bend even if uniform residual stress is applied to the entire surface of the turbine rotor blade. Residual stress is generated because plastic deformation is applied to the turbine blade pole surface at this time.

第1実施形態のタービン動翼のひずみ推定方法を用いたタービン動翼の評価方法は、上述した要因で変形したタービン動翼における非弾性ひずみの発生量を推定し、この推定結果からタービン動翼を評価するものである。   In the turbine blade evaluation method using the turbine blade distortion estimation method of the first embodiment, the amount of inelastic strain generated in the turbine blade deformed due to the above-described factors is estimated. Is to evaluate.

図1は、第1実施形態のタービン動翼のひずみ推定装置を表すブロック構成図である。   FIG. 1 is a block configuration diagram illustrating a turbine rotor blade strain estimation apparatus according to a first embodiment.

図1に示すように、第1実施形態のタービン動翼のひずみ推定装置10は、解析装置11を備えており、この解析装置11は、計測部12と、表示部13と、記憶部14とが接続されている。また、解析装置11は、演算部21と、比較部22と、判定部23とを備えている。   As shown in FIG. 1, the turbine blade distortion estimation device 10 of the first embodiment includes an analysis device 11, which includes a measurement unit 12, a display unit 13, a storage unit 14, and the like. Is connected. The analysis device 11 includes a calculation unit 21, a comparison unit 22, and a determination unit 23.

計測部12は、製造されてガスタービンに装着される前(変形前)や変形処理前のタービン動翼(以下、動翼)31(図5参照)、使用後や変形処理後の動翼33を3次元(3D)計測するものである。解析装置11は、所定の非弾性ひずみを設定して順解析(FEM解析)を行うことで動翼31の解析モデルを形成し、動翼31の解析モデルと変形後の動翼31とを比較し、逆解析を行うことで動翼31における非弾性ひずみの発生量を予測するものである。表示部13は、解析装置11が推定した動翼31における非弾性ひずみの発生量を表示し、記憶部14は、解析装置11が逐次演算して求めたデータを記憶する。   The measuring section 12 is a turbine blade (hereinafter referred to as a moving blade) 31 (see FIG. 5) before being manufactured and mounted on a gas turbine (before deformation) or before deformation processing (see FIG. 5), and a moving blade 33 after use or after deformation processing. Is measured three-dimensionally (3D). The analysis device 11 sets a predetermined inelastic strain and performs forward analysis (FEM analysis) to form an analysis model of the moving blade 31, and compares the analysis model of the moving blade 31 with the deformed moving blade 31. Then, the amount of inelastic strain generated in the moving blade 31 is predicted by performing inverse analysis. The display unit 13 displays the amount of inelastic strain generated in the moving blade 31 estimated by the analysis device 11, and the storage unit 14 stores data obtained by the sequential calculation by the analysis device 11.

また、演算部21は、計測部12が計測した変形前の動翼31の表面及び内部が分割された複数のメッシュを所定の大きさの単位領域に区画し、解析モデルを用いて単位領域の各メッシュに対して一様な非弾性ひずみを入力して動翼31の解析モデル形状を求めるものである。比較部22は、動翼31の解析モデル形状と変形後の動翼31の形状とを比較するものである。判定部23は、比較部22の比較結果が収束条件を満たさないときに演算部21の処理を繰り返す一方、比較部22の比較結果が収束条件を満たすきに解析モデルに入力する非弾性ひずみを確定するものである。   In addition, the calculation unit 21 divides a plurality of meshes obtained by dividing the surface and the inside of the moving blade 31 before deformation measured by the measurement unit 12 into unit regions of a predetermined size, and uses the analysis model to determine the unit regions. A uniform inelastic strain is input to each mesh to obtain an analysis model shape of the moving blade 31. The comparison unit 22 compares the analysis model shape of the moving blade 31 with the deformed shape of the moving blade 31. The determination unit 23 repeats the processing of the calculation unit 21 when the comparison result of the comparison unit 22 does not satisfy the convergence condition, while the inelastic strain input to the analysis model when the comparison result of the comparison unit 22 satisfies the convergence condition. It will be confirmed.

まず、変形後のタービン動翼の評価方法を説明する。図2は、使用後のタービン動翼の評価方法を表すフローチャート、図3は、変形処理後のタービン動翼の評価方法を表すフローチャートである。   First, a method for evaluating a turbine blade after deformation will be described. FIG. 2 is a flowchart showing a turbine blade evaluation method after use, and FIG. 3 is a flowchart showing a turbine blade evaluation method after deformation processing.

図2に示すように、ステップS11にて、動翼31を製造し、ステップS12にて、計測部12により変形前の動翼31の3次元形状を計測する。ステップS13にて、動翼31をガスタービンに装着し、このガスタービンを運転することで使用する。ステップS14にて、ガスタービンの運転を停止し、動翼31をガスタービンから取外し、計測部12により変形後の動翼31の3次元形状を計測する。ステップS15にて、動翼31における非弾性ひずみのFEM解析(有限要素解析)を実施し、ステップS16にて、動翼31のクリープひずみ量(ひずみの発生量)を推定する。そして、ステップS17にて、推定したクリープひずみ量に基づいて動翼31を評価する。この動翼31の評価とは、推定したクリープひずみ量と予め設定された規定値とを比較することで、動翼31の継続した使用を判定するものである。   As shown in FIG. 2, the moving blade 31 is manufactured in step S11, and the three-dimensional shape of the moving blade 31 before deformation is measured by the measuring unit 12 in step S12. In step S13, the rotor blade 31 is mounted on a gas turbine and used by operating the gas turbine. In step S14, the operation of the gas turbine is stopped, the moving blade 31 is removed from the gas turbine, and the three-dimensional shape of the deformed moving blade 31 is measured by the measuring unit 12. In step S15, FEM analysis (finite element analysis) of inelastic strain in the moving blade 31 is performed, and in step S16, the creep strain amount (strain generation amount) of the moving blade 31 is estimated. In step S17, the rotor blade 31 is evaluated based on the estimated creep strain amount. The evaluation of the moving blade 31 is to determine the continued use of the moving blade 31 by comparing the estimated creep strain amount with a preset specified value.

また、図3に示すように、ステップS21にて、動翼31を製造(製造中)し、ステップS22にて、計測部12により変形前の動翼31の3次元形状を計測する。なお、ここで、計測部12による変形前の動翼31の3次元形状の計測が実施できなかったとき、設計モデルを代わりに使用してもよい。ステップS23にて、動翼31に対して変形処理を実施する。この変形処理とは、上述した形状修正のために行う曲げ加工や塑性変形による加工硬化や残留応力の付与を図るために行うショットピーニングなどである。ステップS24にて、変形処理を行った変形後の動翼31を計測部12により3次元形状を計測する。ステップS25にて、動翼31における非弾性ひずみのFEM解析(有限要素解析)を実施し、ステップS26にて、動翼31の塑性ひずみ量(塑性ひずみの発生量)を推定する。そして、ステップS27にて、推定した塑性ひずみ量に基づいて動翼31を評価する。   As shown in FIG. 3, the moving blade 31 is manufactured (during manufacture) in step S <b> 21, and the three-dimensional shape of the moving blade 31 before deformation is measured by the measuring unit 12 in step S <b> 22. Here, when the measurement unit 12 cannot measure the three-dimensional shape of the moving blade 31 before deformation, the design model may be used instead. In step S23, deformation processing is performed on the moving blade 31. This deformation process includes a bending process performed for the above-described shape correction, work hardening by plastic deformation, and shot peening performed for imparting residual stress. In step S <b> 24, the measurement unit 12 measures the three-dimensional shape of the deformed moving blade 31 that has undergone the deformation process. In step S25, FEM analysis (finite element analysis) of inelastic strain in the moving blade 31 is performed, and in step S26, the plastic strain amount (the amount of plastic strain generated) of the moving blade 31 is estimated. In step S27, the rotor blade 31 is evaluated based on the estimated amount of plastic strain.

次に、タービン動翼のひずみ推定方法を説明する。具体的には、図2及び図3にて、ステップS12,S22〜ステップS16,S26の処理の説明である。図4は、第1実施形態のタービン動翼のひずみ推定方法を表すフローチャート、図5は、タービン動翼のひずみ推定方法を説明するための概略図である。   Next, a method for estimating the distortion of the turbine blade will be described. Specifically, in FIG. 2 and FIG. 3, the processing of steps S12 and S22 to steps S16 and S26 will be described. FIG. 4 is a flowchart illustrating a turbine blade distortion estimation method according to the first embodiment, and FIG. 5 is a schematic diagram for explaining the turbine blade distortion estimation method.

図4及び図5に示すように、ステップS31にて、計測部12は、変形前の動翼31の3次元形状を計測し、ステップS32にて、演算部21は、動翼31の解析モデルを作成する。そして、演算部21は、ステップS33にて、変形前の動翼31の表面及び内部が分割された複数のメッシュを所定の大きさの単位領域に区画し、ステップS34にて、解析モデルを用いて単位領域の複数のメッシュに対して一様な非弾性ひずみを入力する。   As shown in FIGS. 4 and 5, in step S <b> 31, the measurement unit 12 measures the three-dimensional shape of the moving blade 31 before deformation, and in step S <b> 32, the calculation unit 21 analyzes the analysis model of the moving blade 31. Create Then, in step S33, the calculation unit 21 partitions a plurality of meshes obtained by dividing the surface and the interior of the moving blade 31 before deformation into unit regions of a predetermined size, and uses the analysis model in step S34. The uniform inelastic strain is input to a plurality of meshes in the unit area.

即ち、変形前の動翼31の表面及び内部を複数のメッシュで分割することで空間を離散化し、複数のメッシュで形成される空間を所定の大きさの単位領域として区画する。例えば、領域Sを拡大したもので説明すると、複数のメッシュを所定の個数で区画し、単位領域A1,A2,A3・・・を設定する。この場合、予め実験や順解析などにより求めた非弾性ひずみの発生量が大きい領域では、1つの単位領域の面積を小さくして単位領域を細かく分割することが望ましい。また、予め実験や順解析などにより求めた亀裂などが発生しやすい領域では、1つの単位領域の面積を小さくして単位領域を細かく分割することが望ましい。更に、予め実験や順解析などにより求めた非弾性ひずみの発生量の差(勾配)が大きい領域では、1つの単位領域の面積を小さくして単位領域を細かく分割することが望ましい。   That is, the space and the inside of the moving blade 31 before deformation are divided by a plurality of meshes to discretize the space, and the space formed by the plurality of meshes is partitioned as a unit region of a predetermined size. For example, when the region S is enlarged, a plurality of meshes are divided by a predetermined number, and unit regions A1, A2, A3,... Are set. In this case, in a region where the amount of inelastic strain generated by experiments or forward analysis is large, it is desirable that the unit region is divided finely by reducing the area of one unit region. Further, in a region where cracks or the like that are obtained in advance by experiments or forward analysis are likely to occur, it is desirable to reduce the area of one unit region and to divide the unit region finely. Furthermore, in a region where the difference (gradient) in the generation amount of inelastic strain obtained in advance by experiments or forward analysis is large, it is desirable to make the unit region finely divided by reducing the area of one unit region.

動翼31の運転中のクリープひずみの逆解析を行う場合は、翼面表面及び内部の全体に単位領域を設ける必要があるが、一方、歪修正の塑性ひずみの逆解析を行う場合は、通常プレス冶具の接する周辺のみ細かく分割すればよい、また、ショットピーニングは極表面にのみ塑性変形が発生するので、単位領域も極表面にのみ設定すればよい。   When performing reverse analysis of creep strain during operation of the rotor blade 31, it is necessary to provide a unit region on the entire surface and inside of the blade surface. On the other hand, when performing reverse analysis of plastic strain for strain correction, It is only necessary to finely divide only the periphery where the press jig comes into contact. Also, since shot peening causes plastic deformation only on the pole surface, the unit region may be set only on the pole surface.

そして、複数の単位領域の各メッシュに対して一様な非弾性ひずみを入力する。ここで、一様な非弾性ひずみとは、所定の条件に対応して設定された大きさが同一の非弾性ひずみである。また、一様な非弾性ひずみとは、単位領域の一方の角部から他方の角部(コントロールポイント)に向けて線形(直性状)に変化(増加または減少)する非弾性ひずみである。この場合、非弾性ひずみの大きさは、逆解析により求める。   Then, uniform inelastic strain is input to each mesh of a plurality of unit regions. Here, the uniform inelastic strain is an inelastic strain having the same magnitude set corresponding to a predetermined condition. The uniform inelastic strain is an inelastic strain that changes (increases or decreases) linearly (directly) from one corner of the unit region to the other corner (control point). In this case, the magnitude of inelastic strain is obtained by inverse analysis.

また、複数のメッシュで形成される空間を所定の大きさの単位領域として区画することで、動翼31の表面及び内部に対して複数の単位領域が設定され、この複数の単位領域ごとに異なる大きさの非弾性ひずみを入力することとなる。そして、ステップS34にて、解析モデルを用いて複数の単位領域における各メッシュに所定の大きさの非弾性ひずみが入力されると、演算部21は、動翼31の解析モデル形状32を求める。   Further, by dividing a space formed by a plurality of meshes as unit areas of a predetermined size, a plurality of unit areas are set for the surface and the inside of the moving blade 31, and each unit area is different. The magnitude of inelastic strain will be input. In step S <b> 34, when an inelastic strain having a predetermined size is input to each mesh in the plurality of unit regions using the analysis model, the calculation unit 21 obtains the analysis model shape 32 of the moving blade 31.

ステップS35にて、計測部12は、変形後の動翼33の3次元形状を計測し、ステップS36にて、比較部22は、動翼31の解析モデル形状32と変形後の動翼33の形状とを比較する。そして、ステップS37にて、判定部23は、動翼31の解析モデル形状32と変形後の動翼33の形状との比較結果が収束条件を満たすかどうかを判定する。   In step S35, the measurement unit 12 measures the three-dimensional shape of the deformed moving blade 33. In step S36, the comparison unit 22 determines the analysis model shape 32 of the moving blade 31 and the deformed moving blade 33. Compare the shape. In step S <b> 37, the determination unit 23 determines whether the comparison result between the analysis model shape 32 of the moving blade 31 and the shape of the deformed moving blade 33 satisfies the convergence condition.

ここで、収束条件とは、動翼31の解析モデル形状32と変形後の動翼33の形状との変形差が予め設定された規定値以内にあることである。具体的に、この変形差とは、動翼31の解析モデル形状32の体積形状と、3次元計測した変形後の動翼33の体積形状とのずれ体積量である。なお、ずれ体積量は、ずれ長さの自乗和でもよい。動翼31の解析モデル形状32の3次元形状と変形後の動翼33の3次元形状とを重ね合わせたとき、両者が重ならない体積である。   Here, the convergence condition is that the deformation difference between the analysis model shape 32 of the moving blade 31 and the shape of the deformed moving blade 33 is within a predetermined value set in advance. Specifically, the deformation difference is a deviation volume amount between the volume shape of the analysis model shape 32 of the moving blade 31 and the volume shape of the deformed moving blade 33 measured in three dimensions. The displacement volume may be a square sum of the displacement lengths. When the three-dimensional shape of the analysis model shape 32 of the moving blade 31 and the three-dimensional shape of the deformed moving blade 33 are superimposed, the volume does not overlap.

ステップS37にて、動翼31の解析モデル形状32と変形後の動翼33の形状との変形差が規定値以内にないと判定(No)されると、ステップS34に戻り、解析モデルに入力する非弾性ひずみの大きさを変更(増加または減少)し、ステップS5,S36,S37の処理を繰り返す。一方、ステップS37にて、動翼31の解析モデル形状32と変形後の動翼33の形状との変形差が規定値以内であると判定(Yes)されると、ステップS38にて、解析モデルに入力した非弾性ひずみの大きさを確定し、推定値として出力する。   If it is determined in step S37 that the deformation difference between the analysis model shape 32 of the moving blade 31 and the shape of the deformed moving blade 33 is not within the specified value (No), the process returns to step S34 and input to the analysis model. The magnitude of the inelastic strain to be changed is changed (increased or decreased), and the processes of steps S5, S36, and S37 are repeated. On the other hand, if it is determined in step S37 that the deformation difference between the analysis model shape 32 of the moving blade 31 and the shape of the deformed moving blade 33 is within the specified value (Yes), the analysis model is determined in step S38. The magnitude of the inelastic strain input to is determined and output as an estimated value.

このように第1実施形態のタービンコンポーネントのひずみ推定方法にあっては、変形前の動翼31を3次元計測するステップと、動翼31の解析モデルを作成するステップと、変形前の動翼31の表面及び内部が分割された複数のメッシュを所定の大きさの単位領域に区画するステップと、解析モデルを用いて単位領域の複数のメッシュに対して一様な非弾性ひずみを入力して動翼31の解析モデル形状32を形成するステップと、変形後の動翼31を3次元計測するステップと、動翼31の解析モデル形状32と変形後の動翼33の形状とを比較するステップと、動翼31の解析モデル形状32と変形後の動翼33の形状との比較結果が収束条件を満たさないときに解析モデルに入力する非弾性ひずみを変更するステップと、動翼31の解析モデル形状32と変形後の動翼33の形状との比較結果が収束条件を満たすときに解析モデルに入力した非弾性ひずみを確定するステップとを有する。   As described above, in the method for estimating the distortion of the turbine component according to the first embodiment, the step of three-dimensionally measuring the moving blade 31 before deformation, the step of creating an analysis model of the moving blade 31, and the moving blade before deformation. A step of dividing a plurality of meshes of which the surface and the interior of 31 are divided into unit areas of a predetermined size, and inputting a uniform inelastic strain to the plurality of meshes of the unit areas using an analysis model The step of forming the analysis model shape 32 of the moving blade 31, the step of three-dimensionally measuring the deformed moving blade 31, and the step of comparing the analysis model shape 32 of the moving blade 31 and the shape of the deformed moving blade 33. Changing the inelastic strain input to the analysis model when the comparison result between the analysis model shape 32 of the moving blade 31 and the shape of the deformed moving blade 33 does not satisfy the convergence condition; And a step of comparison between the shape of the blade 33 after deformation analysis model shape 32 is fixed strain inelastic entered into the analysis model when the convergence condition is satisfied.

従って、単位領域に一様な非弾性ひずみを入力することで、逆解析時における変数(パラメータ)の数が減少して発散を防止することができると共に、自由に非弾性ひずみを入力して解析を行うことで動翼31に発生した非弾性ひずみを精度良く推定することができる。   Therefore, by inputting a uniform inelastic strain in the unit area, the number of variables (parameters) at the time of inverse analysis can be reduced and divergence can be prevented. By performing the above, it is possible to accurately estimate the inelastic strain generated in the moving blade 31.

第1実施形態のタービンコンポーネントのひずみ推定方法では、単位領域を動翼31の表面及び内部に対して複数設定し、複数の単位領域ごとに異なる非弾性ひずみを入力する。従って、動翼31の変形に応じた適正な解析モデル形状32を形成することができる。   In the turbine component strain estimation method of the first embodiment, a plurality of unit regions are set with respect to the surface and the inside of the moving blade 31, and different inelastic strains are input for each of the plurality of unit regions. Therefore, an appropriate analysis model shape 32 corresponding to the deformation of the moving blade 31 can be formed.

第1実施形態のタービンコンポーネントのひずみ推定方法では、解析モデルを用いて入力する一様な非弾性ひずみとして、同一の非弾性ひずみを設定している。従って、変数(パラメータ)の数が減少して発散を防止することができる。   In the turbine component strain estimation method of the first embodiment, the same inelastic strain is set as the uniform inelastic strain input using the analysis model. Therefore, the number of variables (parameters) can be reduced and divergence can be prevented.

第1実施形態のタービンコンポーネントのひずみ推定方法では、解析モデルを用いて入力する一様な非弾性ひずみとして、線形に変化する非弾性ひずみを設定している。従って、変数(パラメータ)の数が減少して発散を防止することができる。   In the turbine component strain estimation method of the first embodiment, linearly changing inelastic strain is set as the uniform inelastic strain input using the analysis model. Therefore, the number of variables (parameters) can be reduced and divergence can be prevented.

第1実施形態のタービンコンポーネントのひずみ推定方法では、収束条件として、動翼31の解析モデル形状32と変形後の動翼33の形状との変形差が予め設定された規定値以内にあることを設定している。従って、解析モデル形状32と変形後の動翼33の形状との変形差が規定値以内にあるときに、収束条件を満たして解析モデルに入力した非弾性ひずみを確定することで、動翼31にマークを付けたり、加工を施す必要がなく、容易に変形差を比較することができる。   In the turbine component strain estimation method according to the first embodiment, as a convergence condition, the deformation difference between the analysis model shape 32 of the moving blade 31 and the shape of the deformed moving blade 33 is within a preset specified value. It is set. Accordingly, when the deformation difference between the analysis model shape 32 and the shape of the deformed moving blade 33 is within a specified value, the inelastic strain that satisfies the convergence condition and is input to the analysis model is determined, thereby moving the moving blade 31. It is not necessary to mark or process the pattern, and the difference in deformation can be easily compared.

第1実施形態のタービンコンポーネントのひずみ推定方法では、変形差として、動翼31の解析モデルの体積形状と変形後の動翼33の体積形状とのずれ体積量を設定している。従って、高精度に変形差を求めることができる。   In the turbine component strain estimation method of the first embodiment, the displacement volume between the volume shape of the analysis model of the moving blade 31 and the volume shape of the deformed moving blade 33 is set as the deformation difference. Therefore, the deformation difference can be obtained with high accuracy.

また、第1実施形態のタービンコンポーネントのひずみ推定装置にあっては、動翼31を3次元計測する計測部12と、計測部12が計測した変形前の動翼31の表面及び内部が分割された複数のメッシュを所定の大きさの単位領域に区画し、解析モデルを用いて単位領域の複数のメッシュに対して一様な非弾性ひずみを入力して動翼31の解析モデル形状32を求める演算部21と、動翼31の解析モデル形状32と計測部12が計測した変形後の動翼33の形状とを比較する比較部22と、比較部22の比較結果が収束条件を満たさないときに演算部21の処理を繰り返す一方、比較部22の比較結果が収束条件を満たすときに解析モデルに入力する非弾性ひずみを確定する判定部23とを備える。   In the turbine component strain estimation apparatus of the first embodiment, the measurement unit 12 that three-dimensionally measures the moving blade 31 and the surface and the inside of the moving blade 31 before deformation measured by the measuring unit 12 are divided. A plurality of meshes are divided into unit regions of a predetermined size, and uniform inelastic strain is input to the plurality of meshes in the unit region using the analysis model to obtain the analysis model shape 32 of the moving blade 31. When the comparison result of the calculation unit 21, the analysis model shape 32 of the moving blade 31 and the shape of the deformed moving blade 33 measured by the measuring unit 12 does not satisfy the convergence condition And a determination unit 23 that determines the inelastic strain to be input to the analysis model when the comparison result of the comparison unit 22 satisfies the convergence condition.

従って、単位領域に一様な非弾性ひずみを入力することで、逆解析時における変数(パラメータ)の数が減少して発散を防止することができると共に、自由に非弾性ひずみを入力して解析を行うことで動翼31に発生した非弾性ひずみを精度良く推定することができる。   Therefore, by inputting a uniform inelastic strain in the unit area, the number of variables (parameters) at the time of inverse analysis can be reduced and divergence can be prevented. By performing the above, it is possible to accurately estimate the inelastic strain generated in the moving blade 31.

また、第1実施形態のタービンコンポーネントの評価方法にあっては、タービンコンポーネントのひずみ推定方法により推定した非弾性ひずみを用いて動翼31を評価する。従って、逆解析時における変数(パラメータ)の数が減少して発散を防止することができると共に、自由に非弾性ひずみを入力して解析を行うことで動翼31に発生した非弾性ひずみを精度良く推定することができ、その結果、動翼31の評価を精度良く実施することができる。   In the turbine component evaluation method of the first embodiment, the moving blade 31 is evaluated using the inelastic strain estimated by the turbine component strain estimation method. Accordingly, the number of variables (parameters) at the time of inverse analysis can be reduced to prevent divergence, and the inelastic strain generated in the moving blade 31 can be accurately analyzed by freely inputting inelastic strain. As a result, the rotor blade 31 can be evaluated with high accuracy.

[第2実施形態]
図6は、第2実施形態のタービン動翼のひずみ推定方法を表すフローチャートである。なお、本実施形態の基本的な構成は、上述した第1実施形態とほぼ同様の構成であり、図1を用いて説明すると共に、上述した第1実施形態と同様の機能を有する部材には、同一の符号を付して詳細な説明は省略する。
[Second Embodiment]
FIG. 6 is a flowchart showing a turbine blade distortion estimation method according to the second embodiment. The basic configuration of the present embodiment is substantially the same as that of the first embodiment described above, and a member having the same function as that of the first embodiment described above will be described with reference to FIG. The same reference numerals are attached and detailed description is omitted.

第2実施形態において、図1に示すように、第2実施形態のタービン動翼の非弾性ひずみ推定装置10は、解析装置11に計測部12と表示部13と記憶部14が接続されて構成され、解析装置11は、演算部21と、比較部22と、判定部23とを備えている。   In the second embodiment, as shown in FIG. 1, the turbine blade inelastic strain estimation device 10 of the second embodiment is configured by connecting a measurement unit 12, a display unit 13, and a storage unit 14 to an analysis device 11. The analysis device 11 includes a calculation unit 21, a comparison unit 22, and a determination unit 23.

計測部12は、使用前(変形前)の動翼31や使用後(変形後)の動翼33を3次元計測するものである。解析装置11は、所定の非弾性ひずみを設定して順解析(FEM解析)を行うことで動翼31の複数の解析モデル形状を形成し、動翼31の解析モデルと変形後の動翼31とを比較することで動翼31における非弾性ひずみの発生量を予測するものである。表示部13は、解析装置11が推定した動翼31における非弾性ひずみの発生量を表示し、記憶部14は、解析装置11が演算して求めた複数の解析モデル形状を記憶する。   The measuring unit 12 performs three-dimensional measurement of the moving blade 31 before use (before deformation) and the moving blade 33 after use (after deformation). The analysis device 11 sets a predetermined inelastic strain and performs forward analysis (FEM analysis) to form a plurality of analysis model shapes of the moving blade 31, and the analysis model of the moving blade 31 and the deformed moving blade 31. Is used to predict the amount of inelastic strain generated in the rotor blade 31. The display unit 13 displays the amount of inelastic strain generated in the moving blade 31 estimated by the analysis device 11, and the storage unit 14 stores a plurality of analysis model shapes calculated by the analysis device 11.

ここで、タービン動翼のひずみ推定方法を説明する。図6に示すように、ステップS41にて、計測部12は、変形前の動翼31の3次元形状を計測し、ステップS42にて、演算部21は、動翼31の解析モデルを作成する。そして、ステップS43にて、演算部21は、解析モデルに複数種類の変形条件(変形モード)を入力して解析することで、動翼31の複数の解析モデル形状を形成し、記憶部14に記憶してデータベースを作成する。ここで、複数種類の変形条件(変形モード)とは、例えば、動翼31に作用する応力、温度などである。   Here, a method for estimating the distortion of the turbine rotor blade will be described. As shown in FIG. 6, in step S41, the measurement unit 12 measures the three-dimensional shape of the moving blade 31 before deformation, and in step S42, the calculation unit 21 creates an analysis model of the moving blade 31. . In step S43, the calculation unit 21 inputs and analyzes a plurality of types of deformation conditions (deformation modes) in the analysis model, thereby forming a plurality of analysis model shapes of the moving blade 31 in the storage unit 14. Remember and create a database. Here, the multiple types of deformation conditions (deformation modes) include, for example, stress and temperature acting on the moving blade 31.

ステップS44にて、計測部12は、変形後の動翼33の3次元形状を計測し、ステップS45にて、比較部22は、記憶部14のデータベースに記憶された動翼31の複数の解析モデル形状と変形後の動翼31の形状とを比較する。そして、ステップS46にて、判定部23は、複数の解析モデル形状と変形後の動翼31の形状との変形差が最小となる解析モデル形状に入力した非弾性ひずみの大きさを確定し、推定値として出力する。   In step S44, the measurement unit 12 measures the three-dimensional shape of the deformed moving blade 33. In step S45, the comparison unit 22 performs a plurality of analyzes of the moving blade 31 stored in the database of the storage unit 14. The model shape is compared with the shape of the rotor blade 31 after deformation. In step S46, the determination unit 23 determines the magnitude of the inelastic strain input to the analysis model shape that minimizes the deformation difference between the plurality of analysis model shapes and the shape of the deformed rotor blade 31. Output as an estimated value.

このように第2実施形態のタービンコンポーネントのひずみ推定方法にあっては、変形前の動翼31を3次元計測するステップと、動翼31の解析モデルを作成するステップと、解析モデルに複数種類の変形条件を入力して動翼31の複数の解析モデル形状を形成するステップと、変形後の動翼31を3次元計測するステップと、動翼31の複数の解析モデル形状と変形後の動翼31の形状とを比較するステップと、動翼31の複数の解析モデル形状と変形後の動翼31の形状との変形差が最小となる解析モデル形状に入力した非弾性ひずみを確定するステップとを有する。   As described above, in the turbine component strain estimation method according to the second embodiment, the step of three-dimensionally measuring the blade 31 before deformation, the step of creating an analysis model of the blade 31, and a plurality of types of analysis models are included. The step of forming a plurality of analysis model shapes of the moving blade 31 by inputting the deformation conditions of the above, a step of three-dimensionally measuring the moving blade 31 after deformation, the plurality of analysis model shapes of the moving blade 31 and the movement after deformation The step of comparing the shape of the blade 31 and the step of determining the inelastic strain input to the analysis model shape that minimizes the deformation difference between the plurality of analysis model shapes of the blade 31 and the shape of the blade 31 after deformation. And have.

従って、解析モデルに複数種類の変形条件を入力して動翼31の複数の解析モデル形状を形成し、変形後の動翼31の形状と比較する。ここで、動翼31の複数の解析モデル形状と変形後の動翼31の形状との変形差が最小となる非弾性ひずみを確定して出力する。そのため、逆解析時における変数(パラメータ)の数が減少して発散を防止することができると共に、自由に非弾性ひずみを入力して解析を行うことで動翼31に発生した非弾性ひずみを精度良く推定することができる。   Accordingly, a plurality of types of deformation conditions are input to the analysis model to form a plurality of analysis model shapes of the moving blade 31 and compared with the shape of the moving blade 31 after deformation. Here, the inelastic strain that minimizes the deformation difference between the plurality of analysis model shapes of the moving blade 31 and the shape of the deformed moving blade 31 is determined and output. Therefore, the number of variables (parameters) at the time of inverse analysis can be reduced and divergence can be prevented, and inelastic strain generated in the moving blade 31 can be accurately analyzed by inputting inelastic strain freely and performing analysis. It can be estimated well.

なお、上述した実施形態では、タービンコンポーネントとしてタービン動翼31を用いて説明したが、タービン動翼31に限らず、圧縮機の動翼や他の回転体、その他の部材でもよく、特に、3次元形状をなす部材に有効的である。   In the above-described embodiment, the turbine blade 31 is used as a turbine component. However, the turbine blade 31 is not limited to the turbine blade 31, and may be a compressor blade, another rotating body, or other member. It is effective for members having a dimensional shape.

10 タービン動翼のひずみ推定装置
11 解析装置
12 計測部
13 表示部
14 記憶部
21 演算部
22 比較部
23 判定部
31,33 タービン動翼(タービンコンポーネント)
32 解析モデル形状
DESCRIPTION OF SYMBOLS 10 Turbine blade distortion estimation device 11 Analysis device 12 Measurement unit 13 Display unit 14 Storage unit 21 Calculation unit 22 Comparison unit 23 Judgment unit 31, 33 Turbine blade (turbine component)
32 Analysis model shape

Claims (9)

変形前のタービンコンポーネントを3次元計測するステップと、
前記タービンコンポーネントの解析モデルを作成するステップと、
前記変形前のタービンコンポーネントが分割された複数のメッシュを所定の大きさの単位領域に区画するステップと、
前記解析モデルを用いて前記単位領域の前記複数のメッシュに対して一様な非弾性ひずみを入力して前記タービンコンポーネントの解析モデル形状を形成するステップと、
変形後のタービンコンポーネントを3次元計測するステップと、
前記タービンコンポーネントの解析モデル形状と変形後の前記タービンコンポーネントの形状とを比較するステップと、
前記タービンコンポーネントの解析モデル形状と変形後の前記タービンコンポーネントの形状との比較結果が収束条件を満たさないときに前記解析モデルに入力する非弾性ひずみを変更して再度前記タービンコンポーネントの解析モデル形状を形成するステップと、
前記タービンコンポーネントの解析モデル形状と変形後の前記タービンコンポーネントの形状との比較結果が収束条件を満たすときに前記解析モデルに入力した非弾性ひずみを確定するステップと、
を有することを特徴とするタービンコンポーネントのひずみ推定方法。
Three-dimensional measurement of turbine components before deformation;
Creating an analytical model of the turbine component;
Partitioning a plurality of meshes obtained by dividing the turbine component before deformation into unit regions of a predetermined size;
Inputting uniform inelastic strain to the plurality of meshes in the unit region using the analysis model to form an analysis model shape of the turbine component;
Three-dimensional measurement of the turbine component after deformation;
Comparing the analytical model shape of the turbine component with the deformed shape of the turbine component;
When the result of comparison between the analytical model shape of the turbine component and the deformed turbine component shape does not satisfy the convergence condition, the inelastic strain input to the analytical model is changed and the analytical model shape of the turbine component is changed again. Forming step;
Determining the inelastic strain input to the analytical model when a comparison result between the analytical model shape of the turbine component and the shape of the turbine component after deformation satisfies a convergence condition;
A strain estimation method for a turbine component, comprising:
前記単位領域は、前記タービンコンポーネントに対して複数設定し、複数の前記単位領域ごとに異なる非弾性ひずみを入力することを特徴とする請求項1に記載のタービンコンポーネントのひずみ推定方法。   2. The turbine component strain estimation method according to claim 1, wherein a plurality of unit regions are set for the turbine component, and different inelastic strains are input to the plurality of unit regions. 3. 前記解析モデルを用いて入力する一様な非弾性ひずみとは、同一の非弾性ひずみであることを特徴とする請求項1または請求項2に記載のタービンコンポーネントのひずみ推定方法。   The turbine component strain estimation method according to claim 1, wherein the uniform inelastic strain input using the analysis model is the same inelastic strain. 前記解析モデルを用いて入力する一様な非弾性ひずみとは、線形に変化する非弾性ひずみであることを特徴とする請求項1または請求項2に記載のタービンコンポーネントのひずみ推定方法。   The turbine component strain estimation method according to claim 1, wherein the uniform inelastic strain input using the analysis model is a linearly varying inelastic strain. 前記収束条件は、前記タービンコンポーネントの解析モデル形状と3次元計測した変形後の前記タービンコンポーネントの形状との変形差が予め設定された規定値以内にあることを特徴とする請求項1から請求項4のいずれか一項に記載のタービンコンポーネントのひずみ推定方法。   2. The convergence condition according to claim 1, wherein a deformation difference between an analysis model shape of the turbine component and a shape of the turbine component after deformation measured in three dimensions is within a predetermined value set in advance. The turbine component strain estimation method according to claim 4. 前記変形差とは、前記タービンコンポーネントの解析モデルの体積形状と3次元計測した変形後の前記タービンコンポーネントの体積形状とのずれ体積量であることを特徴とする請求項5に記載のタービンコンポーネントのひずみ推定方法。   6. The turbine component according to claim 5, wherein the deformation difference is a deviation volume amount between a volume shape of an analysis model of the turbine component and a volume shape of the turbine component after deformation measured three-dimensionally. Strain estimation method. 変形前のタービンコンポーネントを3次元計測するステップと、
前記タービンコンポーネントの解析モデルを作成するステップと、
前記解析モデルに複数種類の変形条件を入力して前記タービンコンポーネントの複数の解析モデル形状を形成するステップと、
変形後のタービンコンポーネントを3次元計測するステップと、
前記タービンコンポーネントの複数の解析モデル形状と変形後の前記タービンコンポーネントの形状とを比較するステップと、
前記タービンコンポーネントの複数の解析モデル形状と変形後の前記タービンコンポーネントの形状との変形差が最小となる前記解析モデル形状に入力した非弾性ひずみを確定するステップと、
を有することを特徴とするタービンコンポーネントのひずみ推定方法。
Three-dimensional measurement of turbine components before deformation;
Creating an analytical model of the turbine component;
Inputting a plurality of types of deformation conditions into the analysis model to form a plurality of analysis model shapes of the turbine component;
Three-dimensional measurement of the turbine component after deformation;
Comparing a plurality of analytical model shapes of the turbine component with shapes of the turbine component after deformation;
Determining an inelastic strain input to the analytical model shape that minimizes a deformation difference between the plurality of analytical model shapes of the turbine component and the shape of the turbine component after deformation;
A strain estimation method for a turbine component, comprising:
タービンコンポーネントを3次元計測する計測部と、
前記計測部が計測した変形前の前記タービンコンポーネントが分割された複数のメッシュを所定の大きさの単位領域に区画し、解析モデルを用いて前記単位領域の前記複数のメッシュに対して一様な非弾性ひずみを入力して前記タービンコンポーネントの解析モデル形状を求める演算部と、
前記タービンコンポーネントの解析モデル形状と前記計測部が計測した変形後の前記タービンコンポーネントの形状とを比較する比較部と、
前記比較部の比較結果が収束条件を満たさないときに前記演算部の処理を繰り返す一方、前記比較部の比較結果が収束条件を満たすときに解析モデルに入力する非弾性ひずみを確定する判定部と、
を備えることを特徴とするタービンコンポーネントのひずみ推定装置。
A measurement unit for three-dimensional measurement of turbine components;
A plurality of meshes obtained by dividing the turbine component before deformation measured by the measurement unit are partitioned into unit regions of a predetermined size, and uniform with respect to the plurality of meshes in the unit regions using an analysis model A calculation unit for inputting an inelastic strain to obtain an analysis model shape of the turbine component;
A comparison unit for comparing the analysis model shape of the turbine component and the shape of the turbine component after deformation measured by the measurement unit;
A determination unit that determines the inelastic strain to be input to the analysis model when the comparison result of the comparison unit does not satisfy a convergence condition and repeats the processing of the calculation unit, while the comparison result of the comparison unit satisfies a convergence condition; ,
A strain estimation device for a turbine component, comprising:
請求項1から請求項7のいずれか一項に記載のタービンコンポーネントのひずみ推定方法により推定した非弾性ひずみを用いて前記タービンコンポーネントを評価することを特徴とするタービンコンポーネントの評価方法。   A turbine component evaluation method, wherein the turbine component is evaluated using inelastic strain estimated by the turbine component strain estimation method according to any one of claims 1 to 7.
JP2016122021A 2016-06-20 2016-06-20 Turbine component strain estimation method and apparatus, and turbine component evaluation method Active JP6578250B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2016122021A JP6578250B2 (en) 2016-06-20 2016-06-20 Turbine component strain estimation method and apparatus, and turbine component evaluation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016122021A JP6578250B2 (en) 2016-06-20 2016-06-20 Turbine component strain estimation method and apparatus, and turbine component evaluation method

Publications (2)

Publication Number Publication Date
JP2017227137A JP2017227137A (en) 2017-12-28
JP6578250B2 true JP6578250B2 (en) 2019-09-18

Family

ID=60891222

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016122021A Active JP6578250B2 (en) 2016-06-20 2016-06-20 Turbine component strain estimation method and apparatus, and turbine component evaluation method

Country Status (1)

Country Link
JP (1) JP6578250B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109595040A (en) * 2018-12-21 2019-04-09 杭州汽轮动力集团有限公司 A kind of regular design method of gas turbine blades distortion

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0861086A (en) * 1994-08-22 1996-03-05 Toshiba Corp Foreign matter collision damage evaluation method in gas turbine blade
JP3788901B2 (en) * 2000-09-27 2006-06-21 株式会社日立製作所 Damage diagnosis device for power generation facilities
JP4631335B2 (en) * 2004-07-21 2011-02-16 Jfeスチール株式会社 Method for determining the shape of a press mold
JP4716171B2 (en) * 2005-05-16 2011-07-06 有限会社エムアンドエムリサーチ Die design support program and method for press working
JP2010224648A (en) * 2009-03-19 2010-10-07 Ihi Corp Analytic model creation device, analytic model creation method, and analytic model creation program

Also Published As

Publication number Publication date
JP2017227137A (en) 2017-12-28

Similar Documents

Publication Publication Date Title
US10330564B2 (en) System and method for predicting distortion of a workpiece resulting from a peening machine process
US9483605B2 (en) Probabilistic high cycle fatigue (HCF) design optimization process
US7870792B2 (en) Forming limit strain analysis
Booysen et al. Fatigue life assessment of a low pressure steam turbine blade during transient resonant conditions using a probabilistic approach
CN103886125B (en) A kind of titanium alloy hot combined shaping method for numerical simulation
WO2008084920A1 (en) Curved surface forming method of a metal plate
CN108763839B (en) A kind of centrifugal impeller boss crack propagation characteristic simulation part design method
KR100783417B1 (en) Curved surface forming method of a steel plate for a ship using a multi-point press
Santos et al. On the determination of flow stress using bulge test and mechanical measurement
CN111339703A (en) Virtual prediction method for material hardening behavior under large strain condition
JP6339969B2 (en) Deformation resistance identification method for thin inspection materials
JP6578250B2 (en) Turbine component strain estimation method and apparatus, and turbine component evaluation method
JP2005227065A (en) Display method for estimated life duration of structure, method of creating life duration distribution, and image display
Ilg et al. Application of a full-field calibration concept for parameter identification of HS-steel with LS-OPT®
JP6446740B2 (en) Simulation apparatus, simulation method, and program
CN117309341A (en) Cylinder liner fatigue testing device and cylinder liner testing method
JP4931511B2 (en) Evaluation Method of Creep Remaining Life of High Temperature Member by Inverse Analysis Considering Stress Relaxation
Reuter et al. Probabilistic Structure-Mechanical Assessment of Rotor Discs Considering Geometry Variations
JP2017090238A (en) Prediction method
KR101869699B1 (en) Method and apparatus for fan simulation by flow analysis
Abdalla et al. Validation of the Nonlinear Superposition Method (NSM) for elastic shakedown limit pressures via comparison with experimental test results of spherical vessels with radial and oblique nozzles
WO2021149170A1 (en) Crack estimation device, crack estimation method, crack inspection method, and failure diagnosis method
Mashiachidi et al. Fatigue life prediction of mistuned steam turbine blades subjected to deviations in blade geometry
Wróbel et al. Numerical models of polymeric composite to simulate fatigue and ageing processes
dos Reis Constitutive parameters identification of metal sheets using circular bulge tests

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20181002

TRDD Decision of grant or rejection written
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20190724

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20190730

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20190826

R150 Certificate of patent or registration of utility model

Ref document number: 6578250

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150