Evaluation of Uncertainties in the Design Process of Complex Mechanical Systems

In this paper, the problem of the evaluation of the uncertainties that originate in the complex design process of a new system is analyzed, paying particular attention to multibody mechanical systems. To this end, the Wiener-Shannon’s axioms are extended to non-probabilistic events and a theory of information for non-repetitive events is used as a measure of the reliability of data. The selection of the solutions consistent with the values of the design constraints is performed by analyzing the complexity of the relation matrix and using the idea of information in the metric space. Comparing the alternatives in terms of the amount of entropy resulting from the various distribution, this method is capable of finding the optimal solution that can be obtained with the available resources. In the paper, the algorithmic steps of the proposed method are discussed and an illustrative numerical example is provided.


Introduction
In industrial applications, the final design solution is almost always an engineering approximation, thus intrinsically prone to uncertainties.A viable approach for the analysis of the approximation in a design solution is based on the use of the non-probabilistic entropy.In the classic theories of Fisher and of Wiener-Shannon, the idea of information is a measure of probabilistic and repetitive events only.However, the concept is applicable beyond the probability field.In this respect, the Wiener-Shannon's axioms can be extended to non-probabilistic events so much so that it is possible to introduce a theory of information for non-repetitive events as a measure of the reliability of data.To this end, solutions consistent with the values of the design constraints may be selected by analyzing the complexity of the relation matrix and using the idea of information in the metric space.The final solution gives an entropic measure of epistemic uncertainties which can be then used in multibody system models analyzed by means of the axiomatic design.Multibody systems are the class of mechanical systems composed of rigid and/or flexible bodies interconnected by kinematic joints [1][2][3][4][5][6][7][8][9][10].They are subjected to nonlinear external force fields and undergo large, finite translations and rotations [11][12][13][14][15]. Examples of multibody systems can be found in several engineering applications [16][17][18][19][20].The complexity of the dynamic behavior of this family of constrained mechanical systems requires the development of advanced analysis and modeling tools to perform virtual prototyping in the multibody framework [21][22][23][24][25][26][27].The problem is particularly challenging when control strategies for flexible multibody systems need to be designed [28][29][30][31][32][33].In general, the design process gives the structure necessary to transform qualitative needs, often stated in non-engineering terms, into real products.In computer simulation of real problems, multiplicity of a data source may lead to a high level of error and, consequently, the final result can be considered only partially satisfactory.This process is carried out through the application of scientific knowledge to the problems.Based on a database of previous designs, the current design process generates several alternatives to be frequently evaluated.Usually, the design process is subdivided into phases with given specifications.Evaluations are required to determine whether the designer needs to go back to one or more phase.In the development of new industrial technology, the designer has the task of selecting among alternative solutions for complex systems.In this general process, besides a set of objective technical requirements, the designer must take into account some subjective factors that only he is able to quantify.The designer must also assess to what extent the selected alternatives fulfill the quantified targets.In particular, in the study of the dynamic behavior of multibody mechanical systems, it is necessary to develop a mathematical model with an appropriate objective function, which expresses the values of the combinations of possible outcomes as a single measure in terms of both cost and effectiveness.From a system design point of view, each subsystem is a system in its own right, i.e., it requires the fulfillment of a set of goals, organized following a hierarchic pattern, in which a single weight is assigned to each target according to the design relevance.In most common situations, a certain number of resources is not correctly allocated because the necessary values are not known or are estimated only in intervals.Consequently, imprecise data and imperfect knowledge may lead to uncertainties, and, therefore, there is the mathematical problem of the division of a quantity Q in a set of n fair values.If the data on distribution are not well-known, it is possible to use the principle of maximum entropy to achieve a consistent solution to the problem.

Complexity
In design research, the importance of explicitly considering uncertainty in the formulation of design is of fundamental importance for producing acceptable results.Multiple alternatives can yield uncertainties because of the lack of information.In this process, a good model helps to decide which alternative to take in order to achieve a final result that is acceptable.The best tools for obtaining good results are those that give a low quantity of uncertainties.In general, a set of interconnected components, interacting among each other in an organized way to achieve the same objective, constitute a system in a broad sense.Multibody mechanical systems represent a classical example of this general idea.Each system can be part of a super-system, i.e., a set of interacting systems.Therefore, to obtain a fast and good design necessitates analyzing the system complexity.Some projects may not need multiple sequences, while a large design may need to be analyzed in more layers.Each project may have its own organizational structure and subsystems.Complex systems are formed by subsystems which must operate in a coordinated way to achieve the desired goals.From a systems design point of view, each subsystem is a system in its own right and, therefore, each subsystem needs a set of requirements necessary to achieve the desired goals.Each requirement can be organized following a hierarchical scale, and its relevance is set by the designers assigning a weight in the hierarchy.The assignment of the weight to each element or system in the super-system is critical to establish a fair allocation of resources.An estimation of consistency of the individual resources must be conducted in order to give an overall assessment of the distribution.

Uncertainty
Uncertainty is a constant in engineering design.This is particularly true for multibody mechanical systems.Actually, conventional design methods suffer from the presence of uncertainty, so much so that new methods and new strategies have been recently introduced.By means of extensive theoretical research and virtual prototyping, a significant improvement can be obtained in the results within the framework of the new methods.The design process carried out under uncertainty is a difficult task because of both the lack of data and, more importantly, the lack of knowledge of the system under analysis.The effects of errors in input data and on the uncertainty of the model may be distinguished and quantified by the Bayesian inference.Choices made during the design with the Bayesian models are based on the theory that allows one to make optimal decisions under uncertainty.The goal of this design process is to maximize the expected utility of the project results, the being latter defined in terms related to the optimal experiment design, depending on the particular objective of the chosen product [34][35][36].

Fair Distribution of Resources
In most common situations, the belief is of no crisp because either its value is unknown or it is only known its range of existence [37].In mathematical terms, the quantity Q has to be divided in a set {q i } of an fair values with every element q i ≥ 0. If the distribution is fair, with v i ≥ 0, which is divided by the weight of the nth share, one can write: and When the information on distribution is not completely well-known, it is possible to use the maximum entropy principle to obtain a solution of the problem.The model of the maximum entropy principle is inborn in Western philosophy and is the mathematical transcription of the ancient principle of wisdom: in the absence of information on events, the same belief holds for any of them.This statement is known as the Principle of Insufficient Reason of Laplace and Bernoulli: if there is no reason to believe that out of a set of possible, mutually exclusive events, no one event is more likely to occur than any other, then one should assume that all events are equally probable.This is a decision principle because, in the absence of a priori reason, one can assign an equal degree of belief to all the events.In presence of incomplete information, the maximum entropy principle assigns and completes the information.Using Shannon's definition, the entropy in a uniform (equally likely) distribution of probability assumes the maximum value.Therefore, one can write: On the basis of the Laplace-Bernoulli principle, one can assign degrees of belief to finite collections of events.In a probability finite space Ω the measure of probability must observe Kolmogorov's axioms, which can be summarized as follows: If in the probability space (Ω, A, P) Ω is a finite non-empty set; A is an σ − algebra off the subset in Ω of the events; P is a probability measure on A, Then the applications of P : A → [0, 1] must be under axioms: Now, if we analyze our problem in Kolmogorov's mathematical formalism [38], we can use his philosophy of probability and how he proposed to relate the mathematical formalism to the real world.With this formalism, we can use the tool of Shannon's entropy to measure the uncertainty of the division.
When one knows only that the belief is included in intervals, an algorithm is required for the division of a quantity.Indicating with n the number of intervals, a set of {v i } fair values under the condition n ∑ i v i = 1 can be found.For the generic v i ≥ 0, one has v i ∈ v i , v i , where v i is a lower value and v i is the upper value of the interval of v i , and therefore: The value v i is the lower possible of its interval value of uncertainty and must be distributed.The fair partition depends on the weight of intervals d i .
The weight W of all intervals is n ∑ The values of v i are given by: In order to have consistent problems, this formula can be used for the distribution of a quantity included in intervals in an uncertain position.If one has a set of n = (m + k + u) elements with k values well-known, m values included in intervals in uncertain position, and u unknown values, in the distribution with i ∈ [1, n] and one cannot assign the Laplace-Bernoulli values because all of the information is not known.Only the values of v k are well-known.The values (k + m) can be divided into two subsets, S = S ∪ S, with the subset S indicating the lower values and the subset S denoting the upper values: One thus obtains the data, ∆ and ∆, for the distribution from the subset S and S, while the sub-rules are: b.
If ∆ > 0 then the distribution is consistent.c.
If ∆ > 0 and ∆ > 0 then the value of the upper set ∆ determines the lower distributions and the one of the lower set ∆ determines the higher distribution (Table 1).

