Topological optimization method and device for anisotropic material structure, equipment and medium

Abstract

The invention discloses a topological optimization method and device for an anisotropic material structure, equipment and a medium, and relates to the technical field of structure topological optimization. The method comprises the steps that a geometric model of a target structure is determined, and a target design area of the geometric model is determined; performing grid division on the target design area, defining design variables and elastic modulus of a target grid in at least two directions, and determining a target function of the target grid; determining constraint conditions of the target grids, and performing iterative optimization on the target functions of the target grids based on the constraint conditions to obtain target density values of the target grids in different directions; and determining a completion instruction in response to the target density values of all the grids of the target design area in different directions, and determining the material distribution of the target design area based on each target density value. According to the scheme, topological optimization can be quickly and accurately carried out on the structure of the anisotropic material, and help is provided for effectively improving the mechanical property of the composite material.

Description

Technical Field

[0001] The present application relates to the technical field of structural topology optimization, and in particular to a topology optimization method, device, equipment and medium for anisotropic material structures. Background Art

[0002] With the continuous development of computer technology, the use of computers to optimize the structure of materials to obtain the optimal material distribution in the design space has been widely used in the industry.

[0003] At present, SIMP (Solid Isotropic Material with Penalization) topology optimization method is mainly used to perform topological optimization of material structure by assuming that there is a fixed relationship between the elastic model of the material and the design variables. However, this method is only applicable to isotropic materials, that is, the materials have the same mechanical properties in all directions. It is difficult to perform structural topological optimization on anisotropic materials such as composite materials or multilayer materials using the existing SIMP topology optimization method.

[0004] How to quickly and accurately perform topological optimization on the structure of anisotropic materials to effectively improve the mechanical properties of composite materials is a key research issue in the industry. Summary of the invention

[0005] The present application provides a topology optimization method, device, equipment and medium for anisotropic material structures, so as to quickly and accurately perform topology optimization on the structures of anisotropic materials, thereby providing assistance for effectively improving the mechanical properties of composite materials.

[0006] According to one aspect of the present application, a topology optimization method for an anisotropic material structure is provided, the method comprising:

[0007] Determining a geometric model of a target structure and determining a target design region of the geometric model; the target structure is composed of anisotropic materials;

[0008] Meshing the target design area, defining design variables and elastic moduli of the target mesh in at least two directions, and determining an objective function of the target mesh based on the design variables and the elastic moduli;

[0009] Determine the constraint conditions of the target grid, and iteratively optimize the objective function of each target grid based on the constraint conditions to obtain target density values ​​of the target grid in different directions;

[0010] In response to an instruction to complete determination of target density values ​​in different directions of all grids of the target design region, material distribution of the target design region is determined based on the target density values.

[0011] According to another aspect of the present application, a topology optimization device for an anisotropic material structure is provided, the device comprising:

[0012] A geometric model determination module, used to determine a geometric model of a target structure and a target design region of the geometric model; the target structure is made of anisotropic material;

[0013] An objective function determination module, used to mesh the target design area, define design variables and elastic moduli of the target mesh in at least two directions, and determine the objective function of the target mesh based on the design variables and the elastic moduli;

[0014] A target density value determination module is used to determine the constraint conditions of the target grid, and iteratively optimize the objective function of each target grid based on each constraint condition to obtain the target density values ​​of the target grid in different directions;

[0015] The material distribution determination module is used to respond to the target density value determination completion instruction of all grids in different directions of the target design area, and determine the material distribution of the target design area based on each target density value.

[0016] According to another aspect of the present application, an electronic device is provided, the electronic device comprising:

[0017] at least one processor; and

[0018] a memory communicatively connected to the at least one processor; wherein,

[0019] The memory stores a computer program that can be executed by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the topology optimization method of the anisotropic material structure described in any embodiment of the present application.

[0020] According to another aspect of the present application, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores computer instructions, and the computer instructions are used to implement the topology optimization method of the anisotropic material structure described in any embodiment of the present application when executed by a processor.

[0021] According to another aspect of the present application, a computer program product is provided, including a computer program, which, when executed by a processor, implements the topology optimization method of the anisotropic material structure described in any embodiment of the present application.

[0022] The technical solution of the embodiment of the present application is as follows: a geometric model of a target structure is determined, and a target design area of ​​the geometric model is determined; the target structure is composed of anisotropic materials; the target design area is meshed, design variables and elastic moduli of the target mesh in at least two directions are defined, and an objective function of the target mesh is determined based on each of the design variables and each of the elastic moduli; constraints of the target mesh are determined, and the objective function of each of the target meshes is iteratively optimized based on each of the constraints to obtain target density values ​​of the target mesh in different directions; a completion instruction is determined in response to target density values ​​of all meshes in different directions of the target design area, and material distribution of the target design area is determined based on each of the target density values, so that topological optimization of the structure of the anisotropic material can be performed quickly and accurately, providing assistance for effectively improving the mechanical properties of composite materials.

[0023] It should be understood that the content described in this section is not intended to identify the key or important features of the embodiments of the present application, nor is it intended to limit the scope of the present application. Other features of the present application will become easily understood through the following description. BRIEF DESCRIPTION OF THE DRAWINGS

[0024] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0025] Figure 1 is a flow chart of a topology optimization method for an anisotropic material structure provided according to Embodiment 1 of the present application;

[0026] Figure 2 is a schematic diagram of a geometric model of a target structure provided according to Embodiment 1 of the present application;

[0027] Figure 3 is a flow chart of a topology optimization method for an anisotropic material structure provided according to Embodiment 2 of the present application;

[0028] Figure 4 is a schematic structural diagram of a topology optimization device for an anisotropic material structure provided in accordance with Embodiment 3 of the present application;

[0029] Figure 5 It is a structural schematic diagram of an electronic device that implements the topology optimization method of the anisotropic material structure of an embodiment of the present application. DETAILED DESCRIPTION

[0030] In order to enable those skilled in the art to better understand the solution of the present application, the technical solution in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of the present application.

[0031] It should be noted that the terms "first", "second", etc. in the specification and claims of the present application and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchangeable where appropriate, so that the embodiments of the present application described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device comprising a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0032] Embodiment 1

[0033] Figure 1 This is a flowchart of a topology optimization method for anisotropic material structures provided in accordance with the first embodiment of the present application. This embodiment is applicable to the case of topology optimization of anisotropic material structures. The method can be performed by a topology optimization device for anisotropic material structures. The topology optimization device for anisotropic material structures can be implemented in the form of hardware and / or software. The topology optimization device for anisotropic material structures can be configured in electronic devices such as computers, servers, or tablet computers. Figure 1 As shown, the method includes:

[0034] Step 110: Determine a geometric model of the target structure and determine a target design region of the geometric model.

[0035] Wherein, the target structure is composed of anisotropic materials; anisotropic materials refer to materials that exhibit different physical or mechanical properties in different directions. Exemplarily, the anisotropic materials involved in this embodiment can be fiber-reinforced composite materials, layered composite materials or single crystal materials, etc., which are not limited in this embodiment; correspondingly, the target structure can be an aircraft shell, a car shell, an engine blade, a pressure vessel, a packaging material or a film, etc., which are not limited in this embodiment.

[0036] Optionally, in this embodiment, after acquiring the target structure, the geometric model of the target structure may be further determined. It is understandable that the geometric model of the target structure may be a three-dimensional model of the target structure, for example, a three-dimensional model of a car engine cover, or a three-dimensional model of an engine blade, etc. Exemplarily, Figure 2 It is a schematic diagram of a geometric model of a target structure provided according to Example 1 of the present application, which shows the three-dimensional structure of the target structure.

[0037] Optionally, in this embodiment, after the geometric model of the target structure is determined, the target design area of ​​the geometric model can be further determined, wherein the target design area can be the entire geometric model or a key part of the geometric model, such as a turning point, a load-bearing part or an intermediate layer of the structure, which is not limited in this embodiment.

[0038] Optionally, in this embodiment, a geometric model of a target structure is determined, and a target design area of ​​the geometric model is determined, including: creating a geometric model of the target structure based on target design software; the geometric model is used to display shape information and size information of the target structure; obtaining design requirements of the target structure, and determining the target design area based on the design requirements, and marking the target design area in the geometric model.

[0039] The target design software may be computer-aided design software, three-dimensional modeling software, or other software that can obtain the three-dimensional structure of the target structure, which is not limited in this embodiment.

[0040] In an optional implementation of the present embodiment, shape information, structural information or size information of the target structure may be input into computer-aided design software, and the computer-aided design software processes the input information to output a geometric model of the target structure; further, the design requirements of the target structure may be obtained. Exemplarily, the design requirements of the target structure may be minimizing the strain energy and mass of the target structure or maximizing the impact resistance of the target structure, which is not limited in the present embodiment.

[0041] In this embodiment, after obtaining the design requirements of the target structure, the target design area can be further determined based on the obtained design requirements, and the target design area can be marked in the geometric model; illustratively, if the design requirements of the target structure are to minimize the mass of the target structure, then the target design area can be the entire geometric model; if the design requirements of the target structure are to maximize the impact resistance of the target structure, then the target design area can be the impacted area in the geometric model, or a buffer material layer (for example, a middle layer of the skin material), etc., which is not limited in this embodiment.

[0042] Step 120 , meshing the target design region, defining design variables and elastic moduli of the target mesh in at least two directions, and determining an objective function of the target mesh based on the design variables and the elastic moduli.

[0043] The target grid may be any grid among a plurality of grids obtained by dividing the target design area, which is not limited in this embodiment.

[0044] In this embodiment, after the target design area is determined in the geometric model, the target design area can be further divided into grids to obtain multiple grids. For example, the target design area can be divided into multiple grids of the same area according to the geometric characteristics of the target design area (for example, area or length, etc.); the target design area can also be divided into multiple grids according to the properties or material distribution of the target design area, which is not limited in this embodiment.

[0045] Furthermore, the design variables and elasticity of the target grid in at least two directions can be defined; illustratively, the initial elastic modulus of the target grid in the fiber direction, the initial density in the fiber direction, the initial elastic modulus in the direction perpendicular to the fiber direction, and the initial density in the direction perpendicular to the fiber direction can be defined; the initial elastic modulus of the target grid in the horizontal direction, the initial density in the horizontal direction, the initial elastic modulus in the vertical direction, the initial density in the vertical direction, the initial elastic modulus in the direction with an angle of 30 degrees to the horizontal direction, and the initial density in the direction with an angle of 30 degrees to the horizontal direction can also be defined; in this embodiment, the value of the initial density in different directions can be 1.

[0046] Furthermore, the objective function of the target grid may be determined based on the defined design variables and elastic modulus of the target grid in at least two directions; in this embodiment, the objective function may be constructed based on the design variables and elastic moduli in different directions.

[0047] Step 130: determine the constraint conditions of the target grid, and iteratively optimize the objective function of each target grid based on the constraint conditions to obtain the target density values ​​of the target grid in different directions.

[0048] Optionally, in this embodiment, the constraint conditions may also be determined according to the design requirements of the target structure; in this embodiment, the constraint conditions may include: load, boundary conditions, volume constraint, stress constraint or displacement constraint, etc.;

[0049] Among them, load refers to the external force or energy input applied to the structure, which can cause structural deformation, stress concentration and possible damage. The target grid may or may not include impact loads; boundary conditions refer to the restrictions on the motion state at the boundary of the target structure; it can be understood that if the target grid is a boundary grid of the non-target design area, then the constraint conditions of the target grid do not include boundary conditions; volume constraints can limit the total volume or mass of the target structure to prevent the optimization result from causing the structure to be too heavy; stress constraints can ensure that the maximum stress of the target structure material during the impact process does not exceed the yield strength of the material to avoid structural damage; displacement constraints can control the deformation of the target structure to avoid exceeding the design allowable range.

[0050] In an optional implementation of this embodiment, the constraint conditions of the target design area may be determined according to the design requirements of the target structure, and further, the constraint conditions of the target grid may be determined according to the grid division situation of the target design area.

[0051] Optionally, in this embodiment, after determining each constraint condition, different variable values ​​can be substituted into the objective function for solving to obtain each output result. Further, it can be determined whether each output result satisfies the above constraints. If the target output result satisfies the above constraints, then the target output result can be determined as the optimal solution of the objective function.

[0052] In another optional implementation of this embodiment, after each constraint condition is determined, iterative calculation can be performed based on each constraint condition and each variable value, and when a set number of iterations is reached, a target output result is output.

[0053] Step 140: In response to an instruction to complete determination of target density values ​​of all grids in different directions of the target design area, determine material distribution of the target design area based on the target density values.

[0054] Optionally, in this embodiment, after determining the target density values ​​of all grids in different directions of the target design area, the material distribution of the target design area can be further determined based on each target density value, and further, the target structure can be manufactured based on the material distribution of the target design area.

[0055] The technical solution of this embodiment is as follows: a geometric model of a target structure is determined, and a target design area of ​​the geometric model is determined; the target structure is composed of anisotropic materials; the target design area is meshed, design variables and elastic moduli of the target mesh in at least two directions are defined, and an objective function of the target mesh is determined based on the design variables and elastic moduli; constraints of the target mesh are determined, and the objective function of each target mesh is iteratively optimized based on the constraints to obtain target density values ​​of the target mesh in different directions; a completion instruction is determined in response to target density values ​​of all meshes in the target design area in different directions, and material distribution of the target design area is determined based on the target density values, so that topological optimization of the structure of anisotropic materials can be performed quickly and accurately, thereby providing assistance for effectively improving the mechanical properties of composite materials.

[0056] Embodiment 2

[0057] Figure 3 This is a flowchart of a topology optimization method for anisotropic material structure provided in Example 2 of the present application. This embodiment is a further refinement of the above technical solution. The technical solution in this embodiment can be combined with various optional solutions in one or more of the above embodiments. Figure 3 As shown, the method includes:

[0058] Step 310: Determine a geometric model of the target structure, and determine a target design region of the geometric model.

[0059] Step 320: input the target geometric model corresponding to the target design area into the target analysis software, and divide the target design area into a plurality of grids; and respectively determine the first design variable and the first elastic modulus of the target grid in the first direction, and the second design variable and the second elastic modulus of the target grid in the second direction.

[0060] Among them, the target analysis software can be computer-aided design software or finite element analysis software, which is not limited in this embodiment; the target grid can be any grid among the grids obtained by dividing the target design area; the first direction can be the horizontal direction of the target grid, or the fiber direction, or any other direction, which is not limited in this embodiment; the second direction can be the vertical direction of the target grid or the direction perpendicular to the fiber direction, or any other direction except the first direction, which is also not limited in this embodiment.

[0061] In an optional implementation of this embodiment, the first direction and the second direction may be determined based on the design requirements of the target structure, for example, each major direction in the design requirements is determined as the first direction, and each minor direction in the design requirements is determined as the second direction. It should be noted that in this embodiment, the number of the first direction and the second direction may be one or more, for example, three first directions and three second directions may be determined.

[0062] Optionally, in this embodiment, after determining the geometric model of the target structure and marking the target design area in the geometric model, the target geometric model corresponding to the target design area can be further input into the target analysis software, thereby dividing the target design area into multiple grids through the target analysis software.

[0063] Furthermore, the first design variable and the first elastic modulus of the target grid in the first direction can be determined; as well as the second design variable and the second elastic modulus of the target grid in the second direction can be determined; wherein the first design variable can be a density value in the first direction, and its initial value can be 1; the second design variable can be a density value in the second direction, and its initial value can be 1.

[0064] In an optional implementation of this embodiment, the first design variable, the first elastic modulus, the second design variable, and the second elastic modulus satisfy the following relationship:

[0065] , ;

[0066] in, is the first elastic modulus; is the first design variable, representing the density of the target grid in the first direction, and its initial value is 1; is the second elastic modulus; is the second design variable, representing the density of the target grid in the second direction, and its initial value is 1; is the penalty factor.

[0067] It should be noted that, in this embodiment, the design variables in each direction not only represent the density of the anisotropic material, but also include the directional information of the anisotropic material, that is, the performance weights of the anisotropic material in different directions.

[0068] The solution of this embodiment defines the relationship between the design variables and the elastic modulus in different directions respectively, which can more accurately reflect the mechanical properties of anisotropic materials and provide a basis for independently adjusting the material properties in different directions during the optimization process.

[0069] Step 330: Determine the objective function of the target grid based on each of the design variables and each of the elastic moduli.

[0070] Optionally, in this embodiment, after the relationship between the design variables and the elastic modulus of the target grid in different directions is defined based on the above steps, the objective function of the target grid can be further determined according to the relationship between the design variables and the elastic modulus in different directions.

[0071] Optionally, in this embodiment, the objective function may be determined based on the following formula:

[0072] ;

[0073] in, is the first design variable, with an initial value of 1; is the penalty factor; is the second design variable, with an initial value of 1; is the elastic modulus of the target mesh in the first direction, is the strain of the target mesh in the first direction; is the elastic modulus of the target mesh in the second direction, is the strain of the target mesh in the second direction, is the volume of the target mesh.

[0074] Correspondingly, in this embodiment, the solution goal that matches the objective function may be minimizing the objective function.

[0075] Step 340: determine the constraint conditions of the target grid, and iteratively optimize the objective function of each target grid based on the constraint conditions to obtain the target density values ​​of the target grid in different directions.

[0076] The constraints of the target mesh may include loads, boundary conditions, volume constraints, stress constraints, and displacement constraints.

[0077] Optionally, in this embodiment, the constraints of the target grid are determined, and the objective function of each target grid is iteratively optimized based on the constraints to obtain target density values ​​of the target grid in different directions. This may include: determining multiple variable values ​​and penalty factor combinations, and substituting each combination into the objective function for solving to obtain multiple density value combinations; the density value combinations include the density values ​​of the target grid in different directions; and when it is determined that the target density value combination satisfies the load, boundary condition, volume constraint, stress constraint or displacement constraint, outputting the target density value.

[0078] In an optional implementation of the present embodiment, different elastic modulus values, strain values, and penalty factor values ​​can be determined, and the different values ​​can be substituted into the above-mentioned objective function respectively, so as to output different density values, for example, multiple density values ​​in the first direction, and multiple density values ​​in the second direction; further, it can be determined whether each density value satisfies the above-mentioned constraints, and the density value in the first direction and the density value in the second direction that satisfy all the constraints are determined as the target density value combination, that is, the target density value finally output.

[0079] In another optional implementation of this embodiment, different elastic modulus values, strain values, and penalty factor values ​​can be determined, and the different values ​​can be substituted into the above-mentioned objective function for iterative calculation. When the iteration stop condition is met, for example, the number of iterations is reached, and the output density values ​​meet the above-mentioned constraints, the target density value can be output.

[0080] Step 350: In response to the target density value determination completion instruction for all meshes in different directions of the target design area, determine the material property of each mesh based on each target density value.

[0081] The material properties may include at least one of the following: elastic modulus, strain, stress, bulk modulus, and Poisson's ratio.

[0082] Optionally, in this embodiment, after determining the density values ​​of each grid in different directions, the material properties such as elastic modulus, strain, stress, bulk modulus or Poisson's ratio in each direction can be further determined based on the density values ​​in each direction, so as to subsequently more accurately determine the material distribution of the target design area.

[0083] Step 360: Determine the material distribution of the target design area.

[0084] Optionally, in this embodiment, after determining the density values ​​and material properties of each grid in the target design area in different directions, the material distribution of the target design area can be further determined based on the density values ​​and material properties of each grid in different directions.

[0085] Step 370: Determine the mechanical properties of the target design area based on the material distribution of the target design area; and perform a manufacturing feasibility assessment on the target structure based on the mechanical properties.

[0086] Optionally, in this embodiment, after determining the material distribution of the target design area, the mechanical properties of the target design area, that is, the mechanical properties of the target structure, can be further determined based on the material distribution of the target design area; further, the manufacturing feasibility of the target structure can be evaluated based on the mechanical properties of the target structure, that is, it can be determined whether the target structure can be manufactured based on the current production environment, and how to manufacture the target structure.

[0087] The solution of this embodiment, after determining the material distribution of the target design area, can further determine the mechanical properties of the target design area based on the material distribution of the target design area; further, the manufacturing feasibility of the target structure can be evaluated based on the mechanical properties of the target structure to determine whether the target structure can be manufactured and how to manufacture the target structure.

[0088] In order to better understand the topology optimization method of the anisotropic material structure involved in this embodiment, a specific example is used below to illustrate it, which mainly includes:

[0089] 1. Geometric model and mesh division

[0090] In this embodiment, a geometric model of the structure may be created and the design area may be meshed.

[0091] 2. Material properties and initial design

[0092] In this embodiment, the elastic modulus of the anisotropic material can be defined and the initial design variables can be set. and is 1.

[0093] 3. Loads and boundary conditions

[0094] In this embodiment, the loads and boundary conditions on the structure can be set.

[0095] 4. Objective function and constraints

[0096] In this embodiment, the objective function can be defined according to the optimization goal, and constraints such as volume and stress can be set.

[0097] 5. Optimize iteration

[0098] In this embodiment, an optimization process may be set in the target design software, and specifically, may include:

[0099] Create a new project: Create a new optimization project in the target design software and add the FEA software as a task component;

[0100] Parameter definition: Define design variables (such as material density in each direction) in the target design software. and ), and set initial values ​​and constraints; design variables can be directly associated with material properties or other model parameters in finite element analysis software.

[0101] In the target design software, define the relationship between material density and elastic modulus by customizing expressions or writing Python scripts; for example, set the elastic modulus to , and the penalty factor p in each iteration can be controlled by the target design software; further, the objective function of the optimization can be defined, such as minimizing the strain energy or mass of the structure, and setting related constraints (e.g., volume constraints, stress limits, etc.).

[0102] Dynamic Adjustment: The target design software can be set to dynamically adjust the design variables (material density) in each iteration and calculate new material property values; the material properties are associated with density through scripts and passed to the FEA software.

[0103] Iterative calculation: Repeat the above steps until the optimization process converges or reaches the preset number of iterations.

[0104] 6. Result analysis and verification

[0105] In this embodiment, the mechanical properties of the optimization results can be verified by finite element analysis, and the manufacturing feasibility can be evaluated.

[0106] The solution of the embodiment of the present application can more accurately reflect the actual performance of anisotropic materials; the optimization result is more in line with the actual working conditions, can effectively improve the mechanical properties of composite materials, and is widely applicable to structural design in the fields of aerospace, automotive industry, civil engineering, etc.

[0107] Embodiment 3

[0108] Figure 4 Schematic diagram of a topology optimization device for anisotropic material structure provided in accordance with the third embodiment of the present application. Figure 4 As shown, the device includes: a geometric model determination module 410, an objective function determination module 420, a target density value determination module 430 and a material distribution determination module 440.

[0109] The geometric model determination module 410 is used to determine the geometric model of the target structure and determine the target design area of ​​the geometric model; the target structure is made of anisotropic material;

[0110] An objective function determination module 420 is used to mesh the target design area, define design variables and elastic moduli of the target mesh in at least two directions, and determine an objective function of the target mesh based on the design variables and the elastic moduli;

[0111] A target density value determination module 430 is used to determine the constraint conditions of the target grid, and iteratively optimize the objective function of each target grid based on each constraint condition to obtain the target density values ​​of the target grid in different directions;

[0112] The material distribution determination module 440 is used to determine the material distribution of the target design area based on each target density value in response to the target density value determination completion instruction of all grids in different directions of the target design area.

[0113] The scheme of this embodiment is to determine the geometric model of the target structure through a geometric model determination module, and determine the target design area of ​​the geometric model; the target structure is composed of anisotropic materials; the target design area is meshed through an objective function determination module, the design variables and elastic moduli of the target mesh in at least two directions are defined, and the objective function of the target mesh is determined based on each of the design variables and each of the elastic moduli; the constraint conditions of the target mesh are determined through a target density value determination module, and the objective function of each of the target meshes is iteratively optimized based on each of the constraint conditions to obtain the target density values ​​of the target mesh in different directions; the material distribution determination module responds to the target density value determination completion instruction of all meshes in the target design area in different directions, and determines the material distribution of the target design area based on each of the target density values, so as to quickly and accurately perform topological optimization on the structure of the anisotropic material, and provide help for effectively improving the mechanical properties of composite materials.

[0114] In an optional implementation of this embodiment, the geometric model determination module 410 is specifically used to create a geometric model of the target structure based on the target design software; the geometric model is used to display the shape information and size information of the target structure;

[0115] The design requirements of the target structure are obtained, and the target design area is determined based on the design requirements, and the target design area is marked in the geometric model.

[0116] In an optional implementation of this embodiment, the target function determination module 420 is specifically used to input the target geometric model corresponding to the target design area into the target analysis software, and divide the target design area into a plurality of grids;

[0117] respectively determining a first design variable and a first elastic modulus of the target grid in a first direction, and a second design variable and a second elastic modulus of the target grid in a second direction;

[0118] The first design variable, the first elastic modulus, the second design variable, and the second elastic modulus satisfy the following relationship:

[0119] , ;

[0120] in, is the first elastic modulus; is the first design variable, representing the density of the target grid in the first direction, and its initial value is 1; is the second elastic modulus; is the second design variable, representing the density of the target grid in the second direction, and its initial value is 1; is the penalty factor.

[0121] In an optional implementation of this embodiment, the objective function determination module 420 is further specifically configured to:

[0122] The objective function is determined based on the following formula:

[0123] ;

[0124] in, is the first design variable, with an initial value of 1; is the penalty factor; is the second design variable, with an initial value of 1; is the elastic modulus of the target mesh in the first direction, is the strain of the target mesh in the first direction; is the elastic modulus of the target mesh in the second direction, is the strain of the target mesh in the second direction, is the volume of the target mesh.

[0125] In an optional implementation of this embodiment, the constraint condition includes at least one of the following: load, boundary condition, volume constraint, stress constraint, and displacement constraint;

[0126] Accordingly, the target density value determination module 430 is specifically used to determine a plurality of variable values ​​and penalty factor combinations, and substitute each of the combinations into the objective function for solving, to obtain a plurality of density value combinations; the density value combinations include the density values ​​of the target grid in different directions;

[0127] When it is determined that the target density value combination satisfies the load, boundary condition, volume constraint, stress constraint or displacement constraint, the target density value is output.

[0128] In an optional implementation of this embodiment, the topology optimization device for anisotropic material structures further includes: a material property determination module, configured to determine the material property of each grid based on each of the target density values;

[0129] The material properties include at least one of the following: elastic modulus, strain, stress, bulk modulus, and Poisson's ratio.

[0130] In an optional implementation of this embodiment, the topology optimization device for anisotropic material structures further includes: a manufacturing feasibility evaluation module for determining the mechanical properties of the target design area based on the material distribution of the target design area;

[0131] The manufacturing feasibility of the target structure is evaluated based on the mechanical properties.

[0132] The topology optimization device for anisotropic material structures provided in the embodiments of the present application can execute the topology optimization method for anisotropic material structures provided in any embodiment of the present application, and has functional modules and beneficial effects corresponding to the execution method.

[0133] Embodiment 4

[0134] Figure 5 A schematic diagram of an electronic device 10 that can be used to implement an embodiment of the present application is shown. The electronic device is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as personal digital processing, cellular phones, smart phones, wearable devices (such as helmets, glasses, watches, etc.) and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples and are not intended to limit the implementation of the present application described and / or required herein.

[0135] like Figure 5As shown, the electronic device 10 includes at least one processor 11, and a memory connected to the at least one processor 11, such as a read-only memory (ROM) 12, a random access memory (RAM) 13, etc., wherein the memory stores a computer program that can be executed by at least one processor, and the processor 11 can perform various appropriate actions and processes according to the computer program stored in the read-only memory (ROM) 12 or the computer program loaded from the storage unit 18 to the random access memory (RAM) 13. In the RAM 13, various programs and data required for the operation of the electronic device 10 can also be stored. The processor 11, the ROM 12, and the RAM 13 are connected to each other through a bus 14. The input / output (I / O) interface 15 is also connected to the bus 14.

[0136] A number of components in the electronic device 10 are connected to the I / O interface 15, including: an input unit 16, such as a keyboard, a mouse, etc.; an output unit 17, such as various types of displays, speakers, etc.; a storage unit 18, such as a disk, an optical disk, etc.; and a communication unit 19, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 19 allows the electronic device 10 to exchange information / data with other devices through a computer network such as the Internet and / or various telecommunication networks.

[0137] The processor 11 may be a variety of general and / or special processing components with processing and computing capabilities. Some examples of the processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various dedicated artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any appropriate processor, controller, microcontroller, etc. The processor 11 executes the various methods and processes described above, such as a topology optimization method for anisotropic material structures, the method comprising: determining a geometric model of a target structure, and determining a target design region of the geometric model; the target structure is composed of anisotropic materials; meshing the target design region, defining design variables and elastic moduli of the target mesh in at least two directions, and determining an objective function of the target mesh based on each of the design variables and each of the elastic moduli; determining constraints of the target mesh, and iteratively optimizing the objective function of each of the target meshes based on each of the constraints to obtain target density values ​​of the target mesh in different directions; determining a completion instruction in response to target density values ​​of all meshes in different directions of the target design region, and determining the material distribution of the target design region based on each of the target density values.

[0138] In some embodiments, the topology optimization method of anisotropic material structure can be implemented as a computer program, which is tangibly contained in a computer-readable storage medium, such as a storage unit 18. In some embodiments, part or all of the computer program can be loaded and / or installed on the electronic device 10 via the ROM 12 and / or the communication unit 19. When the computer program is loaded into the RAM 13 and executed by the processor 11, one or more steps of the topology optimization method of the anisotropic material structure described above can be performed. Alternatively, in other embodiments, the processor 11 can be configured to perform the topology optimization method of the anisotropic material structure in any other appropriate manner (for example, by means of firmware).

[0139] Various implementations of the systems and techniques described above herein can be implemented in digital electronic circuit systems, integrated circuit systems, field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), application specific standard products (ASSPs), systems on chips (SOCs), load programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various implementations can include: being implemented in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which can be a special purpose or general purpose programmable processor that can receive data and instructions from a storage system, at least one input device, and at least one output device, and transmit data and instructions to the storage system, the at least one input device, and the at least one output device.

[0140] The computer programs for implementing the methods of the present application may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, so that when the computer programs are executed by the processor, the functions / operations specified in the flow charts and / or block diagrams are implemented. The computer programs may be executed entirely on the machine, partially on the machine, partially on the machine and partially on a remote machine as a stand-alone software package, or entirely on a remote machine or server.

[0141] In the context of the present application, a computer-readable storage medium may be a tangible medium that may contain or store a computer program for use by or in conjunction with an instruction execution system, device, or equipment. A computer-readable storage medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, devices, or equipment, or any suitable combination of the foregoing. Alternatively, a computer-readable storage medium may be a machine-readable signal medium. A more specific example of a machine-readable storage medium may include an electrical connection based on one or more lines, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0142] To provide interaction with a user, the systems and techniques described herein may be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and a pointing device (e.g., a mouse or trackball) through which the user can provide input to the electronic device. Other types of devices may also be used to provide interaction with the user; for example, the feedback provided to the user may be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user may be received in any form (including acoustic input, voice input, or tactile input).

[0143] The systems and techniques described herein may be implemented in a computing system that includes backend components (e.g., as a data server), or a computing system that includes middleware components (e.g., an application server), or a computing system that includes frontend components (e.g., a user computer with a graphical user interface or a web browser through which a user can interact with implementations of the systems and techniques described herein), or a computing system that includes any combination of such backend components, middleware components, or frontend components. The components of the system may be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: a local area network (LAN), a wide area network (WAN), a blockchain network, and the Internet.

[0144] A computing system may include a client and a server. The client and the server are generally remote from each other and usually interact through a communication network. The client and server relationship is generated by computer programs running on the corresponding computers and having a client-server relationship with each other. The server may be a cloud server, also known as a cloud computing server or cloud host, which is a host product in the cloud computing service system to solve the defects of difficult management and weak business scalability in traditional physical hosts and VPS services.

[0145] It should be understood that the various forms of processes shown above can be used to reorder, add or delete steps. For example, the steps recorded in this application can be executed in parallel, sequentially or in different orders, as long as the expected results of the technical solution of this application can be achieved, and this document is not limited here.

[0146] The above specific implementations do not constitute a limitation on the protection scope of this application. It should be understood by those skilled in the art that various modifications, combinations, sub-combinations and substitutions can be made according to design requirements and other factors. Any modifications, equivalent substitutions and improvements made within the spirit and principles of this application should be included in the protection scope of this application.

Claims

1. A topology optimization method for anisotropic material structure, characterized in that: include: Determining a geometric model of a target structure, and determining a target design region of the geometric model; The target structure is composed of anisotropic material; Meshing the target design area, defining design variables and elastic moduli of the target mesh in at least two directions, and determining an objective function of the target mesh based on the design variables and the elastic moduli; Determine the constraint conditions of the target grid, and iteratively optimize the objective function of each target grid based on the constraint conditions to obtain target density values ​​of the target grid in different directions; In response to an instruction to complete determination of target density values ​​in different directions of all grids of the target design region, material distribution of the target design region is determined based on the target density values.

2. The topology optimization method of anisotropic material structure according to claim 1, characterized in that: The step of determining a geometric model of a target structure and determining a target design region of the geometric model includes: Creating a geometric model of the target structure based on target design software; the geometric model is used to display shape information and size information of the target structure; The design requirements of the target structure are obtained, and the target design area is determined based on the design requirements, and the target design area is marked in the geometric model.

3. The topology optimization method of anisotropic material structure according to claim 1, characterized in that: The step of meshing the target design area and defining design variables and elastic moduli of the target mesh in at least two directions includes: Inputting a target geometric model corresponding to the target design area into target analysis software, and dividing the target design area into a plurality of grids; respectively determining a first design variable and a first elastic modulus of the target grid in a first direction, and a second design variable and a second elastic modulus of the target grid in a second direction; The first design variable, the first elastic modulus, the second design variable, and the second elastic modulus satisfy the following relationship: ， ； in, is the first elastic modulus; is the first design variable, representing the density of the target grid in the first direction, and its initial value is 1; is the second elastic modulus; is the second design variable, representing the density of the target grid in the second direction, and its initial value is 1; is the penalty factor.

4. The topology optimization method of anisotropic material structure according to claim 3, characterized in that: The objective function of determining the target grid based on each of the design variables and each of the elastic moduli includes: The objective function is determined based on the following formula: ； in, is the first design variable, with an initial value of 1; is the penalty factor; is the second design variable, with an initial value of 1; is the elastic modulus of the target mesh in the first direction, is the strain of the target mesh in the first direction; is the elastic modulus of the target mesh in the second direction, is the strain of the target mesh in the second direction, is the volume of the target mesh.

5. The topology optimization method of anisotropic material structure according to claim 4, characterized in that: The constraint condition includes at least one of the following: load, boundary condition, volume constraint, stress constraint and displacement constraint; Accordingly, the determining of the constraint conditions of the target grid and iteratively optimizing the objective function of each target grid based on each constraint condition to obtain the target density values ​​of the target grid in different directions include: Determine multiple variable value and penalty factor combinations, and substitute each of the combinations into the objective function for solving, to obtain multiple density value combinations; the density value combinations include density values ​​of the target grid in different directions; When it is determined that the target density value combination satisfies the load, boundary condition, volume constraint, stress constraint or displacement constraint, the target density value is output.

6. The topology optimization method of anisotropic material structure according to claim 1, characterized in that: Before determining the material distribution of the target design area based on each of the target density values, the method further includes: Determining material properties of each grid based on each of the target density values; The material properties include at least one of the following: elastic modulus, strain, stress, bulk modulus, and Poisson's ratio.

7. The topology optimization method of anisotropic material structure according to claim 1, characterized in that: After determining the material distribution of the target design area based on each of the target density values, the method further includes: Determining the mechanical properties of the target design area based on the material distribution of the target design area; The manufacturing feasibility of the target structure is evaluated based on the mechanical properties.

8. A topology optimization device for anisotropic material structure, characterized in that: include: A geometric model determination module, used to determine a geometric model of a target structure and determine a target design region of the geometric model; The target structure is composed of anisotropic material; An objective function determination module, used to mesh the target design area, define design variables and elastic moduli of the target mesh in at least two directions, and determine the objective function of the target mesh based on the design variables and the elastic moduli; A target density value determination module is used to determine the constraint conditions of the target grid, and iteratively optimize the objective function of each target grid based on each constraint condition to obtain the target density values ​​of the target grid in different directions; The material distribution determination module is used to respond to the target density value determination completion instruction of all grids in different directions of the target design area, and determine the material distribution of the target design area based on each target density value.

9. An electronic device, characterized in that: The electronic device comprises: at least one processor; and a memory communicatively connected to the at least one processor; wherein, The memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can perform the topology optimization method for anisotropic material structure according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the topology optimization method of the anisotropic material structure according to any one of claims 1 to 7 when executed.

Images