线性结构的地震反应

地

震

工

程

学第七章线性结构的地震反应分析翟

永

梅同

济

大

学土木工程学院上

海

防

灾

救

灾

研

究

所同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention第七章

线性结构的地震反应分析‹

概述‹

动力方程的建立‹

时域分析方法‹

频域分析方法‹

振型迭加法‹

反应谱理论2同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.1概述动荷载的定义——大小、方向和作用点随时间变化；在其作用下，结构上的惯性力与外荷比不可忽视的荷载。自重、缓慢变化的荷载，其惯性力与外荷比很小，分析时仍视作静荷载。静荷只与作用位置有关，而动荷是坐标和时间的函数。动荷载的分类简谐荷载周期非简谐荷载冲击荷载确定非周期

突加荷载动荷载其他确定规律的动荷载风荷载地震荷载不确定其他无法确定变化规律的荷载同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.1概述研究内容和任务研究地震荷载作用下工程结构地震反应规律及抗震减灾理论。研究内容：第一类问题：反应分析（结构动力计算）

-----正问题输入（地震荷载）结构（系统）输出（动力反应）第二类问题：参数（或称系统）识别-----反问题输入（地震荷载）结构（系统）输出（动力反应）同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.1概述第三类问题：荷载识别。-----反问题输入（地震荷载）结构（系统）输出（动力反应）第四类问题：控制问题-----控制问题输入（地震荷载）结构（系统）输出（动力反应）控制系统（装置、能量）同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.1概述动力方程的建立动力方程的求解线性结构地震反应分析结构抗震试验方法结构动力特性及其模型化构件、结构的动力性能恢复力曲线、系统识别理论弹塑性动力分析的一般过程弹塑性结构地震反应分析多维地震波作用下的平-扭耦联系统同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.1概述建筑抗震设计隔震与减震建筑物震害的工程控制场地地震小区划地震灾害预测建筑物震害预测地下管网的震害预测同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.1概述线性结构——结构变形与外力保持线性关系，外力去除后，结构变形能完全恢复到原始状态的结构。时域分析法频域分析法确定性反应分析随机反应分析振型迭加法（结构最大反应——反应谱法）（非比例阻尼问题——复模态法）线性结构地震反应分析同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.2

动力方程的建立„

结构的离散化方法„

建立动力平衡方程的基本方法„

一维地震动输入时的动力方程„

多维地震动输入时的动力方程„

多点地震动输入时的动力方程同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.2

动力方程的建立„

结构的离散化方法实际结构质量沿结构几何形状连续分布，属于无限多自由度体系。这不仅导致分析困难，而且从工程角度也没必要。通过离散化方法将无限自由度体系转换为有限自由度体系。常用简化方法有：‡

广义坐标法（形函数法、里兹法）‡

有限单元法‡

集中质量法同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.2

动力方程的建立‡

广义坐标法（形函数法、里兹法）„

广义坐标法假定结构的位移可以用一个有限级数表示，级数的每一项都是一个随时间变化的函数和一个指定的形状函数的乘积。mυ(x)∞∑υ(x,t)

=

ϕ

(x)y

(t)

y

(t)

---广义坐标iiiϕ(x)i=1---形状函数广义坐标个数即为自由度个数inϕ

(0)

=ϕ

(l)

=

0∑(x,t)

(x)y

(t)υ

≈

ϕiiiii=1同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.2

动力方程的建立‡

广义坐标法（形函数法、里兹法）„

通过广义坐标的变换,将无限自由度体系转化成有限自由度体系。广义坐标法为结构离散化提供了基础。然而，对于杆系结构，往往不宜直接采用广义坐标法。其原因在于：1）各根杆件的形状函数不统一；2）广义坐标没有直观的物理含义，不易于由单根杆件集成为结构；3）按广义坐标法形成的系数矩阵往往不规则。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.2

动力方程的建立‡

有限单元法m和静力问题一样，可通过将实际结构离散化为有限个单元的集合，将无限自由度问题化为有限自由度来解决。结点位移个数即为自由度个数„有限单元法是以节点的位移作为结构的广义坐标，并统一规定了各杆件单元共享的形状函数，使广义坐标获得了直观的物理背景及统一的计算格式。„由于任一节点位移仅影响相临单元，结构的向量方程耦联程度较小，质量矩阵、刚度矩阵等将表现出带状特征，方便计算。„有限单元法可以认为是将体系的刚度、质量、荷载、阻尼“等效”地集中于节点处。实现这种等效的过程是通过形状函数完成各点之间的连续性要求。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.2

动力方程的建立‡

集中质量法m将实际结构的质量看成（按一定规则）集中在某些几何点上，除这些点之外物体是无质量的。这样就将无限自由度系统变成一有限自由度系统。„

集中质量法是最早提出来的离散化方法。这一方法人为地把质量集中于一些点处，与之相对应，结构的刚度特性、阻尼特性、荷载特征则被集中于质量的平移自由度方面„

不适当地集中质量可能导致较大的计算误差。应附加动能等效原则（集中质量前后体系的动能不发生显著变化）。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.2

动力方程的建立mm+αm柱m>>m梁m+

mαI

厂房排架水平振动I2I梁时的计算简图(屋架质量远大于柱子质量)三个自由度体系单自由度体系同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPreventionv(t)u(t)θ(t)三个自由度三个自由度水平振动时的计算体系构架式基础顶板简化成刚性块多自由度体系复杂体系可通过加支杆限制质量运动的办法确定体系的自由度同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.2

动力方程的建立•

结构的动力特性结构受动荷载作用，它的反应不仅和动荷载有关，而且还和结构本身固有的特性（包括结构阻尼、频率谱和振型等）有关。结构的固有特性能确定动荷下的反应程度，因此将它们称作结构的动力特性。¾自振频率和频率谱外界干扰消除后，系统在平衡位置附近所产生的振动，称作自由振动（无外荷作用的振动）。自由振动的频率称自振频率。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention•

结构的动力特性实际结构有小于等于（一般等于）自由度数的自振频率，将其按从小到达依次排列，此排列称作频率谱。频率谱中最小的频率称作基本频率，简称基频。其后依次称为第二、三频率等。它们可以通过计算和试验得到。不同结构频率谱的分布是不同的。象单跨梁、不计扭转振动的房屋等，相邻两频率间隔较大，这样的频谱称稀疏型的。对于空间结构、考虑扭转振动的房屋等，频谱中存在密集区，这样的频谱称密集型的。结构的动力反应和它的频谱有密切关系。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention•

结构的动力特性¾

结构的振型当在一定条件下结构按频谱中某一频率振动时，在任意时刻各质量的位移都保持同一比例，也即变形形状是固定的。这一变形形式称作此频率对应的振型。与基频对应的振型称第一振型或基本振型，其他依次称第二、第三振型等等。振型也可通过计算或实验得到，在多自由度体系分析时，它是重要的工具。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention•

结构的动力特性¾

结构的阻尼实际结构的自由振动都是衰减的，经一定时间后将仍处于平衡。这说明振动过程有能量耗散，这种能量耗散作用称作阻尼。产生能量耗散的原因很多，如材料的内摩擦、周围介质对能量的吸收等等。至今为止，阻尼机理仍然是没有解决的问题。为了在动力分析中考虑阻尼的影响，使分析更符合实际，人们提出了种种关于阻尼的假定。这些假定统称作阻尼理论。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention•

结构的动力特性结构分析中常用“等效粘滞”阻尼。所谓等效粘滞阻尼是假设：导致能量耗散是由于存在阻尼力，它和运动的速度成正比，方向和速度方向相反。这比例系数称阻尼系数，其数值由试验确定。阻尼系数根据这一理论，单自由度的阻尼力为速度c

y&同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.3建立动力平衡方程的基本方法„

动平衡法[M

]{u&&}+[C]{u&}+[K]{u}

={P}„

拉格朗日方程法∂T∂∂

yT∂∂

yV∂W∂

yd()

−+=

p

−cdt

∂

y&„

积分变分原理同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.2动力平衡方程的建立„

动平衡法将惯性力作为一种等效荷载，按静力平衡原理建立平衡方程。在质点运动的每一瞬时，作用在质点上的所有外力（荷载与约束力）与假想地加在质点上的惯性力互相平衡，可利用静力学的处理方法建立结构的运动方程。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention单自由度体系运动方程例-1)试建立图示结构的运动方程。P(t)

mh解：由于横梁刚度无穷大，结构只能产生水平位移。EI又设横梁（质量m）位移为u，以它为隔离

mu&&P(t)体，受力如图所示。图中F

和F

可由位移法得到Fs1

cu&

Fs2s1s212EI

kF

=

F

=

u

=

uu

us1s23h

