Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Fractional Series Solution Construction for Nonlinear Fractional Reaction-Diffusion Brusselator Model Utilizing Laplace Residual Power Series
Previous Article in Journal
Transmission Properties of Electromagnetic Waves in Magneto-Electro-Elastic Piezoelectric Electromagnetic Metamaterials
Previous Article in Special Issue
A Hybrid Particle Swarm Optimization-Based Wavelet Threshold Denoising Algorithm for Acoustic Emission Signals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Symmetry Importance in a Relative Entropy Analysis for Some Engineering Problems

Head of Civil Engineering and Transportation, Department of Structural Mechanics, Faculty of Civil Engineering, Architecture and Environmental Engineering, Łódź University of Technology, 90924 Lodz, Poland
Symmetry 2022, 14(9), 1945; https://doi.org/10.3390/sym14091945
Submission received: 29 July 2022 / Revised: 5 September 2022 / Accepted: 13 September 2022 / Published: 18 September 2022

Abstract

:
This paper aims at certain theoretical studies and additional computational analysis on symmetry and its lack in Kullback-Leibler and Jeffreys probabilistic divergences related to some engineering applications. As it is known, the Kullback-Leibler distance in between two different uncertainty sources exhibits a lack of symmetry, while the Jeffreys model represents its symmetrization. The basic probabilistic computational implementation has been delivered in the computer algebra system MAPLE 2019®, whereas engineering illustrations have been prepared with the use of the Finite Element Method systems Autodesk ROBOT® & ABAQUS®. Determination of the first two probabilistic moments fundamental in the calculation of both relative entropies has been made (i) analytically, using a semi-analytical approach (based upon the series of the FEM experiments), and (ii) the iterative generalized stochastic perturbation technique, where some reference solutions have been delivered using (iii) Monte-Carlo simulation. Numerical analysis proves the fundamental role of computer algebra systems in probabilistic entropy determination and shows remarkable differences obtained with the two aforementioned relative entropy models, which, in some specific cases, may be neglected. As it is demonstrated in this work, a lack of symmetry in probabilistic divergence may have a decisive role in engineering reliability, where extreme and admissible responses cannot be simply replaced with each other in any case.

1. Introduction

Uncertainty analysis is one of the most important aspects of a solution to real engineering problems [1], in which geometrical, and material imperfections as well as a statistical scattering of the environmental actions may play a decisive role. Such an analysis may concern the theoretical or numerical determination of probabilistic characteristics of the structural response (expectations, variances, etc.) or finding a difference (distance) of maximum (extreme) responses to their admissible counterparts. This second problem is fundamental for reliability assessment [2], in which engineering structures or systems are verified in the context of deformations, vibrations, stresses, temperatures, or also other physical quantities. Both extreme responses and also their admissible counterparts may have an uncertain nature so that a realistic mathematical problem would be to find the probabilistic distance between two different probability density functions, where one is discrete, and the opposite–most frequently–has continuous character. Therefore, it would be quite natural to apply the probabilistic divergence [3] (relative entropy) to measure the safety margin of both existing as well as newly designed structures and systems. This idea has been explored yet in the literature a few times and may also concern some other research issues like detecting a difference in-between the experimental and numerical results, which is a subject of the specific calibration procedures in deterministic analyses. Some other applications may be seen in eigenfrequency and the induced frequency in forced vibrations of the given system, the extreme and phase change temperatures in the given material, etc. One of the very specific ideas is looking for a distance between original material characteristics and their effective counterparts resulting from some homogenization procedures.
Relative entropy [4] (probabilistic distance) applicable in this context may be defined for univariate and multivariate random sources but always demands a precise definition of two probability distributions (or at least their histograms biased with relatively small numerical errors). Most mathematical models in this area show full symmetry, which means that a choice of primary and target distribution makes no difference. The only exception is the so-called Kullback-Leibler distance [3], which is affected by this choice and which exhibits a certain lack of symmetry and which could be very sensitive to contrasted distribution types and parameters. Mathematical proofs and considerations have been demonstrated in the literature; nevertheless, the importance of real engineering systems and their reliability has not been precisely studied. Such a lack of symmetry may result in some underestimation (or overestimation) of structural safety, which leads to non-optimal design decisions or even to designing over the realistic limits and demands; such an effect is an unwanted situation in both cases.
This asymmetry is the main research problem studied in this work, in which a two-fold methodology has been proposed. The first approach is based upon symbolic integration inherent in the Kullback-Leibler entropy (KL) definition assuming for a brevity of presentation two different Gaussian probability distributions with given intervals of their parameter variations. It is very rare in real engineering problems that the given system response may be described with a specific probability density function and its parameters. That is why the second method of the KL entropy determination is the most common application of the Monte-Carlo simulation as well as the generalized stochastic perturbation technique together with the Finite Element Method (FEM) [5] program. Such a two-fold methodology guarantees a huge variety of possible applications, especially because all probabilistic procedures employed in this study have been programmed by the author in the numerical environment of the system MAPLE. The few computational experiments attached here confirm an importance of the KL entropy asymmetry in both engineering and applied sciences problems while studying a distance between extreme and admissible structural response functions.

