b1461930.006592.001
b1461930.006592.001
b1461930.006592.001
by
Umut Utku CELEP
B.S., Civil Engineering, Dokuz Eylül University, 1998
M.S., Earthquake Engineering, Boğaziçi University, 2001
ACKNOWLEDGEMENTS
I would like to express my deepest gratitude to Professor Nuray Aydınoğlu for his
invaluable guidance and help during the preparation of this dissertation, who had been a
father to me as well as a wise mentor when I need. I would like to thank for his patience
and understanding for the times I had to split focus to my job as a practicing engineer.
I would like to express my sincere gratefulness to my dear wife Beyza Celep for her
support and encouragement. This dissertation would not have been completed without her
presence.
I would also like to mention my special gratitude to my boss and mentor Osman
Müyesser for his support and understanding.
Finally, special thanks to my old friend and colleague Göktürk Önem for his valuable
discussions and feedback.
v
ABSTRACT
Previous research indicates that, shear force demand in yielding walls are not
proportional to the design moments calculated by code procedures and higher shear force
demands develop along the wall with respect to code predictions as a consequence of the
higher mode effects after the plastic hinge formation at the base of the wall. 2007 version
of the Turkish Seismic Design Code takes the dynamic shear amplification phenomena into
account with a constant base shear amplification factor of 1.5 regardless of the first mode
period and ductility level of the wall. However, results obtained from extensive non-linear
time history analyses performed on generic walls in this study indicate that dynamic shear
amplifications increase with increasing first mode period, Strength Reduction Factor (R)
and ground motion intensity. A dynamic base shear amplification relationship as a function
of the first mode period and strength reduction factor has been proposed for the Turkish
Seismic Design Code (2007), based on the regression analysis of the non-linear time
history analysis. A story shear force profile has been suggested for the Turkish Seismic
Design Code (2007), which is intended for not only preventing shear failures at the base
but also along the height of the wall. As a side product of the non-linear time history
analyses, a moment profile has also been proposed for use in the Turkish Seismic Design
Code (2007). A modal decomposition technique is presented in this study for
demonstrating the effects of the higher modes on the dynamic shear amplification
phenomenon.
vi
ÖZET
Geçmiş çalışmalar, perde tabanında plastik kesit oluşması sonrasında ortaya çıkan
yüksek mod etkileri nedeniyle kesme kuvveti talebinin yönetmelik yöntemleri ile elde
edilen tasarım kesme kuvvetlerinden büyük olduğunu ve tasarım momentleri ile uyumsuz
olduğunu göstermiştir. 2007 Türk Deprem Yönetmeliği dinamik kesme kuvveti
büyütmesini periyot ve deprem yükü azaltma katsayısından (R) bağımsız olarak tasarım
taban kesme kuvvetlerine uygulanan 1.5 katsayısı ile dikkate almaktadır. Fakat temsili
perde duvarlar üzerinde yapılan zaman tanım alanında doğrusal olmayan analizlerinin
sonucunda, dinamik kesme kuvveti büyütmelerinin artan birinci mod periyodu ve deprem
yükü azaltma katsayısı ile arttığını göstermektedir. Zaman tanım alanında doğrusal
olmayan analizlerden elde edilen sonuçlara yapılan regresyon analizlerinden; 2007 Türk
Deprem Yönetmeliği’nde kullanılabilecek, periyot ve deprem yükü azaltma katsayısının
fonksiyonu olan bir dinamik taban kesme kuvveti büyütmesi ilişkisi önerilmiştir. Bunun
yanında, 2007 Türk Deprem Yönetmeliği’nde kullanılabilecek perde yüksekliği boyunca
kesme göçmelerini de engellemek amacıyla bir kat kesme kuvveti dağılımı önerilmiştir.
Yapılan zaman tanım alanında doğrusal olmayan analizlerin bir yan ürünü olarak, 2007
Türk Deprem Yönetmeliği için bir perde moment profili önerilmiştir. Çalışmada, yüksek
modların dinamik kesme kuvveti büyütmesi üzerindeki etkisi bir modal ayrıştırma tekniği
ile gösterilmiştir.
vii
TABLE OF CONTENTS
ACKNOWLEDGEMENTS..................................................................................................iv
ABSTRACT...........................................................................................................................v
ÖZET ....................................................................................................................................vi
LIST OF FIGURES ..............................................................................................................ix
LIST OF TABLES..............................................................................................................xiii
LIST OF SYMBOLS ..........................................................................................................xiv
1. INTRODUCTION .............................................................................................................1
1.1. Motivation of the Study .............................................................................................2
1.2. Scope of Work ...........................................................................................................2
2. LITERATURE REVIEW ..................................................................................................4
3. DESIGN OF THE GENERIC STRUCTURAL WALLS................................................11
3.1. Pre-Design of Generic Structural Walls...................................................................12
3.2. Design of Generic Structural Walls .........................................................................14
4. A NOVEL MODAL DECOMPOSITION TECHNIQUE FOR IDENTIFYING
HIGHER MODE EFFECTS IN DYNAMIC SHEAR AMPLIFICATION....................17
4.1. Linear and Non-Linear Modes of Vibration ............................................................18
4.2. Modal Decomposition Technique ............................................................................23
4.2.1. Modal Expressions in The Linear Response ..................................................23
4.2.2. Incremental Form of Modal Expressions in The Non-Linear Response........24
5. NON-LINEAR TIME HISTORY ANALYSES OF GENERIC WALLS.......................37
5.1. Group 1 Analyses and Results .................................................................................39
5.2. Group 2 Analyses and Results .................................................................................45
5.3. Group 3 Analyses and Results .................................................................................49
5.4. Group 4 Analyses and Results .................................................................................53
5.5. Remarks on Results Obtained from Non-Linear Time History Analyses ...............57
6. PROPOSED DYNAMIC SHEAR AMPLIFICATION FACTOR
RELATIONSHIP AND STORY SHEAR PROFILE.....................................................58
6.1. Comparison of the Proposed Amplification Relationship with the Other
Relationships............................................................................................................60
6.2. Proposed Story Shear Force Profile Along the Wall Height ...................................65
viii
LIST OF FIGURES
Figure 2.1. Inertia forces before and after yielding at the base of the wall ......................4
Figure 2.2. Base shear amplification factors suggested by Derecho and Corley
(1984) .............................................................................................................6
Figure 2.3. Base shear amplification factors suggested by Seneviratna and
Krawinkler (1994)..........................................................................................8
Figure 3.1. Illustration of generic walls considered in this study...................................11
Figure 3.2. Design acceleration response spectrum used in the study ...........................14
Figure 3.3. Modelling representation of structural walls................................................15
Figure 4.1. Linear and non-linear 1st mode shapes of the 16 story generic wall ...........19
Figure 4.2. Linear and non-linear 2nd mode shapes of the 16 story generic wall..........20
Figure 4.3. Linear and non-linear 3rd mode shapes of the 16 story generic wall ..........20
Figure 4.4. Identification of linear and non-linear phases of the wall through plastic
hinge rotations at the base ............................................................................25
Figure 4.5. Plastic rotation history at the base of 16 story generic wall (R=6)
obtained from non-linear time history analysis (Whittier Narrows, 1987-
Brea Dam Record)........................................................................................28
Figure 4.6. Modal and superposed base shear history of the 16 story generic wall
designed with R=6........................................................................................29
Figure 4.7. Agreement between the base shear of the 16 story generic wall obtained
from mode superposition and Group 2 non-linear time history analysis .....30
Figure 4.8. Agreement between the overturning moment of the 16 story generic wall
obtained from mode superposition and Group 2 non-linear time history
analysis.........................................................................................................30
Figure 4.9. Increment of modal seismic forces of the 16 story generic wall designed
with R=6 (Immediately before the peak response instant)...........................31
Figure 4.10. Modal seismic forces of the 16 story generic wall designed with R=6 (at
the peak response instant) ............................................................................32
Figure 4.11. Modal and superposed story shear forces of the 16 story generic wall
designed with R=6 (at the peak response instant)........................................33
x
Figure 4.12. Increment of modal moments of the 16 story generic wall designed with
R=6 (immediately before the peak response instant)...................................34
Figure 4.13. First ten seconds of the modal and superposed dynamic base shear
amplification factor history of the 16 story generic wall designed with
R=6...............................................................................................................35
Figure 5.1. Illustration of damping ratio assumption in time history analyses ..............38
Figure 5.2. Illustration of analytical model for Group 1 analyses..................................39
Figure 5.3. Story moment diagrams of 16 story wall with strength level R=6 (Group
1 time history analyses)................................................................................40
Figure 5.4. Story shear force diagrams of 16 story wall with strength level R=6
(Group 1 time history analyses) ...................................................................41
Figure 5.5. Mean story moment diagram of 16 story wall with strength level R=6
(Group 1 time history analyses) ...................................................................42
Figure 5.6. Mean story shear force diagram of 16 story wall with strength level R=6
(Group 1 time history analyses) ...................................................................42
Figure 5.7. Mean story shear force amplification factors computed for the 16 story
wall with strength level R=6 (Group 1 time history analyses).....................43
Figure 5.8. Mean base shear amplification factors obtained from group 1 time history
analyses ........................................................................................................44
Figure 5.9. Illustration of analytical model for Group 2 analyses..................................45
Figure 5.10. Story moment diagrams of 16 story wall with strength level R=6 (Group
2 time history analyses)................................................................................46
Figure 5.11. Story shear force diagrams of 16 story wall with strength level R=6
(Group 2 time history analyses) ...................................................................46
Figure 5.12. Mean story moment diagram of 16 story wall with strength level R=6
(Group 2 time history analyses) ...................................................................47
Figure 5.13. Mean story moment diagram of 16 story wall with strength level R=6
(Group 2 time history analyses) ...................................................................47
Figure 5.14. Mean story shear force amplification factors computed for the 16 story
wall with strength level R=6 (Group 2 time history analyses).....................48
Figure 5.15. Mean base shear amplification factors obtained from Group 3 time
history analyses ............................................................................................48
Figure 5.16. Illustration of analytical model for Group 3 analyses..................................49
xi
Figure 5.17. Story moment diagrams of 16 story wall with strength level R=6 (Group
3 time history analyses)................................................................................50
Figure 5.18. Story shear force diagrams of 16 story wall with strength level R=6
(Group 3 time history analyses) ...................................................................50
Figure 5.19. Mean story moment diagram of 16 story wall with strength level R=6
(Group 3 time history analyses) ...................................................................51
Figure 5.20. Mean story shear force diagram of 16 story wall with strength level R=6
(Group 3 time history analyses) ...................................................................51
Figure 5.21. Mean story shear force amplification factors computed for the 16 story
wall with strength level R=6 (Group 3 time history analyses).....................52
Figure 5.22. Mean base shear amplification factors obtained from Group 3 time
history analyses ............................................................................................52
Figure 5.23. Illustration of analytical model for Group 3 analyses..................................53
Figure 5.24. Story moment diagrams of 16 story wall with strength level R=6 (Group
4 time history analyses)................................................................................54
Figure 5.25. Story shear force diagrams of 16 story wall with strength level R=6
(Group 4 time history analyses) ...................................................................54
Figure 5.26. Mean story moment diagram of 16 story wall with strength level R=6
(Group 4 time history analyses) ...................................................................55
Figure 5.27. Mean story shear force diagram of 16 story wall with strength level R=6
(Group 4 time history analyses) ...................................................................55
Figure 5.28. Mean story shear force amplification factors computed for the 16 story
wall with strength level R=6 (Group 4 time history analyses).....................56
Figure 5.29. Mean base shear amplification factors obtained from Group 4 time
history analyses ............................................................................................56
Figure 6.1. Agreement between the proposed relationship and mean base shear
amplification factors.....................................................................................59
Figure 6.2. Base shear amplification factors computed by using Ghosh and
Markevicius (1990) relationship ..................................................................61
Figure 6.3. Ratio of the base shear amplification factors between the relationship
suggested by Ghosh and Markevicius (1990) and this study.......................61
Figure 6.4. Base shear amplification factors computed by using Eibl and Keintzel
(1988) relationship .......................................................................................62
xii
Figure 6.5. Ratio of the base shear amplification factors between the relationship
suggested by Eibl and Keintzel (1988) and this study .................................63
Figure 6.6. Base shear amplification factors computed by using Rutenberg and Nsieri
(2006) relationship .......................................................................................64
Figure 6.7. Ratio of the base shear amplification factors between the relationship
suggested by Rutenberg and Nsieri (2006) and this study...........................64
Figure 6.8. Observed trend of the shear force profile obtained from non-linear time
history analyses of generic walls..................................................................66
Figure 6.9. Proposed story shear force profile................................................................66
Figure 6.10. Eurocode 8.1 (CEN, 2004) story shear force profile ...................................67
Figure 6.11. Shear force profile proposed by Rutenberg and Nsieri (2006) ....................68
Figure 6.12. Proposed moment profile .............................................................................72
Figure 6.13. Agreement between proposed moment profile and non-linear time history
analysis of the generic walls designed with R=2 .........................................73
Figure 6.14. Agreement between proposed moment profile and non-linear time history
analysis of the generic walls designed with R=4 .........................................74
Figure 6.15. Agreement between proposed moment profile and non-linear time history
analysis of the generic walls designed with R=6 .........................................74
Figure 6.16. Agreement Between proposed moment profile and non-linear time
history analysis of the generic walls designed with R=8 .............................75
Figure 6.17. Agreement between proposed moment profile and non-linear time history
analysis of the generic walls designed with R=10 .......................................75
xiii
LIST OF TABLES
LIST OF SYMBOLS
1. INTRODUCTION
Structural walls systems constitute a noteworthy share in high rise building stock
built in seismic zones. Structural walls offer advantages not only to structural engineers
such as their considerable lateral stiffness and strength and consequent reduction of seismic
deformation demands; but also to architects by whom the structural system can also be
used as partitioning elements with reduced disturbance.
On the other hand structural walls exhibit the disadvantage of redundancy problems
under seismic attack unless they are designed to respond fully elastic, which would be
impractical and uneconomical if not impossible. Thus, general design trend is to
incorporate Strength Reduction Factors (hereinafter will be called as R factors) to reduce
seismic demands on the walls by taking the risk of an acceptable damage that is deemed to
develop under design ground motion levels. Traditionally strength reduction factors
uniformly apply to all design seismic effects including wall moments and shear forces, that
is, wall design shear forces are assumed to be proportional to design moments. However,
previous nonlinear time history analyses by various researchers revealed that, this is not the
case for wall structures designed to respond the seismic attack in the nonlinear range and
shear force demands are greater than the demands expected in the usual design process.
This amplification with respect to design shear forces have been attributed to the higher
mode effects after yielding occurring at the base of the wall; however quantified
explanation of this phenomenon and the range of amplification is still an open field of
research.
Current rationale of most seismic design codes is to inhibit brittle shear failure modes
while attaining the flexural capacity at the base of the wall, which stems from capacity
design principles. However, amplified shear forces, particularly at the base of the wall,
may cause unexpected premature brittle shear failure modes of diagonal tension, diagonal
compression and sliding before the plastic hinge development if required measures are not
taken. This may cause destructive consequences for the structural wall systems designed
for high ductility demands since adequate reinforcement and detailing might not be
provided for an extra shear force demand.
2
Despite the fact that the effect of higher modes are held responsible for shear force
amplification after the plastic hinge formation at the base of the wall, analytical
demonstration of the problem is rather limited in the literature.
This work comprises a parametric study performed on generic walls for deriving a
dynamic base shear amplification factor relationship, which could be proposed for in lieu
of the constant base shear amplification in the Turkish Seismic Design Code (2007). A
modal decomposition technique has been presented to demonstrate the responsibility of
higher modes in dynamic shear amplifications.
2. LITERATURE REVIEW
LOWERED CENTROID
VE1
VE2
heff
heff
My My
VE2 > VE1
INERTIA FORCE DISTRIBUTION AT INERTIA FORCE DISTRIBUTION
FIRST MODE RESPONSE AT THE ONSET AFTER YIELDING AT THE BASE
OF YIELDING AT THE BASE
Figure 2.1. Inertia forces before and after yielding at the base of the wall
derived from conventional code procedures ( Veb ) to obtain the expected base shear demand
Dynamic base shear amplification has not been recognized until the work of
Blakeley et al. (1975), who pioneered the idea that higher mode effects amplify the base
shear with respect to the base shear obtained from conventional code procedures after the
formation of a plastic hinge at the base of the wall. They stated that after a plastic hinge
formation at the base of the wall, the center of lateral inertial loading becomes lower or
higher with respect to the center predicted by the code distribution, effectively causing
higher shear forces than predicted by the code. Blakely et al. (1975) performed non-linear
time history analyses on various structural wall models and recommended base shear
amplification factors as a function of number of stories (n) of the wall structure. Based on
these recommendations the provisions of New Zealand seismic code NZS3101 (Standards
New Zealand 1982) have changed in 1982 and formed the basis of the New Zealand code
perspective for the base shear amplification factors as given in Equation (2.2) which is still
in use (Standards New Zealand 2006).
n
ωn = 0.9 + n≤6
10
(2.2)
n
ωn = 1.3 + ≤ 1.8 n>6
30
It is seen from Equation (2.2) that, New Zealand perspective on the base shear
amplification is dependent on the first mode period (fundamental period) of the wall,
which is implicitly related to the number of stories of the wall (n).
A parametric investigation study carried out by Derecho and Corley (1984) revealed
the increasing dependency of base shear amplification not only to the first mode period but
also to the ductility level of the wall. The authors presented their base shear amplification
results in the form of a chart format as shown in Figure 2.2. Derecho and Corley also stated
6
that there is a linear relationship between the base shear amplification factor and intensity
of the ground motion.
Figure 2.2. Base shear amplification factors suggested by Derecho and Corley (1984)
in which Vn is the base shear capacity for the structure calculated by ultimate limit state
analysis assuming a triangular force distribution, c represents a coefficient that stands for
the higher mode effects varying between 0.25 and 0.30, W is the seismic weight of the wall
and PGA is the design peak ground acceleration. Kabeyasawa concluded that maximum
dynamic shear in a wall under seismic attack is higher than the statically computed shear
7
due to the higher mode effects and the additional shear force demand in the non-linear
range is a function of the peak ground acceleration.
A study by Eibl and Keintzel (1988) has also shown that base shear amplification is
both functions of the first mode period and expected ductility level. The authors have
proposed a base shear amplification factor relationship as given in Equation (2.4), which
later formed the basis of the Eurocode 8 EN 1998-1 (CEN, 2004) provisions for base shear
amplifications of wall systems designed for high ductility. Eibl and Keintzel have stated in
their study that the first term under the square root denotes the contribution the first mode
that is reduced by the strength reduction factor and the second term implicitly denotes the
importance of the second mode in yielding systems.
2 2
⎛γ M ⎞ ⎛ S (T ) ⎞
ε = q ⎜ rd rd ⎟ + 0.1⎜ e C ⎟ ≤ q (2.4)
⎝ q M Ed ⎠ ⎝ Se (T1 ) ⎠
in which ε is the base shear amplification factor, q is the strength reduction factor, γRd the
overstrength factor, MRd the design flexural strength at the wall base, MEd the design
moment at the wall base, TC the upper limit period of the spectral acceleration plateau in
the design spectrum and Se(T) is the ordinate of the elastic response spectrum.
Ghosh and Markevicius (1990) has concluded in their parametric study that no
dependence with the first mode period could be established and claimed that dynamic
shear amplification is likely to be a function of ground motion intensity. The authors has
suggested a relationship for the maximum base shear demand after yielding at the base of
the wall, in which the first term is a function of the intensity of the ground motion and the
second term is the base shear corresponding to plastic hinge formation at the base wall as
presented in Equation (2.5).
PGA My
Vmax = 0.25 × W × + (2.5)
g 0.67 H
where W is the seismic weight of the wall, PGA the design peak ground acceleration, My
the yield moment capacity at the base of the wall and H is the total wall height.
8
Seneviratna and Krawinkler (1994) have emphasized the dependency of the shear
amplification of structural walls on both strength level and the first mode period of the
structure as shown in Figure 2.3. The authors stated that shear amplifications increase for
longer periods and lower strength levels, which can be attributed to the higher mode
vibrations, particularly the second mode, of the elastic portion of wall above the plastic
hinge at the base. They also compared their results with the relationship proposed by
Ghosh and Markevicius (1990) and concluded that Ghosh and Markevicius relationship
predicts the base shear amplifications in good agreement with their results for medium-
high walls with relatively low strength levels, however underestimates amplifications for
tall structures due to the lack of a term related to the period of the structure.
Figure 2.3. Base shear amplification factors suggested by Seneviratna and Krawinkler
(1994)
9
Filiatrault et al. (1994) have studied the problem from the Canadian code approach
and concluded that brittle shear failures are possible for ductile walls designed by the
Canadian seismic code provisions. Based on a parametric study of walls designed
according the Canadian code, they suggested seismic zone dependent strength reduction
factors of 1.0 or 1.5 for shear design to account for base shear amplifications.
Rutenberg and Nsieri (2006) have studied the problem through Eurocode 8
perspective and claimed that current provisions of Eurocode 8 are in need of revision since
they do not predict the shear amplifications realistically for the medium and high ductility
classes of structural walls. Based on a parametric study by using non-linear time history
analyses, the authors recommended the relationship given in Equation (2.6) in lieu of the
current Eurocode relationship for dynamic base shear amplifications.
in which q is the strength reduction factor and T is first mode period of the structure.
Rutenberg and Nsieri (2006) concluded that modification is necessary to Eurocode 8 base
shear amplification relationship, which should increase with increasing period and strength
reduction factors and should take into account the distribution of plasticity along the wall
height.
Krawinkler (2006) expressed that “Once a plastic hinge forms at the base of a wall
structure, the dynamic response characteristics of the structure above the base change
radically. Different ‘modes’ become predominant and the effective lever arm (M/V ratio)
becomes smaller, and for relatively tall and slender structures it becomes much smaller.
Correspondingly, the dynamic shear force that can be generated may become much larger
than that indicated by a static M/V concept in which the resultant lateral force is applied
approximately at 2/3 of the structure height (or higher).” The author emphasized that
despite the detrimental reality that base shear demand is higher than predicted by
traditional code provisions, higher actual shear strength of walls than recommended by
code equations may be compensating this negative effect. The author also pointed out a
very critical issue of wall behavior based on the recent test results, where wall response is
10
often controlled by sliding at the base and shear force distribution may be totally different
than predicted by analytical methods.
It is interesting to note that USA has been rather hesitant to apply dynamic shear
amplification factors in the national seismic design codes, despite several researchers in the
USA have stressed the fact on shear amplifications in their studies. As an exception,
Commentary to the Recommended Lateral Force Requirements of the Structural Engineers
Association of California (SEAOC, 1999) recommended adopting the amplification factors
advocated by the New Zealand code.
Turkish Seismic Design Code (1997) did not take into account the dynamic base
shear amplifications in structural walls at any level of ductility, and the design for shear
was based on the shear forces corresponding to the lateral force distribution used for
flexural design. Present version of the code (2007) takes the base shear amplification
phenomena into account with a constant base shear amplification factor of 1.5 regardless of
the first mode period and ductility level of the wall.
Seismic design codes could prefer utilizing base shear amplification factors
suggested by researchers, as long as they are practical to implement in design practice.
However they have a serious drawback of being related to the intensity of ground motion
records they are derived from. Research on shear amplification of structural walls is still
continuing and there is limited global consensus between the researchers (also between the
codes) on the degree of base shear amplifications, despite the fact that the source of this
amplification is recognized as the higher mode effects initiating after the plastic hinge
formation at the base of the wall.
11
The numerical and statistical findings in this study are based on the results obtained
from the analysis of generic structural walls that are intended to represent the structural
wall systems in the building stock. Thus, generic structural walls with 8, 12, 16, 20 and 30
stories and typical story heights of 3.0m have been considered as illustrated in Figure 3.1,
in order to cover a broad range of wall characteristics in terms of flexibility and strength.
m
m
m
m
m
m
m
m
m
m
m m
m m
m m
m m
30 @ 3m = 90 m
m m m
m m m
m m m
m m m
m m m m
20 @ 3m = 60 m
m m m m
m m m m
16 @ 3m = 48 m
m m m m
m m m m m
12 @ 3m = 36 m
m m m m m
m m m m m
8 @ 3m = 24 m
m m m m m
m m m m m
m m m m m
m m m m m
m m m m m
where possible so as not to deviate from realistic results. Details of the pre-design and
design steps are presented in the following sections of this chapter.
Analysis and design effects obtained herein have been computed by using SAP 2000
v11 (CSI, 2006) structural analysis software licensed to the Department of Earthquake
Engineering of Boğaziçi University.
A simple approach has been utilized to determine the section dimensions of the
generic walls. It has been assumed that the first mode period corresponding to gross
(uncracked) stiffness shall approximately satisfy an old empirical relation correlated with
the total building height (H) as given in Equation (3.1) (Turkish Seismic Design Code,
1997).
3
T1 = 0.05 H w 4 (3.1)
Thus, by assigning building heights to Equation (3.1), first mode periods have been
determined for each generic wall in concern. Assuming that the wall mass and stiffness are
uniformly distributed along the wall height, the natural vibration frequency of the system
can be estimated by Equation (3.2) (Chopra 2006);
3.516 EI g
ω1 = (3.2)
H w2 m( y )
where Ig is the moment of inertia of the gross section properties and m(y) is the mass per
unit length of the wall.
2π
T1 = (3.3)
ω1
13
3 2π
0.05 H w 4 =
3.516 EI g (3.4)
H w2 m( y )
Solution of Equation (3.4) for the generic walls resulted in the following first mode
period, section dimensions and typical story masses, as given in Table 3.1.
Design basis acceleration response spectrum given in Turkish Seismic Design Code
(2007) corresponding to 10 percent probability of exceedence in 50 years has been used for
the representation of ground motion (Figure 3.2). The design acceleration response
spectrum has been assumed to be based on local site conditions with site class Z3 (corner
periods of TA = 0.15sec and TB = 0.60 sec) and Seismic Zone 1 (Effective peak ground
acceleration of 0.40g), respectively.
12
10
8
A (m/s )
2
0
0.00 0.50 1.00 1.50 2.00 2.50 3.00
Period (sec)
Generic structural walls have been modeled with frame (beam) elements in terms of
stiffness and lumped masses at each story level as shown in Figure 3.3. Flexural stiffness
of the walls has been reduced by a factor 0.50 to obtain cracked section stiffness,
approximating the cracking and softening of linear stiffness under seismic attack.
The boundary conditions at the base of the walls have been assumed to be fixed, thus
rotations arising from the foundation or soil stiffness have been neglected.
15
Response spectrum analyses have been performed for each generic wall in order to
obtain the design story moments and shear forces along the wall height. All modes of the
structure exhibiting modal amplitudes in the lateral direction, that is equal to the number of
stories of the wall in concern, have been taken into account during response spectrum
analysis, and critical damping ratio has been assumed as 5 percent.
Each generic wall has been analyzed by taking into account Strength Reduction
Factors R = 2, 4, 6, 8 and 10 in order to provide a strength range. These Strength Reduction
Factors have been incorporated into response spectrum analyses by reducing the design
spectrum as per the requirements of the Turkish Seismic Design Code 2007. It can be
argued that realistic wall buildings in practice are not designed for strength reduction
factors of 8 or 10, which would result in drastically under-designed structures. However, it
is intended that high values of R may have some merit in cases of seismic assessment of
old structural walls designed with lesser seismic demands and weaker capacities.
Modal response quantities derived from response spectrum analyses have been
combined with the CQC modal combination method to obtain the design effects, namely
the story moments and shear forces, which are presented in Appendix A for each generic
wall with their corresponding strength reduction factors.
16
Design moments obtained at the base of the wall from response spectrum analyses
with the corresponding strength reduction factor have been hypothetically assumed as
equal to the design moment capacity or yield moment strength of the section at the base of
the generic walls. Thus, story moment and shear force profiles along the wall height have
also been supposed to be equal to the profiles obtained from the response spectrum
analyses as the yield moment is attained at the base.
17
As with systems exhibiting cantilever like behavior, structural wall systems are
statically unstable but dynamically stable structures when flexural yielding occurs at their
bases. During seismic excitation, as the system alters from linear to nonlinear behavior
mode, a plastic hinge or a series of hinges develop at the base of the wall and the structural
wall loses its redundancy for a time interval until it is stabilized back by inertia forces.
During this time interval, a new system with a completely unique dynamic characteristics
prevail until the system velocity becomes zero, and the linear phase initiates with
unloading. Status of the structural wall before and after the plastic development at the base
of the wall shall be called hereinafter as “the linear phase” and “the non-linear phase”,
respectively, for the ease of understanding.
It has been observed in this study and by various researchers described in Chapter 2
that greater shear forces are obtained from nonlinear time history analyses than the shear
forces expected from design. It is a general consensus that higher mode effects in the non-
linear phases causes this amplification of shear forces with respect to the ones obtained
from linear design, however studies demonstrating individual contribution of higher mode
effects to dynamic amplification are very limited in the literature. Kabeyasawa (1987) has
called the response of second mode effects as “residual fluctuating forces”. Based on a
triangular first mode shape, he predicted the second mode shape from orthogonality
condition with the first mode and extracted the second mode response by eliminating the
triangular first mode effects from the non-linear time history analyses.
the technique. The technique has been demonstrated on the non-linear time history results
of 16 story generic wall designed with a Strength Reduction Factor of R=6 which have
been obtained from the Whittier Narrows 1987 - Brea Dam (Downstream) record in Group
2 analysis as explained in Chapter 5. (single plastic hinge at the base of the wall with
elasto-plastic moment-curvature hysteretic characteristics-see Section 5.2).
Modal characteristics of the wall in the linear and non–linear phases should be
determined prior to the implementation of the modal decomposition technique. Vibration
modes before the plastic hinge development will be denoted as the linear modes, whereas
vibration modes in the non-linear phases of the plastic hinge will be referred to as the non-
linear modes. Eigenvalues and mode shapes both for linear and non-linear phases of the
walls have been obtained by using a built-in feature of MATLAB. This feature utilizes the
Jacobi Method (Bathe, 1996) that converge accurate eigenvalues regardless of the stiffness
matrix conditions, as is the case for the wall systems in the non-linear phase.
General expression for the dynamic eigenvalue problem can be derived from the
equation of motion of an undamped multi degree of freedom system having a mass matrix
of M and a stiffness matrix of K as:
KΦ n = ωn2 MΦ n (4.1)
where, Φn and ωn2 denote the eigenvector and eigenvalue of the n’th mode of the system,
respectively.
Solution of Equation (4.1) is straightforward for the linear phase of the wall, however
the system is theoretically unstable in the nonlinear phase due the plastic hinge formation
and thus eigenvalues can not be obtained through classical methods.
It is interesting to note that the Jacobi Method is able to provide the eigenvalues of
the system even in the non-linear phase since the solution of the eigenvalue problem by the
Jacobi Method does not contain any division by any number. Hence, the method always
19
converges and yields an accurate solution for positive, zero or negative eigenvalues.
(Bathe, 1996; Wilson, 2006)
Figure 4.1 to Figure 4.3 shows the first three of the computed 16 linear and non-
linear mode shapes (eigenvectors) of the 16 story generic wall, respectively.
12
10
STORY
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
MODAL AMPLITUDE
Figure 4.1. Linear and non-Linear 1st mode shapes of the 16 story generic wall
20
12
10
STORY
0
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8
MODAL AMPLITUDE
Figure 4.2. Linear and non-linear 2nd mode shapes of the 16 story generic wall
14
12
10
STORY
2 LIN. PHASE
NONLIN. PHASE
0
-0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
MODAL AMPLITUDE
Figure 4.3. Linear and non-linear 3rd mode shapes of the 16 story generic wall
21
Figure 4.1 reveals a transformation of the first mode shape pattern from parabolic to
straight line as a consequence of the plastic hinge development at the base of the wall in
the non-linear phase. Period of vibration in the first mode has been calculated as 1.30 sec
in the linear phase, whereas in the non-linear phase it is infinite (ω1=0). The zero natural
frequency indicates that modal stiffness for the first mode becomes effectively zero after
yielding takes place provided that no hardening is assumed for the moment-curvature
relationship at plastic hinge forming at the base. Thus, non-linear first mode can be called
as the unstable mode of the system, which would be stabilized by the inertia forces during
seismic excitation.
However, the non-linear first mode of the system has a definite Modal Participation
Factor, since eigenvectors are definite
Φ1*MTx
Γ1* = (4.2)
Φ1*T MΦ1*
It is also interesting to note that, non-linear first mode of the system satisfy the
orthogonality conditions with the non-linear second mode as reported by Kabeyasawa
(1987):
Modification of the second and third mode shapes in the non-linear phase has been
observed as a shift of the bulk of the mode shape downwards as shown in Figure 4.2 and
Figure 4.3, due to the hinge development at the base of the wall. Similar behavior has been
observed for higher modes of the wall above the third non-linear mode. In order to
demonstrate the eigenvectors of the other generic walls buildings with 8, 12, 16, 20 and 30
stories, first two mode shapes, both in the linear and non-linear phases, are presented in
Appendix B.
22
Non-zero periods of vibration have been computed for the higher modes as shown in
Table 4.1. Those modes can be called as stable modes regardless of the plastic hinge
formation at the base of the wall.
Table 4.1. Linear and non-linear higher mode periods of the 16 story generic wall
LINEAR NON-LINEAR
PHASE PHASE
T (sec) T (sec)
MODE 2 0.2145 0.3064
MODE 3 0.0767 0.0947
MODE 4 0.0392 0.0454
MODE 5 0.0237 0.0266
MODE 6 0.0159 0.0174
MODE 7 0.0114 0.0123
MODE 8 0.0086 0.0092
MODE 9 0.0067 0.0071
MODE 10 0.0054 0.0057
MODE 11 0.0045 0.0047
MODE 12 0.0038 0.0039
MODE 13 0.0033 0.0034
MODE 14 0.0030 0.0030
MODE 15 0.0027 0.0028
MODE 16 0.0026 0.0026
23
f m (t ) = M Φ m Γ m am (t ) (4.4)
ΦTm M 1
Γm = (4.5)
ΦTm M Φ m
Φm = Φm Γm
(4.6)
Φ = [Φ1 Φ 2 ..... Φ m ..... Φ n ]
where n refers to total number of degrees of freedom in the direction of seismic excitation.
Effective mass matrix can be defined as:
M = MΦ (4.7)
f (t ) = M a(t ) (4.8)
T
a(t ) = a1 (t ) a2 (t ) .... am (t ) .... an (t ) (4.9)
Equation (4.4) can be written in incremental form at the (i)’th step representing the
response increment in between the time stations ti and ti+1 :
m Γ m ∆am
∆f m(i) = M Φ(i) (i) (i)
(4.10)
prevailing in the (i)’th step, depending on whether the base plastic hinge response is along
loading/unloading lines or on the yielding plateau of the response hysteresis. Now ∆am (t )
refers to m’th mode pseudo-acceleration response increment in the (i)’th step.
where,
Thus, in nonlinear response, two different n by n effective mass matrices, M (i) , are
defined, one each corresponding to linear and non-linear modes, respectively.
From plastic rotation time-history response at the base of the wall, the time steps
associated with loading/unloading and the non-linear phases of the hysteretic response of
the base plastic hinge can be detected as illustrated in Figure 4.4.
ZERO ROTATIONAL
MOMENT VELOCITY
(LINEAR MODES)
LINEAR PHASE
PLASTIC
ROTATION (rad)
ZERO ROTATIONAL
VELOCITY NON-LINEAR PHASE
(NON-LINEAR MODES)
Figure 4.4. Identification of linear and non-linear phases of the wall through plastic hinge
rotations at the base
Thus modal increments of seismic forces, ∆f m(i) , can be readily obtained from
Equation (4.10). Further, base shear and base bending moment modal increments can be
obtained as:
Adding to those at the previous time station ti, m’th mode base shear and base
bending moment at time station ti+1 can be obtained as:
V b (ti+1 ) = V1b (ti+1 ) + V2b (ti+1 ) + .... + Vmb (ti+1 ) + .... + Vnb (ti+1 ) (4.21)
M b (ti+1 ) = M 1b (ti+1 ) + M 2b (ti+1 ) + .... + M mb (ti+1 ) + .... + M nb (ti+1 ) (4.22)
which should be exactly equal to those directly obtained from the nonlinear response
history analysis of the wall.
( )
2
m Γ m ∆d m ω m
∆f m(i) = M Φ(i) (i) (i) (i)
(4.23)
Increments of modal seismic forces for the first non-linear mode should be obtained
as zero from Equation (4.23) since the frequency associated with the first non-linear mode
is obtained as zero provided that there is no strain hardening at the base plastic hinge.
Thus, increments of modal base shear and modal overturning moment associated with the
first non-linear mode should also be zero from Equations (4.17) and (4.18). This finding
indicates that first mode of the wall in the non-linear phase do not contribute to any force
related response. However, non-linear first mode contribute to increments of modal
deformations as the associated incremental displacements will be non-zero as shown in
Equation (4.24).
On the contrary, higher modes in the non-linear phases will maintain their
contribution to base shear and deformation response as their frequencies and incremental
27
displacements are non-zero. On the other hand, it is interesting to note that increments of
modal overturning moments associated with non-linear higher modes shall also be zero at
the base since the product
( )
2
m Γ m ∆d m ωm
y Φ(i) (i) (i) (i)
(4.25)
is calculated as zero for these modes. Thus, it can be stated that non-linear higher modes
initiates non-zero incremental shear forces, which also satisfy zero increments of modal
overturning moments at the base of the wall. Hence, during time intervals of non-linear
phase, incremental shear forces shall continue to develop on the structural wall until the
commencement of unloading in the plastic hinge, which modify the proportionality
assumed between the base shear forces and overturning moments by code procedures.
Figure 4.5 shows the plastic rotation history of the 16 story generic wall designed
with R=6, which has been arranged to demonstrate only the non-linear phases of the wall
with non-zero plastic rotations (linear phases of the wall are shown with zero rotations).
The wall has exhibited the first non-linear response at time t = 4.62 seconds of the record
and sustained this phase for 0.04 seconds with plastic rotations in the positive direction. A
total of 40 non-linear phases has been noted throughout the response history of the generic
wall with 38 of them are associated with plastic rotations in the negative direction which
indicates a response around a residual deformation developed after t = 5.0 seconds.
28
-3
x 10
1
-1
PLASTIC ROTATION (rad)
-2
-3
-4
-5
-6
-7
0 5 10 15 20 25 30
TIME (sec)
Figure 4.5. Plastic rotation history at the base of 16 story generic wall (R=6) obtained from
non-linear time history analysis (Whittier Narrows, 1987-Brea Dam Record)
Modal base shear history of the 16 story generic wall, decomposed with the
technique explained above, is given in Figure 4.6 for the first five modes of the wall. Base
shear response combined with mode superposition method is also provided in last pane of
Figure 4.6. With reference to Figure 4.5, maximum base shear response of the 16 story
generic wall has been observed in the first non-linear phase of the wall which occurred at t
= 4.65 seconds (hereinafter called “the peak response instant”) of the response history as
indicated in Figure 4.6 with a dark circle. This non-linear phase has been observed to
initiate at t = 4.62 seconds and continued through 0.04 seconds to switch to a linear phase
with unloading in the plastic hinge. It is evident from Figure 4.6 that, modal base shear of
the second mode also reaches its first peak value at the peak response instant, where the
contribution of the first mode is comparatively lower. Effect of third and higher modes has
been observed to be negligible at peak response instant as shown in Figure 4.6, indicating
the significant role of the second mode on the base shear response.
Base shear and overturning moment response histories of the 16 story generic wall
computed by the superposition of decomposed modal responses have been compared with
the actual base shear and overturning moment response histories obtained from the Group
29
2 non-linear time history analysis in order to verify the accuracy of the technique as shown
in Figure 4.7. and Figure 4.8, respectively. Both Figure 4.7 and Figure 4.8 demonstrate an
almost perfect match between the responses obtained from mode superposition and non-
linear time history analysis, which indicates that non-linear response can be decomposed
into modal responses by using modal characteristics of the structural walls in an
incremental fashion.
BASE SHEAR (kN)
2000
MODE 1
1000
0
-1000
-2000
0 5 10 15 20 25
BASE SHEAR (kN)
2000
MODE 2
1000
0
-1000
-2000
0 5 10 15 20 25
BASE SHEAR (kN) BASE SHEAR (kN)
2000
MODE 3
1000
0
-1000
-2000
0 5 10 15 20 25
2000
MODE 4
1000
0
-1000
-2000
0 5 10 15 20 25
BASE SHEAR (kN)
2000
MODE 5
1000
0
-1000
-2000
0 5 10 15 20 25
BASE SHEAR (kN)
2000
MODAL SUP.
1000
0
-1000
-2000
0 5 10 15 20 25
TIME (sec)
Figure 4.6. Modal and superposed base shear history of the 16 story generic wall designed
with R=6
30
2000
1500
1000
500
BASE SHEAR (kN)
-500
-1000
-1500
MODAL SUP.
-2000 NL. TIME HIST.
-2500
0 5 10 15 20 25 30
TIME (sec)
Figure 4.7. Agreement between the base shear of the 16 story generic wall obtained from
mode superposition and Group 2 non-linear time history analysis
4
x 10
2.5
2
MODAL SUP.
1.5 NL. TIME HIST.
OVERTUNING MOMENT(kNm)
0.5
-0.5
-1
-1.5
-2
-2.5
0 5 10 15 20 25 30
TIME (sec)
Figure 4.8. Agreement between the overturning moment of the 16 story generic wall
obtained from mode superposition and Group 2 non-linear time history analysis
31
Modal increment of seismic forces associated with the first three modes of the 16
story generic wall obtained immediately before the peak response instant are shown in
Figure 4.9.
12
10
STORY
0
-8 -6 -4 -2 0 2 4 6
MODAL SEISMIC FORCE INCREMENT (kN)
Figure 4.9. Increment of modal seismic forces of the 16 story generic wall designed with
R=6 (Immediately before the peak response instant)
Modal seismic force pattern of the wall at the peak response instant is shown in
Figure 4.10, which emphasizes considerably higher modal seismic forces provided by the
non-linear second mode with respect to other modal contributions.
32
14
12
10
STORY
MODE 1
4
MODE 2
MODE 3
2 MODE 4
MODE 5
0
-300 -200 -100 0 100 200 300
MODAL SEISMIC FORCE (kN)
Figure 4.10. Modal seismic forces of the 16 story generic wall designed with R=6 (at the
peak response instant)
Figure 4.11 demonstrates the significant effect of the non-linear second mode to
story shear forces combined with mode superposition at the peak response instant. Story
shear force quantities are given as absolute values in Figure 4.11 for the ease of
interpretation. It is evident from Figure 4.11 that majority of story shear forces in the non-
linear phases of the wall are provided by the non-linear second mode. Story shear forces
associated with the first mode at the peak response instant shown in Figure 4.11 are those
retained from the prior linear phase, since increments of seismic forces associated with the
first mode are zero in this phase as demonstrated in Figure 4.9. Contribution of the third
and higher modes to story shear force response are virtually negligible based on the
observation that modal increments of seismic forces associated with these modes tend to
reduce the modal shear response in the non-linear phase contrary to the case for the second
non-linear mode.
33
0
0 200 400 600 800 1000 1200 1400 1600 1800 2000
MODAL STORY SHEAR FORCE (kN)
Figure 4.11. Modal and superposed story shear forces of the 16 story generic wall designed
with R=6 (at the peak response instant)
Effect of the non-linear second mode has also been observed on the moment
response of the 16 story generic wall as shown in Figure 4.12 in terms of modal increment
of moments obtained immediately before the peak response instant. Increment of modal
moments associated with the first non-linear mode has been obtained as zero due to zero
increment of seismic forces associated with this mode. However, second non-linear mode,
with zero modal moment increment at the base, has been observed to induce moments at
upper story levels. This effect may contribute to the spread of plasticity at the upper levels
of the wall in the non-linear phases as observed in Group 1 and Group 3 non-linear time
history analysis (See Chapter 5).
34
10
STORY
0
-350 -300 -250 -200 -150 -100 -50 0 50 100
MODAL MOMENT INCREMENT (kNm)
Figure 4.12. Increment of modal moments of the 16 story generic wall designed with R=6
(Immediately before the peak response instant)
Dynamic base shear amplification factors can also be presented in terms of modal
contributions. Figure 4.13 shows first ten seconds of modal dynamic amplification factor
history of the 16 story generic wall, obtained by dividing each modal base shear to the
design base shear of the wall. Dynamic base shear amplification factor combined with
mode superposition method is also provided in last pane of Figure 4.13 and dynamic
amplification factor at the peak response instant has been marked with a dark dot for the
ease of understanding.
35
3
2 MODE 1
1
0
-1
-2
-3
0 1 2 3 4 5 6 7 8 9 10
DYN. AMP. FACT.
3
2 MODE 2
1
0
-1
-2
-3
0 1 2 3 4 5 6 7 8 9 10
DYN. AMP. FACT.
3
2 MODE 3
1
0
-1
-2
-3
0 1 2 3 4 5 6 7 8 9 10
DYN. AMP. FACT.
3
2 MODE 4
1
0
-1
-2
-3
0 1 2 3 4 5 6 7 8 9 10
DYN. AMP. FACT.N
3
2 MODE 5
1
0
-1
-2
-3
0 1 2 3 4 5 6 7 8 9 10
DYN. AMP. FACT.
3
2 MODAL SUP.
1
0
-1
-2
-3
0 1 2 3 4 5 6 7 8 9 10
TIME (sec)
Figure 4.13. First ten seconds of the modal and superposed dynamic base shear
amplification factor history of the 16 story generic wall designed with R=6
Figure 4.13 indicates that, non-linear second mode is chiefly responsible for the
dynamic base shear amplifications, where effect of higher modes above the second mode
has been found to be negligible at the peak response instant. Dynamic base shear
amplification factor associated with the first mode has been computed greater than unity
between 6.30 seconds and 6.40 seconds of the response history (Figure 4.13), due to slight
numerical errors in the decomposition process.
Significant role of the second mode on dynamic shear amplifications has been
demonstrated in this study by utilizing a modal decomposition technique. It has been
shown that, response of yielding structural walls can be decomposed into modal response
contributions by using the modal decomposition technique explained above.
37
In this chapter, base shear amplification factors and their relation with the first mode
period of walls and the strength reduction factors have been investigated through a
parametric study based on non-linear time history analyses. Base shear amplification
factors are proposed based on regression analyses of the parametric time history analysis
results, which are indented to be used in Turkish Seismic Design Code (2007) for shear
design of walls. Furthermore, a story shear force envelope along the wall height is
recommended and interrelated with the amplified base shear.
The generic walls with geometric characteristics given in Chapter 3 have been
analyzed with non-linear time history analyses by using step-by-step integration procedure
in the time domain. Each wall has been consecutively analyzed for the strength levels
corresponding to strength reduction factors of R=2, 4, 6, 8 and 10 with 20 strong ground
motion acceleration records scaled to fit the Turkish Seismic Design Code (2007)
Acceleration Response Spectrum as explained in Appendix C.
Newmark’s average acceleration method (γ=1/2 and β=1/4) has been used as an
unconditionally stable algorithm in step-by-step integration of non-linear time history
analyses. In order to avoid the accuracy problems, time steps has been chosen as 0.005 sec,
which is below 1/10 of the third mode period of all generic walls.
Damping matrix of the system has been constructed with Rayleigh method that is,
the damping matrix has been assumed to be proportional to the stiffness and mass matrices
of the walls. Damping ratios for the first and the third modes have been chosen to be 1 per
38
cent, so as not to over emphasize the damping effects in the second mode as illustrated in
Figure 5.1
ξ 1= ξ3= 0.01
ξ2
ω1 ω2 ω3 ω
Figure 5.1. Illustration of damping ratio assumption in time history analyses
Non-linear time history analyses has been executed under four groups (hereinafter
called Group 1, Group 2, Group 3 and Group 4) in order to investigate the sensitivity of the
base shear amplification factors to the distribution of plastic hinges along the wall height
and the hysteresis characteristics of the plastic hinges. Modeling details and base shear
amplification results for each group are presented in the subsequent sections. Due to vast
number of analyses performed in this study, only the moment and story shear force
diagrams along the wall height, that have been obtained from the generic 16 story wall
with strength reduction factor of R=6 has been presented.
Group 1 and Group 2 non-linear time history analyses have been performed by using
SAP 2000 v11 (CSI, 2006) structural analysis software, whereas Group 3 and Group 4
non-linear time history analyses have been performed by using RUAUMOKO (Carr, 2000)
39
structural analysis software, which are both licensed to the Department of Earthquake
Engineering of Boğaziçi University.
Group1 analyses involve the generic wall analytical models of multiple plastic hinge
distribution along the wall height with elastoplastic moment-curvature hysteresis relations
as shown in Figure 5.2. Other structural modeling assumptions described in Chapter 3 are
also valid for the Group 1 models.
STORY j MASS
PLASTIC HINGE
STORY j-1 MASS MOMENT
PLASTIC HINGE
MASS
PLASTIC HINGE
MASS
CURVATURE
PLASTIC HINGE
STORY 3 MASS
PLASTIC HINGE
STORY 2 MASS
ELASTOPLASTIC HYSTERESIS
PLASTIC HINGE RELATION FOR PLASTIC HINGES
STORY 1 MASS
PLASTIC HINGE
It has been assumed that yield moment capacity is uniform along the wall, being
equal to base design moment corresponding to the strength reduction factors as described
in Chapter 3. Actually, reinforced concrete wall sections do not exhibit perfectly
elastoplastic hysteresis characteristics during flexure as reversals due to seismic attack
become significant, and some degree of stiffness degradation should be expected. Yet,
elastoplastic plastic hinges are customarily used in most non-linear time history analyses
and it is believed that reliable results can be obtained for the benchmarking of studies with
more complex hysteresis relationships.
Non-linear time history analyses have been performed for generic wall models of 8,
12, 16, 20 and 30 stories with strength levels corresponding to R = 2, 4, 6, 8 and 10. For
40
each generic wall model with the aforementioned strength level, maximum story moments
and shear forces have been obtained for each of 20 ground motion record as shown in
Figure 5.3 and Figure 5.4 for 16 story wall with strength level of R = 6.
16
14
12
10
STORY
0
0 0.5 1 1.5 2 2.5
MOMENT (kNm) 4
x 10
Figure 5.3. Story moment diagrams of 16 story wall with strength level R=6 (Group 1 time
history analyses)
41
16
14
12
10
STORY
0
0 500 1000 1500 2000 2500 3000
STORY SHEAR FORCE (kN)
Figure 5.4. Story shear force diagrams of 16 story wall with strength level R=6 (Group 1
time history analyses)
Subsequently story moment and shear forces have been averaged to obtain mean
story moment and story shear force diagrams, as shown in Figure 5.5 and Figure 5.6 for 16
story wall with R = 6 strength level with comparison to the design moment and shear force
diagrams.
42
16
MEAN (T.HIST)
14 DESIGN
12
10
STORY
0
0 0.5 1 1.5 2 2.5
MOMENT (kNm) 4
x 10
Figure 5.5. Mean story moment diagram of 16 story wall with strength level R=6 (Group 1
time history analyses)
16
MEAN (T.HIST)
14 DESIGN
12
10
STORY
0
0 500 1000 1500 2000 2500
STORY SHEAR FORCE (kN)
Figure 5.6. Mean story shear force diagram of 16 story wall with strength level R=6
(Group 1 time history analyses)
Consequently mean story shear forces obtained from time history analyses have been
divided to design story shear forces in order to obtain story shear force amplification
factors as shown in Figure 5.7 for the 16 story wall with R = 6 strength level.
43
16
14
12
10
STORY
0
1.5 2 2.5 3 3.5 4
STORY SHEAR FORCE AMPLIFICATION
Figure 5.7. Mean story shear force amplification factors computed for the 16 story wall
with strength level R=6 (Group 1 time history analyses)
As shown in Figure 5.4, there is a scatter of story shear forces due to the different
characteristics of the ground motion in spite of scaling has been applied to fit a target
response spectrum. This scatter has been observed at all generic walls and corresponding
strength levels and is obviously reflected to the base shear amplification factors; however it
is believed that averaging time history results from 20 ground motion records is justifiable
due to the sufficient number of samples in the group.
Mean base shear amplification factors derived from Group 1 non-linear time history
analyses have been collected in a single diagram (βb – T – R diagram) related with the first
mode periods of the wall corresponding the cracked section stiffness and strength reduction
factors as shown in Figure 5.8
44
GROUP 1 ANALYSES
4
R=2
R=4
2.5
1.5
1
0.5 1 1.5 2 2.5
CRACKED SECTION ELASTIC FIRST MODE PERIOD (sec)
Figure 5.8. Mean base shear amplification factors obtained from Group 1 time history
analyses
45
Group 2 analyses include the generic wall analytical models with a single plastic
hinge at the base of the wall with elastoplastic moment-curvature hysteresis relation as
shown in Figure 5.9. Other structural modeling assumptions described in Chapter 3 are
also valid for the Group 2 models.
STORY j MASS
MASS
MASS
CURVATURE
STORY 3 MASS
STORY 2 MASS
ELASTOPLASTIC HYSTERESIS
RELATION FOR PLASTIC HINGE
STORY 1 MASS
PLASTIC HINGE
It is expected in real wall structures that a single plastic hinge formation at the base
of the wall is not possible and the plasticity spreads upwards from the base with increasing
seismic demand. However, it is intended that Group 2 analyses shall provide insight on the
effect of spread of plasticity on structural response with respect to Group 1 analyses.
Non-linear time history analyses have been performed for Group 2 generic wall
models similar to the Group 1 analyses. Similarly, for each generic wall model with the
corresponding strength level, maximum story moments and shear forces have been
obtained for each of 20 ground motion record. Subsequently mean story moments, story
shear forces and base shear amplification factors have been computed as shown in Figure
5.9 to Figure 5.14 for the 16 story wall with R = 6 strength level.
46
16
14
12
10
STORY
0
0 0.5 1 1.5 2 2.5 3
MOMENT (kNm) 4
x 10
Figure 5.10. Story moment diagrams of 16 story wall with strength level R=6 (Group 2
time history analyses)
16
14
12
10
STORY
0
0 500 1000 1500 2000 2500
STORY SHEAR FORCE (kN)
Figure 5.11. Story shear force diagrams of 16 story wall with strength level R=6 (Group 2
time history analyses)
47
16
MEAN (T.HIST)
14 DESIGN
12
10
STORY
0
0 0.5 1 1.5 2 2.5
MOMENT (kNm) 4
x 10
Figure 5.12. Mean story moment diagram of 16 story wall with strength level R=6 (Group
2 time history analyses)
16
MEAN (T.HIST)
14 DESIGN
12
10
STORY
0
0 500 1000 1500 2000 2500
STORY SHEAR FORCE (kN)
Figure 5.13. Mean story moment diagram of 16 story wall with strength level R=6 (Group
2 time history analyses)
48
16
14
12
10
STORY
0
1.5 2 2.5 3 3.5 4
STORY SHEAR FORCE AMPLIFICATION
Figure 5.14. Mean story shear force amplification factors computed for the 16 story wall
with strength level R=6 (Group 2 time history analyses)
4.5 R=4
R=6
4 R=8
R=10
3.5
2.5
1.5
1
0.5 1 1.5 2 2.5
CRACKED SECTION ELASTIC FIRST MODE PERIOD (sec)
Figure 5.15. Mean base shear amplification factors obtained from Group 3 time history
analyses
49
Group 3 analyses involve the generic wall analytical models with multiple plastic
hinge distribution along the wall height with Takeda (1970) degrading stiffness moment-
curvature hysteretic relationships as shown in Figure 5.16. Other structural modeling
assumptions depicted in Chapter 3 are also valid for the Group 3 models.
STORY j MASS
PLASTIC HINGE
STORY j-1 MASS MOMENT
PLASTIC HINGE
MASS
PLASTIC HINGE
MASS
CURVATURE
PLASTIC HINGE
STORY 3 MASS
PLASTIC HINGE
STORY 2 MASS
TAKEDA HYSTERESIS RELATION
PLASTIC HINGE FOR PLASTIC HINGES
STORY 1 MASS
PLASTIC HINGE
It has been intended that Group 3 analyses shall reflect the effect of degrading
stiffness moment-curvature hysteretic relationship with respect to the Group 1 analyses
that were performed with elastoplastic hysteretic relationship.
Non-linear time history analyses have been performed for the generic wall models in
the Group 3 analyses similar to Group 1 and Group 2 analyses. For each generic wall
model with the corresponding strength level, maximum story moments and shear forces
have been obtained for each of 20 ground motion records. Subsequently mean story
moments, story shear forces and base shear amplification factors have been computed as
shown in Figure 5.17 to Figure 5.21 for the 16 story wall with R = 6 strength level.
50
16
14
12
10
STORY
0
0 0.5 1 1.5 2 2.5
MOMENT (kNm) 4
x 10
Figure 5.17. Story moment diagrams of 16 story wall with strength level R=6 (Group 3
time history analyses)
16
14
12
10
STORY
0
0 500 1000 1500 2000 2500 3000
STORY SHEAR FORCE (kN)
Figure 5.18. Story shear force diagrams of 16 story wall with strength level R=6 (Group 3
time history analyses)
51
16
MEAN (T.HIST)
14 DESIGN
12
10
STORY
0
0 0.5 1 1.5 2 2.5
MOMENT (kNm) 4
x 10
Figure 5.19. Mean story moment diagram of 16 story wall with strength level R=6 (Group
3 time history analyses)
16
MEAN (T.HIST)
14 DESIGN
12
10
STORY
0
0 200 400 600 800 1000 1200 1400 1600 1800 2000
STORY SHEAR FORCE (kN)
Figure 5.20. Mean story shear force diagram of 16 story wall with strength level R=6
(Group 3 time history analyses)
52
16
14
12
10
STORY
0
1.5 2 2.5 3 3.5 4 4.5
STORY SHEAR FORCE AMPLIFICATION
Figure 5.21. Mean story shear force amplification factors computed for the 16 story wall
with strength level R=6 (Group 3 time history analyses)
R=6
3.5
R=8
R=10
2.5
1.5
1
0.5 1 1.5 2 2.5
CRACKED SECTION ELASTIC FIRST MODE PERIOD (sec)
Figure 5.22. Mean base shear amplification factors obtained from Group 3 time history
analyses
53
Group 4 analyses constitute the generic wall analytical models with a single plastic
hinge at the base of the wall height Takeda (1970) degrading stiffness moment-curvature
hysteresis relations as shown in Figure 5.23. Other structural modeling assumptions
described in Chapter 3 are also valid for the Group 4 models.
STORY j MASS
MASS
MASS
CURVATURE
STORY 3 MASS
PLASTIC HINGE
In spite of an unrealistic behavior expected with a single plastic hinge at the base of
the wall, it has been anticipated that Group 4 analyses shall provide insight on the effect of
spread of plasticity on structural response with respect to Group 3 analyses.
Similar to the other analysis groups, non-linear time history analyses have been
performed for the generic wall models in the Group 4 analyses. For each generic wall
model with the corresponding strength level, maximum story moments and shear forces
have been obtained for each of 20 ground motion record. Subsequently mean story
moments, story shear forces and base shear amplification factors have been computed as
shown in Figure 5.24 to Figure 5.28 for the 16 story wall with R = 6 strength level.
54
16
14
12
10
STORY
0
0 0.5 1 1.5 2 2.5 3
MOMENT (kNm) 4
x 10
Figure 5.24. Story moment diagrams of 16 story wall with strength level R=6 (Group 4
time history analyses)
16
14
12
10
STORY
0
0 500 1000 1500 2000 2500
STORY SHEAR FORCE (kN)
Figure 5.25. Story shear force diagrams of 16 story wall with strength level R=6 (Group 4
time history analyses)
55
16
MEAN (T.HIST)
14 DESIGN
12
10
STORY
0
0 0.5 1 1.5 2 2.5
MOMENT (kNm) 4
x 10
Figure 5.26. Mean story moment diagram of 16 story wall with strength level R=6 (Group
4 time history analyses)
16
MEAN (T.HIST)
14 DESIGN
12
10
STORY
0
0 200 400 600 800 1000 1200 1400 1600 1800 2000
STORY SHEAR FORCE (kN)
Figure 5.27. Mean story shear force diagram of 16 story wall with strength level R=6
(Group 4 time history analyses)
56
16
14
12
10
STORY
0
1.5 2 2.5 3 3.5 4
STORY SHEAR FORCE AMPLIFICATION
Figure 5.28. Mean story shear force amplification factors computed for the 16 story wall
with strength level R=6 (Group 4 time history analyses)
R=2
MEAN BASE SHEAR AMPLIFICATION FACTOR
4.5 R=4
R=6
4 R=8
R=10
3.5
2.5
1.5
1
0.5 1 1.5 2 2.5
CRACKED SECTION ELASTIC FIRST MODE PERIOD (sec)
Figure 5.29. Mean base shear amplification factors obtained from Group 4 time history
analyses
57
Based on the information gathered from the Group 1, Group 2, Group 3 and Group 4
non-linear time history analyses results, the following points have been observed;
Based on the remarks stated in Chapter 5, it has been concluded that multiple hinge
distribution along the wall better present the expected wall behavior during seismic attack.
Moreover, it is believed that Takeda (1970) degrading stiffness hysteresis model is more
suitable to represent realistic reinforced concrete flexural behavior. Thus, regression
analysis have been performed for the base shear amplification factors derived from Group
3 non-linear time history analyses in order to propose a relationship between the base shear
dynamic amplification factor, first mode period of the wall corresponding to cracked
section stiffness and the strength reduction factor. It is suggested that the relationship
proposed herein can be used for the design and evaluation of structural wall systems in
accordance with Turkish Seismic Design Code (2007) provisions.
which formed the proposed expression given in Equation (6.1) for the dynamic base shear
amplification factors as presented in Equation (6.2);
The coefficient of multiple determination of the regression (R2) has been computed
as 0.9743 over 1.0000, which suggested a successful fit of the proposed relationship to the
samples. The agreement between the obtained mean dynamic base shear amplification
59
factors and the proposed dynamic base shear amplification relationship is shown in Figure
6.1.
4.5
R=2
R=4
4
BASE SHEAR AMPLIFICATION FACTOR
R=6
R=8
3.5 R=10
2.5
1.5
1
0.5 1 1.5 2 2.5
CRACKED SECTION ELASTIC FIRST MODE PERIOD (sec)
Figure 6.1. Agreement between the proposed relationship and mean base shear
amplification factors
Figure 6.1 suggests that the proposed relationship is in good agreement with the
mean amplification factors for strength reduction factors between 2 and 6, and
overestimates the amplification factors for strength reduction factors 8 and 10 in the high
period range that is believed to be acceptable in terms of a conservative approach.
60
Proposed dynamic base shear amplification relationship given in Equation (6.2) have
been compared with the relationships proposed by other researchers mentioned in Chapter
2. It should be noted that such a comparison may not be appropriate since the base shear
amplification factors are highly affected by the ground motion intensity considered,
number of ground motion records considered and modeling assumptions. However, it is
found useful to compare the trends between the proposed relationship and those proposed
by other researchers.
Figure 6.2 shows the base shear amplification factors computed for the generic walls
considered in this study by using the relationship proposed by Ghosh and Markevicius
(1990). In order to illustrate the deviation between the two relationships base shear
amplification factors computed with Ghosh and Markevicius relationship is divided by the
amplification factors obtained from the proposed relationship as shown in Figure 6.3
Figure 6.3 reveals that proposed relationship in this study provides greater base shear
amplification ratios particularly in the high period range with respect to the relationship
proposed by Ghosh and Markevicius, which may be attributed to lack of a period term in
the Ghosh and Markevicius relationship.
61
3.5
R=2
R=4
2.5
1.5
1
0.5 1 1.5 2 2.5
CRACKED SECTION ELASTIC FIRST MODE PERIOD (sec)
Figure 6.2. Base shear amplification factors computed by using Ghosh and Markevicius
(1990) relationship
0.95 R=2
R=4
0.9 R=6
R=8
RATIO GHOSH / SUGGESTED
0.85 R=10
0.8
0.75
0.7
0.65
0.6
0.55
0.5
0.5 1 1.5 2 2.5
CRACKED SECTION ELASTIC FIRST MODE PERIOD (sec)
Figure 6.3. Ratio of the base shear amplification factors between the relationship suggested
by Ghosh and Markevicius (1990) and this study
62
Figure 6.4 shows the base shear amplification factors by using the relationship
proposed by Eibl and Keintzel (1988). Similarly base shear amplification factors computed
with Eibl and Keintzel relationship is divided by the amplification factors obtained from
the proposed relationship as shown in Figure 6.5 for comparison.
5
R=2
4.5 R=4
BASE SHEAR AMPLIFICATION FACTOR
R=6
R=8
4
R=10
3.5
2.5
1.5
1
0.5 1 1.5 2 2.5
CRACKED SECTION ELASTIC FIRST MODE PERIOD (sec)
Figure 6.4. Base shear amplification factors computed by using Eibl and Keintzel (1988)
relationship
63
2
R=2
R=4
R=6
R=8
R=10
RATIO EIBL / SUGGESTED
1.5
0.5
0.5 1 1.5 2 2.5
CRACKED SECTION ELASTIC FIRST MODE PERIOD (sec)
Figure 6.5. Ratio of the base shear amplification factors between the relationship suggested
by Eibl and Keintzel (1988) and this study
Figure 6.5 indicates the lower amplification factors computed by Eibl and Keintzel
relationship where parallel-like trends of the ratios with respect to strength reduction
factors also suggests the idea that amplification factors have been affected by the different
intensities of ground motion taken into account between the two studies.
Base shear amplification factors by using the relationship proposed by Rutenberg and
Nsieri (2006) is illustrated in Figure 6.6, whereas the ratio between the base shear
amplification factors computed with Rutenberg and Nsieri relationship and the proposed
relationship is shown in Figure 6.7 for comparison.
64
6
R=2
5.5 R=4
3.5
2.5
1.5
1
0.5 1 1.5 2 2.5
CRACKED SECTION ELASTIC FIRST MODE PERIOD (sec)
Figure 6.6. Base shear amplification factors computed by using Rutenberg and Nsieri
(2006) relationship
2
R=2
R=4
R=6
RATIO RUTENBERG / SUGGESTED
R=8
R=10
1.5
0.5
0.5 1 1.5 2 2.5
CRACKED SECTION ELASTIC FIRST MODE PERIOD (sec)
Figure 6.7. Ratio of the base shear amplification factors between the relationship suggested
by Rutenberg and Nsieri (2006) and this study
65
Figure 6.7 depicts slightly higher amplification factors computed by Rutenberg and
Nsieri in which the ratio trends are almost in a linear fashion with increasing strength
reduction factor and period. This result is believed to be related to the higher intensity SAC
ground motion suite (Somerville et al. 1997) used by Rutenberg and Nsieri in their studies.
It should be emphasized that despite the lack of consensus between the researchers
on the magnitude of base shear amplification factors, the above-given results address the
fact that base shear amplification factors increase with increasing strength reduction factor
and increasing first mode period of structural walls. It is believed that suggested base shear
amplification relationship in this study is suitable for the design of structural walls against
shear failures according to Turkish Seismic Design Code (2007).
6.2. Proposed Story Shear Force Profile Along the Wall Height
Mean shear force distribution patterns derived from the non-linear time history
analyses reveal that amplified shear force magnitudes drop gradually to half value at about
40 per cent of the total wall height and follow a practically uniform path up to the tip of the
wall excluding the top one or two stories as shown in Figure 6.8. It is interesting to note
that, normalized story shear force profiles shown in Figure 6.8 encompass all generic walls
designed with the strength reduction factors considered in this study (R=2, 4, 6, 8 and 10),
which indicates that amplified story shear force pattern is neither a function of first mode
period nor the strength reduction factor.
It is believed that a constant shear force zone just above the base of the wall, equal to
the critical wall height as per Turkish Seismic Design Code (2007) will be rational
considering the tension shift effects (Paulay and Priestley, 1992) in the plastic hinge zone.
Above this zone, it is suggested that the story shear forces shall linearly descend to half
magnitude at the 40 per cent of the total wall height and remain constant along the rest of
the wall as shown in Figure 6.9. In cases where 40 per cent of the wall height do not
coincide with a story level, upper end of the linearly descending zone point may be
rounded up to the nearest upper story level.
66
1
OBSERVED TREND
0.9 NL. TIME HIST.
0.8
0.7
NORMALIZED HEIGHT
0.6
0.5
0.4
0.3
0.2
0.1
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
NORMALIZED SHEAR FORCE
Figure 6.8. Observed trend of the shear force profile obtained from non-linear time history
analyses of generic walls
60% Hw
40% Hw
max(l w , Hw /6)
0.5 1.0
Shear force profile proposed in this study is similar to that given in Eurocode 8 EN
1998-1 (CEN, 2004) in most aspects, as shown in Figure 6.10. Major notable differences
67
between the two profiles are the level where the shear forces drop approximately to 50 per
cent of the base shear is defined as 1/3 (33 per cent) of the wall height and no constant
shear in the tension shift zone at the base of the wall is taken into account.
67% H w
33% Hw
0.5 1.0
Figure 6.10. Eurocode 8.1 (CEN, 2004) story shear force profile
Rutenberg and Nsieri (2006) also recommended a similar story shear force profile
based on the non-linear time history results as illustrated in Figure 6.11. The plastic hinge
zone above the base has been recommended as 10 per cent of the total height by Rutenberg
and Nsieri (2006) and the level where the shear forces drop 50 per cent of the base shear
has been related with first mode period of the wall with the expression given in Equation
(6.4)
It is believed that shear force profile proposed by Rutenberg and Nsieri over punishes
the structures in the medium period range, since ξ corresponds to values in the order of 70
per cent and 85 per cent of the total wall height for a structural wall with 1.0 sec and 0.50
sec period, respectively.
68
ξH w
0.10H w
0.5 1.0
Figure 6.11. Shear force profile proposed by Rutenberg and Nsieri (2006)
few non-linear time history analyses should be performed to support the proposed
relationship if the structure is complex in geometry and/or configuration or the importance
level of the structure is high. It should also be noted that the proposed relationship has been
derived from non-linear time history analyses of generic structural walls covering a limited
period range. Thus, it is recommended that the proposed relationship should not be used
for structural wall systems above 30 stories.
In majority of the cases of design practice, even flexural strength of the wall
computed with minimum reinforcement ratios can well exceed the flexural strength
required by analysis. This provided flexural overstrength effectively reduce the strength
reduction factor used in design since the plastic hinge at the base of the wall will tend to
yield at a higher moment demand than required by analysis. Thus, Equation (6.2) can be
modified to take into account the aforementioned flexural overstrength as shown in
Equation (6.5)
R
β b = 1.0 + (0.281× T1-cr + 0.394) × ( − 1.5)0.553 (6.5)
Ψ o
M rb
Ψo = (6.6)
M eb
In Equation (6.6), M eb denotes the moment demand at the base of the wall obtained
from code procedures and M rb denotes the provided moment strength of the section at the
base of the wall computed with code requirements.
2. Multiply the design base shear ( Veb ) with base shear amplification factor ( β b )
computed by using Equation (6.5) to obtain amplified design base shear ( Vab ) as
shown in Equation (6.7).
3. Multiply the amplified design base shear ( Vab ) with the unitless story shear profile
proposed in Section 6.2 to obtain the amplified story shear forces along the wall.
Even though the proposed dynamic base shear amplification relationship is intended
for use in design, it can be employed in the seismic assessment of structural wall systems
by using linear and non-linear analysis procedures as well. Strength reduction factor
obtained by using single mode pushover analysis procedures can be implemented in
Equation (6.5) in order to compute the amplified base shear forces. Implementation of a
dynamic base shear amplification factor is not necessary in multi-mode pushover analysis
(Aydınoğlu, 2003) since the analysis procedure takes the effects of higher modes into
account. For the seismic assessment of structural wall systems by using linear analysis
procedures, third term in Equation (6.5) may be modified as shown in Equation (6.8).
M Eb
Ry = (6.9)
M nb
In Equation (6.9), M Eb denotes the elastic moment demand at the base of the wall
obtained from code procedures and M nb denotes the nominal moment strength of the
section at base of the wall associated with characteristic material strengths.
71
Some degree of sliding shear in the plastic hinge zone of walls is inevitable
particularly at high ductility demands. Recent tests of walls performed on E-Defense
shaking table in Japan (NIED, 2006) and UC San Diego shaking table in the US (UCSD,
2005), indicate the dominance of sliding shear in the plastic hinge zone on the wall
behavior (Krawinkler, 2006). The severity of sliding shear increase with increasing
ductility demand and this phenomenon causes strength deterioration as the number of
cycles experienced by the plastic hinge. Thus, most of the sliding shear resistance will be
provided by the dowel action of the longitudinal reinforcement as a consequence of open
cracks and deteriorated concrete strength in the plastic hinge. Hence, it may be rational to
define an upper bound for base shear amplifications where the behavior will be primarily
governed by sliding at the base of the wall after attaining sliding shear strength.
Based on the recommendations of Paulay and Priestley (1992) for the sliding shear
resistance of longitudinal reinforcement in squat walls, sliding shear resistance at the base
of a of cantilever wall can be predicted by Equation (6.10);
in which, Vfb is the sliding resistance provided by the longitudinal reinforcement, Asw is the
total longitudinal reinforcement area in the web of the section at the base of the wall and fyk
is the characteristic yield strength of reinforcement. Thus, upper bound of the base shear
amplification factor can be expressed as given in Equation (6.11)
0.25 Asw f yk
β max
b
= (6.11)
Veb
It is believed that the upper bound limitation of base shear amplifications as given in
Equation (6.11) will provide rational results for design and prevent unnecessary overdesign
in practice.
72
As a side product of this study, a unitless moment profile along the wall height is
suggested for the Turkish Seismic Design Code (2007), based on the regression of
normalized moments obtained from non-linear time history analyses of the generic walls. It
has been observed in this study that moment profile is strongly dependent on the strength
reduction factor utilized, due to the increasing spread of plasticity above the base of the
wall with increasing strength reduction factor. However, such a dependency was not
observed for the first mode period. A tri-linear moment profile is suggested as shown in
Figure 6.12.
Hw
ηH w
δHw
λ 1.0
Agreement of the proposed moment profile with those obtained from non-linear time
history analyses are shown in Figure 6.13 to Figure 6.17 with respect to the strength
reduction factor utilized.
1
PROPOSED PROFILE
0.9 NL. TIME HIST.
0.8
0.7
NORMALIZED HEIGHT
0.6
0.5
0.4
0.3
0.2
0.1
0
0 0.2 0.4 0.6 0.8 1
NORMALIZED MOMENT
Figure 6.13. Agreement between proposed moment profile and non-linear time history
analysis of the generic walls designed with R=2
74
1
PROPOSED PROFILE
0.9 NL. TIME HIST.
0.8
0.7
NORMALIZED HEIGHT
0.6
0.5
0.4
0.3
0.2
0.1
0
0 0.2 0.4 0.6 0.8 1
NORMALIZED MOMENT
Figure 6.14. Agreement between proposed moment profile and non-linear time history
analysis of the generic walls designed with R=4
1
PROPOSED PROFILE
0.9 NL. TIME HIST.
0.8
0.7
NORMALIZED HEIGHT
0.6
0.5
0.4
0.3
0.2
0.1
0
0 0.2 0.4 0.6 0.8 1
NORMALIZED MOMENT
Figure 6.15. Agreement between proposed moment profile and non-linear time history
analysis of the generic walls designed with R=6
75
1
PROPOSED PROFILE
0.9 NL. TIME HIST.
0.8
0.7
NORMALIZED HEIGHT
0.6
0.5
0.4
0.3
0.2
0.1
0
0 0.2 0.4 0.6 0.8 1
NORMALIZED MOMENT
Figure 6.16. Agreement between proposed moment profile and non-linear time history
analysis of the generic walls designed with R=8
1
PROPOSED PROFILE
0.9 NL. TIME HIST.
0.8
0.7
NORMALIZED HEIGHT
0.6
0.5
0.4
0.3
0.2
0.1
0
0 0.2 0.4 0.6 0.8 1
NORMALIZED MOMENT
Figure 6.17. Agreement between proposed moment profile and non-linear time history
analysis of the generic walls designed with R=10
76
It should be noted that tension shift effect (Paulay and Priestley 1992) on the
proposed moment profile has not been taken into account. Equation (6.13) and Equation
(6.14) should be modified as shown below to incorporate the tension shift effect in
accordance with the Turkish Seismic Design Code (2007):
lw
δ = 0.003 × ( R − 2) 2 + 0.025 × ( R − 2) + max( , 0.167) (6.15)
Hw
l
η = 0.125 × ( R − 2) − 0.005 × ( R − 2) 2 + max( w , 0.167) (6.16)
Hw
77
7. CONCLUSIONS
Amplified shear forces, particularly at the base of the wall, are probable to cause
unexpected premature brittle shear failure modes of diagonal tension, diagonal
compression and sliding before the intended plastic hinge development if the magnified
shear forces are not taken into account. Practically dynamic shear amplification is taken
into account by means of base shear amplification factors applied to the design base shear
calculated with the code procedures. These amplification factors are presented principally
as a function of first mode period and/or strength reduction factors. 2007 version of the
Turkish Seismic Design Code considers the base shear amplification phenomena with a
constant base shear amplification factor of 1.5 regardless of the first mode period and
ductility level of the wall. It is believed that this must be revised with base shear
amplification factors as functions of strength reduction factor and first mode period of the
wall. It is also believed that a shear force profile along the height of the wall is also
necessary in the Turkish Seismic Design Code (2007) not only for preventing shear
failures at the base but also along the height of the structural wall.
Extensive parametric non-linear time history analyses have been performed for code
designed generic structural walls, in order to suggest a base shear amplification
relationship and a shear force profile for the Turkish Seismic Design Code (2007).
Distribution and hysteretic characteristics of the plastic hinges have been considered in the
analyses models of the parametric study, in order to investigate the effect of modeling
assumptions to dynamic shear amplifications. It has been observed that dynamic shear
amplifications are not dependent on the hysteretic characteristics of the plastic hinge.
However, models with single plastic hinge at the base of the wall resulted in comparatively
78
higher story shear force amplification factors and incorrect moment patterns with respect to
distributed plastic hinge models. Results of the non-linear time history analysis revealed
that dynamic base shear amplification factors increase with increasing first mode period as
well as with strength reduction factor. A base shear amplification relationship as functions
of first mode period and strength reduction factor has been proposed based on a regression
analysis of the non-linear time history analysis results. It was suggested that dynamic base
shear amplifications at the base of the wall should be limited with the amplifications
corresponding to the sliding shear strength during load reversals for a better
characterization of the wall behavior under seismic attack. Moreover, a practical story
shear force profile has been suggested, based on the non-linear time history analysis
results, to be used in tandem with the Turkish Seismic Design Code (2007). It has been
observed that story shear force profile is neither a function of first mode period nor the
strength reduction factor. A moment profile has also been proposed for use in the Turkish
Seismic Design Code (2007) as a side product of the parametric non-linear time history
analyses performed in this study. Moment profiles obtained from non-linear time history
analyses indicated a strong dependency on the strength reduction factor utilized but not the
first mode period of the wall. A novel modal decomposition technique is presented in this
study for demonstrating the individual effects of the higher modes on the dynamic shear
amplification phenomenon. The technique is based on inverse extraction of modal
responses from non-linear time history analyses results by using the mode superposition
method. It was shown that non-linear second mode is chiefly responsible for the dynamic
base shear amplifications and the effect of higher modes above the second mode has been
found to be negligible. It is believed that the modal decomposition technique presented in
this study demonstrates the individual contribution of higher mode effects to dynamic
amplification where related studies on this subject are very limited in the literature.
real dynamic behavior of yielding structural walls as more full-scale tests and supporting
numerical research simulating the complex reinforced concrete behavior are carried out.
80
8
R=2
7 R=4
R=6
R=8
6
R=10
5
STORY
0
0 0.5 1 1.5 2 2.5
DESIGN MOMENT (kNm) 4
x 10
8
R=2
7 R=4
R=6
R=8
6
R=10
5
STORY
0
0 200 400 600 800 1000 1200 1400
DESIGN STORY SHEAR FORCE (kN)
12
R=2
R=4
10 R=6
R=8
R=10
8
STORY
0
0 0.5 1 1.5 2 2.5 3 3.5 4
DESIGN MOMENT (kNm) 4
x 10
12
R=2
R=4
10 R=6
R=8
R=10
8
STORY
0
0 200 400 600 800 1000 1200 1400 1600 1800
DESIGN STORY SHEAR FORCE (kN)
16
R=2
14 R=4
R=6
R=8
12
R=10
10
STORY
0
0 1 2 3 4 5 6 7
DESIGN MOMENT (kNm) 4
x 10
16
R=2
14 R=4
R=6
R=8
12
R=10
10
STORY
0
0 500 1000 1500 2000 2500
DESIGN STORY SHEAR FORCE (kN)
20
R=2
18 R=4
R=6
16
R=8
R=10
14
12
STORY
10
0
0 1 2 3 4 5 6 7 8 9 10
DESIGN MOMENT (kNm) 4
x 10
20
R=2
18 R=4
R=6
16
R=8
R=10
14
12
STORY
10
0
0 500 1000 1500 2000 2500 3000
DESIGN STORY SHEAR FORCE (kN)
30
R=2
R=4
25 R=6
R=8
R=10
20
STORY
15
10
0
0 0.5 1 1.5 2 2.5
DESIGN MOMENT (kNm) 5
x 10
30
R=2
R=4
25 R=6
R=8
R=10
20
STORY
15
10
0
0 500 1000 1500 2000 2500 3000 3500 4000 4500
DESIGN STORY SHEAR FORCE (kN)
5
STORY
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
MODAL AMPLITUDE
Figure B.1. Linear and non-linear 1st mode shapes of the 8 story wall building
5
STORY
0
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
MODAL AMPLITUDE
Figure B.2. Linear and non-linear 2nd mode shapes of the 8 story wall building
86
8
STORY
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
MODAL AMPLITUDE
Figure B.3. Linear and non-linear 1st mode shapes of the 12 story wall building
8
STORY
0
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
MODAL AMPLITUDE
Figure B.4. Linear and non-linear 2nd mode shapes of the 12 story wall building
87
12
10
STORY
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
MODAL AMPLITUDE
Figure B.5. Linear and non-linear 1st mode shapes of the 16 story wall building
12
10
STORY
0
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8
MODAL AMPLITUDE
Figure B.6. Linear and non-linear 2nd mode shapes of the 16 story wall building
88
16
14
12
STORY
10
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
MODAL AMPLITUDE
Figure B.7. Linear and non-linear 1st mode shapes of the 20 story wall building
16
14
12
STORY
10
0
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8
MODAL AMPLITUDE
Figure B.8. Linear and non-linear 2nd mode shapes of the 20 story wall building
89
20
STORY
15
10
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
MODAL AMPLITUDE
Figure B.9. Linear and non-linear 1st mode shapes of the 30 story wall building
20
STORY
15
10
0
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8
MODAL AMPLITUDE
Figure B.10. Linear and non-linear 2nd mode shapes of the 30 story wall building
90
Non-linear time history analysis is considered as the most advanced and most
complicated analysis method with respect to the simplified analysis methods used in
engineering practice. Simplified analysis methods have the advantage of specified ground
motion characteristics mostly in the form of response spectra either in the codified or site
specific forms. This eliminates the effort required by the designer to deal with the
variations arising from the ground motion side of the problem and modeling assumptions
and representations becomes the major concerns to deal with. However, reliable results
expected from a time history analysis are also dependent on appropriately selected strong
ground motion records since the ground motion characteristics will also affect the
response.
Two major issues are confronted with the strong ground motion records to be
incorporated in the non-linear time history analyses. First issue particularly comprises the
composition of a record suite with sufficient number of ground motions records having
comparable characteristics. It is evident that ground motion characteristics are primarily
affected by distance to fault, magnitude of the earthquake and local site conditions. Thus,
selected ground motion records should be representative of these three criteria where the
structure shall be designed. Theoretically, the more the number of records incorporated in
analyses, the more reliable results are obtained and an increase of number of samples allow
for working with average analysis results. It is obvious that finding sufficient number of
locally recorded accelerograms satisfying the three criteria is not an easy task, thus either
accelerograms with similar conditions can be selected or synthetically simulated
accelerograms can be developed.
It is intended that selected ground motion records suitably represent the expected
ground motion characteristics for Seismic Zone 1 defined in Turkish Seismic Design Code
(2007) in terms distance and magnitude. Pseudo-acceleration response spectra of the
selected records are given in Figure C.1 with respect to the design response spectrum
constructed for respective Site Class Z3 and the Seismic Zone 1.
92
16
DESIGN
12
10
0
0 0.5 1 1.5 2 2.5 3
PERIOD (sec)
Standardization of a suite of records is one of the key points for non-linear time
history analyses. Standardization is principally based on scaling of real accelerograms to fit
a target design ground motion level, which can be defined through design peak ground
acceleration, intensity or the response spectrum itself. Although details of this subject are
out of the scope of the present study, it should be reminded that any manipulation for the
sake of standardization has the risk of losing the natural characteristics of the real records.
Standardization of the selected records in this study has been achieved by scaling the
records to fit the target design response spectrum in all period ranges. An iterative
procedure has been applied by using a program coded in MATLAB with steps given
below:
1. Obtain target pseudo velocity response spectrum, for 5 per cent damping
ratio( S tpv ( f , ξ = 0.05) )
2. Obtain pseudo-velocity response spectrum of the original record, for 5 per cent
damping ratio ( S pv ( f , ξ = 0.05) ).
93
3. Divide the velocity response spectrum of the original record by the target pseudo-
velocity response spectrum to obtain the scale factors ( α ( f ) ) corresponding to
each frequency.
S pv ( f , ξ = 0.05)
α( f ) =
S tpv ( f , ξ = 0.05)
______
FFT (a(t )) = α ( f ) FFT (a(t ))
__ ______
a (t ) = inv[ FFT (a(t ))]
7. Repeat steps between 2 and 6 by replacing the Scaled Record in Step 6 with the
original record in step 2 until a satisfactory convergence is achieved between the
pseudo acceleration response spectrum of the record and the target pseudo
acceleration response spectrum.
Figure C.2 shows the pseudo-acceleration response spectrum of the original record
with respect to target acceleration response spectrum before iteration. Figure C.3, Figure
C.4 and Figure C.5 show the spectra at iteration steps 5, 10 and 20, respectively.
94
15
RECORD
10
0
0 0.5 1 1.5 2 2.5 3
PERIOD (sec)
Figure C.2 . Comparison of the record spectrum with target spectrum before iteration
15
RECORD
PSEUDO SPECTRAL ACCELERATION (m/sec2)
TARGET
10
0
0 0.5 1 1.5 2 2.5 3
PERIOD (sec)
Figure C.3. Comparison of the record spectrum with target spectrum at iteration step 5
95
15
RECORD
10
0
0 0.5 1 1.5 2 2.5 3
PERIOD (sec)
Figure C.4. Comparison of the record spectrum with target spectrum at iteration step 10
15
RECORD
PSEUDO SPECTRAL ACCELERATION (m/sec2)
TARGET
10
0
0 0.5 1 1.5 2 2.5 3
PERIOD (sec)
Figure C.5. Comparison of the record spectrum with target spectrum at iteration step 20
ORIGINAL RECORD
0.4
ACCELERATION (g)
0.2
-0.2
-0.4
0 5 10 15 20 25 30 35 40
time (sec)
SCALED RECORD
0.4
ACCELERATION (g)
0.2
-0.2
-0.4
0 5 10 15 20 25 30 35 40
time (sec)
Pseudo acceleration response spectra of the selected records after scaling are given in
Figure C.7 with respect to design acceleration response spectra. Accelerograms scaled by
using the methodology above has been used as the ground motion data in the non-linear
time history analyses of the generic walls as explained in Chapter 5.
97
15
PSEUDO SPECTRAL ACCELERATION (m/sec2) TARGET
SCALED REC.
10
0
0 0.5 1 1.5 2 2.5 3
PERIOD (sec)
Figure C.7. Comparison of the acceleration response spectrum of the records with the
target spectrum after scaling
98
Shear Design of the 16 Story Generic Wall Designed With a Strength Reduction
Factor of R=6
The proposed shear design procedure is presented below for a 16 story generic wall
designed with strength reduction factor of R = 6.
Key characteristics of the wall are summarized below;
N= 16 stories
b = 0.30 m (section width in flexure)
lw = 7.50 m (section depth in flexure)
h = 3.00 m (Typical story height)
Hw = 48.00 m (Total wall height measured from the base)
T1-cr = 1.30 sec (First mode period of the wall based on cracked section stiffness
properties)
Design shear forces have been obtained in Chapter 3 by using the response spectrum
method as shown in Figure D.1 below;
16
14
12
10
STORY
0
0 100 200 300 400 500 600 700 800
DESIGN STORY SHEAR FORCE (kN)
Veb = 741 kN
The base shear amplification factor is calculated by using the proposed relation given
in Equation 6.5 by assuming flexural overstrength of 1.0
M eb
= 1.0
M rb
The amplified design base shear is obtained from Equation 6.7 will give;
The story shear force diagram is obtained by using the suggested profile given
Section 6.2
• Determine the height of the critical wall height of constant shear right above the base
of wall
• Determine the level where the story shear forces drop to 50 per cent of the base shear
Amplified story design shear force diagram obtained as above is given in Figure D.2
where it is compared with the quantities obtained from design and Group 3 non-linear time
history analyses.
100
16
AMPLIFIED
14 DESIGN
T. HISTORY
12
10
STORY
0
0 500 1000 1500 2000 2500
STORY SHEAR FORCE (kN)
Figure D.2. Amplified design shear forces for the 16 story wall (R = 6)
101
Proposed Moment Diagram of the 16 Story Generic Wall Designed With a Strength
Reduction Factor of R=6
Design moment at the base of the 16 story generic wall has been obtained in Chapter
3 by using the response spectrum method.
M eb = 14677.34 kNm
Intermediate points on the proposed moment profile are obtained as per Section 6.4
by considering the tension shift effects at the base of the wall.
16
PROPOSED PROFILE
14 NL. TIME HIST.
DESIGN
12
10
STORY
0
0 5000 10000 15000
MOMENT
REFERENCES
Aydınoğlu, M.N., 2003, “An Incremental Response Spectrum Analysis Procedure Based
on Inelastic Spectral Displacements for Multi-Mode Performance Evaluation”,
Bulletin of Earthquake Engineering, Vol. 1, pp. 3-36.
Blakeley, G., C. Cooney and M. Megget, 1975, “Seismic Shear Loading at Flexural
Capacity in Cantilever Wall Structures”, Bulletin of the New Zealand National
Society for Earthquake Engineering, Vol. 8, No.4, pp. 278-290.
Bathe, K-J., 1996, Finite Element Procedures, Prentice Hall, Upper Saddle River, NJ,
USA.
Chopra, A.K., 2007, Dynamics of Structures, Prentice Hall, Upper Saddle River, NJ, USA.
Carr, A.J., 2000, “RUAUMOKO, Program for Inelastic Dynamic Analysis”, Department
of Civil Engineering, University of Canterbury, Christchurch, New Zealand.
CEN, 2004, Eurocode EC8, “Eurocode (EC) 8 Design of Structures for Earthquake
Resistance-Part I General Rules, Seismic Actions and Rules For Buildings (EN
1998-1)”, Comité Européen de Normalisation, Brussels.
CSI, 2006, “SAP 2000, Static and Dynamic Finite Element of Structures”, Computers and
Structures Inc., Berkeley, California, USA.
Eberhard, M. and M. Sözen, 1993, “Behavior Based Method to Determine Design Shear in
Earthquake-Resistant Walls”, Journal of the Structural Division, American Society of
Civil Engineers (ASCE), Vol. 119(2), pp. 619-640.
104
Eibl, J. and E. Keintzel, 1988, “Seismic Shear Forces in RC Cantilever Shear Walls”,
Proceedings of Ninth World Conference on Earthquake Engineering, Kyoto, Japan, 2
August-9 August 1991,Vol. VI, pp. VI/5-VI/10.
Filiatrault, A., D. D’Aronco and R. Tinawi, 1994, “Seismic Shear of Ductile Cantilever
Walls: A Canadian Code Perspective”, Canadian Journal of Civil Engineering, Vol.
21, pp. 363-376.
Krawinkler, H., 2006, “Importance of Good Nonlinear Analysis”, The Structural Design of
Tall and Special Buildings, Vol. 15, pp. 515-531.
NIED, 2006, “Test of a Six Story R/C Building”, Hyogo Earthquake Research Center, E-
Defense Facility, National Research Institute for Earth Science and Disaster
Prevention, Japan, www.bosai.go.jp.
Paulay, T. and M.J.N. Priestley, 1992, “Seismic Design of Reinforced Concrete and
Masonry Buildings”, John Wiley & Sons, New York, USA.
Rutenberg, A. and E. Nsieri, 2006, “The Seismic Shear of Ductile Cantilever Wall Systems
and the EC8 Provisions”, Bulletin of Earthquake Engineering, Vol. 4, pp.1-21.
Standards New Zealand, 1982, “NZS 4203 Code of Practice for Structural Design and
Design Loadings for Buildings”, Wellington, New Zealand.
Standards New Zealand, 2006, “Concrete Structures Standard: Part I-The Design of
Concrete Structures, Part II-Commentary on the Design of Concrete Structures”,
Wellington, New Zealand.
Takeda, T., M. Sözen and N. Nielsen, 1970, “Reinforced Concrete Response to Simulated
Earthquakes”, Journal of the Structural Division, American Society of Civil
Engineers (ASCE), Vol. 96(12), pp 101-113.
The MathWorks Inc., 2004, MATLAB, The Language of Technical Computing, USA.
Turkish Seismic Design Code, 2007, “Specification for Buildings to be Built on Seismic
Areas”, in Turkish, Ministry of Public Works and Settlement, Ankara, Turkey.
Turkish Seismic Design Code, 1997, “Specification for Structures to be Built on Disaster
Areas”, in English and Turkish, Ministry of Public Works and Settlement, Ankara,
Turkey.
UCSD, 2005, “Shake Table Test of a Full-Scale Seven Story Structural Wall”, Englekirk
Center, Powell Structural Research Laboratories, University of California, San
Diego, USA, www.jacobsschool.ucsd.edu/Englekirk.
106
Clough , R.W. and J. Penzien, 1993, Dynamics of Structures, 2nd International Edition,
McGraw Hill Inc., Singapore.
Kappos, A.J. and P. Antoniadis, 2007, “A Contribution to Seismic Shear Design of R/C
Walls in Dual Structures”, Bulletin of Earthquake Engineering, Vol. 5, pp. 443-466.
Linde, P., 1998, “Evaluation of Structural Walls Designed According to Eurocode 8 and
SIA 160”, Earthquake Engineering and Structural Dynamics, Vol. 27, pp. 793-809.
Priestley, M.J.N., 2003, “Myths and Fallacies in Earthquake Engineering, Revisited”, The
Mallet Milne Lecture, IUSS Press, Pavia, ITALY.
Riva, P., A. Meda and E. Giuriani, 2003, “Cyclic Behavior of a Full Scale RC Structural
Wall”, Engineering Structures, Vol. 25, pp. 835-845.