Analysis of Global Sensitivity of Landing Variables on Landing Loads and Extreme Values of the Loads in Carrier-Based Aircrafts

When a carrier-based aircraft is in arrested landing on deck, the impact loads on landing gears and airframe are closely related to landing states. The distribution and extreme values of the landing loads obtained during life-cycle analysis provide an important basis for buffering parameter design and fatigue design. In this paper, the effect of the multivariate distribution was studied based on military standards and guides. By establishment of a virtual prototype, the extended Fourier amplitude sensitivity test (EFAST) method is applied on sensitivity analysis of landing variables. The results show that sinking speed and rolling angle are the main influencing factors on the landing gear’s course load and vertical load; sinking speed, rolling angle, and yawing angle are the main influencing factors on the landing gear’s lateral load; and sinking speed is the main influencing factor on the barycenter overload. The extreme values of loads show that the typical condition design in the structural strength analysis is safe. The maximum difference value of the vertical load of the main landing gear is 12.0%. This research may provide some reference for structure design of landing gears and compilation of load spectrum for carrier-based aircrafts.


Introduction
Carrier-based aircraft is the most important method for aircraft carrier strike group to control sea supremacy, and it is also an indispensable power for modern navies.Due to limited area in deck landing zone and the demand for bolting and go-around, carrier-based aircraft usually lands on deck via impact method under high sinking speed and high engaging speed along a fixed glide path angle [1].The impact load, braking load of arresting cable [2,3], and other loads at the moment when the aircraft touches deck put forward higher requirements for design and analysis of landing gears and airframe structure, especially for the structures closely related to landing [4].
At the primary design stage of an aircraft landing gear, the landing loads and the most severe landing conditions of the landing gear under different landing conditions should be determined according to design outline, which will be adopted in the parameter design of landing gear buffers and structural strength analysis.At the detailed design stage, optimization design will be carried out to balance performance and structure of the landing gear according to the relation between landing conditions and landing loads.Finally, the maneuvering envelope will be determined and the design of service load spectrum will be compiled for the assessment of fatigue life [5,6].
At present, the research literatures and reports on definition of arrested deck landing conditions are mainly military standards and guides [7,8].MIL-A-8863C (AS) provides the landing variables, their distribution form, and empirical formula of mean value and standard deviation that should be considered in deck landing.Micklos [9] provided a measurement report on the landing variables of all kinds of carrier-based aircrafts landed on the Enterprise Aircraft Carrier in both day and night.The measurement of landing variables and loads is a costly and tedious task, thus a large amount of researches have been carried out on the simulation and analysis of landing dynamics to analyze different landing conditions and loads in recent years.Zhang [10] introduced deck motion into a dynamics model through the wheel-deck coordinate system and simulated three landing situations to verify the model.Mikhaluk [11] and Lihua et al. [12] used finite element model to analyze forces on cables under different initial velocities and masses.Sati et al. [13] built a detailed aircraft arresting system using bond graph approach and studied the landing performance at different engagement speeds and masses.
The studies mentioned above give a simple trend of the relationship between landing variables and loads under a few landing conditions.However, the conditions empirically selected are unrepresentative.During the compilation of design service load spectrum, representative values of typical landing conditions can be obtained by discretizing the interval of landing parameters through sensitivity analysis.Incorrect structural load analysis may lead to excessive design, which would induce large weight and high cost at the primary stage of structure design.A further determination of the relationship between landing variables as well as the extreme value of the loads is required.Sebastian [14] analyzed the condition of free flight engagement (FFE), and the results show that the most severe loading condition of the landing gear takes place under the condition of free flight engagement at small sinking speed.Chester [15] researched the response of the main and nose gears via simulation of landing impact considering pitching and heaving degrees of freedom of the aircraft motion.The simulation indicates that the maximum vertical loads of main gears are almost linearly dependent on the sinking speed, while the response of nose gear is very sensitive to the initial values of pitch angle and pitching inertia.Yunwen et al. [16] explored the effect of different landing variables on sinking velocity based on the landing data measured from an E-2C.The results show that aircraft path angle and deck pitch angle are highly correlated with sinking velocity.Although they have studied part of the relevance between landing variables and loads, there is no direct conclusion on method accurate enough for determination of the representative values that can reflect typical landing conditions.Therefore, sensitivity and extreme value analysis need to be studied and it provides an important basis for the parameter design of landing gear buffer, structural optimization design [17], and fatigue analysis.However, theoretical basis and analysis method still need to be further investigated.
The extended Fourier amplitude sensitive test (EFAST) analysis method [18] is adopted in this paper to study the coupling and sensitivity between landing variables and landing loads.The method has been well applied in the field of hydrology [19], physical model [20], and so on.Based on military standards and guides, the effects of the multivariate distribution on the probability density function of single landing variable were analyzed and the distribution of single landing variable was fitted.By using the simplified landing virtual prototype and analysis of a large amount of landing conditions with the multicondition automatic simulation technology, the first-order and global sensitivity coefficients of landing variables on landing loads were calculated quantitatively, and the extreme value conditions and the frequency curves of the landing load were obtained.

Studies on Landing Variables in Military Standards
Many factors would affect the carrier-based aircrafts during the arrested deck landing.The landing variables and their relationships are mainly stipulated by related military standards and specifications.
2.1.Landing Variables and Multivariate Distribution.The military standard MIL-A-8863C (AS) stipulates the landing variables and their distributions.The condition of the landing variables needs to satisfy the demands of the multivariate distributions.The joint probability of the eight landing variables P T is calculated by where V TD is the approaching speed, V E is the engaging speed, V V is the sinking speed, θ P is the aircraft pitching angle, θ R is the aircraft rolling angle, θ R is the aircraft rolling rate, θ Y is the aircraft yawing angle, and d is the off-center arresting distance.The subscript i denotes the initial value of each landing variable.The symbol / < is chosen to be either > or < according to the initial value of the variable.
When the initial value of any landing variable is greater than the average one, the symbol > is chosen and it represents that the cumulative probability of that landing variable is greater than the initial value.Conversely, when the symbol < is chosen, it represents that the cumulative probability of that landing variable is less than the initial value.The evaluation of the landing variables is determined by the extreme value of the cumulative probability P 0 , namely, where x max i and x min i are the maximum and minimum value of the variables.In this paper, only the touch-and-go and arrested landing conditions are considered and the landing variables follow a normal distribution.The cumulative probability P 0 = 1 × 10 −3 , and the multivariate joint probability P T = 7 8125 × 10 −6 2.2.The Effects of the Multivariate Distribution on the Landing Variables.The relationships among landing variables are according to the multivariate distribution.The product of the joint probability of the landing variables does not reflect the effects on the distribution of the single landing variable directly.Therefore, the distribution of landing variables is investigated according to the multivariate distribution.Take the approaching speed V TD as an example.There exists a relationship: where μ V TD is the average of the carrier engaging speeds.Take V TD and V E as an example.The extreme value of the 2 International Journal of Aerospace Engineering joint probability of any two landing variables is given by inserting (3) into (1).
The distribution of the landing variable is mapped to a standard normal distribution in order to unify the range of the variables.The mutual effects among the value ranges of the landing variables are determined by the multivariate distribution [21].In Figure 1, there are six equiprobable design envelopes of any two landing variables.The probability of each equiprobable envelope is shown in the figure.
The equiprobable envelope 6 represents the condition that the joint probability of two landing variables reaches the minimum, the value of which is P T ′ = 0 0005 The equiprobable envelope 1 degenerates to a point which represents the condition that the joint probability of two landing variables reaches the maximum.
The evaluation of the eight landing variables of four groups on the equiprobable envelopes must ensure that the sum of their sequence number is nine, satisfying the demand of the joint probability.For example, when one group of variables is the equiprobable envelope 5, then other three groups must be the equiprobable envelope 1, 1 and 2.
It is difficult to make analytical calculations on the effects of the multivariate distribution on the single landing variable.Monte Carlo method is adopted to sample simulation for the complex probability constraints.Before sampling, the distribution of each variable should be mapped to the standard normal distribution.10 5 simulate samples are extracted according to the distribution of each landing variable.The simulating sample screening is according to the multivariate distribution determined by (1).Figures 2 and 3 show the Monte Carlo sampling condition and the condition of applying multivariate distribution.
From the screening results, it can be seen that only 29,194 out of 10 5 samples remain after adopting multivariate distribution and the 71% of the samples cannot satisfy the requirements of the joint probability.By comparing Figure 2 with Figure 3, it can be known that the multivariate distribution deletes lots of samples that landing variables approaching the extreme value and keeps the samples that landing variables closing to the average.In Figure 3, the equiprobable envelope 6 is the envelope curve of the minimum value of the multivariate distribution shown in Figure 1.
The normal test is carried out by Kolmogorov-Smirnov method.D = 0 005, P = 0 372 > 0 05, and the landing variables after screening obey the normal distribution.The comparisons between distributions of landing variables before and after screening as well as the fitting conditions are shown in Figure 4.It is seen that after adopting the multivariate distribution, the probability density function of the single landing variable changes a lot.
From the analysis above, it is seen that the requirement of the multivariate distribution restricts the optional values   3 International Journal of Aerospace Engineering among the variables and has a large effect on the standard deviation of the single landing variable.The results of the sensitivity analysis are closely related to the probability density function of the input variables.The adoption of the distribution of the landing variables before their multivariate distribution would overestimate the probability that landing variables reach the maximum or minimum values, which may cause the errors of the global sensitivity analysis of the landing load.In this section, based on the requirements of the multivariate distribution, the equiprobable design envelopes of the landing variables are obtained.The distribution of the single landing variable for global sensitivity analysis is studied by Monte Carlo sampling method.

The Landing Virtual Prototype of the Carrier-Based Aircraft and the Multicondition Automatic Simulation
The attitudes of the aircrafts differ a lot at the moment of deck landing.Under the effects of the deck reacting force, the arresting force of the arresting hook, the aerodynamic force, and other forces, the load-time history of the landing gear is rather complicated during the landing process.The establishment of the accurate landing simulation model of the aircraft is the basis of the sensitivity analysis of the landing load.The accuracy of the sensitivity analysis is related with the sample quantity.In order to improve the reliability of the sensitivity analysis, it is necessary to consider the computation efficiency during simulating different landing conditions.The constraints and loads of the aircrafts during landing are analyzed, and the landing virtual prototype of the carrier-based aircraft based on LMS Virtual.Lab is established.For the computation efficiency problem of the massive landing conditions, the secondary development technology is applied to realize the multicondition automatic simulation and data processing of the landing conditions.

Hypothesis of the Aircraft Virtual Prototype and Model
Establishment.The aircraft virtual prototype consists of nose landing gear subsystem, main landing gear subsystem, arresting hook system, airframe model, and deck model.The external forces on the aircraft during landing are body load which includes aircraft gravity and engine thrust, aerodynamic force, acting force between tires, and deck and arresting force.
The virtual prototype has the following reasonable assumptions: (1) Ignore the motion of the aircraft carrier and offcenter landing.
(2) Assume that all components of the subsystems are rigid body, and the mass concentrates on the centroid.
(3) Ignore the effects of yawing landing on the arresting force curve.
According to the real motion of each component of the subsystems during the aircraft landing, the relating motion pairs such as swirl joint, cylinder joint, and spherical joint in subsystems are established.The constraints and load conditions of the subsystems are shown in Figure 5.
The initial condition of the virtual prototype is the moment that the aircraft lands on the deck, and the landing gears just touch the deck.During landing, the engine thrust maintains the maximum (twin engine).Arresting forcedistance curve refers to military standards [8].MK7-3 arresting machine is selected.The curve is obtained by weight interpolation and velocity compensation, shown in Figure 6.
The virtual prototype is established according to E-2 airborne early warning aircraft, shown in Figure 7.

Landing Multicondition Automatic Simulation.
In order to improve the computation efficiency of the massive landing conditions, the automation technology is applied to realize the multicondition automatic simulation and data processing of the landing conditions.The auto-set, renewal, and analysis of the virtual prototype parameters are carried out by LMS    The virtual prototype of LMS Virtual.Lab motion can be divided into the product model and analysis model by func-tion.The secondary development of the virtual prototype computes the initial landing attitude and location of the aircraft.The flow chart of the multicondition automatic simulation of the virtual prototype is shown in Figure 8.
The initial landing attitude of the carrier-based aircraft is determined by the multivariable pitching angle θ P , rolling angle θ R , and yawing angle θ Y The attitude transition matrix from airframe coordinate system to ground coordinate system is calculated by Tait-Bryan angles [22].The transition matrix is used to update the aircraft attitude of the product model as follows: where the pitching angle θ = θ P , rolling angle φ = θ R , and yawing angle ψ = θ Y When the aircraft attitude is determined, the aircraft location can be determined just with the constraints between the aircraft and the deck.When the rolling angle θ R is positive, the aircraft tilts to the right and the right main landing gear just touches the deck, otherwise the opposite.The initial condition of the carrier-based aircraft includes the carrier engaging speed V E , the sinking speed V V , and the aircraft rolling rate θ R These variables can be set by analyzing the model.In the eight landing variables stipulated by military standards and guides, V TD is not the variable at the moment of landing, so it cannot be considered in the virtual prototype.Besides, off-center condition is not taken into consideration.Therefore, both the ground critical velocity V TD and off-center arresting distance d are simplified by using their average ones.
The sinking speed of the carrier-based aircraft is large during landing, and the landing gear bears rather large impact load.Referring to the analysis of the land-based landing load of the aircraft [23], the three directional loads of the landing gear and barycenter overload are selected as landing load to carry out the study.The course load and vertical load of the landing gear act on the center of the tire axle, while the lateral load acts on the contact point between the tires and the ground.The bounce landing is not permitted, so the maximum values of the three directional load and barycenter overload are selected as landing loads for analysis.After the BDF solver finishing solving the landing process, the extreme values of the three directional load and barycenter overload under this condition are output by calling the analysis model.

EFAST Global Sensitivity Analysis Method
During the aircraft landing, there are many landing variables and landing loads.The local sensitivity analysis method is limited by the variable selection and model, so the analytical results are usually not comprehensive or accurate enough.The variable range of the global sensitivity analysis can be extended to the entire variable domain, which can provide rather integrated quantitative analysis results.EFAST method is used to study the global sensitivity of the landing variables to the landing loads.
EFAST is the quantitative global sensitivity analysis method which is developed by Saltelli et al. [18].Based on classical Fourier amplitude sensitivity test (FAST), Saltelli et al. combined Sobol variance decomposition thought to develop the method.This method maps the multidimensional variable to independent variable by appropriate searching curves.Then the search within the multidimensional variable  Assume that the landing analytical model is described as y = f X , where X = X x 1 , x 2 , … , x n , n = 6 is the combination of the landing variables.Before performing the global sensitivity analysis, all variables have to be mapped to the interval [0, 1] to constitute the n-dimensional unit super cube space.
Define an independent variable s and introduce the same mapping s → X, so where s −∞ < s < +∞ is the scalar variable, ω i ∈ Ζ = −∞, … , − 1, 0, 1, … , + ∞ , and its value changes to suitable value according to x i With the change of s, all x i change in the value space K n The oscillation frequency is ω i G i is the searching curve and recommended by Sati et al. [13] as where φ is the random value in 0, 2π Expand f X by Fourier series. where , and the amplitude of the Fourier series is From the properties of A j and B j , it is known that A −j = A j , B −j = −B j , and Λ −j = Λ j By calculating all frequencies, the total variance caused by each variable x i can be obtained, namely, where   International Journal of Aerospace Engineering The total variance of the model is According to Sobol variance decomposition thought, the model can be decomposed into functions of single variable and combined variable, namely, where V i = V E y|x i which represents the variance of x i , V i,j = V E y|x i , x j − V i − V j which represents the variance of x i contributed by x j , and V i,j,k ~V1,2,…,n are by analogy.
The first-order sensitivity coefficient S x i of the variable x i is defined as below.
The total effects of the variable x i include the single effect and interacting effects on the output of the model.The global sensitivity coefficient of the x i is defined as The difference between global sensitivity coefficient S T x i and first-order sensitivity coefficient S x i reflects the degree of interaction between x i and other variables.

Landing Variable Sensitivity and Extreme Value Analysis
When using EFAST method to analyze the global sensitivity, the input of probability distribution of the variable and the output of the model corresponding to the inputting samples are required [18].The landing variables given by standards [7] should be modified according to practical situation.The type of the landing variable distribution maintains constant, then the average and standard deviations of the landing variables need to be modified to satisfy the design requirements.
Table 1 shows the distribution of the six landing variables under the given design requirements.Due to the effects of the deck layout, the value range of the rolling angle and yawing angle are not bilateral symmetric.The input parameters of EFAST analysis method use the distribution of the landing variable and its value range obtained after adopting multivariate distribution presented in Section 2.2.The EFAST sampling is carried out by using (8), then 6438 samples are obtained.The samples in the generated interval [0, 1] are mapped back to the distribution of the landing variable, used to input the variables of the virtual prototype.For all landing conditions, the virtual prototype and its secondary development technology are used to analyze the landing loads automatically.According to the analytical results of the landing load, the first-order and global sensitivity coefficients of landing variables to landing loads are calculated.

The Sensitivity Analysis of Main Landing Gear Load.
When the carrier-based aircraft lands on the deck, main landing gear touches the deck first in most cases.The vertical impact load is the main load that acts on the main landing gear.The sensitivity analysis results of the vertical load of the main landing gear are shown in Figure 9.
The analysis results show that sinking speed and rolling angle are the main affecting variables of the vertical load of the main landing gear.The total contribution of the two variables is more than 82%, while the remaining landing variables only have little coupling effects on the vertical load.The sinking speed has larger first-order sensitivity coefficient in the analysis of the left main landing gear while the rolling angle has larger first-order sensitivity coefficient in the analysis of the right main landing gear.When the right main landing gear touches the deck first, the rolling angle is large and the landing load is relatively severe.When the left main landing gear touches the deck first, the rolling angle is small and the sinking speed is the main effect of the vertical load.
Figure 9 also shows the distribution of the cumulative frequency of the main landing gear's vertical load.The largest vertical load (853 kN) condition of the left main gear is V V = 4 76 m/s and θ R = −1 99 °The largest vertical load (1278 kN) condition of the right main gear is V V = 4 77 m/ s and θ R = 7 09 °The extreme value of the vertical load of the right landing gear is 1.5 times of the one of the left landing gear.
The sensitivity analysis results of the main landing gear's lateral load are shown in Figure 10.It is seen that the sinking speed, rolling angle, and yawing angle are the main affecting factors of the main landing gear's lateral load.The total contribution of the three variables is over 91% while remaining variables barely have direct effects on the lateral load.
From the cumulative frequency curve of the lateral load shown in Figure 10, it is seen that the extreme value of the lateral load of the left main landing gear is nearly the same as that of the right main landing gear.However, the average load of the right main landing gear is larger than that of the left main landing gear.Under the condition that the largest lateral load of the left main landing gear is 120 kN, the sinking speed V V = 4 62 m/s, rolling angle θ R = 7 36 °, and yawing angle θ Y = 0 84 °Under the condition that the largest lateral load of the right main landing gear is 122 kN, the sinking  11.The analysis results show that the sinking speed and rolling angle are the main affecting variables of the main landing gear's course load.The two variables' total contribution is 76% while the remaining landing variables have little coupling effects on the course load.Known from the cumulative frequency curve, the load distribution and extreme load of the left main landing gear is almost the same as that of the right main landing gear.Under the condition that the largest course load of the left main landing gear is 609 kN, the sinking speed V V = 4 59 m/s and rolling angle θ R = −2 04 °Under the condition that the largest course load of the right main landing gear is 721 kN, the sinking speed V V = 4 76 m/s and rolling angle θ R = 7 09 °The course loads of the left and right main landing gears have little difference.The extreme value of the course load also occurs under the condition of large sinking speed and large rolling angle.
Figures 9 and 11 show that the sensitivity analysis results of the course load are similar to that of the vertical load.The coefficients of sinking speed and rolling angle are nearly the same.Besides, the course loads of the right main landing gear are all greater than that of the left main landing gear.The course load has rather high correlation with the vertical load.

