JP2008128698A - Method and apparatus for estimating necessary life difference of significant difference determination in accelerated test - Google Patents

Method and apparatus for estimating necessary life difference of significant difference determination in accelerated test Download PDF

Info

Publication number
JP2008128698A
JP2008128698A JP2006311138A JP2006311138A JP2008128698A JP 2008128698 A JP2008128698 A JP 2008128698A JP 2006311138 A JP2006311138 A JP 2006311138A JP 2006311138 A JP2006311138 A JP 2006311138A JP 2008128698 A JP2008128698 A JP 2008128698A
Authority
JP
Japan
Prior art keywords
life
difference
procedure
magnification
level
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2006311138A
Other languages
Japanese (ja)
Inventor
Takumi Fujita
工 藤田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NTN Corp
Original Assignee
NTN Corp
NTN Toyo Bearing Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by NTN Corp, NTN Toyo Bearing Co Ltd filed Critical NTN Corp
Priority to JP2006311138A priority Critical patent/JP2008128698A/en
Publication of JP2008128698A publication Critical patent/JP2008128698A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To provide a method for determining the presence of significant difference from the life between two lots obtained in an accelerated test and calculating life difference due to magnification. <P>SOLUTION: The random numbers corresponding to the test number of a lot of a level 1 are produced from a certain Weibull distribution, the random number corresponding to the test number of a lot of a level 2 are produced from the same Weibull distribution and the life calculated from the produced random numbers is operated to calculate one set of a life ratio. This processing is repeated the set number of times to calculate the life ration of the number of sets of the set number of times to form the probability distribution and cumulated probability distribution of this life ratio. The maximum and minimum life ratios in set reliable width are read to determine the presence of significant difference. When the presence of the significant difference is determined, the Weibull distribution is changed to a different Weibull distribution to repeat the analyzing step of the cumulated probability distribution. A graph showing the relationship between the set magnification obtained by the repetition of the analyzing step and the maximum life in the reliable width is created and the value of the set magnification corresponding to the life ratio of the lives of the levels 1 and 2 in input data is read to be set as the life difference capable of determining at least the read value. <P>COPYRIGHT: (C)2008,JPO&INPIT

Description

この発明は、軸受等の機械部品または試験片からなる試験対象品を使用環境よりも厳しい所定の環境条件におき、破損が発生した時間から寿命を算出する加速試験において、試験対象品の2つのロット間で有為差が有るか否かの判定、およびある場合にどの程度の有為差があるかを見積もる有為差有無判定・有為寿命差見積もり方法,装置、プログラムに関する。   In the accelerated test for calculating the life from the time when the failure occurs, the test object consisting of a mechanical part such as a bearing or a test object is placed in a predetermined environmental condition that is severer than the use environment. The present invention relates to a method, apparatus, and program for determining whether there is a significant difference between lots, and for determining whether there is a significant difference in some cases, and for determining whether there is a significant difference or estimating a significant life difference.

寿命試験は、軸受等の機械部品の性能を評価するために欠かせない試験の1つである。寿命試験には、大きく分けて(1) 実機の使用環境に近い条件で試験を行う実機試験と、(2) 比較的過酷な条件で寿命試験を行う加速試験がある。前者は、製品が有限時間内に破損するケースが極めて少ないため、ある目標時間まで破損することなく試験が継続すれば、寿命は問題ないと判断する試験である(以下、このような試験を「打切り試験」と呼ぶ)。一方、後者は、比較的短時間で破損が発生するので、ワイブルプロットで寿命が算出でき(例えば非特許文献1)、その算出寿命から性能の優劣を判定する試験である(以下、このような試験を「加速試験」と呼ぶ)。   The life test is one of the tests indispensable for evaluating the performance of mechanical parts such as bearings. Life tests can be broadly divided into (1) actual machine tests that test under conditions close to the actual operating environment and (2) accelerated tests that perform life tests under relatively severe conditions. The former is a test that determines that there is no problem in the service life if the test continues without damage until a certain target time because the product is rarely damaged within a finite time (hereinafter referred to as " Called “censored test”). On the other hand, since the latter occurs in a relatively short time, the life can be calculated by Weibull plot (for example, Non-Patent Document 1), and is a test for determining the superiority or inferiority of the performance from the calculated life (hereinafter, such as this The test is called “accelerated test”).

従来より、寿命試験は経験を積んだ熟練者が行っており、試験条件や試験個数を決める寿命試験の設計と寿命試験結果の解釈に対して経験的に確からしい判断ができたと考えられる。
図14に、従来から行われてきた寿命試験の設計と寿命試験結果の解釈の手順を、打切り試験と加速試験ごとに示す。
また、現在、寿命試験において経験的に判断されているものの詳細を、表1に示す。
Conventionally, the experienced life test has been conducted by experienced experts, and it is considered that he was able to make empirical judgments on the design of the life test and the interpretation of the life test results to determine the test conditions and the number of tests.
FIG. 14 shows a procedure for interpreting a life test design and a life test result, which has been conventionally performed, for each of an abort test and an acceleration test.
The details of what is currently determined empirically in the life test are shown in Table 1.

Figure 2008128698
Figure 2008128698

なお、ワイブル分布を機械部品の寿命判断に用いるものは、種々の特許文献,非特許文献に提案されている。
特開2006−040203号公報 特開2002−277382号公報 特開2005−226829号公報 真壁肇著、信頼性工学入門79、1991年発行
In addition, what uses Weibull distribution for the lifetime judgment of a machine component is proposed by various patent documents and nonpatent literature.
JP 2006-040203 A JP 2002-277382 A JP 2005-226829 A Author Makabe Satoshi, Introduction to Reliability Engineering 79, published in 1991

加速試験は、上記のように破損時間から寿命を算出し、その算出寿命から性能の優劣を判定する試験であるが、2ロットの寿命試験結果から求めた算出寿命の間に、本当に有為差があるかどうかを確かめたい状況が多くある。
従来より、算出寿命の有為差判定は信頼幅という概念を使って行われてきた。しかし、この判定方法には、次のようにいくつかの問題がある。
In the accelerated test, the life is calculated from the failure time as described above, and the superiority or inferiority of the performance is determined from the calculated life. However, there is a significant difference between the calculated life obtained from the life test results of 2 lots. There are many situations where you want to see if there is any.
Conventionally, the significant difference determination of the calculated life has been performed using the concept of confidence width. However, this determination method has several problems as follows.

まず、信頼幅は、2水準間での寿命の有意差を判断できないという問題がある。その理由は、どの程度信頼幅が重なっていたら寿命に優位差があるのかということを定量的に判断できないためである。信頼幅を使って2水準間での寿命の有意差を判定する場合、次のような手順が必要になる。以下、その手順を図15を使って説明する。
1) L10寿命におけるメディアンランクの信頼水準の分布(図中Aの分布)を求める。ここで、この分布はF値表とファイサー(Fisher)の近似式を用いて求める。
2) L10寿命におけるメディアンランクの信頼水準の分布とワイブルスロープからL10のばらつき頻度の分布(図中Bの分布)を求める。
3) このばらつき分布f(x)を使って、図中の式の積分を取り、寿命倍率nが同じ寿命分布から発生する頻度を調べる。これは寿命比のばらつき分布である。
4) 3)の分布の累積確率分布を作成し、そのグラフから90%以上の確率で有意差有と判定できる寿命倍率を計算する。
First, the reliability range has a problem that a significant difference in lifetime between two levels cannot be determined. The reason is that it is impossible to quantitatively determine how much the reliability range overlaps and there is a difference in life. In order to determine a significant difference in life between two levels using the confidence width, the following procedure is required. The procedure will be described below with reference to FIG.
1) Obtain the median rank confidence level distribution (distribution A in the figure) in the L10 life. Here, this distribution is obtained using an F-number table and an approximate expression of Fisher.
2) The distribution of the distribution frequency of L10 (distribution B in the figure) is obtained from the distribution of the confidence level of the median rank in the L10 life and the Weibull slope.
3) Using this variation distribution f (x), the integral of the equation in the figure is taken, and the frequency at which the life multiplication factor n occurs from the same life distribution is examined. This is a variation distribution of the life ratio.
4) The cumulative probability distribution of the distribution of 3) is created, and the life multiplication factor that can be determined as having a significant difference with a probability of 90% or more is calculated from the graph.

上のような手順で、従来の信頼幅を使って寿命の有意差検定や設計が可能になる。しかし、このような手順は複雑であり、実用化が難しい。このため、信頼幅を使った寿命の有意差検定が今まで行われてこなかったのは当然であるといえる。   The above procedure makes it possible to test and design a significant difference in lifetime using the conventional confidence range. However, such a procedure is complicated and difficult to put into practical use. For this reason, it can be said that the test of significant difference in life using the confidence range has not been performed so far.

レオナード・ジー・ジョンソン(LEONARD G. JOHNSON)は、ワイブルスロープごと、試験個数ごとに、2水準間の平均寿命とL10 寿命の有意差を判定する方法を提案している)。しかし、その方法も、実際には使われていないのが現状である。これは、手軽な形で優位差検定を行うツールを残していないことに原因がある。   LEONARD G. JOHNSON proposes a method to determine the significant difference between the average life and L10 life between two levels for each Weibull slope and number of specimens). However, the method is not actually used at present. This is due to the fact that there is no tool to perform the dominant difference test in an easy way.

また、従来は、有為差有りと判定を行っても、少なくともどの程度の寿命差があると言えるのかを、定量的に示す適切な方法がなかった。   Conventionally, there has been no appropriate method for quantitatively indicating at least how much life difference can be said even if it is determined that there is a significant difference.

この発明の目的は、2つのロット間での有為差有無の判定を定量的に適切に行え、かつ有為差があると判定される場合に、少なくともどの程度の寿命差があると言えるのかを、簡単かつ迅速に試算することができ、かつ定量的に求められて信頼性の高いものとでき、熟練者でなくても必要な寿命差を適切に試算することのできる方法、装置、およびコンピュータプログラムを提供することである。   The purpose of this invention is to determine the existence of a significant difference between two lots quantitatively and appropriately, and when it is determined that there is a significant difference, at least how much life difference can be said to be Can be easily and quickly calculated, can be quantitatively obtained and reliable, and can be used to appropriately calculate the required life difference even if not skilled. To provide a computer program.

この発明の加速試験における有為差有無判定・有為寿命差見積もり方法は、軸受等の機械部品または試験片からなる試験対象品を使用環境よりも厳しい所定の環境条件におき、破損が発生した時間から寿命を算出する加速試験において、試験対象品の2つのロットの寿命から、有為差有無の判定、および少なくとも断定できる、倍率による寿命差を算出する方法であって、
コンピュータに対し、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命を入力する過程(M1)と、
上記コンピュータに、上記有為差有無の判定、および少なくとも断定できる寿命差を演算させ演算結果を表示装置の画面に表示させるコンピュータ演算処理過程(M2)とを含む。
上記寿命は、例えばL10寿命(90%の信頼度の寿命)や、L50寿命(50%の信頼度の寿命)等の所定信頼度の寿命ある。
In the accelerated test of the present invention, the presence / absence of significant difference / estimated life difference estimation method is determined by placing a test object consisting of a mechanical part such as a bearing or a test piece under predetermined environmental conditions that are severer than the operating environment, resulting in damage. In the accelerated test for calculating the life from time, a method for calculating the difference in life by magnification, which can be determined and at least determined from the life of two lots of the test object,
A process (M1) of inputting, as input information, the value of the Weibull slope of the Weibull distribution of the test object, the number of each test of the two lots of Level 1 and Level 2 to be compared, and the life of the test result as input information;
And a computer calculation process (M2) for causing the computer to determine whether or not there is a significant difference and to calculate at least a life difference that can be determined and to display a calculation result on a screen of a display device.
The lifetime is a lifetime with a predetermined reliability such as an L10 lifetime (90% reliability lifetime) or an L50 lifetime (50% reliability lifetime).

