JP2023005230A - Damage behavior evaluation method of micro-structure of metallic material - Google Patents

Damage behavior evaluation method of micro-structure of metallic material Download PDF

Info

Publication number
JP2023005230A
JP2023005230A JP2021107010A JP2021107010A JP2023005230A JP 2023005230 A JP2023005230 A JP 2023005230A JP 2021107010 A JP2021107010 A JP 2021107010A JP 2021107010 A JP2021107010 A JP 2021107010A JP 2023005230 A JP2023005230 A JP 2023005230A
Authority
JP
Japan
Prior art keywords
strain
microstructure
stress
plastic strain
fem analysis
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2021107010A
Other languages
Japanese (ja)
Other versions
JP7474225B2 (en
Inventor
翔太 椿
Shota Tsubaki
瑛介 黒澤
Eisuke KUROSAWA
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Kobe Steel Ltd
Original Assignee
Kobe Steel Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Kobe Steel Ltd filed Critical Kobe Steel Ltd
Priority to JP2021107010A priority Critical patent/JP7474225B2/en
Publication of JP2023005230A publication Critical patent/JP2023005230A/en
Application granted granted Critical
Publication of JP7474225B2 publication Critical patent/JP7474225B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

To provide an evaluation method of a damage behavior of a micro-structure of a metallic material in conformity to a real condition targeting the micro-structure of the metallic material to be evaluated without using a complicated numerical analysis model.SOLUTION: A method for evaluating a damage behavior of a micro-structure of a metallic material using a computer executes the steps of: causing a computer to acquire a micro-structure image of the metallic material; creating an FEM analysis model simulating the micro-structure image; acquiring a material property as a parameter; setting a history of the stress or strain applied to each element of the FEM analysis model as a boundary condition; outputting equivalent plastic strain and stress tri-axial degree of each element by executing the FEM analysis; obtaining the distribution of the output equivalent plastic strain and stress tri-axial degree; and evaluating the damage behavior of the micro-structure from the distribution.SELECTED DRAWING: Figure 1

Description

本発明は、金属材料のミクロ組織の損傷挙動評価方法に関する。 The present invention relates to a method for evaluating damage behavior of microstructures of metallic materials.

金属材料の多くは、例えば製造過程の塑性加工時、該金属材料で形成された構造体を使用時に、大きな荷重を受けうる。大きな荷重を受けて、上記金属材料の破断ひずみ量を超える大きなひずみが発生した場合、金属材料は破断し、その結果、製造過程では製造不良となり、構造体として使用時にはその機能が失われる。よって、金属材料の破断特性を事前に把握すること、例えばあるひずみ発生下で材料が破断するかどうかを評価することは非常に重要である。 Many metal materials can receive a large load, for example, during plastic working in the manufacturing process or when using a structure formed of the metal material. When a large strain exceeding the breaking strain amount of the metal material is generated by receiving a large load, the metal material breaks, resulting in a manufacturing defect in the manufacturing process and loss of its function when used as a structure. Therefore, it is very important to know the rupture properties of metallic materials in advance, for example, to evaluate whether the material will rupture under the occurrence of a certain strain.

延性破壊を示す金属材料は、上記ひずみで変形中にミクロ組織内でマイクロボイドが発生し、それが成長および連結することによって最終的に破断に至る。本明細書では、上記マイクロボイドの発生(成長も含む)を「金属材料の損傷」ということがある。そのため、マイクロボイド発生・成長挙動を評価することができれば、その材料の、例えば破断伸び、局所延性等の破断特性を評価できると考えられる。上記マイクロボイドの発生には、材料に生じる塑性ひずみと静水圧応力が関係することが知られており、例えば非特許文献1および非特許文献2では、材料の破断時の塑性ひずみ量と応力三軸度(=静水圧応力/相当塑性ひずみ)の関係を定義して、有限要素法(FEM)等のシミュレーション上でその材料の破断を表現する、GISSMOモデルが使用されている。 Metallic materials exhibiting ductile fracture generate microvoids in the microstructure during deformation due to the above strain, which eventually lead to fracture due to their growth and coalescence. In this specification, the generation (including growth) of microvoids is sometimes referred to as "damage to metal materials." Therefore, if microvoid generation/growth behavior can be evaluated, the rupture properties of the material, such as elongation at break and local ductility, can be evaluated. It is known that the generation of microvoids is related to the plastic strain and hydrostatic stress generated in the material. A GISSMO model is used that defines the relationship of axiality (=hydrostatic stress/equivalent plastic strain) and expresses the fracture of the material on simulations such as the finite element method (FEM).

ところで上記金属材料の開発では、所望の機械的特性が得られるよう、合金成分の調整や熱処理による工夫等により、ミクロ組織の種類、組織形状、各組織の相分率等が様々な金属材料が製造される。しかし、上記非特許文献1および非特許文献2の技術では、上記様々なミクロ組織に応じて、マイクロボイド発生・成長挙動を明確に把握することが難しかった。 By the way, in the development of the above metal materials, metal materials with various types of microstructures, structure shapes, phase fractions of each structure, etc. have been developed by adjusting the alloy components and devising heat treatments so as to obtain the desired mechanical properties. manufactured. However, with the techniques of Non-Patent Document 1 and Non-Patent Document 2, it was difficult to clearly grasp the generation and growth behavior of microvoids according to the various microstructures.

一方、材料のマイクロボイド発生挙動を直接評価する方法が提案されている。例えば非特許文献3には、連続体損傷力学を適用したミクロ組織の解析モデルを用いて、応力解析により材料のマイクロボイド発生挙動を直接評価することが提案されている。 On the other hand, a method for directly evaluating the microvoid generation behavior of materials has been proposed. For example, Non-Patent Document 3 proposes direct evaluation of the microvoid generation behavior of a material by stress analysis using a microstructure analysis model to which continuum damage mechanics is applied.

F.Neukamm et al., 7th European LS-DYNA Conference,2009F. Neukamm et al. , 7th European LS-DYNA Conference, 2009 鎮西将太ら,R&D神戸製鋼技報,Vol.66,No.2(2017),page76-81Shota Chinzei et al., R&D Kobe Steel Technical Report, Vol. 66, No. 2 (2017), pages76-81 米村繁ら,新日鐵住金技報,第410号(2018),page47-56Shigeru Yonemura et al., Nippon Steel & Sumitomo Metal Technical Report, No. 410 (2018), page 47-56

従来、ミクロ組織の損傷挙動の評価には、非特許文献3の連続体損傷力学モデルのような複雑で計算負荷の高い計算モデルを使用する必要があった。しかし、該連続体損傷力学モデルは、決定の困難な材料パラメータを多数含んでいる。よって例えば複相鋼の場合、各構成相について該材料パラメータを全て決定することは難しく、経験的に想定して上記材料パラメータを設定せざるを得ず、実態に即した破断特性の評価は難しかった。 Conventionally, evaluation of microstructure damage behavior required the use of a complex and computationally intensive computational model such as the continuum damage mechanics model of Non-Patent Document 3. However, the continuum damage mechanics model contains many difficult-to-determine material parameters. Therefore, for example, in the case of dual-phase steel, it is difficult to determine all the material parameters for each constituent phase, and the above material parameters have to be set based on empirical assumptions, making it difficult to evaluate the rupture properties in line with the actual situation. rice field.

本発明は、上記事情に鑑みてなされたものであって、複雑な数値解析モデルを使用することなく、かつ評価する金属材料のミクロ組織を対象に、実態に即した金属材料のミクロ組織損傷挙動の評価方法を提供することにある。 The present invention has been made in view of the above circumstances, and does not use a complicated numerical analysis model, and targets the microstructure of the metal material to be evaluated. is to provide an evaluation method for

本発明の態様1は、
コンピュータを用いて金属材料のミクロ組織の損傷挙動を評価する方法であって、コンピュータが、
前記金属材料のミクロ組織画像を取得するステップ1と、
前記ミクロ組織画像を複数の要素に分割し、有限要素法(FEM)によって前記ミクロ組織画像を模擬したFEM解析モデルを作成するステップ2と、
前記FEM解析モデルのパラメータとして前記金属材料の材料特性を取得するステップ3と、
境界条件として、前記FEM解析モデルの各要素に付与する、前記金属材料の破断までに生じる応力若しくはひずみ、または前記応力若しくは前記ひずみの履歴を設定するステップ4と、
前記境界条件でFEM解析を実行し、前記各要素の相当塑性ひずみと応力三軸度を出力するステップ5と、
前記FEM解析において前記FEM解析モデルに与えた応力またはひずみ別に、出力された前記各要素の相当塑性ひずみと応力三軸度の分布を求めるステップ6と、
前記相当塑性ひずみと前記応力三軸度の分布から、ミクロ組織の損傷挙動を評価するステップ7と
を実行する、金属材料のミクロ組織の損傷挙動評価方法である。
Aspect 1 of the present invention is
A computer-aided method for evaluating microstructural damage behavior of a metallic material, the computer comprising:
a step 1 of obtaining a microstructure image of the metallic material;
step 2 of dividing the microstructure image into a plurality of elements and creating an FEM analysis model simulating the microstructure image by a finite element method (FEM);
Step 3 of obtaining material properties of the metallic material as parameters of the FEM analysis model;
Step 4 of setting, as boundary conditions, the stress or strain that occurs until the metal material fractures, or the history of the stress or strain that is applied to each element of the FEM analysis model;
Step 5 of performing an FEM analysis under the boundary conditions and outputting the equivalent plastic strain and stress triaxiality of each element;
Step 6 for obtaining the distribution of the equivalent plastic strain and the stress triaxiality of each of the output elements for each stress or strain applied to the FEM analysis model in the FEM analysis;
A method for evaluating damage behavior of a microstructure of a metal material, comprising: performing step 7 of evaluating damage behavior of a microstructure from the distribution of the equivalent plastic strain and the stress triaxiality.

本発明の態様2は、
前記ステップ4は、各要素の積分点に付与する、前記応力若しくは前記ひずみ、または前記応力若しくは前記ひずみの前記履歴を設定し、
前記ステップ5は、各要素の積分点における前記相当塑性ひずみと応力三軸度を出力する、態様1に記載の金属材料のミクロ組織の損傷挙動評価方法である。
Aspect 2 of the present invention is
The step 4 sets the stress or the strain, or the history of the stress or the strain, to be applied to the integration point of each element,
The step 5 is the damage behavior evaluation method of the microstructure of the metal material according to the aspect 1, wherein the equivalent plastic strain and the stress triaxiality at the integration point of each element are output.

本発明の態様3は、
前記金属材料は、2種類以上の構成相を有し、
前記ステップ5~7を構成相ごとに実施する、態様1または2に記載の金属材料のミクロ組織の損傷挙動評価方法である。
Aspect 3 of the present invention is
The metal material has two or more constituent phases,
3. A method for evaluating damage behavior of a microstructure of a metal material according to aspect 1 or 2, wherein steps 5 to 7 are performed for each constituent phase.

本発明の態様4は、
予め用意された、前記2種類以上の構成相それぞれの応力-ひずみ曲線を、前記ステップ3にて前記材料特性として取得することを含む、態様3に記載の金属材料のミクロ組織の損傷挙動評価方法である。
Aspect 4 of the present invention is
A method for evaluating damage behavior of a microstructure of a metallic material according to aspect 3, comprising obtaining the stress-strain curves of each of the two or more constituent phases prepared in advance as the material properties in the step 3. is.

本発明の態様5は、
前記ステップ7で、予め定められた手順に従って損傷クライテリアを決定し、決定された前記損傷クライテリアを用いてミクロ組織の損傷有無を判断する、態様1から4のいずれか一項に記載の金属材料のミクロ組織の損傷挙動評価方法である。
Aspect 5 of the present invention is
5. The metal material according to any one of aspects 1 to 4, wherein in step 7, a damage criterion is determined according to a predetermined procedure, and the presence or absence of microstructure damage is determined using the determined damage criterion. It is a damage behavior evaluation method of microstructure.

本発明の態様6は、
前記損傷クライテリアの決定を、
前記金属材料からなる試験片を用意し、該試験片に、連続的にひずみを付与しながら、該試験片の同一視野のミクロ組織画像を複数取得するステップ11と、
前記複数のミクロ組織画像から、マイクロボイド発生位置と、マイクロボイド発生時の塑性ひずみ量を決定するステップ12と、
前記ミクロ組織画像を複数の要素に分割した、前記ミクロ組織画像を模擬したFEM解析モデルを作成するステップ13と、
前記FEM解析モデルのパラメータとして前記金属材料の材料特性を取得するステップ14と、
境界条件として、前記FEM解析モデルの各要素に付与する、ひずみ量を設定するステップ15と、
前記境界条件でFEM解析を実行し、前記マイクロボイド発生時の各塑性ひずみ量における前記各要素の相当塑性ひずみと応力三軸度を出力するステップ16と、
出力された前記各要素の相当塑性ひずみと応力三軸度のデータの中から、前記マイクロボイド発生位置における相当塑性ひずみと応力三軸度のデータを抽出するステップ17と、
抽出された各位置の相当塑性ひずみと応力三軸度の分布を求めるステップ18と、
前記分布において、マイクロボイドが発生する相当塑性ひずみと応力三軸度の関係から、損傷クライテリアを決定するステップ19と
を含む工程で行う、態様5に記載の金属材料のミクロ組織の損傷挙動評価方法である。
Aspect 6 of the present invention is
determining the damage criteria,
A step 11 of preparing a test piece made of the metal material and acquiring a plurality of microstructure images of the same field of view of the test piece while continuously applying strain to the test piece;
Step 12 of determining the microvoid generation position and the amount of plastic strain at the time of microvoid generation from the plurality of microstructure images;
step 13 of creating an FEM analysis model simulating the microstructure image by dividing the microstructure image into a plurality of elements;
a step 14 of obtaining material properties of the metallic material as parameters of the FEM analysis model;
Step 15 of setting a strain amount to be applied to each element of the FEM analysis model as a boundary condition;
A step 16 of performing FEM analysis under the boundary conditions and outputting the equivalent plastic strain and stress triaxiality of each element at each plastic strain amount at the time of microvoid generation;
a step 17 of extracting the data of the equivalent plastic strain and the stress triaxiality at the microvoid generation position from the output data of the equivalent plastic strain and the stress triaxiality of each element;
a step 18 of obtaining distributions of equivalent plastic strain and stress triaxiality at each extracted position;
A method for evaluating the damage behavior of a microstructure of a metallic material according to aspect 5, which is performed in the step 19 of determining the damage criteria from the relationship between the equivalent plastic strain at which microvoids are generated and the stress triaxiality in the distribution. is.

本発明の態様7は、
前記金属材料は、2種類以上の構成相を有し、
前記損傷クライテリアの決定でのステップ17~19を構成相ごとに実施する、請求項6に記載の金属材料のミクロ組織の損傷挙動評価方法である。
Aspect 7 of the present invention is
The metal material has two or more constituent phases,
7. A method for evaluating damage behavior of a microstructure of a metallic material according to claim 6, wherein steps 17 to 19 in determining the damage criteria are performed for each constituent phase.

本発明によれば、複雑な数値解析モデルを使用することなく、かつ評価する金属材料のミクロ組織を対象に、実態に即した金属材料のミクロ組織損傷挙動の評価方法を提供することができる。 According to the present invention, it is possible to provide a method for evaluating the microstructure damage behavior of a metal material that is in line with the actual situation, without using a complicated numerical analysis model, and targeting the microstructure of the metal material to be evaluated.

本実施形態に係る評価システムのブロック図である。1 is a block diagram of an evaluation system according to this embodiment; FIG. 本実施形態に係るミクロ組織の損傷挙動を評価する方法をステップ順に示すフロー図である。FIG. 2 is a flow chart showing step by step a method for evaluating damage behavior of a microstructure according to the present embodiment. 本実施形態に係るミクロ組織を模擬した解析モデルを作成する手順の一例を示す図である。It is a figure which shows an example of the procedure of creating the analysis model which simulated the micro structure which concerns on this embodiment. 本実施形態に係るミクロ組織の損傷クライテリアを決定する方法をステップ順に示すフロー図である。FIG. 2 is a flow diagram showing the steps of a method for determining microstructural damage criteria according to the present embodiment. 実施例で用いた材料の顕微鏡写真である。1 is a micrograph of materials used in Examples. 実施例における相当塑性ひずみと応力三軸度のグラフである。4 is a graph of equivalent plastic strain and stress triaxiality in Examples. 実施例における相当塑性ひずみと応力三軸度の別のグラフである。It is another graph of equivalent plastic strain and stress triaxiality in the example. 実施例における相当塑性ひずみと応力三軸度の別のグラフである。It is another graph of equivalent plastic strain and stress triaxiality in the example. 実施例における相当塑性ひずみと応力三軸度の別のグラフである。It is another graph of equivalent plastic strain and stress triaxiality in the example. 実施例で用いた引張試験片の形状を示す図である。It is a figure which shows the shape of the tensile test piece used in the Example. 実施例で用いた準連続ミクロ組織観察-引張試験装置の側面図である。It is a side view of a quasi-continuous microstructure observation-tensile test apparatus used in Examples. 実施例で用いた引張試験片における観察位置を示す図である。It is a figure which shows the observation position in the tensile test piece used in the Example. 実施例において引張試験中の顕微鏡観察時期を示す図である。It is a figure which shows the microscope observation time during a tensile test in an Example. 実施例におけるミクロ観察結果の一部を示す顕微鏡写真である。It is a microscope photograph which shows a part of the microscopic observation result in an Example. 実施例におけるミクロ観察結果の別の一部を示す顕微鏡写真である。It is a microscope photograph which shows another part of the microscopic observation result in an Example. 実施例におけるミクロ観察結果の別の一部を示す顕微鏡写真である。It is a microscope photograph which shows another part of the microscopic observation result in an Example. 実施例におけるミクロ観察結果の別の一部を示す顕微鏡写真である。It is a microscope photograph which shows another part of the microscopic observation result in an Example. 実施例におけるミクロ観察結果の別の一部を示す顕微鏡写真である。It is a microscope photograph which shows another part of the microscopic observation result in an Example. 実施例におけるミクロ観察結果の別の一部を示す顕微鏡写真である。It is a microscope photograph which shows another part of the microscopic observation result in an Example. 実施例におけるミクロ観察結果の別の一部を示す顕微鏡写真である。It is a microscope photograph which shows another part of the microscopic observation result in an Example. 実施例におけるフェライト相の相当塑性ひずみと応力三軸度の関係と損傷クライテリアを示すグラフである。4 is a graph showing the relationship between the equivalent plastic strain of the ferrite phase and the stress triaxiality and the damage criteria in Examples. 実施例におけるマルテンサイト相の相当塑性ひずみと応力三軸度の関係と損傷クライテリアを示すグラフである。4 is a graph showing the relationship between the equivalent plastic strain of the martensite phase and the stress triaxiality and the damage criteria in Examples.

本発明者らは、従来技術の様に複雑な計算モデルを使用しなくとも、評価するミクロ組織を対象に、実態に即した損傷挙動の評価を行うことのできる方法を実現すべく、鋭意研究を行った。その結果、金属組織の数値解析で得られる相当塑性ひずみと応力三軸度を、ミクロ組織の損傷挙動の評価に適用すればよいこと、具体的には、ミクロ組織画像から作成した数値解析モデルにおいて、金属材料の破断までに生じる応力またはひずみまたはそれらの履歴を境界条件として設定して構造計算を行い、得られた相当塑性ひずみ(εeq)と応力三軸度(η)の分布により、ミクロ組織中の、マイクロボイドの発生位置、発生するタイミング、発生量等のマイクロボイド発生挙動に適用すればよいことを見出した。 The present inventors have conducted intensive research to realize a method that can evaluate the damage behavior of the microstructure to be evaluated in accordance with the actual situation without using a complicated calculation model as in the conventional technology. did As a result, the equivalent plastic strain and stress triaxiality obtained by numerical analysis of the metal structure should be applied to the evaluation of the damage behavior of the microstructure. , a structural calculation is performed by setting the stress or strain that occurs until the fracture of the metal material or their history as boundary conditions, and the distribution of the equivalent plastic strain (εeq P ) and stress triaxiality (η) obtained from the microscopic It was found that the method can be applied to microvoid generation behavior such as microvoid generation position, generation timing, generation amount, etc. in the tissue.

本実施形態に係る金属材料のミクロ組織損傷挙動を評価する方法を実行する一例である、評価システム1について説明する。図1は、本実施形態に係る評価システム1のブロック図である。評価システム1は、コンピュータシステムであって、試験部10と、制御部20と、記憶部30と、表示部40と、入力部50とを備える。少なくとも制御部20および記憶部30により、コンピュータが構成される。 An evaluation system 1, which is an example of executing a method for evaluating microstructure damage behavior of a metal material according to the present embodiment, will be described. FIG. 1 is a block diagram of an evaluation system 1 according to this embodiment. The evaluation system 1 is a computer system and includes a testing section 10 , a control section 20 , a storage section 30 , a display section 40 and an input section 50 . At least the control unit 20 and the storage unit 30 constitute a computer.

試験部10は、応力/ひずみ付与部11と、撮影部12とを有する。応力/ひずみ付与部11は、金属材料を保持し、所定の応力またはひずみを付与することができる実験装置である(詳細は後述する図11参照)。応力/ひずみ付与部11は、金属片の伸び計測用の撮影装置(例えば、CCD、CMOS等の撮像素子を備えたビデオカメラ)を備え、応力またはひずみの付与中に金属材料の変形を測定することができる。金属材料の詳細は、後述する。応力/ひずみ付与部11は、金属材料全体が画角に十分収まるように撮影を開始し、金属材料の変形を時間的に連続して撮影する。応力/ひずみ付与部11は、金属材料の左右方向または上下方向の端部間の距離を測定することで、変形を測定できる。また、後述するように金属材料に目印となるマーク(例えばシール)を取り付け、当該マーク間の距離を測定することで、変形を測定してもよい。 The testing section 10 has a stress/strain imparting section 11 and an imaging section 12 . The stress/strain applying unit 11 is an experimental apparatus capable of holding a metal material and applying a predetermined stress or strain (see FIG. 11 described later for details). The stress/strain imparting unit 11 includes an imaging device for measuring elongation of a metal piece (for example, a video camera equipped with an imaging device such as a CCD, CMOS, etc.), and measures deformation of the metal material during the application of stress or strain. be able to. Details of the metal material will be described later. The stress/strain imparting unit 11 starts photographing so that the entire metal material is sufficiently within the angle of view, and continuously photographs the deformation of the metal material over time. The stress/strain imparting unit 11 can measure deformation by measuring the distance between the ends of the metal material in the left-right direction or the up-down direction. Further, as will be described later, the deformation may be measured by attaching a mark (for example, a sticker) to the metal material and measuring the distance between the marks.

撮影部12は、例えば走査型電子顕微鏡(SEM)等であり、応力若しくはひずみを付与される前、および応力/ひずみ付与部で応力若しくはひずみを付与された金属材料の少なくとも一部を部分的に拡大したミクロ組織画像を撮影することができる。撮影部12は、撮影したミクロ組織画像のデータを制御部20に送信する。後述するように制御部20は、受け取ったミクロ組織画像のデータを記憶部30に保存する。 The imaging unit 12 is, for example, a scanning electron microscope (SEM) or the like, and partially captures at least a portion of the metal material before stress or strain is applied and to which stress or strain is applied in the stress/strain applying unit. A magnified microstructure image can be taken. The photographing unit 12 transmits data of the photographed microstructure image to the control unit 20 . As will be described later, the control unit 20 stores the received microstructure image data in the storage unit 30 .

制御部20は、モデル作成処理21、パラメータ設定処理22、解析処理23および分布作成処理24を実行する。制御部20は、プログラムを実行することにより所定の機能を実現するCPUまたはMPUのような汎用プロセッサを含む。制御部20は、例えば記憶部30に格納された演算プログラム等を呼び出して実行することにより、評価システム1における各種の処理を実現する。制御部20は、記憶部30に格納された演算プログラム等を呼び出して実行することにより、モデル作成処理21、パラメータ設定処理22、解析処理23および分布作成処理24を実現できる。制御部20は、ハードウェアとソフトウェアの協働により所定の機能を実現するものに限定されず、所定の機能を実現する専用に設計されたハードウェア回路でもよい。すなわち、制御部20は、CPU、MPU以外にも、GPU、FPGA、DSP、ASIC等、種々のプロセッサで実現され得る。このような制御部20は、例えば、半導体集積回路である信号処理回路で構成され得る。 The control unit 20 executes model creation processing 21 , parameter setting processing 22 , analysis processing 23 and distribution creation processing 24 . The control unit 20 includes a general-purpose processor such as a CPU or MPU that implements predetermined functions by executing programs. The control unit 20 realizes various types of processing in the evaluation system 1 by calling and executing a calculation program or the like stored in the storage unit 30, for example. The control unit 20 can implement a model creation process 21 , a parameter setting process 22 , an analysis process 23 and a distribution creation process 24 by calling up and executing a calculation program or the like stored in the storage unit 30 . The control unit 20 is not limited to one that realizes a predetermined function through the cooperation of hardware and software, and may be a hardware circuit designed exclusively for realizing a predetermined function. That is, the control unit 20 can be realized by various processors such as GPU, FPGA, DSP, ASIC, etc., in addition to CPU and MPU. Such a control unit 20 can be configured by, for example, a signal processing circuit that is a semiconductor integrated circuit.

制御部20は、モデル作成処理21により、撮影部12で撮影された金属材料のミクロ組織画像から、FEM解析で使用するためのFEM解析モデルを作成する。FEM解析モデルを作成するためのプログラム、パラメータ等は、例えば、記憶部30に格納されており、必要に応じて制御部20が読み出し、モデル作成に使用する。後述するように、制御部20は、例えば、イメージベースモデリングによりFEM解析モデルを作成することができる。イメージベースモデリングを行う場合、制御部20は、ミクロ組織画像の少なくとも一部の領域に対して任意の画像処理フィルタを適用し、金属材料の当該領域のミクロ組織形状に基づいた形状データを形成する。そして、制御部20は、後述するパラメータ設定処理22により設定されたFEM解析モデルの各パラメータ(例えば、作成されるモデルの要素の種類、要素数、節点数など)に基づいて、上記領域を模擬したFEM解析モデルを作成できる。作成されたFEM解析モデルは、記憶部30へと出力され、格納される。 The control unit 20 creates an FEM analysis model for use in FEM analysis from the microstructure image of the metal material captured by the image capturing unit 12 by model creation processing 21 . Programs, parameters, and the like for creating the FEM analysis model are stored, for example, in the storage unit 30, and are read by the control unit 20 as necessary and used for model creation. As will be described later, the control unit 20 can create an FEM analysis model by image-based modeling, for example. When performing image-based modeling, the control unit 20 applies any image processing filter to at least a partial region of the microstructure image to form shape data based on the microstructure shape of the region of the metal material. . Then, the control unit 20 simulates the region based on each parameter of the FEM analysis model (for example, the type of element of the model to be created, the number of elements, the number of nodes, etc.) set by the parameter setting process 22 described later. It is possible to create a FEM analysis model with The created FEM analysis model is output to the storage unit 30 and stored.

制御部20は、パラメータ設定処理22により、各種のパラメータを設定する。当該パラメータには、例えば、応力/ひずみ付与部11が付与する応力またはひずみ、制御部20がモデル作成処理21により作成するFEM解析モデルの要素の種類、要素数および節点数が含まれてもよい。また、当該パラメータには、FEM解析によって解析したい金属材料の特性(金属材料のヤング率、降伏応力および耐力、加工硬化指数並びに結晶方位など)、解析処理23によりFEM解析する際の境界条件等が含まれてもよい。当該パラメータは、記憶部30へと出力され、格納される。そして、必要に応じて、各構成要素は記憶部30から当該パラメータを読み込み、使用する。 The control unit 20 sets various parameters by the parameter setting process 22 . The parameters may include, for example, the stress or strain applied by the stress/strain applying unit 11, the type of elements of the FEM analysis model created by the control unit 20 by the model creation process 21, the number of elements, and the number of nodes. . In addition, the parameters include the characteristics of the metal material to be analyzed by FEM analysis (Young's modulus, yield stress and yield strength, work hardening index and crystal orientation of the metal material, etc.), boundary conditions when performing FEM analysis by analysis processing 23, etc. may be included. The parameters are output to and stored in the storage unit 30 . Then, each component reads the parameter from the storage unit 30 and uses it as necessary.

制御部20は、解析処理23により、作成されたFEM解析モデルを用いて、金属材料に応力またはひずみが付与された際の形状変形をシミュレーションするためにFEM解析を行うことができる。解析に使用するFEM解析モデル、シミュレーションを行うためのプログラム、境界条件等は、例えば記憶部30に格納されており、必要に応じて制御部20が読み出し、解析に使用する。本実施形態に係る評価システム1の制御部20は、解析処理23により、FEM解析モデルに与えた応力またはひずみ別に、FEM解析モデルの各要素の積分点における応力三軸度および相当塑性ひずみを算出し、出力することができる。応力三軸度および相当塑性ひずみは、表示部40へと出力されて、表示部40で表示されてもよいし、記憶部30に出力され、格納されてもよいし、後述する出力部60へと出力され、別の制御装置へと伝達されてもよい。 The control unit 20 can use the created FEM analysis model by the analysis processing 23 to perform FEM analysis for simulating shape deformation when stress or strain is applied to the metal material. The FEM analysis model used for the analysis, the program for performing the simulation, the boundary conditions, etc. are stored, for example, in the storage unit 30, and are read by the control unit 20 as necessary and used for the analysis. The control unit 20 of the evaluation system 1 according to the present embodiment calculates the stress triaxiality and the equivalent plastic strain at the integration point of each element of the FEM analysis model for each stress or strain applied to the FEM analysis model by the analysis processing 23. can be output. The stress triaxiality and the equivalent plastic strain may be output to the display unit 40 and displayed on the display unit 40, may be output to the storage unit 30 and stored, or may be output to the output unit 60 described later. and may be transmitted to another control device.

制御部20は、例えばグラフ作成処理等の分布作成処理24により、付与した応力またはひずみ別に、解析処理23により算出し出力したFEM解析モデルの各積分点における応力三軸度および相当塑性ひずみの分布を求めることができる。具体的に例えば、それらを軸とするグラフにプロットすることができる。また、使用者が例えば入力部50を使って積分点を選択し、選択した積分点の応力三軸度および相当塑性ひずみをグラフにプロットすることもできる。グラフを作成するためのプログラム、パラメータ等は、例えば、記憶部30に格納されており、必要に応じて制御部20が読み出し、グラフ作成に使用する。 The control unit 20, for example, by distribution creation processing 24 such as graph creation processing, for each applied stress or strain, distribution of stress triaxiality and equivalent plastic strain at each integration point of the FEM analysis model calculated and output by analysis processing 23 can be asked for. Specifically, for example, they can be plotted on a graph with these as axes. Alternatively, the user can select an integration point using the input unit 50, for example, and plot the stress triaxiality and equivalent plastic strain of the selected integration point on the graph. Programs, parameters, and the like for creating the graph are stored, for example, in the storage unit 30, and are read by the control unit 20 as necessary and used for creating the graph.

記憶部30は、種々の情報を記録できる記録媒体である。記憶部30は、例えば、DRAMやSRAM、フラッシュメモリ等のメモリ、HDD、SSD、その他の記憶デバイスまたはそれらを適宜組み合わせて実現される。記憶部30は、撮影部12で撮影された画像、モデル作成処理21で作成されたFEM解析モデル等のデータを格納することができる。また、記憶部30は、パラメータ設定処理22で設定された、解析に使用する材料特性や境界条件等のパラメータ、および解析処理23で算出し出力された解析結果のデータ等を格納することができる。さらに、必要に応じて、記憶部30は、制御部20、表示部40等の各構成要素へと格納している情報を送信することができる。各構成要素は、記憶部30を介さずに各構成要素間で直接、情報を送受信してもよい。 The storage unit 30 is a recording medium capable of recording various information. The storage unit 30 is realized by, for example, a memory such as a DRAM, an SRAM, a flash memory, an HDD, an SSD, other storage devices, or an appropriate combination thereof. The storage unit 30 can store data such as the image captured by the imaging unit 12 and the FEM analysis model created by the model creating process 21 . In addition, the storage unit 30 can store parameters such as material properties and boundary conditions used for analysis set in the parameter setting process 22, data of analysis results calculated and output in the analysis process 23, and the like. . Furthermore, the storage unit 30 can transmit the stored information to each component such as the control unit 20 and the display unit 40 as necessary. Each component may transmit and receive information directly between the components without going through the storage unit 30 .

表示部40は、任意のデータを表示できるディスプレイであり、金属材料のミクロ組織画像、FEM解析モデル、解析に使用する材料特性や境界条件等のパラメータ、算出された解析結果、作成されたグラフ等、任意の情報を表示することができる。使用者は、表示部40に表示された情報を基に入力部50を操作して制御部20のパラメータ設定処理22によってパラメータを設定したり、撮影部12で撮影されたミクロ組織画像と、制御部20が解析処理23で算出した解析結果とを対比させたりすることができる。 The display unit 40 is a display that can display arbitrary data, such as microstructure images of metal materials, FEM analysis models, parameters such as material properties and boundary conditions used for analysis, calculated analysis results, created graphs, and the like. , can display arbitrary information. The user operates the input unit 50 based on the information displayed on the display unit 40 to set the parameters by the parameter setting processing 22 of the control unit 20, or the microstructure image captured by the imaging unit 12 and the control The analysis result calculated by the analysis processing 23 by the unit 20 can be compared.

入力部50は、例えば、キーボードやマウスなどのコンピュータの一般的な入力装置であり、制御部20のパラメータ設定処理22にパラメータ等を入力したり、表示部40に表示された制御部20の解析処理23による解析結果に対して位置を指示したりできる。また、表示部40にタッチパネルを用いた場合など、表示部40が、使用者が触れた位置を検知する機構のような入力部50に相当する機能を併せ持ってもよい。 The input unit 50 is, for example, a general computer input device such as a keyboard and a mouse, and inputs parameters and the like to the parameter setting process 22 of the control unit 20, and analyzes the control unit 20 displayed on the display unit 40. A position can be indicated with respect to the analysis result by the processing 23 . Further, when a touch panel is used for the display section 40, the display section 40 may also have a function corresponding to the input section 50, such as a mechanism for detecting a position touched by the user.

また、評価システム1は、他の装置へと取得した情報や算出した情報等を伝達することができる出力部60等を備えてもよい。出力部60は、例えば、USB端子またはイーサネット端子である。また、無線通信を行うための通信装置であってもよい。そのような通信装置として、例えば2.4GHz/5.2GHz/5.3GHz/5.6GHz等の周波数を利用して無線通信を行う、Wi-Fi(登録商標)規格に準拠した規格が採用され得る。 The evaluation system 1 may also include an output unit 60 or the like capable of transmitting acquired information, calculated information, and the like to other devices. The output unit 60 is, for example, a USB terminal or an Ethernet terminal. Moreover, it may be a communication device for performing wireless communication. As such a communication device, a standard conforming to the Wi-Fi (registered trademark) standard, which performs wireless communication using frequencies such as 2.4 GHz/5.2 GHz/5.3 GHz/5.6 GHz, is adopted. obtain.

なお、評価システム1は少なくとも制御部20と記憶部30とを有していればよい。例えば、制御部20は、試験部10から画像を受け取るための受信端子、及び、入力部50から使用者の指示を受け付ける受信端子を有し、表示部40および/または出力部60に信号を出力するための出力端子を有していればよい。このとき、試験部10、表示部40、入力部50及び出力部60の一部または全部は、評価システム1に包含されなくてもよい。 Note that the evaluation system 1 only needs to have at least the control unit 20 and the storage unit 30 . For example, the control section 20 has a receiving terminal for receiving an image from the testing section 10 and a receiving terminal for receiving a user's instruction from the input section 50, and outputs a signal to the display section 40 and/or the output section 60. It is only necessary to have an output terminal for At this time, some or all of the testing section 10, the display section 40, the input section 50, and the output section 60 may not be included in the evaluation system 1.

本実施形態に係る金属材料のミクロ組織損傷挙動を評価する方法について、図2のミクロ組織の損傷挙動評価フロー図を例に説明する。 A method for evaluating the microstructure damage behavior of a metal material according to the present embodiment will be described with reference to the microstructure damage behavior evaluation flowchart of FIG. 2 as an example.

[ステップ1:金属材料のミクロ組織画像の取得]
まず図2のステップ1(S1)の通り、金属材料のミクロ組織画像を取得する。
[Step 1: Acquisition of microstructure image of metal material]
First, as in step 1 (S1) of FIG. 2, a microstructure image of the metal material is acquired.

ミクロ組織画像は、例えば、走査型電子顕微鏡(SEM)、光学顕微鏡(OM)、電子線後方散乱回折(EBSD)、コンピュータ断層撮影(CT)、超音波探傷検査(UT)等で観察した画像を用いればよい。ミクロ組織を構成する構成相の種類によっては、エッチング等の前処理をあらかじめ行ってもよい。上記ミクロ組織画像のサイズ、倍率等は、評価領域、評価位置等に応じて適宜設定すればよい。 The microstructure image is, for example, an image observed with a scanning electron microscope (SEM), an optical microscope (OM), electron beam backscatter diffraction (EBSD), computed tomography (CT), ultrasonic inspection (UT), or the like. You can use it. Pretreatment such as etching may be performed in advance depending on the type of constituent phases that constitute the microstructure. The size, magnification, etc. of the microstructure image may be appropriately set according to the evaluation region, evaluation position, and the like.

対象とする金属材料の鋼種、強度レベル等は問わない。金属材料として、鉄鋼に限らず、アルミニウム系、チタン系、銅系、マグネシウム系の純金属または合金が挙げられる。金属材料が鉄鋼である場合、ミクロ組織の種類、組織形状、相分率は限定されない。ミクロ組織を構成する構成相として、例えば、フェライト、マルテンサイト、ベイナイト、パーライト、残留オーステナイト、MAのうちの1以上を含み、その他セメンタイト等の析出物を含みうる。前記組織形状は、結晶の形状をいい、例えば枝状、粒状、針状等が挙げられる。前記鉄鋼の成分組成も限定されない。前記金属材料の製造方法も限定されず、所望のミクロ組織が得られるように、例えば熱間圧延条件、熱処理条件等の製造条件が適宜設定されうる。複数の金属材料を評価、例えば複数の金属材料間で対比する場合は、例えばミクロ組織の構成相の種類、組織形状、相分率のうちの1以上が異なる複数の金属材料を用意することが挙げられる。簡便に説明するため、以下に記載する評価方法において、金属材料は、軟質相としてフェライトおよび硬質相としてマルテンサイトから構成される二相鋼を、金属材料のミクロ組織損傷挙動を評価する対象としている。 The steel type, strength level, etc. of the target metal material do not matter. The metal material is not limited to steel, and includes aluminum, titanium, copper, and magnesium pure metals or alloys. When the metal material is steel, the type of microstructure, structure shape, and phase fraction are not limited. Constituent phases constituting the microstructure include, for example, one or more of ferrite, martensite, bainite, pearlite, retained austenite, and MA, and may also include precipitates such as cementite. The structure shape refers to the shape of crystals, and examples thereof include branch-like, granular, needle-like, and the like. The chemical composition of the steel is also not limited. The method for producing the metal material is not limited, either, and production conditions such as hot rolling conditions and heat treatment conditions can be appropriately set so as to obtain a desired microstructure. When evaluating a plurality of metal materials, for example, when comparing a plurality of metal materials, it is possible to prepare a plurality of metal materials that differ in one or more of, for example, the type of constituent phase of the microstructure, the shape of the structure, and the phase fraction. mentioned. For the sake of simplicity, in the evaluation method described below, the metallic material is a duplex steel composed of ferrite as the soft phase and martensite as the hard phase, and the microstructural damage behavior of the metallic material is evaluated. .

[ステップ2:FEM解析モデルの作成]
図2のステップ2(S2)の通り、制御部20は、モデル作成処理21にて、FEM解析に使用するFEM解析モデルとしてメッシュモデルを作成する。当該解析モデルは、例えば、金属材料の取得したミクロ組織画像を用いてイメージベースモデリングにより作成できる。モデル作成方法の一例を図3に記載する。図3Aは取得されたミクロ組織画像を示す。制御部20は、このようなミクロ組織画像の画素の濃度または輝度に対して所定の閾値を設定し、任意の画素位置の濃度等が当該閾値を超えるか否かによって当該画素位置の値を設定する二値化処理を行うことで二値画像を作成する。図3Bは、図3Aのミクロ組織画像から作成された二値画像である。これにより、マルテンサイト相を抽出し、それ以外の相と分離することができる。
[Step 2: Creation of FEM analysis model]
As shown in step 2 (S2) of FIG. 2, the control unit 20 creates a mesh model as an FEM analysis model to be used for FEM analysis in the model creation process 21. FIG. The analysis model can be created, for example, by image-based modeling using acquired microstructure images of the metal material. An example of the model creation method is described in FIG. FIG. 3A shows the acquired microstructure image. The control unit 20 sets a predetermined threshold value for the density or brightness of pixels in such a microstructure image, and sets a value for a given pixel position depending on whether the density or the like at an arbitrary pixel position exceeds the threshold value. A binary image is created by performing binarization processing. FIG. 3B is a binary image created from the microstructure image of FIG. 3A. This allows the martensite phase to be extracted and separated from the other phases.

次に、制御部20は、二値画像からマルテンサイト相に対応する領域の輪郭を取得する。当該輪郭は、例えばエッジを検出する画像処理フィルタ(例えば、ソーベルフィルタ、ラプラシアンフィルタ等)を用いることで取得することができる。図3Cは、二値画像から取得された、マルテンサイト相(硬質相)に対応する領域の輪郭を表す。 Next, the control unit 20 acquires the contour of the region corresponding to the martensite phase from the binary image. The contour can be acquired, for example, by using an image processing filter (eg, Sobel filter, Laplacian filter, etc.) that detects edges. FIG. 3C represents the contour of the region corresponding to the martensite phase (hard phase) obtained from the binary image.

次に、制御部20は、上記の輪郭から、マルテンサイト相に対応する領域の形状データと、フェライト相に対応する領域の形状データとを作成する。図3Dは、マルテンサイト相に対応する領域Mと、フェライト相に対応する領域Fを作成した形状データである。その後、制御部20は、当該形状データを所定の条件に基づいて設定されるメッシュに分割した2次元のメッシュモデルを作成する。図3Eは、このようにして作成されたメッシュ領域の一例である。本実施形態に係るメッシュモデルでは、四辺形要素のメッシュに分割している。図3Eには、右上の一部の領域を部分的に拡大し、メッシュの形状を例示した。四辺形要素は長方形である必要はない。また四辺形であることも必須ではなく、三角形以上の多角形であればよい。 Next, the control unit 20 creates shape data of the region corresponding to the martensite phase and shape data of the region corresponding to the ferrite phase from the contours. FIG. 3D shows shape data for creating a region M corresponding to the martensite phase and a region F corresponding to the ferrite phase. After that, the control unit 20 creates a two-dimensional mesh model by dividing the shape data into meshes set based on predetermined conditions. FIG. 3E is an example of a mesh region created in this way. The mesh model according to this embodiment is divided into meshes of quadrilateral elements. FIG. 3E shows an example of the shape of the mesh by partially enlarging a portion of the upper right region. A quadrilateral element need not be rectangular. Moreover, it is not essential to be a quadrilateral, and it is sufficient if the polygon is a triangle or more.

本実施形態におけるFEM解析モデルとして、上記のメッシュモデルは、4節点アイソパラメトリック要素を設定して、総要素数が40,279、総節点数が40,507の2Dモデルで作成する。そして、当該2Dモデルを画像の奥行き方向へ1層分拡張し、8節点アイソパラメトリック要素により総要素数が40,279、総節点数が81,014の3Dモデルを作成し、当該3Dモデルを用いてFEM解析を行う。モデルの条件は上記に限定されるものではなく、例えば、要素の種類に関して、構成節点数が異なる要素を使用してもよい。また、総要素数や総節点数が異なる3Dモデルを作成して、FEM解析を行ってもよい。また、3Dモデルではなく、2Dモデルを作成してFEM解析モデルとして使用してもよい。 As the FEM analysis model in this embodiment, the above mesh model is created as a 2D model with a total number of elements of 40,279 and a total number of nodes of 40,507 by setting 4-node isoparametric elements. Then, the 2D model is expanded by one layer in the depth direction of the image, and a 3D model with a total number of elements of 40,279 and a total number of nodes of 81,014 is created by 8-node isoparametric elements, and the 3D model is used. FEM analysis is performed by The conditions of the model are not limited to those described above, and for example, elements with different numbers of constituent nodes may be used with respect to the types of elements. In addition, FEM analysis may be performed by creating 3D models with different total number of elements and total number of nodes. Also, instead of the 3D model, a 2D model may be created and used as the FEM analysis model.

[ステップ3:材料特性の取得]
図2のステップ3(S3)の通り、制御部20は、パラメータ設定処理22により、FEM解析を行うために二相鋼を構成する各相に関する各材料の特性を取得し、パラメータとして設定する。ここで、取得する材料特性とは、例えば、フェライト相とマルテンサイト相それぞれのヤング率、降伏応力および耐力、加工硬化指数並びに結晶方位などである。また、制御部20は、金属材料を構成する各材料に関する応力-ひずみ曲線(本実施形態においては、フェライトに関する応力-ひずみ曲線とマルテンサイトに関する応力-ひずみ曲線)をそれぞれ材料特性として取得する。これにより、制御部20が解析処理23にて、後述するFEM解析を行った際、FEM解析モデルに与える応力またはひずみに対して、当該モデルを構成する各相の形状および特性に基づいて、各積分点に与えられる応力またはひずみが算出される。金属材料を構成する各材料に関する応力-ひずみ曲線は、例えば、当該材料で構成されている試験片に対して引張試験を行うことで、取得することができる。
[Step 3: Acquisition of material properties]
As shown in step 3 (S3) of FIG. 2, the control unit 20 acquires the properties of each material related to each phase that constitutes the duplex steel and sets them as parameters by the parameter setting process 22 in order to perform the FEM analysis. Here, the material properties to be obtained are, for example, Young's modulus, yield stress and yield strength, work hardening index and crystal orientation of each of the ferrite phase and the martensite phase. In addition, the control unit 20 acquires the stress-strain curve of each material constituting the metal material (in this embodiment, the stress-strain curve of ferrite and the stress-strain curve of martensite) as material characteristics. As a result, when the control unit 20 performs the FEM analysis described later in the analysis processing 23, each A stress or strain applied to the integration point is calculated. A stress-strain curve for each material that constitutes a metal material can be obtained, for example, by performing a tensile test on a test piece that is composed of the material.

[ステップ4:境界条件の設定]
図2のステップ4(S4)の通り、制御部20は、パラメータ設定処理22により、FEM解析モデルに与える境界条件を設定する。当該境界条件は、金属材料を構成する構造体が破壊されるまでに実際に材料に生じると考えられる応力またはひずみの履歴を与える。当該応力またはひずみの履歴は、構造体の構造解析によりあらかじめ推定することができる。また、実際に構造物から、例えばひずみゲージを用いることで計測してもよいし、任意の条件を設定してもよい。境界条件は、応力またはひずみの履歴に限定されず、例えば、応力またはひずみを与えてもよい。この場合、制御部20はパラメータ設定処理22により、当該境界条件に基づいて所定の増加幅で、FEM解析を行うためのパラメータを設定する。
[Step 4: Setting Boundary Conditions]
As in step 4 (S4) of FIG. 2, the control unit 20 sets boundary conditions to be given to the FEM analysis model by the parameter setting process 22. FIG. The boundary conditions provide the stress or strain history that is actually considered to occur in the material until the structure that constitutes the metal material fails. The history of stress or strain can be estimated in advance by structural analysis of the structure. Moreover, it may be measured by actually using, for example, a strain gauge from the structure, or arbitrary conditions may be set. Boundary conditions are not limited to stress or strain histories and may, for example, impose stress or strain. In this case, the control unit 20 sets the parameters for performing the FEM analysis by the parameter setting process 22 with a predetermined increase width based on the boundary conditions.

[ステップ5:FEM解析の実行(相当塑性ひずみと応力三軸度の出力)]
図2のステップ5(S5)の通り、制御部20は、解析処理23により、作成したFEM解析モデルを用いてFEM解析を実行する。FEM解析では、作成したFEM解析モデルに対して境界条件として設定した応力またはひずみの履歴を与えることで、金属材料に応力またはひずみが加えられた際にどのように変形するかシミュレーションにより解析することができる。解析により、FEM解析モデルの各要素に与えた応力またはひずみ別に、FEM解析モデルの各要素における「応力三軸度」および「相当塑性ひずみ」を算出し、出力することができる。FEM解析モデルの各要素には、解析時の積分を行うための積分点が予め設定されている。「応力三軸度」および「相当塑性ひずみ」は、積分点に与えられたに与えた応力またはひずみ別に算出される。なお、FEM解析モデルの各要素の積分点以外に、例えば積分点について算出された解から計算される要素解または節点解を用いて「応力三軸度」および「相当塑性ひずみ」を算出してもよい。
[Step 5: Execution of FEM analysis (output of equivalent plastic strain and stress triaxiality)]
As shown in step 5 (S5) of FIG. 2, the control unit 20 executes the FEM analysis using the created FEM analysis model by the analysis processing 23. FIG. In FEM analysis, by giving the history of stress or strain set as boundary conditions to the created FEM analysis model, it is possible to analyze by simulation how the metal material deforms when stress or strain is applied. can be done. Through the analysis, the "stress triaxiality" and "equivalent plastic strain" in each element of the FEM analysis model can be calculated and output for each stress or strain applied to each element of the FEM analysis model. Each element of the FEM analysis model has preset integration points for performing integration during analysis. "Stress triaxiality" and "equivalent plastic strain" are calculated for each stress or strain applied to the integration points. In addition to the integration point of each element of the FEM analysis model, for example, the element solution or nodal solution calculated from the solution calculated for the integration point is used to calculate the "stress triaxiality" and "equivalent plastic strain". good too.

[ステップ6:相当塑性ひずみと応力三軸度の分布作成]
図2のステップ6(S6)の通り、制御部20は、例えばグラフ作成処理等の分布作成処理24により、相当塑性ひずみと応力三軸度の分布を求めることができる。具体的に例えば、相当塑性ひずみと応力三軸度のグラフを作成する。詳細には、前記境界条件として設定した応力またはひずみ別に、出力された前記各位置の相当塑性ひずみと応力三軸度を、これらを軸とするグラフにプロットする。また、制御部20は、算出された相当塑性ひずみと応力三軸度を所定の時間刻みごとに(例えば、FEM解析において1回の解析あたりに進める時間間隔ごとに)当該グラフにプロットしてもよい。前記グラフは、例えば縦軸を相当塑性ひずみ、横軸を応力三軸度とすることができる。後記の実施例では、全ての積分点の相当塑性ひずみと応力三軸度をプロットしているが、これに限られず、前記位置として、前記モデルにおける一部の積分点の、相当塑性ひずみと応力三軸度をプロットしてもよい。
[Step 6: Creation of distribution of equivalent plastic strain and stress triaxiality]
As shown in step 6 (S6) in FIG. 2, the control unit 20 can obtain the distribution of the equivalent plastic strain and the stress triaxiality by the distribution creating process 24 such as the graph creating process. Specifically, for example, a graph of equivalent plastic strain and stress triaxiality is created. Specifically, for each stress or strain set as the boundary condition, the output equivalent plastic strain and stress triaxiality at each position are plotted on a graph with these as axes. In addition, the control unit 20 plots the calculated equivalent plastic strain and stress triaxiality on the graph for each predetermined time interval (for example, for each time interval advanced per analysis in FEM analysis). good. In the graph, for example, the vertical axis can be the equivalent plastic strain, and the horizontal axis can be the stress triaxiality. In the examples below, the equivalent plastic strain and stress triaxiality of all the integration points are plotted, but not limited to this, the equivalent plastic strain and stress of some integration points in the model as the position Triaxiality may be plotted.

[ステップ7:ミクロ組織の損傷挙動の評価]
図2のステップ7(S7)の通り、前記相当塑性ひずみと応力三軸度の分布から、ミクロ組織の損傷挙動を評価する。評価方法の一例として、前記境界条件として設定した応力またはひずみ別に、前記プロット点の分布形態を観察することが挙げられる。または、前記境界条件として設定した応力またはひずみの異なる、プロット点の2以上の分布を対比することによって、変形過程すなわち塑性ひずみ(εmicro)量の変化に伴うプロット点の推移を観察し、マイクロボイド発生挙動の評価に用いることができる。本実施形態によれば、従来技術の様に複雑な連続体損傷モデルを用いる必要がなく、金属材料のマイクロボイド発生挙動を簡便に評価することが可能となる。
[Step 7: Evaluation of microstructural damage behavior]
As shown in step 7 (S7) of FIG. 2, the damage behavior of the microstructure is evaluated from the distribution of the equivalent plastic strain and the stress triaxiality. An example of the evaluation method is to observe the distribution form of the plot points for each stress or strain set as the boundary condition. Alternatively, by comparing two or more distributions of plot points with different stresses or strains set as the boundary conditions, observe the transition of the plot points accompanying the deformation process, that is, the change in the amount of plastic strain (εmicro), and microvoids It can be used to evaluate developmental behavior. According to this embodiment, it is not necessary to use a complicated continuum damage model as in the prior art, and it is possible to easily evaluate the microvoid generation behavior of a metal material.

本実施形態に係る評価方法は、前述の通り、対象とする金属材料を特に限定しない。前記金属材料は、構成相が1種類であってもよいし、2種類以上であってもよい。前記金属材料の構成相が2種類以上である場合、前記ステップ5~7を構成相ごとに実施することができる。この場合、前記ステップ7で、構成相別に示した相当塑性ひずみと応力三軸度のプロット点の分布形態を対比して、マイクロボイド発生のしやすさを評価することが挙げられる。 As described above, the evaluation method according to this embodiment does not particularly limit the target metal material. The metal material may have one type of constituent phase, or may have two or more types of constituent phases. When the metal material has two or more constituent phases, steps 5 to 7 can be performed for each constituent phase. In this case, in step 7, the susceptibility to microvoid generation may be evaluated by comparing the plot point distributions of the equivalent plastic strain and the stress triaxiality shown for each constituent phase.

前記ステップ7で、別途用意した損傷クライテリアを用い、ミクロ組織の損傷有無を判断してもよい。損傷クライテリアとして、下記の通り、実験的に決定した損傷クライテリアを用いることが好ましい。該損傷クライテリアは、まず、ミクロ組織の損傷であるマイクロボイドの発生位置とマイクロボイド発生時の塑性ひずみ量を簡便な実験により求め、これらの情報と、金属材料のFEM解析で得られる、マイクロボイド発生位置の相当塑性ひずみおよび応力三軸度とを、関連付けることで得られる。該損傷クライテリアを用いれば、実態に即したミクロ組織の損傷挙動の評価を行うことができる。なお、本実施例では、該損傷クライテリアを導出するために複数のプロット点をグラフ上に表示したが、複数のプロット点を表示することは必須ではない。複数のプロット点の表示に代えて、導出された該損傷クライテリアを表示してもよい。 In the above step 7, damage criteria prepared separately may be used to determine whether or not the microstructure is damaged. As the damage criteria, experimentally determined damage criteria are preferably used as described below. The damage criteria are first obtained by a simple experiment to determine the location of microvoid generation, which is damage to the microstructure, and the amount of plastic strain at the time of microvoid generation. It is obtained by relating the equivalent plastic strain and the stress triaxiality at the point of origin. By using the damage criteria, it is possible to evaluate the damage behavior of microstructures in line with actual conditions. In this embodiment, a plurality of plot points are displayed on the graph to derive the damage criteria, but it is not essential to display a plurality of plot points. Instead of displaying a plurality of plot points, the derived damage criteria may be displayed.

以下、損傷クライテリアを実験的に決定する方法について、図4のミクロ組織の損傷クライテリア決定フロー図を例に説明する。 The method for experimentally determining the damage criteria will be described below with reference to the microstructure damage criteria determination flow diagram of FIG. 4 as an example.

[ステップ11:塑性変形中の複数のミクロ組織画像の取得]
図4のステップ11(S11)の通り、撮影部12は、塑性変形中の複数のミクロ組織画像を取得する。詳細には、金属材料からなる試験片を用意し、該試験片に、連続的にひずみを付与しながら、該試験片の同一視野のミクロ組織画像を複数取得する。
前記ひずみを付与する方法は特に限定されず、連続的にひずみを付与できればよい。前記方法として、例えば引張試験が挙げられ、例えば一軸方向の引張試験を行うことが挙げられる。上記引張試験を行う場合、後述する実施例に記載の通り、所定の塑性ひずみ量になるまで試験片を引っ張り、引張試験機から試験片を取り外して顕微鏡観察し、再度、試験片を引張試験機に設置して引っ張ることを繰り返し行う、準連続ミクロ組織観察-引張試験が挙げられる。しかしこれに限定されず、種々の塑性ひずみ量を付与時の複数枚のミクロ組織画像を取得できればよい。例えば非特許文献3に示されたようなミクロ組織観察と引張試験の方法等を採用してもよい。
[Step 11: Acquisition of multiple microstructure images during plastic deformation]
As in step 11 (S11) of FIG. 4, the imaging unit 12 acquires a plurality of microstructure images during plastic deformation. Specifically, a test piece made of a metal material is prepared, and a plurality of microstructure images of the same field of view of the test piece are acquired while continuously applying strain to the test piece.
The method of applying the strain is not particularly limited as long as the strain can be applied continuously. The method includes, for example, a tensile test, for example, performing a uniaxial tensile test. When performing the above tensile test, as described in the examples described later, the test piece is pulled until a predetermined amount of plastic strain is reached, the test piece is removed from the tensile tester and observed under a microscope, and the test piece is placed in the tensile tester again. A quasi-continuous microstructure observation-tensile test, in which the specimen is placed in a slab and repeatedly pulled, can be mentioned. However, the present invention is not limited to this, as long as it is possible to acquire a plurality of microstructure images when applying various amounts of plastic strain. For example, the method of microstructure observation and tensile test as shown in Non-Patent Document 3 may be adopted.

前記引張試験に用いる引張試験片は、前記ミクロ組織画像を取得するための面が、鏡面研磨、さらにはミクロ組織の種類により必要に応じて、エッチング等の前処理があらかじめ施されていてもよい。前記ミクロ組織画像は、前記ステップ1と同様に、走査型電子顕微鏡(SEM)等で観察した画像を用いることができる。上記ミクロ組織画像のサイズ、倍率等は、評価領域、評価位置等に応じて適宜設定すればよい。 In the tensile test piece used in the tensile test, the surface for obtaining the microstructure image may be mirror-polished and, if necessary, subjected to pretreatment such as etching depending on the type of microstructure. . For the microstructure image, an image observed with a scanning electron microscope (SEM) or the like can be used, as in step 1 above. The size, magnification, etc. of the microstructure image may be appropriately set according to the evaluation region, evaluation position, and the like.

前記引張試験による引張試験片の公称ひずみ量は、後述する実施例に記載の通り、試験片に伸びを計測するためのマーク(シール)を付した状態で、試験片の変形を撮影し、前記マークの位置から求めることができる。また、試験片の変形前後のミクロ組織写真に表示された金属材料の組織形状等を対比することで、同一の領域を特定し、当該領域の変形量からの推定により、当該写真が撮影された時点でのミクロ組織の局所的な領域の塑性ひずみ量(εmicro)を求めることができる。なお公称ひずみ量および塑性ひずみ量は、例えば、変形前の試験片の形状を100としたときの変位の程度(%)で表される。一例として、本発明者は、公称ひずみεnominal=13.3%、及び公称ひずみεnominal=18.2%の時点において、塑性ひずみεmicro=17.1%、及び塑性ひずみεmicro=35.5%のミクロ組織画像を取得した。 The nominal strain amount of the tensile test piece in the tensile test is obtained by photographing the deformation of the test piece with a mark (sticker) attached to the test piece for measuring elongation, as described in the examples described later. It can be obtained from the position of the mark. In addition, by comparing the structure shape of the metal material displayed in the microstructure photographs before and after deformation of the test piece, the same region was identified, and the photograph was taken by estimating the deformation amount of the region. The amount of plastic strain (ε micro ) in a localized region of the microstructure at the point in time can be determined. Note that the nominal strain amount and the plastic strain amount are expressed, for example, by degrees of displacement (%) when the shape of the test piece before deformation is taken as 100. As an example, we find that at nominal strain ε nominal =13.3% and nominal strain ε nominal =18.2%, plastic strain ε micro =17.1% and plastic strain ε micro =35.2%. 5% microstructural images were acquired.

[ステップ12:ミクロ組織画像からのマイクロボイド発生位置とマイクロボイド発生時の塑性ひずみの把握]
図4のステップ12(S12)の通り、前記複数のミクロ組織画像から、マイクロボイド発生位置と、マイクロボイド発生時の塑性ひずみ量(εmicro)を把握する。マイクロボイドが発生した位置の確認は、前記ミクロ組織画像において目視でまたは画像解析で行うことができる。マイクロボイド発生位置は、金属材料を構成する構成相が2種類以上である場合、いずれの構成相で発生したかも把握する。例えば、ミクロ組織画像に、マイクロボイド発生位置と、マイクロボイド発生時の金属材料に付与した塑性ひずみ量(εmicro)を記録することが挙げられる。
[Step 12: Grasping the microvoid generation position from the microstructure image and the plastic strain at the time of microvoid generation]
As shown in step 12 (S12) of FIG. 4, microvoid generation positions and plastic strain amounts (ε micro ) at the time of microvoid generation are grasped from the plurality of microstructure images. The positions where microvoids are generated can be confirmed visually or by image analysis in the microstructure image. In the case where there are two or more types of constituent phases constituting the metal material, the location of occurrence of microvoids is determined in which constituent phase. For example, a microstructure image may be used to record the location of microvoid generation and the amount of plastic strain (ε micro ) imparted to the metal material when microvoids were generated.

[ステップ13:FEM解析モデルの作成]
図4のステップ13(S13)の通り、制御部20は、モデル作成処理21により、FEM解析に使用するFEM解析モデルとして、メッシュモデルを作成する。FEM解析モデルは、上述のステップ2と同様の方法で作成することができる。
[Step 13: Creation of FEM analysis model]
As shown in step 13 (S13) of FIG. 4, the control unit 20 creates a mesh model as an FEM analysis model to be used for FEM analysis by the model creating process 21. FIG. The FEM analysis model can be created in a manner similar to step 2 above.

[ステップ14:材料特性の取得]
図4のステップ14(S14)の通り、制御部20は、パラメータ設定処理22により、ステップ11で試験した金属材料を構成する各材料の特性を取得する。ここで、取得する材料特性とは、ステップ3で記載しているように、例えば、構成する材料それぞれのヤング率、降伏応力および耐力、加工硬化指数、結晶方位、並びに応力-ひずみ曲線などである。
[Step 14: Acquisition of material properties]
As shown in step 14 (S14) of FIG. 4, the control section 20 acquires the properties of each material constituting the metal material tested in step 11 by the parameter setting process 22. FIG. Here, the material properties to be acquired are, for example, the Young's modulus, yield stress and yield strength, work hardening index, crystal orientation, and stress-strain curve of each constituent material, as described in step 3. .

[ステップ15:境界条件の設定]
図4のステップ15(S15)の通り、制御部20は、パラメータ設定処理22により、FEM解析モデルに与える境界条件を設定する。当該境界条件は、上記のステップ4と同様、金属材料を構成する構造体が破壊されるまでに実際に材料に生じると考えられる応力またはひずみ量の履歴を与える。
[Step 15: Setting Boundary Conditions]
As shown in step 15 (S15) of FIG. 4, the control unit 20 sets boundary conditions to be given to the FEM analysis model by the parameter setting process 22. FIG. Similar to step 4 above, the boundary condition gives the history of the amount of stress or strain that is actually considered to occur in the material until the structure that constitutes the metal material is destroyed.

[ステップ16:FEM解析の実行(相当塑性ひずみと応力三軸度の出力)]
図4のステップ16(S16)の通り、制御部20は、解析処理23により、作成したFEM解析モデルを用いてFEM解析を実行する。上記のステップ5と同様、作成したFEM解析モデルに対して境界条件として設定した応力またはひずみの履歴を与えることで、金属材料に応力またはひずみが加えられた際にどのように変形するかシミュレーションにより解析することができる。解析により、FEM解析モデルに与えた応力またはひずみ別に、FEM解析モデルの各位置における応力三軸度および相当塑性ひずみを算出し、出力することができる。
[Step 16: Execution of FEM analysis (output of equivalent plastic strain and stress triaxiality)]
As shown in step 16 (S16) of FIG. 4, the control unit 20 performs the FEM analysis using the created FEM analysis model by the analysis processing 23. FIG. As in step 5 above, by giving the stress or strain history set as a boundary condition to the created FEM analysis model, we can simulate how the metal material will deform when stress or strain is applied. can be analyzed. By the analysis, the stress triaxiality and the equivalent plastic strain at each position of the FEM analysis model can be calculated and output for each stress or strain applied to the FEM analysis model.

[ステップ17:マイクロボイド発生位置の相当塑性ひずみと応力三軸度の抽出]
図4のステップ17(S17)の通り、制御部20は、解析処理23により、出力された前記各位置の相当塑性ひずみと応力三軸度のデータの中から、マイクロボイド発生位置における相当塑性ひずみと応力三軸度のデータを抽出する。当該抽出データは、ステップ12にて記録した当該マイクロボイドが発生した際の塑性ひずみεmicroを、FEM解析モデルに与えてどのように変形するか解析することで算出した値を抽出する。したがって、マイクロボイドが発生した際の塑性ひずみεmicroに対応した相当塑性ひずみおよび応力三軸度のデータを抽出することができる。
[Step 17: Extraction of equivalent plastic strain and stress triaxiality at microvoid occurrence position]
As shown in step 17 (S17) in FIG. 4, the control unit 20 selects the equivalent plastic strain at the microvoid occurrence position from the data of the equivalent plastic strain and the stress triaxiality output at each position by the analysis processing 23. and extract the stress triaxiality data. The extracted data is a value calculated by applying the plastic strain ε micro when the microvoid is generated recorded in step 12 to the FEM analysis model and analyzing how it is deformed. Therefore, it is possible to extract the data of the equivalent plastic strain and the stress triaxiality corresponding to the plastic strain ε micro when microvoids are generated.

[ステップ18:抽出されたマイクロボイド発生位置の相当塑性ひずみと応力三軸度の分布作成]
図4のステップ18(S18)の通り、制御部20は、例えばグラフ作成処理等の分布作成処理24により、抽出された前記各位置の相当塑性ひずみと応力三軸度の分布を求める。具体的に例えば、これらを軸とするグラフにプロットする。該分布を求めることで、対象の構成相における、マイクロボイドが発生する、相当塑性ひずみと応力三軸度の関係(相当塑性ひずみと応力三軸度の組み合わせ値)が明らかになる。分布としてグラフを作成する場合、前記グラフは、例えば縦軸を相当塑性ひずみ(εeq)、横軸を応力三軸度(η)とすることができる。
[Step 18: Creation of distribution of equivalent plastic strain and stress triaxiality at the extracted microvoid occurrence position]
As shown in step 18 (S18) of FIG. 4, the control unit 20 obtains the distribution of the equivalent plastic strain and the stress triaxiality at each of the extracted positions by the distribution creation process 24 such as the graph creation process. Specifically, for example, they are plotted on a graph with these as axes. Obtaining the distribution clarifies the relationship between the equivalent plastic strain and the stress triaxiality (combined value of the equivalent plastic strain and the stress triaxiality) that causes microvoids in the target constituent phase. When creating a graph as a distribution, the graph can have, for example, the equivalent plastic strain (εeq P ) on the vertical axis and the stress triaxiality (η) on the horizontal axis.

[ステップ19:分布における損傷クライテリアの決定]
図4のステップ19(S19)の通り、前記グラフ等の分布において、マイクロボイドが発生する相当塑性ひずみと応力三軸度の関係から、損傷クライテリアを決定する。本決定方法の分布の一例によれば、マイクロボイドが発生する、相当塑性ひずみと応力三軸度の組み合わせ値がプロット点としてグラフに示される。また前述のとおり、相当塑性ひずみと応力三軸度が大きいほど損傷が生じやすい傾向にある。言い換えると、相当塑性ひずみと応力三軸度の小さい領域はマイクロボイドが生じにくいといえる。よって、損傷クライテリアは、上記マイクロボイドの生じにくい領域と、マイクロボイドが発生した上記プロット点の集合域との境界線として示される。前記プロット点は、多ければ多いほど精度が高まることから、例えば3点以上、更には5点以上示すことが好ましい。制御部20が、分布作成処理24により、プロットされた結果から、損傷クライテリアをグラフ上に表示してもよい。代替的に、制御部20は、分布作成処理24により、FEM解析においてFEM解析モデルに与えた応力またはひずみ別に、出力された各要素の相当塑性ひずみと応力三軸度とから、損傷クライテリアをグラフ上に表示してもよい。
[Step 19: Determine Damage Criteria in Distribution]
As shown in step 19 (S19) of FIG. 4, the damage criteria are determined from the relationship between the equivalent plastic strain at which microvoids are generated and the stress triaxiality in the distributions such as the graphs. According to one example of the distribution of this determination method, the combined values of equivalent plastic strain and stress triaxiality at which microvoids occur are plotted on the graph. In addition, as described above, the greater the equivalent plastic strain and stress triaxiality, the more likely damage is to occur. In other words, microvoids are less likely to occur in regions where the equivalent plastic strain and stress triaxiality are small. Therefore, the damage criterion is indicated as a boundary line between the region where microvoids are unlikely to occur and the clustered region of the plot points where microvoids occur. As the number of plotted points increases, the accuracy increases. Therefore, it is preferable to indicate, for example, 3 points or more, more preferably 5 points or more. The control unit 20 may display the damage criteria on a graph from the results plotted by the distribution creation process 24 . Alternatively, the control unit 20 uses the distribution creation processing 24 to graph the damage criteria from the output equivalent plastic strain and stress triaxiality of each element for each stress or strain applied to the FEM analysis model in the FEM analysis. may be displayed above.

後記する実施例では、相当塑性ひずみと応力三軸度の関係のクライテリアとして反比例の近似式を示しているがこれに限定されない。データに応じて近似式の係数も変わりうる。前記クライテリアの式は、得られたデータに応じて直線近似線等であってもよい。 In the examples described later, an inverse proportional approximation formula is shown as a criterion for the relationship between the equivalent plastic strain and the stress triaxiality, but the present invention is not limited to this. The coefficients of the approximation formula may also change depending on the data. The criterion formula may be a linear approximation line or the like depending on the obtained data.

金属材料が複相鋼であり、構成相が2種類以上である場合、損傷クライテリアは、各構成相で求める。よって複相鋼の場合、損傷クライテリアは、構成相ごとに、上記ステップ17~19を実施して決定する。例えば構成相がフェライト相とマルテンサイト相である場合、後記の実施例に例示の通り、これらフェライト相とマルテンサイト相のそれぞれにおいて、損傷クライテリアを決定することが挙げられる。 If the metal material is a dual phase steel and there are two or more constituent phases, the damage criteria are determined for each constituent phase. Therefore, in the case of dual-phase steel, the damage criteria are determined by performing steps 17 to 19 above for each constituent phase. For example, when the constituent phases are a ferrite phase and a martensite phase, the damage criteria may be determined for each of these ferrite phase and martensite phase as illustrated in the examples below.

前述の通り、従来の連続体損傷モデルを用いると、マイクロボイド発生挙動を直接評価するにあたり、決定の困難なパラメータが複数あるのに対し、本発明の損傷クライテリアの決定方法によれば、上記マイクロボイド発生の実態と、相当塑性ひずみと応力三軸度とを突き合わせて損傷クライテリアを決定しているため、簡易にかつ実態に即した閾値を決定でき、かつこの損傷クライテリアを用いて、実態に即した損傷の評価を行うことができる。 As described above, when using the conventional continuum damage model, there are multiple parameters that are difficult to determine when directly evaluating the microvoid generation behavior. Since the damage criteria are determined by comparing the actual state of void generation with the equivalent plastic strain and the stress triaxiality, it is possible to determine the threshold value easily and in line with the actual state, and using this damage criterion, it is possible to quickly determine the actual state. damage assessment can be performed.

以下、実施例を挙げて本発明をより具体的に説明する。本発明は以下の実施例によって制限を受けるものではなく、前述および後述する趣旨に合致し得る範囲で、適宜変更を加えて実施することも可能であり、それらはいずれも本発明の技術的範囲に包含される。すなわち、以下では軟質相であるフェライト相と硬質相であるマルテンサイト相の二相で構成される二相鋼を用いた実施例を示しているが、これに限定されない。また以下では、応力三軸度および相当塑性ひずみの分布の一例としてグラフの作成を行っているが、これに限定されない。 EXAMPLES Hereinafter, the present invention will be described more specifically with reference to Examples. The present invention is not limited by the following examples, and can be implemented with appropriate modifications within the scope that can match the spirit described above and below. subsumed in That is, although the examples using the duplex steel composed of the two phases of the ferrite phase as the soft phase and the martensite phase as the hard phase are shown below, the present invention is not limited to this. In addition, although the graph is created below as an example of the distribution of the stress triaxiality and the equivalent plastic strain, the present invention is not limited to this.

本実施例では、軟質層と硬質相の二相で構成される二相鋼であって、該二相鋼を構成する相の種類が同じで組織形状と相分率の異なる2種類の鋼板を対象とし、この2種類の鋼板のマイクロボイド発生挙動を、下記の通り評価した。 In this embodiment, two types of steel sheets, which are a duplex steel composed of two phases of a soft layer and a hard phase, have the same type of phases that constitute the duplex steel, but have different structural shapes and phase fractions. The microvoid generation behavior of these two types of steel sheets was evaluated as follows.

1.金属材料の準備
材料1として、C:0.063質量%、Si:0.50質量%、およびMn:1.46質量%(実績)を含む鋼を溶製し、インゴットを得てから鍛造し、熱処理を施して、マルテンサイトとフェライトの二相鋼を用意した。また材料2として、硬質相であるマルテンサイトの面積分率が材料1とは異なる、マルテンサイトとフェライトの二相鋼を用意した。上記材料1と材料2を用いて、下記の手順によりミクロ組織の損傷挙動の評価を行った。
1. Preparation of metal material As material 1, steel containing 0.063% by mass of C, 0.50% by mass of Si, and 1.46% by mass of Mn (actual result) was melted, obtained an ingot, and then forged. , heat-treated to prepare a dual-phase steel of martensite and ferrite. As material 2, a dual-phase steel of martensite and ferrite was prepared in which the area fraction of martensite, which is a hard phase, was different from that of material 1. Using the above materials 1 and 2, microstructure damage behavior was evaluated according to the following procedure.

2.金属材料のミクロ組織画像の取得
上記各材料を、顕微鏡観察面がおおよそ10mm×10mm程度のサイズに切り出し、ナイタールでエッチングしてから、走査型電子顕微鏡にて倍率1000倍で顕微鏡写真を撮影し、ミクロ組織の画像を取得した。材料1の顕微鏡写真を図5A、材料2の顕微鏡写真を図5Bに示す。図5Aにおいて、白色部分は硬質相Qであるマルテンサイト相を示し、グレー部分は軟質相Rであるフェライト相を示す。また、図5Bにおいて、薄いグレー部分は硬質相Qであるマルテンサイト相を示し、濃いグレー部分は軟質相Rであるフェライト相を示す。図5Aおよび図5Bにおいて、該顕微鏡写真の画像解析を行ったところ、硬質相Qであるマルテンサイトの面積分率は、材料1では約15%であり、材料2では約50%であった。
2. Acquisition of microstructure image of metal material Each of the above materials is cut into a size of about 10 mm x 10 mm for microscopic observation, etched with nital, and then photographed with a scanning electron microscope at a magnification of 1000. Microstructural images were acquired. A photomicrograph of material 1 is shown in FIG. 5A and a photomicrograph of material 2 is shown in FIG. 5B. In FIG. 5A, the white portion indicates the hard phase Q, which is the martensite phase, and the gray portion indicates the soft phase R, which is the ferrite phase. In FIG. 5B, the light gray portion indicates the martensite phase, which is the hard phase Q, and the dark gray portion indicates the ferrite phase, which is the soft phase R. 5A and 5B, image analysis of the micrographs revealed that the area fraction of martensite, which is the hard phase Q, was about 15% for material 1 and about 50% for material 2.

3.FEM解析モデルの作成
次に、評価システム1の制御部20にモデル作成処理21を実行させて、FEM解析に使用するFEM解析モデルとして上記材料1と材料2それぞれのミクロ組織の画像を模したメッシュモデルを作成した。FEM解析モデルは、上述のステップ2と同様の方法で作成され得る。本実施例では、ステップ2に記載したFEM解析モデルと同様、メッシュモデルは、4節点アイソパラメトリック要素を設定して、総要素数が40,279、総節点数が40,507の2Dモデルで作成されている。制御部20にモデル作成処理21を実行させることで、当該2Dモデルに対して、画像の奥行き方向へ1層分拡張し、8節点アイソパラメトリック要素により総要素数が40,279、総節点数が81,014の3Dモデルを作成する。その後、制御部20に解析処理23を実行させて、当該3Dモデルを用いてFEM解析を行わせた。
3. Creation of FEM analysis model Next, the control unit 20 of the evaluation system 1 is caused to execute the model creation process 21, and the mesh imitating the image of the microstructure of each of the materials 1 and 2 is used as the FEM analysis model to be used for the FEM analysis. created the model. A FEM analysis model can be created in a manner similar to step 2 above. In this embodiment, similar to the FEM analysis model described in step 2, the mesh model is created as a 2D model with a total number of elements of 40,279 and a total number of nodes of 40,507 by setting 4-node isoparametric elements. It is By causing the control unit 20 to execute the model creation process 21, the 2D model is expanded by one layer in the depth direction of the image, and the total number of elements is 40,279 and the total number of nodes is 40,279 by the 8-node isoparametric element. Create a 3D model of 81,014. After that, the controller 20 was caused to execute the analysis processing 23 to perform FEM analysis using the 3D model.

作成された2Dのメッシュモデルは、硬質相であるマルテンサイトに対応する領域の面積分率が、材料1では約15%となるように、材料2では約50%となるように作成されている。同様に、3Dのメッシュモデルは、硬質相であるマルテンサイトに対応する領域の体積分率が、材料1では約15%となるように、材料2では約50%となるように作成されている。 The created 2D mesh model is created so that the area fraction of the region corresponding to the hard phase martensite is about 15% for material 1 and about 50% for material 2. . Similarly, the 3D mesh model is made such that the volume fraction of the region corresponding to the hard phase martensite is about 15% for material 1 and about 50% for material 2. .

4.材料特性の取得
金属材料を構成する材料であるフェライトおよびマルテンサイトの特性を取得して制御部20のパラメータ設定処理22として入力し、FEM解析モデルにおけるフェライトに対応する領域およびマルテンサイトに対応する領域のパラメータとして設定した。ここで、取得した材料特性とは、ステップ3で記載しているように、例えば、フェライトとマルテンサイトそれぞれのヤング率、降伏応力および耐力、加工硬化指数並びに結晶方位などである。また、金属材料を構成するフェライトとマルテンサイトそれぞれに関する応力-ひずみ曲線を取得し、材料特性として制御部20に設定した。それらの応力-ひずみ曲線は、フェライトで構成されている単相鋼とマルテンサイトで構成されている単相鋼それぞれに対して引張試験を行い、取得した。
4. Acquisition of material properties The properties of ferrite and martensite, which are materials constituting the metal material, are acquired and input as parameter setting processing 22 of the control unit 20, and the regions corresponding to ferrite and martensite in the FEM analysis model was set as a parameter of Here, the acquired material properties are, for example, Young's modulus, yield stress and yield strength, work hardening index and crystal orientation of ferrite and martensite, respectively, as described in step 3. In addition, stress-strain curves for ferrite and martensite, which constitute the metal material, were obtained and set in the control unit 20 as material properties. These stress-strain curves were obtained by performing tensile tests on the single-phase steel composed of ferrite and the single-phase steel composed of martensite, respectively.

5.境界条件の設定
次に、各FEM解析モデルに与える境界条件を制御部20に設定した。当該境界条件は、上記のステップ4と同様、金属材料を構成する構造体が破壊されるまでに実際に材料に生じると考えられるひずみ量の履歴を与えた。本実施例では、後述するように一軸方向の引張試験により引張試験片に対してひずみを生じさせる。したがって、FEM解析モデルに与える境界条件として、引張試験における引っ張り方向に対応する方向をx軸(図5Aおよび図5Bにおける左右方向)、x軸に直交する方向をyおよびz軸(それぞれ図5Aおよび図5Bにおける上下方向および奥行き方向)とすると、ε=-νε=-νεとなるようにひずみ増分を与えた。ここで、εはx軸方向のひずみ増分、εはy軸方向のひずみ増分、εは、z軸方向のひずみ増分である。また、νは、ポアソン比である。本実施例では、ν=0.5として設定した。したがって、各軸の方向のひずみ増分は、x軸方向:y軸方向:z軸方向=1:-0.5:-0.5となるように設定した。
5. Setting of Boundary Conditions Next, the boundary conditions given to each FEM analysis model were set in the controller 20 . Similar to step 4 above, the boundary conditions provided a history of the amount of strain that would actually occur in the material until the structure constituting the metal material was destroyed. In this example, strain is generated in a tensile test piece by a uniaxial tensile test, as will be described later. Therefore, as boundary conditions given to the FEM analysis model, the direction corresponding to the tensile direction in the tensile test is the x axis (horizontal direction in FIGS. 5A and 5B), and the direction orthogonal to the x axis is the y and z axes (FIGS. 5A and 5B, respectively). 5B), the strain increment was given so that ε x =−νε y =−νε z . where ε x is the strain increment along the x-axis, ε y is the strain increment along the y-axis, and ε z is the strain increment along the z-axis. Also, ν is Poisson's ratio. In this embodiment, ν is set to 0.5. Therefore, the strain increment in the direction of each axis was set to x-axis direction:y-axis direction:z-axis direction=1:-0.5:-0.5.

6.FEM解析の実行(相当塑性ひずみと応力三軸度の出力)
次に、制御部20に解析処理23を実行させて、作成したFEM解析モデルを用いてFEM解析を実行させる。FEM解析では、上記のような境界条件で設定された塑性ひずみがFEM解析モデルに与えられた際、解析モデルがどのように変形するかシミュレーションにより解析することができる。それにより、解析モデルでの各位置における変位を算出し、応力三軸度および相当塑性ひずみを算出し、出力することができる。本実施例では、各位置として、FEM解析モデルの各要素の積分点における値を用いた。上記FEM解析は、境界条件として設定する塑性ひずみを変更することで、任意の塑性ひずみに対する解析を行うことができる。
6. Execution of FEM analysis (output of equivalent plastic strain and stress triaxiality)
Next, the controller 20 is caused to execute the analysis processing 23 to execute the FEM analysis using the created FEM analysis model. In the FEM analysis, when the plastic strain set by the boundary conditions as described above is applied to the FEM analysis model, it is possible to analyze by simulation how the analysis model deforms. Thereby, the displacement at each position in the analysis model can be calculated, and the stress triaxiality and the equivalent plastic strain can be calculated and output. In this embodiment, the values at the integration points of the elements of the FEM analysis model are used as the respective positions. By changing the plastic strain set as a boundary condition, the above FEM analysis can analyze any plastic strain.

7.相当塑性ひずみと応力三軸度のグラフ作成
制御部20に分布作成処理24を実行させて、上記計算結果から得られた、ミクロ組織内の各位置、すなわち全積分点の応力三軸度と相当塑性ひずみの値を、境界条件として与えたひずみ値ごとに、縦軸を相当塑性ひずみ(εeq)、横軸を応力三軸度(η)とするグラフにプロットさせた。上記プロットしたグラフとして、x軸方向の塑性ひずみが5%のときの材料1と材料2の軟質相のグラフをそれぞれ図6Aと図6B、x軸方向の塑性ひずみが5%のときの材料1と材料2の硬質相のグラフをそれぞれ図7Aと図7B、x軸方向の塑性ひずみが20%のときの材料1と材料2の軟質相のグラフをそれぞれ図8Aと図8B、x軸方向の塑性ひずみが20%のときの材料1と材料2の硬質相の2軸グラフをそれぞれ図9Aと図9Bに示す。
7. Graph creation of equivalent plastic strain and stress triaxiality The control unit 20 is caused to execute the distribution creation process 24, and each position in the microstructure obtained from the above calculation result, that is, the stress triaxiality and the equivalent The values of plastic strain were plotted on a graph with equivalent plastic strain (εeq P ) on the vertical axis and stress triaxiality (η) on the horizontal axis for each strain value given as a boundary condition. As the graphs plotted above, the graphs of the soft phases of material 1 and material 2 when the plastic strain in the x-axis direction is 5% are shown in FIGS. 6A and 6B, respectively. 7A and 7B are graphs of the hard phase of material 2 and material 2, respectively, and FIGS. Biaxial graphs of the hard phases of Material 1 and Material 2 at a plastic strain of 20% are shown in FIGS. 9A and 9B, respectively.

8.ミクロ組織の損傷挙動の評価
上記ひずみ量別、構成相別に、材料1(硬質相が少ないミクロ組織)と材料2(硬質相が多いミクロ組織)を比較した。その結果、材料1と材料2の、軟質相と硬質相の両方において、同じひずみを与えたときに、材料2の方が、応力三軸度も相当塑性ひずみもより高い値を示す傾向にあることがわかる。またこの傾向は、ひずみ5%よりもひずみ20%の方がより顕著であることがわかる。前述のとおり、相当塑性ひずみと応力三軸度が大きいほど損傷が生じやすくなる傾向にあることから、応力三軸度と相当塑性ひずみが材料1よりも高い値を示す材料2は、材料1よりも、より早くミクロ組織内の微視的な損傷が生じやすいことがこれらのグラフの対比から判断できる。
8. Evaluation of Microstructure Damage Behavior Material 1 (microstructure with less hard phase) and material 2 (microstructure with more hard phase) were compared for each strain amount and each constituent phase. As a result, in both the soft and hard phases of Material 1 and Material 2, when the same strain is applied, Material 2 tends to show higher values for stress triaxiality and equivalent plastic strain. I understand. Moreover, it can be seen that this tendency is more pronounced at a strain of 20% than at a strain of 5%. As mentioned above, the greater the equivalent plastic strain and stress triaxiality, the more likely damage is to occur. From the comparison of these graphs, it can be judged that microscopic damage in the microstructure tends to occur more quickly.

9.損傷クライテリアの決定
下記手順に沿って損傷クライテリアを決定した。
9. Determination of damage criteria Damage criteria were determined according to the following procedure.

(1)塑性変形中の複数のミクロ組織画像の取得(準連続ミクロ組織観察-引張試験)
引張試験途中の複数の段階でミクロ組織観察を行って、引張試験中の応力/ひずみと、該応力/ひずみを受けて変形したミクロ組織の撮影を行った。
詳細には、図10に示す形状の引張試験片を用意した。前記引張試験片はマルテンサイト面積率が13.7%である材料で作製した。そして、該引張試験片の一方の平面をミクロ観察用に鏡面研磨し(ナイタールエッチング有り)、他方の平面にビデオ伸び計用のシールを貼付し、上記した応力/ひずみ付与部11の一例である図11の準連続ミクロ組織観察-引張試験装置70の側面図の通り、(引張)試験片71を、チャック72で固定し、ビデオ伸び計用のシール73を貼付した面側にビデオカメラ74を配置した。そしてビデオカメラで撮影しながら引張試験を行った。引張試験の試験条件は、室温環境下、クロスヘッド変位速度1mm/minとした。(引張)試験片71のシール73貼付面とは反対側の面には、ミクロ組織観察用の研磨面75を有している。
(1) Acquisition of multiple microstructure images during plastic deformation (quasi-continuous microstructure observation-tensile test)
Microstructural observations were performed at multiple stages during the tensile test to photograph the stress/strain during the tensile test and the microstructure deformed under the stress/strain.
Specifically, a tensile test piece having the shape shown in FIG. 10 was prepared. The tensile test piece was made of a material having a martensite area ratio of 13.7%. Then, one plane of the tensile test piece is mirror-polished (with nital etching) for microscopic observation, and a video extensometer seal is attached to the other plane, and an example of the stress/strain imparting unit 11 As shown in the side view of a quasi-continuous microstructure observation-tensile test apparatus 70 in FIG. 11, a (tensile) test piece 71 is fixed with a chuck 72, and a video camera 74 was placed. A tensile test was performed while filming with a video camera. The test conditions for the tensile test were a crosshead displacement rate of 1 mm/min under a room temperature environment. The surface of the (tensile) test piece 71 opposite to the surface to which the seal 73 is applied has a polished surface 75 for microstructure observation.

試験では、まずSEM観察し、ひずみを与える前の状態のミクロ組織写真を取得した。そして、引張試験機に引張試験片71を設置して引っ張ることでひずみを与え、所定の伸び量となったところで引っ張りを止め、引張試験機から引張試験片71を取り出して、引張試験片71のミクロ組織観察用の研磨面75をSEM観察し、その後、引張試験機に引張試験片71を再度設置し、所定の伸び量となるまで引っ張ることを繰り返した。この繰り返しを行うことで、あるタイミング(ひずみ量)のときのミクロ組織を把握できる。観察位置は、図12に示す通り、引張試験片における矢印の位置とした。ミクロ組織の観察は走査型電子顕微鏡を用いて行った。ミクロ組織の観察は、同一視野を、図13に示す通り、測定開始時(無ひずみ)と、くびれ発生の前後を含めた破断直前までの複数時期に行った。上記各SEM観察では、引張試験片に各応力が付加された時点で、引張試験装置から引張試験片を一旦取り外し、試験片の研磨面のミクロ観察を行った。 In the test, SEM observation was first performed to obtain a microstructure photograph of the state before strain was applied. Then, the tensile test piece 71 is placed in a tensile tester and pulled to apply strain, stop pulling when a predetermined amount of elongation is reached, remove the tensile test piece 71 from the tensile tester, and remove the tensile test piece 71. The polished surface 75 for microstructure observation was observed with an SEM, and then the tensile test piece 71 was placed again in the tensile tester and pulled repeatedly until a predetermined amount of elongation was reached. By repeating this process, it is possible to grasp the microstructure at a certain timing (amount of strain). The observation position was the position of the arrow in the tensile test piece as shown in FIG. Observation of the microstructure was performed using a scanning electron microscope. The microstructure was observed in the same visual field at a plurality of times from the start of measurement (no strain) to immediately before fracture including before and after occurrence of constriction, as shown in FIG. 13 . In each of the above SEM observations, when each stress was applied to the tensile test piece, the tensile test piece was temporarily removed from the tensile tester, and the polished surface of the test piece was microscopically observed.

(2)ミクロ組織画像からのマイクロボイド発生位置とマイクロボイド発生時の塑性ひずみの把握
上記ミクロ観察結果の一部を図14A~図14Cに示す。図14Aは、非ひずみ時、すなわち前記図13の〔1〕の黒丸の時点のミクロ組織写真、図14Bは、前記図13の〔2〕の黒丸の時点であって、引張試験による引張方向と同一の軸方向(すなわち、x軸方向)に関する塑性ひずみ(ミクロ組織写真に表示された金属材料の組織形状等をもとに視野の変形量から推定した、ミクロ組織写真に表示されている所定領域のひずみ値)εmicroが17.1%(ビデオ伸び計用のシールの変位量から推定した公称ひずみεnominalが13.3%)時のミクロ組織写真、図14Cは、前記図13の〔3〕の黒丸の時点であって、塑性ひずみεmicroが35.5%(公称ひずみεnominalが18.2%)時のミクロ組織写真である。図14A~図14Cに示す破線で囲われた領域は、上記所定領域を表しており、横方向(x軸方向)に伸び、縦方向(y軸方向)に縮んでいることが分かる。この準連続ミクロ組織観察-引張試験では、マクロひずみ量であるx軸方向の塑性ひずみ量(εmicro)、当該塑性ひずみ量でのマイクロボイド発生位置、マイクロボイド発生起点となった構成相の種類を評価項目とした。
(2) Ascertainment of Microvoid Occurrence Location and Plastic Strain at Microvoid Occurrence from Microstructure Images Some of the microscopic observation results are shown in FIGS. 14A to 14C. 14A is a microstructure photograph at the time of no strain, that is, at the time of the black circle in [1] of FIG. 13, and FIG. 14B is the time of the black circle in [2] in FIG. Plastic strain in the same axial direction (that is, x-axis direction) (a predetermined area displayed in the microstructure photograph, which is estimated from the amount of deformation in the field of view based on the structure shape of the metal material displayed in the microstructure photograph Strain value) ε micro is 17.1% (nominal strain ε nominal estimated from the amount of displacement of the seal for the video extensometer is 13.3%). ] and a microstructure photograph when the plastic strain ε micro is 35.5% (nominal strain ε nominal is 18.2%). 14A to 14C, the area surrounded by broken lines represents the predetermined area, and it can be seen that it extends in the horizontal direction (x-axis direction) and shrinks in the vertical direction (y-axis direction). In this quasi-continuous microstructure observation-tensile test, the amount of plastic strain in the x-axis direction (ε micro ), which is the amount of macro strain, the position where microvoids are generated at the amount of plastic strain, and the type of constituent phase that is the origin of microvoids was used as an evaluation item.

準連続ミクロ組織観察-引張試験でのマイクロボイドの発生は、上記ミクロ組織写真において、図15Aの矢印または図15Bの丸囲みで示されるようなボイドを、目視で確認した。その確認結果の一例を図16~図19に示す。図16および図17は、フェライト相で発生したマイクロボイド発生位置の写真であり、図18および図19は、マルテンサイト相で発生したマイクロボイド発生位置の写真である。図16A~図16Cは引張試験片にくびれが発生する前の、フェライト相におけるマイクロボイドの発生写真であり、図16D、図16E、図17A~図17Dは引張試験片にくびれが発生した後の、フェライト相におけるマイクロボイドの発生写真である。また、図18Aが、引張試験片にくびれが発生する前の、マルテンサイト相におけるマイクロボイドの発生写真であり、図18B~図18D、図19Aおよび図19Bが、引張試験片にくびれが発生した後の、マルテンサイト相におけるマイクロボイドの発生写真である。図16~図19のいずれの写真においても、矢印で示すマイクロボイド発生位置を含む左側の写真と、該写真を拡大し、マイクロボイド発生位置を丸で示した右側の写真を示す。 Observation of quasi-continuous microstructure—The occurrence of microvoids in the tensile test was confirmed visually in the above microstructure photographs, such as voids indicated by arrows in FIG. 15A or circles in FIG. 15B. An example of the confirmation result is shown in FIGS. 16 to 19. FIG. 16 and 17 are photographs of locations where microvoids are generated in the ferrite phase, and FIGS. 18 and 19 are photographs of locations where microvoids are generated in the martensite phase. FIGS. 16A to 16C are photographs of microvoid generation in the ferrite phase before constriction occurs in the tensile test piece, and FIGS. , is a photograph of microvoid generation in the ferrite phase. In addition, FIG. 18A is a photograph of microvoid generation in the martensite phase before constriction occurs in the tensile test piece, and FIGS. 18B to 18D, FIGS. It is a photograph of generation of microvoids in the martensite phase afterward. 16 to 19, the photograph on the left including the locations where microvoids are generated indicated by arrows and the enlarged photograph on the right where the locations where microvoids are generated are indicated by circles.

