Research Article  Open Access
Nonlinear Earthquake Analysis of Reinforced Concrete Frames with Fiber and BernoulliEuler BeamColumn Element
Abstract
A beamcolumn element based on the EulerBernoulli beam theory is researched for nonlinear dynamic analysis of reinforced concrete (RC) structural element. Stiffness matrix of this element is obtained by using rigidity method. A solution technique that included nonlinear dynamic substructure procedure is developed for dynamic analyses of RC frames. A predictedcorrected form of the Bossakα method is applied for dynamic integration scheme. A comparison of experimental data of a RC column element with numerical results, obtained from proposed solution technique, is studied for verification the numerical solutions. Furthermore, nonlinear cyclic analysis results of a portal reinforced concrete frame are achieved for comparing the proposed solution technique with Fibre element, based on flexibility method. However, seismic damage analyses of an 8story RC frame structure with softstory are investigated for cases of lumped/distributed mass and load. Damage region, propagation, and intensities according to both approaches are researched.
1. Introduction
The realistic modeling of the nonlinear static or dynamic behavior of RC structures is a more sophisticated problem due to the inelastic behavior of concrete, plasticity of reinforcement, interface debonding of these materials, and so on. In the last thirty years, researchers have made more effort for the numerical modeling of RC structures and most of the stateoftheart on this problem deals with two main approaches: lumped plasticity modeling [1–3] and distributedinelasticity modeling (i.e., the socalled fibre beamcolumn elements, FBCE) [4–8]. In the first approach, nonlinear springs based on momentrotation and forcedisplacement curves are used and nonlinear volume is assumed to be lumped in the specific location of the element. Mohr et al. [1] were used a series of polynomial shape functions for strain distributions of the vertical and shear on crosssection. Li et al. [2] was proposed a simplified lumped hinge for determination of failure mode of statically indeterminate structure. Reshotkina and Lau [3] was used a concentrated plasticity approach for axialflexureshear interactions on the inelastic behavior of reinforced concrete members under the seismic loads. In the second approach, the nonlinear behaviors of concrete and reinforcement materials are separately calculated. This method divided into two as approaches was based on flexibility and rigidity. Taucer et al. [9] were developed a flexibility method for nonlinear dynamic analysis of reinforced concrete element. This method based on BernoulliEuler hypothesis included biaxial bending and axial force conditions. Nonlinear behavior of concrete and steel were computed by using uniaxial stressstrain relationships. Forcedisplacement interpolation functions are used for the obtaining of element flexibility matrix. Element stiffness matrix is achieved by inverse of flexibility matrix. Furthermore, Ceresa et al. [10] were developed a beamcolumn element based on flexibility matrix by using to the Timoshenko beam theory under cyclic loading. Damages in the element are taken into account with coaxial rotating crack model. Lu et al. [4] were used multilayered shell element and fibre beamcolumn element based on flexibility method for modeling of the different frame and frame with shear wall. However, an important research report about methods based on rigidity was performed by Kawano et al. [6]. In the report, nonlinear analyses of reinforced concrete space frames with multistory and multibay are obtained under the earthquake loads. Cubic Hermitian polynomials for transversal displacements and the linear shape functions for axial displacements are used for obtaining of stiffness matrix. However, the GaussLobatto integration rule was used for obtaining element mass and stiffness matrices. Element damage location was simply presented and element was not divided into subelement called as “segment.” Some calculation problems in heavily damage location was shown and they are reported to arise from flexural deformation. Furthermore, strain rate effects in response to reinforced concrete frames was investigated by Iribarren [7]. Used formulations are based on rigidity method and a strain rate dependent material formulation is developed for both the concrete and steel constitutive response. Brum [8] was obtained nonlinear dynamic analyses of frame and masonry structures by using a method which was based on stiffness matrix. In the research, a uniaxial constitutive model for concrete and masonry is proposed for compression and tension regions of the materials under the cyclic loadings.
In this study, a beamcolumn element based on the EulerBernoulli beam theory is presented for the fiber RC element. Element stiffness matrices are obtained by using rigidity method. The beam or column element is divided into subelements called “segment” [11]. Furthermore, the internal freedoms of this segment are dynamically condensed to external freedoms at the end of the element. Thus, nonlinear dynamic analysis of high RC building can be obtained within a short time. This procedure requires that nested loops of the element and structure are obtained at the same time (the socalled nonlinear dynamic substructure). However, this condensation procedure is not used in the modeling of the fibre element approach (FEA) of the RC element. In addition, uniform or trapezoidal loads of the segment are assumed to be zero or condensed to the external freedoms at the end of the element in the FEA [11]. This case is not preferred due to the required redistribution of loading. However, in this study, external loading of the segment is taken into account by considering the damage occurred in the element.
The present research is organized as follows: (i) presenting of fiber BernoulliEuler beamcolumn element based on rigidity method, (ii) developing of nonlinear dynamic substructure technique of RC frames, (iii) verification of the constitutive model with respect to experimental result of a column structural element and dynamic analysis results of fibre RC element of a portal frame, (iv) obtaining of seismic damage analyses of an 8story RC frame with softstory for distributed/lumped mass and load case, and (v) results.
2. Fiber and BernoulliEuler Approach (FBEA) for Reinforced Concrete Beam Column Element
In this section, theory of Fiber and BernoulliEuler element which is based on stiffness matrix is firstly presented for reinforced concrete section. Elements are divided into subelements, called segment, and they are assumed as substructures of element. Furthermore, crosssection of the each segment is also subdivided into a number of fibers/layers. In the next section, nonlinear dynamic substructures method is also mentioned to be applied to freedoms in element ends.
2.1. Obtaining Segment Stiffness and Mass Matrices
The strain distribution on a crosssection of beam is assumed to be uniform due to axial forces and linear due to only bending according to the BernoulliEuler approach. Furthermore, if a plane section before bending is plane after bending and if the strain is assumed to be small and shear stresses are omitted, the strain of a point on a crosssection in the axial direction can be written as where and are the displacement of the axial and vertical directions of element, respectively (Figure 1). However, if crosssection of element is divided into fibers/layers, this equation for each fiber/layer in local axis can be rewritten as, where is incremental axial strain of a fiber/layer on a segment in local axis direction. Therefore, if the cubic Hermitian polynomials for transversal displacements and the linear shape functions for axial displacements are used [12], (2) can be written as where is displacement vector, including displacement and rotation of a segment. is straindisplacement matrix. Thus, incremental stress in each fiber/layer can obtained as where , determined by using uniaxial stressstrain relationship of using materials, is tangent elasticity modulus of each fiber. Therefore, total strain energy of a segment subelement is obtained as where and are expressed to be length and crosssection area of segment, respectively. is also area of a fiber/layer. Thus, element stiffness matrix can be obtained by using minimum potential energy principle. Stiffness matrix of a segment can be written as The segment stiffness matrix is obtained by using areas, coordinates, and tangent elasticity modulus of fiber/layer for reinforced concrete sections. In this study, linear superposition rule is used for different material properties of the fiber/layer in the section (Figure 1). Furthermore, a fiber concrete or reinforced bar on the crosssection may be cracked or damaged due to external loading. For this reason, a relationship between the damaged and undamaged cases must be obtained for the solutions. Thus, damage in each concrete/reinforcement fiber can be written as [13] where and are undamaged and damaged/tangent elasticity module of the th fiber, respectively. , damage intensities are obtained separately under tensile and compressive stress for each fiber.
(a)
(b)
However, total kinetic energy of particle velocities on the crosssection of an element throughout its neutral axis can be written as where is mass density and is particle velocity. If fiber/layer elements are used for (8), the mass matrices of an element on the local axis are obtained as where is the element shape functions matrix on the local axis and is also mass density of the th fiber. Stiffness and mass matrices of a segment, obtained by local axis, are transformed to global axes by using transformation matrix.
2.2. Obtaining with Nonlinear Dynamic Substructures Technique of Element Stiffness and Mass Matrices
In FBEA, cubic Hermitian shape functions for rotation and shear strain and linear shape functions for the axial deformation are used, respectively. The GaussLobatto integration rule is generally used for obtaining element mass and stiffness matrices along with local axis [6–8]. However, substructure procedures are generally preferred for solutions of high buildings. However, this technique in the previous study is not used due to numerical difficultly. Element stiffness/flexibility matrix is obtained by using crosssection properties on the integration points. Therefore, effects of shape functions on the solution are very important. In this study, an element is divided into subelements which are called “segment” to remove disadvantage effects of the shape functions. The freedoms of segment are condensed to the end freedoms at the ends of the element. Thus, if freedoms at ends and internal regions of element are called as external and internal freedoms, respectively (Figure 2), these external and internal freedoms for stiffness matrices can be written as,where and are stiffness matrices which include external and internal freedoms, respectively [11]. and are the displacement vectors referred by these freedom; and are also the external loads of these freedoms. When applying the substructure procedure to (10a) and (10b), the external load and stiffness matrices of a frame element can be rewritten asHowever, if the Rayleigh method may be used to obtain element damping matrices, this matrix is defined as where and are Rayleigh damping coefficients with respect to mass and stiffness matrices, respectively [14]. However, if substructure procedure in (11a), (11b), and (11c) is applied to mass and damping matrices [15],mass and damping matrices can be obtained for the external freedom. Global stiffness, damping, and mass matrices may be achieved by using the stiffness and mass matrices belonging to external freedom. Thus, global stiffness, mass matrices, and external load vector can be calculated as Therefore, the dynamic equilibrium equations of the structure may be given as where subscripts and indicate that the quantity is related to ground acceleration and static loads.
2.3. Bossakα Form of the Equation of Motion
The equation of motion for the RC frame is given in (15). In this study, the Bossakα integration method, presented by Wood et al. [16], is used for the solution of the equation in the time domain. Integration scheme of the method retains the Newmark method. Furthermore, the Bossakα integration method is required to be modified to (15) in the time domain as follows: where is the Bossak parameter, used for controlling the numerical dissipation. The Bossak parameter should be chosen as in (17) for unconditional stability and secondorder accuracy: In this study, is selected as −0.10. To solve the nonlinear dynamic equation of motion for the RC frame, the NewtonRaphson method is used in conjunction with the predictorcorrector technique. The Bossakα time integration algorithm is given by Wood et al. [16]. Predicted displacement and velocity vectors for the time step are obtained by using displacement and velocity vectors at the time step which is known. Thus, they can be calculated aswhere and are Newmark’s coefficients. The displacement and velocity vectors [17] of the RC frame can be written in terms of the predicted vectors shown in (18a) and (18b). The vectors may be defined asThese relations can be substituted into (16) and a time marching algorithm can be applied to this equation as given in the Appendix.
3. Numerical Applications
3.1. Comparison of Experimental and Numerical Analysis of a RC Column
In this section, experimental cyclic test results of a RC column with numerical solutions obtained from proposed solution method are compared. The experimental test result includes response under uniaxial static and lateral cyclic forces. Loading and material properties of the experimental study are given by Takahashi [18]. The experimental setup is shown in Figure 3. Exponential decreasing functions for the softening region of tensile and compressive strengths of the concrete for the constitutive model of FBEA are used. These functions are shown in Figure 4(a). The bilinear kinematic hardening rule is used for nonlinear behavior of the steel for the two approaches (Figure 4(b)). Static loads are converted to masses which are condensed to the element ends for all solutions and displacement values obtained due to the loads being considered as the initial condition. Acceleration data, a sinus wave shape, and time varying are used for all dynamic solutions. This dynamic load is applied to the horizontal direction at 1280 mm height of RC column. Tangent stiffness matrix is used for the solution and damping matrix is also assumed to be proportional to stiffness matrix.
(a)
(b)
Loaddisplacement curves of experimental and numerical results are given in Figure 5. Maximum and minimum values of cyclic displacement obtained from numerical results are approximately between −0.055 and 0.055 m. Numerical response of RC element which is obtained with FBEA is shown as similar to envelope curve of experimental results. All values of the horizontal displacements are on the envelope curve of the experimental result. It is said that this solution technique and material models of concrete and steel can be used for the solution of the RC structural element under the dynamic loading.
3.2. Nonlinear Dynamic Analyses of a Portal RC Frame Structure
In this section, nonlinear dynamic analyses of a portal RC frame are obtained for the comparing of proposed solution technique with model based on flexibility. SeismoStructure program [19] is used for the flexibility based solutions (with fibre element method). Finite element meshes are shown in Figure 6 for both approaches. Finite element mesh of proposed method is divided into segment for comparing with analysis results of SeismoStructure Program. ACI 31802 [20] code is used for the material properties of concrete and steel. Crosssection and material properties of the beam and column of selected portal RC frame are given in Table 1. The tensile softening region in the fibre element method is not taken into account while the compressive behavior of the concrete is used. In SeismoStructure program, if the tensile stresses obtained from the solutions reach the tensile strength of the concrete, tensile strength suddenly drops to zero. In this case, a sudden collapse of the structure occurs.

(a)
(b)
In this study, exponential decreasing functions in the softening region of tensile and compressive strengths of the concrete for constitutive model of FBEA are used. These functions are shown in Figure 4(a). The bilinear kinematic hardening rule is used for nonlinear behavior of the steel for the two approaches (Figure 4(b)). Static loads are converted to masses which are condensed to the element ends for all solutions and displacement values obtained due to the loads being considered as the initial condition. Acceleration data and a sinus wave shape are used for all dynamic solutions (Figure 7). This dynamic load is applied to the horizontal direction.
Displacement time history graphs of node 3 obtained from FBEA and FEA are shown in Figure 8. Displacement time amplitude values are approximately similar until time of 0.75 sec and, after this time, displacement amplitude values obtained from the FEA are bigger than those obtained from the FBEA. This case arises from not taking into account the concrete tensile softening region for the fibre element method. However, failure of structure is not seen for two approaches in all times.
Accumulated tensile damage cases obtained for both approaches are given in Figure 9. The first damage zones are obtained at the whole of end regions of the beams and at outside surfaces of bottom regions of the columns for the FBEA. Damages are shown at the whole crosssection of bottom regions of the columns for the FEA. Furthermore, obtained damage zone regions in the beam and column for the FBEA and FEA are similar at sec. However, some differences between both approaches are seen from the point of intensities of the damage at this time. However, damage intensities at end regions of the beam according to FBEA are extended and damage regions are propagated to middle region of the beam at the sec. Damage intensities at end regions of the column are also extended and damage regions are propagated to middle region from upper region of the column at the sec. Additionally, damage intensities at end regions of the beam according to FEA are not extended at this time. Damage intensities at end regions of the column are also extended and damage regions are propagated to middle region from upper region of the column; after this time, damage zones are extended but propagation of damages are not seen according to both approaches at the sec.
3.3. Seismic Damage Analyses of an 8Story RC Frame Structure with SoftStory
In this numerical application, nonlinear dynamic analyses of an 8story RC frame structure with softstory are investigated for cases of lumped/distributed mass and load. ACI 31802 [20] code is used for the material properties of concrete and steel. Crosssection and material properties of the beam and column of selected RC frame are given in Table 2. Finite element mesh and gravity loading case are shown in Figure 10. Uniaxial stressstain relationships of the concrete and steel are seen in Figures 4(a) and 4(b), respectively. Static loads and displacements are considered as initial conditions in all solutions. Spectrum acceleration curve, given by Z1 type soil in the Turkish Regulation Code on Building in Disaster Areas [20], is selected as the target spectrum curve. Synthetic earthquake acceleration data for maximum amplitude, 0.3 g, are produced. This synthetic accelerationtime graph is shown in Figure 11 and it is affected on the horizontal direction of the RC frame structure. Tangent stiffness matrix is used for the all solution and damping matrix is also assumed to be proportional to stiffness matrix.

Displacement time history graphs of node 9 obtained from nonlinear dynamic analyses results for lumped/distributed mass and load cases are shown in Figure 12. Displacement amplitude values are approximately similar until time of 1.08 sec and, after this time, displacement amplitude values obtained from the distributed mass and load case are bigger than those obtained from the lumped case. Absolute maximum horizontal displacements for the lumped and distributed cases are obtained as 22.6 and 41.4 mm, respectively. Thus, absolute maximum displacement value according to distributed case is approximately 83% of ratio bigger than that of lumped case. This result arises from the redistribution of loading in the internal regions of the elements for the distributed approach. However, the redistribution of loading in these internal regions is only obtained at the end regions of elements for the lumped approach. The redistribution procedure for loading in these internal regions requires more iterative steps for a distributed approach. However, numerical dissipation is not shown for both solutions and NewtonRaphson procedures for all element solutions which are converged for the two approaches. The Bossakα dynamic integration algorithm is applied successfully to the nonlinear dynamic solutions.
Accumulated tensile damage regions obtained for both mass/load case are given in Figures 13 and 14. Damage zones in the lumped case are obtained at the end of the beam of all floors and at the whole of the beam crosssection for the time = 1.00 sec. Damage zones occurred at upper parts of the beams at the beamcolumn join region and at lower parts of the mid region of the beams until this time for the distributed case. Furthermore, damage zones seen at lower parts of firstfloor columns occurred for both approaches, but intensities and regions of these damages according to distributed approach are bigger than those of lumped mass/load case. The damage intensities are some increased and damage propagations are remained at the same regions for both approaches until sec; after this time, damage intensities according to two approaches are increased at between and 8.00 sec, and no change for the damage zones is obtained. However, increases in the damage intensities achieved a minimal level between and 10.00 sec. It is said that the damage zones obtained according to the distributed/lumped case had important differences.
(a) sec
(b) sec
(c) sec
(d) sec
(e) sec
(f) sec
(g) sec
(h) sec
(a) sec
(b) sec
(c) sec
(d) sec
(e) sec
(f) sec
(g) sec
(h) sec
4. Conclusions
In this study, a beamcolumn element based on the EulerBernoulli beam theory is researched for nonlinear dynamic analysis of reinforced concrete (RC) structural element. Stiffness matrix of this element is obtained by using rigidity method. A solution technique that included nonlinear dynamic substructure procedure is developed for dynamic analyses of RC frames. A predictedcorrected form of the Bossakα method is applied to dynamic integration scheme. A comparison of experimental data of a RC column element with numerical results, obtained from proposed solution technique, is studied for verification of the numerical solutions. Furthermore, nonlinear cyclic analysis results of a portal reinforced concrete frame are achieved for comparing of the proposed solution technique with fibre element, based on flexibility method. However, seismic damage analyses of an 8story RC frame structure with softstory are investigated for cases of lumped/distributed mass and load. Damage region, propagation, and intensities according to both approaches are researched. Results obtained from this study are presented to be itemized as follow.(i)Constitutive model based on rigidity method is used in obtaining the stiffness matrix. The beam or column element is divided into a subelement called the segment. The internal freedoms of this segment are dynamically condensed to the external freedoms at the ends of the element. FBEA requires less freedom than the FEA due to the nonlinear dynamic substructure. Thus, nonlinear dynamic analysis of high RC building can be obtained within short times. This condensation procedure in previous study is not used for the modeling of the FBEA of RC element.(ii)Numerical response of RC element which obtained with FBEA is similar to envelope curve of experimental results. All values of the horizontal displacements are on the envelope curve of the experimental result. It is said that this solution technique and material models of concrete and steel can be used for the solution of the RC structural element under the dynamic loading.(iii)Accumulated tensile damage cases obtained for both FBEA and FEA approaches are similar, but some differents for displacement values obtained according to both methods are seen. These differents are to be small for little values of tensile damage intensities and propagation regions of damage and are to be big depending on increasing of the damage intensities and propagation regions of damage. However, tensile damage regions are obtained in detail according to the FBEA.(iv)When the results obtained from lumped and distributed approaches are investigated, absolute maximum displacement values for distributed approach are seen to be bigger than those for lumped approach. This case arises from the redistribution of loads in the internal regions of the elements for the distributed approach. However, the redistribution of loads in these internal regions is not used for the lumped approach. The redistribution procedure for loading in these internal regions requires more iterative steps.(v)Obtained damage zones for the lumped mass/load case are seen at the end of the beam of all floors and at the whole of the beam crosssection, but obtained damage regions for the distributed mass/load case showed up at upper parts of the beams at the beamcolumn join region and at lower parts of the mid region of the beams. Furthermore, additional damage zones at lower parts of firstfloor columns occurred for both approaches, but intensities and regions of these damages according to distributed approach are bigger than those of lumped approach.(vi)However, numerical dissipation is not shown for all solutions and NewtonRaphson procedures for all element solutions are converged for the two approaches. The Bossak time integration algorithm and proposed solution technique are successfully applied to the nonlinear dynamic solutions.
Appendix
Time Marching Algorithm
The algorithm applied to (16) is as follows.(1)Compute integration parameters (2), , and are known; set global iteration counter, .(3)Predict response at ,(4)Set element counter, .(5)Set element iteration counter, . Subtract incremental displacement vector of external joints from global displacement vectors, (6)Compute the element stiffness, mass, damping matrices, and element external and internal force vectors.(7)Compute iterative and incremental displacement vectors of internal joints,(8)Check for convergence of iteration process (Euclidian norm of the unbalanced internal displacement vector at time step and element iteration ) using an element displacement tolerance (Tolelem).(a)If convergence is achieved. Set(b)If global convergence is not achieved, set and return to Step 5.(9)Compute the effective dynamic stiffness matrix and the vector of unbalanced forces,where and are the tangent stiffness matrix and the incremental restoring force vector of external joints on the global axis of the structure, respectively, which are assembled from the element contributions.(10)Solve incremental displacements in the global axis, (11)Update displacement, velocity, and acceleration vectors,(12)Check for convergence of the iteration process (Euclidian norm of the unbalanced force vector at time step and global iteration ) using a force tolerance (Tolglo).(a)If , convergence is achieved. Set , , and .(b)If global convergence is not achieved, set and return to Step 4.
Abbreviations
FEA:  Fibre element approach 
FBEA:  Fiber and BernoulliEuler approach 
RC:  Reinforced concrete. 
:  Area of the th fiber/layer 
:  Acceleration vector of structure 
:  Crosssection area of segment 
:  Straindisplacement matrix 
:  Damping matrix for external freedoms 
:  Frame damping matrix for external freedoms 
:  Damping matrix for externalinternal freedoms 
:  Damping matrix for internalexternal freedoms 
:  Damping matrix for internal freedoms 
:  Initial elasticity module of the th fiber/layer 
:  Tangent elasticity module of the th fiber/layer 
:  External load vector for external freedoms 
:  External load vector for internal freedoms 
:  Dynamic external force vector for internal freedoms 
:  Dynamic restoring force vector for internal freedoms 
:  Dynamic external force vector of structure 
:  Dynamic restoring force vector of structure 
:  Static external force vector of structure 
:  Stiffness matrix for external freedoms 
:  Frame stiffness matrix for external freedoms 
:  Stiffness matrix for externalinternal freedoms 
:  Stiffness matrix for internalexternal freedoms 
:  Stiffness matrix for internal freedoms 
:  Stiffness matrix of the segment 
:  Length of segment 
:  Mass matrix for external freedoms 
Frame mass matrix for external freedoms  
:  Mass matrix for externalinternal freedoms 
:  Mass matrix for internalexternal freedoms 
:  Mass matrix for internal freedoms 
:  Mass matrix of the segment 
:  Number of total fiber/layer on the fiber crosssection 
:  The element shape functions matrix 
:  Displacement vector 
:  Total kinetic energy of particle velocities on the crosssection of an element throughout its neutral axis 
Tolelem:  Tolerance value of element for NewtonRaphson procedure 
Tolglo:  Tolerance value of structure for NewtonRaphson procedure 
:  Displacement in the axial directions of element local axis 
:  Displacement vector for external freedoms 
:  Displacement vector for internal freedoms 
:  Displacement vector of structure 
:  Predicted displacement vector of structure 
:  Displacement in the vertical directions of element local axis 
Velocity vector of structure  
:  Predicted velocity vector of structure 
:  Particle velocity 
:  Bossak parameter 
:  Rayleigh damping coefficient for mass matrix 
:  First Newmark’s coefficient 
:  Rayleigh damping coefficient for stiffness matrix 
:  Strain of a point on a crosssection in the axial direction of element 
:  Incremental axial strain of the th fiber/layer in local axis direction of a segment 
:  Second Newmark’s coefficient 
:  Vertical coordinates of the th fiber/layer in local axis direction of a segment 
:  Mass density 
:  Total strain energy of a segment. 
Conflict of interests
The author declares that there is no conflict of interests regarding the publication of this paper.
References
 S. Mohr, J. M. Bairán, and A. R. Marí, “A frame element model for the analysis of reinforced concrete structures under shear and bending,” Engineering Structures, vol. 32, no. 12, pp. 3936–3954, 2010. View at: Publisher Site  Google Scholar
 Y. Li, X. Lu, H. Guan, and L. Ye, “An improved tie force method for progressive collapse resistance design of reinforced concrete frame structures,” Engineering Structures, vol. 33, no. 10, pp. 2931–2942, 2011. View at: Publisher Site  Google Scholar
 S. S. Reshotkina and D. T. Lau, “Modeling damagebased degradations in stiffness and strength in the postpeak behaviour in seismic progressive collapse of reinforced concrete structures,” in Proceedings of the 15th World Conference on Earthquake Engineering, Lisboa, Portugal, September 2012. View at: Google Scholar
 X. Lu, X. Lu, H. Guan, and L. Ye, “Collapse simulation of reinforced concrete highrise building induced by extreme earthquakes,” Earthquake Engineering and Structural Dynamics, vol. 42, no. 5, pp. 705–7723, 2012. View at: Google Scholar
 OpenSees, “Open system for earthquake engineering simulation,” Pacific Earthquake Engineering Research Center, University of California, Berkeley, Calif, USA, 2013, http://opensees.berkeley.edu. View at: Google Scholar
 A. Kawano, M. C. Griffith, H. R. Joshi, and R. F. Warner, “Analysis of behaviour and collapse of concrete frames subjected to severe ground motion,” Research Report No. R163, Department of Civil and Environmental Engineering, Adelaide University, South Australia, Australia, 1998. View at: Google Scholar
 B. S. Iribarren, Progressive collapse simulation of reinforced concrete structures: influence of design and material parameters and investigation of the strain rate effects, [Ph.D. thesis], Polytechnic Faculty, Faculty of Applied Sciences, Bruxelles Royal Military Academy , Université Libre de, 2010.
 J. F. S. Brum, A model for the non linear dynamic analysis of reinforced concrete and masonry framed structures [Ph.D. thesis], Universitat Politècnica de Catalunya, 2010.
 F. F. Taucer, E. Spacone, and F. C. Filippou, “A Fibre beamcolumn element for seismic response analysis of reinforced concrete structures,” EERC Report91/17, Earthquake Engineering Research Center, Berkeley, Calif, USA, 1991. View at: Google Scholar
 P. Ceresa, L. Petrini, R. Pinho, and R. Sousa, “A fibre flexureshear model for seismic analysis of RCframed structures,” Earthquake Engineering and Structural Dynamics, vol. 38, no. 5, pp. 565–586, 2009. View at: Publisher Site  Google Scholar
 S. Krishnan, “Threedimensional nonlinear analysis of tall irregular steel buildings subject to strong ground motion,” EERL200301, Earthquake Engineering Research Laboratory, California Institute of Technology, Berkeley, Calif, USA, 2003. View at: Google Scholar
 T. R. Chandrupatla and A. D. Belegundu, Introduction to Finite Elements in Engineering, Prentice Hall, New Jersey, NJ, USA, 2012.
 F. Légeron, P. Paultre, and J. Mazars, “Damage mechanics modeling of nonlinear seismic behavior of concrete structures,” Journal of Structural Engineering, vol. 131, no. 6, pp. 946–955, 2005. View at: Publisher Site  Google Scholar
 K. J. Bathe, Finite Element Procedures in Engineering Analysis, Prentice hall, Englewood Cliffs, NJ, USA, 1982.
 R. J. Guyan, “Reduction of stiffness and mass matrices,” AIAA Journal, vol. 3, no. 2, article 380, 1965. View at: Google Scholar
 W. L. Wood, M. Bossak, and O. C. Zienkiewicz, “A alpha modification of Newmark’s method,” International Journal for Numerical Methods in Engineering, vol. 15, pp. 1562–1566, 1981. View at: Google Scholar
 I. Miranda, R. M. Ferencz, and T. J. R. Hughes, “An improved implicitexplicit time integration method for structural dynamics,” Earthquake Engineering and Structural Dynamics, vol. 18, no. 5, pp. 643–653, 1989. View at: Google Scholar
 Y. Takahashi, Development of high seismic performance RC piers with objectoriented structural analysis [Ph.D. thesis], Faculty of Engineering of Kyoto University, 2002.
 “Earthquake engineering software solution”, SeismoStruc Ver 6, 2013, http://www.seismosoft.com.
 ACI, “Building code requirements for structural concrete,” ACI 31802, American Concrete Institute, Farmington Hills, Mich, USA, 1995. View at: Google Scholar
Copyright
Copyright © 2014 Muhammet Karaton. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.