2. Historical and Mathematical Background

Probabilistic divergence (also known as the relative entropy [6]) has directly followed the foundations of thermodynamics, where Gibbs (1878) entropy of a given thermodynamical system has been introduced in the statistical context as
S = k B i p i ln ( p i )
where k B is the Boltzmann constant (1872), while p i denotes the probabilities of all the microstates. Johann von Neumann (1927) modified this definition to the following form:
S = k B T r ( ρ   ln ( ρ ) )
with ρ being the density matrix of the given quantum mechanical system. Finally, Shannon (1948) [7] proposed this entropy as
H = i p i ln ( p i )
which, with some small modifications and extensions toward continuous probability distributions, has been used until now. Lower index i counts here and further the given populations of random accidents. Quite independently, probabilistic divergence itself has been invented and developed more than a hundred years ago (Hellinger, 1905 [8]) to measure a distance in-between two different probability distributions. It has received a lot of attention in mathematical studies, but even more in applied sciences and engineering applications. Some of them concern quantum mechanics, fracture analysis and some other issues exhibiting statistical scattering. One of the most famous models is the so-called Kullback-Leibler (KL) model [9], whose principal deficiency is a lack of symmetry while analyzing a divergence of two given distributions p and q according to its discrete version.
D K L ( p , q ) = i p i log ( p i q i )
In general cases, such as D K L ( p , q ) D K L ( q , p ) , Jeffreys [10] have proposed a symmetrization of this divergence by simply adding
D J ( p , q ) = D K L ( p , q ) + D K L ( q , p )
A lack of symmetry in the KL divergence is of paramount importance while analyzing engineering reliability problems, where p may denote structural effort and q means structural strength, for instance. Now, the area of relative entropy includes more advanced models adjacent to both univariate and multivariate PDFs as Jensen-Shannon [11,12], Rényi [13], Bhattacharyya [14], Mahalanobis [15], Tsallis [16], and others have proposed. Nevertheless, symmetry, its possible lack, and its consequences for the probabilistic divergence remain one of the most important theoretical and numerical topics.
It is well known that KL divergence has a huge number of well-documented applications in various branches of engineering starting from energy analysis [17,18], marine engineering [19], mechanical problems of failure [20,21,22] as well as in civil engineering analyses [23]. Some recent mathematical works concern the maximum entropy principle [24] or multivariate distributions [25], while physicists are interested in i.e., non-equilibrium random motion [26]. A comparison of various probabilistic divergence models is also available in the literature and is documented for Hellinger and Kullback-Leibler theories [27], Jensen-Shannon and KL distances [28] as well as for a relation of the latter model of Tsallis statistics [29].
The main aim of this work is to compare two closely related probabilistic relative entropies proposed by Kullback-Leibler and Jeffreys [30], and this is done to study an impact of a lack of symmetry in the first theory. The numerical study proposed here concerns three different examples, and it starts from the analytical symbolic determination of both divergences in the case of two different combinations of neighboring Gaussian distributions. The first computational experiments is based upon symbolic integration inherent in the definitions of two relative entropies under consideration. Further, a reliability study relevant to some steel Pratt trusses is presented, in which KL and Jeffreys divergence is used to detect a distance in between the PDFs of extreme and admissible structural deformations. The extreme deformations are determined using the stochastic finite element method approach realized with (i) probabilistic extension of traditional deterministic designing procedures, (ii) the Monte-Carlo simulation analysis and (iii) the iterative generalized stochastic perturbation method, while the admissible deformations are given due to the engineering codes.
The final example deals with the homogenization procedure of some particulate composite, in which divergence between real and homogenized composite characteristics is under investigation. The computational analysis presented in this paper exhibits that a lack of symmetry in the KL divergence may play a remarkable role in some engineering problems. As has been demonstrated here, it is highly important whether the Kullback-Leibler distance is measured in between the admissible and extreme deformations or vice versa. This aspect is completely absent in reliability analysis, for instance, where an absolute value of the difference of admissible and extreme structural behavior mean values (scaled by its variance) is analyzed. A more remarkable effect is seen while approximating the analogous distance between the original and homogenized composite materials characteristics. Some further idea could be a discussion of probabilistic distance between the designed and real deformations or stresses in the erected civil engineering structures or mechanical systems. This aspect has not been discussed until now in the literature, but it may have some importance, especially in highly nonlinear problems with uncertainty when reliability assessment is decisive for structural safety and durability.

3. Relative Entropies and Probabilistic Numerical Methods

Relative entropy equations applied and discussed here follow probabilistic divergence models invented by Kullback and Leibler as well as by Jeffreys [30]. Their definitions for two different probability distributions p(x) and q(x) are formulated in turn as follows [9]:
D K L ( p , q ) = p ( x ) log ( q ( x ) ) d x + p ( x ) log ( p ( x ) ) d x
D J ( p , q ) = D K L ( p , q ) + D K L ( q , p )
As a specific case, one may consider two different Gaussian probability distributions: p ( x ) N ( μ 1 , σ 1 ) q ( x ) N ( μ 2 , σ 2 ) . Then, the KL relative entropy equals to
D K L ( p , q ) = log ( σ 2 σ 1 ) + σ 1 2 + ( μ 1 μ 2 ) 2 2 σ 2 2 1 2
whereas the Jeffreys entropy is calculated as
D J ( p , q ) = log ( σ 2 σ 1 ) + σ 1 2 + ( μ 1 μ 2 ) 2 2 σ 2 2 + log ( σ 1 σ 2 ) + σ 2 2 + ( μ 2 μ 1 ) 2 2 σ 1 2 1
Quite similarly, one may derive the relative entropies for two Weibull distributions having the parameters p ( x ) W ( k 1 , l 1 ) and q ( x ) W ( k 2 , l 2 )
D K L ( p , q ) = log ( k 1 l 1 k 1 ) log ( k 2 l 2 k 2 ) + ( k 1 k 2 ) ( log ( l 1 ) γ k 1 ) + ( l 1 l 2 ) k 2 Γ ( k 2 k 1 + 1 ) 1
As it is seen, reliable determination of the first two probabilistic moments is necessary to study the probabilistic divergence of two Gaussian distributions. This is the reason to engage not only the Monte-Carlo simulation (MCS) approach but also some faster techniques, such as the stochastic semi-analytical method (SAM) and the iterative stochastic perturbation method (SPT) [31,32,33]. These two last methodologies enable remarkable reduction of a solution of the problems with uncertainty, whose computational effort becomes closer to the deterministic origin rather than the MCS time consumption. A computational implementation may be delivered with the use of the MCS approach, in which statistical estimation usually follows the well-known statistical estimators introduced as
E [ y ( b ) ] = 1 M m = 1 M y m ( b ) = 1 M m = 1 M p = 0 P A p b m p
V a r ( y ( b ) ) = 1 M 1 m = 1 M { y m ( b ) E [ y ( b ) ] } 2 = 1 M 1 m = 1 M p = 0 P { A p b m p E [ y ( b ) ] } 2
where b is the input uncertainty source, bm means its mth random numerical realization, and y(b) is the given structural output. Alternatively, the iterative generalized stochastic perturbation technique is also employed. Its fundamental step is the following nth order Taylor expansion of the same response function y(b) [31,32]:
y ( b ) = y 0 ( b 0 ) + i = 1 n ε i i ! i y ( b ) b i | b = b 0 ( Δ b ) i
It concerns the function f of the given random input parameter b about its mean value b0 so that Δb is the first variation of this parameter about its given mean. This expansion is embedded next into the classical integral definition of the kth central probabilistic moment as
μ k ( y ( b ) ) = + ( y ( b ) E [ y ( b ) ] ) k   p b ( x ) d x = + { y 0 ( b 0 ) + i = 1 n ε i i ! i y ( b ) b i | b = b 0 ( Δ b ) i E [ y ( b ) ] } k p b ( x ) d x
Further transforms are possible after a choice of the specific PDF type as well as its upper and lower bounds. It should be mentioned here that the perturbation method is still frequently applied in many mathematical studies focused on dynamical systems and nonlinear phenomena [34,35].
The semi-analytical approach directly uses the definitions of the central probabilistic moments enriched with the polynomial bases as
μ k ( y ( b ) ) = + ( y ( b ) E [ y ( b ) ] ) k   p b ( x ) d x = + { i = 1 N A i ( b 0 ) i E [ i = 1 N A i ( b 0 ) i ] } k p b ( x ) d x
The purpose of the first numerical experiment was to investigate analytically probabilistic divergence for two different Gaussian distributions N 1 ( E 1 , σ 1 )   N 2 ( E 2 , σ 2 ) , and it was entirely done in the computer algebra system MAPLE 2019.2. This purpose has been achieved by introducing two additional parameters being a difference between the coefficients of variation of these two distributions, α 1 and α 2 , as well as between their expectations, accordingly
Δ ( α ) = α 2 α 1 = σ 2 E 2 ,   Δ ( E ) = E 2 E 1
Two different combinations of the input data have been adopted to check some initial parametric sensitivity of both Kullback-Leibler and Jeffreys relative entropies. The first test corresponds to the following data: E1 = 2.0 & α1 = 0.01, whereas the second one uses E1 = 1.0 and α1 = 0.01, respectively. The resulting entropy surfaces have been collected in Figure 1 (the first input combination) as well as in Figure 2 (the second test).
The first general observation is that the Jeffreys model returns a few times larger entropy than the Kullback-Leibler theory so a lack of symmetry matters a lot. The second model exhibits significant sensitivity to the input data, while in the Jeffreys model, this effect is almost negligible. The KL approach seems to be more sensitive to a difference in-between coefficients of variation of two random inputs and then, to their expectations. Symmetrized Jeffreys theory shows even higher changes of the entropy while modifying expectations difference, while sensitivity to statistical scattering of the input data is either very small or may be just postponed. Nevertheless, as one could expect, the larger the difference in input expectations, the higher the resulting relative entropy (a distance between both PDFs). It directly follows a deterministic case of two functions and their mean values. Quite a different observation can be made for two input coefficients of variation–a distance for the PDFs under consideration becomes smaller while increasing their difference. It is because the concentration of the PDF about its mean value is smaller while increasing this coefficient so that the PDF splashes and two different bell-shaped curves are effectively closer. Finally, one may conclude that symmetry (and its possible lack) in the relative entropy computations may play a very important role in applied sciences and cannot be simply neglected.

4. Reliability Structural Analysis and Probabilistic Divergence of the Elastic Pratt Truss Structure

4.1. Numerical Model

The next computational experiment has documented the usage of two aforementioned relative entropies to study probabilistic divergence between the extreme (E) and admissible (R) deformation of some civil engineering structures. This divergence is the basis of the so-called limit function g(R,E) = RE widely employed in the reliability theory [36]. A specific case is investigated, where both R and E have random nature. A simple plane steel Pratt truss with 12 segments is analyzed for this purpose (Figure 3), whose height has been adopted as h = 1.0 m, and its span l is treated as the Gaussian variable having the expectation equal to 18.00 m. The stochastic finite element method (SFEM) analysis of this structure in the context of Bhattacharyya entropy has been analyzed in [37]. The first two probabilistic moments of structural displacements have been determined here via the series of traditional FEM experiments [38], where polynomial bases connecting extreme deformations with the input uncertainty source have been recovered. All these experiments have been carried out here in the civil engineering system Autodesk ROBOT, while the probabilistic part has been completed in the computer algebra system MAPLE 2019.2. This truss FEM model has been uploaded with the constant load q redistributed throughout the upper chord nodal points as the concentrated forces. Its characteristic value has been set as qk = 15 kN/m, while the design load has been proposed for an illustration as qd = 20 kN/m; a dead load of the truss is taken into account separately. All the cross-sections of structural components have been assumed as hot-finished square-hollow steel profiles made of structural steel S235, whose material parameters have been assumed according to the design code Eurocode 3 [39]. The following profiles have been found here as the most optimal solution: (i) upper chord–SHS 140 × 140 × 8, (ii) diagonals–SHS 90 × 90 × 8, and also (iii) lower chord–SHS 120 × 120 × 8. The FEM discretization has been completed using 49 linear two-noded truss finite elements, where horizontal and vertical displacements have been fixed at lower left edge; a vertical displacement has been blocked at the right lower corner. A single finite element represents here a column and a diagonal, while upper and lower chords have been represented by the set of single finite elements from node to node, respectively.
Numerical experiments have been completed assuming geometrical nonlinearities and p-delta effect in the structural behavior. The incremental Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm has been employed for this purpose with the following parameters:
  • load increment number has been fixed as 5;
  • maximum iteration number for a single increment equals 40;
  • increment length reduction number was set as 3;
  • increment length reduction factor was equal to 0.5;
  • the maximum number of line searches equals 0;
  • control parameter for the line search method is set as 0.5;
  • the maximum number of the BFGS corrections is set as 10;
  • relative tolerance for the residual forces and displacements is taken as 0.0001.

4.2. Numerical Results and Discussion

The basic results of numerical studies have been collected in Table 1 below as the functions of an increasing coefficient of variation of the uncertain truss length ( α ( l ) [ 0.00 , 0.15 ] ). They are obtained with three different numerical methods described briefly in the preceding section and they include (i) the resulting coefficient of variation of the extreme vertical displacement α ( u ) obtained on the symmetry axis of the structure, (ii) Kullback-Leibler (DKL) relative entropy, and (iii) Jeffreys probabilistic divergence (DJ). First of all, it is seen that a lack of symmetry in probabilistic divergence measured here as a difference DJ−DKL naturally increases together with an increase of the input coefficient α. Nevertheless, it is rather not larger than 1% of the DJ value at the very upper end of the input α range, which is rather negligible.
The very important result is that all three stochastic numerical methods coincide here almost perfectly, which enables significant reduction of computational time by replacing the MCS with SPT or SAM techniques. Further, both entropies exponentially decrease together with a linear decrease of the input COV, which is quite an expected effect. As one may see, particular numerical values of these entropies cannot be directly used in structural reliability assessment as they are with a few ranges larger than the reliability indices listed in Eurocode 0. Let us mention that some rescaling procedure has been proposed for this purpose in [37].

5. Probabilistic Divergence Application to the Homogenization of some Particulate Composite

5.1. Composite Material Numerical Model

The last numerical case study was prepared to study a probabilistic divergence between the random distribution of some mechanical characteristics of the given composite material and its homogenized characteristics. This experiment has been focused on a particulate composite consisting of a single spherical particle centrally embedded into a cubic matrix volume, see Figure 4. A perfect interface between the particle and the matrix is assumed here, and initial deformation has been completely neglected. Its discretization [40] was made in the preprocessor of the ABAQUS system and is shown in Figure 5.
The homogenization procedure performed and analyzed here is a widely known mathematical-numerical approach invented for both classical composite materials as well as for some multiscale heterogeneous structures. The so-called effective characteristics for these complex materials/structures are determined via analytical calculus or in some FEM computer analyses. All such approximations are based upon an assumption that the deformation energy of the real and of the homogenized media are equal to each other. Closed-form equations leading to these effective characteristics are derived from the specific boundary conditions imposed on the composite representative volume element (RVE) [40].
A composite consisting of the rubber particle embedded into the polymeric matrix is analyzed, whose elastic parameters of the particle have been taken as Ep = 1.0 MPa, νp = 0.4888, while for the polymer matrix as equal to Em = 4.0 GPa and νm = 0.34. Both Young moduli of this composite have been adopted here in turn as Gaussian variables with given expectations and some variability interval for their coefficients of variation to see an influence of the input uncertainty ( α [ 0.00 , 0.20 ] ). Polynomial approximating functions relating the effective elasticity tensor components and Young moduli of the particle and of the matrix have been recovered numerically via certain series of the FEM experiments and the least squares method implemented in the computer algebra system MAPLE 2019.2 as before.

5.2. Numerical Results and Discussion

Numerical results of such analysis have been contained in Figure 6, Figure 7 and Figure 8–for C 1111 ( e f f ) , C 1122 ( e f f ) and C 1212 ( e f f ) while randomizing particle Young modulus (left graphs) and, separately, matrix Young modulus (the right column). Kullback-Leibler and Jeffreys probabilistic divergences have been compared here, and a general observation one can make is that a lack of symmetry is more remarkable in the right graphs, where the absolute extreme distance of the real and homogenized characteristics is a few orders smaller. The opposite case–particle uncertainty–shows practically the same results for two entropy models.
The next observation of a quite general character is that all relative entropies decay very fast while increasing the input coefficient of variation; they are close to 0 for about α 0.20 . The extreme values marked on the vertical axes allow us to compare the distance of the PDF for particle and for matrix Young moduli to the given effective tensor components. As one may expect, the elasticity modulus of the matrix is many times closer to any component of the tensor C i j k l ( e f f ) than the elastic modulus of the rubber particle. Looking for the expected values of the effective tensor and comparing them with Young moduli of both components, it is seen that particle contribution to the overall elastic characteristics is marginal in this context. It is further seen that a relative entropy could be used in composites engineering (not only in the context of the homogenization method) to find out the most distant random parameters from the given limit state, whose uncertainty may be simply postponed in the model. Other parameters having especially the smallest distance to that limit state may be decisive in reliability assessment, for instance, for the given composite material.
The last conclusion is that various relative entropies may result in totally different numerical values; however, a general character of the function H = H(α) remains the same (very similar to engineering reliability indices). The extreme values result in these graphs from the Jeffreys model, and then from the Kullback-Leibler model, entropy is almost the same or slightly smaller for the given COV value.

6. Concluding Remarks

Symmetrical and non-symmetrical probabilistic divergencies have been presented in this work and contrasted with each other to check the possible lack of symmetry in relative entropy and its possible significance in numerical simulation and specific engineering problems. It has been found that this lack of symmetry may have a remarkable influence on relative entropy values in some specific mathematical analyses, but it seems to be less important in engineering analyses, especially when completed thanks to the stochastic finite element methods. Such a conclusion may be drawn for two Gaussian distributions, and further analyses should be provided for other distributions and their pairs having more important and frequent applications in applied sciences and engineering.
Another important research finding reported in this work is that a lack of symmetry in Kullback-Leibler probabilistic divergence may lead to some imprecision while analyzing the reliability of even relatively simple engineering structures. This imprecision may be higher for large-scale structures, and it deserves further numerical studies. The homogenization analysis presented above shows that Kullback-Leibler divergence asymmetry may be very large while computing a distance between the original and homogenized elastic characteristics. This means that a definition proposed by Equation (6) returns completely different results while counting the distance between the original and homogenized tensors components or vice versa. Symmetric probabilistic divergence models are recommended in this context for further studies.
It should be noted that engineering practice also welcomes time-dependent probabilistic divergence analysis, where two distributions of various natures representing extreme effort and the corresponding limit change their parameters in the exploitation time. Then, this divergence may additionally depend on time in some unknown way, which deserves further studies, both practical and numerical. Let us note that some interesting applications would be a determination of the relative probabilistic entropy in conjunction with the probability transformation method (PTM) presented in [41].
Another interesting research avenue would be to contrast the statistics of the deformation computed via the SFEM models with those resulting from laboratory experiments or some in-situ observations. Then, a relative probabilistic entropy could serve as some calibration measure to improve numerical models and methods themselves. The same observation concerns effective material characteristics, which can be determined also in an experimental manner via uniaxial, and biaxial tensile tests. There is no doubt that modern advanced engineering of solid continua, such as optical fiber [42] behavior or specifically coupled and uncoupled fluid flow phenomena [43,44], may need a determination of probabilistic entropy and also probabilistic divergence. Such an application may be very attractive due to the final single graph of the relative entropy allowing for full identification of the resulting uncertainty. An application of the methodology presented in computer analysis of some acoustics and hydrodynamic phenomena [45] may find its importance in biomedical engineering and it deserves further research attention.

Funding

This paper has been written in the framework of the research grant OPUS no. 2021/41/B/ST8/02432 “Probabilistic entropy in engineering computations” sponsored by the National Science Center in Cracow, Poland.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Haldar, A.; Mahadevan, S. Probability, Reliability and Statistical Methods in Engineering Design; Wiley: New York, NY, USA, 2000. [Google Scholar]
  2. Hohenbichler, M.; Rackwitz, R. First-order concepts in systems reliability. Struct. Safety 1987, 1, 177–188. [Google Scholar] [CrossRef]
  3. Kullback, S. Information Theory and Statistics; Wiley: New York, NY, USA, 1959. [Google Scholar]
  4. Leinster, T. Entropy and Diversity. The Axiomatic Approach; Cambridge University Press: Cambridge, UK, 2021. [Google Scholar]
  5. Zienkiewicz, O.C.; Taylor, R.L.; Zhu, J.Z. The Finite Element Method: Its Basis & Fundamentals; Butterworth-Heinemann: Oxford, UK, 2013. [Google Scholar]
  6. Donald, M.J. On the relative entropy. Comm. Math. Phys. 1986, 105, 13–34. [Google Scholar] [CrossRef]
  7. Shannon, C.E. A mathematical theory of communication. Part I & II. Bell Sys. Tech. J. 1948, 27, 379–423, 623–656. [Google Scholar]
  8. Hellinger, E. Neue begrűndung der theorie quadratischer formen von unendlichvielen veränderlichen. J. Reine Angew. Math. 1909, 136, 210–271. [Google Scholar] [CrossRef]
  9. Kullback, S.; Leibler, R.A. On information and sufficiency. The Annals Math. Stat. 1951, 22, 79–86. [Google Scholar] [CrossRef]
  10. Jeffreys, H. An invariant form for the prior probability in estimation problems. Proc. Roy. Soc. London. Ser. A Math. Phys. Sci. 1946, 186, 453–461. [Google Scholar]
  11. Nielsen, F. On a variational definition for the Jensen-Shannon symmetrization of distances based on the information radius. Entropy 2021, 23, 464. [Google Scholar] [CrossRef]
  12. Nielsen, F. On the Jensen-Shannon symmetrization of distances relying on abstract means. Entropy 2020, 21, 485. [Google Scholar] [CrossRef]
  13. Harremoës, P. Interpretation of Rényi entropies and divergences. Phys. A: Stat. Mech. Its Appl. 2006, 365, 57–62. [Google Scholar] [CrossRef]
  14. Bhattacharyya, A. On a measure of divergence between two statistical populations defined by their probability distributions. Bull. Calcutta Math. Soc. 1943, 35, 99–109. [Google Scholar]
  15. Mahalanobis, P.C. On the generalized distance in statistics. Proc. Nat. Inst. Sci. India 1936, 2, 49–55. [Google Scholar]
  16. Kamiński, M. Tsallis entropy in dual homogenization of random composites using the Stochastic Finite Element Method. Int. J. Num. Meth. Engrg. 2018, 113, 834–857. [Google Scholar] [CrossRef]
  17. Fukui, K.; Sato, K.; Mizusaki, J.; Numao, M. Kullback-Leibler Divergence Based Kernel SOM for Visualization of Damage Process on Fuel Cells. In Proceedings of the 22nd IEEE International Conference on Tools with Artificial Intelligence, Arras, France, 27–29 October 2010; pp. 233–240. [Google Scholar] [CrossRef]
  18. Bilionis, I.; Koutsourelakis, P.S. Free energy computations by minimization of Kullback–Leibler divergence: An efficient adaptive biasing potential method for sparse representations. J. Comput. Phys. 2012, 231, 3849–3870. [Google Scholar] [CrossRef] [Green Version]
  19. Jiang, M.; Song, S.; Tang, F.; Li, Y.; Liu, J.; Feng, X. Scan registration for underwater mechanical scanning imaging sonar using symmetrical Kullback–Leibler divergence. J. Electron. Imag. 2019, 28, 013026. [Google Scholar] [CrossRef]
  20. Wei, S.; Wang, D.; Peng, Z. Principal component analysis based Kullback-Leibler divergence for die cracks detection. In Proceedings of the 2020 International Conference on Sensing, Measurement & Data Analytics in the era of Artificial Intelligence (ICSMD), Xi’an, China, 15–17 October 2020; pp. 224–228. [Google Scholar] [CrossRef]
  21. Filippatos, A.; Langkamp, A.; Kostka, P.; Gude, M. A Sequence-Based Damage Identification Method for Composite Rotors by Applying the Kullback–Leibler Divergence, a Two-Sample Kolmogorov–Smirnov Test and a Statistical Hidden Markov Model. Entropy 2019, 21, 690. [Google Scholar] [CrossRef]
  22. Kumar, P.S.; Kumaraswamidhas, L.A.; Laha, S.K. Bearing degradation assessment and remaining useful life estimation based on Kullback-Leibler divergence and Gaussian processes regression. Measurement 2021, 174, 108948. [Google Scholar] [CrossRef]
  23. Zhu, X.Q.; Hao, H.; Uy, B.; Xia, Y.; Mirza, O. Dynamic assessment of shear connection conditions in slab-girder bridges by Kullback-Leibler distance. Adv. Struct. Engrg. 2012, 15, 771–780. [Google Scholar] [CrossRef]
  24. Abbas, A.E.H.; Cadenbach, A.; Salimi, E. A Kullback–Leibler View of Maximum Entropy and Maximum Log-Probability Methods. Entropy 2017, 19, 232. [Google Scholar] [CrossRef]
  25. Contreras-Reyes, J.E. Asymptotic form of the Kullback–Leibler divergence for multivariate asymmetric heavy-tailed distributions. Phys. A Stat. Mech. Its Appl. 2014, 395, 200–208. [Google Scholar] [CrossRef]
  26. Garbaczewski, P. Shannon versus Kullback–Leibler entropies in nonequilibrium random motion. Phys. Lett. A 2005, 341, 33–38. [Google Scholar] [CrossRef]
  27. Ferrante, A.; Pavon, M.; Ramponi, F. Hellinger versus Kullback–Leibler multivariable spectrum approximation. IEEE Trans. Autom. Control 2008, 53, 954–967. [Google Scholar] [CrossRef]
  28. Martín, A.L.; López-Rosa, S.; Angulo, J.; Antolín, J. Jensen–Shannon and Kullback–Leibler divergences as quantifiers of relativistic effects in neutral atoms. Chem. Phys. Letters 2015, 635, 75–79. [Google Scholar] [CrossRef]
  29. Huang, J.; Yong, W.A.; Hong, L. Generalization of the Kullback–Leibler divergence in the Tsallis statistics. J. Math. Anal. Appl. 2016, 436, 501–512. [Google Scholar] [CrossRef]
  30. Nielsen, F. Fast approximations of the Jeffreys divergence between univariate Gaussian mixtures via mixture conversions to exponential-polynomial distributions. Entropy 2021, 23, 1417. [Google Scholar] [CrossRef]
  31. Kamiński, M. The Stochastic Perturbation Method for Computational Mechanics; Wiley: Chichester, UK, 2013. [Google Scholar]
  32. Kamiński, M. On iterative scheme in determination of the probabilistic moments of the structural response in the Stochastic perturbation-based Finite Element Method. Int. J. Num. Meth. Engrg. 2015, 104, 1038–1060. [Google Scholar] [CrossRef]
  33. Kamiński, M. Probabilistic entropy in homogenization of the periodic fiber-reinforced composites with random elastic parameters. Int. J. Num. Meth. Engrg. 2012, 90, 939–954. [Google Scholar] [CrossRef]
  34. D’Ambrosio, R.; Giordano, G.; Paternoster, B.; Ventola, A. Perturbative analysis of stochastic Hamiltonian problems under time discretizations. Appl. Math. Letters 2021, 120, 107223. [Google Scholar] [CrossRef]
  35. Jiang, D.; Zhou, B.; Han, B. Ergodic stationary distribution and extinction of a n-species Gilpin-Ayala competition system with nonlinear random perturbations. Appl. Math. Letters 2021, 120, 107273. [Google Scholar] [CrossRef]
  36. Eurocode 0, PN-EN-1990-1-1; Structural design basis. European Union: Brussels, Belgium, 2002.
  37. Kamiński, M.; Błoński, R. Analytical and numerical reliability analysis of certain Pratt steel truss. Appl. Sci. 2022, 12, 2901. [Google Scholar] [CrossRef]
  38. Oden, J.T. A general theory of finite elements. II. Applications. Int. J. Num. Meth. Engrg. 1969, 1, 247–259. [Google Scholar] [CrossRef]
  39. Eurocode 3, PN-EN-1993-1-1; Designing of steel structures—general rules. European Union: Brussels, Belgium, 2008.
  40. Kamiński, M. Design sensitivity analysis for the homogenized elasticity tensor of a polymer filled with rubber particles. Int. J. Sol. Struct. 2014, 51, 612–621. [Google Scholar] [CrossRef]
  41. Falsone, G.; Laudani, R. A probability transformation method (PTM) for the dynamic stochastic response of structures with non-Gaussian excitations. Engrg. Comput. 2018, 35, 1978–1997. [Google Scholar] [CrossRef]
  42. Yang, D.Y.; Tian, B.; Qu, Q.X.; Zhang, C.R.; Chen, S.S.; Wei, C.C. Lax pair, conservation law, Darboux transformation and localized waves of a variable-coefficient coupled Hirota system in an inhomogeneous optical fiber. Chaos Solitons Fractals 2021, 150, 110487. [Google Scholar] [CrossRef]
  43. Shen, Y.; Tian, B. Bilinear auto-Bäcklund transformations and soliton solutions of a (3+1) – dimensional generalized nonlinear evolution equation for the shallow water waves. Appl. Math. Lett. 2021, 122, 107301. [Google Scholar] [CrossRef]
  44. Zhou, T.Y.; Tian, B.; Chen, Y.Q.; Shen, Y. Painlevé analysis, auto-Bäcklund transformation and analytic solutions of a (2+1) – dimensional generalized Burgers system with the variable coefficients in a fluid. Nonlinear Dynamics 2022, 108, 2417–2428. [Google Scholar] [CrossRef]
  45. Gao, X.Y.; Guo, Y.J.; Shan, W.R. Taking into consideration an extended coupled (2+1) – dimensional Burgers system in oceanography, acoustics and hydrodynamics. Chaos Solitons Fractals 2022, 161, 112293. [Google Scholar] [CrossRef]