評価の一例として、上記図14Cのミクロ組織写真において、上記図16~19に示すように項目を評価した結果を図20に示す。図20において、実線または点線で囲んだ部分はボイドが発生した箇所を示しており、実線の囲みはフェライト相で発生したボイドであることを示し、点線の囲みはマルテンサイト相で発生したボイドであることを示す。図20は公称ひずみεnominalが18.2%時の結果であるが、この図20と同様に、他の公称ひずみεnominal時の、ボイド発生位置、ボイド発生起点となった構成相の種類も、図16~図19の通り観察し、図20に示された様な評価結果を得た。 As an example of evaluation, FIG. 20 shows the results of evaluating the items shown in FIGS. 16 to 19 in the microstructure photograph of FIG. 14C. In FIG. 20, the portions surrounded by solid lines or dotted lines indicate the locations where voids are generated, the solid line surroundings indicate voids generated in the ferrite phase, and the dotted line surroundings indicate voids generated in the martensite phase. indicates that there is FIG. 20 shows the results when the nominal strain ε nominal is 18.2%. Similar to this FIG. , and the evaluation results shown in FIG. 20 were obtained.

(3)上記引張試験に対応したひずみでの解析結果の出力
FEM解析を行うために、制御部20にモデル作成処理21を実行させて、上記と同様の方法にてFEM解析モデルを作成した。当該モデルは、図14Aにて点線で囲われた領域を用いてイメージベースモデリングにより作成されている。試験片のマルテンサイト相の面積分率は13.7%であるため、マルテンサイト相に相当する領域の体積分率が13.7%であるFEM解析モデルを作成した。また、金属材料を構成する材料であるフェライトおよびマルテンサイトの特性をパラメータとして制御部20に設定した。
(3) Output of Analysis Results at Strains Corresponding to the Tensile Test In order to perform FEM analysis, the controller 20 was caused to execute the model creation process 21 to create an FEM analysis model in the same manner as described above. The model was created by image-based modeling using the area enclosed by the dotted line in FIG. 14A. Since the area fraction of the martensite phase in the test piece is 13.7%, an FEM analysis model was created in which the volume fraction of the region corresponding to the martensite phase is 13.7%. Also, the characteristics of ferrite and martensite, which are materials constituting the metal material, are set in the control unit 20 as parameters.

