WO2014033927A1 - 亀裂進展推定方法、及び情報処理装置 - Google Patents
亀裂進展推定方法、及び情報処理装置 Download PDFInfo
- Publication number
- WO2014033927A1 WO2014033927A1 PCT/JP2012/072232 JP2012072232W WO2014033927A1 WO 2014033927 A1 WO2014033927 A1 WO 2014033927A1 JP 2012072232 W JP2012072232 W JP 2012072232W WO 2014033927 A1 WO2014033927 A1 WO 2014033927A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- crack
- depth
- stress
- information processing
- processing apparatus
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M5/00—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
- G01M5/0033—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by determining damage, crack or wear
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/006—Crack, flaws, fracture or rupture
- G01N2203/0062—Crack or flaws
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/006—Crack, flaws, fracture or rupture
- G01N2203/0062—Crack or flaws
- G01N2203/0066—Propagation of crack
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Definitions
- the present invention relates to a crack propagation estimation method and an information processing apparatus, and more particularly to a technique for accurately and simply estimating the progress of a crack generated in a member.
- the plastic strain increment ⁇ p is calculated based on the strain generated in the variable load state of the equipment, and the creep strain is calculated based on the increase in strain generated in the steady load state of the equipment. It is described that the incremental ⁇ c is calculated, the fatigue damage ⁇ p is calculated using ⁇ p, the creep damage ⁇ c is calculated using ⁇ c, and the life of the device is evaluated.
- Patent Document 2 in creep crack growth evaluation, temperature and stress are analyzed from input information and temperature / stress analysis data in a database, and creep cracks are analyzed from analysis information, non-destructive data, and creep crack growth life analysis data. It describes that the progress life is analyzed and the replacement time of the high-temperature equipment target article is determined from the analysis information.
- the present invention has been made in view of such a background, and an object thereof is to provide a crack progress estimation method and an information processing apparatus capable of accurately and easily estimating the progress of a crack generated in a member.
- One of the present invention for solving the above problem is a method of estimating the progress of a crack generated in a member, and the information processing apparatus Stress distribution ⁇ (a) in the depth direction when no crack occurs for each part of the member, the relationship between the depth of the crack that propagates and the creep contribution, and the parameters C and m of the creep contribution and the Paris law
- Stress distribution ⁇ (a) in the depth direction when no crack occurs for each part of the member, the relationship between the depth of the crack that propagates and the creep contribution, and the parameters C and m of the creep contribution and the Paris law
- the creep contribution in the depth of the crack that develops for the predetermined part is obtained from the relationship between the crack depth and the creep contribution stored in the predetermined part,
- the parameters C and m corresponding to the acquired creep contribution are acquired from the relationship between the creep contribution and the Paris rule parameters C and m stored for the predetermined part,
- the information processing apparatus previously stores the stress distribution ⁇ (a) in the depth direction in the case where no crack is generated for each part of the member, the relationship between the depth of the progressing crack and the creep contribution, and the creep contribution and the Paris.
- the relationship with the parameters C and m of the law is stored, and the stress distribution ⁇ (a) in the depth direction, the creep contribution in the depth of the developing crack, and the creep contribution for a predetermined part of the member received from the user Since the corresponding parameters C and m are automatically acquired and the progress of the crack at the predetermined site is estimated, the user can obtain accurate information on the progress of the crack at the predetermined site easily and quickly.
- the information processing apparatus uses the relationship between the crack depth and the creep contribution, for example, a time-series change in stress in the depth of the crack that the member propagates, creep rupture characteristics, and cracks at the site. It is calculated based on the actual measurement value of the number of occurrences of repeated stress until it occurs.
- the information processing apparatus obtains the length of a crack generated on the surface of the predetermined portion from a user, and the depth of the crack generated in the predetermined portion based on the length of the crack.
- the crack is estimated, and the estimated depth of the crack is used as an initial value used for estimating the progress of the crack at the predetermined portion.
- the information processing apparatus automatically estimates the depth of the crack generated in the predetermined portion based on the length of the crack specified by the user. Therefore, even if the user does not input the crack depth, the crack is generated. You can get information on the progress of In consideration of empirical rules and safety of linear fracture mechanics, the information processing apparatus, for example, sets the depth of the crack to 1/3 of the length of the surface crack.
- the information processing apparatus applies a repeated stress to a curve indicating the relationship between the number N of repeated stresses and the crack length a obtained by estimating the progress of cracks in the predetermined portion.
- a tangent line is drawn from the origin to a portion that protrudes upward among the portions that change sharply. Modify the curve. Thereby, a good estimation result close to the result of numerical calculation (for example, ⁇ Jc) can be obtained.
- the numerical analysis is performed only on the stress distribution ⁇ (0) of the part in the case where the crack is not generated in the part where the crack is to be estimated, and the Paris law is used for the subsequent crack development. Therefore, it is possible to estimate the progress of a crack easily and quickly without enormous numerical calculation.
- the crack progress is estimated based on the Paris law by appropriately selecting the parameters C and m in accordance with the creep contribution in the depth of the progressing crack, it is possible to improve the estimation accuracy of the crack progress. Can do.
- the information processing apparatus obtains an actual measurement value of the time series change of the stress at the depth of the crack in which the member propagates, the creep rupture characteristics, and the number of times the repeated stress is generated until the crack is generated. Based on these, the creep contribution in the depth of the crack that develops is calculated.
- the information processing apparatus stores, for example, a relationship between a Paris law constant and a holding time, and the creep contribution degree of the relation coincides with the calculated creep contribution degree in the progressing crack depth.
- the parameters C and m are determined based on the specified relationship, and the crack propagation at the depth is estimated using the determined parameters C and m.
- FIG. 2 is a diagram illustrating a main hardware configuration of an information processing apparatus 100.
- FIG. 2 is a diagram illustrating main functions of the information processing apparatus 100.
- FIG. It is a flowchart explaining crack progress estimation processing S300. This is a screen displayed by the information processing apparatus 100 when accepting specification (input) of information from the user.
- It is an example of stress-strain characteristics.
- It is an example of stress distribution (DELTA) (sigma) (a) of the depth direction in an object site
- DELTA stress distribution
- a of the depth direction in an object site
- It is an example of the time-series change 700 of the stress in the depth of the crack which progresses.
- It is an example of the creep rupture characteristic 800.
- It is an example of the relationship 900 of a Paris law constant and retention time.
- FIG. 1 shows a main hardware configuration of an information processing apparatus 100 used for realizing a crack growth analysis system described as an embodiment.
- the crack growth analysis system of this embodiment is, for example, a high-temperature crack or fatigue generated in a member such as a structure or various equipment used at a high temperature such as a steam turbine or boiler of a power plant (thermal power plant, nuclear power plant, etc.). Used for crack analysis and diagnosis.
- the information processing apparatus 100 includes a central processing unit 101 (CPU, MPU, etc.), a main storage device 102 (ROM, RAM, NVRAM, etc.), and an auxiliary storage device 103 (hard disk drive, magneto-optical disk drive, It comprises an SSD (Solid State Drive, etc.), an input device 104 (keyboard, mouse, touch panel, etc.), an output device 105 (liquid crystal monitor, organic EL panel, etc.), and a communication device 106 (NIC (Network Interface Card, etc.)).
- CPU central processing unit
- main storage device 102 ROM, RAM, NVRAM, etc.
- auxiliary storage device 103 hard disk drive, magneto-optical disk drive, It comprises an SSD (Solid State Drive, etc.)
- an input device 104 keyboard, mouse, touch panel, etc.
- an output device 105 liquid crystal monitor, organic EL panel, etc.
- NIC Network Interface Card, etc.
- FIG. 2 shows the main functions provided by the information processing apparatus 100.
- the information processing apparatus 100 includes a stress-strain characteristic storage unit 201, a crack occurrence number storage unit 202, a stress-strain characteristic selection unit 203, a stress distribution ⁇ (0) calculation unit 204, a stress distribution ⁇ ( a)
- the functions of the calculation unit 205, the crack propagation estimation processing unit 206, the creep contribution calculation unit 207, the Paris law constant determination unit 208, the creep rupture property storage unit 209, and the Paris law constant-retention time storage unit 210 are provided. .
- the stress-strain characteristic storage unit 201 has a plurality of stresses for each stress range for one or more parts (predetermined positions of the structure) of the member for which the progress of the crack is to be estimated. Memorize strain characteristics.
- the crack occurrence number storage unit 202 acquires one or more parts (for example, a steam turbine, a boiler, which will be referred to as a target member hereinafter) that are acquired via the input device 104 and the like to estimate the progress of the crack ( For example, the number of occurrences of cracks in each of the R portion, the cutout portion, the outer peripheral portion, etc.) is stored.
- parts for example, a steam turbine, a boiler, which will be referred to as a target member hereinafter
- the stress-strain characteristic selection unit 203 performs, for each stress range stored in the stress-strain characteristic storage unit 201, for a predetermined part (hereinafter referred to as a target part) for which the progress of the crack of the target member is to be estimated.
- a target part a predetermined part for which the progress of the crack of the target member is to be estimated.
- N f is the number of cracks generated in the target region
- N cp is the number of cracks of cp type (tensile creep strain + compression plastic strain) in the strain range division method
- N pp is the pp type (tensile plastic strain in the strain range division method).
- + Compression plastic strain ⁇ cp is the cp-type strain range in the strain range division method
- ⁇ pp is the pp-type strain range in the strain range division method
- A1, A2, ⁇ 1, ⁇ 2 are all experimentally Required constants (for example, Reference 1 “Damage Analysis and Material Evaluation of Turbine Stop Valves / Control Valves in Ultra High Pressure and High Temperature Plants”, Thermal Nuclear Power Generation, Vol. 35, No. 11, Nov. 1984, Masashi Nakashiro et al. , Page 48 ").
- the stress distribution ⁇ (0) calculation unit 204 is based on the stress-strain characteristics selected by the stress-strain characteristic selection unit 203, and the stress distribution ⁇ (0) ( (Surface stress, surface stress range) is obtained and stored.
- the stress distribution ⁇ (a) calculation unit 205 performs numerical analysis (for example, analysis by the finite element method (FEM)) using ⁇ (0) obtained by the stress distribution ⁇ (0) calculation unit 204, The stress distribution ⁇ (a) in the depth direction at the target site in the case where no crack has occurred in FIG. 5 and FIG. 17 is obtained.
- FEM finite element method
- the crack propagation estimation processing unit 206 estimates the crack progress in the target portion.
- da / dN C ⁇ ( ⁇ K) m
- ⁇ K ⁇ (a) ⁇ ( ⁇ ⁇ a) 1/2
- a is the crack depth
- N is the number of occurrences of repeated stress
- C and m are constants determined according to the target member
- ⁇ K is the stress intensity factor range.
- the creep contribution calculation unit 207 obtains the relationship between the depth of the crack that propagates and the creep contribution (FIG. 18) in the above estimation performed by the crack propagation estimation processing unit 206.
- the Paris law constant determining unit 208 obtains the relationship between the depth of the crack that propagates and the creep contribution obtained by the creep contribution calculating unit 207 (FIG. 18), the creep contribution and the Paris law constants C and m described later. On the basis of the relationship (FIG. 19), the crack growth estimation processing unit 206 determines the constants C and m of the Paris rule used when estimating the growth of the crack at the depth of the crack.
- the creep rupture characteristic storage unit 209 stores a creep rupture characteristic 700 (for example, a creep rupture characteristic obtained by an experiment for the same material as the target member) described later (FIG. 6).
- the Paris law constant-holding time storage unit 210 stores a relationship 800 (FIG. 8) between the Paris law constants C and m and the holding time.
- FIG. 3 is a flowchart for explaining a process (hereinafter referred to as a crack progress estimation process S300) performed by the information processing apparatus 100 when a user estimates crack progress using a crack progress analysis system.
- a crack progress estimation process S300 a process performed by the information processing apparatus 100 when a user estimates crack progress using a crack progress analysis system.
- the crack growth estimation process S300 will be described with reference to FIG.
- the information processing apparatus 100 receives an analysis target (target member, target part, material, etc.), operation conditions (steam temperature, steam pressure, etc.), information on cracks (such as crack occurrence status) from the user via the input device 104. ) Is received (S311).
- an analysis target target member, target part, material, etc.
- operation conditions steam temperature, steam pressure, etc.
- information on cracks such as crack occurrence status
- FIG. 4 shows an example of a screen displayed on the output device 105 by the information processing apparatus 100 when receiving the above.
- the information processing apparatus 100 includes the material and the generation site as an analysis target, the steam temperature and the steam pressure as operation conditions, and the number of times of starting and stopping at the time of occurrence of a crack as information on a crack, The present surface crack length and crack depth (initial crack depth) are accepted.
- the user can omit designation (input) of the stress distribution ⁇ (0) in the target region.
- the information processing apparatus 100 automatically obtains the stress distribution ⁇ (0) by the strain range division method (the stress distribution ⁇ (0) calculation unit 204 described above). function).
- the user can omit designation (input) for the crack depth (initial crack depth).
- the information processing apparatus 100 determines the crack depth from the designated (input) length of the surface crack.
- the information processing apparatus sets the crack depth to 1/3 of the surface crack in consideration of, for example, an empirical rule of linear fracture mechanics and safety.
- the information processing apparatus 100 acquires the number N f of crack occurrences and a plurality of stress-strain characteristics for each stress range stored for the received target part (S312 and S313).
- Fig. 5 shows an example of stress-strain characteristics (extracted from Reference 1).
- the stress-strain characteristic storage unit 201 stores a plurality of stress-strain characteristics as shown in FIG.
- the information processing apparatus 100 selects a plurality of stress-strain characteristics for each stress range that satisfy the above-described relations of Expressions 1 to 3 (S314). Then, the information processing apparatus 100 obtains the stress distribution ⁇ (0) in the target part when there is no crack in the target part based on the selected stress-strain characteristics (S315).
- the information processing apparatus 100 performs numerical analysis (for example, analysis by a finite element method) using the obtained ⁇ (0), and in the case where there is no crack in the target part, A stress distribution ⁇ (a) is obtained (S316).
- FIG. 6 shows an example of the stress distribution ⁇ (a) obtained at this time.
- ⁇ (a) is obtained based on the stress distribution on the straight line in the crack depth direction at the target site, which is obtained by the numerical analysis.
- the absolute value of the stress distribution ⁇ A (a) at the time (t A ) at which the compression is maximum in the numerical analysis and the stress distribution ⁇ B at the time (t B ) at which the tension is maximum in the numerical analysis.
- the sum of the absolute value of (a) is obtained as the stress distribution ⁇ (a) of the target part when no crack occurs in the target part.
- the information processing apparatus 100 determines the progress of cracks (the number N of repeated stress occurrences) in the target region based on the obtained stress distribution ⁇ (a) and the above-described Paris rule (Equations 3 and 4). (Relationship with crack depth) is estimated (S317 to S319).
- the shape of the target member, the nature of the material, and the environment (temperature, pressure, etc.) in which the target member is placed vary depending on the depth of the crack that develops. For this reason, in the above estimation, it is considered that the degree of creep contribution also changes according to the depth of the crack that develops due to repeated stress.
- the information processing apparatus 100 obtains the creep contribution at the depth of the crack, and determines the Paris law according to the obtained creep contribution. Constants C and m are determined (selected).
- the creep contribution in the depth of the crack that propagates is the stress that actually acts on the target part (for example, it acts on a predetermined part of the steam turbine when the steam turbine is started, operated, or stopped).
- the stress rupture characteristics 800 illustrated in FIG. 8 (refer to Reference 2 “NIMS / CDS / No. 31B /”). (Extracted from p.10 of 1994) is substituted into the following equation.
- t is the time
- Nr is the number of occurrences of repeated stress applied until the target site cracks
- t r ( ⁇ (t)) is the rupture time when the stress ⁇ (t) is applied.
- ⁇ (t B ) 180 MPa
- ⁇ (t C ) 20 MPa
- ⁇ (t D ) 140 MPa
- ⁇ (t B ) 180MPa the time t a - continued 3 hours' time t B '
- ⁇ (t C ) 20MPa the time t B' - continued 150 hours of time t C
- ⁇ (t D) 140MPa the time t C -Assume that the time t E continues for 3 hours.
- the rupture time (time to rupture) corresponding to 180 MPa at time t B is 20000 hours
- the rupture time corresponding to 20 MPa at time t C is 1.0E + 07 (exponential notation) time. Since the fracture time corresponding to 140 MPa at time t D can be acquired as 50000 hours, for example, if the number of occurrences Nr of the repeated stress is 617, the creep contribution degree in the depth of the crack that develops from Equation 6 Is It becomes.
- ⁇ (t B ) 140 MPa
- ⁇ (t C ) 15 MPa
- ⁇ (t D ) 110 MPa in the time series change 700 of FIG.
- ⁇ (t C ) 15 MPa continues for 150 hours from time t B ′ ⁇ time t C
- the time to rupture corresponding to 140 MPa at time t B is 4.0E + 04 hours
- the time rupture corresponding to 20 MPa at time t C is 3.0E + 07 hours
- time Since the fracture time corresponding to 110 MPa at t D can be acquired as 1.0E + 05 hours, for example, assuming that the number of occurrences Nr of repeated stress is 617 times, the creep contribution degree in the depth of the crack that develops from Equation 6 Is It becomes.
- ⁇ (t B ) 100 MPa
- ⁇ (t C ) 10 MPa
- ⁇ (t D ) 80 MPa
- ⁇ (t B ) 100 MPa continues for 3 hours from time t A ′ ⁇ time t B ′
- ⁇ (t C ) 10 MPa continues for 150 hours from time t B ′ ⁇ time t C
- the time to rupture corresponding to 100 MPa at time t B is 1.5E + 05 hours
- the time to break corresponding to 10 MPa at time t C is 1.0E + 08 hours
- time Since the fracture time corresponding to 80 MPa at t D can be obtained as 2.5E + 05 hours, for example, if Nr is 617 times, from Equation 6, the creep contribution in the depth of the crack that develops is It becomes.
- ⁇ (t B ) 70 MPa
- ⁇ (t C ) 7 MPa
- ⁇ (t C ) 7 MPa continues for 150 hours from time t B ′ ⁇ time t C
- ⁇ (t D ) It is assumed that 56 MPa continues for 3 hours from time t C to time t E.
- the rupture time (time to rupture) corresponding to 70 MPa at time t B is 3.5E + 05 hours
- the rupture time corresponding to 7 MPa at time t C is 3.0E + 10 (exponential notation) ) time
- the breaking time corresponding to 56MPa at time t D can be acquired as 6.0E + 05 hours, for example, if Nr is 617 times, from equation 6, the creep contribution at a depth of a crack to progress , It becomes.
- the information processing apparatus 100 determines Paris constants C and m to be used when estimating the progress of a crack according to the creep contribution in the depth of the crack that has been determined in this way.
- FIG. 9 is an example of a relationship 900 between the Paris law constant and the retention time stored in the information processing apparatus 100.
- test specimens made of the same material (Cr-Mo-V cast steel) as the actual machine (for example, a steam turbine reheat steam valve bent part of a thermal power plant) are in a specified environment. (Temperature: 550 ° C., stress: ⁇ 220 MPa).
- the data held for 10,000 minutes is estimated based on the graph for holding for 1 minute to 1000 minutes.
- the information processing apparatus 100 contrasts each creep contribution degree (FIG. 9) of each graph obtained as described above with the creep contribution degree obtained from Expression 6, and among the data shown in FIG. The data corresponding to the creep contribution determined from 6 and the value is specified. In the above comparison, logarithmic interpolation is performed as necessary.
- the information processing apparatus 100 determines (selects) the constants C and m obtained from the data specified in this way as constants of the Paris rule used when estimating the progress of the crack at the depth of the crack that propagates.
- FIG. 10 shows an example in which the constants C and m of the Paris rule are determined in accordance with the degree of creep contribution at the depth of the crack that develops.
- the figure shows the Paris law constant C, m based on the 10,000 minute hold data when the crack depth is 1 mm or more and less than 3 mm, and the Paris rule constant based on the 1000 minute hold data when the crack depth is 3 mm or more and less than 5 mm.
- C and m are 5 mm or more and less than 20 mm
- the Paris rule constants C and m based on the data held for 100 minutes are used as the Paris law constants used when estimating the progress of cracks at the depth of the cracks. This is the case.
- the information processing apparatus 100 finishes estimating the progress of the crack (S319: YES), the information processing apparatus 100 outputs the result to the output apparatus 105 (S320). For example, the information processing apparatus 100 ends the estimation of the progress of a crack when the number N of occurrences of repeated stress exceeds a predetermined number.
- FIG. 11 to FIG. 13 show an example of crack propagation estimation results output by the information processing apparatus 100.
- FIG. 11 is a result of estimating the progress of a crack when the crack depth is between 1 mm and 3 mm (accuracy is not guaranteed in the range exceeding 3 mm) using Paris constants C and m based on data held for 10,000 minutes. It is. From the figure, it can be seen that, for example, when the repeated stress is generated 50 times, the depth of the crack reaches 3 mm (3.0E-03 m).
- FIG. 12 uses the Paris rule constants C and m based on the 1000 minute hold data, and the crack depth between 3 mm and 5 mm (accuracy is not guaranteed in the range of less than 3 mm and in the range of more than 5 mm). It is the result of estimating the progress. From the figure, it can be seen that, for example, the crack depth reaches 5 mm (5.0E-03 m) when the stress is generated 19 times after the crack depth exceeds 3 mm.
- FIG. 13 uses the Paris rule constants C and m based on the 100-minute hold data, and the crack depth between 5 mm and 20 mm (accuracy is not guaranteed in the range of less than 5 mm and in the range of more than 20 mm). It is the result of estimating the progress. From the figure, for example, it can be seen that the crack depth reaches 20.9 mm (20.9E-03 m) when the stress is generated 500 times after the crack depth exceeds 5 mm. Further, from the figure, it can be estimated that the growth of the crack stops when the depth of the crack is about 30.0 mm (30.0E-03 m).
- the information processing apparatus 100 sets a tangent line from the origin to the upwardly convex portion of the portion where the sharp change is seen. Then, this tangent (a line connecting the vicinity of the origin (the origin or its vicinity) and the contact point (a line indicated by a broken line in FIG. 14)) is replaced with the original progress curve to obtain a progress curve. For example, when the progress curve shown in FIG. 14 is obtained, the information processing apparatus 100 corrects the progress curve as shown in FIG. By performing such correction, it is possible to obtain a progress curve close to the actual crack progress (or the result of numerical calculation such as ⁇ Jc).
- the crack growth is estimated by appropriately selecting the Paris law constants C and m in accordance with the degree of creep contribution in the depth of the crack that progresses. It can be estimated with accuracy. Further, the Paris law constants C and m are regarded as functions C (t) and m (t) of time, and the Paris law constants C and m to be adopted in the depth of cracks that develop by combining the creep contribution and the holding time. Therefore, it is possible to estimate the progress of the crack easily and accurately using the relationship between the measured Paris law constants C and m and the retention time (the relationship between ⁇ K and da / dN).
- the information processing apparatus 100 simplifies the stress that actually acts on the member, and the crack progresses based on the time-series change of the stress at the depth of the crack that progresses. Since the relationship between the depth of the steel and the creep contribution is calculated, the time-series changes in the stress actually acting on the member can be reflected in the calculation results, and the depth of the crack that progresses in a manner close to the actual situation. The relationship with the creep contribution can be calculated. In addition, it simplifies the stress that actually acts on the member and simulates time-series changes in stress at the depth of the crack that develops, making it easier than considering the complex stress that actually acts on the member. In addition, it is possible to calculate the relationship between the depth of cracks that progress rapidly and the degree of creep contribution.
- the information processing apparatus 100 automatically obtains the stress distribution ⁇ (0) by the strain range division method when the user omits the designation of the stress distribution ⁇ (0). In addition, when the user omits the designation of the crack depth, the information processing apparatus 100 obtains the crack depth from the designated surface crack length.
- Reference numeral 152 shows the estimation results of crack propagation for each of the cases where designation of stress distribution ⁇ (0) and crack depth is omitted (reference numeral 153) and numerical calculation ( ⁇ Jc) (reference numeral 154). Show. As shown in the figure, it can be seen that even when the user omits the designation of the stress distribution ⁇ (0) and the crack depth, the progress of the crack can be estimated with high accuracy.
- the information processing apparatus 100 is obtained in advance and stored (database) the relationship between the depth of the image and the creep contribution (determined from Equation 6), and the relationship between the creep contribution and the constants C and m of the Paris rule.
- the user can more easily and quickly obtain the result of crack propagation prediction.
- the information processing apparatus 100 stores in FIG. 17 an example of the stress distribution in the depth direction when no crack is generated, and FIG. 18 shows an example of the relationship between the depth of the developing crack and the creep contribution.
- FIG. 17 An example of the relationship between the creep contribution and the Paris rule constants C and m is shown in FIG.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Artificial Intelligence (AREA)
- Data Mining & Analysis (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Computational Linguistics (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
Abstract
Description
部材の部位ごとの、亀裂が生じていない場合における深さ方向の応力分布Δσ(a)、進展する亀裂の深さとクリープ寄与度との関係、及び前記クリープ寄与度とパリス則のパラメータC,mとの関係、を記憶し、
ユーザから部材の所定の部位の指定を受け付け、
前記所定部位についての前記深さ方向の応力分布Δσ(a)を取得し、
前記所定部位についての進展する亀裂の深さにおけるクリープ寄与度を、前記所定部位について記憶している、亀裂の深さとクリープ寄与度との前記関係から取得し、
取得した前記クリープ寄与度に対応する前記パラメータC,mを、前記所定部位について記憶している、クリープ寄与度とパリス則のパラメータC,mとの前記関係から取得し、
次の関係
da/dN=C・(ΔK)m
ΔK=Δσ(a)・(π・a)1/2
(但し、aは亀裂深さ、Nは繰り返し応力の発生回数、C,mは前記部材に応じて定まる定数、ΔKは応力拡大係数範囲)
に基づき、前記所定部位における亀裂の進展を推定するものである。
前記部材の部位ごとの、亀裂が生じていない場合におけるその深さ方向の前記応力分布Δσ(a)は、
既知の複数の応力-ひずみ特性のうち次の関係を満たすものを選出し、
1/Nf=1/Npp+1/Ncp
Δεcp=A2・Ncp -α2
Δεpp=A1・Npp -α1
(但し、Nfは前記部位における既知の亀裂発生回数、Ncpはひずみ範囲分割法におけるcp型(引張クリープひずみ+圧縮塑性ひずみ)の亀裂発生回数、Nppはひずみ範囲分割法におけるpp型(引張塑性ひずみ+圧縮塑性ひずみ)の亀裂発生回数、Δεcpはひずみ範囲分割法におけるcp型のひずみ範囲、Δεppはひずみ範囲分割法におけるpp型のひずみ範囲、A1,A2,α1,α2はいずれも実験的に求められる定数)
選出した前記応力-ひずみ特性に基づき、前記部位に亀裂が生じていない場合における当該部位の応力分布Δσ(0)を求め、前記Δσ(0)に基づく数値解析を実施することにより求めたものであり、
前記情報処理装置は、
求めた前記応力分布Δσ(a)と、次の関係
da/dN=C・(ΔK)m
ΔK=Δσ(a)・(π・a)1/2
(但し、aは亀裂深さ、Nは繰り返し応力の発生回数、C,mは前記部材に応じて定まる定数、ΔKは応力拡大係数範囲)
とに基づき、前記部位における亀裂の進展を推定する。
1/Nf=1/Npp+1/Ncp ・・・ 式1
Δεcp=A2・Ncp -α2 ・・・ 式2
Δεpp=A1・Npp -α1 ・・・ 式3
da/dN=C・(ΔK)m ・・・ 式4
ΔK=Δσ(a)・(π・a)1/2 ・・・ 式5
ここでtは時刻、Nrは対象部位に亀裂が発生するまでに作用した繰り返し応力の発生回数、tr(Δσ(t))は応力Δσ(t)を作用させた場合における破断時間である。
となる。
となる。
となる。
となる。
クリープ寄与度(10000分保持)=50×(167/50000)=17%
となる。
クリープ寄与度(1000分保持)=110×(16.7/50000)=3.7%
となる。
クリープ寄与度(100分保持)=110×(1.67/50000)=1.4%
となる。
クリープ寄与度(10分保持)=1025×(0.167/50000)=0.34%
となる。
クリープ寄与度(1分保持)=2300×(0.0167/50000)=0.34%
となる。
201 応力-ひずみ特性記憶部
202 亀裂発生回数記憶部
203 応力-ひずみ特性選出部
204 応力分布Δσ(0)算出部
205 応力分布Δσ(a)算出部
206 亀裂進展推定処理部
207 クリープ寄与度算出部
208 パリス則パラメータ選択部
209 クリープ破断特性記憶部
210 パリス則定数-保持時間記憶部
S300 亀裂進展推定処理
500 進展する亀裂の深さにおける応力の時系列変化
600 クリープ破断特性
700 パリス則定数と保持時間の関係
Claims (12)
- 部材に生じる亀裂の進展を推定する方法であって、
情報処理装置が、
部材の部位ごとの、亀裂が生じていない場合における深さ方向の応力分布Δσ(a)、進展する亀裂の深さとクリープ寄与度との関係、及び前記クリープ寄与度とパリス則のパラメータC,mとの関係、を記憶し、
ユーザから部材の所定の部位の指定を受け付け、
前記所定部位についての前記深さ方向の応力分布Δσ(a)を取得し、
前記所定部位についての進展する亀裂の深さにおけるクリープ寄与度を、前記所定部位について記憶している、亀裂の深さとクリープ寄与度との前記関係から取得し、
取得した前記クリープ寄与度に対応する前記パラメータC,mを、前記所定部位について記憶している、クリープ寄与度とパリス則のパラメータC,mとの前記関係から取得し、
次の関係
da/dN=C・(ΔK)m
ΔK=Δσ(a)・(π・a)1/2
(但し、aは亀裂深さ、Nは繰り返し応力の発生回数、C,mは前記部材に応じて定まる定数、ΔKは応力拡大係数範囲)
に基づき、前記所定部位における亀裂の進展を推定する
亀裂進展推定方法。 - 請求項1に記載の亀裂進展推定方法であって、
前記情報処理装置が記憶する、亀裂の深さとクリープ寄与度との前記関係は、前記部材の進展する亀裂の深さにおける応力の時系列的な変化、クリープ破断特性、及び前記部材の部位に亀裂が発生するまでの繰り返し応力の発生回数の実測値、に基づいて算出したものである
亀裂進展推定方法。 - 請求項1又は2に記載の亀裂進展推定方法であって、
前記情報処理装置が、
ユーザから前記所定部位の表面に生じている亀裂の長さを取得し、
前記亀裂の長さに基づき前記所定部位に生じている前記亀裂の深さを推定し、
推定した前記亀裂の深さを、前記所定部位における亀裂の進展の推定に際して用いる初期値とする
亀裂進展推定方法。 - 請求項3に記載の亀裂進展推定方法であって、
前記情報処理装置が、取得した前記所定部位の表面に生じている前記亀裂の長さに1/3を乗じた値を、前記所定部位に生じている前記亀裂の深さとして推定する
亀裂進展推定方法。 - 請求項1乃至4のいずれか一項に記載の亀裂進展推定方法であって、
前記情報処理装置は、前記所定部位における亀裂の進展を推定することにより得られる、繰り返し応力の回数Nと亀裂の長さaとの関係を示す曲線に、繰り返し応力の回数Nの変化に対して亀裂の長さaが急峻に変化する部分が含まれている場合、前記急峻に変化する部分のうち上に凸になっている部分に原点近傍から接線を引くことにより前記曲線を修正する
亀裂進展推定方法。 - 請求項1乃至5のいずれか一項に記載の亀裂進展推定方法であって、
前記部材の部位ごとの、亀裂が生じていない場合におけるその深さ方向の前記応力分布Δσ(a)は、
既知の複数の応力-ひずみ特性のうち次の関係を満たすものを選出し、
1/Nf=1/Npp+1/Ncp
Δεcp=A2・Ncp -α2
Δεpp=A1・Npp -α1
(但し、Nfは前記部位における既知の亀裂発生回数、Ncpはひずみ範囲分割法におけるcp型(引張クリープひずみ+圧縮塑性ひずみ)の亀裂発生回数、Nppはひずみ範囲分割法におけるpp型(引張塑性ひずみ+圧縮塑性ひずみ)の亀裂発生回数、Δεcpはひずみ範囲分割法におけるcp型のひずみ範囲、Δεppはひずみ範囲分割法におけるpp型のひずみ範囲、A1,A2,α1,α2はいずれも実験的に求められる定数)
選出した前記応力-ひずみ特性に基づき、前記部位に亀裂が生じていない場合における当該部位の応力分布Δσ(0)を求め、前記Δσ(0)に基づく数値解析を実施することにより求めたものであり、
前記情報処理装置は、
求めた前記応力分布Δσ(a)と、次の関係
da/dN=C・(ΔK)m
ΔK=Δσ(a)・(π・a)1/2
(但し、aは亀裂深さ、Nは繰り返し応力の発生回数、C,mは前記部材に応じて定まる定数、ΔKは応力拡大係数範囲)
とに基づき、前記部位における亀裂の進展を推定する
亀裂進展推定方法。 - CPU及びメモリを供える情報処理装置であって、
部材の部位ごとの、亀裂が生じていない場合における深さ方向の応力分布Δσ(a)、進展する亀裂の深さとクリープ寄与度との関係、及び前記クリープ寄与度とパリス則のパラメータC,mとの関係、を記憶する手段と、
ユーザから部材の所定の部位の指定を受け付ける手段と、
前記所定部位についての前記深さ方向の応力分布Δσ(a)を取得する手段と、
前記所定部位についての進展する亀裂の深さにおけるクリープ寄与度を、前記所定部位について記憶している、亀裂の深さとクリープ寄与度との前記関係から取得する手段と、
取得した前記クリープ寄与度に対応する前記パラメータC,mを、前記所定部位について記憶している、クリープ寄与度とパリス則のパラメータC,mとの前記関係から取得する手段と、
次の関係
da/dN=C・(ΔK)m
ΔK=Δσ(a)・(π・a)1/2
(但し、aは亀裂深さ、Nは繰り返し応力の発生回数、C,mは前記部材に応じて定まる定数、ΔKは応力拡大係数範囲)
に基づき、前記所定部位における亀裂の進展を推定する手段と
を備える情報処理装置。 - 請求項7に記載の情報処理装置であって、
亀裂の深さとクリープ寄与度との前記関係は、前記部材の進展する亀裂の深さにおける応力の時系列的な変化、クリープ破断特性、及び前記部位に亀裂が発生するまでの繰り返し応力の発生回数の実測値、に基づいて算出したものである
情報処理装置。 - 請求項7又は8に記載の情報処理装置であって、
ユーザから前記所定部位の表面に生じている亀裂の長さを取得する手段と、
前記亀裂の長さに基づき前記所定部位に生じている前記亀裂の深さを推定する手段と、
推定した前記亀裂の深さを、前記所定部位における亀裂の進展の推定に際して用いる初期値とする手段と
を備える情報処理装置。 - 請求項9に記載の情報処理装置であって、
取得した前記所定部位の表面に生じている前記亀裂の長さに1/3を乗じた値を、前記所定部位に生じている前記亀裂の深さとして推定する
情報処理装置。 - 請求項7乃至10のいずれか一項に記載の情報処理装置であって、
前記所定部位における亀裂の進展を推定することにより得られる、繰り返し応力の回数Nと亀裂の長さaとの関係を示す曲線に、繰り返し応力の回数Nの変化に対して亀裂の長さaが急峻に変化する部分が含まれている場合、前記急峻に変化する部分のうち上に凸になっている部分に原点近傍から接線を引くことにより前記曲線を修正する手段を備える
情報処理装置。 - 請求項7乃至11のいずれか一項に記載の情報処理装置であって、
前記部材の部位ごとの、亀裂が生じていない場合におけるその深さ方向の前記応力分布Δσ(a)を、
既知の複数の応力-ひずみ特性のうち次の関係を満たすものを選出し、
1/Nf=1/Npp+1/Ncp
Δεcp=A2・Ncp -α2
Δεpp=A1・Npp -α1
(但し、Nfは前記部位における既知の亀裂発生回数、Ncpはひずみ範囲分割法におけるcp型(引張クリープひずみ+圧縮塑性ひずみ)の亀裂発生回数、Nppはひずみ範囲分割法におけるpp型(引張塑性ひずみ+圧縮塑性ひずみ)の亀裂発生回数、Δεcpはひずみ範囲分割法におけるcp型のひずみ範囲、Δεppはひずみ範囲分割法におけるpp型のひずみ範囲、A1,A2,α1,α2はいずれも実験的に求められる定数)
選出した前記応力-ひずみ特性に基づき、前記部位に亀裂が生じていない場合における当該部位の応力分布Δσ(0)を求め、前記Δσ(0)に基づく数値解析を実施することにより求める手段と、
求めた前記応力分布Δσ(a)と、次の関係(パリス則)
da/dN=C・(ΔK)m
ΔK=Δσ(a)・(π・a)1/2
(但し、aは亀裂深さ、Nは繰り返し応力の発生回数、C,mは前記部材に応じて定まる定数、ΔKは応力拡大係数範囲)
とに基づき、前記部位における亀裂の進展を推定する手段と
を備える情報処理装置。
Priority Applications (10)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2012/072232 WO2014033927A1 (ja) | 2012-08-31 | 2012-08-31 | 亀裂進展推定方法、及び情報処理装置 |
EP12883660.8A EP2891875B1 (en) | 2012-08-31 | 2012-08-31 | Fissure progress estimation method and information processing device |
JP2014505434A JP5567233B1 (ja) | 2012-08-31 | 2012-08-31 | 亀裂進展推定方法、及び情報処理装置 |
CA2885718A CA2885718A1 (en) | 2012-08-31 | 2012-08-31 | Method of predicting crack growth and information processing device therefor |
US14/424,076 US10275546B2 (en) | 2012-08-31 | 2012-08-31 | Method of predicting crack growth and information processing device therefor |
BR112015004558A BR112015004558A2 (pt) | 2012-08-31 | 2012-08-31 | Método de previsão de crescimento de fissura e dispositivo de processamento de informação do mesmo |
KR1020157005845A KR20150054794A (ko) | 2012-08-31 | 2012-08-31 | 균열 진전 추정 방법 및 정보 처리 장치 |
MX2015002685A MX345386B (es) | 2012-08-31 | 2012-08-31 | Metodo de procesamiento de información y dispositivo de procesamiento de informacion. |
CN201280075936.XA CN104685337A (zh) | 2012-08-31 | 2012-08-31 | 龟裂进展推定方法以及信息处理装置 |
PL12883660T PL2891875T3 (pl) | 2012-08-31 | 2012-08-31 | Sposób oceny postępu pęknięcia i urządzenie przetwarzające informacje |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2012/072232 WO2014033927A1 (ja) | 2012-08-31 | 2012-08-31 | 亀裂進展推定方法、及び情報処理装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2014033927A1 true WO2014033927A1 (ja) | 2014-03-06 |
Family
ID=50182773
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2012/072232 WO2014033927A1 (ja) | 2012-08-31 | 2012-08-31 | 亀裂進展推定方法、及び情報処理装置 |
Country Status (10)
Country | Link |
---|---|
US (1) | US10275546B2 (ja) |
EP (1) | EP2891875B1 (ja) |
JP (1) | JP5567233B1 (ja) |
KR (1) | KR20150054794A (ja) |
CN (1) | CN104685337A (ja) |
BR (1) | BR112015004558A2 (ja) |
CA (1) | CA2885718A1 (ja) |
MX (1) | MX345386B (ja) |
PL (1) | PL2891875T3 (ja) |
WO (1) | WO2014033927A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016194433A (ja) * | 2015-03-31 | 2016-11-17 | みずほ情報総研株式会社 | 損傷評価システム及び損傷評価プログラム |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5567233B1 (ja) | 2012-08-31 | 2014-08-06 | 中国電力株式会社 | 亀裂進展推定方法、及び情報処理装置 |
US9792555B2 (en) * | 2013-01-04 | 2017-10-17 | Siemens Energy, Inc. | Probabilistic modeling and sizing of embedded flaws in ultrasonic nondestructive inspections for fatigue damage prognostics and structural integrity assessment |
WO2017051552A1 (ja) * | 2015-09-25 | 2017-03-30 | 富士フイルム株式会社 | 点検計画立案支援システム、方法およびプログラム |
JP6974984B2 (ja) * | 2017-08-30 | 2021-12-01 | 三菱パワー株式会社 | 余寿命評価方法及び保守管理方法 |
JP6979307B2 (ja) * | 2017-08-30 | 2021-12-08 | 三菱パワー株式会社 | 亀裂評価基準策定方法、内部探傷検査による亀裂評価方法及び保守管理方法 |
EP3734547B1 (en) * | 2017-12-25 | 2022-01-05 | Fujitsu Limited | Image processing program, image processing method, and image processing device |
CN108415966B (zh) * | 2018-02-07 | 2021-09-24 | 中国地质大学(武汉) | 一种基于水箱模型的岩溶流域水文过程模拟方法 |
CN110411833B (zh) * | 2019-08-12 | 2021-09-14 | 武汉钢铁有限公司 | 一种海水腐蚀环境下不同频率裂纹扩展速率预测方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006105673A (ja) * | 2004-10-01 | 2006-04-20 | Toshiba Corp | 表面き裂の進展解析方法および装置 |
JP2007256042A (ja) * | 2006-03-23 | 2007-10-04 | Mitsubishi Heavy Ind Ltd | ガスタービン高温部品のき裂進展予測方法及びこの方法を用いたき裂進展予測装置 |
JP2008003009A (ja) | 2006-06-23 | 2008-01-10 | Toshiba Corp | 高温機器の寿命診断装置、高温機器の寿命診断方法およびプログラム |
JP2009031124A (ja) * | 2007-07-27 | 2009-02-12 | Mitsubishi Heavy Ind Ltd | き裂進展予測方法及びプログラム |
JP2010216983A (ja) | 2009-03-17 | 2010-09-30 | Toshiba Corp | 機器寿命評価システムおよび機器寿命評価方法 |
JP2011095178A (ja) * | 2009-10-30 | 2011-05-12 | Itochu Techno-Solutions Corp | き裂進展解析方法及びプログラム |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7684936B2 (en) | 2008-05-07 | 2010-03-23 | Simmonds Precision Products, Inc. | Method, apparatus and computer program product for determining a future time of a component |
EP2386090B1 (en) * | 2009-01-08 | 2019-07-03 | Battelle Memorial Institute | Path-dependent cycle counting and multi-axial fatigue evaluation of engineering structures |
US9581570B2 (en) * | 2011-02-10 | 2017-02-28 | University Of South Carolina | Determination of the remaining life of a structural system based on acoustic emission signals |
JP5567233B1 (ja) | 2012-08-31 | 2014-08-06 | 中国電力株式会社 | 亀裂進展推定方法、及び情報処理装置 |
-
2012
- 2012-08-31 JP JP2014505434A patent/JP5567233B1/ja not_active Expired - Fee Related
- 2012-08-31 US US14/424,076 patent/US10275546B2/en active Active
- 2012-08-31 MX MX2015002685A patent/MX345386B/es active IP Right Grant
- 2012-08-31 CN CN201280075936.XA patent/CN104685337A/zh active Pending
- 2012-08-31 EP EP12883660.8A patent/EP2891875B1/en active Active
- 2012-08-31 WO PCT/JP2012/072232 patent/WO2014033927A1/ja active Application Filing
- 2012-08-31 PL PL12883660T patent/PL2891875T3/pl unknown
- 2012-08-31 CA CA2885718A patent/CA2885718A1/en not_active Abandoned
- 2012-08-31 BR BR112015004558A patent/BR112015004558A2/pt not_active IP Right Cessation
- 2012-08-31 KR KR1020157005845A patent/KR20150054794A/ko active IP Right Grant
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006105673A (ja) * | 2004-10-01 | 2006-04-20 | Toshiba Corp | 表面き裂の進展解析方法および装置 |
JP2007256042A (ja) * | 2006-03-23 | 2007-10-04 | Mitsubishi Heavy Ind Ltd | ガスタービン高温部品のき裂進展予測方法及びこの方法を用いたき裂進展予測装置 |
JP2008003009A (ja) | 2006-06-23 | 2008-01-10 | Toshiba Corp | 高温機器の寿命診断装置、高温機器の寿命診断方法およびプログラム |
JP2009031124A (ja) * | 2007-07-27 | 2009-02-12 | Mitsubishi Heavy Ind Ltd | き裂進展予測方法及びプログラム |
JP2010216983A (ja) | 2009-03-17 | 2010-09-30 | Toshiba Corp | 機器寿命評価システムおよび機器寿命評価方法 |
JP2011095178A (ja) * | 2009-10-30 | 2011-05-12 | Itochu Techno-Solutions Corp | き裂進展解析方法及びプログラム |
Non-Patent Citations (3)
Title |
---|
KATSUYUKI SHIBATA ET AL.: "Development of RPV Reliability Analysis Code Based on Probabilistic Fracture Mechanics Methodology", JOURNAL OF THE ATOMIC ENERGY SOCIETY OF JAPAN, vol. 43, no. 4, 2001, pages 387 - 396, XP055253801 * |
MASASHI NAKASHIRO ET AL.: "Failure analysis and metallurgical evaluation of turbine stop and control valves used at super high temperature and pressure plant", THERM. NUCL. POWER, vol. 35, no. 11, November 1984 (1984-11-01), pages 48 |
See also references of EP2891875A4 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016194433A (ja) * | 2015-03-31 | 2016-11-17 | みずほ情報総研株式会社 | 損傷評価システム及び損傷評価プログラム |
Also Published As
Publication number | Publication date |
---|---|
JPWO2014033927A1 (ja) | 2016-08-08 |
US10275546B2 (en) | 2019-04-30 |
EP2891875A1 (en) | 2015-07-08 |
PL2891875T3 (pl) | 2018-10-31 |
MX345386B (es) | 2017-01-26 |
BR112015004558A2 (pt) | 2017-08-22 |
KR20150054794A (ko) | 2015-05-20 |
MX2015002685A (es) | 2015-11-23 |
CA2885718A1 (en) | 2014-03-06 |
JP5567233B1 (ja) | 2014-08-06 |
US20150324697A1 (en) | 2015-11-12 |
CN104685337A (zh) | 2015-06-03 |
EP2891875A4 (en) | 2016-04-06 |
EP2891875B1 (en) | 2018-06-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5567233B1 (ja) | 亀裂進展推定方法、及び情報処理装置 | |
US8532941B2 (en) | Fatigue life estimation method and system | |
Zhu et al. | Bayesian framework for probabilistic low cycle fatigue life prediction and uncertainty modeling of aircraft turbine disk alloys | |
Agius et al. | Sensitivity and optimisation of the Chaboche plasticity model parameters in strain-life fatigue predictions | |
JP2008003009A (ja) | 高温機器の寿命診断装置、高温機器の寿命診断方法およびプログラム | |
JP2015532430A (ja) | 確率論的疲労亀裂寿命推定のための方法およびシステム | |
Herz et al. | Numerical simulation of plasticity induced fatigue crack opening and closure for autofrettaged intersecting holes | |
Šeruga et al. | Unification of the most commonly used time–temperature creep parameters | |
EP2724141A2 (en) | System and method for generating a combined model for isothermal and anisothermal fatigue life | |
US20130298691A1 (en) | Part fatigue fracture evaluating apparatus, part fatigue fracture evaluating method, and computer program | |
Guo et al. | A damage coupled elastic-plastic constitutive model and its application on low cycle fatigue life prediction of turbine blade | |
EP3819806A1 (en) | System, apparatus and method for estimating life of components | |
Shen et al. | Analysis on the fatigue damage evolution of notched specimens with consideration of cyclic plasticity | |
Semenov et al. | Thermal fatigue fracture criteria of single crystal heat-resistant alloys and methods for identification of their parameters | |
Shlyannikov et al. | Creep‐fracture resistance parameters determination based on stress and ductility damage models | |
Beesley et al. | A novel simulation for the design of a low cycle fatigue experimental testing programme | |
Tao et al. | Life prediction based on weight‐averaged maximum shear strain range plane under multiaxial variable amplitude loading | |
Hu et al. | Investigation of low cycle fatigue crack propagation behavior of 316H steel at 550℃ based on cyclic response and damage accumulation: experiment and modelling | |
JP5121593B2 (ja) | 機器寿命評価方法および機器寿命評価システム | |
JP5567081B2 (ja) | 亀裂進展推定方法、及び情報処理装置 | |
WO2013114604A1 (ja) | 亀裂進展推定方法、及び情報処理装置 | |
Guan | Probabilistic modeling of threshold stress intensity factor for fatigue endurance reliability prediction | |
JP2013096862A (ja) | き裂進展挙動予測方法 | |
JP5507777B2 (ja) | 亀裂進展推定方法、及び情報処理装置 | |
Cai et al. | A combined fatigue-ductile damage model with optimized parameters for predicting the mechanical response of high-Cr steels at high stress levels |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
ENP | Entry into the national phase |
Ref document number: 2014505434 Country of ref document: JP Kind code of ref document: A |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 12883660 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: MX/A/2015/002685 Country of ref document: MX |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
ENP | Entry into the national phase |
Ref document number: 20157005845 Country of ref document: KR Kind code of ref document: A |
|
ENP | Entry into the national phase |
Ref document number: 2885718 Country of ref document: CA |
|
WWE | Wipo information: entry into national phase |
Ref document number: 14424076 Country of ref document: US |
|
REG | Reference to national code |
Ref country code: BR Ref legal event code: B01A Ref document number: 112015004558 Country of ref document: BR |
|
ENP | Entry into the national phase |
Ref document number: 112015004558 Country of ref document: BR Kind code of ref document: A2 Effective date: 20150226 |