JP2014070976A - Piping stress evaluation device and method - Google Patents

Piping stress evaluation device and method Download PDF

Info

Publication number
JP2014070976A
JP2014070976A JP2012216401A JP2012216401A JP2014070976A JP 2014070976 A JP2014070976 A JP 2014070976A JP 2012216401 A JP2012216401 A JP 2012216401A JP 2012216401 A JP2012216401 A JP 2012216401A JP 2014070976 A JP2014070976 A JP 2014070976A
Authority
JP
Japan
Prior art keywords
branch pipe
stress
spectrum
pipe
acceleration
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.)
Granted
Application number
JP2012216401A
Other languages
Japanese (ja)
Other versions
JP5946064B2 (en
Inventor
Yoshitsugu Nekomoto
善続 猫本
Moriyoshi Suzuki
盛喜 鈴木
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 JP2012216401A priority Critical patent/JP5946064B2/en
Publication of JP2014070976A publication Critical patent/JP2014070976A/en
Application granted granted Critical
Publication of JP5946064B2 publication Critical patent/JP5946064B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

PROBLEM TO BE SOLVED: To provide a piping stress evaluation device and a method for calculating stress with high accuracy and performing stress evaluation at an arbitrary evaluation point by using easy measurement.SOLUTION: A piping stress evaluation device has: a vibration frequency specification section 32 which specifies a first vibration frequency of a branch pipe based on a time history waveform of acceleration when arbitrary vibration is given to the branch pipe; an analytic model determination section 33 which calculates a second vibration frequency of the branch pipe based on a first analytic model obtained by modeling a header by rotational rigidity, horizontal rigidity, and concentrated mass, modeling the branch pipe by a beam element, and modeling a connection part between the header and the branch pipe by the rotational rigidity, and determines a second analytic model in which the second vibration frequency matches to the first vibration frequency; a steady response measurement section 34 which calculates a third acceleration spectrum based on the time history waveform of the acceleration in a steady operation of a plant; and a generated stress evaluation section 35 which calculates peak stress based on the second analytic model and the third acceleration spectrum, compares and evaluates allowable stress of the peak stress and the branch pipe.

Description

本発明は、配管応力評価装置及び方法に関する。   The present invention relates to a piping stress evaluation apparatus and method.

配管の応力評価を行う場合、例えば、図21に示すように、母管51に設けた分岐管52が弁53を有しており、分岐管52が固定壁54にサポート55により固定されている構成において、分岐管52の応力評価を行う場合には、従来は、分岐管52の評価対象部位に歪ゲージ56を貼り付け、歪ゲージ56の出力信号を、アンプ57、A/D変換器58を介して、コンピュータ59に取り込み、図22に示すように、歪と応力の換算により、歪の時刻歴波形(図22(a)参照)を応力の時刻歴波形(図22(b)参照)に換算し、発生応力を算出して、疲労評価を実施していた。   When performing stress evaluation of piping, for example, as shown in FIG. 21, a branch pipe 52 provided in a mother pipe 51 has a valve 53, and the branch pipe 52 is fixed to a fixed wall 54 by a support 55. In the configuration, when stress evaluation of the branch pipe 52 is performed, conventionally, a strain gauge 56 is attached to the evaluation target portion of the branch pipe 52, and an output signal of the strain gauge 56 is supplied to an amplifier 57 and an A / D converter 58. As shown in FIG. 22, the strain time history waveform (see FIG. 22 (a)) is converted into the stress time history waveform (see FIG. 22 (b)). The fatigue stress was evaluated by calculating the generated stress in terms of

又、図23に示すように、分岐管52に設置した加速度計61、62から変位を算出して、配管諸元(部材諸元)に基づき、変位と応力の関係から発生応力を評価したり、図24に示すように、母管51及び分岐管52の解析モデル70を用い、加速度計測点P1、P2での計測データにより応力を算出したりしていた。   23, the displacement is calculated from accelerometers 61 and 62 installed in the branch pipe 52, and the generated stress is evaluated from the relationship between the displacement and the stress based on the piping specifications (member specifications). As shown in FIG. 24, using an analysis model 70 of the mother pipe 51 and the branch pipe 52, stress is calculated from measurement data at the acceleration measurement points P1 and P2.

特開2001−153719号公報JP 2001-153719 A 特開2002−162298号公報JP 2002-162298 A 特開平11−014782号公報Japanese Patent Laid-Open No. 11-014782

しかしながら、図21、図22で説明した歪ゲージ56を用いる場合には、計測が手間であり、貼り付け箇所のみしか発生応力を把握できなかった。又、図23、図24で説明したように、比較的容易に計測を行うため、加速度計61、62の計測データより応力を算定する手法もあるが、解析モデル70が不十分であれば、精度の高い疲労評価は困難であった。   However, when the strain gauge 56 described with reference to FIGS. 21 and 22 is used, the measurement is troublesome, and the generated stress can be grasped only at the pasted portion. As described with reference to FIGS. 23 and 24, there is a method for calculating stress from the measurement data of the accelerometers 61 and 62 in order to perform measurement relatively easily. However, if the analysis model 70 is insufficient, Highly accurate fatigue evaluation was difficult.

本発明は上記課題に鑑みなされたもので、容易な計測を用いて、任意の評価点において、精度の高い応力を算出して、応力評価を行う配管応力評価装置及び方法を提供することを目的とする。   The present invention has been made in view of the above problems, and an object of the present invention is to provide a piping stress evaluation apparatus and method for performing stress evaluation by calculating high-precision stress at an arbitrary evaluation point using easy measurement. And

上記課題を解決する第1の発明に係る配管応力評価装置は、
プラントの母管から分岐した分岐管に、或いは、前記母管及び前記分岐管に取り付ける複数の加速度計と、
前記加速度計で計測した計測値を用いて、前記分岐管に発生する応力を算出し、評価を行う応力評価手段とを備えた配管応力評価装置において、
前記応力評価手段は、
任意の振動を前記分岐管に与えたときの前記複数の加速度計における加速度の時刻歴波形から、周波数分析を行って、前記母管から遠い方の前記加速度計での第1の加速度スペクトルと前記母管に近い方又は前記母管の前記加速度計での第2の加速度スペクトルとを算出し、前記第1の加速度スペクトルを前記第2の加速度スペクトルで除算して、第1の伝達関数を算出し、前記第1の伝達関数のスペクトルピークから前記分岐管の第1の振動数を特定する振動数特定部と、
前記母管を回転剛性、水平剛性及び集中質量でモデル化し、前記分岐管を梁要素でモデル化し、前記母管と前記分岐管の接続部分を回転剛性でモデル化した第1の解析モデルを作成し、前記第1の解析モデルに基づいて、固有値解析を行って、前記分岐管の第2の振動数を算出し、前記第2の振動数が前記第1の振動数に一致するように、前記母管の回転剛性、水平剛性の設定を変更して、前記第2の振動数が前記第1の振動数に一致する第2の解析モデルを決定する解析モデル決定部と、
前記プラントの定常運転時において、前記分岐管の任意の計測点に配置した前記加速度計における加速度の時刻歴波形から、周波数分析を行って、第3の加速度スペクトルを算出する定常応答計測部と、
前記第2の解析モデルを用い、当該第2の解析モデルの周波数応答解析を行って、前記計測点における第4の加速度スペクトルと前記分岐管の任意の評価点における第1の応力スペクトルとを算出し、前記第1の応力スペクトルを前記第4の加速度スペクトルで除算して、第2の伝達関数を算出し、前記第2の伝達関数に前記第3の加速度スペクトルを乗算して、前記評価点における第2の応力スペクトルを算出し、前記第2の応力スペクトルを周波数積分してピーク応力を算出し、前記ピーク応力と前記分岐管の許容応力とを比較することにより、前記分岐管に発生する応力の評価を行う発生応力評価部とを有することを特徴とする。
A piping stress evaluation apparatus according to the first invention for solving the above-mentioned problems is as follows.
A plurality of accelerometers attached to the branch pipe branched from the mother pipe of the plant, or to the mother pipe and the branch pipe;
In the pipe stress evaluation apparatus provided with stress evaluation means for calculating and evaluating the stress generated in the branch pipe using the measurement value measured by the accelerometer,
The stress evaluation means includes
The frequency analysis is performed from the time history waveforms of acceleration in the plurality of accelerometers when arbitrary vibration is applied to the branch pipe, and the first acceleration spectrum in the accelerometer far from the mother pipe and the first acceleration spectrum A first acceleration function is calculated by calculating a second acceleration spectrum that is closer to the mother tube or the accelerometer of the mother tube, and dividing the first acceleration spectrum by the second acceleration spectrum. A frequency specifying unit for specifying the first frequency of the branch pipe from the spectrum peak of the first transfer function;
The first pipe is modeled with rotational rigidity, horizontal rigidity, and concentrated mass, the branch pipe is modeled with a beam element, and a first analysis model is created in which the connection between the mother pipe and the branch pipe is modeled with rotational rigidity. Then, based on the first analysis model, eigenvalue analysis is performed to calculate a second frequency of the branch pipe, so that the second frequency matches the first frequency, An analysis model determination unit that determines a second analysis model in which the second frequency matches the first frequency by changing the rotational rigidity and horizontal rigidity of the mother pipe;
During steady operation of the plant, from a time history waveform of acceleration in the accelerometer arranged at an arbitrary measurement point of the branch pipe, by performing frequency analysis, a steady response measurement unit that calculates a third acceleration spectrum;
Using the second analysis model, frequency response analysis of the second analysis model is performed to calculate a fourth acceleration spectrum at the measurement point and a first stress spectrum at an arbitrary evaluation point of the branch pipe. And dividing the first stress spectrum by the fourth acceleration spectrum to calculate a second transfer function, multiplying the second transfer function by the third acceleration spectrum, and Is generated in the branch pipe by calculating a second stress spectrum at, calculating a peak stress by frequency integrating the second stress spectrum, and comparing the peak stress with an allowable stress of the branch pipe And a generated stress evaluation unit that performs stress evaluation.

上記課題を解決する第2の発明に係る配管応力評価装置は、
上記第1の発明に記載の配管応力評価装置において、
前記解析モデル決定部は、前記分岐管に弁が設けられている場合には、前記弁を集中質量でモデル化して、前記第1の解析モデルを作成することを特徴とする。
A piping stress evaluation apparatus according to the second invention for solving the above-mentioned problems is as follows.
In the piping stress evaluation apparatus according to the first invention,
When the branch pipe is provided with a valve, the analysis model determination unit models the valve with a concentrated mass to create the first analysis model.

上記課題を解決する第3の発明に係る配管応力評価装置は、
上記第1の発明に記載の配管応力評価装置において、
前記解析モデル決定部は、前記分岐管に弁が設けられている場合には、前記弁のハンドルを含めて、前記弁全体を分布質量でモデル化して、前記第1の解析モデルを作成することを特徴とする。
A piping stress evaluation apparatus according to a third invention for solving the above-mentioned problem is as follows.
In the piping stress evaluation apparatus according to the first invention,
The analysis model determination unit creates the first analysis model by modeling the entire valve including a handle of the valve with a distributed mass when the branch pipe is provided with a valve. It is characterized by.

上記課題を解決する第4の発明に係る配管応力評価方法は、
プラントの母管から分岐した分岐管に、或いは、前記母管及び前記分岐管に取り付ける複数の加速度計を用いて、任意の振動を前記分岐管に与えたときの前記複数の加速度計における加速度の時刻歴波形から、周波数分析を行って、前記母管から遠い方の前記加速度計での第1の加速度スペクトルと前記母管に近い方又は前記母管の前記加速度計での第2の加速度スペクトルとを算出し、前記第1の加速度スペクトルを前記第2の加速度スペクトルで除算して、第1の伝達関数を算出し、前記第1の伝達関数のスペクトルピークから前記分岐管の第1の振動数を特定し、
前記母管を回転剛性、水平剛性及び集中質量でモデル化し、前記分岐管を梁要素でモデル化し、前記母管と前記分岐管の接続部分を回転剛性でモデル化した第1の解析モデルを作成し、前記第1の解析モデルに基づいて、固有値解析を行って、前記分岐管の第2の振動数を算出し、前記第2の振動数が前記第1の振動数に一致するように、前記母管の回転剛性、水平剛性の設定を変更して、前記第2の振動数が前記第1の振動数に一致する第2の解析モデルを決定し、
前記プラントの定常運転時において、前記分岐管の任意の計測点に配置した前記加速度計における加速度の時刻歴波形から、周波数分析を行って、第3の加速度スペクトルを算出し、
前記第2の解析モデルを用い、当該第2の解析モデルの周波数応答解析を行って、前記計測点における第4の加速度スペクトルと前記分岐管の任意の評価点における第1の応力スペクトルとを算出し、前記第1の応力スペクトルを前記第4の加速度スペクトルで除算して、第2の伝達関数を算出し、前記第2の伝達関数に前記第3の加速度スペクトルを乗算して、前記評価点における第2の応力スペクトルを算出し、前記第2の応力スペクトルを周波数積分してピーク応力を算出し、前記ピーク応力と前記分岐管の許容応力とを比較することにより、前記分岐管に発生する応力の評価を行うことを特徴とする。
The pipe stress evaluation method according to the fourth invention for solving the above-mentioned problem is as follows.
The acceleration of the plurality of accelerometers when an arbitrary vibration is applied to the branch pipe, using a branch pipe branched from the main pipe of the plant or a plurality of accelerometers attached to the main pipe and the branch pipe. A frequency analysis is performed from the time history waveform, and a first acceleration spectrum at the accelerometer far from the mother tube and a second acceleration spectrum at the accelerometer near or near the mother tube. The first acceleration spectrum is divided by the second acceleration spectrum to calculate a first transfer function, and the first vibration of the branch pipe is calculated from the spectrum peak of the first transfer function. Identify the number,
The first pipe is modeled with rotational rigidity, horizontal rigidity, and concentrated mass, the branch pipe is modeled with a beam element, and a first analysis model is created in which the connection between the mother pipe and the branch pipe is modeled with rotational rigidity. Then, based on the first analysis model, eigenvalue analysis is performed to calculate a second frequency of the branch pipe, so that the second frequency matches the first frequency, Changing the setting of rotational stiffness and horizontal stiffness of the mother pipe to determine a second analysis model in which the second frequency matches the first frequency;
During steady operation of the plant, from the time history waveform of acceleration in the accelerometer arranged at an arbitrary measurement point of the branch pipe, frequency analysis is performed to calculate a third acceleration spectrum,
Using the second analysis model, frequency response analysis of the second analysis model is performed to calculate a fourth acceleration spectrum at the measurement point and a first stress spectrum at an arbitrary evaluation point of the branch pipe. And dividing the first stress spectrum by the fourth acceleration spectrum to calculate a second transfer function, multiplying the second transfer function by the third acceleration spectrum, and Is generated in the branch pipe by calculating a second stress spectrum at, calculating a peak stress by frequency integrating the second stress spectrum, and comparing the peak stress with an allowable stress of the branch pipe It is characterized by evaluating stress.

上記課題を解決する第5の発明に係る配管応力評価方法は、
上記第4の発明に記載の配管応力評価方法において、
前記分岐管に弁が設けられている場合には、前記弁を集中質量でモデル化して、前記第1の解析モデルを作成することを特徴とする。
The pipe stress evaluation method according to the fifth invention for solving the above-mentioned problem is as follows.
In the piping stress evaluation method according to the fourth invention,
When the branch pipe is provided with a valve, the valve is modeled by a concentrated mass to create the first analysis model.

上記課題を解決する第6の発明に係る配管応力評価方法は、
上記第4の発明に記載の配管応力評価方法において、
前記分岐管に弁が設けられている場合には、前記弁のハンドルを含めて、前記弁全体を分布質量でモデル化して、前記第1の解析モデルを作成することを特徴とする。
The pipe stress evaluation method according to the sixth invention for solving the above-mentioned problems is as follows.
In the piping stress evaluation method according to the fourth invention,
When the branch pipe is provided with a valve, the first analysis model is created by modeling the entire valve including the valve handle with a distributed mass.

本発明によれば、加速度計から取得される加速度データを用いて評価を行うので、現場での計測の手間を軽減することができる。又、精度の良い解析モデルを用いて応力を算出するので、信頼性の高い疲労評価を行うことができる。又、解析モデルを用いているので、分岐管の基部のみならず、任意の評価点においても評価を行うことができる。   According to the present invention, since the evaluation is performed using the acceleration data acquired from the accelerometer, it is possible to reduce the labor of measurement at the site. In addition, since the stress is calculated using an accurate analysis model, a highly reliable fatigue evaluation can be performed. Since the analysis model is used, the evaluation can be performed not only at the base of the branch pipe but also at any evaluation point.

本発明に係る配管応力評価装置を説明する概略図である。It is the schematic explaining the piping stress evaluation apparatus which concerns on this invention. 図1に示した配管応力評価装置のブロック図である。It is a block diagram of the piping stress evaluation apparatus shown in FIG. 本発明に係る配管応力評価方法を説明するフローチャートであり、図2における振動数特定部を説明するものである。It is a flowchart explaining the piping stress evaluation method which concerns on this invention, and demonstrates the frequency specific | specification part in FIG. 本発明に係る配管応力評価方法を説明するフローチャートであり、図2における解析モデル決定部を説明するものである。It is a flowchart explaining the piping stress evaluation method which concerns on this invention, and explains the analysis model determination part in FIG. 本発明に係る配管応力評価方法を説明するフローチャートであり、図2における定常応答計測部を説明するものである。It is a flowchart explaining the piping stress evaluation method based on this invention, and demonstrates the steady response measuring part in FIG. 本発明に係る配管応力評価方法を説明するフローチャートであり、図2における発生応力評価部を説明するものである。It is a flowchart explaining the piping stress evaluation method concerning this invention, and demonstrates the generated stress evaluation part in FIG. タッピング試験で取得された加速度の時間歴波形を示すグラフであり、(a)は分岐管、(b)は母管側の計測データである。It is a graph which shows the time history waveform of the acceleration acquired by the tapping test, (a) is a branch pipe, (b) is measurement data on the mother pipe side. 図7(a)、(b)に示した加速度の時間歴波形を周波数分析したグラフであり、(a)は分岐管、(b)は母管側の加速度スペクトルである。FIGS. 7A and 7B are graphs obtained by frequency analysis of the acceleration time history waveform shown in FIGS. 7A and 7B, where FIG. 7A shows a branch pipe and FIG. 図8(a)、(b)に示した加速度スペクトルから求めた母管と分岐管の伝達関数を示すグラフである。It is a graph which shows the transfer function of the main pipe and the branch pipe calculated | required from the acceleration spectrum shown to Fig.8 (a), (b). 本発明に係る配管応力評価方法で用いる解析モデルを示す図である。It is a figure which shows the analysis model used with the piping stress evaluation method which concerns on this invention. 図10に示した解析モデルにおいて、母管のモデル化を説明する図である。It is a figure explaining modeling of a mother pipe in the analysis model shown in FIG. 図10に示した解析モデルにおいて、母管と分岐管との取り合い部の局部剛性を説明する図である。In the analysis model shown in FIG. 10, it is a figure explaining the local rigidity of the connection part of a main pipe and a branch pipe. 図12で説明した局部剛性を求めるために用いるグラフである。13 is a graph used for obtaining the local rigidity described in FIG. 12. 分岐管における振動数と振動モードを説明する図であり、(a)は1次モード、(b)は2次モードである。It is a figure explaining the frequency and vibration mode in a branch pipe, (a) is a primary mode, (b) is a secondary mode. (a)は、定常運転時における分岐管の加速度の時間歴波形を示すグラフであり、(b)は、(a)に示した加速度の時間歴波形を周波数分析した加速度スペクトルのグラフである。(A) is a graph which shows the time history waveform of the acceleration of a branch pipe at the time of steady operation, (b) is a graph of the acceleration spectrum which frequency-analyzed the time history waveform of the acceleration shown to (a). 解析モデルを用いて算出した計測点での加速度スペクトルと評価点での応力スペクトルを説明する図である。It is a figure explaining the acceleration spectrum in the measurement point calculated using the analysis model, and the stress spectrum in an evaluation point. 図16に示した計測点での加速度スペクトルと評価点での応力スペクトルから求めた応力−加速度の伝達関数を示すグラフである。It is a graph which shows the transfer function of the stress-acceleration calculated | required from the acceleration spectrum in the measurement point shown in FIG. 16, and the stress spectrum in an evaluation point. 図15(b)に示した定常運転時における分岐管の加速度スペクトルと、図17に示した応力−加速度の伝達関数とから求めた応力スペクトルを示すグラフである。It is a graph which shows the stress spectrum calculated | required from the acceleration spectrum of the branch pipe at the time of the steady operation shown in FIG.15 (b), and the transfer function of the stress-acceleration shown in FIG. 本発明に係る配管応力評価方法で用いる解析モデルの他の一例を示す図であり、弁の部分のみを示すモデルである。It is a figure which shows another example of the analysis model used with the piping stress evaluation method which concerns on this invention, and is a model which shows only the part of a valve. 図19に示したモデルを用いた解析モデル全体を示す図である。It is a figure which shows the whole analysis model using the model shown in FIG. 歪ゲージを用いた従来の応力評価を説明する概略図である。It is the schematic explaining the conventional stress evaluation using a strain gauge. (a)は、歪の時刻歴波形を示すグラフであり、(b)は、歪から換算した応力の時刻歴波形を示すグラフである。(A) is a graph showing a time history waveform of strain, and (b) is a graph showing a time history waveform of stress converted from strain. 加速度計を用いた従来の応力評価を説明する概略図である。It is the schematic explaining the conventional stress evaluation using an accelerometer. 配管の解析モデルを示す図である。It is a figure which shows the analysis model of piping.

以下、図1〜図20を参照して、本発明に係る配管応力評価装置及び方法の実施形態を説明する。なお、以降の説明では、プラントの配管の構造において、半径の大きい配管を母管と呼び、当該母管から分岐して設けられ、母管より小さい半径を有する配管を分岐管と呼び、振動応力の影響が大きい分岐管を評価対象として説明を行う。   Hereinafter, with reference to FIGS. 1-20, embodiment of the piping stress evaluation apparatus and method which concern on this invention is described. In the following description, in the plant piping structure, a pipe having a large radius is referred to as a mother pipe, a pipe provided by branching from the mother pipe and having a radius smaller than the mother pipe is referred to as a branch pipe, and vibration stress A branch pipe that has a large influence will be described as an evaluation target.

(実施例1)
最初に、図1を参照して、本実施例の配管応力評価装置を説明する。ここで、図1は、本実施例の配管応力評価装置を説明する概略図である。
Example 1
Initially, with reference to FIG. 1, the piping stress evaluation apparatus of a present Example is demonstrated. Here, FIG. 1 is a schematic diagram illustrating the piping stress evaluation apparatus of the present embodiment.

本実施例の配管応力評価装置20は、実際に配置されている配管において、母管11に設けられた分岐管12に発生する応力を算定して、その疲労評価を行うものであり、その応力の算定、評価のため、分岐管12の複数の箇所に(少なくとも2箇所に)加速度計21、22を取り付けて、加速度を計測している。この加速度計21、22は、タッピングハンマー15を用いた加振により生じた分岐管12の加速度やプラントの定常運転時における分岐管12の加速度などを計測することになる。   The pipe stress evaluation apparatus 20 of the present embodiment calculates the stress generated in the branch pipe 12 provided in the mother pipe 11 in the actually arranged pipe, and performs the fatigue evaluation thereof. For calculation and evaluation, accelerometers 21 and 22 are attached to a plurality of locations (at least two locations) of the branch pipe 12 to measure acceleration. The accelerometers 21 and 22 measure the acceleration of the branch pipe 12 caused by the vibration using the tapping hammer 15, the acceleration of the branch pipe 12 during the steady operation of the plant, and the like.

そして、加速度計21で計測された計測データは、計測データの信号を増幅するアンプ23、計測データの信号をA/D(アナログ/デジタル)変換するA/D変換器25を介して、コンピュータ27に入力されて、記録される。同様に、加速度計22で計測された計測データは、計測データの信号を増幅するアンプ24、計測データの信号をA/D変換するA/D変換器26を介して、コンピュータ27に入力されて、記録される。   The measurement data measured by the accelerometer 21 is supplied to a computer 27 via an amplifier 23 that amplifies the measurement data signal and an A / D converter 25 that performs A / D (analog / digital) conversion of the measurement data signal. To be recorded. Similarly, measurement data measured by the accelerometer 22 is input to a computer 27 via an amplifier 24 that amplifies the measurement data signal and an A / D converter 26 that A / D converts the measurement data signal. Recorded.

なお、分岐管12としては、図1に示すように、弁13を有する構成でもよいし、更に、図21で示したように、分岐管12をサポートで固定壁に固定した構成でもよい。又、弁13などが無い分岐管12のみの構成でもよい。   As shown in FIG. 1, the branch pipe 12 may have a configuration having a valve 13, or may further have a configuration in which the branch pipe 12 is fixed to a fixed wall with a support as shown in FIG. 21. Moreover, the structure of only the branch pipe 12 without the valve 13 etc. may be sufficient.

又、分岐管12に設置した加速度計21、22の少なくとも一方は、母管11又は分岐管12の母管11に近い位置に設置することが望ましく、少なくとも他方は、一方の方より母管11から離れた位置に設置することが望ましい。例えば、本実施例では、母管11と接続している分岐管12の基部に加速度計21を設置し、母管11から離れた位置に加速度計22を設置している。   Moreover, it is desirable to install at least one of the accelerometers 21 and 22 installed in the branch pipe 12 at a position closer to the mother pipe 11 of the mother pipe 11 or the branch pipe 12, and at least the other is the mother pipe 11 from one side. It is desirable to install it at a position away from For example, in this embodiment, the accelerometer 21 is installed at the base of the branch pipe 12 connected to the mother pipe 11, and the accelerometer 22 is installed at a position away from the mother pipe 11.

次に、図2〜図18を参照して、本実施例の配管応力評価装置における配管応力評価方法を説明する。ここで、図2は、図1に示した本実施例の配管応力評価装置を説明するブロック図であり、図3〜図6は、図2に示した各ブロックでの手順を説明するフローチャートである。又、図7〜図18は、図3〜図6に示したフローチャートを説明する図である。   Next, with reference to FIGS. 2-18, the piping stress evaluation method in the piping stress evaluation apparatus of a present Example is demonstrated. Here, FIG. 2 is a block diagram for explaining the pipe stress evaluation apparatus of the present embodiment shown in FIG. 1, and FIGS. 3 to 6 are flowcharts for explaining the procedure in each block shown in FIG. is there. 7 to 18 are diagrams for explaining the flowcharts shown in FIGS.

図2に示すように、本実施例の配管応力評価装置20は、主に、振動計測部31、振動特定部32、応力算定部33、定常応答計測部34及び発生応力評価部35から構成されている。本実施例では、コンピュータ27(応力評価手段)が、振動特定部32、応力算定部33、定常応答計測部34及び発生応力評価部35を有している。   As shown in FIG. 2, the pipe stress evaluation apparatus 20 of the present embodiment mainly includes a vibration measurement unit 31, a vibration identification unit 32, a stress calculation unit 33, a steady response measurement unit 34, and a generated stress evaluation unit 35. ing. In the present embodiment, the computer 27 (stress evaluation means) includes a vibration specifying unit 32, a stress calculation unit 33, a steady response measurement unit 34, and a generated stress evaluation unit 35.

[振動計測部]
振動計測部31は、図1に示した加速度計21、22、アンプ23、24、A/D変換器25、26から構成され、振動計測部31で計測した加速度の時刻歴波形がコンピュータ27に入力、記録され、当該計測データが振動計測部31、定常応答計測部34で用いられることになる。
[Vibration measurement unit]
The vibration measurement unit 31 includes the accelerometers 21 and 22, amplifiers 23 and 24, and A / D converters 25 and 26 illustrated in FIG. 1, and the time history waveform of acceleration measured by the vibration measurement unit 31 is stored in the computer 27. The measurement data is input and recorded, and the measurement data is used by the vibration measurement unit 31 and the steady response measurement unit 34.

[振動数特定部]
振動数特定部32では、本実施例の配管応力評価装置20を用いて、タッピング試験を行い、評価対象の配管(ここでは、分岐管12)の振動数ν1(固有振動数)を特定している。振動数特定部32について、図1〜図3及び図7〜図9を参照して説明する。
[Frequency specific part]
The frequency specifying unit 32 performs a tapping test using the pipe stress evaluation device 20 of the present embodiment, and specifies the frequency ν1 (natural frequency) of the pipe to be evaluated (here, the branch pipe 12). Yes. The frequency identification unit 32 will be described with reference to FIGS. 1 to 3 and FIGS. 7 to 9.

まず、図1に示すように、評価対象の分岐管12に加速度計21、22を設置する。ここでは、加速度計21が母管11に近い方に配置され、加速度計22が加速度計21より母管11から遠い方に配置される。そして、タッピングハンマー15を用いて、任意の振動となるタッピング加振を分岐管12に与え、加速度計21、22、アンプ23、24、A/D変換器25、26を用いて、加速度の時刻歴波形の計測データを取得する(ステップS1;タッピング試験)。取得した計測データから、図7(a)、(b)に示すように、分岐管12及び分岐管12の母管11側について、タッピング加振による加速度の時刻歴波形を取得することになる。   First, as shown in FIG. 1, accelerometers 21 and 22 are installed in the branch pipe 12 to be evaluated. Here, the accelerometer 21 is arranged closer to the mother tube 11, and the accelerometer 22 is arranged farther from the mother tube 11 than the accelerometer 21. The tapping hammer 15 is used to give the tapping excitation, which is an arbitrary vibration, to the branch pipe 12, and the acceleration time is measured using the accelerometers 21 and 22, amplifiers 23 and 24, and A / D converters 25 and 26. History waveform measurement data is acquired (step S1; tapping test). From the acquired measurement data, as shown in FIGS. 7A and 7B, a time history waveform of acceleration due to tapping excitation is acquired for the branch pipe 12 and the parent pipe 11 side of the branch pipe 12.

そして、図7(a)、(b)に示す加速度の時刻歴波形の波形を、周波数分析することにより、図8(a)、(b)に示すように、分岐管12及び分岐管12の母管11側について、加速度スペクトルを算出する(ステップS2;周波数分析)。ここで、図8(a)に示す分岐管12の加速度スペクトルを加速度スペクトルA(第1の加速度スペクトル)とし、図8(b)に示す分岐管12の母管11側の加速度スペクトルを加速度スペクトルB(第2の加速度スペクトル)とする。   Then, by analyzing the frequency of the acceleration time history waveform shown in FIGS. 7A and 7B, as shown in FIGS. 8A and 8B, the branch pipe 12 and the branch pipe 12 An acceleration spectrum is calculated for the mother tube 11 side (step S2; frequency analysis). Here, the acceleration spectrum of the branch pipe 12 shown in FIG. 8A is an acceleration spectrum A (first acceleration spectrum), and the acceleration spectrum of the branch pipe 12 shown in FIG. Let B (second acceleration spectrum).

算出した分岐管12の加速度スペクトルAと、分岐管12の母管11側の加速度スペクトルBを用いて、母管11と分岐管12の伝達関数F1(第1の伝達関数)を算出する。具体的には、分岐管12の加速度スペクトルAを、分岐管12の母管11側の加速度スペクトルBを用いて除算することにより、図9に示す伝達関数F1、即ち、周波数に対する加速度応答倍率の関係を求めることができる。図9に示す伝達関数F1のスペクトルピークに基づいて、分岐管12の振動モードと振動数ν1(第1の振動数)を特定することになる(ステップS3;特定)。ここで特定した振動モード及び振動数ν1が後述するステップ14で使用されることになる。   Using the calculated acceleration spectrum A of the branch pipe 12 and the acceleration spectrum B of the branch pipe 12 on the side of the mother pipe 11, a transfer function F1 (first transfer function) between the mother pipe 11 and the branch pipe 12 is calculated. Specifically, by dividing the acceleration spectrum A of the branch pipe 12 by using the acceleration spectrum B on the side of the mother pipe 11 of the branch pipe 12, the transfer function F1, shown in FIG. A relationship can be sought. Based on the spectrum peak of the transfer function F1 shown in FIG. 9, the vibration mode and the frequency ν1 (first frequency) of the branch pipe 12 are specified (step S3; specification). The vibration mode and the frequency ν1 specified here are used in step 14 described later.

[解析モデル決定部]
解析モデル決定部33では、モデル諸元の設定により配管の解析モデルを作成すると共に、振動数特定部32で特定した振動数ν1に解析モデルを用いて算出した振動数ν2が一致するように、解析モデルの諸元のチューニングを行い、最適な解析モデルを決定している。解析モデル決定部33について、図1、図2、図4及び図10〜図14を参照して説明する。
[Analysis model determination unit]
The analysis model determination unit 33 creates an analysis model of the pipe by setting the model specifications, and the frequency ν2 calculated using the analysis model matches the frequency ν1 specified by the frequency specification unit 32. The analysis model is tuned to determine the optimal analysis model. The analysis model determination unit 33 will be described with reference to FIGS. 1, 2, 4, and 10 to 14.

まず、解析モデルM1(第1の解析モデル)を作成するためのモデル諸元、例えば、配管、弁の諸元、局部剛性、サポートの諸元、計測点情報等を設定する(ステップS11;モデル諸元設定)。   First, model specifications for creating the analysis model M1 (first analysis model), for example, piping, valve specifications, local stiffness, support specifications, measurement point information, etc. are set (step S11; model). Specification settings).

