Structural Analysis for PerformanceBased Earthquake Engineering

Basic modeling concepts Nonlinear static pushover analysis Nonlinear dynamic response history analysis Incremental nonlinear analysis Probabilistic approaches

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 1

In performance-based engineering it is necessary to obtain realistic estimates of inelastic deformations in structures so that these deformations may be checked against deformation limits as established in the appropriate performance criteria. Two basic methods are available for determining these inelastic deformations: Nonlinear static pushover analysis and Nonlinear Dynamic Response History analysis. Pushover analysis is the subject of the next several slides.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 1

Nonlinear Static Pushover Analysis

Why pushover analysis? Basic overview of method Details of various steps Discussion of assumptions Improved methods

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 2

These are the basic subtopics discussed in the section on pushover analysis.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 2

Why Pushover Analysis?

Performance-based methods require

reasonable estimates of inelastic deformation or damage in structures.

Elastic Analysis is not capable of providing

this information.

Nonlinear dynamic response history analysis is

capable of providing the required information, but may be very time-consuming.
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 3

The use of pushover analysis may simply be the lesser of all evils. Elastic analysis does not have the capability to compute inelastic deformations, hence it is out. Nonlinear response history analysis (NRHA) is certainly viable but is very time consuming. Also, NHRA may produce a very wide range of responses for a system subjected to a suite of appropriately scaled ground motions. Computed deformation demands can easily range by an order of magnitude (or more) making it difficult to make engineering decisions. Hence, we are left with Nonlinear Static Pushover Analysis (NSPA) as a reasonable alternative.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 3

Why Pushover Analysis?

Nonlinear static pushover analysis may

provide reasonable estimates of location of inelastic behavior.

Pushover analysis alone is not capable of

providing estimates of maximum deformation. Additional analysis must be performed for this purpose. The fundamental issue is How Far to Push?
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 4

NSPA, in addition to providing estimates of deformation demands, provides some useful insight into the pattern of inelastic deformation that may occur. This is very important when assessing desirable behaviors such as strongcolumn weak-beam behavior. In NSPA an inelastic model is developed and is subjected to gravity load followed by a monotonically increasing static lateral load. While the load pattern is defined, the magnitude of the load is not. The fundamental question in pushover analysis is how far to push? Other computational tools, such as the Capacity Spectrum Approach must be used in concert with NSPA to determine how far to push.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 4

Why Pushover Analysis?

It is important to recognize that the purpose

of pushover analysis is not to predict the actual response of a structure to an earthquake. (It is unlikely that nonlinear dynamic analysis can predict the response.)

The minimum requirement for any method

of analysis, including pushover, is that it must be good enough for design.
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 5

It is very important to note that the purpose of NSPA is not to predict the actual performance of a structure. It is doubtful that even NRHA can do this. The purpose of NSPA is to provide information which may used to assess the adequacy of a design of a new or existing building.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 5

Basic Overview of Method

Development of Capacity Curve Prediction of Target Displacement

Capacity-Spectrum Approach (ATC 40) Simplified Approach (FEMA 273, NEHRP) Uncoupled Modal Response History Modal Pushover

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 6

A pushover analysis consists of two parts. First, the pushover or Capacity Curve is determined through application of incremental static loads to an inelastic model of the structure. Second, this curve is used with some other Demand tool to determine the target displacement. A variety of demand tools are available, four of which are presented on this slide. In this course emphasis is placed on the first two approaches.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 6

Development of the Capacity Curve (ATC 40 Approach)

1. Develop Analytical Model of Structure Including: Gravity loads Known sources of inelastic behavior P-Delta Effects 2. Compute Modal Properties: Periods and Mode Shapes Modal Participation Factors Effective Modal Mass 3. Assume Lateral Inertial Force Distribution 4. Construct Pushover Curve 5. Transform Pushover Curve to 1st Mode Capacity Curve 6. Simplify Capacity Curve (Use bilinear approximation)
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 7

The first approach covered is the so-called Demand Capacity Spectrum approach. This method is described in detail in the ATC 40 document.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 7

Development of the Capacity Curve

Pushover Curve Base Shear Capacity Curve Modal Acceleration

Roof Displacement

Modal Displacement
Advanced Analysis 15-5b - 8

Instructional Material Complementing FEMA 451, Design Examples

The nonlinear static analysis of the structure produces a pushover curve as shown at the left. The symbol above the curve indicates that for this curve the lateral load pattern was upper triangular. Other load patterns, such as uniform or proportional to first mode shape will produce different pushover curves. The curve at the right is a simplified first mode bilinear version of the pushover curve. This curve is called a Capacity Curve, or Capacity Spectrum. Note that the quantities on the X and Y axes of the capacity curve are modal acceleration and modal displacement. Details on the development of the Capacity Curve are provided later.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 8

Development of the Demand Curve

1. 2. 3. 4. 5.

Assume Seismic Hazard Level (e.g 2% in 50 years) Develop 5% Damped ELASTIC Response Spectrum Modify for Site Effects Modify for Expected Performance and Equivalent Damping Convert to Displacement-Acceleration Format

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 9

The next step in the analysis is to compute the Demand Curve. This is basically an elastic response spectrum that has been modified for expected performance and equivalent viscous damping. The modifications are HIGHLY EMPIRICAL. The various steps in the development of the demand curve are given here. Details are provided later.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 9

Elastic Spectrum Based Target Displacement

Base Shear/Weight or Pseudoacceleration (g) Elastic Spectrum based demand curve for X% equivalent viscous damping

Point on capacity curve representing X% equivalent viscous damping.

Target Displacement

Spectral Displacement
Advanced Analysis 15-5b - 10

Instructional Material Complementing FEMA 451, Design Examples

