Military Reliability ModelingWilliam P. Fox, Steven B. HortonIntroductionYouare an infantry rifle platoon leader.Yourplatoon is occupying a battleposition and has been ordered to establishan observationpost (OP)on ahilltopapproximatelyonekilometerforwardofyourposition.TheOPwill beoccupiedbythreesoldiersfor24hours.Hourlysituationreports mustbe madebyradio.All necessaryrations,equipment,andsuppliesforthe24hourperiod mustbecarried withthem. TheOPisineffectiveunlessitcancommunicatewithyouin a timelymanner.Therefore,radiocommunicationsmustbereliable.Theradiohasseveral componentswhichaffectitsreliabilityanessentialonebeingthebatteryBatterieshave a useful life whichis notdeterministic(wedo notknowexactlyhowlongabatterywill lastwhenweinstallit).Itslifetimeisavariablewhichmaydependonprevioususe,manufacturingdefects,weather,etc.Thebatterythatisinstalled intheradiopriortoleavingfortheOP could lastonlyafewminutesorfortheentire24hours.Sincecommunicationsaresoimportanttothismission,we are interested in modelingand analyzing the reliability of the battery.We will use the followingdefinition forreliability:If Tisthetimetofailure ofa component ofa system,and f(t)is theprobabilitydistributionfunction ofT,thenthe components'reliabilityattimetisR(t) = P(T >t) = 1 - F(t).R(t)iscalledthereliabilityfunctionandF(t)isthecumulativedistributionfunctionof f(t).Ameasure of this reliabilityis theprobabilitythata givenbattery will last morethan24hours.Ifweknowtheprobabilitydistributionforthebatterylife,wecanuseourknowledgeofprobabilitytheorytodeterminethereliability.Ifthebatteryreliability is below acceptable standards,one solution is to have the soldiers carryspares.Clearly,themore sparesthey carry,theless likelythereistobeafailureincommunicationsduetobatteries.Ofcoursethebatteryisonlyonecomponentof theradio.Others include the antenna,handset, etc.Failure of any one of theessential components causesthesystemtofail.153
153 Military Reliability Modeling William P. Fox, Steven B. Horton Introduction You are an infantry rifle platoon leader. Your platoon is occupying a battle position and has been ordered to establish an observation post (OP) on a hilltop approximately one kilometer forward of your position. The OP will be occupied by three soldiers for 24 hours. Hourly situation reports must be made by radio. All necessary rations, equipment, and supplies for the 24 hour period must be carried with them. The OP is ineffective unless it can communicate with you in a timely manner. Therefore, radio communications must be reliable. The radio has several components which affect its reliability, an essential one being the battery. Batteries have a useful life which is not deterministic (we do not know exactly how long a battery will last when we install it). Its lifetime is a variable which may depend on previous use, manufacturing defects, weather, etc. The battery that is installed in the radio prior to leaving for the OP could last only a few minutes or for the entire 24 hours. Since communications are so important to this mission, we are interested in modeling and analyzing the reliability of the battery. We will use the following definition for reliability: If T is the time to failure of a component of a system, and f(t) is the probability distribution function of T, then the components’ reliability at time t is R(t) = P (T >t) = 1 - F(t). R(t) is called the reliability function and F(t) is the cumulative distribution function of f(t). A measure of this reliability is the probability that a given battery will last more than 24 hours. If we know the probability distribution for the battery life, we can use our knowledge of probability theory to determine the reliability. If the battery reliability is below acceptable standards, one solution is to have the soldiers carry spares. Clearly, the more spares they carry, the less likely there is to be a failure in communications due to batteries. Of course the battery is only one component of the radio. Others include the antenna, handset, etc. Failure of any one of the essential components causes the system to fail
Thisisarelativelysimpleexampleofoneofmanymilitaryapplicationsofreliability.Thischapterwill showwecanuseelementaryprobabilitytogeneratemodels that can be used todeterminethe reliability of military equipmentModeling ComponentReliabilityInthis section,wewill discuss howtomodel component reliability.Recallthatthe reliability function, R(t) is defined as:R(t)= P(T > t) = P(component fails after time t).This canalsobe stated, using Tas thecomponent failure time,asR(t)= P(T> t) = 1- P(T ≤t)= 1- [f(x)dx = 1- F(t) Thus,if weknowtheprobabilitydensityfunction f(t)ofthetimetofailureT,wecan use probability theoryto determinethereliability function R(t).Wenormallythinkofthesefunctions asbeingtimedependent;however,thisisnotalways thecase.Thefunctionmightbe discrete such asthe lifetime ofa cannontube.It isdependentonthenumberofroundsfiredthroughit(adiscreterandomvariable)Auseful probabilitydistribution in reliabilityis theexponential distribution.Recallthat its densityfunction isgivenby[ae-st>0f(t) :otherwise'owhere the parameter a is such that Y, equals the mean of the random variableT. If T denotes the time to failure of a piece of equipment or a system, then isthemeantimetofailurewhichis expressed in units oftime.For applications ofreliability,we will usetheparameter α.Since/is themeantimetofailure,2istheaveragenumberoffailures perunit time orthefailurerate.Forexample, ifa light bulb has a time to failure thatfollows an exponential distribution with ameantimetofailureof50hours,thenitsfailurerateis1lightbulbper50hoursor1/50perhour,sointhiscase=0.02perhour.NotethatthemeanofT,themean time to failure of the component, is Ya.Example1:Let's considertheexamplepresented inthe introduction.Lettherandom variable T bedefined as follows:154
154 This is a relatively simple example of one of many military applications of reliability. This chapter will show we can use elementary probability to generate models that can be used to determine the reliability of military equipment. Modeling Component Reliability In this section, we will discuss how to model component reliability. Recall that the reliability function, R(t) is defined as: R(t) = P(T > t) = P(component fails after time t). This can also be stated, using T as the component failure time, as R(t) P(T t) P(T t) f ( x)dx F(t) t = > = - £ = - = - - ¥ ò 1 1 1 . Thus, if we know the probability density function f (t) of the time to failure T, we can use probability theory to determine the reliability function R(t). We normally think of these functions as being time dependent; however, this is not always the case. The function might be discrete such as the lifetime of a cannon tube. It is dependent on the number of rounds fired through it (a discrete random variable). A useful probability distribution in reliability is the exponential distribution. Recall that its density function is given by f t e t t ( ) = ì í î > - l l 0 0 otherwise , where the parameter l is such that 1 l equals the mean of the random variable T. If T denotes the time to failure of a piece of equipment or a system, then 1 l is the mean time to failure which is expressed in units of time. For applications of reliability, we will use the parameter l . Since 1 l is the mean time to failure, l is the average number of failures per unit time or the failure rate. For example, if a light bulb has a time to failure that follows an exponential distribution with a mean time to failure of 50 hours, then its failure rate is 1 light bulb per 50 hours or 1/50 per hour, so in this case l = 0.02 per hour. Note that the mean of T, the mean time to failure of the component, is 1 l . Example 1: Let's consider the example presented in the introduction. Let the random variable T be defined as follows:
T=timeuntilarandomlyselectedbatteryfailsSuppose radio batteries havea timetofailure that is exponentially distributedwithameanof3ohours.Inthiscase,wecouldwrite1T ~EXP =301Therefore,perhour,sothat301α930f(t) :t > 0 and F(t)dx-3030F(t),theCDFoftheexponential distribution,canbe integratedto obtain-1F(t) =1-e 30.t>oNowwecancomputethereliabilityfunctionforabattery:R(t) = 1- F(t) = 1-e30t>o1-Recall that intheearlierexample,thesoldiersmustoccupytheOPfor24hoursThereliabilityofthebatteryfor24hoursisR(24) = e%(24) =0.4493,so the probability that the battery lasts more than 24 hours is 0.4493.Example2:Wehavetheoptiontopurchaseanewnickelcadmiumbatteryforouroperation.Testinghas shownthatthedistribution ofthetimetofailurecanbe modeled using a parabolic function:XX0≤x≤4848384f(x):[0otherwiseLettherandomvariableTbedefined asfollows:155
155 T = time until a randomly selected battery fails. Suppose radio batteries have a time to failure that is exponentially distributed with a mean of 30 hours. In this case, we could write T ~ EXP = 1 30 l æ è ç ö ø ÷ . Therefore, l = 1 30 per hour, so that f(t) e t = - 1 30 1 30 , t > 0 and F(t) e dx x t = - ò 1 30 1 30 0 . F(t), the CDF of the exponential distribution, can be integrated to obtain F(t) e t = - - 1 1 30 , t > 0. Now we can compute the reliability function for a battery: R(t) F(t) e e t t = - = - - æ è ç ö ø ÷ = - - 1 1 1 1 30 1 30 , t > 0 . Recall that in the earlier example, the soldiers must occupy the OP for 24 hours. The reliability of the battery for 24 hours is ( ) ( ) R 24 e 0 4493 1 30 24 = = - . , so the probability that the battery lasts more than 24 hours is 0.4493. Example 2: We have the option to purchase a new nickel cadmium battery for our operation. Testing has shown that the distribution of the time to failure can be modeled using a parabolic function: f x x x x ( ) = æ è ç ö ø ÷ - æ è ç ö ø ÷ £ £ ì í ï ï î ï ï 384 1 48 0 48 0 otherwise Let the random variable T be defined as follows:
T=time until a randomly selected batteryfailsInthiscase,wecouldwritef(t)=(t / 384)(1-t / 48), 0 ≤t ≤ 48 and F(t)= [(x/384)(1-x/48)dxRecall that intheearlierexamplethesoldiersmustmantheOPfor24hoursThereliabilityofthebatteryfor24hoursistherefore24R(24) = 1- F(24) = 1- [(t / 384)(1-t / 48)dx = 0.500Cwhichisanimprovementoverthebatteriesfromexample1ModelingSeriesSystemsNowweconsiderisasystemwithncomponentsC,,C,...,C,whereeachoftheindividual componentsmustwork inorderforthesystemtofunction.Amodelofthis type of system is shown in figure 1.ccc...Figure 1:Series systemIf we assume these components are mutually independent, the reliability of thistype of system iseasyto compute.Wedenotethe reliability of componentiattime t by R,(t). In other words, R,(t) is simply the probability that component iwill function continuously fromtime O through until time t.We are interested inthe reliability of the entire system of n components, but since these componentsaremutuallyindependent,thesystemreliabilityisR(t)= R,(t).R(t). ...-R,(t) .Example 3: Our radio has several components. Let us assume that there arefourmajor components--they are (in order)thehandset,the battery,thereceiver-transmitter,andthe antenna.Sincetheyall mustfunction properlyfortheradioto operate, we can model the radio with the diagram shown in figure2.156
156 T = time until a randomly selected battery fails. In this case, we could write f (t) = (t / 384)(1- t / 48), 0 £ t £ 48 and F(t) ( x )( x )dx t ò = - 0 / 384 1 / 48 . Recall that in the earlier example the soldiers must man the OP for 24 hours. The reliability of the battery for 24 hours is therefore R(24) 1 F(24) 1 (t 384)(1 t 48)dx 0 500 0 24 = - = - - = ò / / . which is an improvement over the batteries from example 1. Modeling Series Systems Now we consider is a system with n components C1 C2 Cn , ,., where each of the individual components must work in order for the system to function. A model of this type of system is shown in figure 1. C . 1 C2 C3 Cn Figure 1: Series system If we assume these components are mutually independent, the reliability of this type of system is easy to compute. We denote the reliability of component i at time t by Ri (t). In other words, Ri (t) is simply the probability that component i will function continuously from time 0 through until time t. We are interested in the reliability of the entire system of n components, but since these components are mutually independent, the system reliability is R(t) = R1 (t)×R2 (t)× ××× ×Rn (t) . Example 3: Our radio has several components. Let us assume that there are four major components - they are (in order) the handset, the battery, the receiver-transmitter, and the antenna. Since they all must function properly for the radio to operate, we can model the radio with the diagram shown in figure 2
CCFigure2:RadioSystemSuppose we know that the probability that the handset will work for at least 24hoursis0.6703,andthereliabilitiesfortheothercomponentsare0.44930.7261,and0.9531,respectively.Ifweassumethatthecomponentsworkindependentlyofeachother,thentheprobabilitythat theentiresystemworksfor24 hours is:R(24) = R;(24)·R2(24)·R(24)·R4(24) = (.6703)(.4493)(.7261)(.9531) = 0.2084Recall that two events A and B are independent if P(A|B)= P(A)ModelingParallelSystems(TwoComponents)Nowweconsiderasystemwithtwocomponentswhereonlyoneofthecomponentsmustworkforthesystemtofunction.Asystemofthistypeisdepicted in figure 3.Figure3:ParallelSystemofTwoComponentsNotice that in this situation the two components are both put in operation at timeO;theyareboth subjecttofailurethroughout theperiod of interest.Onlywhenbothcomponentsfail beforetimetdoesthesystemfail.Againwealsoassumethat thecomponentsare independent.Thereliability of thistype of system canbe found using the following well known model:P(AU B) = P(A) + P(B)- P(An B)Inthis case,Aistheeventthatthefirst componentfunctionsforlongerthansometime,t,andBistheeventthatthesecond componentfunctions longerthanthe sametime,t.Since reliabilities are probabilities,we can translate theaboveformula into the following:R(t) = R,(t)+R,(t)- R,(t)R2(t)157
157 C1 C2 C3 C4 Figure 2: Radio System Suppose we know that the probability that the handset will work for at least 24 hours is 0.6703, and the reliabilities for the other components are 0.4493, 0.7261, and 0.9531, respectively. If we assume that the components work independently of each other, then the probability that the entire system works for 24 hours is: R(24) = R1 ( 24)×R2 (24)×R3 (24)×R4 (24) = (.6703)( .4493)(.7261)( .9531) = 0.2084 . Recall that two events A and B are independent if P( A|B) = P(A). Modeling Parallel Systems (Two Components) Now we consider a system with two components where only one of the components must work for the system to function. A system of this type is depicted in figure 3. C1 C2 Figure 3: Parallel System of Two Components Notice that in this situation the two components are both put in operation at time 0; they are both subject to failure throughout the period of interest. Only when both components fail before time t does the system fail. Again we also assume that the components are independent. The reliability of this type of system can be found using the following well known model: P( AÈ B) = P(A) + P(B) - P( AÇ B). In this case, A is the event that the first component functions for longer than some time, t, and B is the event that the second component functions longer than the same time, t. Since reliabilities are probabilities, we can translate the above formula into the following: R(t) = R1 (t) + R2 (t) - R1 (t)R2 (t)