[go: up one dir, main page]

CN114564878A - 一种基于图神经网络的结构动力学模型生成方法 - Google Patents

一种基于图神经网络的结构动力学模型生成方法 Download PDF

Info

Publication number
CN114564878A
CN114564878A CN202111674213.9A CN202111674213A CN114564878A CN 114564878 A CN114564878 A CN 114564878A CN 202111674213 A CN202111674213 A CN 202111674213A CN 114564878 A CN114564878 A CN 114564878A
Authority
CN
China
Prior art keywords
matrix
structural
neural network
generating
data
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.)
Granted
Application number
CN202111674213.9A
Other languages
English (en)
Other versions
CN114564878B (zh
Inventor
王迎
张俊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology Shenzhen
Original Assignee
Harbin Institute of Technology Shenzhen
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology Shenzhen filed Critical Harbin Institute of Technology Shenzhen
Priority to CN202111674213.9A priority Critical patent/CN114564878B/zh
Publication of CN114564878A publication Critical patent/CN114564878A/zh
Application granted granted Critical
Publication of CN114564878B publication Critical patent/CN114564878B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • 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
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/08—Learning methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Software Systems (AREA)
  • Artificial Intelligence (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Mathematical Physics (AREA)
  • Computational Linguistics (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Data Mining & Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于图神经网络的结构动力学模型的生成方法，所述方法包括：确定目标结构，获取所述目标结构的结构响应数据和拓扑信息；根据所述结构响应数据和所述拓扑信息对预设的人工神经网络进行训练；将训练后的所述人工神经网络作为所述目标结构对应的结构动力学模型。应用本发明生成的结构动力学模型可以计算(预测)结构的动力学响应。本发明通过将目标结构的拓扑信息引入深度学习网络的框架之中，可以在不失去深度学习网络对于目标结构的物理特征的提取能力的同时，一定程度上提高模型的泛化能力，实现人工神经网络模型在不同拓扑结构中的转换。解决了现有技术中用于模拟结构动力学的深度学习网络泛化能力差的问题。

Description

一种基于图神经网络的结构动力学模型生成方法
技术领域
本发明涉及深度学习领域，尤其涉及的是一种基于图神经网络的结构动力学模型生成方法。
背景技术
传统的结构动力学模拟方法(如有限元法)的目标是基于数学模型建立的，然而复杂结构的数学模型是很难获得的，因此传统方法是不能反映结构真实的物理特性的。为了克服传统模拟方法的局限性，机器学习方法，如SVM[1]，和深度学习方法，如ANN[2,3]、CNN[4]、RNN[5]、LSTM[6,7]等已经被引入到结构动力学的模拟过程中。但是，人工智能模型需要大量的实测数据用来训练，且训练好的模型泛化能力往往很差。因此，基于人工智能算法的结构动力学模拟方法也具有其自身的局限性。
因此，现有技术还有待改进和发展。
发明内容
本发明要解决的技术问题在于，针对现有技术的上述缺陷，提供一种基于图神经网络的结构动力学模型生成方法，旨在解决现有技术中用于模拟结构动力学的深度学习网络泛化能力差的问题。
本发明解决问题所采用的技术方案如下：
第一方面，本发明实施例提供一种基于图神经网络的结构动力学模型生成方法，其中，所述方法包括：
确定目标结构，获取所述目标结构的结构响应数据和拓扑信息；
根据所述结构响应数据和所述拓扑信息对预设的人工神经网络进行训练；
将训练后的所述人工神经网络作为所述目标结构对应的结构动力学模型。
在一种实施方式中，所述获取所述目标结构的结构响应数据，包括：
根据预设的抽样时间点采集目标结构上若干自由度分别对应的位移响应数据、速度响应数据以及加速度响应数据；
将所述位移响应数据、所述速度响应数据以及所述加速度响应数据作为所述结构响应数据。
在一种实施方式中，所述根据所述结构响应数据和所述拓扑信息对预设的人工神经网络进行训练，包括：
根据所述拓扑信息生成发送矩阵和接收矩阵，其中，所述发送矩阵和所述接收矩阵的组合可以反映所述目标结构的若干自由度之间的相互作用拓扑信息；
根据若干所述自由度分别对应的所述位移响应数据生成结构位移向量，根据若干所述自由度分别对应的所述速度响应数据生成结构速度向量，根据若干所述自由度分别对应的所述加速度响应数据生成训练标签数据；
根据所述发送矩阵、所述接收矩阵、所述结构位移向量以及所述结构速度向量和所述训练标签数据对所述人工神经网络进行训练。
在一种实施方式中，所述人工神经网络包括第一多层感知机、第二多层感知机、第三多层感知机、第四多层感知机，所述根据所述发送矩阵、所述接收矩阵、所述结构位移向量以及所述结构速度向量和所述训练标签数据对所述人工神经网络进行训练，包括：
根据所述发送矩阵、所述接收矩阵以及所述结构位移向量生成恢复力计算矩阵，将所述恢复力计算矩阵输入所述第一多层感知机，得到恢复力预测矩阵；
根据所述发送矩阵、所述接收矩阵以及所述结构速度向量生成阻尼力计算矩阵，将所述阻尼力计算矩阵输入所述第二多层感知机，得到阻尼力预测矩阵；
根据所述恢复力预测矩阵和所述阻尼力预测矩阵得到合力的计算矩阵，将所述合力的计算矩阵输入所述第三多层感知机，得到合力预测矩阵，其中，所述合力预测矩阵用于反映恢复力与阻尼力相加的矢量和；
获取外部激励荷载数据，根据所述接收矩阵、所述合力预测矩阵以及所述外部激励荷载数据生成加速度响应计算矩阵，将所述加速度响应计算矩阵输入所述第四多层感知机，得到加速度响应预测数据；
根据所述加速度响应预测数据和所述训练标签数据对所述人工神经网络进行训练。
在一种实施方式中，所述根据所述发送矩阵、所述接收矩阵以及所述结构位移向量生成恢复力计算矩阵，包括：
将所述发送矩阵与所述结构位移向量相乘，得到第一矩阵；
将所述接收矩阵与所述结构位移向量相乘，得到第二矩阵；
将所述第一矩阵与所述第二矩阵按列合并，得到所述恢复力计算矩阵。
在一种实施方式中，所述根据所述发送矩阵、所述接收矩阵以及所述结构速度向量生成阻尼力计算矩阵，包括：
将所述发送矩阵与所述结构速度向量相乘，得到第三矩阵；
将所述接收矩阵与所述结构速度向量相乘，得到第四矩阵；
将所述第三矩阵与所述第四矩阵按列合并，得到所述阻尼力计算矩阵。
在一种实施方式中，所述根据所述恢复力预测矩阵和所述阻尼力预测矩阵得到合力的计算矩阵，包括：
将所述恢复力预测矩阵和所述阻尼力预测矩阵按列合并，得到所述合力的计算矩阵。
在一种实施方式中，所述根据所述接收矩阵、所述合力预测矩阵以及所述外部激励荷载数据生成加速度响应计算矩阵，包括：
将所述合力预测矩阵左乘以所述接收矩阵的转置，得到第五矩阵；
将所述第五矩阵与所述外部激励荷载数据按列合并，得到所述加速度响应计算矩阵。
第二方面，本发明实施例还提供一种基于图神经网络的结构动力学模型生成系统，其中，所述系统包括：
数据获取模块，用于确定目标结构，获取所述目标结构的结构响应数据和拓扑信息；
网络训练模块，用于根据所述结构响应数据和所述拓扑信息对预设的人工神经网络进行训练，将训练后的所述人工神经网络作为所述目标结构对应的结构动力学模型。
第三方面，本发明实施例还提供一种基于图神经网络的结构动力学模型，其中，所述基于图神经网络的结构动力学模型采用上述任意所述的基于图神经网络的结构动力学模型生成方法生成。
本发明的有益效果：本发明实施例通过确定目标结构，获取所述目标结构的结构响应数据和拓扑信息；根据所述结构响应数据和所述拓扑信息对预设的人工神经网络进行训练；将训练后的所述人工神经网络作为所述目标结构对应的结构动力学模型。应用本发明生成的结构动力学模型可以计算(预测)结构的动力学响应。本发明通过将目标结构的拓扑信息引入深度学习网络的框架之中，可以在不失去深度学习网络对于目标结构的物理特征的提取能力的同时，一定程度上提高模型的泛化能力，实现人工神经网络模型在不同拓扑结构中的转换。解决了现有技术中用于模拟结构动力学的深度学习网络泛化能力差的问题。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的基于图神经网络的结构动力学模型生成方法的流程示意图。
图2是本发明实施例提供的弹簧-质量块系统的示意图。
图3是本发明实施例提供的人工神经网络的训练过程示意图。
图4是本发明实施例提供的地震加速度与时间的关系示意图。
图5是本发明实施例提供的三自由度系统中第一个自由度的预测结果示意图。
图6是本发明实施例提供的的三自由度系统中第二个自由度的预测结果示意图。
图7是本发明实施例提供的的三自由度系统中第三个自由度的预测结果示意图。
图8是本发明实施例提供的的四自由度系统中第一个自由度的预测结果示意图。
图9是本发明实施例提供的四自由度系统中第二个自由度的预测结果示意图。
图10是本发明实施例提供的四自由度系统中第三个自由度的预测结果示意图。
图11是本发明实施例提供的四自由度系统中第四个自由度的预测结果示意图。
图12是本发明实施例提供的基于图神经网络的结构动力学模型生成系统的内部模块示意图。
图13是本发明实施例提供的终端的原理框图。
具体实施方式
本发明公开了一种基于图神经网络的结构动力学模型生成方法，为使本发明的目的、技术方案及效果更加清楚、明确，以下参照附图并举实施例对本发明进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
本技术领域技术人员可以理解，除非特意声明，这里使用的单数形式“一”、“一个”、“所述”和“该”也可包括复数形式。应该进一步理解的是，本发明的说明书中使用的措辞“包括”是指存在所述特征、整数、步骤、操作、元件和/或组件，但是并不排除存在或添加一个或多个其他特征、整数、步骤、操作、元件、组件和/或它们的组。应该理解，当我们称元件被“连接”或“耦接”到另一元件时，它可以直接连接或耦接到其他元件，或者也可以存在中间元件。此外，这里使用的“连接”或“耦接”可以包括无线连接或无线耦接。这里使用的措辞“和/或”包括一个或更多个相关联的列出项的全部或任一单元和全部组合。
本技术领域技术人员可以理解，除非另外定义，这里使用的所有术语(包括技术术语和科学术语)，具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是，诸如通用字典中定义的那些术语，应该被理解为具有与现有技术的上下文中的意义一致的意义，并且除非像这里一样被特定定义，否则不会用理想化或过于正式的含义来解释。
传统的结构动力学模拟方法(如有限元法)的目标是基于数学模型建立的，然而复杂结构的数学模型是很难获得的，因此传统方法是不能反映结构真实的物理特性的。为了克服传统模拟方法的局限性，机器学习方法，如SVM[1]，和深度学习方法，如ANN[2,3]、CNN[4]、RNN[5]、LSTM[6,7]等已经被引入到结构动力学的模拟过程中。但是，人工智能模型需要大量的实测数据用来训练，且训练好的模型泛化能力往往很差。因此，基于人工智能算法的结构动力学模拟方法也具有其自身的局限性。
针对现有技术的上述缺陷，本发明提供一种基于图神经网络的结构动力学模型生成方法，所述方法包括：确定目标结构，获取所述目标结构的结构响应数据和拓扑信息；根据所述结构响应数据和所述拓扑信息对预设的人工神经网络进行训练；将训练后的所述人工神经网络作为所述目标结构对应的结构动力学模型。本发明通过将目标结构的拓扑信息引入深度学习网络的训练过程。应用本发明生成的结构动力学模型可以计算(预测)结构的动力学响应。本发明通过将目标结构的拓扑信息引入深度学习网络的框架之中，可以在不失去深度学习网络对于目标结构的物理特征的提取能力的同时，一定程度上提高模型的泛化能力，实现人工神经网络模型在不同拓扑结构中的转换。解决了现有技术中用于模拟结构动力学的深度学习网络泛化能力差的问题。
如图1所示，所述方法包括如下步骤：
步骤S100、确定目标结构，获取所述目标结构的结构响应数据和拓扑信息。
具体地，本实施例中的目标结构可以为任意一个需要进行结构动力学模拟过程的复杂结构。为了生成目标结构对应的结构动力学模型，本实施例需要获取目标结构的结构响应数据和拓扑信息，其中，结构响应数据可以反映目标结构各节点在外部激励下的响应，拓扑信息可以反映目标结构的物理特征。
在一种实现方式中，所述步骤S100具体包括：
步骤S101、根据预设的抽样时间点采集目标结构上若干自由度分别对应的位移响应数据、速度响应数据以及加速度响应数据；
步骤S102、将所述位移响应数据、所述速度响应数据以及所述加速度响应数据作为所述结构响应数据。
具体地，为了得到目标结构的结构响应数据，本实施例需要在预先设定的抽样时间点采集目标结构在各个自由度上的位移、速度以及加速度响应，得到位移响应数据、速度响应数据以及加速度响应数据。然后将这三种数据作为目标结构的结构响应数据。
例如，以图2所示的弹簧-质量块系统为例，该模型属于一维动力学系统，即每个节点只有一个平动自由度，常被应用于计算考虑层间剪切的建筑结构，是建筑结构抗震设计计算中重要的理论基础。图2所示结构的结构响应数据可通过下述公式(1)计算获得，其中，
Figure BDA0003450393440000081
Y分别为各自由度的加速度响应、速度响应和位移响应；k、c、m、f分别为刚度、阻尼、质量和外部激励；i∈[1,ndof]表示自由度编号。
Figure BDA0003450393440000082
其中，
Figure BDA0003450393440000091
Figure BDA0003450393440000092
Figure BDA0003450393440000093
公式(1)还可以写成公式(2)的形式：
Figure BDA0003450393440000094
如图1所示，所述方法还包括：
步骤S200、根据所述结构响应数据和所述拓扑信息对预设的人工神经网络进行训练。
具体地，为了提高人工神经网络的泛化能力，本实施例需要将结构的拓扑信息输入到人工神经网络的框架之中。由于拓扑信息可以反映目标结构的物理特征，因此经过结构响应数据和拓扑信息训练后得到的人工神经网络就可以在不失去深度学习网络的特征提取能力的同时，提高模型的泛化能力，实现人工神经网络模型在不同拓扑结构中的转换。
在一种实现方式中，所述步骤S200具体包括：
步骤S201、根据所述拓扑信息生成发送矩阵和接收矩阵，其中，所述发送矩阵和所述接收矩阵的组合可以反映所述目标结构的若干自由度之间的相互作用拓扑信息；
步骤S202、根据若干所述自由度分别对应的所述位移响应数据生成结构位移向量，根据若干所述自由度分别对应的所述速度响应数据生成结构速度向量，根据若干所述自由度分别对应的所述加速度响应数据生成训练标签数据；
步骤S203、根据所述发送矩阵、所述接收矩阵、所述结构位移向量以及所述结构速度向量和所述训练标签数据对所述人工神经网络进行训练。
具体地，本实施例中基于结构拓扑信息生成一对关系矩阵，即发送矩阵和接收矩阵，通过修改这一对关系矩阵的内部元素可以用来表示所有复杂结构的相互作用拓扑信息，它们组合在一起表示的物理意义与有限元刚度矩阵中元素位置的物理意义相似。然后根据各自由度的位移响应数据生成结构位移向量、速度响应数据生成结构速度向量、加速度响应数据生成训练标签数据(用于训练和测试使用)。最后根据发送矩阵、接收矩阵、结构位移向量以及结构速度向量生成人工神经网络的训练数据，并采用训练标签数据评估人工神经网络输出的结果与真实结果之间差距，进而以两者的差距为导向对人工神经网络中的参数进行迭代更新，直至人工神经网络输出的结果与真实结果之间差距达到预设的训练目标。
举例说明，如图2所示的三自由度系统的动力响应可由下述公式(3)获得，其中的·代表向量按位相乘，若两个向量之间没有符号则表示正常的矩阵相乘。其中，Msr与Mrr组成一对关系矩阵，可以用来表示所有复杂结构的相互作用拓扑信息，Msr称为发送矩阵，Mrr称为接收矩阵，两个矩阵组合在一起表示的物理意义与有限元刚度矩阵中元素位置的物理意义相似。
Figure BDA0003450393440000111
Figure BDA0003450393440000112
其中，等号右边由三部分组成，分别表示动力学系统的：恢复力、阻尼力和外部激励。
在一种实现方式中，所述人工神经网络包括第一多层感知机、第二多层感知机、第三多层感知机、第四多层感知机，所述步骤S203，具体包括：
步骤S2031、根据所述发送矩阵、所述接收矩阵以及所述结构位移向量生成恢复力计算矩阵，将所述恢复力计算矩阵输入所述第一多层感知机，得到恢复力预测矩阵；
步骤S2032、根据所述发送矩阵、所述接收矩阵以及所述结构速度向量生成阻尼力计算矩阵，将所述阻尼力计算矩阵输入所述第二多层感知机，得到阻尼力预测矩阵；
步骤S2033、根据所述恢复力预测矩阵和所述阻尼力预测矩阵得到合力的计算矩阵，将所述合力的计算矩阵输入所述第三多层感知机，得到合力预测矩阵，其中，所述合力预测矩阵用于反映恢复力与阻尼力相加的矢量和；
步骤S2034、获取外部激励荷载数据，根据所述接收矩阵、所述合力预测矩阵以及所述外部激励荷载数据生成加速度响应计算矩阵，将所述加速度响应计算矩阵输入所述第四多层感知机，得到加速度响应预测数据；
步骤S2035、根据所述加速度响应预测数据和所述训练标签数据对所述人工神经网络进行训练。
需要理解的是，如上述公式(3)所示，即使外部激励fi是已知的，但由于结构的刚度和阻尼是无法测量得到的，因此恢复力和阻尼力未知。除此之外，传统的弹簧-质量块系统假设每个节点的质量是由结构的总质量按照比例分配到每个节点上的，这种假设缺乏依据。因此本实施例中通过四个多层感知机(MLP)来间接求得以上三个未知参数。其中，第一多层感知机用于预测恢复力，第二多层感知机用于预测阻尼力，第三多层感知机用于预测恢复力和阻尼力的合力，第四多层感知机用于预测加速度响应。具体地，由于发送矩阵和接收矩阵是基于目标结构的拓扑信息生成的，可以反映目标结构的物理特征，因此根据发送矩阵、接收矩阵以及结构位移向量生成第一多层感知机的输入数据，即恢复力计算矩阵，第一多层感知机基于输入的数据会输出一个恢复力预测矩阵。同理，根据发送矩阵、接收矩阵以及结构速度向量生成第二多层感知机的输入数据，即阻尼力计算矩阵，第二多层感知机基于输入的数据会输出一个阻尼力预测矩阵。然后，基于预测出的恢复力预测矩阵和阻尼力预测矩阵生成第三多层感知机的输入数据，即合力的计算矩阵，第三多层感知机基于输入的数据会输出一个合力预测矩阵。之后，基于已知的外部激励荷载数据、接收矩阵以及合力预测矩阵生成第四多层感知机的输入数据，即加速度响应计算矩阵，第四多层感知机基于输入的数据会输出一个加速度响应预测数据，并将该加速度响应预测数据作为人工神经网络最终的输出数据，由于训练标签数据为加速度响应真实数据，因此通过比较加速度响应预测数据和训练标签数据可以计算出人工神经网路的损失值，该损失值可以反映人工神经网络输出的预测数据和真实数据之间的差距，进而以该损失值为导向对人工神经网络的网络参数进行更新，直至该损失值收敛至训练目标。
在一种实现方式中，所述步骤S2031，具体包括：
步骤S20311、将所述发送矩阵与所述结构位移向量相乘，得到第一矩阵；
步骤S20312、将所述接收矩阵与所述结构位移向量相乘，得到第二矩阵；
步骤S20313、将所述第一矩阵与所述第二矩阵按列合并，得到所述恢复力计算矩阵。
具体地，将发送矩阵和接收矩阵分别乘以结构位移向量，得到第一矩阵和第二矩阵，再将两个矩阵按列合并为一个2列的矩阵，即得到恢复力计算矩阵，之后将其作为第一多层感知机的输入数据(如图3所示)。
在一种实现方式中，所述步骤S2032，具体包括：
步骤S20321、将所述发送矩阵与所述结构速度向量相乘，得到第三矩阵；
步骤S20322、将所述接收矩阵与所述结构速度向量相乘，得到第四矩阵；
步骤S20323、将所述第三矩阵与所述第四矩阵按列合并，得到所述阻尼力计算矩阵。
具体地，将发送矩阵和接收矩阵分别乘以结构速度向量，得到第三矩阵和第四矩阵，再将两个矩阵按列合并为一个2列的矩阵，即得到阻尼力计算矩阵，之后其作为第二多层感知机的输入数据(如图3所示)。
在一种实现方式中，所述步骤S2033，具体包括：
步骤S20331、将所述恢复力预测矩阵和所述阻尼力预测矩阵按列合并，得到所述合力的计算矩阵。
具体地，将获得的恢复力预测矩阵和阻尼力预测矩阵按列合并为一个2列的矩阵，即得到合力的计算矩阵，之后将其作为第三多层感知机的输入数据(如图3所示)。
在一种实现方式中，所述步骤S2034，具体包括：
步骤S20341、将所述合力预测矩阵左乘以所述接收矩阵的转置，得到第五矩阵；
步骤S20342、将所述第五矩阵与所述外部激励荷载数据按列合并，得到所述加速度响应计算矩阵。
具体地，将合力预测矩阵左乘接收矩阵Mrr的转置，得到第五矩阵，再将第五矩阵与外部激励荷载数据按列合并为一个2列的矩阵，即得到加速度响应计算矩阵，之后再将其作为第四多层感知机的输入数据(如图3所示)。
为了证明本发明的技术效果，发明人做了以下实验：
以图2所示的三自由度系统作为训练模型，根据公式(1)生成数据，并训练神经网络，然后应用训练好的神经网络模型预测四自由度系统的加速度响应(设置系统各自由度的初始速度和位移分别为0m/s2和0m，只需通过输入各抽样时间点对应的外部激励数值，就可预测结构的时域响应加速度)。三自由度系统和四自由度系统的参数如表1和表2所示。选择一组地震加速度(BAJA_020787_251)作为外部激励，该地震加速度的波动形式如图4所示。为了计算时间序列，需要对加速度响应进行数值积分来求解节点自由度的速度和位移，数值积分方法如式(4)所示。该算例采用的抽样频率为1000Hz,训练结果和测试结果分别如图5、6、7和图8、9、10、11所示。由结果图可以看出，该方法能够在结构参数未知的情况下，准确的预测结构的加速度响应。且由于拓扑信息的输入，使得神经网络模型的泛化能力得到了大大的提升。
表1.三自由度弹簧-质量块系统参数
Figure BDA0003450393440000151
表2.四自由度弹簧-质量块系统参数
Figure BDA0003450393440000152
Figure BDA0003450393440000153
Figure BDA0003450393440000154
基于上述实施例，本发明还提供了一种基于图神经网络的结构动力学模型生成系统，如图12所示，所述系统包括：
数据获取模块01，用于确定目标结构，获取所述目标结构的结构响应数据和拓扑信息；
网络训练模块02，用于根据所述结构响应数据和所述拓扑信息对预设的人工神经网络进行训练，将训练后的所述人工神经网络作为所述目标结构对应的结构动力学模型。
基于上述实施例，本发明还提供了一种终端，其原理框图可以如图13所示。该终端包括通过系统总线连接的处理器、存储器、网络接口、显示屏。其中，该终端的处理器用于提供计算和控制能力。该终端的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该终端的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现基于图神经网络的结构动力学模型生成方法。该终端的显示屏可以是液晶显示屏或者电子墨水显示屏。
本领域技术人员可以理解，图13中示出的原理框图，仅仅是与本发明方案相关的部分结构的框图，并不构成对本发明方案所应用于其上的终端的限定，具体的终端可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
在一种实现方式中，所述终端的存储器中存储有一个或者一个以上的程序，且经配置以由一个或者一个以上处理器执行所述一个或者一个以上程序包含用于进行基于图神经网络的结构动力学模型生成方法的指令。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本发明所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
综上所述，本发明公开了一种基于图神经网络的结构动力学模型生成方法，所述方法包括：确定目标结构，获取所述目标结构的结构响应数据和拓扑信息；根据所述结构响应数据和所述拓扑信息对预设的人工神经网络进行训练；将训练后的所述人工神经网络作为所述目标结构对应的结构动力学模型。本发明通过将目标结构的拓扑信息引入深度学习网络的框架之中。应用本发明生成的结构动力学模型可以计算(预测)结构的动力学响应。本发明通过将目标结构的拓扑信息引入深度学习网络的框架之中，可以在不失去深度学习网络对于目标结构的物理特征的提取能力的同时，一定程度上提高模型的泛化能力，实现人工神经网络模型在不同拓扑结构中的转换。解决了现有技术中用于模拟结构动力学的深度学习网络泛化能力差的问题。
应当理解的是，本发明的应用不限于上述的举例，对本领域普通技术人员来说，可以根据上述说明加以改进或变换，所有这些改进和变换都应属于本发明所附权利要求的保护范围。

Claims (10)

1.一种基于图神经网络的结构动力学模型生成方法，其特征在于，所述方法包括：
确定目标结构，获取所述目标结构的结构响应数据和拓扑信息；
根据所述结构响应数据和所述拓扑信息对预设的人工神经网络进行训练；
将训练后的所述人工神经网络作为所述目标结构对应的结构动力学模型。
2.根据权利要求1所述的基于图神经网络的结构动力学模型生成方法，其特征在于，所述获取所述目标结构的结构响应数据，包括：
根据预设的抽样时间点采集目标结构上若干自由度分别对应的位移响应数据、速度响应数据以及加速度响应数据；
将所述位移响应数据、所述速度响应数据以及所述加速度响应数据作为所述结构响应数据。
3.根据权利要求2所述的基于图神经网络的结构动力学模型生成方法，其特征在于，所述根据所述结构响应数据和所述拓扑信息对预设的人工神经网络进行训练，包括：
根据所述拓扑信息生成发送矩阵和接收矩阵，其中，所述发送矩阵和所述接收矩阵的组合可以反映所述目标结构的若干自由度之间的相互作用拓扑信息；
根据若干所述自由度分别对应的所述位移响应数据生成结构位移向量，根据若干所述自由度分别对应的所述速度响应数据生成结构速度向量，根据若干所述自由度分别对应的所述加速度响应数据生成训练标签数据；
根据所述发送矩阵、所述接收矩阵、所述结构位移向量以及所述结构速度向量和所述训练标签数据对所述人工神经网络进行训练。
4.根据权利要求3所述的基于图神经网络的结构动力学模型生成方法，其特征在于，所述人工神经网络包括第一多层感知机、第二多层感知机、第三多层感知机、第四多层感知机，所述根据所述发送矩阵、所述接收矩阵、所述结构位移向量以及所述结构速度向量和所述训练标签数据对所述人工神经网络进行训练，包括：
根据所述发送矩阵、所述接收矩阵以及所述结构位移向量生成恢复力计算矩阵，将所述恢复力计算矩阵输入所述第一多层感知机，得到恢复力预测矩阵；
根据所述发送矩阵、所述接收矩阵以及所述结构速度向量生成阻尼力计算矩阵，将所述阻尼力计算矩阵输入所述第二多层感知机，得到阻尼力预测矩阵；
根据所述恢复力预测矩阵和所述阻尼力预测矩阵得到合力的计算矩阵，将所述合力的计算矩阵输入所述第三多层感知机，得到合力预测矩阵，其中，所述合力预测矩阵用于反映恢复力与阻尼力相加的矢量和；
获取外部激励荷载数据，根据所述接收矩阵、所述合力预测矩阵以及所述外部激励荷载数据生成加速度响应计算矩阵，将所述加速度响应计算矩阵输入所述第四多层感知机，得到加速度响应预测数据；
根据所述加速度响应预测数据和所述训练标签数据对所述人工神经网络进行训练。
5.根据权利要求4所述的基于图神经网络的结构动力学模型生成方法，其特征在于，所述根据所述发送矩阵、所述接收矩阵以及所述结构位移向量生成恢复力计算矩阵，包括：
将所述发送矩阵与所述结构位移向量相乘，得到第一矩阵；
将所述接收矩阵与所述结构位移向量相乘，得到第二矩阵；
将所述第一矩阵与所述第二矩阵按列合并，得到所述恢复力计算矩阵。
6.根据权利要求4所述的基于图神经网络的结构动力学模型生成方法，其特征在于，所述根据所述发送矩阵、所述接收矩阵以及所述结构速度向量生成阻尼力计算矩阵，包括：
将所述发送矩阵与所述结构速度向量相乘，得到第三矩阵；
将所述接收矩阵与所述结构速度向量相乘，得到第四矩阵；
将所述第三矩阵与所述第四矩阵按列合并，得到所述阻尼力计算矩阵。
7.根据权利要求4所述的基于图神经网络的结构动力学模型生成方法，其特征在于，所述根据所述恢复力预测矩阵和所述阻尼力预测矩阵得到合力的计算矩阵，包括：
将所述恢复力预测矩阵和所述阻尼力预测矩阵按列合并，得到所述合力的计算矩阵。
8.根据权利要求4所述的基于图神经网络的结构动力学模型生成方法，其特征在于，所述根据所述接收矩阵、所述合力预测矩阵以及所述外部激励荷载数据生成加速度响应计算矩阵，包括：
将所述合力预测矩阵左乘以所述接收矩阵的转置，得到第五矩阵；
将所述第五矩阵与所述外部激励荷载数据按列合并，得到所述加速度响应计算矩阵。
9.一种基于图神经网络的结构动力学模型生成系统，其特征在于，所述系统包括：
数据获取模块，用于确定目标结构，获取所述目标结构的结构响应数据和拓扑信息；
网络训练模块，用于根据所述结构响应数据和所述拓扑信息对预设的人工神经网络进行训练，将训练后的所述人工神经网络作为所述目标结构对应的结构动力学模型。
10.一种基于图神经网络的结构动力学模型，其特征在于，所述结构动力学模型采用上述权利要求1-8任意所述的基于图神经网络的结构动力学模型生成方法生成。
CN202111674213.9A 2021-12-31 2021-12-31 一种基于图神经网络的结构动力学模型生成方法 Active CN114564878B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111674213.9A CN114564878B (zh) 2021-12-31 2021-12-31 一种基于图神经网络的结构动力学模型生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111674213.9A CN114564878B (zh) 2021-12-31 2021-12-31 一种基于图神经网络的结构动力学模型生成方法

Publications (2)

Publication Number Publication Date
CN114564878A true CN114564878A (zh) 2022-05-31
CN114564878B CN114564878B (zh) 2022-10-28

Family

ID=81711195

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111674213.9A Active CN114564878B (zh) 2021-12-31 2021-12-31 一种基于图神经网络的结构动力学模型生成方法

Country Status (1)

Country Link
CN (1) CN114564878B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116305461A (zh) * 2023-03-13 2023-06-23 清华大学 结构响应计算方法、装置、电子设备及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150220311A1 (en) * 2014-02-03 2015-08-06 Richard Salter Computer implemented modeling system and method
CN104898414A (zh) * 2015-03-31 2015-09-09 西南交通大学 一种高速列车动力学性能设计关键参数识别方法
US20210049314A1 (en) * 2019-08-15 2021-02-18 Global Energy Interconnection Research Institute Co. Ltd Systems and Methods of Composite Load Modeling for Electric Power Systems
CN113158543A (zh) * 2021-02-02 2021-07-23 浙江工商大学 一种软件定义网络性能智能预测方法
CN113177626A (zh) * 2021-06-29 2021-07-27 中国科学院自动化研究所 数据驱动的复杂系统机理自动学习方法、系统及设备

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150220311A1 (en) * 2014-02-03 2015-08-06 Richard Salter Computer implemented modeling system and method
CN104898414A (zh) * 2015-03-31 2015-09-09 西南交通大学 一种高速列车动力学性能设计关键参数识别方法
US20210049314A1 (en) * 2019-08-15 2021-02-18 Global Energy Interconnection Research Institute Co. Ltd Systems and Methods of Composite Load Modeling for Electric Power Systems
CN113158543A (zh) * 2021-02-02 2021-07-23 浙江工商大学 一种软件定义网络性能智能预测方法
CN113177626A (zh) * 2021-06-29 2021-07-27 中国科学院自动化研究所 数据驱动的复杂系统机理自动学习方法、系统及设备

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
TOLGA BAKLACIOGLU等: ""Dynamic modeling of exergy efficiency of turboprop engine components using hybrid genetic algorithm-artificial neural networks"", 《ENERGY》 *
俞阿龙: "机器人系统中操作环境动力学模型研究", 《计量学报》 *
张介明: ""利用神经网络建立轮胎动力学模型"", 《上海汽车》 *
袁军社等: "基于神经网络的火箭发动机结构动力学优化", 《火箭推进》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116305461A (zh) * 2023-03-13 2023-06-23 清华大学 结构响应计算方法、装置、电子设备及存储介质
CN116305461B (zh) * 2023-03-13 2023-10-13 清华大学 结构响应计算方法、装置、电子设备及存储介质

Also Published As

Publication number Publication date
CN114564878B (zh) 2022-10-28

Similar Documents

Publication Publication Date Title
Zhang et al. Physics-informed multi-LSTM networks for metamodeling of nonlinear structures
Xu et al. Learning viscoelasticity models from indirect data using deep neural networks
Naveros et al. A spiking neural simulator integrating event-driven and time-driven computation schemes using parallel CPU-GPU co-processing: a case study
Nouri et al. A digital neuromorphic realization of the 2-D Wilson neuron model
Li et al. Structural reliability calculation method based on the dual neural network and direct integration method
Taghizadeh et al. Turbulence closure modeling with data-driven techniques: Investigation of generalizable deep neural networks
CN106650918A (zh) 构建系统模型的方法和装置
Lu et al. Bifidelity data-assisted neural networks in nonintrusive reduced-order modeling
CN114564878B (zh) 一种基于图神经网络的结构动力学模型生成方法
Liu et al. Total solution for structural mechanics problems
Dey et al. Mapping and validating a point neuron model on Intel's neuromorphic hardware Loihi
CN114997036B (zh) 基于深度学习的网络拓扑重构方法、装置和设备
CN114386367B (zh) 一种非线性延迟电路系统的仿真方法、系统及介质
CN114048544B (zh) 飞行器结构响应特性预示方法、装置、设备和介质
Danciu A CNN-based approach for a class of non-standard hyperbolic partial differential equations modeling distributed parameters (nonlinear) control systems
Zhu et al. A bootstrap based virtual sample generation method for improving the accuracy of modeling complex chemical processes using small datasets
CN106021880A (zh) 基于bp神经网络的导管架平台结构响应计算方法
Simpson et al. On the use of nonlinear normal modes for nonlinear reduced order modelling
CN115392594B (zh) 一种基于神经网络和特征筛选的用电负荷模型训练方法
CN113779498B (zh) 离散傅里叶矩阵重构方法、装置、设备和存储介质
EP4109351A1 (en) Inference program, inference method, and information processing apparatus
KR102624710B1 (ko) Gru 기반 구조물 시계열 응답 예측 방법
Taghizadeh Toward turbulence closure modeling with data–driven techniques
Edem et al. One-point quadrature ANS solid-shell element based on a displacement variational formulation Part I–Geometrically linear assessment
CN106803233A (zh) 图像透视变换的优化方法

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