Dynamic Reliability Design of Multicomponent Structure with Improved Weighted Regression Distributed Collaborative Surrogate Model Method

For dynamic reliability design of complex structures with multiple components, an improved weighted regression distributed collaborative surrogate model method (IWRDCSMM) is developed from the extremum response surface method (ERSM), decomposed-coordinated thought, and improved weighted regression principle. .e ERSM is used to address the dynamic reliability and sensitivity analyses of multicomponent structures and enhance the computing efficiency. .e decomposedcoordinated thought is applied to handle the relationship among multiple components. .e improved weighted regression method is used to find the efficient samples with smaller errors to improve the modeling accuracy. .e proposed method is first introduced for dynamic probabilistic analysis (including reliability analysis and sensitivity analysis) of multicomponent structures. .e method is then mathematically modeled by adopting the efficient samples selected based on the improved weighted regression method. Finally, the radial deformation dynamic probabilistic analysis of an aeroengine turbine blisk assembled by blade and disk is accomplished, in respect of the IWRDCSMM, fluid-thermal-structure interaction, and the randomness of input parameters within the time domain [0, T]. .e results illustrate that the reliability degree of turbine blisk radial deformation is 0.9951 when the allowable value uallow is 2.30×10 m, and all the input parameters affecting the turbine blisk radial deformation are gas temperature, angular speed, inlet velocity, outlet pressure, material density, and inlet pressure, successively. As revealed by the comparison of different methods, the IWRDCSMM has high fitting speed and simulation efficiency with the guarantee of accuracy. .e efforts of this study provide a promising dynamic probabilistic analysis technique for complex structures with multiple components and enrich mechanical reliability theory.


Introduction
Mechanical system generally comprises multiple complex structures, and the complex structure is assembled by multiple components. For instance, turbine blisk in an aeroengine contains blade and disk [1]. Each structure plays a key role and is subjected to rigorous workloads during operation. e failure of complex structures easily results in the disorder of the whole mechanical system and even causes a catastrophic accident in the working process [2]. erefore, it is required to investigate the probabilistic analysis of the complex structure to improve the performance and safety of the mechanical system with the consideration of the effects of dynamic loads. Probabilistic analysis generally includes reliability analysis and sensitivity analysis [3].
Lots of works have emerged on the probabilistic analysis. Chowdhury et al. [4] applied Monte Carlo (MC) simulation and finite element (FE) method to evaluate the reliability of the cracked structure. Naess et al. [5] employed MC simulation to estimate the reliability of the complicated structure considering the randomness of inputs. Tomasson and Soder [6] discussed the reliability evaluation of composite power systems using an improved way of applying MC simulation with the cross-entropy method. Jahani et al. [7] utilized interval MC simulation to perform the reliability analysis of the complex structure. Gao et al. [8] presented a mixed perturbation MC method for structural reliability analysis. Zhang et al. [9] studied multiobject reliability analysis of turbine blisk with multidiscipline under multiphysical field interaction by the MC method. Shayanfar et al. [10] sought a correction for the first-order reliability method by employing MC simulation to calculate the failure probability of a structure. Zhao and Ang [11] evaluated the reliability of the mechanical structure by the first-order third-moment method and random parameters. Alibrandi and Koh [12] used the first-order reliability method and secant hyperplane method to estimate the reliability of floating production systems. Lu et al. [13] used the secondorder fourth-moment method to derive the structural reliability based on the first-order fourth-moment method [14]. Lin et al. [15] proposed the score function method with kernel density estimation for structural stochastic sensitivity analysis. All the above efforts complete the probabilistic analyses of the complex structure by only adopting direct simulation methods such as the MC method and first-order reliability method. Although it is precise for these methods, the computational burden is rather unbearable for most of the computing platforms.
To address the computational efficiency with the acceptable accuracy, the surrogate model method (also known as the response surface method (RSM)) was developed for probabilistic analysis. Cheng and Li [16] employed RSM to solve the inverse structural reliability problems with implicit functions. Allaix and Carbone [17] used an improved RSM to evaluate structural reliability. Fang and Tee [18] proposed a novel method for structural reliability analysis based on RSM and improved genetic algorithm. Li et al. [19] applied the RSM to fulfill the seismic reliability evaluation of steeltimber hybrid shear wall systems. Cheng and Liu [20] performed the reliability analysis of steel cable-stayed bridges based on the RSM. Zhang and Bai [21] presented extremum RSM (ERSM) to derive the dynamic reliability analysis of a two-link flexible robot manipulator. Lu et al. [22] developed multiple ERSM to finish the dynamic probabilistic analysis of an aeroengine turbine blisk assembled by disk and blade. Fei et al. [23][24][25] derived the ERSM-based least-squares support vector machine and demonstrated it to be efficient in the nonlinear dynamic probabilistic analysis of turbine components. Kaymaz and Gaspar et al. applied the Kriging model method to study the structural reliability analysis [26,27]. Chojaczyk et al. [28] applied artificial neural network model to study the structural reliability analysis of steel structures. Lehky and Novak [29] used artificial neural network model to solve the inverse structural reliability problem. Cheng [30] integrated the advantages of the artificial neural network method, first-order reliability method, and importance sampling updating method for the serviceability reliability analysis of prestressed concrete bridges. e above efforts can accomplish the reliability and sensitivity analyses with acceptable efficiency, but the problem of unacceptable computational precision still exists since the limitation of quadratic polynomials is regarded as a surrogate model. Besides, these methods are only applicable to the probabilistic analysis of one component and one failure mode. erefore, there is an urgent need for the dynamic probabilistic analysis of the multicomponent structure to develop an efficient method.
is paper is oriented to develop an efficient and accurate analytical method for the dynamic reliability design of the multicomponent structure by integrating the ERSM [23], decomposed-coordinated thought [31,32], and the developed improved weighted regression method. is method is called improved weighted regression distributed collaborative surrogate model method (hereafter, IWRDCSMM). e ERSM is used to address the transient problem in the dynamic probabilistic analysis of the multicomponent structure to enhance the computing efficiency and accuracy. e decomposed-coordinated thought is applied to process the relationship among multiple components. e improved weighted regression method is used to find the samples with smaller errors or larger weight to improve the modeling accuracy. e model is adopted for the dynamic probabilistic analysis of the multicomponent structure within the time e dynamic probabilistic analysis of the radial deformation of the blisk is completed regarding fluidthermal-structural interaction, to demonstrate the feasibility and reasonability of the proposed method. e remaining of this paper is organized as follows. In Section 2, we elaborate the proposed IWRDCSMM. Section 3 applies this developed technique to the dynamic probabilistic analysis for the blisk radial running deformation, by considering the randomness of input variables of inlet velocity, inlet pressure, outlet pressure, gas temperature, material density, and angular speed. Section 4 verifies the feasibility and effectiveness of the method from the accuracy and efficiency of modeling and simulation by the comparison of the IWRDCSMM with other methods. Some conclusions on this study are given in Section 5.

Basic Theory
is section mainly introduces the basic theory of the IWRDCSMM. We first describe the procedure of the dynamic probabilistic analysis for multicomponent structure with the method. e details of ERSM and IWRDCSMM are then given in next two subsections. We finally discuss the thoughts of structural dynamic reliability and sensitivity analyses.
As illustrated in Figure 1, the procedure of multicomponent structural dynamic probabilistic analysis with IWRDCSMM mainly includes dynamic deterministic analysis, mathematical modeling, and dynamic probabilistic analysis. e details of this technique are summarized as follows: (1) Define objective structure with multicomponent, ensure analytical time domain [0, T], and set dynamic loads and boundary conditions. (2) Select influential factors as random input variables, ascertain their numerical features, and then derive the dynamic deterministic analysis of multicomponent structure. (3) Find the extreme values of analytical object and determine the time point and obtain samples of input parameters and output responses with linkage sampling technique [9]. e linkage sampling method fuses the extremum selection method [33] and Latin hypercube sampling method [34,35], which is defined as using a group input sample to obtain the output responses of multiple objectives at the same time. e linkage sampling method with one input parameter x and two output responses y 1 and y 2 is shown in Figure 2.     [21,23,36]. e mathematical model of this method is built in (1) with quadratic polynomial function, in which y QP (x) indicates the extreme value of output response and x expresses the vector of input variables.
in which a is the unknown coefficient of constant, b is the vector of unknown coefficients of linear term, and c is the matrix of unknown coefficients of the quadratic term. Besides, b, c, and x can be written as follows: where k is the number of input parameters. Equation (1) is rewritten as quadratic polynomial function on basis of (2), that is, in which a, b i , and c i are the unknown coefficients of the constant term, linear term, and quadratic term, respectively. To obtain the unknown coefficients a, b i , and c i in (3), the least square method is used to fit the extremum response surface model. Its theory is denoted by where w presents the vector of unknown coefficients, d expresses the matrix of input samples, and y indicates the vector of output response. w, d, and y are shown as follows: where n is the number of samples of input parameters, and x ik is the kth input parameter of the ith sampling date.

Improved Weighted Regression Distributed Collaborative
Surrogate Model Method. ERSM can efficiently address the dynamic probabilistic analysis of single structure. However, it is difficult for this technique to resolve the dynamic reliability and sensitivity analyses of the multicomponent structure. Accordingly, the decomposed-coordinated thought is introduced into the ERSM, which is called as distributed collaborative surrogate model method (DCSMM) [37]. e multicomponent structure always refers to many subcomponents. We consider the structure with m subcomponents to investigate the DCSMM, and then the mathematical model of the DCSMM is established with respect to the least square method. e coordinated function can be written by in which y (j) QP (x) is the decomposed function of the jth subcomponent and is expressed as where a (j) is the unknown coefficient of the constant term of the jth subcomponent, b Despite the DCSMM reasonably handles the dynamic probabilistic analysis of the complex structure [37], the least square method is difficult to handle high-nonlinear function and large-scale samples in mathematical modeling. e conventional weighted regression principle is used to establish the mathematical model with few samples relative to the least square method. e weighted values are obtained by the ratio of the minimum output response to the output responses [38,39]. e selected minimum output response may not be the optimal value, so that the fitting accuracy of mathematical model is unprecise. In this case, the improved weighted regression (IWR) method is developed to compute the weighted values by the ratio of the minimum error to the errors, and the efficient samples with smaller errors or larger weights are selected as the training samples of the fitting model. e errors are the absolute difference between the real value and the computing value. e minimum error is the minimum value of all obtained errors. e efficient samples are the samples with smaller errors or larger weighted values, selected from all samples. e IWR model is established using the selected efficient samples, the fitting performance of which is superiority to the least square method. e detailed validation procedure is stated in Section 4.
In the light of (6) and (7), the formula of the mathematical model of IWRDCSMM is 4 Advances in Materials Science and Engineering where x is the vector of efficient input samples, x is the value of the ith efficient input sample of the jth subcomponent, Y IWR (x) represents the coordinated function that established the IWR, y i are, respectively, the unknown coefficients of the constant term, linear term, and quadratic term, which can be obtained via the IWR, that is, in which w (j) is the vector of unknown coefficients of the jth subcomponent for the IWRDCSMM, d (j) is the matrix of input samples of the jth subcomponent, v (j) is the weighted matrix of efficient samples, and y , and y (j) (10), and the weighted matrix v (j) is defined by (11).
where n is the number of efficient samples.
in which i � 1, 2, . . . , n; e obj is the objective value of the jth subcomponent, namely, the minimum error, (y (j) i (x)) QP is the ith efficient output response of jth subcomponent by using the least square method, and (y (j) i (x)) QP is the ith efficient output response of jth subcomponent via the FE method.
Based on the above theories, the unknown coefficients in (8) are acquired, and the mathematical model of IWRDCSMM is established.

Dynamic Probabilistic Analysis.
e IWRDCSMM is utilized to implement the dynamic probabilistic analysis (reliability analysis and sensitivity analysis) of complex structure with multiple components. In respect of the mathematical model of the IWRDCSMM, the limit state function of analytical object as Y(x) is where Y allow is the allowable value of analytical object, x is the vector of extracting input samples, and x (j) i is the value of the ith input variable of the jth subcomponent. Besides, the safety of complex structure is proved by Y(x) ≥ 0, and Y(x) < 0 indicates that the multicomponent structure failure.
To evaluate the dynamic reliability of multicomponent structure, the MC method is employed to achieve enough samples to obtain the values of limit state function with regard to (12) [40]. By the Chebyshev theorem and Bemoulli theorem [41] in probability theory, the structural failure and reliability probabilities are in which N and N f refer to the total number of extracting samples and the number of failures (as defined in (12) with a given Y allow ), respectively. Sensitivity analysis of multicomponent structure can reflect the effect of the variation of input variables on failure probability or reliability probability, which help us determine the important degrees of input parameters on the output response of analytical object and then provide guidance for the design parameters of complex structures [1,23,40]. Sensitivity analysis typically contains sensitivity degree (S d ) analysis and impact probability (I d ) analysis [25]. Sensitivity degree reflects the influence level of an input variable on the output response, and impact probability indicates how much an input variable affects the output response relative to the other input parameters.
e sensitivity degree S d is defined by the partial derivative of failure or reliability probability function with respect to the mean of input parameters, that is,

Advances in Materials Science and Engineering
where χ is the reliability index, the distribution features of input variables obey normal distribution and are mutually independent. erefore, zP f /zχ can be expressed by in which μ Y and σ Y are the mean value and standard deviation of the limit state function, respectively. Equation (14) is rewritten as follows: e impact probability I d is defined by the ratio of the absolute sensitivity degree of one input parameter to the sum of the absolute sensitivity degrees of all input variables, that is, In addition, the relationship of the impact probability of all input variables is In the light of the above analyses, the reliability degree and sensitivity indexes (sensitivity degree and impact probability) of random input variables on failure probability are acquired to evaluate the design quality of the complex structure.

Case Study: Dynamic Probabilistic Analysis of Turbine Blisk
In this section, the high-pressure turbine blisk of an aeroengine, which comprises blade and disk, is selected to demonstrate the proposed method in this paper. e dynamic probabilistic analysis of turbine blisk radial running deformation is accomplished to verify the feasibility and applicability in engineering practice, by considering the fluid-thermal-structural interaction.

Finite Element Model and Input Variables.
As an important structure of an aeroengine, high-pressure turbine blisk suffers high gas temperature and pressure from main combustion chamber and high centrifugal force due to angular speed. e schematic diagram of a blisk is shown in Figure 3. From Figure 3, it is obvious that the turbine blisk is a typical cyclic and symmetric structure. To save modeling burden and improve simulation time, 1/48 of the blisk model is selected as the object of study. e simplified threedimensional models of the blisk and flow field are drawn in Figures 4 and 5, respectively. In Figure 5, the FTSI surface indicates the fluid-thermal-structural interaction surface which is to transfer fluid loads and thermal loads to the structure.
In view of Figures 4 and 5, the FE/FV models of turbine blisk and flow field are established in Figures 6 and 7, respectively. As illustrated in Figures 6 and 7, the FE models of turbine blade and disk are consisted of tetrahedrons with 48,551 elements, 75,097 nodes and 33,885 elements, 57,213 nodes, and the FV model of flow field has the tetrahedrons with 338,917 elements and 472,930 nodes.
To study the influence of fluid-thermal-structural interaction on the radial deformation of turbine blisk, nickelbased superalloy GH4133 [42] is selected as the material of turbine blisk. e variation of angular speed within a time domain [0 s, 215 s] is shown in Figure 8. Here, 12 points are sampled as the calculation points [43]. e inlet velocity, inlet pressure, outlet pressure, and gas temperature are regarded as input variables to perform the deterministic analysis, and their means are assumed to be 160 m/s, 2 × 10 6 Pa, 5.88 × 10 5 Pa, and 1,200 K, respectively.
Substantially, the input variables process uncertainties and affect the blisk radial deformation, so that the randomness of input parameters needs to be considered. erefore, the above parameters are regarded as random variables. Assuming that all random input variables obey normal distribution and are independent mutually; the distribution characteristics of the random variables are listed in Table 1.

Dynamic Deterministic Analysis.
Based on the determined values of input variables, the dynamic deterministic analysis of blisk radial deformation is carried out by the closecoupling analysis technique [22] with fluid-thermal-structural interaction. e procedure of dynamic deterministic analysis of blisk radial deformation contains fluid analysis, thermal analysis, structural analysis, and iterative solution. e fluid analysis is based on the standard k-ε turbulence model with respect to the boundary conditions of the flow field and finite volume method. With regard to the thermal analysis, the law of energy conservation is utilized to resolve Inlet Outlet FTSI surface       Advances in Materials Science and Engineering the property of turbine blisk flow flied, and then the pressure and gas temperature are transferred to turbine blisk through the FTSI surface. e structural analysis is performed by means of the FE method, in which the tetrahedron shape equation and geometric equation are applied to perform blisk radial deformation analysis. e dynamic deterministic analysis of radial deformation is finished by using the iterative solution.
According to the theory of the close-coupling analysis method, the dynamic deterministic analysis is accomplished, and the changing curves of the radial deformations of blade and disk with time are shown in Figure 9.
As shown in Figure 9, it is obvious that the radial deformation of turbine blisk increases with the increase of angular speed and the maximum emerges in the time domain [165 s, 200 s]. T � 175 s is chosen as the computing point. At this moment, the distributions of pressure on FTSI surface, temperature on turbine blisk, and the nephograms of blade and disk radial deformations are displayed in Figures 10-12, in which P indicates pressure, t expresses temperature, and u b and u d represent the radial deformations of blade and disk, respectively.
As illustrated in Figure 12, the maximum of radial deformation of turbine blisk occurs at the tip of blade and disk. Hence, the risk section responses of turbine blade and disk tip are selected as the output responses of the blisk and are also the maximum output responses of the blisk. e blisk is safe as the risk section is secure. erefore, the output responses of turbine blade and disk tip are treated as the computational points in the blisk dynamic probabilistic analyses.

Derivation of Mathematical Model.
To derive the IWRDCSMM model for the blisk radial deformation, the linkage sampling method is applied to extract 60 samples for random input variables following the numerical characteristics in Table 1. And then 60 samples for the corresponding output response are obtained by dynamic deterministic analyses. From this pool of 60 samples, 40 samples (as shown in Table 2) are selected to establish the initial mathematical model of blade and disk radial running deformations.
e two models are displayed in (19). e remaining 20 samples are treated as the testing samples for the purpose of accuracy assessment.
Based on (19) and (20), samples with smaller errors underlined and bold in Table 2 are chosen from 40 samples and are used to derive the IWRDCSMM model of blade and disk in (20). Advances in Materials Science and Engineering

Dynamic Probabilistic Analysis.
In view of the IWRDCSMM model in (20), the limit state function of the blisk radial deformation is where u allow is the allowable value of radial deformation and v, p in , p out , t, ρ, and w are the extracting samples of input variables. For the dynamic reliability analysis of blisk radial deformation, the allowable value u allow is 2.30 × 10 −3 m in respect of engineering experience. Based on the value and the limit state function, the MC method is applied to obtain 10,000 samples of input variables. e simulation history and histogram distribution of the blisk radial deformation are shown in Figures 13 and 14, respectively, in which u max indicates the maximum of blisk radial deformation.
From Figures 13 and 14, it is obvious that the radial deformation of turbine blisk follows a normal distribution with the mean μ � 2.19 × 10 −3 m and the variance σ 2 � 1.87 × 10 −9 m. And the reliability degree R � 0.9951 is acquired, which satisfies the requirement of engineering design.
Additionally, the sensitivity analysis of blisk radial running deformation with respect to all input parameters is implemented.
e analytical results including sensitivity degree and impact probability are obtained as depicted in Figure 15 and Table 3, respectively.
As demonstrated in Figure 15 and Table 3, the sensitivity degrees of all input parameters have either positive or negative values. e signs are related to whether the variation of an input variable has a positive or negative correlation with the output response [23,33]. Hence, the blisk radial deformation increases with the increase of inlet pressure, outlet pressure, gas temperature, material density, and angular speed, while with the decrease of the inlet velocity. Besides, the impact probability reveals the importance of an input variable relative to other variables. e importance of all input parameters affecting the blisk radial running deformation are successively gas temperature, angular speed, inlet velocity, outlet pressure, material density, and inlet pressure.
In view of the above analysis, the effects of all input parameters on the blisk radial running deformation are acquired.
is information might be used to guide the design of turbine blisk as well as update the model of the complex structure.

IWRDCSMM Verification
is section is to validate the effectiveness and feasibility of the developed IWRDCSMM. e fitting properties of both fitting efficiency and computing accuracy are firstly focused to validate the IWRDCSMM by comparing with distributed collaborative surrogate model method-based least square method (LSDCSMM) and distributed collaborative surrogate model method-based weighted regression principle (WRDCSMM). And then, these models and MC simulation are utilized to implement the dynamic reliability analysis of turbine blisk radial deformation to verify the superiority of the proposed method.

Fitting Property.
is subsection is to verify the IWRDCSMM from the fitting property (including fitting efficiency and accuracy) perspective. e three models (LSDCSMM model, WRDCSMM model, and IWRDCSMM model) are built by the samples in Table 2. erein, the forms of LSDCSMM and IWRDCSMM models are shown in (19) and (20). e formula of the WRDCSMM model is modeled in (22) with 20 samples by using conventional weighted regression principle [44,45].

Advances in Materials Science and Engineering
To support the IWRDCSMM model, the remaining 20 samples are selected as the testing samples. ese models ((19), (20), and (22)) are verified by the below errors.
in which i � 1, 2, . . . , 20; u i,test is the testing value of radial deformation for LSDCSMM, WRDCSMM, or IWRDCSMM corresponding to the ith sample; and u i,true is the true value of radial deformation based on the FE model for the ith sample. e analytical results including the fitting time and accuracy are listed in Table 4.
As revealed in Table 4, the fitting time of the IWRDCSMM model is less than that of the LSDCSMM model and is equal to that of the WRDCSMM model. e reason is that the number of samples for the IWRDCSMM is less than that of LSDCSMM and is equal to that of the WRDCSMM in mathematical modeling. e maximum absolute error, minimum absolute error, and average absolute error for the IWRDCSMM are only 0.155, 0.001, 0.026 and 0.904, 0.500, 0.807 for these of the LSDCSMM and WRDCSMM, respectively. e IWRDCSMM model holds less error and higher precision than the LSDCSMM and WRDCSMM.

Computational Property of Dynamic Reliability Analysis.
To further validate the effectiveness and feasibility of the developed method, MC simulation, LSDCSMM, WRDCSMM, and IWRDCSMM are employed to perform the dynamic reliability analyses of turbine blisk radial running deformation based on the same input variables in Table 1 and computational environment. e computing time and precision of different simulations for different methods are listed in Tables 5 and 6.
As shown in Table 5, in the computing time the IWRDCSMM is far less than the MC simulation and is equal to the LSDCSMM and WRDCSMM. In the computational efficiency, the IWRDCSMM becomes higher with the increasing simulations. As demonstrated in Table 6, the analytical precision of the IWRDCSMM is almost consistent with that of the MC simulation and is higher than those of LSDCSMM and WRDCSMM. As revealed in Tables 5 and 6, the proposed method can effectively accomplish the dynamic reliability analysis of complex structure with multiple components, while the MC simulation cannot when the number of simulation is over 10 3 .
In short, it is fully illustrated that the IWRDCSMM holds high precision and efficiency in modeling and simulation for the nonlinear transient reliability analysis of turbine blisk. e promising approach in this work is testified by the dynamic reliability and sensitivity analyses of complex structure with multiple components.

Conclusions
e main aim of this paper is to develop an improved weighted regression distributed collaborative surrogate model method (IWRDCSMM) with respect to the extremum response surface method (ERSM), decomposed-coordinated thought, and improved weighted regression (IWR), for the dynamic probabilistic analysis of complex structure with multiple components. e dynamic reliability and sensitivity analyses of the radial running deformation for an aeroengine turbine blisk comprising blade and disk are implemented by considering the influence of fluidthermal-structural interaction to verify the effectiveness and feasibility of the proposed method. Some conclusions are acquired as follows: (1) rough the dynamic deterministic analysis of aeroengine turbine blisk, the changing rules of the blade and disk radial running deformations from start to cruise for aircraft are gained by adopting the close-coupling analysis technique within the time domain [0, T]. e turbine blisk radial running deformation reaches the maximum at t � 175 s.  (2) e reliability degree of radial running deformation of turbine blisk is R � 0.9951 based on the dynamic reliability analysis of IWRDCSMM, when the allowable value of turbine blisk radial running deformation u allow is 2.30 × 10 −3 m. (3) From sensitivity analysis, the importance of all input parameters affecting the blisk radial deformation are successively gas temperature, angular speed, inlet velocity, outlet pressure, material density, and inlet pressure. (4) In fitting property, the IWRDCSMM needs smaller samples in modeling than other surrogate models and holds higher fitting efficiency and accuracy. (5) By the comparison of different methods, the proposed IWRDCSMM is promising to enhance computing efficiency with the guarantee of acceptable analytical precision in the dynamic probabilistic analysis of multicomponent structure.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this article.