上記コンピュータ演算処理過程(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)とを含む。
As the computer calculation process (M2),
A first random number generation procedure (N21) for generating a random number according to the Weibull distribution determined from the value of the Weibull slope of the input information of the test target product and the temporarily set life for the number of test pieces of the level 1 lot;
A second random number generation procedure (N22) for generating random numbers according to the same Weibull distribution as the first random number generation procedure (N21) for the test number of lots of level 2.
From the set number of random numbers generated in the first random number generation procedure (N21) by the predetermined lifetime calculation method and the set number of random numbers generated in the second random number generation procedure (N22), the predetermined lifetime calculation method. A life ratio calculation procedure (N23) for calculating a life ratio of a set of level 1 and level 2 lifespans from the obtained lifespan,
The first random number generation procedure (N21), the second random number generation procedure (N23), and the life ratio calculation procedure (N24) are repeated a set number of times to obtain the life ratio of the set number of sets, and the probability distribution of the life ratio A cumulative probability distribution analysis procedure (N24) for creating a cumulative probability distribution and reading the maximum life ratio within the set confidence range from the cumulative probability distribution;
It is judged that there is a significant difference when the life ratio of the life of Level 1 and Level 2 in the above input information exceeds the maximum life ratio obtained by the cumulative probability distribution analysis procedure (N24), and is valid when it is within the maximum life ratio. And a significant difference presence / absence determination procedure (N25) for determining that there is no difference.

また、上記コンピュータ演算処理過程(N2)として、さらに、
上記有為差有無判定手順(N25)で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手順を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す倍率変更繰り返し手順(N26)と、
この倍率変更繰り返し手順(N26)で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する寿命倍率関係グラフ作成手順(N27)と、
上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする有為寿命差読み取り手順(N28)と、
上記有為差有無判定手順(N25)で判定した結果および上記有為寿命差読み取り手順(N28)で読み取った、少なくとも断定できる寿命差を表示装置に表示させる結果出力手順(N29)と、を含む。
In addition, as the computer calculation process (N2),
The cumulative probability distribution analysis is executed when it is determined that there is a significant difference in the significant difference presence / absence determination procedure (N25), and the Weibull distribution is sequentially changed to a Weibull distribution having a setting magnification with a different provisional set life for each repetition. A magnification change repetition procedure (N26) for repeating the procedure until the set magnification reaches the set maximum magnification from the set minimum magnification,
A life magnification relationship graph creation procedure (N27) for creating a life magnification relationship graph, which is a graph showing the relationship between the set magnification repeatedly changed and the maximum life ratio within the reliability range, obtained in this magnification change repetition procedure (N26). )When,
Reading the value of the set magnification corresponding to the life ratio of the life of Level 1 and Level 2 in the input information from the life magnification relationship graph and reading the significant life difference as a life difference that can be determined at least. (N28),
And a result output procedure (N29) for displaying on the display device at least a determinable life difference read in the significant difference presence / absence determination procedure (N25) and the significant life difference reading procedure (N28). .

上記第1,第2乱数発生手順(N21,N22)において、乱数発生に用いるワイブル分布を特定するための寿命は、適宜想定した寿命または実績値を設定すれば良い。
また、寿命比算出手順(N23)において、乱数発生手順で発生させた試験個数の乱数から寿命を求める所定の寿命計算方法は、加速試験において従来から用いられている中の適宜の寿命計算方法を用いれば良い。加速試験においては、破損の発生した時間からワイブルプロットで寿命を求めることが行われており、このような破損時間から寿命を求める計算方法を、上記寿命比算出手順における寿命計算方法として用いる。
ワイブル分布は、次式、
In the first and second random number generation procedures (N21, N22), as the life for specifying the Weibull distribution used for random number generation, a life or actual value estimated as appropriate may be set.
In addition, in the life ratio calculation procedure (N23), the predetermined life calculation method for obtaining the life from the random number of the test pieces generated in the random number generation procedure is an appropriate life calculation method conventionally used in the accelerated test. Use it. In the accelerated test, the life is obtained by a Weibull plot from the time when the breakage occurs, and a calculation method for obtaining the life from the breakage time is used as the life calculation method in the life ratio calculation procedure.
The Weibull distribution is given by

Figure 2008128698
Figure 2008128698

ただし、m:ワイブルスロープ、α:尺度因子、γ:最小寿命、
によって特定される。
Where m: Weibull slope, α: scale factor, γ: minimum life,
Specified by.

軸受等の機械部品の寿命は、ワイブル分布に従うとされている。ワイブル分布は、ワイブルスロープm、尺度因子α、最小寿命γの3つのパラメータを持っており、ワイブルスロープmによって指数分布、対数正規分布、正規分布を表現できる万能分布として知られている。量産される軸受等では、ワイブルスロープは実績値が既知である場合が多く、この発明方法において、ワイブルスロープには、試験対象品の実績値を用いることが好ましい。実績値がない場合は、適宜の方法で見積もったワイブルスロープを用いてもよい。最小寿命γは、種々の規格、例えばISO等によって計算方法が定められており、そのように定められたいずれかの計算方法を用いることが好ましい。尺度因子αは、ワイブルスロープの値、要求寿命の信頼度、要求寿命の値、および上記最小寿命γから一義的に決定される演算式があり、その演算式を用いて特定しても良い。   The life of mechanical parts such as bearings is said to follow the Weibull distribution. The Weibull distribution has three parameters, a Weibull slope m, a scale factor α, and a minimum life γ, and is known as a universal distribution that can express an exponential distribution, a lognormal distribution, and a normal distribution by the Weibull slope m. In mass-produced bearings and the like, the actual value of the Weibull slope is often known, and in the method of the present invention, it is preferable to use the actual value of the product to be tested as the Weibull slope. If there is no actual value, a Weibull slope estimated by an appropriate method may be used. The calculation method of the minimum life γ is determined by various standards such as ISO, and it is preferable to use any one of the calculation methods determined as such. The scale factor α has an arithmetic expression that is uniquely determined from the value of the Weibull slope, the reliability of the required life, the value of the required life, and the minimum life γ, and may be specified using the arithmetic expression.

この方法において、累積分布作成手順(N24)で行う処理は、同じ寿命分布を持つロットから験個数分の試験対象品を2組抜きとり、寿命試験を行った結果、どの程度の寿命差が出るのかを設定回数確認することに対応している。同じ寿命分布を持つロットから試験対象品を抜きとった試験においても、寿命はばらついてる。その寿命比のばらつきを累積確率分布の適宜の信頼区間、例えば5%と95%の区間(90%信頼区間)として定義すると、その寿命ばらつきは、上記累積確率分布から算出できる。この結果は、同じ寿命分布を持つロットから試験個数の試験対象品を2組抜きとって寿命比較を行っても、90%のものは、累積確率分布に示された信頼区間の範囲で寿命比がばらついてしまうことを示している。 逆にいうと上記信頼区間の最大の寿命差があれば、そのロット間の寿命分布が異なっている可能性が高いということになる。
したがって、上記有為差有無判定手順(N25)として、上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手順で得た最大寿命比を超える場合に有為差有りと判定し、最大寿命比以内の場合に有為差無しと判定することで、定量的に適切に有為差有無を判定することができる。
In this method, in the cumulative distribution creation procedure (N24), two lifespans are obtained as a result of performing a life test after extracting two sets of test target products from a lot having the same life distribution. It corresponds to checking the number of times set. Even in a test in which a test object is extracted from a lot having the same life distribution, the life varies. If the variation in the life ratio is defined as an appropriate confidence interval of the cumulative probability distribution, for example, a 5% and 95% interval (90% confidence interval), the life variation can be calculated from the cumulative probability distribution. This result shows that even if two sets of test target products are extracted from lots with the same life distribution, the life ratio is 90% within the confidence interval shown in the cumulative probability distribution. It shows that it will vary. In other words, if there is a maximum life difference in the confidence interval, there is a high possibility that the life distribution between the lots is different.
Therefore, as the significant difference presence / absence determination procedure (N25), when the life ratio of the level 1 and level 2 lives in the input information exceeds the maximum life ratio obtained by the cumulative probability distribution analysis procedure, there is a significant difference. By determining and determining that there is no significant difference within the maximum life ratio, the presence or absence of the significant difference can be appropriately determined quantitatively.

また、寿命が異なるワイブル分布との組み合わせで、上記同様の累積確率分布の計算を行い、例えば、ワイブル分布と寿命が1.1 、1.2 …50倍異なるワイブル分布の組み合わせで計算を行うと、その組み合わせが1.1 、1.2 …50倍と変化するごとに、寿命差が1.1 、1.2 …50倍を中心にしてばらつくはずである。上記倍率変更繰り返し手順(N26)は、このような累積確率分布の計算、およびその確率分布についての上記確率分布分析手順(N24)を行う。
この倍率変更繰り返し手順(N26)で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成し、このグラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取ると、その読み取った値が少なくとも断定できる寿命差となる。
Also, the cumulative probability distribution similar to the above is calculated in combination with the Weibull distribution with different lifetimes. For example, when the calculation is performed with the combination of the Weibull distribution and the Weibull distribution with different lifetimes of 1.1, 1.2, ... Every time it changes 1.1, 1.2 ... 50 times, the life difference should vary around 1.1, 1.2 ... 50 times. The magnification change repetition procedure (N26) performs the calculation of the cumulative probability distribution and the probability distribution analysis procedure (N24) for the probability distribution.
A life magnification relationship graph, which is a graph showing the relationship between the set magnification repeatedly changed and the maximum life ratio within the reliability range, obtained in the magnification change repetition procedure (N26), is created, and from the graph, the input information When the value of the set magnification corresponding to the life ratio of the life of level 1 and level 2 is read, the read value is at least a life difference that can be determined.

この発明方法は、このようにワイブル分布に基づく乱数シミュレーションにより有為差有無を判定し、また少なくとも断定できる寿命差を求めるため、これら有為差有無および少なくとも断定できる寿命差を、簡単かつ迅速に試算することができる。また、定量的に求められて信頼性の高いものとでき、熟練者でなくても試算することができる。この乱数シミュレーションに用いるプログラムは、条件入力を行えば勝手に有為差有無を判定および少なくとも断定できる寿命差の検定を行ってくれる手軽なツールであることが優れた点である。   In the method of the present invention, the presence or absence of a significant difference is determined by random simulation based on the Weibull distribution, and at least the life difference that can be determined is obtained. You can estimate. Further, it can be quantitatively obtained and can be highly reliable, and can be calculated even if it is not an expert. The program used for this random number simulation is excellent in that it is a handy tool that can test the existence of a significant difference and at least determine the life difference that can be determined by inputting conditions.

この発明の加速試験における有為差有無判定・有為寿命差見積もり装置は、軸受等の機械部品または試験片からなる試験対象品を使用環境よりも厳しい所定の環境条件におき、破損が発生した時間から寿命を算出する加速試験において、試験対象品の2つのロットの寿命から、有為差有無の判定、および少なくとも断定できる、倍率による寿命差を算出する装置であって、
演算処理装置(1)と、この演算処理装置(1)の出力を画面に表示する表示装置(2)と、上記演算処理装置(1)に入力を行う入力手段(3)とを備える。
In the accelerated test according to the present invention, the existence / non-existence of a significant difference / estimated life difference estimation device is caused by placing a test object consisting of a mechanical part such as a bearing or a test piece under a predetermined environmental condition that is severer than the use environment. In the accelerated test for calculating the life from time, it is a device for calculating the difference in life by magnification, which can be determined and at least determined from the life of two lots of the test object,
An arithmetic processing device (1), a display device (2) for displaying the output of the arithmetic processing device (1) on a screen, and an input means (3) for inputting to the arithmetic processing device (1) are provided.

上記演算処理装置(1)は、上記表示装置(2)の画面に、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命の入力を促す表示を行わせる促し画面出力手段(7M)と、実行命令に応答して上記有為差有無の判定および少なくとも断定できる寿命差の演算を行いその演算結果を上記表示装置(2)の画面に出力する有為差判定・有為寿命差演算手段(52)とを備える。   The arithmetic processing unit (1) displays, on the screen of the display unit (2), as input information, the value of the Weibull distribution of the test target product, each test of the two lots of level 1 and level 2 to be compared. Prompt screen output means (7M) that prompts the input of the number and test result life, and in response to the execution command, the presence / absence of the significant difference is calculated and at least the life difference that can be determined is calculated, and the calculation result is obtained. It is provided with a significant difference determination / significant life difference calculating means (52) for outputting to the screen of the display device (2).

上記有為差判定・有為寿命差演算手段(52)は、第1乱数発生手段(53)、第2乱数発生手段(54)、寿命比算出手段(55)、累積確率分布分析手段(56)、有為差有無判定手段(57)、倍率変更繰り返し手段(58)、寿命倍率関係グラフ作成手段(59)、有為寿命差読み取り手段(60)、および結果出力手段(61)を有する。   The significant difference determination / significant life difference calculating means (52) includes a first random number generating means (53), a second random number generating means (54), a life ratio calculating means (55), and a cumulative probability distribution analyzing means (56). ), Significant difference presence / absence determination means (57), magnification change repetition means (58), life magnification relation graph creation means (59), significant life difference reading means (60), and result output means (61).

第1乱数発生手段(53)は、上記入力情報のワイブルスロープの値および仮設定寿命から定まるワイブル分布に従った乱数を水準1のロットの試験個数分発生させる手段である。
第2乱数発生手段(54)は、第1乱数発生手段(53)と同じワイブル分布に従った乱数を水準2のロットの試験個数分発生させる手段である。
寿命比算出手段(55)は、第1乱数発生手段(53)で発生させた設定個数の乱数から所定の寿命計算方法で求まる寿命と、第2乱数発生手段(54)で発生させた設定個数の乱数から上記所定の寿命計算方法で求まる寿命とから、1組の水準1,水準2の寿命の寿命比を算出する手段である。
The first random number generation means (53) is a means for generating random numbers according to the Weibull distribution determined from the value of the Weibull slope of the input information and the temporarily set lifetime for the number of test pieces of the level 1 lot.
The second random number generation means (54) is means for generating random numbers according to the same Weibull distribution as the first random number generation means (53) for the number of test pieces of the level 2 lot.
The life ratio calculation means (55) has a life determined by a predetermined life calculation method from a set number of random numbers generated by the first random number generation means (53) and a set number generated by the second random number generation means (54). Is a means for calculating the life ratio of a set of level 1 and level 2 lifespans from the lifespan obtained by the predetermined lifespan calculation method from the random number.

累積確率分布分析手段(56)は、これら第1乱数発生手段(53),第2乱数発生手段(54),および寿命比算出手段(55)の処理を設定回数繰り返し、上記設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成して累積確率分布から、設定信頼幅内の最大の寿命比を読み取る手段である。   The cumulative probability distribution analysis means (56) repeats the processes of the first random number generation means (53), the second random number generation means (54), and the life ratio calculation means (55) a set number of times, and the number of sets of the set number of times The life ratio is obtained, a probability distribution of the life ratio and a cumulative probability distribution are created, and the maximum life ratio within the set reliability range is read from the cumulative probability distribution.

有為差有無判定手段(57)は、上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手段(56)で得た最大寿命比を超える場合に有為差有りと判定し、最大寿命比以内の場合に有為差無しと判定する手段である。   The significant difference presence / absence determining means (57) determines that there is a significant difference when the life ratio of level 1 and level 2 life in the input information exceeds the maximum life ratio obtained by the cumulative probability distribution analysis means (56). It is a means for determining that there is no significant difference when it is within the maximum life ratio.

倍率変更繰り返し手段(58)は、有為差有無判定手段(57)で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手段(56)の処理を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す手段である。
寿命倍率関係グラフ作成手段(59)は、倍率変更繰り返し手段(58)で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する手段である。
有為寿命差読み取り手段(60)は、上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする手段である。
The magnification change repetition means (58) is executed when there is a significant difference in the significant difference presence / absence determination means (57), and the Weibull distribution is sequentially changed to a Weibull distribution having a set magnification with a temporarily set lifetime different for each repetition. In other words, the cumulative probability distribution analysis means (56) repeats the processing until the set magnification becomes the set maximum magnification from the set minimum magnification.
The life magnification relationship graph creating means (59) is a life magnification relationship graph, which is a graph showing the relationship between the setting magnification repeatedly changed and the maximum life ratio within the reliability range, obtained by the magnification change repetition means (58). It is a means to create.
The significant life difference reading means (60) reads the value of the set magnification corresponding to the life ratio of the life of level 1 and level 2 in the input information from the life magnification relationship relation graph, and at least determines the read value. It is a means to make the life difference possible.

結果出力手段(61)は、上記有為差有無判定手段(57)で判定した結果、および上記有為寿命差読み取り手段(60)で読み取った、少なくとも断定できる寿命差を表示装置(2)に表示させる手段である。   The result output means (61) displays the result determined by the significant difference presence / absence determining means (57) and the life difference that can be determined at least determined by the significant life difference reading means (60) on the display device (2). It is a means to display.

この構成の加速試験における有為差有無判定・有為寿命差見積もり装置によると、この発明の有為差有無判定・有為寿命差見積もり方法に使用されて、有為差の有無の判定、および有為差があると判定される場合に、少なくともどの程度の寿命差があると言えるのかを、簡単かつ迅速に試算することができ、かつ定量的に求められて信頼性の高いものとでき、熟練者でなくても少なくとも断定できる寿命差を適切に試算することができる。   According to the significant difference presence / absence determination / significant life difference estimation apparatus in the accelerated test of this configuration, the presence / absence of significant difference is used in the significant difference existence determination / significant life difference estimation method of the present invention, and When it is determined that there is a significant difference, at least how much life difference can be said can be easily and quickly calculated, and quantitatively obtained and reliable. Even if it is not an expert, the life difference which can be determined at least can be estimated appropriately.

この発明の加速試験における有為差有無判定・有為寿命差見積もりプログラム(51)は、コンピュータで実行可能なプログラムであって、
上記表示装置(2)の画面に、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命の入力を促す表示を行わせる促し画面出力手順(N1)と、
実行命令に応答して上記有為差有無の判定および少なくとも断定できる寿命差の演算を行いその演算結果を上記表示装置の画面に出力する有為差判定・有為寿命差演算手順(N2)とを含む。
The significant difference presence / absence judgment / significant life difference estimation program (51) in the acceleration test of the present invention is a computer-executable program,
On the screen of the display device (2), as input information, the value of the Weibull slope of the Weibull distribution of the product under test, the number of each test in the two lots of Level 1 and Level 2 to be compared, and the life of the test result are entered. A prompt screen output procedure (N1) for prompting display;
In response to the execution command, the determination of the presence / absence of the significant difference and the calculation of the life difference that can be determined at least are performed, and the calculation result is output to the screen of the display device. including.

この有為差判定・有為寿命差演算手順(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)とを備える。
This significant difference judgment and significant life difference calculation procedure (N2)
A first random number generation procedure (N21) for generating a random number according to the Weibull distribution determined from the value of the Weibull slope of the input information and the temporarily set life, for the number of test pieces of the level 1 lot;
A second random number generation procedure (N22) for generating random numbers according to the same Weibull distribution as the first random number generation procedure (N21) for the test number of lots of level 2.
From the set number of random numbers generated in the first random number generation procedure (N21) by the predetermined lifetime calculation method and the set number of random numbers generated in the second random number generation procedure (N22), the predetermined lifetime calculation method. A life ratio calculation procedure (N23) for calculating a life ratio of a set of level 1 and level 2 lifespans from the obtained lifespan,
The first random number generation procedure (N21), the second random number generation procedure (N22), and the life ratio calculation procedure (N23) are repeated a set number of times to determine the life ratio of the set number of sets, and the probability distribution of the life ratio A cumulative probability distribution analysis procedure (N24) for creating a cumulative probability distribution and reading the maximum life ratio within the set confidence range from the cumulative probability distribution;
It is judged that there is a significant difference when the life ratio of the life of Level 1 and Level 2 in the above input information exceeds the maximum life ratio obtained by the cumulative probability distribution analysis procedure, and there is no significant difference when it is within the maximum life ratio. And a significant difference presence / absence determination procedure (N25).

有為差判定・有為寿命差演算手順(N2)は、さらに、有為差有無判定手順(N25)で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手順を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す倍率変更繰り返し手順(N26)と、
この倍率変更繰り返し手順(N26)で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する寿命倍率関係グラフ作成手順(N27)と、
上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする有為寿命差読み取り手順(28)と、
上記有為差有無判定手順で判定した結果および上記有為寿命差読み取り手順で読み取った、少なくとも断定できる寿命差を表示装置に表示させる結果出力手順(29)と、
を含む。
The significant difference determination / significant life difference calculation procedure (N2) is further executed when there is a significant difference in the significant difference presence / absence determination procedure (N25). And a cumulative change distribution analysis procedure that repeats the cumulative probability distribution analysis procedure until the set magnification reaches the set maximum magnification from the set minimum magnification,
A life magnification relationship graph creation procedure (N27) for creating a life magnification relationship graph, which is a graph showing the relationship between the set magnification repeatedly changed and the maximum life ratio within the reliability range, obtained in this magnification change repetition procedure (N26). )When,
Reading the value of the set magnification corresponding to the life ratio of the life of Level 1 and Level 2 in the input information from the life magnification relationship graph and reading the significant life difference as a life difference that can be determined at least. (28) and
A result output procedure (29) for displaying on the display device the life difference that can be determined at least, the result determined in the significant difference presence / absence determination procedure and the significant life difference reading procedure;
including.

この構成の加速試験における有為差有無判定・有為寿命差見積もりプログラムは、この発明の有為差有無判定・有為寿命差見積もりの実施に使用され、有為差の有無の判定、および有為差があると判定される場合に、少なくともどの程度の寿命差があると言えるのかを、簡単かつ迅速に試算することができ、かつ定量的に求められて信頼性の高いものとでき、熟練者でなくても必要な寿命差を適切に試算することができる。   The program for determining the presence / absence of significant difference / estimated life difference estimation in the accelerated test of this configuration is used to implement the existence / non-existence of significant difference / estimated life difference estimation of this invention. When it is determined that there is a difference, at least how much life difference can be estimated can be easily and quickly calculated, and it can be quantitatively obtained and reliable. Even if it is not a person, it can estimate the required life difference appropriately.

この発明の加速試験における有為差有無判定・有為寿命差見積もり方法、装置、およびプログラムによると、コンピュータシミュレーションにより、試験対象品のワイブル分布から乱数を発生させ、同じワイブル分布から乱数を発生させ、その発生させた乱数から求まる寿命を演算して1組の寿命比を算出し、この処理を設定回数繰り返して設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成し、上記累積確率分布から設定信頼幅内の最大および最小の寿命比を読み取り、有為差有無の判定を行うため、有為差有無の判定を、簡単かつ迅速に行うことができ、かつ定量的に求められて信頼性の高いものとでき、熟練者でなくても判定することができる。
また、有為差有りと判定した場合に、寿命が異なるワイブル分布に繰り返し毎に変えて、上記累積確率分布分析手順を繰り返し、繰り返しにより得られた設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成し、このグラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とするため、有為差がある場合に、少なくとも断定できる寿命差を、簡単かつ迅速に求めることができ、かつ定量的に求められて信頼性の高いものとでき、熟練者でなくても適切に有為寿命差を求めることができる。
According to the method, apparatus, and program for determining whether there is a significant difference in the accelerated test of the present invention, the random number is generated from the Weibull distribution of the test object by computer simulation, and the random number is generated from the same Weibull distribution. The life obtained from the generated random number is calculated to calculate one set of life ratio, this process is repeated a set number of times to obtain the life ratio of the set number of sets, and the probability distribution of this life ratio and the cumulative probability distribution And read the maximum and minimum life ratio within the set confidence range from the cumulative probability distribution above, and determine the presence or absence of a significant difference, so the presence or absence of a significant difference can be easily and quickly determined, In addition, it can be obtained quantitatively and can be highly reliable, and can be determined even if it is not an expert.
In addition, when it is determined that there is a significant difference, the cumulative probability distribution analysis procedure is repeated by changing the Weibull distribution with different lifetimes every time, and the set ratio obtained by the iterations and the maximum lifetime ratio within the confidence range are repeated. From the graph, the value of the set magnification corresponding to the life ratio of the life of level 1 and level 2 in the input information is read, and the read value is at least In order to make a life difference that can be determined, if there is a significant difference, at least the life difference that can be determined can be determined easily and quickly, and it can be quantitatively determined and highly reliable. Even if it is not, a significant difference in life can be obtained.

この発明の実施形態を説明する。この加速試験における有為差有無判定・有為寿命差見積もり方法は、軸受等の機械部品または試験片からなる試験対象品を使用環境よりも厳しい所定の環境条件におき、破損が発生した時間から寿命を算出する加速試験において、試験対象品の2つのロットの寿命から、有為差有無の判定、および少なくとも断定できる、倍率による寿命差を算出する方法である。   An embodiment of the present invention will be described. In this accelerated test, the existence of a significant difference and the estimation method of the significant life difference are determined by placing a test object consisting of a mechanical part such as a bearing or a test piece under specified environmental conditions that are severer than the operating environment, and starting from the time when the damage occurred. In the accelerated test for calculating the life, this is a method for calculating the life difference by magnification that can be determined and at least determined from the life of two lots of the test object.

以下、この実施形態を図面と共に説明する。この加速試験における有為差有無判定・有為寿命差見積もり方法は、図1に示すコンピュータ1に、乱数シミュレーションプログラムである加速試験における有為差有無判定・有為寿命差見積もりプログラム41を実行させることで行う。コンピュータ1はパーソナルコンピュータ等からなり、中央処理装置4およびメモリ5を有し、所定のオペレーションシステムによって動作するものである。コンピュータ1には、液晶表示装置等の画面によって表示可能な表示装置2と、キーボードやマウス等の入力装置3が接続され、あるいは付属して設けられている。コンピュータ1、表示装置2、入力装置3、および必要寿命差見積もりプログラム41により、図2に各機能達成手段をブロックで示した有為差有無判定・有為寿命差見積もり装置が構成される。同図の有為差有無判定・有為寿命差見積もり装置の構成については、後に説明する。
有為差有無判定・有為寿命差見積もりプログラム51は、図4および図5に流れ図で示す手順を備えるものである。同図の内容は、後に説明する。
Hereinafter, this embodiment will be described with reference to the drawings. In the acceleration test, the computer 1 shown in FIG. 1 executes a significant difference presence / absence determination / significant life difference estimation program 41 in the acceleration test, which is a random number simulation program. Do that. The computer 1 is composed of a personal computer or the like, has a central processing unit 4 and a memory 5, and operates by a predetermined operation system. The computer 1 is provided with a display device 2 that can be displayed on a screen such as a liquid crystal display device and an input device 3 such as a keyboard and a mouse. The computer 1, the display device 2, the input device 3, and the necessary life difference estimation program 41 constitute a significant difference presence / absence determination / significant life difference estimation device in which each function achievement means is shown as a block in FIG. 2. The configuration of the significant difference presence / absence determination / significant life difference estimation apparatus shown in FIG.
The significant difference presence / absence determination / significant life difference estimation program 51 includes the procedures shown in the flowcharts of FIGS. 4 and 5. The contents of this figure will be described later.

この必要寿命差見積もり方法は、図3に示すように、コンピュータ1に対して所定の情報を入力する入力過程M1と、コンピュータ1で演算処理を行って演算結果を出力するコンピュータ演算処理過程M2とからなる。   As shown in FIG. 3, the required life difference estimation method includes an input process M1 for inputting predetermined information to the computer 1, and a computer calculation process M2 for performing calculation processing on the computer 1 and outputting a calculation result. Consists of.

入力過程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の入力、および試験結果の寿命を促す表示が行われ、その個数および寿命を入力する記入用表示が行われる。オペレータは、この記入用表示に対して該当する情報を入力する。
In the input process M1, as shown in FIG. 6, an input screen 2a for prompting input of predetermined input information is displayed on the display device 2 by the output of the computer 1, and a display for prompting predetermined input is displayed on the input screen 2a. Is called.
The input screen 2a is displayed by switching to a plurality of screens in the example of FIG. Of the input screens, the first screen (FIG. 6A) displays whether to select a significant difference between the L10 life and the L50 life. By selecting the display portion consisting of, either the L10 life or the L50 life can be selected.
On the second screen (FIG. 6B), as input information, a display for prompting input of the value of the Weibull slope of the Weibull distribution of the test object and a display for filling in the input are performed. The operator enters a corresponding value in this entry display.
On the second screen (FIG. 6 (C)), the test numbers n1 and n2 of the two lots of level 1 and level 2 to be compared are input and a display for prompting the test result life is performed. A display for entry is entered. The operator inputs information corresponding to the entry display.

ワイブルスロープの値には、試験の実績値を入力することが好ましい。実績値は10個以上の試験で得た結果を用いることが望ましく、より好ましくは20個以上の試験結果である。また、試験個数および試験結果は、実際に行った試験における個数および結果を入力する。なお、試験個数と寿命は、水準2の欄に長寿命側の入力をしなければないらいとする入力規則が、有為差有無判定・有為寿命差見積もりプログラム51に手順として設定されている。この手順は、例えば促し画面出力手順N1の一部として設けられる。以上の点に注意して、所定の条件を入力し、入力画面2aのOKボタンを押すと計算が開始される。   It is preferable to input the actual value of the test as the value of the Weibull slope. As the actual value, it is desirable to use a result obtained by 10 or more tests, and more preferably 20 or more test results. As the number of tests and test results, the number and results of tests actually performed are input. Note that the number of tests and the lifespan are set in the Level 2 column as a procedure in the program for judging whether there is a significant difference or a significant lifespan difference estimation program 51 in which the long life side must be entered. . This procedure is provided as a part of the prompt screen output procedure N1, for example. With the above points in mind, calculation is started when a predetermined condition is input and the OK button on the input screen 2a is pressed.

図3のコンピュータ演算処理過程M2では、入力されたワイブルスロープの値、および2つのロットの各試験個数n1,n2、試験結果の寿命から、有為差有無の判定に必要な寿命差を計算して、有為差有無を判定し、有為差有りの場合には、少なくとも断定できる寿命差を計算する。   In the computer processing step M2 in FIG. 3, the life difference necessary for determining the existence of a significant difference is calculated from the input Weibull slope value, the test numbers n1, n2 of the two lots, and the life of the test result. Whether there is a significant difference or not, if there is a significant difference, calculate at least the life difference that can be determined.

計算が終了すると、図7に示す出力画面に結果が表示される。有為差がない場合は、同図(A)のように、有為差無しという表示と、有為差有りと判断するための必要寿命差が表示される。
有為差がある場合は、同図(B)のように、有為差ありという表示と、少なくとも断定できる寿命差(倍率)が表示される。
When the calculation is completed, the result is displayed on the output screen shown in FIG. When there is no significant difference, as shown in FIG. 6A, a display indicating that there is no significant difference and a required life difference for determining that there is a significant difference are displayed.
When there is a significant difference, as shown in FIG. 5B, a display that there is a significant difference and a life difference (magnification) that can be determined at least are displayed.

図1の有為差有無判定・有為寿命差見積もりプログラム51は、コンピュータ1で実行可能なプログラムであって、図4,図5に流れ図で示す手順を備える。図4に示すように、有為差有無判定・有為寿命差見積もりプログラム51は、促し画面出力手順N1と必要寿命差演算手順N2とでなり、促し画面出力手順N1では、図6と共に前述した入力画面2aを出力する。この入力画面2aに対して、上記各入力情報が入力手段3から入力され、かつ入力画面2aのOKキーのクリック等によって実行命令が入力手段3から入力されると、有為差有無判定・有為寿命差見積もり手順N2が実行される。同図の入力画面に対して入力する過程が、図3の入力過程M1であり、同図のコンピュータ演算処理過程M2は図4の有為差有無判定・有為寿命差見積もり手順N2を実行する過程である。   The significant difference presence / absence determination / significant life difference estimation program 51 of FIG. 1 is a program that can be executed by the computer 1 and includes the procedures shown in the flowcharts of FIGS. As shown in FIG. 4, the significant difference presence / absence determination / significant life difference estimation program 51 includes a prompt screen output procedure N1 and a required life difference calculation procedure N2. The prompt screen output procedure N1 is described above with reference to FIG. The input screen 2a is output. When the input information is input to the input screen 2a from the input means 3 and an execution command is input from the input means 3 by clicking the OK key on the input screen 2a, etc., it is determined whether there is a significant difference. Therefore, the life difference estimation procedure N2 is executed. The input process M1 in FIG. 3 is the input process M1 in FIG. 3, and the computer calculation process M2 in FIG. 4 executes the significant difference presence / absence / significant life difference estimation procedure N2 in FIG. It is a process.

有為差有無判定・有為寿命差見積もり手順N2は、図5に流れ図で示す各手順で構成される。この流れ図には各手順M21〜M29毎の具体的な処理例を注釈として併記してある。
理解の容易のため、具体的数値を上げ、同図の具体的処理例を参照して、各手順N21〜N29を説明する。
The significant difference presence / absence determination / significant life difference estimation procedure N2 includes the respective steps shown in the flowchart of FIG. In this flowchart, a specific processing example for each of the procedures M21 to M29 is shown as an annotation.
For easy understanding, each of the procedures N21 to N29 will be described with reference to a specific processing example in FIG.

今、ある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)。
Assume that two lots of tests are performed in an accelerated test with a Weibull distribution having a certain L10 life (FIG. 8A) and a Weibull slope of 4. Assume that the number of tests n1 in the first lot (level 1) is 3, and the number of tests n2 in the second lot (level 2) is also 3.
First, according to the Weibull distribution, three random numbers corresponding to the test number n1 of the first lot are generated, and L10 life and L50 life are calculated from the three data (N21). As a calculation method, an appropriate method used for calculating the life from the result of the conventional acceleration test is used. Although only the L10 life or only the L50 life may be calculated, in this embodiment, in order to calculate the necessary life difference for each of the L10 life and the L50 life, both the L10 life and the L50 life are calculated. ing.
Next, from the Weibull distribution with the same L10 life, three random numbers that are the test number n2 of the second lot are generated, and the L10 life and L50 life are calculated from the three data by the same method as above. (N22).

次に、得られた1組のL10寿命とL50寿命の寿命比をそれぞれ算出する。すなわちL10寿命同士の寿命比、およびL50寿命同士の寿命比を算出する(N23)。   Next, the respective life ratios of the obtained L10 life and L50 life are calculated. That is, the life ratio between the L10 lives and the life ratio between the L50 lives are calculated (N23).

次に、これらの作業を設定回数(例えば1000回)繰り返し、設定回数組(1000組)のL10およびL50寿命の寿命比の確率分布(図8(B)と累積確率分布(図8(C)を作成する(N24)。
これらは、同じ寿命分布を持つロットから3個の試験片を2組抜きとり、寿命試験を行った結果、どの程度の寿命差が出るのかを設定回数(1000回)確認することに対応している。
Next, these operations are repeated a set number of times (for example, 1000 times), and the probability distribution (FIG. 8B) and the cumulative probability distribution (FIG. 8C) of the life ratio of the L10 and L50 life of the set number of sets (1000 sets). Is created (N24).
These correspond to checking the set number of times (1000 times) to see how much difference in life results when two sets of three test pieces are extracted from the lot with the same life distribution and the life test is conducted. Yes.

度数分布の図(図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倍以上であるといえる。
From the frequency distribution diagram (FIG. 8B), it can be seen that the life varies even though the test pieces are extracted from the lots having the same life distribution. If the variation in the life ratio is defined as a predetermined reliability interval, for example, a 5% and 95% interval (90% confidence interval) of the cumulative probability distribution, the life variation is 0 as shown in FIG. It can be calculated to be 39 to 2.57 times (N243).
As a result, even when two sets of three test pieces were extracted from lots having the same life distribution, the life ratio varied between 0.39 and 2.57 times for 90%. It shows that it will end. In other words, if there is a life difference of 2.57 times or more, which is the maximum value within the variation, it is highly possible that the life distribution between the lots is different. When the number of tests is three, it can be said that the life difference that can determine the significant difference is 2.57 times or more.

そこで、上記ばらつき内の最大値である2.57倍と、入力された2つのロット(水準1,水準2)の寿命差(倍率)、すなわち〔(水準2の寿命)/(水準1の寿命)〕とを比較し(N25)、入力された2つのロットの寿命の寿命差の方が小さいと、有為差無しと反対してその判定結果と上記の有為差を判断できる寿命差の値(上記の例では2.57場合)とを出力画面に表示させ(N29)、コンピュータ1による処理を終了する。   Therefore, the maximum value within the above-mentioned variation is 2.57 times and the life difference (magnification) between the two input lots (level 1 and level 2), that is, [(life of level 2) / (life of level 1). )] (N25), and if the difference in the lifespan of the two input lots is smaller, the difference in lifespan that can be used to judge the difference between the judgment result and the above-mentioned The value (2.57 in the above example) is displayed on the output screen (N29), and the processing by the computer 1 is terminated.

入力された2つのロットの寿命の寿命差の方が大きい場合は、次の手順N26〜N28に進み、少なくとも断定できる寿命差の演算を行う。
手順N26では、寿命分布が異なるワイブル分布との組み合わせで、上述と同様の計算を行う。具体的には、ワイブル分布と寿命が1.1 、1.2 …50倍異なるワイブル分布の組み合わせで計算を行う。このようにして計算を行うと、組み合わせが1.1 、1.2 …50倍と変化するごとに、寿命差が1.1 、1.2 …50倍を中心にしてばらつくはずである。
When the life difference between the two input lots is larger, the process proceeds to the next steps N26 to N28, and at least the life difference that can be determined is calculated.
In procedure N26, the same calculation as described above is performed in combination with a Weibull distribution having a different lifetime distribution. Specifically, the calculation is performed with a combination of the Weibull distribution and the Weibull distribution whose lifespan is 1.1, 1.2... 50 times different. When the calculation is performed in this way, every time the combination changes 1.1, 1.2... 50 times, the life difference should vary around 1.1, 1.2... 50 times.

最後に、組み合わせの倍率を横軸にして、それぞれの累積確率分布の90%区間を縦軸にとったグラフ(図8(D))を作成する(N27)。この図は、寿命がある倍率で異なる2ロットの試験片群から3個の試験片を1組ずつ抜き取って寿命比較を行った場合、どの程度寿命差がばらつくかを示した図である。例えば、4倍の2ロットの試験片群から2個の試験片を1組ずつ抜き取り、寿命比較を行った場合、寿命差は4倍になるケースが最も多いものの、7.5倍の寿命差になるケースも10%の確率で起こりえることがこの図から読み取れる。   Finally, a graph (FIG. 8D) is created (N27) with the combination magnification as the horizontal axis and the 90% interval of each cumulative probability distribution as the vertical axis. This figure shows how the life difference varies when three sets of test pieces are extracted one by one from a group of two different test pieces at different magnifications and the life comparison is performed. For example, when two pairs of test pieces are extracted from a group of 2 lots of test specimens 4 times each and the life comparison is performed, the life difference is most often 4 times, but the life difference is 7.5 times. It can be seen from this figure that the case that becomes can occur with a probability of 10%.

今、問題としている2ロット間の寿命差は3倍である。図から、3倍の寿命差は寿命差が1.7倍以上ある2ロット間でしか起こり得ない寿命差だということが読み取れるので、少なくともいえる寿命差は1.7倍ということができる。このように読み取った値を、少なくとも断定できる寿命差と定める(N28)。
この定めた寿命差の値は、有為差有りという表示と共に、表示装置2の出力画面2b(図7(A))に出力する(N29)。
以上が2ロットの寿命から少なくともいえる寿命差を算出する手順である。
Now, the difference in life between the two lots in question is three times. From the figure, it can be read that the life difference of 3 times is a life difference that can only occur between two lots having a life difference of 1.7 times or more, so at least the life difference that can be said can be 1.7 times. The value read in this way is determined as at least a life difference that can be determined (N28).
The determined life difference value is output to the output screen 2b (FIG. 7A) of the display device 2 together with a display that there is a significant difference (N29).
The above is the procedure for calculating at least the life difference that can be said from the life of two lots.

図4,図5に示し有為差有無判定・有為寿命差見積もりプログラム41についての上記の説明は、具体的に数値を例にとって説明したが、この有為差有無判定・有為寿命差見積もりプログラム41は、整理すると、次の手順により構成される。   The above description of the significant difference presence / absence determination / significant life difference estimation program 41 shown in FIGS. 4 and 5 has been described with reference to specific numerical values. This significant difference presence / absence determination / significant life difference estimation When organized, the program 41 is constituted by the following procedure.

この発明の加速試験における有為差有無判定・有為寿命差見積もりプログラム(51)は、コンピュータで実行可能なプログラムであって、
上記表示装置(2)の画面に、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命の入力を促す表示を行わせる促し画面出力手順(N1)と、
実行命令に応答して上記有為差有無の判定および少なくとも断定できる寿命差の演算を行いその演算結果を上記表示装置の画面に出力する有為差判定・有為寿命差演算手順(N2)とを含む。
The significant difference presence / absence judgment / significant life difference estimation program (51) in the acceleration test of the present invention is a computer-executable program,
On the screen of the display device (2), as input information, the value of the Weibull slope of the Weibull distribution of the product under test, the number of each test in the two lots of Level 1 and Level 2 to be compared, and the life of the test result are entered. A prompt screen output procedure (N1) for prompting display;
In response to the execution command, the determination of the presence / absence of the significant difference and the calculation of the life difference that can be determined at least are performed, and the calculation result is output to the screen of the display device. including.

この有為差判定・有為寿命差演算手順(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)とを備える。
This significant difference judgment and significant life difference calculation procedure (N2)
A first random number generation procedure (N21) for generating a random number according to the Weibull distribution determined from the value of the Weibull slope of the input information and the temporarily set life, for the number of test pieces of the level 1 lot;
A second random number generation procedure (N22) for generating random numbers according to the same Weibull distribution as the first random number generation procedure (N21) for the test number of lots of level 2.
From the set number of random numbers generated in the first random number generation procedure (N21) by the predetermined lifetime calculation method and the set number of random numbers generated in the second random number generation procedure (N22), the predetermined lifetime calculation method. A life ratio calculation procedure (N23) for calculating a life ratio of a set of level 1 and level 2 lifespans from the obtained lifespan,
The first random number generation procedure (N21), the second random number generation procedure (N22), and the life ratio calculation procedure (N23) are repeated a set number of times to determine the life ratio of the set number of sets, and the probability distribution of the life ratio A cumulative probability distribution analysis procedure (N24) for creating a cumulative probability distribution and reading the maximum life ratio within the set confidence range from the cumulative probability distribution;
It is determined that there is a significant difference when the life ratio of the life of Level 1 and Level 2 in the above input information exceeds the maximum life ratio obtained by the cumulative probability distribution analysis procedure, and there is no significant difference when it is within the maximum life ratio. And a significant difference presence / absence determination procedure (N25) for determination.

有為差判定・有為寿命差演算手順(N2)は、さらに、有為差有無判定手順(N25)で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手順を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す倍率変更繰り返し手順(N26)と、
この倍率変更繰り返し手順(N26)で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する寿命倍率関係グラフ作成手順(N27)と、
上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする有為寿命差読み取り手順(28)と、
上記有為差有無判定手順(N28)で判定した結果および上記有為寿命差読み取り手順で読み取った、少なくとも断定できる寿命差を表示装置に表示させる結果出力手順(29)と、を含む。
The significant difference determination / significant life difference calculation procedure (N2) is further executed when there is a significant difference in the significant difference presence / absence determination procedure (N25). And a cumulative change distribution analysis procedure that repeats the cumulative probability distribution analysis procedure until the set magnification reaches the set maximum magnification from the set minimum magnification,
A life magnification relationship graph creation procedure (N27) for creating a life magnification relationship graph, which is a graph showing the relationship between the set magnification repeatedly changed and the maximum life ratio within the reliability range, obtained in this magnification change repetition procedure (N26). )When,
Reading the value of the set magnification corresponding to the life ratio of the life of Level 1 and Level 2 in the input information from the life magnification relationship graph and reading the significant life difference as a life difference that can be determined at least. (28) and
And a result output procedure (29) for displaying on the display device at least a determinable life difference read by the significant difference presence / absence determination procedure (N28) and the significant life difference reading procedure.

上記累積確率分布分析手順(N24)は、手順N21〜N23を設定回数繰り返させる手順N241と、その繰り返し過程で記憶したおいた寿命比から確率分布および累積確率分布を演算する手順N242と、この累積確率分布から90%範囲の大寿命比を読み取る手順N243とからなる。   The cumulative probability distribution analysis procedure (N24) includes a procedure N241 for repeating steps N21 to N23 a set number of times, a procedure N242 for calculating a probability distribution and a cumulative probability distribution from the life ratio stored in the repetition process, And a procedure N243 for reading a long life ratio in the 90% range from the probability distribution.

上記倍率変更繰り返し手順(N26)は、設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す手順N261と、繰り返し毎に異なる設定倍率に変更する倍率変更手順N231とからなる。   The magnification change repetition procedure (N26) includes a procedure N261 that is repeated until the set magnification is changed from the set minimum magnification to the set maximum magnification, and a magnification change procedure N231 that is changed to a different set magnification for each repetition.

各乱数発生手順(N21,N22)の詳細について説明する。これらの手順N21,N22は、ワイブル分布を特定し、その特定したワイブル分布に従ってワイブル乱数を発生させる。
一般に軸受の寿命分布は次式1)のワイブル分布に従うと言われている。
Details of each random number generation procedure (N21, N22) will be described. These procedures N21 and N22 specify a Weibull distribution and generate a Weibull random number according to the specified Weibull distribution.
Generally, it is said that the bearing life distribution follows the Weibull distribution of the following equation 1).

Figure 2008128698
Figure 2008128698

ただし、m:ワイブルスロープ、α:尺度因子、γ:最小寿命、
ワイブル分布は、3つのパラメータを持っており、ワイブルスロープmによって指数分布、対数正規分布、正規分布を表現できる万能分布として知られている。参考として、図10に各種パラメータを変化させた時のワイブル分布の変化を示す。ワイブルスロープmは、分布の形状を支配するパラメータであり、この値が小さいほどばらつきの大きい分布ということができる。尺度因子αは、横軸(寿命)のスケールを変化させるもので、この値が大きいほど寿命は相対的に長くなる。最小寿命γは、寿命分布の横軸(寿命)を単にシフトさせるものである。
Where m: Weibull slope, α: scale factor, γ: minimum life,
The Weibull distribution has three parameters and is known as a universal distribution that can express an exponential distribution, a lognormal distribution, and a normal distribution by the Weibull slope m. For reference, FIG. 10 shows changes in the Weibull distribution when various parameters are changed. The Weibull slope m is a parameter that governs the shape of the distribution, and it can be said that the smaller the value, the larger the variation. The scale factor α changes the scale of the horizontal axis (lifetime), and the larger the value, the longer the life. The minimum life γ simply shifts the horizontal axis (life) of the life distribution.

この実施形態では、ワイブル乱数を発生させるが、この乱数を発生させるためにはワイブル分布の3つのパラメータを決定する必要がある。決め方の手順は、例えば以下のようになる。
1) ワイブルスロープmを実績から決定する。
2) 乱数を発生させたい分布の信頼度(例えばL10寿命であるか、あるいはL50寿 命であるか)を決定する。
3) 信頼度から求めたワイブルスロープmから、最小寿命γを所定の数式を使って決定 する。例えば、L10寿命またはL50寿命から求めた尺度因子αから、
最小寿命γを、例えば、以下の2)式を使って決定する。
この式は、1990年制定のISOの最小寿命であり、実験値からの回帰式である。
In this embodiment, a Weibull random number is generated. In order to generate this random number, it is necessary to determine three parameters of the Weibull distribution. The procedure for deciding is, for example, as follows.
1) The Weibull slope m is determined from the results.
2) Determine the reliability of the distribution for which random numbers are to be generated (for example, whether the life is L10 or L50).
3) From the Weibull slope m obtained from the reliability, the minimum life γ is determined using a predetermined formula. For example, from the scale factor α obtained from the L10 life or L50 life,
The minimum lifetime γ is determined using, for example, the following equation 2).
This equation is the minimum life of ISO established in 1990, and is a regression equation from experimental values.

Figure 2008128698
Figure 2008128698

これは、R≦10の値で、R=0(L10寿命でのa1)のとき、この式は1になるという式である。過去のISOの最少寿命考慮の式では、L10寿命以下の寿命は、この式にL10寿命を書けた値ということで定義されている。Rは信頼度に対応する値(100−Rが信頼度となる値)である。
なお、最小寿命の定め方については、各種の規格(例えばISO)において、時代と共に変更される場合があるが、規格の変更に伴い、実施時の規格に応じた定め方を採用すれば良い。また、最小寿命は、材料試験条件によっても変化するのでより一般的な式で記述するほうが良いとの主張もあり、適宜の値を用いれば良い。
This is an expression such that this value becomes 1 when R ≦ 10 and R = 0 (a1 with L10 life). In a formula that takes into account the minimum lifetime of ISO in the past, a lifetime that is less than or equal to the L10 lifetime is defined as a value that can write the L10 lifetime in this formula. R is a value corresponding to the reliability (a value at which 100-R is the reliability).
Note that the method of determining the minimum life may be changed with the times in various standards (for example, ISO), but a method according to the standard at the time of implementation may be adopted along with the change of the standard. In addition, there is a claim that it is better to describe the minimum life by a more general formula because it changes depending on the material test conditions, and an appropriate value may be used.

