[Editor's Note: This article has been updated since its original publication to reflect a more recent version of the software interface.] The demand for products with high reliability and low manufacturing costs is an ever increasing one in recent times. Inception of a product that beats competition in terms of both reliability and cost requires a design strategy that steers the engineering design process toward higher reliability from the very beginning. Engineering Design By Reliability (EDBR) is a methodology used to produce components that meet a certain reliability goal by designing reliability directly into the components. In this methodology, all design parameters are taken to be random variables. All parameters that define the failure governing stress acting on a component during its mission are taken to be distributed variables and all parameters defining the failure governing strength exhibited by the component during its mission are also taken as distributed. As a result, the failure governing stress and strength are obtained as distributions using various techniques of synthesizing these distributions. These are then coupled mathematically to obtain the reliability of the component, which is then compared to the target reliability. If the obtained reliability is equal to or greater than the target reliability then the design is accepted; otherwise, the design is iterated until the specified reliability requirement is met. Advantages of EDBR
Over Conventional Design Methods Obtaining the
Distribution of a Function of Random Variables The binary synthesis of distributions method assumes a normal distribution for all random variables of the desired function and also assumes a normal distribution for the resulting function distribution. The method uses the mean and standard deviation of all the random variables in a function to arrive at the mean and standard deviation of the function. This method may not provide accurate results for cases where a normal distribution assumption does not apply. The generation of system moments method involves obtaining the first four moments of the function whose distribution is sought. Then the Pearson distribution approximation is used to arrive at the distribution. The Monte Carlo simulation method involves the use of simulation to arrive at the desired distribution of the function of random variables. The accuracy of results obtained here depends largely on the number of simulations used. Nonrepeatability of results is another minor drawback of this method. However, with increased complexity of problems and availability of software tools, simulation is emerging as the method of choice for many analysts. Failure Governing
Stress Distribution It must be realized that distributions of loads, geometric dimensions and any modifying factors used to obtain the stress distribution are componentspecific and have to be determined through observation and experimentation as most of the time general use data for these may be either not available or not applicable. The distribution of loads, moments and torques acting on a component may be determined by attaching transducers to critical locations of the component where stress is expected to be the highest. Distributions of geometric dimensions may be determined by measurements on a sample of components and fitting the most appropriate distribution to the recorded observations. If no such observations are obtainable but a tolerance value is specified and a normal distribution can be assumed for the dimension of the component, then the base dimension can be taken as the mean of the normal distribution, and a sixth of the tolerance as the standard deviation, because the total tolerance is taken to equal six standard deviations. Multiplication factors used should be studied separately to obtain their distributions. Failure Governing
Strength Distribution Estimation of
Component Reliability
Where:
f*(s) is called the failure function. It can be seen that Q is the cumulative function of f*(s). However f*(s) is not a probability density function but simply a failure function and the area under f*(s) is generally less than 1. It should be noted that the unreliability of the component is not given by all of the area of interference of the stress and strength distributions. The area representing the unreliability may be a part of the area of the interference and at times may lie outside the overlap of the two distributions. For instance, if a stress distribution is normal with a mean value of 50 units and standard deviation of 15 units, while the corresponding failure governing strength distribution is normal with a mean of 75 units and standard deviation of 5 units, then the area bound by the failure function representing the unreliability of the component is as shown in Figure 1.
Figure 1: StressStrength distribution interference The reliability R of the component can be obtained as the difference between 1 and the unreliability value obtained from the above expression of Q. Alternatively, the reliability of the component can be calculated as the probability that strength S is greater than all possible values of stress s, which in mathematical terms is expressed as:
Where:
f' *(s) is called the survival function. Once the reliability of the component is known, it can be compared to the target reliability. If the reliability value equals or exceeds the reliability target then the component design is a successful one. If the value of R obtained above is lower than the target value then the design of the component is modified and the design process is repeated after incorporating changes in the design parameters that lead to a higher probability of mission success for the component. The procedure is iterated until the reliability goal is met. In the remainder of this article, a simple example illustrates the EDBR methodology using ReliaSoft’s Weibull++ software Example
Table 1: Typical load values (in lbs) for the application
Table 2: Ultimate tensile strength data (in psi) of 50 specimens of AISI 4340 steel rods. Ref [1]
Stress Formulation Stress formulation in the EDBR methodology is the same as in conventional design methods. For the present example, stress s is simply:
where F represents the tensile load on the rod and d represents the rod's diameter. Both of these variables are distributed. The appropriate distributions and the respective parameters of the distributions for both variables need to be determined. Determination of the Stress Distribution Assuming that the diameters of the rods are normally distributed, the parameters of the corresponding normal distribution can be determined from the given base dimension and tolerance value. Thus the mean of the assumed normal distribution, which is taken to be the base dimension, is 0.5. The standard deviation, which is taken to be a sixth of the total tolerance value, is 0.005. The distribution for the rod diameter d can be represented as N(0.5, 0.005). The available load data can be used to determine the distribution of F using probability plotting, least squares or maximum likelihood estimation. The data can be sent to the Weibull++ software and the Distribution Wizard can be used to decide on the best distribution and calculate the parameters. Weibull++ gives the twoparameter Weibull distribution as the best fit distribution to the present load data. The parameters of this distribution are calculated as beta = 3.3435 and eta = 14278. Thus the distribution of load F can be represented as W(3.3435, 14278). Once the distributions of all the design variables in the failure governing stress formulation are known, the distribution of the failure governing stress s can be obtained using any of the three methods of synthesizing distributions discussed previously. The Monte Carlo simulation method is illustrated here. Monte Carlo Simulation Weibull++ can be used to generate random numbers based on the distribution of the diameter d and the load F. These random numbers can then be used to calculate random stress values using equation Eqn.(5). A distribution can be fitted to these stress values to obtain the failure governing stress distribution. For this example, a set of 1000 random numbers are generated based on the distribution N(0.5, 0.005) of diameter d. A set of 1000 random numbers are also generated based on the load F distribution W(3.343451,14278.157884). These random numbers from the d and F distributions are substituted into Eqn. (5) to obtain a set of 1000 stress values. Figure 2 shows the generated diameter and load values, as well the stress calculation. A distribution is then fitted to the calculated stress values to obtain the failure governing stress distribution. The generalized gamma distribution G(11.195576, 0.262346, 1.063886) is obtained as the best fit distribution. Figure 2: The general spreadsheet (shown in background) with the values obtained from the Monte Carlo utility and resulting stress calculation. The stress values are then transferred to a data folio and analyzed (foreground). Determination of the Strength Distribution A distribution is fitted to the ultimate tensile strength data of Table 2. As expected, the normal distribution gives the best fit to these strength data. The probability plot for the data is shown in Figure 3. The fitted failure governing strength distribution is obtained as N(103421.079990, 2395.106091).
Figure 3: Normal probability plot of the ultimate tensile strength data of Table 2 Calculation of Reliability The Weibull++ stressstrength folio is used to estimate the reliability of the steel rods once the failure governing stress distribution G(11.195576, 0.262346, 1.063886) and the failure governing strength distribution N(103421.079990, 2395.106091) have been obtained. The reliability value is obtained as 97.92%. Figure 4 shows the stressstrength plot of the distributions. It is clear that the desired goal of no more than 10% failures is met by the present design.
Figure 4: StressStrength plot
Alternative Approach Using Monte Carlo Simulation
After a number of random values representing the difference between strength and stress are obtained, a normal distribution is fitted to these points and the desired reliability is obtained as the probability that any of these values would be greater than 0. This approach is demonstrated in Figures 5 and 6. The result, 96.95%, is comparable to the one obtained with the other method.
Figure 5: The userdefined function in the Monte Carlo utility is used to generate values that represent the difference between S and s.
Figure 6: A normal distribution is fitted to the generated values and the QCP is used to calculate the probability that the generated S  s values 0. Conclusion
References
