Next Article in Journal
Constrained Parameterized Differential Dynamic Programming for Waypoint-Trajectory Optimization
Previous Article in Journal
Numerical Investigation of the Vortex Ring Phenomena in Rotorcraft
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rigid–Flexible Coupling Dynamics Analysis of Coordination Arm and Application of a New Directional Subinterval Uncertainty Analysis Method

School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China
*
Author to whom correspondence should be addressed.
Aerospace 2024, 11(6), 419; https://doi.org/10.3390/aerospace11060419
Submission received: 28 March 2024 / Revised: 10 May 2024 / Accepted: 20 May 2024 / Published: 22 May 2024

Abstract

:
Cartridge delivery systems are commonly employed in aerospace engineering for the transportation of cylindrical projectiles. The coordination mechanism plays a pivotal role in ensuring reliable cartridge conveying, with its positioning accuracy being of utmost importance. However, accurately depicting the nonlinear relationship between input parameters and output response is challenging due to the involvement of numerous complex, uncertain factors during the movement process of the coordination mechanism. To address this issue, this study proposes a dynamics model that incorporates hinged gaps to represent rigid–flexible coupling within the coordination mechanism. Experimental validation confirms its effectiveness, while computational efficiency is enhanced through the utilization of a deep learning neural network surrogate model. Furthermore, an improved method for the uncertainty analysis of directional subintervals is introduced and applied to analyze uncertainty in coordination mechanisms, yielding results that demonstrate superior efficiency compared to other approaches.

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.

2. Dynamics Modeling of Coordinating Mechanism

A schematic diagram of the coordination mechanism structure of the coordination drug delivery device is shown in Figure 1 below. The coordination mechanism is driven by a hydraulic cylinder to rotate around the lug shaft, enabling the coordination function of the flipping and drug delivery mechanisms at various angles of attack. This complex and irregular mechanical structure primarily consists of rectangular tubular trusses welded together. Prior to dynamic modeling, it is necessary to simplify its structure and make reasonable assumptions. The rotary wheel hub serves as a pivotal component attached to the lug shaft, connected through high-precision bearings. For dynamic modeling purposes, it can be treated as a rigid body while disregarding any gaps between the lug shafts. The flexible coordination arm acts as the main load-bearing component, linked to the flipping mechanism via a rigid spline shaft and initially secured in position with pre-pressure. The drug delivery mechanism and flipping mechanism are firmly connected using bolts without any movement during coordination; thus, they can be simplified as concentrated masses affixed at the end of the flexible coordination arm. Pin shafts hinge both the upper fulcrums of the coordination cylinder and the flexible coordination arm, which are further connected to the drug delivery device frame through fixed pin shafts. Considering that uncertainties in processing technology may result in unavoidable hinge gaps, these gaps need thorough consideration during dynamics modeling.

2.1. Mathematical Model of Flexible Cantilever Beam

