WO2019208061A1 - 残留応力算出方法 - Google Patents

残留応力算出方法 Download PDF

Info

Publication number
WO2019208061A1
WO2019208061A1 PCT/JP2019/012534 JP2019012534W WO2019208061A1 WO 2019208061 A1 WO2019208061 A1 WO 2019208061A1 JP 2019012534 W JP2019012534 W JP 2019012534W WO 2019208061 A1 WO2019208061 A1 WO 2019208061A1
Authority
WO
WIPO (PCT)
Prior art keywords
residual stress
metal material
finite element
strain
calculated value
Prior art date
Application number
PCT/JP2019/012534
Other languages
English (en)
French (fr)
Inventor
真理子 松田
圭介 沖田
知和 中川
Original Assignee
株式会社神戸製鋼所
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 株式会社神戸製鋼所 filed Critical 株式会社神戸製鋼所
Priority to CN201980027485.4A priority Critical patent/CN112005088B/zh
Priority to FIEP19791476.5T priority patent/FI3786600T3/fi
Priority to EP19791476.5A priority patent/EP3786600B1/en
Priority to ES19791476T priority patent/ES2941730T3/es
Priority to KR1020207033511A priority patent/KR102475949B1/ko
Publication of WO2019208061A1 publication Critical patent/WO2019208061A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/20—Metals
    • G01N33/204—Structure thereof, e.g. crystal structure
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L1/00—Measuring force or stress, in general
    • G01L1/25—Measuring force or stress, in general using wave or particle radiation, e.g. X-rays, microwaves, neutrons
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L5/00—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
    • G01L5/0047—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes measuring forces due to residual stresses
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/20—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by using diffraction of the radiation by the materials, e.g. for investigating crystal structure; by using scattering of the radiation by the materials, e.g. for investigating non-crystalline materials; by using reflection of the radiation by the materials
    • G01N23/2055—Analysing diffraction patterns
    • 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
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/60—Specific applications or type of materials
    • G01N2223/607—Specific applications or type of materials strain

Definitions

  • the present invention relates to a residual stress calculation method.
  • Patent Document 1 describes a residual stress measurement method that can analyze a residual stress for a flat joint or an axisymmetric joint.
  • the residual stress measurement method described in Patent Document 1 is a residual stress analysis method using a general intrinsic strain method.
  • the residual stress is calculated by the following procedure. First, by cutting a structure to be analyzed for residual stress, elastic release strains in three directions are measured respectively. A coefficient matrix indicating the relationship between the unit value of the intrinsic strain and the elastic strain is obtained by finite element analysis. The most probable value of the inherent strain in the three directions is obtained using the measured elastic strain in the three directions and the obtained coefficient matrix. Then, the most probable value of the residual stress is calculated using the most probable value of the inherent strain and the coefficient matrix.
  • the inherent strain method it is necessary to expose the inside of the structure as a cutting surface by cutting the structure in order to measure residual stress or elastic strain generated in the structure.
  • measurement errors always occur in the measurement of residual stress or elastic strain on the surface of the structure and the cutting surface.
  • the number of measurement points of residual stress or elastic strain in the inherent strain method is very large, and is usually several tens to several hundreds. For this reason, it is difficult to appropriately consider measurement errors.
  • the present invention has been made based on the above-described circumstances, and a residual stress calculation method that can easily calculate a residual stress without obtaining a coefficient matrix representing a relationship between intrinsic strain and elastic strain by a finite element method.
  • the purpose is to provide.
  • a first aspect of the present invention made to solve the above-described problem is a method for calculating a residual stress of a metal material using an inherent strain method, which cuts a plurality of locations on the surface of the metal material and the metal material.
  • a step of preparing measured values of residual stress at each of a plurality of locations on the cutting surface formed by the step, and a plurality of finite elements having the same shape as the metal material after cutting, including a part of the surface and the cutting surface A step of creating a model, a step of setting a parameter capable of imparting an inherent strain for each of the plurality of finite element models created in the creation step, and a plurality of the parameters set in the setting step For each of the finite element models, a step of deriving a calculated value of residual stress at each of the plurality of locations based on finite element analysis, and a residual stress derived in the derivation step A difference between the calculated value and the measured value of the residual stress prepared in the preparing step for each of the plurality of locations, and
  • This residual stress calculation method derives a calculated value of residual stress based on finite element analysis for each of a plurality of finite element models having the same shape as the metal material after cutting. Then, the residual stress calculation method finds the optimal solution for the parameter that gives the inherent strain by solving the optimization problem that minimizes the difference between the calculated value and the measured value of the residual stress, and the parameter that becomes this optimal solution Is used to calculate the residual stress in the finite element model of the same shape as the metal material before cutting. Therefore, the residual stress calculation method can easily calculate the residual stress without obtaining a coefficient matrix representing the relationship between the inherent strain and the elastic strain by the finite element method.
  • a second aspect of the present invention made to solve the above-described problem is a method for calculating a residual stress of a metal material using an inherent strain method, in which a plurality of locations on the surface of the metal material and the metal material are cut.
  • a step of preparing measured values of elastic strain at each of a plurality of locations on the cutting surface formed by a plurality of finite elements having the same shape as the metal material after cutting, including a part of the surface and the cutting surface A step of creating a model, a step of setting a parameter capable of imparting an inherent strain for each of the plurality of finite element models created in the creation step, and a plurality of the parameters set in the setting step
  • For each of the finite element models, a step of deriving a calculated value of elastic strain at each of the plurality of locations based on finite element analysis, and an elasticity derived in the derivation step Calculating the difference between the calculated value of the strain and the measured value of the elastic strain prepared in the preparing step for each of the plurality
  • This residual stress calculation method derives a calculated value of elastic strain based on a finite element analysis for each of a plurality of finite element models having the same shape as the metal material after cutting. Then, the residual stress calculation method finds the optimal solution for the parameter that gives the inherent strain by solving the optimization problem that minimizes the difference between the calculated value of the elastic strain and the measured value, and the parameter that becomes this optimal solution Is used to calculate the residual stress in the finite element model of the same shape as the metal material before cutting. Therefore, the residual stress calculation method can easily calculate the residual stress without obtaining a coefficient matrix representing the relationship between the inherent strain and the elastic strain by the finite element method.
  • the metal material includes a cylindrical shaft portion and a plate-like portion protruding from the shaft portion in the entire circumference and in the radial direction, and the surface of the metal material is formed at a connection portion between the shaft portion and the plate-like portion. It is preferable to include a fillet surface.
  • projecting from the shaft part in the entire circumference and the radial direction means projecting in the radial direction over the entire circumference from the outer peripheral surface of the shaft part.
  • the residual stress calculation method uses the range including the fillet surface as a residual stress calculation target, the residual stress generated in the metal material can be appropriately evaluated.
  • the fillet surface includes a region formed in an arc shape, and an extended surface of the cutting surface passes through a center position of the arc of the fillet surface.
  • the residual stress calculation method can efficiently perform the finite element analysis by defining the center position of the arc of the fillet surface as the center of the local coordinate system.
  • the parameters are initially set based on at least one of a half-value width of diffraction X-rays in the X-ray stress measurement for the metal material, a hardness of the metal material, a similar case, and a finite element analysis. Good.
  • the half-width of diffracted X-rays, the hardness of the metal material, the similar case, and the finite element analysis in the X-ray stress measurement for the metal material make it possible to estimate the plastic strain existing inside the metal material.
  • the convergence of the design variable to the optimal solution becomes faster as the initial value of the design variable is closer to the optimal solution. Since the residual stress calculation method can select a value close to the optimal solution as the initial value of the parameter by estimating plastic strain, the optimization problem can be solved efficiently.
  • the adjustment range of the parameter when solving the optimization problem is at least one of the half-value width of diffracted X-rays in the X-ray stress measurement for the metal material, the hardness of the metal material, similar cases, and finite element analysis. It may be set according to the constraint condition based on it. Since the residual stress calculation method employs reasonable constraint conditions for the parameters, the optimization problem can be solved efficiently.
  • the residual stress calculation method when the calculated value falls within the error range of the measured value in the calculating step, the difference between the calculated value and the measured value is set to zero, and the calculated value is an error of the measured value.
  • the difference from the calculated value to the boundary of the error range may be the difference between the calculated value and the measured value. Since the residual stress calculation method calculates the difference between the calculated value and the measured value based on the error range of the measured value, the measurement error can be appropriately taken into account even when the number of measurement points is very large.
  • the residual stress calculation method it is preferable to use a condition that the volume of the metal material does not change due to generation of the inherent strain as one of the constraint conditions when solving the optimization problem. Thereby, the residual stress calculation method can reduce design variables in the optimization problem.
  • the above parameter is preferably a coefficient of intrinsic strain or intrinsic strain approximation function. Thereby, the residual stress calculation method can directly calculate the optimum solution of the inherent strain.
  • the residual stress calculation method of the present invention can easily calculate the residual stress without obtaining a coefficient matrix representing the relationship between the inherent strain and the elastic strain by the finite element method.
  • FIG. 1 is a flowchart showing a residual stress calculation method according to the first embodiment of the present invention.
  • FIG. 2 is a schematic longitudinal sectional view of a metal material used in the residual stress calculation method of FIG.
  • FIG. 3 is a schematic perspective sectional view showing a part of the metal material of FIG.
  • FIG. 4 is an explanatory diagram showing details of a procedure for obtaining the optimum solution of FIG.
  • FIG. 5 is a graph showing the relationship between the X-ray half width and the equivalent intrinsic strain.
  • the residual stress calculation method shown in FIG. 1 is a residual stress calculation method for the metal material 1 using the inherent strain method, and is formed by cutting a plurality of locations on the surface of the metal material 1 and the metal material 1.
  • a preparation step S1 for preparing measured values of residual stress at each of a plurality of locations on the cutting surface, and a plurality of finite element models having the same shape as the metal material 1 after cutting including a part of the surface of the metal material 1 and the cutting surface.
  • a setting step S3 for setting a parameter capable of imparting an inherent strain for each of a plurality of finite element models created in the creating step S2, and a plurality of parameters set in the setting step S3.
  • a derivation step S4 for deriving a calculated value of the residual stress at each of a plurality of locations based on finite element analysis, and a residual stress derived in the derivation step S4 And a calculation step S5 for calculating for each of the difference between the measured value of the residual stress that is prepared by the calculated value and the preparing step S1 in a plurality of locations.
  • the residual stress calculation method repeatedly executes the setting step S3, the derivation step S4, and the calculation step S5 while adjusting the parameters as an optimization problem that minimizes the difference between the calculated value and the measured value calculated in the calculation step S5. After obtaining the optimal solution for the parameter, set the optimal parameter to the finite element model of the same shape as the metal material 1 before cutting, and calculate the residual stress for this finite element model based on the finite element analysis. To do.
  • the load vector P is expressed by the following equation. If the inherent strain ⁇ 0 is known, the displacement u can be obtained using equations (4) to (6). By substituting the obtained displacement u into equation (3) and further using equation (1), total strain ⁇ and residual stress ⁇ can be obtained.
  • N residual stresses determined by measurement are expressed as ⁇ m
  • N residual stresses calculated when the inherent strain is ⁇ 0 are expressed as ⁇ c .
  • the difference R between ⁇ m and ⁇ c is defined by the following equation.
  • the inherent strain can be approximately expressed by a distribution function.
  • the inherent strain ⁇ 0 at an arbitrary point is expressed by the following linear function with M distribution function parameters a.
  • M is a function of coordinates, and may be non-linear with respect to the coordinates.
  • ⁇ c can be obtained using equations (1) to (6). Therefore, ⁇ c can be expressed by the following linear function.
  • H is a coefficient matrix representing the relationship between intrinsic strain and elastic strain.
  • the residual stress calculation method a parameter capable of applying the inherent strain ⁇ 0 is set in a finite element model, and a calculated value of the residual stress is derived by finite element analysis for the finite element model in which the parameter is set. Then, in the residual stress calculation method, the most probable value of the parameter is obtained by solving an optimization problem that minimizes the difference R between the calculated value and the measured value of the residual stress using Equation (7). That is, the residual stress calculation method calculates the most probable value of the inherent strain ⁇ 0 without obtaining the coefficient matrix H indicating the relationship between the inherent strain and the elastic strain.
  • the metal material 1 used in the residual stress calculation method includes a columnar shaft portion 2 and a plate-like portion 3 protruding from the shaft portion 2 in the entire circumference and in the radial direction. Further, an annular fillet surface X is formed at a connection portion between the shaft portion 2 and the plate-like portion 3 so as to go around the shaft portion 2. A plurality of locations on the fillet surface X serve as measurement areas for residual stress. That is, the surface of the metal material 1 that is the measurement target region of the residual stress includes the fillet surface X formed at the connection portion between the shaft portion 2 and the plate-like portion 3. Examples of the material of the metal material 1 include metals such as forged steel, cast steel, and titanium. In addition, the number of the plate-shaped parts 3 is not limited to four shown in FIG.
  • the fillet surface X is formed in a substantially arc shape by cold working. Specifically, as shown in FIG. 3, the metal material 1 is formed so as to include a region in which the fillet surface X is formed in an arc shape in a longitudinal section (a divided surface passing through the center line of the shaft portion 2). Yes.
  • the fillet surface X only needs to be formed so as to include an arc-shaped region having a predetermined diameter (fillet diameter), and may partially include a region that is not arc-shaped.
  • the fillet surface X includes a region formed in an arc shape and is formed in an annular shape so as to go around the shaft portion 2, and thus constitutes a part of the annular surface as a whole. Therefore, the center position of the arc of the fillet surface X draws a circular orbit centered on the center line of the shaft portion 2. Therefore, in the longitudinal section of the metal material 1, the center position of the arc is the origin of the local coordinate system, and the moving radius from this origin is defined as r.
  • a perpendicular direction from the center position of the arc of the fillet surface X toward the center line of the shaft portion 2 is defined as a reference direction
  • a declination from the reference direction in the longitudinal section of the metal material 1 is defined as ⁇ .
  • the depth from the fillet surface X can be expressed by r
  • the position of the fillet surface X in the arc can be expressed by ⁇ .
  • the declination angle is defined as ⁇ with respect to the tangential direction of the circular orbit drawn by the center position of the arc of the fillet surface X.
  • an information processing apparatus including at least a control unit, a storage unit, a display unit, and an operation unit is used.
  • This information processing apparatus is, for example, a general-purpose computer, and operates so as to realize the residual stress calculation method when the control unit reads a program stored in the storage unit.
  • the information processing apparatus further includes an external information acquisition unit for taking in information from outside and storing it in the storage unit. Examples of the external information acquisition unit include a communication unit that transmits and receives information via a communication line or a communication line connected to another device, a portable medium reading unit that reads information from a portable medium that stores information, and the like. .
  • the preparation step S1 is a step of preparing measured values of residual stresses at a plurality of locations on the surface of the metal material 1 and at a plurality of locations on the cutting surface formed by cutting the metal material 1.
  • the information processing apparatus takes in residual stress measurement data from the external information acquisition unit and stores it in the storage unit.
  • the communication unit receives measurement data transmitted from the measurement device and stores the measurement data in the storage unit.
  • the measurement data includes at least data that associates the measurement points of the residual stress with the measurement values of the residual stress. Since the metal material 1 is a metal structure in which the fillet surface X is formed by cold working, an inherent strain is generated in the vicinity of the fillet surface X. That is, it is presumed that the elastic strain derived from the inherent strain is generated in the range from the fillet surface X to a certain depth. Therefore, the measurement location of the residual stress in the metal material 1 is selected over a range in which it is estimated that the elastic strain is generated due to the inherent strain.
  • the measurement data may include residual stress measurement direction data and shape data of the metal material 1 after cutting.
  • the metal material 1 is cut so as to have various cutting surfaces in order to measure the residual stress in various directions.
  • the metal material 1 is cut in consideration of the origin in the local coordinate system.
  • the metal material 1 is cut so that the extended surface of the cutting surface of the metal material 1 passes through the center position of the arc of the fillet surface X.
  • the creation step S2 is a step of creating a plurality of finite element models having the same shape as the metal material 1 after cutting, including a part of the surface of the metal material 1 and the cutting surface.
  • the information processing apparatus creates a finite element model having the same shape as the metal material 1 when the residual stress is measured based on the operation of the operation unit.
  • Commercially available finite element model creation software may be used for creating the finite element model.
  • the creation step S2 Since the metal material 1 is cut so as to have various cutting surfaces in the measurement of the residual stress, in the creation step S2, a part of the surface of the metal material 1 and the metal material 1 after cutting including the cutting surface have the same shape. Multiple finite element models are created. In addition, when the shape data of the metal material 1 after cutting is included in the measurement data of the residual stress, a finite element model may be created based on this shape data.
  • the finite element model is formed by a plurality of elements having unique identification information.
  • the number of divisions of the finite element model is arbitrarily set according to the balance between the calculation capability of the information processing apparatus and the calculation accuracy of the residual stress. However, from the viewpoint of obtaining sufficient calculation accuracy, it is preferable that the division size of the finite element model in the range where the elastic strain is estimated to be derived from the inherent strain is smaller than the division size in the other ranges.
  • the finite element model is preferably divided according to the radius r, the deflection angle ⁇ , and the deflection angle ⁇ in the local coordinate system. When the finite element model is divided in this way, finite element analysis becomes easy.
  • the setting step S3 is a step of setting a parameter that can give an inherent strain to each of the plurality of finite element models created in the creation step S2.
  • the information processing apparatus sets predetermined parameters for each element of the finite element model for each of the plurality of finite element models.
  • the parameters set for each element are parameters that can give inherent strain to the finite element model.
  • parameters that can be given an inherent strain include intrinsic strain, a coefficient of an intrinsic strain approximation function, temperature, energy, and the like. It is preferable that this parameter is a coefficient of an inherent strain or an approximate function of an inherent strain since an optimal solution of an inherent strain can be directly calculated in an analysis of an optimization problem.
  • the coefficient of the inherent strain approximation function is selected as a parameter, the unknown strain is aggregated in the distribution function parameter by approximating the inherent strain with the distribution function. Residual stress can be derived with high accuracy.
  • the half-value width of the diffracted X-rays in the X-ray stress measurement for the metal material 1, the hardness of the metal material 1, similar cases, and the finite element analysis make it possible to estimate the plastic strain existing inside the metal material 1.
  • the parameter which can give an inherent strain is initially set based on at least one of the half-value width of the diffraction X-ray in the X-ray stress measurement for the metal material 1, the hardness of the metal material 1, a similar case, and the finite element analysis. It is preferred if When parameters are initially set by plastic strain estimation, the efficiency of solving optimization problems is improved.
  • the parameters initially set in each element may all have the same value, or some or all may have different values.
  • the parameter may be initially set to an arbitrary value or may not be initially set by estimation of plastic strain.
  • the derivation step S4 is a step of deriving calculated values of residual stresses at a plurality of locations based on the finite element analysis for each of the plurality of finite element models for which parameters are set in the setting step S3.
  • the information processing apparatus derives, for each of the plurality of finite element models, the calculated value of the residual stress at each of a plurality of locations corresponding to the measurement location of the residual stress in the metal material 1 based on the finite element analysis.
  • a well-known method can be used for this finite element analysis. For this reason, the information processing apparatus can perform finite element analysis based on, for example, commercially available finite element analysis software.
  • the calculation step S5 is a step of calculating the difference between the calculated value of the residual stress derived in the derivation step S4 and the measured value of the residual stress prepared in the preparation step S1 for each of a plurality of locations.
  • the information processing apparatus calculates the residual stress calculated using the finite element model and the residual stress measurement stored in the storage unit for each of the measurement points of the residual stress in the metal material 1. Calculate the difference.
  • the difference between the calculated value and the measured value is calculated based on the error range of the measured value. Specifically, when the information processing apparatus calculates the difference between the calculated value and the measured value, and the calculated residual stress value is within the error range of the measured residual stress value, the difference between the calculated value and the measured value is zero. If the calculated value falls outside the error range of the measured value, the calculation is executed so that the difference between the calculated value and the boundary of the error range of the measured value is the difference between the calculated value and the measured value. For this reason, even when the number of residual stress measurement points is very large, a calculation result in which measurement errors are appropriately taken into consideration can be obtained. In addition, as an error range of a measured value, it is good to use the standard deviation for every measurement location, for example.
  • the information processing apparatus sets the objective function to minimize the difference between the calculated value of the residual stress and the measured value of the residual stress, and adjusts the parameter by setting the parameter as a setting variable. While calculating the difference repeatedly based on the finite element model, an optimal parameter solution that minimizes the difference is obtained.
  • a well-known method such as a Hooke-Jeves method which is a discrete search method is used. For this reason, the information processing apparatus can execute an analysis of an optimization problem based on commercially available optimization software.
  • the objective function for example, the sum of squares of the difference between the calculated value of the residual stress and the measured value of the residual stress may be used, or the maximum value of the difference between the calculated value of the residual stress and the measured value of the residual stress. May be used.
  • the determination of whether or not the optimal solution for the parameter has been obtained is performed by determining whether or not a certain period of time has elapsed without significant change in the above-mentioned difference in the optimization problem analysis. In this case, it may be determined that an optimal parameter solution has been obtained.
  • the parameter adjustment range for solving the optimization problem is a constraint based on at least one of the half-value width of the diffracted X-ray in the X-ray stress measurement for the metal material 1, the hardness of the metal material 1, similar cases, and finite element analysis. It may be set according to conditions. For example, when there are similar study cases and the variation of the inherent strain distribution can be grasped, the parameter adjustment range may be set according to the constraint condition based on the variation range.
  • the measurement result of the hardness of the metal material 1, the half value of the diffracted X-ray in the X-ray stress measurement for the metal material may be predicted based on any of the measurement results of the width, and the parameter adjustment range may be set according to the constraint condition according to the prediction result.
  • the parameter adjustment range may be set with at least one of a maximum value and a minimum value depending on a constraint condition.
  • the residual stress calculation method may use a condition that the volume of the metal material 1 does not change due to the generation of the inherent strain, as one of the constraint conditions when solving the optimization problem.
  • the residual stress of the metal material 1 is generated by plastic strain, it can be considered that the volume of the metal material 1 is not changed although the shape of the plastically deformed portion changes. Therefore, the residual stress calculation method can reduce design variables in the optimization problem by using this constraint condition.
  • N finite element models are used as shown in FIG. 4 (N is a natural number of 2 or more).
  • an inherent strain value is set as a parameter for each of the N finite element models.
  • a derivation step S4 for each of the N finite element models, a calculated value of the residual stress at each of a plurality of locations corresponding to the measurement location of the residual stress in the metal material 1 is derived based on the finite element analysis. Since the residual stress is calculated for each of a plurality of locations, the calculated value of the residual stress is associated as a residual stress group for each finite element model.
  • calculation step S5 the difference between the calculated value and the measured value of the residual stress is calculated for each of the plurality of locations.
  • the calculated difference is associated as a difference group for each finite element model.
  • the residual stress calculation method After obtaining an optimum parameter solution by solving an optimization problem, the parameter that becomes the optimum solution is set to a finite element model having the same shape as that of the metal material 1 before cutting. Residual stress is calculated for the model based on finite element analysis. Specifically, after the information processing apparatus obtains an optimum parameter solution by solving the above optimization problem, the parameter that becomes the optimum solution is set to a finite element model having the same shape as the metal material 1 before cutting. To do. The information processing apparatus calculates the residual stress at an arbitrary position based on the finite element analysis for this finite element model. Thereby, the said residual stress calculation method can calculate the estimated value of the residual stress in the arbitrary positions of the metal material 1.
  • the residual stress calculation method derives a calculated value of the residual stress based on the finite element analysis for each of a plurality of finite element models having the same shape as the metal material 1 after cutting. Then, the residual stress calculation method finds the optimal solution for the parameter that gives the inherent strain by solving the optimization problem that minimizes the difference between the calculated value and the measured value of the residual stress, and the parameter that becomes this optimal solution Is used to calculate the residual stress in the finite element model having the same shape as that of the metal material 1 before cutting. Therefore, the residual stress calculation method can easily calculate the residual stress without obtaining a coefficient matrix representing the relationship between the inherent strain and the elastic strain by the finite element method.
  • the residual stress calculation method uses the range including the fillet surface X of the metal material 1 as a target for calculating the residual stress, the residual stress generated in the metal material 1 can be appropriately evaluated.
  • the fillet surface X includes a region formed in an arc shape and the extended surface of the cutting surface of the metal material 1 passes through the center position of the arc of the fillet surface X, the residual stress calculation method is the arc of the fillet surface X.
  • the calculation step S5 if the calculated value of the residual stress is within the error range of the measured value, the difference between the calculated value and the measured value is zero, and the calculated value is outside the error range of the measured value. In this case, since the difference from the calculated value to the boundary of the error range is the difference between the calculated value and the measured value, the measurement error can be appropriately considered even when the number of measurement points is very large.
  • the calculated value of the residual stress is derived based on the finite element analysis for each of the plurality of finite element models having the same shape as the metal material 1 after cutting, and the difference between the calculated value of the residual stress and the measured value is calculated.
  • the finite element analysis is performed for each of the plurality of finite element models having the same shape as the metal material 1 after cutting. Even if the calculated value of elastic strain is derived based on the above, and the optimization problem that minimizes the difference between the calculated value of elastic strain and the measured value is solved, Good.
  • the residual stress calculation method is a residual stress calculation method of the metal material 1 using the inherent strain method, and a plurality of locations on the surface of the metal material 1 and a cutting surface formed by cutting the metal material 1.
  • a setup process that sets parameters that can give inherent strain to each of the multiple finite element models created in the creation process, and a plurality of finite element models for which parameters are set in the setup process
  • the derivation process for deriving the calculated value of elastic strain at each of multiple locations based on finite element analysis, the calculated value of elastic strain derived in the derivation process, and the preparation process A calculation process for calculating the difference between the measured values of the elastic strain for each of a plurality of locations, and adjusting parameters as an optimization problem to minimize the difference between the calculated value and the measurement value calculated in the calculation process.
  • the parameter serving as the optimal solution is set to a finite element model having the same shape as that of the metal material 1 before cutting.
  • the residual stress may be calculated for the model based on finite element analysis.
  • the metal material 1 is provided with the column-shaped axial part 2 and the four plate-shaped parts 3 which protrude from the axial part 2 to a perimeter and radial direction, and the axial part 2 and the plate-shaped part 3 are the following.
  • the metal material used for the residual-stress calculation method is not limited to the thing of such a structure.
  • the metal material may not have a fillet surface.
  • the extension surface of the cutting surface of the metal material 1 is an arc of the fillet surface X
  • the extended surface of the cutting surface of the metal material 1 may not pass through the center position of the arc of the fillet surface X.
  • the difference between the calculated value and the measured value is set to zero, and the calculated value is outside the error range of the measured value.
  • the calculation is performed so that the difference from the calculated value to the boundary of the error range of the measured value is the difference between the calculated value and the measured value.
  • the difference between the measured values may be calculated. However, it is preferable to calculate the difference between the calculated value and the measured value based on the error range of the measured value from the viewpoint of obtaining a calculation result that appropriately takes the measurement error into consideration.
  • the preparation step S1 is executed before the creation step S2, but the preparation step S1 may be executed at an arbitrary timing before the calculation step S5.
  • Cutting pieces of three types of shapes, F piece, C piece, and T piece were prepared from a metal material so that all the residual stresses in the three directions could be measured.
  • Each of the F piece, the C piece, and the T piece is formed such that the extended surface of the cutting surface passes through the origin of the local coordinate system.
  • the T piece has a shape of a thin piece obtained by slicing a metal material by slightly changing ⁇ so as to pass through the center line of the shaft portion. Residual stress in the r direction and ⁇ direction inside the fillet surface while changing ⁇ . Can be measured.
  • the residual stress at a depth of 0.3 mm and 0.5 mm from the fillet surface was measured using an X-ray stress measuring device (MSF-3M manufactured by Rigaku Corporation).
  • the sin 2 ⁇ method was used to measure the residual stress. Measuring ⁇ direction stress while changing ⁇ every 10 degrees in the range of 80 degrees to 100 degrees, measuring ⁇ direction stress changing ⁇ every 10 degrees from 0 degrees to 100 degrees, totaling 28 points The measured value was obtained.
  • the residual stress at a depth of 18 mm to 40 mm from the fillet surface was measured using an X-ray stress measuring device (MSF-3M manufactured by Rigaku Corporation).
  • the sin 2 ⁇ method was used to measure the residual stress.
  • the r-direction stress and the ⁇ -direction stress are measured 10 to 16 points at a time while changing ⁇ by 20 degrees in the range of 0 to 100 degrees by cutting, and for the Cb piece, ⁇ is 10 degrees by cutting.
  • the stress in the r direction and the stress in the ⁇ direction were measured at 10 to 16 points while changing by 20 degrees in the range of 110 to 110 degrees, and a total of 308 measurement values were obtained.
  • the residual stress at a depth of 18 mm to 40 mm from the fillet surface was measured using an X-ray residual stress measuring device (Pulstec ⁇ -360).
  • the cos ⁇ method was used to measure the residual stress.
  • the r-direction stress and the ⁇ -direction stress were measured at 10 to 16 points while changing ⁇ every 10 degrees in the range of 0 to 100 degrees, and a total of 222 measurement values were obtained.
  • finite element models having the same shape as the F piece, C piece, and T piece were created using finite element method analysis software.
  • C piece since ⁇ was changed by cutting, 12 types of finite element models having the same shape as the shape when the residual stress was measured were created.
  • the inherent strain was set by the following distribution function for the created finite element model.
  • R indicates the r direction
  • ⁇ indicates the ⁇ direction.
  • R 0 and ⁇ 0 indicate the starting position of the intrinsic strain calculation
  • ⁇ R and ⁇ indicate incremental values.
  • M and n represent the orders in the r direction and the ⁇ direction, respectively
  • a ij represents a coefficient of the approximate function.
  • R 0 , ⁇ 0 , ⁇ R, and ⁇ are parameters that determine the intrinsic strain calculation region. If these set values are different, a difference occurs in the obtained inherent strain distribution even if the same residual stress data is used. Similarly, the obtained inherent strain distribution varies depending on the values of the orders m and n of the inherent strain distribution function.
  • the coefficients A ij used for calculating the distortions of these components are a ij , b ij , c ij , and d ij , respectively.
  • the equivalent plastic strain has a correlation with the half width (X-ray half width) in the X-ray residual stress measurement. Further, if it is considered that the inherent strain is generated only by the plastic strain, the inherent strain has a correlation with the X-ray half width. Therefore, the initial value of the coefficient of the approximate function of the inherent strain was determined based on the measurement result of the X-ray half width. If m is 2 and n is 8, each of a ij , b ij , c ij and d ij is 16 coefficients. Table 1 shows the initial values of the coefficients of the approximate function determined based on the measurement result of the X-ray half width under these conditions.
  • the residual stress at the position corresponding to the position where the residual stress was actually measured was calculated by finite element analysis using a finite element model in which the inherent strain was set by the distribution function.
  • the square of the difference between the measured value and the calculated value of the residual stress was calculated, and these were all added up.
  • repeat these calculations coefficients a ij, b ij, c ij while adjusting the values of d ij by Hooke-Jeeves method such that the sum of the squares of the differences is minimized.
  • the sum of the squares of the differences did not change significantly in the calculation for a certain period of time, it was determined that the sum of the squares of the differences was minimized.
  • a hole is provided so that the measurement site is exposed using a drill, a strain gauge is attached to the bottom of the hole, and the metal material is cut into the smallest pieces to include the portion with the strain gauge attached.
  • the residual stress is grasped by measuring the elastic release strain before and after cutting.
  • the difference between the measured value and the calculated value of residual stress if the calculated value is within the error range of the measured value, the difference between the calculated value and the measured value is zero, and the calculated value is within the error range of the measured value. If it is outside, the calculation considering the error range in which the difference between the calculated value and the boundary of the error range is the difference between the calculated value and the measured value was also examined. As the error range of the measurement value, the standard deviation for each measurement location was used.
  • FIG. 6 shows the results (Me) of the residual stress measured by (C3) and the cutting method, respectively.
  • the residual stress calculation method of the present invention can easily calculate the residual stress without obtaining a coefficient matrix representing the relationship between the inherent strain and the elastic strain by the finite element method.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Analytical Chemistry (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Toxicology (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Heat Treatment Of Articles (AREA)
  • Control Of Heat Treatment Processes (AREA)

Abstract

本発明の一態様に係る残留応力算出方法は、金属材の複数箇所のそれぞれにおける残留応力の測定値を用意する工程と、金属材と同一形状の複数の有限要素モデルを作成する工程と、複数の有限要素モデルのそれぞれに対し、固有ひずみを付与可能なパラメータを設定する工程と、パラメータを設定された複数の有限要素モデルのそれぞれについて、有限要素解析に基づき複数箇所のそれぞれにおける残留応力の計算値を導出する工程と、残留応力の計算値及び測定値の差を算出する工程とを備え、計算値と測定値との差を最小化する最適化問題として、パラメータの最適解を求めた後、切削前の金属材と同一形状の有限要素モデルについて有限要素解析に基づき残留応力を計算する。

Description

残留応力算出方法
　本発明は、残留応力算出方法に関する。
　金属製構造物の内部には、構造物の製造条件に由来して熱ひずみ、塑性ひずみ等の固有ひずみが生成される。生成された固有ひずみは弾性ひずみを誘起する。このため、構造物には、固有ひずみに由来する残留応力が生じる。
　構造物の内部に生成された固有ひずみは、弾性ひずみを解放するように構造物が分断されても変化しない。そこで、構造物を適切に分断することにより固有ひずみに由来する残留応力又は弾性解放ひずみを測定し、この測定値を用いた逆解析によって構造物内に存在する固有ひずみを適切に把握し、把握した固有ひずみから構造物内の残留応力を算出する固有ひずみ法が知られている。例えば特許文献１には、平板継手又は軸対称継手について残留応力を解析可能な残留応力測定方法が説明されている。
　特許文献１に記載された残留応力測定方法は、一般的な固有ひずみ法を用いた残留応力の解析法である。この方法では、次の手順により残留応力が算出される。まず、残留応力の解析対象となる構造物を切断することにより、３つの方向の弾性解放ひずみがそれぞれ測定される。固有ひずみの単位値と弾性ひずみとの関係を示す係数マトリックスが有限要素解析により求められる。測定された３方向の弾性ひずみと求められた係数マトリックスとを用いて３方向の固有ひずみの最確値が求められる。そして、固有ひずみの最確値と係数マトリックスとを用いて残留応力の最確値が算出される。
　特許文献１の残留応力測定方法では、固有ひずみと弾性ひずみとの関係を示す係数マトリックスを求めることになるが、この計算には膨大な時間と手間が必要である。そこで、計算の手間を省くために、有限要素解析により係数マトリックスを自動的に計算する専用プログラムを製作することが考えられる。しかしながら、この専用プログラムの製作には、有限要素法に関する専門知識が必要となる。
　また、固有ひずみ法では、構造物の内部に生じている残留応力又は弾性ひずみの測定のために、構造物を切削することによって構造物の内部を切削面として露出させる必要がある。しかしながら、構造物の表面及び切削面における残留応力又は弾性ひずみの測定では必ず測定誤差が生じる。また、固有ひずみ法における残留応力又は弾性ひずみの測定点の数は非常に多く、通常、数十から数百である。このため、適切に測定誤差を考慮することは困難である。
特開２００５－１８１１７２号公報
　本発明は、上述のような事情に基づいてなされたものであり、固有ひずみと弾性ひずみとの関係を表す係数マトリックスを有限要素法で求めることなく、簡単に残留応力を算出できる残留応力算出方法を提供することを目的とする。
　上記課題を解決するためになされた本発明の第１の態様は、固有ひずみ法を用いた金属材の残留応力算出方法であって、金属材の表面上の複数箇所及びこの金属材を切削することにより形成される切削面上の複数箇所のそれぞれにおける残留応力の測定値を用意する工程と、上記表面の一部及び上記切削面を含む切削後の上記金属材と同一形状の複数の有限要素モデルを作成する工程と、上記作成工程で作成された複数の上記有限要素モデルのそれぞれに対し、固有ひずみを付与可能なパラメータを設定する工程と、上記設定工程で上記パラメータを設定された複数の上記有限要素モデルのそれぞれについて、有限要素解析に基づき上記複数箇所のそれぞれにおける残留応力の計算値を導出する工程と、上記導出工程で導出された残留応力の上記計算値及び上記用意工程で用意された残留応力の上記測定値の差を上記複数箇所のそれぞれについて算出する工程とを備え、上記算出工程で算出される上記計算値と上記測定値との差を最小化する最適化問題として、上記パラメータを調整しつつ上記設定工程、上記導出工程及び上記算出工程を繰り返し実行することにより上記パラメータの最適解を求めた後、この最適解となるパラメータを切削前の上記金属材と同一形状の有限要素モデルに設定し、この有限要素モデルについて有限要素解析に基づき残留応力を計算する。
　当該残留応力算出方法は、切削後の金属材と同一形状の複数の有限要素モデルのそれぞれについて、有限要素解析に基づき残留応力の計算値を導出する。そして、当該残留応力算出方法は、残留応力の計算値と測定値との差を最小化する最適化問題を解くことにより、固有ひずみを付与するパラメータの最適解を求め、この最適解となるパラメータを用いて切削前の金属材と同一形状の有限要素モデルにおける残留応力を計算する。したがって、当該残留応力算出方法は、固有ひずみと弾性ひずみとの関係を表す係数マトリックスを有限要素法で求めることなく、簡単に残留応力を算出できる。
　上記課題を解決するためになされた本発明の第２の態様は、固有ひずみ法を用いた金属材の残留応力算出方法であって、金属材の表面上の複数箇所及びこの金属材を切削することにより形成される切削面上の複数箇所のそれぞれにおける弾性ひずみの測定値を用意する工程と、上記表面の一部及び上記切削面を含む切削後の上記金属材と同一形状の複数の有限要素モデルを作成する工程と、上記作成工程で作成された複数の上記有限要素モデルのそれぞれに対し、固有ひずみを付与可能なパラメータを設定する工程と、上記設定工程で上記パラメータを設定された複数の上記有限要素モデルのそれぞれについて、有限要素解析に基づき上記複数箇所のそれぞれにおける弾性ひずみの計算値を導出する工程と、上記導出工程で導出された弾性ひずみの上記計算値及び上記用意工程で用意された弾性ひずみの上記測定値の差を上記複数箇所のそれぞれについて算出する工程とを備え、上記算出工程で算出される上記計算値と上記測定値との差を最小化する最適化問題として、上記パラメータを調整しつつ上記設定工程、上記導出工程及び上記算出工程を繰り返し実行することにより上記パラメータの最適解を求めた後、この最適解となるパラメータを切削前の上記金属材と同一形状の有限要素モデルに設定し、この有限要素モデルについて有限要素解析に基づき残留応力を計算する。
　当該残留応力算出方法は、切削後の金属材と同一形状の複数の有限要素モデルのそれぞれについて、有限要素解析に基づき弾性ひずみの計算値を導出する。そして、当該残留応力算出方法は、弾性ひずみの計算値と測定値との差を最小化する最適化問題を解くことにより、固有ひずみを付与するパラメータの最適解を求め、この最適解となるパラメータを用いて切削前の金属材と同一形状の有限要素モデルにおける残留応力を計算する。したがって、当該残留応力算出方法は、固有ひずみと弾性ひずみとの関係を表す係数マトリックスを有限要素法で求めることなく、簡単に残留応力を算出できる。
　上記金属材が円柱状の軸部とこの軸部から全周かつ径方向に突出する板状部とを備え、上記金属材の上記表面が上記軸部と上記板状部との接続部分に形成されるフィレット面を含むとよい。ここで、上述の「軸部から全周かつ径方向に突出する」とは、軸部の外周面から全周に亘って径方向に突出する、ということである。軸部と板状部との接続部分に形成されるフィレット面には、製造時に固有ひずみが生成されやすい。当該残留応力算出方法は、このフィレット面を含む範囲を残留応力の算出対象とするので、金属材の内部に生じる残留応力を適切に評価できる。
　上記フィレット面が円弧状に形成された領域を含み、上記切削面の延長面が上記フィレット面の円弧の中心位置を通るとよい。これにより、当該残留応力算出方法は、フィレット面の円弧の中心位置を局所座標系の中心として定義することで、有限要素解析を効率よく実施できる。
　上記設定工程で、上記パラメータが、上記金属材に対するＸ線応力測定における回折Ｘ線の半価幅、上記金属材の硬さ、類似事例及び有限要素解析の少なくともいずれかに基づき初期設定されるとよい。金属材に対するＸ線応力測定における回折Ｘ線の半価幅、金属材の硬さ、類似事例及び有限要素解析は、金属材の内部に存在する塑性ひずみを推測可能とする。また、最適化問題では、設計変数の初期値が最適解に近いほど設計変数の最適解への収束が速くなる。当該残留応力算出方法は、塑性ひずみの推測によりパラメータの初期値として最適解に近い値を選択できるので、最適化問題を効率よく解くことができる。
　上記最適化問題を解く際の上記パラメータの調整範囲が、上記金属材に対するＸ線応力測定における回折Ｘ線の半価幅、上記金属材の硬さ、類似事例及び有限要素解析の少なくともいずれかに基づく制約条件により設定されるとよい。当該残留応力算出方法は、パラメータに合理的な制約条件を採用するので、最適化問題を効率よく解くことができる。
　当該残留応力算出方法は、上記算出工程で、上記計算値が上記測定値の誤差範囲内となる場合は、上記計算値及び上記測定値の差をゼロとし、上記計算値が上記測定値の誤差範囲外となる場合は、上記計算値から上記誤差範囲の境界までの差を上記計算値及び上記測定値の差とするとよい。当該残留応力算出方法は、測定値の誤差範囲に基づき計算値及び測定値の差を算出するので、測定点の数が非常に多い場合であっても、適切に測定誤差を考慮できる。
　当該残留応力算出方法は、上記最適化問題を解く際の制約条件の一つとして、金属材の体積が固有ひずみの生成によって変化しないという条件を用いるとよい。これにより、当該残留応力算出方法は、最適化問題における設計変数を減らすことができる。
　上記パラメータが、固有ひずみ又は固有ひずみ近似関数の係数であるとよい。これにより、当該残留応力算出方法は、固有ひずみの最適解を直接的に算出できる。
　本発明の残留応力算出方法は、固有ひずみと弾性ひずみとの関係を表す係数マトリックスを有限要素法で求めることなく、簡単に残留応力を算出できる。
図１は、本発明の第一実施形態の残留応力算出方法を示すフローチャートである。 図２は、図１の残留応力算出方法に用いられる金属材の模式的縦断面図である。 図３は、図３の金属材の一部を示す模式的斜視断面図である。 図４は、図１の最適解を求める手順の詳細を示す説明図である。 図５は、Ｘ線半価幅と相当固有ひずみとの関係を示すグラフである。 図６は、θ＝３０度におけるフィレット面からの深さと残留応力との関係を示すグラフである。
　以下、本発明の残留応力算出方法の実施形態について、図を参照しつつ詳説する。
　図１に示す残留応力算出方法は、固有ひずみ法を用いた金属材１の残留応力算出方法であって、金属材１の表面上の複数箇所及びこの金属材１を切削することにより形成される切削面上の複数箇所のそれぞれにおける残留応力の測定値を用意する用意工程Ｓ１と、金属材１の表面の一部及び切削面を含む切削後の金属材１と同一形状の複数の有限要素モデルを作成する作成工程Ｓ２と、作成工程Ｓ２で作成された複数の有限要素モデルのそれぞれに対し、固有ひずみを付与可能なパラメータを設定する設定工程Ｓ３と、設定工程Ｓ３でパラメータを設定された複数の有限要素モデルのそれぞれについて、有限要素解析に基づき複数箇所のそれぞれにおける残留応力の計算値を導出する導出工程Ｓ４と、導出工程Ｓ４で導出された残留応力の計算値及び用意工程Ｓ１で用意された残留応力の測定値の差を複数箇所のそれぞれについて算出する算出工程Ｓ５とを備えている。当該残留応力算出方法は、算出工程Ｓ５で算出される計算値と測定値との差を最小化する最適化問題として、パラメータを調整しつつ設定工程Ｓ３、導出工程Ｓ４及び算出工程Ｓ５を繰り返し実行することによりパラメータの最適解を求めた後、この最適解となるパラメータを切削前の金属材１と同一形状の有限要素モデルに設定し、この有限要素モデルについて有限要素解析に基づき残留応力を計算する。
　当該残留応力算出方法を具体的に説明する前に、残留応力算出の原理、固有ひずみ算出の原理、当該残留応力算出方法に用いられる金属材及び情報処理装置について説明する。
〔残留応力算出の原理〕
　まず、残留応力算出の原理について説明する。固有ひずみをε０とすると、残留応力σは次式で表される。
Figure JPOXMLDOC01-appb-M000001
 　ここで、Ｄは弾性係数マトリックスであり、εは次式の関係を満たす全ひずみである。
Figure JPOXMLDOC01-appb-M000002
 　体積積分の範囲は解析対象となる範囲全体である。また、Ｂは変位ｕとεを関係付ける係数マトリックスであり、それらの関係は次式で表される。
Figure JPOXMLDOC01-appb-M000003
 　一方、荷重と変位との関係は次式で表される。
Figure JPOXMLDOC01-appb-M000004
 　ここで、Ｋは剛性マトリックス、Ｐは荷重ベクトルを表す。有限要素法の基礎理論である仮想仕事の原理より、剛性マトリックスＫは次式で表される。
Figure JPOXMLDOC01-appb-M000005
 　荷重は固有ひずみのみによって生じると仮定すれば荷重ベクトルＰは次式で表される。
Figure JPOXMLDOC01-appb-M000006
 　固有ひずみε０が既知であれば、式（４）から式（６）を用いて変位ｕを求めることができる。求めた変位ｕを式（３）に代入し、さらに式（１）を用いることで、全ひずみε及び残留応力σが得られる。
〔固有ひずみ算出の原理〕
　測定によって決定されたＮ個の残留応力をσｍと表し、固有ひずみをε０とした場合に算出されるＮ個の残留応力をσｃと表す。σｍとσｃとの差Ｒを次式で定義する。
Figure JPOXMLDOC01-appb-M000007
　固有ひずみは、分布関数によって近似的に表すことができる。任意点の固有ひずみε０は、Ｍ個の分布関数パラメータａによって次式の線形関数で表される。
Figure JPOXMLDOC01-appb-M000008
 　ここで、Ｍは座標の関数であり、座標に関して非線形であっても良い。式（８）によって固有ひずみε０が表されるならば、式（１）から式（６）を用いてσｃを求めることができる。したがって、σｃは次式の線形関数で表すことができる。
Figure JPOXMLDOC01-appb-M000009
 　ここで、Ｈは固有ひずみと弾性ひずみとの関係を表す係数マトリックスである。
　一般的な固有ひずみ法では、固有ひずみε０の成分ごとに、Ｍ個の分布関数パラメータａの一つの成分に単位値を与え、その他の成分をゼロとして残留応力を求める手順を繰り返すことによって、係数マトリックスＭ及びＨを求める。そして、式（９）を式（７）に代入し、Ｒが最小となるようにＭ個の分布関数パラメータａを決定することで、固有ひずみε０の最確値を算出する。しかしながら、固有ひずみと弾性ひずみとの関係を示す係数マトリックスＨを求める計算には膨大な時間と手間が必要である。
　一方、当該残留応力算出方法は、固有ひずみε０を付与可能なパラメータを有限要素モデルに設定し、パラメータが設定された有限要素モデルについて、有限要素解析により残留応力の計算値を導出する。そして、当該残留応力算出方法は、式（７）を用いて残留応力の計算値及び測定値の差Ｒを最小化する最適化問題を解くことにより、パラメータの最確値を求める。つまり、当該残留応力算出方法は、固有ひずみと弾性ひずみとの関係を示す係数マトリックスＨを求めることなく、固有ひずみε０の最確値を算出する。
〔金属材〕
　図２に示すように、当該残留応力算出方法に用いられる金属材１は、円柱状の軸部２とこの軸部２から全周かつ径方向に突出する板状部３とを備えている。また、軸部２と板状部３との接続部分には軸部２を一周するように環状のフィレット面Ｘが形成されている。このフィレット面Ｘ上の複数箇所は、残留応力の測定対象領域となっている。つまり、残留応力の測定対象領域となる金属材１の表面が軸部２と板状部３との接続部分に形成されるフィレット面Ｘを含んでいる。金属材１の材料としては、例えば鍛鋼、鋳鋼、チタン等の金属が挙げられる。なお、板状部３の数は、図２に示す４つに限定されない。
（フィレット面）
　フィレット面Ｘは、冷間加工により略円弧状に形成されている。具体的には、図３に示すように、金属材１の縦断面（軸部２の中心線を通る分割面）において、フィレット面Ｘが円弧状に形成された領域を含むように形成されている。なお、フィレット面Ｘは、所定の直径（フィレット径）の円弧状の領域を含むように形成されればよく、円弧状でない領域を一部に含んでいてもよい。
（局所座標系）
　ここで、後述する有限要素解析を効率よく実施するために局所座標系を定義しておく。フィレット面Ｘは、円弧状に形成された領域を含み、かつ軸部２を一周するように環状に形成されているので、全体として円環面の一部を構成する。したがって、フィレット面Ｘの円弧の中心位置は、軸部２の中心線を中心とした円軌道を描く。そこで、金属材１の縦断面内において、円弧の中心位置を局所座標系の原点とし、この原点からの動径をｒと定義する。また、フィレット面Ｘの円弧の中心位置から軸部２の中心線へ向かう垂線の方向を基準方向とし、金属材１の縦断面内における基準方向からの偏角をθと定義する。これにより、フィレット面Ｘからの深さをｒで表現し、かつフィレット面Ｘの円弧における位置をθで表現することができる。さらに、フィレット面Ｘの円弧の中心位置が描く円軌道の接線方向について、偏角をφと定義する。このような局所座標系が採用されると、軸部２の周方向には均一に残留応力が生じていると仮定することにより偏角φを考慮する必要がなくなるため、フィレット面Ｘ近傍の有限要素解析が容易となる。
〔情報処理装置〕
　当該残留応力算出方法には、制御部、記憶部、表示部及び操作部を少なくとも備える情報処理装置が用いられる。この情報処理装置は、例えば汎用のコンピュータであり、記憶部に記憶されたプログラムを制御部が読み込むことによって当該残留応力算出方法を実現するように動作する。また、この情報処理装置は、外部から情報を取り込んで記憶部に記憶するための外部情報取得部をさらに備えている。外部情報取得部としては、例えば通信回線又は他の機器に接続される通信線を介して情報を送受信する通信部、情報を記憶した可搬媒体から情報を読み取る可搬媒体読取部等が挙げられる。
　以下、当該残留応力算出方法を具体的に説明する。
＜用意工程Ｓ１＞
　用意工程Ｓ１は、金属材１の表面上の複数箇所及びこの金属材１を切削することにより形成される切削面上の複数箇所のそれぞれにおける残留応力の測定値を用意する工程である。用意工程Ｓ１では、情報処理装置が、外部情報取得部から残留応力の測定データを取り込んで記憶部に記憶する。なお、情報処理装置が、残留応力の測定用装置を兼ねる構成である場合は、測定機器から送信された測定データを通信部で受信し、この測定データを記憶部に記憶する。
（測定データ）
　測定データは、残留応力の測定箇所及び残留応力の測定値を対応付けたデータを少なくとも含んで構成される。金属材１は、冷間加工によりフィレット面Ｘを形成した金属製構造物であるので、フィレット面Ｘの近傍に固有ひずみが生成される。つまり、フィレット面Ｘから一定の深さまでの範囲に固有ひずみに由来する弾性ひずみが生じていると推測される。したがって、金属材１における残留応力の測定箇所は、固有ひずみに由来して弾性ひずみが生じていると推測される範囲にわたって選択される。なお、測定データには、残留応力の測定方向のデータや切削後の金属材１の形状データが含まれていてもよい。
　また、金属材１は、残留応力の測定において、様々な方向の残留応力を測定するために様々な切削面を有するように切削されるが、特に、局所座標系における原点を考慮して切削される。具体的には、金属材１の切削面の延長面がフィレット面Ｘの円弧の中心位置を通るように金属材１は切削される。これにより、当該残留応力算出方法は、金属材１における残留応力の測定箇所を動径ｒ及び偏角θで表現できるため、有限要素解析を効率よく実施できる。
＜作成工程Ｓ２＞
　作成工程Ｓ２は、金属材１の表面の一部及び切削面を含む切削後の金属材１と同一形状の複数の有限要素モデルを作成する工程である。作成工程Ｓ２では、情報処理装置が、操作部の操作に基づいて、残留応力を測定した際の金属材１と同一形状の有限要素モデルを作成する。有限要素モデルの作成には、市販の有限要素モデル作成ソフトが用いられるとよい。
　金属材１は残留応力の測定において様々な切削面を有するように切削されるため、作成工程Ｓ２では、金属材１の表面の一部及び切削面を含む切削後の金属材１と同一形状の複数の有限要素モデルが作成される。なお、残留応力の測定データに切削後の金属材１の形状データが含まれている場合には、この形状データに基づき有限要素モデルが作成されるとよい。
（有限要素モデル）
　有限要素モデルは、固有の識別情報を有する複数の要素により形成される。有限要素モデルの分割数は、情報処理装置の計算能力と残留応力の計算精度とのバランスに応じて任意に設定される。ただし、十分な計算精度を得る観点から、固有ひずみに由来して弾性ひずみが生じていると推測される範囲における有限要素モデルの分割サイズは、他の範囲における分割サイズより小さいと好ましい。また、有限要素モデルは、局所座標系における動径ｒ、偏角θ及び偏角φに応じて分割されていると好ましい。このように有限要素モデルが分割されていると、有限要素解析が容易となる。
＜設定工程Ｓ３＞
　設定工程Ｓ３は、作成工程Ｓ２で作成された複数の有限要素モデルのそれぞれに対し、固有ひずみを付与可能なパラメータを設定する工程である。設定工程Ｓ３では、情報処理装置が、複数の有限要素モデルのそれぞれに対し、予め定められたパラメータを有限要素モデルの要素毎に設定する。
（パラメータ）
　各要素に設定されるパラメータは、有限要素モデルに対して固有ひずみを付与可能なパラメータである。固有ひずみを付与可能なパラメータとしては、固有ひずみ、固有ひずみ近似関数の係数、温度、エネルギー等が挙げられる。このパラメータは、固有ひずみ又は固有ひずみ近似関数の係数であると、最適化問題の解析において固有ひずみの最適解を直接的に算出できるので好ましい。特に、パラメータとして固有ひずみ近似関数の係数が選択されると、固有ひずみを分布関数で近似することにより、固有ひずみの未知数が分布関数パラメータに集約されるため、残留応力の測定値が少なくても精度よく残留応力を導出できる。
　また、金属材１に対するＸ線応力測定における回折Ｘ線の半価幅、金属材１の硬さ、類似事例及び有限要素解析は、金属材１の内部に存在する塑性ひずみを推測可能とする。最適化問題では、設計変数の初期値が最適解に近いほど設計変数の最適解への収束が速くなる。このため、固有ひずみを付与可能なパラメータは、金属材１に対するＸ線応力測定における回折Ｘ線の半価幅、金属材１の硬さ、類似事例及び有限要素解析の少なくともいずれかに基づき初期設定されると好ましい。パラメータが塑性ひずみの推測によって初期設定されると、最適化問題を解く効率が向上する。
　なお、各要素に初期設定されるパラメータは、全て同じ値であってもよいし、一部又は全部が異なる値であってもよい。また、パラメータは、任意の値に初期設定されてもよく、塑性ひずみの推測によって初期設定されなくてもよい。
＜導出工程Ｓ４＞
　導出工程Ｓ４は、設定工程Ｓ３でパラメータを設定された複数の有限要素モデルのそれぞれについて、有限要素解析に基づき複数箇所のそれぞれにおける残留応力の計算値を導出する工程である。導出工程Ｓ４では、情報処理装置が、複数の有限要素モデルのそれぞれについて、金属材１における残留応力の測定箇所に対応する複数箇所のそれぞれにおける残留応力の計算値を有限要素解析に基づき導出する。この有限要素解析には、周知の手法を用いることができる。このため、情報処理装置は、例えば市販の有限要素解析ソフトウェアに基づき有限要素解析を実行できる。
＜算出工程Ｓ５＞
　算出工程Ｓ５は、導出工程Ｓ４で導出された残留応力の計算値及び用意工程Ｓ１で用意された残留応力の測定値の差を複数箇所のそれぞれについて算出する工程である。算出工程Ｓ５では、情報処理装置が、金属材１における残留応力の測定箇所のそれぞれについて、有限要素モデルを用いて導出した残留応力の計算値と記憶部に記憶された残留応力の測定値との差を算出する。
　算出工程Ｓ５では、測定値の誤差範囲に基づき計算値及び測定値の差が算出される。具体的には、情報処理装置が、計算値及び測定値の差の算出において、残留応力の計算値が残留応力の測定値の誤差範囲内となる場合は、計算値及び測定値の差をゼロとし、計算値が測定値の誤差範囲外となる場合は、計算値から測定値の誤差範囲の境界までの差を計算値及び測定値の差とするように計算を実行する。このため、残留応力の測定点の数が非常に多い場合であっても、適切に測定誤差を考慮した計算結果が得られる。なお、測定値の誤差範囲としては、例えば測定箇所毎の標準偏差が用いられるとよい。
〔最適化問題〕
　当該残留応力算出方法は、算出工程Ｓ５で算出される残留応力の計算値と残留応力の測定値との差を最小化する最適化問題として、パラメータを調整しつつ設定工程Ｓ３、導出工程Ｓ４及び算出工程Ｓ５を繰り返し実行することによりパラメータの最適解を求める。具体的には、情報処理装置が、残留応力の計算値と残留応力の測定値との差を最小化することを目的関数として設定し、パラメータを設定変数として設定することにより、パラメータを調整しつつ有限要素モデルに基づき上記差を繰り返し算出し、上記差を最小化するようなパラメータの最適解を求める。最適化探索には、例えば離散型探索手法であるＨｏｏｋｅ－Ｊｅｅｖｅｓ法等の周知の手法が用いられる。このため、情報処理装置は、市販の最適化ソフトウェアに基づき最適化問題の解析を実行できる。
　なお、目的関数としては、例えば残留応力の計算値と残留応力の測定値との差の二乗和が用いられてもよいし、残留応力の計算値と残留応力の測定値との差の最大値が用いられてもよい。また、パラメータの最適解が得られたか否かの判定については、最適化問題の解析において、例えば上記差が大きく変化しないまま一定時間経過したか否かを判定することにより実施し、一定時間経過した場合にパラメータの最適解が得られたと判定するようにしてもよい。
　最適化問題を解く際のパラメータの調整範囲は、金属材１に対するＸ線応力測定における回折Ｘ線の半価幅、金属材１の硬さ、類似事例及び有限要素解析の少なくともいずれかに基づく制約条件により設定されるとよい。例えば、類似の検討事例があり、固有ひずみ分布のばらつきが把握できている場合、そのばらつき範囲に基づく制約条件によりパラメータの調整範囲を設定してもよい。また例えば、金属材１における残留応力の発生起因となるプロセスを有限要素解析にて再現した結果、金属材１の硬さの測定結果、金属材１に対するＸ線応力測定における回折Ｘ線の半価幅の測定結果のいずれかに基づき固有ひずみの最大値及び最小値を予測し、予測結果に応じた制約条件によりパラメータの調整範囲を設定してもよい。また、パラメータの調整範囲は、最大値及び最小値の少なくともいずれかが制約条件により設定されるとよい。パラメータに合理的な制約条件が課されると、最適化問題を解く効率が向上する。
　また、当該残留応力算出方法は、最適化問題を解く際の制約条件の一つとして、金属材１の体積が固有ひずみの生成によって変化しないという条件を用いるとよい。金属材１の残留応力が塑性ひずみによって生じる場合、塑性変形した部分の形状は変化するものの金属材１の体積は変化していないとみなせる。したがって、当該残留応力算出方法は、この制約条件を用いることで、最適化問題における設計変数を減らすことができる。
　最適化問題の解析には、図４に示すようにＮ個の有限要素モデルが用いられる（Ｎは２以上の自然数）。上述したように、まず設定工程Ｓ３で、Ｎ個の有限要素モデルのそれぞれに対して、例えば固有ひずみ値がパラメータとして設定される。次に導出工程Ｓ４で、Ｎ個の有限要素モデルのそれぞれについて、金属材１における残留応力の測定箇所に対応する複数箇所のそれぞれにおける残留応力の計算値が有限要素解析に基づき導出される。残留応力は複数箇所のそれぞれについて計算されるため、残留応力の計算値は有限要素モデル毎に残留応力群として対応付けられる。次に算出工程Ｓ５で、複数箇所のそれぞれについて、残留応力の計算値及び測定値の差が算出される。算出された差は有限要素モデル毎に差群として対応付けられる。次に、算出された全ての差の合計が最小であるか否かが判定される。差の合計が最小でない場合又は算出工程Ｓ５の実行が初回である場合は、探索手法に基づき固有ひずみ値を調整した後に設定工程Ｓ３に戻り、調整された固有ひずみ値がＮ個の有限要素モデルに再設定される。そして差の合計が最小となるまでこれら処理が繰り返される。
＜残留応力の算出Ｓ６＞
　当該残留応力算出方法は、最適化問題を解くことによりパラメータの最適解を求めた後、この最適解となるパラメータを切削前の金属材１と同一形状の有限要素モデルに設定し、この有限要素モデルについて有限要素解析に基づき残留応力を計算する。具体的には、情報処理装置が、上述の最適化問題を解くことによりパラメータの最適解を求めた後、この最適解となるパラメータを切削前の金属材１と同一形状の有限要素モデルに設定する。そして情報処理装置が、この有限要素モデルについて有限要素解析に基づき任意の位置における残留応力を計算する。これにより、当該残留応力算出方法は、金属材１の任意の位置における残留応力の推測値を算出できる。
（利点）
　当該残留応力算出方法は、切削後の金属材１と同一形状の複数の有限要素モデルのそれぞれについて、有限要素解析に基づき残留応力の計算値を導出する。そして、当該残留応力算出方法は、残留応力の計算値と測定値との差を最小化する最適化問題を解くことにより、固有ひずみを付与するパラメータの最適解を求め、この最適解となるパラメータを用いて切削前の金属材１と同一形状の有限要素モデルにおける残留応力を計算する。したがって、当該残留応力算出方法は、固有ひずみと弾性ひずみとの関係を表す係数マトリックスを有限要素法で求めることなく、簡単に残留応力を算出できる。
　当該残留応力算出方法は、金属材１のフィレット面Ｘを含む範囲を残留応力の算出対象とするので、金属材１の内部に生じる残留応力を適切に評価できる。また、フィレット面Ｘが円弧状に形成された領域を含み、金属材１の切削面の延長面がフィレット面Ｘの円弧の中心位置を通るので、当該残留応力算出方法は、フィレット面Ｘの円弧の中心位置を局所座標系の中心として定義することで、有限要素解析を効率よく実施できる。
　当該残留応力算出方法は、算出工程Ｓ５で、残留応力の計算値が測定値の誤差範囲内となる場合は、計算値及び測定値の差をゼロとし、計算値が測定値の誤差範囲外となる場合は、計算値から誤差範囲の境界までの差を計算値及び測定値の差とするので、測定点の数が非常に多い場合であっても、適切に測定誤差を考慮できる。
＜その他の実施形態＞
　今回開示された実施の形態はすべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は、上記実施形態の構成に限定されるものではなく、請求の範囲によって示され、請求の範囲と均等の意味及び範囲内での全ての変更が含まれることが意図される。
　上記実施形態では、切削後の金属材１と同一形状の複数の有限要素モデルのそれぞれについて、有限要素解析に基づき残留応力の計算値を導出し、残留応力の計算値と測定値との差を最小化する最適化問題を解くことにより、固有ひずみを付与するパラメータの最適解を求めるものについて説明したが、切削後の金属材１と同一形状の複数の有限要素モデルのそれぞれについて、有限要素解析に基づき弾性ひずみの計算値を導出し、弾性ひずみの計算値と測定値との差を最小化する最適化問題を解くことにより、固有ひずみを付与するパラメータの最適解を求めるものであってもよい。つまり、残留応力算出方法は、固有ひずみ法を用いた金属材１の残留応力算出方法であって、金属材１の表面上の複数箇所及びこの金属材１を切削することにより形成される切削面上の複数箇所のそれぞれにおける弾性ひずみの測定値を用意する用意工程と、金属材１の表面の一部及び切削面を含む切削後の金属材１と同一形状の複数の有限要素モデルを作成する作成工程と、作成工程で作成された複数の有限要素モデルのそれぞれに対し、固有ひずみを付与可能なパラメータを設定する設定工程と、設定工程でパラメータを設定された複数の有限要素モデルのそれぞれについて、有限要素解析に基づき複数箇所のそれぞれにおける弾性ひずみの計算値を導出する導出工程と、導出工程で導出された弾性ひずみの計算値及び用意工程で用意された弾性ひずみの測定値の差を複数箇所のそれぞれについて算出する算出工程とを備え、算出工程で算出される計算値と測定値との差を最小化する最適化問題として、パラメータを調整しつつ設定工程、導出工程及び算出工程を繰り返し実行することによりパラメータの最適解を求めた後、この最適解となるパラメータを切削前の金属材１と同一形状の有限要素モデルに設定し、この有限要素モデルについて有限要素解析に基づき残留応力を計算するものであってもよい。
　上記実施形態では、金属材１が、円柱状の軸部２とこの軸部２から全周かつ径方向に突出する４つの板状部３とを備え、軸部２と板状部３との接続部分にフィレット面Ｘを有するものについて説明したが、残留応力算出方法に用いられる金属材は、このような構造のものに限定されない。例えば、金属材がフィレット面を有していなくてもよい。
　上記実施形態では、金属材１のフィレット面Ｘが円弧状に形成された領域を含み、金属材１の残留応力を測定する際に、金属材１の切削面の延長面がフィレット面Ｘの円弧の中心位置を通るように金属材１が切削されるものについて説明したが、金属材１の切削面の延長面がフィレット面Ｘの円弧の中心位置を通らないものであってもよい。ただし、有限要素解析を効率よく実施する観点から、金属材１の切削面の延長面がフィレット面Ｘの円弧の中心位置を通るように金属材１が切削されると好ましい。
　上記実施形態では、算出工程Ｓ５で、残留応力の計算値が残留応力の測定値の誤差範囲内となる場合は、計算値及び測定値の差をゼロとし、計算値が測定値の誤差範囲外となる場合は、計算値から測定値の誤差範囲の境界までの差を計算値及び測定値の差とするように計算が実行されるものについて説明したが、測定誤差を考慮せずに計算値及び測定値の差を算出してもよい。ただし、適切に測定誤差を考慮した計算結果を得る観点から、測定値の誤差範囲に基づき計算値及び測定値の差を算出すると好ましい。
　上記実施形態では、作成工程Ｓ２の前に用意工程Ｓ１を実行するものについて説明したが、用意工程Ｓ１は算出工程Ｓ５より前の任意のタイミングで実行されればよい。
　以下、実施例によって本発明をさらに詳細に説明するが、本発明はこれらの実施例に限定されるものではない。
［残留応力の測定］
　図２に示す金属材について、フィレット面及びフィレット面より内部の残留応力を測定した。金属材の材料としては鋼材を用いた。４つのフィレット面には、同様に冷間加工が施されており、フィレット面は、直径２２ｍｍの円弧状に形成された領域を含んでいる。簡略化のために、４つのフィレット面は同等であり、かつ軸部の周方向には均一に残留応力が生じていると仮定する。また、直径２２ｍｍの円弧の中心位置を局地座標系の原点とする。
　３方向の残留応力を全て測定できるように、金属材からＦ片、Ｃ片、Ｔ片の３種類の形状の切削片を用意した。Ｆ片、Ｃ片、Ｔ片は、いずれも切削面の延長面が局地座標系の原点を通るように形成される。Ｆ片は、θ＝０度で板状部側から軸部を切り落とした形状であり、θを変更しながらのフィレット面の表面付近のθ方向及びφ方向の残留応力の測定が可能である。Ｃ片は、θ＝０度からθ＝１１０度の範囲で軸部側から板状部を円錐形状等に切り落とした形状であり、θを変更しながらのフィレット面より内部のｒ方向及びφ方向の残留応力の測定が可能である。Ｔ片は、軸部の中心線を通るようにφをわずかに変更して金属材をスライスした薄片の形状であり、θを変更しながらのフィレット面より内部のｒ方向及びθ方向の残留応力の測定が可能である。Ｆ片及びＴ片については１つずつ用意し、Ｃ片についてはθ＝０度から２０度ずつ変更してθ＝１００度までを測定するもの（Ｃａ片）と、θ＝１０度から２０度ずつ変更してθ＝１１０度まで測定するもの（Ｃｂ片）の２つを用意した。なお、残留応力の測定位置には電解研磨処理を施した。
　Ｆ片については、Ｘ線応力測定装置（リガク社のＭＳＦ－３Ｍ）を用いてフィレット表面から深さ０．３ｍｍ及び０．５ｍｍの位置の残留応力を測定した。残留応力の測定にはｓｉｎ２ψ法を用いた。θを８０度から１００度の範囲で１０度ずつ変更しながらθ方向応力を測定し、θを０度から１００度の範囲で１０度ずつ変更しながらφ方向応力を測定して、合計２８点の測定値を得た。
　２つのＣ片については、Ｘ線応力測定装置（リガク社のＭＳＦ－３Ｍ）を用いてフィレット表面から深さ１８ｍｍから４０ｍｍの位置の残留応力を測定した。残留応力の測定にはｓｉｎ２ψ法を用いた。Ｃａ片では、切削によりθを０度から１００度の範囲で２０度ずつ変更しながら、ｒ方向応力とφ方向応力とを１０から１６点ずつ測定し、Ｃｂ片では、切削によりθを１０度から１１０度の範囲で２０度ずつ変更しながら、ｒ方向応力とφ方向応力とを１０から１６点ずつ測定して、合計３０８点の測定値を得た。
　Ｔ片については、Ｘ線残留応力測定装置（パルステック社のμ－３６０）を用いてフィレット表面から深さ１８ｍｍから４０ｍｍの位置の残留応力を測定した。残留応力の測定にはｃｏｓα法を用いた。θを０度から１００度の範囲で１０度ずつ変更しながらｒ方向応力とθ方向応力とを１０から１６点ずつ測定して、合計２２２点の測定値を得た。
［有限要素モデルの準備］
　情報処理装置において、有限要素法解析ソフトを用いてＦ片、Ｃ片、Ｔ片と同一形状の有限要素モデルをそれぞれ作成した。なお、Ｃ片については、切削によりθを変更しているので、残留応力を測定した際の形状とそれぞれ同一形状の１２種類の有限要素モデルを作成した。
　固有ひずみを分布関数で近似する手法を採用するため、作成した有限要素モデルに対し、次式の分布関数によって固有ひずみを設定した。
Figure JPOXMLDOC01-appb-M000010
Figure JPOXMLDOC01-appb-M000011
 　ここで、Ｒはｒ方向を示し、αはθ方向を示す。Ｒ０及びα０は固有ひずみ計算の開始位置を示し、ΔＲ及びΔαは増分値を示す。また、ｍ及びｎは、それぞれｒ方向及びθ方向の次数を示し、Ａｉｊは近似関数の係数を示す。Ｒ０、α０、ΔＲ及びΔαは、固有ひずみ計算領域を決定するパラメータであり、これらの設定値が異なれば、同じ残留応力データを用いても、得られる固有ひずみ分布に差異が生じる。同様に、固有ひずみ分布関数の次数ｍ及びｎの値によっても、得られる固有ひずみ分布は異なる。
　図２に示す金属材は軸対称であるので、φ方向のせん断ひずみは０となる。したがって、４つの成分のひずみのみを考慮すればよい。これらの成分のひずみを計算するために用いる係数Ａｉｊをそれぞれａｉｊ、ｂｉｊ、ｃｉｊ、ｄｉｊとした。
　相当塑性ひずみは、図５に示すように、Ｘ線残留応力測定における半価幅（Ｘ線半価幅）との間で相関を持つことがわかっている。また、固有ひずみが塑性ひずみのみによって生じているとみなすと、固有ひずみはＸ線半価幅との間で相関を持つ。そこで、固有ひずみの近似関数の係数について、Ｘ線半価幅の測定結果に基づき初期値を決定した。ｍを２とし、ｎを８とすると、ａｉｊ、ｂｉｊ、ｃｉｊ、ｄｉｊのそれぞれが１６個の係数となる。この条件においてＸ線半価幅の測定結果に基づき決定した近似関数の係数の初期値を表１に示す。
Figure JPOXMLDOC01-appb-T000012
［最適化問題の解析］
　設計変数を固有ひずみの近似関数の係数ａｉｊ、ｂｉｊ、ｃｉｊ、ｄｉｊとし、目的関数を残留応力の測定値と計算値との差を最小化することとして、最適化問題を解析した。最適化問題における制約条件としては係数ａｉｊ、ｂｉｊ、ｃｉｊ、ｄｉｊの絶対値を０．１未満とした。また、最適化計算には汎用の解析ソフト（ダッソー・システムズ社のＩｓｉｇｈｔ）を用い、最適化探索にはＨｏｏｋｅ－Ｊｅｅｖｅｓ法を採用した。
　まず、分布関数によって固有ひずみを設定した有限要素モデルを用いて、実際に残留応力を測定した位置に対応する位置の残留応力を有限要素解析により計算した。次に、各有限要素モデルについて、位置毎に、残留応力の測定値と計算値との差の二乗を算出し、これらを全て合計した。そして、差の二乗の合計が最小となるようにＨｏｏｋｅ－Ｊｅｅｖｅｓ法により係数ａｉｊ、ｂｉｊ、ｃｉｊ、ｄｉｊの値を調整しつつこれらの計算を繰り返した。差の二乗の合計が一定時間の計算において大きく変化しなくなった時点で、差の二乗の合計が最小となったと判定した。
［残留応力の算出］
　切削前の図２の金属材と同一形状の有限要素モデルに対し、係数ａｉｊ、ｂｉｊ、ｃｉｊ、ｄｉｊの最適解を設定し、この有限要素モデルについて有限要素解析に基づきθ＝３０度における深さに対する残留応力を計算した。また、比較のために、ひずみゲージを用いた切断法により図２の金属材のθ＝３０度における深さに対する残留応力を測定した。残留応力の測定方向はθ方向とした。この切断法では、ドリルを用いて測定部位が露出されるように孔を設け、この孔の底部にひずみゲージを貼り、このひずみゲージを添付した部分を含むように金属材を最小片に切断し、切断前後の弾性解放ひずみを測定することにより残留応力を把握する。
　また、残留応力の測定値と計算値との差の算出において、計算値が測定値の誤差範囲内となる場合は、計算値及び測定値の差をゼロとし、計算値が測定値の誤差範囲外となる場合は、計算値から誤差範囲の境界までの差を計算値及び測定値の差とする誤差範囲を考慮した算出についても検討した。測定値の誤差範囲としては、測定箇所毎の標準偏差を用いた。
　さらに、最適化問題を解く際の制約条件の一つとして、金属材の体積が固有ひずみの生成によって変化しないという条件を用いる解析手法についても検討した。体積が変化しないという条件により、ａｉｊ、ｂｉｊ、ｃｉｊの和が０となるので、設計変数を減らすことができた。
　誤差範囲を考慮せずに残留応力を計算した結果（Ｃ１）、誤差範囲を考慮して残留応力を計算した結果（Ｃ２）、体積が変化しないという制約条件を考慮して残留応力を計算した結果（Ｃ３）及び切断法により残留応力を測定した結果（Ｍｅ）をそれぞれ図６に示す。
　図６に示すように、残留応力の計算結果は、いずれも切断法による残留応力の測定結果とよく一致することが確認された。特に、誤差範囲を考慮して残留応力を計算した結果が、残留応力の測定結果と最も一致しているといえる。また、体積が変化しないという制約条件を考慮した残留応力の計算は、設計変数を減らしているにもかかわらず、残留応力の測定結果と十分に一致しているといえる。
　本発明の残留応力算出方法は、固有ひずみと弾性ひずみとの関係を表す係数マトリックスを有限要素法で求めることなく、簡単に残留応力を算出できる。
　１　金属材
　２　軸部
　３　板状部
　Ｘ　フィレット面

Claims (10)

  1. 　固有ひずみ法を用いた金属材の残留応力算出方法であって、
    　金属材の表面上の複数箇所及びこの金属材を切削することにより形成される切削面上の複数箇所のそれぞれにおける残留応力の測定値を用意する工程と、
    　上記表面の一部及び上記切削面を含む切削後の上記金属材と同一形状の複数の有限要素モデルを作成する工程と、
    　上記作成工程で作成された複数の上記有限要素モデルのそれぞれに対し、固有ひずみを付与可能なパラメータを設定する工程と、
    　上記設定工程で上記パラメータを設定された複数の上記有限要素モデルのそれぞれについて、有限要素解析に基づき上記複数箇所のそれぞれにおける残留応力の計算値を導出する工程と、
    　上記導出工程で導出された残留応力の上記計算値及び上記用意工程で用意された残留応力の上記測定値の差を上記複数箇所のそれぞれについて算出する工程と
    　を備え、
    　上記算出工程で算出される上記計算値と上記測定値との差を最小化する最適化問題として、上記パラメータを調整しつつ上記設定工程、上記導出工程及び上記算出工程を繰り返し実行することにより上記パラメータの最適解を求めた後、この最適解となるパラメータを切削前の上記金属材と同一形状の有限要素モデルに設定し、この有限要素モデルについて有限要素解析に基づき残留応力を計算する残留応力算出方法。
  2. 　固有ひずみ法を用いた金属材の残留応力算出方法であって、
    　金属材の表面上の複数箇所及びこの金属材を切削することにより形成される切削面上の複数箇所のそれぞれにおける弾性ひずみの測定値を用意する工程と、
    　上記表面の一部及び上記切削面を含む切削後の上記金属材と同一形状の複数の有限要素モデルを作成する工程と、
    　上記作成工程で作成された複数の上記有限要素モデルのそれぞれに対し、固有ひずみを付与可能なパラメータを設定する工程と、
    　上記設定工程で上記パラメータを設定された複数の上記有限要素モデルのそれぞれについて、有限要素解析に基づき上記複数箇所のそれぞれにおける弾性ひずみの計算値を導出する工程と、
    　上記導出工程で導出された弾性ひずみの上記計算値及び上記用意工程で用意された弾性ひずみの上記測定値の差を上記複数箇所のそれぞれについて算出する工程と
    　を備え、
    　上記算出工程で算出される上記計算値と上記測定値との差を最小化する最適化問題として、上記パラメータを調整しつつ上記設定工程、上記導出工程及び上記算出工程を繰り返し実行することにより上記パラメータの最適解を求めた後、この最適解となるパラメータを切削前の上記金属材と同一形状の有限要素モデルに設定し、この有限要素モデルについて有限要素解析に基づき残留応力を計算する残留応力算出方法。
  3. 　上記金属材が円柱状の軸部とこの軸部から全周かつ径方向に突出する板状部とを備え、
    　上記金属材の上記表面が上記軸部と上記板状部との接続部分に形成されるフィレット面を含む請求項１に記載の残留応力算出方法。
  4. 　上記金属材が円柱状の軸部とこの軸部から全周かつ径方向に突出する板状部とを備え、
    　上記金属材の上記表面が上記軸部と上記板状部との接続部分に形成されるフィレット面を含む請求項２に記載の残留応力算出方法。
  5. 　上記フィレット面が円弧状に形成された領域を含み、
    　上記切削面の延長面が上記フィレット面の円弧の中心位置を通る請求項３又は請求項４に記載の残留応力算出方法。
  6. 　上記設定工程で、上記パラメータが、上記金属材に対するＸ線応力測定における回折Ｘ線の半価幅、上記金属材の硬さ、類似事例及び有限要素解析の少なくともいずれかに基づき初期設定される請求項１から請求項４のいずれか１項に記載の残留応力算出方法。
  7. 　上記最適化問題を解く際の上記パラメータの調整範囲が、上記金属材に対するＸ線応力測定における回折Ｘ線の半価幅、上記金属材の硬さ、類似事例及び有限要素解析の少なくともいずれかに基づく制約条件により設定される請求項１から請求項４のいずれか１項に記載の残留応力算出方法。
  8. 　上記算出工程で、上記計算値が上記測定値の誤差範囲内となる場合は、上記計算値及び上記測定値の差をゼロとし、上記計算値が上記測定値の誤差範囲外となる場合は、上記計算値から上記誤差範囲の境界までの差を上記計算値及び上記測定値の差とする請求項１から請求項４のいずれか１項に記載の残留応力算出方法。
  9. 　上記最適化問題を解く際の制約条件の一つとして、金属材の体積が固有ひずみの生成によって変化しないという条件を用いる請求項１から請求項４のいずれか１項に記載の残留応力算出方法。
  10. 　上記パラメータが、固有ひずみ又は固有ひずみ近似関数の係数である請求項１から請求項４のいずれか１項に記載の残留応力算出方法。
PCT/JP2019/012534 2018-04-25 2019-03-25 残留応力算出方法 WO2019208061A1 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
CN201980027485.4A CN112005088B (zh) 2018-04-25 2019-03-25 残留应力算出方法
FIEP19791476.5T FI3786600T3 (fi) 2018-04-25 2019-03-25 Menetelmä jäännösjännityksen laskemiseksi
EP19791476.5A EP3786600B1 (en) 2018-04-25 2019-03-25 Method for computing residual stress
ES19791476T ES2941730T3 (es) 2018-04-25 2019-03-25 Método para calcular la tensión residual
KR1020207033511A KR102475949B1 (ko) 2018-04-25 2019-03-25 잔류 응력 산출 방법

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2018083986A JP6867329B2 (ja) 2018-04-25 2018-04-25 残留応力算出方法
JP2018-083986 2018-04-25

Publications (1)

Publication Number Publication Date
WO2019208061A1 true WO2019208061A1 (ja) 2019-10-31

Family

ID=68294996

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2019/012534 WO2019208061A1 (ja) 2018-04-25 2019-03-25 残留応力算出方法

Country Status (7)

Country Link
EP (1) EP3786600B1 (ja)
JP (1) JP6867329B2 (ja)
KR (1) KR102475949B1 (ja)
CN (1) CN112005088B (ja)
ES (1) ES2941730T3 (ja)
FI (1) FI3786600T3 (ja)
WO (1) WO2019208061A1 (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111141437A (zh) * 2020-01-02 2020-05-12 中车青岛四方机车车辆股份有限公司 电阻点焊接头内部残余应力测量方法
CN113127999A (zh) * 2021-05-13 2021-07-16 江南造船(集团)有限责任公司 双壁管法兰的应力和强度评估方法
CN114878046A (zh) * 2022-04-22 2022-08-09 北京工业大学 一种测量厚板焊接件内部残余应力的方法
CN115392075A (zh) * 2022-08-11 2022-11-25 北京航空航天大学 一种基于多次走刀的车削残余应力二维仿真方法
EP4155703A4 (en) * 2020-06-02 2023-11-01 Kabushiki Kaisha Kobe Seiko Sho (Kobe Steel, Ltd.) RESIDUAL STRESS MEASUREMENT METHOD

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7566698B2 (ja) * 2020-09-10 2024-10-15 株式会社神戸製鋼所 測定システム及び測定方法
CN114659687B (zh) * 2020-12-23 2024-07-05 财团法人金属工业研究发展中心 曲面块材的残留应力量测方法
CN113297772B (zh) * 2021-06-07 2022-04-01 南昌航空大学 一种使用有限元软件模拟复合材料残余应力分布的方法
CN113378442B (zh) * 2021-06-30 2022-05-10 湖南大学 一种铝合金锻件残余应力的表征方法
CN113504435B (zh) * 2021-07-20 2022-07-08 中南大学 三电平逆变器开路故障诊断方法及系统
CN115479962B (zh) * 2022-09-21 2024-08-30 福州大学 快速准确推算焊接接头内部与内表面残余应力的方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005181172A (ja) 2003-12-22 2005-07-07 Mitsubishi Heavy Ind Ltd 残留応力測定方法及び装置
JP2008058179A (ja) * 2006-08-31 2008-03-13 Tokyo Institute Of Technology 残留応力の評価方法
JP2015145834A (ja) * 2014-02-03 2015-08-13 国立大学法人横浜国立大学 残留応力推定方法、残留応力推定システムおよびプログラム
WO2016140094A1 (ja) * 2015-03-05 2016-09-09 株式会社神戸製鋼所 残留応力推定方法及び残留応力推定装置
WO2016140093A1 (ja) * 2015-03-05 2016-09-09 株式会社神戸製鋼所 残留応力推定方法及び残留応力推定装置
WO2016140091A1 (ja) * 2015-03-05 2016-09-09 株式会社神戸製鋼所 残留応力推定方法及び残留応力推定装置

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102005008281B4 (de) * 2005-02-23 2013-04-04 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Verfahren zur Erfassung von lokalen Eigenspannungen in Festkörperobjekten mit einer Ionenstrahltechnik
JP5356894B2 (ja) * 2009-04-06 2013-12-04 ポリプラスチックス株式会社 残留応力算出方法及び残留応力分布導出方法
JP6060591B2 (ja) * 2012-09-26 2017-01-18 Ｊｆｅスチール株式会社 プレス成形品のスプリングバック要因特定方法および装置
CN103853899B (zh) * 2013-05-22 2017-02-15 中国能建集团装备有限公司南京技术中心 轴类零件疲劳寿命计算方法
JP5955301B2 (ja) * 2013-11-14 2016-07-20 株式会社神戸製鋼所 残留応力算出方法
CN107273630B (zh) * 2017-06-28 2020-05-19 华中科技大学 一种通过参数反演计算调控残余应力加工参数的方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005181172A (ja) 2003-12-22 2005-07-07 Mitsubishi Heavy Ind Ltd 残留応力測定方法及び装置
JP2008058179A (ja) * 2006-08-31 2008-03-13 Tokyo Institute Of Technology 残留応力の評価方法
JP2015145834A (ja) * 2014-02-03 2015-08-13 国立大学法人横浜国立大学 残留応力推定方法、残留応力推定システムおよびプログラム
WO2016140094A1 (ja) * 2015-03-05 2016-09-09 株式会社神戸製鋼所 残留応力推定方法及び残留応力推定装置
WO2016140093A1 (ja) * 2015-03-05 2016-09-09 株式会社神戸製鋼所 残留応力推定方法及び残留応力推定装置
WO2016140091A1 (ja) * 2015-03-05 2016-09-09 株式会社神戸製鋼所 残留応力推定方法及び残留応力推定装置

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111141437A (zh) * 2020-01-02 2020-05-12 中车青岛四方机车车辆股份有限公司 电阻点焊接头内部残余应力测量方法
CN111141437B (zh) * 2020-01-02 2021-12-21 中车青岛四方机车车辆股份有限公司 电阻点焊接头内部残余应力测量方法
EP4155703A4 (en) * 2020-06-02 2023-11-01 Kabushiki Kaisha Kobe Seiko Sho (Kobe Steel, Ltd.) RESIDUAL STRESS MEASUREMENT METHOD
CN113127999A (zh) * 2021-05-13 2021-07-16 江南造船(集团)有限责任公司 双壁管法兰的应力和强度评估方法
CN113127999B (zh) * 2021-05-13 2023-03-17 江南造船(集团)有限责任公司 双壁管法兰的应力和强度评估方法
CN114878046A (zh) * 2022-04-22 2022-08-09 北京工业大学 一种测量厚板焊接件内部残余应力的方法
CN114878046B (zh) * 2022-04-22 2023-11-24 北京工业大学 一种测量厚板焊接件内部残余应力的方法
CN115392075A (zh) * 2022-08-11 2022-11-25 北京航空航天大学 一种基于多次走刀的车削残余应力二维仿真方法

Also Published As

Publication number Publication date
EP3786600A1 (en) 2021-03-03
ES2941730T3 (es) 2023-05-25
CN112005088B (zh) 2021-11-26
JP6867329B2 (ja) 2021-04-28
KR102475949B1 (ko) 2022-12-08
FI3786600T3 (fi) 2023-03-22
JP2019190990A (ja) 2019-10-31
KR20210002595A (ko) 2021-01-08
EP3786600B1 (en) 2023-03-15
CN112005088A (zh) 2020-11-27

Similar Documents

Publication Publication Date Title
WO2019208061A1 (ja) 残留応力算出方法
JP6345617B2 (ja) 残留応力推定方法及び残留応力推定装置
CN107430637B (zh) 残余应力推定方法及残余应力推定装置
EP3070447B1 (en) Residual stress measuring method
KR101949070B1 (ko) 잔류 응력 추정 방법 및 잔류 응력 추정 장치
JP6560998B2 (ja) 残留応力算出方法
Sinha et al. Vibratory parameters of blades from coordinate measurement machine data
CN105136361A (zh) 一种用x射线衍射测定立方单晶体材料应力的方法
JP4533621B2 (ja) 残留応力測定方法及び装置
JP2008058179A (ja) 残留応力の評価方法
JPH05223661A (ja) 残留応力の測定方法
JP2021072121A (ja) Ｃｎｃ機械加工のための適応性経路生成
CN116735060A (zh) 一种航空铝合金t形锻件残余应力的测试方法
KR20080096873A (ko) 1차원 등가 열팽창 계수를 이용한 곡 빌트업 자동 역 설계방법
Todorovic et al. Comparative model analysis of two types of clamping elements in dynamic conditions
JP6376081B2 (ja) 被検体の形状算出方法、鋼板のプレス矯正方法および鋼板の製造方法
JP4269678B2 (ja) 計測データ比較装置および形状品質評価システム
US20240033801A1 (en) Working Method of Material, and Process Design Computer and Program of the Same
Park et al. Application of localization for the fabrication of large curved plates in shipbuilding
PAVLASEK et al. SOFTWARE TOOL FOR CONVERSION OF HISTORICAL TEMPERATURE SCALES
BR102020009604A2 (pt) métodos para analisar um projeto de peça
CN118195551A (zh) 一种地下装配式结构施工安全预警方法、系统及介质
Todorovic et al. Usporedna analiza modela dva tipa elemenata za stezanje u dinamičkim uvjetima

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 19791476

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 20207033511

Country of ref document: KR

Kind code of ref document: A

ENP Entry into the national phase

Ref document number: 2019791476

Country of ref document: EP

Effective date: 20201125