一种基于静载试验识别连续梁桥实际刚度的方法技术

本发明专利技术涉及一种基于静载试验识别连续梁桥实际刚度的方法，包括步骤：依照荷载试验加载工况建立桥梁有限元计算模型；计算结构位移，将位移向量按照测量值和未测量值分开，构建目标函数；选择全部单元刚度作为待识别参数；采用遗传算法，设置待识别参数的上下限进行初步优化计算；将遗传算法计算结果作为初始值，采用L‑M算法进行第二轮优化计算；根据测点布置结合灵敏度分析判断刚度识别的有效区域。本发明专利技术依据静载试验数据进行有限元模型修正，不需要额外布置测点和加载方案；不需要进行有限元模型缩聚，可避免矩阵奇异；选用全部单元刚度作为待识别参数，可避免参数选择不合理造成的误差；结合遗传算法和L‑M算法的特点，提高了优化计算的精度。

【技术实现步骤摘要】
一种基于静载试验识别连续梁桥实际刚度的方法
本专利技术属于路桥工程
，具体来说，涉及到一种基于静载试验识别连续梁桥实际刚度的方法。
技术介绍
桥梁结构刚度状况直接关系着结构的受力性能，是结构计算中非常重要的力学参数。目前结构刚度均采用理论计算值。实际工程中，混凝土弹性模量、截面几何尺寸、截面开裂、预应力损失等均会导致实际刚度与理论刚度存在偏差，从而影响结构受力。对于运营桥梁，荷载试验仍是桥梁承载力评定的主要方法。荷载试验主要通过实测值与理论值的比值，即校验系数进行判断。该方法简单直观，能够得出桥梁结构承载力是否满足，但不能掌握桥梁结构的实际刚度状况。有限元模型修正技术，主要依据实测结果对结构力学参数进行识别，能够精准识别桥梁的实际受力状况，为桥梁承载力的评定提供准确的数据。其中静载试验测试数据准确，加载效率高，能够真实反映桥梁的受力状况，但测点较少造成有限元模型修正困难。本专利技术结合静载试验特点，基于有限元模型修正，提出一种有效识别在役桥梁实际刚度的方法。
技术实现思路
本专利技术的目的是解决在役桥梁实际刚度获取的问题，通过桥梁荷载试验，利用有限元模型修正技术，得出识别桥梁结构实际刚度的实用方法，具体过程如下：步骤一：依照荷载试验加载工况建立桥梁有限元计算模型；步骤二：计算结构位移，将位移向量按照测量值和未测量值分开，构建目标函数；步骤三：选择全部单元刚度作为待识别参数；步骤四：采用遗传算法，设置待识别参数的上下限进行初步优化计算；步骤五：将遗传算法计算结果作为初始值，采用L-M算法进行第二轮优化计算；步骤六：根据测点布置结合灵敏度分析判断刚度识别的有效区域。进一步的，步骤一，桥梁有限元计算模型的力学方程为：[k][δ]＝[F](1)式中[k]为结构的总刚度矩阵，[δ]为结构的位移矩阵，[F]为荷载矩阵，当桥梁静载试验涉及多个加载工况，将不同工况组成荷载矩阵，每列表示不同的工况，每行表示不同的节点号。进一步的，步骤二，计算出结构的位移矩阵,具体如下：[δ]＝[F][K]-1(2)计算出的位移矩阵[δ]为模型所有节点的位移，根据测点布置，从计算的位移矩阵中提取出测点的计算值，当实测值较少时，将[δ]按照测量值和未测量值分开表示，则上式表示为：式中δa为测点的位移，δb为未测点的位移，得出：[δa]＝[F][[K]-1a](4)通过计算与实测值的差值，构建出目标函数：[G]＝[δa]-[Ua]＝[F][[K]-1a]-[Ua](5)式中[Ua]为实测位移矩阵，求得一组参数使得公式(5)计算结果接近于零，根据式(5)计算的[G]为位移残差矩阵，行代表各测点，列代表各工况；将矩阵[G]变换为一维向量A，通过计算向量A的模，当|A|越小，说明残差矩阵越接近于零，使得有限元模型修正转换为求|A|最小值问题，采用优化算法进行计算。进一步的，步骤四至五种，根据步骤二确定目标函数，求|A|的最小值，采用遗传算法进行第一轮优化计算，将遗传算法计算得出的各参数值作为初始值，采用L-M算法进一步优化计算，搜索步长dk计算如下：式中Jk为目标函数对各参数的偏导数，采用有限差分法计算Jk，μk为阻尼因子，I为单位矩阵，Fk为目标函数；式中参数μk的计算如下：μk初始值取较小值，μk初值取0.0001～0.00001，ηk为目标函数与迭代点二次函数增量的比值。附图说明图1是某连续梁桥实测挠度曲线图。图2是L-M算法计算结果图。图3是本专利技术提供方法的计算结果图。图4是本专利技术提供方法的流程框图。图5是桥型布置图。图6是跨中断面布置图。图7是挠度测点布置图。图8和图9分别是工况一和工况二的加载方案布置图。图10是具体实施方式中的迭代过程图。图11是具体实施方式中的迭代过程图。具体实施方式下面结合具体的实施例对本专利技术所述的基于静载试验识别连续梁桥实际刚度的方法做进一步说明，但是本专利技术的保护范围并不限于此。本专利技术的目的是解决在役桥梁实际刚度获取的问题，通过桥梁荷载试验，利用有限元模型修正技术，总结出识别桥梁结构实际刚度的实用方法，具体过程如下：步骤一：依照荷载试验加载工况建立桥梁有限元计算模型；步骤二：计算结构位移，将位移向量按照测量值和未测量值分开，构建目标函数；步骤三：选择全部单元刚度作为待识别参数；步骤四：采用遗传算法，设置待识别参数的上下限进行初步优化计算；步骤五：将遗传算法计算结果作为初始值，采用L-M算法进行第二轮优化计算；步骤六：根据测点布置结合灵敏度分析判断刚度识别的有效区域。例如，某连续梁桥荷载试验，在边跨最大正弯矩加载工况下的实测挠度和计算挠度结果见图1。从图1可以看出，第二跨实测挠度与理论计算挠度存在着较大的差别，采用上述实测数据，按照L-M方法进行结构刚度优化，单元刚度初始值采用计算值，优化计算结果见图2。优化后的结果的并不理想，优化后的计算值与实测值偏差为2.7mm，模型优化的收敛性较差。采用本方法进行优化计算，首先采用遗传算法进行计算，设置参数的上下界限，将识别的刚度结果作为初始值，然后采用L-M优化算法再进行优化计算，计算结果见图3。从图3可看出，本方法优化后的计算值与实测值最大偏差0.9mm，模型优化的收敛误差在可接受范围内。说明本方法通过有限元模型修正，计算值与理论值基本吻合，能够有效地识别出桥梁结构的实际刚度。结合图4，说明本方法的实施方式，具体如下：步骤一：依照荷载试验加载工况建立桥梁有限元计算模型；对步骤一的进一步说明，桥梁有限元计算模型是基于位移法的计算原理，其力学方程为：[k][δ]＝[F](1)式中[k]为结构的总刚度矩阵，[δ]为结构的位移矩阵，[F]为荷载矩阵。按照上式编写相关程序代码，桥梁静载试验一般涉及多个加载工况，可将不同工况组成荷载矩阵，每列表示不同的工况，每行表示不同的节点号。步骤二：计算结构位移，将位移向量按照测量值和未测量值分开，构建目标函数；对步骤二的进一步说明，计算出结构的位移矩阵,具体如下：[δ]＝[F][K]-1(2)计算出的位移矩阵[δ]为模型所有节点的位移，实际桥梁荷载试验中不可能对所有位移均进行测试，因此，根据测点布置，从计算的位移矩阵中提取出测点的计算值。由于实测值较少，将[δ]按照测量值和未测量值分开表示，则上式可以表示为：式中δa为测点的位移，δb为未测点的位移，这样便得出：[δa]＝[F][[K]-1a](4)通过计算与实测值的差值，从而构建出目标函数。[G]＝[δa]-[Ua]＝[F][[K]-1a]-[Ua](5)式中[Ua]为实测位移矩阵。有限元模型修正的核心在于求得一组参数使得公式(5)计算结果为零。事实上由于计算模型误差和传感器测试误差等原因，式(5)不可能绝对为零，只是接近于零。根据式(5)计算的[G]为位移残差矩阵，行代表各测点，列代表各工况。将矩阵[G]变换为一维向量A，通过计算向量A的模，当|A|越小，说明残差矩阵越接近于零。这样有限元模型修正便转换为求|A|最小值问题，可采用各种优化算法进行计算。步骤三：选择全部单元刚度作为待识别参数；对步骤三的进一步说明，从式(4)的表达式可以看出，位移的计算公式为[δa]＝[F][[K]-1a]，其影响因素主要有[F]和[K]。对于桥梁静载试验，加载重量是经过准确称重的，即矩阵[F]是确定。则本文档来自技高网...

