一种重载铁路桥梁有限元模型的位移预测方法

本发明属于重载铁路有限元模型修正领域，具体涉及一种重载铁路桥梁有限元模型的修正方法。

背景技术：

1、近年来，随着我国经济快速发展，交通运输量极大地提高，作为交通运输关键点的桥梁在其中承担着极其重要的角色。为满足日益庞大的货载需要，因此采用建设重载铁路桥梁的方法，以期望解决国家运量大的难题。有限元模型精确与否是当前土建领域极其重要的一环，与项目实施效率、质量息息相关，虽然基于普通桥梁有限元模型修正经验，重载铁路桥梁模型修正有一定实践基础，但是其轴重大、牵引质量大的特点导致该领域仍需进一步探索。

2、一般桥梁有限元模型修正技术基本成熟。已有人针对大跨度连续刚构桥有限元模型修正问题，提出一种基于改进粒子群算法优化bp神经网络的有限元模型修正方法。该方法优化bp神经网络的代理模型使其具有更高的效率和拟合精度，且有效预测了实际钢构桥挠度变化。已有人基于径向神经网络，结合遗传算法，通过精准有限元模型建立分析，完成了对大跨度悬索桥结构行为进行精确预测与分析。已有人将神经网络方法结合matlab用于处理大跨度预应力混凝土连续粱拆施工控制中的箱粱线形控制的误差预测和调整，通过精确有限元模型修正，实现预测和实践的统一。综上所述，有限元模型修正方法基本成熟且进行了一定实际运用，但是基于桥梁位移形变相关方面尚且不足。

3、重载铁路桥梁比之普通桥梁荷载更大，对于力学结构更为严谨，因此一般桥梁有限元模型修正方法尚不适配重载铁路桥梁。此外相关领域基于充分融合观测信息和先验信息的贝叶斯方法，构造包含动态权重系数的目标函数，结合能替代有限元模型进行计算的响应面方法对重载铁路桥梁进行了一定探索。该方法基于重载铁路黄河特大桥实例，获取了主梁控制截面竖向位移相对误差1％以下的成果，但桥梁跨中等位移误差仍然较大，故在重载铁路桥梁位移领域仍需进一步研究，相关模型修正方法仍需研发。

技术实现思路

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

2、本发明的目的在于提供一种重载铁路桥梁有限元模型的修正方法，本方法能够有效解决现有技术中存在的问题，尤其是现有的重载铁路桥梁有限元模型修正方法的修正效率和修正精度较低的问题。

3、本发明提供了一种重载铁路桥梁有限元模型的位移预测方法，包括以下步骤：

4、s1，建立目标桥梁的实体有限元模型；

5、s2，确定有限元模型变量参数x，建立目标函数q(x)；

6、s3，采用响应面法构建响应面模型；

7、s4，根据步骤s3的响应面模型，采用多岛遗传算法对所述有限元模型变量参数进行修正；

8、s5，基于s4的修正后的有限元模型变量参数修正所述实体有限元模型，导入有限元分析软件进行不同荷载位移预测。

9、进一步地，步骤s1中，根据工程数据建立ansys重载铁路桥梁模型，并且计算出全桥关键点位移数据。

10、进一步地，步骤s2所述变量参数x＝[e,ρ,μ]，其中，e为混凝土弹性模量，ρ为混凝土密度，μ为泊松比。

11、进一步地，步骤s2所述目标函数公式：

12、

13、其中zq为有限元模型计算出的桥梁q点位移，为实测q点位移，s表示位移点个数。

14、进一步地，采用响应面法构建响应面模型具体包括以下步骤：

15、s31，基于s1步骤选取i个关键点的初始参数变量xi＝[ei,ρi,μi](i＝1,2，...，n)，以及响应面目标函数y的样本点数值y’＝{y1,y2,...,yn}；

16、s32，根据响应面法原理采用最小二乘法回归分析获得二次多项式的回归系数，得到目标函数y即响应面模型；

17、s33，采用响应面目标函数y替代q(x)，通过数学优化方法对响应面目标函数进行优化，得响应面目标函数y与误差项函数最小时的参数变量取值，基于此参数变量取值的二次多项式即为目标函数q(x)的响应面模型；

18、s34，根据r2分析，若r2>0.95，则转入s35，否则进入步骤s31；

19、s35，输出响应面模型y。

20、进一步地，根据权利要求5所述的一种重载铁路桥梁有限元模型的位移预测方法，其特征在于，步骤s31提出的响应面目标函数表示为：

21、

22、其中ei、ρi、μi表示第i次下参数取值，其中e为混凝土弹性模量，ρ为混凝土密度，μ为泊松比；αa(a＝0，1，...，9)是待定系数。

23、进一步地，进一步地，将选取n个样本点数值y’＝{y1,y2,...,yn}，代入公式2中，得误差项函数s的计算公式：

24、

25、其中yt为样本点数值，为含有待定系数的样本点数值，t＝0，1，...，9。

26、进一步地，根据最小二乘法原理，对误差项函数s求最小值，即得到函数待定最小二乘估计值，求出公式3中待定系数值。

27、进一步地，步骤s34提出的r2表示为：

28、

29、其中y为实际位移值，为平均位移值，为位移预测值。

30、进一步地，基于多岛遗传算法对该有限元模型参数进行模型修正包括以下步骤：

31、s41，初始化群体，将一个大的种群分成若干个子种群，且每个子种群为岛，在每个岛上运用miga算法进行子种群进化；

32、s42，明确群体上每个个体适应度值；

33、s43，评价个体适应度值同时选择将进入下一代的个体；

34、s44，按概率进行交叉操作；

35、s45，按概率进行突变操作；

36、s46，对子代进行规则筛，若没有满足停止条件，则转步骤s42，否则进入步骤s47；

37、s47，输出种群中适应度值最优子代作为问题最优解。

38、进一步地，基于步骤s47输出值建立的有限元模型即为修改后的精确模型，对该模型输入不同荷载即可完成该桥梁全桥位移预测。

39、本发明的优点在于：本发明公开了一种重载铁路桥梁有限元模型的位移预测方法，根据多岛遗传算法结合响应面法的基础上完成模型参数修正，再通过输入不同荷载达到对全桥位移的预测。此方法对于有限元计算次数少，而且能够得到高精度的修正结果，避免了大型有限元模型修正过程中反复迭代计算耗时过长的弊端，提高了大型有限元模型修正工作的效率；基于多岛遗传方法，克服了一般的迭代方法容易陷入局部极小的陷阱而出现“死循环”现象，使迭代顺利进行；此外，响应面方法作为一种试验设计与数理统计相结合的数学分析方法，运用于大型复杂桥梁结构的有限元模型修正中是可行的，能够同时保证较高的修正效率和修正精度。同时该方法对于重载铁路桥梁位移预测精度起到了至关重要的重要。该发明为重载铁路桥梁模型修正提供了一定借鉴意义，推动该领域进一步发展。

技术特征：

1.一种重载铁路桥梁有限元模型的位移预测方法，其特征在于，包括以下步骤：

2.根据权利要求1所述的一种重载铁路桥梁有限元模型的位移预测方法，其特征在于，步骤s1中，根据工程数据建立ansys重载铁路桥梁模型，并且计算出全桥关键点位移数据。

3.根据权利要求1所述的一种重载铁路桥梁有限元模型的位移预测方法，其特征在于，步骤s2所述变量参数x＝[e,ρ,μ]，其中，e为混凝土弹性模量，ρ为混凝土密度，μ为泊松比。

4.根据权利要求1所述的一种重载铁路桥梁有限元模型的位移预测方法，其特征在于，步骤s2所述目标函数公式：

5.根据权利要求3所述的一种重载铁路桥梁有限元模型的位移预测方法，其特征在于，采用响应面法构建响应面模型具体包括以下步骤：

6.根据权利要求5所述的一种重载铁路桥梁有限元模型的位移预测方法，其特征在于，步骤s31提出的响应面目标函数表示为：

7.根据权利要求6所述的一种重载铁路桥梁有限元模型的位移预测方法，其特征在于，将选取n个样本点数值y’＝{y1,y2,...,yn}，代入公式2中，得误差项函数s的计算公式：

8.根据权利要求7所述的一种重载铁路桥梁有限元模型的位移预测方法,其特征在于，根据最小二乘法原理，对误差项函数s求最小值，即得到函数待定最小二乘估计值，求出公式3中待定系数值。

9.根据权利要求5所述的一种重载铁路桥梁有限元模型的位移预测方法，其特征在于，步骤s34提出的r2表示为：

10.根据权利要求1所述的一种重载铁路桥梁有限元模型的位移预测方法，其特征在于，基于多岛遗传算法对该有限元模型参数进行模型修正包括以下步骤：

技术总结
本发明公开了一种重载铁路桥梁有限元模型的位移预测方法，根据多岛遗传算法结合响应面法的基础上完成模型参数修正，再通过输入不同荷载达到对全桥位移的预测。此方法对于有限元计算次数少，而且能够得到高精度的修正结果，避免了大型有限元模型修正过程中反复迭代计算耗时过长的弊端，提高了大型有限元模型修正工作的效率；基于多岛遗传方法，克服了一般的迭代方法容易陷入局部极小的陷阱而出现“死循环”现象，使迭代顺利进行；此外，响应面方法作为一种试验设计与数理统计相结合的数学分析方法，运用于大型复杂桥梁结构的有限元模型修正中是可行的，能够同时保证较高的修正效率和修正精度。同时该方法对于重载铁路桥梁位移预测精度起到了至关重要的重要。

技术研发人员：徐向阳,杨浩,向阳
受保护的技术使用者：苏州大学
技术研发日：
技术公布日：2024/2/29