1. Introduction
Cartridge delivery systems are often used in aerospace engineering to transport cylindrical projectiles. The coordination mechanism is a crucial component of a cartridge conveying system, primarily responsible for receiving paper cartridges from external sources, storing and aligning them to the desired angle, and subsequently propelling them into the pipeline. Its primary function involves the reciprocating movement of the load between the drug receiving position and drug delivery position. In addition to ensuring its normal operation, it must also possess the ability to swiftly and accurately coordinate in place while maintaining a lightweight design under high loads. Based on actual test records, various uncertain factors, such as flexible vibrations within the mechanism and clearances between components, can significantly impact the stability and reliability of the coordination mechanism. Therefore, conducting an uncertainty analysis of this coordination mechanism holds immense significance in advancing research on cartridge conveying systems.
There are various quantitative analysis methods for uncertainty in mechanism dynamics systems. The central idea behind these methods is to clarify different sources of uncertainty in order to fully capture the stochastic behavior of dynamics systems so as to quantify the uncertainty consequences through numerical or experimental methods [
1]. At present, mainstream uncertainty propagation methods mainly include Monte Carlo simulation (MCS), the local expansion method, the orthogonal function expansion method, the numerical integration method, etc. The Monte Carlo method is also known as the random sampling method or statistical test method. With an increase in simulation times, its calculation results will gradually approach an accurate solution. Thanks to the development of computer technology, the Monte Carlo simulation method has been widely used in many engineering fields, and it is usually used as a quasi-accurate calculation method to verify the accuracy of other approximate methods [
2]. However, as this method requires a large number of simulation tests to ensure the accuracy of the calculation results, it has low calculation efficiency and high calculation costs, which make it difficult to apply to many practical engineering problems, especially complex large-scale problems. By improving sampling technology, several variants of the direct Monte Carlo method, such as the important sampling method [
3,
4], the subset simulation method [
5,
6,
7], and the direction sampling method [
8], have been proposed and developed, which can greatly improve the solving efficiency of the Monte Carlo method under certain conditions.
The above methods were all developed based on probability theory. There are two main problems in applying probabilistic methods to uncertainty analysis. First, obtaining an accurate trend of the probability distribution of uncertain variables depends on an extensive statistical analysis of a large number of test samples. This method involves huge computational complexity and may not be feasible or cost-effective for many practical engineering problems. Second, because it is difficult to obtain accurate parameter distributions, subjective assumptions are often made when using probabilistic methods for uncertainty analysis. However, even if there is only a small deviation in the parameter distribution, the results of an uncertainty analysis will produce significant errors [
9]. Non-probabilistic methods can effectively solve these shortcomings. Non-probabilistic models do not require a detailed understanding of the sample characteristics of uncertain variables but focus on determining the size or boundary of variable uncertainty, which is relatively easy. Second, even if the data are limited, the size or upper and lower boundaries of variable uncertainty can be more accurately determined, thus improving the accuracy of uncertainty analysis.
Since the middle and late 1980s, in order to overcome the dependence of classical probability methods on large sample sizes, the interval model has been introduced into the field of structural engineering to describe the fluctuation range of uncertain parameters. In the classical interval model [
10], a single uncertain parameter is described as a fluctuation-bounded interval variable, while the uncertain domain of multiple independent parameters constitutes a “multi-dimensional box”. At present, in the uncertain analysis problem, commonly used interval analysis methods can be divided into the vertex method (VM) [
11,
12], the configuration method [
13,
14,
15,
16], the Taylor expansion method [
17,
18,
19], etc. The vertex method assumes that the structural response obtains its response boundary at the vertex of the uncertain domain. However, the calculation speed of the vertex method increases exponentially with an increase in the interval parameters, which is unacceptable for large structures. The configuration method constructs a reasonable surrogate model by setting a certain number of sample points in the uncertain parameter space so as to obtain a range of structural responses. In order to meet certain accuracy requirements, a large number of sample points are usually required. Especially when there are many interval parameters, the calculation cost increases exponentially and then gets into a dimension disaster. The Taylor expansion method is used to solve the uncertain static response problem. This method involves carrying out the first- or second-order Taylor series expansion of the structural response function at the midpoint of the interval to obtain the value range or interval boundary of the structural response. The biggest advantage of this method is that it has high solving efficiency, and the calculation cost increases linearly with an increase in the interval parameters. The Taylor expansion method only considers the lower-order terms of the Taylor series expansion, ignoring the higher-order terms. When the interval parameters are largely uncertain, it is difficult to use the first-order Taylor expansion method to obtain high-precision response results, and it even leads to a distortion of the results due to linear approximation [
20].
For the interval analysis problem with a high degree of uncertainty or nonlinearity, Qiu and Elishakoff [
21] proposed a subinterval prediction method, which divided the original interval difference into several continuous subintervals, performed an interval analysis on these subintervals, respectively, and finally combined the interval analysis results. Zhou et al. [
22] used the subinterval prediction method to solve the structural analysis problem with large uncertain parameters. Chen et al. [
23] proposed a subinterval homogenization method to identify the elasticity characteristics of periodic microstructure problems with interval uncertainty and compared the results obtained with the MCS method to verify the accuracy of the method. Fu et al. [
24] proposed a dimension reduction subinterval prediction method to combine the dimension reduction method and the subinterval prediction method to deal with strong nonlinear and high-dimensional problems. The subinterval prediction method has good calculation accuracy, but there are still problems in determining the number of subintervals and calculation efficiency, especially for high-dimensional problems, in which the subinterval prediction method will produce a “dimensionality explosion”, resulting in low calculation efficiency.
The proxy model is a response function approximation constructed through an experimental design which serves as an approximate technique for addressing multivariable problems. It significantly enhances computational efficiency while ensuring the fitting accuracy of the proxy model and can be easily integrated with other methods, making it widely employed in the field of engineering design at present. Currently, commonly used proxy models primarily encompass machine learning [
25,
26] and deep learning neural network models. Among them, the deep learning neural network model consists of multiple hidden layers and exhibits exceptional fitting capability for complex, nonlinear relationships. Henceforth, this paper employs the DNN proxy model as a substitute for the high-precision dynamic model for simulation calculations in order to reduce computational costs.
The present paper establishes a rigid–flexible coupling dynamics model of the articulated gap to investigate the coordination mechanism of an ammunition automatic loading system. Experimental tests are conducted to validate the accuracy and effectiveness of this dynamics model. Subsequently, in order to reduce computational complexity during the uncertainty analysis of the coordination mechanism, a deep neural network surrogate model is employed for simulation calculations as a substitute for the high-precision dynamics model. Furthermore, building upon the classical subinterval prediction method based on Taylor expansion, we propose an improved directional subinterval prediction approach that incorporates gradient information within specific intervals. An adaptive convergence method is also designed to enhance calculation efficiency. Finally, we apply this methodology to analyze positioning accuracy uncertainties in the coordination mechanism and compare it with Monte Carlo simulations (MCSs) and other analysis methods, demonstrating its superior efficiency and accuracy.
4. Uncertainty Propagation Analysis Method Based on Subinterval Prediction
In this section, a subinterval prediction-based uncertainty propagation analysis method is proposed and applied to the uncertainty analysis of the coordinating mechanism. This method can accurately predict the upper and lower boundaries of the uncertainty output, providing a theoretical basis for subsequent parameter optimization.
4.1. Directional Subinterval Prediction Method
Taking the two-dimensional uncertain input vector,
, interval analysis problem as an example, for convenience, assume that both input parameter intervals are divided into
subintervals and stipulate that
is
, where
is a non-negative integer. One of the difficulties of the classical subinterval prediction method is that it needs to perform a first-order Taylor approximation on all the subinterval combinations, and in the two-dimensional problem, it needs to perform a first-order Taylor approximation on the midpoint of all the subinterval combinations. To solve this problem, this paper proposes a directional subinterval prediction method, which can reduce the number of points to be calculated. Firstly, select any subinterval,
, in the subinterval domain of the input vector
, and according to the classical subinterval prediction method, conduct a first-order Taylor expansion at the midpoint,
, of the subinterval with an amplitude of
:
where
. The upper and lower bounds of
can be expressed as follows:
To avoid performing first-order Taylor expansions at each central point of the subinterval domain, a reasonable approach is to find an expansion path and predict the upper bound of
by applying first-order Taylor expansions at multiple points on the path. When predicting the upper bound of
, the expansion path of extension points should move in the direction of the increasing output response function. Therefore, according to the gradient of
relative to
at
, the next expansion point,
, required to predict the upper bound can be expressed as follows:
The first-order Taylor expansion of
with an amplitude of
at
is as follows:
where the subinterval
. The upper and lower bounds of
obtained from (26) and (27) can be expressed as follows:
When the above steps are performed
times in the same way, all the expansion points,
, used to predict the upper bound can be uniformly expressed by the following equation:
At the above point, the output response is a first-order Taylor expansion with an amplitude of
, and the upper and lower bounds of all the expansion points’ output responses are obtained as follows:
where the subinterval
, and the upper bound of
can be obtained by interval fusion:
Similar to the method used for calculating the upper bound of the output response value interval, when predicting the lower bound of
, the expansion path of the extension point should move along the direction of the decreasing output response function. Therefore, according to the gradient of
relative to
at
, all the expansion points,
, required for predicting the lower bound can be uniformly expressed as follows:
According to (33) and (34), the upper and lower bounds of the output responses of all the expansion points are as follows:
where the subinterval
. According to Equation (35), the lower bound of
can be obtained by interval fusion:
The directional subinterval prediction is based on gradient analysis by constructing two expansion paths from to and and selecting two subinterval sets along different expansion paths. Therefore, when predicting the response interval, only a small number of points on the paths with increasing or decreasing gradients are selected for the first-order Taylor expansion. Compared with the classical subinterval analysis method, this method can greatly reduce the amount of calculation needed. In order to simplify the analysis process, in the above description of the two-dimensional problem, each input parameter is divided into the same number of subintervals, while different input parameters can be set to different numbers of subintervals according to the needs of the actual problem.
4.2. Adaptive Strategy for Subinterval Partitioning
In the subinterval prediction analysis method, the accuracy of the prediction results is closely related to the number of subintervals. When the input parameter interval is divided into a sufficient number of subintervals, both classical subinterval prediction and improved directional subinterval prediction can achieve sufficient accuracy. However, an increase in the number of subintervals means an increase in the number of calculations and a decrease in calculation efficiency. Therefore, on the premise of ensuring accuracy, it is an important measure to reduce the number of subintervals divided by the input parameter interval as much as possible to improve the calculation efficiency.
According to the directional subinterval prediction method in the previous section, the output response is first expanded by the first-order Taylor series at the midpoint of the input parameter interval, and the upper and lower bounds of the initial interval,
and
, can be obtained by Equations (30) and (31). If the initial interval cannot meet the prediction accuracy of the output response, a subinterval division of the input parameter interval is needed. Therefore, after obtaining the initial interval,
, the parameter interval is divided into two equal subintervals by the midpoint of the interval, where the first-order Taylor series is first expanded. Then, the first-order Taylor series is, respectively, expanded at the midpoints of the two subintervals, and the new upper and lower bounds of the interval
can be obtained by Equations (35) and (39). The adaptive strategy for subinterval division is shown in
Figure 11 below.
Global convergence is guaranteed on the premise of ensuring accuracy. The convergence index is defined as follows:
where
is the number of adaptive iterations for subinterval division;
and
are the upper and lower bounds of the output response value interval calculated from
iterations;
and
are the upper and lower bounds calculated through (
) iterations; and
and
are the defined convergence indexes, namely, the accuracy indexes representing the output response prediction.
4.3. Computing Process
To sum up, the main flow chart of the directional subinterval prediction method is shown in
Figure 12, and the specific implementation steps are as follows:
Step 1: Summarize the main influencing factors according to the mechanism failure mode, list the uncertain input parameters, and determine the uncertain input parameters by referring to the previous structural design, specified indicators, and existing experimental data.
Step 2: Determine the uncertain output response according to the mechanism’s uncertainty problem and construct the output response function, . Let .
Step 3: Find the appropriate expansion path for predicting the upper and lower bounds of the output response according to the gradient analysis of .
Step 4: Calculate the upper and lower bounds of the value interval of and , respectively, according to Equations (33) and (34) and Equations (37) and (38). Predict the upper and lower bounds of the interval according to Equations (35) and (39) through interval fusion.
Step 5: Test whether the upper and lower bounds of obtained through the iterations meet the specified convergence index. If they do meet it, the output response prediction interval is obtained; if not, proceed to Step 6.
Step 6: Divide the value interval of the input parameters into two equal subintervals with the midpoint of the interval as the dividing line, take the midpoints of the two new subintervals as the expansion point, and return to Step 3, with .
5. Analysis of Uncertainty Propagation in Coordination Mechanism Based on Subinterval Prediction
Based on the DNN surrogate model of the coordinating mechanism established in the previous section, the uncertainty of the positioning accuracy of the coordinating mechanism is quantitatively analyzed by using the directional subinterval prediction method and the subinterval division and convergence strategy. According to the theoretical design and engineering practice of the coordinating mechanism, the value interval and uncertainty level of the uncertainty input parameters of the coordinating mechanism are shown in
Table 2.
The uncertain input parameter interval is brought into the DNN surrogate model of the coordinating mechanism, and the directional subinterval prediction method is used to calculate the interval boundary of the output response of the coordinating mechanism, where the convergence index is set to
. In addition, in order to determine the accuracy of the proposed method, the results of the MCS method are used as the accurate solution, where the number of iterations in the MCS method is
, and the results of the first-order Taylor expansion method and the classical subinterval prediction method are used as the reference solution. The results are shown in
Figure 13 and
Table 3.
The output response boundary results of the coordinating mechanisms using the directional subinterval prediction are shown in
Figure 13a. As can be seen from the figure, the convergence results of the output response boundary are obtained only through two adaptive iterations, in which each uncertain input parameter interval is divided into two subintervals. The upper and lower boundaries of the output response of the coordinating mechanism are very close to the exact solution calculated by the MCS method, with relative errors of only 0.72‱ and 1.10‱, and the accuracy can be improved again by increasing the number of subintervals. The results of the first-order Taylor expansion method and the classical subinterval prediction method are shown in
Table 3. From the table, the relative errors of the lower boundary of the output response obtained by the first-order Taylor expansion method and the classical subinterval prediction method are 10.19‱ and 6.85‱, respectively. The results of the directional subinterval prediction method and the classical subinterval prediction method are in good agreement with the exact solution. However, the directional subinterval prediction method only needs 27 iterations of the surrogate model, while the classical subinterval prediction method needs a subinterval analysis of all the possible subinterval combinations, so it needs 2304 iterations. Therefore, the directional subinterval prediction method proposed in this paper can greatly reduce the calculation cost and improve the calculation efficiency under the premise of ensuring accuracy, indicating that this method can be used in uncertainty analysis and follow-up research on coordinating mechanisms.
6. Conclusions
This paper investigates the coordination mechanism of a cartridge conveying system by constructing a rigid–flexible coupling dynamics model of the coordination mechanism with a hinged gap, based on models of shaft–hole hinged gaps and flexible cantilever beams. Through conducting experiments on the principle prototype platform and analyzing recorded data, we compare the output response under actual conditions with that obtained from simulation using the dynamics model to verify its accuracy and effectiveness. However, due to the numerous calculation models required in uncertainty analysis, employing a high-precision dynamics model would result in significant computational costs. Therefore, we utilize a deep neural network surrogate model as an alternative for the simulation calculations instead. Subsequently, we propose a directional subinterval prediction method based on classical subinterval prediction methods and employ an adaptive subinterval partition method to accurately approximate output response interval boundaries while reducing the use of the model. This approach is then applied to predict the positioning accuracy of the coordination mechanism under uncertain input parameters, obtaining an output response interval influenced by these uncertainties. The results are compared with those obtained from Monte Carlo simulations (MCS), the first-order Taylor expansion method, and classical subinterval prediction methods. Our findings demonstrate close agreement between the outcomes derived from our proposed method and the exact solutions derived from the MCS simulations. Moreover, compared to a MCS, which requires function calls in order to obtain results, only 27 function calls are needed when using our proposed method. In contrast to conventional analysis methods like MCSs, not only does our proposed method exhibit higher accuracy, but it also demonstrates superior efficiency.