Enhancing the Performance and Durability of Commercial Vehicle Cargo Box Frames through Modal Analysis

: The cargo box frame (CBF) is the main structure of a commercial vehicle designed to handle loads and components during travel. The chassis is subject to vibrations caused by rough roads and the components mounted on it. This study proposes a procedure for analyzing and validating CBF structures using a combination of non-destructive modal analysis and ﬁnite element analysis to investigate the vibration characteristics of the four-wheel CBF, including its natural frequency and mode shapes. The CBF’s response to various load conditions, including stress distribution and displacement, was analyzed. The results show that the actuation frequency can affect a truck’s chassis due to the CBF’s natural frequency falling within the excitation range. The resulting mode shape can improve CBF strength, reduce weight, identify defective welds, and determine optimal mounting locations based on the center of gravity (CG) for components such as side-swing doors and cold room panels.


Introduction
The consideration of actual welding processes can have an impact on the complexities and uncertainties involved in the Finite Element Analysis (FEA) of structures.FEA utilizes welding models to approximate the effects of the welding process.However, these models may not accurately replicate the actual welding process, resulting in discrepancies between FEA results and the actual welded structure.Determining accurate welding parameters for FEA analysis can be challenging, especially in the case of complex welded structures.Welded joints may possess different material properties compared to the base structure.FEA typically assumes constant material properties, which may not accurately represent the variations in material properties caused by welding.The welding of multiple components can influence the stiffness and natural frequency of the structure.To address these limitations, it is crucial to refine the welding models employed in FEA and validate them using testing data.Additionally, conducting iterative testing and adjusting parameters based on experimental data can enhance the accuracy of FEA analysis for structures with multiple welded components.
Numerous vehicle owners opt to modify or alter their registered vehicles.However, in the absence of proper control or supervision to ensure compliance with automotive engineering principles and appropriate safety standards, this can pose a threat to both the vehicle owner and other road users.Additionally, such modifications may create opportunities for illegal activities.Therefore, it is essential to establish specific criteria to permit vehicle modifications or alterations to ensure compliance with international standards and guarantee safety.The Land Transport Act of B.E. 2522 (1979) is among the regulations that Thailand has utilized to define the following terms [1]:

•
Land transportation refers to the transportation of people, animals, or goods by road using vehicles.

•
Regular transportation refers to transportation for hire on routes designated by the Transport Committee.

•
Small-scale transportation refers to transportation of people or goods, or a combination of both, for hire on routes designated by the transport committee, using vehicles with a combined weight of no more than 4000 kg.

•
Personal transportation refers to transportation of goods for one's own business using vehicles that weigh more than 2200 kg.

•
A pickup truck is defined as a vehicle that has a solid, box-shaped section for carrying cargo with a permanent roof.The cargo area and the driver's compartment may be integrated together or separated, and there may be doors for loading on the side or rear.It is classified as the second type of vehicle used for transporting animals or goods from among the nine types.
The regulations cited in this passage include the Vehicle Act (No. 3) of B.E. 2525 (1982), as well as Sections 12 and 14 of the Vehicle Act, B.E. 2522 (1979), and Section 4 of the Land Transport Department Regulation on Criteria for Requesting and Granting Permits for Modified Vehicles under the Vehicle Act, B.E. 2562 (2019).According to these regulations, "vehicle" refers to personal cars and motorcycles, "chassis" refers to the original manufacturer's framework that supports the weight and cargo of the vehicle, and "body" pertains to the driver's seat and passenger compartment, but excludes the cargo area.Vehicles that have undergone reinforcement, modification, or installation of additional parts are considered to be unlawfully altered, specifically pickup trucks.Nonetheless, this can be rectified by informing the Department of Land Transport where the vehicle is registered about the modification.Failing to report such modifications would be considered a breach of the Land Transport Act of B.E. 2522.As of September 30, 2022, the total number of trucks in Thailand was 1,217,179.Compared to September 30, 2021, there was an increase of 397,122 trucks that are not fixed-route, representing 32.63% (an increase of 5.57%).The type of truck that saw the highest increase was pickup trucks at 6.96%, followed by semitrailers with a total of 98,208 vehicles (representing 24.73% of the total) and pickup trucks with a total of 94,085 vehicles (representing 23.69% of the total) [2].
The cargo box trucks (CBTs) of the single cab type are designed for chassis applications in Thailand, with a focus on providing a dependable service for an extended period.These trucks are adaptable to meet various chassis requirements due to their diverse cargo box sizes and configurations.Moreover, the CBFs possess the flexibility and capacity to produce highly customized truck bodies that can satisfy the most challenging applications.The Body Builder's Guide for HILUX 2020 by Toyota Motor Corporation [3] states that all bodybuilders must prove that their body modifications do not affect the fundamental functions of the base vehicle.Additionally, they must conduct comprehensive research to ensure that any changes made to the standard parts do not create functional or safety issues from both technical and safety perspectives.There are discussions on safety factors in the structural integrity assessment of components with defects [4], focusing on the structural analysis of ladder chassis for higher strength [5], and examining the calculation of safety factors according to the Eurocodes [6].Gulvanessian, Calgaro, and Holicky provide a designer's guide to EN 1990, EUROCODE, which serves as the basis for structural design [7].
The body modifications or alterations should not obstruct the forward area of view, damage the chassis frame, or cause an imbalance in weight between the right and left wheels.To avoid the local concentration of load on the chassis frame, the body-building job must be conducted in a way that distributes the load evenly over the frames, ensuring that all the wheels are positioned on the same plane without distorting the frame.Additionally, the materials and parts used in the body-building or alteration work should be designed and fabricated to facilitate the inspection and maintenance of the chassis parts once they are mounted on the vehicles.There are specific limitations on the length, width, height, and weight of the parts installed according to each specific base vehicle, and it is necessary to comply with any additional restrictions that may be enforced in each country.
The research articles cover various aspects of modal analysis and its applications in different fields.Modal analysis is explored in various contexts within the collection of articles.It encompasses powertrain modal analysis for vehicle drivability studies [8] and modal analysis applied to a conceptual microsatellite design featuring perforated structural components for mass reduction [9].The utilization of the digital image correlation technique is presented for modal analysis [10], along with an exploration of vibration fatigue analysis of carbon steel coil springs under different road excitations [11].Furthermore, novel modal testing methods designed for structures rotating in water are assessed [12].Additionally, the research focuses on the operational modal analysis of an axial compressor rotor and casing system, aiming to identify a digital twin online [13], as well as the modal analysis of lithium-ion batteries for electric vehicles [14].Strain modal testing using fiber Bragg gratings is investigated specifically for automotive applications [15].
In the automotive industry, modal analysis technology plays a significant role.It is extensively applied in studying the dynamic response analysis of vehicle suspension systems, load reverse engineering [16], and vibration fatigue analysis for designing rapid box wagons [17].The articles also delve into the vibration characteristics of a full-side open boxcar [18] and the effects of adding an auxiliary chassis to a 6-ton truck [19].Overall, this research highlights the importance of modal analysis and vibration analysis in various engineering fields.It showcases a wide range of techniques and applications that aid in identifying potential structural weaknesses, and enhancing durability, reliability, and performance.Specifically, the presented article introduces an analytical modeling approach to enhance structural analysis by understanding acoustic emission signals in thin-walled objects [20].It proposes a novel intelligent method for bearing fault diagnosis that employs EEMD permutation entropy and GG clustering to improve accuracy and efficiency [21].Furthermore, it introduces a new online operational modal analysis method for vibration control in linear time-varying structures, offering real-time monitoring and adjustment capabilities [22].Additionally, the article investigates the application of empirical mode decomposition and artificial neural networks for automatic bearing fault diagnosis based on vibration signals, aiming to achieve reliable and efficient fault detection [23].
Understanding the modal behavior of a car chassis is crucial for analyzing its dynamic behavior.The vibrational characteristics of a vehicle depend on its stiffness and mass distribution.Global bending and torsional vibration mode frequencies are often used to evaluate the structural performance of a vehicle.Bending and torsional stiffness are critical for the vibrational behavior of the structure, especially its first natural frequency.Therefore, to enhance the performance and durability of the entire CBF, a combined simulation and test-based method is used for conducting dynamics analysis of the CBF, based on simulation analysis and vibration test results.
This paper is structured as follows: Section 2 presents the research methodology and materials used.Section 3 analyzes the results and provides validation for both simulation and experimental performance analysis.Section 4 interprets the differences found in the requirement expression characteristic analysis.Lastly, Section 5 presents the conclusions derived from the proposed technique.

Non-Destructive Modal Analysis
Modal analysis is the process of determining modal parameters of mechanical systems that describe their dynamic behavior.It can be approached theoretically, numerically, or experimentally.Theoretical modal analysis involves solving the eigenvalue problem of the mathematical model of a mechanical system.The equations that describe the linearized motion of a general form in a particular configuration are provided by where M, D and K are the mass matrix, the damping matrix and the stiffness matrix, respectively.The vector q(t) contains the degrees of freedom (DOF) of the system.Equation ( 1) is a set of differential equations in matrix form for the dynamic response of a structure modelled with a finite number of degrees of freedom.An eigenvalue analysis, commonly known as the most prevalent form of dynamic analysis, is often performed.Apart from determining the frequencies, it is also essential to investigate the mode shapes that occur at the natural frequencies.These mode shapes represent the undamped free vibration response of the structure resulting from an initial disturbance from its static equilibrium position.To obtain this solution, the damping and applied force terms in the general equation are neglected.Subsequently, it is assumed that each node experiences sinusoidal functions (q = A sin(ωt)) with peak amplitudes specific to that node.The number of eigenvalues or natural frequencies in the model is equal to the total number of DOFs.Each eigenvalue or frequency is associated with a corresponding eigenvector or mode shape.Since none of the eigenvectors can be zero vectors (A = 0), the equation that needs to be solved is as follows: It is important to mention that our typical focus lies on the first few eigenvalues of the model.This is because the finite element model, being an approximation of the structure, tends to yield inaccurate higher eigenvalues and vectors.According to the theoretical solution, the structure exhibits indefinite vibrations in various mode shapes.However, due to the presence of damping in all structures, these vibrations ultimately decay over time.
Modal parameters are determined through measurements, which is an experimental approach.This involves recording the time histories of excitation force and corresponding response of the system.The system's dynamic behavior can be described in the frequency domain using the frequency response function (FRF) or in the time domain using the impulse response function (IRF).These functions are calculated from measurement data using the Inverse Fourier Transform.Impact testing has become the most popular method for experimental modal analysis as FRF measurements can be computed using a FFT analyzer.Impact testing is a fast, convenient, and low-cost way of finding the modes of machines and structures.
The FRF is defined as the ratio of output signal to input signal.There are three primary types of FRFs based on the type of response parameter being analyzed.These response parameters can be displacement, velocity, or acceleration in relation to frequency.A set of all the FRFs is called frequency response function matrix.Furthermore, the Laplace transform (s = iω) and neglecting initial conditions in the frequency domain of Equation ( 1) are given by where −1 is called the transfer function matrix.The definition of one element of FRF matrix for receptance is as follows: where λ r is eigenvalue of the r-th mode, φ jr is the j-th element of the r-th natural shapes vector, and n is number of modes.The vibration modes are determined from the FRF or IRF by applying special estimation procedures.Each mode is characterized by the natural frequency, damping ratio and mode shape.These characteristics, commonly called modal parameters, contribute to the understanding and utilization of proper orthogonal decomposition in the analysis and modeling of vibrational behavior in mechanical and structural systems [24,25].

Finite Element Modeling
This methodology involves the FEA of CBF followed by design.The pre-processing stage in the simulation involves CAD modelling, meshing and applying loads and boundary conditions.The solution stage involves element matrix formulations, and the assemblage of global stiffness elements followed by inversions and multiplications.The final stage is post-processing, which involves viewing and editing the results.
The CAD model of the CBF was created in the SolidWorks referring to the frame drawing and body attachment holes of a single cab, as shown in Figure 1.This model is simulated using the ANSYS software, V18.1 version.The overall view orientation is shown in Figure 2. Limitations are set for the length (L), width (W), height (H) and weight of the parts mounted according to the particular base vehicle.Any restrictions which may be imposed in each country should be complied with.The simulation used in this research is the modal analysis for model validation to compare with the EMA, after which this CAD model is used to evaluate the effect of static structural analysis.
r jr vector, and n is number of modes.The vibration modes are determined from the FRF or IRF by applying special estimation procedures.Each mode is characterized by the natural frequency, damping ratio and mode shape.These characteristics, commonly called modal parameters, contribute to the understanding and utilization of proper orthogonal decomposition in the analysis and modeling of vibrational behavior in mechanical and structural systems [24,25].

Finite Element Modeling
This methodology involves the FEA of CBF followed by design.The pre-processing stage in the simulation involves CAD modelling, meshing and applying loads and boundary conditions.The solution stage involves element matrix formulations, and the assemblage of global stiffness elements followed by inversions and multiplications.The final stage is post-processing, which involves viewing and editing the results.
The CAD model of the CBF was created in the SolidWorks referring to the frame drawing and body attachment holes of a single cab, as shown in Figure 1.This model is simulated using the ANSYS software, V18.1 version.The overall view orientation is shown in Figure 2. Limitations are set for the length (L), width (W), height (H) and weight of the parts mounted according to the particular base vehicle.Any restrictions which may be imposed in each country should be complied with.The simulation used in this research is the modal analysis for model validation to compare with the EMA, after which this CAD model is used to evaluate the effect of static structural analysis.The physical and mechanical properties are shown in Table 1.The software is used to determine the dynamic characteristics such as natural frequencies and the corresponding mode shapes of the CBF.The hot rolled structural steel (SS400) material is defined in the property of the static structural and modal modules in ANSYS software.The following in Table 1 shows the Thai Industrial Standard (TIS) of SS400 material mechanical properties, such as the yield point or yield strength (yield stress), tensile strength and elongation.Convergence studies about element size were preformed and The physical and mechanical properties are shown in Table 1.The software is used to determine the dynamic characteristics such as natural frequencies and the corresponding mode shapes of the CBF.The hot rolled structural steel (SS400) material is defined in the property of the static structural and modal modules in ANSYS software.The following in Table 1 shows the Thai Industrial Standard (TIS) of SS400 material mechanical properties, such as the yield point or yield strength (yield stress), tensile strength and elongation.Convergence studies about element size were preformed and optimum element sizes were selected as 10 mm.In the finite element model of the CBF, 552,400 elements of hybrid grid type were used.Hybrid grids combine the advantages of structured (quadrilateral) and unstructured (triangular) elements, enabling efficient meshing in areas where structured grids are suitable, as well as handling complex regions with unstructured elements.They are particularly useful for the intricate or highly distorted geometries of the CBF.The initial testing was subjected to a free-free boundary condition, where it was supported by a chain hoist.This condition was then modeled in Ansys Workbench, simulating a free-free condition without any imposed restrictions on the boundary conditions.Specifically, none of the DOFs under the body attachment holes of the CBF were constrained or selected.

Securing Basic Performance and Safety Requirements
For the purpose of ensuring safe operation, an even distribution of the cargo weight across the entire vehicle is assumed.The ratio provided below should be employed to ascertain the load exerted on the front axle: where Fr is the front axle mass and T is the total vehicle mass.The mass checks are measured the curb mass of the built or altered vehicle.Take a measurement of the front axle mass (Fr) and the rear axle mass (Re) separately.Each value must not be more than their respective tolerances.The relation between the overall vehicle mass and maximum allowable axle weight is as follows: where the A defines the special equipment, accessories and permanent attachment.The B is the total mass of passengers, cargo and baggage.

Experimental Investigation and Results
The following are the steps for analyzing and designing the structure: • Modal Test and Model Validation: these are the steps for modeling and analyzing the CBF structure: a. Measure various dimensions of the commercial vehicle and create a threedimensional structure of CBF using SolidWorks software. b.
Test the main structure characteristics using a modal testing technique by classifying the behavior of different frequencies.c.
Simulate vibration behavior at different frequencies using a finite element method in Ansys software.d.
Compare the accuracy of test results with the model to use in statistical analysis.e.
Simulate damage analysis with Ansys software under various conditions.Change the supporting metal structure for the refrigerated cabinet.

•
The secure running safeties: compare the results before and after improvements.

Modal Test and Model Validation
This section describes the experiments required to perform a modal analysis on the CBF.Non-parametric system identification is applied to the input and output spectra to obtain a set of FRFs.The first step in the procedure is to measure the impulse response of a system, by means of an impact test.This is performed using an impact hammer and tri-axial acceleration sensors.A data acquisition module is utilized to measure the input and output signals between each excitation point and measured output.Analog signals are converted into discrete time signals by the module, which are then communicated to the computer at fixed time intervals specified by the sampling frequency.These input and output signals are transformed to the frequency domain by means of the Discrete Fourier Transform.These transformed signals are used to calculate the FRFs.
The following equipment is required to perform an impact test: 1.
The free boundary condition is subject to limitations imposed by the natural frequencies of the rigid body modes.To investigate the dynamic behavior of a CBF, it is suspended from a chain hoist, creating one or more rigid body modes based on the stiffness of the supporting materials and the total mass of the structure.If the natural frequencies of these rigid body modes are significantly different from the first natural frequency of the structure, the measured FRF data should not be affected by this boundary condition.

2.
The 3D motion at each test point is desired in the resulting mode shapes.A roving triaxial accelerometer is used and the CBF is impacted at a fixed DOF with the hammer.The locations in the structural testing refer to the 3D wireframe model (208 points) of CBF, shown in Figure 3.The tri-axial accelerometer must be simultaneously sampled together with the force data; a four-channel FFT analyzer is required instead of a two-channel analyzer.

3.
The data-acquisition device DT9837B is used with the IEPE signal conditioning unit.This device is a high-accuracy dynamic signal acquisition module that is ideal for vibration measurements.It contains four synchronized 24-bit sensor inputs that provide a data stream that is matched in time, which can be used for signal analysis.During the tests, only one SIMO transfer is measured, using one impulse force hammer (type 9722A500) and one tri-axial accelerometer (type 8763B100) from KISTLER.
3. The data-acquisition device DT9837B is used with the IEPE signal conditioning unit.This device is a high-accuracy dynamic signal acquisition module that is ideal for vibration measurements.It contains four synchronized 24-bit sensor inputs that provide a data stream that is matched in time, which can be used for signal analysis.
During the tests, only one SIMO transfer is measured, using one impulse force hammer (type 9722A500) and one tri-axial accelerometer (type 8763B100) from KISTLER.