2列x方向全部力的平衡方程，即可得结构的运动方程为hmu&&+

cu&

+

ku

=

P(t)同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.2动力平衡方程的建立„

动平衡法影响系数质量、阻尼等。aij

——由j坐标单位物理量在i坐标方向上所引起的力，可以是刚度、所有坐标j处的物理量（包括i处）与相应于坐标i处的影响系数乘积之和即为i坐标方向所受到的力值，如：n∑ffusicikij==弹性力jj=1n∑u&cij阻尼力惯性力jj=1n∑fu&&jmij=mij=1同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.2动力平衡方程的建立„

动平衡法根据动平衡法，上述各力之和等于i坐标处作用的外力pi，即n∑pm

u

+c

u

+

k

u

=(

&&

&)i=1，2…nijjijjijjij则全部n个坐标的运动方程即构成以矩阵形式表示的动力方程：[M

]{u}+[C]

{u}+[K]{u}

=

{P}&&&同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.2动力平衡方程的建立„

拉格朗日方程法将动静法与虚位移原理结合，就得到了动力学普遍方程：受有理想约束的质点系在运动过程中，其上所受的主动力和惯性力在质点系的任何虚位移上所做的虚功之和为零。在质量上考虑惯性力、阻尼力的作用，则在任意瞬时质量应该处于“动平衡”状态，因此根据虚位移原理，外力（动荷载、惯性力、阻尼力）的总虚功应恒等于总虚变形功。只要功、能函数可以用广义坐标表示，就可以由下式导出动力方程，其应用范围不仅包括线性体系，也可以包含非线性体系：∂

∂

∂d

T

T

V∂W(

)−

+

=

p−dt

∂y&

∂y

∂yc∂y同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.2动力平衡方程的建立„

拉格朗日方程法拉氏方程的特点（优点）：•是一个二阶微分方程组，方程个数与体系的自由度相同。形式简洁、结构紧凑。而且无论选取什么参数作广义坐标，方程形式不变。•方程中不出现约束反力，因而在建立体系的方程时，只需分析已知的主动力，不必考虑未知的约束反力。体系越复杂，约束条件越多，自由度越少，方程个数也越少，问题也就越简单。•拉氏方程是从能量的角度来描述动力学规律的，能量是整个物理学的基本物理量而且是标量，因此拉氏方程为把力学规律推广到其他物理学领域开辟了可能性，成为力学与其他物理学分支相联系的桥梁。拉氏方程的价值拉氏方程在理论上、方法上、形式上和应用上用高度统一的规律，描述了力学系统的动力学规律，为解决体系的动力学问题提供了统一的程序化的方法，不仅在力学范畴有重要的理论意义和实用价值，而且为研究近代物理学提供了必要的物理思想和数学技巧。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.2动力平衡方程的建立„

一维地震动输入时的动力方程定参考系和相对于定参考系作匀速直线运动的参考系都是惯性参考系，否则为非惯性参考系。参考系o’x’y’为惯性参考系，oxy为非惯性参考系，只能针对o’x’y’写出质点惯性表达式：fmu&&='i=1,2…Ni

imi相对于定参考系的总位移为=

+u

uu'iigf

(u&&m

u&&=+)则igmii同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.2动力平衡方程的建立„

一维地震动输入时的动力方程弹性力及阻尼力仅与相对位移和相对速度有关，故可在定参考系o’x’y’中应用动平衡法，得到一维地震输入时的动力方程：[]ij

j

ij

jn+

+

+∑

(

)m

u&&

u&&

c

u&

k

u=

0i=1，2…ni

i

gj=1a

u记

=

,

则gg[M

]

{u}+[C]{u}+[K]

{u}

=

−[M

]{I}a&&&g结构反应量是针对动参考系（非惯性参考系）的相对位移、相对速度、相对加速度反应。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention•

建立运动方程的基本步骤以下讨论中一律认为系统的阻尼是等效粘滞阻尼。动平衡法列方程的一般步骤为：1)确定体系的自由度——质量独立位移数；2)建立坐标系，确定未知位移（坐标正向为正）；3)根据阻尼理论确定质量所受的阻尼力；4)根据达朗泊尔原理在质量上假想作用有惯性力；5)取质量为隔离体并作受力图；6)根据动平衡原理列每一质量的瞬时动力平衡方程，此方程就是运动（微分）方程。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention•