Well-Known Values Values Included in Intervals
Values of the Laplace-Bernuille Distribution Using the principle of maximum entropy, the fair distribution can be determined.In this paper, the use of the method by a very simple application reported in Table 2 is introduced.
The proposed method is new.It can be used in the distribution of a quantity in the presence of uncertainty in probabilistic and non-probabilistic problems.

Measure of the Information of a Distribution of Imprecise Data
In the theories of Fisher and Wiener-Shannon, the idea of information is not only a measure only of probabilistic and repetitiveness events, but can also be used as a measure of the fairness of the complex distributions [39,40].The use of the information axioms allows for the development of useful and applicable information models.For every event t is possible to have a measure of information using the mathematical expression: This definition of information has a natural application in the metric space.Norbert Wiener in Cybernetics gave the clearest definition of entropy: The information can be evaluated by the probability and by the non-probabilistic measures of diameters.Thus, it is possible to have the measure of the information from non-probabilistic data.
For the non-probabilistic information, instead of probability, it is possible to utilize a finite number of appropriate proportions, subject to a set of constraints that add up to 1. Let d 1 , d 2 , . . ., d n be non-negative real numbers, and let: One can use the following relation as a measure of information: Thus, J(ρ) is maximum when ρ 1 = ρ 2 = . . .= ρ n , and J(ρ) is minimum when, for every element i, only one number is different from zero.In the final result, there is the final distribution that also takes into account the uncertainty of data, with a fair distribution of imprecise data and a global measure of uncertainty.For example, one can write: If the final distribution is uniform, then the result would be a maximum number of bits because the uniform distribution involves the maximum entropy value: Thus, it is possible to measure the entropy of the two distributions and to choose the best solution for our design.

Conclusions
The principal goal of the authors' research is the development of new methods for the design and the analysis of complex mechanical systems [41][42][43][44][45][46][47][48][49][50][51].In general, the designer of complex systems is often called upon to select a suitable solution among several alternatives.For industrial engineering applications, this is particularly true in the design process of multibody mechanical systems composed of several bodies constrained by kinematic joints.In this challenging process, some subjective factors that are difficult to quantify analytically must be taken into account in order to judge which alternative will effectively fulfill the desired targets.When possible, in this complex process a mathematical model which includes an objective function can be developed order to facilitate the design process using a single quantity expressing the performance of the combinations of all the possible outcomes in terms of cost and effectiveness.From a system design point of view, each subsystem is a system in its own right, i.e., it must achieve a set of goals organized in a hierarchical way by assigning a weight according to their relevance.The assignment of the weight for each element or each system in the super-system can be obtained using the fair distribution model.Using this methodology, the challenge is to establish the right distribution for a given problem.It is thus necessary to devise some mathematical tools to compare and evaluate two or more distributions.For this purpose, Wiener-Shannon's axioms can be extended to non-probabilistic events, so that it is possible to introduce a theory of information for non-repetitive events as a measure of the reliability of data.To this end, one can select solutions consistent with the values of the design constraints by analyzing the complexity of the relation matrix and using the idea of information in the metric space.Analyzing the alternatives by the amount of entropy resulting from the various distribution, this method allows one to find the optimal solution that can be obtained with the available resources.As shown in the paper, the principal algorithmic steps for a correct analysis are three, namely, the creation of the alternative solutions to meet the functional requirements of the project, the evaluation of alternatives based on the resource budget, and the entropic classification of the alternative solutions.In the paper, a simple numerical example is illustrated in order to demonstrate the effectiveness of the proposed methodology.