多軸疲労寿命評価方法
【課題】多軸負荷による疲労損傷を高精度に予測可能な多軸疲労寿命評価方法を提供する。
【解決手段】評価対象の構造物の各時刻での応力またはひずみ状態を基に、各時刻での主応力または主ひずみに対して直交する面である主応力・主ひずみ面をそれぞれ求め、求めた各時刻での主応力・主ひずみ面のそれぞれに対して、当該主応力・主ひずみ面に作用する各時刻での主応力または主ひずみの垂直成分の値をそれぞれ算出すると共に、当該算出した値と、予め求めておいた単軸負荷時のＳＮ線図とに基づき、各主応力・主ひずみ面での累積損傷評価を行い、当該累積損傷評価で最も損傷が大きいと評価された主応力・主ひずみ面を評価対象の主応力・主ひずみ面として、疲労寿命の評価を行う。
【解決手段】評価対象の構造物の各時刻での応力またはひずみ状態を基に、各時刻での主応力または主ひずみに対して直交する面である主応力・主ひずみ面をそれぞれ求め、求めた各時刻での主応力・主ひずみ面のそれぞれに対して、当該主応力・主ひずみ面に作用する各時刻での主応力または主ひずみの垂直成分の値をそれぞれ算出すると共に、当該算出した値と、予め求めておいた単軸負荷時のＳＮ線図とに基づき、各主応力・主ひずみ面での累積損傷評価を行い、当該累積損傷評価で最も損傷が大きいと評価された主応力・主ひずみ面を評価対象の主応力・主ひずみ面として、疲労寿命の評価を行う。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、航空機エンジン、過給器などの回転機械や、圧力容器などを含む一般構造物の疲労寿命を評価する際に用いられる多軸疲労寿命評価方法に関するものである。
【背景技術】
【０００２】
実機部材の受ける荷重は、時刻と共に変化するランダム荷重である。したがって、主応力や主ひずみの方向（主軸という）が変化することも想定される。主軸方向を変化させた多軸疲労試験結果によれば、せん断ひずみと軸ひずみの位相あるいはせん断応力と軸応力の位相がずれるに従い、疲労寿命は単軸の寿命（等価ひずみあるいは等価応力を用いて求めた疲労寿命）よりも低下していくことが報告されている。
【０００３】
こうした問題点を解決するために、負荷経路（主軸方向変化を含む負荷経路）を考慮して、主軸の変化にも対応できる多軸疲労寿命評価方法がいくつか提案されている。
【０００４】
なかでも、非特許文献１，２では、評価対象とする期間中で主応力・主ひずみが最も大きな値となる主応力・主ひずみ面に着目して、主軸の変化量を定量化する方法、および負荷経路（主軸方向変化を含む負荷経路）に基づいた多軸疲労寿命評価方法（ＩＳ（Itoh-Sakane）法）が提案されている。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】伊藤隆基、「非比例多軸低サイクル疲労寿命評価モデル」、材料、社団法人日本材料学会、２００１年１２月１５日、ｖｏｌ．５０、Ｎｏ．１２、ｐｐ．１３１７−１３２２
【非特許文献２】Takamoto Itoh、Tomohiko Ozaki、Toru Amaya、and Masao Sakane、「DETERMINATION OF STRESS AND STRAIN RANGES UNDER NON-PROPORTIONAL CYCLIC LOADING」、8th International Conference on Multiaxial Fatigue & Fracture、２００７年
【非特許文献３】Takamoto Itoh、Masao Sakane、Takahiro Hata、Naomi Hamada、「A design procedure for assessing low cycle fatigue life under proportional and non-proportional loading」、International Journal of Fatigue 28、２００６年、ｐｐ．４５９−４６６
【発明の概要】
【発明が解決しようとする課題】
【０００６】
上述のように、従来のＩＳ法では、評価対象とする期間中（例えば１サイクル中）で主応力・主ひずみが最大となる面に着目して、疲労損傷を評価してきた。つまり、従来のＩＳ法では、主応力・主ひずみが最大となる面が、最も疲労損傷が大きいとして評価を行ってきた。
【０００７】
しかしながら、実機に作用する応力・ひずみの状態によっては、主応力・主ひずみが最大となる面が、最も疲労損傷が大きいとは限らないという問題がある。
【０００８】
例えば、１サイクル中に大きな軸負荷がｎ回作用し、小さいねじり負荷がＮ回作用する場合について考えると、ｎ≒Ｎであれば、主応力・主ひずみが最大となる面が最も疲労損傷が大きい面となる。しかし、ｎ＜＜Ｎである場合、すなわち、大きい軸負荷が１回作用する間に、多数回の小さいねじり負荷が作用するような場合においては、回数の多いねじり負荷による疲労損傷が大きくなり、主応力・主ひずみが最大となる面が最も疲労損傷が大きい面とならない場合もある。
【０００９】
このような場合に、従来通り主応力・主ひずみが最大となる面を評価対象として疲労寿命の評価を行うと、最も疲労損傷が大きい面でない面を評価対象としていることとなり、危険側の評価となってしまう。
【００１０】
そこで、本発明の目的は、上記課題を解決し、多軸負荷による疲労損傷を高精度に予測可能な多軸疲労寿命評価方法を提供することにある。
【課題を解決するための手段】
【００１１】
本発明は上記目的を達成するために創案されたものであり、評価対象の構造物の各時刻での応力またはひずみ状態を基に、各時刻での主応力または主ひずみに対して直交する面である主応力・主ひずみ面をそれぞれ求め、求めた各時刻での主応力・主ひずみ面のそれぞれに対して、当該主応力・主ひずみ面に作用する各時刻での主応力または主ひずみの垂直成分の値をそれぞれ算出し、当該算出した値と、予め求めておいた単軸負荷時のＳＮ線図とに基づき、各主応力・主ひずみ面での累積損傷評価を行い、当該累積損傷評価で最も疲労損傷が大きいと評価された主応力・主ひずみ面を評価対象の主応力・主ひずみ面として、疲労寿命の評価を行う多軸疲労寿命評価方法である。
【００１２】
前記疲労寿命の評価は、前記評価対象の主応力・主ひずみ面での主応力・主ひずみ軸を基準軸とし、［数１］に示す式（１）
【００１３】
【数１】
【００１４】
により、非比例負荷係数ｆNPを求め、求めた非比例負荷係数ｆNPに基づき、下式（２ａ）または下式（２ｂ）
Ｎｆ_NP＝｛（１＋αｆNP）Δε／Ａ｝1/B ・・・（２ａ）
Ｎｆ_NP＝｛（１＋αｆNP）Δσ／Ａ｝1/B ・・・（２ｂ）
但し、Ｎｆ_NP：疲労寿命となるサイクル数
Ａ，Ｂ：材料定数
Δε：基準軸に投影して算出されるひずみ範囲
Δσ：基準軸に投影して算出される応力範囲
α：非比例負荷での追硬化の程度あるいは寿命低下の程度
を表す係数
により寿命評価を行うとよい。
【発明の効果】
【００１５】
本発明によれば、多軸負荷による疲労損傷を高精度に予測可能な多軸疲労寿命評価方法を提供できる。
【図面の簡単な説明】
【００１６】
【図１】本発明の一実施の形態に係る多軸疲労寿命評価方法の手順を示すフローチャートである。
【図２】（ａ），（ｂ）は、図１の多軸疲労寿命評価方法において、ある時刻ｔnの主応力・主ひずみ面に作用する別の時刻ｔiの主応力・主ひずみＳ（ｔi）の垂直成分について説明する図である。
【発明を実施するための形態】
【００１７】
以下、本発明の実施の形態を添付図面にしたがって説明する。
【００１８】
まず、本実施の形態に係る多軸疲労寿命評価方法で用いるＩＳ法について説明する。なお、ＩＳ法については、非特許文献１，２に詳細に記載されているため、ここでは概略のみを説明する。
【００１９】
ＩＳ法は、主応力（ひずみ）基準のクライテリアであり、下式（３）で定義される。
ΔＳＩNP＝（１＋αｆNP）ΔＳＩ ・・・（３）
【００２０】
式（３）において、ΔＳＩは評価対象となる主応力・主ひずみ面に作用する垂直成分として算出される応力（ひずみ）範囲、すなわち、評価対象の主応力・主ひずみ面での主応力・主ひずみ軸を基準軸としたとき、その基準軸に投影して算出される応力（ひずみ）範囲である。また、αは非比例負荷（応力またはひずみの主軸方向が変化する負荷）での追硬化の程度あるいは寿命低下の程度を表す係数であり、非比例負荷による時間強度の低下率からも求めることができる。αは材料定数であるため、予め材料データとして取得しておくとよい。ｆNPは、負荷経路の非比例の程度を表すパラメータであり、非比例負荷係数と呼称する。
【００２１】
非比例負荷係数ｆNPは、［数２］に示す式（１）より得られる。
【００２２】
【数２】
【００２３】
式（１）における積分記号の中の「｜ｅ1×ｅRＳＩ（ｔ）｜」は、主応力（主ひずみ）の方向変化を含む値の変化量を表しており、ｄｓは応力（ひずみ）経路、すなわち負荷経路を表している。負荷経路の全経路に沿って、主応力（ひずみ）の変化量を経路積分することにより、非比例負荷係数ｆNPを求めることができる。なお、式（１）において積分値にπ／（２ＳＩmax・Ｌpath）を乗じているのは、ｆNPを基準化するためである。軸ひずみεまたは軸応力σとせん断ひずみγまたはせん断応力τが９０度の位相差で負荷される円形負荷では、ｆNPは１となる。
【００２４】
また、ｅ1は最大主応力（最大主ひずみ）の方向を示す単位ベクトルであり、ｅRは、ある時刻ｔ（ｅ1を基準（時刻０）とした時刻）での主応力（主ひずみ）の方向を示す単位ベクトルである。よって、ｅ1とｅRの外積は、ｅ1とｅRのなす角度ξ（ｔ）を用いて、ｓｉｎξ（ｔ）で表すことができる。よって、式（１）の積分記号の中の「｜ｅ1×ｅRＳＩ（ｔ）｜」は、「ＳＩ（ｔ）｜ｓｉｎξ（ｔ）｜」と表すこともできる。
【００２５】
以下、本実施の形態に係る多軸疲労寿命評価方法を図１を用いて説明する。
【００２６】
図１に示すように、本実施の形態に係る多軸疲労寿命評価方法では、まず、ステップＳ１にて、評価対象の構造物の各時刻での応力またはひずみ状態を基に、各時刻での主応力または主ひずみを計算する。これにより、各時刻での主応力または主ひずみに対して直交する面である各時刻での主応力・主ひずみ面も得られることになる。
【００２７】
なお、ステップＳ１で計算に用いる各時刻での応力またはひずみ状態としては、例えば、疲労寿命の評価を行う部位がセンサを配置した測定点である場合は、そのセンサによる測定データを用いることができる。また、測定点でない部位の疲労寿命の評価を行う際には、応力またはひずみ状態を任意の方法で解析した解析結果を用いることができる。また、評価を行う期間は、例えば、評価対象の構造物の起動（運用開始）から停止（運用終了）までとすればよい。具体的には、例えば、評価対象の構造物が航空機である場合には、離陸から着陸までの１サイクル、評価対象の構造物がプラントである場合は、前回の定期検査から今回の定期検査までを、評価対象の期間とすればよい。
【００２８】
その後、ステップＳ２にて、ある時刻ｔnの主応力・主ひずみ面に着目し、その面に作用する他の各時刻での主応力または主ひずみの垂直成分の値をそれぞれ算出する。
【００２９】
図２（ａ）に示すように、ある時刻ｔnの主応力・主ひずみがＳ（ｔn）であり、時刻ｔnの主応力・主ひずみ面が符号２１で示される面であったとする。このとき、図２（ｂ）に示すように、別の時刻ｔiの主応力・主ひずみがＳ（ｔi）であったとすると、時刻ｔnの主応力・主ひずみ面２１に作用する別の時刻ｔiの主応力・主ひずみＳ（ｔi）の垂直成分は、図中のＳ（ｔi）’のようになる。同様にして、ある時刻ｔnの主応力・主ひずみ面２１に作用する全ての時刻での主応力・主ひずみの垂直成分を算出する。
【００３０】
その後、ステップＳ３にて、予め求めておいた単軸負荷時のＳＮ線図に基づき、時刻ｔnの主応力・主ひずみ面について累積損傷評価を行う。
【００３１】
累積損傷評価では、マイナー則や修正マイナー則を用いて評価を行うとよい。具体的には、ステップＳ２で算出した任意の時刻ｔiの主応力・主ひずみの垂直成分Ｓ（ｔi）’に対応するサイクル数Ｎiを単軸負荷時のＳＮ線図より求めると共に、そのＳ（ｔi）’が負荷された回数ｎiを求め、下式（４）
Ｄn＝Σ（ｎi／Ｎi） ・・・（４）
により、時刻ｔnの主応力・主ひずみ面での疲労損傷度Ｄnを求める。
【００３２】
その後、ステップＳ４にて、全てのｎについて計算を行ったかを判断し、ＮＯと判断されればステップＳ２に戻る。つまり、ステップＳ２，Ｓ３を全ての時刻について実行し、すべての時刻の主応力・主ひずみ面での疲労損傷度をそれぞれ求める。なお、評価対象とする時刻の間隔については、解析結果を用いる場合には解析を行った設定（ステップ）をそのまま採用すればよく、測定結果を用いる場合には、適宜な間隔に設定すればよい。
【００３３】
その後、ステップＳ５にて、ステップＳ３で求めた疲労損傷度Ｄnが最大となる主応力・主ひずみ面、すなわち、最大の疲労損傷となる主応力・主ひずみ面を選択し、選択した主応力・主ひずみ面での主応力・主ひずみ軸（つまり主応力・主ひずみ面に垂直な軸）を基準軸とし、選択した主応力・主ひずみ面を評価面に決定する。
【００３４】
その後、ステップＳ６にて、決定した基本軸と評価面を用い、ＩＳ法により疲労寿命を評価する。
【００３５】
具体的には、上述の式（１）により、非比例負荷係数ｆNPを求め、求めた非比例負荷係数ｆNPに基づき、下式（２ａ）または下式（２ｂ）
Ｎｆ_NP＝｛（１＋αｆNP）Δε／Ａ｝1/B ・・・（２ａ）
Ｎｆ_NP＝｛（１＋αｆNP）Δσ／Ａ｝1/B ・・・（２ｂ）
但し、Ｎｆ_NP：疲労寿命となるサイクル数
Ａ，Ｂ：材料定数
Δε：基準軸に投影して算出されるひずみ範囲
Δσ：基準軸に投影して算出される応力範囲
α：非比例負荷での追硬化の程度あるいは寿命低下の程度
を表す係数
により寿命評価を行う。なお、式（２ａ），（２ｂ）におけるＡ，Ｂ，αは、いずれも材料定数であるので、予め材料データとして取得しておくとよい。
【００３６】
以上説明したように、本実施の形態に係る多軸疲労寿命評価方法では、評価対象の構造物の各時刻での応力またはひずみ状態を基に、各時刻での主応力または主ひずみに対して直交する面である主応力・主ひずみ面をそれぞれ求め、求めた各時刻での主応力・主ひずみ面のそれぞれに対して、当該主応力・主ひずみ面に作用する各時刻での主応力または主ひずみの垂直成分の値をそれぞれ算出すると共に、当該算出した値と、予め求めておいた単軸負荷時のＳＮ線図とに基づき、各主応力・主ひずみ面での累積損傷評価を行い、当該累積損傷評価で最も損傷が大きいと評価された主応力・主ひずみ面を評価対象の主応力・主ひずみ面として、疲労寿命の評価を行うようにしている。
【００３７】
つまり、本実施の形態に係る多軸疲労寿命評価方法では、評価対象とする損傷面を決定するために、評価対象の応力・ひずみの時刻歴から全ての時刻について主応力・主ひずみを計算し累積損傷評価を行い、最大の疲労損傷を受ける面を評価対象としている。
【００３８】
本発明では、従来方法のように主応力・主ひずみが最大となる面を評価するのではなく、最も疲労損傷が大きい主応力・主ひずみ面を評価対象として疲労寿命の評価を行っているため、小さい負荷が多数回負荷されるような応力・ひずみ状態であっても、危険側の評価となってしまうことを抑制でき、主軸の大きさと方向が複雑に変化する多軸負荷による疲労損傷を高精度に予測することが可能になる。
【００３９】
なお、本実施の形態では、クライテリアとして主応力・主ひずみを用いたが、これに限らず、クライテリアとしてミーゼス型やトレスカ型などの相当応力（ひずみ）を採用することも勿論可能である。つまり、本発明によれば、いずれのクライテリアを採用した場合であっても、疲労損傷が最も大きくなる面を評価対象とし、高精度な疲労寿命の評価を行うことが可能になる。
【技術分野】
【０００１】
本発明は、航空機エンジン、過給器などの回転機械や、圧力容器などを含む一般構造物の疲労寿命を評価する際に用いられる多軸疲労寿命評価方法に関するものである。
【背景技術】
【０００２】
実機部材の受ける荷重は、時刻と共に変化するランダム荷重である。したがって、主応力や主ひずみの方向（主軸という）が変化することも想定される。主軸方向を変化させた多軸疲労試験結果によれば、せん断ひずみと軸ひずみの位相あるいはせん断応力と軸応力の位相がずれるに従い、疲労寿命は単軸の寿命（等価ひずみあるいは等価応力を用いて求めた疲労寿命）よりも低下していくことが報告されている。
【０００３】
こうした問題点を解決するために、負荷経路（主軸方向変化を含む負荷経路）を考慮して、主軸の変化にも対応できる多軸疲労寿命評価方法がいくつか提案されている。
【０００４】
なかでも、非特許文献１，２では、評価対象とする期間中で主応力・主ひずみが最も大きな値となる主応力・主ひずみ面に着目して、主軸の変化量を定量化する方法、および負荷経路（主軸方向変化を含む負荷経路）に基づいた多軸疲労寿命評価方法（ＩＳ（Itoh-Sakane）法）が提案されている。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】伊藤隆基、「非比例多軸低サイクル疲労寿命評価モデル」、材料、社団法人日本材料学会、２００１年１２月１５日、ｖｏｌ．５０、Ｎｏ．１２、ｐｐ．１３１７−１３２２
【非特許文献２】Takamoto Itoh、Tomohiko Ozaki、Toru Amaya、and Masao Sakane、「DETERMINATION OF STRESS AND STRAIN RANGES UNDER NON-PROPORTIONAL CYCLIC LOADING」、8th International Conference on Multiaxial Fatigue & Fracture、２００７年
【非特許文献３】Takamoto Itoh、Masao Sakane、Takahiro Hata、Naomi Hamada、「A design procedure for assessing low cycle fatigue life under proportional and non-proportional loading」、International Journal of Fatigue 28、２００６年、ｐｐ．４５９−４６６
【発明の概要】
【発明が解決しようとする課題】
【０００６】
上述のように、従来のＩＳ法では、評価対象とする期間中（例えば１サイクル中）で主応力・主ひずみが最大となる面に着目して、疲労損傷を評価してきた。つまり、従来のＩＳ法では、主応力・主ひずみが最大となる面が、最も疲労損傷が大きいとして評価を行ってきた。
【０００７】
しかしながら、実機に作用する応力・ひずみの状態によっては、主応力・主ひずみが最大となる面が、最も疲労損傷が大きいとは限らないという問題がある。
【０００８】
例えば、１サイクル中に大きな軸負荷がｎ回作用し、小さいねじり負荷がＮ回作用する場合について考えると、ｎ≒Ｎであれば、主応力・主ひずみが最大となる面が最も疲労損傷が大きい面となる。しかし、ｎ＜＜Ｎである場合、すなわち、大きい軸負荷が１回作用する間に、多数回の小さいねじり負荷が作用するような場合においては、回数の多いねじり負荷による疲労損傷が大きくなり、主応力・主ひずみが最大となる面が最も疲労損傷が大きい面とならない場合もある。
【０００９】
このような場合に、従来通り主応力・主ひずみが最大となる面を評価対象として疲労寿命の評価を行うと、最も疲労損傷が大きい面でない面を評価対象としていることとなり、危険側の評価となってしまう。
【００１０】
そこで、本発明の目的は、上記課題を解決し、多軸負荷による疲労損傷を高精度に予測可能な多軸疲労寿命評価方法を提供することにある。
【課題を解決するための手段】
【００１１】
本発明は上記目的を達成するために創案されたものであり、評価対象の構造物の各時刻での応力またはひずみ状態を基に、各時刻での主応力または主ひずみに対して直交する面である主応力・主ひずみ面をそれぞれ求め、求めた各時刻での主応力・主ひずみ面のそれぞれに対して、当該主応力・主ひずみ面に作用する各時刻での主応力または主ひずみの垂直成分の値をそれぞれ算出し、当該算出した値と、予め求めておいた単軸負荷時のＳＮ線図とに基づき、各主応力・主ひずみ面での累積損傷評価を行い、当該累積損傷評価で最も疲労損傷が大きいと評価された主応力・主ひずみ面を評価対象の主応力・主ひずみ面として、疲労寿命の評価を行う多軸疲労寿命評価方法である。
【００１２】
前記疲労寿命の評価は、前記評価対象の主応力・主ひずみ面での主応力・主ひずみ軸を基準軸とし、［数１］に示す式（１）
【００１３】
【数１】
【００１４】
により、非比例負荷係数ｆNPを求め、求めた非比例負荷係数ｆNPに基づき、下式（２ａ）または下式（２ｂ）
Ｎｆ_NP＝｛（１＋αｆNP）Δε／Ａ｝1/B ・・・（２ａ）
Ｎｆ_NP＝｛（１＋αｆNP）Δσ／Ａ｝1/B ・・・（２ｂ）
但し、Ｎｆ_NP：疲労寿命となるサイクル数
Ａ，Ｂ：材料定数
Δε：基準軸に投影して算出されるひずみ範囲
Δσ：基準軸に投影して算出される応力範囲
α：非比例負荷での追硬化の程度あるいは寿命低下の程度
を表す係数
により寿命評価を行うとよい。
【発明の効果】
【００１５】
本発明によれば、多軸負荷による疲労損傷を高精度に予測可能な多軸疲労寿命評価方法を提供できる。
【図面の簡単な説明】
【００１６】
【図１】本発明の一実施の形態に係る多軸疲労寿命評価方法の手順を示すフローチャートである。
【図２】（ａ），（ｂ）は、図１の多軸疲労寿命評価方法において、ある時刻ｔnの主応力・主ひずみ面に作用する別の時刻ｔiの主応力・主ひずみＳ（ｔi）の垂直成分について説明する図である。
【発明を実施するための形態】
【００１７】
以下、本発明の実施の形態を添付図面にしたがって説明する。
【００１８】
まず、本実施の形態に係る多軸疲労寿命評価方法で用いるＩＳ法について説明する。なお、ＩＳ法については、非特許文献１，２に詳細に記載されているため、ここでは概略のみを説明する。
【００１９】
ＩＳ法は、主応力（ひずみ）基準のクライテリアであり、下式（３）で定義される。
ΔＳＩNP＝（１＋αｆNP）ΔＳＩ ・・・（３）
【００２０】
式（３）において、ΔＳＩは評価対象となる主応力・主ひずみ面に作用する垂直成分として算出される応力（ひずみ）範囲、すなわち、評価対象の主応力・主ひずみ面での主応力・主ひずみ軸を基準軸としたとき、その基準軸に投影して算出される応力（ひずみ）範囲である。また、αは非比例負荷（応力またはひずみの主軸方向が変化する負荷）での追硬化の程度あるいは寿命低下の程度を表す係数であり、非比例負荷による時間強度の低下率からも求めることができる。αは材料定数であるため、予め材料データとして取得しておくとよい。ｆNPは、負荷経路の非比例の程度を表すパラメータであり、非比例負荷係数と呼称する。
【００２１】
非比例負荷係数ｆNPは、［数２］に示す式（１）より得られる。
【００２２】
【数２】
【００２３】
式（１）における積分記号の中の「｜ｅ1×ｅRＳＩ（ｔ）｜」は、主応力（主ひずみ）の方向変化を含む値の変化量を表しており、ｄｓは応力（ひずみ）経路、すなわち負荷経路を表している。負荷経路の全経路に沿って、主応力（ひずみ）の変化量を経路積分することにより、非比例負荷係数ｆNPを求めることができる。なお、式（１）において積分値にπ／（２ＳＩmax・Ｌpath）を乗じているのは、ｆNPを基準化するためである。軸ひずみεまたは軸応力σとせん断ひずみγまたはせん断応力τが９０度の位相差で負荷される円形負荷では、ｆNPは１となる。
【００２４】
また、ｅ1は最大主応力（最大主ひずみ）の方向を示す単位ベクトルであり、ｅRは、ある時刻ｔ（ｅ1を基準（時刻０）とした時刻）での主応力（主ひずみ）の方向を示す単位ベクトルである。よって、ｅ1とｅRの外積は、ｅ1とｅRのなす角度ξ（ｔ）を用いて、ｓｉｎξ（ｔ）で表すことができる。よって、式（１）の積分記号の中の「｜ｅ1×ｅRＳＩ（ｔ）｜」は、「ＳＩ（ｔ）｜ｓｉｎξ（ｔ）｜」と表すこともできる。
【００２５】
以下、本実施の形態に係る多軸疲労寿命評価方法を図１を用いて説明する。
【００２６】
図１に示すように、本実施の形態に係る多軸疲労寿命評価方法では、まず、ステップＳ１にて、評価対象の構造物の各時刻での応力またはひずみ状態を基に、各時刻での主応力または主ひずみを計算する。これにより、各時刻での主応力または主ひずみに対して直交する面である各時刻での主応力・主ひずみ面も得られることになる。
【００２７】
なお、ステップＳ１で計算に用いる各時刻での応力またはひずみ状態としては、例えば、疲労寿命の評価を行う部位がセンサを配置した測定点である場合は、そのセンサによる測定データを用いることができる。また、測定点でない部位の疲労寿命の評価を行う際には、応力またはひずみ状態を任意の方法で解析した解析結果を用いることができる。また、評価を行う期間は、例えば、評価対象の構造物の起動（運用開始）から停止（運用終了）までとすればよい。具体的には、例えば、評価対象の構造物が航空機である場合には、離陸から着陸までの１サイクル、評価対象の構造物がプラントである場合は、前回の定期検査から今回の定期検査までを、評価対象の期間とすればよい。
【００２８】
その後、ステップＳ２にて、ある時刻ｔnの主応力・主ひずみ面に着目し、その面に作用する他の各時刻での主応力または主ひずみの垂直成分の値をそれぞれ算出する。
【００２９】
図２（ａ）に示すように、ある時刻ｔnの主応力・主ひずみがＳ（ｔn）であり、時刻ｔnの主応力・主ひずみ面が符号２１で示される面であったとする。このとき、図２（ｂ）に示すように、別の時刻ｔiの主応力・主ひずみがＳ（ｔi）であったとすると、時刻ｔnの主応力・主ひずみ面２１に作用する別の時刻ｔiの主応力・主ひずみＳ（ｔi）の垂直成分は、図中のＳ（ｔi）’のようになる。同様にして、ある時刻ｔnの主応力・主ひずみ面２１に作用する全ての時刻での主応力・主ひずみの垂直成分を算出する。
【００３０】
その後、ステップＳ３にて、予め求めておいた単軸負荷時のＳＮ線図に基づき、時刻ｔnの主応力・主ひずみ面について累積損傷評価を行う。
【００３１】
累積損傷評価では、マイナー則や修正マイナー則を用いて評価を行うとよい。具体的には、ステップＳ２で算出した任意の時刻ｔiの主応力・主ひずみの垂直成分Ｓ（ｔi）’に対応するサイクル数Ｎiを単軸負荷時のＳＮ線図より求めると共に、そのＳ（ｔi）’が負荷された回数ｎiを求め、下式（４）
Ｄn＝Σ（ｎi／Ｎi） ・・・（４）
により、時刻ｔnの主応力・主ひずみ面での疲労損傷度Ｄnを求める。
【００３２】
その後、ステップＳ４にて、全てのｎについて計算を行ったかを判断し、ＮＯと判断されればステップＳ２に戻る。つまり、ステップＳ２，Ｓ３を全ての時刻について実行し、すべての時刻の主応力・主ひずみ面での疲労損傷度をそれぞれ求める。なお、評価対象とする時刻の間隔については、解析結果を用いる場合には解析を行った設定（ステップ）をそのまま採用すればよく、測定結果を用いる場合には、適宜な間隔に設定すればよい。
【００３３】
その後、ステップＳ５にて、ステップＳ３で求めた疲労損傷度Ｄnが最大となる主応力・主ひずみ面、すなわち、最大の疲労損傷となる主応力・主ひずみ面を選択し、選択した主応力・主ひずみ面での主応力・主ひずみ軸（つまり主応力・主ひずみ面に垂直な軸）を基準軸とし、選択した主応力・主ひずみ面を評価面に決定する。
【００３４】
その後、ステップＳ６にて、決定した基本軸と評価面を用い、ＩＳ法により疲労寿命を評価する。
【００３５】
具体的には、上述の式（１）により、非比例負荷係数ｆNPを求め、求めた非比例負荷係数ｆNPに基づき、下式（２ａ）または下式（２ｂ）
Ｎｆ_NP＝｛（１＋αｆNP）Δε／Ａ｝1/B ・・・（２ａ）
Ｎｆ_NP＝｛（１＋αｆNP）Δσ／Ａ｝1/B ・・・（２ｂ）
但し、Ｎｆ_NP：疲労寿命となるサイクル数
Ａ，Ｂ：材料定数
Δε：基準軸に投影して算出されるひずみ範囲
Δσ：基準軸に投影して算出される応力範囲
α：非比例負荷での追硬化の程度あるいは寿命低下の程度
を表す係数
により寿命評価を行う。なお、式（２ａ），（２ｂ）におけるＡ，Ｂ，αは、いずれも材料定数であるので、予め材料データとして取得しておくとよい。
【００３６】
以上説明したように、本実施の形態に係る多軸疲労寿命評価方法では、評価対象の構造物の各時刻での応力またはひずみ状態を基に、各時刻での主応力または主ひずみに対して直交する面である主応力・主ひずみ面をそれぞれ求め、求めた各時刻での主応力・主ひずみ面のそれぞれに対して、当該主応力・主ひずみ面に作用する各時刻での主応力または主ひずみの垂直成分の値をそれぞれ算出すると共に、当該算出した値と、予め求めておいた単軸負荷時のＳＮ線図とに基づき、各主応力・主ひずみ面での累積損傷評価を行い、当該累積損傷評価で最も損傷が大きいと評価された主応力・主ひずみ面を評価対象の主応力・主ひずみ面として、疲労寿命の評価を行うようにしている。
【００３７】
つまり、本実施の形態に係る多軸疲労寿命評価方法では、評価対象とする損傷面を決定するために、評価対象の応力・ひずみの時刻歴から全ての時刻について主応力・主ひずみを計算し累積損傷評価を行い、最大の疲労損傷を受ける面を評価対象としている。
【００３８】
本発明では、従来方法のように主応力・主ひずみが最大となる面を評価するのではなく、最も疲労損傷が大きい主応力・主ひずみ面を評価対象として疲労寿命の評価を行っているため、小さい負荷が多数回負荷されるような応力・ひずみ状態であっても、危険側の評価となってしまうことを抑制でき、主軸の大きさと方向が複雑に変化する多軸負荷による疲労損傷を高精度に予測することが可能になる。
【００３９】
なお、本実施の形態では、クライテリアとして主応力・主ひずみを用いたが、これに限らず、クライテリアとしてミーゼス型やトレスカ型などの相当応力（ひずみ）を採用することも勿論可能である。つまり、本発明によれば、いずれのクライテリアを採用した場合であっても、疲労損傷が最も大きくなる面を評価対象とし、高精度な疲労寿命の評価を行うことが可能になる。
【特許請求の範囲】
【請求項１】
評価対象の構造物の各時刻での応力またはひずみ状態を基に、各時刻での主応力または主ひずみに対して直交する面である主応力・主ひずみ面をそれぞれ求め、
求めた各時刻での主応力・主ひずみ面のそれぞれに対して、当該主応力・主ひずみ面に作用する各時刻での主応力または主ひずみの垂直成分の値をそれぞれ算出すると共に、
当該算出した値と、予め求めておいた単軸負荷時のＳＮ線図とに基づき、各主応力・主ひずみ面での累積損傷評価を行い、
当該累積損傷評価で最も損傷が大きいと評価された主応力・主ひずみ面を評価対象の主応力・主ひずみ面として、疲労寿命の評価を行う
ことを特徴とする多軸疲労寿命評価方法。
【請求項２】
前記疲労寿命の評価は、
前記評価対象の主応力・主ひずみ面での主応力・主ひずみ軸を基準軸とし、［数１］に示す式（１）
【数１】
により、非比例負荷係数ｆNPを求め、
求めた非比例負荷係数ｆNPに基づき、下式（２ａ）または下式（２ｂ）
Ｎｆ_NP＝｛（１＋αｆNP）Δε／Ａ｝1/B ・・・（２ａ）
Ｎｆ_NP＝｛（１＋αｆNP）Δσ／Ａ｝1/B ・・・（２ｂ）
但し、Ｎｆ_NP：疲労寿命となるサイクル数
Ａ，Ｂ：材料定数
Δε：基準軸に投影して算出されるひずみ範囲
Δσ：基準軸に投影して算出される応力範囲
α：非比例負荷での追硬化の程度あるいは寿命低下の程度
を表す係数
により寿命評価を行う
請求項１記載の多軸疲労寿命評価方法。
【請求項１】
評価対象の構造物の各時刻での応力またはひずみ状態を基に、各時刻での主応力または主ひずみに対して直交する面である主応力・主ひずみ面をそれぞれ求め、
求めた各時刻での主応力・主ひずみ面のそれぞれに対して、当該主応力・主ひずみ面に作用する各時刻での主応力または主ひずみの垂直成分の値をそれぞれ算出すると共に、
当該算出した値と、予め求めておいた単軸負荷時のＳＮ線図とに基づき、各主応力・主ひずみ面での累積損傷評価を行い、
当該累積損傷評価で最も損傷が大きいと評価された主応力・主ひずみ面を評価対象の主応力・主ひずみ面として、疲労寿命の評価を行う
ことを特徴とする多軸疲労寿命評価方法。
【請求項２】
前記疲労寿命の評価は、
前記評価対象の主応力・主ひずみ面での主応力・主ひずみ軸を基準軸とし、［数１］に示す式（１）
【数１】
により、非比例負荷係数ｆNPを求め、
求めた非比例負荷係数ｆNPに基づき、下式（２ａ）または下式（２ｂ）
Ｎｆ_NP＝｛（１＋αｆNP）Δε／Ａ｝1/B ・・・（２ａ）
Ｎｆ_NP＝｛（１＋αｆNP）Δσ／Ａ｝1/B ・・・（２ｂ）
但し、Ｎｆ_NP：疲労寿命となるサイクル数
Ａ，Ｂ：材料定数
Δε：基準軸に投影して算出されるひずみ範囲
Δσ：基準軸に投影して算出される応力範囲
α：非比例負荷での追硬化の程度あるいは寿命低下の程度
を表す係数
により寿命評価を行う
請求項１記載の多軸疲労寿命評価方法。
【図１】
【図２】
【図２】
【公開番号】特開２０１３−４４６６６（Ｐ２０１３−４４６６６Ａ）
【公開日】平成２５年３月４日（２０１３．３．４）
【国際特許分類】
【出願番号】特願２０１１−１８３６３７（Ｐ２０１１−１８３６３７）
【出願日】平成２３年８月２５日（２０１１．８．２５）
【出願人】（０００００００９９）株式会社ＩＨＩ (5,014)
【出願人】（５９３００６６３０）学校法人立命館 (359)
【出願人】（５０４１４５３２０）国立大学法人福井大学 (287)
【Ｆターム（参考）】
【公開日】平成２５年３月４日（２０１３．３．４）
【国際特許分類】
【出願日】平成２３年８月２５日（２０１１．８．２５）
【出願人】（０００００００９９）株式会社ＩＨＩ (5,014)
【出願人】（５９３００６６３０）学校法人立命館 (359)
【出願人】（５０４１４５３２０）国立大学法人福井大学 (287)
【Ｆターム（参考）】
[ Back to top ]