次に、各FEM解析モデルに与える境界条件をパラメータとして制御部20に設定した。当該境界条件は、上記と同様、金属材料を構成する構造体が破壊されるまでに実際に材料に生じると考えられるひずみ量の履歴を与えている。また、ε=-νε=-νεとなるようにひずみ増分を与えており、ν=0.5として設定し、各軸の方向のひずみ増分は、x軸方向:y軸方向:z軸方向=1:-0.5:-0.5となるように設定した。 Next, the boundary conditions given to each FEM analysis model were set in the controller 20 as parameters. Similar to the above, the boundary condition gives the history of the amount of strain that is actually considered to occur in the material until the structure that constitutes the metal material is destroyed. In addition, the strain increment is given so that ε x = -νε y = -νε z , and ν is set to 0.5, and the strain increment in the direction of each axis is x-axis direction: y-axis direction: z It was set so that the axial direction=1:-0.5:-0.5.

そして、制御部20に解析処理23を実行させてFEM解析を行い、上記ミクロ組織の写真から推定された各塑性ひずみεmicroと同等の塑性ひずみを上記FEM解析モデルに対して引張方向に与えた場合の、FEM解析モデルの変形を解析した。すなわち、本実施例では、FEM解析モデルの各積分点に、x軸方向に対して17.1%、35.5%およびその他の推定した塑性ひずみεmicroが与えられた際のそれぞれの塑性ひずみεmicroでのFEM解析モデルの変形を解析した。そして、解析モデルでの各位置における変位を算出し、応力三軸度および相当塑性ひずみを算出し、記憶部30に出力した。すなわち、当該塑性ひずみεmicroが試験片のミクロ組織写真の所定領域にて発生した場合における、所定領域内の任意の位置での応力三軸度および相当塑性ひずみを出力している。 Then, the control unit 20 is caused to perform the analysis processing 23 to perform FEM analysis, and a plastic strain equivalent to each plastic strain ε micro estimated from the microstructure photograph is applied to the FEM analysis model in the tensile direction. We analyzed the deformation of the FEM analysis model in the case. That is, in this example, each plastic strain when 17.1%, 35.5% and other estimated plastic strains ε micro in the x-axis direction are given to each integration point of the FEM analysis model We analyzed the deformation of the FEM analysis model at ε micro . Then, the displacement at each position in the analysis model was calculated, the stress triaxiality and the equivalent plastic strain were calculated, and output to the storage unit 30 . That is, when the plastic strain ε micro occurs in a predetermined region of the microstructure photograph of the test piece, the stress triaxiality and the equivalent plastic strain at an arbitrary position within the predetermined region are output.

