JP5403256B2 - 多結晶金属材の疲労強度評価装置と方法 - Google Patents
多結晶金属材の疲労強度評価装置と方法 Download PDFInfo
- Publication number
- JP5403256B2 JP5403256B2 JP2009263090A JP2009263090A JP5403256B2 JP 5403256 B2 JP5403256 B2 JP 5403256B2 JP 2009263090 A JP2009263090 A JP 2009263090A JP 2009263090 A JP2009263090 A JP 2009263090A JP 5403256 B2 JP5403256 B2 JP 5403256B2
- Authority
- JP
- Japan
- Prior art keywords
- crystal
- crystal grain
- stress
- fatigue strength
- modulus
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Description
そのような多結晶金属材を対象とする疲労強度評価法は、例えば、特許文献1、非特許文献1,2等に開示されている。
非特許文献1,2は、多結晶金属の結晶粒と結晶方位を方位像顕微鏡で特定した後に、高サイクル疲労試験を実施してすべり挙動を観察し、さらに有限要素法(FEM)を用いてすべり発生部における応力を解析したものである。
一方、非特許文献1,2により、多数の結晶粒間相互の変形拘束が強固になっている場合において、高サイクル疲労で粒界近傍に現われる特有なすべりは、その拘束に起因するせん断応力の上昇によるものであることが示された。
しかし、これらの従来技術では、多結晶金属材の疲労強度評価手段が提案されていないため、その確立が要望されていた。
すなわち、本発明の目的は、実在する多結晶金属材について異方性を加味した疲労強度を評価する装置と方法を提供することにある。
前記結晶粒マップと各結晶粒の結晶方位から有限要素モデルを作成するモデル作成手段と、
前記有限要素モデルを用いて各結晶粒に作用する応力分布を算出するFEM解析手段と、
前記応力分布と結晶方位から各結晶粒の結晶すべり面のせん断応力を算出するせん断応力算出手段と、
前記応力分布と前記せん断応力から多結晶金属の疲労強度を評価する疲労強度評価手段とを備え、
前記疲労強度評価手段により、
各結晶粒に対して、多結晶のヤング率で規格化した長手方向の規格化ヤング率と、シュミット因子とを求め、
各結晶粒に対して、前記規格化ヤング率とシュミット因子の積をクラック発生指数として算出し、
前記クラック発生指数が全体に対して大きい結晶粒をクラック予測結晶として疲労強度を評価する、ことを特徴とする多結晶金属材の疲労強度評価装置が提供される。
モデル作成手段により、前記結晶粒マップと各結晶粒の結晶方位から有限要素モデルを作成し、
FEM解析手段により、前記有限要素モデルを用いて各結晶粒に作用する応力分布を算出し、
せん断応力算出手段により、前記応力分布と結晶方位から各結晶粒の結晶すべり面のせん断応力を算出し、
疲労強度評価手段により、前記応力分布と前記せん断応力から、
各結晶粒に対して、多結晶のヤング率で規格化した長手方向の規格化ヤング率と、シュミット因子とを求め、
各結晶粒に対して、前記規格化ヤング率とシュミット因子の積をクラック発生指数として算出し、
前記クラック発生指数が全体に対して大きい結晶粒をクラック予測結晶として疲労強度を評価する、ことを特徴とする多結晶金属材の疲労強度評価方法が提供される。
単一の結晶にある応力が作用するとき、前述のすべり系に対してすべりを起こさせようとするせん断応力(以下「すべり応力」と呼ぶ)は、結晶と応力の幾何学的関係を用いることで算出できる。
単軸応力状態、つまり垂直応力のうちひとつがσ、他の応力成分5つが全てゼロの場合には図1(A)の模式図で示したようなすべり面1(法線n)、すべり方向(方向ベクトルd)と荷重F、断面積Aの関係において、数1の式(1)ですべり応力τが算出できる。ここでσは公称応力である。cosθ・cosφはシュミット因子と呼ばれ、最大で0.5の値を持つ。
なおここですべり系座標系とは、すべり面法線n、すべり方向ベクトルd、t(=n×d)の3方向を基準軸に持つ座標系のことである。この場合のイメージ図を図1(B)に示す。図中では応力は矢印で表示しているが、実際にはテンソルのためベクトルのように内積で投影することはできず、座標変換を行なう必要がある。
この図において、本発明の疲労強度評価装置は、結晶方位測定装置10、結晶方位解析プログラム12、有限要素解析プログラム14、すべり面せん断応力算出プログラム16、及びプログラム14,16をインストールしたコンピュータ18からなる。
結晶方位測定装置10は、例えばEBSP解析装置であり、結晶方位解析プログラム12を内蔵する結晶方位解析手段に相当する。
プログラム14,16をインストールしたコンピュータ18は、モデル作成手段、FEM解析手段、せん断応力算出手段及び疲労強度評価手段に相当する。
結晶方位取得ステップS1では、結晶方位解析手段により、実在する多結晶金属材の解析対象面の結晶粒マップと各結晶粒の結晶方位を取得する。
モデル作成ステップS2では、モデル作成手段により、前記結晶粒マップと各結晶粒の結晶方位から有限要素モデルを作成する。
応力分布算出ステップS3では、FEM解析手段により、前記有限要素モデルを用いて各結晶粒に作用する応力分布を算出する。
せん断応力算出ステップS4では、せん断応力算出手段により、前記応力分布と結晶方位から各結晶粒の結晶すべり面のせん断応力を算出する。
疲労強度評価ステップS5では、疲労強度評価手段により、前記応力分布と前記せん断応力から多結晶金属の疲労強度を評価し出力する。
図4は、単結晶FEMモデルの説明図である。この図において、単結晶FEMモデルは1辺が10mmの立方体であり、底面にZ拘束、原点Oの節点にXY拘束、A節点にY拘束、B節点にX拘束を与えている。また応力はZ方向に+100MPaを付与している。
図5は単結晶FEMモデルの結晶方位模式図である。この図において、(A)を結晶方位1、(B)を結晶方位2、(C)を結晶方位3と呼ぶ。
この結晶方位1〜3に対し、所定のミラー指数で結晶方位を与えて、FEM解析を実施した。表1は、FEM解析結果とその理論値を示している。
また、12個のすべり系のすべり応力を算出したところ、FEM解析結果から座標変換して得られたすべり応力は、単軸応力下でシュミット因子により生じるすべり応力の理論値と高い精度で一致した。従って、上述した単結晶FEM解析は、弾性応力解析と同様、すべり応力の算出においても高い精度が得られることがわかる。
単結晶の場合、表1から明らかなように、結晶方位が異なるとヤング率および変形挙動が大きく異なる。そのため、異なる方位の結晶粒が隣り合う場合、その境界には特異な応力分布が発生することが予想され、隣り合った結晶粒の方位が大きく異なる場合に顕著な応力分布変化が起こると考えられる。
そこで、結晶方位の異なる単結晶要素集合を2個結合して、仮想2結晶の弾性応力解析を行なった。
また、図7は、2結晶FEMモデルの結晶方位模式図である。この図において、(A)は結晶方位1,2を組み合わせた2結晶FEMモデルA、(B)は結晶方位1,3を組み合わせた2結晶FEMモデルB、(C)は結晶方位2,3を組み合わせた2結晶FEMモデルCである。
これらのモデルA,B,Cに対し、FEM解析を実施し、結晶粒界の応力分布とすべり応力分布を理論値と比較した。その結果を以下に説明する。
FEM解析で得られた長手方向応力分布は、Z方向の公称応力は100MPaであるにも関わらず、結晶粒界付近には113〜129MPaの長手応力が発生していた。これは結晶粒が相互に拘束しあうことで応力分布の不均一が生じ、公称応力よりも高い部分が生じたと考えられる。
また単純な2結晶FEMモデルの検討で応力が約30%上乗せされることが明らかになったが、実際の多結晶体では3次元的に拘束が生じることから、さらに大きな応力の不均一が生じ、より高い応力が発生することが予測される。
FEM解析で得られたすべり応力分布は、2結晶FEMモデルAでは、シュミット因子で算出される理論的最大値(50MPa)を越える51MPaという大きなすべり応力が発生した。また、2結晶FEMモデルBでは、理論すべり応力値が27MPaに対して42MPa(156%)という大きなすべり応力が発生した。さらに2結晶FEMモデルCでも理論値より大きなすべり応力が発生しており、結晶粒界の一部では理論値よりも大きなすべり応力が発生することが確かめられた。
以上の結果から、複数の結晶粒が隣り合っている場合の応力分布とすべり応力分布について、上述したFEM解析は、単結晶FEMモデルと同様、2結晶FEMモデルにおいても高い精度が得られることが確認された。
上述した結晶方位解析手段により、実在する多結晶金属材の解析対象面の結晶粒マップと各結晶粒の結晶方位を取得した。
図8は、EBSP測定で得られた結晶粒マップである。この図では、結晶粒界を図中の線で描いている。従って、線で囲まれた領域は各結晶粒に相当する。また、同時にFEM解析により各結晶粒の長手方向応力分布を得た。
なお、多結晶金属材(試験片)は、Ni合金を用い、試験片形状が正方形断面4mm×4mm、長さは40mmのものを使用した。
この図において、結晶粒番号は解析対象領域以外の番号も含んでいることから、7〜20という中途半端な値である。また、クラックは、後述する疲労試験で発生した位置に太線で示している。
結晶粒11と19は長手方向のヤング率が低く、応力の負担が少ないと予想され、FEM解析による結果はそれに合致する結果となった。そしてその周囲の結晶粒は応力負担が大きくなっており、特に結晶粒12番と20番はヤング率が高いため、高い応力が発生していた。
上述した試験片を用いて4点曲げ疲労試験を、5kNモジュール試験機を用いて荷重制御で実施した。
疲労試験後、光学顕微鏡で観察したところ、特定の場所に固執すべり帯(PSB:Parsistent Slip Band)が発生し、これがクラックに成長する様子が確認できた。
試験後のEBSP観察によりクラックの分布と結晶粒分布の位置関係を調査した(図9参照)ところ、ほとんどのクラックが結晶粒界手前あるいは粒界を過ぎた位置で停滞している様子が確かめられた。これは結晶粒が異なるとすべり変形の幾何学的が変化してしまうため、大きな抵抗になると考えられる。図9でCとDのクラックは連結して成長し破断に至ったが、いずれも結晶粒内でクラックが成長し、疲労試験の最終段階で連結した。
図9のA、B、C、Dの位置のクラックは、いずれもFEM解析で応力が高いところに発生し、FEM解析と試験結果の整合性が得られた。
結晶粒12番、13番、14番、20番は長手方向のヤング率が高く、FEM解析でも長手方向に高い応力が発生していた。12番以外の粒ではクラックが発生しており、十分な転位の運動が生じたと考えられる。特に20番は隣り合った19番のヤング率が低いため、応力負担が高く、疲労試験の早い段階で複数のクラックが発生していた。
以下では結晶の原子配列に起因する異方性、すなわち結晶学的なすべり応力の分布と、疲労クラックの対応を検討する。
結晶粒9番、10番、13番、14番、17番、20番が高く、そのうち9番、10番以外にはクラックが発生しており、「規格化ヤング率×シュミット因子」の値がクラック発生領域の目安になると言える。
複数発生したクラックの中で、CとDの位置(それぞれ結晶粒14番と13番の粒内に対応)のクラックは結晶粒界を横切って連結、成長し、破断に至った。このように図10(C)で示した「長手方向規格化ヤング率×シュミット因子」の値が大きい粒が並んだ領域は他になく、すべり変形が生じやすい結晶粒が並んでいるため、他のクラックよりも優先的に成長が進み、破壊に至ったと考えられる。
これにより、クラックの発生しやすい結晶粒を特定し、容易に疲労強度を評価することができる。
これにより、単純応力ではない場合でも実在する多結晶金属材について異方性を加味した疲労強度を評価することができる。
12 結晶方位解析プログラム、
14 有限要素解析プログラム、
16 すべり面せん断応力算出プログラム、
18 コンピュータ
Claims (2)
- 実在する多結晶金属材の解析対象面の結晶粒マップと各結晶粒の結晶方位を取得する結晶方位解析手段と、
前記結晶粒マップと各結晶粒の結晶方位から有限要素モデルを作成するモデル作成手段と、
前記有限要素モデルを用いて各結晶粒に作用する応力分布を算出するFEM解析手段と、
前記応力分布と結晶方位から各結晶粒の結晶すべり面のせん断応力を算出するせん断応力算出手段と、
前記応力分布と前記せん断応力から多結晶金属の疲労強度を評価する疲労強度評価手段とを備え、
前記疲労強度評価手段により、
各結晶粒に対して、多結晶のヤング率で規格化した長手方向の規格化ヤング率と、シュミット因子とを求め、
各結晶粒に対して、前記規格化ヤング率とシュミット因子の積をクラック発生指数として算出し、
前記クラック発生指数が全体に対して大きい結晶粒をクラック予測結晶として疲労強度を評価する、ことを特徴とする多結晶金属材の疲労強度評価装置。 - 結晶方位解析手段により、実在する多結晶金属材の解析対象面の結晶粒マップと各結晶粒の結晶方位を取得し、
モデル作成手段により、前記結晶粒マップと各結晶粒の結晶方位から有限要素モデルを作成し、
FEM解析手段により、前記有限要素モデルを用いて各結晶粒に作用する応力分布を算出し、
せん断応力算出手段により、前記応力分布と結晶方位から各結晶粒の結晶すべり面のせん断応力を算出し、
疲労強度評価手段により、前記応力分布と前記せん断応力から、
各結晶粒に対して、多結晶のヤング率で規格化した長手方向の規格化ヤング率と、シュミット因子とを求め、
各結晶粒に対して、前記規格化ヤング率とシュミット因子の積をクラック発生指数として算出し、
前記クラック発生指数が全体に対して大きい結晶粒をクラック予測結晶として疲労強度を評価する、ことを特徴とする多結晶金属材の疲労強度評価方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009263090A JP5403256B2 (ja) | 2009-11-18 | 2009-11-18 | 多結晶金属材の疲労強度評価装置と方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009263090A JP5403256B2 (ja) | 2009-11-18 | 2009-11-18 | 多結晶金属材の疲労強度評価装置と方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2011108032A JP2011108032A (ja) | 2011-06-02 |
JP5403256B2 true JP5403256B2 (ja) | 2014-01-29 |
Family
ID=44231421
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2009263090A Active JP5403256B2 (ja) | 2009-11-18 | 2009-11-18 | 多結晶金属材の疲労強度評価装置と方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5403256B2 (ja) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5974222B2 (ja) * | 2012-08-16 | 2016-08-23 | 国立研究開発法人産業技術総合研究所 | 顕微分光測定シミュレーション法 |
JP6586842B2 (ja) * | 2015-09-24 | 2019-10-09 | 日本製鉄株式会社 | 金属材料の疲労強度評価方法 |
JP2017161265A (ja) * | 2016-03-07 | 2017-09-14 | 新日鐵住金株式会社 | 単結晶材の材料パラメータ導出装置及び導出方法 |
CN110727989B (zh) * | 2018-06-28 | 2023-08-04 | 中车大同电力机车有限公司 | 结构疲劳强度分析方法、装置以及计算机可读存储介质 |
CN109271713B (zh) * | 2018-09-25 | 2023-02-28 | 重庆大学 | 考虑晶体微结构力学的齿轮接触疲劳分析方法 |
CN109948216B (zh) * | 2019-03-12 | 2023-01-03 | 华东理工大学 | 总应变能密度修正的缺口件低周疲劳预测方法 |
CN110705137B (zh) * | 2019-08-22 | 2023-04-21 | 中车青岛四方机车车辆股份有限公司 | 一种应力幅值和均值的确定方法以及装置 |
CN112199780B (zh) * | 2020-10-28 | 2022-04-15 | 中国舰船研究设计中心 | 一种过渡锻环等厚缩比强度及疲劳联合试验模型设计方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08209273A (ja) * | 1995-02-02 | 1996-08-13 | Hitachi Ltd | 耐熱材料用Ni基合金 |
AU2003304096A1 (en) * | 2003-05-12 | 2004-11-26 | Nihon University | Method for predicting fatigue life of spot-welded structure |
JP2008197852A (ja) * | 2007-02-10 | 2008-08-28 | Phifit Kk | 塑性加工における被加工材の組織発展の解析装置、解析システムおよび記録媒体 |
-
2009
- 2009-11-18 JP JP2009263090A patent/JP5403256B2/ja active Active
Also Published As
Publication number | Publication date |
---|---|
JP2011108032A (ja) | 2011-06-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5403256B2 (ja) | 多結晶金属材の疲労強度評価装置と方法 | |
Zhang et al. | AELAS: Automatic ELAStic property derivations via high-throughput first-principles computation | |
Kalidindi et al. | Macroscopic shape change and evolution of crystallographic texture in pre-textured FCC metals | |
Chen et al. | Active slip system identification in polycrystalline metals by digital image correlation (DIC) | |
Kysar et al. | High strain gradient plasticity associated with wedge indentation into face-centered cubic single crystals: geometrically necessary dislocation densities | |
Anahid et al. | Dwell fatigue crack nucleation model based on crystal plasticity finite element simulations of polycrystalline titanium alloys | |
Xue et al. | The influence of indenter tip radius on the micro-indentation hardness | |
Lin et al. | Crystal plasticity modeling of cyclic deformation for a polycrystalline nickel-based superalloy at high temperature | |
Shi et al. | Constitutive modeling and failure mechanisms of anisotropic tensile and creep behaviors of nickel-base directionally solidified superalloy | |
Zhang et al. | Crystallographic analysis for fatigue small crack growth behaviors of a nickel-based single crystal by in situ SEM observation | |
CN112883602B (zh) | 一种多尺度疲劳裂纹萌生寿命仿真预测方法 | |
Rabbolini et al. | Fatigue crack growth in Haynes 230 single crystals: an analysis with digital image correlation | |
Venkataramani et al. | Crystal plasticity based FE model for understanding microstructural effects on creep and dwell fatigue in Ti-6242 | |
JP2016162381A (ja) | 残留応力推定方法及び残留応力推定装置 | |
JPWO2016002880A1 (ja) | 伸びフランジ割れ予測方法、伸びフランジ割れ予測装置、コンピュータープログラム、及び記録媒体 | |
Bugat et al. | Microstructure and damage initiation in duplex stainless steels | |
P. Jirandehi et al. | General quantification of fatigue damage with provision for microstructure: A review | |
Biswas et al. | Interaction between a notch and cylindrical voids in aluminum single crystals: experimental observations and numerical simulations | |
Yeratapally et al. | Investigation of fatigue crack incubation and growth in cast MAR-M247 subjected to low cycle fatigue at room temperature | |
US10429357B2 (en) | Model creation method and device, and inspection device using the same | |
JP2009264797A (ja) | 多結晶材料の変形特性予測方法及び装置、並びにプログラム及び記録媒体 | |
Dawson et al. | Mechanical metrics of virtual polycrystals (MechMet) | |
Xue et al. | Strain energy density method for estimating fracture toughness from indentation test of 0Cr12Mn5Ni4Mo3Al steel with Berkovich indenter | |
Brommesson et al. | 3D grain structure modelling of intergranular fracture in forged Haynes 282 | |
Deka et al. | Effect of constraint and latent hardening ratio on the plastic flow around a crack tip in a hardening FCC single crystal |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20120926 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20130516 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20130522 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20130709 |
|
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: 20131002 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20131015 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 5403256 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |