桥梁结构地震反应分析学习教案.pptx

《桥梁结构地震反应分析学习教案.pptx》由会员分享，可在线阅读，更多相关《桥梁结构地震反应分析学习教案.pptx（81页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、会计学1桥梁结构地震反应桥梁结构地震反应(fnyng)分析分析第一页，共81页。4.1 引言引言(ynyn)n n位于地震区的公路桥梁，在遭遇地震作用时，结位于地震区的公路桥梁，在遭遇地震作用时，结构将发生构将发生(fshng)(fshng)振动，从而使结构产生随时振动，从而使结构产生随时间变化的位移、速度、加速度、内力和变形等，间变化的位移、速度、加速度、内力和变形等，通称为结构的地震反应。通称为结构的地震反应。n n桥梁结构地震反应的大小，除了与地震地面运动桥梁结构地震反应的大小，除了与地震地面运动有关外，还与结构自身的动力特性（自振周期与有关外，还与结构自身的动力特性（自振周期与阻尼）有

2、关。阻尼）有关。n n所谓桥梁结构地震反应分析，主要是计算桥梁结所谓桥梁结构地震反应分析，主要是计算桥梁结构在地震地面运动作用下的内力和变形。构在地震地面运动作用下的内力和变形。第1页/共81页第二页，共81页。n n桥梁抗震设计，首先要解决桥梁结构地震反应的计算问题。地震地桥梁抗震设计，首先要解决桥梁结构地震反应的计算问题。地震地面运动作为动态作用，其引起的桥梁结构反应遵循一般的结构动力面运动作为动态作用，其引起的桥梁结构反应遵循一般的结构动力学原理。但由于地震地面运动有别于一般的动力荷载学原理。但由于地震地面运动有别于一般的动力荷载(hzi)(hzi)，地，地震地面运动作用下桥梁结构的动力

3、反应分析必然有它的特殊之处。震地面运动作用下桥梁结构的动力反应分析必然有它的特殊之处。n n桥梁结构地震反应分析方法可分为解析法和数值法两类。解析法建桥梁结构地震反应分析方法可分为解析法和数值法两类。解析法建立在对结构充分简化的基础之上，从目前规定的分析方法来看，普立在对结构充分简化的基础之上，从目前规定的分析方法来看，普通简支梁桥和拱桥的地震力计算方法仍是基于解析法。通简支梁桥和拱桥的地震力计算方法仍是基于解析法。第2页/共81页第三页，共81页。n n解析解：一种包含分式、三角函数、指数、对数等基本解析解：一种包含分式、三角函数、指数、对数等基本函数的解的形式。用来求解析解的方法称为解析法

4、函数的解的形式。用来求解析解的方法称为解析法analytic techniquesanalytic techniques、analytic methodsanalytic methods。解析法即常见。解析法即常见的微积分技巧，如分离的微积分技巧，如分离(fnl)(fnl)变量法等。解析解为一封变量法等。解析解为一封闭形式闭形式closed-formclosed-form的函数，对任一独立变量，将其带的函数，对任一独立变量，将其带入解析函数可求出精确的相依变量。入解析函数可求出精确的相依变量。n n 数值解：当无法由微积分技巧求出解析解时，只能利用数值解：当无法由微积分技巧求出解析解时，只能利

5、用数值分析方式求解。在数值分析过程中，首先将原方程数值分析方式求解。在数值分析过程中，首先将原方程式加以简化，如先将微分符号改为差分符号等；再用传式加以简化，如先将微分符号改为差分符号等；再用传统的代数方法将原方程式改写成另一方便求解形式。统的代数方法将原方程式改写成另一方便求解形式。n n 求解步骤：求解步骤：n n 将一独立变量带入，求得相依变量的近似解。求得的将一独立变量带入，求得相依变量的近似解。求得的相依变量为一个个分离相依变量为一个个分离(fnl)(fnl)的数值的数值discrete valuesdiscrete values，不像解析解为一连续的分布，经过上述简化，其精确性不像

6、解析解为一连续的分布，经过上述简化，其精确性不如解析法。不如解析法。第3页/共81页第四页，共81页。第4页/共81页第五页，共81页。第5页/共81页第六页，共81页。4.2 4.2 桥梁结构地震桥梁结构地震桥梁结构地震桥梁结构地震(dzhn)(dzhn)反应分析方法简介反应分析方法简介反应分析方法简介反应分析方法简介分析方法的演变过程分析方法的演变过程自从自从1899 1899 年日本学者大森房吉首次提出用于抗震设计的静力法以来，年日本学者大森房吉首次提出用于抗震设计的静力法以来，桥梁结构地震反应桥梁结构地震反应(f(f nyng)nyng)分析方法历经了从静力法到动力的反分析方法历经了从

7、静力法到动力的反应应(f(f nyng)nyng)谱法和动态时程分析法的演变过程。谱法和动态时程分析法的演变过程。依据输入地震动的特点，桥梁结构地震反应依据输入地震动的特点，桥梁结构地震反应(f(f nyng)nyng)分析方法可以分析方法可以分为两大类：确定性方法和随机振动方法。分为两大类：确定性方法和随机振动方法。第6页/共81页第七页，共81页。n n确定性方法使用天然地震地面运动记录或人工模拟确定性方法使用天然地震地面运动记录或人工模拟地震地面运动作为地震动输入，求出结构的反应。地震地面运动作为地震动输入，求出结构的反应。确定性方法可进一步分为确定性方法可进一步分为(fn wi)(fn

8、 wi)静力法、反应谱静力法、反应谱法和动态时程分析法。法和动态时程分析法。n n随机振动方法把地震视为随机过程，把具有统计性随机振动方法把地震视为随机过程，把具有统计性质的地震动作用在结构上，求出结构的反应。质的地震动作用在结构上，求出结构的反应。n n迄今为止，绝大多数国家现行的结构抗震设计规范迄今为止，绝大多数国家现行的结构抗震设计规范均采用确定性方法。只有欧洲规范（均采用确定性方法。只有欧洲规范（Eurocode 8Eurocode 8）允许使用确定性方法或随机振动方法。允许使用确定性方法或随机振动方法。第7页/共81页第八页，共81页。n n随机振动：未来任一给定时刻的瞬时值不能预先

9、确定的振动，无法随机振动：未来任一给定时刻的瞬时值不能预先确定的振动，无法用确定性函数而须用概率统计方法定量描述其运动规律。用确定性函数而须用概率统计方法定量描述其运动规律。如：车辆如：车辆在高低不平路面上行驶、桥梁在高低不平路面上行驶、桥梁/高层建筑在阵风或地震作用下发生高层建筑在阵风或地震作用下发生的振动就是随机振动。的振动就是随机振动。n n随机振动的单次试验结果有不确定性、不可预估性和不重复性，但随机振动的单次试验结果有不确定性、不可预估性和不重复性，但相同条件下的多次试验结果却有内在的统计规律。将每次试验结果相同条件下的多次试验结果却有内在的统计规律。将每次试验结果都看作一个样本，则

10、它们全体（集合）构成一个随机过程，用以表都看作一个样本，则它们全体（集合）构成一个随机过程，用以表示随机振动的响应。对一随机过程示随机振动的响应。对一随机过程 ，任意时刻的样本取值是随机的，任意时刻的样本取值是随机的，称随机变量。称随机变量。n n用概率统计方法可得到随机过程的各种信息，如用数学用概率统计方法可得到随机过程的各种信息，如用数学(shxu)(shxu)期期望表示随机变量的平均值；用均方值表示随机变量平方的平均值；望表示随机变量的平均值；用均方值表示随机变量平方的平均值；用标准差表示随机变量偏离数学用标准差表示随机变量偏离数学(shxu)(shxu)期望的程度；用概率密度期望的程度

11、；用概率密度函数或概率分布函数表示随机变量在不同范围取值的概率等等，从函数或概率分布函数表示随机变量在不同范围取值的概率等等，从而全面描述随机振动的激励和响应。而全面描述随机振动的激励和响应。第8页/共81页第九页，共81页。静力法静力法一、弹性静力法一、弹性静力法日本学者大森房吉在日本学者大森房吉在1899 1899 年提出年提出(t ch)(t ch)。它假设结构各。它假设结构各个部分与地震动具有相同的振动，结构因地震作用引个部分与地震动具有相同的振动，结构因地震作用引起的惯性力起的惯性力地震力就等于地面运动加速度与结构总地震力就等于地面运动加速度与结构总质量的乘积；再将地震力视为静力作用

12、在结构上，进质量的乘积；再将地震力视为静力作用在结构上，进行结构线弹性静力分析计算。地震力计算公式：行结构线弹性静力分析计算。地震力计算公式：-地震地震(dzhn)(dzhn)系数系数将将F F作为静荷载，按静力计算作为静荷载，按静力计算(j sun)(j sun)方法计算方法计算(j sun)(j sun)结构结构的地震效应。的地震效应。第9页/共81页第十页，共81页。二、二、二、二、Pushover Pushover 法法法法塑性塑性塑性塑性(sxng)(sxng)倒塌模态分析方法（了解）倒塌模态分析方法（了解）倒塌模态分析方法（了解）倒塌模态分析方法（了解）n nPushover Pu

13、shover 法早在法早在20 20 世纪世纪(shj)60(shj)60 年代末就已经提出，在近年里得到年代末就已经提出，在近年里得到很大的发展和应用。从严格意义上看，很大的发展和应用。从严格意义上看，Pushover Pushover 法不能算作一种结法不能算作一种结构地震反应分析的方法，但它提供了一个评估结构地震反应、尤其构地震反应分析的方法，但它提供了一个评估结构地震反应、尤其是非线性地震反应的简单而有效的方法。是非线性地震反应的简单而有效的方法。n nPushover Pushover 法能够追踪结构从屈服直到极限状态的整个非弹性变形过法能够追踪结构从屈服直到极限状态的整个非弹性变形

14、过程。实际进行的程。实际进行的PushoverPushover分析过程，是一种纯粹的非线性静力分析分析过程，是一种纯粹的非线性静力分析过程，因此它与通常的非线性静力分析在计算方法上没有什么不同。过程，因此它与通常的非线性静力分析在计算方法上没有什么不同。第10页/共81页第十一页，共81页。Pushover Pushover 分析与常规非线性静力分析的主要分析与常规非线性静力分析的主要分析与常规非线性静力分析的主要分析与常规非线性静力分析的主要(zh(zh yo)yo)差别：差别：差别：差别：（1 1）Pushover Pushover 分析需要预先假定一个荷载分布模式，而常规的分析需要预先假

15、定一个荷载分布模式，而常规的非线性静力分析外加荷载是确定的。非线性静力分析外加荷载是确定的。（2 2）Pushover Pushover 分析需要预先确定与结构性能目标相对应的位移分析需要预先确定与结构性能目标相对应的位移限值，如屈服位移、倒塌破坏极限位移等，而常规的非线性限值，如屈服位移、倒塌破坏极限位移等，而常规的非线性静力分析无此要求静力分析无此要求(yoqi)(yoqi)。（3 3）Pushover Pushover 分析最终得到一条分析最终得到一条Pushover Pushover 曲线。对桥梁结构，曲线。对桥梁结构，该曲线通常为墩底剪力与上部结构质量中心处的位移之间的该曲线通常为墩

16、底剪力与上部结构质量中心处的位移之间的关系曲线，称为能力曲线；分析过程通常还计算总的结构能关系曲线，称为能力曲线；分析过程通常还计算总的结构能量耗散及等效弹性刚度，并利用单振型反应谱法计算力效应量耗散及等效弹性刚度，并利用单振型反应谱法计算力效应和位移效应和位移效应即所谓需求分析，常规的非线性静力分析则无即所谓需求分析，常规的非线性静力分析则无此过程。此过程。（4 4）Pushover Pushover 分析进行需求分析进行需求/能力比计算，以评估结构的抗震能力比计算，以评估结构的抗震性能，常规的非线性静力分析无此过程。性能，常规的非线性静力分析无此过程。第11页/共81页第十二页，共81页。

17、Pushover Pushover 法的计算法的计算法的计算法的计算(j sun)(j sun)步骤如下：步骤如下：步骤如下：步骤如下：假定一个适当的、沿高度分布的侧向荷载模式；假定一个适当的、沿高度分布的侧向荷载模式；按荷载增量法进行结构非线性分析按荷载增量法进行结构非线性分析(fnx)(fnx)，直至结构到达最终位移限，直至结构到达最终位移限值。增量形式的非线性平衡方程可以写成：值。增量形式的非线性平衡方程可以写成：第12页/共81页第十三页，共81页。计算等效单自由度系统的等效刚度和等效粘滞阻尼比；计算等效单自由度系统的等效刚度和等效粘滞阻尼比；利用反应谱方法计算结构利用反应谱方法计算结

18、构(jigu)(jigu)特征力效应和特征位移效应特征力效应和特征位移效应需求分析；需求分析；进行需求进行需求/能力比计算，评估结构能力比计算，评估结构(jigu)(jigu)的抗震性能。的抗震性能。第13页/共81页第十四页，共81页。反应谱法反应谱法反应谱法反应谱法(同时考虑地面同时考虑地面同时考虑地面同时考虑地面(dmin)(dmin)运动和结构的动力特性运动和结构的动力特性运动和结构的动力特性运动和结构的动力特性)一、单振型反应谱法一、单振型反应谱法对可以近似视为单自由度体系的结构，在计算对可以近似视为单自由度体系的结构，在计算(j sun)(j sun)出结构的出结构的振动周期之后，

19、其最大地震惯性力就可以利用规范反应谱曲振动周期之后，其最大地震惯性力就可以利用规范反应谱曲线求出。线求出。1940年，美国(mi u)皮奥特提出。第14页/共81页第十五页，共81页。第15页/共81页第十六页，共81页。第16页/共81页第十七页，共81页。回顾回顾(hug)(hug)：弹性体系的地震反：弹性体系的地震反应分析（引出地震反应谱概念）应分析（引出地震反应谱概念）一、地震作用下单自由度体系一、地震作用下单自由度体系(tx)(tx)的运动方程的运动方程质点质点(zhdin)(zhdin)位移位移质点加速度质点加速度惯性力惯性力弹性恢复力弹性恢复力阻尼力阻尼力运动方程运动方程第17页

