15第13章_线性动力分析_李永双

上一章我们已经谈到使用反应谱分析方法计算结构地震作用具有很多优点，但它本质上仍然是一种拟动力方法，并不是结构真实的动力响应分析，只是对于结构动力响应最大值进行估算的近似方法。但是如果仅限于线弹性范围内，反应谱分析法仍然是被认为高效而且合理的方法。

时程分析是结构更加真实的动力分析，而且随着计算机速度的提高，最初时程分析在计算时间方面的困难已经被解决。另外，最主要的是结构设计领域的不断发展，结构分析早已经超出了线弹性的范围，面对更多非线性问题，反应谱法动力分析是没有任何意义的，这时就更需要动力时程分析的发展和应用。

SAP2000为我们提供了线性时程分析和非线性时程分析功能，而且已经可以使用线性时程分析结果直接对结构构件进行设计。本章主要向大家介绍线性时程分析的内容。

13.1 线性时程分析的基本理论

线性时程分析本质上仍然是通过对结构基本动力微分方程（12.2）的求解，来得到结构在动力荷载作用下结构的基本响应大小。与反应谱分析方法不同的是，时程分析方法是将动力作用以时间函数的形式引入微分方程，并通过相应的积分方式得到结构每一个时刻的响应以及响应的变化情况。因此时程分析的积分方式是此类型分析的关键问题之一。

在SAP2000程序中，时程分析可以使用的积分方式为模态积分和直接积分。模态积分方式是以结构的模态分析结果为基础，通过结构不同模态的积分求解来得到结构总的响应值，在程序中如果选择使用模态积分求解方式，需要选择积分求解所基于的模态分析工况，因此在定义该工况前需要先定义相关模态分析工况。一般情况下，如果结构需要考虑高阶振型的影响，就应该在该模态分析工况中考虑更多的振型数目。SAP2000中所使用的模态积分方法又被称为FNA（快速非线性分析方法），我们将在后续章节非线性分析中介绍这一内容，本章主要介绍直接积分方法。

13.1.1时间积分方式

直接积分法又被称为逐步求解法，它是最常规的动力分析求解方法，其本质是在一系列时间间隔范围内求解平衡方程。在SAP2000中时间积分方法的选择是通过其它参数区域时间积分参数的定义来完成的，这里提供了Newmark法、Wilson法、排列法、Hiber-Huges-Taytor法和Chung and Hulbert法等方法可供选择，这些方法都是隐式方法，并且大部分都是无条件稳定的方法。正如Wilson教授所提出的：“只有单步的、隐式的、无条件稳定的方法才应该用于实际结构的逐步地震分析。” 由于本书的篇幅所限，本节只介绍在SAP2000中推荐使用的几种直接积分方法。

1. Newmark法

Nemark法是由Nemark在1959年提出的，并且在随后的发展过程中，许多其他研究者对其也进行了修改和改进。为了说明Newmark积分方法，考虑动力微分方程（11.2）式写成如下的形式：

t + Cut + Kut = Ft Mu直接使用泰勒级数获得下列两个额外方程：

(13.1)

t-tut=ut-t+tut2t3t-t t-t...... +uu26t2t-t ...... +u2

(13.2a)

t=ut-t+tut-tu (13.2b)

Newmark截短了这些方程并以下列形式来表达这些方程：

t-tut=ut-t+tut2t-t t3 +uu2

(13.3c)

t=ut-t+tut-t+t2 uu (13.3d)

如果假定加速度在时间步长内是线性的，可得下列方程：

ut  ut t)(u

t (13.4)

将方程(13.4)代入方程(13.3c和13.3d)，可产生标准形式的Newmark方程：

1t-t+()t2ut-t t2ut ut=ut-t+tu2

(13.5a) (13.5b)

t=ut-t(1)tut-ttut u

Newmark为结构系统中每个位移DOF的每个时间步长，通过迭代来求解方程 (13.5a、

t项。 13.5b和13.1)。通过用DOF相关的质量除以方程从方程(13.1)中获得u在SAP2000现在的版本中所使用的Newmark方法并不是上面所述，而是1962年Wilson

改进的Newmark方法。Wilson将Newmark法写成矩阵形式，并添加了刚度与质量比例阻尼，同时通过在每个时间步长对方程直接求解来消除了对迭代的需要。这要求方程 (13.5a和 13.5b) 以下列形式重写：

tb1(ututt)b2uttb3utt utb4(ututt)b5uttb6utt u其中常量b1到b6为：b1

(13.6a) (13.6b)

111b；；；b4tb1；b3222ttb51tb2；b6t(1+b3)。

将方程(13.6a和13.6b)代入方程(13.1)，可在时间“t”处按未知节点位移ut的形式写出

系统动力平衡方程。即：

ttb3utt) (b1M +b4C + K)ut = FtM(b1uttb2u

ttb6utt)  C(b4uttb5u2. Wilson法

(13.7)

2

Newmark方法本身并不是无条件稳定的，1973年Wilson在常规的Newmark方法基础上通过引入一个系数θ达到了无条件的稳定，这就是现在SAP2000程序中所使用的Wilson积分方法。 系数θ的目的是通过修改时间步长，来消除真实解附近不稳定解的振荡。如果θ等于1.0，就是Newmark方法。Wilson方法先将系数θ引入到时间步长公式：

t=t

以及一个由以下定义的荷载：

(13.8a)

Rt=Rt-t(RtRt-t)

(13.8b)

t向量之后，节其中1.0。使用Newmark方法在积分时间步长t计算了加速度u点加速度值、速度值以及位移值就由下列基本方程计算：

t=ut-tu1tut-t) (u

(13.9a) (13.9b)

t=ut-t(1)tut-ttut ut-tut=ut-t+tut2(12)t-t t2ut +u2 (13.10c)

系数θ的使用有助于在系统的高阶振型中去除数值阻尼。然而，对于高阶振型响应是很重要的问题，引入的误差可能比较大。除此之外，这一方法并不能在时间t处精确满足动力平衡方程。因此，SAP2000程序并不推荐使用Wilson方法。

3. Hiber-Huges-Taytor（HHT）法

HHT法本质上仍然是Newmark方法的发展，这一方法是Hughes在1987年提出的。HHT法将系数引入并修改结构动力方程，如公式（13.11）所示，并使用Newmark方法进行求解。

t+(1+) Cut +(1+) Kut Mu

tt+Kutt = (1+)FtFt+Cu （13.11）

该α参数取值为0到-1/3之间。对于α=0，该方法等同于α=0.5和β=0.25的Newmark

方法，也等同于平均加速度法，而且当α=0可以得到各种方法中最高的精度，但是此时会导致过多的高频振型振动，这里的高频振型具有与时间步长同阶或低阶的周期。

更多的时候可能选择α为负值，此时高频振型被阻尼严重地耗散。这不是物理阻尼，因为当使用更小的时间步长时它会减小。在非线性分析中，经常需要使用一个α负值来确保结果的收敛性。

在使用HHT方法时，为了获得最好的结果，一般使用实际的最小的时间步长，然后选择尽可能接近于零的α值。尝试使用不同的α和时间步长来确保结果不过分依赖于这些参数。

当α等于零时，HHT方法还原为常量加速度方法。它在高阶振型中产生数值能量损耗，但是它不能象用于刚度比例阻尼的阻尼比一样被预测。同时，它也不能在时间t处求解基本平衡方程。虽然存在这些不足之处，但是这一方法在计算程序的使用过程中却得到了与使用刚度比例阻尼法的效率很相近的结果。在SAP2000程序中，HHT方法也是程序默认的方法，SAP2000也建议用户使用默认的HHT方法，除非用户对其它方法有特定的需要。

3

13.1.2阻尼参数设置

当结构遭遇地震作用时，即使结构主体构件保持弹性变形状态，结构次要构件的永久变形将耗散一定的能量，从理论上讲这一部分能量是很难估计的。在结构动力分析中，这一部分能量耗散是通过阻尼来定义和实现的。对于数值计算本身，为获得稳定解，多数增量求解法也需加入一定的人工阻尼或者数值阻尼。SAP2000在时程分析中允许工程师使用几种方法定义结构在数值分析中的阻尼参数，其中涉及了模态阻尼比以及质量和刚度比例阻尼的定义方式，并且在分析过程中可以直接考虑连接单元的阻尼属性。

1. 模态阻尼比

结构反应谱分析和时程分析都需要考虑结构阻尼的影响。一般情况下，结构阻尼是通过模态阻尼比来进行定义的，也就是前面模态方程（12.7）和（12.8）中ξn项，我们定义ξn为振型n的阻尼比，它是该模态阻尼与临界阻尼的比值。

一般情况下，在结构反应谱分析和时程分析过程中，混凝土结构的模态阻尼比一般选为0.05，此时表示振动两个相邻极大值之间的衰减比为0.73，而且每个周期的应变能耗损为46.7%。钢结构的模态阻尼比一般选为0.02，此时表示振动两个相邻极大值之间的衰减比为0.88，而且每个周期的应变能耗损为22.7%。因此阻尼比设置的大小对于结构的影响是非常关键的。

振型阻尼比的值应该是在0到1之间。在过去的结构动力分析中，一般情况下各振型所采用的是相同的阻尼比，但是实测数据表明，结构高振型的阻尼比一般大于低振型的阻尼比。使用SAP2000程序时，可以为所有振型指定一个统一的阻尼比，也可以为不同振型指定不同的阻尼比，或根据周期和频率值指定对应阻尼比并默认进行插值，各振型的阻尼比之间是不相关的。

2. 质量和刚度比例阻尼

另外，在SAP2000中还经常用到的另一种阻尼定义——质量和刚度比例阻尼，它也经常用于结构非线性增量分析中。这一阻尼类型也被称为Rayleigh阻尼，它是假设阻尼矩阵与质量矩阵和刚度矩阵成正比。即：

CMK

(13.12)

在振型叠加分析中，阻尼矩阵必须具有下列属性才能分解模态方程：

TT2nnnCnnMnTnKn

T0nCm当nm

(13.13a) (13.13b)

由于质量矩阵和刚度矩阵的正交属性，该方程可重写为：

2 2nnn (13.14)

（13.14）式也可以表达为下式：

n12nn2

(13.15)

很显然，只在两个频率i与j下，就可准确指定模态阻尼，在下列方程中求解η与δ：

4

1i1i1j2j ji

(13.16)

对于一般情况，阻尼在这两个频率中被设定是相等的，因此ξi=ξj=ξ并且比例因数计算

为：

2

ij (13.17a)

ij

(13.17b)

从物理意义上讲，质量比例阻尼的假设意味着存在有外部支承的阻尼器，而使用刚度比例阻尼对结构高阶振型具有阻尼增加效应。虽然Rayleigh型阻尼没有经过物理论证，而且它的使用对大多数结构来说是难以解释的，但是使用这一阻尼方式可以用较大的时间积分步长获得稳定的数值结果，因此在SAP2000程序中仍然使用它。

3. 其它单元的阻尼考虑

目前，结构中经常使用阻尼器、隔振器等非线性连接单元，这些连接单元与一般的结构构件不同，它的目的是主动耗散结构的应变能或削弱能量传输，从力学模型来看，这类连接单元本身具备较大的阻尼值。

当模型中包含这类连接单元时，SAP2000将在动力分析过程中考虑这些阻尼的影响，并把连接单元属性中指定的线性有效阻尼系数或者非线性阻尼自动转换为振型阻尼，这一过程中将忽略振型间的交叉阻尼。这些振型阻尼值对每个振型一般是不同的，他们取决于每一振型在连接单元中引发的变形。在后面非线性连接单元定义部分，将谈到隔振、阻尼等非线性连接单元将需要同时定义线性阻尼属性和非线性阻尼属性，在线性动力分析时，转换的是线性有效阻尼；在非线性分析时，转换的是非线性阻尼。

13.2 时程曲线的输入

在SAP2000中，当需要使用时程分析法进行结构地震作用计算时，需要先将地震波引入到程序中。地震波可以是实际地震记录波，也可以是人工模拟加速度时程曲线。

图13-1 时程函数引入对话框

5

13.2.1地震时程曲线的选取

关于地震波的选取，我国2002版新规范有明确的说明：“采用时程分析法时，应按建筑场地类型和设计地震分组选用不少于二组的实际强震记录和一组人工模拟的加速度时程曲线。（《建筑抗震设计规范（GB 50011-2001）》5.1.2-3）”因此对于时程地震波的选取首先要基于建筑结构的场地类别和地震分组，另外需要选择不同的三种波进行计算，对此规范中也有相关说明。

前面我们在讲述反应谱分析中也谈到，时程分析相对于反应谱分析的一项弱点就是反应谱分析涵盖了一个地区地震作用的整体统计影响，而时程分析只能够反映一次地震作用的影响，而且即使是相同场地类型和设计地震分组的地震波也可能带来结构不同的响应结果。规范所采用的是根据小样本容量下的计算结果来估计地震效应值的方法。其理论依据是通过大量的地震加速度记录输入不同结构类型进行时程分析结果的统计分析，若选用不少于二条实际记录和一条人工模拟的加速度时程曲线作为输入，计算的平均地震效应值不小于大样本容量平均值的保证率在85%以上，而且一般也不会偏大很多。

图13-2 时程曲线定义对话框

《建筑抗震设计规范（GB 50011-2001）》中5.1.2-3还规定所选取的时程曲线“其平均地震影响系数曲线应与振型分解反应谱法所采用的地震影响系数曲线在统计意义上相符”。对于这一点规范也有所解释，“统计意义上相符”指的是时程曲线平均地震影响系数曲线与振型分解反应谱法所用的地震影响系数曲线相比，在各周期点上相差不大于20%。计算结果的平均底部剪力一般不会小于振型分解反应谱法计算结果的80%，每条地震波输入的计算结果不会小于65%。因此，《建筑抗震设计规范（GB 50011-2001）》中5.1.2-3也作了以下

6

规定：“弹性时程分析时，每条时程曲线计算所得结构底部剪力不应小于振型分解反应谱计算结果的65%，多条时程曲线计算所得结构底部剪力不应小于振型分解反应谱计算结果的80%”。在这里需要注意的是在“弹性时程分析时”，因为前面我们提到过反应谱分析仅限于线弹性范围，而使用时程分析除了可以完成线性分析外，也可以进行罕遇地震作用下的弹塑性分析。

在SAP2000程序中，联机的技术资料为工程师提供了国际上所常用的地震时程曲线和国内按场地类型不同归纳的一部分常用的地震时程曲线。对于安装单机版的用户，可以通过访问本机上的安装目录下Time History Functions文件夹来读取这些数据，如果是安装了网络版客户端的用户，可以通过访问网络安装SAP2000程序的Server下相同的目录来调用这些数据。

SAP2000所使用的时程曲线都是文本格式的，对于联机带有的地震时程曲线，可以通过文本格式对于曲线数据进行全面读取，也可以进行所需要的修改。工程师可以根据需要使用文本构建时程曲线数据，也可以使用文本编辑人工模拟加速度曲线。SAP2000读取文本格式的时程曲线时可以分别使用两种方式：一种是时间和函数值对应的数据对格式，另一种是等时间间隔的函数值关系。使用第一种时数据格式为两列，左列为时间参量，右列为对应时间参量下的加速度函数值；第二种格式所有数据均为加速度函数值，时间为默认的等间隔时间，在读取函数前需要确定等时间间隔，对于联机带有的地震波该信息已经在前面文本部分注明，并将该时间间隔值输入到SAP2000对应对话框内。

13.2.2任意动荷载时程曲线的输入

使用SAP2000可以对结构进行地震作用下的动力时程分析，但时程分析的应用范围不仅仅限于此，可以使用时程分析来完成风荷载、工业动力荷载、爆炸荷载等任意形式的动力荷载作用分析。后面在SAP2000的时程工况定义中也可以看到，在荷载类型中不仅可以选择几个方向的加速度荷载，还可以选择任意已经定义的静力荷载工况。

在使用SAP2000来模拟任意动力荷载时，必须具有该荷载的时程曲线数据。任意动力荷载也可以通过文本格式进行构建和编辑，其方法与地震时程曲线定义方法是相同的，唯一不同之处在于地震时程曲线中定义的是加速度大小，而在其它类型动力荷载中定义的是荷载值函数随时间的变化曲线。

在使用SAP2000来模拟任意动力荷载时需要注意，首先除了加速度作用模拟地震作用以外，其它类型的荷载必须是直接作用在结构对象上的荷载，比如说风荷载作用在结构迎风面上、工业冲击荷载作用在设备连接的结构结点上等，对于爆炸荷载，也是从冲击波作用在结构构件开始进行分析的，对于爆炸从爆源通过空间传播到结构上的这一过程是不进行考虑的，换句话说SAP2000并不处理波动传递和场作用问题。

13.3 线性时程分析工况的定义

定义完成时程函数曲线之后，需要定义时程分析工况。当选择定义>分析工况命令后将弹出分析工况定义对话框，然后选择添加新工况并在分析工况类型下拉菜单中选择“Time History”，也就是时程分析工况类型，便可以弹出时程分析工况定义对话框，时程分析工况全部定义都可以在这里完成。

如果需要定义的是线性时程分析，首先需要在分析类型选项中选择线性分析类型，反之选择非线性分析类型，关于非线性时程分析的相关概念和定义过程我们将在后续章节进行介绍。

7

图13-3 时程工况定义对话框

13.3.1时程类型

分析类型定义区域的右边是时程类型的定义，也就是我们前面所谈到的时程积分类型的选择，包含振型（模态）积分方法和直接积分方法两种，当选择直接积分选项时，在对话框下面的其它参数定义区域将自动弹出时间积分方式的选择，点击修改/显示按钮可以修改积分类型，其中所包含的选项就是我们前面所谈到的几种常用的积分方式，选择的积分方式将直接在其它参数区域中显示出来，默认选项是HHT积分方法，它也将在该区域显示。对于各种积分方式的基本理论前面已经进行了简述，需要再次指出的是，SAP2000程序默认积分方法是HHT法，HHT方法一般情况下会得到比较理想的结果，因此程序推荐用户使用默认的方法进行时程分析直接积分，除非用户对其它方法有特定的需要。

13.3.2时程运动类型

在线性时程分析工况下，当我们选择使用振型（模态）积分类型时，该选项下面的时程运动类型将被激活，（当选择直接积分时或选择非线性振型积分类型时，这一项将一直选择为瞬态并且是不可以修改）。时程运动类型是用来定义该时程工况运动类型的选项，此处有三个选项：瞬态、周期和静力，目前静力选项是不可用的。当时程类型定义为振型积分时，我们可以选择的运动类型可以为瞬态和周期两种。

瞬态运动类型是时程分析常用的情况，在这里结构从静止开始，时程分析定义的目的是为结构运动过程中每一个瞬态时刻指定加速度作用或所承受指定的荷载。周期运动类型是用于指定的荷载为周期性的时程分析；也就是说，这些时程作用在给定的分析长度内无限次重复，同时，时程作用效应结果也是周期性的。程序自动调节分析开始时的位移和速度，以使

8

之等于分析结束时的这些值。该结果是结构对一个任意时间变化中循环施加的荷载的稳态反应。对于线性瞬态分析，其初始状态总是零。对于一个线性周期分析，其初始状态设置为和分析结束时相同的状态。

13.3.3初始条件

在时程分析工况定义过程中，可以为本时程工况指定相关的初始条件。默认情况下，初始条件为零，也就是时程分析从结构构件零应力刚度开始。在静力分析工况中也存在初始初始条件指定功能，这些分析工况可以是考虑P-△效应或大位移效应等影响。例如，如果我们需要在水平作用分析过程中考虑由重力引起的P-△效应，可以先定义一个施加重力作用的非线性静力分析工况，选中要考虑P-△效应。然后再定义其它水平分析工况，并选择使用已经定义的非线性工况作为初始条件（或使用的刚度），这些水平分析工况就能够考虑重力P-△效应的影响。

在时程分析中，模态积分的线性时程分析工况是不能够指定初始条件的，原因是该分析使用源于模态分析工况的振型信息作为初始条件。模态积分的非线性时程分析工况和直接时间积分的线性和非线性时程分析工况都可以根据需要选择相应的初始条件。并不是所有的分析工况都可以作为初始条件指定给以上分析工况，对于不同类型的分析工况是不同的，这也是分析连续性的基本要求。只有已经定义的模态积分的非线性时程分析工况才能够作为模态积分的非线性时程分析工况的初始条件，这一对应选择关系是唯一的。对于直接积分的线性和非线性工况，可以作为初始条件的分析工况有：非线性静力工况、非线性施工顺序工况和非线性直接积分时程分析工况。

需要注意的是，初始工况对于所定义的工况的影响范围和深度是不同的，换句话说基于本工况的性质，它对于初始条件分析工况所承接的结果内容是不同的。这主要决定于本分析工况是线性分析工况还是非线性分析工况。如果本工况是线性分析工况，那么它所获得初始条件分析工况信息仅仅是刚度信息，也就是说初始工况中荷载作用下的结构效应不会被带到该工况，初始工况结构构件的内力、应力及变形效应也不会对该工况有任何影响。如果本分析工况为非线性分析工况，那么它从初始条件分析工况所获得的信息是全面的，包括刚度信息，以及全部结构效应信息（内力、应力及变形等），以及先前加载历史导致的材料非线性效应信息等。

对于复杂的土木结构，我们希望在一个单独的分析工况中去解决一个方面的全部问题，但这一点往往是不能够实现的，特别是对于施工和使用过程中经过一系列复杂过程的结构，或者结构本身需要考虑比较多的因素影响，这些因素存在共同作用且从理论上又不能够进行简单叠加。这时候就需要根据程序所提供的初始条件的指定功能来实现以上的分析，有时候可能需要指定几层初始条件来完成一个连续的分析过程。

在动力时程分析过程中，一般情况下建议使用非线性直接积分类型工况去完成连续的时程分析，这也就是我们所说的完整非线性分析。原因就是对线性时程工况指定初始分析工况时，该工况会考虑前面工况非线性刚度影响，但不会考虑后续附加的非线性影响。例如，在线性分析工况中不考虑施加的荷载引起的附加P-△效应。类似地，如果非线性工况已引起移动的单拉支撑进入受压状态，因此取消它们的刚度，则线性荷载将不能使它们重新回到受拉状态。因此在确认与非线性刚度情况相比，非线性效应很小时，可以使用线性直接积分时程分析工况。否则，一般建议进行完整的非线性分析。

13.3.4荷载作用施加

线性时程工况定义的主体部分是在施加的荷载处进行定义的，在这一区域可以定义施加荷载的类型，分为一般荷载类型和加速度荷载类型。当选择加速度荷载类型时，在后面荷载名称位置可以选择三个方向的加速度和绕三个轴的角加速度；当选择为一般荷载类型时，在后面荷载名称位置可以选择已经在静力荷载工况中定义的任意荷载工况名称。当定义完成荷

9

载类型和荷载名称的定义后，需要在函数下拉菜单选择已经定义好的函数名称（也就是前面所提到的时程曲线），以及比例系数项，此处比例系数的意义是将时程曲线中所有函数值乘以一个放大系数或折减系数，来控制时程函数整个时间范围的函数值大小。比如《建筑抗震设计规范（GB 50011-2001）》表5.1.2-2给出了结构时程分析所用地震加速度的时程曲线的最大值，这一最大值的获得就是通过时程曲线中峰值与比例系数的乘积来达到的，并且整个时间范围内函数值都会按比例进行所需的调整。

图13-4 时程工况高级荷载参数定义

上面是对于时程工况荷载定义的基本选项，如果需要对其进行更进一步的定义，可以勾选本区域左下角的显示高级荷载参数选项，将弹出荷载定义的其它选项，包括到达时间、采用坐标系和输入角度。到达时间定义了某一荷载作用波到达结构的时间，比如我们都知道地震作用包含了横波和纵波成分，纵波的传播速度要比横波快，因此结构在同一次地震作用下，横波要晚于纵波到达结构，这时就可以通过到达时间的定义来实现。到达时间一般用于加速度荷载，但是也可以用于其它荷载形式。坐标系和输入角度的定义只针对于加速度荷载，一般情况下地震加速度作用默认加载是参考全局坐标轴来完成的，在此处坐标系选项中，我们可以选择其它用户自定义坐标系统来施加地震加速度，并且可以在“输入角度”项中输入加速度与坐标轴的夹角，来模拟更为一般的地震作用输入。

13.3.5时间步数据

时程分析工况定义对话框中时间步数据区域主要是定义输出步参量的选项。虽然程序会计算每个时间步的时程函数荷载作用结构的效应，但不默认将所有结果都进行输出。最终分

10

析输出的总时间是这里所定义的输出时间步数和输出时间步长的乘积。输出时间步长需要不小于时程曲线输入的步长，否则更小的步长将没有什么意义，一般建议取输出时间步长或其整数倍。输出时间是从零时刻零步数开始算起的，当定义的输出总时间小于总的时程作用时间，程序将仅输出这一部分输出时间。

13.3.6阻尼参数

时程工况定义中另一个主要的部分是阻尼的定义，关于动力分析中阻尼设置的相关理论在前面已经进行了简述，这里只讲述程序中进行线性时程分析阻尼参数的定义方式。对于时程分析采用两种积分方式，阻尼的定义也不同，这里我们分别进行阐述。

1. 振型（模态）积分阻尼设置

当选择的是振型积分类型时，对话框中其它参数区域将显示振型阻尼选项，并显示振型阻尼定义的默认方式，点击修改/显示按钮将弹出振型阻尼定义对话框：

图13-5 振型阻尼定义对话框（常量阻尼参数定义）

此处振型阻尼类型定义有三种方式： (1)所有振型常数阻尼选项

如名称所表示的那样，使用此项是对所有振型的阻尼定义为常数。当选中该选项时，对话框显示一个所有振型常数阻尼的编辑框。例如，要施加5%的阻尼，在该编辑框中输入 0.05 即可。该对话框同时显示一个振型阻尼覆盖项区域，这一项的意义在于当对所有振型定义相同的常数阻尼后，可以选择若干振型定义阻尼覆盖项，对于定义覆盖项的振型，程序将使用覆盖项的阻尼值取代统一定义的阻尼常数值。可以使用添加、修改或删除命令对阻尼覆盖项进行相关的编辑。另外需要注意的是此项输入的阻尼值仍然是阻尼比的值。

11

图13-6 振型阻尼定义对话框（周期或频率差值阻尼参数定义）

(2)周期或频率插值阻尼选项

此项是通过定义频率或周期点系列的阻尼值，并进行插值来确定结构的振型阻尼。使用此项进行定义时，需要定义相关周期或频率值点处的阻尼值。程序在确定振型阻尼时，对于指定点处，阻尼值取定义值；对于所定义的周期或频率最大值和最小值之间的振型通过线性插值确定阻尼值；在指定范围以外，阻尼值为常数，该常数值为离指定点最近的点处的值。另外需要注意的是此项输入的阻尼值仍然是阻尼比的值。

选中该选项时，对话框显示如下选项：

按周期插值选项。如果选中该选项，则周期与阻尼编辑显示区域将变为可操作的，此时阻尼的定义是根据周期值进行插值计算的。可使用这些区域添加、修改或删除各周期的阻尼。要添加一个定义，输入一个值，点击添加按钮。要修改一个值，加亮显示要修改的定义，进行修改，然后点击修改按钮。要删除一个定义，加亮显示要删除的定义，然后点击删除按钮。

按频率插值选项。如果选中该选项，则频率与阻尼编辑显示区域将变为可操作的，此时阻尼的定义是根据频率值进行插值计算的。同样可以使用这些区域添加、修改或删除等命令对于定义进行编辑

该对话框也将显示一个振型阻尼覆盖项区域，可以选择若干振型定义阻尼覆盖项，对于定义覆盖项的振型，程序将使用覆盖项的阻尼值取代前面定义或插值计算的阻尼值。可以使用添加、修改或删除命令对阻尼覆盖项进行相关的编辑。

在使用周期插值或频率插值方法定义振型阻尼的过程中，可以随时切换定义插值定义的方式，当从一种方式转换到另一种方式时，已经定义的阻尼值将随着周期与频率的相互转换保留下来，并按照由小到大的顺序进行列表。

12

图13-7 振型阻尼定义对话框（质量和刚度振型阻尼参数定义）

(3)系数质量和刚度阻尼选项

这一项是使用由刚度矩阵乘以指定的系数和质量矩阵乘以指定的另一个系数的组合计算得到阻尼矩阵的方式。这两个系数可以直接指定，也可以通过在两个不同周期或频率处指定阻尼值等量比例值来计算，这里定义的阻尼值仍然是阻尼比的值。从数量关系上看，刚度比例化阻尼与频率成线性正比，而质量比例化阻尼与周期成线性正比。

选中该选项时，对话框显示如下选项：

直接指定选项。当选中该选项时，直接输入质量比例化系数和刚度比例化系数。所定义的比例系数的意义可以查看前面关于质量和刚度比例阻尼的理论简述。

按周期指定振型阻尼选项。当选中该选项时，在第一、第二周期和相关的阻尼编辑框中输入数值。当在这些编辑框中输入数值后，SAP2000将计算并在相应的显示栏中显示质量比例化系数和刚度比例化系数，这一参数的转换关系可以查看本章中（13.51）式。如果在阻尼编辑框中输入数值时，质量比例化系数和刚度比例化系数却没有发生变化，那么可以点击重新计算系数按钮进行计算。

按频率指定振型阻尼选项。当选中该选项时，在频率和阻尼编辑框中输入数值。当在这些编辑框中输入数值后，SAP2000将计算并在相应的显示栏中显示质量比例化系数和刚度比例化系数，这一参数的转换关系可以查看本章中（13.51）式。如果在阻尼编辑框中输入数值时，质量比例化系数和刚度比例化系数却没有发生变化，那么可以点击重新计算系数按钮进行计算。

该对话框也将显示一个振型阻尼覆盖项区域，可以选择若干振型定义阻尼覆盖项，对于定义覆盖项的振型，程序将使用覆盖项的阻尼值取代前面使用质量和刚度比例系数定义的阻尼值。可以使用添加、修改或删除命令对阻尼覆盖项进行相关的编辑。

13

另外需要注意的是与第二项周期或频率插值阻尼定义中不同，此时使用周期和频率计算指定阻尼时，第一二周期阻尼值的定义二者是完全的。质量和刚度阻尼比例系数取决于使用哪种方式及其参数的设置。

2. 直接积分阻尼设置

图13-8 直接积分阻尼定义对话框

当选择的是直接积分类型时，对话框中其它参数区域将显示阻尼选项，并显示振型阻尼定义的默认方式，点击修改/显示按钮将弹出振型阻尼定义对话框：

因为直接积分方式并不是基于模态分析，也不会进行模态分析，因此它的阻尼设置方式中不包含振型（模态）阻尼的定义方式，其阻尼的设置方式是按照质量和刚度比例阻尼直接定义的。此处质量和刚度比例阻尼的定义方式和参数意义与振型积分阻尼设置的第三种方法——质量和刚度比例阻尼系数选项设置基本上是完全相同的。唯一的不同之处是这里不包含振型阻尼覆盖项的设置，原因仍然是直接积分不会进行模态分析，因此不需要也不能够进行该项的设置。

13.4 线性时程分析结果的读取与输出

当一个线性时程分析工况被定义并运行之后，可以选择对该分析结果进行输出，时程分析的变形结果和内力、应力结果可以使用图形和文本的形式进行输出。针对于时程分析工况，可以输出时程轨迹和节点反应谱曲线。

当使用图形和文本输出时程分析的结果时，其方法与结构其它分析工况的结果输出基本相同，这在本书前面的章节已有过全面的介绍，因此在本节就不再赘述了。唯一不同的一点需要注意，当选择时程工况进行输出时，因为时程分析本质上是一系列连续时间点的全面分析，因此输出变形和内力结果时，需要我们确定是输出每个时间步（一步一步）的结果，还是包络的结果。当选择包络的结果时，程序会输出在所有时间点位置结构效应的包络值，当选择输出一步一步的结果时，程序将可以输出时程分析的每一步结果，这时可以在显示结果设置对话框中选择要输出的时间点的值，如果是图形输出，对话框下边将显示向左向右的箭头，可以通过选择向左向右箭头来连续显示时程分析其它时刻的效应结果。

针对时程分析结果的特性，程序给出了时程轨迹和节点反应谱曲线的输出，本节我们将主要介绍这一部分时程结果的输出方式和使用。

14

13.4.1时程轨迹的显示

高规4.6节给出了结构水平位移限值和舒适度的要求，其中4.6.6条规定，高度超过150米的高层建筑要满足相关条件下结构舒适度要求的限值，也就是表4.6.6结构顶点最大加速度限值的要求。这一最大加速度可以通过风洞试验来确定，也可以通过荷载规范规定的10年一遇的风荷载来计算。如果对结构进行常规的规范等效风荷载计算，是不能够输出结构节点加速度信息的。但是如果使用时程分析进行风荷载的动力分析，就可以通过时程轨迹来输出节点的加速荷载。

点击显示>显示绘图函数命令，可打开如图13-9所示的对话框。

图13-9 时程显示定义对话框

 时程工况

分析工况（多步工况） 选择一个要显示其结果的时程工况的名称。  选择绘图函数

定义绘图函数（按钮） 函数列表 垂直函数 水平函数 显示时程函数对话框。 已定义的函数。Input Energy指输入能量。 要显示其轨迹的函数。 水平坐标轴的函数。  时间范围

从„到„ 恢复默认值（按钮） 输入数值，定义一个显示结果的时间段。如果只想绘制一部分时程轨迹时，该命令比较有用。 用于将时间范围恢复为默认值。  更新轴的范围

水平

用于改变默认的水平坐标轴的范围。可控制时程轨迹的出图范围。 15

垂直 用于改变默认的竖直坐标轴的范围。可控制时程轨迹的出图范围。  轴标记

水平 垂直 定义时程轨迹水平轴名称。 定义时程轨迹竖直轴名称。  选择线对象绘图函数

实线、虚线、原点线 垂直比例系数 线颜色 定义所绘制的时程轨迹曲线线型。 定义所绘制的数值调整系数。 定义所绘制的时程轨迹曲线颜色。  其它

轴网叠加 保存命名集（按钮） 选择该选项时，时程轨迹图中将显示等间距的轴网网格。 把用户现在分析输出的设置保存为名称列表文件，它们可以在以后进行访问。 显示命名集（按钮） 显示（按钮） 完成（按钮） 点击这一按钮访问使用保存命名集命令已经定义过的输出设置文件。 点击此按钮，将显示时程轨迹。 点击此按钮，将保存对时程显示对话框的设定，并关闭时程显示定义对话框。 需要说明的是对于时程输出有时候需要设置许多参量才能够获得满意的输出效果，如果经常使用同样的数据输出方式，而每一次都重复这些设置效率是比较低的。因此程序提供可以使用保存命令集的功能，在一个或更多检查框被选中后, 把用户现在分析输出的设置保存为名称列表文件，当在以后的工程中使用相同的输出设置时，可以点击显示命令集选项，访问已经定义过的输出设置文件，结果将自动按照相应的形式进行输出。

SAP2000可以输出结构能量、构件内力、节点运动形态等一系列时间函数信息。可以通过选中构件后点击定义函数信息来定义要输出的函数信息。当选择定义绘图函数选项后，将弹出图13-10对话框。

图13-10 时程函数对话框

点击该对话框的Add Input Functions下拉列表，可选择该列表中的某一类型的函数。表13-1列出了各类型函数的含义。

表13-1 各类函数的含义 Add Load Functions 添加输入函数。将添加一个输入时程函数。 16

Add Energy Functions Add Base Functions Add Joint Disps/Forces Add Frame Forces 添加能量函数。 添加基底反力函数。 添加节点位移/力函数。 添加线对象内力函数。 在函数区域选择一个节点编号后，在下拉菜单中选择Add Joint Disps/Forces并点击修改/显示时程函数按钮，弹出如图13-11所示的时程点函数对话框。

图13-11 时程点函数对话框

 时程点函数对话框

绘图函数名称 Joint ID 向量类型 振型数 分量 定义函数的名称。 显示该点的编号，也可以填入新的节点编号， 从列出的向量类型中选择一个。 选择是包括全部振型还是只考虑一个振型。 从列出的六个向量分量中选择一个。UX、UY、UZ表示三个整体坐标上的平动，RX、RY、RZ表示绕三个整体坐标轴的旋转。 在函数区域选择一个框架编号后，在下拉菜单中选择Add Frame Forces并点击修改/显示时程函数按钮，弹出如图13-12所示的时程框架函数对话框。

图13-12 框架绘图函数对话框

17

 时程框架函数对话框

绘图函数名称 Frame ID 向量类型 振型数 位置在最近的输出站 定义函数的名称。 显示该框架的编号，也可以填入新的框架编号。 从列出的向量类型中选择一个。 选择是包括全部振型还是只考虑一个振型。 使用相对距离或绝对距离定义输出位置相对于最近测站距离。 为了结果输出对比的方便，SAP2000允许在一个图形中输出不同节点或框架某一效应的结果，并且为了在图形中能够清晰地区分，在对话框中定义每一条曲线的颜色。此时可以先点击选中垂直函数区域栏中的某个轨迹函数名称，然后在该区域为其选择某种显示线型或显示颜色。当然，对于黑白打印机而言，仅仅在颜色上区分结果曲线是不能够实现的，因此考虑到这一点，程序中可以在“选择线对象绘图函数”位置定义选择结果输出的曲线类型。

图13-13 时程函数显示对话框

 文件

打印设置 打印图形 打印表格 获取图片 打开打印设置对话框并设置打印参数。 将图形输出到打印机。 将表格输出到打印机。 获取整个表格、无标题的表格或仅图形部分，将其保存为一个图形文件。  图例

图例 列出了各颜色的线型对应的时程轨迹曲线，以及总的最大值和最小值。 13.4.2反应谱曲线的输出 SAP2000对模型中各点反应谱曲线的输出。使用显示>显示反应谱曲线，可打开反应谱生成对话框，如图13-14～13-18所示。通过该对话框，可设置要显示的反应谱参数并显示反应谱。该对话框共分为5个页面，下面分别对各页面进行了解释。

18

图13-14 反应谱生成对话框-定义页面

 时程工况

时程工况 选择要显示其分析结果的时程分析工况。  选择一个点

选择一个点 选择节点。1STORY1表示第一层的节点1。  向量方向

X方向 Y方向 Z方向 输出X方向的分析结果。 输出Y方向的分析结果。 输出Z方向的分析结果。

图13-15 反应谱生成对话框-坐标轴页面

 横坐标

频率/周期 反应谱的水平坐标轴。既可以是频率，也可以是周期。 19

 坐标

SD SV PSV SA PSA 反应谱的垂直坐标轴。SD为位移谱。 速度谱。 拟速度谱。 加速度谱。 拟加速度谱。

图13-16 反应谱生成对话框-选项页面

 横坐标

算术的 对数的 反应谱展宽 以算术方法表示反应谱的水平坐标。 以对数坐标表示反应谱的水平坐标。 以百分比表示的反应谱峰值宽度。  坐标

算术的 对数的 比例系数 以算术方法表示反应谱的垂直坐标。 以对数方法表示反应谱的垂直坐标。 对输出结果的调整系数。  轴网叠加

轴网叠加 选择该选项时，时程轨迹图中将显示等间距的轴网网格。如不选择此项，则不显示。 20

图13-17 反应谱生成对话框-频率/周期页面

 包括频率/周期

默认 结构 用户 默认频率或周期，视坐标轴页面中所作的选择而定。 由SAP2000根据结构模型计算得到。 由用户指定。在下面的包括频率/周期中输入。

图13-18 反应谱生成对话框-阻尼页面

 阻尼值

阻尼值 可指定多达6个的阻尼值，SAP2000为每个阻尼值生成一条反应谱曲线。 对上述输出参数进行设置之后，点击该对话框下面的显示按钮，弹出如图13-19所示的反应谱曲线对话框。

21

图13-19 反应谱图形显示对话框

 文件

打印设置 打印图形 打印表格 获取图片 打开打印设置对话框并设置打印参数。 将图形输出到打印机。 将表格输出到打印机。 获取整个表格、无标题的表格或仅图形部分，将其保存为一个图形文件。  图例

图例 列出了各颜色的线型对应的阻尼比。 该对话框右侧的图例中，Direction X表示当前显示的是X方向的谱，Factor 1.00指比例系数为1。Widening0.00%指谱的峰值展宽为零。该对话框右下角的坐标值是光标处于图中时，光标所在位置的坐标。图例中以6种不同颜色的文字表示图中不同颜色曲线对应的阻尼水平。

除此以外，在上面的对话框中，同样存在着保存命令集、显示命令集的定义和访问选项，其意义和使用方法与前面时程轨迹结果的输出相同，在此不再赘述。

22