4.
Sampling frequency, also known as sampling rate, refers to the number of samples taken per unit of time during the process of digitizing an analog signal.It represents the rate at which the continuous analog signal is discretized into discrete samples.The choice of an appropriate sampling frequency depends on the characteristics of the signal being sampled and the specific requirements of the application.Higher sampling frequencies capture more detail and allow for the accurate representation of high-frequency components, but they also require more storage and computational resources.Let us define the sampling frequency for our testing based on the number of samples and frequency span in the FRF setup as 4096 lines and 1024 Hz.According to Nyquist's theorem, we can only accurately measure frequencies up to half of the sampling rate.Therefore, the sampling frequency should be set at least twice the maximum frequency component present in the signal.So, our spectral line resolution is calculated as half of the sample rate divided by the number of lines, resulting in 0.5 Hz.The window calculation time is 2 s.Each experiment is repeated at least three times, and the best three measurements are used to determine the averaged FRF for the corresponding SIMO transfer.
Setting up the signal processing: During the measurement, choose the frequency domain and select the FRF (transfer function) for calculation.Use the linear method for averaging and set the number of repetitions to three times.Regarding sampling, set the number of samples to 2048 in the time domain.In the frequency domain, set the maximum frequency to 1280 Hz.Other values will be calculated based on these two parameters.For the trigger configuration, select the trigger type and use Force Channel 1 (connected to the acquisition device, DT9837B) as the reference port.Choose the +slope of the trigger and set the trigger lines and lock lines to observe and adjust the time signal from the hammer impact, as shown in Figure 4a.In this test, use 0.4% of the channel voltage and set a pre-trigger delay of 10 samples.In the setup for the channels spreadsheet, activate four channels and select AC Coupling for one of the piezoelectric force hammers and for the other three channels, connect them to accelerometers.Set the force unit as lbf and the acceleration unit as g.Enter the sensitivity constant of 50 mV/unit for all the transducers (using the instrument's datasheet).During impact testing, use the force window for the hammer channel and the Hanning window for the acceleration channel.
trigger delay of 10 samples.In the setup for the channels spreadsheet, activate four channels and select AC Coupling for one of the piezoelectric force hammers and for the other three channels, connect them to accelerometers.Set the force unit as lbf and the acceleration unit as g.Enter the sensitivity constant of 50 mV/unit for all the transducers (using the instrument's datasheet).During impact testing, use the force window for the hammer channel and the Hanning window for the acceleration channel.Precautions to be taken during signal acquisition: It is essential to carefully examine the response results each time an impact is applied to ensure that the signal patterns are suitable for subsequent signal processing, as shown in Figure 4b.This is particularly important because 208 data points are collected for the FRF, and all three directions of response are measured using accelerometers, resulting in a total of 624 functions.Furthermore, it is not feasible to collect response data in the time domain.Therefore, it is necessary to observe the response results only during the data acquisition process.For instance, let us consider the measurement at point 80 in the -Y direction (denoted by a negative sign, indicating the opposite direction of the reference installation of the accelerometer) compared to point 131 in the +X direction.The response results in the time domain are depicted in Figure 4b.In this specific test, point 80 in the -Y direction is selected as the impact location using a hammer, and roving accelerometers are employed to measure the response at the designated position, which is point 131 in the +X direction on the structure.
5. The responses of the CBF are acquired in the tri-axial accelerometer and hammer data acquisition system and then transferred into MEscope software.For parameter estimation of the CBF, the experimental modal analysis (EMA) software [26] is used.The CBF tests are conducted under predetermined environmental conditions.The frequency span of 0-25 Hz is chosen for the model validation.Precautions to be taken during signal acquisition: It is essential to carefully examine the response results each time an impact is applied to ensure that the signal patterns are suitable for subsequent signal processing, as shown in Figure 4b.This is particularly important because 208 data points are collected for the FRF, and all three directions of response are measured using accelerometers, resulting in a total of 624 functions.Furthermore, it is not feasible to collect response data in the time domain.Therefore, it is necessary to observe the response results only during the data acquisition process.For instance, let us consider the measurement at point 80 in the -Y direction (denoted by a negative sign, indicating the opposite direction of the reference installation of the accelerometer) compared to point 131 in the +X direction.The response results in the time domain are depicted in Figure 4b.In this specific test, point 80 in the -Y direction is selected as the impact location using a hammer, and roving accelerometers are employed to measure the response at the designated position, which is point 131 in the +X direction on the structure.

5.
The responses of the CBF are acquired in the tri-axial accelerometer and hammer data acquisition system and then transferred into MEscope software.For parameter estimation of the CBF, the experimental modal analysis (EMA) software [26] is used.The CBF tests are conducted under predetermined environmental conditions.The frequency span of 0-25 Hz is chosen for the model validation.

6.
The modal identification is performed using a technique of extraction in the peak picking method from the power spectral densities.This technique provides estimations of the natural frequencies and mode shapes.The first six mode shapes obtained from this method are shown in the next part to compare with the simulation, and Figures 5-7 illustrates the demonstration of the motion snapshot of the CBF at 4.31 Hz to 23.8 Hz, respectively.If a CBF is excited near one of its frequencies, the corresponding mode shape will dominate.However, there will still be a small contribution from all other modes present in this response.
The road excitation from the transmission system has typical values varying from 0 to 100 Hz.In practice, diesel engines are known to have an operating speed of about 8 to 33 Hz [19].In the low-speed idling condition, the excitation frequencies are about 8 to 10 Hz.The main excitation is at low speeds, when the truck is in the first gear during the idling condition of engine.At higher gears or speed, the excitation to the chassis is much less.
estimations of the natural frequencies and mode shapes.The first six mode shapes obtained from this method are shown in the next part to compare with the simulation, and Figures 5-7 illustrates the demonstration of the motion snapshot of the CBF at 4.31 Hz to 23.8 Hz, respectively.If a CBF is excited near one of its frequencies, the corresponding mode shape will dominate.However, there will still be a small contribution from all other modes present in this response.The natural frequencies of the CBF were obtained with a range between 0 and 25 Hz.The first six modes were examined.It can be seen from Figures 5-8 that there was a good agreement between the mode shapes in FEA and EMA. Figure 8 illustrates the motion snapshot of the CBF.Identified frequencies from the FEA and corresponding mode shapes obtained from EMA are summarized in Table 2.It obtained similar results, with differences ranging from 0.55 to 1.87 for the first six modes.There were some differences between the results obtained from the FEA and EMA.The data analysis yielded a correlation coefficient of R = 0.9875, a coefficient of determination (r-square) of 0.9752, and a standard error of 1.1148.These statistical measures assess the relationship between a simulated model using FEA as the independent variable and the test results obtained from EMA as the dependent variable in modal analysis.
The road excitation from the transmission system has typical values varying fro to 100 Hz.In practice, diesel engines are known to have an operating speed of about 33 Hz [19].In the low-speed idling condition, the excitation frequencies are about 8 to Hz.The main excitation is at low speeds, when the truck is in the first gear during idling condition of engine.At higher gears or speed, the excitation to the chassis is m less.
The natural frequencies of the CBF were obtained with a range between 0 and 25 The first six modes were examined.It can be seen from Figures 5-8 that there was a g agreement between the mode shapes in FEA and EMA. Figure 8 illustrates the mot snapshot of the CBF.Identified frequencies from the FEA and corresponding mode sha obtained from EMA are summarized in Table 2.It obtained similar results, w differences ranging from 0.55 to 1.87 for the first six modes.There were some differen between the results obtained from the FEA and EMA.The data analysis yielde correlation coefficient of R = 0.9875, a coefficient of determination (r-square) of 0.9752, a standard error of 1.1148.These statistical measures assess the relationship betwee simulated model using FEA as the independent variable and the test results obtained fr EMA as the dependent variable in modal analysis.It is thought that these differences resulted from the uncertainties in the struct steel weldment (boundary conditions) of types and material properties of steel.In FEA, the properties are selected as a density 2.45 g/cm 3 and a Young's modulus of GPa, respectively.It is thought that these differences resulted from the uncertainties in the structural steel weldment (boundary conditions) of types and material properties of steel.In the FEA, the properties are selected as a density 2.45 g/cm 3 and a Young's modulus of 88.7 GPa, respectively.
In the comparison of natural frequencies in the ANSYS modal and experimental modal, the dominant mode was a lateral bending which occurred in the 4-5 Hz range with maximum translation experienced in alternate between front or rear of the CBF.The second to fifth modes were a torsion in the 8-15 Hz range and a 2-nodal-point lateral bending in the 21-24 Hz range in the sixth mode.Thus, the CBF may experience structural resonance at an idling condition.An increase in mass will reduce the natural frequency, while an increase in the stiffness will increase the natural frequency.Now, the CAD model and material properties can be used to evaluate other analytical models if the mechanical system meets its specifications.

Refinement and Enhancement
The process of performing a static structural analysis involves the following main steps.

•
Model Preparation: Create a model of the structure to be analyzed using SolidWorks software.The details are shown in Figures 1 and 2.

•
Meshing: Divide the model into smaller finite elements through a process known as meshing.Convergence studies were conducted to determine the optimal element size, and it was found that an element size of 10 mm yielded the best results.The finite element model of the CBF utilized a hybrid grid type consisting of 552,400 elements.This modeling approach is similar to the modal analysis simulation.

•
Material Definition: Specify the material properties for each element in the model.This includes mechanical properties such as the Young's modulus of 88.7 GPa and density of 2.45 g/cm 3 , which define how the material will respond to applied loads.

•
Boundary Conditions: Apply appropriate boundary conditions to the model, including constraints, loads, and supports.These boundary conditions mimic the real-world environment in which the structure operates.The structure of the CBF is attached to the chassis of a pickup truck with M12 bolts for fastening at eight points, as shown in Figure 9.The structure is designed to carry a total mass of 2000 kg and is uniformly distributed throughout the structure of the refrigerated container floor, calculated from the maximum number of ice loads of 100 bags, with each bag massing 20 kg.The main structure for carrying the load is the upper beam, which is made of seven square tubes, stacked on the C-channel beams.The square tubes of the structure that carries the weight has a size of 3 × 1 1 /5 inches.

•
Analysis Solution: Run the analysis using the ANSYS solver, which solves the governing equations of static equilibrium to determine the response of the structure under the applied loads.The solver calculates the Total deformation, Equivalent Von-Mises stress, Equivalent Von-Mises strain, Factor of Safety (FOS), and other relevant results for the structure.

•
Analyze and interpret the results obtained from the analysis.If necessary, refine the model, mesh, or analysis settings based on the results obtained.This iterative process allows for fine-tuning the analysis to achieve more accurate and reliable results.The simulation results suggest that the central part of the structure is the most vulnerable to damage when exposed to a load, with a maximum deformation of 0.41 mm.The structure experienced a maximum stress of 85 MPa at the connection point between the C-channel beam of the structure and the chassis.However, this area had an FOS value of 2.95.The location with the highest strain of the structure was at the same point as the location with the maximum deformation, with a size of 4.59 × 10 −4 mm/mm.
results for the structure.

•
Analyze and interpret the results obtained from the analysis.If necessary, refine the model, mesh, or analysis settings based on the results obtained.This iterative process allows for fine-tuning the analysis to achieve more accurate and reliable results.The simulation results suggest that the central part of the structure is the most vulnerable to damage when exposed to a load, with a maximum deformation of 0.41 mm.The structure experienced a maximum stress of 85 MPa at the connection point between the C-channel beam of the structure and the chassis.However, this area had an FOS value of 2.95.The location with the highest strain of the structure was at the same point as the location with the maximum deformation, with a size of 4.59 × 10 −4 mm/mm.To improve the handling performance of the vehicle, a model of the vehicle's structure was created, and the impact of the mass of the side-swing doors on the structure of the vehicle was considered.The force generated by the mass of the side-swing doors, which is about 30 kg, combined with the total mass of the structure, is equal to 325 kg.After considering the weight exerted on all eight legs of the structure, it was observed that the weight on each leg was not equal.When viewed from above and facing the same direction as the vehicle, the shape of the mass distribution on the legs was tilted towards the left-bottom side of the vehicle, where the side-swing doors are installed.Leg 6 weighed To improve the handling performance of the vehicle, a model of the vehicle's structure was created, and the impact of the mass of the side-swing doors on the structure of the vehicle was considered.The force generated by the mass of the side-swing doors, which is about 30 kg, combined with the total mass of the structure, is equal to 325 kg.After considering the weight exerted on all eight legs of the structure, it was observed that the weight on each leg was not equal.When viewed from above and facing the same direction as the vehicle, the shape of the mass distribution on the legs was tilted towards the left-bottom side of the vehicle, where the side-swing doors are installed.Leg 6 weighed 81 kg, which was more than leg 3 on the opposite side that weighed 49 kg.To improve the position of the CG of the structure, it was necessary to move the center of gravity towards the front as much as possible, by shifting the side-swing doors to the front of the vehicle.
As a guideline for improving the position of the CG of the structure, it is advisable to shift the side-swing door installation towards the front of the vehicle and evaluate the subsequent variation in the center of mass, as shown in Figure 10.Move it three levels from the original position which are 5, 10, and 20 cm, respectively.According to the simulation, the position of the CG of the structure and the forces generated on each leg changed, as shown in Table 3.This causes the mass that each leg bears to shift more towards the front-left direction, which is the direction that improves the balance of the structure.The side-swing door installation position was changed by moving it to the front of the car.It was found that the mass at the position where the leg supports the force decreased at positions G and H by approximately 30%, causing the position of the CG of the structure to change, as shown in Figure 11.The displacement of the center of mass of the structure is directly proportional to the extent to which the side-swing door is moved.Based on the simulation results, it can be concluded that moving the side-swing door installation position 20 cm to the front shifts the position of the center of mass of the structure forward by a distance of 32 cm.Positions A to H refer to Figure 11b.Due to the results of the simulation of the prototype structure in the previous topic, it was found that the structure still had a very high FOS.It is necessary to reduce the size of some parts of the structure to reduce the weight of the structure and save production costs.Therefore, it was decided to reduce the size of the additional structural steel by changing the size of the square tube (ST) to a C-light lip channel (CL) and using an angle steel bar to fasten it to the structure, as shown in Figure 12.Adjusting the size caused it to lesser affect the load behavior compared to the original model, as shown in Table 3. Due to the results of the simulation of the prototype structure in the previous topic, it was found that the structure still had a very high FOS.It is necessary to reduce the size of some parts of the structure to reduce the weight of the structure and save production costs.Therefore, it was decided to reduce the size of the additional structural steel by changing the size of the square tube (ST) to a C-light lip channel (CL) and using an angle steel bar to fasten it to the structure, as shown in Figure 12.Adjusting the size caused it to lesser affect the load behavior compared to the original model, as shown in Table 3.
of some parts of the structure to reduce the weight of the structure and s costs.Therefore, it was decided to reduce the size of the additional stru changing the size of the square tube (ST) to a C-light lip channel (CL) and steel bar to fasten it to the structure, as shown in Figure 12.Adjusting the s lesser affect the load behavior compared to the original model, as shown i The structure received a load, along with the side-swing doors that w moving the installation position forward by 20 cm, and we found that th collapsed from the load over a distance of 3.1 mm, as shown in Figure 13 was found that the maximum stress and strain values were 363 MPa a mm/mm, respectively.The position where the maximum strain and stres on the steel beam used to attach the C-light lip channel beam structu structure, as shown in Figure 14.The aforementioned strain value is in a ra the steel to undergo permanent deformation, which poses a significa structure.In order to prevent damage to the angle steel bar (flat steel) whe The structure received a load, along with the side-swing doors that were adjusted by moving the installation position forward by 20 cm, and we found that the structure had collapsed from the load over a distance of 3.1 mm, as shown in Figure 13.In addition, it was found that the maximum stress and strain values were 363 MPa and 1.05 × 10 −3 mm/mm, respectively.The position where the maximum strain and stress occurred was on the steel beam used to attach the C-light lip channel beam structure to the main structure, as shown in Figure 14.The aforementioned strain value is in a range that causes the steel to undergo permanent deformation, which poses a significant risk to the structure.In order to prevent damage to the angle steel bar (flat steel) when clamping the top steel beam to the main structure, the thickness of the flat steel used for clamping the beam was simulated to be 1/8 inch or 3.0 mm.The results of the simulation show the maximum force behavior that occurs on the structure, which is presented in Table 4.
Appl.Sci.2023, 13, x FOR PEER REVIEW 18 top steel beam to the main structure, the thickness of the flat steel used for clampin beam was simulated to be 1/8 inch or 3.0 mm.The results of the simulation show maximum force behavior that occurs on the structure, which is presented in Table 4     The simulation results indicate that altering the size of the flat ste structure's susceptibility to damage, and it was concluded that a flat steel t less than 3 mm is necessary to prevent damage to the structure.However, thickness of the flat steel to more than 3.0 mm does not help to make the stru because the structure will still tear at the joints with the C-light lip chan thinness and increased deformation of the C-light lip channel.During mod 01 and 02, the total mass of the CBF main structure and side-swing door w  The simulation results indicate that altering the size of the flat steel impacts the structure's susceptibility to damage, and it was concluded that a flat steel thickness of no less than 3 mm is necessary to prevent damage to the structure.However, increasing the thickness of the flat steel to more than 3.0 mm does not help to make the structure stronger, because the structure will still tear at the joints with the C-light lip channel due to the thinness and increased deformation of the C-light lip channel.During model simulations 01 and 02, the total mass of the CBF main structure and side-swing door was 325 kg.But when the C-light lip channel was modified in model simulations 03 and 04, the total mass decreased to 319 kg, a reduction of 1.8%.Although this modification caused the FOS value to decrease by approximately 53% from the previous value, the FOS value remained within the acceptable range of the safety standards.In most cases, structural designs include a higher factor of safety, with the exception of the aerospace and automobile industries, where minimizing weight is critical to performance and cost.To achieve this, safety factors are deliberately kept low, typically around 1.25 to 2.0, for the design of structures that must withstand static loads with high confidence, and all available design data are used to ensure their integrity.

Secure Running Safety Comparison
The load acting on the structure during the installation of the side-swing doors of the refrigerated truck container is significant.A side-swing door was installed on the left side of the CBF structure for the ease of transporting ice bags out of the refrigerated container.The installation of the side-swing door leads to an imbalance in the total weight of the truck.Once the modified truck structure was complete, it was weighed to check the weight distribution on each side of the wheels.The weight indicator in the rack version of the Dini Argeo, complete with a connector for linking the platforms, is shown in Figure 15.It was found that the weight on each side of the wheels was not balanced, with the weight falling on the left wheel being greater than the right.The details of the topic, Securing Basic Performance and Safety Requirements, will be shown last.Factors affecting the left side structure of the vehicle were the position of the fuel tank combined with the weight of the side-swing doors, and the structure that supports the door force causing the center of mass of the structure to be skewed to the left and located at the rear of the vehicle.This feature was not beneficial for the stability and safety of the vehicle while driving.Therefore, it was necessary to analyze the center of mass of the structure and adjust the appropriate installation position of the side-swing doors.

Secure Running Safety Comparison
The load acting on the structure during the installation of the side-swing doors of the refrigerated truck container is significant.A side-swing door was installed on the left side of the CBF structure for the ease of transporting ice bags out of the refrigerated container.The installation of the side-swing door leads to an imbalance in the total weight of the truck.Once the modified truck structure was complete, it was weighed to check the weight distribution on each side of the wheels.The weight indicator in the rack version of the Dini Argeo, complete with a connector for linking the platforms, is shown in Figure 15.It was found that the weight on each side of the wheels was not balanced, with the weight falling on the left wheel being greater than the right.The details of the topic, Securing Basic Performance and Safety Requirements, will be shown last.Factors affecting the left side structure of the vehicle were the position of the fuel tank combined with the weight of the side-swing doors, and the structure that supports the door force causing the center of mass of the structure to be skewed to the left and located at the rear of the vehicle.This feature was not beneficial for the stability and safety of the vehicle while driving.Therefore, it was necessary to analyze the center of mass of the structure and adjust the appropriate installation position of the side-swing doors.The maximum allowable load was reduced when the mass of the finished vehicle (C.W. + A = T) increased, while keeping the G.V.W. fixed.Table 5 shows the chassis and cab curb mass of the base vehicle (C.W.) and the gross vehicle mass of the finished vehicle (G.V.W.) for the models GUN122R-BTFXYT3 and -BTFLXT3.The maximum allowable load was reduced when the mass of the finished vehicle (C.W. + A = T) increased, while keeping the G.V.W. fixed.Table 5 shows the chassis and cab curb mass of the base vehicle (C.W.) and the gross vehicle mass of the finished vehicle (G.V.W.) for the models GUN122R-BTFXYT3 and -BTFLXT3.

Discussion
The high correlation coefficient (R = 0.9875) signified a robust positive linear relationship between the simulated model derived from FEA and the test results obtained through EMA.Moreover, the r-square value indicated that approximately 97.52% of the variability in the test results can be accounted for by the simulated FEA model.The standard error of 1.1148 represents the average deviation between the observed EMA test results and the predicted values obtained from the simulated FEA model.A lower standard error signifies a more accurate fit of the model in predicting the EMA test results.These findings demonstrate that the FEA model provides a substantial explanation for the observed variation in the EMA test results.Therefore, in the analysis and validation for assessing the failure behavior, the properties were chosen as a density of 2.45 g/cm 3 and a Young's modulus of 88.7 GPa.
The selection of a CL steel section over a ST section can be attributed to two reasons.Firstly, the CL section offers cost advantages as it requires less material compared to a box-shaped section while still providing structural support.Secondly, the CL section offers flexibility in design and installation, allowing for easier customization and adaptation to different construction requirements.However, it is important to note that the CL steel section's open design reduces the cross-sectional area available for carrying loads, resulting in reduced strength and stiffness.This characteristic makes the CL section more susceptible to bending and deformation compared to the ST section.In contrast, the ST steel section is stronger than the CL steel section due to its geometric properties.The ST section has equal thickness on all sides, enabling a uniform distribution of stress and load-carrying capacity.This uniform stress distribution enhances the box-shaped section's ability to resist bending and deformation more effectively.
The results of the comparison of the mass of a CBF four-wheel vehicle before and after modifying the structure and shifting the CG position of the side-swing door towards the front of the car by 20 cm are presented.According to the manufacturer's standards for the GUN122R-BTFXYT3 and -BTFLXT3 vehicle models, the G.V.W. does not exceed the specified values for the front wheel (1165 kg) and rear wheel (1685 kg), and the total G.V.W. is equal to 2850 kg.From the ratio criteria value, it can be seen that the driving safety is 45.56%, which exceeds the requirement of being greater than or equal to 30%.However, the mass affecting the front wheel before modification differs from that affecting the left rear wheel by 65 kg, while the mass of the left and right front wheels differs only by 15 kg, and the mass of the front and rear wheels differs by 50 kg.This could lead to problems with wheel locking during driving.
Regarding the structure after modifying and shifting the CG position of the side-swing door to the front of the car by 20 cm, the mass affecting the left rear wheel decreased, and the difference between the rear wheels reduced by 40 kg (a decrease of 61.5%).Meanwhile, the difference in mass between the left and right front wheels increased by 20 kg.Additionally, the difference in mass between the front and rear wheels decreased from 190 to 145 kg (a decrease of 31%).Although the total vehicle mass before and after the changes differed by only 5 kg, the ratio criteria for driving safety increased to 46.6% without suspension upgrade, ensuring safe driving performance.The total vehicle mass (T) is 2140 kg, which is less than the total G.V.W. limit of 2850 kg, and the maximum mass of passengers, cargo, and baggage (B value) can reach up to 710 kg.However, if loaded to its full capacity of 2000 kg, combined with the total mass before loading of 2135 kg, the total vehicle mass will exceed the total G.V.W. value designed by the manufacturer, which is only 2850 kg.Therefore, if the total G.V.W. exceeds 2200 kg but is not more than 4000 kg, according to the law, an examination must be conducted to evaluate its impact on the reliability, stability, and safety of operation to ensure compliance with engineering principles.The examination must include plans and calculations, along with details of the inspection in the modified or adjusted parts, conducted by a licensed professional engineer to issue a certification letter.
The structure of the CBF has a maximum strain of 1.05 × 10 −3 mm/mm, a maximum tensile strength of 197 MPa, and an FOS value of 1.27.The CG was moved forward by 32 mm.Based on the analysis results, the following recommendations are suggested for improving the CBF, although it will remain safe to use even when loaded to its full capacity of 2000 kg that supports cold room panels:

•
Use flat steel with a thickness of at least 3.0 mm or 1/8 inch to make the angle steel bar for attaching the channels.

•
Move the side-swing door installation position 20 cm forward from the original position, adjust the spacing, and reduce the lateral channels to a suitable distribution range.

•
Control the welding process to ensure completeness at every point of the weld.

Conclusions
This paper presents the steps involved in conducting a structural analysis to evaluate the strength of a car box frame (CBF) in a four-wheel vehicle.It also focuses on determining the optimal installation location for side-swing doors and cold room panels.To ensure vehicle stability and safety during cargo transportation, the analysis takes into consideration the optimum center of gravity (CG) position.Non-destructive modal methods and finite element analysis were employed to analyze CBF four-wheeled vehicles across various dimensions.The MEscope software was utilized to test the characteristics of the main structure, providing a description of the required equipment for conducting impact operations.This equipment included a conditional motion three-axis accelerometer and a data acquisition module.The analysis process involved several steps, including data collection, processing, mode identification, and result interpretation.By simulating vibrational behavior at different frequencies, a comparison was made using the finite element method in Ansys software to conduct damage analysis under various conditions.Furthermore, the combination of real test results and simulations in both techniques enhanced the accuracy of the analysis and verification process.This paper contributes to the design and manufacture of large CBFs with a significant number of welds, ultimately leading to modified CBF designs that reduce the production costs for the competitive CBF assembly industry.

Figure 2 .
Figure 2. View orientation of CBT and CBF: (a) perspective views of commercial vehicle; (b) CAD modelling.

•
Refinement and Enhancement: Adjust the local center of gravity position of the supporting structure for the refrigerated cabinet.Analyze the results of the model structure to propose solutions: a. Adjust the local center of gravity position of the supporting structure for the refrigerated cabinet.b.

Figure 3 .
Figure 3. Structural testing refer of the 3D wireframe model of CBF.4.Sampling frequency, also known as sampling rate, refers to the number of samples taken per unit of time during the process of digitizing an analog signal.It represents the rate at which the continuous analog signal is discretized into discrete samples.The choice of an appropriate sampling frequency depends on the characteristics of the signal being sampled and the specific requirements of the application.Higher sampling frequencies capture more detail and allow for the accurate representation of high-frequency components, but they also require more storage and computational resources.Let us define the sampling frequency for our testing based on the number of samples and frequency span in the FRF setup as 4096 lines and 1024 Hz.According to Nyquist's theorem, we can only accurately measure frequencies up to half of the sampling rate.Therefore, the sampling frequency should be set at least twice the maximum frequency component present in the signal.So, our spectral line resolution is calculated as half of the sample rate divided by the number of lines, resulting in 0.5 Hz.The window calculation time is 2 s.Each experiment is repeated at least three times, and the best three measurements are used to determine the averaged FRF for the corresponding SIMO transfer.

Figure 3 .
Figure 3. Structural testing refer of the 3D wireframe model of CBF.

Figure 4 .
Figure 4. (a) Setting up the processing using acquisition; (b) force and acceleration signal.

Figure 4 .
Figure 4. (a) Setting up the processing using acquisition; (b) force and acceleration signal.

Figure 6 .
Figure 6.Motion snapshot from relative maximum deviation behavior of EMA: (a) at 10.6 Hz; (b) at 10.9 Hz.Figure 6. Motion snapshot from relative maximum deviation behavior of EMA: (a) at 10.6 Hz; (b) at 10.9 Hz.

Figure 8 .
Figure 8. Motion snapshot from relative maximum deviation behavior of FEA.

Figure 8 .
Figure 8. Motion snapshot from relative maximum deviation behavior of FEA.

Figure 10 .
Figure 10.The positions of the side-swing doors' sliding installation design.

Figure 10 .
Figure 10.The positions of the side-swing doors' sliding installation design.

Figure 11 .
Figure 11.Changes in (a) percent mass at leg support, and (b) CG position at different distances.

Figure 11 .
Figure 11.Changes in (a) percent mass at leg support, and (b) CG position at different distances.

Figure 13 .
Figure 13.The maximum deformation of the C-light lip channel structural frame.

Figure 13 .
Figure 13.The maximum deformation of the C-light lip channel structural frame.

Figure 13 .
Figure 13.The maximum deformation of the C-light lip channel structural frame.

Figure 14 .
Figure 14.The position of maximum stress and strain.

Figure 14 .
Figure 14.The position of maximum stress and strain.
Number 1 to 5 refers to Figure 2b.

Table 2 .
Comparison of natural frequencies found in ANSYS modal and experimental modal.

Table 3 .
The weight at the leg of the CBF structure affecting the chassis (kg).

Support Position Distance of the Forward Shift in the Installation of the Side-Swing Doors Form Original Original Shift 5 cm. Shift 10 cm. Shift 20 cm.
Positions A to H refer to Figure 11b.Appl.Sci.2023, 13, x FOR PEER REVIEW 17 of 23

Table 4 .
Comparing simulation results of reducing the size of the box beam.

Table 4 .
Comparing simulation results of reducing the size of the box beam.

Table 5 .
Comparison of wheel mass with accepted basic performance and safety requirements.

Table 5 .
Comparison of wheel mass with accepted basic performance and safety requirements.