Научная статья на тему 'Methods and algorithms for evaluating unknown parameters of operation processes of complex technical systems'

Methods and algorithms for evaluating unknown parameters of operation processes of complex technical systems Текст научной статьи по специальности «Математика»

CC BY
78
9
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук

Аннотация научной статьи по математике, автор научной работы — K. Kolowrocki, J. Soszynska

The paper objectives are to present the methods and tools useful in the statistical identifying unknown parameters of the operation models of complex technical systems and to apply them in the maritime industry. There are presented statistical methods of determining unknown parameters of the semi-markov model of the complex technical system operation processes. There is also presented the chi-square goodness-of-fit test applied to verifying the distributions of the system operation process conditional sojourn times in the particular operation states. Applications of these tools to identifying and predicting the operation characteristics of a ferry operating at the Baltic Sea waters are presented as well

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Methods and algorithms for evaluating unknown parameters of operation processes of complex technical systems»

METHODS AND ALGORITHMS FOR EVALUATING UNKNOWN PARAMETERS OF OPERATION PROCESSES OF COMPLEX TECHNICAL SYSTEMS

K. Kolowrocki, J. Soszynska.

Gdynia Maritime University, Gdynia, Poland e-mail: [email protected], [email protected]

ABSTRACT

The paper objectives are to present the methods and tools useful in the statistical identifying unknown parameters of the operation models of complex technical systems and to apply them in the maritime industry. There are presented statistical methods of determining unknown parameters of the semi-markov model of the complex technical system operation processes. There is also presented the chi-square goodness-of-fit test applied to verifying the distributions of the system operation process conditional sojourn times in the particular operation states. Applications of these tools to identifying and predicting the operation characteristics of a ferry operating at the Baltic Sea waters are presented as well.

1 INTRODUCTION

Many real transportation systems belong to the class of complex systems. It is concerned with the large numbers of components and subsystems they are built and with their operating complexity. Modeling the complicated system operation processes, is difficult because of the large number of the operation states, impossibility of their precise defining and because of the impossibility of the exact describing the transitions between these states. The changes of the operation states of the system operations processes cause the changes of these systems reliability structures and also the changes of their components reliability functions (Blokus-Roszkowska et all. 2008b). The models of various multistate complex systems are considered in (Blokus-Roszkowska et all. 2008a). The general joint models linking these system reliability models with the models of their operation processes (Kolowrocki, Soszynska 2008), allowing us for the evaluation of the reliability and safety of the complex technical systems in variable operations conditions, are constructed in (Blokus-Roszkowska et all. 2008b).

In order to be able to apply these general models practically in the evaluation and prediction of the reliability of real complex technical it is necessary to elaborate the statistical methods concerned with determining the unknown parameters of the proposed models, namely the probabilities of the initials system operation states, the probabilities of transitions between the system operation states and the distributions of the sojourn times of the system operation process in the particular operation states and also the unknown parameters of the conditional multistate reliability functions of the system components in various operation states. It is also necessary the elaborating the methods of testing the hypotheses concerned with the conditional sojourn times of the system operations process in particular operations states and the hypotheses concerned with the conditional multistate reliability functions of the system components in the system various operation states. The model of the operation process of the complex technical system with the distinguished their operation states is proposed in (Kolowrocki, Soszynska 2008). The semi-markov process is used to construct a general probabilistic model of the considered complex industrial system operation process. To construct this model there were defined the vector of the probabilities of the system initial operation states, the

matrix of the probabilities of transitions between the operation states, the matrix of the distribution functions and the matrix of the density functions of the conditional sojourn times in the particular operation states. To describe the system operation process conditional sojourn times in the particular operation states the uniform distribution, the triangular distribution, the double trapezium distribution, the quasi-trapezium distribution, the exponential distribution, the Weibull's distribution, the normal distribution and the chimney distribution are suggested in (Kolowrocki, Soszynska 2008). In this paper, the formulae estimating unknown parameters of these distributions are given and the chi-square test is applied to verifying the hypotheses about these distributions validity. Moreover these tools are applied to unknown parameters estimation and characteristics prediction of the Stena Baltica ferry operation process.

2 IDENTYFICATION OF THE OPERATION PROCESS OF THE COMPLEX TECHNICAL SYSTEM

2.1. Estimation of unknown parameters of the semi-markov model of the operation process

