(57)【要約】【課題】効率よく短時間でウォーカー指数を取得可能な平均応力評価パラメータの算出方法を提供する。【解決手段】平均応力としてウォーカー応力σwを用い、そのウォーカー応力σwから疲労寿命などの疲労強度を評価するに際して、前記ウォーカー応力σwのパラメータとなるウォーカー指数ｍを算出する方法において、応力比Ｒ、疲労寿命、および最大応力σmaxを制御因子とし、かつ温度を信号因子として、実験計画法により最適疲労試験条件を選定し、その最適疲労試験条件で疲労試験を行い、その疲労試験の結果から前記ウォーカー指数ｍを算出する。【選択図】図１
2010112783.tif 000002
【特許請求の範囲】
【請求項１】
平均応力として、下式（１）
σw＝２σa（１−Ｒ）m-1
＝σmax（１−Ｒ）m              …（１）
但し、σa：応力振幅
σmax：最大応力
Ｒ：応力比
ｍ：ウォーカー指数
で表されるウォーカー応力σwを用い、前記ウォーカー応力σwのパラメータとなるウォーカー指数ｍを算出する平均応力評価パラメータの算出方法において、
応力比Ｒ、疲労寿命、および最大応力σmaxを制御因子とし、かつ温度を信号因子として、実験計画法により最適疲労試験条件を選定し、その最適疲労試験条件で疲労試験を行い、その疲労試験の結果から前記ウォーカー指数ｍを算出することを特徴とする平均応力評価パラメータの算出方法。
【請求項２】
前記制御因子および前記信号因子の水準をそれぞれ設定すると共に、疲労強度を評価する材料の疲労強度線図を予め取得しておき、その疲労強度線図を用いて、前記制御因子の各水準と前記信号因子の各水準の組み合わせごとに前記ウォーカー指数ｍを算出し、前記制御因子の各水準、前記信号因子の各水準、および算出した前記ウォーカー指数ｍを直交表に割付けて分散分析を行って、前記制御因子の各水準のＳＮ比を算出し、前記制御因子の前記ＳＮ比が最大となる水準を組み合わせたものを前記最適疲労試験条件として選定する請求項１記載の平均応力評価パラメータの算出方法。
【請求項３】
前記最適疲労試験条件で疲労試験を行って疲労強度線図を求め、その疲労強度線図および式（１）を用いて前記ウォーカー指数ｍを算出する請求項１または２記載の平均応力評価パラメータの算出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、航空機エンジンのシャフト、ディスク、ブレードなどの疲労強度を評価するに際して、平均応力を求めるために用いる平均応力評価パラメータの算出方法に係り、特に、平均応力としてウォーカー応力σwを用い、そのパラメータであるウォーカー指数ｍを算出する方法に関するものである。
【背景技術】
【０００２】
航空機エンジンでは、吸い込み空気がファンにおいて加圧され、その空気の一部がコアエンジンに導かれ、圧縮機で加圧された後、燃焼器内で燃料と混合され、点火されて、高圧タービン、低圧タービンに順次排出される。高圧タービンは圧縮機を駆動し、低圧タービンはファンを駆動する。
【０００３】
航空機エンジンのファンやタービンに用いられるシャフト、ディスク、ブレードなどの部材には、圧縮空気や燃焼空気による応力が常にかかるため、金属疲労によりき裂が発生する場合がある。そのため、これらの部材の疲労寿命などの疲労強度を予め求めておき、き裂が発生する前に取り替える必要がある。
【０００４】
航空機エンジンのシャフト、ディスク、ブレードなどの疲労強度を評価するに際して、平均応力が疲労特性（疲労寿命など）に影響を及ぼすことが知られている。
【０００５】
図３に示すように、応力比Ｒ＞０である場合、応力比Ｒが増えるにつれ平均応力も増加し、疲労寿命が短くなり、疲労強度が低下する傾向がある。応力比Ｒは、最大応力をσmax、最小応力をσminとすると、Ｒ＝σmin／σmaxで表される。
【０００６】
一方、応力比Ｒ＜０である場合は、σminが負となり材料に圧縮の荷重が付加されることになる。この場合、疲労寿命が延び、疲労強度が増加する傾向がある。
【０００７】
応力比Ｒおよび平均応力が疲労特性に及ぼす影響を評価する方法として、ウォーカー（Walker）応力（あるいはひずみ）がある。ウォーカー応力σwは、下式（１）
σw＝２σa（１−Ｒ）m-1
＝σmax（１−Ｒ）m              …（１）
但し、σa：応力振幅
σmax：最大応力
Ｒ：応力比
ｍ：ウォーカー指数
で表される。式（１）において、応力振幅σaは片振幅であり、σa＝（σmax−σmin）／２である。
【０００８】
なお、この出願の発明に関連する先行技術文献情報としては、次のものがある。
【０００９】
【非特許文献１】Ｗａｌｋｅｒ，Ｋ．、「Ｔｈｅ  Ｅｆｆｅｃｔ  ｏｆ  Ｓｔｒｅｓｓ  Ｒａｔｉｏ  Ｄｕｒｉｎｇ  Ｃｒａｃｋ  Ｐｒｏｐａｇａｔｉｏｎ  ａｎｄ  Ｆａｔｉｇｕｅ  ｆｏｒ  ２０２４−Ｔ３  ａｎｄ  ７０７５−Ｔ６  Ａｌｕｍｉｎｕｍ」、ＡＳＴＭ  ＳＴＰ  ４６２、Ａｍｅｒｉｃａｎ  Ｓｏｃｉｅｔｙ  ｆｏｒ  Ｔｅｓｔｉｎｇ  ａｎｄ  Ｍａｔｅｒｉａｌｓ、１９７０年、ｐｐ．１−１４
【発明の開示】
【発明が解決しようとする課題】
【００１０】
疲労強度を評価する際の平均応力としてウォーカー応力σwを用いる場合、そのパラメータとなるウォーカー指数ｍを求める必要がある。
【００１１】
しかしながら、ウォーカー指数ｍは材料、試験片形状、あるいは温度などの疲労試験条件の影響を受けるため、ウォーカー指数ｍを精度よく求めるためには、試験条件を変えて数多くの疲労試験を実施する必要がある。そのため、疲労試験数が非常に多くなり、莫大な時間とコストがかかってしまうという問題があった。
【００１２】
そこで、本発明の目的は、上記課題を解決し、効率よく短時間でウォーカー指数を取得可能な平均応力評価パラメータの算出方法を提供することにある。
【課題を解決するための手段】
【００１３】
本発明は上記目的を達成するために創案されたものであり、請求項１の発明は、平均応力として、下式（１）
σw＝２σa（１−Ｒ）m-1
＝σmax（１−Ｒ）m              …（１）
但し、σa：応力振幅
σmax：最大応力
Ｒ：応力比
ｍ：ウォーカー指数
で表されるウォーカー応力σwを用い、前記ウォーカー応力σwのパラメータとなるウォーカー指数ｍを算出する平均応力評価パラメータの算出方法において、応力比Ｒ、疲労寿命、および最大応力σmaxを制御因子とし、かつ温度を信号因子として、実験計画法により最適疲労試験条件を選定し、その最適疲労試験条件で疲労試験を行い、その疲労試験の結果から前記ウォーカー指数ｍを算出する平均応力評価パラメータの算出方法である。
【００１４】
請求項２の発明は、前記制御因子および前記信号因子の水準をそれぞれ設定すると共に、疲労強度を評価する材料の疲労強度線図を予め取得しておき、その疲労強度線図を用いて、前記制御因子の各水準と前記信号因子の各水準の組み合わせごとに前記ウォーカー指数ｍを算出し、前記制御因子の各水準、前記信号因子の各水準、および算出した前記ウォーカー指数ｍを直交表に割付けて分散分析を行って、前記制御因子の各水準のＳＮ比を算出し、前記制御因子の前記ＳＮ比が最大となる水準を組み合わせたものを前記最適疲労試験条件として選定する請求項１記載の平均応力評価パラメータの算出方法である。
【００１５】
請求項３の発明は、前記最適疲労試験条件で疲労試験を行って疲労強度線図を求め、その疲労強度線図および式（１）を用いて前記ウォーカー指数ｍを算出する請求項１または２記載の平均応力評価パラメータの算出方法である。
【発明の効果】
【００１６】
本発明によれば、最適な疲労試験条件を選定することができるので、疲労強度評価に用いるウォーカー指数を取得するための疲労試験を効率よく、かつ短時間で行うことが可能となり、コストを抑制できる。
【発明を実施するための最良の形態】
【００１７】
以下、本発明の好適な実施の形態を添付図面にしたがって説明する。
【００１８】
本発明の平均応力評価パラメータの算出方法は、例えば、航空機エンジンのシャフト、ディスク、ブレードなどの部材の疲労強度を評価するに際して、平均応力を求めるために用いるウォーカー指数を算出する方法である。すなわち、本発明における平均応力評価パラメータは、ウォーカー指数である。
【００１９】
本発明では、疲労強度を評価する際の平均応力として、下式（１）
σw＝２σa（１−Ｒ）m-1
＝σmax（１−Ｒ）m              …（１）
但し、σa：応力振幅
σmax：最大応力
Ｒ：応力比
ｍ：ウォーカー指数
で表されるウォーカー応力σwを用いる。
【００２０】
ウォーカー応力σwのパラメータとなるウォーカー指数ｍは、材料、試験片形状、あるいは温度などの疲労試験条件の影響を受けるため、ウォーカー指数ｍを精度よく求めるためには、数多くの疲労試験を実施する必要があった。
【００２１】
そこで、本発明者らは、ウォーカー指数ｍを効率よく、かつ精度よく算出できる方法について鋭意検討を行い、その結果、実験計画法を用いて予め最適疲労試験条件を選定しておき、その最適疲労試験条件で疲労試験を行ってウォーカー指数ｍを求めれば、ウォーカー指数ｍを効率よく短時間で、かつ精度よく算出できることを見出した。
【００２２】
図１は、本実施形態においてウォーカー指数ｍを算出する際のフローチャートである。
【００２３】
図１に示すように、ウォーカー指数ｍを算出する際には、まず、制御因子および信号因子を設定する（ステップＳ１）。
【００２４】
本発明では、最適な疲労試験条件を選定するため、制御因子および信号因子として疲労試験条件を採用する。具体的には、制御因子Ａとしてノッチの有無、制御因子Ｂとして応力比Ｒ、制御因子Ｃとして疲労寿命、制御因子Ｄとして最大応力σmaxを採用し、信号因子として温度を採用する。
【００２５】
また、各制御因子および信号因子の水準を設定する。本実施形態で設定した各制御因子の水準を表１に、信号因子の水準を表２に示す。
【００２６】
【表１】
2010112783.tif 000003
【００２７】
【表２】
2010112783.tif 000004
【００２８】
表１において、σyは降伏応力である。また、表２において、信号因子の水準αはＲＴ（Room Temperature；室温）であり、６８Ｆ（３７．８℃）とした。
【００２９】
さらに、疲労強度を評価する材料の疲労強度線図を予め取得しておく。本実施形態では、疲労強度を評価する材料として、ニッケル基超合金であるインコネル（登録商標）７１８を用いた場合を説明する。材料としてインコネル７１８を用いた場合の疲労強度線図を図２に示す。
【００３０】
その後、ステップＳ１で設定した制御因子の水準と信号因子の水準の組み合わせごとにウォーカー指数ｍを算出し、これらを直交表に割付ける（ステップＳ２，Ｓ３）。
【００３１】
本実施形態で用いる直交表を表３に示す。
【００３２】
【表３】
2010112783.tif 000005
【００３３】
表３において、制御因子の列にある数字（１，２，３）は、各制御因子の水準を表し、（α，β，γ）は信号因子の水準を表す。
【００３４】
表１、２で設定した制御因子および信号因子の各水準を表３の直交表に割り当てた一覧を表４に示す。
【００３５】
【表４】
2010112783.tif 000006
【００３６】
表４において、制御因子Ｄ（最大応力）については、第三水準が存在しないため、第三水準を第一水準に置換している。
【００３７】
表４に示すように、例えば、実験番号１では、以下の疲労試験条件が割り当てられる。
【００３８】
制御因子Ａ（ノッチの有無）：Unnotched；試験片の切り欠きがない
制御因子Ｂ（応力比）：All data；図２のデータを用いる場合、応力比Ｒ＝−１．００，−０．５０，０．１０，０．２０，０．５０
制御因子Ｃ（疲労寿命）：All data：図２中の全疲労寿命域（１０4〜１０8）
制御因子Ｄ（最大応力）：All data：図２中の全ての応力条件
これらの条件に合致する疲労試験結果を図２から抽出し、ウォーカー指数ｍを計算する。実験番号１の場合、ＲＴ（室温）でのウォーカー指数ｍは、ｍ＝０．５４となる。疲労強度線図からウォーカー指数ｍを算出する方法については、公知であるため説明を省略するが、例えば、最小二乗法などを用いて算出するとよい。
【００３９】
同様にして、制御因子の水準と信号因子の水準の組み合わせごとにウォーカー指数ｍを算出して、直交表に割り付ける。ウォーカー指数ｍの計算結果を表４に併せて示す。
【００４０】
その後、ステップＳ３で作成した直交表を用いて分散分析を行い、制御因子の各水準のＳＮ比を算出する（ステップＳ４）。
【００４１】
以下、具体例として、制御因子Ａ（ノッチの有無）の第一水準（Unnotched）のＳＮ比ηを算出する手順を説明する。
【００４２】
表４から抜粋した制御因子Ａの第一水準におけるウォーカー指数ｍの計算結果を表５に示す。
【００４３】
【表５】
2010112783.tif 000007
【００４４】
また、説明を簡略化するため、表５の数値を記号化したものを表６に示す。
【００４５】
【表６】
2010112783.tif 000008
【００４６】
以下の計算式では、表６の記号を用いて説明する。
【００４７】
まず、表５のデータから、修正項Ｓmを［数１］に示す式（２）
【００４８】
【数１】
2010112783.tif 000009
【００４９】
より求める。式（２）において、実験数ｎは制御因子Ａを第一水準として設定した実験の数であり、ここでは、ｎ＝９である。また、信号数ｋは信号因子の水準の数であり、ここでは、ｋ＝２である。
【００５０】
表５の値を式（２）に代入すると、修正項Ｓmは、
Ｓm＝（０．５４＋０．６１＋…＋０．５９＋０．５５）／（９×２）
＝５．４１
となる。
【００５１】
また、データの二乗和Ｓtを下式（３）
Ｓt＝（ｙ11）2＋（ｙ12）2＋…＋（ｙ91）2＋（ｙ92）2        …（３）
より求める。表５の値を式（３）に代入すると、データの二乗和Ｓtは、
Ｓt＝０．５４2＋０．６１2＋…＋０．５９2＋０．５５2
＝５．４８
となる。
【００５２】
さらに、信号因子の平均Ｍを、下式（４）
Ｍ＝（Ｍ１＋Ｍ２）／ｋ                              …（４）
より求める。表５の値を式（４）に代入すると、信号因子の平均Ｍは、
Ｍ＝（６８＋１０００）／２
＝５３４
となる。
【００５３】
その後、有効除数ｒを、下式（５）
ｒ＝ｎ［（Ｍ1−Ｍ）2＋（Ｍ2−Ｍ）2］                  …（５）
より求める。表５の値および式（４）で求めた信号因子の平均Ｍの値を式（５）に代入すると、有効除数ｒは、
ｒ＝９［（６８−５３４）2＋（１０００−５３４）2］
＝３９０８８０８
となる。
【００５４】
また、信号因子の一次回帰の効果Ｓβを、［数２］に示す式（６）
【００５５】
【数２】
2010112783.tif 000010
【００５６】
より求める。表５の値と、式（４）で求めた信号因子の平均Ｍの値と、式（５）で求めた有効除数ｒの値とを式（６）に代入すると、信号因子の一次回帰の効果Ｓβは、
Ｓβ＝［（６８−５３４）×４．５５＋（１０００−５３４）×５．３２］／３９０８８０８
＝０．０３３
となる。
【００５７】
さらに、誤差平均の二乗和ＳNを、［数３］に示す式（７）
【００５８】
【数３】
2010112783.tif 000011
【００５９】
より求める。表５の値および式（２）で求めた修正項Ｓmの値を式（７）に代入すると、誤差平均の二乗和ＳNは、
ＳN＝（１．４８2＋１．４５2＋…＋１．１２2）／２−５．４１
＝０．０２１
となる。
【００６０】
また、誤差変動Ｓeを下式（８）
Ｓe＝Ｓt−Ｓm−Ｓβ−ＳN                …（８）
より算出する。式（２）で求めた修正項Ｓmの値と、式（３）で求めたデータの二乗和Ｓtの値と、式（６）で求めた信号因子の一次回帰の効果Ｓβの値と、式（７）で求めた誤差平均の二乗和ＳNの値とを式（８）に代入すると、誤差変動Ｓeは、
Ｓe＝５．４８−５．４１−０．０３３−０．０２１
＝０．０２３
となる。
【００６１】
その後、誤差分散Ｖeを、［数４］に示す式（９）
【００６２】
【数４】
2010112783.tif 000012
【００６３】
より求める。式（８）で求めた誤差変動Ｓeの値と、表５の値とを式（９）に代入すると、誤差分散Ｖeは、
Ｖe＝０．０２３／（９×２−１−９）
＝０．００２９
となる。
【００６４】
最後に、ＳＮ比ηを［数５］に示す式（１０）
【００６５】
【数５】
2010112783.tif 000013
【００６６】
より求める。式（５）で求めた有効除数ｒの値と、式（６）で求めた信号因子の一次回帰の効果Ｓβの値と、式（９）で求めた誤差分散Ｖeの値とを式（１０）に代入すると、ＳＮ比ηは、
η＝１０ｌｏｇ10［０．０３３／（３９０８８０８×０．００２９）］
＝−５５．２７
となる。以上により、制御因子Ａの第一水準のＳＮ比ηが求められる。
【００６７】
同様にして、制御因子Ａの第二水準、および制御因子Ｂ〜Ｄの各水準のＳＮ比ηを算出する。算出した各制御因子Ａ〜ＤのＳＮ比ηの一覧を表７に示す。
【００６８】
【表７】
2010112783.tif 000014
【００６９】
ステップＳ４で各制御因子Ａ〜ＤのＳＮ比ηを算出した後、制御因子Ａ〜Ｄの各水準のＳＮ比ηを比較し、制御因子のＳＮ比ηが最大となる水準を組み合わせたものを最適疲労試験条件として選定する（ステップＳ５）。
【００７０】
表７より、制御因子Ａ（ノッチの有無）では第一水準（Unnotched）、制御因子Ｂ（応力比）では第二水準（Ｒ＝−０．５，０．１，０．５）、制御因子Ｃ（疲労寿命）では第二水準（Ｎ＜１０7）、制御因子Ｄ（最大応力）では第一水準（All data）のＳＮ比ηが最大となり、これらを組み合わせたものが最適疲労試験条件として選定される。
【００７１】
その後、選定した最適疲労試験条件で実際に疲労試験を行って疲労強度線図を求め、その疲労強度線図を用いて、ウォーカー指数ｍを算出する（ステップＳ６）。
【００７２】
以上により、ウォーカー応力σwのパラメータであるウォーカー指数ｍが得られる。
【００７３】
算出したウォーカー指数ｍを用いてウォーカー応力σwを求めることで、ウォーカー応力σwを用いて強度評価を行うことが可能となる。
【００７４】
以上説明したように、本実施形態では、応力比Ｒ、疲労寿命、および最大応力σmaxを制御因子とし、かつ温度を信号因子として、実験計画法により最適疲労試験条件を選定し、その最適疲労試験条件で疲労試験を行い、その疲労試験の結果からウォーカー指数ｍを算出している。
【００７５】
これにより、最適疲労試験条件を選定できるため、無駄な疲労試験を行う必要がなくなり、効率よく、かつ短時間で精度よくウォーカー指数ｍを算出することが可能となり、コストを抑制することができる。
【００７６】
上記実施形態では、制御因子としてノッチの有無、応力比Ｒ、疲労寿命および最大応力σmaxを用い、信号因子として温度を用いたが、制御因子および信号因子はこれらに限定されない。
【図面の簡単な説明】
【００７７】
【図１】本発明においてウォーカー指数を算出する際のフローチャートである。
【図２】本発明においてウォーカー指数を算出する際に用いた、疲労強度を評価する材料の疲労強度線図である。
【図３】疲労強度が応力比および平均応力の影響を受けることを説明する図である。
【図１】
2010112783.tif 000015
【図２】
2010112783.tif 000016
【図３】
2010112783.tif 000017