1、第 4章 静力弹塑性分析原理4.1 概要4.1.1 非线性分析的目的非线性抗震分析方法可分为非线性 静力 分析方法和非线性 动力 分析方法。非线性动力分析方法可以认为是比较准确的方法，但是因为分析时间较长并对技术人员理论水准有较高的要求，所以在实际工程上的普及应用受到了限制。相反静力分析方法虽然在反映结构动力特性方面有所不足，但是因为计算效率较高和操作简单、理论概念清晰等原因被广大设计人员所普遍使用。静力弹塑性分析又被称为 Pushover分析，是 基于性能的抗震设计 (Performance-Based Seismic Design, PBSD)中最具代表性的分析方法。所谓基于性能的抗震设计
2、是以 某种目标性能 (target performance)为设计控制目标，而不是单纯的满足规范要求的极限承载能力的设计方法 。其步骤是先按照规范要求进行抗震分析和构件设计，然后通过 Pushover分析获得结构的极限承载能力，最后通过非线性位移结果评价结构是否满足目标性能要求。目前规范中推荐的基底剪力法和反应谱分析方法均为弹性分析方法，其评价标准是地震作用下的抗力不小于地震作用下产生的内力，这些方法也被称为基于荷载的设计方法。而 基于性能的设计方法则是使用与结构损伤直接相关的位移来评价结构的变形能力 (耗能能力 )，所以又被称为基于位移的设计 (displacement-based desi
3、gn)方法。通过 Pushover分析可得如图 4.1.1所示的荷载 -位移关系曲线 (能力谱 )，根据结构耗能情况可得到非线性需求谱。 能力谱与需求谱的交点就是结构对于地震作用的性能点 (performance point)。 性能点意味着结构对于地震作用所拥有的最大的非线性承载力和最大位移 ，该点在控制目标性能范围内则表示该结构满足了性能要求。通过非线性分析可以了解结构具有的的极限承载能力和安全度。5 % S m o o t h E l a s t i c D e s i g n S p e c t r aP e r f o r m a n c e P o i n tD e m a n d
4、 S p e c t r u mC a p a c i t y S p e c t r u mS p e c t r a l D i s p l a c e m e n tSpectralAccelerationSaSdDd e s i g n图4.1.1 基于位移设计法的结构抗震性能评价4.1.2 静 力 弹 塑 性 分 析 的 抗 震 设 计 原 理基于性能的抗震分析方法有下列四种。 线性静力分析法(Linear Static Procedure, LSP) 线性动力分析法(Linear Dynamic Procedure, LDP) 非线性静力分析法(Nonlinear Static P
5、rocedure, NSP) 非线性动力分析法(Nonlinear Dynamic Procedure, NDP) 其中Pushover分析方法属于非线性静力分析法，又被称为塑性铰分析法。该分析方法主要P u s h o v e rA n a l y s i sC a p a c i t y S p e c t r u m dS D O F S y s t e mrofrofC a p a c i t y C u r v eM D O F S y s t e mt r a n s f o r mFb a s eVb a s eV被应用于受高阶振型和动力特性影响较小的结构。Pushover分析就
6、是按照指定的加载模式逐渐加载至控制目标并获得结构的荷载- 位移曲线(capacity curve)，然后将其转换为单自由度体系的加速度响应和位移响应的能力谱，同时将加速度-周期格式的加速度反应谱转换为加速度-位移格式(ADRS，Acceleration-Displacement Response Spectrum)的需求谱(demand spectrum)，将需求谱和能力谱反映在同一个坐标系中，两条谱曲线的交点(性能点)就是满足该水准地震作用的极限承载能力和变形能力点。因此可通过定义不同的需求谱(小震、中震、大震)，通过验算不同性能水准下的承载力和变形，实现“小震不坏、中震可修、大震不倒”的三
7、水准(performance level)抗震设计原则。结构大师中使用了ATC-40(1996)和FEMA-273(1997)等报告中的能力谱法(Capacity Spectrum Method, CSM)和推荐的参数对构件的抗震性能进行评价。能力谱法的原理如图4.1.2所示。(a) 结构的能力曲线 (capacity curve)和能力谱 (capacity spectrum)(b) 需求谱 (demand spectrum) (c) 性能点 (performance point)图4.1.2 能力谱法(Capacity Spectrum Method, CSM)Pushover分析的目的是
8、要了解结构具有的承载能力和变形能力，钢筋砼结构在进行Pushover分析前必须先进行线弹性分析和构件设计以获得结构的配筋结果。然后才能进一步进行非线性分析。Pushover分析的优点如下： 可获得结构屈服后的响应和极限承载能力 可获得结构耗能能力和位移需求 可获得结构构件的出铰顺序 在维修加固工程中事先了解需要加固的构件4.1.3 静 力 弹 塑 性 分 析 方 法如下图4.1.3所示，结构在横向荷载作用的初期处于弹性状态，当内力超过构件的开裂或屈服内力时部分构件将发生开裂或屈服，构件和结构的刚度和阻尼都将发生变化，荷载和位移的相关关系显示非线性特性。由弹性进入屈服阶段的点A被称为弹性极限，部
9、分构件屈服后随着荷载的增加结构的位移会显著增加，到达B点后较小的外力增量也会发生较大的位移，最后在 C点后即使不再增加外力位移也会增加，C点被称为极限承载能力点。 D e m a n d S p e c t r u mn,2Tad4n,1nTdaaSR e s p o n s e S p e c t r u mt r a n s f o r m5 % E l a s t i cS p e c t r u mP e r f o r m a n c e P o i n tD e m a n d S p e c t r u mC a p a c i t y S p e c t r u mdaxAIn
10、ternalForceD i s p l a c e m e n tAE l a s t i c R a n g eCI n e l a s t i c R a n g eBP l a s t i c H i n g e图4.1.3 内力和位移的关系点C是通过荷载增量进行分析的荷载控制法所能得到稳定解的极限点，要想获得C 点之后的曲线只能通过位移增量进行分析，即采用位移控制法。结构大师中既提供荷载控制法又提供位移控制法。4.1.4 迭代分析方法Pushover分析中由于发生裂缝和屈服造成结构的刚度变化，在分析过程中会产生不平衡力也叫残余力(Residual Force)，为了消除不平衡力需要进
11、行迭代计算使不平衡力达到可以忽略的程度(满足收敛条件)。结构大师的迭代计算方法使用了完全牛顿-拉普森法(Full Newton-Raphson Method)，该方法具有收敛速度快的特点。 使用完全牛顿-拉普森法的非线性分析过程如下图所示，分析过程如下。 (1)nF(0)nK1nU(3)nU(1)n(2)4; R e f e r e n c e L o a d; D i s p l a c e m e n t; R e s i d u a l F o r c e0Pn10Pn (2)n(1)n0 (1)(2)()i1n(2)(3)n AB; I n t e r n a l F o r c e(
12、) C图4.1.4 完全牛顿-拉普森法(1). 在当前步骤(n)增加荷载向量 可得图4.1.4所示的A点，此时的平衡方程式如下。 0nP(4.1.4a)10nnKUF其中，: 当前步骤(n)的结构切线刚度矩阵nK: 当前步骤(n) 的位移增量U: 前次步骤(n-1) 的内力向量1nF: 当前步骤(n) 的加载系数: 荷载向量0P: 当前步骤(n)的荷载向量n可将式(4.1.1)用增量形式表达如下， 01nnnKUPF(4.1.4b)其中，: 当前步骤(n)的荷载增量向量 0nP解式(4.1.2)得位移增量 。nU(2). 利用位移增量 计算各单元的切线刚度和内力，将各单元的内力组合构成切线刚度
13、n矩阵 。将各单元的内力与节点力组合构成内力向量 。此时结构的内力和位移()inK ()inF的关系满足图4.1.4点B上的平衡条件。 (3). 荷载增加 时如果单元发生屈服则单元产生残余力 ，可通过下面的迭代计算0nP()inR消除残余力。(4.1.4c)() ()0()ii innniKUPFR其中，: 当前步骤(n)内的第i次迭代计算时的切线刚度矩阵 ()inK: 当前步骤(n) 内的第i次迭代计算时的位移向量iU: 当前步骤(n) 内的第i次迭代计算时的内力向量()inF: 当前步骤(n) 内的第i次迭代计算时的残余力iR解方程(4.1.3)得位移向量 。计算各单元的内力和切线刚度后可
14、得残余力 ，重()inU ()inR复(1)(3)步骤直到满足收敛条件。 (4). 满足收敛条件时(在点C)将进行下一个增量步骤的分析。 残余力和收敛计算塑形铰的出现造成了单元刚度的变化，单元刚度的变化又引起了单元内力的变化，从而使外力和单元内力之间产生了不平衡力(残余力)。程序中消除残余力的方法如下： 进行收敛迭代计算时(在时程荷载工况对话框中勾选了迭代计算选项)使用完全牛顿-拉普森法进行迭代收敛计算直至满足收敛条件。但是仍有下面的残余力累计到下一个增分步骤的外力中。当最大迭代次数输入1时等同于不进行迭代计算。a. 到最大迭代次数时仍未满足收敛条件时的残余力b. 满足了收敛条件但仍残留的不平
15、衡力 不进行收敛迭代计算时(在时程荷载工况对话框中未勾选迭代计算选项)各增量步骤的残余力将累计到下一个增量步骤中的外力中。因此即便是某个增量步骤中没有收敛只要下一个步骤中收敛时，可以认为最终分析结果收敛。 收敛判断条件 因为不可能完全消除残余力，所以为了既满足计算结果的精确度又保证计算效率，需要设置适当的收敛判断条件。迭代计算的收敛判断采用范数标准，有位移范数、荷载范数、能量范数，可选择其一也可多选作为收敛判断标准。 位移范数(4.1.4d)()()iTinDnU 荷载范数(4.1.4e)()()iTinnF 能量范数(4.1.4f)()()iTinnEU其中， : 位移范数D: 荷载范数 F
16、: 能量范数E: 当前步骤(n) 内的第i次迭代计算累计的位移增量向量()inU: 当前步骤(n) 内的第i次迭代计算的位移向量i: 当前步骤(n) 内的第i次迭代计算的累计内力增量向量()in: 当前步骤(n) 内的第i次迭代计算的内力向量iF4.1.5 初始荷载因为地震作用前结构的竖向荷载是始终存在的，所以有Pushover分析有必要考虑竖向荷载作用下的初始内力状态，这样计算的杆件的内力才是接近真实的。特别是考虑轴力和弯矩相关的柱构件在计算屈服面时需要考虑竖向荷载引起的轴力。 结构大师的Pushover分析对初始荷载也进行非线性分析，以获得更接近于实际情况的初始内力。4.1.6 加载模式P
17、ushover的 横 向 荷 载 应 该 能 相 对 准 确 地 反 映 实 际 地 震 作 用 ， 即 实 际 地 震 力 在 各 楼 层 的 惯 性 力分 布 状 态 ， 这 样 才 能 保 证 分 析 结 果 更 接 近 于 实 际 状 态 。一 般 来 说 ， 在 Pushover分 析 中 推 荐 使 用 两 种 以 上 的 横 荷 载 分 布 模 式 进 行 分 析 ， 通 过 比 较 取 不利 的 结 果 进 行 判 断 。 在 结 构 大 师 中 提 供 了 四 种 横 向 荷 载 加 载 模 式 ， 特 别 是 其 中 的 层 剪 力 模 式 使用 了 反 应 谱 分 析 得
18、 到 的 层 剪 力 分 布 模 式 ， 更 接 近 于 实 际 的 地 震 力 分 布 。 静力荷载模式：按用户定义的静力荷载分布加载。 振型模式：可按振型形状分布模式加载，也可以将几个振型线性组合。 按各楼层的质量分布比例加载 按反应谱分析的层剪力分布模式加载4.1.7 P-Delta效应结构大师的Pushover分析支持梁、柱、支撑的P-Delta效应，非线性墙单元 目前不支持P-Delta效应。P-Delta分析属于几何非线性分析，单元刚度矩阵中要考虑几何刚度，考虑P-delta效应的平衡方程如下。 10GnnKUFP(4.1.7a)其中， : 结构的弹性刚度矩阵K: 结构的几何刚度矩
19、阵G在Pushover分析过程中刚度矩阵行列式应大于零，当刚度矩阵的行列式为零或负值时将忽略几何刚度的影响。(4.1.7b)0GK刚度矩阵的行列式为零或负值的情况如下。 FEMA类型铰的位移超过极限值时(图4.1.4的C点之后)多折线类型铰在出现塑形铰以后新的单元弹性刚度矩阵加上几何刚度矩阵时，对角线上的刚度成分出现0或负值4.2 非线性分析原理4.2.1 荷载控制法和位移控制法（1）荷载控制法荷载控制法是指预设一个最大荷载后逐渐加载至最大荷载的分析方法。结构大师中的荷载控制法不是事先预设最大荷载，而是在逐渐加载后达到结构极限承载力时终止继续加载的方法，即在无法得到稳定解时终止分析的方法。 （
20、2）位移控制法位移控制法是指预设一个控制位移后逐渐加载至最大位移的方法。一般来说控制位移是指控制某个节点的位移，但是在分析过程中发生最大位移的节点可能会发生变化。结构大师中既可以将某个节点作为位移控制节点，也可以在分析过程中更新控制位移的节点，即可以控制所有节点的位移不超过某个限值。 结构大师中控制位移的方法有控制节点的最大位移、控制结构的最大位移与结构高度的比值(分为 1/120、1/100、1/50)，也可以控制结构的层间位移角。4.2.2 步长控制方法结构大师中提供下面三种步长控制方法。 自动调整步长方法 等步长方法 自定义步长函数方法（1）自动调整步长自动调整步长就是在非线性不是很明显
21、的阶段加大步长间距，在非线性比较明显的阶段自动减小步长间距的方法。具体操作步骤如下： 第一阶段：计算弹性极限( )，加载水平荷载并计算各构件的内力和构件屈服内力的比1n值，从而获得达到屈服时的荷载增量。 第二阶段：根据收敛情况自动调整分析步长( )。Lastn（2）等步长将控制荷载或控制位移按总步骤数等分作为分析步长。(a) 等步长方法 (b) 自动调整步长方法图4.2.2a 不同步长控制方法的分析结果（3）自定义步长函数方法按照用户输入的步长函数计算各步长的荷载系数。定义步长函数的方法如下：- 在Pushover荷载工况对话框中设定总步骤数。- 选择步长控制函数并输入函数值(参见下图)。其中
22、步骤数与总步骤数无关，仅是用于 定义函数的X轴的参数(X轴的最大值对应的是总步骤数，输入与总步骤数不同的数值时将按比例调整步骤数)；函数值为荷载系数值。例如当总步骤数为 10而按下面对话框中定义步长函数时，表示第2个步骤的荷载系数为0.6，第10个步骤的荷载系数为1.0。图4.2.2b 不同步长控制方法的分析结果4.2.3 子步骤功能非线性分析中为了能快速的收敛，在每个增量步骤中会采用迭代计算。每个迭代计算的过程也被称为子步骤。结构大师中的子步骤具有自动调整步长的功能。 在当前步骤计算中如果按照最大迭代次数计算没有收敛，程序将启动子步骤功能自动将荷载增量减少一半重新进行分析，如果依然没有收敛将
23、再次减少荷载增量进行分析。程序并不是所有情况都启动子步骤进行迭代计算。根据步长控制方法启动子步骤的方法如下。 (1) 自动调整步长自动调整步长的控制方法本身具有自动调整步长的功能，但是因为无限制的调整步长将严重影响计算效率，因此程序内部设置了三个等级分别对应50次、100次、200次的调整步长次数，超过这些调整次数数依然没有收敛时将启动子步骤功能第二次细分步长。(2) 等步长、步长控制函数先按固定的步长进行分析(包括迭代计算)，没有收敛时将启动子步骤功能。 4.2.4 当前刚度与初始刚度的比值结构大师中使用了当前刚度与初始刚度的比值的概念，用于评价结构目前所处的状态。结构大师中用当前刚度与初始
24、刚度的比值(Current Stiffness Ratio) 判断结构的状态的标准如下： 弹性状态 : 10.%Cs 弹性状态到极限承载能力阶段 : 0. 极限承载能力状态 : . 极限承载能力后的状态 : 0s当前步骤中Cs小于0时程序将自动回退到前一步骤并终止分析。 9 0 %8 0 %7 0 %5 0 %3 0 %0 %LoadD i s p l a c e m e n tCs=sC=1 0 0 %s( L i m i t P o i n t ):. LinearSt Cs=10blRg 16.25 0.513ypiiad25 0.67类型 B(一般的滞回特性) 25 .468ypiia
25、d类型 C(较差的滞回特性) 所有值 0.33(3) 有效周期(4.3.3c)其 中 Tn为 结 构 周 期 ， 可 取 地 震 作 用 方 向 的 第 一 振 型 ， 为 延 性 比 ， =（ api-ay） /(dpi-dy)。063.7()55ypiyieqad063.7()55ypiyieqdaR e d u c e d r e s p o n s e s p e c t r u mdSaSAVE l a s t i c r e s p o n s e s p e c t r u m ( 5 % )A c c e l e r a t i o n r a n g eV e l o c i
26、 t yr a n g e(4) 非弹性需求谱下面介绍使用有效阻尼计算非弹性需求谱的方法。使用有效阻尼系数分别计算加速度影响区段和速度影响区段的谱折减系数(spectrum reduction factor, SR)(参见图4.3.2b) 。谱折减系数可使用Newmark和Hall(1982)推荐的地面运动放大系数计算，加速度区段的谱折减系数(SR A)和速度区段的谱折减系数(SR V)的计算式见(4.3.2d)。ATC-40中对不同结构响应类型规定了谱折减系数的下限值。图4.3.2b 使用谱折减系数对弹性谱进行折减 63.73.2108ln50.342.1.6ypiyiAadforTypeA
27、SR BfC(1.4.9) 63.72.3104ln50.61.5.7ypiyiVadforTypeASR BfC表4.3.2b 结构响应类型对应的谱折减系数的下限值结构响应类型 SRA SRVType A (完整滞回特性) 1.00 0.33 0.50Type B (一般滞回特性) 0.67 0.44 0.56Type C (较差滞回特性) 0.33 0.56 0.67通过对弹性反应谱的折减可得弹塑性需求谱，该弹塑性需求谱是阻尼比为5%的弹塑性需求谱，计算性能点还需要通过计算有效阻尼对应的需求谱。下面介绍计算性能点的方法。 (4) 计算性能点的方法结 构 大 师 中 按 照 ATC-40的
28、推 荐 的 两 种 方 法 计 算 性 能 点 ， 两 种 方 法 都 是 通 过 计 算 有 效 阻 尼 获得 需 求 谱 ， 并 计 算 需 求 谱 与 能 力 谱 的 交 点 的 方 法 。Procedure-A方 法是ATC-40中提供的基本方法，其步骤如下(参见图4.3.2c)。- 首先获得能力谱的初始切线刚度直线与阻尼比为5%的弹性需求谱的交点，将该交点作为初始性能点。- 计算初始性能点上的等效阻尼以及有效阻尼，利用有效阻尼计算弹塑性需求谱，并获得弹塑性需求谱与能力谱的交点，即获得新的性能点。- 反复计算上述过程当性能点上的响应位移和响应加速度满足程序内部设置的误差范围时，将该步骤
29、的性能点作为最终性能点。图4.3.2c 使用 Procedure-A方法计算性能点的示意图(ATC-40)Procedure-B方 法计算性能点的操作步骤如下(参见图4.3.2d)：- 首先假设一个位移延性系数，使用位移延性系数计算结构的有效周期和有效阻尼，然后获得有效周期线(直线)与弹塑性需求谱(使用有效阻尼计算) 的交点。- 这样一系列的交点上的连线与能力谱的交点就是最终的性能点。图4.3.2d 使用Procedure-B方法计算性能点的示意图(ATC-40)Procedure-B方法具有收敛性能好的优点，并且只需要计算一系列不同有效阻尼比和有效周期对应的加速度谱值和位移谱值的点即可，因此
30、计算方法也相对简单。结构大师中使用两种方法计算的性能点的对话框如下。图4.3.2e 使用Procedure-A 方法计算性能点的示意图(结构大师)图4.3.2f 使用Procedure- B方法计算性能点的示意图( 结构大师)4.3.3 对结构和构件的抗震性能评价(1) 对结构抗震性能的评价目前规范提出了“小震不坏、中震可修、大震不倒”的三水准(performance level)抗震设防要求，但是没有明确给出中震可修和大震不倒的具体标准。目前建议的方法如下：- 使用Pushover分析对“小震不坏”的判断方法：根据小震的基底剪力大小找出Pushover分析的步骤，查看该步骤构件是否出现塑性铰
31、，如果没有出现则可判断为结构“小震不坏”。- 使用Pushover分析对“中震可修”的判断方法：可参考FEMA或ATC-40 上对结构生命安全的要求目标，即位移与结构总高度比不超过2%作为控制目标进行验算 (FEMA中对于立即使用要求不大于1%、防止倒塌要求不大于4%)。即计算中震性能点上的位移/总高度比值不超过2%时可判断为结构“中震可修”。- 使用Pushover分析对“大震不倒”的判断方法：当性能点上的弹塑性层间位移角小于抗规表5.5.5条的要求时可认为“大震不倒”。另外还可验算在性能点上是否满足墙柱弱梁、强剪弱弯的要求。(2) 对结构构件的抗震性能的评价对构件的性能评价除了可判断各构件
32、是否达到开裂、屈服、极限状态外，可使用FEMA-273或ATC-40中推荐的方法设置下面不同水准的要求。结构大师在后处理中提供构件在各分析步骤包括性能点上的状态。IO = 立即使用极限状态(Immediate Occupancy)LS = 生命安全极限状态(Life Safety)CP = 防止坍塌极限状态(Collapse Prevention)图4.3.3 构件的使用性能标准4.4 结构大师中的 Pushover分析操作步骤在结构大师中提供两种流程的操作步骤，一种是基于默认设置的完全自动化的操作，一种是用户按常规方法通过自定义方式的操作步骤。 自动化的操作步骤1. 运行自动生成分析数据命令
33、：可自动生成初始荷载工况、静力弹塑性荷载工况、铰类型、分配给各构件的铰特性值、分析控制条件等信息。菜单路径：主菜单 静力弹塑性分析 自动生成 自动生成分析数据。2. 检查数据：检查荷载数据和生成的铰数据的合理性。该步骤也可跳过。菜单路径：主菜单 静力弹塑性分析 检查数据 荷载数据 (或铰数据)。3. 运行分析菜单路径：主菜单 静力弹塑性分析 运行 静力弹塑性分析。4. 查看分析结果：参见使用手册的说明。 用户手动操作步骤1. 定义静力弹塑性分析荷载工况、初始荷载工况、分析控制条件等信息。菜单路径：主菜单 静力弹塑性分析 静力弹塑性荷载 荷载工况。2. 定义梁、柱、支撑的铰类型和墙纤维的本构关系。菜单路径：主菜单 静力弹塑性分析 静力弹塑性铰 定义铰特性值。3. 分别铰特性值：将铰类型分配给构件。菜单路径：主菜单 静力弹塑性分析 静力弹塑性铰 分配铰特性值。4. 检查数据：检查荷载数据和生成的铰数据的合理性。该步骤也可跳过。菜单路径：主菜单 静力弹塑性分析 检查数据 荷载数据 (或铰数据)。5. 运行分析菜单路径：主菜单 静力弹塑性分析 运行 静力弹塑性分析。6. 查看分析结果：参见使用手册的说明。