We assume, similarly as in (Blokus-Roszkowska et all 2008b, Kolowrocki, Soszynska 2008) that a system during its operation at the fixed moment t, t e <0, + » >, may be in one of v, v e N, different operations states zb, b = 1,2, ..., v. Next, we mark by Z(t), t e <0, +<»>, the system operation process, that is a function of a continuous variable t, taking discrete values in the set Z = {z1,z2,...zv } of the operation states. We assume a semi-markov model (Blokus-Roszkowska et all 2008b, Kolowrocki, Soszynska 2008) of the system operation process Z(t) and we mark by 6bl its random conditional sojourn times at the operation states zb when its next operation state is z,, b, l = 1,2,..., v, b * l.

Under these assumption, the operation process may be described by the vector [pb (0)]1xv of probabilities of the system operation process staying in particular operations states at the initial moment t = 0, the matrix [pbl (t)]vxv of the probabilities of the system operation process transitions

between the operation states and the matrix [Hbl (t)]vxv of the distribution functions of the conditional sojourn times 9bl of the system operation process at the operation states or equivalently by the matrix [hbl (t)]vxv of the distribution functions of the conditional sojourn times 9bl of the system operation process at the operation states.

To estimate the unknown parameters of the system operations process, firstly during the experiment, we should collect necessary statistical data performing the following steps:

i) to analyze the system operation process and either to fix or to define its following general parameters:

- the number of the operation states o the system operation process v,

- the operation states o the system operation process z1, z2, ..., zv,

- the duration time of the experiment 0 ;

ii) to fix and to collect the following statistical data necessary to evaluating the probabilities of the initial states of the system operations process:

- the number of the investigated (observed) realizations of the system operation process n(0),

- the numbers of staying of the operation process respectively in the operations states z1, z2, ..., zv, at the initials moment t = 0 of all n(0) observed realizations of the system operation process

n1 (0), «2(0), «v(0);

iii) to fix and to collect the following statistical data necessary to evaluating the probabilities of transitions between the system operation states:

- the numbers nbl, b, l = 1,2,...,v, b^l, of the transitions of the system operation process from the operation state zb to the operation state zl during all observed realizations of the system operation process;

- the numbers nb, b = 1,2,...,v, of departures of the system operation process from the operation states zb ;

iv) to fix and to collect the following statistical data necessary to evaluating the unknown parameters of the distributions of the conditional sojourn times of the system operation process in the particular operation states:

- the realizations 9kbI, k = 1,2, ..., nbl, for each b, l = 1,2,...,v, b * l of the conditional sojourn times 9bI of the system operations process at the operation state zb when the next transition is to the operation state zl during the observation time;

After collecting the above statistical data it is possible to estimate the unknown parameters of the system operation process performing the following steps: i) to determine the vector

[ p(0)] = [ pY (0), p2(0),..., pv (0)],

of the realizations of the probabilities pb (0), b = 1,2,..., v, of the initial states of the system operation process, according to the formula

(0)=

«(0)

for b = 1,2,...,

where n(0) = £ nb (0), is the number of the realizations of the system operation process starting at the

4=1

initial moment t = 0; ii) to determine the matrix

[ Pi ] =

Pll Pl2 ••• Plv

P21 p22 . . . p2v

Pvl Pv2 ... Pvv

of the realizations of the probabilities pu, b, l = 1,2,..., v, of the system operations process transitions from the operations state zb to the operations state zl during the experiment time according to the formula

n

pbl =-iL for b,l = 1,2,..., v, b * l, pbb = 0 for b = 1,2,..., v,

where nb = £ nbl, b = 1,2,..., v, is the realization of the total number of the system operations process

b*l

departures from the operations state zb during the experiment time ©;

iii) to determine the following empirical characteristics of the realizations of the conditional sojourn time of the system operation process in the particular operation states:

V

b

- the realizations of the mean values 9bl of the conditional sojourn times 9bl of the system operations process at the operations state Hbl (t) when the next transition is to the operation state 9bl , according to the formula

— 1 "bl , , , ebl = — i ekbl b,i = 1,2,..., v, b * I,

- the number 9kM, of the disjoint intervals k = 1,2,...,nbl,, 9bl, that include the realizations 9^, k = 1,2,...,nbl, of the conditional sojourn times 9bl at the operation state zb when the next transition is to the operation state zl, according to the formula

the length d of the intervals Ij =< aj,bj), j = 1,2,..., r , according to the formula

R

d = , where R = max 6^ — min 6^ .

r — 1 1<k <nbl 1Sk Snbl

the ends aj, bj, of the intervals Ij =< aJbl,bj), j = 1,2,..., r , according to the formulae

ab = min 6 — d, K = ah + jd, j = 1,2,...,r , abi = bb— , j = 2,3,..., r,

1<k <nbl 2

in the way such that

/1 uI2 u... uIr =< ah,bbi

and

It n Ij = 0 for all i * j, i, j e {1,2,..., r},

- the numbers nj of the realizations 9^ in particular intervals Ij, j = 1,2,...,r , according to the formula

nj = #{k : 9k e Ij, k e {1,2,...,"bl}}, j = 1,2,..., r ,

where I nj = nbl, whereas the symbol # means the number of elements of the set;

j=1

iv) to estimate the parameters of the distributions of the conditional sojourn times of the system operation process in the particular operation states for the following distinguished distributions respectively in the following way: - the uniform distribution with a density function

nbl k=1

r = Vnbl

hbl (t) =

t < x.

bl

1

> xbl < t < ybi

ybl xbl 0, t > ybi,

where 0 < xbl < ybl <

the estimates of the unknown parameters of this distribution are:

xbl = abl , ybl = Xbl + rd ; - the triangular distribution with a density function

i0,

hbl (t) =

t < x

2

t-xb

ybl xbl zbl xbl 2 ybl-t

ybl xbl ybl zbl 0,

bl

xbl < t < zbl

zbl < t < ybl

t > ybl,

where 0 < xbl < zbl < ybl <+<»,

the estimates of the unknown parameters of this distribution are:

xbl = abi, ybl = xbl + rd, i = r;

- the double trapezium distribution with a density function

f0,

hu (t) =

t < x

bl

qu +

wb, +

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Cbl -qbl

zbl ■ xbl

Cbl -Wbl

ybl -zbl

(t-xblX xb, < t < Zbl -(yu-tX zu < t < ybl

t > ybl :

where

Cbl =

2- qu (zu -xbl)- wbl(yu-zbi) ybl - xbl

, 0 < Xbl < Zbl < ybi < +<», 0 < qu < +», 0 < Wbl < +<»,

0 < qu (zbi ~xbt) + wbl(ybl ~zbl) < 2;

the estimates of the unknown parameters of this distribution are:

n n

xbl = ah, ybi = xbl +rd, qbl = , =, = Ob

nbld ' nbld'

-bl "bl'

0

0

- the quasi-trapezium distribution with a density function

hbl (t) =

qbl + (t~xbl),

zbl xbl

Wbl +

ybl zbl

t < xbl xbl < t < zh

1 ^ -f ^ 2 zbl < t < zbl

2

zbl < t < ybl

t > ybl,

where

A _

2 qbl (zbl xbl ) wbl (ybl zbl ) 0 < X < z 1 < z 2 < y < 0 < q <

-T-z---,0 < Xbl < zbl < zbl < ybl > 0 < qbl <+'x>,

zbl zbl + ybl Xbl

0 < Wbl < 0 < qbl(4 ~xbl) + wbl (zh -ybl) < 2 ,

the estimates of the unknown parameters of this distribution are:

xbl = abl , ybl = xbl + rd, qbl =

nbl w _ nbl _ 1 _ Q\ _ 2 _ 02

- wbl _-7' zbl _ 0bl, zbl _

nbld

nbld'

where

— 1 1 n(me) . — 2 °bi _- E , _

nbl n,, + 1

E0. , nme) _ [^];

n(me) 1 _1 nbl - n(me) J_n(me) +1 2

- the exponential distribution with a density function

0,

hbl (t) _\ '

t < Xbl, Xbl > 0

abl exp[-abl (t-xbl)], t > xbl,

where 0 < abl < +<x>, 0 < xbl = aj,,

the estimates of the unknown parameters of this distribution are:

1

Xbl _ abl , abl _ —

0bl Xbl

- the Weibull's distribution with a density function

h (t) _j°' t < Xbl'Xbl > 0' " \außu (t - Xbl)ßbl-1 exp[-abl (t - Xbl)ßbl ], t > Xbl,

where 0 < abl < +<x>, 0 < /3bl < +<x>, 0 < xbl = albl,,

the estimates of the unknown parameters of this distribution are:

0

bl

0

1

«m «bl .

bl- + X Inidj-Xbr )

Xbl = abl, abl = «

abl =

ßbl i=1

X (ßj )

j=l

j \ßbl

nbl X

j=l

Xßl)ßbl lnßbl -Xbl)

- the normal distribution with a density function

hbl (t) =

