CN117763926B - 一种数字孪生信息驱动的高可信结构变形监测方法 - Google Patents
一种数字孪生信息驱动的高可信结构变形监测方法 Download PDFInfo
- Publication number
- CN117763926B CN117763926B CN202410195544.1A CN202410195544A CN117763926B CN 117763926 B CN117763926 B CN 117763926B CN 202410195544 A CN202410195544 A CN 202410195544A CN 117763926 B CN117763926 B CN 117763926B
- Authority
- CN
- China
- Prior art keywords
- deformation
- strain
- field
- training
- model
- 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 77
- 238000012544 monitoring process Methods 0.000 title claims abstract description 37
- 238000012549 training Methods 0.000 claims abstract description 56
- 238000013528 artificial neural network Methods 0.000 claims abstract description 42
- 238000013507 mapping Methods 0.000 claims abstract description 42
- 238000006073 displacement reaction Methods 0.000 claims abstract description 11
- 238000005070 sampling Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims abstract 2
- 238000004458 analytical method Methods 0.000 claims description 36
- 230000006870 function Effects 0.000 claims description 21
- 238000004088 simulation Methods 0.000 claims description 19
- 230000000694 effects Effects 0.000 claims description 13
- 238000010200 validation analysis Methods 0.000 claims description 9
- 239000013598 vector Substances 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 7
- 238000011156 evaluation Methods 0.000 claims description 6
- 230000004913 activation Effects 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 230000003068 static effect Effects 0.000 claims description 5
- 239000000463 material Substances 0.000 claims description 4
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 claims description 2
- 238000010276 construction Methods 0.000 claims description 2
- 230000008878 coupling Effects 0.000 claims description 2
- 238000010168 coupling process Methods 0.000 claims description 2
- 238000005859 coupling reaction Methods 0.000 claims description 2
- 238000002790 cross-validation Methods 0.000 claims description 2
- 238000012795 verification Methods 0.000 claims 3
- 238000013508 migration Methods 0.000 claims 1
- 230000005012 migration Effects 0.000 claims 1
- 230000001105 regulatory effect Effects 0.000 claims 1
- 238000012163 sequencing technique Methods 0.000 claims 1
- 238000005259 measurement Methods 0.000 abstract description 8
- 238000000691 measurement method Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 8
- 238000013527 convolutional neural network Methods 0.000 description 6
- 101100012902 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) FIG2 gene Proteins 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 241000237942 Conidae Species 0.000 description 2
- 101001121408 Homo sapiens L-amino-acid oxidase Proteins 0.000 description 2
- 102100026388 L-amino-acid oxidase Human genes 0.000 description 2
- 238000013526 transfer learning Methods 0.000 description 2
- 101100233916 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) KAR5 gene Proteins 0.000 description 1
- 229910001069 Ti alloy Inorganic materials 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
Classifications
-
- 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
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
一种数字孪生信息驱动的高可信结构变形监测方法,属于工程结构变形监测领域,步骤如下:建立结构的有限元模型进行采样,获得全场应变和变形样本数据,增加样本数据的多样性;以张量或矩阵形式存储应变和变形场样本数据;基于神经网络训练应变场到变形场样本数据的映射关系;实时获取应变传感器数据,并建立应变场的数字孪生模型;将数字孪生模型预测的全场应变数据信息作为神经网络的输入,预测的全场变形数据作为输出,通过预测结果实现全场变形监测。本发明无需进行位移计布置、变形监测可信度高,可解决传统变形测量方法易受结构振动影响,难以针对工程结构进行随动测量以及传统变形监测方法可信度低、稳定性差等难题。
Description
技术领域
本发明属于工程结构变形监测领域,涉及一种数字孪生信息驱动的高可信结构变形监测方法。
背景技术
目前,对工程结构进行实时高可信度结构变形监测,对保障结构安全和功能实现具有重大意义。针对在轨卫星而言,高可信度结构变形监测能够确保卫星的安全运行和实现对地目标的高可信度稳定监测。由于在轨卫星运行环境复杂且承载功能装备数量较多,因此难以利用大量传感器实现对在轨卫星整体结构的变形监测。同时传统直接测量方法,比如摄影测量法,需安装高清相机,易受结构振动影响产生模糊;传统的激光测量法也难以针对结构三维复杂曲面进行随动测量。随着应变传感器的轻型化、小型化、智能一体化发展趋势,应变数据的测量日益便捷。因此通过传感器实测应变数据实现对结构变形的间接测量是一种有效的解决途径。
现有的基于传感器间接测量的变形监测技术包括模态叠加法(Modal parameterestimation of a simulated non-linear structure,Foss, 1999)、逆有限元法(Beamshape Sensing Using Inverse Finite Element Method: Theory and ExperimentalValidation, Gherlone, 2011)、KO理论及曲率变形监测法(Application of FiniteElastic Theory to the Deformation of Rubbery Materials, Blatz, 1962)等,其中模态叠加法受模型误差和模态截断误差的影响较大,逆有限元法难以构造复杂结构单元形函数,KO理论和曲率重构法对复杂结构的变形监测可信度低和稳定性差。由于现阶段,基于应变对结构变形进行监测的技术普遍存在着可信度低和稳定性较差的缺陷,因此需要实现一种更为简单有效、便于实施、并且满足可信度和稳定性需求的变形监测方法。此外,应变传感器数据只有离散的测点,仿真数据难以反应结构变形真实状态,因此利用应变场预测变形场仍存在困难。
而数字孪生技术作为连接物理世界和数字世界的纽带,通过建立物理系统的数字模型、实时监测系统状态并驱动模型动态更新,可实现系统行为更准确的描述和预测。为全场高可信度应变获取及利用全场高可信度应变预测全场变形、实现高可信度结构变形监测提供了有效途径。
发明内容
本发明主要解决传统变形监测方法计算可信度低和稳定性差的问题,提出一种数字孪生信息驱动的高可信结构变形监测方法。
本发明采用的技术方案为:
一种数字孪生信息驱动的高可信结构变形监测方法,所述的高可信结构变形监测方法,首先,利用有限元分析获得具有多样性的应变场和变形场样本数据。其次,通过神经网络训练结构应变场到变形场的映射关系。再次,利用应变传感器实测应变数据和仿真分析结果构建数字孪生模型获得高可信孪生应变场。最后,将孪生应变场带入神经网络训练得到的应变场到变形场映射关系,实现结构变形场的高可信度预测。包括以下步骤:
第一步,根据结构真实状况建立有限元模型,进行仿真分析并采样,获得具有多样性的全场应变和变形样本数据。具体的:
步骤1.1,首先,基于结构真实受载荷状况确定应变传感器布置方案,布置方案包含传感器数目、种类、布置的位置和方向。
步骤1.2,其次,构建结构的有限元模型并进行网格划分,获得有限元网格节点数目为M。其中,有限元模型的材料属性、边界条件根据结构真实情况设置。并根据1.1)中应变传感器布置方案,建立局部坐标系。
步骤1.3,通过仿真分析方法获得有限元模型的应变和变形。其中,仿真分析方法包括静力分析、模态分析、热力耦合分析(为本领域已知内容)。通过设置不同的仿真分析方法及载荷施加方式,获得N组有限元分析结果作为样本。对每一组有限元分析得到的样本,提取有限元网格模型全部节点的应变和变形作为结构应变场数据和变形场数据,最终获得应变场样本数据和变形场样本数据/>。对每一组样本数据,其中/>、、/>、/>均为大小为/>的列向量,/>、/>、/>分别对应应变传感器测量的应变方向,/>对应需要监测的变形方向。
第二步,存储应变场样本数据和变形场样本数据,存储数据的方式包括张量形式或矩阵形式。此处以张量存储进行说明,具体的:
步骤2.1,提取第一步建立的有限元网格模型节点位置坐标,并依次对三个方向的坐标的数值按由大到小的顺序进行排序,分别获得长度分别为m、n、l的向量X、Y、Z。
步骤2.2,对应每一组应变场和变形场样本数据,建立四个大小为m*n*l的张量、/>、/>、/>,分别存储应变三个方向应变场/>及变形场/>。对于每一个张量中的元素有关系式:
其中,Xi、Yj、Zk分别为向量X、Y、Z中的第i、j、k个元素,i、j、k分别为小于m、n、l的正整数,构成张量中元素的坐标。/>为模型网格节点坐标/>构成的集合,为张量中的第i、j、k个元素值,/>为坐标/>位置处的应变值或变形值。
第三步,基于神经网络训练应变场到变形场的映射关系。采用数据驱动方法,利用应变场样本张量和变形场样本张量/>,训练结构应变场到变形场的映射关系。其中,应变场张量/>包含三个方向的应变场/>。在训练映射关系时,所使用的神经网络包括卷积神经网络、自编码器神经网络、对抗生成神经网络等,具体的:
步骤3.1,将第二步获得的样本数据划分为训练集和验证集,划分方法包括留出法、交叉验证法和自举法。其中训练集用于映射关系的训练、验证集用于神经网络训练效果的评估。
步骤3.2,利用训练集构建神经网络,训练应变场到变形场的映射关系,训练过程中,将应变场样本张量/>作为神经网络的输入,变形场样本张量/>作为神经网络的训练目标。
所述的神经网络的构建包含激活函数、损失函数和优化器的设置。其中激活函数包括Sigmoid函数、ReLU函数,损失函数包括相对均方根误差(RRMSE)、均方误差(MSE)、均方根误差(RMSE),优化器包括自适应矩估计(Adam)、随机梯度下降(SGD)算法。
步骤3.3,利用验证集对训练效果进行评估,并根据评估结果调整神经网络的超参数包括学习率、隐层数目等,以获得效果最好的训练结果,最终获得应变场到变形场的最优映射关系。其中神经网络超参数的调整方法包括人工经验、随机、贝叶斯超参数优化,训练效果的评估方法包括但不限于决定系数(R2)和相对均方根误差(RRMSE),其具体公式如下:
其中,为验证集的数据量,/>表示神经网络训练的目标值,/>表示神经网络的预测值,/>表示神经网络训练的目标值的平均值。R2越接近于1表示变形训练效果越好,RRMSE越接近于0表示训练效果越好。
第四步,实时获取结构处于实际加载环境中,预先布置的应变传感器采集的应变数据,构建数字孪生模型,获得结构的高可信孪生应变场。具体的:
步骤4.1,依据结构工作环境中所受载荷状况,进行试验加载,并实时获取应变传感器数据。
步骤4.2,根据结构实际所受加载状况,设置仿真分析方法和载荷,其余条件的设置与第一步保持一致,得到与结构在真实工作环境中所受载荷一致的仿真分析结果。提取有限元网格模型全部节点的应变作为应变场,并利用第二步提取的有限元网格模型节点位置坐标/>,训练有限元网格模型节点位置坐标/>到应变场的映射关系/>,其中映射关系的训练包括径向基函数代理模型、高斯过程回归、深度神经网络。
步骤4.3,确定应变传感器位置在有限元网格模型节点位置坐标中的坐标/>,读取实时获取的应变传感器测量的应变数据/>。将/>和带入构建的有限元网格模型节点位置坐标到应变的映射关系/>,并对映射关系/>进行修调,构建得到数字孪生模型/>。其中对映射关系/>进行修调并构建数字孪生模型的方法包括加法标度函数、协同克里金、迁移学习等。该数字孪生模型的输入为模型节点坐标/>、输出为孪生应变场/>,同样以张量形式存储孪生应变场/>,最终获得结构的高可信孪生应变场张量/>。
第五步,结合数字孪生模型输出的孪生应变场和神经网络训练得到的应变场到变形场映射关系,获得结构变形场的预测结果。具体的:
将第四步获得的孪生应变场张量作为第三步神经网络训练得到的映射关系/>的输入,得到变形场预测结果/>,实现结构变形场监测。
第六步,为验证变形场预测可信度,在试验中选取变形测点确定其坐标并布设位移计,位移计布设方向为需要监测的变形方向,与第一步中提取的应变场方向保持一致。将位移计测量的变形值与第四步输出的变形场中对应坐标/>的变形预测值进行对比,对变形预测可信度进行评估。可信度评估方法包括但不限于决定系数(R2)和相对均方根误差(RRMSE)。
本发明的有益效果为:
(1)本发明针对传统变形反演方法变形监测可信度低、稳定性差的问题,利用仿真分析获得具有多样性的应变场和变形场样本,基于神经网络训练得到应变场到变形场的映射关系。并基于真实加载状态下应变传感器采集的应变和仿真分析结果,构建结构的数字孪生模型并获得高可信孪生应变场。将孪生应变场输入应变场到变形场的映射关系,从而获得高可信的预测变形场。
(2)本发明对结构变形预测可信度高,适用于难以使用位移计等传统变形监测手段进行变形监测的场景。
(3)本发明的创新点为:1)通过张量或矩阵等形式存储节点应变和变形数据能够较好的保留其空间位置信息,便于神经网络模型的训练;2)基于数据驱动思想,利用样本数据并通过神经网络训练获得了结构应变场到变形场映射关系,解决了复杂结构难以解析推导其几何方程的难题;3)通过数字孪生构建了结构的高可信孪生应变场,能够实现结构应变场的实时高可信度输出,并为结构变形场预测提供更精确的应变场信息输入。
附图说明
图1为一种数字孪生信息驱动的高可信结构变形监测方法的实现流程图。
图2为内压舱体的有限元模型示意图;图2(a)为模型整体示意图;图2(b)为模型的正视图;图2(c)为模型的俯视图。
图3为神经网络构建的结构应变场到变形场映射关系示意图。
图4为内压舱体真实变形场和通过不同方法预测的变形场示意图;图4(a)为真实变形场示意图,图4(b)为模态叠加法预测的变形场示意图,图4(c)为逆有限元法预测的变形场示意图,图4(d)为本方法预测的变形场示意图。
图5为本方法与模态叠加法和逆有限元法对变形预测的可信度对比图。
具体实施方式
为使本发明解决的方法问题、采用的方法方案和达到的方法效果更加清楚,下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部内容。
图1为本发明实施例提供的一种数字孪生信息驱动的高可信结构变形监测方法的实现流程图。如图1所示,本发明实施例提供的一种数字孪生信息驱动的高可信结构变形监测方法,具体包括以下步骤:
第一步,建立结构的有限元模型,进行有限元分析并采样,获得具有多样性的全场应变和变形样本数据。具体的:
步骤1.1,根据结构的真实受载状况制定应变传感器布置方案,在结构外表面布置应变传感器,应变传感器数目为50个,类型为三向应变花,其中应变花内敏感栅夹角为,确定每个应变传感器的方向并记录其位置。
步骤1.2,建立内压舱体的有限元模型,如图2(a)、(b)、(c)所示。模型通过壳体建立,其中模型尺寸a1=1900 mm,为模型的整体高度;a2=100 mm、a3=700 mm、a4=450 mm,分别为模型上下柱壳段的高度、模型中段圆柱壳的高度、模型斜锥壳部分的高度;a5=2000 mm,为模型的整体直径长度;a6=500 mm、a7=350 mm,分别为模型斜锥壳和底部平面宽度。设置模型的材料为高温钛合金,弹性模量E=114 GPa,泊松比v=0.31。通过四边形壳体单元对有限元模型进行划分,共计165325个网格节点。设定边界条件为底端固支,即固定所有方向的自由度。根据1.1)传感器布置方向,建立局部柱坐标系,坐标系原点选择在模型的几何中心,柱坐标系轴向为模型的竖直方向、径向为圆柱壳平面内沿半径方向、周向为与径向垂直的方向。
步骤1.3,采用模态分析方法进行仿真分析,并设置频率分析步计算内压舱体模型前20阶模态,得到每一阶模态对应的模型应变和变形结果。提取每一阶模态对应的沿局部坐标系轴向和周向的应变和/>、轴向与周向平面内的切应变/>以及模型整体变形/>。其中,对每一阶模态下的应变/>和变形/>,均为165325*1的列向量。最终获得应变场样本数据/>和变形场样本数据/>。
第二步,采用张量形式存储第一步获得的应变及变形样本数据,具体的:
步骤2.1,从局部坐标系提取模型网格节点位置坐标构建集合/>,并依次对径向坐标x、周向坐标y、轴向坐标z进行由大到小排序,最终得到大小分别为60、60、57的向量X、Y、Z。
步骤2.2,穷举三个向量中元素的组合,依据判别关系式,即判断每一个组合构成的坐标是否属于集合,如果是则将应变场/>和变形场/>在该坐标位置对应的应变和变形分别存储到四个张量对应的元素中,如果不是则存储0元素,最终获得20组每组4个分别存储应变场/>和变形场/>的60*60*57的张量。
第三步,基于第二步获得的20组应变场和变形场张量数据样本建立卷积神经网络训练应变场到变形场的映射关系。如图3所示,在构建神经网络进行映射关系的训练时,对于每一组样本,将应变场张量同时作为输入、变形场张量/>作为训练目标。具体的:
步骤3.1,将每组数据样本按8:2的方式随机划分为训练集和验证集,其中训练集用于映射关系的训练、验证集用于神经网络训练效果的评估。
步骤3.2,通过训练集构建卷积神经网络训练映射关系,卷积神经网络的激活函数选择ReLU函数,优化器选择Adam优化器,并定义损失函数为多输出回归的RMSE,获得初步的训练结果。
步骤3.3,通过验证集和决定系数(R2)对训练的效果进行评估,并依据评估结果采用贝叶斯优化方法对卷积神经网络的学习率和隐层数目等超参数进行优化,最终获得结构三个方向应变场到变形场的最优映射关系。
第四步,基于结构实际所受载荷状况进行试验加载和仿真分析的设置,具体的:
步骤4.1,根据结构在实际工作环境所受载荷状况进行内压静力试验加载设置,并实时获得每一个应变传感器数据。
步骤4.2,修改第一步有限元分析过程,对模型施加0.12 MPa的内压载荷,其余条件与第一步保持一致,设置分析步为静力分析步计算内压舱体模型静力作用下的应变和变形。提取有限元网格模型全部节点的应变作为结构应变场,并利用第二步提取的有限元网格模型节点位置坐标/>,利用径向基函数代理模型构建有限元网格模型节点位置坐标/>到应变场/>的映射关系/>。
步骤4.3,根据应变传感器位置确定其在有限元网格模型节点位置坐标中的坐标/>,并实时读取应变传感器测量的应变数据/>。将/>和结合带入的映射关系/>,利用迁移学习方法构建数字孪生模型获得孪生应变场并进行张量形式存储,得到高可信的孪生应变场张量/>。
第五步,将第四步获得的孪生应变场张量输入到第三步通过卷积神经网络训练完成的应变场到变形场的映射关系/>,得到高可信的结构变形场预测结果,并与真实变形场和使用传统模态叠加法、传统逆有限元法得到的结构变形场预测结果的对比如图4所示,其中,颜色的深度表示变形值的大小,颜色越深表示变形越小,颜色越浅表示变形越大,变形值的单位为毫米。从图4中可以看出,本方法对结构变形的预测结果与真实结果最为接近,变形预测精度显著高于传统模态叠加法和传统逆有限元法。
第六步,在试验中选取变形测点确定其坐标并布设位移计,位移计数量为20个。位移计布设方向为需要监测的变形方向,并与第一步提取的变形方向保持一致。通过变形测点实测变形值和第四步输出的变形场中对应坐标/>的变形预测值进行对比,通过决定系数(R2)和相对均方根误差(RRMSE)对结构变形场预测可信度进行评估,R2越接近于1表示变形预测结果可信度越高,RRMSE越接近于0表示变形预测结果可信度越高。
利用位移计测量20个测点的变形数据,通过R2和RRMSE对变形场预测的可信度进行评估,计算得到R2为0.9658,RRMSE为0.1865。如图5所示,将本方法得到的全场变形预测结果的可信度,与基于传统模态叠加法和传统逆有限元法得到的变形预测结果的可信度进行对比,本方法的预测可信度分别提高18.1%和10.4%。
表1变形预测可信度对比
最后应说明的是:以上各实施例仅用以说明本发明的方法方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通方法人员应当理解:其对前述各实施例所记载的方法方案进行修改,或者对其中部分或者全部方法特征进行等同替换,并不使相应方法方案的本质脱离本发明各实施例方法方案的范围。
Claims (7)
1.一种数字孪生信息驱动的高可信结构变形监测方法,其特征在于,所述的高可信结构变形监测方法,首先,利用有限元分析获得具有多样性的应变场和变形场样本数据;其次,通过神经网络训练结构应变场到变形场的映射关系;再次,利用应变传感器实测应变数据和仿真分析结果构建数字孪生模型获得高可信孪生应变场;最后,将孪生应变场带入神经网络训练得到的应变场到变形场映射关系,实现结构变形场的高可信度预测;具体包括以下步骤:
第一步,根据结构真实状况建立有限元模型,进行仿真分析并采样,获得具有多样性的全场应变和变形样本数据;具体的:
步骤1.1,基于结构真实受载荷状况确定应变传感器布置方案,布置方案包含传感器数目、种类、布置的位置和方向;
步骤1.2,构建结构的有限元模型并进行网格划分,获得有限元网格节点数目为M;其中,有限元模型的材料属性、边界条件根据结构真实情况设置;并根据1.1)中应变传感器布置方案,建立局部坐标系;
步骤1.3,通过仿真分析方法获得有限元模型的应变和变形,通过设置不同的仿真分析方法及载荷施加方式,获得N组有限元分析结果作为样本;对每一组有限元分析得到的样本,提取有限元网格模型全部节点的应变和变形作为结构应变场数据和变形场数据,最终获得应变场样本数据{εx,εy,εxy}N和变形场样本数据{w}N;对每一组样本数据,其中εx、εy、εxy、w均为大小为M*1的列向量,εx,εy,εxy分别对应应变传感器测量的应变方向,w对应需要监测的变形方向;
第二步,存储应变场样本数据和变形场样本数据,存储数据的方式包括张量形式或矩阵形式;此处以张量存储进行说明,具体的:
步骤2.1,提取第一步建立的有限元网格模型节点位置坐标{x,y,z},并依次对三个方向的坐标的数值按由大到小的顺序进行排序,分别获得长度分别为m、n、l的向量X、Y、Z;
步骤2.2,对应每一组应变场和变形场样本数据,建立四个大小为m*n*l的张量Φx、Φy、Φxy、ψ,分别存储应变三个方向应变场εx、εy、εxy及变形场w;对于每一个张量中的元素有关系式:
if(Xi,Yj,Zk)∈Ω
ei,j,k=E(Xi,Yj,Zk)
else
ei,j,k=0
其中,Xi、Yj、Zk分别为向量X、Y、Z中的第i、j、k个元素,i、j、k分别为小于m、n、l的正整数,(Xi,Yj,Zk)构成张量中元素的坐标;Ω为模型网格节点坐标{x,y,z}构成的集合,ei,j,k为张量中的第i、j、k个元素值,E(Xi,Yj,Zk)为坐标(Xi,Yj,Zk)位置处的应变值或变形值;
第三步,基于神经网络训练应变场到变形场的映射关系;采用数据驱动方法,利用应变场样本张量{Φ}N和变形场样本张量{ψ}N,训练结构应变场到变形场的映射关系其中,应变场张量Φ包含三个方向的应变场{Φx,Φv,Φxy};具体的:
步骤3.1,将第二步获得的样本数据划分为训练集和验证集,其中训练集用于映射关系的训练、验证集用于神经网络训练效果的评估;
步骤3.2,利用训练集构建神经网络,训练应变场到变形场的映射关系f,训练过程中,将应变场样本张量{Φ}N作为神经网络的输入,变形场样本张量{ψ}N作为神经网络的训练目标;
步骤3.3,利用验证集对训练效果进行评估,并根据评估结果调整神经网络的超参数,获得训练结果,最终获得应变场到变形场的最优映射关系f;
第四步,实时获取结构处于实际加载环境中,预先布置的应变传感器采集的应变数据,构建数字孪生模型,获得结构的高可信孪生应变场;具体的:
步骤4.1,依据结构工作环境中所受载荷状况,进行试验加载,并实时获取应变传感器数据;
步骤4.2,根据结构实际所受加载状况,设置仿真分析方法和载荷,其余条件的设置与第一步保持一致,得到与结构在真实工作环境中所受载荷一致的仿真分析结果;提取有限元网格模型全部节点的应变作为应变场{ε′x,ε′y,ε′xy},并利用第二步提取的有限元网格模型节点位置坐标{x,y,z},训练有限元网格模型节点位置坐标{x,y,z}到应变场{ε′x,ε′y,ε′xy}的映射关系其中映射关系的训练包括径向基函数代理模型、高斯过程回归、深度神经网络;
步骤4.3,确定应变传感器位置在有限元网格模型节点位置坐标{x,y,z}中的坐标读取实时获取的应变传感器测量的应变数据/>将/>和带入构建的有限元网格模型节点位置坐标到应变的映射关系ξ,并对映射关系ξ进行修调,构建得到数字孪生模型/>该数字孪生模型的输入为模型节点坐标{x,y,z}、输出为孪生应变场/>同样以张量形式存储孪生应变场最终获得结构的高可信孪生应变场张量Φ′={Φ′x,Φ′y,Φ′xy};
第五步,结合数字孪生模型输出的孪生应变场和神经网络训练得到的应变场到变形场映射关系f,获得结构变形场的预测结果;具体的:
将第四步获得的孪生应变场张量Φ′作为第三步神经网络训练得到的映射关系f的输入,得到变形场预测结果ψ′,实现结构变形场监测;
第六步,为验证变形场预测可信度,在试验中选取变形测点确定其坐标并布设位移计,位移计布设方向为需要监测的变形方向,与第一步中提取的应变场方向保持一致;将位移计测量的变形值与第四步输出的变形场中对应坐标/>的变形预测值进行对比,对变形预测可信度进行评估。
2.根据权利要求1所述的一种数字孪生信息驱动的高可信结构变形监测方法,其特征在于,所述的步骤1.3中,仿真分析方法包括静力分析、模态分析、热力耦合分析。
3.根据权利要求1所述的一种数字孪生信息驱动的高可信结构变形监测方法,其特征在于,所述的步骤3.1中,划分样本数据的方法包括留出法、交叉验证法和自举法。
4.根据权利要求1所述的一种数字孪生信息驱动的高可信结构变形监测方法,其特征在于,所述的步骤3.2中,神经网络的构建包含激活函数、损失函数和优化器的设置;其中激活函数包括Sigmoid函数、ReLU函数,损失函数包括相对均方根误差RRMSE、均方误差MSE、均方根误差RMSE,优化器包括自适应矩估计算法Adam、随机梯度下降算法SGD。
5.根据权利要求1所述的一种数字孪生信息驱动的高可信结构变形监测方法,其特征在于,所述的步骤3.3中,神经网络超参数的调整方法包括人工经验、随机、贝叶斯超参数优化,训练效果的评估方法包括决定系数R2和相对均方根误差RRMSE,其具体公式如下:
其中,K为验证集的数据量,mi表示神经网络训练的目标值,表示神经网络的预测值,表示神经网络训练的目标值的平均值R2越接近于1表示变形训练效果越好,RRMSE越接近于0表示训练效果越好。
6.根据权利要求1所述的一种数字孪生信息驱动的高可信结构变形监测方法,其特征在于,所述的步骤4.3中,对映射关系ξ进行修调并构建数字孪生模型的方法包括加法标度函数、协同克里金、迁移学习。
7.根据权利要求1所述的一种数字孪生信息驱动的高可信结构变形监测方法,其特征在于,所述的第六步中,可信度评估方法包括决定系数R2和相对均方根误差RRMSE。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410195544.1A CN117763926B (zh) | 2024-02-22 | 2024-02-22 | 一种数字孪生信息驱动的高可信结构变形监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410195544.1A CN117763926B (zh) | 2024-02-22 | 2024-02-22 | 一种数字孪生信息驱动的高可信结构变形监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117763926A CN117763926A (zh) | 2024-03-26 |
CN117763926B true CN117763926B (zh) | 2024-05-17 |
Family
ID=90314820
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410195544.1A Active CN117763926B (zh) | 2024-02-22 | 2024-02-22 | 一种数字孪生信息驱动的高可信结构变形监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117763926B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117973158B (zh) * | 2024-04-01 | 2024-08-09 | 大连理工大学 | 一种基于数字孪生的大型结构模态全域监测方法 |
CN118052114B (zh) * | 2024-04-16 | 2024-06-28 | 大连理工大学 | 基于“算测融合”的导弹数字孪生建模方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112906136A (zh) * | 2021-03-16 | 2021-06-04 | 上海交通大学 | 船体外板的激光热成形变形预测方法和系统 |
CN112926152A (zh) * | 2021-02-15 | 2021-06-08 | 西北工业大学 | 一种数字孪生驱动的薄壁件装夹力精准控制与优化方法 |
CN113932761A (zh) * | 2021-10-13 | 2022-01-14 | 山东大学 | 基于数字孪生的高速列车承载结构变形监测方法及系统 |
CN114782658A (zh) * | 2022-04-19 | 2022-07-22 | 哈尔滨工业大学 | 基于数字孪生与增强现实的结构试验监测预警系统及方法 |
CN115392077A (zh) * | 2022-08-16 | 2022-11-25 | 中国人民解放军军事科学院国防科技创新研究院 | 基于深度学习的卫星物理场级数字孪生模型构建方法 |
CN116579209A (zh) * | 2023-05-11 | 2023-08-11 | 南京航空航天大学 | 基于数字孪生和神经网络的机翼结构变形实时预测方法 |
CN116911115A (zh) * | 2023-07-06 | 2023-10-20 | 华南理工大学 | 一种数字孪生隧道与智能建造方法及系统 |
CN117034430A (zh) * | 2023-09-11 | 2023-11-10 | 湘潭大学 | 基于深度学习及数字孪生的闸门健康监测方法、系统及计算机可读存储介质 |
CN117387559A (zh) * | 2023-12-12 | 2024-01-12 | 华东交通大学 | 基于数字孪生的混凝土桥监测系统及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10636142B2 (en) * | 2018-04-20 | 2020-04-28 | Siemens Healthcare Gmbh | Real-time and accurate soft tissue deformation prediction |
-
2024
- 2024-02-22 CN CN202410195544.1A patent/CN117763926B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112926152A (zh) * | 2021-02-15 | 2021-06-08 | 西北工业大学 | 一种数字孪生驱动的薄壁件装夹力精准控制与优化方法 |
CN112906136A (zh) * | 2021-03-16 | 2021-06-04 | 上海交通大学 | 船体外板的激光热成形变形预测方法和系统 |
CN113932761A (zh) * | 2021-10-13 | 2022-01-14 | 山东大学 | 基于数字孪生的高速列车承载结构变形监测方法及系统 |
CN114782658A (zh) * | 2022-04-19 | 2022-07-22 | 哈尔滨工业大学 | 基于数字孪生与增强现实的结构试验监测预警系统及方法 |
CN115392077A (zh) * | 2022-08-16 | 2022-11-25 | 中国人民解放军军事科学院国防科技创新研究院 | 基于深度学习的卫星物理场级数字孪生模型构建方法 |
CN116579209A (zh) * | 2023-05-11 | 2023-08-11 | 南京航空航天大学 | 基于数字孪生和神经网络的机翼结构变形实时预测方法 |
CN116911115A (zh) * | 2023-07-06 | 2023-10-20 | 华南理工大学 | 一种数字孪生隧道与智能建造方法及系统 |
CN117034430A (zh) * | 2023-09-11 | 2023-11-10 | 湘潭大学 | 基于深度学习及数字孪生的闸门健康监测方法、系统及计算机可读存储介质 |
CN117387559A (zh) * | 2023-12-12 | 2024-01-12 | 华东交通大学 | 基于数字孪生的混凝土桥监测系统及方法 |
Non-Patent Citations (1)
Title |
---|
《数据驱动的曲面构件形状-拓扑协同优化方法》;高天贺等;《航空学报》;20240125;第45卷(第2期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN117763926A (zh) | 2024-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN117763926B (zh) | 一种数字孪生信息驱动的高可信结构变形监测方法 | |
CN111209708B (zh) | 基于机器学习的桩-土相互作用预测分析方法 | |
Ghaboussi et al. | New nested adaptive neural networks (NANN) for constitutive modeling | |
CN111914481B (zh) | 基于改进的pso-rbf算法的边坡稳定性预测方法 | |
CN114117617B (zh) | 一种三跨内廊式rc框架建筑地震响应的快速预测方法 | |
Johari et al. | Modelling the mechanical behaviour of unsaturated soils using a genetic algorithm-based neural network | |
CN109885916B (zh) | 一种基于lssvm的混合试验在线模型更新方法 | |
CN115017833B (zh) | 基于深度神经网络的高地应力软岩体地应力计算方法 | |
CN109614580A (zh) | 基于在线Xgboost算法的抗震混合试验模型更新方法 | |
CN116306178B (zh) | 基于自适应形函数和等效中性层的结构应变反演方法 | |
CN116522774A (zh) | 基于迁移学习算法的土质边坡稳定性预测方法及预测平台 | |
Fu et al. | Integration of laboratory testing and constitutive modeling of soils | |
CN114117840A (zh) | 一种基于仿真和试验数据混合驱动的结构性能预测方法 | |
CN115455793A (zh) | 基于多尺度模型修正的高层结构复杂构件受力分析方法 | |
CN114896672A (zh) | 一种基于CSO-BiLSTM网络的桥梁动态位移重构方法 | |
CN110631792A (zh) | 基于卷积神经网络的抗震混合试验模型更新方法 | |
CN114996995A (zh) | 一种超材料隔振单元性能预报方法和系统 | |
CN116956440A (zh) | 一种复杂环境下混凝土性态监测点优化布置方法 | |
Ding et al. | Jaya-based long short-term memory neural network for structural damage identification with consideration of measurement uncertainties | |
CN114861519B (zh) | 复杂地质条件下初始地应力场加速优化反演方法 | |
CN110362902B (zh) | 一种基于区间逐维分析的单源动载荷识别方法 | |
CN111177966B (zh) | 一种基于贝叶斯理论的导弹结构不确定载荷区间重构方法 | |
CN118866211B (zh) | 基于数据驱动的颗粒增强金属基复合材料的本构建模方法及应力预测方法 | |
CN114154725A (zh) | 一种基于机器学习算法的地基承载力特征值预测方法 | |
CN117787067B (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 |