Next Article in Journal
A Theoretical Framework for Optimality Conditions of Nonlinear Type-2 Interval-Valued Unconstrained and Constrained Optimization Problems Using Type-2 Interval Order Relations
Next Article in Special Issue
Design Optimization of 3-DOF Redundant Planar Parallel Kinematic Mechanism Based Finishing Cut Stage for Improving Surface Roughness of FDM 3D Printed Sculptures
Previous Article in Journal
Coping with Challenges in Teaching Foreign Languages to Children with Mild Intellectual Disabilities: Stakeholders’ Perspectives
Previous Article in Special Issue
Identification of Inertial Parameters for Position and Force Control of Surgical Assistance Robots
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multidisciplinary Lightweight Optimization for Front Impact Structure of Body Frame Based on Active and Passive Safety

School of Mechanical and Power Engineering, Zhengzhou University, Zhengzhou 450000, China
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(8), 907; https://doi.org/10.3390/math9080907
Submission received: 22 March 2021 / Revised: 14 April 2021 / Accepted: 16 April 2021 / Published: 19 April 2021
(This article belongs to the Special Issue Mathematical Problems in Mechanical Engineering)

Abstract

:
The Analytic Target Cascading (ATC) is an effective method for solving hierarchical Multidisciplinary Design Optimization (MDO) problems. At the same time, this method suffers from poor convergence and low accuracy, which is caused by the inconsistency of system constraints. In this paper, a novel ATC method based on dynamic relaxation factor is proposed. The dynamic relaxation factor of consistency constraint is added in the system level and is adjusted by the deviation of the linking variables between the levels to ensure the feasible region of the design space. The effectiveness and accuracy of this method are verified by a mathematical example. This method is used to solve the lightweight problem of the trussed front part of the vehicle body frame based on active and passive safety to achieve the collaborative optimization of lightweight trussed frame, crash safety, and aerodynamic characteristics. The important value of the novel ATC method based on dynamic relaxation factor in engineering applications is proven.

1. Introduction

Formula Student Electric China (FSEC) is a competition which consists of static and dynamic events. Five items are included in the dynamic competition: linear acceleration test, 8-word detour test, high-speed obstacle avoidance test, durability test, and efficiency test. Therefore, it is necessary to comprehensively evaluate the acceleration, durability, and body weight of the racing car. Under the premise of ensuring the personal safety of the driver, the overall performance of the racing car’s acceleration, driving range, handling stability, braking, and lightweight is required to be optimal.
The electric racing car in the FSEC is a typical trussed frame structure [1,2,3]. At present, most scholars only consider the collision safety of the body for lightweight research. In [4], Xiong studied the multi-objective lightweight optimization problem of the front impact structure of the automobile body. The crashworthiness model of the vehicle was established by taking the thickness of the automobile body and the materials as the design variables. In [5], Xie used the thickness of the front cabin of electric automobile as the design variable to design the vehicle lightweight, and verified the safety of the automobile by the simulation of the collision. However, for automobiles, the research of lightweight should not only include passive safety; active safety is also an important factor that cannot be ignored [6,7]. The main problem of active safety is how to improve the handling stability, in addition to being improved by the electronic control safety system [8,9,10,11]. In [12], Yu pointed out that the vertical load of the tire on the ground is reduced by the aerodynamic lift on the automobile, and the change of the vertical load has a significant impact on the wheel cornering characteristics, thus affecting the handling stability, as well as the dynamic performance and braking performance of the automobile. Therefore, the research on active safety of racing car can be considered to include aerodynamic characteristics. In [13], the author analyzed the aerodynamic characteristics of Ford F-150 race truck. In order to change the aerodynamic lift coefficient of the whole vehicle, the modeling design of the truck was adopted, so as to increase the downward pressure of the tire during driving and improve the handling stability of the truck during driving. Aerodynamic characteristics are closely related to the shape of the front of the vehicle [14,15,16], and are affected by the position of the key nodes in the front of the vehicle. In [17], the author used the key control points of the head shape surface to establish 12 spline curves according to the streamlined head shape of a high-speed train, and then generated seven spline curves for parametric modeling of the head. For the truss body structure, when the node changes, the racing car weight will be affected by the length of the pipe with the change of the node position. Therefore, FSEC racing car is improved in the process of lightweight, which is conducive to improving the active safety and affecting the weight of the racing car.
In general, active and passive safety collaborative optimization is needed to optimize the overall performance of FSEC racing car. The combination of active and passive safety is an effective way to improve vehicle safety. Although many companies and scholars began to pay attention to the research on the combination of active and passive system safety [18,19], the improvement of handling stability in these studies was carried out from the electronic control system. Study of the body structure is ignored, which increases the crash safety and improves the aerodynamic characteristics while lightening, and realizes the integration of active and passive safety. Therefore, this paper studies the abovementioned complex system with multiple couplings from the perspective of structural considerations. The system needs to consider the influence of many disciplines such as the structure, aerodynamic characteristics, and braking performance of the entire vehicle. It is a typical MDO problem, which was formally proposed by the multidisciplinary design optimization technical committee [20]. It fully considers the coordination mechanism between various disciplines in the process of complex product design. It can shorten the design cycle and obtain better design schemes by obtaining the optimal value of conflicting objectives [21,22,23,24,25].
In the study of multidisciplinary optimization engineering problems, MDO optimization method is the key part of MDO research content. It is used to describe the engineering problems involved mathematically. The selection of different types of optimization methods has a great impact on the complexity, accuracy, and calculation time of solution. So far, MDO optimization method has been developed from single-level optimization method to multi-level optimization method. The single-level MDO method is integrated into the system level for design optimization by the design variables and constraints of various sub disciplines, which results in the complexity of the system level. It is not optimized for large-scale engineering cases [26,27,28]. Therefore, multi-level MDO methods are often used in MDO. It can be decomposed into a number of relatively simple suboptimization problems according to disciplines. The system-level control of global design variables is adopted, and the local design variables are controlled by the subsystem. The difference in optimization results is coordinated by the system level; for example [29,30,31], Concurrent Subspace Optimization (CSSO), Collaborative Optimization (CO), Bilevel Integrated System Synthesis (BISS), etc. Many new MDO methods have been developed by scholars on the basis of existing MDO optimization methods. In [32], Michelena proposed an ATC method based on component hierarchical decomposition for large-scale system problems according to their structure. This theory provided a reliable method for solving large and complex problems. In [33], Yao considered structural layout, shape, and size optimization, and proposed a Concurrent Subspace Collaborative Optimization method (CSCO). Many design variables do not need to be optimized simultaneously using this method, which significantly reduces the degree of nonlinearity in structural analysis. In [34], Shin and Park proposed the MDO Based on Independent Subspaces (MDOIS) algorithm, which uses a system analysis method to coordinate the coupling relationship between subsystems.
Among them, ATC is used to solve the problem of lightweight and active passive safety comprehensive optimization design of truss body. System decomposition is used to divide a complex system into smaller subsystems and determine the interaction between them. After decomposition, the different disciplines in the subsystem level are connected with each other by linking design variables, so system coordination must be performed to ensure the convergence of all linking design variables, and the inconsistency of variables between disciplines needs to be considered [35]. In view of the inconsistency between subsystem level, Brauhn [36] used the inequality constraint form to relax the system-level consistency constraints in the CO method, which improved the convergence effect of the CO method. In [37], Li proposed a dynamic relaxation factor method based on the relaxation factor constraint, so that the size of the relaxation factor can be adjusted according to the dynamic inconsistency between disciplines. His theory provides a powerful idea for the research in this paper and is a very useful algorithm. It is used to ensure that the relaxation factor is not only large, but also small, which leads to the failure to meet the relaxation requirements among disciplines and leads to no feasible solution.
In recent years, the intermixing of algorithms has provided a good reference for constructing new algorithms and improving the efficiency of the original ones. In [38], new metaheuristic algorithms are proposed by inserting the model F1 in Slime Mold Algorithm (SMA) and other representative algorithms. In [39], Zapata built a 3D structure using Self-Assembly Algorithm With Particle Swarm Optimization (SAPSO) which is a combination of two emerging techniques (self-assembly and PSO) with the purpose of improving the process of construction of 3D structures and the convergence speed based on decreasing the number of iterations for building any structure. In [40], Roman proposed to mix the Compact Form Dynamic Linearization–Proportional–Derivative Takagi–Sugeno Fuzzy Algorithm (CFDL-PDTSFA) and Virtual Reference Feedback Tuning VRFT to get a new data-driven algorithm of CFDL-PDTSFA-VRFT.
Therefore, in this paper, the idea of dynamic relaxation factor in the CO model is introduced into the ATC model, and the consistency of the linking design variables between the same levels of ATC is used as the dynamic constraint condition. The size of the relaxation factor is adjusted according to the inconsistency of the linking design variables, taking into account the iterative calculation in the optimization process. At the same time, the problem of choosing the size of the relaxation factor is well avoided. However, as the number of iterations increases, the inconsistent information of the linking design variables may be too small and cause the problem to converge incorrectly, so a fixed term is introduced to constrain the lower limit of convergence to ensure the solution feasibility. Then, the dynamic relaxation factor constraint ATC method is used to construct the active and passive safety MDO design platform for the front impact structure of the racing car body, and the polynomial response surface model is used to obtain the analysis models of each level [41] to reduce the calculation time, realize the lightweight of the body frame, and improve the active and passive safety of the whole racing car. It is confirmed that the dynamic relaxation factor ATC method is suitable for practical engineering problems.

2. Materials and Methods

2.1. Multidisciplinary Optimization Model

The CO model is widely used in the MDO method. The main idea of CO is to decompose the complex engineering project system design problem into several subsystem level problems. The design optimization of each subsystem is completed independently without considering the influence of the design results of other subsystems on the system. The incompatibility and inconsistency between subsystems are coordinated by the top system level [31]. It has the characteristics of parallel computing and reduces the difficulty of solving high-order nonlinear equations. However, the algorithm framework is a two-level structure, as shown in Figure 1, which cannot be widely used in all kinds of multi-layer engineering problems, such as automobile manufacturing.
The ATC method originated from the design of automotive products [32]. The design targets are continuously decomposed from the upper level to the lower level. At the same time, the responses of all levels are continuously fed back from the lower level to upper level. Each unit problem is solved independently; the upper and lower levels overlap for optimization. The optimization process stops when the convergence conditions are satisfied. Therefore, ATC can be applied to infinitely divided i (i = 1, 2, ..., N) level problems. Figure 2 shows the three-level decomposition system.

2.2. ATC Model Description

For the automotive industry, a large number of design variables of complete vehicles, systems, subsystems, and parts need to be optimized. The original design problem (vehicle level) can be interpreted as finding a design scheme to minimize the deviation between all design targets and responses on the premise of meeting all constraints, and its mathematical description is shown in Equation (1) [32]:
M i n P R                           W h e r e       R = r x                                           s . t .       g i x , y 0     i = 1 , 2 , , n                     h i x , y = 0     j = 1 , 2 , , m   x m i n x x m a x y m i n y y m a x
The optimized objective function in the Equation (1) is defined as the deviation between the required target P and the system response R obtained by the analysis model. r x , g i x , y , and h i x , y are inequality and equation constraints, x is design variable,   y is a linking design variable.
The abovementioned initial design problem is decomposed, and the original problem is divided into multiple hierarchical structures. When using ATC decomposition, the initial targets and constraints in the system level will be decomposed into multiple subsystem spaces. The optimization target P i j at level i is the minimum response deviation calculated by the upper system ( i − 1) and the lower system ( i + 1). The optimization target P i j and the analysis model r exist in the target cascading structure. The design variables and parameters, as well as the lower level response, are occupied by the analysis model and return the response of the design problem. Response and linking variables are introduced to describe the interaction between the vertical and horizontal subsystem problems. The response is defined as the output of an analysis model, and the linking variable is defined as a common design that exists between two or more design models. Therefore, the general flowchart of ATC is described [36], as shown in Figure 3.
According to the data transfer flow chart shown in Figure 3, the general ATC model is proposed in this article, as shown in Equation (2).
Level i:
{ P i j :     m i n ω i j R R i j R i j U + ω i j y y i j y i j U + ε R + ε y r e s p e c t   t o   x ˜ i j ,   y i j ,   y i + 1 j ,   R i + 1 j ,   ε R ,   ε y w h e r e ,   R i j = r i j R i + 1 j , x ˜ i j , y i j s . t .         ω i j R i = 1 n j = 1 c i + 1 k | | R i + 1 j R i + 1 j L | | 2 2 ε R         ω i j y i = 1 n j = 1 c i + 1 k | | y i + 1 j y i + 1 j L | | 2 2 ε y g i j R i j , x i j , y i j 0 h i j R i j , x i j , y i j = 0 x ˜ i j m i n x ˜ i j x ˜ i j   m a x   , y i j m i n y i j y i j m a x
where P i j is the j-th subsystem optimization problem obtained by the upper level of i   i = 1 ,   2 ,   , n ,   P i j = k 1 , k 2 , , k P i j . R is defined as the response value of the analysis system at this level. R U , R L are the target values passed by the upper and lower systems, respectively. y   is   defined   as   the linking variable between the same level. y U ,     y L are the linking variables transmitted by the upper and lower systems, respectively. ω   represents the weight coefficient value of the response or the linking variable and the matching target importance. ε R is obtained by the deviation between the target value and the response value. The number of subproblems in the ( i + 1)th level is defined as c i + 1 k (k = 1, 2, ..., m). ε y is obtained by the deviation between the linking variable and the calculated value. r represents the analysis model of the level. x ˜ i j m i n , x ˜ i j m a x are used as the lower limit and upper limit of the local design variable x ˜ i j .

2.3. The Introduction of Dynamic Relaxation Factor in ATC

In the standard ATC optimization method, the deviation between the target and the response is not as small as possible. When the deviation is too small, the Kuhn–Tucker condition may not be met during the calculation process, causing convergence problems. Furthermore, too large will lead to the constraints between the systems being out of control. In [37], Li proposed a dynamic relaxation factor method for the CO algorithm multidisciplinary model based on the fixed relaxation factor, so that the size of the relaxation factor can be self-adjusted according to the inconsistency between disciplines, avoiding excessively large or small case. This mathematical method has brought significant contributions to multidisciplinary optimization, as shown in Equation (3).
ε = λ × Δ 2 Δ = | | x c h 1 x c h 2 | | 2   ,   0.5 < λ < 1
where ∆ represents the inconsistency between sublevel 1 and sublevel 2. x c h 1 ,   x c h 2 are the design variables optimization solutions from the two subsystems in the optimization process. λ is a constant coefficient, which is used for the system-level optimization problem that determines a feasible region.
In this method, the relaxation factor is adaptively adjusted according to the inconsistent information of the linking variables in each iteration. However, its applicability is limited and only applies when the linking variables are completely consistent.
To improve the consistency of the ATC method, by referring to the method of Equation (3), the tolerance of dynamic linking variables between i + 1 system level is constructed as the dynamic relaxation factor constraint of ATC method. In this method, the inconsistency of linking variables among subsystems is taken as the constraint condition. The inconsistency between systems is represented by the deviation of the linking variable γ between levels after each iteration, as shown in Equation (4).
γ = n   j , k = 1 | | y i + 1 j y i + 1 k | | 2 2
In order to ensure that the optimization problem has a reasonable calculation domain, the coefficient ρ is added to γ, 0.5 ρ 1 . In the early stage of the optimization process, when the inconsistency between the linking variables is large, γ plays a major role in improving the efficiency of the global search in the system-level optimization process. However, in the later stage, when the amount of inconsistency between disciplines γ gradually decreases, then a fixed relaxation term φ is added to the relaxation factor to limit the lower limit of the relaxation factor. It is used to ensure that the system-level optimization problem always has a certain size of feasible region during the dynamic change of the relaxation factor, as shown in Equation (5).
ε y * = ρ 2 γ + φ
where φ 0.01 ,   0.00001 , the value of φ depends on the order of magnitude of the variable. For example, when the magnitude of the linking variable is about 10 1 and 10 0 , we can choose 0.001 as its value, as used in the mathematical example of Equation (9) in this paper. Of course, it is also feasible to choose a smaller value. However, if the value is too small, non-convergence may occur, and appropriate expansion is needed at this time.
The dynamic relaxation factor ATC model is constructed by introducing ε y * in the original ATC model. In Equation (6), the dynamic relaxation factor ATC model consisting of system level and subsystem level is shown.
{ S y s t e m   l e v e l :             m i n ω i j R R i j R i j U + ω i j y y i j y i j U + ε R + ε y *   r e s p e c t   t o   x ˜ i j ,   y i j ,   y i + 1 j ,   R i + 1 j ,   ε R ,   ε y * w h e r e         R i j = r i j R i + 1 j , x ˜ i j , y i j ε y * = ρ 2 n j , k = 1 | | y i j y i k | | 2 2 + φ s . t .         ω i j R i = 1 n j = 1 c i + 1 k | | R i j R i j L | | 2 2 ε R               ω i j y i = 1 n j = 1 c i + 1 k | | y i j y i j L | | 2 2 ε y * g i j R i j , x i j , y i j 0 h i j R i j , x i j , y i j = 0 x ˜ i j m i n x ˜ i j x ˜ i j   m a x   , y i j m i n y i j y i j m a x S u b s y s t e m   l e v e l     : m i n ω i j R R i + 1 j R i + 1 j U + ω i j y y i + 1 j y i + 1 j U + ε R + ε y *   r e s p e c t   t o   x ˜ i + 1 j ,   y i + 1 j w h e r e         R i + 1 j = r i + 1 j x ˜ i + 1 j , y i + 1 j s . t .       g i j R i + 1 j , x i + 1 j , y i + 1 j 0           h i j R i + 1 j , x i + 1 j , y i + 1 j = 0 x ˜ i + 1 j m i n x ˜ i + 1 j x ˜ i + 1 j   m a x   , y i + 1 j m i n y i + 1 j y i + 1 j m a x
where the parameters φ and ρ are selected according to the actual situation in order to ensure the convergence of the algorithm.   ε R is controlled by the system. In this paper, the default weight coefficient ω is 1, which can be adjusted according to the actual weight for different problems.
The calculation process of the dynamic relaxation factor method is shown in Figure 4:

2.4. Math Problem Application

In order to better prove the ATC method of inter-system linking variable as relaxation factor constraint, a simple mathematical example, Equation (7) [36], is introduced to compare the dynamic relaxation factor and fixed relaxation factor of 0.1, 0.001, and 0.00001, respectively. This example can be decomposed to make it exist only as related variables, which reduces the difficulty of calculation and can be solved by ATC algorithm.
M i n         F = x 1 2 + x 2 2 s . t .     g 1 = x 1 + 0.1 x 2 4 0                   g 2 = 2 0.1 x 1 x 2 0   0 < x 1 < 5 ,   0 < x 2 < 5
The theoretically optimal solution of this problem is x * = 0.19802 , 1.9802 , f x * = 3.9604 . Aiming at this problem, the traditional fixed relaxation factor ATC optimization method is used as Equation (8):
{ S y s t e m   l e v e l :             M i n       F = x 01 2 + x 02 2 + ε 1 + ε 2 W h e r e       R 1 = r 1 = 0.5 × x 01 2 + x 02 2 ,   R 2 = r 2 = 0.5 × x 01 2 + x 02 2 s . t .       x 01 x 11 2 + x 02 x 12 2 ε 1 x 01 x 21 2 + x 02 x 22 2 ε 2 S u b s y s t e m   l e v e l   1                       M i n       j 1 = x 11 x 01 u 2 + x 12 x 02 u 2 s . t .       g 1 = x 11 + 0.1 x 12 4 0 S u b s y s t e m   l e v e l   2                       M i n       j 2 = x 21 x 01 u 2 + x 21 x 02 u 2   s . t .       g 2 = 2 0.1 x 21 x 22 0
The inter-discipline consistency constraint is used as the response deviation to construct the dynamic relaxation factor, as shown in Equation (9).
{ S y s t e m   l e v e l :                 M i n       F = x 01 2 + x 02 2 + ε * W h e r e       R 1 = r 1 = 0.5 × x 01 2 + x 02 2 ,   R 2 = r 2 = 0.5 × x 01 2 + x 02 2 s . t .       x 01 x 11 2 + x 02 x 12 2 ε * x 01 x 21 2 + x 02 x 22 2 ε * ε * = ρ 2 γ + φ = ρ 2 [ x 11 x 21 2 + x 12 x 22 2 ] + φ S u b s y s t e m   l e v e l   1                       M i n       j 1 = x 11 x 01 u 2 + x 12 x 02 u 2 s . t .       g 1 = x 11 + 0.1 x 12 4 0 S u b s y s t e m   l e v e l   2                       M i n       j 2 = x 21 x 01 u 2 + x 21 x 02 u 2   s . t .       g 2 = 2 0.1 x 21 x 22 0
The multidisciplinary optimization platform is built into the Isight software; the above problems are solved respectively. The Modified Method of Feasible Directions (MMFD) algorithm is used in the system level, and the Non-Linear Programming by Quadratic Lagrangian (NLPQL) algorithm is used in the subsystem level. The initial values of the above questions are all set to (1,1). The value of ρ as long as it meets the requirements of the value range and the value of φ have been introduced above.
As shown in Figure 5a–c: (a) is the iterative graph of dynamic relaxation factor constraints, (b,c) are the iterative graphs when the relaxation factor is constrained to 0.1 and 0.001. When the relaxation factor is 0.001, there is no feasible solution due to overconstraint.
Table 1 shows the comparison between traditional ATC and ATC based on dynamic relaxation factor. As shown in Table 1, when the relaxation factor is 0.1 for optimization, the constraint strength is insufficient, which leads to premature convergence, with a large error of 28.79%; when the relaxation factor is 0.001, the optimal value is close to the theoretical solution and the difference is only 2.04%, but there is still a 0.83% gap compared with the dynamic relaxation factor; when the relaxation factor is 0.00001, the relaxation factor constraint is too small, resulting in no feasible solution. In addition, when the dynamic relaxation factor is used for iteration, the number of iterations is only 17 times. Compared with the relaxation factor of 0.1, the number of iteration steps is reduced by three steps, and compared with the relaxation factor of 0.001, it is reduced by 17 steps.
Therefore, dynamic relaxation factor constraints are adopted to avoid constraints and insufficient constraints when selecting relaxation factor constraints, and to improve the efficiency and accuracy of the solution.

3. Engineering Problem Solving

Lightweight optimization for front impact structure of body frame based on active and passive safety is a typical ATC problem. Before the problem is stratified, the optimization objective is lightweight, and the constraints are collision safety and air lift. The research object of this problem is the original truss body, as shown in Figure 6. It is a typical MDO problem to study the lightweight of a truss body while ensuring active and passive safety, which involves many design variables and constraints. Since the linking variables of the two disciplines, crash safety and air lift, are the length and wall thickness of critical pipe fittings, which have a direct impact on lightweighting, the above unstratified problem can be decoupled into a two-level ATC problem.
Therefore, the ATC method is used to decompose the complex engineering system into several subsystems for optimization, which can eliminate the complex system analysis and optimize each subsystem in parallel, reducing the complexity of the whole system. However, auxiliary variables are added in the process of optimization, which also increases the calculation of the whole system to a certain extent, so it is very important to apply efficient optimization algorithm to each subsystem of ATC for optimization calculation. At the same time, the response surface model is used to overcome the problems of spending a lot of time and energy on simulation calculations and low efficiency in finite element simulation analysis.

3.1. Passive Safety

The truss body established in this paper is based on a certain electric formula, combined with the design requirements in the competition rules and created by SolidWorks software. Considering the truss body’s passive safety, materials, processing technology, cost, and other aspects, the pipe fittings of truss body are divided into seven specifications, as shown in Figure 7.
Since (0) the main ring and the front ring are important structures of the truss body, the size specifications and material requirements are relatively strict. They are not used as design variables to ensure the strength and rigidity of the frame structure. The thickness of the pipe fittings (1)–(6) is shown in Table 2 as the design variable t i (i = 1, 2, …, 6).
During the collision, the deformation of the driver’s cockpit should be reduced as far as possible to ensure that the driver has enough living space; the deformation of the front cabin and the motor transmission cabin should be allowed to absorb the kinetic energy transmitted by the body when the collision occurs, so as to reduce the collision acceleration of the cockpit and reduce the impact damage to the driver caused by the collision impact. For passenger cars, a major cause for the injury of the passenger’s leg is that there is body front panel intrusion in the frontal rigid collision. Therefore, the deformation of body front panel which is in contact with the front end of the driver should be reduced as much as possible. This will lead to the driver’s leg injury being small and easy escape [5,42].
Therefore, in the simulation analysis of the frontal crash safety of the racing car, the intrusion at body front panel and the peak acceleration of the cockpit are two important evaluation indicators.

3.2. Active Safety

In the process of vehicle driving, the vertical load of tire is often changed. The cornering characteristic of tire is affected by the change of vertical load, which leads to the change of handling and stability. As shown in Figure 8, there is the influence of different vertical loads on the cornering stiffness of A-70 series tires, B-high performance 70 series tires, and C-60 series tires [12].
Therefore, the vertical load of the air lift tire can be changed by the air lift, which affects the handling and stability of the vehicle to a certain extent. At present, many rail vehicles have been designed for the front part of the vehicle [43,44] to increase the air pressure and reduce the air resistance. In order to reduce the amount of calculation, this paper only designs the head of the truss car body. The front model is established and the aerodynamic characteristics of the front part are optimized to increase the negative pressure during the driving of the front part.
The vertical component of the vehicle becomes the aerodynamic lift during driving. The aerodynamic lift is proportional to the square of the incoming flow velocity, the windward area of the racing car, and the air lift coefficient C l . The air lift coefficient expression is shown in Equation (10):
C l = F z 2 ρ ν 2 S
where S, ν , ρ are defined as windward area, incoming flow velocity, and air density.
The head of the truss body is shown in Figure 9. Pipe fitting 1, pipe fitting 2, and pipe fitting 3 are the longitudinal pipe of body front panel, the middle pipe under the front diaphragm, and the transverse pipe of body front panel. Considering the influence of the body structure on active safety, the length change Δ l j   ( j = 1 ,   2 ,   3 ) of three pipe fittings is taken as the design variable and the length of other members is changed according to the position relationship.
As shown in Table 3, the negative and positive signs are represented as length decrease and length increase. According to the position of the front pipe fittings, the parametric design of the front shape is shown in Figure 10. According to the Computational Fluid Dynamics (CFD) simulation analysis, the aerodynamic characteristics are greatly affected by the concave and convex of the curve at the top of the vehicle head, and the height change value ∆dy at the midpoint of the longitudinal symmetry plane of the vehicle head is taken as the design variable.
The change range curve at the midpoint of the longitudinal symmetry plane is shown in Figure 11. For the points between the longitudinal midpoint and the two ends, the coordinate increment value is changed according to the linear law. The deformation is mainly the change of the unevenness of the curve. The two end points of the control line need to be fixed. The change from the endpoints at both ends to the midpoint of the control line is increasing. When Δ dy < 0, the curve is concave inward, and when Δ dy > 0, the curve is convex outward. The air lift coefficients at different locations are shown in Figure 12.
To make the incoming flow uniform and stable, the wake flow is fully developed, and the calculation domain is too small to avoid backflow. This paper sets the calculation domain of the cuboid as 4 times the head length of the front part, 8 times the head length of the rear part, 6 times the head height of the upper part, and 6 times the head width of the left and right parts, as shown in Figure 13.
ICEM software is used to treat the geometric surface, and tetrahedral mesh elements are selected. Encrypt the grid near the front of the racing car to ensure calculation accuracy [2]. The basic parameter settings are shown in Table 4:

3.3. Engineering Problem ATC Model

To sum up, in order to meet the lightweight needs of the automotive industry, as shown in Equation (11):
Find :   t i , Δ l j ,   Δ dy Objective :   M i n   M t i , Δ l j s . t .         G A t i , Δ l j = A 60 0 G D t i , Δ l j = D 71.6 0 G C l Δ l j ,   Δ dy = m i n C l Δ l j ,   Δ dy t i m i n t i t i m a x Δ l j m i n Δ l j Δ l j m a x Δ dy m i n Δ dy Δ dy m a x i = 1 , 2 ,   ,   6 ; j = 1 , 2 , 3
where A is the peak acceleration during the collision, D is the intrusion at body front panel during the collision, M is the weight of the truss body, C l is the air lift coefficient, t i is the pipe thickness, l i is the pipe length, ∆dy is the change value of midpoint height of longitudinal symmetry plane. According to the calculation framework of the ATC method, the active and passive safety car body crashworthiness structure optimization problem is decomposed into a system-level optimization (lightweight system) and three subsystems (intrusion subsystem, peak acceleration subsystem, and air lift coefficient subsystem) optimization problems. At the same time, the dynamic relaxation factor model is introduced into the ATC method framework. The mathematical description is as follows in Equations (12)–(15):
System level:
{ M i n     P = M t i , Δ l j + ε R C l + ε A + ε D + ε C l ω R C l | | C C l C C l L 2 2 )   ε R C l ω A | | y A y A U 2 2 ε A ω D ( | | y D y D U 2 2 ε D ω C l | | y C l y C l U 2 2 ε C l ε A = ρ A 2 × t A t D 2 2 + Δ l A Δ l D 2 2 + φ A ε D = ρ D 2 × t A t D 2 2 + Δ l A Δ l D 2 2 + φ D ε C l   = ρ C l 2 × Δ l A Δ l C l | | 2 2 + Δ l B Δ l C l | | 2 2 + φ C l
Subsystem level:
Analysis model of peak acceleration at impact:
M i n       J 1 = ω A | | y A y A L | | 2 2 s . t .     G A = A 60 0 y A = Δ l j , t i Δ l i m i n Δ l i Δ l i m a x Δ l j m i n Δ l j Δ l j m a x
Analysis model of intrusion at body front panel:
M i n       J 2 = ω D ( | | y D y D L | | 2 2 ) s . t .     G D = D 71.6 0 y D = Δ l i , t i t i m i n t i t i m a x Δ l j m i n Δ l j Δ l j m a x
Analysis model of air lift coefficient:
M i n       J 3 = ω C l y C l y C l U | | 2 2 + ω R C l C C l C C l L | | 2 2 s . t .     G C l = m i n C C l C C l = R Δ dy , Δ l j         y c l = Δ l j Δ l j m i n Δ l j Δ l j m a x Δ dy m i n Δ d y Δ dy m a x
where y represents the linking variables of different models, and the superscripts U and L represent the system level and the subsystem level. C is the analysis model, which is obtained through the fluid subsystem analysis module. ε D ,   ε A ,   ε C l are defined as the inter-system linking deviation of intrusion, peak acceleration, and air lift coefficient.

3.4. Multidisciplinary Optimization

In order to accurately establish a Response Surface Methodology (RSM), a Design of Experiment (DOE) is usually used to sample the design space. Among the many DOE methods, the Optimal Latin Hypercube Design (Opt LHD) improves the uniformity of the random Latin Hypercube Design, making the fitting of factors and responses more accurate and true. The Opt LHD makes all the test points evenly distributed in the design space as much as possible, with excellent space filling and balance.
Generally speaking, the first-order RSM has a simple structure and requires relatively little test data, but the accuracy of the RSM is low. The combined effect is poor and cannot reflect the interaction between the independent variables. Second-order polynomial RSM is used to make up for this defect as Equation (16) is generally used [41]:
Y ^ = b 0 + i = 1 n b i x i + i = 1 n j = 1 n b i j x i
where b 0 , b i ,   b i j are the undetermined constant term, the first coefficient term, and the quadratic coefficient of the polynomial. Y ^ is the fitting response of the polynomial response surface. Y is the real function response, and x i   ( i = 1 , 2 , n ) is the design variable.
The Opt LHD is used to extract 16 groups of sample points from 10 design variables, t i ,   Δ l j ,   Δ dy . The vehicle weight ( M ), intrusion ( D ), peak acceleration ( A ), and air lift coefficient ( C l ) are taken as the target response. M and   D are solved by LS-DYNA software,   C l is solved by Fluent software, the results are shown in Table 5.
The RSM of the vehicle mass (M), the intrusion (D), the peak acceleration ( A ), and the air lift coefficient ( C l ) are obtained, following Equations (17)–(20):
M = 50.779 0.102 * Δ l 1 0.041 * Δ l 2 + 0.032 * Δ l 3 0.0248 * t 1 0.002 * t 2 0.014 * t 3 0.015 * t 4 + 0.003 * t 5 0.001 * t 6 + 0.001 * Δ l 1 * Δ l 2 + 0.001 * Δ l 1 * Δ l 3 + 0.412 * Δ l 1 * t 1 0.234 * Δ l 1 * t 2 0.04 * Δ l 1 * t 3 + 0.047 * Δ l 1 * t 4 0.286 * Δ l 1 * t 5 + 0.189 * Δ l 1 * t 6 + 0.01 * Δ l 2 * Δ l 2 + 0.001 * Δ l 2 * Δ l 3 + 0.049 * Δ l 2 * t 1 0.106 * Δ l 2 * t 2 + 0.018 *   Δ l 2 * t 3 + 0.076 * Δ l 2 * t 4 0.071 * Δ l 2 * t 5 0.084 * Δ l 2 * t 6 + 0.006 * Δ l 3 * Δ l 3 + 0.33 * Δ l 3 * t 1 0.303 * Δ l 3 * t 2 + 0.057 *   Δ l 3 * t 3 0.096 * Δ l 3 * t 4 0.046 * Δ l 3 * t 5 + 0.06 * Δ l 3 * t 6 0.087 * t 1 * t 1 0.046 * t 1 * t 2 0.122 * t 1 * t 3 0.116 * t 1 * t 4 0.038 * t 1 * t 5 0.038 * t 1 * t 6 0.006 * t 2 * t 2 0.024 * t 2 * t 3 0.033 * t 2 * t 4 + 0.004 * t 2 * t 5 0.003 * t 2 * t 6 0.123 * t 3 * t 3 0.122 * t 3 * t 4 0.011 * t 3 * t 5 0.022 * t 3 * t 6 0.114 * t 4 * t 4 0.012 * t 4 * t 5 0.022 * t 4 * t 6 + 0.009 * t 5 * t 5 + 0.003 * t 5 * t 6 0.001 * t 6 * t 6
D = 24.433 + 0.474 * Δ l 1 + 0.121 * Δ l 2 0.169 * Δ l 3 + 0.075 * t 1 + 0.012 * t 2 + 0.061 * t 3 + 0.044 * t 4 0.017 * t 5 0.001 * t 6 + 0.004 * Δ l 1 * Δ l 1 + 0.003 * Δ l 1 * Δ l 2 0.001 * Δ l 1 * Δ l 3 1.211 * Δ l 1 * t 1 + 1.493 * Δ l 1 * t 2 0.186 * Δ l 1 * t 3 0.175 * Δ l 1 * t 4 + 0.678 * Δ l 1 * t 5 0.56 * Δ l 1 * t 6 + 0.011 * Δ l 2 * Δ l 2 0.005 * Δ l 2 * Δ l 3 0.144 * Δ l 2 * t 1 + 0.444 * Δ l 2 * t 2 0.172 * Δ l 2 * t 3 0.193 * Δ l 2 * t 4 + 0.245 * Δ l 2 * t 5 + 0.272 * Δ l 2 * t 6 1.416 * Δ l 3 * t 1 + 1.09 * Δ l 3 * t 2 + 0.145 * Δ l 3 * t 3 + 0.207 * Δ l 3 * t 4 + 0.157 * Δ l 3 * t 5 0.3 * Δ l 3 * t 6 + 0.272 * t 1 * t 1 + 0.154 * t 1 * t 2 + 0.419 * t 1 * t 3 + 0.349 * t 1 * t 4 + 0.104 * t 1 * t 5 + 0.113 * t 1 * t 6 + 0.043 * t 2 * t 2 + 0.13 * t 2 * t 3 + 0.121 * t 2 * t 4 0.015 * t 2 * t 5 + 0.017 * t 2 * t 6 + 0.547 * t 3 * t 3 + 0.442 *   t 3 * t 4 + 0.035 * t 3 * t 5 + 0.087 * t 3 * t 6 + 0.322 * t 4 * t 4 0.005 * t 4 * t 5 + 0.053 * t 4 * t 6 0.058 * t 5 * t 5 0.027 * t 5 * t 6 0.004 * t 6 * t 6
A = 43.549 + 0.008 * Δ l 1 + 0.012 * Δ l 2 0.009 * Δ l 3 + 0.01 * t 1 + 0.002 * t 2 + 0.001 * t 3 + 0.008 * t 4 0.002 * t 5 0.158 * Δ l 1 * t 1 + 0.049 * Δ l 1 * t 2 + 0.03 * Δ l 1 * t 3 + 0.042 * Δ l 1 * t 4 + 0.045 * Δ l 1 * t 5 0.105 * Δ l 1 * t 6 + 0.006 * Δ l 2 * Δ l 2 + 0.002 * Δ l 2 * Δ l 3 0.055 * Δ l 2 * t 1 + 0.01 * Δ l 2 * t 2 + 0.068 * Δ l 2 * t 3 0.055 * Δ l 2 * t 4 + 0.019 * Δ l 2 * t 5 + 0.025 * Δ l 2 *   t 6 + 0.001 * Δ l 3 * Δ l 3 0.075 * Δ l 3 * t 1 + 0.111 * Δ l 3 * t 2 0.055 * Δ l 3 * t 3 + 0.055 * Δ l 3 * t 4 0.013 * Δ l 3 * t 5 0.025 * Δ l 3 * t 6 + 0.035 * t 1 * t 1 + 0.021 * t 1 * t 2 + 0.042 * t 1 * t 3 *   + 0.051 * t 1 * t 4 + 0.015 * t 1 * t 5 + 0.015 * t 1 * t 6 + 0.008 * t 2 * t 2 + 0.009 * t 2 * t 3 + 0.024 * t 2 * t 4 + 0.001 * t 2 * t 5 + 0.003 * t 2 * t 6 + 0.012 * t 3 * t 3 + 0.04 * t 3 * t 4 0.006 * t 3 * t 5 + 0.001 * t 3 * t 6 + 0.064 * t 4 * t 4 + 0.007 * t 4 * t 5 + 0.012 * t 4 * t 6 0.005 * t 5 * t 5 0.002 * t 5 * t 6
C l = 0.0011 * Δ l 1 + 0.0017 * Δ l 2 0.0015 * Δ l 3 0.0008 * Δ dy + 0.0002 * Δ l 2 * Δ l 2 + 0.0001 * Δ l 3 * Δ l 3 + 0.0002 * Δ d y * Δ d y
In summary, the dynamic relaxation factor constraint ATC method is used to build a multidisciplinary optimization design platform for the engineering problem and optimize it. The Non-dominated Sorting Genetic Algorithm II (NSGA-II) algorithm is used for the system level and the NLPQL algorithm is used for each subsystem. The specific process is shown in Figure 14.
The optimized iterative diagram of the vehicle mass M, the intrusion D, the peak acceleration A , and the air lift coefficient C l by numerical calculation is shown in Figure 15a–d.
After solving by MDO, compare the solved data with the initial results, as shown in Table 6. The reduction effect is represented as “-” in the table.
It can be concluded from the table that before and after optimization, the body weight is reduced by 6.50%, while the peak acceleration is reduced by 12.09%, the intrusion is reduced by 29.18%, and the air lift coefficient is reduced by 12.50%. This confirms the value of the dynamic relaxation factor ATC method in practical application of engineering problems.

4. Conclusions

In this paper, the shortcomings of ATC method such as poor convergence and long iteration time are researched. The novel ATC method based on dynamic relaxation factor is proposed by adding a dynamic relaxation factor that is adjusted by the deviation of the linking variables between the levels to ensure the feasible region of the design space. The optimization framework of this method is as follows: the polynomial response surface model is used to construct the approximate model, the dynamic relaxation factor constraint is applied to the tolerance of the linking variables between levels, and the NSGA- II algorithm and the NLPQL algorithm are used for the system level and subsystem level, respectively, for the complex problems. When the dynamic relaxation factor constraint ATC method is used to solve the problem, its convergence is obviously better than that of the fixed relaxation factor constraint, which provides a method for avoiding the convergence difficulty of ATC method.
In order to ensure the application value of the ATC method in engineering practice, this paper researches the lightweight of the trussed body front impact structure based on active and passive safety. The peak acceleration and the intrusion are used as the response values to evaluate the passive safety and the air lift coefficient is used as the response value to evaluate the active safety. Compared with the traditional lightweight design that only considers passive safety, this paper reasonably associates active safety and passive safety through the body structure. After that, the platform of dynamic relaxation factor constrained ATC method is built, each response model is fitted by RSM, the NSGA-II algorithm and the NLPQL algorithm are used for the system level and subsystem level, respectively. Compared with the initial value, the body weight is reduced by 6.50%, the peak acceleration is reduced by 12.09%, the intrusion is reduced by 29.18%, and the air lift coefficient is reduced by 12.50%. It is useful for studying the lightweight of the body structure while improving the active and passive safety. It is proved that the dynamic relaxation factor ATC method is effective in engineering practice.
However, the accuracy of the response surface approximation model method used in this paper is relatively low. Therefore, higher accuracy approximate model methods, such as reinforcement learning and deep learning, will be used in the future to improve the optimization results.

Author Contributions

Methodology, M.W.; software, M.W.; validation, T.W. and M.W.; formal analysis, X.L.; writing—original draft preparation, M.W.; writing—review and editing, T.W.; visualization, T.W.; supervision, D.Q.; project administration, T.W.; funding acquisition, T.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China—Research on the integrated design method of structure-material-performance of automobile body based on big data (51705468).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Coppola, L.; Marco, B.D.; Niola, V. Impact Attenuator Optimum Design for a FSAE Racing Car by Numerical and Experimental Crash Analysis. Int. J. Automot. Technol. 2020, 21, 1339–1348. [Google Scholar] [CrossRef]
  2. Hetawal, S.; Gophane, M.; Ajay, B.K. Aerodynamic Study of Formula SAE Car. In Proceedings of the 12th Global Congress on Manufacturing and Management, Vellore, India, 8–10 December 2014. [Google Scholar]
  3. Boria, S. Behaviour of an Impact Attenuator for Formula SAE Car under Dynamic Loading. Int. J. Veh. Struct. Sys. 2010, 2, 45–53. [Google Scholar] [CrossRef]
  4. Xiong, F.; Wang, D.; Ma, Z. Lightweight optimization of the front end structure of an automobile body using entropy-based grey relational analysis. Proc. Inst. Mech. Eng. D-J. Aut. 2019, 233, 917–934. [Google Scholar] [CrossRef] [Green Version]
  5. Xie, H.; Wang, P.; Wang, H. Lightweight design for an electric vehicle front cabin. J. Plast Environ. Eng. 2018, 25, 1–7. [Google Scholar]
  6. Mages, M.; Seyffert, M.; Class, U. Analysis of the pre-crash benefit of reversible belt pre- pretensioning in different accident scenarios. In Proceedings of the 22nd International Technical Conference on the Enhanced Safety of Vehicles (ESV) National Highway Traffic Safety Administration, Washington, DC, USA, 20–22 April 2011. [Google Scholar]
  7. Woitsch, G.; Sinz, W. Influences of pre-crash braking induced dummy—Forward displacements on dummy behaviour during EuroNCAP frontal crashtest. Acc. Anal. Prev. 2014, 62, 268–275. [Google Scholar] [CrossRef]
  8. Jeongmin, C.; Kunsoo, H. Active Front Steering for Driver’s Steering Comfort and Vehicle Driving Stability. Int. J. Automot. Technol. 2019, 20, 589–596. [Google Scholar]
  9. Wu, W.; Zou, D.; Ou, J. Adaptive Cruise Control Strategy Design with Optimized Active Braking Control Algorithm. Math. Probl. Eng. 2020, 2020, 1–10. [Google Scholar] [CrossRef]
  10. Narjes, A.; Alireza, K.; Pouria, S. Adaptive yaw stability control by coordination of active steering and braking with an optimized lower-level controller. Int. J. Adapt. Control 2020, 34, 1242–1258. [Google Scholar]
  11. Sujuan, S.; Dong, J.; Chuan, B.O. Multiobjective Optimization of Nonlinear Active Suspension System with Time-Delayed Feedback. Math. Probl. Eng. 2020, 2020, 1–11. [Google Scholar]
  12. Yu, Z. The Theory of Automobile, 5th ed.; China Machine Press: Beijing, China, 2009; pp. 138–139. [Google Scholar]
  13. Roulo, D.; Ptasienski, Z.; Mccumber, B. NASCAR Truck Aerodynamic Analysis and Improvement. In Proceedings of the ASME 2017 International Mechanical Engineering Congress and Exposition, Tampa, FL, USA, 3–9 November 2017. [Google Scholar]
  14. Xiang, Z.; Zhi, J.; Huang, J. A systematic approach for streamlined head form design and evaluation of Chinese high-speed train. Int. J. Rail Transp. 2019, 7, 117–139. [Google Scholar] [CrossRef]
  15. Chen, Z.; Liu, T.; Jiang, Z.; Guo, Z. Comparative analysis of the effect of different nose lengths on train aerodynamic performance under crosswind. J. Fluid Struct. 2018, 78, 69–85. [Google Scholar] [CrossRef]
  16. Li, R.; Xu, P.; Peng, Y. Multi-objective optimization of a high-speed train head based on the FFD method. J. Wend Eng. Ind. Aerod. 2016, 152, 41–49. [Google Scholar] [CrossRef]
  17. Yu, M.; Pan, Z.; Jiang, R. Multi-objective optimization design of the head shape of high-speed train based on approximate model. J. Mech. Eng. 2019, 55, 178–186. [Google Scholar]
  18. Tijssens, M.; Bosma, F.; Kietlinski, K. A Methodology and Tool Chain to Develop Integrated Safety Systems. In Proceedings of the 24th International Technical Conference on the Enhanced Safety of Vehicles (ESV), Gothenburg, Sweden, 8–11 June 2015. [Google Scholar]
  19. Schoeneburg, R.; Breitling, T. Enhancement of Active and Passive Safety by Future PRE-SAFE Systems. In Proceedings of the International Technical Conference on the Enhanced Safety of Vehicles, Washington, DC, USA, 6–9 June 2015. [Google Scholar]
  20. Giesing, J.P.; Barthelemy, J.M. A Summary of Industry MDO Application and Needs; NASA: Washington, DC, USA, 1998.
  21. Dexl, F.; Hauffe, A.; Wolf, K. Multidisciplinary multi-objective design optimization of an active morphing wing section. Struct. Multidiscipl. Optim. 2020, 62, 2423–2440. [Google Scholar] [CrossRef]
  22. Anselma, P.G.; Niutta, C.B.; Mainini, L. Multidisciplinary design optimization for hybrid electric vehicles: Component sizing and multi-fidelity frontal crashworthiness. Struct. Multidiscipl. Optim. 2020, 62, 2149–2166. [Google Scholar] [CrossRef]
  23. Hirschler, T.; Bouclier, R.; Duval, A. Isogeometric sizing and shape optimization of thin structures with a solid-shell approach. Struct. Multidiscipl. Optim. 2019, 59, 767–785. [Google Scholar] [CrossRef] [Green Version]
  24. Martin, I.; Hartwig, L.; Bestle, D. A multi-objective optimization framework for robust axial compressor airfoil design. Struct. Multidiscipl. Optim. 2019, 59, 1935–1947. [Google Scholar] [CrossRef]
  25. Denis, V.; Chirkov, A.S.; Ankudinova, A.E. Skorospelov. Multi-objective shape optimization of a hydraulic turbine runner using efficiency, strength and weight criteria. Struct. Multidiscipl. Optim. 2018, 58, 627–640. [Google Scholar]
  26. Cramer, E.J.; Dennis, J.E.; Frank, P.D. Problem formulation for multidisciplinary optimization. SIAM J. Optim. 1994, 4, 754–776. [Google Scholar] [CrossRef]
  27. Haftka, R.T. Simultaneous analysis and design. AIAA J. 1985, 7, 1099–1103. [Google Scholar] [CrossRef]
  28. Balling, R.J.; Sobieszczanski-Sobieski, J. Optimiza-tion of coupled systems: A critical overview of approaches. AIAA J. 1996, 34, 6–17. [Google Scholar] [CrossRef]
  29. Sobieszczanski-Sobieski, J. Optimization by Decomposition: A Step Form Hierarchic to Non-Hierarchic Systems; NASA: Washington, DC, USA, 1989.
  30. Kroo, J.; Altus, S.; Sobieszczanski-Sobieski, J. Multidisciplinary Optimization Methods for Aircraft Prelimin-Ary Design; NASA: Washington, DC, USA, 1994.
  31. Sobieszczanski-Sobieski, J.; Agte, J.; Robert, S.R. Bi-Level Integrated System Synthesis (Bliss); NASA: Washington, DC, USA, 1998.
  32. Michelena, N.; Kim, H.M.; Papalambros, P.Y. A System Partitioning and Optimization Approach to Target Cascading. In Proceedings of the 12th International Conference on Engineering Design, Munich, Germany, 24–26 August 1999. [Google Scholar]
  33. Yao, W.; Wu, Y.; Chen, F. A concurrent subspace collaborative optimization architecture to structural synthetical optimization design. Struct. Multidiscipl. Optim. 2016, 53, 1197–1207. [Google Scholar] [CrossRef]
  34. Shin, M.K.; Park, G.J. Multidisciplinary design optimization based on independent subspaces. Int. J. Numer. Meth. Engng. 2005, 64, 599–617. [Google Scholar] [CrossRef]
  35. Kim, H.M.; Rideout, D.G.; Papalambros, P.Y. Analytical target cascading in automotive vehicle design. J. Mech. Des. 2003, 125, 481–489. [Google Scholar] [CrossRef] [Green Version]
  36. Braun, R.; Gage, P.; Kroo, I. Implementation and Performance Issues in Collaborative Optimization; American Institute of Aeronautics and Astronautics: Reston, VA, USA, 1996. [Google Scholar]
  37. Li, X.; Li, W.; Liu, C. Geometric analysis of collaborative optimization. Struct. Multidiscipl. Optim. 2008, 35, 301–313. [Google Scholar] [CrossRef]
  38. Precup, R.; David, R.; Roman, R. Slime Mould Algorithm-Based Tuning of Cost-Effective Fuzzy Controllers for Servo Systems. Int. J. Comput. Int. Syst. 2021, 14, 1042–1052. [Google Scholar] [CrossRef]
  39. Zapata, H.; Perozo, N.; Angulo, W. A Hybrid Swarm Algorithm for Collective Construction of 3D Structures. Int. J. Artif. Intell. 2020, 18, 1–18. [Google Scholar]
  40. Roman, R.; Precup, R.; Bojan-Dragos, C. Combined Model-Free Adaptive Control with Fuzzy Component by Virtual Reference. Feedback Tuning for Tower Crane Systems. Procedia Comput. Sci. 2019, 162, 267–274. [Google Scholar] [CrossRef]
  41. Fatma, G.; Islem, L.; Chiraz, H. Optimization of boron removal from water by electrodialysis using response surface methodology. Water Sci. Technol. 2020, 81, 293–300. [Google Scholar]
  42. Srinivas, G.R.; Deb, A.; Sanketh, R. An Enhanced Methodology for Lightweighting a Vehicle Design Considering Front Crashworthiness and Pedestrian Impact Safety Requirements. Adv. Mat. Res. 2017, 173, 623–630. [Google Scholar] [CrossRef]
  43. Xiong, J.; Li, T.; Zhang, J. Shape optimization of high-speed trains under multi running conditions. Sci. China Technol. Sci. 2016, 46, 313–322. [Google Scholar]
  44. Muñoz-Paniagua, J.; García, J.; Crespo, A. Genetically aerodynamic optimization of the nose shape of a high-speed train entering a tunnel. J. Wind Eng. Ind. Aerod. 2014, 130, 48–61. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Collaborative optimization method hierarchy.
Figure 1. Collaborative optimization method hierarchy.
Mathematics 09 00907 g001
Figure 2. ATC level delivery.
Figure 2. ATC level delivery.
Mathematics 09 00907 g002
Figure 3. ATC model data flow level transmission design problem.
Figure 3. ATC model data flow level transmission design problem.
Mathematics 09 00907 g003
Figure 4. ATC dynamic relaxation factor calculation process.
Figure 4. ATC dynamic relaxation factor calculation process.
Mathematics 09 00907 g004
Figure 5. Different relaxation factor iteration; (a) dynamic relaxation factor constraints; (b) the relaxation factor is constrained to be 0.1; (c) the relaxation factor is constrained to be 0.001.
Figure 5. Different relaxation factor iteration; (a) dynamic relaxation factor constraints; (b) the relaxation factor is constrained to be 0.1; (c) the relaxation factor is constrained to be 0.001.
Mathematics 09 00907 g005
Figure 6. Electric formula racing.
Figure 6. Electric formula racing.
Mathematics 09 00907 g006
Figure 7. Grouping of pipe fittings.
Figure 7. Grouping of pipe fittings.
Mathematics 09 00907 g007
Figure 8. Cornering stiffness of different series of tires.
Figure 8. Cornering stiffness of different series of tires.
Mathematics 09 00907 g008
Figure 9. Truss body head pipe fittings.
Figure 9. Truss body head pipe fittings.
Mathematics 09 00907 g009
Figure 10. Shape design of the front part.
Figure 10. Shape design of the front part.
Mathematics 09 00907 g010
Figure 11. The change range of the midpoint of longitudinal symmetry plane of the vehicle head.
Figure 11. The change range of the midpoint of longitudinal symmetry plane of the vehicle head.
Mathematics 09 00907 g011
Figure 12. Air lift coefficients at different locations of the midpoint.
Figure 12. Air lift coefficients at different locations of the midpoint.
Mathematics 09 00907 g012
Figure 13. Setting of computational domain for CFD simulation of vehicle head.
Figure 13. Setting of computational domain for CFD simulation of vehicle head.
Mathematics 09 00907 g013
Figure 14. Construction of multidisciplinary optimization platform for engineering problems.
Figure 14. Construction of multidisciplinary optimization platform for engineering problems.
Mathematics 09 00907 g014
Figure 15. The optimized iterative diagram; (a) the vehicle mass M; (b) the peak acceleration A ; (c) the intrusion D; (d) the air lift coefficient C l .
Figure 15. The optimized iterative diagram; (a) the vehicle mass M; (b) the peak acceleration A ; (c) the intrusion D; (d) the air lift coefficient C l .
Mathematics 09 00907 g015
Table 1. Comparison between traditional ATC and dynamic relaxation factor ATC.
Table 1. Comparison between traditional ATC and dynamic relaxation factor ATC.
Solve under Different ConditionsOptimal SolutionError AccuracyNumber of Iterations
Theoretical solution3.9604————
Solution with a relaxation factor of 0.12.820428.79%20
Solution with a relaxation factor of 0.0013.87962.04%34
Solution with a relaxation factor of 0.00001No feasible solution————
Dynamic relaxation factor3.91251.21%17
Table 2. Main parameters of racing rods.
Table 2. Main parameters of racing rods.
Main Components (Design Variables)Initial Value/mmDesign Upper Limit/mmDesign Lower Limit/mmMaterial
(0) Front ring and main ring1.6--4130 steel
(1) Front baffle, roll cage diagonal brace, main ring diagonal brace, side anti-collision structure, battery protection structure1.62.01.24130 steel
(2) The upper and lower supports of the front clapboard and the upper and lower longitudinal pipes of the motor drive cabin3.54.53.06061 aluminum
(3) Truss support pipe for front cabin and motor drive cabin3.54.53.06061 aluminum
(4) Front partition support and main ring support1.51.61.54130 steel
(5) Cockpit and battery protection structure truss support pipe1.62.01.24130 steel
(6) Mounting cross bar for motor, differential, and drive1.21.60.94130 steel
Table 3. Selection of design variable parameters for front part.
Table 3. Selection of design variable parameters for front part.
Design VariableInitial
Value/mm
Lower
Limit/mm
Upper
Limit/mm
Change value of length of tube 10−5.030.0
Change value of length of tube 20−50.075.0
Change value of length of tube 30−25.015.0
Change value of longitudinal symmetry plane of the vehicle head0−20.025.0
Table 4. Basic parameters of CFD simulation.
Table 4. Basic parameters of CFD simulation.
Calculation Domain Boundary TypeParameter Setting
VELOCITY INLET ν = 11.1 m/s
PRESSURE OUTLETFree pressure outlet, hydraulic diameter: 2.04 m, turbulence intensity: 5%
GROUND Slip   wall ,   ν = 11.1   m/s
WALL
FRONT SURFACE Immovable avoidance surface
Table 5. Sampling data and simulation response value.
Table 5. Sampling data and simulation response value.
Design Variable12345678
  t 1 / mm 1.681.651.812.001.841.921.891.87
  t 2 / mm 3.974.034.303.833.774.173.573.90
  t 3 / mm 4.233.504.304.434.373.773.903.83
  t 4 /mm1.511.501.601.541.581.531.511.53
  t 5 / mm 1.951.761.921.871.602.001.711.63
  t 6 / mm 1.951.811.791.891.731.711.682.00
Δ l 1 / mm −14.00−6.67−21.334.3326.33−17.6730.00−3.00
Δ l 2 / mm 0.00−33.3341.6716.6725.0058.33−8.3375.00
Δ l 3 / mm 15.00−3.67−9.00−25.007.00−19.67−14.339.67
Δ dy / mm −2.00−17.00−14.00−11.0016.0013.001.00−5.00
M / kg 46.43 45.77 45.85 47.79 49.22 48.81 47.52 43.84
D / mm 50.19 49.57 50.13 62.90 50.16 62.26 64.71 63.83
A / g       48.42 45.75 49.84 50.93 48.25 50.00 47.76 48.94
C l       0.26 0.30 0.25 0.25 0.23 0.23 0.24 0.21
Design Variable910111213141516
  t 1 / mm 1.711.951.631.791.731.971.761.60
  t 2 / mm 4.103.704.434.233.504.504.373.63
  t 3 / mm 3.703.634.174.504.103.973.574.03
  t 4 /mm1.571.591.551.521.561.551.571.59
  t 5 / mm 1.731.841.651.811.971.681.891.79
  t 6 / mm 1.601.841.871.631.651.761.971.92
Δ l 1 / mm 15.338.0011.6719.00−10.33−25.0022.670.67
Δ l 2 / mm 66.67−41.678.3333.3350.00−16.67−25.00−50.00
Δ l 3 / mm 11.671.67−17.00−1.004.33−6.3312.33−22.33
Δ dy / mm 4.0019.0025.00−20.0022.0010.00−8.007.00
M / kg 46.66 39.19 43.04 43.08 45.20 41.43 45.85 51.78
D / mm 65.43 67.00 70.84 66.51 69.77 78.54 71.04 65.10
A / g       46.86 48.49 47.23 49.10 47.85 49.05 48.71 47.59
C l       0.20 0.29 0.24 0.24 0.23 0.28 0.27 0.32
Table 6. Data comparison before and after optimization.
Table 6. Data comparison before and after optimization.
Response ValueInitial ValueOptimized ValueImprovement Effect
M/kg46.7643.71−6.50%
A/g53.4246.96−12.09%
D/mm70.6250.02−29.18%
C l 0.240.21−12.50%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, T.; Wang, M.; Li, X.; Qin, D. Multidisciplinary Lightweight Optimization for Front Impact Structure of Body Frame Based on Active and Passive Safety. Mathematics 2021, 9, 907. https://doi.org/10.3390/math9080907

AMA Style

Wang T, Wang M, Li X, Qin D. Multidisciplinary Lightweight Optimization for Front Impact Structure of Body Frame Based on Active and Passive Safety. Mathematics. 2021; 9(8):907. https://doi.org/10.3390/math9080907

Chicago/Turabian Style

Wang, Tingting, Mengjian Wang, Xia Li, and Dongchen Qin. 2021. "Multidisciplinary Lightweight Optimization for Front Impact Structure of Body Frame Based on Active and Passive Safety" Mathematics 9, no. 8: 907. https://doi.org/10.3390/math9080907

APA Style

Wang, T., Wang, M., Li, X., & Qin, D. (2021). Multidisciplinary Lightweight Optimization for Front Impact Structure of Body Frame Based on Active and Passive Safety. Mathematics, 9(8), 907. https://doi.org/10.3390/math9080907

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop