2 桥梁构件的数学与力学基础
一辆车开上桥面,荷载经桥面板传给主梁,再经支座、墩台传向地基。桥梁由相互连接的构件组成;取出其中一片梁、一段拱或一根索,就能分析支座反力、截面应力、构件变形与稳定,以及线形与索力的关系。
本章按“计算模型—数学基础—梁板—受压构件—拱与索”展开。数学基础集中在一节;平衡、微分方程、特征值和优化随不同构件再次出现,每次都对应一个看得见的对象。计算工具与AI可辅助重复运算、绘图和数据整理,本章着重呈现决定结果的对象、假定、平衡与协调、材料关系及边界条件,并将物理量的含义和核验依据放在相应图式与算例旁。后续各桥型章节再把这些构件组织成完整体系。
阅读层次。主线是取出自由体、列平衡,再沿荷载—内力—应力—挠度读梁,继而比较偏心受压、稳定及拱索平衡。数学基础供对应概念的复习;随机过程和调索优化属于向后续章节延伸的内容,可按链接回看。阅读时先确定图中的物体、输入与输出,再对照公式中的每个量。
图中上部辨认梁板、受压构件、拱与索,下部连接共同的分析过程。梁板主要体现弯曲与分配,受压构件体现偏心和稳定,拱与索则对照曲线中的受压与受拉。每个模型都回到同一件事:作用怎样经连接传递,构件怎样响应。
完成本章后,学生能够:
- 从整桥中辨认主梁、桥面板、墩柱、拱肋与拉索,取出一个构件,画出荷载和支承作用,并用合力与合矩核对自由体图。
- 从桥面位移、荷载时程和振动图中读出位移、速度、加速度、频率与振型,说明频率及阻尼怎样改变响应,并区分支承约束、初始状态和作用的随机变化。
- 沿同一片梁指出荷载、剪力、弯矩、截面应力和位移所在的位置,计算简单梁的响应,选择梁、索、板壳或实体单元,并说明截面布置、预应力与施工顺序的影响。
- 从受压墩柱或杆件的偏心与侧向弯曲,解释截面两侧的应力差别和二阶效应,区分材料强度问题与构件失稳问题。
- 对照同跨径的拱与索,说明水平力、曲线形状及受压受拉的关系,计算理想模型的内力,从索形或频率理解索力识别,并把斜拉桥调索要求写成有目标与约束的优化问题。
2.1 从实桥到计算模型
2.1.1 研究对象先于计算公式
研究目标决定模型保留的内容。如果要算一片主梁的纵向弯矩,可以把它简化为具有截面刚度的轴线;如果要算车轮下桥面板的局部弯曲,就必须保留板的宽度、厚度和支承;如果要研究索锚处的应力,锚固范围和连接构造也不能省去。细节多少本身不能判断模型好坏,模型必须能表达所研究的机制。
下面从同一座三维教学桥中选取同一片主梁,依次观察“桥中的构件—取出的构件—作用替代连接”。三种状态沿用相同的几何位置和观察方向,因而可以直接追踪梁端支座、作用点和坐标。教学模型由既有T梁几何派生,选取一跨按独立简支跨处理;画面中的尺寸与计算跨径以支座中心位置为准,不由外观推定实桥参数。
单独打开:从整桥到一片主梁。三种显示状态分别呈现主梁在桥中的位置、取出的主梁,以及用外部作用替代连接后的自由体。梁端接触位置与反力位置对应,向下箭头表示荷载。状态切换与旋转保持同一根梁及其观察方向。
把构件从周围结构中分离,切断连接,并在切断处补上外部对它的作用,得到自由体图。对取出的主梁,支座反力是外力;若研究对象扩大为“梁、支座和桥墩”,梁与支座之间的接触力就成为对象内部的力,基础边界上的作用才需要补画。外力和内力的区别取决于所选对象。
主梁跨中被假想截开时,两部分通过截面传递轴力、剪力、弯矩和扭矩。保留左半梁,截面作用由右半梁施加;保留右半梁,作用方向相反。这里的“截开”是分析手段,真实构件仍保持连续。做图时既要看切断了哪一处,也要看留下哪一侧。
2.1.2 连接和支承怎样进入模型
桥梁构造通过允许或限制运动来传递作用。平面梁节点可用顺梁位移、竖向位移和转角三个自由度描述;空间梁节点还需横向位移及另外两个转角。约束一个自由度,便可能出现与它对应的反力或反力矩。
| 理想化连接或支承 | 允许的运动与传递的作用 | 观察实体时要辨认什么 |
|---|---|---|
| 平面铰支承 | 允许转动,限制两个平动方向;可传两个力分量 | 转动在哪里发生,平移由什么构造限制 |
| 沿梁活动支承 | 允许顺梁移动和转动,限制竖向位移 | 滑移方向、接触面及限位位置 |
| 固定端或刚接 | 限制相应平动和转动;可以传力和力矩 | 构件与节点能否提供所假定的约束 |
| 弹性支承 | 反力随相对位移变化 | 支座、墩和基础的柔度是否需要共同考虑 |
| 仅承压接触 | 受压时接触,分离后不能继续传拉力 | 普通支座是否可能脱空,是否设有抗拔构造 |
“活动”不是所有方向都自由,“铰接”也不是不传水平力。简支梁的一端仍需有沿梁方向的约束以防止整体平移;仅受竖向荷载时,这一约束的水平反力可以为零。支座摩阻、间隙和有限刚度,在温度、制动或地震问题中可能成为主要因素,此时需要修改最初的理想支承。
连接还决定水平力在哪里闭合。拱脚可能把推力交给桥台,也可能由系杆将其闭合在上部结构中;索的拉力通过锚固、桥塔、主梁或锚碇继续传递。辨认完桥型,应当沿每个力分量追到能承受它的构件和基础。梁、拱、索的差别由此转化为可以计算的连接与受力问题。
2.1.3 简化要围绕问题作出
梁理论并不把实体梁变成没有厚度的细线;它用轴线位移与截面量代表一类主要运动。数学上舍去的自由度,必须对应物理上可以暂时忽略的机制。
| 本章常用假定 | 采用的近似与暂不考虑的内容 | 什么情况需要重新考虑 |
|---|---|---|
| 静力 | 加速度引起的惯性力 | 车辆冲击、地震、明显振动 |
| 小变形 | 位移与转角较小,可对几何关系作相应展开;一阶分析另在变形前几何上列平衡 | 轴压显著时,即使侧移较小也需检查二阶效应;大转角或显著索形变化需进一步更新几何 |
| 线弹性 | 应力与应变成正比,暂不计塑性及损伤;相同线性边界和状态下可叠加 | 开裂、屈服、接触变化 |
| 平截面 | 截面变形后仍保持平面,以少量截面运动表示纵向应变;暂不计相应翘曲 | 深梁、锚固区、突变截面和局部承压 |
| 理想柔索 | 忽略抗弯刚度,由张力维持形状 | 短索弯曲刚度、松弛、减振器和复杂锚端 |
静定结构的反力和内力可由平衡求出;超静定结构还要满足变形协调和材料关系。两根平行受拉杆若共同连接于保持水平的刚性横梁,伸长相同,则有
由此得 。刚度比为3∶1时,荷载分担也是3∶1;仅有总力平衡不能得出这一分配。主梁之间、墩梁之间以及多根索之间的共同工作,都需要相应的协调条件。
因此,建立计算关系要同时考虑平衡、几何协调、材料关系和边界条件。时间问题再补初始条件,施工问题再补已有状态。模型能够求解,只说明方程在当前设置下得到了解,并不说明这些设置已经对应实际对象。
2.1.4 用平衡检查第一张计算图
平面静力问题的整体平衡为
力的单位为N或kN,力矩的单位为N·m或kN·m。三维问题需分别检查三个方向的合力和合矩。两个向上的反力之和正确,而荷载位置填错,仍可能不满足力矩平衡。
例2-1 移动荷载,反力为何改变? 另设一片计算跨径 的简支梁,集中力 作用于距左支点 处,忽略自重。对左支点取矩,并列竖向平衡:
所以 、。荷载移至跨中,两端才各承担90 kN。若原位置也平均分配,对左支点的合矩为 ,并未平衡。
约束是否足以防止刚体运动要另行判断。未知反力数恰好等于平衡方程数,仍可能因约束方向平行或作用线布置不当形成可动体系。力能算出来、构件有材料强度、体系不发生刚体运动,是三个需要分别确认的问题。
2.1.5 荷载、内力、应力与位移
本章总体坐标取 顺桥、 横桥、 向上。竖向梁挠度另用 表示,取向下为正;正弯矩使梁下缘受拉。每个局部模型均应显示相应坐标,不能凭图的上下或软件颜色猜正负号。
图中外力作用于梁和边界,内力属于截面,应力属于截面内的材料点;原轴线与变形线之间的距离表示竖向位移。沿梁的静态挠曲线记为 ,同一测点随时间变化的位移在下一节记为 。挠度是竖向位移的常用名称。
| 量 | 它回答的问题 | 读数所在的位置与单位 |
|---|---|---|
| 荷载 | 每米梁受到多大的分布力 | 沿梁位置,kN/m |
| 反力 | 周围支承对所选构件施加多大的力 | 支座或边界,kN |
| 内力 | 截面两侧通过切口传递什么 | 某一截面;kN或kN·m |
| 应力 | 截面上一个小面积承受多大的力密度 | 材料点,MPa |
| 位移 、转角 | 一个点移动多少,截面旋转多少 | 点或截面;mm、rad |
| 刚度 | 材料和截面怎样抵抗相应变形 | 构件;不能与应力或抗力混称 |
| 质量 | 同一加速度对应多大的惯性项 | 集中质量kg或t;应与刚度采用协调单位 |
同一式中必须采用协调的单位。例如 ,而 ;长度的换算因次会直接进入截面刚度。自重是力,质量通过重力加速度换成重量,不能把t直接当作kN。
2.2 桥梁分析的数学基础
桥面受到作用会移动,移动的快慢会变化,不同部位还会一起运动。微积分用来描述这些变化,微分方程把作用与运动联系起来,特征值和特征向量给出结构自身的振动频率与形态。概率描述作用的不确定性,优化则从允许的方案中选出符合目标的状态。下面把这些数学内容放回桥梁中的具体对象。
2.2.1 位移、速度与加速度
图中青色测点沿梁的位置用 表示;它相对原位置的竖向位置变化是竖向位移。把该测点固定下来,其竖向位移随时间的记录写成 ,取向下为正。位移变化的快慢是速度,速度变化的快慢是加速度。下面在同一片梁上分别读取这三个量。
直接拖动梁上的青色测点,就能沿梁读取不同位置的位移;拖动时程图上的游标,便能回到某一时刻。上方的荷载 滑块控制荷载大小,箭头与挠度一起变化;位置和时间也保留滑块,便于精确调整。静态竖向位移沿梁的分布称为挠曲线,本书构件计算记作 ;在某个固定位置随时间变化时用 。挠度就是竖向位移, 则是位置坐标,二者不是同一个量。图中位移有显示放大,读数保持实际数值。
图后的关系写成
三个量的单位依次为m、m/s、m/s²。字母上一个点、两个点分别表示对时间求一次、两次导数。加速度决定惯性项;它不是“位移比较大”的另一种说法。简谐运动中,位移到达端点时速度为零,加速度却指向平衡位置。
微积分同样描述荷载与内力。图中均布荷载每米为 ,长度 上的累计荷载就是 ,对应荷载图下的面积。一般分布荷载的合力为 ;这就是积分在受力图中的直接含义。静力梁中 ,所以剪力为零的位置可以对应弯矩极值。梁板一节再把这些关系连成完整算例。
2.2.2 一个位移坐标的运动方程
一片桥面按某一种主要形态上下弯曲时,可以用一个测点的位移 表示这一形态的幅度。桥的质量、抗弯与支承刚度、振动耗能,分别进入等效质量 、等效刚度 和等效阻尼 ;它们都必须按同一个位移坐标、同一种形态换算。图中始终是桥梁本身。
这个模型的核心关系只有一行:
随加速度变化, 随速度变化, 随位移变化。三项的单位均为力。桥梁的弹性恢复作用与该坐标的位移方向相反,黏性阻尼近似下的阻尼作用与速度方向相反;移至方程左边后写成正号。位移从自重作用后的静平衡位置量起,自重已经由该状态承担。
撤去附加外力后,结构仍可振动。无阻尼时,它自身的固有频率为
刚度增大,振动通常更快;质量增大,振动通常更慢。 的单位是Hz,即每秒振动次数;圆频率 的单位是rad/s。阻尼比 表示相对于临界阻尼的大小。频率和阻尼共同决定动力响应。13
2.2.3 频率与响应的联系
周期荷载的大小和重复快慢是两个独立条件。荷载幅值相同,接近某阶固有频率时,响应仍可能明显放大;能否激发这一阶,还与荷载位置和分布有关。这里先把“输入节奏”和“结构自身节奏”分开,完整的同幅荷载动画、幅频曲线、阻尼比较和算例集中在第9章共振实验。
本章用运动方程理解导数的物理含义;第9章再回答振幅有多大、怎样测量和怎样控制。频率比 无量纲,不是时刻或风速。
2.2.4 支承约束与运动起点
梁端装什么支座,决定它能否移动、能否转动;这就是边界条件。振动从哪个位置、以多大速度开始,则是初始条件。两者在图上分别对应“构件的端部”和“运动的起点”。
对竖向弯曲问题,简支端限制竖向位移而允许转动,因而 、;固定端限制位移与转角,因而 、。没有端部外力与力偶的自由端满足 、,它的位移与转角并不必为零。端部施加力或力矩时,相应内力条件随之改变。
对同一座桥的同一振动形态,初始条件写作 、。桥面具有初始挠曲形态后静止释放,是 ;在平衡位置给一个初速度,是 。它们可以有相同的固有频率,却从不同状态开始运动。实际启动过程通常包含与初始条件有关的过渡响应;第9章共振图显示的是过渡衰减后的稳态部分。
支承、初始位移和初始速度确定以后,运动方程才对应一个具体过程。“从静止开始”只说明初速度为零,不能据此把已有变形也置零。预应力、合龙和拆架同样需要保留前一步已经形成的状态。
2.2.5 特征向量是形态,特征值决定频率
沿桥取三个点,如果它们始终按 的比例往复,这列数描述一种相对形态,即振型向量。中点当前下移10 mm时,两侧点各下移约7.07 mm;比例本身不带毫米。把比例整体乘2,同时把代表幅度除以2,实际运动不变。
有限元模型把节点位移与转角排成向量。忽略外力和阻尼,寻找能够保持比例自由运动的形态,得到
这给线性代数中的特征问题一个具体用途:特征向量对应空间比例,特征值对应圆频率的平方。 刚度矩阵单独的特征值不能直接当作桥梁频率,还需质量、边界和单位一致。14
均匀简支Euler—Bernoulli梁有 、;若 Hz,前四阶为1、4、9、16 Hz。这是特定模型的数学例子,不能通用于斜拉桥和悬索桥。第9章振型与测点直接展示这些形态、节点和时间变化,再解释多模态叠加与实测识别。
特征向量与特征值的几何解释可作数学补充。压力改变构件稳定性时也会出现特征问题,但求得的是临界荷载与屈曲形态;后面的压杆一节再作对照。
2.2.6 连续构件与有限个坐标的联系
一个坐标描述一个主要运动;多个坐标描述多个部位共同运动;整根梁连续的位移则写成 。当位置和时间都需要保留时,连续模型会用到偏微分方程。梁与索的连续模型分别由它们的平衡关系得到。将连续结构离散成节点,或用若干振型表示其运动,就能回到质量、阻尼、刚度与位移组成的矩阵方程。这是本节与第9章动力学、第10章数字建模的衔接。2
2.2.7 概率与随机过程:分布与时程的区别
桥侧风速仪即使一直放在同一位置,读数也会变化。先把这些读数按大小排开,就能看见哪些值常见、哪些值少见;再把它们按时间连起来,才看见风的变化节奏。分布描述取值的机会,过程还保留先后次序。
2.2.7.1 从观测次数到概率密度
下图采用一个明确设定的教学例子:风速的中心取15 m/s,标准差取3 m/s,用正态曲线示意一次观测的分布。这两个数是为了演示而选的参数,并非本桥实测结果,也不是规范规定。直方图把样本按风速区间计数;曲线下的一块面积,表示风速落在这个区间的概率。
拖动阈值,右侧阴影面积随之变化。阈值越靠右,超过它的机会越小。概率密度是曲线的高度,概率是曲线下的面积;不能把某处的高度直接读成百分比。直方图按“次数÷总样本数÷组距”换算成密度后,才能与密度曲线使用同一纵轴。图中的有限样本起伏也不等于理论曲线本身。19
图后的关系是
是一次观测的风速, 是密度, 是累计概率。若把风速减去中心15 m/s,再除以标准差3 m/s,得到无量纲的“离中心几个标准差”;反过来写成 ,只是在还原量纲。例如 对应18 m/s。这不是风速的动力方程,也没有说明相邻时刻怎样关联。 实际采用哪种分布、中心和离散程度多大,需由观测资料及适用条件确定;相同均值和标准差也可能有不同的尾部。
一次取值的分布不保留时间次序;风振还要知道一阵风持续多久、快慢波动怎样组合。完整记录、协方差与弱平稳的对照见下面的选读;年超越概率与事件年率的区别见概率尺度专题,风的时程怎样产生响应见第9章。
选读:共同发生、时间相关与超越次数
2.2.7.2 两种较大作用,同时出现的机会
把同一时间窗口的风速与车流荷载配成一对,每个点就代表一种共同出现的状态。下图使用合成的百分位点示意,并非风与交通的实测联合分布。图中分别把两者最大的10%标为“较大”;右上角同时满足两个条件。这里的10%是便于读图的示例阈值,不能当作工程设计极值。
在独立示例中,两个10%的条件一起满足,概率为 。切换同步程度后,两边各自的10%保持不变,右上角的机会却改变;完全同步的极限示例达到10%。因此,各项都取很大值的组合可能很少见,但不能不看相关性就把概率相乘。真实风与交通是否相关,需要同一时间窗口的数据来判断。
两个事件的共同概率一般写成 ;只有独立时,才简化成 。这也不能代替规范中的作用组合:组合还要考虑作用引起的结构效应、持续时间和设计要求。
2.2.7.3 同样的取值范围,可以有不同的变化节奏
回到风速仪,把读数按先后次序接起来,得到一条时程。同一条件下可能出现许多不同的记录,这些可能记录合起来描述一个随机过程。下面沿青线看一次记录;固定一个时刻,再看其他记录在这个时刻的取值。记录编号不是时间,也不是沿桥的位置。
默认60 s是一分钟,20 Hz表示每0.05 s记录一点,含首尾共1201点;600 s是十分钟,共12001点。记录多长,与每秒记录多少点,是两个量。 降到0.5 Hz时,每2 s才留一点,灰色参考线显示漏掉的快速变化。实际降采样还需先滤除高频成分以避免混叠;密集采样点也不等于彼此独立的样本。
本图用不同随机相位生成24条教学记录,设定平均风速15 m/s、过程标准差3 m/s,按Kaimal型谱分配快慢不同的波动,保留0.002~4 Hz分量后归一化方差;不是桥位实测风,也不是年最大风速模型。15 图上读到的是m/s,计算风力还需空气密度、气动力系数、受风面积等,不能把风速数值直接当成kN。
2.2.7.4 隔一段时间,两次偏离怎样配对
两个过程即使中心相同、波动幅度相同,也可能一个缓慢起伏,一个快速交替。取两次读数,各自减去平均值:都高于平均值或都低于平均值,两个偏差的乘积为正;一高一低,乘积为负。把这种乘积平均,就得到协方差。下图另用两条合成的周期记录专门展示配对,因而相关性可以随节奏再次升高;它不是上图Kaimal型风速记录的统计结果。
图中时间间隔变大,配对关系随之改变;右侧同号与异号区域帮助辨认正、负贡献。协方差的单位是风速单位的平方,即 。再除以两个时刻风速的标准差乘积,得到无量纲的相关系数;对同一个过程在不同时刻配对,称为自相关。零相关只表示这种线性关联为零,不能一般地推断两次读数独立。6
图后的定义为
表示对可能结果取平均, 是相应时刻的平均风速。图中由一段有限记录算出的平均值是样本估计,需要与模型的理论关系区分。分布告诉我们风可能多大,时间相关告诉我们一阵风会怎样持续或交替,二者一起影响桥梁响应。
有时可以在一个合适时间窗口内作近似:平均值不随时刻改变、方差有限,且协方差只随两次读数的间隔改变。这样的二阶统计描述称为弱平稳,它不表示每条风速曲线不变。若暴风发展使平均值明显升高,就不能直接沿用这一假定。本章先读懂“大小”和“节奏”;频谱及随机风响应在第9章继续展开。
2.2.7.5 一段时间内:零次、恰好一次与至少一次
先选定统计单位。例如,把每天的最大风速与同一阈值比较,每天只有“超越”或“不超越”两种结果;也可以按年取最大值。每日概率与年度概率不能混用。 下图可切换“天”和“年”,切换单位只改变示例的时间单位,不会把每日概率自动换算成年概率。
设每个统计单位的超越概率为 ,各单位独立且概率保持不变,共观察 个单位。图中的三个区域互不重叠,合起来覆盖所有可能;它们来自二项分布的次数计算:20
| 图中区域 | 先用文字读 | 再写成概率 |
|---|---|---|
| 零次 | 每一次都不超越 | |
| 恰好一次 | 任选一次超越,其余全部不超越 | |
| 至少两次 | 剩下的所有情况 |
至少一次包含恰好一次、两次、三次…… 因而只需从全部可能中扣除“零次”:
这个相减求的是“至少一次”,不是“只超越一次”。默认按年计算,年超越概率1%、50年内:零次约60.50%,恰好一次约30.56%,至少两次约8.94%,所以至少一次约39.50%。这里计数的是超越年份,一个超越年份内仍可能有多次超过阈值的过程。若按天,则公式中的 必须是每天最大值的超越概率, 必须是天数;实际相邻天可能相关,乘法关系要先检查适用条件。
年超越概率 对应重现期 年。“百年一遇”表示年超越概率1%,不是每隔100年必定发生一次;它描述作用阈值被超过的机会,也不是桥梁失效概率。5
同一示例中,超越年份数的期望为 个年,与“至少一次的概率39.50%”是不同的量。期望表示次数的平均,不必为整数;次数的期望相加不要求独立,而上述“全不超越”概率相乘需要独立条件。
若把阈值写成 ,单次不超越的概率记为 ,那么独立、同分布的 次观测的最大值不超过 ,等价于每次都不超过 ,故 。这仍是图中“零次超越”的同一个道理。连续风速的每个采样点不能直接当成一次独立抽样,年最大值也不能直接套用瞬时风速的分布;进一步内容见统计、重现期与风险。
2.2.8 优化的变量、目标与约束
斜拉桥调索把每根索的张力作为变量,把主梁弯矩或挠度作为目标,并限制索力与支承反力。先在后面的索张力、分力与主梁平衡中看清作用关系,再比较两个优化目标。完整交互只在该处展开;数学上对应“变量—目标—约束”的二次优化。
描述空间分布, 描述一点的运动,振型给出相对比例,随机过程描述随时间变化的不确定性,优化选择符合目标和约束的状态。后面的构件分析继续给这些数学对象以具体位置。
2.3 梁与板:从荷载看内力、应力和变形
梁桥主梁、拱桥桥面梁、斜拉桥主梁和悬索桥加劲梁,都通过梁的弯曲、剪切和扭转参与工作。桥面板则把局部车轮作用传向纵梁与横梁。本节沿同一条关系分析:荷载作用在何处,截面传递什么,材料怎样受力,构件怎样变形。
2.3.1 同一片梁中的四种量
梁受到向下荷载时,支座反力在边界托住它。截取其中一段,截面剪力平衡竖向力,弯矩平衡这些力的转动作用。截面上又通过分布应力实现这两个合量:上部纵向纤维受压、下部纤维受拉,形成内力偶;腹板和相邻材料通过剪应力传递剪力。
内力作用于整个截面,应力属于截面内的材料点。以轴力为例,;弯矩是正应力对相应轴的力矩积分。相同弯矩下,加大截面高度可改变应力分布,不能把弯矩数值直接当作应力大小。位移又是沿构件变形累积的结果,截面弯矩较大不等于该处位移必然最大。
单独打开:梁中的物理量。剖切位置可沿梁移动,截面旁显示对应的剪力、弯矩和受拉受压区,变形线标示该点位移。跨中剪力为零而弯矩与挠度较大;靠近简支端,弯矩趋于零而剪力较大。截面合力、材料点应力和轴线位移分别对应三个观察尺度。
这里采用等截面、均质、未开裂、线弹性简支梁;两端支承中心距30 m,材料与截面参数在例2-3中给出。为了显示变形,图形可以放大,实际位移以读数为准。三维形体帮助定位构件和物理量,计算仍按已声明的一维梁理论进行。
2.3.2 有限元:把桥梁分成能共同工作的单元
当主梁连续多跨、桥面板有多种支承,或梁塔索连在一起时,难以用一条现成公式求出全部响应。有限单元法把对象分成有限个单元,用节点的位移和转角描述运动,再把各单元的平衡与连接协调组装起来。 网格中的线、面、体代表不同的力学近似。
| 单元与几何表示 | 在桥梁中用来算什么 | 必须保留或检查什么 |
|---|---|---|
| 梁单元:以轴线为几何 | 主梁、墩柱、拱肋的整体轴力、弯矩、剪力、扭转及位移 | 截面面积、惯性矩、扭转刚度、方向和端部释放;轴线虽细,截面并未消失 |
| 杆、索单元:沿连接线 | 拉索、吊杆、系杆的轴向作用与线形 | 拉索需考虑初拉力、只受拉和几何变化;普通双向受力杆不能自动代表柔索 |
| 板、壳单元:以中面划网格 | 桥面板的双向弯曲、箱梁板件与局部屈曲 | 厚度、法向、边界和面内连接;板强调弯曲,壳还表达面内作用 |
| 实体单元:在体积内划网格 | 锚固区、支座附近、节点局部三维应力 | 体积材料与接触边界,足够的厚度方向分辨率;全桥全用实体不一定更适合问题 |
单元内部用插值把有限个节点值连成近似位移场;材料关系给出应力或截面内力,平衡给出节点力。共用节点或明确的连接约束使相邻单元协调,在线性静力模型中,组装后得到
这里 是待求节点位移及适用的转角, 是与之配对的节点力和力矩。施加支承条件后求 ,再回到各单元计算内力、应力和反力。加入质量与阻尼就得到动力方程;在无阻尼自由振动下寻找能保持比例的运动,就接回前面的振型特征问题。16
有限元可以回答复杂构件怎样分担荷载、哪里变形较大、哪些模态容易参与、调索怎样改变内力。它也保留建模者的假定:错连了支座、把铰变成刚接、遗漏索初拉力,网格再密仍在求错误的问题。实际计算先核对整体平衡和边界,再与简单梁等已知解对照,并逐步加密观察关心的结果是否稳定。局部应力峰还要区分物理应力集中与理想尖角造成的数值奇异。第10章再展开模型建立、求解与核验。
2.3.3 材料放在哪里:面积、惯性矩和截面模量
面积 表示用了多少材料;截面二次矩 表示材料离所讨论弯曲轴多远。相同面积,把材料移离形心轴,通常可以提高该方向的抗弯刚度。但腹板还要传剪、连接上下翼缘并维持形状,不能因为它靠近中性轴就全部移走。
对同一种材料组成的分块截面,先由 求形心,再用平行轴定理
累加各块对总形心轴的贡献。 是两条平行轴的距离。面积以m²计,惯性矩以m⁴计;形心位置、弯曲轴或长度单位错了,应力与挠度会一起出错。
以 向上、受拉应力为正、正弯矩 使下缘受拉,均质线弹性截面的单向弯曲应力为
设形心至上下缘的距离为 ,截面模量为 、,则上下缘应力分别为 、。非对称T形截面的 通常不同,不能只用同一个“半梁高”计算两缘应力。截面开裂、不同材料组合或局部扰动明显时,需要改变截面模型。
例2-2 同一根矩形梁换个方向。 截面为 。以0.60 m为竖向高度时,、;转90°后,、。两者面积都为0.18 m²,自重相同。
保持材料、跨度和竖向荷载相同,竖放梁的抗弯刚度是横放的4倍,弯曲挠度为横放的四分之一;但截面模量只增加2倍,边缘应力仅减为二分之一。惯性矩解释变形,截面模量解释给定弯矩下的边缘应力。 两者来自同一截面几何,却回答不同问题。
图中仅比较弹性弯曲。实际主梁还须具有相应的剪切、扭转和横向稳定能力,梁高也受建筑高度、运输和施工条件限制。T梁、工字梁与箱梁的材料分布应随这些任务选择。
2.3.4 从弯矩到挠度:保持同一个算例
弯矩 经截面刚度 确定曲率;曲率沿梁累积,并满足支承条件,才确定挠曲线。对于相同材料和截面,跨度或支承改变,构件的力—位移关系也会改变。因此 是截面抗弯刚度,某一点的等效平动刚度 则还取决于体系与加载位置。
例2-3 一片30 m简支梁的连续计算。 取 、、、。材料重度25 kN/m³,形成自重100 kN/m;另加均布荷载20 kN/m,总线荷载 。采用本节三维梁的同一组参数与假定。
由整体平衡及截面平衡,
以m计, 为kN, 为kN·m。三处读数如下。
| 位置 | 剪力 | 弯矩 | 先看见的关系 |
|---|---|---|---|
| 左支点内侧 | 1800 kN | 0 | 简支端不传弯矩,梁段仍传剪力 |
| 四分点 m | 900 kN | 10125 kN·m | 剪力逐渐减小,弯矩继续增加 |
| 跨中 m | 0 | 13500 kN·m | 弯矩的变化率为零,达到最大值 |
跨中 ,上下缘应力分别为 MPa和 MPa。将弯矩关系代入 ,两次积分,并用 确定常数,可得
这个结果可以直接核对三维图:两端位移为零;中点坡度为零;两端反力之和为3600 kN;跨中上下缘应力绝对值相等。对于本例的静定梁,若只把 加倍,反力、剪力和弯矩不变,挠度减半。荷载加倍则使反力、内力、应力和挠度同时加倍。
只把计算跨径从30 m增到36 m,保持 不变,最大弯矩变为1.44倍,挠度变为2.0736倍,约24.121 mm。这不是“跨径增加20%,所有响应增加20%”。如果同时加高梁截面,就要重算惯性矩与自重,不能继续沿用这一倍率。
2.3.5 用能量理解整根梁怎样共同贡献位移
一小段梁弯曲时储存应变能,全梁的弯曲应变能为
它说明高弯矩区、低刚度区会产生较大变形贡献。要找某一点的位移,可以在该处沿目标方向虚设单位力,得到单位力弯矩图 ,再用单位荷载法
其中 表示单位力引起的弯矩系数,不能与动力学质量混淆。真实荷载没有因此改变;虚设的单位力只用来选择想要读取的位移。变截面梁按各段 累积,就能看到不同部位刚度对同一位移的贡献。
对例2-3,在跨中设向下单位力,左半跨 。利用对称性,
这个积分与微分方程解一致,说明11.633 mm由整根梁的曲率分布共同形成。若构件很短或很深,剪切变形不能忽略,总位移还需计入剪切项;杆系中的轴向、弯曲、剪切和扭转贡献可概括为
小写量均来自相应的单位荷载状态; 为已包含截面分布修正的有效剪切面积, 为自由扭转常数。选择保留哪些项,应由构件的主要变形和所需精度决定。
2.3.6 剪切、扭转与板的共同工作
弯矩和剪力由相同荷载产生,但在截面上通过不同应力传递。均质线弹性梁中,常用
表示弯曲剪应力。 为计算位置一侧面积对中性轴的静矩, 为该处有效宽度。薄壁构件常写成剪流 ,单位为N/m,它表示沿薄壁单位长度传递的剪切合力。翼缘与腹板之间的连接因此不能只按外形连接,还必须传递相应的剪流。
荷载偏离剪心时,会产生附加扭矩。以大小为 的竖向力、横向偏心 为例,外加扭矩为 。梁内各段扭矩还由支承抗扭条件决定,不能把外加集中扭矩当成全梁每个截面的内扭矩。
对于单室闭口薄壁截面,按Saint-Venant自由扭转近似,
为壁厚中线围成的面积。环向剪流依赖闭合路径;将箱梁顶板沿纵向切开,原有闭合条件便消失,不能保留原来的 。扭转常数也不能随意用极惯性矩代替。开口截面的翘曲、箱形截面的畸变,是需要分别辨认的变形。
桥面板受轮载时,沿纵横两个方向同时弯曲,板的扭曲也参与平衡。把板画成许多条带有助于辨认方向,但相邻条带仍相互连接,其变形需要协调。板模型适合描述轮载扩散、横向弯曲与板边约束;梁轴线模型只能回答它保留的整体截面问题。支座、锚固、横隔板和开孔附近,还需按局部构造处理作用扩散。
钢正交异性桥面板的纵肋和横梁使两个方向的结构刚度不同,这一特征来自构造布置;不能据此把钢材本身直接看成正交异性材料。箱梁的整体扭转是截面整体转动,畸变是截面形状发生变化,横向板弯曲又是局部板的运动。看三维图时先辨认哪一部分在怎样变形,再决定需要梁、板还是壳的表示。
2.3.7 预应力:在车辆到来以前建立状态
普通钢筋与混凝土在共同变形中承担内力;预应力则通过张拉和锚固,先在梁中引入受控的压力与弯曲。张拉后梁可能上拱,截面已有应力,后续自重、铺装和车辆作用都从这个状态继续增加。预应力没有改变平衡原则,而是改变了结构的起始状态。
仍取 向上、受拉应力为正。有效预压力 作用在形心下方,偏心距 ;正弯矩 使下缘受拉。对均质线弹性完整截面,
第一项是均匀压应力,第二项使下缘预压增大、上缘预压减小,第三项为外荷载弯曲应力。在例2-3跨中,若取 kN、 m,则均匀压应力为3.0 MPa,偏心项的边缘幅值为1.65 MPa。预应力单独作用时,上下缘分别为 MPa、 MPa;叠加本例荷载后,上缘约 MPa,下缘约 MPa。
下缘几乎不再受拉,只是这一截面、这一有效预压力和这一荷载状态的结果。传力时混凝土龄期不同,使用时预应力会损失,连续体系还会出现约束引起的次内力,均需按相应阶段分析。配筋与承载力计算在后续梁桥章节展开。
曲线束的张力方向沿梁改变,因此对梁产生横向作用。以钢束向下偏心 和等效荷载向下为正,在小斜率且张力近似恒定时,
束形二阶导数在这里重新出现:钢束改变方向越快,转向处需要传递的横向作用越大。锚端或折点同时存在集中作用,必须与分布作用一起保留。
例2-4 抛物线束的上托作用。 取 m、跨中向下偏心 m、两端位于形心、有效预压力 kN。束形为
负号表示向上的均布作用。每端切线斜率大小约0.10,端部对梁的向下分力约800 kN;全跨上托力 kN,恰与两端向下力平衡。只画上托力而删去端部作用,会把结构内部预应力误算成额外向上的外力。
2.3.8 施工顺序怎样留在内力中
例2-5 同一成桥外形,两种加载历史。 两跨等跨等刚度梁,每跨20 m。每跨先以简支状态承受自重 kN/m,随后在已有变形状态无强迫地连接成两跨连续梁,再施加两跨均布荷载 kN/m。支座无沉降,忽略接头强迫变形和长期时变效应。
简支阶段的跨中弯矩为2500 kN·m,中支点弯矩为零;后加荷载在连续体系中产生跨中弯矩500 kN·m、中支点弯矩 kN·m。最终跨中为3000 kN·m,中支点为 kN·m。
若把全部70 kN/m一次加到最终连续体系,得到的跨中弯矩为1750 kN·m、中支点为 kN·m。这是另一个加载过程,不能代替先简支后连续。合龙顶升、强迫转角、收缩徐变与预应力还会继续改变状态,具体过程见第4章的阶段施工。
回到梁的观测,计算挠度偏小可能来自刚度偏大,也可能来自支座约束、加载位置或测点不同。应先对应同一构件、同一状态、同一位置与方向,再比较反力、应变和挠度;若再测频率,还需对应质量与振型。网格加密可以减小离散误差,却不能改正单位、支座或阶段设置错误。核对贯穿每个算例,不需要另找一套脱离构件的判断规则。
2.4 受压构件:偏心、二阶效应与稳定
桥墩、索塔和拱肋常在承担较大压力的同时受到弯矩。斜拉桥的主梁也可能因索的水平分力形成压力。这些构件不能仅按“纯受压”理解:轴力与弯矩共同决定截面应力,压力又会放大已有横向变形。分析因此从截面走向整根构件,再走向相互连接的体系。
2.4.1 偏心压力怎样成为轴力与弯矩
压力作用线偏离截面形心时,可以等效为形心处的轴压力 与弯矩 。本小节以压力和压应力为正, 沿所讨论方向的偏心距。对关于所讨论弯曲轴对称的均质、完整、线弹性截面,两侧边缘的截面模量相同,边缘应力为
应力图由“均匀压应力”与“弯曲应力”叠加而成。增加偏心,一边压应力增加,另一边减小,合压力仍保持 。轴力受压不等于每一个材料点都受压;弯矩足够大时,某一侧可以出现拉应力。
例2-6 偏心增大,哪一侧先受拉? 矩形截面宽1.0 m、高1.2 m,压力2400 kN沿高度偏心0.15 m。、,两缘压应力为
即3.5 MPa和0.5 MPa,全截面受压。偏心增为0.25 m时,两缘为4.5 MPa和 MPa,负值表示线弹性模型出现拉应力。
要求全截面不出现拉应力,有 。对本例矩形截面,核心距 m。双向偏心时,矩形核心为菱形,满足
截面核心回答的是完整线弹性截面的无拉应力条件,不能直接用来区分钢筋混凝土构件的大、小偏心受压破坏,也不能代替开裂后截面抗力计算。圬工若不能承受计算拉应力,需要重建受压区和接触状态;钢筋混凝土则要考虑钢筋及混凝土的应变协调。
2.4.2 变形以后,压力的力臂也变了
一阶分析在变形前几何上列平衡,偏心压力产生弯矩 。构件向远离荷载作用线一侧弯曲后,压力对中间截面的力臂变成 ,弯矩随之增为
新增的 是二阶效应。它来自“已有压力乘以横向位移”,不是把荷载再加一遍。高墩、索塔、受压主梁、拱肋和施工中尚未形成完整横向联系的构件,都可能需要考虑这一机制。工程中常用 – 指体系侧移效应,用 – 指构件自身弯曲效应;两者都应由所取对象与实际变形解释。
采用等截面、线弹性、无初弯曲的偏心受压直杆模型,设两端压力在同一侧保持偏心 ,杆轴端点横向位移为零,端部允许转动。压力通过偏心连接传入轴线,所以杆端截面可有 弯矩。以增大压力力臂方向的横向位移 为正,曲率关系与变形后平衡连成
这是关于杆轴坐标的常微分方程。令 ,解得
这就是经典割线关系的位移和弯矩形式,可在马里兰大学力学教材的偏心压杆部分查阅其推导。7 它保留压力对横向变形的放大,同时仍采用线弹性、小曲率及指定端部条件;不能无限外推到大变形、屈服或接触改变后的状态。
单独打开:压力怎样放大弯曲。图中同时标出力的作用线与变形后的杆轴线。压力增大时,横向位移增加力臂,二阶弯矩比一阶弯矩 增长更快。页面可调整压力和偏心距,并显示对应位移与弯矩;端部偏心连接属于这一计算模型的边界。
例2-7 压力尚未达到欧拉值,弯矩已明显放大。 取 GPa、、 m、 mm。对应理想两端铰支直杆的欧拉力为2467.4 kN。压力取其0.45倍,即1110.3 kN;有 。
一阶弯矩为22.207 kN·m;割线放大系数约2.023,二阶最大弯矩约44.922 kN·m,跨中附加横向位移约20.458 mm。此处仅展示弹性放大机制。已知 而未给截面面积、边缘距离及材料强度,尚不能完成截面应力与承载力判断。
2.4.3 欧拉屈曲为什么也是特征问题
若直杆完全无初始弯曲、轴压力严格通过形心,直线状态在任意荷载下形式上都满足横向平衡。问题在于:这个平衡状态能否抵抗微小侧向扰动。对理想两端铰支细长弹性杆,非零侧向形态满足
边界条件只允许特定荷载值出现非零解:
最小正特征值给出基阶欧拉临界力 。形态中的幅值 不能由线性特征问题确定,与振型一样只表示相对形状;实际缺陷、材料非线性和加载过程决定实际位移与承载过程。
对于可用有效长度表述的理想杆,写成 。例2-7中只把有效长度从8 m降为4 m,理论临界力提高为4倍,即9869.6 kN;前提是新增支撑确实限制所考察方向侧移,并有足够的刚度、强度与可靠连接。另一个弯曲方向、受压板件局部屈曲和侧扭失稳仍需分别检查。
将整体杆系离散,常见的屈曲形式为
这里把 定义为参考压荷载对应的几何软化矩阵, 为参考荷载的临界倍数。采用其他符号定义时式中的正负号会改变。它与自由振动的
形式相似,但右侧对应的物理量不同:前者由压力改变平衡稳定性,后者由质量与恢复力决定运动频率。屈曲模态不是振型,荷载倍数不是频率。 两者都需要先声明基态、边界和所采用的线性化。
实际压杆有初弯曲、残余应力和荷载偏心,往往在理论分岔前已有侧移;材料也可能先屈服。因此欧拉值和线性特征屈曲值用于理解与估计理想稳定性,不能直接当作真实构件的设计承载力。施工时尚未安装的横撑,不能提前计入构件约束。
用于压弯判断的 与 应来自同一截面、同一荷载状态。分别取不同状态的最大轴力与最大弯矩,不能自动组成一个真实受力状态。得到弹性应力后,再按材料与构件特点检查强度、开裂、稳定和变形;这些判断始终附着在具体构件上。
2.5 拱与索:曲线传力与索力调整
同样的竖向荷载,可以经向上拱起的构件传递,也可以经向下悬垂的索传递。两者都依靠水平力与曲线几何共同工作:拱主要承受压力,柔索承受拉力。支承必须能承担水平分力;曲线较平时,水平力会明显增大。
单独打开:索形、索力与拱推力。默认 m、 m、 kN/m,水平分力500 kN,端部索力640.3 kN。跨径与荷载不变时,垂度减小使水平分力增大;同形拱具有相同水平力大小,但主要承担轴向压力。页面可调整垂度并对照力的方向,读数均来自按水平投影均布荷载的理想模型。
2.5.1 拱的曲线怎样改变弯矩
把一根水平梁改为向上拱起,并使两端能够形成水平推力,竖向荷载引起的一部分弯矩便可由推力的力矩抵消。取同跨径、同竖向荷载的简支梁弯矩为 ,拱轴相对拱脚连线的竖向坐标为 ,水平推力为 ,有
它并未保证拱处处纯压。只有相应荷载下的压力线与拱轴相合,并满足所取模型的条件时,弯矩才可能消失;荷载分布、拱轴、支承和温度改变,都可能重新产生压弯。
例2-8 抛物线三铰拱的推力。 两拱脚等高,跨径 m,拱顶矢高 m,承受沿水平投影均布的竖向荷载 kN/m。由对称与整体平衡,每端竖向反力400 kN;拱顶铰的弯矩为零,因此
抛物线拱轴 恰与该荷载压力线一致,理想模型中 。拱脚轴向压力大小为 kN。若矢高减半、其他条件不变,水平推力加倍;基础或系杆必须能闭合它。
本节三维对照使用同一跨径、同一均布荷载与同一矢高绝对值,将向上的拱与向下的柔索并列观察。两者水平力公式相同,轴向力的性质相反;拱要特别关心受压稳定,索要保持受拉状态。
拱轴偏离压力线时,以相应截面压轴力 和弯矩 表示法向偏心,,便回到上一节偏心受压的截面应力分析。图上压力线与拱轴的竖向距离为 ,它一般不等于法向偏心。对倾斜截面, 与 也一般不同,不能拿竖向量图值直接代入核心条件。
2.5.2 拱轴怎样找到:把受拉索形倒置
给定支点和竖向荷载,柔索自然沿能够传递拉力的方向形成平衡线;把这条形状倒置,可以得到相应的理想受压线。对桥梁,先看桥面荷载怎样传到拱,再确定压力线,最后判断拟采用的拱轴是否接近它。18
单独打开:拱轴找形。全跨均布荷载下,抛物线与压力线相合;只加载半跨或改为集中力,适合的形状随之变化。保持原抛物线时,图中的两条线分开,拱便同时承担弯矩。把中铰高度调低,水平推力增大,基础需要承担的水平作用也随之改变。
这涉及三类设计思路。拓扑先回答哪些节点通过哪些构件连接,例如有无系杆、推力是否交给桥台;找形在给定连接、荷载和边界下确定合适曲线;拓扑优化则允许进一步改变材料分布或传力路径。当前图保持三铰拱连接,只改变荷载和轴线几何。三者相互关联,但不是同一个计算步骤。
2.5.3 柔索的受力特点
理想柔索主要承受拉力,不能像梁那样用弯矩维持任意形状。索的方向随位置改变,张力的竖向分量因此能够托住荷载;要改变垂度,所需水平分力也会改变。索的线形与索力共同满足平衡,给定其中一个,另一个通常就受到约束。
2.5.4 水平分力与竖向分力从哪里来
取两支点等高的柔索,以左支点为 , 表示相对支点连线向下的垂距。忽略索的抗弯刚度,所有分布荷载竖向作用,没有沿索分布的水平外力,因此张力水平分力 沿索保持不变。索斜率决定竖向分力,局部平衡给出
此处 按单位水平投影长度计。索力方向沿切线,张力 的大小通常沿索变化;只有水平分力 为常量。把整根索的张力都等同于跨中水平分力,会低估较陡位置的索力。
与梁的 对照,柔索依靠已有张力与曲率平衡横向荷载,梁依靠抗弯刚度与曲率的空间变化传力。二者都从微段平衡出发,控制关系却由构件的主要机制决定。
2.5.5 抛物线与悬链线取决于荷载怎样分布
若 沿水平投影均布,且 ,积分得
令跨中垂度为 ,便有
这是抛物线柔索。用许多竖向吊索承托近似均匀桥面荷载时,可以由它建立最初的力学认识;实际离散吊点、索自重、弹性伸长与不均匀荷载还会使索形改变。
若荷载来自索自身重量,并且每米实际索长的重量为常量 ,一段水平投影 上的重量为 ,因此
这不再是前一个常量右端的方程。令 ,两等高支点间的解可写为
它是悬链线,跨中垂度满足 。已知跨径、自重和垂度时,通常先求参数 ,再由 得水平力。MIT的悬链线静力讲义展示了这类由张力与自重确定线形的关系。8
小垂度条件下,悬链线可近似为抛物线,但二者的荷载定义不能交换。前者按实际索长均布,后者按水平投影均布;同时存在桥面作用与索自重时,应把各项按同一坐标正确组合。给定无应力索长、支点位置与材料刚度时,还需把索的弹性伸长纳入协调,不能同时任意指定索长、垂度和索力。
2.5.6 从索形和振动反推索力
索力不只可以由设计给定,也可以从测量反推。激光扫描先得到索的空间点云,提取中心线、锚点和垂度,再用索的平衡关系拟合;振动测量先得到索的运动记录,识别频率,再用索力与频率的关系估计张力。两条路线测到的都是响应或几何,需要力学模型才能变成索力。
激光扫描的原始点云并非一条干净曲线:遮挡、点密度不均、表面反射以及扫描期间风和车辆引起的运动,都会影响拟合。需先识别同一根索、统一坐标、确认锚点,再结合线重和适用的悬链线或索单元模型估计张力。已有桥梁研究将点云提取、力学模型拟合和频率法对照组成完整路线。17
对两等高支点、按水平投影均布荷载的理想索,测得垂度 后可以直接看到
相同跨度与荷载下,垂度小,水平分力大。这个式子中的 是水平分力;较陡位置的总张力仍要由索的斜率计算。斜拉索的两锚点不等高,其垂度通常相对锚点连线测量,不能直接把这个等高索公式中的垂度换成斜索图上的任意距离。
振动法先用长度 、恒定张力 、线质量 的理想直弦说明。两端横向位移为零,忽略垂度、抗弯刚度和阻尼,有
频率越高并不自动意味着索力越大,比较时还要保持索长、线质量、阶次及边界一致。 用m、 用kg/m、 用Hz,所得 为N。例如 m、 kg/m,测得基频1.0 Hz时,估计张力2000 kN;同条件基频提高到1.1 Hz,估计张力变为2420 kN,增加21%。
这一关系来自索微段的横向平衡:张力方向随索的弯曲而变,产生恢复作用,平衡微段惯性。将空间形态与时间运动分开,其表达为
这里 是沿索位置, 是横向位移;偏导只是分别考察“位置变化”和“时间变化”。前面的桥梁主要由抗弯刚度恢复,理想弦主要由已有张力恢复,所以本例 ,而均匀简支梁竖弯的 ,两种频率比不能混用。
实际短索的弯曲刚度、长索的垂度与伸长、锚端柔度、有效振动长度和减振器均会影响结果。3 应比较多个阶次的 ,并与线形、张拉记录或其他测量相互核对。某根索的局部频率也不能直接充当全桥竖弯频率。实际测量的关键链条是:测到什么—采用什么模型—估计哪个力—用什么独立证据核对。
2.5.7 调斜拉索,让主梁的弯曲和线形改变
在索上调张力,在主梁上看效果。 图中的青色内索与棕色外索分别标出每根张力 。拖动索上的圆点,同组左右两根索同步增减张力;选择一根索,锚点旁的箭头显示它怎样拉住主梁。向上的分力 托梁,指向塔的分力 压梁,主梁的线形与弯矩随之改变。
这座教学桥为60 m+60 m,跨中设塔,两端及塔位设竖向支承;两对对称拉索分别由内、外索力控制,默认每根1800 kN,允许范围0~4000 kN。手动增力时若出现承压支承负反力,应重新检查支承状态。
图中的量各有位置:输入是斜拉索张力 ,传给主梁的是锚点分力 、,输出是主梁位移 、弯矩 和轴力 。 灰色参考线保留未变形位置,线形显示放大300倍,实际位移以毫米读数为准;力的倾角仍按未变形的索计算。切换到弯矩时,弯矩图的零线贴着主梁轴线。
点击“相同竖向托力”,每根索都向上托住1500 kN。外索与水平线夹角45°,张力需要2121 kN;内索较陡,夹角约63.4°,张力只需1677 kN。同样向上托住主梁,越平的索需要越大的张力。 两索的竖向箭头等长,沿索箭头却不同;水平分力也分别为1500 kN和750 kN。再切到“相同张力”,各索张力都是1800 kN,竖向托力却不再相同。
“相同竖向托力”用于看清倾角的作用,并不就是最优索力。选择“减小全梁弯矩”或“减小全梁挠度”,计算还要考虑各个锚点对全梁的影响,所得内、外索力以及各索竖向分力都不必相同。实际调索通过张拉或索长调整实现,本图用张力作为直接输入,先呈现从索到梁的作用关系。
这些箭头对应的方程很短。设索与水平线的夹角为 ,则
向上托梁,直接调整竖向弯矩和挠度; 沿梁指向塔,使相应区段承受轴向压力。塔端同时受到方向相反的索力。仅画向上力会漏掉主梁受压与桥塔受力,这也是图中需要保持真实斜索方向的原因。
例2-9 同一座斜拉桥,两种调索目标。 取全长 m,塔位在60 m,塔顶比索锚基准高40 m;外索锚点在20、100 m,内索锚点在40、80 m。外索倾角45°,内索倾角约63.435°。主梁承受80 kN/m均布荷载, kN·m²,在0、60、120 m竖向支承,塔位主梁保持弯曲连续。
用两端简支梁作为计算基准,再以向上的塔位反力 保证跨中位移为零。 表示该基准梁在 处施加单位向下力时, 处的向下位移柔度。于是
这一步让“多一个支承”体现为明确的变形协调条件,而不是任意分配反力。对同一基准梁, 表示单位向下力的弯矩影响,故
两端反力仍由整体平衡求得:
左半桥任一截面的轴向压力,来自该截面外侧各拉索水平分力的累加;右半桥对称。若轴力以受拉为正,则该压力记为负值。在邻近塔的区段,两种索都参与压梁;最外索以外的区段则没有这部分索水平力。图中阶梯状轴力分布由此产生。
| 状态 | 每根内索 / kN | 每根外索 / kN | 两端主梁支承反力各 / kN | 塔位主梁支承反力 / kN |
|---|---|---|---|---|
| 默认调索状态 | 1800.000 | 1800.000 | 901.520 | 2031.438 |
| 弯矩平方积分最小 | 1720.052 | 2567.342 | 630.769 | 1630.769 |
图后的数学表达只需分清三项:
| 优化内容 | 这个桥上的对应对象 |
|---|---|
| 设计变量 | 两对对称拉索的每根张力 |
| 目标 | 减小全梁弯矩平方积分,或减小竖向位移平方积分 |
| 约束 | 每根索的力限值;两端及塔位的主梁承压支承反力不小于零 |
以减小弯矩为目标,可写为
固定几何、线弹性和相同支承条件下, 与 都随两组索力线性变化,平方积分便形成带约束的二次优化问题。4 本页目标是平方积分,不是直接最小化绝对最大弯矩;两种提法应当区分。
改为挠度平方积分最小,内索约1679.832 kN、外索约2623.748 kN,目标改变,解也改变。索力调节幅度不是全部同增同减,必须通过结构影响关系判断。表中的塔位反力是主梁支承处的反力;塔基础还要承担塔端拉索向下分力及实际塔自重等,不能把这两个反力混为一谈。
这个教学模型固定塔顶位置和索方向,把索力作为可调作用,暂不计塔柔度、索伸长协调与主梁二阶压弯。真实斜拉桥需在全桥平衡状态附近建立影响关系,通过索长或张拉控制量实施调整,重算梁塔索共同平衡,再检查各施工阶段。成桥控制在斜拉桥与悬索桥中继续展开。
从梁到压杆、拱和索,构件的物理机制逐步改变:梁以弯曲刚度抵抗横向作用,压杆的轴力又会放大弯曲,拱用曲线与水平推力调整内力,柔索靠拉力维持线形。数学工具反复出现,但每次都需要由对象、荷载与连接重新确定具体方程。
2.6 本章小结与练习
一张有效的计算图,应当让人看清研究对象、连接、作用和所作假定。数学基础把这些内容写成有条件的方程;构件分析再赋予每个量以具体位置:荷载在构件外部,内力在截面,应力在材料点,位移沿构件累积,振型表示相对运动。
同一套平衡、协调和材料关系,分别得到梁的弯曲、压杆的二阶效应、拱的推力与索的线形。特征值既能求振动频率,也能求理想屈曲荷载,但两者的物理意义不同。概率需要明确随机对象与时间窗口,调索优化需要明确目标和约束。读懂这些量以后,才能把构件组织成后续各章的桥梁体系。
| 本章用到的思路 | 在桥梁上怎样落地 |
|---|---|
| 分离与平衡 | 从整桥取出主梁,在切断处补作用,核对合力与合矩 |
| 简化与离散 | 按关心的机制选梁、索、板壳或实体;用节点连接共同工作 |
| 微分与累积 | 从位移读速度、加速度;从荷载到内力、曲率与挠度 |
| 特征分析 | 看清各阶振型及频率,区分振动与屈曲特征问题 |
| 随机描述 | 区分一条风速时程、同一时刻的分布、采样和相关性 |
| 找形、反演与优化 | 由荷载找拱轴,由索形或频率估索力,按主梁目标调整索力 |
AI与脚本可以加快计算、绘图和数据拟合。读者仍要能够在桥上指出每个量的位置,说明模型保留了什么,并判断参数变化为什么引起图中变化。图让变化可见,方程给变化以依据,平衡与独立测量用于检验结果。
2.6.1 本章练习
1. 从主梁到自由体(对应目标:1)。 另设24 m简支梁,180 kN集中力距左支点18 m,忽略自重。画出自由体,说明两端分别限制哪些自由度,求左右竖向反力,并用合力与合矩检查。若研究对象扩大到梁、支座和桥墩,原来的支座接触力应怎样处理?
2. 从受力到运动(对应目标:1、2)。 对 ,在桥梁振动图上指出位移、惯性、恢复与耗能所对应的对象。荷载幅值不变,频率从0.5 Hz移至接近1 Hz时,第9章共振算例的响应怎样变化?保持1 Hz并增大阻尼比时又怎样变化?再用简支端与固定端的图示,说明边界条件和 所表示的初始条件有何不同。
3. 振型不是振幅(对应目标:2)。 均匀简支桥的一阶振型在四分点、跨中、四分之三点的比例为 ,某时刻跨中下移10 mm,求三个点的位移。若振型整体乘2,模态坐标怎样改变才能保持实际运动?已知本例一阶频率1 Hz,列出前四阶竖弯频率,并说明此比例为何不能通用于所有桥型。
4. 概率、次数与时间相关(对应目标:2)。 某阈值的年超越概率为1%,各年独立且概率不变。求50年内至少超越一次的概率与超越年份数的期望,解释二者的区别及其为何不等于结构失效概率。两条风速过程有相同均值和方差,能否据此认定它们引起相同的振动?
5. 梁的截面与已有状态(对应目标:3)。 0.30 m×0.60 m矩形简支梁由竖放改为横放,保持材料、跨度和竖向外荷载相同,求横放与竖放的弯曲挠度比和边缘应力比。再对例2-4的20 m抛物线预应力束,仅把有效预压力减小15%,求均布上托力和每端竖向分力,检查合力是否平衡。说明这两组计算分别改变了什么。若分别计算全梁整体弯矩、桥面板局部弯曲、锚固区三维应力,应优先考虑哪类有限元单元?
6. 截面受拉与构件失稳(对应目标:4)。 矩形截面宽0.8 m、高1.5 m,压力1800 kN沿高度偏心0.30 m。按完整均质线弹性截面、压应力为正,求减压侧边缘应力。另有一根理想轴压杆,保持 不变,有效长度降为原来的0.8倍,求欧拉临界力倍率。解释为何前一问与后一问不能互相代替。
7. 索形、索力与频率(对应目标:5)。 两等高支点间的柔索跨度40 m,承受沿水平投影均布的20 kN/m竖向荷载,跨中垂度由8 m改为4 m。求改变前后的水平分力和最大张力。另设有效长度100 m、线质量50 kg/m的理想直弦,测得基频1.0 Hz,估算张力;指出实际拉索至少两项可能使该估算偏离的因素。
8. 写清楚“最优索力”(对应目标:1、5)。 使用本章60 m+60 m斜拉桥算例,说明两组变量指的是实际索张力还是竖向分力,写出减小弯矩的目标及页面约束。对比弯矩与挠度两个目标下的索力,说明为何会不同。选一个可行状态,用竖向合力与合矩核对反力,再指出主梁轴向压力来自哪里。
参考答案与评价要点
第1题。 kN、 kN; kN, kN·m。平面铰支端限制两平动,活动端限制竖向位移、释放沿梁位移,两端允许转动。扩大对象后,梁与支座之间的接触力成为内部作用,新的外边界作用仍须补画。
第2题。 是所选坐标的位移,、、 分别是惯性、阻尼和弹性恢复项,右侧为外荷载。例中固有频率1 Hz,频率由0.5 Hz接近1 Hz时,稳态幅值明显增大;5%阻尼比、1 Hz处约25.330 mm。阻尼比增至10%,同处降为12.665 mm,荷载幅值仍为10 kN。简支端限制竖向移动、允许转动,固定端同时限制位移与转角;这些约束在运动过程中持续成立。初始位移与初始速度规定运动从什么状态开始,不改变支承本身。
第3题。 三点位移约为 mm。一阶模态坐标为10 mm;振型整体乘2,坐标减半为5 mm,实际位移与频率不变。本例前四阶频率1、4、9、16 Hz,来自均匀简支梁的 ;质量、刚度、支承或结构体系改变,频率比就可能改变。
第4题。 至少超越一次的概率 ;超越年份数期望为 ,即平均0.5个超越年。一个表示是否发生过的机会,一个表示次数的平均值。结构失效还需由作用计算响应并与抗力或性能边界比较。同均值、同方差的过程仍可能有不同时间相关和频谱,不能据此认定动态响应相同。
第5题。 竖放、横放的惯性矩分别为0.00540与0.00135 m⁴,截面模量为0.0180与0.0090 m³;横放/竖放挠度比为4,边缘应力绝对值比为2。预压力降为6800 kN,均布上托力为68 kN/m,每端向下分力680 kN;。前一组只改变截面方向,后一组只改变已有预压力,不能将两种原因混在一起。整体弯矩可选梁单元,桥面板局部弯曲可选板壳单元,锚固区三维应力可选实体单元;单元选择仍需与连接、荷载和结果精度要求配套。
第6题。 m², m³,减压侧应力 MPa,表示线弹性计算出现拉应力。欧拉临界力倍率为 。前一问是给定截面内力下的材料点应力,后一问是理想构件维持平衡的稳定问题。实际设计还需结合材料、缺陷和边界,核心条件不能代替稳定检查。
第7题。 垂度8 m时 kN,最大张力640.3 kN;垂度4 m时 kN,最大张力1077.0 kN。每端竖向分力始终为400 kN,因此索变平时水平力显著增加。独立直弦例的张力为 N,即2000 kN;垂度、弯曲刚度、有效长度、锚端柔度或减振器等均可能造成偏离。
第8题。 变量是每根实际斜索的张力 ,竖向分力为 。目标 ;每根力限0~4000 kN,且三个主梁支承反力非负。弯矩目标得内索1720.052 kN、外索2567.342 kN,挠度目标得内索1679.832 kN、外索2623.748 kN。两个目标衡量的响应不同,所以解不同。核对 kN及 kN·m。主梁压力来自拉索水平分力 ,由每个截面外侧力的累加确定。
AI可以辅助编写计算脚本、绘制内力与变形图、整理测量数据。计算文件应同时保留对象、参数、单位、假定和所用关系,使图形与数值能够追溯到同一模型。结论的依据是平衡、边界、量纲及独立算例的一致性;三维图清楚展示物理量,计算记录说明结果为何成立。
本章公式用于建立构件分析基础。桥梁作用与组合、钢筋混凝土与预应力截面、钢构件稳定的工程设计,分别需结合项目采用的规范与构造要求。规范阅读按本章问题给出查阅入口。
本章参考文献
UNIVERSITY OF MANCHESTER. Seeing and Touching Structural Concepts[EB/OL]. [2026-10-19]. 项目目录. 参考范围:通过对象、变形与受力观察建立结构概念。
MASSACHUSETTS INSTITUTE OF TECHNOLOGY. 2.002 Mechanics and Materials II: Euler–Bernoulli Beam Vibration[EB/OL]. 2004[2026-10-19]. 课程讲义. 参考范围:梁的动力方程、边界与分离变量。
FEDERAL HIGHWAY ADMINISTRATION. Mitigation of Wind-Induced Vibration of Stay Cables: Numerical Simulations and Evaluations: FHWA-HRT-14-049[R/OL]. 2014[2026-10-19]. 第2章:理论背景. 参考范围:理想弦、弯曲刚度和垂度对索振动的影响。
BOYD S, VANDENBERGHE L. Convex Optimization[M/OL]. Cambridge: Cambridge University Press, 2004[2026-10-19]. 作者公开版本. 参考范围:约束最小二乘与二次优化;本章调索算例为独立编制的教学模型。
U.S. GEOLOGICAL SURVEY. The 100-Year Flood[EB/OL]. 2018[2026-10-19]. 年超越概率与重现期.
NIST/SEMATECH. e-Handbook of Statistical Methods: Autocorrelation[EB/OL]. [2026-10-19]. 自相关说明. 极值资料的抽样与时间相关另见NIST极端风速统计说明。
DALLY J W, BONENBERGER R J, FOURNEY W L. Mechanics I: Statics+++[M/OL]. 2016[2026-10-19]. University of Maryland课程版本. 参考范围:第7章偏心受压直杆与割线关系。
MASSACHUSETTS INSTITUTE OF TECHNOLOGY. 2.019 Design of Ocean Systems: Mooring Dynamics, Continued[EB/OL]. 2011[2026-10-19]. 悬链线静力讲义. 参考范围:沿实际索长均布自重的悬链线关系。
肖汝诚, 等. 桥梁结构体系[M]. 北京:人民交通出版社, 2013. 参考范围:以构件、连接与约束认识桥梁体系;本章算例均另列假定与参数。
中华人民共和国交通运输部. 公路桥涵设计通用规范:JTG D60—2015[S/OL]. 2015. 官方发布页.
中华人民共和国交通运输部. 公路钢筋混凝土及预应力混凝土桥涵设计规范:JTG 3362—2018[S/OL]. 2018. 官方发布页.
中华人民共和国交通运输部. 公路钢结构桥梁设计规范:JTG D64—2015[S/OL]. 2015. 官方发布页.
MASSACHUSETTS INSTITUTE OF TECHNOLOGY. 2.003SC Engineering Dynamics: Linear System Modeling, a Single Degree of Freedom Oscillator[EB/OL]. 2011[2026-10-19]. 课程演示与讲义. 参考范围:单自由度模型、初始响应、阻尼与简谐稳态响应。
DELFT UNIVERSITY OF TECHNOLOGY. Computational Modelling: Modal Analysis[EB/OL]. [2026-10-19]. 交互教材. 参考范围:质量与刚度矩阵、广义特征值及模态叠加。
NREL. TurbSim User’s Guide, v2.00[R/OL]. [2026-10-19]. 官方项目手册. 参考范围:谱表示法与式(23)Kaimal型谱;本章另设教学参数及有限频带。
DELFT UNIVERSITY OF TECHNOLOGY. Computational Modelling: 2D Frame Analysis[EB/OL]. [2026-10-19]. 杆梁单元、节点自由度与组装. 单元类型另见Dassault Systèmes Abaqus单元选择文档。
3D Laser Scanning-Based Tension Assessment for Bridge Cables Considering Point Cloud Density[J/OL]. Structural Control and Health Monitoring, 2024: 8094924. 论文原文. 参考范围:桥梁索点云、线形提取、力学拟合及与频率法的比较。
ETH ZÜRICH, BLOCK RESEARCH GROUP. eQUILIBRIUM: Funicular for Vertical Forces[EB/OL]. [2026-10-19]. 竖向力的平衡索形交互. 参考范围:荷载、水平推力、边界与平衡线的关系;本章三铰拱数值由所列模型独立计算。
NIST/SEMATECH. What is a Probability Distribution[EB/OL]. [2026-10-19]. 概率密度与曲线下面积。参考范围:连续变量的区间概率、密度与累计分布。
NIST/SEMATECH. Binomial Distribution[EB/OL]. [2026-10-19]. 二项分布。参考范围:独立、同概率试验的成功次数、零次与恰好一次概率及次数期望。