CN113505880B - 基于特征金字塔网络的拓扑最优结构预测方法 - Google Patents
基于特征金字塔网络的拓扑最优结构预测方法 Download PDFInfo
- Publication number
- CN113505880B CN113505880B CN202110806305.1A CN202110806305A CN113505880B CN 113505880 B CN113505880 B CN 113505880B CN 202110806305 A CN202110806305 A CN 202110806305A CN 113505880 B CN113505880 B CN 113505880B
- Authority
- CN
- China
- Prior art keywords
- design
- loss function
- training data
- neural network
- optimal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000013461 design Methods 0.000 claims abstract description 113
- 238000012549 training Methods 0.000 claims abstract description 59
- 239000000463 material Substances 0.000 claims abstract description 47
- 238000005457 optimization Methods 0.000 claims abstract description 34
- 238000013528 artificial neural network Methods 0.000 claims abstract description 32
- 238000009826 distribution Methods 0.000 claims abstract description 32
- 238000013507 mapping Methods 0.000 claims abstract description 8
- 238000007781 pre-processing Methods 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 41
- 230000008569 process Effects 0.000 claims description 8
- 230000009471 action Effects 0.000 claims description 5
- 238000000265 homogenisation Methods 0.000 claims description 4
- 230000000750 progressive effect Effects 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 15
- 238000006073 displacement reaction Methods 0.000 description 8
- 238000011161 development Methods 0.000 description 3
- 238000003062 neural network model Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005429 filling process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Business, Economics & Management (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Computational Linguistics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Mathematical Physics (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Data Mining & Analysis (AREA)
- Economics (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于特征金字塔网络的拓扑最优结构预测方法，包括：生成训练数据；对训练数据进行预处理；根据预处理后的训练数据构建并训练深度神经网络以拟合在给定约束条件下的设计主体与设计主体的最优材料分布的映射关系，深度神经网络采用特征金字塔网络；构建损失函数，根据损失函数利用反向传播训练更新深度神经网络，损失函数为平均绝对误差损失函数，均方误差损失函数或交叉熵损失函数；利用训练更新后的深度神经网络，预测力学性能指标最优的材料分布。本发明通过利用已有数据训练特征金字塔网络，能够利用训练后的网络模型实现拓扑优化结构的高效快速预测，节约计算时间，降低计算成本。
Description
技术领域
本发明涉及拓扑结构优化技术领域，具体涉及一种基于特征金字塔网络的拓扑最优结构预测方法。
背景技术
随着智能制造的快速发展，对拓扑优化设计提出了缩短设计周期、降低设计成本、保证设计可靠性以及标准化、系列化、通用化等目标，这就要求在给定设计主体、边界条件、负载情况、约束条件和设计目标确定之后，设计一套高效合理的方法对拓扑结构实现省时、降费、提质的优化设计。
拓扑优化结构的设计方案一直是航天领域与机械领域的重要分支，拓扑优化结构的智能化设计也是未来智能制造的重要组成部分。对于拓扑优化结构的设计，需要综合考虑刚度、强度、耐用性、重量比、导热率等多个指标的设计要求，是一个典型的复杂多目标优化设计问题。在传统的拓扑优化设计方法中，主要包括变密度法(SIMP)、均匀化方法(HDM)、渐进结构优化方法(ESO)和水平集方法(LST)。上述几种方法所使用的数学原理是不同的，但都基于梯度算法和有限元计算，它们作为有效的工具在几个领域发挥了极其重要的作用，可以准确地获得不同任务中的最优材料分布。但在拓扑优化设计方法中，直接调用有限元方法来计算刚度、导热率等目标函数，会造成优化时间成本大幅增加，在有限计算资源的条件下甚至无法完成优化求解。并且为了保证设计结构接近最优的结构，大多数时候必须增加设计迭代和变量的数量，而随着优化结构的尺度与维度不断变大，计算时间与计算量将会大幅度提升，导致计算成本增加，研制周期增长。
发明内容
为解决上述现有技术中存在的部分或全部技术问题，本发明提供一种基于特征金字塔网络的拓扑最优结构预测方法。
本发明的技术方案如下：
提供了一种基于特征金字塔网络的拓扑最优结构预测方法，所述方法包括：
生成训练数据，其中，训练数据包括给定的设计主体的约束条件，以及包括在给定约束条件下的设计主体的最优材料分布；
对训练数据进行预处理，以使深度神经网络的输入数据和输出数据的维度相同；
根据预处理后的训练数据构建并训练深度神经网络以拟合在给定约束条件下的设计主体与设计主体的最优材料分布的映射关系，其中，深度神经网络采用特征金字塔网络；
构建损失函数，根据损失函数利用反向传播训练更新深度神经网络，其中，损失函数为平均绝对误差损失函数，均方误差损失函数和交叉熵损失函数中的任一种；
利用训练更新后的深度神经网络，预测力学性能指标最优的材料分布。
在一些可能的实现方式中，所述设计主体的约束条件包括：边界条件，负载条件，体积约束条件和过滤半径约束条件。
在一些可能的实现方式中，所述生成训练数据包括：
对设计主体的设计区域进行网格划分，以使每个网格为一个设计单元，每个设计单元对应一个设计变量，选取设计区域的一个方向为横向，选取与横向垂直的方向为纵向，根据边界条件确定横向边界约束矩阵X1和纵向边界约束矩阵X2，根据负载条件确定横向负载矩阵X3和纵向负载矩阵X4，根据体积约束条件确定体积约束矩阵X5，根据过滤半径约束条件确定过滤半径约束矩阵X6，计算对应的最优材料分布Y1，得到一组训练数据(X1,X2,X3,X4,X5,X6,Y1)；
重复多次上述训练数据的获取过程，直至得到预设数量的训练数据。
在一些可能的实现方式中，所述生成训练数据中，利用传统拓扑优化方法计算最优材料分布，其中，传统拓扑优化方法包括：变密度法、均匀化方法、渐进结构优化方法和水平集方法中的任一种。
在一些可能的实现方式中，设计区域为方形线性设计区域，将方形线性设计区域划分为M1×M2个网格。
在一些可能的实现方式中，所述对训练数据进行预处理中，利用卷积对训练数据(X1,X2,X3,X4,X5,Y1)中的矩阵X1，X2，X3，X4，X5和X6进行处理，以使深度神经网络的输入矩阵X1，X2，X3，X4，X5和X6的维度与深度神经网络的输出矩阵相同。
在一些可能的实现方式中，设计单元对应的设计变量的值设置在[0,1]之间，0表示对应的设计单元内没有材料填充，1表示对应的设计单元内填满材料，0和1之间的中间值表示对应的设计单元内填充有部分材料，且部分材料与填满时的材料的数量比值为当前中间值。
在一些可能的实现方式中，损失函数为均方误差损失函数。
在一些可能的实现方式中，构建均方误差损失函数为：
loss＝MSE
其中，loss表示损失函数，MSE表示材料分布的平均绝对误差。
本发明技术方案的主要优点如下：
本发明的基于特征金字塔网络的拓扑最优结构预测方法通过利用一定数量的训练数据样本训练构建的特征金字塔网络，以学习给定约束条件的设计主体与其拓扑优化结构之间的映射关系，能够利用训练后的特征金字塔网络模型实现拓扑优化结构的高效快速预测，节约计算时间，降低计算成本，且在约束条件发生变化后，仍能够利用训练后的特征金字塔网络模型实现拓扑优化结构的高效快速预测。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明一实施例的基于特征金字塔网络的拓扑最优结构预测方法的流程图；
图2为本发明一实施例的一种设计区域离散示意图；
图3为本发明一实施例的一种设计主体的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明具体实施例及相应的附图对本发明技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本发明的一部分实施例，而不是全部的实施例。基于本发明的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
以下结合附图，详细说明本发明实施例提供的技术方案。
在一般性拓扑优化设计问题中，设计主体分成很多个设计单元，每一个设计单元都是一个设计变量，设计变量的值在[0,1]之间，0表示对应的设计单元内没有材料填充，即密度为0，1表示对应的设计单元内填满材料，即密度为1，0和1之间的中间值表示对应的设计单元内填充有部分材料，且部分材料与填满时的材料的数量比值为当前中间值，即该设计单元为中间密度单元。在一定的区域中需要根据设计填充一定数量材料的设计单元，通过相关力学性质的有限元计算，得出目标函数的大小，使得设计主体的柔度尽可能小。
考虑一个二维平面内的拓扑优化设计问题，其数学模型可以表述为：
其中，C(x)表示优化的目标函数，优化设计的目标是最小化目标函数，x表示设计变量，U表示全局位移矢量，K表示刚度矩阵，xe表示第e个设计变量，ue表示第e个设计单元的单元位移向量，ke表示单元刚度矩阵，p表示惩罚因子，N表示总设计单元数，F表示力矢量，V(x)表示材料体积，V0表示设计区域总体积，Vf表示给定的体积分数。
通过上述公式一可以看出，若要求解该拓扑优化设计问题，必须要实时求解目标方程，得到设计主体的对应的柔度值。该目标函数可以通过有限元的方法进行求解，从而得到准确的目标函数值。但是通常来说，随着设计变量的增多或有限元分析网格的加密，其计算时间也会越来越长，导致计算成本和研制周期增加。
为此，参见图1，本发明一实施例提供了一种基于特征金字塔网络的拓扑最优结构预测方法，该方法包括以下步骤：
步骤100，生成训练数据，其中，训练数据包括给定的设计主体的约束条件，以及包括在给定约束条件下的设计主体的最优材料分布；
步骤200，对训练数据进行预处理，以使深度神经网络的输入数据和输出数据的维度相同；
步骤300，根据预处理后的训练数据构建并训练深度神经网络以拟合在给定约束条件下的设计主体与设计主体的最优材料分布的映射关系，其中，深度神经网络采用特征金字塔网络；
步骤400，构建损失函数，根据损失函数利用反向传播训练更新深度神经网络，其中，损失函数为平均绝对误差损失函数，均方误差损失函数和交叉熵损失函数中的任一种；
步骤500，利用训练更新后的深度神经网络，预测力学性能指标最优的材料分布。
本发明一实施例提供的基于特征金字塔网络的拓扑最优结构预测方法通过利用一定数量的训练数据样本训练构建的特征金字塔网络，以学习给定约束条件的设计主体与其拓扑优化结构之间的映射关系，能够利用训练后的特征金字塔网络模型实现拓扑优化结构的高效快速预测，节约计算时间，降低计算成本，且在约束条件发生变化后，仍能够利用训练后的特征金字塔网络模型实现拓扑优化结构的高效快速预测。
以下对本发明一实施例提供的基于特征金字塔网络的拓扑最优结构预测方法的步骤及原理进行具体说明。
步骤100，生成训练数据。
本发明一实施例中，约束条件包括：边界条件，负载条件，体积约束条件和过滤半径约束条件。
具体地，基于上述设定的设计主体的约束条件，采用以下方式生成训练数据：
对设计主体的设计区域进行网格划分，以使每个网格为一个设计单元，每个设计单元对应一个设计变量，选取设计区域的一个方向为横向，选取与横向垂直的方向为纵向，根据边界条件确定横向边界约束矩阵X1和纵向边界约束矩阵X2，根据负载条件确定横向负载矩阵X3和纵向负载矩阵X4，根据体积约束条件确定体积约束矩阵X5，根据过滤半径约束条件确定过滤半径约束矩阵X6，并计算对应的最优材料分布Y1，得到一组训练数据(X1,X2,X3,X4,X5,X6,Y1)；
重复多次上述训练数据的获取过程，直至得到预设数量的训练数据。
本发明一实施例中，在上述生成训练数据的过程中，可以利用传统拓扑优化方法计算最优材料分布Y1，其中，传统拓扑优化方法包括：变密度法、均匀化方法、渐进结构优化方法和水平集方法中的任一种。
参见图2，设计主体的设计区域可以设计为方形线性设计区域，将方形线性设计区域划分为M1×M2个网格，其中，每个网格刚好是一个设计变量。
参见图3，以图3所示的受力梁为例对步骤100进行示例说明，该受力梁一端横向固定在基础结构上，即有横向约束，且一端的一节点承受力F的作用，另一端的一节点纵向固定在基础结构上，即有纵向约束。
针对上述的受力梁，将受力梁的设计区域设计为方形线性设计区域，对方形线性设计区域进行网格划分，可划分成M1×M2个网格，例如为64×32的网格单元，使得每个网格刚好是一个设计变量xe。在边界条件约束上，设定横向方向的位移约束为ux，设定纵向方向的位移约束为uy，当某节点横向位移有约束，则该节点的ux＝1，若横向位移没有约束，则ux＝0，当某节点纵向位移有约束，则该节点的uy＝1，若纵向位移没有约束，则uy＝0，在上述的受力梁中，由于受力梁左边一端横向全部固支，所以横向边界约束矩阵是左边一列元素都为1，其余元素均为0的(M1+1)×(M2+1)矩阵X1，同理，由于受力梁左边一端没有约束，仅右下角的节点有纵向的约束，所以纵向边界约束矩阵是右下角节点为1，其余元素均为0的(M1+1)×(M2+1)矩阵X2。在负载条件约束上，将承受的力F分解成横向方向px＝Fcosθ，与纵向方向py＝Fsinθ，受力节点值为Fcosθ，其余没有受力的则为0，得到横向负载矩阵是一个受力节点元素为px＝Fcosθ，其余元素均为0的(M1+1)×(M2+1)矩阵X3，同理得到纵向负载矩阵是一个受力节点元素为py＝Fsinθ，其余元素均为0的(M1+1)×(M2+1)矩阵X4。在体积约束条件上，在M1×M2个网格单元中，设定每个xe＝Vf，Vf表示体积分数，生成一个M1×M2的体积约束矩阵X5。在过滤半径约束条件上，在M1×M2个网格单元中，设定每个网格单元中数值r＝Rmin，Rmin表示过滤半径，生成一个M1×M2的过滤半径约束矩阵X6。基于上述确定的约束条件，通过传统拓扑优化方法计算悬臂梁的最优材料分布Y1，获得一组训练数据(X1,X2,X3,X4,X5,X6,Y1)。重复上述过程直至获得足够多的训练数据。
其中，力的作用范围可以在设计主体的任意网格里面，可以根据实际问题的设计要求设定作用范围；体积分数的具体数值范围可以根据实际问题的设计要求进行设定，每个设计单元对应的体积分数可以为设定数值范围中的任一个数值；过滤半径的具体数值范围可以根据实际问题的设计要求进行设定，每个设计单元对应的过滤半径可以为设定数值范围中的任一个数值。
由于训练数据量越多，训练出来的用于预测材料分布的深度神经网络模型的预测精度越高，可选的，本发明一实施例中，可随机生成40000组的训练数据。
步骤200，对训练数据进行预处理，以使深度神经网络的输入数据和输出数据的维度相同。
以上述的受力梁为例，6个维度矩阵(X1,X2,X3,X4,X5,X6)可以看成6个输入到神经网络中的通道，为保证输入通道的一致性，把X1,X2,X3,X4四个通道先做一个卷积核为2×2的卷积，使得四个通道由(M1+1)×(M2+1)矩阵变换为M1×M2矩阵，使得神经网络的六个输入通道与一个输出通道都为M1×M2矩阵，即神经网络的输入矩阵与输出矩阵的维度相同，确保在神经网络中，上采样与下采样的维度变换不容易出错，便于计算。
步骤300，根据预处理后的训练数据构建并训练深度神经网络以拟合在给定约束条件下的设计主体与设计主体的最优材料分布的映射关系。
本发明一实施例中，采用特征金字塔网络(Feature Pyramid Network，FPN)作为深度神经网络模型来拟合在给定约束条件下的设计主体与设计主体的最优材料分布的映射关系。通过利用上述获取的训练数据对特征金字塔网络进行训练和测试，能够使特征金字塔网络学习到训练数据中的物理规律，从而得到一个具有较强泛化能力的神经网络，可以实现对任意给定约束条件下的设计主体的最优材料分布的高精度预测。
从本质上来说，该特征金字塔网络即为一个代理模型，可记为
步骤400，构建损失函数，根据损失函数利用反向传播训练更新深度神经网络。
由于特征金字塔网络模型的输入是经过数据预处理后，六个维度为M1×M2的约束条件对应的矩阵X1,X2,X3,X4,X5,X6，输出是一个维度为M1×M2的预测材料分布可以选择平均绝对误差损失函数，均方误差损失函数和交叉熵损失函数中的任一种损失函数作为损失函数进行反向传播，迭代优化得到最终的特征金字塔网络模型，使得利用得到的特征金字塔网络模型获取的预测材料分布/>可以接近采用传统拓扑优化方法得到的最优材料分布Y1。
可选的，本发明一实施例中，选取均方误差损失函数作为特征金字塔网络模型在训练时的损失函数；通过利用均方误差损失函数在神经网络模型训练过程中进行反向传播，不断迭代优化，进行参数调节，所得到的代理模型相比于采用平均绝对误差损失函数和交叉熵损失函数时的效果更好，预测精度更高。
具体地，构建均方误差损失函数为：
loss＝MSE公式二
其中，loss表示损失函数，MSE表示材料分布的均方误差。
本发明一实施例中，材料分布的均方误差MSE利用以下公式三计算。
其中，M表示总的单元数，表示标签中第i个单元的密度值，/>表示预测中第i个单元的值。
可以看出这是对数据拟合的优化问题。故选取合适的优化算法求解该问题，则可在有限的时间内求出力学性能指标最优的拓扑优化设计方案。
步骤500，利用训练更新后的深度神经网络，预测力学性能指标最优的材料分布。
当完成特征金字塔网络模型的训练更新后，将特征金字塔网络模型作为预测模型，将给定约束条件的设计主体作为模型输入，能够得到对应的设计主体的材料分布，实现在任意给定约束条件的设计主体的拓扑优化结构的高效快速预测。
需要说明的是，在本文中，诸如“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。此外，本文中“前”、“后”、“左”、“右”、“上”、“下”均以附图中表示的放置状态为参照(若存在)。
最后应说明的是：以上实施例仅用于说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (5)
1.一种基于特征金字塔网络的拓扑最优结构预测方法，其特征在于，包括：
生成训练数据，其中，训练数据包括给定的设计主体的约束条件，以及包括在给定约束条件下的设计主体的最优材料分布；
对训练数据进行预处理，以使深度神经网络的输入数据和输出数据的维度相同；
根据预处理后的训练数据构建并训练深度神经网络以拟合在给定约束条件下的设计主体与设计主体的最优材料分布的映射关系，其中，深度神经网络采用特征金字塔网络；
构建损失函数，根据损失函数利用反向传播训练更新深度神经网络，其中，损失函数为平均绝对误差损失函数，均方误差损失函数和交叉熵损失函数中的任一种；
利用训练更新后的深度神经网络，预测力学性能指标最优的材料分布；
所述设计主体的约束条件包括：边界条件，负载条件，体积约束条件和过滤半径约束条件；
所述设计主体为受力梁，所述受力梁一端横向固定在基础结构上，且一端的一节点承受力F的作用，另一端的一节点纵向固定在基础结构上；
所述生成训练数据包括：
将受力梁的设计区域设计为方形线性设计区域，对方形线性设计区域进行网格划分，划分成M1×M2个网格，使得每个网格刚好是一个设计变量xe，在边界条件约束上，获取左边一列元素都为1，其余元素均为0的(M1+1)×(M2+1)横向边界约束矩阵X1，获取右下角节点为1，其余元素均为0的(M1+1)×(M2+1)纵向边界约束矩阵X2，在负载条件约束上，将承受的力F分解成横向方向px＝Fcosθ，与纵向方向py＝Fsinθ，受力节点值为Fcosθ，其余没有受力的则为0，获取一个受力节点元素为px＝Fcosθ，其余元素均为0的(M1+1)×(M2+1)横向负载矩阵X3，获取一个受力节点元素为py＝Fsinθ，其余元素均为0的(M1+1)×(M2+1)纵向负载矩阵X4，在体积约束条件上，在M1×M2个网格单元中，设定每个xe＝Vf，Vf表示体积分数，生成一个M1×M2的体积约束矩阵X5，在过滤半径约束条件上，在M1×M2个网格单元中，设定每个网格单元中数值r＝Rmin，Rmin表示过滤半径，生成一个M1×M2的过滤半径约束矩阵X6，基于上述确定的约束条件，计算受力梁的最优材料分布Y1，获得一组训练数据(X1,X2,X3,X4,X5,X6,Y1)；
重复多次上述训练数据的获取过程，直至得到预设数量的训练数据；
将6个维度矩阵(X1,X2,X3,X4,X5,X6)看成6个输入到深度神经网络中的通道，把X1,X2,X3,X4四个通道先做一个卷积核为2×2的卷积，使得四个通道由(M1+1)×(M2+1)矩阵变换为M1×M2矩阵，使得神经网络的六个输入通道与一个输出通道都为M1×M2矩阵。
2.根据权利要求1所述的基于特征金字塔网络的拓扑最优结构预测方法，其特征在于，所述生成训练数据中，利用传统拓扑优化方法计算最优材料分布，其中，传统拓扑优化方法包括：变密度法、均匀化方法、渐进结构优化方法和水平集方法中的任一种。
3.根据权利要求2所述的基于特征金字塔网络的拓扑最优结构预测方法，其特征在于，设计单元对应的设计变量的值设置在[0,1]之间，0表示对应的设计单元内没有材料填充，1表示对应的设计单元内填满材料，0和1之间的中间值表示对应的设计单元内填充有部分材料，且部分材料与填满时的材料的数量比值为当前中间值。
4.根据权利要求1所述的基于特征金字塔网络的拓扑最优结构预测方法，其特征在于，损失函数为均方误差损失函数。
5.根据权利要求4所述的基于特征金字塔网络的拓扑最优结构预测方法，其特征在于，构建均方误差损失函数为：
loss＝MSE
其中，loss表示损失函数，MSE表示材料分布的均方误差。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110806305.1A CN113505880B (zh) | 2021-07-16 | 2021-07-16 | 基于特征金字塔网络的拓扑最优结构预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110806305.1A CN113505880B (zh) | 2021-07-16 | 2021-07-16 | 基于特征金字塔网络的拓扑最优结构预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113505880A CN113505880A (zh) | 2021-10-15 |
| CN113505880B true CN113505880B (zh) | 2024-04-16 |
Family
ID=78013121
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110806305.1A Active CN113505880B (zh) | 2021-07-16 | 2021-07-16 | 基于特征金字塔网络的拓扑最优结构预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113505880B (zh) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102073761A (zh) * | 2010-12-31 | 2011-05-25 | 锋电能源技术有限公司 | 一种风力发电机组主轴强度计算方法 |
| CN105874510A (zh) * | 2014-01-09 | 2016-08-17 | 西门子产品生命周期管理软件公司 | 用于增材制造的晶格结构的结构保留拓扑优化的方法 |
| CN110414165A (zh) * | 2019-08-01 | 2019-11-05 | 华东交通大学 | 一种基于全局应力约束的多相材料柔顺机构拓扑优化方法 |
| CN111507032A (zh) * | 2020-03-31 | 2020-08-07 | 中国人民解放军军事科学院国防科技创新研究院 | 基于深度学习技术预测温度分布的组件布局优化设计方法 |
| CN111832204A (zh) * | 2020-07-02 | 2020-10-27 | 西安交通大学 | 一种热流耦合结构的非迭代式拓扑优化方法 |
| CN111898730A (zh) * | 2020-06-17 | 2020-11-06 | 西安交通大学 | 一种利用图卷积神经网络结构加速的结构优化设计方法 |
| CN112069623A (zh) * | 2020-09-09 | 2020-12-11 | 大连理工大学 | 一种航空发动机风扇叶片点阵填充方法 |
| CN112100774A (zh) * | 2020-09-16 | 2020-12-18 | 哈尔滨理工大学 | 一种基于变密度法的应力和应变能双约束的拓扑优化方法 |
-
2021
- 2021-07-16 CN CN202110806305.1A patent/CN113505880B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102073761A (zh) * | 2010-12-31 | 2011-05-25 | 锋电能源技术有限公司 | 一种风力发电机组主轴强度计算方法 |
| CN105874510A (zh) * | 2014-01-09 | 2016-08-17 | 西门子产品生命周期管理软件公司 | 用于增材制造的晶格结构的结构保留拓扑优化的方法 |
| CN110414165A (zh) * | 2019-08-01 | 2019-11-05 | 华东交通大学 | 一种基于全局应力约束的多相材料柔顺机构拓扑优化方法 |
| CN111507032A (zh) * | 2020-03-31 | 2020-08-07 | 中国人民解放军军事科学院国防科技创新研究院 | 基于深度学习技术预测温度分布的组件布局优化设计方法 |
| CN111898730A (zh) * | 2020-06-17 | 2020-11-06 | 西安交通大学 | 一种利用图卷积神经网络结构加速的结构优化设计方法 |
| CN111832204A (zh) * | 2020-07-02 | 2020-10-27 | 西安交通大学 | 一种热流耦合结构的非迭代式拓扑优化方法 |
| CN112069623A (zh) * | 2020-09-09 | 2020-12-11 | 大连理工大学 | 一种航空发动机风扇叶片点阵填充方法 |
| CN112100774A (zh) * | 2020-09-16 | 2020-12-18 | 哈尔滨理工大学 | 一种基于变密度法的应力和应变能双约束的拓扑优化方法 |
Non-Patent Citations (2)
| Title |
|---|
| Experimental investigation of strut-and-tie layouts in deep RC beams designed with hybrid bi-linear topology optimization;Jackson L. Jewett 等;《Engineering Structures》;20191231;第1-10页 * |
| 连续体结构拓扑优化敏度过滤研究;张国锋 等;《组合机床与自动化加工技术》;20210630(第6期);第29-32页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113505880A (zh) | 2021-10-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110110413B (zh) | 一种基于材料场缩减级数展开的结构拓扑优化方法 | |
| Jiang et al. | Surrogate-model-based design and optimization | |
| CN112836411B (zh) | 加筋板壳结构的优化方法、装置、计算机设备和存储介质 | |
| CN109344524B (zh) | 一种薄板结构加强筋分布优化方法 | |
| CN110795873B (zh) | 一种考虑尺寸控制的跨尺度拓扑优化方法 | |
| CN111027110B (zh) | 一种连续体结构拓扑与形状尺寸综合优化方法 | |
| US20100058257A1 (en) | Topology optimization method using equivalent static loads | |
| CN112818470B (zh) | 一种基结构的优化方法、装置、计算机设备和存储介质 | |
| CN114741753B (zh) | 一种薄壁加筋结构优化方法、装置、电子设备和存储介质 | |
| Cheng et al. | Interval multiobjective optimization of structures based on radial basis function, interval analysis, and NSGA-II | |
| CN113065268A (zh) | 一种基于代理模型快速多目标空间反射镜优化设计方法 | |
| CN110990944A (zh) | 基于权重比计算的车架多目标拓扑优化方法 | |
| CN113505929B (zh) | 基于嵌入物理约束深度学习技术的拓扑最优结构预测方法 | |
| CN114999591A (zh) | 一种多构型点阵结构的拓扑优化方法 | |
| CN115203997A (zh) | 一种基于多变量设计的点阵-实体复合结构拓扑优化方法 | |
| CN110751729A (zh) | 基于腐蚀-扩散算子的准周期层级结构拓扑优化方法 | |
| CN103065015B (zh) | 一种基于内力路径几何形态的承载结构低碳节材设计方法 | |
| CN113505880B (zh) | 基于特征金字塔网络的拓扑最优结构预测方法 | |
| CN111274624B (zh) | 一种基于rbf代理模型的多工况异形节点拓扑优化设计方法 | |
| CN116362079B (zh) | 一种基于新型插值模型的多材料结构拓扑优化方法 | |
| Zobeiry et al. | Theory-guided machine learning composites processing modelling for manufacturability assessment in preliminary design | |
| CN113139313B (zh) | 基于有限元理论的等效轻量化结构分析方法 | |
| Cunha et al. | Verification of a genetic algorithm for the optimization of stiffened plates through the constructal design method | |
| CN115408914B (zh) | 二维结构的问题无关机器学习拓扑优化方法、介质及产品 | |
| CN116306167B (zh) | 一种基于t样条的复杂机械结构稳健拓扑优化方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |