MAE 430 Project I
1
Project I
11.27.2014
Group 8 Saeb Moosavi, Shakil Bin Zaman MAE 430 Introduction To Reliability In Mechanical Engineering Design Professor Ji Ho Song
Project I 11.27.2014 Group 8 Saeb Moosavi, Shakil Bin Zaman MAE - - PowerPoint PPT Presentation
Project I 11.27.2014 Group 8 Saeb Moosavi, Shakil Bin Zaman MAE 430 Introduction To Reliability In Mechanical Engineering Design Professor Ji Ho Song MAE 430 Project I 1 Contents Data set 1 (Saeb) Linearity Test K-S test
MAE 430 Project I
1
11.27.2014
Group 8 Saeb Moosavi, Shakil Bin Zaman MAE 430 Introduction To Reliability In Mechanical Engineering Design Professor Ji Ho Song
MAE 430 Project I
2
MAE 430 Project I
3
600 38 128 270 376 66 120 47 134 382 187 656 266 661 337 57 681 384 565 53 632 373 141 7 378 253
MAE 430 Project I
4
Step1 Cumulative distribution function estimation Step2 Probability estimation using probability paper Step3 Goodness of fit test (Kolomogorov-Smirnov test)
600 38 128 270 376 66 120 47 134 382 187 656 266 661 337 57 681 384 565 53 632 373 141 7 378 253
26 Data
Target: select the best distribution of given data
MAE 430 Project I
5
F(xj) = j / n
F(xj) = ( j - 0.5) / n
F(xj) = j /(n +1)
F(xj) = ( j - 0.3) /(n + 0.4)
F(xj) = ( j -1) /(n -1)
F(xj) = ( j - 0.375) /(n + 0.25)
Simple C.D. & Mode rank can’t be use because F(x) = 0 or 1
C.D.F. estimation method
MAE 430 Project I
6
Probability distribution
MAE 430 Project I
7
Symmetric simple cumulative distribution
100 200 300 400 500 600 700
1 2
Norm A
Equation y = a + b*
0.91477 Value Standard Err Norm Intercept
0.09701 Norm Slope 0.0043 2.61958E-4
Normal LogNormal Weibull Biexponential
2 3 4 5 6 7
1 2
C A
Equation y = a + b*
0.9666 Value Standard Erro C Intercept
0.21831 C Slope 1.0894 0.04047
2 3 4 5 6 7
1 2
B A
Equation y = a + b*
0.88426 Value Standard Erro B Intercept
0.32344 B Slope 0.83085 0.05996
100 200 300 400 500 600 700
1 2
C A
Equation y = a + b
0.8072 Value Standard Er C Intercept
0.18333 C Slope 0.0050 4.95067E-4
MAE 430 Project I
8
Normal LogNormal Weibull Biexponential
100 200 300 400 500 600 700
0.0 0.5 1.0 1.5 2.0
C A
Equation y = a + b*
0.92916 Value Standard Err C Intercept
0.08076 C Slope 0.0039 2.18083E-4 100 200 300 400 500 600 700
1 2
C A
Equation y = a + b
0.84288 Value Standard Er C Intercept -1.923 0.14794 C Slope 0.0046 3.9949E-4
Mean Rank
2 3 4 5 6 7
0.0 0.5 1.0 1.5 2.0
C A
Equation y = a + b
0.88285 Value Standard Err C Intercept
0.29714 C Slope 0.7581 0.05509
2 3 4 5 6 7
1 2
C A
Equation y = a + b
0.95881 Value Standard Err C Intercept
0.21671 C Slope 0.97004 0.04018
MAE 430 Project I
9
Normal LogNormal Weibull Biexponential
100 200 300 400 500 600 700
0.0 0.5 1.0 1.5 2.0
C A
Equation y = a + b*
0.92212 Value Standard Err C Intercept
0.08907 C Slope 0.0041 2.40535E-4
Median Rank
100 200 300 400 500 600 700
1 2
C A
Equation y = a + b*x
0.84288 Value Standard Error C Intercept
0.14794 C Slope 0.00464 3.9949E-4
2 3 4 5 6 7
0.0 0.5 1.0 1.5 2.0
C A
Equation y = a + b
0.88386 Value Standard Er C Intercept
0.31123 C Slope 0.7979 0.0577 2 3 4 5 6 7
1 2
C A
Equation y = a +
0.9635 Value Standard Er C Intercep
0.21666 C Slope 1.0339 0.04017
MAE 430 Project I
10
Normal LogNormal Weibull Biexponential The Rest Method
100 200 300 400 500 600 700
1 2
C A
Equation y = a + b*
0.91972 Value Standard Err C Intercept
0.09174 C Slope 0.0042 2.4773E-4
100 200 300 400 500 600 700
1 2
C A
Equation y = a + b*
0.81911 Value Standard Err C Intercept
0.17181 C Slope 0.0049 4.63975E-4 2 3 4 5 6 7
1 2
C A
Equation y = a +
0.96477 Value Standard Er C Intercep
0.21693 C Slope 1.0530 0.04022 2 3 4 5 6 7
1 2
C A
Equation y = a +
Value Standard Er C Intercep
0.31543 C Slope 0.8094 0.05848
MAE 430 Project I
11
Symm.S.C Mean Rank Median Rank Rest R2 R R2 R R2 R R2 R Normal 0.9181 0.95822 0.9319 0.9654 0.9252 0.96189 0.9229 0.96069 Log-Normal 0.8888 0.94281 0.7667 0.87565 0.8884 0.9426 0.974 0.98694 Weibull 0.9679 0.98384 0.756 0.8695 0.8954 0.9463 0.787 0.88717 Bi-exponential 0.8149 0.90272 0.8491 0.9215 0.8491 0.9215 0.8263 0.90903
MAE 430 Project I
12
α = 0.198
Normal & LogNormal
α = 0.204
Weibull & Extreme Value
α = 0.170
Normal & LogNormal
α = 0.175
Weibull & Extreme Value
MAE 430 Project I
13
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Biexponential a = 0.01, D = 0.204 a = 0.05, D = 0.175
B A
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Weibull a = 0.01, D = 0.204 a = 0.05, D = 0.175
B A
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
LogNormal a = 0.01, D = 0.198 a = 0.05, D = 0.170
B A
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Normal a = 0.01, D = 0.198 a = 0.05, D = 0.170
B A
σ = 232.56 μ = 299.53 σ = 1.203 μ = 5.279 m = 1.0894 ξ = 329.96 ξ = 200 X0 = 418.2 Symmetric simple cumulative distribution
MAE 430 Project I
14
Mean Rank σ = 256.41 μ = 303.92 σ = 1.319 μ = 5.279 m = 0.97 ξ = 339.54 ξ = 217.39 X0 = 418.04
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Normal a = 0.01, D = 0.198 a = 0.05, D = 0.170
B A
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
LogNormal a = 0.01, D = 0.198 a = 0.05, D = 0.170
B A
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Weibull a = 0.01, D = 0.204 a = 0.05, D = 0.175
B A
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Biexponential a = 0.01, D = 0.204 a = 0.05, D = 0.175
B A
MAE 430 Project I
15
Median Rank σ = 243.9 μ = 303 σ = 1.253 μ = 5.279 m = 1.0339 ξ = 333.95 ξ = 215.52 X0 = 414.59
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Normal a = 0.01, D = 0.198 a = 0.05, D = 0.170
B A
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
LogNormal a = 0.01, D = 0.198 a = 0.05, D = 0.170
B A
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Weibull a = 0.01, D = 0.204 a = 0.05, D = 0.175
B A
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Biexponential a = 0.01, D = 0.204 a = 0.05, D = 0.175
B A
MAE 430 Project I
16
The Rest Method σ = 238.09 μ = 299.66 σ = 1.235 μ = 5.279 m = 1.053 ξ = 332.69 ξ = 204.08 X0 = 416.53
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Normal a = 0.01, D = 0.198 a = 0.05, D = 0.170
B A
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
LogNormal a = 0.01, D = 0.198 a = 0.05, D = 0.170
B A
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Weibull a = 0.01, D = 0.204 a = 0.05, D = 0.175
B A
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Biexponential a = 0.01, D = 0.204 a = 0.05, D = 0.175
B A
MAE 430 Project I
17
According to K-S test results the best distribution function for this data set is Weibull distribution because data are more concentrated in the center of limit line in this distribution.
MAE 430 Project I
18
262 500 328 91 608 119 15 164 211 89 668 278 319 116 330 449 74 128 622 223 98
MAE 430 Project I
19
100 200 300 400 500 600 700
0.0 0.5 1.0 1.5 2.0
C Linear Fit of C C A
Equation y = a + b*x
0.85896 Value Standard Error C Intercept
0.14318 C Slope 0.00452 4.41996E-4
Normal
0.0 0.2 0.4 0.6 0.8 1.0 1.2
0.0 0.5 1.0 1.5 2.0
D Linear Fit of D D A
Equation y = a + b*xLogNormal 0.0 0.2 0.4 0.6 0.8 1.0 1.2
1
E Linear Fit of E E A
Equation y = a + b*xWeibull
100 200 300 400 500 600 700
1 2
K Linear Fit of K K A
Equation y = a + b*xBiexponential
MAE 430 Project I
20
100 200 300 400 500 600 700
0.0 0.5 1.0 1.5 2.0
E Linear Fit of E E A
Equation y = a + b*x
0.86499 Value Standard Erro E Intercept
0.12459 E Slope 0.00403 3.8459E-4
Normal
0.0 0.2 0.4 0.6 0.8 1.0 1.2
0.0 0.5 1.0 1.5 2.0
G Linear Fit of G G A
Equation y = a + b*xLogNormal
0.0 0.2 0.4 0.6 0.8 1.0 1.2
1
H Linear Fit of H H A
Equation y = a + b*xWeibull
100 200 300 400 500 600 700
1
H Linear Fit of H H A
Equation y = a + b*xBiexponential
MAE 430 Project I
21
100 200 300 400 500 600 700
0.0 0.5 1.0 1.5 2.0
L Linear Fit of L L A
Equation y = a + b*xNormal
0.0 0.2 0.4 0.6 0.8 1.0 1.2
0.0 0.5 1.0 1.5 2.0
J Linear Fit of J J A
Equation y = a + b*xLogNormal
0.0 0.2 0.4 0.6 0.8 1.0 1.2
1 2
K Linear Fit of K K A
Equation y = a + b*xWeibull
100 200 300 400 500 600 700
1 2
O Linear Fit of O O A
Equation y = a + b*xBiexponential
MAE 430 Project I
22
100 200 300 400 500 600 700
0.0 0.5 1.0 1.5 2.0
S Linear Fit of S S A
Equation y = a + b*Normal
0.0 0.2 0.4 0.6 0.8 1.0 1.2
0.0 0.5 1.0 1.5 2.0
M Linear Fit of M M A
Equation y = a + bLogNormal
0.0 0.2 0.4 0.6 0.8 1.0 1.2
1 2
N Linear Fit of N N A
Equation y = a + b*xWeibull
100 200 300 400 500 600 700
1 2
V Linear Fit of V V A
Equation y = a + b*xBiexponential
MAE 430 Project I
23
R SD Bie 0.85355 0.6664 Log 0.91496 0.36748 Nor 0.93126 0.37306 Wei 0.94773 0.35471
Median
R SD Bie 0.86261 0.60716 Log 0.91847 0.3844 Nor 0.9330 0.34977 Wei 0.95487 0.35649
The rest
R SD Bie 0.85963 0.62703 Log 0.91954 0.38908 Nor 0.93247 0.35765 Wei 0.95706 0.35578
Mean
R SD Bie 0.87155 0.54506 Log 0.91496 0.36748 Nor 0.93431 0.32461 Wei 0.94773 0.35471
MAE 430 Project I
24
Symm. S.C. Mean Median The rest Normal O O O O LogNormal O O O O Weibull O O O O Bi-exponential X X X X
MAE 430 Project I
25
n
D
= 0.207
n
D
= 0.202
n
D
= 0.178
n
D
= 0.175
MAE 430 Project I
26
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
B B A Normal
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
B B A Weibull
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
B B A Biexponential
Red : alpha = 0.05 Blue : alpha = 0.15
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
B A LogNormal
µ= 5.29 σ=1.886 µ= 255.61 σ=312.88 m= 0.8369 ξ= 202.48 x0= 364.5 ξ= 193.42
MAE 430 Project I
27
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
B B A Normal
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
B B A Weibull
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
B B A Biexponential
Red : alpha = 0.05 Blue : alpha = 0.15
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
B A LogNormal
µ= 5.29 σ=1.886 x0= 364.3 ξ= 198.42 µ= 255.78 σ=322.88 m= 0.8469 ξ= 212.48
MAE 430 Project I
28
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
D D A Normal
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
D D A Weibull
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
B B A Biexponential
Red : alpha = 0.05 Blue : alpha = 0.15
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
B A LogNormal
µ= 5.29 σ=1.886 x0= 362.5 ξ= 191.52 µ= 252.61 σ=314.22 m= 0.8769 ξ= 232.48
MAE 430 Project I
29
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
F F A Normal
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
F F A Weibull
Red : alpha = 0.05 Blue : alpha = 0.15
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
B B A Biexponential
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
B A LogNormal
µ= 5.29 σ=1.886 x0= 364.5 ξ= 193.42 µ= 254.21 σ=311.86 m= 0.8369 ξ= 202.48
MAE 430 Project I
30
– Weibull distribution makes the best distribution of the given data followed by Normal distribution (according to R-value and Standard deviation).(O) – It is imperative that the Bi-exponential distribution is skewed from the straight line. (X)
– Normal distribution gives a better distribution (all values inside the significance level) with no outliers.(O) – LogNormal distribution is rejected undoubtedly. – Weibull and Bi-exponential distributions show a scattered distribution as far as the significance value is concerned. (X) – Both distributions contain some outliers outside the 5% significance level.
MAE 430 Project I
31
262 500 328 91 608 119 15 164 211 89 668 278 319 116 330 449 74 128 622 223 98
600 38 128 270 376 66 120 47 134 382 187 656 266 661 337 57 681 384 565 53 632 373 141 7 378 253
MAE 430 Project I
32
Symmetrical Simple Cumulative Distribution
100 200 300 400 500 600 700
1 2 3
ND Linear Fit of ND ND A
Equation y = a + b*x
0.91505 Value Standard Error ND Intercept
0.07248 ND Slope 0.00457 2.04908E-4
Normal
2 3 4 5 6 7
1 2 3
LN Linear Fit of C LN A
Equation y = a + b*xLogNormal
2 3 4 5 6 7
1 2
W Linear Fit of D W A
Equation y = a + b*x
0.98428 Value Standard Erro D Intercept
0.12215 D Slope 1.2187 0.0227
Weibull
100 200 300 400 500 600 700
1 2
BI Linear Fit of BI BI A
Equation y = a + b*xBi-exponential
MAE 430 Project I
33
100 200 300 400 500 600 700
1 2
E Linear Fit of E E A
Equation y = a + b*xNormal
2 3 4 5 6 7
1 2
F Linear Fit of F F A
Equation y = a + b*xLogNormal
2 3 4 5 6 7
1 2
G Linear Fit of G G A
Equation y = a + b*xWeibull
100 200 300 400 500 600 700
1 2
F Linear Fit of F F A
Equation y = a + b*xBi-exponential
MAE 430 Project I
34
100 200 300 400 500 600 700
1 2 3
D Linear Fit of D D A
Equation y = a + b*xNormal
2 3 4 5 6 7
1 2 3
I Linear Fit of I I A
Equation y = a + b*xLogNormal
2 3 4 5 6 7
1 2
C Linear Fit of C C A
Equation y = a + b*x
0.98304 Value Standard Error C Intercept
0.12254 C Slope 1.1763 0.02278
Weibull
100 200 300 400 500 600 700
1 2
J Linear Fit of J J A
Equation y = a + b*xBi-exponential
MAE 430 Project I
35
100 200 300 400 500 600 700
1 2 3
M Linear Fit of M M A
Equation y = a + b*xNormal
2 3 4 5 6 7
1 2 3
J Linear Fit of J J A
Equation y = a + b*xLog Normal
2 3 4 5 6 7
1 2
K Linear Fit of K K A
Equation y = a + b*xWeibull
100 200 300 400 500 600 700
1 2
N Linear Fit of N N A
Equation y = a + b*xBi-exponential
MAE 430 Project I
36
R SD Bie 0.8904 0.58066 Log 0.9551 0.2961 Nor 0.9151 0.2961 Wei 0.9928 0.1582
Median
R SD Bie 0.8984 0.2737 Log 0.9557 0.2890 Nor 0.9604 0.2737 Wei 0.9917 0.1587
The rest
R SD Bie 0.8957 0.5546 Log 0.9558 0.2914 Nor 0.9595 0.2794 Wei 0.9920 0.1579
Mean
R SD Bie 0.9065 0.4986 Log 0.9553 0.2810 Nor 0.9632 0.2554 Wei 0.9903 0.1643
MAE 430 Project I
37
C. Mean Median The rest Normal O O O O LogNormal X X O O Weibull O O O O Biexponential X X X X
MAE 430 Project I
38
α=0.01
α=0.05
n
D
= 0.154
n
D
= 0.143
n
D
= 0.149
n
D
= 0.128
MAE 430 Project I
39
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
SSC SSC A Bi-exponential
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
SSC SSC A Weibull
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
SSC SSC A LogNormal
Red : alpha = 0.05 Blue : alpha = 0.01
µ= 5.29 σ=1.886 µ= 218.8 σ=286.7 m= 0.9369 ξ= 302.48 x0= 364.5 ξ= 193.42
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
SSC SSC A Normal
MAE 430 Project I
40
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Mean Mean A Bi-exponential
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Mean Mean A Weibull
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Mean Mean A LogNormal
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Mean Mean A Normal
Red : alpha = 0.05 Blue : alpha = 0.01
µ= 5.31 σ=1.686 x0= 354.5 ξ= 192.62 µ= 276.7 σ=316.9 m= 0.8969 ξ= 312.48
MAE 430 Project I
41
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Median Median A Weibull
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Median Median A LogNormal
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Median Median A Normal
Red : alpha = 0.05 Blue : alpha = 0.01
µ= 5.32 σ=1.686 x0= 364.5 ξ= 193.42 µ= 276.7 σ=216.9 m= 0.8969 ξ= 312.48
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Median Median A Bi-exponential
MAE 430 Project I
42
Red : alpha = 0.05 Blue : alpha = 0.01
µ= 5.31 σ=1.686 x0= 364.5 ξ= 193.42 µ= 276.7 σ=216.9 m= 0.8969 ξ= 312.48
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Rest Rest A Normal
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Median Median A LogNormal
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Rest Rest A Weibull
100 200 300 400 500 600 700 0.0 0.2 0.4 0.6 0.8 1.0
Rest Rest A Bi-exponential
MAE 430 Project I
43
– Bi-exponential and Lognormal distributions are not linear in comparison to that of Weibull and Normal. – Weibull shows a better distribution of the linearity test with smallest R-value and SD.
– Weibull distribution has the best fit distribution with all points inside the 5% significance range. – LogNormal and Bi-exponential distributions have values outside the 5% significance level and is therefore rejected. (X) – Normal distribution shows the poorest distribution of K-S test and is therefore rejected.(X) – Normal distribution has points outside the significance level of 5% which is unacceptable.
MAE 430 Project I
44