Air flow simulation-based stake car wind resistance dynamic optimization method and system

Abstract

The invention discloses a stake car wind resistance dynamic optimization method and system based on airflow simulation, and belongs to the field of automobile pneumatic optimization. A stake car three-dimensional structure model is constructed, and spatial distribution parameters and surface roughness parameters of a stake component are extracted; establishing a multi-working-condition external flow field simulation model, and obtaining a unit area wind resistance coefficient distribution diagram; constructing a wind resistance sensitivity matrix, and identifying a high-sensitivity component; performing disturbance modeling on the high-sensitivity component, generating a local deformation configuration sample and establishing a simulation model set; parallel simulation calculation is carried out, and wind resistance coefficient and aerodynamic stress data are extracted; based on a multi-objective optimization algorithm, determining an optimal disturbance configuration; outputting a structure adjustment suggestion for guiding dynamic design optimization; according to the method, collaborative optimization of wind resistance reduction and structure safety can be achieved, and the pneumatic performance and the energy-saving level of the whole stake car are improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of automobile aerodynamic optimization, in particular to a method and system for dynamic optimization of wind resistance of a warehouse van based on airflow simulation. BACKGROUND

[0002] A warehouse van is a widely used freight vehicle in the logistics transportation process, which is characterized by a fence type fence on the side, and is suitable for transporting agricultural and sideline products, livestock, building materials and other goods that need to be ventilated or irregularly stacked. However, due to the non-closed structure of the warehouse van fence, the airflow penetrates irregularly during driving, the turbulence is severe, the overall wind resistance coefficient of the vehicle is higher than that of a van, the fuel consumption is high, the driving stability is poor, and the nonlinear aerodynamic load is generated on the local structure of the vehicle.

[0003] Currently, the research on wind resistance optimization of warehouse vans mainly focuses on the average wind resistance optimization of fixed overall shape, which is difficult to cope with the dynamic changes of wind resistance under different airflow states during high-speed driving, especially in mountainous areas or crosswind sections. The fence structure in the local area is prone to form turbulent wake vortex, which may cause vehicle yaw, sudden increase in fuel consumption, local structure fatigue and even vehicle compartment instability.

[0004] In addition, the existing technology generally lacks a wind resistance dynamic evaluation mechanism combining real-time airflow field simulation and vehicle structure parameter linkage, which cannot identify and optimize the weak position of the aerodynamic performance in the warehouse fence structure during the design stage, resulting in high design redundancy and difficulty in reducing wind resistance without sacrificing the cargo carrying capacity. SUMMARY

[0005] The purpose of the present application is to provide a method and system for dynamic optimization of wind resistance of a warehouse van based on airflow simulation to solve the problems in the background art.

[0006] In order to achieve the above purpose, the present application provides the following technical scheme: a method for dynamic optimization of wind resistance of a warehouse van based on airflow simulation, comprising: constructing a three-dimensional structure model of the target warehouse van, and obtaining the spatial distribution parameters and surface roughness parameters of the fence components; establishing an external flow field simulation model corresponding to the vehicle, simulating the airflow flow state under multiple typical driving conditions, and obtaining the unit area wind resistance coefficient distribution map acting on each fence component; constructing a fence component wind resistance sensitivity matrix according to the obtained unit area wind resistance coefficient distribution map, and extracting the component number set of the high sensitivity area of wind resistance; modeling the structure disturbance of each component in the component number set, generating a plurality of local deformation configuration samples, and constructing a set of disturbed aerodynamic simulation models; Parallelly simulating the aerodynamic simulation model set to obtain total wind resistance coefficients and local aerodynamic stress distributions of each perturbation sample under each working condition, and form an optimization evaluation function group; Performing multi-objective optimization based on the optimization evaluation function group to determine an optimal perturbation configuration, wherein a component number set corresponding to the optimal perturbation configuration satisfies that a total wind resistance coefficient decreases by more than a preset threshold, and local aerodynamic stress is within an allowable range; According to the optimal perturbation configuration, outputting a structure adjustment suggestion for guiding dynamic design optimization of the fence components of the fence vehicle.

[0007] Preferably, the three-dimensional structure model of the target fence vehicle comprises: Based on laser point cloud scanning of the fence vehicle, high-precision point cloud data of the external structure of the vehicle body is obtained; The point cloud data is denoised, registered and surface fitted by using a point cloud reconstruction algorithm to generate a three-dimensional digital geometric model containing detailed components; The fence components in the three-dimensional model are numbered by region, and the spatial position, length, height and thickness information of each component are extracted to form a spatial distribution parameter matrix; Through surface texture image acquisition and gray texture distribution analysis, combined with a standard roughness modeling library, the surface roughness grade of each fence component is estimated and assigned.

[0008] Preferably, the establishment of the corresponding vehicle external flow field simulation model comprises: Based on the three-dimensional digital geometric model, a computational fluid dynamics simulation domain containing the vehicle model and the external air domain is constructed, and boundary conditions and inlet flow velocity parameters are set; Selecting multiple typical driving conditions, including high-speed straight driving, side wind interference driving and long downhill sliding operation scenarios, and respectively constructing a simulation model group corresponding to each working condition; Using large eddy simulation or Reynolds time-averaged equation solving method to perform steady or transient simulation on each group of working conditions to obtain the wind resistance distribution per unit area of the fence component surface; Projecting the simulation output results onto the surface grid of each component to calculate the ratio of wind resistance to area, and generating a unit area wind resistance coefficient distribution map.

[0009] Preferably, the construction of the fence component wind resistance sensitivity matrix comprises: The unit area wind resistance coefficient distribution map is clustered and integrated according to the component number to calculate the average wind resistance coefficient value of each fence component, and a wind resistance coefficient initial vector is constructed; Taking the component size factor and the surface roughness grade as correction weights, the wind resistance coefficient initial vector is normalized to generate a standardized wind resistance response vector; Based on multi-condition wind resistance response data, the variance and maximum gradient value of wind resistance change of each fence component under different conditions are calculated as wind resistance sensitivity index, and a wind resistance sensitivity matrix is ​​constructed. Set a preset sensitivity threshold, and filter the component numbers in the sensitivity matrix whose corresponding values ​​are greater than the sensitivity threshold to form a set of component numbers in the high wind resistance sensitive area.