(4)解析結果からの、マイクロボイド発生位置の相当塑性ひずみと応力三軸度の抽出と、抽出された相当塑性ひずみと応力三軸度のグラフ作成
制御部20に分布作成処理24を実行させて、出力された前記各位置の相当塑性ひずみと応力三軸度のデータの中から、マイクロボイド発生位置における相当塑性ひずみと応力三軸度のデータを抽出した。当該抽出データは、当該マイクロボイドが発生した際の塑性ひずみεmicroをFEM解析モデルに与えて算出した値を抽出している。そして各位置が対応する構成相別に、抽出された前記各位置の相当塑性ひずみと応力三軸度を、縦軸を相当塑性ひずみ(εeq)、横軸を応力三軸度(η)としたグラフにプロットし、各構成相において、マイクロボイドが発生する相当塑性ひずみと応力三軸度の関係を求めた。その結果を、フェライト相については図21、マルテンサイト相については図22に示す。上記図21および図22から、相当塑性ひずみと応力三軸度の値が、ひずみの増加により、早期に高くなる、すなわちグラフ上で早期にグラフ右上に推移するほど、マイクロボイドが早期に発生すると言える。
(4) Extraction of the equivalent plastic strain and stress triaxiality at the microvoid occurrence position from the analysis results, and creation of a graph of the extracted equivalent plastic strain and stress triaxiality. Then, the data of the equivalent plastic strain and the stress triaxiality at the microvoid occurrence position were extracted from the output data of the equivalent plastic strain and the stress triaxiality at each position. The extracted data is a value calculated by applying the plastic strain ε micro when the microvoid is generated to the FEM analysis model. Then, for each constituent phase corresponding to each position, the equivalent plastic strain and stress triaxiality of each extracted position were set to the equivalent plastic strain (εeq P ) on the vertical axis and the stress triaxiality (η) on the horizontal axis. Plotted on a graph, the relationship between the equivalent plastic strain at which microvoids occur and the stress triaxiality was obtained for each constituent phase. The results are shown in FIG. 21 for the ferrite phase and FIG. 22 for the martensite phase. From FIG. 21 and FIG. 22 above, the values of the equivalent plastic strain and the stress triaxiality become higher early with an increase in strain, that is, the earlier the microvoids are generated, the earlier it shifts to the upper right of the graph. I can say

(5)グラフでの損傷クライテリアの決定
図21、図22のそれぞれにおいて、プロット点の領域と、相当塑性ひずみと応力三軸度が小さくなるグラフの左斜め下方向の領域との境界線として、反比例型の近似曲線を算出した。該近似曲線は、各構成相での損傷が発生する閾値曲線であり、該曲線よりも上は損傷が生じる領域であることを示す。
(5) Determination of damage criteria in graphs In each of FIGS. 21 and 22, as a boundary line between the area of the plot points and the area in the lower left direction of the graph where the equivalent plastic strain and stress triaxiality decrease, An inversely proportional approximation curve was calculated. The approximation curve is a threshold curve at which damage occurs in each constituent phase, and indicates the area above the curve where damage occurs.

10.損傷クライテリアを用いたミクロ組織の損傷挙動の評価
FEM解析で得られた、各構成相の応力三軸度と相当塑性ひずみの関係を、上記求めた損傷クライテリアで判断することによって、各構成相の損傷挙動をより明確に判断することができる。詳細には、前記図6、図7、図8および図9のいずれにおいても、材料1よりも材料2のほうがグラフ右上のプロット点が多かった。このことから、材料2のほうが早期にマイクロボイドが発生すると考えられ、材料2のほうが破断特性は劣ると推察される。
10. Evaluation of damage behavior of microstructure using damage criteria Damage behavior can be determined more clearly. More specifically, in all of FIGS. 6, 7, 8 and 9, there were more plotted points in the upper right corner of the graph for material 2 than for material 1. FIG. From this, it is considered that microvoids are generated earlier in material 2, and it is inferred that material 2 is inferior in rupture properties.