ワイブル乱数の発生につき説明する。乱数とは、定性的にはでたらめな数列であって、発生頻度が均一(等確率)で、その発生に規則性がない(無規則性)というものであるが、完全な乱数を発生させることは不可能である。そこで、コンピュータで発生させることのできる疑似乱数を使う。簡易な乱数発生アルゴリズムでは、例えば10進法で20桁ぐらいの周期性が見られるが、周期性が6千桁以上の周期性となるものもあり、このような周期性の少ない乱数発生アルゴリズムを用いることが好ましい。   The generation of Weibull random numbers will be described. Random numbers are qualitatively random sequences that have a uniform frequency of occurrence (equal probability) and no regularity in their generation (irregularity). Is impossible. Therefore, pseudorandom numbers that can be generated by a computer are used. A simple random number generation algorithm, for example, has a periodicity of about 20 digits in decimal notation, but some periodicity has a periodicity of 6,000 digits or more. It is preferable to use it.

この実施形態では、一様な乱数ではなく、ワイブル分布に従った乱数であるワイブル乱数を発生させる。このため発生方法には工夫が必要になる。確率密度関数が複雑な場合、その分布に従う乱数を発生するには棄却法と呼ばれる方法を用いればよく、この実施形態においても、棄却法を用いる。
確率密度関数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を採用し、確率密度分布に従う乱数の数列を作っていく方法を棄却法という。この方法は、条件に合わない乱数を捨てることになるので乱数発生法としては効率がよくないが、よい一様乱数さえ得られれば原理的に正しい数列が得られる方法である。
In this embodiment, a Weibull random number that is not a uniform random number but a random number according to the Weibull distribution is generated. For this reason, a device is required for the generation method. When the probability density function is complicated, a method called a rejection method may be used to generate random numbers according to the distribution, and the rejection method is also used in this embodiment.
As shown in FIG. 11, it is assumed that the domain of the probability density function f (x) is in the range of 0 to X0, and the maximum value of f (x) in the domain is M. If RN is a uniform pseudorandom number in the interval [0,1], a uniform pseudorandom number xi in the interval [0, x0] can be generated by X0 · RN. Similarly, uniform pseudorandom numbers yi in the interval [0, M] can be generated by M · RN. Therefore, when the generated random numbers xi, yi satisfy the condition that f (xi)> yi, the random number xi is adopted as following the given probability density distribution. The random number xi is not adopted. A method of repeating this work, adopting random numbers xi with a probability according to the probability density distribution, and creating a sequence of random numbers according to the probability density distribution is called a rejection method. This method is not efficient as a random number generation method because random numbers that do not meet the conditions are discarded. However, in principle, if a good uniform random number is obtained, a correct number sequence can be obtained.

図5の乱数発生手順N21,N22において、ワイブル分布の特定のための寿命(L10寿命)は、適宜想定した値を、これらの手順N21,N22の計算式等に設定しておいても良く、また入力手段3からの入力によって可変としても良い。
また、寿命比較算出手順N23において、発生したワイブル乱数からのL10寿命を求める演算は、加速試験において従来から用いられている中の適宜の寿命計算方法を用いれば良い。加速試験においては、破損の発生した時間からワイブルプロットで寿命を求めることが行われており、このような破損時間から寿命を求める計算方法を、上記寿命比算出手順における寿命計算方法として用いる。
In the random number generation procedures N21 and N22 in FIG. 5, the life for specifying the Weibull distribution (L10 life) may be set as an appropriately assumed value in the calculation formula of these procedures N21 and N22, Further, it may be variable by input from the input means 3.
Further, in the life comparison calculation procedure N23, the calculation for obtaining the L10 life from the generated Weibull random number may use an appropriate life time calculation method conventionally used in the acceleration test. In the accelerated test, the life is obtained by a Weibull plot from the time when the breakage occurs, and a calculation method for obtaining the life from the breakage time is used as the life calculation method in the life ratio calculation procedure.

この寿命計算には、例えば次の方法が採用できる。
(1) 寿命試験を実施する。
(2) 得られたデータ(破損した時間あるいは破損した負荷回数)を昇順に並び替える。
(3) これらデータを図13のグラフ(ワイブル確率紙)にプロットする(縦軸:累積破損確率、横軸:寿命)。
(4) 図13の紙にプロットしたデータの最適直線を最小二乗法で引く。このとき、L10 寿命以下の位置に最小寿命があるということになるので、L10 寿命の値を10分割し(何分割でも良いがフィッティングでの計算時間が妥当な時間になるように設定する)、累積確率0 %の位置にプロットを加える。10通りの最適曲線で最もデータがフィットする最適直線を採用する。
(5) そうすると、ワイブルスロープがこの線の傾き、最小寿命は、L10 寿命の値を10分割のいずれかの値、L10 寿命(ワイブルスロープが累積確率10%交わる寿命)と尺度因子αの関係からαを決定できる。
For this lifetime calculation, for example, the following method can be employed.
(1) Conduct a life test.
(2) Sort the obtained data (damaged time or number of broken loads) in ascending order.
(3) These data are plotted on the graph (Weibull probability paper) in FIG. 13 (vertical axis: cumulative failure probability, horizontal axis: life).
(4) The optimal straight line of the data plotted on the paper of FIG. 13 is drawn by the method of least squares. At this time, since there is a minimum life at a position below the L10 life, the value of the L10 life is divided into 10 (any number is possible, but the calculation time for fitting is set to be a reasonable time) A plot is added at the position where the cumulative probability is 0%. The optimum straight line that best fits the data is adopted among the 10 optimum curves.
(5) Then, the Weibull slope is the slope of this line, and the minimum life is based on the relationship between the L10 life value in any of 10 divisions, the L10 life (lifetime where the Weibull slope has a cumulative probability of 10%) and the scale factor α. α can be determined.

この打切り時間見積もり装置プログラム51では、有為差のある寿命差有りと判断するための必要寿命差の関係を求めている。図9に寿命差有りと判断するための必要寿命差の関係を示す。この結果は、事前に入力するワイブルスロープと試験個数によって変化するので、図にはワイブルスロープと試験個数を変化させた結果を示している。試験個数が増えるほど、あるいはワイブルスロープが大きいほど、傾きが小さくなることが分かる。これらの結果は、(1) 試験個数が増えるほど寿命差有りと判断するための必要寿命差を小さくできること、(2) ワイブルスロープが大きい試験条件ほど寿命差有りと判断するための必要寿命差を小さくできることを示している。   The censoring time estimation apparatus program 51 obtains the relationship of the necessary life difference for determining that there is a life difference with a significant difference. FIG. 9 shows the relationship of the required life difference for determining that there is a life difference. Since this result changes depending on the Weibull slope and the number of tests input in advance, the figure shows the result of changing the Weibull slope and the number of tests. It can be seen that the inclination decreases as the number of tests increases or the Weibull slope increases. These results show that (1) the difference in life required to determine that there is a difference in life can be reduced as the number of tests increases, and (2) the difference in life required to determine that there is a difference in life for test conditions with a larger Weibull slope. It shows that it can be made smaller.

図2と共に加速試験における有為差有無判定・有為寿命差見積もり装置につきする。この有為差有無判定・有為寿命差見積もり装置は、演算処理装置1と、この演算処理装置1の出力を画面に表示する表示装置2と、上記演算処理装置1に入力を行う入力手段3とを備える。   Together with FIG. 2, a device for determining whether there is a significant difference in the acceleration test and estimating a significant life difference will be described. This significant difference presence / absence determination / significant life difference estimation device includes an arithmetic processing device 1, a display device 2 for displaying the output of the arithmetic processing device 1 on a screen, and an input means 3 for inputting to the arithmetic processing device 1. With.

上記演算処理装置1は、上記表示装置2の画面に、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命の入力を促す表示を行わせる促し画面出力手段7Mと、実行命令に応答して上記有為差有無の判定および少なくとも断定できる寿命差の演算を行いその演算結果を上記表示装置2の画面に出力する有為差判定・有為寿命差演算手段51とを備える。   The arithmetic processing device 1 displays, on the screen of the display device 2, as input information, the value of the Weibull slope of the Weibull distribution of the test object, the number of tests in each of the two lots of level 1 and level 2 to be compared, and the test results. A prompt screen output means 7M for performing a display prompting the input of the life of the user, and in response to the execution command, the presence / absence of the significant difference is determined, and at least the life difference that can be determined is calculated, and the calculation result is displayed on the screen of the display device 2 And a significant life difference calculating / meaning life difference calculating means 51.

上記有為差判定・有為寿命差演算手段52は、第1乱数発生手段53、第2乱数発生手段54、寿命比算出手段55、累積確率分布分析手段56、有為差有無判定手段57、倍率変更繰り返し手段58、寿命倍率関係グラフ作成手段59、有為寿命差読み取り手段60、および結果出力手段61を有する。   The significant difference determination / significant life difference calculation means 52 includes a first random number generation means 53, a second random number generation means 54, a life ratio calculation means 55, a cumulative probability distribution analysis means 56, a significant difference presence / absence determination means 57, A magnification change repeating unit 58, a life magnification relationship graph creating unit 59, a significant life difference reading unit 60, and a result output unit 61 are provided.

第1乱数発生手段53は、上記入力情報のワイブルスロープの値および仮設定寿命から定まるワイブル分布に従った乱数を水準1のロットの試験個数分発生させる手段であり、図5の手順N21で説明した処理を行う。
第2乱数発生手段54は、第1乱数発生手段53と同じワイブル分布に従った乱数を水準2のロットの試験個数分発生させる手段であり、手順N22で説明した処理を行う。
寿命比算出手段55は、第1乱数発生手段53で発生させた設定個数の乱数から所定の寿命計算方法で求まる寿命と、第2乱数発生手段54で発生させた設定個数の乱数から上記所定の寿命計算方法で求まる寿命とから、1組の水準1,水準2の寿命の寿命比を算出する手段であり、手順N23で説明した処理を行う。
The first random number generating means 53 is means for generating random numbers according to the Weibull distribution determined from the value of the Weibull slope of the input information and the temporarily set life for the number of test pieces of the level 1 lot, and will be described in the procedure N21 of FIG. Perform the process.
The second random number generation means 54 is means for generating random numbers according to the same Weibull distribution as the first random number generation means 53 for the test number of lots of level 2, and performs the processing described in step N22.
The life ratio calculating means 55 is configured to calculate the predetermined time from the life determined by a predetermined life calculation method from the set number of random numbers generated by the first random number generating means 53 and the set number of random numbers generated by the second random number generating means 54. It is a means for calculating the life ratio of a set of level 1 and level 2 lifespans from the lifespan determined by the lifespan calculation method, and performs the processing described in step N23.

累積確率分布分析手段56は、これら第1乱数発生手段53,第2乱数発生手段54,および寿命比算出手段55の処理を設定回数繰り返し、上記設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成して累積確率分布から、設定信頼幅内の最大の寿命比を読み取る手段であり、手順N24で説明した処理を行う。   The cumulative probability distribution analysis means 56 repeats the processes of the first random number generation means 53, the second random number generation means 54, and the life ratio calculation means 55 a set number of times, obtains the life ratio of the set number of times, and this life This is means for creating a probability distribution of ratios and a cumulative probability distribution and reading the maximum life ratio within the set confidence width from the cumulative probability distribution, and performs the processing described in step N24.

有為差有無判定手段57は、上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手段56で得た最大寿命比を超える場合に有為差有りと判定し最大寿命比以内の場合に有為差無しと判定する手段であり、手順N25で説明した処理を行う。   The significant difference presence / absence determination means 57 determines that there is a significant difference when the life ratio of the life of the level 1 and level 2 in the input information exceeds the maximum life ratio obtained by the cumulative probability distribution analysis means 56, and determines the maximum life. A means for determining that there is no significant difference when the ratio is within the ratio, and performs the processing described in step N25.

倍率変更繰り返し手段58は、有為差有無判定手段57で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手段56の処理を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す手段であり、手順N26で説明した処理を行う。
寿命倍率関係グラフ作成手段59は、倍率変更繰り返し手段58で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する手段であり、手順N27で説明した処理を行う。
有為寿命差読み取り手段60は、上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする手段であり、手順N28で説明した処理を行う。
The magnification change repetition unit 58 is executed when the significant difference presence / absence determination unit 57 determines that there is a significant difference, and sequentially changes the Weibull distribution to a Weibull distribution having a set magnification with a different temporary setting life for each repetition. The process of the cumulative probability distribution analysis unit 56 is a unit that repeats the above-described setting magnification from the set minimum magnification to the set maximum magnification, and performs the processing described in step N26.
The life magnification relationship graph creating means 59 is a means for creating a life magnification relationship graph which is a graph showing the relationship between the set magnification repeatedly changed and the maximum life ratio within the reliability range obtained by the magnification change repetition means 58. Yes, the process described in step N27 is performed.
The significant life difference reading means 60 reads the value of the set magnification corresponding to the life ratio of the life of Level 1 and Level 2 in the input information from the life magnification relationship relation graph, and the life that can at least determine the read value. It is a means for making a difference, and the processing described in step N28 is performed.

結果出力手段61は、上記有為差有無判定手段57で判定した結果、および上記有為寿命差読み取り手段60で読み取った、少なくとも断定できる寿命差を表示装置2に表示させる手段であり、手順N29で説明した処理を行う。   The result output means 61 is a means for causing the display device 2 to display the result determined by the significant difference presence / absence determining means 57 and the life difference that can be determined at least read by the significant life difference reading means 60, and the procedure N29. Perform the process described in.

この発明の一実施形態に係る有為差有無判定・有為寿命差見積もり装置の概略ブロック図である。1 is a schematic block diagram of a significant difference presence / absence determination / significant life difference estimation apparatus according to an embodiment of the present invention. FIG. 同有為差有無判定・有為寿命差見積もり装置の概念構成を示すブロック図である。It is a block diagram which shows the conceptual structure of the same significant difference existence determination / significant lifetime difference estimation apparatus. 同有為差有無判定・有為寿命差見積もり装置を用いた必要寿命差見積もり方法の概略流れ図である。It is a general | schematic flowchart of the required lifetime difference estimation method using the same significant difference existence determination / significant lifetime difference estimation apparatus. 同有為差有無判定・有為寿命差見積もり方法を実施する有為差有無判定・有為寿命差見積もりプログラムの概略の流れ図である。It is a general | schematic flowchart of the significant difference existence determination / significant lifetime difference estimation program which implements the same significant difference presence determination / significant lifetime difference estimation method. 同プログラムにおける有為差有無判定・有為寿命差見積もり手順の詳細を示す流れ図である。It is a flowchart which shows the detail of the significant difference existence determination / significant lifetime difference estimation procedure in the program. 図2の有為差有無判定・有為寿命差見積もり装置における入力画面例の説明図である。It is explanatory drawing of the example of an input screen in the significant difference existence presence determination / significant lifetime difference estimation apparatus of FIG. 図2の有為差有無判定・有為寿命差見積もり装置における出力画面例の説明図である。It is explanatory drawing of the example of an output screen in the significant difference presence / absence determination / significant lifetime difference estimation apparatus of FIG. (A)はワイブル分布の例のグラフ、(B)は頻度と寿命比の関係例を示す確率分布のグラフ、(C)は累積確率と寿命比の関係例を示すグラフ、(D)は少なくとも断定できる寿命差(倍率)と寿命差との関係を示すグラフである。(A) is a graph of an example of the Weibull distribution, (B) is a graph of a probability distribution showing an example of the relationship between the frequency and the life ratio, (C) is a graph showing an example of the relation between the cumulative probability and the life ratio, and (D) is at least It is a graph which shows the relationship between the life difference (magnification) which can be determined, and a life difference. (A)はワイブルスロープ別の寿命倍率と必要寿命差の関係を示すグラフ、(B)は試験個数別の寿命倍率と必要寿命差の関係を示すグラフである。(A) is a graph showing the relationship between the life magnification for each Weibull slope and the required life difference, and (B) is a graph showing the relationship between the life magnification for each test number and the required life difference. ワイブル分布の各パラメータの影響例を示すグラフである。It is a graph which shows the example of influence of each parameter of a Weibull distribution. ワイブル分布の定め方を示すグラフである。It is a graph which shows how to define a Weibull distribution. 実施形態の方法と従来の信頼幅を使った方法の寿命ばらつきの分布比較図である。It is distribution comparison figure of the lifetime dispersion | variation of the method of embodiment, and the method using the conventional reliability width. ワイブル確率紙の説明図である。It is explanatory drawing of a Weibull probability paper. 従来の打切りおよび加速試験の手順を示す流れ図である。It is a flowchart which shows the procedure of the conventional truncation and acceleration test. 従来の信頼幅を使った寿命の有為差判定方法の説明図である。It is explanatory drawing of the lifetime significant difference determination method using the conventional reliability range.

符号の説明Explanation of symbols

1…コンピュータ(演算処理手段)
2…表示装置
3…入力装置
7M…促し画面出力手段
51…有為差有無判定・有為寿命差見積もりプログラム
52…有為差有無判定・有為寿命差見積もり手段
53…第1乱数発生手段
54…第2乱数発生手段
55…寿命比算出手段
56…累積分布作成手段
57…有為差有無判定手段
58…倍率変更繰り返し手段
59…結果出力手段
1 Computer (arithmetic processing means)
2. Display device 3 ... Input device 7M ... Prompt screen output means 51 ... Significant difference presence / absence determination / significant life difference estimation program 52 ... Significant difference presence / absence determination / significant life difference estimation means 53 ... First random number generation means 54 ... second random number generating means 55 ... life ratio calculating means 56 ... cumulative distribution creating means 57 ... significant difference presence / absence judging means 58 ... magnification change repeating means 59 ... result output means

Claims (3)

