[发明专利]一种基于贝叶斯和GARCH的预应力预测方法

权利要求书

1.一种基于贝叶斯和GARCH的预应力预测方法，其特征在于，所述预应力预测方法包括：

步骤一，选取同规格的预应力构件一和预应力构件二，测量预应力构件一在张拉过程中的应力应变数据，获取一组预应力构件一的应力数据f_i和应变数据ε_i，通过蒙特卡洛方法，若测量精度为90％数据服从正态分布，则对正态分布N(ε_i,0.1f_i/E)进行抽样模拟，得到预应力构件一的应力-应变的先验概率分布P(ε_ik)；

步骤二，选取预定数量测试点，对预应力构件二在张拉过程中的应力-应变数据进行测量，获得对应不同应力的应变集；

步骤三，将步骤一和步骤二中得到的数据作为基础，并根据贝叶斯方法预测预应力构件二后段的应力-应变数据，获得应力-应变的后验概率分布，具体操作为：将步骤二中测得的应变集作为试验样本O_M，根据贝叶斯方法，可得到应力-应变的后验概率分布P'(ε_ik)：

P'(ε_ik)＝N×L(O_M|ε_jk)×P(ε_ik)；

其中，ε_ik为第f_i应力第k次蒙特卡洛试验得到的预应力构件的应变；O_M为对应不同应力的预应力构件的应变集；N为归一化常数；L(O_M|ε_jk)代表给定ε_jk的分布时出现试验样本O_M的概率，满足：

将上述公式代入后验概率分布P'(ε_ik)的公式中，则在第k次蒙特卡洛试验，对应力f_i时，应变ε_ik的后验分布P'(ε_ik)为：

步骤四，将步骤二中测得的应变集作为基础数据，并根据GARCH方法预测预应力构件二后段预设步数以内的应力-应变值；

步骤五，预应力构件二应力-应变预测值的采信：后段预设步数以内采信上述步骤四的预测值为实际预测值，预设步数以上采信步骤三的预测值为实际预测值。

2.根据权利要求1所述的基于贝叶斯和GARCH的预应力预测方法，其特征在于，在所述步骤四中，将步骤二中测得的应变集作为基础数据，设ε_i服从GARCH(m,s)模型，则：

ε_i＝σ_iη_i

其中η_i是均值为0，方差为1的独立分布随机变量序列，对与预测预应力构件二的应力-应变分布可简单取m＝1,s＝1，此时GARCH(1,1)在预设L步内预测为：