Reddy J.N. - Solution Manual For Introduction To Finite Element Analysis-MGH (2005)
Reddy J.N. - Solution Manual For Introduction To Finite Element Analysis-MGH (2005)
net/publication/278726673
CITATIONS READS
4 436
4 authors:
All content following this page was uploaded by Sarah Bobby on 16 August 2015.
ABSTRACT
Topology optimization has recently been explored as a tool for the conceptual design
of efficient structures during the initial building design stage. Although these
techniques were originally developed in a deterministic setting most practical
applications are subject to inherent uncertainties in the design problem (e.g. stochastic
loads, model idealization) and thus there has recently been an effort to incorporate the
effect of these uncertainties into the topology optimization procedure. Within this
context the authors have developed an approach that is able to efficiently account for
the uncertainties in the problem at hand using concepts from Performance-Based
Design. This paper presents a general overview of this method and illustrates its
applicability for the topology design of structures under wind and seismic hazards.
INTRODUCTION
λ ( dm ) = ∫ ∫ G ( dm | edp ) ⋅ dG ( edp | im ) ⋅ dλ ( im )
edp im
=η ⋅ P ( DM > dm )
where dm denotes a damage measure indicating damage, edp indicates an engineering
demand parameter characterizing the structural response, im indicates an intensity
measure describing the intensity of the hazard at the site of interest, λ (a) is the mean
annual rate of exceedance of event A=a, where capital letters indicate random
variables and lowercase letters indicate their realizations, G(a|b) = P(A > a | B = b)
denotes the complementary cumulative distribution function of random variables A
given B = b, f(a) is the probability density function of a, and η is the mean annual rate
of arrival of significant events. As can be seen in Eq. (2), the various terms are
independent of each other; these can be obtained through separate analyses (see
Figure 1): hazard analysis for λ(im), structural analysis for G(edp|im), and damage
analysis for G(dm|edp).
CASE STUDIES
Several case studies will now be given to illustrate the use of the PBTO methodology
in the design of novel, efficient layouts for conceptual building design. Two case
studies will be presented in this section. The objective of both studies will be to
obtain an efficient layout for the bracing of a planar lateral load resisting system
envisioned as part of a 3D building. The Method of Moving Asymptotes was used to
solve the sub-problems of the PBTO method (Svanberg 1987).
Case Study 1. The aim of the first case study is the topology design of the bracing for
a tall building located on the eastern coast of the United States near the city of Miami,
Florida. The objective is to minimize the volume of material used for the bracing
while meeting probabilistic performance constraints describing nonstructural damage.
The dominant hazard in this area is characterized by hurricane events; as such the
models described by Minciarelli et al. (2001) and Diniz et al. (2004) are used to
model the wind hazard. The external loads applied to the structure were determined
by integrating and appropriately scaling pressure measurements obtained from wind
tunnel tests performed at the boundary layer wind tunnel at the Inter-University
Research Centre on Building Aerodynamics and Wind Engineering (CRIACIV)
(Bernardini et al. 2012, Spence et al. 2011) and applying the model described in
Spence and Kareem (2014). The engineering demand parameters of interest were
taken to be the interstory drifts at all floors of the structure and the damage parameter
was taken as that corresponding to the occurrence of non-structural damage. The
nominal capacity for this value was taken as 1/400th of the floor-to-floor height. The
target damage measure was taken to be dmj = 1 for all floors. Two trials are
presented: Trial 1 considers a target mean annual rate of exceedance of λ0j =0.02 and
Trial 2 considers a target mean annual exceedance rate of λ0j =0.002. The initial
design domain is shown by the shaded gray area of Figure 3(a), where the horizontal
and vertical black lines represent frame elements modeling the underlying secondary
system (e.g. floor system) that is used to transfer local loads to the main lateral load
resisting system. Figures 3(b) and 3(c) show the optimal bracings for Trials 1 and 2
respectively. Although a similar high-waisted X-bracing pattern is noticed for both
trials, Trial 2 shows noticeably thicker bracing due to the stricter requirement on the
target mean annual exceedance rate. Figure 4 illustrates trends regarding the
objective and constraint functions. The material volume history given in Figure 4(a)
shows rapid and steady convergence for both trials and illustrates the efficiency of the
method. Figure 4(b) plots the maximum mean annual exceedance rate over all
constraints for each design cycle, and it can be seen that at the final design cycle the
required target mean annual exceedance rate is met for both trials. Finally, Figure
4(c) gives the damage measures with target mean annual exceedance rate for the
initial and final structures for both trials. A damage measure of dm≤1 is shown for
the final structures at all stories, indicating that all constraints were met.
Figure 3: Case Study 1: (a) Initial design domain; Optimal topologies for
(b) Trial 1, (c) Trial 2
Figure 4: Case Study 1 (a) Volume history; (b) Maximum mean annual
exceedance rate over all constraints; (c) Damage measure with target mean
annual exceedance rate
Case Study 2. The second case study considers the topology design of a bracing
system for a three story building under earthquake loading. It is assumed that the
structure is located in California and the earthquake hazard is represented using the
models suggested by Boore (2003) and Vetter and Taflanidis (2012, 2014). As in the
first case study, the engineering demand parameters of interest are taken to be the
interstory drifts at all floors of the structure and the damage parameter was taken as
that corresponding to the occurrence of non-structural damage with a nominal
capacity of 1/400th of the floor-to-floor height. The damage measure was again taken
to be dmj = 1 for all floors. The target mean annual rate of exceedance was taken to
be λ0j =0.01, and earthquake events with magnitudes between 5 and 8 were
considered. The initial design domain is shown by the shaded gray area of Figure
5(a), where the secondary system of the structure is illustrated using horizontal and
vertical black lines representing frame elements. Figure 5(b) shows the optimal
bracing scheme after the problem is solved using the PBTO solution strategy. The
PBTO methodology produced an X-bracing scheme, which is encouraging because an
X-bracing scheme was also determined by Chun et al. (2012) for a structure having a
design domain with the same dimensions. Figure 6 illustrates trends regarding the
objective function and constraints for Case Study 2. The material volume history
given in Figure 6(a) shows rapid and steady convergence, which was achieved after
only 6 design cycles and illustrates the efficiency of the PBTO solution strategy. The
maximum mean annual exceedance rate over all constraints is given in Figure 6(b).
As can be seen, the target mean annual exceedance rate was met for the final design
at design cycle 27. Finally, Figure 6(c) gives the damage measures with target mean
annual exceedance rate for the initial and final structures of Case Study 2. The final
structure displays a damage measure of 1 for the constraints at all stories.
Figure 5: Case Study 2: (a) Initial design domain; (b) Optimal Topology
Figure 6: Case Study 2 (a) Volume history; (b) Maximum mean annual
exceedance rate over all constraints; (c) Damage measure with target mean
annual exceedance rate
CONCLUSIONS
ACKNOWLEDGEMENTS
Support for this work was provided by the NSF Grant No. CMMI-1301008. The
authors also thank Krister Svanberg for providing the MMA algorithm, which was
used as the optimization algorithm in this research.
REFERENCES