軸受等の機械部品または試験片からなる試験対象品を使用環境よりも厳しい所定の環境条件におき、破損が発生した時間から寿命を算出する加速試験において、試験対象品の2つのロットの寿命から、有為差有無の判定、および少なくとも断定できる、倍率による寿命差を算出する方法であって、
コンピュータに対し、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命を入力する過程と、
上記コンピュータに、上記有為差有無の判定、および少なくとも断定できる寿命差を演算させ演算結果を表示装置の画面に表示させるコンピュータ演算処理過程とを含み、
上記コンピュータ演算処理過程として、
試験対象品の上記入力情報のワイブルスロープの値および仮設定寿命から定まるワイブル分布に従った乱数を水準1のロットの試験個数分発生させる第1乱数発生手順と、
第1乱数発生手順と同じワイブル分布に従った乱数を水準2のロットの試験個数分発生させる第2乱数発生手順と、
第1乱数発生手順で発生させた設定個数の乱数から所定の寿命計算方法で求まる寿命と第2乱数発生手順で発生させた設定個数の乱数から上記所定の寿命計算方法で求まる寿命とから、1組の水準1,水準2の寿命の寿命比を算出する寿命比算出手順と、
これら第1乱数発生手順,第2乱数発生手順,および寿命比算出手順を設定回数繰り返し、上記設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成して累積確率分布から、設定信頼幅内の最大の寿命比を読み取る累積確率分布分析手順と、 上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手順で得た最大寿命比を超える場合に有為差有りと判定し、最大寿命比以内の場合に有為差無しと判定する有為差有無判定手順と、
この手順で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手順を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す倍率変更繰り返し手順と、
この倍率変更繰り返し手順で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する寿命倍率関係グラフ作成手順と、
上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする有為寿命差読み取り手順と、
上記有為差有無判定手順で判定した結果および上記有為寿命差読み取り手順で読み取った、少なくとも断定できる寿命差を表示装置に表示させる結果出力手順と、
を含む、加速試験における有為差有無判定・有為寿命差見積もり方法。
In an accelerated test that calculates the life from the time at which damage occurred, placing the test object consisting of mechanical parts such as bearings or test pieces under specified environmental conditions that are severer than the operating environment, the life of the two lots of the test object is calculated. , A method for calculating the difference in life according to magnification, which can be determined whether there is a significant difference, and at least determined,
The process of inputting the value of the Weibull slope of the Weibull distribution of the test object, the number of tests in each of the two lots of Level 1 and Level 2 to be compared, and the life of the test results as input information to the computer,
A computer calculation process for causing the computer to determine the presence or absence of the significant difference and to calculate at least the life difference that can be determined and to display the calculation result on the screen of the display device;
As the above computer processing process,
A first random number generation procedure for generating a random number according to the Weibull distribution determined from the value of the Weibull slope of the input information of the test object and the temporarily set life for the number of test pieces of the level 1 lot;
A second random number generation procedure for generating random numbers according to the same Weibull distribution as the first random number generation procedure for the number of test pieces of the level 2 lot;
From the lifetime obtained by the predetermined lifetime calculation method from the set number of random numbers generated in the first random number generation procedure and the lifetime determined by the predetermined lifetime calculation method from the set number of random numbers generated in the second random number generation procedure, 1 A life ratio calculation procedure for calculating the life ratio of the life of the pair of level 1 and level 2,
The first random number generation procedure, the second random number generation procedure, and the life ratio calculation procedure are repeated a set number of times, the life ratio of the set number of times is obtained, and a probability distribution and a cumulative probability distribution of the life ratio are created and accumulated. The cumulative probability distribution analysis procedure that reads the maximum life ratio within the set confidence range from the probability distribution, and the life ratio of the level 1 and level 2 lives in the input information is the maximum life ratio obtained by the cumulative probability distribution analysis procedure. If it exceeds, it is determined that there is a significant difference, and if it is within the maximum life ratio, there is no significant difference.
This procedure is executed when there is a significant difference, and the Weibull distribution is changed to a Weibull distribution with a different set magnification with a temporary set life, and the cumulative probability distribution analysis procedure is set for the cumulative probability distribution analysis procedure. Repeat magnification change procedure that repeats from the minimum magnification to the set maximum magnification,
A life magnification relationship graph creation procedure for creating a life magnification relationship graph, which is a graph showing the relationship between the set magnification repeatedly changed and the maximum life ratio within the reliability range obtained in this magnification change repetition procedure,
Reading the value of the set magnification corresponding to the life ratio of the life of Level 1 and Level 2 in the input information from the life magnification relationship graph and reading the significant life difference as a life difference that can be determined at least. When,
A result output procedure for displaying on the display device a result determined in the significant difference presence / absence determination procedure and a life difference that can be determined at least read in the significant life difference reading procedure;
Including a significant difference in acceleration tests, and a method for estimating a significant life difference.
軸受等の機械部品または試験片からなる試験対象品を使用環境よりも厳しい所定の環境条件におき、破損が発生した時間から寿命を算出する加速試験において、試験対象品の2つのロットの寿命から、有為差有無の判定、および少なくとも断定できる、倍率による寿命差を算出する装置であって、
演算処理装置と、この演算処理装置の出力を画面に表示する表示装置と、上記演算処理装置に入力を行う入力手段とを備え、
上記演算処理装置は、
上記表示装置の画面に、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命の入力を促す表示を行わせる促し画面出力手段と、
実行命令に応答して上記有為差有無の判定および少なくとも断定できる寿命差の演算を行いその演算結果を上記表示装置の画面に出力する有為差判定・有為寿命差演算手段とを備え、
上記有為差判定・有為寿命差演算手段は、
上記入力情報のワイブルスロープの値および仮設定寿命から定まるワイブル分布に従った乱数を水準1のロットの試験個数分発生させる第1乱数発生手段と、
第1乱数発生手段と同じワイブル分布に従った乱数を水準2のロットの試験個数分発生させる第2乱数発生手段と、
第1乱数発生手段で発生させた設定個数の乱数から所定の寿命計算方法で求まる寿命と第2乱数発生手段で発生させた設定個数の乱数から上記所定の寿命計算方法で求まる寿命とから、1組の水準1,水準2の寿命の寿命比を算出する寿命比算出手段と、
これら第1乱数発生手段,第2乱数発生手段,および寿命比算出手段の処理を設定回数繰り返し、上記設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成して累積確率分布から、設定信頼幅内の最大の寿命比を読み取る累積確率分布分析手段と、
上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手段で得た最大寿命比を超える場合に有為差有りと判定し、最大寿命比以内の場合に有為差無しと判定する有為差有無判定手段と、
この手段で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手段の処理を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す倍率変更繰り返し手段と、
この倍率変更繰り返し手段で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する寿命倍率関係グラフ作成手段と、
上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする有為寿命差読み取り手段と、
上記有為差有無判定手段で判定した結果および上記有為寿命差読み取り手段で読み取った、少なくとも断定できる寿命差を表示装置に表示させる結果出力手段と、
を含む、加速試験における有為差有無判定・有為寿命差見積もり装置。
In an accelerated test that calculates the life from the time at which damage occurred, placing the test object consisting of mechanical parts such as bearings or test pieces under specified environmental conditions that are severer than the operating environment, the life of the two lots of the test object is calculated. , A device for calculating the difference in life by magnification, which can be determined whether there is a significant difference, and at least determined,
An arithmetic processing device, a display device for displaying the output of the arithmetic processing device on a screen, and input means for inputting to the arithmetic processing device,
The arithmetic processing unit is
On the screen of the display device, as input information, a display prompting the user to input the value of the Weibull slope of the Weibull distribution of the product under test, the number of each test in the two lots of Level 1 and Level 2 to be compared, and the life of the test results Prompt screen output means to perform,
In response to the execution command, comprises the presence / absence of the significant difference and the calculation of the life difference that can be determined at least and outputs the calculation result to the screen of the display device.
The above significant difference judgment and significant life difference calculation means are:
First random number generating means for generating random numbers according to the Weibull distribution determined from the value of the Weibull slope of the input information and the temporarily set life for the number of test pieces of the level 1 lot;
Second random number generating means for generating random numbers according to the same Weibull distribution as the first random number generating means for the test number of lots of level 2;
From the lifetime determined by the predetermined lifetime calculation method from the set number of random numbers generated by the first random number generator and the lifetime determined by the predetermined lifetime calculation method from the set number of random numbers generated by the second random number generator. A life ratio calculating means for calculating the life ratio of the life of the pair of level 1 and level 2;
The processing of the first random number generation means, the second random number generation means, and the life ratio calculation means is repeated a set number of times, the life ratio of the set number of times is obtained, and a probability distribution and a cumulative probability distribution of the life ratio are created. A cumulative probability distribution analysis means for reading the maximum life ratio within the set reliability range from the cumulative probability distribution;
It is judged that there is a significant difference when the life ratio of the life of Level 1 and Level 2 in the input information exceeds the maximum life ratio obtained by the cumulative probability distribution analysis means, and there is no significant difference when it is within the maximum life ratio. Means for determining whether there is a significant difference, and
This means is executed when there is a significant difference, and the Weibull distribution is sequentially changed to a Weibull distribution having a different set magnification with a temporarily set life, and the process of the cumulative probability distribution analysis means is performed with the set magnification. Magnification change repetition means that repeats until the set minimum magnification reaches the set maximum magnification,
A life magnification relationship graph creating means for creating a life magnification relationship graph, which is a graph showing the relationship between the set magnification repeatedly changed and the maximum life ratio within the reliability range obtained by the magnification change repetition means,
Meaning life difference reading means for reading the value of the set magnification corresponding to the life ratio of the life of level 1 and level 2 in the input information from the life magnification relation graph and making the read value at least a life difference that can be determined When,
A result output means for displaying a result determined by the significant difference presence / absence determining means and a life difference that can be determined at least, read by the significant life difference reading means, on a display device;
Including significant difference in acceleration test, and useful life difference estimation device.
コンピュータで実行可能なプログラムであって、
表示装置の画面に、入力情報として、試験対象品の持つワイブル分布のワイブルスロープの値、比較する水準1,水準2の2つのロットの各試験個数および試験結果の寿命の入力を促す表示を行わせる促し画面出力手順と、
実行命令に応答して上記有為差有無の判定および少なくとも断定できる寿命差の演算を行いその演算結果を上記表示装置の画面に出力する有為差判定・有為寿命差演算手順とを含み、
この有為差判定・有為寿命差演算手順は、
上記入力情報のワイブルスロープの値および仮設定寿命から定まるワイブル分布に従った乱数を水準1のロットの試験個数分発生させる第1乱数発生手順と、
第1乱数発生手順と同じワイブル分布に従った乱数を水準2のロットの試験個数分発生させる第2乱数発生手順と、
第1乱数発生手順で発生させた設定個数の乱数から所定の寿命計算方法で求まる寿命と第2乱数発生手順で発生させた設定個数の乱数から上記所定の寿命計算方法で求まる寿命とから、1組の水準1,水準2の寿命の寿命比を算出する寿命比算出手順と、
これら第1乱数発生手順,第2乱数発生手順,および寿命比算出手順を設定回数繰り返し、上記設定回数の組数の寿命比を求め、この寿命比の確率分布と累積確率分布を作成して累積確率分布から、設定信頼幅内の最大の寿命比を読み取る累積確率分布分析手順と、 上記入力情報における水準1,水準2の寿命の寿命比が上記累積確率分布分析手順で得た最大寿命比を超える場合に有為差有りと判定し、最大寿命比以内の場合に有為差無しと判定する有為差有無判定手順と、
この手順で有為差有りとした場合に実行され、上記ワイブル分布を仮設定寿命が異なる設定倍率となるワイブル分布に繰り返し毎に順次変えて、上記累積確率分布分析手順を、上記設定倍率が設定最小倍率から設定最大倍率となるまで繰り返す倍率変更繰り返し手順と、
この倍率変更繰り返し手順で得られた、繰り返し変更した設定倍率と上記信頼幅内の最大の寿命比の関係を示すグラフである寿命倍率関係グラフを作成する寿命倍率関係グラフ作成手順と、
上記寿命倍率関係グラフから、上記入力情報における水準1,水準2の寿命の寿命比に対応する上記設定倍率の値を読み取り、その読み取った値を少なくとも断定できる寿命差とする有為寿命差読み取り手順と、
上記有為差有無判定手順で判定した結果および上記有為寿命差読み取り手順で読み取った、少なくとも断定できる寿命差を表示装置に表示させる結果出力手順と、
を含む、加速試験における有為差有無判定・有為寿命差見積もりプログラム。
A program executable on a computer,
On the screen of the display device, as input information, a display that prompts the user to input the Weibull slope value of the Weibull distribution of the test object, the number of each test in the two lots of Level 1 and Level 2 to be compared, and the life of the test results Prompt screen output procedure,
In response to the execution command, includes the determination of the presence or absence of a significant difference and the calculation of at least the life difference that can be determined, and the calculation result is output to the screen of the display device.
This significant difference judgment and significant life difference calculation procedure is:
A first random number generation procedure for generating a random number according to the Weibull distribution determined from the value of the Weibull slope of the input information and the temporarily set life for the number of test pieces of the level 1 lot;
A second random number generation procedure for generating random numbers according to the same Weibull distribution as the first random number generation procedure for the number of test pieces of the level 2 lot;
From the lifetime obtained by the predetermined lifetime calculation method from the set number of random numbers generated in the first random number generation procedure and the lifetime determined by the predetermined lifetime calculation method from the set number of random numbers generated in the second random number generation procedure, 1 A life ratio calculation procedure for calculating the life ratio of the life of the pair of level 1 and level 2,
The first random number generation procedure, the second random number generation procedure, and the life ratio calculation procedure are repeated a set number of times, the life ratio of the set number of times is obtained, and a probability distribution and a cumulative probability distribution of the life ratio are created and accumulated. The cumulative probability distribution analysis procedure that reads the maximum life ratio within the set confidence range from the probability distribution, and the life ratio of the level 1 and level 2 lives in the input information is the maximum life ratio obtained by the cumulative probability distribution analysis procedure. If it exceeds, it is determined that there is a significant difference, and if it is within the maximum life ratio, there is no significant difference.
This procedure is executed when there is a significant difference, and the Weibull distribution is changed to a Weibull distribution with a different set magnification with a temporary set life, and the cumulative probability distribution analysis procedure is set for the cumulative probability distribution analysis procedure. Repeat magnification change procedure that repeats from the minimum magnification to the set maximum magnification,
A life magnification relationship graph creation procedure for creating a life magnification relationship graph, which is a graph showing the relationship between the set magnification repeatedly changed and the maximum life ratio within the reliability range obtained in this magnification change repetition procedure,
Reading the value of the set magnification corresponding to the life ratio of the life of Level 1 and Level 2 in the input information from the life magnification relationship graph and reading the significant life difference as a life difference that can be determined at least. When,
A result output procedure for displaying on the display device a result determined in the significant difference presence / absence determination procedure and a life difference that can be determined at least read in the significant life difference reading procedure;
A program for determining the presence / absence of a significant difference and estimating a significant life difference in an accelerated test.
JP2006311138A 2006-11-17 2006-11-17 Method and apparatus for estimating necessary life difference of significant difference determination in accelerated test Pending JP2008128698A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2006311138A JP2008128698A (en) 2006-11-17 2006-11-17 Method and apparatus for estimating necessary life difference of significant difference determination in accelerated test

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2006311138A JP2008128698A (en) 2006-11-17 2006-11-17 Method and apparatus for estimating necessary life difference of significant difference determination in accelerated test

Publications (1)

Publication Number Publication Date
JP2008128698A true JP2008128698A (en) 2008-06-05

Family

ID=39554687

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006311138A Pending JP2008128698A (en) 2006-11-17 2006-11-17 Method and apparatus for estimating necessary life difference of significant difference determination in accelerated test

Country Status (1)

Country Link
JP (1) JP2008128698A (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009115676A (en) * 2007-11-08 2009-05-28 Ntn Corp Weibull slope estimation method and device in life test
CN102252898A (en) * 2011-03-09 2011-11-23 北京航空航天大学 Method for testing accelerated life of electronic product based on life-stress model
JP2014067098A (en) * 2012-09-24 2014-04-17 Yahoo Japan Corp Work evaluation index generation device
CN112557037A (en) * 2020-11-30 2021-03-26 南京工业大学 Accelerated life test method for wind power main shaft bearing model
CN113762981A (en) * 2021-03-30 2021-12-07 中国人民解放军国防科技大学 Product credibility calculation method based on exponential distribution

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009115676A (en) * 2007-11-08 2009-05-28 Ntn Corp Weibull slope estimation method and device in life test
CN102252898A (en) * 2011-03-09 2011-11-23 北京航空航天大学 Method for testing accelerated life of electronic product based on life-stress model
JP2014067098A (en) * 2012-09-24 2014-04-17 Yahoo Japan Corp Work evaluation index generation device
CN112557037A (en) * 2020-11-30 2021-03-26 南京工业大学 Accelerated life test method for wind power main shaft bearing model
CN112557037B (en) * 2020-11-30 2022-03-11 南京工业大学 Accelerated life test method for wind power main shaft bearing model
CN113762981A (en) * 2021-03-30 2021-12-07 中国人民解放军国防科技大学 Product credibility calculation method based on exponential distribution
CN113762981B (en) * 2021-03-30 2023-11-28 中国人民解放军国防科技大学 Product credibility calculation method based on index distribution

Similar Documents

Publication Publication Date Title
Williams et al. A practical method for statistical analysis of strain–life fatigue data
JP5669630B2 (en) Test case generation method, program and system
JP2008128699A (en) Weibull slope estimation method and device of lifetime test
Agius et al. Sensitivity and optimisation of the Chaboche plasticity model parameters in strain-life fatigue predictions
US7657408B2 (en) Structural analysis apparatus, structural analysis method, and structural analysis program
JP4859634B2 (en) Life test design / judgment method
JP2008128698A (en) Method and apparatus for estimating necessary life difference of significant difference determination in accelerated test
WO2022149297A1 (en) Fatigue life prediction method, fatigue life prediction device, fatigue life prediction program, and recording medium
JP2007108843A (en) Semiconductor device design support method, semiconductor device design support system and semiconductor device design support program
WO2014033927A1 (en) Fissure progress estimation method and information processing device
JP2008145414A (en) Life estimation method and system from truncated life test
JP2017090071A (en) Hardness tester and hardness testing method
JP2008128690A (en) Method and apparatus for estimating ending time in life test and test stop reference
JP2008128696A (en) Method and apparatus for estimating necessary number of tests in accelerated test
JP2008128697A (en) Method and apparatus for estimating necessary life difference of significant difference determination in accelerated test
Guzmán et al. Damage characterization in a ferritic steel sheet: Experimental tests, parameter identification and numerical modeling
Pham et al. Application of a graphical method on estimating forming limit curve of automotive sheet metals
JP5895358B2 (en) Destructive inspection system, destructive inspection method, data processing program, and program recording medium
JP6482779B2 (en) Life test interpretation and design methods, equipment, and programs
JP2007232698A (en) Young&#39;s modulus estimation method, young&#39;s modulus estimation program and young&#39;s modulus estimation device
Oliveira Filho et al. Analysis of sensor placement in beams for crack identification
CN111124918B (en) Test data prediction method and device and processing equipment
JPWO2015037067A1 (en) Ductile fracture evaluation system and ductile fracture evaluation method
Ahmad et al. The probabilistic analysis of fatigue crack effect based on magnetic flux leakage
JP2017102055A (en) Methods for evaluating and calculating compaction characteristics of reinforced concrete, calculation device, and method for designing reinforced concrete