Topological optimization design method and device, storage medium and electronic device

Abstract

The invention discloses a topological optimization design method and device, a storage medium and an electronic device.The method comprises the steps that an optimization finite element model corresponding to a target object is established, loads and boundary conditions are applied to the optimization finite element model, and the loads comprise the load corresponding to netting thrust; constructing an optimization column formula corresponding to the optimization finite element model after the load and the boundary condition are applied, and determining a target penalty factor corresponding to a topological optimization algorithm; and solving the optimization column formula according to the target penalty factor and the topological optimization algorithm so as to carry out topological optimization design on the target object. Through the method, the problem that the topological optimization design of the target object under the condition of considering the netting load cannot be effectively processed in the related technology can be solved.

Description

Technical Field

[0001] The present application relates to the field of conductor frame structure design, and more specifically, to a topology optimization design method and device, a storage medium, and an electronic device. Background Art

[0002] The combination of offshore wind power and aquaculture cages is a typical model of integrated development, which can realize additional economic value. The integrated development of offshore wind power and marine ranches is an effective way to improve the marine economy and increase the rate of return per unit of sea area. Unlike ordinary conductor frames, conductor frames with cages installed are subject to greater wave loads, even up to twice the load of ordinary conductor frame wind turbines. Therefore, how to give a conductor frame structure suitable for the combination of aquaculture cages and conductor frames is a difficult problem that needs to be solved.

[0003] There is no effective solution yet for the problem that related technologies cannot effectively handle the topology optimization design of the target object under the consideration of the net load. Summary of the invention

[0004] The embodiments of the present application relate to a topology optimization design method and device, a storage medium and an electronic device. Through the embodiments of the present application, the problem that the related art cannot effectively handle the topology optimization design of the target object under the consideration of the net load can be solved.

[0005] According to one embodiment of the embodiments of the present application, a topology optimization design method is provided, including: establishing an optimized finite element model corresponding to a target object, and applying loads and boundary conditions to the optimized finite element model, wherein the loads include: loads corresponding to the thrust of a mesh; constructing an optimization formula corresponding to the optimized finite element model after applying the loads and the boundary conditions, and determining a target penalty factor corresponding to a topology optimization algorithm; solving the optimization formula according to the target penalty factor and the topology optimization algorithm to perform topology optimization design on the target object.

[0006] In an exemplary embodiment, establishing an optimized finite element model corresponding to a target object includes: establishing a geometric model corresponding to the target object; importing the geometric model into target software so that the target software meshes the geometric model; dividing the geometric model into a plurality of nodes and a plurality of units according to a meshing result, and determining the geometric model including the plurality of nodes and the plurality of units as the optimized finite element model.

[0007] In an exemplary embodiment, applying a load to the optimized finite element model includes: determining a design domain and a non-design domain in the optimized finite element model; determining a first position for applying the load in the optimized finite element model according to the design domain and the non-design domain; performing a force analysis on the target object, and determining load parameters and load type of the load according to the force analysis results and preset design requirements, wherein the load type includes at least one of the following: a line load, a concentrated load and / or a dynamic load; and applying the load to the first position of the optimized finite element model according to the load parameters and the load type.

[0008] In an exemplary embodiment, applying boundary conditions to the optimized finite element model includes: determining a design domain and a non-design domain in the optimized finite element model; determining a second position in the optimized finite element model to apply the boundary condition based on the design domain and the non-design domain; determining a boundary condition type corresponding to the boundary condition, wherein the boundary condition type includes at least one of the following: a fixed constraint, a sliding constraint, a hinge constraint and / or an environmental boundary condition; and applying the boundary condition of the boundary condition type to the second position of the optimized finite actor model.

[0009] In an exemplary embodiment, constructing an optimization formula corresponding to the optimized finite element model after applying the load and the boundary conditions includes: performing a preliminary analysis on the optimized finite element model after applying the load and the boundary conditions to determine a preliminary analysis result, wherein the preliminary analysis result includes: a load path and a force distribution corresponding to the optimized finite element model; determining an objective function according to the preliminary analysis result; and establishing the optimization formula according to the objective function and preset constraints.

[0010] In an exemplary embodiment, determining a target penalty factor corresponding to a topology optimization algorithm includes: obtaining an initial penalty factor; adjusting the initial penalty factor according to a preset design requirement, and determining the adjusted penalty factor as the target penalty factor.

[0011] According to another embodiment of the embodiment of the present application, a topology optimization design device is also provided, including: an establishment module, used to establish an optimized finite element model corresponding to a target object, and apply loads and boundary conditions to the optimized finite element model, wherein the loads include: loads corresponding to the thrust of the net; a construction module, used to construct an optimization column formula corresponding to the optimized finite element model after applying the loads and the boundary conditions, and determine a target penalty factor corresponding to a topology optimization algorithm; a calculation module, used to solve the optimization column formula according to the target penalty factor and the topology optimization algorithm, so as to perform topology optimization design on the target object.

[0012] According to another aspect of the embodiments of the present application, a computer-readable storage medium is provided, in which a computer program is stored, wherein the computer program is configured to execute the above method when running.

[0013] According to another aspect of the embodiments of the present application, an electronic device is provided, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the method through the computer program.

[0014] According to another aspect of the embodiments of the present application, a computer program product is provided, including a computer program, wherein the computer program is executed by a processor to perform the method.

[0015] In an embodiment of the present application, an optimized finite element model corresponding to the target object is established, and loads and boundary conditions including mesh thrust are applied to the optimized finite element model. An optimization formula corresponding to the optimized finite element model after applying the loads and boundary conditions is constructed, and the target penalty factor corresponding to the topology optimization algorithm is determined; the optimization formula is solved according to the target penalty factor and the topology optimization algorithm to perform topology optimization design on the target object. That is to say, after establishing the optimized finite element model, the embodiment of the present application takes into account the effect of mesh thrust when applying loads and boundary conditions to the optimized finite element model, and then constructs an optimization formula corresponding to the optimized finite element model after applying the loads and boundary conditions, and solves the optimization disadvantage to perform topology optimization design on the target object. Through the embodiment of the present application, the problem that the topology optimization design of the target object under the consideration of mesh load cannot be effectively handled in the related art can be solved, thereby realizing the design of the catheter frame under the consideration of mesh load. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] The drawings described herein are used to provide a further understanding of the present application and constitute a part of the present application. The illustrative embodiments of the present application and their descriptions are used to explain the present application and do not constitute an improper limitation on the present application. In the drawings:

[0017] Figure 1 It is a hardware structure block diagram of a computer terminal of a topology optimization design method according to an embodiment of the present application;

[0018] Figure 2 is a flow chart of a topology optimization design method according to an embodiment of the present application;

[0019] Figure 3 It is a schematic diagram of the relationship between the elastic modulus and density variables of the material;

[0020] Figure 4 It is a flow chart of a method for topological optimization design of a pipe cage aquaculture cage according to an optional embodiment of the present application;

[0021] Figure 5 It is a schematic diagram of the optimized finite element model structure in an optional embodiment of the present application;

[0022] Figure 6 is a schematic diagram of a non-design domain according to an optional embodiment of the present application;

[0023] Figure 7 is a schematic diagram of topology optimization results when the thrust of the net clothing in an optional embodiment of the present application is considered as a concentrated load;

[0024] Figure 8 is a schematic diagram of topology optimization results when the thrust of the net is considered as a line load according to an optional embodiment of the present application;

[0025] Fig. 9 It is a structural block diagram of a topology optimization design device according to an embodiment of the present application. DETAILED DESCRIPTION

[0026] 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.

[0027] 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 interchanged 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 that includes a series of steps or units is not necessarily limited to those steps or units that are clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products, or devices; "plurality" means two or more.

[0028] The method embodiments provided in the embodiments of the present application can be executed in a computer terminal or similar computing device or cloud platform or independent physical server or software platform, wherein the above software platform is run through one or more servers. Taking running on a computer terminal as an example, Figure 1 1 is a hardware structure block diagram of a computer terminal of a topology optimization design method according to an embodiment of the present application. Figure 1As shown, the computer terminal may include one or more ( Figure 1 Only one is shown in the figure) processor 102 and memory 104 for storing data. In an exemplary embodiment, the computer terminal may also include a transmission device 106 for communication functions and an input and output device 108, wherein the processor 102 may include but is not limited to a processing device such as a microprocessor MCU or a programmable logic device FPGA. A person skilled in the art can understand that Figure 1 The structure shown is only for illustration and does not limit the structure of the above-mentioned computer terminal. Figure 1 More or fewer components as shown, or with Figure 1 Equivalent functions or comparisons shown Figure 1 A different configuration with more features is shown.

[0029] The memory 104 can be used to store computer programs, for example, software programs and modules of application software, such as computer programs corresponding to the methods in the embodiments of the present application. The processor 102 executes various functional applications and data processing by running the computer programs stored in the memory 104, that is, the above-mentioned methods are implemented. The memory 104 may include a high-speed random access memory, and may also include a non-volatile memory, such as one or more magnetic storage devices, flash memory, or other non-volatile solid-state memory. In some examples, the memory 104 may further include a memory remotely arranged relative to the processor 102, and these remote memories may be connected to the computer terminal via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0030] The transmission device 106 is used to receive or send data via a network. The specific example of the above network may include a wireless network provided by a communication provider of a computer terminal. In one example, the transmission device 106 includes a network adapter (Network Interface Controller, referred to as NIC), which can be connected to other network devices through a base station so as to communicate with the Internet. In one example, the transmission device 106 can be a radio frequency (RF) module, which is used to communicate with the Internet wirelessly.

[0031] In this embodiment, a topology optimization design method is provided, which is applied to the above-mentioned computer terminal. Figure 2 : is a flowchart of a topology optimization design method according to an embodiment of the present application, and the process includes the following steps:

[0032] Step S202, establishing an optimized finite element model corresponding to the target object, and applying loads and boundary conditions to the optimized finite element model, wherein the loads include: loads corresponding to the thrust of the net;

[0033] Wherein, the target object may be: a catheter rack.

[0034] Step S204, constructing an optimization formula corresponding to the optimized finite element model after applying the load and the boundary condition, and determining a target penalty factor corresponding to the topology optimization algorithm;

[0035] Step S206: solving the optimization formula according to the target penalty factor and the topology optimization algorithm to perform topology optimization design on the target object.

