Failure behavior modeling of slender reinforced concrete columns subjected to eccentric load 1

This work presents a numerical model to simulate the failure behavior of slender reinforced concrete columns subjected to eccentric compression loads. Due to the significant influence of the lateral displacements on the loading state provided by an eccentric load, geometric nonlinearity is considered. The responses of the concrete in tension and compression are described by two scalar damage variables that reduce, respectively, the positive and negative effective stress tensors, which lead to two different damage surfaces that control the dimension of the elastic domain. To describe the behavior of the reinforcements, truss finite elements with elastoplastic material model are employed. Interaction between the steel bars and concrete is modeled through the use of interface finite elements with high aspect ratio and a damage model designed to describe the bond-slip behavior. The results showed that the numerical model is able to represent the nonlinear behavior of slender concrete columns with good accuracy, taking into account: formation of cracks; steel yielding; crushing of the concrete in the compressive region; and interaction between rebars and concrete.


Latin American Journal of Solids and
After the advent of the high strength concrete (HSC) a large number of slender RC columns has been constructed, and more researchers have given special attention to this subject (Pallarés et al., 2008;Kim and Yang, 1995).Consequently, the development of numerical tools to support the analyses has increased and different numerical models to simulate the failure mechanism have been proposed.For instance, in order to investigate the behavior of RC columns under eccentric compression load, with consideration of cracks, Majewski et al. (2008) proposed an approach based on elasto-plasticity with non-local softening model.Recently, Carmona et al. (2013) developed a fracture mechanics model to describe the buckling behavior of lightly concrete columns, wherein the structural response is a result of the propagation of one single crack at the column's central region.Using the finite difference method, Bouchaboub and Samai (2013) proposed an approach to simulate the behavior of slender high strength RC columns under combined biaxial and axial compression.For all the aforementioned works, the numerical tools developed were very useful to support the nonlinear analyses and to simulate the failure mechanism presented by RC columns.
Therefore, seeking to contribute for a better understanding of the failure mechanisms of RC columns, a numerical model composed by some features which may be easily implemented in an existent finite element code is proposed.These features are made up by material models based on the Continuum Damage Mechanics Theory (CDMT) to represent the concrete and steel-concrete interface and an elastic-perfectly plastic model to simulate the steel reinforcement behavior.Moreover, the applicability of the Interface Finite Element (IFE) with very high aspect ratio proposed by Manzoli et al. (2012) is investigated for the type of problem addressed in this paper.The IFE is responsible to connect the steel rebars on the concrete matrix and, in combination with the use of an appropriate constitutive model, it is able to describe the complex bond-slip phenomenon.
Initially, to calibrate the parameters of the material models, simple compression, tension and pullout tests are carried out.After that, to assess the effectiveness of the proposed numerical approach, a number of columns subjected to eccentric compressive load are numerically analyzed.In this study, the obtained results are compared to the experimental and numerical results reported by Kim and Yang (1995); Majewski et al. (2008), respectively.The considered variables were the slenderness ratio ( ) and the concrete strength ( 25.5MPa , 63.5MPa and 86.2MPa ) for the same steel ratio of 1.98% ρ = .This paper is organized in five main sections.In section 2 the constitutive models adopted to represent the concrete, steel reinforcements and steel-concrete interface are described.In section 3, the scheme adopted to connect the steel bars on the concrete matrix through the use of interface finite elements is presented.Later, in section 4, the features of the numerical model proposed for modeling the failure process of RC columns with different geometrical and mechanical properties are assessed.And finally, some concluding remarks are discussed in section 5.

CONSTITUTIVE MODELS
In this section the constitutive models adopted to represent the behavior of the materials involved in the proposed approach are described.For the concrete and interface steel-concrete, two constitutive models available in the literature based on the CDMT were implemented.For this class of models, the mechanical behavior of a damaged material is usually described by using the notion of the effective stress, together with the hypothesis of mechanical equivalence between the damage and Latin American Journal of Solids and Structures 12 (2015) 520-541 the undamaged material (Murakami 2012).In this paper, for the constitutive models based on the CDMT, the effective stress tensor, σ , will assume the following form: where C is the fourth order linear-elastic constitutive tensor, ε is the second order strain tensor and (:) is the tensorial product contracted in two indices.
The reinforcement bars are modeled using the classical elastoplastic model.

A continuum isotropic damage model for concrete
To describe the nonlinear behavior of the concrete, the rate-independent version of the constitutive model based on the CDMT presented by Cervera et al. (1996) was implemented.An important feature of this model is its capability to describe different responses when subject to tensile and compression.To do this, two different scalar damage variables, d + and d − , are defined for tensile and compression, respectively, which take values from zero (undamaged material) to one (fully damaged material).In addition, the effective stress tensor (equation 1) is split into tensile ( ) − σ and compressive ( ) With the aid of the above definitions, the nominal stress tensor for this model is obtained by reducing each part of the effective stress tensor, according to its respective damage variable in tension and compression, ( ) ( ) The equivalent effective tensile and compression norms are defined by the equations 4 and 5, respectively.These positive scalar variables are used to define concepts as loading, unloading and reloading for general stress states, The parameter is a material property that depends on the relation between the uniaxial and biaxial compressive strength of the concrete, β .According Cervera et al. (1996), typical values for concrete are: To describe the evolution of the damage variables d + and d − , two exponential functions to represent the material degradation process in tension and compression are defined.These damage variables are able to determine the softening behavior in tension and the hardening/softening in compression, after reaching the initial elastic limit in tension and compression, respectively, In an uniaxial tensile test, the expression shown in equation 8 is able to represent the softening in the stress-strain curve as a process degradation of quasi-brittle materials.To satisfy the mesh objectivity condition, the energy dissipated by the material in tension must to be properly related to the fracture energy of the material.Therefore, the softening parameter A + is derived from the ratio between the material fracture energy and the geometric factor, ch l , termed characteristic length, which corresponds to the width zone where the degradation concentrates, such that: is written in terms of the tensile strength, t f , the elastic modulus, E , and the (tensile) fracture energy of the material, f G + .The characteristic length, depends on the spatial discretization, and in this paper, is assumed to be the square root of the finite element area.
The parameters A − and B − are defined so that the stress-strain curve of the concrete satisfies two previously selected points of an uniaxial experimental test.
The algorithm of the damage model described in this section was implemented as shown in Table 1.More details about the formulation can be found in Cervera et al. (1996).

Numerical example for validation
To illustrate the behavior of the constitutive model adopted for the concrete, uniaxial load simulations for tensile and compressive states were carried out by imposing horizontal displacement in-crements, in a total of 100 for each branch (tension and compression).The Figure 1 shows the test setup of a fictitious concrete block numerically analyzed with a square cross section ( 2 100mm ) discretized with two three-node triangular finite elements in plane stress condition, with an out-ofplane thickness of 100mm .Table 1: Algorithm of the damage model proposed by Cervera et al. (1996).The parameters adopted in the analyses are: Young's modulus GPa 23.6 Figure 1 shows the normalized stress versus strain curves obtained in the numerical analyses.As can be noted, for both compression and tension tests, the damage model employed is able to describe the main characteristics of the mechanical behavior of the concrete under static load.

Elastoplastic model for steel reinforcements
To simulate the behavior of the steel bars, an elastic-perfectly plastic model was employed.The Figure 2 illustrates the stress-strain relation for this model characterized by the Young's modulus s E and the yield stress y σ .This model is widely described in the literature and a detailed description can be found in Simó and Hughes (1998).

J2 continuum damage model for steel-concrete interface
To represent the bond-slip behavior between the steel bars and concrete, a constitutive damage model consistent with the Continuous Strong Discontinuity Approach (CSDA) is adopted, as proposed by Manzoli et al. (2012).This model is formulated on the effective stress field associated with a scalar damage variable, [0,1] J d ∈ , which degrades only the deviatoric part of the effective stress.Thus, the constitutive relation of the J2 damage model can be written as where S and V σ are the deviatoric and the volumetric components of the effective stress tensor (equation 1), such that: The damage criterion can be expressed in terms of the effective stress as: Latin American Journal of Solids and Structures 12 (2015) 520-541 where J r is the current damage threshold, and J τ is the equivalent stress, given by: The evolution of the damage threshold can be expressed in a closed form, always using the highest value reached by J τ during the loading process, i.e., 0 max( , ) In this paper the evolution of the damage variable assumes the following form: where ( ) J q r is a function (hardening/softening law) that can be calibrated from experimental tests or analytical expressions available in standard codes to describe the steel-concrete interface.
Table 2 illustrates the algorithm of the J2 damage model implemented.More details about this constitutive model can be found in Manzoli et al. (2012).
Calculate and :

Verify if true :
false : OUTPUT: , Table 2: Algorithm of the J2 damage model.

INTERFACE FINITE ELEMENT
In the methodology applied in this work for modeling reinforced concrete columns, the interface steel-concrete is discretized using interface finite elements (IFEs), as proposed by Manzoli et al. (2012).The IFE is the standard triangular (three-node) finite element with very high aspect ratio.
As described by the authors, when the aspect ratio increases (ratio of the largest to the smallest dimension), the element's strains also increase approaching the kinematics of the strong discontinuity, as in the case of the CSDA (Simó et al., 1993;Oliver et al., 1999;Oliver, 2000).Therefore, based on the same principles as those of CSDA, it can be stated that bounded stresses can be obtained from unbounded strains by means of a continuum constitutive relation, which tends toward a discrete constitutive relation as the aspect ratio increases.
Latin American Journal of Solids and Structures 12 (2015) 520-541 The main advantages of the use of IFEs can be summarized as follows: • it is not necessary a fine mesh in the interface steel-concrete, since the aspect ratio of this element can be very high; • discrete constitutive relations are no longer necessary and contact elements (zero-thickness elements) or specific variational formulations are not required to introduce the effects of the cohesive stresses between elements; • a continuum constitutive relation can be used to described the bond-slip mechanism.Regarding the main features of the use of IFEs, it is also important to remark that: • the matrix remains continuous and the displacement compatibility between elements is maintained; • the IFEs overlap the concrete matrix; • when the IFE height is much smaller than the length of the element base, only slip between the steel and concrete is allowed.In this paper, the J2 continuum damage model described in section 2.3 is applied on the IFE to describe the bond degradation.The function ( ) q r of equation ( 15) can be adjusted to reproduce any given relation law between the bond shear stress ( ) s τ and the interface slip s .This paper uses the interface model proposed by CEB-FIP model code 90 (1993) depicted in Figure 3.This model is expressed in terms of the parameters: max τ (maximum shear strength), f τ (residual shear strength), 1 s , 2 s , 3 s and α , which depend on the concrete strength, the bar geometry (ribbed or smooth), the confining situation (confined or unconfined) and bond condition (good or poor).Thus, the hardening/softening law for this model becomes: and 0 0 q = .The parameters G and h are the shear modulus and the height of the IFE, respectively.The features described above aimed to give an idea about the use of the IFE.formulation can be found in Manzoli et al. (2012)

