Microsoft Word - numero_58_art_02_3086 R.N. da Cunha et alii, Frattura ed Integrità Strutturale, 58 (2021) 21-32; DOI: 10.3221/IGF-ESIS.58.02 21 Lumped damage mechanics as a diagnosis tool of reinforced concrete structures in service: case studies of a former bridge arch and a balcony slab Rafael Nunes da Cunha Federal University of Alagoas, Maceió, Brazil rafael.cunha@ctec.ufal.br, http://orcid.org/0000-0003-2503-6758 Camila de Sousa Vieira Laboratory of Mathematical Modelling in Civil Engineering, Federal University of Sergipe, São Cristóvão, Brazil milavieira@gmail.com, http://orcid.org/0000-0002-5371-191X David Leonardo Nascimento de Figueiredo Amorim Laboratory of Mathematical Modelling in Civil Engineering, Federal University of Sergipe, São Cristóvão, Brazil Federal University of Alagoas, Maceió, Brazil davidnf@academico.ufs.br, http://orcid.org/0000-0002-9233-3114 ABSTRACT. Reinforced concrete structures may need repair in order to ensure the designed durability. The assessment of the structural condition usually is made in loco, but sometimes numerical analyses are required as a low cost and effective preliminary diagnosis. Lumped damage mechanics (LDM) may be used as an alternative diagnosis procedure. In LDM approaches for reinforced concrete structures, the concrete cracking in modelled by a damage variable and the reinforcement yielding is represented by plastic deformations. Therefore, this paper presents LDM as a diagnosis tool to analyse actual structures. Then, the procedure is applied in two case studies. The first one is a former bridge arch tested in laboratory. Such bridge was deactivated after almost three decades in service and some damage were observable. The second case study deals with a reinforced concrete balcony that collapsed after fifteen years of service. In the moment of collapse, only the dead load was acting in the balcony. The results show that LDM numerical response of those structures are quite close to laboratory observations, for the former bridge arch, and to in loco measurements, for the balcony. In this sense, according to the analysed structures, LDM may be used as a diagnosis tool. KEYWORDS. Lumped damage mechanics; Reinforced concrete; Diagnosis; Bridge arch; Balcony. Citation: Cunha, R. N., Vieira, C. S., Amorim, D. L. N. F., Lumped damage mechanics as a diagnosis tool of reinforced concrete structures in service: case studies of a former bridge arch and a balcony slab, Frattura ed Integrità Strutturale, 58 (2021) 21-32. Received: 29.04.2021 Accepted: 03.07.2021 Published: 01.10.2021 Copyright: © 2021 This is an open access article under the terms of the CC-BY 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. https://youtu.be/C0VKIqEbuyg R.N. da Cunha et alii, Frattura ed Integrità Strutturale, 58 (2021) 21-32; DOI: 10.3221/IGF-ESIS.58.02 22 INTRODUCTION einforced concrete (RC) structures require periodical inspection, maintenance and repairs to ensure the designed durability. Inspections can identify damages that, if were not corrected, can lead to partial or total collapse of the structure. When a structure shows signs of damage it is necessary to investigate its causes, which results in a diagnosis that serves as the basis for adopting intervention measures. The assessment of structural damage may occur in loco using non-destructive methods such as visual inspection [1], acoustic emission [2] and ground penetrating radar [3], or destructive methods as explosions [4]. This analysis can also take place in laboratory, using samples of an original structure [5] or with the combination of experimental and numerical analysis, aiming to reproduce the experimental behaviour, as the stress/strain and fatigue [6] or the damage and dynamic properties [7]. Other diagnosis techniques are based only on numerical analysis of the structures such as nonlinear dynamic analysis [8] or on artificial neural networks [9]. One of the most used computational method for modelling the occurrence of damage and monitoring the cracking process in reinforced concrete structures is the finite element method (FEM) [10-17]. However, a disadvantage of these procedures using FEM is the significant computational cost. An efficient alternative to such procedures in order to assess damage in reinforced concrete structures is the Lumped Damage Mechanics (LDM). This theory combines the concepts of the classic damage and fracture mechanics in plastic hinges, which defines the concrete cracking as the damage variable. LDM was firstly proposed to analyse RC structures submitted to seismic loads [18-22]. Later, LDM was expanded for several other conditions, as impact loads [23-25] and tunnel linings [26-28]. One of the great advantages of using LDM to evaluate damage of reinforced concrete structures is the use of parameters easily obtained by civil engineers in practice, such as cracking moment and the load conditions. In order to ensure that LDM may be applied in actual engineering problems, Flórez-López et al. [29] defined structure integrity according to the damage level of the structural elements. For reinforced concrete elements, a damage value around 0.3 or 0.4 represents the beginning of the reinforcement yield, and a damage of approximately 0.6 means the bearing capacity of the element [29]. In the light of the foregoing, the objective of this paper is the application of LDM as a diagnosis tool of RC structures, using a former bridge arch and a balcony that suffered collapse as case studies. The first case is an arch that was removed from an actual bridge in China [13] and tested in laboratory and the second one is a collapsed balcony of an actual building in Brazil [30]. Then, the damage variable is used as a diagnosis tool, presenting satisfactory results for both analysed structures. Figure 1: Relations among variables in a structural problem. REVIEW ON LUMPED DAMAGE MECHANICS FOR FRAMES AND ARCHES [26-29] he variables of a structural problem can be categorised in three sets: kinematic, static and internal ones (Fig 1). Such variables are related by three sets of equations, which can be called as kinematic equations, equilibrium relations and constitutive law (Fig 1). Classic finite element analysis usually relates directly the matrix of generalised displacements {U} and matrix of generalised external forces {P}. Note that the approach presented in this paper presents more matrices, defined in this section. Then, considering a structure composed by frame and arch elements in a XZ Cartesian reference [26], such quantities are described by: R T R.N. da Cunha et alii, Frattura ed Integrità Strutturale, 58 (2021) 21-32; DOI: 10.3221/IGF-ESIS.58.02 23        1 1 1 ... ... ... T i i i j j j n n nu w u w u w u wU (1)        1 1 1 ... ... ... T u w ui wi i uj wj j un wn nP P P P P P P P P P P PP (2) where u1 is the node 1 displacement parallel to the X-axis, w1 is the node 1 displacement parallel to the Z-axis and θ1 is the node 1 rotation on the XZ plane. Analogously, Pun is the applied force on node n parallel to the X-axis, Pwn is the applied force on node n parallel to the Z-axis and Pθn is the applied moment on node n on the XZ plane. The structural equilibrium is given by:          1 m T b b b P B M (3) being [B]b the transformation matrix and {M}b the generalised stresses of a finite element b, both given as follows:     Ti j ib M M NM (4)                                   Frame element: sin cos sin cos 0 0 0 ... 1 ... 0 ... 0 0 0 sin cos sin cos 0 0 0 ... 0 ... 1 ... 0 0 0 0 0 0 ... cos sin 0 ... cos sin 0 ... 0 0 0 1 b b b b b b b b b b b b b b b b b b b b b L L L L L L L L node node i node j node n B (5)                                                                                       Arch element: 1 sin cos 0 0 0 ... 1 ... sin sin sin cos 0 0 0 ... 0 ... sin sin cos 1 sin cos 1 cos 0 0 0 ... cos sin 0 ... sin sin cos cos sin ... sin b b b b b b b b b b b b b b b b b b b b b b b b b b b b b node node i R R R R R R B                                                                                                   sin cos cos 0 ... 0 0 0 sin cos cos sin sin cos cos ... 1 ... 0 0 0 sin sin cos 1 sin cos 1 cos ... cos sin 0 ... 0 0 0 sin sin b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b R R R R R R node j node n (6) R.N. da Cunha et alii, Frattura ed Integrità Strutturale, 58 (2021) 21-32; DOI: 10.3221/IGF-ESIS.58.02 24 where Mi and Mj are bending moments on the nodes i and j of the element, respectively, Ni is the axial force on the node i of the element, αb is the orientation of the frame element and Lb its length; βb, χb and Rb geometrically define the circular arch element (Fig 2). Note that the second member of (3) is the sum of the internal forces of each finite element. Figure 2: Generalised stresses and geometric definition of (a) arch and (b) frame finite elements. The kinematic relations for a finite element are given by:          bbΦ B U U (7) being [B(U)]b the same transformation matrix previously presented in (5-6) considering large displacements and {Φ} the matrix of generalised deformations, which is conjugated to {M} and written as:       Ti j ibΦ (8) where φi and φj are flexural rotations due to the bending moments Mi and Mj, respectively, and δi is the axial deformation related to Ni. If small displacements and small deformations are assumed, (7) can be rewritten as follows:       0b bΦ B U (9) being [B0] the transformation matrix for the initial condition. The elastic relations of the element are obtained by assuming the deformation equivalence hypothesis, by applying the Castigliano’s Theorem [26-27, 29] and considering the inelastic phenomena are concentrated at hinges located at the edges of the element (Fig 3):          p bbbΦ Φ F D M (10) where {Φp}b is the matrix of generalised plastic deformations, that accounts for reinforcement yielding, and [F(D)]b is the damaged flexibility matrix, both written as:      0 Tp p pi jbΦ (11)                      0 0 0 11 12 13 0 0 0 12 22 23 0 0 0 13 23 33 1 1 1 1 i b j F F F d F F F d F F F F D (12) R.N. da Cunha et alii, Frattura ed Integrità Strutturale, 58 (2021) 21-32; DOI: 10.3221/IGF-ESIS.58.02 25 being di and dj the damage variables at hinges i and j, respectively, which account for concrete cracking. The terms Fij0 are described in the Appendix for frame and arch elements. Figure 3: Lumped damage modelling for reinforced concrete members [26-29]. In order to complete the constitutive law, the evolution laws of internal variables (damage and plastic rotation) must be considered. Thus, consider the complementary energy of the finite element:     1 2 T p b b W M Φ Φ (13) The evolution law of the damage variables are obtained throughout the Griffith criterion i.e.                                       2 0 11 2 2 0 22 2 0 0 2 1 0 0 2 1 i i i i i ii i i i j j j j j jj j j j d G Y d M FW G dG Y d d d d G Y d M FW G dG Y d d d (14) where Gi and Gj are the damage driving moments for the hinges i and j, respectively, and Y(di) and Y(dj) are the cracking resistance functions [28]. The cracking resistance function for reinforced concrete elements [29] is given by:       0 ln 1 1 i i i d Y d Y q d (15) being Y0 and q model parameters. The plastic evolution laws for each hinge of the element are given by:                                 0 0 0 0 0 10 0 0 0 0 10 0 p i i pi i ip ii i p j j j p j jp jj j f M f C k df f M f C k df (16) R.N. da Cunha et alii, Frattura ed Integrità Strutturale, 58 (2021) 21-32; DOI: 10.3221/IGF-ESIS.58.02 26 being fi and fj the yield functions for the hinges i and j, respectively, and C and k0 model parameters. Note that the model parameters Y0, q, C and k0 can be directly related to the classic theory of reinforced concrete structures, as presented by [20, 27-29]. CASE STUDIES Bridge Concrete Arch [13] he first case study is the failure analysis of two reinforced concrete arch ribs that were part of a 28 years old bridge in China. The geometry of the arches is described in Tab. 1. The failure experiments were made in laboratory where the arches were submitted to static loads [13]. Fig. 4 presents the location of the applied loads (P). Span (L) Height (H) Radius Inclined angle at arch foot Cross section width (T) Cross section height (h) 19.72 m 3.08 m 17.89 m 35.5° 20 cm 28 cm Table 1: Arch geometric parameters [13]. Figure 4: Lumped damage modelling for the arch tested by [13]. The loads were applied at middle (point E), at the left quarter plus 27 cm (point C) and at the right quarter of the arch (point G), which were applied in eighteen load steps. The two end sections of the arches were fixed. The reinforcement is located near the upper and lower sides of the arch cross section. Three bars with 16 mm of diameter composed each layer. Tab. 2 presents the elastic constants of the materials. Material Concrete Steel Elastic modulus E (GPa) 9.8 210 Poisson’s ratio ν 0.17 0.27 Table 2: Material elastic constants [13]. Due the long usage time and problems related to the removal and transportation of the arches, both of them suffered damages at the two ends (nodes A and I). Moreover, at 1/8 distance of the left and right supports (nodes B and H) the arches were also damaged, when they were being dismantling in the field, resulting in two notches. It was obtained that the collapse load (Pmax) was equal to 68 kN and 65 kN for the tested arches. Tab. 3 presents the experimental data of the vertical displacements at the measured location, for different load ratios. T R.N. da Cunha et alii, Frattura ed Integrità Strutturale, 58 (2021) 21-32; DOI: 10.3221/IGF-ESIS.58.02 27 Load P/Pmax (%) Vertical displacement (mm) Node E 0 0 6.76 - 8.46 - 10.15 4.43 11.84 5.02 13.53 5.77 16.91 7.24 27.99 12.63 39.06 17.39 47.91 22.35 52.35 24.96 54.57 26.53 59.41 30.07 68.97 36.4 73.09 39.49 77.5 43.45 90.88 56.43 Table 3: Experimental data of vertical displacements for Pmax = 68kN [13]. Figure 5: Inclined balcony slab [30]. Collapse of a reinforced concrete balcony slab [30] A balcony of a fifteen years old building suddenly collapsed. The balcony structure severely cracked at the clamped end of the cantilevered slab, rotating approximately 17 degrees, but remained connected to the rest of the structure of the building, as is shown in Fig. 5. The balcony floor structure is a reinforced concrete slab with 12 cm of thickness, 9.32 m long and 2.03 m wide (Fig. 6). The cantilevered slab is clamped in the outer edge of the living room slab, whose thickness is 15 cm, and both are simply supported on the facade beam. The bottom of the two slabs are on the same level, so the difference in thickness between them results in a 3 cm step at the top of the cantilevered slab, located next to the beam (Fig. 7). R.N. da Cunha et alii, Frattura ed Integrità Strutturale, 58 (2021) 21-32; DOI: 10.3221/IGF-ESIS.58.02 28 The visual inspection of the fractured slab showed steel bars with 9.5 mm of diameter spaced 150 mm. Analysis of samples of concrete and steel bars resulted in characteristic compressive strength (fck) equal to 30 MPa and yield stress of 500 MPa, respectively. The permanent vertical loads on the slab are concrete dead load (3.0 kN/m2), weight of the bottom coating (0.4 kN/m2), weight of the laying mortar (1.0 kN/m2) and weight of the marble floor (0.56 kN/m2). These loads generate a bending moment in the clamped section equal to -10.55 kNm/m. The flexural design considering only the permanent vertical loads and concrete with fck of 14 MPa, value adopted in the building design, results in a steel area equal to 3.4 cm2/m. This value is lower than the steel area actually utilized on the balcony slab, which is 4.72 cm2/m. Figure 6: Dimensions of the living room and the balcony in cm. Figure 7: Cross-section A-A of the slabs (dimensions in cm). Thus, it was concluded that the accident was caused by the abrupt rectification of the steel bars of the negative reinforcement that were installed with double curvature in the form of a step, due to the rupture of the surface layer of concrete by the action of the straightened bars (Fig. 8). Figure 8: Detail of the steel bars before the accident [30]. RESULTS Bridge Concrete Arch [13] he arch was divided in eight finite elements, as shown in Fig. 4. Although a complete nonlinear analysis is possible [14-29], since the goal of this paper is to present a diagnosis of the analysed structure, a single measurement of the tested arch [13] is taken into account. Then, the procedure was carried out for the step of 10.15% of the collapse T R.N. da Cunha et alii, Frattura ed Integrità Strutturale, 58 (2021) 21-32; DOI: 10.3221/IGF-ESIS.58.02 29 load i.e., experimental value of 6.902 kN, since the collapse load is 68 kN. The vertical displacement at the mid-span of the tested arch is equal to 4.43 mm. The LDM analysis reached applied load of 6.802kN (about 10% of the collapse load) and vertical displacement of 4.334mm for the initial damage values in nodes A, B, H and I equal to 0.25, 0.23, 0.22 and 0.12, respectively. Compared to the experimental values, the absolute and relative displacement errors were equal to 0.096 mm and 2.16%. For the estimated initial damage values in the arch, it is possible to suggest minor repairs at nodes A, B, H and I, considering what is indicated by Flórez-López et al. [29] (Tab. 4), if the arch was not tested up to collapse and then returned to service. Performance level Maximum expected damage in beams Maximum expected damage in columns Description 1 0.30 0.10 The elements do not require any intervention after the event 2 0.40 0.30 Some minor repairs may be needed 3 0.50 0.40 The element requires reparation at reasonably costs 4 0.60 0.50 The structure requires a major rehabilitation process 5 > 0.70 > 0.60 Inacceptable structural behaviour Table 4: Performance level for damage [29]. Collapse of a reinforced concrete balcony slab [30] Since the acting permanent load is 4.96kN/m² during the accident, it was verified that the damage evolution law (Eqn. 14) that G < Y(d), i.e. the loads did not cause the damage in the clamped section. This is corroborated by comparing the cracking moment (15.73kNm/m) with the acting bending moment in the clamped section (10.22kNm/m). To achieve the cracking moment (15.73kNm/m) and the ultimate moment (23.69 kNm/m), the load in the balcony should be 7.63kN/m² and 11.5kN/m², respectively. Based on that and in the conditions observed in loco, the accident occurred due to the action of the straightened bars. In the ultimate condition of the slab section, was observed that the residual thickness was approximately 2 cm, as shown in Fig. 9. The simulation, considering large displacements, reached a damage value of 0.9995 and the horizontal and vertical displacement in the edge of the slab are equal to 13.57cm and 58.99cm, respectively. These values lead to an angle approximately equal to 17.30°, which confirms the in loco observations. Note that the damage obtained in the analysis characterises an inacceptable structural behaviour [29] (Tab. 4), denoting that almost there is no cross section at the clamped end. Assuming that the damage penalises directly the inertia moment [27], the remaining thickness at the clamped end is estimated as 0.95cm, quite close to the field observations. Figure 9: Geometry and typical reinforcement after the accident [25-26]. R.N. da Cunha et alii, Frattura ed Integrità Strutturale, 58 (2021) 21-32; DOI: 10.3221/IGF-ESIS.58.02 30 CONCLUSION his paper presents lumped damage mechanics (LDM) as a diagnosis tool to evaluate damage levels (concrete cracking) in reinforced concrete (RC) structures. Differently from all previous studies, where LDM was used in nonlinear analyses, this paper showed that the LDM formulation can be easily applied in practical engineering problems. Therewith, LDM was applied as a diagnosis tool in two different structures: a former bridge arch in China that was tested in laboratory and a RC balcony slab in Brazil that collapsed after 15 years of service. For the former bridge arch, LDM response was quite close to the experimental one and then initial damage values in four different points were estimate around 0.20, concluding that minor repairs might be needed if the arch was returned to service. On the other hand, the RC slab balcony collapsed with only permanent loads where, the reinforcement straightened up due to an incorrect design. By applying LDM, it was noticed that the analysis reached large displacements as the actual balcony, with a rotation near 17º, quite close to in loco observations. The calculated damage was 0.9995, characterising collapse. Also, the residual thickness is estimated near to 1.0 cm, quite close to in loco measure (2.0 cm). Finally, regarding the analysed engineering problems in this paper, it is possible to observe that the LDM may be adequately a diagnosis tool to evaluate the levels of damage in structures, being able to determine the distribution of initial/current damage values and to verify its acceptability and safety for use. Moreover, it is feasible tool for technical analysis of actual structural accidents, such as the case of the RC balcony slab in Brazil. ACKNOWLEDGEMENTS he first author acknowledges CAPES for the financial support of his D.Sc. studies. The second author acknowledges PNPD/CAPES for the financial support of her Postdoc studies. APPENDIX ust for convenience, the flexibility matrix of a damaged member (12) is reproduced here:                      0 0 0 11 12 13 0 0 0 12 22 23 0 0 0 13 23 33 1 1 1 1 i b j F F F d F F F d F F F F D (A.1) For frame members such matrix is denoted as:                        1 0 1 3 6 1 0 6 1 3 0 0 b b i b b b j b L L d EI EI L L EI d EI L AE F D (A.2) For arch elements the terms of the flexibility matrix are given by: T T J R.N. da Cunha et alii, Frattura ed Integrità Strutturale, 58 (2021) 21-32; DOI: 10.3221/IGF-ESIS.58.02 31                                                         2 0 11 2 2 0 12 2 2 2 2 0 13 2 4 sin 3 2 cos 3sin cos sin cos1 1 2 2sin sin 2 sin sin cos sin cos1 1 2 2sin sin 5sin 5 cos cos 3 2 cos1 2 sin 1 2 b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b R F EI R AE R F EI R AE R F EI                                                         2 0 22 2 2 2 0 23 2 2 3 2 0 33 sin sin cos cos sin sin cos sin cos1 1 2 2sin sin sin cos sin cos1 2 sin sin sin cos cos1 2 sin 3sin cos b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b AE R F EI R AE R F EI AE R F                          2 2 cos 2 3sin cos sin sin sin cos cos sin b b b b b b b b b b b b b b EI R AE (A.3) REFERENCES [1] Rehman, S.K.U., Ibrahim, Z., Memom, S.A. and Jameel, M. (2016). Nondestructive test methods for concrete bridges: a review, Constr. Build. Mater., 107, pp. 58-86. DOI: 10.1016/j.conbuildmat.2015.12.011. [2] Noorsuhada, M.N. (2016). An overview on fatigue damage assessment of reinforced concrete structures with the aid of acoustic emission technique, Constr. Build. Mater., 112, pp. 424-439. DOI: 10.1016/j.conbuildmat.2016.02.206. [3] Solla, M., Lorenzo, H., Novo, A. and Caamaño, J. (2012). Structural analysis of the Roman Bibei bridge (Spain) based on GPR data and numerical modelling, Automat. Constr., 22, pp. 334–339. DOI: 10.1016/j.autcon.2011.09.010. [4] Ruggiero, A., Bonora, N., Curiale, G., De Muro, S., Iannitti, G., Marfia, S., Sacco, E., Scafati, S. and Testa, G. (2019). Full scale experimental tests and numerical model validation of reinforced concrete slab subjected to direct contact explosion, Int. J. Impact Eng., 132, pp. 103309. DOI: 10.1016/j.ijimpeng.2019.05.023. [5] Jesus, A.M.P. de., Silva, A.L.L. da., Figueiredo, M. V., Correia, J.A.F.O., Ribeiro, A.S. and Fernandes, A.A. (2011). Strain-life and crack propagation fatigue data from several Portuguese old metallic riveted bridges, Eng. Fail. Anal., 18(1), pp. 148–163. DOI: 10.1016/j.engfailanal.2010.08.016. [6] Liu, Z., Correia, J., Carvalho, H., Mourão, A., Jesus, A., Calçada, R. and Berto, F. (2019). Global‐local fatigue assessment of an ancient riveted metallic bridge based on submodelling of the critical detail, Fatigue Fract. Eng. Mater. Struct., 42(2), pp. 546–560. DOI: 10.1111/ffe.12930. [7] Masciotta, M.G., Pellegrini, D., Brigante, D., Barontini, A., Lourenço, P.B., Girardi, M., Padovani, C. and Fabbrocino, G. (2019). Dynamic characterization of progressively damaged segmental masonry arches with one settled support: experimental and numerical analyses, Frat. Ed Integrità Strutt., 14(51), pp. 423–441. DOI: 10.3221/IGF-ESIS.51.31. [8] Vásques, J.A., Junemann, R., de la Llera, J.C., M.Eeri, Hube, M.A., M.Eeri, and Chacón, M.F. (2020). Three-dimensional nonlinear response history analyses for earthquake damage assessment: A reinforced concrete wall building case study. Earthq. Spectra, 37(1), pp. 235-261. DOI: 10.1177/8755293020944180. [9] Ukrainczyk, N., Ukrainczyk, V. (2008). A neural network method for analysing concrete durability, Mag. Concr. Res., 60(7), pp. 475–486, DOI: 10.1680/macr.2007.00016. R.N. da Cunha et alii, Frattura ed Integrità Strutturale, 58 (2021) 21-32; DOI: 10.3221/IGF-ESIS.58.02 32 [10] Wu, J., Zhou, Y.M., Zhang, R., Liu, C.B. and Zhang, Z.C. (2020). Numerical simulation of reinforced concrete slab subjected to blast loading and the structural damage assessment, Eng. Fail. Anal., 118, pp. 104926. DOI: 10.1016/j.engfailanal.2020.104926. [11] Abedeni, M. and Zhang, C.W. (2021). Dynamic vulnerability assessment and damage prediction of RC columns subjected to severe impulsive loading, Struct. Eng. Mech., 77(4), pp. 441-461. DOI: 10.12989/sem.2021.77.4.441. [12] Tenorio-Monteiro, E. and Juárez-Luna, G. (2021). Beam-column finite element with embedded discontinuities for modelling damage in reinforced concrete prismatic elements, Structures, 29, pp. 1934-1953. DOI: 10.1016/j.istruc.2020.12.055. [13] Tang, X.S., Zhang, J.R., Li, C.X., Xu, F.H. and Pan, J. (2005). Damage analysis and numerical simulation for failure process of a reinforced concrete arch structure, Comput. Struct., 83, pp. 2609–2631. DOI: 10.1016/j.compstruc.2005.03.017. [14] Ooi, E.T., Yang, Z.J. (2011). Modelling crack propagation in reinforced concrete using a hybrid finite element–scaled boundary finite element method, Eng. Fract. Mech., 78(2), pp. 252–73, Doi: 10.1016/j.engfracmech.2010.08.002. [15] Bakir, P.G., Erdogan, Y.S. (2013). Damage identification in reinforced concrete beams by finite element model updating using parallel and hybrid genetic algorithms, Int. J. Comput. Methods, 10(03), pp. 1350010, DOI: 10.1142/S0219876213500102. [16] Shardakov, I.N., Shestakov, A.P., Glot, I.O., Bykov, A.A. (2016). Process of cracking in reinforced concrete beams (simulation and experiment), Frat. Ed Integrità Strutt., 10(38), pp. 339–50, DOI: 10.3221/IGF-ESIS.38.44. [17] Alañón, A., Cerro-Prada, E., Vázquez-Gallo, M.J., Santos, A.P. (2018). Mesh size effect on finite-element modeling of blast-loaded reinforced concrete slab, Eng. Comput., 34(4), pp. 649–58, DOI: 10.1007/s00366-017-0564-4. [18] Flórez-López J. (1993) Modelos de daño concentrado para la simulación del colapso de pórticos planos, Rev. Int. Mét. Num. Cálc. Dis. Ing., 9(2), pp. 123–139. [19] Cipollina, A., López-Inojosa, A. and Flórez-López, J. (1995). A simplified damage mechanics approach to nonlinear analysis of frames, Comput. Struct., 54(6), pp. 1113–26. DOI: 10.1016/0045-7949(94)00394-I. [20] Perdomo, M.E., Ramírez, A. and Flórez-López, J. (1999). Simulation of damage in RC frames with variable axial forces, Earthq. Eng. Struct. Dyn., 28(3), pp. 311–28. DOI: 10.1002/(sici)1096-9845(199903)28:3<311::aid-eqe819>3.0.co;2-d. [21] Rajasankar, J., Iyer, N.R. and Prasad, A.P. (2009). Modelling inelastic hinges using CDM for nonlinear analysis of reinforced concrete frame structures. Comput. Concr.,6(4), pp. 319–41. DOI: 10.12989/cac.2009.6.4.319. [22] Teles, D.V.C., Cunha, R.N., Amorim, D.L.N.F, Picón, R. and López, J.F. (2021). Parametric study of dynamic behaviour of RC dual system design with the Brazilian Standard Code using the lumped damage model, J. Brazilian Soc. Mech. Sci. Eng., 43(5), pp. 246. DOI: 10.1007/s40430-021-02977-8. [23] Yang, T. -S. and Wang, J. -L. (2010). Damage analysis of three-dimensional frame structure suffering from impact, J. Vib. Shock, 29(12), pp. 177–180. [24] Teles, D.V.C., Oliveira, M.C. and Amorim, D.L.N.F. (2020). A simplified lumped damage model for reinforced concrete beams under impact loads, Eng. Struct., 205. DOI: 10.1016/j.engstruct.2019.110070. [25] Oliveira, M.C., Teles, D.V.C. and Amorim, D.L.N.F. (2020). Shear behaviour of reinforced concrete beams under impact loads by the Lumped Damage framework, Frat. ed Integrita Strutt., 53, pp. 13-25. DOI: 10.3221/IGF-ESIS.53.02. [26] Amorim, D.L.N.F., Proença, S.P.B. and Flórez-López, J. (2013). A model of fracture in reinforced concrete arches based on lumped damage mechanics, Int. J. Solids Struct., 50(24), pp. 4070–4079. DOI: 10.1016/j.ijsolstr.2013.08.012. [27] Amorim, D.L.N.F., Proença, S.P.B. and Flórez-López, J. (2014). Simplified modeling of cracking in concrete: Application in tunnel linings, Eng. Struct., 70, pp. 23–35. DOI: 10.1016/j.engstruct.2014.03.031. [28] Brito, T.I.J., Santos, D.M., Santos, F.A.S., Cunha, R.N. and Amorim, D.L.N.F. (2020). On the lumped damage modelling of reinforced concrete beams and arches, Frat. ed Integrita Strutt., 54, pp. 1-20. DOI: 10.3221/IGF-ESIS.54.01. [29] Flórez-López, J., Marante, M.E., Picón, R. (2015). Fracture and damage mechanics for structural engineering of frames: State-of-the-art industrial applications, Engineering Science Reference. ISBN: 978-1466663794. [30] Lima, N.A. (1996). Empuxo no vazio provoca o colapso da estrutura de uma varanda, In: Acidentes Estruturais na Construção Civil, São Paulo, PINI, 1, pp. 149-154. ISBN: 978-8572660617.