Uncertainty Analysis of Aircraft Center of Gravity Deviation and Passenger Seat Allocation Optimization

: The traditional method of allocating passenger seats based on compartments does not effectively manage an aircraft’s center of gravity (CG), resulting in a notable divergence from the desired target CG (TCG). In this work, the Boeing B737-800 aircraft was employed as a case study, and row-based and compartment-based integer programming models for passenger allocation were examined and constructed with the aim of addressing the current situation. The accuracy of CG control was evaluated by comparing the row-based and compartment-based allocation techniques, taking into account different bodyweights and numbers of passengers. The key contribution of this research is to broaden the range of the mobilizable set for the aviation weight and balance (AWB) model, resulting in a significant reduction in the range of deviations in the center of gravity outcomes by a factor of around 6 to 16. The effectiveness of the row-based allocation approach and the impact of passenger weight randomness on the deviation of an airplane’s CG were also investigated in this study. The Monte Carlo method was utilized to quantify the uncertainty associated with passenger weight, resulting in the generation of the posterior distribution of the aircraft’s center of gravity (CG) deviation. The outcome of the row-based model test is the determination of the range of passenger numbers that can be effectively allocated under different TCG conditions.


Introduction
The process of aviation weight and balance planning involves the development of a strategic plan to allocate passengers, carriages, baggage, and mail to various cabin compartments.The aim of this plan is to maintain the equilibrium of the aircraft throughout operation, taking into consideration factors such as load capacity and the position of the center of gravity.The primary goal of this system is to maintain the aircraft's center of gravity (CG) and overall weight within a safe range, as this is a crucial factor in ensuring the safety of flight operations.
Aircraft exhibit distinct characteristics compared to land or sea transportation vehicles, as depicted in Figure 1.Land and water transport vehicles have the ability to bear weight by means of wheels on the ground or hulls on water surfaces, thereby providing multi-point support and ensuring high stability.Their load distribution balancing criteria are often less stringent, with differential weight constraints imposed between the left and right sides as well as the front and rear of the body.Nevertheless, during flight, the aerodynamic force exerted by the weight is primarily concentrated on the wings, resulting in a singular point of force.Additionally, the dynamic nature of the atmospheric environment necessitates the maintenance of the aircraft's operational stability.Hence, in contrast to land and marine transportation, the balance requirements of air transportation are exceedingly stringent, thus limiting the CG within a specific range and mandating that it is close to a predetermined target CG (TCG).
athematics 2024, 12, x FOR PEER REVIEW 2 of 1 exceedingly stringent, thus limiting the CG within a specific range and mandating that i is close to a predetermined target CG (TCG).The aircraft CG position is the basis of the aircraft's balance during operation, and it analysis and evaluation have attracted the attention of many scholars.Idan, Iosilevskii and Nazarov (2003) used artificial neural networks (NNs) to construct an aircraft weigh and CG position [1].Zhang, Yang, and Shen (2012) proposed a hybrid CG estimation and sliding mode adaptive control approach and introduced an online CG estimation method based on adaptive weight data fusion to guarantee smaller dynamic inversion errors [2] Yang et al. ( 2016) proposed a Gaussian process regression model for the estimation of the CG location of fixed-wing aircraft [3].Their model is a data-based approach explicitly tackling uncertainties caused by the quality and quantity of the data.In the work of Juma and Abdullah (2018), a baggage load planning system for a Boeing 737-800 aircraft is con sidered to obtain a good CG position for improved fuel consumption [4].Chaves, Silves tre, and Gamboa (2018) presented an approach to onboard in-flight weight and balance estimation systems for the identification of weight and CG position using cruise angle o attack, Mach number, and elevator deflection values [5].Zhao et al. (2021) constructed a mixed integer programming model to maximize the total payload and minimize the CG deviation and proposed a new set of CG envelope constraints [6].
The CG of the AWB plan is not exactly equal to the airlines' TCG.Its deviation i uncertain because of the randomness of individual passenger weight and seat allocation Normally, airline AWB planning assigns passengers based on aircraft compartments, no seats in operation.This causes little diversity in allocation plans, and passengers' seat al locations are uncertain.Moreover, not many studies have been conducted on cabin as signments, as their focus has been on cabin seat comfort and passenger seat selection pref erences as well as other relevant factors.Liu, YX et al. ( 2023) improved cabin seating dis tribution using the Wells-Riley model to reduce the probability of passengers contracting respiratory illnesses and to improve the economic efficiency of airlines [21].From the per Lift Supportiveness Supportiveness tiveness The aircraft CG position is the basis of the aircraft's balance during operation, and its analysis and evaluation have attracted the attention of many scholars.Idan, Iosilevskii, and Nazarov (2003) used artificial neural networks (NNs) to construct an aircraft weight and CG position [1].Zhang, Yang, and Shen (2012) proposed a hybrid CG estimation and sliding mode adaptive control approach and introduced an online CG estimation method based on adaptive weight data fusion to guarantee smaller dynamic inversion errors [2].Yang et al. (2016) proposed a Gaussian process regression model for the estimation of the CG location of fixed-wing aircraft [3].Their model is a data-based approach explicitly tackling uncertainties caused by the quality and quantity of the data.In the work of Jumat and Abdullah (2018), a baggage load planning system for a Boeing 737-800 aircraft is considered to obtain a good CG position for improved fuel consumption [4].Chaves, Silvestre, and Gamboa (2018) presented an approach to onboard in-flight weight and balance estimation systems for the identification of weight and CG position using cruise angle of attack, Mach number, and elevator deflection values [5].Zhao et al. (2021) constructed a mixed integer programming model to maximize the total payload and minimize the CG deviation and proposed a new set of CG envelope constraints [6].
The position of the aircraft's center of gravity (CG) during operation is primarily influenced by factors such as the aircraft's empty weight CG; fuel quantity; and the distribution of passengers, cargo, baggage, and other relevant elements.The fuel quantity for an AWB plan is typically estimated based on the aircraft's empty weight CG for a particular flight.The primary function of stowage is the precise placement and distribution of passengers, freight, luggage, and other related items.In terms of air cargo transportation, several mathematical models and optimization methods have been developed to address AWB (Larsen & Mikkelsen, 1980 [7]; Brosh, 1981 [8]; Amiouny et al., 1992 [9]; Wodziak & Fadel, 1994 [10]; Mathur, 1998 [11]; Heidelberg et al., 1998 [12]; Dahmani & Krichen, 2013 [13]).However, as Brandt and Nickel (2019) have mentioned, the problems with real-world scales and practical constraints have yet to be solved [14].Regarding air passenger transportation, the literature on AWB planning mainly concernspassenger seat allocation, with several real-world operation requests (Schultz and Soolaki 2021 [15]; Ren, Pan, and Jiang 2022 [16]; Notomista et al. 2016 [17]; Birolini et al. 2022 [18]; Wong et al. 2009 [19]; Ma et al. 2023 [20]).
The CG of the AWB plan is not exactly equal to the airlines' TCG.Its deviation is uncertain because of the randomness of individual passenger weight and seat allocation.Normally, airline AWB planning assigns passengers based on aircraft compartments, not seats in operation.This causes little diversity in allocation plans, and passengers' seat allocations are uncertain.Moreover, not many studies have been conducted on cabin assignments, as their focus has been on cabin seat comfort and passenger seat selection preferences as well as other relevant factors.Liu, YX et al. (2023) improved cabin seating distribution using the Wells-Riley model to reduce the probability of passengers contracting respiratory illnesses and to improve the economic efficiency of airlines [21].From the perspective of an airline company, J Qiu et al. (2017) collected and analyzed passengers' content and operating habits regarding cabin multimedia devices, built a model of passengers' personalities and preferences, predicted passengers' likely positive opinions of each other, and intelligently assigned them to adjacent seats [22].Castro, J et al. (2021) proposed a seat allocation model considering the future sales of the airline [23].Milne, R. J. et al. (2024) used a greedy algorithm to assign seats to passengers in order to separate infectious passengers from susceptible ones, reducing the rate of disease transmission [24].However, all of these methods can cause the CG of the AWB plan to vary considerably from the TCG.Moreover, passengers vary greatly in terms of individual weight and cabin carry-on baggage.For example, in 2023, a large number of South Korean passengers secretly took heavy grain on airplanes from China back to their country, where rice is quite expensive.It was notably reported that the weight of the cabin carry-on baggage of these "grain customers" was 5-8 times that of an ordinary passenger.However, this factor was ignored in AWB planning and caused other CG errors.Additional ballasts are adopted by most airlines to reduce CG errors and limit CG deviations within reasonable bounds.This makes the aircraft heavier, increasing fuel consumption and carbon emissions.There are many studies in the research literature on uncertainty, such as Birge and Louveaux 2011 [25]; Hannan et al. 2020 [26]; Crainic, Djeumou Fomeni, and Rei 2021 [27]; and Crainic et al. 2024 [28].However, to the best of our knowledge, there is little research on uncertainty in AWB planning.
The current body of research on passenger loading mostly concentrates on the operational economics and implementation of loading systems, with loading being considered as an additional instrument.Nevertheless, there is a notable lack of research on the control intensity of the CG considering the results of passenger allocation and the uncertainty associated with passenger weight.In addressing this research gap, the focus of this study is on maximizing the model's utility by investigating how the allocation of passengers and cargo can improve control of the aircraft's CG, as well as setting an appropriate TCG.
In order to mitigate the CG deviation and enhance loading accuracy, a row-based linear integer programming model was devised to optimize the allocation of passenger seats.In contrast to the compartment-based model, the proposed approach has the potential to effectively reduce CG deviation and achieve optimal control over the aircraft's CG configuration.Furthermore, by employing the Monte Carlo methodology, the posterior distribution of the aircraft's CG was derived based on the varying weight distributions of passengers.The influence of passenger weight uncertainty on airplane CG was quantified by calculating confidence intervals at different levels of the CG deviation distribution.The optimal placement of the TCG was ascertained through a comparative analysis of the confidence intervals for a range of appropriate deviation distributions.This improves the usefulness of the scheduling model.
The sections of this article are organized as follows: Section 1 provides the introduction and introduces the background and topic; Section 2 describes the compartment-based and row-based mathematical models; Section 3 compares the results of the two AWB models; Section 4 describes the analysis of the effect of passenger weight uncertainty on the center of gravity of an aircraft using the row-based AWB model; and Section 5 summarizes the conclusions.

