A Topological Optimization Method and System for Precise Microstructure Filling

Abstract

The present invention discloses a topology optimization method and system for precise microstructure filling. The method includes two aspects: a synthetic microstructure interpolation method based on an initial microstructure and an update method for the unit filling microstructure during the optimization process. The synthetic microstructure based on the initial microstructure interpolates the implicit function expressions of different microstructures to obtain the implicit function expression of the synthetic microstructure, and a corresponding synthetic microstructure geometric model can be established. The update method for the unit filling microstructure during the optimization process combines the distribution of different initial microstructure filling domains within the unit with the synthetic microstructure interpolation method to obtain the microstructure of each unit, and updates are performed in each iteration to clarify the filling microstructure of each unit after each iteration. The geometric configuration of the optimization result of the present invention is clear, without intermediate density units and unclear unit filling microstructures, and has good manufacturability.

Description

Technical Field

[0001] The present invention relates to the technical field of topology optimization, and particularly relates to a topology optimization method and system for precise microstructure filling. Background Art

[0002] Topology optimization is a structural design method that finds the optimal topological configuration within a design domain based on given constraints. With the development of additive manufacturing technology, topology optimization now has more development directions, and multi-material topology optimization and multi-microstructure topology optimization have thus developed rapidly.

[0003] An important problem in multi-material topology optimization and multi-microstructure topology optimization is the material description problem of the overlapping parts of multiple materials. In traditional topology optimization, different material interpolation models are often used to approximately replace the material properties of the overlapping part cells. However, most traditional methods obtain material properties from material numerical values, and when the optimized results are actually manufactured, it will be difficult to manufacture due to the unclear definition of the material structure in the overlapping area. In order to better apply multi-material / multi-microstructure topology optimization in actual engineering, it is necessary to clearly define the material description of the overlapping parts.

[0004] Therefore, in multi-material / multi-microstructure topology optimization, clearly defining the material of the overlapping parts is very helpful for improving the manufacturability of the topology optimization results. Summary of the Invention

[0005] In order to overcome the defects and deficiencies of the existing technology, the present invention provides a topology optimization method and system for precise microstructure filling to solve the problem of unclear overlapping part structures in traditional topology optimization methods when dealing with multi-material / multi-micro-material problems. It includes a synthetic microstructure structure interpolation method based on an initial microstructure and an update method for filling the microstructure of cells during the optimization process. The geometric configuration of its optimization results is clear, without intermediate density cells and unclear cell filling microstructures, and has good manufacturability.

[0006] In order to achieve the above object, the present invention adopts the following technical solutions:

[0007] The present invention provides a topology optimization method for precise microstructure filling, including the following steps:

[0008] Set the basic parameters of topology optimization;

[0009] Set the type and parameters of the triply periodic minimal surface of the initial microstructure;

[0010] Set the initial configuration within the design domain and the grid scale of the design domain;

[0011] Interpolate the implicit function expression of the triply periodic minimal surface of the initial microstructure to obtain the function expression of the synthesized triply periodic minimal surface;

[0012] Construct the corresponding triply periodic minimal surface porous structure using the function expression of the surface;

[0013] Adopt the homogenization method to calculate the macroscopic equivalent mechanical properties of the constructed porous structure, and fit the macroscopic equivalent mechanical properties of different porous structures into a performance function related to the weight value of the triply periodic minimal surface of the initial density in the synthesized surface expression;

[0014] Iterate in a loop to traverse each element in the design domain;

[0015] Obtain the equivalent mechanical properties of the element according to the volume fraction of the microstructure filling domain based on the initial triply periodic minimal surface in the element and the performance function;

[0016] Perform finite element solution;

[0017] Update the component distribution, shape and size in the design domain;

[0018] Judge whether it converges. If it is determined not to converge, return to the iterative loop step. If it is determined to converge, output the final topology optimization data.

[0019] As a preferred technical solution, the basic parameters of the topology optimization are set, specifically including boundary conditions, volume constraints and objective functions.