According to the schematic diagram of the coordination mechanism in Figure 1, the rigid–flexible coupling system of the coordination mechanism constitutes a typical class I rigid–flexible coupling problem [27], in which the coordination arm can be regarded as a flexible cantilever beam, and the overturning and dosing mechanisms can be regarded as fixed loads. Based on this, the dynamic modeling of the coordination arm is conducted, and its deformation model is shown in Figure 2. The coordinate system encompasses the inertial coordinate system, O X Y , and the floating coordinate system, o x y , established based on the flexible beam.
It can be seen from Figure 2 that the vector of the origin of the floating coordinate system with respect to the inertial coordinate system is r O . Take any point, P , on the flexible cantilever beam; its vector in the floating coordinate system before deformation is r , and its coordinates in the floating coordinate system are [ x , 0 ] T . After the deformation of the flexible cantilever beam, the point reaches point P , and then the vector of point P with respect to the inertial coordinate system can be expressed as follows:
r p = r O + C ( r + R R ) ,
where C is the coordinate transformation matrix of the floating coordinate system relative to the inertial coordinate system, R R is the deformation displacement vector from point P to point P before and after deformation, and its coordinates in the floating coordinate system are as follows:
r u = [ u ( x , t ) v ( x , t ) ] = [ w 1 ( x , t ) + w c ( x , t ) w 2 ( x , t ) ] = [ w 1 ( x , t ) 1 2 0 x ( w 2 ( ε , t ) ε ) 2 d ε w 2 ( x , t ) ] ,
where w 1 ( x , t ) is the x -axis directional extension in the floating coordinate system, w 2 ( x , t ) is the lateral bending deformation, and w c ( x , t ) = 1 2 0 x ( w 2 ( ε , t ) ε ) 2 d ε is the longitudinal deformation caused by the bending deformation of the flexible cantilever beam, which is called the nonlinear coupling deformation [1]. When the flexible cantilever beam works at a relatively high speed, the influence of u c on its dynamic characteristics cannot be ignored.
The finite element method and the hypothetical mode method are usually adopted to solve the vibration model of a flexible cantilever beam. Because the dynamics model established by the finite element method has many degrees of freedom, which is not conducive to active control design, the structural control of the whole system is actually realized by controlling several major modes. Therefore, the discretization method of the hypothetical mode method is more convenient for subsequent control research to establish the dynamics model. The specific expression for the deformation of the flexible cantilever beam is as follows:
{ w 1 ( x , t ) = N 1 ( x ) q 1 ( t ) w 2 ( x , t ) = N 2 ( x ) q 2 ( t ) ,
where N 1 ( x ) and N 2 ( x ) are, respectively, the eigenfunction row vectors of longitudinal and transverse vibrations and q 1 ( t ) and q 2 ( t ) are, respectively, the eigenfunction column vectors of longitudinal and transverse vibrations, which can be expressed as follows:
{ N i ( x ) = ( N 1 ( i ) ( x ) , N 2 ( i ) ( x ) , N n ( i ) ( x ) ) q i = ( q 1 ( i ) , q 2 ( i ) , q n ( i ) ) T ,
where, for i = 1 , 2 , n is the number of hypothetical modes. The eigenfunction can be expressed as follows:
N n ( 1 ) ( x ) = sin ( 2 j 1 ) 2 L x ,   j = 1 , 2 n N n ( 2 ) ( x ) = cos β j x cosh β j x + γ j ( sin β j x sinh β j x ) ,   j = 1 , 2 n γ i = ( cos β i x + cosh β i x ) ( sin β i x + sinh β i x )
Then, the deformation displacement vector, R R , from point P to point P before and after deformation can be expressed as follows:
R R = [ w 1 1 2 0 x ( w 2 ( ε , t ) ε ) 2 d ε w 2 ] = [ N 1 ( x ) q 1 ( t ) 1 2 q 2 T S ( x ) q 2 ( t ) N 2 ( x ) q 2 ( t ) ]
In summary, the generalized coordinates of a flexible cantilever beam can be expressed as follows:
Q f = [ r O , θ , q 1 T , q 2 T ] T ,
where r O and θ are the reference coordinates of the floating coordinate system. According to the second type of Lagrange equations, the rigid–flexible coupling dynamic equation of the flexible cantilever beam can be derived as follows:
M Q ¨ f + C Q ˙ f + K Q f = F ,
where M is the generalized mass matrix, C is the generalized damping matrix, K is the generalized stiffness matrix, and F is the generalized force matrix.

2.2. Mathematical Model of Contact Force in Clearance of Shaft–Hole Joint

When the two bodies represented by the shaft and hub undergo a curved surface contact collision, the contact collision force, F C , can be expressed by the combined force of the normal contact force, F N , and the tangential contact force, F T :
F C = F N + F T
In the classical dynamic modeling of contact collision force, the Hertz contact model is usually used to represent the contact of shaft–hole articulated clearance, but this model assumes that the contact collision between the shaft and the hub is completely elastic, ignoring the energy loss in the process. However, in practical engineering problems, more factors should be considered in the contact collision force, such as the relative contact velocity, the geometric parameters of the shaft and hub, material properties, etc. In this paper, a mixed contact force model is adopted to construct the normal contact collision force of shaft–hole articulated clearance [27], which can be expressed as follows:
F N = K n δ n + D m o d δ ˙ ,
where D m o d = 3 K n ( 1 - c r 2 ) e 2 ( 1 - c r ) δ n 4 δ ˙ - is the improved nonlinear damping coefficient. K n = 1 8 π E 2 δ ( 3 ( R h - R a ) + 2 δ ) 2 ( R h - R a + δ ) 3 is the improved nonlinear stiffness coefficient in the mixed contact force model, where E = E a ( 1 v h 2 ) + E h ( 1 v a 2 ) E h E a is the composite elastic modulus, E a and E h are the elastic moduli of the shaft and hub, v a and v h are the Poisson’s ratios of the shaft and hub, δ is the initial collision velocity, c r is the recovery coefficient, and n is the exponential coefficient.
The tangential contact force of shaft–hole articulated clearance is usually replaced by a tangential friction force. In this paper, the Ambrosio coulomb friction model is adopted, which is a continuous function of the tangential relative velocity and can guarantee the stability and continuity of numerical integration [27]. Its expression form is shown as follows:
F T = f k c d F N sgn ( V t ) ,
where f k is the friction coefficient and c d is the dynamic correction coefficient related to tangential velocity, which can be expressed as follows:
c d = { 0 V t v 0 V t v 0 v 1 v 0     v 0 V t v 1 1 V t v 1 ,
where v 0 and v 1 are the boundary values of the given tangential velocity interval.

2.3. Experimental Study and Model Verification of Coordination Mechanism

In order to verify the accuracy and effectiveness of the dynamics model of the coordinating mechanism, this section will conduct a comparative verification of the dynamics model based on a principle prototype platform built in the laboratory. The principle prototype platform is shown in Figure 3. The test process of the coordination mechanism is simplified as follows: the hydraulic oil source provides power to the coordination driving cylinder, while the PLC control box regulates the opening degree of the hydraulic servo valve based on instructions, thereby enabling the completion of coordinated actions by driving the coordination cylinder. During this process, pressure sensors installed on the coordination driving cylinder measure and record data regarding both rodless cavity and rod cavity pressures. These measured values serve as inputs for a dynamic simulation model of the coordination mechanism. Additionally, an angle encoder is positioned at the hinge connecting the coordination arm and ear shaft to monitor real-time changes in angle relative to their initial positions. This information serves as output for evaluating the actual conditions within the coordination mechanism. All aforementioned data are recorded in real time using a data acquisition instrument for subsequent model comparison and verification.
This paper mainly conducts tests and records for the condition that the coordination mechanism is fully charged and the coordination arm is coordinated from the 3° drug contact position on the horizontal surface to the 60° angle position below the horizontal surface, namely a complete coordination angle of 63°. As shown in Figure 4, Figure 5 and Figure 6 below, Figure 4 and Figure 5 are the pressure smooth curves of the rodless cavity and the rod-loaded cavity of the coordination driving cylinder measured by the pressure sensor in real time during the coordination mechanism’s coordination to the 60° angle action, respectively, and Figure 6 is the cylinder force curve obtained by converting the pressure data of the rodless cavity and the rod-loaded cavity of the coordination driving cylinder.
In order to verify the validity of the dynamics simulation model, Test Data 1 were selected as the dynamic input into the dynamics simulation model of the coordinating mechanism, and the angular displacement of the coordinating mechanism coordinated to a 60° incident angle under the full charge state was compared with the real-time measured results of the angle encoder when the coordinating mechanism was coordinated to a 60° incident angle under the actual working conditions, as shown in Figure 7. Through comparative analysis, the simulation results were found to be close to the test results, indicating the rationality and validity of the dynamics simulation model of the coordinating mechanism, and it can be considered that the simulation model can be applied to the subsequent dynamics characteristic analysis, uncertainty analysis, and optimization.

3. Construction of Coordination Agency Agent Model

The structure of the ammunition automatic loading system is complex and has obvious structural characteristics. The structural dynamics model takes into account a variety of factors, resulting in a strong nonlinear relationship between input and output. In the process of the uncertainty analysis and optimization of the coordination mechanism, the dynamics model needs to be called upon for calculation many times. With an improvement in the accuracy of the dynamics model, the cost of a single calculation increases significantly. Directly calling upon the original dynamics model will lead to unbearable calculation costs. Therefore, it is a common solution to use a high-precision surrogate model instead of the dynamics simulation model. The main advantage of the surrogate model method is the ability to use an implicit approximation function to approximately replace the original dynamics model so as to carry out calculation analysis or optimization and reduce calculation costs. In this section, a surrogate model of the coordination mechanism is constructed based on a deep neural network (DNN).

3.1. Deep Neural Network Architecture

The neural network method is widely used in surrogate model research due to its excellent nonlinear fitting ability. It is usually composed of a multilayer system, consisting of an input layer, an output layer, and a hidden layer, which can transform input information into useful output information. In research on multi-hidden neural network algorithms, many scholars have found that with an increase in the number of hidden layers, the ability of neural network models to analyze complex patterns is also improved. Multi-hidden neural networks are actually deep versions of shallow neural networks which attempt to express features by using more neurons, which is typical of deep neural networks (DNNs). Compared with the classical structure of artificial neural network models, the main difference between these and DNNs is that they have more hidden layers. As long as there are more than three hidden layers, it can be defined as a DNN, and its basic architecture is shown in Figure 8. DNNs can be described as having a layout consisting of multilayer neurons, which are interconnected by weights and adjusted iteratively through an optimization process. The initial data are input into the input layer, and then the input layer distributes the data to each neuron in the first hidden layer. The output of the first hidden layer is then passed to the next hidden layer until it reaches the output layer, where data are distributed to the output neurons, which are equal in number to the number of categories. During training, the DNN uses a backpropagation algorithm. It compares the network output with the expected output to calculate the error and propagates the error back. The random gradient descent method, combined with a loss function, is used to update the weights. This process is iterated and repeated, updating all training data on the entire network in a training round. The single-neuron model of a DNN is shown in Figure 9.
The construction of a deep neural network (DNN) agent model typically involves two stages: forward propagation and backward propagation. During the forward propagation stage, the distance from the input layer to the output layer is computed using training samples in order to estimate the output result of the test samples. Subsequently, network parameters are adjusted to minimize the discrepancy between the calculated outputs for training samples and their target values. In the backward propagation stage, the gradient descent algorithm is employed to optimize the loss function and determine weight matrices and bias terms that yield minimal loss values as optimal parameters. The mean-squared error function (MSE) serves as the utilized loss function in this study:
J ( W , B , x , y ) = 1 n i = 1 n y i f ( x i ) 2
In this equation, W and B represent the input weight matrix and bias matrix, respectively, while y i and f ( x i ) represent the theoretical value and predicted value of the output of layer i , respectively. represents the Euclidean norm. For each layer, W and B satisfy the following relationship:
f ( x i ) = σ ( z i ) = σ ( W i f ( x i 1 ) + B i )
where z i is the input of layer i , and σ ( ) is the activation function, which is adopted as the ReLU activation function, whose nonlinear ability can make the DNN perform better in nonlinear regression, not only accelerating the convergence speed but also improving the model’s predictive ability. Then, the loss function in Equation (13) can be rewritten as follows:
J ( W , B , x , y ) = 1 n i = 1 n y i σ ( W i f ( x i 1 ) + B i ) 2
From the above equation, we can obtain the gradients of W and B as follows:
J ( W , B , x , y ) W i = [ ( y i f ( x i ) ) σ ( z i ) ] [ f ( x i 1 ) ] T
J ( W , B , x , y ) B i = ( y i f ( x i ) ) σ ( z i )
Then, the gradient for any layer j is as follows:
δ j = J ( W , B , x , y ) z j = ( z i z i 1 z i 1 z i 2 z j + 1 z j ) J ( W , B , x , y ) z i
According to the forward propagation algorithm, W j and B j satisfy the following relationship:
z j = W j f ( x j 1 ) + B j
Then, the gradients of W j and B j in layer j are as follows:
J ( W , B , x , y ) W j = δ j [ f ( x j 1 ) ] T
J ( W , B , x , y ) B j = δ j
From the recursive formula above, we can derive the following:
δ j = ( W j + 1 ) T δ j + 1 σ ( z j )
If the hidden layer j contains m neurons, then the updates for W j and B j of layer j are represented as follows:
W j = W j α p = 1 m δ p , j [ f ( x p , j 1 ) ] T
B j = B j α p = 1 m δ p , j
In this formula, α represents the learning efficiency.
During the forward propagation process of constructing a DNN, the weights and biases of each layer’s neuron nodes are initialized randomly, resulting in a stochastic prediction, f ( x i ) . Subsequently, employing the backpropagation algorithm, the neural network parameters undergo fine-tuning based on the loss function, Equation (13), to progressively converge towards their theoretical values. After multiple iterations of weight and bias updates, the trained DNN surrogate model is obtained.

3.2. Construction of Proxy Model for Coordination Mechanism Based on DNN

According to the construction method of the DNN surrogate model, training samples for training the DNN surrogate model are first obtained, and multiple values are extracted and combined in the subintervals of each uncertain input parameter interval to form the input samples for training the DNN surrogate model. The input samples are taken into the coordination mechanism dynamics model for solving to obtain corresponding output responses to form the training set, validation set, and test set of the surrogate model, and then the input samples and output responses are used for training the DNN surrogate model. Among them, 1000 groups of input samples are selected, and the sampling intervals of uncertain input parameters of the input samples are shown in Table 1.
The regression performance of the coordinating agency DNN surrogate model on the training set, validation set, and test set is shown in Figure 10. As can be seen from the training results, the complex correlation coefficients, R , of the training set, validation set, and test set are all above 0.9, indicating the accuracy and effectiveness of the trained coordinating agency DNN surrogate model, which can be used for subsequent uncertainty analysis and interval parameter optimization.

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, X = ( X 1 , X 2 ) , interval analysis problem as an example, for convenience, assume that both input parameter intervals are divided into m subintervals and stipulate that m is 2 j , where j 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, X 0 I , in the subinterval domain of the input vector X , and according to the classical subinterval prediction method, conduct a first-order Taylor expansion at the midpoint, X 0 c , of the subinterval with an amplitude of Δ X / m :
f ( X 0 I ) = f ( X 0 c ) + f ( X ) X | X = X 0 c Δ X m ,
where X 0 I = [ X 0 c X r / m , X 0 c + X r / m ] . The upper and lower bounds of f ( X 0 I ) can be expressed as follows:
f l ( X 0 I ) = f ( X 0 c ) | f ( X ) X | X = X 0 c | X r m
f u ( X 0 I ) = f ( X 0 c ) + | f ( X ) X | X = X 0 c | X r m
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 f ( X ) by applying first-order Taylor expansions at multiple points on the path. When predicting the upper bound of f ( X ) , the expansion path of extension points should move in the direction of the increasing output response function. Therefore, according to the gradient of f ( X ) relative to X at X 0 c , the next expansion point, X u 1 c , required to predict the upper bound can be expressed as follows:
X u 1 c = X 0 c + s i g n ( f ( X ) X | X = X 0 c ) X r m
The first-order Taylor expansion of f ( X ) with an amplitude of Δ X / m at X u 1 c is as follows:
f ( X u 1 I ) = f ( X u 1 c ) + f ( X ) X | X = X u 1 c Δ X m
where the subinterval X u 1 I = [ X u 1 c X r / m , X u 1 c + X r / m ] . The upper and lower bounds of f ( X u 1 ) obtained from (26) and (27) can be expressed as follows:
f l ( X u 1 I ) = f ( X u 1 c ) | f ( X ) X | X = X u 1 c | X r m
f u ( X u 1 I ) = f ( X u 1 c ) + | f ( X ) X | X = X u 1 c | X r m
When the above steps are performed n 1 times in the same way, all the expansion points, X u i c , used to predict the upper bound can be uniformly expressed by the following equation:
X u i c = X u ( i 1 ) c + s i g n ( f ( X ) X | X = X u ( i 1 ) c ) X r m ,   i = 1 , 2 , , m 1
At the above point, the output response is a first-order Taylor expansion with an amplitude of Δ X / m , and the upper and lower bounds of all the expansion points’ output responses are obtained as follows:
f l ( X u i I ) = f ( X u i c ) | f ( X ) X | X = X u i c | X r m ,   i = 1 , 2 , , m 1
f u ( X u i I ) = f ( X u i c ) + | f ( X ) X | X = X u i c | X r m ,   i = 1 , 2 , , m 1
where the subinterval X u i I = [ X u i c X r / m , X u i c + X r / m ] , and the upper bound of f ( X u 1 I , X u 2 I , , X u i I ) can be obtained by interval fusion:
f u ( X u 1 I , X u 2 I , , X u i I ) = i = 1 , 2 , , m 1 f u ( X u i I ) = max i = 1 , 2 , , m 1 f u ( X u i I )
Similar to the method used for calculating the upper bound of the output response value interval, when predicting the lower bound of f ( X ) , the expansion path of the extension point should move along the direction of the decreasing output response function. Therefore, according to the gradient of f ( X ) relative to X at X 0 c , all the expansion points, X l i c , required for predicting the lower bound can be uniformly expressed as follows:
X l i c = X l ( i 1 ) c s i g n ( f ( X ) X | X = X l ( i 1 ) c ) X r m ,   i = 1 , 2 , , m 1
According to (33) and (34), the upper and lower bounds of the output responses of all the expansion points are as follows:
f l ( X l i I ) = f ( X l i c ) | f ( X ) X | X = X l i c | X r m ,   i = 1 , 2 , , m 1
f u ( X l i I ) = f ( X l i c ) + | f ( X ) X | X = X l i c | X r m ,   i = 1 , 2 , , m 1
where the subinterval X l i I = [ X l i c X r / m , X l i c + X r / m ] . According to Equation (35), the lower bound of f ( X l 1 I , X l 2 I , , X l i I ) can be obtained by interval fusion:
f l ( X l 1 I , X l 2 I , , X l i I ) = i = 1 , 2 , , m 1 f l ( X l i I ) = min i = 1 , 2 , , m 1 f l ( X l i I )
The directional subinterval prediction is based on gradient analysis by constructing two expansion paths from X 0 to X u i c and X l i c 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, f 0 l ( X I ) and f 0 u ( X I ) , 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, f ( X ) , 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 f ( X ) 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:
| f q u ( X I ) f q 1 u ( X I ) f q u ( X I ) | ε 1 | f q l ( X I ) f q 1 l ( X I ) f q l ( X I ) | ε 2 ,
where q is the number of adaptive iterations for subinterval division; f q l ( X I ) and f q u ( X I ) are the upper and lower bounds of the output response value interval calculated from q iterations; f q 1 l ( X I ) and f q 1 u ( X I ) are the upper and lower bounds calculated through ( q 1 ) iterations; and ε 1 and ε 2 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, f ( X I ) . Let q = 0 .
Step 3: Find the appropriate expansion path for predicting the upper and lower bounds of the output response according to the gradient analysis of f ( X I ) .
Step 4: Calculate the upper and lower bounds of the value interval of f ( X u i I ) and f ( X l i I ) , respectively, according to Equations (33) and (34) and Equations (37) and (38). Predict the upper and lower bounds of the interval f ( X I ) according to Equations (35) and (39) through interval fusion.
Step 5: Test whether the upper and lower bounds of f ( X I ) obtained through the q 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 q = q + 1 .

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 ε 1 = ε 2 = 0.01 . 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 10 6 , 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 10 6 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.

Author Contributions

Writing—original draft, X.G.; Writing—review & editing, L.C. and J.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhou, K.; Wang, Z.; Gao, Q.; Yuan, S.; Tang, J. Recent advances in uncertainty quantification in structural response characterization and system identification. Probabilistic Eng. Mech. 2023, 74, 103507. [Google Scholar] [CrossRef]
  2. Elishakoff, I.; Ren, Y. The bird’s eye view on finite element method for structures with large stochastic variations. Comput. Methods Appl. Mech. Eng. 1999, 168, 51–61. [Google Scholar] [CrossRef]
  3. Ye, N.; Lu, Z. Single-loop cross-entropy-based directional importance sampling for efficient estimation of failure probability function. Aerosp. Sci. Technol. 2024, 146, 108962. [Google Scholar] [CrossRef]
  4. Ting, Y.; Zhenzhou, L. A novel single-loop Kriging importance sampling method for estimating failure probability upper bound under random-interval mixed uncertainties. Aerosp. Sci. Technol. 2023, 143, 108722. [Google Scholar]
  5. Liao, Z.; Xia, W.; He, X. An investigation into Markov chain Monte Carlo algorithms for Subset simulation: Emphasizing uncertainty analysis. Comput. Struct. 2024, 294, 107268. [Google Scholar] [CrossRef]
  6. Prabhu, S.; Ehrett, C.; Javanbarg, M.; Brown, D.A.; Lehmann, M.; Atamturktur, S. Uncertainty Quantification in Fault Tree Analysis: Estimating Business Interruption due to Seismic Hazard. Nat. Hazards Rev. 2020, 21, 04020015. [Google Scholar] [CrossRef]
  7. Saberi, S.; Abdollahi, A.; Friswell, M.I. Probability analysis of bistable composite laminates using the subset simulation method. Compos. Struct. 2021, 271, 114120. [Google Scholar] [CrossRef]
  8. Poikonen, T.; Blattner, P.; Kärhä, P.; Ikonen, E. Uncertainty analysis of photometer directional response index f SUB2/SUB using Monte Carlo simulation. Metrologia 2012, 49, 727–736. [Google Scholar] [CrossRef]
  9. Ben-Haim, Y.; Elishakoff, I. Convex Models of Uncertainty in Applied Mechanics; Elsevier: Amsterdam, The Netherlands, 2013. [Google Scholar]
  10. Zhou, Y.; Huang, G.; Wang, S.; Zhai, Y.; Xin, X. Water resources management under dual uncertainties: A factorial fuzzy two-stage stochastic programming approach. Stoch. Environ. Res. Risk Assess. 2016, 30, 795–811. [Google Scholar] [CrossRef]
  11. Qiu, Z.; Lv, Z. The vertex solution theorem and its coupled framework for static analysis of structures with interval parameters. Int. J. Numer. Methods Eng. 2017, 112, 711–736. [Google Scholar] [CrossRef]
  12. Ma, Y.; Liang, Z.; Chen, M.; Hong, J. Interval analysis of rotor dynamic response with uncertain parameters. J. Sound Vib. 2013, 332, 3869–3880. [Google Scholar] [CrossRef]
  13. Liu, Y.S.; Wang, X.J.; Wang, L. Interval uncertainty analysis for static response of structures using radial basis functions. Appl. Math. Model. 2019, 69, 425–440. [Google Scholar] [CrossRef]
  14. Liu, Y.S.; Wang, X.J.; Wang, L.; Zheng, L. A Bayesian collocation method for static analysis of structures with unknown-but-bounded uncertainties. Comput. Methods Appl. Mech. Eng. 2019, 346, 727–745. [Google Scholar] [CrossRef]
  15. Wang, L.; Chen, Z.; Yang, G. An interval uncertainty analysis method for structural response bounds using feedforward neural network differentiation. Appl. Math. Model. 2020, 82, 449–468. [Google Scholar] [CrossRef]
  16. Wei, T.H.; Meng, G.W.; Li, F.; Cong, Y.B. Interval analysis of structure based on dimension reduction method and Chebyshev polynomial. Huazhong Univ. Sci. Technol. (Nat. Sci. Ed.) 2021, 41, 14–19. [Google Scholar]
  17. Chen, S.H.; Lian, H.D.; Yang, X.W. Interval static displacement analysis for structures with interval parameters. Int. J. Numer. Methods Eng. 2002, 53, 393–407. [Google Scholar] [CrossRef]
  18. Ma, K.; Li, B.H.; Yang, K.; Liu, Q.L. First and second-order sensitivity method of structure static displacement. J. Jilin Univ. (Eng. Technol. Ed.) 2021, 51, 472–477. [Google Scholar]
  19. Li, J.P.; Chen, J.J.; Zhu, Z.Q.; Li, J.P.; Chen, J.J.; Zhu, Z.Q. A Method for Solving the Structural Interval Finite Element Equations. Eng. Mech. 2010, 27, 79–83. [Google Scholar]
  20. Wu, P.G.; Ni, B.Y.; Jiang, C. An Interval Finite Element Method Based on the Neumann Series Expansion. Chin. J. Theor. Appl. Mech. 2020, 52, 1431–1442. [Google Scholar]
  21. Qiu, Z.; Elishakoff, I. Antioptimization of structures with large uncertain-but-non-random arameters via interval analysis. Comput. Methods Appl. Mech. Eng. 1998, 152, 361–372. [Google Scholar] [CrossRef]
  22. Zhou, Y.T.; Jiang, C.; Han, X. Interval and subinterval analysis methods of the structural analysis and their error estimations. Int. J. Comput. Methods 2006, 3, 229–244. [Google Scholar] [CrossRef]
  23. Chen, N.; Yu, D.; Xia, B.; Liu, J.; Ma, Z. Interval and subinterval homogenization-based method for determining the effective elastic roperties of periodic microstructure with interval parameters. Int. J. Solids Struct. 2017, 106–107, 174–182. [Google Scholar] [CrossRef]
  24. Fu, C.M.; Cao, L.X.; Tang, J.C.; Long, X.Y. A subinterval decomposition analysis method for uncertain structures with large uncertainty parameters. Comput. Struct. 2018, 197, 58–69. [Google Scholar] [CrossRef]
  25. Aloisio, A.; Contento, A.; Alaggio, R.; Quaranta, G. Physics-based models, surrogate models and experimental assessment of the vehicle–bridge interaction in braking conditions. Mech. Syst. Signal Process. 2023, 194, 110276. [Google Scholar] [CrossRef]
  26. Cunha, B.Z.; Droz, C.; Zine, A.M.; Foulard, S.; Ichchou, M. A review of machine learning methods applied to structural dynamics and vibroacoustic. Mech. Syst. Signal Process. 2023, 200, 110535. [Google Scholar] [CrossRef]
  27. Chen, Z.Q.; Qian, L.F.; Chen, G.S.; Nie, S.C.; Yin, Q. Dynamics of luffing motion of a hydraulically driven shell manipulator with revolute clearance joints. Def. Technol. 2022, 18, 689–708. [Google Scholar]
Figure 1. Schematic diagram of coordination mechanism.
Figure 1. Schematic diagram of coordination mechanism.
Aerospace 11 00419 g001
Figure 2. Deformation of flexible cantilever beam.
Figure 2. Deformation of flexible cantilever beam.
Aerospace 11 00419 g002
Figure 3. Principle prototype platform composition diagram.
Figure 3. Principle prototype platform composition diagram.
Aerospace 11 00419 g003
Figure 4. Coordination cylinder’s rodless chamber pressure test data curve.
Figure 4. Coordination cylinder’s rodless chamber pressure test data curve.
Aerospace 11 00419 g004
Figure 5. Coordinated cylinder’s pressure test data curve with rod cavity.
Figure 5. Coordinated cylinder’s pressure test data curve with rod cavity.
Aerospace 11 00419 g005
Figure 6. Coordinated cylinder’s resultant pressure test data curve.
Figure 6. Coordinated cylinder’s resultant pressure test data curve.
Aerospace 11 00419 g006
Figure 7. Comparison of angular displacement test data and simulation data of coordinate mechanism.
Figure 7. Comparison of angular displacement test data and simulation data of coordinate mechanism.
Aerospace 11 00419 g007
Figure 8. DNN architecture.
Figure 8. DNN architecture.
Aerospace 11 00419 g008
Figure 9. Single-neuron model.
Figure 9. Single-neuron model.
Aerospace 11 00419 g009
Figure 10. Performance of DNN proxy model for coordination agencies.
Figure 10. Performance of DNN proxy model for coordination agencies.
Aerospace 11 00419 g010
Figure 11. Schematic diagram of adaptive strategy for subinterval partition.
Figure 11. Schematic diagram of adaptive strategy for subinterval partition.
Aerospace 11 00419 g011
Figure 12. Main flow chart of directional subinterval prediction method.
Figure 12. Main flow chart of directional subinterval prediction method.
Aerospace 11 00419 g012
Figure 13. Output response boundary and error comparison of coordination mechanisms obtained by directional subinterval prediction method and MCS method.
Figure 13. Output response boundary and error comparison of coordination mechanisms obtained by directional subinterval prediction method and MCS method.
Aerospace 11 00419 g013
Table 1. Input parameter training set sampling interval.
Table 1. Input parameter training set sampling interval.
SymbolUncertain ParametersValue RangeTheoretical Design Value
M X T Loading (without drug) mass of drug-coordinating arm/kg[101.14, 105.73]102.06
m y Module drug weight (6 blocks)/kg[20.2, 20.4]20.24
R 1 Upper cylinder fulcrum pin radius/mm[12.4, 12.5]12.5
R 2 Upper cylinder fulcrum hole radius/mm[12.5, 12.6]12.5
R 3 Radius of coordinator fulcrum hole/mm[12.5, 12.6]12.5
R 4 Lower fulcrum shaft radius of cylinder/mm[14.9, 15]15
R 5 Lower fulcrum hole radius of cylinder/mm[15, 15.1]15
μ Coefficient of friction[0.25, 0.35]0.3
Table 2. Uncertainty input parameter interval and uncertainty level.
Table 2. Uncertainty input parameter interval and uncertainty level.
IndexUncertain ParametersValue RangeTheoretical Design Value
1Loading (without drug) mass of drug-coordinating arm/kg[101.14, 105.73]0.0222
2Module drug weight (6 blocks)/kg[20.2, 20.4]0.0049
3Upper cylinder fulcrum pin radius/mm[12.4, 12.5]0.004
4Upper cylinder fulcrum hole radius/mm[12.5, 12.6]0.004
5Radius of coordinator fulcrum hole/mm[12.5, 12.6]0.004
6Lower fulcrum shaft radius of cylinder/mm[14.9, 15]0.0033
7Lower fulcrum hole radius of cylinder/mm[15, 15.1]0.0033
8Coefficient of friction[0.25, 0.35]0.1667
Table 3. Upper and lower boundaries of output responses of coordination mechanisms obtained by different analytical methods.
Table 3. Upper and lower boundaries of output responses of coordination mechanisms obtained by different analytical methods.
MethodsLower Bound of Output ResponseUpper Bound of Output ResponseNumber of Calls
ValueErrorValueError
MCS26.291727.2592 10 6
First-order Taylor expansion26.264910.19‱27.308918.23‱9
Classical subinterval prediction26.27376.85‱27.27495.76‱2304
Directional subinterval prediction26.28980.72‱27.25621.10‱27
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gao, X.; Chen, L.; Tang, J. Rigid–Flexible Coupling Dynamics Analysis of Coordination Arm and Application of a New Directional Subinterval Uncertainty Analysis Method. Aerospace 2024, 11, 419. https://doi.org/10.3390/aerospace11060419

AMA Style

Gao X, Chen L, Tang J. Rigid–Flexible Coupling Dynamics Analysis of Coordination Arm and Application of a New Directional Subinterval Uncertainty Analysis Method. Aerospace. 2024; 11(6):419. https://doi.org/10.3390/aerospace11060419

Chicago/Turabian Style

Gao, Xuan, Longmiao Chen, and Jingsong Tang. 2024. "Rigid–Flexible Coupling Dynamics Analysis of Coordination Arm and Application of a New Directional Subinterval Uncertainty Analysis Method" Aerospace 11, no. 6: 419. https://doi.org/10.3390/aerospace11060419

APA Style

Gao, X., Chen, L., & Tang, J. (2024). Rigid–Flexible Coupling Dynamics Analysis of Coordination Arm and Application of a New Directional Subinterval Uncertainty Analysis Method. Aerospace, 11(6), 419. https://doi.org/10.3390/aerospace11060419

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