[0010] Preferably, the structural perturbation modeling for each component in the component number set includes: Based on the spatial geometric information and boundary constraints of each component in the component number set, a set of disturbance parameters is set, including disturbance amplitude, disturbance direction and disturbance type; Different combinations of perturbation parameters are applied to each component, and multiple local deformation configuration samples are generated by shape parameter transformation. The deformation types include deflection, torsion and thickness deformation. Each local deformation configuration is embedded into the original three-dimensional model of the vehicle to construct the disturbed geometric model of the vehicle, while keeping the structure of the remaining components unchanged; For all perturbation configuration models, a set of corresponding computational fluid dynamics simulation models is generated in batches, and the simulation boundary conditions are uniformly set as the data input set for solving the optimization evaluation function.

[0011] Preferably, parallel simulation calculations of the aerodynamic simulation model set include: For each disturbance configuration model, a unified set of multi-condition simulation parameters is set, including flow velocity, inflow angle, ground sliding conditions and turbulence model, and a condition configuration file is constructed. The perturbation configuration model set is combined with the corresponding operating condition configuration in batches, and parallel simulation is performed using a high-performance computing platform or cloud computing environment. After the simulation, the total drag coefficient and aerodynamic load distribution on the surface of the fence components for each disturbance configuration under each working condition are extracted and stored in vector form. Based on the changes in drag coefficient and local aerodynamic stress of components under different working conditions, a drag optimization objective function and a stress constraint function are constructed to form a multi-objective optimization evaluation function set.

[0012] Preferably, the multi-objective optimization based on the set of optimization evaluation functions includes: Based on the sample set of disturbance configurations and the corresponding wind resistance objective function and stress constraint function, the design variable space and objective function space are constructed, and the total wind resistance coefficient reduction threshold T1 and the maximum allowable stress threshold T2 are defined. A multi-objective evolutionary algorithm is used to iteratively optimize and search the perturbation configuration samples, resulting in a Pareto optimal solution set that balances wind resistance performance and aerodynamic strength. In the Pareto optimal solution set, a solution set is selected where the wind resistance decreases by more than the threshold T1 and the local stress of all disturbing components does not exceed the threshold T2. The set of component numbers involved in the optimal disturbance configuration that meets the screening criteria is denoted as M′, and the wind resistance optimization adjustment suggestions and component deformation parameters corresponding to the configuration are output.

[0013] Preferably, the output structure adjustment suggestions based on the optimal perturbation configuration include: The perturbation parameters of each component in the optimal perturbation configuration are analyzed, including perturbation amplitude, perturbation direction and perturbation type, and the difference analysis is performed with the original component geometric parameters. Based on the difference analysis results, a structural adjustment mapping table is constructed to clarify the deformation operation type and numerical range of each disturbed component; Based on the spatial positioning relationship of components in the whole vehicle, evaluate the geometric compatibility of the disturbance configuration with adjacent components and the overall vehicle shape. If the assembly constraints are not met, propose adjustment suggestions.

[0014] This invention also provides a dynamic optimization system for the wind resistance of a cargo box truck based on airflow simulation, comprising: The 3D structure modeling module constructs a 3D structural model of the target cargo truck and obtains the spatial distribution parameters and surface roughness parameters of its fence components. The working condition simulation module establishes a simulation model of the external flow field of the corresponding vehicle, simulates the airflow state under various typical driving conditions, and obtains the distribution map of the wind resistance coefficient per unit area acting on each fence component. The component identification module constructs a wind resistance sensitivity matrix for fence components based on the obtained wind resistance coefficient distribution map per unit area, and extracts the set of component numbers for areas with high wind resistance sensitivity. The configuration sample generation module performs structural perturbation modeling on each component in the component number set, generates multiple local deformation configuration samples, and constructs a set of perturbed aerodynamic simulation models. The data extraction module performs parallel simulation calculations on the aerodynamic simulation model set to obtain the total drag coefficient and local aerodynamic stress distribution of each disturbance sample under each working condition, and forms an optimization evaluation function set. The configuration selection module performs multi-objective optimization based on the optimization evaluation function set to determine the optimal disturbance configuration. The component number set corresponding to the optimal disturbance configuration satisfies that the total drag coefficient reduction is greater than the preset threshold and the local aerodynamic stress is within the allowable range. The structural adjustment suggestion output module outputs structural adjustment suggestions based on the optimal disturbance configuration to guide the dynamic design optimization of the fence components of the cargo box vehicle.

[0015] The technical effects and advantages provided by the present invention in the above technical solution are as follows: 1. This invention proposes a systematic, quantifiable, and iterative method for optimizing wind resistance by constructing a three-dimensional structural model of a staked vehicle and a multi-condition aerodynamic simulation environment, combined with wind resistance sensitive area identification, disturbance modeling, and parallel computing. This method can not only accurately identify key component areas affecting the overall vehicle's wind resistance, but also effectively select the optimal disturbance configuration that balances wind resistance reduction and structural safety through a multi-objective optimization algorithm, thereby significantly improving the aerodynamic performance and energy efficiency of the staked vehicle under different driving conditions.

[0016] 2. The dynamic optimization process of this invention supports structural-level perturbation modeling and sensitivity analysis, resulting in higher simulation accuracy and design guidance value. Through modular modeling, perturbation sample generation, parallel simulation calculation, and parameterized output, this invention can achieve an automated closed loop for wind resistance optimization and structural adjustment suggestions, providing an efficient, practical, and scalable engineering solution for lightweight design, energy saving, and wind farm safety of cargo vehicles. Attached Figure Description

[0017] To more clearly illustrate the technical solutions in the embodiments of this application or the prior art, the drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments recorded in this invention. For those skilled in the art, other drawings can be obtained based on these drawings.

[0018] Figure 1 This is a flowchart of the method of the present invention.

[0019] Figure 2 This is a flowchart of the system modules of the present invention. Detailed Implementation

[0020] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0021] Example 1, please refer to Figure 1 As shown in this embodiment, the dynamic optimization method for wind resistance of a cargo box vehicle based on airflow simulation includes: Construct a three-dimensional structural model of the target cargo truck and obtain the spatial distribution parameters and surface roughness parameters of its fence components; Establish a simulation model of the external flow field of the corresponding vehicle to simulate the airflow state under various typical driving conditions and obtain the distribution map of the wind resistance coefficient per unit area acting on each fence component. Based on the obtained distribution map of wind resistance coefficient per unit area, a wind resistance sensitivity matrix of fence components is constructed, and the set of component numbers in the high wind resistance sensitive area is extracted. Structural perturbation modeling is performed on each component in the component number set to generate multiple local deformation configuration samples, and a set of perturbed aerodynamic simulation models is constructed. Parallel simulation calculations are performed on the aerodynamic simulation model set to obtain the total drag coefficient and local aerodynamic stress distribution of each disturbance sample under each working condition, forming an optimization evaluation function set. Multi-objective optimization is performed based on the set of optimization evaluation functions to determine the optimal disturbance configuration. The set of component numbers corresponding to the optimal disturbance configuration satisfies that the decrease in total drag coefficient is greater than a preset threshold and the local aerodynamic stress is within the allowable range. Based on the optimal perturbation configuration, structural adjustment suggestions are output to guide the dynamic design optimization of the fence components of the cargo box vehicle.

[0022] Select a physical prototype of the target cargo van model, ensuring the vehicle surface is clean and unobstructed, and place it in an unobstructed open space or laboratory scanning platform.

[0023] A multi-station laser scanning device (such as a FARO Focus S350 or Leica RTC360 3D laser scanner) is used to set up scanning points around the cargo truck at 360 degrees to complete a multi-view scan of the vehicle body. The scanning device should be set to high-precision mode, with a point spacing accuracy better than 1 mm, and the scanning coverage should include at least the entire vehicle shell, outline, and all fence components.

[0024] The raw point cloud data obtained from each scanning perspective is imported into point cloud processing software (such as Geomagic, CloudCompare, or PolyWorks) for initial stitching and merging.

[0025] The merged point cloud data is subjected to noise point removal, outlier repair and boundary identification. Specifically, the RANSAC (Random Sample Consensus) algorithm can be used to identify noise points and the data can be smoothed by median filtering.

[0026] Global alignment is achieved using point cloud registration technology to ensure consistency between data from different perspectives. This embodiment employs the Iterative Closest Point Registration (ICP) algorithm, whose core idea is to achieve rigid body transformation alignment between point clouds by minimizing the Euclidean distance error function between two sets of point clouds.

[0027] The processed point cloud data is converted into a surface mesh model based on a fitting algorithm. The Poisson surface reconstruction algorithm is preferably used to generate continuous surfaces, forming a 3D digital geometric model containing the complete vehicle body outline and structural details. The generated model retains the spatial information of multi-level components such as fences, posts, and base plates, and is output in STL or OBJ standard formats.

[0028] In the generated 3D digital model, fence components are distinguished from the whole vehicle model through semantic segmentation or manual annotation. Each fence component is identified as an independent geometric object and assigned a unique number.

[0029] For each numbered component, its spatial center coordinates (x, y, z) in the global coordinate system, along with its corresponding length, width, and height dimensions, are extracted. This invention organizes the spatial parameters of all components into a spatial distribution parameter matrix P, where each row represents a component, and the column elements are, in order, the center coordinates and dimension values.

[0030] For example, the i-th row of matrix P is represented as: Pi = [xi, yi, zi, li, wi, hi]; where xi, yi, and zi are the geometric center coordinates of the i-th fence component; li, wi, and hi are the length, width, and height of the component, respectively. This parameter matrix provides the basis for subsequent regional correlation analysis of wind drag coefficients.

[0031] High-resolution surface texture acquisition of the target cargo truck can be achieved by using an industrial-grade structured light camera (such as GOMATOS) combined with a light source to capture surface texture images of each fence component at a fixed angle. The image resolution should be better than 100 pixels / cm.

[0032] Image grayscale analysis was performed on the acquired surface texture images to extract characteristic indicators such as the root mean square error of the grayscale distribution, gradient direction distribution density, and local peak change rate. Based on the ISO 25178 roughness standard, a mapping relationship was established between the above image texture features and roughness levels.

[0033] In this embodiment, the roughness level is defined as five levels (R1~R5), where R1 is the mirror smooth level and R5 is the rough processing level. A KNN (K-nearest neighbor) classification model is constructed using the training set to map the texture features corresponding to each component to the roughness level, and this level information is appended to the component attributes.

[0034] For example: if the mean square error of the grayscale distribution is less than 5 and the local peak density is lower than the threshold of 1, it is judged as R1; if the gradient change rate is higher than 15% and the grayscale fluctuation exceeds 12, it is judged as R5.

[0035] Ultimately, each fence component includes the following parameters: spatial location (x, y, z), dimensional information (length, width, height), and surface roughness grade (R1~R5). All component attributes are uniformly stored in the component attribute library, providing accurate input conditions for subsequent airflow simulation load calculation, wind resistance sensitive area extraction, and disturbance modeling.

[0036] Based on the three-dimensional digital geometric model generated in the previous stage, import it into a CFD (Computational Fluid Dynamics) simulation platform (such as ANSYS Fluent, Simcenter STAR-CCM+, or OpenFOAM) to construct a simulation domain that includes the vehicle model and the air region.

[0037] The outer boundary of the air region is a rectangular envelope, with its length (X-axis) being 6 times the vehicle's length, its width (Y-axis) being 4 times the vehicle's width, and its height (Z-axis) being 3 times the vehicle's height. The vehicle model is placed at the front 1 / 3 of the simulation domain, with its bottom in contact with the ground of the simulation domain to simulate actual road contact conditions.

[0038] The boundary conditions are set as follows: The front-end boundary is set as the velocity inlet, and the input wind speed value is set according to the corresponding operating conditions; The rear boundary is set as a pressure outlet, and the pressure is set to atmospheric pressure. The vehicle surface is set to an impenetrable wall boundary condition; The ground is set to a sliding surface to simulate the movement of a vehicle relative to the ground. Air is assumed to be an incompressible fluid with a density of 1.225 kg / m³ and a viscosity of [missing value]. Pa·second.

[0039] To ensure the simulation model has universality and scenario relevance, multiple simulation sub-models for different typical driving conditions need to be constructed, mainly including the following three categories: High-speed straight-line constant speed driving condition: The input wind speed is set to 100 kilometers per hour to simulate the situation of a vehicle driving at a stable speed on a straight section of a highway.

[0040] Crosswind interference condition: Set the prevailing wind direction to a 30-degree angle with the vehicle's direction of travel and the wind speed to 60 kilometers per hour to simulate coastal areas, mountain passes, and other areas susceptible to crosswinds.

[0041] Long downhill coasting condition: The simulation domain is set with a 5-degree tilt angle, the vehicle is kept in a coasting state without power, and the wind speed is set to 80 kilometers per hour to evaluate the impact of changes in the wake behind the vehicle on the load on the fence.

[0042] A separate simulation sub-model is created for each operating condition, keeping the boundary settings consistent with the air domain size, and only adjusting the flow direction angle and velocity parameters.

[0043] To obtain more detailed data on local airflow variations, this embodiment selects Large Eddy Simulation (LES) as the primary solution method to capture transient turbulent behavior between complex components. When computational resources are limited, the Reynolds-Averaged Navier-Stokes (RANS) model can be used for steady-state conditions, with the SST k-ω model being preferred, balancing computational efficiency and accuracy.

[0044] The solution process is as follows: Mesh generation: A structured and refined mesh is used, with the mesh size refined to within 2 mm on the surface of the fence components to ensure accurate simulation of the turbulent boundary layer; Initial conditions: The initial air temperature is set at 20 degrees Celsius, and the initial turbulence intensity is 5%. Time step setting: The LES method uses a time step of 0.001 seconds and a total simulation time of 2 seconds; Convergence criterion: Residuals less than As a convergence criterion; The solver is set to a pressure-based solver, and the SIMPLE algorithm is used for pressure-velocity coupling.

[0045] After the simulation calculations are completed, the force data for each surface element of the fence component is exported, including the aerodynamic force vector per unit area. The wind resistance coefficient per unit area is defined as follows: The wind resistance coefficient per unit area = local wind resistance / component unit area; where local wind resistance is a composite of wall pressure and shear force output by the simulation platform, and its direction of action is consistent with the prevailing wind direction. The surface of each fence component is divided into multiple grid units, and after calculating the wind resistance coefficient of each unit, a wind resistance coefficient distribution map of the component is formed.

[0046] Taking component number Mi as an example, its surface is composed of n grid cells, and the wind resistance coefficient corresponding to cell j is denoted as Cij. Then the wind resistance coefficient spectrum per unit area of ​​component Mi is the set {Ci1, Ci2, ..., Cin}.

[0047] This atlas can be viewed as a two-dimensional matrix data structure and supports mapping onto a three-dimensional model surface for visualization analysis. Simulation results are exported as CSV or VTK files for subsequent use in wind resistance sensitivity matrix construction and disturbance optimization algorithms.

[0048] The wind drag coefficient distribution map per unit area was clustered and integrated according to the component number. The wind drag coefficient distribution map is the result output by the CFD simulation platform under various typical working conditions, which is represented as the local wind drag (per unit area) value on the surface grid cell of each fence component, in Newtons per square meter.

[0049] The average drag coefficient of a component is obtained by averaging the drag coefficient values ​​of all grid cells under the same component number. Let the component number be Mi, and the drag coefficients of its corresponding n grid cells be C1, C2, ..., Cn. Then the average drag coefficient Ci_avg of component Mi can be expressed as: Ci_avg = (C1 + C2 + ... + Cn) / n; this average value reflects the aerodynamic drag level of the component under the corresponding operating conditions.

[0050] After performing the above processing on all components, an initial wind resistance coefficient vector R0 is constructed, where: R0 = [C1_avg, C2_avg, ..., Cm_avg]; m is the total number of fence components.

[0051] To improve the objectivity and comparability of wind resistance response, this invention introduces two component physical characteristics as weighting coefficients to normalize the initial wind resistance vector R0: Component size factor Si: This represents the influence of component size on wind resistance. Its value is the ratio of the component's projected area to the average projected area, i.e.: Si = Ai / A_avg; where Ai is the projected area of ​​component Mi facing the airflow direction, and A_avg is the average projected area of ​​all components.

[0052] Surface roughness grade factor Qi: In this embodiment, roughness grades R1~R5 are converted into numerical grades 1~5. The larger the value, the rougher the surface and the greater the impact on wind resistance. After normalization, the value of Qi is between [0.2, 1.0].

[0053] The formula for calculating the standardized wind resistance response value Ei of the final component is: Ei = Ci_avg × Si × Qi; Ei is normalized for all components to generate a standard wind resistance response vector R*, whose element range is [0, 1], which is used to construct subsequent matrices with multi-condition data.

[0054] To characterize the stability and sensitivity of a component's wind resistance response under different operating conditions, it is necessary to integrate wind resistance response data from multiple operating conditions and calculate the volatility and response gradient of each component.

[0055] Suppose n typical operating conditions are established (such as high-speed straight road, crosswind, downhill, etc.), and the standardized wind resistance response value of each component under each operating condition is denoted as Ei1, Ei2, ..., Ein. Then the variance of the wind resistance response Vi of the component is defined as: ; where μi is the mean from Ei1 to Ein.

[0056] In addition, the maximum variation gradient Gi is introduced to measure the degree of abrupt change in the component's drag response under different operating conditions. Its calculation formula is as follows: , where j ≠ k, and j, k ∈ [1, n]; that is, take the maximum difference between the response values ​​under any two working conditions.

[0057] Finally, Vi and Gi are linearly weighted and combined to form the sensitivity index Si: ; where α and β are weighting parameters, set according to the simulation scenario, generally taken as α = 0.5, β = 0.5. Alternatively, α = 0.4, β = 0.6 can be set according to the specific scenario to emphasize the impact of abrupt changes. The Si values ​​of all components are combined to form the wind resistance sensitivity matrix S: S = [S1, S2, ..., Sm].

[0058] Each element in the matrix represents the overall response sensitivity of the corresponding component to changes in wind resistance under different operating conditions. The higher the value, the more significant its impact on the aerodynamic performance of the whole vehicle.

[0059] To identify aerodynamically weak areas in a structure, a selected threshold T is set in the sensitivity matrix S. If Si ≥ T for a certain component, then the component is identified as a highly wind-resistant component.

[0060] In the experiment, the present invention sets T to 0.8 (normalized value), that is, only components within the highest sensitivity of 20% are extracted.

[0061] Construct a set of component numbers M, which contains all component numbers that satisfy Si ≥ T: M = {Mi | Si ≥ T, i∈ [1, m]}; set M is the target region for subsequent disturbance modeling and structural optimization, representing the group of components with the strongest wind resistance response and the greatest impact on the aerodynamic performance of the whole vehicle.

[0062] After determining the set M of highly wind-sensitive components, structural disturbance modeling is prepared for each component in the set. First, a set of disturbance parameters D needs to be established. This parameter set describes the control dimensions of the component during structural deformation, and mainly includes: Disturbance amplitude Δ: Represents the intensity of change in the local geometric features of a component, set as a percentage range relative to the original size. For example, the value range of disturbance amplitude Δ can be set to ±5%, ±10%, or ±15%.

[0063] Disturbance direction Θ: Defines the directionality of the disturbance, mainly including: along the prevailing wind direction (X-axis); the normal direction perpendicular to the prevailing wind direction (Y-axis); and local torsional direction (rotation around the Z-axis). Disturbance type τ: Indicates the deformation mode adopted by the structural disturbance, including but not limited to: linear stretching (increase in member length); local contraction (decrease in member height); central depression (depression within the member's central area); edge upturning or downturning (deflection at the member's ends); and torsional deformation (rotation of the cross-section). The parameter set D can be formalized as a triple: D = (Δ, Θ, τ).

[0064] Each component can select multiple combination schemes from the parameter set D according to the preset disturbance strategy, and generate corresponding local deformation configurations respectively.

[0065] For each component Mi, based on the perturbation parameter set D, multiple local perturbation deformation configuration samples are constructed using the shape parameter-driven modeling method. The specific operation steps are as follows: Import the original geometry of the component into a CAD environment (such as Siemens NX, SolidWorks, or Rhinoceros) and build a parametric modeling sketch.

[0066] By using each combination of disturbance parameter set D as input variables, the control points, sketch constraints, or curve fitting relationships are adjusted to achieve geometric deformation.

[0067] The shape-driven modeling module generates multiple sets of different deformation configuration samples. For example, for component Mi, deformation samples Mi1, Mi2, Mi3...Mik are generated, and each sample corresponds to a set of perturbation parameters.

[0068] Ensure that each deformable sample does not overlap or break at geometric continuity and connection surface boundaries to maintain its embedding compatibility in the whole vehicle model.

[0069] A typical example of a perturbation deformation sample is shown below: Component number: M17; Perturbation parameter combination: Δ = +10%, Θ = Y-axis, τ = central depression; Deformation sample number generated: M17_3; Description: The component sinks by 10% in the central region in the vertical wind direction to simulate the buckling state under wind pressure.

[0070] After constructing the deformable configuration sample, it needs to be embedded into the full vehicle 3D structural model. The operation steps are as follows: Load the original 3D model of the vehicle and locate the geometric position, connection points, and boundary relationships of component Mi.

[0071] Delete the original component Mi and replace it with the perturbation sample Mi_j, ensuring interface alignment and keeping adjacent components unchanged.

[0072] Verify the integrity of the model to ensure that the replacement of disturbed components does not affect the overall shape closure of the vehicle body.

[0073] Each disturbance configuration replaces only one component, while the remaining components remain unchanged, in order to analyze the impact of a single component disturbance on the overall wind resistance.

[0074] Finally, for each sensitive component Mi, k disturbance configuration samples are used to generate k vehicle-level 3D models, forming a set of vehicle models after disturbance.

[0075] Import all perturbation configuration models into the CFD simulation platform to automatically generate corresponding simulation calculation models, including: Geometric import and flow field construction: Each disturbance model is imported into the simulation software, maintaining the same boundary dimensions, inlet wind speed, outlet pressure, and other conditions as the original simulation field.

[0076] Mesh generation and refinement: Local mesh refinement is performed in the disturbed component area to ensure simulation accuracy; the original simulation model settings are used in the remaining areas.

[0077] Physics field settings and solution configuration: Use the same solver settings as the original model (such as RANS or LES). Set consistent parameters such as time step and convergence criteria; Maintain consistent wind direction, temperature, and turbulence intensity to ensure comparability of comparison results.

[0078] Model naming and organization: Each disturbance simulation model is named with the component number and disturbance combination number, such as M17_3_FlowModel, forming a structured model library.

[0079] Finally, a multidimensional disturbance simulation sample set is obtained, denoted as: F = {F1, F2, ..., Fk}; where each Fi represents an aerodynamic simulation model corresponding to a specific disturbance configuration.

[0080] First, a unified set of simulation parameters needs to be established for each disturbance configuration. To ensure comparability, this invention standardizes and models multiple typical driving conditions, specifically including: High-speed straight-line driving condition: Simulated vehicle speed of 100 km / h, wind direction is the same as vehicle direction; Crosswind interference condition: Set the wind direction to be at a 30-degree angle to the vehicle direction, and the wind speed to 60 kilometers per hour; Downhill skiing conditions: Set a slope of 5 degrees and a wind speed of 80 kilometers per hour.

[0081] The above operating conditions correspond to different boundary condition settings in the simulation platform, including: Velocity inlet boundary conditions: specify wind speed and flow direction; Sliding ground boundary conditions: Simulating relative driving conditions; Turbulence model selection: This invention preferably uses the SST k-ω model in the RANS (Reynolds-Averaged Navier–Stokes) model, which balances simulation efficiency and accuracy; Other parameters: air density is 1.225 kg / m³, dynamic viscosity is... Pa·second.

[0082] The above parameters are uniformly encapsulated into a working condition configuration file (e.g., JSON or XML format) for subsequent simulation calls and automated deployment.

[0083] After completing the correspondence between configuration samples and operating parameters, the task is uploaded to a high-performance computing platform or cloud computing environment (such as ANSYS HPC cluster, AWS EC2 CFD computing service, etc.), and a parallel simulation scheduling mechanism is adopted to significantly improve computing efficiency.

[0084] The specific steps of parallel scheduling are as follows: Task generation: Each perturbation configuration model is combined with all preset operating conditions to form a perturbation operating condition task pair. For example, the simulation task for configuration Mi_j under operating condition Ck is T(i,j,k).

[0085] Resource allocation strategy: Use task pool scheduling algorithms (such as static block allocation or dynamic load balancing mechanism) to distribute tasks to multiple computing nodes for parallel execution.

[0086] Automated solution process: The solution process for each task is controlled by scripts, including mesh generation, initialization, iteration control, and convergence determination (residual less than 1 / 3). ) and export the results.

[0087] Automatic result archiving: After each simulation task is completed, the output results are automatically compressed and archived to the server database by naming the configuration number + operating condition number.

[0088] This parallel mechanism can complete hundreds of simulation tasks within 48 hours, significantly improving the efficiency of system evaluation.

[0089] After the simulation is completed, key evaluation indicators need to be extracted from each set of perturbation models, mainly including: Total drag coefficient (Cd): Automatically calculated by CFD software, it is defined as the ratio of aerodynamic drag to air dynamic pressure acting on the vehicle, expressed as: Where: Fd is the total wind resistance in Newtons; ρ is the air density; V is the relative wind speed; and A is the windward area. This index is used as the main optimization objective, and the smaller the expected value, the better.

