Battery cooling plate design method fusing topological optimization and bionic design

Abstract

The invention discloses a battery cooling plate design method fusing topological optimization and bionic design, which comprises the following steps: firstly, carrying out dimension reduction treatment, setting boundary conditions by taking a cooling plate neutral layer as a topological optimization design domain, and establishing a two-dimensional topological model; a Brinkman model is adopted to describe fluid flow characteristics, a topological optimization model with a multi-objective optimization function and a volume fraction as constraint conditions is constructed, optimization design is carried out on a two-dimensional topological model, and a topological optimization cooling plate is output; analyzing a flow channel of the topological optimization cooling plate; filling the solid area with a feather-shaped structure to obtain the bionic cooling plate; a box dimension method is adopted to calculate flow channels of the two, and if a box dimension difference value is smaller than a set threshold value, a corresponding bionic cooling plate is output; and if the box dimension difference value is not smaller than the set threshold value, returning to the replacement link for correction until the box dimension difference value is smaller than the set threshold value. According to the invention, the structural stability and manufacturability are improved while the heat dissipation performance is ensured.

Description

Technical Field

[0001] This invention relates to the technical field of battery cooling plates, and more particularly to a battery cooling plate design method that integrates topology optimization and biomimetic design. Background Technology

[0002] Energy storage systems are crucial for power supply, but they generate significant heat during charge-discharge cycles, which, if not managed properly, can damage battery performance and lifespan. The core task of battery cooling plates is to manage the heat generated during battery operation, ensuring the battery remains within a safe and efficient temperature range. Existing cooling plate optimizations often focus on optimizing channel dimensions and prioritize maximum heat transfer rate, leading to excessive localized pressure drops. To balance heat transfer rate and hydrodynamic performance, topology optimization has become an effective design method. However, existing topology optimization structures typically contain numerous small and complex channels, failing to adequately consider fabrication feasibility.

[0003] Therefore, the above problems urgently need to be solved. Summary of the Invention

[0004] Purpose of the invention: The purpose of this invention is to provide a battery cooling plate design method that integrates topology optimization and biomimetic design. This invention improves structural stability and manufacturability while ensuring heat dissipation performance.

[0005] Technical Solution: To achieve the above objectives, this invention discloses a battery cooling plate design method integrating topology optimization and biomimetic design, comprising the following steps: S1. First, the three-dimensional flow and heat transfer problem of the battery cooling plate is simplified into a two-dimensional problem. The neutral layer of the cooling plate is used as the topology optimization design domain. Boundary conditions are set, and the corresponding two-dimensional topology model is established. S2. Based on the two-dimensional topology model, the Brinkman model is used to describe the flow characteristics of fluid in porous media. A multi-objective optimization function with the maximum heat transfer rate and minimum flow dissipation as the optimization objectives and a topology optimization model with volume fraction as the constraint is constructed. The two-dimensional topology model is optimized and designed to output a topology-optimized cooling plate that satisfies the maximum heat transfer rate and minimum flow dissipation. S3. Analyze the flow channel distribution of the topology-optimized cooling plate to obtain the main channel and branch channels, and then extract the flow channel contour; arrange the feather-like structure according to the flow channel contour constraint, replace the original solid area to align the feather shaft along the extension direction of the main channel, adjust the arrangement position to ensure the connectivity of the branch channels, and ensure that the volume fraction of the flow channel in the cooling plate remains unchanged, thus obtaining the biomimetic cooling plate. The box-count method is used to calculate the flow channel structure of the topology-optimized cooling plate and the biomimetic cooling plate. If the difference in box-count is less than a set threshold, the corresponding biomimetic cooling plate is output. If the difference in box-count is not less than the set threshold, the process returns to the replacement stage for correction, reducing the number of feather-like structures and adjusting the feather axis direction to ensure that the main flow channel is consistent with the main flow channel direction. The minimum width of the branch channel is constrained to make the width of the branch channel similar to that of the topology-optimized cooling plate until the difference in box-count is less than the set threshold.

[0006] Optionally, in step S1, the overall structure of the cooling plate is symmetrical about the central plane in terms of geometry and boundary conditions. Half of the neutral layer of the cooling plate is selected as the topology optimization design domain, and the lower boundary is set as a symmetrical boundary condition.

[0007] Optionally, step S2 specifically includes the following steps: S21. To achieve an accurate characterization of flow resistance in porous media regions, volume force is introduced as a damping term in the incompressible laminar Navier–Stokes equations, thus obtaining the Brinkman model suitable for topology optimization calculations. S22. Based on Darcy's law, construct an expression relating volumetric force to reverse osmosis rate and flow velocity; S23. Introduce design variables to model the solid and fluid regions in a unified manner. The relationship between reverse osmosis and design variables can be expressed by Darcy interpolation equations. S24. Since there are solid and fluid regions within the design domain, and the heat transfer mechanisms in the two regions are different, design variables are used to unify the heat transfer equations for the solid and fluid regions. S25. In order to obtain a flow channel structure with maximum heat dissipation capacity and low flow power consumption, the heat transfer rate and power dissipation in the design domain are calculated, normalized and combined by weight to construct a multi-objective optimization function. S26. Analyzing the thermal conductivity, fluid flow properties, and design variables in topology optimization, the topology optimization problem is expressed as finding design variables that minimize the multi-objective optimization function under constraints. The constraints are the range of values ​​and volume fraction of the design variables. Thus, the expression for the topology optimization problem is obtained. S27. During the topology optimization solution process, density filtering and projection processing are performed on the design variables to obtain the final design variables, and the topology-optimized cooling plate that satisfies the maximum heat transfer and minimum flow dissipation is output.

[0008] Optionally, in step S21, the Brinkman model is represented by the equation describing the law of conservation of mass and the equation describing the momentum conservation of fluid motion, specifically in the form of: , , in For gradient operators, For fluid density, For fluid velocity, Dynamic viscosity; The pressure of the fluid; For velocity gradient; It is a volume force.

[0009] Optionally, in step S22, the volume force With reverse osmosis and flow rate The relational expression is: .

[0010] Optionally, the Darcy interpolation equation in step S23 is: , , In the formula, For solid permeability, For liquid permeability, The Reynolds number is... For Darcy's number, This is the penalty factor for the Darcy interpolation model. For design variables.

[0011] Optionally, the heat transfer equation in step S24 is: , In the formula, This refers to the isobaric heat capacity of a fluid. For the design domain temperature; The thermal conductivity of the solid region is... The thermal conductivity of the fluid region. Heat generated in the solid region; when When, the heat transfer equation is transformed into the heat conduction control equation for the solid region; when When the heat transfer equation is transformed into the convection-diffusion equation of the fluid region, when When the heat transfer equation is in this state, it represents a smooth transition between the two mechanisms, reflecting the conjugate heat transfer characteristics at the solid-fluid interface.

[0012] Optionally, the expression for the multi-objective optimization function in step S25 is: , , , , , in It is the heat transfer rate within the design domain. It is the power dissipation of the flow within the design domain, for and Normalization process is performed to obtain and ,right and Weighted combination yields a multi-objective optimization function , and These are the minimum and maximum heat transfer rates; and For minimum and maximum flow power dissipation; The coefficient of an ideal heat source; As a weighting factor; For reference temperature, For designing the domain area.

[0013] Optionally, the expression for the topology optimization problem in step S26 is: , In the formula, This represents the volume fraction.

[0014] Optionally, in step S27, a Helmholtz filter is used for density filtering. The expression for the Helmholtz filter is: , In the formula, These are the design variables before filtering; For the filtered design variables; R min The filter radius; Hyperbolic tangent projection is used for projection processing. The expression for hyperbolic tangent projection is: , In the formula, These are the projected design variables; For the projection point; The projection slope.

[0015] Beneficial Effects: Compared with existing technologies, this invention has the following significant advantages: This invention solves the problems of complex structures and poor manufacturability of existing cooling plates by optimizing the flow channel structure and introducing biomimetic design principles, thereby improving structural stability and manufacturability while ensuring heat dissipation performance; This invention significantly improves the overall performance of the cooling plate by designing a cooling plate that balances heat transfer efficiency and flow resistance through topology optimization; This invention simplifies the complex structure obtained from topology optimization by combining biomimetic principles and drawing on the distribution characteristics of bird feathers, not only preserving flow advantages but also improving the manufacturability and structural stability of the cooling plate; This invention uses the box-counting method to determine the inheritance and similarity between the topological structure and the biomimetic structure, continuously revising the biomimetic structure, making the design process scientifically based and universally applicable. Attached Figure Description

[0016] Figure 1 This is a flowchart of the design method of the present invention; Figure 2 This is a simplified schematic diagram of the three-dimensional model of the present invention; Figure 3 This is a schematic diagram of the topology-optimized cooling plate structure of the present invention; Figure 4 This is a schematic diagram of the biomimetic cooling plate design of the present invention; Figure 5 This is a diagram showing the box dimension calculation results of this invention. Detailed Implementation

[0017] The technical solution of the present invention will be further described below with reference to the accompanying drawings.

[0018] like Figure 1 As shown, this invention discloses a battery cooling plate design method that integrates topology optimization and biomimetic design, comprising the following steps: like Figure 2 As shown in Figure S1, the three-dimensional structure of the battery cooling plate is first simplified and analyzed. Considering that the cooling plate is relatively thin and the fluid flow and temperature distribution change weakly along the thickness direction, the three-dimensional flow and heat transfer problem is reduced to a two-dimensional problem. The neutral layer of the cooling plate is used as the topology optimization design domain, and the corresponding two-dimensional topology model is established.

[0019] The overall structure of the cooling plate is geometrically and boundary-conditionally symmetrical about the central plane. Therefore, half of the neutral layer of the cooling plate is selected as the design domain, and the lower boundary is set as a symmetrical boundary condition to ensure the accuracy of calculation while significantly reducing the amount of calculation. This step S1 not only maintains the integrity and symmetry of the physical problem, but also avoids redundant calculations and improves the efficiency and convergence stability of topology optimization.

[0020] S2. Based on the two-dimensional topology model, the Brinkman model is used to describe the flow characteristics of fluid in porous media. A topology optimization model is constructed with the goal of maximizing heat transfer rate and minimizing flow dissipation as the optimization objectives and volume fraction as the constraint condition. The two-dimensional topology model is then optimized.

[0021] Step S2 includes the following specific steps: S21. To accurately characterize the flow resistance in porous media regions, volume force is introduced as a damping term into the incompressible laminar Navier–Stokes equations, thus obtaining the Brinkman model suitable for topology optimization calculations. This yields the equations describing the law of mass conservation and the equations describing the momentum conservation in fluid motion, with the specific equation forms as follows: , , in For gradient operators, For fluid density, For fluid velocity, Dynamic viscosity; The pressure of the fluid; For velocity gradient; It is a volume force.

[0022] S22. According to Darcy's law, volume force With reverse osmosis and flow rate The specific formula is as follows:

[0023] The smaller the reverse osmosis rate α, the less resistance the fluid experiences, representing a fluid; conversely, the larger the reverse osmosis rate α, the greater the force it experiences, representing a solid.

[0024] S23. To achieve unified modeling across different regions, design variables are introduced. To unify the two regions, reverse osmosis rate With design variables The relationship can be expressed by the Darcy interpolation equation, specifically: , , In the formula, For solid permeability, For liquid permeability, The Reynolds number is... For Darcy's number, This is the penalty factor for the Darcy interpolation model.

[0025] S24. Since there are solid and fluid regions within the design domain, and the heat transfer mechanisms in the two regions are different, design variables are required. To unify the heat transfer equations for these two regions, the heat transfer equations are as follows: , In the formula, This refers to the isobaric heat capacity of a fluid. For the design domain temperature; The thermal conductivity of the solid region is... The thermal conductivity of the fluid region. Heat generated in the solid region; when When the heat transfer equation is transformed into the heat conduction control equation for the solid region, heat transfer mainly relies on heat conduction; when When the heat transfer equation is transformed into the convection-diffusion equation of the fluid region, energy transfer is mainly dominated by fluid convection; when When the heat transfer equation is applied, it shows a smooth transition between the two mechanisms and can truly reflect the conjugate heat transfer characteristics at the interface between the solid and the fluid.

[0026] S25. To obtain a flow channel structure with maximum heat dissipation capacity and low flow power consumption, the heat transfer rate and power dissipation within the design domain are calculated, normalized, and combined by weights to construct a multi-objective optimization function. This multi-objective optimization function can obtain a cooling channel topology layout that better meets actual engineering needs while considering both heat dissipation performance and energy loss. The expression of the multi-objective optimization function is: , , , , , in It is the heat transfer rate within the design domain. It is the power dissipation of the flow within the design domain, for and Normalization process is performed to obtain and ,right and Weighted combination yields a multi-objective optimization function , and These are the minimum and maximum heat transfer rates; and For minimum and maximum flow power dissipation; The coefficient of an ideal heat source; As a weighting factor; For reference temperature, For designing the domain area.

[0027] This invention uses a uniform and constant ideal heat source coefficient to represent the actual heating process in the solid region of the battery. To provide an equivalent description; in the mathematical model, introduce This option selects the solid region; when When the material is in a fluid state, it does not generate heat, therefore The integral is guaranteed to be 0; when When , it indicates that the material region is solid and fully participates in the heat dissipation calculation; when At this time, the material region is a material that lies between solid and fluid. This item acts as a weighting factor, used to calculate the proportion of solid material in the material region.

[0028] S26. Analysis of thermal conductivity, fluid flow properties, and design variables in topology optimization; the topology optimization problem is formulated as finding design variables. This makes the multi-objective optimization function Minimum, set the constraint condition as For values ​​between 0 and 1, the volume fraction is 0.5, expressed as: , In the formula, It is the volume fraction. Set to 0.5.

[0029] S27. In the process of topology optimization, it is inevitable that phenomena such as checkerboard patterns, gray-scale cells and mesh dependencies will appear in the design domain. In order to obtain optimization results with clear boundaries, continuous structure and easy subsequent processing, it is necessary to perform density filtering and projection processing on the design variables during the calculation process. Helmholtz filter and hyperbolic tangent projection are used for density filtering and projection processing.

[0030] The expression for the Helmholtz filter is: , In the formula, These are the design variables before filtering; For the filtered design variables; R min The filter radius; Hyperbolic tangent projection expression: , In the formula, These are the projected design variables; For the projection point; The projection slope.

[0031] The local density field is smoothed by using filters to reduce the checkerboard effect; the design variables can be further made closer to 0 or 1 by using projection functions, making the boundary between the solid and fluid regions clearer.

[0032] After the topology optimization solution is completed, the output is the design variables after density filtering and hyperbolic tangent projection. As the final optimization result, design variables It can effectively reflect the clear boundary between the solid and fluid regions, thus outputting a topology-optimized cooling plate that satisfies maximum heat transfer and minimum flow dissipation. Topology optimization of the two-dimensional model is performed using COMSOL software, such as... Figure 3 As shown.

[0033] like Figure 4 As shown in Figure S3, the flow channel distribution of the topology-optimized cooling plate is analyzed. The channel that undertakes the main flow transport, has a wide cross-section and continuous characteristics is defined as the main channel, and the channel width of the main channel is greater than 0.5 mm. The channel that is responsible for local heat dissipation and has a narrow cross-section is defined as the branch channel, and the channel width of the branch channel is less than 0.5 mm. Then the flow channel contour is extracted.

[0034] In nature, the shape of feathers helps to efficiently guide airflow. Therefore, this invention uses feather shapes to fill solid areas. The geometric features of feathers are extracted to characterize their flow guiding properties, so that the main direction of the feather shaft matches the main direction of the flow channel. The tapered curvature is optimized to give it good flow guiding performance. The feather-like structure is adaptively arranged according to the flow channel contour constraints, replacing the original solid areas so that the feather shaft is aligned with the extension direction of the main flow channel. The arrangement position is adjusted to ensure the connectivity of the branch channels, and the volume fraction of the flow channel in the cooling plate remains unchanged at 50%. This ensures that the biomimetic structure can fill the irregular solid areas in the topological flow channel and maintain the continuity of the flow field through the natural flow guiding shape of the feathers. Finally, a biomimetic cooling plate that combines the advantages of topological optimization flow and the manufacturability of the biomimetic structure is obtained.

[0035] To verify the inheritance of the morphological features of the biomimetic cooling plate from the topology-optimized cooling plate, the box dimension method was used to calculate the fractal dimension of the flow channel structure of the two cooling plates. The box dimension can quantitatively describe the complexity and distribution characteristics of the flow channel structure. The closer the dimension values ​​are, the higher the similarity between the two in terms of geometric shape.

[0036] The flow channel structure of the two cooling plates was calculated using the box-dimensional method, and the formula is as follows: , In the formula, The side length of the box represents the scale at which the image is observed; For a side length of The minimum number of boxes required to cover the structure within the given grid.