The Demand Curve is used in concert with the Capacity Curve to predict the target displacement. A trial-an-error procedure is typically used to compute the target displacement.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 10

Review of MDOF Dynamics (1)

Original Equations of Motion:

&& + Cu & + Ku = MRu &&g Mu

K = M 2

Transformation to Modal Coordinates:

1 1 R= . 1

y1 y y = 2 = [1 2 3 ... n ] . yn & + C y & + Ky = MRu &&g M& y

Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 11

u = y

Recall that a main step in the NSPA procedure is the conversion of the pushover curve (in the force vs. displacement domain) to the capacity curve (in the spectral acceleration vs. spectral displacement domain). To facilitate an explanation of this conversion, a review of MDOF dynamics is provided. Here the MDOF equations are shown. Terminology follows that in the textbook by Clough and Penzien. The first step is to transform from natural coordinates (displacements at the various DOF) to modal coordinates (amplitudes of mode shapes).

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 11

Review of MDOF Dynamics (2)

Use of Orthogonality Relationships:

& + T C y & + T Ky = T MRu &&g T M& y

T M = M * T C = C *
T K = K *

i T Mi = mi * i T Ci = ci * i T K i = k i *
T

SDOF equation in Mode i :

&i + ci* y & i + ki* yi = i MRu &&g mi* & y

Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 12

The orthogonality conditions are used to decouple the equations, resulting in one equation for each mode.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 12

Review of MDOF Dynamics (3)

Simplify by dividing through by mi and noting
*

ci* = 2 i i mi*

ki* = i2 * mi
2 i

i T MR & &i + 2 i i y & i + yi = T &&g = i u &&g y u i Mi

Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 13

Dividing through by the generalized mass in each mode produces the standard modal equation as shown. Note that this is identical to the standard SDOF equation except for the presence of the Gamma term which is referred to as the modal participation factor of the mode.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 13

Review of MDOF Dynamics (4)

T MR & &i + 2 i i y & i + 2 yi = iT &&g = i u &&g y u i Mi

Modal Participation Factor:

i T MR i = T i Mi
Important Note: i depends on mode shape scaling
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 14

It is important to note that the amplitude of the modal participation factor is dependent on the (arbitrary) modal scaling factor. This is evident from the fact that one appears in the numerator and two terms appear in the denominator.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 14

Variation of First Mode Participation Factor with First Mode Shape 1=1.0
1.0

1=1.4
1.0

1=1.6
1.0

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 15

This slide shows how the modal participation factor is dependent on the shape of the mode (which is independent of the scale factor). Note that the mode shapes have been normalized such that the top level displacement is 1.0.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 15

Review of MDOF Dynamics (5)

Any Mode of MDOF system

& &i + 2 i i y & i + i 2 yi = i u &&g y

SDOF system

&& + 2 D & + 2 D = u &&g D i i i i i i

If we obtain the displacement Di(t) from the response of a SDOF we must multiply by 1 to obtain the modal amplitude response yi(t). history

y1 (t ) = 1 Di (t )
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 16

Note that the only real difference between the a single mode of the MDOF system and the SDOF system is the modal participation factor on the RHS of the individual mode of the MDOF. Code based response spectra (used in determining the target displacement) DO NOT have the modal participation factor built in.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 16

Review of MDOF Dynamics (6)

If we run a SDOF Response history analysis:

yi (t ) = i Di (t )
If we use a response spectrum:

yi ,max = i Di ,max

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 17

The response history or response spectrum ordinate of a single mode of a MDOF system is easily obtained from the equivalent SDOF system.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 17

Review of MDOF Dynamics (7)

In general

yi (t ) = i Di (t )
Recalling

ui (t ) = i yi (t )
Substituting

ui (t ) = ii Di (t )
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 18

The natural displacement vector (e.g. nodal displacements) in any mode is given by the lower equation. This is obtained by simple algebraic manipulation of two previous equations.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 18

Review of MDOF Dynamics (8)

Applied static forces required to produce ui(t):

Fi (t ) = Kui (t ) = i Ki Di (t )
Recall

Ki = i2 Mi
where

Fi (t ) = i Mi i2 Di (t ) = i Mi ai (t )
Fi ( t ) = S i ai ( t ) Si = i Mi
Advanced Analysis 15-5b - 19

Instructional Material Complementing FEMA 451, Design Examples

Given the displacements and the elastic stiffness K, the equivalent static forces F required to produce the displacements can be obtained for any mode. The equation is manipulated to obtain the equivalent static forces in terms of pseudoacceleration and a force distribution vector S.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 19

Review of MDOF Dynamics (9)

Total shear in mode:

Vi (t ) = i (Mi ) Rai (t ) = iiT MRai (t )

T

Vi = FiT R

a (t ) Vi ( t ) = M i i
Effective Modal Mass:

= [i MR] M i iT Mi
T

Important Note: M i does NOT depend on mode shape scaling

Advanced Analysis 15-5b - 20

Instructional Material Complementing FEMA 451, Design Examples

The total shear in each mode is obtained as shown in the top equation. Some algebraic manipulation results in the effective modal mass for each mode. Note that this quantity is NOT dependent on mode shape scaling (as a pair of s appear in the numerator and the denominator). Though not evident from this slide the sum of the effective mass in all of the modes is equal to the total mass of the system.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 20

Variation of First Mode Effective Mass with First Mode Shape

M 1 = 1.0 M Total M 1 = 0.9 M Total M 1 = 0.8 M Total
1.0

1.0

1.0

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 21

This slide shows how the effective mass is dependent on mode shape. Again the modes have been normalized to a value of 1.0 at the top level. It should be noted that the first case is actually impossible for an MDOF system as all of the effective mass is in the first mode (leaving none for the higher modes).

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 21

Review of MDOF Dynamics (10)