[0090] Local aerodynamic stress distribution: Extract the pressure and shear force per unit area on the surface of the disturbed component, calculate the maximum stress σmax and its location, and use it to determine the structural safety and fatigue risk of the component.

[0091] Data vector format organization: The results of each disturbance configuration under multiple operating conditions are organized into vector groups. For example: The drag vector of configuration Mi_j is: [Cd1, Cd2, Cd3]; The maximum stress vector of configuration Mi_j is: [σmax1, σmax2, σmax3].

[0092] All perturbation configuration data are organized in matrix form and input into the optimization function module.

[0093] Based on the data extracted from the simulation, a set of multi-objective optimization evaluation functions is established, including at least one main objective function and one constraint function: Objective function F1: Minimize the weighted average of the total drag coefficients: F1(Mi_j) = (Cd1 + Cd2 + Cd3) / 3; where Cd1~Cd3 are the total drag coefficients of configuration Mi_j under three operating conditions. The objective is to find the disturbance configuration that minimizes F1.

[0094] Constraint function F2: Maximum stress constraint: F2(Mi_j) = max(σmax1, σmax2, σmax3) ≤ σ_allow; where σ_allow is the safe stress threshold of the material, which is set to the range of 120 MPa to 250 MPa according to the safety standards of commonly used structural components of aluminum alloy or steel, depending on the material of the component.

[0095] The final optimized evaluation function set is expressed as follows: MinimizeF1(Mi_j); Subject toF2(Mi_j) ≤ σ_allow; This function set supports subsequent optimization of perturbation configurations using genetic algorithms (such as NSGA-II) or multi-objective particle swarm optimization (MOPSO).

[0096] We construct a perturbation optimization design space and objective function space, which serve as the input parameter set for a multi-objective evolutionary algorithm.

[0097] The disturbance design variables are denoted as a triple: D = (Δ, Θ, τ); where: Δ is the disturbance amplitude, defined as the percentage change in the geometric dimensions of the component, with a value range of ±5%, ±10%, ±15%; Θ is the disturbance direction, including along the wind direction (X-axis), perpendicular to the wind direction (Y-axis), and torsion around the Z-axis; τ is the disturbance type, including deformation modes such as stretching, compression, indentation, torsion, and edge warping.

[0098] Each set of perturbation variables defines a unique perturbation configuration sample, corresponding to the simulation model and its output data.

[0099] The objective function has a multi-objective structure, including: The objective function for wind resistance, F1, is to minimize the weighted average of the overall drag coefficient Cd of the vehicle under various operating conditions. Let the drag coefficients under n typical operating conditions be Cd1 to Cdn, then: F1 = (Cd1 + Cd2 + … + Cdn) ÷ n; The stress constraint function F2 constrains the maximum aerodynamic stress σ_max of each disturbed component to not exceed the allowable stress σ_allow of the component material. Specifically, it is defined as: F2 = max(σ1, σ2, ..., σm) ≤ σ_allow; where σi represents the maximum aerodynamic stress response of the i-th disturbed component under multiple operating conditions, and m is the number of disturbed components.

[0100] Drag reduction threshold T1: The percentage decrease in drag coefficient relative to the original undisturbed model, preferably set to 8%; Stress tolerance threshold T2: Determined based on the mechanical properties of the material. In this embodiment, 150 MPa is used as σ_allow for the aluminum alloy fence component.

[0101] After defining the variable space and objective function, the Non-Dominated Sorting Genetic Algorithm II (NSGA-II) is used as the optimization solution tool in this invention.

[0102] The core process of the NSGA-II algorithm is as follows: Population initialization: Randomly select perturbation combinations from the perturbation configuration sample set as the initial population. The population size is set to 100 individuals, and each individual corresponds to a perturbation configuration.

[0103] Fitness assessment: For each individual, calculate its value under the objective function F1 and the constraint function F2, and classify it according to the non-dominated ranking principle.

[0104] Selection and crossover operations: Elite individuals are selected by crowding distance sorting to participate in crossover and mutation, generating the next generation of perturbation combinations.

[0105] Iterative optimization: Repeat the evolutionary process of crossover, mutation, and fitness evaluation, with the number of iterations set to 50 rounds, or until the objective function changes converge.

[0106] Output Pareto optimal solution set: All non-dominated solutions constitute the Pareto front, representing the optimal trade-off between wind resistance performance and structural safety.

[0107] This algorithm has global optimization capabilities and is suitable for complex nonlinear multi-objective problems.

[0108] The Pareto optimal solution set obtained from the optimization algorithm is further screened to ensure that the final optimization result is feasible and safe in practical engineering.

[0109] The filtering criteria are as follows: The wind resistance objective function value F1 in the solution set must be greater than the wind resistance reduction threshold T1 (e.g., 8%) relative to the unperturbed model. The maximum aerodynamic stress σ_max of all disturbed components under multiple operating conditions must be less than the allowable stress threshold T2 (e.g., 150 MPa).

[0110] The solutions that satisfy the above two conditions constitute the final set of optimized feasible solutions, from which the configuration with the best performance can be selected for implementation.

[0111] For the finally determined optimal perturbation configuration, the component numbers involved are extracted, an optimized component number set M′ is generated, and structural adjustment suggestions are generated simultaneously.

[0112] The specific steps are as follows: Traverse the structural model of the optimal perturbation configuration, extract the component numbers of all perturbations, such as M17, M22, M38, etc., and form a set M′; For each disturbed component, record its disturbance parameters (Δ, Θ, τ) as a structural deformation suggestion. For example: M17: +10% stretch, Y-axis direction; M22: -5% thickness, Z-axis direction; M38: +15% torsion, 12 degrees angle.

[0113] The above set of numbers and parameter suggestions can be output as an engineering adjustment manual, which can be directly connected to CAD modeling or production line manufacturing instructions.

[0114] For all disturbed components in the final determined optimal disturbance configuration, disturbance parameters are analyzed, and their deformation data, including disturbance amplitude, disturbance direction and disturbance type, are extracted.

[0115] The disturbance parameter triplet is denoted as: D(i) = (Δi, Θi, τi); where: Δi represents the disturbance amplitude of component i, defined as the percentage change relative to the original size, in % %. Θi represents the direction of disturbance, with values ​​including prevailing wind direction (X-axis), vertical wind direction (Y-axis), and vertical or torsional direction (Z-axis); τi represents the type of disturbance, supporting five basic deformations: linear stretching, local contraction, central depression, edge warping, and cross-sectional torsion.