exp[-(t mb\)2], t e (-œ, œ), 2a,,

where - <x> < mu < +<x>, 0 < abi < +<x>,

the estimates of the unknown parameters of this distribution are:

— 2 _2 l nbl . 2

mbi = ßbl, abl = abl =-X (ß -mbl ) ,

- the chimney distribution with a density function

hbl (t) =

zbl Xbl

,2 1

db

ybl zbl 0,

t < X

bl

Xbl < t < zbl

1 ^ ^ 2 zbl < t < zbl

zh < t < ybl

t > ybl,

where 0 < xbl < z\, < zl, < ybl < 0 < qbl < 0 < wbl < +M, abl > cbl > 0, dbl > 0, abl + cbl + dbl = b

The estimates of the unknown parameters of this distribution are:

xbl = abl , ybl = xbl + rd,

and moreover, if

nbl = max{njl}, i = j, where j e{12,..., r},

1< j<r

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

is such a number of the interval for which nj = nbl, then: for i = 1

« ii

zh = Xbl + (i - 1)d, 4 = Xbl +id, abl = 0, cbl =—, dbl = bl

«b+1 +...+«b

n

bl

1

0

a

bl

c

bl

n

bl

bl

n

when n'b+1 _ 0 or n'b+1 * 0 and > 3,

n' + nM ri+2 + + ri

zlbl _ Xbl + (i - 1)d, 4 _ Xl + (i + 1)d, abl _ 0, cbl _ -, dbl _-»—■■ b'

n

bl

when rib+1 * 0 and < 3,

bl n'+l

nbl

for i _ 2,3,..., r-1

nl, +... + nb-1 ni , nb+1 +... + nb

zu _ Xbl + (i 1)d, zbl _ Xbl +id, abl _ , cbl _ , dbl

nbl nbl nbl

n

when n'b-1 _ 0 or n'b-1 * 0 and —^ > 3 and

nb-

nb +... + rib-1 n[, + ni+1 nb+2 +... + nb

zu _ Xb, + (i - 1)d, 4 _ XU + (i + 1)d, abl _ b' , cu -^, du _ b'

nbl nbl nbl

when n'bl 1 _ 0 or ribl 1 * 0 and —-- > 3 and

n'bl

4 _ Xbl + (i - 2)d, zh _ Xbl +id, abl _

n\, +... + n'u2 ^ n'b- + nl , _nb+1 +... + nh

bl ' "bl ' "bl nbl nbl

nn when nb - * 0 and -b- < 3 and when n'b+1 _ 0 or ri^1 * 0 and > 3,

nb- nb+

nb +... + nb-1 nbl1 + ni, + n':1 , nb,2 +... + nb

„2 „ , /-,. , i\ j bl ' ■■■ ' "bl "bl ' ,lbl ' /7 ,lbl 1 ••• 1 ,lbl

zu _ Xbl + (i- 2)d, zu _ Xbl + (i + 1)d, abl _-, cbl _-, dbl _■

nbl nbl nbl

when n'bi1 * 0 and —-- < 3 and when n'b+1 * 0 and —+1 < 3,

nb, nb+

for i _ r

n\, +... + nb1

4 _ Xbl + (i - 1)d, zh _ Xbl + id, abl _ -—, cbl _—, dbl

nn

nbl

n

when nb- _ 0 or nb- * 0 and -JbL- > 3 ,

n'bl

= xbl + (i - 2)d,

Zbl = xbl +

nhl +... + n

bl

nb-1 + n

Cu, =■

bl

dM = 0,

i-2

abl =

z

n

n

bl

bl

n

when nl-1 ^ 0 and < 3

nbl

2.2. Identification of distributions of conditional sojourn times in operation states

To formulate and next to verify the non-parametric hypothesis concerning the form of the distribution function Hbl (t) of the system conditional sojourn time 9bl at the operation state zb

when the next transition is to the operation state zl, on the basis of its realizations 9kbi, k = 1,2,...,nbl, it is necessary to proceed according to the following scheme:

- to construct and to plot the realization of the histogram of the system conditional sojourn time 9bl at the operation states, defined by the following formula

— nJ

K (t) = -f for t elj,

nbi

- to analyze the realization of the histogram, comparing it with the graphs of the density functions hbl (t) of the previously distinguished distributions, to select one of them and to formulate the null hypothesis H0 and the alternative hypothesis HA, concerning the unknown form of the distribution function Hbl(t) of the conditional sojourn time 9bl in the following form:

H0 : The system conditional sojourn time 9bl at the operation state zb when the next transition is to the operations state zl, has the distribution function Hbl (t),

HA : The system conditional sojourn time 9bl at the operation state zb when the next transition is to the operations state zl, has the distribution function different from Hbl (t),

- to join each of the intervals Ij that has the number nj of realizations is less than 4 either with the neighbour interval Ij+1 or with the neighbour interval Ij-1 this way that the numbers of realizations

in all intervals are not less than 4,

- to fix a new number of intervals r ,

- to determine new intervals Ij =< aj,bj), j = 1,2,.., r,

- to fix the numbers nj of realizations in new intervals ^, j = 1,2,.., r,

- to calculate the hypothetical probabilities that the variable 9bl takes values from the interval Ij, under the assumption that the hypothesis H0 is true, i.e. the probabilities

Pj = P(9bl e Ij) = P(aj < 9bl < bj ) = Hbl (bbj) -Hbl (aj), j = 1,2,..., f,

where Hbl (bbj) and Hbl (aj) are the values of the distribution function Hbl (t) of the random variable 9bl defined in the null hypothesis H0,

- to calculate the realization of the %2 (chi-square)-Pearson's statistics Unu, according to the formula

= £ (nbl -nblPj f J=1 n., p j

u

- to assume the significance level a (a = 0.01, a = 0.02, a = 0.05 or a = 0.10) of the test,

- to fix the number f-l-1 of degrees of freedom, substituting for l for the distinguished distributions respectively the following values: l = 0 for the uniform, triangular, double trapezium, quasi-trapezium and chimney distributions, l = 1 for the exponential distribution, l = 2 for the Weibull's and normal distributions,

- to read from the Tables of the %2 - Pearson's distribution the value ua for the fixed values of the significance level a and the number of degrees of freedom f-l -1 such that the following equality holds

P(Unu > Ua) = 1 -a,

and next to determine the critical domain in the form of the interval (ua,+<x>) and the acceptance domain in the form of the interval < 0, ua >,

- to compare the obtained value unu of the realization of the statistics U with the red from the Tables critical value ua of the chi-square random variable and to verify previously formulated the null hypothesis H0 in the following way: if the value unbi does not belong to the critical domain, i.e. when unu < ua,then we do not reject the hypothesis H0, otherwise if the value unbi belongs to the critical domain, i.e. when unbi > ua, then we reject the hypothesis H0 in favor of the hypothesis HA .

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

3 APPLICATION IN MARITIME TRANSPORT

3.1. The Stena Baltica ferry operation process and its statistical identification

Taking into account the operation process of the considered ferry we distinguish the following as its eighteen operation states:

• an operation state z1 - loading at Gdynia Port,

• an operation state z2 -unmooring operations at Gdynia Port,

• an operation state z3 - leaving Gdynia Port and navigation to "GD" buoy,

• an operation statez4-navigation at restricted waters from "GD" buoy to the end of Traffic Separation Scheme,

• an operation state z5 - navigation at open waters from the end of Traffic Separation Scheme to "Angoring" buoy,

• an operation state z6 - navigation at restricted waters from "Angoring" buoy to "Verko" Berth at Karlskrona,

• an operation state z7 - mooring operations at Karlskrona Port,

• an operation state z8 - unloading at Karlskrona Port,

• an operation state z9 - loading at Karlskrona Port,

• an operation state z10 - unmooring operations at Karlskrone Port,

• an operation state z11 - ship turning at Karlskrone Port,

• an operation state z12 - leaving Karlskrone Port and navigation at restricted waters to "Angoring" buoy,

• an operation state z13 - navigation at open waters from "Angoring" buoy to the entering Traffic Separation Scheme,

• an operation state z14 -navigation at restricted waters from the entering Traffic Separation Scheme to "GD" buoy,

• an operation state z15 - navigation from "GD" buoy to turning area,

• an operation state z16 - ship turning at Gdynia Port,

• an operation state z17 - mooring operations at Gdynia Port,

• an operation state z18 - unloading at Gdynia Port.

To identify all parameters of Stena Baltica ferry operation process the statistical data about this process is needed. The statistical data that has been collected up to now is given in Tables 1-7 in (Soszynska et all 2009, Appendix 5A). In the Tables 1-7 there are presented the realizations 9kbI, k = 1,2,...,42, for each b = 1,2,...,17, l = b +1 and b = 18, l = 1 of the ship operation process conditional sojourn times 9bl, b = 1,2,...,17, l = b +1 and b = 18, l = 1 in the state zb while the next transition is the state zt during the experiment time 0 = 42 days.

These statistical data allow us, applying the methods and procedures given in the section 2, to formulate and to verify the hypotheses about the conditional distribution functions Hbl (t) of the Stena Baltica ferry operation process sojourn times 9bl, b = 1,2,...,17, l = b +1 and b = 18, l = 1 in the state zb while the next transition is to the state zl on the base of their realizations 9kbI, k = 1,2,...,42. On the basis of the statistical data, given in the Appendix 5A in [5], the vector of the probabilities of the system initial operation states was evaluated in the following form

[ Pb (0)] = [1, 0, 0,..., 0, 0].

The matrix of the probabilities pbl of transitions from the operation state zb into the operation state zl were evaluated as well. Their evaluation are given in the matrix below

[ Pbl ]1

010...00 001...00

000...01 100...00

Next the matrix [hbl (t)]18118 of conditional density functions of the system operation process Z(t)

conditional sojourn times 9bl (Soszynska et all 2009, Appendix 5A) were evaluated.

The results of the distributions unknown parameters estimation and the hypotheses testing are as

follows:

- the conditional sojourn time 912 have a triangular distribution with the density function

¿12 (t) = \

0, t < 7,

0.000441- 0.0031, 7 < t < 54, 0.044 -0.000431; 54 < t < 103, 0, t > 103;

- the conditional sojourn time 923 have an exponential distribution with the density function

Mt) =

0,

t < 1.6

1.03exp[-1.03(t-1.6)], t > 1.6;

- the conditional sojourn time d34 have a steep chimney distribution with the density function

MO =

0, t < 29,

0.0278, 29 < t < 35,

0.1984, 35 < t < 38,

0.0266, 38 < t < 47,

0, t > 47;

- the conditional sojourn time 6>45 have a chimney distribution with the density function

MO =

0, t < 41,

0.0095, 41 < t < 46,

0.0762, 46 < t < 56,

0.0127, 56 < t < 71,

0, t > 71;

- the conditional sojourn time d56 have a double trapezium distribution with the density function

MO =

0, t < 467.8,

- 0.00004t + 0.0277, 467.8 < t < 525.95,

- 0.00006t + 0.0397, 525.95 < t < 650.2, 0, t < 650.2;

- the conditional sojourn time 067 have a double trapezium distribution with the density function

Mt) =

0, t < 31.9,

0.0067t-0.1747, 31.9 < t < 37.17, 0.0031t-0.0395, 37.7 < t < 45.1, 0, t > 45.1;

- the conditional sojourn time 6>78 have a double trapezium distribution with the density function

K(t) =

0, t < 4.5,

- 0.0183t + 0.2922, 4.5 < t < 7.02,

- 0.0069t + 0.2122, 7.02 < t < 10.5, 0, t > 10.5;

- the conditional sojourn time 6>89 have a triangular distribution with the density function

h»(t) =

0, t < 0,

0.0021t, 0 < t < 21.4

- 0.002t + 0.087, 21.4 < t < 44.4,

0, t > 44.4;

- the conditional sojourn time 6>910 have a double trapezium distribution with the density function

h910(t) =

0, t < 14.6,

0.0001t + 0.0109, 14.6 < t < 52.26,

0.0062, 52.2 < t < 127.4,

0, t > 127.4;

- the conditional sojourn time 6>1011 have a double trapezium distribution with the density function

h1011(t)

0, t < 1.6, - 0.3071t +1.0014, 1.6 < t < 2.93, 0.0398t -0.0145, 2.93 < t < 6.4, 0, t > 6.4;

the conditional sojourn time 9nn have a quasi-trapezium distribution with the density function

h1112(t) =

0, t < 3.8,

-148.899t + 567.4863, 3.8 < t < 3.81, 0.181, 3.81 < t < 4.48,

0.3773t -1.5094, 4.48 < t < 6.2, 0, t > 6.2;

- the conditional sojourn time 01213 have a triangular distribution with the density function

h1213(t) =

0, t < 18.7,

0.025t -0.4675, 18.7 < t < 23.86,

- 0.012t + 0.4116, 23.86 < t < 34.3,

0, t > 34.3;

- the conditional sojourn time 6>1314 have a chimney distribution with the density function

h1314(t) =

0, t < 410,

0.0017, 410 < t < 478,

0.0189, 478 < t < 512,

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

0.0024, 512 < t < 614,

0, t > 614;

- the conditional sojourn time 6>1415 have a double trapezium distribution with the density function

h1415(t) =

0, t < 36.8,

-0.0006t + 0.0518, 36.8 < t < 50.14, 0.0003t + 0.0062, 50.14 < t < 75.2, 0, t > 75.2;

- the conditional sojourn time 91516 have a chimney distribution with the density function

h1516(t) =

0, t < 30,

0.0317, 30 < t < 33,

0.2698, 33 < t < 36,

0.0084, 36 < t < 48,

0, t > 48;

the conditional sojourn time 91617 have a triangular distribution with the density function

10,

h1617(t) =

t < 2.7,

0.305t -0.823, 2.7 < t < 4.52, -0.313t +1.9719, 4.52 < t < 6.3, 0, t > 6.3;

- the conditional sojourn time 91718 have a double trapezium distribution with the density function

h1718(t) =

0, t < 2.3,

- 0.1134t + 0.6707, 2.3 < t < 5.62, 0.0071t - 0.0063, 5.62 < t < 10.7, 0, t > 10.7;

- the conditional sojourn time 9181 have a triangular distribution with the density function

h181(t) =

0, t < 0,

0.0023t, 0 < t < 18.74,

-0.0016t + 0.0729, 18.74 < t < 45.59,

0, t > 45.59.

3.2. The Stena Baltica ferry operation process prediction

On the basis of the previous section, the mean values Mbl = E[9bl ], b,l = 1,2,...,18, b ^ l, (12) in (Kolowrocki, Soszynska 2008) of the system operation process Z(t) conditional sojourn times in particular operation states were determined and there are given by:

M12 = 54.33, M23 = 2.57, M34 = 36.57,

M45 = 52.5, M56 = 525.95, M67 = 37.16,

M78 = 7.02, M89 = 21.43, M910 = 53.69,

M1011 = 2.93, M1112 = 4.38, M1213 = 23.86,

M1314 = 509.69, M1415 = 50.14, M1516 = 34.28,

M1617 = 4.52, M1718 = 5.62, M181 = 18.74.

Hence, by (21) in (Kolowrocki, Soszynska 2008), the unconditional mean sojourn time in the particular operation states are given by:

M1 = E[91 ] = p12 M12 = 1 • 54.33 = 54.33,

M 2 = E[92 ] = p23M 23 = 1 • 2.57 = 2.57,

M 3 = E[93] = p34 M 34 = 1 • 36.57 = 36.57,

M 4 = E [94] = p45M45 = 1 • 52.5 = 52.5,

M 5 = E[95 ] = p56 M 56 = 1 • 525.95 = 525.95,

M 6 = E9] = P67 M 67 = 1 • 37.16 = 37.16,

M 7 = E[97] = p78 M 78 = 1 • 7.02 = 7.02,

M 8 = E9] = pggMs9 = 1 • 21.43 = 21.43,

M 9 = E[99] = P910 M- 910 = 1 • 53.69 = 53.69,

M!0 = E[9j0] = pV)nMvm = 1 • 2.93 = 2.93,

M11 = E[9n] = pm2 M1112 = 1 • 4.38 = 4.38,

M12 = E[912] = Pj213 M1213 = 1 • 23.86 = 23.86,

M13 = E[913] = p1314M1314 = 1 • 509.69 = 509.69,

M14 = E[914] = p1415 M1415 = 1 • 50.14 = 50.14,

M !5 = E[9^5 ] = p1516 M ,5,6 = 1 • 34.28 = 34.28,

M16 = E[9|6 ] = P1617 1617 = 1 • 4.52 = 4.52,

M17 = E[917] = p1718 M1718 = 1 • 5.62 = 5.62,

M18 = E[918] = P181M181 = 1 • 18.74 = 18.74.

Since from the system of equations below (23) in (Kolowrocki, Soszynska 2008) that takes the form

n1 = n2 = ... = n18 = 0.056,

then the limit values of the transient probabilities (the portions of time of a week, as the operation process is periodic) pb (t) at the operational states zb, according to (22) in (Kolowrocki, Soszynska 2008), are given by

p1 = 0.037, p2 = 0.002, p3 = 0.025,

p4 = 0.036, p5 = 0.364, p6 = 0.025,

p7 = 0.005, p8 = 0.014, p9 = 0.037,

p10 = 0.002, p11 = 0.003, p12 = 0.017,

p13 = 0.354, p14 = 0.035, p15 = 0.024,

