(57)【要約】【課題】摩擦攪拌接合のシミュレーション解析に必要な広域の温度条件と広域のひずみ速度条件に適用可能な材料構成則において、複数の定数を精度よく且つ簡易に同定できるようにする。【解決手段】所定の対象材料の特性を示す材料構成則の関係式に含まれる複数の定数を同定する摩擦攪拌接合解析における材料構成則の定数同定方法において、前記対象材料でなる試験片１，２を用いて所定の条件下で摩擦攪拌接合試験を行い、所定の物理量Ｔａを計測し、前記複数の定数に仮定の値が与えられた前記関係式を用いて、条件を前記摩擦攪拌接合試験と整合させて有限要素解析を行い、前記物理量Ｔｂを算出し、計測された物理量Ｔａと算出された物理量Ｔｂとの一致度を評価し、所定の同定完了条件が成立するまで、上記の評価に基づいて、前記仮定の値を変更して上記の有限要素解析と評価を繰り返し行う。【選択図】図６
2017187343.tif 000002
【特許請求の範囲】
【請求項１】
所定の対象材料の特性を示す材料構成則の関係式に含まれる複数の定数を同定する摩擦攪拌接合解析における材料構成則の定数同定方法であって、
前記対象材料でなる試験片を用いて所定の条件下で摩擦攪拌接合試験を行い、所定の物理量を計測する物理量計測工程と、
前記複数の定数に仮定の値が与えられた前記関係式を用いて、条件を前記摩擦攪拌接合試験と整合させて有限要素解析を行い、前記物理量を算出する物理量算出工程と、
前記物理量計測工程で計測された物理量と、前記物理量算出工程で算出された物理量との一致度を評価する評価工程と、を有し、
所定の同定完了条件が成立するまで、前記評価工程の評価に基づいて、前記仮定の値を変更して前記物理量算出工程と前記評価工程を繰り返し行うことを特徴とする摩擦攪拌接合解析における材料構成則の定数同定方法。
【請求項２】
前記物理量は、摩擦攪拌接合が行われているときの前記対象材料の温度であることを特徴とする請求項１に記載の摩擦攪拌接合解析における材料構成則の定数同定方法。
【請求項３】
前記摩擦攪拌接合試験は、前記試験片の被接合部に、高速で回転するツールの一端側を前記試験片の表面側から押し込むことで行われるものであり、
前記物理量計測工程では、前記表面側から見たときの前記被接合部からの距離が異なる複数箇所、及び、前記表面からの深さが異なる複数箇所で前記試験片の温度を計測することを特徴とする請求項２に記載の摩擦攪拌接合解析における材料構成則の定数同定方法。
【請求項４】
所定の対象材料の特性を示す材料構成則の関係式に含まれる複数の定数を同定する摩擦攪拌接合解析における材料構成則の定数同定システムであって、
前記対象材料でなる試験片を用いて所定の条件下で行われた摩擦攪拌接合試験で計測された前記試験片の所定の物理量の情報を取得する情報取得手段と、
前記複数の定数に仮定の値が与えられた前記関係式を用いて、前記摩擦攪拌接合試験と条件を整合させた有限要素解析を行い、前記物理量を算出する物理量算出手段と、
前記情報取得手段によって取得された物理量と、前記物理量算出手段によって算出された物理量との一致度を評価する評価手段と、
所定の同定完了条件が成立するまで、前記評価手段による評価に基づいて、前記仮定の値を変更して、前記物理量算出手段による前記物理量の算出と、前記評価手段による評価とを繰り返し行う定数同定手段と、を備えたことを特徴とする摩擦攪拌接合解析における材料構成則の定数同定システム。
【請求項５】
前記物理量は、摩擦攪拌接合が行われているときの前記対象材料の温度であることを特徴とする請求項４に記載の摩擦攪拌接合解析における材料構成則の定数同定システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、摩擦攪拌接合のシミュレーション解析に用いられる材料構成則の定数同定方法及び定数同定システムに関する。
【背景技術】
【０００２】
部材間を接合する技術の一つとして、特許文献１，２に開示された摩擦攪拌接合が知られている。図１１及び図１２に示すように、摩擦攪拌接合では、先端部にショルダ３２と該ショルダ３２から突出したプローブ３１とを備えたツール３０が用いられ、該ツール３０は、周方向Ａに高速で回転しながら、プローブ３１とショルダ３２が材料２１，２２間の突き合わせ部２４に押し込まれるように軸方向Ｂの力が加えられ、該突き合わせ部２４に沿って横方向Ｃに材料２１，２２に対して相対的に並進移動される。
【０００３】
これによって、材料２１，２２は、ツール３０との間に生じる摩擦熱によってツール３０の周囲において粘塑性状態になると共に、該粘塑性材料がツール３０の高速回転によって攪拌されることで材料発熱が発生する。この結果、材料２１，２２の突き合わせ部２４及びその周辺部では、固相線温度には達しない程度に温度上昇した固相状態で再結晶が生じ、これにより、材料２１，２２間が突き合わせ部２４に沿って接合される。これにより、材料２１，２２間には、微細化された金属組織からなる機械的性質に優れた接合部が形成され得る。
【０００４】
上記のような摩擦攪拌接合は、一部のアルミニウム合金において実用化が果たされているが、その他のアルミニウム合金や、アルミニウム合金以外の金属、或いは樹脂などにも適用可能な技術であり、今後、種々の材料への応用が期待されている。
【０００５】
摩擦攪拌接合において、ツールの回転速度や駆動トルク、材料に対するツールの並進移動速度などの条件が適切でなければ、材料の内部に欠陥が生じることがあり、これに起因して、継手強度が不足して接合部として十分な性能を有しない可能性がある。そこで、適切な条件で摩擦攪拌接合を行うためには、ツールの耐久性、ツール駆動装置の剛性、治具による材料の保持強度などを考慮しながら、接合機材の設計を行う必要がある。
【０００６】
現在、摩擦攪拌接合用の接合機材の設計や接合条件の設定にあたっては、接合機材の実機を試作して、この実機によって摩擦攪拌接合を適切に実行可能であるか否かを確認する試験を行っているのが実情である。
【０００７】
ところが、このように実機を用いた試験を行うためには、実機の試作及び試験に多大なコストと時間がかかってしまうことから、有限要素法を用いたシミュレーション解析システムを利用することで、コストの削減と時間の短縮を図ることが検討されている。
【０００８】
この場合、有限要素解析には、摩擦攪拌接合が行われる材料の特性を適切に示し得る材料構成則の関係式が用いられる必要がある。特に、摩擦攪拌接合が行われる材料部分では、粘塑性状態での材料発熱が生じるため、この現象を解析において高精度で予測し得る材料構成則の関係式が必要になる。
【０００９】
また、粘塑性状態で攪拌される材料では、ひずみ速度や温度などが部位によって大きく異なると共に、各部位においても大きく変動する。具体的には、摩擦攪拌接合が行われる材料部分において、そのひずみ速度域は例えば１０−３〜１０２ｓ−１程度の広域に亘り、温度域は例えば室温以上、固相線温度（アルミニウム合金の場合は６５０℃程度）以下の広域に亘る。
【００１０】
したがって、摩擦攪拌接合の解析には、一定温度及び一定ひずみ速度で行われる一般的な材料試験で得られるような材料構成則の関係式（所定温度での応力とひずみの関係式など）を用いることは適当でなく、温度及びひずみ速度に対する依存性が大きく、広範囲の温度域及びひずみ速度域での適用が可能な材料構成則の関係式を用いることが好ましい。そのような構成則の関係式の例として、下記の数式１に示されるJohnson-cook則の関係式、数式２及び数式３に示されるNorton-hoff則の関係式などが知られている。なお、数式２において速度依存性を示す定数「ｍ」には、数式３において温度依存性を示す定数「ｍ０、ｍ１、ｍ２」を含む右項全体が代入される。
【００１１】
【数１】
2017187343.tif 000003
【００１２】
【数２】
2017187343.tif 000004
【００１３】
数式１に示されるJohnson-cook則の関係式と、数式２，３に示されるNorton-hoff則の関係式のいずれにおいても、左項が応力σで構成され、右項が、材料の変形抵抗を示す項、材料の温度に依存する項、及び、材料のひずみ速度に依存する項を乗じた数式で構成されている。このように構成された材料構成則の関係式は、広範囲の温度域及びひずみ速度域で適用され得ることから、摩擦攪拌接合中における粘塑性材料の特性を示すものとして好適に用いることができる。
【００１４】
上記のJohnson-cook則の関係式（数式１）やNorton-hoff則の関係式（数式２，３）には複数の定数が含まれており、これらの定数は材料毎に決められるものである。したがって、このような関係式を摩擦攪拌接合の解析に用いる際には、材料構成則の関係式に含まれる定数を予め同定しておく必要がある。
【先行技術文献】
【特許文献】
【００１５】
【特許文献１】特許第２７１２８３８号公報
【特許文献２】特許第２７９２２３３号公報
【発明の概要】
【発明が解決しようとする課題】
【００１６】
ところが、上記のように材料の変形抵抗を示す項、材料の温度に依存する項及び材料のひずみ速度に依存する項を有する材料構成則の関係式の構成は複雑であることから、対象材料についての熱間ねじり試験や熱間引張試験などの材料試験によって得られる応力や、ひずみ、ひずみ速度などの値を上記の関係式に代入することで複数の方程式を立てて、複数の定数の解を求めようとしても、単純な計算によっては解を得ることができない。
【００１７】
そこで、対象材料についての材料試験によって得られる応力とひずみ速度などの計測値を用いると共に、例えば遺伝的アルゴリズムやシミュレーティッドアニーリングなどをベースとした何らかの最適化手法によって、材料構成則の関係式に含まれる複数の定数に仮の値を設定する工程、材料試験と条件を整合させた有限要素解析を行って、応力とひずみ速度などの算出値を得る工程、及び、材料試験で得られる計測値と有限要素解析で得られる算出値との一致度を評価する工程からなるサイクルを繰り返しながら、最適な定数の値の組合せを探索することで、定数の同定を行うことが考えられる。
【００１８】
しかしながら、上述したように摩擦攪拌接合によって粘塑性状態となる材料部分の温度及びひずみ速度はいずれも広範囲に亘ることから、摩擦攪拌接合の解析に用いられる材料構成則の定数を精度よく同定するためには、広域に亘る複数のひずみ速度条件及び複数の温度条件の下で材料試験を行い、これらの試験結果との対比を行うために、材料試験と同じ複数の条件下での有限要素解析を行う必要がある。
【００１９】
例えば、ひずみ速度及び温度のそれぞれについて４つ〜７つ程度の条件が設定される場合には、これらの条件の組合せはおよそ２０〜４０パターンとなり、該パターン毎に材料試験と有限要素解析を行う必要がある。そして、材料試験での計測値と解析での算出値との一致度を評価するためには、２０〜４０パターンずつの値同士を対比させる必要がある。
【００２０】
このような対比による一致度の評価を行いながら上記のJohnson-cook則の関係式（数式１）やNorton-hoff則の関係式（数式２，３）などに含まれる複数の定数の組合せを最適化手法によって同定するためには、同定が完了するまでに、定数の値の仮設定、解析及び評価からなるサイクルを数百回程度繰り返す必要がある。
【００２１】
そして、定数の値を仮設定し直す度に２０〜４０パターンの解析をやり直すことになるため、１つの対象材料についての材料構成則の定数を同定するまでに、最終的に１万回程度の解析が行われることになり、膨大な時間とコストを要してしまう。
【００２２】
さらに、以上のような定数同定方法において解析結果との対比が行われる材料試験に関しても、熱間試験中に材料の温度が設定温度を超えて上昇したり、熱間引張試験中におけるネッキングにより応力状態が変化したりする問題がある。そのため、必ずしも精度が高くない材料試験の結果を用いた評価が行われることから、材料構成則の定数を精度よく同定できないという問題もある。
【００２３】
そこで、本発明は、摩擦攪拌接合のシミュレーション解析に必要な広域の温度条件と広域のひずみ速度条件に適用可能な材料構成則において、複数の定数を精度よく且つ簡易に同定できるようにすることを課題とする。
【課題を解決するための手段】
【００２４】
前記課題を解決するため、本発明に係る摩擦攪拌接合解析における材料構成則の定数同定方法及び定数同定システムは、次のように構成したことを特徴とする。
【００２５】
まず、本願の請求項１に記載の発明は、所定の対象材料の特性を示す材料構成則の関係式に含まれる複数の定数を同定する摩擦攪拌接合解析における材料構成則の定数同定方法であって、
前記対象材料でなる試験片を用いて所定の条件下で摩擦攪拌接合試験を行い、所定の物理量を計測する物理量計測工程と、
前記複数の定数に仮定の値が与えられた前記関係式を用いて、条件を前記摩擦攪拌接合試験と整合させて有限要素解析を行い、前記物理量を算出する物理量算出工程と、
前記物理量計測工程で計測された物理量と、前記物理量算出工程で算出された物理量との一致度を評価する評価工程と、を有し、
所定の同定完了条件が成立するまで、前記評価工程の評価に基づいて、前記仮定の値を変更して前記物理量算出工程と前記評価工程を繰り返し行うことを特徴とする。
【００２６】
なお、本願の発明において、「同定完了条件」の具体的な条件は特に限定されるものでないが、例えば、同定されるべき複数の定数に関しての最適な値の組合せを探索するために所定の最適化手法が用いられる場合、該最適化手法における所定の「探索完了条件」を「同定完了条件」としてもよい。
【００２７】
ここでいう「探索完了条件」の具体例としては、複数の定数に与えられる仮定の値の変更回数が、所定の最適化手法における規定回数に達すること、又は、直前の「評価工程」において、計測された物理量と算出された物理量との一致度が所定条件を満たすと評価されることなどが挙げられる。
【００２８】
請求項２に記載の発明に係る摩擦攪拌接合解析における材料構成則の定数同定方法は、前記請求項１に記載の発明において、
前記物理量は、摩擦攪拌接合が行われているときの前記対象材料の温度であることを特徴とする。
【００２９】
請求項３に記載の発明に係る摩擦攪拌接合解析における材料構成則の定数同定方法は、前記請求項２に記載の発明において、
前記摩擦攪拌接合試験は、前記試験片の被接合部に、高速で回転するツールの一端側を前記試験片の表面側から押し込むことで行われるものであり、
前記物理量計測工程では、前記表面側から見たときの前記被接合部からの距離が異なる複数箇所、及び、前記表面からの深さが異なる複数箇所で前記試験片の温度を計測することを特徴とする。
【００３０】
また、本願の請求項４に記載の発明は、所定の対象材料の特性を示す材料構成則の関係式に含まれる複数の定数を同定する摩擦攪拌接合解析における材料構成則の定数同定システムであって、
前記対象材料でなる試験片を用いて所定の条件下で行われた摩擦攪拌接合試験で計測された前記試験片の所定の物理量の情報を取得する情報取得手段と、
前記複数の定数に仮定の値が与えられた前記関係式を用いて、前記摩擦攪拌接合試験と条件を整合させた有限要素解析を行い、前記物理量を算出する物理量算出手段と、
前記情報取得手段によって取得された物理量と、前記物理量算出手段によって算出された物理量との一致度を評価する評価手段と、
所定の同定完了条件が成立するまで、前記評価手段による評価に基づいて、前記仮定の値を変更して、前記物理量算出手段による前記物理量の算出と、前記評価手段による評価とを繰り返し行う定数同定手段と、を備えたことを特徴とする。
【００３１】
請求項５に記載の発明に係る摩擦攪拌接合解析における材料構成則の定数同定システムは、前記請求項４に記載の発明において、
前記物理量は、摩擦攪拌接合が行われているときの前記対象材料の温度であることを特徴とする。
【発明の効果】
【００３２】
本願の請求項１に記載の発明に係る摩擦攪拌接合解析における材料構成則の定数同定方法によれば、対象材料でなる試験片を用いた摩擦攪拌接合試験が実際に行われると共に、複数の定数に仮定の値が与えられた関係式を用いて、摩擦攪拌接合試験と条件が整合された有限要素解析が行われた上で、１つの条件下での摩擦攪拌接合試験で得られる物理量の計測値と、同じく１つの条件下での解析で得られる同物理量の算出値との一致度の評価に基づいて、対象材料の温度域及びひずみ速度域が広範に亘る摩擦攪拌接合に対応した材料構成則の定数を同定することができる。
【００３３】
上記の物理量の一致度の評価においては、１回の試験結果と１回の解析結果とを対比させるだけで済むため、試験の回数が１回のみで済むと共に、仮に、定数の同定に至るまでに定数の仮定の値の変更、解析及び評価からなるサイクルが数百回繰り返されたとしても、最終的な解析の回数がそれを超えることがない。
【００３４】
よって、材料試験の結果を利用する従来の同定方法では、広範な温度域及びひずみ速度域に対応させるために数十パターンの材料試験を要すると共に、これに対応する数十パターンの解析が上記のサイクル毎に行われることで最終的に一万回程度の解析が必要になるのに対して、本発明によれば、上記のように試験及び解析の回数が少なくて済むため、時間とコストを効果的に短縮できる。したがって、本発明によれば、摩擦攪拌接合の有限要素解析に適した材料構成則の定数を簡易且つ効率的に同定できる。
【００３５】
また、従来の同定方法に用いられる材料試験では、温度やひずみ速度、応力等の値を必ずしも精度よく得られない問題があるのに対して、本発明によれば、実際の摩擦攪拌接合が行われるときの計測値が用いられることで、定数の同定を高精度で行うことができる。したがって、このように高精度に同定された定数を有する材料構成則の関係式を用いて、摩擦攪拌接合の有限要素解析を精度よく行うことができる。
【００３６】
本願の請求項２に記載の発明に係る摩擦攪拌接合解析における材料構成則の定数同定方法によれば、摩擦攪拌接合試験において計測された対象材料の温度は、摩擦攪拌接合によって生じる材料発熱に対応した物理量であるため、当該計測温度と、解析で算出された温度との一致度の評価に基づいて、材料発熱が生じる摩擦攪拌接合の解析に適した材料構成則の定数を精度よく同定できる。
【００３７】
本願の請求項３に記載の発明に係る摩擦攪拌接合解析における材料構成則の定数同定方法によれば、試験片の表面側から見たときの被接合部からの距離が異なる複数箇所、及び、試験片の表面からの深さが異なる複数箇所で計測された温度と、これらに対応する各箇所で算出された温度との一致度の評価に基づいて、対象材料における被接合部からの距離や表面からの深さの違いによる温度差が生じる摩擦攪拌接合の解析に適した材料構成則の定数を精度よく同定できる。
【００３８】
本願の請求項４に記載の発明に係る摩擦攪拌接合解析における材料構成則の定数同定システムによれば、対象材料でなる試験片を用いて実際に行われた１つの条件下での摩擦攪拌接合試験で得られた物理量の計測値と、複数の定数に仮定の値が与えられた関係式が用いられると共に、摩擦攪拌接合試験と条件が整合されて行われた１つの条件下での有限要素解析で得られた同物理量の算出値との一致度の評価に基づいて、対象材料の温度域及びひずみ速度域が広範に亘る摩擦攪拌接合の解析に適した材料構成則の定数を同定することができる。そのため、請求項１に係る発明と同様の効果を得ることができる。
【００３９】
本願の請求項５に記載の発明に係る摩擦攪拌接合解析における材料構成則の定数同定システムによれば、摩擦攪拌接合試験において計測された対象材料の温度は、摩擦攪拌接合によって生じる材料発熱に対応した物理量であるため、請求項２に係る発明と同様の効果を得ることができる。
【図面の簡単な説明】
【００４０】
【図１】本発明の一実施形態に係る摩擦攪拌接合解析における材料構成則の定数同定方法における摩擦攪拌接合試験に用いられる試験片とツールを示す平面図である。
【図２】同試験片と同ツールを示す図１のII−II線断面図である。
【図３】本発明の一実施形態に係る摩擦攪拌接合解析における材料構成則の定数同定システムの全体構成を示す図である。
【図４】同システムの記憶装置の構成を示す図である。
【図５】比較例に係る定数同定動作を示すフローチャートである。
【図６】実施例に係る定数同定動作を示すフローチャートである。
【図７】定数初期値を設定する方法を説明するためのグラフである。
【図８】摩擦攪拌接合試験と条件を整合させた有限要素解析で算出された温度分布の一例を示すコンター図である。
【図９】摩擦攪拌接合試験で計測された温度と有限要素解析で算出された温度の一例を示すグラフである。
【図１０】定数データベースの一例を示すテーブルである。
【図１１】摩擦攪拌接合を説明するための斜視図である。
【図１２】摩擦攪拌接合を説明するための断面図である。
【発明を実施するための形態】
【００４１】
以下、本発明の実施形態に係る摩擦攪拌接合解析における材料構成則の定数同定方法及び定数同定システムについて、添付図面を参照しながら詳細に説明する。
【００４２】
本実施形態に係る摩擦攪拌接合解析における材料構成則の定数同定方法は、所定の対象材料の特性を示す材料構成則の関係式に含まれる複数の定数を同定する方法である。具体的に、本実施形態では、上記の数式２及び数式３に示されるNorton-hoff則の関係式に含まれる５つの定数（Ｋ、β、ｍ０、ｍ１、ｍ２）を同定する方法について説明する。
【００４３】
ただし、本発明は、Norton-hoff則の関係式に限られず、例えば上記の数式１に示されるJohnson-cook則の関係式など、摩擦攪拌接合の有限要素解析に利用可能なその他の材料構成則の関係式にも適用可能である。
【００４４】
本実施形態では、対象材料の特性を示す材料構成則の関係式に関して、以下に説明する方法によって定数が同定されることで、当該関係式を用いて、対象材料についての摩擦攪拌接合の有限要素解析を行うことが可能になる。したがって、対象材料の摩擦攪拌接合に用いられる接合機材の設計や接合条件の設定等において、有限要素法を用いたシミュレーション解析システムを利用することが可能になるため、設計等にかかるコストと時間を効果的に削減することができる。
【００４５】
［摩擦攪拌接合試験］
図１及び図２に示すように、本実施形態に係る摩擦攪拌接合解析における材料構成則の定数同定方法では、先ず、摩擦攪拌接合によって第１試験片１と第２試験片２を相互に接合させる摩擦攪拌接合試験が行われる。第１及び第２試験片１，２は、同一の対象材料からなる。
【００４６】
なお、摩擦攪拌接合試験では、一方の試験片１の材料のみを対象材料として、他方の試験片２として、異なる材料からなるものを用いてもよい。
【００４７】
第１及び第２試験片１，２は、例えば長方形の板材であり、互いに同じ大きさ及び厚みを有する。第１及び第２試験片１，２は、例えば長辺側の側面同士が付き合わされた状態で水平面上に載置され、所定の治具（図示せず）によって位置決めされた状態で、試験片１，２間の突き合わせ部３において互いに接合される。
【００４８】
摩擦攪拌接合試験に用いられる接合機材は、上下方向に延びる例えば円柱状のツール６を備えている。ツール６には、ツール６の下端面からツール６の軸心に沿って下方に突出したプローブ７と、ツール６の下端面においてプローブ７の周囲に形成されたショルダ８とが設けられている。
【００４９】
ツール６は、その軸心周りの周方向Ａに回転駆動されながら軸方向Ｂに下降されることで、試験片１，２間の突き合わせ部３にプローブ７とショルダ８とが押し込まれる。その後、ツール６の上下方向位置が維持されるように、ツール６による試験片１，２に対する軸方向Ｂ下側への押し付けと、ツール６の周方向Ａの回転駆動とが継続されながら、突き合わせ部３に沿った横方向Ｃへのツール６の並進移動が行われる。
【００５０】
ツール６の回転速度、ツール６を回転駆動するための駆動トルク、試験片１，２に対してツール６を軸方向Ｂに押し付けるためにツール６に作用させる軸力、ツール６の並進移動速度、治具による試験片１，２の位置決め強度等の各種条件は、一般的な摩擦攪拌接合で採用されるような条件に設定される。
【００５１】
以上の条件下で摩擦攪拌接合が行われるとき、試験片１，２における複数の温度計測ポイントＣＨ１〜ＣＨ６において温度が計測される。本実施形態では、第１試験片１に設定された第１〜第５ポイントＣＨ１〜ＣＨ５と、第２試験片２に設定された第６ポイントＣＨ６にそれぞれ温度センサ１８（図３参照）が取り付けられ、各温度センサ１８によって各ポイントＣＨ１〜ＣＨ６の温度が計測される。
【００５２】
図１に示すように、第１試験片１における第１〜第５ポイントＣＨ１〜ＣＨ５に関して、平面視における突き合わせ部３からの距離は、第１ポイントＣＨ１及び第３ポイントＣＨ３が第１距離Ｌ１、第２ポイントＣＨ２及び第４ポイントＣＨ４が第１距離Ｌ１よりも大きな第２距離Ｌ２、第５ポイントＣＨ５が第２距離Ｌ２よりも大きな第３距離Ｌ３とされている。また、第２試験片２における第６ポイントＣＨ６は、平面視における突き合わせ部３からの距離が第３距離Ｌ３よりも僅かに大きな第４距離Ｌ４とされている。
【００５３】
ツール６の並進移動方向Ｃにおいて、第１及び第２ポイントＣＨ１，ＣＨ２は同じ位置とされ、第３及び第４ポイントＣＨ３，ＣＨ４は、第１及び第２ポイントＣＨ１，ＣＨ２よりも下流側の同じ位置とされ、第５及び第６ポイントＣＨ５，ＣＨ６は、第３及び第４ポイントＣＨ３，ＣＨ４よりも更に下流側の同じ位置とされている。
【００５４】
また、図２に示すように、試験片１，２の表面からの深さは、第１及び第２ポイントＣＨ１，ＣＨ２がＤ１とされており、第３及び第４ポイントＣＨ３，ＣＨ４がＤ１よりも大きなＤ２とされており、第５及び第６ポイントＣＨ５，ＣＨ６は、試験片１，２の表面に位置付けられている。
【００５５】
なお、温度計測ポイントＣＨ１〜ＣＨ６の設定は、以上の構成に限られるものでなく、複数箇所に設定されるのであれば種々の変更が可能である。ただし、温度計測ポイントＣＨ１〜ＣＨ６は、平面視における突き合わせ部３からの距離が異なる複数箇所、及び、試験片１，２の表面からの深さが異なる複数箇所に設定されることが好ましく、これにより、粘塑性状態の材料部分が攪拌される摩擦攪拌接合時において、突き合わせ部３からの距離や表面からの深さに応じて異なる試験片１，２の温度を計測することができる。
【００５６】
各温度計測ポイントＣＨ１〜ＣＨ６において計測される温度は、並進移動方向Ｃにおけるツール６の位置に応じて変化するが、本実施形態では、並進移動方向Ｃにおいて、ツール６のプローブ７の軸心が各温度計測ポイントＣＨ１〜ＣＨ６と同じ位置に位置するときに計測された温度のデータが、材料構成則の定数の同定に用いられる。
【００５７】
以上のように摩擦攪拌接合試験が行われた後、該試験で計測された温度のデータを用いて、次に説明する材料構成則の定数同定システムによる定数同定動作が行われる。
【００５８】
［摩擦攪拌接合解析における材料構成則の定数同定システムの概要］
図３は、摩擦攪拌接合解析における材料構成則の定数同定システムの中心となるコンピュータ１０の構成を示す図である。このコンピュータ１０は、ＣＰＵ等の中央演算装置１１と、メモリまたはハードディスク等の記憶装置１２と、キーボード、マウス等の入力装置１３と、ディスプレイ、プリンタ等の出力装置１４とを有する。
【００５９】
中央演算装置１１は、記憶装置１２にアクセス可能に構成されていると共に、入力装置１３から送られる信号、及び、上述の摩擦攪拌接合試験において試験片１，２（図１及び図２参照）の温度を検知する温度センサ１８から送られる信号等を入力可能に構成されている。
【００６０】
中央演算装置１１は、記憶装置１２に記憶されている情報と、入力装置１３等を介して入力された情報とを用いて、各種演算処理を実行したり、演算処理によって得られた情報を記憶装置１２に保存したりする。また、中央演算装置１１は、必要に応じて、出力装置１４に制御信号を出力して、ディスプレイの表示や印刷などを行う。
【００６１】
なお、温度センサ１８は、必ずしも中央演算装置１１に接続されなくてもよく、この場合、摩擦攪拌接合試験で計測された温度のデータは、入力装置１３等を介して記憶装置１２に保存されてもよい。
【００６２】
［記憶装置］
図４に示すように、記憶装置１２は、材料構成則の定数の同定に用いられる各種プログラムを記憶するプログラム記憶部１２ａと、各種プログラムの動作に用いられる各種データを記憶するデータ記憶部１２ｂとを有している。
【００６３】
プログラム記憶部１２ａには、定数同定動作を全般的に行う定数同定プログラム、摩擦攪拌接合試験と条件を整合させた有限要素解析を行う摩擦攪拌接合解析プログラム、摩擦攪拌接合解析プログラムによる有限要素解析に用いられる材料構成則の関係式の定数に与えられる仮定の値（以下、「定数仮定値」という）の初期値（以下、「定数初期値」という）を設定する定数初期値設定プログラム、及び、定数同定プログラムによって同定された定数を所定の定数データベースに登録する定数登録プログラムなどが記憶されている。
【００６４】
データ記憶部１２ｂには、上記の摩擦攪拌接合試験で計測された温度のデータが記憶される計測温度データファイル、摩擦攪拌接合解析プログラムによる有限要素解析によって算出された温度のデータが記憶される算出温度データファイル、摩擦攪拌接合解析プログラムによる有限要素解析に用いられる定数仮定値のデータが記憶される定数仮定値データファイル、定数初期値設定プログラムによる定数初期値の設定に用いられる各種参考データが記憶される参考データファイル、及び、既に材料構成則の定数同定が完了している各種材料に関する定数のデータが登録された定数データベースなどが設けられている。
【００６５】
［摩擦攪拌接合解析における材料構成則の定数同定システムの動作］
次に、材料構成則の定数同定システムの動作例について、図６に示す定数同定動作が行われる実施例と、図５に示す定数同定動作が行われる比較例とを対比させながら説明する。実施例と比較例のいずれにおいても、上記の数式２及び数式３に示されるNorton-hoff則の関係式に含まれる５つの定数（Ｋ、β、ｍ０、ｍ１、ｍ２）を同定するための動作について説明する。
【００６６】
［比較例］
比較例に係る摩擦攪拌接合解析における材料構成則の定数同定方法では、図５に示す定数同定動作が行われる前に、予め、対象材料について、熱間ねじり試験や熱間引張試験などの材料試験が行われ、この材料試験によって、応力とひずみの関係を示すデータ、及び、応力とひずみ速度の関係を示すデータが得られる。
【００６７】
この材料試験におけるひずみ速度の条件としては、摩擦攪拌接合によって粘塑性状態となる材料部分に生じ得る広範囲のひずみ速度に対応させるべく、例えば１．０×１０−３〜１．０×１０２ｓ−１程度の広域に亘る少なくとも４つ〜７つ程度の条件が設定される。なお、設定されたひずみ速度条件を実現するためには、材料試験で用いられるアクチュエータを、対象材料のひずみ速度に対応する速度で駆動させればよい。
【００６８】
また、温度の条件としては、摩擦攪拌接合時において、摩擦熱や材料発熱によって固相線温度（アルミニウム合金の場合は６５０℃程度）付近まで温度上昇する被接合部の温度と、被接合部から遠い材料部分の温度との差が大きくなることに対応させるべく、例えば３００〜６００℃程度の広域に亘る少なくとも４つ〜７つ程度の条件が設定される。
【００６９】
したがって、合計２０〜４０程度の条件での材料試験が行われることになり、温度条件毎に、応力とひずみの関係を示すデータと、応力とひずみ速度の関係を示すデータとが得られ、これらのデータは、コンピュータによる図５に示す定数同定動作に用いられる。
【００７０】
図５に示す比較例の制御動作には、主として、例えば遺伝的アルゴリズムをベースとした公知の最適化手法によって定数の値の最適な組合せを探索するように構成されたプログラムが用いられる。
【００７１】
この制御動作では、先ず、ステップＳ１０１において、上述の２０〜４０パターンの材料試験で得られたデータが取得される。具体的には、温度条件毎に得られた応力とひずみの関係を示すデータ、及び、温度条件毎に得られた応力とひずみ速度の関係を示すデータが取得される。
【００７２】
続くステップＳ１０２では、ステップＳ１０１で取得されたデータのうち、応力とひずみ速度の関係を示すデータに基づいて、同じく応力とひずみ速度の関係を示すNorton-hoff則の関係式（数式２，３）に含まれる５つの定数についての定数初期値が設定される。
【００７３】
具体的には、材料試験で得られた応力とひずみ速度の関係（ステップＳ１０１）と、Norton-hoff則の関係式（数式２，３）で示される応力とひずみ速度の関係との一致度が最も高くなるような定数の値の組合せが最小二乗近似等によって求められ、該組合せを構成する各定数の値が定数初期値（定数仮定値の初期値）に設定される。
【００７４】
次のステップＳ１０３では、ステップＳ１０２で設定された定数初期値を含むNorton-hoff則の関係式（数式２，３）を用いて、上述した２０〜４０パターンの材料試験に条件を整合させた有限要素解析が行われる。すなわち、材料試験と同じく、温度とひずみ速度の条件が異なる２０〜４０パターンの有限要素解析が行われる。これらの解析によって、１つのパターン毎に、応力とひずみの関係が算出される。
【００７５】
続くステップＳ１０４では、材料試験で得られた応力とひずみの関係（ステップＳ１０１）と、有限要素解析で得られた応力とひずみの関係（ステップＳ１０３）とが１つのパターン毎に対比される。つまり、ステップＳ１０４では、試験結果と解析結果との対比が２０〜４０パターン行われる。これらの対比によって、試験結果と解析結果との一致度が最小二乗近似等によって評価される。
【００７６】
続くステップＳ１０５では、図５の制御動作を行うプログラムで採用された最適化手法における所定の探索完了条件が成立したか否かが判定される。探索完了条件の具体例としては、例えば、後述するステップＳ１０６（定数仮定値の変更）の実行回数が、前記最適化手法で規定された所定回数に達することが挙げられる。
【００７７】
ただし、探索完了条件は、定数同定動作を行うプログラムに採用される最適化手法の種類に応じて異なるものであり、探索完了条件の別の例としては、直前のステップＳ１０４において、試験結果と解析結果との一致度が所定条件を満たすと評価されることなどが挙げられる。
【００７８】
ステップＳ１０５の判定の結果、探索完了条件が成立していない場合は、ステップＳ１０６において定数仮定値が変更された上で、ステップＳ１０３の有限要素解析がやり直されて、該有限要素解析の結果について、材料試験の結果との一致度が再び評価される（ステップＳ１０４）。
【００７９】
このような定数仮定値の変更（ステップＳ１０６）、有限要素解析（ステップＳ１０３）、及び、試験結果と解析結果との一致度の評価（ステップＳ１０４）のサイクルは、ステップＳ１０５の判定において探索完了条件が成立するまで繰り返される。
【００８０】
このサイクルにおける各ステップＳ１０６では、これ以前に設定された定数仮定値とこれに対応するステップＳ１０４の評価とに基づいて、上記の最適化手法に従って決められた値が、新たな定数仮定値として設定される。
【００８１】
また、このサイクルにおいて、各ステップＳ１０３では、直前のステップＳ１０６で再設定された定数仮定値を含むNorton-hoff則の関係式（数式２，３）を用いて、材料試験の条件と整合した２０〜４０パターンの有限要素解析が行われ、各ステップＳ１０４では、直前のステップＳ１０３で得られた２０〜４０パターンの解析結果が材料試験の結果と対比されて、両結果の一致度が評価される。
【００８２】
ステップＳ１０５の判定の結果、上記の最適化手法における探索完了条件が成立すると、ステップＳ１０７において、Norton-hoff則の関係式（数式２，３）に含まれる５つの定数（Ｋ、β、ｍ０、ｍ１、ｍ２）の値の組合せが、例えば、上記のサイクルの繰り返しの中で最も高い一致度の評価（ステップＳ１０４）が得られたときの定数仮定値の組合せに同定されて、定数同定動作が終了する。
【００８３】
以上で説明した比較例の方法において、２０〜４０パターンの全ての条件について試験結果と解析結果との一致度が高くなるような定数（Ｋ、β、ｍ０、ｍ１、ｍ２）の最適な組合せを見つけ出すためには、上述した定数仮定値の変更（ステップＳ１０６）、有限要素解析（ステップＳ１０３）、及び、試験結果と解析結果との一致度の評価（ステップＳ１０４）のサイクルを数百回程度繰り返す必要がある。
【００８４】
また、比較例の方法では、ステップＳ１０６において定数仮定値を再設定する度に、ステップＳ１０３において２０〜４０パターンの解析を行うことになるため、最終的にステップＳ１０７において定数の同定がなされるまでに、１万回程度の解析が行われることになる。仮に１回の解析に１時間かかる場合、全ての解析が終わるまでに１３ヶ月程度かかることになり、膨大な時間とコストを要してしまう。
【００８５】
さらに、比較例では、熱間ねじり試験や熱間引張試験などの材料試験で得られるデータが用いられるが、熱間試験中には材料の温度が設定温度を超えて上昇しやすい問題があると共に、熱間引張試験中にはネッキングにより応力状態が変化しやすい問題がある。そのため、材料試験で得られるデータの精度は必ずしも高くなく、このように精度が得られ難いデータを用いた評価に基づいて定数の同定が行われる比較例の方法では、定数を必ずしも精度よく同定できないという問題もある。
【００８６】
［実施例］
これに対して、本実施形態では、図６に示す実施例に係る制御動作が上述のコンピュータ１０（図３参照）によって行われる。
【００８７】
実施例に係る制御動作には、主として、記憶装置１２のプログラム記憶部１２ａ（図４参照）に記憶された定数同定プログラムが用いられ、この他にも、定数初期値設定プログラムや摩擦攪拌接合解析プログラム、定数登録プログラムなどが適宜併せて用いられる。
【００８８】
定数同定プログラムは、例えば遺伝的アルゴリズムをベースとした公知の最適化手法によって、Norton-hoff則の関係式（数式２，３）に含まれる５つの定数（Ｋ、β、ｍ０、ｍ１、ｍ２）の値の最適な組合せを探索するように構成されている。また、定数同定プログラムは、定数の同定に際して、定数初期値設定プログラムや摩擦攪拌接合解析プログラム、定数登録プログラムなど、他のプログラムを必要に応じて動作させて、これらのプログラムと連係して動作する。
【００８９】
図６に示す制御動作では、先ず、ステップＳ１において、記憶装置１２のデータ記憶部１２ｂ（図４参照）に記憶された各種データが取得される。
【００９０】
ステップＳ１で取得されるデータには、上記の摩擦攪拌接合試験における計測温度のデータ、及び、次のステップＳ２における定数初期値の設定に用いられる各種参考データが含まれる。
【００９１】
ステップＳ１において、計測温度のデータは、データ記憶部１２ｂに記憶された計測温度データファイル（図４参照）から取得される。ここで取得される計測温度のデータは、図１に示す並進移動方向Ｃにおいて、プローブ７の軸心が各温度計測ポイントＣＨ１〜ＣＨ６と同じ位置に位置するときに各温度計測ポイントＣＨ１〜ＣＨ６で計測された温度のデータである。
【００９２】
ステップＳ１において、各種参考データは、データ記憶部１２ｂに記憶された参考データファイル（図４参照）から取得される。ここで取得される参考データは、Norton-hoff則の関係式（数式２，３）に含まれる５つの定数（Ｋ、β、ｍ０、ｍ１、ｍ２）について定数初期値を設定するために用いられる各種情報のデータである。
【００９３】
具体的には、例えば、定数Ｋ以外の４つの定数（β、ｍ０、ｍ１、ｍ２）の値となり得る各範囲に関して、経験的に知られている対象材料又はこれに類似する材料の温度依存性及び速度依存性の変動幅から定められる範囲の上限値及び下限値のデータと、後述のように定数Ｋの定数初期値の設定に用いられるJIS Z 2241の引張試験（室温、低速）の結果（応力とひずみ速度の関係）のデータとが、上記の参考データとして用いられる。
【００９４】
ただし、参考データの構成はこれに限定されるものでなく、例えば、既に同定が完了している種々の材料の定数のデータが、上記の参考データとして用いられてもよい。
【００９５】
続くステップＳ２では、ステップＳ１で取得された参考データに基づいて、Norton-hoff則の関係式（数式２，３）に含まれる５つの定数（Ｋ、β、ｍ０、ｍ１、ｍ２）についての定数初期値が設定される。
【００９６】
例えば、４つの定数（β、ｍ０、ｍ１、ｍ２）の各定数初期値は、ステップＳ１で取得された上限値と下限値との平均値に設定される。このようにしてＫ以外の定数の値が仮に設定されたNorton-hoff則の関係式（数式２，３）に従えば、例えば図７における領域Ｐで示される応力とひずみ速度の関係が得られる。
【００９７】
そして、このようにＫ以外の定数の値が仮に設定されたNorton-hoff則の関係式（数式２，３）に、ステップＳ１で取得された引張試験の結果における応力とひずみ速度の値（図７における点Ｑ）を代入することで、残りの定数Ｋの定数初期値が算出される。このようにして５つの定数（Ｋ、β、ｍ０、ｍ１、ｍ２）の値の全てが仮に設定されたNorton-hoff則の関係式（数式２，３）に従えば、例えば図７における曲線Ｒで示される応力とひずみ速度の関係が得られる。
【００９８】
ただし、ステップＳ２において、定数初期値の設定方法は特に上記の方法に限定されるものでなく、例えば、対象材料に類似する材料に関して、既に定数（Ｋ、β、ｍ０、ｍ１、ｍ２）の同定が完了している材料が複数存在する場合は、これらの材料の定数の値に基づく計算によって得られた値（例えば、平均値又は中央値）を、対象材料についての定数初期値として設定してもよい。
【００９９】
ステップＳ２で設定された定数初期値は、データ記憶部１２ｂにおける定数仮定値データファイル（図４参照）に記憶される。
【０１００】
次のステップＳ３では、プログラム記憶部１２ａに記憶された摩擦攪拌接合解析プログラム（図４参照）によって、有限要素解析が行われる。ステップＳ３の有限要素解析は、ステップＳ２で設定された定数初期値を含むNorton-hoff則の関係式（数式２，３）を用いて、上述の摩擦攪拌接合試験に条件を整合させて行われる。
【０１０１】
ステップＳ３の有限要素解析では、例えば、図８（ａ）に示される平面視でのコンター図、及び、図８（ｂ）に示される仮想ツール１０６の並進移動方向から見たコンター図のように、並進移動中の仮想ツール１０６の周辺における仮想対象材料１０１の温度が部位毎に算出される。これによって得られる図８（ａ）及び図８（ｂ）に示すような温度分布のデータは、データ記憶部１２ｂにおける算出温度データファイル（図４参照）に記憶される。
【０１０２】
続くステップＳ４では、ステップＳ１で取得された摩擦攪拌接合試験における計測温度のデータと、ステップＳ３の有限要素解析による算出温度のデータとが対比され、両者の一致度が評価される。
【０１０３】
具体的には、例えば図９に示すように、上述した摩擦攪拌接合試験における試験片１，２の第１〜第６ポイントＣＨ１〜ＣＨ６（図１，２参照）での計測温度Ｔａと、ステップＳ３の有限要素解析で得られた温度分布（図８参照）における各ポイントＣＨ１〜ＣＨ６に相当する箇所での算出温度Ｔｂとが対比され、これらの各箇所での対比に基づいて、総合的な計測温度Ｔａと算出温度Ｔｂとの一致度が最小二乗近似等によって評価される。
【０１０４】
続くステップＳ５では、定数同定プログラムで用いられる最適化手法における所定の探索完了条件が成立したか否かが判定される。具体的な探索完了条件は、例えば、後述するステップＳ６（定数仮定値の変更）の実行回数が、前記最適化手法で規定された所定回数に達することとされる。
【０１０５】
ただし、具体的な探索完了条件は、定数同定プログラムで採用される最適化手法の種類に応じて異なるものであり、探索完了条件の別の例としては、例えば、直前のステップＳ４において、計測温度Ｔａと算出温度Ｔｂとの一致度が所定条件を満たすと評価されることなどが挙げられる。
【０１０６】
ステップＳ５の判定の結果、探索完了条件が成立していない場合は、ステップＳ６において定数仮定値が変更された上で、ステップＳ３において有限要素解析がやり直されると共に、続くステップＳ４において、やり直された有限要素解析の結果について、摩擦攪拌接合試験の結果との一致度が再び評価される（ステップＳ４）。
【０１０７】
このような定数仮定値の変更（ステップＳ６）、有限要素解析（ステップＳ３）、及び、計測温度Ｔａと算出温度Ｔｂとの一致度の評価（ステップＳ４）のサイクルは、ステップＳ５の判定において探索完了条件が成立するまで繰り返される。
【０１０８】
このサイクルにおける各ステップＳ６では、これ以前に設定された定数仮定値とこれに対応するステップＳ４の評価とに基づいて、上記の最適化手法に従って決められた値が、新たな定数仮定値として設定される。
【０１０９】
また、このサイクルにおける各ステップＳ３では、直前のステップＳ６で再設定された定数仮定値を含むNorton-hoff則の関係式（数式２，３）を用いて、摩擦攪拌接合試験と条件を整合させた有限要素解析が行われ、各ステップＳ４では、直前のステップＳ３の有限要素解析で得られた算出温度Ｔｂが摩擦攪拌接合試験での計測温度Ｔａと対比されて、両温度Ｔａ，Ｔｂの一致度が評価される。
【０１１０】
ステップＳ５の判定の結果、上記の最適化手法における探索完了条件が成立すると、ステップＳ７において、Norton-hoff則の関係式（数式２，３）に含まれる５つの定数（Ｋ、β、ｍ０、ｍ１、ｍ２）の値の組合せが、例えば、上記のサイクル（ステップＳ６、ステップＳ３〜ステップＳ５）の繰り返しの中で最も高い一致度の評価（ステップＳ４）が得られたときの定数仮定値の組合せに同定されて、定数同定動作が終了する。
【０１１１】
ただし、ステップＳ７における定数の同定の仕方は上記のものに限られるものでなく、定数同定プログラムで採用される最適化手法の種類に応じて異なる。ステップＳ７におけるその他の同定の仕方として、例えば、最後のステップＳ６で設定された最終の定数仮定値の組合せが定数の値として同定されてもよい。
【０１１２】
ステップＳ７で同定された定数（Ｋ、β、ｍ０、ｍ１、ｍ２）の値は、プログラム記憶部１２ａに記憶された定数登録プログラム（図４参照）によって、データ記憶部１２ｂに記憶された定数データベース（図４参照）に、対象材料の特性を示す材料構成則の定数として新たに登録されて、定数同定動作が終了する。
【０１１３】
例えば図１０に示すように、定数データベースには、Norton-hoff則の関係式（数式２，３）に含まれる５つの定数（Ｋ、β、ｍ０、ｍ１、ｍ２）について同定された値が材料毎に登録される。したがって、定数の同定を終えた材料については、定数データベースに登録された定数の値を用いることで、摩擦攪拌接合の有限要素解析を精度よく行うことができる。
【０１１４】
以上で説明した実施例の方法によれば、材料試験の結果を利用する比較例の方法（図５参照）とは異なり、対象材料でなる試験片１，２（図１及び図２参照）を用いた摩擦攪拌接合試験が実際に行われると共に、５つの定数（Ｋ、β、ｍ０、ｍ１、ｍ２）に定数仮定値が与えられたNorton-hoff則の関係式（数式２，３）を用いて、摩擦攪拌接合試験と条件が整合された有限要素解析（図６のステップＳ３）が行われた上で、１つの条件下での摩擦攪拌接合試験で得られる計測温度Ｔａと、同じく１つの条件下での有限要素解析で得られる算出温度Ｔｂとの一致度の評価（図６のステップＳ４）に基づいて、対象材料の温度域及びひずみ速度域が広範に亘る摩擦攪拌接合に対応したNorton-hoff則の関係式（数式２，３）の定数（Ｋ、β、ｍ０、ｍ１、ｍ２）を同定することができる。
【０１１５】
計測温度Ｔａと算出温度Ｔｂとの一致度の評価（図６のステップＳ４）においては、１回の試験結果と１回の解析結果とを対比させるだけで済むため、図６に示す定数同定動作の前に予め行われる試験の回数が１回のみで済むと共に、図６に示す定数同定動作において、仮に、定数の同定（図６のステップＳ７）に至るまでに定数仮定値の変更（図６のステップＳ６）、解析（図６のステップＳ３）及び評価（図６のステップＳ４）からなるサイクルが数百回繰り返されたとしても、最終的な解析の回数がそれを超えることがない。
【０１１６】
しかも、上記のサイクルにおいて、１回の評価（図６のステップＳ４）に際しては、１パターンの試験結果と解析結果とを対比させるだけで済むため、数十パターンの試験結果と解析結果とを対比させる比較例に比べて、定数の同定（図６のステップＳ７）に至るまでに必要な上記サイクルの回数を低減しやすくなり、例えば百回程度のサイクルによって高精度の同定をなし得る。
【０１１７】
よって、上述の比較例の方法では、数十パターンの材料試験を要すると共に、これに対応する数十パターンの解析（図５のステップＳ１０３）が上記のようなサイクル毎に行われることで最終的に一万回程度の解析が必要になるのに対して、実施例の方法によれば、上記のように試験及び解析の回数が少なくて済むため、時間とコストを効果的に短縮できる。
【０１１８】
したがって、本実施例によれば、摩擦攪拌接合の有限要素解析に適切に用いられ得るNorton-hoff則の関係式（数式２，３）の定数（Ｋ、β、ｍ０、ｍ１、ｍ２）を簡易且つ効率的に同定できる。
【０１１９】
また、比較例の方法に用いられる材料試験では、温度やひずみ速度、応力等の値を必ずしも精度よく得られないのに対して、実施例の方法によれば、摩擦攪拌接合が行われるときに対象材料に生じる摩擦熱や材料発熱に応じて上昇する温度を用いた評価（図６のステップＳ４）に基づいて、摩擦熱や材料発熱が生じる摩擦攪拌接合の解析に適したNorton-hoff則の関係式（数式２，３）の定数（Ｋ、β、ｍ０、ｍ１、ｍ２）の同定を高精度で行うことができる。
【０１２０】
しかも、実施例では、図１及び図２に示す摩擦攪拌接合試験において、試験片１，２の表面側から見たときの突き合わせ部３からの距離Ｌ１，Ｌ２，Ｌ３，Ｌ４が異なる複数箇所ＣＨ１〜ＣＨ６、及び、試験片１，２の表面からの深さＤ１，Ｄ２が異なる複数箇所ＣＨ１〜ＣＨ６で計測された温度Ｔａと、有限要素解析における対応箇所について算出された温度Ｔｂとの一致度の評価に基づいて、対象材料における被接合部からの距離や表面からの深さの違いによる温度差が生じる摩擦攪拌接合の解析に適したNorton-hoff則の関係式（数式２，３）の定数（Ｋ、β、ｍ０、ｍ１、ｍ２）を精度よく同定できる。
【０１２１】
したがって、本実施形態によれば、以上のように高精度に同定された定数（Ｋ、β、ｍ０、ｍ１、ｍ２）を有するNorton-hoff則の関係式（数式２，３）を用いて、摩擦攪拌接合の有限要素解析を精度よく行うことができる。
【０１２２】
以上、上述の実施形態を挙げて本発明を説明したが、本発明は上述の実施形態に限定されるものではない。
【０１２３】
例えば、上述の実施形態では、Norton-hoff則の関係式（数式２，３）の定数（Ｋ、β、ｍ０、ｍ１、ｍ２）を同定する例を説明したが、本発明が適用される具体的な材料構成則はこれに限定されるものでなく、例えば、Johnson-cook則の関係式（数式１）の定数（Ａ，Ｂ，ｎ，Ｃ，ｍ）を同定する場合にも本発明を好適に適用可能である。
【０１２４】
また、上述の実施形態では、図６に示す定数同定動作において、遺伝的アルゴリズムをベースにした最適化手法が用いられる例を説明したが、本発明において、定数同定動作に用いられる最適化手法はこれに限られず、その他の種々の最適化手法を用いることができる。
【０１２５】
さらに、上述の実施形態では、摩擦攪拌接合試験での計測結果と、有限要素解析での算出結果との対比において、温度の一致度を評価する例（図６のステップＳ４）を説明したが、本発明において、一致度の評価がなされる物理量は「温度」に限られるものでなく、例えば、ツール６（図１及び図２参照）に作用する回転トルク、軸方向Ｂにおいてツール６に作用する軸力、又は、並進移動方向Ｃにおいてツール６に作用する抵抗などであってもよい。
【０１２６】
これら３つの変形例のうち、ツール６に作用する回転トルクの一致度が評価される第１の例では、ツール６の回転速度が一定に維持されながら摩擦攪拌接合が行われるときにツール６に作用する回転トルクについて、摩擦攪拌接合試験での計測値と、これと条件を整合させた有限要素解析での算出値との一致度が評価されるようにすればよい。この場合、摩擦攪拌接合試験における回転トルクの計測及び有限要素解析における回転トルクの算出は、例えば、材料に対するツール６の押し込み深さが異なる複数のタイミング、及び、並進移動方向Ｃにおけるツール６の位置が異なる複数のタイミングで行われることが好ましく、これによって得られる複数の計測値と複数の算出値との対比に基づいて一致度が評価されることで、定数の同定を精度よく行うことができる。
【０１２７】
ツール６に作用する軸力の一致度が評価される第２の例では、材料に対してツール６が所定深さまで押し込まれているときにツール６に軸方向Ｂに作用する軸力について、摩擦攪拌接合試験での計測値と、これと条件を整合させた有限要素解析での算出値との一致度が評価されるようにすればよい。この場合、摩擦攪拌接合試験における軸力の計測及び有限要素解析における軸力の算出は、例えば、材料に対するツール６の押し込み深さが一定に維持された状態において、並進移動方向Ｃにおけるツール６の位置が異なる複数のタイミングで行われることが好ましく、これによって得られる複数の計測値と複数の算出値との対比に基づいて一致度が評価されることで、定数の同定を精度よく行うことができる。
【０１２８】
並進移動方向Ｃにおいてツール６に作用する抵抗の一致度が評価される第３の例では、材料に対するツール６の押し込み深さが一定に維持されると共に、ツール６の回転速度が一定に維持されながら、並進移動方向Ｃにおけるツール６の並進移動が行われているときに、並進移動方向Ｃにおいてツール６に作用する抵抗について、摩擦攪拌接合試験での計測値と、これと条件を整合させた有限要素解析での算出値との一致度が評価されるようにすればよい。この場合、摩擦攪拌接合試験における抵抗の計測及び有限要素解析における抵抗の算出は、例えば、並進移動方向Ｃにおけるツール６の位置が異なる複数のタイミングで行われることが好ましく、これによって得られる複数の計測値と複数の算出値との対比に基づいて一致度が評価されることで、定数の同定を精度よく行うことができる。
【産業上の利用可能性】
【０１２９】
以上のように、本発明によれば、摩擦攪拌接合のシミュレーション解析に必要な広域の温度条件と広域のひずみ速度条件に適用可能な材料構成則において、複数の定数を精度よく且つ簡易に同定できるようにすることが可能となるから、摩擦攪拌接合のシミュレーション解析分野において好適に利用される可能性がある。
【符号の説明】
【０１３０】
１，２    試験片
３        突き合わせ部
６        ツール
７        プローブ
８        ショルダ
１０      コンピュータ
１１      中央演算装置
１２      記憶装置
１２ａ    プログラム記憶部
１２ｂ    データ記憶部
１３      入力装置
１４      出力装置
１８      温度センサ
１０１    仮想対象材料
１０６    仮想ツール
ＣＨ１    第１ポイント
ＣＨ２    第２ポイント
ＣＨ３    第３ポイント
ＣＨ４    第４ポイント
ＣＨ５    第５ポイント
ＣＨ６    第６ポイント
Ｔａ      計測温度
Ｔｂ      算出温度
【図１】
2017187343.tif 000005
【図２】
2017187343.tif 000006
【図３】
2017187343.tif 000007
【図４】
2017187343.tif 000008
【図５】
2017187343.tif 000009
【図６】
2017187343.tif 000010
【図７】
2017187343.tif 000011
【図８】
2017187343.tif 000012
【図９】
2017187343.tif 000013
【図１０】
2017187343.tif 000014
【図１１】
2017187343.tif 000015
【図１２】
2017187343.tif 000016