运动方程建立举例两自由度体系运动方程Py(t)Px(t)m例)试建立图示结构的运动方程。解：沿位移正向加限制位移的支座如图所示。由位移法或弯矩分配法可做出支座单位位移的弯矩图如图示。图中1M118EI30EIM1

==MM2M2M31227l27l12EI6EIM3

=M4

=227l7lMM41同济大学Tongji

University•

运动方程建立举例上海防灾救灾研究所ShanghaiInstituteofDisasterPreventionPy(t)Px(t)两自由度体系运动方程图中m18EI30EIM1

==M227l27l12EI1k11M16EIM3

=M4

=22M7l7lM21k21由此可求得图示反力(刚度)系数kij−18EI

M348EIk12===k11k2137l37lMM4k2−18EI12EI221k12k

=22337l7l同济大学Tongji

University•

运动方程建立举例上海防灾救灾研究所ShanghaiInstituteofDisasterPreventionfIyfdy两自由度体系运动方程取质量为隔离体，加惯性力fIx、fdxf

，阻尼力f

、

f

和弹性恢复力f

、IydxdyexfIxf

Pex

xf

。由达朗泊尔原理、阻尼理论和上述feyeyPy结果可得f

=

−mu&&;

=

−

−f

c

u&

c

v&;

f

k

u

k

v=

−

−Ixdx1112ex1112f

=

−mv&&;

f

c

u

c

v;

f

k

u

k

v=

−

−=

−

−&

&Iydy2122ey2122列平衡方程并以矩阵方程表示，则得运动方程如下⎧

⎫⎡

⎤⎧

⎫

⎡c

c

⎤⎧

⎫

⎡k

k

⎤⎧

⎫

Pm

0

u&&u&ux+1112+1112=⎨

⎬⎢

⎥

⎢⎨

⎬⎥

⎢⎨

⎬

⎨

⎬⎥P0

m

v&&

c

c

v&

k

k

v⎣

⎦⎩

⎭

⎣⎦⎩

⎭

⎣⎦⎩

⎭⎩

⎭y21222122记作[k]称刚度阵同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention•

体系运动方程的一般形式在单自由度和两自由度的基础上，不难推广得到n个自由度体系的情况。在记[M]—质量阵、[C]—阻尼阵、[K]—刚度阵、[P]eq—等效荷载阵；[d]、[v]、[a]—为位移、速度、加速度阵；

[Δ]P—荷载位移阵情况下动力方程列式结果[M][a]+[C][v]+[K][d]=[P]eq同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.2动力平衡方程的建立„

多维地震动输入时的动力方程„实际地震动的地面运动，包括六个分量：三个平动分量和三个转动分量。因此，地震时结构的反应是针对六维非惯性参考系的反应。&&

+&

+=

−[

M

]({U&&

}

+

[

X

]{θ

})&&[

M

]{U

}

[C

]{U

}

[

K

]{U

}gθ&&„目前对地面转动分量的观测资料很少，大部分研究很少考虑转动的影响，而采用更为简单的动力方程：&&

+&

+=

−&&[M

]{U

}

[C

]{U

}

[K

]{U

}

[M

]{U

}g同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.2动力平衡方程的建立„

多点地震动输入时的动力方程对于长跨桥梁、坝体、管道网路等较大范围内的结构，地震波在结构基础面上的传播要经历一定的时间。即在同一时刻，结构各支撑点所承受的地面运动是不同的，这就是所谓的多点地震输入问题。这时必须考虑各支撑点间相对运动引起的结构内的拟静力应力。&&&−1&&[M

]

{U}

[C]

{U}

[K]

{U}

[M

][K]

[K

]{U

}++=ggl&&U[Kg]为因支座相对运动所产生的弹性耦合矩阵；为支座加速gl度过程，应注意，对不同的点，在同一时刻的值是不同的，具体过程由支撑输入地震波的情况决定。当仅考虑单方向水平波的输入时，不同的

可只考虑相位差的变化，可引用错时输入的办&&Ug法来处理地震反应的计算过程。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.2动力平衡方程的建立„

多点地震动输入时的动力方程当地面各点运动加速度过程完全相同，即各支撑点处的相对位移为零，则−1=

−[K]

[K

]

[I]g运动方程变为：&&

&&&[M]{U}+[C]{U}+[K]{U}=−[M]{U

}g同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention7.3

时域分析方法„

线性加速度方法„

纽马克β法和威尔逊θ法„

时域分析方法的收敛性与稳定性„

阻尼的处理同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

线性加速度方法•运动增量方程：Δ

&&

+[M

]{

U}

[C]{

U}

[K]{

U}

[M

]{

U

}Δ

&

+Δ=

−Δ

&&jjjgj[M]不随时间变化，如果已知t时刻[C]、[K]、位移、速度、加速度（称状态向量），设法从增量方程求得位移、速度、加速度的增量，则显然可以求得t+Δt时刻状态向量，重复这一过程即可求得非线性问题的数值解答。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

线性加速度方法•运动增量方程：Δ

&&

+[M

]{

U}

[C]{

U}

[K]{

U}Δ

&

+Δ=

−[M

]{

U

}Δ

&&jjjgjt

t

tΔ

=

−线性加速度算法假定，在时段内，结构加j+1j速度反应是关于时间的线性函数，即取：&&U−

U&&Δ

U&&jU&&&

=j+

1j==

常量jΔ

tΔ

t将位移{U}按泰勒级数在tj附近展开：&&&&&&{U}{U}{U}U

t

+τ

=

U

+

τ

+

τ2+

τ3+...{

(

)}

{

}jjjjj1!2!3!τ对

求导，U

t

+τ

=

U

+

U

τ

+

U

τ&&&&{

}

...&&&+2{

(

)}

{

}

{

}12jjjj同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

线性加速度方法{U(tj

+τ)}={U}当τ=∆t时，j+1&&&&&&U

−

U

=

U

Δt

+

U

Δt{

}2+

U

Δt{

}

...3+{

}

{

}

{

}1216j+1jjjj&

&

&&&&&{

}

{

}

{

}1

{

}

...U

−

U

=

U

Δt

+

U

Δt2

+j+1jjj2U&&

−U&&ΔU&&将假定+U&&&

=j

1j=j=

常量

带入上两式，则jΔtΔt&&&&&Δ

=

Δ

+{

U}

{U}

t{U}

tΔ2+

Δ

Δ{

U}

t21216(1)(2)jjjj&

&&{ΔU}

={U}

Δt

+&&{ΔU}

Δt12jjj同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

线性加速度方法&&{

U}Δ{U}Δt&&{ΔU}

=

6

−

6

−3{U}1）可得jj由（jjΔ

2t{ΔU}j将上式带入（2）可得&&&&{Δ

}

=3U−3{

}

Δ

−

{

}

ΔU

t

U

t12jΔtjj将上两式带入运动增量方程，即可将原来的增量微分方程化为关于ΔU的代数方程，并令[

]

[

]

[

]

[

]K

=M

+C

+

K6Δt23Δtj{U}

+

3{U}

)

+[C](

3{U}

+

Δt

{U}

)&&{ΔP}

=