p16 = 0.003, p17 = 0.004, p18 = 0.013.

Hence by (26) in (Kolowrocki, Soszynska 2008), the mean values of the system operation process total sojourn times 9b in the particular operation states zb, for the operation time 0 = 1 month = 720 hours are approximately given by

we get

E[9J = pl9 = 26.64, E[92] = p29 = 1.44,

E[93] = p9 = 18.00, E[94] = p9 = 25.92,

E[95] = p9 = 262.08, E[96] = p69 = 18.00,

E[97] = p9 = 3.6, E[98] = p9 = 10.08,

E[99] = p99 = 26.64, E[910] = pw9 = 1.44,

E[9„] = pn9 = 2.16, E[912] = pn9 = 12.24,

E[9j3] = p139 = 25.88, E[9,4] = pu9 = 25.20,

E[915] = p159 = 17.28, E[916] = p169 = 2.16,

E[0„] = pxie = 2.88, E[0J = pKd = 9.36.

4 CONCLUSION

The statistical methods and algorithms for the unknown parameters of the operation process of complex technical systems in variable operation conditions are proposed. Next, these methods are applied to estimating the operation process of Stena Baltica ferry operating between Gdynia Port in Poland and Karsklone Port in Sweden. The proposed methods other very wide applications to port and shipyard transportation systems operation processes characteristics evaluation are obvious. The results are expected to be the basis to the reliability and safety of complex technical systems optimization and their operation processes effectiveness and cost analysis

Acknowledgements

The paper describes part of the work in the Poland-Singapore Joint Research Project titled "Safety and Reliability of Complex Industrial Systems and Processes" supported by grants from the Poland's Ministry of Science and Higher Education (MSHE grant No. 63/N-Singapore/2007/0) and the Agency for Science, Technology and Research of Singapore (A*STAR SERC grant No. 072 1340050).

5 REFERENCES

Blokus-Roszkowska, A., Guze, S., Kolowrocki, K., Kwiatuszewska-Sarnecka, B., Soszynska, J. 2008a. Models of safety, reliability and availability evaluation of complex technical systems related to their operation processes. WP 4 - Task 3.1 - 31.05.2008. Poland-Singapore Joint Research Project, 20072010.

Blokus-Roszkowska, A., Guze, S., Kolowrocki, K., Kwiatuszewska-Sarnecka, B., Soszynska, J. 2008b.

Models of safety, reliability and availability evaluation of complex technical systems related to their operation processes. WP 4 - Task 4.1 - 31.05.2008. Poland-Singapore Joint Research Project, 20072010.

Blokus-Roszkowska, A. Guze, S. Kolowrocki, K., Kwiatuszewska-Sarnecka, B., Milczek B., Soszynska, J. 2009. Data Mining for Identification and Prediction of Safety and Reliability Characteristics of Complex Industrial Systems and Processes. WP6 - Task 6.1. Methods-algorithms for evaluating unknown parameters of system operation processes, reliability and safety models. WP6 - Task 6.1.2 -16.03.2009. Poland-Singapore Joint Research Project, 2007-2010. Grabski, F. 2002. Semi-Markov Models of Systems Reliability and Operations Analysis. Warsaw: Systems

Research Institute, Polish Academy of Science. Kolowrocki, K., Soszynska, J. 2008. A general model of technical systems operation processes related to their environment and infrastructure. WP 2 - Task 2.1 - 31.05.2008. Poland-Singapore Joint Research Project, 2007-2010.

Soszynska, J., Kolowrocki, K., Kaminski, P., Jurdzinski, M., Guze, S., Milczek, B., Golik, P. 2009. Data Mining for Identification and Prediction of Safety and Reliability Characteristics of Complex Industrial Systems and Processes. WP6- Task 6.2. Preliminary safety and risk analysis and evaluation of a Stena Baltica ferry in variable operation conditions. Sub-Task 6.2.5 - Sub-Task 6.2.5 - 16.03.2009. Poland-Singapore Joint Research Project, 2007-2010. Kolowrocki, K., Kwiatuszewska- Sarnecka, B. 2009. Data Mining for Identification and Prediction of Safety and Reliability Characteristics of Complex Industrial Systems and Processes. WP6 - Task 6.1. Methods-algorithms for evaluating unknown parameters of system operation processes, reliability and safety models - a preliminary study. WP6 - Task 6.1.1 - 16.03.2009. Poland-Singapore Joint Research Project, 2007-2010.

i Надоели баннеры? Вы всегда можете отключить рекламу.