[0037] This invention uses Python to calculate the box dimension, such as... Figure 5 As shown, the calculation results indicate that the box dimension of the topology-optimized cooling plate is... The box dimension of the biomimetic cooling plate is The difference between the two is only 0.069, which is less than 0.1, indicating that the biomimetic structure is highly consistent with the topology optimization result in terms of fractal characteristics. This result verifies that the biomimetic cooling plate improves the manufacturability of the structure while maintaining the topological flow field distribution law, and has good morphological inheritance and flow guiding performance. If the difference in box dimension between the two structures after replacing the topology cooling plate structure with the biomimetic structure is not less than 0.1, it indicates that the biomimetic structure generated by the replacement deviates from the flow characteristics of the topology optimization in terms of local geometric complexity or channel distribution. Corrections need to be made in the biomimetic replacement process to reduce the number of feather-like structures in the biomimetic structure, avoiding increased structural complexity due to excessive flow channels. Simultaneously, the direction of the feather shafts should be adjusted to ensure the dominant flow channel aligns with the main topological flow channel direction. While maintaining a fluid volume fraction of 50%, the minimum width of the branch channels should be constrained to make it similar to the width of the topology-optimized cooling plate, reducing the increase in structural complexity caused by excessively narrow local branch channels. These corrections effectively reduce the geometric complexity of the biomimetic structure, making its fractal box dimension closer to the topology optimization result, thereby improving the consistency between the two in overall shape and flow distribution.

[0038] This invention further exports the COMSOL 2D model, models it in Solidworks, and uses the professional fluid dynamics software Fluent for thermal simulation. This embodiment takes the side cooling plate of a 280Ah battery module (1P5S) as the research object. The cooling plate is a cuboid with a length of 356.19mm, a width of 207.2mm, and a thickness of 6mm. It has a 20mm long and 4mm wide inlet and outlet at the center of both sides. Simplifying it to a rectangle with a length of 356.19mm and a width of 207.2mm, and considering the symmetrical structure of the design domain, subsequent design is carried out after symmetrical division, such as... Figure 1As shown in the figure, aluminum was set as the solid material and water as the coolant in the topology optimization. The design domain was optimized with the maximum heat transfer rate and minimum flow dissipation as optimization objectives and a volume fraction of 50% as a constraint. The flow channel distribution of the topology-optimized structure was then analyzed, and the flow channel contour was extracted. The geometric features of feathers were extracted to characterize their flow guiding properties, and the feather-like structure was replaced with the irregular solid regions in the topology flow channel to obtain a manufacturable biomimetic cooling plate. The fractal box dimension was calculated to prove the similarity between the two structures. Finally, a thermal management system for a 280Ah 1P5S battery was established, and the heat dissipation performance of the traditional cooling plate, the topology-optimized cooling plate, and the biomimetic cooling plate were compared and analyzed. The results show that the biomimetic cooling plate structure has better heat dissipation performance than the traditional cooling plate and better manufacturability than the topology-optimized cooling plate, providing an important reference method for the optimized design of battery cooling plates.

Claims

1. A battery cooling plate design method integrating topology optimization and biomimetic design, characterized in that, Includes the following steps: S1. First, the three-dimensional flow and heat transfer problem of the battery cooling plate is simplified into a two-dimensional problem. The neutral layer of the cooling plate is used as the topology optimization design domain. Boundary conditions are set, and the corresponding two-dimensional topology model is established. S2. Based on the two-dimensional topology model, the Brinkman model is used to describe the flow characteristics of fluid in porous media. A multi-objective optimization function with the maximum heat transfer rate and minimum flow dissipation as the optimization objectives and a topology optimization model with volume fraction as the constraint is constructed. The two-dimensional topology model is optimized and designed to output a topology-optimized cooling plate that satisfies the maximum heat transfer rate and minimum flow dissipation. S3. Analyze the flow channel distribution of the topology-optimized cooling plate to obtain the main channel and branch channels, and then extract the flow channel contour; arrange the feather-like structure according to the flow channel contour constraint, replace the original solid area to align the feather shaft along the extension direction of the main channel, adjust the arrangement position to ensure the connectivity of the branch channels, and ensure that the volume fraction of the flow channel in the cooling plate remains unchanged, thus obtaining the biomimetic cooling plate. The flow channel structure of the topology-optimized cooling plate and the bionic cooling plate is calculated using the box dimension method. If the difference in box dimension is less than a set threshold, the corresponding bionic cooling plate is output. If the difference in box dimension is not less than the set threshold, the process returns to the replacement stage for correction, reducing the number of feather-like structures and adjusting the direction of the feather shaft to ensure that the main flow channel is aligned with the main flow channel. The minimum width of the branch channel is constrained to make the width of the branch channel similar to that of the topology-optimized cooling plate, until the difference in box dimension between the two is less than the set threshold.

2. The battery cooling plate design method integrating topology optimization and biomimetic design according to claim 1, characterized in that, In step S1, the overall structure of the cooling plate is symmetrical about the central plane in terms of geometry and boundary conditions. Half of the neutral layer of the cooling plate is selected as the topology optimization design domain, and the lower boundary is set as the symmetrical boundary condition.

3. The battery cooling plate design method integrating topology optimization and biomimetic design according to claim 1, characterized in that, Step S2 specifically includes the following steps: S21. To achieve an accurate characterization of flow resistance in porous media regions, volume force is introduced as a damping term in the incompressible laminar Navier–Stokes equations, thus obtaining the Brinkman model suitable for topology optimization calculations. S22. Based on Darcy's law, construct an expression relating volumetric force to reverse osmosis rate and flow velocity; S23. Introduce design variables to model the solid and fluid regions in a unified manner. The relationship between reverse osmosis and design variables can be expressed by Darcy interpolation equations. S24. Since there are solid and fluid regions within the design domain, and the heat transfer mechanisms of the two regions are different, design variables are used to unify the heat transfer equations of the solid and fluid regions. S25. In order to obtain a flow channel structure with maximum heat dissipation capacity and low flow power consumption, the heat transfer rate and power dissipation in the design domain are calculated, normalized and combined by weight to construct a multi-objective optimization function. S26. Analyzing the thermal conductivity, fluid flow properties, and design variables in topology optimization, the topology optimization problem is expressed as finding design variables that minimize the multi-objective optimization function under constraints. The constraints are the range of values ​​and volume fraction of the design variables. Thus, the expression for the topology optimization problem is obtained. S27. During the topology optimization solution process, density filtering and projection processing are performed on the design variables to obtain the final design variables, and the topology-optimized cooling plate that satisfies the maximum heat transfer and minimum flow dissipation is output.

4. The battery cooling plate design method integrating topology optimization and biomimetic design according to claim 3, characterized in that, In step S21, the Brinkman model is represented by the law of conservation of mass and the law of conservation of momentum describing fluid motion. The specific equations are as follows: ， ， in For gradient operators, For fluid density, For fluid velocity, Dynamic viscosity; The pressure of the fluid; For velocity gradient; It is a volume force.

5. The battery cooling plate design method integrating topology optimization and biomimetic design according to claim 4, characterized in that, The volume force in step S22 With reverse osmosis and flow rate The relational expression is: 。 6. The battery cooling plate design method integrating topology optimization and biomimetic design according to claim 5, characterized in that, The Darcy interpolation equation in step S23 is: ， ， In the formula, For solid permeability, For liquid permeability, Let Reynolds number be 1. For Darcy's number, This is the penalty factor for the Darcy interpolation model. For design variables.

7. The battery cooling plate design method integrating topology optimization and biomimetic design according to claim 6, characterized in that, The heat transfer equation in step S24 is: ， In the formula, This refers to the isobaric heat capacity of a fluid. For the design domain temperature; The thermal conductivity of the solid region is... The thermal conductivity of the fluid region. Heat generated in the solid region; when When, the heat transfer equation is transformed into the heat conduction control equation for the solid region; when When the heat transfer equation is transformed into the convection-diffusion equation of the fluid region, when When the heat transfer equation is in this state, it represents a smooth transition between the two mechanisms, reflecting the conjugate heat transfer characteristics at the solid-fluid interface.

8. The battery cooling plate design method integrating topology optimization and biomimetic design according to claim 7, characterized in that, The expression for the multi-objective optimization function in step S25 is: ， ， ， ， ， in It is the heat transfer rate within the design domain. It is the power dissipation of the flow within the design domain, for and Normalization process is performed to obtain and ,right and Weighted combination yields a multi-objective optimization function , and These are the minimum and maximum heat transfer rates; and For minimum and maximum flow power dissipation; The coefficient of an ideal heat source; As a weighting factor; For reference temperature, For designing the domain area.

9. A battery cooling plate design method integrating topology optimization and biomimetic design according to claim 8, characterized in that, The expression for the topology optimization problem in step S26 is: ， In the formula, This represents the volume fraction.

10. A battery cooling plate design method integrating topology optimization and biomimetic design according to claim 9, characterized in that, In step S27, a Helmholtz filter is used for density filtering. The expression for the Helmholtz filter is: ， In the formula, These are the design variables before filtering; For the filtered design variables; R min The filter radius; Hyperbolic tangent projection is used for projection processing. The expression for hyperbolic tangent projection is: ， In the formula, These are the projected design variables; The projection point; The projection slope.

Images