CN111737835B - Three-period minimum curved surface-based three-dimensional porous heat dissipation structure design and optimization method - Google Patents
Three-period minimum curved surface-based three-dimensional porous heat dissipation structure design and optimization method Download PDFInfo
- Publication number
- CN111737835B CN111737835B CN202010597487.1A CN202010597487A CN111737835B CN 111737835 B CN111737835 B CN 111737835B CN 202010597487 A CN202010597487 A CN 202010597487A CN 111737835 B CN111737835 B CN 111737835B
- Authority
- CN
- China
- Prior art keywords
- function
- optimization
- wall thickness
- porous
- porous structure
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 83
- 230000017525 heat dissipation Effects 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000013461 design Methods 0.000 title claims abstract description 28
- 230000000737 periodic effect Effects 0.000 claims abstract description 33
- 230000008569 process Effects 0.000 claims abstract description 15
- 239000000463 material Substances 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims description 13
- 230000004907 flux Effects 0.000 claims description 12
- 239000013598 vector Substances 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 7
- 239000011343 solid material Substances 0.000 claims description 6
- 239000011148 porous material Substances 0.000 claims description 5
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 108010074864 Factor XI Proteins 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- HDRXZJPWHTXQRI-BHDTVMLSSA-N diltiazem hydrochloride Chemical compound [Cl-].C1=CC(OC)=CC=C1[C@H]1[C@@H](OC(C)=O)C(=O)N(CC[NH+](C)C)C2=CC=CC=C2S1 HDRXZJPWHTXQRI-BHDTVMLSSA-N 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims 1
- 230000035945 sensitivity Effects 0.000 claims 1
- 238000004519 manufacturing process Methods 0.000 abstract description 9
- 238000007639 printing Methods 0.000 abstract description 6
- 238000010146 3D printing Methods 0.000 abstract description 4
- 238000011960 computer-aided design Methods 0.000 abstract description 2
- 238000002474 experimental method Methods 0.000 description 3
- 238000010206 sensitivity analysis Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/10—Additive manufacturing, e.g. 3D printing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
A three-cycle minimum curved surface-based design and optimization method for a three-dimensional porous heat dissipation structure belongs to the field of computer aided design. First, a porous structure is built by an implicit function representation of a three-cycle infinitesimal surface. Secondly, according to a steady-state heat conduction equation, the heat dissipation problem is converted into a minimization problem of heat dissipation weakness under a given constraint condition. Then, a global-local interpolation method is used for directly solving the parameter function. And finally, carrying out periodic optimization and wall thickness optimization on the modeling problem to obtain the optimized porous sheath-like structure with smooth period and wall thickness variation. The porous structure of the present invention greatly improves heat dissipation performance, heat conduction efficiency and efficiency. The porous structure designed by the invention has the characteristics of smoothness, full connectivity, controllability, quasi-self-supporting property and the like, the applicability and manufacturability of the structure are ensured by the characteristics, the porous structure is suitable for a 3D printing manufacturing method, the internal structure in the printing process does not need to be additionally supported, and the printing time and the printing material are saved.
Description
Technical Field
The invention belongs to the field of engineering design and manufacture, and relates to a design and optimization method of a three-dimensional porous heat dissipation structure, which is suitable for heat dissipation structures of various large-scale engineering machinery, radiators of automobiles and the like and related parts of combustion appliances.
Background
How to construct a light and efficient heat dissipation structure has received wide attention in various engineering fields. The conventional radiator structure cannot achieve high heat conduction efficiency. Through the research on the porous structure, the heat-conducting property of the structure can be effectively improved, however, the representation and optimization method becomes a technical bottleneck restricting the further development of the structure.
The traditional heat dissipation structure design depends on the basic theory and practical experience of calorifics, and the heat dissipation problem of a complex structure is difficult to solve. Subsequently, topology optimization methods using porous structures such as tree topologies, truss/frame structures, microstructures, etc. have emerged to address the above issues, and these porous structures can be used to calculate high-degree-of-freedom topologies for cooling channels. However, these methods have a common problem in that they require a large number of design variables and are expensive to optimize due to time-consuming re-gridding.
In recent years, porous structures based on three-cycle extremely small curved surfaces have been widely used in the fields of tissue engineering, lightweight manufacturing, biomedicine, and the like. The porous structure based on the three-cycle extremely-small curved surface has the advantages of good connectivity, easiness in control, high specific strength, rigidity and the like. The porous structure based on the three-cycle infinitesimal surface (especially for large and complex porous structures) is very time-consuming and memory-consuming to represent by a polyhedral mesh (tetrahedron or hexahedron), and the traditional finite element-based porous structure processing method is almost heuristic without effective optimization, so that the porous structure based on the three-cycle infinitesimal surface is less researched in heat dissipation.
In view of the above, an efficient representation and optimization method is proposed to obtain the period and wall thickness of a three-cycle minimum curved porous sheath structure suitable for heat dissipation. The main optimization processes include period optimization and wall thickness optimization. The former is a coarse adjustment of the structure and the latter is a fine adjustment of the structure. Firstly, expressing a porous sheath-like structure by using an implicit function, wherein the implicit function is controlled by a periodic parameter function and a wall thickness parameter function; on the basis, a steady-state heat conduction equation with boundary conditions can be conveniently and directly established into a mathematical model by using a function expression, and then the optimization problem of the model is converted into the solution of the two continuous parameter functions; and finally, from the perspective of discretization, the two parameter functions can be effectively calculated by utilizing implicit function representation and radial basis interpolation without re-meshing, and finally the optimized heat dissipation porous structure with smooth period and wall thickness variation is obtained.
Disclosure of Invention
The invention provides an effective expression and optimization method of a porous heat dissipation structure based on a three-cycle extremely-small curved surface. First, a porous structure is built by an implicit function representation of a three-cycle infinitesimal surface. Secondly, according to a steady-state heat conduction equation, the heat dissipation problem is converted into a minimization problem of heat dissipation weakness under a given constraint condition. Then, a global-local interpolation method is used for directly solving the parameter function. Finally, we perform periodic optimization and wall thickness optimization on the modeling problem to obtain an optimized porous sheath-like structure with smooth period and wall thickness variation.
The technical scheme adopted by the invention is as follows:
a three-cycle extremely-small-curved-surface-based design and optimization method for a three-dimensional porous heat dissipation structure comprises the following steps:
method for expressing porous structure
The commonly used three-cycle minimum curved surface mostly has implicit function representation, taking the P minimum curved surface as an example:
wherein r is a three-dimensional vector, and x, y and z are corresponding coordinates of the three-dimensional vector.
The periodic parameter function P (r) >0 can be directly added into the function representation of the three-period extremely-small curved surface, and in order to make the distance scale of the directional distance field in the periodic variation process basically unchanged, the implicit function representation is improved as follows:
wherein P (r) controls the continuous change of the pore period, and a porous curved surface with smooth transition on the space is constructed; other types of three-cycle pole faces were treated in the same manner.
The porous structure with thickness based on the three-cycle extremely-small curved surface can be used for the improved implicit function curved surface by controlling the parameter function W (r) of the wall thicknessTwo offset surfaces are obtained after the offset to two sides, and are expressed as:
finally, the three-cycle minimal surface-based porous sheath-like structure is represented by an intersection operator continuous function:
in the above definition, the parameter functions p (r) >0 for controlling the period and the parameter functions w (r) >0 for controlling the wall thickness are introduced to realize the control of the shape and the periodic pores of the porous structure, and the porous structure with the wall thickness meeting the requirement is finally generated by optimizing the parameter functions p (r) and w (r).
The porous structure defined by the function inherits the excellent characteristics of the three-period extremely-small curved surface, such as high surface area-volume ratio, full connectivity, high smoothness and controllability. High surface to volume ratio and full connectivity facilitate heat dissipation from the structure. The structural function provides a computationally optimized method based on the high degree of controllability of the period and wall thickness. Good smoothness and connectivity facilitates 3D printing manufacturing, ensures manufacturing accuracy, and can remove excess material (e.g., excess liquid in SLA) in 3D manufacturing.
(II) optimization procedure for heat dissipation problem
The invention mainly focuses on the heat dissipation problem under the condition of steady-state heat conduction, after a model heat source and boundary conditions are given, the internal space of the model is filled by using the constructed porous sheath-like structure, and the optimal distribution of the porous structure period and the wall thickness is solved under the conditions of given material volume constraint and periodic function gradient constraint.
1. Problem modeling
Based on the above purpose, the heat dissipation problem model is established as follows:
such that:
where C is the heat dissipation weakness, T is the temperature field, Ω is the given design area, Φ is the functional representation of the porous sheath-like structure given above, Q is the heat flux density of the internal heat source, Q is the heat flux density of the internal heat sourcesIs the Neumann boundary gammaQThe heat flux in the upper direction of the normal line,is a given temperature on the Directilet boundary, λ is the thermal conductivity;is a vector differential operator which is a function of,x, Y, Z represent unit vectors in the positive direction of the three coordinate axes x, y, and z, respectively;is the corresponding test function and is,Sob1is a first order Sobelov space, V is the volume of the porous structure,is corresponding volume constraint, and adds gradient constraint of periodic variation to avoid the damage of porous structure caused by severe variation of periodic functionAnd formula for calculation of a mode having a gradientH (x) is the Heaviside function, H (x) is 0 when x is negative, otherwise 1, and H (x) is defined as the continuous function H in order to minimize the optimization problem and avoid checkerboard phenomenaη(x) The function is defined as:
where η is a regularization parameter used to control the number of non-zero elements in the global stiffness matrix, and is typically taken to be 10-3An interval of intermediate values is defined. Furthermore, the thermal conductivity λ of the material of the porous structure is calculated from the structure function Φ and should be set as:ξ=H(Φ) Is the volume ratio of the solid material, lambdaSAnd λDRepresenting the thermal conductivity of the solid material and the pore portion, respectively.
2. Discretization
The method is characterized in that a double-scale grid is adopted in the discretization process, the three-dimensional design space of a design domain is divided into uniform hexahedron finite elements called coarse units firstly, the coarse units are used for generating a temperature field, and the number n of the coarse unitssIs determined by the volume of the design space; each coarse cell is then further subdivided into smaller hexahedral cells, called fine cells, which are used for more precise geometric calculations of volume, etc., where the number n of fine cells in each coarse cellbThe default setting is 27. A discrete form of the optimization problem (1.6-1.9) is obtained:
such that:
KT=Q (1.12)
where T is the temperature field, Q is the heat source and heat flux terms, K is the stiffness matrix, V is the volume of the porous structure,is a corresponding volume constraint, Nb=nb×nsIs the total number of the thin cells,is the function value of phi of the l-th node in the j-th cell, vbIs the volume of the fine grid cell, G is the global gradient constraint of the structure, | Ω | is the volume of Ω, nlIs made by designingThe number of sub-regions in the domain,is the ith sub-region omegaiThe number of the thin units in the inner,is the gradient of the periodic function at the ith point in the s-th cell,is the ith sub-region omegaiThe local gradient within the constraint value of the local gradient,is the volume of the ith coarse grid cell, p>0 is a penalty factor for the global gradient constraint and has:
3. global-local interpolation
The optimization of a period parameter function and a thickness parameter function is converted into the optimization of a limited number of design variables by adopting a global-local radial basis interpolation algorithm, and the key idea is to decompose a large coefficient matrix into a weighted small coefficient matrix for solving.
Taking the periodic parameter function as an example, Ω is first divided into nlSub-regionsThe radial basis interpolation in the local ellipsoid (containing the corresponding sub-region) yields the local periodic parameter function:
wherein psik(r) is a number fromk(r) defined weight parameter, dk(r)=‖r-Ck‖2Is interpolated to the ellipsoid central point CkDistance of (1) ()+Is that the truncation function satisfies x>Time 0 (x)+X, otherwise (x)+=0,Rk(r) is a length function with respect to the radius, Pk(r) is the subregion ΩkA corresponding local periodic parameter function within the local ellipsoid, and defined as:
wherein R iski(r)=(r-Oki)2log(|r-Oki|) is the sheet radial basis function,is the sub-region omegakControl points in the corresponding local ellipsoid, qki(r) is a primary term of coordinates x, y, z, akiAnd bkjThe coefficients to be found are the quadratic term and the primary term, respectively, and m is the number of primary terms (m is 4 by default).
Global-local radial basis interpolation can be simplified to the following form:
wherein n istIs the total number of control points in the design field omega (generally 400), Ni(r) is the corresponding calculable coefficient function,is a function of the period of the control point. The proposed global local interpolation method can improve the calculation efficiency and make the structure change smoothly.
4. Modeling problem optimization
The three-dimensional heat dissipation optimization method based on the constructed optimization problem comprises two parts of period optimization and wall thickness optimization. The period and the wall thickness of the porous structure based on the three-period extremely-small curved surface are independently controlled by a period function P (r) and a wall thickness function W (r), the period optimization is coarse adjustment of the structure, the wall thickness optimization is fine adjustment, and the specific optimization process is as follows:
step 1: optimizing the period; firstly, converting function optimization into optimization of interpolation basis function parameters by using a radial basis interpolation method; randomly selecting n in the solution domaintA base point of interpolationThen there is the interpolation form:
thus, the optimization problem is transformed into a parametric variableThe optimization problem of (2); finally, the optimization variables are derived by the objective function and the constraint function as follows:
wherein,the objective function, the volume constraint and the gradient constraint respectively to the parametric variable PiThe equation for the partial derivative is solved,is an intermediate equation to be calculated in the process of gradient partial derivative calculation; n issIs the number of coarse cells, nbIs the number of thin cells in each thick cell;is phi function value of the l node in the k coarse unit and the j fine unit; from the formula of thermal conductivityXi is knownkjIs the heat conductivity lambda of the k-th coarse cell and the j-th fine cellkjA corresponding parameter factor xi; k0Is the initial stiffness matrix. In MMA solver, byAn optimized porous structure with a smoothly varying period can be obtained. Since the wall thickness function w (r) is fixed, the structural porosity as a whole increases with increasing periodic function p (r), and therefore a period-optimized convergence is easy to achieve. In our experiment, the cycle optimization converged to 70 iterations.
Step 2: optimizing the wall thickness; similarly, the control point based on W (r) (variable is) And constructing a wall thickness function W (r) by adopting a radial basis interpolation method, wherein the corresponding sensitivity analysis is as follows:
wherein,respectively, an objective function, a volume constraint, and a parametric variable WiSolving an equation of partial derivative;is phi function value of the l node in the k coarse unit and the j fine unit; xikjIs the heat conductivity lambda of the k-th coarse cell and the j-th fine cellkjThe corresponding parameter factor. Since the wall thickness variation is smoother than the periodic variation, the gradient constraint of w (r) is no longer required. Finally, selection is made in the MMA solverAndan optimized porous structure with both smooth periodicity and wall thickness variation can be obtained. Since the optimization period function p (r) is fixed and the structural porosity monotonically increases with the increase of the wall thickness function w (r), the convergence of the wall thickness optimization is also better achieved. In the experiment, the wall thickness optimization converged to 30 iterations.
The invention discloses a design and optimization system of a porous sheath-shaped heat dissipation structure for 3D printing, and belongs to the field of computer-aided design and industrial design and manufacture. The proposed porous structure is expressed in an implicit function form, and has good connectivity, controllability, mechanical property, thermal property, higher surface area to volume ratio and smoothness. The porous structure is applied to the three-dimensional heat dissipation problem, and an optimized porous structure with continuous geometric change and smooth topological change is obtained. Compared with the existing traditional heat dissipation structure, the porous structure greatly improves the heat dissipation performance, the heat conduction efficiency and the efficiency. The porous structure designed by the invention has the characteristics of smoothness, full connectivity, quasi-self-supporting property and the like, the applicability and manufacturability of the structure are ensured, the porous structure is suitable for a common 3D printing manufacturing method, the internal structure in the printing process does not need to be additionally supported, and the printing time and the printing material can be saved.
Drawings
Fig. 1 is a flow chart of the design and optimization of a three-dimensional porous heat dissipation structure based on three-cycle extremely-small curved surfaces.
Fig. 2 is a diagram of the design and optimization results of a three-cycle minimal curved surface-based three-dimensional porous heat dissipation structure, and a, b, and c are diagrams of the optimization results of three different three-cycle minimal curved surfaces.
Detailed Description
The following further describes a specific embodiment of the present invention with reference to the drawings and technical solutions.
The implementation of the invention can be divided into a plurality of main steps of porous sheath structure function representation, establishment of a heat dissipation problem optimization model and discretization thereof, flow optimization and the like:
method for expressing porous shell-like structure
Firstly, establishing an improved implicit numerical surface:
wherein r is a three-dimensional vector, x, y and z are respectively corresponding coordinates of the three-dimensional vector, and P (r) controls the continuous change of the hole period, so that a hole curved surface with smooth transition in space is constructed.
Further, a multi-scale porous sheath-like structure with thickness is constructed: the porous structure with the thickness based on the three-cycle extremely-small curved surface can be obtained by shifting the improved implicit function curved surface to two sides by utilizing a parameter function W (r) for controlling the wall thickness, wherein the two shifted curved surfaces are expressed as follows:
finally obtaining a porous sheath-like structure based on the three-cycle extremely-small curved surface:
to sum up, the value intervals of P (r) (P curved surface are [0.5,2], G curved surface are [0.37,2], D curved surface are [0.5,2], IWP curved surface is [0.48,2]) control the period of the porous structure, w (r) (P curved surface is [0.02,0.95], G curved surface is [0.02,1.35], D curved surface is [0.02,0.7], IWP curved surface is [0.02,2.95]) control the wall thickness of the porous structure.
(II) modeling and optimization based on porous shell-like structure
1. Modeling of heat dissipation problems
The problem of minimizing the weakness of heat dissipation is used here to create a porous structure optimization problem. The inner space of the model is filled by using the constructed porous sheath-like structure with the aim of minimizing the average temperature of the structure and with the constraint of the model volume and boundary conditions, so that the period and the wall thickness of the porous structure have optimal distribution under the condition of the constraint of the material volume.
Based on the above objective, the problem model is established as follows:
such that:
where C is the heat dissipation weakness, T is the temperature field, Ω is the given design area, Φ is the functional representation of the porous sheath-like structure given above, Q is the heat flux density of the internal heat source, Q is the heat flux density of the internal heat sourcesIs the Neumann boundary gammaQThe heat flux in the upper direction of the normal line,is a given temperature on the Directilet boundary, λ is the thermal conductivity;is a vector differential operator which is a function of,x, Y, Z represent unit vectors in the positive direction of the three coordinate axes x, y, and z, respectively;is the corresponding test function and is,Sob1is a first order Sobelov space, V is the volume of the porous structure,is corresponding volume constraint, and adds gradient constraint of periodic variation to avoid the damage of porous structure caused by severe variation of periodic functionAnd formula for calculation of a mode having a gradientH (x) is the Heaviside function, H (x) is 0 when x is negative, otherwise 1, and H (x) is defined as the continuous function H in order to minimize the optimization problem and avoid checkerboard phenomenaη(x) The function is defined as:
where η is a regularization parameter used to control the number of non-zero elements in the global stiffness matrix, and is generally taken as 10-3An interval of intermediate values is defined. Furthermore, the thermal conductivity λ of the material of the porous structure is calculated from the structure function Φ and should be set as:xi H (Φ) is the volume ratio of the solid material, λSAnd λDRepresenting the thermal conductivity of the solid material and the pore portion, respectively.
2. Discretization of optimization problems
In the discretization process, a solving area is subdivided into two sets of uniform grids with different precisions: the temperature field is interpolated with a coarse grid, the model is described with a fine grid and an integral calculation is performed.
Obtaining a discrete form of the optimization problem:
such that:
KT=Q (2.11)
where T is the temperature field, Q is the heat source and heat flux terms, K is the stiffness matrix, V is the volume of the porous structure,is a corresponding volume constraint, Nb=nb×nsIs the total number of the thin cells,is the function value of phi of the l-th node in the j-th cell, vbIs the volume of the fine grid cell, G is the global gradient constraint of the structure, | Ω | is the volume of Ω, nlIs the number of sub-regions in the design domain,is the ith sub-region omegaiThe number of the thin units in the inner,is the gradient of the periodic function at the ith point in the s-th cell,is the ith sub-region omegaiThe local gradient within the constraint value of the local gradient,is the volume of the ith coarse grid cell, p>0 is a penalty factor for the global gradient constraint and has:
the optimization of the period parameter function and the thickness parameter function is converted into the optimization of a limited number of design variables by adopting a global-local RBF interpolation algorithm, and the global-local radial basis interpolation can be simplified into the following form:
wherein n istIs the total number of control points in the design field omega (value 400), Ni(r) is the corresponding calculable coefficient function,is the value of the periodic function of the control point,is the value of the wall thickness function at the control point. Since the position of the control points is not changed during the optimization process, the coefficient function Ni(r) may be calculated prior to optimization.
3. Modeling problem optimization
Only two unknown parameter functions p (r) and w (r) need to be optimized here. The specific optimization process is as follows:
step 1: optimizing the period; firstly, converting function optimization into optimization of interpolation basis function parameters by using a radial basis interpolation method; randomly selecting n in the solution domaintA base point of interpolationThere is an interpolation form of (2.14). Thus, the optimization problem is transformed into a parametric variableThe optimization problem of (2); finally, the optimization variables are derived by the objective function and the constraint function as follows:
wherein,the objective function, the volume constraint and the gradient constraint respectively to the parametric variable PiThe equation for the partial derivative is solved,is an intermediate equation to be calculated in the process of gradient partial derivative calculation; n issIs the number of coarse cells, nbIs the number of thin cells in each thick cell;is phi function value of the l node in the k coarse unit and the j fine unit; from the formula of thermal conductivityXi is knownkjIs the heat conductivity lambda of the k-th coarse cell and the j-th fine cellkjA corresponding parameter factor xi; k0Is the initial stiffness matrix. Given the variable sensitivity analysis, the optimized periodic parameter function can be obtained by the well-known MMA method, thus yielding a periodically optimized structure and serving as the initial structure for wall thickness optimization.
Step 2: optimizing the wall thickness; similarly, the control point based on W (r) (variable is) And constructing a wall thickness function W (r) by adopting a radial basis interpolation method, wherein the corresponding sensitivity analysis is as follows:
wherein,respectively, an objective function, a volume constraint, and a parametric variable WiSolving an equation of partial derivative;is phi function value of the l node in the k coarse unit and the j fine unit; xikjIs the heat conductivity lambda of the k-th coarse cell and the j-th fine cellkjThe corresponding parameter factor. Substituting into MMA algorithm to obtain the solution of optimization problem.
And the design and optimization of a heat dissipation structure based on the porous structure representation of the three-period extremely-small curved surface are provided. The porous structure is expressed in an implicit function form, and has good connectivity, controllability, higher surface-to-volume ratio, higher smoothness and good mechanical and thermal properties. Various experiments show that the proposed porous structure greatly improves heat dissipation performance, efficiency of heat conduction and efficiency. To achieve high thermal conductivity efficiency, there is a balance between period and wall thickness for a given volume constraint, and the optimized structure period and wall thickness variation is smooth and natural, which facilitates structural stress and fabrication. Compared with the traditional heat dissipation structure and the grid structure, the optimized porous structure has higher heat dissipation efficiency (lower heat dissipation weakness).
Claims (1)
1. A three-cycle extremely-small-curved-surface-based design and optimization method for a three-dimensional porous heat dissipation structure is characterized by comprising the following steps:
method for expressing porous structure
The three-cycle infinitesimal surfaces all have implicit function representations, and the implicit function representation of the P infinitesimal surface is as follows:
wherein r is a three-dimensional vector, and x, y and z are corresponding coordinates of the three-dimensional vector respectively;
the periodic parameter function P (r) >0 is directly added into the function representation of the three-period extremely-small curved surface, and in order to enable the distance scale of the directional distance field to be unchanged in the periodic change process, the latent function is improved and represented as follows:
wherein P (r) controls the continuous change of the pore period, and a porous curved surface with smooth transition on the space is constructed;
the three-cycle minimum curved surface-based porous structure with thickness is subjected to the improved implicit function curved surface by controlling the parameter function W (r) of the wall thicknessTwo offset curved surfaces are obtained after the offset to the two sides, and are expressed as:
finally, the three-cycle minimal surface-based porous sheath-like structure is represented by an intersection operator continuous function:
introducing a parameter function P (r) >0 for controlling the period and a parameter function W (r) >0 for controlling the wall thickness to realize the control of the shape and the periodic holes of the porous structure, and finally generating the porous structure with the wall thickness meeting the requirement by optimizing the parameter functions P (r) and W (r);
(II) optimization procedure for heat dissipation problem
In the heat dissipation problem under the condition of steady-state heat conduction, after a model heat source and boundary conditions are given, the internal space of the model is filled by using the porous sheath-like structure, and the optimal distribution of the period and the wall thickness of the porous structure is solved under the conditions of volume constraint and periodic function gradient constraint of given materials;
(1) problem modeling
The heat dissipation problem model is established as follows:
such that:
where C is the heat dissipation weakness, T is the temperature field, Ω is the given design area, Φ is the functional representation of the porous sheath-like structure given above, Q is the heat flux density of the internal heat source, Q is the heat flux density of the internal heat sourcesIs the Neumann boundary gammaQThe heat flux in the upper direction of the normal line,is a given temperature on the Directilet boundary, λ is the thermal conductivity;is a vector differential operator which is a function of,x, Y, Z represent unit vectors in the positive direction of the three coordinate axes x, y, and z, respectively;is the corresponding test function and is,Sob1is a first order Sobelov space, V is the volume of the porous structure,is corresponding volume constraint, and adds gradient constraint of periodic variation to avoid the damage of porous structure caused by severe variation of periodic functionAnd formula for calculation of a mode having a gradientH (x) is the Heaviside function, H (x) is 0 when x is negative, otherwise 1, and H (x) is defined as the continuous function H in order to minimize the optimization problem and avoid checkerboard phenomenaη(x) The function is defined as:
wherein eta is a regularization parameter for controlling the number of non-zero elements in the global stiffness matrix, and the parameter eta is 10-3Defining an interval of intermediate values; furthermore, the material thermal conductivity λ of the porous structure is calculated from the structural function Φ, set as:xi H (Φ) is the volume ratio of the solid material, λSAnd λDRespectively representing the thermal conductivity of the solid material and the hole part;
(2) discretization
In the discretization process, a solving area is subdivided into two sets of uniform grids with different precisions:using a coarse grid to interpolate a temperature field, using a fine grid to describe a model and performing integral calculation; the number of coarse cells is nsNumber of fine cells n in each coarse cellbDefault setting is 27; a discrete form of the optimization problem (1.6-1.9) is obtained:
such that:
KT=Q(1.12)
where T is the temperature field, Q is the heat source and heat flux terms, K is the stiffness matrix, V is the volume of the porous structure,is a corresponding volume constraint, Nb=nb×nsIs the total number of the thin cells,is the function value of phi of the l-th node in the j-th cell, vbIs the volume of the fine grid cell, G is the global gradient constraint of the structure, | | Ω | | is the volume of Ω, nlIs the number of sub-regions in the design domain,is the ith sub-region omegaiThe number of the thin units in the inner,is that the periodic function is in the s-th detailThe gradient of the ith point within the element,is the ith sub-region omegaiThe local gradient within the constraint value of the local gradient,is the volume of the ith coarse grid cell, p>0 is a penalty factor for the global gradient constraint and has:
(3) global-local interpolation
The optimization of a periodic parameter function and a thickness parameter function is converted into the optimization of a limited number of design variables by adopting a global-local radial basis interpolation algorithm, and the key idea is to decompose a large coefficient matrix into a weighted small coefficient matrix for solving;
periodic parametric function, first dividing Ω into nlSub-regionsAnd (3) obtaining a local periodic parameter function by radial basis interpolation in a local ellipsoid containing the corresponding sub-region:
wherein psik(r) is a number fromk(r) defined weight parameter, dk(r)=||r-Ck||2Is interpolated to the ellipsoid central point CkDistance of (1) ()+Is that the truncation function satisfies x>Time 0 (x)+X, otherwise (x)+=0,Rk(r) is a length function with respect to the radius, Pk(r) is the subregion ΩkA corresponding local periodic parameter function within the local ellipsoid, and defined as:
wherein R iski(r)=(r-Oki)2log(|r-Oki|) is the sheet radial basis function,is the sub-region omegakControl points in the corresponding local ellipsoid, qki(r) is a primary term of coordinates x, y, z, akiAnd bkjThe coefficients to be solved of the secondary term and the primary term are respectively, m is the number of the primary term, and m is 4;
global-local radial basis interpolation is simplified to the form:
wherein n istIs the total number of control points in the design domain omega, and takes 400, Ni(r) is a function of the corresponding coefficients,is the period function value of the control point;
(4) modeling problem optimization
The three-dimensional heat dissipation optimization method based on the constructed optimization problem comprises a period optimization part and a wall thickness optimization part; the period and the wall thickness of the porous structure based on the three-period extremely-small curved surface are independently controlled by a period function P (r) and a wall thickness function W (r), the period optimization is coarse adjustment of the structure, the wall thickness optimization is fine adjustment, and the specific optimization process is as follows:
step 1: optimizing the period; firstly, optimizing the function by utilizing a radial basis interpolation methodConverting the parameters into interpolation basis function parameters for optimization; randomly selecting n in the solution domaintA base point of interpolationThen there is the interpolation form:
conversion of periodic optimization problem into parametric variablesThe optimization problem of (2); finally, the optimization variables are derived by the objective function and the constraint function as follows:
wherein,the objective function, the volume constraint and the gradient constraint respectively to the parametric variable PiThe equation for the partial derivative is solved,is an intermediate equation to be calculated in the process of gradient partial derivative calculation; n issIs the number of coarse cells, nbIs the number of thin cells in each thick cell;is phi function value of the l node in the k coarse unit and the j fine unit; from the formula of thermal conductivityξkjIs the heat conductivity lambda of the k-th coarse cell and the j-th fine cellkjA corresponding parameter factor xi; k0Is an initial stiffness matrix; in MMA solver, byObtaining an optimized porous structure with a stable period;
step 2: optimizing the wall thickness; similarly, based on the control point of W (r), the variable isAnd (3) constructing a wall thickness function W (r) by adopting a radial basis interpolation method, and analyzing the corresponding sensitivity as follows:
wherein,respectively, an objective function, a volume constraint, and a parametric variable WiSolving an equation of partial derivative;is phi function value of the l node in the k coarse unit and the j fine unit; xikjIs the heat conductivity lambda of the k-th coarse cell and the j-th fine cellkjCorresponding parameter factors; since the wall thickness variation is smoother than the periodic variation, the gradient constraint of w (r) is no longer required; finally, selection is made in the MMA solverAndan optimized porous structure is obtained with both smooth periodicity and wall thickness variation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010597487.1A CN111737835B (en) | 2020-06-28 | 2020-06-28 | Three-period minimum curved surface-based three-dimensional porous heat dissipation structure design and optimization method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010597487.1A CN111737835B (en) | 2020-06-28 | 2020-06-28 | Three-period minimum curved surface-based three-dimensional porous heat dissipation structure design and optimization method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111737835A CN111737835A (en) | 2020-10-02 |
CN111737835B true CN111737835B (en) | 2021-10-15 |
Family
ID=72651390
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010597487.1A Active CN111737835B (en) | 2020-06-28 | 2020-06-28 | Three-period minimum curved surface-based three-dimensional porous heat dissipation structure design and optimization method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111737835B (en) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112191849B (en) * | 2020-10-10 | 2023-03-28 | 重庆邮电大学 | Gradient porous heat dissipation device design and material increase manufacturing method based on temperature distribution |
CN112395746B (en) * | 2020-11-03 | 2023-12-08 | 山东大学 | Method, microstructure, system and medium for calculating microstructure family equivalent material property |
CN114621013B (en) * | 2020-12-11 | 2023-02-07 | 上海微电子装备(集团)股份有限公司 | Periodic structure, wafer bearing table and manufacturing method thereof |
CN113158526B (en) * | 2021-05-07 | 2024-03-19 | 杭州电子科技大学 | TPMS-based controllable pore structure design method |
CN113392570B (en) * | 2021-06-08 | 2022-10-18 | 哈尔滨工业大学 | Method for evaluating homogeneity degree of pore structure of cement-based material particle stacking system |
CN113722779A (en) * | 2021-09-08 | 2021-11-30 | 大连理工大学 | Parameterized carving design method based on thin shell structure |
CN114037811B (en) * | 2021-11-18 | 2022-05-27 | 北京优锘科技有限公司 | 3D temperature field graph rendering method, apparatus, medium, and device based on directed distance field |
CN114589783B (en) * | 2022-03-11 | 2023-04-07 | 山东大学 | 3D printing porous communicated reef generation method and system |
CN114623705B (en) * | 2022-03-15 | 2022-10-18 | 大连理工大学 | Heat exchanger based on GD type hybridization minimum curved surface disturbance structure |
CN115017560A (en) * | 2022-03-30 | 2022-09-06 | 西安交通大学 | 3D printing-oriented multi-scale multi-material self-adaptive design method |
CN114970252A (en) * | 2022-05-12 | 2022-08-30 | 重庆大学 | Variable porosity three-cycle minimum curved surface porous catalyst carrier and establishment method thereof |
CN117911651A (en) * | 2024-01-29 | 2024-04-19 | 上海慕灿信息科技有限公司 | Three-dimensional structured sub-grid dividing method |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017192508A1 (en) * | 2016-05-03 | 2017-11-09 | Masdar Institute Of Science And Technology | Catalytic converter substrates comprising triply periodic minimal surfaces |
CN107885945A (en) * | 2017-11-22 | 2018-04-06 | 大连理工大学 | Towards the multiple dimensioned pore space structure lightweight modeling method of 3D printing |
CN109145427A (en) * | 2018-08-14 | 2019-01-04 | 大连理工大学 | A kind of porous structure design and optimization method based on three period minimal surfaces |
CN110377960A (en) * | 2019-06-24 | 2019-10-25 | 浙江大学 | The building storage method of the non-homogeneous porous support of biological tissue based on B-spline body |
CN111062166A (en) * | 2019-12-16 | 2020-04-24 | 浙江大学 | Three-period minimum curved surface porous structure topology optimization method based on variable density method |
CN111159903A (en) * | 2019-12-31 | 2020-05-15 | 重庆邮电大学 | Design and manufacturing method of compact type multi-channel multi-fluid heat exchange device |
-
2020
- 2020-06-28 CN CN202010597487.1A patent/CN111737835B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017192508A1 (en) * | 2016-05-03 | 2017-11-09 | Masdar Institute Of Science And Technology | Catalytic converter substrates comprising triply periodic minimal surfaces |
CN107885945A (en) * | 2017-11-22 | 2018-04-06 | 大连理工大学 | Towards the multiple dimensioned pore space structure lightweight modeling method of 3D printing |
CN109145427A (en) * | 2018-08-14 | 2019-01-04 | 大连理工大学 | A kind of porous structure design and optimization method based on three period minimal surfaces |
CN110377960A (en) * | 2019-06-24 | 2019-10-25 | 浙江大学 | The building storage method of the non-homogeneous porous support of biological tissue based on B-spline body |
CN111062166A (en) * | 2019-12-16 | 2020-04-24 | 浙江大学 | Three-period minimum curved surface porous structure topology optimization method based on variable density method |
CN111159903A (en) * | 2019-12-31 | 2020-05-15 | 重庆邮电大学 | Design and manufacturing method of compact type multi-channel multi-fluid heat exchange device |
Non-Patent Citations (2)
Title |
---|
基于三周期极小曲面的模型空间剖分研究;周丽超;《中国优秀硕士学位论文全文数据库 信息科技辑》;20200116;全文 * |
基于三周期极小曲面的模型轻量化方法;冯紫鑫;《中国优秀硕士学位论文全文数据库 信息科技辑》;20180316;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111737835A (en) | 2020-10-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111737835B (en) | Three-period minimum curved surface-based three-dimensional porous heat dissipation structure design and optimization method | |
WO2022000132A1 (en) | Method for designing and optimizing three-dimensional porous heat dissipation structure on basis of three-cycle minimal curved surface | |
CN109145427B (en) | Porous structure design and optimization method based on three-cycle minimum curved surface | |
CN111709171B (en) | Isogeometric solving and heat dissipation topology generation method for heat flow strong coupling problem | |
CN111859693B (en) | Efficient porous structure representation and optimization method | |
CN110069800B (en) | Three-dimensional structure topology optimization design method and equipment with smooth boundary expression | |
CN109657284B (en) | Metamaterial-oriented equal-geometry topology optimization method | |
CN109376497B (en) | Method for obtaining continuous gradient porous structure of extremely-small curved surface | |
CN113204906B (en) | Multiphase material topology optimization design method and system considering structural stability | |
CN112765856B (en) | Mixed level set method for topological optimization of functionally graded porous structure | |
CN110795873A (en) | Cross-scale topology optimization method considering size control | |
CN113094944A (en) | Micro-channel radiator and microscopic scale topology optimization design method thereof | |
CN109670207A (en) | A kind of dynamics integrated design method towards a variety of porous material structures | |
CN114254408A (en) | Gradient lattice isogeometric topology optimization method based on proxy model | |
CN114254409B (en) | Multi-scale topology optimization method based on isogeometric analysis | |
CN109783950A (en) | The method of topological optimization design of connectivity structure in increasing material manufacturing | |
CN110245410B (en) | Multi-parametric variable-based topological optimization design method for thermoelastic structure of multiphase material | |
CN117373579B (en) | Multi-microstructure multi-scale parallel topology optimization method and system under action of dynamic load in time domain | |
CN113191016B (en) | Body expression model-based multi-material product modeling and analyzing integrated method | |
Liu et al. | An efficient data-driven optimization framework for designing graded cellular structures | |
CN116187073A (en) | Topology optimization method for anisotropic material transient heat transfer structure based on grid-free EFGM | |
CN113515824B (en) | Topology optimization design method for cooperation of rib layout and substrate shape | |
CN111859482B (en) | Lightweight design method for muscle type thin-wall structure | |
CN113722779A (en) | Parameterized carving design method based on thin shell structure | |
CN116579151B (en) | Non-uniform lattice structure optimization design method based on MMC framework |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |