Determination of Johnson–Cook Material and Failure Model Constants for High-Tensile-Strength Tendon Steel in Post-Tensioned Concrete Members

: The estimation of damage in steel tendons is important for evaluating the remaining capacity of existing tensioned members. This research focuses on calculating Johnson–Cook (LC) model and damage parameters of high-strength steel material through quasi-static and dynamic uniaxial tests. Finite element analysis is used to replicate the experimental procedure, and through dynamic image correlation analysis, the numerical results accuracy is veriﬁed. In this investigation, it is found that the JC model can accurately replicate deformation and stress concentration under different strain rates and triaxiality conditions and, thus, can be used for fracture analysis of prestressed concrete members. equivalent plastic strain used to model softening behavior of the material in this work under reference conditions.


Introduction
Clear understanding of damage and deformation caused in high-strength steel tendons due to plastic deformation, strain rate and stress concentration effects is of critical importance for the accurate evaluation of existing tensioned concrete members. During regular operation and function of tensioned members, tendons should remain elastic or well within their serviceability state. Despite that, in recent history, due to corrosion, damage on tensioned members of bridge girders has been reported by several researchers [1][2][3][4][5][6]. Corrosion-induced failure of tendons is not thoroughly included in typical design practices of prestressed members and may cause the collapse of bridges in which they are most commonly installed. In current practice, manual inspection along the length of tendons where corrosion is suspected to have occurred is necessary. Corrosion of tendons in these applications commonly occurs due to two primarily two reasons. On the one hand, cracks are formed in surrounding concrete or due to the prestressed member being in close proximity to chloride-contaminated water. On the other hand, tendon rupture at inadequate grouting positions in the sheathing occurs due to the ingress of water from anchoring points. Once corrosion has occurred in a tendon, depending on the severity, it can cause significant reduction in the cross-sectional area of the member [7] and thus induce an unfavorable stress concentration condition. Furthermore, in cases of severe corrosion, the remaining cross-sectional area might be inefficient for bearing the imposed service dead loads and might yield or fracture ( Figure 1). In the latter case, depending on the condition of the surrounding tendons, complete failure might occur if during fracture and dynamic loading, redistribution capacity of the remaining tendons is exceeded. That is why the investigation of fracture performance of high-strength tendons under high strain rates and stress concentration conditions is important to properly evaluate the remaining capacity of a tensioned member. Once the state of corrosion in a tensioned member is verified, risk assessment analysis can be performed, and thus, the remaining capacity can be simulated with the aid of finite element (FE) tools. One of the most common tools used by researchers [8][9][10][11] to evaluate the performance and characteristics of metals under coupled stress and high strain-rate conditions is the Johnson-Cook (JC) material model with damage [12] and it is currently one of the most widely incorporated models in commercial FE software packages, due to its ability to predict material behavior with accuracy and speed, and because it couples a flow stress model with strain rates, elevated temperatures as well as stress concentration conditions. K. Xu et al. [13] performed an experimental investigation of seven high-strength steels in an effort to propose a modification to the traditional JC constitutive model. In their research, BH300, HSLA350, 440 W, HSS590, TRIP590, DP600 and DP800 steel materials were used in uniaxial tensile tests with strain rates ranging from 0.005 s −1 to 1000 s −1 , at normal environmental temperatures, as the heat-related material softening was not the primary objective of the research. Results were used to calibrate material parameters of the traditional JC constitutive model as well as for the evaluation of resulting differences between their proposed model, the traditional JC model, and the experimental data. From their research, it was found that from strain ranges of 2-15% and a tensile strength of 450-850 MPa, the proposed model had an average error of 2%, which is acceptable.
K. Vedantam et al. [14] investigated the mechanical response of two types of steel, Mild and DP590, in tension, at room temperature, using quasi-static and split Hopkinson bar techniques at strain rates ranging from 0.001 s −1 to 1800 s −1 , and the resulting data were used to calculate the JC model parameters. It was found that for increasing strain rates, fracture strain as well as ultimate stress values increased in a similar manner approaching ultimate stress values of 1000 MPa. Finally, detailed JC material parameters were presented.
From the performed literature investigation, it was clearly identified that to accurately model the fracture behavior of high-strength steels, proper material definition and model calibration is required. The data necessary for the aforementioned FE modeling need to be obtained through expensive and time-consuming experimental effort under both high strain rates and high-stress concentration conditions to accurately consider both damage initiation and progression parameters. For high-tensile-strength tendon material commonly used in Japanese infrastructure in general, readily available model constants are not available. In this work the overall behavior, including plastic deformation and the fracture characteristics of medium-carbon high-strength steel used in tensioned members, has been studied through extensive experimental analysis under quasi-static and medium strain-rate loading conditions as well as stress concentration through the implementation of the tensile testing of notched specimens. Failure parameters and material constants for the JC model under room temperature have been calculated through the analysis of experimental data. Damage growth parameters are also introduced and proposed for accurate modeling of necking and fracture of tensile specimens under uniaxial loading conditions. The constants have been evaluated through numerical modeling of dog-bonetype tensile specimens under uniaxial loading conditions for strain rates similar to the experimental configuration as well as implementing digital image correlation (DIC) for the verification of strain propagation during different strain-rate and stress concentration conditions.

Materials and Experimental Procedure
In this research, the SBPR 930/1080 Type B No. 1 medium carbon high-strength steel was investigated, and its chemical composition is presented in Table 1 (in wt%). The material used for the manufacturing of the specimens was supplied by a local company according to JIS G 3109 [15], in 450 mm × 32 mm cylindrical pieces out of which the dog-bone type specimens were manufactured using a manual lathe. Geometrical details can be seen in Figure 2, for cylindrical specimens used throughout this research for both quasi-static and dynamic tensile tests, with smooth as well as notched gauge lengths. For the tensile tests, MTS 244.11 servo-hydraulic actuator ( Figure 3) was used. Acceleration, velocity, displacement, and excitation frequency characteristics are represented in Figure 4. Despite the dynamic characteristics of the experimental procedure, experimental parameters were well within the capability envelope of the utilized actuator. The actuator was mounted on a loading frame using ball-joints for both fixed and extendable part of the piston to allow for increased mobility and flexibility under various testing conditions and requirements. For this research, both specimens and actuator were positioned and fixed in a vertical orientation to ensure an inline application of pulling force ( Figure 5). Furthermore, to ensure constant pulling rate, the loading end of the specimens were constructed in such a way as to allow for initial retraction of the actuator's piston without exerting force onto the specimen. Once required velocity is achieved, and after that stage, contact between the top of the specimen and the mounting fixture (plate) at the movable end of the actuator is initiated, transferring the resulting load onto the specimen body. As can be seen in Table 2, for dynamic loading cases, three specimens were tested for each required strain rate, and the resulting force-displacement data were converted into true stress-strain data using standard equations for uniaxial tensile tests. To further aid the calibration of the initial part of flow stress-strain curves, 2 mm strain gauges were attached using adhesive onto the specimens, and strain data were captured for values up to 20,000 µm, before the adhesive failure resulting in detachment of the strain gauges positioned in the middle of the respective gauge lengths.
From the flow stress-strain data of smooth (no notch) specimens (Figure 6a), it can be seen that for the same strain rate, similar stress-strain curves were obtained. For quasi-static strain rate (0.001 s −1 ), this specific material showed the most ductile behavior fracturing, on average at 0.127 strain. In uniaxial dynamic tensile (Figure 6b) tests, it can be seen that for the case of 0.5 s −1 strain rate, material showed a more brittle behavior while fracturing, on average at 0.092 strain. In comparison to the 0.5 s −1 strain rate case, specimens tested at 1 s −1 showed an increase in the ductility fracturing, on average at 0.099 strain, as well as showing the highest ultimate stress values overall. Finally, for the case of 2 s −1 strain rate, specimens presented most brittle behavior out of all cases, while not showing any significant increase in their ultimate stress value when compared to the slowest dynamic case. In Figure 7, engineering stress-strain (initial cross-sectional area is used for the strain calculations) and converted true stress-strain (strain calculations are performed based on the actual cross-sectional area, which changes with time) data can be seen for quasi-static  Table 2. Similarly to the case of variable strain rate data, in the case of changing the notch size and radius, good agreement is observed between similar cases. Smooth specimens were the most ductile, and the fracture strain progressively reduced as the minimum specimen radius and notch radius decreased. Appl        From the flow stress-strain data of smooth (no notch) specimens (Figure 6a), it can be seen that for the same strain rate, similar stress-strain curves were obtained. For quasistatic strain rate (0.001 s −1 ), this specific material showed the most ductile behavior  est dynamic case. In Figure 7, engineering stress-strain (initial cross-sectional area is used for the strain calculations) and converted true stress-strain (strain calculations are performed based on the actual cross-sectional area, which changes with time) data can be seen for quasi-static testing of notched specimens listed in Table 2. Similarly to the case of variable strain rate data, in the case of changing the notch size and radius, good agreement is observed between similar cases. Smooth specimens were the most ductile, and the fracture strain progressively reduced as the minimum specimen radius and notch radius decreased. Notable reduction in both ultimate stress and failure strain for the C7_0.001R0.4 and C8_0.001R0.4, was observed, resulting in a considerably brittle fracture behavior in comparison to the smooth cases.  est dynamic case. In Figure 7, engineering stress-strain (initial cross-sectional area is used for the strain calculations) and converted true stress-strain (strain calculations are performed based on the actual cross-sectional area, which changes with time) data can be seen for quasi-static testing of notched specimens listed in Table 2. Similarly to the case of variable strain rate data, in the case of changing the notch size and radius, good agreement is observed between similar cases. Smooth specimens were the most ductile, and the fracture strain progressively reduced as the minimum specimen radius and notch radius decreased. Notable reduction in both ultimate stress and failure strain for the C7_0.001R0.4 and C8_0.001R0.4, was observed, resulting in a considerably brittle fracture behavior in comparison to the smooth cases.  Notable reduction in both ultimate stress and failure strain for the C7_0.001R0.4 and C8_0.001R0.4, was observed, resulting in a considerably brittle fracture behavior in comparison to the smooth cases.

Johnson-Cook Model
The Johnson-Cook model is able to accurately analyze and predict stress-strain behavior for ductile materials, and its applicability as well as accuracy has been studied thoroughly in the literature for steel or aluminum alloys under combined conditions of large deformation, high strain rate as well as elevated temperatures focusing on metal forming or impact performance [16][17][18][19]. The JC stress model is expressed in Equation (1).
in which σ represents von Mises or equivalent stress, A is the yield stress of tested material under reference conditions (strain rate and temperature), B is the strain hardening constant, n is the strain hardening coefficient, C is a coefficient of strain rate resulting in postyield strengthening of the material, . ε * is the dimensionless strain rate, T * homologous  (2) and (3).
where . ε p is the accumulated plastic strain and . ε 0 is the reference strain rate which, in this work, was taken as 0.001 s −1 . T m is defined as the melting temperature of the material and T ref is the reference temperature. For the scope of this research, performance of the highstrength medium carbon steel material was investigated under quasi-static and medium dynamic strain rates as well as varying stress concentration conditions, but the temperature factor was not considered since the primary failure factor of tensioned members is usually due to corrosion, as mentioned in previous sections.

Determination of Material Constants A, B, n
ε 0 and T = T re f in Equation (1), the second and third parentheses are omitted, since the effects of strain rate strengthening and thermal softening are neglected. Modifying the remaining terms by taking the natural logarithm on both sides and using the averaged true stress-strain data from cases C1_0.001_NR, C2_0.001_NR and C3_0.001NR and linearly plotting the ln(σ − A) term with of Equation (4), a linear regression model was used to fit the data points, as can be seen in Figure 8.

Determination of Material Constant C
For the purpose of this work and while not considering thermal softening effects, Equation (1) can be modified as: To obtain the C parameter, stress-strain data at four different strain rates (0.001 s −1 , 0.5 s −1 , 1 s −1 , 2 s −1 ) were used to plot Figure 9, while utilizing A, B, n constants that were calculated in the previous section and substituted in Equation (5). Afterwards, first-order linear fitting was performed using a vertical axis intercept value of 1, since Equation (1) is in the form of y a bx = + (C parameter's calculation sensitivity analysis is presented in Appendix B). Similarly to Figure 8, from the slope of the linear regression fit, the C parameter was calculated, which can also be seen in Table 3.  The A parameter was also calculated under reference strain conditions using the 0.2% offset method. For the linear fitting presented in Figure 8, an R 2 factor of more than 97.5% was achieved, resulting in a good accuracy of the regression model. As a result of the latter, material constants A, B and n were calculated from the slope and intercept of the curve, which can be seen in Table 3.

Determination of Material Constant C
For the purpose of this work and while not considering thermal softening effects, Equation (1) can be modified as: To obtain the C parameter, stress-strain data at four different strain rates (0.001 s −1 , 0.5 s −1 , 1 s −1 , 2 s −1 ) were used to plot Figure 9, while utilizing A, B, n constants that were calculated in the previous section and substituted in Equation (5). Afterwards, first-order linear fitting was performed using a vertical axis intercept value of 1, since Equation (1) is in the form of y = a + bx (C parameter's calculation sensitivity analysis is presented in Appendix B). Similarly to Figure 8, from the slope of the linear regression fit, the C parameter was calculated, which can also be seen in Table 3. To obtain the C parameter, stress-strain data at four different strain rates (0.001 s −1 , 0.5 s −1 , 1 s −1 , 2 s −1 ) were used to plot Figure 9, while utilizing A, B, n constants that were calculated in the previous section and substituted in Equation (5). Afterwards, first-order linear fitting was performed using a vertical axis intercept value of 1, since Equation (1) is in the form of y a bx = + (C parameter's calculation sensitivity analysis is presented in Appendix B). Similarly to Figure 8, from the slope of the linear regression fit, the C parameter was calculated, which can also be seen in Table 3. The material constants that were calculated from the aforementioned constitutive equations for the JC model are summarized in Table 3.  The material constants that were calculated from the aforementioned constitutive equations for the JC model are summarized in Table 3.

Johnson-Cook Damage Model Parameters
Substituting the material constants from Table 3 into Equation (1), the following relationships can be formed according to the JC model for stress, strain and strain deformation rate, as can be seen in Equation (6). When comparing experimental data with analytical prediction from Equation (6), good accuracy can be observed until the onset of damage and necking of the tensile specimen ( Figure 10).
(MPa) (6) mation rate, as can be seen in Equation (6). When comparing experimental data with analytical prediction from Equation (6), good accuracy can be observed until the onset of damage and necking of the tensile specimen ( Figure 10).  To accurately simulate damage that occurs in the material model with regard to the JC damage parameter setting, it is important to define at which point damage is calculated. In this study, after careful consideration (Appendix A), the authors decided to use the Damage Initiation point of Figure 10. With that, the JC damage model is used to relate fracture strain with stress triaxiality ratio, strain rate as well as temperature [11,20], and it is expressed in Equation (7).
D1 to D5 represent damage constants for the JC model, m σ is the mean stress (hydrostatic) and eq σ is the equivalent stress (von Mises). As damage occurs in an element governed by JC damage model, it is accumulated based on a damage law and can be represented by Equation (8) [21]. When damage occurs, during high levels of deformation, material strength is reduced [11] and the resulting relation for stress during this damage evolution step is presented in Equation (9). To accurately simulate damage that occurs in the material model with regard to the JC damage parameter setting, it is important to define at which point damage is calculated. In this study, after careful consideration (Appendix A), the authors decided to use the Damage Initiation point of Figure 10. With that, the JC damage model is used to relate fracture strain with stress triaxiality ratio, strain rate as well as temperature [11,20], and it is expressed in Equation (7).
D 1 to D 5 represent damage constants for the JC model, σ m is the mean stress (hydrostatic) and σ eq is the equivalent stress (von Mises). As damage occurs in an element governed by JC damage model, it is accumulated based on a damage law and can be represented by Equation (8) [21]. When damage occurs, during high levels of deformation, material strength is reduced [11] and the resulting relation for stress during this damage evolution step is presented in Equation (9).
where ∆ε is the equivalent plastic strain increment and ε f is the equivalent strain to fracture under certain stress, strain rate and temperature conditions.
In Equation (9), σ D is the resulting stress after damage in an element has been initiated, and D is a damage parameter with the following conditions (0 ≤ D ≤ 1). In Equation (7), σ m σ eq can also be defined as stress triaxiality ratio η * [22,23], and along with equivalent stress can be obtained from undamaged material, while considering plastic deformation up until the onset of necking. According to the work of Bridgman [24], stress triaxiality values can be estimated from uniaxial tests of round specimens according to the analytical model which is presented in Equation (10).
In the model, η * is the stress triaxiality state value, R represents the radius of notch that the specimen is manufactured according to and α represents the minimum cross sec-tion's radius. Triaxialities calculated according to Bridgman's model for different notched specimens can be seen in Table 2.
Neglecting the effects of strain rate and temperature, Equation (7) can be simplified representing fracture strain in terms of the aforementioned D 1 to D 3 damage parameters and stress triaxiality ratio effects. When plotting the fracture strain-stress triaxiality ratio (Figure 11), experimental tension data [25] is used, with 0.001 s −1 strain rate for smooth and notched specimens from Figure 7b and Table 2, in the form of y = A + B · exp(R0 · x); D 1 to D 3 damage parameters can be calculated from the exponential coefficients of the equation similar in principal to the derivation of Equation (4). Afterwards, the D 4 strainrate-dependent parameter was calculated by rewriting Equation (7) complete with the previously calculated D 1 to D 3 damage parameters according to Equation (11). In detail, when plotting Figure 12) and using a linear regression fitting equation intercepting the vertical axis at a value of 1.0, from the slope of the resulting equation, the final JC damage parameter D 4 was calculated.  The calculated JC damage model parameters can be seen summarized in Table 4 and The calculated JC damage model parameters can be seen summarized in Table 4 and can be used in FE software to simulate yield and fracture of high-strength tendons in prestressed concrete applications.  Figure 11. Fracture strain and stress triaxiality relationship from uniaxial tensile test data. The calculated JC damage model parameters can be seen summarized in Table 4 and can be used in FE software to simulate yield and fracture of high-strength tendons in prestressed concrete applications.

Numerical Simulation
The purpose of the numerical modeling was to verify the reproducibility of the experimental results in the commercially available finite element software Abaqus Explicit [26], as part of a broader research work aimed at modeling the dynamic fracture behavior of prestressed concrete members. To accurately reproduce the dynamic fracture effects of high-tensile-strength tendons, non-linear dynamic analysis was used throughout the

Numerical Simulation
The purpose of the numerical modeling was to verify the reproducibility of the experimental results in the commercially available finite element software Abaqus Explicit [26], as part of a broader research work aimed at modeling the dynamic fracture behavior of prestressed concrete members. To accurately reproduce the dynamic fracture effects of high-tensile-strength tendons, non-linear dynamic analysis was used throughout the simulation process of uniaxial tensile tests. Full-scale, three-dimensional models were created, accurately reproducing the geometrical properties of the constructed specimens shown in Figure 2. In both single element and full-scale analysis, geometric nonlinearity related effects were taken into account by using the "NLgeom" option available in Abaqus Explicit.

Numerical Simulation of Singular Finite Element
Initially, to verify the accuracy of the analytical model in FE simulations, a single 8-noded cubical-shaped C3D8R element measuring 1 mm 3 was modeled, and suitable boundary conditions were applied to simulate biaxial symmetry in the two axes perpendicular to the loading direction. To maintain a stress triaxiality ratio η * of 0.333 throughout the tensile test, the bottom four nodes were restrained in the direction of applied force, and four nodes on each of the two faces perpendicular to the loading axis had their movement in the two orthogonal axes restrained, respectively, as can be seen in Figure 13a. The top four nodes were free to move in the direction of loading, and to achieve that a velocity-based loading condition was applied.
Similarly to the experimental procedure, loading speeds of 0.1-200 mm/s were applied to the top four nodes simultaneously. To reduce inertia-related effects at the beginning of the simulation, velocity amplitude was smoothly applied to the simulation for the first 1/10th of the overall step's duration and then kept constant until the completion of each test (Figure 13b). Duration of the tensile phase of the simulation was adjusted each time according to the required strain rate in order to allow for sufficient simulation time and up until the complete damage being registered at the tested element.
In Abaqus, several ductile material models are available that can accurately capture the deformation of steel materials, but in this work, the JC flow stress model and, correspondingly, the JC damage model was utilized using material parameters that were calculated in previous sections. Along with the data presented in Tables 3 and 4, for this material, according to the manufacturer's specifications, Young's Modulus E = 210 GPa and Poisson's ratio of ν = 0.28, and furthermore, the density was set as ρ = 7.85E − 09 tonnes/mm 3 .
Similarly to the experimental procedure, loading speeds of 0.1-200 mm/s were applied to the top four nodes simultaneously. To reduce inertia-related effects at the beginning of the simulation, velocity amplitude was smoothly applied to the simulation for the first 1/10th of the overall step's duration and then kept constant until the completion of each test (Figure 13b). Duration of the tensile phase of the simulation was adjusted each time according to the required strain rate in order to allow for sufficient simulation time and up until the complete damage being registered at the tested element.  In Figure 14, a comparison between the experimental true stress-strain curves for different strain rates of Table 2 with the results obtained from a singular finite element are presented. Overall, good accuracy was achieved between experimental and numerical results with the exception of 1 s −1 strain rate results in which although the failure strain was similar, the ultimate stress value had a difference of approximately 8%. It is believed that this is due to the non-monotonic nature of failure strain as well as ultimate stress that was observed during the experimental procedure (Figure 6b).
To accurately model the material softening behavior past the damage initiation point [27] represented in Figure 15a, in which σ y0 is the yield stress, ε pl 0 is the equivalent plastic strain at the damage initiation point, and ε pl f is the equivalent plastic strain at failure when the scalar damage parameter D = 1 (Equation (8)), damage and strain are correlated. For modeling of post-damage-initiation softening until element failure data, post-peak stress σ D is calculated based on the difference between experimental and JC analytical prediction (Figure 10). Based on the difference between relevant stress value σ D and from Equation (9), sets of D − ε were calculated. Their respective correlation is presented in Figure 15b, showing the correlation between damage parameter and equivalent plastic strain used to model softening behavior of the material in this work under reference conditions. In Figure 14, a comparison between the experimental true stress-strain curves for different strain rates of Table 2 with the results obtained from a singular finite element are presented. Overall, good accuracy was achieved between experimental and numerical results with the exception of 1 s −1 strain rate results in which although the failure strain was similar, the ultimate stress value had a difference of approximately 8%. It is believed that this is due to the non-monotonic nature of failure strain as well as ultimate stress that was observed during the experimental procedure (Figure 6b). To accurately model the material softening behavior past the damage initiation point [27] represented in Figure 15a

Numerical Simulation of Full-Scale Tensile Specimens
To simulate the ductile failure of the dog-bone-shaped tensile specimens, a threedimensional model was constructed, replicating in detail the geometrical properties of the manufactured tensile specimens. After some initial mesh sensitivity analysis (Appendix C), the maximum size of C3D8R elements was chosen as equal to 3 mm × 2 mm × 2 mm (coarse mesh), and the minimum size at the working length was set as 1 mm × 1 mm × 1 mm (fine mesh), which resulted in a total of 25,344 elements ( Figure 16). For the material modeling, parameters stated in Section 3 were used, and boundary conditions were utilized in accordance with the experimental setup. Similarly to the experimental procedure, quasi-static 0.001 s −1 and 0.5 s −1 strain rate uniaxial tensile simulations were performed, and the results can be seen in Figure 17. Similarly to the experimental process, in the case of numerical simulations, load and displacement were monitored throughout the uniaxial tensile test. Afterwards, the obtained load-displacement data were transformed into engineering stress-strain curves and, subsequently, to true stress-strain curves. Overall, good accuracy is observed between experimental and numerical results for both cases. In