Pullout test
In order to investigate the capability of the J2 damage model (described in the the interface finite elements (presented in section 3 test with the same parameters that will be used fo carried out.
Figure 4 shows the test setup, including the boundary conditions and finite element mesh.steel bar of diameter 10mm is embedded in a concrete block with a rectangular cross section of 75mm 100mm × .The concrete specimen is discretized with ments in plane stress condition, with an out of plane thickness of tized by 6 two-node (one-dimensional) truss elements with cross 2 78.54mm .The bar elements are connected to the concrete elements by ments with thickness of 1mm h = and out-of-plane thickness of the perimeter of the cross-section of the bar in contact with the concrete.has been carried out by imposing a total displacement of 2mm (see Figure 4).
The concrete is assumed to be linearly elastic, with a Young's modulus of Poisson's ratio of Figure 5 shows the average bond stress (corresponding area of the contact surface between the bar and concrete ent stages of the bond degradation process can be observed The features described above aimed to give an idea about the use of the IFE.Details about its Manzoli et al. (2012).
In order to investigate the capability of the J2 damage model (described in the subsection 2.3) with presented in section 3) in representing the interface behavior, a pullout test with the same parameters that will be used for modeling the reinforced concrete columns is , including the boundary conditions and finite element mesh.A is embedded in a concrete block with a rectangular cross section The concrete specimen is discretized with 96 three-node triangular finite elements in plane stress condition, with an out of plane thickness of 75mm .The steel bar is discredimensional) truss elements with cross-section area of s A = To investigate the complex behavior of reinforced concrete columns subjected to eccentric compre sive load, a series of experimental tests was carried out by variables like concrete strength, the slenderness ratio and the longitudinal steel ratio on the ult mate load capacity were investigated by the researchers.A total of 30 tied reinforced columns with a square cross heights, and L three slenderness ratio were defined: λ = crete compressive strengths, 25.5MPa (normal or conventional strength concrete (middle strength concrete -MSC) and 86.2MPa dinal steel ratios of 1.98% ρ = and 3.95% ρ = considered a load eccentricity of 24mm and the rate of loading was controlled by a constant incr ment rate of the vertical displacement.

Finite element model
In this paper, some of the RC columns experimentally tested by cally investigated using the ingredients of the developed methodology.

Description of the columns experimentally tested by Kim and Yang (1995)
To investigate the complex behavior of reinforced concrete columns subjected to eccentric compre sive load, a series of experimental tests was carried out by Kim and Yang (1995).The influence variables like concrete strength, the slenderness ratio and the longitudinal steel ratio on the ult mate load capacity were investigated by the researchers.A total of 30 tied reinforced columns with a square cross-section of .Moreover, three different co (normal or conventional strength concrete -NSC), 63.5MPa 86.2MPa (high strength concrete -HSC), and two longit 3.95% were used.Figure 6 illustrates the test setup.It was and the rate of loading was controlled by a constant incr columns experimentally tested by Kim and Yang (1995) are numer cally investigated using the ingredients of the developed methodology.Table 3 shows the geome according to the parameters of the interface model suggested by .The nonlinear behavior observed in the first stage of the curve corinduces a micro expression .The stage II .Finally, the horizontal To investigate the complex behavior of reinforced concrete columns subjected to eccentric compreshe influences of variables like concrete strength, the slenderness ratio and the longitudinal steel ratio on the ultiand three different were experimentally analyzed.Consequently, .Moreover, three different con-

63.5MPa
HSC), and two longituillustrates the test setup.It was and the rate of loading was controlled by a constant increare numerishows the geomet-rical and mechanical properties of the nine columns analyzed with a reinforcement ratio of and different values of the concrete strength and slenderness ratio.For all the numerical models, the same boundary conditions have been applied.Figure 7 illustrates the boundary conditions for the columns with 10 λ = .An incremental vertical compression load was applied until failure.The nonlinear solution is obtained using displacement control method, in which the load process is controlled by the lateral deflection of a left-side middle point of the columns (see Figure 7).
Table 4 shows the total number of nodes and the number of each type of finite elements of the numerical analyses with slenderness ratio of  Due to the rigid behavior of the two parallel steel plates, a linear elastic material with a high modulus of elasticity is employed.For the steel bars, the elastic in subsection 2.2 is assumed with Young' modulus of 387MPa y σ = .The interface behavior is described by the J2 damage model described in tion 2.3 with the same parameters employed in the pullout analysis performed in the For the three types of concrete, NSC, MSC and HSC, the damage model proposed by al. (1996), described in subsection 2.1 is used.The parameters employed are listed in Detailed data of the finite element mesh built for the numerical models.
Due to the rigid behavior of the two parallel steel plates, a linear elastic material with a high modulus of elasticity is employed.For the steel bars, the elastic-perfectly plastic material described is assumed with Young' modulus of 210GPa s E = and yield strength of .The interface behavior is described by the J2 damage model described in subse with the same parameters employed in the pullout analysis performed in the subsection 3.1 For the three types of concrete, NSC, MSC and HSC, the damage model proposed by Cervera et is used.The parameters employed are listed in Table 5.The same parameters utilized to obtain the stress-strain curve shown in Figure 1 were used for the analyses of columns with NSC (C1, C4 and C7).For the columns with MSC (C2, C5 and C8) and HSC (C3, C6 and C9), the parameters as the Young's modulus and tensile strength have been calculated by the expressions suggested by Kim and Yang (1995), 3, 320 6900 respectively.Moreover, for these columns (MSC and HSC) the compression stress thresholds were obtained by imposing that the stress-strain curve for an uniaxial compressive test (as illustrated in Figure 1) reaches the respective compressive strength for the MSC and HSC.Table 5 presents all the mechanical properties adopted for the damage model proposed by Cervera et al. (1996).Table 5: Mechanical properties adopted for the damage model proposed by Cervera et al. (1996).

Comparison between numerical and experimental results
In this section the results obtained in the numerical simulations using the proposed methodology are compared with the experimental results obtained by Kim and Yang (1995).In addition, some results are also compared with the numerical analyses performed by Majewski et al. (2008).

NSC-Columns
This section presents the results obtained in the numerical simulations of the NSC-Columns with three different slenderness ratios: (column C7), as given in Table 3.
First, the numerical results have been compared with those results obtained by Majewski et al. (2008) in terms of axial force versus axial displacement (Figure 8).Then, they also have been compared against the experimental data obtained by Kim and Yang (1995) in terms of axial force versus lateral deflection, as shown in Figure 9.It can be noted that the results obtained with the proposed methodology are in very agreement with the results available in the literature.These graphics show that the ultimate load decreases with increasing slenderness ratio.
Figures 10, 12 and 14 show the tensile damage distributions along the columns C1, C4 and C7, respectively.In turn, the compressive damage distributions are depicted in Figures 11, 13 and 15, for these same columns.For the column C1, with the smallest slenderness ratio, the failure process occurs by the co pressive strain in concrete (Figure 11).Figures 8 and 9 show the small axial displacement and la eral deflection presented by this column.

Latin American Journal of Solids and Structures
On the other hand, for the columns C4 and C7 with the highest slenderness ratios, a large lateral and axial deflection are observed.For these cases, the failure process is driven by the increasing in tensile and compressive stress in a small region at the middle height of the columns, with a pr nounced plastic yielding of the longitudinal steel bars under normal tensile stress, forming hinge mechanism or local buckling, as illustrated in Figures 1 structural curves shown in Figures 8 and 9.For the column C1, with the smallest slenderness ratio, the failure process occurs by the com-

MSC-and HSC-Columns
). Figures 8 and 9 show the small axial displacement and lat-On the other hand, for the columns C4 and C7 with the highest slenderness ratios, a large lateral For these cases, the failure process is driven by the increasing in tensile and compressive stress in a small region at the middle height of the columns, with a pronounced plastic yielding of the longitudinal steel bars under normal tensile stress, forming a plastic hinge mechanism or local buckling, as illustrated in Figures 12(b   good agreement with the experimental results obtained by Kim and Yang (1995).Once again, the responses show that the ultimate load decreases with increasing slenderness ratio.The tensile and compressive damage distributions for these columns can be seen from Figure 18 to 23.In general, the failure process of these columns is very similar of those obtained for the NSC-Columns, except by its magnitude.Due to the small lateral deflection presented by the columns C2 and C3 (see Figures 16 and 17), with a small slenderness ratio, the process failure of these columns is conducted by crushing of the concrete in the compressed region of the columns, as illustrated in Figure 21 for the HSC-Column.
The columns C5 and C6, C8 and C9 with slenderness ratio of , respectively, present a large lateral deflections (Figures 16 and 17).For this reason, these columns present a degradation process driven by the increase of tensile stress providing the formation, growth and coalescence of cracks, with a pronounced yielding of the tensile longitudinal steel bars, which culminates with the strain localization in a small region at the middle height of the columns, exhibiting a plastic hinge mechanism (see Figures 18(b   To conclude, it is also important to note that the methodology employed was able to capture the influence of the concrete strength.The results showed that for columns with slenderness ratio 10 λ = , the load capacity increases, with the increasing of compressive strength.The same behavior does not occur to the columns with higher slenderness ( ), because the yielding of the longitudinal steel bar in tension is the factor that limits the ultimate load capacity.The curves illustrated in Figures 9, 16 and 17, for slender columns, shows that the increasing of concrete strength does not provide a significant gain in terms of the structural strength.
It should also be seen in Figures 18, 20 and 22, that the value of tensile damage is intensified in almost horizontal alignments of finite elements, from the tensioned face, with almost constant spacing.These lines of elements represent the tensile cracks, and the space between the cracks is due to the bond degradation between steel bar and concrete.

CONCLUSIONS
In the present work the authors sought to contribute for a better understanding of the failure behavior of reinforced concrete columns subject to eccentric compressive loads.With that purpose, a numerical approach has been proposed and its features have been implemented into an in-house finite element code.To represent the concrete behavior, the rate-independent version of the damage constitutive model proposed by Cervera et al. (1996) with different damage variables for tension and compression was implemented.Firstly, this damage constitutive model was assessed through the numerical analysis of a simple concrete block under tension and then under compression.The results showed that the damage constitutive model is able to represent the structural nonlinearity resulting from the fracture process of concrete in the tensioned region and the crushing failure of concrete in the compressed region.The complex bond-slip phenomenon on the steel-concrete interface was represented by the use of an interface finite element with very high aspect ratio and a J2 damage constitutive model, which degrades only the deviatoric part of the effective stress tensor, as described by Manzoli et al. (2012).A pullout test was performed initially to calibrate the parameters of the interface model for posterior analysis of the columns.To represent the steel bars, twonode linear finite elements (truss elements) and an elastic-perfectly plastic model, which were available in the finite element code, were used.
In the sequence, using the ingredients described above, the failure behavior of some of the columns experimentally tested by Kim and Yang (1995) were investigated.In addition, some results were also compared with those obtained numerically by Majewski et al. (2008).The variables considered in the analysis were the slenderness ratio ( ) and the concrete strength ( 25.5MPa , 63.5MPa and 86.2MPa ) for a steel ratio of 1.98% in all the cases.
The numerical results exhibited good agreements with both, the experimental and the numerical results available in the literature, which can be confirmed by analyzing the curves obtained in terms of axial force versus axial displacement (Figure 8) and axial force versus lateral deflection (Figures 9,16 and 17).In general, the methodology was able to reproduce the influence of the slenderness ratio and the concrete strength in the same way that was reported by Kim and Yang (1995) and by Majewski et al. (2008).The responses show that the ultimate load decreases with increasing slenderness ratio.Similar influence of the concrete strength on the failure process of NSC, MSC and HSC columns was observed.Moreover, the methodology was also able to reproduce the main fea-Latin American Journal of Solids and Structures 12 (2015) 520-541 tures of the failure process of the columns.Two distinct failure modes were observed.For less slender columns, the concrete plays an important role in the compression region.For them, the ultimate load occurs due to crushing of concrete.For the slender columns, the degradation process is driven by the increase of tensile stress providing the formation, growth and coalescence of cracks, with a pronounced yielding of the tensile longitudinal steel bars, which culminates with the strain localization in a small region at the middle height of the columns, exhibiting a plastic hinge mechanism.
Therefore, it is the authors' conclusion that the presented methodology is able to represent efficiently, qualitatively and quantitatively, the failure mechanism of reinforced concrete columns under eccentric compression load, reproducing efficiently the structural degradation of slender reinforced concrete columns.The main advantage obtained by the use of the methodology proposed was its easy application and calibration of the parameters of the assumed material models, including the complex interface behavior between the steel and concrete.In future investigations, the authors intend to extend the analyses for 3D problems and to consider the influence of the steel ratio in the failure process.Also, a more deep analysis about the crack pattern should be made.
oct σ − and oct τ − are the octahedral normal and shear stresses, respectively, obtained from − σ .Two independent damage criteria, one for tension and other for compression, are expressed as: r + and r − are the current damage thresholds, which are updated continuously to control the size of the expanding damage surface.The boundary damage surfaces for the effective stresses are expressed threshold for damage.The evolution of the damage thresholds can be expressed in a closed form, always using the highest values reached by τ + and τ − , during the loading process, i.e.,

Figure 1 :
Figure 1: Normalized stress-strain curves for uniaxial tension and compression loads.

Figure 2 :
Figure 2: Constitutive model adopted for the steel bars.

Figure 4 :
Figure 4: Pullout test: (a) detail of the interface elements; (b) geometry, finite element mesh and boundary cond tions.

Failure
behavior modeling of slender reinforced concrete columns subjected to eccentric load 528 Latin American Journal of Solids and Structures 12 (2015) 520-541

Failure
behavior modeling of slender reinforced concrete columns subjected to eccentric load Latin American Journal of Solids and Structures 12 (2015) 520-541 according to the parameters of the interface model suggested by .The nonlinear behavior observed in the first stage of the curve co responds to the rebar ribs penetration into the mortar matrix phenomenon, which induces a micro cracking formation and a local crushing.This stage is represented by the expression FIP model code 90 (1993), for the interval 1 0 s s ≤ ≤ .The stage II illustrates the linear softening behavior of the bond stress for 2 3 s s s ≤ ≤ .Finally, the horizontal (stage III) represents the residual bond capacity, f τ .slipcurve obtained with the J2 damage model.

Figure 6 :
Figure 6: Geometry and boundary condition of reinforced concrete columns experimentally tested byKim and Yang (1995).
node triangular Latin American Journal of Solids and Structures finite elements were used for the concrete and for the steel plates.The steel bars are represented by two-node linear finite elements (one-dimensional truss elements), which are connected to the nodes of the concrete elements by interface finite elements, as depicted in

Figure 7 :
Figure 7: Numerical model built for the analysis of the

Failure
behavior modeling of slender reinforced concrete columns subjected to eccentric load Latin American Journal of Solidsand Structures 12 (2015) 520-541were used for the concrete and for the steel plates.The steel bars are represented by dimensional truss elements), which are connected to the nodes s by interface finite elements, as depicted in Figure7.Numerical model built for the analysis of the columns with 10 were used for the concrete and for the steel plates.The steel bars are represented by dimensional truss elements), which are connected to the nodes Due to the rigid behavior of the two parallel steel plates, a linear elastic material with a

Figure 11 :
Figure 11: Compressive damage distribution for the NSC and (b) at the ultimate load.

Figure 12 :
Figure 12: Tensile damage distribution for the NSC-Column with at the ultimate load.

Figures
Figures 16 and 17 show the numerical results obtained for MSC spectively, for the three different slenderness ratios considered.For both cases the results are in Latin American Journal of Solids and Structures 12 (2015) 520-541 Figures 16 and 17 show the numerical results obtained for MSC-Columns and HSC-Columns, respectively, for the three different slenderness ratios considered.For both cases the results are in

Figure 13 :Figure 14 :
Figure 13: Compressive damage distribution for the NSC

Figure 15 :
Figure 15: Compressive damage distribution for the NSC

Figure 18 :
Figure 18: Tensile damage distribution for the MSC (b) at the ultimate load.

Figure 22 :Figure 23 :
Figure 22: Tensile damage distribution for the HSC-Column with (b) at the ultimate

Table 3 :
Geometrical and mechanical properties of the columns analyzed.