AbstractThis study validates the use of available tools for nonlinear finite element dynamic analysis (NLFEDA) modeling of unbonded post-tensioned reinforced concrete (PT-RC) members under drop weight impact and is the first attempt at low-velocity impact to the authors’ knowledge. The study was performed in two parts. The first portion of the study determined the appropriate mesh and material properties of the testing setup and specimens that were generated by the authors. The second portion of the study then established the level of convergence of the model based on a wide range of design and loading parameters. The NLEFDA was validated using experimental data obtained from laboratory investigations on the behavior of shear-critical and flexural PT-RC members. The interaction between the prestress, as applied from the PT tendon, and the surrounding concrete was well modeled in the presence of large shear failures of the concrete. The modeling techniques, coupled with experimental studies, can be helpful in expanding the range of understanding of the behavior of unbonded PT-RC members.IntroductionPrestressed concrete construction is commonly found in structures where horizontal and vertical deformation is a critical aspect of safety and design (Kang and Ibrahim Ary 2012; Lee et al. 2015; Shin et al. 2020). Developed in the late nineteenth century by Eugene Freyssinet, the use of prestressed concrete provided an economical means to help reduce material usage by aiding in the design of longer floor spans and more slender structures (Bondy and Allred 2018). The pre-compressive stresses applied to the concrete are provided by prestressing longitudinal reinforcement, referred to as tendons, which are comprised of individual or bundled hard-drawn wires and bars of high strength alloy steel. In this manner, the strong compressive properties of concrete are further utilized in conjunction with the prestressing tendon and other longitudinal reinforcement. In the 1950s, the development of unbonded tendons for post-tensioned reinforced concrete (PT-RC) design helped to extend the life span of prestressed concrete structures. As a duct is used for the insertion of the PT tendon, this allowed for additional economical aspects with the capability to replace or re-tension the tendon, common for structures such as nuclear containment facilities and bridge piers which must be regularly maintained for structural integrity.Inherently, as the use of unbonded PT-RC construction has become widely adopted, the need to understand member behavior under all loading scenarios is important. This is especially true for dynamic loading conditions as structures such as bridge piers can be subjected to potentially large impacting masses on a recurring basis. To study these scenarios, impact tests on concrete members have been performed under low velocity impact using a drop weight apparatus, where low velocity was defined as 10 m/s or less (Nghiem et al. 2021). This testing is based on the use of a hammer of significant mass dropped from a controlled height to fail the specimen in a three-point bending configuration. When impacted in such a manner, the main concern for concrete members pertains to: (1) large inertial effects that lead to increasing shear failure modes; (2) the influence of stress wave propagation effects from the impacting hammer on reaction forces; and (3) localized damage at the loading point.Many drop weight laboratory investigations have been performed on traditional reinforced concrete (RC) type members to address both flexural- (Kishi et al. 2001; Fujikake et al. 2009; Tachibana et al. 2010; Kishi and Mikami 2012; Nghiem and Kang 2020) and shear-critical members (Kishi et al. 2002a, b; Bhatti et al. 2009; Saatci and Vecchio 2009; Nghiem and Kang 2020). Unfortunately, in comparison, little information has been previously collected on the impact behavior of unbonded PT-RC beams, as only a few experimental studies were readily available (Ishikawa et al. 1998, 2000; Wu et al. 2016). To populate this previously unavailable lab data, a recent comprehensive study was performed on both flexural- and shear-type unbonded PT-RC members (see Nghiem et al. 2021). The specimens in that study had identical longitudinal reinforcement, but varying shear reinforcement ratios were intended to be both statically flexural- and shear-critical (Vn/Pn near or below 1.0), according to the ACI code. The main variables investigated included: increasing levels of effective PT force (Pse) and the application of two different impact energies per each specimen type (F = flexural, S = shear). From that study, it was concluded that: (1) shear failure mechanisms are still important and present even for statically flexural members; (2) modes of shear failure increase with increasing impact velocity; (3) the engagement of the PT tendon increases member resistance and helps restrain cracking at impact as observed with the increased values of the reaction force; and (4) increased values of initial PT forces seem to exhibit increased shear failure mechanisms in the form of increased brittle failure modes and arching mechanisms from the support point towards the impacting point (Nghiem et al. 2021).However, experimental studies are expensive and time consuming, which lead to many limiting factors (Kang and Wallace 2008). These restrictions are especially evident in drop weight testing with physical constraints such as: (1) allowable span lengths due to testing facilities; (2) allowable drop height limited by the capacity of the testing frame and reaction floor; and (3) dimensions of the member cross-section as limited by the guiding rails used for positioning of the drop hammer. Additionally, the overall number of tested specimens capable of being tested is a large economical restraint. To further evaluate the behavior of concrete members under impact without those limitations, the use of nonlinear finite element dynamic analysis (NLFEDA) has provided additional insights. However, proper validations must be performed to prove its consistency with observed real-life phenomena.The purpose of this numerical investigation is to exhibit the capability of using NLFEDA techniques for the analysis of unbonded PT-RC beams under impact. First, a brief review was done to discuss past successful attempts at the modeling of concrete members, namely RC. This is then followed by the development of an NLFEDA model of unbonded PT-RC members. Validations were performed in two key steps: (1) a sensitivity study to determine appropriate member geometries (with respect to mesh size) and material properties; and (2) establishing the overall level of accuracy of the resulting model considering a wide range of design and loading parameters. The experimental study completed by Nghiem et al. (2021) was used for all validation purposes and was the basis for the development of the NLFEDA mesh.Research SignificanceThe intent of this study is to present a suitable and easy to replicate method for the modeling of unbonded PT-RC beams under impact. As prior studies have relied on indirect methods in which temperature-induced shrinkage is applied to the tendon prior to loading, a direct method for the application of the pre-compressive tendon forces is performed. Additionally, past modeling has focused on members with bonded tendons; in contrast, this study focuses on addressing members with unbonded PT tendons, which, in some cases, can be more applicable for structures such as nuclear containment facilities and other high importance structures. Beyond the scope of this study, the modeling methods can potentially be used for predicting or assessing the dynamic behavior of other unbonded PT structures under other similar loading conditions for future code or guide development. Forensic investigation is one of the areas that requires a nonlinear model.Review of Previous Modeling MethodsThe use of specialized NLFEDA software has become increasingly popular. Seen within the previous decade, while still relatively few in number, an increasing amount of NLFEDA studies has been performed on drop weight testing of concrete members (Bhatti et al. 2009, 2011; Cotsovos 2010; Jiang et al. 2012; Adhikary et al. 2012; Jiang and Chorzepa 2015; Yi et al. 2016; Pham and Hao 2017; Zhao et al. 2017). This is due to significant advances in computational power and wider accessibility. However, one downside is that the obtained results are not readily accepted and are always susceptible to criticism without adequate justification. Commercial products, as those provided by LS-DYNA (2007), have been shown to provide accurate results. Developed by the former Livermore Software Technology Corporation (LSTC) with origins beginning in the 1970s (currently owned by ANSYS), LS-DYNA began as the first 3-D capable software for the simulation of structural elements under impact and dynamic events. Its ability to be utilized under a number of different materials and loading conditions has led to its widespread use in industries ranging from automobile, aerospace, construction, and civil engineering, to manufacturing and bioengineering applications. With respect to the low velocity impact of concrete members, its successful use has been validated by authors such as Jiang et al. (2012), Jiang and Chorzepa (2015), and Zhao et al. (2017). However, depending upon the complexity of the model, computational time can still pose a problem and balance must be made in terms of (1) ease of use, (2) required computational power, and (3) model accuracy. Sensitivity is commonly tuned to the geometric mesh size and properties of the materials used. For the modeling of large scale concrete members, two mesh sizes have been commonly incorporated: a relatively coarser mesh in the range of 20–25 mm (Bhatti et al. 2009, 2011; Jiang and Chorzepa 2015; Zhao et al. 2017) and finer meshes measuring around 10 mm (Jiang et al. 2012; Pham and Hao 2017). It was noted that anything smaller than 5 mm only added to the computational time and exhibited no large benefits in model accuracy (Pham and Hao 2017).A critical aspect when performing the numerical modeling of concrete members is the consideration of how to represent member cracking and potential large brittle failure effects. How this is represented varies greatly based on the materials models used in the NLFEDA environment and is largely subjective to each researcher on how they want to interpret member cracking. A commonly applied method is “element erosion.” In this case, the erosion or removal of elements is done to remove exceedingly distorted Lagrangian elements to prevent numerical problems with respect to model convergence. When compared to real world scenarios, erosion can be seen as a way to visualize the effects of target fragmentation, spalling, scabbing of concrete, or penetration during impact. Typical methods for the assessment of element erosion are based on the following physical quantities at integration points of each element: the equivalent plastic shear strain, maximum principal strain, minimum pressure (positive in compression), damage index, and equivalent stress. Normally, the maximum principle strain has been an effective criteria for erosion criteria. Values of erosion are set at some value in reference to when the principle strain exceeds the set threshold, element deletion will occur (i.e., a value of 1.01 would mean 1% erosion over the principle strain). As erosion is a non-physical parameter, the amount of erosion is then adjusted to some measured output, such as reaction force or displacement (as based on the experimental study used for validation). As will be observed, many successful modeling attempts have been performed and the modeling methods can vary significantly from author to author. A brief review of recent modeling is provided below in chronological order.Early Modeling Methods of RC Members by Bhatti et al. (2009, 2011)The study performed by Bhatti et al. (2009) was an early use of LS-DYNA for the simulation of RC members under impact. RC beams were simulated to gain additional insight on the behavior of the shear failure phenomena related to drop weight impact. The model of Bhatti et al. (2009) was capable of predicting impact force and reaction at the supports with an error of 20%. In contrast, displacements were better evaluated with an error margin of 10%. However, as has been observed from other authors, the post-peak behavior of the displacement is not easily replicated, as Bhatti et al. (2009) stated that the free vibrational response was not aligned with the experimental results. Visually, the crack patterns at the surface were easily replicated within their numerical model when compared to the experimental findings used for their validation. An additional study with a similar mesh was performed later by Bhatti et al. (2011). This new study focused on the behavior of lightweight aggregate concrete under drop weight impact. Both of these early studies helped to show the validity of the use of LS-DYNA on modeling the behavior of concrete elements upon impact.Numerical Bonding Methods of Steel Reinforcement to Concrete by Jiang et al. (2012)A thorough sensitivity study was done in which Jiang et al. (2012) looked carefully at two methods for the bonding of steel to concrete. This study was done to more accurately model strain rate effects to investigate the influence of several important behaviors of the concrete material of RC members under impact loads. Their model was validated against the laboratory findings of Fujikake et al. (2009). With respect to the bond between steel and concrete, two methods were studied independently for comparison. Bonding was performed through (1) the use of a common node between the two elements; and (2) by use of the keyword “CONSTRAINED_LAGRANGE_IN_SOLID” to bond the steel reinforcement to the concrete elements.From their results, impact force and displacement correlated well with the experimental findings. However, member deformation behavior was divided into two phases of a pre- and post-peak. Simulations during pre-peak were good but lesser agreement was observed in the post-peak response. This same behavior was observed from the study by Bhatti et al. (2009) and later noted by Zhao et al. (2017). The authors noted that this observation is due to the release of the elastic deformation of RC during impact. In terms of accuracy, Jiang et al. (2012) noted lesser errors for the specimens tested at higher impacts. Comparison of the two constraint methods for the reinforcement to the concrete showed that common nodes provided a somewhat better result than the use of constraining the reinforcement per certain keywords; however, the detail and difficulty to create the mesh is not worth the increased amount of effort. The authors stated “[that] the coupled method can be used as a substitute for the common nodes method [as the] the small difference regarding impact response of RC structures using two methods [can be] neglected.”Numerical Modeling of Prestressed Girders by Jiang and Chorzepa (2015)Jiang and Chorzepa (2015) modeled pre-tensioned prestressed concrete (PC) girders with bonded tendons under impact using LS-DYNA. As of the publication date (2015), few numerical studies on the PC members were documented. This may have been due to the difficulties involved in modeling the prestressing force in an explicit finite element analysis program. This coupled with the inherent challenges in producing an effective model make the modeling of PC members additionally difficult. One of the objectives of Jiang and Chorzepa (2015) was to find an efficient material model to use in terms of strain rate, prestress loss on impact response, and failure modes of the PC members. The authors looked carefully at determining an effective value of erosion for the material model (considering principal strain). Values of erosion were varied from 1.06 to 1.14. They noted that element erosion has a significant influence on the peak impact force. This is due to the fact that the eroded elements directly reduce the contact area and stiffness around the impact zone. In comparison, it was observed that the erosion values have relatively negligible influence on the midspan displacement. They stated that this is attributed to the fact that the displacement is primarily controlled by the overall response of the beam, whereas the peak impact force may be affected by the local response (e.g., spalling) of the beam. For their particular model, an erosion value of 1.06 provided the most accurate results when compared to the experimental findings.Modeling of Shear Members by Zhao et al. (2017)A study by Zhao et al. (2017) was performed to address the lack of studies on shear-deficient RC members. Compared with flexural-critical members, the behavior of shear beams is not well understood, due to a fundamental lack of research on the shear failure and associated cracking mechanisms. Their study included both experimental testing and numerical modeling. Generally, their results were in good agreement with experimental results. However, they noted that models of flexural type beams were better converged than shear failure specimens. Much like previous authors, two phases of the member deformation were noted with high deviation from the experimental results during the second phase or after peak displacements were obtained. This is due to the fact that “shear cracks develop and the beam is basically [separated into distinct parts] even during early stages of impact, therefore the later stage of the response of the beam is actually a non-continuum mechanism problem [introduced using] NLFEDA [methods].” However, ultimately the results were still in good agreement. The authors noted higher accuracy with flexural RC members than when compared to shear-critical members. With respect to the post-peak response, this phase could be less important from a practical or design perspective as it is a post-failure behavior.Modeling MethodsAs discussed above, commercial software LS-DYNA provides an easy-to-use method capable at delivering adequate results and was used for all numerical modeling. The approach used for the generation of the mesh for the specimens and the testing setup are provided in the following subsections. A main focus for this study is then devoted to determining the appropriate model properties, which is in reference to mesh geometries (element size) and value of material erosion for concrete elements. Using the results of Nghiem et al. (2021), validations were performed in two key steps: (1) the first step involved determining the key parameters as based on two key laboratory specimens; and (2) then the simulation of all laboratory specimens with the determined mesh and material properties to find the level of convergence of the NLFEDA model with that of the experimental results.Testing SetupModel geometries were constructed to simulate the testing setup and boundary conditions as established per the drop weight experiment (Nghiem et al. 2021), in which both shear and flexural-critical type unbonded PT-RC beams were arranged in a three-point bending configuration. For this study, shear specimens (type S) were defined as members with minimal transverse reinforcement (required only for construction purposes during the laboratory investigation; Nghiem et al. 2021) with shear to flexural capacities (Vn/Pn) near 1.0 and below. Accordingly, flexural (type F) specimens were constructed and modeled with transverse reinforcement at an on-center spacing of 100 mm and values of Vn/Pn in the range of 4.0–5.3, well above unity. A mock-up of the experimental testing frame is shown in Fig. 1. Specific to this experimental study were the established boundary conditions at the left and right supports, in which the left support was allowed to rotate while rotation and displacement along the horizontal direction were allowed for the right support. Per laboratory testing, this was done to minimize possible effects of catenary action. Additional details of the experimental testing setup and specimens can be found in Appendix A.To replicate similar conditions in an NLFEDA environment, the support was constructed utilizing the following mesh as shown in Fig. 2. The use of 20 mm thick rigid cross plates measuring 60×220 mm at the top and bottom of the member are used to restrict uplifting of the beam. Cross plates were tied through use of rigid beam elements creating a pseudo box for placement of the concrete beam member using keyword “TIED_NODES_TO_SURFACE_OFFSET” to perfectly connect the two elements. Similarly, the bottom cross plates were then tied to cylindrical elements where boundary conditions were applied (using keyword “TIED_SURFACE_TO_SURFACE_OFFSET), defined by allowance for rotation at the left support and rotation plus movement along the horizontal axis at the right support. Lastly, the drop weight hammer was replicated by a 220-mm diameter cylindrical rigid mass equal to 420 kg [see Fig. 2]. The impacting velocity of the hammer was controlled by setting an initial velocity at the start of simulations. The full testing setup and the beam model are given in Fig. 3; please see Figs. B1 and B2 (in the Appendix B) for additional details.Specimen GeometryTo generate the beam specimens, two typical sizes of mesh were generated for evaluation in the sensitivity study, defined herein as a coarse mesh and fine mesh. Figs. B3 and B4 of Appendix B show representative meshes of the concrete. Per Fig. B3 of the Appendix, it can be seen that to form the PT duct, similarly to the experimental specimens, two circular voids were formed at 260 mm offset from the top compressive fiber. Note the variance that exists in element shape and size near the duct area, which was done to account for the PT duct area only. The maximum dimension along any axis of the concrete mesh element was equal to 20 mm and 10 mm for coarse and fine geometries, respectively. The elevation of the mesh, also shown in Fig. 3, was formed using default Type 1 eight node solid elements. All steel reinforcement was formed using 20 or 10 mm beam elements. For the bonded reinforcement, default Type 1 Hughes-Liu beam elements were used and bonded using the keyword “CONSTRAINED_LAGRANGE_IN_SOLID”. This keyword provided the coupling mechanism of the beam elements (rebar, shear stirrups) to the solid concrete elements. Unbonded PT reinforcement was modeled using Type 9 Spotweld beam elements (based on Hughes-Liu formulation). In lieu of any constraint to the solid concrete elements, unbonded conditions were provided using keyword “AUTOMATIC_NODES_TO_SURFACE” to model the unbonded interface or interaction of the tendon to the PT duct.Concrete Material ModelLS-DYNA material “CSCM_CONCRETE” (ID# 159) was used to define the concrete properties. The main consideration for its implementation was due to its ease of use and large amount of user documentation available. This material model has been established by the Federal Highway Administration (FHWA) in its applicable use in dynamic situations and has been previously and successfully used by other authors (Zhao et al. 2017). This model has five main features: (1) use of isotropic constitutive equations; (2) three stress invariant yield surfaces with translation for pre-peak hardening; (3) hardening cap that expands and contracts; (4) damage-based softening with erosion and modulus reduction; and (5) rate effects for increasing strength in high-strain rate applications. A general shape of the model is given in Fig. B5 (of the Appendix). More details on the model can be found elsewhere (LS-DYNA 2007). For this study, details were provided to determine an appropriate threshold for element erosion (as based on principal strain). As will be discussed in the provided sensitivity study, erosion was varied from 1.01 to 1.15 and then correlated with the experimental findings (Nghiem et al. 2021). See Table 1 for input material strengths. The necessary amount of concrete erosion was calibrated while default values were utilized for all other parameters.Table 1. Material properties for concrete and steel (units: MPa)Table 1. Material properties for concrete and steel (units: MPa)Concreteϕ10 (trans.)ϕ20 (long.)ϕ25 (PT bar)fc′fytfyfyp21656464973Steel Reinforcement Material ModelsLS_DYNA material “PLASTIC_KINEMATIC” (ID# 003) was used to define the conventional bonded steel reinforcement, while “SPOTWELD” (ID# 100)” was used for the unbonded PT tendons, where both materials provide isotropic hardening plasticity. “SPOTWELD” material was used due to the ease of applying initial PT forces through the user keyword “INITIAL_AXIAL_FORCE_BEAM”. Initial PT forces were applied during the dynamic relaxation period prior to the application of the impact hammer or dynamic analysis. PT tendons were anchored to the member through the use of rigid plates at both the left and right of the beam member. Connection of the PT tendons to the rigid plates was provided by the keyword “TIED_NODES_TO_SURFACE_OFFSET”, in which the nodes at the left and right were constrained to the surface of the rigid plate. The input material strengths can be seen in Table 1 and are based on the measured lab properties per Nghiem et al. (2021). Typical constitutive behaviors of the two materials are observed in Fig. B6 (of the Appendix).Sensitivity StudyTo determine suitable modeling parameters, a deep sensitivity study was performed based on the two mesh sizes (coarse=20 mm, fine=10 mm) and thresholds of element erosion for concrete (in range from 1.01 to 1.15). These two variables were the basis of the analysis, which ultimately determined the final characteristics of the model. The results were validated using two key specimens from the laboratory investigation (Nghiem et al. 2021), member ID SMPT170 and FMPT180. These specimens were chosen as ideal candidates due to the similarities in the impact energy and applied post-tensioning force; however, both belong to different specimen series, shear and flexural (S & F). Table 2 and Figs. 4 and 5 provide the specimen design and loading conditions. Additional experimental details can be found in Appendix A.Table 2. Design and testing parameters of key specimens used for the sensitivity studyTable 2. Design and testing parameters of key specimens used for the sensitivity studyIDPpt_se (kN)fse (MPa)Vn (kN)Pn (kN)(Vn/Pn)m (kg)v (m/s)Ek (kJ)SMPT170(N)1703463062551.204207.9813.36FMPT180(N)1703671,0362584.024208.2514.31Numerical analyses were performed for both specimens, at coarse and fine meshes (20 mm & 10 mm) with values of erosion equal to 1.01, 1.02. 1.03, 1.04, 1.05, 1.075, 1.10, and 1.15. The resulting plastic strain is shown in Figs. 6 and 7 for both specimens SMPT170(N) and FMPT180(N), respectively. Note the notation “(N)” is used to designate simulated NLFEDA members in lieu of those tested in the lab. Figures refer to models tested with meshes equal to 10 mm. The resulting failure is a representative case for both mesh dimensions (20 and 10 mm). Both shear and flexural models show similar failure patterns in agreement with laboratory testing. Visually, the erosion and plastic strain can be seen as a direct form of the observed laboratory failures.The erosion sensitivity of the NLFEDA model to shear failure is extremely noticeable, as evidenced in the simulation of SMPT170(N). For erosion equal to 1.01 and 1.02, the formation of the shear plug (as shown from distinct diagonal cracking originating from the top compressive fiber) is provided by removing elements, in which a clear triangular pattern is formed. When increasing the tolerance (erosion = 1.03, 1.04 1.05… 1.15) the threshold for deletion increases, thereby showing the failure as large strains in lieu of the deleted elements. When comparing the resulting visual failure to that of FMPT180(N), the differences are clear. This is especially notable at values of 1.01 and 1.02. However, when erosion is greater than 1.02, the two series of specimens (shear and flexural) are not as visually dissimilar.In determining model convergence with the laboratory results, emphasis was placed on member resistance against impacting load. This was taken quantitatively as measured by the reaction force (taken at summation of the left and right supports). This was done as Zhao et al. (2017) and other authors mentioned the ability of LS-DYNA to model member parameters, such as force, much better than member deformation, especially during the initial stages of loading toward the peak response.Table 3 shows the peak reaction values for specimens SMPT170(N) and FMPT180(N) at mesh sizes of both 20 and 10 mm, respectively, and at all cases of erosion (1.01 to 1.15). For the shear model of SMPT170(N), the obtained values at 20 mm mesh were unconservative with values overcalculated when compared to the laboratory results. The resulting maximum reaction force for SMPT170 was equal to 773 kN. For NLFEDA results, this value was in a range of 766 kN at high erosion (=1.15) and around 850 kN (erosion=1.02 and 1.03), with resulting values approximately 10% over the lab-obtained value. The amount of error greater than the lab value is given to the right peak value shown in Table 3. No direct trend was observed when looking at the reaction force quantitatively. The same unconservative behavior was also observed for FMPT180(N), with a range of values from around 900 to 920 kN, compared to the laboratory result of 813 kN. The same amount of error was also observed, around 10% for all erosion values.Table 3. Values of maximum reaction values corresponding to mesh size and allowable thresholds of erosionTable 3. Values of maximum reaction values corresponding to mesh size and allowable thresholds of erosionMesh sizeErosionRmax (kN)Error (%)Rmax (kN)Error (%)20 mm1.01835.738.12914.9312.541.02858.3511.04917.2612.821.03852.8110.32919.2013.061.04833.707.85896.6010.281.05833.937.88903.6211.151.075830.807.48903.4411.121.10833.117.78904.2511.221.15766.76−0.81900.2710.7310 mm1.01632.87−18.13772.38−5.001.02663.50−14.17766.15−5.761.03696.68−9.87641.23−21.131.04593.95−23.16642.89−20.921.05595.12−23.01631.14−22.371.075626.46−18.96641.92−21.041.10610.88−20.97641.76−21.061.15607.30−21.44625.84−23.02In contrast, the use of 10 mm typical mesh dimension provided conservative calculations with results in good approximation with the laboratory study. For SMPT170(N), values were in agreement with −10% to −20% error with larger errors more associated with high values of erosion (>1.04). The same trends are observed from flexural specimen FMPT180(N) with values in the range of −5% and −25% errors. The lower values of error were associated with lower erosion thresholds with values of 1.01 and 1.02.In considering visually observed failures, the use of an erosion threshold equal to 1.02 appeared to be the ideal choice for this particular drop weight test setup. Visually, a threshold of 1.02 showed differences of failure in the two specimens. However, for specimen SMPT170(N), a large number of elements were deleted. This behavior was characteristic of the lab testing, in which high brittle failure occurred. This was in the form of an immediate shear plug at the contact of the hammer, followed by brittle failures at the top and bottom fiber of the left span as detailed by Nghiem et al. (2021). Please note that the use of low erosion may be atypical for members under cyclic loadings due to repeated openings and closings of the cracks. The immediate erosion for this specific case is suitable, as member failure was relatively instantaneous as noted by Nghiem et al. (2021). With a mesh equal to 10 mm, quantitatively, the use of erosion at 1.02 was able to produce conservative calculations with errors of −15% and 5% for SMPT170(N) and FMPT180(N), respectively (see Table 3). As noted by previous authors, further reduction in the mesh size would most likely not lead to more positive results and only increase calculation time and space required for data storage.Model ConvergenceTo access the full characteristics of the model under the governing mesh size and erosion parameters, the NLFEDA model was run for the full matrix of laboratory testing (Nghiem et al. 2021). Convergence of the NLFEDA is determined below, in which the range of errors associated with a wider range of measured parameters is evaluated. Calibration of the erosion was done relative to the peak reaction of the beam. In addition to the two laboratory specimens, SMPT170 and FMPT180, the full testing matrix expanded on the design parameters with increasing levels of initial PT force and impacting energies. Additional details on the testing matrix are given in Appendix A.Shear SpecimensFor shear specimens, two impacting velocities were utilized: approximately 5.5 and 8.0 m/s, pertaining to low and medium velocity (SL & SM), respectively. The resulting failures during the loading period are given in Fig. 8, where the progressive damage of SMPT170(N) is provided as a representative case. Time starts at the time of impact (t=0 ms) to 15 ms after initial impact (time the hammer). Damage of the model can be observed in two items: (1) propagation of the plastic deformation; and (2) progressive erosion of the elements during the propagation of the impact force to the reaction supports. From Fig. 8, typical behavior is governed by: (1) development of the shear plug at the impact point at 0.25 ms after hammer contact; (2) full development of the shear plug at 1.00 ms after impact; (3) arching type shear failure developing at the support points at 2.00 ms after impact; and (4) increasing damage leading up to the peak response at around 15.00 ms after impact. Comparing shear members at lower impact velocity [SLPT172(N) (see Fig. C4 of the Appendix) with that at higher velocity SMPT170(N)], a slightly larger shear plug is observed and the qualitative comparisons are analogous to that as observed in the lab. Note that the slight asymmetry of the damage pattern is due to the asymmetrical testing setup which was also accounted for in the model, in which the left support allowed for rotation while the right support allowed for both rotation and movement along the longitudinal axis. Also, a slight eccentricity between the beam’s midspan and impact location might have existed; however, it was not considered in the numerical modeling.A representative case for comparisons of the time histories for all measured forces is given in Fig. 9, per SMPT170(N). These are of the impact force, reaction force, and PT force. Note that behavior of the impact force is provided in plots of “Peak of Hammer” and “Hammer”; with the former used for observing the short duration and high amplitude of the initial phase, and the latter providing the entirety of the loading. Qualitatively, all parameters show good behavior when compared to the laboratory testing. The impact force shows two phases as established in Nghiem et al. (2021), with a short duration high amplitude phase followed by a more drawn out loading period. However, one thing to note is the slight delay of the peak impact force during the first phase, evident in Fig. 9. This characteristic seems minor when looking at the overall loading history as it is only about a 0.05 ms offset of the peak of NLFEDA and laboratory calculations. Stress wave propagation effects were also well simulated as evidenced in the delay of the initiation of reaction forces from the time of impact (t=0 ms) to the development of support forces (at approximately t=1.5 ms). Behavior of the PT force is also very good but trails off after the post-peak response, noting that the model was more susceptible to large concrete failures or the ability to maintain prestress after impact for shear-critical members. However, member deformation was not as well modeled, as seen in Fig. C21 of the Appendix for SMPT170(N). In contrast, for similar shear members tested at lower impact velocity, SLPT172(N), member deformations were much more in line with lab results (see Fig. 10). The obtained errors of all peak responses of all shear members can be seen in Table 4. Detailed time histories of all members are given in Appendix C.Table 4. Validation of shear-type unbonded members, SL & SM (units: error %)Table 4. Validation of shear-type unbonded members, SL & SM (units: error %)SeriesIDPpeak (%)Rmax (%)Ppt_max (%)δmax (%)SLSLPT120(N)40.46−17.63−4.86−8.01SLPT172(N)49.12−22.33−2.282.34SMSMRC(N)13.99−21.55N/A−44.37SMPT0(N)19.31−14.78−15.82−18.87SMPT57(N)11.73−4.99−25.92N/ASMPT118(N)13.28−14.745.07−43.62SMPT170(N)9.18−14.17−3.93−28.60As seen from Table 4, when comparing the peak values of the model with those of the laboratory, a broad range of agreement was observed across all parameters. Particularly, shear specimens at lower impact velocity (series SL) registered noticeably larger impact forces compared to that in the lab. This amount of error was in the range of 50%. However, the reaction forces and PT forces were both agreed with errors of approximately −20% to −5% and −2% to −10% of the experimental values. When increasing the impact velocity, modeling of series SM registered peak impact forces much more in line with those of the experimental values. This behavior was also observed per Jiang et al. (2012). However, larger differences in the calculated member deformations were observed in the range of approximately 20%–40%. As will be established below, when compared to the modeling of flexural type members, obtaining a high level of accuracy for shear members in an NLFEDA setting is much more difficult with shear members.Flexural SpecimensA representative case of the propagation of failure for flexural specimens is given in Fig. 11, in which the flexural specimen at medium impact velocity, FMPT180(N), is shown. The initial development of the main modes of failure is typically the same as observed for shear models, except that a reduced amount of elemental erosion is largely observed. Additionally, some type of flexural cracking can be seen occurring at 1.00 ms for FMPT180(N) in the form of two vertical cracks from the bottom fiber up towards the impact point. When increasing the impact velocity, model FHPT200(N) showed larger deleted elements or shear failure (see Fig. C45 of the Appendix), which was also characteristic of the lab setting, in which increased failure modes with increased impact velocity were noted.The representative case for comparisons of the time histories for all measured quantities are given in Fig. 12, per FMPT180(N). These are of the impact force, reaction force, and PT force. More similar to the noted observations of shear specimens at medium impact (series SM), both flexural series at medium and high impact velocities, FM and FH, show a broad range of agreement with that of the laboratory results. However, performance of flexural specimens was generally better than shear specimens. This is given by better accuracy of the member deformation. The amount of error was around 20% to 30% of the peak midspan deformation. In the case of some members, errors were in the range of 3%–6%, noting a great level of accuracy. This is in agreement with the results stated by Zhao et al. (2017). While LS-DYNA is capable of producing good results of concrete members, flexural members are more accurately modeled and the sensitivity of the model to shear is high. The range of errors of all peak values can be seen in Table 5. Detailed time histories of all members are given in Appendix C.Table 5. Validation of flexural type unbonded members, FM & FH (units: error %)Table 5. Validation of flexural type unbonded members, FM & FH (units: error %)SeriesIDPpeak (%)Rmax (%)Ppt_max (%)δmax (%)FMFMRC(N)26.01−13.04N/A−35.87FMPT0(N)28.77−13.68−29.72−29.53FMPT60(N)18.88−3.31−12.89−20.16FMPT109(N)34.80−20.32−4.05−18.75FMPT180(N)19.64−5.76−5.99N/AFHFHPT0(N)7.18−4.39N/A−25.36FHPT122(N)22.90−5.36−10.41−3.18FHPT200(N)8.70−3.27−5.31−6.79DiscussionThere are few examples for the numerical modeling of unbonded PT-RC members under low velocity drop weight impact. The study done here was to substantiate the use of readily available commercial software and material models, provided here in the use of LS-DYNA. From the results, the use of LS-DYNA provided a practical method of NLFEDA modeling of unbonded PT-RC members, both at replicating the prestressing force provided by the tendon and also the impacting conditions for both shear- and flexural-critical members. Qualitatively, all measured parameters, such as impact and reaction forces, PT forces, and member deformation. were in good agreement with those of the experimental values, with only minor differences in the distribution. Additionally, all major failure modes were observed in the form of element erosion, noting the occurrence of shear failure common to concrete members under this particular loading. However, some differences were noted with the peak response of all parameters with a broad range of errors dependent upon various design or loading conditions.Nevertheless, the observed differences of the NLFEDA modeling with that of the experimental agreed with those, as noted by previous authors. These details include the: (1) ability to model flexural specimens more accurately than shear specimens; and (2) difficulty at modeling deformation, especially during the post-peak phase. Both of these attributes are associated with the non-continuum problem as stated by Zhao et al. (2017), in which it is noted that the model basically begins to be separated into two or more distinct pieces upon large shear failures. This problem is difficult to avoid, as one key element of the lab study was the extent of shear failure especially for members without transverse reinforcement and the increased brittle nature observed at higher levels of prestress. This led to a particular importance in capturing these characteristics. By introducing low thresholds of element erosion (=1.02) to represent the brittle nature of the unbonded PT-RC members under impact, conservative results were obtained, and the interaction of the PT tendon with the surrounding concrete in presence of large failures was met with similar convergence between both shear and flexural members when compared to the experimental testing.A critical modeled aspect was the loss of residual PT tendon forces after impact. This was one observation of the lab study (Nghiem et al. 2021), in which the brittle nature of shear members directly led to a loss in the ability to maintain PT force after impact. Per Nghiem et al. (2021), it was noted that: (1) increases in impact velocity increased shear failure mechanisms which directly led to inabilities of shear members to maintain prestress; and (2) increases in the value of initial PT potentially led to increased brittle failure mechanisms thereby increasing the potential inability to maintain prestress after impact. For this study, the use of erosion to model the loss of strength with respect to the ability of the member to maintain prestress after impact was on the conservative side. This was observed in the modeling of members SLPT172(N), SMPT118(N), and SMPT170(N), in which erosion of elements contained within the impacting zone (shear plug) leads to the failure of the surrounding concrete to hold PT forces, thereby resulting in PT forces after impact being less than the initial values. Additionally, as noted in the lab setting, this behavior seems to also be attributed to increased levels of prestress. For the modeling of flexural members, this was observed for FHPT200(N), a member with relatively higher levels of initial PT forces and impacted at a higher impact velocity. Per the NLFEDA model, this interaction between the PT tendon and the surrounding concrete in the presence of large failures is a crucial aspect even for flexural members that should be studied more thoroughly.ConclusionsThe use of LS-DYNA for NLFEDA modeling of unbonded PT-RC members under low velocity drop weight impact was validated under a broad range of design and loading conditions. When considering the lab results of Nghiem et al. (2021), a model with a typical mesh size of 10 mm and low thresholds for element erosion (=1.02) of concrete were suitable at providing conservative and practical results. Calculations made per NLFEDA were able to visually show important failures during the loading process and maintain a broad range of validation with respect to impact and reaction forces and member deformation. Importance was placed on the ability to capture the distinct failure modes of shear- and flexural-critical members, with shear models (SL and SM) being able to reflect the extreme brittle nature observed in the experimental testing through the large number of eroded elements in the form of a shear plug at the loading point. Conversely, for flexural models (FM and FH), erosion was less severe but still showed aspects of shear failure combined with flexural failure modes throughout the loading process. Additionally, the interaction between the prestress as applied from the PT tendon was well modeled in the presence of large shear failures of the concrete. With respect to the observed phenomenon and inertial effects of low velocity impact loading, this study substantiated the use of commercial software and material models. The modeling performed in this study can provide beneficial and additional methods that can easily replicate the behavior of unbonded members under low velocity impact. These modeling tools, coupled with experimental studies, can be helpful in expanding the range of understanding of the behavior of such members by capturing data not easily obtained in a lab setting, either due to testing complexity or economic constraints.Data Availability StatementAll data, models, and code generated or used during the study appear in the published article.AcknowledgmentsThis work was supported by the Nuclear Safety Research Program through the Korea Foundation Of Nuclear Safety (KoFONS) using the financial resource granted by the Nuclear Safety and Security Commission (NSSC) of the Republic of Korea (No. 2003007). Partial support by the Institute of Construction and Environmental Engineering at Seoul National University is also appreciated. The views expressed are those of the authors and do not necessarily represent those of the sponsors.References Adhikary, S. D., B. Li, and F. Fujikake. 2012. “Dynamic behavior of reinforced concrete beams under varying rates of concentrated loading.” Int. J. Impact Eng. 47 (Sep): 24–38. https://doi.org/10.1016/j.ijimpeng.2012.02.001. Bhatti, A. Q., N. Kishi, and H. Mikami. 2011. “An applicability of dynamic response analysis of shear-failure type RC beams with lightweight aggregate concrete under falling-weight impact loading.” Mater. Struct. 44 (1): 221–231. https://doi.org/10.1617/s11527-010-9621-9. Bhatti, A. Q., N. Kishi, H. Mikami, and T. Ando. 2009. “Elasto-plastic impact response analysis of shear-failure-type RC beams with shear rebars.” Mater. Des. 30 (3): 502–510. https://doi.org/10.1016/j.matdes.2008.05.068. Bondy, K. D., and B. Allred. 2018. Post-tensioned concrete principles and practice. Lexington, KY: Lulu Press. Cotsovos, D. M. 2010. “A simplified approach for assessing the load-carrying capacity of reinforced concrete beams under concentrated load applied at high rates.” Int. J. Impact Eng. 37 (8): 907–917. https://doi.org/10.1016/j.ijimpeng.2010.01.005. Ishikawa, N., H. Enrin, S. Katsuki, and T. Ohta. 1998. “Dynamic behavior of prestressed concrete beams under rapid speed loading.” In Structures under shock and impact V, 35. Southampton, UK: Computational Mechanics Publications. Ishikawa, N., S. Katsuki, and K. Takemoto. 2000. “Dynamic analysis of prestressed concrete beams under impact and high speed loadings.” In Structures under shock and impact VI, edited by N. Jones, C. A. Brebbia, and A. M. Rejendran, 48. Southampton, UK: WIT Press. Kang, T. H.-K., and M. Ibrahim Ary. 2012. “Shear-strengthening of reinforced & prestressed concrete beams using FRP: Part II—Experimental investigation.” Int. J. Concr. Struct. Mater. 6 (1): 49–57. https://doi.org/10.1007/s40069-012-0005-0. Kang, T. H.-K., and J. W. Wallace. 2008. “Seismic performance of reinforced concrete slab-column connections with thin plate stirrups.” ACI Struct. J. 105 (5): 617. Kishi, N., and H. Mikami. 2012. “Empirical formulas for designing reinforced concrete beams under impact loading.” ACI Struct. J. 109 (4): 509–519. Kishi, N., H. Mikami, and T. Ando. 2002a. “Impact-resistance behavior of shear-failure-type RC beams under falling-weight impact loading.” In Structures under shock and impact VII, 63. Southampton, UK: WIT Press. Kishi, N., O. Nakano, K. G. Matsuoka, and T. Ando. 2001. “Experimental study on ultimate strength of flexural-failure-type RC beams under impact loading.” In Proc., Transactions of the 16th Int. Conf. on Structural Mechanics in Reactor Technology. Berlin: International & American Associations for Structural Mechanics in Reactor Technology (IASMiRT). Lee, S.-H., K.-J. Shin, and T. H.-K. Kang. 2015. “Flexural strengthening of continuous concrete beams using externally prestressed bars.” PCI J. 60 (1): 68–86. https://doi.org/10.15554/pcij60.1-01. LS-DYNA. 2007. Keyword user’s manual, version 971. Livermore, CA: LST Corp. Nghiem, A., C. Demartino, Y. Xiao, and T. H.-K. Kang. 2021. “Impact behavior of unbonded post-tensioned concrete beams.” ACI Struct. J. 118 (1): 201–214. Nghiem, A., and T. H.-K. Kang. 2020. “Drop-weight testing on concrete beams and ACI design equations for maximum and residual deflections under low-velocity impact.” ACI Struct. J. 117 (2): 199–210. Saatci, S., and F. J. Vecchio. 2009. “Effects of shear mechanisms on impact behavior of reinforced concrete beams.” ACI Struct. J. 106 (1): 78–86. Shin, H., T. H.-K. Kang, and J. H. Park. 2020. “Grouted extruded-strand tendons: Friction coefficients and differential individual strand forces.” ACI Struct. J. 117 (3): 223–234. Tachibana, S., H. Masuya, and S. Nakamura. 2010. “Performance based design of reinforced concrete beams under impact.” Nat. Hazards Earth Syst. Sci. 10 (6): 1069–1078. https://doi.org/10.5194/nhess-10-1069-2010. Yi, W.-J., D.-B. Zhao, and S. K. Kunnath. 2016. “Simplified approach for assessing shear resistance of reinforced concrete beams under impact loads.” ACI Struct. J. 113 (4): 747–756. https://doi.org/10.14359/51688617. Zhao, D.-B., W.-J. Yi, and S. K. Kunnath. 2017. “Shear mechanisms in reinforced concrete beams under impact loading.” ACI Struct. J. 143 (9): 04017089.