[0116] The original geometric parameters of component i are denoted as: G0(i) = (Li, Hi, Ti, αi); where: Li is the length, Hi is the height, Ti is the thickness, and αi is the original curvature or angle parameter.

[0117] For component i in the optimal configuration, the deformed parameter G1(i) is extracted, and the difference ΔG(i) is calculated: ΔG(i) = G1(i) − G0(i); this difference is used to quantitatively describe the deformation intensity. For example, if ΔG(i).Li = +12 mm, it means that the length of the component increases by 12 mm after the disturbance; if αi changes from 0° to 15°, it means that the component experiences 15 degrees of edge warping. After unifying the units of all difference data, they are recorded in an intermediate data structure for subsequent suggestion generation.

[0118] Based on the above difference analysis results, this invention constructs a structural adjustment mapping table for engineering implementation, which clearly indicates the deformation type, numerical range and operation direction to be implemented for each disturbed component.

[0119] Table 1 Structural Adjustment Mapping Table

[0120] As shown in Table 1, the component number is M17; the disturbance type τ = central depression; the disturbance direction Θ = Y-axis (perpendicular to the wind direction); the original height Hi = 800 mm, and the deformed height H'i = 720 mm; therefore, Δi = −10%.

[0121] Based on the disturbance type and the magnitude of the difference quantization, and combined with the standard processing methods and process specifications in the engineering database, adjustment suggestions are automatically generated. For example: Tensile disturbance → Increase the length of the rod, which can be achieved by welding or module extension; Torsional disturbance → Adjust the angle of the connecting piece or use an eccentric support; Thickness disturbance → Use thicker plates or surface overlay welding.

[0122] Considering that although the disturbance configuration is optimal in terms of wind resistance performance, it may cause assembly interference, local interference or geometric overrun problems in the actual vehicle structure, a space compatibility assessment is required to ensure that the optimal design is feasible for assembly.

[0123] Extract the 3D positioning data of each component in the vehicle model: P(i) = (Xi, Yi, Zi); where Xi, Yi, and Zi represent the center position of component i in the vehicle coordinate system.

[0124] The AABB (Axis-Aligned Bounding Box) algorithm is used to quickly detect collisions between the deformed geometric boundary of the disturbed component and the boundaries of its neighboring components to determine whether there are assembly conflicts or overlaps.

[0125] If a spatial conflict is detected, the system will automatically adjust the disturbance amplitude Δi or direction Θi based on the degree of overlap of the conflict boundaries, generating a "feasible deformation alternative." For example, the original suggestion was +15% torsion, but the conflict area was too large; after fine-tuning, the suggestion is +10% torsion, eliminating the conflict. This adjustment process ensures that design optimization will not disrupt the original assembly logic and structural stability.

[0126] Example 2, please refer to Figure 2 As shown in this embodiment, the dynamic optimization system for wind resistance of a cargo box vehicle based on airflow simulation includes: The 3D structure modeling module constructs a 3D structural model of the target cargo truck and obtains the spatial distribution parameters and surface roughness parameters of its fence components. The working condition simulation module establishes a simulation model of the external flow field of the corresponding vehicle, simulates the airflow state under various typical driving conditions, and obtains the distribution map of the wind resistance coefficient per unit area acting on each fence component. The component identification module constructs a wind resistance sensitivity matrix for fence components based on the obtained wind resistance coefficient distribution map per unit area, and extracts the set of component numbers for areas with high wind resistance sensitivity. The configuration sample generation module performs structural perturbation modeling on each component in the component number set, generates multiple local deformation configuration samples, and constructs a set of perturbed aerodynamic simulation models. The data extraction module performs parallel simulation calculations on the aerodynamic simulation model set to obtain the total drag coefficient and local aerodynamic stress distribution of each disturbance sample under each working condition, and forms an optimization evaluation function set. The configuration selection module performs multi-objective optimization based on the optimization evaluation function set to determine the optimal disturbance configuration. The component number set corresponding to the optimal disturbance configuration satisfies that the total drag coefficient reduction is greater than the preset threshold and the local aerodynamic stress is within the allowable range. The structural adjustment suggestion output module outputs structural adjustment suggestions based on the optimal disturbance configuration to guide the dynamic design optimization of the fence components of the cargo box vehicle.

[0127] The above description is merely a specific embodiment of this application, but the scope of protection of this application is not limited thereto. Any changes or substitutions that can be easily conceived by those skilled in the art within the scope of the technology disclosed in this application should be included within the scope of protection of this application.

Claims

1. A method for dynamic optimization of wind resistance of a gondola car based on airflow simulation, characterized in that: The application relates to a method for dynamic design optimization of a warehouse fence vehicle. The method comprises the following steps: a three-dimensional structural model of the target warehouse fence vehicle is constructed, and the spatial distribution parameters and surface roughness parameters of the fence components are obtained; a corresponding vehicle external flow field simulation model is established, the airflow flow state under multiple typical driving conditions is simulated, and the unit area wind resistance coefficient distribution diagram acting on each fence component is obtained; a fence component wind resistance sensitivity matrix is constructed according to the obtained unit area wind resistance coefficient distribution diagram, and the component number set of the wind resistance high sensitivity area is extracted; structure perturbation modeling is performed on each component in the component number set, a plurality of local deformation configuration samples are generated, and a set of perturbed aerodynamic simulation models is constructed; parallel simulation calculation is performed on the set of aerodynamic simulation models, the total wind resistance coefficient and the local aerodynamic stress distribution of each perturbed sample under each condition are obtained, and a set of optimization evaluation functions is formed; multi-objective optimization is performed based on the set of optimization evaluation functions, and the optimal perturbed configuration is determined, wherein the component number set corresponding to the optimal perturbed configuration satisfies that the total wind resistance coefficient decreases by more than a preset threshold, and the local aerodynamic stress is within an allowable range; 2. The method of claim 1, wherein: structure adjustment suggestions are output according to the optimal perturbed configuration, which are used for guiding the dynamic design optimization of the warehouse fence vehicle fence components. The construction of the three-dimensional structural model of the target warehouse fence vehicle comprises the following steps: high-precision point cloud data of the external structure of the vehicle body is obtained based on laser point cloud scanning of the warehouse fence vehicle; point cloud data is denoised, registered and surface fitted by using a point cloud reconstruction algorithm to generate a three-dimensional digital geometric model containing detailed components; the fence components in the three-dimensional model are numbered by region, the spatial position, length, height and thickness information of each component are extracted, and a spatial distribution parameter matrix is formed; 3. The method of claim 2, wherein: the surface roughness grade of each fence component is estimated and assigned through surface texture image acquisition and gray texture distribution analysis combined with a standard roughness modeling library. The establishment of the corresponding vehicle external flow field simulation model comprises the following steps: a computational fluid dynamics simulation domain containing the vehicle model and the external air domain is constructed based on the three-dimensional digital geometric model, and boundary conditions and inlet flow velocity parameters are set; a plurality of typical driving conditions are selected, including high-speed straight driving, side wind interference driving and long downhill sliding operation scenes, and simulation model groups corresponding to the conditions are constructed respectively; steady-state or transient simulation is performed on each group of conditions by using a large eddy simulation or Reynolds time-averaged equation solving method, and the unit area wind resistance distribution of the fence component surface is obtained; 4. The method of claim 3, wherein: the simulation output results are projected onto the component surface grid, the ratio of wind resistance to area is calculated, and the unit area wind resistance coefficient distribution diagram is generated. The construction of the fence component wind resistance sensitivity matrix comprises the following steps: the unit area wind resistance coefficient distribution diagram is clustered and integrated according to the component number, the average wind resistance coefficient value of each fence component is calculated, and a wind resistance coefficient initial vector is constructed; the component size factor and the surface roughness grade are taken as correction weights, the wind resistance coefficient initial vector is normalized, and a standardized wind resistance response vector is generated; based on the multi-condition wind resistance response data, the variance and the maximum gradient value of the wind resistance change of each fence component under different conditions are calculated as wind resistance sensitivity indexes, and a wind resistance sensitivity matrix is constructed. A preset sensitivity threshold is set, and component numbers corresponding to values greater than the sensitivity threshold in the sensitivity matrix are screened to form a component number set of the wind resistance high sensitivity area.

5. The method of claim 4, wherein: The structure perturbation modeling of each component in the component number set comprises: Based on the spatial geometric information and boundary constraint conditions of each component in the component number set, a set of perturbation parameters is set, including perturbation amplitude, perturbation direction and perturbation type; Different perturbation parameter combinations are applied to each component respectively, and a plurality of local deformation configuration samples are generated by a shape parameter transformation method, and the deformation types include flexure, torsion and thickness deformation; Each local deformation configuration is embedded in the original whole vehicle three-dimensional model to construct a perturbed whole vehicle geometric model, and the structures of the remaining components remain unchanged; For all perturbed configuration models, a corresponding computational fluid dynamics simulation model set is generated in batches, the simulation boundary conditions are uniformly set, and the data input set for solving the optimization evaluation function is obtained.

6. The method of claim 5, wherein: The parallel simulation calculation of the aerodynamic simulation model set comprises: For each perturbed configuration model, a set of uniform multi-working condition simulation parameters is set, including flow velocity, flow angle, ground sliding condition and turbulence model, and a working condition configuration file is constructed; The perturbed configuration model set is combined with the corresponding working condition configuration in batches, and parallel simulation solving is performed by using a high-performance computing platform or a cloud computing environment; After the simulation is completed, the total wind resistance coefficient and the aerodynamic load distribution of the surface of the fence component of each perturbed configuration under each working condition are extracted and stored in vector form; According to the wind resistance coefficient variation and the local aerodynamic stress variation value of the component under different working conditions, a wind resistance optimization objective function and a stress constraint function are constructed to form a multi-objective optimization evaluation function group.

7. The method of claim 6, wherein: The multi-objective optimization based on the optimization evaluation function group comprises: Based on the perturbed configuration sample set and the corresponding wind resistance objective function and stress constraint function, a design variable space and a target function space are constructed, and a total wind resistance coefficient drop threshold T1 and a maximum allowable stress threshold T2 are defined; A multi-objective evolutionary algorithm is used to iteratively optimize and search the perturbed configuration sample to obtain a set of Pareto optimal solutions balanced between wind resistance performance and aerodynamic strength; In the Pareto optimal solution set, a solution set is screened out, in which the wind resistance drop amplitude exceeds the threshold T1 and the local stress of all perturbed components does not exceed the threshold T2; The component number set involved in the optimal perturbed configuration that meets the screening condition is denoted as M', and the wind resistance optimization adjustment suggestion corresponding to the configuration and the component deformation parameter are output.

8. The method of claim 7, wherein: The structure adjustment suggestion output according to the optimal perturbed configuration comprises: The perturbation parameters of each component in the optimal perturbed configuration are analyzed, including perturbation amplitude, perturbation direction and perturbation type, and difference analysis is performed on the original component geometric parameters; Based on the difference analysis result, a structure adjustment mapping table is constructed to clearly define the deformation operation type and numerical range of each perturbed component; Combined with the spatial positioning relationship of the component in the whole vehicle, the geometric compatibility of the perturbed configuration to the adjacent components and the whole vehicle shape is evaluated, and if the assembly constraint is not met, an adjustment suggestion is proposed.

9. A system for dynamic optimization of wind resistance of a gondola car based on airflow simulation, for implementing the method for dynamic optimization of wind resistance of a gondola car based on airflow simulation according to any one of claims 1 to 8, characterized in that: It comprises: A three-dimensional structure modeling module constructs a three-dimensional structure model of the target warehouse fence vehicle and obtains spatial distribution parameters and surface roughness parameters of the fence components of the vehicle; The working condition simulation module establishes an external flow field simulation model corresponding to the vehicle, simulates the airflow flow state under multiple typical driving conditions, and obtains a unit area wind resistance coefficient distribution map acting on each fence component; The component identification module constructs a fence component wind resistance sensitivity matrix according to the obtained unit area wind resistance coefficient distribution map, and extracts a component number set of a wind resistance high sensitivity area; The configuration sample generation module performs structure disturbance modeling on each component in the component number set, generates multiple local deformation configuration samples, and constructs a set of aerodynamic simulation models after disturbance; The data extraction module performs parallel simulation calculation on the set of aerodynamic simulation models, obtains the total wind resistance coefficient and the local aerodynamic stress distribution of each disturbance sample under each working condition, and forms an optimization evaluation function group; The configuration screening module performs multi-objective optimization based on the optimization evaluation function group to determine an optimal disturbance configuration, wherein the component number set corresponding to the optimal disturbance configuration satisfies that the total wind resistance coefficient decreases by more than a preset threshold, and the local aerodynamic stress is within an allowable range; The structure adjustment suggestion output module outputs a structure adjustment suggestion according to the optimal disturbance configuration, which is used to guide the dynamic design optimization of the fence components of the warehouse fence vehicle.

Images