近年、金属材料の成分組成と製造プロセスが複雑となりつつあることに伴い、金属材料のミクロ組織のスケールでは、組織や強度の分布(不均質性)が可避的に発現しやすく、該組織等の不均質性はマクロな強度特性に大きく影響する。更に、上記金属材料で形成の構造体に要求される強度特性は近年ますます厳しくなっている。しかし、本実施形態に係る方法によれば、上記組織等の不均質性の影響を正確に考慮して、実用スケールの構造体(部材)のマクロ強度特性を精度高く評価することができる。 In recent years, as the component composition and manufacturing process of metallic materials have become more complex, on the scale of the microstructure of metallic materials, the distribution (heterogeneity) of the structure and strength tends to unavoidably appear, and the structure etc. The inhomogeneity of steel greatly affects the macroscopic strength properties. Furthermore, the strength characteristics required for structures formed from the above metal materials have become more and more severe in recent years. However, according to the method according to the present embodiment, it is possible to accurately evaluate the macro-strength characteristics of a practical-scale structure (member) by accurately considering the influence of the inhomogeneity of the structure and the like.

本開示の請求項に記載の金属材料のミクロ組織の損傷挙動評価方法は、ハードウェア資源、例えば、プロセッサ、メモリ、と、ソフトウェア資源(コンピュータプログラム)との協働などによって実現される。 The damage behavior evaluation method of the microstructure of the metallic material according to the claims of the present disclosure is realized by cooperation of hardware resources, such as processors, memories, and software resources (computer programs).

本実施形態に係る方法は、より厳しい強度特性の求められる、例えば建築物、船舶、海洋構造物、橋梁、タンク等の構造体を構成する金属材料の開発に有効に活用することができる。 The method according to the present embodiment can be effectively utilized in the development of metal materials that constitute structures such as buildings, ships, offshore structures, bridges, tanks, etc., which require stricter strength characteristics.

1 評価システム
10 試験部
11 応力/ひずみ付与部
12 撮影部
20 制御部
21 モデル作成処理
22 パラメータ設定処理
23 解析処理
24 分布作成処理
30 記憶部
40 表示部
50 入力部
60 出力部
70 準連続ミクロ組織観察-引張試験装置
71 引張試験片
72 固定用チャック
73 ビデオ伸び計用のシール
74 ビデオカメラ
75 ミクロ組織観察用の研磨面
1 evaluation system 10 test unit 11 stress/strain applying unit 12 imaging unit 20 control unit 21 model creation processing 22 parameter setting processing 23 analysis processing 24 distribution creation processing 30 storage unit 40 display unit 50 input unit 60 output unit 70 quasi-continuous microstructure Observation-Tensile Testing Apparatus 71 Tensile Specimen 72 Clamping Chuck 73 Seal for Video Extensometer 74 Video Camera 75 Polished Surface for Microstructural Observation

Claims (7)

コンピュータを用いて金属材料のミクロ組織の損傷挙動を評価する方法であって、コンピュータが、
前記金属材料のミクロ組織画像を取得するステップ1と、
前記ミクロ組織画像を複数の要素に分割し、有限要素法(FEM)によって前記ミクロ組織画像を模擬したFEM解析モデルを作成するステップ2と、
前記FEM解析モデルのパラメータとして前記金属材料の材料特性を取得するステップ3と、
境界条件として、前記FEM解析モデルの各要素に付与する、前記金属材料の破断までに生じる応力若しくはひずみ、または前記応力若しくは前記ひずみの履歴を設定するステップ4と、
前記境界条件でFEM解析を実行し、前記各要素の相当塑性ひずみと応力三軸度を出力するステップ5と、
前記FEM解析において前記FEM解析モデルに与えた応力またはひずみ別に、出力された前記各要素の相当塑性ひずみと応力三軸度の分布を求めるステップ6と、
前記相当塑性ひずみと前記応力三軸度の分布から、ミクロ組織の損傷挙動を評価するステップ7と
を実行する、金属材料のミクロ組織の損傷挙動評価方法。
A computer-aided method for evaluating microstructural damage behavior of a metallic material, the computer comprising:
a step 1 of obtaining a microstructure image of the metallic material;
step 2 of dividing the microstructure image into a plurality of elements and creating an FEM analysis model simulating the microstructure image by a finite element method (FEM);
Step 3 of obtaining material properties of the metallic material as parameters of the FEM analysis model;
Step 4 of setting, as boundary conditions, the stress or strain that occurs until the metal material fractures, or the history of the stress or strain that is applied to each element of the FEM analysis model;
Step 5 of performing an FEM analysis under the boundary conditions and outputting the equivalent plastic strain and stress triaxiality of each element;
Step 6 for obtaining the distribution of the equivalent plastic strain and the stress triaxiality of each of the output elements for each stress or strain applied to the FEM analysis model in the FEM analysis;
A method for evaluating damage behavior of a microstructure of a metal material, comprising: step 7 of evaluating damage behavior of a microstructure from the distribution of the equivalent plastic strain and the stress triaxiality.
前記ステップ4は、各要素の積分点に付与する、前記応力若しくは前記ひずみ、または前記応力若しくは前記ひずみの前記履歴を設定し、
前記ステップ5は、各要素の積分点における前記相当塑性ひずみと応力三軸度を出力する、請求項1に記載の金属材料のミクロ組織の損傷挙動評価方法。
The step 4 sets the stress or the strain, or the history of the stress or the strain, to be applied to the integration point of each element,
2. The method for evaluating damage behavior of a microstructure of a metallic material according to claim 1, wherein said step 5 outputs said equivalent plastic strain and stress triaxiality at integration points of each element.
前記金属材料は、2種類以上の構成相を有し、
前記ステップ5~7を構成相ごとに実施する、請求項1または2に記載の金属材料のミクロ組織の損傷挙動評価方法。
The metal material has two or more constituent phases,
3. The method for evaluating damage behavior of a microstructure of a metallic material according to claim 1, wherein steps 5 to 7 are performed for each constituent phase.
予め用意された、前記2種類以上の構成相それぞれの応力-ひずみ曲線を、前記ステップ3にて前記材料特性として取得することを含む、請求項3に記載の金属材料のミクロ組織の損傷挙動評価方法。 Damage behavior evaluation of the microstructure of the metallic material according to claim 3, comprising obtaining the stress-strain curves of each of the two or more constituent phases prepared in advance as the material properties in the step 3. Method. 前記ステップ7で、予め定められた手順に従って損傷クライテリアを決定し、決定された前記損傷クライテリアを用いてミクロ組織の損傷有無を判断する、請求項1から4のいずれか一項に記載の金属材料のミクロ組織の損傷挙動評価方法。 5. The metal material according to any one of claims 1 to 4, wherein in step 7, a damage criterion is determined according to a predetermined procedure, and the presence or absence of microstructure damage is determined using the determined damage criterion. damage behavior evaluation method of the microstructure of. 前記損傷クライテリアの決定を、
前記金属材料からなる試験片を用意し、該試験片に、連続的にひずみを付与しながら、該試験片の同一視野のミクロ組織画像を複数取得するステップ11と、
前記複数のミクロ組織画像から、マイクロボイド発生位置と、マイクロボイド発生時の塑性ひずみ量を決定するステップ12と、
前記ミクロ組織画像を複数の要素に分割した、前記ミクロ組織画像を模擬したFEM解析モデルを作成するステップ13と、
前記FEM解析モデルのパラメータとして前記金属材料の材料特性を取得するステップ14と、
境界条件として、前記FEM解析モデルの各要素に付与する、ひずみ量を設定するステップ15と、
前記境界条件でFEM解析を実行し、前記マイクロボイド発生時の各塑性ひずみ量における前記各要素の相当塑性ひずみと応力三軸度を出力するステップ16と、
出力された前記各要素の相当塑性ひずみと応力三軸度のデータの中から、前記マイクロボイド発生位置における相当塑性ひずみと応力三軸度のデータを抽出するステップ17と、
抽出された各位置の相当塑性ひずみと応力三軸度の分布を求めるステップ18と、
前記分布において、マイクロボイドが発生する相当塑性ひずみと応力三軸度の関係から、損傷クライテリアを決定するステップ19と
を含む工程で行う、請求項5に記載の金属材料のミクロ組織の損傷挙動評価方法。
determining the damage criteria,
A step 11 of preparing a test piece made of the metal material and acquiring a plurality of microstructure images of the same field of view of the test piece while continuously applying strain to the test piece;
Step 12 of determining the microvoid generation position and the amount of plastic strain at the time of microvoid generation from the plurality of microstructure images;
step 13 of creating an FEM analysis model simulating the microstructure image by dividing the microstructure image into a plurality of elements;
a step 14 of obtaining material properties of the metallic material as parameters of the FEM analysis model;
Step 15 of setting a strain amount to be applied to each element of the FEM analysis model as a boundary condition;
A step 16 of performing FEM analysis under the boundary conditions and outputting the equivalent plastic strain and stress triaxiality of each element at each plastic strain amount at the time of microvoid generation;
a step 17 of extracting the data of the equivalent plastic strain and the stress triaxiality at the microvoid generation position from the output data of the equivalent plastic strain and the stress triaxiality of each element;
a step 18 of obtaining distributions of equivalent plastic strain and stress triaxiality at each extracted position;
6. Damage behavior evaluation of the microstructure of the metal material according to claim 5, which is performed in the process including step 19 of determining the damage criteria from the relationship between the equivalent plastic strain at which microvoids are generated and the stress triaxiality in the distribution. Method.
前記金属材料は、2種類以上の構成相を有し、
前記損傷クライテリアの決定でのステップ17~19を構成相ごとに実施する、請求項6に記載の金属材料のミクロ組織の損傷挙動評価方法。
The metal material has two or more constituent phases,
7. The method for evaluating damage behavior of a microstructure of a metal material according to claim 6, wherein steps 17 to 19 in determining the damage criteria are performed for each constituent phase.
JP2021107010A 2021-06-28 2021-06-28 Damage behavior evaluation method for microstructure of metallic materials Active JP7474225B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2021107010A JP7474225B2 (en) 2021-06-28 2021-06-28 Damage behavior evaluation method for microstructure of metallic materials

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2021107010A JP7474225B2 (en) 2021-06-28 2021-06-28 Damage behavior evaluation method for microstructure of metallic materials

Publications (2)

Publication Number Publication Date
JP2023005230A true JP2023005230A (en) 2023-01-18
JP7474225B2 JP7474225B2 (en) 2024-04-24

Family

ID=85107179

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021107010A Active JP7474225B2 (en) 2021-06-28 2021-06-28 Damage behavior evaluation method for microstructure of metallic materials

Country Status (1)

Country Link
JP (1) JP7474225B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117034430A (en) * 2023-09-11 2023-11-10 湘潭大学 Gate health monitoring method, system and computer readable storage medium based on deep learning and digital twin

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7360442B2 (en) 2005-12-29 2008-04-22 Accellent, Inc. Method for measuring and calculating tensile elongation of ductile metals
JP2012104042A (en) 2010-11-12 2012-05-31 Kobe Steel Ltd Prediction method of uniform extension and prediction program for uniform extension
JP6159551B2 (en) 2013-03-29 2017-07-05 株式会社Jsol Part breakage prediction system and part breakage prediction program
JP6232934B2 (en) 2013-10-31 2017-11-22 新日鐵住金株式会社 Material characteristic acquisition method and apparatus, program, and recording medium

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117034430A (en) * 2023-09-11 2023-11-10 湘潭大学 Gate health monitoring method, system and computer readable storage medium based on deep learning and digital twin
CN117034430B (en) * 2023-09-11 2024-01-09 湘潭大学 Gate health monitoring method, system and computer readable storage medium based on deep learning and digital twin

Also Published As

Publication number Publication date
JP7474225B2 (en) 2024-04-24

Similar Documents

Publication Publication Date Title
Nicoletto et al. Characterization of microshrinkage casting defects of Al–Si alloys by X-ray computed tomography and metallography
Tarigopula et al. A study of large plastic deformations in dual phase steel using digital image correlation and FE analysis
Baudoin et al. Influence of critical resolved shear stress ratios on the response of a commercially pure titanium oligocrystal: crystal plasticity simulations and experiment
Heripre et al. Coupling between experimental measurements and polycrystal finite element calculations for micromechanical study of metallic materials
Abbassi et al. Failure analysis based on microvoid growth for sheet metal during uniaxial and biaxial tensile tests
Ghosh et al. Crystal plasticity FE modeling of Ti alloys for a range of strain-rates. Part II: Image-based model with experimental validation
Olofsson et al. Characterisation and investigation of local variations in mechanical behaviour in cast aluminium using gradient solidification, Digital Image Correlation and finite element simulation
Li et al. On the evolution of lattice deformation in austenitic stainless steels—The role of work hardening at finite strains
JP5967321B2 (en) Stretch flange crack prediction method, stretch flange crack prediction apparatus, computer program, and recording medium
Wilson-Heid et al. Anisotropic multiaxial plasticity model for laser powder bed fusion additively manufactured Ti-6Al-4V
Fagerholt et al. Experimental and numerical investigation of fracture in a cast aluminium alloy
Shade et al. Exploring new links between crystal plasticity models and high-energy X-ray diffraction microscopy
Khan et al. Low cycle lifetime assessment of Al2024 alloy
Sharaf et al. Modeling the microstructure influence on fatigue life variability in structural steels
Gerke et al. Experiments and numerical simulation of damage and fracture of the X0-specimen under non-proportional loading paths
Cheong et al. Predicting fatigue crack initiation through image-based micromechanical modeling
Andani et al. Estimation of micro-Hall-Petch coefficients for prismatic slip system in Mg-4Al as a function of grain boundary parameters
Tang et al. On the inhomogeneous deformation behavior of magnesium alloy beam subjected to bending
JP5009222B2 (en) Method and apparatus for predicting deformation characteristics of polycrystalline material, program and recording medium
US11047812B2 (en) Macrotexture map visualizing texture heterogeneity in polycrystalline parts
Nowell et al. Investigation of fatigue crack models by micro-scale measurement of crack tip deformation
JP7474225B2 (en) Damage behavior evaluation method for microstructure of metallic materials
Farahani et al. A GTN failure analysis of an AA6061-T6 Bi-failure specimen
Balusu et al. Investigating the relationship between grain orientation and surface height changes in nickel polycrystals under tensile plastic deformation
Evrard et al. Influence of boundary conditions on bi-phased polycrystal microstructure calculation

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230901

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20240318

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20240409

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20240412

R150 Certificate of patent or registration of utility model

Ref document number: 7474225

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150