S1 + S 2 + ... + S n = MR

S
k =1

i ,k

=M i

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 22

In may be shown that the sum of the S vectors is equal to the product of M and R. The sum of the entries in each row of each S vector is the effective modal mass in that mode. Note that i is an index over modes and k is an index over DOF.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 22

Simple Numerical Example

0 50 50 K = 50 110 60 0 60 130 1.267 S1 = 1.060 0.600 0 1.0 0 M = 0 1.1 0 0 0 1 . 2 0.071 S3 = 0.183 0.172

0.338 S 2 = 0.223 0.428

S1,k = 2.927
k =1

S 2,k = 0.313
k =1

S
k =1

3, k

= 0.060

1.0 S1 + S 2 + S3 = 1.1 1.2

Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 23

This example illustrates some of the properties of the S vectors.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 23

Review of MDOF Dynamics (11)

Displacement Response in single mode:

ui (t ) = ii Di (t )
From Response-History or Response Spectrum Analysis

Total shear in single mode:

a (t ) Vi ( t ) = M i i
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 24

We are now ready to make the appropriate transformations. The quantities Di and ai would come from a linear RHA or response spectrum analysis of the equivalent SDOF system for the i-th mode. The ui and Vi terms are the equivalent structural displacements and forces in the i-th mode. If the structural forces and displacements are known (as from a pushover analysis), the modal equivalents, Di and ai, may be determined.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 24

First Mode Response as Function of System Response

Modal Displacement:

D1 (t ) =

u1,roof (t ) 11,roof

Modal Acceleration:

a1 (t ) =

V1 (t ) M 1
Advanced Analysis 15-5b - 25

Instructional Material Complementing FEMA 451, Design Examples

Here are the final equations used to make the transformation from the pushover curve (in the base shear vs roof displacement domain) to the capacity curve (in the modal acceleration vs modal displacement domain).

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 25

Converting Pushover Curve to Capacity Curve

u1,roof (t ) = 11,roof D1 (t )
D1(t)

V
First Mode SDOF System (modal coords)

V
First Mode System (natural coords)
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 26

Here the transformation from first mode system natural displacement coordinates to first mode modal amplitude coordinates is shown.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 26

Converting Pushover Curve to Capacity Curve

Base Shear
V (t ) Modal a1 ( t ) = 1 Acceleration M 1

Pushover Curve Roof Displacement

Capacity Curve

Modal Displacement u (t ) D1 (t ) = 1 11,roof

Advanced Analysis 15-5b - 27

Instructional Material Complementing FEMA 451, Design Examples

Finally, the first mode capacity curve is obtained from the pushover curve through the use of the transformation equations determined on the last several slides. We will get back to the use of the capacity curve later.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 27

Development of Pushover Curve

Potential Plastic Hinge Location (Must be predicted and possibly corrected)

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 28

Now it is necessary to discus the development of the pushover curve itself. In the development of the curve it is first necessary to develop a realistic nonlinear model of the system. All possible sources of inelastic deformation should be included in the analytical model. If it is found during analysis that sections that were not modeled inelastically develop forces or moments in excess of yield capacity the model should be modified to include such behavior and the analysis should be rerun.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 28

Event-to-Event Pushover Analysis

Create Mathematical Model Apply gravity load to determine initial nodal displacements and member forces Apply lateral load sufficient to produce single yield event Update nodal displacements and member forces Modify structural stiffness to represent yielding

Continue Until Sufficient Load or Displacement is obtained.

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 29

This is the basic flowchart for event-to-event pushover analysis. Each step will be explained in more detail in later slides. Note that the analysis may be performed under force control or under displacement control. Displacement control is required if the tangent stiffness matrix of the structure is not positive definite at any step (usually the latter steps). Note that this sequence assumes that no yielding occurs under gravity load. (If it does, the structure should be redesigned!)

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 29

Initial Gravity Load Analysis

Moment

MG

A,B
Rotation

Moments plotted on tension side.

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 30

The first step in any pushover analysis is to run a gravity analysis. It is extremely rare that yielding will occur in the gravity analysis, however the pattern of moment and forces that develop in the individual structural components will have an effect on the location of and sequencing of hinges in the lateral load phase of the analysis. The gravity load analysis will also cause gravity related P-Delta effects to be activated (if such effects are explicitly included in the analytical model).

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 30

Analysis 1: Gravity Analysis

Base Shear

Roof Displacement
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 31

This slide shows the state of the structure just after gravity loads are applied but before any lateral load has been applied.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 31

Lateral Load Analysis (Acting Alone)

ML
Moment

Total Load = V

A B

B A Rotation

Moments plotted on tension side.

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 32

Now the lateral load is applied. The idealized moments in two potential hinging regions are shown for the lateral load only. Insufficient lateral load has been applied to cause yielding.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 32

Combined Load Analysis Including Total Load V

MG+ML
Moment

B Total Load = V
A B

Rotation

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 33

If the member forces from gravity load are added to the member forces from the lateral loads it is seen that the moment computed at the right span, right hinge is well in excess of the capacity. The program performing the analysis will then compute the fraction of the lateral load, that when added to the gravity load, causes first yielding in the structure.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 33

Analysis 2a First Lateral Analysis

Base Shear

Roof Displacement
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 34

Here the total load V is applied to the structure which has not yet yielded.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 34

Combined Load Analysis:

Determine amount of Lateral Load Required to Produce First Yield MG+ ML
Moment

Total Load = V
Rotation

For all potential hinges (i) find i such that

M G ,i + i M L ,i = M P , i
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 35

We have applied too much lateral load. Hence, we want to compute the portion of load, V, that just causes the first yielding.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 35

Analysis 2b Adjust Load to First Yield

Base Shear
Old Tangent Stiffness

New Tangent Stiffness

1 1

Roof Displacement
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 36

The pushover curve is not at the state shown, with only one hinge present.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 36

Analysis 3a
Modify System Stiffness Apply Remainder of Load

Base Shear

VR=V(1-1)
1 1

Roof Displacement
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 37

We now apply the remainder of the load VR = V(1-). We will want to determine how much of the remaining load causes the next hinge to form.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 37

Determine amount of Lateral Load Required to Produce Second Yield

Total Load = (1-1)V

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 38

The next hinge will occur at the right of the second story girder of the right bay.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 38

Analysis 3b Adjust Load to Second Yield

Base Shear
Old Tangent Stiffness New Tangent Stiffness

2 1 2 1

Roof Displacement
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 39

The second hinge is formed and the stiffness is changed.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 39

Analysis 4a
Modify System Stiffness Apply Remainder of Load
Base Shear

VR
2 1 2 1

Roof Displacement
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 40

The remaining load is applied, and the next hinge location is found.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 40

Analysis 4b Adjust Load to Third Yield

Base Shear

Old Tangent Stiffness New Tangent Stiffness

3 2 1 3 2 1

Roof Displacement
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 41

It appears that the next hinge will form at the right hand side of the first story girder in the left bay.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 41

Analysis 5a..
Base Shear

VR
2 1

2 3 1

Roof Displacement
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 42

The procedure is continued until adequate displacement has been obtained. A maximum expected displacement would be 3% of the height of the structure (as this is in excess of the seismic drift limit in most codes). The Provisions and a few other documents require that the pushover curve be extended to 1.5 times the Target Displacement, where the Target Displacement is determined empirically. The empirical expressions for computing Target Displacement are discussed later.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 42

Convert Pushover Curve to Capacity Curve

Base Shear
V (t ) Modal a1 ( t ) = 1 Acceleration M 1

Pushover Curve Roof Displacement

Capacity Curve
ACTUAL SIMPLIFIED

Modal Displacement u (t ) D1 ( t ) = 1

11,roof

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 43

In the capacity-spectrum approach it is necessary to transform the pushover curve (in Force-Displacement space) into a Capacity Curve (in Modal Acceleration-Modal Displacement Space). The transformation equations developed earlier are used for this purpose. In many cases it is convenient to replace the capacity curve by a simplified bilinear version as shown. We will use the bilinear version in subsequent discussions.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 43

Equivalent Viscous Damping

Base Shear/Weight or Pseudoacceleration (g) Elastic Spectrum based demand curve for X% equivalent viscous damping

Point on capacity curve representing X% Equivalent Viscous Damping.

Target Displacement

Spectral Displacement
Advanced Analysis 15-5b - 44

Instructional Material Complementing FEMA 451, Design Examples

The next step is to determine the target displacement. Here, we are using the Capacity Curve in association with a Demand Curve which is an elastic response spectrum modified for site effects, and then modified for an amount of equivalent viscous damping, X, which is consistent with an amount of hysteretic energy dissipated by the system. It is important to note here that the demand spectrum is plotted as pseudoacceleration vs displacement, not pseudoacceleration vs period as is traditional. We will get back to this later. Before doing so it is instructive to discuss the concept of equivalent viscous damping which is a key (yet dubious) ingredient in the procedure.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 44

Computing Damping Ratio from Damping Energy and Strain Energy

F FD Area=ED u

F FS
K

Area=ES u

E D = FD u
= Cu 2 = 2m 2u 2

ES = 0.5Fs u = 0.5 Ku 2

= 0.5m 2u 2
ED 4ES
Advanced Analysis 15-5b - 45

Instructional Material Complementing FEMA 451, Design Examples

This is the derivation of damping ratio as computed on the basis of damping energy and strain energy. It is applicable only to systems under harmonic resonance.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 45

Computing Damping Ratio from Damping Force and Elastic Force

F F FD Area=ED u FS Area=ES u

E D = FD u

ES =

1 Fs u 2

ED F = D 4ES 2 FS
Advanced Analysis 15-5b - 46

Instructional Material Complementing FEMA 451, Design Examples

A manipulation of terms in the previous slide shows that for a system under harmonic resonance, the damping ratio may be expressed as 0.5 times the ratio of the damping force to the elastic spring force. This equation holds at any time in the response history after the transients have fully damped out.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 46

Computing True Viscous Damping Ratio from Damping Energy and Strain Energy

ED F = D 4ES 2 FS

Note: System must be in steady state harmonic RESONANT response for this equation to work.

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 47

This slide is here to emphasize the key restriction on the previous derivations.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 47

Harmonic Resonant Response from NONLIN

15000

10000

5000 Force, Kips

-5000

-10000

Spring Damper

-15000 -6 -4 -2 0 2 4 6 Displacement, Inches

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 48

This is an example of a system response at harmonic resonance. The damping energy is computed from the perimeter of the ellipse (shown in black).

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 48

Harmonic Non-Resonant Response from NONLIN

2000

1500

1000

500 Force, Kips

-500

-1000

-1500

Spring Damper

-2000 -1.00

-0.50

0.00 Displacement, Inches

0.50

1.00

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 49

This is a similar plot for a nonresonant system. There is no true elliptical steady state response.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 49

Results from NONLIN Using: =

System Period = 0.75 seconds Harmonic Loading Target Damping Ratio 5% Critical Loading Period (sec) 0.50 0.75 1.00 Damping Force (k) 118 984 197 Spring Force (k) 787 9828 2251

ED F = D 4ES 2 FS

Damping Ratio % 7.50 5.00 3.75

Resonant

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 50

Here is a numerical example for a system with 5% damping at resonance. The use of the upper equation produces various results, depending on the ratio of the loading period to the system period. Only the resonant response produces the correct damping ratio.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 50

Results from NONLIN Using: =

System Period = 0.75 seconds Harmonic Loading Target Damping Ratio 20% Critical Loading Period (sec) 0.50 0.75 1.00 Damping Force (k) 430 999 1888 Spring Force (k) 717 2498 5666

ED F = D 4ES 2 FS

Damping Ratio % 30.0 20.0 16.7

Resonant

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 51

Here, a similar analysis was performed, except the target damping ratio was 20% of critical. Again, only the resonant response produces the correct damping ratio.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 51

Computing Equivalent Viscous Damping Ratio from Yield-Based Hysteretic Energy and Strain Energy
ES ED u EH u ES

Viscous System

Yielding System

EH 4ES
Advanced Analysis 15-5b - 52

Instructional Material Complementing FEMA 451, Design Examples

Now, an attempt is made to compute an equivalent viscous damping ratio for a system whose energy dissipation is hysteretic, rather than viscous. It is very important to note that in the viscous system the elastic energy is based on the initial stiffness of the system, whereas in the yielding system the secant stiffness is used.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 52

Actual Yielding System

System Stiffness and Energy Dissipation
Rigid Rigid

Equivalent Elastic System

System Stiffness

System Energy Dissipation

ES EH u ED

ES u
Advanced Analysis 15-5b - 53

Instructional Material Complementing FEMA 451, Design Examples

This slide compares the conceptual differences between equivalent viscous damping for hysteretic systems, and equivalent viscous damping in elastic systems.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 53

Original Yielding System

m
Initial Stiffness k
Initial Frequency:

Posin(Inelt) =

k m

Resonant Frequency:

Inelastic =
F Fy k uy ksec umax u

2 ( 0.5 sin 2 )
uy umax )

= cos 1 (1 2

Maximum Steady State Response (loaded at Inelastic):

umax =

4u y P 4 o ku y
Advanced Analysis 15-5b - 54

Instructional Material Complementing FEMA 451, Design Examples

Here we show some actual properties of the yielding system. The formulas shown here are from the 1968 article by Jennings.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 54

Equivalent Elastic System

m Initial Stiffness k Posin(t) Resonant Frequency:

sec =

ksec m

Maximum Steady State Resonant Response: F Fy ksec uy umax u

u max = Po 2 sec ksec

sec

Equivalent Damping: uy 1 = 0.637(1 ) = 0.637(1 ) umax

Advanced Analysis 15-5b - 55

Instructional Material Complementing FEMA 451, Design Examples

These are the properties of the equivalent system. Note the relationship between equivalent damping and ductility demand.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 55

Equivalent Elastic System when Strain Hardening is Included

F

K
K u uy Fmax Fy ksec umax

sec 0.637

( Fy umax Fmax u y )

Fmax umax
Advanced Analysis 15-5b - 56

Instructional Material Complementing FEMA 451, Design Examples

For a system with strain hardening, the equivalent viscous damping is computed as shown here. Note that the quantity is the strain hardening ratio.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 56

Equivalent Elastic System when Strain Hardening is Included

F

K
K u uy Fmax Fy ksec umax

Equiv 0.637

uy Fy 1 1 = 0 . 637 Fmax u max ( 1) + 1

Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 57

This is a recast of the previous formula in terms of ductility and strain hardening ratio.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 57

Effect of Secondary Stiffness On Equivalent Viscous Damping

0.6

=0.0 Equivalent Viscous Damping Ratio

0.5

0.4

=0.05 =0.10 =0.15 =0.20

0.3

0.2

0.1

0 1 2 3 4 5 6 7 8 9

Ductility Demand
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 58

This is a plot of equivalent viscous damping ratios versus ductility demand for systems with different strain hardening ratios. Note that systems with high strain hardening ratios actually have reduced equivalent damping ratios at larger ductility demands.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 58

Reduction in Ductility Demand with Strain Hardening Ratio

(W = 11250 k, K = 918 k/in., T=1.0 sec, El Centro Ground Motion)
8
Yield Strength

7 6

600 800 1000

Ductility Demand

5 4 3 2 1 0 0.00 0.05 0.10 0.15 0.20 0.25

Strain Hardening Ratio

Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 59

In practical cases the reduction in damping with larger ductility ratios and strain hardening ratios is not likely a problem due to the fact that (as expected) the ductility demand reduces with strain hardening ratio. For systems with high strain hardening it is unlikely that the ductility demand will be high enough to indicate decreasing damping with increased ductility demand.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 59

Total System Damping (% Critical)

Total = 5 + Equiv

Robust Shaking Duration Short Long

Moderately Robust

Pinched Or Brittle

=1 = .7

= .7 = .33

= .7 = .33
Advanced Analysis 15-5b - 60

See ATC 40 for Exact Values

Instructional Material Complementing FEMA 451, Design Examples

In the previous derivations it was assumed that the inelastic hysteretic behavior is robust. For systems with less robust behavior the energy dissipated per cycle will be reduced. Also, the longer the duration of strong shaking, the more the likelihood of reduced stiffness and strength. The kappa term of ATC 40 compensates for this effect. The total damping in a system is the 5% inherent damping plus kappa times the additional equivalent viscous damping from inelastic energy dissipation.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 60

Equivalent Viscous Damping Values for EPP System

(Values Shown are Percent Critical)
1.2

1.0

5 10 20

30

40

0.8

F/Fy

0.6

0.4

0.2

0.0 0.0 0.5 1.0 1.5 2.0 2.5

U/Uy
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 61

This plot shows the equivalent viscous damping of 5 to 40 percent superimposed in a normalized bilinear capacity curve. This follows the formula that damping=5+63.7(1-1/mu). Note that the equivalent damping at yield is 5% (as expected for an elastic system). Note also that equivalent damping increases very rapidly with ductility demand. At a ductility demand of only 2, equivalent viscous damping is about 35% of critical.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 61

Equivalent Viscous Damping Values for System

With 5% Strain Hardening Ratio
(Values Shown are Percent Critical)
1.2 1.0 0.8 F/Fy 0.6 0.4 0.2 0.0 0.0 0.5 1.0 1.5 u/uy
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 62

5 10 20

30

40

2.0

2.5

3.0

This is similar to the previous plot, but strain hardening is 5% of the initial stiffness. Note that equivalent viscous damping is somewhat less than 35% critical when the ductility demand is 2.0.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 62

Pseudoacceleration Spectrum in Traditional Format

1.2 1.0

SDS

Pseudoacceleration, g

0.8

0.6

0.4

SD1/T

0.2

0.0 0.0 1.0 2.0 Period, Seconds 3.0 4.0

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 63

Now we are ready to discuss the demand curve. Shown is a the NEHRP response spectrum normalized to a maximum value of 1.0 g. Recall that this spectrum includes site effects as well as the 2/3 factor to account for expected good behavior. In the western U.S. this is equivalent to a 10% in 50 year earthquake. In the eastern U.S. it is closer to a 5% in 50 year earthquake. Note that the response spectrum is plotted in the traditional manner of pseudoacceleration vs period of vibration.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 63

Pseudoacceleration (Demand) Spectrum in ADRS Format (5% Damping)

1.2
T=.10 T=.50 T=1.0

1.0 Pseudoacceleration, g 0.8 0.6 0.4

T=1.5

T=2.0

T=3.0

0.2
T=4.0

0.0 0.0 5.0 10.0 15.0 20.0 Spectral Displacement, Inches

Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 64

In this plot the response spectrum is plotted in so-called AccelerationDisplacement space, hence the name Acceleration Displacement Response Spectrum (ADRS). Here, period values are shown as diagonal lines.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 64

Spectral Reduction Factors for Increased Equivalent Damping

1.2 Robust Degrading Severely Degrading 0.8 1.0

Spectral Reduction Factor

0.6

0.4

0.2

0.0 0

W E I V E R
10 20

E D I L S this
30 40 50

Equivalent Viscous Damping, Percent Critical

Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 65

The previous spectrum was computed for a system with 5% damping. For higher levels of viscous damping the resulting displacements and accelerations will be lower. It is unclear why different curves are provided for systems with degrading and severely degrading response a response spectrum is by definition based on a linear elastic analysis of a system with a certain level of damping. Also, the effect of hysteretic behavior on damping is already included in the kappa factor shown on a previous slide.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 65

Demand Spectra for Various Damping Values

1.2
T=.50 5% Damping T=1.0

1.0
T=1.5

Pseudoacceleration, g

10%

0.8
20%

0.6

30% 40% T=2.0

0.4
T=3.0

0.2
T=4.0

0.0 0 5 10 15 20 Spectral Displacement, Inches

Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 66

Here we show a family of demand spectra for various damping values ranging from 5 to 40% of critical. At any period value the displacement and the pseudoaccelerations are significantly reduced as damping increases.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 66

Combined Capacity-Demand Spectra

1.2
5% Damping

1.0 Pseudoacceleration, g
10%

0.8
20% 5% 10% 20% 30% 40%

0.6

30% 40%

0.4

0.2

0.0 0 5 10 15 20 Spectral Displacement, Inches

Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 67

Now the demand spectrum and the capacity spectrum are plotted on the same graph. This is the advantage of making the initial pushover curve transformation from force-displacement to modal acceleration- modal displacement. Here, diagonal lines are used to label the damping values on the capacity plot. We will zoom in on the area in the circle on the next slide.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 67

Finding the Target Displacement

0.4
40% 30% 20% 10% 5% Damping

5% 10%

Pseudoacceleration, g

0.3

20%

30%
0.2

40%
0.1

0.0 0 2 4 6

Target Disp = 6.2 in.

8 10 Spectral Displacement, Inches
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 68

Finally, we are able to determine the target displacement. It is the displacement corresponding to the point where the X percent demand spectra and the X percent damping point on the capacity spectrum intersect. This point is found graphically on this plot. It may also be found by iteration or for simple capacity and demand curves, a closed-form solution may be found. Recall that this displacement corresponds to the first mode SDOF system and must therefore be transformed back to the displacement within the MDOF system.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 68

You are Not Done Yet!

Note: The target displacement from the Capacity-Demand
diagram corresponds to a first mode SDOF system. It must be multiplied by the first mode modal participation factor and the modal amplitude of the first mode mode shape at the roof to determine displacements or deformations in the original system. Hinge rotations may then be obtained for comparison with performance criterion.

Knowing the target displacement, the base shear can

be found from the original pushover curve.

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 69

These are the final steps in the process.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 69

There is sometimes cause to fear that scientific technique, that proud servant of engineering arts, is trying to swallow its master
Professor Hardy Cross

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 70

It has taken a lot of work to find the one target displacement. Such an analysis would need to be repeated for different lateral load patterns, and possibly for various perturbations in hysteretic behavior. Fortunately, several commercial programs (e.g SAP2000, RAM PERFORM) make the process relatively simple. However, there are many many assumptions and simplifications involved in the process, and one might wonder if a more simple approach could be used without tremendous loss in accuracy. Two such simplified approaches are given in the next several slides.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 70

Simplified Pushover Approaches: 2003 NEHRP Recommended Provisions

Procedure is presented in Appendix to Chapter 5 Gravity Loads include 25% of live load (but Provisions are not specific on P-Delta Modeling Requirements) Lateral Loads Applied in a First Mode Pattern Structure is pushed to 150% of target displacement Target displacement is assumed equal to the displacement computed from a first mode response spectrum analysis, multiplied by the factor Ci Ci adjusts for error in equal displacement theory when structural period is low

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 71

The first of the simplified procedures is given by the 2003 NEHRP Provisions. A list of the basic assumptions is presented on this slide. The use of the first mode displacement follows the equal displacement observation that is the displacements predicted from an elastic and inelastic analysis of the same structure are approximately the same. The NEHRP response spectrum is used in the displacement part of the analysis.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 71

Simplified Pushover Approaches: 2003 NEHRP Provisions (2)

Ci =

(1 Ts / T1 ) + (Ts / T1 ) Rd
Sa

Ci=1 if Ts/T1 <1

TS = S D1/ S DS

Rd =

1.5 R 0

Ts

T
Advanced Analysis 15-5b - 72

Instructional Material Complementing FEMA 451, Design Examples

Ci is a correction factor for very short period structures for which it has been observed that the equal displacement approach is not particularly reliable.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 72

Simplified Pushover Approaches: 2003 NEHRP Provisions (3)

Member strengths need not be evaluated Component deformation acceptance based on

laboratory tests Maximum story drift may be as high as 1.25 times standard limit Nonlinear Analysis must be Peer Reviewed

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 73

The first statement is a testament to the old seismic design adage strength is essential but otherwise unimportant. The Provisions do not provide acceptance criteria for component deformations. It is suggested by the Provisions that such limits be based on available test data. Instead, performance criteria from ATC40 or FEMA 356 may be used. An Advantage to performing pushover analysis (or nonlinear response history analysis) is that the allowable story drifts are increased by 25%.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 73

Simplified Pushover Approaches: FEMA 356*. (Also used in FEMA 350)

Procedure presented in Chapter 3 More detailed (thoughtful) treatment than in NEHRP Recommended Provisions Principal Differences: > Apply 25% of unreduced Gravity Load > Use of two different lateral load patterns > P-Delta effects included > Consideration of Hysteretic Behavior
* FEMA 273 in Prestandard Format

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 74

FEMA 356 gives a similar but somewhat more detailed procedure. The principal differences are shown here.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 74

Simplified Pushover Approaches: FEMA 356 (2)

T t = C0C1C2C3 S a e 2 g 4
t = Target Displacement

Spectral Displacement

C0 = Modification factor to relate roof displacement to first mode spectral displacement. C1 = Modification factor to relate expected maximum inelastic displacement to displacement calculated from elastic response (similar to NEHRP Provisions Ci)

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 75

In FEMA 356, the target displacement is computed from a product of multipliers on the spectral displacement. The first multiplier, Co, is effectively the first mode participation factor times the ordinate of the mode shape at the roof of the structure. The second multiplier, C1, accounts for errors in the equal displacement concept for low period buildings.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 75

Simplified Pushover Approaches: FEMA 356 (3)

T t = C0C1C2C3 S a e 2 g 4
C2 = Modification factor to represent effect of pinched hysteretic loop, stiffness degradation, and strength loss. C3 = Modification factor to represent increased displacements due to dynamic P-Delta effect

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 76

Multiplier C2 is an adjustment for hysteretic behavior, and C3 is a modifier for P-Delta effects.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 76

Discussion of Assumptions
1. Dynamic effects are ignored 2. Duration effects are ignored 3. Choice of lateral load pattern 4. Only first mode response included 5. Use of elastic response spectrum 6. Use of equivalent viscous damping 7. Modification of response spectrum for higher damping

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 77

This is a list of the most pertinent (glaring) assumptions in pushover analysis. The list is quite long, and because of these issues, many engineers and researchers believe that pushover analysis is not a viable analysis/design tool. This is difficult to argue with. However, a pushover analysis does provide more information than does a purely elastic analysis. In particular it is beneficial to know the sequencing of hinging and to develop a true estimate of overstrength. Pushover analysis can be considered as a useful evaluation tool, not to be used alone, but used in concert with other tools to assess the likely performance of a structure. This being said, some of the potential problems with pushover analysis are illustrated in the following slides.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 77

True Acceleration vs Pseudoacceleration 30% Critical Damping

200 True Total Acceleration Pseudo Acceleration Acceleration, in/sec/sec 150

100

50

0 0.0 1.0 2.0 3.0 Undamped Period, Seconds

Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 78

The main issue with pushover analysis is the use of the equivalent viscous damping to predict the response of a yielding system. One reason for concern is that we use a pseudo acceleration spectrum as the basis for predicting displacements in a SDOF system. Recall that the pseudoacceleration spectrum is derived from the true displacement spectrum by dividing each displacement value by the frequency squared. For low damping values, there is negligible error in this assumption. For higher damping values the error can be quite large, particularly when the system has a long period of vibration. This plot shows the comparison of a pseudoacceleration spectrum and a true acceleration spectrum for a system with 30% critical damping. At the higher period values the error is 30%, on the unconservative side.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 78

Relative Error Between True Acceleration and Pseudoacceleration

50% 5% Damping 10% Damping 15% Damping 20% Damping 25% Damping 30% Damping

40%

Error, Percent

30%

20%

10%

0% 0.00

1.00

2.00

3.00

Period of Vibration, Seconds

Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 79

This plot shoes the error in computed pseudoacceleration and true acceleration for a range of different damping values. In general, the higher the damping value the larger the error, and the longer the period of vibration, the higher the error.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 79

Equivalent Elastic System

m Initial Stiffness k Posin(t) Resonant Frequency:

sec =

ksec m

F Fy ksec umax u

Maximum Steady State Resonant Response: k umax = sec 2 sec Po

Equivalent Damping: u sec = 0.637(1 y ) umax
Advanced Analysis 15-5b - 80

Instructional Material Complementing FEMA 451, Design Examples

Another big concern is the use of equivalent viscous damping itself. This is a review of how the equivalent damping is computed.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 80

These systems have the same hysteretic Energy Dissipation, the same AVERAGE (+/-) displacement, but considerably DIFFERENT maximum displacement.

umin

umax

umin

umax

umin

umax

The equivalent viscous damping (see previous slide) is good at predicting the AVERAGE displacement, but CAN NOT predict the true maximum displacement.
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 81

Here, three systems that have the same hysteretic energy have considerably different maximum displacements due to residual deformation. This effect can occur for a variety of reasons, some having to do with ground motion, and others with structural properties. Jennings pointed this out in his original paper that forwarded the use of equivalent damping. To date, this effect is not included in pushover analysis. It is included automatically in response history analysis.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 81

Improved Pushover Methods

Use of Inelastic Response Spectrum Adaptive Load Patterns Use of SDOF Response History Analysis Inclusion of Higher Mode Effects

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 82

These are several of the ideas that have been used to improve pushover analysis. In most cases the desire is to produce a predicted response that more closely matched that from response history analysis.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 82

Elastic Spectrum Based Target Displacement

Base Shear/Weight or Pseudoacceleration (g)
Elastic Spectrum based demand curve for X% equivalent viscous damping

Point on capacity curve representing X% equivalent viscous damping.

Target Displacement

Spectral Displacement
Advanced Analysis 15-5b - 83

Instructional Material Complementing FEMA 451, Design Examples

One approach is to use ductility, not equivalent damping, to determine the target displacement. Here the concept of using equivalent damping is reviewed.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 83

Inelastic Response Spectrum Based Target Displacement

Base Shear/Weight or Pseudoacceleration (g)

Inelastic Spectrum based Demand Curve for ductility demand of X.

Point on capacity curve representing ductility demand of X.

Target Displacement

Spectral Displacement
Advanced Analysis 15-5b - 84

Instructional Material Complementing FEMA 451, Design Examples

Instead of using equivalent viscous damping, one can use ductility as a modifier for the elastic spectrum. The target displacement is now found as the projection on the displacement axis of the point at which a spectrum with a ductility of X meets the ductility of X on the capacity curve. It has been shown by Chopra and others that this produces more consistent results when compared to nonlinear response history analysis.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 84

Inelastic Spectrum Based Target Displacement

Gives the same results as the equal displacement

theory for (longer period) EPP systems

When compared to inelastic response history

analysis, the use of inelastic spectra gives better results than ATC 40 procedure.

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 85

This is a summary of the advantages of using the inelastic spectrum.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 85

Computing Target Displacements from Response History Analysis of SDOF Systems

Method called Uncoupled Modal Response History Analysis

(UMRHA) is described by Chopra and Goel. See, for example, Appendix A of PEER Report 2001/03, entitled Modal Pushover Analysis Procedure to Estimate Seismic Demands for Buildings.

In the UMHRA method, the undamped mode shapes are

used to determine a static load pattern for each mode.

Using these static lateral loads, a series of pushover curves

and corresponding bilinear capacity curves are obtained for the first few modes. This is done using the procedures described earlier for the ATC 40 approach.

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 86

Another approach which uses inelastic response history analysis of uncoupled SDOF inelastic systems has been suggested by Chopra. Again, this procedure appears to produce better results that do the more basic approaches. The steps in the procedure are outlined in this and the next slide.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 86

Computing Target Displacements from Response History Analysis of SDOF Systems (2)

Using an appropriate ground motion, a nonlinear dynamic

response history analysis is computed for each modal bilinear system. This may be accomplished using NONLIN or NONLIN-Pro.

The modal response histories are transformed to system

coordinates and displacement (and deformation) response histories are obtained for each mode.

The modal response histories are added algebraically to

determine the final displacement (deformations). In the Modal Pushover approach, the individual response histories are combined using SRSS.
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 87

Continuation of previous series.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 87

Computing Target Displacements from Response History Analysis of SDOF Systems (3)

Results from such an analysis are detailed in PEER Report

2001/16, entitled Statistics of SDF-System Estimate of Roof Displacement for Pushover Analysis of Buildings. Conclusions from above report (paraphrased by F. Charney): For larger ductility demands the SDOF method, using only the first mode, overestimates roof displacements and the bias increases for longer period buildings. For small ductility demand systems, the SDOF system, using only the first mode, underestimates displacement, and the bias increases for longer period systems.
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 88

These are the conclusions from the report describing the method.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 88

Conclusions (continued)
First mode SDOF estimates of roof displacements due to individual ground motions can be alarmingly small (as low as 0.31 to 0.82 times exact) to surprisingly large (1.45 to 2.15 times exact). Errors increase when P-Delta effects are included. (Note: the method includes P-Delta effects only in the first mode). The large errors arise because for individual ground motions the first mode SDOF system may underestimate or overestimate the residual deformation due to yield-induced permanent drift. The error is not improved significantly by including higher mode contributions. However, the dispersion is reduced when elastic or nearly elastic systems are considered.
Instructional Material Complementing FEMA 451, Design Examples Advanced Analysis 15-5b - 89

Continuation of the conclusions.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 89

Computing Target Displacements from Response History Analysis of SDOF Systems Problems with the method:

No rational basis
Does not include P-Delta effects in higher modes Can not consider differences in hysteretic behavior of
individual components

No reduction in effort compared to full time-history analysis Problem of ground motion selection and scaling still exists

Instructional Material Complementing FEMA 451, Design Examples

Advanced Analysis 15-5b - 90

The problem with the improved approach is that it has absolutely no rational basis. If one is going to go through all the effort indicated, one might as well use response history analysis which is the subject of the following set of slides. Proceed to Topic 15-5c.

FEMA 451B Topic 15-5b Notes

Advanced Analysis 15-5b - 90