[0020] As a preferred technical solution, the type and parameters of the triply periodic minimal surface of the initial microstructure are set. One triply periodic minimal surface corresponds to one filling microstructure, and the parameter values adopt the variable parameter values in the Fourier implicit expression of the triply periodic minimal surface.

[0021] As a preferred technical solution, set the initial configuration in the design domain, including the number of components in the design domain, the type, size, shape and position distribution of the filling microstructure of each component.

[0022] As a preferred technical solution, the interpolation of the implicit function expression of the triply periodic minimal surface of the initial microstructure specifically includes the following steps:

[0023]

[0024] Among them, F e represents the Fourier expression of the triply periodic minimal surface filled in the element, e represents the number of the current element, and F i represents the Fourier expression of the initial triply periodic minimal surface, and λ irepresents the weight value of the three - periodic minimal surface of the initial material within the element, \(i\) represents the number of the initial three - periodic minimal surface, and \(n\) represents the number of three - periodic minimal surfaces.

[0025] As a preferred technical solution, the construction of the corresponding three - periodic minimal - surface porous structure using the functional expression of the surface specifically includes the following steps:

[0026] Obtain the corresponding patches by synthesizing the functional expression of the three - periodic minimal surface;

[0027] Offset the patches of the synthesized surface, and the offset amount is the wall thickness of the porous structure;

[0028] Enclose the patches of the initial surface and the offset surface to solidify them.

[0029] As a preferred technical solution, obtaining the equivalent mechanical properties of the element according to the volume fraction of the micro - structure filling domain based on the initial three - periodic minimal surface within the element and the performance function specifically includes the following steps:

[0030] Calculate the volume fraction of the micro - structure filling domain based on the initial three - periodic minimal surface within the element;

[0031] Perform a weighted average on the volume fractions of all micro - structure filling domains in the element, and then take the value obtained by the weighted average of the volume fraction of each micro - structure filling domain as the weight value in the corresponding micro - structure expression within the element, where the weight value is the proportion of the volume fraction of different micro - structure filling domains in the sum of the volume fractions of all micro - structure filling domains within the element;

[0032] Combine the weight value of the initial three - periodic minimal surface within the element with the interpolation method to obtain the Fourier implicit function expression of the three - periodic minimal surface filled within the element. Combine the functional relationship between the weight value and the weight value of the initial structure and the performance of the porous structure to obtain its performance parameters, and substitute the performance parameters into the micro - structure model for calculation to obtain the actual equivalent mechanical properties of the element.

[0033] As a preferred technical solution, the specific calculation formula for the volume fraction of the micro - structure filling domain based on the initial three - periodic minimal surface within the element is:

[0034]

[0035]

[0036]

[0037] Among them, represents the topological description value of the micro - structure filling domain within the design domain, represents the weight value of the micro - structure at the element nodes, represents the volume fraction of the microstructure filling domain based on the initial triply periodic minimal surface within the unit, M represents the number of nodes of the unit, and N represents the number of microstructures occupying the nodes of the unit. represents the area occupied by the microstructure within the design domain, and H represents the Heaviside function;

[0038] The specific calculation formula for the weight value of the initial microstructure within the unit is expressed as:

[0039]

[0040] Among them, represents the weight value of the initial triply periodic minimal surface within the unit, and ρ represents the volume fraction of the microstructure filling domain based on the initial triply periodic minimal surface within the unit;

[0041] The specific calculation formula for the equivalent mechanical properties is expressed as:

[0042]

[0043] Among them, P e represents the actual mechanical property parameter of the unit, represents the mechanical property parameter of the microstructure filled within the unit.

[0044] As a preferred technical solution, the moving asymptote method is used to update the component distribution and shape dimensions within the design domain.

[0045] The present invention also provides a topological optimization system for precise microstructure filling, including: a setting module, an interpolation calculation module, a porous structure construction module, a porous structure equivalent mechanical property calculation module, a performance function fitting module, a cyclic iteration module, a unit equivalent mechanical property calculation module, a finite element solution module, an update module, a convergence judgment module, and an output module;

[0046] The setting module is used to set the basic parameters of the topological optimization, set the type and parameters of the initial triply periodic minimal surface of the initial microstructure, and set the initial configuration within the design domain and the grid scale of the design domain;

[0047] The interpolation calculation module is used to interpolate the implicit function expression of the initial triply periodic minimal surface of the initial microstructure to obtain the function expression of the synthesized triply periodic minimal surface;

[0048] The porous structure construction module is used to construct the corresponding triply periodic minimal surface porous structure by using the function expression of the surface;

[0049] The porous structure equivalent mechanical property calculation module is used to calculate the macroscopic equivalent mechanical properties of the constructed porous structure by using the homogenization method;

[0050] The performance function fitting module is used to fit the macroscopic equivalent mechanical properties of different porous structures into a performance function related to the weight value of the triply periodic minimal surface of the initial density in the synthetic surface expression;

[0051] The cyclic iteration module is used to perform iterative cycles and traverse each triply periodic minimal surface;

[0052] The element equivalent mechanical property calculation module is used to obtain the equivalent mechanical property of the element according to the volume fraction of the microstructure filling domain based on the initial triply periodic minimal surface in the element and the performance function;

[0053] The finite element solving module is used to perform finite element solving;

[0054] The update module is used to update the component distribution and shape dimensions in the design domain;

[0055] The convergence judgment module is used to judge whether convergence occurs. If it is determined that there is no convergence, it returns to the iterative cycle step. If it is determined that convergence occurs, it outputs data;

[0056] The output module is used to output the final topology optimization data.

[0057] Compared with the prior art, the present invention has the following advantages and beneficial effects:

[0058] (1) In each unit within the design domain of the present invention, there is a clearly defined microstructure filling, especially there will be no phenomenon of uncertain microstructure at the overlapping part of multiple materials.

[0059] (2) The transition of the connection between different structures of the present invention is smooth, and there will be no disconnection or stress concentration phenomenon.

[0060] (3) The geometric configuration of the optimization result of the present invention is clear, without intermediate density units and unclear unit filling microstructures, and has good manufacturability. Description of the Drawings

[0061] Figure 1 is a flowchart of the topology optimization method for precise microstructure filling of the present invention;

[0062] Figure 2(a) is a schematic diagram of the initial material distribution within the design domain of the present invention;

[0063] Figure 2(b) is a schematic diagram of the final optimization result of the present invention;

[0064] Figure 3(a) is a schematic diagram of the final optimization result model of the traditional multi-material topology optimization method;

[0065] Figure 3(b) is a schematic diagram of the final optimization result model of the topology optimization method for precise microstructure filling of the present invention. Detailed implementation manners

[0066] In order to make the objectives, technical solutions and advantages of the present invention more clear and understandable, the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0067] Embodiment 1

[0068] As Figure 1 shown, this embodiment provides a topology optimization method for precise microstructure filling, based on the topology optimization framework of the moving deformable component method, including two aspects: the synthetic microstructure interpolation method based on the initial microstructure and the update method of the unit filling microstructure during the optimization process. The synthetic microstructure based on the initial microstructure interpolates the implicit function expressions of different microstructures to obtain the implicit function expression of the synthetic microstructure, and then a corresponding synthetic microstructure geometric model can be established; the update method of the unit filling microstructure during the optimization process is specifically to combine the distribution of different initial microstructure filling domains in the unit with the synthetic microstructure interpolation method to obtain the microstructure of each unit, and update it every iteration to clarify the filling microstructure of each unit after each iteration.

[0069] Specifically, it includes the following steps:

[0070] S1: Set the basic parameters of the topology optimization;

[0071] In this embodiment, the basic parameters include boundary conditions such as the size, load, and constraints of the design domain, and set the volume constraint and the objective function.

[0072] S2: Set the type and parameters of the three-period minimal surface of the initial microstructure;