【技术保护点】
一种基于静载试验识别连续梁桥实际刚度的方法，其特征在于，包括如下步骤：步骤一：依照荷载试验加载工况建立桥梁有限元计算模型；步骤二：计算结构位移，将位移向量按照测量值和未测量值分开，构建目标函数；步骤三：选择全部单元刚度作为待识别参数；步骤四：采用遗传算法，设置待识别参数的上下限进行初步优化计算；步骤五：将遗传算法计算结果作为初始值，采用L‑M算法进行第二轮优化计算；步骤六：根据测点布置结合灵敏度分析判断刚度识别的有效区域。

【技术特征摘要】
1.一种基于静载试验识别连续梁桥实际刚度的方法，其特征在于，包括如下步骤：步骤一：依照荷载试验加载工况建立桥梁有限元计算模型；步骤二：计算结构位移，将位移向量按照测量值和未测量值分开，构建目标函数；步骤三：选择全部单元刚度作为待识别参数；步骤四：采用遗传算法，设置待识别参数的上下限进行初步优化计算；步骤五：将遗传算法计算结果作为初始值，采用L-M算法进行第二轮优化计算；步骤六：根据测点布置结合灵敏度分析判断刚度识别的有效区域。2.根据权利要求1所述的基于静载试验识别连续梁桥实际刚度的方法，其特征在于，步骤一，桥梁有限元计算模型的力学方程为：[k][δ]＝[F](1)式中[k]为结构的总刚度矩阵，[δ]为结构的位移矩阵，[F]为荷载矩阵，当桥梁静载试验涉及多个加载工况，将不同工况组成荷载矩阵，每列表示不同的工况，每行表示不同的节点号。3.根据权利要求2所述的基于静载试验识别连续梁桥实际刚度的方法，其特征在于，步骤二，计算出结构的位移矩阵，具体如下：[δ]＝[F][K]-1(2)计算出的位移矩阵[δ]为模型所有节点的位移，根据测点布置，从计算的位移矩阵中提取出测点的计算值，当实测值较少时，将[δ]按照测量值和未测量值分开表示，则上式表示为：式中δa为测...

【专利技术属性】
技术研发人员：谢海龙，张俊成，曹蔚枝，郭茂泉，姬兵亮，宋茂林，罗鹏，王彬，何鑫，
申请(专利权)人：山西省交通科学研究院，
类型：发明
国别省市：山西,14