[0036] In an embodiment of the present application, an optimized finite element model corresponding to the target object is established, and loads and boundary conditions including mesh thrust are applied to the optimized finite element model. An optimization formula corresponding to the optimized finite element model after applying the loads and boundary conditions is constructed, and the target penalty factor corresponding to the topology optimization algorithm is determined; the optimization formula is solved according to the target penalty factor and the topology optimization algorithm to perform topology optimization design on the target object. That is to say, after establishing the optimized finite element model, the embodiment of the present application takes into account the effect of mesh thrust when applying loads and boundary conditions to the optimized finite element model, and then constructs an optimization formula corresponding to the optimized finite element model after applying the loads and boundary conditions, and solves the optimization disadvantage to perform topology optimization design on the target object. Through the embodiment of the present application, the problem that the topology optimization design of the target object under the consideration of mesh load cannot be effectively handled in the related art can be solved, thereby realizing the design of the catheter frame under the consideration of mesh load.

[0037] Optionally, the above-mentioned step S202 of establishing an optimized finite element model corresponding to the target object includes: establishing a geometric model corresponding to the target object; importing the geometric model into the target software so that the target software meshes the geometric model; dividing the geometric model into a plurality of nodes and a plurality of units according to a result of the meshing, and determining the geometric model including the plurality of nodes and the plurality of units as the optimized finite element model.

[0038] It is understandable that establishing an optimized finite element model corresponding to the target object is a systematic engineering process, involving the creation of a detailed mesh model from the initial geometry concept to the optimization analysis of the model. Specifically:

[0039] (1) Establish a geometric model corresponding to the target object.

[0040] An accurate geometric model of the target object can be created in three-dimensional computer-aided software (CAD software). The geometric model should reflect the shape, size and structural characteristics of the target object in detail. For example, for a pipe cage aquaculture cage, the geometric model needs to include: the pipe rack's columns, beams, connection nodes, as well as the shape of the cage and the connection method between the net and the pipe rack.

[0041] (2) Meshing the geometric model.

[0042] The geometric model is imported into the finite element analysis software. In the finite element analysis software, the geometric model is converted into a mesh model, that is, meshing is performed. Meshing is the discretization of continuous geometric shapes into many small, regular units, each of which can be regarded as a tiny part of the structure.

[0043] The quality of meshing directly affects the accuracy and computational efficiency of finite element analysis. The selection of unit size is a trade-off process. Smaller units can provide more detailed analysis results, but will increase computing time and resource requirements; larger units may simplify the calculation, but at the expense of the accuracy of the results. In the embodiment of the present application, more than 30,000 nodes and units can be selected to ensure the accuracy of the optimized finite element model.

[0044] (3) Divide the geometric model into multiple nodes and multiple units.

[0045] After meshing, the geometric model is discretized into a series of nodes and units. Nodes are connection points in the mesh model, and units are small volumes or areas composed of nodes, such as tetrahedrons, hexahedrons, triangles or quadrilaterals, etc. Nodes and units represent points and volumes in the finite element model, respectively. Each unit is assigned specific properties, such as material properties (elastic modulus, Poisson's ratio, density, etc.) and geometric properties (shape, size). These properties will be used to calculate stress, strain and displacement within the unit.

[0046] Optionally, the above-mentioned step S202 of applying a load to the optimized finite element model includes: determining a design domain and a non-design domain in the optimized finite element model; determining a first position for applying the load in the optimized finite element model according to the design domain and the non-design domain; performing a force analysis on the target object, and determining load parameters and load type of the load according to the force analysis results and preset design requirements, wherein the load type includes at least one of the following: line load, concentrated load and / or dynamic load; and applying the load to the first position of the optimized finite element model according to the load parameters and the load type.

[0047] It can be understood that applying loads to the optimized finite element model can ensure that the results of the topology optimization design can reflect the structural performance under actual working conditions. Specifically:

[0048] (1) Determine the design domain and non-design domain in the optimized finite element model.

[0049] Design domain: the area where material distribution or structural layout can be freely changed during optimization;

[0050] Non-design domain: The area in the structure where the original material distribution and layout must remain unchanged to ensure the critical performance or function of the structure.

[0051] In the embodiment of the present application, the design domain can be: the non-main leg part of the catheter frame (i.e., the target object) or the area where the cage is connected to the catheter frame, and the non-design domain can be: the four main legs of the catheter frame, and the material distribution of these main legs must remain unchanged during the optimization process to ensure the stability of the catheter frame.

[0052] (2) Determine the first position where the load is applied.

[0053] The location where the loads are applied is determined based on the design domain and non-design domain and the actual use scenario of the structure. In the case of the jacket cage, the loads (e.g. wind, wave, net thrust) mainly act within the design domain because the material distribution in these areas is optimized based on the loads. For example, net thrust may act on a specific point or line where the jacket and the cage meet, while the top wind load acts on the key nodes on the top of the jacket.

[0054] (3) Perform stress analysis on the target object and determine the load parameters and load type based on the stress analysis results and preset design requirements.

[0055] Force analysis is to predict the response of the structure under different load conditions by optimizing the finite element model, including stress, strain and displacement. Based on the preset design requirements (for example, structural strength, stability, cost constraints, etc.) and the analysis results, the type of load to be applied (the type of load can be: line load, concentrated load or dynamic load) and specific parameters (load size, direction) can be determined. For example, in the embodiment of the present application, the net thrust may be set as a concentrated load or a line load, depending on the way the net is connected to the conductor frame.

[0056] (4) Applying the load at a first position of the optimized finite element model.

[0057] You can select the load type (e.g., concentrated load or line load), specify the load size and direction, and select the load application point or area (i.e., the first position) in the optimized finite element model. In the embodiment of the present application, a wind load (FX=1N, MY=120N·m) may be applied to the top node, and a net thrust (FX=2.5N or 0.139N) may be applied to the point or line where the net box is connected to the conductor frame.

[0058] Correct setting of loads is the key to ensuring that the optimization analysis results can accurately reflect the structural performance.

[0059] Optionally, the above-mentioned step S202 of applying boundary conditions to the optimized finite element model includes: determining a design domain and a non-design domain in the optimized finite element model; determining a second position for applying the boundary condition in the optimized finite element model according to the design domain and the non-design domain; determining a boundary condition type corresponding to the boundary condition, wherein the boundary condition type includes at least one of the following: a fixed constraint, a sliding constraint, a hinge constraint and / or an environmental boundary condition; and applying the boundary condition of the boundary condition type to the second position of the optimized finite actor model.

[0060] It can be understood that applying boundary conditions to the optimized finite element model can ensure that the optimized finite element model can accurately reflect the constraint state of the structure when simulating real-world loads. Specifically:

[0061] (1) Determine the design domain and non-design domain.

[0062] The determination of the design domain and the non-design domain is consistent with the technical solution for determining the design domain and the non-design domain corresponding to the above-mentioned “applying a load to the optimized finite element model”.

[0063] (2) Determine the second location where the boundary condition is applied.

[0064] The position where the boundary conditions are applied should reflect the constraint state of the structure in the actual working environment. For example, in the embodiment of the present application, all nodes at the bottom of the target object (i.e., the jacket) are fixed, which usually means that under water, the bottom of the jacket is fixed to the seabed or seabed, and no displacement or rotation is allowed. Therefore, these nodes are set to the second position to apply fixed constraint boundary conditions.

[0065] (3) Determine the boundary condition type corresponding to the boundary condition.

[0066] The type of boundary condition depends on the characteristics of the structure and the support conditions. Common types include:

[0067] 1) Fixed constraint: It limits all displacements and rotations of the structure at specific points. The fixed nodes in the model are usually set to this type of constraint to simulate the fixed support state.

[0068] 2) Sliding constraint: allows the structure to slide within a certain plane, but restricts displacement or rotation in other directions. This type of constraint is used when relative sliding is possible between contact surfaces.

[0069] 3) Hinge constraint: allows the structure to rotate at the hinge point, but limits displacement. This constraint is usually used to simulate the working condition of hinge connection.

[0070] 4) Environmental boundary conditions: reflect the impact of the external environment on the structure, such as water flow dynamic pressure, temperature changes, etc. In underwater structure optimization, the static pressure of water can be regarded as an environmental boundary condition.

[0071] In the optimization process of the conductor frame box in the embodiment of the present application, the fixed constraint type is mainly used to cover all the bottom nodes to ensure that the model reflects the actual fixed support state in the optimization iteration.

[0072] (4) Apply boundary conditions to the second position of the optimized finite element model.

[0073] For example, for fixed constraints, you can choose to optimize all nodes at the bottom of the finite element model and impose zero displacement and rotation constraints on these nodes to ensure that they remain fixed during subsequent load application and optimization iterations, simulating actual support conditions.

[0074] Optionally, the above-mentioned step S204 of constructing an optimization formula corresponding to the optimized finite element model after applying the load and the boundary conditions includes: performing a preliminary analysis on the optimized finite element model after applying the load and the boundary conditions to determine a preliminary analysis result, wherein the preliminary analysis result includes: a load path and a force distribution corresponding to the optimized finite element model; determining an objective function according to the preliminary analysis result; and establishing the optimization formula according to the objective function and preset constraints.

[0075] It is understood that the steps of determining the optimized column formula include:

[0076] (1) Perform a preliminary analysis of the optimized finite element model after applying loads and boundary conditions.

[0077] Pre-analysis is the process of making a preliminary assessment of the response of the model under given loads and boundary conditions before formal optimization. Pre-analysis can include calculations of rigid body displacement, stress, strain, and deformation of the target object. The pre-analysis results provide information on the mechanical behavior of the structure (i.e., the target object) under different working conditions. Key information can be the load path (i.e., how the load is transmitted in the structure) and the force distribution (stress and displacement distribution of each part of the structure).

[0078] In the embodiment of the present application, the preliminary analysis simulates how the mesh thrust and top load affect the entire jacket structure and how these forces are distributed in the structure, including possible hot spots (ie, stress concentration areas) and deformation modes.

[0079] (2) Determine the objective function based on the preliminary analysis results.

[0080] The objective function is a mathematical expression of the optimal design, which defines the goal of the optimization, such as minimizing the weight of the structure, maximizing the stiffness of the structure, minimizing the cost, or optimizing the performance of the structure. Based on the preliminary analysis results, the key performance indicators of the structure can be identified and converted into objective functions.

[0081] For example, in the present application embodiment, the objective function may be defined as minimizing the total weight or cost of the structure while ensuring that the structural stiffness meets a specific safety standard. In topology optimization, the objective function usually involves the flexibility of the structure (i.e., the inverse stiffness of the structure) or the total mass.

[0082] (3) Establish an optimized column formula.

[0083] The optimization formula combines the objective function and constraints to form a mathematical problem, which describes how to adjust the design variables (such as the density of the unit) to optimize the objective function while satisfying a series of constraints (such as volume constraints, manufacturing constraints, stability constraints, etc.).

[0084] In the embodiments of the present application, the optimization formula may describe how to change the distribution of materials in the design domain of the jacket under the action of the net thrust and the top load to simultaneously meet the stiffness maximization and volume percentage constraints.

[0085] When constructing the optimization formula, you can use the continuum structural topology optimization method, such as the Solid Isotropic Material with Penalization (SIMP model), which promotes the polarization of material distribution through the penalty factor p (which can be taken as 2.8), making the optimized structure closer to the ideal binary (with or without material) configuration, thereby improving the efficiency and performance of the structure.

[0086] In the SIMP model, it is assumed that the relative cell density x e and the unit elastic modulus E e There is a power function relationship between them, and the mathematical expression is: Wherein, E0: elastic modulus of the entity material; p: penalty factor, generally p≥1, in the embodiment of the present application, p can be taken as: p=2.8.

[0087] In addition, in order to avoid the occurrence of singularity in the stiffness matrix during the calculation process, the relationship between the relative unit density and the unit elastic modulus can be: Among them, E min : Elastic modulus of the hollow part of the material, take E min =10 -9 E0; Figure 3 The relationship between the elastic modulus of the material and the density variable is shown in the formula. The relationship between the elastic modulus of the material (Ee) and the density variable (E0) is as follows: Figure 3 As shown, when x e When they are equal to 0 and 1 respectively, the element represents the hollow part and solid part of the material, and the elastic modulus of the material corresponds to E min and E0. The penalty factor p is assigned to x e Materials between 0 and 1 have lower energy efficiency, making x e The optimization process is polarized, making the SIMP model closer to the discrete variable topology optimization model.

[0088] Optionally, the above step S204 of determining a target penalty factor corresponding to the topology optimization algorithm includes: obtaining an initial penalty factor; adjusting the initial penalty factor according to preset design requirements, and determining the adjusted penalty factor as the target penalty factor.

[0089] It is understood that the steps of determining the target penalty factor include:

[0090] (1) Obtain the initial penalty factor.

[0091] Before starting topology optimization, an initial penalty factor (p) needs to be selected. The penalty factor plays a key role in the above SIMP model, which controls the material density variable (x e ) and unit elastic modulus (E e ). A higher penalty factor pushes the design variable toward 0 (indicating material removal) or 1 (indicating material presence), prompting the optimization process to produce cleaner material distribution results.

[0092] (2) Adjust the initial penalty factor.

[0093] The initial penalty factor may need to be adjusted according to specific design requirements and optimization goals to achieve the best optimization effect. Design requirements may include but are not limited to: structural stability, manufacturing feasibility, cost control, and specific indicators of structural performance. In the embodiment of the present application, it is necessary to ensure that the optimized design can not only withstand the net thrust and top load.

[0094] The purpose of adjusting the penalty factor is to optimize the material distribution so that the design meets the performance requirements while avoiding unnecessary material waste, thereby finding a balance between structural efficiency and cost. In this case, by considering the impact of net thrust as a line load or concentrated load, the goal of the optimization algorithm is to find a jacket design that can meet the structural stiffness requirements while minimizing the use of materials.

[0095] (3) The adjusted penalty factor is determined as the target penalty factor.

[0096] After adjusting and optimizing the initial penalty factor, a penalty factor that can balance the design requirements and performance goals was obtained. This adjusted penalty factor was determined as the target penalty factor (p = 2.8) for subsequent topology optimization calculations. During the optimization iteration process, the target penalty factor will guide the algorithm on how to adjust the density of the unit to maximize the structural stiffness while considering volume constraints and manufacturing feasibility.

[0097] In order to better understand the process of the above-mentioned topology optimization design method, the above-mentioned topology optimization design method is described below in combination with an optional embodiment, but it is not used to limit the technical solution of the embodiment of the present application.

[0098] The optional embodiment of the present application provides a method for topological optimization design of a pipe cage aquaculture cage. The method used in the optional embodiment of the present application includes:

[0099] (1) Continuum structure topology optimization method.

[0100] Assume that each element e in the design domain is assigned a variable relative density variable The value range is [0, 1], and the relative density variable x e There is a one-to-one mapping relationship between the overall variable matrix x and the overall flexibility c of the structure. By solving the minimum value of the overall flexibility c of the continuum, the maximum value of the overall stiffness of the continuum can be obtained, and then the density distribution of each part of the continuum structure can be obtained. Therefore, the topological optimization formula for maximizing the stiffness of the continuum structure is:

[0101] find:x={x1,x2,x3,...,x N}

[0102]

[0103] subject to:V(x) / V0=κ

[0104] KU=F

[0105] 0≤x e ≤1(e=1,2,3,...,N).

[0106] Where, c: compliance; U: overall displacement matrix; K: overall stiffness matrix; E e : elastic modulus of unit e; u e : displacement matrix of unit e; k0: unit stiffness matrix when the density value is 1; V: structure volume; V0: design domain volume; κ: volume percentage.

[0107] (2) SIMP difference.

[0108] (3) Objective function sensitivity.

[0109] The optimization solution based on gradient information needs to derive the sensitivity of the objective function and the constraint function, introduce the Lagrangian operator, and the flexibility is expressed as: c(x) = U T F+λ T (KU-F).

[0110] On both sides of the above equation, e Taking partial derivatives we get:

[0111]

[0112] Since F is the external load and the structural quantity x e No relation, so Therefore, we can get:

[0113]

[0114] Let F T +λ T K=0, we can get: T K=-F T .

[0115] Since K is a symmetric matrix, that is, K = K T , for λ T K=-F T Transposing both sides yields: λ=-U.

[0116] therefore, in,

[0117] Then derive:

[0118] The final result

[0119] (4) Optimization criteria method.

[0120] The common topology optimization formula uses mathematical programming or optimization criterion method OC for optimization solution under the premise of obtaining sensitivity. The mathematical programming method has the advantages of high accuracy and wide application range, but its disadvantage is that the calculation speed is too slow when there are many optimization variables in practice. The optimization criterion method has high convergence iteration efficiency, and the number of iterations will not increase with the increase of the number of design variables. The mathematical formula for OC optimization solution is:

[0121]

[0122] Where, m: motion limit; η: numerical damping coefficient.

[0123] B e Defined as:

[0124] Among them, v e : The unit volume when the density value is 1.

[0125] The Lagrangian operator γ can be determined by bisection, so that the updated volume ratio constraint strictly satisfies the set conditions and the density value of each round of optimization iteration can satisfy the volume constraint conditions.

[0126] Figure 4 Flow chart of a topology optimization design method for a pipe cage aquaculture cage according to an optional embodiment of the present application. Figure 4 As shown:

[0127] Step S401, determining an optimized finite element model.

[0128] To establish a topology optimization finite element model, the number of nodes is generally required to be greater than 30,000 and the number of elements is generally required to be greater than 30,000.

[0129] Figure 5 is a schematic diagram of the optimized finite element model structure in an optional embodiment of the present application, such as Figure 5 As shown, an optional embodiment of the present application may include 31926 nodes and 31681 units.

[0130] Step S402, determining the design domain and the non-design domain.

[0131] Figure 6 is a schematic diagram of a non-design domain according to an optional embodiment of the present application, such as Figure 6 As shown, the optional embodiment of the present application has a jacket height of 52.6m, a top opening of 16m, and a bottom opening of 30m. The four main legs are undesignable areas (red areas), that is, the four main legs must be guaranteed to have material distribution.

[0132] Step S403: applying a load.

[0133] In an optional embodiment of the present application, a rigid element is established and six top loads are applied. Six directional loads are applied to the top of the jacket foundation (the negative sign indicates a downward direction), and topology optimization is performed without considering the non-design domain and with considering the non-design domain.

[0134] Among them, top load: FX=1N, MY=120N.

[0135] Net load: The total horizontal thrust of the net is 10N. The net load can be considered as a four-point loading method, and the average thrust at each point is FX=2.5N.

[0136] Step S404: applying boundary conditions.

[0137] Apply boundary conditions to fix all nodes at the bottom.

[0138] Step S405, determining the optimized column formula.

[0139] Step S406, determining the interpolation penalty factor (ie, penalty factor) p.

[0140] Generally speaking, p≥1. In an optional embodiment of the present application, the value may be: p=2.8.

[0141] Step S407, obtaining optimization results.

[0142] Figure 7 1 is a schematic diagram of the topology optimization result of considering the thrust of the net as a concentrated load in an optional embodiment of the present application. Figure 7 As shown:

[0143] In the working condition setting, the net thrust (acting on four points), bending moment and horizontal thrust are considered.

[0144] Setup requirements:

[0145] The walls around the conductor frame are the design area; the effect of the net load is considered; the conductor frame height range is -30 to 22.6m, and the net range is -10 to 4m; the horizontal thrust of the net and the top horizontal thrust are 10:1, and the net load can consider the four-corner point / two-line loading method; the ratio of the top horizontal thrust to the bending moment is 1:120; the stiffness maximization model under the volume ratio constraint; set rotational symmetry + plane symmetry to ensure that the material layout of the four faces is exactly the same; set the minimum size to eliminate the checkerboard phenomenon.

[0146] Figure 8 is a schematic diagram of the topology optimization result of considering the thrust of the net as a line load according to an optional embodiment of the present application, such as Figure 8 As shown: the net load is distributed as a line load (two lines), and the size is evenly divided. The thrust of the net line load is evenly distributed on 72 nodes, and FX on each node is 0.139N.

[0147] The results show that the structure contains two layers of X-shaped cross bars but they are all distributed from the top of the net to the bottom of the catheter frame. More materials can be arranged at this position, so the structure can be strengthened.

[0148] Through the description of the above implementation methods, those skilled in the art can clearly understand that the method according to the above embodiment can be implemented by means of software plus a necessary general hardware platform, and of course by hardware, but in many cases the former is a better implementation method. Based on such an understanding, the technical solution of the present application, or the part that contributes to the prior art, can be embodied in the form of a software product, which is stored in a storage medium (such as ROM / RAM, magnetic disk, optical disk), and includes a number of instructions for a terminal device (which can be a mobile phone, computer, server, or network device, etc.) to execute the methods of each embodiment of the present application.

[0149] The present application also provides a structural block diagram of a topology optimization design device. Fig. 9 is a structural block diagram of a topology optimization design device according to an embodiment of the present application; Fig. 9 As shown, including:

[0150] Establishing module 92, used to establish an optimized finite element model corresponding to the target object, and apply loads and boundary conditions to the optimized finite element model, wherein the loads include: loads corresponding to the thrust of the net;

[0151] A construction module 94 is used to construct an optimization formula corresponding to the optimized finite element model after applying the load and the boundary condition, and determine a target penalty factor corresponding to the topology optimization algorithm;

[0152] The calculation module 96 is used to solve the optimization formula according to the target penalty factor and the topology optimization algorithm to perform topology optimization design on the target object.

[0153] Through the above-mentioned device, an optimized finite element model corresponding to the target object is established, and loads and boundary conditions including mesh thrust are applied to the optimized finite element model. The optimization formula corresponding to the optimized finite element model after the loads and boundary conditions are applied is constructed, and the target penalty factor corresponding to the topology optimization algorithm is determined; the optimization formula is solved according to the target penalty factor and the topology optimization algorithm to perform topology optimization design on the target object. That is to say, after the optimized finite element model is established, when the loads and boundary conditions are applied to the optimized finite element model, the effect of the mesh thrust is taken into account in the embodiment of the present application, and then the optimization formula corresponding to the optimized finite element model after the loads and boundary conditions are applied is constructed, and the optimization disadvantage is solved to perform topology optimization design on the target object. Through the embodiment of the present application, the problem that the topology optimization design of the target object under the mesh load cannot be effectively handled in the related technology can be solved, thereby realizing the design of the catheter frame under the mesh load.

[0154] In an exemplary embodiment, the establishment module 92 is also used to establish a geometric model corresponding to the target object; import the geometric model into the target software so that the target software meshes the geometric model; divide the geometric model into multiple nodes and multiple units according to the mesh division result, and determine the geometric model including the multiple nodes and the multiple units as the optimized finite element model.

[0155] In an exemplary embodiment, the establishment module 92 is also used to determine the design domain and the non-design domain in the optimized finite element model; determine the first position for applying the load in the optimized finite element model according to the design domain and the non-design domain; perform a force analysis on the target object, and determine the load parameters and load type of the load according to the force analysis results and preset design requirements, wherein the load type includes at least one of the following: line load, concentrated load and / or dynamic load; apply the load to the first position of the optimized finite element model according to the load parameters and the load type.

[0156] In an exemplary embodiment, the establishment module 92 is also used to determine the design domain and the non-design domain in the optimized finite element model; determine a second position in the optimized finite element model to apply the boundary condition based on the design domain and the non-design domain; determine the boundary condition type corresponding to the boundary condition, wherein the boundary condition type includes at least one of the following: fixed constraints, sliding constraints, hinged constraints and / or environmental boundary conditions; and apply the boundary condition of the boundary condition type to the second position of the optimized finite actor model.

[0157] In an exemplary embodiment, the construction module 94 is also used to perform a preliminary analysis on the optimized finite element model after applying the load and the boundary conditions to determine the preliminary analysis results, wherein the preliminary analysis results include: the load path and force distribution corresponding to the optimized finite element model; determining the objective function according to the preliminary analysis results; and establishing the optimization formula according to the objective function and preset constraints.

[0158] In an exemplary embodiment, the construction module 94 is further configured to obtain an initial penalty factor; adjust the initial penalty factor according to a preset design requirement, and determine the adjusted penalty factor as the target penalty factor.

[0159] An embodiment of the present application further provides a storage medium, which includes a stored program, wherein the program executes any of the above methods when it is run.

[0160] Optionally, in this embodiment, the storage medium may be configured to store program codes for executing the following steps:

[0161] S1, establishing an optimized finite element model corresponding to the target object, and applying loads and boundary conditions to the optimized finite element model, wherein the loads include: loads corresponding to the thrust of the net clothing;

[0162] S2, constructing an optimization column formula corresponding to the optimized finite element model after applying the load and the boundary condition, and determining a target penalty factor corresponding to the topology optimization algorithm;

[0163] S3, solving the optimization formula according to the target penalty factor and the topology optimization algorithm to perform topology optimization design on the target object.

[0164] An embodiment of the present application further provides an electronic device, including a memory and a processor, wherein a computer program is stored in the memory, and the processor is configured to run the computer program to execute the steps in any one of the above method embodiments.

[0165] Optionally, the electronic device may further include a transmission device and an input / output device, wherein the transmission device is connected to the processor, and the input / output device is connected to the processor.

[0166] Optionally, in this embodiment, the processor may be configured to perform the following steps through a computer program:

[0167] S1, establishing an optimized finite element model corresponding to the target object, and applying loads and boundary conditions to the optimized finite element model, wherein the loads include: loads corresponding to the thrust of the net clothing;

[0168] S2, constructing an optimization column formula corresponding to the optimized finite element model after applying the load and the boundary condition, and determining a target penalty factor corresponding to the topology optimization algorithm;

[0169] S3, solving the optimization formula according to the target penalty factor and the topology optimization algorithm to perform topology optimization design on the target object.

[0170] An embodiment of the present application further provides a computer program product, including a computer program, which is used by a processor to execute the steps in any of the above method embodiments.

[0171] Optionally, in this embodiment, the computer program product may be executed by a processor to perform the following steps:

[0172] S1, establishing an optimized finite element model corresponding to the target object, and applying loads and boundary conditions to the optimized finite element model, wherein the loads include: loads corresponding to the thrust of the net clothing;

[0173] S2, constructing an optimization column formula corresponding to the optimized finite element model after applying the load and the boundary condition, and determining a target penalty factor corresponding to the topology optimization algorithm;

[0174] S3, solving the optimization formula according to the target penalty factor and the topology optimization algorithm to perform topology optimization design on the target object.

[0175] Optionally, in this embodiment, the above-mentioned storage medium may include but is not limited to: a USB flash drive, a read-only memory (ROM), a random access memory (RAM), a mobile hard disk, a magnetic disk or an optical disk, and other media that can store program codes.

[0176] Optionally, the specific examples in this embodiment may refer to the examples described in the above embodiments and optional implementation modes, and this embodiment will not be described in detail here.

[0177] Obviously, those skilled in the art should understand that the above modules or steps of the present application can be implemented by a general computing device, they can be concentrated on a single computing device, or distributed on a network composed of multiple computing devices, and optionally, they can be implemented by a program code executable by a computing device, so that they can be stored in a storage device and executed by the computing device, and in some cases, the steps shown or described can be executed in a different order from that herein, or they can be made into individual integrated circuit modules, or multiple modules or steps therein can be made into a single integrated circuit module for implementation. Thus, the present application is not limited to any specific combination of hardware and software.

[0178] The above description is only the preferred embodiment of the present application and is not intended to limit the present application. For those skilled in the art, the present application may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the principles of the present application shall be included in the protection scope of the present application.

Claims

1. A topology optimization design method, characterized in that: include: Establishing an optimized finite element model corresponding to the target object, and applying loads and boundary conditions to the optimized finite element model, wherein the loads include: loads corresponding to the thrust of the net; Constructing an optimization formula corresponding to the optimized finite element model after applying the load and the boundary conditions, and determining a target penalty factor corresponding to the topology optimization algorithm; The optimization formula is solved according to the target penalty factor and the topology optimization algorithm to perform topology optimization design on the target object.

2. The topology optimization design method according to claim 1, characterized in that: Establish an optimized finite element model corresponding to the target object, including: Establishing a geometric model corresponding to the target object; Importing the geometric model into the target software so that the target software performs mesh division on the geometric model; The geometric model is divided into a plurality of nodes and a plurality of units according to a division result of the mesh division, and the geometric model including the plurality of nodes and the plurality of units is determined as the optimized finite element model.

3. The topology optimization design method according to claim 1, characterized in that: Applying a load to the optimized finite element model includes: Determining a design domain and a non-design domain in the optimized finite element model; determining a first position for applying the load in the optimized finite element model according to the design domain and the non-design domain; Performing a force analysis on the target object, and determining the load parameters and load type of the load according to the force analysis results and preset design requirements, wherein the load type includes at least one of the following: linear load, concentrated load and / or dynamic load; The load is applied at a first location of the optimized finite element model according to the load parameter and the load type.

4. The topology optimization design method according to claim 1, characterized in that: Applying boundary conditions to the optimized finite element model includes: Determining a design domain and a non-design domain in the optimized finite element model; determining a second position for applying the boundary condition in the optimized finite element model according to the design domain and the non-design domain; Determine a boundary condition type corresponding to the boundary condition, wherein the boundary condition type includes at least one of the following: a fixed constraint, a sliding constraint, a hinge constraint and / or an environmental boundary condition; and apply the boundary condition of the boundary condition type to a second position of the optimized finite actor model.

5. The topology optimization design method according to claim 1, characterized in that: Constructing an optimization formula corresponding to the optimized finite element model after applying the load and the boundary condition, including: Performing a pre-analysis on the optimized finite element model after applying the load and the boundary condition to determine a pre-analysis result, wherein the pre-analysis result includes: a load path and a force distribution corresponding to the optimized finite element model; Determine the objective function according to the pre-analysis result; The optimization formula is established according to the objective function and preset constraints.

6. The topology optimization design method according to claim 1, characterized in that: Determine the target penalty factor corresponding to the topology optimization algorithm, including: Get the initial penalty factor; The initial penalty factor is adjusted according to preset design requirements, and the adjusted penalty factor is determined as the target penalty factor.

7. A topology optimization design device, characterized in that: include: An establishment module is used to establish an optimized finite element model corresponding to the target object, and apply loads and boundary conditions to the optimized finite element model, wherein the loads include: loads corresponding to the thrust of the net clothing; A construction module is used to construct an optimization formula corresponding to the optimized finite element model after applying the load and the boundary condition, and determine a target penalty factor corresponding to the topology optimization algorithm; A calculation module is used to solve the optimization formula according to the target penalty factor and the topology optimization algorithm to perform topology optimization design on the target object.

8. A computer-readable storage medium, characterized in that: The computer-readable storage medium includes a stored program, wherein the program executes the method described in any one of claims 1 to 6 when executed.

9. An electronic device comprising a memory and a processor, characterized in that: A computer program is stored in the memory, and the processor is configured to execute the method according to any one of claims 1 to 6 through the computer program.

10. A computer program product, comprising a computer program, characterized in that The computer program is used by a processor to execute the steps of the method described in any one of claims 1 to 6.

Images