Distribusi Probabilitas … (Suyono) 165

DISTRIBUSI PROBABILITAS DARI TOTAL BIAYA PERAWATAN KOMPONEN

Suyono1

ABSTRACT

Article discusses the probability distribution of the total maintenance cost of a component over a finite time horizon. The maintenance cost is assumed to be a function of the component lifetime. The marginal distribution of the total maintenance cost is presented in the form of Laplace transform. Article also presents the mean and second moment of the total maintenance cost. Keywords: Laplace transform, renewal process, total maintenance cost.

ABSTRAK

Artikel membahas distribusi probabilitas dari total biaya perawatan komponen pada interval waktu berhingga. Biaya perawatan dianggap sebagai fungsi dari waktu hidup komponen. Distribusi probabilitas dari total biaya perawatan komponen, beserta mean dan momen keduanya disajikan dalam bentuk transformasi Laplace. Kata Kunci: transformasi Laplace, proses renewal, biaya maintenance total

1 Jurusan Matematika, FMIPA, Universitas Negeri Jakarta, Jl.Pemuda Rawamangun, Jakarta

synjkt@yahoo.com

Jurnal Mat Stat, Vol. 8 No. 2 Juli 2008: 165-170 166

PENDAHULUAN

Perhatikan sebuah komponen yang merupakan elemen dari sebuah sistem. Agar sistem dapat bekerja dengan baik maka komponen tersebut harus dirawat. Perawatan dapat dilakukan secara preventif, yakni perawatan yang dilakukan sebelum komponen tersebut rusak, atau dilakukan penggantian ketika komponen sudah rusak. Dalam makalah ini dianggap kedua macam perawatan ini akan membuat komponen dapat bekerja kembali sebagaimana komponen yang baru. Terkait dengan perawatan komponen, salah satu kuantitas yang penting untuk dianalisa adalah total biaya perawatan (life-cycle cost) pada suatu interval waktu, lihat Wagner (1975).

Beberapa peneliti telah membahas total biaya perawatan komponen. Harga harapan rata-rata biaya per satuan waktu dalam interval waktu [0,∞) telah dibahas oleh Noortwijk dan Frangopol (2004). Dalam makalah yang lain Noortwijk (2003) membahas sifat asimtotik dari variansi total biaya perawatan. Dalam makalah ini akan dibahas total biaya perawatan pada interval waktu berhingga [0,t], t ≥ 0. Susunan makalah ini adalah sebagai berikut. Di Bagian 2 dibahas model matematika untuk total biaya perawatan komponen. Selanjutnya di Bagian 3 dibahas distribusi probabilitas total biaya perawatan komponen yang disajikan dalam bentuk transformasi Laplace. Mean dan momen kedua total biaya perawatan disajikan pada Bagian 4. Bagian terakhir berisi sebuah contoh aplikasi.

MODEL TOTAL BIAYA PERAWATAN KOMPONEN

Misalkan 0 = S0 < S1 < S2 < … adalah waktu-waktu dilakukannya perawatan (baik perawatan preventif maupun penggantian) suatu komponen dari sebuah sistem. Perawatan dilakukan sedemikian hingga akan menjadikan komponen dapat bekerja kembali sebagaimana komponen yang baru. Dengan asumsi ini barisan waktu-waktu perawatan dari komponen, yakni (Sj, j ≥ 1), dapat dimodelkan dengan proses renewal. Ini berarti bahwa waktu-waktu antar perawatan, yakni

Tj = Sj – Sj-1, j = 1, 2, 3, …,

merupakan variabel-variabel acak non-negatif, saling independen dan berdistribusi identik. Banyaknya perbaikan komponen yang dilakukan pada interval waktu [0,t], t ≥ 0, adalah

N(t) = max{j ≥ 0 : Sj ≤ t}.

Notasikan dengan Cj biaya perawatan yang dilakukan pada waktu Sj. Maka total biaya perawatan komponen pada interval waktu [0,t] dapat dinyatakan sebagai

⎪⎩

⎪⎨⎧

>≤= ∑ =

tTtTCtK

tN

j j

1

)(

1 1

jika,0jika,)( (1)

Selanjutnya akan diasumsikan ((Tj, Cj), j ≥ 0) merupakan barisan vektor acak non-negatif, saling independen dan berdistribusi identik dengan fungsi distribusi kumulatif bersama

H(x,y) = P(Tj ≤ x, Cj ≤ y).

Distribusi Probabilitas … (Suyono) 167

Secara khusus, jika barisan-barisan (Tj) dan (Sj) saling independen maka (K(t), t ≥ 0), merupakan proses renewal reward, lihat Ross (2003).

Notasikan F dan G fungsi-fungsi distribusi kumulatif dari Tj dan Cj , yakni

F(x) = H(x,∞) = P(Tj ≤ x) dan G(y) = H(∞,y) =P(Cj ≤ y).

Transformasi Laplace Stieltjes dari F dan H masing-masing akan dinotasikan dengan

][)()(0

* jTx eExdFeF βββ −∞

− == ∫

dan

∫ ∫∞

−−∞

−− ==0 0

* ][),(),( jj CTyx eEyxdHeH αβαβαβ

dimana α, β ≥ 0.

DISTRIBUSI PROBABILITAS TOTAL BIAYA PERAWATAN

Distribusi probabilitas dari total biaya perawatan K(t) yang disajikan pada persamaan (1) diberikan pada teorema berikut dalam bentuk transformasi Laplace. Teorema 1. Untuk α, β > 0

)],(1[)(1][ *

*

0

)(

αββββα

HFdteeE ttK

−−

=∫∞

−− .

Bukti: Dengan mengkondisikan pada kejadian {T1 = s, Y1 = y} diperoleh, untuk α >0,

)](1[),(][

),(],|[

),(],|[][

0 0

)]([

011

)(

0 011

)()(

tFysdHeE

ysdHyYsTeE

ysdHyYsTeEeE

y

t

s

stKy

y ts

tK

y

t

s

tKtK

−+=

==+

===

∫ ∫

∫ ∫

∫ ∫

∞

= =

−+−

∞

=

∞

=

−

∞

= =

−−

α

α

αα

Dengan mengambil transformasi Laplace pada kedua ruas persamaan di atas diperoleh, untuk β > 0,

Jurnal Mat Stat, Vol. 8 No. 2 Juli 2008: 165-170 168

)](1[1][),(

)](1[1][),(

])(1[1),(][

)(1),(][

)](1[),(][][

*

0

)(*

*

0 0 0

)(

00 0 0

)(][

00 0

)(

00 0 0

)]([

0

)(

ββ

αβ

ββ

β

β

βα

βαβα

βαβα

βαβα

ββαβα

FdueeEH

FdueeEysdHe

tdFeysdudHeEee

dtetFysdtdHeEee

dtetFdtysdHeeEdteeE

u

uuK

y s u

uuKsy

t

y s u

uKusy

t

y s st

stKty

t

t

t y

t

s

tstKyttK

−+=

−+=

−+=

−+=

−+=

∫

∫ ∫ ∫

∫∫ ∫ ∫

∫∫ ∫ ∫

∫∫ ∫ ∫∫

∞

=

−−

∞

=

∞

=

∞

=

−−−−

∞−

∞

=

∞

=

∞

=

−+−−

∞−

∞

=

∞

=

∞

=

−−−−

∞

=

−∞

=

∞

= =

−−+−∞

−−

Sebagai akibatnya

)],(1[)(1][ *

*

0

)(

αββββα

HFdteeE ttK

−−

=∫∞

−− . ●

Untuk mendapatkan distribusi probabilitas dari K(t) transformasi Laplace di atas pada

umumnya harus diinversi secara numerik yang algoritmanya dapat dilihat pada, misalnya Abate and Whitt (1992).

Jika biaya perawatan Cj merupakan fungsi dari waktu hidup komponen, Tj, yakni jika Cj = c(Tj) untuk suatu fungsi Borel c, maka diperoleh teorema berikut. Teorema 2. Jika Cj = c(Tj) untuk suatu fungsi Borel c maka, untuk α, β > 0, diperoleh

])(1[

)(1][

0

)(

*

0

)(

∫∫ ∞ −−

∞−−

−

−=

sdFe

FdteeEssc

ttK

βα

βα

β

β . (2)

Bukti: Langkah-langkah pembuktian teorema ini serupa dengan bukti Teorema 1, hanya saja pada langkah pertama untuk menghitung transformasi Laplace dari K(t) dikondisikan pada kejadian {T1 = s} yang menghasilkan, untuk α >0,

)](1[)(][][0

)]()([)( tFsdFeEeEt

stKsctK −+= ∫ −+−− αα .

Selanjutnya dengan mengambil transformasi Laplace pada kedua ruas persamaan di atas bukti selesai. ●

MEAN DAN MOMEN KEDUA DARI TOTAL BIAYA PERAWATAN KOMPONEN

Mean dari total biaya perawatan komponen K(t) sebagaimana disajikan pada persamaan (1) dengan biaya perawatan Cj = c(Tj) untuk suatu fungsi Borel c disajikan dalam teorema berikut. Teorema 3. Untuk β > 0,

)](1[

)()()]([ *

0

0 ββ

β

β

F

sdFescdtetKE

s

t

−= ∫∫

∞ −∞−

Distribusi Probabilitas … (Suyono) 169

Bukti: Dengan mengambil turunan pertama terhadap α pada kedua ruas persamaan (2) diperoleh

2

0

)(

0

)(*

0

)(

])(1[

])()()][(1[])([

∫∫

∫ ∞ −−

∞ −−∞−−

−

−−=−

sdFe

sdFescFdteetKE

ssc

ssc

ttK

βα

βα

βα

β

β.

Bukti selesai dengan mengambil α = 0 pada persamaan di atas. ● Transformasi Laplace untuk momen kedua dari K(t) diberikan dalam teorema berikut. Teorema 4. Untuk β > 0,

2*0 0

22*

0

2

)](1[

)]()([2)()()](1[)]([

ββ

β ββ

β

F

sdFescsdFescFdtetKE

ss

t

−

+−= ∫ ∫

∫∞ ∞ −−∞

−

Bukti: Ambil turunan kedua terhadap α pada kedua ruas persamaan (2) dan selanjutnya ambil α = 0. ●

APLIKASI

Anggap waktu-waktu antar perawatan komponen, Tj, saling independen dan berdistribusi identik dengan mean 1/2, yakni

0,1)( 2 ≥−= − xexF x .

Di sini barisan waktu-waktu perawatan (Sj, j ≥ 0) merupakan proses Poisson dengan laju (rate) 2. Untuk β ≥ 0,

22)(*

+=β

βF .

Anggap biaya perawatan komponen Cj = c(Tj) dengan c(s) = es, s ≥ 0. Maka

122)()(

0

2

0 +== ∫∫

∞ −−∞ −

βββ dseeesdFesc ssss

Dan

βββ 22)()(

0

22

0

2 == ∫∫∞ −−∞ − dseeesdFesc ssss .

Dengan menggunakan hasil-hasil pada Teorema 3 dan Teorema 4 diperoleh

)1()2(2)]([ 2

0 ++

=∫∞

−

ββββ dtetKE t

dan

323

2

0

2 )2(2)1(

)2(8)]([ββ

ββββ +

+++

=∫∞

− dtetKE t

Jurnal Mat Stat, Vol. 8 No. 2 Juli 2008: 165-170 170

Kedua transformasi Laplace ini dapat diinversi secara analitik dan masing-masing menghasilkan )1(24)]([ tettKE −−−=

dan tt teetttKE −− +−++= 83863616)]([ 22 .

Dengan menggunakan Teorema 2 diperoleh

1

1

1

0

)(

!)()(1)(][

−∞

=

−∞

−−⎥⎦

⎤⎢⎣

⎡−+

−−+= ∑∫

i

ittK

iiAdteeE

λβαλλββα .

Untuk mendapatkan distribusi probabilitas dari total biaya perawatan komponen transformasi Laplace ini harus diinversi secara numerik.

PENUTUP Dalam makalah ini telah dibahas model matematika untuk total biaya perawatan komponen pada interval waktu berhingga [0,t], t ≥ 0. Waktu-waktu perawatan komponen dimodelkan dengan proses renewal. Biaya perawatan komponen dianggap hanya bergantung pada lama waktu hidup komponen. Dengan asumsi-asumsi tersebut telah diperoleh hasil berupa distribusi probabilitas total biaya perawatan dalam bentuk transformasi Laplace ganda. Selain itu dalam makalah ini juga disajikan transformasi Laplace dari mean dan momen kedua total biaya perawatan. Sebagai kelanjutan dari hasil-hasil yang telah diperoleh dalam makalah ini kiranya menarik dan penting untuk dikaji sifat-sifat asimtotik dari total biaya perawatan yang sampai sejauh ini belum penulis dapatkan.

DAFTAR PUSTAKA

Abate, J. and Whitt, W. (1992), The Fourier-series method for inverting transforms of probability distributions, Queuing Systems 10, pp. 5 – 88.

Ross, H.M. (2000), Introduction to Probability Models; Seventh Edition, Academic Press, San Diego. Van Noortwijk, J.M. and Frangopol, D.M. (2004), Two probabilistic life-cycle maintenance models

for deteriorating civil infrastructures, Prob. Eng. Mech. 19, pp. 354-359. Van Noortwijk, J.M. (2003), Explicit formulas for the variance of discounted life-cycle cost,

Reliability Engineering and System Safety 80, pp. 185 – 195. Wagner, H.M. (1975), Principle of Operations Research; Second Edition, Englewood Cliffs, New

Jersey.