20、/共81页第十八页，共81页。n n阻尼（阻尼（dampingdamping）：使振动能量随时间或距离逐步耗散）：使振动能量随时间或距离逐步耗散(ho sn)/(ho sn)/损损的因素，如振动系统内部质点间相对运动的阻碍、外部介质摩擦等。的因素，如振动系统内部质点间相对运动的阻碍、外部介质摩擦等。n n换句话说，换句话说，阻尼指任何振动系统在振动中，由于外界作用和阻尼指任何振动系统在振动中，由于外界作用和/或系或系统本身固有的原因引起的振动幅度逐渐下降的特性，以及此一特性统本身固有的原因引起的振动幅度逐渐下降的特性，以及此一特性的量化表征。的量化表征。第18页/共81页第十九页，共81页。二

21、、单自由度体系二、单自由度体系(tx)(tx)动力学分析动力学分析1.1.单自由度体系单自由度体系(tx)(tx)自由振动自由振动（1 1）无阻尼）无阻尼(zn)(zn)时时时时（2 2）有阻尼时）有阻尼时阻尼比阻尼比无阻尼圆频率无阻尼圆频率第19页/共81页第二十页，共81页。m m 将荷载看成是连续作用的一系列冲量将荷载看成是连续作用的一系列冲量(chngling)(chngling)，求出每个冲量，求出每个冲量(chngling)(chngling)引引起的位移后将这些位移相加即起的位移后将这些位移相加即为动荷载引起的位移。为动荷载引起的位移。2.2.单自由度体系单自由度体系(tx)(t

22、x)受迫振动受迫振动-冲量法，杜哈美积分冲量法，杜哈美积分(jfn)(jfn)求解求解一个恒力的冲量指的是这个力与其作用时间的乘积。第20页/共81页第二十一页，共81页。杜哈美积分杜哈美积分杜哈美积分杜哈美积分(jfn)(jfn)n n在振动学中著名的杜哈美积分在振动学中著名的杜哈美积分(Duhamels Integral)(Duhamels Integral)：对于受迫振动：对于受迫振动（强迫振动），可以将强迫力时程分解为一系列的脉冲的叠加，如（强迫振动），可以将强迫力时程分解为一系列的脉冲的叠加，如果已知系统在单个脉冲下的响应，并注意到果已知系统在单个脉冲下的响应，并注意到 s s 时刻

23、时刻(shk)(shk)的脉冲的脉冲只对时间只对时间 t s t s 的响应有影响，那么整个系统在的响应有影响，那么整个系统在 t t 时刻时刻(shk)(shk)的响的响应就等于所有应就等于所有t t时刻时刻(shk)(shk)以前的脉冲各自单独作用下的叠加。以前的脉冲各自单独作用下的叠加。n n用用 h(u)h(u)表示系统在单位脉冲作用下表示系统在单位脉冲作用下 u u 时刻时刻(shk)(shk)的响应。的响应。s s 时刻时刻(shk)(shk)的脉冲在系统的脉冲在系统 t(t s)t(t s)时刻时刻(shk)(shk)产生的影响就等于产生的影响就等于 h(t-s)h(t-s)，将

24、所有，将所有 s(0t)s(0t)加起来，得到整个系统在加起来，得到整个系统在 t t 时刻时刻(shk)(shk)的响应。的响应。n n对于离散时间，就是相加；对于连续时间，变成积分。对于离散时间，就是相加；对于连续时间，变成积分。m第21页/共81页第二十二页，共81页。动荷载的位移反应动荷载的位移反应(fnyng),(fnyng),不计不计阻尼阻尼m-杜哈美积分杜哈美积分(jfn)(jfn)计阻尼计阻尼(zn)(zn)时时第22页/共81页第二十三页，共81页。三、单自由度体系地震作用三、单自由度体系地震作用(zuyng)(zuyng)分分析析运动运动(yndng)(yndng)方程方程

25、或或其中其中(qzhng)(qzhng)由由Duhamel积分可得零初始条件下质点相对于地面的位移为积分可得零初始条件下质点相对于地面的位移为最大位最大位移反应移反应第23页/共81页第二十四页，共81页。质点相对于地面质点相对于地面(dmin)(dmin)的速度为的速度为质点相对于地面的最大速度质点相对于地面的最大速度(sd)(sd)反应为反应为第24页/共81页第二十五页，共81页。质点质点(zhdin)(zhdin)的绝对加速度为的绝对加速度为质点相对质点相对(xingdu)(xingdu)于地面的最大加速度反应为于地面的最大加速度反应为第25页/共81页第二十六页，共81页。地震地震(

26、dzhn)(dzhn)反应谱反应谱最大相对速度最大相对速度(xin du s d)(xin du s d)最大加速度最大加速度最大反应(fnyng)之间的关系在阻尼比、地面运动确定后，最大反应只是在阻尼比、地面运动确定后，最大反应只是结构周期结构周期的函数。的函数。单自由度体系在给定的地震作用下某个最大反应与体系单自由度体系在给定的地震作用下某个最大反应与体系自振周期的关系曲线称为该反应的地震反应谱。自振周期的关系曲线称为该反应的地震反应谱。最大相对位移最大相对位移第26页/共81页第二十七页，共81页。地震地震(dzhn)(dzhn)反应谱反应谱第27页/共81页第二十八页，共81页。位移位

27、移(wiy)(wiy)反应谱反应谱Elcentro 1940（N-S）地震记录第28页/共81页第二十九页，共81页。相对速度相对速度(xin du s d)(xin du s d)反应谱反应谱Elcentro 1940（N-S）地震记录第29页/共81页第三十页，共81页。绝对绝对(judu)(judu)加速度反应谱加速度反应谱Elcentro 1940（N-S）地震记录第30页/共81页第三十一页，共81页。相对位移反应谱相对位移反应谱,记为记为SDSD绝对加速度反应谱，绝对加速度反应谱，记为记为PSA相对速度反应谱，相对速度反应谱，记为记为PSV地震地震(dzhn)(dzhn)反应谱的特

28、点反应谱的特点1.1.阻尼比对反应阻尼比对反应(fnyng)(fnyng)谱影响很大。谱影响很大。2.2.对于加速度反应谱，当结构周期小对于加速度反应谱，当结构周期小 于某个值时幅值随周期急剧于某个值时幅值随周期急剧(jj)(jj)增大，增大，大于某个值时，快速下降。大于某个值时，快速下降。3.3.对于速度反应谱，当结构周期小于某对于速度反应谱，当结构周期小于某 个值时幅值随周期增大，随后趋于常数。个值时幅值随周期增大，随后趋于常数。4.4.对于位移反应谱，幅值随周期增大。对于位移反应谱，幅值随周期增大。gal,gal,加速度单位，伽利加速度单位，伽利略，略，1gal=1cm/s2,1gal=

29、1cm/s2,1g=980gal,g1g=980gal,g为重力加为重力加速度速度第31页/共81页第三十二页，共81页。不同场地条件对反应不同场地条件对反应(fnyng)(fnyng)谱的影响谱的影响 由于由于(yuy)(yuy)诸多随机因素影响，由不同记录得到的加速度反应谱诸多随机因素影响，由不同记录得到的加速度反应谱具有很大的随机性。在大量地震加速度记录输入后绘制得到众具有很大的随机性。在大量地震加速度记录输入后绘制得到众多反应谱曲线的基础上，再经过平均与光滑化，得到供设计使多反应谱曲线的基础上，再经过平均与光滑化，得到供设计使用的规范反应谱曲线，即平均加速度反应谱。用的规范反应谱曲线，

30、即平均加速度反应谱。周期（周期（s)s)岩石岩石坚硬场地坚硬场地厚的无粘性土层厚的无粘性土层软土层软土层结构的阻尼比和场地条件对反应结构的阻尼比和场地条件对反应(fnyng)(fnyng)谱有很大影响。谱有很大影响。第32页/共81页第三十三页，共81页。-重力荷载代表值重力荷载代表值-地震系数（反映震级、震中距、地基等的影响）地震系数（反映震级、震中距、地基等的影响）-动力系数动力系数(反映结构的特性反映结构的特性,如周期、阻尼等的影响如周期、阻尼等的影响)按静力计算方法计算结构的地震按静力计算方法计算结构的地震(dzhn)(dzhn)效应效应单振型反应单振型反应(fnyng)(fnyng)

31、谱法的地震力计算公式（谱法的地震力计算公式（19891989规范）：规范）：第33页/共81页第三十四页，共81页。反应反应(fnyng)谱常表示为动力放大系数形式，动力系数谱常表示为动力放大系数形式，动力系数 查表得查表得到到19891989规范规范(gufn)(gufn)第34页/共81页第三十五页，共81页。第35页/共81页第三十六页，共81页。20082008规范：水平规范：水平规范：水平规范：水平(shu(shu png)png)设计加速度反应谱设计加速度反应谱设计加速度反应谱设计加速度反应谱第36页/共81页第三十七页，共81页。第37页/共81页第三十八页，共81页。第38页/