配管、弁の諸元としては、具体的には、母管11、分岐管12、弁13の質量、剛性などを入力すると共に、母管11、分岐管12のサポート(固定箇所)の剛性などを設定する。例えば、図11に示すように、母管11が多スパン母管モデルとして、複数のサポート14に支持、固定されている場合には、母管11は、回転剛性及び水平剛性としてモデル化できるので、回転バネ41及び水平バネ42、43を用いてモデル化すると共に、母管11の質量を母管集中質量45としてモデル化する。母管集中質量45は、分岐管12を挟む2つのサポート14の間の母管11の質量の50%を目処とする。   As the specifications of the piping and valves, specifically, the masses, rigidity, etc. of the mother pipe 11, the branch pipe 12, and the valve 13 are input, and the rigidity of the support (fixed location) of the mother pipe 11, the branch pipe 12 and the like. Set. For example, as shown in FIG. 11, when the mother pipe 11 is supported and fixed to a plurality of supports 14 as a multi-span mother pipe model, the mother pipe 11 can be modeled as rotational rigidity and horizontal rigidity. Modeling is performed using the rotary spring 41 and the horizontal springs 42 and 43, and the mass of the mother pipe 11 is modeled as a mother pipe concentrated mass 45. The mother pipe concentrated mass 45 is targeted at 50% of the mass of the mother pipe 11 between the two supports 14 sandwiching the branch pipe 12.

又、分岐管12については、例えば、図10に示すように、梁要素J1としてモデル化し、弁13は、2つの弁集中質量46、47としてモデル化する。弁集中質量46は弁13本体を、弁集中質量47は弁13のハンドル13aをモデル化したものである。なお、形状線L1は、後述する図14においてハンドル13aの動きを示すために図示したものである。   For example, as shown in FIG. 10, the branch pipe 12 is modeled as a beam element J1, and the valve 13 is modeled as two valve concentrated masses 46 and 47. The valve concentrated mass 46 models the main body of the valve 13, and the valve concentrated mass 47 models the handle 13 a of the valve 13. The shape line L1 is illustrated to show the movement of the handle 13a in FIG. 14 described later.

又、例えば、図12に示すように、母管11と分岐管12の取り合い部(接続部分)の局部剛性、つまり、分岐管12が接続される母管11の表面11aの部分の局部剛性は、回転剛性としてモデル化できるので、回転バネ44を用いてモデル化する。この際、母管11の半径と肉厚の比、母管11と分岐管12の半径比が判れば、図13に示すグラフから、局部剛性(回転剛性)の大きさを求めることができ、局部剛性(回転剛性)の大きさを設定することができる。   For example, as shown in FIG. 12, the local rigidity of the joint portion (connection portion) of the mother pipe 11 and the branch pipe 12, that is, the local rigidity of the surface 11a portion of the mother pipe 11 to which the branch pipe 12 is connected is Since it can be modeled as rotational rigidity, the rotational spring 44 is used for modeling. At this time, if the ratio of the radius and thickness of the mother pipe 11 and the radius ratio of the mother pipe 11 and the branch pipe 12 are known, the magnitude of local rigidity (rotational rigidity) can be obtained from the graph shown in FIG. The size of the local rigidity (rotational rigidity) can be set.

又、加速度計21、22を設置した計測点の情報、例えば、設置位置の情報等を設定する。   In addition, information on measurement points where the accelerometers 21 and 22 are installed, for example, information on installation positions is set.

以上の設定を行うことにより、評価対象の配管(分岐管12)に関して、図10に示す解析モデルM1を作成することになる(ステップS12;解析モデル作成)。   By performing the above settings, the analysis model M1 shown in FIG. 10 is created for the pipe to be evaluated (branch pipe 12) (step S12; analysis model creation).

作成した解析モデルM1を用いて、固有値解析を実施し、振動数ν2(第2の振動数)、振動モードを算出する(ステップS13;固有値解析)。例えば、図1で示した配管を対象にして、分岐管12及び弁13の固有値解析を行うと、図14(a)に示すように、1次モードとして、62.9Hzの振動数、又、図14(b)に示すように、2次モードとして、72.1Hzの振動数を算出することができる。なお、図14(a)、(b)において、点線は基本の形状を表し、実線はモード形状(振動時の形状)を表している。又、固有値解析としては、公知の解析方法を用いればよい。   Using the created analysis model M1, eigenvalue analysis is performed to calculate the frequency ν2 (second frequency) and the vibration mode (step S13; eigenvalue analysis). For example, when the eigenvalue analysis of the branch pipe 12 and the valve 13 is performed on the pipe shown in FIG. 1, the frequency of 62.9 Hz is obtained as the primary mode as shown in FIG. As shown in FIG. 14B, the frequency of 72.1 Hz can be calculated as the secondary mode. In FIGS. 14A and 14B, the dotted line represents the basic shape, and the solid line represents the mode shape (shape during vibration). As eigenvalue analysis, a known analysis method may be used.

そして、ステップS3で特定した振動数ν1とステップS13において算出した振動数ν2とを比較し、一致する場合には、後述するステップS31へ進み、一致しない場合には、ステップS15へ進む(ステップS14;比較)。ここでは、ステップS12で作成した解析モデルM1の精度を確認することになる。   Then, the frequency ν1 specified in step S3 is compared with the frequency ν2 calculated in step S13. If they match, the process proceeds to step S31 described later, and if they do not match, the process proceeds to step S15 (step S14). ; Comparison). Here, the accuracy of the analysis model M1 created in step S12 is confirmed.

振動数ν1と振動数ν2が一致しない場合には、配管、サポートの諸元を見直し、修正する(ステップS15;修正)。このとき、分岐管12については、その長さや肉厚などから、剛性、質量分布を正確に設定できるので、その諸元は変更せず、又、母管11と分岐管12の取り合い部の局部剛性も、図13で説明したように、正確に設定できるので、その諸元は変更しないが、モデル諸元の中で不確定要素となりやすい母管11の諸元(主に、母管11の水平剛性、回転剛性)について、特に、サポート14の形状(門型、片持ち梁型など)に影響を受け、正確に評価することが難しい、サポート14による水平剛性、回転剛性について、修正を行い、更に、場合によっては、母管11の母管集中質量45について、修正を行い、振動数ν1と振動数ν2が一致するまで、ステップS11〜S15を繰り返して、チューニングを行う。このような手順で解析モデルM2(第2の解析モデル)が決定される。決定された解析モデルM2は最適化されたものであり、解析モデルとしての精度が向上することになる。   If the frequency ν1 and the frequency ν2 do not match, the specifications of the piping and the support are reviewed and corrected (step S15; correction). At this time, since the rigidity and mass distribution of the branch pipe 12 can be accurately set based on the length and thickness of the branch pipe 12, the specifications thereof are not changed, and the local part of the joint portion between the mother pipe 11 and the branch pipe 12 is not changed. As described with reference to FIG. 13, the rigidity can be set accurately, so the specifications thereof are not changed, but the specifications of the mother pipe 11 that is likely to be an uncertain factor in the model specifications (mainly, the specifications of the mother pipe 11). Regarding horizontal rigidity and rotational rigidity, the horizontal rigidity and rotational rigidity of the support 14 are corrected, especially affected by the shape of the support 14 (gate type, cantilever type, etc.) and difficult to evaluate accurately. In some cases, the mother tube concentrated mass 45 of the mother tube 11 is corrected, and tuning is performed by repeating steps S11 to S15 until the frequency ν1 and the frequency ν2 coincide. The analysis model M2 (second analysis model) is determined by such a procedure. The determined analysis model M2 is optimized, and the accuracy as the analysis model is improved.

[定常応答計測部]
定常応答計測部34では、本実施例の配管応力評価装置20を用いて、実際の配管において、プラントの定常運転時の分岐管12の加速度計測を実施し、定常運転時のおける加速度スペクトルを算出している。定常応答計測部34について、図1、図2、図5及び図15を参照して説明する。
[Steady response measurement unit]
The steady response measuring unit 34 uses the pipe stress evaluation device 20 of the present embodiment to measure the acceleration of the branch pipe 12 during the steady operation of the plant in the actual pipe and calculate the acceleration spectrum during the steady operation. doing. The steady response measuring unit 34 will be described with reference to FIGS. 1, 2, 5, and 15.

まず、分岐管12に加速度計21、22を設置し、設置した加速度計21、22、アンプ23、24、A/D変換器25、26を用いて、プラントの定常運転時における加速度の時刻歴波形の計測データを取得する(ステップS21;計測)。加速度計21、22を設置する位置は、振動数特定部32におけるタッピング試験(ステップS1参照)と同じ位置でもよいが、計測点の設置位置の情報が明らかであれば、分岐管12の任意の位置でも良く、又、加速度計21、22のいずれか一方のみの計測でもよい。加速度計21、22を用いて取得した計測データから、図15(a)に示すように、分岐管12における加速度の時刻歴波形を得ることになる。   First, accelerometers 21 and 22 are installed in the branch pipe 12, and using the installed accelerometers 21 and 22, amplifiers 23 and 24, and A / D converters 25 and 26, a time history of acceleration during steady operation of the plant. Waveform measurement data is acquired (step S21; measurement). The position where the accelerometers 21 and 22 are installed may be the same position as the tapping test (see step S1) in the frequency identification unit 32. However, if information on the installation position of the measurement point is clear, any position of the branch pipe 12 may be set. The position may be measured, or only one of the accelerometers 21 and 22 may be measured. From the measurement data acquired using the accelerometers 21 and 22, a time history waveform of acceleration in the branch pipe 12 is obtained as shown in FIG.

そして、図15(a)に示す加速度の時刻歴波形の波形を、周波数分析することにより、図15(b)に示すように、分岐管12における加速度スペクトルを算出する(ステップS22;分析)。ここで、図15(b)に示す分岐管12の加速度スペクトルを加速度スペクトルC(第3の加速度スペクトル)とする。ここで分析した加速度スペクトルCが後述するステップ33で使用されることになる。   Then, by analyzing the frequency of the acceleration time history waveform shown in FIG. 15A, an acceleration spectrum in the branch pipe 12 is calculated as shown in FIG. 15B (step S22; analysis). Here, the acceleration spectrum of the branch pipe 12 shown in FIG. 15B is assumed to be an acceleration spectrum C (third acceleration spectrum). The acceleration spectrum C analyzed here is used in step 33 described later.

[発生応力評価部]
発生応力評価部35では、決定した解析モデルM2を用い、ステップS22で分析された加速度スペクトルCに基づいて、評価点で発生する応力を算出し、算出した応力から疲労評価を行っている。発生応力評価部35について、図1、図2、図6及び図16〜図18を参照して説明する。なお、以降の説明では、一例として、加速度計21を設置した分岐管12の基部の位置を評価点としているが、以降の手順では、分岐管12の任意の位置を評価点としても、当該評価点で発生する応力を算出して、疲労評価を行うことができる。
[Generated stress evaluation section]
The generated stress evaluation unit 35 uses the determined analysis model M2 to calculate the stress generated at the evaluation point based on the acceleration spectrum C analyzed in step S22, and performs fatigue evaluation from the calculated stress. The generated stress evaluation unit 35 will be described with reference to FIGS. 1, 2, 6 and 16 to 18. In the following description, as an example, the position of the base of the branch pipe 12 where the accelerometer 21 is installed is used as an evaluation point. However, in the following procedure, the evaluation is performed even if an arbitrary position of the branch pipe 12 is used as an evaluation point. Fatigue evaluation can be performed by calculating the stress generated at a point.

ステップS14で決定した解析モデルM2を用いて、周波数応答解析を実施する。例えば、図16に示すように、計測点(図16では、加速度計22を設置した位置)における加速度の周波数応答と、評価点(図16では、加速度計21を設置した位置)における応力の周波数応答を算出する(ステップS31;周波数応答解析)。ここで、図16に示す計測点での加速度スペクトルを加速度スペクトルD(第4の加速度スペクトル)とし、評価点での応力スペクトルを応力スペクトルE(第1の応力スペクトル)とする。   A frequency response analysis is performed using the analysis model M2 determined in step S14. For example, as shown in FIG. 16, the frequency response of acceleration at the measurement point (in FIG. 16, the position where the accelerometer 22 is installed) and the frequency of stress at the evaluation point (the position where the accelerometer 21 is installed in FIG. 16). A response is calculated (step S31; frequency response analysis). Here, the acceleration spectrum at the measurement point shown in FIG. 16 is the acceleration spectrum D (fourth acceleration spectrum), and the stress spectrum at the evaluation point is the stress spectrum E (first stress spectrum).

ステップS31で求めた計測点での加速度スペクトルDと、評価点での応力スペクトルEとを用いて、応力と加速度の伝達関数を算出する。具体的には、応力スペクトルEを、加速度スペクトルDを用いて除算することにより、図17に示す伝達関数、即ち、周波数に対する応力−加速度応答倍率の関係を求めることができる(ステップS32;伝達関数計算)。ここで、図17に示す伝達関数を伝達関数F2(第2の伝達関数)とする。   Using the acceleration spectrum D at the measurement point obtained in step S31 and the stress spectrum E at the evaluation point, a transfer function of stress and acceleration is calculated. Specifically, by dividing the stress spectrum E by using the acceleration spectrum D, the transfer function shown in FIG. 17, that is, the relationship between the stress-acceleration response magnification with respect to the frequency can be obtained (step S32; transfer function). Calculation). Here, the transfer function shown in FIG. 17 is defined as a transfer function F2 (second transfer function).

ステップS32で求めた伝達関数F2に、ステップS22で分析した加速度スペクトルCを乗算する。これにより、評価点における応力スペクトルを求めることができる。例えば、図17に示した伝達関数F2のグラフと図15(b)に示した加速度スペクトルCのグラフを乗算することにより、図18に示す周波数対応の応力スペクトルを求めることになる(ステップS33;応力スペクトル計算)。ここで、図18に示す応力スペクトルを応力スペクトルG(第2の応力スペクトル)とする。   The transfer function F2 obtained in step S32 is multiplied by the acceleration spectrum C analyzed in step S22. Thereby, the stress spectrum in an evaluation point can be calculated | required. For example, the stress spectrum corresponding to the frequency shown in FIG. 18 is obtained by multiplying the graph of the transfer function F2 shown in FIG. 17 by the graph of the acceleration spectrum C shown in FIG. 15B (step S33; Stress spectrum calculation). Here, let the stress spectrum shown in FIG. 18 be a stress spectrum G (second stress spectrum).

ステップS33で求めた応力スペクトルGを周波数積分することにより、応力の標準偏差σを算出する。標準偏差を3倍した値は最大値に近い値となるので、ここでは、算出した標準偏差σを3倍した値、即ち、[3σ]をピーク応力として算出する(ステップS34;ピーク応力算出)。   The stress standard deviation σ is calculated by frequency integrating the stress spectrum G obtained in step S33. Since the value obtained by multiplying the standard deviation by three is close to the maximum value, the value obtained by multiplying the calculated standard deviation σ by three, that is, [3σ] is calculated as the peak stress (step S34; peak stress calculation). .

ステップS34で算出したピーク応力3σを、分岐管12の許容応力と比較することにより、対象配管である分岐管12の応力評価を行う(ステップS35;応力評価)。このようにして、ピーク応力3σが分岐管12の許容応力の範囲内であるかどうかを判断することにより、高い信頼性で配管の疲労評価を行うことになる。   By comparing the peak stress 3σ calculated in step S34 with the allowable stress of the branch pipe 12, the stress of the branch pipe 12 that is the target pipe is evaluated (step S35; stress evaluation). Thus, by judging whether or not the peak stress 3σ is within the allowable stress range of the branch pipe 12, the fatigue evaluation of the pipe is performed with high reliability.

以上説明してきたように、本実施例の配管応力評価装置20では、加速度計21、22により取得した加速度の時刻歴波形のデータを用いているため、現場での計測の手間を軽減することができる。   As described above, in the piping stress evaluation apparatus 20 of the present embodiment, the time history waveform data of acceleration acquired by the accelerometers 21 and 22 is used. it can.

又、発生応力評価部35において、解析モデルM2は、分岐管12のモデルと共に、母管11のモデルまで含まれており、更に、母管11と分岐管12の取り合い部の局部剛性もモデル化されているので、解析モデルとしての精度が格段に向上すると共に、当該解析モデルM2から算定する応力の精度も向上し、算定した応力、そして、その疲労評価の信頼性が高いものとなる。なお、母管11と分岐管12の取り合い部の局部剛性は、図13で説明したように、配管の条件でグラフ化されており、容易に解析モデルM2に反映することができる。   Further, in the generated stress evaluation unit 35, the analysis model M2 includes not only the model of the branch pipe 12 but also the model of the mother pipe 11, and further models the local rigidity of the joint between the mother pipe 11 and the branch pipe 12. Therefore, the accuracy as the analysis model is remarkably improved, the accuracy of the stress calculated from the analysis model M2 is also improved, and the calculated stress and the reliability of the fatigue evaluation become high. Note that the local rigidity of the joint between the mother pipe 11 and the branch pipe 12 is graphed under the piping conditions as described with reference to FIG. 13, and can be easily reflected in the analysis model M2.

加えて、発生応力評価部35で用いる解析モデルM2は、タッピング試験により特定した振動数に一致するように、チューニング(最適化)されているので、実際の配管に適応した精度の良い解析モデルとして構築される。このようにして最適化された精度の良い解析モデルM2を用い、更に、実際の定常運転時の加速度スペクトルCを用いて、評価点の応力を算出するので、実際の定常運転時に応じた応力を精度良く算出することになり、信頼性の高い疲労評価を行うことができる。   In addition, the analysis model M2 used in the generated stress evaluation unit 35 is tuned (optimized) so as to match the frequency specified by the tapping test. Built. The stress at the evaluation point is calculated by using the optimized analysis model M2 optimized in this way, and further using the acceleration spectrum C at the actual steady operation, so the stress corresponding to the actual steady operation is calculated. The calculation is performed with high accuracy, and a highly reliable fatigue evaluation can be performed.

又、解析モデルM2を用いて、応力を算出しているので、分岐管12の基部に限らず、分岐管12の任意の評価点においても、応力を算出することができ、その疲労評価を行うことができる。   Moreover, since the stress is calculated using the analysis model M2, the stress can be calculated not only at the base of the branch pipe 12, but also at any evaluation point of the branch pipe 12, and the fatigue evaluation is performed. be able to.

(実施例2)
本実施例は、実施例1で説明した配管応力評価装置20(図1、図2参照)を用いると共に、実施例1で説明した配管応力評価方法(図3〜図6参照)と基本的に同じ手順を実施するものであるが、分岐管12に取り付けてある弁13の解析モデルM1が実施例1(図10参照)とは相違する。ここでは、解析モデルM1の相違点について説明し、同等の点については、その説明を省略する。
(Example 2)
This embodiment uses the pipe stress evaluation apparatus 20 (see FIGS. 1 and 2) described in the first embodiment and is basically the same as the pipe stress evaluation method (see FIGS. 3 to 6) described in the first embodiment. Although the same procedure is performed, the analysis model M1 of the valve 13 attached to the branch pipe 12 is different from that of the first embodiment (see FIG. 10). Here, differences between the analysis models M1 will be described, and descriptions of equivalent points will be omitted.

実施例1では、図10で説明したように、分岐管12を梁要素J1としてモデル化すると共に、弁13を弁集中質量46、47としてモデル化した。一方、本実施例では、弁13を集中質量としてモデル化するのではなく、弁分布質量48、49としてモデル化している。例えば、図19では、弁13全体の質量の50%を弁13の重心に、その10%をハンドル13aの重心に、残りの40%を弁13全体に分布させることで、弁分布質量48、49としてモデル化している。このように、本実施例では、ハンドル13aの部分も弁分布質量49としてモデル化の対象としている。このような弁分布質量48、49は、弁13の形状に応じて任意に変更すればよい。   In Example 1, the branch pipe 12 was modeled as the beam element J1 and the valve 13 was modeled as the valve concentrated masses 46 and 47 as described with reference to FIG. On the other hand, in this embodiment, the valve 13 is not modeled as a concentrated mass, but is modeled as valve distribution masses 48 and 49. For example, in FIG. 19, 50% of the mass of the entire valve 13 is distributed at the center of gravity of the valve 13, 10% is distributed at the center of gravity of the handle 13a, and the remaining 40% is distributed throughout the valve 13. 49 as a model. As described above, in the present embodiment, the handle 13 a is also modeled as the valve distribution mass 49. Such valve distribution masses 48 and 49 may be arbitrarily changed according to the shape of the valve 13.

このように、分岐管12の弁13を弁分布質量48、49でモデル化することにより、分岐管12自体の振動特性の特定精度が向上する。その結果、より信頼性の高い応力を算出することができ、その疲労評価の信頼性もより高いものになる。   Thus, by specifying the valve 13 of the branch pipe 12 with the valve distribution masses 48 and 49, the accuracy of specifying the vibration characteristics of the branch pipe 12 itself is improved. As a result, stress with higher reliability can be calculated, and the reliability of the fatigue evaluation becomes higher.

本発明は、プラント等における配管の分岐管の疲労評価に好適なものである。   The present invention is suitable for fatigue evaluation of a branch pipe of a pipe in a plant or the like.

11 母管
12 分岐管
13 弁
15 タッピングハンマー
21、22 加速度計
23、24 アンプ
25、26 A/D変換器
27 コンピュータ
31 振動計測部
32 振動数特定部
33 解析モデル決定部
34 定常応答計測部
35 発生応力評価部
41 回転バネ(母管回転剛性)
42、43 水平バネ(母管水平剛性)
44 回転バネ(母管局部剛性(回転剛性))
45 母管集中質量
46、47 弁集中質量
48、49 弁分布質量
J1 梁要素
DESCRIPTION OF SYMBOLS 11 Mother pipe 12 Branch pipe 13 Valve 15 Tapping hammer 21, 22 Accelerometer 23, 24 Amplifier 25, 26 A / D converter 27 Computer 31 Vibration measurement part 32 Frequency identification part 33 Analysis model determination part 34 Steady response measurement part 35 Generated stress evaluation unit 41 Rotating spring (Rolling rigidity of mother pipe)
42, 43 Horizontal spring (base pipe horizontal rigidity)
44 Rotating spring (Home pipe local rigidity (rotational rigidity))
45 Concentrated mass of mother pipe 46, 47 Concentrated mass of valve 48, 49 Distributed mass of valve J1 Beam element

Claims (6)

プラントの母管から分岐した分岐管に、或いは、前記母管及び前記分岐管に取り付ける複数の加速度計と、
前記加速度計で計測した計測値を用いて、前記分岐管に発生する応力を算出し、評価を行う応力評価手段とを備えた配管応力評価装置において、
前記応力評価手段は、
任意の振動を前記分岐管に与えたときの前記複数の加速度計における加速度の時刻歴波形から、周波数分析を行って、前記母管から遠い方の前記加速度計での第1の加速度スペクトルと前記母管に近い方又は前記母管の前記加速度計での第2の加速度スペクトルとを算出し、前記第1の加速度スペクトルを前記第2の加速度スペクトルで除算して、第1の伝達関数を算出し、前記第1の伝達関数のスペクトルピークから前記分岐管の第1の振動数を特定する振動数特定部と、
前記母管を回転剛性、水平剛性及び集中質量でモデル化し、前記分岐管を梁要素でモデル化し、前記母管と前記分岐管の接続部分を回転剛性でモデル化した第1の解析モデルを作成し、前記第1の解析モデルに基づいて、固有値解析を行って、前記分岐管の第2の振動数を算出し、前記第2の振動数が前記第1の振動数に一致するように、前記母管の回転剛性、水平剛性の設定を変更して、前記第2の振動数が前記第1の振動数に一致する第2の解析モデルを決定する解析モデル決定部と、
前記プラントの定常運転時において、前記分岐管の任意の計測点に配置した前記加速度計における加速度の時刻歴波形から、周波数分析を行って、第3の加速度スペクトルを算出する定常応答計測部と、
前記第2の解析モデルを用い、当該第2の解析モデルの周波数応答解析を行って、前記計測点における第4の加速度スペクトルと前記分岐管の任意の評価点における第1の応力スペクトルとを算出し、前記第1の応力スペクトルを前記第4の加速度スペクトルで除算して、第2の伝達関数を算出し、前記第2の伝達関数に前記第3の加速度スペクトルを乗算して、前記評価点における第2の応力スペクトルを算出し、前記第2の応力スペクトルを周波数積分してピーク応力を算出し、前記ピーク応力と前記分岐管の許容応力とを比較することにより、前記分岐管に発生する応力の評価を行う発生応力評価部とを有することを特徴とする配管応力評価装置。
A plurality of accelerometers attached to the branch pipe branched from the mother pipe of the plant, or to the mother pipe and the branch pipe;
In the pipe stress evaluation apparatus provided with stress evaluation means for calculating and evaluating the stress generated in the branch pipe using the measurement value measured by the accelerometer,
The stress evaluation means includes
The frequency analysis is performed from the time history waveforms of acceleration in the plurality of accelerometers when arbitrary vibration is applied to the branch pipe, and the first acceleration spectrum in the accelerometer far from the mother pipe and the first acceleration spectrum A first acceleration function is calculated by calculating a second acceleration spectrum that is closer to the mother tube or the accelerometer of the mother tube, and dividing the first acceleration spectrum by the second acceleration spectrum. A frequency specifying unit for specifying the first frequency of the branch pipe from the spectrum peak of the first transfer function;
The first pipe is modeled with rotational rigidity, horizontal rigidity, and concentrated mass, the branch pipe is modeled with a beam element, and a first analysis model is created in which the connection between the mother pipe and the branch pipe is modeled with rotational rigidity. Then, based on the first analysis model, eigenvalue analysis is performed to calculate a second frequency of the branch pipe, so that the second frequency matches the first frequency, An analysis model determination unit that determines a second analysis model in which the second frequency matches the first frequency by changing the rotational rigidity and horizontal rigidity of the mother pipe;
During steady operation of the plant, from a time history waveform of acceleration in the accelerometer arranged at an arbitrary measurement point of the branch pipe, by performing frequency analysis, a steady response measurement unit that calculates a third acceleration spectrum;
Using the second analysis model, frequency response analysis of the second analysis model is performed to calculate a fourth acceleration spectrum at the measurement point and a first stress spectrum at an arbitrary evaluation point of the branch pipe. And dividing the first stress spectrum by the fourth acceleration spectrum to calculate a second transfer function, multiplying the second transfer function by the third acceleration spectrum, and Is generated in the branch pipe by calculating a second stress spectrum at, calculating a peak stress by frequency integrating the second stress spectrum, and comparing the peak stress with an allowable stress of the branch pipe A piping stress evaluation apparatus, comprising: a generated stress evaluation unit that performs stress evaluation.
請求項1に記載の配管応力評価装置において、
前記解析モデル決定部は、前記分岐管に弁が設けられている場合には、前記弁を集中質量でモデル化して、前記第1の解析モデルを作成することを特徴とする配管応力評価装置。
In the piping stress evaluation apparatus according to claim 1,
The said analysis model determination part models the said valve by concentrated mass, and produces the said 1st analysis model, when the valve is provided in the said branch pipe, The piping stress evaluation apparatus characterized by the above-mentioned.
請求項1に記載の配管応力評価装置において、
前記解析モデル決定部は、前記分岐管に弁が設けられている場合には、前記弁のハンドルを含めて、前記弁全体を分布質量でモデル化して、前記第1の解析モデルを作成することを特徴とする配管応力評価装置。
In the piping stress evaluation apparatus according to claim 1,
The analysis model determination unit creates the first analysis model by modeling the entire valve including a handle of the valve with a distributed mass when the branch pipe is provided with a valve. Piping stress evaluation device characterized by this.
プラントの母管から分岐した分岐管に、或いは、前記母管及び前記分岐管に取り付ける複数の加速度計を用いて、任意の振動を前記分岐管に与えたときの前記複数の加速度計における加速度の時刻歴波形から、周波数分析を行って、前記母管から遠い方の前記加速度計での第1の加速度スペクトルと前記母管に近い方又は前記母管の前記加速度計での第2の加速度スペクトルとを算出し、前記第1の加速度スペクトルを前記第2の加速度スペクトルで除算して、第1の伝達関数を算出し、前記第1の伝達関数のスペクトルピークから前記分岐管の第1の振動数を特定し、
前記母管を回転剛性、水平剛性及び集中質量でモデル化し、前記分岐管を梁要素でモデル化し、前記母管と前記分岐管の接続部分を回転剛性でモデル化した第1の解析モデルを作成し、前記第1の解析モデルに基づいて、固有値解析を行って、前記分岐管の第2の振動数を算出し、前記第2の振動数が前記第1の振動数に一致するように、前記母管の回転剛性、水平剛性の設定を変更して、前記第2の振動数が前記第1の振動数に一致する第2の解析モデルを決定し、
前記プラントの定常運転時において、前記分岐管の任意の計測点に配置した前記加速度計における加速度の時刻歴波形から、周波数分析を行って、第3の加速度スペクトルを算出し、
前記第2の解析モデルを用い、当該第2の解析モデルの周波数応答解析を行って、前記計測点における第4の加速度スペクトルと前記分岐管の任意の評価点における第1の応力スペクトルとを算出し、前記第1の応力スペクトルを前記第4の加速度スペクトルで除算して、第2の伝達関数を算出し、前記第2の伝達関数に前記第3の加速度スペクトルを乗算して、前記評価点における第2の応力スペクトルを算出し、前記第2の応力スペクトルを周波数積分してピーク応力を算出し、前記ピーク応力と前記分岐管の許容応力とを比較することにより、前記分岐管に発生する応力の評価を行うことを特徴とする配管応力評価方法。
The acceleration of the plurality of accelerometers when an arbitrary vibration is applied to the branch pipe, using a branch pipe branched from the main pipe of the plant or a plurality of accelerometers attached to the main pipe and the branch pipe. A frequency analysis is performed from the time history waveform, and a first acceleration spectrum at the accelerometer far from the mother tube and a second acceleration spectrum at the accelerometer near or near the mother tube. The first acceleration spectrum is divided by the second acceleration spectrum to calculate a first transfer function, and the first vibration of the branch pipe is calculated from the spectrum peak of the first transfer function. Identify the number,
The first pipe is modeled with rotational rigidity, horizontal rigidity, and concentrated mass, the branch pipe is modeled with a beam element, and a first analysis model is created in which the connection between the mother pipe and the branch pipe is modeled with rotational rigidity. Then, based on the first analysis model, eigenvalue analysis is performed to calculate a second frequency of the branch pipe, so that the second frequency matches the first frequency, Changing the setting of rotational stiffness and horizontal stiffness of the mother pipe to determine a second analysis model in which the second frequency matches the first frequency;
During steady operation of the plant, from the time history waveform of acceleration in the accelerometer arranged at an arbitrary measurement point of the branch pipe, frequency analysis is performed to calculate a third acceleration spectrum,
Using the second analysis model, frequency response analysis of the second analysis model is performed to calculate a fourth acceleration spectrum at the measurement point and a first stress spectrum at an arbitrary evaluation point of the branch pipe. And dividing the first stress spectrum by the fourth acceleration spectrum to calculate a second transfer function, multiplying the second transfer function by the third acceleration spectrum, and Is generated in the branch pipe by calculating a second stress spectrum at, calculating a peak stress by frequency integrating the second stress spectrum, and comparing the peak stress with an allowable stress of the branch pipe A pipe stress evaluation method characterized by performing stress evaluation.
請求項4に記載の配管応力評価方法において、
前記分岐管に弁が設けられている場合には、前記弁を集中質量でモデル化して、前記第1の解析モデルを作成することを特徴とする配管応力評価方法。
In the piping stress evaluation method according to claim 4,
A pipe stress evaluation method characterized in that, when a valve is provided in the branch pipe, the first analysis model is created by modeling the valve with a concentrated mass.
請求項4に記載の配管応力評価方法において、
前記分岐管に弁が設けられている場合には、前記弁のハンドルを含めて、前記弁全体を分布質量でモデル化して、前記第1の解析モデルを作成することを特徴とする配管応力評価方法。
In the piping stress evaluation method according to claim 4,
When the branch pipe is provided with a valve, pipe stress evaluation is characterized in that the first analysis model is created by modeling the entire valve including a handle of the valve with a distributed mass. Method.
JP2012216401A 2012-09-28 2012-09-28 Pipe stress evaluation apparatus and method Active JP5946064B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2012216401A JP5946064B2 (en) 2012-09-28 2012-09-28 Pipe stress evaluation apparatus and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2012216401A JP5946064B2 (en) 2012-09-28 2012-09-28 Pipe stress evaluation apparatus and method

Publications (2)

Publication Number Publication Date
JP2014070976A true JP2014070976A (en) 2014-04-21
JP5946064B2 JP5946064B2 (en) 2016-07-05

Family

ID=50746318

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012216401A Active JP5946064B2 (en) 2012-09-28 2012-09-28 Pipe stress evaluation apparatus and method

Country Status (1)

Country Link
JP (1) JP5946064B2 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2550192A (en) * 2016-05-12 2017-11-15 Xodus Group Ltd Pipework Fatigue Lifetime Measurement
KR101862729B1 (en) * 2016-12-13 2018-05-31 사단법인 캠틱종합기술원 Apparatus for analyzing vibration characteristics of tool using ultrasonic excitation
CN112735237A (en) * 2020-12-30 2021-04-30 许昌学院 Mechanical vibration test device
CN116680848A (en) * 2023-06-14 2023-09-01 西南石油大学 Pipeline suspending section safety evaluation system, device and medium
US12055453B2 (en) 2021-09-15 2024-08-06 Kabushiki Kaisha Toshiba Stress estimation method for machine structure and monitoring method for machine structure

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09145530A (en) * 1995-11-20 1997-06-06 Mitsubishi Heavy Ind Ltd Fatigue evaluation apparatus for piping system
JP2001153719A (en) * 1999-11-24 2001-06-08 Toshiba Corp Analysis system for vibration characteristic

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09145530A (en) * 1995-11-20 1997-06-06 Mitsubishi Heavy Ind Ltd Fatigue evaluation apparatus for piping system
JP2001153719A (en) * 1999-11-24 2001-06-08 Toshiba Corp Analysis system for vibration characteristic

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2550192A (en) * 2016-05-12 2017-11-15 Xodus Group Ltd Pipework Fatigue Lifetime Measurement
KR101862729B1 (en) * 2016-12-13 2018-05-31 사단법인 캠틱종합기술원 Apparatus for analyzing vibration characteristics of tool using ultrasonic excitation
CN112735237A (en) * 2020-12-30 2021-04-30 许昌学院 Mechanical vibration test device
CN112735237B (en) * 2020-12-30 2023-01-17 许昌学院 Mechanical vibration test device
US12055453B2 (en) 2021-09-15 2024-08-06 Kabushiki Kaisha Toshiba Stress estimation method for machine structure and monitoring method for machine structure
CN116680848A (en) * 2023-06-14 2023-09-01 西南石油大学 Pipeline suspending section safety evaluation system, device and medium
CN116680848B (en) * 2023-06-14 2023-12-19 西南石油大学 Pipeline suspending section safety evaluation system, device and medium

Also Published As

Publication number Publication date
JP5946064B2 (en) 2016-07-05

Similar Documents

Publication Publication Date Title
CN108195535B (en) Bolt joint looseness detection method and system based on nonlinear excitation characteristics
Lees et al. Model-based identification of rotating machines
JP5946064B2 (en) Pipe stress evaluation apparatus and method
CN107387333B (en) Method and device for determining fatigue damage of tower component of wind generating set
Liu et al. Weighted contribution rate of nonlinear output frequency response functions and its application to rotor system fault diagnosis
KR101541978B1 (en) Apparatus for detectng the tube wall thinning and method thereof
Siringoringo et al. Noncontact operational modal analysis of structural members by laser Doppler vibrometer
Rolfes et al. Integral SHM-system for offshore wind turbines using smart wireless sensors
Foti Dynamic identification techniques to numerically detect the structural damage
CA2824714A1 (en) Method for improving determination of mode shapes for a mechanical structure and applications hereof
Hong et al. Comprehensive comparison of macro-strain mode and displacement mode based on different sensing technologies
CN112162034B (en) Steel-concrete combined beam damage identification method applying structural noise
KR101653116B1 (en) Strain estimation system, global and local safety evaluation system, and method for structure through impact hammer test
EP2244081A2 (en) Structural integrity monitoring system
Teloli et al. Good practices for designing and experimental testing of dynamically excited jointed structures: The Orion beam
Lukasiewicz et al. Calculation of strains in dents based on high resolution in-line caliper survey
JP6075543B2 (en) Nonlinear response diagnostic apparatus and diagnostic method
JP3145625B2 (en) Piping system fatigue evaluation device
KR20150004127A (en) Method of Estimating Displacement of a Structure based on Acceleration and Strain
EP3708992B1 (en) Estimating device, estimating method, and program storing medium
EP3764074A1 (en) Pipeline diagnosing device, asset management device, pipeline diagnosing method, and recording medium
Pandiripalli Repeatability and reproducibility studies: A comparison of techniques
JP5390447B2 (en) Method for estimating noise countermeasure effect
EP4372348A2 (en) Tension in post-tensioned rods
JP2019100873A (en) Analysis device, analysis system, and analysis method

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20150818

TRDD Decision of grant or rejection written
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20160420

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

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20160426

RD03 Notification of appointment of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7423

Effective date: 20160426

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160524

R151 Written notification of patent or utility model registration

Ref document number: 5946064

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151