[M

](−{ΔU

}+&&&&&&6Δtjgjjj2j同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

线性加速度方法则运动增量方程可写为（拟静力增量方程）：[

]K

{ΔU}

={ΔP}jjj[

]K{ΔP}为拟静力刚度矩阵，

为拟静力荷载向量状态数j为0,1,2,….N，N为计算时程的离散时段数。根据微分方程的初始条件，在任一时刻tj+1,&&&{U}{U}jj均可得知，故可象求解静力位移那样求解上述拟静力增量方程。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

线性加速度方法特别注意：¾拟静力刚度矩阵不仅与[K]有关，而且与质量和阻尼有关；拟静力荷载向量不仅取决于地震加速度增量，而且取决于前¾一时刻的计算反应值&&&{U}{U}jj¾这样将使动力分析计算中的误差逐级积累，严重时甚至导致结果的发散¾为了减少误差积累，可以应用加速度平衡校正算法，根据&&{ΔU}增量动平衡方程求解Δ

&&

=

−

Δ

&&

−{

U}

{

U

}

[M

]

([C]{

U}

[K]{

U}

)Δ

&

+

Δ−1jg

jjj同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

线性加速度方法•线性加速度方法一般步骤：[M

][K][C]、阻尼阵1、生成质量阵、刚度阵2、计算拟静力刚度矩阵；&&&{U}{U}3、从初始条件、开始计算拟静力荷载向量004、求解拟静力增量方程，得到相对位移增量；5、计算相对速度增量，分别迭加相对位移增量与相对速度增量得到本步末的计算位移与计算速度；6、计算本步加速度增量，迭加加速度增量得本步末加速度计算值7、以本步末的速度和加速度作为初始状态，返回到第3步继续下一时段的计算。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

线性加速度方法‹线性加速度方法特点：‡

计算和理论分析表明，为使计算有足够的精度，积分步长应小于系统周期的十分之一。‡

不仅能用于分析线性结构，而且能用于分析非线性结构。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

纽马克β法和威尔逊θ法1、纽马克β法是一种将线性加速度方法普遍化的方法。其基本方程为：1{U

}

{U

}

{U

}

t

(=+

&

Δ

+

−

β

&&

Δ+

β

&&){U}

(

t)

{U

}(

t)Δ22j+1jjjj+12&&&&&&{U

}

=

{U

}+

(1−α){U}

Δt+

α{U

}Δtj+1jjj+1显然，线性加速度法是纽马克β法在α＝1/2，β＝1/6时的特例。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

纽马克β法和威尔逊θ法2、威尔逊θ法是线性加速度法的变形。两者的区别在于，线性加速度方法是在时刻t＋Δt使用运动方程，而θ法则将运动方程应用于更后一点时刻t＋θ·Δt（θ>1）。其基本方程为：θΔ2(

t)θΔ2(

t){U

(t

+

θΔt}

=

{U

(t)}

+

{U

}θ

⋅

Δt

+&{U&&(t)}

+{U&&(t

+

θ

⋅

Δt}36θΔtθΔt{U

(t+

θΔt}

=

{U

(t)}+

{U

(t)}+

{U

(t+θ

⋅Δt}&&&&&&22而在t＋θ·Δt时的运动方程为：&&&&&[M]{U(t+θ

⋅Δt)}+[C]{U(t+θ

⋅Δt)}+[K]{U(t+θ

⋅Δt)

}=

−[M]{U

(t+θ

⋅Δt)}g同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

纽马克β法和威尔逊θ法应用上述基本公式求得t＋θ·Δt时刻的反应值，再在t与t＋θ·Δt之间作线性内差求得t＋Δt

时刻的反应量作为下一时段计算的初始状态，继续下一步计算。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

时域分析方法的收敛性与稳定性1、时域分析的收敛性，是指当步长Δt趋于无穷小时，逐步积分法解的误差也趋于无穷小。以上方法都是收敛的。2、稳定性，是指差分法在任意时间步长上所得到的解是否会因为初始条件或计算过程中舍入误差的扩散而导致无限的增长或振荡问题。若结果不受时间步长的影响，则算法是无条件稳定的，反之为有条件稳定。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

时域分析方法的收敛性与稳定性纽马克β法，当α＝1/2，β＝1/4时（平均常加速度值）为无条件稳定的。α＝1/2，β＝1/6时（线性加速度算法）为有条件稳定的（稳定步长的临界值为Δt/T

<0.55,T

为指定结构最minmin小周期）。对威尔逊θ法，采用全量形式求解，θ≥1.37时为无条件稳定的。3、具有收敛性和稳定性的算法并不能保证解的误差在可接受的范围内，控制计算误差在容许误差范围内的方法与积分步长的选取密切相关。对于线性结构，Δt的选取应满足：Δt＝min{1/6T

,1/2T

,

ΔT

}se0T

为结构最低周期,T

为地震波最低周期分量，ΔT

为地震波时se0程的数值化时间间隔同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

阻尼的处理对于均质材料或单一类型的材料，可以采用瑞利比例阻尼的假定：

[C]=a[M]+b[K]比例常数可以根据振型分解方法由选定的两个振型阻尼比和相应的自振频率表示。ζζa

=2(

−

)/(

1

−

1

)jjiω

ω22ω

ωiij2=

ζ

ω

−ζ

ω

ω

−ω2b

2()/(

)j

ji

ijiωζi、

i

为第i振型的阻尼比和自振频率。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention第三节

频域分析方法一、什么是频域分析方法二、频域传递函数三、线性单自由度体系的地震反应四、线性多自由度体系的地震反应同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

什么是频域分析方法？•线性系统的动力响应，在时域内表现为振幅反应时程随时间的变化，在频域内则表现为系统能量在各频段内的分布。•时域分析方法的基本思路是将时间过程离散化，在每个小时段内把动力问题化为拟静力问题求解，然后迭加得到总体反应。••频域分析的基本思路是将频域离散化，针对每个小频段内的动力问题运用频域传递函数求解，然后迭加得到总体反应。广泛应用于土结构相互作用的确定性地震反应分析和线性结构的随机地震反应分析等领域。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

频域传递函数„„传递函数描述系统输入与输出量之间的一般关系。频域传递函数描述线性系统输出量与输入量在频域内的传递关系。设一般的线性系统为：

ωH(i

)a

y

a

y+a2

yL

a

y

b

x

b

x(n−2)

+

+

=+L

b

x(m−1)

+

+(n)(n−1)+(m)01n01m式中，a，b表示与时间无关的常系数；y为输出量；x为输入量，（.）表示关于时间t的导数的阶数。对上式作关于t的傅立叶变换，可得一般项：∞∫(n

−

j

)

−

iω

t=

ω

ωa

(

i

)

Y

(

)n

−

ja

j

yedtj−

∞∞∫(m−

j)

−iωt

=

ω

m−

j

ωb

x

e

dt

b

(

i)

X(

)jj−∞同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

频域传递函数n∑A

(iω)

=

a

(ωi)n−

j令：njj=0m∑B

(iω)

=

b

(ωi)m−

jmjj=0A

iω

Y

ω

=

B

iω

X

ω(

)

(

)则方程化为：

(

)

(

)nmBm

(iω)令H(iω)

=An

(iω)Y(ω)

=

H(iω)X(ω)则方程变为：H(iω)称为系统的频域传递函数。对于线性系统，计算传递函数时，只要把原微分方程中的微分算子带之以

(iω)n

并做简单代数运算即可。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

频域传递函数频域传递函数表示输入为单位谐和激励时的系统稳态反应。是一个复数，H(iω)可以写成模与幅角的形式：H(i

)

H(i

)eω

=

ω

−ϕ(ω)iH(iω)式中，模表示系统反应与激励在频域内的幅值ϕ(ω)表示反应与激励之间的比，又称为增益因子；幅角相位差，又叫相位因子同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

线性单自由度体系的地震反应线性单自由度体系在一维地震动输入时的动力方程为：mu&&

+

cu&

+

ku

=

−mu&&g令ω0

=

k

/

m,ξ

=

c/

2mω0u&&

+

2ξω

u&

+

u

=

−u&&ω2则动力方程变为：00gu&&g由于输入量是地面加速度而不是地面位移，所以关于输入量的算子D(n)的阶数是零。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

线性单自由度体系的地震反应此体系的频域传递函数为：1H

(iω

)

=ω2−

ω2+

2ζω

0ω

i0−

u&&通过傅立叶变换，可以将输入转换为频域内的函数：g∞∫(ω

)=

−

u

t

e−

i

ω

tdtP&&

(

)g−

∞则线性单自由度体系的频域地震反应解为：P

(ω

)ω

2

−

ω

2

+

2ξω

ω

iu

(ω

)

=00同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

线性单自由度体系的地震反应通过傅立叶逆变换，可将解答u(w)转化为时域反应：∞∫u(t)

=

1

u(ω)e

ω

dωi

t2π−∞同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention•频域方法进行线性体系地震反应分析的一般步骤H(iω)1、根据动力方程求出频域传递函数P(ω)u(ω)2、用快速傅立叶变换将输入变换为频域内的函数3、应用频率传递函数对每一频率分量求出反应量u(t)4、采用快速傅立叶逆变换将反应量转化为时域反应同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

线性多自由度体系的地震反应多自由度体系一般是多输入、多输出系统，因此，频域传递函数有交叉性。交叉性频域传递函数计算的原则是：逐个输入，分离输出。为此，先定义广义频域传递函数的概念H

(iω

)广义频域传递函数

是指在第k个自由度处输入单位谐lk和激励时所引起的第l个自由度的输出反应值。因此，当第k个自由度处输入不是单位谐和激励，而是一般的

x

k

(

t

)

时，则第l自Y

ω

=

H

iω

X

ω由度处关于频率的反应值为：

(

)(

)

(

)lklkk设线性多自由度体系具有输入的自由度数为m，则所有输入在l端产生的反应值可由迭加原理给出：m∑Y

l

(ω

)

=H

(iω

)

X

(ω

)lkkk

=

1同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

线性多自由度体系的地震反应采用广义频域传递函数概念进行多自由度体系频域分析的优点在于它提供了一个一般的理论框架，就具体计算而言，这种分析是相当繁琐的。因此实际应用中往往利用振型分解法。即先将多自由度体系转化为一系列等效单自由度体系，然后利用单自由度体系的频域分析方法进行计算，最后应用振型迭加原理给出总体反应。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention第四节

振型迭加法一、振型迭加法的引入二、振型时域分析法三、振型频域分析法四、振型分解反应谱法同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

振型迭加法的引入‹

线性体系的时域分析法与频域分析法都是从对输入的离散化着手进行体系动力反应的分析。‹

振型迭加法则通过对结构振动特征的离散化来实现体系动力反应的离散化，然后根据实际需要，选取部分主导特征反应，运用迭加原理求取结构体系的动力反应。‹

对振动主导特征反应量的计算则可以采用时域分析的方法，也可以采用频域分析的方法。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention„

振型迭加法的引入‹

利用振型分解原理，将耦合的动力方程化为解耦的等效单自由度方程分别求解，然后将各振型反应迭加起来，获取体系的总动力反应，即为振型迭加法。‹

在结构体系的运动过程中，前几个低阶振型的运动在总运动中依次占据主导地位。对于高层建筑或动力自由度较多的体系，一般选取前9～15个振型，对于大量的较低的一般性建筑，可以取前1～3个振型分析。同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention•运动方程线性多自由度体系输入一维水平地震动的运动方程mmmmNiI

=

m

(&x&

+

&x&

）惯性力xiiigi弹性恢复力=++S

k

x

k

x

Lk

xii1

1i2

2in

n21=++阻尼力

R

c

x

c

x&

Lc

x&&ii1

1i2

2in

nnn∑

∑运动方程+c

x&+=

−k

x

m

&x&m

&x&i

iij

iij

ii

gxg(t)j=1j=1i=1,2,LN[m]{x}+[c]{x}+[k]{x}=

−[m]{I}x

(t)&&&&&g同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention•运动方程的解运动方程的解：采用结构动力学中的振型分解法，多自由度线性体系的振动位移x(t)可以表示为各振型下位移反应的叠加（线性组合）。∑x

(t

)

=

X

q

(t

)ijij3(

)1

32

33

3++2(

)1

22

23

21(t)2(t)31(

)1

12

13

1同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention•运动方程的解——续以两个自由度线性体系为例将质点

m

和

m

在地震作用下任一时刻

的位移

x

(t)和

x

(t)用1212其两个振型的线性组合

来表示，即：x

(t)

=

q

(t)X

+

q

(t)X

⎫1111221⎬x

(t)

=

q

(t)X

+

q

(t)X2112222⎭q

(t)

、q

(t)是时间的函数，称为广

义坐标，表示在质点任

一12时刻的变位中第一振型

与第二振型所占的分量

。[m]{x}+[c]{x}+[k]{x}=

−[m]{I}x

(t)&&&&&代入运动方程g同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention•运动方程的解——续假定阻尼矩阵是质量矩阻尼矩阵亦能满足正交[c

]=

α

[m

]+

α

[k

]阵和刚度矩阵的线性组条件，以消除振型之间合，从而使的耦合，即令：12式中

α

、

α

−

−

−

比例常数故得：12[m

][X

]{q}+

(α

[m

]+

α

[k

])[X

]{q}+

[k

][X

]{q}=

−[m

][x

]&&&&&120将上式等号两边各乘以{X

}

，得：Tj{X

}

[m

][X

]{q}+

{X

}

(α

[m

]+

α

[k

])[X

]{q}+

{X

}

[k

][X

]{q}=

−{X

}

[m

][x

]TjTjTjTj&&&&&120根据振型对质量矩阵和刚度矩阵的正交性，整理后得

;q&&

+

(α

+

α

ω)q&

j

+

ωq

=

−γ

&x&0(

j

=

1,

2,

L

n

)2j2jj12jjn∑m

i

Xm

i

X{X

}

[m

]{1}Tjji式中γ

j

==i

1={X

}

[m

]{X

}Tjn∑2jiji

=1γ

j

−

−

−

体系在地震反应中第即当各质点位移x

=

x

=

L

x

=

L

x

=

1时的

q

值。j

振型的振型参与系数。12jnj同济大学Tongji

University上海防灾救灾研究所ShanghaiInstituteofDisasterPrevention•运动方程的解——续则运动方程变成q&&

+

2ζ

ω

q&

+ω

q

=

−γ

&x&

(

j

=1,

2,

Ln)2jj

j

jj

jj

0式中ζ

j

−

−

−

对应于

j振型的阻尼比，ωj为相应于第j振型的原频率；γ

j为第j振型的振型参与系数