32、共81页第三十九页，共81页。第39页/共81页第四十页，共81页。第40页/共81页第四十一页，共81页。第41页/共81页第四十二页，共81页。竖向设计竖向设计竖向设计竖向设计(shj)(shj)加速度反应谱加速度反应谱加速度反应谱加速度反应谱第42页/共81页第四十三页，共81页。规则规则规则规则(guz)(guz)桥梁与非规则桥梁与非规则桥梁与非规则桥梁与非规则(guz)(guz)桥梁定义桥梁定义桥梁定义桥梁定义第43页/共81页第四十四页，共81页。如柱式墩地震如柱式墩地震如柱式墩地震如柱式墩地震(dzhn)(dzhn)力计算力计算力计算力计算第44页/共81页第四十五页，共81页。

33、第45页/共81页第四十六页，共81页。二、振型分解反应谱法二、振型分解反应谱法二、振型分解反应谱法二、振型分解反应谱法(多自由度弹性体系多自由度弹性体系多自由度弹性体系多自由度弹性体系(t(t x)x)的地震的地震的地震的地震反应分析反应分析反应分析反应分析)n n对不能简化为单自由度系统的复杂结构，可通过振型对不能简化为单自由度系统的复杂结构，可通过振型分解法，化为类似于单自由度系统进行计算。分解法，化为类似于单自由度系统进行计算。n n振型（模态或模态振型）振型（模态或模态振型）:n n 结构系统按其某一自振周期振动时的变形模式。结构系统按其某一自振周期振动时的变形模式。n n振型是指体

34、系的一种固有的特性。它与固有频率相对振型是指体系的一种固有的特性。它与固有频率相对(xingdu)(xingdu)应，即为对应固有频率体系自身振动的形应，即为对应固有频率体系自身振动的形态。每一阶固有频率都对应一种振型。态。每一阶固有频率都对应一种振型。第46页/共81页第四十七页，共81页。譬如一个垂直弹簧，上端固定，下端连接一重物，组成一个简单系统。这个重物可以上下振动，可以摇摆振动，还可以绕弹簧轴扭动，这就是譬如一个垂直弹簧，上端固定，下端连接一重物，组成一个简单系统。这个重物可以上下振动，可以摇摆振动，还可以绕弹簧轴扭动，这就是该系统的三个基本振型，对应这三个振型有三个固有频率。还可能

35、有其他的振型和对应频率，通过模态分析可以计算出来。该系统的三个基本振型，对应这三个振型有三个固有频率。还可能有其他的振型和对应频率，通过模态分析可以计算出来。模态分析：指的是对一个系统（也可以是一个部件）振动特性所作的分析，以得到该系统的一系列的振型和对应各个振型的频率，所有的这些模态分析：指的是对一个系统（也可以是一个部件）振动特性所作的分析，以得到该系统的一系列的振型和对应各个振型的频率，所有的这些(zhxi)(zhxi)频率都叫做固有频率。频率都叫做固有频率。第47页/共81页第四十八页，共81页。n振型与体系实际的振动(zhndng)形态不一定相同。nn通过模态分析计算出振型和所对应的

36、频率后，按照频率从低到高的排列，来说第一振型，第二振型等等。实际结构的振动(zhndng)形态并不是一个规则的形状，而是各阶振型相叠加的结果。第48页/共81页第四十九页，共81页。n n计算系统的振型和所对应的频率意义：计算系统的振型和所对应的频率意义：n n 一个系统本身是不会振起来，但如果有外部的激一个系统本身是不会振起来，但如果有外部的激励（周期性的力，位移，速度和加速度），系统就励（周期性的力，位移，速度和加速度），系统就会发生振动，这就是会发生振动，这就是(jish)(jish)所谓的强迫振动。所谓的强迫振动。n n 通过振动理论分析可知，如果一个外部激励的方通过振动理论分析可知，

37、如果一个外部激励的方式和频率正好和该系统的某个振型及其所对应的频式和频率正好和该系统的某个振型及其所对应的频率一致，那就会发生共振现象。一旦发生了共振现率一致，那就会发生共振现象。一旦发生了共振现象，那么该系统就会发生非常激烈的振动。在没有象，那么该系统就会发生非常激烈的振动。在没有阻尼的情况下，发生共振时该系统的振幅会是无限阻尼的情况下，发生共振时该系统的振幅会是无限大。而一般情况下系统多少都会有些阻尼，包括材大。而一般情况下系统多少都会有些阻尼，包括材料的内部阻尼，此时系统的振幅不会无限大，但也料的内部阻尼，此时系统的振幅不会无限大，但也是非常大，这是不允许的，因为这个系统可能很快是非常大

38、，这是不允许的，因为这个系统可能很快会破坏，或者发出极大的噪音。计算系统的振型和会破坏，或者发出极大的噪音。计算系统的振型和所对应的频率，可以大致评估一下该系统的抗振刚所对应的频率，可以大致评估一下该系统的抗振刚度。度。第49页/共81页第五十页，共81页。n n理论上讲，有限元系统有多少个自由度，就有多少个振型。但通常忽略高阶理论上讲，有限元系统有多少个自由度，就有多少个振型。但通常忽略高阶振型（振型（3000 Hz-6000 Hz3000 Hz-6000 Hz，视具体情况而定）。因为高阶振型在有系统阻尼作，视具体情况而定）。因为高阶振型在有系统阻尼作用下，影响不大。模态分析时，给一个频率上

39、限用下，影响不大。模态分析时，给一个频率上限(或者所要计算或者所要计算(j sun)(j sun)频率频率的个数的个数)。大于该频率上限的振型，不用计算。大于该频率上限的振型，不用计算(j sun)(j sun)。n n在给定的频率范围内，振型越少越好。振型少，发生共振的可能性小。一阶在给定的频率范围内，振型越少越好。振型少，发生共振的可能性小。一阶频率越高越好，这也是为了避免共振。频率越高越好，这也是为了避免共振。第50页/共81页第五十一页，共81页。前八阶模态分别为纵向、前八阶模态分别为纵向、竖向、竖向、扭转、扭转、竖向、竖向、扭转、扭转、竖向、侧向加扭转以及竖向、侧向加扭转以及竖向方向

40、竖向方向(fngxing)的振动。大跨斜的振动。大跨斜拉桥作为一种具有较高柔拉桥作为一种具有较高柔性的结构形式，模态分布性的结构形式，模态分布较为密集，如频率小于较为密集，如频率小于1的模态共有的模态共有17阶。阶。第51页/共81页第五十二页，共81页。n在一般的有限元分析中，由于系统的自由度很多，同时在研究系统的响应时，往往只需要了解少数较低的特征值及相应的特征向量，因此在有限元分析中，发展了一些适应上述特点而效率较高的解法（子空间迭代法）。n midas Civil中除了(ch le)提供精确的特征向量法分析外，还提供了与荷载相关的Ritz向量分析法。多重Ritz向量能用于线性和非线性结

41、构的动力分析。与精确特征向量法相比，多重Ritz向量法用更少的时间可产生更精确的结果。第52页/共81页第五十三页，共81页。第53页/共81页第五十四页，共81页。1.子空间子空间(kngjin)迭代法迭代法n子空间迭代法是求解大型矩阵特征值问题的最常用最有效的方法之一，它适合于求解部分特征值解，被广泛应用于结构动力学的有限元分析中。n 子空间迭代法是假设r个起始向量（采用移频法，通过特征值的移动和已收敛的特征向量的移出，使r保持在较小的数值，从而显著提高计算效率和改进收敛速度(sd)），同时进行迭代（通过求解减缩广义特征值问题）以求得矩阵的前p（r）个特征值和特征向量。（如果r不是足够大，

42、一方面可能漏掉可能激起的振型；另一方面又可能引入不可能激起的振型）。第54页/共81页第五十五页，共81页。2.Lanczos方法方法(fngf)nLanczos方法和Ritz向量法的共同特点是直接生成一组Lanczos向量或者Ritz向量，对运动方程进行缩减，然后求解缩减了的运动方程的特征值问题，避免了迭代步骤（采用直接叠加法），从而具有更高的计算效率。n Lanczos法和Ritz向量法本质上一致，但是在实际计算中，由于计算机的截断误差和舍入误差，导致数值上的不稳定性（例如虚假的多重特征值现象），因此妨碍了Lanczos方法的实际应用。n20世纪70年代以后，很多研究工作者提出了不少Lan

43、czos向量的重正交技术(jsh)以调高其算法的稳定性，Ritz向量法从这个意义上说可以是这种，但由于他改变了成Lanczos向量的算法公式，导致以后求解的不是对角矩阵的特征值问题，而是一般矩阵特征值问题。第55页/共81页第五十六页，共81页。3.多重多重Ritz向量法（求解的是一般向量法（求解的是一般(ybn)矩阵特征值问题）矩阵特征值问题）n多重Ritz向量法认可结构动态响应是空间荷载分布的函数，考虑(kol)动力荷载的空间分布（当定义了初始向量后，第一个向量块的静态响应就来源于该初始荷载向量），可以避免漏掉可能激起的振型和引入不可能激起的振型，能够显著提高计算效率。n在计算结构动力响应

44、时，用相同数目的振型进行叠加，Ritz向量直接叠加法可以有比子空间迭代法更高的精度，这是由于后者包含了实际不被激起的振型。n子空间迭代法求出结构的前r阶振型，而Ritz向量直接叠加法求出的是和激发荷载向量直接相关的振型。n 第56页/共81页第五十七页，共81页。n因此用振型分解反应谱法和振型叠加法进行结构动力分析时，一般建议采用Ritz向量法进行结构的振型分析。程序允许在三个自由方向的加速度、静力荷载工况和内嵌非线性连接荷载中选择初始向量，可以指定任意数量的初始向量。n对于非线性连接荷载，当模型中有非线性连接单元(dnyun)时，该项自动选择为非线性连接的单元(dnyun)数量。如果没有非线

45、性连接单元(dnyun)，则该项为0。如果没有非线性连接单元(dnyun)，各非线性连接的初始荷载工况包含的初始向量数量的数值将不发生作用。第57页/共81页第五十八页，共81页。4.子空间子空间(kngjin)迭代法、迭代法、Lanczos法和多重法和多重Ritz向量法算例比较向量法算例比较第58页/共81页第五十九页，共81页。由图可以由图可以(ky)清楚的发现：清楚的发现：多重Ritz向量法和另外2种分析方法得到的振型图是大不相同的。通过模态查看，可以发现多重Ritz向量法所有振型都是对称的（荷载作用是对称的），因为它考虑了空间荷载分布状态及动力贡献，所以他忽略(hl)了所有反对称振型。

46、对于反对称振型，并不是由荷载激发的，荷载在这些振型的动力贡献为零。第59页/共81页第六十页，共81页。一一)多自由度弹性体系多自由度弹性体系(tx)(tx)动力分析回顾动力分析回顾自由振动自由振动(zhndng)(zhndng)分析分析运动运动(yndng)(yndng)方程方程设方程的特解为设方程的特解为m m1m m2-频率方程频率方程-振型方程振型方程第60页/共81页第六十一页，共81页。二二).).振型分解振型分解(fnji)(fnji)法法(不计阻尼不计阻尼)运动运动(yndng)(yndng)方程方程设设代入运动代入运动(yndng)(yndng)方程，得方程，得方程两端左乘方

47、程两端左乘第61页/共81页第六十二页，共81页。折算体系折算体系-j-j振型广义质量振型广义质量-j-j振型广义荷载振型广义荷载-j-j振型广义刚度振型广义刚度第62页/共81页第六十三页，共81页。计算计算(j sun)(j sun)步骤步骤:2.2.求广义质量求广义质量(zhling)(zhling)、广义荷载、广义荷载;3.3.求组合求组合(zh)(zh)系数系数;4.4.按下式求位移：按下式求位移：1.1.求振型、频率：求振型、频率：折算体系折算体系第63页/共81页第六十四页，共81页。三三).).振型分解振型分解(fnji)(fnji)法法(计阻尼计阻尼)阻尼力阻尼力-阻尼阻尼(

48、zn)(zn)矩阵矩阵-当质点当质点j j有单位速度有单位速度 ,其余质点速度为其余质点速度为0 0时时,质点质点i i上的阻尼力上的阻尼力.若下式成立若下式成立(chngl)(chngl)则将则将 称作正交阻尼矩阵称作正交阻尼矩阵,称作振型称作振型j j的广义阻尼系数的广义阻尼系数.第64页/共81页第六十五页，共81页。运动运动(yndng)(yndng)方程方程设设令令-第第j j振型阻尼比振型阻尼比(由试验确定由试验确定).).计算计算(j sun)(j sun)步骤步骤:1.1.求振型、频率求振型、频率(pnl);(pnl);2.2.求广义质量、广义荷载求广义质量、广义荷载;4.4.

49、求组合系数求组合系数;5.5.求位移求位移;3.3.确定振型阻尼比确定振型阻尼比;第65页/共81页第六十六页，共81页。四四).).正交阻尼正交阻尼(zn)(zn)矩阵的构成矩阵的构成其中，其中，a 0 a 0、a1 a1由试验由试验(shyn)(shyn)确确定。定。通过实测获得两个振型阻尼比通过实测获得两个振型阻尼比 和和 。同理同理-瑞利阻尼瑞利阻尼(zn)(zn)矩矩阵阵第66页/共81页第六十七页，共81页。五五)、计算、计算(j sun)(j sun)水平地震作用的振型分解反应谱法水平地震作用的振型分解反应谱法作用作用(zuyng)(zuyng)于于i i质点上的力有质点上的力有

50、m1m2mimNxixg(t)惯性力惯性力弹性恢复力弹性恢复力阻尼力阻尼力运动方程运动方程第67页/共81页第六十八页，共81页。设设代入运动代入运动(yndng)(yndng)方程，得方程，得方程两端左乘方程两端左乘第68页/共81页第六十九页，共81页。-j-j振型广义振型广义(gungy)(gungy)质量质量-j-j振型广义阻尼振型广义阻尼(zn)(zn)系数系数-j-j振型广义振型广义(gungy)(gungy)刚度刚度第69页/共81页第七十页，共81页。-j-j振型的振型参与振型的振型参与(cny)(cny)系数系数第70页/共81页第七十一页，共81页。对于对于(duy)(du