[0073] In step S2, one type of three-period minimal surface corresponds to one filling microstructure. In this method, the types of microstructures need to be greater than 1, and the parameter values are the variable parameter values in the Fourier implicit expression of the three-period minimal surface. In this embodiment, the primitive and gyroid surfaces are selected as the types of the three-period minimal surfaces of the initial material.

[0074] S3: Set the initial configuration within the design domain and the mesh scale of the design domain;

[0075] The initial configuration in step S3 is the relevant parameters of the components within the design domain, including the number of components within the design domain, the type, size, shape, and position distribution of the filling microstructure of each component.

[0076] As shown in Fig. 2(a), it is the initial configuration of this embodiment. Among them, the part enclosed by the dashed line is the component filled with the microstructure based on the primitive surface, and the part enclosed by the solid line is the component filled with the microstructure based on the gyroid surface. In the initial design of this embodiment, there are a total of 8 components in the design domain, among which 6 components are filled with the microstructure based on the primitive surface, and the remaining ones are filled with the microstructure based on the gyroid surface. All components have the same size and shape in the design domain.

[0077] S4: Interpolate the implicit function expression of the triply periodic minimal surface of the initial microstructure to obtain the function expression of the synthesized triply periodic minimal surface;

[0078] The interpolation calculation method of the Fourier implicit expression of different triply periodic minimal surfaces in step S4 is as follows:

[0079]

[0080] Among them, F e represents the Fourier expression of the triply periodic minimal surface filled in the unit, e represents the number of the current unit, and F i represents the Fourier expression of the initial triply periodic minimal surface, and λ i represents the weight value of the triply periodic minimal surface of the initial material in the unit. At this time, it refers to the ratio of the volume fraction of the microstructure filling domain based on the initial triply periodic minimal surface in the unit to the volume fraction of the solid microstructure filling domain in the unit, and there is i represents the number of the initial triply periodic minimal surface, and n represents the number of triply periodic minimal surfaces.

[0081] S5: Use the function expression of the surface to construct the corresponding triply periodic minimal surface porous structure.

[0082] The construction of the porous structure model in step S5 is essentially based on the initial and the synthesized triply periodic minimal surfaces in step S4 to construct the porous structure. The construction of the porous structure model is realized through the following steps: (1) Obtain the corresponding patches through the function expression of the synthesized triply periodic minimal surface; (2) Offset the patches of the synthesized surface. The offset amounts and offset directions of all involved surfaces need to be consistent, and this offset amount is the wall thickness of the porous structure; (3) Close the patches of the initial surface and the offset surface of the surface and solidify them.

[0083] S6: Use the homogenization method to calculate the macroscopic equivalent mechanical properties of the constructed porous structure, and fit the macroscopic equivalent mechanical properties of different porous structures into a performance function related to the weight value of the initial density triply periodic minimal surface in the synthesized surface expression.

[0084] Step S6: The performance calculation of the porous structure is to calculate the macroscopic performance of a unit cell of the porous structure. In a triply periodic minimal surface, a unit cell is a complete microstructure of the periodic structure body.

[0085] The performance function in step S6 is a function between the mechanical performance of the synthesized surface and the weight value of the initial triply periodic minimal surface in its surface function expression.

[0086] S7: Iterative loop i = 1, 2, … n.

[0087] S8: Obtain the equivalent mechanical performance of the unit according to the volume fraction of different microstructures in the unit and the performance function. The specific steps include:

[0088] (1) The volume fraction of different microstructures in step S8 refers to the volume fraction of the microstructure filling domain based on the initial triply periodic minimal surface in the unit. The specific calculation method is as follows:

[0089]

[0090]

[0091]

[0092] Among them, represents the topological description value of the microstructure filling domain in the design domain, represents the weight value of the microstructure at the unit node, represents the volume fraction of the microstructure filling domain based on the initial triply periodic minimal surface in the unit. M represents the number of nodes of the unit, and N represents the number of microstructures occupying the unit nodes. represents the area occupied by the microstructure in the design domain, and H represents the Heaviside function.

[0093] (2) The specific calculation method of the method for obtaining the unit performance by the distribution of microstructures in the unit in step S8 is to perform a weighted average on the volume fractions of all microstructure filling domains in the unit, and then take the value obtained by the weighted average of the volume fraction of each microstructure filling domain as the weight value in the corresponding microstructure expression in the unit. The weight value is the proportion of the volume fraction of the different microstructure filling domains in the unit in the sum of the volume fractions of all microstructure filling domains in the unit. The weight value of the initial microstructure in the unit is calculated by the following formula:

[0094]

[0095] Among them, represents the weight value of the initial triply periodic minimal surface in the unit, and ρ represents the volume fraction of the microstructure filling domain based on the initial triply periodic minimal surface in the unit.

[0096] (3) The weight value of the initial triply periodic minimal surface within the element, combined with the interpolation method in step S4, yields the Fourier implicit function expression of the specific triply periodic minimal surface filled within the element. The function relationship between the weight value and the weight value of the initial structure obtained in step S6 and the performance of the porous structure can be combined to obtain its performance parameters. Substituting the performance parameters into the microstructure model for calculation, the actual equivalent mechanical performance of the element is obtained and calculated through the following formula:

[0097]

[0098]

[0099] where P e represents the actual mechanical performance parameters of the element, represents the mechanical performance parameters of the microstructure filled within the element.

[0100] S9: Conduct finite element solution.

[0101] S10: Update the component distribution within the design domain.

[0102] The update method in step S10 is to use the moving asymptote method to update the shape size and distribution (spatial position) of the components within the design domain.

[0103] S11: Determine whether to converge. If not, jump to step S7; if convergent, output the data.

[0104] In this embodiment, the convergence criterion can be set according to actual requirements. Available convergence criteria include judging the change amount of the shape and position of all components within the design domain before and after iteration. When the change amount is less than the threshold, it converges.

[0105] The output data in step S11 includes the number, position distribution, and shape size of the components within the design domain.

[0106] As shown in Fig. 2(b), it is the final optimized result corresponding to the initial configuration. Among them, the black part is the overlapping part of the two materials, and there are various microstructures filled inside, all of which are microstructures based on synthetic surfaces. The part enclosed by the dashed line is the component filled with the microstructure based on the primitive surface, and the part enclosed by the solid line is the component filled with the microstructure based on the gyroid surface. As shown in Fig. 3(a) and Fig. 3(b), they are respectively the final optimized result model of the traditional multi-material topology optimization method and the final optimized result model of the topology optimization method with precise microstructure filling. It can be seen that the topology optimization method with precise microstructure filling has a smooth and clear transition between different materials, avoiding the problems of non-smooth connection and non-connection between different materials, and at the same time avoiding the problem of unclear structure in the overlapping part of multiple materials during the manufacturing process, which can better improve the manufacturability of the multi-material topology optimization result.

[0107] Example 2

[0108] This embodiment provides a topology optimization system with precise microstructure filling, including: a setting module, an interpolation calculation module, a porous structure construction module, a porous structure equivalent mechanical property calculation module, a performance function fitting module, a cyclic iteration module, a unit equivalent mechanical property calculation module, a finite element solution module, an update module, a convergence judgment module, and an output module;

[0109] In this embodiment, the setting module is used to set the basic parameters of the topology optimization, set the type and parameters of the three-period minimal surface of the initial microstructure, set the initial configuration within the design domain, and the grid scale of the design domain;

[0110] In this embodiment, the interpolation calculation module is used to interpolate the implicit function expression of the three-period minimal surface of the initial microstructure to obtain the function expression of the synthesized three-period minimal surface;

[0111] In this embodiment, the porous structure construction module is used to construct the corresponding three-period minimal surface porous structure by using the function expression of the surface;

[0112] In this embodiment, the porous structure equivalent mechanical property calculation module is used to calculate the macroscopic equivalent mechanical properties of the constructed porous structure by using the homogenization method;

[0113] In this embodiment, the performance function fitting module is used to fit the macroscopic equivalent mechanical properties of different porous structures into a performance function related to the weight value of the three-period minimal surface of the initial density in the synthesized surface expression;

[0114] In this embodiment, the cyclic iteration module is used to perform iterative loops to traverse each unit within the design domain;

[0115] In this embodiment, the unit equivalent mechanical property calculation module is used to obtain the equivalent mechanical property of the unit according to the volume fraction of the microstructure filling domain based on the initial triply periodic minimal surface in the unit and the performance function;

[0116] In this embodiment, the finite element solution module is used to perform finite element solution;

[0117] In this embodiment, the update module is used to update the component distribution and shape size in the design domain;

[0118] In this embodiment, the convergence judgment module is used to judge whether convergence occurs. If it is determined that there is no convergence, it returns to the iterative loop step. If it is determined that convergence occurs, it outputs data;

[0119] In this embodiment, the output module is used to output the final topology optimization data.

[0120] The above embodiments are preferred embodiments of the present invention, but the embodiments of the present invention are not limited to the above embodiments. Any other changes, modifications, substitutions, combinations, and simplifications made without departing from the spirit and principle of the present invention shall be equivalent replacement methods and are all included in the protection scope of the present invention.

Claims

1. A topological optimization method for precise micro-structure filling, characterized in that It includes the following steps: Set the basic parameters of topology optimization; Set the type and parameters of the triply periodic minimal surface of the initial microstructure; Set the initial configuration within the design domain and the mesh scale of the design domain; Interpolate the implicit function expression of the triply periodic minimal surface of the initial microstructure to obtain the function expression of the synthesized triply periodic minimal surface; Construct the corresponding triply periodic minimal surface porous structure using the function expression of the surface; Adopt the homogenization method to calculate the macroscopic equivalent mechanical properties of the constructed porous structure, and fit the macroscopic equivalent mechanical properties of different porous structures into a performance function related to the weight value of the triply periodic minimal surface of the initial density in the synthesized surface expression; Iterate in a loop to traverse each element within the design domain; Obtain the equivalent mechanical properties of the element according to the volume fraction of the microstructure filling domain based on the initial triply periodic minimal surface within the element and the performance function. The specific steps include: Calculate the volume fraction of the microstructure filling domain based on the initial triply periodic minimal surface within the element; Perform a weighted average of the volume fractions of all microstructure filling domains in the element, and then take the value after the weighted average of the volume fraction of each microstructure filling domain as the weight value in the corresponding microstructure expression within the element, where the weight value is the proportion of the volume fraction of different microstructure filling domains within the element in the sum of the volume fractions of all microstructure filling domains within the element; Combine the weight value of the initial triply periodic minimal surface within the element with the interpolation method to obtain the Fourier implicit function expression of the triply periodic minimal surface filled within the element. Combine the functional relationship between the weight value and the weight value of the initial structure and the performance of the porous structure to obtain its performance parameters, and substitute the performance parameters into the microstructure model for calculation to obtain the actual equivalent mechanical properties of the element; The specific calculation formula for the volume fraction of the microstructure filling domain based on the initial triply periodic minimal surface within the element is: Among them, represents the topological description value of the microstructure filling domain within the design domain, represents the weight value of the microstructure at the unit nodes, represents the volume fraction of the microstructure filling domain based on the initial triply periodic minimal surface within the unit. M represents the number of nodes of the unit, and N represents the number of microstructures occupying the unit nodes, represents the area occupied by the microstructure within the design domain, and H represents the Heaviside function; Perform finite element solution; Update the component distribution and shape dimensions within the design domain; Judge whether it converges. If it is determined not to converge, return to the iterative loop step. If it is determined to converge, output the final topology optimization data.

2. The topology optimization method for precise microstructure filling according to claim 1, wherein The setting of the basic parameters of topology optimization specifically includes boundary conditions, volume constraints, and objective functions.

3. The topology optimization method for precise microstructure filling according to claim 1, characterized in that, The setting of the type and parameters of the triply periodic minimal surface of the initial microstructure. One triply periodic minimal surface corresponds to one filling microstructure, and the parameter values adopt the variable parameter values in the Fourier implicit expression of the triply periodic minimal surface.

4. The topology optimization method for precise microstructure filling according to claim 1, characterized in that Set the initial configuration within the design domain, including the number of components within the design domain, the type, size, shape, and position distribution of the filling microstructure of each component.

5. The topology optimization method for precise microstructure filling according to claim 1, characterized in that The interpolation of the implicit function expression of the triply periodic minimal surface of the initial microstructure specifically includes the following steps: Among them, F e represents the Fourier expression of the triply periodic minimal surface filled in the unit, e represents the number of the current unit, F i represents the Fourier expression of the initial triply periodic minimal surface, λ i represents the weight value of the triply periodic minimal surface of the initial material in the unit, i represents the number of the initial triply periodic minimal surface, and n represents the number of triply periodic minimal surfaces.

6. The topology optimization method for precise microstructure filling according to claim 1, characterized in that The construction of the corresponding triply periodic minimal surface porous structure using the function expression of the surface specifically includes the following steps: Obtain the corresponding patches through the function expression of the synthesized triply periodic minimal surface; Offset the patches of the synthesized surface, and the offset amount is the wall thickness of the porous structure; Close the patches of the initial surface and the offset surface of the surface and solidify them.

7. According to the topology optimization method with precise microstructure filling described in claim 1, wherein The specific calculation formula for the weight value of the initial microstructure within the element is expressed as: Among them, represents the weight value of the initial three-period minimal surface within the element, and e represents the number of the current element; The specific calculation formula for equivalent mechanical properties is expressed as: Among them, P e represents the actual mechanical property parameters of the unit, and represents the mechanical property parameters of the microstructures filled in the unit.

8. The topology optimization method for precise microstructure filling according to claim 1, characterized in that, The moving asymptote method is used to update the component distribution and shape dimensions within the design domain.

9. A topological optimization system for precise microstructure filling, characterized in that, A topology optimization method for realizing the precise microstructure filling described in any one of claims 1-8 above, comprising: a setting module, an interpolation calculation module, a porous structure construction module, a porous structure equivalent mechanical property calculation module, a performance function fitting module, a loop iteration module, an element equivalent mechanical property calculation module, a finite element solution module, an update module, a convergence judgment module, and an output module; The setting module is used to set the basic parameters of topology optimization, set the type and parameters of the three-period minimal surface of the initial microstructure, set the initial configuration within the design domain, and the grid scale of the design domain; The interpolation calculation module is used to interpolate the implicit function expression of the three-period minimal surface of the initial microstructure to obtain the function expression of the synthesized three-period minimal surface; The porous structure construction module is used to construct the corresponding three-period minimal surface porous structure using the function expression of the surface; The porous structure equivalent mechanical property calculation module is used to calculate the macroscopic equivalent mechanical properties of the constructed porous structure using the homogenization method; The performance function fitting module is used to fit the macroscopic equivalent mechanical properties of different porous structures into a performance function related to the weight value of the three-period minimal surface of the initial density in the synthesized surface expression; The loop iteration module is used to perform iterative loops and traverse each three-period minimal surface; The element equivalent mechanical property calculation module is used to obtain the equivalent mechanical property of the element according to the volume fraction of the microstructure filling domain based on the initial three-period minimal surface within the element and the performance function; The finite element solution module is used to perform finite element solution; The update module is used to update the component distribution and shape dimensions within the design domain; The convergence judgment module is used to judge whether convergence occurs. If it is determined that there is no convergence, it returns to the iterative loop step. If it is determined that convergence occurs, it outputs data; The output module is used to output the final topology optimization data.

Images