JP5251911B2 - Residual stress calculation device, residual stress measurement device, residual stress calculation method, residual stress measurement method, and program - Google Patents
Residual stress calculation device, residual stress measurement device, residual stress calculation method, residual stress measurement method, and program Download PDFInfo
- Publication number
- JP5251911B2 JP5251911B2 JP2010066987A JP2010066987A JP5251911B2 JP 5251911 B2 JP5251911 B2 JP 5251911B2 JP 2010066987 A JP2010066987 A JP 2010066987A JP 2010066987 A JP2010066987 A JP 2010066987A JP 5251911 B2 JP5251911 B2 JP 5251911B2
- Authority
- JP
- Japan
- Prior art keywords
- plastic deformation
- metal plate
- residual stress
- parameter
- sound
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Description
本発明は、残留応力算出装置、残留応力測定装置、残留応力算出方法、残留応力測定方法およびプログラムに関する。 The present invention relates to a residual stress calculation device, a residual stress measurement device, a residual stress calculation method, a residual stress measurement method, and a program.
金属板等の被検査体の内部に残留している応力（残留応力）は、複雑で不均一な応力分布を生じ、被検査体の機械的な性質に大きな影響を与えるものであるため、残留応力の大きさを正確に測定することは、被検査体の性能を維持する上で重要な技術である。そのため、金属板等の被検査体に残留している応力を測定するための技術について、様々な検討がなされており、各種の方法が提案されている。 Residual stress (residual stress) remaining inside the object to be inspected, such as a metal plate, creates a complex and uneven stress distribution and greatly affects the mechanical properties of the object to be inspected. Measuring the magnitude of the stress accurately is an important technique for maintaining the performance of the object to be inspected. For this reason, various studies have been made on techniques for measuring the stress remaining in an object to be inspected, such as a metal plate, and various methods have been proposed.
例えば、以下の特許文献１及び特許文献２には、Ｘ線又は中性子線を被検査体に照射し、得られた回折パターンから結晶の格子間隔を測定して、被検査体の応力を測定する方法が開示されている。また、以下の非特許文献１には、被検査体中を伝播する縦波及び横波の音速から音弾性パラメータを測定し、予め測定しておいた主応力和と主応力差の関係から主応力を算出する方法が開示されている。ここで、上記音弾性パラメータとは、主応力和、又は、主応力差と相関のある複数の超音波音速から算出されるパラメータのことである。
For example, in
しかしながら、上記特許文献１に記載の方法では、被検査体表面の残留応力しか測定することができないという問題がある。また、上記特許文献２に記載の方法では、被検査体内部の残留応力も測定可能ではあるが、測定装置が大規模になってしまうという問題がある。また、上記特許文献１及び特許文献２に記載の方法は、Ｘ線又は中性子線を利用する方法であるため、被爆の危険性があるという問題も存在する。
However, the method described in
他方、上記非特許文献１の方法は、比較的簡単な測定装置で被検査体内部の残留応力を測定可能ではあるが、被検査体が塑性変形している場合には、測定誤差が大きくなるという問題があった。
On the other hand, although the method of
そこで、本発明は、上記問題に鑑みてなされたものであり、本発明の目的とするところは、塑性変形している被検査体の内部に残留している応力を、安全かつ正確に測定することが可能な、残留応力算出装置、残留応力測定装置、残留応力算出方法、残留応力測定方法およびプログラムを提供することにある。 Therefore, the present invention has been made in view of the above problems, and an object of the present invention is to safely and accurately measure the stress remaining in the plastically deformed object. Another object is to provide a residual stress calculation device, a residual stress measurement device, a residual stress calculation method, a residual stress measurement method, and a program.
上記課題を解決するために、本発明のある観点によれば、塑性変形した金属板中に複数の超音波を発生させることで測定された金属板中の複数の超音波音速に関する音速情報を取得する音速情報取得部と、前記金属板の塑性変形量又は塑性変形量に換算可能な値である塑性変形特徴量を算出する塑性変形特徴量算出部と、前記音速情報取得部が取得した前記音速情報に基づいて、前記金属板中の音弾性パラメータを算出する音弾性パラメータ算出部と、前記塑性変形特徴量算出部により算出された前記塑性変形特徴量を利用して、前記音弾性パラメータの応力による変化を表す係数の補正を行う音弾性パラメータ係数補正部と、前記音弾性パラメータ算出部が算出した音弾性パラメータと、前記音弾性パラメータ係数補正部が補正した係数とに基づいて、前記金属板に残留している応力を算出する残留応力算出部と、を備える残留応力算出装置が提供される。 In order to solve the above-described problem, according to one aspect of the present invention, sound speed information related to a plurality of ultrasonic sound speeds in a metal plate measured by generating a plurality of ultrasonic waves in a plastically deformed metal plate is obtained. The sound velocity information acquisition unit, the plastic deformation amount of the metal plate or a plastic deformation feature amount calculation unit that calculates a plastic deformation feature amount that is a value that can be converted into the plastic deformation amount, and the sound speed acquired by the sound speed information acquisition unit. Based on the information, the acoustoelastic parameter calculating unit for calculating the acoustoelastic parameter in the metal plate, and the plastic deformation feature amount calculated by the plastic deformation feature amount calculating unit, the stress of the acoustoelastic parameter A sound elasticity parameter coefficient correction unit that corrects a coefficient representing a change caused by the sound elasticity parameter, a sound elasticity parameter calculated by the sound elasticity parameter calculation unit, a coefficient corrected by the sound elasticity parameter coefficient correction unit, Based on the residual stress calculation unit for calculating a stress remaining in the metal plate, the residual stress calculating device with a provided.
前記残留応力算出装置は、前記係数と前記塑性変形特徴量との関係を表す相関情報が格納された記憶部を更に備え、前記音弾性パラメータ係数補正部は、前記記憶部に格納された前記相関情報と、前記塑性変形特徴量算出部が算出した前記塑性変形特徴量とを利用して、前記係数を決定することが好ましい。 The residual stress calculation device further includes a storage unit storing correlation information representing a relationship between the coefficient and the plastic deformation feature amount, and the acoustoelastic parameter coefficient correction unit is the correlation stored in the storage unit. Preferably, the coefficient is determined using information and the plastic deformation feature amount calculated by the plastic deformation feature amount calculation unit.
前記音弾性パラメータ算出部は、前記音速情報を利用して、前記音弾性パラメータとして音速比Ｒ及び音響複屈折Ｂを算出し、前記音弾性パラメータ係数補正部は、前記塑性変形特徴量算出部が算出した前記塑性変形特徴量に基づいて、以下の式１における係数ＣＲ及びＲ０、並びに、以下の式２における係数ＣＡ及びＢ０を決定することが好ましい。
The acoustoelastic parameter calculation unit calculates the sound velocity ratio R and the acoustic birefringence B as the acoustoelastic parameters using the sound velocity information, and the acoustoelastic parameter coefficient correction unit includes the plastic deformation feature amount calculation unit. calculated based on the plastic deformation characteristic quantity, the coefficient C R and R 0 in
ここで、下記式１及び式２において、σ１及びσ２は、金属板中に残留している主応力を表す。
Here, in the following
前記残留応力算出部は、前記音弾性パラメータ算出部が算出した音速比Ｒ及び音響複屈折Ｂと、前記音弾性パラメータ係数補正部が決定した前記係数Ｒ０、ＣＲ、Ｂ０及びＣＡと、を利用し、前記式１及び前記式２に基づいて金属板に残留している主応力を算出することが好ましい。
The residual stress calculation unit includes a sound speed ratio R and an acoustic birefringence B calculated by the acoustoelastic parameter calculation unit, and the coefficients R 0 , C R , B 0 and C A determined by the acoustoelastic parameter coefficient correction unit. It is preferable to calculate the principal stress remaining on the metal plate based on the
前記塑性変形特徴量は、前記金属板に生じた塑性歪であってもよい。
また、前記塑性変形特徴量は、塑性変形に伴う前記金属板の板厚の変化の割合である減肉比であってもよい。
The plastic deformation feature amount may be a plastic strain generated in the metal plate.
The plastic deformation feature amount may be a thickness reduction ratio that is a rate of change in the thickness of the metal plate accompanying plastic deformation.
また、上記課題を解決するために、本発明の別の観点によれば、塑性変形した金属板中に複数の超音波を発生させることで、当該金属板中を伝播する前記複数の超音波の音速を測定する音速測定装置と、前記音速測定装置が測定した音速の測定結果を利用して、前記塑性変形した金属板中に残留する応力を算出する残留応力算出装置と、を備え、前記残留応力算出装置は、前記音速測定装置から、前記塑性変形した金属板中を伝播する複数の超音波音速に関する測定結果である音速情報を取得する音速情報取得部と、前記金属板の塑性変形量又は塑性変形量に換算可能な値である塑性変形特徴量を算出する塑性変形特徴量算出部と、前記音速情報取得部が取得した前記音速情報に基づいて、前記金属板中の音弾性パラメータを算出する音弾性パラメータ算出部と、前記塑性変形特徴量算出部により算出された前記塑性変形特徴量を利用して、前記音弾性パラメータの応力による変化を表す係数の補正を行う音弾性パラメータ係数補正部と、前記音弾性パラメータ算出部が算出した音弾性パラメータと、前記音弾性パラメータ係数補正部が補正した係数とに基づいて、前記金属板に残留している応力を算出する残留応力算出部と、を有する残留応力測定装置が提供される。 In order to solve the above problem, according to another aspect of the present invention, by generating a plurality of ultrasonic waves in a plastically deformed metal plate, the plurality of ultrasonic waves propagating in the metal plate are generated. A sonic velocity measuring device for measuring the sonic velocity, and a residual stress calculating device for calculating a residual stress in the plastically deformed metal plate using a measurement result of the sonic velocity measured by the sonic velocity measuring device. The stress calculation device includes a sound speed information acquisition unit that acquires sound speed information that is a measurement result regarding a plurality of ultrasonic sound speeds propagating in the plastically deformed metal plate from the sound speed measurement device, and a plastic deformation amount of the metal plate or Based on the sound velocity information acquired by the plastic deformation feature amount calculating unit that calculates a plastic deformation feature amount that is a value that can be converted into a plastic deformation amount, and calculating the acoustic elasticity parameter in the metal plate Sound elasticity A parameter calculation unit, a acoustoelastic parameter coefficient correction unit that corrects a coefficient representing a change due to stress in the acoustoelastic parameter using the plastic deformation feature value calculated by the plastic deformation feature value calculation unit, and A residual stress calculating unit that calculates a residual stress on the metal plate based on the acoustoelastic parameter calculated by the acoustoelastic parameter calculating unit and the coefficient corrected by the acoustoelastic parameter coefficient correcting unit. A stress measurement device is provided.
また、上記課題を解決するために、本発明の更に別の観点によれば、塑性変形した金属板中に複数の超音波を発生させることで測定された金属板中の複数の超音波音速に関する音速情報を取得するステップと、前記金属板の塑性変形量又は塑性変形量に換算可能な値である塑性変形特徴量を算出するとともに、取得した前記音速情報に基づいて、前記金属板中の音弾性パラメータを算出するステップと、算出した前記塑性変形特徴量を利用して、前記音弾性パラメータの応力による変化を表す係数の補正を行うステップと、算出された音弾性パラメータと、補正された係数とに基づいて、前記金属板に残留している応力を算出するステップと、を含む残留応力算出方法が提供される。 In order to solve the above problems, according to still another aspect of the present invention, the present invention relates to a plurality of ultrasonic sound speeds in a metal plate measured by generating a plurality of ultrasonic waves in a plastically deformed metal plate. Obtaining a sound velocity information; calculating a plastic deformation amount of the metal plate or a plastic deformation feature amount that is a value that can be converted into a plastic deformation amount; and, based on the obtained sound velocity information, a sound in the metal plate. A step of calculating an elastic parameter; a step of correcting a coefficient representing a change of the acoustoelastic parameter due to stress using the calculated plastic deformation feature; a calculated acoustoelastic parameter; and a corrected coefficient And calculating a residual stress on the metal plate, based on the above.
また、上記課題を解決するために、本発明の更に別の観点によれば、塑性変形した金属板中に複数の超音波を発生させることで、当該金属板中を伝播する前記複数の超音波の音速を測定するステップと、前記金属板中を伝播する複数の超音波音速に関する測定結果である音速情報を取得するステップと、前記金属板の塑性変形量又は塑性変形量に換算可能な値である塑性変形特徴量を算出するとともに、取得した前記音速情報に基づいて、前記金属板中の音弾性パラメータを算出するステップと、算出した前記塑性変形特徴量を利用して、前記音弾性パラメータの応力による変化を表す係数の補正を行うステップと、算出された音弾性パラメータと、補正された係数とに基づいて、前記金属板に残留している応力を算出するステップと、を含む残留応力測定方法が提供される。 In order to solve the above problems, according to still another aspect of the present invention, the plurality of ultrasonic waves propagating in the metal plate by generating the plurality of ultrasonic waves in the plastically deformed metal plate. A step of measuring the sound velocity of the metal plate, a step of obtaining sound velocity information that is a measurement result of a plurality of ultrasonic sound velocities propagating in the metal plate, and a plastic deformation amount of the metal plate or a value that can be converted into a plastic deformation amount. While calculating a certain plastic deformation feature amount, calculating the acoustoelastic parameter in the metal plate based on the acquired sound velocity information, and using the calculated plastic deformation feature amount, Correcting a coefficient representing a change due to stress, and calculating a residual stress on the metal plate based on the calculated acoustoelastic parameter and the corrected coefficient. Residual stress measurement method is provided.
また、上記課題を解決するために、本発明の更に別の観点によれば、コンピュータに、塑性変形した金属板中に複数の超音波を発生させることで測定された金属板中の複数の超音波音速に関する音速情報を取得する音速情報取得機能と、前記金属板の塑性変形量又は塑性変形量に換算可能な値である塑性変形特徴量を算出する塑性変形特徴量算出機能と、前記音速情報取得機能が取得した前記音速情報に基づいて、前記金属板中の音弾性パラメータを算出する音弾性パラメータ算出機能と、前記塑性変形特徴量算出機能が算出した前記塑性変形特徴量を利用して、前記音弾性パラメータの応力による変化を表す係数の補正を行う音弾性パラメータ係数補正機能と、前記音弾性パラメータ算出機能が算出した音弾性パラメータと、前記音弾性パラメータ係数補正機能が補正した係数とに基づいて、前記金属板に残留している応力を算出する残留応力算出機能と、を実現させるためのプログラムが提供される。 In order to solve the above problems, according to still another aspect of the present invention, a plurality of ultrasonic waves in a metal plate measured by causing a computer to generate a plurality of ultrasonic waves in a plastically deformed metal plate. A sound speed information acquisition function for acquiring sound speed information related to the sound wave speed, a plastic deformation feature amount calculation function for calculating a plastic deformation feature amount that is a value that can be converted into a plastic deformation amount or a plastic deformation amount of the metal plate, and the sound speed information. Based on the sound velocity information acquired by the acquisition function, using the acoustoelastic parameter calculation function for calculating the acoustoelastic parameter in the metal plate, and the plastic deformation feature amount calculated by the plastic deformation feature amount calculation function, A acoustoelastic parameter coefficient correction function for correcting a coefficient representing a change due to stress of the acoustoelastic parameter, an acoustoelastic parameter calculated by the acoustoelastic parameter calculation function, and the acoustoelastic parameter. Based on the coefficient meter factor correction function is corrected, program for implementing a residual stress calculating function of calculating a stress remaining in the metal plate is provided.
以上説明したように本発明によれば、金属板内部を伝播する音波の音速の関係から算出した音弾性パラメータと、塑性変形が音弾性パラメータに及ぼす影響を補正するための補正値とを利用して、金属板中の残留応力を算出するため、塑性変形している被検査体の内部に残留している応力を、安全かつ正確に測定することが可能となる。 As described above, according to the present invention, the acoustoelastic parameter calculated from the relationship of the sound velocity of the sound wave propagating through the metal plate and the correction value for correcting the influence of plastic deformation on the acoustoelastic parameter are used. Thus, since the residual stress in the metal plate is calculated, it is possible to safely and accurately measure the stress remaining in the plastically deformed test object.
以下に添付図面を参照しながら、本発明の好適な実施の形態について詳細に説明する。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。 Exemplary embodiments of the present invention will be described below in detail with reference to the accompanying drawings. In addition, in this specification and drawing, about the component which has the substantially same function structure, duplication description is abbreviate | omitted by attaching | subjecting the same code | symbol.
（目的）
まず、本発明の実施形態について説明するに先立ち、図１９〜図２５を参照しながら従来の音弾性効果を利用した残留応力の測定方法について説明し、従来の残留応力測定方法の問題点を明らかにする。また、この従来の残留応力測定方法の問題点を踏まえたうえでの本発明が目的とするところについて、説明を行う。図１９は、従来の音弾性効果を利用した残留応力測定方法の流れを示した流れ図である。図２０及び図２１は、塑性歪がないときの音弾性パラメータの測定結果を示したグラフ図である。図２２及び図２３は、塑性歪があるときの音弾性パラメータの測定結果を示したグラフ図である。図２４及び図２５は、塑性歪を考慮しないときの応力測定結果を示したグラフ図である。
(the purpose)
First, prior to describing the embodiment of the present invention, a conventional method for measuring residual stress using the acoustoelastic effect will be described with reference to FIGS. 19 to 25, and the problems of the conventional residual stress measuring method will be clarified. To. Further, the purpose of the present invention based on the problems of the conventional residual stress measurement method will be described. FIG. 19 is a flowchart showing the flow of a conventional residual stress measurement method using the acoustoelastic effect. 20 and 21 are graphs showing measurement results of acoustoelastic parameters when there is no plastic strain. 22 and 23 are graphs showing measurement results of acoustoelastic parameters when there is plastic strain. 24 and 25 are graphs showing stress measurement results when plastic strain is not considered.
まず、図１９を参照しながら、従来の音弾性（ａｃｏｕｓｔｏｅｌａｓｔｉｃｉｔｙ）効果を利用した残留応力測定方法の流れを簡単に説明する。 First, with reference to FIG. 19, the flow of a conventional residual stress measurement method using the acoustic elasticity effect will be briefly described.
従来の音弾性効果を利用した残留応力測定方法では、まず、所定の音速測定装置を利用して、被検査体（例えば、金属板等）の内部を伝播する音波（超音波）の音速を測定する（ステップＳ１１）。具体的には、音速測定装置を利用して、被検査体の内部に縦波超音波及び横波超音波を発生させて、多重反射エコー間隔を測定したり、共振周波数を測定したりする。これにより、縦波超音波及び横波超音波について、例えば以下のような物理量を測定することができる。 In the conventional residual stress measurement method using the acoustoelastic effect, first, the sound velocity of the sound wave (ultrasonic wave) propagating inside the object to be inspected (for example, a metal plate) is measured using a predetermined sound velocity measuring device. (Step S11). Specifically, longitudinal sound waves and transverse wave ultrasonic waves are generated inside the object to be inspected using a sound velocity measuring device, and the multiple reflection echo interval is measured or the resonance frequency is measured. Thereby, for example, the following physical quantities can be measured for the longitudinal ultrasonic waves and the transverse ultrasonic waves.
・縦波について
エコー間隔 ＴＬ
共振周波数 ｆＬｎ（ｎは、次数である。）
・横波について
エコー間隔 ＴＳ１，ＴＳ２
共振周波数 ｆＳ１ｎ，ｆＳ２ｎ（ｎは、次数である。）
・ About longitudinal waves Echo interval TL
Resonant frequency f Ln (n is the order)
・ About transverse waves Echo interval T S1 , T S2
Resonant frequencies f S1n , f S2n (n is the order)
ここで、音速とエコー間隔の間、及び、音速と共振周波数の間には、超音波の伝播距離（例えば、被検査体の厚み）ｔを用いて、以下の式１１及び式１２で表される関係が成立する。なお、以下の式１１及び式１２において、Ｃは音速を表す。また、ｎは、共振周波数の次数である。 Here, between the sound velocity and the echo interval, and between the sound velocity and the resonance frequency, the propagation distance of the ultrasonic wave (for example, the thickness of the object to be inspected) t is used, and is expressed by the following equations 11 and 12. The relationship is established. In the following formulas 11 and 12, C represents the speed of sound. N is the order of the resonance frequency.
Ｃ＝２ｔ／Ｔ ・・・（式１１）
Ｃ＝２ｆｎｔ／ｎ ・・・（式１２）
C = 2t / T (Formula 11)
C = 2f n t / n (Expression 12)
ここで、上記式１１において、Ｔは、縦波のエコー間隔ＴＬ及び横波のエコー間隔ＴＳ１，ＴＳ２を指し、上記式１２において、ｆｎは、縦波の共振周波数ｆＬｎ及び横波の共振周波数ｆＳ１ｎ，ｆＳ２ｎを指すものである。 In Equation 11, T indicates the longitudinal wave echo interval T L and the transverse wave echo intervals T S1 and T S2 , and in Equation 12, f n represents the longitudinal wave resonance frequency f Ln and the transverse wave frequency. The resonance frequencies f S1n and f S2n are indicated.
よって、縦波超音波から得られた物理量（エコー間隔又は共振周波数）を利用することで、縦波超音波の音速ＣＬを算出することができ、横波超音波から得られた物理量を利用することで、横波超音波の音速ＣＳ１，ＣＳ２を算出することができる。 Thus, by utilizing the physical quantity obtained from the longitudinal ultrasonic waves (the echo interval or a resonance frequency), it is possible to calculate the sound velocity C L of longitudinal ultrasonic waves, using a physical quantity obtained from the transverse ultrasonic waves Thus, the sound speeds C S1 and C S2 of the transverse ultrasonic waves can be calculated.
続いて、従来の音弾性効果を利用した残留応力測定方法では、測定した音速を利用して、各種の音弾性パラメータを算出する（ステップＳ１３）。このような音弾性パラメータの例として、以下の式１３で表される音速比Ｒと、以下の式１４で表される音響複屈折Ｂと、を挙げることができる。 Subsequently, in the conventional residual stress measurement method using the acoustoelastic effect, various acoustoelastic parameters are calculated using the measured sound speed (step S13). Examples of such acoustoelastic parameters include a sound speed ratio R expressed by the following formula 13 and an acoustic birefringence B expressed by the following formula 14.
ここで、音弾性効果を利用した残留応力測定方法では、以下で説明するように音速を測定することが目的ではなく、上記式１３及び式１４に示したような音弾性パラメータを用いて、被検査体に残留している主応力を測定することが目的である。上記式１３及び式１４の内容から明らかなように、音弾性パラメータを表す式の分子及び分母それぞれには、音速が同一の次数で含まれているため、上記式１１及び式１２におけるｔ（超音波の伝播距離）は、分子と分母で打ち消しあうこととなる。従って、従来の音弾性効果を利用した残留応力測定方法では、エコー間隔や共振周波数等の測定値を利用するようにすれば、超音波の伝播距離（例えば、被検査体の厚み等）は、測定する必要はない。 Here, the residual stress measurement method using the acoustoelastic effect is not intended to measure the speed of sound, as will be described below, but using the acoustoelastic parameters as shown in the above equations 13 and 14, The purpose is to measure the main stress remaining in the test object. As is clear from the contents of Equation 13 and Equation 14 above, the numerator and denominator of the equation representing the acoustoelastic parameter include the same order of sound velocity. The propagation distance of the sound wave is canceled by the numerator and the denominator. Therefore, in the conventional residual stress measurement method using the acoustoelastic effect, if measurement values such as the echo interval and the resonance frequency are used, the propagation distance of the ultrasonic wave (for example, the thickness of the object to be inspected) is There is no need to measure.
続いて、従来の残留応力測定方法では、算出した音弾性パラメータを利用して、被検査体に残留している応力を算出する（ステップＳ１５）。 Subsequently, in the conventional residual stress measurement method, the stress remaining in the object to be inspected is calculated using the calculated acoustoelastic parameter (step S15).
ここで、音速比Ｒ及び音響複屈折Ｂは、被検査体に残留している主応力σ１及びσ２との間に以下の式１５及び式１６で表されるような関係が成立することが知られている。 Here, the relationship between the sound velocity ratio R and the acoustic birefringence B is expressed by the following equations 15 and 16 between the main stresses σ 1 and σ 2 remaining in the object to be inspected. It has been known.
ここで、上記式１５において、ＣＲは主応力和変化による音速比の変化を表す係数であり、定数項Ｒ０は無応力時の音速比である。また、上記式１６において、ＣＡは主応力差変化による音響複屈折の変化を表す係数であり、定数項Ｂ０は被検査体の無応力時の組織異方性を表す。 Here, in the above formula 15, C R is a coefficient representing a change in the sound velocity ratio by the main stress sum changes, the constant term R 0 is the speed of sound ratio under no stress. In the above equation 16, C A is a coefficient representing the change in acoustic birefringence due to the change in main stress difference, and the constant term B 0 represents the tissue anisotropy when the test object is not stressed.
そこで、ある被検査体について、音速比Ｒと主応力和（σ１＋σ２）との関係、及び、音響複屈折Ｂと主応力差（σ１−σ２）との関係を予め測定しておき、算出した音弾性パラメータから主応力和及び主応力差を算出して、被検査体に残留している主応力σ１及びσ２を算出する。 Therefore, the relationship between the sound speed ratio R and the principal stress sum (σ 1 + σ 2 ) and the relationship between the acoustic birefringence B and the principal stress difference (σ 1 −σ 2 ) are measured in advance for a certain object to be inspected. Then, the main stress sum and the main stress difference are calculated from the calculated acoustoelastic parameters, and the main stresses σ 1 and σ 2 remaining in the object to be inspected are calculated.
ここで、被検査体として、引張強度の異なる２種類の鋼板（サンプルＡ：板厚１ｍｍの引張強度９８０ＭＰａ級高張力鋼板、サンプルＢ：板厚１．４ｍｍの引張強度５９０ＭＰａ級高張力鋼板）を用意し、引張試験機により荷重を加えながら音弾性パラメータを測定した。準備した２種類の鋼板には、塑性変形は生じていない。得られた測定結果を、図２０及び図２１に示す。 Here, two types of steel plates having different tensile strengths (Sample A: tensile strength 980 MPa class high strength steel plate with a thickness of 1 mm, sample B: tensile strength 590 MPa class high strength steel plate with a thickness of 1.4 mm) are used as test objects. The acoustoelastic parameter was measured while applying a load with a tensile tester. The two types of prepared steel plates are not plastically deformed. The obtained measurement results are shown in FIGS.
図２０及び図２１において、上側のグラフ図が音速比Ｒと応力（すなわち、鋼板に加えた荷重の大きさであり、主応力和に相当するもの）との関係を示したグラフ図である。また、下側のグラフ図が音響複屈折Ｂと応力（すなわち、鋼板に加えた荷重の大きさであり、主応力差に相当するもの）との関係を示したグラフ図である。 20 and 21, the upper graph is a graph showing the relationship between the sound speed ratio R and the stress (that is, the magnitude of the load applied to the steel sheet, which corresponds to the principal stress sum). Further, the lower graph is a graph showing the relationship between the acoustic birefringence B and the stress (that is, the magnitude of the load applied to the steel sheet and corresponding to the main stress difference).
図２０及び図２１から明らかなように、被検査体である鋼板に塑性変形が生じていない場合、加えた荷重（応力）と音弾性パラメータとの間には、線形関係があることがわかる。従って、この線形関係を利用することで、算出した音弾性パラメータから主応力和及び主応力差を算出することができるため、被検査体に残留している主応力σ１及びσ２を特定することができる。 As apparent from FIGS. 20 and 21, it is understood that there is a linear relationship between the applied load (stress) and the acoustoelastic parameter when plastic deformation does not occur in the steel plate as the object to be inspected. Therefore, since the principal stress sum and the principal stress difference can be calculated from the calculated acoustoelastic parameters by using this linear relationship, the principal stresses σ 1 and σ 2 remaining in the inspection object are specified. be able to.
他方、被検査体に大きな歪が加えられており、被検査体が塑性歪を受けているときの音弾性パラメータの測定結果について説明する。上述のサンプルＡ及びサンプルＢに対して予め所定の塑性歪を与えておき、これら塑性歪が与えられたサンプルに対して、引張試験機により引張応力を加えながら、音弾性パラメータを測定した。 On the other hand, the measurement result of the acoustoelastic parameter when a large strain is applied to the inspection object and the inspection object is subjected to plastic strain will be described. Predetermined plastic strain was given to the above-mentioned sample A and sample B, and the acoustoelastic parameters were measured while applying tensile stress to the sample given these plastic strains with a tensile tester.
得られた結果を、図２２及び図２３に示した。
なお、鋼板に対して与えた塑性歪の大きさは、鋼板が伸びた割合で示している。例えば、各グラフ図において２％という表記が意味するものは、鋼板に引張応力を加えることで鋼板の長さが２％伸びる程度の引張応力を加えたということである。各グラフ図から明らかなように、塑性歪の存在により、グラフ図の傾きや切片の位置が変化している。グラフ図の傾きや切片の位置が異なるということは、塑性歪が生じることで同一応力条件下であっても音弾性パラメータが変化しているということを示すものである。
The obtained results are shown in FIG. 22 and FIG.
In addition, the magnitude | size of the plastic strain given with respect to the steel plate is shown by the ratio which the steel plate extended. For example, the notation of 2% in each graph means that a tensile stress that increases the length of the steel sheet by 2% is applied by applying a tensile stress to the steel sheet. As is apparent from each graph, the inclination of the graph and the position of the intercept change due to the presence of plastic strain. The fact that the inclination of the graph and the position of the intercept are different indicate that the acoustoelastic parameter is changed even under the same stress condition due to the occurrence of plastic strain.
図２０及び図２１に示した、塑性歪が存在しないときの音弾性パラメータと応力との関係を利用して、図２２及び図２３に示した音弾性パラメータに基づいて残留応力を算出した結果を、図２４及び図２５に示す。図２４が、サンプルＡについての残留応力測定結果であり、図２５が、サンプルＢについての残留応力測定結果である。 The results of calculating the residual stress based on the acoustoelastic parameters shown in FIGS. 22 and 23 using the relationship between the acoustoelastic parameters and the stress when no plastic strain exists, as shown in FIGS. 24 and 25. FIG. 24 shows the residual stress measurement results for sample A, and FIG. 25 shows the residual stress measurement results for sample B.
行った測定は、被検査体である鋼板に対して一軸引張試験を行った結果に基づくものである。そのため、引っ張り方向と平行な方向の主応力σ１については引張応力に応じて増加する線形関係が成り立つとともに、引っ張り方向に直交する方向の主応力σ２については、引張応力によらず０近傍で一定の値となることが予想される。図２４及び図２５をみると、塑性歪が大きくなるにつれて、塑性歪がない場合（０％）に比べて測定結果の誤差が大きくなっていることがわかる。 The measurement performed is based on the result of a uniaxial tensile test performed on a steel plate as an object to be inspected. Therefore, a linear relationship that increases according to the tensile stress is established for the principal stress σ 1 in the direction parallel to the tensile direction, and the principal stress σ 2 in the direction orthogonal to the tensile direction is near 0 regardless of the tensile stress. A constant value is expected. 24 and 25, it can be seen that as the plastic strain increases, the error in the measurement result increases as compared to the case where there is no plastic strain (0%).
このように、従来の音弾性効果を利用した残留応力測定方法では、弾性変形領域ではなく塑性変形領域に相当する歪みが被検査体に加えられた場合には誤差が大きくなって、被検査体に残留している応力の大きさを正確に測定できないという問題があった。 As described above, in the conventional residual stress measurement method using the acoustoelastic effect, when a strain corresponding to the plastic deformation region, not the elastic deformation region, is applied to the inspection object, the error becomes large, and the inspection object There is a problem that the magnitude of the residual stress cannot be measured accurately.
そこで、本発明者らは、塑性変形した被検査体であっても、被検査体に残留している応力の大きさを正確に算出することが可能な残留応力の算出方法を見いだすことを目的として、鋭意検討を行った。その結果、以下で説明するような、本発明の実施形態に係る残留応力の算出方法及び残留応力の測定方法に想到した。 Therefore, the present inventors have aimed to find a method for calculating a residual stress that can accurately calculate the magnitude of the stress remaining in the object to be inspected even if the object is plastically deformed. As a result, we conducted intensive studies. As a result, the inventors have arrived at a residual stress calculation method and a residual stress measurement method according to an embodiment of the present invention as described below.
（第１の実施形態）
＜残留応力測定装置の構成について＞
まず、図１を参照しながら、本発明の第１の実施形態に係る残留応力測定装置の構成について、説明を行う。図１は、本実施形態に係る残留応力測定装置の構成を示した説明図である。
(First embodiment)
<About the configuration of the residual stress measurement device>
First, the configuration of the residual stress measurement apparatus according to the first embodiment of the present invention will be described with reference to FIG. FIG. 1 is an explanatory view showing a configuration of a residual stress measuring apparatus according to this embodiment.
本実施形態に係る残留応力測定装置１０は、被検査体である金属板サンプルＳに残留する応力を測定する装置である。この残留応力測定装置１０は、図１に例示したように、音速測定装置１００と、残留応力算出装置２００と、を主に備える。
The residual
また、本実施形態に係る残留応力測定装置１０は、図１に示した装置のほかに、表面粗度を測定する装置や、歪ゲージなどの歪みを測定する装置や、被検査体である金属板の板厚を測定する各種の装置等を備えていてもよい。
In addition to the apparatus shown in FIG. 1, the residual
［音速測定装置について］
音速測定装置１００は、被検査体である金属板サンプルＳに超音波を発生させて、金属板Ｓ中を伝播する超音波の音速を測定する装置である。この音速測定装置１００は、シングアラウンド法、パルスエコーオーバーラップ法、パルス重畳法、共振法、時間間隔平均法等といったものであってもよい。以下では、例えば図２に示したような、多重反射エコー間隔を利用した音速測定装置１００、又は、共振周波数を利用した音速測定装置１００を使用する場合について説明するが、本発明に係る音速測定装置１００が図２に示した例に限定されるわけではない。
[Sonic velocity measuring device]
The sound
図２は、本実施形態に係る音速測定装置１００の構成を示した説明図である。
まず、図２（ａ）を参照しながら、多重反射エコー間隔を利用した音速測定装置の構成について、簡単に説明する。
FIG. 2 is an explanatory diagram showing the configuration of the sound
First, with reference to FIG. 2 (a), the configuration of a sound speed measuring apparatus using multiple reflection echo intervals will be briefly described.
多重反射エコー間隔を利用した音速測定装置１００は、図２（ａ）に示したように、パルサーレシーバー１０１と、プローブ１０３と、を主に備える。
As shown in FIG. 2A, the sound
パルサーレシーバー１０１は、いわゆる超音波受発信装置であり、後述するプローブ１０３を介して、発生させた超音波を被検査体である金属板Ｓに発信するとともに、金属板Ｓ内で反射した反射波を受信する装置である。
The
また、プローブ１０３は、パルサーレシーバー１０１が発生させた超音波を金属板Ｓに送信するとともに、金属板Ｓ内で反射した反射波を受信する。
The
多重反射エコー間隔を利用した音速測定法では、図２（ａ）右側に示したグラフ図のように、第１反射波、第２反射波・・・といった多重反射波を観測し、この反射波の間隔（エコー間隔、図２（ａ）のグラフ図における時間Ｔ）を測定する。多重反射波のエコー間隔は、金属板Ｓの板厚ｔをパラメータとすると、金属板中の音速Ｃとの間に、上記式１１に示したような関係を有するため、この多重反射波のエコー間隔Ｔから金属板Ｓ中の音速を算出することが可能である。 In the sound velocity measurement method using the multiple reflection echo interval, multiple reflected waves such as a first reflected wave, a second reflected wave,... Are observed as shown in the graph on the right side of FIG. Interval (echo interval, time T in the graph of FIG. 2A) is measured. The echo interval of the multiple reflected waves has a relationship as shown in the above equation 11 with the sound velocity C in the metal plate, where the thickness t of the metal plate S is a parameter. The speed of sound in the metal plate S can be calculated from the interval T.
次に、図２（ｂ）を参照しながら、共振周波数を利用した音速測定装置の構成について、簡単に説明する。 Next, with reference to FIG. 2 (b), the configuration of the sound speed measuring device using the resonance frequency will be briefly described.
共振周波数を利用した音速測定装置１００は、図２（ｂ）に示したように、正弦波発生装置１１１と、バースト波発生装置１１３と、電磁超音波変換素子（Ｅｌｅｃｔｒｏ−Ｍａｇｎｅｔｉｃ Ａｃｏｕｓｔｉｃ Ｔｒａｎｓｄｕｃｅｒ：ＥＭＡＴ）１１５と、積分器１１７と、を主に備える。
As shown in FIG. 2B, the sound
正弦波発生装置１１１は、任意の周波数を有する正弦波を発生させる装置である。正弦波発生装置１１１が発生させた正弦波は、後述するバースト波発生装置１１３がバースト波を発生させる際の搬送波として利用される。
The
バースト波発生装置１１３は、正弦波発生装置１１１が発生させた正弦波を搬送波として利用し、バースト波を発生させる装置である。バースト波発生装置１１３が発生させたバースト波は、後述するＥＭＡＴ１１５を介して金属板Ｓに発信される。また、バースト波発生装置１１３は、ＥＭＡＴ１１５を介して受信した金属板Ｓ中で反射された反射波を、後述する積分器１１７に伝送する。
The
電磁超音波変換素子（ＥＭＡＴ）１１５は、バースト波発生装置１１３が発生させたバースト波を振動（すなわち、超音波）に変換して、被検査体である金属板Ｓに伝達させる。また、ＥＭＡＴ１１５は、金属板Ｓ中で反射された反射波を受信して電気信号へと変換し、バースト波発生装置１１３に伝送する。このＥＭＡＴ１１５を利用することで、水や油などの接触媒質が不要となり、金属板Ｓ中を伝播する超音波の音速を非接触で測定することが可能となる。
The electromagnetic ultrasonic transducer (EMAT) 115 converts the burst wave generated by the
積分器１１７は、ＥＭＡＴ１１５が受信した超音波（金属板Ｓ中で反射された反射波）に対応する電気信号を積分する装置である。この積分器１１７により、金属板Ｓ中を伝播した超音波の共振周波数を特定することが可能となる。
The
共振周波数を利用した音速測定法では、図２（ｂ）右側に示したグラフ図のように、共振周波数ｆ０，ｆ１，ｆ２，・・・を観測し、これらの共振周波数ｆｎを測定する。共振周波数は、金属板Ｓの板厚ｔをパラメータとすると、金属板中の音速Ｃとの間に、上記式１２に示したような関係を有するため、この共振周波数から金属板Ｓ中の音速を算出することが可能である。 The speed-of-sound measurement method using the resonance frequency, as shown in FIG. 2 (b) graph shown on the right side, the resonance frequency f 0, f 1, f 2 , observes ..., these resonant frequencies f n taking measurement. The resonance frequency has a relationship as shown in the above equation 12 with the sound speed C in the metal plate when the plate thickness t of the metal plate S is used as a parameter. Therefore, the sound speed in the metal plate S is calculated from this resonance frequency. Can be calculated.
なお、図２に示した音速測定装置１００は、被検査体である金属板Ｓに対して縦波超音波及び横波超音波を発生させることが好ましいが、音速測定装置１００が発生させる超音波の種類は、上述の例に限定されるわけではない。音速測定装置１００は、ＳＨ波などの各種の超音波を発生させて、被検査体中を伝播する超音波の音速を測定してもよい。
Note that the sonic
例えば、表面ＳＨ波を利用して音速を測定する場合には、上記式１５と組み合わせる式として、上記式１６のかわりに以下の式１７に示す音響異方性ΦＳを用いて、音速を測定することが可能である。 For example, when the sound speed is measured using surface SH waves, the sound speed is measured using acoustic anisotropy ΦS shown in the following expression 17 instead of the above expression 16 as an expression combined with the above expression 15. It is possible.
ここで、上記式１７において、ＶＳ１、ＶＳ２は主応力方向に偏向した表面ＳＨ波の音速であり、ＣＳは、主応力差変化による音響異方性の変化を表す係数である。 Here, in Expression 17, V S1 and V S2 are the sound speeds of the surface SH waves deflected in the main stress direction, and C S is a coefficient representing a change in acoustic anisotropy due to a change in main stress difference.
このような音速測定装置１００で測定された金属板Ｓ中を伝播する超音波の音速に関する測定結果は、音速に関する音速情報として、後述する残留応力算出装置２００に通知される。
The measurement result relating to the sound velocity of the ultrasonic wave propagating through the metal plate S measured by such a sound
以上、本実施形態に係る音速測定装置１００について説明した。
続いて、図３を参照しながら、本実施形態に係る残留応力測定装置１０が備える残留応力算出装置２００の構成について説明する。
The sound
Next, the configuration of the residual
［残留応力算出装置について］
図３は、本実施形態に係る残留応力算出装置２００の構成を示したブロック図である。本実施形態に係る残留応力算出装置２００は、図３に示したように、音速情報取得部２０１と、音弾性パラメータ算出部２０３と、塑性変形特徴量算出部２０５と、音弾性パラメータ係数補正部２０７と、記憶部２０９と、残留応力算出部２１１と、を主に備える。また、残留応力算出装置２００は、更に、表示制御部２１３と、表示部２１５とを備えていてもよい。
[Residual stress calculation device]
FIG. 3 is a block diagram showing the configuration of the residual
音速情報取得部２０１は、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、ＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）、通信装置等により実現される。音速情報取得部２０１は、音速測定装置１００が塑性変形した金属板中に超音波を発生させることで測定した金属板中の音速に関する情報（音速情報）を取得する。この音速情報は、縦波超音波を用いて測定された縦波の音速に関する情報と、２種類の横波超音波を用いて測定された横波の音速に関する情報と、を含むことが好ましい。この音速情報は、音速そのものを表す情報であってもよく、エコー間隔や共振周波数などの音速を算出するために利用可能な情報であってもよい。音速情報取得部２０１は、取得した音速情報（すなわち、音速測定値に関する情報）を、後述する音弾性パラメータ算出部２０３に出力する。
The sound speed
音弾性パラメータ算出部２０３は、例えば、ＣＰＵ、ＲＯＭ、ＲＡＭ等により実現される。音弾性パラメータ算出部２０３は、音速情報取得部２０１から通知された音速情報に基づいて、被検査体である金属板Ｓに関する音弾性パラメータを算出する。具体的には、音弾性パラメータ算出部２０３は、音速情報取得部２０１から通知された音速情報に基づいて、金属板中に残留する主応力の和又は差と相関のある物理量を、音弾性パラメータとして算出する。このような、主応力の和又は差と相関のある物理量の例として、例えば、音速比Ｒ及び音響複屈折Ｂを挙げることができる。なお、以下では、音弾性パラメータ算出部２０３が、音弾性パラメータとして、音速比Ｒ及び音響複屈折Ｂを算出する場合について、説明を行うものとする。
The acoustic elasticity
音弾性パラメータ算出部２０３は、音速情報取得部２０１から通知された音速情報を利用して、下記式１０１に基づいて音速比Ｒを算出する。また、音弾性パラメータ算出部２０３は、音速情報取得部２０１から通知された音速情報を利用して、下記式１０２に基づいて音響複屈折Ｂを算出する。音弾性パラメータ算出部２０３は、これら音速比Ｒ及び音響複屈折Ｂの算出が終了すると、算出したこれらの音弾性パラメータを、後述する残留応力算出部２１１に出力する。
The sound elasticity
塑性変形特徴量算出部２０５は、例えば、ＣＰＵ、ＲＯＭ、ＲＡＭ、通信装置等により実現される。塑性変形特徴量算出部２０５は、金属板の塑性変形の大きさを表す塑性変形特徴量を算出する。この塑性変形特徴量は、歪ゲージ等で測定された塑性変形量そのものであってもよく、塑性変形量に換算可能な、変形前後の板厚から算出される減肉比であってもよい。ここで、変形前後の板厚から算出される減肉比は、以下の式により定義される値である。また、塑性変形量に換算可能な値としては、他に表面粗度や硬度、ＦＥＭ解析結果から推定した塑性変形量等がある。
The plastic deformation feature
減肉比＝（ｔ−ｔ０）／ｔ０ Thinning ratio = (t−t 0 ) / t 0
ここで、上記減肉比の定義式において、ｔ０は、塑性変形が加えられる前の金属板の板厚を示しており、ｔは、塑性変形が加えられた後の金属板の板厚を示している。 Here, in the definition formula of the thickness reduction ratio, t 0 indicates the plate thickness of the metal plate before plastic deformation is applied, and t indicates the plate thickness of the metal plate after plastic deformation is applied. Show.
塑性変形特徴量算出部２０５は、取得した塑性変形量そのもの、又は、塑性変形量に換算可能な値である塑性変形特徴量を、後述する音弾性パラメータ係数補正部２０７に出力する。
The plastic deformation feature
なお、塑性変形特徴量算出部２０５は、塑性変形特徴量の算出に際して、後述する記憶部２０９等に格納されている各種の検量線、近似式、ルックアップテーブル、データベース等を参照してもよい。
The plastic deformation feature
音弾性パラメータ係数補正部２０７は、例えば、ＣＰＵ、ＲＯＭ、ＲＡＭ等により実現される。音弾性パラメータ係数補正部２０７は、前記塑性変形特徴量算出部により算出された前記塑性変形特徴量を利用して、前記音弾性パラメータの応力による変化を表す係数の補正を行う。以下、より詳細に、音弾性パラメータ係数補正部２０７が実施する係数の補正処理について説明する。
The acoustoelastic parameter
先に挙げた式１５及び式１６に示したように、音弾性パラメータ算出部２０３によって算出される音弾性パラメータは、主応力の和及び差を用いて定式化される量である。そこで、本実施形態に係る残留応力算出装置２００では、塑性変形による音弾性パラメータの変化を補正するために、式１５及び式１６等の理論式に含まれる係数の値を補正する。
As shown in Equations 15 and 16, the acoustoelastic parameter calculated by the acoustoelastic
本実施形態に係る残留応力算出装置２００では、塑性変形特徴量算出部２０５により算出される塑性変形特徴量と、式１５及び式１６等に記載されているような、理論式の係数との関係を表す相関情報を、予め生成しておく。
In the residual
例えば、音弾性パラメータ算出部２０３の算出する音弾性パラメータが音速比Ｒ及び音響複屈折Ｂであり、塑性変形特徴量算出部２０５の算出する塑性変形特徴量が塑性歪に関するものである場合を考える。この場合、以下に示す相関情報に対応する検量線、近似式、ルックアップテーブル、データベース等が予め作成され、残留応力算出装置２００の記憶部２０９等に格納される。
For example, consider a case where the acoustoelastic parameters calculated by the acoustoelastic
−塑性歪と式１５に示した係数ＣＲとの関係を表す相関情報
−塑性歪と式１５に示した係数Ｒ０との関係を表す相関情報
−塑性歪と式１６に示した係数ＣＡとの関係を表す相関情報
−塑性歪と式１６に示した係数Ｂ０との関係を表す相関情報
- correlation information expressing the relationship between the coefficient C R shown in plastic strain and expression 15 - plastic strain and correlation information expressing the relationship between the coefficient R 0 shown in Formula 15 - coefficient shown in plastic strain and Formula 16 C A Information representing the relationship between the correlation strain and the correlation information representing the relationship between the plastic strain and the coefficient B 0 shown in Equation 16
音弾性パラメータ係数補正部２０７は、これらの相関情報に対応する各種の検量線等と、塑性変形特徴量算出部２０５から通知された塑性変形特徴量とを利用して、通知された塑性変形特徴量に対応する係数の値を算出する。これにより、式１５及び式１６に記載されている各係数が補正されたこととなる。
The acoustoelastic parameter
音弾性パラメータ係数補正部２０７は、補正後の係数を、後述する残留応力算出部２１１に出力する。
The acoustoelastic parameter
記憶部２０９は、本実施形態に係る残留応力算出装置２００が備えるストレージ装置の一例である。記憶部２０９には、音弾性パラメータ係数補正部２０７が音弾性パラメータの係数の補正に際して利用する、塑性変形特徴量と比例係数との相関関係及び塑性変形特徴量と定数項との相関関係を示す各種の検量線、近似式、ルックアップテーブル、データベース等が格納されている。また、記憶部２０９には、塑性変形特徴量算出部２０５が塑性変形特徴量を算出する際に利用する各種の検量線、近似式、ルックアップテーブル、データベース等が格納されていてもよい。更に、記憶部２０９には、本実施形態に係る残留応力算出装置２００が、何らかの処理を行う際に保存する必要が生じた様々なパラメータや処理の途中経過等、または、各種のデータベース等が、適宜記録される。この記憶部２０９は、残留応力算出装置２００が備える各処理部が、自由に読み書きを行うことが可能である。
The
残留応力算出部２１１は、例えば、ＣＰＵ、ＲＯＭ、ＲＡＭ等により実現される。残留応力算出部２１１は、音弾性パラメータ算出部２０３が算出した音弾性パラメータと、音弾性パラメータ係数補正部２０７が補正した係数とを利用して、被検査体である金属板Ｓに残留している主応力の大きさを算出する。以下、残留応力算出部２１１が実行する残留応力の算出手順について、詳細に説明する。
The residual
残留応力算出部２１１には、音弾性パラメータ算出部２０３が算出した、主応力の和を用いて定式化される音弾性パラメータと、主応力の差を用いて定式化される音弾性パラメータとが入力される。残留応力算出部２１１は、この２種類の音弾性パラメータの値を、音弾性パラメータ係数補正部２０７から通知された係数を利用して補正する。その後、残留応力算出部２１１は、２種類の補正後の音弾性パラメータを利用して、主応力σ１及び主応力σ２の大きさを算出する。
The residual
ここで、音弾性パラメータ算出部２０３により算出される音弾性パラメータが、音速比Ｒ（主応力の和を用いて定式化される音弾性パラメータ）及び音響複屈折Ｂ（主応力の差を用いて定式化される音弾性パラメータ）であるものとする。
Here, the acoustoelastic parameters calculated by the acoustoelastic
残留応力算出部２１１は、音弾性パラメータ算出部２０３より通知された音速比Ｒ及び音響複屈折Ｂを、音弾性パラメータ係数補正部２０７から通知された係数Ｒ０，ＣＲ，Ｂ０，ＣＡを用いて、以下に示す式１０３及び式１０４の右辺に示したように補正する。ここで、以下に示した式１０３及び式１０４は、音速比Ｒ及び音響複屈折Ｂの定義式（式１５及び式１６）を変形したものである。
The residual
補正の結果得られる値は、式１０３及び式１０４の左辺から明らかなように、主応力の和及び主応力の差に対応する値である。そこで、残留応力算出部２１１は、以下に示す式１０５及び式１０６と、補正の結果得られる値とを利用して、金属板Ｓに残留している主応力σ１及びσ２を算出する。
The values obtained as a result of the correction are values corresponding to the sum of the main stresses and the difference between the main stresses, as is apparent from the left side of the
このようにして算出される主応力σ１及びσ２は、塑性変形が音弾性パラメータに及ぼす影響を考慮した上で算出されたものであり、塑性変形が生じている被検査体であっても、従来の音弾性効果を利用した残留応力測定方法に比べて精度の良い結果を得ることができる。 The principal stresses σ 1 and σ 2 calculated in this way are calculated in consideration of the effect of plastic deformation on the acoustoelastic parameter, and even if the object to be inspected is undergoing plastic deformation. Thus, a more accurate result can be obtained as compared with the conventional residual stress measurement method using the acoustoelastic effect.
残留応力算出部２１１は、このようにして算出した被検査体に残留している残留応力の算出結果を、後述する表示制御部２１３に出力する。また、残留応力算出部２１１は、残留応力の算出結果を、残留応力を算出した日時に関する情報や、被検査体に関する情報等と関連付けて、履歴情報として記憶部２０９等に記録してもよい。また、残留応力算出部２１１は、算出した残留応力を、残留応力算出装置２００の外部に設けられた各種の装置に出力してもよい。この算出結果の出力は、所定の通信網を介して直接行われてもよく、各種の記録媒体を介して間接的に行われてもよい。
The residual
表示制御部２１３は、例えば、ＣＰＵ、ＲＯＭ、ＲＡＭ等により実現される。表示制御部２１３は、残留応力算出部２１１が算出した被検査体である金属板Ｓの残留応力の大きさを残留応力算出装置２００が備える表示部２１５に表示する際の表示制御を行う。また、表示制御部２１３は、残留応力の算出結果以外にも、算出した音弾性パラメータの値や、補正された係数の値や、塑性変形特徴量と音弾性パラメータの係数との間の相関関係を表すグラフ図等など、各種の情報を表示部２１５に表示させることができる。
The
表示部２１５は、例えば、ＣＰＵ、ＲＯＭ、ＲＡＭ、出力装置等により実現される。表示部２１５は、表示制御部２１３によって制御されながら、残留応力算出装置２００が備えるディスプレイ等の出力装置に、残留応力の算出結果を含む各種の情報を表示する。これにより、残留応力算出装置２００の利用者は、被検査体である金属板に残留している残留応力の大きさを、その場で把握することができる。
The
以上、本実施形態に係る残留応力算出装置２００の機能の一例を示した。上記の各構成要素は、汎用的な部材や回路を用いて構成されていてもよいし、各構成要素の機能に特化したハードウェアにより構成されていてもよい。また、各構成要素の機能を、ＣＰＵ等が全て行ってもよい。従って、本実施形態を実施する時々の技術レベルに応じて、適宜、利用する構成を変更することが可能である。
Heretofore, an example of the function of the residual
なお、上述のような本実施形態に係る残留応力算出装置の各機能を実現するためのコンピュータプログラムを作製し、パーソナルコンピュータ等に実装することが可能である。また、このようなコンピュータプログラムが格納された、コンピュータで読み取り可能な記録媒体も提供することができる。記録媒体は、例えば、磁気ディスク、光ディスク、光磁気ディスク、フラッシュメモリなどである。また、上記のコンピュータプログラムは、記録媒体を用いずに、例えばネットワークを介して配信してもよい。 A computer program for realizing each function of the residual stress calculation apparatus according to the present embodiment as described above can be produced and installed in a personal computer or the like. In addition, a computer-readable recording medium storing such a computer program can be provided. The recording medium is, for example, a magnetic disk, an optical disk, a magneto-optical disk, a flash memory, or the like. Further, the above computer program may be distributed via a network, for example, without using a recording medium.
＜残留応力算出方法の具体例について＞
続いて、本実施形態に係る残留応力算出装置２００が実施する残留応力算出方法の具体例について、図４Ａ〜図１６を参照しながら説明する。
<Specific example of residual stress calculation method>
Subsequently, a specific example of a residual stress calculation method performed by the residual
［塑性変形特徴量として塑性歪を利用する場合について］
まず、図４Ａ〜図９を参照しながら、塑性変形特徴量として塑性歪を利用する場合の例について、具体的に説明する。図４Ａ〜図５Ｂは、塑性歪と音弾性パラメータの係数との関係を示したグラフ図である。図６及び図７は、塑性歪による補正をしたときの応力測定結果を示したグラフ図である。図８及び図９は、塑性歪による補正をしたときの測定誤差の変化を示したグラフ図である。
[When plastic strain is used as a plastic deformation feature]
First, an example in which plastic strain is used as the plastic deformation feature amount will be specifically described with reference to FIGS. 4A to 9. 4A to 5B are graphs showing the relationship between the plastic strain and the coefficient of the acoustoelastic parameter. 6 and 7 are graphs showing stress measurement results when correction is made by plastic strain. 8 and 9 are graphs showing changes in measurement error when correction is made by plastic strain.
以下に示す具体例は、被検査体として引張強度の異なる２種類の鋼板（サンプルＡ：板厚１ｍｍの引張強度９８０ＭＰａ級高張力鋼板、サンプルＢ：板厚１．４ｍｍの引張強度５９０ＭＰａ級高張力鋼板）を用意して、測定を行った結果である。 Specific examples shown below are two types of steel plates having different tensile strengths as samples to be inspected (sample A: tensile strength 980 MPa class high strength steel plate with a thickness of 1 mm, sample B: tensile strength 590 MPa class high tension with a thickness of 1.4 mm). It is the result of having prepared and measured it.
まず、図４Ａ〜図５Ｂに着目する。図４Ａ及び図４Ｂは、サンプルＡについて、塑性歪と音弾性パラメータの係数との関係を示したグラフ図であり、図５Ａ及び図５Ｂは、サンプルＢについて、塑性歪と音弾性パラメータの係数との関係を示したグラフ図である。また、図４Ａ及び図５Ａは、各サンプルの音速比Ｒに関する測定結果を示しており、図４Ｂ及び図５Ｂは、各サンプルの音響複屈折Ｂに関する測定結果を示している。 First, attention is focused on FIGS. 4A to 5B. 4A and 4B are graphs showing the relationship between the plastic strain and the coefficient of the acoustoelastic parameter for the sample A. FIGS. 5A and 5B show the coefficient of the plastic strain and the acoustoelastic parameter for the sample B. It is the graph which showed this relationship. 4A and 5A show measurement results regarding the sound speed ratio R of each sample, and FIGS. 4B and 5B show measurement results regarding the acoustic birefringence B of each sample.
ここで、図４Ａ〜図５Ｂについて、横軸は、被検査体に予め与えられた塑性歪（％）を表しており、塑性歪の大きさは、鋼板が伸びた割合で示している。また、各グラフ図の縦軸は、係数の大きさを示している。各係数は、各塑性歪の値に対応する塑性変形を有するサンプルそれぞれについて引張試験機を用いた測定を行って、例えば図２０及び図２１に示したようなグラフ図を作成し、各グラフ図の傾き及び切片の値を求めることで算出した。また、各グラフ図には、塑性歪と係数との関係を表す近似直線と、その近似式が記してある。 Here, in FIG. 4A to FIG. 5B, the horizontal axis represents the plastic strain (%) given in advance to the object to be inspected, and the magnitude of the plastic strain is shown by the rate at which the steel plate is stretched. Moreover, the vertical axis | shaft of each graph figure has shown the magnitude | size of the coefficient. Each coefficient is measured by using a tensile tester for each sample having plastic deformation corresponding to each plastic strain value, and for example, graphs as shown in FIGS. 20 and 21 are created. It calculated by calculating | requiring the inclination and intercept value. Each graph shows an approximate straight line representing the relationship between the plastic strain and the coefficient, and an approximate expression thereof.
各グラフ図から明らかなように、係数と塑性歪との間には相関関係が認められ、両者の関係は、例えば線形近似が可能であることがわかった。そこで、各グラフ図から得られた近似式を、塑性歪と係数との相関関係を表す情報として利用し、サンプルＡ及びサンプルＢの残留応力を算出した。 As is clear from each graph, it was found that there was a correlation between the coefficient and the plastic strain, and the relationship between the two could be linearly approximated, for example. Therefore, the approximate expression obtained from each graph was used as information representing the correlation between the plastic strain and the coefficient, and the residual stress of Sample A and Sample B was calculated.
すなわち、塑性変形特徴量算出部２０５は、塑性変形が加えられる前後の金属板の長さの変化に関する情報を取得して長さの変化の割合を算出し、得られた値を塑性変形特徴量（すなわち、図４Ａ〜図５Ｂにおける塑性歪）とする。また、音弾性パラメータ係数補正部２０７は、該当するサンプルの各近似式に算出された塑性歪の大きさを代入して、各係数を決定する。その後、残留応力算出部２１１は、音弾性パラメータ係数補正部２０７が算出した係数と、音弾性パラメータ算出部２０３が算出した音弾性パラメータとを利用して、サンプルＡ及びサンプルＢの残留応力を算出した。
That is, the plastic deformation feature
得られた結果（残留応力の測定結果）を、図６及び図７に示す。図６及び図７において、横軸は応力（引張応力、単位ＭＰａ）を示し、縦軸は、主応力の大きさ（ＭＰａ）を示している。図６及び図７について、引張応力が加えられた方向に対して平行な方向の主応力σ１に着目すると、金属板に加えられた塑性変形の大きさが変化しても、得られる主応力σ１の大きさには、図２４及び図２５に示したような大きな変化は見られない。また、引っ張り応力が加えられた方向に直交する方向の主応力σ２に着目すると、図６及び図７では、値が０近傍でほぼ一定となっていることがわかる。 The obtained results (residual stress measurement results) are shown in FIGS. 6 and 7, the horizontal axis represents stress (tensile stress, unit MPa), and the vertical axis represents the magnitude of main stress (MPa). 6 and 7, when attention is paid to the principal stress σ 1 in a direction parallel to the direction in which the tensile stress is applied, the principal stress obtained even if the magnitude of the plastic deformation applied to the metal plate is changed. There is no significant change in the magnitude of σ 1 as shown in FIGS. Further, when attention is paid to the principal stress σ 2 in the direction orthogonal to the direction in which the tensile stress is applied, it can be seen that the values are almost constant in the vicinity of 0 in FIGS.
また、図８及び図９は、塑性歪による補正をしたときの測定誤差の変化の様子を示している。図８及び図９は、実際に加えられた引張応力の大きさと算出した主応力との差をプロットしたものである。図８及び図９から明らかなように、塑性歪による補正をしなかった場合（図８（ａ）及び図９（ａ））には、塑性歪の大きさが大きくなるにつれて誤差も大きくなっていることがわかる。また、塑性歪による補正を行った場合（図８（ｂ）及び図９（ｂ））には、塑性歪の大きさが大きくなっても、誤差は大きくなっておらず、また、誤差の絶対値も補正をしなかった場合に比べて小さいことがわかる。 8 and 9 show how the measurement error changes when correction is made by plastic strain. 8 and 9 plot the difference between the magnitude of the tensile stress actually applied and the calculated principal stress. As is apparent from FIGS. 8 and 9, when correction by plastic strain is not performed (FIGS. 8A and 9A), the error increases as the plastic strain increases. I understand that. In addition, when the correction by plastic strain is performed (FIGS. 8B and 9B), the error does not increase even when the plastic strain increases, and the absolute value of the error does not increase. It can be seen that the value is also smaller than when no correction is made.
このように、塑性歪による補正を行った上で残留応力の算出を行うことで、被検査体である金属板に塑性変形が加わっている場合であっても、残留応力を精度良く算出することが可能となる。 In this way, by calculating the residual stress after correcting by plastic strain, it is possible to accurately calculate the residual stress even when plastic deformation is applied to the metal plate that is the object to be inspected. Is possible.
［塑性変形特徴量として減肉比を利用する場合について］
続いて、図１０〜図１６を参照しながら、塑性変形特徴量として減肉比を利用する場合の例について、具体的に説明する。図１０は、塑性歪と減肉比との関係を示したグラフ図である。図１１Ａ〜図１２Ｂは、減肉比と音弾性パラメータの係数との関係を示したグラフ図である。図１３及び図１４は、減肉比による補正をしたときの応力測定結果を示したグラフ図である。図１５及び図１６は、減肉比による補正をしたときの測定誤差の変化を示したグラフ図である。
[When using the reduction ratio as a plastic deformation feature value]
Subsequently, an example in the case of using the thickness reduction ratio as the plastic deformation feature amount will be specifically described with reference to FIGS. FIG. 10 is a graph showing the relationship between the plastic strain and the thinning ratio. 11A to 12B are graphs showing the relationship between the thinning ratio and the coefficient of the acoustoelastic parameter. FIG. 13 and FIG. 14 are graphs showing the stress measurement results when correction is performed by the thinning ratio. FIG. 15 and FIG. 16 are graphs showing changes in measurement error when correction is performed using the thinning ratio.
以下では、塑性変形特徴量として、以下の式１０７で表される減肉比を利用した場合の例について、説明を行う。なお、以下の式１０７について、ｔ０は、塑性変形が加えられる前の金属板の板厚を示しており、ｔは、塑性変形が加えられた後の金属板の板厚を示している。下記式１０７からも明らかなように、減肉比は、塑性変形による板厚の変化の割合を示すパラメータである。 Below, the example at the time of utilizing the thickness reduction ratio represented by the following formula | equation 107 as a plastic deformation feature-value is demonstrated. In the following expression 107, t 0 indicates the plate thickness of the metal plate before plastic deformation is applied, and t indicates the plate thickness of the metal plate after plastic deformation is applied. As is clear from the following formula 107, the thickness reduction ratio is a parameter indicating the rate of change in plate thickness due to plastic deformation.
減肉比＝（ｔ−ｔ０）／ｔ０ ・・・（式１０７） Thinning ratio = (t−t 0 ) / t 0 (Equation 107)
ここで、被検査体である金属板の板厚は、マイクロメータ等による測定結果を利用することも可能であり、ＦＥＭ解析による推定結果を利用することも可能であるが、以下に示す例では、マイクロメータによる測定結果を利用した。 Here, the thickness of the metal plate as the object to be inspected can use the measurement result by a micrometer or the like, and can also use the estimation result by FEM analysis. In the example shown below, The measurement result with a micrometer was used.
図１０は、上述のサンプルＡ及びサンプルＢについて実際に塑性歪及び減肉比を算出し、その対応関係を示したものである。図１０から明らかなように、塑性歪と減肉比との間には相関関係が認められるため、塑性歪に換えて減肉比を利用することでも、被検査体に残留する残留応力を算出することが可能であると考えられる。 FIG. 10 shows the correlation between the plastic strain and the thinning ratio actually calculated for Sample A and Sample B described above. As is clear from FIG. 10, since a correlation is recognized between the plastic strain and the thinning ratio, the residual stress remaining in the object to be inspected can be calculated by using the thinning ratio instead of the plastic strain. It is considered possible to do.
そこで、サンプルＡ及びサンプルＢについて、塑性変形による減肉比と係数との関係を測定した。得られた結果が、図１１Ａ〜図１２Ｂである。図１１Ａ及び図１１Ｂは、サンプルＡについて、減肉比と音弾性パラメータの係数との関係を示したグラフ図であり、図１２Ａ及び図１２Ｂは、サンプルＢについて、減肉比と音弾性パラメータの係数との関係を示したグラフ図である。また、図１１Ａ及び図１２Ａは、各サンプルの音速比Ｒに関する測定結果を示しており、図１１Ｂ及び図１２Ｂは、各サンプルの音響複屈折Ｂに関する測定結果を示している。 Therefore, for sample A and sample B, the relationship between the thinning ratio due to plastic deformation and the coefficient was measured. The obtained result is FIG. 11A-FIG. 12B. 11A and 11B are graphs showing the relationship between the thinning ratio and the coefficient of acoustoelastic parameter for sample A. FIGS. 12A and 12B are the graphs of the thinning ratio and acoustoelastic parameter for sample B. It is the graph which showed the relationship with a coefficient. Moreover, FIG. 11A and FIG. 12A have shown the measurement result regarding the sound speed ratio R of each sample, and FIG. 11B and FIG. 12B have shown the measurement result regarding the acoustic birefringence B of each sample.
なお、図１１Ａ等において、「８．０Ｅ−０６」のような記載がなされているが、これは、「８．０×１０−６」という値に対応している。 In FIG. 11A and the like, a description such as “8.0E-06” is made, which corresponds to a value of “8.0 × 10 −6 ”.
各グラフ図から明らかなように、音弾性パラメータの係数と減肉比との間には相関関係が認められ、両者の関係は、例えば線形近似が可能であることがわかった。そこで、各グラフ図から得られた近似式を、減肉比と音弾性パラメータの係数との相関関係を表す情報として利用し、サンプルＡ及びサンプルＢの残留応力を算出した。 As is apparent from the graphs, a correlation was recognized between the coefficient of the acoustoelastic parameter and the thinning ratio, and it was found that the relationship between the two can be linearly approximated, for example. Accordingly, the approximate expression obtained from each graph was used as information indicating the correlation between the thinning ratio and the coefficient of the acoustoelastic parameter, and the residual stress of Sample A and Sample B was calculated.
得られた結果（残留応力の測定結果）を、図１３及び図１４に示す。図１３及び図１４において、横軸は応力（引張応力、単位ＭＰａ）を示し、縦軸は、主応力の大きさ（ＭＰａ）を示している。図１３及び図１４について、引張応力が加えられた方向に対して平行な方向の主応力σ１に着目すると、金属板に加えられた塑性変形の大きさが変化しても、得られる主応力σ１の大きさには、図２４及び図２５に示したような大きな変化は見られない。また、引っ張り応力が加えられた方向に直交する方向の主応力σ２に着目すると、図１３及び図１４では、値が０近傍でほぼ一定となっていることがわかる。 The obtained results (residual stress measurement results) are shown in FIGS. 13 and 14, the horizontal axis represents stress (tensile stress, unit MPa), and the vertical axis represents the magnitude of main stress (MPa). 13 and 14, when attention is paid to the principal stress σ 1 in the direction parallel to the direction in which the tensile stress is applied, the principal stress obtained even if the magnitude of the plastic deformation applied to the metal plate changes. There is no significant change in the magnitude of σ 1 as shown in FIGS. Further, when attention is paid to the main stress σ 2 in the direction orthogonal to the direction in which the tensile stress is applied, it can be seen that the values are almost constant in the vicinity of 0 in FIGS.
また、図１５及び図１６は、減肉比による補正をしたときの測定誤差の変化の様子を示している。図１５及び図１６から明らかなように、減肉比による補正をしなかった場合（図１５（ａ）及び図１６（ａ））には、減肉比の絶対値が大きくなるにつれて誤差も大きくなっていることがわかる。また、減肉比による補正を行った場合（図１５（ｂ）及び図１６（ｂ））には、減肉比の絶対値が大きくなっても、誤差は大きくなっておらず、また、誤差の絶対値も補正をしなかった場合に比べて小さいことがわかる。 FIG. 15 and FIG. 16 show how the measurement error changes when correction is performed using the thinning ratio. As is apparent from FIGS. 15 and 16, when the correction by the thinning ratio is not performed (FIGS. 15A and 16A), the error increases as the absolute value of the thinning ratio increases. You can see that Further, when the correction by the thinning ratio is performed (FIGS. 15B and 16B), the error does not increase even if the absolute value of the thinning ratio increases, and the error It can be seen that the absolute value of is smaller than the case where no correction is made.
このように、塑性歪そのものではなく塑性歪と相関のある特徴量を利用して補正を行った上で残留応力の算出を行うことで、被検査体である金属板に塑性変形が加わっている場合であっても、残留応力を精度良く算出することが可能となる。 In this way, the plastic deformation is applied to the metal plate, which is the object to be inspected, by calculating the residual stress after performing correction using the feature quantity correlated with the plastic strain, not the plastic strain itself. Even in this case, the residual stress can be calculated with high accuracy.
なお、図４Ａ〜図５Ｂ、及び、図１１Ａ〜図１２Ｂでは、音弾性パラメータの係数と塑性変形特徴量との間の相関関係を線形近似で近似しているが、２次関数、３次関数・・・等の多次関数、指数関数、対数関数等の各種の曲線で近似を行ってもよい。 In FIGS. 4A to 5B and FIGS. 11A to 12B, the correlation between the coefficient of the acoustoelastic parameter and the plastic deformation feature amount is approximated by linear approximation. Approximation may be performed using various curves such as multi-order functions such as... Exponential function, logarithmic function.
＜残留応力算出方法の流れについて＞
続いて、図１７を参照しながら、本実施形態に係る残留応力算出装置２００で実施される残留応力算出方法の流れについて、簡単に説明する。図１７は、本実施形態に係る残留応力算出方法の流れを示した流れ図である。
<Residual stress calculation method flow>
Subsequently, a flow of a residual stress calculation method performed by the residual
なお、以下の説明に先立ち、残留応力算出装置２００の記憶部２０９等には、塑性変形特徴量と音弾性パラメータの係数との相関関係を表す情報が予め格納されているものとする。
Prior to the following description, it is assumed that information indicating the correlation between the plastic deformation feature quantity and the coefficient of the acoustoelastic parameter is stored in advance in the
まず、残留応力算出装置２００の音速情報取得部２０１は、例えば音速測定装置１００等から、被検査体である金属板についての音速情報（音速の測定結果）を取得する（ステップＳ１０１）。音速情報取得部２０１は、被検査体についての音速情報を取得すると、取得した音速情報を、音弾性パラメータ算出部２０３に出力する。
First, the sound speed
次に、音弾性パラメータ算出部２０３は、音速情報取得部２０１から通知された音速情報を利用して、金属板の音弾性パラメータを算出する（ステップＳ１０３）。算出される音弾性パラメータは、例えば、主応力の和を用いて定式化される音弾性パラメータ、及び、主応力の差を用いて定式化される音弾性パラメータであり、このような音弾性パラメータの一例として、音速比Ｒ及び音響複屈折Ｂを挙げることができる。音弾性パラメータ算出部２０３は、算出した音弾性パラメータを、残留応力算出部２１１に出力する。
Next, the sound elasticity
他方、残留応力算出装置２００の塑性変形特徴量算出部２０５は、被検査体である金属板の塑性歪、板厚、長さ等の情報を取得して、塑性変形の度合いを示す特徴量である塑性変形特徴量を算出する（ステップＳ１０５）。このような塑性変形特徴量の一例として、塑性歪や減肉比等を挙げることができる。塑性変形特徴量算出部２０５は、算出した塑性変形特徴量を、音弾性パラメータ係数補正部２０７に出力する。
On the other hand, the plastic deformation feature
音弾性パラメータ係数補正部２０７は、塑性変形特徴量算出部２０５から通知された塑性変形特徴量と、記憶部２０９等に格納されている相関関係を表す情報とを利用して、音弾性パラメータの係数を補正する（ステップＳ１０７）。この音弾性パラメータの係数は、例えば、音弾性パラメータを主応力の和又は差を用いて定式化する際に式中に表れる係数である。音弾性パラメータ係数補正部２０７は、音弾性パラメータの係数を補正すると、補正した音弾性パラメータの係数を残留応力算出部２１１に出力する。
The acoustoelastic parameter
次に、残留応力算出部２１１は、音弾性パラメータ算出部２０３から通知された音弾性パラメータを、音弾性パラメータ係数補正部２０７から通知された係数を利用して補正する。これにより、塑性変形による影響を考慮した音弾性パラメータの値を得ることができる。その後、残留応力算出部２１１は、補正後の音弾性パラメータを利用して、被検査体である金属板に残留している残留応力を算出する（ステップＳ１０９）。
Next, the residual
以上説明したような流れで処理が行われることにより、本実施形態に係る残留応力算出方法では、被検査体である金属板が塑性変形している場合であっても、精度良く残留応力を算出することが可能となる。 By performing the process according to the flow described above, the residual stress calculation method according to the present embodiment calculates the residual stress with high accuracy even when the metal plate as the object to be inspected is plastically deformed. It becomes possible to do.
なお、図１７では、音弾性パラメータを算出するための一連の流れと、塑性変形特徴量を算出するための一連の流れとが平行して実施される場合について示しているが、音弾性パラメータの算出に関する流れと塑性変形特徴量の算出に関する流れとは、平行して行われなくとも良い。 FIG. 17 shows a case where a series of flows for calculating the acoustoelastic parameter and a series of flows for calculating the plastic deformation feature amount are performed in parallel. The flow relating to the calculation and the flow relating to the calculation of the plastic deformation feature amount do not have to be performed in parallel.
（ハードウェア構成について）
次に、図１８を参照しながら、本発明の実施形態に係る残留応力算出装置２００のハードウェア構成について、詳細に説明する。図１８は、本発明の実施形態に係る残留応力算出装置２００のハードウェア構成を説明するためのブロック図である。
(About hardware configuration)
Next, the hardware configuration of the residual
残留応力算出装置２００は、主に、ＣＰＵ９０１と、ＲＯＭ９０３と、ＲＡＭ９０５と、を備える。また、残留応力算出装置２００は、更に、バス９０７と、入力装置９０９と、出力装置９１１と、ストレージ装置９１３と、ドライブ９１５と、接続ポート９１７と、通信装置９１９とを備える。
The residual
ＣＰＵ９０１は、演算処理装置および制御装置として機能し、ＲＯＭ９０３、ＲＡＭ９０５、ストレージ装置９１３、またはリムーバブル記録媒体９２１に記録された各種プログラムに従って、残留応力算出装置２００内の動作全般またはその一部を制御する。ＲＯＭ９０３は、ＣＰＵ９０１が使用するプログラムや演算パラメータ等を記憶する。ＲＡＭ９０５は、ＣＰＵ９０１の実行において使用するプログラムや、その実行において適宜変化するパラメータ等を一次記憶する。これらはＣＰＵバス等の内部バスにより構成されるバス９０７により相互に接続されている。
The
バス９０７は、ブリッジを介して、ＰＣＩ（Ｐｅｒｉｐｈｅｒａｌ Ｃｏｍｐｏｎｅｎｔ Ｉｎｔｅｒｃｏｎｎｅｃｔ／Ｉｎｔｅｒｆａｃｅ）バスなどの外部バスに接続されている。
The
入力装置９０９は、例えば、マウス、キーボード、タッチパネル、ボタン、スイッチおよびレバーなどユーザが操作する操作手段である。また、入力装置９０９は、例えば、赤外線やその他の電波を利用したリモートコントロール手段（いわゆる、リモコン）であってもよいし、残留応力算出装置２００の操作に対応したＰＤＡ等の外部接続機器９２３であってもよい。さらに、入力装置９０９は、例えば、上記の操作手段を用いてユーザにより入力された情報に基づいて入力信号を生成し、ＣＰＵ９０１に出力する入力制御回路などから構成されている。残留応力算出装置２００のユーザは、この入力装置９０９を操作することにより、残留応力算出装置２００に対して各種のデータを入力したり処理動作を指示したりすることができる。
The
出力装置９１１は、取得した情報をユーザに対して視覚的または聴覚的に通知することが可能な装置で構成される。このような装置として、ＣＲＴディスプレイ装置、液晶ディスプレイ装置、プラズマディスプレイ装置、ＥＬディスプレイ装置およびランプなどの表示装置や、スピーカおよびヘッドホンなどの音声出力装置や、プリンタ装置、携帯電話、ファクシミリなどがある。出力装置９１１は、例えば、残留応力算出装置２００が行った各種処理により得られた結果を出力する。具体的には、表示装置は、残留応力算出装置２００が行った各種処理により得られた結果を、テキストまたはイメージで表示する。他方、音声出力装置は、再生された音声データや音響データ等からなるオーディオ信号をアナログ信号に変換して出力する。
The
ストレージ装置９１３は、残留応力算出装置２００の記憶部の一例として構成されたデータ格納用の装置である。ストレージ装置９１３は、例えば、ＨＤＤ（Ｈａｒｄ Ｄｉｓｋ Ｄｒｉｖｅ）等の磁気記憶部デバイス、半導体記憶デバイス、光記憶デバイス、または光磁気記憶デバイス等により構成される。このストレージ装置９１３は、ＣＰＵ９０１が実行するプログラムや各種データ、および外部から取得した各種のデータなどを格納する。
The
ドライブ９１５は、記録媒体用リーダライタであり、残留応力算出装置２００に内蔵、あるいは外付けされる。ドライブ９１５は、装着されている磁気ディスク、光ディスク、光磁気ディスク、または半導体メモリ等のリムーバブル記録媒体９２１に記録されている情報を読み出して、ＲＡＭ９０５に出力する。また、ドライブ９１５は、装着されている磁気ディスク、光ディスク、光磁気ディスク、または半導体メモリ等のリムーバブル記録媒体９２１に記録を書き込むことも可能である。リムーバブル記録媒体９２１は、例えば、ＣＤメディア、ＤＶＤメディア、Ｂｌｕ−ｒａｙメディア等である。また、リムーバブル記録媒体９２１は、コンパクトフラッシュ（登録商標）（ＣｏｍｐａｃｔＦｌａｓｈ：ＣＦ）、フラッシュメモリ、または、ＳＤメモリカード（Ｓｅｃｕｒｅ Ｄｉｇｉｔａｌ ｍｅｍｏｒｙ ｃａｒｄ）等であってもよい。また、リムーバブル記録媒体９２１は、例えば、非接触型ＩＣチップを搭載したＩＣカード（Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ ｃａｒｄ）または電子機器等であってもよい。
The
接続ポート９１７は、機器を残留応力算出装置２００に直接接続するためのポートである。接続ポート９１７の一例として、ＵＳＢ（Ｕｎｉｖｅｒｓａｌ Ｓｅｒｉａｌ Ｂｕｓ）ポート、ＩＥＥＥ１３９４ポート、ＳＣＳＩ（Ｓｍａｌｌ Ｃｏｍｐｕｔｅｒ Ｓｙｓｔｅｍ Ｉｎｔｅｒｆａｃｅ）ポート、ＲＳ−２３２Ｃポート等がある。この接続ポート９１７に外部接続機器９２３を接続することで、残留応力算出装置２００は、外部接続機器９２３から直接各種のデータを取得したり、外部接続機器９２３に各種のデータを提供したりする。
The
通信装置９１９は、例えば、通信網９２５に接続するための通信デバイス等で構成された通信インターフェースである。通信装置９１９は、例えば、有線または無線ＬＡＮ（Ｌｏｃａｌ Ａｒｅａ Ｎｅｔｗｏｒｋ）、Ｂｌｕｅｔｏｏｔｈ（登録商標）、またはＷＵＳＢ（Ｗｉｒｅｌｅｓｓ ＵＳＢ）用の通信カード等である。また、通信装置９１９は、光通信用のルータ、ＡＤＳＬ（Ａｓｙｍｍｅｔｒｉｃ Ｄｉｇｉｔａｌ Ｓｕｂｓｃｒｉｂｅｒ Ｌｉｎｅ）用のルータ、または、各種通信用のモデム等であってもよい。この通信装置９１９は、例えば、インターネットや他の通信機器との間で、例えばＴＣＰ／ＩＰ等の所定のプロトコルに則して信号等を送受信することができる。また、通信装置９１９に接続される通信網９２５は、有線または無線によって接続されたネットワーク等により構成され、例えば、インターネット、家庭内ＬＡＮ、赤外線通信、ラジオ波通信または衛星通信等であってもよい。
The
以上、本発明の実施形態に係る残留応力算出装置２００の機能を実現可能なハードウェア構成の一例を示した。上記の各構成要素は、汎用的な部材を用いて構成されていてもよいし、各構成要素の機能に特化したハードウェアにより構成されていてもよい。従って、本実施形態を実施する時々の技術レベルに応じて、適宜、利用するハードウェア構成を変更することが可能である。
Heretofore, an example of the hardware configuration capable of realizing the function of the residual
以上、添付図面を参照しながら本発明の好適な実施形態について詳細に説明したが、本発明はかかる例に限定されない。本発明の属する技術の分野における通常の知識を有する者であれば、特許請求の範囲に記載された技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、これらについても、当然に本発明の技術的範囲に属するものと了解される。 The preferred embodiments of the present invention have been described in detail above with reference to the accompanying drawings, but the present invention is not limited to such examples. It is obvious that a person having ordinary knowledge in the technical field to which the present invention pertains can come up with various changes or modifications within the scope of the technical idea described in the claims. Of course, it is understood that these also belong to the technical scope of the present invention.
例えば、上述の実施形態では、音弾性パラメータとして音速比Ｒ及び音響複屈折Ｂを例に挙げたが、本発明の実施形態に係る残留応力算出方法では、音弾性パラメータは上述の例に限定されるわけではない。例えば、音響複屈折Ｂの替わりに、表面ＳＨ波の音速異方性に関する音弾性パラメータを利用しても良い。 For example, in the above-described embodiment, the acoustic velocity ratio R and the acoustic birefringence B are exemplified as the acoustoelastic parameters. However, in the residual stress calculation method according to the embodiment of the present invention, the acoustoelastic parameters are limited to the above-described example. I don't mean. For example, instead of the acoustic birefringence B, a acoustoelastic parameter relating to the sound velocity anisotropy of the surface SH wave may be used.
１０ 残留応力測定装置
１００ 音速測定装置
１０１ パルサーレシーバー
１０３ プローブ
１１１ 正弦波発生装置
１１３ バースト波発生装置
１１５ 電磁超音波変換素子（ＥＭＡＴ）
１１７ 積分器
２００ 残量応力算出装置
２０１ 音速情報取得部
２０３ 音弾性パラメータ算出部
２０５ 塑性変形特徴量算出部
２０７ 音弾性パラメータ係数補正部
２０９ 記憶部
２１１ 残留応力算出部
２１３ 表示制御部
２１５ 表示部
Ｓ 金属板
DESCRIPTION OF
117
Claims (10)
前記金属板の塑性変形量又は塑性変形量に換算可能な値である塑性変形特徴量を算出する塑性変形特徴量算出部と、
前記音速情報取得部が取得した前記音速情報に基づいて、前記金属板中の複数の超音波音速から算出される、主応力和又は主応力差と相関のあるパラメータである音弾性パラメータを算出する音弾性パラメータ算出部と、
前記塑性変形特徴量算出部により算出された前記塑性変形特徴量を利用して、前記音弾性パラメータの応力による変化を表す係数の補正を行う音弾性パラメータ係数補正部と、
前記音弾性パラメータ算出部が算出した音弾性パラメータと、前記音弾性パラメータ係数補正部が補正した係数とに基づいて、前記金属板に残留している応力を算出する残留応力算出部と、
を備えることを特徴とする、残留応力算出装置。 A sound speed information acquisition unit for acquiring sound speed information related to a plurality of ultrasonic sound speeds in the metal plate measured by generating a plurality of ultrasonic waves in the plastically deformed metal plate;
A plastic deformation feature amount calculating unit that calculates a plastic deformation feature amount that is a value that can be converted into a plastic deformation amount or a plastic deformation amount of the metal plate;
Based on the sound speed information acquired by the sound speed information acquisition unit, a sound elasticity parameter, which is a parameter correlated with a main stress sum or a main stress difference, is calculated from a plurality of ultrasonic sound speeds in the metal plate. An acoustoelastic parameter calculator;
An acoustoelastic parameter coefficient correction unit that corrects a coefficient representing a change of the acoustoelastic parameter due to stress using the plastic deformation feature value calculated by the plastic deformation feature value calculating unit;
Based on the acoustoelastic parameter calculated by the acoustoelastic parameter calculating unit and the coefficient corrected by the acoustoelastic parameter coefficient correcting unit, a residual stress calculating unit that calculates the residual stress on the metal plate;
A residual stress calculation device comprising:
前記音弾性パラメータ係数補正部は、前記記憶部に格納された前記相関情報と、前記塑性変形特徴量算出部が算出した前記塑性変形特徴量とを利用して、前記係数を決定することを特徴とする、請求項１に記載の残留応力算出装置。 The residual stress calculation device further includes a storage unit storing correlation information representing a relationship between the coefficient and the plastic deformation feature amount,
The acoustoelastic parameter coefficient correction unit determines the coefficient using the correlation information stored in the storage unit and the plastic deformation feature amount calculated by the plastic deformation feature amount calculation unit. The residual stress calculation apparatus according to claim 1.
前記音弾性パラメータ係数補正部は、前記塑性変形特徴量算出部が算出した前記塑性変形特徴量に基づいて、以下の式１における係数ＣＲ及びＲ０、並びに、以下の式２における係数ＣＡ及びＢ０を決定することを特徴とする、請求項２に記載の残留応力算出装置。
ここで、下記式１及び式２において、σ１及びσ２は、金属板中に残留している主応力を表す。
The acoustoelastic parameter coefficient correction unit, based on the plastic deformation feature amount calculated by the plastic deformation feature amount calculation unit, coefficients C R and R 0 in the following formula 1 and coefficient C A in the following formula 2 And B 0 are determined, The residual stress calculation apparatus according to claim 2, wherein:
Here, in the following formulas 1 and 2, σ 1 and σ 2 represent main stresses remaining in the metal plate.
前記音速測定装置が測定した音速の測定結果を利用して、前記塑性変形した金属板中に残留する応力を算出する残留応力算出装置と、
を備え、
前記残留応力算出装置は、
前記音速測定装置から、前記塑性変形した金属板中を伝播する複数の超音波音速に関する測定結果である音速情報を取得する音速情報取得部と、
前記金属板の塑性変形量又は塑性変形量に換算可能な値である塑性変形特徴量を算出する塑性変形特徴量算出部と、
前記音速情報取得部が取得した前記音速情報に基づいて、前記金属板中の複数の超音波音速から算出される、主応力和又は主応力差と相関のあるパラメータである音弾性パラメータを算出する音弾性パラメータ算出部と、
前記塑性変形特徴量算出部により算出された前記塑性変形特徴量を利用して、前記音弾性パラメータの応力による変化を表す係数の補正を行う音弾性パラメータ係数補正部と、
前記音弾性パラメータ算出部が算出した音弾性パラメータと、前記音弾性パラメータ係数補正部が補正した係数とに基づいて、前記金属板に残留している応力を算出する残留応力算出部と、
を有することを特徴とする、残留応力測定装置。 A sound speed measuring device that measures the speed of sound of the plurality of ultrasonic waves propagating in the metal plate by generating a plurality of ultrasonic waves in the plastically deformed metal plate;
A residual stress calculation device that calculates a residual stress in the plastically deformed metal plate using a measurement result of the sound velocity measured by the sound velocity measurement device;
With
The residual stress calculation device is:
A sound speed information acquisition unit that acquires sound speed information, which is a measurement result of a plurality of ultrasonic sound speeds propagating in the plastically deformed metal plate, from the sound speed measurement device;
A plastic deformation feature amount calculating unit that calculates a plastic deformation feature amount that is a value that can be converted into a plastic deformation amount or a plastic deformation amount of the metal plate;
Based on the sound speed information acquired by the sound speed information acquisition unit, a sound elasticity parameter, which is a parameter correlated with a main stress sum or a main stress difference, is calculated from a plurality of ultrasonic sound speeds in the metal plate. An acoustoelastic parameter calculator;
An acoustoelastic parameter coefficient correction unit that corrects a coefficient representing a change of the acoustoelastic parameter due to stress using the plastic deformation feature value calculated by the plastic deformation feature value calculating unit;
Based on the acoustoelastic parameter calculated by the acoustoelastic parameter calculating unit and the coefficient corrected by the acoustoelastic parameter coefficient correcting unit, a residual stress calculating unit that calculates the residual stress on the metal plate;
A residual stress measuring device comprising:
前記金属板の塑性変形量又は塑性変形量に換算可能な値である塑性変形特徴量を算出するとともに、取得した前記音速情報に基づいて、前記金属板中の複数の超音波音速から算出される、主応力和又は主応力差と相関のあるパラメータである音弾性パラメータを算出するステップと、
算出した前記塑性変形特徴量を利用して、前記音弾性パラメータの応力による変化を表す係数の補正を行うステップと、
算出された音弾性パラメータと、補正された係数とに基づいて、前記金属板に残留している応力を算出するステップと、
を含むことを特徴とする、残留応力算出方法。 Obtaining sound velocity information regarding a plurality of ultrasonic sound speeds in the metal plate measured by generating a plurality of ultrasonic waves in the plastically deformed metal plate;
The plastic deformation amount of the metal plate or a plastic deformation feature value that is a value that can be converted into the plastic deformation amount is calculated, and is calculated from a plurality of ultrasonic sound speeds in the metal plate based on the acquired sound speed information. Calculating an acoustoelastic parameter which is a parameter correlated with the principal stress sum or principal stress difference;
Using the calculated plastic deformation feature amount, correcting a coefficient representing a change due to stress of the acoustoelastic parameter;
Calculating a residual stress on the metal plate based on the calculated acoustoelastic parameter and the corrected coefficient;
The residual stress calculation method characterized by including.
前記金属板中を伝播する複数の超音波音速に関する測定結果である音速情報を取得するステップと、
前記金属板の塑性変形量又は塑性変形量に換算可能な値である塑性変形特徴量を算出するとともに、取得した前記音速情報に基づいて、前記金属板中の複数の超音波音速から算出される、主応力和又は主応力差と相関のあるパラメータである音弾性パラメータを算出するステップと、
算出した前記塑性変形特徴量を利用して、前記音弾性パラメータの応力による変化を表す係数の補正を行うステップと、
算出された音弾性パラメータと、補正された係数とに基づいて、前記金属板に残留している応力を算出するステップと、
を含むことを特徴とする、残留応力測定方法。 Measuring the speed of sound of the plurality of ultrasonic waves propagating through the metal plate by generating a plurality of ultrasonic waves in the plastically deformed metal plate;
Obtaining sound speed information that is a measurement result of a plurality of ultrasonic sound speeds propagating in the metal plate;
The plastic deformation amount of the metal plate or a plastic deformation feature value that is a value that can be converted into the plastic deformation amount is calculated, and is calculated from a plurality of ultrasonic sound speeds in the metal plate based on the acquired sound speed information. Calculating an acoustoelastic parameter which is a parameter correlated with the principal stress sum or principal stress difference;
Using the calculated plastic deformation feature amount, correcting a coefficient representing a change due to stress of the acoustoelastic parameter;
Calculating a residual stress on the metal plate based on the calculated acoustoelastic parameter and the corrected coefficient;
A residual stress measurement method comprising:
塑性変形した金属板中に複数の超音波を発生させることで測定された金属板中の複数の超音波音速に関する音速情報を取得する音速情報取得機能と、
前記金属板の塑性変形量又は塑性変形量に換算可能な値である塑性変形特徴量を算出する塑性変形特徴量算出機能と、
前記音速情報取得機能が取得した前記音速情報に基づいて、前記金属板中の複数の超音波音速から算出される、主応力和又は主応力差と相関のあるパラメータである音弾性パラメータを算出する音弾性パラメータ算出機能と、
前記塑性変形特徴量算出機能が算出した前記塑性変形特徴量を利用して、前記音弾性パラメータの応力による変化を表す係数の補正を行う音弾性パラメータ係数補正機能と、
前記音弾性パラメータ算出機能が算出した音弾性パラメータと、前記音弾性パラメータ係数補正機能が補正した係数とに基づいて、前記金属板に残留している応力を算出する残留応力算出機能と、
を実現させるためのプログラム。
On the computer,
A sound speed information acquisition function for acquiring sound speed information related to a plurality of ultrasonic sound speeds in a metal plate measured by generating a plurality of ultrasonic waves in a plastically deformed metal plate,
A plastic deformation feature amount calculation function for calculating a plastic deformation feature amount that is a value that can be converted into a plastic deformation amount or a plastic deformation amount of the metal plate;
Based on the sound speed information acquired by the sound speed information acquisition function, a sound elasticity parameter, which is a parameter correlated with a main stress sum or a main stress difference, is calculated from a plurality of ultrasonic sound speeds in the metal plate. A sound elasticity parameter calculation function;
Acoustoelastic parameter coefficient correction function for correcting a coefficient representing a change due to stress of the acoustoelastic parameter using the plastic deformation feature value calculated by the plastic deformation feature value calculating function;
A residual stress calculation function for calculating a stress remaining in the metal plate based on the acoustoelastic parameter calculated by the acoustoelastic parameter calculation function and the coefficient corrected by the acoustoelastic parameter coefficient correction function;
A program to realize
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2010066987A JP5251911B2 (en) | 2010-03-23 | 2010-03-23 | Residual stress calculation device, residual stress measurement device, residual stress calculation method, residual stress measurement method, and program |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2010066987A JP5251911B2 (en) | 2010-03-23 | 2010-03-23 | Residual stress calculation device, residual stress measurement device, residual stress calculation method, residual stress measurement method, and program |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2011196953A JP2011196953A (en) | 2011-10-06 |
| JP5251911B2 true JP5251911B2 (en) | 2013-07-31 |
Family
ID=44875352
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2010066987A Expired - Fee Related JP5251911B2 (en) | 2010-03-23 | 2010-03-23 | Residual stress calculation device, residual stress measurement device, residual stress calculation method, residual stress measurement method, and program |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP5251911B2 (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102818763B (en) * | 2012-07-30 | 2014-11-19 | 首钢总公司 | Hot-rolled steel plate residual stress calculating method suitable to production field |
| JP6362533B2 (en) * | 2013-12-24 | 2018-07-25 | 株式会社神戸製鋼所 | Residual stress evaluation method and residual stress evaluation apparatus |
| RU2556331C1 (en) * | 2014-02-28 | 2015-07-10 | Федеральное государственное бюджетное учреждение науки Физико-технический институт Уральского отделения Российской академии наук (ФТИ УрО РАН) | Apparatus for determining acoustoelastic communication coefficient |
| RU2589219C2 (en) * | 2014-08-13 | 2016-07-10 | Общество с ограниченной ответственностью "ИНКОТЕС" (ООО "ИНКОТЕС") | Method of evaluating stress-strain state of elements of complex structures |
| GB2529484A (en) * | 2014-08-22 | 2016-02-24 | Univ Sheffield | Deriving contact stress or contact loadusing ultrasound data |
| JP6529887B2 (en) * | 2015-10-22 | 2019-06-12 | 株式会社神戸製鋼所 | Residual stress evaluation method |
| CN106770402B (en) * | 2017-01-11 | 2023-08-04 | 中国工程物理研究院核物理与化学研究所 | Three-dimensional calibration measuring device for neutron diffraction stress analysis |
| JP6981521B1 (en) | 2020-12-08 | 2021-12-15 | Ｊｆｅスチール株式会社 | Residual stress calculation method |
| CN113074849A (en) * | 2021-03-26 | 2021-07-06 | 重庆交通大学 | Concrete surface absolute stress measuring method based on laser ultrasonic technology |
| CN115144112A (en) * | 2022-06-23 | 2022-10-04 | 武汉中科创新技术股份有限公司 | Residual stress detection method and device, electronic equipment and storage medium |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH07286916A (en) * | 1994-04-15 | 1995-10-31 | Sumitomo Metal Ind Ltd | Method for measuring residual stress |
| JPH11264775A (en) * | 1998-03-18 | 1999-09-28 | Osaka Gas Co Ltd | Pipe measuring and analyzing method and pipe analyzing device |
| JP4031711B2 (en) * | 2002-06-17 | 2008-01-09 | 新日本製鐵株式会社 | Residual stress distribution measuring apparatus and residual stress distribution measuring method |
| JP4087730B2 (en) * | 2003-03-26 | 2008-05-21 | 株式会社酒井鉄工所 | Stress anisotropy detection method and stress diagnosis method |
| JP2008076387A (en) * | 2006-08-24 | 2008-04-03 | Toshiba Corp | Ultrasonic stress measuring method and instrument |
| JP2009025093A (en) * | 2007-07-18 | 2009-02-05 | Nichizou Tec:Kk | Electromagnetic ultrasonic measuring device, and measuring method of plate thickness and stress using electromagnetic ultrasonic wave |
-
2010
- 2010-03-23 JP JP2010066987A patent/JP5251911B2/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2011196953A (en) | 2011-10-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5251911B2 (en) | Residual stress calculation device, residual stress measurement device, residual stress calculation method, residual stress measurement method, and program | |
| JP4954729B2 (en) | Concrete pile soundness evaluation support device, soundness evaluation support method, and soundness evaluation support program | |
| WO2017071605A1 (en) | Elasticity detection method and device | |
| EP3592227A1 (en) | Methods and systems to measure and evaluate stability of medical implants | |
| WO2015159378A1 (en) | Ultrasonic inspection device and ultrasonic inspection method | |
| EP2426490B1 (en) | Ultrasonic testing method | |
| JP2019070627A (en) | Nondestructive inspection system | |
| JP4764921B2 (en) | Ultrasonic exploration method using resonance phenomenon | |
| JP2011013150A (en) | Method and device for measuring crack state of concrete structure | |
| CN103075981A (en) | Ultrasonic thickness measuring method | |
| Santos et al. | Metrological validation of a measurement procedure for the characterization of a biological ultrasound tissue-mimicking material | |
| JP2006200901A (en) | Ultrasonic inspection method and device | |
| JP2011242332A (en) | Ultrasonic nondestructive measurement method, ultrasonic nondestructive measurement instrument, and program | |
| JP2010038696A (en) | Non-destructive evaluation method of degree of metal fatigue damage and ultrasonic metal fatigue damage degree measuring instrument | |
| JP2007017300A (en) | Surface inspection device and surface inspection method | |
| JP2009025093A (en) | Electromagnetic ultrasonic measuring device, and measuring method of plate thickness and stress using electromagnetic ultrasonic wave | |
| JP2015078910A (en) | Ultrasonic measuring device and calibration method thereof | |
| JP6529853B2 (en) | Residual stress evaluation method | |
| JP5447323B2 (en) | Sonic velocity distribution measuring apparatus, sonic velocity distribution measuring method and program | |
| JP2001337077A (en) | Device for non-destructively inspecting exfoliation in concrete structure | |
| JP4745083B2 (en) | Stress measuring apparatus and measuring method thereof | |
| Zhang et al. | A new nondestructive technique for measuring pressure in vessels by surface waves | |
| JP2005326238A (en) | Crack development monitoring method and its device by ultrasonic tofd method | |
| JP5240218B2 (en) | Crystal grain size measuring apparatus, method and program | |
| JP5630790B2 (en) | Method for estimating the rate of thinning |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20120209 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20130319 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20130321 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20130401 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 5251911 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20160426 Year of fee payment: 3 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| LAPS | Cancellation because of no payment of annual fees |