Microsoft Word - numero_73_art_02_5324.docx V. Bomfim et alii, Fracture and Structural Integrity, 73 (2025) 12-22; DOI: 10.3221/IGF-ESIS.73.02 12 A simplified nonlinear model for bamboo-reinforced concrete beams based on lumped damage mechanics Vitória Maria Souza Bomfim, Sergio Luiz de Sousa Nunes, Gilberto Messias dos Santos Júnior Disaster Research Institute, Federal University of Sergipe, São Cristóvão, Brazil vitoriamaria@academico.ufs.br, http://orcid.org/0009-0006-8547-9584 snunes@academico.ufs.br, http://orcid.org/0009-0001-8051-3417 gilbertomessias@academico.ufs.br, http://orcid.org/0009-0009-5107-6148 Camila de Sousa Vieira Campus do Sertão, Federal University of Alagoas, Delmiro Gouveia, Brazil camila.vieira@delmiro.ufal.br, http://orcid.org/0000-0002-5371-191X David Leonardo Nascimento de Figueiredo Amorim Disaster Research Institute, Federal University of Sergipe, São Cristóvão, Brazil Post-Graduate Programme in Civil Engineering, Federal University of Alagoas, Maceió, Brazil davidnf@academico.ufs.br, http://orcid.org/0000-0002-9233-3114 KEYWORDS. Bamboo-reinforced concrete, Nonlinear model, Lumped damage mechanics. Citation: Bomfim, V.M.S., Nunes, S.L.S., Santos Júnior, G.M., Vieira, C.S., Amorim, D.L.N.F., A simplified nonlinear model for bamboo-reinforced concrete beams based on lumped damage mechanics, Fracture and Structural Integrity, 73 (2025) 12-22. Received: 08.01.2025 Accepted: 05.04.2025 Published: 06.04.2025 Issue: 07.2025 Copyright: © 2025 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/kwfzIxqwjK0 V. Bomfim et alii, Fracture and Structural Integrity, 73 (2025) 12-22; DOI: 10.3221/IGF-ESIS.73.02 13 INTRODUCTION amboo is a renewable and sustainable material with high tensile strength that has been used as an alternative to steel as reinforcement of concrete elements [1]. The study of flexural performance of bamboo-reinforced concrete (BRC) beams in terms of load capacity, deflection and failure showed its feasibility [2-3]. However, to design BRC beams, it is necessary to consider the intrinsic characteristics of bamboo, it such as its reduced bond with concrete [4]. The behaviour and strength of BRC beams depend on the bond strength between bamboo and cement, the reinforcement ratio, the application of confinement, the presence of admixtures, and the strength of concrete [5]. Ibrahim et al. [6] tested experimentally bamboo-reinforced concrete beams subjected to flexural loads. The study verified the influence of bamboo’s cross-sectional area on the beam’s load-carrying capacity and the influence of its ultimate tensile strength on deflection. Besides experiments, the behaviour of BRC elements is also analysed through numerical methods. Awoyera et al. [7] validated the experimental evaluation of flexural behaviour of large-scale BRC beams with finite element modelling performed using ABAQUS® software. They demonstrated that members reinforced with 50% bamboo, although with about 14% lesser strength but with minimal deformation and crack propagation, can also be a sustainable alternative for construction. Besides real-life experiments, Mondal et al. [8] utilised finite element numerical experiments to develop a load and resistance factor design framework for BRC beams. They showed that a strength reduction factor to consider the slippage of bamboo inside the concrete could be utilised in the design equation. These papers usually apply two- or three- dimensional finite element analysis, which might require considerable computational effort. Alternatively, lumped damage mechanics (LDM) can be an interesting tool for approaching bamboo-reinforced concrete structures since it is based on key concepts of classic fracture [9-11] and damage mechanics [12]. LDM was originally developed for seismic analysis of conventional reinforced concrete frames [13]. Later, it was developed for different materials and load conditions [14-21]. Recently, LDM was extended to two-dimensional continuum media [22-23] and reinforced concrete slabs [24]. Note that other approaches might also be helpful in analysing reinforced (steel or bamboo) concrete structures, such as continuum damage and cohesive fracture approaches. Regarding continuum damage modelling, concrete damage plasticity (CDP) modelling is quite effective in analysing reinforced concrete structures under different load conditions [25-27]. Another option is to analyse complex concrete structures by cohesive fracture models [28-30]. Regardless of the accuracy of such approaches, lumped damage models may present more efficient simulations. According to Bosse et al. [31], when compared to CDP, lumped damage modelling of reinforced concrete structures demands computational resources approximately 10,000 times lower. Therefore, this paper proposes a novel lumped damage model for bamboo-reinforced concrete beams. The proposed model is easy to implement and feasible for practical applications, especially if several numerical analyses are required, e.g., Monte Carlo simulations on structural reliability. PROPOSED LUMPED DAMAGE MODEL FOR BAMBOO-REINFORCED CONCRETE BEAMS Strain equivalence hypothesis and its application in bamboo-reinforced concrete beams rom classic damage mechanics, the first main concept to analyse is the effective stress. For the sake of simplicity, consider a uniaxial case. If the applied Cauchy stress () implies in damaged material, the effective stress can be defined as:    1  (1) where  is the damage variable. Then, the strain equivalence hypothesis states that an undamaged material can replace the damaged material submitted to a Cauchy stress state with the same strain state submitted by the effective stress (Fig. 1). Therefore, the elasticity law (Hooke’s law) is rewritten using the effective stress, i.e.                       e dE E E E1 1  (2) B F V. Bomfim et alii, Fracture and Structural Integrity, 73 (2025) 12-22; DOI: 10.3221/IGF-ESIS.73.02 14 being E the elasticity modulus and  the total strain. Note that the total strain can be divided into an elastic ( e) and a damaged (d) part. Figure 1: Strain equivalence hypothesis. Now, considering a BRC beam element (Fig. 2), its deformed shape can be described by two relative rotations (i and j) at the elements’ ends. Such relative rotations are assembled in the generalised deformations matrix as follows:      T i jΦ (3) where the superscript T means ‘transpose of’. To consider the aforementioned concepts from classic damage mechanics, the generalised deformations matrix is described as a sum of two parts: elastic {e} and damaged one {d}, i.e.       e dΦ Φ Φ (4) Note that the elastic part of the generalised deformations {e} represents the elastic behaviour of the beam element. On the other hand, the damaged part {d} takes into account the concrete cracking of the BRC beam element (Fig. 2), which is considered lumped at the element’s ends, mathematically represented by inelastic hinges and calculated by damage variables (di and dj). Note that no plastic deformations were considered since bamboo reinforcement does not yield. Figure 2: Generalised deformations (left) and stresses (right) of a beam element and the lumped damage of the inelastic hinges. Therefore, considering an analogous relation as the one presented in (2), the elastic law for a beam element is: V. Bomfim et alii, Fracture and Structural Integrity, 73 (2025) 12-22; DOI: 10.3221/IGF-ESIS.73.02 15                                                                               i i j j i j L dL L EI dEI EI L dL L EI EI EI d L L EI d EI L L EI EI d 0 0 0 3 13 6where , , 0 6 3 3 1 3 1 6 and 6 3 1 e dΦ Φ Φ F M C D M F D M F C D F D (5) being [F0] the elastic flexibility matrix, [C(D)] the compliance matrix due to concrete cracking in the inelastic hinges, [F(D)] the flexibility matrix, L the element’s length, EI the element’s flexure stiffness, and {M} the generalised stress matrix, which contains two bending moments located at the elements’ ends (Mi and Mj):     T i jM MM (6) Note that [C(D)] is null if both damage values are zero, and its main diagonal tends to infinity if both damage values tend to one, which reproduces perfect hinges. Both terms in the main diagonal are obtained directly from the strain equivalence hypothesis (2). Proposed lumped damage model The complementary energy (W) of the BRC beam element is given by:               T TW 1 1 2 2 M Φ M F D M (7) The damage driving moments (Gi and Gj) are obtained by differentiating the complementary energy with respect to the damage variables, i.e.           i i i i j j j j M LW G d EI d M LW G d EI d 2 2 2 2 6 (1 ) 6 (1 ) (8) Since the damage variables represent concrete cracking, the damage evolution laws for each inelastic hinge are defined by a Griffith generalised criterion, such as follows:                   i i i i i i j j j j j j d G R d G R d d d G R d G R d d 0 ( ) ( ) 0 0 ( ) ( ) 0 (9) where R(di) and R(di) are the cracking resistance functions for both inelastic hinges. The cracking resistance function must account for this behaviour since BRC beams tend to be more deformable than conventional ones (steel). Therefore, the proposed cracking resistance function for BRC beams is: V. Bomfim et alii, Fracture and Structural Integrity, 73 (2025) 12-22; DOI: 10.3221/IGF-ESIS.73.02 16              d R d R k d d 5 0 ln 1 exp 1 (10) being R0 and k model parameters. Note that the model parameters present clear physical meaning since if there is no concrete cracking, i.e. d = 0, the parameter k exerts no influence on the cracking resistance; then R0 is defined as the initial cracking resistance. On the other hand, while concrete cracking propagates, i.e. d > 0, the parameter k inserts a second term in the cracking resistance function, which is responsible for increasing the cracking resistance due to the bamboo reinforcement. Finally, the exponential term introduces the necessary deformability for the BRC beams. This characteristic can be observed in a bending moment vs. damage curve, as the one presented in Fig. 3. It is worth noting that there is an ultimate damage value (du) related to the maximum bending moment (Fig. 3). Numerically the bending moment vs. damage curve goes up to d → 1.00; however, the ultimate damage value (du) is considered as the maximum damage value, which is related to the bamboo slippage i.e. collapse. A better explanation about the model parameters (R0 and k) is presented in the Appendix of this paper. Figure 3: Bending moment vs. damage for a bamboo-reinforced concrete beam. Notwithstanding, the model parameters are directly related to BRC properties. For any BRC beam, it is possible to calculate the first cracking moment (Mr) and the ultimate bending moment (Mu). First, by equalling (8) and (10), we have the cracking propagation criterion:            dM L R k d dEI d 2 5 02 ln 1 exp 16 (1 ) (11) Considering the beginning of concrete cracking, M = Mr and d = 0.0 in the previous equation. Then, the initial cracking resistance (R0) is defined as: rM L R EI 2 06 (12) The ultimate bending moment (Mu) is reached for the load-bearing capacity, and the damage presents an ultimate value du (Fig. 3).            uu r u uu dM L M L k d EI dEI d 2 2 5 2 ln 1 exp 6 16 (1 ) (13) Since the previous equation presents two unknown variables (k and du), the second equation is: V. Bomfim et alii, Fracture and Structural Integrity, 73 (2025) 12-22; DOI: 10.3221/IGF-ESIS.73.02 17                           u r u u u u u d d M LM d k d d d d d EI 2 5 5 1 5 4 exp ln 1 exp 0 3 (14) Then, by solving the nonlinear system composed of equations (13-14), the parameters k and du are found. Finite element analysis Regarding the lumped damage framework, the local analysis in the finite element programming usually splits the classic stiffness-based formulation into three matrix equations, as illustrated in Fig. 4. Note that the generalised displacements {U} are obtained throughout usual finite element analysis for any load step. Then, the generalised deformations matrix is calculated by the kinematic relation, expressed as follows:                              L L L L L L L L sin cos sin cos ; where sin cos sin cos Φ B U B (15) where α is the beam’s inclination. The generalised stresses matrix as well as the damage variables can be calculated by the elasticity law (5), being the damage evolution (9) calculated by a simple prediction-correction algorithm. Finally, the external forces {P} are balanced by the element’s internal forces i.e.       T P B M (16) Figure 4: Local finite element analysis. If convergence is achieved, the global analysis moves to the next load step, where the damage variables penalise the element’s stiffness matrix, obtained by substituting (5) and (15) into (16). RESULTS AND DISCUSSION his section compares the proposed model to experiments on BRC beams to analyse its accuracy. BRC beams from Mondal et al. [8] Mondal et al. [8] tested three sets of BRC beams with different reinforcement ratios: 1.5% (BRCB-1), 2.5% (BRCB-2), and 3.7% (BRCB-3). The characteristic concrete compressive strength was equal to 20MPa for all beams. All beams present a T V. Bomfim et alii, Fracture and Structural Integrity, 73 (2025) 12-22; DOI: 10.3221/IGF-ESIS.73.02 18 cross-section of 20cm × 25cm and were submitted to a four-point bending test, where the effective span is 2.7m (Fig. 5). The distance between two consecutive load points is the same, i.e. 90cm. The displacements of all beams were measured at mid-span by an LVDT. By taking advantage of the problem’s symmetry, only two finite elements were needed (Fig. 5). For the tested beams [8], the first cracking moment (Mr) and the ultimate bending moment (Mu), as well as the model parameters R0 and k, are given in Tab. 1. Note that the finite element between the support and the applied load is defined as ‘element #1’ and the finite element between the applied load and the symmetry support is named ‘element #2’. Both elements presented the same ultimate damage value for each beam’s set: 0.8899 for BRCB-1, 0.8900 for BRCB-2, and 0.8901 for BRCB-3. The comparison between the numerical and experimental results is depicted in Fig. 5. For the three sets, the numerical results present good agreement with the experimental envelope, especially considering the BRC beams’ load- bearing capacity. The numerical failure, i.e. when the damage reaches its ultimate value, was quite close to the experiments. Figure 5: Test set-up [8], mathematical model and numerical results compared to the experimental ones. Beams Mr (kN.mm) Mu (kN.mm) Element #1 Element #2 R0 (kN.mm) k (kN.mm) R0 (kN.mm) k (kN.mm) BRCB-1 13,500.0 34,200.0 5.826 -1.795 2.913 -0.897 BRCB-2 15,007.5 46,107.0 7.199 -3.264 3.600 -1.632 BRCB-3 15,606.0 51,115.5 7.785 -4.013 3.892 -2.006 Table 1: Model parameters for the finite element analysis. BRC beams from Mali and Datta [32] Mali and Datta [32] tested several BRC beams under different reinforcement ratios. Only three of those beams showed flexure failure. Such beams have a cross-section of 14cm × 15cm and were submitted to a four-point bending test, where the effective span is L = 1.10m (Fig. 6). By taking advantage of the problem’s symmetry, only two finite elements were used (Fig. 6). The longitudinal reinforcement ratio is 3.8%. For the tested beams [32], the first cracking moment (Mr) and the ultimate bending moment (Mu) are 3,300.00kN.mm and 8,121.67kN.mm, respectively. Therefore, the model parameters are R0 = 0.939kN.mm and k = -0.273kN.mm for the finite element between the support and the applied load, and R0 = 0.469kN.mm and k = -0.137kN.mm for the finite element between the applied load and the symmetry support. Both V. Bomfim et alii, Fracture and Structural Integrity, 73 (2025) 12-22; DOI: 10.3221/IGF-ESIS.73.02 19 elements presented the ultimate damage value of 0.8899. The comparison between the numerical and experimental results is depicted in Fig. 6. Note that the proposed model is well-fitted to the experimental envelope. On the other hand, numerical failure, i.e., when the damage reaches its ultimate value, occurs later than the experimental behaviour. Figure 6: Test set-up [32], mathematical model and numerical results compared to the experimental ones. CONCLUSIONS his paper proposes a novel lumped damage model for bamboo-reinforced concrete (BRC) beams. Its accuracy was verified by comparing it with experimental results of BRC beams submitted to four-point bending tests available in the technical literature. The experimental load vs. displacement curves were compared with the results obtained with the proposed model, where the good agreement of the curves demonstrates the model’s accuracy. The numerical analyses ended when the damage variable reached its ultimate value, characterising the bamboo slippage. The ultimate damage value of each beam was calculated, and it was quite close to the experiments regarding load vs. displacement curves. Note that only two finite elements were necessary to achieve a good agreement with the experiments, showing the model’s feasibility for practical applications. Such characteristics are advantageous when considering structural reliability analysis, such as Monte Carlo. REFERENCES [1] Azuwa, S.B. (2024). A review on experimental observation on structural performance of bamboo reinforced concrete beam, Heliyon, 10, e24628. DOI: 10.1016/j.heliyon.2024.e24628. [2] Govindan, B., Ramasamy, V., Panneerselvam, B. and Rajan D. (2022). Performance assessment on bamboo reinforced concrete beams, Innov. Infrastruct. Solut., 7(16). DOI: 10.1007/s41062-021-00616-8. [3] Bala, A. and Gupta, S. (2023). Engineered bamboo and bamboo-reinforced concrete elements as sustainable building materials: A review, Constr. Build. Mater., 394, e132116. DOI: 10.1016/j.conbuildmat.2023.132116. [4] Mondal, B., Maity, D. and Patra, P.K. (2022). Bond behavior between bamboo and normal-strength concrete: experimental and numerical investigation, Pract. Period. Struct. Des. Constr., 27(3), 04022037. DOI: 10.1061/(ASCE)SC.1943-5576.0000715. [5] Sayed, U., Dauletbek, A., Xin, X., Lorenzo, R. and Li, H. (2022). A review on the mechanical behaviour of bamboo reinforced concrete beams, J. Renew. Mater., 10(12). DOI: 10.32604/jrm.2022.022624. [6] Ibrahim, W., El-Fattah, W., Hassan, H.A., Ehab, A. and Elkarem, A. (2024). Flexural behavior of bamboo concrete beams, Innov. Infrastruct. Solut., 9(371). DOI: 10.1007/s41062-024-01672-6. [7] Awoyera, P.O., Karthik, S., Rao, P.R.M. and Gobinath, R. (2019). Experimental and numerical analysis of large-scale bamboo-reinforced concrete beams containing crushed sand, Innov. Infrastruct. Solut., 4(41). DOI: 10.1007/s41062-019-0228-x. [8] Mondal, B., Maity, D. and Patra, K.P. (2023). Load and resistance factor design for bamboo reinforced concrete beam in ultimate flexural limit state, Struct. Saf., 102, 102323. DOI: 10.1016/j.strusafe.2023.102323. [9] Broek, D. (1984). Elementary engineering fracture mechanics, Dordrecht, Martinus Nijhoff Publishers. T V. Bomfim et alii, Fracture and Structural Integrity, 73 (2025) 12-22; DOI: 10.3221/IGF-ESIS.73.02 20 DOI: 10.1007/978-94-009-4333-9. [10] Tvergaard, V. (2007). Numerical modelling in non linear fracture mechanics. Fracture and Structural Integrity, 1(1), pp. 25–28. DOI: 10.3221/IGF-ESIS.01.03. [11] Taylor, D. (2014). Fracture Mechanics: Inspirations from Nature. Fracture and Structural Integrity, 8(30), pp. 1–6. DOI: 10.3221/IGF-ESIS.30.01. [12] Lemaitre, J. and Chaboche, J.L. (1985). Mécaniques des matériaux solides, 1. ed., Paris, Dunod. [13] 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–1126. DOI: 10.1016/0045-7949(94)00394-I. [14] Faleiro, J., Oller, S. and Barbat, A.H. (2010). Plastic-damage analysis of reinforced concrete frames. Eng. Comput., 27(1), pp. 57–83. DOI: 10.1108/02644401011008522. [15] Alva, G.M.S. and El Debs, A.L.H.C. (2010). Application of lumped dissipation model in nonlinear analysis of reinforced concrete structures, Eng. Struct., 132(4), pp. 974-981. DOI: 10.1016/j.engstruct.2009.12.024. [16] Perdomo, M.E., Picón, R., Marante, M.E., Hild, F., Roux, S. and Flórez-López, J. (2013). Experimental analysis and mathematical modeling of fracture in RC elements with any aspect ratio, Eng. Struct., 46, pp. 407–416, DOI: 10.1016/j.engstruct.2012.07.005. [17] Toi, Y. and Hasegawa, K.H. (2011). Element-size independent, elasto-plastic damage analysis of framed structures using the adaptively shifted integration technique, Comput. Struct., 89(23-24), pp. 2162-2168. DOI: 10.1016/j.compstruc.2011.09.002. [18] Bai, Y., Kurata, M., Flórez-López, J. and Nakashima, M. (2016). Macromodeling of crack damage in steel beams subjected to nonstationary low cycle fatigue. J. Struct. Eng., 142(10), pp. 1–13. DOI: 10.1061/(ASCE)ST.1943-541X.0001536. [19] 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. Fracture and Structural Integrity, 14(54), pp. 1-20. DOI: 10.3221/IGF-ESIS.54.01. [20] Teles, D.V.C., Cunha, R.N., Amorim, D.L.N.F., Picón, R.A. and Flórez-López, J. (2021). Parametric study of dynamic behaviour of RC dual system design with the Brazilian Standard Code using the lumped damage model. J. Braz. Soc. Mech. Sci. Eng., 43, 246. DOI: 10.1007/s40430-021-02977-8. [21] Oliveira, J.M.J., Vieira, C.S., Silva, M.F. and Amorim, D.L.N.F. (2023). Fracture modelling of steel fibre reinforced concrete structures by the lumped damage mechanics: Application in precast tunnel segments. Eng. Struct., 278, 115487. DOI: 10.1016/j.engstruct.2022.115487. [22] Picón, R.A., Santos, D.M., Teles, D.V.C., Amorim, D.L.N.F., Zhou, X., Bai, Y., Proença, S.P.B. and Flórez-López, J. (2021). Modeling of localization using Nash variational formulations: The extended damage mechanics. Eng. Fract. Mech., 258, 108083. DOI: 10.1016/j.engstruct.2024.117993. [23] Nardi, D.C. and Leonel, E.D. (2024). An extended lumped damage mechanics IGABEM formulation for quasi-brittle material failure. Eng. Anal. Bound. Elem., 169(A), 105955. DOI: 10.1016/j.enganabound.2024.105955. [24] Cunha, R.N., Amorim, D.L.N.F., Proença, S.P.B. and Flórez-López, J. (2024). Modeling the initiation and propagation of complex networks of cracks in reinforced concrete plates. Eng. Struct., 308, 117993. DOI: 10.1016/j.engstruct.2024.117993. [25] Costa, P.O.B., Bosse, R.M. and Gidrão, G.M.S. (2022). Behavior assessment of asymmetrical building with concrete damage plasticity (CDP) under seismic load, Fracture and Structural Integrity, 16(61), pp. 108-118. DOI: 10.3221/IGF-ESIS.61.07. [26] Nguyen, Q., and Livaoğlu, R. (2023). Degradation of the first frequency of an RC frame with damage levels. Fracture and Structural Integrity, 17(64), pp. 1–10. DOI: 10.3221/IGF-ESIS.64.01. [27] Correia, V.C., and Barboza, A.S.R. (2024). Numerical Analysis of the Beam-Column Joint on Steel Fiber-Reinforced Concrete Structures. International Journal of Advances in Engineering & Technology, 17(5), pp. 473-487. DOI: 10.5281/zenodo.14172950. [28] Tarasovs, S., Kruminš, J., and Tamužs, V. (2015). Modelling of the fracture toughness anisotropy in fiber reinforced concrete. Fracture and Structural Integrity, 10(35), pp. 271–277. DOI: IGF-ESIS.35.31. [29] Almeida, L.P.R., de Lima Junior, E.T. and Barbirato, J.C.C. (2022). Probabilistic dipole BEM model for cohesive crack propagation analysis. J Braz. Soc. Mech. Sci. Eng., 44, 485. DOI: 10.1007/s40430-022-03765-8. [30] Amorim, D.L.N.F., Picón, R., Vieira, C.S., and Flórez-López, J. (2024). Intra-element versus inter-element crack propagation: the numerical extensometer approach. J Braz. Soc. Mech. Sci. Eng. 46, 360. DOI: 10.1007/s40430-024-04951-6. V. Bomfim et alii, Fracture and Structural Integrity, 73 (2025) 12-22; DOI: 10.3221/IGF-ESIS.73.02 21 [31] Bosse, R.M., Flórez-López, J., Gidrão, G.M.S., Rodrigues, I.D., and Beck, A.T. (2024). Collapse mechanisms and fragility curves based on Lumped Damage Mechanics for RC frames subjected to earthquakes. Eng. Struct. 311, pp.118115. DOI: 10.1016/j.engstruct.2024.118115. [32] Mali, P.R. and Datta, D. (2020). Experimental evaluation of bamboo reinforced concrete beams. J. Build. Eng., 28, pp.101071. DOI: 10.1016/j.jobe.2019.101071. APPENDIX A umped damage models for conventional (steel) reinforced concrete elements usually uses the following relation for the cracking resistance function:         d R d R q d0 ln 1 ( ) 1 (A.1) where R0 and q are model parameters (see e.g. [20] for a brief review). On the other hand, the novel lumped damage model for bamboo-reinforced concrete beams presented in this paper uses (10) as the cracking resistance function, written again as follows just for convenience:              d R d R k d d 5 0 ln 1 exp 1 (A.2) where R0 and k are model parameters. Firstly, note that R0 is the same term for both equations, which can be defined as the initial cracking resistance of the gross cross-section, i.e. R0 is related to the first cracking moment (Mr). The terms q and k from the previous equations are related to the additional cracking resistance introduced by both types of reinforcement. Note that both terms, q and k, are related to reinforcement ratio and resistance. Now, focusing on (A.2), it is well-known that the bamboo reinforcement does not yield and does not reach its ultimate strength since slippage and debonding from concrete usually occur first. Therefore, the parameter k is also related to the ultimate bending moment (Mu) of the bamboo-reinforced cross-section, which accounts for slippage and debonding. Finally, yet importantly, since bamboo-reinforced concrete beams usually present more ductility than equivalent beams with conventional (steel) reinforcement [8], an exponential term based on the damage variable is added to penalise the reinforcement parameter k to reproduce the ductility behaviour of bamboo-reinforced concrete beams [8]. Therefore, for a simple comparison, assume that a bamboo-reinforced concrete beam presents exactly the first cracking moment (Mr) and ultimate bending moment (Mu) of another beam with conventional reinforcement. Once R0, q, and k are properly obtained, the ultimate damage for the conventional reinforced concrete beam is around 0.6-0.65, while the bamboo-reinforced concrete beam presents the ultimate damage at around 0.8-0.9. The initial damage propagation for the bamboo-reinforced concrete beam explains this difference (Fig. A.1). Figure A.1: Bending moment vs. damage for steel- and bamboo-reinforced concrete beams. L V. Bomfim et alii, Fracture and Structural Integrity, 73 (2025) 12-22; DOI: 10.3221/IGF-ESIS.73.02 22 Now, in order to present a better interpretation of the parameter k, the following relation is obtained by rearranging (13):                u r u u u u ML M EI d k d d d 2 2 2 5 6 (1 ) ln 1 exp 1 (A.3) Note that the parameter k is related to Mu, Mr and du. The numerator presents a difference between the square of the effective ultimate bending moment and the square of the first cracking moment. Since du is related to the neutral axis position, both concrete compressive resistance and reinforcement ratio also play a role in parameter k. Finally, differently from the lumped damage models for conventional reinforced concrete elements, the proposed model for bamboo-reinforced concrete beams does not consider plastic deformations. This assumption was made to avoid unnecessary uncertainties since mono-sign cyclic tests were not found, however, it introduces limitations to the proposed model. Nevertheless, the proposed methodology may be useful for practical applications, where the actual behaviour of bamboo-reinforced concrete beams can be estimated with satisfactory accuracy. << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles true /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Error /CompatibilityLevel 1.4 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.0000 /ColorConversionStrategy /CMYK /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams false /MaxSubsetPct 100 /Optimize true /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments true /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile () /AlwaysEmbed [ true ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 300 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.50000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile () /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /CZE /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /HUN /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice