SEAU Practical Uses of Nonlinear Pushover Analysis
SEAU Practical Uses of Nonlinear Pushover Analysis
Disclaimer
This presentation is meant to neither explicitly endorse nor discourage
the use of the nonlinear pushover analysis method for any given design
scenario. As with any method, nonlinear pushover has its benefits and
limitations. Whether nonlinear pushover analysis is appropriate lies at
the discretion of the engineer in responsible charge.
1
2/3/2017
Presentation Outline
1. What is Nonlinear Analysis?
2. What is Nonlinear Pushover Analysis?
3. Why Use Nonlinear Pushover Analysis?
4. What are the Common Perceptions?
5. Embracing Innovation.
6. Embracing Nonlinearity.
7. Nonlinear Pushover Analysis – A Good First Step.
8. Nonlinear Element Modeling
9. Examples
10. Nonlinear Pushover Limitations
11. Time History – Quick Summary
12. Software
Wikipedia:
In physical sciences, a nonlinear system is a system in which the output is
not directly proportional to the input.
2
2/3/2017
Types of Nonlinearity
• Geometric Nonlinearity – Also known as P-Delta.
• Rate Dependent Nonlinearity – Properties change as a function of rate
(e.g. viscous damper)
• History Dependent Nonlinearity – Properties change as a function of
repeated load (cumulative ductility or fatigue).
• Cyclical Dependent Nonlinearity – Change of hysteretic properties
(e.g. tension only braces).
• Contact Nonlinearity (e.g. footing/soil interface, pounding)
• Material Nonlinearity – e.g. yielding of steel.
Geometric Nonlinearity
P P
θ2 u1 2 θ3
This example is for illustrative F Axial loads in members (1, 3)
purposes and assumes equal change as a result of load F.
member lengths and section 1 3
properties
Note: The P1, P2 and P3 forces are the developed frame forces which are resolved through iterative
processes.
3
2/3/2017
Stocky
Concentric Braces
4
2/3/2017
Contact Nonlinearity
Footing Deflection
8
0
0 5 10 15 20 25 30 Pounding
-1
Material Nonlinearity
Some materials are more ductile than others….
Some materials behave better when confined….
5
2/3/2017
6
2/3/2017
Element Nonlinearity
Repeated cycles of stable (material) nonlinearity are preferred and even encouraged in
the codes. This is reflected in the R factors. Let’s observe…
F F F F
∆ ∆ ∆ ∆
R RR R R R R R R R R R ??
BRB/SRMF SCSW OCBF
7
2/3/2017
1000
800
Base Shear
600
Nonlinear
rxn rxn 400 Linear
200
8
2/3/2017
9
2/3/2017
10
2/3/2017
• More than 75% of respondents said that guidance for most of these topics [nonlinear
analysis and design] was ambiguous.
• Claims of inadequate software (21%), too complicated (29%), not practical/time
consuming (61%), lack of research (22%), lack of guidelines (43%).
• More guidance is needed for: Modeling NL Elements (42%), NL Procedures (18%),
Benchmark problems with solutions (35%).
11
2/3/2017
“Nonlinear analysis can be used when owners request ways to reduce costs (for new
construction) by optimizing material use, more likely though as demonstrating a building
retrofit is perhaps not even necessary (or if it is, that only minor changes/systems are
needed rather than what the code would require), or even as a way of quantifying
performance for owners, insurance and risk managers that may look at inventories of
structures. So, the major question to the owner is whether they’d prefer to pay now or pay
later?”
12
2/3/2017
5. Embracing Innovation
13
2/3/2017
A Tip Calculator?
14
2/3/2017
What do we do…fundamentally?
Design Buildings?
How?
Mathematical Models…
15
2/3/2017
Mathematical Models:
16
2/3/2017
Mathematical Models:
Mn=FyZx Mn=Asfy(d-a/2)
17
2/3/2017
This law is now used in the semiconductor industry to guide long-term planning and to set
targets for research and development.
What does this say about the advancement of computer power over the last 20 years?
18
2/3/2017
Moore’s Law:
https://en.wikipedia.org/wiki/Moore%27s_law
Embracing Nonlinearity
19
2/3/2017
6. Embracing Nonlinearity
Lumped Plasticity
Models
Embracing Nonlinearity
Professor Ed Wilson:
20
2/3/2017
Embracing Nonlinearity
21
2/3/2017
91% Mass
Activation
22
2/3/2017
The simple calculation of story drifts is not appropriate when using dynamic
response spectrum analyses because combined modal displacements can be
less than the combined modal story drift (signs are lost). This is a
consequence of simplified approaches.
SEAU: 2017 Education Conference
In other
words, what
0.38 Seconds
happens to
0.28 Seconds
these mode
shapes when
yielding
occurs?
91% Mass
Activation
23
2/3/2017
A Nonlinear Model:
24
2/3/2017
A Nonlinear Model
Nonlinear models explicitly include nonlinear mechanisms or nonlinear material
behaviors that become activated upon reaching a theoretical load or displacement
threshold.
Nonlinear
Linear:
Linear
σ = εE
Nonlinear:
σ = ??
SEAU: 2017 Education Conference
Stress Diagrams
Fb=M/S Fy Fy
Mn=FyZx
25
2/3/2017
Stress Diagrams
σc σc 0.85f’c
Mn=Asfy(d-a/2)
fy
σs fy
26
2/3/2017
Stress Diagrams
εc εc εcu=0.003
εs εt>0.00207 εt>0.00207
Elastic Elasto-Plastic Plastic
Strain Diagrams
SEAU: 2017 Education Conference
N.A.
d/2
Strain(ε)
φ
Curvature: φ = ε/(d/2)
SEAU: 2017 Education Conference
27
2/3/2017
200
Moment (k-ft)
150
100
50
0
0 0.0002 0.0004 0.0006 0.0008 0.001 0.0012
Curvature (1/in)
250
200
50
0
0 0.0002 0.0004 0.0006 0.0008 0.001 0.0012
Curvature (1/in)
28
2/3/2017
250
200
150
Moment (k-ft)
0
0 0.0002 0.0004 0.0006 0.0008 0.001 0.0012
Curvature (1/in)
29
2/3/2017
450
• Time 15000
Base Force (kips)
400
350 • Complexity 10000
300
• Project Needs
Beam Force (kip-in)
5000
250
200
150
• $$ -0.05 -0.04 -0.03 -0.02 -0.01
0
0 0.01 0.02 0.03 0.04 0.05
50
-10000
0
0 10 20 30 40 50 60 -15000
30
2/3/2017
Nonlinear Analysis
Fy 0.85f’c
Mn=FyZx Mn=Asfy(d-a/2)
fy
These very familiar models represent the ultimate strength limit state of these beams.
SEAU: 2017 Education Conference
31
2/3/2017
W14x159
Force
32
2/3/2017
Nonlinear Pushover:
∆
Force Pattern
Force
Reaction
Nonlinear Pushover:
33
2/3/2017
Nonlinear Pushover:
Nonlinear Pushover:
34
2/3/2017
Nonlinear Pushover:
Distributed Plasticity:
Reinforcement
Unconfined Concrete
Confined Concrete
35
2/3/2017
36
2/3/2017
37
2/3/2017
Comparison:
Distributed Plasticity Lumped Plasticity
Strain Hardening?
Distributed Lumped
SEAU: 2017 Education Conference
38
2/3/2017
Comparison
Distributed Lumped
SEAU: 2017 Education Conference
Key Differentiators:
• Time
• Complexity
• Project Needs
• $$
• Accuracy
• Software
39
2/3/2017
Nonlinear Pushover
40
2/3/2017
Nonlinear Pushover
41
2/3/2017
Roof
Displacement
Base
Reaction
42
2/3/2017
43
2/3/2017
44
2/3/2017
Knowing the
order of hinge
occurrence may
or may not be
important, but
knowing that
beam hinging is
likely to dominate
the nonlinear
behavior is
important
500
450
400
Base Force (kips)
350
300
250
200 4
150
100
50
0
0 10 20 30 40 50 60
Roof Displacement (inches)
45
2/3/2017
500
450
Prevention
Occupancy
400
Base Force (kips)
Immediate
Collapse
350
Damage
300
Control
Safety
Life
250
200
150
100
50
0
0 10 20 30 40 50 60
Roof Displacement (inches)
Nonlinear Pushover
46
2/3/2017
Nonlinear Pushover
47
2/3/2017
Frame 1
Frame 1
48
2/3/2017
Frame 1
Frame 1
49
2/3/2017
Frame 1
Frame 1
50
2/3/2017
Frame 2
Frame 2
51
2/3/2017
Frame 2
Frame 2
52
2/3/2017
Frame 2
Frame 3
53
2/3/2017
Frame 3
Frame 3
54
2/3/2017
Frame 3
Frame 3
55
2/3/2017
700
600
500
Base Shear (kips)
400 Frame 1
2
Frame 2
Frame 3
300
200
100
1
0
0 2 4 6 8 10 12 14
Rooftop Displacement (inches)
SEAU: 2017 Education Conference
A Caveat…
While this example corroborates the requirements of AISC 341, research by Sen,
Roeder, Lehman & Berman at University of Washington shows much better performance
in braced frames when quasistatically loaded under a fully reversed increasing amplitude
cyclical protocol…
See Structure
Magazine – July 2015
SEAU: 2017 Education Conference
56
2/3/2017
Nonlinear Pushover
Why?
1 3
What happens if the geometric nonlinear stiffness matrix overpowers the initial
matrix?
24EI/L3 6EI/L2 6EI/L2 -6P1/5L-6P3/5L -P1/10 -P3/10
6EI/L2 8EI/L 2EI/L + -P1/10 -2P1L/15-2P2L/15 -P2L/30
6EI/L2 2EI/L 8EI/L -P3/10 -P2L/30 -2P2L/15-2P3L/15
The effects of geometric nonlinearity can become greater than the effects of
strain hardening.
SEAU: 2017 Education Conference
57
2/3/2017
Pushover Analysis:
Roof Disp.
350
300
Base Shear
250
200
150
100
50
0
0 2 4 6 8 10
Reaction
P 250
θ2 2 θ3 200
F W/O P-delta
150
W/ P-delta
1 3 100
50 Geometric
0 nonlinearity
0 2 4 6 8 10 controls over strain
Rooftop Displacement hardening
58
2/3/2017
Base Shear
250
200
W/O P-delta
150
W/ P-delta
100
50
0
0 2 4 6 8 10
Rooftop Displacement
Base
Reaction
Nonlinear Pushover
59
2/3/2017
North-South Pushover
North-South Pushover
60
2/3/2017
North-South Pushover
North-South Pushover
61
2/3/2017
North-South Pushover
North-South Pushover
62
2/3/2017
North-South Pushover
Failure occurs at 38 inches rooftop displacement
Pushover Curve
1400
1200
1000
Base Shear (kip)
200 4
0
0 5 10 15 20 25 30 35 40
Displacement (in)
63
2/3/2017
Nonlinear Pushover
Nonlinear Pushover
64
2/3/2017
Nonlinear Pushover
Nonlinear Pushover
65
2/3/2017
Nonlinear Pushover
Brace Hysteretic Backbone
A typical brace hysteretic backbone: 600
500
400
300
0
-2 -1 0 1 2 3 4 5
-100
-200
-300
-400
Axial Deformation (in)
Nonlinear Pushover
300
200
100
0
-2 -100 0 2 4 6
-200
-300
-400
Axial Deformation (in)
Concept courtesy of
Brent Maxfield, SE
66
2/3/2017
Nonlinear Pushover
Nonlinear Pushover
67
2/3/2017
68
2/3/2017
80
80 Fmax=68.8 ksi
Fmax=62.8 ksi
60 Fye=55 ksi
60 Fye=51.7 ksi
40 40
TENSION TENSION
27.5 ksi
20
20
Stress (ksi)
Stress (ksi)
0
-3 -2 -1 0 1 2 3 0
-1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 3 3.5
-20
COMPRESSION
-40 COMPRESSION
Fye=-51.7 ksi
-40
-42.24 ksi
-60
Fmax=-65.2 ksi
-60
-80 Strain (%)
Strain (%)
tension Intermediate
Local buckling less critical
displacement
than inelastic buckling
typical
buckling
brace
unbonded
brace
compression
69
2/3/2017
70
2/3/2017
71
2/3/2017
72
2/3/2017
73
2/3/2017
Nonlinear Pushover
Bennett Federal Building
N-S Push-over V vs Roof Displacement
35000
30000
25000
20000
BRBF
15000
10000 BRBF Braces
5000 Begin to Yield
0
0 5 10 15
Roof Displacement - Inches
74
2/3/2017
Nonlinear Pushover
Horizontal Combinations
In the same direction, along the same line - Shear Wall and
Moment frame:
w14x120
20’ 20’
SEAU: 2017 Education Conference
75
2/3/2017
Horizontal Combinations
In the same direction, along the same line - Shear Wall and Heavy
Moment frame:
w24x370
12” thick wall
12’ f’c=4,000psi
I = 0.35Ig
7.69k 7.69k 84.62k
20’ 20’
SEAU: 2017 Education Conference
Horizontal Combinations
Unless a stiffness compatibility exists, logic may predicate that the less
rigid system be abandoned entirely.
76
2/3/2017
Horizontal Combinations
Dual System?
Horizontal Combinations
77
2/3/2017
Horizontal Combinations
Horizontal Combinations
78
2/3/2017
Horizontal Combinations
Horizontal Combinations
79
2/3/2017
Horizontal Combinations
Horizontal Combinations
80
2/3/2017
Horizontal Combinations
Horizontal Combinations
81
2/3/2017
Horizontal Combinations
Horizontal Combinations
Pushover Curve
600
500
Base Shear (kip)
400
300
200
100
0
0 1 2 3 4 5 6 7
Displacement (in)
SEAU: 2017 Education Conference
82
2/3/2017
Horizontal Combinations
Unless a stiffness compatibility exists, logic may predicate that the less
rigid system be abandoned entirely.
Horizontal Combinations
Another Approach?
83
2/3/2017
Horizontal Combinations
concrete
12” thick wall with
ρvert = 0.0015
reinforcement
Horizontal Combinations
Stress in
Reinforcement
Stress in
Concrete
84
2/3/2017
Horizontal Combinations
Stress in
Reinforcement
Stress in
Concrete
Horizontal Combinations
Stress in
Reinforcement
Stress in
Concrete
85
2/3/2017
Horizontal Combinations
Stress in
Reinforcement
Stress in
Concrete
Horizontal Combinations
Stress in
Reinforcement
Stress in
Concrete
86
2/3/2017
Horizontal Combinations
Pushover Curves
400
350
300
200
150
100
50
0
0 1 2 3 4 5 6 7
Displacement (in)
Nonlinear Pushover
87
2/3/2017
Nonlinear Pushover
Seismic Retrofit
Nonlinear Pushover
Seismic Retrofit
88
2/3/2017
89
2/3/2017
Nonlinear Pushover
Seismic Retrofit
90
2/3/2017
0.8
0.6
1
0.8
0.6
0.4
Acceleration(g)
Acceleration (g)
0.4
0.2
0.2 0
0 -0.2
-0.4
-0.2
-0.6
-0.4
-0.8
-0.6 -1
0 10 20 30 40 50 60 0 10 20 30 40 50 60 70 80
Time (sec) Time (sec)
91
2/3/2017
15000
10000
Beam Force (kip-in)
5000
0
-0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05
-5000
-10000
-15000
-20000
Beam Rotation (rad)
92
2/3/2017
Keeping it Simple?
Successful & efficient
nonlinear modeling
requires a targeted
approach.
93
2/3/2017
• Time
• $$
• Ground Motions
• Owner Objectives
• Computing Power
• Review
Software
• Hand 1.0 (by Ron Hamburger)
• SAP 2000 Ultimate
• ETABS Ultimate
• Perform 3D
• STAAD Pro
• Opensees (Open System for Earthquake Engineering Simulation,
NEES)
• ANSYS
• LS Dyna
• ANSR (proprietary)
• RISA?
94
2/3/2017
The End
95