CN112016143A - 一种梁式桥跨中竖向位移随温度变化的计算方法 - Google Patents

一种梁式桥跨中竖向位移随温度变化的计算方法 Download PDF

Info

Publication number
CN112016143A
CN112016143A CN202010814388.4A CN202010814388A CN112016143A CN 112016143 A CN112016143 A CN 112016143A CN 202010814388 A CN202010814388 A CN 202010814388A CN 112016143 A CN112016143 A CN 112016143A
Authority
CN
China
Prior art keywords
span
bridge
vertical displacement
main beam
node
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
CN202010814388.4A
Other languages
English (en)
Other versions
CN112016143B (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.)
University of Science and Technology Beijing USTB
Original Assignee
University of Science and Technology Beijing USTB
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 University of Science and Technology Beijing USTB filed Critical University of Science and Technology Beijing USTB
Priority to CN202010814388.4A priority Critical patent/CN112016143B/zh
Priority to US17/098,475 priority patent/US11579038B2/en
Publication of CN112016143A publication Critical patent/CN112016143A/zh
Application granted granted Critical
Publication of CN112016143B publication Critical patent/CN112016143B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M5/00—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
    • G01M5/0008—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings of bridges
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M5/00—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
    • G01M5/0041—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by determining deflection or stress
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Civil Engineering (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Architecture (AREA)
  • Bridges Or Land Bridges (AREA)

Abstract

本发明提供一种梁式桥跨中竖向位移随温度变化的计算方法，属于桥梁结构分析与监测技术领域。该方法将主梁在各个支座处的结点转角设为未知量，根据结点的弯矩平衡条件，建立方程，然后引入数列，建立各未知量之间的数量关系，将关系式带入方程，解出各个结点的转角解析公式，再通过叠加原理求出各个结点的弯矩解析公式，最后通过虚功原理求出各跨主梁的跨中竖向位移的解析公式。该方法针对边跨长度相等且跨数任意的等截面梁式桥，给出的解析公式属于精确解，计算方便，通用性强，而且从机理层面揭示了梁式桥在温度变化下的行为规律，可用于指导桥梁监测系统的测点布设，并为温度变形基准模型的建立提供先验知识。

Description

一种梁式桥跨中竖向位移随温度变化的计算方法
技术领域
本发明涉及桥梁结构分析与监测技术领域，特别是指一种梁式桥跨中竖向位移随温度变化的计算方法。
背景技术
简支梁和连续梁是交通运输网络中最常使用的桥梁形式，本发明将其统称为梁式桥。梁式桥的结构变形是桥梁运营期结构健康监测关注的重点，通常由主梁梁端伸缩位移和主梁跨中竖向位移表征。梁式桥的主梁伸缩位移已经得到了很好的研究，可通过一维热胀冷缩公式进行估计。相比之下，针对主梁跨中竖向位移随温度变化的研究要少得多。现场监测表明，桥梁顶面因接受阳光的直接照射，在晴朗的天气条件下会与桥梁底面之间产生明显的竖向温度差，从而影响桥梁在竖直方向的变形甚至改变结构的受力状态。为正确理解梁式桥在温度变化下的行为规律，并为桥梁结构健康监测系统的传感器布设提供理论依据，有必要研究主梁顶、底面之间的温度差与主梁跨中竖向位移的关系。
目前计算梁式桥温度变形的方法主要有两类：一类是基于结构力学原理的理论分析，另一类是基于有限元计算软件的数值模拟。第一类方法通常按照三弯矩方程进行计算。三弯矩方程的应用比较麻烦，而且对于不同跨径布置的结构需要分别进行分析，通用性差。数值模拟也是一个费时费力的过程，不易揭示现象背后的机理，且有一定的技术门槛，难以满足工程师现场应用的需要。关于梁式桥在顶面与底面温度差作用下的跨中竖向位移，至今尚未看到通用性强、形式简单、参数关系清晰的统一计算公式。
发明内容
本发明要解决的技术问题是提供一种梁式桥跨中竖向位移随温度变化的计算方法。
该方法适用于各跨主梁为同种材质的等截面直线桥，且主梁顶面与底面之间的温度差沿桥梁走向处处相等。若桥梁总跨数为奇数，则最中间一跨为主跨，其余为边跨；若桥梁总跨数为偶数，则最中间2跨为主跨，其余为边跨。桥梁主跨长度可以与边跨不同，但所有边跨的长度相等。该方法的计算步骤如下：
(1)将主梁在各个支座处的结点转角设为未知量，根据结点的弯矩平衡条件，建立方程：
当梁式桥的总跨数为奇数时：
Figure BDA0002632157060000021
当梁式桥的总跨数为偶数时：
Figure BDA0002632157060000022
其中：zk、zk+1、zk-1是结点转角未知量，以顺时针转动为正，k是变量的下标，取k＝1,2,…,u，u是未知量的个数，取
Figure BDA0002632157060000023
即取不小于n2的整数中最小的一个，n是梁式桥的总跨数；因为结构和温度变化具有对称性，所以独立的结点转角未知量个数为u个。
当总跨数n为奇数时，设位于中央的第u跨为主跨，跨长为l1；其余各跨为边跨，长度均为l0；当总跨数n为偶数时，设位于中央的第u和u+1跨为主跨，跨长均为l1；其余各跨为边跨，长度均为l0；
i0和i1分别是边跨和主跨主梁的抗弯线刚度，定义为i0＝(E·I)/l0和i1＝(E·I)/l1，其中，E为主梁材料的弹性模量，I为主梁截面的抗弯惯性矩；
MT是由主梁顶面与底面温度差引起的两端固定梁的固端弯矩，以主梁下缘受拉为正，MT＝(α·ΔT·E·I)/h，其中，α是主梁材料的线膨胀系数；ΔT为主梁顶面与底面之间的温度差，以顶面温度高于底面温度为正；h为截面高度；
在上面的公式中，当u＝1时，则只有式(A1)、(B1)有效，且此时的梁式桥只有主跨而没有边跨，可令l0＝l1；当u＝2时，则只有式(A1)、(A3)、(B1)、(B3)有效；当u≥3时，式(A1)、(A2)、(A3)、(B1)、(B2)、(B3)均有效。
(2)引入数列，建立各个结点转角未知量之间的数量关系：
引入数列{ak}，其通项为
Figure BDA0002632157060000031
其中，常数
Figure BDA0002632157060000032
常数λ对奇数跨梁式桥取
Figure BDA0002632157060000033
对偶数跨梁式桥取
Figure BDA0002632157060000034
ξ为边跨与主跨的长度比，即ξ＝l0/l1；数列{ak}的下标k＝1,2,…,u-1。
借助数列{ak}，zk与zk-1之间的关系可以表示为：
zk＝-au-k+1·zk-1
其中：下标k＝2,3,…,u；
(3)将步骤(2)中变量间的关系式代入步骤(1)中的方程中，解出各个结点转角的解析公式：
Figure BDA0002632157060000035
其中，k＝1,2,…,u；对于奇数跨梁式桥，λ取λ1，且根据对称性可得，z2u+1-k＝-zk；对于偶数跨梁式桥，λ取λ2，且根据对称性可得，z2u+2-k＝-zk，zu+1＝0；
(4)利用步骤(3)的结果，通过叠加原理求出各个结点的弯矩解析公式：
若以使主梁下缘受拉的弯矩为正弯矩，则当梁式桥的总跨数为奇数时，第k个结点的弯矩Mk为
Figure BDA0002632157060000036
此时，k＝1,2,…,u-1；
第u个结点的弯矩Mu为
Figure BDA0002632157060000037
根据对称性，第2u+1-k个结点的弯矩等于第k个结点，即
M2u+1-k＝Mk，此时，k＝1,2,…,u；
当梁式桥的总跨数为偶数时，第k个结点的弯矩Mk为
Figure BDA0002632157060000041
此时，k＝1,2,…,u；
第u+1个结点的弯矩Mu+1为
Figure BDA0002632157060000042
根据对称性，第2u+2-k个结点的弯矩等于第k个结点，即
M2u+2-k＝Mk，此时，k＝1,2,…,u；
(5)利用步骤(4)的结果，通过虚功原理求出各跨主梁的跨中竖向位移的解析公式：
若规定主梁跨中向下运动对应于正位移，则当梁式桥的总跨数为奇数时，第k跨主梁的跨中竖向位移Δk为
Figure BDA0002632157060000043
此时，k＝1,2,…,u-1；
第u跨主梁的跨中竖向位移Δu为
Figure BDA0002632157060000044
根据对称性，第2u-k跨主梁的跨中竖向位移Δ2u-k为
Δ2u-k＝Δk，此时，k＝1,2,…,u-1；
当梁式桥的总跨数为偶数时，第k跨主梁的跨中竖向位移Δk为
Figure BDA0002632157060000045
此时，k＝1,2,…,u-1；
第u跨主梁的跨中竖向位移Δu为
Figure BDA0002632157060000046
根据对称性，第2u+1-k跨主梁的跨中竖向位移Δ2u+1-k为
Δ2u+1-k＝Δk，此时，k＝1,2,…,u。
上述梁式桥为同种材质的等截面直线桥；主梁顶面与底面之间的温度差沿桥梁走向处处相等；若桥梁总跨数为奇数，则最中间一跨为主跨，其余为边跨；若桥梁总跨数为偶数，则最中间2跨为主跨，其余为边跨；桥梁所有边跨的长度相等。
当梁式桥各跨长度相等时，
奇数跨梁式桥第k跨主梁的跨中竖向位移Δk统一写为
Figure BDA0002632157060000051
此时，k＝1,2,…,u；
根据对称性，第2u-k跨主梁的跨中竖向位移Δ2u-k为
Δ2u-k＝Δk，此时，k＝1,2,…,u-1；
偶数跨梁式桥第k跨主梁的跨中竖向位移Δk统一写为
Figure BDA0002632157060000052
此时，k＝1,2,…,u；
根据对称性，第2u+1-k跨主梁的跨中竖向位移Δ2u+1-k为
Δ2u+1-k＝Δk，此时，k＝1,2,…,u。
当梁式桥各跨长度相等时，
由主梁顶面与底面之间温度差引起的最外侧边跨主梁的跨中竖向位移的变化幅度大于其余各跨，当总跨数n增多时，最外侧边跨主梁的跨中竖向位移趋近于极限值
Figure BDA0002632157060000053
Figure BDA0002632157060000054
由主梁顶面与底面之间温度差引起的中间跨主梁的跨中竖向位移的变化幅度小于其余各跨，当总跨数n增多时，中间跨主梁的跨中竖向位移趋近于0。
本发明的上述技术方案的有益效果如下：
上述方案中，给出了任意跨梁式桥跨中竖向位移随主梁顶面与底面温度差变化的统一计算公式，利用计算器便可方便地得到精确计算结果，通用性强，易于进行参数分析。本方法避免了繁琐的结构力学分析，也避免了有限元数值模拟，非常适用于现场匡算和测试结果检核，具有明显的实际应用价值。更重要的是，基于本方法得到的计算公式，可以得到梁式桥内力和变形随温度变化的普遍性规律，并可发现温度变形随梁式桥跨数增多时的极限值，具有明确的理论意义。另外，本方法所得结果可用于指导梁式桥结构健康监测系统的测点布设，并为温度变形基准模型的建立提供先验知识。
附图说明
图1为本发明实施例中奇数跨梁式桥的分析模型；
图2为本发明实施例中偶数跨梁式桥的分析模型。
具体实施方式
为使本发明要解决的技术问题、技术方案和优点更加清楚，下面将结合附图及具体实施例进行详细描述。
本发明提供一种梁式桥跨中竖向位移随温度变化的计算方法。
该方法包括步骤如下：
(1)将主梁在各个支座处的结点转角设为未知量，根据结点的弯矩平衡条件，建立方程：
当梁式桥的总跨数为奇数时：
Figure BDA0002632157060000061
当梁式桥的总跨数为偶数时：
Figure BDA0002632157060000062
其中：zk、zk+1、zk-1是结点转角未知量，以顺时针转动为正，k是变量的下标，取k＝1,2,…,u，u是未知量的个数，取
Figure BDA0002632157060000063
即取不小于n2的整数中最小的一个，n是梁式桥的总跨数；
当总跨数n为奇数时，设位于中央的第u跨为主跨，跨长为l1；其余各跨为边跨，长度均为l0；当总跨数n为偶数时，设位于中央的第u和u+1跨为主跨，跨长均为l1；其余各跨为边跨，长度均为l0；
i0和i1分别是边跨和主跨主梁的抗弯线刚度，定义为i0＝(E·I)/l0和i1＝(E·I)/l1，其中，E为主梁材料的弹性模量，I为主梁截面的抗弯惯性矩；
MT是由主梁顶面与底面温度差引起的两端固定梁的固端弯矩，以主梁下缘受拉为正，MT＝(α·ΔT·E·I)/h，其中，α是主梁材料的线膨胀系数；ΔT为主梁顶面与底面之间的温度差，以顶面温度高于底面温度为正；h为截面高度；
(2)引入数列，建立各个结点转角未知量之间的数量关系：
引入数列{ak}，其通项为
Figure BDA0002632157060000071
其中，常数
Figure BDA0002632157060000072
常数λ对奇数跨梁式桥取
Figure BDA0002632157060000073
对偶数跨梁式桥取
Figure BDA0002632157060000074
ξ为边跨与主跨的长度比，即ξ＝l0/l1；数列{ak}的下标k＝1,2,…,u-1。
借助数列{ak}，zk与zk-1之间的关系可以表示为：
zk＝-au-k+1·zk-1
其中：下标k＝2,3,…,u；
(3)将步骤(2)中变量间的关系式代入步骤(1)中的方程中，解出各个结点转角的解析公式：
Figure BDA0002632157060000075
其中，k＝1,2,…,u；对于奇数跨梁式桥，λ取λ1，且根据对称性可得，z2u+1-k＝-zk；对于偶数跨梁式桥，λ取λ2，且根据对称性可得，z2u+2-k＝-zk，zu+1＝0；
(4)利用步骤(3)的结果，通过叠加原理求出各个结点的弯矩解析公式：
若以使主梁下缘受拉的弯矩为正弯矩，则当梁式桥的总跨数为奇数时，第k个结点的弯矩Mk为
Figure BDA0002632157060000076
此时，k＝1,2,…,u-1；
第u个结点的弯矩Mu为
Figure BDA0002632157060000077
根据对称性，第2u+1-k个结点的弯矩等于第k个结点，即
M2u+1-k＝Mk，此时，k＝1,2,…,u；
当梁式桥的总跨数为偶数时，第k个结点的弯矩Mk为
Figure BDA0002632157060000081
此时，k＝1,2,…,u；
第u+1个结点的弯矩Mu+1为
Figure BDA0002632157060000082
根据对称性，第2u+2-k个结点的弯矩等于第k个结点，即
M2u+2-k＝Mk，此时，k＝1,2,…,u；
(5)利用步骤(4)的结果，通过虚功原理求出各跨主梁的跨中竖向位移的解析公式：
若规定主梁跨中向下运动对应于正位移，则当梁式桥的总跨数为奇数时，第k跨主梁的跨中竖向位移Δk为
Figure BDA0002632157060000083
此时，k＝1,2,…,u-1；
第u跨主梁的跨中竖向位移Δu为
Figure BDA0002632157060000084
根据对称性，第2u-k跨主梁的跨中竖向位移Δ2u-k为
Δ2u-k＝Δk，此时，k＝1,2,…,u-1；
当梁式桥的总跨数为偶数时，第k跨主梁的跨中竖向位移Δk为
Figure BDA0002632157060000085
此时，k＝1,2,…,u-1；
第u跨主梁的跨中竖向位移Δu为
Figure BDA0002632157060000086
根据对称性，第2u+1-k跨主梁的跨中竖向位移Δ2u+1-k为
Δ2u+1-k＝Δk，此时，k＝1,2,…,u。
在具体应用中，步骤(1)中以主梁在各个支座处的结点转角为未知量，建立结点弯矩平衡方程的推导具体为：
对于附图1中的奇数跨梁式桥分析模型，总跨数为n＝2·u-1，其中整数u≥1，并可记为
Figure BDA0002632157060000091
(即取不小于n2的整数中最小的一个)。设第u跨为主跨，长度为l1，其余各跨为边跨，长度均为l0。记边跨和主跨主梁的抗弯线刚度分别为i0和i1，其中i0＝(E·I)/l0，i1＝(E·I)/l1。其中，E为主梁材料的弹性模量，I为主梁截面的抗弯惯性矩。主梁截面的高度为h，主梁顶面与底面之间的温度差为ΔT，主梁材料的线膨胀系数为α。若引入边跨与主跨的长度比ξ＝l0/l1，则i1＝i0ξ。
根据结构力学的转角位移方程，可以写出主梁在每个结点处的弯矩平衡方程。对于结点1，有：
4i0·z1+2i0·z2+MT＝0 (1)
其中MT为两端固定单跨梁由顶面与底面温度差引起的固端弯矩，若以主梁下缘受拉为正，则查表可得
Figure BDA0002632157060000092
对于结点k＝2,3,…,u-1，同样可以列出弯矩平衡方程：
2i0·zk-1+8i0·zk+2i0·zk+1＝0 (3)
对于结点u，弯矩平衡方程为：
2i0·zu-1+(4i0+4i1)·zu+2i1·zu+1＝0 (4)
由于结构和温度作用存在对称性，所以zk＝-z2u+1-k，其中k＝1,2,…,u，因此式(4)成为：
2i0·zu-1+(4i0+2i1)·zu＝0 (5)
同时，对结点k＝u+1,u+2，…,2u处列写的弯矩平衡方程与式(1)、(3)、(5)等效，因此独立的结点转角未知量个数为u个。
同理，对于附图2中的偶数跨梁式桥分析模型，总跨数为n＝2·u，其中整数u≥1。设位于中央的第u和u+1跨为主跨，跨长均为l1；其余各跨为边跨，长度均为l0。根据结构力学的转角位移方程，结点1处主梁的弯矩平衡方程与式(1)相同，结点k＝2,3,…,u-1处的弯矩平衡方程与式(3)相同，结点u处的弯矩平衡方程为：
2i0·zu-1+(4i0+4i1)·zu+2i1·zu+1＝0 (6)
由于结构和温度作用存在对称性，所以zk＝-z2u+2-k，其中k＝1,2,…,u，且zu+1＝0。因此式(6)成为：
2i0·zu-1+(4i0+4i1)·zu＝0 (7)
同时，对结点k＝u+1处列写的弯矩平衡方程为恒等式，对结点k＝u+2,u+3,…,2u+1处列写的弯矩平衡方程与对结点k＝1,2,…,u处列写的平衡方程等效，因此独立的结点转角未知量个数仍为u个。
在上面的推导中，当u＝1时，则只有式(1)有效，且此时只有主跨而没有边跨，可令l0＝l1；当u＝2时，则只有式(1)、(5)或式(1)、(7)有效；当u≥3时，以上各式均有效。
步骤(2)中引入数列，建立各个结点转角未知量之间的数量关系的推导具体为：
对奇数跨梁式桥的情况，由式(5)可得
Figure BDA0002632157060000101
令a1＝1(2+ξ)，并将zu＝-a1·zu-1代入式(3)，可得zk与zk-1之间的关系(k＝u-1,u-2,…,2)：
zk＝-au-k+1·zk-1 (9)
其中数列ak(k＝2,3,…,u-1)满足：
Figure BDA0002632157060000102
式(10)两边同时减去常数
Figure BDA0002632157060000103
即
Figure BDA0002632157060000104
对式(11)右边做恒等变形：
Figure BDA0002632157060000105
令分子中
Figure BDA0002632157060000106
可得
Figure BDA0002632157060000107
根据式(11)可得到如下两式：
Figure BDA0002632157060000108
Figure BDA0002632157060000111
上面两式相除可得：
Figure BDA0002632157060000112
因此数列
Figure BDA0002632157060000113
是以
Figure BDA0002632157060000114
为公比的等比数列，首项为
Figure BDA0002632157060000115
因此
Figure BDA0002632157060000116
在已知a1＝1/(2+ξ)时，可根据式(16)求出ak的表达式(k＝1,2,…u-1)：
Figure BDA0002632157060000117
式中常数
Figure BDA0002632157060000118
对偶数跨梁式桥的情况，由式(7)可得
Figure BDA0002632157060000119
此时令a1＝1(2+2ξ)，其余推导过程与奇数跨梁式桥相同，则zk与zk-1之间仍然满足式(9)的关系式，其中ak的表达式(k＝1,2,…,u-1)式(17)仍然适用，只需将λ＝λ1修改为
Figure BDA00026321570600001110
步骤(3)将步骤(2)中变量间的关系式代入步骤(1)中的方程中，解出各个结点转角的解析公式。具体为：
根据式(17)中ak的通项可知：
Figure BDA00026321570600001111
将z2＝-au-1·z1代入式(1)可以解出z1：
Figure BDA00026321570600001112
根据式(9)可以进一步得到其余未知量的解：
Figure BDA0002632157060000121
其中，式(21)适用于k＝1,2,…,u。对于奇数跨梁式桥，λ取λ1，且根据对称性可得，z2u+1-k＝-zk。对于偶数跨梁式桥，λ取λ2，且根据对称性可得，z2u+2-k＝-zk，zu+1＝0。
步骤(4)中利用步骤(3)的结果，通过叠加原理求出各个结点的弯矩解析公式。具体为：
规定以使主梁下缘受拉的弯矩为正弯矩。对于奇数跨梁式桥，第k个(k＝1,2,…,u-1)结点的弯矩Mk等于
Mk＝4i0·zk+2i0·zk+1+MT (22)
将式(2)和式(21)代入式(22)，化简可得：
Figure BDA0002632157060000122
由于对称性(zu+1＝-zu)，所以第u个结点的弯矩为
Figure BDA0002632157060000123
根据对称性，第2u+1-k个结点(k＝1,2,…,u)的弯矩为
M2u+1-k＝Mk (25)
对于偶数跨梁式桥，第k个(k＝1,2,…,u)结点的弯矩Mk的表达式的形式与式(22)相同，但与奇数跨梁式桥不同的是，zk的表达式中应用λ2代替λ1。由此可得偶数跨梁式桥第k个(k＝1,2,…,u)结点的弯矩Mk为
Figure BDA0002632157060000124
由于zu+1＝0，所以第u+1个结点的弯矩Mu+1为：
Figure BDA0002632157060000125
根据对称性，第2u+2-k个(k＝1,2,…,u)结点的弯矩为
M2u+2-k＝Mk (28)
步骤(5)中利用步骤(4)的结果，通过虚功原理求出各跨主梁的跨中竖向位移的解析公式。具体为：
规定主梁跨中向下运动对应于正位移。当梁式桥的总跨数为奇数时，第k跨(k＝1,2,…,u-1)主梁的跨中竖向位移可通过虚功原理并结合图乘法求得：
Figure BDA0002632157060000131
同理可得第u跨主梁的跨中竖向位移为
Figure BDA0002632157060000132
根据对称性，第2u-k跨(k＝1,2,…,u-1)主梁的跨中竖向位移为
Δ2u-k＝Δk (31)
当梁式桥的总跨数为偶数时，第k跨(k＝1,2,…,u-1)主梁的跨中竖向位移为：
Figure BDA0002632157060000133
同理可得第u跨主梁的跨中竖向位移为
Figure BDA0002632157060000134
根据对称性，第2u+1-k跨(k＝1,2,…,u)主梁的跨中竖向位移为
Δ2u+1-k＝Δk (34)
当梁式桥各跨长度相等时，有ξ＝1。对奇数跨梁式桥来说，
Figure BDA0002632157060000141
此时用于计算各跨主梁跨中竖向位移的式(29)和式(30)可统一写为：
Figure BDA0002632157060000142
其中，k＝1,2,…,u。根据对称性，第2u-k跨(k＝1,2,…,u-1)主梁的跨中竖向位移仍按式(31)计算。
当梁式桥各跨长度相等时，对偶数跨梁式桥来说，
Figure BDA0002632157060000143
此时用于计算各跨主梁跨中竖向位移的式(32)和式(33)可统一写为：
Figure BDA0002632157060000144
其中，k＝1,2,…,u。根据对称性，第2u+1-k跨(k＝1,2,…,u)主梁的跨中竖向位移仍按式(34)计算。
根据式(35)可知，对于固定u的奇数跨梁式桥来说，各跨主梁的跨中竖向位移的幅度|Δk|正比于
Figure BDA0002632157060000145
由于
Figure BDA0002632157060000146
所以可引入两个函数f1(x)＝x+1/x和
Figure BDA0002632157060000147
则通过复合函数有
Figure BDA0002632157060000148
由于k＝1,2,…,u，所以g1(k)≥1且随k的增加而减小；另一方面，根据导数f1′(x)＝1-1/x2可知，当x≥1时，f1(x)将单调增加。因此，f1(g1(k))将随k的增加而单调减小，即各跨主梁的跨中竖向位移的幅度以第k＝1跨为最大，第k＝u跨为最小。
同理，对于固定u的偶数跨梁式桥来说，由式(36)可知，各跨主梁的跨中竖向位移的幅度|Δk|正比于
Figure BDA0002632157060000149
引入函数
Figure BDA00026321570600001410
则通过复合函数可得
Figure BDA00026321570600001411
因k＝1,2,…,u，g1(k)≥1并随k的增加而减小；同时，导数
Figure BDA00026321570600001412
所以f2(g1(k))仍然随k的增加而单调减小。因此，对偶数跨梁式桥来说，各跨主梁的跨中竖向位移的幅度以第k＝1跨为最大，第k＝u和k＝u+1跨为最小。
综合以上可知，无论是奇数跨或偶数跨梁式桥，由主梁顶面与底面之间的温度差引起的各跨主梁跨中竖向位移的变化幅度以最外侧边跨为最大，以位于中央的主跨为最小。
对第k＝1跨主梁，当u趋向于无穷大时，式(35)和式(36)均趋向于常数
Figure BDA0002632157060000151
Figure BDA0002632157060000152
因此，梁式桥最外侧边跨主梁的跨中竖向位移随桥梁跨数的增加存在极限值
Figure BDA0002632157060000153
同理，对第k＝u跨主梁，当u趋向于无穷大时，式(35)和式(36)均趋向于0。因此，梁式桥位于中央的主跨主梁的跨中竖向位移随桥梁跨数的增加也存在极限值，该极限值为0。
下面结合具体实施例予以说明。
实施例1
某梁式桥为6跨连续梁结构，每跨跨径为110m，主梁材料为钢材，梁高4.5m。由于阳光直射，主梁顶面与底面的温度差为20℃。需确定第1跨和第3跨跨中桥面的竖向位移。
由上述条件可知，跨数n＝6，为偶数跨连续梁的情况，独立的结点转角未知量个数u＝3。主跨与边跨的跨径相等，l0＝l1＝110m，主梁线膨胀系数为α＝1.2×10-5℃-1，梁高h＝4.5m，主梁顶面与底面的温度差为ΔT＝20℃。将上述参数代入式(36)可得，当k＝1时，Δ1＝-0.0295m；当k＝3时，Δ3＝-0.0016m。可见第1跨主梁跨中桥面的竖向位移幅度大于第3跨，位移数值前的负号表示跨中桥面向上运动。
实施例2
设梁式桥的跨数为n，每跨跨径均为110m，主梁材料为钢材，梁高4.5m，主梁顶面与底面的温度差为20℃。计算n＝1,2,…,10时，梁式桥最外侧边跨主梁的跨中竖向位移。
由上述条件可知以下参数：l0＝l1＝110m，α＝1.2×10-5℃-1，h＝4.5m，ΔT＝20℃，k＝1。根据n是奇数或偶数，分别将这些参数代入式(35)或式(36)，并令
Figure BDA0002632157060000154
(即取不小于n2的整数中最小的一个)，从而得到不同总跨数n的梁式桥最外侧边跨主梁的跨中竖向位移(以跨数增加为序)：-0.0807m，-0.0202m，-0.0323m，-0.0288m，-0.0297m，-0.0295m，-0.0295m，-0.0295m，-0.0295m，-0.0295m。不难发现从n＝6开始，最外侧边跨主梁的跨中竖向位移就很接近-0.0295m了，该值与式(37)计算出的极限值一致。这一算例清楚地表明，当连续梁跨数较多时，由主梁顶面与底面的温度差引起的主梁跨中竖向位移较小，在现有仪器设备的测试精度下不容易测准。
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明所述原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (4)

1.一种梁式桥跨中竖向位移随温度变化的计算方法，其特征在于：包括步骤如下：
(1)将主梁在各个支座处的结点转角设为未知量，根据结点的弯矩平衡条件，建立方程：
当梁式桥的总跨数为奇数时：
Figure FDA0002632157050000011
当梁式桥的总跨数为偶数时：
Figure FDA0002632157050000012
其中：zk、zk+1、zk-1是结点转角未知量，以顺时针转动为正，k是变量的下标，取k＝1,2,…,u，u是未知量的个数，取
Figure FDA0002632157050000013
即取不小于n/2的整数中最小的一个，n是梁式桥的总跨数；
当总跨数n为奇数时，设位于中央的第u跨为主跨，跨长为l1；其余各跨为边跨，长度均为l0；当总跨数n为偶数时，设位于中央的第u和u+1跨为主跨，跨长均为l1；其余各跨为边跨，长度均为l0；
i0和i1分别是边跨和主跨主梁的抗弯线刚度，定义为i0＝(E·I)/l0和i1＝(E·I)/l1，其中，E为主梁材料的弹性模量，I为主梁截面的抗弯惯性矩；
MT是由主梁顶面与底面温度差引起的两端固定梁的固端弯矩，以主梁下缘受拉为正，MT＝(α·ΔT·E·I)/h，其中，α是主梁材料的线膨胀系数；ΔT为主梁顶面与底面之间的温度差，以顶面温度高于底面温度为正；h为截面高度；
(2)引入数列，建立各个结点转角未知量之间的数量关系：
引入数列{ak}，其通项为
Figure FDA0002632157050000014
其中，常数
Figure FDA0002632157050000021
常数λ对奇数跨梁式桥取
Figure FDA0002632157050000022
对偶数跨梁式桥取
Figure FDA0002632157050000023
ξ为边跨与主跨的长度比，即ξ＝l0/l1；数列{ak}的下标k＝1,2,…,u-1；
借助数列{ak}，zk与zk-1之间的关系定义为：
zk＝-au-k+1·zk-1
其中：下标k＝2,3,…,u；
(3)将步骤(2)中变量间的关系式代入步骤(1)中的方程中，解出各个结点转角的解析公式：
Figure FDA0002632157050000024
其中，k＝1,2,…,u；对于奇数跨梁式桥，λ取λ1，且根据对称性可得，z2u+1-k＝-zk；对于偶数跨梁式桥，λ取λ2，且根据对称性可得，z2u+2-k＝-zk，zu+1＝0；
(4)利用步骤(3)的结果，通过叠加原理求出各个结点的弯矩解析公式：
若以使主梁下缘受拉的弯矩为正弯矩，则当梁式桥的总跨数为奇数时，第k个结点的弯矩Mk为
Figure FDA0002632157050000025
此时，k＝1,2,…,u-1；
第u个结点的弯矩Mu为
Figure FDA0002632157050000026
根据对称性，第2u+1-k个结点的弯矩等于第k个结点，即M2u+1-k＝Mk，此时，k＝1,2,…,u；
当梁式桥的总跨数为偶数时，第k个结点的弯矩Mk为
Figure FDA0002632157050000027
此时，k＝1,2,…,u；
第u+1个结点的弯矩Mu+1为
Figure FDA0002632157050000028
根据对称性，第2u+2-k个结点的弯矩等于第k个结点，即
M2u+2-k＝Mk，此时，k＝1,2,…,u；
(5)利用步骤(4)的结果，通过虚功原理求出各跨主梁的跨中竖向位移的解析公式：
若规定主梁跨中向下运动对应于正位移，则当梁式桥的总跨数为奇数时，第k跨主梁的跨中竖向位移Δk为
Figure FDA0002632157050000031
此时，k＝1,2,…,u-1；
第u跨主梁的跨中竖向位移Δu为
Figure FDA0002632157050000032
根据对称性，第2u-k跨主梁的跨中竖向位移Δ2u-k为
Δ2u-k＝Δk，此时，k＝1,2,…,u-1；
当梁式桥的总跨数为偶数时，第k跨主梁的跨中竖向位移Δk为
Figure FDA0002632157050000033
此时，k＝1,2,…,u-1；
第u跨主梁的跨中竖向位移Δu为
Figure FDA0002632157050000034
根据对称性，第2u+1-k跨主梁的跨中竖向位移Δ2u+1-k为
Δ2u+1-k＝Δk，此时，k＝1,2,…,u。
2.根据权利要求1所述的梁式桥跨中竖向位移随温度变化的计算方法，其特征在于：所述梁式桥为同种材质的等截面直线桥；主梁顶面与底面之间的温度差沿桥梁走向处处相等；若桥梁总跨数为奇数，则最中间一跨为主跨，其余为边跨；若桥梁总跨数为偶数，则最中间2跨为主跨，其余为边跨；桥梁所有边跨的长度相等。
3.根据权利要求1所述的梁式桥跨中竖向位移随温度变化的计算方法，其特征在于：当梁式桥各跨长度相等时，
奇数跨梁式桥第k跨主梁的跨中竖向位移Δk统一写为
Figure FDA0002632157050000041
此时，k＝1,2,…,u；
根据对称性，第2u-k跨主梁的跨中竖向位移Δ2u-k为
Δ2u-k＝Δk，此时，k＝1,2,…,u-1；
偶数跨梁式桥第k跨主梁的跨中竖向位移Δk统一写为
Figure FDA0002632157050000042
此时，k＝1,2,…,u；
根据对称性，第2u+1-k跨主梁的跨中竖向位移Δ2u+1-k为
Δ2u+1-k＝Δk，此时，k＝1,2,…,u。
4.根据权利要求1所述的梁式桥跨中竖向位移随温度变化的计算方法，其特征在于：当梁式桥各跨长度相等时，
由主梁顶面与底面之间温度差引起的最外侧边跨主梁的跨中竖向位移的变化幅度大于其余各跨，当总跨数n增多时，最外侧边跨主梁的跨中竖向位移趋近于极限值
Figure FDA0002632157050000043
Figure FDA0002632157050000044
由主梁顶面与底面之间温度差引起的中间跨主梁的跨中竖向位移的变化幅度小于其余各跨，当总跨数n增多时，中间跨主梁的跨中竖向位移趋近于0。
CN202010814388.4A 2020-08-13 2020-08-13 一种梁式桥跨中竖向位移随温度变化的计算方法 Active CN112016143B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202010814388.4A CN112016143B (zh) 2020-08-13 2020-08-13 一种梁式桥跨中竖向位移随温度变化的计算方法
US17/098,475 US11579038B2 (en) 2020-08-13 2020-11-16 Method for calculating temperature-dependent mid-span vertical displacement of girder bridge

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010814388.4A CN112016143B (zh) 2020-08-13 2020-08-13 一种梁式桥跨中竖向位移随温度变化的计算方法

Publications (2)

Publication Number Publication Date
CN112016143A true CN112016143A (zh) 2020-12-01
CN112016143B CN112016143B (zh) 2021-09-07

Family

ID=73504324

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010814388.4A Active CN112016143B (zh) 2020-08-13 2020-08-13 一种梁式桥跨中竖向位移随温度变化的计算方法

Country Status (2)

Country Link
US (1) US11579038B2 (zh)
CN (1) CN112016143B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023087399A1 (zh) * 2021-11-19 2023-05-25 中车长春轨道客车股份有限公司 一种无源自运动结构设计方法及低温温控开关

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117592252B (zh) * 2023-10-26 2024-09-17 桐乡市苏台高速公路投资开发有限公司 π型钢-混连续刚构桥准静力位移影响线快速计算方法
CN117951771B (zh) * 2023-11-29 2024-07-16 中铁五局集团电务工程有限责任公司 面向最优受力的连续钢箱梁桥施工支架位置快速确定方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150286751A1 (en) * 2014-04-08 2015-10-08 Gci, Inc. Load rating of bridges, including substructure
CN109101745A (zh) * 2018-08-28 2018-12-28 北京科技大学 一种斜拉桥跨中挠度随温度变化的估算方法
CN109933746A (zh) * 2019-03-26 2019-06-25 北京科技大学 悬索桥中跨主缆跨中挠度和高程随温度变化的估算方法
CN110941872A (zh) * 2019-12-16 2020-03-31 北京科技大学 悬索桥主缆垂度和塔顶水平位移随温度变化的确定方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150286751A1 (en) * 2014-04-08 2015-10-08 Gci, Inc. Load rating of bridges, including substructure
CN109101745A (zh) * 2018-08-28 2018-12-28 北京科技大学 一种斜拉桥跨中挠度随温度变化的估算方法
CN109933746A (zh) * 2019-03-26 2019-06-25 北京科技大学 悬索桥中跨主缆跨中挠度和高程随温度变化的估算方法
CN110941872A (zh) * 2019-12-16 2020-03-31 北京科技大学 悬索桥主缆垂度和塔顶水平位移随温度变化的确定方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
周毅等: "斜拉桥跨中竖向位移的温度灵敏度系数研究", 《工程力学》 *
周毅等: "温度对斜拉桥跨中竖向位移的作用机理研究", 《工程力学》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023087399A1 (zh) * 2021-11-19 2023-05-25 中车长春轨道客车股份有限公司 一种无源自运动结构设计方法及低温温控开关

Also Published As

Publication number Publication date
US20220050008A1 (en) 2022-02-17
CN112016143B (zh) 2021-09-07
US11579038B2 (en) 2023-02-14

Similar Documents

Publication Publication Date Title
CN112016143B (zh) 一种梁式桥跨中竖向位移随温度变化的计算方法
Karina Calculation of deformations of a cantilever-frame planar truss model with an arbitrary number of panels
CN106896156B (zh) 逐跨均匀荷载面曲率差的梁结构损伤识别方法
Han et al. Exact dynamic characteristic analysis of a double-beam system interconnected by a viscoelastic layer
CN108121865B (zh) 一种基于伴随变量法的含有多不确定参数的结构区间响应传播分析方法
Chen et al. Comparison of regularization methods for moving force identification with ill-posed problems
CN104951603B (zh) 一种基于悬链线的变电站软导线下料长度数值求解方法
CN109101745B (zh) 一种斜拉桥跨中挠度随温度变化的估算方法
CN111783201B (zh) 一种三跨自锚式悬索桥动力特性的快速分析方法
Cong et al. Planar multimodal 1: 2: 2 internal resonance analysis of cable-stayed bridge
Ruocco et al. Hencky bar-chain model for buckling analysis of non-uniform columns
CN111308979B (zh) 基于多率迟延状态空间模型的u控制系统
CN106096257A (zh) 一种非线性索单元分析方法及系统
CN111324949A (zh) 一种考虑噪声影响的工程结构柔度识别方法
CN111783199A (zh) 一种多段式索缆结构动力特性的精细化快速求解方法
CN115828759A (zh) 一种多约束拱桥施工期最优扣索索力快速确定方法和装置
CN108763667B (zh) 大曲率曲线钢-混凝土组合箱梁桥简化设计方法
Jamshidi-Ghaleh et al. Diffusion coefficient measurements of transparent liquid solutions using Moiré deflectometry
CN112287441B (zh) 一种等跨等截面连续梁影响线的计算方法
He et al. A multi-scale wavelet finite element model for damage detection of beams under a moving load
Liu et al. Determination of asymptotic field coefficients for V-notches by the weak form quadrature element method
Ghosh et al. Tusas: A fully implicit parallel approach for coupled phase-field equations
Corradi et al. Inelastic analysis of reinforced concrete frames
Życzkowski Post-buckling analysis of non-prismatic columns under general behaviour of loading
CN111859655B (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