(57)【要約】      （修正有）【課題】使用後の軸受から、ラジアル荷重、アキシアル荷重、および回転回数を合理的に推定することのできる使用条件推定方法を提案する。【解決手段】使用後の転がり軸受の回転輪と固定輪に対してＸ線分析を行い、このＸ線分析の結果から回転輪の最大転動体荷重の推定値を得ると共に、固定輪の任意の１点以上の位置での転動体荷重の推定値を得る（Ｓ１）。得られた回転輪の最大転動体荷重の推定値と固定輪の任意の１点以上の位置での転動体荷重の推定値とから軸受の負荷分布を推定する（Ｓ２）。この負荷分布と軸受における転動体と内外輪との接触角とから、ラジアル荷重とアキシアル荷重とを推定する（Ｓ３）。負荷分布と、負荷回数Ｎ、繰り返し応力Ｓ、Ｘ線分析で求まる半価幅ｗ(°）の関係を求めておいた結果とから、使用された回転回数を推定する（Ｓ４）。【選択図】図１
2011069684.tif 000002
【特許請求の範囲】
【請求項１】
使用後の転がり軸受の回転輪と固定輪に対してＸ線分析を行い、このＸ線分析の結果から回転輪の最大転動体荷重の推定値を得ると共に、固定輪の任意の１点以上の位置での転動体荷重の推定値を得る過程と、
この過程で得られた回転輪の最大転動体荷重の推定値と固定輪の任意の１点以上の位置での転動体荷重の推定値とから軸受の負荷分布を推定する過程と、
この推定した軸受の負荷分布と軸受における転動体と内外輪との接触角とから、前記軸受の使用された使用条件であるラジアル荷重とアキシアル荷重とを推定する過程と、
を含む転がり軸受の使用条件推定方法。
【請求項２】
請求項１おいて、前記軸受における転動体と内外輪との接触角を、転がり軸受の固定輪の転走面から推定する転がり軸受の使用条件推定方法。
【請求項３】
請求項１おいて、前記軸受における転動体と内外輪との接触角を、転がり軸受の設計値から決定する転がり軸受の使用条件推定方法。
【請求項４】
請求項１ないし請求項３のいずれか１項において、前記使用条件を推定する転がり軸受がラジアル軸受である転がり軸受の使用条件推定方法。
【請求項５】
請求項１ないし請求項３のいずれか１項において、前記使用条件を推定する転がり軸受がスラスト軸受である転がり軸受の使用条件推定方法。
【請求項６】
使用後の転がり軸受の回転輪と固定輪に対してＸ線分析を行い、このＸ線分析の結果から回転輪の最大転動体荷重の推定値を得ると共に固定輪の任意の１点以上の位置での転動体荷重とを得る過程と、
この過程で得られた前記回転輪の最大転動体荷重の推定値と固定輪の任意の１点以上の位置での転動体荷重とから軸受の負荷分布を推定する過程と、
この推定した軸受の負荷分布と、負荷回数Ｎ、繰り返し応力Ｓ、Ｘ線分析で求まる半価幅ｗ( °) の関係を求めておいた結果とから、軸受の使用された回転回数を推定する過程と、
を含む転がり軸受の使用条件推定方法。
【請求項７】
請求項６において、負荷回数Ｎ、繰り返し応力Ｓ、Ｘ線分析で求まる半価幅ｗ( °) の関係を下式に当てはめる余寿命推定方法。
【数１】
2011069684.tif 000021
ただし、ｗ0 は未疲労での半価幅、σY は降伏応力、f 、g 、h 、k は正の定数である。
【請求項８】
請求項７において、前記繰り返し応力Ｓは相当応力σe である余寿命推定方法。
【請求項９】
請求項６ないし請求項８のいずれか１項において、前記軸受における転動体と内外輪との接触角を、転がり軸受の固定輪の転走面から推定する転がり軸受の使用条件推定方法。
【請求項１０】
請求項６ないし請求項８のいずれか１項において、前記軸受における転動体と内外輪との接触角を、転がり軸受の設計値から決定する転がり軸受の使用条件推定方法。
【請求項１１】
請求項６ないし請求項１０のいずれか１項において、前記使用された回転回数を推定する転がり軸受がラジアル軸受である転がり軸受の使用条件推定方法。
【請求項１２】
請求項６ないし請求項１０のいずれか１項において、前記使用された回転回数を推定する転がり軸受がスラスト軸受である転がり軸受の使用条件推定方法。
【請求項１３】
使用後の転がり軸受の回転輪と固定輪に対してＸ線分析を行い、このＸ線分析の結果から回転輪の最大転動体荷重の推定値を得ると共に、固定輪の任意の１点以上の位置での転動体荷重の推定値を得る過程と、
この過程で得られた回転輪の最大転動体荷重の推定値と固定輪の任意の１点以上の位置での転動体荷重の推定値とから軸受の負荷分布を推定する過程と、
この推定した軸受の負荷分布と軸受における転動体と内外輪との接触角とから、前記軸受の使用された使用条件であるラジアル荷重とアキシアル荷重とを推定する過程と、
前記の推定した軸受の負荷分布と、負荷回数Ｎ、繰り返し応力Ｓ、Ｘ線分析で求まる半価幅ｗ( °) の関係を求めておいた結果とから、軸受の使用された回転回数を推定する過程と、
を含む転がり軸受の使用条件推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、転がり軸受の使用された使用条件の推定方法に関し、より具体的にはＸ線分析により軸受の使用された使用条件であるラジアル荷重、アキシアル荷重、回転回数等を推定する方法に関する。
【背景技術】
【０００２】
転がり軸受の寿命は、使用荷重、潤滑条件、材料等に依存することが知られている。従来より、軸受の寿命予測は、使用荷重、潤滑条件、材料等を考慮して作成された寿命計算式を使って行われている（非特許文献１）。この計算式は、転がり軸受をある条件で使用する際にどのくらいの期間使用できるかを見積もるため、あるいは、要求される使用期間で軸受が破損しないためにどのような条件で転がり軸受を使用すればよいかを見積もるために使用されている。
【０００３】
一般に、軸受は寿命計算式に基づいて設定した使用条件で使用される。したがって、通常の条件で軸受が使用されるかぎりは、軸受の寿命が問題になることはないはずである。しかしながら、軸受の寿命が市場で問題となる状況がしばしば生じる。これは、実際の軸受の使用条件が設計した条件と異なっている場合があることが一因である。
【０００４】
設計した寿命よりも早期に破損した軸受では、その破損原因を推定するために、使用条件を推定する調査が行われる。軸受の使用条件の推定方法としては、(1) 使用温度推定（非特許文献２）、(2) Ｘ線分析による使用面圧推定（非特許文献２）、(3) 潤滑条件推定（特許文献１，２）、(4) 荷重推定（特許文献３，４）等がある。これらは、いずれも軸受の破損原因の推定や余寿命推定などに用いられるものであるが、この中で軸受にとって最も基本であり、最も重要な使用条件の推定は荷重推定である。これは転がり軸受の寿命を決める最も基本的な因子が動等価荷重と呼ばれる荷重に関係する因子であるためである。
【０００５】
式(1) に転がり軸受の寿命計算式を示す。
【数１】
2011069684.tif 000003
ラジアル軸受の場合の動等価荷重：Ｐr ＝ＸＦr ＋ＹＦa   ・・・(2) 
スラスト軸受の場合の動等価荷重：Ｐa ＝Ｆa ＋１．２Ｆr ・・・(3) 
【０００６】
ここで、Ｃは静定格荷重で既知の値(kgf) 、Ｐは動等価荷重(kgf) 、ｐは玉軸受で３、コロ軸受で10/3、Ｘはラジアル荷重係数で既知の値（非特許文献３）、Ｙはアキシアル荷重係数で既知の値（非特許文献３）、Ｆr はラジアル荷重(kgf) 、Ｆa はアキシアル荷重(kgf) である。
【０００７】
動等価荷重は、ラジアル荷重Ｆr (kgf) とアキシアル荷重Ｆa (kgf) から求められる。しかし、使用後の軸受から、ラジアル荷重Ｆr (kgf) とアキシアル荷重Ｆa (kgf) を合理的に推定する方法は無かった。
一方、軸受の使用条件ではないものの、軸受の破損原因の推定には、軸受がどのくらいの時間使用されていたか（＝軸受の回転回数）を知ることも重要である。しかし、軸受の回転回数を使用後の軸受から推定する方法は存在しない。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開２００５−３４５１３２号公報
【特許文献２】特開２００４−２０３７８号公報
【特許文献３】特開２００１−１２４６６５公報
【特許文献４】特開２００２−２５７７９７号公報
【非特許文献】
【０００９】
【非特許文献１】岡本純三著, ころがり軸受・ころ軸受の動的負荷容量−ルンドベルグとパルムグレン（Lundberg-Palmgren ）理論の詳解−，千葉大学工学部機械工学科機械要素講座，(1988)
【非特許文献２】対馬全之, 前田喜久男共著, ベアリングエンジニア, 48 (1984) 1-17.
【非特許文献３】ＮＴＮ社発行，ＮＴＮ転がり軸受総合カタログ, CAT. No202- ＶII/J, (2002).
【非特許文献４】Ｔ．Ａ．ハリス等（T. A. Harris et. al.）著, 転がり軸受の解析（Rolling Bearing Analysis） 5th ed., CPC Press, (2006) , 106p.
【非特許文献５】Ｋ．Ｔ．ジョンソン（K. L. Johnson ）著 ,接触理論（ Contact Mechanics）, (1989), 102p.
【非特許文献６】Ｘ線応力測定法標準(2002 年度版)-鉄鋼編, 日本材料学会, (2002)
【発明の概要】
【発明が解決しようとする課題】
【００１０】
転がり軸受のラジアル荷重Ｆr (kgf) 、アキシアル荷重Ｆa (kgf) 、軸受の回転回数の推定は、その破損原因の推定にとって重要なものであるが、その推定方法として合理的な方法が無かった。
【００１１】
この発明の目的は、使用後の軸受から、その使用条件であるラジアル荷重およびアキシアル荷重を合理的に推定することのできる転がり軸受の使用条件推定方法を提案することにある。
この発明の他の目的は、使用後の軸受から、その使用条件である軸受の回転回数を推定することのできる転がり軸受の使用条件推定方法を提案することにある。
【課題を解決するための手段】
【００１２】
この発明における第１の転がり軸受の使用条件推定方法は、使用後の転がり軸受の回転輪と固定輪に対してＸ線分析を行い、このＸ線分析の結果から回転輪の最大転動体荷重の推定値を得ると共に、固定輪の任意の１点以上の位置での転動体荷重の推定値を得る過程（Ｓ１）と、
この過程（Ｓ１）で得られた回転輪の最大転動体荷重の推定値と固定輪の任意の１点以上の位置での転動体荷重の推定値とから軸受の負荷分布を推定する過程（Ｓ２）と、
この推定した軸受の負荷分布と軸受における転動体と内外輪との接触角とから、前記軸受の使用された使用条件であるラジアル荷重とアキシアル荷重とを推定する過程（Ｓ３）とを含む。
【００１３】
この方法によると、使用後の転がり軸受のラジアル荷重およびアキシアル荷重を、Ｘ線分析を使って合理的に調べることができる。そのため、軸受の破損原因の正確な推定の一助となり、軸受の破損原因の調査をより詳細に実施できるようになる。
なお、上記Ｘ線分析では、Ｘ線照射により得られたデータから、表面からの深さに対応する残留応力と半価幅の分布を測定する。この残留応力あるいは半価幅の分布の結果から、回転輪の最大転動体荷重の推定値を得る。以下で言うＸ線分析においても、これと同様に残留応力と半価幅分布を測定する。
【００１４】
この発明において、前記軸受における転動体と内外輪との接触角は、転がり軸受の固定輪の転走面から推定しても良い。前記転動体と内外輪との接触角は、転がり軸受の設計値から決定しても良い。
【００１５】
この発明方法で使用条件を推定する転がり軸受は、ラジアル軸受であっても良く、またスラスト軸受であっても良い。
【００１６】
この発明における第２の転がり軸受の使用条件推定方法は、使用後の転がり軸受の回転輪と固定輪に対してＸ線分析を行い、このＸ線分析の結果から回転輪の最大転動体荷重の推定値を得ると共に固定輪の任意の１点以上の位置での転動体荷重とを得る過程（Ｓ１）と、
この過程（Ｓ１）で得られた前記回転輪の最大転動体荷重の推定値と固定輪の任意の１点以上の位置での転動体荷重とから軸受の負荷分布を推定する過程（Ｓ２）と、
この推定した軸受の負荷分布と、あらかじめ負荷回数Ｎ、繰り返し応力Ｓ、Ｘ線分析で求まる半価幅ｗ( °) の関係を求めておいた結果とから、軸受の使用された回転回数を推定する過程（Ｓ４）とを含む。
【００１７】
この方法によると、使用後の転がり軸受の回転回数を、Ｘ線分析を使って合理的に調べることができる。そのため、軸受の破損原因の正確な推定の一助となり、軸受の破損原因の調査をより詳細に実施できるようになる。
【００１８】
第２の転がり軸受の使用条件推定方法において、負荷回数Ｎ、繰り返し応力Ｓ、Ｘ線分析で求まる半価幅ｗ( °) の関係を下式に当てはめても良い。
【数２】
2011069684.tif 000004
ただし、ｗ0 は未疲労での半価幅、σY は降伏応力、f 、g 、h 、k は正の定数である。
この場合に、前記繰り返し応力Ｓは相当応力σe であっても良い。
【００１９】
第２の転がり軸受の使用条件推定方法における前記各方法において、前記軸受における転動体と内外輪との接触角を、転がり軸受の固定輪の転走面から推定しても良い。また、前記軸受における転動体と内外輪との接触角を、転がり軸受の設計値から決定しても良い。
【００２０】
第２の転がり軸受の使用条件推定方法において、前記使用された回転回数を推定する転がり軸受は、ラジアル軸受であっても良く、またスラスト軸受であっても良い。
【００２１】
この発明の第３の転がり軸受の使用条件推定方法は、使用後の転がり軸受の回転輪と固定輪に対してＸ線分析を行い、このＸ線分析の結果から回転輪の最大転動体荷重の推定値を得ると共に、固定輪の任意の１点以上の位置での転動体荷重の推定値を得る過程（Ｓ１）と、
この過程（Ｓ１）で得られた回転輪の最大転動体荷重の推定値と固定輪の任意の１点以上の位置での転動体荷重の推定値とから軸受の負荷分布を推定する過程（Ｓ２）と、
この推定した軸受の負荷分布と軸受における転動体と内外輪との接触角とから、前記軸受の使用された使用条件であるラジアル荷重とアキシアル荷重とを推定する過程（Ｓ３）と、
前記の推定した軸受の負荷分布と、あらかじめ負荷回数Ｎ、繰り返し応力Ｓ、Ｘ線分析で求まる半価幅ｗ( °) の関係を求めておいた結果とから、軸受の使用された回転回数（Ｓ４）を推定する過程と、を含む。
【００２２】
この方法によると、使用後の転がり軸受のラジアル荷重、アキシアル荷重、および回転回数を、Ｘ線分析を使って合理的に調べることができる。そのため、軸受の破損原因の正確な推定の一助となり、軸受の破損原因の調査をより詳細に実施できるようになる。
【発明の効果】
【００２３】
この発明の第１の転がり軸受の使用条件推定方法は、使用後の転がり軸受の回転輪と固定輪に対してＸ線分析を行い、このＸ線分析の結果から回転輪の最大転動体荷重の推定値を得ると共に、固定輪の任意の１点以上の位置での転動体荷重の推定値を得る過程と、この過程で得られた回転輪の最大転動体荷重の推定値と固定輪の任意の１点以上の位置での転動体荷重の推定値とから軸受の負荷分布を推定する過程と、この推定した軸受の負荷分布と軸受における転動体と内外輪との接触角とから、前記軸受の使用された使用条件であるラジアル荷重とアキシアル荷重とを推定する過程とを含む方法であるため、使用後の転がり軸受のラジアル荷重およびアキシアル荷重を、Ｘ線分析を使って合理的に調べることができる。そのため、軸受の破損原因の正確な推定の一助となり、軸受の破損原因の調査をより詳細に実施できるようになる。
【００２４】
この発明の第２の転がり軸受の使用条件推定方法は、使用後の転がり軸受の回転輪と固定輪に対してＸ線分析を行い、このＸ線分析の結果から回転輪の最大転動体荷重の推定値を得ると共に固定輪の任意の１点以上の位置での転動体荷重とを得る過程と、この過程で得られた前記回転輪の最大転動体荷重の推定値と固定輪の任意の１点以上の位置での転動体荷重とから軸受の負荷分布を推定する過程と、この推定した軸受の負荷分布と、負荷回数Ｎ、繰り返し応力Ｓ、Ｘ線分析で求まる半価幅ｗ( °) の関係を求めておいた結果とから、軸受の使用された回転回数を推定する過程とを含む方法であるため、使用後の転がり軸受の回転回数を、Ｘ線分析を使って合理的に調べることができる。そのため、軸受の破損原因の正確な推定の一助となり、軸受の破損原因の調査をより詳細に実施できるようになる。
【００２５】
この発明の第３の転がり軸受の使用条件推定方法は、使用後の転がり軸受の回転輪と固定輪に対してＸ線分析を行い、このＸ線分析の結果から回転輪の最大転動体荷重の推定値を得ると共に、固定輪の任意の１点以上の位置での転動体荷重の推定値を得る過程と、この過程で得られた回転輪の最大転動体荷重の推定値と固定輪の任意の１点以上の位置での転動体荷重の推定値とから軸受の負荷分布を推定する過程と、この推定した軸受の負荷分布と軸受における転動体と内外輪との接触角とから、前記軸受の使用された使用条件であるラジアル荷重とアキシアル荷重とを推定する過程と、前記の推定した軸受の負荷分布と、負荷回数Ｎ、繰り返し応力Ｓ、Ｘ線分析で求まる半価幅ｗ( °) の関係を求めておいた結果とから、軸受の使用された回転回数を推定する過程とを含む方法であるため、使用後の転がり軸受のラジアル荷重、アキシアル荷重、および回転回数を、Ｘ線分析を使って合理的に調べることができる。そのため、軸受の破損原因の正確な推定の一助となり、軸受の破損原因の調査をより詳細に実施できるようになる。
【図面の簡単な説明】
【００２６】
【図１】この発明の一実施形態に係る転がり軸受の使用条件推定方法を示す流れ図である。
【図２】ラジアル軸受にラジアル荷重Ｆr がのみが作用している場合の負荷帯の模式図である。
【図３】時刻ｔ＝０の時に転動体に最大の荷重が作用している状況の説明図である。
【図４】Ｘ分析の結果となる残留応力分布図の一例を示す説明図である。
【図５】接触面下表面下応力とＺΣρ値との関係を示すグラフである。
【図６】Ｘ分析の測定原理の説明図である。
【図７】この実施形態の使用条件推定方法の推定対象とする各種ラジアル軸受の例を示す断面図である。
【図８】この実施形態の使用条件推定方法の推定対象とするスラスト軸受の一例を示す断面図である。
【発明を実施するための形態】
【００２７】
この発明の一実施形態を図面と共に説明する。この転がり軸受の使用条件推定方法は、使用後の転がり軸受の回転輪と固定輪に対してＸ線分析を行い、このＸ線分析の結果から回転輪の最大転動体荷重の推定値を得ると共に、固定輪の任意の１点以上の位置での転動体荷重の推定値を得る過程（Ｓ１）と、この過程（Ｓ１）で得られた回転輪の最大転動体荷重の推定値と固定輪の任意の１点以上の位置での転動体荷重の推定値とから軸受の負荷分布を推定する過程（Ｓ２）と、この推定した軸受の負荷分布と軸受における転動体と内外輪との接触角とから、前記転がり軸受の使用された使用条件であるラジアル荷重とアキシアル荷重とを推定する過程（Ｓ３）と、前記の推定した軸受の負荷分布と、あらかじめ負荷回数Ｎ、繰り返し応力Ｓ、Ｘ線分析で求まる半価幅ｗ( °) の関係を求めておいた結果とから、軸受の使用された回転回数を推定する過程（Ｓ４）とを含む。
【００２８】
この方法によると、ラジアル荷重、アキシアル荷重、および軸受の回転回数を、Ｘ線分析を使って推定することができるので、軸受の破損原因の正確な推定の一助となり、軸受の破損原因の調査をより詳細に実施できるようになる。
【００２９】
以下では、Ｘ線を使った転がり軸受のラジアル荷重Ｆr  (kgf)、アキシアル荷重Ｆa  (kgf)、軸受の回転回数の推定方法について、具体的に説明していく。
まずはじめに、ラジアル荷重Ｆr  (kgf)とアキシアル荷重Ｆa  (kgf)を求める方法を説明する。転がり軸受のラジアル荷重Ｆr  (kgf)とＦa アキシアル荷重 (kgf)は以下の式(4),(5) で求められる（非特許文献１）。
【００３０】
【数３】
2011069684.tif 000005
ここで、Ｊr は負荷率εできまる定数でラジアル積分値、Ｊa は負荷率εで決まる定数でアキシアル積分値、Ｚは転動体個数、Ｑmax は最大転動体荷重(kgf) 、αは接触角(rad) である。
【００３１】
これより、ラジアル荷重Ｆr  (kgf)とアキシアル荷重Ｆa  (kgf)は、ε、Ｑmax  (kgf)、α(rad) の３つが求まれば推定できることになる。まず、従来のＸ線分析による負荷推定（非特許文献２）を軸受の回転輪で実施し、その軸受の最大接触面圧Ｐmax(kgf/mm2 を推定する。その後、最大接触面圧Ｐmax(kgf/mm2)から最大転動体荷重Ｑmax (kgf) を従来の方法で求める（非特許文献４）。また、α(rad) は軸受の転走跡の跡から、あるいは軸受の設計の諸元から推測できる。
【００３２】
上記Ｘ線分析では、Ｘ線照射により得られたデータから、表面からの深さに対応する主せん断応力の残留応力の分布を測定する。上記Ｘ線分析による負荷推定は、この残留応力分布の測定結果から、回転輪の最大転動体荷重の推定値を得る推定方法であり、次のように行う。
軸受に作用するヘルツの最大接触面圧Ｐmax と軸受の接触楕円短軸半径とには、ある一定の関係が存在する。なお、接触楕円は、荷重が加わったときに接点のまわりに生じる楕円状の接触面を言う。また、軸受使用時に作用する主せん断応力分布のピーク位置と軸受の接触楕円短軸半径とにも、ある一定の関係が存在する。したがって、軸受使用時に生成する主せん断応力分布のピーク位置が測定できれば、軸受に作用するヘルツの最大接触面圧Ｐmax を求めることができる。負荷荷重の推定は、軸受使用時に作用する主せん断応力分布のピーク位置を残留応力分布測定によって予測し、軸受に作用するヘルツの最大接触面圧Ｐmax を求める分析方法である。
【００３３】
負荷荷重推定の具体的方法を説明する。
高い接触面圧の下では転がり接触により最大せん断応力が作用する位置に対応する深さＺ45°に圧縮応力のピークを持つ残留応力分布が生成される（図４−(1) ピーク例）。このピーク位置を用いて接触面圧Ｐmax を推定することができる。実機で使用された軸受は、異物の噛み込みや温度上昇などにより、残留応力分布にピークが認められない場合がある（図１−(2) 生成深さ例）。このような場合は、残留応力の生成深さから接触面圧Ｐmax を推定する。上記Ｘ分析では、図４の残留応力分布図を求める。
〔圧縮残留応力のピーク位置から負荷荷重を推定する場合〕
図５に示す接触面下のせん断応力分布の横軸（Σρ×Ｚ）の深さＺにピーク位置（mm）を当てはめ、計算されたΣρ×Ｚ（Σρ：接触する物体間の曲率和）値（図５中、Ａ）を点線上にたどり、直線Ｐと交わったところを読み取る。
〔圧縮残留応力の生成深さ位置から負荷荷重を推定する場合〕
同じく、図５の横軸のＺに生成深さを当てはめ、計算されたΣρ×Ｚ値（図５中、Ｂ）を点線上にたどり、τC ＝600 ＭＰａ（臨界せん断応力：材料内部に残留応力が生成するために必要な最低限のせん断応力）の線と交わったところを読み取る。
【００３４】
また、上記の最大接触面圧Ｐmax(kgf/mm2)から最大転動体荷重Ｑmax (kgf) を求める方法では、次のように求める。
今、Ｘ線分析でＰmax が得られたとすれば、以下の式から転動体荷重Q 、接触楕円の長軸半径a,短軸半径ｂをそれぞれ求めることができる。
【数４】
2011069684.tif 000006
ここで、Ｒx1，Ｒx2，Ｒy1，Ｒy2は、ｘ，ｙ方向の２物体１，２（図示せず）の曲率半径である。最大接触面圧Ｐmax はＸ線分析で既に得られているので、結果、転動体荷重Ｑが求まることになる。
【００３５】
次に、負荷率εを求める手順を説明する。一般に、ラジアル軸受にラジアル荷重Ｆr(kgf)とアキシアル荷重Ｆa  (kgf)の合成荷重が加わると、個々の転動体に加わる荷重は均一ではなくなる。ルンドベルグとパルムグレン（LundbergとPalmgren）（非特許文献１）は、この問題について検討し、任意の角度 (rad)ごとに転動体荷重の関係である式(6) を求めた。
【００３６】
【数５】
2011069684.tif 000007
ここで、Ｑmax は最大転動体荷重、αは接触角、ｔはヘルツの接触理論から点接触では１．１、線接触では１．５である。
【００３７】
今、回転輪の負荷推定の結果からＱmax (kgf) は既知であるから、εを求めるためには、固定輪のある角度に作用していた転動体荷重Ｑ(kgf) を求める必要がある。一般には、固定輪の正確な負荷位置を求めることは難しいので、固定輪の２点で負荷推定を行うことになる。具体的には、固定輪の負荷域中の任意の２点ψ1  (rad)、ψ2 ＝ψ1 ＋φ (rad)に作用していた荷重を負荷推定で求め、以下の連立方程式を解けばεを求めることができる。ここで、ψ1  (rad)からずらした角度φ (rad)は記録しておく必要がある。
【数６】
2011069684.tif 000008
【００３８】
この式であれば、1 変数非線形方程式であるので容易に解くことができる。この式が解ければ、ψ1  (rad)は求められるから、その結果を式(7) か式(8) に代入すればεを求めることができる。その後、式(4) と式(5) から、ラジアル荷重Ｆr  (kgf)とスラスト荷重Ｆa  (kgf)を求める（必要であれば、式(2) と式(3) から動等価荷重(kgf) も計算できる）。
【００３９】
以下では、軸受の回転回数が測定されていない場合でも適用できるＸ線分析による回転回数の推定方法について説明していく。この方法は２段階である。まず、負荷推定で得た結果から最大転動体荷重Ｑmax  (kgf)を求め、その値と軸受の諸元から、回転輪が回転したときに固定輪に作用する内部応力を転動体が通過するごとに求める。次に、内部応力が作用したときの半価幅の低下量を実験で求めた半価幅、応力、負荷回数の関係から見積もり、負荷を受ける前の半価幅から引く。この計算を転動体が通過するごとに繰り返し、内外輪で起こる半価幅の低下をシミュレートする。最後に、この半価幅の低下が測定値と一致したときの軸受の回転回数を求める。
【００４０】
まずはじめに、回転輪（内輪）が回転したときに固定輪（外輪）に作用する内部応力を転動体が通過するごとに求める方法について説明する。図２に、ラジアル軸受にラジアル荷重Ｆr  (kgf)のみが作用しており、軸受の使用すきまが０である場合の負荷帯を示す（すきま＝０、ε＝０．５）。
【００４１】
このとき、図２の角度ψ (rad)ごとの荷重分布Ｑ(kgf) は式(6) で表される。負荷率εと最大転動体荷重Ｑmax  (kgf)は既に分析によって分かっている。したがって、式(6) より転動体荷重Ｑ(kgf) は分かるので、それぞれの角度ψでの最大接触面圧Ｐmax(kgf/mm2)と接触楕円の長・短軸半径ａ，ｂ (mm) は、従来の方法で推定できる（非特許文献４）。すなわち、接触楕円の長・短軸半径ａ，ｂは、段落〔００３４〕で前述した式により求まる。
【００４２】
一方、内外輪に作用する内部応力は、平面ひずみ状態と仮定することができるので、接触中心直下の深さｚ(mm)に対する各種応力成分σx ，σy ，σz ，τxy，τyz，τzx(kgf/mm2) は以下の式から計算できる（非特許文献５）。
【００４３】
【数７】
2011069684.tif 000009
ここで、σx ，σy ，σz  (kgf/mm2)は垂直応力成分、τxy，τyz，τzx (kgf/mm2)はせん断応力成分、νはポアソン比で０．３である。
【００４４】
したがって、接触中心直下の深さｚ(mm)における相当応力は式(16)から計算できる。
【数８】
2011069684.tif 000010
【００４５】
以上から、深さごとの相当応力は軸受内外輪の各位置ψ (rad)で計算することができる。しかし、この応力が軸受の回転時にどのような状況で繰り返されるかについては不明である。そこで、軸受内に作用する応力がどのような状況で繰り返されるかについて考える。
【００４６】
今、軸受の内輪がｎi (min-1) の速度で回転する場合を考える。また、時間ｔ＝０ (min)の時、転動体は、最大の転動体荷重Ｑmax  (kgf)が生じる位置にあるとする。その状況を図３に示す。
【００４７】
図中の内輪のＡ点は、時間ｔ＝０ (min)においてＱmax  (kgf)の負荷を受けている。この状態から、内輪を回転させていくと、内輪の回転速度ｎi (min-1) は式(17)で与えられる転動体の公転速度ｎe (min-1) よりも速いので、内輪のＡ点は次の転動体に追いつき、負荷を受けることになる。次の転動体にＡ点が負荷を受ける時間ｔ1 (min) は、内輪の回転位置と転動体の公転位置＋位相差が一致する点であるから、ｔ1 ｄmin)は式(18)で求めることができる。
【００４８】
【数９】
2011069684.tif 000011
ここで、Ｚは転動体個数、Ｄa は転動体直径(mm)、ｄp は転動体ピッチ円径(mm)、αは接触角(rad) 、ｎ0 は外輪の回転速度(min-1) である。
【００４９】
これより、内輪のＡ点が１個目の転動体から負荷を受ける時間ｔ1  (min)は式(20)になる。
【数１０】
2011069684.tif 000012
【００５０】
同様に、内輪のＡ点がａ個目の転動体から負荷を受けるｔa 時間（ min-1）とそのときのＡ点の角度ψa  (rad)は式(21)と式(22)から計算できる。
【数１１】
2011069684.tif 000013
【００５１】
また、内輪のＡ点が個目の転動体から負荷を受ける荷重Ｑ（ψa ） (kgf)は式(23)になる。
【数１２】
2011069684.tif 000014
【００５２】
以上より、内輪Ａ点がａ個目の転動体から負荷を受ける角度ψa  (rad)が計算でき、そのときの転動体荷重Ｑ（ψa ） (kgf)が計算できるので、最終的に、角度（ψa ） (rad)での内部応力を求めることができる。ここで、以上の計算は内輪の１点のみに着目しているが、軸受の回転回数が多くなれば、負荷がランダムになり、どの位置においても同じ負荷を受けると考えることができる。したがって、位相の異なる他の領域について負荷の状況を別途考える必要は無い（非特許文献１）。また、この計算は内輪回転での計算であるが、外輪回転でも同様な計算ができる。
【００５３】
次に、内部応力から半価幅ｗ( °) の変化がどのように起こっていくかを計算する手順を説明する。鋼の疲労の程度は、繰り返し応力σe  (kgf/mm2)の大きさと負荷回数Ｎによって決まる。また、半価幅ｗ( °) の変化は鋼の疲労の程度を表す。したがって、半価幅ｗ( °) 、繰り返し応力σe  (kgf/mm2)、負荷回数Ｎの関係が既知であれば、繰り返し応力と半価幅ｗ( °) から負荷回数N を見積もることができる。今、半価幅ｗ( °) 、繰り返し応力σe  (kgf/mm2)、負荷回数Ｎの関係は式(24)で表せると仮定する。
【００５４】
【数１３】
2011069684.tif 000015
ここで、ｗ0 は未疲労での半価幅ｗの値( °) 、ｆ，ｇ，ｈ，ｋは正の定数、
σe は繰り返し応力で相当応力の値(kgf/mm2) 、σY は降伏応力106kgf/mm2である。
【００５５】
この式は、全く物理的な意味を持たないが、半価幅ｗ( °) が負荷回数Ｎと繰り返し応力σe  (kgf/mm2)の増加に対して単純減少する形で、負荷回数Ｎと繰り返し応力σe  (kgf/mm2)の増加に対しては線形、非線形どちらの変化であっても実験結果を適合させることができる。また、繰り返し応力σe  (kgf/mm2)には塑性変形の程度を表す相当応力σe  (kgf/mm2)を適用することとし、降伏応力σY  (kgf/mm2)以下では、半価幅ｗ( °) の低下は起こらないので、式の形はそれを考慮した形になっている。ここで、降伏応力σY (kgf/mm2) としては、600 ×√3 ≒106kgf/mm2を採用した。ｆ，ｇ，ｈ，ｋの定数は、実験結果（実験で得た半価幅ｗ( °) 、繰り返し応力σe  (kgf/mm2)、負荷回数Ｎの関係）を式(24)で非線形重回帰分析することによって決定する。これより、負荷回数Ｎは任意の繰り返し応力σe  (kgf/mm2)と半価幅ｗ( °) に対して求めることができるようになる。
【００５６】
次に、半価幅ｗ( °) の低下を計算していく手順を具体例で説明していく。今、６２０６玉軸受にラジアル荷重700kgfが作用しており、時間ｔ＝０min において、軸受が図２の状態から、内輪が１min -1で回転し始めるとする。まず、ｔ＝０ minにおける内輪Ａ点直下の最大転動体荷重Ｑmax (kgf) を計算する。今、時間ｔ＝０min 、荷重；Ｆr ＝７００kgf 、転動体数Ｚ＝９、接触角α＝０ (rad)、転動体ピッチ直径ｄP ＝４５．５mm、転動体直径ＤP ＝９．５２５mm、玉軸受のラジアル積分値Ｊr ＝０．２２８８（ε＝０．５）、内輪回転数ｎi (min-1) であるから、Ａ点の転動体荷重は以下のように求められる。
【００５７】
【数１４】
2011069684.tif 000016
【００５８】
次に、得られた転動体荷重から、接触楕円の長軸半径ａ (mm) と短軸半径ｂ (mm) 及び最大接触面圧Ｐmax  (kgf/mm2)を計算する。今、各曲率半径は６２０６玉軸受の諸元から、Ｒx1mm、Ｒx2 mm 、Ｒy1 mm 、Ｒy2 mm であるから、ａ (mm) 、ｂ (mm) 、Ｐmax (kgf/mm2) は以下になる。
ａ＝２．７０５３５２mm
ｂ＝０．１９３２８３mm
Ｐmax ＝３０４．０４８kgf/mm2 
【００５９】
次に、得られたａ (mm) 、ｂ (mm) 、Ｐmax  (kgf/mm2 ）から、内部応力を計算する。今、例として、ｚ＝0.2mm の深さにおける内部応力を計算する。
【００６０】
【数１５】
2011069684.tif 000017
【００６１】
次に、上述の計算で得た繰り返し応力σe  (kgf/mm2)が１回加わった後の半価幅ｗ( °) を計算する。仮に、実験から半価幅ｗ( °) 、繰り返し応力σe  (kgf/mm2)、負荷回数Ｎの関係として式(26)が得られたとする。
【数１６】
2011069684.tif 000018
【００６２】
未疲労の半価幅ｗ0(°) が７°( 実際には実験で求められる) であるとすると、前述の計算で得た繰り返し応力σe が１回加わった後の半価幅ｗ( °) は以下のように計算できる。
【数１７】
2011069684.tif 000019
【００６３】
これは、７００kgf で荷重を受けた内輪のＡ点が時間t=0minにおいて、1 回の負荷を受けたとき、深さｚ＝０．２mmにおける半価幅ｗ( °) が1.44E-07( °) 低下することを示している。式(26)の形から分かるように、2 回目以降の負荷に対しては、半価幅ｗ( °) の変化は非線形性が出てくるので、負荷が2 回目以降の半価幅ｗ( °) の低下量を見積もる式は、この式を負荷回数Ｎに対して偏微分した式(27)になる。
【００６４】
【数１８】
2011069684.tif 000020
【００６５】
以上のように、半価幅ｗ( °) の低下の挙動は、軸受が回転するときの変化を１回の負荷ごとに計算して求めることができる。この計算を繰り返すと、半価幅ｗ( °) は次第に低下し、最終的には、実験で得た半価幅（この例では深さ位置ｚ＝０．２mmで測定した値）と一致するときがくる。このとき、軸受の回転角度ψa は同時に計算しているから、この角度を２πで割れば、軸受の回転回数が計算できる。以上のように軸受の回転回数は、Ｘ線分析により求められる。
【００６６】
次に、上記Ｘ線分析過程Ｓ１ａで行うＸ分析の測定原理を図６と共に説明する（非特許文献６を引用）。物質にある波長を持ったＸ線をある角度で入射すると、入射Ｘ線は、図６（Ａ）に示すように、物質に応じて特定の角度に回折される。この回折Ｘ線の角度と強度は、物質中の原子の種類とその構成に依存する。したがって、物質にＸ線を入射し、反射してきたＸ線の角度と強度を測定すれば、その物質の状態（相）を同定できる。この実施形態におけるＸ線分析はこのＸ線回折を利用したもので、本来は結晶の状態（相）を同定する分析になる。その回折角度と物質との相互作用は次のブラッグの式で表される（λは入射Ｘ線の波長で、n は整数）。
ブラッグ反射の条件式：２ｄ・sin θ＝ｎλ
この式中に格子面の間隔ｄが入っている。もし、圧縮や引張の応力働いていたら、物質中の面間隔は若干変化する。この面間隔の変化を測定しているのが、Ｘ線分析による残留応力の測定になる。一方、半価幅は面間隔のばらつきを示したものになる。焼入された物質は急速に起こる変態で格子面が揃っていない。したがって面間隔のばらつきは、ばらつきが大きくなっている。転動等の疲労を受けると、格子が揺さぶられるので、焼入れ前の鉄本来の面間隔に緩和していき、結果、Ｘ線分析値の半価幅が小さくなっていく。
【００６７】
なお、図７は、この実施形態の使用条件推定方法の推定対象とする各種ラジアル軸受の例を示し、図８は、実施形態の使用条件推定方法の推定対象とするスラスト軸受の一例を示す。
【符号の説明】
【００６８】
Ｓ１：Ｘ線分析・転動体荷重推定過程
Ｓ２：負荷分布推定過程
Ｓ３：ラジアル荷重．アキシアル荷重推定過程
Ｓ４：回転回数推定過程
【図１】
2011069684.tif 000022
【図２】
2011069684.tif 000023
【図３】
2011069684.tif 000024
【図４】
2011069684.tif 000025
【図５】
2011069684.tif 000026
【図６】
2011069684.tif 000027
【図７】
2011069684.tif 000028
【図８】
2011069684.tif 000029