■キーワード：■結晶塑性解析，その場回折試験，格子ひずみ■■ 1m =)=1g )++=)))(() ここで，(0は，それぞれすべり系αの分解せん断応力，ひずみ速度依存性指数，および参照せん断0sgn)(−)は加工硬化挙動を1exp−，および(1)を次式で仮定する． １．研究の目的と背景■計算機の高速化やメモリの大容量化に伴い，結晶塑性構成式を導入した数値解析手法の利用が拡大している．特に塑性異方性が顕著な材料においては，微視的変形機構の活動度を適切に考慮することで，集合組織の発達や巨視的変形挙動の異方性を合理的に表現可能であることから，今後も応用が拡大するものと考えられる．また，材料科学分野においては，結晶塑性有限要素法(CPFEM)１，２）や結晶塑性高速フーリエ変換(CPFFT)３～５）を用いた不均一変形挙動の解析が盛んに行われている． 結晶塑性解析を実施する上で問題になることの■つとして，材料パラメータ同定の困難さがある．特に等価ではない複数の変形機構が存在する場合や複相合金の場合は，異なる材料パラメータの組合せで同一の巨視的変形挙動が表現できてしまう場合もあり，付加的なフィッティング対象を用いた系統的なパラメータ同定手法の構築が求められる．付加的フィッティング対象としては，デジタル画像相関法(DIC)等により取得される不均一変形に関するデータや集合組織発達データ，様々な負荷条件で系統的に実施した材料試験データ，さらに放射光や中性子線を利用したその場回折試験により得られる格子ひずみ発達に関するデータ等がある６～９）．特に格子ひずみ発達は，結晶配向ごとに得られる変形中の塑性変形機構の活動が反映された平均情報であるため，塑性変形に関連する材料パラメータの同定に用いられることが多い．しかしながら，格子ひずみ発達を再現するためのパラメータ同定は未だに試行錯誤的に行われており，多大な時間と労力を要している． 本研究では，より現実的な結晶塑性解析を実施するために必要な材料パラメータ同定を効率的に実施する手法の確立を目的として，格子ひずみ発達をフィッティング対象としたパラメータ同定手法の構築を試みた．具体的には，パラメータ同定における試行錯誤的なプロセスをできるだけ自動化することを目的として，勾配法およびメタヒューリスティックアプローチの1つである遺伝的アルゴリズム(GA)による最適化の導入を検討した． ■２．解析方法■■ ・■■結晶塑性構成式■■本研究ではPeirceら１）により提案された速度依存型結晶塑性構成則を導入した大変形有限要素法を用いた．すべり系αのせん断ひずみ速度(熊本大学■大学院先端科学研究部■（ ■ ■年度■一般研究開発助成■■■■ ■ ■■ ■■■■）■教授■眞山■剛■()()g)，m，および()ˆdd ここで，()はすべり系間の相互作用を表現する係数であるが，本研究ではすべての組合せが等価であると仮定して，()ˆは，全すべり系の累積すべり量Γの関数である次式Voce則を用いた． ()()1()ˆ0()ひずみ速度を表している．速度依存性が十分に小さいとき，式(1)中()gはすべり系αの臨界分解せん断応力(CRSS)と近似的に一致する．すべり系レベルでの加工硬化挙動を表現する()gの発展則として，本研究では下記の硬化則を用いた． ()()本研究では解析対象を複相マグネシウム(Mg)合金である長周期積層(LPSO)型Mg-Zn-Y合金１０，１１）とした．本合金はα-Mg相とLPSO相からなる二相合金であり，α-Mg相は商用Mg合金と同様，底面すべり系，柱面すべり系，二次錐面
元のページ ../index.html#147