Nose Landing Gear Load and Barycenter Overload
Analysis.The sensitivity analysis results of the nose landing  gear's three directional loads are shown in Figure 12.It can be known that the sinking speed and rolling angle are the main affecting variables of the nose landing gear's course load and their first-order sensitivity coefficients are 0.59 and 0.14, respectively.The pitching angle is the secondary affecting variable, and its first-order sensitivity coefficient is 0.004.Other variables have a few coupling effects on the course load.
For the lateral load of the nose landing gear, the rolling angle is the main affecting variable and its first-order sensitivity coefficient is 0.6.The sinking speed, pitching angle, and yawing angle have a few effects on the lateral load.For the vertical load of the nose landing gear, the pitching angle, rolling angle, and sinking speed are the main affecting variables and their first-order sensitivity coefficients are 0.53, 0.14, and 0.07, respectively.Other variables do not have appreciated effects on the vertical load.
From the cumulative frequency curve of the nose landing gear's three directional loads, the sinking speed V V = 5 41 m/s, the pitching angle θ P = 0 62 °, and the rolling angle θ R = 3 93 °under the condition that the largest course load of the nose landing gear is 133 kN.Under the condition that the largest lateral load of the nose landing gear is 93 kN, the sinking speed V V = 4 77 m/s, the pitching angle θ P = 3 09 °, and the rolling angle θ R = 7 09 °, which is the large rolling angle condition.Under the condition that the largest vertical load of the nose landing gear is 196 kN, the sinking speed V V = 3 44 m/s, the pitching angle θ P = 6 23 °, and the rolling  The sensitivity analysis results of the barycenter overload of the airframe are shown in Figure 13.It is seen that sinking speed and rolling angle are the main affecting variables and their first-order sensitivity coefficients are 0.74 and 0.18, respectively.Under the condition that the largest barycenter overload is 5.11, the sinking speed V V = 4 77 m/s and the rolling angle θ R = 7 09 °, which is the condition of large rolling angle and relatively large sinking speed.The barycenter overload of the airframe mainly concentrates between 2.0 and 4.0.

Sensitivity Analysis of the Landing
Load of the Carrier-Based Aircraft.According to the sensitivity analysis results of the main landing gear, nose landing gear, and barycenter overload, the influence degree of the landing variables on landing loads can be classified as shown in Table 2.The sinking speed and pitching angle are the main affecting factors of the landing loads.The pitching angle and yawing angle have significant effects on some loads.The ground critical velocity and rolling rate have little effect on the landing loads, and they are influenced a lot by the change of the arresting cable load and the motion of the deck due to the restrictions of the simplified condition of the virtual prototype.

Extreme Value Analysis of the Landing Loads of the
Carrier-Based Aircraft.In the structure design stage, it is difficult to obtain the most severe load condition of the airframe during its service.When determining the most severe parameter combination, some conservative conditions are usually selected.For example, a certain or several extreme values of landing variables and the average of remaining landing variables are selected as extreme value conditions of the load for the researches.The conditions represent different landing conditions such as tail sinking, yawing landing, and landing with single landing gear.Some conditions will not happen during the service period, which leads to structural overdesign.
There are many landing samples selected in the sensitivity analysis of the landing load of the carrier-based aircraft.Approximately, the extreme load conditions of these samples can be used as extreme load conditions of all landing conditions to carry out the comparative analysis.From the computed results of the 6438 samples in EFAST sensitivity analysis method, the extreme value conditions of the three directional loads of the main and nose landing gear and barycenter overload are usually the condition that the sinking speed, rolling angle, or pitching angle are large and the landing variables do not pass the maximum of the value range.For further study in the rationality and security of the conservative condition design, typical condition design is selected for the comparative calculation.
Table 3 shows eight kinds of typical design conditions selected to compare with samples of EFAST sensitivity analysis in extreme loads of each direction.From the comparison, it is seen that the maximum gravity overload is slightly larger than that of the EFAST samples and the difference value is 23.5% approximately.
The extreme vertical and course loads of the landing gear of the typical design condition are nearly the same as the extreme values of EFAST samples.The maximum course load of the left and right main landing gear of the typical design condition is almost the same as the maximum load of the EFAST samples, and the percentage differences are about 6.2% and 16.0%, respectively.The maximum vertical load of the left and right main landing gear of the typical design condition is almost the same as the maximum load of the EFAST samples, and the percentage differences are about 2.5% and 12.0%, respectively.The maximum course and vertical loads of the nose landing gear of the typical design condition are quite different from the maximum loads of the EFAST samples, and the percentage differences are about 26.3% and 28.6%, respectively.
For the lateral load which is the secondary load, the maximum load of the typical design condition is quite different from that of the EFAST samples and the percentage difference is about 35.2%.From the overall analysis, the gravity overload and extreme load in each direction of the eight    International Journal of Aerospace Engineering kinds of typical design conditions are nearly the same as those of the EFAST samples, and the key loads of the typical design conditions are slightly larger.The extreme loads obtained by the conservative condition design are safe to be used in structure strength analysis and would not cause overdesign.

Conclusions
Based on the results reported, several conclusions can be drawn.Firstly, the multivariate distribution has a large influence on the distribution of the landing variables.After multivariate distribution processing, the distribution of the landing variables is still close to the normal distribution, and the mean values of the variables maintain constant while the standard deviation decreases.
Secondly, the results of global sensitivity analysis show that all the six landing variables have effects on landing loads.The sinking speed and rolling angle almost influence all kinds of landing loads significantly and are the main affecting factors.The pitching angle and yawing angle may only influence some landing variables a lot, and they are the secondary affecting factors.The results of extreme value analysis show that it is safe to adopt typical condition design in the structural strength analysis.

Figure 4 :
Figure 4: Comparisons of sample distributions before and after screening.

Figure 8 :
Figure 8: Flow chart of the virtual prototype multicondition automatic simulation.

6
Sensitivity factor of vertical load of the left main gear Frequency curve of vertical load of the left main gear Frequency curve of vertical load of the right main gear

Figure 9 :
Figure 9: Sensitivity and cumulative frequency curve of vertical load of the main landing gear.
Sensitivity factor of lateral load of the left main gear Frequency curve of lateral load of the left main gear Sensitivity factor of lateral load of the right main gear Frequency curve of lateral load of the right main gear

Figure 10 :
Figure 10: Sensitivity and cumulative frequency curve of lateral load of the main landing gear.

First 5 (
Sensitivity factor of course load of the left main gear Frequency curve of course load of the left main gear d) Frequency curve of course load of the right main gear

Figure 11 :
Figure 11: Sensitivity and cumulative frequency curve of course load of the main landing gear.

Figure 12 :
Figure 12: Sensitivity and cumulative frequency curve of three directional loads of the nose landing gear.
Frequency curve of barycenter overload

Figure 13 :
Figure 13: Sensitivity and cumulative frequency curve of the barycenter overload.

Table 1 :
Aircraft landing variables distribution.

Table 2 :
Sensibility of landing variables to landing loads.