基于kriging代理模型的桥梁数字孪生体构建方法与流程

本发明属于桥梁结构监测领域，涉及一种桥梁运营安全监测数字孪生体构建方法，具体涉及一种基于kriging代理模型的桥梁运营安全监测数字孪生体构建方法。

背景技术：

1、实际桥梁结构其结构参数受运营环境、施工误差等因素影响，其结果往往和设计参数值存在差异，导致在采用设计参数值进行结构分析过程中，桥梁结构响应与实际传感器监测分析值产生误差，为了保证数字孪生体的准确构建，需要构建实桥在运营环境中所受影响在虚拟模型中的实现方法。

2、针对上述问题，在运营环境下有监测系统的桥梁结构，引入桥梁结构代理模型概念，以结构随机有限元为核心，构建桥梁结构概率有限元基准模型算法以建立能够表征结构刚度参数偏移概率分布。考虑结构参数之间灵敏度问题，采用基于矩阵摄动法的代理模型灵敏度计算分析方法，为修正参数的选择提供依据。通过基于贝叶斯后验分布的最优参数值以及考虑协方差的最优估计共同构建代理模型目标函数，以实际桥梁的历史监测数据分析计算频率为频率先验分布，常规设计参数值及其经验分布范围为修正参数先验分布，经过多次迭代从而实现结构刚度参数偏移概率分布算法，由于修正参数组合试验样本较多且有限元分析计算时间较长，导致代理模型的计算量较大，基于桥梁结构监测数字孪生特点，基于正交试验设计法的方法设计概率种群，构建修正参数与频率之间的kriging模型以提高效率，结合结构刚度参数偏移概率分布算法，从而实现桥梁结构代理模型的完整构建。

技术实现思路

1、针对上述既有结构监测分析值存在的不足，本发明提供了一种基于kriging代理模型的桥梁数字孪生体构建方法，该方法可大幅度提高计算效率和准确率。

2、本发明的目的是通过以下技术方案实现的：

3、一种基于kriging代理模型的桥梁数字孪生体构建方法，包括如下步骤：

4、步骤一：根据所监测的桥梁结构采集到的实时车辆数据，利用两个正态分布函数加权和的双峰型概率分布对全年累计车辆数据按类统计，基于证据理论的多源信息融合，结合数值仿真和监测数据分析，得到结构响应特征和预测结果，构建桥梁安全阈值等级标准，具体步骤如下：

5、步骤一一：构建两个正态分布函数加权和的双峰型概率分布：

6、

7、

8、式中，σ1、σ2为两峰值各自正态方差；μ1、μ2为两峰值各自正态均值；α为第一正态分布权重；

9、步骤一二：基于证据理论的多源信息融合，结合数值仿真和监测数据分析，构建桥梁安全阈值等级标准：

10、(1)根据桥梁结构数值仿真和监测数据，确定辨识框架θ，在辨识框架θ下，确定n条证据的基本可信度分配函数m1，m2，m3···mn，若集函数m：2θ→[0,1]，则满足：

11、

12、

13、式中，m为框架θ上的基本可信度分配；m(a)为a的基本可信数；

14、(2)明确每个基本信度分配函数的焦元，并计算冲突系数k：

15、如果m：2θ→[0,1]为框架θ上的基本可信度分配，则：

16、

17、若m(a)＞0，则a为信度函数bel的焦元；

18、假设bel1和bel2是同一个识别框架θ下的两个信度函数，则冲击系数k：

19、

20、式中，m1和m2分别是他们所对应的基本可信度分配，焦元分别为a1，a2，a3···ak，k为焦元的个数；

21、(3)利用dempster合成法则对n条证据的基本可信度进行合成：

22、

23、(4)根据桥梁运营期结构特点将评估阈值等级分为三级：

24、1)第一级：考虑日常运营期在车流量最大阶段时桥梁运营安全，即采用日常运营中的高峰期评估值；

25、2)第二级：考虑全年运营期在重车荷载作用下的桥梁运营安全，即采用全年中的重车评估值；

26、3)第三级：基于物理模型的桥梁结构监测评估方法，即采用ansys计算出来的结构评估值；

27、步骤二：根据步骤一中得到的桥梁安全阈值对应情况，构建交通荷载概率分布算法，计算不同安全阈值等级对应的桥梁结构交通荷载概率分布；利用所构建的交通荷载概率分布算法和采用蒙特卡洛模拟构建的桥梁车辆荷载计算模型，计算对应安全阈值下车辆概率分布的桥梁瞬态分析响应；具体步骤如下：

28、步骤二一：假设某类车辆荷载车重数据分布服从概率密度函数为f(x|θ1,θ2,...,θm)，其中x为车重数据，m为荷载个数，θ为函数参数样本空间，θ1,θ2,...,θm∈θ为函数参数个数，则车重数据的联合概率密度函数为：

29、

30、若有使得

31、

32、则得到极大似然估计量为即得到该车重数据概率密度函数参数；

33、设某一类车辆荷载样本的分布函数为fwn(x)，根据经验确定车辆荷载概率分布函数为fw(x)，设假设h0为经验确定的车辆荷载概率分布和样本概率分布一致，由此可得检验统计量d定义如下：

34、

35、步骤二二：基于步骤二一中获得的若干车辆荷载概率模型分布，根据各车型荷载模型、车辆间距和车轴等信息参数，采用蒙特卡洛模拟产生随机车辆，按照车流量在高峰期间各车型比例随机产生，模式如下所示：

36、将桥梁点位纵向分为k段，分别位于各跨中边跨位置，横向根据车道数量分为n段，以模拟车辆来回过程，则桥梁荷载布设施加位置如下式所示：

37、

38、由此得到车辆荷载概率分布通过蒙特卡洛模拟产生车辆概率分布荷载为：

39、

40、式中，z1,z2,...,zg为模拟车辆荷载质量；g为模拟车辆荷载质量个数；t1,t2,...,tg为形状参数；

41、进而通过模拟车辆概率分布荷载m在x矩阵中的随机分散分布，得到概率分布荷载工况，其中对应分布荷载位置如下所示：

42、pij＝[zh th]t

43、式中，zh为第h处的模拟车辆荷载质量，th为第h处的形状参数；

44、其他未分布荷载值为0，实现实际车辆车流的随机性位置分布，通过全年车流车速分析，确定对应时间段车辆数概率分布，根据车辆位置进行足量分析计算，得到有限元车辆荷载仿真空间分布，进而实现计算不同安全阈值等级对应的桥梁结构交通荷载概率分布和对应安全阈值下车辆概率分布的桥梁瞬态分析响应；

45、步骤三：依据步骤二中计算的不同安全阈值等级下交通荷载概率分布和桥梁瞬态分析响应结果，采用pod算法对原始样本数据进行降维，利用pce模型对kriging模型进行改进重构，结合实桥响应，构建实桥响应与荷载组合的基于代理模型的实桥数字孪生体pod-pce-kriging模型；提出概率加权混合加点准则，通过k-fold交叉验证方法计算误差，对kriging代理模型进行修正优化设计，实现基于代理模型的实桥数字孪生体kriging模型的建立，具体步骤如下：

46、步骤三一：构建pod-pce-kriging模型代替结构分析计算，具体步骤如下：

47、(1)通过pod算法对响应数据空间进行降维；

48、1)设n维场量空间{xn∈ω}的一组“最优”正交基{u1,u2,u3···un}，假设xn可由若干正交基近似，则：

49、

50、式中，ui为xn的特征向量或基向量；αi为ui的相关系数；r为降维后的维度数；

51、2)对桥梁结构响应数据进行中心化处理，则得到新的中心化响应数据集：

52、

53、式中，为中心化样本数据集；为针对一维度的样本均值；

54、3)对新得到的中心化响应数据集进行奇异值分解：

55、

56、式中，λ＝{λ1,λ2,λ3···λn}为协方差矩阵由大到小排列的特征值组成的向量，λ为u的特征值，n为特征值的个数；

57、(2)利用pce模型对kriging模型进行改进重构：

58、1)kriging代理模型中某一点未知点x的预测响应为：

59、

60、式中，f(x)为回归模型；μ为回归系数；g(x)为代理模型的估计误差；

61、2)采用pce模型替代原本kriging模型中常用的回归部分，对于任一样本，表示为以下pce：

62、

63、式中，{βa,a∈nn}为多变量标准正交多项式fa(x)的系数；其中fai(i)为第i个变量的第ai阶多项式；

64、对于s个响应数据，回归部分f(β,x)表示为：

65、f(β,x)＝[f(β,x1),f(β,x2),f(β,x3),...,f(β,xs)]t

66、3)g(x)表示平稳随机过程，代表回归模型近似目标函数的系统偏差，其满足以下条件：

67、

68、式中，σ2代表方差；r(θ,xi,xj)是假定的关于超参数的相关系数；

69、4)根据最优线性无偏估计要求，求得：

70、

71、5)通过求解无约束最优化问题得到超参数θ：

72、

73、6)某一未知点x的估计响应值为：

74、

75、r(x)＝[r(θ,x1,x),r(θ,x2,x),r(θ,x3,x)···r(θ,xn,x),]

76、选择相关函数r为：

77、

78、步骤三二：基于概率加权混合加点准则进行代理模型序列优化设计，具体步骤如下：

79、(1)在设计空间中选取少量初始样本x0，循环次数i＝0；

80、(2)进行数值计算或者结构分析，获取样本xi对应的实际响应f(xi)；

81、(3)基于已有样本xi及其响应值f(xi)，运用kriging代理模型技术，构建近似函数

82、(4)采用概率加权混合加点准则的优化算法对代理模型进行优化，得到第i次循环的最优解并计算实际响应值fimin；

83、(5)核查优化是否收敛：若优化结果已满足收敛条件，则输出最优设计{xmin,fmin}；否则i＝i+1，继续步骤(6)；

84、(6)根据局部加点准则可计算第一个候选样本点xinew1为：

85、

86、式中，“arg”表示为自变量与因变量之间的映射关系，返回值为当前模型最小值对应的输入参数；

87、假设d维设计空间中存在一组样本点x＝[x1,x2,x3···xn]，根据这组样本将设计空间分割为n个泰森多边形：

88、ri＝{x∈ω|||x-xi||2＜||x-xj||2}j≠i

89、式中：ω代表整个设计空间；x为空间中任意一向量点；

90、在k-fold交叉验证中，将已有的xi-1随机地分为k组互斥的子集di＝[d1,d2,d3···dn]，样本的预测误差e为：

91、

92、式中，f(di)为di集合中所有样本的实际响应值；是由除di组以外样本构建的代理模型；表示di中样本的预测响应值；

93、将样本的交叉验证误差进行排序，选出预测误差最大的点xc，根据全局加点准则筛选样本，计算第二个候选样本点xinew2：

94、xnew2＝argmax(||x-xc||2)

95、采用概率加权混合加点准则计算x混：

96、

97、(7)将x混加入总样本中，得到集合xi，转至步骤2)，继续进行迭代。

98、相比于现有技术，本发明具有如下优点：

99、1、本发明利用交通荷载概率分布算法，构建交通荷载概率分布及模型，构建实桥响应与荷载组合kriging模型，实现基于监测系统的数字孪生体响应仿真。

100、2、本发明可以有效地提升数字孪生体的构建和更新效率，大幅度提高模型精度以及所生成数字孪生体响应仿真准确率，可以达到所构建pod-pce-kriging模型具备足够精度能代替结构分析过程、提高计算效率和预测精度的目的，适用于运营周期内桥梁与隧道的诊断评估及监测。

101、3、本发明解决了现有基于代理模型的桥梁数字孪生体构建方法对修正参数组合试验样本较多且有限元分析计算时间较长，导致代理模型的计算量较大的问题。