# 3D Fiber-Based Frame Element with Multiaxial Stress Interaction for RC Structures.

1. IntroductionThe increasing interest towards performance-based engineering requires robust and accurate numerical models for strength and ductility assessment of reinforced concrete (RC) structures. A great majority of RC structures, such as bridges and frame buildings, can be idealized with one-dimensional frame elements for global nonlinear analysis. However, these elements usually rely on the classical Euler-Bernoulli beam assumption, that is, plane sections remain plane and orthogonal to the element axis, and uniaxial stress conditions, which is only valid in the so-called B-regions. As an example, Figure 1 shows the results from such an element applied to the modeling of a floor subassembly subjected to a vertical concentrated load in the middle of the floor beam [1]. This creates multidirectional shear flow in the spandrel beam coupled with normal stresses arising from bending. The frame element is unable to reproduce the cracked torsional stiffness of the spandrel beam, thus significantly overestimating the overall cracked stiffness and strength of the test. In the same figure, the results from 3D nonlinear finite element analysis (FEA) have been included, demonstrating much better agreement [2].

From the above, it is clear that standard flexural frame elements, commonly used in commercial software, present limitations under multiaxial loading conditions, which have motivated the development of more advanced frame elements for global nonlinear static and dynamic analysis. Several formulations have been proposed in 2D addressing the problem of axial-flexure-shear interaction. Some authors used an equivalent truss model in shear combined with the Euler-Bernoulli beam assumption and uniaxial stress response in flexure [3]. Further improvements were obtained when shear deformations were included in the element kinematics (Timoshenko beam theory) [4], which also allowed using two-dimensional constitutive models at each individual fiber based on either fixed [5-7] or variable shear strain profile [8]. Different constitutive models have been used, ranging from orthotropic smeared crack models (e.g., modified compression field theory (MCFT) [9] and fixed angle softened truss models (FA-STM) [10]) to softened plasticity damage models [11, 12], and microplane models [13], among others. A detailed review can be found in [14].

Research devoted to RC 3D frame elements has been more limited mostly due to complexities in the development of robust and accurate triaxial material models for concrete and their implementation into fiber frame elements. Mazars et al. [15] developed a 3D displacement-based Timoshenko beam element with a 3D constitutive model for concrete based on continuum damage mechanics. The Saint-Venant torsion was introduced by means of the warping function derived from linear elastic analysis under pure torsion. Gregori et al. [16] used a similar element to that used by Mazars et al. but with a parabolic shear strain profile and Coulomb theory for torsion. The constitutive model was based on the MCFT. The cross section was subdivided into three regions based on the dominant stress conditions: 1D, 2D, and 3D, where the corresponding constitutive models were applied. Mullapudi and Ayoub [17] formulated a 3D flexibility-based fiber element using the FA-STM. The beam kinematics was based on the Timoshenko beam theory in conjunction with the Coulomb torsion, that is, no warping was considered.

The abovementioned work was based on the classical Timoshenko beam theory with either uniform or parabolic shear profiles. This approach does not satisfy interfiber sectional equilibrium; hence, more advanced (higher-order) beam theories have been recently investigated which satisfy interfiber sectional equilibrium in an average sense [18-21]. However, only a limited number have been successfully applied to 3D analysis of RC frame structures [22].

In the present work, a nonlinear frame element based on the Timoshenko beam assumption is developed. The cross section is discretized into several fibers, which allows direct consideration of multiaxial stress interaction at the material level. A triaxial smeared crack constitutive model is implemented for concrete, whereas reinforcement is treated as smeared within the concrete with uniaxial stress-strain behavior. The Saint-Venant theory is used for modeling torsion, which is suitable for noncircular solid cross sections with warping effects. The element formulation has been implemented into a force-based frame element and applied to monotonic and cyclic analysis of RC columns and beams. Comparison against experimental results shows similar levels of accuracy to those obtained from detailed 3D FEA with solid elements.

2. Frame Element Formulation

In the Timoshenko beam theory (TBT), cross sections are not constrained to remain orthogonal to the beam axis, as in the Euler-Bernoulli beam theory, due to shear distortions of an infinitesimal beam element. The kinematics of any point of the beam is completely described by three independent displacements and rotations of the beam axis as follows (Figure 2):

d = ([u.sub.o] [v.sub.o] [w.sub.o] [[theta].sub.x] [[theta].sub.y] [[theta].sub.z]. (1)

The total displacement at any given point of the section is related to the beam axis displacements as

u (x, y, z) = [u.sub.o] (x) - [[theta].sub.z](x)y + [[theta].sub.y](x)z + [omega] (y, z)[alpha] (x), (2)

v(x,y,z) = [v.sub.o] (x) - [[theta].sub.x] (x)z, (3)

w(x, y, z) = [w.sub.o] (x) + [[theta].sub.x] (x)y, (4)

where u, v, and w are the axial and transverse displacements of a generic point of the section with coordinates (x, y, z), [omega] is the warping function, and a is the warping parameter equal to [partial derivative][[theta].sub.x]/[partial derivative]x according to the Saint-Venant theory. Assuming small displacements, axial and shear strains at a given section fiber are given as

[mathematical expression not reproducible], (5)

where [[epsilon].sub.o] is the average axial strain, [X.sub.y] and [X.sub.z] the bending curvatures, [[gamma].sub.oy] and [[gamma].sub.oz] the section shear deformations, and [alpha] is the derivative of the torsional rotation [[theta].sub.x]. The above relationships can be expressed in matrix as [epsilon] = [B.sub.[epsilon]]e, where [B.sub.[epsilon]] is the strain-displacement compatibility matrix, [epsilon] = ([[epsilon].sub.x], [[gamma].sub.xy], [[gamma].sub.xz]), and e = ([[epsilon].sub.o], [[gamma].sub.oy], [[gamma].sub.oz]).

The section forces consistent with the section strains are obtained from integration of the axial and shear stress distributions over the cross section as follows:

[mathematical expression not reproducible]. (6)

The relation between section forces and section strains is given by the section stiffness matrix. In the nonlinear case, incremental section forces and strains are related by the tangent stiffness matrix of the section as follows:

[mathematical expression not reproducible], (7)

where [bar.D] is the condensed (3 x 3) tangent material stiffness matrix. Note that [B.sub.[sigma]] and [B.sub.[epsilon]] are different; hence, [K.sub.s] will not be symmetric in the general case, with the exception of the circular cross section for which the warping function is zero. The [bar.D] matrix is obtained from condensation of the (6 x 6) material stiffness matrix assuming that the total stresses [[sigma].sub.y], [[sigma].sub.z], and [[tau].sub.yz] are zero in the beam element as follows:

[mathematical expression not reproducible], (8)

where [[sigma].sub.c,i] and [[sigma].sub.s,i] are the stresses in concrete and steel, respectively, and [[rho].sub.si] (i = y, z) is the reinforcement ratio at each individual section fiber. This process is performed iteratively at each section fiber until convergence, which requires implementing an additional loop within the fiber state-determination procedure.

The warping function in (2) is obtained from the solution of an elastic prismatic beam under pure torsion based on the Saint-Venant theory of torsion, which is defined by the following boundary value problem:

[[partial derivative].sup.2][omega]/[partial derivative][y.sup.2] + [d.sup.2][omega]/[partial derivative][z.sup.2] = 0 in [OMEGA], (9)

[partial derivative][omega]/[partial derivative]y [n.sub.y] [partial derivative][omega]/[partial derivative]z = z[n.sub.y] - y[n.sub.z] on [partial derivative][OMEGA], (10)

where [OMEGA] and [partial derivative][OMEGA] represent the cross-sectional domain and its boundary, respectively, and [n.sub.y] and [n.sub.z] are the components of the vector normal to the boundary. Equations (9) and (10) represent a boundary value problem that can be solved using the finite element technique. In the present case, four-node plane isoparametric elements were used to discretize the cross section and calculate the elastic warping function w. Figure 3 shows the calculated warping functions for the rectangular and hollow cross sections used for verification in Section 5.

3. Constitutive Model

A 3D orthotropic smeared crack material model is formulated for cracked reinforced concrete, which represents an extension of the 2D model presented in [23] for membrane elements. A similar approach was also adopted in [17, 24]. The model is formulated in the material directions (Figure 4), which are coincident with the crack directions determined at first cracking as

min([[epsilon].sub.1] + [[epsilon].sub.cp1], [e.sub.2] + [[epsilon].sub.cp2], [[epsilon].sub.3] + [[epsilon].sub.cp3]) (11)

where [[epsilon].sub.cpi] is the plastic compressive strain (positive) and [[epsilon].sub.i] is the principal strain (negative for tension). [[epsilon].sub.i] can be obtained upon solving the eigenvalue problem for the strain tensor in the xyz reference system, [[epsilon].sub.xyz]. After the first cracking, the material directions are kept constant throughout the analysis, assuming that secondary cracks are orthogonal to the initial ones. Total strains in the material directions, [[epsilon].sub.123], are obtained from [[epsilon].sub.xyz] as

[[epsilon].sub.123] = T[[epsilon].sub.xyz], (12)

where T is the standard (6 x 6) transformation matrix. The corresponding concrete stresses, [[sigma].sub.123], are determined from equivalent uniaxial constitutive laws for axial and shear stresses (Figure 5(a)). In compression, the monotonic envelope is defined by the Popovics curve scaled by the lateral confining factor and a softening coefficient due to orthogonal tensile strains [9] given as

[[beta].sub.[sigma]] = 1 / 0.80 - 0.34 ([[epsilon].sub.av] / [[epsilon].sub.o]), [[epsilon].sub.av] = [[epsilon].sub.j] + [[epsilon].sub.k] / 2, (13)

where [[epsilon].sub.j] and [[epsilon].sub.k] are the total tensile strains in the orthogonal directions j and k.

In tension, a physically based approach is used such that average concrete stresses are determined using equilibrium, compatibility, and constitutive relations applied on a cracked concrete element with three different regions as follows: an elastic, fully bonded region (b), a crack region (cr), and a bond-slip region (sl) where the slip between concrete and reinforcement occurs (Figure 5(b)). Thus, average concrete stress in the ith material direction can be obtained as

[[sigma].sup.av.sub.ci] = (1 - [[lambda].sub.i])[s.sub.i] [[sigma].sup.b.sub.ci] + ([[lambda].sub.i] - [[epsilon].sub.i])[s.sub.i] [[sigma].sup.sl.sub.ci] + [s.sub.i][[epsilon].sub.i] [[sigma].sup.cr.sub.ci] (14)

where [[sigma].sup.b.sub.ci], [[sigma].sup.sl.sub.ci], and [[sigma].sup.cr.sub.ci] are concrete stresses in the (b), (sl), and (cr) regions, respectively. [s.sub.i], [[lambda].sub.i], and [[epsilon].sub.i] are the average spacing, unbonded length parameter, and tensile strain in the ith direction, respectively. This approach allows full characterization of the tension stiffening, crack-closing, and crack-opening regimes based on principles of mechanics. Further details can be found elsewhere [23, 25].

Shear stresses, [[tau].sub.ij], arising from deviation between principal directions and crack directions are related to the shear strains, [[gamma].sub.ij], as

[[tau].sub.ij] = [G.sub.ij][[gamma].sub.ij] < [v.sub.ci,max], [G.sub.ij] = 1/2 [[sigma].sub.i] - [[sigma].sub.j] / [[epsilon].sub.i] - [[epsilon].sub.j], (15)

where [[sigma].sub.i], [[sigma].sub.j] and [[epsilon].sub.i], [[epsilon].sub.j] are the normal stresses and strains in the ith and jth crack directions, respectively. Note that, in this approach, dowel action is not directly accounted for. Shear failure due to aggregate interlock is included by setting an upper limit on the maximum shear stress transferred along cracks [9] as follows:

[v.sub.ci,max] = 0.18 [square root of ([f'.sub.c])] / 0.31 + (24[w.sub.ij] / 16 + [a.sub.g]), (16)

where [a.sub.g] is the maximum aggregate size and [w.sub.ij] is the average crack width (mm). Assuming smeared cracks, the average crack width is given as

[w.sub.ij] = [s.sub.ij][[epsilon].sub.n], (17)

where [[epsilon].sub.n] is the maximum (in absolute value) tensile strain normal to the crack. The average crack spacing, [s.sub.ij], is taken as the average between i and j directions as follows:

[s.sub.ij] = [s.sub.i] + [s.sub.j] / 2, (18)

where the average crack spacing in the ith direction is determined from those in the x, y, and z directions as follows:

[s.sub.i] = 1 / ([n.sub.1i]/[s.sub.x]) + ([n.sub.2i]/[s.sub.y]) + ([n.sub.3i]/[s.sub.z]), (19)

where [n.sub.ji] are the direction cosines between xyz and the crack directions.

Finally, steel stresses are determined based on the uniaxial Menegotto-Pinto hysteretic model [26] with modified isotropic hardening [27]. For the ith reinforcement group, which orientation is defined by the unitary vector [n.sub.si], the total axial strain can be obtained upon transforming [[epsilon].sub.xyz] to the corresponding direction. Then, steel stresses are obtained using the abovementioned model and transformed back to the xyz reference system as

[[sigma].sub.si,xyz] = [T.sup.T.sub.si] [([[rho].sub.si] [[sigma].sub.si] 0 0 0 0 0).sup.T], (20)

where [[sigma].sub.si] is the uniaxial steel stress, [[rho].sub.si] is the reinforcement ratio for the ith reinforcement group, and [T.sub.si] is the corresponding transformation matrix between reference systems. [[sigma].sub.si,xyz] stresses are added to the concrete stresses, [[sigma].sub.c,xyz], in order to obtain the total stresses, [[sigma].sub.xyz], which are used for the calculation of the vector of internal element forces. Similarly, the ith reinforcement contribution to the material stiffness matrix, [D.sub.si,xyz], is added to that provided by concrete, [D.sub.c,xyz], such that the total material stiffness matrix is given as

[D.sub.xyz] = [T.sup.T] [D.sup.c,123] T + [T.sup.T.sub.si] [D.sub.si] [T.sub.si], (21)

where

[mathematical expression not reproducible], (22)

where [E.sub.ci] and [E.sub.si] are the tangent moduli for concrete and steel and [G.sub.ij] is the shear stiffness.

4. Frame Element Implementation

The proposed frame model was implemented using Matlab[C] into a three-dimensional, force-based fiber element based on the state-determination formulation described by Spacone et al. [28]. In this formulation, once the global displacements are known, an iterative procedure is invoked at the element level in order to obtain the corresponding section strains at each integration point. This is performed in the basic element configuration upon removal of the rigid body modes as follows (Figure 6):

[DELTA]q = R[DELTA][d.sub.E], (23)

where [DELTA]q is the increment of element displacements in the basic configuration, [DELTA][d.sub.E] is the increment of global displacements, at a given load step, and R is the rigid body transformation matrix given as

[mathematical expression not reproducible]. (24)

The element forces, [DELTA]Q, can be obtained from the element flexibility matrix, [F.sub.e], which is known from the previous converged load step as

[DELTA]Q = [(Fe).sup.-1] [DELTA]q. (25)

Equation (25) is only an approximation since [F.sub.e] is changing within each load step due to material nonlinearities. Section forces, [DELTA]S, are interpolated exactly from element forces. At the kth integration point, this is given as

[DELTA][S.sub.k] = [b.sub.k][DELTA]Q, (26)

where [b.sub.k] is the force interpolation matrix at the kth integration point with coordinate [x.sub.k], which is given as

[mathematical expression not reproducible]. (27)

An estimate of the section deformations at the kth integration point, [mathematical expression not reproducible][e.sub.k], can be made based on the section flexibility matrix, [F.sub.s,k], which is obtained from the inverse of the section stiffness matrix given in (7):

[DELTA][e.sub.k] = [F.sub.s,k] [DELTA][S.sub.k]. (28)

At this point, the total strains at each individual fiber can be determined, and the corresponding stresses can be obtained using the constitutive models presented in Section 3. Integration of the total stresses over the cross section yields the section forces, [DELTA][S.sub.int,k], and an updated section flexibility matrix. The unbalanced section forces, [DELTA][S.sub.u,k], given by the difference between [DELTA][S.sub.k] and [DELTA][S.sub.int,k], are used to update the section deformations as [DELTA][e.sub.k] = [DELTA][e.sub.k] + [F.sub.s,k][DELTA][S.sub.u,k], which are further used to correct the element forces by the following quantity:

[DELTA][Q.sub.u] = [F.sup.-1.sub.e] [NIP.summation over (k=1)] [w.sub.k] [b.sup.T.sub.k] [DELTA][e.sub.k], (29)

where [w.sub.k] are the integration weights at the kth integration point and [F.sub.e] is the updated element flexibility matrix given as

[F.sub.e] = [NIP.summation over (k=1)] [w.sub.k] [b.sup.T.sub.k] [F.sub.s,k] [b.sub.k]. (30)

The procedure is repeated with the new element forces until satisfying element convergence, that is, when [parallel][DELTA][S.sub.k] - [DELTA][S.sub.int,k][parallel] < tol at each kth integration point. Once element convergence is achieved, the element forces obtained in the basic configuration are transformed back to the global configuration and assembled within the vector of internal forces of the structure. Global equilibrium between the vector of applied external forces and internal forces is checked at this point. If this is not satisfied, the residuum between external and internal forces and the tangent stiffness matrix of the structure are used in a subsequent iteration to compute a new increment of nodal displacements. A flow chart summarizing the numerical procedure is shown in Figure 7.

5. Verification Examples

5.1. Biaxial Bending Test. Bousias et al. [29] tested flexure-dominated RC cantilever columns in biaxial cyclic bending. Different paths of imposed lateral displacement were investigated, ranging from uncoupled bidirectional to circular patterns. All columns were 250 mm x 250 mm x 1490 mm rectangular specimens fixed to a stiff base foundation. The shear span-to-depth ratio was 6. The columns were well reinforced in shear with a double hoop pattern of 8 mm diameter stirrups at 70 mm spacing. Longitudinal reinforcement consisted of 16 mm bars uniformly distributed around the perimeter. Material properties are summarized in Table 1. A constant axial load was applied at the top of each column. Two tests are reproduced here, S1 and S7, with the displacement paths shown in Figure 8.

The columns were modeled with one single FB element with 5 IPs. The cross section was discretized into 10 x 10 fibers, with the longitudinal and transverse reinforcement smeared as shown in Figure 9. The contribution of the inclined hoops was decomposed in the x and z local directions and added to the hoops aligned with the axes. In any case, the influence of transverse reinforcement in the response of the member is rather small since the behavior is dominated by flexure. The main effect of the transverse reinforcement is confining of concrete, which was accounted for with an estimated confinement factor of 1.1.

The hysteretic lateral force-displacement response is shown in Figure 10. The numerical results agree well with the experimental data for both S1 and S7 specimens. Wide hysteretic loops can be observed for specimen S7 due to the nature of the imposed displacement. This aspect is well captured by the analytical model. For specimen S1, analytical results slightly underestimate the strength of the member, which might be related to material modeling of steel, although good agreement is observed in terms of hysteretic pinching and residual displacements. Figure 11 also reports the measured horizontal reactions in each direction. It is noted that the horizontal reactions are not totally decoupled for the S1 column, as it happens with displacements, presumably due to biaxial yielding. It shall be mentioned that testing and measuring the response of biaxially loaded columns involves several difficulties, such as misalignment of the cross-sectional principal axes with those of the actuator, and inclusion of P-delta effects, among others.

5.2. Torsion and Bending. The series of RC beams tested by Onsongo [30] in combined torsion and flexure are analyzed herein. The test campaign consisted of two series of five beams each, being the series TBO over-reinforced and designed to fail in concrete crushing without yielding of neither longitudinal nor transverse reinforcement, whereas the series TBU was under-reinforced, with either longitudinal or transverse reinforcement yielding prior to failure. Within each series, the ratio of applied torque to applied moment was varied from 0 to 5.059, being the rest of the test variables essentially constant. The test setup consisted of adjustable steel level arms at the ends of the specimen, where an eccentric jack load was applied, resulting in constant bending moment and torque along the specimen. Beams TBO3 and TBU3 were selected for analysis, both having the same cross section. Material and loading properties are summarized in Table 2.

Half of the beam was modeled using one single element with 5 IPs (Figure 12). Rigid elements were connected to the end of the beam in order to reproduce the support conditions and the eccentric load application. The torsional and bending rotations around z-axis were restrained at node 1. At node 4, the beam was free to rotate, and only the vertical displacement was restrained. Additionally, out-of-plane degrees of freedom were restrained in order to create an isostatic system and avoid stiffness matrix singularities. The cross section was discretized into 64 fibers, with both longitudinal and transverse reinforcements smeared within the web and flange regions as shown in Figure 12. The warping function of the hollow section is already presented in Figure 3. It was calculated as the difference between the warping functions of two solid rectangular sections, one with the outer dimensions and the other one with the dimensions of the inner core.

Figures 13 and 14 compare analytical and experimental results in terms of applied bending moment-curvature and applied torque-twist for the two beams. In general, good agreement can be found with the experiment. The moment-curvature response is better captured for the over-reinforced beam (TBO3) since it essentially behaves elastically in bending, whereas some differences can be observed for the under-reinforced beam (TBU3).

Results for the hoop strains at middepth are shown in Figure 15 for both beams. These are somewhat overestimated in the postcracking range, especially for beam TBO3, due to overestimation of the twist. In both beams, yielding of hoops occurred in both y and z reinforcements at several cross-sectional fibers.

From inspection of local strains in the longitudinal reinforcement, it was found that this did not yield for beam TBO3 but did yield for beam TBU3, with the values of 1.7[per thousand] for TBO3 and 2.5[per thousand] for TBU3 at the bottom flange. Concrete crushing occurred for the beam TBO3, reaching principal compressive strains of about 2.5[per thousand] and corresponding stresses of 18 MPa, slightly less than [f'.sub.c] due to concrete softening at orthogonal tensile strains. For TBU3, principal compressive strains remained below 2[per thousand].

5.3. Bending, Shear, and Torsion. The test of the floor subassembly by Collins and Lampert [1] presented in Figure 1 is considered here for verification. This test setup produces torsion, bending, and shear in the spandrel beam, and bending and shear in the floor beam. The deflection of the floor beam depends on the postcracked torsional and bending stiffness of the spandrel beam. Both beams were reinforced with longitudinal and transverse reinforcements, with the material properties summarized in Table 3.

The T specimen was modeled with 4 elements, with 5 IPs each (Figure 16). A mesh of 10 x 10 fibers was used to discretize both cross sections. Longitudinal reinforcement was smeared in the top three layers and in the four bottom layers for the floor beam. For the spandrel beam, it was smeared in the top and bottom three layers. Transverse reinforcement was smeared over the entire section domain for both beams. The torsional rotation was fixed at nodes 1 and 3, being the rest of rotations free. At node 5, the beam was simply supported. A vertical concentrated load was applied at node 4 under displacement control with 2 mm displacement increments.

Figure 17 shows the results for the applied load-deflection at node 4 and for the load-twist response of the spandrel beam. Furthermore, the results from 3D FEA presented by Valipour and Foster [2] are shown. The authors used brick elements with a 3D nonlinear cementitious material model, which is a fracture-plastic constitutive law for cracked concrete, and 1D elements for the reinforcement with a linear-hardening stress-strain relationship. Good correlation can be observed with the experiments in both cases. At 80 mm vertical deflection, the twist predicted by the model is 0.010 rad/m, whereas in the experiment, FEA is about 0.014 rad/m. In the frame element analysis, yielding of longitudinal top and bottom reinforcement occurred in the floor beam and only of the transverse reinforcement in the spandrel beam. In the experiment, yielding of the bottom longitudinal reinforcement also occurred in the spandrel beam.

6. Conclusions

Results from a 3D fiber-based frame element with multiaxial stress interaction have been presented. The proposed element is intended for nonlinear analysis of 3D frame structures such as those in the field of earthquake engineering, when consideration of multiaxial effects due to the presence of disturbed regions is important. Established frame analysis theories only consider the interaction between bending and normal forces with some "adjustments" for flexure-shear interaction. The proposed frame element formulation considers full coupling between bending, shear, and torsion in a consistent manner based on the fiber-section discretization approach and 3D constitutive models. Despite the conceptual simplicity of the implemented constitutive models and beam kinematics, the element accuracy was comparable to that from refined finite element analysis, at least for the investigated case study involving shear, bending, and torsion coupling. Further verification studies, for instance, regarding cyclic shear and torsion, are expected in the future depending on the availability of experimental data. Potential enhancement of the constitutive models will be the object of future research as well. The possibility of replacing traditional flexural frame elements, based on the Euler-Bernoulli beam assumption, for global nonlinear analysis of framed structures should be further addressed. An important aspect, not considered in the present work, is the computational time, which obviously lies in between flexural frame elements and solid finite elements.

https://doi.org/10.1155/2018/8596970

Data Availability

The authors declare that the data that support the findings of this study are available from the corresponding author on request. Furthermore, the data related to the verification examples are available within the article and in the papers cited in the references.

Conflicts of Interest

The authors declare that there are not conflicts of interest regarding the publication of this paper.

Acknowledgments

The Ph.D. scholarship support provided by the University School for Advanced Studies IUSS Pavia, Pavia, is greatly acknowledged.

References

[1] M. P. Collins and P. Lampert, "Redistribution of moments at cracking. The key to simpler torsion design?," ACI Structural Journal, vol. 35, pp. 343-384, 1973.

[2] H. R. Valipour and S. J. Foster, "Nonlinear reinforced concrete frame element with torsion," Engineering Structures, vol. 32, no. 4, pp. 998-1002, 2010.

[3] J. Guedes, P. Pegon, and A. A. Pinto, Fibre/Timoshenko Beam Element in CASTEM 2000, Special Publication Nr. I.94.31 Applied Mechanics Unit, Safety Technology Institute, Commission of the European Communities, Joint Research Centre, Ispra Establishment, Ispra, Italy, 1994.

[4] A. Marini and E. Spacone, "Analysis of reinforced concrete elements including shear effects," ACI Structural Journal, vol. 103, no. 5, pp. 645-655, 2006.

[5] P. Ceresa, L. Petrini, R. Pinho, and R. Sousa, "A fibre flexure-shear model for seismic analysis of RC-framed structures," Earthquake Engineering and Structural Dynamics, vol. 38, no. 5, pp. 565-586, 2009.

[6] S. Guner and F. J. Vecchio, "Analysis of shear-critical reinforced concrete plane frame elements under cyclic loading," Journal of Structural Engineering, vol. 137, no. 8, pp. 834-843, 2011.

[7] D. Ferreira, J. M. Bairan, and A. Marl, "Shear strengthening of reinforced concrete beams by means of vertical prestressed reinforcement," Structure and Infrastructure Engineering: Maintenance, Management, vol. 12, no. 3, pp. 394-410, 2015.

[8] A. Kagermanov and P. Ceresa, "An exact shear strain approach for RC frame elements with axial-shear interaction," in Proceedings of the VII European Congress on Computational Methods in Applied Sciences and Engineering, ECCOMAS, pp. 5-10, Crete, Greece, June 2016.

[9] F. J. Vecchio and M. P. Collins, "The modified compression field theory for reinforced concrete elements subjected to shear," ACI Journal, vol. 83, no. 2, pp. 219-231, 1986.

[10] X. Pang and T. T. C. Hsu, "Fixed angle softened truss model for reinforced concrete," ACI Structural Journal, vol. 93, no. 2, pp. 196-208, 1996.

[11] F. C. Filippou and A. Saritas, "A beam finite element for shear-critical RC beams," ACI Special Publications, vol. 237, pp. 295-310, 2006.

[12] D. C. Feng, X. D. Ren, and J. Li, "Softened damage-plasticity model for analysis of cracked reinforced concrete structures," Journal of Structural Engineering, vol. 144, no. 6, 2018.

[13] M. Petrangeli, P. E. Pinto, and V. Ciampi, "Fibre element for cyclic bending and shear of RC structures. I: theory," Journal of Engineering Mechanics, vol. 125, no. 9, pp. 994-1001, 1999.

[14] P. Ceresa, L. Petrini, and R. Pinho, "Flexure-shear fiber beam-column elements for modeling frame structures under seismic loading--state of the art," Journal of Earthquake Engineering, vol. 11, no. 1, pp. 46-88, 2007.

[15] J. Mazars, P. Kotronis, F. Ragueneau, and G. Casaux, "Using multifiber beams to account for shear and torsion. Applications to concrete structural elements," Computers Methods in Applied Mechanics and Engineering, vol. 195, no. 52, pp. 7264-7281, 2006.

[16] J. N. Gregori, P. F. Miguel, M. A. Fernandez, and J. R. Marti-Vargas, "A theoretical model for including the effect of monotonic shear loading in the analysis of reinforced concrete beams," Engineering Structures, vol. 52, pp. 257-272, 2013.

[17] T. R. S. Mullapudi and A. Ayoub, "Analysis of reinforced concrete columns subjected to combined axial, flexure, shear, and torsional loads," Journal of Structural Engineering, vol. 139, no. 4, pp. 561-573, 2013.

[18] A. Correia, J. Pacheco, and R. Pinho, "Force-based higher-order beam element with flexural-shear-torsional interaction in 3D frames," Part I: Theory. Engineering Structures, vol. 89, pp. 204-217, 2015.

[19] J. M. B Garcia and A. R. Mari Bernat, "Coupled model for the nonlinear analysis of anisotropic section subjected to general 3D loading. Part 1: theoretical formulation," Computers and Structures, vol. 84, no. 31-32, pp. 2254-2263, 2006.

[20] S. Mohr, J. M. Bairan, and A. R. Mari, "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.

[21] P. Di Re, D. Addessi, and F. C. Filippou, "Mixed 3D beam element with damage plasticity for the analysis of RC members under warping torsion," Journal of Structural Engineering, vol. 144, no. 6, p. 04018064, 2018.

[22] A. Kagermanov and P. Ceresa, "RC fiber-based beam-column element with flexure-shear-torsion interaction," in High Tech Concrete: Where Technology and Engineering Meet, D. Hordijk and M. Lukovic, Eds., Springer, Cham, Switzerland, 2018.

[23] A. Kagermanov and P. Ceresa, "Physically-based cyclic tensile model for RC membrane elements," Journal of Structural Engineering, vol. 142, no. 12, article 04016118, 2016.

[24] F.J. Vecchio and R. G. Selby, "Toward compression-field analysis of reinforced concrete solids," Journal of Structural Engineering, vol. 117, no. 6, 1991.

[25] A. Kagermanov and P. Ceresa, "Fiber-section model with an exact shear strain profile for two-dimensional RC frame structures," Journal of Structural Engineering, vol. 143, no. 10, article 04017132, 2017.

[26] M. Menegotto and P. E. Pinto, "Method of analysis of cyclically loaded reinforced concrete plane frames including changes in geometry and non-elastic behavior of elements under combined normal force and bending," in Proceedings of the IABSE Symposium, Resistance and Ultimate Deformability of Structures Acted on by Well-Defined Repeated Loads, IABSE, Zurich, Switzerland, 1973.

[27] F. C. Filippou, E. G. Popov, and V. V. Bertero, "Effects of bond deterioration on hysteretic behavior of reinforced concrete joints," Tech. Rep. UCB/EERC-83/19, EERC, University of California, Berkeley, CA, USA, 1983.

[28] E. Spacone, F. C. Filippou, and F. F. Taucer, "Fibre beam column model for non-linear analysis of R/C frames: Part I," Formulation. Earthquake Engineering and Structural Dynamics, vol. 25, no. 7, pp. 71-125, 1996.

[29] N. S. Bousias, G. Verzeletti, M. N. Fardis, and E. Gutierrez, "Load-path effects in column biaxial bending with axial force," Journal of Engineering Mechanics, vol. 121, no. 5, 1995.

[30] W. Onsongo, The diagonal compression field theory for reinforced concrete beams subjected to combined torsion, flexure and axial load, Ph.D. thesis, University of Toronto, Toronto, ON, Canada, 1978.

Alexander Kagermanov (1) and Paola Ceresa [ID] (1,2)

(1) University School for Advanced Studies IUSS Pavia, 27100 Pavia, Italy

(2) RED Risk Engineering + Development, 27100 Pavia, Italy

Correspondence should be addressed to Paola Ceresa; paola.ceresa@iusspavia.it

Received 31 March 2018; Accepted 2 July 2018; Published 15 August 2018

Academic Editor: Fabrizio Paolacci

Caption: FIGURE 1: Beam test in torsion, shear, and flexure [1, 2].

Caption: FIGURE 2: Beam kinematics and section discretization.

Caption: FIGURE 3: Calculated warping functions for rectangular and hollow-core section.

Caption: FIGURE 4: 3D constitutive model reference systems.

Caption: FIGURE 5: Cyclic response of concrete (a) and cracked concrete element for tensile model (b).

Caption: FIGURE 6: Frame element degrees of freedom in the global (a) and basic (b) configurations.

Caption: FIGURE 7: Flow chart summarizing the force-based element iteration procedure.

Caption: FIGURE 8: Imposed top displacement on the columns.

Caption: FIGURE 9: Numerical model of the Bousias column and section discretization.

Caption: FIGURE 10: Lateral force-displacement response of the columns [29].

Caption: FIGURE 11: Horizontal reactions in the columns [29].

Caption: FIGURE 12: Numerical model of the Onsongo beam and section discretization.

Caption: FIGURE 13: Moment-curvature and torque-twist response for beam TBO3 [30].

Caption: FIGURE 14: Moment-curvature and torque-twist response for beam TBU3 [30].

Caption: FIGURE 15: Torque-hoop strains for beams TBO3 and TBU3 [30].

Caption: FIGURE 16: Numerical model and section discretization of the floor subassembly.

Caption: FIGURE 17: Load-deflection and load-twist response of the floor subassembly [1, 2].

TABLE 1: Material properties of the columns [29]. Test unit [f'.sub.c] (MPa) [[epsilon].sub.o] [f.sub.y] (MPa) S1 29 0.002 460 S7 28.1 0.002 460 Test unit [f.sub.u] [[epsilon].sub.u] P/[f'.sub.c] * [A.sub.c] (MPa) ([per thousand]) S1 710 11 0.12 S7 710 11 0.12 TABLE 2: Material properties of the beam in bending and torsion [30]. Beam [f.sub.c] (MPa) [[epsilon].sub.o] T/M TBO3 19.1 0.0024 0.701 TBU3 34.8 0.0031 0.701 TABLE 3: Material properties of the beam in bending, shear, and torsion [1]. Longitudinal steel [phi] 16 Concrete [f'.sub.c] [[epsilon].sub.o] [E.sub.c] [f.sub.y] [E.sub.s] (MPa) (MPa) (MPa) (MPa) 24.1 0.0022 24,000 434 200,000 Transverse steel [f'.sub.c] [phi] (mm) [f.sub.y] (MPa) [E.sub.s] (MPa) (MPa) 24.1 9.5 294 200,000 6.5 310 200,000

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Kagermanov, Alexander; Ceresa, Paola |

Publication: | Advances in Civil Engineering |

Date: | Jan 1, 2018 |

Words: | 5964 |

Previous Article: | Investigations on Engineering Properties of Solidified/Stabilized Pb-Contaminated Soil Based on Alkaline Residue. |

Next Article: | An Analytical Study of Square CFT Columns in Bracing Connection Subjected to Axial Loading. |

Topics: |