Figure 1. Kullback-Leibler (left) versus Jeffreys (right) entropy in the first test.
Figure 1. Kullback-Leibler (left) versus Jeffreys (right) entropy in the first test.
Symmetry 14 01945 g001
Figure 2. Kullback-Leibler (left) versus Jeffreys (right) entropy in the second test.
Figure 2. Kullback-Leibler (left) versus Jeffreys (right) entropy in the second test.
Symmetry 14 01945 g002
Figure 3. Static scheme of the given Pratt truss structure–Autodesk ROBOT 2021.
Figure 3. Static scheme of the given Pratt truss structure–Autodesk ROBOT 2021.
Symmetry 14 01945 g003
Figure 4. Spatial idealization of the RVE for periodic particulate composite.
Figure 4. Spatial idealization of the RVE for periodic particulate composite.
Symmetry 14 01945 g004
Figure 5. ABAQUS model of the particulate composites.
Figure 5. ABAQUS model of the particulate composites.
Symmetry 14 01945 g005
Figure 6. Relative probabilistic entropies of C 1111 and C 1111 ( e f f ) for the particulate composite; particle Young modulus uncertainty–left graph; matrix Young modulus uncertainty–right graph.
Figure 6. Relative probabilistic entropies of C 1111 and C 1111 ( e f f ) for the particulate composite; particle Young modulus uncertainty–left graph; matrix Young modulus uncertainty–right graph.
Symmetry 14 01945 g006
Figure 7. Relative probabilistic entropies of C 1122 and C 1122 ( e f f ) for the particulate composite; particle Young modulus uncertainty—(left) graph; matrix Young modulus uncertainty—(right) graph.
Figure 7. Relative probabilistic entropies of C 1122 and C 1122 ( e f f ) for the particulate composite; particle Young modulus uncertainty—(left) graph; matrix Young modulus uncertainty—(right) graph.
Symmetry 14 01945 g007
Figure 8. Relative probabilistic entropies of C 1212 and C 1212 ( e f f ) for the particulate composite; particle Young modulus uncertainty—(left) graph; matrix Young modulus uncertainty—(right) graph.
Figure 8. Relative probabilistic entropies of C 1212 and C 1212 ( e f f ) for the particulate composite; particle Young modulus uncertainty—(left) graph; matrix Young modulus uncertainty—(right) graph.
Symmetry 14 01945 g008
Table 1. Relative entropies versus output uncertainty in the Pratt truss structure.
Table 1. Relative entropies versus output uncertainty in the Pratt truss structure.
Numerical
Method
Data Typeα = 0.025α = 0.050α = 0.075α = 0.100α = 0.125α = 0.150
Probabilistic analytical
method
α(u)0.09140.18270.27360.36390.45330.5418
DKL6369.2781706.486846.114548.691415.147347.052
DJ6403.8831716.821851.966552.988418.740350.277
Stochastic perturbation techniqueα(u)0.09140.18270.27360.36390.45330.5417
DKL6369.2781706.486846.114548.691415.147347.052
DJ6403.8841716.821851.966552.988418.740350.277
Monte-Carlo simulation
approach
α(u)0.09140.18180.27220.36180.45070.5385
DKL6367.9281702.557842.733545.431411.813343.520
DJ6402.5421712.961848.615549.743415.412346.747
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kamiński, M. On the Symmetry Importance in a Relative Entropy Analysis for Some Engineering Problems. Symmetry 2022, 14, 1945. https://doi.org/10.3390/sym14091945

AMA Style

Kamiński M. On the Symmetry Importance in a Relative Entropy Analysis for Some Engineering Problems. Symmetry. 2022; 14(9):1945. https://doi.org/10.3390/sym14091945

Chicago/Turabian Style

Kamiński, Marcin. 2022. "On the Symmetry Importance in a Relative Entropy Analysis for Some Engineering Problems" Symmetry 14, no. 9: 1945. https://doi.org/10.3390/sym14091945

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop