CN113722961A - 一种基于广义协同高斯过程模型的结构不确定性量化解析方法 - Google Patents

一种基于广义协同高斯过程模型的结构不确定性量化解析方法 Download PDF

Info

Publication number
CN113722961A
CN113722961A CN202111024605.0A CN202111024605A CN113722961A CN 113722961 A CN113722961 A CN 113722961A CN 202111024605 A CN202111024605 A CN 202111024605A CN 113722961 A CN113722961 A CN 113722961A
Authority
CN
China
Prior art keywords
precision
low
gaussian process
process model
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.)
Granted
Application number
CN202111024605.0A
Other languages
English (en)
Other versions
CN113722961B (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202111024605.0A priority Critical patent/CN113722961B/zh
Publication of CN113722961A publication Critical patent/CN113722961A/zh
Application granted granted Critical
Publication of CN113722961B publication Critical patent/CN113722961B/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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/08—Probabilistic or stochastic CAD

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开了一种基于广义协同高斯过程模型的结构不确定性量化解析方法，包括以下步骤：A.获得高、低精度的训练样本集；B.利用所述的低精度训练样本集建立低精度高斯过程模型；C.利用所述的高精度训练样本集建立差值高斯过程模型；D.利用所述低精度高斯过程模型和差值高斯过程模型共同构成广义协同高斯过程模型；E.在所述广义协同高斯过程模型框架里解析地计算结构响应的均值和方差。本发明整合了高、低精度的样本数据，用数学模型代替结构物理模型，能够以较低成本达到较高的计算精度；同时本发明所提解析方法能够快速直接地得到均值和方差的解析结果，准确量化结构的不确定性大小。

Description

一种基于广义协同高斯过程模型的结构不确定性量化解析 方法
技术领域
本发明涉及工程结构不确定性量化领域，具体为基于广义协同高斯过程模型的结构不确定性量化解析方法。
背景技术
土木工程的结构参数不可避免存在不确定性，可大致归纳为以下几个方面：结构材料物理特性的不确定性(如材料的弹性模量、质量密度的不确定性)；几何尺寸的不确定性(如构件的长度和截面高、宽度的不确定性)；结构边界条件的不确定性(如构件之间连接的不确定性)。引起结构参数不确定性的因素很多，比如加工容差、装配磨损、环境侵蚀、参数自身固有的随机性。传统的结构设计分析方法(如安全系数法)，一般将相关的不确定性参数当作确定量来考虑，再根据经验或以往资料确定安全系数的值，最终得到确定的计算结果。但参数不确定性必然导致结构响应具有不确定性，若简单地忽略这些不确定性，会造成对结构的认识不足，导致结构设计无法准确地反映结构的实际情况，可能会使结构处于一种不安全的状态。因此，为了提高工程结构的可靠度，为设计人员进行工程结构设计与决策提供更准确的信息，需要量化参数不确定性传递给结构的不确定性大小。
蒙特卡洛法(MCS)是常用的不确定性量化方法，需要对不确定性参数进行大量采样，然后进行相应的有限元模型计算，再对模型计算结果进行统计分析。MCS方法具有适用范围广、稳定性好、易实现的优点，但大量次数的有限元模型分析会导致计算成本很高，难以应用于大型复杂结构的不确定性量化。为克服MCS法计算效率低的不足，代理模型方法采用数学模型近似代替结构物理模型，后续不确定性量化无需原始物理模型，大大降低了计算成本。其中，一种基于高斯过程模型(GPM)的结构不确定性量化方法，将均值和方差的复杂高维积分问题转化为简单一维积分问题，得到了统计矩的解析结果。
GPM的建立涉及到训练样本，训练样本精度越高，建立的模型越准确。高精度样本数据的获取需要建立复杂的有限元模型，这在一定程度上也会增加计算成本，导致采用高精度样本建模效率较低。
综上所述，要通过较低的计算成本达到满意的建模精度，需要一种采用较少高精度样本数据就能获得良好预测精度的模型，以较低成本得到较高精度的统计矩计算结果，准确地量化结构不确定性大小。
发明内容
本发明针对代理模型用来不确定性量化，存在高精度样本计算成本高和低精度样本建模精度低的问题，提出一种基于广义协同高斯过程模型(GC-GPM)的结构不确定性量化解析方法。
本发明的目的是通过以下技术方案来实现的：一种基于广义协同高斯过程模型的结构不确定性量化解析方法，包括以下步骤：
A.获得高、低精度的训练样本集
A1.采样概率分布服从正态或均匀分布的结构不确定性参数，对于结构不确定性参数的概率分布不服从正态或均匀分布时，根据概率相等的原则将其转化为服从正态或均匀分布的参数ui后再采样，转化表达式如下：
Figure BDA0003241956340000021
其中，
Figure BDA0003241956340000022
是参数xi的概率分布函数，G表示正态分布或均匀分布函数；然后通过拉丁超立方抽样获得结构不确定性参数的两组样本点，数量分别为n1和n2，将n1个样本点作为低精度输入值X1，n2个样本点作为高精度输入值X2，其中，n1＞n2；
A2.在ANSYS中建立低精度和高精度的结构有限元模型，分别利用所述的低精度输入值X1和高精度输入值X2进行对应的有限元模型的循环计算，得到低精度观测值Y1和高精度观测值Y2，分别组成低精度训练样本集
Figure BDA0003241956340000031
和高精度训练样本集
Figure BDA0003241956340000032
B.利用所述的低精度训练样本集建立低精度高斯过程模型；具体为：
将低精度训练样本集
Figure BDA0003241956340000033
代入高斯过程模型进行建模，假定低精度高斯过程模型超参数为
Figure BDA0003241956340000034
Figure BDA0003241956340000035
表示低精度高斯过程模型的过程方差；l1，k表示当低精度输入值从x1，k到x′1，k时低精度输出值y1的变化速度，l1，k越大，速度越慢，k表示向量的第k个维度，共有d维。根据最大化边缘似然函数MLE得到超参数Θ1，获得训练完成后的低精度高斯过程模型y1：
Figure BDA0003241956340000036
式中，
Figure BDA0003241956340000037
表示高斯过程，
Figure BDA0003241956340000038
和
Figure BDA0003241956340000039
分别表示低精度高斯过程模型的预测值和预测方差，可写为：
Figure BDA00032419563400000310
Figure BDA00032419563400000311
式中，
Figure BDA00032419563400000312
C1＝C(X1，X1)是样本点集X1的协方差函数矩阵，表示低精度样本点之间的相关性；
Figure BDA00032419563400000314
表示待预测点x*与低精度输入值X1之间的相关性，x1，i(i＝1，2，...，n1)表示组成低精度输入值X1的第i个样本点；
Figure BDA00032419563400000313
e1表示长度为n1的单位向量。
C.利用所述的高精度训练样本集建立差值高斯过程模型
C1.假设最终预测值可由所述低精度高斯过程模型的预测值和一个高斯误差线性表示：
Figure BDA0003241956340000041
其中，y1(x)是低精度模型的真实值；y2(x)是高精度模型的真实值；δ2(x)表示高斯过程误差；⊥表示相互独立；ρ1是一个比例系数，通过最大化边缘似然函数(MLE)获得；
C2.考虑步骤B中所述低精度高斯过程模型的预测值误差，预测值与真实值之间的关系可表示为：
Figure BDA0003241956340000042
其中，ε1(x)表示低精度高斯过程模型的预测值误差，对于两个不同的输入值，ε1(x)是相互独立的，即对任何x≠x′，有Cov(ε(x)，ε(x′))＝0；y1(x)服从正态分布：
Figure BDA0003241956340000043
C3.联立步骤C1和步骤C2中公式，可将y2(x)与ρ1y1(x)的差值δ2(x)表示为：
Figure BDA0003241956340000044
其中，δ2(x)服从正态分布：
Figure BDA0003241956340000045
C4.将步骤C3中所述的
Figure BDA0003241956340000046
记为δ′2，将高精度训练样本集
Figure BDA0003241956340000047
代入高斯过程模型的建立过程，假定差值高斯过程模型超参数为
Figure BDA0003241956340000048
Figure BDA0003241956340000049
表示差值高斯过程模型的过程方差；l2，k表示当高精度输入值从x2，k到x′2，k时差值δ′2的变化速度，l2，k越大，速度越慢，k表示向量的第k个维度，共有d维。根据MLE得到超参数Θ2，获得训练完成后的差值高斯过程模型：
Figure BDA00032419563400000410
式中，
Figure BDA00032419563400000411
表示高斯过程，
Figure BDA00032419563400000412
和
Figure BDA00032419563400000413
分别表示差值高斯过程模型的预测值和预测方差，可写为：
Figure BDA0003241956340000051
Figure BDA0003241956340000052
式中，
Figure BDA0003241956340000053
C2＝C(X2，X2)是X2的协方差函数矩阵；
Figure BDA0003241956340000054
表示低精度高斯过程模型在X2处的预测方差；
Figure BDA0003241956340000055
表示待预测点x*与X2之间的相关性，x2，i(i＝1，2，...，n2)表示组成高精度输入值X2的第i个样本点；
Figure BDA0003241956340000056
e2表示长度为n2的单位向量。
D.利用低精度高斯过程模型和差值高斯过程模型共同构成广义协同高斯过程模型；具体为：
将步骤B和步骤C所述的低精度高斯过程模型和差值高斯过程模型线性组合，共同构成广义协同高斯过程模型，得到最终的高精度预测值
Figure BDA0003241956340000057
和预测方差
Figure BDA0003241956340000058
Figure BDA0003241956340000059
Figure BDA00032419563400000510
E.通过步骤D构建的广义协同高斯过程模型解析地计算结构响应的均值和方差；
E1.根据统计原理，均值
Figure BDA00032419563400000511
和方差
Figure BDA00032419563400000512
的表达式为：
Figure BDA00032419563400000513
Figure BDA00032419563400000514
式中，p(·)表示概率密度函数；x表示输入值；
Figure BDA0003241956340000061
表示预测值，vy表示预测方差。
通过将协方差函数分离为正态分布密度函数的形式，可将步骤D中所述
Figure BDA0003241956340000062
和
Figure BDA0003241956340000063
的不可积分形式转化为可积分形式：
Figure BDA0003241956340000064
Figure BDA0003241956340000065
式中，
Figure BDA0003241956340000066
α1(2)，i(j)表示α1(2)的第i(j)个元素；Aji表示
Figure BDA0003241956340000067
第j行第i列的元素；Bji表示
Figure BDA0003241956340000068
第j行第i列的元素；
Figure BDA0003241956340000069
Figure BDA00032419563400000610
p1(2)，i表示p1(2)的第i个元素；q1＝C1e1，q2＝Λ2e2，q1(2)，j表示q1(2)的第j个元素；d表示向量的维度；
Figure BDA00032419563400000611
xk表示x*的第k个元素；
Figure BDA00032419563400000612
表示样本点集X1(2)第i(j)行、第k列的元素。
E2.将步骤E1中所述的
Figure BDA00032419563400000613
和
Figure BDA00032419563400000614
代入
Figure BDA00032419563400000615
和
Figure BDA00032419563400000616
中，基于广义协同高斯过程模型，将均值和方差的高维积分转化为一维积分：
Figure BDA0003241956340000071
Figure BDA0003241956340000072
E3.E2中所述的一维积分
Figure BDA0003241956340000073
和
Figure BDA0003241956340000074
可统一地表达为：
Figure BDA0003241956340000075
式中，
Figure BDA0003241956340000076
表示均值为ζ，方差为
Figure BDA0003241956340000077
的正态概率密度函数；p(xk)表示参数xk的概率密度函数。
当参数为正态分布或均匀分布时，通过解析地计算一维积分，最终可以直接得到均值和方差的解析结果，量化结构不确定性大小，解析结果为：
Figure BDA0003241956340000081
式中，
Figure BDA0003241956340000082
表示参数xk服从均值为ξ，方差为θ2的正态分布；
Figure BDA0003241956340000083
表示参数xk服从上下限分别为
Figure BDA0003241956340000084
和x的均匀分布；Φ(·)表示标准正态概率分布函数。
进一步地，在步骤C中，由于考虑了所述低精度高斯过程模型在所述高精度输入值X2处的预测误差ε1(x)，因此不需要满足嵌套样本数据
Figure BDA0003241956340000085
的限制条件。
进一步地，步骤A2中，根据结构的复杂程度，在ANSYS能够计算的条件下，将结构几何模型划分较少的单元作为低精度有限元模型，将划分单元数为低精度有限元模型的4-5倍作为高精度有限元模型；或者将复杂原始模型作为高精度有限元模型，通过结构受力方式及自重相等进行等效简化的简单模型作为低精度有限元模型。通过低精度模型计算得到的低精度观测值与真实值偏差较大，因此称为低精度模型。
进一步地，步骤B和步骤C4中，高斯过程模型的建立过程包括以下步骤：
高斯过程模型表达式如下：
Figure BDA0003241956340000086
式中，
Figure BDA0003241956340000087
表示高斯过程；m(x)为均值函数，采用常数形式μ；C(x，x′)为平方指数协方差函数，表示为：
Figure BDA0003241956340000088
式中，x和x′是输入值；xk与x′k分别是x和x′的第k个元素，d表示输入值的维度；高斯过程模型的超参数定义为Θ＝{l1，l2，...，lk，...，ld，η2}，η2表示过程方差；lk表示输出值在输入值从xk到x′k时的变化速度，lk越大，速度越慢；
假设有n个观测值的训练样本集
Figure BDA00032419563400000911
其中，
Figure BDA0003241956340000091
Figure BDA00032419563400000912
根据模型输出Y服从高斯分布的先验假设，则有：
Figure BDA0003241956340000092
同样的，待预测点x*处的预测值y*与Y也服从高斯分布：
Figure BDA0003241956340000093
根据贝叶斯原理，预测值y*的后验分布为：
Figure BDA0003241956340000094
将所述P(Y)和P(Y，y*)代入后验分布的计算公式，可得：
Figure BDA0003241956340000095
其中，
Figure BDA0003241956340000096
Figure BDA0003241956340000097
式中，
Figure BDA00032419563400000913
C＝C(X，X)；α＝C-1(Y-eμ*)；
Figure BDA00032419563400000914
e表示长度为n的单位向量。
其中，高斯过程模型的超参数Θ＝{l1，l2，...，lk，...，ld，η2}可通过最大化边缘似然函数(MLE)求得，即最小化负对数边缘似然函数
Figure BDA0003241956340000098
Figure BDA0003241956340000099
Figure BDA00032419563400000910
本发明的优势是：
1.相较于传统的蒙特卡洛方法，本发明提出的基于广义协同高斯过程模型(GC-GPM)的结构不确定性量化解析方法采用数学模型代替结构物理模型，省去了大量原始模型计算，大大降低了计算成本，同时保证了计算结果有良好的精度；
2.本发明采用整合了高、低精度样本数据的GC-GPM，低精度样本用来拟合原始模型的整体趋势，较少的高精度样本用来对该趋势进行修正，以较低成本达到满意的建模精度，与一般的高斯过程模型(GPM)方法相比具有高效率的优势；
3.本发明基于的GC-GPM框架考虑了所述低精度高斯过程模型的预测值误差，相较于一般的多精度高斯过程模型，本发明适用范围更广，不仅适用于嵌套样本数据
Figure BDA0003241956340000101
还适用于非嵌套样本数据
Figure BDA0003241956340000102
4.本发明是基于GC-GPM推导的结构不确定性量化解析方法，相较于使用一般的高斯过程模型和其他多精度高斯过程模型直接预测再统计预测结果，本发明所提解析方法能够快速直接地得到均值和方差的解析结果，准确量化结构的不确定性大小。
附图说明
图1：基于GC-GPM的结构不确定性量化解析方法流程图；
图2：GC-GPM的建模步骤流程图；
图3：单层球面壳的高低精度有限元模型。
具体实施方式
下面结合一单层球面壳的固有频率统计矩计算进一步说明本发明。假定钢管直径、钢材密度和弹性模量为不确定性参数，分别服从均匀、正态和对数正态分布，计算不确定参数下网壳结构固有频率的均值和方差。利用本发明所提解析方法进行计算，具体包括以下步骤：
A.获得高、低精度的训练样本集。
A1.结构不确定性参数的概率分布包括但不限于正态分布、均匀分布、对数正态分布、威布尔分布、贝塔分布、伽马分布，将不服从正态或均匀分布的单层球面壳的结构不确定性参数(钢管直径、钢材密度和弹性模量)根据概率相等的原则转化为正态或均匀分布，再通过拉丁超立方抽样获得所述结构不确定性参数的两组样本点，数量分别为n1、n2，将n1个样本点作为低精度输入值X1，n2个样本点作为高精度输入值X2，其中，n1＞n2。本算例中，n1＝30，n2＝12；
A2.利用所述的低精度输入值X1和高精度输入值X2进行对应的低精度有限元模型的循环计算，得到单层网壳固有频率的低精度数据Y1和高精度数据Y2，分别组成低精度训练样本集
Figure BDA0003241956340000111
和高精度训练样本集
Figure BDA0003241956340000112
本算例中，采用ANSYS软件建立该网壳的有限元模型，所有杆件均采用BEAM188梁单元模拟，根据结构的复杂程度，在ANSYS能够计算的条件下，将结构几何模型划分较少的单元作为低精度有限元模型，将划分单元数为低精度有限元模型的4-5倍作为高精度有限元模型；或者将复杂原始模型作为高精度有限元模型，通过结构受力方式及自重相等进行等效简化的简单模型作为低精度有限元模型。通过低精度模型计算得到的低精度观测值与真实值偏差较大，因此称为低精度模型。高精度有限元模型将每根梁单元划分5个单元，共375个单元；低精度有限元模型将每个梁单元划分1个单元，共75个单元。高、低精度有限元模型如图3所示。
B.利用所述的低精度训练样本集
Figure BDA0003241956340000113
建立低精度高斯过程模型；具体为：
高斯过程模型的建立过程包括以下步骤：
高斯过程模型表达式如下：
Figure BDA0003241956340000121
式中，
Figure BDA0003241956340000122
表示高斯过程；m(x)为均值函数，采用常数形式μ；C(x，x′)为平方指数协方差函数，表示为：
Figure BDA0003241956340000123
式中，x和x′是输入值；xk与x′k分别是x和x′的第k个元素，d表示输入值的维度；高斯过程模型的超参数定义为Θ＝{l1，l2，...，lk，...，ld，η2}，η2表示过程方差；lk表示输出值在输入值从xk到x′k时的变化速度，lk越大，速度越慢；
假设有n个观测值的训练样本集
Figure BDA0003241956340000124
其中，
Figure BDA0003241956340000125
Figure BDA00032419563400001212
根据模型输出Y服从高斯分布的先验假设，则有：
Figure BDA0003241956340000126
同样的，待预测点x*处的预测值y*与Y也服从高斯分布：
Figure BDA0003241956340000127
根据贝叶斯原理，预测值y*的后验分布为：
Figure BDA0003241956340000128
将所述P(Y)和P(Y，y*)代入后验分布的计算公式，可得：
Figure BDA0003241956340000129
其中，
Figure BDA00032419563400001210
Figure BDA00032419563400001211
式中，
Figure BDA00032419563400001213
C＝C(X，X)；α＝C-1(Y-eμ*)；
Figure BDA00032419563400001214
e表示长度为n的单位向量。
其中，高斯过程模型的超参数Θ＝{l1，l2，...，lk，...，ld，η2}可通过最大化边缘似然函数(MLE)求得，即最小化负对数边缘似然函数
Figure BDA0003241956340000131
Figure BDA0003241956340000132
Figure BDA0003241956340000133
将低精度训练样本集
Figure BDA0003241956340000134
代入高斯过程模型进行建模，假定低精度高斯过程模型超参数为
Figure BDA0003241956340000135
Figure BDA0003241956340000136
表示低精度高斯过程模型的过程方差；l1，k表示当低精度输入值从x1，k到x′1，k时低精度输出值y1的变化速度，l1，k越大，速度越慢，k表示向量的第k个维度，共有d维。根据最大化边缘似然函数MLE得到超参数Θ1获得训练完成后的网壳固有频率的低精度高斯过程模型y1：
Figure BDA0003241956340000137
式中，
Figure BDA0003241956340000138
表示高斯过程，
Figure BDA0003241956340000139
和
Figure BDA00032419563400001310
分别表示低精度高斯过程模型的预测值和预测方差，可写为：
Figure BDA00032419563400001311
Figure BDA00032419563400001312
式中，
Figure BDA00032419563400001313
C1＝C(X1，X1)是样本点集X1的协方差函数矩阵，表示低精度样本点之间的相关性；
Figure BDA00032419563400001314
表示待预测点x*与低精度输入值X1之间的相关性，x1，i(i＝1，2，...，n1)表示组成低精度输入值X1的第i个样本点；
Figure BDA00032419563400001315
e1表示长度为n1的单位向量。
C.利用所述的高精度训练样本集建立差值高斯过程模型
C1.假设最终预测值可由所述低精度高斯过程模型的预测值和一个高斯误差线性表示：
Figure BDA0003241956340000141
其中，y1(x)是网壳固有频率低精度模型的真实值；y2(x)是网壳固有频率高精度模的型真实值；δ2(x)表示高斯过程误差；⊥表示相互独立；ρ1是一个比例系数，通过最大化边缘似然函数(MLE)获得；
C2.考虑步骤B中所述网壳固有频率的低精度高斯过程模型的预测值误差，其预测值与真实值之间的关系可表示为：
Figure BDA0003241956340000142
其中，ε1(x)表示低精度高斯过程模型的预测值误差，对于两个不同的输入值，ε1(x)是相互独立的，即对任何x≠x′，有Cov(ε(x)，ε(x′))＝0；y1(x)服从正态分布：
Figure BDA0003241956340000143
考虑了所述低精度高斯过程模型在所述高精度输入值X2处的预测误差ε1(x)，不需要满足嵌套样本数据
Figure BDA0003241956340000144
的限制条件。
C3.联立步骤Cl和步骤C2中公式，可将y2(x)与ρ1y1(x)的差值δ2(x)表示为：
Figure BDA0003241956340000145
其中，δ2(x)服从正态分布：
Figure BDA0003241956340000146
C4.将步骤C3中所述的
Figure BDA0003241956340000147
记为δ′2，将高精度训练样本集
Figure BDA0003241956340000148
代入高斯过程模型的建立过程，假定差值高斯过程模型超参数为
Figure BDA0003241956340000149
Figure BDA00032419563400001410
表示差值高斯过程模型的过程方差；l2，k表示当高精度输入值从x2，k到x′2，k时差值δ′2的变化速度，l2，k越大，速度越慢，k表示向量的第k个维度，共有d维。根据MLE得到超参数Θ2，获得训练完成后的差值高斯过程模型：
Figure BDA00032419563400001411
式中，
Figure BDA0003241956340000151
表示高斯过程，
Figure BDA0003241956340000152
和
Figure BDA0003241956340000153
分别表示差值高斯过程模型的预测值和预测方差，可写为：
Figure BDA0003241956340000154
Figure BDA0003241956340000155
式中，
Figure BDA0003241956340000156
C2＝C(X2，X2)是X2的协方差函数矩阵；
Figure BDA0003241956340000157
表示低精度高斯过程模型在X2处的预测方差；
Figure BDA0003241956340000158
表示待预测点x*与X2之间的相关性，x2，i(i＝1，2，...，n2)表示组成高精度输入值X2的第i个样本点；
Figure BDA0003241956340000159
e2表示长度为n2的单位向量。
D.利用低精度高斯过程模型和差值高斯过程模型共同构成广义协同高斯过程模型；具体为：
将步骤B和步骤C所述的网壳固有频率的低精度高斯过程模型和差值高斯过程模型线性组合，共同构成广义协同高斯过程模型，得到最终单层球面壳固有频率的高精度预测值
Figure BDA00032419563400001510
和预测方差
Figure BDA00032419563400001511
Figure BDA00032419563400001512
Figure BDA00032419563400001513
E.通过步骤D构建的单层球面壳固有频率的广义协同高斯过程模型解析地计算该网壳固有频率的均值和方差；
E1.根据统计原理，均值
Figure BDA00032419563400001514
和方差
Figure BDA00032419563400001515
的表达式为：
Figure BDA00032419563400001516
Figure BDA0003241956340000161
式中，p(·)表示概率密度函数；x表示输入值；
Figure BDA0003241956340000162
表示预测值，vy表示预测方差。
通过将协方差函数分离为正态分布密度函数的形式，可将步骤D中所述
Figure BDA0003241956340000163
和
Figure BDA0003241956340000164
的不可积分形式转化为可积分形式：
Figure BDA0003241956340000165
Figure BDA0003241956340000166
式中，
Figure BDA0003241956340000167
α1(2)，i(j)表示α1(2)的第i(j)个元素；Aji表示
Figure BDA0003241956340000168
第j行第i列的元素；Bji表示
Figure BDA0003241956340000169
第j行第i列的元素；
Figure BDA00032419563400001610
Figure BDA00032419563400001611
p1(2)，i表示p1(2)的第i个元素；q1＝C1e1，q2＝Λ2e2，q1(2)，j表示q1(2)的第j个元素；d表示向量的维度；
Figure BDA00032419563400001612
xk表示x*的第k个元素；
Figure BDA00032419563400001613
表示样本点集X1(2)第i(j)行、第k列的元素。
E2.将步骤E1中所述的
Figure BDA0003241956340000171
和
Figure BDA0003241956340000172
代入
Figure BDA0003241956340000173
和
Figure BDA0003241956340000174
和，基于广义协同高斯过程模型，将网壳固有频率均值和方差的高维积分转化为一维积分：
Figure BDA0003241956340000175
Figure BDA0003241956340000176
Figure BDA0003241956340000181
E3.E2中所述的一维积分
Figure BDA0003241956340000182
和
Figure BDA0003241956340000183
可统一地表达为：
Figure BDA0003241956340000184
式中，
Figure BDA0003241956340000185
表示均值为ζ，方差为
Figure BDA0003241956340000186
的正态概率密度函数；p(xk)表示参数xk的概率密度函数。
当参数为正态分布或均匀分布时，通过解析地计算一维积分，最终可以直接得到均值和方差的解析结果，量化结构不确定性大小，解析结果为：
Figure BDA0003241956340000191
式中，
Figure BDA0003241956340000192
表示参数xk服从均值为ξ，方差为θ2的正态分布；
Figure BDA0003241956340000193
表示参数xk服从上下限分别为
Figure BDA0003241956340000194
和x的均匀分布；Φ(·)表示标准正态概率分布函数。
在本算例中，得到该网壳结构固有频率的均值和方差的解析结果，并与MCS法对比，计算结果见表1。由表1可知，GC-GPM法与MCS法的计算结果非常吻合，均值的最大误差仅为0.0070％，方差的最大误差仅为0.6295％。因此，GC-GPM法对不确定性量化具有较高的计算精度。
表1 GC-GPM前5阶固有频率统计矩的计算结果及相对误差
Figure BDA0003241956340000196
为进一步证明本发明所提解析方法的优势，将其与传统高斯过程模型(GPM)方法进行对比。分别采用15、25、35、45个高精度样本数据建立GPM，并计算该网壳前5阶固有频率的均值和方差，同时与MCS法对比，计算结果列于表2。由表1、2可知，要达到与12个高精度样本点的GC-GPM法的计算精度相当，传统GPM法需要35个高精度样本点，而获取高精度样本的计算成本相对较高。以上分析表明，GC-GPM法具有高精度低成本的优势。
表2：GPM前5阶固有频率统计矩的计算结果及相对误差
Figure BDA0003241956340000195
Figure BDA0003241956340000201
注：表中15h、25h、35h、45h分别表示15、25、35、45个高精度样本
本说明书实施案例所述的内容仅仅是对发明构思的实现形式的列举，本发明的保护范围不应当被视为仅限于实施案例所陈述的具体形式，本发明的保护范围也包括本领域技术人员根据本发明构思所能够想到的等同技术手段。

Claims (4)

1.一种基于广义协同高斯过程模型的结构不确定性量化解析方法，其特征在于，包括以下步骤：
A.获得高、低精度的训练样本集
A1.采样概率分布服从正态或均匀分布的结构不确定性参数，对于结构不确定性参数的概率分布不服从正态或均匀分布时，根据概率相等的原则将其转化为服从正态或均匀分布的参数ui后再采样，转化表达式如下：
Figure FDA0003241956330000011
其中，
Figure FDA0003241956330000012
是参数xi的概率分布函数，G表示正态分布或均匀分布函数；然后通过拉丁超立方抽样获得结构不确定性参数的两组样本点，数量分别为n1和n2，将n1个样本点作为低精度输入值X1，n2个样本点作为高精度输入值X2，其中，n1＞n2；
A2.在ANSYS中建立低精度和高精度的结构有限元模型，分别利用所述的低精度输入值X1和高精度输入值X2进行对应的有限元模型的循环计算，得到低精度观测值Y1和高精度观测值Y2，分别组成低精度训练样本集
Figure FDA0003241956330000013
和高精度训练样本集
Figure FDA0003241956330000014
B.利用所述的低精度训练样本集建立低精度高斯过程模型；具体为：
将低精度训练样本集
Figure FDA0003241956330000015
代入高斯过程模型进行建模，假定低精度高斯过程模型超参数为
Figure FDA0003241956330000016
Figure FDA0003241956330000017
表示低精度高斯过程模型的过程方差；l1，k表示当低精度输入值从x1，k到x′1，k时低精度输出值y1的变化速度，l1，k越大，速度越慢，k表示向量的第k个维度，共有d维。根据最大化边缘似然函数MLE得到超参数Θ1，获得训练完成后的低精度高斯过程模型y1：
Figure FDA0003241956330000018
式中，
Figure FDA0003241956330000019
表示高斯过程，
Figure FDA00032419563300000110
和
Figure FDA00032419563300000111
分别表示低精度高斯过程模型的预测值和预测方差，可写为：
Figure FDA0003241956330000021
Figure FDA0003241956330000022
式中，
Figure FDA0003241956330000023
C1＝C(X1，X1)是样本点集X1的协方差函数矩阵，表示低精度样本点之间的相关性；
Figure FDA0003241956330000024
表示待预测点x*与低精度输入值X1之间的相关性，x1，i(i＝1，2，...，n1)表示组成低精度输入值X1的第i个样本点；
Figure FDA0003241956330000025
e1表示长度为n1的单位向量。
C.利用所述的高精度训练样本集建立差值高斯过程模型
C1.假设最终预测值可由所述低精度高斯过程模型的预测值和一个高斯误差线性表示：
Figure FDA0003241956330000026
其中，y1(x)是低精度模型的真实值；y2(x)是高精度模型的真实值；δ2(x)表示高斯过程误差；⊥表示相互独立；ρ1是一个比例系数，通过最大化边缘似然函数(MLE)获得；
C2.考虑步骤B中所述低精度高斯过程模型的预测值误差，预测值与真实值之间的关系可表示为：
Figure FDA0003241956330000027
其中，ε1(x)表示低精度高斯过程模型的预测值误差，对于两个不同的输入值，ε1(x)是相互独立的，即对任何x≠x′，有Cov(ε(x)，ε(x′))＝0；y1(x)服从正态分布：
Figure FDA0003241956330000028
C3.联立步骤C1和步骤C2中的公式，可将y2(x)与ρ1y1(x)的差值δ2(x)表示为：
Figure FDA0003241956330000031
其中，δ2(x)服从正态分布：
Figure FDA0003241956330000032
C4.将步骤C3中所述的
Figure FDA0003241956330000033
记为δ′2，将高精度训练样本集
Figure FDA0003241956330000034
代入高斯过程模型的建立过程，假定差值高斯过程模型超参数为
Figure FDA0003241956330000035
Figure FDA0003241956330000036
表示差值高斯过程模型的过程方差；l2，k表示当高精度输入值从x2，k到x′2，k时差值δ′2的变化速度，l2，k越大，速度越慢，k表示向量的第k个维度，共有d维。根据MLE得到超参数Θ2，获得训练完成后的差值高斯过程模型：
Figure FDA0003241956330000037
式中，
Figure FDA0003241956330000038
表示高斯过程，
Figure FDA0003241956330000039
和
Figure FDA00032419563300000310
分别表示差值高斯过程模型的预测值和预测方差，可写为：
Figure FDA00032419563300000311
Figure FDA00032419563300000312
式中，
Figure FDA00032419563300000313
C2＝C(X2，X2)是X2的协方差函数矩阵；
Figure FDA00032419563300000314
表示低精度高斯过程模型在X2处的预测方差；
Figure FDA00032419563300000315
表示待预测点x*与X2之间的相关性，x2，i(i＝1，2，...，n2)表示组成高精度输入值X2的第i个样本占.
Figure FDA00032419563300000316
e2表示长度为n2的单位向量。
D.利用低精度高斯过程模型和差值高斯过程模型共同构成广义协同高斯过程模型；具体为：
将步骤B和步骤C所述的低精度高斯过程模型和差值高斯过程模型线性组合，共同构成广义协同高斯过程模型，得到最终的高精度预测值
Figure FDA00032419563300000317
和预测方差
Figure FDA00032419563300000318
Figure FDA0003241956330000041
Figure FDA0003241956330000042
E.通过步骤D构建的广义协同高斯过程模型解析地计算结构响应的均值和方差；
E1.根据统计原理，均值
Figure FDA0003241956330000043
和方差
Figure FDA0003241956330000044
的表达式为：
Figure FDA0003241956330000045
Figure FDA0003241956330000046
式中，p(·)表示概率密度函数；x表示输入值；
Figure FDA00032419563300000411
表示预测值，vy表示预测方差。
通过将协方差函数分离为正态分布概率密度函数的形式，可将步骤D中所述
Figure FDA0003241956330000047
和
Figure FDA0003241956330000048
的不可积分形式转化为可积分形式：
Figure FDA0003241956330000049
Figure FDA00032419563300000410
式中，
Figure FDA0003241956330000051
α1(2)，i(j)表示α1(2)的第i(j)个元素；Aji表示
Figure FDA0003241956330000052
第j行第i列的元素；Bji表示
Figure FDA00032419563300000512
第j行第i列的元素；
Figure FDA0003241956330000053
Figure FDA0003241956330000054
p1(2)，i表示p1(2)的第i个元素；q1＝C1e1，q2＝Λ2e2，q1(2)，j表示q1(2)的第j个元素；d表示向量的维度；
Figure FDA0003241956330000055
xk表示x*的第k个元素；
Figure FDA0003241956330000056
表示样本点集X1(2)第i(j)行、第k列的元素。
E2.将步骤E1中所述的
Figure FDA0003241956330000057
和
Figure FDA00032419563300000511
代入
Figure FDA0003241956330000058
和
Figure FDA0003241956330000059
中，基于广义协同高斯过程模型，将均值和方差的高维积分转化为一维积分：
Figure FDA00032419563300000510
Figure FDA0003241956330000061
E3.E2中所述的一维积分
Figure FDA0003241956330000062
和
Figure FDA0003241956330000063
可统一地表达为：
Figure FDA0003241956330000064
式中，
Figure FDA0003241956330000065
表示均值为ζ，方差为
Figure FDA0003241956330000066
的正态分布概率密度函数；p(xk)表示参数xk的概率密度函数。
当参数为正态分布或均匀分布时，通过解析地计算一维积分，最终可以直接得到均值和方差的解析结果，量化结构不确定性大小，解析结果为：
Figure FDA0003241956330000071
式中，
Figure FDA0003241956330000072
表示参数xk服从均值为ξ，方差为θ2的正态分布；
Figure FDA0003241956330000073
表示参数xk服从上下限分别为
Figure FDA0003241956330000074
和x的均匀分布；Φ(·)表示标准正态概率分布函数。
2.根据权利要求1所述的基于广义协同高斯过程模型的结构不确定性量化解析方法，其特征在于，在步骤C中，考虑了所述低精度高斯过程模型在所述高精度输入值X2处的预测误差ε1(x)，不需要满足嵌套样本数据
Figure FDA0003241956330000075
的限制条件。
3.根据权利要求1所述的基于广义协同高斯过程模型的结构不确定性量化解析方法，其特征在于，步骤A2中，根据结构的复杂程度，在ANSYS能够计算的条件下，将结构几何模型划分较少的单元作为低精度有限元模型，将划分单元数为低精度有限元模型的4-5倍作为高精度有限元模型；或者将复杂原始模型作为高精度有限元模型，通过结构受力方式及自重等相同进行等效简化的简单模型作为低精度有限元模型。通过低精度模型计算得到的低精度观测值与真实值偏差较大，因此称为低精度模型。
4.根据权利要求1所述的基于广义协同高斯过程模型的结构不确定性量化解析方法，其特征在于，步骤B和步骤C4中，高斯过程模型的建立过程包括以下步骤：
高斯过程模型表达式如下：
Figure FDA0003241956330000076
式中，
Figure FDA0003241956330000081
表示高斯过程；m(x)为均值函数，采用常数形式μ；C(x，x′)为平方指数协方差函数，表示为：
Figure FDA0003241956330000082
式中，x和x′是输入值；xk与x′k分别是x和x′的第k个元素，d表示输入值的维度；高斯过程模型的超参数定义为Θ＝{l1，l2，...，lk，...，ld，η2}，η2表示过程方差；lk表示输出值在输入值从xk到x′k时的变化速度，lk越大，速度越慢；
假设有n个观测值的训练样本集
Figure FDA0003241956330000083
其中，
Figure FDA0003241956330000084
Figure FDA00032419563300000812
根据模型输出Y服从高斯分布的先验假设，则有：
Figure FDA0003241956330000085
同样的，待预测点x*处的预测值y*与Y也服从高斯分布：
Figure FDA0003241956330000086
根据贝叶斯原理，预测值y*的后验分布为：
Figure FDA0003241956330000087
将所述P(Y)和P(Y，y*)代入后验分布的计算公式，可得：
Figure FDA0003241956330000088
其中，
Figure FDA0003241956330000089
Figure FDA00032419563300000810
式中，
Figure FDA00032419563300000813
C＝C(X，X)；α＝C-1(Y-eμ*)；
Figure FDA00032419563300000814
e表示长度为n的单位向量。
其中，高斯过程模型的超参数Θ＝{l1，l2，...，lk，...，ld，η2}可通过最大化边缘似然函数(MLE)求得，即最小化负对数边缘似然函数
Figure FDA00032419563300000811
Figure FDA0003241956330000091
Figure FDA0003241956330000092
CN202111024605.0A 2021-09-01 2021-09-01 一种基于广义协同高斯过程模型的结构不确定性量化解析方法 Active CN113722961B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111024605.0A CN113722961B (zh) 2021-09-01 2021-09-01 一种基于广义协同高斯过程模型的结构不确定性量化解析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111024605.0A CN113722961B (zh) 2021-09-01 2021-09-01 一种基于广义协同高斯过程模型的结构不确定性量化解析方法

Publications (2)

Publication Number Publication Date
CN113722961A true CN113722961A (zh) 2021-11-30
CN113722961B CN113722961B (zh) 2024-02-13

Family

ID=78680815

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111024605.0A Active CN113722961B (zh) 2021-09-01 2021-09-01 一种基于广义协同高斯过程模型的结构不确定性量化解析方法

Country Status (1)

Country Link
CN (1) CN113722961B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114996995A (zh) * 2022-05-20 2022-09-02 华中科技大学 一种超材料隔振单元性能预报方法和系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2437034A2 (en) * 2010-09-30 2012-04-04 Siemens Corporation System and method for conditional multi-output regression for machine condition monitoring
CN106056127A (zh) * 2016-04-07 2016-10-26 江南大学 一种带模型更新的gpr在线软测量方法
CN110197286A (zh) * 2019-05-10 2019-09-03 武汉理工大学 一种基于混合高斯模型和稀疏贝叶斯的主动学习分类方法
WO2021063136A1 (zh) * 2019-09-30 2021-04-08 江苏大学 一种数据驱动的高精度组合导航数据融合方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2437034A2 (en) * 2010-09-30 2012-04-04 Siemens Corporation System and method for conditional multi-output regression for machine condition monitoring
CN106056127A (zh) * 2016-04-07 2016-10-26 江南大学 一种带模型更新的gpr在线软测量方法
CN110197286A (zh) * 2019-05-10 2019-09-03 武汉理工大学 一种基于混合高斯模型和稀疏贝叶斯的主动学习分类方法
WO2021063136A1 (zh) * 2019-09-30 2021-04-08 江苏大学 一种数据驱动的高精度组合导航数据融合方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
万华平;任伟新: "复杂结构不确定性量化分析和敏感性分析的高斯过程模型框架", 第七届全国结构抗振控制与健康监测学术会议 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114996995A (zh) * 2022-05-20 2022-09-02 华中科技大学 一种超材料隔振单元性能预报方法和系统
CN114996995B (zh) * 2022-05-20 2024-07-02 华中科技大学 一种超材料隔振单元性能预报方法和系统

Also Published As

Publication number Publication date
CN113722961B (zh) 2024-02-13

Similar Documents

Publication Publication Date Title
CN105608263B (zh) 一种面向涡轮叶盘结构寿命概率分析的自适应处理方法
CN113627098B (zh) 一种cfd模型确认方法及产品设计方法
CN109726437B (zh) 一种舱门气动载荷等效节点力处理方法
CN112580228B (zh) 一种涡扇发动机混合结构风扇叶片结构设计优化方法
CN109545290B (zh) 一种基于Voronoi分形技术的非晶合金自由体积检测方法
CN109599866B (zh) 一种预测辅助的电力系统状态估计方法
CN113722961A (zh) 一种基于广义协同高斯过程模型的结构不确定性量化解析方法
CN112434380A (zh) 一种飞机机翼重量估算方法
CN105046041A (zh) 一种压气机可控扩散叶型稳健优化设计方法
CN115495965A (zh) 复杂航空结构在混合不确定性下时变可靠性的分析方法
CN112307536A (zh) 一种大坝渗流参数反演方法
CN106407601B (zh) 一种基于数据挖掘技术的气动特性数据处理方法
Chen et al. Towards high-accuracy deep learning inference of compressible turbulent flows over aerofoils
CN108694472B (zh) 预测误差极值分析方法、装置、计算机设备和存储介质
CN111931272B (zh) 任意尺度边坡等精度安全系数计算方法及网格划分方法
CN107808021B (zh) 基于cfd的流体设备的阻力计算方法
CN110348701B (zh) 一种水库群防洪调度风险传递规律分析方法
CN116933386A (zh) 一种基于mcok代理模型的飞行器气动数据融合方法
CN112287484A (zh) 一种基于矢量代理模型的复杂工程系统可靠性设计方法
CN115293069B (zh) 一种用于飞行器外流场仿真控制参数智能优化的系统
Chau et al. A structural reliability analysis method based on radial basis function
CN116011071A (zh) 基于主动学习的空中造楼机结构可靠性分析方法及系统
CN106980675B (zh) 一种高效的桥梁结构健康预警系统
CN111754027B (zh) 一种基于概率语言术语集的pm2.5排放权分配方法
CN106342305B (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