加速試験における有為差有無判定・有為寿命差見積もり方法および装置
【課題】加速試験で得られた2つのロット間寿命から有為差有無の判定、および倍率による寿命差を算出する方法を提供する。
【解決手段】あるワイブル分布から水準1のロットの試験個数分の乱数を発生させ、同じワイブル分布から水準2のロットの試験個数分の乱数を発生させる。その乱数から求まる寿命を演算し1組の寿命比を算出する。この処理を設定回数繰り返し設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成し、設定信頼幅内の最大および最小の寿命比を読み取り、有為差有無の判定を行う。有為差有りと判定した場合、異なるワイブル分布に変えて、上記累積確率分布分析手順を繰り返す。繰り返しにより得られた設定倍率と信頼幅内の最大の寿命比の関係を示すグラフを作成し、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする。
【解決手段】あるワイブル分布から水準1のロットの試験個数分の乱数を発生させ、同じワイブル分布から水準2のロットの試験個数分の乱数を発生させる。その乱数から求まる寿命を演算し1組の寿命比を算出する。この処理を設定回数繰り返し設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成し、設定信頼幅内の最大および最小の寿命比を読み取り、有為差有無の判定を行う。有為差有りと判定した場合、異なるワイブル分布に変えて、上記累積確率分布分析手順を繰り返す。繰り返しにより得られた設定倍率と信頼幅内の最大の寿命比の関係を示すグラフを作成し、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする。
【発明の詳細な説明】
【技術分野】
【0001】
この発明は、軸受等の機械部品または試験片からなる試験対象品を使用環境よりも厳しい所定の環境条件におき、破損が発生した時間から寿命を算出する加速試験において、試験対象品の2つのロット間で有為差が有るか否かの判定、およびある場合にどの程度の有為差があるかを見積もる有為差有無判定・有為寿命差見積もり方法,装置、プログラムに関する。
【背景技術】
【0002】
寿命試験は、軸受等の機械部品の性能を評価するために欠かせない試験の1つである。寿命試験には、大きく分けて(1) 実機の使用環境に近い条件で試験を行う実機試験と、(2) 比較的過酷な条件で寿命試験を行う加速試験がある。前者は、製品が有限時間内に破損するケースが極めて少ないため、ある目標時間まで破損することなく試験が継続すれば、寿命は問題ないと判断する試験である(以下、このような試験を「打切り試験」と呼ぶ)。一方、後者は、比較的短時間で破損が発生するので、ワイブルプロットで寿命が算出でき(例えば非特許文献1)、その算出寿命から性能の優劣を判定する試験である(以下、このような試験を「加速試験」と呼ぶ)。
【0003】
従来より、寿命試験は経験を積んだ熟練者が行っており、試験条件や試験個数を決める寿命試験の設計と寿命試験結果の解釈に対して経験的に確からしい判断ができたと考えられる。
図14に、従来から行われてきた寿命試験の設計と寿命試験結果の解釈の手順を、打切り試験と加速試験ごとに示す。
また、現在、寿命試験において経験的に判断されているものの詳細を、表1に示す。
【0004】
【表1】
【0005】
なお、ワイブル分布を機械部品の寿命判断に用いるものは、種々の特許文献,非特許文献に提案されている。
【特許文献1】特開2006−040203号公報
【特許文献2】特開2002−277382号公報
【特許文献3】特開2005−226829号公報
【非特許文献1】真壁肇著、信頼性工学入門79、1991年発行
【発明の開示】
【発明が解決しようとする課題】
【0006】
加速試験は、上記のように破損時間から寿命を算出し、その算出寿命から性能の優劣を判定する試験であるが、2ロットの寿命試験結果から求めた算出寿命の間に、本当に有為差があるかどうかを確かめたい状況が多くある。
従来より、算出寿命の有為差判定は信頼幅という概念を使って行われてきた。しかし、この判定方法には、次のようにいくつかの問題がある。
【0007】
まず、信頼幅は、2水準間での寿命の有意差を判断できないという問題がある。その理由は、どの程度信頼幅が重なっていたら寿命に優位差があるのかということを定量的に判断できないためである。信頼幅を使って2水準間での寿命の有意差を判定する場合、次のような手順が必要になる。以下、その手順を図15を使って説明する。
1) L10寿命におけるメディアンランクの信頼水準の分布(図中Aの分布)を求める。ここで、この分布はF値表とファイサー(Fisher)の近似式を用いて求める。
2) L10寿命におけるメディアンランクの信頼水準の分布とワイブルスロープからL10のばらつき頻度の分布(図中Bの分布)を求める。
3) このばらつき分布f(x)を使って、図中の式の積分を取り、寿命倍率nが同じ寿命分布から発生する頻度を調べる。これは寿命比のばらつき分布である。
4) 3)の分布の累積確率分布を作成し、そのグラフから90%以上の確率で有意差有と判定できる寿命倍率を計算する。
【0008】
上のような手順で、従来の信頼幅を使って寿命の有意差検定や設計が可能になる。しかし、このような手順は複雑であり、実用化が難しい。このため、信頼幅を使った寿命の有意差検定が今まで行われてこなかったのは当然であるといえる。
【0009】
レオナード・ジー・ジョンソン(LEONARD G. JOHNSON)は、ワイブルスロープごと、試験個数ごとに、2水準間の平均寿命とL10 寿命の有意差を判定する方法を提案している)。しかし、その方法も、実際には使われていないのが現状である。これは、手軽な形で優位差検定を行うツールを残していないことに原因がある。
【0010】
また、従来は、有為差有りと判定を行っても、少なくともどの程度の寿命差があると言えるのかを、定量的に示す適切な方法がなかった。
【0011】
この発明の目的は、2つのロット間での有為差有無の判定を定量的に適切に行え、かつ有為差があると判定される場合に、少なくともどの程度の寿命差があると言えるのかを、簡単かつ迅速に試算することができ、かつ定量的に求められて信頼性の高いものとでき、熟練者でなくても必要な寿命差を適切に試算することのできる方法、装置、およびコンピュータプログラムを提供することである。
【課題を解決するための手段】
【0012】
この発明の加速試験における有為差有無判定・有為寿命差見積もり方法は、軸受等の機械部品または試験片からなる試験対象品を使用環境よりも厳しい所定の環境条件におき、破損が発生した時間から寿命を算出する加速試験において、試験対象品の2つのロットの寿命から、有為差有無の判定、および少なくとも断定できる、倍率による寿命差を算出する方法であって、
コンピュータに対し、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命を入力する過程(M1)と、
上記コンピュータに、上記有為差有無の判定、および少なくとも断定できる寿命差を演算させ演算結果を表示装置の画面に表示させるコンピュータ演算処理過程(M2)とを含む。
上記寿命は、例えばL10寿命(90%の信頼度の寿命)や、L50寿命(50%の信頼度の寿命)等の所定信頼度の寿命ある。
【0013】
上記コンピュータ演算処理過程(M2)として、
試験対象品の上記入力情報のワイブルスロープの値および仮設定寿命から定まるワイブル分布に従った乱数を水準1のロットの試験個数分発生させる第1乱数発生手順(N21)と、
第1乱数発生手順(N21)と同じワイブル分布に従った乱数を水準2のロットの試験個数分発生させる第2乱数発生手順(N22)と、
第1乱数発生手順(N21)で発生させた設定個数の乱数から所定の寿命計算方法で求まる寿命と第2乱数発生手順(N22)で発生させた設定個数の乱数から上記所定の寿命計算方法で求まる寿命とから、1組の水準1,水準2の寿命の寿命比を算出する寿命比算出手順(N23)と、
これら第1乱数発生手順(N21),第2乱数発生手順(N23),および寿命比算出手順(N24)を設定回数繰り返し、上記設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成して累積確率分布から、設定信頼幅内の最大の寿命比を読み取る累積確率分布分析手順(N24)と、
上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手順(N24)で得た最大寿命比を超える場合に有為差有りと判定し、最大寿命比以内の場合に有為差無しと判定する有為差有無判定手順(N25)とを含む。
【0014】
また、上記コンピュータ演算処理過程(N2)として、さらに、
上記有為差有無判定手順(N25)で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手順を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す倍率変更繰り返し手順(N26)と、
この倍率変更繰り返し手順(N26)で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する寿命倍率関係グラフ作成手順(N27)と、
上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする有為寿命差読み取り手順(N28)と、
上記有為差有無判定手順(N25)で判定した結果および上記有為寿命差読み取り手順(N28)で読み取った、少なくとも断定できる寿命差を表示装置に表示させる結果出力手順(N29)と、を含む。
【0015】
上記第1,第2乱数発生手順(N21,N22)において、乱数発生に用いるワイブル分布を特定するための寿命は、適宜想定した寿命または実績値を設定すれば良い。
また、寿命比算出手順(N23)において、乱数発生手順で発生させた試験個数の乱数から寿命を求める所定の寿命計算方法は、加速試験において従来から用いられている中の適宜の寿命計算方法を用いれば良い。加速試験においては、破損の発生した時間からワイブルプロットで寿命を求めることが行われており、このような破損時間から寿命を求める計算方法を、上記寿命比算出手順における寿命計算方法として用いる。
ワイブル分布は、次式、
【0016】
【数1】
【0017】
ただし、m:ワイブルスロープ、α:尺度因子、γ:最小寿命、
によって特定される。
【0018】
軸受等の機械部品の寿命は、ワイブル分布に従うとされている。ワイブル分布は、ワイブルスロープm、尺度因子α、最小寿命γの3つのパラメータを持っており、ワイブルスロープmによって指数分布、対数正規分布、正規分布を表現できる万能分布として知られている。量産される軸受等では、ワイブルスロープは実績値が既知である場合が多く、この発明方法において、ワイブルスロープには、試験対象品の実績値を用いることが好ましい。実績値がない場合は、適宜の方法で見積もったワイブルスロープを用いてもよい。最小寿命γは、種々の規格、例えばISO等によって計算方法が定められており、そのように定められたいずれかの計算方法を用いることが好ましい。尺度因子αは、ワイブルスロープの値、要求寿命の信頼度、要求寿命の値、および上記最小寿命γから一義的に決定される演算式があり、その演算式を用いて特定しても良い。
【0019】
この方法において、累積分布作成手順(N24)で行う処理は、同じ寿命分布を持つロットから験個数分の試験対象品を2組抜きとり、寿命試験を行った結果、どの程度の寿命差が出るのかを設定回数確認することに対応している。同じ寿命分布を持つロットから試験対象品を抜きとった試験においても、寿命はばらついてる。その寿命比のばらつきを累積確率分布の適宜の信頼区間、例えば5%と95%の区間(90%信頼区間)として定義すると、その寿命ばらつきは、上記累積確率分布から算出できる。この結果は、同じ寿命分布を持つロットから試験個数の試験対象品を2組抜きとって寿命比較を行っても、90%のものは、累積確率分布に示された信頼区間の範囲で寿命比がばらついてしまうことを示している。 逆にいうと上記信頼区間の最大の寿命差があれば、そのロット間の寿命分布が異なっている可能性が高いということになる。
したがって、上記有為差有無判定手順(N25)として、上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手順で得た最大寿命比を超える場合に有為差有りと判定し、最大寿命比以内の場合に有為差無しと判定することで、定量的に適切に有為差有無を判定することができる。
【0020】
また、寿命が異なるワイブル分布との組み合わせで、上記同様の累積確率分布の計算を行い、例えば、ワイブル分布と寿命が1.1 、1.2 …50倍異なるワイブル分布の組み合わせで計算を行うと、その組み合わせが1.1 、1.2 …50倍と変化するごとに、寿命差が1.1 、1.2 …50倍を中心にしてばらつくはずである。上記倍率変更繰り返し手順(N26)は、このような累積確率分布の計算、およびその確率分布についての上記確率分布分析手順(N24)を行う。
この倍率変更繰り返し手順(N26)で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成し、このグラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取ると、その読み取った値が少なくとも断定できる寿命差となる。
【0021】
この発明方法は、このようにワイブル分布に基づく乱数シミュレーションにより有為差有無を判定し、また少なくとも断定できる寿命差を求めるため、これら有為差有無および少なくとも断定できる寿命差を、簡単かつ迅速に試算することができる。また、定量的に求められて信頼性の高いものとでき、熟練者でなくても試算することができる。この乱数シミュレーションに用いるプログラムは、条件入力を行えば勝手に有為差有無を判定および少なくとも断定できる寿命差の検定を行ってくれる手軽なツールであることが優れた点である。
【0022】
この発明の加速試験における有為差有無判定・有為寿命差見積もり装置は、軸受等の機械部品または試験片からなる試験対象品を使用環境よりも厳しい所定の環境条件におき、破損が発生した時間から寿命を算出する加速試験において、試験対象品の2つのロットの寿命から、有為差有無の判定、および少なくとも断定できる、倍率による寿命差を算出する装置であって、
演算処理装置(1)と、この演算処理装置(1)の出力を画面に表示する表示装置(2)と、上記演算処理装置(1)に入力を行う入力手段(3)とを備える。
【0023】
上記演算処理装置(1)は、上記表示装置(2)の画面に、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命の入力を促す表示を行わせる促し画面出力手段(7M)と、実行命令に応答して上記有為差有無の判定および少なくとも断定できる寿命差の演算を行いその演算結果を上記表示装置(2)の画面に出力する有為差判定・有為寿命差演算手段(52)とを備える。
【0024】
上記有為差判定・有為寿命差演算手段(52)は、第1乱数発生手段(53)、第2乱数発生手段(54)、寿命比算出手段(55)、累積確率分布分析手段(56)、有為差有無判定手段(57)、倍率変更繰り返し手段(58)、寿命倍率関係グラフ作成手段(59)、有為寿命差読み取り手段(60)、および結果出力手段(61)を有する。
【0025】
第1乱数発生手段(53)は、上記入力情報のワイブルスロープの値および仮設定寿命から定まるワイブル分布に従った乱数を水準1のロットの試験個数分発生させる手段である。
第2乱数発生手段(54)は、第1乱数発生手段(53)と同じワイブル分布に従った乱数を水準2のロットの試験個数分発生させる手段である。
寿命比算出手段(55)は、第1乱数発生手段(53)で発生させた設定個数の乱数から所定の寿命計算方法で求まる寿命と、第2乱数発生手段(54)で発生させた設定個数の乱数から上記所定の寿命計算方法で求まる寿命とから、1組の水準1,水準2の寿命の寿命比を算出する手段である。
【0026】
累積確率分布分析手段(56)は、これら第1乱数発生手段(53),第2乱数発生手段(54),および寿命比算出手段(55)の処理を設定回数繰り返し、上記設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成して累積確率分布から、設定信頼幅内の最大の寿命比を読み取る手段である。
【0027】
有為差有無判定手段(57)は、上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手段(56)で得た最大寿命比を超える場合に有為差有りと判定し、最大寿命比以内の場合に有為差無しと判定する手段である。
【0028】
倍率変更繰り返し手段(58)は、有為差有無判定手段(57)で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手段(56)の処理を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す手段である。
寿命倍率関係グラフ作成手段(59)は、倍率変更繰り返し手段(58)で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する手段である。
有為寿命差読み取り手段(60)は、上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする手段である。
【0029】
結果出力手段(61)は、上記有為差有無判定手段(57)で判定した結果、および上記有為寿命差読み取り手段(60)で読み取った、少なくとも断定できる寿命差を表示装置(2)に表示させる手段である。
【0030】
この構成の加速試験における有為差有無判定・有為寿命差見積もり装置によると、この発明の有為差有無判定・有為寿命差見積もり方法に使用されて、有為差の有無の判定、および有為差があると判定される場合に、少なくともどの程度の寿命差があると言えるのかを、簡単かつ迅速に試算することができ、かつ定量的に求められて信頼性の高いものとでき、熟練者でなくても少なくとも断定できる寿命差を適切に試算することができる。
【0031】
この発明の加速試験における有為差有無判定・有為寿命差見積もりプログラム(51)は、コンピュータで実行可能なプログラムであって、
上記表示装置(2)の画面に、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命の入力を促す表示を行わせる促し画面出力手順(N1)と、
実行命令に応答して上記有為差有無の判定および少なくとも断定できる寿命差の演算を行いその演算結果を上記表示装置の画面に出力する有為差判定・有為寿命差演算手順(N2)とを含む。
【0032】
この有為差判定・有為寿命差演算手順(N2)は、
上記入力情報のワイブルスロープの値および仮設定寿命から定まるワイブル分布に従った乱数を水準1のロットの試験個数分発生させる第1乱数発生手順(N21)と、
第1乱数発生手順(N21)と同じワイブル分布に従った乱数を水準2のロットの試験個数分発生させる第2乱数発生手順(N22)と、
第1乱数発生手順(N21)で発生させた設定個数の乱数から所定の寿命計算方法で求まる寿命と第2乱数発生手順(N22)で発生させた設定個数の乱数から上記所定の寿命計算方法で求まる寿命とから、1組の水準1,水準2の寿命の寿命比を算出する寿命比算出手順(N23)と、
これら第1乱数発生手順(N21),第2乱数発生手順(N22),および寿命比算出手順(N23)を設定回数繰り返し、上記設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成して累積確率分布から、設定信頼幅内の最大の寿命比を読み取る累積確率分布分析手順(N24)と、
上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手順で得た最大寿命比を超える場合に有為差有りと判定し、最大寿命比以内の場合に有為差無しと判定する有為差有無判定手順(N25)とを備える。
【0033】
有為差判定・有為寿命差演算手順(N2)は、さらに、有為差有無判定手順(N25)で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手順を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す倍率変更繰り返し手順(N26)と、
この倍率変更繰り返し手順(N26)で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する寿命倍率関係グラフ作成手順(N27)と、
上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする有為寿命差読み取り手順(28)と、
上記有為差有無判定手順で判定した結果および上記有為寿命差読み取り手順で読み取った、少なくとも断定できる寿命差を表示装置に表示させる結果出力手順(29)と、
を含む。
【0034】
この構成の加速試験における有為差有無判定・有為寿命差見積もりプログラムは、この発明の有為差有無判定・有為寿命差見積もりの実施に使用され、有為差の有無の判定、および有為差があると判定される場合に、少なくともどの程度の寿命差があると言えるのかを、簡単かつ迅速に試算することができ、かつ定量的に求められて信頼性の高いものとでき、熟練者でなくても必要な寿命差を適切に試算することができる。
【発明の効果】
【0035】
この発明の加速試験における有為差有無判定・有為寿命差見積もり方法、装置、およびプログラムによると、コンピュータシミュレーションにより、試験対象品のワイブル分布から乱数を発生させ、同じワイブル分布から乱数を発生させ、その発生させた乱数から求まる寿命を演算して1組の寿命比を算出し、この処理を設定回数繰り返して設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成し、上記累積確率分布から設定信頼幅内の最大および最小の寿命比を読み取り、有為差有無の判定を行うため、有為差有無の判定を、簡単かつ迅速に行うことができ、かつ定量的に求められて信頼性の高いものとでき、熟練者でなくても判定することができる。
また、有為差有りと判定した場合に、寿命が異なるワイブル分布に繰り返し毎に変えて、上記累積確率分布分析手順を繰り返し、繰り返しにより得られた設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成し、このグラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とするため、有為差がある場合に、少なくとも断定できる寿命差を、簡単かつ迅速に求めることができ、かつ定量的に求められて信頼性の高いものとでき、熟練者でなくても適切に有為寿命差を求めることができる。
【発明を実施するための最良の形態】
【0036】
この発明の実施形態を説明する。この加速試験における有為差有無判定・有為寿命差見積もり方法は、軸受等の機械部品または試験片からなる試験対象品を使用環境よりも厳しい所定の環境条件におき、破損が発生した時間から寿命を算出する加速試験において、試験対象品の2つのロットの寿命から、有為差有無の判定、および少なくとも断定できる、倍率による寿命差を算出する方法である。
【0037】
以下、この実施形態を図面と共に説明する。この加速試験における有為差有無判定・有為寿命差見積もり方法は、図1に示すコンピュータ1に、乱数シミュレーションプログラムである加速試験における有為差有無判定・有為寿命差見積もりプログラム41を実行させることで行う。コンピュータ1はパーソナルコンピュータ等からなり、中央処理装置4およびメモリ5を有し、所定のオペレーションシステムによって動作するものである。コンピュータ1には、液晶表示装置等の画面によって表示可能な表示装置2と、キーボードやマウス等の入力装置3が接続され、あるいは付属して設けられている。コンピュータ1、表示装置2、入力装置3、および必要寿命差見積もりプログラム41により、図2に各機能達成手段をブロックで示した有為差有無判定・有為寿命差見積もり装置が構成される。同図の有為差有無判定・有為寿命差見積もり装置の構成については、後に説明する。
有為差有無判定・有為寿命差見積もりプログラム51は、図4および図5に流れ図で示す手順を備えるものである。同図の内容は、後に説明する。
【0038】
この必要寿命差見積もり方法は、図3に示すように、コンピュータ1に対して所定の情報を入力する入力過程M1と、コンピュータ1で演算処理を行って演算結果を出力するコンピュータ演算処理過程M2とからなる。
【0039】
入力過程M1では、図6に示すように所定の入力情報の入力を促す入力画面2aが、コンピュータ1の出力によって表示装置2に表示され、入力画面2a中に、所定の入力を促す表示が行われる。
入力画面2aは、同図の例では複数の画面に切り替えて表示される。入力画面のうち、第1の画面(図6(A)では、L10寿命とL50寿命とのいずれで有為差を検定するかを選択させる表示が行われる。この画面に対して、ラジオボタン等からなる表示部分を選択することで、L10寿命とL50寿命とのいずれかが選択できる。
第2の画面(図6(B))では、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値の入力を促す表示、および入力を記入する記入用表示が行われる。オペレータは、この記入用表示に、該当する値を入力する。
第2の画面(図6(C))では、比較する水準1および水準2の2つのロットの各試験個数n1,n2の入力、および試験結果の寿命を促す表示が行われ、その個数および寿命を入力する記入用表示が行われる。オペレータは、この記入用表示に対して該当する情報を入力する。
【0040】
ワイブルスロープの値には、試験の実績値を入力することが好ましい。実績値は10個以上の試験で得た結果を用いることが望ましく、より好ましくは20個以上の試験結果である。また、試験個数および試験結果は、実際に行った試験における個数および結果を入力する。なお、試験個数と寿命は、水準2の欄に長寿命側の入力をしなければないらいとする入力規則が、有為差有無判定・有為寿命差見積もりプログラム51に手順として設定されている。この手順は、例えば促し画面出力手順N1の一部として設けられる。以上の点に注意して、所定の条件を入力し、入力画面2aのOKボタンを押すと計算が開始される。
【0041】
図3のコンピュータ演算処理過程M2では、入力されたワイブルスロープの値、および2つのロットの各試験個数n1,n2、試験結果の寿命から、有為差有無の判定に必要な寿命差を計算して、有為差有無を判定し、有為差有りの場合には、少なくとも断定できる寿命差を計算する。
【0042】
計算が終了すると、図7に示す出力画面に結果が表示される。有為差がない場合は、同図(A)のように、有為差無しという表示と、有為差有りと判断するための必要寿命差が表示される。
有為差がある場合は、同図(B)のように、有為差ありという表示と、少なくとも断定できる寿命差(倍率)が表示される。
【0043】
図1の有為差有無判定・有為寿命差見積もりプログラム51は、コンピュータ1で実行可能なプログラムであって、図4,図5に流れ図で示す手順を備える。図4に示すように、有為差有無判定・有為寿命差見積もりプログラム51は、促し画面出力手順N1と必要寿命差演算手順N2とでなり、促し画面出力手順N1では、図6と共に前述した入力画面2aを出力する。この入力画面2aに対して、上記各入力情報が入力手段3から入力され、かつ入力画面2aのOKキーのクリック等によって実行命令が入力手段3から入力されると、有為差有無判定・有為寿命差見積もり手順N2が実行される。同図の入力画面に対して入力する過程が、図3の入力過程M1であり、同図のコンピュータ演算処理過程M2は図4の有為差有無判定・有為寿命差見積もり手順N2を実行する過程である。
【0044】
有為差有無判定・有為寿命差見積もり手順N2は、図5に流れ図で示す各手順で構成される。この流れ図には各手順M21〜M29毎の具体的な処理例を注釈として併記してある。
理解の容易のため、具体的数値を上げ、同図の具体的処理例を参照して、各手順N21〜N29を説明する。
【0045】
今、あるL10寿命を持つワイブル分布(図8(A))で、ワイブルスロープは4としている)である加速試験で、2ロットの試験を行ったとする。1ロット目(水準1)の試験個数n1が3個で、2ロット目(水準2)の試験個数n2も3個であるとする。
まず、そのワイブル分布に従って1ロット目の試験個数n1である3個の乱数を発生させ、その3個のデータからL10寿命とL50寿命を算出する(N21)。算出方法は、従来の加速試験の結果からの寿命算出に用いられている適宜の方法を使用する。なお、L10寿命だけ、またはL50寿命だけを演算しても良いが、この実施形態では、L10寿命とL50寿命のそれぞれについて必要寿命差を計算するために、L10寿命とL50寿命の両方を算出している。
次に、同じL10寿命を持つワイブル分布から、第2ロット目の試験個数n2である3個の乱数を3個発生させ、その3個のデータからL10寿命とL50寿命を上記と同じ方法で算出する(N22)。
【0046】
次に、得られた1組のL10寿命とL50寿命の寿命比をそれぞれ算出する。すなわちL10寿命同士の寿命比、およびL50寿命同士の寿命比を算出する(N23)。
【0047】
次に、これらの作業を設定回数(例えば1000回)繰り返し、設定回数組(1000組)のL10およびL50寿命の寿命比の確率分布(図8(B)と累積確率分布(図8(C)を作成する(N24)。
これらは、同じ寿命分布を持つロットから3個の試験片を2組抜きとり、寿命試験を行った結果、どの程度の寿命差が出るのかを設定回数(1000回)確認することに対応している。
【0048】
度数分布の図(図8(B))から、同じ寿命分布を持つロットから試験片を抜きとって試験を行ったにもかかわらず、寿命はばらついていることが分かる。その寿命比のばらつきを、所定信頼度区間、例えば累積確率分布の5%と95%の区間(90%信頼区間)として定義すると、その寿命ばらつきは、図8(C)に示したように0.39〜2.57倍であると算出できる(N243)。
この結果は、同じ寿命分布を持つロットから3個の試験片を2組抜きとって寿命比較を行っても、90%のものは0.39〜2.57倍の間で寿命比がばらついてしまうことを示している。逆にいうと、ばらつき内の最大値である2.57倍以上の寿命差があれば、そのロット間の寿命分布が異なっている可能性が高いということになるので、2ロットの比較試験で試験個数がそれぞれ3個の場合、有為差を判断できる寿命差は2.57倍以上であるといえる。
【0049】
そこで、上記ばらつき内の最大値である2.57倍と、入力された2つのロット(水準1,水準2)の寿命差(倍率)、すなわち〔(水準2の寿命)/(水準1の寿命)〕とを比較し(N25)、入力された2つのロットの寿命の寿命差の方が小さいと、有為差無しと反対してその判定結果と上記の有為差を判断できる寿命差の値(上記の例では2.57場合)とを出力画面に表示させ(N29)、コンピュータ1による処理を終了する。
【0050】
入力された2つのロットの寿命の寿命差の方が大きい場合は、次の手順N26〜N28に進み、少なくとも断定できる寿命差の演算を行う。
手順N26では、寿命分布が異なるワイブル分布との組み合わせで、上述と同様の計算を行う。具体的には、ワイブル分布と寿命が1.1 、1.2 …50倍異なるワイブル分布の組み合わせで計算を行う。このようにして計算を行うと、組み合わせが1.1 、1.2 …50倍と変化するごとに、寿命差が1.1 、1.2 …50倍を中心にしてばらつくはずである。
【0051】
最後に、組み合わせの倍率を横軸にして、それぞれの累積確率分布の90%区間を縦軸にとったグラフ(図8(D))を作成する(N27)。この図は、寿命がある倍率で異なる2ロットの試験片群から3個の試験片を1組ずつ抜き取って寿命比較を行った場合、どの程度寿命差がばらつくかを示した図である。例えば、4倍の2ロットの試験片群から2個の試験片を1組ずつ抜き取り、寿命比較を行った場合、寿命差は4倍になるケースが最も多いものの、7.5倍の寿命差になるケースも10%の確率で起こりえることがこの図から読み取れる。
【0052】
今、問題としている2ロット間の寿命差は3倍である。図から、3倍の寿命差は寿命差が1.7倍以上ある2ロット間でしか起こり得ない寿命差だということが読み取れるので、少なくともいえる寿命差は1.7倍ということができる。このように読み取った値を、少なくとも断定できる寿命差と定める(N28)。
この定めた寿命差の値は、有為差有りという表示と共に、表示装置2の出力画面2b(図7(A))に出力する(N29)。
以上が2ロットの寿命から少なくともいえる寿命差を算出する手順である。
【0053】
図4,図5に示し有為差有無判定・有為寿命差見積もりプログラム41についての上記の説明は、具体的に数値を例にとって説明したが、この有為差有無判定・有為寿命差見積もりプログラム41は、整理すると、次の手順により構成される。
【0054】
この発明の加速試験における有為差有無判定・有為寿命差見積もりプログラム(51)は、コンピュータで実行可能なプログラムであって、
上記表示装置(2)の画面に、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命の入力を促す表示を行わせる促し画面出力手順(N1)と、
実行命令に応答して上記有為差有無の判定および少なくとも断定できる寿命差の演算を行いその演算結果を上記表示装置の画面に出力する有為差判定・有為寿命差演算手順(N2)とを含む。
【0055】
この有為差判定・有為寿命差演算手順(N2)は、
上記入力情報のワイブルスロープの値および仮設定寿命から定まるワイブル分布に従った乱数を水準1のロットの試験個数分発生させる第1乱数発生手順(N21)と、
第1乱数発生手順(N21)と同じワイブル分布に従った乱数を水準2のロットの試験個数分発生させる第2乱数発生手順(N22)と、
第1乱数発生手順(N21)で発生させた設定個数の乱数から所定の寿命計算方法で求まる寿命と第2乱数発生手順(N22)で発生させた設定個数の乱数から上記所定の寿命計算方法で求まる寿命とから、1組の水準1,水準2の寿命の寿命比を算出する寿命比算出手順(N23)と、
これら第1乱数発生手順(N21),第2乱数発生手順(N22),および寿命比算出手順(N23)を設定回数繰り返し、上記設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成して累積確率分布から、設定信頼幅内の最大の寿命比を読み取る累積確率分布分析手順(N24)と、
上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手順で得た最大寿命比を超える場合に有為差有りと判定し最大寿命比以内の場合に有為差無しと判定する有為差有無判定手順(N25)とを備える。
【0056】
有為差判定・有為寿命差演算手順(N2)は、さらに、有為差有無判定手順(N25)で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手順を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す倍率変更繰り返し手順(N26)と、
この倍率変更繰り返し手順(N26)で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する寿命倍率関係グラフ作成手順(N27)と、
上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする有為寿命差読み取り手順(28)と、
上記有為差有無判定手順(N28)で判定した結果および上記有為寿命差読み取り手順で読み取った、少なくとも断定できる寿命差を表示装置に表示させる結果出力手順(29)と、を含む。
【0057】
上記累積確率分布分析手順(N24)は、手順N21〜N23を設定回数繰り返させる手順N241と、その繰り返し過程で記憶したおいた寿命比から確率分布および累積確率分布を演算する手順N242と、この累積確率分布から90%範囲の大寿命比を読み取る手順N243とからなる。
【0058】
上記倍率変更繰り返し手順(N26)は、設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す手順N261と、繰り返し毎に異なる設定倍率に変更する倍率変更手順N231とからなる。
【0059】
各乱数発生手順(N21,N22)の詳細について説明する。これらの手順N21,N22は、ワイブル分布を特定し、その特定したワイブル分布に従ってワイブル乱数を発生させる。
一般に軸受の寿命分布は次式1)のワイブル分布に従うと言われている。
【0060】
【数2】
【0061】
ただし、m:ワイブルスロープ、α:尺度因子、γ:最小寿命、
ワイブル分布は、3つのパラメータを持っており、ワイブルスロープmによって指数分布、対数正規分布、正規分布を表現できる万能分布として知られている。参考として、図10に各種パラメータを変化させた時のワイブル分布の変化を示す。ワイブルスロープmは、分布の形状を支配するパラメータであり、この値が小さいほどばらつきの大きい分布ということができる。尺度因子αは、横軸(寿命)のスケールを変化させるもので、この値が大きいほど寿命は相対的に長くなる。最小寿命γは、寿命分布の横軸(寿命)を単にシフトさせるものである。
【0062】
この実施形態では、ワイブル乱数を発生させるが、この乱数を発生させるためにはワイブル分布の3つのパラメータを決定する必要がある。決め方の手順は、例えば以下のようになる。
1) ワイブルスロープmを実績から決定する。
2) 乱数を発生させたい分布の信頼度(例えばL10寿命であるか、あるいはL50寿 命であるか)を決定する。
3) 信頼度から求めたワイブルスロープmから、最小寿命γを所定の数式を使って決定 する。例えば、L10寿命またはL50寿命から求めた尺度因子αから、
最小寿命γを、例えば、以下の2)式を使って決定する。
この式は、1990年制定のISOの最小寿命であり、実験値からの回帰式である。
【0063】
【数3】
【0064】
これは、R≦10の値で、R=0(L10寿命でのa1)のとき、この式は1になるという式である。過去のISOの最少寿命考慮の式では、L10寿命以下の寿命は、この式にL10寿命を書けた値ということで定義されている。Rは信頼度に対応する値(100−Rが信頼度となる値)である。
なお、最小寿命の定め方については、各種の規格(例えばISO)において、時代と共に変更される場合があるが、規格の変更に伴い、実施時の規格に応じた定め方を採用すれば良い。また、最小寿命は、材料試験条件によっても変化するのでより一般的な式で記述するほうが良いとの主張もあり、適宜の値を用いれば良い。
【0065】
ワイブル乱数の発生につき説明する。乱数とは、定性的にはでたらめな数列であって、発生頻度が均一(等確率)で、その発生に規則性がない(無規則性)というものであるが、完全な乱数を発生させることは不可能である。そこで、コンピュータで発生させることのできる疑似乱数を使う。簡易な乱数発生アルゴリズムでは、例えば10進法で20桁ぐらいの周期性が見られるが、周期性が6千桁以上の周期性となるものもあり、このような周期性の少ない乱数発生アルゴリズムを用いることが好ましい。
【0066】
この実施形態では、一様な乱数ではなく、ワイブル分布に従った乱数であるワイブル乱数を発生させる。このため発生方法には工夫が必要になる。確率密度関数が複雑な場合、その分布に従う乱数を発生するには棄却法と呼ばれる方法を用いればよく、この実施形態においても、棄却法を用いる。
確率密度関数f(x)の変域が図11のように、0からX0 の範囲にあるとみなされるものとし、その変域内でのf(x)の最大値をMとする。RNを区間〔0,1 〕での一様擬似乱数とするとX0 ・RNにより、区間〔0,x0〕での一様擬似乱数xiを発生することができる。同様にして、M・RNにより、区間〔0,M 〕での一様擬似乱数yiを発生することができる。そこで、このようにして発生させた乱数xi,yiがf(xi)> yi となる条件を満足する場合には、乱数xiは与えられた確率密度分布に従うものとして採用し、満足しなければ、その乱数xiを不採用とする。この作業を繰り返し、確率密度分布に従う確率で乱数xiを採用し、確率密度分布に従う乱数の数列を作っていく方法を棄却法という。この方法は、条件に合わない乱数を捨てることになるので乱数発生法としては効率がよくないが、よい一様乱数さえ得られれば原理的に正しい数列が得られる方法である。
【0067】
図5の乱数発生手順N21,N22において、ワイブル分布の特定のための寿命(L10寿命)は、適宜想定した値を、これらの手順N21,N22の計算式等に設定しておいても良く、また入力手段3からの入力によって可変としても良い。
また、寿命比較算出手順N23において、発生したワイブル乱数からのL10寿命を求める演算は、加速試験において従来から用いられている中の適宜の寿命計算方法を用いれば良い。加速試験においては、破損の発生した時間からワイブルプロットで寿命を求めることが行われており、このような破損時間から寿命を求める計算方法を、上記寿命比算出手順における寿命計算方法として用いる。
【0068】
この寿命計算には、例えば次の方法が採用できる。
(1) 寿命試験を実施する。
(2) 得られたデータ(破損した時間あるいは破損した負荷回数)を昇順に並び替える。
(3) これらデータを図13のグラフ(ワイブル確率紙)にプロットする(縦軸:累積破損確率、横軸:寿命)。
(4) 図13の紙にプロットしたデータの最適直線を最小二乗法で引く。このとき、L10 寿命以下の位置に最小寿命があるということになるので、L10 寿命の値を10分割し(何分割でも良いがフィッティングでの計算時間が妥当な時間になるように設定する)、累積確率0 %の位置にプロットを加える。10通りの最適曲線で最もデータがフィットする最適直線を採用する。
(5) そうすると、ワイブルスロープがこの線の傾き、最小寿命は、L10 寿命の値を10分割のいずれかの値、L10 寿命(ワイブルスロープが累積確率10%交わる寿命)と尺度因子αの関係からαを決定できる。
【0069】
この打切り時間見積もり装置プログラム51では、有為差のある寿命差有りと判断するための必要寿命差の関係を求めている。図9に寿命差有りと判断するための必要寿命差の関係を示す。この結果は、事前に入力するワイブルスロープと試験個数によって変化するので、図にはワイブルスロープと試験個数を変化させた結果を示している。試験個数が増えるほど、あるいはワイブルスロープが大きいほど、傾きが小さくなることが分かる。これらの結果は、(1) 試験個数が増えるほど寿命差有りと判断するための必要寿命差を小さくできること、(2) ワイブルスロープが大きい試験条件ほど寿命差有りと判断するための必要寿命差を小さくできることを示している。
【0070】
図2と共に加速試験における有為差有無判定・有為寿命差見積もり装置につきする。この有為差有無判定・有為寿命差見積もり装置は、演算処理装置1と、この演算処理装置1の出力を画面に表示する表示装置2と、上記演算処理装置1に入力を行う入力手段3とを備える。
【0071】
上記演算処理装置1は、上記表示装置2の画面に、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命の入力を促す表示を行わせる促し画面出力手段7Mと、実行命令に応答して上記有為差有無の判定および少なくとも断定できる寿命差の演算を行いその演算結果を上記表示装置2の画面に出力する有為差判定・有為寿命差演算手段51とを備える。
【0072】
上記有為差判定・有為寿命差演算手段52は、第1乱数発生手段53、第2乱数発生手段54、寿命比算出手段55、累積確率分布分析手段56、有為差有無判定手段57、倍率変更繰り返し手段58、寿命倍率関係グラフ作成手段59、有為寿命差読み取り手段60、および結果出力手段61を有する。
【0073】
第1乱数発生手段53は、上記入力情報のワイブルスロープの値および仮設定寿命から定まるワイブル分布に従った乱数を水準1のロットの試験個数分発生させる手段であり、図5の手順N21で説明した処理を行う。
第2乱数発生手段54は、第1乱数発生手段53と同じワイブル分布に従った乱数を水準2のロットの試験個数分発生させる手段であり、手順N22で説明した処理を行う。
寿命比算出手段55は、第1乱数発生手段53で発生させた設定個数の乱数から所定の寿命計算方法で求まる寿命と、第2乱数発生手段54で発生させた設定個数の乱数から上記所定の寿命計算方法で求まる寿命とから、1組の水準1,水準2の寿命の寿命比を算出する手段であり、手順N23で説明した処理を行う。
【0074】
累積確率分布分析手段56は、これら第1乱数発生手段53,第2乱数発生手段54,および寿命比算出手段55の処理を設定回数繰り返し、上記設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成して累積確率分布から、設定信頼幅内の最大の寿命比を読み取る手段であり、手順N24で説明した処理を行う。
【0075】
有為差有無判定手段57は、上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手段56で得た最大寿命比を超える場合に有為差有りと判定し最大寿命比以内の場合に有為差無しと判定する手段であり、手順N25で説明した処理を行う。
【0076】
倍率変更繰り返し手段58は、有為差有無判定手段57で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手段56の処理を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す手段であり、手順N26で説明した処理を行う。
寿命倍率関係グラフ作成手段59は、倍率変更繰り返し手段58で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する手段であり、手順N27で説明した処理を行う。
有為寿命差読み取り手段60は、上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする手段であり、手順N28で説明した処理を行う。
【0077】
結果出力手段61は、上記有為差有無判定手段57で判定した結果、および上記有為寿命差読み取り手段60で読み取った、少なくとも断定できる寿命差を表示装置2に表示させる手段であり、手順N29で説明した処理を行う。
【図面の簡単な説明】
【0078】
【図1】この発明の一実施形態に係る有為差有無判定・有為寿命差見積もり装置の概略ブロック図である。
【図2】同有為差有無判定・有為寿命差見積もり装置の概念構成を示すブロック図である。
【図3】同有為差有無判定・有為寿命差見積もり装置を用いた必要寿命差見積もり方法の概略流れ図である。
【図4】同有為差有無判定・有為寿命差見積もり方法を実施する有為差有無判定・有為寿命差見積もりプログラムの概略の流れ図である。
【図5】同プログラムにおける有為差有無判定・有為寿命差見積もり手順の詳細を示す流れ図である。
【図6】図2の有為差有無判定・有為寿命差見積もり装置における入力画面例の説明図である。
【図7】図2の有為差有無判定・有為寿命差見積もり装置における出力画面例の説明図である。
【図8】(A)はワイブル分布の例のグラフ、(B)は頻度と寿命比の関係例を示す確率分布のグラフ、(C)は累積確率と寿命比の関係例を示すグラフ、(D)は少なくとも断定できる寿命差(倍率)と寿命差との関係を示すグラフである。
【図9】(A)はワイブルスロープ別の寿命倍率と必要寿命差の関係を示すグラフ、(B)は試験個数別の寿命倍率と必要寿命差の関係を示すグラフである。
【図10】ワイブル分布の各パラメータの影響例を示すグラフである。
【図11】ワイブル分布の定め方を示すグラフである。
【図12】実施形態の方法と従来の信頼幅を使った方法の寿命ばらつきの分布比較図である。
【図13】ワイブル確率紙の説明図である。
【図14】従来の打切りおよび加速試験の手順を示す流れ図である。
【図15】従来の信頼幅を使った寿命の有為差判定方法の説明図である。
【符号の説明】
【0079】
1…コンピュータ(演算処理手段)
2…表示装置
3…入力装置
7M…促し画面出力手段
51…有為差有無判定・有為寿命差見積もりプログラム
52…有為差有無判定・有為寿命差見積もり手段
53…第1乱数発生手段
54…第2乱数発生手段
55…寿命比算出手段
56…累積分布作成手段
57…有為差有無判定手段
58…倍率変更繰り返し手段
59…結果出力手段
【技術分野】
【0001】
この発明は、軸受等の機械部品または試験片からなる試験対象品を使用環境よりも厳しい所定の環境条件におき、破損が発生した時間から寿命を算出する加速試験において、試験対象品の2つのロット間で有為差が有るか否かの判定、およびある場合にどの程度の有為差があるかを見積もる有為差有無判定・有為寿命差見積もり方法,装置、プログラムに関する。
【背景技術】
【0002】
寿命試験は、軸受等の機械部品の性能を評価するために欠かせない試験の1つである。寿命試験には、大きく分けて(1) 実機の使用環境に近い条件で試験を行う実機試験と、(2) 比較的過酷な条件で寿命試験を行う加速試験がある。前者は、製品が有限時間内に破損するケースが極めて少ないため、ある目標時間まで破損することなく試験が継続すれば、寿命は問題ないと判断する試験である(以下、このような試験を「打切り試験」と呼ぶ)。一方、後者は、比較的短時間で破損が発生するので、ワイブルプロットで寿命が算出でき(例えば非特許文献1)、その算出寿命から性能の優劣を判定する試験である(以下、このような試験を「加速試験」と呼ぶ)。
【0003】
従来より、寿命試験は経験を積んだ熟練者が行っており、試験条件や試験個数を決める寿命試験の設計と寿命試験結果の解釈に対して経験的に確からしい判断ができたと考えられる。
図14に、従来から行われてきた寿命試験の設計と寿命試験結果の解釈の手順を、打切り試験と加速試験ごとに示す。
また、現在、寿命試験において経験的に判断されているものの詳細を、表1に示す。
【0004】
【表1】
【0005】
なお、ワイブル分布を機械部品の寿命判断に用いるものは、種々の特許文献,非特許文献に提案されている。
【特許文献1】特開2006−040203号公報
【特許文献2】特開2002−277382号公報
【特許文献3】特開2005−226829号公報
【非特許文献1】真壁肇著、信頼性工学入門79、1991年発行
【発明の開示】
【発明が解決しようとする課題】
【0006】
加速試験は、上記のように破損時間から寿命を算出し、その算出寿命から性能の優劣を判定する試験であるが、2ロットの寿命試験結果から求めた算出寿命の間に、本当に有為差があるかどうかを確かめたい状況が多くある。
従来より、算出寿命の有為差判定は信頼幅という概念を使って行われてきた。しかし、この判定方法には、次のようにいくつかの問題がある。
【0007】
まず、信頼幅は、2水準間での寿命の有意差を判断できないという問題がある。その理由は、どの程度信頼幅が重なっていたら寿命に優位差があるのかということを定量的に判断できないためである。信頼幅を使って2水準間での寿命の有意差を判定する場合、次のような手順が必要になる。以下、その手順を図15を使って説明する。
1) L10寿命におけるメディアンランクの信頼水準の分布(図中Aの分布)を求める。ここで、この分布はF値表とファイサー(Fisher)の近似式を用いて求める。
2) L10寿命におけるメディアンランクの信頼水準の分布とワイブルスロープからL10のばらつき頻度の分布(図中Bの分布)を求める。
3) このばらつき分布f(x)を使って、図中の式の積分を取り、寿命倍率nが同じ寿命分布から発生する頻度を調べる。これは寿命比のばらつき分布である。
4) 3)の分布の累積確率分布を作成し、そのグラフから90%以上の確率で有意差有と判定できる寿命倍率を計算する。
【0008】
上のような手順で、従来の信頼幅を使って寿命の有意差検定や設計が可能になる。しかし、このような手順は複雑であり、実用化が難しい。このため、信頼幅を使った寿命の有意差検定が今まで行われてこなかったのは当然であるといえる。
【0009】
レオナード・ジー・ジョンソン(LEONARD G. JOHNSON)は、ワイブルスロープごと、試験個数ごとに、2水準間の平均寿命とL10 寿命の有意差を判定する方法を提案している)。しかし、その方法も、実際には使われていないのが現状である。これは、手軽な形で優位差検定を行うツールを残していないことに原因がある。
【0010】
また、従来は、有為差有りと判定を行っても、少なくともどの程度の寿命差があると言えるのかを、定量的に示す適切な方法がなかった。
【0011】
この発明の目的は、2つのロット間での有為差有無の判定を定量的に適切に行え、かつ有為差があると判定される場合に、少なくともどの程度の寿命差があると言えるのかを、簡単かつ迅速に試算することができ、かつ定量的に求められて信頼性の高いものとでき、熟練者でなくても必要な寿命差を適切に試算することのできる方法、装置、およびコンピュータプログラムを提供することである。
【課題を解決するための手段】
【0012】
この発明の加速試験における有為差有無判定・有為寿命差見積もり方法は、軸受等の機械部品または試験片からなる試験対象品を使用環境よりも厳しい所定の環境条件におき、破損が発生した時間から寿命を算出する加速試験において、試験対象品の2つのロットの寿命から、有為差有無の判定、および少なくとも断定できる、倍率による寿命差を算出する方法であって、
コンピュータに対し、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命を入力する過程(M1)と、
上記コンピュータに、上記有為差有無の判定、および少なくとも断定できる寿命差を演算させ演算結果を表示装置の画面に表示させるコンピュータ演算処理過程(M2)とを含む。
上記寿命は、例えばL10寿命(90%の信頼度の寿命)や、L50寿命(50%の信頼度の寿命)等の所定信頼度の寿命ある。
【0013】
上記コンピュータ演算処理過程(M2)として、
試験対象品の上記入力情報のワイブルスロープの値および仮設定寿命から定まるワイブル分布に従った乱数を水準1のロットの試験個数分発生させる第1乱数発生手順(N21)と、
第1乱数発生手順(N21)と同じワイブル分布に従った乱数を水準2のロットの試験個数分発生させる第2乱数発生手順(N22)と、
第1乱数発生手順(N21)で発生させた設定個数の乱数から所定の寿命計算方法で求まる寿命と第2乱数発生手順(N22)で発生させた設定個数の乱数から上記所定の寿命計算方法で求まる寿命とから、1組の水準1,水準2の寿命の寿命比を算出する寿命比算出手順(N23)と、
これら第1乱数発生手順(N21),第2乱数発生手順(N23),および寿命比算出手順(N24)を設定回数繰り返し、上記設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成して累積確率分布から、設定信頼幅内の最大の寿命比を読み取る累積確率分布分析手順(N24)と、
上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手順(N24)で得た最大寿命比を超える場合に有為差有りと判定し、最大寿命比以内の場合に有為差無しと判定する有為差有無判定手順(N25)とを含む。
【0014】
また、上記コンピュータ演算処理過程(N2)として、さらに、
上記有為差有無判定手順(N25)で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手順を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す倍率変更繰り返し手順(N26)と、
この倍率変更繰り返し手順(N26)で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する寿命倍率関係グラフ作成手順(N27)と、
上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする有為寿命差読み取り手順(N28)と、
上記有為差有無判定手順(N25)で判定した結果および上記有為寿命差読み取り手順(N28)で読み取った、少なくとも断定できる寿命差を表示装置に表示させる結果出力手順(N29)と、を含む。
【0015】
上記第1,第2乱数発生手順(N21,N22)において、乱数発生に用いるワイブル分布を特定するための寿命は、適宜想定した寿命または実績値を設定すれば良い。
また、寿命比算出手順(N23)において、乱数発生手順で発生させた試験個数の乱数から寿命を求める所定の寿命計算方法は、加速試験において従来から用いられている中の適宜の寿命計算方法を用いれば良い。加速試験においては、破損の発生した時間からワイブルプロットで寿命を求めることが行われており、このような破損時間から寿命を求める計算方法を、上記寿命比算出手順における寿命計算方法として用いる。
ワイブル分布は、次式、
【0016】
【数1】
【0017】
ただし、m:ワイブルスロープ、α:尺度因子、γ:最小寿命、
によって特定される。
【0018】
軸受等の機械部品の寿命は、ワイブル分布に従うとされている。ワイブル分布は、ワイブルスロープm、尺度因子α、最小寿命γの3つのパラメータを持っており、ワイブルスロープmによって指数分布、対数正規分布、正規分布を表現できる万能分布として知られている。量産される軸受等では、ワイブルスロープは実績値が既知である場合が多く、この発明方法において、ワイブルスロープには、試験対象品の実績値を用いることが好ましい。実績値がない場合は、適宜の方法で見積もったワイブルスロープを用いてもよい。最小寿命γは、種々の規格、例えばISO等によって計算方法が定められており、そのように定められたいずれかの計算方法を用いることが好ましい。尺度因子αは、ワイブルスロープの値、要求寿命の信頼度、要求寿命の値、および上記最小寿命γから一義的に決定される演算式があり、その演算式を用いて特定しても良い。
【0019】
この方法において、累積分布作成手順(N24)で行う処理は、同じ寿命分布を持つロットから験個数分の試験対象品を2組抜きとり、寿命試験を行った結果、どの程度の寿命差が出るのかを設定回数確認することに対応している。同じ寿命分布を持つロットから試験対象品を抜きとった試験においても、寿命はばらついてる。その寿命比のばらつきを累積確率分布の適宜の信頼区間、例えば5%と95%の区間(90%信頼区間)として定義すると、その寿命ばらつきは、上記累積確率分布から算出できる。この結果は、同じ寿命分布を持つロットから試験個数の試験対象品を2組抜きとって寿命比較を行っても、90%のものは、累積確率分布に示された信頼区間の範囲で寿命比がばらついてしまうことを示している。 逆にいうと上記信頼区間の最大の寿命差があれば、そのロット間の寿命分布が異なっている可能性が高いということになる。
したがって、上記有為差有無判定手順(N25)として、上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手順で得た最大寿命比を超える場合に有為差有りと判定し、最大寿命比以内の場合に有為差無しと判定することで、定量的に適切に有為差有無を判定することができる。
【0020】
また、寿命が異なるワイブル分布との組み合わせで、上記同様の累積確率分布の計算を行い、例えば、ワイブル分布と寿命が1.1 、1.2 …50倍異なるワイブル分布の組み合わせで計算を行うと、その組み合わせが1.1 、1.2 …50倍と変化するごとに、寿命差が1.1 、1.2 …50倍を中心にしてばらつくはずである。上記倍率変更繰り返し手順(N26)は、このような累積確率分布の計算、およびその確率分布についての上記確率分布分析手順(N24)を行う。
この倍率変更繰り返し手順(N26)で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成し、このグラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取ると、その読み取った値が少なくとも断定できる寿命差となる。
【0021】
この発明方法は、このようにワイブル分布に基づく乱数シミュレーションにより有為差有無を判定し、また少なくとも断定できる寿命差を求めるため、これら有為差有無および少なくとも断定できる寿命差を、簡単かつ迅速に試算することができる。また、定量的に求められて信頼性の高いものとでき、熟練者でなくても試算することができる。この乱数シミュレーションに用いるプログラムは、条件入力を行えば勝手に有為差有無を判定および少なくとも断定できる寿命差の検定を行ってくれる手軽なツールであることが優れた点である。
【0022】
この発明の加速試験における有為差有無判定・有為寿命差見積もり装置は、軸受等の機械部品または試験片からなる試験対象品を使用環境よりも厳しい所定の環境条件におき、破損が発生した時間から寿命を算出する加速試験において、試験対象品の2つのロットの寿命から、有為差有無の判定、および少なくとも断定できる、倍率による寿命差を算出する装置であって、
演算処理装置(1)と、この演算処理装置(1)の出力を画面に表示する表示装置(2)と、上記演算処理装置(1)に入力を行う入力手段(3)とを備える。
【0023】
上記演算処理装置(1)は、上記表示装置(2)の画面に、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命の入力を促す表示を行わせる促し画面出力手段(7M)と、実行命令に応答して上記有為差有無の判定および少なくとも断定できる寿命差の演算を行いその演算結果を上記表示装置(2)の画面に出力する有為差判定・有為寿命差演算手段(52)とを備える。
【0024】
上記有為差判定・有為寿命差演算手段(52)は、第1乱数発生手段(53)、第2乱数発生手段(54)、寿命比算出手段(55)、累積確率分布分析手段(56)、有為差有無判定手段(57)、倍率変更繰り返し手段(58)、寿命倍率関係グラフ作成手段(59)、有為寿命差読み取り手段(60)、および結果出力手段(61)を有する。
【0025】
第1乱数発生手段(53)は、上記入力情報のワイブルスロープの値および仮設定寿命から定まるワイブル分布に従った乱数を水準1のロットの試験個数分発生させる手段である。
第2乱数発生手段(54)は、第1乱数発生手段(53)と同じワイブル分布に従った乱数を水準2のロットの試験個数分発生させる手段である。
寿命比算出手段(55)は、第1乱数発生手段(53)で発生させた設定個数の乱数から所定の寿命計算方法で求まる寿命と、第2乱数発生手段(54)で発生させた設定個数の乱数から上記所定の寿命計算方法で求まる寿命とから、1組の水準1,水準2の寿命の寿命比を算出する手段である。
【0026】
累積確率分布分析手段(56)は、これら第1乱数発生手段(53),第2乱数発生手段(54),および寿命比算出手段(55)の処理を設定回数繰り返し、上記設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成して累積確率分布から、設定信頼幅内の最大の寿命比を読み取る手段である。
【0027】
有為差有無判定手段(57)は、上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手段(56)で得た最大寿命比を超える場合に有為差有りと判定し、最大寿命比以内の場合に有為差無しと判定する手段である。
【0028】
倍率変更繰り返し手段(58)は、有為差有無判定手段(57)で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手段(56)の処理を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す手段である。
寿命倍率関係グラフ作成手段(59)は、倍率変更繰り返し手段(58)で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する手段である。
有為寿命差読み取り手段(60)は、上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする手段である。
【0029】
結果出力手段(61)は、上記有為差有無判定手段(57)で判定した結果、および上記有為寿命差読み取り手段(60)で読み取った、少なくとも断定できる寿命差を表示装置(2)に表示させる手段である。
【0030】
この構成の加速試験における有為差有無判定・有為寿命差見積もり装置によると、この発明の有為差有無判定・有為寿命差見積もり方法に使用されて、有為差の有無の判定、および有為差があると判定される場合に、少なくともどの程度の寿命差があると言えるのかを、簡単かつ迅速に試算することができ、かつ定量的に求められて信頼性の高いものとでき、熟練者でなくても少なくとも断定できる寿命差を適切に試算することができる。
【0031】
この発明の加速試験における有為差有無判定・有為寿命差見積もりプログラム(51)は、コンピュータで実行可能なプログラムであって、
上記表示装置(2)の画面に、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命の入力を促す表示を行わせる促し画面出力手順(N1)と、
実行命令に応答して上記有為差有無の判定および少なくとも断定できる寿命差の演算を行いその演算結果を上記表示装置の画面に出力する有為差判定・有為寿命差演算手順(N2)とを含む。
【0032】
この有為差判定・有為寿命差演算手順(N2)は、
上記入力情報のワイブルスロープの値および仮設定寿命から定まるワイブル分布に従った乱数を水準1のロットの試験個数分発生させる第1乱数発生手順(N21)と、
第1乱数発生手順(N21)と同じワイブル分布に従った乱数を水準2のロットの試験個数分発生させる第2乱数発生手順(N22)と、
第1乱数発生手順(N21)で発生させた設定個数の乱数から所定の寿命計算方法で求まる寿命と第2乱数発生手順(N22)で発生させた設定個数の乱数から上記所定の寿命計算方法で求まる寿命とから、1組の水準1,水準2の寿命の寿命比を算出する寿命比算出手順(N23)と、
これら第1乱数発生手順(N21),第2乱数発生手順(N22),および寿命比算出手順(N23)を設定回数繰り返し、上記設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成して累積確率分布から、設定信頼幅内の最大の寿命比を読み取る累積確率分布分析手順(N24)と、
上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手順で得た最大寿命比を超える場合に有為差有りと判定し、最大寿命比以内の場合に有為差無しと判定する有為差有無判定手順(N25)とを備える。
【0033】
有為差判定・有為寿命差演算手順(N2)は、さらに、有為差有無判定手順(N25)で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手順を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す倍率変更繰り返し手順(N26)と、
この倍率変更繰り返し手順(N26)で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する寿命倍率関係グラフ作成手順(N27)と、
上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする有為寿命差読み取り手順(28)と、
上記有為差有無判定手順で判定した結果および上記有為寿命差読み取り手順で読み取った、少なくとも断定できる寿命差を表示装置に表示させる結果出力手順(29)と、
を含む。
【0034】
この構成の加速試験における有為差有無判定・有為寿命差見積もりプログラムは、この発明の有為差有無判定・有為寿命差見積もりの実施に使用され、有為差の有無の判定、および有為差があると判定される場合に、少なくともどの程度の寿命差があると言えるのかを、簡単かつ迅速に試算することができ、かつ定量的に求められて信頼性の高いものとでき、熟練者でなくても必要な寿命差を適切に試算することができる。
【発明の効果】
【0035】
この発明の加速試験における有為差有無判定・有為寿命差見積もり方法、装置、およびプログラムによると、コンピュータシミュレーションにより、試験対象品のワイブル分布から乱数を発生させ、同じワイブル分布から乱数を発生させ、その発生させた乱数から求まる寿命を演算して1組の寿命比を算出し、この処理を設定回数繰り返して設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成し、上記累積確率分布から設定信頼幅内の最大および最小の寿命比を読み取り、有為差有無の判定を行うため、有為差有無の判定を、簡単かつ迅速に行うことができ、かつ定量的に求められて信頼性の高いものとでき、熟練者でなくても判定することができる。
また、有為差有りと判定した場合に、寿命が異なるワイブル分布に繰り返し毎に変えて、上記累積確率分布分析手順を繰り返し、繰り返しにより得られた設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成し、このグラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とするため、有為差がある場合に、少なくとも断定できる寿命差を、簡単かつ迅速に求めることができ、かつ定量的に求められて信頼性の高いものとでき、熟練者でなくても適切に有為寿命差を求めることができる。
【発明を実施するための最良の形態】
【0036】
この発明の実施形態を説明する。この加速試験における有為差有無判定・有為寿命差見積もり方法は、軸受等の機械部品または試験片からなる試験対象品を使用環境よりも厳しい所定の環境条件におき、破損が発生した時間から寿命を算出する加速試験において、試験対象品の2つのロットの寿命から、有為差有無の判定、および少なくとも断定できる、倍率による寿命差を算出する方法である。
【0037】
以下、この実施形態を図面と共に説明する。この加速試験における有為差有無判定・有為寿命差見積もり方法は、図1に示すコンピュータ1に、乱数シミュレーションプログラムである加速試験における有為差有無判定・有為寿命差見積もりプログラム41を実行させることで行う。コンピュータ1はパーソナルコンピュータ等からなり、中央処理装置4およびメモリ5を有し、所定のオペレーションシステムによって動作するものである。コンピュータ1には、液晶表示装置等の画面によって表示可能な表示装置2と、キーボードやマウス等の入力装置3が接続され、あるいは付属して設けられている。コンピュータ1、表示装置2、入力装置3、および必要寿命差見積もりプログラム41により、図2に各機能達成手段をブロックで示した有為差有無判定・有為寿命差見積もり装置が構成される。同図の有為差有無判定・有為寿命差見積もり装置の構成については、後に説明する。
有為差有無判定・有為寿命差見積もりプログラム51は、図4および図5に流れ図で示す手順を備えるものである。同図の内容は、後に説明する。
【0038】
この必要寿命差見積もり方法は、図3に示すように、コンピュータ1に対して所定の情報を入力する入力過程M1と、コンピュータ1で演算処理を行って演算結果を出力するコンピュータ演算処理過程M2とからなる。
【0039】
入力過程M1では、図6に示すように所定の入力情報の入力を促す入力画面2aが、コンピュータ1の出力によって表示装置2に表示され、入力画面2a中に、所定の入力を促す表示が行われる。
入力画面2aは、同図の例では複数の画面に切り替えて表示される。入力画面のうち、第1の画面(図6(A)では、L10寿命とL50寿命とのいずれで有為差を検定するかを選択させる表示が行われる。この画面に対して、ラジオボタン等からなる表示部分を選択することで、L10寿命とL50寿命とのいずれかが選択できる。
第2の画面(図6(B))では、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値の入力を促す表示、および入力を記入する記入用表示が行われる。オペレータは、この記入用表示に、該当する値を入力する。
第2の画面(図6(C))では、比較する水準1および水準2の2つのロットの各試験個数n1,n2の入力、および試験結果の寿命を促す表示が行われ、その個数および寿命を入力する記入用表示が行われる。オペレータは、この記入用表示に対して該当する情報を入力する。
【0040】
ワイブルスロープの値には、試験の実績値を入力することが好ましい。実績値は10個以上の試験で得た結果を用いることが望ましく、より好ましくは20個以上の試験結果である。また、試験個数および試験結果は、実際に行った試験における個数および結果を入力する。なお、試験個数と寿命は、水準2の欄に長寿命側の入力をしなければないらいとする入力規則が、有為差有無判定・有為寿命差見積もりプログラム51に手順として設定されている。この手順は、例えば促し画面出力手順N1の一部として設けられる。以上の点に注意して、所定の条件を入力し、入力画面2aのOKボタンを押すと計算が開始される。
【0041】
図3のコンピュータ演算処理過程M2では、入力されたワイブルスロープの値、および2つのロットの各試験個数n1,n2、試験結果の寿命から、有為差有無の判定に必要な寿命差を計算して、有為差有無を判定し、有為差有りの場合には、少なくとも断定できる寿命差を計算する。
【0042】
計算が終了すると、図7に示す出力画面に結果が表示される。有為差がない場合は、同図(A)のように、有為差無しという表示と、有為差有りと判断するための必要寿命差が表示される。
有為差がある場合は、同図(B)のように、有為差ありという表示と、少なくとも断定できる寿命差(倍率)が表示される。
【0043】
図1の有為差有無判定・有為寿命差見積もりプログラム51は、コンピュータ1で実行可能なプログラムであって、図4,図5に流れ図で示す手順を備える。図4に示すように、有為差有無判定・有為寿命差見積もりプログラム51は、促し画面出力手順N1と必要寿命差演算手順N2とでなり、促し画面出力手順N1では、図6と共に前述した入力画面2aを出力する。この入力画面2aに対して、上記各入力情報が入力手段3から入力され、かつ入力画面2aのOKキーのクリック等によって実行命令が入力手段3から入力されると、有為差有無判定・有為寿命差見積もり手順N2が実行される。同図の入力画面に対して入力する過程が、図3の入力過程M1であり、同図のコンピュータ演算処理過程M2は図4の有為差有無判定・有為寿命差見積もり手順N2を実行する過程である。
【0044】
有為差有無判定・有為寿命差見積もり手順N2は、図5に流れ図で示す各手順で構成される。この流れ図には各手順M21〜M29毎の具体的な処理例を注釈として併記してある。
理解の容易のため、具体的数値を上げ、同図の具体的処理例を参照して、各手順N21〜N29を説明する。
【0045】
今、あるL10寿命を持つワイブル分布(図8(A))で、ワイブルスロープは4としている)である加速試験で、2ロットの試験を行ったとする。1ロット目(水準1)の試験個数n1が3個で、2ロット目(水準2)の試験個数n2も3個であるとする。
まず、そのワイブル分布に従って1ロット目の試験個数n1である3個の乱数を発生させ、その3個のデータからL10寿命とL50寿命を算出する(N21)。算出方法は、従来の加速試験の結果からの寿命算出に用いられている適宜の方法を使用する。なお、L10寿命だけ、またはL50寿命だけを演算しても良いが、この実施形態では、L10寿命とL50寿命のそれぞれについて必要寿命差を計算するために、L10寿命とL50寿命の両方を算出している。
次に、同じL10寿命を持つワイブル分布から、第2ロット目の試験個数n2である3個の乱数を3個発生させ、その3個のデータからL10寿命とL50寿命を上記と同じ方法で算出する(N22)。
【0046】
次に、得られた1組のL10寿命とL50寿命の寿命比をそれぞれ算出する。すなわちL10寿命同士の寿命比、およびL50寿命同士の寿命比を算出する(N23)。
【0047】
次に、これらの作業を設定回数(例えば1000回)繰り返し、設定回数組(1000組)のL10およびL50寿命の寿命比の確率分布(図8(B)と累積確率分布(図8(C)を作成する(N24)。
これらは、同じ寿命分布を持つロットから3個の試験片を2組抜きとり、寿命試験を行った結果、どの程度の寿命差が出るのかを設定回数(1000回)確認することに対応している。
【0048】
度数分布の図(図8(B))から、同じ寿命分布を持つロットから試験片を抜きとって試験を行ったにもかかわらず、寿命はばらついていることが分かる。その寿命比のばらつきを、所定信頼度区間、例えば累積確率分布の5%と95%の区間(90%信頼区間)として定義すると、その寿命ばらつきは、図8(C)に示したように0.39〜2.57倍であると算出できる(N243)。
この結果は、同じ寿命分布を持つロットから3個の試験片を2組抜きとって寿命比較を行っても、90%のものは0.39〜2.57倍の間で寿命比がばらついてしまうことを示している。逆にいうと、ばらつき内の最大値である2.57倍以上の寿命差があれば、そのロット間の寿命分布が異なっている可能性が高いということになるので、2ロットの比較試験で試験個数がそれぞれ3個の場合、有為差を判断できる寿命差は2.57倍以上であるといえる。
【0049】
そこで、上記ばらつき内の最大値である2.57倍と、入力された2つのロット(水準1,水準2)の寿命差(倍率)、すなわち〔(水準2の寿命)/(水準1の寿命)〕とを比較し(N25)、入力された2つのロットの寿命の寿命差の方が小さいと、有為差無しと反対してその判定結果と上記の有為差を判断できる寿命差の値(上記の例では2.57場合)とを出力画面に表示させ(N29)、コンピュータ1による処理を終了する。
【0050】
入力された2つのロットの寿命の寿命差の方が大きい場合は、次の手順N26〜N28に進み、少なくとも断定できる寿命差の演算を行う。
手順N26では、寿命分布が異なるワイブル分布との組み合わせで、上述と同様の計算を行う。具体的には、ワイブル分布と寿命が1.1 、1.2 …50倍異なるワイブル分布の組み合わせで計算を行う。このようにして計算を行うと、組み合わせが1.1 、1.2 …50倍と変化するごとに、寿命差が1.1 、1.2 …50倍を中心にしてばらつくはずである。
【0051】
最後に、組み合わせの倍率を横軸にして、それぞれの累積確率分布の90%区間を縦軸にとったグラフ(図8(D))を作成する(N27)。この図は、寿命がある倍率で異なる2ロットの試験片群から3個の試験片を1組ずつ抜き取って寿命比較を行った場合、どの程度寿命差がばらつくかを示した図である。例えば、4倍の2ロットの試験片群から2個の試験片を1組ずつ抜き取り、寿命比較を行った場合、寿命差は4倍になるケースが最も多いものの、7.5倍の寿命差になるケースも10%の確率で起こりえることがこの図から読み取れる。
【0052】
今、問題としている2ロット間の寿命差は3倍である。図から、3倍の寿命差は寿命差が1.7倍以上ある2ロット間でしか起こり得ない寿命差だということが読み取れるので、少なくともいえる寿命差は1.7倍ということができる。このように読み取った値を、少なくとも断定できる寿命差と定める(N28)。
この定めた寿命差の値は、有為差有りという表示と共に、表示装置2の出力画面2b(図7(A))に出力する(N29)。
以上が2ロットの寿命から少なくともいえる寿命差を算出する手順である。
【0053】
図4,図5に示し有為差有無判定・有為寿命差見積もりプログラム41についての上記の説明は、具体的に数値を例にとって説明したが、この有為差有無判定・有為寿命差見積もりプログラム41は、整理すると、次の手順により構成される。
【0054】
この発明の加速試験における有為差有無判定・有為寿命差見積もりプログラム(51)は、コンピュータで実行可能なプログラムであって、
上記表示装置(2)の画面に、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命の入力を促す表示を行わせる促し画面出力手順(N1)と、
実行命令に応答して上記有為差有無の判定および少なくとも断定できる寿命差の演算を行いその演算結果を上記表示装置の画面に出力する有為差判定・有為寿命差演算手順(N2)とを含む。
【0055】
この有為差判定・有為寿命差演算手順(N2)は、
上記入力情報のワイブルスロープの値および仮設定寿命から定まるワイブル分布に従った乱数を水準1のロットの試験個数分発生させる第1乱数発生手順(N21)と、
第1乱数発生手順(N21)と同じワイブル分布に従った乱数を水準2のロットの試験個数分発生させる第2乱数発生手順(N22)と、
第1乱数発生手順(N21)で発生させた設定個数の乱数から所定の寿命計算方法で求まる寿命と第2乱数発生手順(N22)で発生させた設定個数の乱数から上記所定の寿命計算方法で求まる寿命とから、1組の水準1,水準2の寿命の寿命比を算出する寿命比算出手順(N23)と、
これら第1乱数発生手順(N21),第2乱数発生手順(N22),および寿命比算出手順(N23)を設定回数繰り返し、上記設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成して累積確率分布から、設定信頼幅内の最大の寿命比を読み取る累積確率分布分析手順(N24)と、
上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手順で得た最大寿命比を超える場合に有為差有りと判定し最大寿命比以内の場合に有為差無しと判定する有為差有無判定手順(N25)とを備える。
【0056】
有為差判定・有為寿命差演算手順(N2)は、さらに、有為差有無判定手順(N25)で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手順を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す倍率変更繰り返し手順(N26)と、
この倍率変更繰り返し手順(N26)で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する寿命倍率関係グラフ作成手順(N27)と、
上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする有為寿命差読み取り手順(28)と、
上記有為差有無判定手順(N28)で判定した結果および上記有為寿命差読み取り手順で読み取った、少なくとも断定できる寿命差を表示装置に表示させる結果出力手順(29)と、を含む。
【0057】
上記累積確率分布分析手順(N24)は、手順N21〜N23を設定回数繰り返させる手順N241と、その繰り返し過程で記憶したおいた寿命比から確率分布および累積確率分布を演算する手順N242と、この累積確率分布から90%範囲の大寿命比を読み取る手順N243とからなる。
【0058】
上記倍率変更繰り返し手順(N26)は、設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す手順N261と、繰り返し毎に異なる設定倍率に変更する倍率変更手順N231とからなる。
【0059】
各乱数発生手順(N21,N22)の詳細について説明する。これらの手順N21,N22は、ワイブル分布を特定し、その特定したワイブル分布に従ってワイブル乱数を発生させる。
一般に軸受の寿命分布は次式1)のワイブル分布に従うと言われている。
【0060】
【数2】
【0061】
ただし、m:ワイブルスロープ、α:尺度因子、γ:最小寿命、
ワイブル分布は、3つのパラメータを持っており、ワイブルスロープmによって指数分布、対数正規分布、正規分布を表現できる万能分布として知られている。参考として、図10に各種パラメータを変化させた時のワイブル分布の変化を示す。ワイブルスロープmは、分布の形状を支配するパラメータであり、この値が小さいほどばらつきの大きい分布ということができる。尺度因子αは、横軸(寿命)のスケールを変化させるもので、この値が大きいほど寿命は相対的に長くなる。最小寿命γは、寿命分布の横軸(寿命)を単にシフトさせるものである。
【0062】
この実施形態では、ワイブル乱数を発生させるが、この乱数を発生させるためにはワイブル分布の3つのパラメータを決定する必要がある。決め方の手順は、例えば以下のようになる。
1) ワイブルスロープmを実績から決定する。
2) 乱数を発生させたい分布の信頼度(例えばL10寿命であるか、あるいはL50寿 命であるか)を決定する。
3) 信頼度から求めたワイブルスロープmから、最小寿命γを所定の数式を使って決定 する。例えば、L10寿命またはL50寿命から求めた尺度因子αから、
最小寿命γを、例えば、以下の2)式を使って決定する。
この式は、1990年制定のISOの最小寿命であり、実験値からの回帰式である。
【0063】
【数3】
【0064】
これは、R≦10の値で、R=0(L10寿命でのa1)のとき、この式は1になるという式である。過去のISOの最少寿命考慮の式では、L10寿命以下の寿命は、この式にL10寿命を書けた値ということで定義されている。Rは信頼度に対応する値(100−Rが信頼度となる値)である。
なお、最小寿命の定め方については、各種の規格(例えばISO)において、時代と共に変更される場合があるが、規格の変更に伴い、実施時の規格に応じた定め方を採用すれば良い。また、最小寿命は、材料試験条件によっても変化するのでより一般的な式で記述するほうが良いとの主張もあり、適宜の値を用いれば良い。
【0065】
ワイブル乱数の発生につき説明する。乱数とは、定性的にはでたらめな数列であって、発生頻度が均一(等確率)で、その発生に規則性がない(無規則性)というものであるが、完全な乱数を発生させることは不可能である。そこで、コンピュータで発生させることのできる疑似乱数を使う。簡易な乱数発生アルゴリズムでは、例えば10進法で20桁ぐらいの周期性が見られるが、周期性が6千桁以上の周期性となるものもあり、このような周期性の少ない乱数発生アルゴリズムを用いることが好ましい。
【0066】
この実施形態では、一様な乱数ではなく、ワイブル分布に従った乱数であるワイブル乱数を発生させる。このため発生方法には工夫が必要になる。確率密度関数が複雑な場合、その分布に従う乱数を発生するには棄却法と呼ばれる方法を用いればよく、この実施形態においても、棄却法を用いる。
確率密度関数f(x)の変域が図11のように、0からX0 の範囲にあるとみなされるものとし、その変域内でのf(x)の最大値をMとする。RNを区間〔0,1 〕での一様擬似乱数とするとX0 ・RNにより、区間〔0,x0〕での一様擬似乱数xiを発生することができる。同様にして、M・RNにより、区間〔0,M 〕での一様擬似乱数yiを発生することができる。そこで、このようにして発生させた乱数xi,yiがf(xi)> yi となる条件を満足する場合には、乱数xiは与えられた確率密度分布に従うものとして採用し、満足しなければ、その乱数xiを不採用とする。この作業を繰り返し、確率密度分布に従う確率で乱数xiを採用し、確率密度分布に従う乱数の数列を作っていく方法を棄却法という。この方法は、条件に合わない乱数を捨てることになるので乱数発生法としては効率がよくないが、よい一様乱数さえ得られれば原理的に正しい数列が得られる方法である。
【0067】
図5の乱数発生手順N21,N22において、ワイブル分布の特定のための寿命(L10寿命)は、適宜想定した値を、これらの手順N21,N22の計算式等に設定しておいても良く、また入力手段3からの入力によって可変としても良い。
また、寿命比較算出手順N23において、発生したワイブル乱数からのL10寿命を求める演算は、加速試験において従来から用いられている中の適宜の寿命計算方法を用いれば良い。加速試験においては、破損の発生した時間からワイブルプロットで寿命を求めることが行われており、このような破損時間から寿命を求める計算方法を、上記寿命比算出手順における寿命計算方法として用いる。
【0068】
この寿命計算には、例えば次の方法が採用できる。
(1) 寿命試験を実施する。
(2) 得られたデータ(破損した時間あるいは破損した負荷回数)を昇順に並び替える。
(3) これらデータを図13のグラフ(ワイブル確率紙)にプロットする(縦軸:累積破損確率、横軸:寿命)。
(4) 図13の紙にプロットしたデータの最適直線を最小二乗法で引く。このとき、L10 寿命以下の位置に最小寿命があるということになるので、L10 寿命の値を10分割し(何分割でも良いがフィッティングでの計算時間が妥当な時間になるように設定する)、累積確率0 %の位置にプロットを加える。10通りの最適曲線で最もデータがフィットする最適直線を採用する。
(5) そうすると、ワイブルスロープがこの線の傾き、最小寿命は、L10 寿命の値を10分割のいずれかの値、L10 寿命(ワイブルスロープが累積確率10%交わる寿命)と尺度因子αの関係からαを決定できる。
【0069】
この打切り時間見積もり装置プログラム51では、有為差のある寿命差有りと判断するための必要寿命差の関係を求めている。図9に寿命差有りと判断するための必要寿命差の関係を示す。この結果は、事前に入力するワイブルスロープと試験個数によって変化するので、図にはワイブルスロープと試験個数を変化させた結果を示している。試験個数が増えるほど、あるいはワイブルスロープが大きいほど、傾きが小さくなることが分かる。これらの結果は、(1) 試験個数が増えるほど寿命差有りと判断するための必要寿命差を小さくできること、(2) ワイブルスロープが大きい試験条件ほど寿命差有りと判断するための必要寿命差を小さくできることを示している。
【0070】
図2と共に加速試験における有為差有無判定・有為寿命差見積もり装置につきする。この有為差有無判定・有為寿命差見積もり装置は、演算処理装置1と、この演算処理装置1の出力を画面に表示する表示装置2と、上記演算処理装置1に入力を行う入力手段3とを備える。
【0071】
上記演算処理装置1は、上記表示装置2の画面に、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命の入力を促す表示を行わせる促し画面出力手段7Mと、実行命令に応答して上記有為差有無の判定および少なくとも断定できる寿命差の演算を行いその演算結果を上記表示装置2の画面に出力する有為差判定・有為寿命差演算手段51とを備える。
【0072】
上記有為差判定・有為寿命差演算手段52は、第1乱数発生手段53、第2乱数発生手段54、寿命比算出手段55、累積確率分布分析手段56、有為差有無判定手段57、倍率変更繰り返し手段58、寿命倍率関係グラフ作成手段59、有為寿命差読み取り手段60、および結果出力手段61を有する。
【0073】
第1乱数発生手段53は、上記入力情報のワイブルスロープの値および仮設定寿命から定まるワイブル分布に従った乱数を水準1のロットの試験個数分発生させる手段であり、図5の手順N21で説明した処理を行う。
第2乱数発生手段54は、第1乱数発生手段53と同じワイブル分布に従った乱数を水準2のロットの試験個数分発生させる手段であり、手順N22で説明した処理を行う。
寿命比算出手段55は、第1乱数発生手段53で発生させた設定個数の乱数から所定の寿命計算方法で求まる寿命と、第2乱数発生手段54で発生させた設定個数の乱数から上記所定の寿命計算方法で求まる寿命とから、1組の水準1,水準2の寿命の寿命比を算出する手段であり、手順N23で説明した処理を行う。
【0074】
累積確率分布分析手段56は、これら第1乱数発生手段53,第2乱数発生手段54,および寿命比算出手段55の処理を設定回数繰り返し、上記設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成して累積確率分布から、設定信頼幅内の最大の寿命比を読み取る手段であり、手順N24で説明した処理を行う。
【0075】
有為差有無判定手段57は、上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手段56で得た最大寿命比を超える場合に有為差有りと判定し最大寿命比以内の場合に有為差無しと判定する手段であり、手順N25で説明した処理を行う。
【0076】
倍率変更繰り返し手段58は、有為差有無判定手段57で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手段56の処理を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す手段であり、手順N26で説明した処理を行う。
寿命倍率関係グラフ作成手段59は、倍率変更繰り返し手段58で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する手段であり、手順N27で説明した処理を行う。
有為寿命差読み取り手段60は、上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする手段であり、手順N28で説明した処理を行う。
【0077】
結果出力手段61は、上記有為差有無判定手段57で判定した結果、および上記有為寿命差読み取り手段60で読み取った、少なくとも断定できる寿命差を表示装置2に表示させる手段であり、手順N29で説明した処理を行う。
【図面の簡単な説明】
【0078】
【図1】この発明の一実施形態に係る有為差有無判定・有為寿命差見積もり装置の概略ブロック図である。
【図2】同有為差有無判定・有為寿命差見積もり装置の概念構成を示すブロック図である。
【図3】同有為差有無判定・有為寿命差見積もり装置を用いた必要寿命差見積もり方法の概略流れ図である。
【図4】同有為差有無判定・有為寿命差見積もり方法を実施する有為差有無判定・有為寿命差見積もりプログラムの概略の流れ図である。
【図5】同プログラムにおける有為差有無判定・有為寿命差見積もり手順の詳細を示す流れ図である。
【図6】図2の有為差有無判定・有為寿命差見積もり装置における入力画面例の説明図である。
【図7】図2の有為差有無判定・有為寿命差見積もり装置における出力画面例の説明図である。
【図8】(A)はワイブル分布の例のグラフ、(B)は頻度と寿命比の関係例を示す確率分布のグラフ、(C)は累積確率と寿命比の関係例を示すグラフ、(D)は少なくとも断定できる寿命差(倍率)と寿命差との関係を示すグラフである。
【図9】(A)はワイブルスロープ別の寿命倍率と必要寿命差の関係を示すグラフ、(B)は試験個数別の寿命倍率と必要寿命差の関係を示すグラフである。
【図10】ワイブル分布の各パラメータの影響例を示すグラフである。
【図11】ワイブル分布の定め方を示すグラフである。
【図12】実施形態の方法と従来の信頼幅を使った方法の寿命ばらつきの分布比較図である。
【図13】ワイブル確率紙の説明図である。
【図14】従来の打切りおよび加速試験の手順を示す流れ図である。
【図15】従来の信頼幅を使った寿命の有為差判定方法の説明図である。
【符号の説明】
【0079】
1…コンピュータ(演算処理手段)
2…表示装置
3…入力装置
7M…促し画面出力手段
51…有為差有無判定・有為寿命差見積もりプログラム
52…有為差有無判定・有為寿命差見積もり手段
53…第1乱数発生手段
54…第2乱数発生手段
55…寿命比算出手段
56…累積分布作成手段
57…有為差有無判定手段
58…倍率変更繰り返し手段
59…結果出力手段
【特許請求の範囲】
【請求項1】
軸受等の機械部品または試験片からなる試験対象品を使用環境よりも厳しい所定の環境条件におき、破損が発生した時間から寿命を算出する加速試験において、試験対象品の2つのロットの寿命から、有為差有無の判定、および少なくとも断定できる、倍率による寿命差を算出する方法であって、
コンピュータに対し、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命を入力する過程と、
上記コンピュータに、上記有為差有無の判定、および少なくとも断定できる寿命差を演算させ演算結果を表示装置の画面に表示させるコンピュータ演算処理過程とを含み、
上記コンピュータ演算処理過程として、
試験対象品の上記入力情報のワイブルスロープの値および仮設定寿命から定まるワイブル分布に従った乱数を水準1のロットの試験個数分発生させる第1乱数発生手順と、
第1乱数発生手順と同じワイブル分布に従った乱数を水準2のロットの試験個数分発生させる第2乱数発生手順と、
第1乱数発生手順で発生させた設定個数の乱数から所定の寿命計算方法で求まる寿命と第2乱数発生手順で発生させた設定個数の乱数から上記所定の寿命計算方法で求まる寿命とから、1組の水準1,水準2の寿命の寿命比を算出する寿命比算出手順と、
これら第1乱数発生手順,第2乱数発生手順,および寿命比算出手順を設定回数繰り返し、上記設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成して累積確率分布から、設定信頼幅内の最大の寿命比を読み取る累積確率分布分析手順と、 上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手順で得た最大寿命比を超える場合に有為差有りと判定し、最大寿命比以内の場合に有為差無しと判定する有為差有無判定手順と、
この手順で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手順を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す倍率変更繰り返し手順と、
この倍率変更繰り返し手順で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する寿命倍率関係グラフ作成手順と、
上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする有為寿命差読み取り手順と、
上記有為差有無判定手順で判定した結果および上記有為寿命差読み取り手順で読み取った、少なくとも断定できる寿命差を表示装置に表示させる結果出力手順と、
を含む、加速試験における有為差有無判定・有為寿命差見積もり方法。
【請求項2】
軸受等の機械部品または試験片からなる試験対象品を使用環境よりも厳しい所定の環境条件におき、破損が発生した時間から寿命を算出する加速試験において、試験対象品の2つのロットの寿命から、有為差有無の判定、および少なくとも断定できる、倍率による寿命差を算出する装置であって、
演算処理装置と、この演算処理装置の出力を画面に表示する表示装置と、上記演算処理装置に入力を行う入力手段とを備え、
上記演算処理装置は、
上記表示装置の画面に、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命の入力を促す表示を行わせる促し画面出力手段と、
実行命令に応答して上記有為差有無の判定および少なくとも断定できる寿命差の演算を行いその演算結果を上記表示装置の画面に出力する有為差判定・有為寿命差演算手段とを備え、
上記有為差判定・有為寿命差演算手段は、
上記入力情報のワイブルスロープの値および仮設定寿命から定まるワイブル分布に従った乱数を水準1のロットの試験個数分発生させる第1乱数発生手段と、
第1乱数発生手段と同じワイブル分布に従った乱数を水準2のロットの試験個数分発生させる第2乱数発生手段と、
第1乱数発生手段で発生させた設定個数の乱数から所定の寿命計算方法で求まる寿命と第2乱数発生手段で発生させた設定個数の乱数から上記所定の寿命計算方法で求まる寿命とから、1組の水準1,水準2の寿命の寿命比を算出する寿命比算出手段と、
これら第1乱数発生手段,第2乱数発生手段,および寿命比算出手段の処理を設定回数繰り返し、上記設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成して累積確率分布から、設定信頼幅内の最大の寿命比を読み取る累積確率分布分析手段と、
上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手段で得た最大寿命比を超える場合に有為差有りと判定し、最大寿命比以内の場合に有為差無しと判定する有為差有無判定手段と、
この手段で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手段の処理を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す倍率変更繰り返し手段と、
この倍率変更繰り返し手段で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する寿命倍率関係グラフ作成手段と、
上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする有為寿命差読み取り手段と、
上記有為差有無判定手段で判定した結果および上記有為寿命差読み取り手段で読み取った、少なくとも断定できる寿命差を表示装置に表示させる結果出力手段と、
を含む、加速試験における有為差有無判定・有為寿命差見積もり装置。
【請求項3】
コンピュータで実行可能なプログラムであって、
表示装置の画面に、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命の入力を促す表示を行わせる促し画面出力手順と、
実行命令に応答して上記有為差有無の判定および少なくとも断定できる寿命差の演算を行いその演算結果を上記表示装置の画面に出力する有為差判定・有為寿命差演算手順とを含み、
この有為差判定・有為寿命差演算手順は、
上記入力情報のワイブルスロープの値および仮設定寿命から定まるワイブル分布に従った乱数を水準1のロットの試験個数分発生させる第1乱数発生手順と、
第1乱数発生手順と同じワイブル分布に従った乱数を水準2のロットの試験個数分発生させる第2乱数発生手順と、
第1乱数発生手順で発生させた設定個数の乱数から所定の寿命計算方法で求まる寿命と第2乱数発生手順で発生させた設定個数の乱数から上記所定の寿命計算方法で求まる寿命とから、1組の水準1,水準2の寿命の寿命比を算出する寿命比算出手順と、
これら第1乱数発生手順,第2乱数発生手順,および寿命比算出手順を設定回数繰り返し、上記設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成して累積確率分布から、設定信頼幅内の最大の寿命比を読み取る累積確率分布分析手順と、 上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手順で得た最大寿命比を超える場合に有為差有りと判定し、最大寿命比以内の場合に有為差無しと判定する有為差有無判定手順と、
この手順で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手順を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す倍率変更繰り返し手順と、
この倍率変更繰り返し手順で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する寿命倍率関係グラフ作成手順と、
上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする有為寿命差読み取り手順と、
上記有為差有無判定手順で判定した結果および上記有為寿命差読み取り手順で読み取った、少なくとも断定できる寿命差を表示装置に表示させる結果出力手順と、
を含む、加速試験における有為差有無判定・有為寿命差見積もりプログラム。
【請求項1】
軸受等の機械部品または試験片からなる試験対象品を使用環境よりも厳しい所定の環境条件におき、破損が発生した時間から寿命を算出する加速試験において、試験対象品の2つのロットの寿命から、有為差有無の判定、および少なくとも断定できる、倍率による寿命差を算出する方法であって、
コンピュータに対し、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命を入力する過程と、
上記コンピュータに、上記有為差有無の判定、および少なくとも断定できる寿命差を演算させ演算結果を表示装置の画面に表示させるコンピュータ演算処理過程とを含み、
上記コンピュータ演算処理過程として、
試験対象品の上記入力情報のワイブルスロープの値および仮設定寿命から定まるワイブル分布に従った乱数を水準1のロットの試験個数分発生させる第1乱数発生手順と、
第1乱数発生手順と同じワイブル分布に従った乱数を水準2のロットの試験個数分発生させる第2乱数発生手順と、
第1乱数発生手順で発生させた設定個数の乱数から所定の寿命計算方法で求まる寿命と第2乱数発生手順で発生させた設定個数の乱数から上記所定の寿命計算方法で求まる寿命とから、1組の水準1,水準2の寿命の寿命比を算出する寿命比算出手順と、
これら第1乱数発生手順,第2乱数発生手順,および寿命比算出手順を設定回数繰り返し、上記設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成して累積確率分布から、設定信頼幅内の最大の寿命比を読み取る累積確率分布分析手順と、 上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手順で得た最大寿命比を超える場合に有為差有りと判定し、最大寿命比以内の場合に有為差無しと判定する有為差有無判定手順と、
この手順で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手順を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す倍率変更繰り返し手順と、
この倍率変更繰り返し手順で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する寿命倍率関係グラフ作成手順と、
上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする有為寿命差読み取り手順と、
上記有為差有無判定手順で判定した結果および上記有為寿命差読み取り手順で読み取った、少なくとも断定できる寿命差を表示装置に表示させる結果出力手順と、
を含む、加速試験における有為差有無判定・有為寿命差見積もり方法。
【請求項2】
軸受等の機械部品または試験片からなる試験対象品を使用環境よりも厳しい所定の環境条件におき、破損が発生した時間から寿命を算出する加速試験において、試験対象品の2つのロットの寿命から、有為差有無の判定、および少なくとも断定できる、倍率による寿命差を算出する装置であって、
演算処理装置と、この演算処理装置の出力を画面に表示する表示装置と、上記演算処理装置に入力を行う入力手段とを備え、
上記演算処理装置は、
上記表示装置の画面に、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命の入力を促す表示を行わせる促し画面出力手段と、
実行命令に応答して上記有為差有無の判定および少なくとも断定できる寿命差の演算を行いその演算結果を上記表示装置の画面に出力する有為差判定・有為寿命差演算手段とを備え、
上記有為差判定・有為寿命差演算手段は、
上記入力情報のワイブルスロープの値および仮設定寿命から定まるワイブル分布に従った乱数を水準1のロットの試験個数分発生させる第1乱数発生手段と、
第1乱数発生手段と同じワイブル分布に従った乱数を水準2のロットの試験個数分発生させる第2乱数発生手段と、
第1乱数発生手段で発生させた設定個数の乱数から所定の寿命計算方法で求まる寿命と第2乱数発生手段で発生させた設定個数の乱数から上記所定の寿命計算方法で求まる寿命とから、1組の水準1,水準2の寿命の寿命比を算出する寿命比算出手段と、
これら第1乱数発生手段,第2乱数発生手段,および寿命比算出手段の処理を設定回数繰り返し、上記設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成して累積確率分布から、設定信頼幅内の最大の寿命比を読み取る累積確率分布分析手段と、
上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手段で得た最大寿命比を超える場合に有為差有りと判定し、最大寿命比以内の場合に有為差無しと判定する有為差有無判定手段と、
この手段で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手段の処理を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す倍率変更繰り返し手段と、
この倍率変更繰り返し手段で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する寿命倍率関係グラフ作成手段と、
上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする有為寿命差読み取り手段と、
上記有為差有無判定手段で判定した結果および上記有為寿命差読み取り手段で読み取った、少なくとも断定できる寿命差を表示装置に表示させる結果出力手段と、
を含む、加速試験における有為差有無判定・有為寿命差見積もり装置。
【請求項3】
コンピュータで実行可能なプログラムであって、
表示装置の画面に、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命の入力を促す表示を行わせる促し画面出力手順と、
実行命令に応答して上記有為差有無の判定および少なくとも断定できる寿命差の演算を行いその演算結果を上記表示装置の画面に出力する有為差判定・有為寿命差演算手順とを含み、
この有為差判定・有為寿命差演算手順は、
上記入力情報のワイブルスロープの値および仮設定寿命から定まるワイブル分布に従った乱数を水準1のロットの試験個数分発生させる第1乱数発生手順と、
第1乱数発生手順と同じワイブル分布に従った乱数を水準2のロットの試験個数分発生させる第2乱数発生手順と、
第1乱数発生手順で発生させた設定個数の乱数から所定の寿命計算方法で求まる寿命と第2乱数発生手順で発生させた設定個数の乱数から上記所定の寿命計算方法で求まる寿命とから、1組の水準1,水準2の寿命の寿命比を算出する寿命比算出手順と、
これら第1乱数発生手順,第2乱数発生手順,および寿命比算出手順を設定回数繰り返し、上記設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成して累積確率分布から、設定信頼幅内の最大の寿命比を読み取る累積確率分布分析手順と、 上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手順で得た最大寿命比を超える場合に有為差有りと判定し、最大寿命比以内の場合に有為差無しと判定する有為差有無判定手順と、
この手順で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手順を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す倍率変更繰り返し手順と、
この倍率変更繰り返し手順で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する寿命倍率関係グラフ作成手順と、
上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする有為寿命差読み取り手順と、
上記有為差有無判定手順で判定した結果および上記有為寿命差読み取り手順で読み取った、少なくとも断定できる寿命差を表示装置に表示させる結果出力手順と、
を含む、加速試験における有為差有無判定・有為寿命差見積もりプログラム。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【公開番号】特開2008−128698(P2008−128698A)
【公開日】平成20年6月5日(2008.6.5)
【国際特許分類】
【出願番号】特願2006−311138(P2006−311138)
【出願日】平成18年11月17日(2006.11.17)
【出願人】(000102692)NTN株式会社 (9,006)
【Fターム(参考)】
【公開日】平成20年6月5日(2008.6.5)
【国際特許分類】
【出願日】平成18年11月17日(2006.11.17)
【出願人】(000102692)NTN株式会社 (9,006)
【Fターム(参考)】
[ Back to top ]