Numerical Simulation of Full-Scale Tensile Specimens
To simulate the ductile failure of the dog-bone-shaped tensile specimens, a threedimensional model was constructed, replicating in detail the geometrical properties of the manufactured tensile specimens. After some initial mesh sensitivity analysis (Appendix C), the maximum size of C3D8R elements was chosen as equal to 3 mm × 2 mm × 2 mm (coarse mesh), and the minimum size at the working length was set as 1 mm × 1 mm × 1 mm (fine mesh), which resulted in a total of 25,344 elements ( Figure 16). For the material modeling, parameters stated in Section 3 were used, and boundary conditions were utilized in accordance with the experimental setup. Similarly to the experimental procedure, quasistatic 0.001 s −1 and 0.5 s −1 strain rate uniaxial tensile simulations were performed, and the results can be seen in Figure 17. Similarly to the experimental process, in the case of numerical simulations, load and displacement were monitored throughout the uniaxial tensile test. Afterwards, the obtained load-displacement data were transformed into engineering stress-strain curves and, subsequently, to true stress-strain curves. Overall, good accuracy is observed between experimental and numerical results for both cases. In particular, for the quasi-static strain rate case, apart from some initial discrepancy in the postyielding stress capacity (numerical results overestimate the experimental case by 2.6%) the stress-strain curve follows closely the experimental results, and a 2.4% difference in ultimate stress is observed. Furthermore, although the fracture strain between the numerical analysis and experimental results is similar, the numerical simulation retains a higher stress state for larger strain values decreasing sharply, leading to element failure. For the case of the 0.5 s −1 strain rate, after yielding, the numerical model underestimates the experimental results by 2.7%, but after that, the numerical stress-strain curve follows closely the experimental one. Finally, a 3.6% difference is observed for the failure strain between numerical and experimental results. The above statements can be seen summarized in Table 5.   (Table 4) JC damage parameters.

Numerical Model Verifications for Smooth Specimens
To further validate the accuracy of the JC flow stress model as well as the JC damage parameters that were calculated in the previous sections, DIC analysis was attempted for specimens subjected to uniaxial tension under different strain rates. To perform the DIC analysis, commercially available software was used called "GOM Correlate". It is a DIC evaluation software program that is extensively used for material research and compound testing purposes. GOM Correlate follows a parametric approach that guarantees reliability for measuring required strains through a parametric approach. While using GOM Correlate, users have to define initial parameters regarding strain surface components for the area of interest, and the software creates facets which are square shaped across the whole range of supplied image data. Based on a high-contrast stochastic pattern that the user has to apply on the area of interest where strain is required to be measured, the software identifies these facets based on the quality of the aforementioned stochastic pattern. An additional parameter that requires adjustment from the user is the distance between adjacent facets on a center-to-center basis. This setting directly influences and it is correlated to the   (Table 4) JC damage parameters.

Numerical Model Verifications for Smooth Specimens
To further validate the accuracy of the JC flow stress model as well as the JC damage parameters that were calculated in the previous sections, DIC analysis was attempted for specimens subjected to uniaxial tension under different strain rates. To perform the DIC analysis, commercially available software was used called "GOM Correlate". It is a DIC evaluation software program that is extensively used for material research and compound testing purposes. GOM Correlate follows a parametric approach that guarantees reliability for measuring required strains through a parametric approach. While using GOM Correlate, users have to define initial parameters regarding strain surface components for the area of interest, and the software creates facets which are square shaped across the whole range of supplied image data. Based on a high-contrast stochastic pattern that the user has to apply on the area of interest where strain is required to be measured, the software identifies these facets based on the quality of the aforementioned stochastic pattern. An additional parameter that requires adjustment from the user is the distance between adjacent facets on a center-to-center basis. This setting directly influences and it is correlated to the  (Table 4) JC damage parameters.

Numerical Model Verifications for Smooth Specimens
To further validate the accuracy of the JC flow stress model as well as the JC damage parameters that were calculated in the previous sections, DIC analysis was attempted for specimens subjected to uniaxial tension under different strain rates. To perform the DIC analysis, commercially available software was used called "GOM Correlate". It is a DIC evaluation software program that is extensively used for material research and compound testing purposes. GOM Correlate follows a parametric approach that guarantees reliability for measuring required strains through a parametric approach. While using GOM Correlate, users have to define initial parameters regarding strain surface components for the area of interest, and the software creates facets which are square shaped across the whole range of supplied image data. Based on a high-contrast stochastic pattern that the user has to apply on the area of interest where strain is required to be measured, the software identifies these facets based on the quality of the aforementioned stochastic pattern. An additional parameter that requires adjustment from the user is the distance between adjacent facets on a center-to-center basis. This setting directly influences and it is correlated to the point density within the area of interest, and by decreased said distance, higher spatial resolution can be obtained by decreasing the distance of adjacent facets [28][29][30].
In this work, a full-frame CMOS camera was used to capture 1920 × 1080 pixelsized image series and videos to be used for the DIC analysis. For DIC, the captured stochastic pattern (Figure 18) was processed using a facet size of 14 and a distance of 9 to evaluate the corresponding strain fields. To compare the experimental results with the numerical modeling, strain was recorded along the axis of the cylindrical specimens and was compared with the resulting strain in central nodes of the FE model. Due to camera limitations regarding video framerate capturing capabilities, as well as applied pattern-related shortcomings, although DIC analysis was attempted for all experimental cases, it was only successful for quasi-static and 0.5 s −1 experimental cases. Namely, for higher strain-rates, to capture a significant amount of image series, 120 fps video recording was attempted, but only a resolution up to 1280 × 720 pixels was available. The lowered resolution coupled with brittleness of the coloring used to create the stochastic speckle pattern severely limited the amount of usable data for the DIC analysis. The resulting strain profile can be seen in Figure 19 for (a) the quasi-static tensile case and for (b) the dynamic case with a strain rate of 0.5 s −1 . In both Figure 19a,b, the horizontal axis data for finite element modeling (FEM) case have been shifted by an amount suitable to each case in order to align the portion of the gauge length in which necking and, correspondingly, fracture occurred. To aid with the visualization of the data, true strain and normalized gauge length data points are being used. For the quasi-static case, due to stochastic pattern degradation, strain data were captured up to approximately 0.2 strain. Appl modeling, strain was recorded along the axis of the cylindrical specimens and was compared with the resulting strain in central nodes of the FE model. Due to camera limitations regarding video framerate capturing capabilities, as well as applied pattern-related shortcomings, although DIC analysis was attempted for all experimental cases, it was only successful for quasi-static and 0.5 s −1 experimental cases. Namely, for higher strain-rates, to capture a significant amount of image series, 120 fps video recording was attempted, but only a resolution up to 1280 × 720 pixels was available. The lowered resolution coupled with brittleness of the coloring used to create the stochastic speckle pattern severely limited the amount of usable data for the DIC analysis. The resulting strain profile can be seen in Figure 19 for (a) the quasi-static tensile case and for (b) the dynamic case with a strain rate of 0.5 s −1 . In both Figure 19a,b, the horizontal axis data for finite element modeling (FEM) case have been shifted by an amount suitable to each case in order to align the portion of the gauge length in which necking and, correspondingly, fracture occurred.
To aid with the visualization of the data, true strain and normalized gauge length data points are being used. For the quasi-static case, due to stochastic pattern degradation, strain data were captured up to approximately 0.2 strain.
(a) (b) Figure 18. Stochastic speckle pattern used for DIC analysis during (a) beginning of the tensile test and (b) one frame before rupture.
When comparing the FEM analysis results with DIC for 0.001 s −1 strain rate, around the area where necking occurs, high strain region is concentrated for 22.6% of the normalized length as opposed to the DIC, which is 26.5%, resulting in the tensile specimen forming a longer necking region by 3.9%. For the case of 0.5 s −1 strain rate, even better accuracy  In Figure 20, the strain map results from DIC analysis and FE simulations are overlapped visualizing the results presented in Figure 19 for smooth specimens at the last captured frame before rupture. In each corresponding figure, FE simulation strain contours have been scaled to the value of DIC analysis, and the overlapping FEM image has been repositioned to align the necking region with the DIC image. As a result, in the case of 0.001 s −1 strain rate, near the center of the necking region, the resulting strain values exceed the visualization boundary limits of 0.0-0.2 strain and, as a result, the region with a strain rate higher than these values are presented in grey color. When comparing the FEM analysis results with DIC for 0.001 s −1 strain rate, around the area where necking occurs, high strain region is concentrated for 22.6% of the normalized length as opposed to the DIC, which is 26.5%, resulting in the tensile specimen forming a longer necking region by 3.9%. For the case of 0.5 s −1 strain rate, even better accuracy is observed with the FEM results underestimating the length where necking occurs by 1.6%. Overall, in both cases, good accuracy is observed between the DIC and FEM results, further reinforcing the suitability of the proposed JC model and damage parameters for SBPR 930/1080 Type B No. 1 tendon high-strength material.
In Figure 20, the strain map results from DIC analysis and FE simulations are overlapped visualizing the results presented in Figure 19 for smooth specimens at the last captured frame before rupture. In each corresponding figure, FE simulation strain contours have been scaled to the value of DIC analysis, and the overlapping FEM image has been repositioned to align the necking region with the DIC image. As a result, in the case of 0.001 s −1 strain rate, near the center of the necking region, the resulting strain values exceed the visualization boundary limits of 0.0-0.2 strain and, as a result, the region with a strain rate higher than these values are presented in grey color. Appl

Conclusions and Recommendations
In this work, numerous tensile tests at room temperature and strain rates 0.001 s −1 -0.5 s −1 were performed in order to calculate the Johnson-Cook model and damage parameters for SBPR 930/1080 Type B No. 1 tendon material aimed at fracture analysis of posttensioned concrete members. Overall, the results obtained in this work, after the calculation of JC parameters, showed good agreement with the experimental data. To verify the agreement between numerical data and experimental tensile data, commercial FEM software was used. The experimental tensile tests were replicated in detail in order to verify the performance of the damage model using both experimental observations, experimental stress-strain data, as well as DIC analysis. The JC model is found to be able to closely predict experimental data with less effort in comparison to other analytical models, but on the other hand, to properly calibrate the related parameters, extensive numerical data are required from several experimental cases. It was found that besides small prediction differences between the FE simulation and experimental results, good accuracy was achieved in predicting the effects of strain concentration and geometrical deformation (necking). Based on these research outcomes, the followed procedure can be applied to closely predict the performance of tested material for fracture analysis of post-tensioned concrete members.
As for recommendations, in this work, all testing was performed at room temperature, and thus, we were not able to identify damage-related parameters with regard to material softening under elevated temperatures; therefore, further testing could be performed to identify these parameters, and dynamic testing could be performed under even higher strain rates, in order to gain a broader picture of strain-related hardening effects. Furthermore, by utilizing the JC model and damage parameters presented in this work, researchers can simulate and estimate the remaining strength of a damaged or corroded PC tendon. By measuring or estimating the remaining cross-section of a tendon, it is possible to accurately estimate stress-concentration states surrounding the damaged or corroded area and construct relevant fragility curves based on the anticipated loading conditions.
Author Contributions: I.G.: conceptualization; methodology; investigation; data analysis and curation; validation; writing-original manuscript. N.C.: supervision, reviewing and editing. All authors have read and agreed to the published version of the manuscript.

Conclusions and Recommendations
In this work, numerous tensile tests at room temperature and strain rates 0.001 s −1 -0.5 s −1 were performed in order to calculate the Johnson-Cook model and damage parameters for SBPR 930/1080 Type B No. 1 tendon material aimed at fracture analysis of post-tensioned concrete members. Overall, the results obtained in this work, after the calculation of JC parameters, showed good agreement with the experimental data. To verify the agreement between numerical data and experimental tensile data, commercial FEM software was used. The experimental tensile tests were replicated in detail in order to verify the performance of the damage model using both experimental observations, experimental stress-strain data, as well as DIC analysis. The JC model is found to be able to closely predict experimental data with less effort in comparison to other analytical models, but on the other hand, to properly calibrate the related parameters, extensive numerical data are required from several experimental cases. It was found that besides small prediction differences between the FE simulation and experimental results, good accuracy was achieved in predicting the effects of strain concentration and geometrical deformation (necking). Based on these research outcomes, the followed procedure can be applied to closely predict the performance of tested material for fracture analysis of post-tensioned concrete members.
As for recommendations, in this work, all testing was performed at room temperature, and thus, we were not able to identify damage-related parameters with regard to material softening under elevated temperatures; therefore, further testing could be performed to identify these parameters, and dynamic testing could be performed under even higher strain rates, in order to gain a broader picture of strain-related hardening effects. Furthermore, by utilizing the JC model and damage parameters presented in this work, researchers can simulate and estimate the remaining strength of a damaged or corroded PC tendon. By measuring or estimating the remaining cross-section of a tendon, it is possible to accurately estimate stress-concentration states surrounding the damaged or corroded area and construct relevant fragility curves based on the anticipated loading conditions. Calculations and parameter estimations in Sections 3 and 4 were performed considering the damage initiation point of Figure 10 and subsequent material softening described in Section 4.1 and Figure 15b. Following that procedure, FEM simulations up to the point of fracture are able to closely match the experimental data both on a singular element basis ( Figure 14) as well as full model simulations ( Figure 17). When using fracture point as the basis for calculating JC damage parameters, the correlation between fracture strain and stress triaxiality from uniaxial tensile test data can be seen in Figure A1. The curve, although exponential shaped, has a distinctively different shape as a result of the D1-D3 parameters. Similarly, when inputting newly calculated parameters into Equation (11), the angle of linear fit equation as well as corresponding data points result in different D4 parameter, as can be seen in Figure A2.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
Calculations and parameter estimations in Sections 3 and 4 were performed considering the damage initiation point of Figure 10 and subsequent material softening described in Section 4.1 and Figure 15b. Following that procedure, FEM simulations up to the point of fracture are able to closely match the experimental data both on a singular element basis ( Figure 14) as well as full model simulations ( Figure 17). When using fracture point as the basis for calculating JC damage parameters, the correlation between fracture strain and stress triaxiality from uniaxial tensile test data can be seen in Figure A1. The curve, although exponential shaped, has a distinctively different shape as a result of the D1-D3 parameters. Similarly, when inputting newly calculated parameters into Equation (11), the angle of linear fit equation as well as corresponding data points result in different D4 parameter, as can be seen in Figure A2. The calculated JC damage model parameters for fracture point-based calculation are summarized in Table A1, and similarly to Section 3.3, they can be used in FE software to simulate yield and fracture of high-strength tendons but with an evident overestimation of their corresponding softening behavior. With the parameters of Table A1, a comparison similar to Figure 14 can be seen between experimental true stress-strain curves for different strain rates of Table 2, with the results obtained from a singular finite element. It is evident that due to the usage of fracture strain instead of corresponding damage initiation value, larger discrepancies are observed between the experimental FE simulation results, especially in the case of the 0.001 s −1 strain rate. The calculated JC damage model parameters for fracture point-based calculation are summarized in Table A1, and similarly to Section 3.3, they can be used in FE software to simulate yield and fracture of high-strength tendons but with an evident overestimation of their corresponding softening behavior.  Table A1, a comparison similar to Figure 14 can be seen between experimental true stress-strain curves for different strain rates of Table 2, with the results obtained from a singular finite element. It is evident that due to the usage of fracture strain instead of corresponding damage initiation value, larger discrepancies are observed between the experimental FE simulation results, especially in the case of the 0.001 s −1 strain rate. Figure A3. Comparison between experimental (Exp, continuous lines) and FE simulation (Sim, dotted lines) results for 0.001 s− 1 -2s− 1 strain rates using JC damage parameters of Table A1.
When performing numerical simulations of full-scale tensile specimens similar to Section 4.2, it can be seen that using damage parameters of Table A1, both ultimate stress state and rupture strain are being overestimated in comparison to the experimental results, and thus, the JC damage parameters of Section 3.3 are recommended for usage in FE software applications ( Figure A4). Exp 0.001s -1 Sim 0.5s -1 Exp 0.5s -1 Sim 1s -1 Exp 1s -1 Sim 2s -1 Exp 2s -1 Figure A3. Comparison between experimental (Exp, continuous lines) and FE simulation (Sim, dotted lines) results for 0.001 s −1 -2 s −1 strain rates using JC damage parameters of Table A1.
When performing numerical simulations of full-scale tensile specimens similar to Section 4.2, it can be seen that using damage parameters of Table A1, both ultimate stress state and rupture strain are being overestimated in comparison to the experimental results, and thus, the JC damage parameters of Section 3.3 are recommended for usage in FE software applications ( Figure A4). Figure A4. True stress-strain plot comparison for 0.001 s− 1 and 0.5 s− 1 strain rates for average experimental and numerical simulation using JC damage parameters of Table A1.

Appendix B
Regarding the calculation procedure for the C parameter followed in Figure 9 and Equation (5), in order to reduce the influence of quasi-static data (ln * = −6.90), 1/5th of the data was removed, and the remaining data were plotted again in Figure A5. As it can be seen from the comparison of Figures 9 and A5, the influence of the amount of data used is minimal for the calculation of the C parameter since the data accumulation trend remains similar.  Table A1.

Appendix B
Regarding the calculation procedure for the C parameter followed in Figure 9 and Equation (5), in order to reduce the influence of quasi-static data (ln ε * = −6.90), 1/5th of the data was removed, and the remaining data were plotted again in Figure A5. As it can be seen from the comparison of Figures 9 and A5, the influence of the amount of data used is minimal for the calculation of the C parameter since the data accumulation trend remains similar.  Table A1.

Appendix B
Regarding the calculation procedure for the C parameter followed in Figure 9 and Equation (5), in order to reduce the influence of quasi-static data (ln * = −6.90), 1/5th of the data was removed, and the remaining data were plotted again in Figure A5. As it can be seen from the comparison of Figures 9 and A5, the influence of the amount of data used is minimal for the calculation of the C parameter since the data accumulation trend remains similar.

Appendix C
To ensure that the results and findings of this work are applicable to other FE model simulations, mesh sensitivity analysis was performed for the case of a full-scale 3D model. According to the Abaqus Analysis User's Manual [27], the stress-strain relationship that is used to define material behavior, can no longer represent material behavior after the onset of material damage (stated as damage initiation point in Figure 10). If the finite element model were to continue following the behavior defined in the stress-strain relationship, a strong mesh dependency would occur based on strain localization. In order to overcome this issue, Abaqus uses a different approach (Damage Evolution Law) to model material softening behavior past the damage initiation point. Specifically, Hillerborg et al. fracture energy proposal [31] is adopted, which decouples mesh dependency from material behavior once damage is initiated. In their proposal, fracture energy is defined according to Equation (A1).
where is the equivalent plastic displacement as work per unit area of crack that has formed. Prior to damage initiation point, is considered as zero, and after that, it is calculated based on Equation (A2).

= ̅ (A2)
where L is defined as the characteristic element length, and for 3D elements used in this work, it is calculated as the ratio of element volume to area of the largest face of the element ( = . . ⁄ ). To illustrate the mesh independency from the JC model and damage parameters (Tables 3 and 4), in Figure A7, a comparison can be seen between average experimental results for the 0.001 s −1 strain rate from Table 2 and FE simulation for average mesh sizes of 1 × 1 × 1 ( = 1) and 2 × 2 × 2 ( = 2) mm that were used in the gauge length region of the full-scale FE model.

Appendix C
To ensure that the results and findings of this work are applicable to other FE model simulations, mesh sensitivity analysis was performed for the case of a full-scale 3D model. According to the Abaqus Analysis User's Manual [27], the stress-strain relationship that is used to define material behavior, can no longer represent material behavior after the onset of material damage (stated as damage initiation point in Figure 10). If the finite element model were to continue following the behavior defined in the stress-strain relationship, a strong mesh dependency would occur based on strain localization. In order to overcome this issue, Abaqus uses a different approach (Damage Evolution Law) to model material softening behavior past the damage initiation point. Specifically, Hillerborg et al. fracture energy proposal [31] is adopted, which decouples mesh dependency from material behavior once damage is initiated. In their proposal, fracture energy is defined according to Equation (A1). where L is defined as the characteristic element length, and for 3D elements used in this work, it is calculated as the ratio of element volume to area of the largest face of the element (L = Vol./L. Area). To illustrate the mesh independency from the JC model and damage parameters (Tables 3 and 4), in Figure A7, a comparison can be seen between average experimental results for the 0.001 s −1 strain rate from Table 2 and FE simulation for average mesh sizes of 1 × 1 × 1 (L = 1) and 2 × 2 × 2 (L = 2) mm that were used in the gauge length region of the full-scale FE model. It can be seen from Figure A7 that the FE analysis results closely match regardless of the mesh size that was utilized. The results of this analysis were expected and can be used to further validate the element damage evolution law followed by Abaqus. In Figure A8, strain map results similar to Figure 20, can be seen for DIC, and average mesh sizes of 1 × 1 × 1 and 2 × 2 × 2 mm are aligned around the necking region for smooth specimens at the last captured frame before rupture. It can be seen that regardless of the mesh size, similar results are obtained, but due to the decreased number of elements along the gauge length, small strain discrepancies are observed.  Figure A8. Comparison between DIC and FE simulation strain map results for smooth specimens at the last captured frame before rupture at 0.001 s −1 strain rate for 1 × 1 × 1 and 2 × 2 × 2 mm average mesh sizes. Figure A7. True stress-strain plot comparison for 0.001 s −1 strain rate for average experimental and numerical simulation using average mesh sizes of 1 × 1 × 1 and 2 × 2 × 2 mm.
It can be seen from Figure A7 that the FE analysis results closely match regardless of the mesh size that was utilized. The results of this analysis were expected and can be used to further validate the element damage evolution law followed by Abaqus. In Figure A8, strain map results similar to Figure 20, can be seen for DIC, and average mesh sizes of 1 × 1 × 1 and 2 × 2 × 2 mm are aligned around the necking region for smooth specimens at the last captured frame before rupture. It can be seen that regardless of the mesh size, similar results are obtained, but due to the decreased number of elements along the gauge length, small strain discrepancies are observed. It can be seen from Figure A7 that the FE analysis results closely match regardless of the mesh size that was utilized. The results of this analysis were expected and can be used to further validate the element damage evolution law followed by Abaqus. In Figure A8, strain map results similar to Figure 20, can be seen for DIC, and average mesh sizes of 1 × 1 × 1 and 2 × 2 × 2 mm are aligned around the necking region for smooth specimens at the last captured frame before rupture. It can be seen that regardless of the mesh size, similar results are obtained, but due to the decreased number of elements along the gauge length, small strain discrepancies are observed.  Figure A8. Comparison between DIC and FE simulation strain map results for smooth specimens at the last captured frame before rupture at 0.001 s −1 strain rate for 1 × 1 × 1 and 2 × 2 × 2 mm average mesh sizes.