Description and Model Formulation
Our methodology involves loading the aircraft with a target center of gravity (PTCG) position.The goal is to obtain a center of gravity (PCG) that closely matches this target value.The center of gravity of an aircraft is calculated by dividing the total moment on the aircraft by the total weight, as described in Equation (1): where z is the index of passenger compartments; B z is the arm of the zth compartment; h is the index of cargo holds; B h is the arm of the hth lower cargo hold; B OEW is the arm of the aircraft's empty operational weight (OEW); B TOF is the arm of the aircraft's takeoff fuel (TOF); p is the index of passengers; w P is the average weight of a passenger; W OEW is the aircraft's empty operational weight; W TOF is the takeoff fuel weight; x z is the number of passengers allocated to the zth compartment; and w h is the weight allocated to the hth lower cargo hold.
In practical applications, the relative position (P CG ) is a more convenient and widely used representation compared to the absolute position (B CG ).It is defined as the ratio of the wing reference distance measured from the leading edge of the mean aerodynamic chord (MAC) to the CG position to the MAC chord length (as illustrated in Figure 2. Therefore, the expression for the P CG can be represented as Equation ( 2).Its unit is %MAC.
where L MAC is the length of the MAC and L EMAC is the distance to the leading edge of the MAC.where z is the index of passenger compartments; Bz is the arm of the zth compartment; h is the index of cargo holds; Bh is the arm of the hth lower cargo hold; BOEW is the arm of the aircraft's empty operational weight (OEW); BTOF is the arm of the aircraft's takeoff fuel (TOF); p is the index of passengers; wP is the average weight of a passenger; WOEW is the aircraft's empty operational weight; WTOF is the takeoff fuel weight; xz is the number of passengers allocated to the zth compartment; and wh is the weight allocated to the hth lower cargo hold.
In practical applications, the relative position (PCG) is a more convenient and widely used representation compared to the absolute position (BCG).It is defined as the ratio of the wing reference distance measured from the leading edge of the mean aerodynamic chord (MAC) to the CG position to the MAC chord length (as illustrated in Figure 2. Therefore, the expression for the PCG can be represented as Equation ( 2).Its unit is %MAC.
where LMAC is the length of the MAC and LEMAC is the distance to the leading edge of the MAC.

Figure 2.
The relationship between the absolute position (BCG) and relative position (PCG).

The Compartment-Based AWB Model
When allocating passengers to seats within an airplane cabin, it is customary to partition the cabin into various compartments based on seat classes, denoted by codes such as OA, OB, OC, and so on (as depicted in Figure 3).The number of passengers in each compartment is determined by calculating the weight balance.In order to determine the CG, the mean arm length inside each compartment is utilized, without taking into account the differences in arm values among seats in various rows within each compartment.The computation of the CG is thereby simplified.

The Compartment-Based AWB Model
When allocating passengers to seats within an airplane cabin, it is customary to partition the cabin into various compartments based on seat classes, denoted by codes such as OA, OB, OC, and so on (as depicted in Figure 3).The number of passengers in each compartment is determined by calculating the weight balance.In order to determine the CG, the mean arm length inside each compartment is utilized, without taking into account the differences in arm values among seats in various rows within each compartment.The computation of the CG is thereby simplified.where z is the index of passenger compartments; Bz is the arm of the zth compartment; h is the index of cargo holds; Bh is the arm of the hth lower cargo hold; BOEW is the arm of the aircraft's empty operational weight (OEW); BTOF is the arm of the aircraft's takeoff fuel (TOF); p is the index of passengers; wP is the average weight of a passenger; WOEW is the aircraft's empty operational weight; WTOF is the takeoff fuel weight; xz is the number of passengers allocated to the zth compartment; and wh is the weight allocated to the hth lower cargo hold.
In practical applications, the relative position (PCG) is a more convenient and widely used representation compared to the absolute position (BCG).It is defined as the ratio of the wing reference distance measured from the leading edge of the mean aerodynamic chord (MAC) to the CG position to the MAC chord length (as illustrated in Figure 2. Therefore, the expression for the PCG can be represented as Equation ( 2).Its unit is %MAC.
where LMAC is the length of the MAC and LEMAC is the distance to the leading edge of the MAC.

Figure 2.
The relationship between the absolute position (BCG) and relative position (PCG).

The Compartment-Based AWB Model
When allocating passengers to seats within an airplane cabin, it is customary to partition the cabin into various compartments based on seat classes, denoted by codes such as OA, OB, OC, and so on (as depicted in Figure 3).The number of passengers in each compartment is determined by calculating the weight balance.In order to determine the CG, the mean arm length inside each compartment is utilized, without taking into account the differences in arm values among seats in various rows within each compartment.The computation of the CG is thereby simplified. s.t.
x z ≤ N P (4) P CG ≥ P FCG (9) The aim of objective ( 1) is to minimize the deviation between the loaded CG and TCG.The aircraft CG is typically expressed in two ways: absolute position (BCG) and relative position (PCG).
Constraint (2) is the total passenger limit; Constraint (3) is the maximum number limit for each compartment; Constraint (4) is the total weight limit for cargo, baggage, and mail; Constraint (5) is the maximum load limit for each cargo hold; Constraint (6) is the total payload limit for the aircraft; and Constraints (7) and ( 8) are the pitch balance constraints for the aircraft, indicating that the CG of the aircraft needs to be within the forward and aft limits.
A somewhat aftward center of gravity (CG) typically results in decreased aerodynamic drag during flight, thereby conserving fuel.In order to adhere to national regulations for "carbon peaking" and "carbon neutrality", airlines frequently opt for a somewhat aftward TCG with the aim of achieving cost savings.During AWB planning, the distribution of passengers and cargo is modified to minimize the difference between the loaded CG and the TCG.Nevertheless, potential variations arise due to the intrinsic limitations associated with the allocation of passengers per compartment.In addition, the variability in passenger and cargo weights, offsets in cargo center of gravity, and differences in fuel weight and density collectively lead to increased loading uncertainty, as depicted in Figure 4.
A somewhat aftward center of gravity (CG) typically results in decreased aerod drag during flight, thereby conserving fuel.In order to adhere to national regulati "carbon peaking" and "carbon neutrality," airlines frequently opt for a somewhat a TCG with the aim of achieving cost savings.During AWB planning, the distribution sengers and cargo is modified to minimize the difference between the loaded CG a TCG.Nevertheless, potential variations arise due to the intrinsic limitations associate the allocation of passengers per compartment.In addition, the variability in passeng cargo weights, offsets in cargo center of gravity, and differences in fuel weight and d collectively lead to increased loading uncertainty, as depicted in Figure 4. .The objective of AWB planning |ΔCG|.When allocating passengers based on compartments, it is challenging to a precise control over the CG.This is because there are limitations on the available w adjustable compartments, and number of passengers.As a result, |ΔCG| is greate As an illustration, if there is just one traveler on the flight and the TCG is 17% assigning the passenger to separate compartments will produce distinct CGs.F B737-800, the values would be 16.9244% in the OA compartment, 17.2177% in the O 17.8844% in the OC, and |ΔCG| is at least 0.0756% MAC.The third row gives the p ger a CG position of 17.0033% MAC and reduces |ΔCG| to 0.0033% MAC if the pas is allocated according to seats.This approach has a noticeably higher degree of acc A smaller value of |CG| indicates a narrower disparity between the loaded a get CG.Airlines have the option to select a more forward-looking TCG to reduce fu sumption, carbon emissions, and expenses as well as a more efficient AWB design t inate the requirement for additional fuel and ballasts in order to enhance strategi ning.Furthermore, considering effective administration of the aircraft's CG can the likelihood of mishaps caused by an imprecise CG, hence improving flight safet sequently, it becomes essential to partition the airplane cabin into each row of se the assigned passengers.

The Row-Based AWB Model
The AWB model for air transportation based on seating rows is described as f The quantity of cargo assigned to each cargo hold and the number of passengers as to each row of seats in the cabin serve as the decision variables.In order to lower and enhance the accuracy of the AWB plan, a linear mixed integer programming was developed to assign passengers by rows while accounting for restrictions s cargo and seat capacities, weight, and balance limits.
Let xr be the integer decision variable representing the number of passeng signed to the rth row, where r = 1, 2, ..., NR, and NR is the number of rows in the cab Let |∆CG| := |P CG − P TCG | and ∆CG := P CG − P TCG .The objective of AWB planning is min |∆CG|.When allocating passengers based on compartments, it is challenging to achieve precise control over the CG.This is because there are limitations on the available weight, adjustable compartments, and number of passengers.As a result, |∆CG| is greater.
As an illustration, if there is just one traveler on the flight and the TCG is 17% MAC, assigning the passenger to separate compartments will produce distinct CGs.For the B737-800, the values would be 16.9244% in the OA compartment, 17.2177% in the OB, and 17.8844% in the OC, and |∆CG| is at least 0.0756% MAC.The third row gives the passenger a CG position of 17.0033% MAC and reduces |∆CG| to 0.0033% MAC if the passenger is allocated according to seats.This approach has a noticeably higher degree of accuracy.
A smaller value of |CG| indicates a narrower disparity between the loaded and target CG.Airlines have the option to select a more forward-looking TCG to reduce fuel consumption, carbon emissions, and expenses as well as a more efficient AWB design to eliminate the requirement for additional fuel and ballasts in order to enhance strategic planning.Furthermore, considering effective administration of the aircraft's CG can reduce the likelihood of mishaps caused by an imprecise CG, hence improving flight safety.Consequently, it becomes essential to partition the airplane cabin into each row of seats for the assigned passengers.

The Row-Based AWB Model
The AWB model for air transportation based on seating rows is described as follows: The quantity of cargo assigned to each cargo hold and the number of passengers assigned to each row of seats in the cabin serve as the decision variables.In order to lower |∆CG| and enhance the accuracy of the AWB plan, a linear mixed integer programming model was developed to assign passengers by rows while accounting for restrictions such as cargo and seat capacities, weight, and balance limits.
Let x r be the integer decision variable representing the number of passengers assigned to the rth row, where r = 1, 2, . .., N R , and N R is the number of rows in the cabin.Let the decision variable w h represent the weight assigned to the hth lower cargo hold, h = 1, 2, . .., NH.The following is the row-based AWB model.
x r ≤ N P (12) x r ≤ n r z = 1, 2, . . ., N R ( 13) P CG ≥ P FCG (17) where n r represents the number of seats in the rth row of the cabin.The remaining terms are consistent with the compartment-based AWB model.

Comparison of the CG Deviation
We used the commercial solver Gourbi to solve our model.The efficiency of the passenger allocation optimization procedure was evaluated by calculating and analyzing the models of the row-based and compartment-based AWB models.We also provide the results of the traditional manual loading method (using manifests for loading).The test results with the TCG at 20% MAC and passenger counts ranging from 30 to 125 are presented in Table 1 and Figure 5.The row-based allocation model has the lowest |CG| values and the smallest gap between the real center of gravity and the ideal center of gravity, as Table 1 and Figure 5 illustrate.The traditional manual loading method has the highest |CG| values and the largest gap between the real center of gravity and the ideal center of gravity, as shown in Table 1 and Figure 5.The results from compartment-based allocation range from 0.001% MAC to 0.03% MAC, manual loading results range from 0.005% MAC to 0.3% MAC, and row-based allocation yields the best results from 0.00006% MAC to 0.005% MAC.The row-based allocation model has the lowest |CG| values and the smallest gap between the real center of gravity and the ideal center of gravity, as Table 1 and Figure 5 illustrate.The traditional manual loading method has the highest |CG| values and the largest gap between the real center of gravity and the ideal center of gravity, as shown in Table 1 and Figure 5.The results from compartment-based allocation range from 0.001% MAC to 0.03% MAC, manual loading results range from 0.005% MAC to 0.3% MAC, and row-based allocation yields the best results from 0.00006% MAC to 0.005% MAC.
Table 2 demonstrates how much faster the compartment-based model produces results than the row-based model.The speed of the row-based model is due to its search space being over ten times larger than that of the compartment-based model.

Comparison of Results with Different TCG
We compared the ∆CGs of compartment-based and row-based models using different TCGs (16%, 19%, 22%, 25%, and 28% MAC) as shown in Figure 6.The allocation methodology utilizing rows results in lower values of ∆CGs compared to the compartment-based method, as depicted in Figure 5.This is due to the fact that the row-based model has 29 rows that can be assigned and about 10 times more sets that can be mobilized than the compartment-based model, which only has three compartments that can be assigned and fewer sets that can be mobilized.The compartment-based model's failure to handle the CG when the number of passengers increases leads to periodic fluctuations in the ∆CGs.The model is limited to cyclically distributing people between several compartments.Furthermore, as the TCG increases in size, the fluctuations become increasingly prominent and exhibit a higher degree of periodicity.It can also be observed from Figure 5

The Effect of Variances in Passenger Weight on CG
Most airlines in China use 75 kg as the normal passenger weight for allocation.However, planners rely on calculations that may not adequately represent the current population demographics, and the exact weight of the passengers is unknown (Melis et al. 2019 [29]).The use of standard passenger weights results in an overestimation of the aircraft's performance characteristics, particularly its change in center of gravity (ΔCG) and fuel cost.
We express the passenger weight distribution as wr~N (75, 15).The passenger weight distribution is represented as a normal distribution with a mean of 75 and a standard deviation of 15.According to the "Aircraft Weight and Balance Control Regulations (CAAC 2019 [30])" and the "2022 China Resident Height and Weight Health Data Report" (as shown in Figure 7), the average weight of passengers and their cabin luggage is 75 kg.The weight fluctuation is 15 kg, considering differences in weight between genders and other aspects such as whether passengers are carrying bags.

The Effect of Variances in Passenger Weight on CG
Most airlines in China use 75 kg as the normal passenger weight for allocation.However, planners rely on calculations that may not adequately represent the current population demographics, and the exact weight of the passengers is unknown (Melis et al. 2019 [29]).The use of standard passenger weights results in an overestimation of the aircraft's performance characteristics, particularly its change in center of gravity (∆CG) and fuel cost.
We express the passenger weight distribution as w r ~N (75, 15).The passenger weight distribution is represented as a normal distribution with a mean of 75 and a standard deviation of 15.According to the "Aircraft Weight and Balance Control Regulations (CAAC 2019 [30])" and the "2022 China Resident Height and Weight Health Data Report" (as shown in Figure 7), the average weight of passengers and their cabin luggage is 75 kg.The weight fluctuation is 15 kg, considering differences in weight between genders and other aspects such as whether passengers are carrying bags.
distribution is represented as a normal distribution with a mean of 75 and a standard deviation of 15.According to the "Aircraft Weight and Balance Control Regulations (CAAC 2019 [30])" and the "2022 China Resident Height and Weight Health Data Report" (as shown in Figure 7), the average weight of passengers and their cabin luggage is 75 kg.The weight fluctuation is 15 kg, considering differences in weight between genders and other aspects such as whether passengers are carrying bags.

The Effect on the Mean
We utilized the Monte Carlo method to evaluate the uncertainty linked to passenger weight.The procedures were as follows: Firstly, the row-based model was utilized to ascertain the number of passengers per seat in the cabin.Secondly, a normal distribution was employed to generate random values (w r ) for Nr passengers.Subsequently, the previously standardized weight of passengers (75 kg) was replaced.Lastly, the computation of the CG was performed for each substitution.By iteratively substituting the passenger weight, we obtained a posterior distribution of the CG that takes into account the passenger's uncertainty.Each of the numbers of passengers (120, 140, and 160) was tested 1000 times with the aim of determining the distributions of CG deviation.The TCG of 24% MAC was used for these tests, as shown in Figure 8.

The Effect on the Mean
We utilized the Monte Carlo method to evaluate the uncertainty linked to passen weight.The procedures were as follows: Firstly, the row-based model was utilized to certain the number of passengers per seat in the cabin.Secondly, a normal distribut was employed to generate random values (wr) for Nr passengers.Subsequently, the p viously standardized weight of passengers (75 kg) was replaced.Lastly, the computat of the CG was performed for each substitution.By iteratively substituting the passeng weight, we obtained a posterior distribution of the CG that takes into account the pass ger's uncertainty.Each of the numbers of passengers (120, 140, and 160) was tested 10 times with the aim of determining the distributions of CG deviation.The TCG of 2 MAC was used for these tests, as shown in Figure 8. Figure 8 demonstrates that the CG deviation follows a normal distribution as a res of the uncertainty in passenger weight.There is a notable change in the average value the absolute difference between the CG when the population size approaches 160.T issue arises when the TCG position is set excessively towards the rear (or front), result in the row-based allocation model saturation of the seats in either the back or front ro due to an excess of passengers.In this scenario, the model can only allocate passengers a sequential manner, either starting from the back and moving towards the front or sta ing from the front and moving towards the back.This results in a significant magnitu of |ΔCG|, leading to a loss of control over the CG.Furthermore, there exists a correlat Figure 8 demonstrates that the CG deviation follows a normal distribution as a result of the uncertainty in passenger weight.There is a notable change in the average value of the absolute difference between the CG when the population size approaches 160.This issue arises when the TCG position is set excessively towards the rear (or front), resulting in the row-based allocation model saturation of the seats in either the back or front rows due to an excess of passengers.In this scenario, the model can only allocate passengers in a sequential manner, either starting from the back and moving towards the front or starting from the front and moving towards the back.This results in a significant magnitude of |∆CG|, leading to a loss of control over the CG.Furthermore, there exists a correlation between the weight of passengers and the quantity of individuals that can be accommodated.In order to incorporate the impact of fluctuations in passenger weight into the CG deviation, it is necessary to consider the effective range of control over the number of passengers.
The bisection approach was used to determine the effective control range of the number of passengers for the model with a TCG of 24% MAC.Two sets of numbers of passengers, specifically (12, 13, 14) and (158, 159, 160), were examined, and the outcome is depicted in Figures 9 and 10.The average value of the distribution of ∆CG deviates from 24% MAC (the average of |∆CG| > 1% MAC) when the number of passengers exceeds 159 or falls below 13.Therefore, the effective control range of the number of passengers for the model with a TCG of 24% MAC is determined to be between 13 and 159.That is, assuming a TCG of 24% MAC, the row-based allocation model can optimally distribute seats when the passenger count is between 14 and 158.The TCG is a significant factor in determining the magnitude of |ΔCG|.Thus, we can determine the effective ranges of numbers of passengers by repeating the preceding operation and considering various TCGs and numbers of passengers, as depicted in Figure 11.The TCG is a significant factor in determining the magnitude of |ΔCG|.Thus, we can determine the effective ranges of numbers of passengers by repeating the preceding operation and considering various TCGs and numbers of passengers, as depicted in Figure 11.The TCG is a significant factor in determining the magnitude of |∆CG|.Thus, we can determine the effective ranges of numbers of passengers by repeating the preceding operation and considering various TCGs and numbers of passengers, as depicted in Figure 11.The TCG is a significant factor in determining the magnitude of |ΔCG determine the effective ranges of numbers of passengers by repeating the p tion and considering various TCGs and numbers of passengers, as depicte

The Effect on the Standard Deviation
We analyzed the standard deviation of the distributions of |∆CG|, which are displayed in Figure 12.This analysis was conducted using a TCG of 22% MAC and a number of passengers ranging from 20 to 120.Although the average distribution remains continuous near the TCG, the variability in the absolute difference in the CG (|∆CG|) is much different.

The Effect on the Standard Deviation
We analyzed the standard deviation of the distributions of |ΔCG|, which are d played in Figure 12.This analysis was conducted using a TCG of 22% MAC and a num of passengers ranging from 20 to 120.Although the average distribution remains cont uous near the TCG, the variability in the absolute difference in the CG (|ΔCG|) is mu different.
Further experiments were carried out in order to monitor the fluctuations in standard deviation values as the number of passengers changed (Figure 13).Research shown that there is a positive correlation between the standard deviation of ΔCG dis butions and the number of passengers.However, this correlation weakens as the num of passengers increases.The correlation between the standard deviation and the num of passengers for other TCGs is comparable to 22% MAC.Further experiments were carried out in order to monitor the fluctuations in the standard deviation values as the number of passengers changed (Figure 13).Research has shown that there is a positive correlation between the standard deviation of ∆CG distributions and the number of passengers.However, this correlation weakens as the number of passengers increases.The correlation between the standard deviation and the number of passengers for other TCGs is comparable to 22% MAC.

The Confidence Interval of the Mean of ΔCG
All other ΔCGs are uniformly regarded as the same distribution when passengers is within the model's effective range.According to the principl entropy, this distribution must encompass all deviations in the center of g senger counts that are within the effective range of the model (Figure 14).distribution was chosen to be a normal distribution with a mean (μ) of −0.2 standard deviation of 1.1% MAC.Subsequently, the deviations of the ce were calculated for confidence levels of 90%, 95%, and 99%.The findings ar Table 3.

The Confidence Interval of the Mean of ∆CG
All other ∆CGs are uniformly regarded as the same distribution when the number of passengers is within the model's effective range.According to the principle of maximum entropy, this distribution must encompass all deviations in the center of gravity for passenger counts that are within the effective range of the model (Figure 14).The maximum distribution was chosen to be a normal distribution with a mean (µ) of −0.2% MAC and a standard deviation of 1.1% MAC.Subsequently, the deviations of the center of gravity were calculated for confidence levels of 90%, 95%, and 99%.The findings are displayed in Table 3.   From Table 3, it can be seen that the deviation of the center of gravity within the valid range of the model does not exceed 3.5% MAC.so according to this result the center of gravity envelope is shrunk accordingly.From Table 3, it can be seen that the deviation of the center of gravity within the valid range of the model does not exceed 3.5% MAC.so according to this result the center of gravity envelope is shrunk accordingly.

The Effect of the Number of Passengers Outside Effective Ranges
While there is a small probability that the number of passengers will fall outside of effective ranges during normal flight operations, the means of the ∆CG distribution show a noticeable change.We calculated the distribution of ∆CG using TCGs ranging from 17% to 24% MAC and passenger counts that were outside of their effective limits.The ∆CG shifts to the left for TCGs ranging from 21 to 24% MAC and to the right for TCGs ranging from 17 to 20% MAC as the number of passengers increases, as depicted in Figure 15.In addition, the means exhibit an arithmetic progression with TCGs for the same number of passengers.More precisely, with every 1% MAC increase in TCGs, the mean varies by approximately 1.01% MAC.When the TCG is constant, the mean shifts by approximately 0.3% of the mean for every additional passenger.When the number of passengers is below the lower limit of the model's effective number of assignments, the average value of ∆CG approaches zero as the number of people increases.This is because the increase in the number of people improves the adaptability of the model's assignments.When the number of passengers achieves the capacity of the cabin, there is variability in the standard deviation within the range of 0.55 to 0.6.When the number of passengers is small, the relationship between the standard deviation and the number of people is more pronounced (Figure 16).Tables 4 and 5 display the confidence intervals that we computed outside the upper and lower bounds of the model's effective allocation range.When the number of passengers achieves the capacity of the cabin, there is variability in the standard deviation within the range of 0.55 to 0.6.When the number of passengers is small, the relationship between the standard deviation and the number of people is more pronounced (Figure 16).Tables 4 and 5  When the number of passengers achieves the capacity of the cabin, there is variability in the standard deviation within the range of 0.55 to 0.6.When the number of passengers is small, the relationship between the standard deviation and the number of people is more pronounced (Figure 16).Tables 4 and 5

The Effect on the Mean
The following findings can be deduced by examining the correlation between the absolute value of the change in center of gravity (|∆CG|) and the uncertainty in passenger weight: When the number of passengers is within the range that the model can effectively regulate, the change in center of gravity (|∆CG|) is insignificant.However, as the number of passengers surpasses this range, the distribution of |∆CG| changes, and the deviation increases.Given that a majority of cases have an occupancy rate of above 80% or 90%, situations that fall below the lower limit of model effectiveness are uncommon.Hence, when addressing the issue of passenger weight uncertainty, the main focus should be on improving the maximum number of passengers that the allocation model can efficiently control.Figure 14 clearly demonstrates that a TCG of 20-21% MAC decreases both the average value of |∆CG| and the range of numbers of passengers that have no effect.Hence, by establishing the TCG to range from 20 to 21% MAC, the efficiency of the row-based allocation model can be optimized.

Discussion
Firstly, the row-based model provides approximately 6 to 16 times greater control of the center of gravity than the compartment-based model and 60 to 80 times greater control than traditional manual loading methods.Therefore, the row-based AWB model can significantly improve the effectiveness of civil passenger transportation allocation.Secondly, the recommended TCG for the row-based model achieves an effective headcount range of 95.29% to 98.23% of the actual range of passengers.This effective range is approximately 4-10 percentage points higher than the least suggested TCG placements (17% MAC effective range at 91.17% and 24% MAC effective range at 86.47%).Moreover, the suggested TCG position has about 2-4 times the center of gravity control of the least recommended TCG position when the number of passengers exceeds the upper boundary of the effective range of the row-based AWB model.However, the level of control that the recommended TCG position has over the least recommended TCG position decreases when the number of passengers falls below the lower boundary of the effective range of the row-based ABW model.Nevertheless, according to the regulations outlined in the relevant legislation, it is highly uncommon for an aircraft with a capacity of 170 individuals to transport so few people.To summarize the above two points, the recommended TCG in the row-based model not only caters to a wider range of numbers of passengers on a flight but also exhibits more effective control over the center of gravity, with general applicability beyond this range.

Conclusions
In this study, a linear integer programming methodology was utilized to develop a precise airplane loading and balancing strategy.Passengers in the cabin are allocated in seat rows, while cargo is distributed in the cargo holds.The row-based allocation approach surpasses the conventional compartment-based strategy in terms of achieving precise control over the aircraft's TCG and having a wider effective control range for the number of passengers.This method enhances the control of the center of gravity (CG) and flexibility of distribution.The row-based allocation model was used in this study to examine the relationship between the ∆CG and the uncertainty in passenger weight.The optimization of the TCG configuration is based on the established upper limits for the number of passengers.The subsequent course of action should involve undertaking further investigation into loading optimization and robust optimization in the presence of uncertainties related to passenger weight, baggage, cargo, mail, and other relevant factors.

Figure 1 .
Figure 1.Force characteristics of different means of transport.

Figure 1 .
Figure 1.Force characteristics of different means of transport.

Figure 2 .
Figure 2. The relationship between the absolute position (BCG) and relative position (PCG).

Figure 5 .
Figure 5.Comparison of allocation results for three methods.

Figure 5 .
Figure 5.Comparison of allocation results for three methods.

Figure 6 .
Figure 6.Comparison of control strength at different center of gravity positions: (a) row-based model; (b) compartment-based model.

Figure 6 .
Figure 6.Comparison of control strength at different center of gravity positions: (a) row-based model; (b) compartment-based model.

Figure 8 .
Figure 8. Passenger weight uncertainty and aircraft center of gravity bias distribution.

Figure 8 .
Figure 8. Passenger weight uncertainty and aircraft center of gravity bias distribution.

Figure 11 .
Figure 11.Effective ranges of numbers of passengers for different TCGs.

Figure 11 .
Figure 11.Effective ranges of numbers of passengers for different TCGs.

Figure 12 .Figure 12 .
Figure 12.Distribution of center of gravity deviations for 20 and 120 passengers.

Figure 12 .
Figure 12.Distribution of center of gravity deviations for 20 and 120 passengers.

Figure 13 .
Figure 13.Variation in standard deviation with the number of passengers.

Figure 13 .
Figure 13.Variation in standard deviation with the number of passengers.

Figure 15 .
Figure 15.Change in the mean value of the distribution concerning the number of persons: (a) numerical upper boundary; (b) numerical lower boundary.

Figure 15 .
Figure 15.Change in the mean value of the distribution concerning the number of persons: (a) numerical upper boundary; (b) numerical lower boundary.

Figure 15 .
Figure 15.Change in the mean value of the distribution concerning the number of persons: (a) numerical upper boundary; (b) numerical lower boundary.

Figure 16 .
Figure 16.Changes in the degree of variance of the distribution with the number of people: (a) numerical upper boundary; (b) numerical lower boundary.

Figure 16 .
Figure 16.Changes in the degree of variance of the distribution with the number of people: (a) numerical upper boundary; (b) numerical lower boundary.

Table 1 .
Comparison of allocation results for three methods with varying numbers of passengers.

Table 2
demonstrates how much faster the compartment-based model produces results than the row-based model.The speed of the row-based model is due to its search space being over ten times larger than that of the compartment-based model.

Table 2 .
Comparison of solution time for two models with varying numbers of passengers.

Table 2 .
Comparison of solution time for two models with varying numbers of passengers.

Table 3 .
Confidence intervals for the center of gravity deviation.Because we only use a set maximum distribution in this table to encompass the other case distributions, it only has the first three rows of values representing the range of center of gravity offsets at different confidence levels. 1

Table 3 .
Confidence intervals for the center of gravity deviation.
display the confidence intervals that we computed outside the upper and lower bounds of the model's effective allocation range.

Table 4 .
Confidence intervals for the center of gravity deviation (numerical upper boundary).

Table 4 .
Confidence intervals for the center of gravity deviation (numerical upper boundary).