JP6640451B2 - Residual stress evaluation method - Google Patents

Residual stress evaluation method Download PDF

Info

Publication number
JP6640451B2
JP6640451B2 JP2015021571A JP2015021571A JP6640451B2 JP 6640451 B2 JP6640451 B2 JP 6640451B2 JP 2015021571 A JP2015021571 A JP 2015021571A JP 2015021571 A JP2015021571 A JP 2015021571A JP 6640451 B2 JP6640451 B2 JP 6640451B2
Authority
JP
Japan
Prior art keywords
construction
impact pressure
water jet
residual stress
value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2015021571A
Other languages
Japanese (ja)
Other versions
JP2016140974A (en
Inventor
直輝 小川
直輝 小川
真弓 越智
真弓 越智
展之 堀
展之 堀
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Mitsubishi Heavy Industries Ltd
Original Assignee
Mitsubishi Heavy Industries Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Mitsubishi Heavy Industries Ltd filed Critical Mitsubishi Heavy Industries Ltd
Priority to JP2015021571A priority Critical patent/JP6640451B2/en
Priority to PCT/JP2015/062096 priority patent/WO2016125319A1/en
Priority to EP15881136.4A priority patent/EP3241644B1/en
Priority to US15/546,077 priority patent/US20180010205A1/en
Publication of JP2016140974A publication Critical patent/JP2016140974A/en
Application granted granted Critical
Publication of JP6640451B2 publication Critical patent/JP6640451B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • CCHEMISTRY; METALLURGY
    • C21METALLURGY OF IRON
    • C21DMODIFYING THE PHYSICAL STRUCTURE OF FERROUS METALS; GENERAL DEVICES FOR HEAT TREATMENT OF FERROUS OR NON-FERROUS METALS OR ALLOYS; MAKING METAL MALLEABLE, e.g. BY DECARBURISATION OR TEMPERING
    • C21D7/00Modifying the physical properties of iron or steel by deformation
    • C21D7/02Modifying the physical properties of iron or steel by deformation by cold working
    • C21D7/04Modifying the physical properties of iron or steel by deformation by cold working of the surface
    • C21D7/06Modifying the physical properties of iron or steel by deformation by cold working of the surface by shot-peening or the like
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B23MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
    • B23PMETAL-WORKING NOT OTHERWISE PROVIDED FOR; COMBINED OPERATIONS; UNIVERSAL MACHINE TOOLS
    • B23P17/00Metal-working operations, not covered by a single other subclass or another group in this subclass
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/32Investigating strength properties of solid materials by application of mechanical stress by applying repeated or pulsating forces
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/40Investigating hardness or rebound hardness
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B24GRINDING; POLISHING
    • B24CABRASIVE OR RELATED BLASTING WITH PARTICULATE MATERIAL
    • B24C1/00Methods for use of abrasive blasting for producing particular effects; Use of auxiliary equipment in connection with such methods
    • B24C1/10Methods for use of abrasive blasting for producing particular effects; Use of auxiliary equipment in connection with such methods for compacting surfaces, e.g. shot-peening
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L5/00Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
    • G01L5/0047Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes measuring forces due to residual stresses
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02Details not specific for a particular testing method
    • G01N2203/0202Control of the test
    • G01N2203/0212Theories, calculations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02Details not specific for a particular testing method
    • G01N2203/0202Control of the test
    • G01N2203/0212Theories, calculations
    • G01N2203/0216Finite elements
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Immunology (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Mechanical Engineering (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Materials Engineering (AREA)
  • Metallurgy (AREA)
  • Organic Chemistry (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Description

本開示は、ウォータジェットピーニングの施工により改善される残留応力の評価方法、この方法を実行するためのプログラム、この方法を実行する装置に関する。   The present disclosure relates to a method for evaluating residual stress improved by performing water jet peening, a program for executing the method, and an apparatus for executing the method.

構造物(金属など)の高温水中における劣化現象の一つとして応力腐食割れ(SCC:Stress Corrosion Cracking)が知られている。このSCCの原因が溶接部などに発生した引張残留応力である場合、その劣化防止対策として、この構造物に対してウォータジェットピーニング(以下、WJPとする)が施工されることがある。このWJPは、構造物(対象物)の表面近傍の残留応力を改善可能な技術であり、水中で施工可能な技術である。より詳細には、WJPでは、対象物の表面に対してノズルから高圧水を噴射し、この噴射により発生する気泡(キャビテーション気泡)の崩壊時の衝撃圧を利用して対象物の表面に塑性変形を起こさせる。そして、こうした塑性変形を通して、対象物の表面近傍の残留応力が低減され、あるいは、その表面近傍の引張残留応力が圧縮残留応力に改善されることでSCCが抑制される。   Stress corrosion cracking (SCC) is known as one of the degradation phenomena of a structure (metal or the like) in high-temperature water. When the cause of the SCC is a residual tensile stress generated in a weld or the like, water jet peening (hereinafter, referred to as WJP) may be applied to the structure as a measure to prevent the deterioration. This WJP is a technology capable of improving residual stress in the vicinity of the surface of a structure (object) and is a technology that can be constructed in water. More specifically, in WJP, high-pressure water is jetted from a nozzle to the surface of an object, and plastic deformation is applied to the surface of the object by utilizing the impact pressure at the time of collapse of bubbles (cavitation bubbles) generated by this injection. Wake up. Then, through such plastic deformation, the residual stress near the surface of the object is reduced, or the tensile residual stress near the surface is improved to the compressive residual stress, thereby suppressing SCC.

一方、WJPを施工後の構造物の残留応力を評価する方法も幾つか提案されており、下記の特許文献1では、施工部位の残留応力を予測する方法が開示されている。具体的には、特許文献1では、WJP施工後の残留応力とキャビテーション気泡の崩壊圧との相関関係と、WJP施工対象表面に作用するキャビテーションエネルギーに基づいてキャビテーション気泡の崩壊圧を求めることができるという知見のもと、ノズルから噴出する噴流を解析して得られるキャビテーション気泡の気泡内圧と気泡数密度を用いてキャビテーションエネルギーを算出し、キャビテーションエネルギーに基づいてキャビテーション気泡の崩壊圧を算出し、この崩壊圧を用いて施工対象表面のWJP施工後の残留応力を算出している。   On the other hand, several methods for evaluating the residual stress of a structure after WJP is applied have been proposed, and Patent Literature 1 below discloses a method for estimating the residual stress of a construction site. Specifically, in Patent Document 1, the collapse pressure of cavitation bubbles can be obtained based on the correlation between the residual stress after WJP construction and the collapse pressure of cavitation bubbles and the cavitation energy acting on the surface to be subjected to WJP construction. Based on the knowledge, the cavitation energy is calculated using the internal pressure and the number density of the cavitation bubbles obtained by analyzing the jet flow ejected from the nozzle, and the collapse pressure of the cavitation bubbles is calculated based on the cavitation energy. The residual stress after the WJP on the construction target surface is calculated using the collapse pressure.

特許第5011416号公報Japanese Patent No. 50111416

上述の通り、特許文献1では、キャビテーションエネルギー(キャビテーション気泡の気泡内圧と気泡数密度)の解析結果に基づいてキャビテーション気泡の崩壊圧を得ている。これに対して、本開示では、WJPのキャビテーション気泡の発生・消滅状況に着目し、キャビテーション気泡により生じる現実の衝撃圧に相関する相関値を解析により求め、この解析された相関値と実験値を対応づけることでキャビテーション気泡により生じる現実の衝撃圧の分布(施工対象表面の各位置における衝撃圧)を予測できることを新たに見出した。そして、この衝撃圧の予測値からWJP施工後の残留応力が評価される。   As described above, in Patent Document 1, the collapse pressure of the cavitation bubble is obtained based on the analysis result of the cavitation energy (the internal pressure of the cavitation bubble and the bubble number density). In contrast, in the present disclosure, attention is paid to the occurrence and disappearance of cavitation bubbles in WJP, a correlation value that correlates to the actual impact pressure generated by cavitation bubbles is obtained by analysis, and the analyzed correlation value and the experimental value are calculated. It has been newly found that it is possible to predict the distribution of the actual impact pressure caused by the cavitation bubbles (impact pressure at each position on the construction target surface) by associating them. Then, the residual stress after WJP construction is evaluated from the predicted value of the impact pressure.

上述の事情に鑑みて、本発明の少なくとも一実施形態は、ウォータジェットピーニング(WJP)におけるキャビテーション気泡の発生・消滅状況の解析に基づいた残留応力の評価方法を提供することを目的とする。   In view of the above circumstances, at least one embodiment of the present invention aims to provide a method for evaluating residual stress based on an analysis of the generation and disappearance of cavitation bubbles in water jet peening (WJP).

(1)本発明の少なくとも一実施形態に係る残留応力評価方法は、
施工対象へのウォータジェットピーニングの施工条件を設定する条件設定工程と、
ノズルモデルから施工対象モデルへ流体を噴出させたときの噴流を前記施工条件に応じて解析して、前記施工対象モデルの表面の各位置での、前記流体の単位体積に含まれる気泡の体積率であるボイド率、及び前記流体の単位体積中で単位時間に崩壊する前記気泡の体積率である崩壊率を求める解析工程と、
前記各位置での前記ボイド率と前記崩壊率との積である衝撃圧相関値を求める衝撃圧相関値算出工程と、
前記施工条件による前記ウォータジェットピーニングによって前記施工対象の表面に作用する衝撃圧の実験値である衝撃圧実験値を取得する実験値取得工程と、
前記施工対象モデルの表面の各位置における前記衝撃圧相関値と、前記施工条件による前記ウォータジェットピーニングによって前記施工対象の表面に作用する前記衝撃圧実験値とを対応づけて、前記施工条件による前記ウォータジェットピーニングによって前記施工対象の表面の各位置に作用する前記衝撃圧の予測値である衝撃圧予測値を取得する予測工程と、
前記予測工程により取得された前記衝撃圧予測値を入力条件として、前記施工条件による前記ウォータジェットピーニングを施工後の前記施工対象における残留応力を算出する残留応力評価工程と、を備える。
(1) The residual stress evaluation method according to at least one embodiment of the present invention includes:
A condition setting step of setting water jet peening construction conditions for the construction object,
Analyzing the jet flow when the fluid is jetted from the nozzle model to the construction target model according to the construction conditions, at each position on the surface of the construction target model, the volume fraction of bubbles contained in the unit volume of the fluid An analysis step to determine the void fraction, and the collapse rate, which is the volume fraction of the bubbles that collapse in a unit time in the unit volume of the fluid,
Impact pressure correlation value calculation step of obtaining an impact pressure correlation value that is a product of the void rate and the collapse rate at each position,
An experimental value acquisition step of acquiring an impact pressure experimental value that is an experimental value of an impact pressure acting on the surface of the construction target by the water jet peening according to the construction condition,
The impact pressure correlation value at each position of the surface of the construction target model, and the impact pressure experimental value acting on the surface of the construction target by the water jet peening according to the construction condition, the impact pressure experimental value according to the construction condition Prediction step of obtaining an impact pressure predicted value that is a predicted value of the impact pressure acting on each position of the construction target surface by water jet peening,
A residual stress evaluation step of calculating a residual stress in the construction target after the water jet peening under the construction conditions, using the impact pressure predicted value acquired in the prediction step as an input condition.

上記(1)の構成によれば、ウォータジェットピーニング(WJP)におけるキャビテーション気泡の発生・消滅状況を解析することで、施工対象表面近傍に生じるWJPによる現実の衝撃圧の分布(施工対象表面の各位置での衝撃圧)を予測し、予測された衝撃圧に基づいてWJP施工後の残留応力を評価することができる。   According to the configuration of the above (1), by analyzing the occurrence and disappearance of cavitation bubbles in water jet peening (WJP), the distribution of the actual impact pressure due to WJP generated near the construction target surface (each of the construction target surface) (Impact pressure at the position), and the residual stress after WJP construction can be evaluated based on the predicted impact pressure.

(2)幾つかの実施形態では、上記(1)の構成において、
前記衝撃圧予測値は、前記施工対象モデルの表面の各位置における前記衝撃圧相関値と前記衝撃圧実験値とを対応づける係数kを決定することにより取得される。
上記(2)の構成によれば、上記の係数kを決定することで、施工対象モデルの表面の各位置における衝撃圧相関値と衝撃圧実験値とを対応づけることができ、これによって、施工対象の表面近傍におけるウォータジェットピーニング(WJP)施工後の現実の残留応力を解析・評価することができる。
(2) In some embodiments, in the configuration of the above (1),
The shock pressure predicted value is obtained by determining a coefficient k that associates the shock pressure correlation value at each position on the surface of the construction target model with the shock pressure experimental value.
According to the configuration of the above (2), by determining the coefficient k, it is possible to associate the shock pressure correlation value at each position on the surface of the construction target model with the shock pressure experiment value. It is possible to analyze and evaluate the actual residual stress after the water jet peening (WJP) in the vicinity of the target surface.

(3)幾つかの実施形態では、上記(1)〜(2)の構成において、
前記残留応力の目標値を設定する目標値設定工程と、
前記残留応力が前記目標値を満たさない場合に前記施工条件を変更する施工条件変更工程と、
前記施工条件変更工程により変更された変更後の施工条件により前記解析工程と、前記衝撃圧相関値算出工程と、前記残留応力評価工程とを行う再評価工程と、
前記再評価工程によって算出される残留応力が前記目標値を満たす場合には、前記再評価工程で用いられた前記施工条件を前記施工対象への前記施工条件として決定する施工条件決定工程と、をさらに備える。
上記(3)の構成によれば、残留応力の目標値と、任意の施工条件によるウォータジェットピーニング(WJP)施工後の残留応力の解析値との比較に基づいて、この目標値を満たす施工条件が決定される。このため、プラントの仕様に合わせて実績のない施工条件によりWJPを施工するような場合でも、WJP施工後の残留応力を解析により評価することができ、プラントの仕様に応じた適切な施工条件を速やかに決定することができる。さらに、このようにして決定された施工条件によりWJPを現実に施工することで、信頼性のあるWJPの施工を行うことができる。また、所望の施工条件による施工実施が可能なWJP施工装置の設計に利用することもできる。
(3) In some embodiments, in the above configurations (1) and (2),
A target value setting step of setting a target value of the residual stress,
An execution condition changing step of changing the execution condition when the residual stress does not satisfy the target value,
The analysis step, the impact pressure correlation value calculation step, and the re-evaluation step of performing the residual stress evaluation step by the changed execution condition changed by the execution condition changing step,
When the residual stress calculated by the re-evaluation step satisfies the target value, a construction condition determination step of determining the construction condition used in the re-evaluation step as the construction condition for the construction target, Further prepare.
According to the configuration of the above (3), on the basis of a comparison between the target value of the residual stress and the analysis value of the residual stress after the water jet peening (WJP) is performed under an arbitrary processing condition, the processing condition satisfying the target value Is determined. For this reason, even when WJP is performed under construction conditions that have no track record in accordance with plant specifications, the residual stress after WJP construction can be evaluated by analysis, and appropriate construction conditions according to plant specifications can be determined. Can be determined promptly. Furthermore, reliable WJP construction can be performed by actually constructing WJP according to the construction conditions determined in this way. Further, the present invention can be used for designing a WJP construction apparatus capable of performing construction under desired construction conditions.

(4)幾つかの実施形態では、上記(3)の構成において、
前記再評価工程によって算出される残留応力が前記目標値を満たさない場合には、前記施工条件変更工程によりさらに前記施工条件を変更し、前記施工条件変更工程によりさらに変更された変更後の施工条件により前記再評価工程を行う。
上記(4)の構成によれば、任意の施工条件によるウォータジェットピーニング(WJP)施工後の残留応力の評価結果がその目標値を満たさない場合には、他の任意の施工条件による同様の評価が行われる。これによって、目標値が満たされる施工条件を決定することができる。
(4) In some embodiments, in the configuration of the above (3),
If the residual stress calculated by the reevaluation step does not satisfy the target value, further change the execution condition by the execution condition change step, the changed execution condition further changed by the execution condition change step Performs the reevaluation step.
According to the configuration of the above (4), when the evaluation result of the residual stress after the water jet peening (WJP) under the arbitrary construction conditions does not satisfy the target value, the similar evaluation under the other arbitrary construction conditions Is performed. As a result, it is possible to determine the construction condition that satisfies the target value.

(5)幾つかの実施形態では、上記(1)〜(4)の構成において、
前記施工条件は、前記ウォータジェットピーニングによるウォータジェットの噴射時間、前記ウォータジェットの噴射速度、前記ウォータジェットの流量、前記ウォータジェットピーニングの施工範囲、前記ウォータジェットの噴射距離、前記気泡の半径、ノズル角度、前記施工対象の表面の傾斜角度のうちの少なくとも1つを含む。
上記(5)の構成によれば、WJP施工後の残留応力に影響する可能性のある各種条件が施工条件に含まれるので、任意の施工条件よるWJP施工後の残留応力を精度よく評価することができる。
(5) In some embodiments, in the above configurations (1) to (4),
The execution conditions are: water jet injection time by the water jet peening, the water jet injection speed, the water jet flow rate, the water jet peening execution range, the water jet injection distance, the bubble radius, and the nozzle. And at least one of an angle and a tilt angle of the surface to be applied.
According to the above configuration (5), since various conditions that may affect the residual stress after WJP are included in the execution conditions, it is necessary to accurately evaluate the residual stress after WJP under any arbitrary conditions. Can be.

(6)本発明の少なくとも一実施形態に係るウォータジェットピーニングの評価プログラムは、
コンピュータに、
施工対象へのウォータジェットピーニングの施工条件を設定する条件設定工程と、
ノズルモデルから施工対象モデルへ流体を噴出させたときの噴流を前記施工条件に応じて解析して、前記施工対象モデルの表面の各位置での、前記流体の単位体積に含まれる気泡の体積率であるボイド率、及び前記流体の単位体積中で単位時間に崩壊する前記気泡の体積率である崩壊率を求める解析工程と、
前記各位置での前記ボイド率と前記崩壊率との積である衝撃圧相関値を求める衝撃圧相関値算出工程と、
前記施工条件による前記ウォータジェットピーニングによって前記施工対象の表面に作用する衝撃圧の実験値である衝撃圧実験値を取得する実験値取得工程と、
前記施工対象モデルの表面の各位置における前記衝撃圧相関値と、前記施工条件による前記ウォータジェットピーニングによって前記施工対象の表面に作用する前記衝撃圧実験値とを対応づけて、前記施工条件による前記ウォータジェットピーニングによって前記施工対象の表面の各位置に作用する前記衝撃圧の予測値である衝撃圧予測値を取得する予測工程と、
前記予測工程により取得された前記衝撃圧予測値を入力条件として、前記施工条件による前記ウォータジェットピーニングを施工後の前記施工対象における残留応力を算出する残留応力評価工程と、を実行させる。
(6) The water jet peening evaluation program according to at least one embodiment of the present invention includes:
On the computer,
A condition setting step of setting water jet peening construction conditions for the construction object,
Analyzing the jet flow when the fluid is jetted from the nozzle model to the construction target model according to the construction conditions, at each position on the surface of the construction target model, the volume fraction of bubbles contained in the unit volume of the fluid An analysis step to determine the void fraction, and the collapse rate, which is the volume fraction of the bubbles that collapse in a unit time in the unit volume of the fluid,
Impact pressure correlation value calculation step of obtaining an impact pressure correlation value that is a product of the void rate and the collapse rate at each position,
An experimental value acquisition step of acquiring an impact pressure experimental value that is an experimental value of an impact pressure acting on the surface of the construction target by the water jet peening according to the construction condition,
The impact pressure correlation value at each position of the surface of the construction target model, and the impact pressure experimental value acting on the surface of the construction target by the water jet peening according to the construction condition, the impact pressure experimental value according to the construction condition Prediction step of obtaining an impact pressure predicted value that is a predicted value of the impact pressure acting on each position of the construction target surface by water jet peening,
A residual stress evaluation step of calculating a residual stress in the object after the water jet peening according to the application condition, using the impact pressure predicted value acquired in the prediction step as an input condition.

上記(6)の構成によれば、ウォータジェットピーニング(WJP)におけるキャビテーション気泡の発生・消滅状況を解析することで、施工対象表面近傍に生じるWJPによる現実の衝撃圧の分布(施工対象表面の各位置での衝撃圧)を予測し、予測された衝撃圧に基づいてWJP施工後の残留応力を評価することができる。   According to the configuration of the above (6), by analyzing the generation and disappearance of cavitation bubbles in water jet peening (WJP), the distribution of the actual impact pressure due to WJP generated near the construction target surface (each of the construction target surface) (Impact pressure at the position), and the residual stress after WJP construction can be evaluated based on the predicted impact pressure.

(7)幾つかの実施形態では、上記(6)の構成において、
前記衝撃圧予測値は、前記施工対象モデルの表面の各位置における前記衝撃圧相関値と前記衝撃圧実験値とを対応づける係数kを決定することにより取得される。
上記(7)の構成によれば、上記の係数kを決定することで、施工対象モデルの表面の各位置における衝撃圧相関値と衝撃圧実験値とを対応づけることができ、これによって、施工対象の表面近傍におけるウォータジェットピーニング(WJP)施工後の現実の残留応力を解析・評価することができる。
(7) In some embodiments, in the configuration of the above (6),
The shock pressure predicted value is obtained by determining a coefficient k that associates the shock pressure correlation value at each position on the surface of the construction target model with the shock pressure experimental value.
According to the configuration of (7), by determining the coefficient k, it is possible to associate the shock pressure correlation value at each position on the surface of the construction target model with the shock pressure experiment value. It is possible to analyze and evaluate the actual residual stress after the water jet peening (WJP) in the vicinity of the target surface.

(8)幾つかの実施形態では、上記(6)〜(7)の構成において、
前記残留応力の目標値を設定する目標値設定工程と、
前記残留応力が前記目標値を満たさない場合に前記施工条件を変更する施工条件変更工程と、
前記施工条件変更工程により変更された変更後の施工条件により前記解析工程と、前記衝撃圧相関値算出工程と、前記残留応力評価工程とを行う再評価工程と、
前記再評価工程によって算出される残留応力が前記目標値を満たす場合には、前記再評価工程で用いられた前記施工条件を前記施工対象への前記施工条件として決定する施工条件決定工程と、をさらに実行させる。
上記(8)の構成によれば、残留応力の目標値と、任意の施工条件によるウォータジェットピーニング(WJP)施工後の残留応力の解析値との比較に基づいて、この目標値を満たす施工条件が決定される。このため、プラントの仕様に合わせて実績のない施工条件によりWJPを施工するような場合でも、WJP施工後の残留応力を解析により評価することができ、プラントの仕様に応じた適切な施工条件を速やかに決定することができる。さらに、このようにして決定された施工条件によりWJPを現実に施工することで、信頼性のあるWJPの施工を行うことができる。また、所望の施工条件による施工実施が可能なWJP施工装置の設計に利用することもできる。
(8) In some embodiments, in the configuration of the above (6) to (7),
A target value setting step of setting a target value of the residual stress,
An execution condition changing step of changing the execution condition when the residual stress does not satisfy the target value,
The analysis step, the impact pressure correlation value calculation step, and the re-evaluation step of performing the residual stress evaluation step by the changed execution condition changed by the execution condition changing step,
When the residual stress calculated by the re-evaluation step satisfies the target value, a construction condition determination step of determining the construction condition used in the re-evaluation step as the construction condition for the construction target, Let it run more.
According to the configuration of the above (8), on the basis of a comparison between the target value of the residual stress and the analysis value of the residual stress after the water jet peening (WJP) under the arbitrary execution conditions, the execution condition satisfying the target value is obtained. Is determined. For this reason, even when WJP is performed under construction conditions that have no track record in accordance with plant specifications, the residual stress after WJP construction can be evaluated by analysis, and appropriate construction conditions according to plant specifications can be determined. Can be determined promptly. Furthermore, reliable WJP construction can be performed by actually constructing WJP according to the construction conditions determined in this way. Further, the present invention can be used for designing a WJP construction apparatus capable of performing construction under desired construction conditions.

(9)幾つかの実施形態では、上記(8)の構成において、
前記再評価工程によって算出される残留応力が前記目標値を満たさない場合には、前記施工条件変更工程によりさらに前記施工条件を変更し、前記施工条件変更工程によりさらに変更された変更後の施工条件により前記再評価工程を行う。
上記(9)の構成によれば、任意の施工条件によるウォータジェットピーニング(WJP)施工後の残留応力の評価結果がその目標値を満たさない場合には、他の任意の施工条件による同様の評価が行われる。これによって、目標値が満たされる施工条件を決定することができる。
(9) In some embodiments, in the configuration of the above (8),
If the residual stress calculated by the reevaluation step does not satisfy the target value, further change the execution condition by the execution condition change step, the changed execution condition further changed by the execution condition change step Performs the reevaluation step.
According to the configuration of the above (9), when the evaluation result of the residual stress after the water jet peening (WJP) under the arbitrary construction conditions does not satisfy the target value, the similar evaluation under the other arbitrary construction conditions Is performed. As a result, it is possible to determine the construction condition that satisfies the target value.

(10)幾つかの実施形態では、上記(6)〜(9)の構成において、
前記施工条件は、前記ウォータジェットピーニングによるウォータジェットの噴射時間、前記ウォータジェットの噴射速度、前記ウォータジェットの流量、前記ウォータジェットピーニングの施工範囲、前記ウォータジェットの噴射距離、前記気泡の半径、ノズル角度、前記施工対象の表面の傾斜角度のうちの少なくとも1つを含む。
上記(10)の構成によれば、WJP施工後の残留応力に影響する可能性のある各種条件が施工条件に含まれるので、任意の施工条件よるWJP施工後の残留応力を精度よく評価することができる。
(10) In some embodiments, in the above configurations (6) to (9),
The execution conditions are: water jet injection time by the water jet peening, the water jet injection speed, the water jet flow rate, the water jet peening execution range, the water jet injection distance, the bubble radius, and the nozzle. And at least one of an angle and a tilt angle of the surface to be applied.
According to the configuration of the above (10), since various conditions that may affect the residual stress after WJP are included in the execution conditions, it is necessary to accurately evaluate the residual stress after WJP under any execution conditions. Can be.

(11)本発明の少なくとも一実施形態に係るウォータジェットピーニングの評価装置は、
施工対象へのウォータジェットピーニングの施工条件を受け付ける条件受付部と、
ノズルモデルから施工対象モデルへ流体を噴出させたときの噴流を前記施工条件に応じて解析して、前記施工対象モデルの表面の各位置での、前記流体の単位体積に含まれる気泡の体積率であるボイド率、及び前記流体の単位体積中で単位時間に崩壊する前記気泡の体積率である崩壊率を求める解析部と、
前記各位置での前記ボイド率と前記崩壊率との積である衝撃圧相関値を求める衝撃圧相関値算出部と、
前記施工条件による前記ウォータジェットピーニングによって前記施工対象の表面に作用する衝撃圧の実験値である衝撃圧実験値を取得する実験値取得部と、
前記施工対象モデルの表面の各位置における前記衝撃圧相関値と、前記施工条件による前記ウォータジェットピーニングによって前記施工対象の表面に作用する前記衝撃圧実験値とを対応づけて、前記施工条件による前記ウォータジェットピーニングによって前記施工対象の表面の各位置に作用する前記衝撃圧の予測値である衝撃圧予測値を取得する予測部と、
前記予測部により取得された前記衝撃圧予測値を入力条件として、前記施工条件による前記ウォータジェットピーニングを施工後の前記施工対象における残留応力を算出する残留応力評価部と、を備える。
(11) The water jet peening evaluation apparatus according to at least one embodiment of the present invention includes:
A condition receiving unit for receiving water jet peening construction conditions for the construction object,
Analyzing the jet flow when the fluid is jetted from the nozzle model to the construction target model according to the construction conditions, at each position on the surface of the construction target model, the volume fraction of bubbles contained in the unit volume of the fluid An analysis unit for determining a void fraction, which is a volume fraction of the bubble that collapses per unit time in a unit volume of the fluid, and
An impact pressure correlation value calculation unit that determines an impact pressure correlation value that is a product of the void rate and the collapse rate at each of the positions,
An experimental value acquisition unit that acquires an impact pressure experimental value that is an experimental value of an impact pressure acting on the surface of the construction target by the water jet peening according to the construction condition.
The impact pressure correlation value at each position of the surface of the construction target model, and the impact pressure experimental value acting on the surface of the construction target by the water jet peening according to the construction condition, the impact pressure experimental value according to the construction condition A prediction unit that acquires an impact pressure predicted value that is a predicted value of the impact pressure acting on each position of the surface of the construction target by water jet peening,
A residual stress evaluation unit that calculates a residual stress in the construction target after the water jet peening under the construction conditions, using the impact pressure predicted value acquired by the prediction unit as an input condition.

上記(11)の構成によれば、ウォータジェットピーニング(WJP)におけるキャビテーション気泡の発生・消滅状況を解析することで、施工対象表面近傍に生じるWJPによる現実の衝撃圧の分布(施工対象表面の各位置での衝撃圧)を予測し、予測された衝撃圧に基づいてWJP施工後の残留応力を評価することができる。   According to the configuration of the above (11), by analyzing the occurrence and disappearance of cavitation bubbles in water jet peening (WJP), the distribution of the actual impact pressure due to WJP generated near the construction target surface (each of the construction target surface) (Impact pressure at the position), and the residual stress after WJP construction can be evaluated based on the predicted impact pressure.

(12)幾つかの実施形態では、上記(11)の構成において、
前記衝撃圧予測値は、前記施工対象モデルの表面の各位置における前記衝撃圧相関値と前記衝撃圧実験値とを対応づける係数kを決定することにより取得される。
上記(12)の構成によれば、上記の係数kを決定することで、施工対象モデルの表面の各位置における衝撃圧相関値と衝撃圧実験値とを対応づけることができ、これによって、施工対象の表面近傍におけるウォータジェットピーニング(WJP)施工後の現実の残留応力を解析・評価することができる。
(12) In some embodiments, in the configuration of the above (11),
The shock pressure predicted value is obtained by determining a coefficient k that associates the shock pressure correlation value at each position on the surface of the construction target model with the shock pressure experimental value.
According to the configuration of the above (12), by determining the coefficient k, it is possible to associate the shock pressure correlation value at each position on the surface of the construction target model with the shock pressure experimental value. It is possible to analyze and evaluate the actual residual stress after the water jet peening (WJP) in the vicinity of the target surface.

(13)幾つかの実施形態では、上記(11)〜(12)の構成において、
前記残留応力の目標値を設定する目標値設定部と、
前記残留応力が前記目標値を満たさない場合に前記施工条件を変更する施工条件変更部と、
前記残留応力が前記目標値を満たす場合に、前記残留応力の算出に用いられた施工条件を前記施工対象への前記施工条件として決定する施工条件決定部と、をさらに備え、
前記施工条件変更部により変更された変更後の施工条件により、前記解析部と前記衝撃圧相関値算出部と前記残留応力評価部とにより再評価が行われると共に、
前記再評価によって算出される残留応力が前記目標値を満たす場合には、前記施工条件決定部は、前記再評価で用いられた前記施工条件を前記施工対象への前記施工条件として決定する。
上記(13)の構成によれば、残留応力の目標値と、任意の施工条件によるウォータジェットピーニング(WJP)施工後の残留応力の解析値との比較に基づいて、この目標値を満たす施工条件が決定される。このため、プラントの仕様に合わせて実績のない施工条件によりWJPを施工するような場合でも、WJP施工後の残留応力を解析により評価することができ、プラントの仕様に応じた適切な施工条件を速やかに決定することができる。さらに、このようにして決定された施工条件によりWJPを現実に施工することで、信頼性のあるWJPの施工を行うことができる。また、所望の施工条件による施工実施が可能なWJP施工装置の設計に利用することもできる。
(13) In some embodiments, in the configuration of the above (11) to (12),
A target value setting unit for setting a target value of the residual stress,
An execution condition changing unit that changes the execution condition when the residual stress does not satisfy the target value,
When the residual stress satisfies the target value, further comprising a construction condition determining unit that determines the construction conditions used for calculating the residual stress as the construction conditions for the construction target,
The re-evaluation is performed by the analysis unit, the impact pressure correlation value calculation unit, and the residual stress evaluation unit according to the changed execution condition changed by the execution condition changing unit,
When the residual stress calculated by the reevaluation satisfies the target value, the construction condition determination unit determines the construction condition used in the reevaluation as the construction condition for the construction target.
According to the configuration of the above (13), on the basis of a comparison between the target value of the residual stress and the analysis value of the residual stress after the water jet peening (WJP) under the arbitrary processing conditions, the execution condition satisfying the target value is obtained. Is determined. For this reason, even when WJP is performed under construction conditions that have no track record in accordance with plant specifications, the residual stress after WJP construction can be evaluated by analysis, and appropriate construction conditions according to plant specifications can be determined. Can be determined promptly. Furthermore, reliable WJP construction can be performed by actually constructing WJP according to the construction conditions determined in this way. Further, the present invention can be used for designing a WJP construction apparatus capable of performing construction under desired construction conditions.

(14)幾つかの実施形態では、上記(13)の構成において、
前記再評価によって算出される残留応力が前記目標値を満たさない場合には、前記施工条件変更部によりさらに前記施工条件を変更し、前記施工条件変更部によりさらに変更された変更後の施工条件により前記再評価を行う。
上記(14)の構成によれば、任意の施工条件によるウォータジェットピーニング(WJP)施工後の残留応力の評価結果がその目標値を満たさない場合には、他の任意の施工条件による同様の評価が行われる。これによって、目標値が満たされる施工条件を決定することができる。
(14) In some embodiments, in the configuration of the above (13),
When the residual stress calculated by the re-evaluation does not satisfy the target value, further change the execution condition by the execution condition changing unit, and by the changed execution condition further changed by the execution condition changing unit. The reevaluation is performed.
According to the configuration of the above (14), when the evaluation result of the residual stress after the water jet peening (WJP) under the arbitrary construction conditions does not satisfy the target value, the similar evaluation under the other arbitrary construction conditions Is performed. As a result, it is possible to determine the construction condition that satisfies the target value.

(15)幾つかの実施形態では、上記(11)〜(14)の構成において、
前記施工条件は、前記ウォータジェットピーニングによるウォータジェットの噴射時間、前記ウォータジェットの噴射速度、前記ウォータジェットの流量、前記ウォータジェットピーニングの施工範囲、前記ウォータジェットの噴射距離、前記気泡の半径、ノズル角度、前記施工対象の表面の傾斜角度のうちの少なくとも1つを含む。
上記(15)の構成によれば、WJP施工後の残留応力に影響する可能性のある各種条件が施工条件に含まれるので、任意の施工条件よるWJP施工後の残留応力を精度よく評価することができる。
(15) In some embodiments, in the above configurations (11) to (14),
The execution conditions are: water jet injection time by the water jet peening, the water jet injection speed, the water jet flow rate, the water jet peening execution range, the water jet injection distance, the bubble radius, and the nozzle. And at least one of an angle and a tilt angle of the surface to be applied.
According to the above configuration (15), since various conditions that may affect the residual stress after WJP are included in the execution conditions, it is necessary to accurately evaluate the residual stress after the WJP is executed under any execution conditions. Can be.

本発明の少なくとも一実施形態によれば、ウォータジェットピーニング(WJP)におけるキャビテーション気泡の発生・消滅状況の解析に基づいた残留応力の評価方法が提供される。   According to at least one embodiment of the present invention, there is provided a method for evaluating residual stress based on analysis of a state of generation and disappearance of cavitation bubbles in water jet peening (WJP).

本発明の一実施形態に係るWJP評価装置の構成を概略的に示す図である。It is a figure showing roughly composition of a WJP evaluation device concerning one embodiment of the present invention. 本発明の一実施形態に係るWJPの評価方法の実施手順を示すフロー図である。It is a flowchart which shows the implementation procedure of the evaluation method of WJP which concerns on one Embodiment of this invention. 本発明の一実施形態に係る施工対象表面と噴射ノズルの関係を説明するための図である。It is a figure for explaining the relation between the construction object surface and an injection nozzle concerning one embodiment of the present invention. 本発明の一実施形態に係る噴射ノズルを示す図である。It is a figure showing an injection nozzle concerning one embodiment of the present invention. 衝撃圧相関値Pcと衝撃圧実験値Pr、衝撃圧予測値Ppの関係を図3に対応して説明するための図である。FIG. 4 is a diagram for explaining a relationship between an impact pressure correlation value Pc, an impact pressure experimental value Pr, and an impact pressure predicted value Pp, corresponding to FIG. 3. WJPによる衝撃圧の圧力分布の例を図3に対応して説明するための図である。FIG. 4 is a diagram for explaining an example of a pressure distribution of an impact pressure by WJP corresponding to FIG. 3. WJPによる衝撃圧の圧力分布の他の例を説明するための図である。It is a figure for explaining other examples of pressure distribution of the impact pressure by WJP. プラントの仕様例を説明するための図である。It is a figure for explaining the example of specification of a plant. プラントの仕様の他の例を説明するための図である。It is a figure for explaining other examples of specifications of a plant. 本発明の他の一実施形態に係るWJP評価装置の構成を概略的に示す図である。It is a figure showing roughly composition of a WJP evaluation device concerning other one embodiment of the present invention. 本発明の他の一実施形態に係るWJPの評価方法の実施手順を示すフロー図である。It is a flowchart which shows the implementation procedure of the WJP evaluation method which concerns on another embodiment of this invention.

以下、添付図面を参照して本発明の幾つかの実施形態について説明する。ただし、実施形態として記載されている又は図面に示されている構成部品の寸法、材質、形状、その相対的配置等は、本発明の範囲をこれに限定する趣旨ではなく、単なる説明例にすぎない。
例えば、「ある方向に」、「ある方向に沿って」、「平行」、「直交」、「中心」、「同心」或いは「同軸」等の相対的或いは絶対的な配置を表す表現は、厳密にそのような配置を表すのみならず、公差、若しくは、同じ機能が得られる程度の角度や距離をもって相対的に変位している状態も表すものとする。
例えば、「同一」、「等しい」及び「均質」等の物事が等しい状態であることを表す表現は、厳密に等しい状態を表すのみならず、公差、若しくは、同じ機能が得られる程度の差が存在している状態も表すものとする。
例えば、四角形状や円筒形状等の形状を表す表現は、幾何学的に厳密な意味での四角形状や円筒形状等の形状を表すのみならず、同じ効果が得られる範囲で、凹凸部や面取り部等を含む形状も表すものとする。
一方、一の構成要素を「備える」、「具える」、「具備する」、「含む」、又は、「有する」という表現は、他の構成要素の存在を除外する排他的な表現ではない。
Hereinafter, some embodiments of the present invention will be described with reference to the accompanying drawings. However, the dimensions, materials, shapes, relative arrangements, and the like of the components described in the embodiments or shown in the drawings are not intended to limit the scope of the present invention thereto, but are merely illustrative examples. Absent.
For example, expressions representing relative or absolute arrangement such as “in a certain direction”, “along a certain direction”, “parallel”, “orthogonal”, “center”, “concentric” or “coaxial” are strictly described. Not only does such an arrangement be shown, but also a state of being relatively displaced by an angle or distance that allows the same function to be obtained.
For example, expressions such as "identical", "equal", and "homogeneous", which indicate that things are in the same state, not only represent exactly the same state, but also have a tolerance or a difference that provides the same function. An existing state shall also be represented.
For example, the expression representing a shape such as a square shape or a cylindrical shape not only represents a shape such as a square shape or a cylindrical shape in a strictly geometrical sense, but also an uneven portion or a chamfer as long as the same effect can be obtained. A shape including a part and the like is also represented.
On the other hand, the expression “comprising”, “comprising”, “including”, “including”, or “having” one component is not an exclusive expression excluding the existence of another component.

図1は、本発明の一実施形態に係るウォータジェットピーニング(以下、WJPとする)の施工後の残留応力を評価するWJPの評価装置10(以下、WJP評価装置10)の機能ブロックを示す図である。
WJP評価装置10は、金属などの構造物(施工対象40)にWJPを施工する前に、WJP施工後における施工対象40の表面近傍の現実の残留応力(WJPの施工結果)を事前に評価することが可能な装置である。このWJP評価装置10は、図1に示される実施形態のように、WJP評価プログラム34を含むコンピュータであっても良い。以下では、WJP評価装置10を、WJP評価プログラム34を含むコンピュータとして説明する。
FIG. 1 is a diagram showing functional blocks of a WJP evaluation apparatus 10 (hereinafter, WJP evaluation apparatus 10) for evaluating residual stress after construction of water jet peening (hereinafter, WJP) according to an embodiment of the present invention. It is.
The WJP evaluation device 10 evaluates in advance the actual residual stress (result of WJP execution) near the surface of the execution target 40 after WJP execution, before executing WJP on a structure such as metal (construction target 40). It is a device that can do it. The WJP evaluation device 10 may be a computer including a WJP evaluation program 34 as in the embodiment shown in FIG. Hereinafter, the WJP evaluation device 10 will be described as a computer including the WJP evaluation program 34.

上述のWJP評価プログラム34を含むコンピュータは、図1に示される実施形態では、各種演算を行うCPU20と、CPU20のワークエリア等になるメモリ13(主記憶装置)と、ハードディスクドライブ装置等の補助記憶装置30とを備えている。また、コンピュータは、図1の例示のように、キーボードやマウス等の入力装置11と、表示装置(出力装置)12と、入力装置11及び表示装置12の入出力インタフェース14と、ネットワークを介して外部と通信するための通信インタフェース15と、ディスク型記憶媒体Mに対してデータの記憶処理や再生処理を行う記憶・再生装置16と、を備えていても良い。   In the embodiment shown in FIG. 1, the computer including the above-described WJP evaluation program 34 includes a CPU 20 for performing various calculations, a memory 13 (main storage device) serving as a work area for the CPU 20, and an auxiliary storage such as a hard disk drive. And an apparatus 30. In addition, as illustrated in FIG. 1, the computer includes an input device 11 such as a keyboard and a mouse, a display device (output device) 12, an input / output interface 14 of the input device 11 and the display device 12, and a network. It may include a communication interface 15 for communicating with the outside, and a storage / reproducing device 16 for performing data storage processing and data reproduction processing on the disk-type storage medium M.

補助記憶装置30には、図1の例示では、施工対象40でのWJPの施工結果を事前に評価可能な上記のWJP評価プログラム34と、OS(Operating System)プログラム37とが予め格納されている。すなわち、本実施形態のWJP評価装置10は、コンピュータにWJP評価プログラム34をインストールしたものである。WJP評価プログラム34は、図1に示される実施形態では、数値流体力学(CFD:Computational Fluid Dynamics)に基づいてWJPの施工で生じる噴流を解析する流れ解析モジュール35と、この流れ解析モジュール35での解析結果に基づいてWJPによる処理範囲を求める残留応力解析モジュール36と、を有する。これらのプログラム34,37は、例えば、記憶・再生装置16を介して、ディスク型記憶媒体Mから補助記憶装置30に取り込まれても良いし、通信インタフェース15を介して外部の装置から補助記憶装置30に取り込まれても良い。なお、WJP評価プログラム34に含まれるモジュール(35、36)の数は任意であり、他の幾つかの実施形態では、WJP評価プログラム34は1以上のモジュールにより構成されても良い。   In the example of FIG. 1, the WJP evaluation program 34 and the OS (Operating System) program 37 capable of preliminarily evaluating the WJP execution result of the execution target 40 are stored in the auxiliary storage device 30 in advance. . That is, the WJP evaluation device 10 of the present embodiment is obtained by installing the WJP evaluation program 34 on a computer. In the embodiment shown in FIG. 1, the WJP evaluation program 34 includes a flow analysis module 35 for analyzing a jet generated in WJP construction based on computational fluid dynamics (CFD), and a flow analysis module 35 for the flow analysis module 35. A residual stress analysis module 36 for obtaining a processing range by WJP based on the analysis result. These programs 34 and 37 may be loaded into the auxiliary storage device 30 from the disk-type storage medium M via the storage / reproducing device 16 or from an external device via the communication interface 15. 30. Note that the number of modules (35, 36) included in the WJP evaluation program 34 is arbitrary, and in some other embodiments, the WJP evaluation program 34 may be composed of one or more modules.

また、補助記憶装置30には、さらに、WJP評価プログラム34の実行過程で用いられる各種データが格納されても良い。図1に示される実施形態では、WJPの施工条件を規定する施工条件データ31(後述)と、施工条件データ31を構成する各種条件によってWJPを施工した場合に生じる噴流の解析結果である流れ解析データ32と、WJPの処理範囲データ33とが格納されている。   Further, the auxiliary storage device 30 may further store various data used in the execution process of the WJP evaluation program 34. In the embodiment shown in FIG. 1, a flow analysis which is an analysis result of a jet generated when WJP is executed under the execution condition data 31 (described later) which defines the execution condition of WJP and various conditions constituting the execution condition data 31. Data 32 and WJP processing range data 33 are stored.

CPU20は、機能的に、流れ解析モジュール35による噴流の解析に必要な各種パラメータを受け付けるパラメータ受付部21と、WJPの施工条件を受け付ける条件受付部22と、施工条件データ31を構成する各種条件に基づいてWJPの施工時に生じる噴流を解析する流れ解析部23と、施工対象に関する目的の処理範囲を受け付ける目的範囲受付部24と、噴流中に生じる気泡(キャビテーション気泡)のボイド率f及び崩壊率η(後述)を求めるボイド率・崩壊率算出部25と、WJPにより生じる気泡によってモデル化された施工対象40である施工対象モデル45の表面に作用すると推定される衝撃圧Pの相関値(衝撃圧相関値Pc(後述))を求める衝撃圧相関値算出部26と、上記の施工条件によるWJPによって施工対象40の表面に作用する衝撃圧Pの実験値である衝撃圧実験値Prを取得する実験値取得部27と、施工対象モデル45の表面の各位置における衝撃圧相関値Pc(衝撃圧相関値Pcの分布)と、施工条件によるWJPによって施工対象40の表面に作用する衝撃圧実験値Prとを対応づけて、施工条件によるWJPによって施工対象40の表面の各位置に作用する衝撃圧Pの予測値である衝撃圧予測値Pp(衝撃圧予測値Ppの分布)を取得する予測部28と、予測部28により取得された衝撃圧予測値Ppを入力条件として、施工条件によるWJPを施工後の施工対象40における残留応力を算出する残留応力評価部29と、を備える。   The CPU 20 functionally includes a parameter receiving unit 21 that receives various parameters necessary for the analysis of the jet by the flow analysis module 35, a condition receiving unit 22 that receives the WJP construction conditions, and various conditions that constitute the construction condition data 31. A flow analysis unit 23 that analyzes a jet generated at the time of WJP construction, a target range receiving unit 24 that receives a target processing range related to a construction target, a void rate f and a collapse rate η of bubbles (cavitation bubbles) generated in the jet. The correlation value (impact pressure) of the void rate / collapse rate calculation unit 25 for determining the later described (impact pressure P) which is estimated to act on the surface of the construction object model 45 which is the construction object 40 modeled by bubbles generated by WJP A shock pressure correlation value calculation unit 26 for calculating a correlation value Pc (described later) and construction by WJP under the above construction conditions An experimental value acquisition unit 27 that acquires an impact pressure experimental value Pr that is an experimental value of the impact pressure P acting on the surface of the elephant 40, and an impact pressure correlation value Pc (impact pressure correlation value) at each position on the surface of the construction target model 45. (Distribution of Pc) and the impact pressure experimental value Pr acting on the surface of the construction subject 40 by WJP according to the construction condition, and the impact pressure P acting on each position on the surface of the construction subject 40 by WJP according to the construction condition. A prediction unit 28 that obtains a shock pressure predicted value Pp (distribution of the shock pressure predicted value Pp), which is a predicted value, and a WJP according to a construction condition, after the impact condition predicted value Pp acquired by the prediction unit 28 is used as an input condition. And a residual stress evaluation unit 29 for calculating a residual stress in the construction target 40 of the above.

CPU20によって実行される上記の各機能部は、いずれも、上述の補助記憶装置30からメモリ13(主記憶装置)にロードされたWJP評価プログラム34をCPU20が実行することで機能する。より具体的には、図1に示される実施形態では、WJP評価プログラム34の流れ解析モジュール35の実行で、パラメータ受付部21、条件受付部22、流れ解析部23が機能する。また、WJP評価プログラム34の残留応力解析モジュール36の実行で、目的範囲受付部24、ボイド率・崩壊率算出部25、衝撃圧相関値算出部26、実験値取得部27、予測部28、残留応力評価部29が機能する。   Each of the functional units executed by the CPU 20 functions when the CPU 20 executes the WJP evaluation program 34 loaded from the auxiliary storage device 30 into the memory 13 (main storage device). More specifically, in the embodiment shown in FIG. 1, the execution of the flow analysis module 35 of the WJP evaluation program 34 causes the parameter receiving unit 21, the condition receiving unit 22, and the flow analyzing unit 23 to function. In addition, by executing the residual stress analysis module 36 of the WJP evaluation program 34, the target range receiving unit 24, the void ratio / collapse ratio calculating unit 25, the impact pressure correlation value calculating unit 26, the experimental value acquiring unit 27, the predicting unit 28, The stress evaluation unit 29 functions.

このような構成を備えるWJP評価装置10によるWJPの評価は、図2に示すフローチャートに従って実行されている。すなわち、図2は、幾つかの実施形態におけるWJPの評価方法の実行手順を示す図である。なお、以下では、説明を簡略化するため、図3に例示されるように、施工対象40の表面は平面であり、この表面の法線上(図3の例示では噴射軸Aiに一致)に位置するノズル50からこの表面上の任意の一点に向かって水を含む流体を法線に沿って噴射するものとして説明する。   The WJP evaluation by the WJP evaluation device 10 having such a configuration is executed according to the flowchart shown in FIG. That is, FIG. 2 is a diagram illustrating an execution procedure of a WJP evaluation method according to some embodiments. In the following, for simplicity of description, as illustrated in FIG. 3, the surface of the construction target 40 is a plane, and is located on a normal line of the surface (coincides with the injection axis Ai in the example of FIG. 3). A description will be given assuming that a fluid containing water is jetted from a nozzle 50 to an arbitrary point on the surface along a normal line.

図2のステップS21において、評価者は、WJP評価装置10の入力装置11を操作して、流れ解析モジュール35の各種パラメータを設定する。言い換えると、WJP評価装置10のパラメータ受付部21は、流れ解析モジュール35の各種パラメータを受け付ける(S21:プログラムのパラメータ設定工程(パラメータ受付工程))。WJP評価装置10のパラメータ受付部21は、これらのパラメータを受け付けると、メモリ13上に展開されている流れ解析モジュール35の該当箇所に、各パラメータを設定する。   In step S21 of FIG. 2, the evaluator operates the input device 11 of the WJP evaluation device 10 to set various parameters of the flow analysis module 35. In other words, the parameter receiving unit 21 of the WJP evaluation device 10 receives various parameters of the flow analysis module 35 (S21: a parameter setting step of a program (a parameter receiving step)). Upon receiving these parameters, the parameter receiving unit 21 of the WJP evaluation device 10 sets each parameter in a corresponding location of the flow analysis module 35 developed on the memory 13.

この流れ解析モジュール35は、図2に示される実施形態では、乱流を数値解析するためのLES(Large Eddy Simulation)モデルと、水と水中に存在する多数の気泡の挙動を数値解析するための2相流モデルと、気泡の発生や消滅等を含む気泡の挙動を数値解析するためのキャビテーションモデルを採用している。そして、このパラメータ設定工程(S21)では、流れ解析モジュール35の各種パラメータとして、例えば、蒸発係数、凝縮係数、気泡核生成部位体積率が設定されている。また、このパラメータ設定工程(S21)では、雰囲気気圧(施工部位の水深D)、水の密度及び粘性、飽和蒸気圧の環境(水温、圧力など)に対する依存性を示すデータが設定されており、WJP施工条件(後述)の設定によりWJPの現実の施工時におけるこれらの値と同じ値が設定されることになる。さらに、このパラメータ設定工程(S21)では蒸気の密度が設定されており、蒸気が例えば理想気体であると仮定した蒸気の密度が設定されても良い。   In the embodiment shown in FIG. 2, the flow analysis module 35 includes an LES (Large Eddy Simulation) model for numerically analyzing turbulence and a numerical analysis for behavior of water and a large number of bubbles existing in the water. A two-phase flow model and a cavitation model for numerically analyzing the behavior of bubbles including generation and disappearance of bubbles are adopted. In the parameter setting step (S21), for example, an evaporation coefficient, a condensation coefficient, and a volume ratio of a bubble nucleation site are set as various parameters of the flow analysis module 35. In this parameter setting step (S21), data indicating the atmospheric pressure (water depth D of the construction site), the density and viscosity of water, and the dependence of the saturated vapor pressure on the environment (water temperature, pressure, etc.) are set. By setting WJP execution conditions (described later), the same values as those values at the time of actual execution of WJP are set. Further, in the parameter setting step (S21), the density of the steam is set, and the density of the steam may be set assuming that the steam is, for example, an ideal gas.

次に、ステップS22において、施工対象40に対するWJPの施工条件が設定される(条件設定工程(条件受付工程))。例えば、WJP評価装置10の入力装置11が評価者によって操作されることで、WJPの施工条件がWJP評価装置10に設定される。言い換えると、WJP評価装置10の条件受付部22が、WJPの施工条件を受け付ける。この条件設定工程(S22)は、モデル設定工程(モデルの数値データ受付工程)(S22−1)と、WJP実施条件設定工程(WJP実施条件受付工程)(S22−2)とを含む。   Next, in step S22, WJP execution conditions for the execution target 40 are set (condition setting step (condition receiving step)). For example, when the input device 11 of the WJP evaluation device 10 is operated by the evaluator, the WJP execution conditions are set in the WJP evaluation device 10. In other words, the condition receiving unit 22 of the WJP evaluation device 10 receives WJP construction conditions. This condition setting step (S22) includes a model setting step (model numerical data receiving step) (S22-1) and a WJP execution condition setting step (WJP execution condition receiving step) (S22-2).

ステップS22−1におけるモデル設定工程では、条件受付部22が、施工対象40のモデルである施工対象モデル45(図3参照)が存在する空間の座標系、この施工対象モデル45を確定する数値デー夕、及び施工対象40に流体(水など)を噴射するノズル50のモデルであるノズルモデル55(図3参照)を確定する数値データを受け付ける。なお、これらの情報は、施工条件データ31の一部として補助記憶装置30に格納されても良い。   In the model setting process in step S22-1, the condition receiving unit 22 determines whether the construction target model 45 (see FIG. 3), which is a model of the construction target 40, exists in a coordinate system of a space, and numerical data for determining the construction target model 45. In the evening, numerical data for determining a nozzle model 55 (see FIG. 3), which is a model of the nozzle 50 that injects a fluid (water or the like) to the construction target 40, is received. These pieces of information may be stored in the auxiliary storage device 30 as a part of the construction condition data 31.

図3の例示では、このモデル設定工程(S22−1)において設定される上記の座標系は、ノズルモデル55の噴射軸Aiと施工対象モデル45の表面とが交差する点を原点とし、原点から噴射軸Aiが延びる方向をZ軸、このZ軸に垂直な軸をX軸、Z軸及びX軸に垂直な軸をY軸とするXYZ座標系である。また、図3では、Z軸(噴射軸Ai)上に位置するノズルモデル55からZ軸(噴射軸Ai)に沿って流体が原点方向に噴射されている。なお、XYZ座標系は、ノズルモデル55の初期位置における噴射軸Ai上の出口51oの位置を原点としてもよい。   In the example of FIG. 3, the coordinate system set in the model setting step (S22-1) has the origin at a point where the injection axis Ai of the nozzle model 55 intersects with the surface of the construction target model 45, and from the origin. This is an XYZ coordinate system in which the direction in which the injection axis Ai extends is the Z axis, the axis perpendicular to the Z axis is the X axis, and the Z axis and the axis perpendicular to the X axis are the Y axis. In FIG. 3, the fluid is ejected from the nozzle model 55 located on the Z axis (ejection axis Ai) in the direction of the origin along the Z axis (ejection axis Ai). In the XYZ coordinate system, the position of the outlet 51o on the injection axis Ai at the initial position of the nozzle model 55 may be set as the origin.

また、施工対象モデル45を確定する上述の数値データは、図3の例示では、施工対象40の表面が平面であるため、施工対象モデル45の表面が平面であることを示す情報が設定されることになる。   In the above-described numerical data for determining the construction target model 45, in the example of FIG. 3, since the surface of the construction target 40 is a plane, information indicating that the surface of the construction target model 45 is a plane is set. Will be.

ノズルモデル55を確定する上述の数値データについては、WJPの現実の施工において用いられる流体を噴射するノズル50を表現することが可能な数値データとなる。例えば、現実のノズル50は、図4に示されるように、円柱形状の部材の一方の端面から他方の端面へ貫く噴射軸Aiに沿って流路が形成されており、ノズル50の一方の端面における流路の開口が入口51iを成し、他方の端面における流路の開口が出口51oを成している。そして、この流路には、出口51o側に向かうに連れて流路径が次第に小さくなる縮径部52と、縮径部52で小さくなった流路径が維持される小径部53と、小径部53から出口に向かうに連れて次第に流路径が大きくなる拡径部54とを有している。   The above-described numerical data that determines the nozzle model 55 is numerical data that can represent the nozzle 50 that ejects a fluid used in actual WJP construction. For example, in the actual nozzle 50, as shown in FIG. 4, a flow path is formed along an injection axis Ai penetrating from one end face of the cylindrical member to the other end face, and one end face of the nozzle 50 is formed. The opening of the flow path at, forms the inlet 51i, and the opening of the flow path at the other end surface forms the outlet 51o. The flow path has a reduced diameter portion 52 whose diameter gradually decreases toward the outlet 51o, a small diameter portion 53 in which the reduced flow diameter is maintained by the reduced diameter portion 52, and a small diameter portion 53. And a diameter-enlarging portion 54 in which the flow path diameter gradually increases from the outlet to the outlet.

このようなノズル50(図4)に対しては、ノズルモデル55を確定する数値データは、例えば、ノズル50の入口51iの流路径di、小径部53の流路径ds、出口51oの流路径do、噴射軸Ai方向における縮径部52の長さNi、同方向における小径部53の長さNs、同方向における拡径部54の長さNoとしても良い。なお、本実施形態では、ノズル50の形状として以上で説明した形状のものを用いているが、他の幾つかの実施形態では、他の形状のノズルモデルを用いてもよく、この他の形状を特定可能なパラメータがモデル設定工程(S22−1)で設定される。   For such a nozzle 50 (FIG. 4), numerical data for determining the nozzle model 55 include, for example, a flow path diameter di of the inlet 51 i of the nozzle 50, a flow path diameter ds of the small diameter portion 53, and a flow path diameter do of the outlet 51 o. The length Ni of the reduced diameter portion 52 in the direction of the injection axis Ai, the length Ns of the small diameter portion 53 in the same direction, and the length No of the enlarged diameter portion 54 in the same direction may be used. In the present embodiment, the shape of the nozzle 50 described above is used as the shape of the nozzle 50. However, in some other embodiments, a nozzle model of another shape may be used. Are set in the model setting step (S22-1).

一方、ステップS22−2におけるWJP実施条件設定工程では、条件受付部22が、WJP実施条件として、WJPの施工結果に影響するよう各条件を受け付ける。このWJP実施条件の各条件は、例えば、ノズル50からの噴流の吐出圧、ノズル50からの噴流の流量、噴流の噴射距離G(ノズル50の出口51oから施工対象40の表面までの距離)(図3参照)、ノズル50の噴射軸Aiと施工対象40の表面の法線とがなす角度θ(施工対象表面の傾斜角度)、施工対象の表面の水深D、水の温度、ノズル50からの流体の噴射時間t、ノズル50からの噴射速度、流速、施工範囲S、気泡半径、施工対象の表面とノズル50の噴射軸Aiとがなす角度であるノズル角度φ、WJPの施工間隔であるピッチ等であり、WJP実施条件設定工程(S22−2)では、これらのうちの少なくとも1つの条件を受け付ける。これらを施工条件データ31の一部として補助記憶装置30に格納されても良い。なお、図2に示される実施形態では、ステップS21で設定されるプログラムのパラメータとステップS22で設定されるWJPの施工条件とを別々に設定しているが、ステップS21で設定されるパラメータもWJPの施工条件の各種条件として設定しても良い。そして、条件設定工程(S22)で設定された条件下での噴流の解析が次のステップS23で実行される(S23:解析工程)。   On the other hand, in the WJP execution condition setting step in step S22-2, the condition receiving unit 22 receives each condition as a WJP execution condition so as to affect the WJP execution result. The WJP execution conditions include, for example, the discharge pressure of the jet from the nozzle 50, the flow rate of the jet from the nozzle 50, and the jet distance G (the distance from the outlet 51o of the nozzle 50 to the surface of the construction target 40) ( (See FIG. 3), the angle θ between the injection axis Ai of the nozzle 50 and the normal of the surface of the construction target 40 (the inclination angle of the construction target surface), the water depth D of the construction target surface, the water temperature, The injection time t of the fluid, the injection speed from the nozzle 50, the flow velocity, the execution range S, the bubble radius, the nozzle angle φ which is the angle formed between the surface to be installed and the injection axis Ai of the nozzle 50, and the pitch which is the execution interval of WJP In the WJP execution condition setting step (S22-2), at least one of these conditions is accepted. These may be stored in the auxiliary storage device 30 as a part of the construction condition data 31. In the embodiment shown in FIG. 2, the program parameters set in step S21 and the WJP execution conditions set in step S22 are separately set, but the parameters set in step S21 are also set in WJP. May be set as various types of construction conditions. Then, the analysis of the jet under the conditions set in the condition setting step (S22) is executed in the next step S23 (S23: analysis step).

ステップS23における解析工程(S23)は、流れ解析工程(S23−1)と、ボイド率f及び崩壊率ηの算出工程(S23−2)を含む。
この流れ解析工程(S23−1)では、流れ解析部23が、WJPの条件設定工程(S22)で設定された条件下での噴流を解析して、施工対象モデル45の表面上の各位置での各時刻における気泡の発生数及び気泡の消滅数を求める。この解析結果は、流れ解析データ32として補助記憶装置30に格納されても良い。より詳細には、この流れ解析工程(S23−1)では、幾つかの実施形態では、施工対象モデル45の表面上で各位置での各時刻における気泡の発生数及び気泡の消滅数を以下のように求める。
The analysis step (S23) in step S23 includes a flow analysis step (S23-1) and a calculation step (S23-2) of the void rate f and the collapse rate η.
In this flow analysis step (S23-1), the flow analysis unit 23 analyzes the jet under the conditions set in the WJP condition setting step (S22), and at each position on the surface of the construction target model 45 The number of generated bubbles and the number of disappeared bubbles at each time are calculated. This analysis result may be stored in the auxiliary storage device 30 as the flow analysis data 32. More specifically, in this flow analysis step (S23-1), in some embodiments, the number of generated bubbles and the number of disappeared bubbles at each time on each position on the surface of the construction target model 45 are as follows. To ask.

すなわち、パラメータ設定工程(S21)で設定したモデル(例えば、上述の2相流モデルおよびキャビテーションモデル)を使用して、条件設定工程(S22)にて設定された条件にてCFD(例えば、非定常LES(LargeEddySimulation))による解析を実施する。解析は、通常は十分な統計情報が得られるまでの時間実施される。この解析は、例えば、ANSYS社の解析コードFLUENT等の計算機能を用いて実施しても良い。また、キャピテーションモデルについては、例えば、Philip J.Zwart,Andrew G.Gerber,Thabet Belamri:“A Two−Phase Flow Model for Predicting Cavitation.Dynamics” ICMF 2004 International Conference on Multiphase Flow,Yokohama, Japan,May 30−June 3,2004 Paper No.152にて述べられているものを使用しても良い。   That is, using the model (for example, the above-described two-phase flow model and cavitation model) set in the parameter setting step (S21), the CFD (for example, the unsteady state) is performed under the conditions set in the condition setting step (S22). Analysis by LES (Large Eddy Simulation) is performed. The analysis is usually performed until sufficient statistical information is obtained. This analysis may be performed using a calculation function such as the analysis code FLUENT of ANSYS. For the caption model, see, for example, Philip J. et al. Zwart, Andrew G .; Gerber, Thabet Belamri: "A Two-Phase Flow Model for Predicting Cavitation. Dynamics", ICMF 2004 International Conference, No. 3, January, March, January, 2010. 152 may be used.

引き続くステップS23−2では、流れ解析工程(S23−1)の完了後、ボイド率・崩壊率算出部25によって、施工対象モデル45の表面上の各位置での気泡に関するボイド率f及び崩壊率ηを算出される(S27:ボイド率f及び崩壊率ηの算出工程)。
ここで、ボイド率fとは、水を含む流体の単位体積に含まれる気泡の体積率であり、崩壊率ηとは、水を含む流体の単位体積中で単位時間に崩壊する気泡の体積率である。ボイド率・崩壊率算出部25は、補助記憶装置30あるいはメモリ13に格納されている流れ解析データ32を用いて、施工対象モデル45の表面上の各位置における流体の単位体積中で、噴射時間中における気泡の単位時間毎の体積率を求める。また、ボイド率・崩壊率算出部25は、補助記憶装置30あるいはメモリ13に格納されている流れ解析データ32を用いて、施工対象モデル45の表面上の各位置における流体の単位体積中で、各時刻における気泡の消滅数から、施工対象モデル45の表面上の各位置における流体の単位体積中で、該当時刻を含む単位時間で崩壊する気泡の体積率を求める。そして、図2に示される実施形態では、ボイド率・崩壊率算出部25は、単位時間毎の体積率の平均値をボイド率fとし、各単位時間で崩壊する気泡の体積率の平均値を崩壊率ηとして算出している。他の幾つかの実施形態では、平均値に限定されず、他の統計的手法により得られる値をボイド率fと崩壊率ηとしても良い。
In the subsequent step S23-2, after the flow analysis step (S23-1) is completed, the void rate / collapse rate calculation unit 25 calculates the void rate f and collapse rate η for bubbles at each position on the surface of the construction target model 45. Is calculated (S27: Step of calculating void ratio f and collapse ratio η).
Here, the void fraction f is the volume fraction of bubbles contained in a unit volume of a fluid containing water, and the decay rate η is the volume fraction of bubbles that collapse in a unit time in a unit volume of a fluid containing water. It is. The void rate / collapse rate calculation unit 25 uses the flow analysis data 32 stored in the auxiliary storage device 30 or the memory 13 to calculate the injection time in the unit volume of the fluid at each position on the surface of the construction target model 45. The volume ratio of air bubbles per unit time is determined. Further, the void rate / collapse rate calculation unit 25 uses the flow analysis data 32 stored in the auxiliary storage device 30 or the memory 13 to calculate, in the unit volume of the fluid at each position on the surface of the construction target model 45, From the number of disappearing bubbles at each time, the volume ratio of bubbles that collapse in a unit time including the time in the unit volume of the fluid at each position on the surface of the construction target model 45 is determined. In the embodiment illustrated in FIG. 2, the void ratio / collapse ratio calculation unit 25 sets the average value of the volume ratio per unit time as the void ratio f, and calculates the average value of the volume ratio of the bubbles that collapse in each unit time. It is calculated as the decay rate η. In some other embodiments, the values are not limited to the average value, and values obtained by other statistical methods may be used as the void rate f and the decay rate η.

ステップS24では、解析工程(S23)の完了後、衝撃圧相関値算出部26が施工対象モデル45の表面上の各位置での衝撃圧相関値をボイド率f及び崩壊率ηに基づいて算出する(S24:衝撃圧相関値算出工程)。詳述すると、WJPによる現実の衝撃圧Pは、経験則として、以下の数式に示すように、気泡の崩壊率ηと、気泡のボイド率fと、水深Dと噴流の流量に依存する係数kの積で表すことができる。
P=k×η×f
In step S24, after completion of the analysis step (S23), the impact pressure correlation value calculation unit 26 calculates an impact pressure correlation value at each position on the surface of the construction target model 45 based on the void rate f and the collapse rate η. (S24: Impact pressure correlation value calculation step). More specifically, the actual impact pressure P by WJP is, as an empirical rule, a coefficient k depending on a bubble collapse rate η, a bubble void rate f, a water depth D, and a jet flow rate as shown in the following equation. It can be expressed by the product of
P = k × η × f

そこで、上記の式中の気泡のボイド率fと気泡の崩壊率ηとの積を衝撃圧相関値Pcとして算出する。すなわち、衝撃圧相関値算出部26は、施工対象モデル45の表面上の各位置でのボイド率fと該当位置での崩壊率ηとを乗算することで、各位置での衝撃圧相関値Pc(衝撃圧相関値Pcの分布)が求められることになる。   Therefore, the product of the void ratio f of the bubble and the collapse ratio η of the bubble in the above equation is calculated as the impact pressure correlation value Pc. That is, the shock pressure correlation value calculation unit 26 multiplies the void rate f at each position on the surface of the construction target model 45 by the collapse rate η at the position to obtain the shock pressure correlation value Pc at each position. (Distribution of impact pressure correlation value Pc) is obtained.

ステップS25において、解析工程(S23)で用いられたのと同じ施工条件によるWJPによって施工対象40の表面に作用する衝撃圧Pの実験値である衝撃圧実験値Prが取得される(S25:実験値取得工程)。衝撃圧実験値Prは、事前に試験片に対してWJPを行うことで得られる衝撃圧Pの測定データ(衝撃圧実験値Pr)をWJP評価プログラム34に取り込むことで取得しても良い。衝撃圧実験値PrのWJP評価プログラム34による取り込みは、評価者によって入力されることにより行っても良いし、補助記憶装置30から読み込むことで行っても良い。
そして、上記の衝撃圧相関値算出工程(S24)と実験値取得工程(S25)により得られるデータに基づいて、予測工程(S26)が行われる。
In step S25, an impact pressure experimental value Pr, which is an experimental value of the impact pressure P acting on the surface of the construction target 40, is obtained by WJP under the same construction conditions used in the analysis step (S23) (S25: experiment). Value acquisition step). The shock pressure experimental value Pr may be acquired by loading measurement data (impact pressure experimental value Pr) of the shock pressure P obtained by performing WJP on the test piece in advance into the WJP evaluation program 34. The loading of the shock pressure experimental value Pr by the WJP evaluation program 34 may be performed by inputting by the evaluator or by reading from the auxiliary storage device 30.
Then, a prediction step (S26) is performed based on the data obtained in the shock pressure correlation value calculation step (S24) and the experimental value acquisition step (S25).

ステップS26における予測工程では、施工対象モデル45の表面の各位置における衝撃圧相関値Pcと、衝撃圧相関値Pcの算出に用いられたのと同じ施工条件によるWJPによって施工対象40の表面に作用する衝撃圧実験値Prとを対応づけて、この施工条件によるWJPによって施工対象40の表面の各位置に作用する衝撃圧Pの予測値である衝撃圧予測値Pp(衝撃圧予測値Ppの分布)が取得される。   In the prediction step in step S26, the impact pressure correlation value Pc at each position on the surface of the construction target model 45 and the effect on the surface of the construction target 40 by WJP under the same construction conditions used to calculate the shock pressure correlation value Pc. The impact pressure experimental value Pr is associated with the impact pressure predicted value Pp (the distribution of the impact pressure predicted value Pp), which is the predicted value of the impact pressure P acting on each position on the surface of the construction target 40 by WJP under the construction conditions. ) Is obtained.

より詳細には、上述した経験側を表す数式の通り、衝撃圧Pは衝撃圧相関値Pcに比例するので、実験により得られる衝撃圧Pと衝撃圧相関値Pcの相関が最も高くなるような上記の係数kを求める。すなわち、衝撃圧予測値Ppは、施工対象モデル45の表面の各位置における衝撃圧相関値Pc(衝撃圧相関値Pcの分布)と衝撃圧実験値Prとを対応づける係数kを決定することにより取得される。より具体的には、衝撃圧実験値Prが得られた施工対象表面の位置と、この位置に対応する衝撃圧相関値Pcとの差を、衝撃圧実験値Prが得られた位置のそれぞれで求め、これらの差がある範囲以内になるような係数kを求めることにより行われる。
例えば、最小二乗法によって係数kを求めても良い。
More specifically, since the impact pressure P is proportional to the impact pressure correlation value Pc according to the above-described mathematical expression representing the empirical side, the correlation between the impact pressure P and the impact pressure correlation value Pc obtained through an experiment is the highest. The above coefficient k is obtained. That is, the impact pressure prediction value Pp is determined by determining a coefficient k that associates the impact pressure correlation value Pc (distribution of the impact pressure correlation value Pc) at each position on the surface of the construction target model 45 with the impact pressure experimental value Pr. Is obtained. More specifically, the difference between the position of the construction target surface at which the shock pressure experimental value Pr was obtained and the shock pressure correlation value Pc corresponding to this position is calculated for each of the positions at which the shock pressure experimental value Pr was obtained. This is done by finding a coefficient k that makes these differences fall within a certain range.
For example, the coefficient k may be obtained by the least square method.

これらの衝撃圧相関値Pc、衝撃圧実験値Pr、衝撃圧予測値Ppの関係を図5により説明する。すなわち、図5は、衝撃圧相関値Pcと衝撃圧実験値Pr、衝撃圧予測値Ppの関係を図3に対応して説明するための図である。図5では、ノズル50の噴射軸Aiの中心(ノズル中心O)は横軸と縦軸の交点であり、ノズル中心Oからの位置は横軸で示されており、各位置における衝撃圧実験値Pr(丸印)は縦軸で示されている。また、図5には、噴射時間t(sec)が2秒(塗りつぶされた丸印)と10秒(白抜きの丸印)である異なる2つの衝撃圧実験値Prが示されており、施工対象表面における同じ位置で比較すると、噴射時間が2秒の場合に比較して10秒の場合の方が衝撃圧実験値Prは大きな値を示している。   The relationship among the shock pressure correlation value Pc, the shock pressure experimental value Pr, and the shock pressure predicted value Pp will be described with reference to FIG. That is, FIG. 5 is a diagram for explaining the relationship among the shock pressure correlation value Pc, the shock pressure experimental value Pr, and the shock pressure predicted value Pp, corresponding to FIG. In FIG. 5, the center of the injection axis Ai of the nozzle 50 (nozzle center O) is the intersection of the horizontal axis and the vertical axis, and the position from the nozzle center O is indicated by the horizontal axis. Pr (circle) is indicated by the vertical axis. FIG. 5 shows two different impact pressure experimental values Pr in which the injection time t (sec) is 2 seconds (filled circle) and 10 seconds (open circle). When compared at the same position on the target surface, the experimental value Pr of the impact pressure shows a larger value when the injection time is 10 seconds than when the injection time is 2 seconds.

一方、図5には、上述の衝撃圧相関値Pcが破線で示されているが、図5の例示では、衝撃圧相関値Pcは、衝撃圧実験値Prと比較すると、各位置において大きな値となっている。衝撃圧相関値Pcは、前述の通り、ボイド率f×崩壊圧ηであり、言い換えると、上述の数式における係数kを1とした時の値と言える。つまり、衝撃圧相関値Pcと衝撃圧実験値Prとが各位置において大きく異なるのは、衝撃圧実験値Prと衝撃圧相関値Pcを対応づけるような係数kが設定されていないことによる。   On the other hand, in FIG. 5, the above-described shock pressure correlation value Pc is indicated by a broken line, but in the example of FIG. 5, the shock pressure correlation value Pc is a large value at each position when compared with the shock pressure experimental value Pr. It has become. The impact pressure correlation value Pc is, as described above, the void ratio f × the collapse pressure η. In other words, the impact pressure correlation value Pc can be said to be a value when the coefficient k in the above equation is set to 1. That is, the reason why the shock pressure correlation value Pc and the shock pressure experimental value Pr are greatly different at each position is that the coefficient k that associates the shock pressure experimental value Pr with the shock pressure correlation value Pc is not set.

そして、上述の予測工程(S26)において両者を対応づけることで取得されたものが衝撃圧予測値Ppであり、図5の例示では、細線で示される衝撃圧予測値Ppが噴射時間2秒の衝撃圧実験値Prに対応し、太線で示される衝撃圧予測値Ppが噴射時間10秒の衝撃圧実験値Prに対応している。この二つの衝撃圧予測値Ppは、両方ともノズル中心Oから左右対称となっている。これは、図3では、平面である施工対象モデル45の法線とノズルモデル55の噴射軸Aiが一致しているなどのため、ノズル中心Oから同心円状に衝撃圧Pの等圧力線Peが分布しているためである(図6A参照)。なお、施工対象モデル45の法線とノズルモデル55が一致していない場合などでは、施工対象表面における衝撃圧Pの分布は図6Aとは異なったものになり、例えば、図6Bでは楕円状に等圧力線Peが分布している。   The shock pressure predicted value Pp obtained by associating the two with each other in the above-described prediction step (S26) is the shock pressure predicted value Pp shown in the example of FIG. Corresponding to the shock pressure experimental value Pr, the shock pressure predicted value Pp indicated by the bold line corresponds to the shock pressure experimental value Pr for the injection time of 10 seconds. Both of these two predicted impact pressure values Pp are symmetric with respect to the center O of the nozzle. This is because, in FIG. 3, since the normal line of the construction target model 45 which is a flat surface coincides with the injection axis Ai of the nozzle model 55, the isopressure line Pe of the impact pressure P is concentrically formed from the nozzle center O. This is because they are distributed (see FIG. 6A). In the case where the normal line of the construction target model 45 does not coincide with the nozzle model 55, the distribution of the impact pressure P on the construction target surface is different from that in FIG. 6A. For example, in FIG. The iso-pressure lines Pe are distributed.

ステップS27では、予測工程(S26)により取得された衝撃圧予測値Ppを入力条件として、施工条件によるWJPを施工後の施工対象40における残留応力が算出される。すなわち、衝撃圧予測値に基づいて、例えばFEM解析を行うことで、残留応力が算出される(S27:残留応力解析工程)。   In step S27, the residual stress in the construction target 40 after WJP is constructed according to the construction conditions is calculated using the impact pressure predicted value Pp acquired in the prediction step (S26) as an input condition. That is, residual stress is calculated by performing, for example, FEM analysis based on the impact pressure predicted value (S27: residual stress analysis step).

上記の実施形態によれば、ウォータジェットピーニング(WJP)におけるキャビテーション気泡の発生・消滅状況を解析することで、施工対象表面近傍に生じるWJPによる現実の衝撃圧を予測し、予測された衝撃圧に基づいてWJP施工後の残留応力を評価することができる。   According to the above embodiment, by analyzing the occurrence and disappearance of cavitation bubbles in water jet peening (WJP), the actual impact pressure due to WJP generated near the construction target surface is predicted, and the predicted impact pressure is calculated. Based on the result, the residual stress after WJP construction can be evaluated.

上記に説明した幾つかの実施形態では、設定された施工条件によるWJPにより改善される施工対象40の表面近傍の残留応力が評価されている。他の幾つかの実施形態では、このようにして評価される残留応力の評価結果に基づいて、WJPの施工条件が決定される。これは、プラントの仕様(例えば、鏡板61に溶接される管台62のサイズなど)が異なる場合には、それぞれのプラントの仕様に対応可能な信頼性のあるWJPの施工条件を迅速に決定しようとするものである。   In some embodiments described above, the residual stress near the surface of the construction target 40, which is improved by WJP under the set construction conditions, is evaluated. In some other embodiments, WJP application conditions are determined based on the evaluation results of the residual stress evaluated in this manner. This means that, when the specifications of the plant (for example, the size of the nozzle 62 to be welded to the end plate 61) are different, the reliable WJP construction conditions that can correspond to the specifications of each plant are quickly determined. It is assumed that.

例えば、図7A〜図7Bの例示では、プラントで用いられる円筒の管台62の仕様が異なっている。具体的には、図7Aに例示される管台62の直径L1は図7Bに例示される管台62の直径L2よりも小さい。また、図7Aの複数の管台62同士の間隔W1は図7Bの管台62同士の間隔W2よりも広く、これに応じてノズル50と施工表面とがなす角度φが異なっている。具体的には、図7Aでの角度φ1は、図7Bでの角度φ2よりも小さい(角度φ1<角度φ2)。また、WJPによる衝撃圧Pの大きさはノズル中心O(図5〜図6B参照)からの距離に依存するので、WJPによる残留応力の改善効果は施工範囲Sの大小に依存する。このため、所望の改善効果を得るために施工範囲Sは適切に設定されている。そして、図7Aでは、管台62の周囲に対してWJPの施工が90度ピッチで行われている。つまり、管台62の周囲に90度毎にノズル50を位置決めし、位置決めされたそれぞれの位置でWJPが施工されており、合計4回のWJPによって管台62の周囲全体が施工されている。これに対して、図7Bでは、管台62の周囲がより広く、管台62の周囲に対してWJPの施工が45度ピッチで行われている。このように施工条件が異なると、ノズル50の1ピッチ分で施工される残留応力低減有効範囲が異なってくるため、ピッチの増加やノズル角度φの調整などを行い、目標とする残留応力が満足されることを確認する必要が生じる。   For example, in the examples of FIGS. 7A and 7B, the specifications of the cylindrical nozzle 62 used in the plant are different. Specifically, the diameter L1 of the nozzle 62 illustrated in FIG. 7A is smaller than the diameter L2 of the nozzle 62 illustrated in FIG. 7B. The interval W1 between the plurality of nozzles 62 in FIG. 7A is wider than the interval W2 between the nozzles 62 in FIG. 7B, and the angle φ between the nozzle 50 and the construction surface is different accordingly. Specifically, the angle φ1 in FIG. 7A is smaller than the angle φ2 in FIG. 7B (angle φ1 <angle φ2). Further, since the magnitude of the impact pressure P by WJP depends on the distance from the nozzle center O (see FIGS. 5 to 6B), the effect of improving the residual stress by WJP depends on the size of the working range S. For this reason, the working range S is appropriately set in order to obtain a desired improvement effect. In FIG. 7A, WJP is performed around the nozzle 62 at a 90-degree pitch. That is, the nozzle 50 is positioned around the nozzle 62 every 90 degrees, and WJP is performed at each of the determined positions, and the entire circumference of the nozzle 62 is performed by a total of four WJPs. On the other hand, in FIG. 7B, the circumference of the nozzle 62 is wider, and the WJP is performed on the circumference of the nozzle 62 at a pitch of 45 degrees. If the application conditions are different as described above, the effective range of the residual stress reduction applied for one pitch of the nozzle 50 is different. Needs to be confirmed.

このため、幾つかの実施形態では図8に示されるように、WJP評価装置10は、図1に示される機能部に加えて、さらに、残留応力の目標値を設定する目標値設定部210と、残留応力が目標値を満たさない場合に施工条件を変更する施工条件変更部211と、残留応力が目標値を満たす場合に、残留応力の算出に用いられた施工条件を施工対象40への施工条件として決定する施工条件決定部212と、をさらに備える。また、これらの機能がコンピュータ(CPU20)によって実行される場合には、いずれの機能部も、上述の補助記憶装置30からメモリ13(主記憶装置)にロードされたWJP評価プログラム34をCPU20が実行することで機能する。より具体的には、図8に示される実施形態では、WJP評価プログラム34の流れ解析モジュール35の実行で、目標値設定部210、施工条件変更部211、施工条件決定部212、が機能している。
以下では、図9に示すフローチャートに従って、本実施形態におけるWJPの評価方法の実行手順について説明する。
For this reason, in some embodiments, as shown in FIG. 8, the WJP evaluation device 10 further includes a target value setting unit 210 that sets a target value of residual stress, in addition to the functional units illustrated in FIG. 1. A construction condition changing unit 211 for changing construction conditions when the residual stress does not satisfy the target value, and constructing the construction conditions used for calculating the residual stress on the construction object 40 when the residual stress satisfies the target value. And a construction condition determining unit 212 that determines the condition. When these functions are executed by the computer (CPU 20), the CPU 20 executes the WJP evaluation program 34 loaded from the auxiliary storage device 30 into the memory 13 (main storage device). It works by doing. More specifically, in the embodiment shown in FIG. 8, the execution of the flow analysis module 35 of the WJP evaluation program 34 causes the target value setting unit 210, the construction condition changing unit 211, and the construction condition determining unit 212 to function. I have.
Hereinafter, the execution procedure of the WJP evaluation method according to the present embodiment will be described with reference to the flowchart shown in FIG.

図9のステップS90において、残留応力の目標値を設定する(S90:目標値設定工程)。目標値は、例えば、施工範囲Sにおいて残留応力が0以下となることであっても良い。また、例えば、評価者により入力された目標値を受付けて、WJP評価プログラム34で利用されるよう設定される。その後、上述の残留応力の評価(図2のステップS21〜ステップS27参照)と同様の評価が行われる。つまり、図9のS91〜S97の各ステップは、それぞれ、図2のステップS21〜ステップS27と同じであるため各ステップの説明を省略する。   In step S90 of FIG. 9, a target value of the residual stress is set (S90: target value setting step). The target value may be, for example, that the residual stress is 0 or less in the construction range S. In addition, for example, the target value input by the evaluator is received and set to be used in the WJP evaluation program 34. Thereafter, the same evaluation as the above-described evaluation of the residual stress (see steps S21 to S27 in FIG. 2) is performed. That is, steps S91 to S97 in FIG. 9 are the same as steps S21 to S27 in FIG. 2, respectively, and thus description of each step is omitted.

引き続くステップS98において、算出された残留応力と、先に定められている目標値(ステップS90)とが比較される。そして、前記残留応力が前記目標値を満たさない場合には、ステップS99においてWJPの施工条件が変更される(S99:施工条件変更工程)。具体的には、WJPの施工条件設定工程(S92)で設定される各種条件のうちの少なくとも1つが変更される。例えば、前記ウォータジェットピーニングによるウォータジェットの噴射時間、前記ウォータジェットの噴射速度、前記ウォータジェットの流量、前記ウォータジェットピーニングの施工範囲S、前記ウォータジェットの噴射距離、前記気泡の半径、ノズル角度φ、前記施工対象表面の傾斜角度θのうちの少なくとも1つであっても良い。そして、施工条件変更工程(S99)によって変更された新たなWJPの施工条件に基づいて残留応力の評価を再度行われる(S91〜S97)。   In the following step S98, the calculated residual stress is compared with the previously determined target value (step S90). If the residual stress does not satisfy the target value, the WJP execution conditions are changed in step S99 (S99: execution condition changing step). Specifically, at least one of the various conditions set in the WJP construction condition setting step (S92) is changed. For example, the jet time of the water jet by the water jet peening, the jet speed of the water jet, the flow rate of the water jet, the working range S of the water jet peening, the jet distance of the water jet, the radius of the bubble, the nozzle angle φ , At least one of the inclination angles θ of the construction target surface. Then, the residual stress is evaluated again based on the new WJP execution conditions changed in the execution condition changing step (S99) (S91 to S97).

また、ステップS98において、算出された残留応力が目標値を満たす場合には、ステップS910において、その残留応力を算出するのに用いられた施工条件が、施工対象40に対する現実の施工条件として決定される。その後、フローを終了する。   If the calculated residual stress satisfies the target value in step S98, the execution condition used to calculate the residual stress is determined as the actual execution condition for the execution target 40 in step S910. You. Then, the flow ends.

なお、施工条件に対応した係数kが予め得られている場合には、実験値取得工程(S95)と、予測工程(S96)を省略し、これらに代えて、衝撃圧相関値算出工程(S94)で取得された衝撃圧相関値Pcと係数kに基づいて取得される衝撃圧予測値Ppを、残留応力解析(S97)の入力条件として直接採用しても良い。また、残留応力の評価結果に基づいて、ノズル50の機能や形状などのWJP施工装置の仕様を検討し、上記の目標値を満たす施工条件による施工を実行可能なWJP施工装置の設計に利用しても良い。   When the coefficient k corresponding to the construction condition has been obtained in advance, the experimental value acquisition step (S95) and the prediction step (S96) are omitted, and instead, the impact pressure correlation value calculation step (S94) The shock pressure prediction value Pp obtained based on the shock pressure correlation value Pc and the coefficient k obtained in the step S) may be directly adopted as an input condition for the residual stress analysis (S97). In addition, based on the evaluation results of the residual stress, the specifications of the WJP construction equipment such as the function and shape of the nozzle 50 are examined, and the design is used to design a WJP construction equipment capable of performing construction under the construction conditions satisfying the above target values. May be.

上記の実施形態によれば、プラントの仕様に合わせて実績のない施工条件によりWJPを施工するような場合でも、WJP施工後の残留応力を解析により評価することができ、プラントの仕様に応じた適切な施工条件を速やかに決定することができる。さらに、このようにして決定された施工条件によりWJPを現実に施工することで、信頼性のあるWJPの施工を行うことができる。また、所望の施工条件による施工実施が可能なWJP施工装置の設計に利用することもできる。   According to the above embodiment, even when WJP is constructed under construction conditions that have no track record in accordance with the specifications of the plant, the residual stress after WJP construction can be evaluated by analysis. Appropriate construction conditions can be determined promptly. Furthermore, reliable WJP construction can be performed by actually constructing WJP according to the construction conditions determined in this way. Further, the present invention can be used for designing a WJP construction apparatus capable of performing construction under desired construction conditions.

他の幾つかの実施形態では、WJPの施工条件に対応づけて、残留応力や、流れ解析データ32、崩壊率η、ボイド率f、衝撃圧相関値Pc、衝撃圧実験値Pr、係数k、衝撃圧予測値Ppなど関連する情報をデータベース化して利用しても良い。これによって、データベースに保持されている施工条件から関連する情報を容易に得ることができる。   In some other embodiments, the residual stress, flow analysis data 32, collapse rate η, void fraction f, impact pressure correlation value Pc, impact pressure experimental value Pr, coefficient k, Related information such as the impact pressure predicted value Pp may be stored in a database and used. This makes it possible to easily obtain related information from the construction conditions stored in the database.

本発明は上述した実施形態に限定されることはなく、上述した実施形態に変形を加えた形態や、これらの形態を適宜組み合わせた形態も含む。
The present invention is not limited to the above-described embodiment, and includes a form in which the above-described embodiment is modified and a form in which these forms are appropriately combined.

10 評価装置
11 入力装置
12 表示装置
13 メモリ
14 入出力インタフェース
15 通信インタフェース
16 記憶・再生装置

20 CPU
21 パラメータ受付部
22 条件受付部
23 流れ解析部
24 目的範囲受付部
25 ボイド率・崩壊率算出部
26 衝撃圧相関値算出部
27 実験値取得部
28 予測部
29 残留応力評価部

30 補助記憶装置
31 施工条件データ
32 流れ解析データ
33 処理範囲データ
34 評価プログラム
35 流れ解析モジュール
36 残留応力解析モジュール
37 OSプログラム

40 施工対象
45 施工対象モデル

50 ノズル
51i 入口
51o 出口
52 縮径部
53 小径部
54 拡径部
55 ノズルモデル

61 鏡板
62 管台

Ai 噴射軸
D 水深
G 噴射距離
M ディスク型記憶媒体
P 衝撃圧
Pc 衝撃圧相関値
Pe 等圧力線
Pp 衝撃圧予測値
Pr 衝撃圧実験値

di 入口流路径
do 出口流路径
ds 流路径
f ボイド率
k 係数
t 噴射時間
O ノズル中心
L1 管台の直径
L2 管台の直径
W1 管台管の間隔
W2 管台管の間隔
S 施工範囲
φ ノズル角度

Reference Signs List 10 evaluation device 11 input device 12 display device 13 memory 14 input / output interface 15 communication interface 16 storage / playback device

20 CPU
21 Parameter Receiving Unit 22 Condition Receiving Unit 23 Flow Analysis Unit 24 Target Range Receiving Unit 25 Void Fraction / Collapse Rate Calculation Unit 26 Impact Pressure Correlation Value Calculation Unit 27 Experimental Value Acquisition Unit 28 Prediction Unit 29 Residual Stress Evaluation Unit

Reference Signs List 30 auxiliary storage device 31 construction condition data 32 flow analysis data 33 processing range data 34 evaluation program 35 flow analysis module 36 residual stress analysis module 37 OS program

40 Construction target 45 Construction target model

50 Nozzle 51i Inlet 51o Outlet 52 Reduced diameter part 53 Small diameter part 54 Large diameter part 55 Nozzle model

61 head plate 62 nozzle

Ai Injection axis D Water depth G Injection distance M Disk-type storage medium P Impact pressure Pc Impact pressure correlation value Pe Isopressure line Pp Impact pressure predicted value Pr Impact pressure experimental value

di Inlet flow path diameter do Outlet flow path diameter ds Flow path diameter f Void fraction k Coefficient t Injection time O Nozzle center L1 Diameter of nozzle L2 Diameter of nozzle W1 Diameter of nozzle Nozzle W2 Distance of nozzle Nozzle S Construction range φ Nozzle angle

Claims (12)

施工対象へのウォータジェットピーニングの施工条件を設定する条件設定工程と、
ノズルモデルから施工対象モデルへ流体を噴出させたときの噴流を前記施工条件に応じて解析して、前記噴流による衝撃圧の分布が生じる前記施工対象モデルの表面の各位置での、前記流体の単位体積に含まれる気泡の体積率であるボイド率、及び前記流体の単位体積中で単位時間に崩壊する前記気泡の体積率である崩壊率を求める解析工程と、
前記各位置での前記ボイド率と前記崩壊率との積である衝撃圧相関値を求める衝撃圧相関値算出工程と、
前記施工条件による前記ウォータジェットピーニングによって現実の衝撃圧の分布が生じている際の前記施工対象または試験片の表面である実表面における少なくとも1つの位置での前記現実の衝撃圧の実験値である衝撃圧実験値を取得する実験値取得工程と、
取得された前記衝撃圧実験値と、前記施工対象モデルの表面上の前記衝撃圧の分布における位置関係が、前記衝撃圧実験値が取得された前記実表面上の位置の前記現実の衝撃圧の分布における位置関係と同じとなるような前記施工対象モデルの表面の位置における前記衝撃圧相関値と、を対応づける係数を決定することにより、前記施工条件による前記ウォータジェットピーニングによって前記施工対象の表面の各位置に作用する前記衝撃圧の予測値である衝撃圧予測値を取得する予測工程と、
前記予測工程により取得された前記衝撃圧予測値を入力条件として、前記施工条件による前記ウォータジェットピーニングを施工後の前記施工対象における残留応力を算出する残留応力評価工程と、を備えることを特徴とする残留応力評価方法。
A condition setting step of setting water jet peening construction conditions for the construction object,
Analyzing the jet flow when the fluid is jetted from the nozzle model to the construction target model according to the construction conditions, at each position on the surface of the construction target model where the distribution of the impact pressure by the jet flow occurs , the flow of the fluid Analysis step to determine the void fraction, which is the volume fraction of bubbles contained in the unit volume, and the collapse rate, which is the volume fraction of the bubbles that collapse in a unit time in the unit volume of the fluid,
Impact pressure correlation value calculation step of obtaining an impact pressure correlation value that is a product of the void rate and the collapse rate at each position,
It is an experimental value of the actual impact pressure at at least one position on the actual surface that is the surface of the application target or the test piece when the actual impact pressure distribution is generated by the water jet peening under the application condition. An experimental value acquiring step of acquiring an impact pressure experimental value,
The acquired impact pressure experimental value and the positional relationship in the distribution of the impact pressure on the surface of the construction target model are the actual impact pressure at the position on the actual surface where the impact pressure experimental value was acquired. by determining said impact pressure correlation value at position of the surface of the working object model, such as the same as the positional relationship in the distribution, Ru the association coefficients, the working object by said water jet peening by the working conditions A prediction step of obtaining a shock pressure prediction value that is a prediction value of the shock pressure acting on each position of the surface of the,
A residual stress evaluation step of calculating a residual stress in the construction target after performing the water jet peening according to the construction condition, using the impact pressure predicted value acquired in the prediction step as an input condition, Residual stress evaluation method.
前記残留応力の目標値を設定する目標値設定工程と、
前記残留応力が前記目標値を満たさない場合に前記施工条件を変更する施工条件変更工程と、
前記施工条件変更工程により変更された変更後の施工条件により前記解析工程と、前記衝撃圧相関値算出工程と、前記残留応力評価工程とを行う再評価工程と、
前記再評価工程によって算出される残留応力が前記目標値を満たす場合には、前記再評価工程で用いられた前記施工条件を前記施工対象への前記施工条件として決定する施工条件決定工程と、をさらに備えることを特徴とする請求項1に記載の残留応力評価方法。
A target value setting step of setting a target value of the residual stress,
An execution condition changing step of changing the execution condition when the residual stress does not satisfy the target value,
The analysis step, the impact pressure correlation value calculation step, and the re-evaluation step of performing the residual stress evaluation step by the changed execution condition changed by the execution condition changing step,
When the residual stress calculated by the re-evaluation step satisfies the target value, a construction condition determination step of determining the construction condition used in the re-evaluation step as the construction condition for the construction target, The residual stress evaluation method according to claim 1, further comprising:
前記再評価工程によって算出される残留応力が前記目標値を満たさない場合には、前記施工条件変更工程によりさらに前記施工条件を変更し、前記施工条件変更工程によりさらに変更された変更後の施工条件により前記再評価工程を行うことを特徴とする請求項に記載の残留応力評価方法。 If the residual stress calculated by the reevaluation step does not satisfy the target value, further change the execution condition by the execution condition change step, the changed execution condition further changed by the execution condition change step The residual stress evaluation method according to claim 2 , wherein the re-evaluation step is performed by: 前記施工条件は、前記ウォータジェットピーニングによるウォータジェットの噴射時間、前記ウォータジェットの噴射速度、前記ウォータジェットの流量、前記ウォータジェットピーニングの施工範囲、前記ウォータジェットの噴射距離、前記気泡の半径、ノズル角度、前記施工対象の表面の傾斜角度のうちの少なくとも1つを含むことを特徴とする請求項1〜のいずれか1項に記載の残留応力評価方法。 The execution conditions are: water jet injection time by the water jet peening, the water jet injection speed, the water jet flow rate, the water jet peening execution range, the water jet injection distance, the bubble radius, and the nozzle. The residual stress evaluation method according to any one of claims 1 to 3 , further comprising at least one of an angle and an inclination angle of the surface to be applied. コンピュータに、
施工対象へのウォータジェットピーニングの施工条件を設定する条件設定工程と、
ノズルモデルから施工対象モデルへ流体を噴出させたときの噴流を前記施工条件に応じて解析して、前記噴流による衝撃圧の分布が生じる前記施工対象モデルの表面の各位置での、前記流体の単位体積に含まれる気泡の体積率であるボイド率、及び前記流体の単位体積中で単位時間に崩壊する前記気泡の体積率である崩壊率を求める解析工程と、
前記各位置での前記ボイド率と前記崩壊率との積である衝撃圧相関値を求める衝撃圧相関値算出工程と、
前記施工条件による前記ウォータジェットピーニングによって現実の衝撃圧の分布が生じている際の前記施工対象または試験片の表面である実表面における少なくとも1つの位置での前記現実の衝撃圧の実験値である衝撃圧実験値を取得する実験値取得工程と、
取得された前記衝撃圧実験値と、前記施工対象モデルの表面上の前記衝撃圧の分布における位置関係が、前記衝撃圧実験値が取得された前記実表面上の位置の前記現実の衝撃圧の分布における位置関係と同じとなるような前記施工対象モデルの表面の位置における前記衝撃圧相関値と、を対応づける係数を決定することにより、前記施工条件による前記ウォータジェットピーニングによって前記施工対象の表面の各位置に作用する前記衝撃圧の予測値である衝撃圧予測値を取得する予測工程と、
前記予測工程により取得された前記衝撃圧予測値を入力条件として、前記施工条件による前記ウォータジェットピーニングを施工後の前記施工対象における残留応力を算出する残留応力評価工程と、を実行させるためのウォータジェットピーニングの評価プログラム。
On the computer,
A condition setting step of setting water jet peening construction conditions for the construction object,
Analyzing the jet flow when the fluid is jetted from the nozzle model to the construction target model according to the construction conditions, at each position on the surface of the construction target model where the distribution of the impact pressure by the jet flow occurs , the flow of the fluid Analysis step to determine the void fraction, which is the volume fraction of bubbles contained in the unit volume, and the collapse rate, which is the volume fraction of the bubbles that collapse in a unit time in the unit volume of the fluid,
Impact pressure correlation value calculation step of obtaining an impact pressure correlation value that is a product of the void rate and the collapse rate at each position,
It is an experimental value of the actual impact pressure at at least one position on the actual surface that is the surface of the application target or the test piece when the actual impact pressure distribution is generated by the water jet peening under the application condition. An experimental value acquiring step of acquiring an impact pressure experimental value,
The acquired impact pressure experimental value and the positional relationship in the distribution of the impact pressure on the surface of the construction target model are the actual impact pressure at the position on the actual surface where the impact pressure experimental value was acquired. by determining said impact pressure correlation value at position of the surface of the working object model, such as the same as the positional relationship in the distribution, Ru the association coefficients, the working object by said water jet peening by the working conditions A prediction step of obtaining a shock pressure prediction value that is a prediction value of the shock pressure acting on each position of the surface of the,
A residual stress evaluation step of calculating a residual stress in the object after the water jet peening according to the application condition, using the impact pressure predicted value obtained in the prediction step as an input condition; and Jet peening evaluation program.
前記残留応力の目標値を設定する目標値設定工程と、
前記残留応力が前記目標値を満たさない場合に前記施工条件を変更する施工条件変更工程と、
前記施工条件変更工程により変更された変更後の施工条件により前記解析工程と、前記衝撃圧相関値算出工程と、前記残留応力評価工程とを行う再評価工程と、
前記再評価工程によって算出される残留応力が前記目標値を満たす場合には、前記再評価工程で用いられた前記施工条件を前記施工対象への前記施工条件として決定する施工条件決定工程と、をさらに実行させるための請求項に記載のウォータジェットピーニングの評価プログラム。
A target value setting step of setting a target value of the residual stress,
An execution condition changing step of changing the execution condition when the residual stress does not satisfy the target value,
The analysis step, the impact pressure correlation value calculation step, and the re-evaluation step of performing the residual stress evaluation step by the changed execution condition changed by the execution condition changing step,
When the residual stress calculated by the re-evaluation step satisfies the target value, a construction condition determination step of determining the construction condition used in the re-evaluation step as the construction condition for the construction target, The water jet peening evaluation program according to claim 5 , which is further executed.
前記再評価工程によって算出される残留応力が前記目標値を満たさない場合には、前記施工条件変更工程によりさらに前記施工条件を変更し、前記施工条件変更工程によりさらに変更された変更後の施工条件により前記再評価工程を行うことを特徴とする請求項に記載のウォータジェットピーニングの評価プログラム。 If the residual stress calculated by the reevaluation step does not satisfy the target value, further change the execution condition by the execution condition change step, the changed execution condition further changed by the execution condition change step 7. The water jet peening evaluation program according to claim 6 , wherein the re-evaluation step is performed by: 前記施工条件は、前記ウォータジェットピーニングによるウォータジェットの噴射時間、前記ウォータジェットの噴射速度、前記ウォータジェットの流量、前記ウォータジェットピーニングの施工範囲、前記ウォータジェットの噴射距離、前記気泡の半径、ノズル角度、前記施工対象の表面の傾斜角度のうちの少なくとも1つを含むことを特徴とする請求項のいずれか1項に記載のウォータジェットピーニングの評価プログラム。 The execution conditions are: water jet injection time by the water jet peening, the water jet injection speed, the water jet flow rate, the water jet peening execution range, the water jet injection distance, the bubble radius, and the nozzle. The water jet peening evaluation program according to any one of claims 5 to 7 , wherein the program includes at least one of an angle and an inclination angle of a surface of the construction target. 施工対象へのウォータジェットピーニングの施工条件を受け付ける条件受付部と、
ノズルモデルから施工対象モデルへ流体を噴出させたときの噴流を前記施工条件に応じて解析して、前記噴流による衝撃圧の分布が生じる前記施工対象モデルの表面の各位置での、前記流体の単位体積に含まれる気泡の体積率であるボイド率、及び前記流体の単位体積中で単位時間に崩壊する前記気泡の体積率である崩壊率を求める解析部と、
前記各位置での前記ボイド率と前記崩壊率との積である衝撃圧相関値を求める衝撃圧相関値算出部と、
前記施工条件による前記ウォータジェットピーニングによって現実の衝撃圧の分布が生じている際の前記施工対象または試験片の表面である実表面における少なくとも1つの位置での前記現実の衝撃圧の実験値である衝撃圧実験値を取得する実験値取得部と、
取得された前記衝撃圧実験値と、前記施工対象モデルの表面上の前記衝撃圧の分布における位置関係が、前記衝撃圧実験値が取得された前記実表面上の位置の前記現実の衝撃圧の分布における位置関係と同じとなるような前記施工対象モデルの表面の位置における前記衝撃圧相関値と、を対応づける係数を決定することにより、前記施工条件による前記ウォータジェットピーニングによって前記施工対象の表面の各位置に作用する前記衝撃圧の予測値である衝撃圧予測値を取得する予測部と、
前記予測部により取得された前記衝撃圧予測値を入力条件として、前記施工条件による前記ウォータジェットピーニングを施工後の前記施工対象における残留応力を算出する残留応力評価部と、を備えることを特徴とするウォータジェットピーニングの評価装置。
A condition receiving unit for receiving water jet peening construction conditions for the construction object,
Analyzing the jet flow when the fluid is jetted from the nozzle model to the construction target model according to the construction conditions, at each position on the surface of the construction target model where the distribution of the impact pressure by the jet flow occurs , the flow of the fluid An analysis unit that determines a void fraction that is a volume fraction of bubbles contained in a unit volume, and a collapse rate that is a volume fraction of the bubbles that collapses in a unit time in a unit volume of the fluid,
An impact pressure correlation value calculation unit that determines an impact pressure correlation value that is a product of the void rate and the collapse rate at each of the positions,
It is an experimental value of the actual impact pressure at at least one position on the actual surface that is the surface of the application target or the test piece when the actual impact pressure distribution is generated by the water jet peening under the application condition. An experimental value acquiring unit for acquiring an impact pressure experimental value;
The acquired impact pressure experimental value and the positional relationship in the distribution of the impact pressure on the surface of the construction target model are the actual impact pressure at the position on the actual surface where the impact pressure experimental value was acquired. by determining said impact pressure correlation value at position of the surface of the working object model, such as the same as the positional relationship in the distribution, Ru the association coefficients, the working object by said water jet peening by the working conditions A prediction unit that obtains an impact pressure predicted value that is a predicted value of the impact pressure acting on each position on the surface of the
With the impact pressure prediction value obtained by the prediction unit as an input condition, a residual stress evaluation unit that calculates a residual stress in the construction target after performing the water jet peening under the construction condition, Water jet peening evaluation device.
前記残留応力の目標値を設定する目標値設定部と、
前記残留応力が前記目標値を満たさない場合に前記施工条件を変更する施工条件変更部と、
前記残留応力が前記目標値を満たす場合に、前記残留応力の算出に用いられた施工条件を前記施工対象への前記施工条件として決定する施工条件決定部と、をさらに備え、
前記施工条件変更部により変更された変更後の施工条件により、前記解析部と前記衝撃圧相関値算出部と前記残留応力評価部とにより再評価が行われると共に、
前記再評価によって算出される残留応力が前記目標値を満たす場合には、前記施工条件決定部は、前記再評価で用いられた前記施工条件を前記施工対象への前記施工条件として決定することを特徴とする請求項に記載のウォータジェットピーニングの評価装置。
A target value setting unit for setting a target value of the residual stress,
An execution condition changing unit that changes the execution condition when the residual stress does not satisfy the target value,
When the residual stress satisfies the target value, further comprising a construction condition determining unit that determines the construction conditions used for calculating the residual stress as the construction conditions for the construction target,
The re-evaluation is performed by the analysis unit, the impact pressure correlation value calculation unit, and the residual stress evaluation unit according to the changed execution condition changed by the execution condition changing unit,
When the residual stress calculated by the reevaluation satisfies the target value, the execution condition determining unit determines the execution condition used in the reevaluation as the execution condition for the execution target. The water jet peening evaluation device according to claim 9 , wherein:
前記再評価によって算出される残留応力が前記目標値を満たさない場合には、前記施工条件変更部によりさらに前記施工条件を変更し、前記施工条件変更部によりさらに変更された変更後の施工条件により前記再評価を行うことを特徴とする請求項10に記載のウォータジェットピーニングの評価装置。 When the residual stress calculated by the re-evaluation does not satisfy the target value, further change the execution condition by the execution condition changing unit, and by the changed execution condition further changed by the execution condition changing unit. The water jet peening evaluation apparatus according to claim 10 , wherein the reevaluation is performed. 前記施工条件は、前記ウォータジェットピーニングによるウォータジェットの噴射時間、前記ウォータジェットの噴射速度、前記ウォータジェットの流量、前記ウォータジェットピーニングの施工範囲、前記ウォータジェットの噴射距離、前記気泡の半径、ノズル角度、前記施工対象の表面の傾斜角度のうちの少なくとも1つを含むことを特徴とする請求項11のいずれか1項に記載のウォータジェットピーニングの評価装置。 The execution conditions are: water jet injection time by the water jet peening, the water jet injection speed, the water jet flow rate, the water jet peening execution range, the water jet injection distance, the bubble radius, and the nozzle. The water jet peening evaluation apparatus according to any one of claims 9 to 11 , wherein the evaluation apparatus includes at least one of an angle and an inclination angle of the surface to be applied.
JP2015021571A 2015-02-05 2015-02-05 Residual stress evaluation method Active JP6640451B2 (en)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2015021571A JP6640451B2 (en) 2015-02-05 2015-02-05 Residual stress evaluation method
PCT/JP2015/062096 WO2016125319A1 (en) 2015-02-05 2015-04-21 Residual stress evaluation method
EP15881136.4A EP3241644B1 (en) 2015-02-05 2015-04-21 Residual stress evaluation method
US15/546,077 US20180010205A1 (en) 2015-02-05 2015-04-21 Residual stress evaluation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2015021571A JP6640451B2 (en) 2015-02-05 2015-02-05 Residual stress evaluation method

Publications (2)

Publication Number Publication Date
JP2016140974A JP2016140974A (en) 2016-08-08
JP6640451B2 true JP6640451B2 (en) 2020-02-05

Family

ID=56563680

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2015021571A Active JP6640451B2 (en) 2015-02-05 2015-02-05 Residual stress evaluation method

Country Status (4)

Country Link
US (1) US20180010205A1 (en)
EP (1) EP3241644B1 (en)
JP (1) JP6640451B2 (en)
WO (1) WO2016125319A1 (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109284570A (en) * 2018-10-19 2019-01-29 燕山大学 A kind of stress assessment method and stress assessment system
CN110095226B (en) * 2019-05-30 2020-11-10 合肥工业大学 Gas impact jet flow pressure measuring device based on intermittent type spiral motion mechanism
CN110263401B (en) * 2019-06-12 2023-04-07 南京毕慕智能建筑科技有限公司 Method for evaluating residual stress relaxation effect of steel box girder top plate-longitudinal rib welding details
CN111220311B (en) * 2020-02-16 2021-05-28 山东开泰智能抛喷丸技术研究院有限公司 Detection apparatus for shot blasting dynamics is thrown to shot-blasting machine bullet
CN112560259A (en) * 2020-12-11 2021-03-26 北京航空航天大学 Transmission shaft shot blasting surface residual stress rapid prediction method based on elastic-plastic collision
TWI788873B (en) * 2021-06-08 2023-01-01 逢甲大學 Method and system for measuring interfacial stress and residual stress in multilayer thin films coated on a substrate
CN113821963B (en) * 2021-11-24 2022-03-18 武汉光谷航天三江激光产业技术研究院有限公司 Compression buckling test method and equipment for laser welding wallboard structure
CN116105905B (en) * 2023-02-03 2023-11-03 保利长大工程有限公司 Construction platform stress checking calculation system based on bridge impact drilling construction system
CN116754166B (en) * 2023-05-11 2023-12-19 中国人民解放军军事科学院系统工程研究院 Combined head and neck and trunk impact calibration test device and method

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2771202B1 (en) * 1997-11-19 2000-01-21 Inst Nat Rech Inf Automat ELECTRONIC IMAGE DATA PROCESSING DEVICE, FOR SIMULATION OF THE DEFORMABLE BEHAVIOR OF AN OBJECT
JP3803734B2 (en) * 1999-01-26 2006-08-02 株式会社日立製作所 Water jet peening equipment
US7065479B2 (en) * 2002-05-28 2006-06-20 General Electric Company Method for determining and compensating for peening-induced distortion
JP4855920B2 (en) * 2006-12-26 2012-01-18 日立Geニュークリア・エナジー株式会社 Residual stress evaluation method for water jet peening surface
US7716961B2 (en) * 2007-08-29 2010-05-18 Hitachi-Ge Nuclear Energy, Ltd. Method for executing water jet peening
JP5072885B2 (en) * 2009-03-04 2012-11-14 三菱重工業株式会社 How to set shot peening conditions
JP2011033582A (en) * 2009-08-05 2011-02-17 Toshiba Corp Residual stress evaluation method and shock load evaluation method of peened part
JP5011416B2 (en) * 2010-03-25 2012-08-29 日立Geニュークリア・エナジー株式会社 Water jet peening residual stress evaluation method and water jet peening construction method
MX348809B (en) * 2012-06-27 2017-06-28 Sintokogio Ltd Shot peening method, shot peening evaluation method, and shot peening evaluation assembly structure.
JP6169957B2 (en) * 2013-11-29 2017-07-26 三菱重工業株式会社 Water jet peening pre-evaluation method, program for executing the method, apparatus for executing the method, and water jet peening construction method

Also Published As

Publication number Publication date
JP2016140974A (en) 2016-08-08
EP3241644A4 (en) 2018-01-31
EP3241644A1 (en) 2017-11-08
EP3241644B1 (en) 2022-10-19
WO2016125319A1 (en) 2016-08-11
US20180010205A1 (en) 2018-01-11

Similar Documents

Publication Publication Date Title
JP6640451B2 (en) Residual stress evaluation method
US9573246B2 (en) Method of pre-evaluating water jet peening, recording medium having program recorded thereon to execute the method, device for executing the method, and water jet peening method
Sedlar et al. Numerical and experimental investigation of three-dimensional cavitating flow around the straight NACA2412 hydrofoil
US7912687B2 (en) Methods of predicting cavitation damage
Nguyen et al. Modeling of the bubble collapse with water jets and pressure loads using a geometrical volume of fluid based simulation method
Wu et al. Vortex Induced Motions of a Column Stabilized Floater Part II: CFD Benchmark and Prediction
Hassanien et al. On the use of surrogate models in reliability-based analysis of dented pipes
Cao et al. Mechanical behaviors of pipeline inspection gauge (pig) in launching process based on Coupled Eulerian-Lagrangian (CEL) method
Asadi et al. A computational study on droplet impingement onto a thin liquid film
Majidi et al. Towards numerical simulations of supersonic liquid jets using ghost fluid method
EP2676211B1 (en) Method of determining mechanical performance of a structure
Kurimoto et al. Response surface modeling of diesel spray parameterized by geometries inside of nozzle
JP2010020452A (en) Erosion prediction/evaluation method and device
Tufisi et al. The effect of a crack near the fixed end on the natural frequencies of a cantilever beam
JP2009128085A (en) Residual stress analysis method by finite element method
CN109635394B (en) Method and device for automatically judging safety of target material in high-energy pipeline rupture range
Peters et al. Numerical modelling and prediction of erosion induced by hydrodynamic cavitation
Kalleklev et al. Design guideline for riser collision
Doroudi Ansys Fluent modelling of an underexpanded supersonic sootblower jet impinging into recovery boiler tube geometries
Shaari et al. Fatigue crack growth analysis on square prismatic with embedded cracks under tension loading
Lefebvre et al. Failure probability assessment using co-kriging surrogate models
Hartung et al. Towards Immersed Boundary Methods For Complex Roughness Structures In Scale-Resolving Simulations.
KR20150054245A (en) Linear approximation of the model within the design space optimal design method
Zhang et al. A coupled peridynamics and DEM-IB-CLBM method for sand erosion prediction in a viscous fluid
KR102248881B1 (en) Real-time flow visualization apparatus using plant operation data and method therefor

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20180202

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20190507

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190705

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: 20191217

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20191226

R150 Certificate of patent or registration of utility model

Ref document number: 6640451

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150